殷彩玉 ，金澤宇 ，諶勇 ，華宏星

1上海交通大學機械系統(tǒng)與振動國家重點實驗室，上海200240

2上海交通大學振動沖擊、噪聲研究所，上海200240

0 引 言

潛艇是關系國家安全的重要威懾與反擊力量，也是敵方防范和打擊的重點目標，因此提高潛艇抗水下爆炸沖擊的能力具有非常重要的意義。遠場水下爆炸將主要破壞潛艇的內(nèi)部設備，水下爆炸沖擊波傳遞到雙層殼體潛艇內(nèi)部設備的主要傳播路徑為：爆炸源→水→輕外殼→舷間水→耐壓殼→設備基座→設備。因此，可以從2個方面著手提升潛艇及艦載設備的抗沖擊能力：一是通過改進或加強潛艇殼體結(jié)構、在設備上加裝性能更好的隔振器來提升潛艇殼體的強度及艦載設備的抗沖擊性能，但在潛艇空間和排水量有限的條件下，通過隔振器來進一步提高艦載設備的抗沖擊性能已難以再取得較大的突破；二是降低沖擊載荷，即通過某種手段（例如在耐壓殼體上敷設一層耐壓抗沖覆蓋層）來降低傳遞到潛艇殼體上的沖擊載荷，從而降低殼體形變及由殼體傳遞到內(nèi)部設備上的載荷。

之所以選擇在潛艇耐壓殼而不是在輕外殼上敷設抗沖覆蓋層，有以下2個原因：一是輕外殼較輕薄，即使敷設抗沖覆蓋層，也會由于輕外殼易變形而導致抗沖覆蓋層的變形有限，最終沖擊波還是很容易傳遞到潛艇內(nèi)部的耐壓殼體上；二是耐壓殼體較厚且不容易變形，在沖擊載荷作用下敷設在耐壓殼體上的抗沖覆蓋層的形變較大，可以有效地吸收沖擊波能量，從而降低傳遞到耐壓殼體上的沖擊載荷。

目前，艦艇抗沖覆蓋層的研究工作主要局限于水面艦船，國內(nèi)學者已開展了大量的理論、仿真及試驗研究，提出在船體濕表面上敷設超彈性多孔橡膠覆蓋層可以緩沖和吸收水下爆炸的沖擊能量。諶勇等［1］針對粘貼于船模上的橡膠夾芯覆蓋層在水下爆炸時的防護作用開展了試驗研究，結(jié)果表明，橡膠夾芯覆蓋層可以有效降低結(jié)構在沖擊波作用下的響應峰值，該覆蓋層相當于一個低通濾波器，對高頻信號的衰減明顯，但對低頻信號的影響不大。章振華等［2-3］采用試驗方法研究了橡膠多孔單胞結(jié)構在低速沖擊作用下的動態(tài)響應，運用Abaqus有限元軟件分析了橡膠多孔覆蓋層連續(xù)體模型對艦船的沖擊防護效果，并與試驗結(jié)果進行了對比分析，結(jié)果表明，橡膠多孔覆蓋層可以降低船體的動態(tài)響應，其防護效果較好。肖鋒等［4-7］采用Abaqus有限元軟件研究了超彈性橡膠分層圓孔蜂窩覆蓋層和六邊形橡膠蜂窩并聯(lián)歐拉桿結(jié)構的動態(tài)壓縮行為，還分析了分層圓孔結(jié)構和橡膠蜂窩結(jié)構在水下爆炸載荷作用下的沖擊防護性能。劉建湖等［8］論述了國內(nèi)外艦艇抗爆抗沖擊技術的現(xiàn)狀與發(fā)展趨勢，分析了制約我國該項技術發(fā)展的瓶頸問題，提出了增強艦艇抗爆抗沖擊能力的建議。

上述分析的橡膠多孔抗沖覆蓋層是針對水面艦船而設計的，然而潛艇在水下作業(yè)時需承受較高的靜水壓力，如果將該類空腔較大的覆蓋層直接應用于潛艇，其覆蓋層空腔在一定的靜水壓力下必然會被壓縮且發(fā)生較大的變形。由于覆蓋層發(fā)揮抗沖擊性能的一個重要機理是基于空腔的變形吸能，故當空腔壓縮變形后將嚴重影響覆蓋層的抗沖擊性能。因此，上述的多孔橡膠抗沖覆蓋層并不適用于潛艇。本文擬針對潛艇抗沖擊環(huán)境的特殊性，提出其適用的耐壓抗沖覆蓋層，分析其耐靜壓特性及抗沖擊特性，用以為潛艇的抗沖擊研究提供技術支撐。

1 耐壓抗沖覆蓋層結(jié)構形式

耐壓抗沖覆蓋層在潛艇上的粘貼方式如圖1所示。在水下爆炸載荷作用下，潛艇將同時受到靜水壓力和沖擊波的作用。因此，潛艇耐壓抗沖覆蓋層需要滿足3個條件：1）耐高靜水壓力；2）隔沖耗能；3）降低傳遞到潛艇耐壓殼體的應力及能量。滿足上述條件的結(jié)構的力學特性如圖2所示，即覆蓋層具有較高的彈性模量，可以承受靜壓且變形較??；應力平臺階段可以吸收沖擊波能量。圖2中：σ為應力，MPa；?為應變。

具備此類力學特性的結(jié)構形式主要有：超彈性多孔橡膠結(jié)構；塑性泡沫結(jié)構；三明治夾層板結(jié)構等。由于三明治夾層板在后屈曲行為中的軟化效應很大，故會限制其吸能效果。因此，本文將主要分析多孔橡膠結(jié)構和塑性泡沫結(jié)構。

1.1 橡膠多孔覆蓋層

針對多孔橡膠覆蓋層，根據(jù)其變形方式可以分為2種：共面變形結(jié)構和異面變形結(jié)構。異面變形結(jié)構主要以蜂窩結(jié)構為主，例如六邊形蜂窩結(jié)構、四邊形蜂窩結(jié)構和三角形蜂窩結(jié)構，如圖3（a）所示。共面變形典型結(jié)構包括手性結(jié)構、蜂窩結(jié)構、圓孔結(jié)構和隧道結(jié)構等，如圖3（b）所示。針對圖3所示的橡膠多孔結(jié)構，取其孔隙率為43.4%，并以1/s的速率壓縮后即可得到異面變形結(jié)構和共面變形結(jié)構的典型變形形式，如圖4所示。不同橡膠多孔結(jié)構的應力應變曲線如圖5所示。由圖4和圖5可知，橡膠多孔抗沖覆蓋層的特點如下：

1）線彈性階段的彈性模量較?。?0 MPa左右），不會超過橡膠基體的彈性模量。

2）應力平臺較低。

3）應力平臺較短。

因此，如果將耐靜壓能力較低的水面艦艇抗沖覆蓋層應用于高靜水壓環(huán)境下的潛艇，其抗沖覆蓋層將產(chǎn)生比較大的變形。在沖擊載荷作用下，其變形空間是有限的，故其抗沖擊性能并不理想。

1.2 塑性泡沫材料

塑性泡沫材料具有密度低、強度高、剛度高、造價低、可設計性強的特點，且該類材料很容易加工成復雜的形狀，已被廣泛應用于抗沖擊設計領域［9-11］。在宏觀力學方面，泡沫材料具有典型的三階段特性：線彈性階段、應力平臺階段及密實化階段。在線彈性階段，應變一般小于0.05，且彈性模量較高；在應力平臺階段，可以發(fā)生較大的塑性變形并吸收相當大的能量，同時其傳遞到結(jié)構的應力較低。由此可知，塑性泡沫結(jié)構在彈性階段的高剛度特性為其耐靜壓能力提供了保障。

1.3 耐壓抗沖覆蓋層

設計潛艇耐壓抗沖擊覆蓋層時，應考慮耐靜水壓力和抗沖擊這2個方面的能力。塑性泡沫材料有較高的比剛度，可以耐受靜水壓力且變形很小。同時，塑性泡沫材料在屈服以后有很寬的應力平臺，其吸能效果較好。考慮覆蓋層的實際粘貼情況，本文將提出潛艇耐壓抗沖覆蓋層的結(jié)構形式，即塑性泡沫與橡膠的復合覆蓋層，其可能的結(jié)構形式如圖6所示。

2 數(shù)值計算模型的有效性驗證

本文擬采用Abaqus/Explicit有限元軟件分析抗沖覆蓋層在水下爆炸載荷下的沖擊防護問題，首先將選用Schiffer等［12］的激波管實驗對數(shù)值計算方法的有效性加以驗證。實驗裝置如圖7所示，覆蓋層由鋁面板和泡沫組成，其中pst為靜水壓力。前、后面板的速度及位移的實驗結(jié)果與仿真結(jié)果對比如圖8所示，其中t為時間。由圖8可知，仿真結(jié)果與實驗結(jié)果吻合較好，驗證了仿真模型及數(shù)值計算方法的有效性。

3 計算模型

3.1 幾何模型

如圖9所示，計算模型包含2個部分：耐壓抗沖覆蓋層和水域，其中潛艇結(jié)構假設為固定邊界。耐壓抗沖覆蓋層由隧道式多孔橡膠覆蓋層和剛性聚氨酯泡沫復合而成，其具體尺寸如圖7所示，模型中水域的長度為2 m。在有限元軟件中，將采用4節(jié)點的平面應變單元CPE4R來模擬耐壓抗沖覆蓋層，采用4節(jié)點聲學單元AC2D4R來模擬水，采用邊界條件Tie來約束泡沫與橡膠、橡膠與水，覆蓋層設周期邊界。

3.2 材料屬性

本文采用Hyperelastic材料模型來模擬橡膠，其單軸應力應變曲線如圖10（a）所示（拉伸為正），泊松比取為0.485。采用Crushable foam材料模型來模擬剛性聚氨酯泡沫材料，其壓縮應力應變曲線如圖10（b）所示：在剛塑性模型中，其屈服強度為1.280 6 MPa，密實化應變?yōu)?.549 1；在線彈性階段，其彈性模量和泊松比分別為30 MPa和0.3；在塑性階段，采用各向同性硬化模型，塑性泊松比設為0，Compression Yield Stress Ratio設為1.732。

水的密度和聲速分別為1000kg/m3和1500 m/s，其體積模量為2.25×109Pa。為了模擬水的空化效應，本文采用Total Wave進行計算，其中水的空化閾值設為0。

3.3 沖擊環(huán)境參數(shù)

水下爆炸沖擊波可以表示為指數(shù)衰減波［13］：

式中：pin為水下爆炸沖擊波；p0為沖擊波峰值；θ為衰減系數(shù)，在實際爆炸中與炸藥當量和爆距有關。

本文主要研究耐壓抗沖覆蓋層在靜水壓力和沖擊載荷共同作用下的沖擊防護效果，并未涉及實艇計算，故假設以下2種不同強度的沖擊波：

1）沖擊波1：p0=10 MPa，θ=1 ms。

2）沖擊波2：p0=15 MPa，θ=1 ms。

考慮4種靜水壓力，即0.5，1，1.5，2 MPa。

4 耐壓抗沖覆蓋層的沖擊防護效果分析

4.1 耐靜壓特性

圖11所示為覆蓋層在不同靜水壓力作用下的變形圖，其中S為覆蓋層在靜水壓力作用下的Mises應力。由圖11可知，當靜水壓力超過覆蓋層的屈服強度后，覆蓋層將出現(xiàn)明顯變形；而當靜水壓力小于覆蓋層的屈服強度時，覆蓋層的變形幾乎可以忽略。因此，可以根據(jù)實際設計中的耐靜壓指標，選取不同屈服強度的覆蓋層。

4.2 抗沖擊特性

圖12所示為耐壓抗沖覆蓋層在2種沖擊波作用下的變形情況，其中PEEQ為覆蓋層的等效塑性變形。由圖12可知，沖擊波越大，等效塑性變形就越大。對于剛性聚氨酯泡沫，其靠近橡膠部分的變形要小于遠離橡膠部分，這是因為橡膠幾乎不可壓縮，從而阻礙了其周圍泡沫的變形。

由于濕表面各節(jié)點的壓力基本相同，故將濕表面各節(jié)點的壓力求和并取平均值，即可得到濕表面的平均壓力pwet，如圖13所示。由圖13可知，在峰值為10 MPa，衰減系數(shù)為1 ms的入射指數(shù)衰減沖擊波作用下，濕表面壓力迅速降低，然后基本保持為應力平臺，該應力平臺的均值與泡沫材料的屈服強度相當；在峰值為15 MPa，衰減系數(shù)為1 ms的入射指數(shù)衰減沖擊波作用下，濕表面壓力在應力平臺階段迅速上升，這主要是因為沖擊波強度較大，覆蓋層的厚度不足以吸收全部的沖擊波能量，導致覆蓋層發(fā)生了密實化。此外，敷設覆蓋層后，濕表面壓力形態(tài)與入射指數(shù)衰減沖擊波的差別很大，即由峰值高、脈寬短的指數(shù)衰減波變?yōu)榉逯档?、持續(xù)時間長的壓力波。

圖14所示為傳遞到潛艇殼體的平均支反力Fr。由圖14可知，在峰值為10 MPa，衰減系數(shù)為1 ms的入射指數(shù)衰減沖擊波作用下，當覆蓋層厚度足以吸收沖擊波能量時，傳遞到殼體的應力由覆蓋層的靜態(tài)屈服強度決定；在峰值為15 MPa，衰減系數(shù)為1 ms的入射指數(shù)衰減沖擊波作用下，當覆蓋層密實化以后，傳遞到殼體的應力也會急劇上升。因此，應避免覆蓋層發(fā)生完全密實化。與入射指數(shù)衰減沖擊波峰值相比，傳遞到殼體的壓力峰值要遠低于沖擊波峰值，這說明耐壓抗沖覆蓋層可以有效地降低傳遞到結(jié)構的沖擊波強度。

5 結(jié) 語

本文針對潛艇所處的特殊沖擊環(huán)境，提出了適用于潛艇的耐壓抗沖覆蓋結(jié)構形式，即塑性泡沫與橡膠的復合耐壓抗沖覆蓋層，并驗證了其抵抗靜水壓力和沖擊載荷的能力，可為潛艇抗沖覆蓋層的設計提供參考。但是，需要注意2個問題：一是如果覆蓋層的屈服強度小于靜水壓力，則覆蓋層在靜水壓力作用下將壓潰，因此在設計中覆蓋層的屈服強度要大于工作環(huán)境的靜水壓力；二是如果沖擊波強度過大，覆蓋層就會發(fā)生密實化，進而影響其抗沖擊效果，故在實際設計中應避免覆蓋層密實化現(xiàn)象的發(fā)生。

［1］諶勇，華宏星，汪玉，等.超彈性夾芯覆蓋層的水下爆炸防護性能［J］. 爆炸與沖擊，2009，29（4）：395-400.CHEN Y，HUA H X，WANG Y，et al.Protective effects of hyper-elastic sandwiches coated onto metal boxes subjected to underwater explosion［J］.Explosion and Shock Waves，2009，29（4）：395-400（in Chinese）.

［2］章振華，諶勇，華宏星，等.超彈性夾芯覆蓋層抗沖擊性能分析及實驗研究［J］.振動與沖擊，2012，31（5）：132-134.ZHANG Z H，CHEN Y，HUA H X，et al.Simulation and experimental research of the hyperelastic sandwich coatings in crush dynamics［J］.Journal of Vibration and Shock，2012，31（5）：132-134（in Chinese）.

［3］章振華，諶勇，肖鋒，等.敷設超彈性覆蓋層艦船水下爆炸沖擊實驗與仿真分析［J］.振動與沖擊，2014，33（10）：106-112.ZHANG Z H，CHEN Y，XIAO F，et al.Underwater explosion research on ships with hyper-elastic coating［J］.Journal of Vibration and Shock，2014，33（10）：106-112（in Chinese）.

［4］肖鋒，諶勇，孫靖雅，等.超彈性橡膠材料分層圓孔蜂窩防護覆蓋層動態(tài)壓縮行為及性能研究［J］.振動與沖擊，2013，32（22）：13-20.XIAO F，CHEN Y，SUN J Y，et al.Dynamic crush behavior and performance of layeredhoneycomb hyperelastic rubber claddings［J］.Journal of Vibration and Shock，2013，32（22）：13-20（in Chinese）.

［5］肖鋒，華宏星，諶勇，等.潛艇濕表面抗沖覆蓋層壓縮特性及抗沖擊性能研究［J］.振動與沖擊，2013，32（18）：126-132.XIAO F，HUA H X，CHEN Y，et al.Performance analysis on compression and shock resistance of anti-shock layer coated onto submarine wetted surface［J］.Journal of Vibration and Shock，2013，32（18）：126-132（in Chinese）.

［6］肖鋒，華宏星，諶勇，等.設計參數(shù)對手性蜂窩橡膠覆蓋層水下爆炸抗沖擊性能的影響［J］.振動與沖擊，2014，33（1）：56-62.XIAO F，HUA H X，CHEN Y，et al.The influence of design parameters on underwater explosion shock resistance of chiral honeycomb rubber cladding［J］.Journal of Vibration and Shock，2014，33（1）：56-62（in Chinese）.

［7］肖鋒，諶勇，朱大巍，等.水下非接觸性爆炸下分層圓孔覆蓋層動響應及抗沖擊性能研究［J］.振動與沖擊，2013，32（19）：111-118.XIAO F，CHEN Y，ZHU D W，et al.Dynamic response and shock resistance of layered claddings subjected to a non-contact underwater explosion［J］.Journal of Vibration and Shock，2013，32（19）：111-118（in Chinese）.

［8］劉建湖，周心桃，潘建強，等.艦艇抗爆抗沖擊技術現(xiàn)狀和發(fā)展途徑［J］.中國艦船研究，2016，11（1）：46-56，71.LIU J H，ZHOU X T，PAN J Q，et al.The state analysis and technical development routes for the anti-explosion and shock technology of naval ships［J］.Chinese Journal of Ship Research，2016，11（1）：46-56，71（in Chinese）.

［9］FLECK N A，DESHPANDE V S.The resistance of clamped sandwich beams to shock loading［J］.Journal of Applied Mechanics，2004，71（3）：386-401.

［10］TILBROOK M T，DESHPANDE V S，F(xiàn)LECK N A.The impulsive response of sandwich beams：analytical and numerical investigation of regimes of behavior［J］.Journal of the Mechanics and Physics of Solids，2006，54（11）：2242-2280.

［11］MCSHANE G J，DESHPANDE V S，F(xiàn)LECK N A.Underwater blast response of free-standing sandwich plates with metallic lattice cores［J］.International Journal of Impact Engineering，2010，37（11）：1138-1149.

［12］SCHIFFER A， TAGARIELLIV L.One-dimensional response of sandwich plates to underwater blast：fluid-structure interaction experiments and simulations［J］.International Journal of Impact Engineering，2014，71：34-49.

［13］COLE R H.Underwater explosions［M］.Princeton，NJ，USA：Princeton University Press，1948.