| 查看: 2632 | 回复: 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 |
回复此楼» 猜你喜欢
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有147人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有17人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
9楼2016-10-17 10:11:53
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13154.3
- 红花: 26
- 帖子: 3250
- 在线: 698.5小时
- 虫号: 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