采用morison方程基于matlab计算大直径波浪力,小尺寸双柱式桥墩波浪力的数值分析...

小尺寸双柱式桥墩波浪力的数值分析

小尺寸双柱式桥墩波浪力的数值分析

柳春光1，2，李璐璐1，张士博1

(1.大连理工大学建设工程学部，辽宁大连116024；

2.大连理工大学海岸和近海工程国家重点实验室，辽宁大连116024)

摘 要：Morison方程适用于小直径孤立柱的波浪力计算，对于双柱式桥墩的波浪力计算已不能直接采用Morison方程。为了能利用Morison方程求解双柱式桥墩的波浪力，定义了双柱式桥墩前柱、后柱及整体结构的波浪力干扰系数。采用线性波浪理论，利用推板式造波方法模拟线性波，将ANSYS Workbench平台作为FLUENT和ANSYS的数据传递平台，对深水环境下双柱式桥墩受到的波浪力以及波浪场作用下双柱式桥墩的结构变形进行单向流固耦合数值模拟。数值模拟结果表明，双柱式桥墩前、后柱波浪力的干扰系数在一定柱距范围内均随柱距的增大近似的成线性增大，当柱距超过一定值后趋于稳定。由于双柱间存在相位差，双柱墩波浪力的总干扰系数以及双柱墩的结构变形却随着双柱中心距的增大而减小。这些变化曲线为实际工程中深水桥梁双柱式桥墩的波浪荷载计算提供了一种更为简便的算法。

关键词：Morison方程；线性波；双柱式桥墩；流固耦合；结构变形

深水桥梁主要分为三种［1］：第一种是建在江、河上的桥梁；第二种是建在水库中的桥梁；第三种就是现在发展迅速的跨海大桥。对于这些深水桥梁，受到的荷载情况比较复杂，例如波浪与水流作用、地震作用以及风载都是主要的动荷载［2］，这些动荷载主要作用在深水桥梁的基础部分。其中，波浪荷载作为一种永久的主要动荷载，在深水桥梁工程设计中引起了人们极大的注意。

对于单柱结构的波浪荷载计算目前仍广泛采用Morison方程，但对于双柱或多柱式结构的波浪力求解则相当复杂。Bushnell M J［3］在振荡水洞中研究了两个桩和3×3组合桩之间的相互干扰；Sarpkaya T等［4］用振荡水流对不同桩柱的两个桩进行了实验；Chakrabarti S K［5］在波浪水槽中测量了并列2个桩、3个桩和5个桩在不同桩距时的波浪力；王爱群等［6］、徐立论等［7］试验研究了单排两桩、三桩、2×2组合桩3×3组合桩、4×4组合桩的波浪力；俞聿修等［8-9］试验研究了不规则波作用下双桩上的波浪力。本文通过对双柱式桥墩进行波浪作用下的数值模拟［10］可以建立简便、有效的双柱式结构波浪力求解公式，为实际桥梁工程设计提供一定的依据。

本文以ANSYS Workbench平台作为计算流体力学软件FLUENT和有限元结构分析软件ANSYS的数据传递平台［1］，对不同柱距的圆形截面双柱式桥墩结构在波浪场中进行单项流固耦合数值模拟［11］，得到了双柱式桥墩的前柱和后柱以及整个双柱式结构在波浪场中受到的波浪力干扰系数的变化规律。

1 理论分析

1.1 波浪理论

在海洋工程中，波浪理论分为线性波浪理论和非线性波浪理论。对于线性波浪理论我们通常采用微幅波理论；非线性波浪理论主要包括Stokes波理论、孤立波理论、椭圆余弦波理论等，本文主要介绍微幅波理论。

微幅波理论又称Airy波理论，是用势函数来研究波浪的运动。微幅波理论忽略了函数中非线性项的影响，由此可以将问题线性化，得到线性波的势函数［12］：

波面方程：

速度势函数：

根据速度势可以求得波浪水质点的水平速度、垂直速度、水平加速度和垂直加速度。

水平速度：

水平加速度：

式中：H为波高；L为波长；T为周期；ω为频率；d为水深；k为波数，

在有限水深的情况下

thkd，可通过迭代计算得到本文波长L=25 m。

1.2 推板式造波

本文中波浪的生成采用推板式造波方法［13-14］，推板式造波是一种仿物理造波，即用数值方法模拟实验室物理造波机进行造波。推板式造波的原理是将水池的波浪入口边界作为动边界，通过模拟动边界做简谐运动，使水池中的水产生波动，进而形成波浪。推波板做简谐运动的速度［13］为：

式中X0是推波板的冲程，其中在t≤2T时间内，由于刚开始造波，所以对推波板的速度进行了关于时间的光滑处理。

已知推波板的速度方程，运用微幅波理论，并且在结果中忽略掉了只产生局部干扰的非传播项，最后得到了推板式造波的波面方程以及速度势函数［14］如下：

</