微幅晃动下考虑地基效应的矩形贮液结构动力响应研究

微幅晃动下考虑地基效应的矩形贮液结构动力响应研究程选生，景 伟*，李国亮，李 德，马 亮（兰州理工大学 土木工程学院 甘肃，兰州730050） 摘 要：一般情况下矩形贮液结构通过底板直接坐落在地基上，有必要在考虑地基效应的基础上研究其动力特性。定义了微幅晃动波高限值，建立了矩形贮液结构的简化模型。采用人工边界模拟地基效应，建立了有限元计算模型，对比了两类模型对应的晃动波高以验证简化模型的合理性，基于势流理论在微幅晃动范围内研究了不同边界条件、不同地震波和不同地震加速度幅值下矩形贮液结构的动力响应。结果表明，两类模型得到的晃动波高趋势一致，且最大晃动波高相差较小。在微幅晃动范围，液体最大晃动波高和地震加速度幅值成线性关系，考虑土-结构相互作用后，壁板水平位移、壁板有效应力和液体晃动压力有较大程度的的减小，而地基效应对液体晃动波高影响很小。 关键词：微幅晃动；土-结构相互作用；矩形贮液结构；粘弹性边界；流-固耦合；动力响应 中图分类号：TU357 文献标识码：A DOI:Dynamic response of rectangular liquid-storage structure considering foundation effect under small amplitude sloshingCHENG Xuan-sheng, JING Wei, LI Guo-liang, LI De, MA Liang(School of Civil Engineering, Lanzhou University of Technology, Lanzhou, 730050, China) Abstract: Rectangular liquid-storage structure is often located on the foundation through the bottom plate, it is necessary to study the dynamic response characteristics on the basis of considering the foundation effect. The limit height of small amplitude sloshing is defined, the simplified model of the rectangular liquid storage structure is established. Using artificial boundary to simulate foundation effect, the finite element calculation model is established, comparing sloshing height corresponding to two kinds of models to verify the rationality of the simplified model. Based on the potential flow theory, the dynamic responses of rectangular liquid-storage structure with different boundary conditions, different seismic waves and different seismic acceleration amplitudes are studied. The results show that sloshing height trend of the two kinds of models is consistent, and the difference of maximum sloshing height is small. In the range of small amplitude sloshing, there is a linear relationship between the maximum sloshing height and the amplitude of earthquake acceleration, after considering soil-structure interaction, tank wall horizontal displacement, tank wall effective stress and liquid sloshing pressure have a greater degree of reduction, but foundation effect has little influence on liquid sloshing height. Keywords: small amplitude sloshing; soil-structure interaction; rectangular liquid-storage structure; viscoelastic boundary; fluid-structure interaction; dynamic response基金项目：国家自然科学基金项目（51368039；51478212）；教育部博士点基金（博导类）（20136201110003）；甘肃省科技支撑计划项目（144GKCA032）收稿日期：2015 -10 -第一作者 程选生，男，工学博士，教授，1972年生通信作者 景 伟，男，博士研究生，1988年生钢筋混凝土矩形贮液结构在给排水等生命线工程和石油化工行业有着广泛的应用，该类特种结构在地震作用下失效不仅会影响人民的日常生活和工业企业的正常生产，而且还会引起液体泄漏导致火灾及环境污染等次生灾害，更甚者会威胁人民的生命安全。贮液结构往往通过底板直接坐落在地基上，有必要在考虑地基效应的基础上研究该类结构的动力响应，关于地基效应对结构动力响应的影响，国内外学者做了大量的研究。夏栋舟等1推导了土-结构相互作用体系的耦合阻尼比公式，结果表明考虑土-结构相互作用后能够增大体系的阻尼性能，大大降低上部结构的地震作用；刘洁平等2提出了土-结构相互体系地震反映分析的简化评估方法；陈少林等3提出了一种分析三维土结构动力相互作用的时域直接方法，该方法可同时得到土体和结构的反应；张之颖等4采用振动台试验研究了软弱基础上土-结构相互作用体系的阻尼特性；刘伟庆等5运用振动台试验研究得到土-结构相互作用会使结构的频率减小，能够改变结构体系的阻尼比；易伟建等6指出考虑土-结构相互作用后，基础变形对结构位移响应有很大的影响；刘毅等7采用整体有限元法分析了土-结构动力相互作用(SSDI)下结构的动力特性，得到SSDI能使基础底面峰值加速度增大5%30%，且地基土越软增大幅度越大。贮液结构会由于流-固耦合的影响，在动力响应方面表现出一定的特殊性。王博等8采用集中弹簧模型考虑土-结构相互作用，结果表明土-结构相互对大型渡槽的模态和地震响应有影响；Larkin9研究了考虑土-结构相互作用的储液池地震响应，得到与固支边界相比峰值加速度可增大2.3倍；程选生等10通过引入无量纲参数建立了弹性地基上矩形贮液结构液-固耦合系统的振动方程；Kianoush等11运用有限元法并考虑土结构相互作用研究了地震频谱特性对混凝土矩形储液池动力响应的影响，得出考虑地基效应后，体系的最大脉冲基底剪力和基底弯矩会随着土体刚度的变化增大或减小；Patel等12研究得出了地基参数对土结构相互作用有很大的影响，而且储液深度对结构响应的影响极大；张华等13研究了渡槽流固耦合体系的动力响应，得到SSI效应对渡槽结构横槽向刚度和主应力有影响。课题组已在贮液结构的耦合振动和液动压力14-16、液-固耦合振动试验17、贮液结构自振特性的计算模型18和贮液结构液-固耦合地震响应19等方面进行了大量的工作。 矩形贮液池结构由于流-固耦合作用表现出一定的特殊性，为了全面掌握该类结构在地震作用下的灾变问题，有必要研究考虑地基效应对其动力响应的影响。本文通过在贮液池底部设置粘弹性人工边界来考虑地基效应，定义了微幅晃动波高限值，基于计算流体动力学，采用势流理论研究了矩形贮液池结构在微幅晃动下的动力响应，通过对比不同边界条件的计算结果，得到了地基效应对贮液池结构动力响应的影响规律，有助于为矩形贮液结构的合理设计提供理论依据。1 微幅晃动 液体晃动波高的大小对贮液结构求解方法的选择及动力响应都有很大的影响，因此，合理地确定微幅和大幅晃动的界限对贮液结构流-固耦合研究显得尤为重要。若液体为微幅晃动（图1），则其求解可采用线性势流理论，从而使得计算分析大大简化。因此有必要以速度势函数f为分析对象，导出液体微幅晃动需要满足的条件。hzx图1 微幅晃动Fig.1 Small sloshing 由速度势函数表示的流体运动Bernoulli方程20可表达为式（1）： （1） 由式（1）得到，当自由液面未扰动时，z=0，并令此时作用的大气压强为P0，则C(t)=P0/r。 由自由液面动力学边界条件可知，在自由表面作用着大气压强P=P0，于是在液体自由表面式（1）可简化为： （2） 若要满足线性化，由式（2）得到速度势函数必须满足： （3） 即 （4） 黄德波20利用分离变量法得到了液体线性晃动对应的速度势函数的近似表达式：（5） 将式（5）带入（4）并进行化简得到： （6）式中，a为波幅，h为贮液高度。 表达式（6）需在流体域的各个位置都满足，则可令z=0，进一步由式（6）得到： （7） 对于微幅晃动，引入高阶无穷小量e，考虑到表达式（7）的任意性，且液体晃动以一阶为主，则可取k=1，由此得到微幅晃动波幅a需要满足的条件： （8）式中，e为高阶无穷小量（10-1）。2 计算模型及动力响应求解在矩形贮液池底部使用人工粘弹性边界模拟土-结构相互作用，将贮液池中的液体按三质点体系模型进行简化，即液体总体分为3部分：刚性质量m0、对流质量mc和脉冲质量mi，由于钢筋混凝土结构本身质量较大，在动力分析中还需要考虑混凝土质量m，由于在模型中假定液体刚性质量和贮液池一起做刚性运动，为了简化并减少自由度，因此可以将钢筋混凝土贮液池的质量m和液体刚性质量m0进行叠加，得到总刚性质量m+m0，最终得到考虑地基效应钢筋混凝土矩形贮液结构的简化力学模型如图2所示21-22。Kc/2K0C0u0(t)ackucuiu0mim0+mCi/2Kc/2Ci/2Kc/2Kc/2Cc/2Cc/2mc图2 简化模型Fig.2 Simplified model根据Hamilton原理，得到该系统的能量平衡方程如式（9）： （9）式中：T为系统的动能；V为系统势能；dWc为系统阻尼耗散能。 （10） （11）（12） 将式（10）、（11）和（12）带入式（9）得到系统的运动方程为：（13）式中：Cc、Ci、C0分别为对流质量、脉冲质量、刚性质量对应的阻尼；Kc、Ki、K0分别为对流质量、脉冲质量和刚性质量对应的刚度；c和k为考虑地基效应的附加阻尼和刚度；uc、ui和u0分别为对流质量、脉冲质量、刚性质量对应的位移；为地面加速度。 在ti时刻位移给定的前提下，可采用Wilson-q方法得到ti+qDt时刻的位移，当q足够大（1.37）时，Wilson-q是无条件稳定的16。采用该种方法求解时，存在逐步积分表达式（14）： （14） 式中等效刚度矩阵和等效荷载向量可表示为： （15） （16）式中：、和分别为体系的质量矩阵、阻尼矩阵和刚度矩阵。 （17） （18） （19） 由式（14）得到u(ti+qDt)，由（17）、（18）和（19）得到、和。3 数值算例3.1 计算模型钢筋混凝土矩形贮液结构尺寸为6m6m4.8m，壁板厚度为0.3m，储存液体高度为2.4m。对于矩形贮液结构，可取e为10-1，由此得到微幅晃动的限值为0.24m。考虑流-固耦合，在贮液池底部设置人工粘弹性边界模拟土-结构相互作用，采用实体单元模拟壁板，采用势流单元模拟液体，液体密度为1000kg/m-3，体积模量为2.3109Pa。运用ADINA建立计算模型如图3所示，混凝土材料参数见表1，粘弹性边界参数见表2。（a）固支边界（b）人工粘弹性边界图3 计算模型Fig.3 Computational model表1 混凝土材料参数Tab.1 Concrete material parameters密度r/(kg/m-3)切线模量E/(Pa)抗拉强度s/（Pa）抗拉屈服强度s/（Pa）25003.3810102.391062.2106抗压强度s/（Pa）压应变e极限抗压强度s/（Pa）极限压应变e-40106-0.002-26.8106-0.004表2 粘弹性边界参数Tab.2 Viscoelastic boundary parameters参数粘弹性边界I粘弹性边界II粘弹性边界IIICBN(Ns/m)1.201073.84781068.21106KBN(N/m)7.441072.02421075.09107CBT(Ns/m)7.581062.05681065.18106KBT(N/m)3.721071.01211072.541073.2 模型验证简化模型动力方程的质量参数m0、mc和mi，阻尼参数Cc、Ci和C0以及刚度参数Kc、Ki和K0参考文献23计算得到，具体参数见表3，地基参数k和c分别为7.44107 N/m和1.20107 Ns/m。运用MATLAB编程求解简化模型对应的动力方程，简化模型对应的晃动波高由式（20）得到22： （20）式中，b为平行于地震作用方向贮液结构的边长，g为重力加速度。表3 简化模型参数Tab.3 Parameters of simplified model质量（kg/m3）mcmim06.60431041.80981042.2595104阻尼（Ns/m）CcCiC01.37891031.65711057.4415106刚度（N/m）KcKiK02.87891051.51721082.3378108运用数值仿真工具ADINA求解计算模型，有限元计算模型粘弹性边界取参数I，假定混凝土为线弹性材料，弹性模量为3.381010 Pa，密度为2500 kg/m3，泊松比为0.2。选取Chichi波为地震作用，PGA为0.4 g。对比简化模型和计算模型对应的晃动波高，以便验证模型的合理性，两类模型晃动波高的计算结果对比见图4。图4 计算结果对比Fig.4 Comparison of calculation results由图4得到，三质点简化模型和有限元计算模型在计算液体晃动波高时得到的趋势基本一致，同时，两者得到的晃动波高最大值分别为0.235 m和0.211 m，计算结果相差较小，从而说明了考虑地基效应混凝土矩形贮液结构三质点简化模型的合理性。3.3 边界条件对贮液结构动力响应的影响在临界波高确定的基础上，基于势流理论研究液体微幅晃动问题。为了研究地基效应对矩形贮液池结构动力响应的影响规律，本文考虑了四种边界条件：（1）底部固支；（2）粘弹性边界I；（3）粘弹性边界II；（4）粘弹性边界III。所选取的用于时程分析的四条地震波Chichi波、Kobe波、Trinidad波和Friuli波对应的反应谱曲线见图5。图5 反应谱Fig.5 Response spectrum基于增量动力分析原理，对选取的四条地震波调幅，使其幅值依次为0.10g、0.15g、0.20g、0.25g、0.30g、0.35g和0.40g。通过对计算模型输入Chichi波、Kobe波、Trinidad波和Friuli波，得到矩形贮液池结构在液体微幅晃动下的动力响应，最大晃动波高和地震幅值的关系见图6，最大晃动波高对应的节点液体晃动时程曲线见图7，池壁水平位移时程曲线见图8，最大液动压力和地震幅值的关系见图9，池壁最大有效应力和地震幅值的关系见图10。限于篇幅，图7图10只列出了Chichi波和Kobe波对应的计算结果。（a）固支边界（b）粘弹性边界I（c）粘弹性边界II（d）粘弹性边界III图6 最大晃动波高与地震幅值关系Fig.6 Relationship between the maximum sloshing height and seismic amplitude由图6可以看出，在液体微幅晃动下，四种边界条件下液体最大晃动波高与地震加速度幅值成线性关系，在地震三要素中持续时间和加速度幅值相等的情况下，Chichi波、Kobe波和Friuli波三条地震波对应的液体最大晃动波高和Trinidad波对应的值有较大的差别，主要原因在于地震动频谱特性对液体晃动波高有很大的影响。该贮液结构的周期在3s左右，由图5地震波反应谱曲线可以看出，3s左右Chichi波、Kobe波和Friuli对应的反应谱曲线非常接近，即结构第一周期对应的加速度反应谱Sa(T1)相差很小，而Trinidad波对应的加速度反应谱最小，使得Trinidad波对应的最大晃动波高在同一地震波幅值下与其它三条波相比其值最小。（a）Chichi波（0.4g）（b）Kobe波（0.4g）图7 四种边界下晃动波高时程曲线Fig.7 Time history curves of sloshing height under four kinds of boundary conditions由图7得到在Chichi波和Kobe波作用下，钢筋混凝土矩形贮液池的液体微幅晃动波高在四种边界条件下相差非常小，即就贮液池液体晃动波高来说，是否考虑地基效应以及地基参数取值对其计算结果的影响很小。（a） Chichi波（0.4g）（b） Kobe波（0.4g）图8 四种边界下壁板水平位移时程曲线Fig.8 Time history curves of wall horizontal displacement under four kinds of boundary conditions由图8得到支座边界条件对池壁水平位移的影响不容忽视，贮液池固支的池壁水平位移计算结果明显大于粘弹性边界，最大水平位移在固支和粘弹性边界下，计算结果最大能相差好几倍。考虑土-结构相互作用时，地基的刚度和阻尼参数取值对池壁水平位移有影响，由于地基阻尼等的影响，使得结构的响应有很大程度的减小。考虑土-结构相互作用后，当地基参数取值合理，在允许的误差范围时，不同参数取值引起的结构动力响应相差较小，可以满足工程的分析精度要求。（a） Chichi波（b） Kobe波图9 最大液动压力与地震幅值关系Fig.9 Relationship between the maximum hydraulic pressure and seismic amplitude由图9得到在四种边界条件下最大液体晃动压力随着地震加速度幅值的增大而增大，考虑地基效应后，晃动压力值相比于固支边界有较大的减小，最大可相差9%左右。另外，考虑地基效应后，地基参数的不同取值对液动压力有影响，但不是很大。（a） Chichi波（b）Kobe波图10 池壁有效应力与地震幅值关系Fig.10 Relationship between effective stress of tank wall and seismic amplitude由图10得到，在Chichi波和Kobe波作用下，四种边界条件对应的池壁有效应力最大值随着地震动幅值的增大而增大，考虑地基效应对池壁有效应力有很大的影响，固支和粘弹性边界相比，有效应力最大值可相差1.5倍左右。另外，地基参数的取值对有效应力计算结果有影响，但总体来看，其影响比较小，在某些地震加速度幅值下，四种粘弹性边界计算结果基本相等。4 结论本文基于势流理论，在结构几何尺寸及水深等参数确定的情况下，定义了微幅晃动临界波高限值，通过在贮液池底部设置人工粘弹性边界来考虑地基效应，研究了矩形贮液结构在微幅晃动下的动力响应，主要结论如下：（1）在微幅晃动范围内，地震作用下液体晃动最大波高和地震加速度幅值基本成线性关系，该结论有助于对贮液结构在液体微幅晃动下，同一地震波而加速度幅值不同的晃动波高进行预测。（2）四种边界条件下，相同位置节点的液体晃动波高时程曲线相差很小，由此得到考虑地基效应与否对晃动波高的计算结果影响不是很大。（3）考虑地基效应后池壁水平位移、液体最大晃动压力和池壁有效应力都有较大程度的减小，经过合理设计，可采用较薄的壁板从而节约资源。四种边界条件下，地震幅值越大，晃动压力和有效应力也越大。（4）考虑土-结构相互作用后，当地基参数取值合理，在允许的误差范围时，不同参数取值引起的结构动力响应相差较小，可以满足工程的分析精度要求。参考文献1 夏栋舟, 何益斌, 刘建华. 土-结构动力相互作用体系阻尼及地震反应分析J. 岩土力学, 2009, 30(10): 2923-2928. Xia Dong-zhou, He Yi-bin, Liu Jian-hua. Study of damping property and seismic action effect for soil-structure dynamic interaction systemJ. Rock and Soil Mechanics, 2009, 30(10): 2923-2928.2 刘洁平, 张令心. 高层建筑土-结构相互作用地震反应分析简化评估方法J. 土木工程学报, 2010, 43(12): 28-34. Liu Jie-ping, Zhang Ling-xin. A simplified assessment method for seismic response analysis ofsoil-structure interaction of high-rise buildingsJ. China Civil Engineering Journal, 2010, 43(12): 28-34.3 陈少林, 唐敢, 刘启方, 等. 四维土结构动力相互作用的一种时域直接分析方法J. 地震工程与工程振动, 2010, 30(2): 24-31. Cheng Shao-lin, Tang Gan, Liu Qi-fang, et al. A direct time-domain method for analysis of three-dimensional soil-structure dynamic interactionJ. Journal of Earthquake Engineering and Engineering Vibration, 2010, 30(2): 24-31.4 张之颖, 赵钟斗, 吕西林, 等. SSI体系阻尼特性振动台模型试验研究J. 土木工程学报, 2010, 43(2): 100-104. Zhang Zhi-ying, Zhao Zhong-dou, Lu Xi-lin, et al. Shaking table tests of the damping behavior of SSI systemsJ. China Civil Engineering Journal, 2010, 43(2): 100-104.5 刘伟庆, 李昌平, 王曙光,等. 不同土性地基上高层隔震结构振动台试验对比研究J. 振动与冲击, 2013, 32(16): 128-133. Liu Wei-qing, Li Chang-ping, Wang Shu-guang, et al. Comparative study on high-rise isolated structure founded on various soil foundations by using shaking table testsJ. Journal of Vibration and Shock, 2013, 32(16): 128-133.6 易伟建, 夏玲琼. 近场地震下考虑桩-土-框架结构相互作用的拟动力试验J. 土木工程学报, 2014, 47(4): 1-8. Yi Wei-jian, Xia Ling-qiong. Pseudo dynamic experimental analysis of pile-soil-frame interaction in near-fault earthquakesJ. China Civil Engineering Journal, 2014, 47(4): 1-8.7 刘毅, 薛素铎, 李雄彦. 土结构相互作用下网架结构动力性能研究J. 振动与冲击, 2014, 33(10): 22-28. Liu Yi, Xue Su-duo, Li Xiong-yan. Grid structure dynamic performance analysis considering soil-structure interactionJ. Journal of Vibration and Shock, 2014, 33(10): 22-28.8 王博, 徐建国. 大型渡槽考虑土-结构相互作用的动力分析J. 世界地震工程, 2000, 16(3): 110-115. Wang Bo, Xu Jian-guo. Dynamic analysis of large aqueducts with soil-structure interactionJ. World Earthquake Engineering, 2000, 16(3): 110-115.9 Larkin T. Seismic response of liquid storage tanks incorporating soil structure interactionJ. Journal of Geotechnical and Geoenvironmental Engineering, 2008, 134: 1804-1814.10 程选生, 杜永峰. 弹性地基上矩形贮液结构的液-固耦合振动特性J. 工程力学, 2011, 28(2): 186-191. Cheng Xuan-sheng, Du Yong-feng. Vibration characteristic analysis of rectangular liquid-storage structures considering liquid-solid coupling on elastic foundationJ. Engineering Mechanics, 2011, 28(2): 186-191.11 Kianoush M R, Ghaemmaghami A R. The effect of earthquake frequency content on the seismic behavior of concrete rectangular liquid tanks using the finite element method incorporating soil-structure interactionJ. Engineering Structures, 2011, 33: 2186-2200.12 Patel C N, Patel H S. Soil-Foundation-Structure Interaction Effects in Seismic behaviour of RC Elevated Water Tank C / Proceedings of International Conference on Advances in Tribology and Engineering Systems, Part of the series Lecture Notes in Mechanical Engineering. 2013: 465-477.13 张华, 纪爱丽, 徐威, 等. 考虑土结构相互作用渡槽流固耦合体风振响应分析J. 中南大学学报(自然科学版), 2015, 46(5): 1845-1850. Zhang Hua, Ji Ai-li, Xu Wei, et al. Analysis of wind-induced vibration response of fluid-structure interaction system for aqueduct considering soil-structure interactionJ. Journal of Central South University (Science and Technology), 2015, 46(5): 1845-1850.14 程选生, 杜永峰, 李慧. 钢筋混凝土矩形贮液结构的液-固耦合晃动弹性底板J. 应用力学学报, 2008,25(4): 622-626. Cheng Xuan-sheng, Du Yong-feng, Li Hui. Liquid-solid coupling sloshing of reinforced concrete rectangular liquid-storage tankselastic soleplateJ. Chinese Journal of Applied Mechanics, 2008,25(4): 622-626.15 程选生. 钢筋混凝土矩形贮液结构的液-固耦合振动J. 煤炭学报, 2009, 34(3): 340-344. Cheng Xuan-sheng. Liquid-solid coupling vibration of reinforced concrete rectangular liquid-storage tanks J. Journal of China