| 查看: 3096 | 回复: 9 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
如何用UDF求解面积分
|
|||||
|
如何用UDF求解面积分？ 我的问题是要给一个UDS变量定义其在边界的值等于另外一个UDS变量在某个面的积分值！ 帮助文档中有部分求体积分，但是面积分怎么求？ |
回复此楼» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有20人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有48人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
dai26520
专家顾问 (著名写手)
-
专家经验: +219 - 仿真EPI: 1
- 应助: 224 (大学生)
- 金币: 4097
- 散金: 4143
- 红花: 107
- 帖子: 1850
- 在线: 529.6小时
- 虫号: 1619980
- 注册: 2012-02-16
- 专业: 高分子合成化学
- 管辖: 仿真模拟
|
#include "udf.h" real NV_VEC(A); real avg_pressure; real sum_P_A=0.0; real sum_A=0.0; real outlet_down_ID=xx; Thread *thread_out; face_t f; Domain *domain; DEFINE_ADJUST(outlet_down,domain) { domain=Get_Domain(1); thread_out=Lookup_Thread(domain,outlet_down_ID);/*outlet_down_ID是出口面的id*/ begin_f_loop(f,thread_out) { F_AREA(A,f,thread_out); sum_A+=NV_MAG(A); sum_P_A+=NV_MAG(A)*F_P(f, thread_out); } end_f_loop(f,thread_out) avg_pressure=sum_P_A/sum_A; } DEFINE_PROFILE(pressure_down,t,i) { face_t f; begin_f_loop(f,t) { F_PROFILE(f,t,i) = avg_pressure; } end_f_loop(f,t) } 给你一个完整的例子参考下，算的平均压力。 |
6楼2016-08-03 15:48:21
mycc
专家顾问 (正式写手)
-
专家经验: +181 - 应助: 230 (大学生)
- 金币: 2589.6
- 红花: 35
- 帖子: 531
- 在线: 58.8小时
- 虫号: 881523
- 注册: 2009-10-23
- 性别: GG
- 专业: 流体力学
- 管辖: 仿真模拟
5楼2016-08-03 10:40:35
7楼2016-08-03 17:47:06
8楼2016-08-03 17:47:20
