卵形消化池基于ALE的流固耦合动力分析.docVIP

卵形消化池基于ALE的流固耦合动力分析 　　摘 要： 本文对卵形消化池结构和流体分别建立有限元模型，应用ALE方法处理流固耦合边界处的网格，分别计算了卵形消化池在天津波和EL Centro波作用下的液体压力和内力响应，研究了卵形消化池的动力特性，结果表明卵形消化池具有良好的抗震性能。 　　关键词：卵形消化池 ALE方法 流固耦合 动力特性 　　中图分类号：C35 文献标识码： A 　　一、前言 　　卵形消化池结构处于正常满水工作状态下，若发生高强度地震时，有可能发生高度的非线性响应，造成破坏。因此，有必要对卵形消化池进行地震分析。 　　本文将利用ADINA，将卵形消化池以及流体作为一个非线性的流固耦合系统来进行分析，采用时程分析法，进行结构动力分析，以期能为卵形消化池的抗震设计提供参考。 　　二、ALE描述 　　处理具有自由液面的流-固动力耦合问题时，采用ALE描述跟踪自由液面和运动边界。使用ALE（Arbitrary Langrange Eluer）坐标系，即任意拉格朗日-欧拉坐标系，该坐标系可以任意速度在空间运动，可处理两相界面的协调及自由表面问题。 　　对于盛在固定区域中的流体，通常采用欧拉描述来进行计算。而对于运动边界的处理，采用ALE描述比较有效 [1]，可使所计算的网格不依赖于流体的运动，使发展的网格能够合理的保持，同时能精确的描述边界。 　　三、工程概况及地震波的选取 　　本文以某一卵形消化池实例为研究对象并进行计算分析。消化池池体采用无粘结预应力混凝土结构，由上部卵形壳体结构和钢筋混凝土基础承台部分组成，单池容积为11，500m3， 卵壳内壁竖向剖面曲线半径21m，池体总高度为42.25m，地面以上部分为壳体部分，高度为31m，最大内径为25.2m，采用强度等级为C30的混凝土，弹性模量为，柏松比为0.16。混凝土密度均取2551千克/m3。壳体截面变厚度，与顶部圆锥交界处厚度为0.45m，与底部基础承台交界处厚度0.9m，由上至下均匀变化。地面以下部分，埋深11.25m，采用C40混凝土，弹性模量为，泊松比为0.2。混凝土抗渗等级均为S8。池壁水平向及垂直向均采用低松驰高强钢绞线施加预应力。图3-1为该消化池池体的竖向剖面图。 　　作用在消化池结构上的静力荷载包括结构自重以及消化池内液体的液压，气压，温度以及预应力荷载。 　　重力加速度取，池内液体重度标准值取，体积模量=2.8e9Pa，粘性系数为。池内工作液面标高为+37.95m。 　　图3-1消化池池体的竖向剖面图 　　根据《建筑抗震设计规范》，采用时程分析法进行地震波作用下的计算。本工程抗震设防烈度为8度，设计基本地震加速度值为0.2g，设计地震分组为第一组，场地土类别为第三类。本文将选取两条典型的地震波记录进行流固耦合计算分析，分别为天津波和el-centro波。考虑中震，将二者最大加速度峰值均调幅至200cm/s2。两种地震波均取南北向水平地震波，施加于模型X轴。 　　天津波持续时间约20秒，南北向水平地面加速度峰值为145.8cm/s2，强度较低。EL Centro波水平地面峰值加速度是341.7cm/s2，震动持续时间大约为30秒。本文将取前20秒时程记录进行分析。在ADINA中，FSI方法只能对位移波进行处理，经过积分将加速度波转换为位移波。 　　四、卵形消化池动力分析 　　对流体与结构分别建立有限元模型，结构上部壳体采用四节点壳单元，下部基础承台采用八节点实体单元。池内液体利用CFD建模，采用27节点流体单元，映射网格。如图4-1、4-2所示： 　　图4-1 结构模型 图4-2 流体模型 　　沿X轴分别输入天津波和EL Centro波的南北向位移荷载，取地震输入方向的内力进行分析。 　　1.液体压力分析 　　尽管两种地震波作用下动水压力的大小不同，但其分布的趋势是基本相同的。在自由液面位置附近静水压力和动水压力都比较小，随着水深的增加而逐渐增大，在地面位置附近达到最大值后，动水压力反而逐渐减小。主要原因是动水压力实际上是池壁附近流体的惯性力造成的，因此与结构的形式有关。而卵形消化池体的中间部分半径较大，上下两端逐渐减小。 　　相比之下，天津波产生的液体压力较大。两条地震波产生的动水压力最大值均发生在壳体与块体交界处附近。相对来说，天津波的频谱范围较为丰富，覆盖了卵形消化池结构的基频，因此液体的动力响应较输入相同加速度峰值的EL Centro波要大，二者的波形也相差较大，但都在静水压力值上下波动。 　　2.内力分析 　　两种波均选用前十秒进行计算，地震波沿笛卡尔坐标系的x轴方向输入。对在二者激励下，处，地面以上部分消化池池壁动力响应过程中的环向力、经向力最大值和最小值进行比较，得出结论：