Method, System, Device and Medium for Online Monitoring of Plane Stress Field without Baseline Data Based on Piezoelectric Transducer Array

Abstract

The disclosure provides a method, system, device and medium for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array. Since Lamb waves have complex multi-mode characteristics, a suitable excitation frequency needs to be selected according to geometric dimensions of the structure to be measured, and then, only low-order mode Lamb waves are excited inside the measured structure to avoid serious waveform aliasing. For isotropic measured objects, anisotropic characteristics will be generated under the action of pre-stresses, that is, the propagation velocities of ultrasonic waves in all directions are different, but there is a linear relationship between velocity changes in different propagation directions and stresses. Therefore, there is still a linear relationship between the difference of velocity changes in different propagation directions and the stress. According to this characteristic, a characterization method of an absolute stress field without baseline data can be implemented. The method of the disclosure can make full use of the low attenuation characteristics of the Lamb waves to realize online monitoring of the plane stress field with a large coverage area.

Claims

1. A method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array, comprising: according to geometric dimensions and material parameters of a measured object, establishing a dispersion curve of Lamb waves, determining a cut-off frequency of a first-order Lamb wave mode, then selecting a suitable frequency within the cut-off frequency to excite piezoelectric wafers, and generating pure S0-mode Lamb waves by means of symmetrical excitation; according to the invariant property of acoustoelastic coefficients of isotropic materials in different coordinate systems, obtaining acoustoelastic equations in different propagation directions, then calculating a difference to obtain a linear relationship between the difference of velocity changes in different directions and the stress, and at the same time, according to a linear superposition relationship between stresses in different directions, calculating the magnitude of a plane stress field of the measured object through the difference of three groups of velocity changes in different directions, thereby finally realizing the accurate characterization of the plane stress field by means of measurement results of the piezoelectric transducer array in different directions; and exciting and receiving Lamb wave signals inside the measured object, performing Hilbert transform on the received signals, extracting signal envelopes, then calculating acoustic time differences in different propagation directions through a cross-correlation algorithm, according to a linear relationship between acoustic time differences and velocity differences, determining that there is still a linear relationship between stresses and acoustic time differences in different propagation directions, and then, substituting three groups of acoustic time differences in different directions into the pre-calibrated acoustoelastic equations to solve the plane stress field of the measured object.

2. The method according to claim 1, wherein according to the geometric dimensions and the material parameters of the measured object, the dispersion curve of Lamb waves is established, and a Lamb wave dispersion equation of a free plate is a Rayleigh-Lamb wave dispersion equation, which meets: $\begin{array}{cc}\frac{\tan \left(\mathrm{qh}\right)}{\tan \left(\mathrm{ph}\right)}=-\frac{4k^2\mathrm{pq}}{{\left(q^2-k^2\right)}^2},& \left(1\right)\end{array}$ $\begin{array}{cc}\frac{\tan \left(\mathrm{qh}\right)}{\tan \left(\mathrm{ph}\right)}=-\frac{{\left(q^2-k^2\right)}^2}{4k^2\mathrm{pq}},& \left(2\right)\end{array}$ wherein p and q are respectively expressed as: $\begin{array}{cc}p=\sqrt{\frac{{\omega }^2}{c_L^2}-k^2},& \left(3\right)\end{array}$ $\begin{array}{cc}q=\sqrt{\frac{{\omega }^2}{c_T^2}-k^2},& \left(4\right)\end{array}$ in the equation, c.sub.L and c.sub.T respectively represent velocities of longitudinal waves and transverse waves, h represents a half of a thickness of the plate, ω represents an angular frequency of ultrasonic waves, and k represents a wave number.

3. The method according to claim 2, wherein modulation signals with a frequency of 500 kHz are selected to excite the Lamb waves.

4. The method according to claim 3, wherein assuming that the measured object is subjected to plane stresses, in a coordinate system x-o-y, normal stresses are σ.sub.11 and σ.sub.22, and shear stresses are σ.sub.12 and σ.sub.21; according to the symmetry of the stresses, the values of σ.sub.12 and σ.sub.21 are equal, therefore, for Lamb waves propagating in an x direction, the acoustoelastic effect of the plane stress field is expressed as:
Δv=K.sub.11σ.sub.11+K.sub.12σ.sub.12+K.sub.22σ.sub.22  (5), in the equation, K.sub.11, K.sub.12 and K.sub.22 respectively represent acoustoelastic coefficients of the Lamb waves in different propagation directions; when the measured object belongs to an isotropic material, the acoustoelastic coefficients in different coordinate systems are exactly the same, therefore, the acoustoelastic effect of the Lamb waves propagating in an x′ direction is expressed as:
Δv=K.sub.11σ′.sub.11+K.sub.12σ′.sub.12+K.sub.22σ′.sub.22  (6), in the equation, σ′.sub.11, σ′.sub.12 and σ′.sub.22 represent normal stresses and shear stresses in a new coordinate system x′-o-y′; according to a coordinate transformation relationship between stresses, the following is obtained: $\begin{array}{cc}{\sigma }_{11}^{\prime }=\frac{{\sigma }_{11}+{\sigma }_{22}}{2}+\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\cos \left(2\theta \right)+{\sigma }_{12}\sin \left(2\theta \right),& \left(7\right)\end{array}$ $\begin{array}{cc}{\sigma }_{22}^{\prime }=\frac{{\sigma }_{11}+{\sigma }_{22}}{2}-\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\cos \left(2\theta \right)-{\sigma }_{12}\sin \left(2\theta \right),& \left(8\right)\end{array}$ $\begin{array}{cc}{\sigma }_{12}^{\prime }=-\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\sin \left(2\theta \right)+{\sigma }_{12}\cos \left(2\theta \right),& \left(9\right)\end{array}$ wherein θ represents an angle; a difference between propagation velocities along angles θ.sub.1 and θ.sub.2 is: $\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =K_{11}[\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\\ +{\sigma }_{12}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)]\\ +K_{22}[\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\left(-\cos \left(2{\theta }_1\right)+\cos \left(2{\theta }_2\right)\right)\\ +{\sigma }_{12}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)]\\ +K_{12}[\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)\\ +{\sigma }_{12}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)]\end{array};& \left(10\right)\end{array}$ in Equation (10) the influence of each component of the plane stress field on dv has a linear superposition relationship, and when the plane stress field has only one component, dv is respectively simplified as: $\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =\left(K_{11}-K_{22}\right)\left[\frac{{\sigma }_{11}}{2}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{11}}{2}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\end{array},& \left(11\right)\end{array}$ $\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =\left(K_{11}-K_{22}\right)\left[\frac{{\sigma }_{22}}{2}\left(-\cos \left(2{\theta }_1\right)+\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{22}}{2}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)\right]\end{array},& \left(12\right)\end{array}$ $\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =\left(K_{11}-K_{22}\right)\left[{\sigma }_{12}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\\ +K_{12}\left[{\sigma }_{12}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\end{array};& \left(13\right)\end{array}$ it is found that to complete the measurement of the plane stress field, it is only necessary to determine the values of K.sub.11−K.sub.22 and K.sub.12 through calibration tests; assuming that there is only σ.sub.11, θ.sub.1 is 0° and θ.sub.2 is 90°, the following is obtained:
dv=Δv(θ.sub.1)−Δv(θ.sub.2)=(K.sub.11−K.sub.22)σ.sub.11  (14); assuming that θ.sub.1 is 0° and θ.sub.2 is 45°, when there is σ.sub.11, dv is expressed as: $\begin{array}{cc}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)=\left(K_{11}-K_{22}\right)\frac{{\sigma }_{11}}{2}+K_{12}\frac{{\sigma }_{11}}{2};& \left(15\right)\end{array}$ two acoustoelastic values of K.sub.11−K.sub.22 and K.sub.12 are determined through Equations (14) and (15); a whole calibration process does not require zero-stress samples, and in a tensile test, it is necessary to measure a Lamb wave group velocity difference between two groups of different propagation directions; the value of K.sub.11−K.sub.22 is set as K, and then, Equation (10) is transformed into: $\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =K\left[\frac{{\sigma }_{11}}{2}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{11}}{2}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\\ +K\left[\frac{{\sigma }_{22}}{2}\left(-\cos \left(2{\theta }_1\right)+\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{22}}{2}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)\right]\\ +K\left[{\sigma }_{12}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\\ +K_{12}\left[{\sigma }_{12}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\end{array};& \left(16\right)\end{array}$ in Equation (16), dv is a measured value, K and K.sub.12 are known quantities, and three components of the plane stress field are unknown quantities, therefore, by solving the equations established by three groups of dv, two normal stresses and one shear stress is obtained; and in three times of dv measurement, it is not necessary to use the data in a zero-stress state as reference data, therefore, the measurement of the plane stress field without baseline data is realized finally through multi-angle acquisition of piezoelectric transducer array sensors.

5. The method according to claim 4, wherein pre-stresses are applied to the measured object in different directions, and then, the influence of the stresses on the S0-mode group velocity is analyzed; a stress of T.sub.33=100 MPa is applied to the measured object to obtain a velocity change in each propagation direction relative to a zero-stress state within a range of 0-360° in a plane perpendicular to a z axis; and 20 MPa is taken as a step value, T.sub.33 is changed from 0 MPa to 100 MPa to obtain velocity difference changes of 0° and 90° as well as 0° and 45°, and acoustoelastic coefficients K and K.sub.12 of the measured object are calculated according to Equations (14) and (15).

6. The method according to claim 5, wherein after the calibration of the acoustoelastic coefficients is completed, actual measurement of the stress state of the measured structure is performed, an arbitrary waveform generator is used to generate low-voltage modulation signals which are then amplified by a power amplifier to excite piezoelectric wafers to generate Lamb waves, S0-mode Lamb waves propagate inside the measured object and are received by three piezoelectric wafer sensors, received signals are mV-order weak signals, the weak signals are amplified and filtered, a lower cut-off frequency and an upper cut-off frequency of a filter are set according to bandwidths of excitation signals, the amplified and filtered signals are input into a high-speed data acquisition system, the signals are encoded by a field programmable gate array (FPGA) chip, and sampled signals are transmitted into a host computer by a PXIE bus and are stored for subsequent data processing.

7. The method according to claim 6, wherein Hilbert transform is performed on the signals received by the piezoelectric wafer sensors; after the Hilbert transform, the amplitude envelopes of the received signals is extracted, and then, a propagation acoustic time difference between S0-mode Lamb waves in different propagation directions is calculated by a cross-correlation algorithm; since the propagation distances between transmitting transducers and receiving transducers are the same, for isotropic materials, the relationship between a propagation time difference and a velocity difference meets: $\begin{array}{cc}\Delta t=t_1-t_2=\frac{L}{c_1}-\frac{L}{c_2}=\frac{L\left(c_2-c_1\right)}{c_1{c}_2}\approx -\frac{L\Delta c}{c^2},& \left(17\right)\end{array}$ in the equation, L represents a distance between receiving and transmitting transducers, and c represents a propagation velocity when there is no stress; and in Equation (17) there is an approximate linear relationship between acoustic time differences in different propagation directions and stresses, therefore, the acoustic time differences obtained by the cross-correlation algorithm are substituted into the calibrated acoustoelastic equations to determine the magnitude of the plane stress field inside the measured object at this time.

8. A system for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array, comprising: a dispersion curve establishing module, wherein by using the dispersion curve establishing module, according to geometric dimensions and material parameters of a measured object, a dispersion curve of Lamb waves is established, a cut-off frequency of a first-order Lamb wave mode is determined, then a suitable frequency is selected within the cut-off frequency to excite piezoelectric wafers, and pure S0-mode Lamb waves are generated by means of symmetrical excitation; a stress field characterization module, wherein by using the stress field characterization module, according to the invariant property of acoustoelastic coefficients of isotropic materials in different coordinate systems, acoustoelastic equations in different propagation directions are obtained, then a difference is calculated to obtain a linear relationship between the difference of velocity changes in different directions and the stress, and at the same time, according to a linear superposition relationship between stresses in different directions, the magnitude of a plane stress field of the measured object is calculated through the difference of three groups of velocity changes in different directions, thereby finally realizing the accurate characterization of the plane stress field by means of measurement results of the piezoelectric transducer array in different directions; and a stress field solving module, wherein by using the stress field solving module, Lamb wave signals are excited and received inside the measured object, Hilbert transform is performed on the received signals, signal envelopes are extracted, then acoustic time differences in different propagation directions are calculated through a cross-correlation algorithm, according to a linear relationship between acoustic time differences and velocity differences, it is determined that there is still a linear relationship between stresses and acoustic time differences in different propagation directions, and then, three groups of acoustic time differences in different directions are substituted into the pre-calibrated acoustoelastic equations to solve the plane stress field of the measured object.

9. An electronic device, comprising a memory and a processor, wherein the memory stores computer programs, and when the processor executes the computer programs, the steps of the methods according to claim 1 are implemented.

10. A computer-readable storage medium for storing computer instructions, wherein when the computer instructions are executed by a processor, the steps of the methods according to claim 1 are implemented.

Description

BRIEF DESCRIPTION OF FIGURES

[0043] FIG. 1 is a schematic diagram of a symmetrical excitation manner of S0-mode Lamb waves, where 1 and 2 respectively represent a piezoelectric wafer exciter, and 3 represents an object to be measured.

[0044] FIG. 2 is a schematic diagram of a group velocity dispersion curve of a 1 mm aluminum plate 6061.

[0045] FIG. 3 is a schematic diagram of a measured object subjected to a plane stress field.

[0046] FIG. 4 is a diagram of group velocity changes in different directions when F=500 kHz and T.sub.33=100 MPa.

[0047] FIG. 5 is a schematic diagram of a difference between velocity changes in 0° and 90° directions when F=500 kHz and T.sub.33=100 MPa.

[0048] FIG. 6 is a schematic diagram of a difference between velocity changes in 0° and 45° directions when F=500 kHz and T.sub.33=100 MPa.

[0049] FIG. 7 is an integrated schematic diagram of a measurement system, where 1 represents a tensiometer, 2 represents an arbitrary waveform generator, 3 represents a power amplifier, 4 represents a high-bandwidth receiving and amplifying device, 5 represents a three-channel high-speed data acquisition system, 6 represents a personal computer (PC), 7 represents a piezoelectric wafer exciter, and 8, 9 and 10 respectively represent a piezoelectric wafer sensor.

DETAILED DESCRIPTION

[0050] The technical solutions in the examples of the disclosure are clearly and completely described below with reference to the drawings in the examples of the disclosure. It is apparent that the described examples are only a part of the examples of the disclosure, and are not all of the examples. Based on the examples of the disclosure, all other examples obtained by a person of ordinary skill in the art without involving any inventive effort all fall within the protection scope of the disclosure.

[0051] Aiming at the defects of the existing methods, in order to meet the urgent need for online monitoring of a plane stress field without baseline data, the disclosure invents a method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array. On the basis of the acoustoelastic effect of Lamb waves, the linear relationship between the stress and the difference of velocity changes of the same-mode Lamb waves in different propagation directions is innovatively used to realize the accurate characterization of the plane stress field without the need of zero-stress standard sample reference data. Since Lamb waves have complex multi-mode characteristics, a suitable excitation frequency needs to be selected according to geometric dimensions of the structure to be measured, and then, only low-order mode Lamb waves are excited inside the measured structure to avoid serious waveform aliasing. For isotropic measured objects, anisotropic characteristics will be generated under the action of pre-stresses, that is, the propagation velocities of ultrasonic waves in all directions are different, but there is a linear relationship between velocity changes in different propagation directions and stresses. Therefore, there is still a linear relationship between the difference of velocity changes in different propagation directions and the stress. According to this characteristic, a characterization method of an absolute stress field without baseline data can be implemented. Furthermore, the method of the disclosure can make full use of the low attenuation characteristics of the Lamb waves to realize online monitoring of the plane stress field with a large coverage area.

[0052] The technical solution adopted by the disclosure is a method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array. First, according to geometric dimensions and material parameters of a measured object, a dispersion curve of Lamb waves is established to ensure that an excitation frequency of Lamb signals is lower than a first-order cut-off frequency, so as to avoid complex signal aliasing caused by the appearance of high-order modes. Previous studies have found that S0-mode Lamb waves have higher sensitivity to stresses, therefore, in the disclosure, the generation of an A0 mode is inhibited by means of symmetrical excitation of upper and lower surfaces, so as to ensure that only the S0 mode is excited inside the measured object. Then, through theoretical and simulated analysis, the linear relationship between S0-mode Lamb waves in different propagation directions and stresses is obtained, so as to further determine the linear relationship between the difference of velocity changes in different propagation directions and the stress, thereby finally realizing the accurate characterization of the plane stress field by means of measurement results of the piezoelectric transducer array in different directions. According to the selected excitation frequency, excitation signals modulated by a Hanning window are loaded inside a signal generator and are amplified by a power amplifier to drive a transmitting probe, the signals of a receiving probe at the other end are amplified by first-order weak signals, then received by a three-channel high-speed acquisition board and transmitted to a host computer through a PXIE bus and stored, and then, data processing is performed. Hilbert transform is performed on the signals, signal envelopes are extracted, then the acoustic time differences between Lamb wave signals in different propagation directions are calculated by a cross-correlation algorithm, and the results are substituted into the pre-calibrated acoustoelastic equations set to determine a mean stress field between the receiving and transmitting probes, thereby realizing the state characterization of the plane stress field of the measured object.

[0053] With reference to FIG. 1 to FIG. 7, the disclosure provides a method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array, including:

[0054] according to geometric dimensions and material parameters of a measured object, a dispersion curve of Lamb waves is established, a cut-off frequency of a first-order Lamb wave mode is determined, then a suitable frequency is selected within the cut-off frequency to excite piezoelectric wafers, and pure S0-mode Lamb waves are generated by means of symmetrical excitation;

[0055] according to the invariant property of acoustoelastic coefficients of isotropic materials in different coordinate systems, acoustoelastic equations in different propagation directions are obtained, then a difference is calculated to obtain a linear relationship between the difference of velocity changes in different directions and the stress, and at the same time, according to a linear superposition relationship between stresses in different directions, the magnitude of a plane stress field of the measured object is calculated through the difference of three groups of velocity changes in different directions, thereby finally realizing the accurate characterization of the plane stress field by means of measurement results of the piezoelectric transducer array in different directions; and

[0056] Lamb wave signals are excited and received inside the measured object, Hilbert transform is performed on the received signals, signal envelopes are extracted, then acoustic time differences in different propagation directions are calculated through a cross-correlation algorithm, according to a linear relationship between acoustic time differences and velocity differences, it is determined that there is still a linear relationship between stresses and acoustic time differences in different propagation directions, and then, three groups of acoustic time differences in different directions are substituted into the pre-calibrated acoustoelastic equations to solve the plane stress field of the measured object.

[0057] According to the geometric dimensions and the material parameters of the measured object, the dispersion curve of Lamb waves is established, and a Lamb wave dispersion equation of a free plate is a Rayleigh-Lamb wave dispersion equation, which meets:

[00009]$\begin{array}{cc}\frac{\tan \left(qh\right)}{\tan \left(ph\right)}=-\frac{4k^{2}pq}{{\left(q^2-k^2\right)}^2},& \left(1\right)\end{array}$$\begin{array}{cc}\frac{\tan \left(qh\right)}{\tan \left(ph\right)}=-\frac{{\left(q^2-k^2\right)}^2}{4k^{2}pq},& \left(2\right)\end{array}$

[0058] where p and q are respectively expressed as:

[00010]$\begin{array}{cc}p=\sqrt{\frac{{\omega }^2}{c_L^2}-k^2},& \left(3\right)\end{array}$$\begin{array}{cc}q=\sqrt{\frac{{\omega }^2}{c_T^2}-k^2},& \left(4\right)\end{array}$

[0059] in the equation, c.sub.L and c.sub.T respectively represent velocities of longitudinal waves and transverse waves, h represents a half of a thickness of the plate, ω represents an angular frequency of ultrasonic waves, and k represents a wave number.

[0060] Assuming that the object to be measured is an aluminum plate 6061 with a thickness of 1 mm, a dispersion curve of the structure to be measured can be obtained by solving Equations (1) and (2). According to the dispersion curve, it can be determined that the cut-off frequency of the first-order Lamb wave mode is 1.6 MHz, therefore, modulation signals with a frequency of 500 kHz are selected to excite the Lamb waves.

[0061] Assuming that the measured object is subjected to plane stresses, in a coordinate system x-o-y, normal stresses are σ.sub.11 and σ.sub.22, and shear stresses are σ.sub.12 and σ.sub.21; according to the symmetry of the stresses, the values of σ.sub.12 and σ.sub.21 are equal, therefore, for Lamb waves propagating in an x direction, the acoustoelastic effect of the plane stress field can be expressed as:

Δv=K.sub.11σ.sub.11+K.sub.12σ.sub.12+K.sub.22σ.sub.22  (5),

[0062] in the equation, K.sub.11, K.sub.12 and K.sub.22 respectively represent acoustoelastic coefficients of the Lamb waves in different propagation directions;

[0063] when the measured object belongs to an isotropic material, the acoustoelastic coefficients in different coordinate systems are exactly the same, therefore, the acoustoelastic effect of the Lamb waves propagating in an x′ direction can be expressed as:

Δv=K.sub.11σ′.sub.11+K.sub.12σ′.sub.12+K.sub.22σ′.sub.22  (6),

[0064] in the equation, σ′.sub.11, σ′.sub.12 and σ′.sub.22 represent normal stresses and shear stresses in a new coordinate system x′-o-y′;

[0065] according to a coordinate transformation relationship between stresses, the following can be obtained:

[00011]$\begin{array}{cc}{\sigma }_{11}^{\prime }=\frac{{\sigma }_{11}+{\sigma }_{22}}{2}+\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\cos \left(2\theta \right)+{\sigma }_{12}\sin \left(2\theta \right),& \left(7\right)\end{array}$$\begin{array}{cc}{\sigma }_{22}^{\prime }=\frac{{\sigma }_{11}+{\sigma }_{22}}{2}-\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\cos \left(2\theta \right)-{\sigma }_{12}\sin \left(2\theta \right),& \left(8\right)\end{array}$$\begin{array}{cc}{\sigma }_{12}^{\prime }=-\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\sin \left(2\theta \right)+{\sigma }_{12}\cos \left(2\theta \right),& \left(9\right)\end{array}$

[0066] where θ represents an angle;

[0067] a difference between propagation velocities along angles θ.sub.1 and θ.sub.2 is:

[00012]$\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =K_{11}[\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\\ +{\sigma }_{12}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)]\\ +K_{22}[\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\left(-\cos \left(2{\theta }_1\right)+\cos \left(2{\theta }_2\right)\right)\\ +{\sigma }_{12}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)]\\ +K_{12}[\frac{{\sigma }_{11}-{\sigma }_{22}}{2}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)\\ +{\sigma }_{12}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)]\end{array};& \left(10\right)\end{array}$

[0068] it can be seen from Equation (10) that the influence of each component of the plane stress field on dv has a linear superposition relationship, and when the plane stress field has only one component, dv can be respectively simplified as:

[00013]$\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =\left(K_{11}-K_{22}\right)\left[\frac{{\sigma }_{11}}{2}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{11}}{2}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\end{array},& \left(11\right)\end{array}$$\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =\left(K_{11}-K_{22}\right)\left[\frac{{\sigma }_{22}}{2}\left(-\cos \left(2{\theta }_1\right)+\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{22}}{2}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)\right]\end{array},& \left(12\right)\end{array}$$\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =\left(K_{11}-K_{22}\right)\left[{\sigma }_{12}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\\ +K_{12}\left[{\sigma }_{12}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\end{array};& \left(13\right)\end{array}$

[0069] it can be found that to complete the measurement of the plane stress field, it is only necessary to determine the values of K.sub.11−K.sub.22 and K.sub.12 through calibration tests;

[0070] assuming that there is only σ.sub.11, θ.sub.1 is 0° and θ.sub.2 is 90°, the following can be obtained:

dv=Δv(θ.sub.1)−Δv(θ.sub.2)=(K.sub.11−K.sub.22)σ.sub.11  (14);

[0071] assuming that θ.sub.1 is 0° and θ.sub.2 is 45°, when there is σ.sub.11, dv is expressed as:

[00014]$\begin{array}{cc}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)=\left(K_{11}-K_{22}\right)\frac{{\sigma }_{11}}{2}+K_{12}\frac{{\sigma }_{11}}{2};& \left(15\right)\end{array}$

[0072] two acoustoelastic values of K.sub.11−K.sub.22 and K.sub.12 are determined through Equations (14) and (15); a whole calibration process does not require zero-stress samples, and in a tensile test, it is necessary to measure a Lamb wave group velocity difference between two groups of different propagation directions;

[0073] the value of K.sub.11−K.sub.22 is set as K, and then, Equation (10) can be transformed into:

[00015]$\begin{array}{cc}\begin{array}{c}\mathrm{dv}=\Delta v\left({\theta }_1\right)-\Delta v\left({\theta }_2\right)\\ =K\left[\frac{{\sigma }_{11}}{2}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{11}}{2}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\\ +K\left[\frac{{\sigma }_{22}}{2}\left(-\cos \left(2{\theta }_1\right)+\cos \left(2{\theta }_2\right)\right)\right]\\ -K_{12}\left[\frac{{\sigma }_{22}}{2}\left(-\sin \left(2{\theta }_1\right)+\sin \left(2{\theta }_2\right)\right)\right]\\ +K\left[{\sigma }_{12}\left(\sin \left(2{\theta }_1\right)-\sin \left(2{\theta }_2\right)\right)\right]\\ +K_{12}\left[{\sigma }_{12}\left(\cos \left(2{\theta }_1\right)-\cos \left(2{\theta }_2\right)\right)\right]\end{array};& \left(16\right)\end{array}$

[0074] in Equation (16), dv is a measured value, K and K.sub.12 are known quantities, and three components of the plane stress field are unknown quantities, therefore, by solving the equations established by three groups of dv, two normal stresses and one shear stress can be obtained; and in three times of dv measurement, it is not necessary to use the data in a zero-stress state as reference data, therefore, the measurement of the plane stress field without baseline data can be realized finally through multi-angle acquisition of piezoelectric transducer array sensors.

[0075] Pre-stresses are applied to the measured object in different directions, and then, the influence of the stresses on the S0-mode group velocity is analyzed; a stress of T.sub.33=100 MPa is applied to the measured object to obtain a velocity change in each propagation direction relative to a zero-stress state within a range of 0-360° in a plane perpendicular to a z axis; and under the action of stresses, an isotropic material exhibits anisotropy, that is, propagation velocities in different directions are different. When the measured object is subjected to a stress in a T.sub.33 direction, the velocity change sensitivity of an S0 mode is better than that of an A0 mode. But under the action of T.sub.33, there is an angle which enables the velocity change of the S0 mode to be zero. Therefore, when measurement is performed with data in a single direction, if the relationship between a stress direction and a measurement direction is unknown, this may cause incorrect results. This problem can be effectively solved by the multi-directional data collected by the piezoelectric array. 20 MPa is taken as a step value, T.sub.33 is changed from 0 MPa to 100 MPa to obtain velocity difference changes of 0° and 90° as well as 0° and 45°, and acoustoelastic coefficients K and K.sub.12 of the measured object are calculated according to Equations (14) and (15).

[0076] After the calibration of the acoustoelastic coefficients is completed, actual measurement of the stress state of the measured structure can be performed, an arbitrary waveform generator is used to generate low-voltage modulation signals which are then amplified by a power amplifier to excite piezoelectric wafers to generate Lamb waves, S0-mode Lamb waves propagate inside the measured object and are received by three piezoelectric wafer sensors, received signals are mV-order weak signals, the weak signals are amplified and filtered, a lower cut-off frequency and an upper cut-off frequency of a filter are set according to bandwidths of excitation signals, the amplified and filtered signals are input into a high-speed data acquisition system, the signals are encoded by an FPGA chip, and sampled signals are transmitted into a host computer by a PXIE bus and are stored for subsequent data processing.

[0077] Hilbert transform is performed on the signals received by the piezoelectric wafer sensors; after the Hilbert transform, the amplitude envelopes of the received signals can be extracted, and then, a propagation acoustic time difference between S0-mode Lamb waves in different propagation directions can be calculated by a cross-correlation algorithm; since the propagation distances between transmitting transducers and receiving transducers are the same, for isotropic materials, the relationship between a propagation time difference and a velocity difference meets:

[00016]$\begin{array}{cc}\Delta t=t_1-t_2=\frac{L}{c_1}-\frac{L}{c_2}=\frac{L\left(c_2-c_1\right)}{c_1{c}_2}\approx -\frac{L\Delta c}{c^2},& \left(17\right)\end{array}$

[0078] in the equation, L represents a distance between receiving and transmitting transducers, and c represents a propagation velocity when there is no stress; and

[0079] it can be seen from Equation (17) that there is an approximate linear relationship between acoustic time differences in different propagation directions and stresses, therefore, the acoustic time differences obtained by the cross-correlation algorithm are substituted into the calibrated acoustoelastic equations to determine the magnitude of the plane stress field inside the measured object at this time.

[0080] The disclosure further provides a system for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array, including:

[0081] a dispersion curve establishing module, where by using the dispersion curve establishing module, according to geometric dimensions and material parameters of a measured object, a dispersion curve of Lamb waves is established, a cut-off frequency of a first-order Lamb wave mode is determined, then a suitable frequency is selected within the cut-off frequency to excite piezoelectric wafers, and pure S0-mode Lamb waves are generated by means of symmetrical excitation;

[0082] a stress field characterization module, where by using the stress field characterization module, according to the invariant property of acoustoelastic coefficients of isotropic materials in different coordinate systems, acoustoelastic equations in different propagation directions are obtained, then a difference is calculated to obtain a linear relationship between the difference of velocity changes in different directions and the stress, and at the same time, according to a linear superposition relationship between stresses in different directions, the magnitude of a plane stress field of the measured object is calculated through the difference of three groups of velocity changes in different directions, thereby finally realizing the accurate characterization of the plane stress field by means of measurement results of the piezoelectric transducer array in different directions; and

[0083] a stress field solving module, where by using the stress field solving module, Lamb wave signals are excited and received inside the measured object, Hilbert transform is performed on the received signals, signal envelopes are extracted, then acoustic time differences in different propagation directions are calculated through a cross-correlation algorithm, according to a linear relationship between acoustic time differences and velocity differences, it is determined that there is still a linear relationship between stresses and acoustic time differences in different propagation directions, and then, three groups of acoustic time differences in different directions are substituted into the pre-calibrated acoustoelastic equations to solve the plane stress field of the measured object.

[0084] The disclosure further provides an electronic device, including a memory and a processor. The memory stores computer programs, and when the processor executes the computer programs, the steps of the method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array are implemented.

[0085] The disclosure further provides a computer-readable storage medium for storing computer instructions. When the computer instructions are executed by the processor, the steps of the method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array are implemented.

[0086] Under normal temperature conditions, the measured object is an aluminum plate 6061 with a thickness of 1 mm, which is an isotropic material, an arbitrary waveform generator is used to generate a 500 kHz pulse waveform modulated by a Hanning window, low-frequency signals generated by a signal generator are subjected to first-stage amplification by an Aigtek power amplifier, piezoelectric wafer exciters are excited by high-voltage signals, trigger signals are generated at the same time, and ultrasonic signals obtained by piezoelectric wafer sensors are acquired by a high-speed data acquisition board. Before the received signals are acquired, the signals are subjected to signal weakening and band-pass filtering first, the signals are amplified to an input voltage range of the data acquisition board, low-frequency noise and high-frequency noise in the signals are filtered out at the same time, the received signals are acquired continuously for 16 times, and the acquired signals are subjected to smooth filtering to filter out a part of electronic noise to improve a signal-to-noise ratio of the received signals. Then, the acoustic time differences of S0-mode Lamb waves in different propagation directions are calculated by means of the Hilbert transform and the cross-correlation algorithm, the acoustic time differences are substituted into the pre-calibrated acoustoelastic equations, and the whole plane stress field of the measured object can be calculated through three groups of acoustic time differences in different propagation directions.

[0087] The method for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array, provided by the disclosure, innovatively considers the combined action of multi-directional Lamb waves, and according to the linear relationship between velocity changes in different propagation directions under the action of stresses and the stresses, linear equations between velocity differences in different propagation directions and stresses are obtained, so that online monitoring of a state of an absolute plane stress field of the measured structure can be realized without any reference baseline data.

[0088] The method, system, device and medium for online monitoring of a plane stress field without baseline data based on a piezoelectric transducer array, provided by the disclosure, have been described in detail above. The principles and implementations of the disclosure are described by specific examples in the text. The description of the above examples is only used to help understand the method and core idea of the disclosure. Meanwhile, for those skilled in the art, according to the idea of the disclosure, there will be changes in specific implementations and application scope. In conclusion, the content of the description shall not be construed as limiting the disclosure.