D.C Electric motor
This example shows some of the capabilities of the LS-DYNA EM solver (Post R13). A standard voltage of 12 V is applied between the two brushes of the motor. Permanent magnets located on each side of the armature cause the branches to start moving. The EM contact algorithm allows the current path to "switch" between the different branches and maintain the movement. The use of LS-DYNA FEM/BEM method to handle the electromagnetic fields as well as its integrated and automatic coupling with the structure makes it a powerful tool for such simulations.
*EM_SOLVER_BEMMAT
*EM_CIRCUIT
*EM_CONTACT
*EM_CONTROL
*EM_CONTROL_CONTACT
*EM_CONTROL_TIMESTEP
*EM_MAT_001
*EM_MAT_002
*EM_OUTPUT
*EM_PERMANENT_MAGNET
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_FEMBEM_MONOLITHIC
*KEYWORD
$ ****************************************************************************
$ electromagnetism
$ ****************************************************************************
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls dt
1 &em_femmtx&em_bemmtx
*EM_CONTROL_TIMESTEP
1,&emdt
$*EM_CONTROL_COUPLING
$,2
*EM_CIRCUIT
1,2,22
2,1,2
*DEFINE_CURVE
22
0,0.
5.0e-4,0.
1.e-3,&volt
1000,&volt
*EM_CONTACT
1,1,30,41,&eps1,&eps2,,&d0
*EM_CONTACT
2,1,30,42,&eps1,&eps2,,&d0
*EM_CONTROL_CONTACT
1,1,2,1
*EM_SOLVER_FEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ reltol maxiter stype precon uselast ncyclfem msglvlFem
1.e-4 10000 1 1 1 0
*EM_SOLVER_BEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ reltol maxiter solvetype precon uselast ncyclbem
&btol 10000 2 2 1
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
1 &ptol
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
2 &qtol
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
3 &wtol
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$ ID Part ID - North South M=0.38/mu0
24 24 1 2 &mag
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$ ID Part ID - North South M=0.38/mu0
25 25 3 4 &mag
*DEFINE_CURVE
100,
0,9.e5
5e-3,9.e5
1.,9.e5
*EM_MAT_001
11,1
*EM_MAT_002
$ Magnet material 1
$ mid mtype sigma eosid murel eosmu
24 4 &sigma2 0 &murel2
$ 25 4 &sigma1 0 &murel1
*EM_MAT_002
$ Magnet material 2
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid murel eosmu
25 4 &sigma2 0 &murel2
*EM_MAT_001
$ Armature
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid
31 4 &sigma3
*EM_MAT_001
$ Armature
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid
32 4 &sigma3
*EM_MAT_001
$ Armature
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid
33 4 &sigma3
*EM_MAT_001
$ Armature
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid
34 4 &sigma3
*EM_MAT_001
$ coil material
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid
41 4 &sigma4
*EM_MAT_001
$ coil material
$--------1---------2---------3---------4---------5---------6---------7---------8
$ mid mtype sigma eosid
42 4 &sigma4
*EM_OUTPUT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matS matF solS solF mesh memory timing d3plotAsc
2 2 2 2 0 0 0
$ mf2 gmv d3plotFor timeHist
$ 0 1 0
*EM_SOLVER_FEMBEM_MONOLITHIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$ COUPLING KSP_TYPE KSP_ATOL KSP_RTOL KSP_MAXIT GMRES_RST FEM_PC DEBUG
0 &monotola &monotolr 100000
$--------1---------2---------3---------4---------5---------6---------7---------8
$ FORCE NEWT_STOL NEWT_ATOL NEWT_RTOL NEWT_MAXI
$ 1.e-9 1.e-9 1.e-5
$--------1---------2---------3---------4---------5---------6---------7---------8
$ LS_ON LS_FTOL LS_GTOL LS_RTOL LS_SAMPL LS_NUM_SP LS_MAXFUN
$ 1 1.e-9 0.01 1.e-5 30000
$--------1---------2---------3---------4---------5---------6---------7---------8
$$USE_ENERG STOP_ERR WIDTH_ERR XTRAPF
*END
