考虑液舱晃荡影响的FLNG运动预报方法技术

本发明专利技术涉及一种考虑液舱晃荡影响的FLNG运动预报方法，包括以下步骤：1)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力；2)将步骤1)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新；3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。与现有技术相比，本发明专利技术具有节约计算时间，大大提高了液舱和FLNG的设计效率等优点。

【技术实现步骤摘要】
考虑液舱晃荡影响的FLNG运动预报方法
本专利技术涉及一种船舶运动预报，尤其是涉及一种考虑液舱晃荡影响的FLNG(大型浮式液化天然气船)运动预报方法。
技术介绍
众所周知，石油资源的短缺以及对能源需求的不断增长，导致海上油气田的大规模开发，同时海洋边际气田的开发也提上了日程。为了经济有效地开发海洋边际油气田，国际上提出了大型浮式液化天然气船(以下简称FLNG)的概念。拥有庞大体积的FLNG长期系泊于作业海域，将经历复杂的海洋环境条件和各种各样的舱内液体装载状况。恶劣的海洋环境条件将引起FLNG的剧烈运动，从而引起舱内液体的晃荡。FLNG的船体运动引起舱内LNG液体晃荡，该晃荡一方面将会对夜舱壁面形成较大的冲击压力，另外一方面反过来也会影响FLNG船体的运动性能。因而，这是一种典型的耦合运动性能。因此，舱内液体晃荡与船体运动耦合响应的预报成为FLNG设计中必须考虑的重要因素之一。目前，数值分析是研究舱内液体晃荡对船体运动影响最主要的研究手段之一。通过数值分析，可以较为全面地获得船体运动与舱内液体之间的耦合响应，影响因素以及响应机理等内容，获得较为可靠的计算结果。计算结果可用来校验实验和数值仿真的结果精度。晃荡问题的研究，对我国FLNG液舱设计和FLNG船自主研发，都有重要的意义。舱内液体晃荡问题的难点在于准确模拟液舱内部的自由液面情况。从现有的研究报道可以看出，对于舱内液体晃荡问题的模拟方法主要为仿真计算。仿真计算虽然可以模拟舱内液体的晃荡，但仍难以考虑船体运动的影响。传统的关于晃荡的研究并没有关注与船体运动之间的耦合影响。
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑液舱晃荡影响的FLNG运动预报方法。本专利技术的目的可以通过以下技术方案来实现：一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，包括以下步骤：1)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力；2)将步骤1)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新；3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。所述的速度势为一个二维矩形液舱在作周期性微幅横荡的情况下由控制方程和各边界条件以及初始条件所唯一确定的时域解。所述的速度势具体计算过程如下：11)假定舱内流体是无粘性、均匀、不可压缩和无旋的，在全局坐标系中，舱内液体的流动速度势的求解是一个混合初边值问题，其中控制方程为拉普拉斯方程，边界条件为自由面的动力学条件和运动学条件以及液舱壁不可穿透条件，初始条件为液面静止；12)将上述方程在全局坐标中的表达转换至随体坐标系中；13)将速度势分解后代入定解条件中，从而得到舱内液体非线性晃荡的定解条件在随体坐标系中的表达形式；14)利用该定解条件可以求得舱内液体晃荡的速度势：应用面元法对速度势进行求解，每一个单元坐标系的转换关系由一个雅可比矩阵完成，将每个单元内的变量方程及雅克比矩阵代入边界积分方程，取各个节点处的速度势进行方程组的建立并求解，可以得到各节点速度势。所述的舱壁的压力计算如下：通过引入加速度势，建立关于的边值问题并求解，并求出流体对舱壁的压力，其中为局部坐标系下舱内流体速度势。所述的舱壁的压力计算详细过程如下：通过以下方程，求得当前时刻的压力：p为所求的舱壁压力，ρ为液体密度，为舱内液体速度势，x和z为液体在x和z轴方向上的坐标，u和ω分别代表FLNG船体在x和z轴方向上运动的速度分量，t为时间变量，η为自由面高度，g为重力加速度；在上述方程中，除了其余所有参数均在步骤1)的定解问题的求解过程中得到，为了求解该未知项，引入加速度势概念，关于该未知项的边值问题表示为如下公式：在流域内；在自由液面上；在液舱壁面上；其中，ax和az分别代表FLNG船体加速度在x和z方向上的分量，为FLNG船体的速度矢量，和分别代表液舱壁面的单位法向量和单位切向量；解得以上定解问题后，即可根据压力公式求得流场内的压力分布，沿壁面压力积分，得到舱壁上的压力。所述的船体运动方程式具体如下：其中mjk,μjk和cjk分别代表浮体的质量、频率无穷大时的附加质量以及静水恢复力系数；k和j代表船体运动模态；ξk为船体位移向；Fj(t)代表作用于FLNG船体上的力，包括外部波浪力以及内部液体晃荡引起的作用力；hjk(t-τ)代表延时函数，表征流体的记忆效应，τ为参数。所述的自由液面的更新具体为：根据半拉格朗日方法更新边界节点的位置及节点上的变量值：假定自由面上的节点仅做垂向运动，则可以得到节点的速度矢量、自由面上变量的全导数，从而得到自由面上速度势的全导数和进行自由液面、速度势更新的方程：为自由面液体速度势，x和z为液体在x和z轴方向上的坐标，u和ω分别代表FLNG船体在x和z轴方向上运动的速度分量，t为时间变量，η为自由面高度，g为重力加速度；通过这两个方程对自由液面处的速度势及波面升高进行更新。与现有技术相比，本专利技术具有以下优点：(1)在时域耦合计算模型中，采用半拉格朗日方法对舱内液体晃荡时的非线性自由液面进行捕捉，使得自由面模型更加精确；(2)基于加速度势概念，建立了关于加速度势的初边值问题，从而更为精确地报舱内液体晃荡对舱壁造成的冲击力；(3)本方法比CFD方法更为节约计算时间，大大提高了液舱和FLNG的设计效率。附图说明图1是液舱晃荡问题的计算域示意图；其中LW为液舱壁，LF为自由表面，X、Z分别为水平及垂直坐标。图2是二维船体运动计算域示意图；其中LB为船体湿表面，LL、LR分别为左右侧边界，LF为自由液面，H为水深，D+、D-为计算域，X、Z分别为水平及垂直坐标。图3是液舱与船体运动耦合示意图。具体实施方式下面结合附图和具体实施例对本专利技术进行详细说明。实施例以下实施例在以本专利技术技术方案为前提下进行实施，给出了详细的实施方式和过程，但本专利技术的保护范围不限于下述的实施例。本专利技术提供一种船舶运动预报方法，其步骤包括：第一、求一个矩形液舱在一个时间步内舱内液体流动的速度势。第二、引入加速度势，建立了关于加速度势的边值问题，求得流体对舱壁的压力。第三，建立船体运动的边值问题，求得速度势，波浪力，阻尼系数和延迟函数，建立运动方程。第四、求解运动方程，获得船舶位置和速度。由采用半拉格朗日方法进行数值模拟，以表征一个时间步后液舱内自由液面的变化。船舶的位置和液舱自由面情况将作为第一步中求解速度势的液舱边界条件。对于步骤一：假定舱内流体是无粘性、均匀、不可压缩和无旋的，在全局坐标系中，舱内液体的流动速度势的求解是一个混合初边值问题：控制方程为拉普拉斯方程，边界条件为自由面的动力学条件和运动学条件以及液舱壁不可穿透条件，初始条件为液面静止。在实际的数值计算中，需要将上述方程在全局坐标中的表达转换至随体坐标系中。将速度势分解后代入定解条件中，从而得到舱内液体非线性晃荡的定解条件在随体坐标系中的表达形式。利用该定解条件可以求得舱内液体晃荡的速度势。应用面元法对速度势进行求解。每一个单元坐标系的转换关本文档来自技高网...

【技术保护点】
一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，包括以下步骤： 1)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力； 2)将步骤1)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新； 3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。

【技术特征摘要】
1.一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，包括以下步骤：1)基于势流理论和面元法，在时域内对舱内液体的非线性晃荡进行描述，求得舱内液体晃荡的速度势，和晃荡所引起的液体对舱壁的压力；2)将步骤1)得到的压力作为船体运动方程式右端的外力输入并考虑外部波浪力的影响，计算船体的运动响应，求得船体的位置，同时进行自由液面的更新；3)将得到的当前时刻船体位置和更新后的自由液面作为舱内液体晃荡问题的边界条件，从而求得下一时刻舱内液体的晃荡问题，如此不断循环，至设定时间后即可停止。2.根据权利要求1所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的速度势为一个二维矩形液舱在作周期性微幅横荡的情况下由控制方程和各边界条件以及初始条件所唯一确定的时域解。3.根据权利要求2所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的速度势具体计算过程如下：11)假定舱内流体是无粘性、均匀、不可压缩和无旋的，在全局坐标系中，舱内液体的流动速度势的求解是一个混合初边值问题，其中控制方程为拉普拉斯方程，边界条件为自由面的动力学条件和运动学条件以及液舱壁不可穿透条件，初始条件为液面静止；12)将上述方程在全局坐标中的表达转换至随体坐标系中；13)将速度势分解后代入定解条件中，从而得到舱内液体非线性晃荡的定解条件在随体坐标系中的表达形式；14)利用该定解条件可以求得舱内液体晃荡的速度势：应用面元法对速度势进行求解，每一个单元坐标系的转换关系由一个雅可比矩阵完成，将每个单元内的变量方程及雅克比矩阵代入边界积分方程，取各个节点处的速度势进行方程组的建立并求解，可以得到各节点速度势。4.根据权利要求2所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的舱壁的压力计算如下：通过引入加速度势，建立关于的边值问题并求解，并求出流体对舱壁的压力，其中为局部坐标系下舱内流体速度势。5.根据权利要求4所述的一种考虑液舱晃荡影响的FLNG运动预报方法，其特征在于，所述的舱壁的压力计算详细过程如下：通过以下方程，求得当前时刻的压力：p为所求的舱壁压力，ρ为液体...

【专利技术属性】
技术研发人员：胡志强，杨建民，赵文华，王舒雅，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：上海;31