考虑液舱晃荡影响的flng运动预报方法
【专利摘要】本发明涉及一种考虑液舱晃荡影响的FLNG运动预报方法，包括以下步骤：1)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力；2)将步骤1)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新；3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。与现有技术相比，本发明具有节约计算时间，大大提高了液舱和FLNG的设计效率等优点。
【专利说明】考虑液舱晃荡影响的FLNG运动预报方法
【技术领域】
[0001]本发明涉及一种船舶运动预报，尤其是涉及一种考虑液舱晃荡影响的FLNG(大型浮式液化天然气船)运动预报方法。
【背景技术】
[0002]众所周知，石油资源的短缺以及对能源需求的不断增长，导致海上油气田的大规模开发，同时海洋边际气田的开发也提上了日程。为了经济有效地开发海洋边际油气田，国际上提出了大型浮式液化天然气船(以下简称FLNG)的概念。拥有庞大体积的FLNG长期系泊于作业海域，将经历复杂的海洋环境条件和各种各样的舱内液体装载状况。恶劣的海洋环境条件将引起FLNG的剧烈运动，从而引起舱内液体的晃荡。FLNG的船体运动引起舱内LNG液体晃荡，该晃荡一方面将会对夜舱壁面形成较大的冲击压力，另外一方面反过来也会影响FLNG船体的运动性能。因而，这是一种典型的耦合运动性能。因此，舱内液体晃荡与船体运动耦合响应的预报成为FLNG设计中必须考虑的重要因素之一。
[0003]目前，数值分析是研究舱内液体晃荡对船体运动影响最主要的研究手段之一。通过数值分析，可以较为全面地获得船体运动与舱内液体之间的耦合响应，影响因素以及响应机理等内容，获得较为可靠的计算结果。计算结果可用来校验实验和数值仿真的结果精度。晃荡问题的研究，对我国FLNG液舱设计和FLNG船自主研发，都有重要的意义。舱内液体晃荡问题的难点在于准确模拟液舱内部的自由液面情况。从现有的研究报道可以看出，对于舱内液体晃荡问题的模拟方法主要为仿真计算。仿真计算虽然可以模拟舱内液体的晃荡，但仍难以考虑船体运动的影响。传统的关于晃荡的研究并没有关注与船体运动之间的率禹合影响。
【发明内容】
[0004]本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑液舱晃荡影响的FLNG运动预报方法。
[0005]本发明的目的可以通过以下技术方案来实现:
[0006]一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，包括以下步骤:
[0007]I)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力；
[0008]2)将步骤I)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新；
[0009]3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。
[0010]所述的速度势为一个二维矩形液舱在作周期性微幅横荡的情况下由控制方程和各边界条件以及初始条件所唯一确定的时域解。
[0011]所述的速度势具体计算过程如下:[0012]11)假定舱内流体是无粘性、均匀、不可压缩和无旋的，在全局坐标系中，舱内液体的流动速度势的求解是一个混合初边值问题，其中控制方程为拉普拉斯方程，边界条件为自由面的动力学条件和运动学条件以及液舱壁不可穿透条件，初始条件为液面静止；
[0013]12)将上述方程在全局坐标中的表达转换至随体坐标系中；
[0014]13)将速度势分解后代入定解条件中，从而得到舱内液体非线性晃荡的定解条件在随体坐标系中的表达形式；
[0015]14)利用该定解条件可以求得舱内液体晃荡的速度势:应用面元法对速度势进行求解，每一个单元坐标系的转换关系由一个雅可比矩阵完成，将每个单元内的变量方程及雅克比矩阵代入边界积分方程，取各个节点处的速度势进行方程组的建立并求解，可以得到各节点速度势。
[0016]所述的舱壁的压力计算如下:
[0017]通过引入加速度势，建立关于
【权利要求】
1.一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，包括以下步骤: 1)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力； 2)将步骤I)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新； 3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。
2.根据权利要求1所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的速度势为一个二维矩形液舱在作周期性微幅横荡的情况下由控制方程和各边界条件以及初始条件所唯一确定的时域解。
3.根据权利要求2所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的速度势具体计算过程如下: 11)假定舱内流体是无粘性、均匀、不可压缩和无旋的，在全局坐标系中，舱内液体的流动速度势的求解是一个混合初边值问题，其中控制方程为拉普拉斯方程，边界条件为自由面的动力学条件和运动学条件以及液舱壁不可穿透条件，初始条件为液面静止； 12)将上述方程在全局坐标中的表达转换至随体坐标系中； 13)将速度势分解后 代 入定解条件中，从而得到舱内液体非线性晃荡的定解条件在随体坐标系中的表达形式； 14)利用该定解条件可以求得舱内液体晃荡的速度势:应用面元法对速度势进行求解，每一个单元坐标系的转换关系由一个雅可比矩阵完成，将每个单元内的变量方程及雅克比矩阵代入边界积分方程，取各个节点处的速度势进行方程组的建立并求解，可以得到各节点速度势。
4.根据权利要求1所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的舱壁的压力计算如下: 通过引入加速度势，建立关于#的边值问题并求解，并求出流体对舱壁的压力，其中I?
dt为局部坐标系下舱内流体速度势。
5.根据权利要求4所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的舱壁的压力计算详细过程如下: 通过以下方程，求得当前时刻的压力:
6.根据权利要求5所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的船体运动方程式具体如下:
7.根据权利要求6所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的自由液面的更新具体为: 根据半拉格朗日方法更新边界节点的位置及节点上的变量值:假定自由面上的节点仅做垂向运动，则可以得到节点的速度矢量、自由面上变量的全导数，从而得到自由面上速度势的全导数和进行自由液面、速度势更新的方程:
【文档编号】G06F17/50GK103984793SQ201410108093
【公开日】2014年8月13日 申请日期:2014年3月21日 优先权日:2014年3月21日
【发明者】胡志强, 杨建民, 赵文华, 王舒雅 申请人:上海交通大学