MPS方法在船舶大幅晃荡中的应用.pdfVIP

第 28卷第4期 荆楚理工学院学报 2013年 8月 Vo1．28No．4 JournalofJingchuUniversityofTechnology Aug．2013 MPS方法在船舶大幅晃荡中的应用 赵龙华 ，孙新蕾2 (1．上海外高桥造船有限公司生产设计所，上海 200137；2．中南财经政法大学武汉学院信息系，湖北 武汉 430079) 摘要：晃荡是一种非线性 自由面流动现象，剧烈时常伴随 自由面的翻卷和破碎等复杂流动现象。文章将 移动粒子法(MovingParticleSemi—implicit)应用于液体大幅晃荡问题，对二维液舱在共振频率附近时的晃荡 问题进行了研究。计算结果表明：MPS方法能够较好地模拟液体的大幅晃荡，数值得到的 自由面形状和拍击 压力能够很好地与实验和其他数值结果相吻合。激励频率对拍击压力的峰值和作用时间都有较大影响 ，当激 励频率低于共振频率时，第一个压力峰值作用时间很短，第二个压力峰值作用时间较长；当激励频率大于共振 频率时，两个压力峰值作用时间接近。此外，当激励频率接近共振频率时压力峰值较大。 关键词：移动粒子法(MrS)；液舱晃荡；拍击压力；自由面流动 中图分类号：U661．1 文献标志码 ：A 文章编号：1008—4657(2013)04—0o67一o6 0 引言 液舱晃荡是船舶海洋工程中常见的流动现象，带有很强的非线性和随机性。晃荡对液舱壁面产生 的拍击压力可能会造成结构的破坏，并且影响船舶航行的稳定性_】I2J。因此，液舱晃荡成为工程中研究 的热点，对晃荡载荷的预测是船舶航行安全性评估的重要内容之一。 由于晃荡问题十分复杂，常带有 自由面的破碎、翻卷、液体飞溅等现象，因此理论研究晃荡问题十分 困难。早期的研究常基于实验，但实验费用高、周期长，提供的信息有限。随着计算流体力学的发展，数 值模拟成为研究晃荡问题的重要工具。早期的学者基于势流理论求解晃荡问题，这方面比较有代表性 的工作如：Fahinsen发展了一种边界元方法来研究非线性晃荡问题 I4J，wu等采用了有限元 (FEM)对 三维晃荡进行了研究 J，Frandsen发展了一个完全非线性的有限差分(FDM)方法 。虽然势流理论的 研究很大程度地揭示了晃荡问题的特点，但势流方法难以处理大变形的自由面问题，如在剧烈晃荡的自 由面翻卷、破碎等现象。近年来，基于粘流的CFD技术发展十分迅速，常被用来研究复杂的晃荡问题。 由于晃荡是一种典型的自由面流动问题，因此准确预测晃荡流动的一个关键因素就是模拟 自由面变形 ， 为此，人们开发了许多数值方法，如 Hu等发展了CIP(ConstraintInterpolationProfile)方法¨ ，并对二 维、三维剧烈晃荡问题进行了研究，Chen等采用Level—set方法分析了液舱在不同填充率和激励频率时 的晃荡现象 J。在国内，端木玉等基于VOF方法分析了不同液舱结构对晃荡的影响¨，郭晓宇等采用 两相 VOF方法计算了低充水液舱晃荡时的气垫效应_l，刘永涛等利用 VOF方法研究了三维大幅晃 荡 ，并探讨了表面张力的影响。 以上这些 CFD方法都属于网格类方法，计算过程依赖于网格来实现。近年来，一种新颖的计算方 法一无网格粒子法，受到了国内外学者的关注。无网格粒子法的求解过程不需要计算网格，流体域采用 自由移动的粒子来表达，这些粒子具有质量、动量等物理信息，基于拉格朗 日描述法进行追踪。粒子法 求解具有几个优点：可以处理动边界大变形问题；对 自由面的追踪不存在数值耗散。鉴于这些特点，粒 子法从一诞生就受到了很大的关注，并被广泛研究，应用于求解复杂的流动问题。 收稿 日期：2013—06—17 作者简介：赵龙华(1982一)，男，黑龙江尚志人，上海外高桥造船有限公司助理工程师。研究方向：船舶设计； 孙新蕾(1985一)，女，黑龙江绥化人，中南财经政法大学武汉学院讲师，硕士。研究方向：数值分析。 67 目前，已有一些关于粒子法在晃荡问题中应用的工作，如 Delorme等利用 SPH(Smoothe