 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 10120 | 回復: 19 | |||
[求助]
fluent中的體熱源加載問題（實際上是熱流密度，不是熱生成率）
|
|||
|
做焊接激光模擬，采用高斯旋轉體熱源。隨時間移動。融化凝固模型。 查詢很多資料和問了若干人后，他們的建議都是加載在能量源項。 /*采用旋轉高斯體熱源*/ /*能量源項*/ #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.)); /*坐標與熱源中心的距離*/ rh=sqrt(log(HH/x[2])/3.)*oq;/*一定熱源高度對應的半徑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; } 條件語句比較復雜，主要怕出問題。呵呵 interpret后沒有錯誤，但加載出不了結果。 希望大家給點意見。 [ 來自科研家族 流體數(shù)值模擬 ] [ Last edited by b07988 on 2012-3-4 at 22:31 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 電弧 |
» 猜你喜歡
311求調(diào)劑
已經(jīng)有4人回復
-
304求調(diào)劑
已經(jīng)有4人回復
-
354求調(diào)劑
已經(jīng)有7人回復
-
321求調(diào)劑
已經(jīng)有3人回復
-
274求調(diào)劑
已經(jīng)有5人回復
-
281求調(diào)劑（0805）
已經(jīng)有15人回復
-
085600材料與化工調(diào)劑 324分
已經(jīng)有8人回復
-
能源材料化學課題組招收碩士研究生8-10名
已經(jīng)有10人回復
-
267一志愿南京工業(yè)大學0817化工求調(diào)劑
已經(jīng)有8人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有10人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求助fluent中速度顯示問題
已經(jīng)有22人回復
- fluent中關于拉瓦爾噴管的問題
已經(jīng)有26人回復
- 求助：fluent中關于壁面?zhèn)鳠釂栴}
已經(jīng)有12人回復
- fluent中的數(shù)據(jù)統(tǒng)計問題
已經(jīng)有9人回復
- 關于fluent中殘差震蕩的問題~誠求各位大俠
已經(jīng)有10人回復
- fluent中壁面熱邊界條件的問題
已經(jīng)有16人回復
- fluent氣液雙相流中的Patch問題
已經(jīng)有15人回復
- 【求助】fluent 液體仿真的時候熱源怎么加載
已經(jīng)有8人回復
- 【求助】Fluent中的對稱問題
已經(jīng)有10人回復
- 【求助】fluent中固體與流體區(qū)域間的interior界面，出現(xiàn)警告提示該怎么解決
已經(jīng)有6人回復
- 【求助】FLUENT三維模型中入口速度的方向設置問題【已完成】
已經(jīng)有17人回復
- 【求助】FLUENT中的內(nèi)熱源怎么設置？在哪里設置？
已經(jīng)有16人回復
- 【求助】有關fluent中DPM模型的模擬問題
已經(jīng)有16人回復
【答案】應助回帖
★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
b07988(金幣+3): ★有幫助 給出了若干建議，謝謝 2012-03-05 12:30:54
xiegangmai(金幣+2): 鼓勵討論交流！ 2012-03-05 20:26:19
感謝參與，應助指數(shù) +1
b07988(金幣+3): ★有幫助 給出了若干建議，謝謝 2012-03-05 12:30:54
xiegangmai(金幣+2): 鼓勵討論交流！ 2012-03-05 20:26:19
|
1, 我記得傳熱學一個默認的假設是熱擾動傳遞速率無限大，但對與激光加熱等加熱時間極短的問題，該假設不在適用，這樣常見的傳熱學中微分方程就要做相應的調(diào)整，不知道你的激光焊接問題村不存在這個問題？ 2 條件語句中有幾個將y坐標和某一值做等值比較，首先C_CENTROID取的是網(wǎng)格中心的坐標值，是一系列離散的值，你等確定你的網(wǎng)格正好有一一個或幾個網(wǎng)格的y坐標值是這個值（個人覺得，除非你自己輸入坐標生成網(wǎng)格，利用網(wǎng)格生成軟件生成的網(wǎng)格是不能保證的）？也就是說你那個幾個有等號的條件判斷語句有很大的可能是遍歷完整個網(wǎng)格也找不到符合條件的，這應該是加載不出結果的原因。所以需要修改判斷條件。 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; /*焊接速度*/，也就是熱源移動的速度，我覺得應該不會有影響。我所寫的代碼中沒有微分方程。除了dS[eqn] =0; 2.這個問題確實沒仔細考慮，但仔細想想這樣寫的等式其實沒什么問題。熱源是一個倒錐形。所涉及的條件是剛好在那個錐尖處。網(wǎng)格比較密，產(chǎn)生的誤差可以忽略。就算關鍵的地方?jīng)]有滿足坐標要求，但仍舊有部分熱源存在的。但實際加載沒有效果。 3.如圖，因為表達式存在1/log(H/z),其中z不能等于0，H不能等于z,所以條件語句比較復雜。呵呵 殘差圖就是直線，基本沒有波動，說明熱源沒有加載上。 熱源模型表達式 中心那根線的熱流密度（處處相等） 整體形狀 |
3樓2012-03-05 12:48:52
6樓2012-03-05 18:09:05
【答案】應助回帖
★ ★ ★ ★ ★ ★
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ù)應該在六位以上，因為是等值比較，只需要最后以為小數(shù)不等，條件就不成立。 你可以在你設的有熱源的條件語句部分加一個Message語句，把該點的坐標值輸出到控制窗口里，根據(jù)這個判斷，而非根據(jù)殘差曲線。 |
4樓2012-03-05 15:30:36
5樓2012-03-05 16:19:29
nuliwo1
銀蟲 (小有名氣)
- 應助: 3 (幼兒園)
- 金幣: 4176.8
- 紅花: 1
- 帖子: 57
- 在線: 35小時
- 蟲號: 1529733
- 注冊: 2011-12-08
- 專業(yè): 金屬基復合材料
7樓2012-03-06 20:36:41
8樓2012-03-07 12:15:11
9樓2012-11-28 11:20:18
10樓2012-11-30 11:23:27
相關版塊跳轉
我要訂閱樓主 b07988 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調(diào)劑 +7 | Liwangman 2026-03-15 | 7/350 |
|
|
|
[考研] 0703化學調(diào)劑，求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 070300化學319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考博] 環(huán)境領域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學666 2026-03-16 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
|
[考研] 本人考085602 化學工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 285化工學碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|