| 查看: 3087 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
如何用UDF求解面积分
|
||||
|
如何用UDF求解面积分？ 我的问题是要给一个UDS变量定义其在边界的值等于另外一个UDS变量在某个面的积分值！ 帮助文档中有部分求体积分，但是面积分怎么求？ |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
9楼2016-08-04 08:10:02
mycc
专家顾问 (正式写手)
-
专家经验: +181 - 应助: 230 (大学生)
- 金币: 2589.6
- 红花: 35
- 帖子: 531
- 在线: 58.8小时
- 虫号: 881523
- 注册: 2009-10-23
- 性别: GG
- 专业: 流体力学
- 管辖: 仿真模拟
5楼2016-08-03 10:40:35
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
7楼2016-08-03 17:47:06
100