一種路基穩(wěn)定可靠度的計算方法
【專利摘要】本發(fā)明公開了一種路基穩(wěn)定可靠度的計算方法，它利用量子粒子群算法優(yōu)化LSSVM參數(shù)，利用LSSVM建立輸入?yún)?shù)與計算結(jié)果間的映射關(guān)系，再進行路基可靠度的計算，以解決蒙特卡羅法計算可靠度耗時過多的限制的問題，為路基穩(wěn)定可靠度計算提供了一條實用途徑。
【專利說明】一種路基穩(wěn)定可靠度的計算方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及最優(yōu)化及可靠度計算領(lǐng)域，具體地指一種路基穩(wěn)定可靠度的計算方 法。

【背景技術(shù)】
[0002] 目前鐵路路基穩(wěn)定的計算采用的是總安全系數(shù)法，將影響安全的各種因素概括為 一個安全系數(shù)，但不能反映工程的實際安全水平，限制了鐵路工程設(shè)計水平的發(fā)展。因此， 將設(shè)計參數(shù)視為隨機變量，用可靠指標(biāo)描述路基穩(wěn)定安全度更加合理。蒙特卡羅法是一種 較精確的可靠度計算方法，其他可靠度計算方法計算結(jié)果常常用蒙特卡洛法進行驗證。但 采用蒙特卡洛法計算可靠度時，需要進行大批量抽樣，抽樣次數(shù)NS100/Pf。如果失效概率 較小，則計算次數(shù)必須足夠大才能保證可靠度的計算精度。復(fù)雜結(jié)構(gòu)一次抽樣計算少則幾 分鐘，多則幾個小時甚至幾天，要進行成千上萬次的抽樣計算在實際工程中是不現(xiàn)實的。


【發(fā)明內(nèi)容】

[0003] 本發(fā)明的目的就是要解決上述【背景技術(shù)】的不足，提供一種路基穩(wěn)定可靠度的計算 方法，以解決蒙特卡羅法計算可靠度耗時過多的問題。
[0004] 本發(fā)明的技術(shù)方案為：一種路基穩(wěn)定可靠度的計算方法，其特征在于，它包括以下 步驟：
[0005]1)、分析實際問題，確定隨機變量（設(shè)隨機變量的個數(shù)為N)，確定計算模式與計算 方法，建立極限狀態(tài)方程Z;
[0006]2)、采用正交試驗設(shè)計方法提供有代表性的訓(xùn)練樣本：
[0007] 2. 1)、在隨機變量的 3 〇 范圍內(nèi)丨/?, -V- 1 ?Sr,講，R+ 1 _<5"/7/, +3cr,'}選擇并 進行N因素五水平正交試驗設(shè)計選擇訓(xùn)練樣本；
[0008]2. 2)、對訓(xùn)練樣本進行尺度變換，使得各隨機變量在訓(xùn)練中具有相同的地位；

【權(quán)利要求】
1. 一種路基穩(wěn)定可靠度的計算方法，其特征在于，它包括以下步驟： 1) 、分析實際問題，確定隨機變量（設(shè)隨機變量的個數(shù)為N)，確定計算模式與計算方 法，建立極限狀態(tài)方程Z; 2) 、采用正交試驗設(shè)計方法提供有代表性的訓(xùn)練樣本： 2. 1)、在隨機變量的3σ范圍內(nèi)：m, - 3cr;，m, - 1 ·5τ,講，柯,+ 1 ·?,w; +3σ,丨選擇并進行 N因素五水平正交試驗設(shè)計選擇訓(xùn)練樣本； 2. 2)、對訓(xùn)練樣本進行尺度變換，使得各隨機變量在訓(xùn)練中具有相同的地位； -- (式 1) Ui 式1中：Xi為隨機變量X的第i個樣本，X'i為尺度變換以后的樣本，μi和σi分別 為隨機變量的均值及其方差； 2. 3)、利用步驟1)中建立的極限狀態(tài)方程計算出訓(xùn)練樣本本Xi對應(yīng)的計算結(jié)果yi(Z >〇輸出I;Z彡O輸出-1); 3) 、采用LSSVM建立樣本輸入值與輸出值之間的映射關(guān)系，采用量子粒子群算法對 LSSVM參數(shù)進行優(yōu)化： 3.1) 、確定LSSVM形狀參數(shù)δ和懲罰系數(shù)γ的取值范圍，然后將取值范圍變換到 [〇, 1]區(qū)間； 3. 2)、將整個尺度變換后的訓(xùn)練樣本fi隨機的分成η組樣本集，取出其中的一組樣 本集為驗證集，用剩下的樣本集作為訓(xùn)練集來設(shè)計預(yù)測模型，設(shè)兌為驗證集通過預(yù)測模型 得到的估計值，計算出相應(yīng)驗證集的預(yù)測誤差E1 ; 3. 3)、再把步驟3. 2)中取出的樣本集放回原訓(xùn)練樣本中，取出另外一組樣本集作為驗 證集，用剩下的樣本集作為訓(xùn)練集來設(shè)計預(yù)測模型，同理計算相應(yīng)的預(yù)測誤差為E2 ; 3. 4)、重復(fù)設(shè)計模型η次，計算出相應(yīng)驗證集的預(yù)測誤差為EpΕ2、Ε3···En ; 3. 5)、將各個相應(yīng)驗證集的預(yù)測誤差求和： MSE = ±Ej(式2) 7=1 式2中，Ej為對應(yīng)驗證集的預(yù)測誤差值； 3. 6)、參數(shù)選擇問題轉(zhuǎn)化為優(yōu)化問題： MinimizeMSE Subjecttoδr ^ [0, 1]Y1 ^ [0, 1] (式 3); 3. 7)、采用量子粒子群算法來求解此優(yōu)化問題，確定LSSVM的形狀參數(shù)δ和懲罰系數(shù) Y; 4) 、計算可靠度： 4.1) 、按照隨機變量各自的概率分布生成隨機抽樣點：首先在區(qū)間（0，1)內(nèi)產(chǎn)生均勻 分布的隨機數(shù)，然后再變換成給定分布的隨機數(shù)； 4. 2)、將隨機抽樣點送入訓(xùn)練好的LSSVM，得到對應(yīng)的輸出量λ; 4. 3)、設(shè)隨機抽樣點總個數(shù)為Ν，統(tǒng)計其中失效的個數(shù)為Nf，則結(jié)構(gòu)的失效概率為Pf = Nf/N； 4. 4)、根據(jù)失效概率與可靠度間的關(guān)系，即可計算得結(jié)構(gòu)的可靠度指標(biāo)。
2. 根據(jù)權(quán)利要求1所述的一種路基穩(wěn)定可靠度的計算方法，其特征在于：步驟3)中， 驗證集的預(yù)測誤差Ej計算公式為： m1 Ej_λΙ (式 4) 式4中，m為每一組樣本個數(shù)，yi為計算輸出值，兌為模型估計值。
3. 根據(jù)權(quán)利要求1所述的一種路基穩(wěn)定可靠度的計算方法，其特征在于：步驟4. 1) 中，采用乘同余法在區(qū)間（〇，1)內(nèi)產(chǎn)生均勻分布的隨機數(shù)。
【文檔編號】G06F19/00GK104318121SQ201410623520
【公開日】2015年1月28日 申請日期:2014年11月7日 優(yōu)先權(quán)日:2014年11月7日
【發(fā)明者】汪瑩?dān)Q, 趙新益, 曾長賢, 劉慶輝, 陳世剛, 管海濤, 楊松, 蔣道君 申請人:中鐵第四勘察設(shè)計院集團有限公司