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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 读入UDF，一初始化FLUENT就闪退
北京石油化工学院2026年研究生招生接收调剂公告
52/1返回列表
查看: 3370  |  回复: 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的回帖

新世界2019

新虫 (小有名气)
引用回帖:
10楼: Originally posted by 预言tq at 2019-05-23 21:10:17
请问问题解决了吗

请问解决了吗？

发自小木虫Android客户端
回复此楼
11楼2020-02-29 12:34:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

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的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 招收生物学/细胞生物学调剂 +4 IceGuo 2026-03-26 5/250 2026-03-29 01:25 by griffith2014
[考研] 070305高分子化学与物理 304分求调剂 +7 c297914 2026-03-28 7/350 2026-03-29 01:17 by 544594351
[考研] 348求调剂 +3 小懒虫不懒了 2026-03-28 3/150 2026-03-29 00:39 by 544594351
[考研] 298求调剂 +3 种圣赐 2026-03-28 3/150 2026-03-28 23:58 by wxiongid
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 311（085601）求调剂 +4 liziyeyeye 2026-03-28 4/200 2026-03-28 18:50 by 535743368
[考研] 一志愿华北电力大学能动专硕，293，求调剂 +3 15537177284 2026-03-23 5/250 2026-03-28 16:11 by xxxsssccc
[考研] 320分，材料与化工专业，求调剂 +9 一定上岸aaa 2026-03-27 13/650 2026-03-28 15:00 by 神马都不懂
[考研] 329求调剂 +6 星野? 2026-03-26 6/300 2026-03-28 14:14 by 唐沐儿
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[考研] 一志愿华东理工大学081700，初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 343求调剂 +4 赠我一本书 2026-03-23 4/200 2026-03-27 00:40 by wxiongid
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 085601求调剂总分293英一数二 +4 钢铁大炮 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
信息提示
关闭
请填处理意见
关闭