| 查看: 2655 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
fluent udf三维力源项问题求助 已有1人参与
|
|||
|
采用由于使用的是二维模型，进行旋转模拟三维流动，分析出受力形式，并表达成柱坐标的方式，因此需要对三方向源项进行加载，新人第一次弄这个，求大神指点 #include"udf.h" #include"math.h" #define SINB 0.707 /* 叶片安装角度，45°是为0.707；30°时为0.5；60°时为0.866； */ #define COSB 0.707 /* 叶片安装角度，45°是为0.707；30°时为0.866；60°时为0.5； */ double value() { double m; if(C_U(c,t)==0) { if(C_V(c,t)>0) { m=1.57; } if(C_U(c,t)<0) { m=-1.57; } else { m=atan(C_V(c,t)/C_U(c,t)) ; } return m; } DEFINE_SOURCE(xmom_source,c,t,dS,eqn) { real source; source = -68.38+68.38*COSB; return source; } DEFINE_SOURCE(swirl_mom_source,c,t,dS,eqn) { real source; source =value()-1.57; return source; } DEFINE_SOURCE(radial_mom_source,c,t,dS,eqn) { real source; source =68.38*SINB; return source; } 现在无法导入，伤心啊 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fluent使用UDF时想要编译速度随时间变化函数
已经有14人回复
- 急，在FLUENT中如何同时调用两个UDF？
已经有11人回复
- FLUENT三维建模后入口速度为抛物线怎么写UDF？
已经有8人回复
- Fluent中 UDF ADJUST 宏链接问题
已经有7人回复
- fluent UDF 修改物性
已经有13人回复
- fluent中UDF编写自定义矢量——UDF大牛请进
已经有7人回复
- fluent中UDF编译的问题！
已经有6人回复
- fluent 加载udf 失败 请高手不吝赐教
已经有8人回复
- fluent udf compile失败
已经有8人回复
- 关于fluentUDF的一个小问题！
已经有5人回复
- 学习fluent udf
已经有5人回复
- fluent udf请教
已经有16人回复
- fluent udf的应用
已经有20人回复
- 走过路过都看看吧，fluent udf 是个神奇的语言
已经有20人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- 求助！！！我想要在fluent中针对一个面进行积分，那么UDF应该用什么宏函数？怎么用？
已经有5人回复
- fluent中用udf编写三维速度边界条件出现错误，大家来指导一下啊
已经有9人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
|
对编程一窍不通。。。。 #include"udf.h" #include"math.h" #define SINB 0.707 /* 叶片安装角度，45°是为0.707；30°时为0.5；60°时为0.866； */ #define COSB 0.707 /* 叶片安装角度，45°是为0.707；30°时为0.866；60°时为0.5； */ real value(cell_t c, Thread *t) { real m; if(C_U(c,t)==0) { if(C_V(c,t)>0) { m=1.57; } if(C_U(c,t)<0) { m=-1.57; } else { m=atan(C_V(c,t)/C_U(c,t)) ; } return m; } DEFINE_SOURCE(xmom_source,c,t,dS,eqn) { real source; source = -68.38+68.38*COSB; return source; } DEFINE_SOURCE(swirl_mom_source,c,t,dS,eqn) { real source; source =value(c,t,)-1.57; dS[eqn]=0; return source; } DEFINE_SOURCE(radial_mom_source,c,t,dS,eqn) { real source; source =68.38*SINB; return source; } |
5楼2014-11-21 14:43:27
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
2楼2014-11-21 13:24:32
其实没看懂，按你的改一下，还是有错误，能否详细帮忙指点一下，先谢谢了 3楼2014-11-21 14:09:17
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
4楼2014-11-21 14:24:00
