| 查看: 559 | 回复: 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 |
回复此楼» 猜你喜欢
0703总分331求调剂
已经有6人回复
-
085600，321分求调剂
已经有7人回复
-
325求调剂
已经有4人回复
-
085600，320分求调剂
已经有12人回复
-
一志愿0817化学工程与技术，求调剂
已经有27人回复
-
一志愿北京化工085600 310分求调剂
已经有8人回复
-
一志愿哈尔滨工业大学085600英一数二337分求调剂
已经有10人回复
-
22408求调剂 354分 可跨专业
已经有3人回复
-
277工科求调剂
已经有4人回复
-
一志愿沪9，求生物学调剂，326分
已经有3人回复
