趙世軍,張志成
(山西省檢驗檢測中心,山西 太原 030012)
起重機是一種省力的值班設備,廣泛應用于建筑、橋梁、水利、工業生產等工程項目[1]。橋式起重機因其結構承載能力強、應用范圍廣,而成為最常見的起重機類型。隨著生產力的提高,為了滿足日益繁重的工程量,橋式起重機常常需要日夜進行起重作業。在這種工況下,橋式起重機結構容易發生故障。這種值班設備的故障可能會造成人員傷亡,甚至引發嚴重的社會問題[2]。因此,為了指導橋式起重機的安全設計,有必要確定橋式起重機的結構可靠性相關變量靈敏度。
起重機是一種重型民用機械,具有質量大、體積大等特點[3],這導致了起重機試驗的巨大費用和勞動力成本。此外,由于相關部門對起重機試驗的投資有限,因此起重機試驗及其數據收集是一項挑戰。使用計算模型(如有限元法)可以模擬和預測起重機的性能。然而,由起重機的復雜結構建立的有限元模型(Finite Element Modeling,FEM)的單次計算時間可以達到10 s 甚至更長;同時,對于可靠性問題,FEM 的計算次數可能達到數萬次,因此在解決工程問題時計算時長是不可接受的[4]。近年來,代理模型得到了廣泛的應用,其主要優點是使用替代工具來預測復雜的計算結果。理論上,如果模型足夠精確,那么模型的最優解與實際結果將會非常近似。除了解決可靠性問題外,在結構安全方面,還需要獲得參數的可靠性靈敏度,以指導安全設計。通過靈敏度分析,可得到不同參數對性能的影響程度。在靈敏度研究中,Sobol[5]提出的基于方差的全局靈敏度引起了廣泛關注。Sobol 靈敏度理論通過比較基于失效概率的不確定參數全局靈敏度指標,可衡量不確定參數對可靠性的影響程度。
針對橋式起重機結構可靠性靈敏度計算的難點,本文提出了一種基于人工神經網絡(Artificial Neural Network,ANN)的仿真計算框架。該仿真計算框架基于橋式起重機結構分析,利用Sobol靈敏度理論和ANN的計算方案,研究失效模式下的兩種危險工況的可靠性靈敏度。結果表明,載荷、截面尺寸和彈性模量等參數的可靠性靈敏度指標較高,由此說明這些參數對起重機結構失效的影響較大。
對于輸入變量x={x1,x2,…,xn}和對應響應Y(x),Sobol[5]提出了一種分解公式。這種分解公式可以把響應按照變量的維度分成不同部分,響應的分解公式的表達式為

在式(1)的基礎上,Sobol 同時提出方差的分解公式的表達式為

式中:Vi為變量xi的一階方差貢獻。一階方差貢獻Vi還可以表示為計算公式為

式中:x-i為除xi之外的輸入變量。
由此,對于變量xi,基于方差貢獻的靈敏度指標Si表示為

ANN包括輸入層、隱藏層和輸出層。輸入節點提供來自外部參數的輸入,以形成輸入層;隱藏節點未連接到外部響應并計算輸入節點的信息,將其傳遞給輸出節點;一組輸出節點構成輸出層,將網絡中的信息傳輸到外部世界。本文利用ANN模型,預測大樣本下訓練樣本的極限狀態響應。ANN的預測模型表示為

式中:Fq、P、PG、PH為橋式起重機的相關載荷參數;t1、t2、t3、t4為橋式起重機的結構設計參數;E為彈性模量。
本文提出的仿真計算框架用于綜合評價與評估工程實際中橋式起重機的可靠性靈敏度。給出了QD 系列吊鉤橋式起重機結構可靠性影響參數,見表1。可靠性靈敏度結果涉及失效模式下的兩種危險工況,即滿載小車位于跨中的工況1 和滿載小車位于跨端的工況2。

表1 橋式起重機結構可靠性靈敏度參數
圖1 為不同參數在兩種危險工況下的可靠性靈敏度計算結果。由圖1 可知,在工況1下,截面尺寸t1的靈敏度指標S5最大。這表明,不確定參數t1對強度破壞的影響最大;截面尺寸t3和集中載荷的靈敏度指標S2和S7,僅次于S5,說明這些參數對失效概率也有顯著影響。在工況2下,采用基于ANN的方法和蒙特卡洛仿真(Monte Carlo Simulation,MCS)方法的計算結果相似;與工況1 相比,在工況2 下的結構可靠性靈敏度計算結果發生了變化,應該更多地關注參數截面尺寸t1、集中載荷P、和均布載荷Fq;在這些不確定參數中,結構參數截面尺寸t1對失效概率的影響最大。

圖1 橋式起重機結構可靠性靈敏度計算結果
結果表明,MCS 方法計算結構可靠性靈敏度效率低下,而基于ANN 的方法計算結構可靠性靈敏度具有高效性和精確性。針對結構可靠性靈敏度參數對橋式起重機結構失效的影響,在安全設計時應該更多地關注結構參數t1和t3,集中載荷P 和均布載荷Fq。如何在更大的樣本量下提高仿真計算框架的計算效率是未來研究的方向。