 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 10159 | 回復(fù): 19 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
fluent中的體熱源加載問題（實際上是熱流密度，不是熱生成率）
|
|||
|
做焊接激光模擬，采用高斯旋轉(zhuǎn)體熱源。隨時間移動。融化凝固模型。 查詢很多資料和問了若干人后，他們的建議都是加載在能量源項。 /*采用旋轉(zhuǎn)高斯體熱源*/ /*能量源項*/ #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;/*一定熱源高度對應(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后沒有錯誤，但加載出不了結(jié)果。 希望大家給點意見。 [ 來自科研家族 流體數(shù)值模擬 ] [ Last edited by b07988 on 2012-3-4 at 22:31 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 電弧 |
» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有12人回復(fù)
-
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑
已經(jīng)有14人回復(fù)
-
085410 一志愿211 22408分數(shù)359求調(diào)劑
已經(jīng)有4人回復(fù)
-
271求調(diào)劑
已經(jīng)有19人回復(fù)
-
385分 生物學(xué)（071000）求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿安徽大學(xué)計算機科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑
已經(jīng)有3人回復(fù)
-
318求調(diào)劑，計算材料方向
已經(jīng)有8人回復(fù)
-
291求調(diào)劑
已經(jīng)有25人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助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īng)有9人回復(fù)
- 關(guān)于fluent中殘差震蕩的問題~誠求各位大俠
已經(jīng)有10人回復(fù)
- fluent中壁面熱邊界條件的問題
已經(jīng)有16人回復(fù)
- fluent氣液雙相流中的Patch問題
已經(jīng)有15人回復(fù)
- 【求助】fluent 液體仿真的時候熱源怎么加載
已經(jīng)有8人回復(fù)
- 【求助】Fluent中的對稱問題
已經(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)助: 0 (幼兒園)
- 金幣: 507
- 帖子: 6
- 在線: 20.3小時
- 蟲號: 2099974
- 注冊: 2012-11-01
- 專業(yè): 金屬材料的力學(xué)行為
9樓2012-11-28 11:20:18
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
b07988(金幣+3): ★有幫助 給出了若干建議，謝謝 2012-03-05 12:30:54
xiegangmai(金幣+2): 鼓勵討論交流！ 2012-03-05 20:26:19
感謝參與，應(yīng)助指數(shù) +1
b07988(金幣+3): ★有幫助 給出了若干建議，謝謝 2012-03-05 12:30:54
xiegangmai(金幣+2): 鼓勵討論交流！ 2012-03-05 20:26:19
|
1, 我記得傳熱學(xué)一個默認的假設(shè)是熱擾動傳遞速率無限大，但對與激光加熱等加熱時間極短的問題，該假設(shè)不在適用，這樣常見的傳熱學(xué)中微分方程就要做相應(yīng)的調(diào)整，不知道你的激光焊接問題村不存在這個問題？ 2 條件語句中有幾個將y坐標(biāo)和某一值做等值比較，首先C_CENTROID取的是網(wǎng)格中心的坐標(biāo)值，是一系列離散的值，你等確定你的網(wǎng)格正好有一一個或幾個網(wǎng)格的y坐標(biāo)值是這個值（個人覺得，除非你自己輸入坐標(biāo)生成網(wǎng)格，利用網(wǎng)格生成軟件生成的網(wǎng)格是不能保證的）？也就是說你那個幾個有等號的條件判斷語句有很大的可能是遍歷完整個網(wǎng)格也找不到符合條件的，這應(yīng)該是加載不出結(jié)果的原因。所以需要修改判斷條件。 3，對于1，請教一下高斯熱源定義。 |
2樓2012-03-04 22:58:13
★ ★
xiegangmai(金幣+2): 鼓勵討論交流！ 2012-03-05 20:26:43
xiegangmai(金幣+2): 鼓勵討論交流！ 2012-03-05 20:26:43
|
1.你所說的加熱時間極短，v0=2.e-3; /*焊接速度*/，也就是熱源移動的速度，我覺得應(yīng)該不會有影響。我所寫的代碼中沒有微分方程。除了dS[eqn] =0; 2.這個問題確實沒仔細考慮，但仔細想想這樣寫的等式其實沒什么問題。熱源是一個倒錐形。所涉及的條件是剛好在那個錐尖處。網(wǎng)格比較密，產(chǎn)生的誤差可以忽略。就算關(guān)鍵的地方?jīng)]有滿足坐標(biāo)要求，但仍舊有部分熱源存在的。但實際加載沒有效果。 3.如圖，因為表達式存在1/log(H/z),其中z不能等于0，H不能等于z,所以條件語句比較復(fù)雜。呵呵 殘差圖就是直線，基本沒有波動，說明熱源沒有加載上。 熱源模型表達式 中心那根線的熱流密度（處處相等） 整體形狀 |
3樓2012-03-05 12:48:52
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★
b07988(金幣+4): ★★★很有幫助 謝謝你的建議，用message做判斷是個好方法 2012-03-05 16:09:04
xiegangmai(金幣+2): 謝謝參與！ 2012-03-05 20:27:08
b07988(金幣+4): ★★★很有幫助 謝謝你的建議，用message做判斷是個好方法 2012-03-05 16:09:04
xiegangmai(金幣+2): 謝謝參與！ 2012-03-05 20:27:08
|
我沒注意到你還有r<=rh這一句。 至于那幾個等式比較的語句，并不是你所說的網(wǎng)格較密，誤差可忽略。fluent數(shù)據(jù)是real類型的，小數(shù)點位數(shù)應(yīng)該在六位以上，因為是等值比較，只需要最后以為小數(shù)不等，條件就不成立。 你可以在你設(shè)的有熱源的條件語句部分加一個Message語句，把該點的坐標(biāo)值輸出到控制窗口里，根據(jù)這個判斷，而非根據(jù)殘差曲線。 |
4樓2012-03-05 15:30:36
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京科技，085601總分305求調(diào)劑 +9 | 半生瓜！ 2026-04-01 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-31 | 7/350 |
|
|
|
[考研] 316求調(diào)劑 +9 | 舟自梗 2026-04-01 | 10/500 |
|
|
|
[考研] 285求調(diào)劑 +7 | AZMK 2026-03-30 | 13/650 |
|
|
|
[考研] 321求調(diào)劑 +8 | y-yh 2026-04-01 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +11 | Amber00 2026-03-31 | 11/550 |
|
|
|
[碩博家園] 博一被送出聯(lián)培感覺不適應(yīng)怎么辦 +3 | 全村的狗 2026-03-31 | 3/150 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +18 | zhubinhao 2026-03-27 | 18/900 |
|
|
|
[考研] 0805求調(diào)劑 +6 | 是水分 2026-03-31 | 6/300 |
|
|
|
[考研] 0856調(diào)劑 +7 | 曲聽筠 2026-03-30 | 7/350 |
|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +6 | 小小機靈蟲 2026-03-29 | 12/600 |
|
|
|
[考研] 材料工程專碩求調(diào)劑 +10 | hyl3153942 2026-03-29 | 10/500 |
|
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|