AIR-COUPLED ULTRASONIC PLANE STRESS DETECTION METHOD FOR COMPOSITE MATERIAL BASED ON DUAL-MODAL SOUND-TIME RATIOS OF LAMB WAVE

Abstract

An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave includes: S1: assembling a detection device; S2: based on the detection device and the dual-modal sound-time ratios of the Lamb wave in S1, using unidirectionally loaded stress to obtain different stress coefficient relations; S3: based on the stress coefficient relations in S2, solving stress coefficients; S4: based on the stress coefficients in S3, acquiring three sound-time ratios; S5: based on the sound-time ratios in S4, describing a stress state of a detection point; and S6: repeating S4 and S5 till completing detection and scanning. The method improves the accuracy of stress coefficient calibration and air-coupled ultrasonic stress representation of a composite material panel greatly.

Claims

1. An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, comprising the following steps: S1: assembling a detection device; S2: based on the detection device and the dual-modal sound-time ratios of the Lamb wave in S1, using unidirectionally loaded stress to obtain different stress coefficient relations; S3: based on the stress coefficient relations in S2, solving stress coefficients; S4: based on the stress coefficients in S3, acquiring three sound-time ratios; S5: based on the three sound-time ratios in S4, describing a stress state of a detection point; and S6: repeating S4 and S5 till completing detection and scanning.

2. The air-coupled ultrasonic plane stress detection method according to claim 1, wherein S1 is implemented as follows: the detection device comprises a function signal generator, a 50-ohm load, a low pass filter, a high speed acquisition board, an anti-aliasing filter, a pre-amplifier, an upper computer, an exciting air-coupled transducer, a receiving air-coupled transducer and a composite material panel; wherein a signal emitted by the function signal generator is emitted to the composite material panel through the 50-ohm load, the low pass filter and the exciting air-coupled transducer, the composite material panel reflects the received signal back to the receiving air-coupled transducer, and the receiving air-coupled transducer transmits the signal to the upper computer through the pre-amplifier, the anti-aliasing filter and the high speed acquisition board.

3. The air-coupled ultrasonic plane stress detection method according to claim 1, wherein the dual-modal sound-time ratios of the Lamb wave in S2 comprise: defining a sound velocity matrix [V.sub.ij]:
V.sub.ij=(V.sub.ij).sub.0+?V.sub.ij,?V.sub.ij=?.sub.ij?.sub.kl (1) wherein ?.sub.ij is a matrix representing an orthotropic stress coefficient, and a form of the matrix is similar to a form of a rigidity matrix: $\begin{matrix} [?_{ij}] = [\begin{matrix} ?_{11} & ?_{12} & ?_{13} \\ ?_{12} & ?_{22} & ?_{23} & 0 \\ ?_{13} & ?_{23} & ?_{33} \\ ?_{44} \\ 0 & ?_{55} \\ ?_{66} \end{matrix}] & (2) \end{matrix}$ for a stress in an XOY plane, ?.sub.zz=?.sub.xz=?.sub.yz=0; then: $\begin{matrix} [? V_{ij}] = [\begin{matrix} ? V_{X} \\ ? V_{Y} \\ ? V_{Z} \end{matrix}] = ? [\begin{matrix} ?_{11} ?_{xx} + ?_{12} ?_{yy} & ?_{66} ?_{xy} / 2 \\ ?_{66} ?_{xy} / 2 & ?_{12} ?_{xx} + ?_{22} ?_{yy} \\ ?_{13} ?_{xx} + ?_{23} ?_{yy} \end{matrix}] & (3) \end{matrix}$ a sound velocity change matrix [?V.sub.ij] has a vector decomposition characteristic similar to a vector decomposition characteristic of the stress, a principle value of [?V.sub.ij] is as follows: $\begin{matrix} ? V_{X, Y} = \frac{1}{2} {[(?_{11} + ?_{12}) ?_{xx} + (?_{12} + ?_{22}) ?_{yy}] ? \sqrt{{[(?_{11} - ?_{12}) ?_{xx} + (?_{12} - ?_{22}) ?_{yy}]}^{2} + {(?_{66} ?_{xy})}^{2}}} & (4) \end{matrix}$ $\begin{matrix} \tan 2 ? = \frac{?_{66} ?_{xy}}{(?_{11} - ?_{12}) ?_{xx} + (?_{12} - ?_{22}) ?_{yy}} & (5) \end{matrix}$ there is the following relation between a stress in any direction and a principal stress ?.sub.1 and ?.sub.2: $\begin{matrix} {\begin{matrix} ?_{xx} = \frac{1}{2} [?_{1} + ?_{2} + (?_{1} - ?_{2}) \cos 2 ?] \\ ?_{yy} = \frac{1}{2} [?_{1} + ?_{2} - (?_{1} - ?_{2}) \cos 2 ?] \\ ?_{xy} = \frac{1}{2} (?_{1} - ?_{2}) \sin 2 ? \end{matrix} & (6) \end{matrix}$ the equation (6) is substituted into the equation (4) to obtain sound velocity change ?V.sub.? in any direction in a state of plane stress as follows: $\begin{matrix} {? V_{?} = ? \frac{1}{2} [\frac{?_{1} + ?_{2}}{2} (?_{11} + ?_{22} + ? 2 ?_{12}) + \frac{?_{1} - ?_{2}}{2} (?_{11} - ?_{22}) \cos 2 ?] + ? ? \frac{1}{2} {{[\begin{matrix} \frac{?_{1} + ?_{2}}{2} (?_{11} - ?_{22}) + \\ \frac{?_{1} - ?_{2}}{2} (?_{11} + ?_{22} - 2 ?_{12}) \cos 2 ? \end{matrix}]}^{2} + ? ? {(\frac{?_{1} - ?_{2}}{2} ?_{66} \sin 2 ?)}^{2}}}^{\frac{1}{2}} \cos 2 (? - ?) & (7) \end{matrix}$ the equation (5) is substituted into the equation (7) to obtain a relation between a velocity change and the principal stress of the plane as follows: $\begin{matrix} ? V_{?} = \frac{?_{11} + ?_{22} + 2 ?_{12}}{4} (?_{1} + ?_{2}) + \frac{?_{11} - ?_{22}}{4} (?_{1} - ?_{2}) \cos 2 ? + \frac{?_{11} - ?_{22}}{4} (?_{1} + ?_{2}) \cos 2 ? + \frac{?_{11} + ?_{22} - 2 ?_{12}}{4} (?_{1} - ?_{2}) \cos 2 ?cos 2 ? + \frac{?_{66}}{4} (?_{1} - ?_{2}) \sin 2 ?sin2? & (8) \end{matrix}$ the following constants are introduced: $\begin{matrix} {\begin{matrix} n_{1} = (?_{11} - ?_{22}) / 4 \\ n_{2} = (?_{11} + ?_{22} + 2 ?_{12}) / 4 \\ n_{3} = (?_{11} + ?_{22} - 2 ?_{12}) / 4 \\ n_{4} = ?_{66} / 4 \end{matrix} & (9) \end{matrix}$ the equation (8) is transformed into:
?V.sub.?=K.sub.1?.sub.1+K.sub.2?.sub.2
K.sub.1=n.sub.1(cos 2?+cos 2?)+n.sub.2+n.sub.3 cos 2? cos 2?+n.sub.4 sin 2? sin 2?
K.sub.2=?n.sub.1(cos 2??cos 2?)+n.sub.2?n.sub.3 cos 2? cos 2??n.sub.4 sin 2? sin 2?(10) wherein K.sub.1 and K.sub.2 are proportionality coefficients between the velocity change and the principal stress ?.sub.1 and ?.sub.2 of the plane, respectively.

4. The air-coupled ultrasonic plane stress detection method according to claim 3, wherein the Lamb wave is propagated in the composite material panel, and phase velocities of a symmetric mode and an asymmetric mode of the Lamb wave are represented as:
V.sub.S0?=K.sub.1?.sub.1+K.sub.2?.sub.2+V.sub.S0
V.sub.A0?=K.sub.1?.sub.1+K.sub.2?.sub.2+V.sub.A0 (11) wherein V.sub.S0? and V.sub.A0? are the phase velocities of the symmetric mode and the asymmetric mode in the stress state, and V.sub.S0 and V.sub.A0 are initial phase velocities of the symmetric mode and the asymmetric mode in an unstress state.

5. The air-coupled ultrasonic plane stress detection method according to claim 4, wherein a ratio of a propagation time of the asymmetric mode to a propagation time of the symmetric mode is taken as follows: $\begin{matrix} R = \frac{T_{A 0 ?}}{T_{S 0 ?}} = \frac{V_{S 0 ?}}{V_{A 0 ?}} = \frac{K_{1} ?_{1} + K_{2} ?_{2} + V_{A 0}}{K_{1} ?_{1} + K_{2} ?_{2} + V_{S 0}} = 1 + \frac{V_{A 0} - V_{S 0}}{K_{1} ?_{1} + K_{2} ?_{2} + V_{S 0}} & (12) \end{matrix}$ the equation (12) is a relation between the ratio of the propagation time of the symmetric mode to the propagation time of the asymmetric mode of the Lamb wave and the principal stress in a plane problem of the composite material; there are three unknown quantities ?.sub.1, ?.sub.2 and ? for an orthotropic panel in the stress state of the plane; wherein the ratios R.sub.?1, R.sub.?2 and R.sub.?3 of the propagation time of the symmetric mode to the propagation time of the asymmetric mode of the Lamb wave in three different directions ?.sub.1, ?.sub.2 and ?.sub.3 are measured, and the ratios are solved through simultaneous equations.

6. The air-coupled ultrasonic plane stress detection method according to claim 5, wherein the unidirectionally loaded stress is ?.sub.1=? and ?.sub.2=0, and the equation (12) is transformed into: $\begin{matrix} R = \frac{T_{A 0 ?}}{T_{S 0 ?}} = \frac{V_{S 0 ?}}{V_{A 0 ?}} = \frac{K_{1} ? + V_{A 0}}{K_{1} ? + V_{S 0}} = 1 + \frac{V_{A 0} - V_{S 0}}{K_{1} ? + V_{S 0}} & (17) \end{matrix}$ $K_{1} = n_{1} (\cos 2 ? + \cos 2 ?) + n_{2} + n_{3} \cos 2 ?cos2? + n_{4} \sin 2 ?sin2?$ ? in the equation (17) is a first included angle between a fiber direction and the stress, ? is a second included angle between a detection direction and the fiber direction, and the first included angle and the second included angle are combined differently to obtain different stress coefficient relations.

7. The air-coupled ultrasonic plane stress detection method according to claim 6, wherein the step of solving the stress coefficients in S3 is implemented as follows: in a case where ?=?=0?,
R=(2n.sub.1+n.sub.2+n.sub.3)?(18) in a case where ?=0?, ?=45?
R=(n.sub.1+n.sub.2)?(19) in a case where ?=?=45?,
R=(n.sub.2+n.sub.4)?(20) in a case where ?=45?, ?=90? $\begin{matrix} R^{} = (n_{2} - n_{1}) ? & (21) \end{matrix}$ $wherein$ $\begin{matrix} R^{} = \frac{V_{A 0} - V_{S 0}}{R - 1} - V_{S 0} & (22) \end{matrix}$ the equations (18)-(22) are simultaneous to obtain n.sub.1, n.sub.2, n.sub.3 and n.sub.4.

8. The air-coupled ultrasonic plane stress detection method according to claim 7, wherein the step of acquiring the three sound-time ratios in S4 comprises: placing an exciting air-coupled transducer and a receiving air-coupled transducer at a dip angle ? on one side of a composite material panel test sample, allowing an included angle between a connection direction of the exciting air-coupled transducer and the receiving air-coupled transducer and a principal direction of a fiber to be ?=0?, and setting a distance between the exciting air-coupled transducer and the receiving air-coupled transducer to be L; performing stretching by using unidirectional stress and superposing the stretching direction with the fiber direction of the composite material panel, wherein the included angle ?=0?; determining an exciting signal cycle as N, determining a center frequency f of the exciting air-coupled transducer and the receiving air-coupled transducer, and extracting and receiving sounds and times of the asymmetric mode A0 and the symmetric mode S0 in an echo signal and substituting the sounds and times into the equation (17) to obtain the sound-time ratios; applying stretching and compressing stress with a certain step length to a composite material panel standard component by stretching and compressing experiments to obtain the sound-time ratios R under different stresses, respectively; in order to eliminate accidental errors, repeating the stretching and compressing experiments for 10 times under a same stress to obtain 10 sound-time ratios, and averaging the 10 sound-time ratios as the sound-time ratio under the stress, and substituting the sound-time ratio into the equation (22) to obtain R; acquiring N groups of data (R.sub.i,?.sub.i), i=1, 2, 3 . . . N, and obtaining a relation expression between the stress and the sound-time ratio through least-squares algorithm linear fitting,
R=k?(23) wherein k is a slope of a straight line obtained by least-squares algorithm linear fitting; and selecting ?=0?, ?=45?; ?=?=45?; ?=45?, ?=90? respectively and repeating the above process for three times to obtain the equations (18)-(21) that are simultaneous to obtain n.sub.1, n.sub.2, n.sub.3 and n.sub.4.

9. The air-coupled ultrasonic plane stress detection method according to claim 8, wherein based on the three sound-time ratios in S4, the step of describing the stress state of the detection point in S5 comprises: performing step-scanning on a to-be-measured region of the composite material panel test sample in x and y directions with ?x and ?y as step pitches, respectively, wherein for each stress detection point P, the included angle ? between the principal stress ?.sub.1 and ?.sub.2 and the fiber direction and the stress direction are unknown, selecting ?.sub.1=?45?, ?.sub.2=0? and ?.sub.3=45? by using the S4 to perform experiments for three times, acquiring the stress state of the point P through the equations (12)-(16) in combination with the stress coefficients n.sub.1, n.sub.2, n.sub.3 and n.sub.4 acquired by calibration, repeating the process, and after completing the scanning, acquiring the stress state of the entire composite material test sample.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0049] FIG. 1 is a schematic diagram of an analytical coordinate system of the present invention.

[0050] FIG. 2 is a schematic diagram of an experimental device of the present invention.

[0051] FIG. 3 is a flow diagram of a method of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0052] The technical solutions in the embodiments of the present invention will be clearly and completely described below in combination with the accompanying drawings in the embodiments of the present invention. Apparently the embodiments described are merely some rather than all of the embodiments of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art without making creative efforts fall within the protection scope of the present invention.

[0053] The center frequency of the air-coupled transducer is determined as f and the dip angle thereof is ? according to a frequency dispersion curve and a Snell law of the Lamb wave. The included angle between a connection direction of the exciting air-coupled transducer and the receiving air-coupled transducer and the principal direction of the fiber is allowed to be ?. The air-coupled transducer is placed as required on one side of the composite material panel, and the propagation distance of the acoustic wave between the exciting air-coupled transducer and the receiving air-coupled transducer is set to be an appropriate value. In order to guarantee an acoustic beam with enough energy, the exciting signal cycle is determined as N, a signal generator is connected to a power amplifier to generate an exciting signal, and subjected to impedance matching of the 50-ohm load and the low pass filter, the exciting signal is applied to the exciting air-coupled transducer. The echo is received by the receiving air-coupled transducer in a receiving position, is amplified and filtered by the pre-amplifier and the anti-aliasing filter, and is collected by the high speed acquisition board and uploaded to the upper computer. The sounds and times of the zero-order symmetrical mode (S0) and the zero-order asymmetric model (A0) are extracted to acquire the sound-time ratio R.sub.?. The included angle ? between a connection direction of the exciting air-coupled transducer and the receiving air-coupled transducer and the principal direction of the fiber of the material is changed, the sound-time ratios R.sub.?1, R.sub.?2 and R.sub.?3 are measured and extracted three times, and the sound-time ratios are simultaneous to obtain the stress state in the corresponding position (x,y). Next, ?x and ?y steps can be performed along the scanning direction to continue the above process till scanning is completed, so as to obtain the stress state of the entire composite material panel.

[0054] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave specifically includes the following steps: [0055] S1: a detection device is assembled; [0056] S2: based on the detection device and the dual-modal sound-time ratios of the Lamb wave in S1, unidirectionally loaded stress is used to obtain different stress coefficient relations; [0057] S3: based on the stress coefficient relations in S2, stress coefficients are solved; [0058] S4: based on the stress coefficients in S3, three sound-time ratios are acquired; [0059] S5: based on the sound-time ratios in S4, a stress state of a detection point is described; and [0060] S6: S4 and S5 are repeated till detection and scanning are completed.

[0061] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where S1 is specifically as follows: the detection device includes a function signal generator, a 50-ohm load, a low pass filter, a high speed acquisition board, an anti-aliasing filter, a pre-amplifier, an upper computer, an exciting air-coupled transducer 1, a receiving air-coupled transducer 2 and a composite material panel 3; [0062] a signal emitted by the function signal generator is emitted to the composite material panel 3 through the 50-ohm load, the low pass filter and the exciting air-coupled transducer 1, the composite material panel 3 reflects the received signal back to the receiving air-coupled transducer 2, and the receiving air-coupled transducer 2 transmits the signal to the upper computer through the pre-amplifier, the anti-aliasing filter and the high speed acquisition board.

[0063] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where the dual-modal sound-time ratios of the Lamb wave in S2 are specifically as follows: the technical principle of the dual-modal sound-time ratios of the Lamb wave are specifically described with reference to FIG. 1. The analytical coordinate system shown in FIG. 1 is established to analyze the propagation law of the Lamb wave in the composite material panel. Three coordinate systems are used in FIG. 1, which are respectively (1) the coordinate system XOY (the fiber coordinate system, usually known) of the principal direction of the material; (2) the coordinate system ?.sub.1O?.sub.2 of the principal stress (known for detection); and (3) the measurement coordinate system XiOYi (actual detection direction, known). With respect to the composite material panel, the paving direction of the fiber is often easily accessible. Therefore, during ultrasonic stress measurement, the fiber direction is taken as angular reference all the time.

[0064] In an orthotropic medium, the ultrasonic wave has different propagation velocities in different directions. The ultrasonic wave propagated in the principal direction of the material has an extreme value of the propagation velocity, that is, the ultrasonic wave corresponds to the principal direction of the ultrasonic wave. In the unstress state, it is assumed that the ultrasonic wave in the XOY plane has two principal propagation directions X.sub.0 and Y.sub.0. At the time, it is in a natural state. When the material bears the principal stress of the plane forming an angle ? with the fiber direction, it corresponds to a deformation state. Relative to the unstress state, the propagation sound velocity and the principal propagation direction affected by the stress will change. It is assumed that the coordinate system of the principal propagation direction at the time is XOY, and the coordinate system forms an angle ? with the coordinate system of the material direction. It can be known from an elastic wave equation in the anisotropic medium that there is an approximately linear relation between the velocity change ?V.sub.ij(i,j=1,2,3) and the ultrasonic wave and the stress ? ?.sub.kl(k,l=1,2,3). For the convenience of analysis, the sound velocity matrix [V.sub.ij] is defined:

V.sub.ij=(V.sub.ij).sub.0+?V.sub.ij,?V.sub.ij=?.sub.ij?.sub.kl (1) [0065] where ?.sub.ij is a matrix representing an orthotropic stress coefficient, and the form of the matrix is similar to that of a rigidity matrix:

[00011] $\begin{matrix} [?_{ij}] = [\begin{matrix} ?_{11} & ?_{12} & ?_{13} \\ ?_{12} & ?_{22} & ?_{23} & 0 \\ ?_{13} & ?_{23} & ?_{33} \\ ?_{44} \\ 0 & ?_{55} \\ ?_{66} \end{matrix}] & (2) \end{matrix}$

[0066] For the stress in an XOY plane, ?.sub.zz=?.sub.xz=?.sub.yz=0;

[0067] Then:

[00012] $\begin{matrix} [? V_{ij}] = [\begin{matrix} ? V_{X} \\ ? V_{Y} \\ ? V_{Z} \end{matrix}] = ? [\begin{matrix} ?_{11} ?_{xx} + ?_{12} ?_{yy} & ?_{66} ?_{xy} / 2 \\ ?_{66} ?_{xy} / 2 & ?_{12} ?_{xx} + ?_{22} ?_{yy} \\ ?_{13} ?_{xx} + ?_{23} ?_{yy} \end{matrix}] & (3) \end{matrix}$

[0068] A sound velocity change matrix [?V.sub.ij] has a vector decomposition characteristic similar to that of the stress, a principle value of [?V.sub.ij] is as follows:

[00013] $\begin{matrix} ? V_{X, Y} = \frac{1}{2} {[(?_{1 1} + ?_{1 2}) ?_{xx} + (?_{1 2} + ?_{2 2}) ?_{yy}] ? \sqrt{{[(?_{1 1} - ?_{1 2}) ?_{x x} + (?_{1 2} - ?_{2 2}) ?_{yy}]}^{2} + {(?_{6 6} ?_{x y})}^{2}}} & (4) \end{matrix}$ $\begin{matrix} \tan 2 ? = \frac{?_{6 6} ?_{x y}}{(?_{1 1} - ?_{1 2}) ?_{x x} + (?_{1 2} - ?_{2 2}) ?_{yy}} & (5) \end{matrix}$ [0069] there is the following relation between the stress in any direction and the principal stress ?.sub.1 and ?.sub.2:

[00014] $\begin{matrix} {\begin{matrix} ?_{x x} = \frac{1}{2} [?_{1} + ?_{2} + (?_{1} - ?_{2}) \cos 2 ?] \\ ?_{yy} = \frac{1}{2} [?_{1} + ?_{2} - (?_{1} - ?_{2}) \cos 2 ?] \\ ?_{x y} = \frac{1}{2} (?_{1} - ?_{2}) \sin 2 ? \end{matrix} & (6) \end{matrix}$ [0070] the equation (6) is substituted into the equation (4) to obtain sound velocity change ?V.sub.? in any direction in the state of plane stress as follows:

[00015] $\begin{matrix} ? V_{?} = \frac{1}{2} [\frac{?_{1} + ?_{2}}{2} (?_{11} + ?_{22} + 2 ?_{12}) + \frac{?_{1} - ?_{2}}{2} (?_{11} - ?_{22}) \cos 2 ?] + \frac{1}{2} {{[\begin{matrix} \frac{?_{1} + ?_{2}}{2} (?_{11} - ?_{22}) + \\ \frac{?_{1} - ?_{2}}{2} (?_{11} + ?_{22} - 2 ?_{12}) \cos 2 ? \end{matrix}]}^{2} + {(\frac{?_{1} - ?_{2}}{2} ?_{66} \sin 2 ?)}^{2}}^{\frac{1}{2}} \cos 2 (? - ?) & (7) \end{matrix}$ [0071] the equation (5) is substituted into the equation (7) to obtain a relation between the velocity change and the principle stress of the plane as follows:

[00016] $\begin{matrix} ? V_{?} = \frac{?_{1 1} + ?_{2 2} + 2 ?_{1 2}}{4} (?_{1} + ?_{2}) + \frac{?_{1 1} - ?_{2 2}}{4} (?_{1} - ?_{2}) \cos 2 ? + \frac{?_{1 1} - ?_{2 2}}{4} (?_{1} + ?_{2}) \cos 2 ? + \frac{?_{1 1} + ?_{2 2} - 2 ?_{1 2}}{4} (?_{1} - ?_{2}) \cos 2 ?cos2? + \frac{?_{6 6}}{4} (?_{1} - ?_{2}) \sin 2 ?sin2? & (8) \end{matrix}$ [0072] the following constants are introduced:

[00017] $\begin{matrix} {\begin{matrix} n_{1} = (?_{1 1} - ?_{2 2}) / 4 \\ n_{2} = (?_{1 1} + ?_{2 2} + 2 ?_{1 2}) / 4 \\ n_{3} = (?_{1 1} + ?_{2 2} - 2 ?_{1 2}) / 4 \\ n_{4} = ?_{6 6} / 4 \end{matrix} & (9) \end{matrix}$ [0073] the equation (8) can be transformed into:

?V.sub.?=K.sub.1?.sub.1+K.sub.2?.sub.2

K.sub.1=n.sub.1(cos 2?+cos 2?)+n.sub.2+n.sub.3 cos 2? cos 2?+n.sub.4 sin 2? sin 2?

K.sub.2=?n.sub.1(cos 2??cos 2?)+n.sub.2?n.sub.3 cos 2? cos 2??n.sub.4 sin 2? sin 2?(10) [0074] where K.sub.1 and K.sub.2 are proportionality coefficients between velocity change and the principal stress ?.sub.1 and ?.sub.2 of the plane, respectively.

[0075] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where the Lamb wave is propagated in the composite material panel the phase velocities of a symmetric mode and an asymmetric mode of the Lamb wave can be represented as:

V.sub.S0?=K.sub.1?.sub.1+K.sub.2?.sub.2+V.sub.S0

V.sub.A0?=K.sub.1?.sub.1+K.sub.2?.sub.2+V.sub.A0 (11) [0076] where V.sub.S0? and V.sub.A0? are phase velocities of the symmetric mode and the asymmetric mode in the stress state, and V.sub.S0 and V.sub.A0 are initial phase velocities of the symmetric mode and the asymmetric mode in the unstress state.

[0077] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where a ratio of propagation time of the asymmetric mode to propagation time of the symmetric mode is taken as follows:

[00018] $\begin{matrix} R = \frac{T_{A 0 ?}}{T_{S 0 ?}} = \frac{V_{S 0 ?}}{V_{A 0 ?}} = \frac{K_{1} ?_{1} + K_{2} ?_{2} + V_{A 0}}{K_{1} ?_{1} + K_{2} ?_{2} + V_{S 0}} = 1 + \frac{V_{A 0} - V_{S 0}}{K_{1} ?_{1} + K_{2} ?_{2} + V_{S 0}} & (12) \end{matrix}$ [0078] the equation (12) is the relation between the ratio of the propagation time of the symmetric mode to the propagation time of the asymmetric mode of the Lamb wave and the principal stress in the plane problem of the composite material. There are three unknown quantities ?.sub.1, ?.sub.2 and ? for the orthotropic panel in the stress state of the plane. Therefore, it needs to measure the ratios R.sub.?1, R.sub.?2 and R.sub.?3 of the propagation time of the symmetric mode to the propagation time of the asymmetric mode of the Lamb wave in three different directions ?.sub.1, ?.sub.2 and ?.sub.3, and the ratios are solved through simultaneous equations.

[0079] To simplify calculation, measurement angles ?.sub.1=?45?, ?.sub.2=0? and ?.sub.3=45? are taken to acquire the ratios of the propagation times of the symmetric mode and the asymmetric mode of the corresponding Lamb wave, R.sub.1, R.sub.2 and R.sub.3:

[00019] $\begin{matrix} {\begin{matrix} R_{1} = \frac{V_{A 0} - V_{S 0}}{R_{? 1} - 1} - V_{s 0} \\ R_{2} = \frac{V_{A 0} - V_{S 0}}{R_{? 2} - 1} - V_{s 0} \\ R_{3} = \frac{V_{A 0} - V_{S 0}}{R_{? 3} - 1} - V_{s 0} \end{matrix} & (13) \end{matrix}$ [0080] for the convenience of calculation and expression, the following variable parameters are introduced:

[00020] $\begin{matrix} {\begin{matrix} L = ?_{1} + ?_{2} \\ M = \sin 2 ? (?_{1} - ?_{2}) \\ N = \cos 2 ? (?_{1} - ?_{2}) \end{matrix} & (14) \end{matrix}$ [0081] a calculation expression of the amplitude and direction of the principal stress can be obtained by solving the equation (14):

[00021] $\begin{matrix} {\begin{matrix} ?_{1} = \frac{1}{2} (\sqrt{M^{2} + N^{2}} + L) \\ ?_{2} = \frac{1}{2} (L - \sqrt{M^{2} + N^{2}}) \\ ? = \frac{1}{2} \arctan (\frac{M}{N}) \end{matrix} & (15) \end{matrix}$ [0082] the three equations in the three directions are simultaneous to obtain:

[00022] $\begin{matrix} {\begin{matrix} L = \frac{(n_{1} + n_{3}) (R_{1} + R_{3}) - 2 n_{1} R_{2}}{2 (n_{2} n_{3} - {n_{1}}^{2})} \\ N = \frac{(n_{1} + n_{2} (R_{1} + R_{3}) - 2 n_{2} R_{2}}{2 ({n_{1}}^{2} - n_{2} n_{3})} \\ M = \frac{R_{3} - R_{1}}{2 n_{4}} \end{matrix} & (16) \end{matrix}$ [0083] the equation (16) is substituted into the equation (15) to obtain the amplitude and direction of the principal stress in the plane of the measured composite material, so as to acquire the stress state in the plane in any direction.

[0084] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where the unidirectionally loaded stress is ?.sub.1=? and ?.sub.2=0, and the equation (12) can be transformed into:

[00023] $\begin{matrix} R = \frac{T_{A 0 ?}}{T_{S 0 ?}} = \frac{V_{S 0 ?}}{V_{A 0 ?}} = \frac{K_{1} ? + V_{A 0}}{K_{1} ? + V_{S 0}} = 1 + \frac{V_{A 0} - V_{S 0}}{K_{1} ? + V_{S 0}} & (17) \end{matrix}$ $K_{1} = n_{1} (\cos 2 ? + \cos 2 ?) + n_{2} + n_{3} \cos 2 ?cos2? + n_{4} \sin 2 ?sin2?$ [0085] ? in the equation (17) is an included angle between a fiber direction and the stress, ? is an included angle between a detection direction and the fiber direction, and the two angles are combined differently to obtain different stress coefficient relations.

[0086] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where the solving the stress coefficients in S3 is specifically as follows: the following angle combination calibration experiment is designed in the present invention to determine the ultrasonic stress coefficients: [0087] in a case where ?=?=0?,

R=(2n.sub.1+n.sub.2+n.sub.3)?(18) [0088] in a case where ?=0?, ?=45?

R=(n.sub.1+n.sub.2)?(19) [0089] in a case where ?=?=45?,

R=(n.sub.2+n.sub.4)?(20) [0090] in a case where ?=45?, ?=90?

[00024] $\begin{matrix} R^{} = (n_{2} - n_{1}) ? & (21) \end{matrix}$ $where$ $\begin{matrix} R^{} = \frac{V_{A 0} - V_{S 0}}{R - 1} - V_{S 0} & (22) \end{matrix}$

[0091] With respect to the composite material panel, V.sub.S0 and V.sub.A0 are easily accessible, so that the equations (18)-(22) are simultaneous to obtain n.sub.1, n.sub.2, n.sub.3 and n.sub.4.

[0092] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where the acquiring the three sound-time ratios in S4 is specifically as follows: the specific calibration process is described below by taking ?=?=0? as an example:

[0093] The Lamb wave has the symmetric mode and the asymmetric mode and dispersion characteristics. Multi-order symmetric modes (S.sub.0, S.sub.1 . . . S.sub.i) and asymmetric modes (A.sub.0, A.sub.1 . . . A.sub.i) can be excited at the same exciting frequency. In order to allow the air-coupled transducer to excite relatively pure modes in the composite material panel, the exciting frequency f is determined according to the frequency dispersion curve of the Lamb wave in combination with the thickness of the to-be-detected component. After the product of frequency and thickness (frequency?thickness of the composite material panel) is determined, the phase velocity of the asymmetric mode A.sub.0 can be known, so that the dip angle ? of the exciting and receiving air-coupled transducers are determined according to the first critical angle of refraction of the Snell law in combination with the propagation velocity in air. The exciting air-coupled transducer is placed at a dip angle ? on one side of a composite material panel test sample, the included angle between the connection direction of the exciting air-coupled transducer and the receiving air-coupled transducer and the principal direction of the fiber is allowed to be ?=0?, and the distance between the receiving air-coupled transducer and the exciting air-coupled transducer is set to be L. Stretching is performed by using unidirectional stress and the stretching direction is superposed with the fiber direction of the composite material panel, i.e., the included angle ?=0?; in order to guarantee the acoustic beam with enough energy, the exciting signal cycle is determined as N, the signal generator is connected to the power amplifier to generate an exciting signal according to the above determined center frequency f of the air-coupled transducer, and subjected to impedance matching of the 50-ohm load and the low pass filter, the exciting signal is applied to the exciting air-coupled transducer. The echo is received by the receiving air-coupled transducer in the receiving position, is amplified and filtered by the pre-amplifier and the anti-aliasing filter, and is collected by the high speed acquisition board and uploaded to the upper computer. Sounds and times of the asymmetric mode A.sub.0 and the symmetric mode S.sub.0 in an echo signal are extracted and received and the sounds and times are substituted into the equation (17) to obtain the sound-time ratios.

[0094] Stretching and compressing stress (it is stipulated that the tensile stress is negative and the compressive stress is positive) is applied with a certain step length to a composite material panel standard component by means of stretching and compressing experiments to obtain the sound-time ratios R under different stress, respectively, in order to eliminate accidental errors, the experiments are repeated for 10 times under the same stress to obtain 10 sound-time ratios, and the sound-time ratios averaged as the sound-time ratio under the stress, and the sound-time ratio is substituted into the equation (22) to obtain R; N groups of data (R.sub.i,?.sub.i) are acquired, i=1, 2, 3 . . . N, and a relation expression between the stress and the sound-time ratio is obtained through least-squares algorithm linear fitting,

R=k?(23) [0095] where k is a slope of a straight line obtained by least-squares algorithm linear fitting; and [0096] ?=0?, ?=45?; ?=?=45?; ?=45?, ?=90? are selected respectively, and the above process is repeated three times to obtain the equations (18)-(21) that are simultaneous to obtain n.sub.1, n.sub.2, n.sub.3 and n.sub.4.

[0097] An air-coupled ultrasonic plane stress detection method for a composite material based on dual-modal sound-time ratios of a Lamb wave, where based on the sound-time ratios in S4, the describing the stress state of the detection point in S5 specifically includes: step-scanning is performed on a to-be-measured region of the composite material panel test sample in x and y directions with ?x and ?y as step pitches, respectively, where for each stress detection point P, the included angle ? between the principal stress ?.sub.1 and ?.sub.2 and the fiber direction and the stress direction are unknown, selecting ?.sub.1=?45?, ?.sub.2=0? and ?.sub.3=45? by using the S4 for experiments three times, acquiring the stress state of the point P through the equations (12)-(16) in combination with the stress coefficients n.sub.1, n.sub.2, n.sub.3 and n.sub.4 acquired by calibration, repeating the process, and after completing the scanning, acquiring the stress state of the entire composite material test sample.

AIR-COUPLED ULTRASONIC PLANE STRESS DETECTION METHOD FOR COMPOSITE MATERIAL BASED ON DUAL-MODAL SOUND-TIME RATIOS OF LAMB WAVE

Assignee

Inventors

Cpc classification

Classification Explorer

G01L1/255

PHYSICS

Classification Explorer

G01N29/07

PHYSICS

Classification Explorer

G01N29/041

PHYSICS

International classification

Classification Explorer

G01N29/04

PHYSICS

Classification Explorer

G01N29/07

PHYSICS

Classification Explorer

G01L1/25

PHYSICS

Abstract

Claims

Description