圆管渡槽抗震计算流体等效简化模型

1、圆管渡槽抗震计算流体等效简化模型：In this paper ， a simplified equivalent model was proposed to analyze the sloshing fluid in the pipeline bridge.The water motion in a pipeline bridge is assumed to be inviscid ，irrotational，compressible ，and linear ( smalldisplacement ) .The transverse sloshing fluid in a pipeline b

2、ridge can be simplified as a fixed rigid mass M0 and a mass-spring system ( M1，K1).According to the equivalence rule ，the actual fluid( computed by the finite elementmodel) and its equivalent model have the same the first sloshing frequency and the same effects on the pipe.The approximate analytical

3、 formulae of the fixed mass M0，mass-spring oscillator (M1，K1)，and its locations in the pipe were acquired by the least-squares and curve fitting algorithms.The simplified equivalent fluid model and fininte element model were used to calculate the first natural frequency and seismic response of a pip

4、eline bridege (fluid-structurecoupling system ).The numericalresults calculated by the two fluid models agreed well ， which confirmed the validity and accuracy of the equivalent model.Therefore ，the simplified model of sloshing fluid can analyze the seismic response of the support structure of pipel

5、ine bridge ， which can provide reference of seismic resistance for the researchers and structural engineers.1 概述我国新修订的水工建筑物抗震设计规范 1 首次增加了 渡槽结构抗震计算与设计的相关条文。 在渡槽抗震计算中， 如何 考虑渡槽内流体晃动以及与结构的相互作用影响是渡槽抗震计 算的关键问题 2 。渡槽抗震分析中，水体可按势流（无旋、无 黏性）考虑 3-4 ，常用的数值分析方法有：（ 1）有限元法，可 分为位移型有限元法 5-6 ，压力型有限元法 7 ，ALE 有限元法 8 ；（2

6、）边界元法 9 ；（ 3）有限体积法及有限差分法 10 等。上部结构可按有限元法建模， 流体与槽身交接面按滑动边界 条件考虑，即：流体与槽体表面的法向位移保持一致，切向位移 不加约束。采用数值方法可满足渡槽体系空间三向地震作用分 析，但流体数值模型以及流一固边界的处理是一件很复杂的工 作，且耗时费钱，不便于实际工程应用。为了简化流体计算，Graham与Rodriguez11首先提出了液 体晃动等效模型的概念， 他们将矩形容器内的液体作用等效为一 个固定质量及一连串的弹簧振子。 Housner12 基于一个物理直 观，提出了一个更为简单的等效模型， 这个模型在土木与水利工 程界得到了广泛应用。最

7、近 Li 与 Wang 1 3针对前述经典解答的 缺陷，提出了一个补充的精确解答，完善了经典解答。李遇春与来明14 在文献 13 的基础上，对固定（脉冲）质量及其位置精 确计算公式进行了数学拟合简化， 提出了表达简单且计算精度高 的建议公式。对于其它截面形状的渡槽，如U 形、梯形等截面的渡槽，Li、Di及Gong15提出了一种半解析/半数值的方法，得 到流体等效模型的近似解。 将渡槽内的流体按等效的简化模型替 代后，抗震计算中的流 - 固耦合动力学问题得到了大大的的简化， 计算精度可满足工程设计要求。 我国新的修订规范给出了矩形与“ U型渡槽中水体的（等效）简化计算方法，然而对于另一种 比较常

8、见的圆形渡槽，规范并未给出其流体简化的计算方法。 圆形（截面）渡槽为封闭的结构，这种结构形式可以避免输水过 程中的蒸发与污染， 这种管线结构在国外的调水工程得到了较为 广泛应用 16 ，比较著名的有法国的 Saint Bachi 渡槽17 等。 圆形（特别是大管径） 渡槽结构在未来的调水工程具有广泛应用 前景。目前关于圆形渡槽（输水管线桥）的地震反应分析研究并 不多见，圆管内晃动流体的（等效）简化模型研究尚未见到文献 报道。本文基于已有的研究成果， 建立圆管内晃动流体的等效力 学模型， 从而建立简化的圆管渡槽抗震计算方法， 为圆管渡槽抗 震设计与研究提供参考。2 圆管内晃动流体的等效力学模型图

9、 1 （ a ）所示的单位长度内圆管内充有静止深度为h 的水体，其中 R 为圆形管道内径，在构建晃动流体等效模型时，基于 以下假定：（ 1）管内的水体为可压缩、无黏、无旋流体的小幅晃动流体，水体的密度及弹性模量分别取定为：p L=1 000 kg/m3, Ev=2.067X 109 Pa；（ 2）管道壁微小变形对流体晃动的影响忽 略不计；（ 3）仅考虑流体晃动的一阶晃动模态，高阶流体晃动 模态忽略不计， 对于一般的流体晃动问题， 其一阶晃动模态起主 要作用，计算中仅考虑一阶晃动模态的影响就可获得较好的计算 精度。于是圆管内晃动的水体可以简化为图1 （b）所示的固定质量（M0和弹簧-质量（M1,

10、 K1）组成的动力系统，其中 M0和 hO分别为脉冲（固定）质量及其作用高度，M1和hl分别为一阶对流质量及其作用高度， K1 为弹簧刚度。根据等效原则：（ 1）流体的原始系统与等效系统具有相同 的自然晃动（振动）频率；（ 2）流体的原始系统与等效系统在 任意水平动力加速度的作用下， 对整个管道的动反力 （包括合力 与合力矩） 相等。采用文献 15 的方法， 借助有限元程序 （ ANSYS code） 15 ， 18 进行计算分析（关于槽内有限元模拟的细节可参 见文献1 5 ，为节省篇幅，本文不再赘述），可以得到上述等效 力学模型的拟合表达式为式中：1为管内流体一阶晃动圆频率；g为重力加速度；

11、 D=2R为管内直径；M为管内（单位长度内）流体总质量。由于管 道内流体动反力为一个汇交力系，这个力系汇交于圆管的中心 点，所以等效模型的质量 M0与M1的位置正好通过圆心，即有： h0=h仁艮图2表示了有限元与式（1）至式（3）结果的比较， 结果显示拟合公式（ 1）至公式（ 3）具有很好的计算精度，可满足工程计算的要求。3 数值算例设有一跨圆管渡槽如图3 （a）所示，其支撑结构为排架形 式，结构尺寸见图3（ a），其中尺寸单位为 mm设管道内径R 为3 000 mm，渡槽跨度为24 m，管道内充液深度为 3 600 mm， 管壁厚为300 mm管道结构钢筋混凝土材料密度为2 500 kg/m

12、3,杨氏模量3.0 X 1010 Pa，泊松比为0.167。图3 （b）所示为圆 管渡槽有限元计算模型， 其中水体采用 Fluid80 单元，水体密度 为 1 000 kg/m3 ，水体与管道壁之间的界面视为滑移边界条件， 即流体和管道壁之间的法向相对位移强制为零， 而释放二者之间 的切向相对位移， 约束水体两端面纵向位移。 管道壁及支撑结构 分别采用Shell181和Beam4单元。图3 （c）所示为圆管渡槽等 效力学计算模型， 其中管道内全部水体沿长度方向分成 24等分， 每一等分（长度为1 m）以固定质量（MQ h0）和质量-弹簧（M1, K1，h1）系统替代，其中固定质量采用Mass2

13、1单元，弹簧采用Combi n14 单元。由拟合简化计算公式（1）得到流体一阶晃动频率为 宀12， 由公式（2）和（3）分别得到M1/M和M0/M以及相应的一阶对 流质量M1和脉冲质量MQ由公式（4）得到弹簧刚度K1，同时 由公式（ 5）分别得到对流质量作用高度 h1 和脉冲质量作用高度 h0。管线桥数值算例中具体计算参数见表 1。首先对包含流体的结构系统进行模态分析，两个系统（图 3(b)与图3 (c)得到的第一阶自振频率分别为 f1=2.481Hz 和 f1=2.486Hz ，两者的相对误差仅为 0.201%，具有高度的一致 性。其次对包含流体的结构系统进行地震反应分析， 地震输入波 选取

14、El-Centro (N-S)波，加速度峰值调整为 0.350 m/s2，地 震波输入时长为10 s。图4为盖梁和连梁端部处水平位移时程 响应曲线，等效模型和有限元模型计算得到的反应曲线吻合很 好。表 2 表示了盖梁与连梁端部水平位移最大值， 两者最大相对 误差盖梁为 4.666%，连梁为 6.211%，等效模型的计算结果满足 工程精度的要求。图 5 为盖梁和连梁端部弯矩时程响应曲线， 等效模型和有限 元模型计算得到的曲线吻合很好。 表3显示了盖梁与连梁端部弯 矩最大值，两者最大误差：盖梁为 6.222%，连梁为 6.221%，计 算结果误差满足工程精度的要求。由于排架支撑结构刚度较小，在地震作用(或风荷载作用) 下容易发生破坏，地震作用(或风荷载)可能成为支撑结构的控 制性荷载， 所以结构的地震反应分析对于排架支撑结构而言尤为 重要。需要说明的是： 管道截面动弯矩及动剪力是抗震设计所需的 荷载，等效简化模型实际上就是动水压力的积分结果， 在抗