见帖的各位专家老师好，我在做一个二维轴对称模型的摩擦生热案例，参考了蓝牙老师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