24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yangjg1991

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 12
帖子: 19
在线: 7.9小时
虫号: 3455960
注册: 2014-10-05
性别: GG
专业: 防灾工程

[求助] 读入UDF，一初始化FLUENT就闪退已有2人参与

我现在做的是FLUENT动网格模拟，想让流体材料的粘度随着外部变量I不断变化，而I是随着时间变化的，就读入I文件，编好了UDF，在FLUENT里面能interpreted，但只要一初始化initialize，FLUENT就闪退，不知道是UDF程序有问题，还是操作步骤有问题，求大神帮助。
UDF程序如下：

/* UDF for Herschel-Bulkley viscosity */
#include "udf.h"
#include "math.h"
#include"stdio.h"
DEFINE_PROPERTY(cell_viscosity,f,t)
{
real vis, s_mag, s_mag_c, sigma_y,n,k,b;
int i;
double a[1001];
FILE* fp=fopen("I.txt","r" ) ;
for(i=0;i<1001;i++)
{
fscanf(fp,"%lf",&b);
a=b;
};
i=0;
begin_f_loop(f, t)
{
real time=CURRENT_TIME;
int nDiv=(int)(time*100);
{
if(time==0)
{
s_mag = C_STRAIN_RATE_MAG(f,t);
k= -1.574*pow(a,3)+7.102*pow(a,2)-10.24*a+4.728;
n= 0.1573*pow(a,2)-0.607*a+0.6781;
sigma_y = -2.145e4*pow(a,2)+7.505*a ;
s_mag_c = 376.1*pow(a,3)-1332*pow(a,2)+1561*a+196.8 ;
if (s_mag < s_mag_c)
{
vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c);
}
else
{
vis = sigma_y / s_mag + k*pow(s_mag/s_mag_c, (n-1));
}
}
else if(time!=0&&nDiv%2==0)
{
i++;
s_mag = C_STRAIN_RATE_MAG(f,t);
k= -1.574*pow(a,3)+7.102*pow(a,2)-10.24*a+4.728;
n= 0.1573*pow(a,2)-0.607*a+0.6781;
sigma_y = -2.145e4*pow(a,2)+7.505*a ;
s_mag_c = 376.1*pow(a,3)-1332*pow(a,2)+1561*a+196.8 ;
if (s_mag < s_mag_c)
{
vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c);
}
else
{
vis = sigma_y / s_mag + k*pow(s_mag/s_mag_c, (n-1));
}
}
else
{
s_mag = C_STRAIN_RATE_MAG(f,t);
k= -1.574*pow(a,3)+7.102*pow(a,2)-10.24*a+4.728;
n= 0.1573*pow(a,2)-0.607*a+0.6781;
sigma_y = -2.145e4*pow(a,2)+7.505*a ;
s_mag_c = 376.1*pow(a,3)-1332*pow(a,2)+1561*a+196.8 ;
if (s_mag < s_mag_c)
{
vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c);
}
else
{
vis = sigma_y / s_mag + k*pow(s_mag/s_mag_c, (n-1));
}
}
}
return vis;
}
end_f_loop(f,t);
}