TEAM 3 problem
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 3 problem is a classic validation test case often studied ([3], [1]) and consists of a conducting ladder, having two holes with a coil above carrying a sinusoidal current [2].
The current in the coil is considered uniform (no eddy currents) while the ladder’s induced current diffuses through its thickness (eddy currents). Several positions of the coil were originally studied and for this test case, we are going to focus on the case where the coil is located above one of the holes.
The main objective of this test case is to study the behavior of the magnetic field along a line A-B that goes from (x=0, y=-55, z=0.5) mm to (x=0, y=55, z=0.5) mm (i.e along the symmetrical axis of the problem in the (x,y) plane, and between the coil and the ladder in the z direction).
[1] Z. A. Razek, New technique for solving three-dimensional multiply connected eddycurrent problems, IEE Proceedings, 137 (1990).
[2] L. R. Turner, Proceedings of the vancouver team workshop, (1988).
[3] D. Zheng, Three dimensional eddy current analysis by the boundary element method,
IEE Transaction of Magnetics, (1997).
*EM_BEM_MATRICES
*EM_BEM_MATRICES
*EM_CIRCUIT_SOURCE
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_DATABASE_POINTOUT
*EM_MAT_001
*EM_MAT_001
*EM_OUTPUT
*EM_POINT_SET
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_FEMBEM
*EM_SOLVER_FEMBEM_MONOLITHIC
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LST-Ansys)
$
$
$X
$X------------------------------------------------------------------------------
$# UNITS: g/mm/s
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
TEAM 3 Problem
$ ****************************************************************************
$ mesh
$ ****************************************************************************
*INCLUDE
mesh_team3_coarse.k
$ ****************************************************************************
$ mechanics
$ ****************************************************************************
$--------1---------2---------3---------4---------5---------6---------7---------8
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
2e-1
*CONTROL_TIMESTEP
3
*DATABASE_BINARY_D3PLOT
5.e-4
*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,7000e-6,200e9,0.3
$# 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,7000e-6,200e9,0.3
$# 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_ELASTIC
5 2687.1e-6 69.e+09 .3
*MAT_ELASTIC
3 8928.6e-6 200.e+09 .3
*CONTROL_STRUCTURED
*HOURGLASS
3
$ ****************************************************************************
$ 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,1.e-5
*EM_SOLVER_FEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ relTol maxIter solveType precond useLast
1.e-6 10000 2 1 1
*EM_SOLVER_BEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ relTol maxIter solveType precond uselast
1.e-6 10000 2 2 1
*EM_SOLVER_FEMBEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$ relTol maxIter
1.e-3 50
*EM_BEM_MATRICES
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
1 1.e-9
*EM_BEM_MATRICES
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
2 1.e-9
*EM_CIRCUIT_SOURCE
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R L C V0
1 11 50 -1.260e6 0.
$ ssidCurr ssidVltin ssidVltOt partID
1 1
$*EM_EXTERNAL_FIELD
$$--------1---------2---------3---------4---------5---------6---------7---------8
$$ fieldId fieldType FieldDef lcidX lcidY lcidZ
$ 1 1 10 21 22 23
*DEFINE_CURVE
21
0.,0.
10.,0.
*DEFINE_CURVE
22
0.,0.
10.,0.
*DEFINE_CURVE
23
0.,0.01
10.,0.01
$*EM_ROGO
$$--------1---------2---------3---------4---------5---------6---------7---------8
$$ rogoId setID setType currType
$ 1 1 1 1
$ 2 2 1 1
$ 3 3 1 1
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
2 2 32.78 0
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
1 1 32.78 0
*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_DATABASE_POINTOUT
1
1
*EM_POINT_SET
1
1,0.0,-55,0.5
2,0.0,-50,0.5
3,0.0,-45,0.5
4,0.0,-40,0.5
5,0.0,-35,0.5
6,0.0,-30,0.5
7,0.0,-25,0.5
8,0.0,-20,0.5
9,0.0,-15,0.5
10,0.0,-10,0.5
11,0.0,-05,0.5
12,0.0,0,0.5
13,0.0,05,0.5
14,0.0,10,0.5
15,0.0,15,0.5
16,0.0,20,0.5
17,0.0,25,0.5
18,0.0,30,0.5
19,0.0,35,0.5
20,0,40,0.5
21,0.0,45,0.5
22,0,50,0
23,0.0,55,0.5
$--------1---------2---------3---------4---------5---------6---------7---------8
$ ****************************************************************************
$ Misc: load curves,...
$ ****************************************************************************
*DEFINE_CURVE
100
0.,9.81e3
10.,9.81e3
*DEFINE_CURVE
3
0.,1.e-5
10.,1.e-5
*END
