| 查看: 1066 | 回复: 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 |
» 猜你喜欢
散金币祈福
已经有187人回复
求助,如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
物理学I论文润色/翻译怎么收费?
已经有91人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后
已经有117人回复
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
求助中医院临床中医中标的青年标书学习,谢谢
已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
【求助】结构优化不收敛
已经有9人回复
【求助】高斯计算过渡态不收敛,请问如何解决?
已经有8人回复
【求助】siesta中transiesta模块的收敛问题
已经有13人回复
【求助】castep几何优化收敛的问题
已经有11人回复
2楼2013-05-16 17:08:17










回复此楼