2.1. Basic Principles of Electromagnetic Ultrasonic Testing
The electromagnetic ultrasonic detection techniques dealt with in this paper all take electrical conductors as the measured objects, and their basic principles are based on electromagnetic induction phenomena, and the fundamental laws of macroscopic electromagnetic phenomena can be expressed very simply by Maxwell's equations [
19], namely:
where
is the electric field strength;
is the magnetic field strength;
is the electric flux density;
is the magnetic flux density;
is the current density;
is the charge density.
In general, the electromagnetic properties of a medium under the action of an electromagnetic field can also be expressed by the following relation:
where
is the conductivity of the medium;
is the dielectric constant of the medium;
is the magnetic permeability of the medium.
Formulas (1) - (7) form the physical basis of the electromagnetic detection technique in this paper, and the solution of the electromagnetic detection problem can be found by combining the above formulas and by setting the boundary conditions of the specific model and applying the force.
2.2. Antenna Wheel-Rail Electromagnetic Ultrasonic Inspection Mechanism
Electromagnetic ultrasonic detection mechanism is: in close proximity to the measured metal surface of the coil with high-frequency alternating current, the measured metal surface skin depth will produce a change in the electromagnetic field, and then in the surface will produce a same frequency of the induced eddy currents, the induced eddy currents in the applied bias magnetic field, resulting in the same frequency of the Lorentz force, which makes the specimen inside the point of the mass of the periodic vibration of the ultrasonic excitation, so as to detect the wheel-rail [
20,
21].
Electromagnetic Acoustic Transducer (EMAT) is a non-contact transducer that can be used to detect damage on the surface of an antenna wheel-rail by means of electromagnetic coupling. EMAT is a non-contact transducer, which excites and receives ultrasonic waves in a conductor by electromagnetic coupling, and can easily excite various types of ultrasonic waves. Among them, the body wave can well detect the internal defects of the specimen, but there is a blind spot for surface and near-surface defects [
22]; while the energy of the surface wave is mainly concentrated in the skinning depth of the specimen surface when it is propagated, which is suitable for the detection of surface opening defects of thick plates [
23,
24]. The application in this paper is the detection of surface and near-surface defects by surface waves. EMAT is generally composed of a combination of a high-frequency coil, a biased magnetic field and a conductor under test. The principle is shown in
Figure 2.
When the coil is passed through an alternating excitation current, a varying magnetic field is generated, and the alternating magnetic field induces eddy currents on the surface of the antenna wheel-rail, namely:
The Lorentz force is formed by the eddy currents in the wheel track under the action of the biased magnetic field
, namely:
where
is the induced eddy current;
is the alternating magnetic field strength;
is the Lorentz force;
is the bias magnetic field strength.
The Lorentz force causes the mass inside the wheel-rail to vibrate, which generates ultrasonic waves. When there are no defects, the ultrasound waves propagate along the surface, and the energy attenuation on the propagation path is very small during acoustic wave detection, it allows the ultrasound guided waves to propagate over long distances on the surface of the antenna wheel-rail or inside the antenna wheel-rail. When there are defects on the surface or near surface of the wheel-rail, the ultrasonic waves encountering the defects will undergo reflection phenomena and transmission phenomena, which involves the fluctuation problem of elastic dynamics, the solution of it is based on the fundamental solution of elastic dynamics. The research models in this paper are all based on 2-D model of the antenna wheel-rail, and the assumption of 2-D plane strain problem is made here. It is assumed that the defects in the detected 2-D model are actually defective grooves in the direction of the perpendicular plate cross-section.
The following boundary integral equation for elastic dynamics can be introduced using the weighted residual method:
where
and
(j=1,2;k=1,2) are the elastic dynamics basis solutions, respectively represents the displacement and the surface force;
is the two-dimensional region boundary;
,
are the displacements in the
and
directions; and
,
are the surface forces of
in the
and
directions. According to the division of the boundary cells, the boundary equations with continuous integration over the whole boundary can be discretized into the following form:
where
is the total number of boundary cells.
Constant element is used in the boundary element program discussed in this paper. This element is geometrically represented as a straight line segment, and the physical quantity is constant on each element. Although the calculation accuracy of constant elements is not as high as that of more complex linear and higher-order elements, according to the research of Cho[
25,
26], the constant elements can also get better results in elastic wave problems.
If constant cells are used, the displacements and surface forces in the integral expressions contained in the boundary integral equations in discrete form can be referred to outside the integral sign, so that what needs to be solved is actually the integral of the underlying solution. The boundary integral in discrete form can be written as:
The unit number
is iterated within the unit number range
(
), and finally
equations are obtained, which are sorted as follows:
where:
The solution of the elements of the coefficient matricesand needs to be discussed in two cases: 1) when, and contain non-singular integrals, which can be solved directly by using the Gaussian product formula machine; 2) when, the integrals in and have singularities at , which have to be dealt with in a special way.
After obtaining all the coefficients of the complex coefficient matrices
and
, it is necessary to introduce the boundary conditions to obtain the system of linear equations for the final required solution. The hybrid boundary element model of a rectangular defective slot is shown in
Figure 3. Ultrasonic waves propagate within the elastic flat plate interface, the upper and lower surfaces of the plate and the defects are free boundaries with a surface force of 0; the so-called virtual boundary is established a certain distance away from the left and right sides of the defects, and the virtual boundary and the free surfaces together form a closed model area; the appropriate frequency thickness product is selected according to the design, so that a single-mode RAM wave is incident from the left virtual boundary into the model area and interacts with the defects, a part of the RAM wave reflects back to the left. The single mode Lamb wave is incident from the left virtual boundary and interacts with the defects, part of the Lamb wave is reflected back to the left virtual boundary, and the other part of the Lamb wave is transmitted to the right virtual boundary.
When a symmetric or antisymmetric mode Lamb wave of a certain frequency thickness and order is input to the left virtual boundary of the model as an incident mode, the reflection and transmission modes obtained after the interaction of the incident mode and the defect may contain any mode that can exist at that frequency thickness, as determined by the dispersion curve. For a linear material, the total displacement of the left boundary should be the sum of the incident mode displacements and the displacements of all possible reflection modes, expressed as follows:
Similarly the right boundary displacement should be the sum of the displacements of all possible transmission modes with the following expression:
where
is the total displacement;
is the reflection mode displacement;
is the transmission mode displacement.