| 查看: 2663 | 回复: 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 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有89人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有27人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有44人回复
8楼2016-10-15 17:08:45
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13419.3
- 红花: 26
- 帖子: 3293
- 在线: 699.2小时
- 虫号: 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
