波浪地震耦合作用下悬浮隧道动力响应分析

罗 刚，张玉龙，潘少康，任 毅

(1.长安大学公路学院，西安710064；2.广西新发展交通集团有限公司，南宁530029；3.绍兴文理学院土木工程学院，绍兴312000)

悬浮隧道又称阿基米德桥是挪威学者提出的一种跨越长大水域的新型结构形式[1]，相比于其它结构由于其悬浮在水中一定深度使其具有受外界环境干扰小、造价低、绿色环保等优点[2]。目前悬浮隧道尚处于研究阶段，包括挪威、意大利、中国、韩国、日本等国学者都对其开展了大量的研究[3-8]。

波浪荷载是悬浮隧道服役期间受到的长期荷载。在波浪荷载的计算中，Morison 方程被认为是一种简便的计算方法，项贻强等[9]通过Stokes波浪理论和Morison 方程研究了线性波浪理论计算悬浮隧道波浪力的适用范围。葛斐等[10-11]通过Airy 线性波浪理论和Morison 方程对悬浮隧道锚索的在波浪荷载下的动力响应和海流作用下的涡激振动进行了研究。麦继婷等[12]通过Morison 方程分析了悬浮隧道放置深度、海流速度、波浪周期等因素对水平波浪力的影响。对于更为复杂的非线性波浪力问题则可以采用势流理论进行分析，Paik 等[13]采用边界元法将悬浮隧道简化为三维梁单元，基于绕射理论求解无粘不可压缩的Navier-Stokes方程对悬浮隧道受到的波浪力进行了研究。Wu 等[14]基于势流理论针对悬浮隧道受到的非线性波浪力进行了研究。除了荷载的计算外，波浪荷载下悬浮隧道的响应问题也是研究的重点，Kunisu[15]采用边界元法对悬浮隧道在波浪荷载下的动力响应进行了计算，分析了隧道尺寸和形状对其响应的影响。Remseth 等[16]提出了一种计算波浪荷载下悬浮隧道随机动力响应的方法，分析了波向、阻尼系数、刚度和自振频率的影响。Faggiano等[17]对波浪力作用下悬浮隧道结构的响应进行详细分析，讨论的悬浮隧道不同方向的受力情况以及随深度的变化规律。Dai等[18]以带锚索的浮筒式悬浮隧道为研究对象，分析了长期波浪荷载作用下悬浮隧道的响应问题。Jin等[19]采用室内试验和数值模拟相结合的方式分析了波浪作用下黏滞阻尼系数、浮重比、锚索刚度对悬浮隧道动力响应的影响。

悬浮隧道的修建场地往往在有地震带的海峡上，这给悬浮隧道带来极大的威胁。在此方面Morita 等[20]以格林函数为基础将悬浮隧道锚索简化为支撑弹簧，进行了竖向地震作用下悬浮隧道的响应分析。Lee等[21]在考虑海水的可压缩性、弹性海床对压力波吸收的条件下，研究了海水对SFT系统地震响应的主要影响。数值模拟是研究地震的重要手段，Leira[22]采用大质量法对多点激励下SFT 的动力响应进行了数值模拟，研究了地震波通道效应对悬浮隧道的动力特性的影响。Fogazzi等[23]采用改进的杆单元来模拟锚索，梁单元模拟管体建立悬浮隧道的有限元模型，进行了地震激励下的动力响应分析。上述研究多是针对单一荷载进行讨论，然而工程实际中荷载往往都不是单一出现的，因此有一些学者开展了多荷载耦合下悬浮隧道的响应研究。Jin 和Kim[24]在波浪荷载的基础上考虑了移动荷载，研究了两者耦合作用下悬浮隧道的特性。DiPilato等[25]采用数值模拟的方式对波浪、地震耦合作用下悬浮隧道的动力特性进行了研究。Wu 等[26]基于三角级数法得到了随机地震荷载，建立了地震作用下的水动力荷载计算公式，对悬浮隧道锚索在地震和流体共同作用下的动力特性进行了研究。此外，一些学者也对其他海洋结构进行了多荷载耦合作用的响应分析[27-28]。

对悬浮隧道进行合理的简化，建立合适的理论模型也是其动力分析的重要部分。孙胜男等[29-30]将悬浮隧道锚索简化为简支梁，分析了锚索在轴向激励和随机激励下的动力响应。惠磊等[31]等建立了冲击作用下悬浮隧道的圆柱壳模型，采用等效质量法进行求解，分析了碰撞中的能量转化和冲击点的最大位移。陈健云等[32-33]通过将管体简化为质量块，将锚索简化为欧拉梁建立了悬浮隧道耦合振动模型，考虑了水流作用下悬浮隧道的耦合振动响应。Lin 等[34]通过将管体和锚索简化为耦合连接的欧拉梁系统，建立了悬浮隧道管体-锚索耦合振动模型，研究了移动荷载下悬浮隧道的动力响应。总体而言，目前关于悬浮隧道研究主要集中在单一荷载、特定构件的研究，对于荷载耦合和结构耦合的研究较少，不能较为真实反映悬浮隧道的实际工况。

为了综合考虑波浪地震荷载以及悬浮隧道结构的耦合效应，本文在前人研究的基础上，通过Stokes波浪方程和三角级数法计算波浪荷载和地震荷载，基于D’Alembert 原理建立波浪地震作用下管体-锚索耦合振动模型，通过四阶Runge-Kutta 法数值积分求解微分方程组，对地震、波浪的荷载参数和系统响应进行分析和讨论。

1 荷载模型

1.1 波浪荷载

采用微幅波理论的一阶近似不能满足波浪荷载计算精度要求，应当考虑非线自由表面对波浪荷载的影响。Stokes波浪理论是计算非线性波浪荷载的常用方法，在保证计算精度的情况下，本文选择Stokes三阶波浪方程来进行计算，波面方程、速度势函数和弥散关系为：

在已知波高H、周期T以及水低深度d的情况下，可以先求波长L和a，进一步则可解出Stokes三阶波浪方程的全部方程。将波浪引起的水质点速度式(2)代入Morison 方程即可计算出作用在悬浮隧道上的波浪荷载，在后文中将进行这部分计算。

1.2 地震荷载

为了对悬浮隧道进行地震响应时域分析，必须建立合适的地震加速度时程曲线。一种可行的方法是通过三角级数法计算出平稳过程条件下地震加速度时程曲线[35]：

式中：S(ω)为设计地震功率谱；Δω为圆频率增量，Δω=2π/Td，Td是地震的持续时间；ωk=kΔω；εk为随机相位改变量，其变化范围为0～2π。

考虑到随机地震激励的非平稳特性，根据功率谱包络理论中的非均匀调制函数对式(3)进行修正[36]：

式中：t0和tn分别为控制平稳段的首末时刻；c为衰减常数。

对于地震谱的选择具有较大任意性，结合实际的场地条件和结构特性采用设计地震加速度响应谱，能够使悬浮隧道的抗震分析更加的精细化和合理。目前中国罕有海洋结构相关的抗震设计标准。因此，本文参考国际标准化组织(ISO)制定的海洋结构抗震设计标准[37]，地震设计反应谱需要根据场地分类、峰值加速度和阻尼比三个参数进行选择，具体如图1所示。

图1 地震设计加速度响应谱Fig.1 Seismic design acceleration response spectrum

通过设计地震响应谱则地震的功率谱可以近似通过下式得到[38]：

式中：ζ 为结构的阻尼比，此处取0.05；Sa(ω)为与Sa(T)对应的圆频率响应谱，可由对后者进行T=2π/ω变量替换得到；P为超越概率常数，在此取0.9。

由于式(5)的转换存在误差，需要对人工地震功率谱进行修正。对于每一个频率点处的人工地震功率谱，以计算反应谱和目标反应谱的比值作为验证精度值，当精度值不满足拟合精度值时，按照比值进行功率谱值的调整，调整后再次生成人工地震动，并用新人工地震动计算相应的反应谱，再与目标反应谱比较，进行精度值检验和调整。如此反复，直至达到精度要求。

通过上述方法，可以建立与目标场地和结构相匹配的加速度时程曲线，则作用在结构上的地震荷载可用下式来表示：

式中：Fg为作用在结构上的地震荷载；m为结构的质量。

1.3 流体作用

如果流体在某一点处的速度和加速度已知，则可以通过Morison 方程即式(7)计算出流体的作用力，方程的第一项是流体作用引起的阻尼力，第二项是流体作用引起的惯性力。

式中：ρw海水的密度；D为悬浮隧道的外径；v为水质点速度；w为管体振动的位移函数；CD和CM分别为流体作用的阻尼力系数和附加质量系数，可通过流体试验测定，在缺乏大量试验的条件下CD可近似取1.0，CM可近似取2.0[9]。

已有研究表明：在相同荷载作用下，悬浮隧道水平方向的响应远小于竖直方向[39]，因此本文仅考虑悬浮隧道的竖向受力及响应。波浪和地震荷载会通过流体间接的作用在悬浮隧道系统上，结构自身的运动也会导致流体作用，此处将同时通过Morison 方程计算三者的流体作用。在计算中锚索受到的流体作用与隧道管体形式相同，但考虑到锚索的直径较小且具有较大的安装水深，在此假设锚索不受到波浪力的作用[11]。

式中：u为锚索的位移；di为锚索的直径；xg为地震引起的水质点位移。

2 隧道模型及方程求解

目前在悬浮隧道理论计算方面，主要有弹性地基梁模型(Beam on elastic foundation)、锚索振动模型(Tether vibration model)和管体-锚索耦合振动模型(Tube-tether coupled model)三种[40-41]。其中，弹性地基梁模型将悬浮隧道锚索简化为弹性支撑，将管体简化为弹性支承梁，能够对管体的动力响应进行分析但忽略了锚索的振动。锚索振动模型将悬浮隧道锚索简化为倾斜的梁，将管体简化为质量块或作用在锚索上的简谐激励，能够对锚索的动力响应进行分析，但忽略了管体的振动。管体-锚索耦合振动模型同时考虑了锚索和管体的振动以及耦合效应能够同时对管体和锚索进行动力响应分析，因此本文将选用该模型进行计算。

2.1 锚索振动

如图2(a)，在悬浮隧道的两侧，假设等间距h布置着Ns个相同的锚索，管体在浮力作用下使锚索处于绷紧的状态，此时锚索具有一定的初张力T0。管体发生w的位移后，在与锚索连接的xi处与其发生位移协调，同时锚索自身也会产生的法向振动，因此，除初张力外锚索还承受附加的张力ΔT。

此外，由于锚索自身的重力和柔度较大，往往表现出一定的垂度，使其受力特性发生改变。为了考虑锚索的垂度效应，采用等效刚度法来计算锚索的弹性模量[42]。

2.2 系统微分方程及求解

悬浮隧道管体受到波浪力、流体作用力以及锚索的张力，通过锚索的张力项隧道管体与锚索之间相互耦合，在锚索振动方程的基础上可以写出悬浮隧道管体-锚索系统的耦合振动方程，如式(11)。

式中：下标i代表第i根锚索，mi为单位长度锚索的质量；ui为锚索的位移；ci为锚索的粘滞阻尼系数；Fdi为第i根锚索受到的流体作用力。

图2 悬浮隧道耦合振动系统Fig.2 Coupled vibration system of SFT

式中：mtb为单位长度管体的质量；ctb为管体的粘滞阻尼系数；EtbItb为管体的弯曲刚度；δ 为狄拉克函数。

上述微分方程组可采用Galerkin 法即式(12)进行化简，将偏微分方程组转换为常微分方程组：由于在锚索的振动中第一阶振型占主导地位[43]，在此取锚索的第一阶模态进行简化计算，化简后的结果为：

上述的时变耦合微分方程组无法求得解析解，对于悬浮隧道管体微分方程而言，计算n阶振型需要联立求解n个微分方程。对于悬浮隧道锚索而言，锚索数量Ns对应着需要求解Ns个微分方程。本文计算对耦合微分方程系统进行降阶处理，然后采用四阶龙格库塔法数值积求解，具体步骤如图3所示。

图3 系统微分方程系统求解步骤Fig.3 Stepsof solving system differential equation

参考国内外待建悬浮隧道的参数，在下文的计算中所用参数如表1。

3 计算与分析

3.1 不同地震谱

为了研究不同场地类别和地震作用下悬浮隧道的动力响应，在此根据已经建立的设计地震谱选择了三条与之匹配的地震记录进行对比分析，如表2所示，同时将四种地震的加速的时程曲线如图4(a)所示。

表1 悬浮隧道系统基本参数Table 1 Primary parametersof SFTsystem

表 2天然地震波参数Table 2 Parameters of natural earthquakes

由图4(a)可以看出，基于三角级数法生成的人工地震波，在时域上具有较大的分布密度和加速度峰值，主要地震作用时间在前40 s。而基于设计谱在太平洋地震工程研究中心(Peer)上选择的地震波具有较大的差异性，主要作用时间及加速度幅值均不相同，但能为悬浮隧道的抗震性能进行多方面的检验。图4(b)为本文中锚索位移的计算结果与Wu 等[26]和Su[41]的结果进行比较，在计算的过程选择了相同的地震荷载参数和锚索结构参数。后者采用的计算模型为锚索振动模型，该模型将悬浮隧道管体的作用力简化为作用在锚索上的简谐荷载，主要考虑锚索的振动作用。由于采用了随机生成的地震波，计算出来的结果必然存在一定的差异。尽管本文所用悬浮隧道模型与他们不同，但计算结果较为相似。值得注意的本文在此对比的是3#锚索的位移，2#锚索因为处于管体的跨中位置，同时受到参数共振的影响产生了较大位移，这在非耦合模型中是没有体现的。

通过多点虚拟激励法将人工地震和自然地震荷载施加在悬浮隧道系统上，同时考虑波浪荷载计算后得到图4(c)～图4(d)动力响应结果。由图4(c)可知，在人工地震作用下悬浮隧道管体的峰值位移为0.30 m，出现在11.7 s。人工地震波的持续时间约为40 s，从图中可以看出前40 s，管体震荡较为剧烈。另外从图中也能看出2#锚索位移最大，为2.15 m，出现在9.9 s，其余两对锚索的最大位移分别为0.57和0.27 m。

由图4(d)可以发现，不同地震荷载输入下结构的响应具有较大的差异性，Helena Montana-01地震的特点是峰值加速度大、持续时间短，因此前20 s内管体的振动最为剧烈，但所造成的位移不大。Imperial Valley-02地震与前者相似，但地震持续时间较长，峰值加速度也更大，其管体最大位移可达0.41 m。Humbolt Bay 地震的水平加速度峰值较小，但持续时间较长，管体的最大位移为0.02 m。总体而言，不同地震作用下管体的位移的幅值和频率均有差别，但主要影响因素还是地震的峰值加速度和持续时间。

3.2 地震方向和大小

在结构抗震分析中，一般将地震分为水平地震作用和竖向地震作用来分析。此处通过海洋结构抗震规范计算出竖向地震下悬浮隧道的响应与水平地震下的响应相比较。图5(a)给出了不同方向地震作用下管体跨中位移时程曲线，由图中可知水平地震的最大位移为0.30 m 方向为正，而竖直地震的最大位移为0.23 m 方向为负，可以看出竖向地震产生的管体位移约为水平地震的76.7%。

图5(b)给出了不同地震方向下锚索的跨中位移时程曲线。由图中可以看出无论是水平还是竖向地震作用下2#锚索的位移都最大。水平地震和竖向地震对应的锚索最大位移分别为1.23 m 和0.34 m。可见水平地震会使锚索产生较大的位移响应，其最大位移约为竖向作用下的3.6倍。比较管体和锚索的位移，可以看出地震对于锚索的影响大于管体，这是因为锚索固定在海床上直接受到地震的作用，而管体则是通过锚索间接受到地震作用。图5(c)给出了不同地震方向下地震峰值加速度(PGA)与隧道管体最大位移之间的关系曲线，由于采用三角级数法生成的地震谱具有一定的随机性，系统响应从而也有了随机性使得PGA 与最大位移之间没有形成明显的函数关系，但对数据进行拟合分析后发现，PGA 与最大位移之间近似呈线性关系。

图4 不同地震波对悬浮隧道响应的影响Fig.4 Impact of different seismic waves on responses of SFT

3.3 波浪作用

为了分析波浪作用对地震荷载作用下悬浮隧道系统响应的影响，本小节分别计算了仅地震、仅波浪、地震波浪作用的三种情况。

图5 不同地震方向及大小对悬浮隧道系统的影响Fig.5 Impact of different seismic directions and PGA on responses of SFT system

图6(a)中给出了三种荷载情况下的管体弯矩时程曲线，可以看出地震作用和波浪作用下隧道响应具有不同的作用特点，仅地震、仅波浪、地震波浪共同作用下最大弯矩分别为12.04×108N·m、0.18×108N·m 和15.58×108N·m，在地震作用上考虑波浪后管体弯矩增大了29.4%。图6(b)给出了三种荷载情况下的2#锚索的张力时程曲线，仅地震、仅波浪、地震波浪作用下锚索张力分别为183.27×105N、1.65×105N、261.09×105N，在地震作用上考虑波浪后锚索张力增大了42.5%，由此可以看出对于波浪荷载对悬浮隧道系统的响应具有一定促进作用。此处仅波浪作用下锚索的位移小于仅地震作用。然而在后续计算中发现不同波浪参数对系统响应影响较大，并非恒小于地震作用。

进一步的，图6(c)～图6(d)对波浪荷载的参数进行了讨论，选择了不同波高、波长和周期进行计算，为了避免地震荷载的干扰，在计算中仅考虑了波浪荷载。图6(c)给出了波高、波长与管体最大位移的关系曲线，由图中可以看出管体最大位移与波长、波高直接近似呈线性关系。进一步的，图6(d)给出了不同波浪周期下管体最大位移和波浪频率，由图6(d)可知随波浪周期的增大管体最大位移出现明显的峰值，最大位移达0.74 m。随波浪周期的增大波浪频率呈反比例下降，在波浪周期较小时(小于10 s)，波浪周期与管体1阶频率相近能够产生共振。

图6 不同波浪参数对悬浮隧道响应的影响Fig.6 Impact of different wave parameters on responses of SFT

4 结论

通过Stokes三阶波浪方程及三角级数法计算出波浪和地震荷载，提出了一种计算波浪地震耦合作用下悬浮隧道系统响应的简化分析模型，对荷载参数和系统参数进行讨论得到以下结论：

(1)本文提出的波浪地震耦合作用下悬浮隧道简化分析模型与锚索振动模型在计算锚索位移时具有较好的一致性，但后者不能考虑管体和锚索的耦合振动效应。

(2)地震的方向对于系统响应具有显著影响，在相同峰值加速度情况下水平地震会使系统产生更大位移响应，同期锚索的位移远大于管体。随地震荷载峰值加速度的增加系统的最大位移响应约呈线性增加趋势。

(3)在地震荷载的基础上，考虑波浪荷载后系统的响应出现一定增长。随波高和波长的增大，系统响应整体呈线性增大趋势。此外当波浪周期较小时(小于10 s)与悬浮隧道管体1阶频率相近容易引发系统的共振。