| 查看: 2593 | 回复: 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 |
回复此楼» 猜你喜欢
美国教授开除中国研究生：我就不该录取你！
已经有6人回复
-
院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度
已经有6人回复
-
某系主任被同校教授举报包养情妇、长期嫖娼
已经有14人回复
-
一审一个审稿人，小修，会怎么样呀？
已经有5人回复
-
特别资助审核状态
已经有6人回复
-
国自然青年基金只有一篇能第一标注，其他都是第二标注能结题吗？
已经有5人回复
-
杭州国企和浙江高校如何选择？
已经有20人回复
-
发个文章，结果被某学报请的外审PUA了
已经有8人回复
-
估计今年青基又没戏
已经有8人回复
-
267求调剂
已经有3人回复
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13377.8
- 红花: 265
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
5楼2018-05-18 07:19:08
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13377.8
- 红花: 265
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
7楼2018-09-16 07:13:07
25