论文部分内容阅读
入水问题作为工程应用中一种比较常见的现象,一直以来都是船舶和海洋工程领域中十分重要的课题。由于涉及到移动物体边界和自由液面,因此能够准确地数值模拟入水过程是具有挑战性的。本文采用自由液面捕捉方法,将自由面作为计算域内的接触间断进行捕捉,建立非均匀不可压缩的欧拉方程作为流场控制方程组。并且利用直角切割网格来离散计算域,通过局部更新物体边界附近的网格信息来解决动边界问题。在此基础上,采用基于单元中心的有限体积法对控制方程进行离散。其中,采用Roe的近似Riemann解算子的复合格式计算流体网格边界处的流通量。边界处的流场变量采用分片线性表达式进行重构,利用最小二乘法求解流场变量梯度,引入限制器来抑制高阶重构方法在间断解附近产生的虚假振动,物体边界处的流通量由Riemann间断理论得出。采用压强分裂算法对压强梯度进行预处理,从而解决由于自由面引起的静压梯度不连续的问题。使用人工压缩法将压强项引入到不可压缩约束方程中,采用双时间推进方法来获得流场的非定常数值解。通过流固耦合算法,将本文方法扩展应用到刚体自由入水问题上。采用二维刚性楔形、圆形剖面常速入水和楔形、船型剖面的自由入水作为实例,验证了本文计算方法的有效性和适用性。采用上述流场计算方法对双圆柱同步入水问题进行了求解,结果表明双圆柱的相互作用对压强分布和所受水作用力的影响较为显著,并且发现双圆柱半径、间距和入水速度对这种相互作用效果有很大的影响。然后,变换这三个状态参量,设计各种不同的双圆柱同步入水数值试验。根据计算结果,在不同入水阶段,对于附加水作用力分别采用类似正态分布和Weibull分布的解析表达式进行数值拟合。并且通过两个任意双圆柱入水模型对估算表达式进行验证,计算所得结果满足工程精度。这说明本文所得的代数表达式能够用于工程中双圆柱入水的附加水作用力估算。这种思想方法同样可以用于多柱状体的入水分析中。将造消波理论引入到本文流场计算方法中建立数值波浪水槽。控制圆柱入水的时间和运动形式来研究不同波浪参数对圆柱入水过程的影响。然后在圆柱运动方程中加入绳索张力的作用,从而实现对悬吊圆柱在波浪中入水过程的数值模拟。最后变化计算模型的圆柱运动和波浪参数,分析不同工况下绳索张力变化和重物运动的特性。将刚架结构有限元理论与本文流场计算方法相结合,并采用单向耦合策略联合计算流场和弹性变形,从而将本文方法扩展应用到水弹性问题的求解中。建立弹性楔形体的自由入水模型,将本文结果与实验数据进行比较,证明了本文方法求解弹性楔形体入水问题的正确性和合理性。然后改变楔形体的入水参数建立不同的模型,分别进行数值计算,并且分析了入水参数对自由入水的弹性楔形体的整体运动和局部变形的影响。最后针对三维船体从船台上纵向重力式下水这一工程实际问题,采用Fluent软件计算船体下水流场,并在此基础上编制的船体下水三维运动方程的VC源代码来进行二次开发,从而达到动态模拟船体下水过程的目的。对渤海造船厂159000t油船下水实例进行数值计算,结果表明采用本文的下水数值方法分析船体下水过程是准确可行的。然后以大连重工集团的298000t油船首部从H3倾斜平台下水为例,来分析船体和船台各状态参数对下水运动的影响,得到一些有用结论用来指导实船下水工艺。