Resistive heating example
This LS-DYNA simulation shows a simple resistive heating problem. The EM resistive heating solver considers the rise time of the current to be slow or non existent (DC current). As such, no inductive and diffusive effects appear which allows the solver to remove the BEM system. The EM solve is consequently very fast and a very high EM timestep can be used.
https://www.dynaexamples.com/em/em_resistive
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Resistive heating example
This LS-DYNA simulation shows a simple resistive heating problem. The EM resistive heating solver considers the rise time of the current to be slow or non existent (DC current). As such, no inductive and diffusive effects appear which allows the solver to remove the BEM system. The EM solve is consequently very fast and a very high EM timestep can be used.
This LS-DYNA simulation shows a simple resistive heating problem. The EM resistive heating solver considers the rise time of the current to be slow or non existent (DC current). As such, no inductive and diffusive effects appear which allows the solver to remove the BEM system. The EM solve is consequently very fast and a very high EM timestep can be used. In the current example, no heat convection or radiation is taken into account for simplification purposes.
Temperature fringes
*CONTROL_SOLUTION *CONTROL_TERMINATION *CONTROL_TIMESTEP *CONTROL_THERMAL_SOLVER *CONTROL_THERMAL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *EM_CIRCUIT *EM_CIRCUIT_ROGO *EM_CONTROL *EM_CONTROL_TIMESTEP *EM_MAT_001 *EM_OUTPUT *EM_SOLVER_FEM *END *INITIAL_TEMPERATURE_SET *KEYWORD *MAT_ELASTIC *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: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM Resistive heating example
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$# emsol
3
*EM_CONTROL_TIMESTEP
$# tstype dtconst
1 &em_dt
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM MAT and EOS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$ em_mid mtype sigma eos
1 2 &em_cond1
*EM_MAT_001
$ em_mid mtype sigma eos
2 2 &em_cond2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CIRCUIT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CIRCUIT
$# circid circtyp lcid r/f l/a c/t0 v0
1 1 10
$# sidcurr sidvin sidvout partid
3 1 2
*EM_CIRCUIT_ROGO
$# rogid setid settype curtyp
1 3 1 1
*DEFINE_CURVE_TITLE
current curve
10
0,0
0.1,5000
0.2,2000
29,2000
30,800
120,0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM SOLVER $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*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