TEAM 2 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, excitation, etc., and each backed by a real laboratory device, on which measurements can be made.
Problem 2 [1] of the TEAM workshops for eddy current code comparison is a conducting cylinder sphere in a uniform magnetic field (in the Y direction). The cylinder is modelled here as a very long 3D part. The field is spatially uniform but varies with a given frequency. The objective is to measure the magnetic flux amplitude and phase inside and outside the cylinder at different instants and compare to the analytical solution and the frequency Eddy current solver is used.
[1]INFINITE CYLINDER IN A UNIFORM SINUSOIDAL FIELD (COMPARISON OF RESULTS, PROBLEM 2). C.R.I. EMSON. COMPEL—The International Journal for Computation and Mathematics in Electrical and Electronic Engineering Vol. 7, No. 1/2, 29-45 © 1988 James & James Science Publishers Ltd
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_EXTERNAL_FIELD
*EM_DATABASE_POINTOUT
*EM_MAT_001
*EM_OUTPUT
*EM_POINT_SET
*EM_SOLVER_FEM
*EM_SOLVER_FEMBEM_MONOLITHIC
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LST-Ansys)
$
$
$X
$X------------------------------------------------------------------------------
$# UNITS: kg/m/s
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R15.0.0 (or higher) with double precision
*KEYWORD
*TITLE
TEAM 2 Problem
$ ****************************************************************************
$ mesh
$ ****************************************************************************
*INCLUDE
mesh.k
*PARAMETER
R Tend 2.e-4
R dtplot 1.e-4
R dt 1.e-5
$
$--- em
R em_cond2.538071e7
R em_dt 1.e-4
R B0 0.1
R Freq 60.
$ *****************************mo***********************************************
$ mechanics
$ ****************************************************************************
$--------1---------2---------3---------4---------5---------6---------7---------8
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
&Tend
*CONTROL_TIMESTEP
3
*DATABASE_BINARY_D3PLOT
&dtplot
*PART
1/4 Tube
$ pid secid mid
1 1 1
*PART
3/4 Tube
$ pid secid mid
2 1 2
*SECTION_SOLID
$ sid elform
1 1
*MAT_ELASTIC
2,7000e,200e9,0.3
*MAT_RIGID
1,7000e,200e9,0.3
$# cmo con1 con2
1.000000 7 7
$ ****************************************************************************
$ electromagnetism
$ ****************************************************************************
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls emdtinit emdtmax emtinit emtend ncyclFem ncyclBem
4 &Freq 5000 5000
*EM_CONTROL_COUPLING
,2
*EM_DATABASE_PARTDATA
1
*EM_CONTROL_TIMESTEP
1,&em_dt
*DEFINE_CURVE
3
0,&dt
1000,&dt
*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
1.e-6 10000 2 4 1
*EM_SOLVER_FEMBEM_MONOLITHIC
,,1.e-16,1.e-5,10000,,
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
1 1.e-16
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
2 1.e-16
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matId relTol
3 1.e-16
*define_curve
189
0.,&B0
1.,&B0
*define_curve
190
0.0,0.0
1.0,0.0
*EM_EXTERNAL_FIELD
$--------1---------2---------3---------4---------5---------6---------7---------8
$ fieldId fieldType FieldDef lcidX lcidY lcidZ
1 1 1 190 189 190
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
1 2 &em_cond
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
2 2 &em_cond
*EM_OUTPUT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ matS matF solS solF mesh memory timing d3plotAsc
3,3,3,3,,
$4,4,4,4
$ mf2 gmv d3plotFor timeHist
0 0 0
*EM_DATABASE_POINTOUT
1
1
*EM_POINT_SET
1
1,0.01,0,0.0
2,0.02,0,0.0
3,0.03,0,0.0
4,0.04,0,0.0
5,0.05,0,0.0
6,0.05842,0,0.0
7,0.06,0,0.0
8,0.065,0,0.0
9,0.06858,0,0.0
10,0.07,0,0.0
11,0.075,0,0.0
12,0.1,0,0.0
13,0.15,0,0.0
21,0.009397,0.003420,0.0
22,0.018794,0.006840,0.0
23,0.028191,0.010261,0.0
24,0.037588,0.013681,0.0
25,0.046985,0.017101,0.0
26,0.054897,0.019981,0.0
27,0.056382,0.020521,0.0
28,0.061080,0.022231,0.0
29,0.064444,0.023456,0.0
30,0.065778,0.023941,0.0
31,0.070477,0.025652,0.0
32,0.093969,0.034202,0.0
33,0.140954,0.051303,0.0
41,0.007071,0.007071,0.0
42,0.014142,0.014142,0.0
43,0.021213,0.021213,0.0
44,0.028284,0.028284,0.0
45,0.035355,0.035355,0.0
46,0.041309,0.041309,0.0
47,0.042426,0.042426,0.0
48,0.045962,0.045962,0.0
49,0.048493,0.048493,0.0
50,0.049497,0.049497,0.0
51,0.053033,0.053033,0.0
52,0.070711,0.070711,0.0
53,0.106066,0.106066,0.0
*END
