QiU Jinhao，TAO Chongcong，Ji Hongli，ZHANG Chao，ZHAO Jinling,3

1.State Key Laboratory of Mechanics and Control of Mechanical Structures，Nanjing University of Aeronautics and Astronautics，Nanjing 210016，P.R.China；

2.College of Aerospace Engineering，Nanjing University of Aeronautics and Astronautics，Nanjing 210016，P.R.China；

3.School of Mechanical and Power Engineering，Nanjing Tech University，Nanjing 211800，P.R.China

Abstract: Laser ultrasonic technique has received increasing attentions in the past decade due to its contactless nature and a wide range of applications have been reported. in this review，applications of laser ultrasonic technique developed at Nanjing University of Aeronautics and Astronautics（NUAA）as well as elsewhere for non-destructive testing in composite laminates are presented. The principles of generating and detecting in a laser ultrasonic system are introduced，three different system configurations are also introduced with each configuration's advantages and disadvantages explained. More importantly，two major applications developed at NUAA for composite laminates are presented including damage detection，stiffness reconstruction and fatigue life prediction. Both applications are realized by a fixed-point PZT sensor and scanning pulse laser based on the linear reciprocal theorem. Analytical method and numerical models are employed and developed to realize the functionalities.

Key words: laser ultrasonic;composite laminates;damage detection;stiffness reconstruction

0 introduction

Fiber reinforced polymers（FRP）are widely used in engineering applications recently due to their superior mechanical properties over traditional metallic materials such as damage tolerance，fatigue and corrosion resistance and especially their high specific stiffness and strength. However，traditional composite structures e.g. composite laminates also come with their weaknesses such as sensitivity to impact damages，non-critical damage developments under relative low loads and composite laminates also show decrease in load-bearing capacity when subjected to cyclic loads. Therefore，various non-destructive testing（NDT）methods have been proposed and tested for damage detections in composite laminates，among which，laser ultrasonic technique shows great potentials due to its contactless nature.

Ultrasonic is a widely used approach for damage detection，the most traditional one being ultrasonic C-scan［1］. One or more probes are usually used to generate and sense bulk ultrasonic waves，which propagate in the structure in the thickness direction as shown in Fig.1（a）. The collected data of transmitted and reflected waves are processed in the time domain or frequency domain to detect the internal damages in structures. However，the detection result，which depends on the experience and skill of inspectors，suffers from the complex interference of multiple reflected and diffracted waves in thin platelike structure. in addition，direct-contact scan and indirect-contact scan through the media（such as water and oil）are costly and labor intensive，especially in large-scale inspection applications.

Apart from ultrasonic C - scan approach，the Lamb wave-based NDE techniques that have been widely used in plate-like structures also show appealing features in detecting damages because Lamb wave has the advantages of high sensitivity to various damage types and the capability of long transmission distance［2-3］. and in general，Lamb wave can be generated and received by piezoelectric sensor pairs［4］that are distributed over the inspection region，as shown in Fig.1（b）. By analyzing the changes of Lamb wave propagation characteristics or using the time-of-flight information of the scattered waves，the damage position can be detected. However，contact-type transducers that are usually surface mounted or embedded on the test structures lead to some limitations of this method in wide applications that：（1）the excitation and sensing positions are fixed at several discrete points. To improve the inspection quality and extend the inspection region，a large amount of piezoelectric sensors are required；（2）transducer installation and cabling can be costly and labor intensive，especially in the large-scale inspection applications；（3）many contact transducers are not applicable under harsh environments such as high temperature and radioactive conditions；and（4）the added weight of the transducers and cables is not cost-effective in practical applications.

Fig.1 Different inspection principles using bulk and lamb waves

Laser ultrasonic technique uses pulse laser to generate ultrasonic waves and overcomes most of the aforementioned disadvantages. The generation of ultrasound with lasers can actually go back to nearly as far as the invention of the first ruby laser，which was used to generate ultrasound or shock waves in materials［5-6］. Ultrasonic can be excited by laser through two different approaches：a thermoelastic way or a vaporized way，the former is truely non-destructive while the latter can cause surface damages. Since the scope of this review is focused on applications in NDT fields，only the thermoelastic approach will be included. The principle of ultrasonic excitation based on thermoelastic mechanism is to heat the material locally to expand，thereby producing a strain and a corresponding stress that becomes the source of wave propagating［7］. The heated point inside the material is a center of expansion to which 3 orthogonal force dipoles can be associated［8-9］，thus the material effectively becomes the emitting transducer of ultrasound. Besides emitting transducer［10-13］，laser can also be used as receiving transducer with an interferometer，which results in three different configurations for a laser ultrasonic system：fixed-point PZT sensor with scanning pulse laser，fixed-point PZT excitation with scanning laser Doppler vibrometer（LDV）sensing and pulse laser excitation with LDV sensing，the first two being a semi-laser system and the third a true laser system. Brief descriptions of the three systems are as follows.

Strategy 1Fixed-point PZT sensor and scanning pulse laser—The system of the first scheme in Fig.2 consists of a Q-switched diode-pumped solidstate laser，a galvanometric laser mirror scanner，an acoustic emission（AE）sensor，a programmable filter，a digitizer，and a computer for hardware control and image processing.A schematic diagram of the entire experimental system configuration is shown in Fig.3. The AE sensor providing ultrasonic reception is fixed and focuses at a particular sensing region on the specimen，whereas the laser beam acting as ultrasonic generation is directed to the targeted specimen with the help of a laser mirror scanner for remote and controllable scanning of the targeted area.The filters，digitizer，and computer hardware are used for data acquisition and image processing for damage evaluation and visualization.All devices are packed in a controller.

Fig.2 Four different laser scanning schemes for wavefield visualization[14]

Fig.3 Experimental setup

During the data acquisition process，the different ultrasonic responses between movable ultrasonic source point A and fixed sensor point B can be obtained experimentally as Case 1 in Fig.4. Based on the linear reciprocal theorem［15-16］the wave propagation from point A to B can be directly converted into the propagation from point B to A. Thus，the obtained data in Case 1 can represent the wave propagation from the fixed laser point B to an inspection region，which is distributed with a series of AE sensors as shown in Case 2 from Fig.4. Then，the responses of each scanning point at the same time t are plotted on an intensity snapshot to represent the wavefield at time point t. Finally，the snapshots are continually displayed in a time series and show the wave propagation in the inspection region as shown in Fig.5［17］.

Fig.4 Schematic diagram of laser ultrasonic scanning technique

Fig.5 Data structure in wave propagation visualization

Strategy 2Fixed-point PZT excitation and scanning LDV sensing—Different from Strategy 1，the data acquisition from this inspection strategy can be directly regarded as the Case 2 in Fig.4. The waveform generated by the piezoelectric transducer can be controlled by the input signal，which is more flexible in inspection application. When the sensing laser beam is perpendicular to the polished target surface，the signal to noise ratio（SNR）of the ultrasonic signal is high due to the majority of the incident laser beam being reflected straight back to the receiver. However，the SNR of the ultrasonic decreases as the incident angle of the laser beam increases. Thus，multiple ultrasonic signals often need to be measured at a single sensing point and averaged to improve the SNR［18］.

Strategy 3LDV sensing and pulse laser excitation—When fixing the LDV sensing point and scanning the pulse laser in the inspection region，the ultrasonic propagation can be constructed sharing the same principle with Strategy 1. On the contrary，the situation that is to scan the LDV sensing point and fix the pulse laser is the same with Strategy 2.Both of them can achieve a complete noncontact inspection.

1 Damage Visualization Based on Laser Ultrasonic imaging

A major application and advantage for laser ultrasonic in the field of NDT is the capability to achieve damage detection and visualization for both metallic and composite materials［19-31］. Lee et al.［20］proposed an ultrasonic propagation imaging for automatic damage visualization for steel pipeline with multiple kinds of damages in nuclear power plant. in their later work，Lee et al.［22］designed a laser ultrasonic system with functions of excitation，sensing，and diagnosis for damage visualization of target located at a long distance. A fixed-point PZT sensor and scanning pulse laser configuration was used.Sohn et al.［23］studied the applicability of non-contact wavefield imaging techniques to detect delamination and debond in composite structures，where a fixedpoint PZT excitation and scanning LDV sensing configuration was used. Martarelli et al.［26］explored the applicability of delamination detection with laser ultrasonic. Chia et al.［24］emended PZT sensors into aircraft's wing structure made of CFRP material and utilized scanning laser to achieve visualization of damages and manufacture defects. An et al.［14］developed a complete noncontact laser ultrasonic scanning system with an LDV sensing and pulse laser excitation configuration，the system was applied for crack visualization and diagnosis in an aluminum plate. Liu et al.［30］also successfully located the fatigue crack in an aluminum plate using laser ultrasonic. in a more recent work，Park et al.［27］applied the noncontact laser ultrasonic system to detect delamination and debond in both a CFRP composite aircraft wing structure and a GFRP composite wind blade structure. Laser ultrasonic can also be used with other non-contact techniques to form a hybrid non-contact inspection system. Cavuto et al.［29］combined laser ultrasonic and air-coupled probe to improve the performance of train axle inspection. Dixon et al.［21］used laser ultrasonic in conjunction with electromagnetic acoustic transducer（EMAT）to detect cracks in metal sheets. in terms of real-world applications，Lee et al.［25］used optical fiber to direct the laser in a pipeline system for nuclear power plants to detect typical pipe defects without baseline data. Bate et al.［31］measured the texture in metal at elevated temperature using laser ultrasonic. Sun and Zhou［28］used laser ultrasonic to detect drill-induced delamination.

in the following，a damage visualization algorithm based on anomalous incident wave（AiW）developed at NUAA will be introduced to help readers better understand the principles of laser ultrasonicbased damage visualization.

1.1 Ultrasonic wave in damaged region

Assuming that the ultrasonic wave passes through a damage as shown in Fig.6. The incident wave wi（x，t）propagates along the positive x-direction and the reflected wave wr（x，t） propagates along the negative x-direction. Generally，the harmonic wave can be given by the following equations［14］

Fig.6 Schematic diagram of the ultrasonic wave interacting with a damage

where ω is the angular frequency，k is the wavenum-ber，x0and l0represent the position and width of the slit，α and β are the phase shifts caused by the damage，A represents the incident wave amplitude while B and C represent the amplitudes of residual incident wave by transmission and reflected wave.

1.2 Damage detection algorithm using AiW energy map

To represent the discontinuous characteristic of the incident wavefield in the damage area，the AiW is defined as the differences between two incident wave signals obtained at adjacent points. Taking two incident wave signals wi（x，i）and wi（x+Δx，i）as an example，the AiW at the position x can be calculated as

where Δx is the minimum spatial resolution of the laser scanning process and the time-lag drmaxcan be calculated through a cross-correlation r（d）.

where drmaxrepresents the time that the incident wave travels from x to x+ Δx to make r（drmax）reach the maximum value. Considering that the incident wave propagates in the continuous medium，two incident waves at the adjacent points are similar with each other. By using Eq.（3），the traveling incident wave in regions 1 and 3 as shown in Fig.6 can be eliminated. As a result，the energy changes in region 2 caused by the damage will be highlighted and the discontinuous characteristic of Eq.（1）is extracted. To evaluate the size and shape of the damage by a digital image，the AiW energy EAiWcan be calculated by the sum of squared AiW［32］

Compared with calculating anomalous wave directly from the original wavefield［17，33］，the AiW extracts the wave propagation changes caused by the damage as well. However，due to the fact that the incident wave is much stronger than reflected wave in terms of signal strength，drmaxonly denotes the time delay of the incident waves. Because the reflected wave and incident wave have the opposite time delays between the adjacent points，the reflected wave is also enlarged by this process. As the reflected wave propagating towards the health area in region 1，the energy of the anomalous wave is also distributed near the damaged area，which reduces the resolution of the damage image. After filtering the reflected wave，the AiW energy eliminates the reflected wave energy in health area that makes the shape of the damage much clearer.

As shown in Fig.6，a part of the incident wave is reflected with the propagation direction reversed.Because the damage is nonpenetrated，the other part of the incident wave passes through the damage and propagates in the original direction. Thus，the one-dimensional wave field w（x，t）near the damage area can be given by combining the incident wave field wi（x，t）and reflected wave field wr（x，t）.

To construct the AiW map，the incident wave field wishould be extracted from the original wave field w. From Eqs.（1，2），the major difference between the incident and reflected waves is the sign in front of the wavenumber k that indicates the wave propagation direction. To obtain the incident wave，twodimensional Fourier transform （2D-FT） converts the wavefield from space-time（x-t） domain into wavenumber-frequency （k-ω） domain to separate the waves with different propagation directions as［34-35］

where W（k，ω）is the wavefield in wavenumberfrequency domain and“j”represents the imaginary unit. The wavefield is divided into two parts：（1）the incident wave part that is in the quadrant with kω＜0 and（2）the reflected wave part that is in the quadrant with kω ＞0. Two - dimensional inverse Fourier transform can then be applied to recover the time domain signals of these two components as follows

where Φi（R）is a window function defined as follows，used to separate the incident and reflected waves.

The principle is explained using 1D wave-field based 2D-FT as above. Similarly，the above 1D wave decomposition method can be extended to 2D case by means of 3D Fourier transform（3D-FT）and three - dimensional inverse Fourier transform（3D-iFT），for obtaining the wavefield by a laser ultrasonic scanning technique. Because the sensor response is collected by digital acquisition，the wavefield takes the form of discrete data. Thus，the Fourier transform can be realized by fast Fourier transform（FFT）.

1.3 Experimental validation

To examine the feasibility of the proposed laser ultrasonic scanning system and image processing technique，the scheme using fixed-point AE sensor and scanning pulse laser as shown in Fig.2 is employed for ultrasonic imaging and damage evaluation.Experiments of a composite structure were conducted.

A symmetric CFRP laminated plate with a stacking configuration of［45°/-45°/0°/90°］sis used to validate the damage evaluation method［32］.The inplane dimensions of the CFRP plate are 480 mm×480 mm and the thickness is 1 mm. To make an artificial damage，a drop-weight impact test was performed according to Test Method D7136［36］. Damage is imparted through out - of -plane，concentrated impact（perpendicular to the plane of the CFRP plate）using a drop weight with a hemispherical striker tip as shown in Fig.7. The impact energy Eimpactcan be calculated by

where m is the mass of the drop weight，g is the acceleration due to gravity and h is the drop-height. in this validation，the impact energy wasset to be 10 J and the impact damage cannot be detected visually.However，delamination，split and crack should exist inside the structure.

Fig.7 Schematic diagram of the drop-weight impact test for a CFRP plate

Fig.8 Picture of the CFRP plate with a damage caused by a 10 J impact

As shown in Fig.8，the impacted side of the CFRP laminated plate is scanned by the movable laser point and the inspection region is 100 mm ×100 mm. The AE sensor is placed on the opposite side of the CFRP plate with a distance of 12 cm away from the centre of the inspection region. The spatial interval of the laser scan is 1 mm and the signal is filtered with a band-pass frequency from 100 kHz to 200 kHz.

Two wavefields at the different time points are shown in Fig.9. Compared with the ultrasonic wave in a metal structure，the ultrasonic wave in the CFRP plate has low SNR due to the large material damping and the wave velocity is also lower. in order to ensure both of the S0and A0waves can pass through the inspection region，the sampling period is set to 200 μs. Same with the phenomenon in a metal plate with a slit，the amplitude changes around the damage region which can be used to detect the position of the damage.

Fig.9 Experimental wavefields at different times

The original wave energy map is constructed as shown in Fig.10. The wave energy decreases as the wave propagates. The anomalous wave energy region hardly shows the impact damage. To evaluate the inspection result，F(xiàn)ig.11 shows the damage image obtained by a conventional immersion ultrasonic C-scan method. Compared with Fig.11，the large amount of the wave energy in the“health”area reduces the resolution of damage image in Fig.10.

After filtering the reflected wave，the AiW energy map can be obtained according to Eq.（5）.Compared with calculating anomalous wave directly from the original wavefield，AiW extracts the wave propagation changes caused by the impact damage，and removes the negative influence from the reflected waves. Fig.12 shows the damage evaluation result using AiW energy map.

Fig.10 Damage evaluation using the original wave energy map

Fig.11 Damage evaluation using immersion ultrasonic C-scan method

Fig.12 Damage evaluation using the AiW energy map

2 Stiffness Reconstruction

identification of material's mechanical properties for composites is fundamental in laminate structures design，and it also provides a nondestructive approach for evaluating composite fatigue accumulated damages［37-38］. The stiffness of laminates may change due to stress or heat induced fatigue damages. Evaluating material properties of the composites periodically thus enables monitoring the health of the structures. Conventional techniques like tensile，compression and shear tests can be destructive in nature and provide only a few of the elastic moduli. Ultrasonic techniques are more advantageous over the conventional methods due to their qualification for nondestructive measurement of the elastic moduli of such materials［39-45］.

Ultrasonic-based stiffness measurement methods come with two categories：the bulk wave and the Lamb wave. Using the bulk wave to measure certain stiffness parameter is well known［46-59］. however，this method is insensitive to the in-plane C12and out-of-plane C66. Hence only seven stiffness parameters out of nine for composite materials can be effectively reconstructed. Another limitation of the bulk wave method is that separate samples are need-ed，whereas Lamb wave-based techniques are in situ.

Lamb wave method has also been investigated to a certain degree［51，60-62］. Jen et al.［61］investigated the inverse problem for isotropic plates but did not conduct any experimental study. Rao［62］reconstructed the elastic moduli of a unidirectional composite plate but the error is unacceptable for C12and C66.Vishnuvardhan et al.［51］reconstructed the nine elastic moduli of composite laminates using a singletransmitter-multiple-receiver configuration. The inversion algorithms，GAs，are used to minimize the error between the measured velocities and the calculated ones. However，the experimentally measured group velocities were considered equal to the phase velocities during the reconstruction process，in which case the inversion results might be incorrect.

At NUAA，a laser ultrasonic system is incorporated to increase the accuracy for Lamb wave velocity measurements and GAs are then applied to invert the measured velocities and conclude to the reconstructed material properties.

2.1 GAs for the inversion of phase velocities

The GA-based inverse method of reconstruction starts with a population of randomly generated candidate values for the stiffness coefficient set，each of that are then evaluated using the error function. A few of them with the lowest error function are preserved，while the others are replaced with new trials. GAs parameters are listed in Table 1.The error function to be minimized，during the inversion of the ultrasonic phase velocities，can be represented as follows

The objective function above describes the relative error between the experimental velocity veandthe theoretical ones vt. When the propagation direction is specified，n denotes the number of the chosen frequencies.

Table 1 GAs parameters used in the inversion procedure

A basic knowledge that the propagation direction of phase velocity is not same as that of the group velocity in an anisotropic plate［63］as shown in Fig.13，should be taken into account in GAs. The snapshot in Fig.13 is based on an FEM simulation of Lamb wave propagation at a given frequency in a unidirectional plate whose properties are listed in Table 2. The direction of group velocity is along the white dot line that implies the sensing line in experiments，while the relevant phase velocity travels perpendicular to the wave front. Because the measured velocity veis obtained from the group angle，the relevant theoretical phase velocity，which is perpendicular to the wave front，should be transformed to the group velocity coordinate by where θ is the angle between the phase velocity and the group velocity.

Fig.13 Schematic of different propagation angle between phase velocity and group velocity in a unidirectional plate (based on an FEM snapshot)

Table 2 Typical stiffness coefficients[65] for two CFRP plates with thickness of 1 mm and density of 1 600 kg/m3 GPa

2.1.1 Sensitivity analysis of phase velocities to the change of stiffness coefficients

in the study［64］，the authors reconstructed the material properties using a series of angles from 0°to 90°（considering the symmetric characteristic of the chosen plate）with an increment of 0.2°，which increases the complexity of the inversion problem by a significant margin. To simplify the reconstruction process and improve the efficiency，a sensitivity analysis should be performed to identify regions where the change of the material properties significantly affects the measured Lamb wave velocities.in addition，the relative sensitivity of unknown parameters（elastic moduli）becomes important in a multiparameter optimization problem as described in this section. if the sensitivities of the nine parameters are remarkably different from each other，then the performance of the GAs may not be reliable，that is，the error of the reconstructed elastic moduli will be larger for some of the elastic moduli，to which the velocities are less sensitive.

To simplify the inversion procedure and improve the reliability of the reconstructed multiple parameters，a sensitivity analysis of the phase velocity to the change of the nine stiffness coefficients along 0° to 90° is conducted at a specified frequency，as illustrated in Fig.14. The horizontal coordinate denotes the group velocity direction and the phase velocities have been transformed to the group velocity direction. The vertical coordinate denotes the relative difference of the phase velocities when each stiffness coefficient is increased by 20%.

During the sensitivity analysis，the material properties of the two plates are given in Table 2 that shows the typical stiffness coefficients of a unidirectional and a cross-ply CFRP plates. However，this does not mean that the experimental specimens are of the same material properties in Table 2.

According to Fig.14（a），for both plates，phase velocities of A0mode in 0° direction are most sensitive to C11and C55. Phase velocities of A0mode in 90° direction are most sensitive to C22and C44，while phase velocities of A0mode in the 20°—60° region are sensitive to C12and C66. However，A0mode velocities are much less sensitive to the left three parameters（C13，C23，and C33）and thus A0mode velocities cannot be used to reconstruct the left three stiffness coefficients efficiently.

Fig.14 The relative difference of phase velocities when varying each stiffness coefficient by 20%

As for S0mode，as shown in Fig.14（b），velocities in 0° direction are more sensitive to C11and C13，velocities in 90° direction are more sensitive to C22，C23，and C33. Velocities in 20° direction are most sensitive to C12and C66for the unidirectional plate，while velocities in 45° direction are most sensitive to C12and C66for the cross-ply one.

Although both A0and S0mode velocities in 0°direction are sensitive to C11，S0mode is chosen to reconstruct C11due to its higher objective function value（higher sensitivity）compared with that of A0mode. Similarly，S0mode velocities are more sensitive to C22，C12，and C66and will be inversed to reconstruct these three parameters.

By comparison and combination of A0mode and S0mode，the reconstructing method regarding the specific combination of modes and directions is illustrated in Tables 3，4. The C11and C22using S0mode velocities are reconstructed first and regarded as known when inverting C44and C55using A0mode velocities. To avoid non-uniqueness of solution for the multiparameter problem，phase velocities at a serious of frequencies（dispersion curves）will be utilized in GAs.

2.1.2 Reconstruction method

The reconstruction of the nine stiffness coefficients was carried out for two different CFRP plates. With MATLAB Optimization Tool，the material properties were identified by inverting the experimental velocities shown in Fig.15. According to the reconstruction method in Tables 3，4，C55and C44are taken from the reconstruction results in 0°and 90° direction of A0mode，respectively. C11and C13are identified using the S0mode velocities in 0°direction. Based on the S0mode velocities in 90° direction，C22，C23，and C33are calculated. As shown in Fig.15（b），phase velocities in the unidirectional plate along 20° direction are sensitive not only to C12and C66but also to others. To assure the accuracy of the reconstructed value of C12and C66，the other seven parameters should be identified first and substituted into GAs for the further reconstruction of C12and C66.

Fig.15 The experimental phase dispersion curves along specific directions

Table 3 Reconstruction method of stiffness coefficients in a unidirectional plate

Table 4 Reconstruction method of stiffness coefficients in a cross-ply plate

2.2 Results and discussion

2.2.1 Reconstructed results

The reconstructed results are listed in Tables 5，6. The relative phase velocities error along 0° and 90oof S0mode is less than 1% and thus the reconstructed results（C11，C13，C22，C23and C33）are of high accuracy according to the sensitivity analysis in Fig.14. The objective function value in 90° direction of A0mode is always the highest，which might result from the inaccuracy and the unsmooth dispersion curve of the experimental velocities in 90° direction of A0mode.

2.2.2 Verification of the reconstructed results

To verify the accuracy of the reconstructed results in Tables 5，6，the nine parameters are substituted into a commercial software ABAQUS，as the material properties，to simulate the propagation characteristics of Lamb waves. A normal-to-plane nodal force，whose time distribution is taken from the experimental signal sensed by the AE transducer when the laser source points at the position of the sensor，as shown in Fig.3，is set as the load condition in the FEM simulation. The model is pinned at the boundaries. The element is of liner solid type and the element size is set as 0.3 mm×0.3 mm.The maximum time step of the explicit scheme is set as 0.01 μs. The simulation results are compared with the experimental images，as shown in Fig.16.The wave curves of the simulated three basic Lamb wave modes coincide well with the experimental ones.

Meanwhile，the calculated and measured velocities curves of S0mode are compared in Fig.17. As illustrated in Figs.17（a），（c），the theoretical and experimental dispersion curves in 0° and 90° directions of S0mode coincide well. There exists difference between the theoretical and experimental dis-persion curves in other directions（20° for the unidirectional plate and 45° for the cross-ply one），which is mainly because that the reconstruction process of these directions is based on the 0° and 90° directions results and the previous error is brought in.

Table 5 Reconstructed parameters and the objective function of the unidirectional plate

Table 6 Reconstructed parameters and the objective function value of the cross-ply plate

Fig.16 Comparison of the snapshots of Lamb waves propagation (Unidirectional plate: (a) — (d); Cross-ply plate:(e)—(f))

information in only three directions is utilized in the inversion procedure and thus velocities in other directions need to be verified. For the unidirectional plate，the experimental velocities of S0mode（f = 500 kHz）in other seven directions（10°，30°，40°，50°，60°，70°，and 80°）were extracted from the laser ultrasonic results and were found to agree well with the theoretical ones，as shown in Fig.17（b）. in Fig.17（d），the measured velocities of S0mode（f = 500 kHz）in the cross-ply plate in other eight directions（10°，20°，30°，40°，50°，60°，70°，and 80°）are taken into account and also match the theoretical ones.

3 Conclusions

Fig.17 Comparison of the experimental and theoretical phase velocities of S0 mode based on the reconstructed results

Laser ultrasonic technique has shown great potentials for NDT applications in composite structures. The utilization of laser ultrasonic can help to ensure the safety of composite structures and prolong the service life. it is extremely useful in the background of composite materials being widely used in various industries. This paper concludes two main applications developed at NUAA that employs laser ultrasonic technique for damage detection and stiffness parameter reconstruction for composite laminates，which covers from the manufacture to various service conditions of composite laminates.For damage detection，an AiW-based damage visualization algorithm is proposed that not only captures the wave propagation changes caused by the impact damage，but also removes the negative influence from the reflected waves resulting in a muchimproved resolution for damage visualization. in terms of stiffness reconstruction，a thorough sensitivity analysis has been conducted to reduce the complexity for the inverse problem，which was solved using GAs，the utilization of laser ultrasonic technique improved the accuracy for velocity extraction leading to a better reconstruction result.

Based on current foundations， ongoing research aims to develop a comprehensive and compact laser inspection system that is capable to realize in-situ，noncontact and nondestructive inspections for both normal and curved composite structures，detecting macroscopic damages e.g. impact damages，major delamination etc.，as well as microscopic damages，e.g.fatigue cracks.