记录一下蓝牙老师的资料，已做备份。考虑侵彻过程中温度的来源主要来自两个方面，塑性功和摩擦热导致温度上升。
先看一下蓝牙老师的仿真结果，如图所示。
热力耦合中涉及到的主要关键字如下：
带温度属性的材料参数，参考如下
*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 热导率 J /m/c=kg m/s^3/k
仿真开始时的温度初始化
*INITIAL_TEMPERATURE_SET
$# nsid temp loc
0 273.0 0
跟温度相关的求解设置
*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
考虑热传导的接触设置
*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
接触的控制