| 5 | 1/1 | 返回列表 |
| 查看: 2981 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
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, !環(huán)境溫度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 !每一個單元循環(huán) i為運算時單元號 *if,esel(i),eq,1,then !如果單元被選擇則執(zhí)行 獲取中心點坐標 *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 |
新蟲 (正式寫手)
木蟲 (著名寫手)

木蟲 (著名寫手)
|
我覺得太大了。因為溫度出現負值那一步應該是收斂的,所以你設最小步長很小沒用。你可以試試把最大步長降到0.01看看會不會出負值 發(fā)自小木蟲IOS客戶端 |

專家顧問 (職業(yè)作家)
![]() |
專家經驗: +809 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 384求調劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 考研化學學碩調劑,一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
[考研] 材料工程專碩 348分求調劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 307求調劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 306求0703調劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 330求調劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 310求調劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 材料專碩英一數二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 材料學求調劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 319求調劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 一志愿吉林大學材料學碩321求調劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 求調劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 288求調劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 收復試調劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 326求調劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|