热电耦合分析不收敛 求助(急!)
我在做一个含孔洞的金属薄片的热电耦合分析,然后根据热结果计算薄片里面的热应力。代码如下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+07CRITERION=2805.
CURRENT CONVERGENCE VALUE= 0.3979E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 2 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.2555E-07
HT FLOW CONVERGENCE VALUE= 0.1283E-02CRITERION= 0.1283E-05
CURRENT CONVERGENCE VALUE= 0.1009E-04CRITERION=2535. <<< CONVERGED
EQUIL ITER 3 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.2856E-10
HT FLOW CONVERGENCE VALUE= 0.6036E-05CRITERION= 0.6253E-08
CURRENT CONVERGENCE VALUE= 0.1323E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 4 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.1922E-14
HT FLOW CONVERGENCE VALUE= 0.6444E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1228E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 5 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.2032E-14
HT FLOW CONVERGENCE VALUE= 0.5674E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1239E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 6 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2199E-14
HT FLOW CONVERGENCE VALUE= 0.5299E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1215E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 7 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.2016E-14
HT FLOW CONVERGENCE VALUE= 0.5079E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1199E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 8 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2218E-14
HT FLOW CONVERGENCE VALUE= 0.5471E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1224E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER 9 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.2418E-14
HT FLOW CONVERGENCE VALUE= 0.6237E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1273E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER10 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2629E-14
HT FLOW CONVERGENCE VALUE= 0.6632E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1249E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER11 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1893E-14
HT FLOW CONVERGENCE VALUE= 0.6151E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1205E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER12 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2080E-14
HT FLOW CONVERGENCE VALUE= 0.5708E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1212E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER13 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1871E-14
HT FLOW CONVERGENCE VALUE= 0.5652E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1215E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER14 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.1995E-14
HT FLOW CONVERGENCE VALUE= 0.5325E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1210E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER15 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1861E-14
HT FLOW CONVERGENCE VALUE= 0.5683E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1195E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER16 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2032E-14
HT FLOW CONVERGENCE VALUE= 0.6178E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1206E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER17 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.2094E-14
HT FLOW CONVERGENCE VALUE= 0.6181E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1240E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER18 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2224E-14
HT FLOW CONVERGENCE VALUE= 0.6205E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1226E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER19 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1954E-14
HT FLOW CONVERGENCE VALUE= 0.5453E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1205E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER20 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.1971E-14
HT FLOW CONVERGENCE VALUE= 0.4814E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1195E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER21 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1800E-14
HT FLOW CONVERGENCE VALUE= 0.4835E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1188E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER22 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.1968E-14
HT FLOW CONVERGENCE VALUE= 0.5453E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1238E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER23 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1873E-14
HT FLOW CONVERGENCE VALUE= 0.5705E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1225E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER24 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC= -0.2136E-14
HT FLOW CONVERGENCE VALUE= 0.6120E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1237E-06CRITERION=2535. <<< CONVERGED
EQUIL ITER25 COMPLETED.NEW TRIANG MATRIX.MAX DOF INC=0.1830E-14
HT FLOW CONVERGENCE VALUE= 0.5487E-08CRITERION= 0.1644E-08
CURRENT CONVERGENCE VALUE= 0.1226E-06CRITERION=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.
请大家帮忙啊!万分感谢!
求指导啊 我看了下,没有发现计算时间,你设置下吧。尽量将载荷步设置成和热分析一样的,或者把载荷步设置小一点,deltim,0.01,0.001,0.1,试试。
希望对你有帮助 schuhi 发表于 2011-11-10 12:42 static/image/common/back.gif
我看了下,没有发现计算时间,你设置下吧。尽量将载荷步设置成和热分析一样的,或者把载荷步设置小一点,de ...
在热分析的时候就不收敛了 我先按照你的试试 谢谢了 schuhi 发表于 2011-11-10 12:42 static/image/common/back.gif
我看了下,没有发现计算时间,你设置下吧。尽量将载荷步设置成和热分析一样的,或者把载荷步设置小一点,de ...
我已经按照你给我说的方法做了一下还是不收敛 您好!那你热分析时打算设置多少时间?假如是3s,建议你这样设置试试,
time,3
deltim,0.02,0.001,0.01,
希望你早日解决问题。 schuhi 发表于 2011-11-12 09:44 static/image/common/back.gif
您好!那你热分析时打算设置多少时间?假如是3s,建议你这样设置试试,
time,3
deltim,0.02,0.001,0.01,
我的计算时间是T,一个设定好的参数 静萧寒 发表于 2011-11-14 09:25 static/image/common/back.gif
我的计算时间是T,一个设定好的参数
呵呵,不好意思啊 电流输入面上的点需要做一下耦合
页:
[1]