| 查看: 2614 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题，作者yunshan216将赠送您 5 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
yunshan216新虫 (正式写手)
|
[求助]
热分析不收敛，出现负温度已有2人参与
|
|
|
在热分析时不收敛，温度还会比常温还低，不知道怎么回事。热原模型是变半径的高斯热源加载在圆柱面上。没有考虑对流。我是哪里设置不对或是应该怎么修改呢？命令流如下： /NOPR KEYW,PR_SET,1 KEYW,PR_STRUC,0 KEYW,PR_THERM,1 KEYW,PR_FLUID,0 KEYW,PR_ELMAG,0 KEYW,MAGNOD,0 KEYW,MAGEDG,0 KEYW,MAGHFE,0 KEYW,MAGELC,0 KEYW,PR_MULTI,0 /GO /COM, Thermal /PREP7 ET,1,SOLID70 ET,2,SURF152 KEYOPT,2,1,0 KEYOPT,2,2,0 KEYOPT,2,3,0 KEYOPT,2,4,1 KEYOPT,2,5,0 KEYOPT,2,6,0 KEYOPT,2,7,0 KEYOPT,2,8,1 KEYOPT,2,9,0 KEYOPT,2,11,0 KEYOPT,2,13,0 units,si TOFFST,273 MP,KXX,1,,15 MP,DENS,1,,7930 MP,C,1,,510 MPTEMP,1,0,1424,1426,2849 MPDATA,ENTH,1,,0,5.75906e9,7.92e9,1.3675e10 *set,rj,30e-6 cylind,rj,,0,-2*rj,-90,90 wprota,,90 vsbw,all wprota,,-90 /pnum,line,1 lplot mshape,0,3d !指定划分单元的形状 mshkey,1 !指定映射网格划分 lsel,s,loc,z,0 lesize,all,,,30 lsel,s,length,,2*rj lesize,all,,,60,0.1 vmesh,all asel,s,,,12 nsla,s,1 ! /VIEW,1,1,2,3 ! /ANG,1 ! /REP,FAST TYPE, 2 MAT, 1 REAL, ESYS, 0 SECNUM, esurf,all alls FINISH /SOLU ANTYPE,4 TRNOPT,FULL LUMPM,0 ALLSEL,ALL FLST,2,5876,1,ORDE,2 FITEM,2,1 FITEM,2,-5876 IC,P51X,TEMP,25, !环境温度25度 *set,u,3 *set,id,3.2 *set,td,20 *set,F,0.6 *set,t_off,500e-6 *set,t_fre,1 !脉冲周期次数 *do,m,1,t_fre,1 *do,n,1,td,1 esel,s,type,,2 *get,nemax_2,elem,,num,max !获取单元最大序号 *get,nemin_2,elem,,num,min !获取单元最小序号 *do,i,nemin_2,nemax_2,1 !每一个单元循环 i为运算时单元号 *if,esel(i),eq,1,then !如果单元被选择则执行 获取中心点坐标 *SET,xsy,centrx(i) *SET,ysy,centry(i) *SET,zsy,centrz(i) *SET,rr,abs(sqrt((xsy/1)*(xsy/1)+(zsy)*(zsy))) *set,rt,14.823*n**0.2648*id**0.3733*1e-6 *if,ysy,gt,0,then *if,rr,le,rt,then *SET,qr,4.55*F*u*id/3.1415926/rt/rt*exp(-4.5*rr*rr/rt/rt) sfe,i,1,hflux,,qr !在i单元表面施加面载荷 *endif *endif *endif *enddo alls TIME,(m-1)*(t_off+td*1e-6)+n*1e-6 !设置时间 AUTOTS,1 DELTIM,0.5e-6, , ,1 KBC,1 TSRES,ERASE TIMINT,1 TINTP,0.005, , ,1,0.5,0.2, OUTRES,ALL,ALL, ALLSEL,ALL LNSRCH,1 ALLSEL,ALL SOLVE *enddo !脉冲结束t=td esel,s,type,,2 sfedele,all,all,hflux !sfa,3,1,conv,10,25 TIME,(m-1)*(t_off+td*1e-6)+t_off+td*1e-6 AUTOTS,1 DELTIM,100e-6, , ,1 KBC,1 TSRES,ERASE TIMINT,1 TINTP,0.005, , ,1,0.5,0.2, OUTRES,ALL,ALL, ALLSEL,ALL LNSRCH,1 ALLSEL,ALL SOLVE *enddo  @wuming524 |
回复此楼» 猜你喜欢
求助：医生想润
已经有3人回复
-
为何一方面国内大学教师无效内卷过劳逝世，而另一方面国家却在硬核科技上被卡脖子深重
已经有8人回复
-
双非院士组还是211普通老师
已经有8人回复
-
青椒的未来在哪里？
已经有4人回复
-
网传南京林业大学加拿大博士首聘期不合格，自杀
已经有31人回复
-
加上“青年”两个字，意义就变了
已经有13人回复
-
寻找高分子领域博导
已经有8人回复
-
连续3年国社科青年本子都没中，今年最后一次青年了，写本子完全浪费时间
已经有12人回复
-
有博导缺学生吗，我可以当牛做马
已经有44人回复
-
北京航空航天大学计算机学院罗洪斌课题组招收2024年学术型博士研究生
已经有7人回复
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
2楼2018-05-17 01:19:43
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
4楼2018-05-17 22:58:16
50