TEAM 13 problem

The TEAM 13 problem is a non linear magnetostatic case where an exciting coil with D.C current is set between two steel channels not aligned with each other and a steel plate inserted between the channels. The applied current is enough to reach steel saturation. This case is an important milestone in the validation of LS-DYNA's FEM/BEM monolithic solver for Eddy currents.

TEAM (Testing Electromagnetic Analysis Methods) represents an open international working group aiming to compare electromagnetic analysis computer codes. TEAM Workshops are meetings of this group. A series of TEAM Workshops was started in 1986 and has been organized in two-year rounds, each comprising a series of ”Regional” workshops and a ”Final” Workshop, as a satellite event of the COMPUMAG Conference. The TEAM problems consist in a list of test-problems, with precisely defined dimensions, constitutive laws of materials, excitations, etc., and each backed by a real laboratory device, on which measurements can be made.

The TEAM 13 problem is a non linear permeability case where an exciting coil with D.C current is set between two steel channels not aligned with each other and a steel plate inserted between the channels. The applied current is enough to reach steel saturation. This case is an important milestone in the validation of LS-DYNA's FEM/BEM monolithic solver for Eddy currents.

The two download links show the two ways the problem can be tackled. Either the results are made time dependent by considering the conductivity of the plate (Eddy current diffusion) or a magnetostatic problem is solved in which case, the first timestep will be the costliest and directly give the saturated solution. For this problem, since the behavior of the magnetic fields is largely driven by the non-linear permeability properties of the steel plates, results should be almost identical between the two configurations.

Reference results can be found in [1].

[1] T. Nakata and K. Fujiwara. Summary of results for benchmark
problem 13 (3-d nonlinear magnetostatic model). COMPEL-The international journal for computation and mathematics in electrical and electronic engineering, 11(3):345–369, 1992.


*EM_CIRCUIT_SOURCE
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_DATABASE_POINTOUT
*EM_MAT_001
*EM_MAT_002
*EM_EOS_PERMEABILITTY
*EM_OUTPUT
*EM_POINT_SET
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_SOLVER_FEM
*EM_SOLVER_FEMBEM_MONOLITHIC

*KEYWORD
*TITLE
TEAM 13
$ Example provided by Iñaki (Ansys LST)
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Recommended LS-DYNA MPP Dev November 2021 (or more recent)
$X Double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (g/mm/s)
$X------------------------------------------------------------------------------
$X
*PARAMETER
R T_end 1.e-4
R dt_plot 5.e-5
$
$--- em
R em_cond 0.0
R em_dt 1.e-5
R curr 1.e6
R rtolbem 1.e-6
R tolbemM 1.e-12
R rtol_ksp 1.e-4
R atol_ksp 1.e-8
R rtol_new 1.e-2
R stol_new 1.e-2
R atol_new 1.e-6
I ls_on 1
$
$
$--- struc
$
R dt_struc 1.e-6
Rrho_struc 7000.e-6
R E_struc 2.e11
R nu_struc 0.3
$ ****************************************************************************
$ mesh
$ ****************************************************************************
*INCLUDE
mesh_team13.k
$ ****************************************************************************
$ mechanics
$ ****************************************************************************
$--------1---------2---------3---------4---------5---------6---------7---------8
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
&T_end
*CONTROL_TIMESTEP
15
*DEFINE_CURVE
15
0.,&dt_struc
1.,&dt_struc
*DATABASE_BINARY_D3PLOT
&dt_plot
*PART
cond
$ pid secid mid eosid hgid grav adpopt tmid
1 1 1
*PART
cond
$ pid secid mid eosid hgid grav adpopt tmid
2 1 2
*SECTION_SOLID
$ sid elform
1 1
*MAT_RIGID
1,&rho_struc,&E_struc,&nu_struc
$# cmo con1 con2
1.000000 7 7
$# lco or a1 a2 a3 v1 v2 v3
0.000 0.000 0.000 0.000 0.000 0.000
*MAT_RIGID
2,&rho_struc,&E_struc,&nu_struc
$# cmo con1 con2
1.000000 7 7
$# lco or a1 a2 a3 v1 v2 v3
0.000 0.000 0.000 0.000 0.000 0.000
$ ****************************************************************************
$ electromagnetism
$ ****************************************************************************
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls emdtinit emdtmax emtinit emtend ncyclFem ncyclBem
1 5000 5000
*EM_CONTROL_TIMESTEP
$--------1---------2---------3---------4---------5---------6---------7---------8
$ TSTYPE DTCONST LCID FACTOR
1 &em_dt 50
*DEFINE_CURVE
50
0.,1.e-5
2.e-3,1.e-5
2.1e-3,2.e-5
5.e-3,2.e-5
5.1e-3,1.e-4
1.,1.e-4
*EM_SOLVER_FEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ relTol maxIter solveType precond useLast
1.e-6 10000 1 1 1
*EM_SOLVER_BEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ relTol maxIter solveType precond uselast
&rtolbem 10000 2 4 1
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
1 &tolbemM
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
2 &tolbemM
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
3 &tolbemM
*EM_CIRCUIT_SOURCE
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R L C V0
1 1 200
$ ssidCurr ssidVltin ssidVltOt partID
1 1 2000
*DEFINE_CURVE
200
0,&curr
1000,&curr
*DEFINE_CURVE
21
0.,0.
10.,0.
*DEFINE_CURVE
22
0.,0.
10.,0.
*DEFINE_CURVE
23
0.,0.01
10.,0.01
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
101 101 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
102 102 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
103 103 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
104 104 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
105 105 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
106 106 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
107 107 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
108 108 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
109 109 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
110 110 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
111 111 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
112 112 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
113 113 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
114 114 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
115 115 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
116 116 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
117 117 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
118 118 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
119 119 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
120 120 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
121 121 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
122 122 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
123 123 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
124 124 1 3
*EM_CIRCUIT_ROGO
$--------1---------2---------3---------4---------5---------6---------7---------8
$ rogoId setID setType currType
125 125 1 3
*EM_MAT_002
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId murel eosMure
$ 2 2 7.505 1.
2 2 &em_cond 1. 10
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
1 1
*EM_EOS_PERMEABILITY
10,1,400
*INCLUDE
BfuncH2.k
*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 0
$ mf2 gmv d3plotFor timeHist
0 0 0
*EM_POINT_SET
1
26,10.,20.,55.0
27,20.,20.,55.0
28,30.,20.,55.0
29,40.,20.,55.0
30,50.,20.,55.0
31,60.,20.,55.0
32,70.,20.,55.0
33,80.,20.,55.0
34,90.,20.,55.0
35,100.,20.,55.0
36,110.,20.,55.0
*EM_DATABASE_POINTOUT
1
1
*EM_SOLVER_FEMBEM_MONOLITHIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$ COUPLING KSP_TYPE KSP_ATOL KSP_RTOL KSP_MAXIT GMRES_RST DEBUG
0 &atol_ksp &rtol_ksp 10000
$--------1---------2---------3---------4---------5---------6---------7---------8
$ NEWT_STOL NEWT_ATOL NEWT_RTOL NEWT_MAXI
0 &stol_new &atol_new &rtol_new 500
$--------1---------2---------3---------4---------5---------6---------7---------8
$ LS_ON LS_FTOL LS_GTOL LS_RTOL LS_SAMPL LS_NUM_SP LS_MAXFUN
0 0.01 0.5 0.01
*END