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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 读入UDF，一初始化FLUENT就闪退
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 3364  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 学硕274求调剂 +8 Li李鱼 2026-03-26 8/400 2026-03-28 16:53 by xuxiang
[考研] 一志愿厦门大学化学学硕307求调剂 +10 y7czhao 2026-03-26 10/500 2026-03-28 14:23 by 唐沐儿
[考研] 一志愿武汉理工，总分321，英一数二，求老师收留。 +6 nnnnnnn5 2026-03-25 6/300 2026-03-28 13:33 by 果果妈咪
[考研] 311求调剂 +4 冬十三 2026-03-24 4/200 2026-03-28 13:17 by 唐沐儿
[考研] 070300求调剂306分 +4 26要上岸 2026-03-27 4/200 2026-03-28 13:06 by 唐沐儿
[考研] 291求调剂 +7 孅華 2026-03-22 7/350 2026-03-28 04:02 by fmesaito
[考研] 考研调剂 +4 Sanmu-124 2026-03-26 4/200 2026-03-27 17:49 by kiokin
[考研] 求调剂，一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 材料与化工328分调剂 +6 。，。，。，。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 296求调剂 +4 汪！？！ 2026-03-25 7/350 2026-03-25 16:41 by 汪！？！
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 277分求调剂，跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
信息提示
关闭
请填处理意见
关闭