发明名称 一种钢悬链式立管纵向运动的分析方法
摘要 本发明公开了一种钢悬链式立管纵向运动的分析方法。该方法对现有钢悬链式立管纵向运动分析技术进行了改进,在考虑了钢悬链式立管横向运动的刚体摆动的基础上建立了钢悬链式立管纵向运动的分析模型。由于增加了钢悬链式立管横向运动的刚体摆动对纵向弯曲运动的影响,使钢悬链式立管的纵向运动分析更加符合实际的结构状态。
申请公布号 CN102708281A 申请公布日期 2012.10.03
申请号 CN201210107701.6 申请日期 2012.04.10
申请人 中国海洋大学 发明人 黄维平;白兴兰;刘娟
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 代理人
主权项 1.一种钢悬链式立管纵向运动的分析方法,其特征在于:建立的钢悬链式立管纵向运动方程如下:公式1:<img file="259136DEST_PATH_FSB00000841291000011.GIF" wi="919" he="109" />公式2:<img file="176277DEST_PATH_FSB00000841291000012.GIF" wi="1217" he="134" />在所述公式1和公式2中:r<sub>b</sub>——弯曲位移;<img file="956014DEST_PATH_FSB00000841291000013.GIF" wi="32" he="49" />——弯曲速度;<img file="351223DEST_PATH_FSB00000841291000014.GIF" wi="33" he="51" />——弯曲加速度;<img file="165595DEST_PATH_FSB00000841291000015.GIF" wi="32" he="50" />——刚体摆动速度,<img file="988058DEST_PATH_FSB00000841291000016.GIF" wi="216" he="92" /><img file="255091DEST_PATH_FSB00000841291000017.GIF" wi="34" he="49" />——刚体摆动加速度,<img file="632153DEST_PATH_FSB00000841291000018.GIF" wi="235" he="85" />α<sub>r</sub>——微元段转动的角位移;<img file="301032DEST_PATH_FSB00000841291000019.GIF" wi="45" he="52" />——微元段转动的角速度;<img file="294396DEST_PATH_FSB000008412910000110.GIF" wi="46" he="51" />——微元段转动的角加速度;s——微元段距转动轴的矢径;s——微元段距转动轴矢径的模;s<sub>1</sub>——微元段距转动轴矢径的模在x轴的投影;s<sub>2</sub>——微元段距转动轴矢径的模在y轴的投影;s<sub>3</sub>——微元段距转动轴矢径的模在z轴的投影;m——钢悬链式立管单位长度的质量;m<sub>a</sub>——钢悬链式立管单位长度的附加质量;c——钢悬链式立管单位长度的结构阻尼系数;c<sub>a</sub>——钢悬链式立管单位长度的附加阻尼系数; k——立管弯曲刚度,包括弹性弯曲刚度和几何刚度;q——作用在钢悬链式立管上的海洋环境荷载;g——重力加速度;基于所述钢悬链式立管纵向运动方程,按以下步骤计算运动响应:1)设定笛卡儿坐标系,设定钢悬链式立管的顶点为坐标原点,设坐标系的y-z平面为悬链线平面;2)根据钢悬链式立管的悬挂点和触地点坐标,按如下公式3计算摆动轴的单位矢量:公式3:<img file="48725DEST_PATH_FSB00000841291000021.GIF" wi="404" he="62" />式中:<img file="785737DEST_PATH_FSB00000841291000022.GIF" wi="34" he="40" />——钢悬链式立管的悬挂点与触地点连线形成的摆动轴的单位矢量<img file="309122DEST_PATH_FSB00000841291000023.GIF" wi="319" he="118" /><img file="738966DEST_PATH_FSB00000841291000024.GIF" wi="305" he="148" /><img file="918275DEST_PATH_FSB00000841291000025.GIF" wi="376" he="128" /><img file="193398DEST_PATH_FSB00000841291000026.GIF" wi="830" he="125" />其中:x<sub>D</sub>,X<sub>O</sub>,y<sub>D</sub>,y<sub>O</sub>,z<sub>D</sub>,z<sub>O</sub>分别为钢悬链式立管触地点和悬挂点的坐标值,在平衡位置时,c<sub>1</sub>=0;3)按如下公式4计算矢量s:公式4:<img file="836869DEST_PATH_FSB00000841291000027.GIF" wi="358" he="84" />式中:s<sub>1</sub>=(x<sub>B</sub>-x<sub>A</sub>),s<sub>2</sub>=(y<sub>B</sub>-y<sub>A</sub>),s<sub>3</sub>=(z<sub>B</sub>-z<sub>A</sub>)在平衡位置时,s<sub>1</sub>=0;4)将所述公式1表示为坐标分量的形式的如下公式5、公式6和公式7:公式5:<img file="172036DEST_PATH_FSB00000841291000028.GIF" wi="984" he="92" />公式6:<img file="900957DEST_PATH_FSB00000841291000029.GIF" wi="993" he="120" />公式7:<img file="979772DEST_PATH_FSB000008412910000210.GIF" wi="1108" he="133" />5)用有限元方法将所述公式5、公式6、公式7和公式2离散后即可得到如下公式8和公式9所示的钢悬链式立管纵向运动的有限元方程:公式8:<img file="477749DEST_PATH_FSB00000841291000031.GIF" wi="652" he="92" />公式9:<img file="249396DEST_PATH_FSB00000841291000032.GIF" wi="724" he="117" />在所述公式8和公式9中:[M]——系统的质量矩阵;[C]——系统的阻尼矩阵;[K]——系统的刚度矩阵;<img file="386985DEST_PATH_FSB00000841291000033.GIF" wi="62" he="65" />——系统的加速度向量;<img file="269491DEST_PATH_FSB00000841291000034.GIF" wi="64" he="63" />——系统的速度向量;{r}——系统的位移向量;{F}——系统的荷载向量;[I]——系统的转动惯量矩阵;[C<sub>α</sub>]——系统的水动力阻尼矩阵;[K<sub>α</sub>]——系统的刚体摆动回复力系数矩阵;<img file="621975DEST_PATH_FSB00000841291000035.GIF" wi="73" he="65" />——系统的刚体摆动角加速度向量;<img file="564523DEST_PATH_FSB00000841291000036.GIF" wi="73" he="65" />——系统的刚体摆动角速度向量;{α}——系统的刚体摆动角位移向量;{M<sub>α</sub>}——系统的刚体摆动外力矩向量;6)采用时程积分法求解公式8和公式9即可得到钢悬链式立管的纵向运动响应;所述时程积分法包括Newmark-β法或Wilson-θ法。 
地址 266000 山东省青岛市崂山区松岭路238号