求助！！！udf自定义采空区孔隙率
我做采空区瓦斯模拟，自定义采空区孔隙率。x为采空区后方，y为工作面平行方向，工作面长度100m，多孔介质80m高。代码：
#include "udf.h"
DEFINE_PROPERTY(mykongxilv, cell, thread)
{
int i;
real kongxilv;
real x;
real xx1,xx2,xx3;
cell_t c;
Thread *t;
face_t f;
begin_f_loop(f,thread)
{
C_CENTROID(x,c,t);
real xx1=x;x向采空区后方，长度150
real xx2=x;y向，工作面方向，工作面长100m
real xx3=x;z向，高，80m
if((xx1<=150)&&(xx1>=0)&&(xx2<=50)&&(xx2>=-50))
{
if((xx3<=80)&&(xx3>=0)&&(xx2>=0))
kongxilv=((0.2*exp(-0.0221*xx1)+0.1)*(exp(-0.15*(50-xx2))+1));
else if((xx3<=80)&&(xx3>=0)&&(xx2<=0))
kongxilv=((0.2*exp(-0.0221*xx1)+0.1).*(exp(-0.15*(50+xx2))+1));
}
else kongxilv=0;
return kongxilv;
}
end_f_loop(f,thread)
}
错误提示：
>
cpp -I"C:\Fluent.Inc\fluent6.3.26/src" -I"C:\Fluent.Inc\fluent6.3.26/cortex/src" -I"C:\Fluent.Inc\fluent6.3.26/client/src" -I"C:\Fluent.Inc\fluent6.3.26/multiport/src" -I. -DUDFCONFIG_H="<udfconfig-host.h>" "D:\GAMBITMODEL\mykongxilv.c"
Error: D:\\GAMBITMODEL\\mykongxilv.c: line 14: parse error.
Error: D:\\GAMBITMODEL\\mykongxilv.c: line 15: parse error.
Error: D:\\GAMBITMODEL\\mykongxilv.c: line 16: parse error.
Error: D:\\GAMBITMODEL\\mykongxilv.c: line 22: parse error.
cpp -I"C:\Fluent.Inc\fluent6.3.26/src" -I"C:\Fluent.Inc\fluent6.3.26/cortex/src" -I"C:\Fluent.Inc\fluent6.3.26/client/src" -I"C:\Fluent.Inc\fluent6.3.26/multiport/src" -I. -DUDFCONFIG_H="<C:/Users/ADMINI~1/AppData/Local/Temp/udfconfig-node0.h>" "C:/User
/ADMINI~1/AppData/Local/Temp/mykongxilv.c.0.c"
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.0.c: line 14: parse error.
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.0.c: line 15: parse error.
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.0.c: line 16: parse error.
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.0.c: line 22: parse error.
cpp -I"C:\Fluent.Inc\fluent6.3.26/src" -I"C:\Fluent.Inc\fluent6.3.26/cortex/src" -I"C:\Fluent.Inc\fluent6.3.26/client/src" -I"C:\Fluent.Inc\fluent6.3.26/multiport/src" -I. -DUDFCONFIG_H="<C:/Users/ADMINI~1/AppData/Local/Temp/udfconfig-node1.h>" "C:/User
/ADMINI~1/AppData/Local/Temp/mykongxilv.c.1.c"
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.1.c: line 14: parse error.
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.1.c: line 15: parse error.
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.1.c: line 16: parse error.
Error: C:/Users/ADMINI~1/AppData/Local/Temp/mykongxilv.c.1.c: line 22: parse error.
回复 1# 龙biscuit 的帖子
我的模型y向是100m，这个附图是160的例子。我的x是150m这个附图也是150m。上面正负号“士”号取值;当y为负值时取正号,当y为正值时取负号回复 2# 龙biscuit 的帖子
你的注释部分要用/* */real xx1=x;/*x向采空区后方，长度150*/
real xx2=x;/*y向，工作面方向，工作面长100m*/
real xx3=x;/*z向，高，80m*/
否则编译不能通过。
回复 3# classic1573 的帖子
哦，多谢，我昨天把后面的注释去掉了就行了回复 3# classic1573 的帖子
大神，在帮我看看吧，我的代码如下/* 下面是对坐标轴在进风巷口的多孔介质孔隙率编程*/
#include "udf.h"
#include "metric.h"
#defineL 100
DEFINE_PROPERTY(mycell_poriness,cell,thread)
{
real myc_p;
real x;
cell_t c;
Thread *t;
real xx1,xx2,xx3;
begin_c_loop(c,thread)
{
C_CENTROID(x,c,t)
xx1=x;/*x轴方向，x指向采空区后方，长150m*/
xx2=x;/*y轴方向，y指向工作面平行方向，及工作面长100m，左右各50m*/
xx3=x;/*z轴方向，z指向采空区高度方向，多孔介质高20m*/
if((xx1<=150)&&(xx1>=0)&&(xx2<=50)&&(xx2>=-50))
{
if((xx3<=20)&&(xx3>=0)&&(xx2>=50))
myc_p=((0.2*exp(-0.0221*xx1)+0.1)*(exp(-0.15*(L-xx2))+1));
else if((xx3<=20)&&(xx3>=0)&&(xx2<=50))
myc_p=((0.2*exp(-0.0221*xx1)+0.1)*(exp(-0.15*(L+xx2))+1));
}
else myc_p=0;
return myc_p;
}
end_c_loop(c,thread)
}
这次用个小的模型实验，发现打开组分运输，一插入udf后就会出现the FLUENT process could not be started。这是什么原因啊，我的case和*.c放在同一个目录了，网格质量也很好。
回复 5# 龙biscuit 的帖子
原因找到了，是因为我的坐标轴方向搞错了。 师兄，您做过的源文件可否分享一下，也是做这个的，UDF、网格都编好了，就是不知道具体在Fluent里怎么设置，QQ是447848206，先谢谢师兄:$龙biscuit，您做过的源文件可否分享一下吗？练练手了。qq7549228
前辈，我最近也在学习编写孔隙率udf, 想请教两个问题：
1. 定义时应该用DEFINE_PROPERTY还是DEFINE_PROFILE呢？
2. udf成功导入fluent后，应该在Porosity中设置孔隙率为多少，我发现Porosity只能设置为constant或者new input parameter，没有user-defined。
页:
[1]