考虑土-储罐-液体相互作用的隔震储罐简化力学模型及地震响应研究.pdf

第 37 卷第 3 期2024 年 3 月振 动 工 程 学 报Journal of Vibration EngineeringVol.37 No.3Mar.2024考虑土储罐液体相互作用的隔震储罐简化力学模型及地震响应研究李想1，2，张书进1，3，孙建刚2，4，王尊策1，徐蕾2，崔利富2，吕远5（1.东北石油大学机械科学与工程学院，黑龙江 大庆 163318；2.大连民族大学土木工程学院，辽宁 大连 116650；3.中国石油大庆油田有限责任公司采油工程研究院，黑龙江 大庆 163453；4.防灾科技学院土木工程学院，北京 101601；5.哈尔滨工业大学土木与环境工程学院，广东 深圳 518055）摘要:为了寻求一种既能减震又能降低成本的储罐结构，提出一种隔震储罐结构体系，根据力的平衡原则推导了滚动隔震装置恢复力模型，得出了复合滚动隔震装置的恢复力模型。基于三质点模型和场地土模型，提出了考虑土储罐液体相互作用（STLI）的隔震储罐的简化力学模型和运动方程，并研究了考虑 STLI效应和不考虑 STLI效应时抗震储罐和隔震储罐在不同场地的地震响应。结果表明：隔震储罐能够有效地降低其基底剪力及倾覆弯矩，但对晃动波高的控制有限，建议在高烈度区，满足晃动波高的前提下，隔震储罐在设计时可以降低烈度。考虑 STLI效应后，抗震储罐的基底剪力和倾覆弯矩明显降低，且从类场地到类场地差异率逐渐增大，软土场地降低最为显著。隔震储罐地震响应受 STLI效应的影响较小，可以有效地隔断上部结构与场地土之间的耦联，弱化 STLI效应对上部结构的影响。关键词:土储罐液体相互作用；隔震储罐；复合滚动隔震；简化力学模型；地震响应中图分类号:TU352.1 文献标志码:A 文章编号:1004-4523（2024）03-0436-12 DOI：10.16385/ki.issn.1004-4523.2024.03.008引言立式石油储罐是石油化工行业重要的存储设施，由于其存储易燃、易爆和有毒介质，一旦发生地震灾害容易引起火灾、爆炸或者造成环境污染等次生灾害，对生态环境和人类生存造成严重的危害。1964年 6月 16日，日本新潟地震中储油罐破坏后发生大火和爆炸，大火持续 15天，烧毁油罐 84座，造成严重的经济损失和空气污染1。1978 年 6 月 12 日，日本宫城地震，成品油罐被破坏，造成了陆地和海洋的大面积污染。1976 年 7 月 28 日，唐山大地震，唐山钢厂和天津化工厂各两座 1000 m3的燃料油罐最下层圈壁发生外屈曲，角焊缝开裂，罐内贮油全部泄漏。立式储罐的震害主要表现为罐壁底部的象足式鼓曲，浮顶上部构件损坏，锚固件、罐底板和罐底贴角焊缝的破坏，管接头及其附件的破坏以及地基液化等。由以上储罐震害可知，立式储罐壁薄、容积大，其地震响应主要涉及提离、罐壁应力、罐壁失稳等。为此国内外学者针对储罐结构的特点，从理论、数值仿真以及试验多维度揭示复杂的地震响应机理，从而提高储罐抗震能力，避免其因为地震作用而产生破坏。就当前的研究成果来看，结构控制思想的引入很好地起到了减小地震破坏作用的效果29。Chalhoub10进行了采用基础隔震技术的储罐振动台试验，结果表明，采用橡胶隔震装置的储罐能有效地降 低 地 震 响 应，但 晃 动 高 度 略 有 增 加。Shrimali等1112研究了单向和双向地震激励下采用叠层橡胶隔震支座的储罐的地震响应，结果表明该方法可以有效地降低储罐动力响应。Cheng等13考虑储罐液固耦合的影响，研究了采用橡胶隔震支座的矩形钢筋混凝土液体储罐的动态响应。结果表明，由于橡胶隔震支座过滤了外部激励的高频分量，因此储罐不会发生高频共振，当外部激发频率接近液体的一阶频率时，会发生显著的共振。孙建刚14从罐壁的柔性和液固耦合运动出发，建立了隔震三质点控制体系简化分析的力学模型，给出了立式储罐基础隔震控制体系的运动方程。目 前，对 储 罐 隔 震 问 题 的 研 究 进 展 十 分 迅速1516，然而大多数研究假设储罐基础是刚性的，忽略了土储罐液体相互作用对地震响应的影响。近几十年来，国内外学者广泛研究了土结构相互作用收稿日期:2022-07-12；修订日期:2022-09-13基金项目:国家自然科学基金资助项目（52178461，51878124，51478090）；广东省重点领域研发计划项目（2019B111102001）。第 3 期李想，等：考虑土储罐液体相互作用的隔震储罐简化力学模型及地震响应研究对各种结构体系地震响应的影响1722。由于液体具有流动性，所以 SSI 效应对储罐的影响正逐渐转化为 更 为 复 杂 的 土 储 罐 液 体 相 互 作 用（STLI）。Haroun 等23研究了考虑土与储罐相互作用的立式储罐地震响应，研究表明在水平地震作用下，土与结构相互作用能放大储罐的动水压力和倾覆弯矩。孙建刚等24采用数值方法研究了考虑土储罐相互作用的大型立式浮顶储罐的地震响应，研究表明，土壤摆动效应对结构的基底剪力和倾覆弯矩有很大的影响，特别是当场地条件为软土时。综上所述，隔震对减小储罐地震响应具有很好的效果，而且 STLI 效应对储罐结构地震响应的影响是不可忽视的。从目前国内外关于储罐隔震体系的研究来看，要实现上述隔震体系，就要使用大量的隔震装置，储罐基础构造也要产生巨大的改变，这必将增加结构体系的建造成本，也使得隔震技术难以大面积推广应用在储罐中。鉴于此，本文在局部改变原结构体系的基础上，提出可以有效降低结构地震响应的隔震储罐结构体系。目前关于土与储罐相互作用的研究仅针对立式储罐，而隔震储罐与立式储罐在结构形式上有很大不同，因此，本文推导了适用于该结构体系的复合滚动隔震装置的恢复力模型，而后提出了基于 STLI 效应的隔震储罐简化力学模型，研究了考虑 STLI效应和忽略 STLI效应的隔震储罐和抗震储罐在不同场地的地震响应，探讨了考虑 STLI效应的必要性。1隔震储罐当前的储罐基础隔震体系如图 1 所示，采用隔震装置将储罐与其基础分隔开。本文所提出的隔震储罐（如图 2 所示）将罐壁和储罐底板分开，并将储罐底板支撑在砂垫层上，罐壁支撑在复合滚动隔震装置上，采用柔性膜将罐壁和储罐底板连接起来，本文所采用的柔性膜是一种新型弹性体材料，该材料具有韧性高、防水、耐磨、强度高的特性，并具有良好的粘结性，可以保证罐壁在水平方向运动，同时防止储罐内液体泄漏。该隔震储罐由于罐壁和储罐底板是断开的，并且通过柔性膜连接，所以当发生地震时砂垫层的弓起作用和底板提离作用的作用力传递路径被切断，可以消除弓起破坏和底板提离。由于复合滚动隔震装置的引入能够很好地降低储罐地震响应，而且隔震装置仅作用在罐壁底部，所以能大大减少隔震装置的布置数量，储罐基础结构形式的变化也不大，能有效地降低建造成本，具有良好的经济效益。2复合滚动隔震装置2.1复合滚动隔震装置的构成复合滚动隔震装置由滚动隔震装置和铅芯阻尼器并联而成，如图 3所示。滚动隔震装置由带有上、下凹槽的盖板和滚球构成，凹槽设计为圆弧形。将铅芯阻尼器插入橡胶管中，而后将其插入盖板，其中上盖板为贯穿孔，下盖板未贯穿。盖板的开孔需比橡胶管尺寸略小，利用橡胶的挤压受力使铅芯阻尼器固定，其中铅芯阻尼器高度略小于隔震装置的整体高度，以消解滚珠滚动时所引起的复合滚动隔震装置竖向位移，并在铅芯阻尼器顶部布置一根弹簧，用螺栓固定。图 1 储罐基础隔震体系Fig.1 Foundation isolation system of tank图 2 本文提出的隔震储罐结构体系Fig.2 Structure system of proposed isolation tank图 3 复合滚动隔震装置Fig.3 Composite rolling isolation device437振 动 工 程 学 报第 37 卷2.2滚动隔震装置恢复力模型参考文献 25，通过对滚动隔震装置的恢复力模型进行计算可知，其自振周期和刚度与单摆类似，具体表达式为：Ti=2Sg=22(R-r)g（1）k=2m（2）式中S为上凹槽质心运动轨迹；R 为凹槽半径；r为滚球半径；为滚动隔震装置圆频率；m 为滚动隔震装置上部质量。从式（1）和（2）中可以看出，滚动隔震装置的周期和刚度仅与凹槽和滚球的半径有关，而与其相对位移无关。这种方法大大简化了恢复力模型的计算过程，当滚动隔震装置发生较大相对位移时，采用上述方法计算会产生较大误差。本文根据力的平衡原则推导出滚动隔震装置变周期和变刚度恢复力模型。以滚动隔震装置中一个滚动单元为研究对象，其上凹槽受力分析如图 4所示。根据力的平衡原则可得滚球与上凹槽接触面的平衡方程：W cos +F sin -N=0（3）W sin -F cos +T=0（4）式中W 为上部结构作用于滚动隔震装置的竖向荷载；F 为滚动隔震装置的恢复力；T 和 N 分别为上凹槽和滚球接触面的切向摩擦力和法向反力；为转动角度。切向摩擦力 T可表示为26：T=sgn(v)Nr=sgn(xcos)Wr cos 11-sgn()xcos rtan （5）根据式（3）（5）可得：F=W tan +sgn(x)Wr1+tan21-sgn()xrtan （6）假设滚球中心的运动方程为：x2+y2=(R-r)2（7）根据式（6）可求得滚动隔震装置恢复力模型为：F=-W4(R-r)2-x2i-12xi+sgn(xi)Wr1+()y21-sgn()xiry=-kxi+Fd（8）式中xi为上凹槽的位移；x 为滚球中心的位移；x为滚球中心的速度。式（8）中第一项为滚动隔震装置弹性恢复力，第二项为滚动隔震装置阻尼力Fd。滚动隔震装置刚度为：k=W4(R-r)2-x2i-12（9）滚动隔震装置周期为：Ti=2mk=24(R-r)2-x2i12g（10）由式（10）可知，滚动隔震装置的隔震周期主要与凹槽半径 R 和滚球半径 r 及其所处的位置有关。在进行隔震设计时，应首先确定滚球的尺寸大小，而后根据所需隔震周期进行凹槽半径的设计。2.3滚动隔震装置恢复力模型验证为了验证滚动隔震装置恢复力模型的正确性，对比分析滚动隔震装置恢复力模型和有限元仿真模型在水平正弦位移激励下的计算结果。选取凹槽半径 R=0.6 m，滚球半径 r=0.075 m。将滚动隔震装置简化为二维平面内的运动，采用 ADINA 有限元软件建立仿真模型，如图 5 所示。单元材料模型选用双线性的弹塑性模型，材料弹性模量为 206 GPa，屈 服 强 度 为 518 MPa，泊 松 比 为 0.3，密 度 为7850 kg/m3。在模型顶部施加 1 N/m 的线荷载以及如图 6 所示的正弦水平位移激励y=A sin(6t)，滚动隔震装置恢复力模型计算结果与有限元仿真模型计算结果如图 7所示。从图 7 中可以看出，滚动隔震恢复力模型和有限元仿真模型的计算结果十分接近，验证了本文所图 4 上凹槽受力分析图Fig.4 Stress analysis of upper groove438第 3 期李想，等：考虑土储罐液体相互作用的隔震储罐简化力学模型及地震响应研究推导的滚动隔震装置恢复力学模型的准确性。2.4复合滚动隔震装置恢复力模型采用复合滚动隔震装置的隔震储罐在遭遇大风天气或者小幅地震时，滚动隔震装置的起滚力和铅芯阻尼器的初始刚度能保证储罐上部结构不发生大幅振动。当遭遇强震时，滚动隔震装置带动铅芯阻尼器进行水平运动，并在往复运动中耗散地震能量。本文铅芯阻尼器的恢复力模型采用 BoucWen光滑型恢复力模型27：Fl=klx+(1-)klz（11）z=Aix-i|x|zi-1z-x|zi（12）式中为屈服后与屈服前的水平刚度之比；kl为铅芯的弹性刚度；z为滞变位移。Ai，i，i等参数可通过参数识别得到。根据式（8），可得复合滚动隔震装置恢复力模型为：F=kfx+Ff+klx+(1-)klz（13）式中kf为复合滚动隔震装置的自复位刚度。3考 虑 STLI 效 应 的 隔 震 储 罐 简 化力学模型3.1场地土简化模型本文将场地土简化为三自由度力学模型2022，分别为平动自由度xH(t)，摆动自由度x(t)和附加自由度x(t)，如图 8所示。图 8 中，ri，e，mf和 If分别为结构基础的半径、埋置深度、质量和质量惯性矩；水平平动刚度系数kH和阻尼系数cH的计算公式分别为2122：kH=82sri2-()1+eri（14）cH=ris()0.68+0.57erikH（15）式中为土的密度；为土的泊松比；s为场地土的等效剪切波速。摆动的转动刚度系数k和阻尼系数c的计算公式分别为2122：k=82sri33(1-)1+2.3eri+0.58()eri3（16）c=ris 0.15631eri-0.08906()eri2-0.00874()eri3kH（17）土体附加的自由度包括质量惯性矩I和阻尼参数c2122，分别表示为：I=()ris2 0.33+0.1()eri2k（18）c=ris 0.4+0.03()eri2k（19）储罐结构的动响应包括液体晃动和液固耦合等复杂的动态过程，目前国内外常用 Housner28和图 5 有限元仿真模型Fig.5 Finite element simulation model图 6 水平位移激励Fig.6 Horizontal displacement excitation图 7 滚动隔震装置恢复力模型和有限元仿真模型滞回曲线Fig.7 Hysteresis curves of restoring force model and finite element simulation model of rolling isolation device图 8 土-结构力学模型Fig.8 Soil-structure mechanical model439振 动 工 程 学 报第 37 卷孙建刚等24提出的刚性理论和弹性理论力学模型来解决储罐复杂的动响应问题，本文基于上述理论推断，考虑土结构相互作用的隔震储罐简化力学模型实际上是土壤模型与势流体理论的有机组合。将图8中的场地土模型与隔震储罐结合，构成图 9中考虑STLI效应的隔震储罐的力学分析模型。3.2隔震储罐简化动力学模型基于弹性理论的立式储罐水平基础隔震力学模型28，将隔震储罐内液体简化为对流质量 mc、液固耦合质量 mi和刚性质量 m0。等效高度分别为 hc，hi和 h0，储液高度为 Hl。液固耦合质量和对流质量通过等效弹簧刚度 kc，ki及阻尼系数 cc，ci与罐壁相连，k0为隔震装置等效刚度，Ff为由摩擦力构成的恢复力，Fl为铅芯阻尼器恢复力。地面运动位移、对流晃动位移、液固耦合位移和隔震装置位移分别为xg(t)，xc(t)，xi(t)和x0(t)。考虑场地土储罐液体相互作用（STLI）的影响，得出考虑 STLI 效应的隔震储罐简化力学模型，如图 10所示。3.3运动方程根据考虑 STLI效应的隔震储罐简化动力学模型，利用 Hamilton 原理得出六自由度体系的运动控制方程：MX+CX+KX=-F（20）其中：M=mc0mcmcmchc00mimimimihi0mcmimi+mc+m0mc+mi+m0mchc+mihi+m0h00mcmimc+mi+m0mc+mi+m0+mfmchc+mihi+m0h00mcHcmiHimcHc+miHi+m0H0mcHc+miHi+m0H0mch2c+mih2i+m0h20+I0000000I，C=cccc0cHc+c-c-cc，K=kckik0kHk0，X=xcxix0 xHxx，F=mcmimc+mi+m0mc+mi+m0+mfmchc+mihi+m0h00 xg+00Ff+Fl000，式中I0为罐体绕中心的转动惯量；Hc，Hi，H0表示3个质点的高度。隔震储罐的基底剪力、倾覆弯矩和晃动波高分别为：Q(t)=-mcxc(t)+x0(t)+xH(t)+hcx(t)+xg(t)-mixi(t)+x0(t)+xH(t)+hix(t)+xg(t)-m0 x0(t)+xH(t)+h0 x(t)+xg(t)（21）M(t)=-mchcxc(t)+x0(t)+xH(t)+hcx(t)+xg(t)-mihixi(t)+x0(t)+xH(t)+hix(t)+xg(t)-m0h0 x0(t)+xH(t)+h0 x(t)+xg(t)（22）hv=-0.837g x0(t)+xH(t)+Hlx(t)+xg(t)+xc(t)Ri（23）式中Ri为储罐半径。图 9 考虑 STLI效应的隔震储罐力学分析模型Fig.9 Mechanical analysis model of isolation tank considering STLI effect图 10 考虑 STLI效应的隔震储罐简化力学模型Fig.10 Simplified mechanical model of isolation tank considering STLI effect440第 3 期李想，等：考虑土储罐液体相互作用的隔震储罐简化力学模型及地震响应研究4简化力学模型验证4.1储罐基本参数及隔震装置设计本文以 3000 m3立式储罐为原型罐设计隔震储罐，其基本参数如表 1所示。设计隔震装置初始隔震周期为 Ti=3 s，滚动摩阻系数取为 0.005 cm，滚球的半径为 r=0.08 m，根据式（10）可算得凹槽半径 R=1.2 m。在隔震储罐罐壁底部均匀布置 6个复合滚动隔震装置。根据预期 隔 震 层 偏 移，选 定 铅 芯 阻 尼 器 有 效 高 度 H=300 mm，直径 D=100 mm，建立铅芯阻尼器有限元仿真模型，并对其进行拟静力数值仿真分析，其中模型选用 3D Solid单元，材料弹性模量为 16.5 GPa，屈服强度为 14 MPa，泊松比为 0.42。有限元模型如图11所示。在铅芯阻尼器顶部施加如图 12 所示的正弦水平位移激励y=A sin(2t)，铅芯阻尼器的滞回曲线如图 13所示，通过对数值仿真模型的滞回曲线进行参数识别，获取 BoucWen 光滑型恢复力模型的参数，如表 2所示。4.2有限元仿真模型根据隔震储罐所处位置的实际地质条件，建立考虑场地土影响的有限元仿真模型。其中对于土壤模型，当其直径 Di与高度 Hi之比大于 10 时29，可基本消除土壤边界对结构的影响。隔震储罐所处位置的土层总高度 Hi为 17 m，本文选择土壤直径 Di为210 m，Di/Hi10，满足要求。表 3 为隔震储罐所处位置的土层物理参数。采用 Shell单元模拟储罐罐壁和柔性膜，材料模型分别选用弹塑性模型和超弹性模型。采用 3D Fluid 单元模拟储罐内液体，液面设置为自由面单元。场地土采用 3D Solid单元，材料选用莫尔库仑材料模型。罐壁底部均匀布置 6组复合滚动隔震装置，采用带有高度的线弹簧单元模拟滚动隔震装置，通过对滚动隔震装置恢复力滞回曲线进行参数识别及数据拟合，即可获得表征滚动隔震装置水平方向力学性能的非线性弹簧单元物理参数，铅芯阻尼器采用 3D Solid 单元模拟，材料模型选用双线性的弹塑性模型。有限元模型如图 14所示。表 1 储罐基本参数Tab.1 Basic parameters of tank储罐半径/m8.25储液高度/m13.07壁厚/m6.1110-3罐壁弹性模量/Pa2.061011液体密度/（kg m-3）1103泊松比0.3高径比1.58基础半径/m8.5基础埋置深度/m1.2图 11 铅芯阻尼器有限元仿真模型Fig.11 Finite element simulation model of lead core damper图 12 水平位移激励Fig.12 Horizontal displacement excitation图 13 滞回曲线Fig.13 Hysteresis curves表 2 BoucWen模型参数Tab.2 Parameters of BoucWen modelAi1.5i17001770i1kl/(N m1)817327900.003441振 动 工 程 学 报第 37 卷4.3场地土分类及地震波选择根据等效剪切波速度和土层覆盖厚度可将场地土分为类、类、类和类。根据规范30可计算得出隔震储罐所处位置土层的等效剪切波速和覆盖层厚度分别为 146.3 m/s 和 17 m，属于类场地。因此，所选的地震波必须符合类场地的特点。本文选择 3 条实际记录的地震波和 1 条人工合成的地震波作为地震动输入，加速度峰值调整为 0.2g，加速度反应谱曲线如图 15所示。4.4计算结果对比分析通过计算得出采用有限元模型和理论模型的考虑 STLI效应的隔震储罐地震响应，表 4为理论模型和有限元模型计算出的地震响应峰值，图 16为 Taft波作用下考虑 STLI 效应的隔震储罐基底剪力、倾覆弯矩和晃动波高的时程曲线。从表 4 和图 16 中的数据可以看出，有限元模型和理论模型得出的地震响应峰值非常接近，而且地震响应时程曲线契合较好。同时从表 4中可以看出有限元模型计算所得地震响应峰值均大于简化力学模型，且平均差异率均在 20%以内，有限元数值仿真结果与简化力学模型计算结果相差不大，两者互为验证，所以本文提出的考虑 STLI 的隔震储罐简化力学模型是可靠的。5算例分析参考 4.1 节储罐基本参数及隔震装置设计进行隔震储罐地震响应研究，储罐基本参数如表 1所示，其中 e/ri=0.14。5.1地震动输入根据规范30要求，本文选择，类场地中 3 条实际记录的地震波和 1 条人工合成的地震波图 15 加速度反应谱Fig.15 Acceleration response spectrum表 3 土层物理参数Tab.3 Physical parameters of soil layer土层素填土淤泥质黏土粉质黏土粗砂覆盖厚度/m2.665.43密度/(kg m3)1600159018401850弹性模量/MPa48.4428.04296.52577.42泊松比0.40.450.40.35内摩擦角/（）12101628黏聚力/kPa88200剪切波速/（m s1）125120156265图 14考虑 STLI效应的隔震储罐有限元模型Fig.14Finite element model of isolation tank considering STLI effect442第 3 期李想，等：考虑土储罐液体相互作用的隔震储罐简化力学模型及地震响应研究作为地震动输入。研究考虑 STLI效应的隔震储罐地震响应，调整加速度时程曲线峰值为 0.4g，四类场地加速度反应谱如图 17 所示。虽然场地土的类型与土壤的软硬程度和土壤厚度相关30，但是场地土简化力学模型只涉及等效剪切波速2122，所以本文以不同剪切波速来表示不同场地土类型。5.2地震响应研究本文采用Newmark法，对比分析了不同场地忽略 STLI 效应和考虑 STLI 效应的抗震储罐和隔震储罐地震响应，具体计算结果如表 5，6所示。从表 5，6数据中可以看出，不同场地地震波作用下，不考虑 STLI效应和考虑 STLI效应的隔震储罐都可以有效降低储罐的基底剪力和倾覆弯矩，但不能有效控制晃动波高，甚至有放大波高的效应。隔震储罐基底剪力和倾覆弯矩与抗震储罐相比，在坚硬场地土（类场地、类场地、类场地）上的减震效果明显好于在软弱场地土（类场地）上的减震效果，其减震率达 67%以上；对于类场地（软土），隔震储罐的基底剪力和倾覆力矩也有明显的减震效果，其减震率也可以达到30%以上。对于晃动波高的减震效果而言，不同场地地震波作用下，隔震储罐的减震效果差异较大，有时隔震对晃动波高有放大效应。在人工波 4（类场地）作用下，隔震储罐较抗震储罐的晃动响应明显增加，隔震后晃动波高较隔震前的晃动波高增加了 60%以上，这主要是因为隔震装置会延长结构的自振周期，而晃动波高属于长周期振动，因此隔震储罐对晃动表 4 理论模型和有限元模型地震响应峰值对比Tab.4 Seismic response comparison of theoretical model and finite element model基底剪力/(106 N)倾覆弯矩/(107 N m)晃动波高/m理论模型有限元模型理论模型有限元模型理论模型有限元模型Taft1.2481.3410.9331.0560.3270.348PEL1.1961.3780.9711.1860.4450.518LWD1.2691.130.9411.190.1370.155人工波 31.5261.651.1771.3470.5440.679平均值1.3101.3751.0061.1950.3630.425平均值差异率/%4.9618.8217注：平均值差异率=|理论模型-有限模型理论模型100%。图 16 地震响应时程曲线Fig.16 Seismic response timehistory curves图 17 加速度反应谱Fig.17 Acceleration response spectrum443振 动 工 程 学 报第 37 卷波高的控制并不明显。而在其他地震波作用时，储罐的晃动波高与抗震储罐的晃动波高相比变化不大。由表 5，6 数据可知，考虑 STLI 效应后，抗震储罐的基底剪力和倾覆弯矩明显降低，且从类场地到类场地差异率逐渐增大，坚硬场地差异率最低，其中类场地土 Qianan 波作用下，考虑 STLI 效应的抗震储罐基底剪力和倾覆弯矩差异率分别为7.77%和 4.99%；软土场地差异率最为显著，其中类场地土 Pasadena 波作用下，考虑 STLI 效应的抗震储罐基底剪力和倾覆弯矩差异率分别为 40.52%和 41.93%；而考虑 STLI效应后抗震储罐的晃动波高变化不明显。这主要是因为考虑 STLI 效应后，相当于在储罐底部植入了柔性层，等效于在抗震储罐底部增加了隔震装置，所以考虑 STLI 效应后基底剪力和倾覆弯矩逐渐降低，由于液体晃动属于长周期运动，所以考虑 STLI 效应对晃动波高的影响不明显，而由于类场地接近刚性，土体变形对结构的反应影响很小，所以类场地土地震波作用下，考虑和忽略 STLI效应对抗震储罐地震响应差异率明显低于其他类型场地土。隔震储罐地震响应受 STLI 效应的影响较小，两者基底剪力、倾覆弯矩和晃动波高的差异率最大值为 13.53%。这是由于场地土的等效侧向刚度远大于隔震储罐隔震层的水平侧向刚度，因此场地土对隔震储罐的自振周期影响很小。不考虑 STLI 效应和考虑 STLI 效应的隔震储罐其自振周期分别为 1.962 和 1.957 s，基本保持一致。说明采用复合滚动隔震装置的隔震储罐可以有效地隔断上部结构与场地土之间的耦联，弱化 STLI效应对上部结构的影响。总之，隔震储罐比传统抗震储罐的基底剪力和倾覆力矩明显降低，由结构抗震知识可知，减震率在75%时，相当于抗震烈度降低 2 度，减震率在 50%时，相当于抗震烈度降低 1度。综上分析，建议在高表 5，类场地地震响应对比Tab.5 Seismic response comparison of sites，计算内容忽略STLI效应考虑STLI效应考虑和忽略STLI效应影响的对比分析基底剪力/(106 N)倾覆弯矩/(107 N m)晃动波高/m基底剪力/(106 N)倾覆弯矩/(107 N m)晃动波高/m基底剪力差异率/%倾覆弯矩差异率/%晃动波高差异率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震抗震隔震抗震隔震类场地Qianan18.6911.40892.4711.8651.01891.420.3080.2944.6417.2381.40791.8311.2731.01291.020.3080.2925.087.770.044.990.6200.58Guangzhou23.0771.11295.1814.3231.28191.060.1560.1550.920.731.12694.5713.1861.36689.640.1560.163-4.2710.171.277.946.70.065.17TH1TG02520.8841.49192.8613.381.30890.220.4440.4312.9318.7111.5291.8811.851.28989.120.4450.4254.4610.411.9511.431.450.231.39人工波 123.3081.74592.5114.3761.24491.350.2060.224-8.7421.2661.71391.9413.0281.19790.810.2070.224-8.658.761.839.383.780.490类场地Lanzou20.4171.55692.3813.2661.07191.930.4040.32719.0915.5712.71582.569.6842.04378.911.2851.42-10.4425.440.1826.743.850.030.61RH220.8842.7286.9813.2191.96785.121.2851.411-9.8113.6351.54188.78.516188.260.4050.33118.2133.220.9735.816.60.11.19TH1TG03519.1971.4492.511.9681.19590.020.5010.4970.8413.7851.47589.38.521.20185.90.5020.4921.9428.192.4528.810.530.11.03人工波 221.4233.03585.8313.9012.04485.31.151.172-1.9117.0483.06582.0210.3122.04680.161.1511.176-2.2220.420.9825.820.110.060.37444第 3 期李想，等：考虑土储罐液体相互作用的隔震储罐简化力学模型及地震响应研究烈度地区，在保证液体晃动波高的要求下，可以采用降烈度设计隔震储罐。6结论（1）根据力的平衡原则推导了滚动隔震装置恢复力模型，并对比分析了滚动隔震恢复力模型和有限元仿真模型在水平正弦位移激励下的计算结果，两者滞回曲线的计算结果十分接近，验证了本文所推导的滚动隔震恢复力学模型的准确性。提出并设计了复合滚动隔震装置，建立了相应的恢复力力学模型。（2）提出了隔震储罐结构体系，依据速度势理论，考虑土储罐液体的相互作用（STLI）的影响，构建了其简化力学模型及相应的运动控制方程。（3）隔震储罐能够有效地降低储罐的基底剪力及倾覆弯矩，但对储液晃动波高的控制有限。考虑STLI效应后，抗震储罐的基底剪力和倾覆弯矩明显降低，且从类场地到类场地差异率逐渐增大，软土场地降低最为显著。隔震储罐地震响应受 STLI效应的影响较小，可以有效地隔断上部结构与场地土之间的耦联，弱化 STLI效应对上部结构的影响。参考文献:1孙建刚.立式储罐地震响应控制研究 D.哈尔滨：中国地震局工程力学研究所，2002.Sun Jiangang.Research on the earthquake response control of the vertical storage tanksD.Harbin：Institute of Engineering Mechanics，CEA，2002.2 Drosos J C，Tsinopoulos S V，Karabalis D L.Seismic retrofit of spherical liquid storage tanks with energy dissipative devices J.Soil Dynamics and Earthquake Engineering，2019，119：158-169.3Miguel Ormeo，Tam Larkin，Nawawi Chouw.Experimental study of the effect of a flexible base on the seismic response of a liquid storage tankJ.Thin-Walled Structures，2019，139：334-346.表 6，类场地地震响应对比Tab.6 Seismic response comparison of sites，计算内容忽略STLI效应考虑STLI效应考虑和忽略STLI效应影响的对比分析基底剪力/(106 N)倾覆弯矩/(107 N m)晃动波高/m基底剪力/(106 N)倾覆弯矩/(107 N m)晃动波高/m基底剪力差异率/%倾覆弯矩差异率/%晃动波高差异率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震减震率/%抗震隔震抗震隔震抗震隔震类场地Taft20.1202.42687.9412.2031.51587.590.6470.714-10.3612.512.50080.057.5651.35782.060.6520.72-10.7637.823.0538.0110.430.821.19PEL24.2172.37690.1915.1911.56989.670.9600.9451.5613.552.42482.118.13461.66679.530.9640.951.5244.042.0246.456.150.420.42LWD29.9651.724394.24518.631.33792.8250.19410.161216.96516.6871.72589.66310.221.15688.6910.1980.16815.244.310.0345.1413.531.753.91人工波 316.7933.07681.6810.7572.19579.601.1171.194-6.9110.7143.0771.356.5842.17167.031.1191.201-7.3236.200.2138.791.090.230.61类场地Pasadena16.2004.62471.469.5892.41474.830.9520.9094.529.6364.7750.465.5692.50155.100.9560.9213.6740.523.2441.933.580.461.35TH415.6752.13686.379.8731.17188.140.3450.382-10.729.1132.2475.455.1021.17676.950.3510.384-9.3541.864.7248.320.441.740.47TH4TG09022.7263.17886.0214.4802.24584.502.4491.621733.7710.6143.2069.816.3712.18265.752.4441.61733.8453.300.8356.002.800.180.28人工波 414.7567.30650.498.8463.86656.300.6471.041-60.9010.6697.4230.416.0883.96834.820.6661.065-60.0427.701.6231.182.642.892.33445振 动 工 程 学 报第 37 卷4Zhang R F，Weng D G，Ren X S.Seismic analysis of a LNG storage tank isolated by a multiple friction pendulum systemJ.Earthquake Engineering and Engineering Vibration，2011，10（2）：253-262.5Colombo J I，Almazn J L.Experimental investigation on the seismic isolation for a legged wine storage tankJ.Jo