| 查看: 2616 | 回复: 11 | ||
[求助]
UDF关于圆球表面力的计算，Compute force and moment ，出现问题，求大神 已有2人参与
|
||
|
大家好： 小弟最近在用Compute_Force_And_Moment这个函数，进行力和力矩的计算。我将函数计算得到的结果与Fluent自身计算的结果进行对比，发现了很奇怪的现象。即二维情况下，x方向的力不正确，y方向力和z方向的力矩是正确的。而在三维情况下，得到的三个力和力矩都是正确的！！！下面是代码和截图！大家帮忙看看，卡了好久了。希望能解答！ #include "udf.h" #include"mem.h" #include"dynamesh_tools.h" #include"math.h" /************************************************************************/ static real vel[3]={0,0,0}; static real omega[3]={0,0,0}; static real body_centroid[3]={0,0,0}; real NV_VEC(force),NV_VEC(moment); /************************************************************************/ DEFINE_EXECUTE_AT_END(comeon) { FILE *fp=NULL; char filename1[]="f_x.txt"; char filename2[]="f_y.txt"; char filename3[]="m_z.txt"; Domain * domain; Thread * tf1; domain=Get_Domain(1); tf1=Lookup_Thread(domain,3); Compute_Force_And_Moment (domain, tf1, body_centroid, force, moment, TRUE); /************************************************************************/ Message ("force_x=%.1f ,force_y=%.1f ,force_z=%.1f,moment_z=%f\n",force[0],force[1],force[2],moment[2]); /************************************************************************/ fp=fopen("f_x.txt","a"  ;fprintf(fp,"%f\n",force[0]); fclose(fp); fp=fopen("f_y.txt","a" ;fprintf(fp,"%f\n",force[1]); fclose(fp); fp=fopen("m_z.txt","a" ;fprintf(fp,"%f\n",moment[2]); fclose(fp); }  图片1.png  图片2.png |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有169人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
8楼2016-10-15 17:08:45
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13069.3
- 红花: 26
- 帖子: 3231
- 在线: 697.9小时
- 虫号: 3071487
- 注册: 2014-03-19
- 专业: 园艺学与植物营养学
2楼2016-10-14 19:35:39
3楼2016-10-15 11:17:09
4楼2016-10-15 13:20:46
5楼2016-10-15 14:41:49
6楼2016-10-15 16:20:10
7楼2016-10-15 16:30:09
9楼2016-10-17 10:11:53
10楼2016-10-17 10:47:43