Structure-thermal coupled analysis I

In this example a steel frame is subjected to a thermal load which results in thermal expansion of the structure. The simulation is set up as a structure-thermal coupled analysis using (static) implicit time integration. For the definition of the thermal load the keyword *BOUNDARY_TEMPERATURE is used and temperature fields are calculated by the thermal solver. The temperature fields are handend over to the structural solver which calculates thermal expansion due to the temperature field. The temperature of the steel frame is initially at 20.0 °C and is ramped up to 200 °C. Thermal material properties are defined with *MAT_THERMAL_ISOTROPIC for the thermal solver. The temperature dependent material properties for the strucural solver, like thermal expansion coefficient, are defined with *MAT_ELASTIC_PLASTIC_THERMAL.

*BOUNDARY_SPC_SET
*BOUNDARY_TEMPERATURE_SET
*CONTROL_ACCURACY
*CONTROL_IMPLICIT_GENERAL
*CONTROL_IMPLICIT_SOLUTION
*CONTROL_SOLUTION
*CONTROL_TERMINATION
*CONTROL_THERMAL_SOLVER
*CONTROL_THERMAL_TIMESTEP
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DATABASE_GLSTAT
*DATABASE_MATSUM
*DEFINE_CURVE
*END
*INCLUDE
*INITIAL_TEMPERATURE_SET
*KEYWORD
*MAT_ELASTIC_PLASTIC_THERMAL
*MAT_THERMAL_ISOTROPIC
*PART
*SECTION_SHELL
*SET_NODE_GENERAL
*SET_NODE_LIST_TITLE
*TITLE

*KEYWORD
$
$=============================CONTROL DEFINITIONS ==============================
$
*TITLE
Steel frame - structure-thermal coupled analysis
*CONTROL_SOLUTION
$     soln       nlq     isnan     lcint     lcacc     ncdcf
         2
*CONTROL_THERMAL_SOLVER
$    atype     ptype    solver     cgtol       gpt    eqheat     fwork       sbc
         1         1        11
$   msglvl    maxitr    abstol    reltol     omega                           tsf

*CONTROL_TIMESTEP
$   dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
                 0.9
*CONTROL_THERMAL_TIMESTEP
$       ts       tip       its      tmin      tmax     dtemp      tscp      lcts
         1
*CONTROL_TERMINATION
$   endtim    endcyc     dtmin    endeng    endmas     nosol
       1.0
*CONTROL_ACCURACY
$      osu       inn    pidosu      iacc
         1         4                   1
$--------1---------2---------3---------4---------5---------6---------7---------8
$ IMPLICIT CONTROL CARDS
$--------1---------2---------3---------4---------5---------6---------7---------8
*CONTROL_IMPLICIT_GENERAL
$   imflag       dt0    imform      nsbs       igs     cnstn      form    zero_v
         1       0.1                            
*CONTROL_IMPLICIT_SOLUTION
$   nsolvr    ilimit    maxref     dctol     ectol     rctol     lstol    abstol
        12                        0.0010   0.01000                       1.0E-20
$    dnorm    diverg     istif   nlprint    nlnorm
         0         0         0         1         4
$
$================================ OUTPUTDATA ===================================
$
*DATABASE_BINARY_D3PLOT
$       dt
       0.1
*DATABASE_GLSTAT
$       dt
       0.1
*DATABASE_MATSUM
$       dt
       0.1
$
$============================== PART DEFINITIONS ===============================
$
*PART
$ title
column_1
$      pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         1         1                                                 1
*PART
$ title
bar
$      pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2         1         1                                                 1
*PART
$ title
column_2
$      pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         1         1                                                 1
$
$============================= SECTION PROPERTIES ==============================
$
*SECTION_SHELL
$      sid    elform      shrf       nip     propt   qr/irid     icomp
         1       -16     0.833       5.0
$       t1        t2        t3        t4
       3.0       3.0       3.0       3.0
$
$====================== MECHANICAL MATERIAL PROPERTIES =========================
$
*MAT_ELASTIC_PLASTIC_THERMAL
$      mid        ro
         1   7.85e-9
$       t1        t2
         0       400
$       e1        e2
    210000    210000
$      pr1       pr2
       0.3       0.3
$   alpha1    alpha2
     12e-6     12e-6
$    sigy1     sigy2
       500       500
$    etan1     etan2
     10000     10000
$
$======================== THERMAL MATERIAL PROPERTIES ==========================
$
*MAT_THERMAL_ISOTROPIC
$     tmid       tro     tgrlc    tgmult      tlat      hlat
         1   7.85e-9
$       hc        tc		 
    4.60e8      40.0
$
$================= THERMAL BOUNDARY AND INITIAL CONDITIONS =====================
$
*INITIAL_TEMPERATURE_SET
$  nid/sid     tlcid       loc
         0        20
$
*BOUNDARY_TEMPERATURE_SET
$  nid/sid     tlcid     tmult       loc    tdeath    tbirth
         1       101       1.0
$
*DEFINE_CURVE
$     lcid
       101
$ time,temp
0,20
1,200
$
$====================== MECHANICAL BOUNDARY CONDITIONS =========================
$
*BOUNDARY_SPC_SET
$ nid/nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
         2                   1         1         1         1         1         1
*BOUNDARY_SPC_SET
$ nid/nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
         3                   1         1         1         1         1         1
$=======================NODE/ELEMENT/SET/SEGMENT DEFINTIONS ====================
*INCLUDE
geo.k
$
*SET_NODE_GENERAL
$      sid
         1
$   option
       ALL
*SET_NODE_LIST_TITLE
fixing_1
$      sid       da1       da2       da3       da4    solver
         2 
$     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         4         1         3         2                                        
*SET_NODE_LIST_TITLE
fixing_2
$      sid       da1       da2       da3       da4    solver
         3 
$     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
       161       163       164       162
$
$--------1---------2---------3---------4---------5---------6---------7---------8
$
*END