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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent udf compile失败
91/1返回列表
查看: 2401  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhoufu1989

银虫 (小有名气)

[求助] fluent udf compile失败

才接触udf不久的新人，希望各位前辈指点迷津！！
模型如图示：，要对红色的圆柱壁面添加热流。
http://bcs.duapp.com/emuchnet/20 ... _1369793765_108.png，热流随角度的关系如图所示，采用分段多项式拟合。
由于我采用的是直角坐标，所以上面拟合的关系式，要将角度（圆柱面半径25mm）转换成x、y坐标，编写的udf如下：
#include "udf.h"
#define A 57.295709
#define R 25.000
DEFINE_PROFILE(heat_flux,thread,index)
{
        real x[ND_ND];
        real a,b,temp1,temp2,temp3,temp4,temp5,temp6,temp7,temp8;
        face_t f;
        begin_f_loop(f,thread)
        {
                F_CENTROID(x,f,thread);
                a=x[0];
                b=x[1];
                if(a>=0&b>=0)
                {
                        temp1=A*asin(b/R);
                        F_PROFILE(f,thread,index)=773.08579+6.5979*temp1-0.01532*temp1*temp1-0.00467*temp1*temp1*temp1+0.000033032*temp1*temp1*temp1*temp1;
                }
                else if(b>0&-16.8897<=a<0)
                {
                        temp2=180-A*asin(b/R);
                        F_PROFILE(f,thread,index)=1310690.000-37867.31033*temp2+353.99224*temp2*temp2-1.0564*temp2*temp2*temp2;
                }
                else if(b>0&-24.0908<=a<-16.8897)
                {
                        temp3=180-A*asin(b/R);
                        F_PROFILE(f,thread,index)=485341.70814-7575.60803*temp3+43.67718*temp3*temp3-0.08533*temp3*temp3*temp3;
                }
                else if(b>=0&-25.0<=a<-24.0908)
                {
                        temp4=180-A*asin(b/R);
                        F_PROFILE(f,thread,index)=-21122000.000+492125.95098*temp4-4243.96404*temp4*temp4+16.09596*temp4*temp4*temp4-0.02267*temp4*temp4*temp4*temp4;
                }
                else if(b>=0&-23.7081<=a<-25.0)
                {
                        temp5=180-A*asin(b/R);
                        F_PROFILE(f,thread,index)=-21122000.000+492125.95098*temp5-4243.96404*temp5*temp5+16.09596*temp5*temp5*temp5-0.02267*temp5*temp5*temp5*temp5;
                }
                else if(b<0&-23.7081<a<=-16.5655)
                {
                        temp6=180-A*asin(b/R);
                        F_PROFILE(f,thread,index)=164115.58656-971.46727*temp6-0.07801*temp6*temp6+0.0094*temp6*temp6*temp6;
                }
                else if(b<0&-16.5655<a<=0)
                {
                        temp7=180-A*asin(b/R);
                        F_PROFILE(f,thread,index)=-15203500.000+187415.08644*temp7-761.84536*temp7*temp7+1.0232*temp7*temp7*temp7;
                }
                else
                {
                        temp8=360+A*asin(b/R);
                        F_PROFILE(f,thread,index)=337883.03869-4342.70635*temp8+20.62194*temp8*temp8-0.04289*temp8*temp8*temp8+0.000033032*temp8*temp8*temp8*temp8;
                }
        }
        end_f_loop(f,thread)
}
interpret不报错，但是算了一下，最后热流没加上去。
compile出现如下错误：
已复制         1 个文件。
(system "copy D:\software\fluent\fluent\Fluent.Inc\fluent6.3.26\src\makefile_nt.udf libudf\ntx86\3d\makefile"
已复制         1 个文件。
(chdir "libudf"()
(chdir "ntx86\3d"()

Done.
"f:/gambit/sample"

Opening library "libudf"...
Error: open_udf_library: 系统找不到指定的文件。

Error Object: ()

谢谢大家了！！
fluent udf compile失败
模型图.png


fluent udf compile失败-1
b.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
早起的鸟儿有虫吃……
1楼 2013-05-29 10:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoufu1989

银虫 (小有名气)
原先只用一个temp的，但是后来提醒有覆盖，所以用了8个temp变量的
赞一下 回复此楼
早起的鸟儿有虫吃……
2楼2013-05-29 10:34:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Scorpio007

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-06-03 21:15:48
如果确定正确的udf是可以编译的，那么else if(b>=0&-23.7081<=a<-25.0)a的范围有问题，改一下试试……
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

zhoufu1989
……
3楼2013-05-29 21:53:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoufu1989

银虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by Scorpio007 at 2013-05-29 21:53:49
如果确定正确的udf是可以编译的，那么else if(b>=0&-23.7081<=a<-25.0)a的范围有问题，改一下试试……

我用ansys fluent能够加载热流了，现在出现新的问题了，热流加的不对，出现了负值，与我期望的不一样。这是怎么回事？希望您能帮帮我
赞一下 回复此楼
早起的鸟儿有虫吃……
4楼2013-05-31 16:19:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-06-02 11:07:11
zhoufu1989: 金币+50, ★★★★★最佳答案, 帮我解决了问题，非常感谢！ 2013-06-02 19:31:26
zhoufu1989: 回帖置顶 2013-06-02 19:31:30
粗略扫了一下程序，没细究，发现你都把逻辑与（&&）写成了按位与（&）,这个对你判断有影响没？
赞一下 回复此楼
5楼2013-05-31 17:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)
★ ★
xiegangmai: 金币+2, 谢谢参与！ 2013-06-02 11:07:17
zhoufu1989: 回帖置顶 2013-06-02 19:31:50
关系写法有问题，C语言和数学表达是不一样的，在你的程序中，举例来说：
else if(b>0&-16.8897<=a<0)  这句，先不管意义，正确语法应该为：
else if(b>0&&-16.8897<=a&&a<0)
   程序中像这样的错误不在少数。
赞一下 回复此楼
6楼2013-05-31 17:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoufu1989

银虫 (小有名气)
引用回帖:
6楼: Originally posted by fangdian at 2013-05-31 17:05:50
关系写法有问题，C语言和数学表达是不一样的，在你的程序中，举例来说：
else if(b>0&-16.8897<=a<0)  这句，先不管意义，正确语法应该为：
else if(b>0&&-16.8897<=a&&a< ...

您说的太对了，&&，这个我知道，-16.8897<=a<0，我把这个与fortran用混了，怪不得老是有问题，现在改过来后能加热流了，从云图上看是正确的。真的太感谢您了
赞一下 回复此楼
早起的鸟儿有虫吃……
7楼2013-06-02 19:28:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anieping

铁虫 (初入文坛)
引用回帖:
7楼: Originally posted by zhoufu1989 at 2013-06-02 19:28:45
您说的太对了，&&，这个我知道，-16.8897<=a<0，我把这个与fortran用混了，怪不得老是有问题，现在改过来后能加热流了，从云图上看是正确的。真的太感谢您了...

请教下楼主   加了热流之后   怎么查看所施加的热流的云图或XY plot呀？谢谢！
赞一下 回复此楼
8楼2013-10-09 10:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

丑小鸭亦潇洒

木虫 (著名写手)

科研小虫
请问楼主：这个UDF代码错在那里了呢？最后调整出来了吧，可否分享学习下啊，求教呢
赞一下 回复此楼
9楼2017-07-10 22:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhoufu1989 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-15 15:35 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 14:14 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:14 by C79jjtjAKjEn
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:14 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
关闭
请填处理意见
关闭