| 查看: 637 | 回复: 0 | ||
[求助]
二维剖面强迫横摇问题
|
|
大家好，第一次在小木虫上发帖求助，金币有点少，但是真心求教，求大家帮帮忙。 我现在做的问题如下： 想要用FLUENT计算一个二维剖面在水中强迫横摇的问题。二维剖面做简谐运动，幅值频率分别是A和φ0。应用动网格中的滑移网格来解决这个问题。下面是我自己编写的UDF。 UDF code: #include"udf.h" #define freq 0.84 //频率 #define ampl 0.175 //幅值 #define density 998.2 //水密度 #define GRAVITY 9.81 //重力加速度 int id=14; //二维剖面的ID DEFINE_CG_MOTION(int_block_roll, dt, vel, omega, time, dtime) { omega[2] = 2*M_PI*freq*ampl*cos(2*M_PI*freq*time); } DEFINE_EXECUTE_AT_END(exe_end) { real x_cg[3],m_glob[3],f_glob[3],sita,time,m_plus; real x[ND_ND],A[ND_ND]; FILE *fp; Domain *domain=Get_Domain(1); face_t f; Thread *tf=Lookup_Thread(domain,id); time = CURRENT_TIME; x_cg[0] = 0.0; x_cg[1] = 0.0; m_plus=0.0; sita = ampl*sin(2*M_PI*freq*time); if(!Data_Valid_P()) return; Compute_Force_And_Moment(domain,tf,x_cg,f_glob,m_glob,TRUE); begin_f_loop(f,tf) { F_CENTROID(x,f,tf); F_AREA(A,f,tf); if(x[1]<0.0) m_plus+=density*GRAVITY*x[1]*fabs(x[1]*A[0]-x[0]*A[1]); } end_f_loop(f,tf) m_glob[2]+=m_plus; #if !RP_NODE fp=fopen("roll.txt","a"  ;fprintf(fp,"%f,%f,%f,%f,%f,%f ",time,sita,f_glob[0],f_glob[1],NV_MAG(f_glob),m_glob[2]); fclose(fp); #endif } 模拟中我使用的二维剖面是根据Ikeda的试验来确定的。网格的划分使用ANSYS ICEM来画的结构网格。参考很多做同样问题的文献，最后是要监测二维剖面受到的力矩的大小。求大家帮忙看看哪里出了问题，谢谢！ 万分感谢！@wuming524 |
回复此楼» 猜你喜欢
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有11人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
20