 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 10119 | 回復(fù): 19 | |||
[求助]
fluent中的體熱源加載問題（實(shí)際上是熱流密度，不是熱生成率）
|
|||
|
做焊接激光模擬，采用高斯旋轉(zhuǎn)體熱源。隨時(shí)間移動。融化凝固模型。 查詢很多資料和問了若干人后，他們的建議都是加載在能量源項(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;/*一定熱源高度對應(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ù)此樓» 收錄本帖的淘帖專輯推薦
 電弧 |
» 猜你喜歡
314求調(diào)劑
已經(jīng)有8人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
化工學(xué)碩306求調(diào)劑
已經(jīng)有10人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦
已經(jīng)有9人回復(fù)
-
311求調(diào)劑
已經(jīng)有12人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求助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中殘差震蕩的問題~誠求各位大俠
已經(jīng)有10人回復(fù)
- fluent中壁面熱邊界條件的問題
已經(jīng)有16人回復(fù)
- fluent氣液雙相流中的Patch問題
已經(jīng)有15人回復(fù)
- 【求助】fluent 液體仿真的時(shí)候熱源怎么加載
已經(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)助回帖
★ ★ ★ ★ ★
感謝參與，應(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è)是熱擾動傳遞速率無限大，但對與激光加熱等加熱時(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è)有等號的條件判斷語句有很大的可能是遍歷完整個(gè)網(wǎng)格也找不到符合條件的，這應(yīng)該是加載不出結(jié)果的原因。所以需要修改判斷條件。 3，對于1，請教一下高斯熱源定義。 |
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; /*焊接速度*/，也就是熱源移動的速度，我覺得應(yīng)該不會有影響。我所寫的代碼中沒有微分方程。除了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á)式 中心那根線的熱流密度（處處相等） 整體形狀 |
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í)
- 蟲號: 1529733
- 注冊: 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í)
- 蟲號: 2099974
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +6 | Liwangman 2026-03-15 | 6/300 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|