版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

yunshan216

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1300
紅花: 4
沙發(fā): 1
帖子: 513
在線: 34.7小時
蟲號: 8588912
注冊: 2018-04-22
專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)

[求助] 熱分析不收斂，出現(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為運(yùn)算時單元號
*if,esel(i),eq,1,then !如果單元被選擇則執(zhí)行獲取中心點(diǎn)坐標(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

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有6人回復(fù)
招08考數(shù)學(xué) 已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有6人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
285求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有7人回復(fù)
318求調(diào)劑已經(jīng)有6人回復(fù)

1樓 2018-05-16 21:18:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kingspin

木蟲 (著名寫手)

應(yīng)助: 254 (大學(xué)生)
金幣: 4472.5
紅花: 74
帖子: 1635
在線: 315.5小時
蟲號: 2585875
注冊: 2013-08-07
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

一般應(yīng)該是increment太大了。瞬態(tài)傅立葉方程對步長有要求

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

歡迎加入Digimat技術(shù)交流討論群366061054，了解復(fù)合材料多尺度仿真技術(shù)

2樓2018-05-17 01:19:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yunshan216

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1300
紅花: 4
沙發(fā): 1
帖子: 513
在線: 34.7小時
蟲號: 8588912
注冊: 2018-04-22
專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)

引用回帖:

2樓: Originally posted by kingspin at 2018-05-17 01:19:43
一般應(yīng)該是increment太大了。瞬態(tài)傅立葉方程對步長有要求

謝謝你。一個子步0.1到1微秒，應(yīng)該不是太長吧。

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

3樓2018-05-17 21:59:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kingspin

木蟲 (著名寫手)

應(yīng)助: 254 (大學(xué)生)
金幣: 4472.5
紅花: 74
帖子: 1635
在線: 315.5小時
蟲號: 2585875
注冊: 2013-08-07
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

引用回帖:

3樓: Originally posted by yunshan216 at 2018-05-17 21:59:12
謝謝你。一個子步0.1到1微秒，應(yīng)該不是太長吧。
...

我覺得太大了。因?yàn)闇囟瘸霈F(xiàn)負(fù)值那一步應(yīng)該是收斂的，所以你設(shè)最小步長很小沒用。你可以試試把最大步長降到0.01看看會不會出負(fù)值

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

歡迎加入Digimat技術(shù)交流討論群366061054，了解復(fù)合材料多尺度仿真技術(shù)

4樓2018-05-17 22:58:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shgao20

專家顧問 (職業(yè)作家)

專家經(jīng)驗(yàn): +809
仿真EPI: 1
應(yīng)助: 1311 (講師)
金幣: 13381.4
紅花: 267
帖子: 3531
在線: 473.4小時
蟲號: 3644138
注冊: 2015-01-14
專業(yè): 動力學(xué)與控制
管轄: 仿真模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

你使用了階躍方式加載，時間步長的影響可能就不大了。
改用斜坡方式加載試試，即：將所有 KBC,1 改為 KBC,0。

贊一下

回復(fù)此樓

5樓2018-05-18 07:19:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yunshan216

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1300
紅花: 4
沙發(fā): 1
帖子: 513
在線: 34.7小時
蟲號: 8588912
注冊: 2018-04-22
專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)

引用回帖:

5樓: Originally posted by shgao20 at 2018-05-18 07:19:08
你使用了階躍方式加載，時間步長的影響可能就不大了。
改用斜坡方式加載試試，即：將所有 KBC,1 改為 KBC,0。

謝謝你，最后我調(diào)節(jié)了材料焓值，就收斂了�？赡苁峭蛔兲炝�，算不出來。還有一個問題，在同一個曲面上施加熱流載荷和對流載荷，使用了表面效應(yīng)單元。但是施加了對流后計(jì)算的結(jié)果是一樣的。是哪里出現(xiàn)問題了么？

贊一下

回復(fù)此樓

6樓2018-09-15 22:17:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shgao20

專家顧問 (職業(yè)作家)

專家經(jīng)驗(yàn): +809
仿真EPI: 1
應(yīng)助: 1311 (講師)
金幣: 13381.4
紅花: 267
帖子: 3531
在線: 473.4小時
蟲號: 3644138
注冊: 2015-01-14
專業(yè): 動力學(xué)與控制
管轄: 仿真模擬

【答案】應(yīng)助回帖

多半是對流載荷有問題，不過，不了解你的具體情況，無法分析問題在哪里。

贊一下

回復(fù)此樓

7樓2018-09-16 07:13:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 yunshan216 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	16/800	2026-03-22 11:20 by 脫穎而出
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5	張vvvv 2026-03-15	7/350	2026-03-21 19:23 by ColorlessPI
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6	自由煎餅果子 2026-03-16	6/300	2026-03-21 17:21 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8	niko- 2026-03-19	9/450	2026-03-20 21:57 by luoyongfeng
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	3/150	2026-03-20 20:35 by 學(xué)員8dgXkO
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..

24小時熱門版塊排行榜

[求助] 熱分析不收斂，出現(xiàn)負(fù)溫度 已有2人參與

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] 熱分析不收斂，出現(xiàn)負(fù)溫度已有2人參與