收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

柔性立管涡激振动时域响应分析

黄智勇  
【摘要】: 地球的表面约70%为海洋,蕴藏在海洋里的资源将成为人类赖以生存的重要保障。就目前国际工业发展趋势而言,对能源的需求量,特别是石油的需求量日益增大,深海石油的开采已成为解决石油能源危机的主要途径。与深海石油工业密切相关的一个国际性的难题便是输油立管的涡激振动问题。 立管的涡激振动是指洋流流经圆柱状细长柔性立管时,由于压力的不均将产生交替脱落的漩涡,由泻涡诱发的非定常流体力会引起立管的振动响应。在一定来流条件下,泻涡将使立管产生强烈的振动,由于泻涡的周期较短,涡激振动成为立管疲劳损伤导致结构破坏的主要因素。涡激振动是一个强非线性的的流固耦合问题,特别是涉及到高雷诺数条件下湍流的研究。 本文作为国家自然科学基金专项基金项目“深海平台的动力特性研究”(项目编号:50323004)和上海市专项基金项目(项目编号:05DJ14001)的组成部分,旨在通过数值计算方法来研究湍流,并在湍流的基础上实现流固耦合的计算,即实现深海立管涡激振动的数值模拟,建立深海立管涡激振动时域预报模型。为了更接近深海立管实时作业环境,数值计算实现了均匀来流和剪切来流条件下细长柔性立管的时域响应预报。 本文参考了国际上有关湍流研究和细长柔性立管涡激振动的最新研究成果,对研究湍流理论的主要数值计算方法进行了比较和讨论,并对细长柔性立管涡激振动的热点问题展开了深入的研究,其主要内容可以概括如下: 1总结了深海立管涡激振动响应的频域和时域预报模型,对深海立管涡激振动响应的实验进行了归纳,指出现有涡激振动预报模型的缺陷和实验研究的主要方向。对研究湍流理论的三种主要数值计算方法进行综合评述,阐述了湍流的本质以及各湍流研究数值计算方法的优缺点,着重对直接数值模拟方法和雷诺平均法计算湍流的主要思想和基本理论进行了详细的介绍。 2采用直接数值模拟对低雷诺数条件下圆柱绕流进行了研究,其中来流条件分别为对均匀来流、横向剪切流和平面剪切流来流。数值计算成功再现了实验中观测到的三维尾涡结构,这在国内数值研究报道中纯属首次。分析和讨论了三维圆柱绕流的物理特性诸如流体动能、雷诺应力、泻涡频率、边界层以及圆柱体表面的压力系数,探究了剪切流对流体流动和泻涡模式的影响。基于直接数值模拟对三维流体流动计算的优越性,用数值实验的方法成功解析了物理实验结果与二维数值计算结果之间存在差异的原因。 3采用雷诺平均方法结合湍流模式对低质量比弹性支撑圆柱体的涡激振动进行了计算。数值模拟考虑了圆柱体横向和流向的振动对横向振幅幅值的影响,并与限制流向运动时的结果进行了比较。数值计算找出了横向和流向自由度条件下,产生较大横向振幅幅值时所对应的圆柱体的临界质量比。本文的数值结果成功呈现了物理实验中观测到的各种尾涡模式,并得出了均匀来流条件下实验中常见的圆柱体的运动轨迹。 4在横向剪切来流条件下,对两向自由度弹性支撑圆柱体的涡激振动进行了计算。数值模拟得到的横向振幅、运动轨迹、升力系数和位移时历曲线的结果表明,对于不同的雷诺数,存在一个剪切率锁定区间。在横向来流剪切率较小时,圆柱体的运动轨迹与均匀来流条件下的运动轨迹一致。随着剪切率的增加,圆柱体的运动形式表现出一种全新的轨迹。 5基于切片理论并结合有限体积法和动力有限元方法,成功建立了深海立管涡激振动响应的时域计算模型。对阶梯状来流以及轴向剪切流条件下深海立管涡激振动响应进行了时域计算,分析了立管在不同来流条件下所激发的振型和模态。本时域预报模型还记录了沿立管轴向分布的不同尾涡结构,重现了物理实验中所观测到的一系列流体流动现象。 6尽管在阶梯状来流条件下所预报的高阶模态响应与实验结果之间还存在一定差异,但本文的时域预报模型在轴向剪切流条件下所得到的低阶模态响应与实验结果吻合较好,且比同类时域模型和频率模型预报的结果更为精确。 7考虑雷诺数和质量阻尼比对自激振动系统横向振幅幅值的影响,并应用受迫振动的实验数据,从理论上建立了弹性支撑圆柱体涡激振动振幅幅值的预报模型。该理论模型所预报的结果能与现有大部分物理实验结果相吻合,同时可为细长立管频率预报模型考虑雷诺数和质量阻尼比的影响提供参考。 以上对本文研究的主要内容与进行了概述,本论文的结论主要有以下几个方面: 1首次用直接数值模拟方法对各种流条件下的圆柱饶流进行了计算。在低雷诺数条件下,成功再现了均匀来流下的三维尾涡模式;如Mode A和Mode B。在横向剪切来流条件下,与来流速度较低一侧对应的尾涡能量在流体运动过程中逐渐衰减;与来流速度较高一侧对应的尾涡能量在流体运动过程中逐渐增强,以至于在尾流较远处只剩下能量较强的尾涡。横向剪切流对圆柱绕流的影响还体现在边界层和表面压强系数的改变。数值计算结果表明,圆柱体表面的边界层和压强系数在横向剪切流作用下沿轴向出现偏转,这与流体流动的物理特性相一致。而在轴向剪切流条件下,泻涡形式沿轴向倾斜,顺时针和逆时针旋转的尾涡不存在能量的衰减或增强,脱落后的尾涡在流向方向依然保持平行。在横向-轴向组合剪切来流条件下,尾涡的运动不仅偏向于来流速度较低的一侧,而且沿轴向倾斜脱落。组合剪切流可看做是单个横向剪切流和轴向剪切流作功的叠加,但不会出现横向剪切流和轴向剪切流之间相互抑制的现象。 2采用直接数值模拟法研究了无滑移边界条件和自由表面边界条件对尾涡的影响,找出了物理实验结果Kiya et al (1980)与近年来数值模拟结果Kang (2007)存在较大差异的主要原因。指出无滑移边界条件将能有效抑制尾涡的脱落,而自由表面边界条件则对尾涡的脱落没有多大的影响。因而,物理实验的实施必须考虑到无滑移边界条件的影响,而数据的采集应尽量避免处在无滑移边界条件影响的区域之内。 3数值计算得出了两向自由度弹性支撑圆柱体产生较大横向振幅的临界质量比* m = 3.5。当圆柱体的质量比高于临界质量比时,限制流向运动对横向振幅幅值的影响不大;当质量比降至临界质量比时,流向和横向的振动将激发出更大的横向振幅幅值,流向运动对横向振幅的影响在此凸现。数值计算还重现了实验中不同横向振幅响应分支对应的尾涡模式,特别是与超上端分支对应的2T尾涡模式。在均匀来流条件下,数值模拟还成功得出了圆柱体“8”型的运动轨迹。 4在横向剪切来流条件下,两向自由度弹性支撑圆柱体的涡激振动存在一个剪切率的锁定区间。剪切率锁定区间与雷诺数的大小成正比。在剪切率的锁定区间内,横向振幅增大;脱离剪切率锁定区间后,横向振幅迅速下降。在横向来流剪切率较小时,圆柱体的运动轨迹与均匀来流条件下的“8”型的运动轨迹一致。随着剪切率的增加,圆柱体“8”型运动轨迹上端消失而表现为“小雨点”的外形。 5在阶梯状来流条件下,立管的响应在时域内并不表现为单一模态,而是在各个模态之间进行跳转。当某一激发模态的权重占绝对主导地位时,立管的变形与该激发模态下的形状一致;而当立管激发各模态的权重并不由唯一的激发模态所主导时,立管的响应形式将在占主导地位的几种激发模态之间跳转。通过记录时域内沿立管轴向分布的流体力,发现流体力沿立管轴向具有分段相关性,这为本文基于切片法所建立的时域预报模型的合理性提供了依据。沿立管轴向分布的附加质量系数存在波动区间,这与立管呈现的多模态响应相关联。当立管的涡激振动表现为多模态响应时,本文提出了附加质量系数的计算应当根据激发模态的权重来考虑各个模态所对应附加质量系数的贡献。数值计算沿立管轴向同时捕捉到了2S和2P尾涡模式,重现了这一经典的物理现象;在不同的立管截面之间,还观测到尾涡泻涡时间的切换现象。 6在轴向剪切来流条件下,本文的时域预报模型所得到的低阶模态响应与实验结果吻合较好。本文采用的时域模型采用了较多的切片数量来考虑流体力的轴向相关性,利用动网格和变形网格来模拟流固耦合运动对流体流动的影响,还同步考虑到立管的变形引起的总体刚度矩阵的非线性性。因而,本文所预报的结果较同类切片模型和频率预报模型的结果更为精确。 7通过沟通受迫振动和自激振动内在联系,建立理论模型成功实现了由受迫振动实验数据来预报自激振动圆柱体的横向振幅幅值。该理论模型可同时考虑到雷诺数和质量阻尼比对横向振幅幅值的影响,同时也再次证明雷诺数是影响涡激振动横向振幅的一个重要参数。本文的理论预报模型可为深海立管频域预报模型的发展和完善提供理论依据。 通过对本文主要内容和结论的阐述,可以反映出该论文的创新点表现在以下几个方面: 1利用直接数值模拟方法首次对横向剪切来流条件下的圆柱绕流进行了计算,对尾涡模式和流体流动的物理特性进行了详细探讨。 2对两向自由度弹性支撑圆柱体在均匀来流和横向剪切流条件下进行了数值模拟,得出了均匀来流条件下产生较大横向振幅所对应的临界质量比为3.5;得出横向剪切流条件下存在剪切流的锁定区间,并发现了一种全新的“小雨点”运动轨迹。 3建立了柔性立管涡激振动时域响应预报模型,在时域内成功预报了立管的模态响应。 4发展和完善了频域预报模型。在受迫振动实验数据基础上,考虑了雷诺数和质量阻尼比得出了自激振动横向振幅峰值的理论预报公式。 总而言之,本文从简单的圆柱绕流开始,逐步考虑到弹性支撑圆柱体和柔性立管的涡激振动响应,并成功建立了柔性立管涡激振动的时域预报模型。本文时域模型的预报结果能准确地反映出物理实验中的诸多现象和结论,同时,从本文时域模型得到的结论中可为完善和发展柔性立管频率预报模型提供思路和参考。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 潘志远;崔维成;刘应中;;阶梯状来流中立管的涡激振动响应预报[J];上海交通大学学报;2006年06期
2 王志国,陈启富,李维扬;圆柱的涡激振动估算[J];哈尔滨工程大学学报;1987年01期
3 徐枫;欧进萍;;低雷诺数下弹性圆柱体涡激振动及影响参数分析[J];计算力学学报;2009年05期
4 黄智勇;崔维成;黄小平;;均匀平面剪切流作用下圆柱体的涡激振动[J];哈尔滨工程大学学报;2007年12期
5 刘兆存,金忠青;二维平行壁面剪切流失稳内部流动结构初探[J];水利水电科技进展;1995年03期
6 王一飞;潘志远;黄小平;崔维成;;深海立管涡激振动疲劳损伤影响因素分析[J];船舶力学;2006年05期
7 王一飞;吴晓源;潘志远;黄小平;崔维成;;弹性模量对深海立管涡激振动疲劳损伤的影响(英文)[J];船舶力学;2007年06期
8 陈伟;宗智;;立管涡激振动的离散涡数值模拟[J];中国海洋平台;2010年01期
9 玄金;黄维平;曹静;;深水顶张力立管涡激振动疲劳损伤分析[J];船海工程;2009年04期
10 陈起富;姚熊亮;李维扬;;并列双柱圆柱体在均匀流场中的涡激振动[J];哈尔滨工程大学学报;1991年01期
11 任铁;付世晓;李润培;杨建民;;基于模型试验的实尺度立管涡激振动响应预报方法研究[J];船舶力学;2011年04期
12 赵秋亮;;海底电力电缆涡激响应的数值求解方法[J];广西质量监督导报;2008年08期
13 潘志远;崔维成;张效慈;;细长海洋结构物涡激振动研究综述(英文)[J];船舶力学;2005年06期
14 黄维平;曹静;张恩勇;唐世振;;大柔性圆柱体两自由度涡激振动试验研究[J];力学学报;2011年02期
15 荣深涛;杨健;徐岩;;科特(Couette)数Pr=1的环隙科特紊流[J];北京交通大学学报;1991年01期
16 雷松;张文首;岳前进;赵岩;;FDPSO张力系统的非线性动力分析[J];中国造船;2010年03期
17 饶志标;杨建民;付世晓;李润培;;剪切流下钢悬链线立管涡激振动响应研究[J];振动与冲击;2010年10期
18 上官丽红;;顶张力对深水刚性立管涡激振动及疲劳损伤的影响[J];船海工程;2011年02期
19 常向东;朱晨辉;;河流穿越管道冲管段涡激振动研究[J];管道技术与设备;2008年02期
20 唐世振;黄维平;刘建军;;深水立管两向自由度涡激振动的数值分析[J];振动与冲击;2010年09期
中国重要会议论文全文数据库 前10条
1 黄维平;孙铭远;;圆柱管束涡激振动实验分析研究[A];第十五届中国海洋(岸)工程学术讨论会论文集(上)[C];2011年
2 林黎明;宋芳;曾晓辉;;应用新的尾流振子模型研究大长径比下立管的涡激振动[A];北京力学会第十六届学术年会论文集[C];2010年
3 王懿;段梦兰;曹静;;深水钢质悬链线立管涡激振动疲劳寿命分析[A];第十五届中国海洋(岸)工程学术讨论会论文集(上)[C];2011年
4 周阳;;螺旋侧板涡激振动控制实验研究[A];第十五届中国海洋(岸)工程学术讨论会论文集(上)[C];2011年
5 葛斐;王雷;陆维;洪友士;;剪切流中大长细比圆柱体的二维涡激振动[A];第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会论文集[C];2009年
6 汪源源;王威琪;;利用超声多普勒技术检测血流的湍流和涡流[A];首届全国医疗器械学术与产业论坛论文集[C];2002年
7 李敢;黄小平;崔维成;;钢悬链线立管涡激疲劳损伤敏感参数研究[A];纪念徐秉汉院士船舶与海洋结构力学学术会议论文集[C];2011年
8 阳祥;李增耀;陶文铨;;离心浮升力对旋转方通道内湍流混合对流影响的数值模拟[A];中国力学学会学术大会'2009论文摘要集[C];2009年
9 魏东泽;黄维平;;深水刚性立管与弹性立管涡激振动分析[A];第十五届中国海洋(岸)工程学术讨论会论文集(上)[C];2011年
10 郑文青;杨和振;;钢悬链立管顶端柔性接头的疲劳分析研究[A];中国钢结构协会海洋钢结构分会2010年学术会议暨第六届理事会第三次会议论文集[C];2010年
中国博士学位论文全文数据库 前10条
1 黄智勇;柔性立管涡激振动时域响应分析[D];上海交通大学;2008年
2 张永波;深海输液立管涡激振动预报及抑振技术研究[D];中国海洋大学;2011年
3 唐国强;立管涡激振动数值模拟方法及物理模型实验[D];大连理工大学;2011年
4 宋吉宁;立管涡激振动的实验研究与离散涡法数值模拟[D];大连理工大学;2012年
5 高云;钢悬链式立管疲劳损伤分析[D];大连理工大学;2011年
6 李朋;深水海洋输流立管涡激振动及干涉试验研究[D];中国海洋大学;2012年
7 孟丹;钢悬链线输流立管非线性振动数值模拟[D];中国海洋大学;2010年
8 戴伟;基于响应面方法的立管结构可靠性研究[D];哈尔滨工程大学;2012年
9 李效民;顶张力立管动力响应数值模拟及其疲劳寿命预测[D];中国海洋大学;2010年
10 王一飞;深海立管涡激振动疲劳损伤预报方法研究[D];上海交通大学;2008年
中国硕士学位论文全文数据库 前10条
1 张友林;海洋立管涡激振动抑制方法研究[D];江苏科技大学;2011年
2 董婧;挠性立管涡激振动的离散涡方法研究[D];大连理工大学;2010年
3 刘庆海;钢悬链线立管涡激振动试验研究[D];中国海洋大学;2012年
4 任铁;实尺度雷诺数复杂来流条件下的立管涡激振动试验研究[D];上海交通大学;2011年
5 李相环;海洋立管涡激振动抑振装置的试验研究[D];中国海洋大学;2012年
6 顾恩凯;考虑流固耦合的大位移柔性立管涡激振动研究[D];中国海洋大学;2012年
7 陈伟;涡激振动的离散涡数值模拟[D];大连理工大学;2009年
8 白长旭;波流耦合作用下的立管涡激振动分析[D];大连理工大学;2008年
9 吴刚;海洋工程结构振动与声学特性计算[D];上海交通大学;2008年
10 苑健康;深水顶张式立管螺旋侧板抑制VIV机理研究[D];中国海洋大学;2010年
中国重要报纸全文数据库 前10条
1 郭志平;研究湍流的后起之秀[N];浙江日报;2000年
2 宋方;全时域升空 全天候作战[N];解放军报;2005年
3 本报记者 连锦添;遭遇湍流[N];人民日报海外版;2002年
4 记者王春;国际湍流论坛再掀漩流[N];科技日报;2002年
5 任一言;佘振苏:一位长江学者的承诺[N];光明日报;2007年
6 上海社会科学院研究员 花建;在文化投资的湍流中[N];中国文化报;2002年
7 石岩;学之大者 融通中西 为国为民 不懈追求[N];人民政协报;2009年
8 鲍勇剑;管理与力学[N];21世纪经济报道;2007年
9 本报记者 宋茜 郑颖 张武明;越过湍流是坦途[N];江西日报;2009年
10 ;溶解氧(DO)[N];中国水利报;2004年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978