基于響應(yīng)面模型的復(fù)合材料軸機(jī)械連接可靠性分析

盧永文 楊朱偉 楊 堯 魏 濤 秦 瑯

(中國船舶集團(tuán)公司第七一一研究所，上海 201108)

1 研究背景

傳統(tǒng)船舶推進(jìn)系統(tǒng)中軸系重量約占整船重量的2%(約100 t～200 t)，使用復(fù)合材料替代金屬可以使軸系減重25%～80%，因此復(fù)合材料軸在船舶推進(jìn)系統(tǒng)的傳動功率密度提高上具有重要的應(yīng)用前景。同時，復(fù)合材料傳動軸還具備隔聲減振效果，是未來船舶推進(jìn)方式的重要發(fā)展方向。目前在美國海軍YP-654巡邏艇、瑞典海軍SMYGE氣墊船以及美國Sacramento級補(bǔ)給艦等艦艇上均使用了復(fù)合材料軸系[1]。復(fù)合軸系中涉及到大量復(fù)合材料與其他部件的連接，連接形式包括了膠接和機(jī)械連接等形式。其中機(jī)械連接在傳遞載荷能力以及可維修性上具有明顯優(yōu)勢，是復(fù)合材料連接結(jié)構(gòu)的重要形式。但是在復(fù)合材料結(jié)構(gòu)上開孔會對復(fù)合材料強(qiáng)度明顯削弱，同時螺栓以及鉚釘?shù)氖褂眠€會增加結(jié)構(gòu)重量。

復(fù)合材料相對金屬材料具有力學(xué)性能分散性大的特點，傳統(tǒng)的安全系數(shù)法進(jìn)行結(jié)構(gòu)強(qiáng)度設(shè)計時是通過安全系數(shù)對各變量不確定性用一個均值進(jìn)行簡化，會導(dǎo)致設(shè)計結(jié)果危險或保守[2]。在美國航空手冊(MIL-HDBK-17-1F)中根據(jù)復(fù)合材料部件是否屬于關(guān)鍵件，定義了復(fù)合材料結(jié)構(gòu)設(shè)計的A基準(zhǔn)值和B基準(zhǔn)值，表示在給定置信度γ條件下，可以期望β×100%的樣本值大于該下限值的問題[3]。因此，復(fù)合材料基準(zhǔn)值設(shè)計法是基于力學(xué)性能的下分位進(jìn)行分析，表明結(jié)構(gòu)設(shè)計保守。結(jié)構(gòu)的概率可靠性設(shè)計方法則是基于概率論與數(shù)理統(tǒng)計理論，將結(jié)構(gòu)設(shè)計中的材料性能的分散性以及外載荷的隨機(jī)性等不確定因素直接進(jìn)行了量化，通過計算得到結(jié)構(gòu)的失效概率充分考慮了材料的分散性[4]。因此對于復(fù)合材料傳動軸機(jī)械連接部位進(jìn)行強(qiáng)度可靠性分析具有重要的意義。

某型碳纖維復(fù)合材料傳動軸額定工況下承擔(dān)著大轉(zhuǎn)矩，目前采用4倍額定工況進(jìn)行確定性設(shè)計，對于材料的隨機(jī)性未加以考量，因此需要對結(jié)構(gòu)進(jìn)行可靠性分析，對結(jié)構(gòu)設(shè)計進(jìn)行指導(dǎo)。

本文首先對結(jié)構(gòu)強(qiáng)度可靠性分析方法進(jìn)行研究，建立了基于響應(yīng)面模型的結(jié)構(gòu)強(qiáng)度可靠度計算方法。隨后對某型復(fù)合材料傳動軸機(jī)械連接結(jié)構(gòu)進(jìn)行了多種額定扭轉(zhuǎn)工況下強(qiáng)度的確定性分析，最后基于確定性計算結(jié)果對復(fù)合材料軸機(jī)械連接部位強(qiáng)度可靠性進(jìn)行了計算分析。

2 基于響應(yīng)面模型的結(jié)構(gòu)可靠性分析

2.1 結(jié)構(gòu)可靠性理論

結(jié)構(gòu)可靠性是指結(jié)構(gòu)或結(jié)構(gòu)構(gòu)件在設(shè)計工作期內(nèi)滿足各項規(guī)定要求的能力,用可靠度來度量可靠性的大小，表示結(jié)構(gòu)在設(shè)計期限內(nèi)，依據(jù)設(shè)計工況使用時，保持結(jié)構(gòu)的安全性、功能性以及耐久性等功能狀態(tài)的能力[5]。

目前，針對結(jié)構(gòu)可靠性的分析模型主要包括概率可靠性分析模型、模糊可靠性分析模型以及非概率可靠性分析模型三種。其中概率可靠性分析模型以及模糊可靠性分析模型都屬于概率可靠性分析模型，概率可靠性分析模型主要針對變量可以依據(jù)大量數(shù)據(jù)統(tǒng)計得到其概率分布，表明了變量的客觀性，而模糊可靠性模型則主要是依據(jù)主觀判斷來量化一些難以通過大量數(shù)據(jù)得到變量概率分布的方法，體現(xiàn)了變量的模糊性。非概率可靠性分析模型則完全不涉及到概率，描述的是對于一些不確定數(shù)據(jù)較少，采用集合模型來描述這些輸入的不確定因素。

結(jié)構(gòu)可靠性采用可靠度來進(jìn)行量化，代表了結(jié)構(gòu)在規(guī)定時間和條件下完成預(yù)定功能的概率，用Pr表示。由于結(jié)構(gòu)的可靠性是與時間以及工作條件密切相關(guān)的，同時預(yù)定功能也是特定結(jié)果，因此在進(jìn)行結(jié)構(gòu)可靠性計算時必須明確。 結(jié)構(gòu)可靠度受到材料性能、結(jié)構(gòu)的幾何尺寸以及外載荷作用等影響，這些影響因素可以用基本變量的形式來表示，并構(gòu)建一個隱式結(jié)構(gòu)功能函數(shù)。設(shè)隨機(jī)變量為X=(X1,X2,X3,…,Xn)，則結(jié)構(gòu)應(yīng)力狀態(tài)包括失效狀態(tài)、臨界狀態(tài)和安全狀態(tài)三種狀態(tài)，不同應(yīng)力狀態(tài)的結(jié)構(gòu)功能函數(shù)如下所示：

其中，Z(X)表示結(jié)構(gòu)功能狀態(tài)函數(shù)，R(X)表示結(jié)構(gòu)所有抗力效應(yīng)，S(X)表示結(jié)構(gòu)載荷效應(yīng)。

結(jié)構(gòu)可靠與結(jié)構(gòu)失效是兩個互不相容事件，因此可以得到結(jié)構(gòu)的可靠度與結(jié)構(gòu)失效概率的大小滿足：

Pr+Pf=1

假設(shè)隨機(jī)變量的失效概率密度函數(shù)為fX(x1,x2,…,xn)，則可以得到結(jié)構(gòu)的失效概率：

由Z(X)=R(X)-S(X)可以得到：

從而得到結(jié)構(gòu)可靠度的大小為：

Pr=1-Pf

2.2 響應(yīng)面模型

目前針對結(jié)構(gòu)可靠性的計算方法包括一次二階矩法和蒙特卡洛法，其中一次二階矩法包括中心點法和驗算點法。中心點法和驗算點法是解析方法，他們的區(qū)別在于，中心點法不考慮變量的分布形式，在平均值處進(jìn)行泰勒展開進(jìn)行結(jié)構(gòu)可靠度的解析求解。驗算點法則通過對變量進(jìn)行當(dāng)量正態(tài)化，將不滿足正態(tài)分布的變量進(jìn)行變換使得其滿足正態(tài)分布，進(jìn)而進(jìn)行結(jié)構(gòu)可靠度的解析求解。

蒙特卡洛法是基于概率論的大數(shù)理論，通過對大量樣本點進(jìn)行統(tǒng)計計算，失效次數(shù)與總的試驗總次數(shù)的比值就是結(jié)構(gòu)的可靠度。由于蒙特卡洛法的理論基礎(chǔ)是大數(shù)理論，因此在使用蒙特卡洛法進(jìn)行可靠性計算時沒有太多的限制。但是如果需要更為可靠的計算結(jié)果則需要進(jìn)行成千上萬次的試驗。因此單純使用蒙特卡洛法進(jìn)行結(jié)構(gòu)的可靠性計算存在計算量過大的問題。針對該問題，目前也提出了一些解決方法，主要是通過將結(jié)構(gòu)隱式狀態(tài)方程進(jìn)行近似擬合，從而實現(xiàn)大大降低計算量大小。

為了解決在可靠性計算中隱式狀態(tài)方程難以求解的問題，研究人員發(fā)展了響應(yīng)面模型、神經(jīng)網(wǎng)絡(luò)模型以及kriging模型等近似模型，通過將目標(biāo)函數(shù)、約束條件與設(shè)計變量之間隱式狀態(tài)函數(shù)使用近似模型進(jìn)行表達(dá)，得到了狀態(tài)方程的顯示表達(dá)式[6]。響應(yīng)面模型具有使用簡單、通用性強(qiáng)、獨立性高以及能與多種商業(yè)軟件結(jié)合等多種優(yōu)點，在本文的計算中選用響應(yīng)面模型進(jìn)行計算。

利用響應(yīng)面模型進(jìn)行可靠性計算時，基于試驗設(shè)計得到計算中的樣本點，并選取適當(dāng)?shù)暮瘮?shù)形式對響應(yīng)面進(jìn)行表示，因此響應(yīng)面的函數(shù)形式會影響最終響應(yīng)面擬合的效率和精度。假設(shè)響應(yīng)面包括隨機(jī)變量X=(X1,X2,X3,…,Xn)，可以得到目標(biāo)函數(shù)與隨機(jī)變量之間的響應(yīng)面函數(shù)關(guān)系如下：

y=β1g1(X1,X2,X3,…,Xn)+…+
βmgm(X1,X2,X3,…,Xn)+ε

其中g(shù)m為所選取的函數(shù)形式，βm為回歸系數(shù)，ε為隨機(jī)誤差[7]。為了評估響應(yīng)面對原始模型的擬合精度，需要對響應(yīng)面函數(shù)進(jìn)行方差分析，計算此時的相關(guān)系數(shù)R2，相關(guān)系數(shù)越接近1表明響應(yīng)面擬合效果越好。

目前應(yīng)用較為廣泛的函數(shù)形式是線性多項式和二次多項式，其中線性多項式函數(shù)形式待定系數(shù)較少，所需樣本點少，計算量小，但缺少交叉項的考慮，對于隱式狀態(tài)方程的非線性行為缺乏考慮。二次多項式在一定程度上克服了線性多項式缺乏對非線性行為考慮的問題，但是待定系數(shù)也顯著增加，使得計算量大大提高?？紤]到本文目的在于建立一套針對復(fù)合材料結(jié)構(gòu)的可靠性計算方法，因此采用了一次多項式進(jìn)行響應(yīng)面擬合，其基本函數(shù)形式如下[8]。

2.3 樣本設(shè)計方法

進(jìn)行響應(yīng)面擬合時，樣本點選取會對擬合速率與精度產(chǎn)生影響，因此采用恰當(dāng)?shù)臉颖驹O(shè)計方法至關(guān)重要。

拉丁超立方設(shè)計方法是一種能夠用樣本點很好的充滿整個空間的試驗設(shè)計方法，其基本原理是假設(shè)在M維空間中，樣本數(shù)量為N，然后將每個設(shè)計變量劃分為N份，這樣就可以將設(shè)計空間分為NM個子空間[9]。對于單個設(shè)計變量而言，通過N次抽樣得到的樣本空間就可以代表這個變量在這個設(shè)計空間中的隨機(jī)分布。其具有在任一維空間中都有與樣本數(shù)量相同的子區(qū)間，并且在子區(qū)間中只有有且僅有一個樣本，樣本在每個子空間可以隨機(jī)選取。但是其不能保證樣本點在空間中分布的均勻性，基于拉丁超立方設(shè)計形成了最優(yōu)拉丁超立方設(shè)計方法，該方法兼顧了樣本點的正交性和空間中分布的均勻性。如圖1所示是九個樣本點在空間中的分布對比，從圖中樣本點的分布可以看出最優(yōu)拉丁超立方設(shè)計方法能夠更好的對空間進(jìn)行描述[10]。因此，本文進(jìn)行響應(yīng)面擬合時選用的樣本設(shè)計方法為最優(yōu)拉丁超立方設(shè)計方法。

圖1 (1)拉丁超立方設(shè)計方法(2)最優(yōu)拉丁超立方設(shè)計方法

3 算例

3.1 復(fù)合材料軸計算模型

3.1.1 計算模型簡化

本文計算模型依據(jù)某型號復(fù)合材料傳動軸進(jìn)行建立，針對螺栓連接結(jié)構(gòu)進(jìn)行了簡化，建立的1/16模型如圖2所示。其中復(fù)合材料軸的厚度為66 mm，鋪層順序為[±45°]，金屬法蘭與螺栓均采用45#鋼，彈性模量210 GPa，泊松比0.3。計算選用Abaqus/Standard模塊，采用C3D8R單元，單元數(shù)為7 804，節(jié)點數(shù)為20 647。對金屬法蘭施加固支約束，螺栓、金屬法蘭以及復(fù)合材料之間的接觸采用罰函數(shù)來模擬摩擦接觸。

圖2 復(fù)合材料軸連接結(jié)構(gòu)有限元模型Fig.2 Finite element model of composite shaft connection structure

3.1.2 強(qiáng)度準(zhǔn)則

碳纖維增強(qiáng)復(fù)合材料層合結(jié)構(gòu)的破壞形式分為首層失效準(zhǔn)則和最終層失效準(zhǔn)則，首層失效準(zhǔn)則是當(dāng)出現(xiàn)某一層失效時就判定結(jié)構(gòu)發(fā)生破壞，最終層失效準(zhǔn)則是當(dāng)所有層都發(fā)生破壞時才判定結(jié)構(gòu)發(fā)生破壞。在實際的工程應(yīng)用中，出于安全的考慮，通常采用保守的設(shè)計方案，因此在本文的計算中采用首層失效準(zhǔn)則進(jìn)行損傷判斷。

復(fù)合材料的失效準(zhǔn)則包括最大應(yīng)力準(zhǔn)則、最大應(yīng)變準(zhǔn)則、Hashin失效準(zhǔn)則以及蔡吳(Tasi-Wu)失效準(zhǔn)則等，目前在工程應(yīng)用領(lǐng)域Tasi-Wu失效準(zhǔn)則應(yīng)用較為廣泛，因此本文的計算中選用Tasi-Wu失效準(zhǔn)則作為材料的失效判據(jù)[11]。其具體形式如下。

式中：

在本文的計算中取I12、I13、I23為0，由此可以得到在平面應(yīng)力狀態(tài)下的Tasi-Wu失效準(zhǔn)則：

當(dāng)Tasi-Wu失效系數(shù)大于1時判定材料發(fā)生破壞。

3.1.3 材料

本文計算的復(fù)合材料軸使用T800纖維增強(qiáng)復(fù)合材料，單層板的材料參數(shù)如表1所示，單層厚度0.3 mm。

表1 復(fù)合材料單層板材料參數(shù)

3.2 復(fù)合材料軸機(jī)械連接部位強(qiáng)度計算

本文計算所用的復(fù)合材料傳動軸額定轉(zhuǎn)矩為4 000 KNm，在4倍額定轉(zhuǎn)矩工況時進(jìn)行強(qiáng)度計算，計算結(jié)果如圖3所示。

圖3 4倍額定轉(zhuǎn)矩工況強(qiáng)度計算結(jié)果

由計算結(jié)果可知，在4倍額定轉(zhuǎn)矩工況下復(fù)合材料傳動軸的最大蔡吳失效系數(shù)為0.888<1，軸的強(qiáng)度滿足設(shè)計要求。隨后計算了復(fù)合材料傳動軸在5倍額定工況下的強(qiáng)度，計算結(jié)果如圖4所示。

圖4 5倍額定轉(zhuǎn)矩工況強(qiáng)度計算結(jié)果Fig.4 Five times the rated torque strength calculation results

由計算結(jié)果可知，在5倍額定轉(zhuǎn)矩工況下復(fù)合材料軸的最大蔡吳失效系數(shù)為1.110>1，強(qiáng)度不滿足設(shè)計要求。

綜合上述強(qiáng)度計算結(jié)果可知，復(fù)合材料傳動軸在4倍額定轉(zhuǎn)矩工況下強(qiáng)度滿足設(shè)計要求，而在5倍額定轉(zhuǎn)矩工況時則不安全。因此需要考慮材料分散性對結(jié)構(gòu)強(qiáng)度的影響，對4倍額定轉(zhuǎn)矩時的結(jié)構(gòu)強(qiáng)度進(jìn)行可靠性分析。

3.3 結(jié)構(gòu)強(qiáng)度可靠性計算

通過3.2中的確定性分析分析結(jié)果，得到了有限元計算的模型輸入文件(.inp)和計算結(jié)果文件(.dat)?；陧憫?yīng)面的結(jié)構(gòu)可靠性計算方法，在Isight中對狀態(tài)方程的響應(yīng)面進(jìn)行擬合并進(jìn)行蒙特卡洛模擬，搭建的計算流程如圖4所示。

圖5 在Isight中搭建的計算流程

根據(jù)擬合的響應(yīng)面的精度的不同，需要的樣本點數(shù)量也存在差異。對4倍額定轉(zhuǎn)矩工況的結(jié)構(gòu)可靠性進(jìn)行分析，選取20個樣本點進(jìn)行線性響應(yīng)面擬合，得到此時響應(yīng)面擬合精度R2=1.00，證明擬合的響應(yīng)面能夠很好的對原始狀態(tài)方程進(jìn)行描述。

隨后對擬合得到的響應(yīng)面進(jìn)行10萬次蒙特卡洛模擬，得到結(jié)構(gòu)可靠度為0.94。根據(jù)航空中對于非關(guān)鍵件可靠性要求的B基準(zhǔn)值0.94，此時結(jié)構(gòu)可靠度0.94>0.90，說明目前結(jié)構(gòu)設(shè)計給予的安全系數(shù)4滿足可靠性要求。由于此時結(jié)構(gòu)強(qiáng)度可靠度高于B基準(zhǔn)值，因此可以降低此時結(jié)構(gòu)可靠度，通過減少目前結(jié)構(gòu)鋪層厚度來提高材料利用率。

如圖6所示，通過分析各輸入變量對輸出結(jié)果的影響，得到影響結(jié)構(gòu)的蔡吳失效系數(shù)最大值的主要輸入變量為X1t、X2t和。其中X1t、X2t與蔡吳失效系數(shù)的大小呈現(xiàn)負(fù)相關(guān)，G23與蔡吳失效系數(shù)的大小則呈現(xiàn)正相關(guān)。因此可以在保持結(jié)構(gòu)厚度不變的情況下改用X1t、X2t更小或G23更大的復(fù)合材料進(jìn)行傳動軸的結(jié)構(gòu)設(shè)計。

圖6 4倍額定轉(zhuǎn)矩時各輸入變量對輸出結(jié)果的影響

4 結(jié)論

通過對結(jié)構(gòu)可靠性計算方法進(jìn)行研究，建立了基于響應(yīng)面模型的結(jié)構(gòu)可靠性分析方法，并利用Isight集成Abaqus，計算了復(fù)合材料軸機(jī)械連接部位強(qiáng)度可靠度大小，獲得的主要結(jié)論如下：

(1)建立了一套基于響應(yīng)面模型的復(fù)合材料結(jié)構(gòu)可靠性分析方法，通過對隱式狀態(tài)方程近似擬合，解決了復(fù)雜結(jié)構(gòu)狀態(tài)方程難以求解的問題，具有重要的工程意義；

(2)根據(jù)復(fù)合材料軸機(jī)械連接部位的結(jié)構(gòu)強(qiáng)度可靠度的計算結(jié)果，在4倍額定轉(zhuǎn)矩工況下的結(jié)構(gòu)可靠度滿足B基準(zhǔn)值要求。因此，在可靠度滿足B基準(zhǔn)值的情況下，可以通過減少復(fù)合材料軸的鋪層厚度來提高材料利用率；

(3)根據(jù)各輸入變量對蔡吳失效系數(shù)的貢獻(xiàn)比例，可以在保持結(jié)構(gòu)鋪層厚度不變的情況下改用、較小或較大的復(fù)合材料來進(jìn)行軸的設(shè)計。