Induction Hardening
Theory and Physics
Principles of Induction Hardening
Professor, how does induction hardening work?
It induces eddy currents on the surface of a steel component using a coil carrying high-frequency current, rapidly heating only the surface layer via the skin effect, followed by rapid quenching to harden it through martensitic transformation.
Heating depth $\approx \delta$:
For steel, $\mu$ changes significantly with temperature. At the Curie temperature (approx. 770°C), $\mu_r \to 1$ and $\delta$ increases sharply.
So the heating depth suddenly increases when it exceeds the Curie temperature.
Guidelines for frequency and heating depth:
| Frequency | Hardened Layer Depth | Application |
|---|---|---|
| 1–10 kHz | 3–10 mm | Large gears, shafts |
| 10–100 kHz | 1–3 mm | Medium-sized parts |
| 100–500 kHz | 0.3–1 mm | Small parts, thin layers |
Summary
- Surface heating via skin effect — Control hardening depth with frequency
- Curie temperature — $\delta$ changes due to abrupt change in $\mu_r$
- Eddy current loss → Joule heating — A coupled electromagnetic-thermal problem
Principles of Induction Hardening—The Skin Effect "Rapidly Heats Only the Skin of Steel"
Induction hardening is a heat treatment method that rapidly heats only the surface layer of a steel component using eddy currents to austenitize it, then cools it to generate martensite (a high-hardness phase). Since the skin depth δ is inversely proportional to the square root of frequency, the heating depth can be controlled by selecting the frequency (high frequency → thin hardened layer, low frequency → deep hardened layer). This design freedom of "determining depth by frequency" makes induction hardening a technology suitable for precision heat treatment of engine components and gears.
Physical Meaning of Each Term
- Electric field term $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$: Faraday's law of electromagnetic induction. A time-varying magnetic flux density generates an electromotive force. [Everyday Example] A bicycle dynamo (generator) generates voltage in a nearby coil by rotating a magnet—a direct application of this law that a changing magnetic field induces an electric field. Induction cooking (IH) heaters also use the same principle, where high-frequency magnetic field changes induce eddy currents in the pot bottom, heating it via Joule heat.
- Magnetic field term $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$: Ampère-Maxwell's law. Current and displacement current generate a magnetic field. [Everyday Example] When current flows through a wire, a magnetic field is created around it—this is Ampère's law. Electromagnets operate on this principle, passing current through a coil to create a strong magnetic field. Smartphone speakers also apply this law: current → magnetic field → force on the diaphragm. At high frequencies (e.g., GHz-band antennas), the displacement current $\partial D/\partial t$ cannot be ignored, describing electromagnetic wave radiation.
- Gauss's law $\nabla \cdot \mathbf{D} = \rho_v$: Indicates that electric charge is the divergence source of electric flux. [Everyday Example] Rubbing hair with a plastic sheet creates static electricity, making hair stand up—the charged sheet (electric charge) radiates electric field lines outward, exerting force on the light hair. Capacitor design calculates the electric field distribution between electrodes using this law. ESD (electrostatic discharge) countermeasures are also based on electric field analysis following Gauss's law.
- Magnetic flux conservation $\nabla \cdot \mathbf{B} = 0$: Expresses that magnetic monopoles do not exist. [Everyday Example] Cutting a bar magnet in half does not create a magnet with only an N pole or only an S pole—N and S poles always exist as a pair. This means magnetic field lines form "closed loops with no start or end points." In numerical analysis, the formulation using vector potential $\mathbf{B} = \nabla \times \mathbf{A}$ is used to satisfy this condition, automatically guaranteeing magnetic flux conservation.
Assumptions and Applicability Limits
- Linear material assumption: Permeability and permittivity are independent of magnetic/electric field strength (nonlinear B-H curve needed in saturation region)
- Quasi-static approximation (low frequency): Displacement current term can be ignored ($\omega \varepsilon \ll \sigma$). Common in eddy current analysis.
- 2D assumption (cross-section analysis): Effective when current direction is uniform and edge effects can be ignored.
- Isotropic assumption: For anisotropic materials (e.g., rolling direction of silicon steel sheets), direction-specific property definitions are needed.
- Non-applicable cases: Additional constitutive laws are needed for plasma (ionized gas), superconductors, nonlinear optical materials.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Magnetic flux density $B$ | T (Tesla) | 1T = 1 Wb/m². Permanent magnets: 0.2–1.4T |
| Magnetic field strength $H$ | A/m | Horizontal axis of B-H curve. Conversion from CGS unit Oe (Oersted): 1 Oe = 79.577 A/m |
| Current density $J$ | A/m² | Calculated from conductor cross-sectional area and total current. Note non-uniform distribution due to skin effect. |
| Permeability $\mu$ | H/m | $\mu = \mu_0 \mu_r$. In vacuum $\mu_0 = 4\pi \times 10^{-7}$ H/m |
| Electrical conductivity $\sigma$ | S/m | Copper: approx. 5.96×10⁷ S/m. Decreases with temperature rise. |
Numerical Methods and Implementation
Electromagnetic-Thermal Coupled FEM
How do you set up a simulation for induction hardening?
Couple two governing equations.
Electromagnetic field: $\nabla \times (\nu \nabla \times \mathbf{A}) + \sigma\partial\mathbf{A}/\partial t = \mathbf{J}_0$
Heat conduction: $\rho c_p \partial T/\partial t = \nabla \cdot (k \nabla T) + Q_{eddy}$
$Q_{eddy} = |\mathbf{J}|^2/\sigma$ is the eddy current heat generation. $\mu$, $\sigma$, $k$, $c_p$ are all temperature-dependent.
Is strong coupling necessary?
Due to the strong temperature dependence of $\mu$, weak coupling (alternating calculation) with updates every few steps yields sufficient accuracy. However, time steps need to be refined near the Curie temperature.
Summary
- Coupled electromagnetic+thermal FEM — Pass $Q_{eddy}$ as a heat source
- Temperature dependence of material properties — $\mu(T)$, $\sigma(T)$ are important
- Weak coupling is practical — Manageable with updates every few steps
FEM for Induction Hardening—Triple Coupling of Electromagnetic, Thermal, and Microstructural Transformation is Required
Numerical analysis of induction hardening requires coupling three phenomena: electromagnetic field analysis (calculation of eddy current heat generation), heat conduction analysis (calculation of temperature distribution), and metallic microstructural transformation (austenitization and martensitic transformation). Particularly, ignoring the latent heat of microstructural transformation (martensitic transformation heat) can cause calculated cooling rates to deviate from measurements by 10–20%, reducing prediction accuracy for hardened layer depth. Since each property (permeability, thermal conductivity, specific heat) depends on temperature and microstructure, nonlinear coupled analysis is essential.
Edge elements (Nedelec elements)
Elements specialized for electromagnetic field analysis. Automatically guarantee continuity of tangential components and eliminate spurious modes. Standard for 3D high-frequency analysis.
Nodal elements
Used for scalar potential formulations. Effective for scalar potential methods in magnetostatics and electrostatic field analysis.
FEM vs BEM (Boundary Element Method)
FEM: Handles nonlinear materials and inhomogeneous media. BEM: Naturally handles infinite domains (open region problems). Hybrid FEM-BEM is also effective.
Nonlinear Convergence (Magnetic saturation)
Nonlinearity of B-H curve handled by Newton-Raphson method. Residual criterion: $||R||/||R_0|| < 10^{-4}$ is typical.
Frequency Domain Analysis
Reduced to a steady-state problem via time-harmonic assumption. Requires complex number operations, but broadband characteristics are obtained via time-domain analysis.
Time Domain Time Step
Time step less than 1/20 of the highest frequency component is required. Implicit time integration allows larger steps but accuracy must be considered.
Choosing Between Frequency Domain and Time Domain
Related Topics
なった
詳しく
報告