载液船舱流体晃荡时域模拟的分析研究.pdfVIP

A Dissertation Submitted to Shanghai Jiao Tong University for Master Degree The simulation and analysis of ship tank sloshing in time domain Major ：Design and Manufacture of Naval Architecture and Ocean Engineering Candidate : He Zhao Supervisor: Zhu Renchuan School of Naval Architecture, Ocean and Civil Engineering Shanghai Jiao Tong University February 2010 2 载液船舱流体晃荡时域模拟研究 摘 要 液舱内流体晃荡问题是水动力学比较经典的问题之一，一直倍受 学术界和工程界的关注。随着 LNG 和 LPG 等特种液货船型的研制和 广泛应用,液体晃荡问题的研究已经成为船舶水动力学中的一个热门 课题。针对于船舶液舱内流体晃荡问题的试验研究日本学者开展得较 多，如 Katsuya Umemoto 和 Takao Yoshikawa 对不同结构的液舱进行 过一系列的实验[1] ；由于存在强烈的非线性自由面现象，理论研究方 面进展并不大，理论研究者也尝试采用了多种理论方法，如边界元法， 有限元法，以及近年发展较快的计算流体方法等。有限元法和计算流 体方法占用大量的计算机资源，计算速度较慢。应用边界元法求解水 动力学问题，可以对研究问题进行降低一维来处理，从而对计算机资 源要求相对较低，计算速度快，常被用来模拟计算伴有自由面现象的 水动力学问题，如 Faltinsen 等对没有隔板晃荡问题开展的理论研究 [2] [3] 并进行了试验验证 , 朱仁传等 采用了多区域边界元方法对有隔板 的时域晃荡问题进行过研究。实际上载液船舱都装有能渗透的隔板， 缪国平等[4, 5]试图对有渗透隔板的液舱内流体晃荡现象进行理论研 究，但是采用的方法是线性频域的，并不能很好地描述流体晃荡问题 的非线性现象。 液舱内的隔板相对较薄，采用边界元方法求解带隔板的液体晃荡 3 问题时，其厚度一般忽略不计，在隔板上布置源点时，会出现有些场 点与源点相距太近，给格林函数带来了奇性，也会使数值计算产生误 差，虽然 Nishino 等通过对边界元公式进行特殊处理[6]，并且引进了 薄板边界元方法，但是该方法处理带有可渗透隔板的情形存在一定的 困难。 针对以上情况，本文提出了利用多域边界元方法对二维和三维带 有可渗透隔板的液舱内流体晃荡的非线性时域问题进行计算模拟，所 得结果与已发表论文中算例模拟结果吻合良好。通常，我们认为二维 的液舱模型等价于宽度趋向于无穷的三维液舱，因为它忽略了沿宽度 方向舱壁对液舱内液体晃动的影响。本文讨论了液舱模型的三维效 应，为其他学者提供了理论参考。另外，我们模拟了不同宽度，隔板 高度和渗透率的液舱模型的晃动现象并且讨论了他们对液舱晃动固 有频率的影响。该方法可以进一步应用于载液船舱、减摇水舱等设置 不同渗透率隔板的水动力问题的分析，在船舶的减摇、减振等工程应 用方面也具有重要意义。 另外本文对液舱内的晃荡和船舶之间的耦合运动也做了初步的 探讨。液舱内液体的晃荡由船舶运动的激励而产生，反过来由于液体 晃荡对舱壁的冲击产生的力又影响到船舶的运动，因此考虑液舱晃荡 和船舶运动的耦合作用在船舶设计或者运动控制中起着至关重要的 作用。