| 7 | 1/1 | 返回列表 |
| 查看: 2980 | 回復(fù): 6 | |||
yunshan216新蟲 (正式寫手)
|
[求助]
熱分析不收斂,出現(xiàn)負(fù)溫度 已有2人參與
|
|
在熱分析時不收斂,溫度還會比常溫還低,不知道怎么回事。熱原模型是變半徑的高斯熱源加載在圓柱面上。沒有考慮對流。我是哪里設(shè)置不對或是應(yīng)該怎么修改呢?命令流如下: /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 !指定映射網(wǎng)格劃分 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 !脈沖周期次數(shù) *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í)行 獲取中心點坐標(biāo) *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 !設(shè)置時間 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 !脈沖結(jié)束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 |
木蟲 (著名寫手)

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

專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗: +809 |
新蟲 (正式寫手)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗: +809 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 8/400 |
|
|
[考研] 333求調(diào)劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|