一個小的案例，

      以石黨勇的《基于ansysls/dyna8.1的顯式動力學》的3D 彈體穿透雙層鋼板的案例為基礎，加上熱力耦合，考慮塑性功轉化為熱能，以及摩擦生熱  單位制 cm-g-us-k

                                                                  溫度云圖

                                     Mises應力云圖

                               彈體上隨機選取的節點的溫度隨時間變化曲線

熱力耦合k文件設置的關鍵關鍵字

1.溫度屬性材料參數

*MAT_THERMAL_ISOTROPIC

$#    tmid       tro     tgrlc    tgmult      tlat      hlat    

         1      7.83       0.0       0.0     830.0       0.0

$#      hc        tc  

4.80000E-64.8000E-12

國際單位制比熱單位為J/kg/c=（m/s)^2/k  熱導率W/m/c=kg m/s^3/k

切換到當前單位制：比熱J/kg/c=10^-8*(cm/us)^2/K , 

                              熱導率W/m/c=kg m/s^3/k=10^3*10^2/10^18=10^(3+2-18) g*cm/us^3/

2.溫度初始化

*INITIAL_TEMPERATURE_SET

$#    nsid      temp       loc   

         0     273.0         0

3,溫度相關的求解設置

*CONTROL_SOLUTION

$#    soln       nlq     isnan     lcint     

         2         0         0         0

*CONTROL_THERMAL_NONLINEAR

$#  refmax       tol       dcp    lumpbc    thlstl    nlthpr    phchpn      

        20       0.0       1.0         0       0.0         0       0.0

*CONTROL_THERMAL_SOLVER

$#   atype     ptype    solver     cgtol       gpt    eqheat     fwork       sbc

         1         1         31.00000E-4         8       1.0       0.9       0.0

$#  msglvl    maxitr    abstol    reltol     omega    unused    unused       tsf

*CONTROL_THERMAL_TIMESTEP

$#      ts       tip       its      tmin      tmax     dtemp      tscp      lcts

         1       1.0      0.05      1E-4         2      60.0       0.5         0

*CONTROL_TIMESTEP

$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st

       0.0       0.9         0       0.0       0.0         0         0         0

4考慮熱傳遞的接觸設置

*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_THERMAL

$#     cid                                                                 title

$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr

        3          1         3         3         0         0         0         0

$#      fs        fd        dc        vc       vdc    penchk        bt        dt

       0.5       0.5           

$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf

       1.0       1.0                          

$#      cf      frad       htc      lmin      lmax    ftoslv    bc_flg      algo

                       4.8E-12      0.01      0.05       0.5

5.接觸控制

*CONTROL_CONTACT

$#  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass

        10       0.0         1         0         1         0         1         0

$#  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj

        

$#   sfric     dfric       edc       vfc        th     th_sf    pen_sf      

       

$#  ignore    frceng   skiprwg    outseg   spotstp   spotdel   spothin       

         0         1         0         0         0         0       0.0

$#    isym    nserod    rwgaps    rwgdth     rwksf      icov    swradf    ithoff

       

$#  shledg    pstiff    ithcnt    tdcnof     ftall    unused    shltrw      

        

$$