基于AK-FORM方法和降維方法的高效時變可靠度分析方法

關(guān)鍵詞：時變可靠度;PHI2方法；一次可靠度方法；Kriging模型;跨越率；降維方法中圖分類號：TB114.3；TU318 文獻(xiàn)標(biāo)志碼：A 文章編號：2096-6717（2025）03-0162-09

An efficient time-varying reliability method based on AK-FORM method and dimension reductionmethod

ZHANGLiang1，TANGYake1，NIUKai'，LIHaoyu ² ，YU Shujun2 （1.State Grid Henan Economic Research Institute， Zhengzhou 45oo52，P.R.China；2.Schoolof Civil Engineering， Chongqing University， Chongqing 40oo45，P.R. China）

Abstract： The PHI2 method is commonly used to perform the time-varying structural reliability analysis，with thecalculation of outcrossing rate being pivotal to its implementation.Achieving suficient accuracy is often necessary to calculate the outcrossng rate at a large number of moments.However，for practical problems involving complex limit state surfaces，the calculation of the outcrossing rate at each moment can be very timeconsuming.To further improve the eficiency of PHI2 method，three strategies are proposed to be introduced in this paper to improve the eficiencyof calculating outcrossing rate.Firstly，the strategy without Cholesky decomposition is used to reduce the number of random variables，while the corrsponding calculation of correlation coeficients is given. Secondly，the improved first-order reliability method based on the adaptive Kriging model（AK-FORM） is introduced to efficiently calculate the reliability index at each moment.Finally， the two-dimensional integral isconverted into aone-dimensional integral byusing the dimensionreduction method.The threeimprovement strategies are combined with thePHI2 method，which forms an eficient timevarying reliability analysis method based on the AK-FORM method and the dimension reduction method，i.e.， the K-PHI2 method.Meanwhile，only the strategy without Cholesky decomposition is combined with the PHI2 method to form the PHI2- method.The calculation results of numerical and engineering examples demonstrate that the PHI2- and K-PHI2 methods proposed in this paper have the same high accuracy as the PHI2 method， and both are better than the PHI2+ method （an improved method based on PHI2） in terms of accuracy; compared with the PHI2 and PHI2 + methods，the PHI2-method has a little improvement in efficiency，while the K-PHI2 method further greatly improves the eficiency of time-varying reliability analysis on this basis.

Keywords： time-dependent reliability；PHI2 method;irst-order reliability method （FORM）；Kriging model; outcrossing rate；dimension reduction method

結(jié)構(gòu)可靠度分析旨在通過考慮結(jié)構(gòu)或荷載的隨機(jī)性確定結(jié)構(gòu)在規(guī)定時間內(nèi)的失效概率。傳統(tǒng)的結(jié)構(gòu)可靠度分析[1-2]一般認(rèn)為結(jié)構(gòu)的模型特性與服役環(huán)境不隨時間而變化。實際上，結(jié)構(gòu)在服役期間的幾何形狀、材料性能等特性會隨時間而退化，所受荷載亦隨時間而變化[3-6]，對結(jié)構(gòu)進(jìn)行時變可靠度分析才能評估其真實狀態(tài)。目前，基于首次超越破壞準(zhǔn)則的時變可靠度分析方法主要分為兩種[7-9]：一種是基于極值的方法，另一種是基于跨越率的方法。

基于極值的方法基本思路是根據(jù)極限狀態(tài)函數(shù)極值的概率特征獲得失效概率，從而將時變問題轉(zhuǎn)化為時不變問題[10]。蒙特卡洛模擬（MCS）方法是一種典型的基于極值的時變可靠度方法，該方法簡單且精確，但在實際應(yīng)用中需消耗大量計算資源[1]。為提高對可靠度問題的分析效率，研究者們提出了兩種改進(jìn)思路：一種是對采樣方法進(jìn)行改進(jìn)，如引入重要性采樣方法[2]和子集模擬方法[13]等，但這種方法在處理小失效概率事件時仍不可避免地需要大量樣本以達(dá)到較高精度；另一種是近似處理極限狀態(tài)函數(shù)極值的概率分布，如引入概率密度演化方法[14]和代理模型方法[2，5-6]等。其中，Kriging代理模型因具有同時提供預(yù)測點期望和標(biāo)準(zhǔn)差的特性而被廣泛應(yīng)用?？紤]到基于一次性抽樣的普通Kriging方法無法保證訓(xùn)練樣本的有效性，將Kriging模型和高效的主動學(xué)習(xí)策略相結(jié)合的AK類方法[15-17]近年成為研究熱點。將AK方法與一次可靠度方法（FORM）相結(jié)合的AK-FORM方法在保證精度的同時有效提高了FORM方法的效率，是一種具有發(fā)展前景的可靠度分析方法。

基于跨越率的方法基本思路是基于泊松過程、馬爾可夫過程等假設(shè)，利用數(shù)值積分方法將任意離散時刻極限狀態(tài)函數(shù)的跨越率映射為失效概率?？缭铰首钤缬蒖ice[18提出，但其概念難以應(yīng)用于實際。隨后，研究者們通過引入首次超越公式求解方法和基于跨越事件發(fā)生服從泊松過程假設(shè)的解析方法等，拓寬了跨越率模型的應(yīng)用范圍。近年來，PHI2方法19成為了一種被廣泛應(yīng)用的時變可靠度分析方法，該方法將跨越率的計算轉(zhuǎn)換為并聯(lián)的靜態(tài)問題，從而實現(xiàn)了跨越率的高效求解。為達(dá)到足夠的計算精度，PHI2方法會將時間周期以較小的時間增量進(jìn)行離散化，即需要計算大量時刻處的跨越率。顯然，跨越率的求解效率是影響PHI2方法效率的關(guān)鍵因素，而跨越率的求解主要涉及到相關(guān)系數(shù)計算、相鄰時刻處的可靠指標(biāo)計算以及二維積分計算。在處理極限狀態(tài)面很復(fù)雜的實際工程問題時，單次跨越率的求解效率有待提高。為提升PHI2方法的效率，研究者們提出了PHI2+[20]、EPHI2[11]等方法，但這些方法大多僅從某一方面改進(jìn)了跨越率的計算效率，并且有可能獲得不準(zhǔn)確的可靠度評估結(jié)果。

筆者從3個方面改進(jìn)PHI2方法，形成一種高效的時變可靠度分析方法。首先，采用無Cholesky分解的策略來減少可靠指標(biāo)求解過程中的隨機(jī)變量數(shù)目，然后采用基于主動學(xué)習(xí)Kriging模型的改進(jìn)一次可靠度（AK-FORM）方法來計算各時刻的可靠指標(biāo)，再利用降維方法將二維積分轉(zhuǎn)化為一維積分，從而完成跨越率的計算。最后通過數(shù)值算例和工程算例驗證建議方法的精度和效率。

1 時變可靠度分析的PHI2方法

1.1 PHI2方法的原理

假設(shè)結(jié)構(gòu)時變可靠度問題的極限狀態(tài)方程為

Z（t）=G（X，Y（t），t）

式中： 為 n-1 維隨機(jī)變量向量；Y（t） 為一維隨機(jī)過程。

基于跨越事件的方法是計算時變可靠性問題最常用的方法之一，將時間周期 [0，t] 內(nèi)跨越事件的數(shù)量定義為 N⁺（0，t） ，結(jié)構(gòu)在時間周期 [0，t] 內(nèi)的累計失效概率可表示為[19-20]

由式（2）可知， P_f，c（0，t） 的上界為[21]

式中： P_t，i（0） 表示初始時刻結(jié)構(gòu)的失效概率； v⁺（τ）為 τ 時刻的跨越率，可由式（4）計算。

式中： Δτ 為時間增量。

若 Y（t） 在 τ 與 τ+Δτ 時刻的截口隨機(jī)變量分別用 Y_k^（1） 與 Y_k^（2） 表示，兩者間的相關(guān)系數(shù)為

若 Y（t） 為平穩(wěn)非高斯隨機(jī)過程，則 Y_k^（1） 與 均為非高斯變量，可引入Nataf變換將其轉(zhuǎn)換為獨立的標(biāo)準(zhǔn)正態(tài)變量。為簡便，僅考慮 Y（t） 為平穩(wěn)高斯隨機(jī)過程的情形，可先通過線性變換將其轉(zhuǎn)化為均值為0、標(biāo)準(zhǔn)差為1的高斯過程 。相應(yīng)地，Y_k^（1） 與 可表示為相應(yīng)的標(biāo)準(zhǔn)正態(tài)變量 與 的線性函數(shù)，且 與 的相關(guān)系數(shù)仍為ρ_Yk（τ，τ+Δτ） ；然后可引入Cholesky分解，將 與 用獨立的標(biāo)準(zhǔn)正態(tài)變量 表示，即

于是，在標(biāo)準(zhǔn)正態(tài)空間中， τ 時刻的極限狀態(tài)面可由超平面近似為

式中： （204u₃…u_n+1 分別為隨機(jī)向量 X 轉(zhuǎn)換至標(biāo)準(zhǔn)正態(tài)空間后對應(yīng)的標(biāo)準(zhǔn)正態(tài)變量； β（τ） 為 τ 時刻的可靠指標(biāo)。同理， τ+Δτ 時刻的極限狀態(tài)面可近似為

式中： （20α_h，3…α_h，n+1）：β（τ+Δτ） 為 τ+Δτ 時刻的可靠指標(biāo)。在PHI2方法中[19]， τ 和 τ+Δτ 時刻的線性化功能函數(shù)的相關(guān)系數(shù)為

式中： α～（τ）=（α_l，1，0，α_l，3…α_l，n+1） 。此時，跨越率的計算可轉(zhuǎn)換為式（10）所示并聯(lián)靜態(tài)問題的求解。

式中： 均可由FORM方法求解； 為二維標(biāo)準(zhǔn)正態(tài)分布函數(shù)，即

式中： φ₂（?） 為二維標(biāo)準(zhǔn)正態(tài)概率密度函數(shù)。

1. 2 PHI2方法的計算性能

利用PHI2方法進(jìn)行時變可靠度分析的關(guān)鍵在于跨越率 v⁺（τ） 的求解，主要包括3個步驟：首先，通過Cholesky分解將隨機(jī)變量 與 獨立化，并推導(dǎo)相關(guān)系數(shù) ρ_G（τ，τ+Δτ） ；然后，利用FORM方法計算各時刻的可靠指標(biāo)；最后，利用二維數(shù)值積分方法求解 Φ₂[β（τ），-β（τ+Δτ），ρ_G（τ，τ+Δτ）] 顯然，上述3個步驟的計算效率對方法的總效率有著至關(guān)重要的影響。

1.2.1Cholesky分解對計算效率的影響

由式（6）可知，通過引入Cholesky分解， τ+Δτ 時刻的隨機(jī)變量 表示為2個隨機(jī)變量的線性組合，因此，計算 τ+Δτ 時刻的可靠指標(biāo)時，所需考慮的隨機(jī)變量數(shù)目由 n 增加至 n+1 。若采用FORM方法計算可靠指標(biāo)，則每一迭代過程所需計算的梯度值亦增加1個；若采用單邊差分法計算梯度值，相應(yīng)功能函數(shù)的調(diào)用次數(shù)亦增加1次。因此，采用Cholesky分解將隨機(jī)變量 獨立化會導(dǎo)致τ+Δτ 時刻可靠指標(biāo)計算的每次迭代都多調(diào)用1次功能函數(shù)，從而降低可靠指標(biāo)的求解效率。

不難發(fā)現(xiàn)，若能不增加 τ+Δτ 時刻隨機(jī)變量的數(shù)目，將有助于提高該時刻可靠指標(biāo)的計算效率。

1.2.2各時刻可靠指標(biāo)計算效率的影響

為達(dá)到足夠的失效概率預(yù)測精度，PHI2方法通常將時間周期以較小的時間增量進(jìn)行離散化，因此，需要計算大量時刻處的跨越率。由式（10）可知，求解任意時刻的跨越率均涉及到 β（τ） 和 β（τ+ Δτ 的計算。因此，若采用FORM方法計算各時刻的可靠指標(biāo)，PHI2方法將涉及大量FORM方法的運行，顯然，F(xiàn)ORM方法的分析效率會在很大程度上影響PHI2方法的分析效率。然而，F(xiàn)ORM方法的每一迭代過程均需通過調(diào)用功能函數(shù)確定迭代點的功能函數(shù)值與梯度值，若能在迭代過程中減少功能函數(shù)的調(diào)用次數(shù)，時變可靠度分析效率將得到極大提高。

1.2.3二維積分的求解

由式（11）可知， Φ₂[β（τ），-β（τ+Δτ），ρ_G（τ，τ+ Δτ） ]本質(zhì)上是二維積分，盡管可采用數(shù)值積分求解，但并不特別方便。尤其需要指出的是，PHI2方法中， .ρ_G（τ，τ+Δτ） 的值較為接近一1[19]，此時采用二維數(shù)值積分時 φ₂（?） 中的大部分積分節(jié)點為無效節(jié)點[22]，導(dǎo)致二維數(shù)值積分方法效率較低。

若能避免直接采用二維數(shù)值積分求解 Φ₂[β（τ） ，-β（τ+Δτ），ρ_G（τ，τ+Δτ）]， 將有助于改善時變可靠度分析的性能。

2 高效的時變可靠度分析方法

為提高PHI2方法的計算效率，引入3個策略對其進(jìn)行改進(jìn)：

1）不引入Cholesky分解將隨機(jī)變量 與 獨立化，以避免各個時刻可靠指標(biāo)計算時隨機(jī)變量數(shù)目的增加。

2）引人基于主動學(xué)習(xí)Kriging模型的改進(jìn)一次可靠度方法計算各時刻的可靠指標(biāo)，以提高FORM方法的計算效率。

3）利用降維方法將式（11）中的二維積分轉(zhuǎn)化為一維積分，以簡化計算

2.1無Cholesky分解時相關(guān)系數(shù)的計算

若不引入Cholesky分解將隨機(jī)變量 與 獨立化，即 τ 時刻與 τ+Δτ 時刻的 與 均為基本隨機(jī)變量，那么 τ 時刻與 τ+Δτ 時刻的功能函數(shù)均只涉及 n 個隨機(jī)變量，這有助于改善 β（τ+Δτ） 的計算效率，但 ρ_G（τ，τ⁺Δτ） 的表達(dá)式（9）將不再適用。

綜合式（6）式（7）不難發(fā)現(xiàn)， τ 時刻的近似極限狀態(tài)面仍可由式（7）表示，但 τ+Δτ 時刻的近似極限狀態(tài)面需修改為

式中 。很明顯， 中僅包含 n 個隨機(jī)變量。定義標(biāo)準(zhǔn)正態(tài)變量

V_L 與 V_H 間的相關(guān)系數(shù)為

ρ_VLVH=E（V_L?V_H）=

當(dāng) i≠j 時，有 E[u_i?u_j]=0 。因此，式（15）可改寫為

由于式（16）中 ，式（16）可被進(jìn)一步化簡成

由式（7）式（12）可知， .ρ_G（τ，τ⁺Δτ） 為 ρ_VLVH 的相反數(shù)，即

值得指出的是，若將無Cholesky分解的策略直接與PHI2方法相結(jié)合，在計算 β（τ+Δτ） 時僅涉及 n 個隨機(jī)變量，可進(jìn)一步提高PHI2方法的分析效率，進(jìn)而形成一個改進(jìn)的PHI2方法，將此方法記為PHI2方法。

2.2基于主動學(xué)習(xí)Kriging模型的改進(jìn)一次可靠度方法

各時刻可靠指標(biāo)的求解效率對時變可靠度分析的效率十分重要。因此，采用高效的基于主動學(xué)習(xí)Kriging模型的改進(jìn)一次可靠度方法（AK-FORM）[]求解 β（τ） 和 β（τ+Δτ） ，以改善各時刻可靠指標(biāo)的計算效率。

該方法在迭代過程中有兩個階段：全局搜索階段和局部搜索階段。若第 q 次迭代點 u^（q） 與第 q-1 次迭代點 的距離 ||u^（q）-u^（q-1）| |大于等于某一較小的閾值 （可取 c=0.1， ，則屬于全局搜索階段，否則屬于局部搜索階段。在全局搜索階段沿用已有FORM方法迭代至 ∥u^（q）-u^（q-1）∥lt;0.1 。在局部搜索階段，首先以全局搜索階段的迭代點和計算梯度值所用差分點為初始訓(xùn)練點 為初始樣本點數(shù)）建立初始Kriging模型，然后根據(jù)Kriging模型在后續(xù)迭代點處的預(yù)測精確性自適應(yīng)地更新模型。假設(shè)第l次更新的Kriging模型為 ，由 確定的迭代點記作 。該迭代點的預(yù)測精度可由 S（u^（M+l）） 衡量。

式中： 和 分別為 對應(yīng)的預(yù)測值和標(biāo)準(zhǔn)差。由Kriging模型的高斯特性可知，S（u^（M+l））≥50 時， 對應(yīng)的預(yù)測響應(yīng) 的誤差小于 6% 的概率大于 99.7% ，說明 具有較高預(yù)測精度，不需要更新 。反之，需要將 加人訓(xùn)練點集 并結(jié)合該點真實的響應(yīng) z^（k+l） 更新 ，即

式中 ：f（U） 為多項式基函數(shù)； ：P^（l+1）?m（U）^（l+1）"分別為利用訓(xùn)練點""及其真實的響應(yīng)建立的回歸系數(shù)向量和方差為 σ²"的零均值高斯過程。在局部搜索階段的每一迭代過程中，均采用更新的Kriging模型式（20）代替功能函數(shù)進(jìn)行計算。正是由于在局部搜索階段引入了自適應(yīng)的更新Kriging模型，AK-FORM方法通常較常規(guī)的FORM方法具有更高的計算效率。AK-FORM方法的流程圖見文獻(xiàn)[17]。

2.3二維標(biāo)準(zhǔn)正態(tài)分布函數(shù)的降維積分計算方法

為避免直接采用二維數(shù)值積分求解 Φ₂[β（τ） -β（τ+Δτ），ρ_G（τ，τ+Δτ）] ，引入降維方法23將二維只分簡化為一維積分，即

式中： 表示一維標(biāo)準(zhǔn)正態(tài)分布函數(shù)。將式（21）代入式（10），此時僅用一維數(shù)值積分方法即可求解τ 時刻的跨越率，相較于二維積分的求解更加簡單高效，并且解決了 ρ_G（τ，τ+Δτ） 取值接近一1所導(dǎo)致的積分節(jié)點利用率低的問題，從而進(jìn)一步提升了PHI2方法的時變可靠度分析效率。需要指出的是，體系可靠度分析研究中，針對 Φ₂[β（τ），-β（τ+Δτ） ρG（τ，τ十△τ）]給出了一些簡化的近似計算方法[24-26]，盡管由于不涉及積分而計算簡單，但其適用范圍與精度往往受到一定影響。

2.4時變可靠度分析及步驟

將上述3個策略與PHI2方法相結(jié)合，即形成了建議的基于AK-FORM方法和降維方法的高效時變可靠度分析方法，簡記為K-PHI2方法。該方法流程圖見圖1，其主要步驟包括

1）將時間段 [0，t] 以dt為間隔進(jìn)行離散，并確定時間增量 Δτ ，使得 ρ_Yk（τ，τ+Δτ） 的取值位于[0.990，0.995]之間。

2）利用AK-FORM方法計算 τ 和 τ+Δτ 時刻的可靠指標(biāo) β（τ），β（τ+Δτ） 和靈敏度系數(shù) Δτ） 。

3）參考無Cholesky分解策略，由式（18）計算 τ 和 τ+Δτ 時刻間的極限狀態(tài)面相關(guān)系數(shù) ρ_G（τ，τ⁺ Δτ， 。

4）將步驟2）和3）中得到的 β（τ），β（τ+Δτ） 和 ρ_G （τ，τ+Δτ） 代入式（21），然后利用一維數(shù)值積分方法求解

5）將 Φ₂[β（τ），-β（τ+Δτ），.._G（τ，τ+Δτ） ]代入式（10），確定 τ 時刻的跨越率 v⁺（τ） 。

6）不斷重復(fù)步驟2） ～5 ）計算新時間 τ+dt 的跨越率，直至得到所有時刻的跨越率。

7）由式（3）計算失效概率 P_f，c（0，t） 。

3 算例分析

首先通過一個設(shè)計的簡單數(shù)值算例說明K-PHI2方法的計算流程，然后用一個經(jīng)典的具有顯式表達(dá)式的工程算例對比各方法的計算性能，最后將驗證后的方法應(yīng)用到了涉及有限元分析的工程算例中，進(jìn)而表明K-PHI2對于各種時變可靠度問題的適用性。各方法的效率以功能函數(shù)調(diào)用次數(shù) N 來評估，精度由可靠指標(biāo)的相對誤差 εr_β 來評估，即

式中： β_mcs 為MCS方法的可靠指標(biāo)計算結(jié)果； β 為其余方法（即PHI2方法[19]、PHI2-方法、 PHI2+ 方法[20]以及K-PHI2方法）的可靠指標(biāo)計算結(jié)果。其中，PHI2 + 方法基于有限差分法對PHI2方法進(jìn)行了改進(jìn)，其跨越率表達(dá)式為

式中： φ（?） 為一維標(biāo)準(zhǔn)正態(tài)概率密度函數(shù)。

3.1 算例1：數(shù)值算例

考察由文獻(xiàn)[17]修改而來的功能函數(shù)

式中：退化系數(shù) r=0.009 ；隨機(jī)變量向量 R 和 G 分別為構(gòu)件的初始抗力和永久荷載效應(yīng)，各隨機(jī)變量相互獨立，其概率信息見表1； Y（t） 為高斯過程，在此為時變的可變荷載效應(yīng) Q（t） ，其間隔時間為 Δτ 的兩截口隨機(jī)變量間的相關(guān)系數(shù)為

式中： λ 為相關(guān)長度，取 λ=1/12a ，時間間隔 Δτ= 0.1λ，此時 ρ_Yk（τ，τ+Δτ）=0.99 。各隨機(jī)變量和隨機(jī)過程的統(tǒng)計信息見表1。

該算例對構(gòu)件5a內(nèi)的時變可靠度進(jìn)行評估，時間段[0，5]被均勻分成125個區(qū)間。由于引入了無Cholesky分解策略，K-PHI2方法每個時刻的跨越率計算僅涉及3個隨機(jī)變量。以 τ=1 時的跨越率v⁺（1） 的求解為例，簡要說明K-PHI2方法的計算流程：首先，利用AK-FORM方法計算可靠指標(biāo) 和靈敏度系數(shù)α（1）=（-0.8158415，0.1039628，0.5688535） 、α（1+Δτ）=（-0.8158360，0.1039643，0.5688611） ，這一過程僅需調(diào)用36次功能函數(shù)。相比較而言，PHI2-方法計算可靠指標(biāo)和靈敏度系數(shù)也僅涉及3個隨機(jī)變量，但由于采用了傳統(tǒng)FORM方法，其功能函數(shù)調(diào)用次數(shù)為88。由于需要Cholesky分解，PHI2方法計算可靠指標(biāo)和靈敏度系數(shù)會涉及4個隨機(jī)變量。此外，PHI2方法還采用傳統(tǒng)FORM方法，其功能函數(shù)調(diào)用次數(shù)為99。然后，K-PHI2方法直接由式（18）得到 τ 和 τ+Δτ 時刻的極限狀態(tài)面相關(guān)系數(shù) ρ_G（1，1+Δτ）=-0.9967801 。接下來，將 與 ρ_G（1，1+Δτ） 帶人由降維策略得到的式（21），即可直接利用一維數(shù)值積分方法求解 10^-6 。最后，將求解結(jié)果代入式（10）即得到 τ 時刻的跨越率 。顯然，K-PHI2在3個方面對PHI2方法進(jìn)行了改進(jìn)，各方法的計算結(jié)果見表2和圖2、圖3。

表2記錄了各方法計算出的第5年時結(jié)構(gòu)的失效概率以及總的功能函數(shù)調(diào)用次數(shù)。圖2、圖3分別展示了各方法計算出的5a內(nèi)結(jié)構(gòu)的時變可靠指標(biāo)和失效概率。由圖可知， PHI2+ 方法精度較差，而提出的PHI2-、K-PHI2方法與PHI2方法精度相當(dāng)，這與表1中的結(jié)果相吻合，表明PHI2-、K-PHI2方法與PHI2方法一樣具有良好精度。此外，相較于PHI2方法，由于不引入Cholesky分解，PHI2-方法減少了 τ+Δτ 時刻的基本隨機(jī)變量數(shù)目，功能函數(shù)的調(diào)用次數(shù)就由12375次減少為11000次，說明建議的不引入Cholesky分解策略有助于提升可靠度分析效率。而在PHI2-方法的基礎(chǔ)上，K-PHI2方法采用更高效的AK-FORM方法代替FORM方法計算 β（τ） 和 β（τ+Δτ） ，從而使功能函數(shù)調(diào)用次數(shù)由11000次降低至4590次，說明K-PHI2方法能進(jìn)一步提升可靠度分析效率。需要指出的是，上述各方法的效率均高于MCS方法。綜上所述，PHI2-方法在一定程度上提高了PHI2方法的效率，K-PHI2方法在其基礎(chǔ)上做出進(jìn)一步改善，在保證精度的同時極大地提高了時變可靠度分析效率。

3.2算例2：工程算例1一顯式表達(dá)式算例

考察如圖4所示的矩形截面簡支腐蝕梁的可靠度問題[1，19-20]，功能函數(shù)為

a（t）=a₀-2rt;b（t）=b₀-2rt

式中： 分別為材料屈服應(yīng)力、初始梁寬和梁，各隨機(jī)變量相互獨立，其概率信息見表 3;Y（t） 為高斯過程，在此為時變的豎向集中荷載 F（t） ，其相關(guān)長度和時間間隔同算例1。 L=5m 為梁的跨長； c=78.5kN/m³ 為鋼密度； r=0.03 mm/a 為腐蝕速率。

表4統(tǒng)計了各方法對該算例的計算結(jié)果，圖5、圖6為各方法所得可靠指標(biāo)和失效概率。結(jié)果表明，PHI2 + 方法精度欠佳；其余方法與MCS方法之間的可靠指標(biāo)相對誤差均在 5% 以內(nèi)。其中，PHI2方法總共需要調(diào)用68150次功能函數(shù)，而PHI2-方法和K-PHI2方法分別僅需調(diào)用61955次和23391次功能函數(shù)。說明在保證精度的前提下，PHI2-方法的效率稍高于PHI2方法和PHI2 + 方法，而K-PHI2方法的效率遠(yuǎn)高于其他各方法。實際工程中存在許多與此算例類似的問題，而建議方法適用于分析此類問題，且相較于PHI2方法和 PHI2+ 方法性能更優(yōu)，能高效準(zhǔn)確地對結(jié)構(gòu)可靠度進(jìn)行分析。

3.3算例3：工程算例2一有限元算例

圖7所示桁架結(jié)構(gòu)27的時變可靠度對應(yīng)的功能函數(shù)為

式中： 分別為桿的橫截面積和彈性模量初始值，各隨機(jī)變量相互獨立，其概率信息見表5； Y（t） 為高斯過程，在此為時變的豎向集中荷載P（t） ，其相關(guān)長度和時間間隔同算例1。 s（P，A，E（t）） 為節(jié)點 O 處的撓度 V ，由有限元分析計算得出，其中， 為退化系數(shù)； s₀ 為s（P，A，E（t） ）的閾值，取 0.105m 。

表6、圖8和圖9中給出了各方法的計算結(jié)果。由于 PHI2+ 會出現(xiàn)相對誤差 gt;100% 的情況，圖表中僅列出了PHI2、PHI2-和K-PHI2方法的計算結(jié)果。結(jié)果表明，K-PHI2方法僅需調(diào)用5016次功能函數(shù)就能達(dá)到與PHI2相當(dāng)?shù)木?，相較于其他各方法效率更高。盡管此算例的功能函數(shù)為隱式函數(shù)且其響應(yīng)涉及有限元求解，但計算結(jié)果表明，K-PHI2方法仍然具有較好的計算性能，反映了其處理實際工程問題的適用性。

4結(jié)論

從3個方面改善各時刻跨越率的計算效率，從而形成一種高效的時變可靠度分析方法（K-PHI2方法），得到如下結(jié)論：

1）PHI2-方法能在一定程度上提高PHI2方法的效率，即不引入Cholesky分解將隨機(jī)變量 與 獨立化，能避免 τ+Δτ 時刻隨機(jī)變量數(shù)目的增多，從而在一定程度上提高時變可靠度分析效率。

2）引入AK-FORM方法計算 τ 和 τ+Δτ 時刻的可靠指標(biāo) 和靈敏度系數(shù) α（τ），α（τ+ Δτ） ，明顯減少了這一過程中功能函數(shù)的調(diào)用次數(shù)；此外，引入降維方法將 τ+Δτ） ]的求解簡化為一維積分問題，進(jìn)一步提高了時變可靠度分析效率。

3算例結(jié)果表明，提出的K-PHI2方法同時適用于數(shù)值算例和工程算例，在保證精確度的情況下，其效率明顯高于PHI2方法、PHI2-方法和PHI2 + 方法。

參考文獻(xiàn)

[1]YANG MD，ZHANGD Q，HANX.New efficient androbust method for structural reliabilityanalysis and itsapplication in reliability-based design optimization [J]. Computer Methods in Applied Mechanics and Engineer ing，2020，366：113018.

[2］范文亮，劉丞，李正良.基于HLRF法與修正對稱秩1 方法的改進(jìn)可靠度方法[J].工程力學(xué)，2022，39（9）：1-9. FANWL，LIUC，LIZL.Improvedreliabilitymethod based on HLRF and modified symmetric rank1 method[J]. EngineeringMechanics，2022，39（9）： 1-9. （in Chinese）

[3]JIANGC，WEI XP，WUB，et al.An improved TRPD method for time-variant reliability analysis [J]. Structural and Multidisciplinary Optimization，2O18，58 （5）：1935-1946.

[4]ZAFART，WANGZL.An efficient method for timedependent reliability prediction using domain adaptation [J].Structural and Multidisciplinary Optimization，2020， 62（5）：2323-2340.

[5]JIANGC，HUZ，LIUYX，etal.Asequential calibration andvalidationframeworkformodeluncertainty quantification and reduction [J]. Computer Methods in Applied Mechanics and Engineering，202O，368： 113172.

[6]WU JH， ZHANG D Q，JIANG C，et al. On reliability analysis method through rotational sparse grid nodes [J]. Mechanical Systems and Signal Processing，2O21，147： 107106.

[7]HAWCHAR L，EL SOUEIDY C P，SCHOEFS F. Principal component analysis and polynomial chaos expansion for time-variant reliability problems[J].Reliability Engineeringamp; System Safety，2O17，167： 406-416.

[8]YU S，ZHANG Y W，LI Y，et al.Time-variant reliability analysis via approximation of the first-crossing PDF [J]. Structural and Multidisciplinary Optimization， 2020，62（5）： 2653-2667.

［9］唐和生，郭雪媛，薛松濤.基于廣義子集模擬和自適應(yīng) Kriging模型的非線性隨機(jī)動力系統(tǒng)的時變可靠性分析 [J].振動與沖擊，2021，40（21）：47-54. TANG H S，GUO X Y，XUE S T. Time-varying reliabilityanalysis of nonlinear stochastic dynamic systems based on generalized subset simulation and adaptive Kriging model[J]. Journal of Vibration and Shock，2O21，40 （21）： 47-54. （in Chinese）

[10]HU Z， DU X P. A sampling approach to extreme value distribution for time-dependent reliability analysis [J]. Journal of Mechanical Design，2013，135（7）： 071003.

[11]LIXW，ZHAO YG，ZHANGXY，et al.Explicit model of outcrossing rate for time-variant reliability [J]. ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems，Part A： Civil Engineering，2O22，8（1）： 04021087.

[12］羅立勝，陳志華.銹蝕鋼構(gòu)件時變可靠度的重要抽樣法 研究[J].工業(yè)建筑，2019，49（8）：1-5，107. LUO L S，CHEN Z H. Time-dependent reliability analysis of corroded steel members based on the important sampling method [J]. Industrial Construction，2O19，49 （8）：1-5，107.（in Chinese）

[13] WANG D P，QIU H B，GAO L， et al. A single-loop Kriging coupled with subset simulation for time-dependent reliability analysis [J].Reliability Engineeringamp;. System Safety，2021，216：107931.

[14］王思文，王憲杰，胡彪，等.基于概率密度演化理論的 LRB隔震支座時變可靠度研究[J].計算力學(xué)學(xué)報， 2020，37（6）： 776-782. WANG SW，WANG XJ，HUB，et al. Study on timevarying reliability of LRB isolated bearings based on probability density evolution method [J]. Chinese Journal of Computational Mechanics，2020，37（6）：776-782. （in Chinese）

[15]WANG YJ，PANH，SHI YN，et al. A new activelearning estimation method for the failure probability of structural reliability based on Kriging model and simple penalty function [J]. Computer Methods in Applied Mechanics and Engineering，2023，41O： 116035.

[16] HONG L X. SHANG R. IIS 7. et al. Portfolin alloca tion strategy for active learning Kriging-based structural reliability analysis[J]. Computer Methods in Applied Me chanics and Engineering，2023，412： 116066.

[17]劉丞，范文亮，余書君，等.基于主動學(xué)習(xí)Kriging模型 的改進(jìn)一次可靠度方法[J].工程力學(xué)，2024，41（2）： 35-42. LIU C，F(xiàn)AN WL，YU S J，et al. Improved first order reliability method based on adaptive kriging model online first[J]. Engineering Mechanics，2024，41（2）： 35-42. （in Chinese）

[18] RICE S O.Mathematical analysis of random noise [J]. The Bell System Technical Journal，1944，23（3）： 282-332.

[19]ANDRIEU-RENAUD C，SUDRET B，LEMAIRE M. The PHI₂ method： A way to compute time-variant reliability [J].Reliability Engineeringamp;.System Safety， 2004，84（1）： 75-86.

[20] SUDRET B. Analytical derivation of the outcrossing rate in time-variant reliability problems [J]. Structure and Infrastructure Engineering，2008，4（5）： 353-362.

[21] SHINOZUKA M. Probability of structural failure under random loading[J]. Journal of the Engineering Mechanics Division，1964， 90（5）： 147-170.

[22］范文亮，楊朋超，李正良.基于Mehler公式的等效相關(guān) 系數(shù)求解技術(shù)[J].同濟(jì)大學(xué)學(xué)報（自然科學(xué)版），2016， 44（6）： 838-844. FAN W L，YANG P C，LI Z L. A technique for solution of equivalent correlation coefficients based on the Mehler’s formula [J]. Journal of Tongji University （Natural Science），2016，44（6）：838-844. （in Chinese）

[23]貢金鑫，趙國藩.二維正態(tài)分布函數(shù)值的一個近似算法 [J].計算結(jié)構(gòu)力學(xué)及其應(yīng)用，1996，13（4）：494-499. GONG JX，ZHAO GF. An approximate algorithm for bivariate normal integral [J]. Chinese Journal of ComputationalMechanics，1996，13（4）： 494-499.（in Chinese）

[24]FENG Y S.A method for computing structural system reliability with high accuracy [J]. Computers amp;. Structures，1989，33（1）： 1-5.

[25]董聰，酈正能.結(jié)構(gòu)系統(tǒng)可靠性精確計算理論[J].強(qiáng)度 與環(huán)境，1995，22（3）：46-51. DONG C，LI Z N. An exact evaluation theory of structural system reliability [J]. Structure amp;.Environment Engineering，1995，22（3）： 46-51. （in Chinese）

[26]姚繼濤，趙國藩，浦聿修.二維標(biāo)準(zhǔn)正態(tài)聯(lián)合概率的計 算[J].建筑結(jié)構(gòu)學(xué)報，1996，17（4）：10-19. YAO JT，ZHAO G F，PU Y X. Calculation of joint probability of two-dimensional standard normal [J]. Journal of Building Structures，1996，17（4）：10-19. （in Chinese）

[27] SUN Z L，WANG J，LIR，et al. LIF： A new Kriging based learning function and its application to structural reliability analysis [J]. Reliability Engineeringamp;. System Safety，2017，157：152-165.

（編輯胡玲）