24小时热门版块排行榜    

查看: 581  |  回复: 2

lekuncui

新虫 (初入文坛)


[交流] udf编程求助

需要获取一个内部边界(interior)的力,下面程序出问题,#include "udf.h"
DEFINE_ON_DEMAND(on_demand_F)
{        
real velocity,flux,F=0;
Thread *t;
face_t f;/* 定义面变量 */
Domain *domain;/* 定义计算区域指针 */
domain=Get_Domain(1);/* 获取指针地址 */
t=Lookup_Thread(domain,10); /* 获取那个内部边界interior面的指针地址,即此面的ID */
begin_f_loop(f, t)
{ velocity=F_U(f,t); /* 获得网格点速度 */
flux=F_FLUX(f,t)*2*3.1415926;/* 获得网格点的流量 */
F+=velocity*flux;/* 计算力 */
}
end_f_loop(f, t)/* interior面循环结束 */
printf("F is %g\n",F);
}
Error:
FLUENT 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
求指教......
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

47425804

木虫 (正式写手)


★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+1, 谢谢参与! 2013-02-21 15:24:27
内部面不存在F_U宏。
2楼2013-02-21 14:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hiqun

禁虫 (正式写手)


★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+3, 谢谢参与! 2013-02-23 19:15:26
begin_f_loop(f, t)
{ velocity=F_U(f,t); /* 获得网格点速度 */
flux=F_FLUX(f,t)*2*3.1415926;/* 获得网格点的流量 */
F+=velocity*flux;/* 计算力 */
}
end_f_loop(f, t)/* interior面循环结束 */

define_demand命令中,计算力的话,是否可采用begin_c_loop,计算每个节点的力,求和,取平均值。具体例子,可参见udf help中的demand中计算温度的平均值得方式。
3楼2013-02-21 17:30:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lekuncui 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见