康玉梅
, 刘子傲, 吴鹏飞
东北大学 资源与土木工程学院, 辽宁 沈阳 110819
收稿日期：2018-08-23
基金项目：国家自然科学基金资助项目(U1602232，51525402，51374049，51574060)。
作者简介：康玉梅(1973-)，女，辽宁灯塔人，东北大学院副教授。
摘要：针对JCSS推荐的荷载效应组合法假设荷载效应之间既相互独立又完全相关的内部矛盾和时段分析法对工程中的非平稳荷载效应过程无法实现有效分析的问题, 对考虑荷载作用顺序的结构可靠度分析方法进行了研究.首先提出了非平稳荷载效应过程下的结构可靠度分析方法; 其次考虑荷载作用顺序的影响, 改进了结构在受多时段荷载作用下的可靠度计算模型; 最后运用跨阈理论分析了结构损伤对可靠度的影响, 推导出综合考虑非平稳荷载过程、荷载作用顺序和结构损伤的可靠度计算公式.算例分析表明:该方法适用性更广泛, 计算结果更真实可靠.
关键词：结构可靠度时段分析法非平稳荷载效应荷载作用顺序荷载组合
Structural Reliability Analysis Method Considering Load Sequence
KANG Yu-mei
, LIU Zi-ao, WU Peng-fei
School of Resources & Civil Engineering, Northeastern University, Shenyang 110819, China
Corresponding author: KANG Yu-mei, E-mail:
kangyumei@mail.neu.edu.cnAbstract: In view of the contradiction which load effects are assumed to be both independent and completely related in load-effect combination method recommended by JCSS, and the ineffectiveness of time-interval analysis method for non-stationary load-effect process in engineering, structural reliability analysis method considering load sequence is performed. Firstly, the structural reliability analysis method under non-stationary load effect process is proposed. Secondly, the structural reliability calculation model under multi-period load considering load sequence is improved. Finally, the influence of structural damage on reliability is analyzed by cross-threshold theory, and the formula for calculating structural reliability considering non-stationary load process, loading sequence and structural damage is derived. The results showed that the method is more applicable and the calculation results are more realistic and reliable.
Key words: structural reliabilitytime period analysis methodnon-stationary load effectload sequenceload combination
结构的可靠性包括安全性、适用性和耐久性, 工程上以可靠度衡量结构可靠性, 即结构在规定的时间、规定的条件下完成预定功能的概率
[1-3].
荷载效应组合在结构可靠度分析和结构设计中占有重要地位
[4].1976年由“国际结构安全度联合委员会(JCSS)”公布了“对各类结构和各种荷载的统一规则”, 即JCSS组合规则, 当作用于结构上的各环境荷载可视为相互独立的随机过程时, 即可用JCSS组合方法进行荷载效应组合, 并将计算二阶可靠指标的有效数值过程
[5]确立为可靠指标的标准算法(JCSS法).但是, 这种方法的假设条件要求荷载效应之间相互独立且完全相关, 这一内部矛盾，使得在计算结构可靠度时结果与真实情况并不完全相符.时段分析法较JCSS法更符合结构作用的真实性
[6].但是, 时段分析法只考虑了平稳荷载效应过程, 未考虑荷载作用时间顺序的影响.在实际工程中, 很多荷载都具有一定的非平稳性, 而且荷载的作用顺序对结构必然存在影响.
1 结构失效概率的时段分析法1.1 结构的失效概率结构的功能函数
Z[7]可表示为
 | (1) |
式中:
R表示结构抗力;
S表示荷载效应.
Z>0时, 结构可靠;
Z < 0时, 结构失效;
Z=0时, 结构处于极限状态.功能函数中的自变量为不同的随机变量, 具有一定的不确定性, 因此用概率来度量结构的可靠性是最好的方式.
Z≥0时的概率
ps为可靠概率, 反之
pf为失效概率, 即
 | (2) |
1.2 JCSS推荐的荷载效应组合方法当结构受两个以上的荷载作用时, 各个荷载同时出现最大值的可能性较小, 因此多载荷作用时需乘以相应的系数进行折减
[8].在JCSS推荐的荷载组合方法中, 荷载效应被假定为等时段的平稳二项矩形波过程, 并且荷载与荷载效应之间有明确的线性关系, 即
 | (3) |
式中:
ki, Δ
Ti,
τij,
Si(
τij)分别为
Si(
t)的时段数、时段长度、第
j个时域以及在该时域的最大值.
目标使用期
τ内的最大荷载组合效应为
 | (4) |
按照JCSS推荐的方法, 最大荷载效应组合为下述
n种情况中的最大值.
 | (5) |
式中，
Si(Δ
Ti),
Si(Δ
Ti－1),
Si(
T)为
Si(
t)在时段Δ
Ti, Δ
Ti－1,
T内的最大值, 但式(5)中第
i种组合是作了如式(3)所示的假定得到的.
 | (6) |
式中:
Sj(
τj－1, l)为
Si(
t)在时段
τj－1, l内的最大值(
l=1, 2, …,
kj－1).按照JCSS法的概率模型,
Sj(
τj1),
Sj(
τj2), …,
Sj(
τjkj)之间以及
Sj(
τj－1, 1),
Sj(
τj－1, 2), …,
Sj(
τj－1, kj－1)之间都是相互独立的, 而JCSS法中式(6)又假定它们完全相关, 这会导致其计算结果与实际结果误差较大, 最大荷载效应组合的分析结果也会因荷载效应之间相关性的增强而偏小, 进行结构设计则偏于不安全.
1.3 时段分析法在时段分析法中, 假定抗力为随机变量, 荷载效应为
n个平稳随机过程.这里以结构受三种荷载效应为例, 其随时间变化的样本函数如
图 1所示.结构在每一时段内的功能函数为
图 1(Fig. 1)
 | (7) |
假设荷载效应为等时段的平稳二项矩形波过程, 则结构每个时段内的失效概率或可靠指标都相同, 结构在目标使用期
τ内的失效概率
[9]为
 | (8) |
式中:
β为功能函数
H(
τ3i)的可靠指标;
ρ为
H(
τ3i)的相关系数.
 | (9) |
式中:
i,
l=1, 2, 3;
H(·)为正态化和线性化后的功能函数;
YR,
YSl为标准正态随机变量;
φ(·)为标准正态密度函数;
Φ-1(·)为标准正态分布函数的反函数;
fx(·)为随机变量
x的概率密度函数;
Fx(·)为随机变量
x的概率分布函数;
R*,
Sl*为设计验算点;
Φn为
n维标准正态函数.
确定可靠指标和相关系数后, 结构的失效概率便可用式(8)进行计算, 它与抗力、荷载效应的概型是完全相符的, 且不必采用近似的方法进行荷载效应组合, 以及寻找最不利的荷载效应组合.
2 考虑荷载作用顺序的可靠度分析2.1 非平稳荷载过程可靠度分析方法假设结构的抗力服从正态分布, 结构的荷载过程为非平稳随机过程, 在不同的时段有着不同的均值和方差.荷载与抗力相互独立, 且在不同时段荷载随机过程有不同的均值与方差.假定荷载效应的均值样本函数如
图 2所示.
图 2(Fig. 2)
将结构承载力极限状态函数当量正态化可得
 | (10) |
其正态化后的等效均值和标准差为
 | (11) |
相应地, 结构抗力与荷载效应可表示为
 | (12) |
式中:
μ,
x*,
β可通过反复迭代求得.当迭代后的‖
x‖满足精度要求时, 可求相应的可靠指标, 即
β=(
β1,
β2, …,
βn)
T.假设荷载随机过程为独立增量过程, 此时不同时段间的相关系数为
[10] | (13) |
将结构使用期内荷载效应分为
n个时段, 结构在使用期内的失效概率为
 | (14) |
当结构各时段失效模式之间的相关系数相同时, 式(14)可简化为
 | (15) |
其中,
φ(
x)是标准正态分布概率密度函数.如果各时段失效模式的可靠指标也相等, 则式(15)可进一步简化为
[6] | (16) |
但是实际工程中很难保证各时段之间有相同的相关系数或可靠指标, 因此求解结构使用期内的失效概率往往只能采用式(17)的方法进行计算:
 | (17) |
在计算过程中涉及到多维积分的求解, 这不但给求解造成较大困难, 而且进行数值积分的精度也难以得到保证.因此, 需要对计算方法进行简化.求解时先求得前两个时段的失效概率, 由于两个时段之间只有一个相关系数, 所以可视为各时段失效模式之间相关系数相同的情况, 按式(18)进行求解.
 | (18) |
则这两时段的可靠指标：
 | (19) |
同理, 可求得
β34,
β56, …这样便将原来的
n个时段缩减到原来的一半.依此类推, 最后可求出结构在使用期内的失效概率和可靠指标, 避免了高维积分的计算困难与精度问题.
2.2 考虑荷载作用顺序的条件可靠度分析方法假设抗力服从正态分布, 不随时间变化, 且与荷载效应相互独立.荷载效应随时间变化服从不同的概率分布, 将荷载效应分为几个相等时段, 假设每个时段内荷载效应的分布函数不变, 荷载随机变量服从联合分布为正态的时间序列, 即
 | (20) |
其中:
R为
n×
n阶协方差矩阵;
,
Aij为
R中第
ij个元素的代数余子式.
计算结构
t时间内可靠度时, 将
n个时段内的荷载效应进行联合求解, 对
n维联合密度函数积分求失效概率, 即
 | (21) |
这种方法虽然在一定程度上简化了计算, 但未考虑时间的先后, 若将某两个时段的荷载效应互换后, 所得的
t时间内的失效概率并不发生变化, 这显然与真实情况有偏差.为了更准确地求解结构的可靠度, 基于正面信息的条件可靠度大于无条件概率, 基于负面信息的条件可靠度小于无条件概率
[11]，这里, 前一时段荷载可视为基于负面信息的条件可靠度, 每一个时段的失效概率都是建立在前一时段的基础上的条件概率.
假设荷载效应在
n个时间序列内呈正态分布, 则其在第一时段内的极限状态方程为
Z=
R－
S1.由于相互独立且服从正态分布的随机变量的线性组合仍服从正态分布, 所以在第一时段内
 | (22) |
求解第二时段的失效概率时, 应在第一时段荷载作用的条件下进行第二次荷载作用, 即
 | (23) |
而不是单纯的
P{
Z(
X1,
X2) < 0}.同理, 之后每个时段的失效概率为
 | (24) |
显然, 过去考虑时间序列上的荷载效应的方法偏于冒进, 且随着时间区段的增加, 产生的误差逐渐加大, 这会导致在考虑时间变量时的可靠指标与实际有较大偏差, 而用本文的方法可以更精确地求得结构在荷载效应随时间变化情况下的可靠度.
2.3 计算简化方法每个时段内荷载的概率密度函数都相同, 将每相邻两时段的荷载效应和失效概率整合, 依此类推, 求得总的失效概率.下面以
图 2所示荷载效应的均值变化情况进行说明.
首先, 求
t1,
t3,
t5,
t7时段的失效概率:
 | (25) |
其次, 分别在
t1,
t3,
t5,
t7时段条件下求
t2,
t4,
t6,
t8时段的失效概率, 每两时段的荷载效应服从二维正态分布.同理, 可求出
t2,
t4,
t6,
t8时段的失效概率
pf2,
pf4,
pf6,
pf8.按同样的方法可以继续对时段进行整合, 最后可求出结构在整个时间
t内的可靠度.
2.4 荷载效应组合假设结构受到三个荷载效应共同作用, 每个时段内荷载效应的概率密度函数相同, 如
图 1所示, 其纵坐标为荷载效应在各时段内的均值.
首先, 对最小荷载时段进行整合, 求出
S3(
τ3i)时段整合后的失效概率:
 | (26) |
其次, 在进行下一步时段整合前, 先对
τ2i上各时段的
S3与
S2进行荷载组合, 求出考虑
S3与
S2两种荷载效应下的
τ2i上各时段失效概率:
 | (27) |
最后, 继续向下一时段整合, 进入下一时段后, 采用同样的方法先对
S1进行荷载效应组合.依此类推, 得出总的失效概率
pf.
该方法考虑了荷载效应随时间变化的先后顺序, 每次计算只用到二维正态分布的概率求解.当时间区段较多时, 计算步骤也较多, 但很容易通过编程实现.此外, 该方法在荷载组合时具有更高的精度.
3 抗力随时间变化的可靠度分析在讨论荷载概率模型时, 常采用基于随机过程的跨阈理论来估计使用期荷载的最大值.但是, 跨阈理论模型不限于估计荷载概率模型的最大值分布, 实际应用时可定义一个阈值
α, 当荷载值超过阈值
α时, 会对结构造成一定损伤, 使抗力
R有一定程度的削减, 而导致可靠指标降低.若无视荷载作用的时间顺序, 会增大计算误差.每次荷载超过阈值
α时, 抗力会产生损伤, 损伤后新的阈值
α也随之减小, 如
图 3所示.因此, 同样大小的荷载效应值
S, 在
t1时段可能不会造成结构的损伤, 但如果发生在
t2时段, 就可能对结构抗力产生影响, 所以这时荷载作用顺序对结构的可靠性影响很大.
图 3(Fig. 3)
本文考虑荷载在每一时段内平稳分布的情况,
图 3中纵坐标表示荷载均值分布情况, 用
vα(
ti)表示单位时间内对跨阈值
α向上跨阈次数的期望.由Rice公式可知
 | (28) |
式中,
为
S的导数, 则
ti时段荷载跨阈次数为
 | (29) |
设
ti时段内荷载超过阈值
α的总时间为
Tαi, 且每个时段内荷载效应都是平稳的, 则
Tαi的均值:
 | (30) |
式中,
Fs(·)为
S的分布函数.设一次跨阈的持续时间为
dα, 时段
ti内总跨阈次数为
Nαti, 显然
 | (31) |
 | (32) |
 | (33) |
如果荷载发生变化的概率与时间无关, 那么荷载发生变化的时点可用泊松过程来描述, 即荷载在时间
t以前发生变化的次数
N(
t)服从泊松分布, 荷载在变化
k次后开始首次超越阈值的概率为
 | (34) |
式中,
λ为单位时间荷载发生变化次数的期望与方差.若荷载超过阈值
α变化
k次以后降到
α以下, 则跨阈持续时间的分布函数和概率密度为
 | (35) |
 | (36) |
则跨阈持续时间的期望：
 | (37) |
又由式(35)可得
 | (38) |
 | (39) |
对于平稳离散型的荷载效应, 每次荷载状态发生变动前荷载大小保持恒定值, 此时上跨域次数为
 | (40) |
每一时段内荷载为平稳过程, 分布函数相同, 则
 | (41) |
 | (42) |
 | (43) |
 | (44) |
为简化模型, 假设每次荷载跨阈时抗力折减都相同, 且每次跨阈时长都取均值
μdαti, 每次荷载超过
α值时抗力
R都乘以折减系数
λ, 则
 | (45) |
式中:
λ(Δ
ti)为每次跨阈持续时间Δ
t下的结构损伤情况, 是与跨阈时长有关的函数;
k是荷载超过阈值
α后变化的次数.设结构受多时段荷载作用, 荷载每次超过阈值后结构都会发生损伤, 荷载与阈值的变化情况如
图 4所示.
图 4(Fig. 4)
假设结构每次损伤后阈值的降低程度相同, 则
 | (46) |
4 算例分析某钢筋混凝土梁承受自重和两种活荷载共同作用, 梁抗力服从对数正态分布, 抗力概率特性如
表 1所示.
表 1(Table 1)
 表 1 结构抗力概率特性Table 1 Probabilistic characteristics of structural esistance | 概率特性 | 1 | 2 | 3 | 4 | 5 | | 均值/kN | 42 | 41.8 | 41.2 | 40.6 | 40.2 | | 标准差 | 5.46 | 5.532 | 5.604 | 5.676 | 5.748 |
| 表 1 结构抗力概率特性 Table 1 Probabilistic characteristics of structural esistance |
假设结构受到的重力服从正态分布, 为平稳随机过程, 其均值为5 kN, 变异系数为0.070, 结构受到的两种活荷载都服从极值Ⅰ型分布, 活荷载的均值为7.16 kN, 方差为13.663, 风荷载的均值为10.1 kN, 方差为28.282.结构的可靠指标计算结果如
图 5中情况1所示.
图 5中情况2表示不考虑抗力劣化; 情况3表示不考虑历史荷载的影响; 情况4表示既不考虑抗力劣化, 也不考虑历史荷载的影响.对比情况3与情况4可知:考虑结构损伤与不考虑结构损伤的分析结果有很大差异; 对比情况2和情况4可知:在结构本身可靠性较高时, 是否考虑加载历史会对计算结果产生影响, 但影响不是十分明显; 对比情况1和情况3可知:当结构可靠性逐渐降低时, 加载历史的影响会逐渐增大.
图 5(Fig. 5)
5 结论1) 对非平稳随机荷载作用下的结构可靠度进行了研究, 给出了结构承受非平稳随机荷载时的计算公式; 针对求解时可能遇到的多维积分问题, 给出了相应的简化方法和结构受非平稳荷载组合时的求解步骤.
2) 鉴于结构受多时段荷载作用时历史荷载对后续时段有影响, 研究了可以综合考虑各时段荷载之间相互影响的结构可靠度分析方法和考虑荷载作用顺序变化对结构可靠度产生的影响.
3) 基于跨阈理论考虑抗力随时间的变化情况, 给出了综合考虑跨阈次数、跨阈持续时间、抗力折减程度的结构可靠度计算方法, 根据实际计算中的不同需求, 代入不同的损伤模型进行计算.
参考文献 | [1] | Yao J T, Cheng K K. Discussion of sensitivity factors and dominating variable for design value method[J].Journal of Civil Engineering, 2017, 21(1): 37–44.
|
| [2] | Tabandeh A, Gardoni P, Murphy C. A reliability-based capability approach[J].Risk Analysis, 2017, 38(2): 410–424.
|
| [3] | Wang C, Li Q W, Pang L, et al. Estimating the time-dependent reliability of aging structures in the presence of incomplete deterioration information[J].Journal of Zhejiang University(Science A), 2016, 17(9): 677–688.DOI:10.1631/jzus.A1500342
|
| [4] | Yao J T, Gu H. Probability model and combination method of crane load of industrial building[J].Journal of Building Structures, 2016, 37(11): 160–166.
|
| [5] | Song S F, Lu Z Z, Qiao H W. Subset simulation for structural reliability sensitivity analysis[J].Reliability Engineering & System Safety, 2017, 94(2): 658–665.
|
| [6] | 姚继涛. 结构可靠度的时段分析方法[J].土木工程学报, 2005, 38(7): 1–5.
( Yao Ji-tao. Time interval analysis method of structural reliability[J].Journal of Civil Engineering, 2005, 38(7): 1–5.DOI:10.3321/j.issn:1000-131X.2005.07.001)
|
| [7] | Han X, Jiang C, Liu L X, et al. Response-surface-based structural reliability analysis with random and interval mixed uncertainties[J].Science China, 2014, 57(7): 1322–1334.DOI:10.1007/s11431-014-5581-6
|
| [8] | Ching J, Phoon K K, Yang J J. Role of redundancy in simplified geotechnical reliability-based design-a quantile value method perspective[J].Structural Safety, 2015, 55(1): 37–48.
|
| [9] | Wang P, Zhang J, Zhai H, et al. A new structural reliability index based on uncertainty theory[J].Chinese Journal of Aeronautics, 2017, 30(4): 1451–1458.DOI:10.1016/j.cja.2017.04.008
|
| [10] | 李全旺, 王草. 荷载随机过程相关性对结构时变可靠度的影响[J].清华大学学报(自然科学版), 2014, 54(10): 1316–1320.
( Li Quan-wang, Wang Cao. The effect of correlation of load random process on time-varying reliability of structures[J].Journal of Tsinghua University(Natural Science Edition), 2014, 54(10): 1316–1320.)
|
| [11] | 张庆功, 张鹏, 胡启国. 关于条件可靠度不等式的证明[J].哈尔滨工业大学学报, 2005, 37(2): 212–213.
( Zhang Qing-gong, Zhang Peng, Hu Qi-guo. Proof of inequality of conditional reliability[J].Journal of Harbin University of Technology, 2005, 37(2): 212–213.DOI:10.3321/j.issn:0367-6234.2005.02.020)
|
