24小时热门版块排行榜

返回列表

lml_sky

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 4151.8
帖子: 659
在线: 40.6小时
虫号: 3048420
注册: 2014-03-13
专业: 生物医学系统建模及仿真

[交流] 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止回阀的动态分析时受了点启发，写出以上代码，只是运行后总出现负体积。那本书中给的例子是平动问题，转动是我自己按常识写的。另外用那书中给的例子的网格文件，我一步步的按着书上说明操作，也会出现负体积，实在不解。敢问有没有遇到相同情况的虫友啊，望同仁能多多讨论，相互启发！

回复此楼