热电耦合分析不收敛 求助（急！）

静萧寒

我在做一个含孔洞的金属薄片的热电耦合分析，然后根据热结果计算薄片里面的热应力。代码如下

FINISH
/CLEAR
/filname,DIANREYINGLI            
/title,ElectricThermalstress   
*DIM,A,,1,1,1
*DIM,TBUCHANG,,1,1,1   
*DIM,V,,1,4                                
*VREAD,A(1),NG,TXT,,JIK,1,1,1
(F5.0)
*VREAD,TBUCHANG(1),TBUCHANG,TXT,,JIK,1,1,1
(F20.6)
*VREAD,V(1,1),V,TXT,,IJK,1,4,1
(4F12.6)
                     
NGA=A(1)                     
VA=V(1,1)
L=V(1,2)
H=V(1,3)*V(1,4)
T=TBUCHANG(1)

/PREP7
*DIM,B,,NGA,2                               !定义二维数组NGA*2
*VREAD,B(1,1),FORANSYS,TXT,,JIK,2,NGA       !从FORANSYS.TXT读入数据给二维数组
(2F14.6)  ！以上都是读取文本文件的数据
*DO,J,1,NGA  
   K,J,B(J,1),B(J,2),0
*ENDDO                                       !创建2*NG-2个关键点
*DO,J,1,NGA-1
LSTR,J,J+1                                  !依据编号依次把创建关键点之间的连线
*ENDDO
LSTR,NGA,1
AL,ALL                                       !形成裂纹面
BLC5,0,0,L,H                                 !形成矩形面
ASBA,2,1                                     !相减形成模型面
ET,1,PLANE67                                !定义单元类型热电耦合
MP,DENS,1,26989                             !定义材料属性定值密度
MP,KXX,1,210                                !热传导率 kxx210
MP,C,1,880                                  !比热容
MP,RSVX,1,,2.65E-8
SMRT,1         
AMESH,ALL                                    !划分网格
DL,NGA+4,,VOLT,VA                            !边界条件
DL,NGA+2,,VOLT,0

/SOLU
TOFFST,273
ANTYPE,4
TUNIF,25,   
TIME,T   
AUTOTS,ON   
SOLVE
FINISH
/POST1                                      !后处理数据的提取
*DIM,TMP,,NGA,1
*DO,J,1,NGA                                 !关键点所对应的节点编号
*GET,TMP(J),NODE,J,TEMP                     !和关键点编号相同故为1-NGA
*ENDDO
     
*CFOPEN,WENDU,TXT           
*VWRITE,TMP(1)
(F10.6)
*CFCLOS
/prep7                                      !重新进入前处理，进行热应力耦合分析
ETCH,TTS                                    !转化单元类型
MP,EX,1,7E10                                !输入材料结构性能参数
MP,ALPX,1,23E-6
MP,NUXY,1,0.26

LDREAD,TEMP,,, , ,'DIANREYINGLI','rst',' '   
/SOLU
SOLVE
/POST1
*DIM,U,,NGA,3
*DO,J,1,NGA                                  !关键点所对应的节点编号
  *GET,U(J,1),NODE,J,S,1                     !和关键点编号相同故为1~NGA
  *GET,U(J,2),NODE,J,S,2                       
  *GET,U(J,3),NODE,J,S,3
*ENDDO     
*CFOPEN,TJIEGUO,TXT
*VWRITE,U(1,1),U(1,2),U(1,3)                 !将主应力写入文本TJIEGUO.TXT中
(3F17.5)
*CFCLOS
结果出现如下错误
EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   12.00   
     HT FLOW CONVERGENCE VALUE= 0.2805E+07  CRITERION=  2805.   
     CURRENT CONVERGENCE VALUE= 0.3979E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   2 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2555E-07
     HT FLOW CONVERGENCE VALUE= 0.1283E-02  CRITERION= 0.1283E-05
     CURRENT CONVERGENCE VALUE= 0.1009E-04  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   3 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2856E-10
     HT FLOW CONVERGENCE VALUE= 0.6036E-05  CRITERION= 0.6253E-08
     CURRENT CONVERGENCE VALUE= 0.1323E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   4 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.1922E-14
     HT FLOW CONVERGENCE VALUE= 0.6444E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1228E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   5 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2032E-14
     HT FLOW CONVERGENCE VALUE= 0.5674E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1239E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   6 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2199E-14
     HT FLOW CONVERGENCE VALUE= 0.5299E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1215E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   7 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2016E-14
     HT FLOW CONVERGENCE VALUE= 0.5079E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1199E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   8 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2218E-14
     HT FLOW CONVERGENCE VALUE= 0.5471E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1224E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER   9 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2418E-14
     HT FLOW CONVERGENCE VALUE= 0.6237E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1273E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  10 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2629E-14
     HT FLOW CONVERGENCE VALUE= 0.6632E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1249E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  11 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1893E-14
     HT FLOW CONVERGENCE VALUE= 0.6151E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1205E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  12 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2080E-14
     HT FLOW CONVERGENCE VALUE= 0.5708E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1212E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  13 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1871E-14
     HT FLOW CONVERGENCE VALUE= 0.5652E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1215E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  14 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.1995E-14
     HT FLOW CONVERGENCE VALUE= 0.5325E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1210E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  15 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1861E-14
     HT FLOW CONVERGENCE VALUE= 0.5683E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1195E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  16 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2032E-14
     HT FLOW CONVERGENCE VALUE= 0.6178E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1206E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  17 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2094E-14
     HT FLOW CONVERGENCE VALUE= 0.6181E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1240E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  18 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2224E-14
     HT FLOW CONVERGENCE VALUE= 0.6205E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1226E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  19 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1954E-14
     HT FLOW CONVERGENCE VALUE= 0.5453E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1205E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  20 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.1971E-14
     HT FLOW CONVERGENCE VALUE= 0.4814E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1195E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  21 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1800E-14
     HT FLOW CONVERGENCE VALUE= 0.4835E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1188E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  22 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.1968E-14
     HT FLOW CONVERGENCE VALUE= 0.5453E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1238E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  23 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1873E-14
     HT FLOW CONVERGENCE VALUE= 0.5705E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1225E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  24 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2136E-14
     HT FLOW CONVERGENCE VALUE= 0.6120E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1237E-06  CRITERION=  2535.     <<< CONVERGED
    EQUIL ITER  25 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1830E-14
     HT FLOW CONVERGENCE VALUE= 0.5487E-08  CRITERION= 0.1644E-08
     CURRENT CONVERGENCE VALUE= 0.1226E-06  CRITERION=  2535.     <<< CONVERGED
** ERROR **                           CP =      10.891   TIME= 14:20:02
Solution not converged at time 1.E-06 (load step 1 substep 1).         
  Run terminated.                                                        
** WARNING **                         CP =      10.891   TIME= 14:20:02
The unconverged solution (identified as time 1.E-06 substep 999999) is  
output for analysis debug purposes.  Results should not be used for     
any other purpose.                                 

请大家帮忙啊！万分感谢！

静萧寒

求指导啊

schuhi

我看了下，没有发现计算时间，你设置下吧。尽量将载荷步设置成和热分析一样的，或者把载荷步设置小一点，deltim,0.01,0.001,0.1,试试。
希望对你有帮助

静萧寒

schuhi 发表于 2011-11-10 12:42
我看了下，没有发现计算时间，你设置下吧。尽量将载荷步设置成和热分析一样的，或者把载荷步设置小一点，de ...

在热分析的时候就不收敛了 我先按照你的试试 谢谢了

静萧寒

schuhi 发表于 2011-11-10 12:42
我看了下，没有发现计算时间，你设置下吧。尽量将载荷步设置成和热分析一样的，或者把载荷步设置小一点，de ...

我已经按照你给我说的方法做了一下  还是不收敛

schuhi

您好！那你热分析时打算设置多少时间？假如是3s，建议你这样设置试试，
time,3
deltim,0.02,0.001,0.01,
希望你早日解决问题。

静萧寒

schuhi 发表于 2011-11-12 09:44
您好！那你热分析时打算设置多少时间？假如是3s，建议你这样设置试试，
time,3
deltim,0.02,0.001,0.01,

我的计算时间是T，一个设定好的参数

schuhi

静萧寒 发表于 2011-11-14 09:25
我的计算时间是T，一个设定好的参数

呵呵，不好意思啊

xhj543240343

电流输入面上的点需要做一下耦合