| 查看: 1825 | 回复: 0 | |||
lml_sky木虫 (正式写手)
|
[交流]
udf计算刚体力矩问题
|
|
最近一直在做，fluent动网格挡板开启的模拟，真的好累心，遇到了能遇到的各种麻烦 直接入主题吧，问题相关的核心代码如下： begin_f_loop (f, t) { F_AREA (A, f, t); F_CENTROID (B,f,t); m_w+=F_P (f, t)*(A[0]*(C_X-B[1])-A[1]*(C_Y-B[0])); } end_f_loop (f, t) 其中线程指针 t 指向的是一个平面二维矩形框，矩形框的周围充满了定向流动的流体，想求矩形框的四个面（也就是四个边）所受的流体的冲击力相对于点（C_X,C_Y）合力矩。其中（C_X,C_Y）在矩形的一条长边上。 不知道代码写的对不对，在网上一直在找求转动问题的udf，一直没找到，在看了《ANSYS流固耦合分析与工程实例》中5.4止回阀的动态分析时受了点启发，写出以上代码，只是运行后总出现负体积。那本书中给的例子是平动问题，转动是我自己按常识写的。另外用那书中给的例子的网格文件，我一步步的按着书上说明操作，也会出现负体积，实在不解。敢问有没有遇到相同情况的虫友啊，望同仁能多多讨论，相互启发！ |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
点击这里搜索更多相关资源