南水北调中线漕河渡槽结构地震动力特性研讨.pdfVIP

第九届全国冲击动力学学术会议论文集 南水北调中线漕河渡槽结构地震动力特性研究 彭辉 雷进生 (三峡大学土木水电学院，宜昌443002) 摘要：基于考虑槽身与槽内在地震荷载下水体流固耦合问题，采用迦辽金法推导了具有压缩性流体与固体 流同耦合问题的动力有限元分析模型，并计算了南水北调中线漕河渡槽结构的地震反应。计算结果表明， 随着水位的上升，渡槽结构地震反应增大，说明槽身与水体的相互耦合对渡槽结构地震反应影响显著，槽 身本身应力、变形也会发生显著的改变。计算成果有利于指导施工和运行期安全监测与管理。 关键词：渡槽，流固耦合，地震动力特性，有限元 1 引言 工程领域中，经常会遇到处于流场中的结构物与流体的相互作用问题，即流同耦合问题。对流固 出的…。郑哲敏研究员在20世纪50年代对平板及悬臂梁与水的耦合振动进行了研究乜’引。在20世 纪60～70年代，国外对流固耦合问题的研究比较活跃，研究对象主要是水丁结构及船舶等在波浪力 的作用下与水的耦合振动i4’，英国学者R．E．D．Bishop等创立了水弹性理论怕7。。现阶段对流【司耦合 系统进行分析有两种方法：一是解析一数值法，即对结构采用有限元离散，对流体则用近似解析描述， 法。对线性流同耦合系统(线弹性结构和理想不可压缩流体)，有限元方法通常有两类：一种是结构和 流体以位移矢量为场变量的位移一位移格式旧’10’，另一种是结构位移矢量和流体场变量的混合型格式， 如位移一压力格式儿川，位移一位移势格式和位移一速度势格式¨引。在有限元方法中，对于不可压缩流体， 最常用的是用附涟水质量解耦【1引。 在考虑流体压缩性时，本文用迦辽金法建立渡槽与水体之间的流固耦合有限元方程。研究结果 表明，用数值方法能够较好地模拟水体与大型挡水结构的耦合作用及地震响应问题。 2计算原理 流固耦合振动问题是一个相当复杂、涉及结构动力学和流体动力学的跨学科的问题。通常可根 据不同类型的流【司耦合问题的具体情况作一些假设，以简化数学模犁，提高计算效率。在文中，既考虑 大体积结构的流固耦合振动，也考虑流体发生小波动时与大体积结构的相互影响。在实际流同耦合计 算中，为了求解方便，尽管可假设流体是无粘、无旋的理想流体，在大多数情况下仍不能求得解析解， 通常只能够用数值解。只有极个别的很有规则形状结构才能用特殊函数，如圆柱体用Bessel函数， 方程)和流体的连续方程综合考虑，经简化后，得出了均匀、无粘、可压缩性流体，在小波动条件下， 采用流【司耦合分析的位移一压力(甜，，p)格式导出有限元方程。 2．1流固耦合系统的动力学模型和基本方程及边界条件”引 假设流体为无粘、可压缩和小扰动的，并假定流体自由液面为小波动(渡槽所在地地震烈度小)。 固体则考虑为线弹性的(冈渡槽结构工作时不允许混凝土开裂，分析中不考虑材料非线性)，其中， K和以分别表示固体域和流体域，鼠代表流同交界面，代表流体刚性固定面界面，S，代表流 体自由表面界面，瓯代表固体位移边界，代表固体的力边界，n，为流体边界单位外法线向鼋， 433 第九届全国冲击动力学学术会议论文集 n，为固体边界单位外法线向量。在流同交界面上任意一点处，n，和n‘厂的方向相反： (1)流体域(巧) 流体场方程为： p埘一与p：o ％ 其中p为流体压力，e。为流体中的卢速。 (2)流体的边界条件 刚性固定边界()上：．兰：o (2) 自由液面(s，)上： 孚+土p：o (3) ’dZ g (3)固体域(K) 固体场方程： 吒√+Z=麒玩 (4) 其中％为固体应力分量，“，为固体位移分量，Z为固体体积力分量，以为I司体质量密度。 (4)固体边界条件 力边界条件(昆)：％％=Z