TEAM 28 Levitation problem
The TEAM 28 problem introduces an axisymmetric electrodynamic levitation device which consists of a conducting plate over two stranded coils. The aim is to determine the dynamic characteristics of the levitating plate (after some damped oscillations, the plate attains a stationary levitation height). A coupled solution of the electromagnetic and the mechanical problem is necessary.
References :
[1] S. K. G. L. W. R. Hans Karl, Joachim Fetzer, Description of team workshop problem 28: An electrodynamic levitation device, Institut fur Theorie der Elektrotechnik, Universitat Stuttgart, Pfa enwaldring 47, 70550 Stuttgart, Germany.
[2] F. Rapetti, An overlapping mortar element approach to coupled magneto-mechanical problems, Mathematics and Computers in Simulation, (2008).
[3] J. F. Stefan Kurz and G. Lehner, Threedimensional transient bem-fem coupled analysis of electrodynamic levitation problems, IEEE TRANSACTIONS ON MAGNETICS, 32.
Lorentz force fringes
*BOUNDARY_SPC_SET_BIRTH_DEATH *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_COORDINATE_SYSTEM *DEFINE_CURVE *DEFINE_TABLE *EM_2DAXI *EM_CIRCUIT_SOURCE *EM_CONTROL *EM_CONTROL_TIMESTEP *EM_MAT_001 *EM_ROTATION_AXIS *END *KEYWORD *MAT_RIGID *PARAMETER *PART *SECTION_SHELL *SECTION_SOLID *TITLE
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R9.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (g/mm/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM An Electrodynamic Levitation Device (TEAM 28)
*PARAMETER
R T_end 1.
R dt_plot 1.e-3
$
$--- em
R em_cond 34.
R em_dt 1.e-5
R em_freq 50.
R em_amp1 -11.52e6
R em_amp2 19.2e6
Iem_femmtx 25
Iem_bemmtx 25
$
$
$--- struc
$
R dt_struc 1.e-5
Rrho_struc 2687e-6
R E_struc 2.e11
R nu_struc 0.3
*EM_CONTROL
$# emsol numls dtinit dtmax t_init t_end ncyclfem ncyclbem
1 0 0.000 0.000 0.000 0.00 &em_femmtx&em_bemmtx
*EM_CONTROL_TIMESTEP
$# tstype dtconst
1 &1.e-5
*EM_MAT_001
$# mid mtype sigma eosid
1 4 &em_cond 0
*EM_MAT_001
$# mid mtype sigma eosid
2 1
*EM_MAT_001
$# mid mtype sigma eosid
3 1
*EM_CIRCUIT_SOURCE
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R L C V0
2 11 &em_freq &em_amp1
$ ssidCurr ssidVltin ssidVltOt partID
8 3
*EM_CIRCUIT_SOURCE
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R L C V0
1 11 &em_freq &em_amp2
$ ssidCurr ssidVltin ssidVltOt partID
7 2
*EM_ROTATION_AXIS
$--------1---------2---------3---------4---------5---------6---------7---------8
$ xpoint ypoint zpoint xaxis yaxis zaxis numSec
0. 0. 0. 0. 0. 1. 32
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$ partId segSetId direction 2Doption
1 -9 3 6 1
$ 1 9 6 3
*include
struc.k
*END