Summary of the invention
The object of the present invention is to provide the pulse eddy current detection methods and detection device of a kind of ferromagnetic metallic material conductivity.
On the basis of the time domain analytic solutions of ferromagnetic plate Pulsed eddy current testing model induced voltage, using induced voltage experiment curv,
The least square problem established between induced voltage time-domain signal measured value and calculated value carrys out inverting and is detected ferromagnetic component
Conductivity and magnetic conductivity;And the amplitude by being stepped up pulse excitation electric current, it obtains electric under different amplitude pulse current excitations
The change curve of conductance inversion result;Finally by the method for curve matching, find the change curve of conductivity inversion result with
The intersection point of axis of ordinates, the conductivity of as tested ferromagnetic component.
Present invention firstly provides a kind of pulse eddy current detection methods of ferromagnetic metallic material conductivity, the specific steps are as follows:
Step 1: the acquisition of the Pulsed eddy current testing signal of tested ferromagnetic component.
Step 2: the inversion method of tested ferromagnetic component electromagnetic parameter.
Step 3: eliminating influence of the magnetic conductivity variation to conductivity testing result using curve matching.
The method advantage that the present invention implements Pulsed eddy current testing to ferromagnetic metallic material conductivity is:
(1) a kind of non-contact, lossless detection method can be provided for the conductivity of ferromagnetic metallic material in industry, filled up
Such technological gap.Can indirectly judgement material hardness, condition of heat treatment, residual stress etc., help to improve material processing work
Skill is horizontal, provides a kind of new nondestructiving detecting means for links such as the productions, maintenance and monitoring of key position ferromagnetic component.
(2) decoupling between conductivity and magnetic conductivity is realized, impulse eddy current method is made to detect ferromagnetic metallic material conductivity
It is possibly realized.Theoretically it can be proved that conductivity can be coupled with magnetic conductivity when eddy current signal can not penetrate tested ferromagnetic component
Together, using eddy current testing signal, it can only be accurately finally inversed by the ratio of the two, can not detect conductivity.In the present invention,
Using impulse eddy current method, eddy current signal can be made effectively to penetrate tested ferromagnetic component, it, can be accurately anti-using time domain induced voltage signal
The product and relative permeability of tested ferromagnetic component conductivity and wall thickness and product two amounts of wall thickness are performed, to realize
The decoupling of conductivity and magnetic conductivity makes it possible that impulse eddy current method detects ferromagnetic metallic material conductivity.
(3) influence of the magnetic conductivity nonlinear change to conductivity testing result is weakened, conductivity testing result is improved
Precision and reliability.For ferromagnetic metallic material, the magnetic conductivity of material is difficult directly to measure, and the size of magnetic conductivity by
The factors such as tested ferromagnetic material microstructure, remanent magnetism, pulse excitation magnetic field strength influence, and will affect conductivity testing result.?
In the present invention, the magnetic conductivity of tested ferromagnetic component is on the one hand set as a unknown parameter, at the signal proposed through the invention
Reason method carrys out the value that inverting determines detection magnetic conductivity every time;On the other hand, by electric under different amplitude pulse current excitations
Conductance inversion result carries out curve fitting, and finds corresponding conductivity inversion result when external magnetic excitation field is zero and ties as detection
Fruit, to weaken the influence caused by testing result of magnetic conductivity nonlinear change, improve conductivity detection precision and can
By property.
Specific embodiment
Below in conjunction with drawings and examples, the present invention is described in further detail.
The present invention provides the pulse eddy current detection method and detection device of a kind of ferromagnetic metallic material conductivity, the detection
Device is as shown in Figure 1, include coil probe, host, DA analog-digital converter, power amplification circuit, sampling resistor, the first AD modulus
Converter and the 2nd AD analog-digital converter, the coil probe includes coil rack, excitation coil and detection coil, and coil is visited
Head (as shown in Figure 1A, 1 is coil rack, and 2 be excitation coil, and 3 be detection coil) is placed in above tested ferromagnetic component, and coil is visited
Head the distance between lower edge and tested ferromagnetic component upper surface (referred to as coil probe lift off) are denoted as lo, under coil is popped one's head in
The wall thickness of the tested ferromagnetic component in side is denoted as d.The coil rack uses air-cored cylindrical coil skeleton, successively twines on coil rack
Around excitation coil and detection coil, wherein the position of excitation coil and detection coil be can be interchanged.The output end of host connects DA
The input terminal of digital analog converter, the input terminal of the output end connection power amplification circuit of DA digital analog converter, power amplification circuit
Output end connection sampling resistor input terminal, the second output terminal of sampling resistor connects the input of the 2nd AD analog-digital converter
End, the second input terminal of the output end connection host of the 2nd AD analog-digital converter, the first output end of sampling resistor connect excitation
The both ends at the both ends of coil, detection coil connect the first AD analog-digital converter input terminal, the output end of the first AD analog-digital converter
The first input end of host is connected, wherein host will connect to DA digital analog converter output drive digital signal, DA digital analog converter
The excitation digital signal of receipts is converted to excitation analog signal and is transferred to power amplification circuit, and power amplification circuit is by received excitation
Pulse excitation electric current is exported after analog signal amplifying power to sampling resistor, supply swashs after sampling resistor receives pulse excitation electric current
Coil is encouraged, detection coil will connect induced voltage analog signal transmission to the first AD analog-digital converter, the first AD analog-digital converter
The induced voltage analog signal of receipts is converted to induced voltage discrete digital signal and is transferred to host progress as the result is shown.
A kind of pulse eddy current detection method of ferromagnetic metallic material conductivity provided by the invention, includes the following steps:
Step 1: the acquisition (SAP) of the Pulsed eddy current testing signal of tested ferromagnetic component:
The obtaining step of pulse eddy current testing signal in step 1 of the present invention are as follows:
Step SAP-1, by coil probe (as shown in Figure 1A, 1 is coil rack, and 2 be excitation coil, and 3 be detection coil)
It is disposed vertically above tested ferromagnetic component, hold-in winding probe lift off lo(0-20mm) is constant, as shown in Figure 1;
Step SAP-2, excitation coil and detection coil are wound on coil rack, and the both ends of excitation coil connect sampling resistor
The first output end, sampling resistor input termination power amplification circuit output end, detection coil two termination the first AD moulds
The input terminal of number converter, the second output terminal of current sampling resistor connect the input terminal of the 2nd AD analog-digital converter;
Step SAP-3 is exported with host programming and is continued pulsewidth as 10~500ms, and the excitation number that amplitude is 0.1~1V is believed
Number;After DA digital analog converter, become to continue pulsewidth to be 10~500ms, amplitude is the excitation analog signal of 0.1~1V, output
To power amplification circuit;After power amplification circuit amplifying power, exports and continue pulsewidth as 10~500ms, amplitude is 0.1~5A
Pulse excitation electric current, after pulse excitation electric current flows through sampling resistor, supply excitation coil, I0For pulse excitation current amplitude;
Step SAP-4, meanwhile, the induced voltage time-domain signal u (t) at the first AD analog-digital converter acquisition testing coil both ends
(unit V), and will be in the induced voltage signal u (t) that collected storage to host;It is acquired and is sampled with the 2nd AD analog-digital converter
Resistance both ends pulse excitation current signal i (t) (unit A), and store into host.
In the present invention, the induced voltage time-domain signal that tested ferromagnetic component is carried out using Pulsed eddy current testing device and
The acquisition of exciting current signal, referred to as Pulsed eddy current testing signal acquisition step, Signal Acquisition Procedure,
SAP。
Step 2: the inversion method (PIP) of tested ferromagnetic component electromagnetic parameter:
After the induced voltage for collecting detection coil both ends according to SAP step, how to be determined by induced voltage and be detected iron
The conductivity and relative permeability of magnetic component are the key that signal processings in the present invention.Author Chen Xingle is equal to 2014, delivers
Paper " Excitation current waveform for eddy on " NDT&E International " volume 66
Current testing on the thickness of ferromagnetic plates ", gives ferromagnetic plate pulse
EDDY CURRENT model, using the ferromagnetic plate Pulsed eddy current testing model, when being passed through pulse excitation electric current i (t) in excitation coil
When, the induced voltage time-domain expression at detection coil both ends are as follows:
π value is 3.14;
E is the bottom of natural logrithm, value 2.72;
σ is the conductivity of tested ferromagnetic component, unit S/m;
D is the wall thickness of tested ferromagnetic component, unit m;
μ0For space permeability, value is 4 π × 10-7H/m;
μrFor the relative permeability for being detected ferromagnetic component;
B is the cutoff boundary of ferromagnetic plate Pulsed eddy current testing model, unit m;
λiIndicate first-order bessel function J1(bλiThe positive root of i-th of)=0;
J0(x) 0 rank Bessel function of the first kind is indicated;
ξkIt is transcendental equationK-th of positive root;
Coefficient
Time constant
T indicates the time;
I (t) is the pulse excitation current signal of sampling resistor acquisition, unit A;
I ' (t) indicates pulse excitation current signal to the derivative of time t;
" * " indicates convolution algorithm
CdFor the coil factor of excitation coil, CpFor the coil factor of detection coil, calculating formula can refer to Chen Xingle and be equal to
2014, " the Excitation current waveform for being published on " NDT&E International " volume 66
An eddy current testing on the thickness of ferromagnetic plates " text.
The present invention is based on the time domain analytic solutions of above-mentioned ferromagnetic plate Pulsed eddy current testing model induced voltage, establish induced electricity
The least square problem between time-domain signal measured value and calculated value is pressed, to the conductivity and magnetic conductivity progress at test point
Inverting, the specific steps are as follows:
Step PIP-1, by coil probe lift-off lo, tested ferromagnetic component wall thickness d and pulse excitation current signal i (t)
It is substituted into formula (1) as known quantity, while by the relative permeability μ of the conductivityσ of tested ferromagnetic component and tested ferromagnetic componentr
It is set as unknown parameter, that is, needs the parameter vector x=(x of inverting1,x2)T=(σ, μr)T;
Step PIP-2, according to SAP step, with the induction at the detection coil both ends that the first AD analog-digital converter collects
Voltage time-domain signal measurement data is (t1,u1),(t2,u2),…,(tm,um), by the reason of its induced voltage calculated with formula (1)
By calculated value u (x, ti) compare, enable induced voltage time-domain signal measured value uiError sum of squares between calculated value is minimum
Carry out inverted parameters x, that is, establish least square problem:
Remember residual function ri(x)=ui-u(x,ti), i=1,2 ..., m, and remember residual functional vector r (x)=(r1(x),r2
(x),...,rm(x))T, R2For two-dimentional real number space, m is the number of measurement, x1And x2Indicate two elements in vector x.
Step PIP-3, using iterative algorithm, solves the optimal solution x of least square problem (2) in host*, iteration
The calculating step of algorithm are as follows:
(1) initial point is given(wherein σ(1)=106~107S/m,Permit
Perhaps error ε > 0 (generally takes ε=10-3), set k=1;
(2) parameter vector for walking kthIt substitutes into calculating formula (1), calculates each of kth step
Time point tiInduced voltage calculated value u (x(k),ti), then with induced voltage time-domain signal measured value uiIt is poor to make, meter
Calculate the residual functional value of kth step
ri(x(k))=ui-u(x(k),ti), i=1,2 ..., m,
And obtain the residual functional vector r of kth step(k);Then induced voltage theoretical curve is further calculated by formula (1)
To the first-order partial derivative of the conductivityσ of tested ferromagnetic component
And induced voltage theoretical curve is to the relative permeability μ of tested ferromagnetic componentrFirst-order partial derivative
Obtain the matrix A of m × 2k=(aij)m×2;J=1,2.
(3) solve system of equation
Acquire direction vector b;
(4) the parameter vector x walked from kth(k)It sets out, along the direction vector b of kth step(k)Make linear search, finds out step-length
αk, so that
I.e. using α as the function f (x of independent variable(k)+αb(k)) when being minimized, the value of independent variable α is αk。
And it enables
x(k+1)=x(k)+αkb(k);
(5) if | | x(k+1)-x(k)| |≤ε then stops calculating, and obtains the optimal solution x of least square problem (2)*=x(k+1);
Otherwise, k:=k+1, return step (2) are set.
Step PIP-4 is solved the optimal solution of least square problem (2) by above-mentioned iterative algorithmAfterwards,
The inversion result σ * for obtaining the conductivity of tested ferromagnetic component, is detected the inversion result of the relative permeability of ferromagnetic component
By the inversion result of the inversion result of the conductivity of tested ferromagnetic component and the relative permeability of tested ferromagnetic component and pulse electricity
Flow excitation amplitude I0It is mapped, is saved in host.
In the present invention, host detects signal analysis and processing to induced voltage, is finally inversed by the conductivity of tested ferromagnetic component
With the process of the relative permeability of tested ferromagnetic component, the referred to as inversion step of parameter, Parameters Inversion
Procedure, PIP.
Step 3: utilizing curve matching, eliminates influence of the magnetic conductivity variation to conductivity testing result:
For ferromagnetic metallic material, the magnetic conductivity of material is difficult directly to measure, and the size of magnetic conductivity is by tested iron
The factors such as magnetic Fine Texture of Material, remanent magnetism, pulse excitation magnetic field strength influence, and therefore, detect every time, are detected the magnetic of ferromagnetic component
Conductance may all change, to influence conductivity testing result.In the present invention, on the one hand by the magnetic of tested ferromagnetic component
Conductance is set as a unknown parameter, and the signal processing method proposed through the invention carrys out the value that inverting determines magnetic conductivity;It is another
Aspect is carried out curve fitting by the inversion result to the conductivity for being detected ferromagnetic component under different amplitude pulse current excitations,
The inversion result of the conductivity of corresponding tested ferromagnetic component when external magnetic excitation field is zero is found as testing result, to cut
Weak magnetic conductance nonlinear change is influenced caused by testing result.Specific implementation step are as follows:
Step 1: obtaining the wall thickness d of tested ferromagnetic component using vernier caliper or audigage measurement, coil is popped one's head in
It is disposed vertically above tested ferromagnetic component, the lift off l of coil probe is gone out using ruler or vernier caliper measuremento;
Step 2: the lift off l of fixed coil probeoIt is constant, by pulse excitation current amplitude I0With 0.2A stepping, from
0.2A is gradually increased 3.0A, and according to SAP step, acquires the corresponding detection coil two of different pulse excitation current amplitudes respectively
The signal for holding induced voltage time-domain signal and sampling resistor both ends pulse excitation electric current, is stored in host;
Step 3: according to PIP step, being finally inversed by the corresponding tested ferromagnetic component of different amplitude pulse exciting currents respectively
The relative permeability of conductivity and tested ferromagnetic component, and draw out tested ferromagnetic component conductivity inversion result with pulse
Exciting current amplitude I0Change curve, as shown in black circle in Fig. 2;
Step 4: utilization index functionTo in Fig. 2 be detected ferromagnetic component conductivity inversion result σ * into
Row curve matching is as illustrated in solid line in figure 2 conductivity matched curve, to obtain factor sigma in exponential function0It is fitted and ties with ε
FruitWithWhereinThe conductivity testing result of as tested ferromagnetic component.