This LS-DYNA simulation shows a simple inductive heating problem. A micro EM timestep is calculated using the circuit’s current period divided by a factor NUMLS. Over a whole period, the full Eddy current problem is solved. An average Joule heating is calculated which is then given to the thermal solver over all the subsequent periods until reaching a time defined by a EM macro timestep. At that point, the matrices are recomputed and the whole procedure is repeated. If the conductor do not move and have a constant conductivity, then the EM macro timestep can be as long as the run. An optional ‘switch.k’ file can be included which turn the EM solver on and off based on a specific node temperature.
Temperature fringes
*BOUNDARY_PRESCRIBED_MOTION_RIGID *CONTROL_SOLUTION *CONTROL_TERMINATION *CONTROL_TIMESTEP *CONTROL_THERMAL_SOLVER *CONTROL_THERMAL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_FUNCTION *DEFINE_CURVE_TITLE *EM_CIRCUIT *EM_CIRCUIT_ROGO *EM_CONTROL *EM_CONTROL_SWITCH *EM_MAT_001 *EM_OUTPUT *EM_SOLVER_FEM *EM_SOLVER_BEM *EM_SOLVER_BEMMAT *END *INITIAL_TEMPERATURE_SET *KEYWORD *MAT_RIGID *MAT_THERMAL_ISOTROPIC *PARAMETER *PART *SECTION_SOLID *TITLE
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (g/mm/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM Inductive heating example
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls emdtinit
2 &numls&em_macrdt
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CIRCUIT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
1 11 025.e3 200.e4
$ ssidCurr ssidVltin ssidVltOt partID
-11 1 2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM MAT and EOS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
1 2 &em_cond1
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
2 4 &em_cond2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM SOLVER $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_SOLVER_BEMMAT
$# matid reltol
1 1e-6
*EM_SOLVER_BEMMAT
$# matid reltol
2 1e-6
*EM_SOLVER_BEM
$# reltol maxit stype precon uselas ncyclbem
1e-6 1000 2 2 1&em_bemmtx
*EM_SOLVER_FEM
$# reltol maxit stype precon uselas ncyclbem
1e-3 1000 1 1 1&em_femmtx
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM OUTPUT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_OUTPUT
$ matS matF solS solF mesh
2 2 2 2 0
*END
