如何用lsdyna做摩擦生热呢？
见帖的各位专家老师好，我在做一个二维轴对称模型的摩擦生热案例，参考了蓝牙老师3d模型的关键字进行了计算，调通了以后不能执行，但是找不到问题，这是热力部分的关键字参数：
*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
*MAT_THERMAL_ISOTROPIC
$# tmid tro tgrlc tgmult tlat hlat
2 2.4 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
*COOL_THERMAL_NONLINEAR
$# refmax tol dcp lumpbc thlstl nlthpr phchpn
20 0.0 1.0 0 0.0 0 0.0
*COOL_THERMAL_SOLVER
$# atype ptype solver cgtol gpt eqheat fwork sbc
1 1 31.00000E-4 1 1.0 0.9 0.0
$# msglvl maxitr abstol reltol omega unused unused tsf
0 5001.0000E-101.00000E-4 1.0 1.0
*COOL_TIMESTEP
$$ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.9 0 0.0 0.0 0 0 0
$# dt2msf dt2mslc imscl unused unused rmscl unused ihdo
0.0 0 0 0.0 0
*COOL_SOLUTION
$# soln nlq isnan lcint
2 0 0 0
*CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE_THERMAL_TITLE
$# cid title
1
$# sids sidm sfact freq fs fd dc
2 1 0.1 2 0.0 0.0 0.0
$# tbirth tdeath sos som nds ndm cof init
0.01.00000E20 1.0 1.0 0 0 0 0
$# k rad h lmin lmax chlm bc_flag
0.0 0.04.8000E-12 0.01 0.05 1.0 0
$# vc vdc ipf slide istiff tiedgap igapcl tietyp
0.0 10.0 0 0 0 0.0 0 0
$# sldsos sldsom tdpen
0.0 0.0 0.0
