24小时热门版块排行榜

返回列表

嗜睡猪崽子

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 558.4
散金: 105
红花: 9
帖子: 161
在线: 65.2小时
虫号: 1831911
注册: 2012-05-23
性别: GG
专业: 高分子材料的加工与成型

[求助] 用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)

} 编译完成后，三个方向的动量源项“不是所有的控件路径都返回值” 这是什么情况啊？而且导入fluent中，计算也不收敛

回复此楼