| 查看: 2682 | 回复: 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 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有166人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
8楼2016-10-15 17:08:45
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13780.3
- 红花: 26
- 帖子: 3312
- 在线: 706小时
- 虫号: 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
5