【技术实现步骤摘要】
一种基于全局形函数的系泊系统求解方法
[0001]本专利技术涉及浮式结构的系泊系统
，具体涉及一种基于全局形函数的系泊系统求解方法，适用于悬链线式系泊系统。
技术介绍
[0002]这里的陈述仅提供与本公开相关的
技术介绍
，而不必然地构成现有技术。
[0003]浮式平台是海上作业的重要设备，是走向深远海的重要依托，其中浮式平台的系泊系统是保障其在作业过程中的安全，因此对系泊系统进行动力响应分析至关重要。系泊系统是长期服役在复杂海况下的水下大柔性结构，具有很强的非线性行为，目前常见的系泊系统分析软件如SESAM采用的是有限元方法， Orcaflex采用的是集中质量法，NREL的FAST则提供了准静态系泊模型、集中质量法和弹性杆理论三种求解方法。然而浮式平台的振动特征和运动特征会给系泊系统带来复杂的非线性边界效应，系泊系统自身的大柔性和大变形增加了动力分析的复杂程度，软件分析的整体耦合分析中，系泊系统的求解占用大部分计算资源，方法采用Newmark
‑
beta法和牛顿法对平衡方程进行数值求解，在整体耦合分析中，占用大部分的计算资源，因此系泊系统动力分析领域需要一种简化数值求解过程、提高分析和计算效率的方法，来有效节省初步设计的时间成本。
技术实现思路
[0004]本专利技术针对考虑系泊系统大变形和非线性的计算效率问题，提供一种基于全局形函数的全新系泊系统求解方法，该方法可以快速求解系泊系统动力响应，适用于系泊系统的响应评估。
[0005]本专利技术的技术方案为：>[0006]本专利技术提供了一种基于全局形函数的系泊系统求解方法，包括以下步骤：
[0007]S1：根据系泊系统的设计参数，将系泊系统分为系泊缆和锚；系泊缆在海水中静平衡状态存在悬浮段和躺底段，采用悬链线方程求解系泊缆悬浮段的初始位置和躺底段与土体的静平衡位置，并将锚简化为弹性支座，建立系泊系统初始状态模型；
[0008]S2：步骤S1中的静力模型求解，为动力模型提供迭代的初始条件，现进行动力模型构建：将系泊缆离散为杆单元模型，构建基于哈密顿原理的能量平衡方程，建立系泊缆的弹性大变形方程，并引入泰勒展开，略去高阶项，构建基于能量平衡方程的变形因子，从而建立考虑系泊系统大变形的刚度矩阵；
[0009]S3：在动力分析中，需要施加边界条件进行约束，以下是边界条件模型：建立系泊缆与海水相互作用体系：针对小尺度柔性结构的波浪力模型，引入莫里森方程，并针对躺底段构建弹性海床模型，通过变分法原理，得出系泊缆在受非线性边界条件、水动力荷载和流荷载作用下的动力学模型；建立系泊缆的非线性边界和变形模型：将海床模拟为线弹性地基，并建立系泊缆与海床摩擦模型；
[0010]S4：由步骤S2和步骤S3建立的考虑系泊缆非线性边界和变形的系泊系统控制方
程，引入状态空间模型，对控制方程进行降阶处理，最后采用数值求解方法对控制方程进行求解。
[0011]进一步地，步骤S1具体实现如下：
[0012]采用悬链线方程进行求解系泊缆静力部分，过程如下：
[0013]根据牛顿第二定律和触底点的边界条件，可得：
[0014][0015]y|
x＝0
＝0
ꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0016][0017]dV＝wds
ꢀꢀꢀꢀꢀꢀꢀ
(4)
[0018]其中，S为系泊缆原长；
[0019]s为其悬浮段长度；
[0020]p为悬浮段水平长度；
[0021]q为躺底段水平长度；
[0022]X为锚泊点与导览孔的水平距离；
[0023]d为水深；
[0024]T为导缆孔处系泊缆张力；
[0025]V为T的垂向分量；
[0026]H为T的水平分量；
[0027]θ为系泊缆单元与水平方向的夹角；
[0028]由悬链线理论可知，导缆孔处系泊缆张力的水平分量H为常数，并由近似条件可知：根据式(1)～(4)求解方程，可得悬浮段方程为：
[0029][0030]其中，
[0031]w为系泊缆单位长度的湿重；
[0032]由躺底段稳定条件可知：
[0033]y＝0(
‑
q<x<0)
ꢀꢀꢀꢀꢀꢀ
(6)
[0034]将式(5)、(6)整合，可得完整的悬链线方程：
[0035][0036]上式中存在三个未知量，a、p和q，同时也存在三个约束方程，分别是：
[0037]y|
x＝p
＝d
ꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0038][0039]p+q＝X
ꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0040]根据上述三个约束方程求得悬链线的超越方程：
[0041][0042]对超越方程进行迭代求解即可以上求出未知量a、p和q。
[0043]进一步地，步骤S2具体实现如下：
[0044]假设系泊缆满足线弹性变形，则其形变方程可以表示如下：
[0045][0046][0047]其中，ds和ds'分别表示单元变形前后的长度；
[0048]dx、dy和dz是单元长度在x、y和z方向的分量；
[0049]i和j表示单元两端的节点数；
[0050]u、v和w表示i节点在x、y和z方向的位移；
[0051]u+du、v+dv和w+dw表示j节点在x、y和z方向的位移；
[0052][0053]其中对形变ε求关于u'、v'和w'的偏微分，可得：
[0054][0055]因此，形变的全微分是：
[0056]δε＝(x'+u')δu'+(y'+v')δv'+(z'+w')δw'
ꢀꢀꢀꢀ
(15)
[0057]由最小作用量原理可知：
[0058][0059]其中，T、V和W分别是系统的动能、势能和非保守力虚功；
[0060]δ是变分符号；
[0061]对动能、势能和非保守力虚功分别进行变分；
[0062][0063][0064][0065]其中，ρ为系泊缆密度，
[0066]A为系泊缆横截面积；
[0067]..表示对时间求二阶导数；
[0068]U为系泊缆的位移；
[0069]T为系泊缆静态张力；
[0070]E为杨氏模量；
[0071]g为重力加速度；
[0072]f
c
为阻尼力；
[0073]将式(17)～(19)代入式(16)，可得：
[0074][0075]采用有限元方法将平衡方程进行离散：
[0076]对于单元e，令U
e
＝{u
e
,v
e
,w
e
}
T
＝N
e
d
e
和d
e
...
【技术保护点】
【技术特征摘要】
1.一种基于全局形函数的系泊系统求解方法，其特征在于，包括以下步骤：S1：根据系泊系统的设计参数，将系泊系统分为系泊缆和锚；系泊缆在海水中静平衡状态存在悬浮段和躺底段，采用悬链线方程求解系泊缆悬浮段的初始位置和躺底段与土体的静平衡位置，并将锚简化为弹性支座，建立系泊系统初始状态模型；S2：步骤S1中的静力模型求解，为动力模型提供迭代的初始条件，现进行动力模型构建：将系泊缆离散为杆单元模型，构建基于哈密顿原理的能量平衡方程，建立系泊缆的弹性大变形方程，并引入泰勒展开，略去高阶项，构建基于能量平衡方程的变形因子，从而建立考虑系泊系统大变形的刚度矩阵；S3：在动力分析中，需要施加边界条件进行约束，以下是边界条件模型：建立系泊缆与海水相互作用体系：针对小尺度柔性结构的波浪力模型，引入莫里森方程，并针对躺底段构建弹性海床模型，通过变分法原理，得出系泊缆在受非线性边界条件、水动力荷载和流荷载作用下的动力学模型；建立系泊缆的非线性边界和变形模型：将海床模拟为线弹性地基，并建立系泊缆与海床摩擦模型；S4：由步骤S2和步骤S3建立的考虑系泊缆非线性边界和变形的系泊系统控制方程，引入状态空间模型，对控制方程进行降阶处理，最后采用数值求解方法对控制方程进行求解。2.根据权利要求1所述的一种基于全局形函数的系泊系统求解方法，其特征在于，步骤S1具体实现如下：采用悬链线方程进行求解系泊缆静力部分，过程如下：根据牛顿第二定律和触底点的边界条件，可得：y|
x＝0
＝0
ꢀꢀꢀꢀꢀꢀ
(2)dV＝wds
ꢀꢀꢀꢀꢀꢀ
(4)其中，S为系泊缆总长；s为其悬浮段长度；p为悬浮段水平长度；q为躺底段水平长度；X为锚泊点与导览孔的水平距离；d为水深；T为导缆孔处系泊缆张力；V为T的垂向分量；H为T的水平分量；θ为系泊缆单元与水平方向的夹角；由悬链线理论可知，导缆孔处系泊缆张力的水平分量H为常数，并由近似条件可知：根据式(1)～(4)求解方程，可得悬浮段方程为：
其中，w为系泊缆单位长度的湿重；由躺底段稳定条件可知：y＝0(
‑
q<x<0)
ꢀꢀꢀꢀꢀꢀ
(6)将式(5)、(6)整合，可得完整的悬链线方程：式(7)中存在三个未知量，a、p和q，同时也存在三个约束方程，分别是：y|
x＝p
＝d
ꢀꢀꢀꢀꢀꢀꢀ
(8)p+q＝X
ꢀꢀꢀꢀꢀꢀꢀ
(10)根据上述三个约束方程求得悬链线的超越方程：对超越方程进行迭代求解即可以上求出未知量a、p和q。3.根据权利要求1所述的一种基于全局形函数的系泊系统求解方法，其特征在于步骤S2具体实现如下：假设系泊缆满足线弹性变形，则其形变方程可以表示如下：假设系泊缆满足线弹性变形，则其形变方程可以表示如下：其中，ds和ds'分别表示单元变形前后的长度；dx、dy和dz是单元长度在x、y和z方向的分量；i和j表示单元两端的节点数；u、v和w表示i节点在x、y和z方向的位移；u+du、v+dv和w+dw表示j节点在x、y和z方向的位移；其中对形变ε求关于u'、v'和w'的偏微分，可得：因此，形变的全微分是：δε＝(x'+u')δu'+(y'+v')δv'+(z'+w')δw'
ꢀꢀꢀꢀ
(15)由最小作用量原理可知：
其中，T、V和W分别...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。