This LS-DYNA simulation shows a simple Tube expansion example using the EM solver. A current is imposed in the coil which produces induced currents in the workpiece. This in turn generates Lorentz forces which cause the tube to expand. This problem is of medium size so running with 4 or more CPUs is recommended.
EM current density 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 *EM_SOLVER_BEM *EM_SOLVER_BEMMAT *END *INITIAL_TEMPERATURE_SET *KEYWORD *MAT_024 *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 Electromagnetic forming problem
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$ emsol
1
*EM_CONTROL_TIMESTEP
$# tstype dtcons lcid
1 &em_dt
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CIRCUIT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CIRCUIT
$ circid circtype lcid R L C V0
1 1 10 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
2 1 3
*EM_CIRCUIT_ROGO
$ rogoId setId setType currType
1 4 1 1
*EM_CIRCUIT_ROGO
$ rogoId setId setType currType
2 2 1 1
$---+----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 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
