Tube expansion problem — Welcome to LS-DYNA Examples

Home » Knowledge Base » EM / Electromagnetics » Tutorial Eddy Currents » Tube expansion problem

This LS-DYNA simulation shows a simple Tube expansion example using the EM solver. A current is imposed in the coil which produces induced currents in the workpiece. This in turn generates Lorentz forces which cause the tube to expand. This problem is of medium size so running with 4 or more CPUs is recommended.

Description

This LS-DYNA simulation shows a simple Tube expansion example using the EM solver. A current is imposed in the coil which produces induced currents in the workpiece. This in turn generates Lorentz forces which cause the tube to expand. This problem is of medium size so running with 4 or more CPUs is recommended.

Animated Result

EM current density fringes

ezgif.comvideotogif.gif

Keywords

                                                               
*CONTROL_SOLUTION
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*CONTROL_THERMAL_SOLVER
*CONTROL_THERMAL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_TITLE
*EM_CIRCUIT
*EM_CIRCUIT_ROGO
*EM_CONTROL
*EM_CONTROL_TIMESTEP
*EM_MAT_001
*EM_OUTPUT
*EM_SOLVER_FEM
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*END
*INITIAL_TEMPERATURE_SET
*KEYWORD
*MAT_024
*MAT_RIGID
*MAT_THERMAL_ISOTROPIC
*PARAMETER
*PART
*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 R8.0.0 (or higher) with double precision 
$X
$X------------------------------------------------------------------------------
$# UNITS: (g/mm/s) 
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM Electromagnetic forming 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         1        10   &em_res   &em_ind   &em_cap    &em_v0
$ ssidCurr ssidVltin ssidVltOt    partID
         2         1         3
*EM_CIRCUIT_ROGO
$   rogoId     setId   setType  currType
         1         4         1         1
*EM_CIRCUIT_ROGO
$   rogoId     setId   setType  currType
         2         2         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_cond1          
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         2         4 &em_cond2    
$---+----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-6 
*EM_SOLVER_BEMMAT
$#   matid                                                                reltol
         2                                                                  1e-6      
*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      
*END

Figures

Cylinder Flow

Download

EM_example_tubeexpansion.zip