| 查看: 2327 | 回复: 5 | |||
2265973243银虫 (小有名气)
|
[求助]
利用UDF控制旋转机械被动旋转问题请教各位前辈 已有3人参与
|
|
给为前辈好，我最近在做一个旋转机械方面的问题，先给定旋转机械一个初始角速度，再通过UDF计算出旋转机械受到的力矩，利用这个力矩通过相应方程计算出在一个时间步内角速度的增量，再把原来的角速度加上这个增量得到一个新的角速度，依次进行下去，我编写的UDF程序如下： #include "udf.h" static real omega_z=0.0,J=200; DEFINE_CG_MOTION(ZL,dt,vel,omega,time,dtime) { face_t f; cell_t c; real f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND],domega_z; Domain *domain = Get_Domain (1); Thread *tf1 = Lookup_Thread (domain,1); int i; NV_S(omega, =, 0.0); if (!Data_Valid_P()) return; for(i=0;i<=ND_ND;i++) { f_glob=0; m_glob=0; } for(i=0;i<ND_ND;i++) x_cg=DT_CG(dt); if(time<=0.0002) { omega[2]=-29.87; Message ("\ntime=%f ,x=%.10lf ,y=%.10lf , force_x=%.1f ,force_y=%.1f ,moment_x=%f,moment_y=%f,moment_z=%f ,omega_z=%f\n", time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2],omega_z); return; } Compute_Force_And_Moment (domain, tf1, x_cg, f_glob, m_glob, TRUE); domega_z=m_glob[2]*0.0002/J;//计算得到角速度增量 omega_z+=domega_z;//原来的角速度加上增量 Message ("\ntime=%f ,x=%.10lf ,y=%.10lf , force_x=%.1f ,force_y=%.1f ,moment_x=%f,moment_y=%f,moment_z=%f ,omega_z=%f\n", time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2],omega_z); omega[2]=omega_z; } 程序可以正常运行，但就是计算出的角速度不对，我的初始角速度是--29.87，但是这里最终得到的角速度好像只有角速度增量的之和； 如果我根据上一时间步计算出了角速度的增量，那么我要得到上一时间步角速度和这个增量之和（也即当前时间步的角速度），请问这个要怎么修改呢？非常感谢，祝生活愉快！！ |
回复此楼» 猜你喜欢
重要来源：本周末出结果
已经有3人回复
-
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
xingfuww
专家顾问 (正式写手)
-
专家经验: +78 - 仿真EPI: 2
- 应助: 46 (小学生)
- 金币: 6637.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 284小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
2楼2016-08-22 08:47:38
xingfuww
专家顾问 (正式写手)
-
专家经验: +78 - 仿真EPI: 2
- 应助: 46 (小学生)
- 金币: 6637.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 284小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
3楼2016-08-22 10:09:29
4楼2016-10-24 17:12:55
5楼2016-10-24 17:56:17
mycc
专家顾问 (正式写手)
-
专家经验: +181 - 应助: 230 (大学生)
- 金币: 2589.6
- 红花: 35
- 帖子: 531
- 在线: 58.8小时
- 虫号: 881523
- 注册: 2009-10-23
- 性别: GG
- 专业: 流体力学
- 管辖: 仿真模拟
6楼2016-10-25 10:31:48
