电磁热耦合，计算第一步就出错，求帮助~

huyang940

第一次做电磁热耦合，钢坯的几何是自己做的，其他设置步骤是根据教材来的，但是提交计算后第一步的时间增量怎么会是10e-30这个数量级的，然后迭代到10e-19后就报错了。  
用的钢坯是截面有阶梯状的圆柱，所以线圈用两个分开的矩形模拟了。以下是命令流。求大虾帮助~  
/FILNAME,11,0   
/REPLOT,RESIZE  
*SET,freq,150000
*SET,pi,4*atan(1)   
*SET,cond,.392e7
*SET,muzero,4e-7*pi
*SET,mur,200
*SET,skind,sqrt(1/(pi*freq*cond*muzero*mur))
*SET,ftime,3
*SET,finc,.05   
*SET,time,0
*SET,delt,.01   
/PREP7  
!*  
!*  
ET,1,PLANE53
!*  
ET,2,PLANE53
!*  
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,5,0
KEYOPT,1,7,0
!*  
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,1
KEYOPT,2,4,0
KEYOPT,2,5,0
KEYOPT,2,7,0
!*  
!*  
EMUNIT,MKS  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,1,,1
MPDE,ALL,1  
TBDE,ALL,1  
MPTEMP,,,,,,,   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,1,,1
MPTEMP,,,,,,,,  
MPTEMP,1,25.5   
MPTEMP,2,160
MPTEMP,3,291.5  
MPTEMP,4,477.6  
MPTEMP,5,635
MPTEMP,6,698
MPTEMP,7,709
MPTEMP,8,720.3  
MPTEMP,9,742
MPTEMP,10,761   
MPTEMP,11,1000  
MPDATA,MURX,2,,200  
MPDATA,MURX,2,,190  
MPDATA,MURX,2,,182  
MPDATA,MURX,2,,161  
MPDATA,MURX,2,,135  
MPDATA,MURX,2,,104  
MPDATA,MURX,2,,84   
MPDATA,MURX,2,,35   
MPDATA,MURX,2,,17   
MPDATA,MURX,2,,1
MPDATA,MURX,2,,1
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,125
MPTEMP,3,250
MPTEMP,4,375
MPTEMP,5,500
MPTEMP,6,625
MPTEMP,7,750
MPTEMP,8,875
MPTEMP,9,1000   
MPDATA,RSVX,2,,1.84e-7  
MPDATA,RSVX,2,,2.72e-7  
MPDATA,RSVX,2,,3.84e-7  
MPDATA,RSVX,2,,5.12e-7  
MPDATA,RSVX,2,,6.56e-7  
MPDATA,RSVX,2,,8.24e-7  
MPDATA,RSVX,2,,1.032e-6
MPDATA,RSVX,2,,1.15e-6  
MPDATA,RSVX,2,,1.2e-6   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,27
MPTEMP,3,127
MPTEMP,4,327
MPTEMP,5,527
MPTEMP,6,727
MPTEMP,7,765
MPTEMP,8,765.001
MPTEMP,9,927
MPDATA,ENTH,2,,0
MPDATA,ENTH,2,,91609056
MPDATA,ENTH,2,,453285756
MPDATA,ENTH,2,,1.2748e9
MPDATA,ENTH,2,,2.2519e9
MPDATA,ENTH,2,,3.3396e9
MPDATA,ENTH,2,,3.548547e9   
MPDATA,ENTH,2,,3.548556e9   
MPDATA,ENTH,2,,4.352e9  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,730
MPTEMP,3,930
MPTEMP,4,1000   
MPDATA,KXX,2,,60.64
MPDATA,KXX,2,,29.5  
MPDATA,KXX,2,,28
MPDATA,KXX,2,,28
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EMIS,2,,0.68
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,3,,1
!*  
/REPLOT,RESIZE  
K,1,0,0,0,             !创建节点
K,2,0.025,0,0,  
K,3,0.025,0.015,0,  
K,4,0.015,0.015,0,  
K,5,0.015,0.04,0,   
K,6,0,0.04,0,   
K,7,0.03,0,0,   
K,8,0.03,0.02,0,
K,9,0.02,0.02,0,
K,10,0.02,0.04,0,   
/REPLOT
LSTR,       1,       2     !连线
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       5  
LSTR,       5,       6  
LSTR,       6,       1  
LSTR,       2,       7  
LSTR,       7,       8  
LSTR,       8,       9  
LSTR,       9,      10  
LSTR,      10,       5  
KPLOT   
/PNUM,KP,1  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
LPLOT   
FLST,2,6,4      !线围成面（钢坯料）
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
FITEM,2,5   
FITEM,2,6   
AL,P51X      
FLST,2,8,4  
FITEM,2,7   
FITEM,2,8   
FITEM,2,9   
FITEM,2,10  
FITEM,2,11  
FITEM,2,4   
FITEM,2,3   
FITEM,2,2   
AL,P51X
RECTNG,0.020,0.022,0.020,0.040,     
RECTNG,0.03,0.032,0,0.02,   
RECTNG,0.022,0.02451,0.02,0.04,
RECTNG,0.032,0.03451,0,0.02,
/FOC,   1,  0.255791421483E-01,  0.162913035094E-01,   0.00000000000
/REPLO  
K,111,0.07,0,0,
/FOC,   1,  0.405964673151E-01,  0.184366356761E-01,   0.00000000000
/REPLO  
KPLOT   
APLOT   
KPLOT   
K,112,0.07,0.04,0,  
LSTR,      24,     111  
LSTR,     111,     112  
LSTR,     112,      21  
LSTR,      21,      20  
LSTR,      20,      25  
LSTR,      25,      24  
FLST,2,6,4  
FITEM,2,28  
FITEM,2,29  
FITEM,2,30  
FITEM,2,21  
FITEM,2,31
FITEM,2,25
AL,P51X
APLOT                       !电磁分析几何创建完毕
FLST,2,7,5,ORDE,2            !布尔操作
FITEM,2,1   
FITEM,2,-7  
AGLUE,P51X  
NUMCMP,AREA                   !重新编号
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
/FOC,   1,  0.369027718125E-01,  0.185463494039E-01,   0.00000000000
/REPLO  
/FOC,   1,  0.358422057771E-01,  0.181440657353E-01,   0.00000000000
/REPLO  
ASEL,S, , ,       7         !对钢坯表面设置单元大小
/REPLOT
FLST,2,1,5,ORDE,1   
FITEM,2,7   
AESIZE,P51X,0.0005,
ASEL,S, , ,       1           !钢坯内部单元大小
/REPLOT
FLST,2,1,5,ORDE,1   
FITEM,2,1   
AESIZE,P51X,0.0015,   
ALLSEL,ALL  
/REPLOT
FLST,5,4,5,ORDE,2            !空气和线圈单元大小
FITEM,5,2   
FITEM,5,-5  
ASEL,S, , ,P51X
/REPLOT
FLST,2,4,5,ORDE,2   
FITEM,2,2   
FITEM,2,-5  
AESIZE,P51X,0.001,
  
ASEL,S, , ,       6          !外部空气单元大小
/REPLOT
FLST,2,1,5,ORDE,1   
FITEM,2,6   
AESIZE,P51X,0.0035,

ALLSEL,ALL  
/REPLOT
FLST,5,2,5,ORDE,2    !单元类型设置
FITEM,5,1   
FITEM,5,7   
ASEL,S, , ,P51X
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,7   
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,4   
FITEM,5,-5  
ASEL,S, , ,P51X
FLST,5,2,5,ORDE,2   
FITEM,5,4   
FITEM,5,-5  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       3, ,   2,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,2   
FITEM,5,-3  
ASEL,S, , ,P51X
/REPLOT
FLST,5,2,5,ORDE,2   
FITEM,5,2   
FITEM,5,-3  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   2,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ASEL,S, , ,       6
CM,_Y,AREA  
ASEL, , , ,       6
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   2,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ALLSEL,ALL  
/REPLOT
CM,_Y,AREA  
ASEL, , , ,       7
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,2
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,5,4,5,ORDE,2   
FITEM,5,2   
FITEM,5,-5  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,2
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,6   
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,2
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2       !网格划分完毕
!*  
SAVE
ALLSEL,ALL  

NSEL,S,LOC,X,0              !电磁边界条件
FLST,2,55,1,ORDE,4  
FITEM,2,5098
FITEM,2,5132
FITEM,2,5152
FITEM,2,-5204   
D,P51X, ,0, , , ,AZ        !电磁边界条件设置完毕
               
ESEL,S,MAT,,3   !电流密度载荷设置
/REPLOT
FLST,2,360,2,ORDE,2
FITEM,2,1263
FITEM,2,-1622   
BFE,P51X,JS,1, , ,79.55e6, ,
ALLSEL,ALL  
/REPLOT
!*  

ET,4,PLANE55         !热分析单元类型
!*  
KEYOPT,4,1,0
KEYOPT,4,3,1
KEYOPT,4,4,0
KEYOPT,4,8,0
KEYOPT,4,9,0
!*  
FLST,3,2,5,ORDE,2     !热分析几何分网
FITEM,3,1   
FITEM,3,7   
AGEN,2,P51X, , , , , , ,1   
FLST,5,9,5,ORDE,2   
FITEM,5,1   
FITEM,5,-9  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   4,       0,    !热分析网格属性
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ASEL,S, , ,       9
FLST,2,1,5,ORDE,1   
FITEM,2,9   
AESIZE,P51X,5e-4,   
ASEL,S, , ,       8
FLST,2,1,5,ORDE,1   
FITEM,2,8   
AESIZE,P51X,1.5e-3,
ASEL,S, , ,       9
CM,_Y,AREA  
ASEL, , , ,       9
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,2
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ASEL,S, , ,       8
CM,_Y,AREA  
ASEL, , , ,       8
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,2
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2       !热分析几何分网完毕
!*  
ALLSEL,ALL  
/REPLOT
ALLSEL,ALL       !热辐射边界条件
/REPLOT
ASEL,S, , ,       9
ALLSEL,BELOW,AREA   
/REPLOT
/FOC,   1,  0.321526596708E-01,  0.182424536315E-01,   0.00000000000
/REPLO  
/FOC,   1,  0.282171438241E-01,  0.206529570876E-01,   0.00000000000
/REPLO  
/FOC,   1,  0.267413253816E-01,  0.225223271148E-01,   0.00000000000
/REPLO  
/FOC,   1,  0.240848521851E-01,  0.222271634263E-01,   0.00000000000
/REPLO  
/REPLOT,RESIZE  
/FOC,   1,  0.232643441263E-01,  0.195675855806E-01,   0.00000000000
/REPLO  
/FOC,   1,  0.232974817292E-01,  0.199983744178E-01,   0.00000000000
/REPLO  
FLST,5,102,1,ORDE,3
FITEM,5,7648
FITEM,5,7658
FITEM,5,-7758   
NSEL,R, , ,P51X
NPLOT   
EPLOT   
FLST,2,102,1,ORDE,3
FITEM,2,7648
FITEM,2,7658
FITEM,2,-7758   
SF,P51X,RDSF,0.68,1,   !热辐射边界条件设置完毕
ALLSEL,ALL  
/REPLOT

STEF,0.567E-07        !辐射选项控制
TOFFST,273  
RADOPT,0.10000000149,0.999999974738E-04,0,1000,0.10000000149,0.10000000149  
SPCTEMP,1,25
!*  
!*  
HEMIOPT,10            !辐射形状系数
V2DOPT,1,20,0,200   
VFOPT,OFF, , , ,BINA
!*  
SAVE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ESEL,S,MAT,,2   !设置热耦合边界
/REPLOT
FLST,2,3106,2,ORDE,6
FITEM,2,1   
FITEM,2,-1022   
FITEM,2,1623
FITEM,2,-2155   
FITEM,2,2508
FITEM,2,-4058   
BFE,P51X,FVIN, ,1,  
FINISH  
/SOL
!*  
!*  
  
!*  
ALLSEL,ALL
MFAN,ON,MFSC   !打开多物理场耦合开关
!*  
!*  
MFEL,1,1,2, , , , , , , ,    !定义电磁和热场控制文件
MFFN,1,
!*  
!*  
MFEL,2,4, , , , , , , , ,   
MFFN,2,
!*   

MFOR,       1,       2    !求解顺序
!*  
MFTI,ftime,                !求解时间控制
MFDT,tinc,0,001,1,0
MFRS,0,
!*  
!*  
MFCO,ALL,0.001,        !定义收敛误差
!*  
ANTYPE,3                 !电磁场求解类型   
HARFRQ,0,100000,         !加热频率
NSUBST, ,   
KBC,0   
!*  
!*  
OUTRES,ALL,ALL,       !输出控制
TUNIF,45,              !温度载荷
MFCM,       1, ,       !写电磁场物分析文件
MFCL,SOLU, ,           !清理分析数据
!*  
ANTYPE,4              !热分析类型
!*   
TRNOPT,FULL            
LUMPM,0
!*  
!*  
NLGEOM,0
SSTIF,0
NROPT,AUTO, ,   
EQSLV, , ,0,
PRECISION,0
MSAVE,0
TOFFST,273,            !温度偏移量设置
!*  
TUNIF,45,              !加载温度载荷
!*  
TIME,3                  !求解设置               
AUTOTS,1
DELTIM,0.01,0.005,0.01,1
KBC,1   
!*  
TSRES,ERASE

!*  
CNVTOL,HEAT,1,0.005,2, ,       !收敛准则
!*  
MFCM,       2, ,                !写温度场分析文件
MFVO,       1,       1,HGEN,       2    !耦合关系设置
MFVO,       1,       2,TEMP,       1
ALLSEL,ALL