UDF非牛顿流体粘度

hanhasiaoqier

下面是我从网上找的模版，然后自己进行了修改，这个UDFinterprete成功了，但是initialize出错了，错误提示是：FLUENT received fatal signal (ACCESS_VIOLATION) 1. Note exact events leading to error. 2. Save case/data under new name. 3. Exit program and restart。。。4.。。

我加载的是粘度UDF ，是从网上找的模版，我自己都没吃透中间四五行程序语句，我按我的粘度模型进行了修改，大家帮我看看这个程序哪有问题可以吗？

#include<stdio.h>
#include "udf.h"
real T,vis, s_mag, s_mag_c, sigma_y,n,k,ky;
real C_1=0.06;
real C_2=1.5;
real C_3=41580;
real C_4=366;
real C_5=6;
int ia ;
DEFINE_PROPERTY(hb_viscosity,c,t)
{
T=C_T(c,t);
s_mag=C_STRAIN_RATE_MAG(c,t);

if (ia==0.0)
{C_1=RP_Get_Real("c_1");
C_2=RP_Get_Real("c_2");
C_3=RP_Get_Real("c_3");
C_4=RP_Get_Real("c_4");
C_5=RP_Get_Real("c_5");
ia=1;}
k=C_1;
n=C_2;
sigma_y=C_3;
s_mag_c=C_4;
ky=C_5;
if(s_mag<s_mag_c)
{vis=ky*pow(s_mag,(n-1));}
else
{vis=sigma_y/s_mag+k*pow(s_mag,(n-1));}
return vis;
}
网上有些人说是要定义内存空间什么的，但是也有人说定义完内存空间，初始是可以通过，但是计算还是会出错，
大家帮帮我吧，，小弟等着毕业呢啊，，谢谢各位大虾啊

needload

同求











上海搬家公司