Disclosure of Invention
The invention aims to provide a wing defect detection device and method for solving the problems in the prior art.
In order to achieve the above purpose, the invention provides a wing defect detection device, comprising:
The lifting platform and the guide rail platform are respectively fixed on the ground and are arranged in parallel;
The rotating assembly comprises a connecting disc and a main shaft, the connecting disc is rotationally connected to the top end of the lifting platform through the main shaft, two turntable joints are fixedly connected to one side, away from the lifting platform, of the connecting disc, and the wing is connected with the turntable joints;
the detection assembly comprises a triaxial mobile station and a laser ultrasonic probe, the triaxial mobile station is arranged on the guide rail platform, and the laser ultrasonic probe is fixedly connected to the triaxial mobile station and used for detecting the wing;
The laser ultrasonic probe is a single-excitation-point double-detection-point laser ultrasonic probe.
Preferably, the two turntable joints are respectively positioned at two sides of the connecting disc and positioned on the same diameter line.
Preferably, the lifting table is fixedly connected with a transmission case, a servo motor is fixedly connected in the transmission case, an output shaft of the servo motor is connected with the main shaft through a speed reducer, a coupler and an absolute photoelectric encoder, and a gap is arranged between the main shaft and the lifting table.
A wing defect detection method comprising the steps of:
s1, mounting the wing on the connecting disc for laser ultrasonic scanning;
s2, performing subsurface crack defect detection on the wing;
S3, detecting surface crack defects of the wing;
s4, performing three-dimensional imaging on the wing to obtain a three-dimensional image of the internal defect of the wing.
Preferably, in the step S1, the wing is kept parallel to the horizontal plane in the initial position, the height of the lifting platform is H J, the distance between the main shaft and the top surface of the lifting platform is H z, and the maximum wing width of the wing is W max, so that H J+hz≥1.5Wmax is satisfied;
The height of the guide rail platform is H, the initial height of the Z axis of the triaxial mobile station is Z pt0, and H J+hz=H+Zpt0 is met.
Preferably, during laser ultrasonic scanning, a global rectangular coordinate system X-Y-Z is firstly established, the X direction points to the wingspan direction of the wing at the initial position, the Y direction points to the wingspan direction vertical to the X direction, the XOY plane is parallel to the horizontal plane, the direction vertical to the XOY plane is the Z direction, then the surface of the wing is scattered into points, laser ultrasonic scanning is carried out along a discretization curve by connecting a plurality of discrete points into a line, a local coordinate system Lm-Ln-Lz is established according to the scanning curve, and a wingspan direction data set p (I 1i,Lm,t),p(I1i, lm, t) and a wingspan direction data set p (J j, ln, t) are obtained and used for imaging defects.
Preferably, the subsurface crack defect is a planar defect, and in the case that a subsurface defect is located in the wingspan section LnOLz corresponding to the wingspan curve J s, the detection of the subsurface defect is realized by utilizing the wingspan data set p (J j, ln, t), and the specific steps of the subsurface crack defect detection process are as follows:
Determining the vertex coordinates Ln 0 of subsurface defects, namely, in the scanning process, only scanning data p (J s, ln, t) on a wing circumferential direction curve J s, generating a defect echo signal RL, and observing the defect echo signal RL of a scanning dataset to determine the wing circumferential section where the subsurface defects are located;
Observing scanning data p (J s, ln, t), and obtaining subsurface defect vertex coordinates Ln 0 when the arrival time t RL of the reflected wave RL obtains the minimum value and the defect vertex is positioned at the midpoint of an excitation point and a detection point II in the Ln direction;
Determination of subsurface defect vertex coordinates Lz 0, wherein the defect vertex coordinates are (Ln 0,Lz0), the coordinates of excitation points are (Ln 1, 0), the coordinates of detection points two are (Ln 1 +l, 0), l is the distance between the excitation points and the detection points two, and c is the wave velocity, then the subsurface crack burial depth Lz 0 can be calculated by the following formula:
Sub-surface defect detection of spanwise cross-sections is performed by scanning data sets p (I 1i, lm, t) and p (I 2i, lm, t).
Preferably, the surface crack defect is a planar defect, and in the case that a surface defect is located in the span section LmOLz corresponding to the span curve I 1s, the specific steps of surface defect detection are as follows:
Sequentially drawing all scanned spanwise cross-section B pictures according to a scanned data set p (I 1i, lm, t);
In the scanning process, the excitation point and the detection point are fixed, when no defect exists between the excitation point and the detection point, the propagation time of the surface acoustic wave from the excitation point to the detection point is the same, and B scanning imaging is carried out to form a continuous straight line;
When a defect exists between the excitation point and the detection point, most waves are reflected when the surface waves are transmitted to the side edge of the defect from the excitation point, the detected signal intensity at the detection point is very low, and the scanning imaging is carried out to form a straight line with a gap, so that the wingspan section I 1s where the position of the surface defect is located is determined;
The notch position in the B scanning chart is the position of the surface defect, and spanwise coordinates Lms1 and Lms2 can be obtained according to the B scanning chart;
surface defect detection of the airfoil circumferential section was performed by p (J s, ln, t).
Preferably, three-dimensional imaging can be achieved in two ways according to the scanned data sets p (I 1i,Lm,t),p(I2i, lm, t) and p (J j, ln, t), and the position and shape of the internal defect are presented;
Extracting a two-dimensional dataset p 1i (Lm, t) corresponding to a spanwise curve I 1i for focused imaging of a two-dimensional spanwise section corresponding to a spanwise curve I 1i, wherein the specific steps include calculating an equivalent velocity in a frequency-wavenumber domain by a numerical calculation method, calculating a phase shift factor according to the equivalent velocity, and calculating a two-dimensional spectrum extrapolation and an imaging matrix, wherein the equivalent velocity can be calculated by the following formula:
Changing imaging depth Lz, repeating the steps until the amplitude value calculation of all imaging points is completed, and obtaining an imaging matrix I (Lm, lz) of a two-dimensional span section corresponding to a span curve I 1i;
Collecting all spanwise cross-section imaging matrixes to obtain a three-dimensional imaging matrix I (I 1i, lm, lz) of the wing in the section L 1;
And carrying out coordinate conversion to obtain a three-dimensional imaging matrix I (x, y, z) under a rectangular coordinate system, and imaging.
Preferably, the specific steps of three-dimensional imaging of the wing based on the wing circumferential scan data set p (J j, ln, t) are as follows:
Dividing the peripheral section of the wing into a plurality of contour sections, expanding the outer contour into a straight line, expanding the peripheral section of the wing in the imaging section into a plane, and performing frequency domain imaging in the imaging section to obtain a plane frequency domain imaging matrix in the contour section;
splicing the curved surface imaging in all the contour sections into a frequency domain imaging matrix I (Ln, lz) of a two-dimensional section of the wing circumference corresponding to a wing circumference curve J j;
And sequentially performing two-dimensional spanwise section imaging corresponding to the peri-wing curve, and collecting all peri-wing section imaging matrixes to obtain a three-dimensional imaging matrix I (J j, ln, lz) of the wing.
Compared with the prior art, the invention has the following advantages and technical effects:
1. According to the wing defect detection device provided by the invention, nondestructive detection is performed on the wing skin defect of the carbon fiber composite material, the wing is arranged on the connecting disc, the main shaft drives the wing to rotate, the triaxial mobile station adjusts the position of the laser ultrasonic probe, so that the omnibearing non-contact detection of the wing is realized, and the tangent line at the excitation point at any moment is kept vertical, thereby facilitating the detection of the laser ultrasonic probe, and further realizing the omnibearing detection of the wing.
2. According to the wing defect detection method provided by the invention, the wing is scattered into a space curve in a discrete mode aiming at the defect of the irregular surface of the wing, which is difficult to detect in a traditional scanning mode, so that the resolution is improved, and the wing defect detection method can go deep into the inside of a multilayer structure of a composite material to carry out comprehensive and accurate detection.
3. According to the invention, the laser ultrasonic probe with single excitation point and double detection points is adopted, ultrasonic signals in the wingspan direction and the wingspan direction can be measured simultaneously, and the wingspan curve scanning data set can be obtained in a simple and convenient manner by matching with the laser ultrasonic scanning mode provided by the invention, so that three-dimensional imaging of defects is performed.
4. The detection method provided by the invention can simultaneously carry out surface defects and subsurface defects in the wingspan direction and the wingspan direction, and avoid missed detection caused by crack orientation.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
Embodiment one:
referring to fig. 1 to 11, the present invention provides an apparatus for detecting wing defects, comprising:
the lifting platform 1 and the guide rail platform 2 are respectively fixed on the ground and are arranged in parallel;
the rotating assembly comprises a connecting disc 3 and a main shaft 4, wherein the connecting disc 3 is rotationally connected to the top end of the lifting platform 1 through the main shaft 4, two turntable joints 5 are fixedly connected to one side, away from the lifting platform 1, of the connecting disc 3, and wings 6 are connected with the turntable joints 5;
the detection assembly comprises a triaxial moving table 7 and a laser ultrasonic probe 8, wherein the triaxial moving table 7 is arranged on the guide rail platform 2, and the laser ultrasonic probe 8 is fixedly connected to the triaxial moving table 7 for detecting the wing 6;
the laser ultrasonic probe 8 is a single-excitation-point double-detection-point laser ultrasonic probe.
Further optimizing scheme, two carousel joints 5 are located respectively and connect disc 3 both sides, and are located same diameter line.
According to a further optimization scheme, a transmission case 9 is fixedly connected to the lifting table 1, a servo motor 10 is fixedly connected to the transmission case 9, an output shaft of the servo motor 10 is connected with a main shaft 4 through a speed reducer 11, a coupling 12 and an absolute photoelectric encoder 13, and a gap is formed between the main shaft 4 and the lifting table 1.
A wing defect detection method comprising the steps of:
S1, mounting a wing 6 on a connecting disc 3 for laser ultrasonic scanning;
S2, performing subsurface crack defect detection on the wing 6;
s3, detecting surface crack defects of the wing 6;
s4, performing three-dimensional imaging on the wing 6 to obtain a three-dimensional image of the internal defect of the wing 6.
In one embodiment of the invention, the defect detection is the detection of surface, subsurface, internal defects of the wing 6.
In the further optimization scheme, in the S1, when the wing 6 is in the initial position, the wing 6 is kept parallel to the horizontal plane, the height of the lifting platform 1 is H J, the distance between the main shaft 4 and the top surface of the lifting platform 1 is H z, and the maximum wing width of the wing 6 is W max, so that H J+hz≥1.5Wmax is satisfied;
The height of the guide rail platform 2 is H, the initial height of the Z axis of the triaxial mobile platform 7 is Z pt0, and H J+hz=H+Zpt0 is satisfied.
In a further optimization scheme, a global rectangular coordinate system X-Y-Z is firstly established during laser ultrasonic scanning, the X direction points to the wingspan direction of the wing 6 at the initial position, the Y direction points to the wingspan direction vertical to the X direction, the XOY plane is parallel to the horizontal plane, the direction vertical to the XOY plane is the Z direction, then the surface of the wing 6 is scattered into points, the laser ultrasonic scanning is carried out along a discretization curve by connecting a plurality of discrete points into a line, a local coordinate system Lm-Ln-Lz is established according to the scanning curve, and a wingspan direction data set p (I 1i,Lm,t),p(I1i, lm, t) and a wingspan direction data set p (J j, ln, t) are obtained and are used for imaging defects.
In one embodiment of the present invention, when a defect is imaged, the surface of the wing 6 is first discretized into points, as shown in fig. 3, by connecting a plurality of discrete points into a line, laser ultrasonic scanning in S1 is facilitated, the spanwise direction of the surface of the wing 6 is M, the circumferential direction of the wing is N, the distance λ between every two discrete points on the left side of the wing 6 in the N direction, the contour line in the direction is discretized into N points, the perimeter c 0 =n×λ of the contour line on the leftmost side is measured once every distance λ along the M direction, the perimeter c n of the contour line on the N direction is measured once every distance λ along the M direction, the perimeter c m of the contour line of the wing 6 at the M λ position on the leftmost side of the wing 6 is recorded, each contour line is discretized into N points, the distance between every two points of the contour line of the wing 6 in the N direction is c mλ/c0, the wing 6 is divided into M segments of contour lines in the M direction, the total length of the wing 6 in the M direction is l=m×x, the wing 6 is divided into two segments of the left and right and the length of the wing 6 is divided into two segments of the m× 1+m2 M (the two discrete points in the m×n directions of the surface of the wing 1=m1×λ,L2=m2);
Establishing a coordinate system, namely establishing a global space rectangular coordinate system X-Y-Z at the end part of the wing 6 by taking the joint center point of the two wings 6 as an origin O, wherein the X direction points to the wingspan direction of the wing 6 at an initial position, the Y direction points to the wingspan direction vertical to the X direction, the XOY plane is parallel to a horizontal plane, the direction vertical to the XOY plane is the Z direction, the coordinates (d, e and f) of each discrete point are determined through a secondary coordinate system, when the wing 6 is at the initial position, the transverse maximum plane of the wing 6 is parallel to the XOY plane, the position of the bottom plate of the triaxial moving platform 7 on the X-Y-Z coordinate is (0, e 0,fpt0)-(L,e0,fpt0), L is the wingspan length of the wing 6, e 0 is the distance between the triaxial moving platform 7 and the wing 6, and e 0≥0.5Wmax+0.5Wpt,Wpt is the width of the guide rail platform 2;
Connecting discrete points into a line, namely connecting the discrete points into curves along the N direction in two profile sections of L 1 and L 2 to form a wing 6 wing circumferential section profile curve, connecting the discrete points along the N direction with the discrete points on the next section profile curve closest to each point along the M direction to form a spanwise curve, jointly forming a reticular space curve, defining the curve along the M direction closest to the laser ultrasonic probe 8 in the initial state as I 10, defining the curve along the N direction as I 1i in the L 1 section, defining the curve along the N direction as I 2i in the L 2 section, defining the curve along the leftmost N direction of the wing 6 in the initial state as J 0, and defining the curve along the M direction as J x lambda as J j;
The method comprises the steps of L 1 section scanning, firstly solving a tangent f-f 0=k0(e-e0 of J 0 at the intersection point (d 0,e0,f0) of J 0 and I 10, calculating the inclination angle of the tangent at the point by k 0, driving the wing 6 to rotate by the connecting disc 3 to keep the tangent vertical, regulating the position of the laser ultrasonic probe 8 by the Z axis to enable a laser path to be vertical to a curve I 10 at the intersection point, scanning the leftmost end of the wing 6 at a speed v along the I 10 to the right, stopping moving when the X-direction scanning position is d=vt=L 1, driving the wing 6 to rotate by the connecting disc 3, enabling the tangent f-f 1=k1(e-e1 of J 1 at the intersection point (d 1,e1,f1) of J 1 and I 11 to be vertical, regulating the Z axis to enable the ultrasonic laser probe path to be vertical to the I 11 at the intersection point, scanning the scanning along the curve I 11 from the right to the left, and re-timing the scanning time, thus finishing a round of scanning, enabling the excitation point of the ith scanning to be (d i,ei,fi), and repeating the steps until all the scanning steps L 1i curve is completed. The laser ultrasound probe 8 operates in the scanning process as shown in fig. 4. When the laser ultrasonic probe 8 works, 1 excitation point and 2 detection points exist on the surface of the wing 6 at the same time, the coordinates of the excitation point are (d i,ei,fi), a spanwise signal is measured by the first detection point, a wingspan direction signal is measured by the second detection point, a local coordinate system Lm-Ln-Lz is established according to a scanning curve, the Lz direction is the normal direction of the surface of the wing 6 at the excitation point (d i,ei,fi), the normal direction points to the inside of the wing 6, the Lm direction points to the scanning direction of the current wingspan curve I 1i, and Ln points to the wingspan direction at the point. By scanning all L 1i curves, the scan dataset p (I 1i, lm, t) on the spanwise I 1i curve and the upper scan dataset p (J j, ln, t) of the L 1 segment spanwise curve J j can be obtained simultaneously;
L 2 section scanning, namely turning the wing 6 back to the initial position when the L 1 section scanning is finished, calculating a tangent line of J 0 at the intersection point of J 0 and L 20, keeping the tangent line vertical, scanning the tangent line in the spanwise direction along the curve L 20 by adjusting the light path of the laser ultrasonic probe 8 to be perpendicular to L 20 at the intersection point, repeating the L 1 section scanning until all the L 2i scanning is finished, and obtaining scanning data sets p (I 2i, lm, t) on the spanwise I 2i curve and upper scanning data sets p (J j, ln, t) on the L 2 section wingspan direction curve J j.
According to a further optimization scheme, subsurface crack defects are plane defects, the subsurface crack defects are detected by utilizing a peri-wing data set p (J j, ln, t) when the subsurface crack defects are located in a wing span section LnOLz corresponding to a peri-wing curve J s, and the subsurface crack defects are detected by the following steps:
Determining the vertex coordinates Ln 0 of subsurface defects, namely, in the scanning process, only scanning data p (J s, ln, t) on a wing circumferential direction curve J s, generating a defect echo signal RL, and observing the defect echo signal RL of a scanning dataset to determine the wing circumferential section where the subsurface defects are located;
Observing scanning data p (J s, ln, t), and obtaining subsurface defect vertex coordinates Ln 0 when the arrival time t RL of the reflected wave RL obtains the minimum value and the defect vertex is positioned at the midpoint of an excitation point and a detection point II in the Ln direction;
Determination of subsurface defect vertex coordinates Lz 0, wherein the defect vertex coordinates are (Ln 0,Lz0), the coordinates of excitation points are (Ln 1, 0), the coordinates of detection points two are (Ln 1 +l, 0), l is the distance between the excitation points and the detection points two, and c is the wave velocity, then the subsurface crack burial depth Lz 0 can be calculated by the following formula:
Sub-surface defect detection of spanwise cross-sections is performed by scanning data sets p (I 1i, lm, t) and p (I 2i, lm, t).
In one embodiment of the invention, S2 subsurface defect detection, wherein subsurface defects are planar defects, and the scanning process can realize subsurface defect detection on the airfoil circumferential section LnOLz and the span section LmOLz. Assuming that a subsurface defect is located in the spanwise section LnOLz corresponding to the airfoil circumferential curve J s, the subsurface defect detection is performed as follows:
In the scanning process of the section L 1 and the section L 2, only a defect echo signal RL appears when scanning data p (J s, ln, t) of a curve J s in the circumferential direction of the wing, and the section of the wing where the subsurface defect is located can be determined by observing a scanning data set p (J j, ln, t);
Observing scanning data p (J s, ln, t), and obtaining subsurface defect vertex coordinates Ln 0 when the arrival time t RL of the reflected wave RL obtains the minimum value and the defect vertex is positioned at the midpoint of an excitation point and a detection point II in the Ln direction;
Determination of subsurface Defect vertex coordinate Lz 0 As shown in FIG. 5, where defect vertex coordinate is (Ln 0,Lz0), excitation point coordinate is (Ln 1, 0), detection point two is (Ln 1 +l, 0), l is distance between excitation point and detection point two, and c is wave velocity, subsurface crack burial depth Lz 0 can be calculated by:
The sub-surface defect vertex position (J s,Ln0,Lz0) of the peripheral section of the J s wing is obtained through the steps, and the defect vertex position (x 0,y0,z0) under a rectangular coordinate system can be obtained through further coordinate transformation. Similarly, spanwise subsurface defect detection can be performed by scanning the data sets p (I 1i, lm, t) and p (I 2i, lm, t).
In a further optimized scheme, the surface crack defect is a plane defect, and when one surface defect is positioned in the spanwise section LmOLz corresponding to the spanwise curve I 1s, the specific steps of surface defect detection are as follows:
Sequentially drawing all scanned spanwise cross-section B pictures according to a scanned data set p (I 1i, lm, t);
In the scanning process, the excitation point and the detection point are fixed, when no defect exists between the excitation point and the detection point, the propagation time of the surface acoustic wave from the excitation point to the detection point is the same, and B scanning imaging is carried out to form a continuous straight line;
When a defect exists between the excitation point and the detection point, most waves are reflected when the surface waves are transmitted to the side edge of the defect from the excitation point, the detected signal intensity at the detection point is very low, and the scanning imaging is carried out to form a straight line with a gap, so that the wingspan section I 1s where the position of the surface defect is located is determined;
The notch position in the B scanning chart is the position of the surface defect, and spanwise coordinates Lms1 and Lms2 can be obtained according to the B scanning chart;
surface defect detection of the airfoil circumferential section was performed by p (J s, ln, t).
In one embodiment of the present invention, S3 surface defect detection is performed, wherein the surface defect is a planar defect. The location of the defect may be determined by presenting the scan dataset in the form of a B-scan. Assuming that a surface defect is located in the spanwise section LmOLz corresponding to the spanwise curve I 1s, the specific steps for surface defect inspection are as follows:
Sequentially drawing all scanned spanwise cross-section B pictures according to a scanned data set p (I 1i, lm, t);
In the scanning process, the excitation point and the detection point are fixed, when no defect exists between the excitation point and the detection point, the propagation time of the surface acoustic wave from the excitation point to the detection point is the same, and B scanning imaging is carried out to form a continuous straight line;
as shown in fig. 6, when a defect exists between the excitation point and the detection point, most of waves are reflected when the surface wave passes from the excitation point to the side edge of the defect, the detected signal intensity at the detection point is very low, and the detected signal intensity is scanned and imaged into a straight line with a gap, so that the spanwise cross section I 1s where the position of the surface defect is located is determined;
and the notch position in the B scanning chart is the position of the surface defect. As shown in fig. 7, spanwise coordinates Lms1 and Lms2 are obtained from the B-scan. In the same way, the surface defect detection of the airfoil circumferential section can be performed by p (J s, ln, t).
According to a further optimization scheme, three-dimensional imaging can be realized in two ways according to the scanned data sets p (I 1i,Lm,t),p(I2i, lm, t) and p (J j, ln, t), and the positions and the shapes of internal defects are presented;
Extracting a two-dimensional dataset p 1i (Lm, t) corresponding to a span curve I 1i for focusing imaging of a two-dimensional span section corresponding to a span curve I 1i, wherein the specific steps include calculating an equivalent speed in a frequency-wave number domain by a numerical calculation method, calculating a phase shift factor according to the equivalent speed, and calculating a two-dimensional spectrum extrapolation and an imaging matrix, wherein the equivalent speed can be calculated by the following formula:
Changing imaging depth Lz, repeating the steps until the amplitude value calculation of all imaging points is completed, and obtaining an imaging matrix I (Lm, lz) of a two-dimensional span section corresponding to a span curve I 1i;
Collecting all spanwise imaging matrixes to obtain a three-dimensional imaging matrix I (I 1i, lm, lz) of the L 1 -section wing 6;
And carrying out coordinate conversion to obtain a three-dimensional imaging matrix I (x, y, z) under a rectangular coordinate system, and imaging.
Further optimizing scheme, the specific steps of three-dimensional imaging of the wing 6 based on the wing circumferential direction scanning data set p (J j, ln, t) are as follows:
Dividing the peripheral section of the wing into a plurality of contour sections, expanding the outer contour into a straight line, expanding the peripheral section of the wing in the imaging section into a plane, and performing frequency domain imaging in the imaging section to obtain a plane frequency domain imaging matrix in the contour section;
splicing the curved surface imaging in all the contour sections into a frequency domain imaging matrix I (Ln, lz) of a two-dimensional section of the wing circumference corresponding to a wing circumference curve J j;
Sequentially performing two-dimensional spanwise cross-section imaging corresponding to the peri-wing curve, and collecting all peri-wing cross-section imaging matrixes to obtain a three-dimensional imaging matrix I (J j, ln, lz) of the wing 6
In one embodiment of the invention, three-dimensional imaging of defects inside the S4 wing can be realized by scanning the data sets p (I 1i,Lm,t),p(I2i, lm, t) and p (J j, ln, t), so that the defects can be more intuitively presented. An example of three-dimensional imaging of an L 1 -section airfoil 6 of the present invention based on a spanwise scan dataset p (I 1i, lm, t) is shown below:
In step s4.1.1, a two-dimensional dataset p 1i (Lm, t) corresponding to the spanwise curve I 1i is extracted for focused imaging of a two-dimensional spanwise section corresponding to the spanwise curve I 1i. The specific implementation steps of the two-dimensional focusing imaging algorithm are as follows:
At step S4.1.2, an equivalent speed c eq is calculated. The calculation principle of the equivalent speed c eq is shown in fig. 8. The coordinates of the target point G are (Lm 0, lz 0), the coordinates of the excitation point are (LmE, 0), the coordinates of the corresponding detection point I are (LmR, 0), the distance exists between the two points, the equivalent detection point M (LmR-l/2, 0) is positioned at the midpoint of the two points, gamma T is the excitation angle, gamma R is the detection angle, gamma M is the equivalent detection angle, and c is the actual wave velocity. The calculation process of the equivalent speed c eq is as follows:
calculating excitation angles gamma T of different target points when detecting the position change, detecting angles gamma R and equivalent detection angles gamma M, and obtaining a functional relation f (gamma T,γR,γM) of gamma T、γR and gamma M;
the wave number domain excitation angle gamma 'T and the wave number domain detection angle gamma' R are calculated in the wave number domain, and the calculation process can be represented by the following formula:
kLm=k(sinγ′T+sinγ′R)
Wherein k Lm is the Lm-direction wave number, k is the wave number, k=ω/c, ω is the frequency, determined by the spanwise sweep step;
The calculated wave number domain excitation angle gamma ' T and wave number domain detection angle gamma ' R calculate a wave number domain equivalent detection angle gamma ' M according to a functional relation f (gamma T,γR,γM);
calculating an equivalent wave speed c eq from the wave number domain excitation angle gamma ' T, the wave number domain detection angle gamma ' R, the wave number domain equivalent detection angle gamma ' M and the actual wave speed c, wherein the calculation process is shown in the following formula:
Performing two-dimensional Fourier transform on a two-dimensional dataset P 1i (Lm, t) corresponding to the spanwise curve I 1i to obtain a two-dimensional frequency spectrum P (k Lm, 0, omega) with the detection depth lz=0;
The two-dimensional spectrum of depth lz=0 is extrapolated to the plane of depth Lz, resulting in the two-dimensional spectrum P (k Lm, 0, ω) of depth Lz, calculated as follows:
P(kLm,Lz,ω)=P(kLm,0,ω)α(kLm,Lz,ω)
Where α (k Lm, lz, ω) is a phase shift factor, which can be calculated by the following formula, where jj represents an imaginary unit:
Performing two-dimensional Fourier transform on Lm and Lz by using a two-dimensional frequency spectrum P (k Lm, lz, omega) of a depth Lz plane to obtain an imaging matrix of the depth Lz, wherein t=0 in the inverse transform process, and the amplitude I (Lm, lz) of an imaging point (Lm, lz) can be calculated by the following formula:
Changing imaging depth Lz, repeating the steps until the amplitude value calculation of all imaging points is completed, obtaining an imaging matrix I (Lm, lz) of a wingspan curve I 1i corresponding to a two-dimensional wingspan section, and drawing a two-dimensional wingspan section focusing imaging as shown in figure 9;
And step S4.1.3, repeating the steps S4.1.1-S4.1.2, and sequentially performing two-dimensional spanwise cross-section imaging corresponding to the spanwise curve. Collecting all spanwise imaging matrixes to obtain a three-dimensional imaging matrix I (I 1i, lm, lz) of the L 1 -section wing 6;
In step S4.1.4, coordinate transformation is performed to obtain a three-dimensional imaging matrix I (x, y, z) under a rectangular coordinate system, and a three-dimensional focusing image of the L 1 -section wing 6 is drawn as shown in fig. 10.
An example of three-dimensional imaging of the wing 6 based on the wing circumferential scan data set p (J j, ln, t) of the present invention is shown below:
S4.2.1, the shape of the wing circumferential section of the wing 6 is irregular, the section is divided into a plurality of contour sections, the outer contour of the wing 6 is unfolded into a straight line, the normal direction is set as the depth direction, the wing circumferential section of the wing 6 in the imaging section is unfolded into a plane, and frequency domain imaging is carried out in the imaging section according to the method in step S4.1.2, so that plane frequency domain imaging in the contour section is obtained;
Step S4.2.2, repeating step S4.2.1 until the curved surface images of all the contour segments are obtained, and splicing the curved surface images in all the contour segments into a frequency domain imaging matrix I (Ln, lz) of a two-dimensional section of the wing circumference corresponding to the wing circumference curve J j, as shown in FIG. 11;
step S4.2.3, repeating the steps S4.2.1-S4.2.2, and sequentially performing two-dimensional spanwise cross-section imaging corresponding to the airfoil circumference curve. Collecting all the imaging matrixes of the wing circumferential sections to obtain a three-dimensional imaging matrix I (J j, ln, lz) of the wing 6 of the L 1 section;
and S4.2.4, carrying out coordinate conversion to obtain a three-dimensional imaging matrix I (x, y, z) under a rectangular coordinate system, and drawing a three-dimensional focusing image of the wing 6.
Embodiment two:
referring to fig. 12-13, the difference between the present embodiment and the first embodiment is that a rotary clamping mechanism is further provided, the rotary clamping mechanism includes a turntable 14, a movable clamping rail 15 is provided on the turntable 14, the turntable 14 is mounted on the movable rail 16 through a supporting frame, and a driving motor 17 is provided on the supporting frame to drive the turntable 14 to rotate.
In one embodiment of the invention, in use, the left end of the wing 6 is mounted on the connecting disc 3, the wing 6 is rotated to the initial position where the maximum section is parallel to the ground, the bottom moving slide rail 16 is moved, the rotary clamping mechanism is aligned with the right end of the wing 6 in the length direction, the turntable 14 is rotated, the clamping slide rail 15 is slid and moved, the right end of the wing 6 is clamped on the moving clamping slide rail 15, and due to the influence of gravity, the right end has a certain deflection deformation in the mounting of the wing 6, the turntable 14 is adjusted again, and the wing 6 is rotated to the initial position where the maximum section is parallel to the ground.
The wing 6 can be clamped in an auxiliary mode through the arranged rotary clamping mechanism, the defect of curvature is reduced, and therefore detection is more accurate.
In the third embodiment, the difference between the present embodiment and the second embodiment is that the moving clamping slide rail 15 and the corresponding turntable 14 are provided with long holes for extending the right end of the wing 6, and the turntable 14 and the laser ultrasonic probe 8 synchronously move along the axial direction, so that the bending deformation of the detection point under gravity is further reduced, and the detection accuracy is further improved.
In the description of the present invention, it should be understood that the terms "longitudinal," "transverse," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," and the like indicate or are based on the orientation or positional relationship shown in the drawings, merely to facilitate description of the present invention, and do not indicate or imply that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus should not be construed as limiting the present invention.
The foregoing embodiments are merely illustrative of the preferred embodiments of the present invention, and the scope of the present invention is not limited thereto, but various modifications and improvements made by those skilled in the art to which the present invention pertains are made without departing from the spirit of the present invention, and all changes and modifications and improvements fall within the scope of the present invention as defined in the appended claims.