24小时热门版块排行榜    

查看: 1677  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 3/150 2026-08-22 14:02 by h4CP7TrQR8Lg
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-22 11:28 by h4CP7TrQR8Lg
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见