Euler-Bernoulli海洋立管涡致强迫振动响应研究.pdf

1958西南石油大学学报（自然科学版）2023 年 8 月 第 45 卷 第 4 期Journal of Southwest Petroleum University（Science&Technology Edition）Vol.45 No.4 Aug.2023DOI：10.11885/j.issn.1674 5086.2021.02.23.01文章编号：1674 5086（2023）04 0133 10中图分类号：TE973；O327文献标志码：AEuler-Bernoulli 海洋立管涡致强迫振动响应研究赵 翔1*，谭 明1，李映辉2，邵永波11.西南石油大学土木工程与测绘学院，四川 成都 6105002.西南交通大学力学与工程学院，四川 成都 610031摘要：针对海洋立管（Pipe-in-pipe，PIP）系统在海水作用下发生的振动问题，开展了对 PIP 系统在涡致强迫振动下的动力学响应研究，分析了在涡致强迫振动下海洋立管外管直径、轴向拉力、外激力频率对海洋立管位移响应的影响规律。基于 Euler-Bernoulli 双梁模型，采用 Lamb-Oseen 涡模型，建立了动力学模型，利用格林函数法求得该强迫振动的稳态响应。结果表明，随着管道直径增加，外激力增加，产生最大力幅值的位置离管道越远；轴向拉力对外部管道的影响较大，对内部管道的影响较小；无因次频率取 0.4 时，外部管道位移超出允许变形极限，内外管壁发生周期碰撞，易对海洋立管造成损伤。关键词：海洋立管；涡致强迫振动；稳态响应；格林函数法；Euler-Bernoulli 双梁Vortex Induced Forced Vibration of Euler-BernoulliPipe-in-pipe SystemZHAO Xiang1*,TAN Ming1,LI Yinghui2,SHAO Yongbo11.School of Civil Engineering and Geomatics,Southwest Petroleum University,Chengdu,Sichuan 610500,China2.School of Applied Mechanics and Engineering,Southwest Jiaotong University,Chengdu,Sichuan 610031,ChinaAbstract:In order to solve the vibration problem of a Pipe-in-pipe(PIP)system under the action of seawater,the dynamicresponse of the vortex induced forced vibration of the PIP system is studied.Then the numerical analysis is carried out to studyinfluences of the outer pipe diameter,axial load,and external excitation frequency on displacement responses of the PIP systemunder the action of vortex.Based on Euler-Bernoulli beams theory,Lamb-Oseen vortex model is used to establish the dynamicmodel the PIP system,and Green function method is used to obtain the steady-state response of the forced vibration of the PIPsystem.Results show that as the diameter of the outer pipe increases,the vortex excitation force increases,and the position ofthe maximum amplitude becomes farther from the outer pipe;the axial tension has a greater impact on the outer pipe,but has asmaller effect on the inner pipe;when the dimensionless external excitation frequency sets to 0.4,the displacement of the outerpipe exceeds the limit of allowable deformation,and periodic collisions between the inner and outer pipe can easily happen,which cause damage to the PIP system.Keywords:pipe-in-pipe;vortex induced forced vibration;steady-state response;Greens function method;Euler-Bernoullidouble-beams网络出版地址：http：/ 翔，谭 明，李映辉，等.Euler-Bernoulli 海洋立管涡致强迫振动响应研究J.西南石油大学学报（自然科学版），2023，45（4）：133 142.ZHAOXiang,TAN Ming,LI Yinghui,et al.Vortex Induced Forced Vibration of Euler-BernoulliPipe-in-pipe SystemJ.Journal of Southwest PetroleumUniversity(Science&Technology Edition),2023,45(4):133142.*收稿日期：2021 02 23网络出版时间：2023 07 05通信作者：赵 翔，E-mail：zhaoxiang_基金项目：国家自然科学基金（12072301，11872319）；四川省科技厅国际合作项目（2020YFH0127）；工程结构安全评估与防灾技术四川省青年科技创新团队（2019JDTD0017）134西南石油大学学报（自然科学版）2023 年引言随着陆地油气资源的不断开采，其储备量逐渐枯竭，包括中国在内的许多国家都将目标转向了海洋。在海洋油气开采的过程中，海洋立管是极其关键的设备之一，主要连接海上悬浮平台与井口，为采集油气提供输送管道。海洋立管在工作时受复杂的工作环境以及洋流、波流等各种动态荷载的影响，会引起立管的振动1，因此，研究海洋立管的振动特性，确保海洋立管的安全性以及稳定性尤为重要。PIP 系统作为一种海洋管道结构，具有较好的抗振效果，近年来得到了广泛应用。PIP 系统通常由外部管道、内部管道和填充在两个管道之间的绝缘层 3 部分构成。通过绝缘层连接外部和内部管道的系统以及由大量舱壁连接的系统是常见的 PIP 系统类型2。内部管道用于输送石油等流体，外部管道涂有防腐材料，并保护内部管道免受冲击。已经有许多研究人员通过分析或者实验的方式研究了PIP 系统的静态特性和动力学问题3 4。PIP 系统在流体作用下发生涡致强迫振动是非常普遍的现象。当涡致强迫振动的频率接近 PIP 系统自身的频率时，会产生大幅度的振动，从而形成较大威胁5。现有研究海洋立管受涡致强迫振动的响应方法有数值方法6 9和实验方法10 12。此前的研究工作中有关轴向拉力和阻尼对振动特性影响的研究较少13。实验方法获得的数据可靠但研究成本较高，而数值方法常用的是经验模型法中的尾流振子模型法，这种模型受到广泛认可。关于抑制立管振动的研究比较多，例如，增加不同截面形状的肋条，来研究其对振动响应的影响等14。尾流振子模型属于非线性振动模型，研究发现基于加速度模型来描述立管尾部流场的尾流振子模型更加适用15。本文采用的 Lamb-Oseen 涡模型属于线性振动模型，这种线性模型在许多压电能量采集器的设计和实验中应用较多，它主要是吸收涡脱落的振动能量，并且在实验中有较好的效果。海洋立管是一种细长结构体，其跨径比较大，并且管道的弯曲变形与剪切变形可以忽略不计，因此，本文采用 Euler-Bernoulli 双梁模型对管道进行模拟16 17。基于 Euler-Bernoulli 双梁模型，Zhang等18研究了在轴向压缩荷载作用下双梁系统自由振动问题。Mao 等19 20采用 Adomian 分解方法研究了 Winkler 层连接的双梁系统的自由振动问题。Abu-hilal 等21通过拉普拉斯变换得到了一般边界条件下 Euler-Bernoulli 双梁模型稳态响应的格林函数。Li 等22建立了具有转动阻尼和平移阻尼的Timoshenko 梁模型，并得到了 Timoshenko 梁稳态响应的格林函数。Zhao 等23研究了多裂纹 Euler-Bernoulli 梁的受迫振动，得到了多裂纹梁的格林函数。他们所采用的梁模型以及格林函数法为本文的研究奠定了理论基础。目前，对 PIP 系统响应的研究主要是计算流体力学分析（Computational Fluid Dynamics，CFD）及实验方法，并未给出确定的解析解，对于 PIP 系统响应的快速预测具有一定的阻碍。基于此，采用 Lamb-Oseen 涡模型对 PIP 系统的外激力建模来研究其涡致强迫振动问题，求解其涡致强迫振动问题的闭形式解析解，并对 PIP 系统的响应进行快速预测，这对研究 PIP 系统涡致强迫振动问题是十分有意义的。本文对海洋立管的研究以横向振动为主，分析了 Euler-Bernoulli 双梁模型在涡致强迫振动作用下的动力学响应。通过建模和拉普拉斯变换推导动力学问题的格林函数，得到了海洋立管的响应。通过数值分析，验证了双梁模型的正确性，并分别对不同管道直径、轴向拉力及外激力频率分析了系统响应的规律。1 PIP 系统的模型和动力学响应图 1 给出了 PIP 的模型示意图。管道系统的长度为 L。外部管道的外径和内径分别为 D1和 D2，内部管道的外径和内径分别为 d1和 d2，U0表示流体的平均速度。PIP 系统通常具有很高的细长比24 26，因此，采用 Euler-Bernoulli 双梁模型，并且将中间绝缘层视为线性弹簧和阻尼器。!#$%&(#$xyzD1D2d1d2o图 1PIP 系统示意图Fig.1Schematic diagram of PIP system第 4 期赵 翔，等：Euler-Bernoulli 海洋立管涡致强迫振动响应研究1351.1 格林函数的解本文考虑了在涡致强迫振动影响下海洋立管的动力学响应，如图 2 所示，建立系统的振动方程为16 17E1I1w1+m34(w1+w2)+m1 w1+k(w1 w2)+c(w1 w2)+F1w1=p(x,t)E2I2w2+m34(w1+w2)+m2 w2 k(w1 w2)c(w1 w2)+F2w2=0(1)式中：E1，E2外部、内部管道的弹性模量，Pa；I1，I2外部、内部管道的惯性矩，m4；w1，w2w1、w2对 x 的四阶导数；m1，m2，m3外部、内部管道以及绝缘层单位长度的质量，kg；w1，w2w1、w2对时间的二阶导数，m/s2；k绝缘层弹性模量，Pa；w1，w2外部、内部管道的横向位移，m，与管道位置 x10和时间 t 相关；c绝缘层单位长度的阻尼，Ns；w1，w2w1、w2对时间的一阶导数，m/s；F1作用在外部管道的轴向拉力，N；F2作用在内部管道的轴向拉力，N；w1，w2w1、w2对 x 的二阶导数；p(x,t)流体产生作用在外部管道上的外激力，N；t时间，s。F1F2F1F2w1()x t，!#$%#$w2()x t，xp()x t，ckL图 2具有黏弹性层的双梁系统Fig.2Double-beam system with a viscoelastic layer假 设 外 激 力 p(x,t)=P(x)eit（式 中：P(x)外激力幅值，N；x管道的长度方向；i虚数单位；漩涡脱落角频率，rad/s，且=2f，f漩涡脱 落 频 率，Hz），相应 的 挠 度 也 可 以表 示为wn(x,t)=Wn(x)eit（式中：Wn管道的横向位移，m，与管道位置 x10相关）。为了简化系统方程，将上述假定代入式（1）中，消去时间变量。因此，可以得到关于 Wn（x）的微分方程E1I1W1m324(W1+W2)m12W1+k(W1 W2)+ic(W1 W2)+F1W1=P(x)E2I2W2m324(W1+W2)m22W2 k(W1 W2)ic(W1 W2)+F2W2=0(2)式中：W1，W2W1、W2对 x 的四阶导数，m3；W1，W2W1、W2对 x 的二阶导数，m1。简化成下列形式W1+a1W1+a2W1+a3W2=CP(x)W2+b1W2+b2W2+b3W1=0(3)其中：a1=F1E1I1，a2=m342 m12+k+icE1I1，a3=m342 k icE1I1，C=1E1I1，b1=F2E2I2，b2=m342 m22+k+icE2I2，b3=m342 k icE2I2格林函数是已经被证实的获得线性系统动态响应的有效方法16。双梁系统稳态强迫振动的格林函数是在单位谐波载荷作用下系统的动态响应。式（3）的格林函数解为W1+a1W1+a2W1+a3W2=C(x x0)W2+b1W2+b2W2+b3W1=0(4)式中：()狄拉克函数。将式（4）中的 x 进行 Laplace 变换和逆变换计算，可以推导出格林函数G1(x;x0)=11(x)W1(0)+21(x)W1(0)+31(x)W1(0)+41(x)W1(0)+51(x)W2(0)+61(x)W2(0)+71(x)W2(0)+81(x)W2(0)+H(x x0)91(x x0)(5)136西南石油大学学报（自然科学版）2023 年G2(x;x0)=12(x)W1(0)+22(x)W1(0)+32(x)W1(0)+42(x)W1(0)+52(x)W2(0)+62(x)W2(0)+72(x)W2(0)+82(x)W2(0)+H(x x0)92(x x0)(6)式中：G1(x;x0)，G2(x;x0)外部管道和内部管道的Green 函数；H(x x0)单位阶跃函数。其中，系数 ij（i=19，j=12）和常数 W1(0)、W2(0)、W1(0)以及 W2(0)均可通过文献16的计算方法，结合边界条件求解得到，从而得到 G1(x;x0)和G2(x;x0)G1(x;x0)=21(x)W1(0)+41(x)W1(0)+61(x)W2(0)+81(x)W2(0)+H(x x0)91(x x0)(7)G2(x;x0)=22(x)W1(0)+42(x)W1(0)+62(x)W2(0)+82(x)W2(0)+H(x x0)92(x x0)(8)当双梁体系受荷载时，控制方程的解为W1(x)W2(x)=L0G1(x;x0)G2(x;x0)P(x0)dx0(9)2 外激力计算图 3 为流体对海洋立管作用的示意图。本文采用的是 Lamb-Oseen 涡流模型。假设涡致强迫振动产生的外激力为简谐力：p(x,t)=P(x)eit，假设振动是周期性的，且管道上的力分布均匀，则 P(x,t)只与时间相关，不随 x 的变化而变化。其中，f 的表达式为27f=St(U0/D)(10)式中：StStrouhal 数，St=0.21；U0平均流体速度，m/s；D钝体直径，m。U0DrdrA y z(,)C(,)y zccyzo图 3涡致强迫振动对管道的作用示意图Fig.3The action diagram of vortex-induced vibration on pipeline注：r漩 涡 中 心 一 点 到 海 洋 立 管 上 一 点 的 距 离，m，r=(y yc)2+(z zc)2；dr漩涡中心与钝体中心水平方向上的相对距离，m。在理论分析中，为了考虑涡流对立管的影响，漩涡产生的距立管中心最小和最大距离分别为0.5D 和 15.0D，超出此范围后漩涡不再对立管产生影响。此外，在理论分析中，外激力频率与涡流脱落频率相等，即可导出外激力幅值 P(x)的表达式。聚焦在一个稳定的涡面上，假设漩涡中心位置C(yc,zc)，A(y,z)为外部管道上一点，速度为 va。则A 点和 C 点之间的压差为27p=12(v2c v2a)(11)式中：流体的密度，kg/m3；vc漩涡中心的平均水平速度，m/s；va外部管道 A 点的速度，m/s。根据 Lamb-Oseen 涡流模型假设和特性，单个涡流的运动可以分解为平均水平速度 vc和旋转速度 c，点 A 相对于涡流中心 C 的切向速度可以建模为vc va=cr(12)式中：c涡核的角速度，rad/s。根据斯托克斯定理，平均转速也可以表示为 c=/2；涡流强度，rad/s。这种压差产生了主要作用在管道上的外激力，表达式如下p=2(vcr 2r2/4)(13)P(x)=lpds(14)式中：l海洋管道截面曲线。将压差 p 代入式（14），得到P(x)=2lvc(y yc)2+(z zc)2 2(y yc)2+(z zc)2/4ds(15)式（15）表明施加在管道上的外激力主要受涡流强度、位置和水平速度的影响。第 4 期赵 翔，等：Euler-Bernoulli 海洋立管涡致强迫振动响应研究137其中=(0/4vet)er20/4vet0(16)式中：0漩涡从圆柱上脱落的初始环流，无因次，0=(U0D)/2St；r0刚性涡核的半径，m；假设 r0=zc；ve漩涡的等效耗散因子，ve=0.002；t0初始漩涡产生时间，s，t0=dr/U0。此外，平均流体速度 U0可以近似为漩涡的平均水平速度 vc。文 献27表 明，当 圆 柱 上 的 漩 涡 稳 定 时，zc=1.3D/2。从而得到外激力关于 yc的关系，最后得到外激力达到最大时的位置 yc。单位长度外激力幅值 P(x)可表示为P(x)=2lvc(y yc)2+(z 1.3D2)2 0.252(y yc)2+(z 1.3D2)2ds(17)再将 P(x)代入式（9），即可得到 PIP 系统的响应。3 数值分析与讨论为了方便计算，管道长度取 10 m，外部管道的外径 D1=0.26 m，内径 D2=0.22 m；内部管道的外径 d1=0.20 m，内径 d2=0.16 m；管道密度取7 850 kg/m3，弹性模量 2.01011Pa，海水密度为1 025 kg/m3。为了方便讨论，引入无因次量X=xLgi(X;X0)=Gi(x;x0)Lk0=kEF=FiPEWi(X)=Wi(X)L=0式中：X无因次坐标；gi(X;X0)无因次挠度；k0无因次弹性模量；F无因次轴向拉力；Wi(X)无因次位移；无因次频率；PEEuler-Bernoulli 双梁的临界屈曲荷载，N，PE=EI2/L2；0简支 Euler-Bernoulli 双梁的一阶固有频率，Hz，0=2(EI/A)0.5/L2。3.1 解的有效性3.1.1 格林函数解的有效性实际上，当绝缘层的弹性模量以及阻尼为零时，忽略绝缘层质量，上述解就是单梁的格林函数。即 g1(X;X0)是 k=0 情况下单梁的格林函数，并且在这种特定情况下，显然可知 g2(X;X0)=0。图 4显示了无因次挠度 gi(X;0.5)和无因次坐标 X的函数。由图4可以看出，双梁系统的无因次挠度g1(X;X0)与单梁的无因次挠度 g(X;X0)完全吻合，双梁系统的无因次挠度 g2(X;X0)为零，验证了本文方法的正确性。00.21070.41070.61070.81071.010700.20.40.60.81.0X!?22#$%&$%&1.2107g X 0.5i();图 4解的有效性验证Fig.4Validation of the solution3.1.2 响应解的有效性通过前面的计算，得到了 PIP 系统涡致强迫振动位移响应，并且借助现有的研究结果与本文的结果进行对比验证，其基本参数如表 1 所示。由于文献28中的立管为单管，所以在这里需要进行一些调整，令绝缘层的弹性模量以及阻尼为零，忽略绝缘层质量，从而退化成单立管模型，以方便数据的分析以及对比。文献28中的数据是由柳军采用 CFD 方法得到的海洋立管位移包络图。由于 Lamb-Oseen 涡模型是将流体对 PIP 系统的作用简化成均匀分布在管道上的简谐力，并且忽略了管道振动对流体的作用，因此，计算结果与文献数据存在一定的差异（图 5）。通过对比，本文的结果和文献结果在同138西南石油大学学报（自然科学版）2023 年一外部流体速度下的各点横向最大位移和振动形态基本一致，可证明采用 Lamb-Oseen 涡模型进行模拟预测是可行的。表 1基本参数Tab.1Foundamental parameters立管长度/m立管材料密度/（kgm3）外径/m流体密度/（kgm3）立管壁厚/m1207 8500.261 0250.019 2流体速度/（ms1）顶张力/kNStrouhal 数弹性模量/GPa0.542000.2200020406080100120-0.8-0.6-0.4-0.200.20.40.60.8&/mW XD11()/!#$%?28图 5PIP 系统位移响应验证Fig.5Displacement response verification of PIP system3.2 外部管道外径的影响如图 6 所示，表示外部管道外径对外激力的影响，且 F1=F2=0，k=5 1010GPa，c=0.059 Ns，=0.2，忽略绝缘层质量。外激力计算采用Lamb-Oseen 涡模型，从文献27中可知，外激力恒大于零。从图 6 中可以看出：当外部管道 D1=0.26 m 时，在 dr=2.0 m 时达到幅值；当外部管道 D1=0.28 m时，在 dr=2.6 m 时达到幅值；当外部管道 D1=0.30 m时，在 dr=3.2 m 时达到幅值；当外部管道 D1=0.32 m时，在 dr=3.8 m 时达到幅值；当外部管道 D1=0.34 m时，在 dr=4.6 m 时达到幅值。由此得出结论：外部管道外径越大，流体产生的外激力幅值越大，并且产生幅值的位置离管道越远。00123451.01.52.02.5P X()/kNdr/mD1=0.26D1=0.28D1=0.30D1=0.32D1=0.340.5图 6外部管道外径对外激力的影响Fig.6The influence of pipeline diameter on vortex induced force图 7 为外部管道外径对系统的影响。从图 7可以看出，由于流速较小只能激发管道的一阶模态，即在 X=0.5 处位移达到最大；在 不变的情况下，随着管道直径的增加，外部管道的响应也增加，内部管道的响应减小，内外管道的位移响应存在较大差距，所以管道直径不宜取太大，从而减小振动响应。00.11030.21030.31030.41030.51030.610300.20.40.60.81.0X0.7103WX1()00.20.40.60.81.0X00.11060.21060.31060.41060.5106WX2()D1=0.26D1=0.28D1=0.30D1=0.32D1=0.34D1=0.26D1=0.28D1=0.30D1=0.32D1=0.34a!#$b%#$图 7外部管道外径对系统的影响Fig.7Effect of pipe diameter on system3.3 无因次轴向拉力的影响图 8 为不同无因次轴向拉力和管道模型位移响应的关系图，且 F1=F2，c=0.059 Ns，=0.2，k=51010GPa。从图8可以看出，当 F在0.20.8变化时，外部管道位移逐渐增大，并且在 0.60.8 变化时位移变化幅度较大；内部管道受无因次轴向拉力的影响较小，位移响应变化较小，当 F=0.4 时，内部管道位移发生较大变化。第 4 期赵 翔，等：Euler-Bernoulli 海洋立管涡致强迫振动响应研究13900.51031.01031.51032.01032.5103F=0.2F=0.4F=0.6F=0.8-2.5105-2.0105-1.5105-1.0105-0.5105000.20.40.60.81.0X00.20.40.60.81.0X0.5105WX2()3.0103WX1()F=0.2F=0.4F=0.6F=0.8a!#$b%#$图 8不同无因次轴向拉力对应的管道位移响应图Fig.8The response diagram of pipe displacement corresponding to different tensile forces3.4 无因次频率的影响流体速度发生改变时，作用在海洋立管上的外激力及频率发生改变，这会对管道模型的位移产生影响。从前面计算可知，流速和无因次频率成正相关。如图 9 所示，绘制了无因次频率和外激力的关系图，在本算例中，k=5 1010GPa，c=0.059 Ns。D1=0.26 m，L=3.9 m。001234dr/m12345P X()/kN=0.2=0.4=0.6=0.8=.1?0=.21图 9无因次频率与外激力的关系图Fig.9Diagram of dimensionless extraneous excitation frequency andvortex induced force从图9可以看出，当=0.2时，dr=2.0m时 P(X)达到幅值；当=0.4 时，dr=3.2 m 时 P(X)达到幅值；当 0.6时，dr=3.9m时 P(X)达到幅值；当=0.6时，P(X)达到最大幅值。由此得出结论：随着无因次频率的增大，外激力幅值先增大后减小，并且产生幅值的位置离管道越来越远。图 10 为无因次频率与外部管道位移的关系图，k=5 1010GPa，c=0.059 Ns，F=0.2。从图 10可以看出，随无因次频率的增大，外部管道位移先增大后减小，在=0.8 时发生反向。当 达到 0.4时，外部管道位移超出允许变形范围极限，允许变化范围极限是中间绝缘层的厚度 0.02 m，此时，可判断内外管壁发生周期碰撞，易对海洋立管造成损伤。020406080100120140WX1()0.40.60.81.0Xa=.2.00?800.2 0=.2 0=.4 0=.6 0=.80-500-400-300-200-100100WX1()00.20.40.60.81.0Xb=.2.1 18 1=.2 1=.4 1=.6 1=.8图 10无因次频率与外部管道位移响应图Fig.10Diagram of dimensionless frequency and outerpipeline response如图 11 所示，绘制了无因次频率与内部管道位移的关系图，绝缘层弹性模量为 5 1010GPa，阻尼为 0.059 Ns，无因次轴向拉力 0.2。从图 11 可以看出，随无因次频率的增大，内部管道位移先增大后减小，在=0.6 时发生反向。相比于外部管道，内部管道的位移较小。随着无因次频率的继续增大，便会激发系统高阶模态。从图 12 中可以看出，随着无因次频率的继续增加，外部管道和内部管道的位移越来越小。140西南石油大学学报（自然科学版）2023 年001510251000.20.40.60.81.0X00.20.40.60.81.0X-6000-4000-20002000WX2()3510WX2()0=.2 0=.4 0=.6 0=.8 1=.2 1=.4 1=.6 1=.8a=.2.00?8b=.2.1 18图 11无因次频率与内部管道响应的关系图Fig.11Diagram of dimensionless frequency and inner pipeline responsea!#$%&()*+,b-#$%&()*+,-32010-22010-120100120102201000.20.40.60.81.0X32010WX2()-252910-202910-152910-102910-5291005291000.20.40.60.81.0X102910WX1()图 12无因次较大频率与管道系统响应的关系图Fig.12Diagram of dimensionless higher frequency and piping system response4 结论1）在其他条件不变的情况下，管道直径增加，外激力增大，外激力产生最大幅值的位置远离管道，外部管道的振动响应增大，内部管道的振动响应减小。2）轴向拉力的增加会增大外部管道的位移响应，影响较大，对内部管道的影响较小。3）随着无因次频率的增大，外部管道和内部管道的位移都先增大后减小。无因次频率取 0.4 时，外部管道位移超出允许变形极限，此时，可判断内外管壁发生周期碰撞，易对海洋立管造成损伤。参考文献1姜峰，郑运虎，梁瑞，等.海洋立管湿模态振动分析J.西南石油大学学报（自然科学版），2015，37（5）：159 166.doi：10.11885/j.issn.1674-5086.2013.11.11.03JIANG Feng,ZHENG Yunhu,LIANG Rui,et al.An ana-lysis of the wet modal vibration of marine riserJ.Jour-nal of Southwest Petroleum University(Science&Tech-nology Edition),2015,37(5):159166.doi:10.11885/j.-issn.1674-5086.2013.11.11.032KARAMPOUR H,ALRSAI M,ALBERMANI F,etal.Propagation buckling in subsea pipe-in-pipe sys-temsJ.Journal of Engineering Mechanics,2017,143(9):04017113.doi:10.1061/(ASCE)EM.1943-7889.00013373LI M,ZHAO X,LI X,et al.Stability analysis of oil-conveying pipes on two-parameter foundations with ge-neralized boundary condition by means of Greens func-tionsJ.Engineering Structures,2018,173:300312.doi:10.1016/j.engstruct.2018.-07.0014ALRSAI M,KARAMPOUR H,ALBERMANI F.On col-lapse of the inner pipe of a pipe-in-pipe system underexternal pressureJ.Engineering Structures,2018,172:614628.doi:10.1016/j.engstruct.2018.06.0575GUO Yang,ZHU Bo,ZHAO Xiang,et al.Dynamic cha-racteristics and stability of pipe-in-pipe system conveyingtwo-phaseflowinthermalenvironmentJ.AppliedOceanResearch,2020,103:102333.doi:10.1016/j.apor.2020.-1023336GAO Yun,FU Shixiao,WANG Jungao,et al.Experimen-tal study of the effects of surface roughness on the vortex-inducedvibrationresponseofaflexiblecylinderJ.Ocean第 4 期赵 翔，等：Euler-Bernoulli 海洋立管涡致强迫振动响应研究141Engineering,2015,103:4054.doi:10.1016/j.oceaneng.-2015.04.0527宋磊建，付世晓，于大鹏，等.剪切流下发生涡激振动的柔性立管阻力特性研究J.力学学报，2016，48（2）：300 306.doi：10.6052/0459-1879-15-309SONG Leijian,FU Shixiao,YU Dapeng,et al.Inves-tigation of drag forces for flexible risersundergoing vor-tex-induced vibration in sheared flowJ.Chinese Journalof Theoretical and Applied Mechanics,2016,48(2):300306.doi:10.6052/0459-1879-15-3098徐万海，马烨璇，罗浩，等.柔性圆柱涡激振动流体力系数识别及其特性J.力学学报，2017，49（4）：818 827.doi：10.6052/0459-1879-16-263XU Wanhai,MA Yexuan,LUO Hao,et al.Identifica-tion and characteristics of hydrodynamic coefficients for aflexible cylinder undergoing vortex-induced vibrationJ.Chinese Journal of Theoretical and Applied Mechanics,2017,49(4):818827.doi:10.6052/0459-1879-16-2639周阳，黄维平，杨斌，等.带螺旋侧板立管两向涡激振动的试验研究J.振动与冲击，2018，37（17）：249 255.doi：10.13465/ki.jvs.2018.17.035ZHOU Yang,HUANG Weiping,YANG Bin,et al.Testsfor vortex-induced vibration in 2-DOF risers with helicalstrakesJ.Journal of Vibration and Shock,2018,37(17):249255.doi:10.13465/ki.jvs.2018.17.03510ZHAO Ming,CHENG Liang,AN Hongwei,et al.Three-dimensional numerical simulation of vortex-induced vi-bration of an elastically mounted rigid circular cylinder insteady currentJ.Journal of Fluids and Structures,2014,50:292311.doi:10.1016/j.jfluidstructs.2014.05.01611唐友刚，樊娟娟，张杰，等.高雷诺数下圆柱顺流向和横向涡激振动分析J.振动与冲击，2013，32（13）：88 92.doi：10.13465/ki.jvs.2013.13.012TANG Yougang,FAN Juanjuan,ZHANG Jie,et al.Inline and trans