冒如權(quán)

（海裝駐上海地區(qū)軍事代表室，上海 200000）

引言

艦船裝備抗沖擊能力是關(guān)系到艦船作戰(zhàn)效能與安全性的重要因素之一，目前國(guó)內(nèi)主要采用沖擊試驗(yàn)或抗沖擊動(dòng)力學(xué)分析方法對(duì)艦船設(shè)備抗沖擊能力進(jìn)行驗(yàn)證，驗(yàn)證結(jié)果僅能回答設(shè)備抗沖擊能力是否滿足要求，結(jié)構(gòu)強(qiáng)度、剛度以及設(shè)備功能性能是否滿足研制要求規(guī)定的技術(shù)指標(biāo)，往往無(wú)法對(duì)設(shè)備本身在生產(chǎn)制造、安裝工藝以及材料特性等環(huán)節(jié)的隨機(jī)特性進(jìn)行分析驗(yàn)證，而載荷、材料、結(jié)構(gòu)以及制造工藝等環(huán)節(jié)的隨機(jī)特性決定了設(shè)備自身抗沖擊性能帶有一定的概率分布特性，本文擬采用DDAM方法與Monte-Carlo法相結(jié)合，考慮沖擊過(guò)程中各影響因素的隨機(jī)特性分布，計(jì)算基于ANSYS-PDS環(huán)境下通海閥箱消聲器的抗沖擊響應(yīng)，評(píng)估艦船設(shè)備的抗沖擊可靠性能力。

1 研究背景

在抗沖擊可靠性計(jì)算方面，海軍裝備研究院的馮麟涵選取艦用增壓鍋爐為研究對(duì)象，采用數(shù)值計(jì)算與基于改進(jìn)神經(jīng)網(wǎng)絡(luò)的設(shè)備沖擊響應(yīng)預(yù)報(bào)模型相結(jié)合的方法統(tǒng)計(jì)分析設(shè)備抗沖擊能力的分布特性[1]；哈爾濱工程大學(xué)的曹貽鵬教授選用動(dòng)態(tài)設(shè)計(jì)分析法（DDAM）對(duì)船舶柴油機(jī)排氣消聲器進(jìn)行抗沖擊分析,研究消聲器結(jié)構(gòu)形式、剛度與質(zhì)量分布對(duì)抗沖擊特性的規(guī)律影響[2]；計(jì)晨[3]采用DDAM方法計(jì)算了柴油機(jī)的沖擊特性響應(yīng)。

2 分析模型

目前艦船設(shè)備抗沖擊仿真普遍采用美軍標(biāo)的DDAM法，DDAM法是基于模態(tài)疊加法的結(jié)構(gòu)振動(dòng)設(shè)計(jì)分析方法，該方法首先對(duì)系統(tǒng)動(dòng)力學(xué)模型進(jìn)行模態(tài)分析，得到系統(tǒng)模態(tài)陣型和模態(tài)質(zhì)量后，根據(jù)設(shè)計(jì)沖擊譜得出各階模態(tài)的模態(tài)位移和應(yīng)力，以此為輸入計(jì)算結(jié)構(gòu)的位移和應(yīng)力[4]。

本文的研究對(duì)象海水消聲器安裝位置為甲板部位，屬于乙類(lèi)設(shè)備，其沖擊輸入如表1所示。

甲板安裝位置，A0、V0表達(dá)式為：

基于表1的沖擊環(huán)境，在ANSYS環(huán)境下編制APDL加載程序，計(jì)算消聲器的沖擊動(dòng)力學(xué)響應(yīng)，計(jì)算結(jié)果位移及應(yīng)力響應(yīng)如圖1～2所示。

3 沖擊可靠性分析

在抗沖擊動(dòng)力學(xué)分析的基礎(chǔ)上，引入概率統(tǒng)計(jì)學(xué)，綜合考慮不確定性因素對(duì)設(shè)備抗沖擊能力的影響，一般考慮設(shè)備結(jié)構(gòu)尺寸、材料屬性、輸入載荷等變量，在一定的樣本量統(tǒng)計(jì)基礎(chǔ)上，確定各隨機(jī)變量分布特點(diǎn)。基于“應(yīng)力-強(qiáng)度干涉模型”建立極限狀態(tài)方程，隨后進(jìn)行Monte-Carlo抽樣?？箾_擊可靠性分析的技術(shù)框架是將影響設(shè)備抗沖擊能力的各物理量考慮為服從一定隨機(jī)分布的隨機(jī)變量，利用Monte-Carlo法在一定量本的抽樣頻率基礎(chǔ)上，確定設(shè)備沖擊失效概率，進(jìn)而獲得可靠度及靈敏度等參數(shù)。

可靠性分析在ANSYS-PDS環(huán)境下進(jìn)行，理論基礎(chǔ)是在確定性有限元分析的基礎(chǔ)上，基于隨機(jī)有限元理論，對(duì)極限狀態(tài)方程進(jìn)行大量本抽樣循環(huán)計(jì)算，通過(guò)Monte-Carlo法計(jì)算得出結(jié)構(gòu)可靠度與靈敏度等參數(shù)[5]?？煽啃苑治龅年P(guān)鍵在于分析文件的編制，分析文件包括前處理模塊、有限元分析過(guò)程及后處理模塊，最后建立極限狀態(tài)方程，將全部分析過(guò)程的APDL命令流以txt格式提取出來(lái)，作為PDS循環(huán)仿真的對(duì)象。具體計(jì)算步驟如下：

1）設(shè)定可靠性分析輸入?yún)?shù)

對(duì)通海閥箱消聲器進(jìn)行抗沖擊可靠性分析，首先需要確定各隨機(jī)變量分布類(lèi)型及分布特性，根據(jù)數(shù)據(jù)統(tǒng)計(jì)，部分隨機(jī)變量分布情況如表2所示。

編制確定性有限元分析過(guò)程的APDL命令流，并添加結(jié)果提取命令，生成可靠性分析的分析文件，隨后進(jìn)入蒙特卡洛仿真循環(huán)，選擇拉丁超立方抽樣方法（LHS），定義循環(huán)次數(shù)1 000次，運(yùn)行結(jié)束后在后處理階段對(duì)數(shù)據(jù)進(jìn)行處理。

2）可靠性分析后處理

查看所有隨機(jī)輸入變量、隨機(jī)輸出變量的抽樣過(guò)程，包括樣本點(diǎn)曲線、樣本點(diǎn)均值曲線、累積分布函數(shù)、累積分布直方圖、結(jié)構(gòu)可靠度、輸出變量靈敏度、散點(diǎn)圖等模擬結(jié)果。

表1 沖擊輸入?yún)?shù)

圖1 Z方向位移云圖

圖2 Z方向應(yīng)力云圖

①查看變量抽樣歷史曲線

利用GUI操作查看樣本均值抽樣歷史曲線，見(jiàn)圖3。

從各隨機(jī)變量的抽樣歷史曲線可以看出，在進(jìn)行了250次LHS抽樣后，輸入變量基本收斂，這說(shuō)明LHS法抽樣次數(shù)足夠多，經(jīng)過(guò)1 000次模擬基本可以達(dá)到較為準(zhǔn)確的模擬效果。從圖中可以看出每個(gè)隨機(jī)變量的各項(xiàng)隨機(jī)指標(biāo)，以法蘭螺栓孔徑R0為例，其均值（MEAN）為11 mm，標(biāo)準(zhǔn)差（STDEV）為0.5 mm，不對(duì)稱系數(shù)（SKEW）為-0.015，斜率（KURT）為0.072，最小值（MIN）9.1 mm，最大值（MAX）12.7 mm。

表2 隨機(jī)輸入變量與輸出變量

②查看概率密度直方圖

調(diào)用命令“PDSHIST，LHSPLATE，XXX，REL”查看隨機(jī)輸入變量的直方圖，見(jiàn)圖4。

從圖4可以看出R0和S概率分布直方圖接近概率函數(shù)曲線，即高斯分布，并且曲線光滑，說(shuō)明抽樣次數(shù)可以滿足要求。

③查看可靠度

可靠度數(shù)值由消聲器極限狀態(tài)函數(shù)Z決定，當(dāng)Z＜0時(shí)的統(tǒng)計(jì)數(shù)據(jù)見(jiàn)圖5。

從圖5中可以看出，經(jīng)過(guò)1 000次抽樣模擬，極限狀態(tài)函數(shù)Z的均值為150 MPa，最大值為204 MPa，最小值為86.5 MPa。在95 %的置信區(qū)間下，Z＜0的概率值為0，即通海閥箱消聲器的失效概率為0，可靠度R=100 %。

④查看靈敏度

調(diào)用“PDSENS，LHSPLATE，XXX，BOTH，RANK，0.025”查看輸出變量靈敏度，見(jiàn)圖6，查看各輸入變量對(duì)輸出變量的影響。

從極限狀態(tài)函數(shù)Z的靈敏度可以看出，對(duì)消聲器結(jié)構(gòu)可靠性影響最大的變量是屈服極限S、法蘭螺栓孔徑R0和材料密度。屈服極限S的靈敏度為正，消聲器可靠性將隨著S增大而增大；法蘭螺栓孔徑R0和材料密度的靈敏度為負(fù)，說(shuō)明消聲器可靠性隨著R0和DENSITY的增大而減小。

圖3 樣本點(diǎn)樣本均值抽樣曲線

圖4 概率密度直方圖

圖5 通海閥箱消聲器結(jié)構(gòu)可靠度

從最大應(yīng)力MAXSTR的靈敏度可以看出，對(duì)消聲器最大應(yīng)力影響最大的因素是法蘭螺栓孔徑R0和材料密度。R0和密度的靈敏度為正，說(shuō)明消聲器最大應(yīng)力將隨著R0和材料密度的增大而減小。

⑤查看消聲器可靠性分析的參數(shù)抽樣散點(diǎn)圖

調(diào)用“PDSCAT，LHSPLATE，X，Y，POLY，1”查看X變量與Y變量的抽樣散點(diǎn)圖以分析二者的相關(guān)性，見(jiàn)圖7。

圖6 輸出變量靈敏度

圖7 散點(diǎn)圖

可以看出屈服極限S和極限狀態(tài)函數(shù)Z符合線性正相關(guān)，法蘭螺栓孔徑R0和最大等效應(yīng)力MAXSTR符合正相關(guān)，這與靈敏度的餅狀圖結(jié)果一致。

3）抗沖擊可靠性分析結(jié)果

通過(guò)1 000次Monte-Carlo抽樣仿真計(jì)算結(jié)果顯示，通海閥箱消聲器抗沖擊結(jié)構(gòu)可靠度為100 %；通過(guò)靈敏度分析，消聲器屈服極限S對(duì)極限狀態(tài)函數(shù)Z的靈敏度為正，說(shuō)明材料屈服極限對(duì)消聲器的可靠性影響最大，且隨著屈服極限的增大，可靠度隨之增大；法蘭螺栓孔徑R0和密度的靈敏度為負(fù)，消聲器的可靠性隨著R0和密度的增大而減小。

4 結(jié)論

本文通過(guò)對(duì)典型艦船設(shè)備通海閥箱消聲器進(jìn)行抗沖擊可靠性計(jì)算，結(jié)合抗沖擊仿真與Monte-Carlo循環(huán)抽樣模擬，在ANSYS-PDS環(huán)境中，在確定性有限元數(shù)值計(jì)算結(jié)果的基礎(chǔ)上，考慮抗沖擊過(guò)程中不確定因素對(duì)計(jì)算結(jié)果的影響，對(duì)結(jié)構(gòu)優(yōu)化改進(jìn)，可靠性定量分析具有一定的參考意義。