Axisymmetric problem — Welcome to LS-DYNA Examples

Home » Knowledge Base » EM / Electromagnetics » Tutorial Eddy Currents » Axisymmetric problem

The R9 version introduced many new advanced capabilities to the axisymmetric solver. For example, it is now possible to simulate pancake or helicoid type of coils by connecting the different turns using EM_CIRCUIT_CONNECT therefore expanding the range of applicable EM forming or bending problems that can be solved using the axisymmetric solver.

Description

The R9 version introduced many new advanced capabilities to the axisymmetric solver. For example, it is now possible to simulate pancake or helicoid type of coils by connecting the different turns using EM_CIRCUIT_CONNECT therefore expanding the range of applicable EM forming or bending problems that can be solved using the axisymmetric solver. Run times are tremendously reduced which allows the set up of advanced optimization studies in combination with LS-Opt.

Animated Result

Temperature fringes

ezgif.comvideotogif.gif

Keywords

*CONTACT_SURFACE_TO_SURFACE    
*CONTROL_CONTACT                                                        
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_FUNCTION
*DEFINE_CURVE_TITLE
*DEFINE_TABLE
*EM_2DAXI
*EM_DATABASE_PARTDATA
*EM_CIRCUIT
*EM_CIRCUIT_CONNECT
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_MAT_001
*EM_OUTPUT
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_ROTATION_AXIS
*END
*HOURGLASS
*KEYWORD
*MAT_ELASTIC
*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 Electromagnetic axisymmetric problem
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM CONTROL                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$    emsol        
         1         
*EM_CONTROL_TIMESTEP
$#  tstype    dtcons      lcid                            
         1    &em_dt
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM CIRCUIT                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CIRCUIT
$   circid  circtype      lcid         R         L         C        V0
         1         3             &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
         6         4         5
*EM_CIRCUIT
$   circid  circtype      lcid         R         L         C        V0
         2         3             &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
         9         7         8
*EM_CIRCUIT
$   circid  circtype      lcid         R         L         C        V0
         3         3             &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
        12        10        11
*EM_CIRCUIT
$   circid  circtype      lcid         R         L         C        V0
         4         3             &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
        15        13        14 
*EM_CIRCUIT
$   circid  circtype      lcid         R         L         C        V0
         5         3             &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
        18        16        17 
*EM_CIRCUIT
$   circid  circtype      lcid         R         L         C        V0
         6         3             &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
        21        19        20         
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   connId  connType  circuit1  circuit2        c1        c2
         1         1         1         2        1.       -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   connId  connType  circuit1  circuit2        c1        c2
         2         1         1         3        1.       -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   connId  connType  circuit1  circuit2        c1        c2
         3         1         1         4        1.       -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   connId  connType  circuit1  circuit2        c1        c2
         4         1         1         5        1.       -1.
*EM_CIRCUIT_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   connId  connType  circuit1  circuit2        c1        c2
         4         1         1         6        1.       -1.
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM MAT and EOS                                   $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         1         2  &em_cond          
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         2         1     
*EM_MAT_001            
$   em_mid     mtype     sigma       eos              
         3         4  &em_cond  
*EM_MAT_001            
$   em_mid     mtype     sigma       eos              
         4         4  &em_cond  
*EM_MAT_001            
$   em_mid     mtype     sigma       eos              
         5         4  &em_cond  
*EM_MAT_001            
$   em_mid     mtype     sigma       eos              
         6         4  &em_cond  
*EM_MAT_001            
$   em_mid     mtype     sigma       eos              
         7         4  &em_cond  
*EM_MAT_001            
$   em_mid     mtype     sigma       eos              
         8         4  &em_cond        
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM 2D AXI                                        $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_ROTATION_AXIS
$   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 startSsid   endSsid    numSec     new2D
         3         3         3         1         1         2        32         1         
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$   partId  segSetId direction  2Doption startSsid   endSsid    numSec     new2D
         1         6         3         1         4         5        32        
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$   partId  segSetId direction  2Doption startSsid   endSsid    numSec     new2D
         4         9         3         1         7         8        32         
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$   partId  segSetId direction  2Doption startSsid   endSsid    numSec     new2D
         5        12         3         1        10        11        32         
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$   partId  segSetId direction  2Doption startSsid   endSsid    numSec     new2D
         6        15         3         1        13        14        32        
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$   partId  segSetId direction  2Doption startSsid   endSsid    numSec     new2D
         7        18         3         1        16        17        32        
*EM_2DAXI
$--------1---------2---------3---------4---------5---------6---------7---------8
$   partId  segSetId direction  2Doption startSsid   endSsid    numSec     new2D
         8        21         3         1        19        20        32         
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM SOLVER                                        $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_SOLVER_BEMMAT
$#   matid                                                                reltol
         1                                                                  1e-9 
*EM_SOLVER_BEMMAT
$#   matid                                                                reltol
         2                                                                  1e-9      
*EM_SOLVER_BEM
$#  reltol     maxit     stype    precon    uselas  ncyclbem
      1e-6      1000         2         2         1&em_bemmtx
*EM_SOLVER_FEM
$#  reltol     maxit     stype    precon    uselas  ncyclbem
      1e-3      1000         1         1         1&em_femmtx
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM OUTPUT                                        $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 
*EM_OUTPUT
$     matS      matF      solS      solF      mesh
         2         2         2         2         0     
*EM_DATABASE_PARTDATA
$#  outlvl
         1 
*END

Figures

Cylinder Flow

Download