Patent ID: 7505853

Claim:
A method for analysis or inspection of a multi-ply composite material with anisotropic layers comprising: (a) determining properties of the multi-ply composite material, wherein the properties comprise the ply lay up, the number of plies in the material, the fiber orientation in the plies and/or the thicknesses of the plies; (b) computing a stiffness matrix in global coordinates for each ply in the multi-ply composite material using a Christoffel equation, and changing from one coordinate system to another using calculations based on a transformed version of Christoffel matrix elements and a Christoffel tensor, and using the stiffness matrix to determine how a wave should travel and how the wave travels in the material; (c) computing a common reference using an algorithm comprising a Christoffel equation; (d) making an initial guess of an ultrasonic wavenormal in the first ply of the composite material and calculating wave phase velocity by assuming that the wavenormal is a straightline path even though the wavenormal actually bends at each interface; (e) determining the angular deviation of the first guess wavenormal from an actual ultrasonic wavenormal by determining how far away the initial guess of step (d) is from the actual wavenormal by (1) determining the actual ultrasonic wavenormal by determining where the initial guess of the wavenormal in the first ply is being sent from, applying the boundary conditions and the reflection and refraction at each of the ply interfaces, and determining where the ultrasonic wavenormal actually ends up; (2) determining how far away the initial guess of step (d) is from where the actual ultrasonic wavenormal as determined in (e)(1) by determining the angular deviation; (f) adjusting the initial guess wavenormal of step (d) using the angular deviation determination of step( e); (g) computing a new wavenormal and its phase velocity by applying Snell's law for reflection at each interface in the composite material to compute the new wavenormal and its phase velocity in each ply layer of the multi-ply composite material until the focus depth is reached; (h) calculating the group velocity and direction in each layer and summing all the vectors; (i) determining if the ray path of the ultrasonic wave goes to the focus depth coordinates; and (j)(1) if the ray path of the ultrasonic wave does go to the focus depth, the transit time to the focus coordinates and back is computed and the method comprises a further step comprising moving the to the next probe or focus coordinates and repeating (iterating) the method at step (g) where Snell's law is utilized to calculate the wavenormal and its phase velocity in each layer below, until the focus depth is reached; or (2) if the ray path of the ultrasonic wave does not go to the focus depth, the wavenormal is adjusted from the initial guess and its phase velocity is calculated and stored and/or presented to a user.