Axisymmetric problem
The R9 version introduced many new advanced capabilities to the axisymmetric solver. For example, it is now possible to simulate pancake or helicoid type of coils by connecting the different turns using EM_CIRCUIT_CONNECT therefore expanding the range of applicable EM forming or bending problems that can be solved using the axisymmetric solver.
https://www.dynaexamples.com/em/eddycurr/em_axi
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Axisymmetric problem
The R9 version introduced many new advanced capabilities to the axisymmetric solver. For example, it is now possible to simulate pancake or helicoid type of coils by connecting the different turns using EM_CIRCUIT_CONNECT therefore expanding the range of applicable EM forming or bending problems that can be solved using the axisymmetric solver.
The R9 version introduced many new advanced capabilities to the axisymmetric solver. For example, it is now possible to simulate pancake or helicoid type of coils by connecting the different turns using EM_CIRCUIT_CONNECT therefore expanding the range of applicable EM forming or bending problems that can be solved using the axisymmetric solver. Run times are tremendously reduced which allows the set up of advanced optimization studies in combination with LS-Opt.
Temperature fringes
*CONTACT_SURFACE_TO_SURFACE *CONTROL_CONTACT *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_FUNCTION *DEFINE_CURVE_TITLE *DEFINE_TABLE *EM_2DAXI *EM_DATABASE_PARTDATA *EM_CIRCUIT *EM_CIRCUIT_CONNECT *EM_CONTROL *EM_CONTROL_TIMESTEP *EM_MAT_001 *EM_OUTPUT *EM_SOLVER_FEM *EM_SOLVER_BEM *EM_SOLVER_BEMMAT *EM_ROTATION_AXIS *END *HOURGLASS *KEYWORD *MAT_ELASTIC *MAT_RIGID *PARAMETER *PART *SECTION_SHELL *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 R9.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (g/mm/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM Electromagnetic axisymmetric 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 3 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
6 4 5
*EM_CIRCUIT
$ circid circtype lcid R L C V0
2 3 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
9 7 8
*EM_CIRCUIT
$ circid circtype lcid R L C V0
3 3 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
12 10 11
*EM_CIRCUIT
$ circid circtype lcid R L C V0
4 3 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
15 13 14
*EM_CIRCUIT
$ circid circtype lcid R L C V0
5 3 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
18 16 17
*EM_CIRCUIT
$ circid circtype lcid R L C V0
6 3 &em_res &em_ind &em_cap &em_v0
$ ssidCurr ssidVltin ssidVltOt partID
21 19 20
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ connId connType circuit1 circuit2 c1 c2
1 1 1 2 1. -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ connId connType circuit1 circuit2 c1 c2
2 1 1 3 1. -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ connId connType circuit1 circuit2 c1 c2
3 1 1 4 1. -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ connId connType circuit1 circuit2 c1 c2
4 1 1 5 1. -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ connId connType circuit1 circuit2 c1 c2
4 1 1 6 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_cond
*EM_MAT_001
$ em_mid mtype sigma eos
2 1
*EM_MAT_001
$ em_mid mtype sigma eos
3 4 &em_cond
*EM_MAT_001
$ em_mid mtype sigma eos
4 4 &em_cond
*EM_MAT_001
$ em_mid mtype sigma eos
5 4 &em_cond
*EM_MAT_001
$ em_mid mtype sigma eos
6 4 &em_cond
*EM_MAT_001
$ em_mid mtype sigma eos
7 4 &em_cond
*EM_MAT_001
$ em_mid mtype sigma eos
8 4 &em_cond
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM 2D AXI $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_ROTATION_AXIS
$ xpoint ypoint zpoint xaxis yaxis zaxis numSec
0. 0. 0. 0. 0. 1. 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
3 3 3 1 1 2 32 1
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
1 6 3 1 4 5 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
4 9 3 1 7 8 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
5 12 3 1 10 11 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
6 15 3 1 13 14 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
7 18 3 1 16 17 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption startSsid endSsid numSec new2D
8 21 3 1 19 20 32
$---+----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-9
*EM_SOLVER_BEMMAT
$# matid reltol
2 1e-9
*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
*EM_DATABASE_PARTDATA
$# outlvl
1
*END