 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 10118 | 回復(fù): 19 | |||
[求助]
fluent中的體熱源加載問題（實(shí)際上是熱流密度，不是熱生成率）
|
|||
|
做焊接激光模擬，采用高斯旋轉(zhuǎn)體熱源。隨時(shí)間移動(dòng)。融化凝固模型。 查詢很多資料和問了若干人后，他們的建議都是加載在能量源項(xiàng)。 /*采用旋轉(zhuǎn)高斯體熱源*/ /*能量源項(xiàng)*/ #include "udf.h" DEFINE_SOURCE(heat_flux, cell, thread, dS, eqn) { real x[ND_ND]; real y,I,U,oq,t0,v0,Q,time,cs,HH; real r,rh; real source_heat; cell_t c; C_CENTROID(x,c,thread); y=0.7; /*有效功率因子*/ U=16.; I=100.; v0=2.e-3; /*焊接速度*/ Q=y*I*U; /*有效功率*/ cs=3./pow(oq,2.); /*形狀因子*/ HH=1.e-3; /*熱源高度*/ oq=2.5e-3; /*熱源半徑*/ time= RP_Get_Real("flow-time"  ;r=sqrt(pow(x[0]-v0*time,2.)+pow(x[1],2.)); /*坐標(biāo)與熱源中心的距離*/ rh=sqrt(log(HH/x[2])/3.)*oq;/*一定熱源高度對(duì)應(yīng)的半徑rh*/ if (x[2]==0&&r<=oq) { source_heat=3.*cs*Q/(3.14*HH*(1.-1./exp(3)))*exp(-3.*pow(r,2.)/pow(oq,2.)); dS[eqn] =0; } else if (x[2]==HH&&r==0.) { source_heat=3.*cs*Q/(3.14*HH*(1.-1./exp(3))); dS[eqn] =0; } else if (x[2]>HH) { source_heat=0; dS[eqn] =0; } else if (r<=rh) { source_heat=3.*cs*Q/(3.14*HH*(1.-1./exp(3)))*exp(-3.*cs/log10(HH/x[2])*pow(r,2.)); dS[eqn] =0; } else { source_heat= dS[eqn] =0; } return source_heat; } 條件語句比較復(fù)雜，主要怕出問題。呵呵 interpret后沒有錯(cuò)誤，但加載出不了結(jié)果。 希望大家給點(diǎn)意見。 [ 來自科研家族 流體數(shù)值模擬 ] [ Last edited by b07988 on 2012-3-4 at 22:31 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 電弧 |
» 猜你喜歡
280求調(diào)劑
已經(jīng)有7人回復(fù)
-
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有5人回復(fù)
-
303求調(diào)劑
已經(jīng)有6人回復(fù)
-
290求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703化學(xué)求調(diào)劑 總分331
已經(jīng)有3人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有8人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
-
材料，紡織，生物（0856、0710），化學(xué)招生啦
已經(jīng)有9人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助fluent中速度顯示問題
已經(jīng)有22人回復(fù)
- fluent中關(guān)于拉瓦爾噴管的問題
已經(jīng)有26人回復(fù)
- 求助：fluent中關(guān)于壁面?zhèn)鳠釂栴}
已經(jīng)有12人回復(fù)
- fluent中的數(shù)據(jù)統(tǒng)計(jì)問題
已經(jīng)有9人回復(fù)
- 關(guān)于fluent中殘差震蕩的問題~誠(chéng)求各位大俠
已經(jīng)有10人回復(fù)
- fluent中壁面熱邊界條件的問題
已經(jīng)有16人回復(fù)
- fluent氣液雙相流中的Patch問題
已經(jīng)有15人回復(fù)
- 【求助】fluent 液體仿真的時(shí)候熱源怎么加載
已經(jīng)有8人回復(fù)
- 【求助】Fluent中的對(duì)稱問題
已經(jīng)有10人回復(fù)
- 【求助】fluent中固體與流體區(qū)域間的interior界面，出現(xiàn)警告提示該怎么解決
已經(jīng)有6人回復(fù)
- 【求助】FLUENT三維模型中入口速度的方向設(shè)置問題【已完成】
已經(jīng)有17人回復(fù)
- 【求助】FLUENT中的內(nèi)熱源怎么設(shè)置？在哪里設(shè)置？
已經(jīng)有16人回復(fù)
- 【求助】有關(guān)fluent中DPM模型的模擬問題
已經(jīng)有16人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
b07988(金幣+3): ★有幫助 給出了若干建議，謝謝 2012-03-05 12:30:54
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2012-03-05 20:26:19
感謝參與，應(yīng)助指數(shù) +1
b07988(金幣+3): ★有幫助 給出了若干建議，謝謝 2012-03-05 12:30:54
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2012-03-05 20:26:19
|
1, 我記得傳熱學(xué)一個(gè)默認(rèn)的假設(shè)是熱擾動(dòng)傳遞速率無限大，但對(duì)與激光加熱等加熱時(shí)間極短的問題，該假設(shè)不在適用，這樣常見的傳熱學(xué)中微分方程就要做相應(yīng)的調(diào)整，不知道你的激光焊接問題村不存在這個(gè)問題？ 2 條件語句中有幾個(gè)將y坐標(biāo)和某一值做等值比較，首先C_CENTROID取的是網(wǎng)格中心的坐標(biāo)值，是一系列離散的值，你等確定你的網(wǎng)格正好有一一個(gè)或幾個(gè)網(wǎng)格的y坐標(biāo)值是這個(gè)值（個(gè)人覺得，除非你自己輸入坐標(biāo)生成網(wǎng)格，利用網(wǎng)格生成軟件生成的網(wǎng)格是不能保證的）？也就是說你那個(gè)幾個(gè)有等號(hào)的條件判斷語句有很大的可能是遍歷完整個(gè)網(wǎng)格也找不到符合條件的，這應(yīng)該是加載不出結(jié)果的原因。所以需要修改判斷條件。 3，對(duì)于1，請(qǐng)教一下高斯熱源定義。 |
2樓2012-03-04 22:58:13
★ ★
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2012-03-05 20:26:43
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2012-03-05 20:26:43
|
1.你所說的加熱時(shí)間極短，v0=2.e-3; /*焊接速度*/，也就是熱源移動(dòng)的速度，我覺得應(yīng)該不會(huì)有影響。我所寫的代碼中沒有微分方程。除了dS[eqn] =0; 2.這個(gè)問題確實(shí)沒仔細(xì)考慮，但仔細(xì)想想這樣寫的等式其實(shí)沒什么問題。熱源是一個(gè)倒錐形。所涉及的條件是剛好在那個(gè)錐尖處。網(wǎng)格比較密，產(chǎn)生的誤差可以忽略。就算關(guān)鍵的地方?jīng)]有滿足坐標(biāo)要求，但仍舊有部分熱源存在的。但實(shí)際加載沒有效果。 3.如圖，因?yàn)楸磉_(dá)式存在1/log(H/z),其中z不能等于0，H不能等于z,所以條件語句比較復(fù)雜。呵呵 殘差圖就是直線，基本沒有波動(dòng)，說明熱源沒有加載上。 熱源模型表達(dá)式 中心那根線的熱流密度（處處相等） 整體形狀 |
3樓2012-03-05 12:48:52
6樓2012-03-05 18:09:05
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★
b07988(金幣+4): ★★★很有幫助 謝謝你的建議，用message做判斷是個(gè)好方法 2012-03-05 16:09:04
xiegangmai(金幣+2): 謝謝參與！ 2012-03-05 20:27:08
b07988(金幣+4): ★★★很有幫助 謝謝你的建議，用message做判斷是個(gè)好方法 2012-03-05 16:09:04
xiegangmai(金幣+2): 謝謝參與！ 2012-03-05 20:27:08
|
我沒注意到你還有r<=rh這一句。 至于那幾個(gè)等式比較的語句，并不是你所說的網(wǎng)格較密，誤差可忽略。fluent數(shù)據(jù)是real類型的，小數(shù)點(diǎn)位數(shù)應(yīng)該在六位以上，因?yàn)槭堑戎当容^，只需要最后以為小數(shù)不等，條件就不成立。 你可以在你設(shè)的有熱源的條件語句部分加一個(gè)Message語句，把該點(diǎn)的坐標(biāo)值輸出到控制窗口里，根據(jù)這個(gè)判斷，而非根據(jù)殘差曲線。 |
4樓2012-03-05 15:30:36
★ ★
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2012-03-05 20:27:27
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2012-03-05 20:27:27
|
在HH高度處，如果此處沒有節(jié)點(diǎn)的話，條件不成立其實(shí)沒關(guān)系的。。呵呵 附近的節(jié)點(diǎn)總有一個(gè)熱流密度賦值給它。 問題是： 能量源項(xiàng)的單位是：w/m3 而這個(gè)熱源是w/m2 體熱源的加載都是用能量源項(xiàng)表示，我就疑惑了~。 但文獻(xiàn)表示這個(gè)熱源就是熱流密度。。。 用message是個(gè)好方法，接觸udf幾個(gè)月， 沒有調(diào)試經(jīng)驗(yàn)，謝謝你的建議。 |
5樓2012-03-05 16:19:29
nuliwo1
銀蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 4176.8
- 紅花: 1
- 帖子: 57
- 在線: 35小時(shí)
- 蟲號(hào): 1529733
- 注冊(cè): 2011-12-08
- 專業(yè): 金屬基復(fù)合材料
7樓2012-03-06 20:36:41
8樓2012-03-07 12:15:11
科研小王子
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 507
- 帖子: 6
- 在線: 20.3小時(shí)
- 蟲號(hào): 2099974
- 注冊(cè): 2012-11-01
- 專業(yè): 金屬材料的力學(xué)行為
9樓2012-11-28 11:20:18
10樓2012-11-30 11:23:27
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 b07988 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 求調(diào)劑（材料與化工327） +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|