| 查看: 1054 | 回复: 1 | |||
[求助]
求助各位,出现不收敛的问题
|
|
在沿着圆柱轴向的一个面施加热,但是出现不收敛的问题,实在不知道是哪里出了问题,大家帮帮忙看一下~ 警告提示: *** WARNING *** CP = 770.083 TIME= 21:42:28 The calculated reference HT FLOW CONVERGENCE VALUE= 1.557378551E-14 is less than the default minimum of 1E-6 (which can be modified using the MINREF field of the CNVTOL command). Check results carefully. 程序命令流如下: FINISH /CLEAR /BATCH /FILNAME,AL-THERMAL,1 /TITLE,THERMAL R0=0.011 !管道外径 Ri=0.01 !管道内径 D=0.001 !管厚 L=0.003 !管长 /PREP7 ET,1,SOLID70 MPTEMP,1,20,100,200,250,300,350 !temp table(`c) 6个点 !材料密度 MP,DENS,1,2.8e3 !density(Kg/m^3) [2800Kg/m^3] !比热 c = 770.045 + 1.671*t - 0.000152*t^2 MPDATA,c,1,1,0.900e3,0.921e3,1.047e3,1.092e3,1.130e3,1.172e3 !热传导系数 MPDATA,KXX,1,1,0.121e3,0.121e3,0.121e3,0.121e3,0.121e3,0.121e3 TOFFST,273 !参考温度,初始温度 PCIRC,R0,RI,-90,90 VOFFST,1,-L,, WPROTA,0,0,90 WPROTA,0,88,0 VSBW,1 WPROTA,0,4,0 VSBW,3 WPROTA,0,-92,0 WPROTA,0,0,-90 NUMMRG,KP NUMCMP,ALL SAVE LESIZE,2,,,25 !径向划分 LESIZE,4,,,25 LESIZE,5,,,25 LESIZE,15,,,25 LESIZE,6,,,25 LESIZE,8,,,25 LESIZE,18,,,25 LESIZE,22,,,25 LESIZE,16,,,116 !其他区域网格划分 LESIZE,17,,,116 LESIZE,7,,,116 LESIZE,21,,,116 LESIZE,1,,,116 LESIZE,3,,,116 LESIZE,13,,,116 LESIZE,14,,,116 LESIZE,25,,,6 !施热区的网格划分 LESIZE,26,,,6 LESIZE,27,,,6 LESIZE,28,,,6 MSHAPE,0,3D LESIZE,23,,,15 VSWEEP,1,1,2 LESIZE,19,,,15 VSWEEP,2,7,8 VSWEEP,3,13,14 ALLSEL,ALL SAVE FINISH /SOLU ANTYPE,4 NROPT,FULL NLGEOM,OFF !关闭大变形 LUMPM,0 TIMINT,ON !打开瞬态效应 OUTPR,BASIC,ALL, OUTRES,ALL,ALL TINTP,,0.25,0.5,1,0.5,0 toffst,0 TUNIF,30 tref,30 asel,s,area,,15 esla *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNCCSYS *SET,_FNCNAME,'load' *SET,_FNCCSYS,1 ! /INPUT,myload.func,,,1 *DIM,%_FNCNAME%,TABLE,6,35,1,,,,%_FNCCSYS% ! ! Begin of equation: 10^20/5.316*{TIME}*exp(-(0.011*{Y})^2/0.003^2)*exp(-{TIME} ! /10^-9) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), 0.0 *SET,%_FNCNAME%(4,0,1), 0.0 *SET,%_FNCNAME%(5,0,1), 0.0 *SET,%_FNCNAME%(6,0,1), 0.0 *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 20, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 1, -3, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 5.316, 0, 0, -3 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 4, -1 *SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 1 *SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 0.011, 0, 0, 3 *SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -2, 3, 3 *SET,%_FNCNAME%(0,12,1), 0.0, -2, 0, 2, 0, 0, -3 *SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -3, 17, -2 *SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -4, 3, -5 *SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 0.003, 0, 0, 0 *SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 2, 0, 0, -3 *SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, -3, 17, -4 *SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -2, 4, -5 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 7, 1, -3, 0, 0 *SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -4, 3, 1 *SET,%_FNCNAME%(0,25,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,26,1), 0.0, -4, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -2, 2, -4 *SET,%_FNCNAME%(0,28,1), 0.0, -2, 0, 9, 0, 0, -5 *SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, -5, 3, -2 *SET,%_FNCNAME%(0,30,1), 0.0, -2, 0, 10, 0, 0, -4 *SET,%_FNCNAME%(0,31,1), 0.0, -5, 0, 1, -2, 17, -4 *SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 4, -5 *SET,%_FNCNAME%(0,33,1), 0.0, -1, 7, 1, -2, 0, 0 *SET,%_FNCNAME%(0,34,1), 0.0, -2, 0, 1, -3, 3, -1 *SET,%_FNCNAME%(0,35,1), 0.0, 99, 0, 1, -2, 0, 0 ! End of equation: 10^20/5.316*{TIME}*exp(-(0.011*{Y})^2/0.003^2)*exp(-{TIME}/ ! 10^-9) !--> SFE,ALL,1,HFLUX,%LOAD% ALLSEL,ALL SBCTRAN TRNOPT,FULL !设置计算方法 *DO,DELTA,1e-9,10e-9,1e-9 !设置循环计算 TIME,DELTA !设置荷载步时间 AUTOTS,0 !关闭自动时间步长 KBC,0 !设置阶梯荷载步 TSRES,ERASE !Defines an array of keytimes at which the time-stepping strategy changes OUTRES,ALL,ALL !输出所有项 SOLVE !开始计算 *ENDDO !终止循环计算 FINISH |
» 猜你喜欢
基金申请
已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有135人回复
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
帮我的英语口语老师找学生
已经有0人回复
什么时候开奖?
已经有13人回复
散金币祈福
已经有94人回复
青基已中
已经有1人回复
散金币祈福
已经有185人回复
求助,如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
【求助】结构优化不收敛
已经有9人回复
【求助】高斯计算过渡态不收敛,请问如何解决?
已经有8人回复
【求助】siesta中transiesta模块的收敛问题
已经有13人回复
【求助】castep几何优化收敛的问题
已经有11人回复
2楼2013-05-16 17:08:17










回复此楼
30