CN102706525A - 一种lng 储罐隔震效应的确定方法 - Google Patents

Abstract

本发明提供一种LNG储罐隔震效应的确定方法，其包括以下步骤：确定LNG储罐的基本信息，计算内罐模型参数，计算外罐模型参数，计算基础参数，建立LNG储罐隔震力学模型，根据储罐所在地的地震反应谱选择地震波，对建立的分析模型进行特定地震波作用下的动力反应分析，从而确定LNG储罐隔震效应。本发明把隔震技术应用到LNG储罐这一固液共同作用的特殊工业结构上，具体地，把储罐的力学分析模型建立到成熟的具有隔震支座恢复力模型的商业软件上，即可方便地进行时程分析，求解出储罐结构的动力反应，从而很好地确定LNG储罐的隔震效应。

Description

一种LNG 储罐隔震效应的确定方法

技术领域

本发明涉及一种LNG储罐的性能分析方法，更具体地涉及一种LNG储罐隔震效应的确定方法。

背景技术

我国是世界上地震灾害最严重的国家之一，比如近年来发生在我国的汶川地震和舟曲地震自然灾害，因此新建和既有建筑在地震发生时是否能够不倒塌是结构设计人员的重要任务，特别是LNG储罐，LNG储罐一旦遭到地震破坏，不但威胁到天然气的存储和影响人们的正常生产生活，而且引发的二次灾害对周围居民生活和环境的危害更是十分巨大。

根据LNG行业设计规范NFPA 59A-2006或EN 14620-1:2006的要求，LNG储罐应进行在操作基准地震(OBE)和安全停车地震(SSE)作用下的承载力极限状态和正常使用极限状态验算，以确保储罐的安全。其中，OBE的定义为50年内超越概率为10%(重现期为475年)，阻尼比为5%的反应谱表示的地震；SSE的定义为50年内超越概率为2%(重现期为2475年)或1%(重现期为4975年)，阻尼比为5%的反应谱表示的地震。

从上述地震计算的定义来看，OBE地震相当于中国规范的基本烈度地震，也即是多遇地震的3倍大小；SSE地震则相当于中国规范的罕遇地震。由此可见，储罐的地震作用设防标准要比普通建筑物的标准高很多，故地震作用力也是非常大的。因此，如何有效降低储罐在地震作用下的响应是LNG储罐设计中的一个重要环节。

基底隔震技术被证明是降低地震作用的有效途径。隔震原理是通过结构底部加入隔震装置（如橡胶支座、摩擦摆支座）延长储罐自振周期，避开场地的特征周期，减小地震力向上部结构的传递，达到降低结构地震响应的目的。现阶段的隔震技术发展得较为成熟了，各国的建筑结构规范也都对常用结构形式的隔震设计提出了具体的要求，但对如何把隔震技术应用到LNG储罐这一固液共同作用的特殊工业结构上，则还是空白。

目前，国内的LNG事业正处于高速发展的阶段，国内其他工程设计单位的设计建造技术还处于起步阶段，还没有独立完成设计建造LNG储罐的例子。

发明内容

本发明目的在于提供一种LNG储罐隔震效应的确定方法，其中把隔震技术应用到LNG储罐这一固液共同作用的特殊工业结构上，具体地，把储罐的力学分析模型建立到成熟的具有隔震支座恢复力模型的商业软件上，即可方便地进行时程分析，求解出储罐结构的动力反应，从而很好地确定LNG储罐的隔震效应。

本发明采用如下技术方案来实现，提供一种LNG储罐隔震效应的确定方法，其特征在于步骤如下：

(1)确定LNG储罐的基本信息：

a)确定钢质内罐的基本参数，利用第一长度测量装置确定内罐的直径D(m)，并确定内罐半径R(m)，内罐罐壁等效厚度t(mm)，内罐材料弹性模量E(MPa)，液体的密度ρ_L(kg/m³)和液体的设计液位H(m)这些参数；

b)确定混凝土外罐的基本参数，利用第二长度测量装置确定罐壁外表面直径D₁(m)、罐壁内表面直径d(m)以及罐壁高度h_w(m)，并确定罐壁质量m_w(kg)，罐顶高度h_r(m)，罐顶质量m_r(kg)，底板厚度h_b(m)，底板质量m_b(kg)和混凝土弹性模量E_c(MPa)；

c)确定单桩的水平刚度k_h(N/m)和桩身竖向刚度k_{V_P}(N/m)；

d)确定底板的桩位布置图及桩根数n_P；

e)确定隔震支座的力学参数，弹性刚度K_u(N/m)、屈服后刚度K_d(N/m)、屈服力F_y(N)和竖向刚度k_{V_B}(N/m)；

f)确定储罐所在场地的地震影响系数曲线，该地震影响系数曲线即一反应谱；

(2)计算内罐模型参数

$m_L=\frac{\mathrm{\π}{D}^2}{4}{\mathrm{\ρ}}_{L}H,$ m_L-为液体总质量(kg)；

$m_c=0.230\frac{D}{H}\tanh \left(\frac{3.67H}{D}\right)m_L,$ m_c—为内罐液体晃动对流部分质量(kg)；

$h_c=[1.0-\frac{\cosh \left(\frac{3.67H}{D}\right)-1.937}{\frac{3.67H}{D}\sinh \left(\frac{3.67H}{D}\right)}]H,$ h_c—为晃动对流质量作用点等效高度；

$T_c=\frac{1.0404\sqrt{D}}{\sqrt{\tanh \left(\frac{3.67H}{D}\right)}},$ T_c—液体的晃动对流周期(s)；

$K_c={\left(\frac{2\mathrm{\π}}{T_c}\right)}^2{m}_c,$ K_c—为晃动对流作用等效刚度(N/m)；

$C_c=2{\mathrm{\ξ}}_c\sqrt{K_c.m_c};$ c_c—为晃动对流液体的阻尼(N.s/m)，ξ_c为液体晃动阻尼比；

$m_i=\left\{\begin{array}{cc}\frac{\tanh \left(0.866\frac{D}{H}\right)m_L}{0.866\frac{D}{H}}& D/H\&GreaterEqual;1.333\\ (1.0-0.218\frac{D}{H})m_L& D/H<1.333\end{array}\right.;$ m_i—为内罐液体冲击部分质量(kg)；

$h_i=0.375[1.0+1.333[\frac{0.866\frac{D}{H}}{\tanh \left(0.866\frac{D}{H}\right)}-1.0\left]\right]H$ 当D/H≥1.333

$h_i=[0.500+0.060\frac{D}{H}]H$ 当D/H＜1.333，h_i—为冲击质量作用点等效高度(考虑了液体对底板的动压力)(m)；

$T_i=C_i\frac{h\sqrt{\mathrm{\&rho;}}}{\sqrt{t/R}\&CenterDot;\sqrt{E}},$ $C_i=\left(\frac{2}{\sqrt{H/R}\&times;(0.46-0.15H/R+0.017{(H+/R)}^2)}\right),$ T_i—液体的冲击周期(s)；

$K_i={\left(\frac{2\mathrm{\&pi;}}{T_i}\right)}^2{m}_i,$ K_i—为冲击作用等效刚度(N/m)；

$C_i=2{\mathrm{\&xi;}}_i\sqrt{K_i.m_i};$ c_i—为冲击液体的阻尼(N.s/m)，ξ_i为液体冲击阻尼比；

(3)计算外罐模型参数

m_t＝m_w+m_r，m_t—罐壁和罐顶的总质量(kg)；

$h_t=\frac{0.5m_w{h}_w+m_r(h_w+\frac{1}{3}{h}_r)}{m_w+m_r},$ h_t—罐壁和罐顶的等效作用高度(m)；

$I=\frac{\mathrm{\&pi;}({D_1}^4-d^4)}{64},$ I—罐壁惯性矩(m⁴)；

$K_t=\frac{3E_{c}I}{h_w^3},$ K_t—为外罐的等效刚度(N/m)；

$C_t=2{\mathrm{\&xi;}}_t\sqrt{K_t.m_t},$ C_t—为外罐的阻尼(N.s/m)，ξ_t为混凝土的阻尼比；

(4)计算基础参数

$I_p=\underset{i=1}{\overset{n_p}{\mathrm{\&Sigma;}}}{x}_i^2,$ I_P—为桩群对底板中心的惯性矩(m²)，x_i—为桩心到底板中心的距离；

K_rot＝K_{V_P}I_p，K_rot—为桩群的转动刚度(N.m/rad)；

K_bv=min(K_{V_B}，K_{V_P})，K_bv—为基础竖向刚度和隔震支座的竖向刚度二者间的较小值；

c_bv=0，c_bv—为桩的竖向阻尼；

(5)建立LNG储罐隔震力学模型

利用上面的模型参数以及基础参数建立力学分析模型，并根据所述力学分析模型在结构软件上建立分析模型：

(6)根据储罐所在地的地震反应谱选择地震波

依据建筑抗震设计规范的规定，正确选择输入的地震加速度时程曲线，该时程曲线要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。

(7)对建立的分析模型进行特定地震波作用下的动力反应分析，从而确定LNG储罐隔震效应。

附图的详细描述

图1是LNG储罐的桩位布置图；

图2示出了LNG储罐隔震力学模型；

图3示出了隔震计算的有限元模型；

图4示出了EI_CENTRO波，Sun_10波，以及一条人工波的波形；

图5示出了图4中的三条波转化成反应谱后与场地反应谱的比较图；

图6示出了EI_CENTRO波作用下支座的剪力-位移滞回曲线，Sun_10波作用下支座的剪力-位移滞回曲线，以及人工波1作用下支座的剪力-位移滞回曲线；

图7是表达式被转换后形成的图表。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明而不是对本发明的限制。

依据本发明的第一具体实施例，利用以下步骤实现LNG储罐隔震效应的确定：确定LNG储罐的基本信息，计算内罐模型参数，计算外罐模型参数，计算基础参数，建立LNG储罐隔震力学模型，根据储罐所在地的地震反应谱选择地震波，对建立的分析模型进行特定地震波作用下的动力反应分析，从而确定LNG储罐隔震效应。

在确定LNG储罐的基本信息中，在进行计算前，需事先确定下列基本信息，作为计算的基本参数。确定钢质内罐的基本参数，内罐的直径D(m)，内罐半径R(m)，内罐罐壁等效厚度t(mm)，内罐材料弹性模量E(MPa)，液体的密度ρ_L(kg/m³)和液体的设计液位H(m)。确定混凝土外罐的基本参数，罐壁外表面直径D₁(m)，罐壁内表面直径d(m)，罐壁高度h_w(m)，罐壁质量m_w(kg)，罐顶高度h_r(m)，罐顶质量m_r(kg)，底板厚度h_b(m)，底板质量m_b(kg)和混凝土弹性模量E_c(MPa)。确定单桩的水平刚度k_h(N/m)和桩身竖向刚度k_{V_P}(N/m)。确定底板的桩位布置图及桩根数n_P。确定隔震支座的力学参数，弹性刚度K_u(N/m)、屈服后刚度K_d(N/m)、屈服力F_y(N)和竖向刚度k_{V_B}(N/m)。确定储罐所在场地的地震影响系数曲线(反应谱)。

之后，计算内罐模型参数：

$m_L=\frac{\mathrm{\&pi;}{D}^2}{4}{\mathrm{\&rho;}}_{L}H,$ m_L-为液体总质量(kg)；

$m_c=0.230\frac{D}{H}\tanh \left(\frac{3.67H}{D}\right)m_L,$ m_c—为内罐液体晃动对流部分质量(kg)；

$h_c=[1.0-\frac{\cosh \left(\frac{3.67H}{D}\right)-1.937}{\frac{3.67H}{D}\sinh \left(\frac{3.67H}{D}\right)}]H,$ h_c—为晃动对流质量作用点等效高度(考虑了液体对底板的动压力)(m)；

$T_c=\frac{1.0404\sqrt{D}}{\sqrt{\tanh \left(\frac{3.67H}{D}\right)}},$ T_c—液体的晃动对流周期(s)；

$K_c={\left(\frac{2\mathrm{\&pi;}}{T_c}\right)}^2{m}_c,$ K_c—为晃动对流作用等效刚度(N/m)；

$C_c=2{\mathrm{\&xi;}}_c\sqrt{K_c.m_c};$ c_c—为晃动对流液体的阻尼(N.s/m)，ξ_c为液体晃动阻尼比；

$m_i=\left\{\begin{array}{cc}\frac{\tanh \left(0.866\frac{D}{H}\right)m_L}{0.866\frac{D}{H}}& D/H\&GreaterEqual;1.333\\ (1.0-0.218\frac{D}{H})m_L& D/H<1.333\end{array}\right.;$ m_i—为内罐液体冲击部分质量(kg)；

$h_i=0.375[1.0+1.333[\frac{0.866\frac{D}{H}}{\tanh \left(0.866\frac{D}{H}\right)}-1.0\left]\right]H$ 当D/H≥1.333

$h_i=[0.500+0.060\frac{D}{H}]H$ 当D/H＜1.333，h_i—为冲击质量作用点等效高度(考虑了液体对底板的动压力)(m)；

$T_i=C_i\frac{h\sqrt{\mathrm{\&rho;}}}{\sqrt{t/R}\&CenterDot;\sqrt{E}},$ $C_i=\left(\frac{2}{\sqrt{H/R}\&times;(0.46-0.15H/R+0.017{(H+/R)}^2)}\right),$ T_i—液体的冲击周期(s)；

$K_i={\left(\frac{2\mathrm{\&pi;}}{T_i}\right)}^2{m}_i,$ K_i—为冲击作用等效刚度(N/m)；

$C_i=2{\mathrm{\&xi;}}_i\sqrt{K_i.m_i};$ c_i—为冲击液体的阻尼(N.s/m)，ξ_i为液体冲击阻尼比；

下一步骤是计算外罐模型参数，具体为：

m_t＝m_w+m_r，m_t—罐壁和罐顶的总质量(kg)；

$h_t=\frac{0.5m_w{h}_w+m_r(h_w+\frac{1}{3}{h}_r)}{m_w+m_r},$ h_t—罐壁和罐顶的等效作用高度(m)；

$I=\frac{\mathrm{\&pi;}({D_1}^4-d^4)}{64},$ I—罐壁惯性矩(m⁴)；

$K_t=\frac{3E_{c}I}{h_w^3},$ K_t—为外罐的等效刚度(N/m)；

$C_t=2{\mathrm{\&xi;}}_t\sqrt{K_t.m_t},$ C_t—为外罐的阻尼(N.s/m)，ξ_t为混凝土的阻尼比；接着的步骤是计算基础参数，具体为：

$I_p=\underset{i=1}{\overset{n_p}{\mathrm{\&Sigma;}}}{x}_i^2,$ I_P—为桩群对底板中心的惯性矩(m²)，x_i—为桩心到底板中心的距离；

K_rot＝K_{V_P}I_p，K_rot—为桩群的转动刚度(N.m/rad)；

K_bv=min(K_{V_B}，K_{V_P})，K_bv—为基础竖向刚度和隔震支座的竖向刚度二者间的较小值；

c_bv=0，c_bv—为桩的竖向阻尼；

然后建立LNG储罐隔震力学模型，具体地利用上面的模型参数以及基础参数建立力学分析模型，并根据所述力学分析模型在结构软件上建立分析模型。

根据储罐所在地的地震反应谱选择地震波，依据建筑抗震设计规范(GB 50011-2001)的5.1.2条的规定，正确选择输入的地震加速度时程曲线，该时程曲线要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。对建立的分析模型进行特定地震波作用下的动力反应分析，得到分析结果。

本发明的第二实施例是依据第一实施例的LNG储罐隔震效应确定方法的一个计算实例。

首先确定LNG储罐的基本信息：

内罐的直径D=80m，内罐半径R=40m，内罐罐壁等效厚度t=20mm，内罐材料弹性模量E=200000MPa，液体的密度ρ_L=480kg/m³和液体的设计液位H=34.26m。

接着确定混凝土外罐的基本参数，罐壁外表面直径D₁=83.6m，罐壁内表面直径d=82m，罐壁高度h_w=38.55m，罐壁质量m_w=2.065E+07kg，罐顶高度h_r=11.375m，罐顶质量m_r=0.685E+07kg，底板厚度h_b=0.9m，底板质量m_b=1.415E+07kg和混凝土弹性模量E_c=3450MPa。

之后确定单桩的水平刚度k_h和桩身竖向刚度k_{V_P}；k_h＝29MN/m，k_{V_P}=600MN/m；并确定底板的桩位布置图及桩根数n_P；n_P=360，桩位布置图如图1所示。

接着确定隔震支座的力学参数，弹性刚度K_u(N/m)、屈服后刚度K_d(N/m)、屈服力F_y(N)和竖向刚度k_{V_B}(N/m)；K_u=14MN/m，K_d=1.4MN/m，F_y=80kN，k_{V_B}=1400MN/m。以及确定储罐所在场地的地震影响系数曲线(反应谱)。该场地的地震影响系数曲线(反应谱)如下式和图表所示：

S_a(T)＝A_max·β(T)

$\mathrm{\&beta;}\left(T\right)=\left\{\begin{array}{cc}1& T\&le;0.04s\\ 1+({\mathrm{\&beta;}}_{\max }-1)\frac{T-0.04}{T_1-0.04}& 0.04s<T\&le;T_1\\ {\mathrm{\&beta;}}_{\max }& T_1<T\&le;T_2\\ {\mathrm{\&beta;}}_{\max }{(T_2/T)}^{\mathrm{\&gamma;}}& T_2<T\&le;8.0s\end{array}\right.$

具体参数如下表

可以把上述表达式转化为图表形式，所述图表在图7中示出。

随后计算内罐模型参数。

$m_L=\frac{\mathrm{\&pi;}{D}^2}{4}{\mathrm{\&rho;}}_{L}H=\frac{\mathrm{\&pi;}\&times;{80}^2}{4}\&times;480\&times;34.26=8.266E+07\mathrm{kg}$

$m_c=0.230\frac{D}{H}\tanh \left(\frac{3.67H}{D}\right)m_L=4.072E+07\mathrm{kg}$

$h_c=[1.0-\frac{\cosh \left(\frac{3.67H}{D}\right)-1.937}{\frac{3.67H}{D}\sinh \left(\frac{3.67H}{D}\right)}]H=28.826m$

$T_c=\frac{1.0404\sqrt{D}}{\sqrt{\tanh \left(\frac{3.67H}{D}\right)}}=9.716s$

$K_c={\left(\frac{2\mathrm{\&pi;}}{T_c}\right)}^2{m}_c=1.701E+07N/m$

$C_c=2{\mathrm{\&xi;}}_c\sqrt{K_c.m_c}=2\&times;0.005\sqrt{1.701E+07\&times;4.072E+07}=2.632E+05N.s/m$

$m_i=\frac{\tanh \left(0.866\frac{D}{H}\right)m_L}{0.866\frac{D}{H}}=3.95E+07\mathrm{kg}$

$h_i=0.375[1.0+1.333[\frac{0.866\frac{D}{H}}{\tanh \left(0.866\frac{D}{H}\right)}-1.0\left]\right]H=31.588m$

$C_i=\left(\frac{2}{\sqrt{H/R}\&times;(0.46-0.15H/R+0.017{(H/R)}^2)}\right)=6.282$

$T_i=C_i\frac{H\sqrt{\mathrm{\&rho;}}}{\sqrt{t/R}\&CenterDot;\sqrt{E}}=0.472s$

$K_i={\left(\frac{2\mathrm{\&pi;}}{T_i}\right)}^2{m}_i=7.00E+09N/m$

$C_i=2{\mathrm{\&xi;}}_i\sqrt{K_i.m_i}=2\&times;0.05\&times;\sqrt{(9.0E+09)\&CenterDot;(3.95E+07)}=5.257E+07N.s/m$

然后计算外罐模型参数。

m_t＝m_w+m_r=(2.065E+07)+(0.685E+07)=2.75E+07kg

$h_t=\frac{0.5m_w{h}_w+m_r(h_w+\frac{1}{3}{h}_r)}{m_w+m_r}=25.02m$

$I=\frac{\mathrm{\&pi;}({D_1}^4-d^4)}{64}=1.84E+05m^4$

$K_t=\frac{3E_{c}I}{h_w^3}=3.323E+10N/m$

$C_t=2{\mathrm{\&xi;}}_t\sqrt{K_t.m_t}=2\&times;0.05\&times;\sqrt{(3.323E+10)\&CenterDot;(2.75E+07)}=9.56E+07N.s/m.$

接着计算基础参数。

$I_p=\underset{i=1}{\overset{n_p}{\mathrm{\&Sigma;}}}{x}_i^2=183473m^2$

K_rot＝K_{V_P}I_p=600MN/mx183473=1.1E+14N.m/rad

K_bv=min(K_{V_B}，K_{V_P})=min(1400MN/m，600MN/m)=600MN/m

c_bv=0

在上述各步骤的基础上建立LNG储罐隔震力学模型，并且在图2中示出了LNG储罐隔震力学模型。本次计算采用商业软件SAP2000来进行，其中的计算模型如图3所示。

其中，节点2，3和4的质量分别为m_t，m_i和m_c，节点1的质量为m_b；节点5和6为构造点，水平位置在节点1的左右侧，竖向高度和节点1一样，用定义节点束缚的方式把节点5，6和1约束在一起；节点7也为构造点，位置在节点1的正下方，与节点1的距离为支座顶面和底板中心线的距离。在节点7处“绘制1节点连接”单元8，单元8的“连接/支座类型”为橡胶支座(rubber isolator)；单元3，4和5为弹簧单元，用来模拟k_t,k_i和k_c。

并根据储罐所在地的地震反应谱选择地震波。根据筛选结果，最后选择的地震波为2条实际发生的地震波：EI_CENTRO和Sun_10波，另外再加一条人工波，具体图形如图4所示。并且这三条波转化成反应谱后与场地反应谱的比较如图5所示。

把人工波1、EI_CENTRO波和Sun_10波输入建立的分析模型，进行特定地震波作用下的动力反应分析，得到分析结果。其中在图6中示出上述三条地震波作用下支座的剪力-位移滞回曲线。

把上面三条地震波作用的计算结果汇总如下表：

非隔震情况下，支座的计算结果如下表所示：

设置隔震支座的计算结果与非隔震的计算结果比较如下表所示：

由上表可知，在设置隔震支座后，桩顶的剪力和弯矩大幅度减小（约80%左右），桩顶的受力得到了很大的改善，储罐的抗震性能得到了很大的提高。

另外，本发明的LNG储罐隔震效应的确定方法在唐山LNG项目中得到了成功的实施和应用。

尽管参照优选实施例已经描述了本发明，本领域熟练技术人员将认识到，可以进行形式和细节上的改变，只要不脱离本发明的精神和范围。本发明试图不局限于被公开的具体实施例，如预期用于实施本发明的最佳模式，相反，本发明将包括落入附加权利要求的范围的全部实施例。

Claims (1)

1.一种LNG储罐隔震效应的确定方法，其特征在于步骤如下：

(1)确定LNG储罐的基本信息：

a)确定钢质内罐的基本参数，利用第一长度测量装置确定内罐的直径D(m)，并确定内罐半径R(m)，内罐罐壁等效厚度t(mm)，内罐材料弹性模量E(MPa)，液体的密度ρ_L(kg/m³)和液体的设计液位H(m)这些参数；

b)确定混凝土外罐的基本参数，利用第二长度测量装置确定罐壁外表面直径D₁(m)、罐壁内表面直径d(m)以及罐壁高度h_w(m)，并确定罐壁质量m_w(kg)，罐顶高度h_r(m)，罐顶质量m_r(kg)，底板厚度h_b(m)，底板质量m_b(kg)和混凝土弹性模量E_c(MPa)；

c)确定单桩的水平刚度k_h(N/m)和桩身竖向刚度k_{V_P}(N/m)；

d)确定底板的桩位布置图及桩根数n_P；

e)确定隔震支座的力学参数，弹性刚度K_u(N/m)、屈服后刚度K_d(N/m)、屈服力F_y(N)和竖向刚度k_{V_B}(N/m)；

f)确定储罐所在场地的地震影响系数曲线，该地震影响系数曲线即一反应谱；

(2)计算内罐模型参数

$m_L=\frac{\mathrm{\&pi;}{D}^2}{4}{\mathrm{\&rho;}}_{L}H,$ m_L—为液体总质量(kg)；

$m_c=0.230\frac{D}{H}\tanh \left(\frac{3.67H}{D}\right)m_L,$ m_c—为内罐液体晃动对流部分质量(kg)；

$h_c=[1.0-\frac{\cosh \left(\frac{3.67H}{D}\right)-1.937}{\frac{3.67H}{D}\sinh \left(\frac{3.67H}{D}\right)}]H,$ h_c—为晃动对流质量作用点等效高度；

$T_c=\frac{1.0404\sqrt{D}}{\sqrt{\tanh \left(\frac{3.67H}{D}\right)}},$ T_c—液体的晃动对流周期(s)；

$K_c={\left(\frac{2\mathrm{\&pi;}}{T_c}\right)}^2{m}_c,$ K_c—为晃动对流作用等效刚度(N/m)；

$C_c=2{\mathrm{\&xi;}}_c\sqrt{K_c.m_c};$ c_c—为晃动对流液体的阻尼(N.s/m)，ξ_c为液体晃动阻尼比；

$m_i=\left\{\begin{array}{cc}\frac{\tanh \left(0.866\frac{D}{H}\right)m_L}{0.866\frac{D}{H}}& D/H\&GreaterEqual;1.333\\ (1.0-0.218\frac{D}{H})m_L& D/H<1.333\end{array}\right.;$ m_i—为内罐液体冲击部分质量(kg)；

$h_i=0.375[1.0+1.333[\frac{0.866\frac{D}{H}}{\tanh \left(0.866\frac{D}{H}\right)}-1.0\left]\right]H$ 当D/H≥1.333

$h_i=[0.500+0.060\frac{D}{H}]H$ 当D/H＜1.333，h_i—为冲击质量作用点等效高度(考虑了液体对底板的动压力)(m)；

$T_i=C_i\frac{h\sqrt{\mathrm{\&rho;}}}{\sqrt{t/R}\&CenterDot;\sqrt{E}},$ $C_i=\left(\frac{2}{\sqrt{H/R}\&times;(0.46-0.15H/R+0.017{(H+/R)}^2)}\right),$ T_i—液体的冲击周期(s)；

$K_i={\left(\frac{2\mathrm{\&pi;}}{T_i}\right)}^2{m}_i,$ K_i—为冲击作用等效刚度(N/m)；

$C_i=2{\mathrm{\&xi;}}_i\sqrt{K_i.m_i};$ c_i—为冲击液体的阻尼(N.s/m)，ξ_i为液体冲击阻尼比；

(3)计算外罐模型参数

m_t＝m_w+m_r，m_t—罐壁和罐顶的总质量(kg)；

$h_t=\frac{0.5m_w{h}_w+m_r(h_w+\frac{1}{3}{h}_r)}{m_w+m_r},$ h_t—罐壁和罐顶的等效作用高度(m)；

$I=\frac{\mathrm{\&pi;}({D_1}^4-d^4)}{64},$ I—罐壁惯性矩(m⁴)；

$K_t=\frac{3E_{c}I}{h_w^3},$ K_t—为外罐的等效刚度(N/m)；

$C_t=2{\mathrm{\&xi;}}_t\sqrt{K_t.m_t},$ C_t—为外罐的阻尼(N.s/m)，ξ_t为混凝土的阻尼比；(4)计算基础参数

$I_p=\underset{i=1}{\overset{n_p}{\mathrm{\&Sigma;}}}{x}_i^2,$ I_P—为桩群对底板中心的惯性矩(m²)，x_i—为桩心到底板中心的距离；

K_rot＝K_{V_P}I_p，K_rot—为桩群的转动刚度(N.m/rad)；

K_bv=min(K_{V_B}，K_{V_P})，K_bv—为基础竖向刚度和隔震支座的竖向刚度二者间的较小值；

c_bv=0，c_bv—为桩的竖向阻尼；

(5)建立LNG储罐隔震力学模型

利用上面的模型参数以及基础参数建立力学分析模型，并根据所述力学分析模型在结构软件上建立分析模型。

(6)根据储罐所在地的地震反应谱选择地震波

依据建筑抗震设计规范的规定，正确选择输入的地震加速度时程曲线，该时程曲线要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。

(7)对建立的分析模型进行特定地震波作用下的动力反应分析，从而确定LNG储罐隔震效应。

Images