TEAM 4 problem — Welcome to LS-DYNA Examples

The TEAM 4 problem is a simple Eddy current code benchmarking problem where a solid brick is placed inside a external magnetic field that decays exponentially with time. The main objective of the problem is to calculate the total circulating current in the brick function of time.

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 4 problem consists of a rectangular aluminium brick [1] with a rectangular hole placed in a uniform magnetic field. The magnetic field is perpendicular to the faces with the hole, and decays exponentially with time. The main objective of the problem is to calculate the total circulating current in the brick function of time.

[1] A. Kameari, Results for benchmark calculations of problem 4 (the felix brick), COMPEL, Vol 7. Num 1 (1988).

*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_EXTERNAL_FIELD
*EM_CIRCUIT_ROGO
*EM_DATABASE_POINTOUT
*EM_MAT_001
*EM_OUTPUT
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_SOLVER_FEMBEM

*KEYWORD
*INCLUDE
curve.k
*DEFINE_CURVE
21
0,0
0.1,0
*EM_CIRCUIT_ROGO
1,1,1,1
*EM_EXTERNAL_FIELD
$——–1———2———3———4———5———6———7———8
$ fieldId fieldType FieldDef lcidX lcidY lcidZ
1 1 1 21 21 20
*EM_CONTROL
$——–1———2———3———4———5———6———7———8
$ emsol numls emdtinit emdtmax emtinit emtend ncyclFem ncyclBem
1
*EM_CONTROL_TIMESTEP
1,2.e-4
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$ em_mid mtype sigma eos skinDepth
1 4 2.538e1 0 0
*EM_OUTPUT
$——–1———2———3———4———5———6———7———8
$ matS matF solS solF mesh
1 1 1 1 0
*MAT_RIGID
1 0.0027136.8900E+10 0.330000 0.000 0.000 0
1,7,7

*CONTROL_HOURGLASS
1 0.100000
*CONTROL_TIMESTEP
,,,,,16
*DEFINE_CURVE
16
0,1.e-4
1000,1.e04
*CONTROL_TERMINATION
0.035 0 0.000 0.000 0.000
*SECTION_SOLID
1 1 0
*DATABASE_BINARY_D3PLOT
1.0000E-3 0 0 0 0
0
*PART

1 1 1 0 0 0 0 0
*INCLUDE
mesh_.k
*END

 

Current Density Vectors :

Total Current in the brick :

– EM_example_TEAM4.zip