Stresses are solved for two parallel cracks in an infinite orthotropic plate during passage of incoming shock stress waves normal to their surfaces. Fourier transformations were used to reduce the boundary conditions with respect to the cracks to two pairs of dual integral equations in the Laplace domain. To solve these equations, the differences in the crack surface displacements were expanded to a series of functions that are zero outside the cracks. The unknown coefficients in the series were solved using the Schmidt method so as to satisfy the conditions inside the cracks. The stress intensity factors were defined in the Laplace domain and were inverted numerically to physical space. Dynamic stress intensity factors were calculated numerically for selected crack configurations.