结构的可靠性包括安全性、适用性和耐久性, 工程上以可靠度衡量结构可靠性, 即结构在规定的时间、规定的条件下完成预定功能的概率[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) |
当结构受两个以上的荷载作用时, 各个荷载同时出现最大值的可能性较小, 因此多载荷作用时需乘以相应的系数进行折减[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所示.结构在每一时段内的功能函数为
(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所示.
将结构承载力极限状态函数当量正态化可得
(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阶协方差矩阵;
计算结构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中纵坐标表示荷载均值分布情况, 用vα(ti)表示单位时间内对跨阈值α向上跨阈次数的期望.由Rice公式可知
(28) |
式中,
(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所示.
假设结构每次损伤后阈值的降低程度相同, 则
(46) |
某钢筋混凝土梁承受自重和两种活荷载共同作用, 梁抗力服从对数正态分布, 抗力概率特性如表 1所示.
假设结构受到的重力服从正态分布, 为平稳随机过程, 其均值为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可知:当结构可靠性逐渐降低时, 加载历史的影响会逐渐增大.
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 ) |