郭 慶

(上海電器科學(xué)研究所（集團(tuán)）有限公司，上海 200063）

0 引言

爆炸載荷作用下的機(jī)械沖擊環(huán)境是船舶生命力的基礎(chǔ)。船用配電板在遭受水下爆炸沖擊時(shí)，內(nèi)部設(shè)備可能產(chǎn)生的變形、損壞或功能喪失等情況，對(duì)船舶動(dòng)力系統(tǒng)影響重大[1]。因此，電氣設(shè)備的抗沖擊能力的提升對(duì)提升船舶生命力具有重大意義。目前，主要采用實(shí)船爆炸、沖擊機(jī)試驗(yàn)以及數(shù)值仿真三種方法對(duì)設(shè)備抗沖擊性能進(jìn)行研究[2]。實(shí)船爆炸方式雖具有高可靠性、高精度，但試驗(yàn)耗費(fèi)大，目前只在美軍中采用。沖擊機(jī)試驗(yàn)方法，雖解決了實(shí)船試驗(yàn)耗費(fèi)大的問題，但對(duì)于參試品，只能判斷試驗(yàn)是否通過，對(duì)未通過產(chǎn)品不能全面徹底分析。與前兩種方法相比較，數(shù)值仿真方法則具有經(jīng)濟(jì)、靈活、不需真實(shí)物理模型和周期短等優(yōu)點(diǎn)，目前正逐步推廣。

李華橋[3]對(duì)核電廠儀控設(shè)備斷路器系統(tǒng)抗沖擊進(jìn)行了研究，利用有限元軟件和PDCA循環(huán)解決了沖擊載荷下柜門斷路器跳閘問題。洪德新[4]以船用大容量配電板為研究對(duì)象，仿真計(jì)算了不同沖擊環(huán)境下的柜體響應(yīng)，但未進(jìn)行相關(guān)試驗(yàn)考核，即仿真計(jì)算結(jié)果只具有一定的參考價(jià)值。周馬俊[5]等基于沖擊譜的時(shí)域信號(hào)計(jì)算了船舶風(fēng)機(jī)雙層隔振裝置的沖擊響應(yīng)，并與試驗(yàn)進(jìn)行了對(duì)比，計(jì)算結(jié)果與試驗(yàn)結(jié)果比較符合，對(duì)實(shí)際工程的抗沖擊性能分析具有指導(dǎo)意義。張毅敏[6]等對(duì)船舶設(shè)備分別用時(shí)域分析法和瞬態(tài)模態(tài)動(dòng)力學(xué)方法進(jìn)行了計(jì)算，結(jié)果表明瞬態(tài)模態(tài)動(dòng)力學(xué)方法是一種有效的船舶設(shè)備抗沖擊性能分析評(píng)估方法。賀少華[7]等從理論方面對(duì)船舶設(shè)備沖擊響應(yīng)計(jì)算方法進(jìn)行了全面評(píng)估，總結(jié)了各方法的適用場(chǎng)合，對(duì)設(shè)備抗沖設(shè)計(jì)提供了理論支撐。姚熊亮[8]等以船舶主汽輪機(jī)組為計(jì)算實(shí)例，通過抗沖擊計(jì)算分析，發(fā)現(xiàn)了設(shè)備的薄弱環(huán)節(jié)和危險(xiǎn)區(qū)域，對(duì)同類產(chǎn)品抗沖擊設(shè)計(jì)及優(yōu)化具有一定的指導(dǎo)意義。

本文針對(duì)某型船，從船舶沖擊計(jì)算方法基礎(chǔ)理論出發(fā)，選取了一種適合某船用配電板沖擊仿真的計(jì)算方法。通過有限元仿真計(jì)算了不同隔振系統(tǒng)在配電板沖擊下的隔振效果，從中選取了最優(yōu)方案并進(jìn)行了相關(guān)沖擊試驗(yàn)，樣機(jī)滿足工程實(shí)際相關(guān)需求。本文仿真計(jì)算了三種隔振系統(tǒng)的隔振效果，計(jì)算結(jié)果表明“結(jié)構(gòu)三”設(shè)計(jì)方案在同等沖擊條件下應(yīng)力及應(yīng)變最小，系統(tǒng)減震效果最好。

1 沖擊計(jì)算方法

船舶設(shè)備的沖擊計(jì)算方法大致有靜態(tài)等效法、動(dòng)態(tài)設(shè)計(jì)分析法（DDAM）和時(shí)域模擬法[9]。靜態(tài)等效法只考慮了受沖擊結(jié)構(gòu)的質(zhì)量效應(yīng)，對(duì)設(shè)備材料及結(jié)構(gòu)特性隨沖擊的變化不予考慮且該方法只考慮了一階模態(tài)響應(yīng),當(dāng)設(shè)備破壞的主要原因?yàn)榈碗A模態(tài)時(shí)適用；動(dòng)態(tài)設(shè)計(jì)法，其本質(zhì)是頻域有限元中的模態(tài)疊加法，對(duì)結(jié)構(gòu)的模態(tài)特性進(jìn)行分析，得到相應(yīng)的模態(tài)振型、模態(tài)質(zhì)量，對(duì)應(yīng)邊界載荷可以得到模態(tài)頻率下的應(yīng)力應(yīng)變；時(shí)域模擬法則采用實(shí)測(cè)的時(shí)間歷程曲線作為載荷輸入，對(duì)設(shè)備進(jìn)行瞬態(tài)動(dòng)響應(yīng)分析，可分析設(shè)備的非線性響應(yīng)。下面對(duì)時(shí)域模擬法與 DDAM 法的算法進(jìn)行解析，瞬態(tài)動(dòng)力學(xué)整體的通用方程如下：

其中，[M]為系統(tǒng)質(zhì)量矩陣；[C]為系統(tǒng)阻尼矩陣；[K]為系統(tǒng)剛度矩陣；分別為系統(tǒng)的彈性力、阻尼力、慣性力；分別為系統(tǒng)的位移、速度、加速度時(shí)間函數(shù)。

對(duì)于DDAM法，線彈性元件系統(tǒng)中，[M]，[C]為常數(shù)；而對(duì)于非線性元器件系統(tǒng)，其真實(shí)的[M]，[C]為時(shí)變矩陣。DDAM的模態(tài)疊加法，是將式（1）由時(shí)域轉(zhuǎn)換為頻域，對(duì){x(t)}做模態(tài)坐標(biāo)轉(zhuǎn)換，在模態(tài)空間中，對(duì)方程截取高階頻率范圍做近似降維解耦。由此可知，最后求得式（1）近似解的方法將不再適用于非線性元器件系統(tǒng)運(yùn)動(dòng)方程的求解。

對(duì)于非線性元器件運(yùn)動(dòng)方程的求解，一般采用隱式積分的newmark積分算法，具體算法如下：對(duì)式(1)任意時(shí)刻，及 t+Δt，系統(tǒng)的增量平衡方程為

由式（13）可求ti+Δt時(shí)刻的位移、速度、加速度表達(dá)式，即時(shí)間過程中任意一點(diǎn)的位移、速度、加速度，從而可以求得系統(tǒng)中任意點(diǎn)的應(yīng)力應(yīng)變。

Δt是時(shí)間步長(zhǎng)積分，也是決定求解精度的決定條件，Δt的步長(zhǎng)取值大小要根據(jù)具體的分析對(duì)象的輸入譜源來具體分析。

綜上所述：當(dāng)設(shè)備主要破壞原因?yàn)榈碗A模態(tài)時(shí)一般采用靜態(tài)等效法；動(dòng)態(tài)設(shè)計(jì)法適用于設(shè)備自身為線彈性材料，而對(duì)于強(qiáng)沖擊環(huán)境的設(shè)備的非線性形變計(jì)算無法精確計(jì)算；相比于動(dòng)態(tài)設(shè)計(jì)法，時(shí)域模擬法方法考慮了材料的非線性問題，可以提供更詳細(xì)的結(jié)果，更精確的外部載荷，對(duì)分析設(shè)備的結(jié)構(gòu)屬性更靈活。

2 抗沖擊設(shè)計(jì)

2.1 柜體結(jié)構(gòu)設(shè)計(jì)

目前，船用配電設(shè)備框架結(jié)構(gòu)主要采用鋼板折彎焊接成形，主體結(jié)構(gòu)分為圍板、頂板和底板三大部分組成，鋼板一般采用2～3 mm優(yōu)質(zhì)鋼板組成，保證有足夠的機(jī)械強(qiáng)度，能達(dá)到抗沖擊、振動(dòng)的船用機(jī)械使用環(huán)境。為充分了解不同型材結(jié)構(gòu)在不同情況下的受力情況，保持型材同等長(zhǎng)度且厚度均為2 mm條件下，將結(jié)構(gòu)一端固定，分別對(duì)型材三個(gè)方向施加F=500 N的正應(yīng)力，如圖1所示。

圖1 不同型材結(jié)構(gòu)示意圖

不同型材結(jié)構(gòu)產(chǎn)生的最大應(yīng)力及應(yīng)變情況見表1，由表1可知：

1）Fx方向：結(jié)構(gòu)Fx方向受力時(shí)，結(jié)構(gòu)A產(chǎn)生的應(yīng)力和應(yīng)變更小，產(chǎn)生的應(yīng)力是其他結(jié)構(gòu)的1/3左右，結(jié)構(gòu)B和結(jié)構(gòu)C受Fx方向產(chǎn)生的應(yīng)力、應(yīng)變基本相同，結(jié)構(gòu)D產(chǎn)生的應(yīng)力及應(yīng)變最大；

2）Fy方向：Fy方向受力時(shí)，結(jié)構(gòu)A和結(jié)構(gòu)B產(chǎn)生的應(yīng)力和應(yīng)變更小，產(chǎn)生的應(yīng)力是結(jié)構(gòu)C的1/3左右，結(jié)構(gòu)D產(chǎn)生的應(yīng)力及應(yīng)變最大；

3）F拉/壓方向：結(jié)構(gòu)A產(chǎn)生的應(yīng)力及應(yīng)變最小，結(jié)構(gòu)C產(chǎn)生的應(yīng)力及應(yīng)變最大。

表1 不同型材結(jié)構(gòu)最大應(yīng)力、應(yīng)變值

綜上所述，結(jié)構(gòu)A在上述幾種結(jié)構(gòu)中受不同方向產(chǎn)生的應(yīng)力、應(yīng)變值相對(duì)較小，故此結(jié)構(gòu)較其他結(jié)構(gòu)更穩(wěn)定。

2.2 有限元建模

以結(jié)構(gòu) A 為基礎(chǔ)建立的船用配電板Solidworks三維模型見圖 2(a)，由于本次計(jì)算主要考核配電板柜體結(jié)構(gòu)強(qiáng)度是否滿足垂向沖擊環(huán)境要求，不用于考核柜體內(nèi)部元器件的結(jié)構(gòu)強(qiáng)度，故將內(nèi)部元器件簡(jiǎn)化為一個(gè)質(zhì)量點(diǎn)來代替，僅保留柜體的外部框架，并將該質(zhì)量點(diǎn)的自由度與柜體底部表面節(jié)點(diǎn)進(jìn)行耦合，有限元模型見圖2(b)。

圖2 配電板建模

2.3 隔振器布置

船用配電板鋼絲繩隔振器選型以抗沖為主要目的時(shí)，需根據(jù)系統(tǒng)輸入環(huán)境要求的嚴(yán)酷度確保系統(tǒng)有足夠的“能容”，且系統(tǒng)峰值響應(yīng)頻率設(shè)計(jì)低頻范圍內(nèi)，系統(tǒng)峰值相應(yīng)頻率可按如下公式進(jìn)行計(jì)算。

式中：fn-系統(tǒng)峰值相應(yīng)頻率，Hz；K-系統(tǒng)靜剛度，N/m；a-動(dòng)剛度系數(shù)；m-系統(tǒng)承載質(zhì)量，kg。

根據(jù)以上選型原則，鋼絲繩隔振器初步設(shè)計(jì)方案如下：

1）三屏配電板柜體頂部各放置一個(gè)鋼絲繩隔振器；

2）轉(zhuǎn)角屏兩側(cè)配電板在柜門及背板底部垂向各布置4個(gè)鋼絲繩隔振器；

3）考慮轉(zhuǎn)角屏內(nèi)部只含主母排，無其他元器件，在底部外直角側(cè)分別布置一個(gè)鋼絲繩隔振器。

2.4 邊界條件

船用配電板抗沖擊計(jì)算初始邊界條件如下：

1）內(nèi)部元器件簡(jiǎn)化為一個(gè)質(zhì)量點(diǎn)，并將該質(zhì)量點(diǎn)的自由度與柜體底部表面節(jié)點(diǎn)進(jìn)行耦合。

2）柜體材料為Q235，其材料參數(shù)見表2。

3）根據(jù)設(shè)備實(shí)際使用環(huán)境，設(shè)備垂向沖擊譜見表3。

表2 材料參數(shù)

表3 設(shè)備設(shè)計(jì)沖擊譜

4）鋼絲繩隔振器參數(shù)采用彈簧單元進(jìn)行模擬替代，其剛度為466 N/m，阻尼為16.5。

5）配電板通過彈簧單元與底部基座相連且底邊默認(rèn)為固定邊界條件。

2.5 數(shù)值計(jì)算結(jié)果

船用配電板在初始設(shè)計(jì)布置方案下垂向沖擊仿真計(jì)算結(jié)果見圖3。由仿真結(jié)果可知：

1）柜體應(yīng)力較大區(qū)域主要在柜體底部，且主要集中在柜體與鋼絲繩隔振器連接處；

2）最大應(yīng)力處發(fā)生在轉(zhuǎn)角屏柜體底部A處，其應(yīng)力為 221 MPa,雖未超過材料許用應(yīng)力 235 MPa，但柜體底部框架出現(xiàn)應(yīng)力集中超過強(qiáng)度的問題，柜體梁可能會(huì)產(chǎn)生一定變形。

3）沖擊載荷下柜體位移響應(yīng)差異較大。轉(zhuǎn)角屏兩側(cè)外部框架底部位移響應(yīng)約為（5～6 mm）為其余部分的響應(yīng)（2～3 mm）的2倍以上。

4）該方案中的鋼絲繩隔振器分布未達(dá)到最佳優(yōu)化狀態(tài)，可以對(duì)底部鋼絲繩隔振器的分布方案進(jìn)行進(jìn)一步的優(yōu)化設(shè)計(jì)，以降低在沖擊載荷作用下的柜體應(yīng)力及位移響應(yīng)。

基于上述仿真計(jì)算結(jié)果，對(duì)底部隔振器布置進(jìn)行優(yōu)化，即在發(fā)生應(yīng)力最大底部框架A處加裝一個(gè)垂向隔振器；同時(shí)考慮柜體位移響應(yīng)差異較大，在保證隔振器數(shù)量不變情況下將原方案設(shè)計(jì)中頂部隔振器改為背部隔振器布置，保持其他邊界條件不變，仿真計(jì)算結(jié)果見圖4、圖5。

由仿真計(jì)算結(jié)果可知：

1）轉(zhuǎn)角屏底部框架 A處加裝隔振器后，應(yīng)力響應(yīng)最大值由221 MPa變?yōu)?14 MPa,應(yīng)變值由6.2 mm變?yōu)?.6 mm，最大值變化不大且應(yīng)力、應(yīng)變分布與原方案基本相同。

2）將頂部隔振器變更為背部隔振器布置后，同等隔振器數(shù)量布置條件下，應(yīng)力響應(yīng)最大值由213.7 MPa變?yōu)?93.3 MPa，較優(yōu)化前降低了10%左右；應(yīng)變響應(yīng)最大值由5.6 mm降低為5.3 mm，較優(yōu)化前降低了5%左右。

在進(jìn)行仿真確定模型滿足沖擊要求后，制作了產(chǎn)品樣機(jī)并在國家標(biāo)準(zhǔn)要求下對(duì)配電板進(jìn)行了沖擊試驗(yàn)，試后設(shè)備無機(jī)械損傷現(xiàn)象，滿足了沖擊試驗(yàn)要求。

圖3 垂向沖擊下應(yīng)力應(yīng)變分布圖

圖4 加裝隔振器后應(yīng)力應(yīng)變分布圖

圖5 背部隔振器布置方案應(yīng)力應(yīng)變分布圖

4 結(jié)論

通過對(duì)現(xiàn)有的抗沖擊計(jì)算方法進(jìn)行了分析，并對(duì)某型號(hào)船用配電板沖擊響應(yīng)進(jìn)行了優(yōu)化，得到結(jié)論如下：

1）時(shí)域模擬法由于考慮了材料的非線性問題，對(duì)分析設(shè)備的機(jī)構(gòu)屬性更靈活，可作為船用設(shè)備抗沖擊性能評(píng)估的一種分析方法。

2）對(duì)船用配電板研制階段進(jìn)行抗沖擊能力分析及優(yōu)化，在同等隔振器布置情況下，應(yīng)力響應(yīng)較優(yōu)化前降低了10%左右，降低了產(chǎn)品在沖擊環(huán)境下的沖擊響應(yīng)，提高了產(chǎn)品可靠性。

3）產(chǎn)品樣機(jī)在國家標(biāo)準(zhǔn)要求的沖擊環(huán)境下進(jìn)行了試驗(yàn)，試驗(yàn)結(jié)果滿足使用要求，為后續(xù)同類產(chǎn)品設(shè)計(jì)及優(yōu)化提供了理論支撐與設(shè)計(jì)參考。