This problem demonstrates using LS-DYNA to solve for the unconstrained expansion of a block due to heating. The model consists of one 8 node brick element at an initial temperature of 10. The brick material is given a volumetric thermal generation rate. Explicit time integration is used for the structural calculations and implicit time integration is used for the thermal calculations. Implicit time integration is unconditionally stable and, therefore, a larger thermal time step can be taken.
*CONTROL_SOLUTION *CONTROL_TERMINATION *CONTROL_THERMAL_SOLVER *CONTROL_THERMAL_TIMESTEP *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_GLSTAT *ELEMENT_SOLID *END *INITIAL_TEMPERATURE_NODE *KEYWORD *MAT_ELASTIC_PLASTIC_THERMAL *MAT_THERMAL_ISOTROPIC *NODE *PART *SECTION_SOLID *TITLE
$ Thermal Expansion of a Block
$
*KEYWORD
$
$=============================CONTROL DEFINITIONS ==============================
$
*TITLE
thermal expansion of a block
*CONTROL_SOLUTION
2
*CONTROL_THERMAL_SOLVER
1 0 1
*CONTROL_TIMESTEP
.01
*CONTROL_THERMAL_TIMESTEP
0 1. .1
*CONTROL_TERMINATION
3.
*DATABASE_BINARY_D3PLOT
.01
*database_glstat
.03
$
$============================== PART DEFINITIONS ===============================
$
*PART
$ PID SECID MID TMID
slab
1 1 1 1
$
$============================= SECTION PROPERTIES ==============================
$
*SECTION_SOLID
$ SECID ELFORM
1 1
$
$====================== MECHANICAL MATERIAL PROPERTIES =========================
$
*MAT_ELASTIC_PLASTIC_THERMAL
1 1.
0. 10. 20. 30. 40. 50.
1.e+10 1.e+10 1.e+10 1.e+10 1.e+10 1.e+10
.3 .3 .3 .3 .3 .3
0.e-06 20.e-07 40.e-07 60.e-07 80.e-07 100.e-07
1.e+20 1.e+20 1.e+20 1.e+20 1.e+20 1.e+20
0. 0. 0. 0. 0. 0.
$
$======================== THERMAL MATERIAL PROPERTIES ==========================
$
*MAT_THERMAL_ISOTROPIC
1 1. 0 10.
1. 1.
$
$============================= NODE DEFINTIONS =================================
$
*NODE
1 0. 0. 0. 7
2 1. 0. 0. 5
3 1. 1. 0. 3
4 0. 1. 0. 6
5 0. 0. 1. 4
6 1. 0. 1. 2
7 1. 1. 1. 0
8 0. 1. 1. 1
$
$============================ ELEMENT DEFINITIONS ==============================
$
*ELEMENT_SOLID
1 1 1 2 3 4 5 6 7 8
$
$============================ ELEMENT DEFINITIONS ==============================
$
*INITIAL_TEMPERATURE_NODE
1 10.
2 10.
3 10.
4 10.
5 10.
6 10.
7 10.
8 10.
*END
