| 1 | 1/1 | 返回列表 |
| 查看: 501 | 回復(fù): 0 | ||
[求助]
用fluent做焊接模擬。程序編譯通過了,計算不收斂,求大神看一下啊,謝謝
|
|
#include "udf.h" #define Mu 4*3.1415927e-7 #define T0 293 #define hc 80 #define dsity 1810 #define g 9.81 #define bata 1e-4 DEFINE_PROFILE(heat_flux,t,i) { real x[ND_ND]; face_t f; real r,m,n,k,qmax,a1,b1,c1,a2,b2,c2,s,e; qmax=8e6; a1=7.648e6; b1=0.5133; c1=0.5986; a2=7.647e6; b2=-0.5133; c2=0.5985; s=5.67e-8; e=0.06; begin_f_loop(f,t) { real T=F_T(f,t); F_CENTROID(x,f,t); m=x[0]; n=x[1]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(n,2)); if(r<=3) F_PROFILE(f,t,i)=a1*exp(-pow((r-b1)/c1,2))+a2*exp(-pow((r-b2)/c2,2))-hc*(T-T0)-s*e*(pow(T,4)-pow(T0,4)); else F_PROFILE(f,t,i)=-hc*(T-T0)-s*e*(pow(T,4)-pow(T0,4)); } end_f_loop(f,t) } DEFINE_SOURCE(x_velocity_source,cell,thread,dS,eqn) { real source; real x[ND_ND]; real r,m,y,z,k,Jr; begin_f_loop(cell,thread) { C_CENTROID(x,cell,thread); m=x[0]; y=x[1]; z=x[2]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(y,2)); Jr=7.473e6*exp(-pow((r-1.211)/1.451,2))+7.473e6*exp(-pow((r+1.211)/1.451,2)); source=-Mu*Jr*(1-Jr/(7.473e6))*(1-z/5)*m/(12*3.1415926*pow(r,2)); dS[eqn]=0; return source; } end_f_loop(cell,cell_thread) } DEFINE_SOURCE(y_velocity_source,cell,thread,dS,eqn) { real source; real x[ND_ND]; real r,m,y,z,k,Jr; begin_f_loop(cell,thread) { C_CENTROID(x,cell,thread); m=x[0]; y=x[1]; z=x[2]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(y,2)); Jr=7.473e6*exp(-pow((r-1.211)/1.451,2))+7.473e6*exp(-pow((r+1.211)/1.451,2)); source=-Mu*Jr*(1-Jr/(7.473e6))*(1-z/5)*y/(12*3.1415926*pow(r,2)); dS[eqn]=0; return source; } end_f_loop(cell,cell_thread) } DEFINE_SOURCE(z_velocity_source,cell,thread,dS,eqn) { real source; real T=C_T(cell,thread); real x[ND_ND]; real r,m,y,z,k,Jr; begin_f_loop(cell,thread) { C_CENTROID(x,cell,thread); m=x[0]; y=x[1]; z=x[2]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(y,2)); Jr=7.473e6*exp(-pow((r-1.211)/1.451,2))+7.473e6*exp(-pow((r+1.211)/1.451,2)); source=Mu*7473000*(1-Jr/7473000)*(1-z/5)/(12*3.1415926*pow(r,2))-dsity*g*bata*(T-800); dS[eqn]=0; return source; } end_f_loop(cell,cell_thread) } 編譯完成后,三個方向的動量源項“不是所有的控件路徑都返回值” 這是什么情況啊?而且導(dǎo)入fluent中,計算也不收斂 |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|