顧 倩，夏 恒，何 軍

(上海交通大學(xué)船舶海洋與建筑工程學(xué)院，上海 200240)

生命線工程系統(tǒng)是對維系城市經(jīng)濟(jì)發(fā)展和社會生活具有重大影響的基礎(chǔ)性工程設(shè)施系統(tǒng)，如交通、通訊、供水、供電、煤氣和輸油等工程系統(tǒng)。根據(jù)網(wǎng)絡(luò)理論，生命線工程系統(tǒng)的網(wǎng)絡(luò)模型由四個要素構(gòu)成，即節(jié)點、線路、流量和賦權(quán)形式。而根據(jù)網(wǎng)絡(luò)賦權(quán)形式的不同，網(wǎng)絡(luò)系統(tǒng)可分為邊權(quán)網(wǎng)絡(luò)(只考慮邊型生命線構(gòu)件或結(jié)構(gòu)的地震失效)、點權(quán)網(wǎng)絡(luò)(只考慮節(jié)點型生命線構(gòu)件或結(jié)構(gòu)的地震失效)和一般賦權(quán)網(wǎng)絡(luò)(同時考慮邊型和節(jié)點型生命線構(gòu)件或結(jié)構(gòu)的地震失效)。

在生命線地震工程中，將系統(tǒng)的抗震可靠度定義為地震作用下生命線系統(tǒng)拓?fù)渚W(wǎng)絡(luò)中源點到終點的連通可靠度。一般的生命線工程系統(tǒng)往往具有多個源點和多個終點，如上海供水管網(wǎng)系統(tǒng)包含多個水廠和多個用戶端。對于這一問題，可通過設(shè)置一個單向聯(lián)結(jié)各源點或各終點的虛擬節(jié)點的辦法將其轉(zhuǎn)化為兩終端問題[1]。此時，兩虛擬節(jié)點的連通可靠度即為原多源多匯系統(tǒng)的連通可靠度。為了簡化分析過程，本文僅考慮單源點到單終點的兩終端生命線系統(tǒng)網(wǎng)絡(luò)連通可靠度問題。生命線工程網(wǎng)絡(luò)模型規(guī)模龐大、類型多樣，導(dǎo)致其具有很大的復(fù)雜性。另外，地震作用下生命線工程還具有失效相關(guān)性，包括兩個方面：第一，地震動的空間相關(guān)性，即同一地震輸入對不同元件的破壞是具有相關(guān)性的；第二，系統(tǒng)結(jié)構(gòu)形式導(dǎo)致的失效相關(guān)性，網(wǎng)絡(luò)系統(tǒng)中節(jié)點的失效會引起此節(jié)點相聯(lián)邊的失效，反之，與一節(jié)點相聯(lián)的所有邊的失效也相當(dāng)于此節(jié)點的失效。因此，要進(jìn)行生命線工程網(wǎng)絡(luò)系統(tǒng)抗震可靠度估計的研究，難點就在于其計算復(fù)雜性和失效相關(guān)性。

一方面，網(wǎng)絡(luò)系統(tǒng)可靠度的計算問題是典型的NP-Hard 問題，計算量呈現(xiàn)非多項式增長的趨勢。在過去的幾十年中，國內(nèi)外學(xué)者圍繞著如何降低網(wǎng)絡(luò)系統(tǒng)可靠度的計算量問題展開了一系列的研究，由此產(chǎn)生的方法主要可分為隨機(jī)模擬法和概率解析法。早期的Monte Carlo 隨機(jī)模擬方法，在大型失效獨立網(wǎng)絡(luò)可靠性分析中，往往能取得較好的效果。但是，這一算法在計算工作量和算法收斂性方面存在一定的局限性[2?3]。概率解析算法是通過應(yīng)用網(wǎng)絡(luò)分析技術(shù)，建立精確或近似的數(shù)學(xué)模型。在概率解析算法中，較為成熟的方法有不交積之和算法[4?5]、因子分解定理算法[6?8]、二分決策圖算法[9?10]以及最小路遞推分解算法[1]。這些算法在應(yīng)用中有各自的適用范圍，且求解效率不盡相同。本文的研究將以遞推分解思想為基礎(chǔ)。該思想最早于1979 年由Dotson 等[11]提出，對降低網(wǎng)絡(luò)算法的空間復(fù)雜性研究起到了巨大的推動作用。后續(xù)又經(jīng)過了一系列的發(fā)展，2002年，李杰和何軍[1,12?13]正式提出了遞推分解算法(RDA)，有效地降低了網(wǎng)絡(luò)算法的復(fù)雜性，并且能夠用于大型網(wǎng)絡(luò)的可靠度估計。RDA 算法的詳細(xì)介紹，將在本文第1 節(jié)中給出。

另一方面，系統(tǒng)中元件失效具有相關(guān)性。在生命線工程系統(tǒng)抗震可靠度的研究中，要考慮元件失效相關(guān)的情況，需要解決的是不交最小路(或不交最小割)事件的聯(lián)合失效概率的計算問題。解決這類問題的思路是，在已知基本隨機(jī)變量數(shù)值特征的基礎(chǔ)上，結(jié)合結(jié)構(gòu)功能函數(shù)，計算相關(guān)失效系統(tǒng)的可靠度。Copula 是建立多維隨機(jī)變量聯(lián)合分布的有效方法[14]，Gumbel Copula 就是工程中常用的一種阿基米德型Copula。1981 年，Ang[15]提出了概率網(wǎng)絡(luò)估算技術(shù)(PNET 方法)，該方法將所有主要失效模式按彼此相關(guān)系數(shù)的大小分成若干組，在每一組中選出一個失效概率最大的失效模式作為代表模式，然后假定各代表失效模式相互獨立，再來估算網(wǎng)絡(luò)系統(tǒng)的失效概率。2000 年，Philippe[16]提出串聯(lián)系統(tǒng)可靠度的一階矩方法，該方法采用Cholesky 分解方法把結(jié)構(gòu)功能函數(shù)的基本隨機(jī)變量進(jìn)行分解，結(jié)果表明Maincon 算法的分析結(jié)果精度較高。2013 年，Won-Hee Kang 等[17]提出了SCM 算法，該算法是一種快速可靠性評價算法，主要由“與”合并和“或”合并兩種運算組成，通過這兩種運算將網(wǎng)絡(luò)化簡直至合并為一個點[18]，從而由最終化簡得到的節(jié)點可靠度估計系統(tǒng)可靠度。

本文基于原RDA 算法，先對原來分別按照邊權(quán)、點權(quán)和一般賦權(quán)，需要運用對應(yīng)程序去計算系統(tǒng)可靠度的情況進(jìn)行整合簡化，發(fā)展出同時適用于三種賦權(quán)形式的統(tǒng)一RDA 算法。進(jìn)一步在統(tǒng)一RDA 算法的基礎(chǔ)上，引入SCM 和Gumbel Copula函數(shù)兩種方法分別計算RDA 算法中不交最小路和不交最小割的聯(lián)合失效概率，發(fā)展出相依失效生命線工程系統(tǒng)抗震可靠度估計的統(tǒng)一RDA 算法。最后，以一假想電力系統(tǒng)為例，將其簡化為一36 節(jié)點格柵型拓?fù)渚W(wǎng)絡(luò)，估計該網(wǎng)絡(luò)系統(tǒng)相依失效的動力抗震可靠度，從而驗證本文方法的有效性。

1 統(tǒng)一RDA 算法

本文以文獻(xiàn)[1]提出的RDA 算法為基礎(chǔ)，該算法通過引入結(jié)構(gòu)函數(shù)遞推分解的概念，采用實時不交化的研究方法，能夠有效降低網(wǎng)絡(luò)算法的復(fù)雜性。然而，對于邊權(quán)、點權(quán)和一般賦權(quán)三種不同賦權(quán)形式的網(wǎng)絡(luò)系統(tǒng)，RDA 算法的推導(dǎo)過程有所區(qū)別，實際計算時也需要對不同形式的網(wǎng)絡(luò)系統(tǒng)采用相應(yīng)的算法程序，過程稍顯繁瑣。本章先對原RDA 算法的推導(dǎo)過程進(jìn)行簡要介紹，在此基礎(chǔ)上，將其三部分程序整合化簡，發(fā)展出統(tǒng)一RDA 算法。最后，通過簡單算例，驗證統(tǒng)一RDA算法的正確性。

1.1 原RDA 算法簡介

一般來說，生命線工程系統(tǒng)的網(wǎng)絡(luò)分析模型是一個兩狀態(tài)關(guān)聯(lián)系統(tǒng)，即元件和系統(tǒng)只有安全和失效兩個狀態(tài)，用系統(tǒng)的結(jié)構(gòu)函數(shù)φ(G)來表示，φ(G)=1 表示系統(tǒng)處于安全狀態(tài)，φ(G)=0則表示系統(tǒng)處于失效狀態(tài)。

1.1.1 邊權(quán)網(wǎng)絡(luò)系統(tǒng)

定義系統(tǒng)G=(V,E)中從一起始點到一終止點的所有最小路中路徑最短的一條為系統(tǒng)G=(V,E)安全的一個基本事件A1=a1,1a1,2···a1,m1，則系統(tǒng)的結(jié)構(gòu)函數(shù)：

式中：Ak為系統(tǒng)的第k條最小路；K為系統(tǒng)的所有最小路數(shù)；a1,m1為A1中第m1條邊。

由布爾運算有：

由互斥和公式：

將式(3)代入式(2)中，由布爾簡化得到：

式中，Gi為從網(wǎng)絡(luò)系統(tǒng)G=(V,E) 中去掉邊a1,i后得到的網(wǎng)絡(luò)系統(tǒng)圖，即由a1,1a1,2···1,i決定的系統(tǒng)子圖。

繼續(xù)對各連通子圖分別尋找最短最小路，并按上述格式分解，運用布爾運算法則進(jìn)行化簡與歸并，則在最終分解完畢，不再存在連通子圖時，有：

式中：ci為第i個連通子圖結(jié)構(gòu)函數(shù)前的系數(shù)；Ai為 第i個 連 通 子 圖 的 一個最短最小路；mnc為系統(tǒng)所有連通子圖數(shù)；Li(i=2,3,···,mnc)為歸并運算后的不交最小路。

若令L1=A1，則式(5)可以寫為更為一般的形式：

1.1.2 點權(quán)網(wǎng)絡(luò)系統(tǒng)

點權(quán)網(wǎng)絡(luò)系統(tǒng)的RDA 算法與邊權(quán)網(wǎng)絡(luò)系統(tǒng)的推導(dǎo)過程大致相同，但在最小路的生成和遞推生成子系統(tǒng)圖時有所區(qū)別。RDA 算法采用寬度優(yōu)先搜索生成系統(tǒng)或各子系統(tǒng)的最小路，此時的最小路為頂點最小路。在點權(quán)網(wǎng)路系統(tǒng)中，此時寬度優(yōu)先搜索生成的最小路就是真實的點最小路，不需要再對其進(jìn)行吸收歸并運算。

遞推分解生成點的不交最小路后，點權(quán)網(wǎng)絡(luò)系統(tǒng)可靠度為：

因此，點權(quán)網(wǎng)絡(luò)系統(tǒng)除了網(wǎng)絡(luò)系統(tǒng)可靠度和網(wǎng)絡(luò)失效概率表達(dá)式與邊權(quán)有所不同，其余均與邊權(quán)網(wǎng)絡(luò)系統(tǒng)一致。

1.1.3 一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)

一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)的RDA 算法比邊權(quán)和點權(quán)更為復(fù)雜。在一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)中，頂點和邊可以同時存在失效狀態(tài)。網(wǎng)絡(luò)系統(tǒng)頂點的失效，將導(dǎo)致與此頂點相聯(lián)的邊失去連接功能，實際上相當(dāng)于此頂點的相連邊也發(fā)生了失效；反過來，與一頂點相連的所有邊的失效也相當(dāng)于此頂點的失效。反映在不交最小路和不交最小割上，頂點和邊的狀態(tài)為相依隨機(jī)事件,不交最小路和不交最小割事件的發(fā)生概率不能再采用元件失效獨立假定下的計算公式，而是需要計算相依隨機(jī)事件的聯(lián)合失效概率。

一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)的RDA 算法采用了Torrieri[19]在1994 年提出的一種計算一般賦權(quán)有向網(wǎng)絡(luò)系統(tǒng)不交最小路發(fā)生概率的方法，該方法把頂點可靠度嵌入邊中，使一般賦權(quán)有向網(wǎng)絡(luò)系統(tǒng)等價變換為獨立失效的邊權(quán)網(wǎng)絡(luò)系統(tǒng)，通過計算獨立失效邊權(quán)有向網(wǎng)絡(luò)系統(tǒng)不交邊最小路的發(fā)生概率，得到一般賦權(quán)有向網(wǎng)絡(luò)系統(tǒng)不交最小路的發(fā)生概率，從而避免了聯(lián)合失效概率的計算。

因此，對于一個一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)，先采用獨立邊權(quán)網(wǎng)絡(luò)系統(tǒng)的遞推分解算法找到所有的不交最小路和不交最小割，再對于每一條不交最小路和不交最小割，用Torrieri 提出的方法計算發(fā)生概率，最后得到系統(tǒng)可靠度，具體過程如下。

1.2 統(tǒng)一RDA 算法

原RDA 算法的推導(dǎo)過程表明：系統(tǒng)的結(jié)構(gòu)函數(shù)不必通過求出系統(tǒng)所有最小路來建立，而是可以通過系統(tǒng)的一個最小路分解，來遞推下一個互斥最小路，并且由過程中產(chǎn)生的不連通子圖，可以給出互補(bǔ)結(jié)構(gòu)函數(shù)，從而利用概率不等式計算給出的系統(tǒng)可靠度的上、下界，實時估計系統(tǒng)連通可靠度的近似值。進(jìn)一步地，結(jié)合計算機(jī)編程，可以大幅降低大型網(wǎng)絡(luò)系統(tǒng)可靠度計算的復(fù)雜性。

生命線工程網(wǎng)絡(luò)系統(tǒng)存在單向邊和雙向邊，區(qū)別在于該邊的兩個頂點是否有向。在網(wǎng)絡(luò)圖建模時，以鄰接向量矩陣作為網(wǎng)絡(luò)圖的存貯方式，將雙向邊拆分成兩條單向邊，此時單向邊或雙向邊就可以在鄰接向量矩陣中進(jìn)行區(qū)分。

由于三種賦權(quán)形式的網(wǎng)絡(luò)系統(tǒng)的RDA 算法推導(dǎo)過程有所不同，原RDA 算法的計算程序也分成三部分，在計算時，三種賦權(quán)形式的網(wǎng)絡(luò)系統(tǒng)需要輸入的條件也不同：邊權(quán)網(wǎng)絡(luò)系統(tǒng)由于頂點永遠(yuǎn)安全，默認(rèn)各頂點失效概率為0，只需要輸入各邊可靠度以及網(wǎng)絡(luò)的鄰接向量矩陣；同理，點權(quán)網(wǎng)絡(luò)系統(tǒng)默認(rèn)各邊失效概率為0，只需要輸入各節(jié)點可靠度以及網(wǎng)絡(luò)的鄰接向量矩陣；一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)則需要輸入各邊可靠度、各節(jié)點可靠度以及網(wǎng)絡(luò)的鄰接向量矩陣。

為了簡化實際運用中的計算過程，希望將三部分程序合為一個。根據(jù)三種賦權(quán)形式網(wǎng)絡(luò)的特點，本文在原RDA 算法的基礎(chǔ)上，發(fā)展出統(tǒng)一RDA 算法。

以一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)的RDA 算法程序為載體，將邊權(quán)和點權(quán)的計算程序歸并入一般賦權(quán)程序中，操作方法如下：在計算邊權(quán)網(wǎng)絡(luò)系統(tǒng)可靠度時，采用一般賦權(quán)程序，其中各邊可靠度和網(wǎng)絡(luò)鄰接向量矩陣的輸入與邊權(quán)的輸入相同，各節(jié)點可靠度均取為1，表示此時各節(jié)點是絕對安全的，這樣也符合邊權(quán)系統(tǒng)的定義。同理，用一般賦權(quán)程序計算點權(quán)系統(tǒng)時，各節(jié)點可靠度和網(wǎng)絡(luò)鄰接向量矩陣的輸入與點權(quán)的輸入相同，各邊可靠度均取為1。

1.3 算例驗證

為了驗證統(tǒng)一RDA 算法的正確性，采用下圖1所示網(wǎng)絡(luò)，分別計算邊權(quán)程序、點權(quán)程序和一般賦權(quán)程序退化的邊權(quán)、點權(quán)系統(tǒng)的可靠度結(jié)果，見表1，對比驗證。

圖1 17 節(jié)點網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.1 17-node network structure

表1 17 節(jié)點網(wǎng)絡(luò)可靠度計算結(jié)果Table 1 Reliability calculation results of 17-node network system

該系統(tǒng)有17 個節(jié)點，32 條邊，系統(tǒng)的源點和終點分別為1 和17，邊或節(jié)點的可靠度均取為0.9。

程序計算結(jié)果表明：由統(tǒng)一RDA 算法退化的邊權(quán)、點權(quán)與原邊權(quán)、點權(quán)算得的可靠度結(jié)果完全一致，但不交最小路數(shù)與不交最小割數(shù)不同，而與一般賦權(quán)一致。其中，統(tǒng)一RDA 算法退化的點權(quán)算得的不交最小路數(shù)與不交最小割數(shù)與原點權(quán)算法相差較大，而統(tǒng)一RDA 算法退化的邊權(quán)算得的不交最小路數(shù)與不交最小割數(shù)與原邊權(quán)算法相近，統(tǒng)一RDA 算法的計算量略有增大。

綜上所述，原RDA 算法以邊權(quán)網(wǎng)絡(luò)系統(tǒng)為主要研究對象，進(jìn)行了詳細(xì)的遞推分解過程的推導(dǎo)，給出完整的可用于計算邊權(quán)網(wǎng)絡(luò)系統(tǒng)可靠度的算法。而點權(quán)網(wǎng)絡(luò)系統(tǒng)和一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)的RDA 算法均以邊權(quán)網(wǎng)絡(luò)系統(tǒng)為基礎(chǔ)，結(jié)合自身特點略有改動。而由于三種賦權(quán)形式網(wǎng)絡(luò)系統(tǒng)RDA算法的區(qū)別，與之對應(yīng)的計算程序也分成三個。為了簡化計算過程，以一般賦權(quán)網(wǎng)絡(luò)系統(tǒng)的RDA算法程序為載體，將邊權(quán)和點權(quán)的計算程序歸并入一般賦權(quán)程序中，又將三部分整合為一體，發(fā)展出統(tǒng)一RDA 算法，不僅使用更加方便，也使RDA 算法程序更加完整緊湊。

2 相依失效工程網(wǎng)絡(luò)系統(tǒng)抗震可靠度估計的RDA 算法

在同一震源下，由于生命線工程結(jié)構(gòu)各設(shè)備間往往存在相互連接的關(guān)系，因此實際上工程網(wǎng)絡(luò)系統(tǒng)是失效相關(guān)的。本章研究相依失效工程網(wǎng)絡(luò)系統(tǒng)抗震可靠度的估計問題，在統(tǒng)一RDA 算法的基礎(chǔ)上，實際上就是研究不交最小路和不交最小割的聯(lián)合概率計算問題，其求解過程包含系統(tǒng)元件(結(jié)構(gòu))極值響應(yīng)邊緣分布模型估計[20?21]和極值響應(yīng)聯(lián)合分布模型估計兩部分。

為了有效解決不交最小路(割)的聯(lián)合概率計算問題，本文先引入移位廣義對數(shù)正態(tài)分布(Shifted Generalized Lognormal Distribution，簡稱SGLD)模型建立結(jié)構(gòu)荷載效應(yīng)的邊緣分布函數(shù)，再分別采用SCM 和Gumbel Copula，建立結(jié)構(gòu)荷載效應(yīng)的聯(lián)合分布函數(shù)，計算聯(lián)合失效概率。最后，將算得的不交最小路和不交最小割的聯(lián)合概率嵌入統(tǒng)一RDA 算法，從而有效估計相依失效生命線工程網(wǎng)絡(luò)系統(tǒng)的可靠度。

2.1 基于SGLD 模型的邊緣分布函數(shù)

SGLD 是一種四階矩概率模型，不僅具有確定的概率密度函數(shù)表達(dá)式和累積分布函數(shù)表達(dá)式，而且具備較廣的偏態(tài)系數(shù)-峰度系數(shù)空間[22]，因此，可用來模擬大多數(shù)的函數(shù)類型?；赟GLD模型，可給出地震作用下結(jié)構(gòu)響應(yīng)邊緣分布的尾部估計。

SGLD 的累積分布函數(shù)的理論表達(dá)式為：

2.2 基于SCM 的不交最小路(割)的聯(lián)合概率計算

SCM 是多維正態(tài)分布的一種新的計算方法，其基本思想是把N維正態(tài)分布函數(shù)看成由N個具有邏輯關(guān)系的事件組成的系統(tǒng)的可靠概率，通過依次把兩個具有邏輯關(guān)系的事件組合成一個復(fù)合事件，最終將整個系統(tǒng)復(fù)合成一個總事件，通過總復(fù)合事件的可靠度指標(biāo)來計算多維正態(tài)分布的聯(lián)合概率函數(shù)值。

SCM 的計算步驟如下：

2.3 基于Gumbel Copula 的不交最小路(割)的聯(lián)合概率計算

本文采用單一參數(shù)多變量Gumbel Copula，主要需要估計其參數(shù)θ。基于兩水準(zhǔn)方法估計參數(shù)時的反應(yīng)樣本，由矩方法的基本思想，讓基于樣本的Kendall’s tau 和基于總體的Kendall’s tau 相等，先得到等級相關(guān)系數(shù)Kendall’s tau，進(jìn)而再計算得到參數(shù)θ。

對于n維Gumbel CopulaCθ(u)，u=(u1,u2,···,un)，其樣本一致性系數(shù)Kendall’s tau 的計算公式為：

2.4 相依失效生命線工程系統(tǒng)抗震可靠度估計的統(tǒng)一RDA 算法

對于相依失效生命線工程系統(tǒng)的可靠度估計問題，先由統(tǒng)一RDA 算法找到所有的不交最小路和不交最小割，再對于每一條不交最小路和不交最小割，考慮相依失效，采用SCM 或Gumbel Copula計算聯(lián)合失效概率，最后回到統(tǒng)一RDA 算法，估計整個系統(tǒng)的抗震可靠度。同時，結(jié)合本文第二章，此方法也可退化應(yīng)用于相依失效的邊權(quán)網(wǎng)絡(luò)系統(tǒng)和點權(quán)網(wǎng)絡(luò)系統(tǒng)。

因此，本文的方法可用于三種賦權(quán)形式的相依失效工程網(wǎng)絡(luò)系統(tǒng)的抗震可靠度估計，使用范圍廣，精確度高，具有廣闊的應(yīng)用前景。

3 算例分析

在各類生命線工程系統(tǒng)中，區(qū)域電力系統(tǒng)結(jié)構(gòu)的地震破壞分析具有代表性。相對于電力系統(tǒng)中的工程結(jié)構(gòu)和高壓電器設(shè)備而言，輸電線路的地震破壞較輕，因此，當(dāng)不考慮電力系統(tǒng)中輸電線路的地震失效時，電力系統(tǒng)被認(rèn)為是一個點權(quán)系統(tǒng)。如果把電力系統(tǒng)中某一結(jié)構(gòu)或設(shè)備作為源點，另一結(jié)構(gòu)或設(shè)備作為終點，那么由源點到終點的整個系統(tǒng)就構(gòu)成了一個拓?fù)渚W(wǎng)絡(luò)，其中各個結(jié)構(gòu)或設(shè)備作為網(wǎng)絡(luò)的節(jié)點，輸電線路作為網(wǎng)絡(luò)的邊。

本章以一假想電力系統(tǒng)為例，將其簡化為下圖2 所示的36 節(jié)點格柵型網(wǎng)絡(luò)，采用第二章提出的算法，估計該網(wǎng)絡(luò)系統(tǒng)相依失效的動力抗震可靠度。先建立系統(tǒng)的動力學(xué)模型，利用SGLD 模型建立結(jié)構(gòu)的邊緣分布函數(shù)，估計元件單體的可靠度。再由元件失效相依，分別采用SCM 和Gumbel Copula，建立結(jié)構(gòu)的聯(lián)合分布函數(shù)，計算聯(lián)合失效概率。最后，采用統(tǒng)一RDA 算法，得到系統(tǒng)的抗震動力可靠度估計結(jié)果。

圖2 格柵型網(wǎng)絡(luò)結(jié)構(gòu)示意圖Fig.2 Schematic diagram of grid network structure

本文采用的相依失效生命線工程系統(tǒng)抗震可靠度估計的統(tǒng)一RDA 算法，不僅可以系統(tǒng)性地考慮網(wǎng)絡(luò)系統(tǒng)結(jié)構(gòu)相依失效的情況，還能在保證精度的前提下大大減少計算量，充分體現(xiàn)該方法的有效性。

3.1 地震作用下的系統(tǒng)動力學(xué)模型

要進(jìn)行36 節(jié)點網(wǎng)絡(luò)系統(tǒng)相依失效下的抗震可靠度評估研究，必須先建立系統(tǒng)的動力學(xué)模型，計算各元件單體的隨機(jī)地震響應(yīng)。由于動力相互作用的存在，有限元分析中的模型變得十分復(fù)雜，計算費時費力，不適用于大量的研究分析工作中。

假設(shè)圖2 所示包含36 個設(shè)備的電力系統(tǒng)，同一列的兩個節(jié)點(如節(jié)點1 和2，節(jié)點3 和4)由硬母線或帶滑移接頭的軟母線連接，因此，需要考慮動力相互作用，而橫排連接母線為一般柔性母線，不引起動力相互作用。

為了簡化電力系統(tǒng)中設(shè)備的非線性隨機(jī)地震動分析計算復(fù)雜的問題，美國學(xué)者Der Kiureghian[24?25]在1999 年最早提出了基于單自由度設(shè)備模型的理想化系統(tǒng)簡化方法，建立了系統(tǒng)的動力學(xué)模型。該方法的基本思路是：用一個假設(shè)的位移形函數(shù)描述設(shè)備的變形，那么每個設(shè)備將會被簡化成一個單自由度振子，這個振子的有效質(zhì)量、有效剛度、有效阻尼以及外部慣性力均與原設(shè)備對應(yīng)。由硬母線或帶滑移接頭軟母線連接的兩個設(shè)備則被理想化為一個雙自由度系統(tǒng)，各參數(shù)仍與設(shè)備自身特性以及連接母線的特性相對應(yīng)，理想化模型示意圖見圖3。

圖3 相連設(shè)備的理想化動力學(xué)模型[26]Fig.3 Idealized dynamic model of connected equipment[26]

隨機(jī)地震作用下，兩設(shè)備間的連接母線進(jìn)行非線性隨機(jī)振動?？紤]動力相互作用，在雙自由度系統(tǒng)內(nèi)建立母線的滯回回復(fù)力模型。

以i、j表示兩個存在相互作用的設(shè)備，具體公式如下：

表2 電器設(shè)備的有效質(zhì)量和剛度Table 2 Values of effective mass and stiffnessof equipment items

本節(jié)建立系統(tǒng)的動力學(xué)模型，用于計算電力系統(tǒng)內(nèi)各單元的隨機(jī)地震響應(yīng)，由響應(yīng)樣本估計單體的邊緣分布和系統(tǒng)的聯(lián)合分布，最終基于RDA算法可得到該電力系統(tǒng)相依失效的動力抗震可靠度。

3.2 結(jié)構(gòu)響應(yīng)隨機(jī)樣本生成

本文利用Mathematica10.0 程序進(jìn)行編程計算，為了保證計算精度，隨機(jī)抽樣生成2000 個地震激勵樣本，作用于該36 節(jié)點網(wǎng)絡(luò)結(jié)構(gòu)上，將結(jié)構(gòu)同一格的兩個節(jié)點理想化為一個雙自由度系統(tǒng)，得到36 個節(jié)點在地震動作用下的2000 組隨機(jī)響應(yīng)數(shù)據(jù)樣本。圖4 為3 組地面地震動加速度時程樣本曲線，圖5 為這3 組地震輸入下，節(jié)點1 對應(yīng)的響應(yīng)曲線。利用這2000 組響應(yīng)樣本，先通過SGLD 模型估計結(jié)構(gòu)隨機(jī)響應(yīng)的邊緣分布，再利用SCM 和Gumbel Copula 函數(shù)分別計算聯(lián)合失效概率，即可得到該結(jié)構(gòu)相依失效下的結(jié)構(gòu)系統(tǒng)抗震可靠度。

圖4 地震譜樣本曲線Fig.4 Sample curves of seismic spectra

圖5 節(jié)點1 對應(yīng)的隨機(jī)響應(yīng)曲線Fig.5 Random response curves corresponding to Node 1

3.3 結(jié)構(gòu)體系邊緣分布分析

基于2000 組隨機(jī)響應(yīng)樣本，本節(jié)分別采用矩方法和兩水準(zhǔn)方法，先識別出SGLD 模型的模型參數(shù)，結(jié)果列于表3 中，再由模型參數(shù)，得到結(jié)構(gòu)體系響應(yīng)尾部分布的近似值。兩水準(zhǔn)方法對應(yīng)的兩個控制點分別取為P1=0.1、P2=0.05。圖6 為其中兩個節(jié)點的超越概率曲線圖，橫軸表示變量X的限值，縱軸表示當(dāng)X取不同限值時，滿足該分布函數(shù)的結(jié)構(gòu)響應(yīng)的超越概率Pf。

圖6 表明：矩方法和兩水準(zhǔn)方法計算得到的超越概率吻合度較高。在分布的中后段，兩種方法給出的超越概率結(jié)果幾乎一致。因此，可以認(rèn)為基于兩水準(zhǔn)方法的SGLD 模型能夠用于本結(jié)構(gòu)響應(yīng)分布的尾部估計，且效率較高。

3.4 結(jié)構(gòu)體系的聯(lián)合概率計算和系統(tǒng)可靠度計算

本節(jié)基于2000 組隨機(jī)響應(yīng)樣本，分別采用SCM 和Gumbel Copula 計算不交最小路和不交最小割的聯(lián)合概率，再將聯(lián)合概率計算結(jié)果嵌入統(tǒng)一RDA 算法，即可算得系統(tǒng)的抗震可靠度。

SCM 需要計算2000 組響應(yīng)樣本的相關(guān)系數(shù)和各元件的可靠度指標(biāo)，再依次把兩個事件組合成一個復(fù)合事件，最終將整個系統(tǒng)復(fù)合成一個總事件，通過總復(fù)合事件的可靠度指標(biāo)來計算多維正態(tài)分布的聯(lián)合概率函數(shù)值。而Gumbel Copula 需要計算Kendall’s tau 和模型參數(shù)θ，結(jié)合邊緣分布結(jié)果來建立Copula 函數(shù)。算得聯(lián)合概率后，將結(jié)果嵌套進(jìn)統(tǒng)一RDA 算法，即可得到結(jié)構(gòu)系統(tǒng)的相依失效可靠度。

由2000 組樣本數(shù)據(jù)算得的相關(guān)系數(shù)為36 階對稱矩陣，因數(shù)據(jù)較多，此處僅列出前6 個節(jié)點之間的相關(guān)系數(shù)矩陣，見表4。算得的Gumbel Copula 所需的τ=0.005185， θ=1.0694。結(jié)構(gòu)系統(tǒng)的可靠度計算結(jié)果見表5。

表3 SGLD 模型參數(shù)Table 3 SGLD model parameters

圖6 各構(gòu)件結(jié)構(gòu)響應(yīng)的超越概率曲線圖Fig.6 Exceedance probability curves of structural response of each component

表4 局部相關(guān)系數(shù)矩陣Table 4 Local correlation coefficient matrix

表5 36 節(jié)點結(jié)構(gòu)網(wǎng)絡(luò)系統(tǒng)可靠度計算結(jié)果Table 5 Reliability calculation results of 36-node network system

圖7 為以Gumbel Copula 可靠度計算結(jié)果為基準(zhǔn)的相對誤差變化曲線，橫坐標(biāo)為RDA 算法中設(shè)定的誤差限值，縱坐標(biāo)為兩種方法可靠度計算結(jié)果的相對誤差。

圖7 兩種方法可靠度計算結(jié)果的相對誤差Fig.7 Relative error of reliability calculation results of two methods

以上結(jié)果表明：一方面，隨著設(shè)定誤差界限的降低，RDA 算法需要計算的不交最小路、不交最小割以及計算時間都在變多，但仍遠(yuǎn)小于傳統(tǒng)遞推分解法，體現(xiàn)了本文計算方法的優(yōu)勢；另一方面，在相同的誤差限值下，SCM 和Gumbel Copula 兩種方法的計算事件數(shù)完全一致，這是因為不交最小路和不交最小割數(shù)均是由RDA 算法確定。兩種方法的結(jié)構(gòu)可靠概率計算結(jié)果非常接近，以Gumbel Copula 為基準(zhǔn)的最大相對誤差不超過0.2‰，但在誤差限值較大的區(qū)域出現(xiàn)較大偏差。同等條件下，Gumbel Copula 的計算時間要低于SCM，且計算過程更為簡潔。

3.5 統(tǒng)一RDA 算法的正確性驗證

當(dāng)不考慮電力系統(tǒng)中輸電線路的地震失效時，本章36 節(jié)點電力系統(tǒng)被認(rèn)為是一個點權(quán)系統(tǒng)，表5 為由統(tǒng)一RDA 算法退化計算的點權(quán)網(wǎng)絡(luò)系統(tǒng)的抗震可靠度結(jié)果。本節(jié)將用SCM 和Gumbel Copula 計算的不交最小路和不交最小割的聯(lián)合概率嵌入點權(quán)的RDA 算法，再次驗證統(tǒng)一RDA 算法的正確性，計算結(jié)果見表6。

表6 點權(quán)網(wǎng)絡(luò)系統(tǒng)的可靠度計算結(jié)果Table 6 Reliability calculation results of node weight network system

計算結(jié)果驗證了統(tǒng)一RDA 算法的正確性，但該網(wǎng)絡(luò)結(jié)構(gòu)由統(tǒng)一算法退化到點權(quán)的不交最小路和不交最小割數(shù)遠(yuǎn)大于原點權(quán)算法，計算量增加。

結(jié)合圖1、圖2 兩個網(wǎng)絡(luò)的可靠度計算結(jié)果，當(dāng)網(wǎng)絡(luò)系統(tǒng)中節(jié)點少、邊多，兩者數(shù)量相差較大，網(wǎng)絡(luò)結(jié)構(gòu)比較復(fù)雜時，統(tǒng)一RDA 算法實用性較高；而當(dāng)網(wǎng)絡(luò)系統(tǒng)中節(jié)點數(shù)與邊數(shù)比較接近，網(wǎng)絡(luò)結(jié)構(gòu)比較簡單時，統(tǒng)一RDA 算法的計算量可能會大大增加，從而導(dǎo)致計算費用增加，不適合應(yīng)用于工程實踐。因此，統(tǒng)一RDA 算法研究的意義主要在于理論簡化，并以此為基礎(chǔ)進(jìn)一步進(jìn)行相依失效情況下的生命線工程系統(tǒng)抗震可靠度估計研究。

綜上所述，雖然統(tǒng)一RDA 算法根據(jù)不同的網(wǎng)絡(luò)系統(tǒng)結(jié)構(gòu)，使用有一定的局限性，但其將原RDA算法整合成一個程序，并且基于統(tǒng)一RDA 算法發(fā)展出的相依失效生命線工程系統(tǒng)抗震可靠度估計的統(tǒng)一RDA 算法，不僅可以較好地考慮元件相依失效，而且可以退化到邊權(quán)和點權(quán)。對于統(tǒng)一RDA 算法計算量增加的問題，計劃后續(xù)將做更深入的研究，對統(tǒng)一RDA 算法中的不交最小路和不交最小割進(jìn)行歸并，從而減小不交最小路數(shù)和不交最小割數(shù)，降低計算量。本文的方法為以后的研究奠定了良好的理論基礎(chǔ)，希望通過進(jìn)一步的研究改進(jìn)，能夠?qū)⑵涓玫貞?yīng)用于工程實踐中。

4 結(jié)論

本文先對原RDA 算法三部分程序進(jìn)行整合簡化，發(fā)展出統(tǒng)一RDA 算法，再在統(tǒng)一RDA 算法的基礎(chǔ)上，引入SCM 和Gumbel Copula 函數(shù)兩種方法，分別計算RDA 算法中不交最小路和不交最小割事件的聯(lián)合失效概率，最后發(fā)展出相依失效生命線工程系統(tǒng)抗震可靠度估計的統(tǒng)一RDA 算法，較好地解決了生命線工程系統(tǒng)中元件失效相關(guān)的可靠度估計問題。算例計算結(jié)果表明：該方法在保證計算精度的前提下，可以用于相依失效生命線網(wǎng)絡(luò)系統(tǒng)的抗震可靠度評估。

本文的主要研究結(jié)論有：

(1)相依失效生命線工程系統(tǒng)抗震可靠度估計的統(tǒng)一RDA 算法不僅可以較好地考慮元件相依失效，而且可以退化到邊權(quán)和點權(quán)，具有良好的適用性；

(2) SCM 和Gumbel Copula 可以用于結(jié)構(gòu)體系可靠度估計中的聯(lián)合概率計算。相比之下，Gumbel Copula 計算過程更加簡潔，計算效率更高；

(3)本文方法適用于大多數(shù)的生命線工程系統(tǒng)抗震可靠度的估計，對于統(tǒng)一RDA 算法退化到邊權(quán)、點權(quán)計算量增加的問題，后續(xù)將做進(jìn)一步的研究。