• ISSN 1673-5722
  • CN 11-5429/P

基于DOE考虑结构参数不确定性的砌体结构地震易损性分析

周强 李长亮 赵文洋

刘洪涛,孔鹏超,王作虎,廖维张,2022. 钢筋混凝土预制拼装柱扭转力学性能数值模拟与参数分析. 震灾防御技术,17(2):372−380. doi:10.11899/zzfy20220217. doi: 10.11899/zzfy20220217
引用本文: 周强,李长亮,赵文洋,2024. 基于DOE考虑结构参数不确定性的砌体结构地震易损性分析. 震灾防御技术,19(3):601−612. doi:10.11899/zzfy20240318. doi: 10.11899/zzfy20240318
Liu Hongtao, Kong Pengchao, Wang Zuohu, Liao Weizhang. Numerical Simulation and Parameter Analysis of Torsion Mechanical Properties of Reinforced Concrete Precast Assembly Columns[J]. Technology for Earthquake Disaster Prevention, 2022, 17(2): 372-380. doi: 10.11899/zzfy20220217
Citation: Zhou Qiang, Li Changliang, Zhao Wenyang. Seismic Vulnerability Analysis of Masonry Structures with Uncertain Parameters Based on DOE[J]. Technology for Earthquake Disaster Prevention, 2024, 19(3): 601-612. doi: 10.11899/zzfy20240318

基于DOE考虑结构参数不确定性的砌体结构地震易损性分析

doi: 10.11899/zzfy20240318
基金项目: 国家自然科学基金(51968047、51608249)
详细信息
    作者简介:

    周强,男,生于1983年。博士,副教授,硕士生导师。主要从事结构抗震及灾场仿真研究。E-mail:zhouqiang@ncu.edu.cn

Seismic Vulnerability Analysis of Masonry Structures with Uncertain Parameters Based on DOE

  • 摘要: 砌体结构参数具有离散性高且其受到地震作用后响应结果非线性强等特点,所以在砌体地震易损性研究中结构参数不确定性对结果产生的影响不容忽视。针对汶川地震中具有代表性的砌体结构,提出基于DOE(试验设计)方法考虑在砌体结构地震易损性研究中结构参数不确定性的影响。首先,使用Plackett-Burman方法开展结构参数的灵敏度分析,筛选出对砌体结构地震响应影响较大的3个结构参数;然后,根据筛选结果及地震参数PGA进行试验设计,进而建立结构参数与最大层间位移角的响应面回归模型;最后,通过进行蒙特卡罗模拟获得地震易损性曲线,并进一步评估结构参数不确定性对砌体地震易损性分析的影响程度。研究结果表明,弹性模量、密度和抗拉强度是对砌体结构易损性影响较大的3个参数;当地震动PGA为0.1 g时,相较于其他结构参数,弹性模量对结构地震响应的影响最显著;当PGA为0.2 g时,3个结构参数的影响由大到小依次为弹性模量、密度和抗拉强度,当PGA为0.4 g时,3个结构参数的影响由大到小依次为弹性模量、抗拉强度与密度。本文所提出的方法计算量少、精度高,可为有效解决砌体结构参数种类众多和材料非线性强等难题提供新思路。
  • 近年来,我国经济高速发展的同时能源需求持续增长,发展核能是确保我国能源供应安全的有力保障,同时全球自然灾害频发也对核能应用安全提出了更高要求,由于我国位于环太平洋地震带与欧亚地震带的交汇地带,核电结构地震安全性评价十分必要,除要求考虑土-结构相互作用外,还需考虑周围邻近结构的影响(Lou等,2011王国波等,2018),即结构-土-结构相互作用。

    相邻结构之间的相互作用是基于土-结构相互作用研究成果的拓展,当介质中存在多个结构时,由于结构响应引起土体扰动,并通过土体传播,土-结构问题演变为多个结构之间的交叉相互作用问题,Whitman(1969)首次提出基础通过土体相互耦合的问题,前期研究土-结构系统模型多采用弹性地基上单质量系统模型。Warburton等(1971)将地基考虑为弹性半空间,相邻上部结构简化为圆柱质量块,基础假设为刚性圆形基础,研究结果表明当1个质量块被外部谐波激发时,第2个质量块通过相对较小的扰动改变受激质量的位移分量。支承在柔性基础上的结构动力响应与支承在刚性基础上的结构动力响应有很大不同(Lehmann等,2001李培振等,2014),刚性基础的假设与震时系统实际运动状态存在一定差距,产生差异的重要原因之一是结构振动产生的能量部分通过柔性基础耗散。田彼得等(1987)基于频域采用子结构法对相邻结构动力相互作用体系动力特性进行了分析,研究结果表明高大结构的存在使低矮结构顶层相对位移频响曲线的第1共振频率显著降低。目前基于频域仅可采用线性方法或等效线性化近似考虑非线性问题,而地震波激励下结构一般会进入非线性状态。柳玉印等(2018)对相邻新、旧重力坝的彼此相互作用进行模拟,结果表明相邻结构的存在对结构峰值响应有一定影响。Luco等(1973)、Murakami等(1977)及Wong等(1975)用简化的二维平面应变模型模拟无限长三维结构,研究结果表明在低频情况下,对于邻近较大结构物的较小结构物而言,结构-土-结构相互作用(Structure-Soil-Structure Interaction,SSSI)的影响较大,而高频情况下可忽略邻近结构,仅考虑单体结构;地震波能量主要集中在低频区域,可见在地震波激励下相邻高层建筑对底层建筑相邻结构的相互作用影响较大。但以上研究简化了上部结构模型,忽略了上部结构尺寸和土体之间相互作用的影响。Ghandil等(2016)利用等效线性模型近似考虑了土体中的大剪切应变,通过研究3个不同高度相邻建筑的相互作用,发现当建筑物间距为建筑基底宽度的一半时,相互作用不可忽视。Aldaikh等(20152016)探讨了地震波激励下SSSI对3个相邻建筑的影响,提出了3个相邻建筑通过土壤耦合的理论公式,并通过振动台试验证明该简单理论模型在描述线弹性半空间上相邻建筑群的动力相互作用是有效的。何涛等(2020)依托振动台试验并结合子结构试验技术,以2个4层钢框架为例,模拟了震时相邻结构的加速度放大系数,变化规律,研究结果表明考虑相邻结构影响时的震时地震响应峰值较单体结构小。Wang等(2022)针对不同尺度的土-结构相互作用(Soil-Structure Interaction,SSI)系统和SSSI系统进行了振动台试验,结果表明SSSI效应对结构反应的影响随着结构间距的增大而减小,随着结构尺寸的增大和地震波激励方向的改变而增大或减小。韩冰等(2019)采用间接边界元法研究了SSSI效应对系统频率的影响,并得出随着结构间距的增大相邻结构相互作用对结构系统频率的影响呈下降趋势的结论。以上研究充分说明了震时相邻结构的存在在一定程度上改变了原有结构的动力特性,且入射波角度的改变可能导致相邻结构系统处于最不利工况状态。

    反应堆厂房通常不是独立于附属结构而存在,受限于核能发电的工艺要求,反应堆厂房周围需建造一系列辅助厂房,以维持整个核电站的正常运作,如果SSSI效应影响较大,应将其作为核电站抗震设计和安全评估的重要特性之一。近年来,关于核电站SSSI的研究发展迅速,Lee等(1973a1973b)采用三维结构模型研究了SSSI效应对相邻反应堆厂房地震响应的影响,并提出了近似解析-数值方法,以求解弹性半空间表面上3个刚性圆形基础的相互作用。Kitada等(1999)对邻近反应堆厂房进行了一系列现场原位试验,研究发现SSSI的影响与外部激励荷载方向和相邻结构空间位置密切相关。Matthees等(1982)对核电站相邻结构在水平地震作用下相互作用的灵敏度进行了研究,结果表明相互作用现象对核电结构响应的影响是不可忽视的。上述研究中相邻结构的基础是分离的,考虑核电结构特殊性及其区别于一般民用建筑的空间布置,核电结构为增强整体性,减少地震作用下邻近基础的相互影响,设计时将核电附属厂房与反应堆厂房共用1个基础,应稼年等(1995)在同一基础上考虑了2个辅助厂房与安全壳在地震波激励下的平-扭耦联动力响应,并指出邻近厂房的存在会增加反应堆厂房震时响应幅值。

    考虑相邻厂房对反应堆厂房的影响,即考虑SSSI效应,计算量较大,因此已有研究多将结构模型简化,若对结构进行较详细地建模分析,需采用高效的计算方法。本文在核电结构中土-结构相互作用分区分析方法(PASSI)的基础上,提出显-隐式单元层计算方法,可在满足高效性的同时,提高稳定性。以我国某核电站为研究对象,其4个辅助厂房同反应堆厂房建造在同个基础上,设计了4个工况模型,分析了周围辅助厂房对反应堆厂房自振频率及地震响应的影响。

    SSI系统分析模型如图1所示。该系统由两部分组成,分别为结构子系统和土体子系统,结构子系统包含上部结构、基础和部分土体,土体子系统为全部土体区域,其与结构子系统存在重叠区域。

    图 1  SSI系统分析模型
    Figure 1.  Soil-structure interaction system

    将土体离散成八节点六面体单元,并将节点划分为内部点和人工边界点,如图1所示。内部点采用集中质量有限元法计算(廖振鹏,2002),动力方程的一般形式可写为:

    $$ {{\boldsymbol{M}}_i}{\ddot {\boldsymbol{u}}_i} + \sum\limits_j^N {{{\boldsymbol{C}}_{ij}}{{\dot {\boldsymbol{u}}}_i}} + \sum\limits_j^N {{{\boldsymbol{K}}_{ij}}{{\boldsymbol{u}}_i}} = {{\boldsymbol{F}}_i} $$ (1)

    式中,$ {{\boldsymbol{M}}_i} $为集中在内部点i上的质量阵,具体形式如式(2)所示;$ {{\boldsymbol{C}}_{ij}} $$ {{\boldsymbol{K}}_{ij}} $分别为内部点i与相邻节点j之间的阻尼阵和刚度阵;$ {\ddot {\boldsymbol{u}}_i} $$ {\dot {\boldsymbol{u}}_i} $$ {{\boldsymbol{u}}_i} $分别为内部点i的加速度向量、速度向量、位移向量;N为与节点i相邻节点总数;$ {F_{{i}}} $为内部点i的外部荷载向量。

    $$ {{\boldsymbol{M}}_i}{\text{ = }}\left[ {\begin{array}{*{20}{c}} {{m_i}}&0&0 \\ 0&{{m_i}}&0 \\ 0&0&{{m_i}} \end{array}} \right] $$ (2)

    对式(1)采用以下时步积分格式离散:

    $$ \dot {\boldsymbol{u}}_i^p = \frac{{{\boldsymbol{u}}_i^p - {\boldsymbol{u}}_i^{p - 1}}}{{\Delta t}} $$ (3)
    $$ \ddot {\boldsymbol{u}}_i^p = \frac{{{\boldsymbol{u}}_i^{p + 1} - 2{\boldsymbol{u}}_i^p + {\boldsymbol{u}}_i^{p - 1}}}{{\Delta {t^2}}} $$ (4)

    式中,$ {\boldsymbol{u}}_i^{p + 1}、{\boldsymbol{u}}_i^p、{\boldsymbol{u}}_i^{p - 1} $分别为i节点$ p + 1、p、p - 1 $时刻的位移;$ \Delta t $为时间步长;$ \dot {\boldsymbol{u}}_i^p、\ddot {\boldsymbol{u}}_i^p $分别为i节点$ p $时刻速度和加速度。

    进而得到:

    $$ {\boldsymbol{u}}_i^{p + 1} = 2{\boldsymbol{u}}_i^p - {\boldsymbol{u}}_i^{p - 1} - \frac{{\Delta t}}{{{{{m}}_i}}}\sum\limits_e {{\boldsymbol{c}}_i^e} ({\boldsymbol{u}}_i^p - {\boldsymbol{u}}_i^{p - 1}) - \frac{{\Delta {t^2}}}{{{m_i}}}\sum\limits_e {{\boldsymbol{k}}_i^e} {\boldsymbol{u}}_i^p + \frac{{\Delta {t^2}}}{{{m_i}}}{{\boldsymbol{F}}_i} $$ (5)

    采用刘晶波等(20052022)提出的黏弹性人工边界,在边界上设置切向、法向弹簧和阻尼,按式(6)、式(7)计算:

    $$ {K_{{\rm{BT}}}} = {\alpha _{\rm{T}}}\frac{G}{R} ,\qquad {C_{{\rm{BT}}}} = \rho {C_{\text{S}}} $$ (6)
    $$ {K_{{\rm{NT}}}} = {\alpha _{\rm{N}}}\frac{G}{R} ,\qquad {C_{{\rm{NT}}}} = \rho {C_{\rm{P}}} $$ (7)

    式中,G为场地的剪切模量;$ {K_{{\rm{BT}}}}、{K_{{\rm{NT}}}} $为别切向和法向弹簧刚度;$ {\alpha _{\rm{T}}}、{\alpha _{\rm{N}}} $分别为粘弹性人工边界的切向和法向修正系数;R为波源点至人工边界的距离;$ {C_{{\rm{BT}}}}、{C_{{\rm{NT}}}} $分别为切向阻尼和法向系数;ρ为介质密度;$ {C_{\rm{S}}}、{C_{\rm{P}}} $分别为介质的剪切波和压缩波波速。

    在人工边界上施加地震动输入的等效节点力为:

    $$ {{\boldsymbol{F}}_b} = ({{\boldsymbol{\sigma}} _b} + c{\dot {\boldsymbol{u}}_b} + k{{\boldsymbol{u}}_b})A $$ (8)

    式中,$ {{\boldsymbol{\sigma}} _b} $为边界点上b点在t时刻的自由场应力;ck分别为式(6)、式(7)中的弹簧和阻尼系数;$ {\dot {\boldsymbol{u}}_b} $$ {{\boldsymbol{u}}_b} $分别为人工边界处自由场速度向量和位移向量;A为边界点影响面积。

    结构子系统运动方程为:

    $$ {\boldsymbol{M}}\ddot u + {\boldsymbol{C}}\dot u + {\boldsymbol{K}}u = {\boldsymbol{F}} $$ (9)

    采用Newmark隐式算法进行瞬态动力直接积分得:

    $$ ({b_1}{\boldsymbol{M}} + {b_4}{\boldsymbol{C}} + {\boldsymbol{K}}){{\boldsymbol{u}}^{p + 1}} = {{\boldsymbol{F}}^{p + 1}} + {\boldsymbol{M}}({b_1}{{\boldsymbol{u}}^p} - {b_2}{\dot {\boldsymbol{u}}^p} - {b_3}{\ddot {\boldsymbol{u}}^p}) + {\boldsymbol{C}}({b_4}{{\boldsymbol{u}}^p} - {b_5}{\dot {\boldsymbol{u}}^p} - {b_6}{\ddot {\boldsymbol{u}}^p}) $$ (10)
    $$ {b_1} = 1/(\beta \Delta {t^2}),{b_2} = 1/(\beta \Delta t),{b_3} = 1 - 1/(2\beta ),{b_4} = \gamma \Delta t{b_1},{b_5} = 1 + \gamma \Delta t{b_2},{b_6} = \Delta t(1 + \gamma {b_3} - \gamma ) $$ (11)

    求解式(10)可得结构上各节点反应。

    在PASSI中,显、隐式域不重叠,采用交互力和位移的方式进行耦合,当应用于柔性基础时,稳定性欠佳。设置显-隐式单元层,该单元层既属于显式单元,又属于隐式单元,为显、隐式单元的重叠区,通过交互位移的方式实现耦合分析(Belytschko等,1978)。结构子系统采用隐式分析,为隐式单元;土体子系统采用显式分析,为显式单元;结构子系统与土体子系统的重叠区域为显-隐式单元层。结构子系统边界点(图1中黑色实心、红色外圆点)同时也是土体子系统的内部节点,同样,土体子系统边界点(图1中红色实心、黑色外圆点)也是结构子系统的内部节点。以一维网格为例,说明显-隐式交替计算的实现,如图2(a)所示,将节点分为2个集合,即显式计算节点和隐式计算节点,分别用黑色和红色实心圆点表示;将单元分为3个集合,即显式单元(E)、隐式单元(I)和显-隐式单元(EI),显-隐式单元为显式单元和隐式单元重叠区,图中红色实心、黑色外圆点既是显式单元区的边界点,又是隐式单元区的内部点,计算中按隐式求解;黑色实心、红色外圆点既是隐式单元区的边界点,又是显式单元区的内部点,计算中按显式求解。显式计算点的网格信息流用黑色实线表示,隐式计算点的信息流用红色实线表示。假设 $ t $ 和($ t - \Delta t $)时刻系统反应已知,可由式(5)将显式计算点(含隐式单元区的边界点,图2(a)中的J−1点)更新到($ t + \Delta t $)时刻,显式内部点J−1的位移传递至隐式边界点J−1,可由式(10)计算得到隐式节点($ t + \Delta t $)时刻的响应,再将隐式单元内部点(图2(a)中的J节点)的反应传递至显式单元区,从而将整个系统 $ t $ 时刻响应更新到($ t + \Delta t $)时刻,重复该过程,即可实现系统的时程分析。

    图 2  显-隐式数据交互示意
    Figure 2.  Explicit-implicit data interaction

    当隐式部分时间步长为显式部分时间步长的NE−I倍时,在隐式时间步长内,信息可在显式网格中跨越NE−I个节点,那么EI单元层须包含NE−I层单元。以图2(b)所示N=5为例(Belytschko等,19791985),需在时间步($ t + 5\Delta t $)内更新J−5之后的所有显式节点,求显式域$ u_{J - 5}^{t + 5\Delta t} $时,$ u_{J - 6}^{t + 4\Delta t} $$ u_{J - 5}^{t + 4\Delta t} $$ u_{J - 4}^{t + 4\Delta t} $$ u_{J - 6}^{t + 3\Delta t} $$ u_{J - 5}^{t + 3\Delta t} $$ u_{J - 4}^{t + 3\Delta t} $需已知,根据图2(b),假设t为0时刻,系统0时刻及前一时刻响应已知,通过式(5)将显式域 $ t $ 时刻响应更新到($ t + \Delta t $)时刻,那么显式域($ t + 2\Delta t $)时刻响应可由($ t + \Delta t $)、$ t $时刻求得,按此递推方法,($ t + 3\Delta t $)和($ t + 4\Delta t $)时刻系统响应可分别由各自时刻的前2个时刻推得,最终可求得显式域($ t + 5\Delta t $)时刻系统响应。首先将显-隐式单元层分配给显式,显式先集成分区,通过J−5更新隐式边界点;然后将隐式单元层利用显式的输入作为边界条件,通过式(10)求解,并将隐式域节点J−4到J节点的信息传输至显式域,因显式域是根据系统2个时刻响应得到下一时刻响应,所以隐式域需连续计算2个时刻,同时隐式域计算每个时刻节点J−4到J节点的信息均需传输至显式域,以保证显式域节点J−4到J节点的2个时刻系统响应已知,作为显式域的初始条件计算显式域的系统响应。重复以上步骤可获得整个系统的响应。

    基于PASSI的显-隐式单元层计算方法计算流程如下:

    (1)设置人工边界,采用上文介绍的黏弹性边界;

    (2)输入脉冲及地震波,采用式(8)等效节点力的方式输入;

    (3)采用集中质量有限元法,由式(5)得到显式域土体节点位移;

    (4)求得显-隐式交界点位移,由式(5)得到交界点(p+1)时刻位移$ \left\{ {U_{Jp + 1}^E} \right\} $,并作为隐式结构子系统的边界位移输入;

    (5)获得结构子系统的边界条件,通过Fortran程序调用ANSYS软件进行动力时程分析,采用Newmark隐式计算方法,由式(10)得到隐式节点(p+1)时刻的位移,并将其更新到显式域节点上;

    (6)重复执行步骤(1)~(5),得到系统任意时刻的响应。

    为验证上述方法的有效性,设计了场地分析和土-结构相互作用分析算例,分别采用本文显-隐式交互算法和ABAQUS软件进行分析。在土体模型底边界垂直入射单位脉冲位移波,其中脉冲宽度为0.25 s,时间步距为0.001 s,计算步数为1 000,计算时长为1.0 s,脉冲波位移时程及频谱如图3所示。

    图 3  脉冲波位移时程和频谱曲线
    Figure 3.  Displacement time history and spectrum of SV/P pulse wave

    首先从三维半无限空间中截取$40 \; {\rm{m}}\; { \times }\; 40\; {\rm{m}}\; { \times }\; 50 \;{\rm{m}}$的有限土体,顶部自由,采用$1\; {\rm{m}}\; { \times }\; 1\; {\rm{m}}\; {\times }\; 1 \;{\rm{m}}$的单元离散,在其底部及4个侧面施加黏弹性人工边界,设置1层土,土体材料参数如表1所示,场地监测点如图4(b)所示分别取场地中心的顶部点a、中部点b和底部点c

    表 1  土体参数
    Table 1.  Soil parameters
    材料厚度/m弹性模量/GPa泊松比密度/kg·m−3剪切波速/m·s−1压缩波速/m·s−1
    软土500.1080.351 000200416
    下载: 导出CSV 
    | 显示表格
    图 4  ABAQUS软件中黏弹性人工边界及监测点示意
    Figure 4.  Viscoelastic artificial boundary and monitoring points in ABAQUS

    结合ABAQUS软件二次开发功能(曹金凤等,2011苏景鹤等,2016),采用前文阐述的黏弹性人工边界理论,在ABAQUS软件中实现黏弹性人工边界(刘晶波等,2007王飞等,2018),ABAQUS软件中施加等效弹簧阻尼系统如图4(a)所示,场地监测点如图4(b)所示。

    SV波和P波垂直入射场地位移时程响应如图5所示,实线为显-隐式单元层计算方法计算结果(Chen等,2022),虚线为ABAQUS软件计算结果。由图5可知,2种方法计算结果基本重合,在一定程度上检验了显-隐式单元层计算方法的准确性。

    图 5  SV波和P波垂直入射场地位移时程曲线
    Figure 5.  Site response of SV and P wave

    结构模型采用$1 \;{\rm{m}}\;{ \times }\; 1 \;{\rm{m}}\;{ \times }\; 10 \;{\rm{m}}$的结构柱,用$ 1\; {\rm{m}}\;{ \times }\; 1 \;{\rm{m}}\;{ \times }\; 1 \;{\rm{m}} $的实体单元进行离散,显-隐式单元层计算方法分为两部分,一部分为土体计算,另一部分为结构计算,考虑到核电结构体量大、单元种类多等特点,采用Fortran程序调用ANSYS软件计算上部结构相应。土体模型采用前文所述材料参数和单元尺寸,输入同样的激励,土-结构相互作用模型如图6所示。

    图 6  土-结构相互作用模型
    Figure 6.  Soil-structure interaction model

    结构底部和顶部点位移响应如图7所示,由图7可知,显-隐式单元层计算方法与ABAQUS软件计算得到的结构受迫振动阶段结果基本相符,但结构自由振动阶段的计算结果有所差异,主要因为2种方法在阻尼设置上不同。

    图 7  结构位移时程曲线
    Figure 7.  Displacement time history of nuclear power plant

    本文算例分析在台式机上完成,采用Windows 10专业版操作系统,计算机处理器为Intel(R) Core(TM)i7-9700 CPU @3.00 Hz 3.00 GHz,内存16.0 GB,硬盘500 GB,计算效率如表2所示。

    表 2  显-隐式单元层计算方法和ABAQUS软件计算效率
    Table 2.  PASSI and ABAQUS calculation efficiency
    算例单元数/个节点数/个自由度数/个
    显-隐式单元层计算方法
    计算时间/min
    ABAQUS软件
    计算时间/min
    场地80 00085 731257 1934106
    土-结构相互作用80 01085 771257 31333106
    注:自由度数=节点数*3(每个节点自由度个数)
    下载: 导出CSV 
    | 显示表格

    表2可知,在仅计算场地模型时,ABAQSU软件计算用时为显-隐式单元层计算方法的25倍左右。当考虑土-结构相互作用时,ABAQUS软件计算用时为显-隐式单元层计算方法的3倍左右。对于本文简单的上部结构模型,其计算量可忽略不计,相比于场地模型,ABAQSU软件计算时间基本不变,而显-隐式单元层计算方法每时步需调用ANSYS软件,需重启动,导致计算效率降低。对于简单的上部结构,由于显-隐式单元层计算方法重启动时间大于结构分析时间,因此显-隐式单元层计算方法计算效率相对于ABAQUS软件计算效率提高不明显,但对于大型复杂的核电结构,显-隐式单元层计算方法重启动时间远小于结构分析时间,效率提高更明显。

    综上所述,基于PASSI的显-隐式单元层计算方法与ABAQUS软件计算精度相当,但计算效率较ABAQUS软件高。因此,采用显-隐式单元层计算方法进行辅助厂房对反应堆厂房的地震响应影响分析。

    华龙一号核电站模型由反应堆厂房、燃料厂房、电器厂房和安全厂房组成(图8),4个辅助厂房呈四面环绕反应堆厂房对称布置,其中反应堆厂房设置了双层安全壳,外层安全壳、电器厂房、燃料厂房连成一体,中间不设缝。安全厂房与外层安全壳不连接,辅助厂房各自独立,所有厂房共用同块混凝土基础底板。

    图 8  核电站平面布置示意
    Figure 8.  Floor plan of nuclear power plant

    反应堆厂房模型如图9所示,反应堆厂房监测点分别取外壳顶部点a、内壳顶部点b、外壳中部点c、内壳中部点d、反应堆底板点f、反应堆上部点e

    图 9  反应堆厂房监测点分布示意
    Figure 9.  Monitoring points of reactor plant

    核电站抗震计算模型采用三维有限元模型,需根据结构形式和受力特点,选择适当的单元进行模拟。墙体和楼板采用板壳单元模拟,板壳单元节点位于墙体和楼板厚度中心处;底板由于厚度较大,采用实体单元模拟,材料参数如表3所示。

    表 3  核电站结构材料参数
    Table 3.  Material parameters of nuclear power plant
    编号材料结构弹性模量/GPa泊松比密度/kg·m−3
    1C30厂房/基础32.50.22 400
    2C40内/外壳36.00.22 450
    下载: 导出CSV 
    | 显示表格

    x轴正向为反应堆厂房指向安全厂房方向,y轴正向为反应堆厂房指向电器厂房方向,z轴正向为竖直向上方向,坐标轴原点为反应堆厂房中心。

    地基为弹性半空间,材料参数按表4所示的地质勘察资料取值,场地范围按规范要求(American Society of Civil Engineers(ASCE),2016;《核电厂抗震设计标准》(GB 50267—2019)(中华人民共和国住房和城乡建设部等,2019))选择$340 \;{\rm{m}}\;{ \times }\; 400 \;{\rm{m}}\;{ \times }\;60 \;{\rm{m}}$,将计算区域离散为六面体八节点实体单元,单元尺寸为$2 \;{\rm{m}}\;{ \times } \;2 \; {\rm{m}}\;{ \times }\;2 \;{\rm{m}}$,在其底部及4个侧面边界施加黏弹性人工边界。

    表 4  核电站厂址参数
    Table 4.  Soil parameters of nuclear power plant site
    材料厚度/m弹性模量/GPa泊松比密度/kg·m−3剪切波速/m·s−1压缩波速/m·s−1
    岩石6046.90.262 6502 6734 639
    下载: 导出CSV 
    | 显示表格

    将控制点处(地基表面)地震动幅值调为0.2 g,根据RG160规范谱人工合成的加速度时程如图10(a)所示,持续时间20.48 s,其反应谱如图10(b)所示。从PEER强震数据库中选出RSN5806_IWATE_55461NS(简称RSN5806)、RSN4872_CHUETSU_65053NS(简称RSN4872)地震波,用于RG160地震动输入结构响应规律验证,RSN5806加速度时程如图10(c)所示,持续时间60.00 s,其反应谱如图10(d)所示。RSN4872加速度时程如图10(e)所示,持续时间60.00 s,其反应谱如图10(f)所示。假定地震波垂直入射,由控制点处的地震动,按传递矩阵方法反演得到基岩半空间的入射波,并得到边界处的自由场,按式(8)得到人工边界点等效地震作用。

    图 10  地震动加速度时程与反应谱曲线
    Figure 10.  Ground motion acceleration time history and response spectrum

    为考察相邻附属厂房对反应堆厂房地震响应的影响,考虑以下4种核电站结构模型:不考虑相邻厂房的反应堆厂房有限元模型,为Case1,如图11(a)所示;安全厂房A-反应堆厂房-安全厂房B有限元模型,为Case2,如图11(b)所示;燃料厂房-反应堆厂房-电器厂房有限元模型,为Case3,如图11(c)所示;核电站整体有限元模型,为Case4,如图11(d)所示。

    图 11  核电站4种工况模型
    Figure 11.  Four working conditions model for nuclear power plant

    对各工况进行模态分析,结果如表5所示。由表5可知,Case1和Case2模型前2阶模态频率相同,由于Case2模型中安全厂房与反应堆厂房分离,说明Case2模型前2阶模态是反应堆厂房的模态,整体而言,反应堆厂房的局部模态频率小于安全厂房。由于燃料厂房、电器厂房与反应堆厂房互相连为整体,增强了结构整体刚度,使Case3模型前35阶自振频率整体上大于Case1模型,且Case3模型第1阶自振频率最大。Case2模型第3~5阶频率基本与Case4模型第1~3阶频率相同,这3阶频率是安全厂房的自振频率。Case3模型第1~3阶频率基本与Case4模型第4~6阶频率相同,这是反应堆厂房、燃料厂房和电器厂房作为整体结构的前3阶频率,Case4模型前3阶频率基本与Case2模型第3~5阶频率相同,这3阶频率是安全厂房的自振频率,经燃料厂房和电器厂房加强后,反应堆厂房局部模态频率大于安全厂房频率。Case1和Case2模型第2阶振型参与系数最大,对应的频率均为3.007 Hz左右;Case3模型第1阶振型参与系数最大,Case4模型第4阶振型参与系数最大,对应的频率均为4.017 Hz左右。4种工况中参与系数最大的振型均与反应堆厂房局部振型相关,主要反映反应堆厂房的自振特性。

    表 5  模态分析
    Table 5.  Modality analysis
    振型Case1Case2Case3Case4
    频率/Hz振型参与系数频率/Hz振型参与系数频率/Hz振型参与系数频率/Hz振型参与系数
    13.006 690.186 263.006 690.186 194.016 961.000 003.726 480.426 64
    23.007 491.000 003.007 491.000 004.068 740.277 893.861 540.413 38
    34.078 570.540 833.725 200.677 044.114 500.094 623.885 430.002 48
    44.096 080.381 573.858 000.683 884.210 240.482 594.017 181.000 00
    55.488 440.002 043.885 280.009 754.739 490.029 304.068 750.275 37
    65.513 010.001 434.078 570.540 695.472 460.169 674.114 530.097 12
    75.540 460.003 314.096 080.381 505.566 560.002 174.210 380.480 55
    85.572 030.004 844.559 210.068 575.610 470.012 114.560 350.045 69
    95.629 190.003 165.439 430.064 355.762 070.080 104.739 620.029 71
    105.658 220.010 715.488 440.002 056.163 020.037 355.441 090.042 47
    116.426 340.004 545.513 010.001 436.470 850.006 885.472 660.167 66
    126.471 340.015 755.540 460.003 316.496 260.389 725.566 560.002 12
    136.496 200.598 365.572 030.004 826.615 100.003 045.610 470.011 96
    146.614 210.001 585.629 190.003 176.875 020.134 965.702 480.044 00
    156.875 130.207 355.658 220.010 727.121 850.000 715.762 180.079 50
    167.117 110.000 875.698 290.070 077.183 560.005 666.163 160.037 35
    177.171 630.005 926.426 340.004 557.186 370.071 546.470 850.006 79
    187.186 590.070 526.471 340.015 757.228 460.359 016.496 260.386 15
    197.218 460.014 676.496 200.598 317.349 150.546 566.615 100.003 03
    207.235 540.000 056.614 210.001 597.778 060.042 946.875 020.133 73
    217.778 050.071 156.875 130.207 347.984 620.026 007.121 850.000 71
    227.984 600.042 407.117 110.000 878.112 550.020 307.183 560.005 58
    238.112 540.033 217.171 630.005 938.241 690.004 937.186 370.070 85
    248.241 670.008 197.186 590.070 518.355 480.012 787.228 550.356 06
    258.355 470.020 327.218 460.014 678.603 980.003 847.349 130.541 36
    268.457 000.005 007.235 540.000 048.879 180.000 607.778 060.042 54
    278.603 970.006 277.778 050.071 158.980 480.011 407.984 620.025 76
    288.977 700.014 417.984 600.042 399.092 810.031 828.112 550.020 11
    299.092 800.051 768.112 540.033 219.161 850.066 658.241 690.004 88
    309.229 790.118 338.241 670.008 199.229 870.071 678.355 480.012 67
    319.414 480.015 028.355 470.020 329.414 480.009 218.603 980.003 80
    329.518 820.077 798.457 000.005 049.450 980.026 218.879 180.000 58
    339.698 920.086 258.603 970.006 279.518 870.045 528.980 480.011 29
    349.776 220.456 368.977 700.014 429.787 150.014 039.092 810.031 53
    359.841 210.024 009.092 800.051 769.841 170.020 969.161 810.065 98
    下载: 导出CSV 
    | 显示表格

    3条地震动输入下反应堆厂房位移响应规律基本一致,因此仅给出RG160地震动输下的反应堆厂房监测点位移响应,如图12所示。由图12可知,Case1和Case2模型对应的反应堆厂房位移响应较接近,Case3和Case4模型对应的反应堆厂房位移响应较接近,且响应较大。对于反应堆底板点f和上部点e,4种工况下的位移基本一致;对于外壳中部点c、内壳中部点d及外壳顶部点a、内壳顶部点b,Case3和Case4模型对应的反应堆厂房位移响应较Case1和Case2模型大,且顶部点位移差异较中部点更明显。综上所述,燃料厂房和电器厂房与外壳连接时,对反应堆厂房位移响应有影响,对反应堆底部和中部位移的影响较小,但对壳顶部位移的影响较大;安全厂房与反应堆厂房分离时,对反应堆厂房位移响应的影响较小。

    图 12  反应堆厂房监测点位移时程曲线
    Figure 12.  Displacement time history of each monitoring point on the reactor plant

    RG160地震动输入下4种工况反应堆厂房监测点加速度反应谱如图13所示。由图13可知,Case1和Case2模型对应的反应堆厂房加速度反应谱较接近,说明安全厂房对反应堆厂房加速度反应谱的影响较小;Case3和Case4模型对应的反应堆厂房加速度反应谱基本一致,且与Case1有差异,说明燃料厂房和电器厂房对反应堆厂房加速度反应谱有影响。对于反应堆底板点f和上部点e,4种工况下的反应堆厂房加速度反应谱基本一致,但安全壳(测点ad)加速度反应谱有所差异,即燃料厂房和电器厂房对反应堆底部和中部加速度反应谱的影响较小,对反应堆厂房安全壳加速度反应谱有明显影响。当频率>10 Hz时,4种工况下测点ad加速度反应谱基本一致;当频率为6~10 Hz时,Case3和Case4模型对应的测点ad加速度响应大于Case1和Case2模型;当频率为2~6 Hz时,Case3和Case4模型对应的加速度响应小于Case1和Case2模型;当频率<2 Hz时,Case3和Case4模型对应的加速度响应大于Case1和Case2模型,即燃料厂房和电器厂房的存在增大了反应堆厂房2 Hz以下的低频响应,因此增大了位移响应,如图13所示。

    图 13  反应堆厂房监测点加速度反应谱
    Figure 13.  Acceleration response spectrum of monitoring point on the reactor plant

    RSN5806和RSN4872地震动输入下4种工况反应堆厂房监测点加速度反应谱如图14所示,图中实线为RSN5806地震动输入下反应堆厂房加速度反应谱,虚线为RSN4872地震动输入下反应堆厂房加速度反应谱。由图14可知,在RSN5806和RSN4872地震动输入下,各工况对应的反应堆厂房加速度响应规律与RG160地震动输入时基本一致,即Case1和Case2模型对应的反应堆厂房加速度反应谱较接近,Case3和Case4模型对应的加速度反应谱基本一致,4中工况下应堆底板点f和上部点e加速度反应谱基本一致。

    图 14  反应堆厂房监测点加速度反应谱
    Figure 14.  Acceleration response spectrum of monitoring point on the reactor plant

    本文基于PASSI算法,提出显-隐式单元层计算方法,实现显、隐式交替计算。通过场地分析算例和土-结构相互作用分析算例,与ABAQUS软件计算结果进行对比,验证计算方法的可行性。以某核电站为研究对象,分析了相邻厂房对反应堆厂房地震响应的影响,得到以下主要结论:

    (1)显-隐式单元层计算方法较PASSI算法稳定性好,但增加了重叠区计算量。

    (2)安全厂房与反应堆厂房共用基础,但不连接,结构整体模态和自振频率为独立厂房模态和自振频率的集合。在基岩场地上,安全厂房通过基础和场地对反应堆厂房的作用较小,对反应堆厂房地震响应的影响较小。

    (3)燃料厂房、电器厂房与反应堆厂房共用基础,并与反应堆厂房相连,增加了反应堆厂房刚度,作为整体结构时,与反应堆厂房相关的自振频率大于单独的反应堆厂房。在本文地震动输入下,燃料厂房和电器厂房加大了反应堆厂房位移峰值,减小了反应堆厂房加速度反应谱峰值,并使反应堆厂房顶部点加速度反应谱峰值向高频移动。

    (4)对于与反应堆厂房共用同一基础的辅助厂房,应与反应堆厂房作为整体进行地震响应分析,至少应将与反应堆厂房相连的辅助厂房作为整体进行分析。

    对于非基岩场地,结构-土-结构相互作用可能更明显,其对反应堆厂房地震响应的影响与输入地震动特性等有关,需另外讨论。

  • 图  1  结构首层平面布置及外观

    Figure  1.  Floor plan and structural appearance

    图  2  结构数值模型及真实震害对比图

    Figure  2.  The damage comparison between numerical structural model and real earthquake

    图  3  地震动加速度反应谱、平均反应谱及规范设计谱

    Figure  3.  Seismic acceleration response spectrum, average response spectrum and code design spectrum

    图  4  基于DOE的结构地震易损性分析方法步骤

    Figure  4.  Flowchart of structural seismic vulnerability analysis method based on DOE

    图  5  标准化效应

    Figure  5.  Standardized effect diagram

    图  6  易损性曲线

    Figure  6.  Comparison of vulnerability curves

    表  1  数值模型信息

    Table  1.   Numerical model information

    建筑构件 单元类型 构件之间的相互作用 材料属性
    悬挑梁 C3D8R Tie 设计规范混凝土本构模型
    现浇楼板 C3D8R Tie 设计规范混凝土本构模型
    钢筋 T3D2 Embedded region 双直线理想弹塑性模型
    墙体 C3D8R Tie 整体式砌体结构本构模型
    下载: 导出CSV

    表  2  模型前2阶自振周期

    Table  2.   The first 2 vibration cycles of the model

    振型自振周期/s振型特征
    一阶0.141横向平动
    二阶0.115纵向平动
    下载: 导出CSV

    表  3  结构参数及其分布

    Table  3.   Structural model parameters and their distribution

    参数 平均值 变异系数 分布类型
    砌体材料参数 阻尼比 0.05 0.3 正态
    弹性模量/MPa 1 350 0.15 正态
    抗压强度/kPa 2 448.29 0.17 对数正态
    抗拉强度/kPa 235.938 0.2 对数正态
    密度/(kg·m−3) 1 600 0.1 正态
    混凝土材料参数 密度/(kg·m−3) 2 400 0.07 正态
    结构几何参数 层高/m 3.3 0.05
    墙厚/mm 240 0.05
    下载: 导出CSV

    表  4  输入参数及标准化

    Table  4.   Input parameters and standardization

    输入参数 下限值 平均值 上限值
    砌体阻尼比 输入值 0.035 0.05 0.065
    标准化 −1 0 1
    砌体弹性模量/MPa 输入值 1 147.5 1 350.0 1 552.5
    标准化 −1 0 1
    砌体抗压强度/kPa 输入值 2 032.08 2 448.29 2 864.50
    标准化 −1 0 1
    砌体抗拉强度/kPa 输入值 188.750 235.938 283.126
    标准化 −1 0 1
    砌体材料密度/(kg·m−3) 输入值 1 440 1 600 1 760
    标准化 −1 0 1
    混凝土密度/(kg·m−3) 输入值 2 328 2 400 2 568
    标准化 −1 0 1
    层高/m 输入值 3.135 3.300 3.465
    标准化 −1 0 1
    墙厚/mm 输入值 228 240 252
    标准化 −1 0 1
    下载: 导出CSV

    表  5  Plackett-Burman试验设计

    Table  5.   A case of Plackett-Burman test design

    项目 序号
    1 2 3 4 5 6 7 8 9 10 11 12
    砌体阻尼比 −1 1 −1 −1 1 1 −1 1 −1 −1 1 1
    砌体密度 −1 −1 1 −1 −1 1 −1 1 1 1 1 −1
    混凝土密度 −1 −1 1 1 1 1 −1 −1 1 −1 −1 1
    砌体抗拉强度 −1 1 1 1 −1 1 1 −1 −1 −1 1 −1
    砌体弹性模量 −1 1 −1 −1 1 1 1 −1 1 1 −1 −1
    墙厚 −1 1 1 −1 −1 −1 1 1 1 −1 −1 1
    层高 −1 −1 1 −1 1 −1 1 −1 −1 1 1 1
    砌体抗压强度 −1 −1 −1 1 −1 1 1 1 −1 1 −1 1
    下载: 导出CSV

    表  6  中心复合试验设计

    Table  6.   A case of central composite test design

    运行序砌体弹性模量砌体密度砌体抗拉强度
    111−1
    2−1−1−1
    3100
    4010
    51−1−1
    6−11−1
    7000
    8111
    9−1−11
    10−111
    11−100
    121−11
    130−10
    14001
    1500−1
    下载: 导出CSV

    表  7  中心复合设计输出结果

    Table  7.   Output results from the central composite design

    运行序最大层间位移角(×10−7
    0.05 g0.1 g0.15 g...0.35 g0.4 g
    $ {\hat y_{u1}} $$ {\hat y_{\sigma 1}} $$ {\hat y_{u2}} $$ {\hat y_{u2}} $$ {\hat y_{u3}} $$ {\hat y_{\sigma 3}} $......$ {\hat y_{u7}} $$ {\hat y_{\sigma 7}} $$ {\hat y_{u8}} $$ {\hat y_{\sigma 8}} $
    13 949.92 209.48 740.63 417.635 78010 027......29 8327 952.235 7791 002
    22 998.81 969.577 020.42 665.228 1836 577......24 3166 219.328 1826 577.2
    34 499.31 736.19 200.13733.539 73912 563......33 4419 797.439 73912 563
    .......................................
    154 235.61 558.18 759.32 412.739 00312 186......32 4729 523.839 00312 186
    下载: 导出CSV

    表  8  响应面模型检验结果

    Table  8.   Verification of response surface model

    响应面标号 $ {R^2} $/% $ R_{\mathrm{A}}^2 $/%
    1 99.59 98.31
    2 98.61 97.52
    3 98.36 97.62
    4 98.76 97.32
    5 97.79 96.80
    6 98.26 97.82
    7 97.12 96.32
    8 95.30 93.25
    下载: 导出CSV

    表  9  极限状态限值

    Table  9.   Limitation value of limited state

    极限状态 轻微破坏(LS1 中等破坏(LS2 严重破坏(LS3 毁坏(LS4
    $ {\theta _{\max }} $ 1/2 000 1/1 600 1/700 1/350
    下载: 导出CSV
  • 陈力波,黄才贵,谷音,2018. 基于改进响应面法的公路简支梁桥地震易损性分析. 工程力学,35(4):208−218. doi: 10.6052/j.issn.1000-4750.2017.06.0512

    Chen L. B., Huang C. G., Gu Y., 2018. Seismic vulnerability analysis of simply supported highway bridges based on an improved response surface method. Engineering Mechanics, 35(4): 208−218. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.06.0512
    黄琼,2014. 基于响应面法的钢筋混凝土框架结构有限元模型修正研究. 兰州:兰州理工大学.

    Huang Q., 2014. Finite element model updating of reinforced concrete frame structure based on response surface method. Lanzhou:Lanzhou University of Technology. (in Chinese)
    李龙师,2020. 松原查干花震区农村房屋抗震性能研究及震害预测. 哈尔滨:中国地震局工程力学研究所.

    Li L. S., 2020. Seismic Capacity Analysis and Earthquake damage prediction of rural houses in Chaganhua earthquake area of Songyuan City. Harbin:Institute of Engineering Mechanics,China Earthquake Administration. (in Chinese)
    吕大刚,刘洋,于晓辉,2019. 第二代基于性能地震工程中的地震易损性模型及正逆概率风险分析. 工程力学,36(9):1−11,24.

    Lü D. G., Liu Y., Yu X. H., 2019. Seismic fragility models and forward-backward probabilistic risk analysis in second-generation performance-based earthquake engineering. Engineering Mechanics, 36(9): 1−11,24. (in Chinese)
    宋帅,钱永久,钱聪,2018. 桥梁地震需求中随机参数的重要性分析方法研究. 工程力学,35(3):106−114. doi: 10.6052/j.issn.1000-4750.2016.11.0856

    Song S., Qian Y. J., Qian C., 2018. Research on methods for importance analysis of random parameter in bridge seismic demand. Engineering Mechanics, 35(3): 106−114. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.11.0856
    孙柏涛,张桂欣,2012. 汶川8.0级地震中各类建筑结构地震易损性统计分析. 土木工程学报,45(5):26−30.

    Sun B. T., Zhang G. X., 2012. Statistical analysis of the seismic vulnerability of various types of building structures in Wenchuan M 8.0 earthquake. China Civil Engineering Journal, 45(5): 26−30. (in Chinese)
    熊立红,吴文博,孙悦,2012. 汶川地震作用下约束砌体房屋的抗震能力分析. 土木工程学报,45(S2):103−108.

    Xiong L. H., Wu W. B., Sun Y., 2012. Seismic performance of confined masonry buildings during the Wenchuan earthquake. China Civil Engineering Journal, 45(S2): 103−108. (in Chinese)
    杨卫忠,2008. 砌体受压本构关系模型. 建筑结构,38(10):80−82.

    Yang W. Z., 2008. Constitutive relationship model for masonry materials in compression. Building Structure, 38(10): 80−82. (in Chinese)
    杨玉成,杨柳,高云学等,1982. 现有多层砖房震害预测的方法及其可靠度. 地震工程与工程振动,2(3):75−86.

    Yang Y. C., Yang L., Gao Y. X., et al., 1982. Method of damage prediction for existing multi story brick buildings and its reliability. Earthquake Engineering and Engineering Vibration, 2(3): 75−86. (in Chinese)
    叶继红,江力强,2021. 考虑多重不确定性的我国多层冷成型钢结构地震风险评估. 土木工程学报,54(2):74−83,126.

    Ye J. H., Jiang L. Q., 2021. Seismic risk assessment of mid-rise cold-formed steel structures in China considering various uncertainties. China Civil Engineering Journal, 54(2): 74−83,126. (in Chinese)
    于晓辉,2012. 钢筋混凝土框架结构的概率地震易损性与风险分析. 哈尔滨:哈尔滨工业大学.

    Yu X. H., 2012. Probabilistic seismic fragility and risk analysis of reinforced concrete frame structures. Harbin:Harbin Institute of Technology. (in Chinese)
    于晓辉,吕大刚,2012. 考虑结构不确定性的地震倒塌易损性分析. 建筑结构学报,33(10):8−14.

    Yu X. H., Lü D. G., 2012. Seismic collapse fragility analysis considering structural uncertainties. Journal of Building Structures, 33(10): 8−14. (in Chinese)
    张令心,李孟达,刘洁平等,2018. 考虑结构不确定性的框剪结构地震易损性分析. 自然灾害学报,27(4):112−118.

    Zhang L. X., Li M. D., Liu J. P., et al., 2018. Seismic fragility analysis of frame-shear wall structures considering structural parameter uncertainty. Journal of Natural Disasters, 27(4): 112−118. (in Chinese)
    张尚荣,谭平,杜永峰等,2014. 基于响应面法的层间隔震结构地震易损性分析. 振动与冲击,33(15):42−48.

    Zhang S. R., Tan P., Du Y. F., et al., 2014. Seismic fragility analysis of inter-story isolation structures based on response surface method. Journal of Vibration and Shock, 33(15): 42−48. (in Chinese)
    周强,2012. 砌体结构抗震试验及弹塑性地震反应分析. 哈尔滨:哈尔滨工程大学.

    Zhou Q., 2012. Seismic test and elastic-plastic seismic response analysis of masonry structure. Harbin:Harbin Engineering University. (in Chinese)
    周强,闵全环,熊拥军等,2018. 江西典型村镇既有房屋抗震性能分析. 自然灾害学报,27(1):96−105.

    Zhou Q., Min Q. H., Xiong Y. J., et al., 2018. Seismic performance analysis of existing typical rural buildings in Jiangxi. Journal of Natural Disasters, 27(1): 96−105. (in Chinese)
    朱震宇,2015. 砌体墙抗震性能分析模型及其工程应用. 上海:上海交通大学.

    Zhu Z. Y., 2015. Analysis model of seismic behavior of masonry walls and its engineering application. Shanghai:Shanghai Jiao Tong University. (in Chinese)
    Choudhury T., Kaushik H. B., 2018. Seismic response sensitivity to uncertain variables in RC frames with infill walls. Journal of Structural Engineering, 144(10): 04018184. doi: 10.1061/(ASCE)ST.1943-541X.0002190
    Demirel İ. O., 2010. A nonlinear equivalent frame model for displacement based analysis of unreinforced brick masonry buildings. Ankara: Middle East Technical University.
    FEMA, 2009. Quantification of building seismic performance factors. Washington: Federal Emergency Management Agency.
    Freeny A., 1988. Empirical model building and response surfaces. Technometrics, 30(2): 229−231. doi: 10.1080/00401706.1988.10488371
    Gokkaya B. U., Baker J. W., Deierlein G. G., 2016. Quantifying the impacts of modeling uncertainties on the seismic drift demands and collapse risk of buildings with implications on seismic design checks. Earthquake Engineering & Structural Dynamics, 45(10): 1661−1683.
    Loukas Y. L., 2001. A Plackett-Burman screening design directs the efficient formulation of multicomponent DRV liposomes. Journal of Pharmaceutical and Biomedical Analysis, 26(2): 255−263. doi: 10.1016/S0731-7085(01)00419-8
    Montgomery D. C., 2001. Design and analysis of experiments. 5th ed. New York: John Wiley & Sons, Inc.
    Parisi F., Augenti N., 2012. Uncertainty in seismic capacity of masonry buildings. Buildings, 2(3): 218−230. doi: 10.3390/buildings2030218
    Towashiraporn P., 2004. Building seismic fragilities using response surface metamodels. Atlanta: Georgia Institute of Technology.
  • 期刊类型引用(0)

    其他类型引用(1)

  • 加载中
图(6) / 表(9)
计量
  • 文章访问数:  122
  • HTML全文浏览量:  14
  • PDF下载量:  23
  • 被引次数: 1
出版历程
  • 收稿日期:  2023-09-06
  • 网络出版日期:  2024-10-15
  • 刊出日期:  2024-09-01

目录

/

返回文章
返回