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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » CFD基本理论 » 帮忙修改一下圆柱绕流动网格udf，急
51/1返回列表
查看: 1378  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liubojia

新虫 (初入文坛)

[求助] 帮忙修改一下圆柱绕流动网格udf，急 已有2人参与

圆柱绕流动网格问题udf
中间位置圆柱的面积要分解成想x，y方向现在这么面积分量不知道怎么在udf中
CODE:
DEFINE_ON_DEMAND(print_force)
{
   real A[ND_ND];
   face_t f;
   real dfx=0,dfy=0;
   int zone_ID=3;
   Thread*thread_name=Lookup_Thread(domain,zone_ID);
   begin_f_loop(f,thread)
   {
         F_AREA(A,f,thread)
         dfx+=F_P(f,thread)*A[1];
         dfy+=F_P(f,thread)*A[0];
     }
     end_f_loop(f,thread)
     printf("x-force=%f y-force=%f,dfx,dfy)

上面A[0] 和A[1]分别是中间圆柱在Y方向和X方向的面积分量    目前就是这个A[0]和A[1]不知道应该怎么定义出来，求大神帮帮忙  ， \"帮忙修改一下圆柱绕流动网格udf，急\"
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-05-11 12:27:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aprilhello

新虫 (初入文坛)
您好，你运行成功了吗，能不能给我发一份udf，谢谢！

发自小木虫Android客户端
赞一下 回复此楼
7楼2017-03-28 10:24:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

gxw13131

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
不明白你的意思，但你的udf里那个thread_name得改成thread,其实就是一个遍历用的指针。单元面积分量就是通过F_AREA()这个宏得到的，放在A[]数组里面，A[0]为x方向，A[1]为y方向。三维的话以此类推。

[ 发自小木虫客户端 ]
赞一下 回复此楼
2楼2015-05-11 12:58:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liubojia

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by gxw13131 at 2015-05-11 12:58:09
不明白你的意思，但你的udf里那个thread_name得改成thread,其实就是一个遍历用的指针。单元面积分量就是通过F_AREA()这个宏得到的，放在A[]数组里面，A为x方向，A为y方向。三维的话以此类推。
...

现在的问题就是如何在F_AREA()中定义A[]这个你会不帮我下出来行不
然后我把金币给你  
F_AREA(A,f,thread)
dfx+=F_P(f,thread)*F_AREA(A,f,thread)[0];
  dfy+=F_P(f,thread)*F_AREA(A,f,thread)[1]
赞一下 回复此楼
3楼2015-05-11 14:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gxw13131

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liubojia: 金币+10, ★有帮助 2015-05-11 22:36:52
引用回帖:
3楼: Originally posted by liubojia at 2015-05-11 14:53:15
现在的问题就是如何在F_AREA()中定义A[]这个你会不帮我下出来行不
然后我把金币给你  
F_AREA(A,f,thread)
dfx+=F_P(f,thread)*F_AREA(A,f,thread);
  dfy+=F_P(f,thread)*F_AREA(A,f,thread)...

A[]你不是在前面已经定义了吗？F_AREA(A,f,thread)这个语句实际上已经把单元各个面积分量输出到A[]里去了。
你这个udf不是想计算圆柱的升力和阻力吗？完全没必要用udf。在monitors里面，残差下面那个create，就可以直接输出升力和阻力了，还可以保存到文件里方便后处理。
赞一下 回复此楼
4楼2015-05-11 19:05:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 303求调剂 +6 元夕元 2026-03-20 7/350 2026-03-25 12:00 by edmund7
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-25 4/200 2026-03-25 10:37 by barlinike
[考研] 求调剂 一志愿 本科 北科大 化学 343 +4 13831862839 2026-03-24 5/250 2026-03-25 09:47 by 无际的草原
[考研] 07化学280分求调剂 +7 722865 2026-03-23 7/350 2026-03-25 09:29 by aa331100
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 一志愿武理材料工程348求调剂 +6 ￣^￣゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 求调剂材料学硕080500，总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 考研调剂 +3 呼呼？~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
信息提示
关闭
请填处理意见
关闭