CN113109447A - Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material - Google Patents

Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material Download PDF

Info

Publication number
CN113109447A
CN113109447A CN202110423348.1A CN202110423348A CN113109447A CN 113109447 A CN113109447 A CN 113109447A CN 202110423348 A CN202110423348 A CN 202110423348A CN 113109447 A CN113109447 A CN 113109447A
Authority
CN
China
Prior art keywords
wave
composite material
dynamic
waves
homogenized
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110423348.1A
Other languages
Chinese (zh)
Inventor
肖海
刘梦龙
李伦
李壮
张耀辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN202110423348.1A priority Critical patent/CN113109447A/en
Publication of CN113109447A publication Critical patent/CN113109447A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654—Imaging
    • G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/0289—Internal structure, e.g. defects, grain size, texture

Abstract

The invention provides a dynamic homogenized ultrasonic full-focusing defect imaging method and a system for a composite material, wherein the method comprises the following steps: step 1: analyzing and modeling the propagation process of sound waves in the periodic layering of the composite laminated board based on a recursive stiffness matrix method; step 2: dynamic homogenization modeling is carried out on the periodic layering units of the composite material laminated plate through a Frokay wave theory, and the homogenization range corresponding to the composite material laminated plate in different layering directions is obtained; and step 3: the energy propagation speed of the quasi-longitudinal wave is obtained according to the dispersion effect of the wave vector in the homogeneous anisotropic material; and 4, step 4: and (3) correcting the arrival time by utilizing the energy propagation speed of the anisotropic composite material in the dynamic homogenization range, so as to perform ultrasonic full-focus defect imaging detection on the composite material laminated plate in various stacking directions. The method and the device realize the visual imaging of the anisotropic composite material of the heterogeneous nature by solving the time compensation problem caused by the sound velocity difference in different propagation directions caused by the composite material layering structure and the center frequency of the excitation signal.

Description

Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material
Technical Field
The invention relates to the field of nondestructive testing, in particular to a dynamic homogenized ultrasonic full-focus defect imaging method and system for a composite material.
Background
The composite material has the advantages of light weight, high specific strength, good corrosion resistance, excellent fatigue durability, strong designability, large-scale integral molding and the like, and is widely applied to the fields of aerospace, automobiles, oceans, biomedicine and the like. As more and more are used in primary load-bearing structures, the thickness of the composite material is also increasing, leading to an increased probability of its defects such as matrix cracking, fiber fracture, interfacial debonding and delamination during manufacturing and service. Nondestructive testing methods for internal defects of materials include X-ray, Computed Tomography (CT), eddy current, thermal imaging, ultrasound, and the like, wherein ultrasound testing has been widely used due to its low cost, safety, no radiation, and convenient testing. The ultrasonic nondestructive detection is to generate ultrasonic pulses in the material, analyze the reflected or transmitted signals of sound waves to obtain defect information, acquire full matrix data through a linear array or area array probe, and perform visual imaging on the defects by combining a full focusing algorithm.
The existing full-focusing imaging algorithm can well detect the defects of isotropic materials, and for anisotropic composite materials, the full-focusing algorithm is applied to defect imaging, so that the arrival time difference caused by the wave speed difference of sound waves propagating along different directions needs to be corrected.
Disclosure of Invention
The invention provides a dynamic homogenized ultrasonic full-focus defect imaging method for a composite material, which comprises the following steps of:
step 1: analyzing and modeling the propagation process of sound waves in the periodic layering of the composite laminated board based on a recursive stiffness matrix method;
step 2: dynamic homogenization modeling is carried out on the periodic layering units of the composite material laminated plate through a Frokay wave theory, and the homogenization range corresponding to the composite material laminated plate in different layering directions is obtained;
and step 3: the energy propagation speed of the quasi-longitudinal wave is obtained according to the dispersion effect of the wave vector in the homogeneous anisotropic material;
and 4, step 4: and (3) correcting the arrival time by utilizing the energy propagation speed of the anisotropic composite material in the dynamic homogenization range, so as to perform ultrasonic full-focus defect imaging detection on the composite material laminated plate in various stacking directions.
As a further improvement of the present invention, in step 1, for a composite material laminated plate structure of an arbitrary ply, a global stiffness matrix of a minimum period unit is obtained according to material parameters and ply information, specifically as follows:
as shown in FIG. 1, the displacement vector u of the j-th layer unit of the composite material laminatejAnd the traction vector σjExpressed as a composite of six sub-waves within the layer:
Figure BDA0003028693770000021
wherein a represents an amplitude, n is 1,2,3 represents a partial wave, ± represents positive and negative directions along a thickness direction z-axis,
Figure BDA0003028693770000022
which represents the vector of the polarization, and,
Figure BDA0003028693770000023
represents wave number, ω represents frequency; definition vector d±nIn relation to wavenumber and polarization vector: component of which
Figure BDA0003028693770000024
Wherein c isijlm(i, j, l, m ═ 1,2,3) is a tensor table of stiffness coefficientsShown in the specification;
the relationship between the displacement and the traction of the upper and lower surfaces of the jth layer is as follows:
Figure BDA0003028693770000025
wherein P is±(3×3)＝[p±1,p±2,p±3]，D±(3×3)＝[d±1,d±2,d±3]，
Figure BDA0003028693770000026
hjIs a single layer thick; layer j stiffness matrix KjComprises the following steps:
Figure BDA0003028693770000027
obtaining global rigidity matrix K of J-layer unit structure by applying recursive algorithmJ：
Figure BDA0003028693770000028
Wherein the content of the first and second substances,
Figure BDA0003028693770000029
Figure BDA00030286937700000210
as a further improvement of the present invention, in step 2, a florkside wave number is obtained by solving a florkside characteristic equation, and then the number of waves stably propagating in the florkside wave, that is, the number of real solutions, is respectively solved by traversing the incident angle and the excitation frequency of the waves, so as to find out the solved homogenization area, which is specifically as follows:
for the smallest unit in a periodically layered composite, the Flokay periodic condition is:
Figure BDA0003028693770000031
where β is the Frokay wavenumber, d is the cell thickness, I (6 × 6) is the identity matrix, subscript 0 denotes the cell upper surface, J denotes the cell lower surface; considering the laminate minimum unit consisting of a J-layer anisotropic monolayer in any direction, its global stiffness matrix Kc(6X 6) layer stiffness matrix KjAccording to the recursive algorithm described above, we obtain:
Figure BDA0003028693770000032
wherein the content of the first and second substances,
Figure BDA0003028693770000033
and (3) combining the two formulas to obtain a Frokay Kve characteristic equation:
Figure BDA0003028693770000034
in the interval-pi < beta d < pi, the Frouka wave characteristic equation has six solutions about Frouka wave number beta, wherein the complex solution represents attenuation waves, the real solution represents stably propagating waves, and the solution of the Frouka wave characteristic equation corresponding to the anisotropic material can be divided into a positive group and a negative group which respectively represent the wave numbers of the Frouka waves propagating along the positive direction and the negative direction of the thickness z axis.
As a further improvement of the present invention, in the step 3, in the homogeneous anisotropic material, the propagation velocities of the acoustic wave in different directions exhibit anisotropy, the actual energy propagation velocity of the acoustic wave can be obtained from a dispersion relation, the dispersion relation between the frequency ω and the wavevector k is represented by G (ω, k) 0 or G (ω, k, θ) 0, and the dominant frequency is represented by ω W (k) or ω W (k, θ), and the energy propagation velocity c is represented by G (ω, k) 0 or G (ω, k), W (k, θ), andeis defined as:
Figure BDA0003028693770000035
as a further improvement of the present invention, in step 4, for any pixel point (x, z), the pixel value I (x, z) is:
Figure BDA0003028693770000036
wherein n is the number of array elements, hijWhen array element i is used as exciting unit, the signal received by array element j is delayed by time delay ti(x, z) and tj(x, z) are the acoustic waves respectively located at xiArray element i at and at xjThe propagation time from the array element j to the pixel point (x, z) is determined by the corresponding distance dt，drEnergy propagation velocity c matched with refraction angle corresponding to sound wave propagation patht，crThe calculation formula is obtained as follows:
Figure BDA0003028693770000041
the invention also discloses a dynamic homogenized ultrasonic full-focus defect imaging system for the composite material, which comprises:
an analytic modeling unit: the recursive stiffness matrix method is used for carrying out analytic modeling on the propagation process of sound waves in the periodic laminated composite material laminated plate.
A homogenization modeling unit: the method is used for carrying out dynamic homogenization modeling on the composite material laminated plate through the Frokay wave theory to obtain the homogenization range corresponding to the composite material laminated plates with different layers.
A solving unit: the method is used for solving the energy propagation speed of the quasi-longitudinal wave according to the dispersion effect of the wave vector in the anisotropic material.
An imaging detection unit: on the basis of a traditional ultrasonic full-focusing algorithm, the arrival time is corrected according to the energy propagation speed curve obtained by the solving unit, so that ultrasonic full-focusing defect imaging detection of the composite material laminated plate in various laminating directions is realized.
In the analysis modeling unit, for a composite material laminated plate structure with any layers, the material parameters and the layer information of single-layer fibers are known, and a global stiffness matrix of a minimum period unit is obtained through a recursive stiffness matrix method and a Flokay wave periodic condition.
In the homogenization modeling unit, the Frokay wave number is obtained by solving a Frokay characteristic equation, then the incident angle and the excitation frequency of the wave are traversed, the number of the waves which are stably propagated in the Frokay wave, namely the number of real number solutions, is respectively solved, and the solved homogenization region is found out.
In the solving unit, in a homogeneous structure after dynamic homogeneous equivalence, the propagation speeds of the sound waves in different directions are anisotropic, the energy propagation speeds in all directions can be obtained through solving of the dispersion effect of wave vectors, and the propagation time of the sound waves on different paths can be obtained according to the energy propagation speeds.
In the imaging detection unit, the propagation time of the sound wave from the excitation unit to the pixel point and from the pixel point to the two-section path of the receiving unit is solved based on the energy propagation speed curve, namely on the basis of the traditional ultrasonic full-focusing imaging algorithm, the propagation time is corrected according to the energy propagation speed curve obtained by dynamic homogenization for the composite materials in different layering directions, and therefore the purpose of improving the full-focusing imaging effect of the heterogeneous composite materials is achieved.
The invention has the beneficial effects that: the invention discloses a dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite materials, which are characterized in that a dynamic homogenized model is constructed by combining a recursive stiffness matrix method and a Flokay wave theory, sound velocities in different propagation directions under a composite material layer structure and excitation signal center frequency are obtained, arrival time differences caused by the sound velocity differences are corrected in a full-focusing algorithm, and visual imaging of damage in a heterogeneous anisotropic composite material is realized.
Drawings
FIG. 1 is a schematic view of a composite laminate panel of the present invention;
FIG. 2 is a schematic representation of a composite laminate coordinate system of the present invention;
FIG. 3 is a plot of the passband and stopband of the Frokay-Kai wave of the present invention (white: 2 stably propagating waves, passband; gray: 1 stably propagating wave; black: attenuated only wave, stopband);
FIG. 4 is a schematic anisotropy plot of phase velocity and energy propagation velocity of acoustic waves of different frequencies propagating in [0/90] ply structure composite material according to the present invention;
FIG. 5 is a flow chart of the dynamic homogenized ultrasonic full focus defect imaging method for composite materials of the present invention.
Detailed Description
As shown in FIG. 5, the invention discloses a dynamic homogenized ultrasonic full-focus defect imaging method and system for composite materials, which realize the visual imaging of anisotropic composite materials through solving the time compensation problem caused by the sound velocity difference in different propagation directions caused by the layering structure of the composite materials and the center frequency of an excitation signal.
Generally, composite materials will take a periodic layup pattern to achieve the desired optimum material properties, while different layup patterns will cause the composite material laminate to become a heterogeneous structure, resulting in complex acoustic reflection and transmission phenomena between the layers. According to the Flokay wave state homogenization theory, from the wave propagation angle, the heterogeneous periodic structure can be regarded as a homogeneous structure in a certain range of the excitation frequency and the propagation angle.
1. Analytical modeling of acoustic wave propagation in periodically laid-up composite laminates
And (3) obtaining a global rigidity matrix by applying a recursive algorithm according to the layer rigidity matrix:
Figure BDA0003028693770000051
the global stiffness matrix can be used to describe the displacement-traction relationship of the upper and lower surfaces of the overall composite laminate. Given the material coefficients of the single layer fibers and the lay-up direction of the laminate, a global stiffness matrix of the composite laminate can be determined by a recursive stiffness matrix method.
2. Solving the Flokay characteristic equation
Figure BDA0003028693770000061
Within the interval- π < β d < π, the equation has six solutions for Froukeley wavenumbers β, where the complex solution represents an evanescent wave and the real solution represents a stably propagating wave. The solution of the anisotropic material corresponding to the characteristic equation can be divided into a positive group and a negative group, which respectively represent the wave numbers of the Frokay waves propagating along the positive direction and the negative direction of the z axis.
3. Solving for dynamic homogenization range
For a composite material laminated plate structure with any laying layer, knowing material parameters of a single layer and laying layer information of the composite material laminated plate structure, the global stiffness matrix of a minimum period unit can be obtained, then solving a Fraoki wave characteristic equation to obtain the number of Fraoki waves, traversing the incidence angle and excitation frequency of the waves, and respectively solving the number of the waves which are stably propagated in the Fraoki waves, namely the number of real number solutions, so that the solved homogenization area can be found out.
For example, the ply orientation is [0/90]]48The minimum periodic unit of the composite material with the thickness of 125um of a single layer, the pass band and the stop band of the Frokay wave are shown in FIG. 3 (considering the x-z plane):
the area shown in fig. 3 can be divided into three parts:
(1) a homogenization region; the incident angle (lower than theta) corresponding to the regioncr1) And at excitation frequencies (below 3MHz), the florokay waves in the minimum period unit of the composite material are all stably propagating waves, so that the composite materials of different ply structures can be regarded as homogeneous materials in corresponding homogenization regions.
(2) A partial homogenization region; corresponding angle of incidence (above theta)cr1But less than thetacr2) And 1 stably propagating wave and 1 evanescent wave in the excitation frequency range, and can therefore also be considered as a homogenized region.
(3) A non-homogenizeable region; there is no stably propagating wave in this region and it cannot be homogenized.
Through dynamic homogenization, the composite material in any layering direction can be regarded as a homogeneous structure within a certain incidence angle and excitation frequency range, so that the propagation process of sound waves in the anisotropic material is simplified, and great convenience is provided for subsequent full-focus imaging.
4. Solving an energy propagation velocity curve for a homogeneous anisotropic structure
In the homogeneous anisotropic material, the energy propagation speeds of the sound waves along different directions are obtained according to the dispersion effect of the wave vector. The phase velocity and energy propagation velocity of sound propagation in the composite material laminate having a laminate structure of [0/90] and a thickness of 125um can be determined by the above-described method, as shown in fig. 4.
5. Ultrasound full focus imaging
The propagation time of a path of an excitation unit, an imaging point and a receiving unit is corrected by utilizing an energy propagation speed curve of the anisotropic composite material in a dynamic homogenization range, so that ultrasonic full-focus defect imaging detection is performed on the composite material laminated plates in various stacking directions.
In summary, for the anisotropic composite laminate, the invention firstly analyzes and models the propagation process of sound waves in the laminate based on a recursive stiffness matrix method, then dynamically homogenizes and models the composite laminate by the Froky wave theory, so as to obtain the homogenizeable range corresponding to the composite materials of different layers, then obtains the propagation velocity of quasi-longitudinal waves according to the dispersion effect of wave vectors in the homogeneous anisotropic materials, and combines the full matrix capture data obtained by experiments to realize the visual detection of the anisotropic composite laminate by improving the full focus imaging algorithm.
The invention has the beneficial effects that:
the invention discloses a dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite materials, which are characterized in that a dynamic homogenized model is constructed by combining a recursive stiffness matrix method and a Flokay wave theory, sound velocities in different propagation directions under a composite material layer structure and excitation signal center frequency are obtained, arrival time differences caused by the sound velocity differences are corrected in a full-focusing algorithm, and visual imaging of damage in a heterogeneous anisotropic composite material is realized.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (10)

1. A dynamic homogenized ultrasonic full-focus defect imaging method for composite materials is characterized by comprising the following steps of:
step 1: analyzing and modeling the propagation process of sound waves in the periodic layering of the composite laminated board based on a recursive stiffness matrix method;
step 2: dynamic homogenization modeling is carried out on the periodic layering units of the composite material laminated plate through a Frokay wave theory, and the homogenization range corresponding to the composite material laminated plate in different layering directions is obtained;
and step 3: the energy propagation speed of the quasi-longitudinal wave is obtained according to the dispersion effect of the wave vector in the homogeneous anisotropic material;
and 4, step 4: and (3) correcting the arrival time by utilizing the energy propagation speed of the anisotropic composite material in the dynamic homogenization range, so as to perform ultrasonic full-focus defect imaging detection on the composite material laminated plate in various stacking directions.
2. The method for imaging defects of composite material through dynamic homogenized ultrasonic full focus according to claim 1, wherein in the step 1, for a composite material laminated plate structure of any ply, a global stiffness matrix of a minimum period unit is obtained according to material parameters and ply information, and the method is as follows:
as shown in FIG. 1, the displacement vector u of the j-th layer unit of the composite material laminatejAnd the traction vector σjExpressed as a composite of six sub-waves within the layer:
Figure FDA0003028693760000011
wherein a represents an amplitude, n is 1,2,3 represents a partial wave, ± represents positive and negative directions along a thickness direction z-axis,
Figure FDA0003028693760000012
which represents the vector of the polarization, and,
Figure FDA0003028693760000013
represents wave number, ω represents frequency; definition vector d±nIn relation to wavenumber and polarization vector: component of which
Figure FDA0003028693760000014
Wherein c isijlm(i, j, l, m ═ 1,2,3) is a tensor representation of the stiffness coefficients;
the relationship between the displacement and the traction of the upper and lower surfaces of the jth layer is as follows:
Figure FDA0003028693760000015
wherein P is±(3×3)＝[p±1,p±2,p±3]，D±(3×3)＝[d±1,d±2,d±3]，
Figure FDA0003028693760000016
hjIs a single layer thick; layer j stiffness matrix KjComprises the following steps:
Figure FDA0003028693760000021
obtaining global rigidity matrix K of J-layer unit structure by applying recursive algorithmJ：
Figure FDA0003028693760000022
Wherein the content of the first and second substances,
Figure FDA0003028693760000023
Figure FDA0003028693760000024
3. the method for imaging defects of dynamic homogenized ultrasonic full focus of composite materials according to claim 1, wherein in the step 2, the number of Frokay waves is obtained by solving a Frokay characteristic equation, then the number of waves stably propagating in the Frokay waves, namely the number of real solutions, is respectively solved by traversing the incidence angle and the excitation frequency of the waves, so as to find out the homogenization region, which is specifically as follows:
for the smallest unit in a periodically layered composite, the Flokay periodic condition is:
Figure FDA0003028693760000025
where β is the Frokay wavenumber, d is the cell thickness, I (6 × 6) is the identity matrix, subscript 0 denotes the cell upper surface, J denotes the cell lower surface; considering the laminate minimum unit consisting of a J-layer anisotropic monolayer in any direction, its global stiffness matrix Kc(6X 6) layer stiffness matrix KjAccording to the recursive algorithm described above, we obtain:
Figure FDA0003028693760000026
wherein the content of the first and second substances,
Figure FDA0003028693760000027
and (3) combining the two formulas to obtain a Frokay Kve characteristic equation:
Figure FDA0003028693760000028
in the interval-pi < beta d < pi, the Frouka wave characteristic equation has six solutions about Frouka wave number beta, wherein the complex solution represents attenuation waves, the real solution represents stably propagating waves, and the solution of the Frouka wave characteristic equation corresponding to the anisotropic material can be divided into a positive group and a negative group which respectively represent the wave numbers of the Frouka waves propagating along the positive direction and the negative direction of the thickness z axis.
4. The method for imaging defects of dynamic homogenized ultrasonic full focus of composite material as claimed in claim 1, wherein in said step 3, the propagation velocity of sound wave in different directions in the homogeneous anisotropic material is anisotropic, and the actual energy propagation velocity can be obtained by wave vector dispersion relation; the dispersion relation between the frequency ω and the wave vector k is represented by G (ω, k) 0 or G (ω, k, θ) 0, the dominant is represented by ω (W) (k) or ω (W, θ), and the energy propagation velocity ceIs defined as:
Figure FDA0003028693760000031
5. the method for imaging defects of full focus of dynamically homogenized ultrasound of composite material according to claim 1, wherein in step 4, for any pixel point (x, z), its pixel value I (x, z) is:
Figure FDA0003028693760000032
wherein n is the number of array elements, hijWhen array element i is used as exciting unit, the signal received by array element j is delayed by time delay ti(x, z) and tj(x, z) are the acoustic waves respectively located at xiArray element i at and at xjThe propagation time from the array element j to the pixel point (x, z) is determined by the corresponding distance dt，drEnergy propagation velocity c matched with refraction angle corresponding to sound wave propagation patht，crThe calculation formula is obtained as follows:
Figure FDA0003028693760000033
6. a dynamic homogenized ultrasonic full focus defect imaging system for composite materials, comprising:
an analytic modeling unit: the method is used for analyzing and modeling the propagation process of sound waves in the periodic laminated composite material laminated plate by a recursive stiffness matrix method;
a homogenization modeling unit: the method is used for carrying out dynamic homogenization modeling on the composite material laminated plate through a Frokay wave theory to obtain the homogenization range corresponding to the composite material laminated plates with different layers;
a solving unit: the energy propagation speed of the quasi-longitudinal wave is obtained according to the dispersion effect of the wave vector in the homogeneous anisotropic material;
an imaging detection unit: on the basis of a traditional ultrasonic full-focusing algorithm, the arrival time is corrected according to the energy propagation speed curve obtained by the solving unit, so that ultrasonic full-focusing defect imaging detection of the composite material laminated plate in various laminating directions is realized.
7. The system for imaging defects of dynamic homogenized ultrasonic full focus of composite materials according to claim 6, wherein in the analytical modeling unit, for a composite material laminated plate structure of any ply, the material parameters of single-layer fibers and ply information are known, and a global stiffness matrix of a minimum period unit can be obtained through a recursive stiffness matrix method and a Flokay wave periodicity condition.
8. The dynamic homogenized ultrasonic full-focus defect imaging system of composite materials according to claim 6, wherein in the homogenization modeling unit, a Frokay wave number is obtained by solving a Frokay characteristic equation, and then the number of waves stably propagating in the Frokay wave, namely the number of real solution, is respectively solved by traversing the incidence angle and the excitation frequency of the waves, so as to find out the solved homogenization area.
9. The system for imaging defects of dynamic homogenized ultrasonic full focus of composite materials according to claim 6, wherein in the solving unit, in the homogenized structure after the dynamic homogenization equivalence, the energy propagation speed can be obtained by solving through the dispersion effect of wave vectors, and the propagation time of the sound waves on different paths can be obtained according to the energy propagation speed.
10. The dynamic homogenized ultrasonic full-focus defect imaging system of composite materials according to claim 6, wherein in the imaging detection unit, the propagation time of the sound wave from the excitation unit to the pixel point and from the pixel point to the two paths of the receiving unit is solved based on the energy propagation speed curve calculated in claim 9, that is, on the basis of a traditional ultrasonic full-focus imaging algorithm, the propagation time is corrected according to the energy propagation speed curve obtained by dynamic homogenization for composite materials in different layering directions, so as to achieve the purpose of improving the full-focus imaging effect of the heterogeneous composite materials.
CN202110423348.1A 2021-04-20 2021-04-20 Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material Pending CN113109447A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110423348.1A CN113109447A (en) 2021-04-20 2021-04-20 Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110423348.1A CN113109447A (en) 2021-04-20 2021-04-20 Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material

Publications (1)

Publication Number Publication Date
CN113109447A true CN113109447A (en) 2021-07-13

Family

ID=76718662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110423348.1A Pending CN113109447A (en) 2021-04-20 2021-04-20 Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material

Country Status (1)

Country Link
CN (1) CN113109447A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114384152A (en) * 2022-01-13 2022-04-22 山东大学 Ultrasonic guided wave damage positioning method and system based on search point matching
CN114428119A (en) * 2022-01-20 2022-05-03 重庆大学 Method for inverting elastic constant of composite material with anisotropic characteristics
CN116577417A (en) * 2023-07-13 2023-08-11 浙江大学 Self-adaptive ultrasonic full-focusing defect imaging method for composite material
CN117147706A (en) * 2023-11-01 2023-12-01 浙江大学 Composite material anisotropic imaging method based on quasi-longitudinal wave group velocity

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085038A (en) * 2017-05-08 2017-08-22 北京工业大学 The total focus imaging method of the amendment of anisotropy weld defect array detection
CN107219305A (en) * 2017-06-02 2017-09-29 北京航空航天大学 A kind of total focus imaging detection method based on annular array transducer
CN111610257A (en) * 2020-05-18 2020-09-01 武汉大学 Array ultrasonic focusing imaging correction method for metal additive manufacturing heterogeneous tissue
CN113552219A (en) * 2021-07-28 2021-10-26 浙江大学 Ultrasonic self-focusing detection method for hole defects of multilayer structure

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085038A (en) * 2017-05-08 2017-08-22 北京工业大学 The total focus imaging method of the amendment of anisotropy weld defect array detection
CN107219305A (en) * 2017-06-02 2017-09-29 北京航空航天大学 A kind of total focus imaging detection method based on annular array transducer
CN111610257A (en) * 2020-05-18 2020-09-01 武汉大学 Array ultrasonic focusing imaging correction method for metal additive manufacturing heterogeneous tissue
CN113552219A (en) * 2021-07-28 2021-10-26 浙江大学 Ultrasonic self-focusing detection method for hole defects of multilayer structure

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
L. WANG ET AL.: "Analysis of Ultrasonic Wave Propagation in Multiply Composites: Homogenization and Effective Anisotropic Media", 《REVIEW OF PROGRESS IN QUANTITATIVE NONDESTRUCTIVE EVALUATION》 *
L. WANG ET AL.: "Floquet wave homogenization of periodic anisotropic media", 《THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA》 *
LI LIN ET AL.: "Total focusing method imaging of multidirectional CFRP laminate with model-based time delay correction", 《NDT AND E INTERNATIONAL》 *
MENGLONG LIU ET AL.: "Floquet wave theory-based time-corrected ultrasonic total focusing method for fiber-reinforced composite laminate", 《ULTRASONICS》 *
S.I. ROKHLIN ET AL.: "Ultrasonic waves in layered anisotropic media: characterization of multidirectional composites", 《INTERNATIONAL JOURNAL OF SOLIDS AND STRUCTURES》 *
葛瑞•马沃可 等: "《岩石物理手册 孔隙介质中地震分析工具》", 31 March 2008 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114384152A (en) * 2022-01-13 2022-04-22 山东大学 Ultrasonic guided wave damage positioning method and system based on search point matching
CN114384152B (en) * 2022-01-13 2023-09-01 山东大学 Ultrasonic guided wave damage positioning method and system based on search point matching
CN114428119A (en) * 2022-01-20 2022-05-03 重庆大学 Method for inverting elastic constant of composite material with anisotropic characteristics
CN114428119B (en) * 2022-01-20 2023-11-21 重庆大学 Method for inverting elastic constant of composite material with anisotropic characteristic
CN116577417A (en) * 2023-07-13 2023-08-11 浙江大学 Self-adaptive ultrasonic full-focusing defect imaging method for composite material
CN116577417B (en) * 2023-07-13 2023-09-29 浙江大学 Self-adaptive ultrasonic full-focusing defect imaging method for composite material
CN117147706A (en) * 2023-11-01 2023-12-01 浙江大学 Composite material anisotropic imaging method based on quasi-longitudinal wave group velocity
CN117147706B (en) * 2023-11-01 2024-01-30 浙江大学 Composite material anisotropic imaging method based on quasi-longitudinal wave group velocity

Similar Documents

Publication Publication Date Title
CN113109447A (en) Dynamic homogenized ultrasonic full-focusing defect imaging method and system for composite material
Yu et al. Damage detection in quasi-isotropic composite bends using ultrasonic feature guided waves
Ramadas et al. Interaction of Lamb mode (Ao) with structural discontinuity and generation of “Turning modes” in a T-joint
Asokkumar et al. Comparison of ultrasonic non-contact air-coupled techniques for characterization of impact-type defects in pultruded GFRP composites
Yu et al. Lamb wave-based damage imaging of CFRP composite structures using autoencoder and delay-and-sum
Jodhani et al. Ultrasonic non-destructive evaluation of composites: A review
Ramadas et al. Numerical and experimental studies on propagation of A0 mode in a composite plate containing semi-infinite delamination: Observation of turning modes
Zhang et al. Visualized characterization of diversified defects in thick aerospace composites using ultrasonic B-scan
Wang et al. A modal decomposition imaging algorithm for ultrasonic detection of delamination defects in carbon fiber composite plates using air-coupled Lamb waves
Jakubczak et al. Non-destructive testing investigation of gaps in thin Glare laminates
Yang et al. Ultrasonic detection methods for mechanical characterization and damage diagnosis of advanced composite materials: A review
Liao et al. Super-resolution ultrasound lamb wave NDE imaging of anisotropic airplane laminates via deconvolutional neural network
Smagulova et al. Novel processing algorithm to improve detectability of disbonds in adhesive dissimilar material joints
Ishii et al. Transmission of ultrasonic waves at oblique incidence to composite laminates with spring-type interlayer interfaces
Zhang et al. Through transmission ultrasonic inspection of fiber waviness for thickness-tapered composites using ultrasound non-reciprocity: Simulation and experiment
Tian et al. Optimized ultrasonic total focusing imaging of diverse and multiple defects in crossply CFRP: Floquet wave theory, numerical simulation, and experimental validation
Chong et al. Full-field ultrasonic inspection for a composite sandwich plate skin-core debonding detection using laser-based ultrasonics
Luo et al. Rapid damage reconstruction imaging of composite plates using non-contact air-coupled Lamb waves
Tong et al. Disbond contour estimation in aluminum/CFRP adhesive joint based on the phase velocity variation of Lamb waves
CN106248009A (en) A kind of ultrasonic thickness-measuring method of highly attenuating fiberglass
Wang et al. Ultrasonic wave interaction with multidirectional composites: Modeling and experiment
Liu Research progress of acoustic emission detection technology based on modal theory
Padiyar M et al. Quantitative characterization of interface delamination in composite T-joint using couplant-free Lamb wave methods
Zou et al. Study on ultrasonic testing of interface bonding state of aluminum matrix composites
Zhong et al. Ultrasonic Testing Techniques for Nondestructive Evaluation of Fiber-Reinforced Composite Structures

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210713