【技术实现步骤摘要】
一种大幅六自由度运动的流固耦合模拟方法
本专利技术涉及的是一种计算机辅助工程(CAE)、计算流体力学(CFD)、流固耦合方法。
技术介绍
随着计算机技术和数值方法的飞速发展，计算流体力学方法(CFD)开始活跃于空气动力学和船舶与海洋工程的各个领域当中，已有越来越多研究工作基于CFD方法开展研究。然而CFD的发展一直受困于多个瓶颈，其中一个关键问题就是动网格技术。在结构物流固耦合运动问题当中，需要处理其六个自由度的运动。结构物在流体中具有六个自由度的运动，并且在流体的作用下，其运动规律极其复杂。因此，要准确地模拟计算结构物的运动，并考虑流固耦合作用，就必须要让网格“动”起来。传统的动网格技术主要有两大类:一种方案是滑移网格方法。滑移网格技术是常用的动网格方法，这种动网格方法能够处理沿直线平移或者绕着固定轴旋转等网格运动，但是只适用于简单的单个自由度运动，并且运动轨迹固定。在船舶与海洋工程中，滑移网格方法在实际应用中仅限于船舶固定、螺旋桨在船后旋转的情况，无法适应于船体的旋转运动，且一旦船舶拥有较大幅度的...
【技术保护点】
1.一种大幅六自由度运动的流固耦合模拟方法，其特征是：/n(1)划分结构物六自由度网格，以结构物的旋转中心为球心，将计算域划分为包含结构物的球形计算域Ω
【技术特征摘要】
1.一种大幅六自由度运动的流固耦合模拟方法，其特征是：
(1)划分结构物六自由度网格，以结构物的旋转中心为球心，将计算域划分为包含结构物的球形计算域ΩS和球外计算域ΩR，并分别生成计算网格τS和τR；定义τS和τR的球面交界面网格为ΠS和ΠR；
(2)设定参数：分别给定结构物的旋转中心和重心在地球坐标系下的初始坐标和结构物的质量m、基于结构物重心的主转动惯量(Ixcg,Iycg,Izcg)；
(3)读入网格、设置边界条件和初场；
(4)重构球面supermesh；
(5)求解流场；
(6)求解各物理量；
(7)由当前流场和网格节点位置，求解结构物在地球坐标系下的位移(x1,x2,x3,φ,θ,ψ)，其中x1,x2,x3和φ,θ,ψ分别为结构物在地球坐标系下的位移和旋转；
(8)由(x1,x2,x3,φ,θ,ψ)更新结构物重心、旋转中心、网格τS所有节点的坐标；
(9)采用弹簧法，即将τR看作一张具有弹性的网，网格节点之间的联系看作是一根独立的弹簧，在球面交界面网格ΠR运动时，由刚度系数重新分布τR网格节点；返回步骤(4)，开始下一个时间步的计算，直至计算完成。
2.根据权利要求1所述的一种大幅六自由度运动的流固耦合模拟方法，其特征是：重构球面supermesh包括以下步骤：
(4a)分别生成ΠR和ΠS的单元列表和n和m分别为ΠR和ΠS的单元和的总数；
(4b)遍历的单元对遍历的单元由AABB-BOX算法在地球坐标系下判断和是否空间重叠；
(4c)若重叠，采用球心投影，将地球坐标系下和的投射至承影平面得到和
(4d)若不重叠，则返回步骤(4b)；
(4e)由逐边裁剪法对和进行相交计算，得到相交部分
(4f)由逆球心投影将转换为三维笛卡尔坐标系下的Kint；
(4g)在地球坐标系下根据Kint和的面积比值判断Kint是否为有效相交；
(4h)若有效相交，则计算Kint的几何信息，并储拓扑关系；将加入supermesh列表
(4i)若无效相交，则返回步骤(4b)。
3.根据权利要求2所述的一种大幅六自由度运动的流固耦合模拟方法，其特征是：步骤(4b)在地球坐标系下由AABB-BOX算法判断和是否重叠包括以下步骤：
(4ba)分别对和建立直平行六面体和所述ω指对于由N个节点xi＝(xi,yi,zi),i＝1,2,...,N组成的空间平面K，分别取点A＝(max{xi},max{yi},max{zi})和点B＝(min{xi},min{yi},min{zi})，i＝1,2,...,N...
【专利技术属性】
技术研发人员：明平剑，肖艺，张文平，张新玉，柳贡民，曹贻鹏，国杰，赵晓臣，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙;23
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。