| 查看: 879 | 回复: 3 | |||
dinosaur1123新虫 (著名写手)
|
[求助]
UDF编程问题，试了好几天还是不对，求大神指导已有1人参与
|
|
想模拟一个管道旋转并且有水平方向的振动，即X方向做简谐运动同时绕管子轴心旋转。用滑移网格写了一个UDF，但是总是调试不对，旋转轴心的坐标不随时间变化，并且算到第11步就报错。 程序如下： #include "udf.h" #include "stdio.h" #define PI 3.141592654 DEFINE_ZONE_MOTION(originx_volecity, omega, axis, origin, velocity, time, dtime) { real f, omegar, theta0, R, rpm; real s; real dx, ts; s=RP_Get_Real("flow-time"  ;ts=s-0.09424100816249847; /* 这是由于前面没有加X方向的振动计算了一段时间0.09424100816249847s */ rpm = 972.683363; /* 转速 */ f = rpm/60; omegar = 2*PI*f; theta0 = PI/2; R = 0.001; /* 幅值 */ velocity[0] = -1.* omegar * R * sin(omegar * ts + theta0); velocity[1] = 0.0; velocity[2] = 0.0; dx=velocity[0]*dtime; origin[0] = origin[0]+dx; origin[1] = 0.0; origin[2] = 0.0; *omega = 2*PI*f; Message("\n originx=%f,dx=%f",origin[0],dx); } DEFINE_PROFILE(inlet_z_velocity, thread, index) /* 这是定义入口速度，和上面没关系 */ { real x[ND_ND]; real r,r0; face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); r0=x[0]*x[0]+x[1]*x[1]; r = sqrt(r0); F_PROFILE(f, thread, index) =0.51969*pow((1.0-r/0.025),0.14286); } end_f_loop(f, thread) } UDF编译完就先迭代了3次，这是什么原因？（见图1  ）然后就是origin[0]的值没有变化，第一步（图1）和第11步（图2  ）都是-0.000017另外，用全局变量把a=origin[0]的值存储再令origin[0]=a+dx也试过，也不行，求解释啊大神  |
回复此楼» 猜你喜欢
申请26博士
已经有5人回复
-
职称评审没过，求安慰
已经有22人回复
-
垃圾破二本职称评审标准
已经有15人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有20人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
dinosaur1123
新虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2261.6
- 散金: 395
- 红花: 11
- 帖子: 1002
- 在线: 155小时
- 虫号: 2192856
- 注册: 2012-12-18
- 性别: GG
- 专业: 污染控制化学
2楼2016-09-27 23:35:39
mycc
专家顾问 (正式写手)
-
专家经验: +181 - 应助: 230 (大学生)
- 金币: 2589.6
- 红花: 35
- 帖子: 531
- 在线: 58.8小时
- 虫号: 881523
- 注册: 2009-10-23
- 性别: GG
- 专业: 流体力学
- 管辖: 仿真模拟
3楼2016-09-28 08:23:02
dinosaur1123
新虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2261.6
- 散金: 395
- 红花: 11
- 帖子: 1002
- 在线: 155小时
- 虫号: 2192856
- 注册: 2012-12-18
- 性别: GG
- 专业: 污染控制化学
4楼2016-09-28 08:26:44