Modeling High Frequency Electromagnetic Skin Depth Variations in Strain Transformed Multi Phase Alloys
Strain transformation collapses high-frequency skin depth by increasing permeability, requiring multi-frequency impedance modeling to decouple phase from defects.

Matrix
Plastic deformation in metastable austenitic stainless steels and transformation-induced plasticity alloys alters the fundamental parameters of electromagnetic wave propagation. When face-centered cubic austenite undergoes mechanical shear, deformation twins and stacking faults accumulate until body-centered cubic or body-centered tetragonal martensite nucleates along shear band intersections. This crystallographic restructuring generates drastic physical divergences in non-destructive testing and inductive sensor measurements.
Austenite exhibits a relative magnetic permeability of approximately 1.003 at room temperature. Martensite formed under plastic strain behaves ferromagnetically, driving local relative permeability values between 15 and 95 under weak interrogating fields.
The electrical conductivity also shifts as lattice dislocations proliferate and phase boundaries multiply. Cold-drawn AISI 304 sheet exhibits an electrical conductivity of 1.38 megasiemens per meter in the annealed state. Following forty percent plastic deformation, the accumulation of point defects, grain elongation, and alpha-prime phase boundaries reduces electrical conductivity to 1.22 megasiemens per meter.
The electromagnetic penetration depth depends mathematically on both properties according to the classical relationship where skin depth equals the inverse square root of pi, excitation frequency, magnetic permeability, and electrical conductivity.
Interrogation frequencies between 100 kilohertz and 20 megahertz compress this penetration into the deformed surface layer. At 1 megahertz, the calculated skin depth for annealed austenite sits at 428 micrometers. Once strain-induced martensite reaches a volume fraction of thirty-five percent, the relative permeability rises to 45 at an excitation magnetic field of 10 microteslas, dropping the local skin depth to 68 micrometers.
The interrogating field encounters an impedance boundary that confines eddy currents entirely within the altered surface profile.
A thirty percent cold-reduction in AISI 304 drives skin depth down from 428 micrometers to 68 micrometers at 1 megahertz.
Phase transformation proceeds non-uniformly through sheet thickness. Surface friction against tooling, bending moments during stamping, and localized shear bands generate steep martensite volume gradients from the outer surface to the core. High-frequency eddy current coils capture an integrated impedance response across this stratified profile rather than a bulk homogeneous material parameter.
- Permeability compression occurs when external mechanical stress biases magnetic domains, shifting the apparent skin depth without phase fraction changes.
- Liftoff suppression decouples probe standoff variations from subsurface magnetic permeability gradients through tuned bridge phase angles.
- Conductivity shadowing masks deep martensitic transformation zones when high surface dislocation density concentrates the eddy current path.
Ignoring the dual shift in permeability and conductivity causes automated eddy current sorting equipment to report false crack indications or misclassify work-hardened pressings as scrap metal.

Diffusion
Planar electromagnetic propagation through an inhomogeneous conductor requires solving the one-dimensional Helmholtz diffusion equation with position-dependent material parameters. In a homogenous half-space, magnetic field intensity decays exponentially with depth. When severe plastic deformation establishes a depth-dependent profile of relative permeability and electrical conductivity, the governing wave equation incorporates spatial derivatives of magnetic permeability alongside the classical diffusion term.
Let z denote depth from the sample surface into the material interior. The differential governing relationship follows the form where the second spatial derivative of the magnetic field intensity minus the logarithmic spatial derivative of magnetic permeability times the first spatial derivative of the magnetic field intensity equals the product of imaginary unit, angular excitation frequency, local magnetic permeability, local electrical conductivity, and the magnetic field intensity. The explicit appearance of the permeability gradient acts as an internal reflection mechanism, scattering electromagnetic energy back toward the coil face.

What Governs the High Frequency Permeability Roll-Off?
Eddy current dampening restricts magnetic domain wall displacement as excitation frequencies pass 500 kilohertz. The effective reversible permeability of strain-induced alpha-prime martensite decays monotonically across high excitation spectra. At quasi-static frequencies below 1 kilohertz, magnetic domain walls oscillate freely under low-amplitude alternating fields, yielding peak relative permeability.
Above 2 megahertz, domain wall motion freezes, leaving spin rotation as the dominant magnetization mechanism, which caps the relative permeability of the transformed phase at values between 4 and 12.
| True Plastic Strain | Martensite Fraction (%) | Conductivity (MS/m) | Relative Permeability at 100 kHz | Relative Permeability at 5 MHz | Calculated Skin Depth at 5 MHz (um) |
|---|---|---|---|---|---|
| 0.00 | 0.0 | 1.38 | 1.003 | 1.002 | 191.4 |
| 0.10 | 4.5 | 1.34 | 3.800 | 2.100 | 134.1 |
| 0.20 | 14.2 | 1.29 | 18.500 | 6.400 | 78.1 |
| 0.30 | 28.7 | 1.25 | 44.000 | 10.200 | 63.1 |
| 0.40 | 46.1 | 1.21 | 72.000 | 14.800 | 53.2 |
Layered discretization models partition the surface-affected zone into discrete planar slabs of uniform physical properties. For a stamped sheet showing a 200-micrometer work-hardened zone, dividing the volume into twenty individual 10-micrometer layers allows transmission line transfer matrix methods to compute input impedance at the air-metal interface. Each layer maintains a characteristic propagation constant and intrinsic electromagnetic wave impedance derived from its local strain-induced properties.
High excitation frequencies isolate the near-surface strata. At 10 megahertz, electromagnetic energy remains confined within the initial 40 micrometers, providing sensitivity to surface shear bands while eliminating interference from substrate thickness tolerances. Conversely, multi-frequency sweeps traversing 100 kilohertz to 10 megahertz provide the phase trajectory signatures required to invert depth profiles.
The boundary condition between adjacent computational strata mandates continuous tangential magnetic and electric fields, capturing the phase shift induced by sharp martensite fronts.
The unresolved physical question centers on whether sub-surface residual stress gradients distort high-frequency domain wall relaxation independently of the actual volume fraction of transformed martensite.

Coil
Inductive sensing elements translate electromagnetic field variations into measurable impedance vectors. An air-cored probe placed above a multi-phase metallic surface experiences reflected impedance changes governed by Dodd and Deeds analytical boundary formulations. When the underlying alloy undergoes localized strain-induced phase transformation, coil inductance increases due to ferromagnetic coupling, while coil resistance rises as eddy currents concentrate inside the compressed skin depth layer.
High-frequency probe architecture demands tight control of parasitic winding capacitance. Self-resonance frequencies must sit at least one decade above the operating measurement band. For a probe operating between 1 megahertz and 15 megahertz, miniature planar spiral coils on polyimide substrates or single-layer ferrite-shielded cylindrical windings deliver clean impedance planes.
Ferrite shielding suppresses magnetic flux leakage, confining interrogation fields to a spot diameter comparable to coil outer radius.

Will Multi Frequency Excitation Resolve Case Hardening?
Simultaneous multi-tone excitation separates the non-linear surface phase fraction from gradual mechanical hardening. Driving an interrogation transducer with superimposed carrier frequencies at 250 kilohertz, 2 megahertz, and 8 megahertz generates three discrete impedance loci. Lower frequencies penetrate through the deformed casing into the paramagnetic core, measuring bulk sheet dimensions and baseline conductivity.
Higher frequencies reflect near-surface martensitic transformation and micro-cracking.
| Transducer Format | Outer Diameter (mm) | Inductance in Air (uH) | Resonance Peak (MHz) | Usable Bandwidth | Spot Resolution (mm) |
|---|---|---|---|---|---|
| Planar Spiral PCB | 1.8 | 0.42 | 84 | 500 kHz to 20 MHz | 0.8 |
| Shielded Ferrite Pot | 3.2 | 1.85 | 38 | 100 kHz to 10 MHz | 1.5 |
| Differential Split-D | 2.5 | 0.95 | 52 | 250 kHz to 15 MHz | 1.1 |
| Air-Core Solenoid | 5.0 | 4.20 | 19 | 50 kHz to 5 MHz | 3.0 |
Signal conditioning hardware balances the probe within an alternating-current bridge circuit feeding an analog demodulator. In-phase and quadrature voltage components correspond directly to real and imaginary impedance deviations. Bridge drift stemming from ambient temperature variations requires continuous compensation, as copper coil resistance changes by 0.39 percent per kelvin.
- Ferrite core saturation occurs under strong permanent magnet bias, distorting impedance linearity during multi-phase characterization.
- Cable capacitance loading suppresses high-frequency bridge sensitivity when line lengths between coil and demodulator exceed two meters.
- Contact liftoff variations introduce large reactive signals that mask small changes in transformed volume fraction.
Thermal resistance shifts inside uncompensated coil windings dwarf subtle martensitic phase signals during continuous inline inspection.
Instrumentation suppliers routinely assert that phase-angle rotation alone eliminates liftoff interference across any arbitrary metallic alloy without recalibration.

Inversion
Extracting continuous phase-depth profiles from discrete impedance data represents an ill-posed inverse problem. As frequency increases, the spatial sensitivity kernel of the interrogation coil concentrates progressively toward the outer boundary. Inverting swept-frequency impedance data reconstructs the spatial variation of magnetic permeability and electrical conductivity through regularized optimization routines.
The forward model computes predicted coil impedance using discretized transfer matrices across twenty or more physical layers. Discrepancies between measured impedance vectors and simulated model points form an objective residual vector. Because experimental noise and instrumentation drift destabilize direct matrix inversion, Tikhonov regularization stabilizes the numerical solution by penalizing unphysical spatial oscillations in reconstructed permeability profiles.
- Baseline calibration on fully annealed standards establishes baseline conductivity and confirms bridge balance across all evaluation frequencies.
- Liftoff trajectory establishment records probe response across calibrated non-conductive shims to establish the phase angle of physical separation.
- Swept frequency measurement across deformed parts acquires real and imaginary impedance values at a minimum of eight logarithmic frequency points between 100 kilohertz and 10 megahertz.
- Regularized inversion execution matches measured impedance loci against a multi-layer forward model to extract the depth-resolved martensite volume fraction.
Magnetoelastic coupling introduces substantial cross-sensitivity during profile reconstruction. Tensile residual stresses expand magnetic domain spacing, enhancing local magnetic permeability parallel to the stress axis, while compressive stresses compress domains and lower apparent permeability. This mechanical effect alters the apparent high-frequency skin depth independently of the crystallographic phase transformation.
A component standard specifies calibration on identical alloy grades processed through matching thermal cycles before stamping.
Separating stress-induced permeability variations from true martensitic phase conversion requires calibrating the inversion kernel against reference coupons subjected to verified tensile strains under vibrating sample magnetometer control. Without this mechanical calibration baseline, mathematical inversion routines interpret residual elastic tension as localized material transformation.
Thick transformed surface zones damp high-frequency eddy currents, leaving deeper microstructural features invisible to electromagnetic interrogation.

Scrap
Automated manufacturing operations processing metastable austenitic stainless steels and high-strength dual-phase sheets face severe economic loss from mismanaged eddy current rejection thresholds. In automotive press shops stamping deep-drawn battery enclosures or structural crossmembers, localized strain transformation naturally occurs along draw radiuses. Inappropriate high-frequency sensor calibration flags these work-hardened zones as cracks, wall-thinning defects, or unapproved material substitutions.
| Inspection Principle | Measurement Speed (parts/min) | Sensitivity to Cracks (um) | Phase Quantification Capability | Liftoff Tolerance (mm) | Relative Capital Cost |
|---|---|---|---|---|---|
| Single-Frequency Eddy Current | 120 | 50 | Poor | 0.2 | 1.0 |
| Multi-Frequency Impedance | 45 | 25 | High | 0.5 | 2.8 |
| Magnetic Flux Leakage | 90 | 100 | Moderate | 1.5 | 1.4 |
| Resonant Ultrasound | 15 | 200 | None | Direct contact | 3.5 |
Single-frequency eddy current instruments configured for standard non-ferrous metals fail immediately on deformed multi-phase alloys. As a stamped part passes under an inspection probe at 500 kilohertz, the abrupt change in skin depth caused by alpha-prime martensite swings the reflected bridge signal across the crack-detection threshold. Production facilities discard structurally sound stampings, driving scrap rates past twelve percent on newly commissioned tooling.
Multi-frequency instrumentation resolves this operational bottleneck by tracking the phase angle trajectory. True physical cracks disrupt current circulation paths, causing an instantaneous resistance increase paired with negligible inductance variation. Strain-induced transformation shifts both inductive and resistive components along a characteristic arc determined by the non-uniform skin depth collapse.
Signal classification algorithms embedded within field-programmable gate arrays isolate crack indications from phase changes within four milliseconds.
Procurement contracts specifying incoming alloy verification require rigorous electromagnetic acceptance standards. Delivery terms according to ISO 21968 and ASTM E1004 mandate that eddy current sorting gates define permissible relative permeability boundaries alongside minimum electrical conductivity values, preventing the wrongful rejection of standard production lots exhibiting normal cold-work transformation.


