一种单跨管道横流向涡激振动的预报方法技术

本发明专利技术公开了一种单跨管道横流向涡激振动的预报方法，该方法包括建立悬跨段管道结构与外界流场的流固耦合作用模型、建立跨肩处管道结构与海床土体的管土耦合作用模型、确定尾流振子模型参数、确定土体支撑刚度和土体阻尼系数等步骤。其优点是：采用尾流振子模型模拟悬跨段的流耦合作用并以试验数据对该模型参数进行标定，考虑了跨肩的管土耦合作用并根据管道运动的方向和位移确定海床土体对运动管道的支持刚度和土体阻尼系数，同时根据管道的运动情况对嵌入土体的最大深度进行实时的更新，很好的模拟了跨肩处管道结构与海床土体之间的相互作用，为海底悬跨管道涡激振动的研究、防范和治理提供了可靠的理论依据。

【技术实现步骤摘要】

【技术保护点】
一种单跨管道横流向涡激振动的预报方法，其特征是，该方法包括以下步骤：(1)建立悬跨段管道结构与外界流场的流固耦合作用模型和跨肩处管道结构与海床土体的管土耦合作用模型：a)悬跨段管道结构与外界流场的流固耦合作用模型：my··+(cs+cw)y·-Ty′′+EIy′′′′=12ρDV2·CL]]>q··+ϵωs(q2-q)q·+ωs2q=ADy··]]>b)跨肩处管道结构与海床土体的管土耦合作用模型：my··+(cs+csoil)y·-Ty′′+EIy′′′′+ksoily=0]]>式中：—管道横向加速度；—管道横向速度；y—管道横向位移；y″—管道横向位移y对管道轴向位移x的二阶偏导；y″″—管道横向位移y对管道轴向位移x的四阶偏导；m—质量项，包括管道结构质量，管内流体质量和管外流场附加质量ma；ma—附加质量，ma＝CaπρD2/4，Ca为附加质量系数；cs—管道结构阻尼系数，cs＝2mωnζ，ζ为结构阻尼比，ωn为管道结构固有圆频率；cw—水动力阻尼系数，cw＝CDρDV/2,CD为平均拖曳力系数；csoil—土体阻尼系数；T—管道受到的轴向拉力；EI—管道结构弯曲刚度；ksoil—土体支撑刚度；ρ—外界流体密度；D—管道直径；V—外界来流速度；CL—瞬时涡激升力系数；—无量纲涡激升力系数对时间t的二阶导数；—无量纲涡激升力系数对时间t的一阶导数；q—无量纲涡激升力系数，CL0为管道静止时的涡激升力系数；ε—尾流振子模型参数；A—尾流振子模型参数；ωs—漩涡脱落频率，ωs＝2πStV/D，St为斯特罗哈数；此为管道微段振动的控制方程，采用有限元法或者有限差分法以连接结点表示整个连续管道，对与结点相关的微分方程在时域上进行求解，得到单跨管道横流向涡激振动的响应；(2)根据试验数据和海床土体性质确定尾流振子模型参数、土体支撑刚度和土体阻尼系数：a)确定尾流振子模型参数：根据试验数据对尾流振子模型中的参数ε和A进行标定：定义约化速度Vr：Vr=2πVωnD]]>当Vr>5,A＝10，当0≤Vr<5，A＝4；参数ε和A满足关系式：CL02(SG+π3St2γ)1+AϵCL04(SG+π3St2γ)=1.12e-1.05SG]]>式中：SG—质量阻尼参数，γ—迟滞系数，γ=CDπ2St.]]>b)确定土体支撑刚度和土体阻尼系数：根据管道的运动情况，对管道嵌入土的最大体深度进行实时更新：如果y>y1,则y1＝y；如果y≤y1,则y1＝y1；式中：y—管道运动位移；y1—管道嵌入土体的最大深度；由管道嵌入土体的最大深度确定最大土体作用力、最大土体吸附力作用点和土体作用力消失点：P1=a(yD)b(Su0+Sugy1)D]]>y2=y1-(1-υ)P1k01+φυ-φ]]>式中：P1—管道嵌入土体的最大深度所对应的土体作用力；y1—嵌入土体最大深度；y2—最大土体吸附力作用点；y3—土体作用力消失点；a—与管沟形状和管道结构表面粗糙度有关的参数；b—与管沟形状和管道结构表面粗糙度有关的参数；Su0—海底泥面土体的不排水抗剪强度；Sug—海底土体不排水抗剪强度垂向增量；υ—海床土体性质参数；k0—海床土体性质参数；φ—海床土体性质参数；海床土体对管道的支撑刚度与管道的运动方向和位置有关：当管道远离土体且管道位置处于y1和y2之间时，土体支撑刚度为：ksoil=(υ-φ)k01+υ;]]>当管道远离土体且管道位置处于y2和y3之间时，土体支撑刚度为：ksoil=-φk0(υ-φ)ψ(1+υ)(1+φ);]]>当管道压缩土体且管道位置处于y1和y3之间时，土体支撑刚度为：ksoil=k0(υ-φ)(1+ψ)(1+υ)(1+φ);]]>其余状态下，土体支撑刚度都为零；土体阻尼系数可以根据下式确定：csoil=ksoilωEDissipated2πEElastic;]]>式中：csoil—土体阻尼系数；ksoil—线性化的土体刚度；ω—响应模态的圆频率；...

【技术特征摘要】

【专利技术属性】
技术研发人员：徐万海，谢武德，袁曦林，崔文瀚，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12