CN112906103B - 主余震作用下rc框架结构抗震评估方法 - Google Patents

主余震作用下rc框架结构抗震评估方法 Download PDF

Info

Publication number
CN112906103B
CN112906103B CN202110069550.9A CN202110069550A CN112906103B CN 112906103 B CN112906103 B CN 112906103B CN 202110069550 A CN202110069550 A CN 202110069550A CN 112906103 B CN112906103 B CN 112906103B
Authority
CN
China
Prior art keywords
energy
model
aftershock
frame structure
main
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110069550.9A
Other languages
English (en)
Other versions
CN112906103A (zh
Inventor
籍多发
翟长海
温卫平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202110069550.9A priority Critical patent/CN112906103B/zh
Publication of CN112906103A publication Critical patent/CN112906103A/zh
Application granted granted Critical
Publication of CN112906103B publication Critical patent/CN112906103B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)

Abstract

主余震作用下RC框架结构抗震评估方法,本发明涉及一种RC框架结构的抗震评估方法,它要解决现有提出的抗震评估方法均是基于单次地震,无法考虑余震作用的问题。抗震评估方法:一、挑选主余震地震动记录;二、建立IMK模型的能量平衡公式,推导能量系数公式;三、调整IMK模型的动力特性,根据步骤一中的主余震地震动记录和步骤二中能量系数计算公式,求解主余震作用下IMK模型的能量系数γ,拟合得到能量系数模型;步骤四、根据能量系数模型计算能量需求谱,计算RC框架结构基于能量的能力谱,确定能力谱和能量需求谱的交点,获得地震响应。本发明与非线性动力分析法的结果吻合程度较高,验证了该方法的正确性。

Description

主余震作用下RC框架结构抗震评估方法
技术领域
本发明涉及一种RC框架结构的抗震评估方法。
背景技术
震害调查资料表明,强余震会进一步加重结构的震害。余震不仅会造成严重的经济损失和人员伤亡,还会导致震后恢复时间延长。例如,2010年新西兰7.4级主震后的一年半内,共发生了上万次余震,严重阻碍了震后修复重建工作的进度,修复工作至今仍未完成。
我国乃至世界上的抗震规范均仅仅考虑主震作用而忽略了主余震对结构造成的累积损伤,因此,关于如何评估这类结构在主余震作用下的抗震性能,目前还缺少科学系统的基础研究,历史地震的教训也昭示了对RC框架结构进行抗震评估的重要性和迫切性。
发明内容
本发明的目的是要解决现有提出的抗震评估方法均是基于单次地震,无法考虑余震作用的问题,而提出的一种基于能量的主余震的抗震评估方法。
本发明主余震作用下RC框架结构抗震评估方法按照以下步骤实现:
步骤一、以场地和余震与主震强度比值为条件,挑选主余震地震动记录;
步骤二、建立IMK模型的能量平衡公式,推导能量系数公式如下:
Figure RE-GDA0002987068440000011
其中:γ代表能量系数,μ代表延性系数,xu代表极限位移,xy代表屈服位移,xc代表峰值位移,Sv代表拟速度,m代表质量,α代表屈服刚度,Ry代表强度折减系数,αh代表软化刚度,Ee代表弹性能量,Ep代表塑性能量;
步骤三、调整IMK模型的动力特性,根据步骤一中的主余震地震动记录和步骤二中能量系数计算公式,求解主余震作用下IMK模型的能量系数γ,拟合得到能量系数模型为:
Figure RE-GDA0002987068440000021
其中:a=a1+a2ξ+a3ξ2,c=c1+c2ξ+c3ξ2,T为IMK模型的周期,a1-a3、b1-b8、c1-c3均为拟合系数;
步骤四、根据步骤三中能量系数模型计算能量需求谱Ea,计算公式如下:
Figure RE-GDA0002987068440000022
计算RC框架结构基于能量的能力谱:
Ek=ΔEk
Figure RE-GDA0002987068440000023
Figure RE-GDA0002987068440000024
其中:Ek为能力谱,ΔEk是第k步RC框架结构的能量增量,ΔEk,i是第k步第i层的能量增量,Fk,i是第k步第i层的水平力,ΔFk,i是第k步第i层的水平力增量,Δδk,i是第k步第i层的水平位移增量;
确定能力谱和能量需求谱的交点,获得RC框架结构的地震响应,通过地震响应实现余震作用下RC框架结构的抗震评估。
本发明利用已记录的壳内和俯冲区域主余震地震动记录,选用能够表征RC框架结构的单自由度模型-IMK模型为例,建立基于能量平衡的抗震评估方法。本发明所述的主余震作用下RC框架结构抗震评估方法主要包括四个部分:第一、对现有主余震地震动进行收集、分类(壳内和俯冲区域),分别计算余震相对于主震的强度。第二、以IMK模型为例,建立能量平衡方程,并推导出能量系数公式。第三、进行参数分析,提出主余震作用IMK 模型的能量系数预测模型。第四、利用能量系数计算出需求谱,并找到与能力谱的交点,最终获得主余震作用下RC框架结构的地震响应。
附图说明
图1为实施例中建筑的平面图,图中连接处为抗弯连接;
图2为实施例中建筑的立面图,图中连接处为单铰连接;
图3为实施例中建筑结构在振型1下的三阶弹性模态图;
图4为实施例中建筑结构在振型2下的三阶弹性模态图;
图5为实施例中建筑结构在振型3下的三阶弹性模态图;
图6为实施例中建筑结构在振型1下的等效单自由度力-位移关系图;
图7为实施例中建筑结构在振型2下的等效单自由度力-位移关系图;
图8为实施例中建筑结构在振型3下的等效单自由度力-位移关系图;
图9为实施例中壳内主余震作用下ΙΜK模型的能量系数图,沿着箭头方向μ的取值依次为2、3、4、5和6;
图10为实施例中俯冲区域主余震作用下ΙΜK模型的能量系数图,沿着箭头方向μ的取值依次为2、3、4、5和6。
图11为实施例中
Figure RE-GDA0002987068440000031
条件下本发明抗震评估方法和非线性动力分析计算的楼层位移对比图,其中□代表非线性动力,○代表实施例中抗震评估方法;
图12为实施例中
Figure RE-GDA0002987068440000032
条件下本发明抗震评估方法和非线性动力分析计算的层间位移角对比图,其中□代表非线性动力,○代表实施例中抗震评估方法;
图13为实施例中
Figure RE-GDA0002987068440000033
条件下本发明抗震评估方法和非线性动力分析计算的楼层位移对比图,其中□代表非线性动力,○代表实施例中抗震评估方法;
图14为实施例中
Figure RE-GDA0002987068440000034
条件下本发明抗震评估方法和非线性动力分析计算的层间位移角对比图,其中□代表非线性动力,○代表实施例中抗震评估方法。
具体实施方式
具体实施方式一:本实施方式主余震作用下RC框架结构抗震评估方法按照以下步骤实现:
步骤一、以场地和余震与主震强度比值为条件,挑选主余震地震动记录;
步骤二、建立IMK模型的能量平衡公式,推导能量系数公式如下:
Figure RE-GDA0002987068440000035
其中:γ代表能量系数,μ代表延性系数,xu代表极限位移,xy代表屈服位移,xc代表峰值位移,Sv代表拟速度,m代表质量,α代表屈服刚度,Ry代表强度折减系数,αh代表软化刚度,Ee代表弹性能量,Ep代表塑性能量;
步骤三、调整IMK模型的动力特性,根据步骤一中的主余震地震动记录和步骤二中能量系数计算公式,求解主余震作用下IMK模型的能量系数γ,拟合得到能量系数模型为:
Figure RE-GDA0002987068440000041
其中:a=a1+a2ξ+a3ξ2,c=c1+c2ξ+c3ξ2,T为IMK模型的周期,a1-a3、b1-b8、c1-c3均为拟合系数;
步骤四、根据步骤三中能量系数模型计算能量需求谱Ea,计算公式如下:
Figure RE-GDA0002987068440000042
计算RC框架结构基于能量的能力谱:
Ek=ΔEk
Figure RE-GDA0002987068440000043
Figure RE-GDA0002987068440000044
其中:Ek为能力谱,ΔEk是第k步RC框架结构的能量增量,ΔEk,i是第k步第i层的能量增量,Fk,i是第k步第i层的水平力,ΔFk,i是第k步第i层的水平力增量,Δδk,i是第k步第i层的水平位移增量;
确定能力谱和能量需求谱的交点,获得RC框架结构的地震响应,通过地震响应实现余震作用下RC框架结构的抗震评估。
具体实施方式二:本实施方式与具体实施方式一不同的是步骤一中选择壳内区域地震和俯冲区域地震。
具体实施方式三:本实施方式与具体实施方式一或二不同的是步骤一中余震与主震强度比值选用
Figure RE-GDA0002987068440000045
来表示,定义如下:
Figure RE-GDA0002987068440000046
其中PGAas为余震的峰值加速度,PGAms为主震的峰值加速度。
具体实施方式四:本实施方式与具体实施方式三不同的是
Figure RE-GDA0002987068440000047
的值选取0.1-0.5、0.5-0.8、0.8-1.0和1.0-2.2。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是步骤三中所述的 IMK模型的动力特性包括IMK模型的周期范围、周期间隔和侧向最大位移。
具体实施方式六:本实施方式与具体实施方式五不同的是IMK模型的周期范围为0.1-3.0s,周期间隔为0.1s。
具体实施方式七:本实施方式与具体实施方式五不同的是IMK模型的侧向最大位移根据延性系数μ确定,μ的取值为2、3、4、5、6。
本实施方式延性系数μ的计算公式为:
Figure RE-GDA0002987068440000051
其中所述xmax为主余震作用下IMK模型的最大位移反应,xy为IMK模型的屈服位移。
实施例:本实施例主余震作用下RC框架结构抗震评估方法按照以下步骤实现:
步骤一、本实施例RC框架结构建筑物的结构图如图1和图2所示,从余震和主震相对强度和场地条件的角度来选取有代表性的主余震地震动,选择壳内地震和俯冲区域地震,选用
Figure RE-GDA0002987068440000052
来表示余震和主震相对强度,定义如下:
Figure RE-GDA0002987068440000053
其中PGAas为余震的峰值加速度,PGAms为主震的峰值加速度,将主余震分为四类,分别对应于
Figure RE-GDA0002987068440000054
在0.1-0.5、0.5-0.8、0.8-1.0和1.0-2.2,根据场地条件,分为两类,场地C和场地D。具体分布见下表1所示;
表1主余震地震动分布
Figure RE-GDA0002987068440000055
步骤二、从能量平衡原理出发,建立IMK模型(图9和图10所示)的能量平衡公式,并给出其适用范围,推导能量系数公式如下:
Figure RE-GDA0002987068440000061
其中:γ代表能量系数,μ代表延性系数,xu代表极限位移,xy代表屈服位移,xc代表峰值位移,Sv代表拟速度,m代表质量,α代表屈服刚度,Ry代表强度折减系数,αh代表软化刚度,Ee代表弹性能量,Ep代表塑性能量;
步骤三、调整IMK模型的动力特性,IMK模型的周期范围为0.1-3.0s,周期间隔为0.1s, IMK模型的侧向最大位移根据延性系数μ确定,根据步骤一中的主余震地震动记录和步骤二中能量系数计算公式,求解主余震作用下IMK模型的能量系数γ,拟合得到能量系数模型为:
Figure RE-GDA0002987068440000062
其中:a=a1+a2ξ+a3ξ2,c=c1+c2ξ+c3ξ2,T为IMK模型的周期,a1-a3、b1-b8、c1-c3为拟合系数,ξ代表阻尼比。见表2所示;
表2拟合系数取值
Figure RE-GDA0002987068440000063
Figure RE-GDA0002987068440000071
步骤四、根据步骤三中能量系数模型计算能量需求谱Ea,计算公式如下:
Figure RE-GDA0002987068440000072
计算RC框架结构基于能量的能力谱:
Ek=ΔEk
Figure RE-GDA0002987068440000073
Figure RE-GDA0002987068440000074
其中:Ek为能力谱,ΔEk是第k步RC框架结构的能量增量,ΔEk,i是第k步第i层的能量增量,Fk,i是第k步第i层的水平力,ΔFk,i是第k步第i层的水平力增量,Δδk,i是第k步第i层的水平位移增量;其中Δδk,i由Pushover程序设定,pushover加载时控制每步位移,当达到加载位移时,此时所用力为Fk,i,两步所用力的差值为ΔFk,i
确定能力谱和能量需求谱的交点,从而获得RC框架结构的地震响应,通过地震响应实现余震作用下RC框架结构的抗震评估。
本实施例选取了某6层建筑为例,来演示本方法评估地震作用的准确性。该结构的前三阶模态如图3-5所示,其等效单自由度如图6-8所示。图11-图14给出了主余震作用下该结构的楼层位移、层间位移角分布情况,从该图可知,随着余震地震动的相对强度增加,该楼层的位移和层间位移角都随之增加了,同时,本发明所提方法与传统非线性动力分析法的结果吻合程度较高,验证了该方法的正确性。

Claims (6)

1.主余震作用下RC框架结构抗震评估方法,其特征在于该抗震评估方法按照以下步骤实现:
步骤一、以场地和余震与主震强度比值为条件,选择壳内区域地震和俯冲区域地震,挑选主余震地震动记录;
步骤二、建立IMK模型的能量平衡公式,推导能量系数公式如下:
Figure FDA0003213839370000011
其中:γ代表能量系数,μ代表延性系数,xu代表极限位移,xy代表屈服位移,xc代表峰值位移,Sv代表拟速度,m代表质量,α代表屈服刚度,Ry代表强度折减系数,αh代表软化刚度,Ee代表弹性能量,Ep代表塑性能量;
步骤三、调整IMK模型的动力特性,根据步骤一中的主余震地震动记录和步骤二中能量系数计算公式,求解主余震作用下IMK模型的能量系数γ,拟合得到能量系数模型为:
Figure FDA0003213839370000012
其中:a=a1+a2ξ+a3ξ2,c=c1+c2ξ+c3ξ2,T为IMK模型的周期,a1-a3、b1-b8、c1-c3均为拟合系数,ξ代表阻尼比;
步骤四、根据步骤三中能量系数模型计算能量需求谱Ea,计算公式如下:
Figure FDA0003213839370000013
计算RC框架结构基于能量的能力谱:
Ek=ΔEk
Figure FDA0003213839370000014
Figure FDA0003213839370000021
其中:Ek为能力谱,ΔEk是第k步RC框架结构的能量增量,ΔEk,i是第k步第i层的能量增量,Fk,i是第k步第i层的水平力,ΔFk,i是第k步第i层的水平力增量,Δδk,i是第k步第i层的水平位移增量;
确定能力谱和能量需求谱的交点,获得RC框架结构的地震响应,通过地震响应实现余震作用下RC框架结构的抗震评估。
2.根据权利要求1所述的主余震作用下RC框架结构抗震评估方法,其特征在于步骤一中余震与主震强度比值选用
Figure FDA0003213839370000022
来表示,定义如下:
Figure FDA0003213839370000023
其中PGAas为余震的峰值加速度,PGAms为主震的峰值加速度。
3.根据权利要求2所述的主余震作用下RC框架结构抗震评估方法,其特征在于▽PGA的值选取0.1-0.5、0.5-0.8、0.8-1.0和1.0-2.2。
4.根据权利要求1所述的主余震作用下RC框架结构抗震评估方法,其特征在于步骤三中所述的IMK模型的动力特性包括IMK模型的周期范围、周期间隔和侧向最大位移。
5.根据权利要求4所述的主余震作用下RC框架结构抗震评估方法,其特征在于IMK模型的周期范围为0.1-3.0s,周期间隔为0.1s。
6.根据权利要求1所述的主余震作用下RC框架结构抗震评估方法,其特征在于IMK模型的侧向最大位移根据延性系数μ确定,μ的取值为2、3、4、5、6。
CN202110069550.9A 2021-01-19 2021-01-19 主余震作用下rc框架结构抗震评估方法 Active CN112906103B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110069550.9A CN112906103B (zh) 2021-01-19 2021-01-19 主余震作用下rc框架结构抗震评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110069550.9A CN112906103B (zh) 2021-01-19 2021-01-19 主余震作用下rc框架结构抗震评估方法

Publications (2)

Publication Number Publication Date
CN112906103A CN112906103A (zh) 2021-06-04
CN112906103B true CN112906103B (zh) 2021-10-15

Family

ID=76115785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110069550.9A Active CN112906103B (zh) 2021-01-19 2021-01-19 主余震作用下rc框架结构抗震评估方法

Country Status (1)

Country Link
CN (1) CN112906103B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105528520A (zh) * 2015-12-09 2016-04-27 哈尔滨工业大学 一种基于主余震序列的损伤指数评估方法
CN108256141A (zh) * 2017-12-11 2018-07-06 哈尔滨工业大学 一种基于Copula理论的主余震联合易损性分析方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105528520A (zh) * 2015-12-09 2016-04-27 哈尔滨工业大学 一种基于主余震序列的损伤指数评估方法
CN108256141A (zh) * 2017-12-11 2018-07-06 哈尔滨工业大学 一种基于Copula理论的主余震联合易损性分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Earthquake risk in reinforced concrete buildings during aftershock sequences based on period elongation and operational earthquake forecasting;Konstantinos Trevlopoulos et al.;《Structural Safety》;20200218;全文 *
主余震地震动特征及结构地震反应分析;籍多发;《中国博士学位论文全文数据库》;20190131(第01期);全文 *
基于改进IMK恢复力模型的填充墙RC框架结构地震易损性分析;张永兵 等;《重庆大学学报》;20180331;第41卷(第3期);全文 *

Also Published As

Publication number Publication date
CN112906103A (zh) 2021-06-04

Similar Documents

Publication Publication Date Title
Michel et al. Full‐scale dynamic response of an RC building under weak seismic motions using earthquake recordings, ambient vibrations and modelling
CN106709182B (zh) 一种地震作用下顺层岩质边坡稳定可靠性安全评价方法
CN108957528B (zh) 设定地震风险的性能抗震设计评估方法
WO2019148950A1 (zh) 城市建筑群地震反应非线性历程分析方法及装置
CN108446444B (zh) 基于性能水准的多模态性能化抗震设计方法
WO2019153798A1 (zh) 基于性态双谱的多模态性能化抗震设计方法
CN112906103B (zh) 主余震作用下rc框架结构抗震评估方法
CN108427833B (zh) 基于三维弹塑性地震反应谱的性能化抗震设计方法
CN114792020A (zh) 一种基于机器学习的建筑抗震韧性快速评估方法及系统
Seyedkazemi et al. Comparison of static pushover analysis and IDA-based probabilistic methods for assessing the seismic performance factors of diagrid structures
Abbasnia et al. An improved displacement-based adaptive pushover procedure based on factor modal combination rule
Sabouniaghdam et al. Probabilistic collapse assessment of steel frame structures considering the effects of soil-structure interaction and height
Zhou et al. Performance evaluation of CWH base isolated building during two major earthquakes in Christchurch
CN102706525B (zh) 一种lng 储罐隔震效应的确定方法
CN101967860A (zh) 最不利设计地震动记录的确定方法
Samanta et al. Duration-dependent inelastic response spectra and effect of ground motion duration
CN110727026B (zh) 一种基于地震-海啸连续作用的非弹性位移比谱模型
CN110990910B (zh) 时程激励下线性耗能结构响应的快速迭代法
Kalkan et al. Lateral load distribution in nonlinear static procedures for seismic design
Sozen Lateral drift of reinforced concrete structures subjected to strong ground motion
Masaeli et al. Effect of Soil and Structure Nonlinear Interaction on the Efficiency of Tuned Mass Damper
Fallah et al. Accuracy evaluation of the modal pushover analysis method in the prediction of seismic response of vertically irregular frames
Chase et al. DAMAGE ACCUMULATION IN A TWO-STORY WOOD-FRAME BUILDING IN SEQUENCES OF INDUCED EARTHQUAKES
Mohammadi et al. Alternative Method to Estimate IDA Curves
Ramirez et al. Simplified analysis for preliminary design of base-isolated structures

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant