一种基于非平稳随机核分布函数的地震可靠度计算方法技术

本发明专利技术公开了一种基于非平稳随机核分布函数的地震可靠度计算方法，包括以下步骤：(1)随机抽样生成结构样本；(2)建立结构的确定性模型，并生成非平稳随机地震动样本；(3)随机地震动样本调幅与循环；(4)将所有地震动水准下的各条地震动波的强度参数IM汇总，并生成IM的一维概率分布曲线；(5)将所有地震动水准下的各条波的强度参数IM和工程需求参数EDP汇总，生成IM

【技术实现步骤摘要】
一种基于非平稳随机核分布函数的地震可靠度计算方法


[0001]本专利技术涉及土木工程概率地震性能评估方法，尤其涉及一种基于非平稳随机核分布函数的地震可靠度计算方法。

技术介绍

[0002]20世纪90年代以来，基于性能的地震工程逐渐成为国际上地震工程领域的主要发展趋势。基于性能的地震工程发展可分为三个阶段，第一阶段以确定性的结构性态反应为性能目标，第二阶段以结构整体可靠度为性能目标，第三阶段即目前以控制地震风险和地震损失为性能目标。现阶段的地震可靠度分析方法为参数型，即通常假设其变量符合对数正态分布。然而，当地震动较大、结构倒塌点较多时，将不再满足对数正态分布假设，基于该假设的可靠度评估结果失真严重。

技术实现思路

[0003]专利技术目的：本专利技术旨在提供一种避免评估结果失真、提高计算精度的基于非平稳随机核分布函数的地震可靠度计算方法。
[0004]技术方案：本专利技术所述的一种基于非平稳随机核分布函数的地震可靠度计算方法，包括以下步骤：
[0005](1)确定结构变量的类型、分布与参数，并随机抽样生成结构样本n1；
[0006](2)建立结构的确定性模型，并生成非平稳随机地震动样本n2；引入非平稳随机核分布函数Sginp(t,β)，生成非平稳随机地震动ginp(θ,t)；
[0007](3)随机地震动样本调幅与循环，开展非线性动力时程响应分析，得到不同地震动水准下的强度参数IM和工程需求参数EDP；
[0008](4)将所有地震动水准下的各条地震动波的强度参数IM汇总，生成强度参数IM的基于非平稳随机核分布函数的一维概率分布曲线；
[0009](5)将所有地震动水准下的各条非平稳随机地震动样本的强度参数IM和工程需求参数EDP汇总，生成强度参数IM和工程需求参数EDP的基于非平稳随机核分布函数的二维联合概率分布曲线；
[0010](6)基于非平稳随机核分布函数的概率函数，计算地震可靠度。
[0011]优选地，所述步骤(1)中所述参数包括结构荷载、结构尺寸和结构材料。
[0012]优选地，所述步骤(2)包括：确定非平稳随机地震动的强度参数IM和工程需求参数EDP类型；确定非平稳随机地震动的强度参数调幅间隔与调幅数量m。
[0013]优选地，所述步骤(2)中非平稳随机核分布函数S
ginp
(t,β)为：
[0014][0015]其中，t表示激励时间，β表示激励频率，β
g
(t)表示场地分类非平稳参数，β
f
(t)表示地震动频率非平稳参数，ξ
g
(t)表示地震动分组非平稳参数，ξ
f
(t)表示地震动折减非平稳参
数，A
amp
(t)表示地震动峰值加速度到达时间的强度参数，S
amp
(t)表示地震动等效峰值强度参数；
[0016]非平稳随机地震动g
inp
(θ,t)为：
[0017][0018]其中，β
k
表示第k项激励频率，β
if
表示间隔频率，N
tr
表示截断项数，Θ1
k
表示地震动随机相位角一，Θ2
k
表示地震动随机相位角二，sin表示正弦函数，cos表示余弦函数。
[0019]优选地，所述步骤(3)包括：
[0020](3.1)设定初始参数y＝1，确定此时的地震动强度参数IM_y，在第y水准下，随机匹配n1和n2，开展非线性动力时程响应分析，得到第y水准下的工程需求参数EDP_
y
；
[0021](3.2)判断y与m的大小，若y小于m，则y＝y+1，将非平稳随机地震动g
inp
(θ,t)样本调幅至第y水准，并返回步骤(3.1)循环；若y大于等于m，则进入步骤(4)。
[0022]优选地，所述步骤(4)中强度参数IM的基于非平稳随机核分布函数的一维概率分布曲线K
X
(x)为
[0023][0024]其中N＝n2×
m，x
i
表示第i条非平稳随机地震动样本的强度参数IM的值；b表示基于非平稳随机核分布函数的一维估计带宽，所述一维估计带宽b为：
[0025]b＝S
ginp
(t，β)
×
σ
×
(n2×
m)
‑
0.4
[0026]其中σ为各条地震动波IM的标准差。
[0027]优选地，所述步骤(5)中强度参数IM和工程需求参数EDP的基于非平稳随机核分布函数的二维联合概率分布曲线f
X
(x)为：
[0028][0029]其中N＝n2×
m，d＝2，X
i
表示二维向量(IM，EDP)的点对形式，T为转置运算，B
‑1为B的逆矩阵；B表示基于非平稳随机核分布函数的二维估计带宽，所述二维估计带宽B为：
[0030][0031]其中，Φ
‑1表示标准正态分布函数的逆函数，π表示圆周率，arc表示求逆，PA1表示二维估计带宽协相关系数，n＝n1
×
n2。
[0032]优选地，所述步骤(6)中基于非平稳随机核分布函数的概率函数E
X
(x)为：
[0033][0034]其中，Q
X
(x)表示强度参数IM和工程需求参数EDP的基于非平稳随机核分布函数的二维联合概率分布曲线对工程需求参数EDP单变量积分。
[0035]优选地，所述Q
X
(x)为
[0036][0037]其中，δc表示可靠度极限状态限值。
[0038]优选地，所述步骤(6)中地震可靠度F(a,δ
c
)为
[0039][0040]其中a和δ为地震可靠度F(a,δ
c
)的变量，分别表示地震动强度和工程需求，IM
i
和Δ
i
分别为IM和EDP的第i组数值，F(a,δ
c
)表示地震动强度为a、可靠度极限状态限值为δ
c
情况下基于非平稳随机核分布函数的地震可靠度计算结果。
[0041]有益效果：与现有技术相比，本专利技术具有如下显著优点：(1)避免了传统的对数正态分布条件假设，同时避免了在地震动较大、结构倒塌点较多时地震可靠度评估结果的失真风险；(2)通过计算强度参数和需求参数间的一维及二维带宽矩阵，实现了在少量样本点情况下的地震可靠度计算，避免了传统蒙特卡洛抽样可靠度方法的计算负担。
附图说明
[0042]图1为本专利技术中一种基于非平稳随机核分布函数的地震可靠度计算方法的流程图；
[0043]图2为本专利技术实施例中结构示意图；
[0044]图3为本专利技术实施例中结构模拟示意图；
[0045]图4为本专利技术实施例中基于本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于非平稳随机核分布函数的地震可靠度计算方法，其特征在于，包括以下步骤：(1)确定结构变量的类型、分布与参数，并随机抽样生成结构样本n1；(2)建立结构的确定性模型，并生成非平稳随机地震动样本n2；引入非平稳随机核分布函数S
ginp
(t,β)，生成非平稳随机地震动g
inp
(θ,t)；(3)随机地震动样本调幅与循环，开展非线性动力时程响应分析，得到不同地震动水准下的强度参数IM和工程需求参数EDP；(4)将所有地震动水准下的各条地震动波的强度参数IM汇总，生成强度参数IM的基于非平稳随机核分布函数的一维概率分布曲线；(5)将所有地震动水准下的各条非平稳随机地震动样本的强度参数IM和工程需求参数EDP汇总，生成强度参数IM和工程需求参数EDP的基于非平稳随机核分布函数的二维联合概率分布曲线；(6)基于非平稳随机核分布函数的概率函数，计算地震可靠度。2.根据权利要求1所述的地震可靠度计算方法，其特征在于，所述步骤(1)中所述参数包括结构荷载、结构尺寸和结构材料。3.根据权利要求1所述的地震可靠度计算方法，其特征在于，所述步骤(2)包括：确定非平稳随机地震动的强度参数IM和工程需求参数EDP类型；确定非平稳随机地震动的强度参数调幅间隔与调幅数量m。4.根据权利要求3所述的地震可靠度计算方法，其特征在于，所述步骤(2)中非平稳随机核分布函数S
ginp
(t,β)为：其中，t表示激励时间，β表示激励频率，β
g
(t)表示场地分类非平稳参数，β
f
(t)表示地震动频率非平稳参数，ξ
g
(t)表示地震动分组非平稳参数，ξ
f
(t)表示地震动折减非平稳参数，A
amp
(t)表示地震动峰值加速度到达时间的强度参数，S
amp
(t)表示地震动等效峰值强度参数；非平稳随机地震动g
inp
(θ,t)为：其中，β
k
表示第k项激励频率，β
if
表示间隔频率，N
tr
表示截断项数，Θ1
k
表示地震动随机相位角一，Θ2
k
表示地震动随机相位角二，sin表示正弦函数，cos表示余弦函数。5.根据权利要求4所述的地震可靠度计算方法，其特征在于，所述步骤(3)包括：(3.1)设定初始参数y＝1，确定此时的地震动强度参数IM_y，在第y水准下，随机匹配n1和n2，开展非线性动力时程响应分析，得到第y水准下的工程需求...

【专利技术属性】
技术研发人员：曹徐阳，沈德建，汪鑫，李明，纵达，马天晓，
申请(专利权)人：河海大学，
类型：发明
国别省市：