Thermal

This example shows the use of *LOAD_THERMAL_VARIABLE keyword in transient analysis to thermally preload a structure.

*CONTROL_CONTACT
*CONTROL_CPU
*CONTROL_ENERGY
*CONTROL_OUTPUT
*CONTROL_SHELL
*CONTROL_STRUCTURED
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DATABASE_EXTENT_BINARY
*DEFINE_CURVE
*ELEMENT_SHELL
*ELEMENT_SOLID
*END
*KEYWORD
*LOAD_THERMAL_CONSTANT
*LOAD_THERMAL_LOAD_CURVE
*LOAD_THERMAL_VARIABLE
*MAT_ELASTIC
*MAT_ELASTIC_PLASTIC_THERMAL
*NODE
*PART
*SECTION_SHELL
*SECTION_SOLID
*TITLE

  *KEYWORD
  $ shell is mat4; solid is mat1
  $  Use *load_thermal_variable in transient analysis to thermally preload 
  $   structure.
  $ ***********************************************************************
  $  3 forms of temperature loading are given below.
  $  (none are suited for thermally prestressing with DRelaxtion
  $    and holding thermal stress during transient)
  $ ***********************************************************************
  $ ********  METHOD 1 *****************
  $ temperature for solids is relative to T = 0 
  $  but temperature for shells is relative to T(t=0) (pre-LS941 only).
  $   BUG in SHELLS pre-LS941  !!!!!!!!!!
  $    itemp = -2
  $*LOAD_THERMAL_constant
  $0
  $25
  $ ***********************************************************************
  $ ********  METHOD 2 *****************
  $  temperature is relative to T(t=0) as documented  (OK)
  $    itemp = 2
  $*load_thermal_load_curve
  $2
  $ ***********************************************************************
  $ ********  METHOD 3 *****************
  $  temperature is relative to T(t=0)   (OK) 
  $    itemp = -9999
  *load_thermal_variable
  0
  1.0, 0, 2
  $ ***********************************************************************
  $  temperature load curve
  *DEFINE_CURVE
           2         
  $  material feels no load if temperature is constant
  $      0.00000000E+00      2.50000000E+01
        0.00000000E+00      0.00000000E+01
        1.00000005E-03      2.50000000E+01
        1.00000000E+02      2.50000000E+01
  *section_solid
  2
  *part
  
  2,2,2
  *element_solid
  $ shared nodes
  1,2,4,9,8,3,59,54,53,58
  $ independent nodes
  $1,2,104,109,108,103,159,154,153,158
  *control_structured
  *TITLE
  tempcyl.vari
  *CONTROL_TERMINATION
   1.00000-3         0 0.0000000         0 0.0000000
  *CONTROL_TIMESTEP
   0.0000000 0.0000000         0 0.0000000 0.0000000         0         0         0
  *CONTROL_SHELL
   0.0000000         0         0         0         0         0         0
  *CONTROL_CONTACT
   0.0000000 0.0000000         0         0         0         0         0
           0         0         0         0 0.0000000         0         0         0
  *CONTROL_OUTPUT
           0         0         0         0 0.0000000         0         0
  *CONTROL_ENERGY
           1         2         1         1
  *DATABASE_BINARY_D3PLOT
   1.00000-5                   0
  *DATABASE_EXTENT_BINARY
           0         0         0         1         1         1         1         0
           0         0         0         0         0         0
  *CONTROL_CPU
   0.0000000
  *mat_elastic
  2,1e-3,1e7,.3
  *MAT_ELASTIC_PLASTIC_THERMAL
           1 1.00000-3
  -1000.0000 0.0000000 1000.0000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
   1.00000+7 1.00000+7 1.00000+7 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
   0.3000000 0.3000000 0.3000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
   2.00000-4 2.00000-4 2.00000-4 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
   1.00000+6 1.00000+6 1.00000+6 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
   1.00000+6 1.00000+6 1.00000+6 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
  *SECTION_SHELL
           1         2 0.0000000 3.0000000 0.0000000 0.0000000
   0.5000000 0.5000000 0.5000000 0.5000000 0.0000000
  *PART
  material type # 4  (Thermo-Elastic-Plastic)                                     
           1         1         1         0         0         0         0         0
  *NODE
  $
  *ELEMENT_SHELL
  $
  *END