| 查看: 2616 | 回复: 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 |
回复此楼» 猜你喜欢
化学B口多少分能上会呀
已经有6人回复
-
答辩
已经有4人回复
-
TsCl保护羟基为什么不反应
已经有21人回复
-
期刊投稿进度
已经有8人回复
-
为什么有的影响因子高的期刊分区不高呢？
已经有5人回复
-
又忍不住想吐槽文科的学术圈
已经有9人回复
-
大扩招
已经有14人回复
-
经常觉得挺累的
已经有7人回复
-
有没有在职教师同时做博后的？
已经有5人回复
-
两个产物点爬小板离得非常近，怎么靠大板提纯
已经有7人回复
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
4楼2018-05-17 22:58:16
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
2楼2018-05-17 01:19:43
yunshan216
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1300
- 红花: 4
- 沙发: 1
- 帖子: 513
- 在线: 34.7小时
- 虫号: 8588912
- 注册: 2018-04-22
- 专业: 机械摩擦学与表面技术
3楼2018-05-17 21:59:12
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13377.8
- 红花: 265
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
5楼2018-05-18 07:19:08
48