| 查看: 2878 | 回复: 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 |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有197人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
8楼2016-10-15 17:08:45
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13285.8
- 红花: 26
- 帖子: 3459
- 在线: 717.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
28