非線性體系動力可靠性分析的等效Duffing體系法

張振浩，楊偉軍

（長沙理工大學 土木與建筑學院，長沙410076）

非線性體系的動力可靠性問題具有重要的理論與實用意義。但考慮了結構非線性因素后，結構體系在隨機激勵下的隨機反應求解變得十分復雜，其動力可靠性分析的難度也就更大［1］。在非線性體系的隨機振動理論中，等效線性化方法是求解非線性系統隨機響應時應用得最廣泛的方法并且目前仍在得到不斷發展［2－6］。與精確解或數值模擬結果的比較表明，等效線性化方法所給出的二階矩精度通常是令人滿意的，但是，等效線性化方法給出的其它統計量，如相關函數、極值等可能是不可靠的。因此，等效線性化法給出的高安全界限時的超越統計量可能是嚴重錯誤的；又如對非線性阻尼系統，用等效線性化得到的首次超越概率與實際值之差可達數個量級［7］。

鑒于上述原因，學者們致力于尋求更好的近似方法。在這當中，等效非線性體系法實質上就曾經被采用過。等效非線性體系法的思想最早是由Caughey提出的［8］，但他提出的方法只適用于原體系是擬線性的情況。朱位秋等［9－10］建立一種適合于求解擬李亞普諾夫系統隨機響應的等效非線性系統法。該方法采用的“最佳”等效原則是使等效系統與原系統具有相同的平均能量變化規律（即具有相同的漂移和擴散系數）。筆者采用的等效非線性體系法將保留原體系的剛度非線性特征，而將非線性阻尼等效線性化，即將原非線性體系等效為具有線性阻尼、非線性剛度的結構體系，而這類具有線性阻尼、非線性剛度的非線性體系是可以通過FPK方程法求得其穩態反應過程的精確概率分布的。該方法對原非線性體系無特殊要求，具有普適性。

1 擬Duffing體系的等效非線性分析

1.1 平穩激勵下2非線性體系間的等效分析

一般情形下的單自由度非線性體系的振動微分方程可表示為式（1），

設與方程（1）等價的非線性體系為式（2）所示的Duffing體系，

式中，ce、ke分別為等效阻尼系數與等效剛度系數；ε為常數，當ε＝0時結構退化為線性體系；β（X）為X的奇函數，且有

體系之間的誤差可以用方程（1）與（2）之間的差值來表示，以e（t）表示：

誤差項e（t）也是一個隨機過程。為使等效體系最優地逼近原體系，等效準則采用使等效體系與原體系之間的絕對偏差為最小，對于隨機過程，等同于使e（t）的平方的均值（即e（t）的均方值）最?。?1］。按此準則來確定等效參數ce與ke。

根據式（3），有

式（4）可以將 E［e2（t）］看作是等效參數ce、ke的二元函數。根據多元函數求極值的方法，可知，要使E［e2（t）］取最小值就相當于要使式（5）成立。

根據式（5），并利用求數學期望與求導運算間的可交換性，最終整理可得式（6）。

將式（6）聯立組成關于ce和ke的方程組，由此求得等效參數ce、ke，見式（7）。

由式（7）可見，要求解等效參數ce、ke，必須知道上式右端的那些數學期望。在不作任何假設的情況下，這些期望值是很難求得的，因為這需要知道反應X（t）和（t）聯合概率分布，而這是未知的。

因為激勵F（t）為平穩過程，若略去反應過程的過渡階段，直接考慮穩態下的反應狀況，此時根據平穩過程與它的均方導數在同一時刻上總是互不相關的這一結論，可知平穩位移反應X（t）和速度反應（t）互不相關，因此有E［X（t）（t）］＝0。于是，式（7）可簡化為式（8）。

在隨機等效分析中，通常用等效體系反應的聯合概率密度來代替原體系反應的聯合概率密度來確定式（7）、（8）或（9）中的數學期望值。由于式（7）、（8）或（9）中的數學期望值是由等效非線性方程（2）得出的，因此這些期望的表達式中將總含有ce和ke。所以，與等效線性化方法類似，等效非線性化方法中，為最終得到ce和ke的具體值，通用的方法是采用迭代法求解：首先假設ce和ke的初值；然后將其代入等效方程（2），由FPK方程法求出反應X（t）和X·（t）的一階矩、二階矩以及二階聯合矩；再由式（7）或（8）或（9）求出第一次迭代的ce1和ke1；如此重復以上步驟，直到求得滿足收斂準則的ce、ke終值；最后，由ce、ke終值代入等效方程（2），將求得的解作為原非線性體系的近似解。

1.2 非平穩激勵情形的討論

對于非平穩隨機激勵的情況，由式（7）、（8）或（9）可明顯看出，由于ce、ke直接與反應的統計矩有關，而非平穩反應的統計矩是時間t的函數［12］，因此，體系的等效參數是隨時間而變化的，即有

此時，等效阻尼和等效剛度以及反應統計矩，就需要從t1＝Δt的離散時刻起按以上步驟迭代計算，直到計算到所需要的時刻tk＝kΔt。

2 結構動力可靠度

通過等效非線性分析將原非線性體系等效為Duffing體系后，可通過FPK方程求得體系的聯合概率密度，于是可以方便地采用經典Poisson過程法求得體系的動力可靠度。

Poisson過程法計算基于首次超越破壞機制的動力可靠度基本公式可表示為式（11）［13］，

式中，b1、－b2為雙側安全界限，T 為時段長，（t））為反應過程與安全界限的交差速率，可由賴斯公式計算見式（12），［13］

3 算 例

考慮van der Pol振子受高斯白噪聲激勵，

式中，W（t）是譜密度為S0的高斯白噪聲，ε′為常參數。試基于首次超越破壞機制分析該非線性體系的動力可靠性。

體系（13）為一復雜非線性體系，首先構造與其等價的非線性體系見式（14），

式中ce、ke為根據2體系間的誤差最小的準則確定的等效參數。在體系（14）中，它們有明確的物理意義，即ce、ke分別為當ε＝0時Duffing體系退化為線性體系時的阻尼系數和剛度系數。

3.1 等效參數的求解

根據式（9）求出ce、ke的表達式。經比較知：

于是，根據數學期望的運算性質，式（9）中的各項期望可求得為

經過化簡計算，等效參數ce與ke最終可表示為反應均方值E（X2）的函數見式（15）和式（16）。

3.2 等效體系的隨機反應求解

根據FPK方程法求解等效非線性體系（14），結 果如表1所示。

表1 等效體系（14）的解析解

表中，常數C由歸一化條件確定，C＝表示，體系（14）當ε＝0時退化為線性體系時（）（）的位移反應X0t和速度反應t的平穩方差，這2個方差由線性隨機振動分析理論可求得，見式（17）。

式（18）即是關于ce和ke的二元方程組，但難以求其解析解，只能通過數值方法求數值解。確定ce和ke值后，再將其反代入表1中的各式即得各項隨機反應結果。

3.3 分析結果及討論

文獻［9］給出了van der Pol振子的近似聯合概率密度函數，為式（19），

其中erfc（·）為余補誤差函數，erfc（x）＝1－erf（x）＝經與數值結果的對比表明，在ε′為小參數時式（19）給出的結果是精度較高的。

以下就將本文的計算結果與式（19）的結果以及等效線性化法分析的結果進行比較。若取參數ε＝0，因為等效非線性體系退化為了線性體系，所以此時可以得到將原非線性體系等效線性化分析的結果。

由表2計算結果對比可以看到，等效非線性法給出的計算結果與文獻［9］的結果甚接近，結果是可靠的。此外，與等效線性化分析結果的比較表明，等效非線性分析所給出的結果精度確有提高，而且有隨著等效體系非線性參數ε的增大，其計算結果精度也有提高的趨勢。因此，等效非線性分析方法是可行的。由表3的計算結果可見，若將原非線性體系等效為線性體系分析其動力可靠性，誤差確實比較大，結果與采用 Monte－Carlo數值模擬法［15］的計算結果吻合較好，等效非線性法的計算結果精度明顯提高。

表2 隨機反應分析結果對比

表3 動力可靠度計算結果對比

4 結 論

非線性體系的動力可靠性分析精度的高低，其關鍵在于體系的隨機反應分析的精度高低。本文提出的基于等效Duffing體系的等效非線性化法，將具有一般普遍性的非線性體系等效為線性阻尼而剛度非線性的這一類可以通過FPK方程求得其精確穩態聯合概率密度函數的非線性Duffing體系。算例分析表明，本文方法的計算結果精度較之等效線性化法的精度要好，提高了非線性體系動力可靠性分析結果的精度。此外，由于所采用的等效非線性體系中含有可控制體系非線性強弱的參數ε，因此改變ε值的大小便可容易獲得將原非線性體系等效為不同強弱非線性體系時的分析結果；特別的，當取ε為零，便可得到等效線性化的分析結果。這對于問題的研究頗為方便。而ε的最佳取值問題，即ε取多大值時能夠獲得精度最高的非線性體系可靠度計算結果，這將是需要進一步開展的研究工作。

最后應指出，分析過程中為了簡化計算，略去了小參數的高階項，同時處理高階反應矩時所采用的正態降階法也采用了一些近似假設，這些將會對計算結果的精度產生一定影響。

［1］楊偉軍，張振浩.基于連續Markov過程首超時間概率分析的結構動力可靠性研究［J］.工程力學，2011，28（7）：124－129.YANG Weijun，ZHANG Zhenhao.Structural dynamic reliability study based on the probability analysis of first－passage time of continuous markov process［J］.Engineering Mechanics，2011，28（7）：124－129.

［2］蘇亮，王毅.等效線性化方法中系統參數求解的優化算法［J］.工程力學，2011，28（9）：23－29.SU Liang，WANG Yi.New optimization algorithm for determining system parameters of equivalent linearization method［J］.Engineering Mechanics，2011，28（9）：23－29.

［3］Guyader A C，Iwan D.Determining equivalent linear parameters for use in a capacity spectrum method of analysis［J］.Journal of Structural Engineering，2006，132（1）：59－67.

［4］Lin Y Y，Miranda E.Noniterative equivalent linear method for evaluation of existing structures ［J］.Journal of Structural Engineering，2008，134（11）：1685－1695.

［5］陳衍茂，劉濟科.一種改進的等效線性化方法［J］.應用力學學報，2008，25（2）：296－298.CHEN Yanmao，LIU Jike.Improved equivalent linearization method［J］.Chinese Journal of Applied Mechanics，2008，25（2）：296－298.

［6］王偉，劉必燈，周正華，等.剛度和阻尼頻率相關的等效線性化方法［J］.巖土力學，2010，31（12）：3928－3933.WANG Wei，LIU Bideng，ZHOU Zhenghua，et al.Equivalent linear method considering frequency dependent stiffness and damping ［J］.Rock and Soil Mechanics，2010，31（12）：3928－3933.

［7］ZHU W Q，HUANG Z L，SUZUKI Y.Equivalent nonlinear system method for stochastically excited and dissipated partially integrable Hamiltonian systems［J］.International Journal of Non－Linear Mechanics，2001，36（5）：773－786.

［8］Caughey T K.On the response of nonlinear oscillators to stochastic excitation.American Society of Mechanical Engineerings，Applied Meehanics Division，1984，65：9－14.

［9］朱位秋，余金壽.預測非線性系統隨機響應的等效非線性系統法［J］.固體力學學報，1989，10（1）：34－44.ZHU Weiqiu，YU Jinshou.Equivalent nonlinear system method for predicting response of nonlinear systems to random excitations［J］.Acta Mechanica Solida Sinica，1989，10（1）：34－44.

［10］朱位秋.非線性隨機動力學與控制［M］.北京：科學出版社，2003.

［11］張振浩.結構動力可靠性理論研究及RC梁橋抗震可靠性分析［D］.長沙：長沙理工大學，2010.

［12］楊偉軍，張振浩，林立.基于破壞指標界限值的結構抗震可靠度分析［J］.地震工程與工程振動，2010，30（1）：77－83.YANG Weijun，ZHANG Zhenhao，LIN Li.Research on seismic reliability of structure based on boundary value of damage index［J］.Journal of Earthquake Engineering and Engineering Vibration，2010，30（1）：77－83.

［13］Ang A H S，Tang W H.Probability concepts in engineering—emphasis on applications to civil and environmental engineering ［M］.John Wiley ＆ Sons Ltd，2007.

［14］李杰，陳建兵.隨機振動理論與應用新進展［M］.上海：同濟大學出版社，2009.

［15］劉佩，姚謙峰.采用重要抽樣法的結構動力可靠度計算［J］.計算力學學報，2009，26（6）：851－855.LIU Pei，YAO Qianfeng.Dynamic reliability calculation based on importance sampling method［J］.Chinese Journal of Computational Mechanics，2009，29（6）：851－855.