24小时热门版块排行榜    

查看: 1505  |  回复: 1

shenjiawei

新虫 (初入文坛)

[求助] 求助求助!Fluent六自由度UDF问题。

本人想模拟波浪流经船体时发生的自由横摇运动。但模拟过程中,发现船舶随重心的旋转为单向,直接翻转,然后出现负体积计算终止,不太符合实际。请问UDF中是否哪里编写有误?
另外,请问如何监测船体横摇角度的变化,如何输出?

附上我编写的UDF,求大神们帮忙看一看。跪谢!!!
[code]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenjiawei

新虫 (初入文坛)

接上楼。UDF代码如下,求各位大神帮忙看看。

#include "udf.h"



#define m      72697.4721; /* 定义全局变量,mass质量通过自主设定的m*求得,我用的是铝的质量 */



#if !RP_NODE   

FILE *result0;
#endif

DEFINE_SDOF_PROPERTIES(m_first0, prop, dt, time, dtime)
{

  real x_cg[ND_ND];
  real f_glob0[ND_ND],m_glob0[ND_ND];

         
                   #if !RP_HOST  
  
                 

                                          
                 

                  Domain *domain= Get_Domain (1);
                  Thread  *tf0= Lookup_Thread (domain, 6);/* fyq1*/               
               

               
               
       

        x_cg[0]=DT_CG(dt)[0];
        x_cg[1]=DT_CG(dt)[1];


          
                  #endif                 
                                          prop[SDOF_MASS]   =72697.4721;

                                   prop[SDOF_IZZ]    =3932.34;
                                  
                                                      #if !RP_HOST             
              Compute_Force_And_Moment (domain, tf0, x_cg, f_glob0, m_glob0, FALSE);
             
              domain= Get_Domain (1);
                              #endif          

             
         node_to_host_real(f_glob0, ND_ND);
     node_to_host_real(m_glob0, ND_ND);


     
   

     prop[SDOF_LOAD_F_X] =0;
     prop[SDOF_LOAD_F_Y] =f_glob0[1];

     prop[SDOF_LOAD_M_Z]= m_glob0[2];  
  
#if !RP_NODE
                result0=fopen("m1_yundong.dat","a";
               


   fprintf(result0," %10.6f  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f  %10.6f \n",time,DT_CG(dt)[0],DT_CG(dt)[1],DT_VEL_CG(dt)[0],DT_VEL_CG(dt)[1],f_glob0[0],f_glob0[1]);
                               
                fclose(result0);       
     #endif          

}
2楼2020-01-14 20:09:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shenjiawei 的主题更新
信息提示
请填处理意见