矩形渡槽大幅晃动有限元模拟.pdfVIP

矩形渡槽大幅晃动有限元模拟 褚嘉怡 河海大学土木工程系，南京（210098 ） E-mail：zhujiayiae@ 摘 要：本文应用大型通用有限元软件ADINA探讨了计入流体的粘性、大幅（非线性）晃 动影响，不可压缩粘性流体与矩形渡槽弹性体的二维动力耦合问题，分析流体的大幅晃动对 矩形渡槽结构所带来的动力影响问题;在水平谐波激励方式下,针对不同的激励振幅和激励频 率，研究渡槽结构和水体的振动反应包括渡槽结构的倾覆力和倾覆力矩，水体自由液面角点 的波高和水体底部角点的动水压力的反应分析。数值分析表明,渡槽结构中水体的晃动具有 明显的非线性特征,且随着激励振幅的增加,水体晃动的非线性增强. 关键词：流体大幅晃动；流固耦合；矩形渡槽； ALE ；ADINA 1. 引言 地震一般指地壳的天然震动，地震发生时产生的地震波引起对地面建筑物的破坏，导致 人员伤亡，造成了地震灾害。地震对建筑物的破坏，主要是由地震力通过地震波起作用的， 即纵波地震力使建筑物上下颠簸，引起建筑物的纵向结构松动，随后横波地震力再使建筑物 发生水平晃动，引起横向结构损坏。当先颠后晃的地震力超过建筑物的承受力时，在几秒钟 内就能使建筑物遭受破坏。 我国地处欧亚板块的东南部，受环太平洋地震带和欧亚地震带的影响，是个多地震的国 家。而我国制订的用来优化配置水资源、解决北方地区缺水问题的重大战略措施——兴建南 水北调工程，但仅中线工程其中大部分施工地点都是位于地震烈度为 7 度以上的地区，有的 还在地震高发区。因此作为南水北调工程中主要的输水建筑物——渡槽，特别是大型渡槽在 地震作用下的安全性问题更加突出,也一直是水工结构抗震领域关注的焦点。 渡槽中水体在外界激励(地震或风) 作用下会产生大幅晃动. 晃动（Sloshing ）是指两种 或两种以上互不相溶流体（一般为气体或液体）在有限空间内的运动，晃动是一种非常复杂 的流体运动现象，其特点是存在可自由移动的自由表（界）面( 自由表面一方面为流域的边 界，它与控制方程及其它条件组成了问题的定解条件; 自由表面的精确位置事先是未知的， 且其位置每时每刻在变化的，必需作为输运方程的一部分进行求解) ，并且与限制其运动范 围的固体壁面具有强烈的交互作用。由于水体的质量和弹性的变化，结构的动力响应将发生 显著变化，且这种变化往往对结构的安全不利，如将增大结构的动位移和动应力等。晃动是 一种特殊的流动，呈现出很强的非线性和随机性。大幅晃动产生的动水压力会对渡槽结构产 生强烈的冲击，成为渡槽结构承受的主要外力作用之一。所以研究在地震作用下，水体大幅 晃动对渡槽结构的影响是很有现实意义的。若不考虑水体晃动的影响,直接计算渡槽的动力 特性是不准确的,这对渡槽的抗震也可能不利。这种渡槽与槽内流体耦合(简称流固耦合)作用 的研究,是近年来研究人员较为关注的热点和难点问题之一。 目前,考虑流固耦合作用计算渡槽动力特性的方法大体上分为两种:一是线性方法,即把 水体当作附加质量作用于渡槽结构上,与渡槽一起进行动力研究,如附加质量、Housner模型; 二是非线性方法,即研究液体非线性晃动对渡槽的影响。对于小幅晃动问题一般作线性化处 理后可解析法或半解析法近似求解，但是对于具有强非线性的大幅晃动问题的有效分析方法 是进行数值求解，即求解具有自由液面边界的Navier -Stokes方程。用数值方法分析流固耦合 - 1 - 问题时，一般对固体结构采用有限元分析处理，对流体部分经常采用的方法有差分法，有限 体积法，有限元法，边界元法等，因此在流固耦合问题的数值分析中，有差分-有限元，有 限体积-有限元，有限元-有限元，边界元-有限元等不同的处理方法，它们之间的差别主要在 于对流体域进行离散的方法上. 近年来ADINA[1] 已经发展成为一套能全面分析固体，流体和流固耦合问题的系统。 ADINA系统是以通用有限元和有限体积法的理论为基础分析固体，水体流动和流固耦合作 用的。这套系统是将K.J. Bathe在计算力学领域所做的三方面理论研究成果[2] ：固体[2-4，12]， ，