储液容器内液体晃荡的非线性动力学分析 numerical analysis on liquid sloshing in storage container by nonlinear dynamics method.pdfVIP

第34卷第1期 爆炸与冲击 V01．34，No．1 2014年1月 EXPLOSIONANDSHOCKWAVES Jan．，2014 文章编号：1001—1455(2014)01—0087—06 储液容器内液体晃荡的非线性动力学分析。 李文盛1，赵友清1，贾善坡1’2”，王凯1，谭继可1 (1．长江大学城市建设学院，湖北荆州434023； 2．山东大学岩土与结构工程研究中心，山东济南250061； 3．中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室，湖北武汉430071) 摘要：基于非线性波动理论模型，求解储液容器内液体晃动的固有频率、模态及动力学响应问题。流体 使用Us-“。状态方程，利用ABAQUS软件的自适应网格技术，建立储液容器液体晃动数学模型，通过施加水 平简谐激励得到液体晃动的固有频率和模态，并与解析解对比，验证了该方法的准确性与可行性。然后，分析 了矩形储液容器在多种激励作用下液体非线性晃动响应特性。 关键词：流体力学；模态；非线性动力学；液体晃动；充液容器；地震响应 中图分类号：0351．3 国标学科代码：1302511文献标志码：A 液体晃荡问题广泛存在石油、核工程、船舶、化工等领域，具有非常复杂的流体运动现象，呈现很强 的非线性和随机性。当外界激励频率接近储液系统的固有频率时，内部储存液体晃动将变得更剧烈并 导致容器破坏，不仅造成重大经济损失及人员伤亡，还造成更严重的次生灾害。因此，对液体晃荡和控 制晃荡的研究引起了广泛关注u‘2j。 目前，对液体晃荡问题的数值研究方法主要有MAC法、VOF法、有限单元法(混合插值FEM、分 方法研究了带有自由液面的大幅晃动问题，但这种方法计算存储量大、稳定性差，V．Armenio[41提出了 永涛等Ⅲ基于VOF法得到液体晃荡运动的数值计算模型；周宏等凹1利用任意拉格朗日欧拉法(ALE) 晃动的理想流体，采用边界元方法研究了三维液体的晃动问题；C．Z．Wang等口叩基于流体速度势理论 对流体非线性晃荡问题进行了有限元分析。 本文中，假设矩形储液容器一个方向的尺寸远大于其他两个方向的尺寸，将三维液体晃荡问题简化 为二维大变形平面应变问题，基于非线性波动理论模型，运用有限元软件ABAQus，流体采用Us-U。状 态方程描述，运用自适应网格技术，建立储液容器液体晃动数学模型，通过对其施加水平简谐激励得到 液体晃动的固有频率和模态，并将其与解析解对比，验证方法的准确性与可行性。本文中考虑的是整个 储液系统在激励作用下的非线性晃动响应特性，并非单纯分析流体。 1 非线性波动理论模型 使用流体动力学非线性波动理论求解储液系统流体自由液面大幅波动问题，将三维模型简化为二 维几何非线性平面应变问题，本文中运用有限元软件ABAQUS[1“，采用状态方程进行求解。在非线性 波动理论模型有限元网格中，流体部分采用4节点平面应变单元CPE4R，每个节点具有2个位移自由 度；刚体部分采用R2D2单元；刚体与流体相互作用使用法向接触，其中刚体表面为主面，流体表面为从 面；储液容器壁高足够大以保证晃动时液体没有溢出。 *收稿日期：2012—06—08；修回日期：2012—10—15 基金项目：国家自然科学基金项目；湖北省教育厅科学研究计划项目( 作者简介：李文盛(1968一)，男，副教授，硕士生导师。 万方数据 88 爆 炸 与 冲 击 第34卷 为了得到稳定的数值解，假设流体材料为几乎不可压缩的，且黏性很小。流体体积响应通过Us-u。 线性状态方程和牛顿黏性剪切模型建模。在储液系统模型中，为确定材料体积强度和压力提供了一个 粒子速度U。的线性关系，有： P—Kf (1) 式中：f一1--po／p，po是扰动