2D Axisymmetric Inductive heating
2D axisymmetric Inductive heating example with periodic recomputation of EM fields. This example uses the monolithic solver to solve the FEM/BEM system which allows to use higher timesteps (See choice of numls).
https://www.dynaexamples.com/em/eddycurr/2dindheat
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
2D Axisymmetric Inductive heating
2D axisymmetric Inductive heating example with periodic recomputation of EM fields. This example uses the monolithic solver to solve the FEM/BEM system which allows to use higher timesteps (See choice of numls).
2D axisymmetric Inductive heating example with periodic recomputation of EM fields. This example uses the monolithic solver to solve the FEM/BEM system which allows to use higher timesteps (See choice of numls).
References :
[1]Approximate analytical solution for induction heating of solid cylinders
Todd A. Jankowski∗, Norma H. Pawley, Lindsey M. Gonzales, Craig A. Ross, James D. Jurney, Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM 87545, USA, Annual Mathematical Modelling 40 (2016) 2770-2782
*EM_2DAXI
*EM_BOUNDARY *EM_CIRCUIT *EM_CONTROL *EM_CONTROL_TIMESTEP
*EM_EOS_TABULATED1 *EM_MAT_001 *EM_ROTATION_AXIS
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_SOLVER_FEM
*EM_SOLVER_FEMBEM_MONOLITHIC *END
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LST-Ansys)
$
$
$X
$X------------------------------------------------------------------------------
$# UNITS: (S.I)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
2D Axisymmetric Inductive heating example
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CONTROL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$ emsol numls
2 &numls
*EM_CONTROL_TIMESTEP
1,&em_macrdt
*EM_SOLVER_FEMBEM_MONOLITHIC
0,1,1.e-12,1.e-3
*EM_ROTATION_AXIS
0,0,0,0,0,1,32
*EM_2DAXI
100,12,,,62,82
*EM_2DAXI
101,1,,,51,71
*EM_2DAXI
102,2,,,52,72
*EM_2DAXI
103,3,,,53,73
*EM_2DAXI
104,4,,,54,74
*EM_2DAXI
105,5,,,55,75
*EM_2DAXI
106,6,,,56,76
*EM_2DAXI
107,7,,,57,77
*EM_2DAXI
108,8,,,58,78
*EM_2DAXI
109,9,,,59,79
*EM_2DAXI
110,10,,,60,80
*EM_2DAXI
111,11,,,61,81
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM CIRCUIT $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
1 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
1 51 71
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
2 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
2 52 72
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
3 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
3 53 73
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
4 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
4 54 74
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
5 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
5 55 75
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
6 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
6 56 76
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
7 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
7 57 77
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
8 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
8 58 78
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
9 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
9 59 79
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
10 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
10 60 80
*EM_CIRCUIT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ circid circtype lcid R/F L/A C/t0 V0
11 11 &freq &
$ ssidCurr ssidVltin ssidVltOt partID
11 61 81
*EM_BOUNDARY
21,9
*EM_BOUNDARY
22,9
*EM_BOUNDARY
23,9
*EM_BOUNDARY
24,9
*EM_BOUNDARY
25,9
*EM_BOUNDARY
26,9
*EM_BOUNDARY
27,9
*EM_BOUNDARY
28,9
*EM_BOUNDARY
29,9
*EM_BOUNDARY
30,9
*EM_BOUNDARY
31,9
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ EM MAT and EOS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
1 2 &em_cond1
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$ em_mid mtype sigma eosId
2 4 &em_cond2 100
$ 2 4 &em_cond2
*EM_EOS_TABULATED1
100,201
$---+----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 5000 2 2 1&em_bemmtx
*EM_SOLVER_FEM
$# reltol maxit stype precon uselas ncyclbem
1e-4 5000 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
