24小时热门版块排行榜    

查看: 1987  |  回复: 6

butter_L

新虫 (小有名气)

[求助] UDF新手,解释没有问题,计算第一步开始报错 已有1人参与

UDF新手。自己在做固液两相流对管内壁的磨损的模拟。用的fluent,写了一小段UDF想试一试,Interpreted没有没有错误,hook之后,开始计算第一步就出现错误

Error:
C:\PROGRA~1\ANSYSI~1\v145\fluent\fluent14.5.0\win64\3ddp\fl1450s.exe received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: #f

按照网上说的,UDM已经分配了,初始化之前导入UDF或者初始化之后导入UDF,将已经计算几步之后的数据导入、再导入UDF,这些方法都已经尝试过,仍然报错。

这是我的UDF和case,大神帮忙看下,不胜感激!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : fluent1.rar
  • 2019-09-20 16:24:21, 9.57 M

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

butter_L

新虫 (小有名气)

#include "udf.h"
#include "math.h"
#include "metric.h"
#include "mem.h"
DEFINE_ADJUST(www,d)
{
        real w;
    int n;
    real densic=3200;
    real c;
    real alpha;
    real A[ND_ND];
    real vel[3];
    real velm;
    real Am;
    real co;
    real velx;
    real vely;
    real velz;
    real qm;
    Thread *tf;
    face_t f;
    c=sqrt(3*480*1000000/densic);
    thread_loop_f(tf,d)
    {
            begin_f_loop(f,tf)
            {
                real denmix=1300;
                velx=F_U(f,tf);
            vely=F_V(f,tf);
            velz=F_W(f,tf);
            vel[1]=velx;
            vel[2]=vely;
            vel[3]=velz;
            qm=velx*denmix*3.1416*0.0045*0.0045;
            F_AREA(A,f,tf); /*面法线向量*/
                velm=NV_MAG(vel);
                Am=NV_MAG(A);
                co=NV_DOT(vel,A)/(velm+0.1)/(Am+0.1); /*面法线向量与速度向量夹角的余弦值*/
                alpha=asin(co)*180/3.1416; /*速度向量与面的夹角度数*/
                w=7/3.1416*qm/denmix*pow((velm/c),2.5)*sin(2*alpha)*sin(alpha); /*单位时间切削量m3*/
                F_UDMI(f,tf,0)=w; /*计算结果储存到0文件*/
            }
            end_f_loop(f,tf)
    }
}
2楼2019-09-20 16:27:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfuww

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
NV_MAG 貌似只能用来计算 内部传递的矢量,你试试自己把求模 用平方开根号来写出来。尽量别用它自带的函数。免得出错
3楼2019-09-20 18:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfuww

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by butter_L at 2019-09-20 16:27:12
#include "udf.h"
#include "math.h"
#include "metric.h"
#include "mem.h"
DEFINE_ADJUST(www,d)
{
        real w;
    int n;
     ...

还有 UDMI 函数需要初始化的,用DEFINE_ON_DEMAND来初始化一下。
4楼2019-09-20 18:13:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfuww

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
butter_L: 金币+20, ★★★很有帮助 2019-09-21 14:18:33
我已经简单测试了,你的udmi不需要初始化也可以,因为你第一步没调用到,你的UDMI 设置返回值为0 的时候是没问题的,程序能跑。说明你的计算公式有问题,如前面的回复,自己把计算方程写进去就好,别用自带的函数,免得出错。
还有你的UDMI 有没有启用啊,在user-defined 里面 memory 里面 number of user defined memory locations 把0 改成1.
5楼2019-09-20 18:24:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfuww

专家顾问 (正式写手)

我看你是遍历面 f? 改成遍历cell 试试。我这边试了用cell 可以,
6楼2019-09-20 18:44:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

butter_L

新虫 (小有名气)

引用回帖:
5楼: Originally posted by xingfuww at 2019-09-20 18:24:38
我已经简单测试了,你的udmi不需要初始化也可以,因为你第一步没调用到,你的UDMI 设置返回值为0 的时候是没问题的,程序能跑。说明你的计算公式有问题,如前面的回复,自己把计算方程写进去就好,别用自带的函数, ...

好的,感谢!我试一下
7楼2019-09-21 08:25:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 butter_L 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 06:29 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 06:13 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 06:00 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 05:34 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
信息提示
请填处理意见