Coupled structural thermal analysis example

Home » Knowledge » Implicit » Guideline Examples » Coupled structural thermal analysis example

A flat panel is fully constrained at one edge, see Figure 34. The initial temperature of the moving block (yellow in Figure 34) is 500 K, while all other parts have an initial temperature of 300 K. First, the moving block is placed on the panel (by gravity) for a duration of 100 s, then a motion is applied to close the contact with the fixed block (blue in Figure 34). For more details on thermal contacts, see Section 6.1.8. The temperature distribution after 1 s is shown in Figure 35 (cut view) and after 500 s in Figure 36. Note that the shell thickness is also visualized in these images, and that the fringe limits differ.

The example keyword file is thermal.key. For further examples of thermal analyses, see also Knowledge-base Thermal .

Example Files: thermal
NOTE! All examples are included in the Complete Guideline

Main Input

*KEYWORD
$===============================================================================
$ Geometry
$===============================================================================
*INCLUDE
geo_thermal002.key
$===============================================================================
$ Loads and boundary conditions
$===============================================================================
*BOUNDARY_SPC_SET
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
       109         0         1         1         1         1         1         1
*BOUNDARY_PRESCRIBED_MOTION_SET
       112         3         2      1001      -1.5                          100.
*DEFINE_CURVE_TITLE
Motion cube
      1001
0.,0.
10.,1.
1000.,1.
*LOAD_BODY_Z
      1002     9810.
*DEFINE_CURVE_TITLE
Gravity
      1002
0.,0.
1.,1.
1000.,1.
*INITIAL_TEMPERATURE_SET
$#    nsid      temp       loc   
       301 300.00000         0
*INITIAL_TEMPERATURE_SET
$ --- for the shells
       302 300.00000        -1
       302 300.00000         0
       302 300.00000         1
*INITIAL_TEMPERATURE_SET
       501 500.00000         0
$===============================================================================
$ Contacts
$===============================================================================
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_THERMAL_ID
       101Heater cube to shell
$#1
         4         1         3         3
$#2     FS
       0.1
$#3

$#4 --- thermal properties
     0.000     0.000 3000.0000 1.0000E-3 1.0000E-2       0.5 
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_THERMAL_ID
       102Support cube to shell
$#1
         7         1         3         3
$#2     FS
       0.1
$#3

$#4 --- thermal properties
     0.000     0.000 3000.0000 1.0000E-3 1.0000E-2       0.5 
$===============================================================================
$ Control cards
$===============================================================================
*INCLUDE
../INCLUDE_FILES/control_cards_nonlin.key
*CONTROL_IMPLICIT_DYNAMICS
$#   imass     gamma      beta    tdybir    tdydth    tdybur     irate     
         1  0.600000  0.380000                  1.        2.         1
*CONTROL_SOLUTION
$#    soln       nlq     isnan     lcint     
         2         0         0      1000
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas      
  500.0000       
*DEFINE_CURVE_TITLE
Implicit timesteps
700,
0.,0.01
1.,1.
100.,1.
110.,1.
500.,1.
$ --- thermal control
*INCLUDE
../INCLUDE_FILES/control_thermal_transient.key
*CONTROL_THERMAL_TIMESTEP
$#      ts       tip       its      tmin      tmax     dtemp      tscp      lcts
         0        1.     1.000
$ ---
*CONTROL_MPP_IO_NODUMP
$===============================================================================
$ Database cards
$===============================================================================
*INCLUDE
../INCLUDE_FILES/database_cards_static.key
*DATABASE_EXTENT_BINARY
         0         0         3         1         1         1         1         1
         0         1         0         0         0         0         0         0
$#3nintsld                                             therm
         1         0        1.         0         0         3                    
$#4           resplt
         0         0
$===============================================================================
$ Simulation title
$===============================================================================
*TITLE
180531, Thermal example with contact
*END