应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 读入UDF，一初始化FLUENT就闪退
111/2返回列表12下一页
查看: 3649  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yangjg1991

铁虫 (初入文坛)

[求助] 读入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);
}
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2016-04-22 14:39:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yangjg1991

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by moons2 at 2016-04-22 20:38:20
不加载udf的情况下初始化会闪退吗

不加载UDF的话，一切运行正常
赞一下(1人) 回复此楼
5楼2016-04-27 16:46:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangjg1991

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by ttzhsz at 2016-04-22 19:12:21
是不是运行库没装好

运行其他不随时间变化的UDF是没有问题的，所以应该不是运行库的问题
赞一下(1人) 回复此楼
7楼2016-04-27 16:50:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ttzhsz

铁虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
是不是运行库没装好
赞一下 回复此楼
2楼2016-04-22 19:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwx1990

金虫 (正式写手)
问题很多，当然不能成功了。比如，除了包含udf.h以外都不要，文本l中读数据错…

发自小木虫Android客户端
赞一下 回复此楼
努力工作
3楼2016-04-22 20:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moons2

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
不加载udf的情况下初始化会闪退吗
赞一下 回复此楼
4楼2016-04-22 20:38:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangjg1991

铁虫 (初入文坛)
引用回帖:
3楼: Originally posted by xwx1990 at 2016-04-22 20:13:22
问题很多，当然不能成功了。比如，除了包含udf.h以外都不要，文本l中读数据错…

你好，请问具体有什么问题，烦请一一列举一下，多谢多谢
赞一下 回复此楼
6楼2016-04-27 16:47:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangjg1991

铁虫 (初入文坛)
自己顶一下，希望能有大神来解答呀
赞一下 回复此楼
8楼2016-05-04 09:45:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

funnywjy

新虫 (初入文坛)
请问楼主还在么？有个小小的请求。我也需要用到HB模型，请问楼主的公式(包括各参数的公式)是在哪里找到的？厚脸皮求一下出处>_<
赞一下 回复此楼
9楼2016-07-05 12:58:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

预言tq

新虫 (著名写手)
请问问题解决了吗

发自小木虫Android客户端
赞一下 回复此楼
10楼2019-05-23 21:10:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangjg1991 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭