TEAM 28 Levitation problem — Welcome to LS-DYNA Examples

Home » Knowledge Base » EM / Electromagnetics » Tutorial Eddy Currents » 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.

Description

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.

Animated Result

Lorentz force fringes

ezgif.comvideotogif.gif

Keywords

*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

Reduced Input

$-----------------------------------------------------------------------------
$
$ 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

Figures

Cylinder Flow

Download