CN108460225A - A kind of emulation mode and system of anisotropy sheet metal induction type magnetic acoustic image - Google Patents
A kind of emulation mode and system of anisotropy sheet metal induction type magnetic acoustic image Download PDFInfo
- Publication number
- CN108460225A CN108460225A CN201810232394.1A CN201810232394A CN108460225A CN 108460225 A CN108460225 A CN 108460225A CN 201810232394 A CN201810232394 A CN 201810232394A CN 108460225 A CN108460225 A CN 108460225A
- Authority
- CN
- China
- Prior art keywords
- surface wave
- wave
- axis
- test specimen
- along
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The present invention discloses a kind of emulation mode and system of anisotropy sheet metal induction type magnetic acoustic image.Method includes：Establish conductivity anisotropy sheet metal simulation architecture model；The conductivity of emulation test specimen different directions is set according to model；Model is placed in magnetostatic field and carries out emulation experiment, obtains inductive loop；According to inductive loop, the static Lorentz force and alternation Lorentz force that test specimen is subject to are determined；Surface wave is obtained according to static Lorentz force and alternation Lorentz force；According to the conductivity in variant direction, velocity of wave of the computational chart surface wave on corresponding direction；According to velocity of wave of the surface wave on different directions, the surface wave displacement on corresponding direction is obtained；According to each surface wave displacement, surface wave Displacements Distribution figure is obtained.Method or system using the present invention by changing the conductivity in all directions, can accurately emulate power distribution and the surface wave displacement of distribution of conductivity anisotropy metal sheet surface.
Description
Technical field
The present invention relates to industrial nondestructive testing fields, more particularly to a kind of anisotropy sheet metal induction type magnetic sound spectrogram
The emulation mode and system of picture.
Background technology
Induction type magnetosonic (Magnetoacoustic Tomography with Magnetic Induction, MAT-MI)
Non-destructive testing technology is a kind of new non-destructive detection method drawn and be tested metal specimen surface conductivity distribution map.Its principle be by
Tested metal specimen is placed in magnetostatic field, introduces alternating magnetic field by coil, the sense with alternating magnetic field same frequency is generated in test specimen
It should be vortexed；Under inductive loop and the collective effect of magnetostatic field, material particle causes week due to the effect by Lorentz force
Phase property is vibrated, and is propagated outward in the form of ultrasonic wave, i.e. magnetoacoustic signals；It is received around test specimen using electromagnet ultrasonic changer
Magnetoacoustic signals are simultaneously sent into computer, and the conductivity spatial distribution of surface of test piece can be reconstructed, and judge lacking for test specimen accordingly
It falls into situation and it is positioned.
Currently, the research for MAT-MI non-destructive testing technologies is built upon measured material distribution of conductivity respectively to same mostly
Under the premise of property.But in a practical situation, metal material is typically anisotropic medium.Conductivity anisotropy is to tested
The distribution of inductive loop density, sound source distribution and acoustic pressure distribution have an impact in part.
Invention content
The object of the present invention is to provide a kind of emulation mode and system of anisotropy sheet metal induction type magnetic acoustic image,
By change all directions on conductivity, accurately emulate distribution of conductivity anisotropy metal sheet surface power distribution and
Surface wave displacement.
To achieve the above object, the present invention provides following schemes：
A kind of emulation mode of anisotropy sheet metal induction type magnetic acoustic image, the emulation mode include：
Conductivity anisotropy sheet metal simulation architecture model is established, the simulation architecture model includes that emulation is tested
Test specimen and emulation broken line coil；
According to the simulation architecture model, the conductivity of the emulation test specimen different directions is set；
The simulation architecture model is placed in magnetostatic field and carries out emulation experiment, obtains inductive loop；
According to the inductive loop, the static Lorentz force and alternation Lorentz force that the test specimen is subject to are determined；
Surface wave is obtained according to the static Lorentz force and the alternation Lorentz force；
According to the conductivity in variant direction, velocity of wave of the surface wave on corresponding direction is calculated；
According to velocity of wave of the surface wave on different directions, the surface wave displacement on corresponding direction is obtained；
According to each surface wave displacement, surface wave Displacements Distribution figure is obtained.
Optionally, described according to the simulation architecture model, the conductivity of the test specimen different directions is set, specifically
Including：
Three-dimensional cartesian rectangular coordinate system is established based on the test specimen, wherein in the surface of the test specimen
The heart is origin, parallel with the test specimen horizontal plane for X-axis, vertical with X-axis for Y-axis, vertical with the faces XOY to be
Z axis；
Conductivity of the test specimen on three X-axis, Y-axis and Z axis directions is set.
Optionally, described that the simulation architecture model is placed in magnetostatic field, calculate the static state that the test specimen is subject to
Lorentz force and alternation Lorentz force, specifically include：
According toObtain the static long-range navigation that test specimen is subject to
Hereby power Fs(h,t)；
Wherein, Fs(h, t) is the static Lorentz force that test specimen is subject to, and σ is conductivity, μ0For space permeability, j is
Imaginary part unit, BsFor static magnetic field strength, z0For the vertical range of broken line coil and test specimen, k is Radial Integrals variable, and I is sharp
The amplitude of electric current is encouraged, ω is the angular frequency of exciting current, J0(k) it is first kind Oth order Bessel function that independent variable is k, J0
(kr) it is first kind Oth order Bessel function that independent variable is kr.
According toObtain test specimen by
The alternation Lorentz force F arrivedd(h,t)；
Wherein, Fd(h, t) is the alternation Lorentz force that test specimen is subject to, μ0For space permeability, I is exciting current
Amplitude, k are Radial Integrals variable, z0For the vertical range of broken line coil and test specimen, j is imaginary part unit, J0(k) it is from change
Amount is the first kind Oth order Bessel function of k, J0(kr) it is first kind Oth order Bessel function that independent variable is kr.
Optionally, described to be specifically included along the velocity of wave of different directions according to surface wave described in the Conductivity Calculation：
The propagation of surface wave is only located at YOZ planes, unrelated with X-axis；
According toObtain the velocity of wave v that the surface wave is propagated along Y-axisRY；
Wherein, vRYFor the velocity of wave that the surface wave is propagated along Y-axis, σyFor the conductivity of Y direction, μmIt is described tested
The Lame constants of part, ρmFor the density of the test specimen；
According toObtain the velocity of wave v that the surface wave is propagated along Z axisRZ；
Wherein, vRZFor the velocity of wave that the surface wave is propagated along Z axis, σzFor the conductivity of Z-direction, μmIt is described tested
The Lame constants of part, ρmFor the density of the test specimen.
Optionally, described that the surface wave displacement is obtained according to the velocity of wave, it specifically includes：
According toCalculate surface wave wave number K along the face that Y direction is propagatedRY；
Wherein, KRYFor the surface wave along the face that Y direction is propagated wave number, ωRYIt is that the surface wave is square along Y-axis
To vibration angular frequency, vRYThe velocity of wave propagated along Y-axis for the surface wave；
According toSurface wave described in t moment is calculated along Y
The in-plane displacement u (y, t) of axis；
Wherein, u (y, t) be test specimen midpoint h at surface wave in moment t along the in-plane displacement of Y-axis, ωRYFor surface
For wave along the vibration angular frequency of Y-axis positive direction, α and β are attenuation coefficient of the surface wave along Z axis positive direction, KTFor shear wave wave number, KRYFor
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen；
According toCalculate surface wave wave number K outside the face that Z-direction is propagatedRZ；
Wherein, KRZFor the surface wave outside the face that Z-direction is propagated wave number, ωRZIt is that the surface wave is square along Z axis
To vibration angular frequency, vRZThe velocity of wave propagated along Z axis for the surface wave；
According toSurface wave described in t moment is calculated along Z
The face outer displacement u (z, t) of axis；
Wherein, u (z, t) be test specimen midpoint h at surface wave in moment t along the face outer displacement of Z axis, ωRZFor surface
For wave along the vibration angular frequency of Z axis positive direction, α and β are surface wave along the attenuation coefficient of Z axis positive direction, KTFor shear wave wave number, KRZIt is
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen.
Optionally, the surface wave Displacements Distribution figure includes surface wave in-plane displacement distribution map and surface corrugated outer displacement point
Butut；
It is described that surface wave Displacements Distribution figure is obtained according to the surface wave displacement, it specifically includes：
According to formula
Meter
Calculate the gray value G at pixel (m, n) in the surface wave in-plane displacement gray-scale mapy(m,n)；
According to each gray value Gy(m, n) obtains surface wave in-plane displacement distribution map Gy；
Wherein, Gy(m, n) is the surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,
T) | y=0,1 ..., W；T=0,1 ..., T } be the surface wave in-plane displacement minimum value, max u (y, t) | y=0,
1,...,W；T=0,1 ..., T } be the surface wave in-plane displacement maximum value, u (y, t) be test specimen midpoint h at table
Surface wave is in moment t along the in-plane displacement of Y-axis；
According to formula
Meter
Calculate the gray value G at pixel (m, n) in the surface corrugated outer displacement gray-scale mapz(m,n)；
According to each gray value Gz(m, n) obtains surface corrugated outer displacement distribution map Gz；
Wherein, Gz(m, n) is the surface corrugated outer displacement gray-scale map GzGray value at middle pixel (m, n)；min{u(z,
T) | z=0,1 ..., H；T=0,1 ..., T } be surface corrugated outer displacement minimum value, max u (z, t) | z=0,
1,...,H；T=0,1 ..., T } be surface corrugated outer displacement maximum value, u (z, t) be test specimen midpoint h at table
Surface wave is in moment t along the face outer displacement of Z axis.
To achieve the above object, the present invention also provides following schemes：
A kind of analogue system of anisotropy sheet metal induction type magnetic acoustic image, the analogue system include：
Simulation architecture model building module, it is described for establishing conductivity anisotropy sheet metal simulation architecture model
Simulation architecture model includes emulation test specimen and emulation broken line coil；
Conductivity setup module, for according to the simulation architecture model, the electricity of the test specimen different directions to be arranged
Conductance；
Inductive loop acquisition module carries out emulation experiment for the simulation architecture model to be placed in magnetostatic field, obtains
Inductive loop；
Lorentz force determining module, for according to the inductive loop, determining the static long-range navigation that the test specimen is subject to
Hereby power and alternation Lorentz force；
Surface wave acquisition module, for obtaining surface wave according to the static Lorentz force and the alternation Lorentz force；
Velocity of wave computing module calculates the surface wave along corresponding direction for the conductivity according to variant direction
Velocity of wave；
Surface wave displacement acquisition module is obtained for the velocity of wave according to surface wave on different directions on corresponding direction
Surface wave displacement；
Distribution map acquisition module, for according to each surface wave displacement, obtaining surface wave Displacements Distribution figure.
Optionally, the conductivity setup module, specifically includes：
Rectangular coordinate system unit establishes three-dimensional cartesian rectangular coordinate system, wherein with institute for being based on the test specimen
The centre of surface for stating test specimen is origin, parallel with the test specimen horizontal plane for X-axis, vertical with X-axis for Y-axis,
It is vertical with the faces XOY for Z axis；
Conductivity setting unit, for conductivity of the test specimen on three X-axis, Y-axis and Z axis directions to be arranged.
Optionally, the surface wave displacement acquisition module, specifically includes：
Wave number acquiring unit in face is used for basisThe surface wave is calculated along the face that Y direction is propagated
Wave number KRY, wherein KRYFor the surface wave along the face that Y direction is propagated wave number, ωRYIt is the surface wave along Y-axis positive direction
Vibration angular frequency, vRYThe velocity of wave propagated along Y-axis for the surface wave；
In-plane displacement acquiring unit is used for basisMeter
Surface wave described in t moment is calculated along the in-plane displacement u (y, t) of Y-axis；
Wherein, u (y, t) be test specimen midpoint h at surface wave in moment t along the in-plane displacement of Y-axis, ωRYFor surface
For wave along the vibration angular frequency of Y-axis positive direction, α and β are attenuation coefficient of the surface wave along Z axis positive direction, KTFor shear wave wave number, KRYFor
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen；
Wave number acquiring unit outside face is used for basisCalculate the face external wave that the surface wave is propagated along Z-direction
Number KRZ, wherein KRZFor the surface wave outside the face that Z-direction is propagated wave number, ωRZIt is the surface wave along Z axis positive direction
Vibrate angular frequency, vRZThe velocity of wave propagated along Z axis for the surface wave；
Face outer displacement acquiring unit is used for basisMeter
Surface wave described in t moment is calculated along the face outer displacement u (z, t) of Z axis；
Wherein, u (z, t) be test specimen midpoint h at surface wave in moment t along the face outer displacement of Z axis, ωRZFor surface
For wave along the vibration angular frequency of Z axis positive direction, α and β are surface wave along the attenuation coefficient of Z axis positive direction, KTFor shear wave wave number, KRZIt is
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen.
Optionally, the surface wave Displacements Distribution figure acquisition module, specifically includes：
First gray value acquiring unit, for according to formula
Meter
Calculate the gray value G at pixel (m, n) in the surface wave in-plane displacement gray-scale mapy(m,n)；
In-plane displacement distribution map acquiring unit, for according to each gray value Gy(m, n) obtains surface wave in-plane displacement
Distribution map Gy；
Wherein, Gy(m, n) is the surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,
T) | y=0,1 ..., W；T=0,1 ..., T } be the surface wave in-plane displacement minimum value, max u (y, t) | y=0,
1,...,W；T=0,1 ..., T } be the surface wave in-plane displacement maximum value, u (y, t) be test specimen midpoint h at table
Surface wave is in moment t along the in-plane displacement of Y-axis；
Second gray value acquiring unit, for according to formula
Calculate the gray value G at pixel (m, n) in the surface corrugated outer displacement gray-scale mapz(m,n)；
Face outer displacement distribution map acquiring unit, for according to each gray value Gz(m, n) obtains surface corrugated outer displacement
Distribution map Gz；
Wherein, Gz(m, n) is the gray value at pixel (m, n) in the surface corrugated outer displacement gray-scale map Gz；min{u
(z, t) | z=0,1 ..., H；T=0,1 ..., T } be surface corrugated outer displacement minimum value, max u (z, t) | z=0,
1,...,H；T=0,1 ..., T } be surface corrugated outer displacement maximum value, u (z, t) be test specimen midpoint h at table
Surface wave is in moment t along the face outer displacement of Z axis.
According to specific embodiment provided by the invention, the invention discloses following technique effects：
The present invention discloses a kind of emulation mode of anisotropy sheet metal induction type magnetic acoustic image, establish conductivity respectively to
Anisotropic sheet metal simulation architecture model；The conductivity of emulation test specimen different directions is set according to model；Model is placed in
Emulation experiment is carried out in magnetostatic field, obtains inductive loop；According to inductive loop, the static Lorentz force that test specimen is subject to is determined
With alternation Lorentz force；Surface wave is obtained according to static Lorentz force and alternation Lorentz force；According to the conductance in variant direction
Rate, velocity of wave of the computational chart surface wave on corresponding direction；According to velocity of wave of the surface wave on different directions, obtain on corresponding direction
Surface wave displacement；According to each surface wave displacement, surface wave Displacements Distribution figure is obtained, method of the invention can be each by changing
Conductivity on direction accurately emulates power distribution and the surface wave displacement of distribution of conductivity anisotropy metal sheet surface.
Description of the drawings
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the present invention
Example, for those of ordinary skill in the art, without having to pay creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is the emulation mode flow chart of anisotropy sheet metal induction type magnetic acoustic image of the embodiment of the present invention；
Fig. 2 is simulation architecture illustraton of model of the embodiment of the present invention；
Fig. 3 is rectangular coordinate system figure of the embodiment of the present invention；
Fig. 4 is the stereogram of test specimen of the embodiment of the present invention and rectangular coordinate system；
Fig. 5 is the analogue system structure chart of anisotropy sheet metal induction type magnetic acoustic image of the embodiment of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of emulation mode and system of anisotropy sheet metal induction type magnetic acoustic image,
By change all directions on conductivity, accurately emulate distribution of conductivity anisotropy metal sheet surface power distribution and
Surface wave displacement.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, below in conjunction with the accompanying drawings and specific real
Applying mode, the present invention is described in further detail.
Fig. 1 is the emulation mode flow chart of anisotropy sheet metal induction type magnetic acoustic image of the embodiment of the present invention.Such as Fig. 1
Shown, a kind of emulation mode of anisotropy sheet metal induction type magnetic acoustic image, the emulation mode includes：
Step 101：Conductivity anisotropy sheet metal simulation architecture model is established, the simulation architecture model includes
Emulate test specimen and emulation broken line coil；
Step 102：According to the simulation architecture model, the conductivity of the test specimen different directions is set；
Step 103：The simulation architecture model is placed in magnetostatic field and carries out emulation experiment, obtains inductive loop；
Step 104：According to the inductive loop, the static Lorentz force and alternation long-range navigation that the test specimen is subject to are determined
Hereby power；
Step 105：Surface wave is obtained according to the static Lorentz force and the alternation Lorentz force；
Step 106：According to the conductivity in variant direction, velocity of wave of the surface wave on corresponding direction is calculated；
Step 107：According to velocity of wave of the surface wave on different directions, the surface wave displacement on corresponding direction is obtained；
Step 108：According to each surface wave displacement, surface wave Displacements Distribution figure is obtained.
Wherein, step 102 specifically includes：
Three-dimensional cartesian rectangular coordinate system is established based on the test specimen, wherein in the surface of the test specimen
The heart is origin, parallel with the test specimen horizontal plane for X-axis, vertical with X-axis for Y-axis, vertical with the faces XOY to be
Z axis；
Conductivity of the test specimen on three X-axis, Y-axis and Z axis directions is set.
Step 103 specifically includes：
According toObtain the static long-range navigation that test specimen is subject to
Hereby power Fs(h,t)；
Wherein, Fs(h, t) is the static Lorentz force that test specimen is subject to, and σ is conductivity, μ0For space permeability, j is
Imaginary part unit, BsFor static magnetic field strength, z0For the vertical range of broken line coil and test specimen, k is Radial Integrals variable, and I is sharp
The amplitude of electric current is encouraged, ω is the angular frequency of exciting current, J0(k) it is first kind Oth order Bessel function that independent variable is k, J0
(kr) it is first kind Oth order Bessel function that independent variable is kr.
According toObtain test specimen by
The alternation Lorentz force F arrivedd(h,t)；
Wherein, FdThe alternation Lorentz force that (h, t) test specimen is subject to, μ0For space permeability, I is the width of exciting current
Value, k are Radial Integrals variable, z0For the vertical range of broken line coil and test specimen, j is imaginary part unit, J0(k) it is independent variable
For the first kind Oth order Bessel function of k, J0(kr) it is first kind Oth order Bessel function that independent variable is kr.
Step 106 specifically includes：
The propagation of surface wave is only located at YOZ planes, unrelated with X-axis；
According toObtain the velocity of wave v that the surface wave is propagated along Y-axisRY；
Wherein, vRYFor the velocity of wave that the surface wave is propagated along Y-axis, σyFor the conductivity of Y direction, μmIt is described tested
The Lame constants of part, ρmFor the density of the test specimen；
According toObtain the velocity of wave v that the surface wave is propagated along Z axisRZ；
Wherein, vRZFor the velocity of wave that the surface wave is propagated along Z axis, σzFor the conductivity of Z-direction, μmIt is described tested
The Lame constants of part, ρmFor the density of the test specimen.
Step 107 specifically includes：
According toCalculate surface wave wave number K along the face that Y direction is propagatedRY；
Wherein, KRYFor the surface wave along the face that Y direction is propagated wave number, ωRYIt is that the surface wave is square along Y-axis
To vibration angular frequency, vRYThe velocity of wave propagated along Y-axis for the surface wave；
According toSurface wave described in t moment is calculated along Y
The in-plane displacement u (y, t) of axis；
Wherein, u (y, t) be test specimen midpoint h at surface wave in moment t along the in-plane displacement of Y-axis, ωRYFor surface
For wave along the vibration angular frequency of Y-axis positive direction, α and β are attenuation coefficient of the surface wave along Z axis positive direction, KTFor shear wave wave number, KRYFor
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen；
According toCalculate surface wave wave number KRZ outside the face that Z-direction is propagated；
Wherein, KRZFor the surface wave outside the face that Z-direction is propagated wave number, ωRZIt is that the surface wave is square along Z axis
To vibration angular frequency, vRZThe velocity of wave propagated along Z axis for the surface wave；
According toSurface wave described in t moment is calculated along Z
The face outer displacement u (z, t) of axis；
Wherein, u (z, t) be test specimen midpoint h at surface wave in moment t along the face outer displacement of Z axis, ωRZFor surface
For wave along the vibration angular frequency of Z axis positive direction, α and β are surface wave along the attenuation coefficient of Z axis positive direction, KTFor shear wave wave number, KRZIt is
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen.
In step 108, the surface wave Displacements Distribution figure includes position outside surface wave in-plane displacement distribution map and surface corrugated
Move distribution map；
According to formula
Calculate the gray value G at pixel (m, n) in the surface wave in-plane displacement gray-scale mapy(m,n)；
According to each gray value Gy(m, n) obtains surface wave in-plane displacement distribution map Gy；
Wherein, Gy(m, n) is the surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,
T) | y=0,1 ..., W；T=0,1 ..., T } be the surface wave in-plane displacement minimum value, max u (y, t) | y=0,
1,...,W；T=0,1 ..., T } be the surface wave in-plane displacement maximum value, u (y, t) be test specimen midpoint h at table
Surface wave is in moment t along the in-plane displacement of Y-axis；
According to formula
Calculate the gray value G at pixel (m, n) in the surface corrugated outer displacement gray-scale mapz(m,n)；
According to each gray value Gz(m, n) obtains surface corrugated outer displacement distribution map Gz；
Wherein, Gz(m, n) is the gray value at pixel (m, n) in the surface corrugated outer displacement gray-scale map Gz；min{u
(z, t) | z=0,1 ..., H；T=0,1 ..., T } be surface corrugated outer displacement minimum value, max u (z, t) | z=0,
1,...,H；T=0,1 ..., T } be surface corrugated outer displacement maximum value, u (z, t) be test specimen midpoint h at table
Surface wave is in moment t along the face outer displacement of Z axis.
Emulation mode provided by the invention, by establishing conductivity anisotropy sheet metal simulation architecture model；According to
The conductivity of model setting emulation test specimen different directions；Model is placed in magnetostatic field and carries out emulation experiment, is incuded
Vortex；According to inductive loop, the static Lorentz force and alternation Lorentz force that test specimen is subject to are determined；According to static Lorentz
Power and alternation Lorentz force obtain surface wave；According to the conductivity in variant direction, wave of the computational chart surface wave on corresponding direction
Speed；According to velocity of wave of the surface wave on different directions, the surface wave displacement on corresponding direction is obtained；According to each surface wave displacement,
Obtain surface wave Displacements Distribution figure, can by change all directions on conductivity, accurately emulate distribution of conductivity respectively to
The power distribution and surface wave displacement of anisotropic metal sheet surface.
Fig. 2 is simulation architecture illustraton of model of the embodiment of the present invention.As shown in Fig. 2, test specimen 1 is cuboid sheet metal,
Its thickness is more than 4 times of ultrasonic wavelength, places it in the magnetostatic field provided by permanent magnet 2 (such as ndfeb magnet), simultaneously will
The broken line coil 3 for being passed through simple sinusoidal alternating current is placed in above test specimen to be measured, and test specimen is sufficiently large relative to coil, to avoid side is generated
Edge effect receives magnetoacoustic signals by electromagnet ultrasonic changer 4 around test specimen.
Fig. 3 is three-dimensional cartesian rectangular coordinate system figure of the embodiment of the present invention.As shown in figure 3, being established to test specimen three-dimensional
Descartes sits straight footmark system OXYZ, and coordinate origin O is located at test specimen centre of surface.The conductivity tensor of test specimen is：
Wherein D is spin matrix：
In formula, γ is that anisotropy moves towards angle, and ψ is anisotropy inclination angle, and χ is anisotropy drift angle.It is test specimen
Reference conductivity rate tensor：
In formula, σx、σy、σzIt is the conductivity on tri- orthogonal directions of X, Y, Z respectively.
Magnetostatic field 5 along Z axis negative direction, intensity B are applied to test specimens.Broken line coil 3 shares n foldings and adjacent leads
Spacing between line is L, and the vertical range between broken line coil 3 and test specimen 1 is z0.Broken line coil 3 is reduced to mutually
Parallel and mutually discrete conducting wire, wherein being passed through the simple sinusoidal alternating current vertical with Z axis as exciting current：
Jc(t)=Iejωt
Wherein, h=(x, y, z) is the position vector of any in three-dimensional cartesian coordinate system OXYZ；T is the time；Jc(t) when being
Between t exciting current intensity；I is the amplitude of exciting current；ω is the angular frequency of exciting current；J is imaginary part unit.Exciting current
Alternating magnetic field is generated around test specimen, the intensity at time t, position h is Bd(h,t)。
The zone of action below coil includes RaAnd RbTwo parts.RaIn point coordinates meet 0 ＜ z ＜ z0, RaIn magnetic arrow
Gesture is mainly caused by the surface eddy of exciting current and test specimen in coil, and calculation formula is
Wherein, Aa(h, t) is RaInterior position h is in the magnetic vector potential of moment t；It is spatial point h in XOY plane
On projection and the distance between coordinate origin；μ is the magnetic conductivity of test specimen；μ0It is space permeability；K is that Radial Integrals become
Amount；J is imaginary part unit；Q is imaginary values；J0() is first kind Oth order Bessel function.
RbIn point coordinates meet z≤0, RbIn magnetic vector potential generated by the surface eddy of test specimen, calculation formula is
Wherein, Ab(h, t) is RbInterior position h is in the magnetic vector potential of moment t.
Single broken line circle at surface of test piece h to be measured moment t generate inductive loop be：
Then n broken lines circle is in the inductive loop of moment t generation at surface of test piece h to be measured：
Je(h, t)=nJ1e(h,t)
Inductive loop JeRespectively static Lorentz force and alternation Lorentz are generated with magnetostatic field and alternating magnetic field interaction
Power, as sound source.Circular is：Je(h, t) and BsInteract the static Lorentz force generated at moment t position h
For：
Je(h, t) and Bd(h, t) at moment t position h axial component interaction generate alternation Lorentz force be：
Dislocation charge in test specimen 1 generates the period under the collective effect of static Lorentz force and alternation Lorentz force
The vibration of property, and propagated outward in the form of surface wave.Surface wave displacement is only located in YOZ planes, unrelated with X.In moment t table
Surface wave is respectively along the in-plane displacement knead dough outer displacement that Y and Z both directions are propagated：
Wherein, u (y, t) is in-plane displacement of the surface wave at 1 midpoint h of test specimen in moment t；U (z, t) is tested
The face outer displacement of surface wave at the h of part midpoint in moment t；It is vibration angular frequency of the surface wave along Y-axis positive direction；It is vibration angular frequency of the surface wave along Z axis positive direction；L is the spacing between adjacent windings；α and β is surface wave edge
The attenuation coefficient of Z axis positive direction：
It is wave number of the surface wave along Y-axis positive direction；It is wave of the surface wave along Z axis positive direction
Number；It is longitudinal wave wave number；It is shear wave wave number；vTIt is transverse wave speed：
vLIt is longitudinal wave velocity：
vRYAnd vRZIt is velocity of wave of the surface wave along Y-axis positive direction and Z axis positive direction respectively：
ρmIt is the density of test specimen；λmAnd μmIt is the Lame constants of test specimen.
By the surface wave in-plane displacement knead dough outer displacement curve at 1 surface each point of test specimen, being converted to grayscale respectively is
256 gray level image GyAnd Gz.The specific method is as follows：
Surface wave in-plane displacement at 1 surface each point of test specimen is u (y, t), wherein y=0,1 ..., W, t=0,
1 ..., T, Fig. 4 are the stereogram of test specimen of the embodiment of the present invention 1 and rectangular coordinate system.As shown in figure 4, W is test specimen
Along the width of Y-direction, T is time of measuring boundary.Then surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n) is
Wherein, Gy(m, n) is surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,t)|y
=0,1 ..., W；T=0,1 ..., T } and max u (y, t) | y=0,1 ..., W；T=0,1 ..., T } it is surface wave respectively
The minimum value and maximum value of in-plane displacement.
Surface corrugated outer displacement at each point of test specimen surface is u (z, t), wherein z=0,1 ..., H, t=0,
1 ..., T, Fig. 4 are the stereogram of test specimen of the embodiment of the present invention and rectangular coordinate system.As shown in figure 4, H is test specimen edge
The thickness of Z-direction, T are time of measuring boundaries.Then the gray value in surface corrugated outer displacement gray-scale map Gz at pixel (m, n) is
Wherein, Gz(m, n) is surface corrugated outer displacement gray-scale map GzGray value at middle pixel (m, n)；min{u(z,t)|z
=0,1 ..., H；T=0,1 ..., T } and max u (z, t) | z=0,1 ..., H；T=0,1 ..., T } it is surface wave respectively
The minimum value and maximum value of face outer displacement.
Fig. 5 is a kind of analogue system structure chart of anisotropy sheet metal induction type magnetic acoustic image of the embodiment of the present invention.
As shown in figure 5, a kind of analogue system of anisotropy sheet metal induction type magnetic acoustic image, the analogue system include：
Simulation architecture model building module 201, for establishing conductivity anisotropy sheet metal simulation architecture model, institute
It includes emulation test specimen and emulation broken line coil to state simulation architecture model；
Conductivity setup module 202, for according to the simulation architecture model, the test specimen different directions to be arranged
Conductivity；
Inductive loop acquisition module 203 carries out emulation experiment for the simulation architecture model to be placed in magnetostatic field, obtains
To inductive loop；
Lorentz force determining module 204, for according to the inductive loop, determining the static Lip river that the test specimen is subject to
Lun Zili and alternation Lorentz force；
Surface wave acquisition module 205, for obtaining surface according to the static Lorentz force and the alternation Lorentz force
Wave；
Velocity of wave computing module 206 calculates the surface wave along corresponding direction for the conductivity according to variant direction
Velocity of wave；
Surface wave displacement acquisition module 207 is obtained for the velocity of wave according to surface wave on different directions on corresponding direction
Surface wave displacement；
Distribution map acquisition module 208, for according to each surface wave displacement, obtaining surface wave Displacements Distribution figure.
Wherein, the conductivity setup module 202, specifically includes：
Rectangular coordinate system unit establishes three-dimensional cartesian rectangular coordinate system, wherein with institute for being based on the test specimen
The centre of surface for stating test specimen is origin, parallel with the test specimen horizontal plane for X-axis, vertical with X-axis for Y-axis,
It is vertical with the faces XOY for Z axis；
Conductivity setting unit, for conductivity of the test specimen on three X-axis, Y-axis and Z axis directions to be arranged.
The surface wave displacement acquisition module 207, specifically includes：
Wave number acquiring unit in face is used for basisThe surface wave is calculated along the face that Y direction is propagated
Wave number KRY, wherein KRYFor the surface wave along the face that Y direction is propagated wave number, ωRYIt is the surface wave along Y-axis positive direction
Vibration angular frequency, vRYThe velocity of wave propagated along Y-axis for the surface wave；
In-plane displacement acquiring unit is used for basisMeter
Surface wave described in t moment is calculated along the in-plane displacement u (y, t) of Y-axis；
Wherein, u (y, t) be test specimen midpoint h at surface wave in moment t along the in-plane displacement of Y-axis, ωRYFor surface
For wave along the vibration angular frequency of Y-axis positive direction, α and β are attenuation coefficient of the surface wave along Z axis positive direction, KTFor shear wave wave number, KRYFor
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen；
Wave number acquiring unit outside face is used for basisCalculate the face external wave that the surface wave is propagated along Z-direction
Number KRZ, wherein KRZFor the surface wave outside the face that Z-direction is propagated wave number, ωRZIt is the surface wave along Z axis positive direction
Vibrate angular frequency, vRZThe velocity of wave propagated along Z axis for the surface wave；
Face outer displacement acquiring unit is used for basisMeter
Surface wave described in t moment is calculated along the face outer displacement u (z, t) of Z axis；
Wherein, u (z, t) be test specimen midpoint h at surface wave in moment t along the face outer displacement of Z axis, ωRZFor surface
For wave along the vibration angular frequency of Z axis positive direction, α and β are surface wave along the attenuation coefficient of Z axis positive direction, KTFor shear wave wave number, KRZIt is
The surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen.
The surface wave Displacements Distribution figure acquisition module 208, specifically includes：
First gray value acquiring unit, for according to formula
Calculate the gray value G at pixel (m, n) in the surface wave in-plane displacement gray-scale mapy(m,n)；
In-plane displacement distribution map acquiring unit, for according to each gray value Gy(m, n) obtains surface wave in-plane displacement
Distribution map Gy；
Wherein, Gy(m, n) is the surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,
T) | y=0,1 ..., W；T=0,1 ..., T } be the surface wave in-plane displacement minimum value, max u (y, t) | y=0,
1,...,W；T=0,1 ..., T } be the surface wave in-plane displacement maximum value, u (y, t) be test specimen midpoint h at table
Surface wave is in moment t along the in-plane displacement of Y-axis；
Second gray value acquiring unit, for according to formula
Calculate the gray value G at pixel (m, n) in the surface corrugated outer displacement gray-scale mapz(m,n)；
Face outer displacement distribution map acquiring unit, for according to each gray value Gz(m, n) obtains surface corrugated outer displacement
Distribution map Gz；
Wherein, Gz(m, n) is the surface corrugated outer displacement gray-scale map GzGray value at middle pixel (m, n)；min{u(z,
T) | z=0,1 ..., H；T=0,1 ..., T } be surface corrugated outer displacement minimum value, max u (z, t) | z=0,
1,...,H；T=0,1 ..., T } be surface corrugated outer displacement maximum value, u (z, t) be test specimen midpoint h at table
Surface wave is in moment t along the face outer displacement of Z axis.
Each embodiment is described by the way of progressive in this specification, the highlights of each of the examples are with other
The difference of embodiment, just to refer each other for identical similar portion between each embodiment.For system disclosed in embodiment
For, since it is corresponded to the methods disclosed in the examples, so description is fairly simple, related place is said referring to method part
It is bright.
Principle and implementation of the present invention are described for specific case used herein, and above example is said
The bright method and its core concept for being merely used to help understand the present invention；Meanwhile for those of ordinary skill in the art, foundation
The thought of the present invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.
Claims (10)
1. a kind of emulation mode of anisotropy sheet metal induction type magnetic acoustic image, which is characterized in that the emulation mode packet
It includes：
Conductivity anisotropy sheet metal simulation architecture model is established, the simulation architecture model includes emulation test specimen
With emulation broken line coil；
According to the simulation architecture model, the conductivity of the emulation test specimen different directions is set；
The simulation architecture model is placed in magnetostatic field and carries out emulation experiment, obtains inductive loop；
According to the inductive loop, the static Lorentz force and alternation Lorentz force that the test specimen is subject to are determined；
Surface wave is obtained according to the static Lorentz force and the alternation Lorentz force；
According to the conductivity in variant direction, velocity of wave of the surface wave on corresponding direction is calculated；
According to velocity of wave of the surface wave on different directions, the surface wave displacement on corresponding direction is obtained；
According to each surface wave displacement, surface wave Displacements Distribution figure is obtained.
2. the emulation mode of anisotropy sheet metal induction type magnetic acoustic image according to claim 1, described according to institute
Simulation architecture model is stated, the conductivity of the test specimen different directions is set, is specifically included：
Three-dimensional cartesian rectangular coordinate system is established based on the test specimen, wherein the centre of surface with the test specimen is
Origin, it is parallel with the test specimen horizontal plane for X-axis, it is vertical with X-axis for Y-axis, it is vertical with the faces XOY for Z axis；
Conductivity of the test specimen on three X-axis, Y-axis and Z axis directions is set.
3. the emulation mode of anisotropy sheet metal induction type magnetic acoustic image according to claim 1, it is described will be described
Simulation architecture model is placed in magnetostatic field, calculates static Lorentz force and alternation Lorentz force that the test specimen is subject to, tool
Body includes：
According toObtain the static Lorentz force F that test specimen is subject tos
(h,t)；
Wherein, Fs(h, t) is the static Lorentz force that test specimen is subject to, and σ is conductivity, μ0For space permeability, j is imaginary part
Unit, BsFor static magnetic field strength, z0For the vertical range of broken line coil and test specimen, k is Radial Integrals variable, and I is excitation electricity
The amplitude of stream, ω are the angular frequency of exciting current, J0(k) it is first kind Oth order Bessel function that independent variable is k, J0(kr) it is
Independent variable is the first kind Oth order Bessel function of kr.
According toObtain what test specimen was subject to
Alternation Lorentz force Fd(h,t)；
Wherein, Fd(h, t) is the alternation Lorentz force that test specimen is subject to, μ0For space permeability, I is the amplitude of exciting current,
K is Radial Integrals variable, z0For the vertical range of broken line coil and test specimen, j is imaginary part unit, J0(k) be independent variable it is k
First kind Oth order Bessel function, J0(kr) it is first kind Oth order Bessel function that independent variable is kr.
4. the emulation mode of anisotropy sheet metal induction type magnetic acoustic image according to claim 2, described according to institute
Surface wave described in Conductivity Calculation is stated along the velocity of wave of different directions, is specifically included：
The propagation of surface wave is only located at YOZ planes, unrelated with X-axis；
According toObtain the velocity of wave v that the surface wave is propagated along Y-axisRY；
Wherein, vRYFor the velocity of wave that the surface wave is propagated along Y-axis, σyFor the conductivity of Y direction, μmFor the test specimen
Lame constants, ρmFor the density of the test specimen；
According toObtain the velocity of wave v that the surface wave is propagated along Z axisRZ；
Wherein, vRZFor the velocity of wave that the surface wave is propagated along Z axis, σzFor the conductivity of Z-direction, μmFor the test specimen
Lame constants, ρmFor the density of the test specimen.
5. the emulation mode of anisotropy sheet metal induction type magnetic acoustic image according to claim 4, described according to institute
It states velocity of wave and obtains the surface wave displacement, specifically include：
According toCalculate surface wave wave number K along the face that Y direction is propagatedRY；
Wherein, KRYFor the surface wave along the face that Y direction is propagated wave number, ωRYFor the surface wave shaking along Y-axis positive direction
Dynamic angular frequency, vRYThe velocity of wave propagated along Y-axis for the surface wave；
According toSurface wave described in t moment is calculated along Y-axis
In-plane displacement u (y, t)；
Wherein, u (y, t) be test specimen midpoint h at surface wave in moment t along the in-plane displacement of Y-axis, ωRYIt is surface wave along Y
The vibration angular frequency of axis positive direction, α and β are attenuation coefficient of the surface wave along Z axis positive direction, KTFor shear wave wave number, KRYIt is described
Surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen；
According toCalculate surface wave wave number K outside the face that Z-direction is propagatedRZ；
Wherein, KRZFor the surface wave outside the face that Z-direction is propagated wave number, ωRZFor the surface wave shaking along Z axis positive direction
Dynamic angular frequency, vRZThe velocity of wave propagated along Z axis for the surface wave；
According toSurface wave described in t moment is calculated along the face of Z axis
Outer displacement u (z, t)；
Wherein, u (z, t) be test specimen midpoint h at surface wave in moment t along the face outer displacement of Z axis, ωRZIt is surface wave along Z
The vibration angular frequency of axis positive direction, α and β are surface wave along the attenuation coefficient of Z axis positive direction, KTFor shear wave wave number, KRZIt is described
Surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen.
6. the emulation mode of anisotropy sheet metal induction type magnetic acoustic image according to claim 4, the surface wave
Displacements Distribution figure includes surface wave in-plane displacement distribution map and surface corrugated outer displacement distribution map；
It is described that surface wave Displacements Distribution figure is obtained according to the surface wave displacement, it specifically includes：
According to formula
Calculate the gray value G at pixel (m, n) in the surface wave in-plane displacement gray-scale mapy(m,n)；
According to each gray value Gy(m, n) obtains surface wave in-plane displacement distribution map Gy；
Wherein, Gy(m, n) is the surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,t)|y
=0,1 ..., W；T=0,1 ..., T } be the surface wave in-plane displacement minimum value, max u (y, t) | y=0,1 ...,
W；T=0,1 ..., T } be the surface wave in-plane displacement maximum value, u (y, t) be test specimen midpoint h at surface wave exist
In-plane displacements of the moment t along Y-axis；
According to formula
Calculate the second gray value G at pixel (m, n) in the surface corrugated outer displacement gray-scale mapz(m,n)；
According to each second gray value Gz(m, n) obtains surface corrugated outer displacement distribution map Gz；
Wherein, Gz(m, n) is the surface corrugated outer displacement gray-scale map GzGray value at middle pixel (m, n)；min{u(z,t)|z
=0,1 ..., H；T=0,1 ..., T } be surface corrugated outer displacement minimum value, max u (z, t) | z=0,1 ...,
H；T=0,1 ..., T } be surface corrugated outer displacement maximum value, u (z, t) be test specimen midpoint h at surface wave exist
Face outer displacements of the moment t along Z axis.
7. a kind of analogue system of anisotropy sheet metal induction type magnetic acoustic image, which is characterized in that the analogue system packet
It includes：
Simulation architecture model building module, for establishing conductivity anisotropy sheet metal simulation architecture model, the emulation
Structural model includes emulation test specimen and emulation broken line coil；
Conductivity setup module, for according to the simulation architecture model, the conductivity of the test specimen different directions to be arranged；
Inductive loop acquisition module carries out emulation experiment for the simulation architecture model to be placed in magnetostatic field, is incuded
Vortex；
Lorentz force determining module, for according to the inductive loop, determining the static Lorentz force that the test specimen is subject to
With alternation Lorentz force；
Surface wave acquisition module, for obtaining surface wave according to the static Lorentz force and the alternation Lorentz force；
Velocity of wave computing module calculates velocity of wave of the surface wave on corresponding direction for the conductivity according to variant direction；
Surface wave displacement acquisition module obtains the surface on corresponding direction for the velocity of wave according to surface wave on different directions
Wave displacement；
Distribution map acquisition module, for according to each surface wave displacement, obtaining surface wave Displacements Distribution figure.
8. the analogue system of anisotropy sheet metal induction type magnetic acoustic image according to claim 7, the conductivity
Setup module specifically includes：
Rectangular coordinate system unit establishes three-dimensional cartesian rectangular coordinate system, wherein with the quilt for being based on the test specimen
The centre of surface of test block is origin, parallel with the test specimen horizontal plane for X-axis, vertical with X-axis for Y-axis, with institute
It is Z axis that it is vertical, which to state the faces XOY,；
Conductivity setting unit, for conductivity of the test specimen on three X-axis, Y-axis and Z axis directions to be arranged.
9. the analogue system of anisotropy sheet metal induction type magnetic acoustic image according to claim 8, the surface wave
Displacement acquisition module, specifically includes：
Wave number acquiring unit in face is used for basisCalculate surface wave wave number along the face that Y direction is propagated
KRY, wherein KRYFor the surface wave along the face that Y direction is propagated wave number, ωRYFor the surface wave shaking along Y-axis positive direction
Dynamic angular frequency, vRYThe velocity of wave propagated along Y-axis for the surface wave；
In-plane displacement acquiring unit is used for basisCalculate t
In-plane displacement u (y, t) of the surface wave described in moment along Y-axis；
Wherein, u (y, t) be test specimen midpoint h at surface wave in moment t along the in-plane displacement of Y-axis, ωRYIt is surface wave along Y
The vibration angular frequency of axis positive direction, α and β are attenuation coefficient of the surface wave along Z axis positive direction, KTFor shear wave wave number, KRYIt is described
Surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen；
Wave number acquiring unit outside face is used for basisCalculate surface wave wave number outside the face that Z-direction is propagated
KRZ, wherein KRZFor the surface wave outside the face that Z-direction is propagated wave number, ωRZFor the surface wave shaking along Z axis positive direction
Dynamic angular frequency, vRZThe velocity of wave propagated along Z axis for the surface wave；
Face outer displacement acquiring unit is used for basisWhen calculating t
Carve face outer displacement u (z, t) of the surface wave along Z axis；
Wherein, u (z, t) be test specimen midpoint h at surface wave in moment t along the face outer displacement of Z axis, ωRZIt is surface wave along Z
The vibration angular frequency of axis positive direction, α and β are surface wave along the attenuation coefficient of Z axis positive direction, KTFor shear wave wave number, KRZIt is described
Surface wave is along the wave number of Y-axis positive direction, μmFor the Lame constants of the test specimen.
10. the analogue system of anisotropy sheet metal induction type magnetic acoustic image according to claim 9, the surface wave
Displacements Distribution figure acquisition module, specifically includes：
First gray value acquiring unit, for according to formula
Calculate the gray value G at pixel (m, n) in the surface wave in-plane displacement gray-scale mapy(m,n)；
In-plane displacement distribution map acquiring unit, for according to each gray value Gy(m, n) obtains surface wave in-plane displacement distribution map
Gy；
Wherein, Gy(m, n) is the surface wave in-plane displacement gray-scale map GyGray value at middle pixel (m, n)；min{u(y,t)|y
=0,1 ..., W；T=0,1 ..., T } be the surface wave in-plane displacement minimum value, max u (y, t) | y=0,1 ...,
W；T=0,1 ..., T } be the surface wave in-plane displacement maximum value, u (y, t) be test specimen midpoint h at surface wave exist
In-plane displacements of the moment t along Y-axis；
Second gray value acquiring unit, for according to formulaMeter
Calculate the gray value G at pixel (m, n) in the surface corrugated outer displacement gray-scale mapz(m,n)；
Face outer displacement distribution map acquiring unit, for according to each gray value Gz(m, n) obtains surface corrugated outer displacement distribution map
Gz；
Wherein, Gz(m, n) is the surface corrugated outer displacement gray-scale map GzGray value at middle pixel (m, n)；min{u(z,t)|z
=0,1 ..., H；T=0,1 ..., T } be surface corrugated outer displacement minimum value, max u (z, t) | z=0,1 ...,
H；T=0,1 ..., T } be surface corrugated outer displacement maximum value, u (z, t) be test specimen midpoint h at surface wave exist
Face outer displacements of the moment t along Z axis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810232394.1A CN108460225B (en) | 2018-03-21 | 2018-03-21 | Simulation method and system of anisotropic metal sheet induction type magnetoacoustic image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810232394.1A CN108460225B (en) | 2018-03-21 | 2018-03-21 | Simulation method and system of anisotropic metal sheet induction type magnetoacoustic image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108460225A true CN108460225A (en) | 2018-08-28 |
| CN108460225B CN108460225B (en) | 2021-08-31 |
Family
ID=63237583
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810232394.1A Active CN108460225B (en) | 2018-03-21 | 2018-03-21 | Simulation method and system of anisotropic metal sheet induction type magnetoacoustic image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108460225B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112287572A (en) * | 2019-07-09 | 2021-01-29 | 中国航发商用航空发动机有限责任公司 | Complex system and method and device for optimizing and verifying lightning stroke direct effect protection of complex system |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE2261500B2 (en) * | 1971-12-17 | 1975-07-03 | Pioneer Electronic Corp., Tokio | Electroacoustic converter |
| CN101791219A (en) * | 2010-03-03 | 2010-08-04 | 中国科学院电工研究所 | Magnetic-acoustic electrical impedance imaging method and device |
| CN102788836A (en) * | 2012-07-26 | 2012-11-21 | 中国科学院电工研究所 | Magneto-acoustic microscopic imaging method and imaging system |
| CN104146708A (en) * | 2014-07-08 | 2014-11-19 | 中国医学科学院生物医学工程研究所 | Magnetoacoustic image reconstruction method and device based on reception property of acoustic transducer |
| CN106023277A (en) * | 2016-05-18 | 2016-10-12 | 华北电力大学（保定） | Magnetic induction magnetoacoustic endoscopic image modeling and simulation method |
| CN106388774A (en) * | 2016-08-28 | 2017-02-15 | 李珣 | Portable sensing type magnetoacoustic two-dimensional conductivity imaging apparatus |
| CN107099695A (en) * | 2017-05-09 | 2017-08-29 | 于泠然 | A kind of deformable permanent-magnet alloy, deformable permanent-magnet alloy strip and preparation method thereof |
-
2018
- 2018-03-21 CN CN201810232394.1A patent/CN108460225B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE2261500B2 (en) * | 1971-12-17 | 1975-07-03 | Pioneer Electronic Corp., Tokio | Electroacoustic converter |
| CN101791219A (en) * | 2010-03-03 | 2010-08-04 | 中国科学院电工研究所 | Magnetic-acoustic electrical impedance imaging method and device |
| CN102788836A (en) * | 2012-07-26 | 2012-11-21 | 中国科学院电工研究所 | Magneto-acoustic microscopic imaging method and imaging system |
| CN104146708A (en) * | 2014-07-08 | 2014-11-19 | 中国医学科学院生物医学工程研究所 | Magnetoacoustic image reconstruction method and device based on reception property of acoustic transducer |
| CN106023277A (en) * | 2016-05-18 | 2016-10-12 | 华北电力大学（保定） | Magnetic induction magnetoacoustic endoscopic image modeling and simulation method |
| CN106388774A (en) * | 2016-08-28 | 2017-02-15 | 李珣 | Portable sensing type magnetoacoustic two-dimensional conductivity imaging apparatus |
| CN107099695A (en) * | 2017-05-09 | 2017-08-29 | 于泠然 | A kind of deformable permanent-magnet alloy, deformable permanent-magnet alloy strip and preparation method thereof |
Non-Patent Citations (1)
| Title |
|---|
| 刘思佳等: "Effect of low-magnitude different-frequency whole-body vibration on subchondral trabecular bone microarchitecture, cartilage degradation, bone/cartilage turnover, and joint pain in rabbits with knee osteoarthritis", 《BMC MUSCULOSKELETAL DISORDERS》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112287572A (en) * | 2019-07-09 | 2021-01-29 | 中国航发商用航空发动机有限责任公司 | Complex system and method and device for optimizing and verifying lightning stroke direct effect protection of complex system |
| CN112287572B (en) * | 2019-07-09 | 2023-10-31 | 中国航发商用航空发动机有限责任公司 | Complex system and lightning stroke direct effect protection optimization and verification method and device thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108460225B (en) | 2021-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pei et al. | A modified meander-line-coil EMAT design for signal amplitude enhancement | |
| Kang et al. | Enhancement of signal amplitude of surface wave EMATs based on 3-D simulation analysis and orthogonal test method | |
| Wang et al. | 3-D modeling and analysis of meander-line-coil surface wave EMATs | |
| Sun et al. | Improvement of unidirectional focusing periodic permanent magnet shear-horizontal wave electromagnetic acoustic transducer by oblique bias magnetic field | |
| KR101061590B1 (en) | Magnetostrictive transducers, structural diagnostic devices and structural diagnostic methods using the same | |
| CN108426948A (en) | A kind of electromagnet ultrasonic changer and its working method of the single mode Lamb wave of excitation | |
| KR101066248B1 (en) | Non-contact type transducer for rod member having multi-loop coil | |
| CN103994923A (en) | Suspension electromagnetic incentive resonant type fatigue test method | |
| Sun et al. | Lamb wave signal selective enhancement by an improved design of meander-coil electromagnetic acoustic transducer | |
| Itaya et al. | Visualization of eddy current distributions for arbitrarily shaped coils parallel to a moving conductor slab | |
| CN112083501A (en) | Seawater ion motion magnetic field simulation method and device for underwater target detection | |
| CN108460225A (en) | A kind of emulation mode and system of anisotropy sheet metal induction type magnetic acoustic image | |
| Wang et al. | Modeling and comparison of three bulk wave EMATs | |
| US20110031966A1 (en) | Non-contact type transducer having multi-loop coil for plate member | |
| Gao et al. | Skin depth and detection ability of magneto-optical imaging for weld defects in alternating magnetic field | |
| Honzík et al. | Acoustic fields in thin fluid layers between vibrating walls and rigid boundaries: integral method | |
| Guan et al. | Accelerometer transverse sensitivity testing with double turntable centrifuge | |
| CN208383812U (en) | A kind of electromagnet ultrasonic changer exciting single mode Lamb wave | |
| Gao et al. | A novel vibration velocity measurement method for electromagnetic vibrators | |
| WO2021039880A1 (en) | Information processing system for nondestructive inspection and nondestructive inspection method | |
| CN114894415A (en) | Method for determining optimal measuring point of sensor based on vibration signal of surface of transformer box body | |
| Wang et al. | A novel method for modeling and analysis of meander-line-coil surface wave EMATs | |
| Boughedda et al. | Electromagnetic acoustic transducer for cracks detection in conductive material | |
| Gebrel et al. | Dynamics of a ring-type macro gyroscope under electromagnetic external actuation forces | |
| Hao et al. | Circuit-field coupled finite element analysis method for an electromagnetic acoustic transducer under pulsed voltage excitation |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |