CN116595821A - 一种考虑建模不确定性的结构地震易损性分析方法 - Google Patents

一种考虑建模不确定性的结构地震易损性分析方法 Download PDF

Info

Publication number
CN116595821A
CN116595821A CN202310413818.5A CN202310413818A CN116595821A CN 116595821 A CN116595821 A CN 116595821A CN 202310413818 A CN202310413818 A CN 202310413818A CN 116595821 A CN116595821 A CN 116595821A
Authority
CN
China
Prior art keywords
modeling
vulnerability
parameters
uncertainty
earthquake
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.)
Pending
Application number
CN202310413818.5A
Other languages
English (en)
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.)
CITIC General Institute of Architectural Design and Research Co Ltd
Original Assignee
CITIC General Institute of Architectural Design and Research Co Ltd
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 CITIC General Institute of Architectural Design and Research Co Ltd filed Critical CITIC General Institute of Architectural Design and Research Co Ltd
Priority to CN202310413818.5A priority Critical patent/CN116595821A/zh
Publication of CN116595821A publication Critical patent/CN116595821A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提出一种考虑建模不确定性的结构地震易损性分析方法,包括不确定性建模参数的选取,精细化有限元模型的建立,基于耐震时程法的地震动输入,基于云图法的地易损性分析,基于龙卷风法的敏感性分析,基于中心点法的不确定性分析。本发明能够对建模随机参数的敏感性进行排序,量化建模不确定性对地震易损性分析的影响。识别出的影响相对较小或不明显的建模参数有助于减少这些参数的计算工作量,确定的关键建模参数对于实现改进的地震易损性评估是很关键的,降低关键建模参数的不确定性可以提高结构的可靠性估计。

Description

一种考虑建模不确定性的结构地震易损性分析方法
技术领域
本发明涉及工程结构抗震性能评估技术领域,具体涉及一种考虑建模不确定性的结构地震易损性分析方法。
背景技术
基于性能抗震设计的出现使得地震工程的研究和应用越来越倾向于根据一组性能化目标进行结构设计。这意味着工程师必须确保结构在整个使用期间满足足够的安全水平。结构的安全性一般通过可靠性分析进行量化,主要目的是考虑所有不确定性的来源及其对结构性能的影响,从而计算失效概率。
地震易损性分析是评估地震作用下结构可靠性的有力工具,被广泛应用于预测结构在地震作用下的破坏程度。然而,由于模型本身及其参数都是不确定的,结构的失效概率与易损性分析中的初步假设密切相关,模型输入参数的微小变化可能导致地震易损性估计的相关可变性。
如果不知道建模参数对传播各种不确定性的敏感性,就很难知道易损性评估的水平。如果忽略重要参数,这可能导致不可靠的易损性估计;相反地如果所有随机参数的不确定性都考虑,则可能会在统计和模拟上花费不必要的精力。在设计和优化过程中提高对建模参数的认识,能够使工程师对目标性能水平做出更好的决策,更关键的结构就会承受更少的损坏和最小的功能损失,结构所承受的危险就更加可控制和预测。
发明内容
本发明的目标是开发有效的方法处理结构地震易损性中建模方面的不确定性,指导未来的不确定性处理、计算资源的管理以及通过现场测试等方式对精确建模参数的估计进行投资,进而提供一种考虑建模不确定性的结构地震易损性分析方法。
本发明提供的技术方案是:一种考虑建模不确定性的结构地震易损性分析方法,包括以下若干步骤:
第一步:不确定性建模参数的选取。
通过对既有研究成果的整理分析,确定结构建模随机参数,假定建模参数的随机分布特征,选择参数分析区间定出建模参数的均值、上界值、下界值。
第二步:精细化有限元模型的建立。
选取合适的材料本构模型、截面模型、单元模型,采用有限元软件建立精细化结构数值模型。根据第一步的建模参数的分析区间分别建立均值模型、上界值模型、下界值模型。
第三步:基于耐震时程法的地震动输入。
选择目标反应谱和目标时间点,根据在某一时间内反应谱的大小与持续时间成线性关系,求解下列方程合成耐震时程曲线。
式中:SaT(T,t)为t时刻的加速度反应谱;T为结构的自振周期;t为任意时间点;tTarget为目标时间点;SaC(T)为目标加速度反应谱;SuT(T,t)为目标位移反应谱;ag为寻求的耐震加时程曲线;Sa(T,t)和Su(T,t)分别为ag在t时刻的加速度反应谱和位移反应谱;α为位移谱的相对权重系数。
第四步:基于云图法的地易损性分析。
假设结构的需求和能力的概率特征由对数正态分布描述,结构需求参数的均值与地震动强度参数的关系满足幂指数关系,采用下列公式建立结构的地震易损性函数。
式中:Pf为结构达到或超过损伤极限状态的概率;Φ[]为标准正态分布函数;为结构需求参数的均值;βSd|IM为地震响应在给定IM范围内的条件对数标准差;/>βc分别为结构能力参数的均值和对数标准差;IM为地震动强度参数;a和b为未知的回归系数;Di为在第i条地震动作用下结构需求响应的峰值;N为进行非线性动力分析的次数。
将第三步得到的耐震时程曲线输入到第二步建立的数值模型中进行非线性时程分析得到结构的地震响应。定义结构的损伤极限状态,选取合适的地震动强度参数和结构需求参数,代入到地震易损性函数中得到易损性曲线。
第五步:基于龙卷风法的敏感性分析。
以结构达到或超过某种损伤极限状态的概率为50%时所对应的易损性中位值为摆幅指标,根据第四步得到的易损性曲线采用下式计算摆幅值,根据摆幅值的大小从高到低进行排序,建立模型输入参数对输出响应相对影响的龙卷风图。
RΔ=|Δμ+μ|+|Δμμ-|
式中:Δμ+、Δμ、Δμ-分别为基于上界值、均值、下界值模型得到的易损性中位值。
第六步:基于中心点法的不确定性分析。
如果非线性函数为X=g(X1,X2,…,Xn),将X在各变量的均值点处展开成泰勒级数,并取线性项,采用差分法近似确定g(X)对Xi的梯度。
式中:Xi为建模随机参数;μXi分别为建模随机参数的均值、上界值、下界值;和/>分别为建模随机参数/>和/>时的易损性中位值。
采用下式求解建模相关的对数标准差βM,反映与结构建模相关的不确定性。
式中:σX为标准差,βMi为第i个建模随机变量的对数标准差。
将βM代入易损性函数得到下式考虑建模不确性的易损性函数。
相比于现有技术,本发明的一种考虑建模不确定性的结构地震易损性分析方法,能够对建模随机参数的敏感性进行排序,量化建模不确定性对地震易损性分析的影响。识别出的影响相对较小或不明显的建模参数有助于减少这些参数的计算工作量,确定的关键建模参数对于实现改进的地震易损性评估是很关键的,降低关键建模参数的不确定性可以提高的结构可靠性估计。
附图说明
图1是本发明一种考虑建模不确定性的结构地震易损性分析方法的流程图;
图2是本发明一实施例中双层框架式桥墩的计算模型图;图中,F—竖向荷载,P—水平荷载,1—节点,2—非线性梁柱单元,3—零长度截面单元;
图3是本发明合成的耐震时程曲线的示意图;图中,(a)为时程曲线,(b)为加速度反应谱;
图4是本发明结构模型在各极限状态下的龙卷风图;图中,(a)为结构模型在轻微损伤状态下的龙卷风图,(b)为结构模型在中度损伤状态下的龙卷风图,(c)为结构模型在重度损伤状态下的龙卷风图,(d)为结构模型在完全损伤状态下的龙卷风图;
图5是本发明考虑与未考虑建模不确定性的地震易损性曲线的对比图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
如图1所示,本实施例的一种考虑建模不确定性的结构地震易损性分析方法,包括:
步骤S1,不确定性建模参数的选取;
步骤S2,精细化有限元模型的建立;
步骤S3,基于耐震时程法的地震动输入;
步骤S4,基于云图法的地易损性分析;
步骤S5,基于龙卷风法的敏感性分析;
步骤S6,基于中心点法的不确定性分析。
该方法的具体步骤如下:
1.不确定性建模参数的选取
通过对既有研究成果的整理分析,确定了9个结构建模随机参数:阻尼比λw、容重系数ξ、混凝土28天抗压强度fc、混凝土峰值强度对应的应变ε0、混凝土残余强度fcu、混凝土残余强度对应的应变ε20、纵向钢筋初始弹性模量E0、纵向钢筋屈服强度fy、纵向钢筋应变硬化率B。随机参数的不确定性通常使用前两阶矩来描述,即均值和标准差或变异系数,假定这两个参数的随机分布特征详见下表,并分别采用第5个和第95个百分位的值作为上界值X+和下界值X-
2.精细化有限元模型的建立
以一座双层框架式桥墩为例,采用OpenSees开放平台建立双层框架式桥墩非线性有限元模型。无约束混凝土和约束混凝土均采用Kent-Park本构模型,钢筋采用Giuffre-Menegotto-Pinto本构模型,钢筋的粘结滑移采用Bond-SP01本构模型。截面采用纤维截面,分为纵向钢筋、非约束混凝土和约束混凝土。立柱和横梁均采用基于刚度法的梁柱单元,并采用零长度截面单元模拟立柱底部纵筋的粘结滑移。计算模型如图2所示。根据第一步的建模参数分析区间分别建立均值模型、上界值模型、下界值模型。
3.基于耐震时程法的地震动输入
初始地震动的输入采用由规范反应谱生成的人工地震波,取设计加速度反应谱最大值Smax=0.5g,特征周期Tg=0.35s,阻尼比ξ=5%。目标时间点为取10s,通过Matlab编程优化算法求解下式,合成总持续时间为30s的耐震时程曲线,如图3所示。
式中:SaT(T,t)为t时刻的加速度反应谱;T为结构的自振周期;t为任意时间点;tTarget为目标时间点;SaC(T)为目标加速度反应谱;SuT(T,t)为目标位移反应谱;ag为寻求的耐震加时程曲线;Sa(T,t)和Su(T,t)分别为ag在t时刻的加速度反应谱和位移反应谱;α为位移谱的相对权重系数。
4.基于云图法的地易损性分析
假设结构的需求和能力的概率特征由对数正态分布描述,结构需求参数的均值与地震动强度参数的关系满足幂指数关系,采用下列公式建立结构的地震易损性函数。
式中:Pf为结构达到或超过损伤极限状态的概率;Φ[]为标准正态分布函数;为结构需求参数的均值;βSd|IM为地震响应在给定IM范围内的条件对数标准差;/>βc分别为结构能力参数的均值和对数标准差;IM为地震动强度参数;a和b为未知的回归系数;Di为在第i条地震动作用下结构需求响应的峰值;N为进行非线性动力分析的次数。
将第三步得到的耐震时程曲线输入到第二步建立的数值模型中进行非线性时程分析得到结构的地震响应。钢筋混凝土桥墩的四种损伤极限状态分别对应为:纵筋首次屈服;混凝土开始剥落;水平承载力达到最大值;水平承载力下降至最大值的85%,定义结构的损伤极限状态如下表。以结构基本周期对应的5%阻尼比加速度反应谱Sa(T1,5%)为地震动强度参数,以位移延性比μΔ为结构需求参数,代入到地震易损性函数中得到易损性曲线。
5.基于龙卷风法的敏感性分析
以结构达到或超过某种损伤极限状态的概率为50%时所对应的易损性中位值为摆幅指标,根据第四步得到的易损性曲线采用下式计算摆幅值,根据摆幅值的大小从高到低进行排序,建立模型输入参数对输出响应相对影响的龙卷风图,如图4所示。
RΔ=|Δμ+μ|+|Δμμ-|
式中:Δμ+、Δμ、Δμ-分别为基于上界值、均值、下界值模型得到的易损性中位值。
6.基于中心点法的不确定性分析
如果非线性函数为X=g(X1,X2,…,Xn),将X在各变量的均值点处展开成泰勒级数,并取线性项,采用差分法近似确定g(X)对Xi的梯度。
式中:Xi为建模随机参数;μXi分别为建模随机参数的均值、上界值、下界值;和/>分别为建模随机参数/>和/>时的易损性中位值。
采用下式求解建模相关的对数标准差βM,反映与结构建模相关的不确定性。
式中:σX为标准差,βMi为第i个建模随机变量的对数标准差。
计算得建模相关的对数标准差βM,对应于轻微损伤、中度损伤、重度损伤、完全损伤分别是0.112、0.179、0.262、0.335,那么考虑建模不确定性后易损性对数标准差将分别提高6.4%,16.1%、13.3、20.1%。
将βM代入易损性函数得到下式考虑建模不确性的易损性函数,考虑与未考虑建模不确定性的易损性曲线的对比如图5所示。
本实施例的一种考虑建模不确定性的结构地震易损性分析方法,能够对建模随机参数的敏感性进行排序,量化建模不确定性对地震易损性分析的影响。识别出的影响相对较小或不明显的建模参数有助于减少这些参数的计算工作量,确定的关键建模参数对于实现改进的地震易损性评估是很关键的,降低关键建模参数的不确定性可以提高结构的可靠性估计。
本发明内容及具体实施例意在证明本发明所提供技术方法的实际应用,不应解释为对本发明保护范围的限定。在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。

Claims (7)

1.一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,包括如下步骤:
步骤1,不确定性建模参数的选取;
步骤2,精细化有限元模型的建立;
步骤3,基于耐震时程法的地震动输入;
步骤4,基于云图法的地易损性分析;
步骤5,基于龙卷风法的敏感性分析;
步骤6,基于中心点法的不确定性分析。
2.如权利要求1所述的一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,所述步骤1具体包括:
通过对既有研究成果的整理分析,确定结构建模随机参数,假定建模参数的随机分布特征,选择分析区间定出建模参数的均值、上界值和下界值。
3.如权利要求2所述的一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,所述步骤2具体包括:
选取合适的材料本构模型、截面模型和单元模型,采用有限元软件建立精细化结构数值模型;根据步骤1的建模参数分析区间分别建立均值模型、上界值模型和下界值模型。
4.如权利要求3所述的一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,所述步骤3具体包括:
选择目标反应谱和目标时间点,根据在某一时间内反应谱的大小与持续时间成线性关系,求解下列方程合成耐震时程曲线
式中:SaT(T,t)为t时刻的加速度反应谱;T为结构的自振周期;t为任意时间点;tTarget为目标时间点;SaC(T)为目标加速度反应谱;SuT(T,t)为目标位移反应谱;ag为寻求的耐震加时程曲线;Sa(T,t)和Su(T,t)分别为ag在t时刻的加速度反应谱和位移反应谱;α为位移谱的相对权重系数。
5.如权利要求4所述的一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,所述步骤4具体包括:
假设结构的需求和能力的概率特征由对数正态分布描述,结构需求参数的均值与地震动强度参数的关系满足幂指数关系,采用下列公式建立结构的地震易损性函数
式中:Pf为结构达到或超过损伤极限状态的概率;Φ[]为标准正态分布函数;为结构需求参数的均值;βSd|IM为地震响应在给定IM范围内的条件对数标准差;/>βc分别为结构能力参数的均值和对数标准差;IM为地震动强度参数;a和b为未知的回归系数;Di为在第i条地震动作用下结构需求响应的峰值;N为进行非线性动力分析的次数;
将步骤三得到的耐震时程曲线输入到步骤二建立的数值模型中进行非线性时程分析得到结构的地震响应。定义结构的损伤极限状态,选取合适的地震动强度参数和结构需求参数,代入到地震易损性函数中得到易损性曲线。
6.如权利要求5所述的一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,所述步骤5具体包括:
以结构达到或超过某种损伤极限状态的概率为50%时所对应的易损性中位值为摆幅指标,根据步骤4得到的易损性曲线采用下式计算摆幅值,根据摆幅值的大小从高到低进行排序,建立模型输入参数对输出响应相对影响的龙卷风图
RΔ=Δμ+μμμ-
式中:Δμ+、Δμ、Δμ-分别为基于上界值、均值、下界值模型得到的易损性中位值。
7.如权利要求6所述的一种考虑建模不确定性的结构地震易损性分析方法,其特征在于,所述步骤6具体包括:
如果非线性函数为X=g(X1,X2,…,Xn),将X在各变量的均值点处展开成泰勒级数,并取线性项,采用差分法近似确定g(X)对Xi的梯度,i=1~n,
式中:Xi为建模随机参数;μXi分别为建模随机参数的均值、上界值、下界值;/>和/>分别为建模随机参数/>和/>时的易损性中位值;
采用下式求解建模相关的对数标准差βM,反映与结构建模相关的不确定性
式中:σX为标准差;βMi为第i个建模随机变量的对数标准差;
将βM代入易损性函数得到下式考虑建模不确性的易损性函数
CN202310413818.5A 2023-04-15 2023-04-15 一种考虑建模不确定性的结构地震易损性分析方法 Pending CN116595821A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310413818.5A CN116595821A (zh) 2023-04-15 2023-04-15 一种考虑建模不确定性的结构地震易损性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310413818.5A CN116595821A (zh) 2023-04-15 2023-04-15 一种考虑建模不确定性的结构地震易损性分析方法

Publications (1)

Publication Number Publication Date
CN116595821A true CN116595821A (zh) 2023-08-15

Family

ID=87603475

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310413818.5A Pending CN116595821A (zh) 2023-04-15 2023-04-15 一种考虑建模不确定性的结构地震易损性分析方法

Country Status (1)

Country Link
CN (1) CN116595821A (zh)

Similar Documents

Publication Publication Date Title
Lignos et al. Development and utilization of structural component databases for performance-based earthquake engineering
Lagaros et al. Fragility assessment of steel frames using neural networks
Alembagheri et al. Seismic performance sensitivity and uncertainty analysis of gravity dams
Faggella et al. Probabilistic seismic response analysis of a 3-D reinforced concrete building
Brozovič et al. Envelope‐based pushover analysis procedure for the approximate seismic response analysis of buildings
Hancilar et al. Fragility functions for code complying RC frames via best correlated IM–EDP pairs
Salvatori et al. Probabilistic seismic performance of masonry towers: General procedure and a simplified implementation
Heo et al. Damage-based seismic performance evaluation of reinforced concrete frames
Hariri-Ardebili et al. Single and multi-hazard capacity functions for concrete dams
CN114792020B (zh) 一种基于机器学习的建筑抗震韧性快速评估方法及系统
CN117195594B (zh) 隧道岩爆等级评估方法、装置、电子设备及存储介质
Feng et al. Fiber-based modeling of circular reinforced concrete bridge columns
Li et al. Efficient structural seismic performance evaluation method using improved endurance time analysis
Li et al. Fast seismic response estimation of tall pier bridges based on deep learning techniques
CN116467789A (zh) 一种基于结构应变能的建筑结构整体损伤量化评估方法
Wahrhaftig et al. Buckling and free vibration analysis of non-prismatic columns using optimized shape functions and Rayleigh method
CN108763164B (zh) 煤与瓦斯突出反演相似度的评价方法
Sheikh et al. Seismic fragility analysis using nonlinear autoregressive neural networks with exogenous input
Aghaeipoor et al. Seismic damage of submerged intake tower under the sequence of mainshocks and aftershocks
CN113011053A (zh) 一种基于性能裕量的结构确信可靠性分析方法及系统
Wang et al. Role of strong motion duration on seismic responses of high concrete faced rockfill dams
CN116595821A (zh) 一种考虑建模不确定性的结构地震易损性分析方法
CN110208087B (zh) 脉动压裂循环载荷下岩石强度计算方法
CN112733395A (zh) 型钢混凝土组合构件从损伤到失效的综合评价分析方法
Fitzjohn et al. A combined SHM/IDA method for assessing collapse capacity and risk in subsequent ground motions

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