孫龍泉，王都亮，李志鵬，劉登科

(哈爾濱工程大學(xué) 船舶工程學(xué)院，哈爾濱 150009)

結(jié)構(gòu)物入水，如飛船返回艙入水[1]、飛機(jī)艙段入水[2]、空投魚雷入水[3]等，都是涉及多相流動(dòng)、自由液面和動(dòng)邊界多場(chǎng)耦合的復(fù)雜過(guò)程，其作用時(shí)間短暫，參數(shù)變化劇烈。尤其航行體高速入水時(shí)，結(jié)構(gòu)會(huì)遭遇到強(qiáng)烈的瞬時(shí)沖擊壓力及過(guò)載，如果不進(jìn)行緩沖降載，容易發(fā)生結(jié)構(gòu)的損壞及內(nèi)部器件的失靈，造成無(wú)法挽回的損失。因此研究高速入水降載技術(shù)意義重大。

在入水降載方面，Hinkley[4]首先提出在航行體頭部加裝頭帽來(lái)保護(hù)航行體結(jié)構(gòu)不受入水沖擊載荷影響過(guò)大，入水沖擊力使頭罩破壞解體，不影響航行體本體正常工作。最初頭帽為一個(gè)整體，經(jīng)過(guò)發(fā)展改進(jìn)，頭帽分為外罩和緩沖材料兩個(gè)部分，外罩起到在空中飛行階段整流作用，緩沖材料起到隔沖降載、緩沖吸能、隔絕碎片的作用。在后續(xù)對(duì)頭帽的研究中，宣建明等[5]通過(guò)入水試驗(yàn)，采用復(fù)合材料作為外罩的主體材料，研究了外罩的破壞形式。雷江濤等[6]通過(guò)有限元分析的方法，研究了頭帽材料、壁厚和緩沖墊對(duì)航行器頭帽分離可靠性的影響。徐新棟等[7]提出了兩種使頭帽入水后更易破裂的改進(jìn)方法，一種為降低頭帽和雷體之間的摩擦力，另一種為降低整流罩強(qiáng)度。并通過(guò)仿真計(jì)算的方式進(jìn)行了驗(yàn)證。錢立新等[8]對(duì)帶頭帽的魚雷的入水過(guò)程進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)驗(yàn)證，得出“撐進(jìn)破壞”是帶頭罩魚雷入水時(shí)頭罩破壞的有效方式。預(yù)制槽弱鏈結(jié)構(gòu)是實(shí)現(xiàn)“撐進(jìn)破壞“方式的有效結(jié)構(gòu)形式的結(jié)論。這些研究主要針對(duì)外罩的破壞模式，為外罩結(jié)構(gòu)的設(shè)計(jì)提供理論基礎(chǔ)，但對(duì)內(nèi)部降載元件的設(shè)計(jì)及高速入水條件下的優(yōu)化未做過(guò)多涉及。錢立新等[9]通過(guò)分析空氣動(dòng)力載荷和入水沖擊載荷，提出了一種用切線法設(shè)計(jì)正切尖拱魚雷頭罩的可行域的方法。該研究對(duì)外罩結(jié)構(gòu)尺寸的設(shè)計(jì)提供了可行參考。而針對(duì)內(nèi)部降載元件的研究中，路龍龍[10]使用有限元軟件ANSYS/LS-DYNA 計(jì)算了小角度入水和垂直入水時(shí)受到的入水載荷，提出了入水緩沖保護(hù)頭帽的設(shè)計(jì)方案，簡(jiǎn)單分析了緩沖材料的質(zhì)量和結(jié)構(gòu)要求。王永虎等[11]采用入水沖擊加載條件下的本構(gòu)關(guān)系式，計(jì)算了帶緩沖頭帽入水的動(dòng)態(tài)緩沖性能。討論了不同入水速度、不同密度泡沫對(duì)緩沖效果的影響。王永虎[12]又引入硬質(zhì)聚氨酯泡沫塑料(rigid polyurethane foams,RPUF)在沖擊條件下的本構(gòu)關(guān)系，并進(jìn)行了相應(yīng)修正，求出入水沖擊響應(yīng)本構(gòu)關(guān)系式，得出了快速確定RPUF控制參數(shù)的方法。這些研究中未涉及到高速入水，并且工況較為單一，所做的工作仍不夠充分，結(jié)論無(wú)法支持起入水工程應(yīng)用研究。

基于上述原因，本文使用歐拉-拉格朗日耦合算法(CEL方法)建立航行體-空氣-水流固耦合模型，分析在不同入水速度(50～150 m/s)下，多種泡沫鋁材料在入水過(guò)程中的降載效果，探討入水速度和最優(yōu)泡沫鋁參數(shù)之間的規(guī)律，為入水沖擊結(jié)構(gòu)降載設(shè)計(jì)提供參考。

1 模型建立

為了降低航行體入水沖擊載荷，在航行體的頭部加上緩沖材料，在航行體入水砰擊到水面的瞬間，依靠緩沖材料壓縮吸收能量，從而起到緩沖作用，降低砰擊載荷對(duì)航行體的影響。

1.1 物理模型設(shè)計(jì)

泡沫鋁是一種功能性很強(qiáng)的降載吸能材料，具有密度小、吸收沖擊能力強(qiáng)、易加工、成形精度高等優(yōu)點(diǎn)可用于抗沖擊結(jié)構(gòu)設(shè)計(jì)。相對(duì)密度μ(μ=ρ*/ρs，ρ*為泡沫鋁表現(xiàn)密度，ρs為致密鋁密度)范圍一般為2%～60%，可由發(fā)泡工藝過(guò)程控制。帶緩沖材料的航行體入水初始狀態(tài)如圖1所示。

圖1 帶緩沖材料的航行體入水初始狀態(tài)Fig.1 Initial state of water entry for a vehicle with buffer material

泡沫鋁在壓縮時(shí)，其應(yīng)力應(yīng)變曲線呈現(xiàn)明顯的三階段特征，即彈性變形階段，塑性屈服平臺(tái)階段，壓縮密實(shí)化階段。泡沫鋁典型的應(yīng)力應(yīng)變曲線如圖2所示。

圖2 泡沫鋁的壓縮應(yīng)力應(yīng)變曲線Fig.2 Compressive stress-strain curve of foamed aluminum

(1)

此方法計(jì)算出泡沫鋁(μ=12.7%)的平臺(tái)應(yīng)力為4 MPa，εd為0.75，與文獻(xiàn)[16]中的實(shí)驗(yàn)數(shù)據(jù)吻合較好。根據(jù)其擬合的經(jīng)驗(yàn)公式，擬合出各相對(duì)密度下的應(yīng)力應(yīng)變曲線，如圖 3所示。

圖3 不同相對(duì)密度泡沫鋁的擬合應(yīng)力應(yīng)變曲線Fig.3 Fitting stress-strain curves of foamed aluminum with different relative densities

1.2 數(shù)值模型設(shè)計(jì)

采用耦合歐拉-拉格朗日方法，進(jìn)行流固耦合求解，航行體和泡沫鋁采用拉格朗日單元，空氣域和水域采用歐拉單元。

選取MK46魚雷作為需要降載的航行體，其長(zhǎng)度為2.67 m，直徑為0.324 m，重量為230 kg。

圖4 MK46魚雷模型頭部形狀尺寸(m)Fig.4 Head shape and size of MK46(m)

航行體建立可變形體，采用剛體(rigid body)約束，模型為空心，網(wǎng)格采用六面體結(jié)構(gòu)化網(wǎng)格，保證航行體截面上至少有30個(gè)網(wǎng)格，網(wǎng)格尺寸設(shè)為0.02，網(wǎng)格數(shù)量為54 496，網(wǎng)格劃分如圖 5所示。材料為鋁合金7075，參數(shù)如表 1所示。

表1 航行體材料參數(shù)Tab.1 Material parameters of vehicle

圖5 航行體模型網(wǎng)格劃分Fig.5 Mesh division of vehicle model

泡沫鋁緩沖結(jié)構(gòu)形狀設(shè)計(jì)為圓臺(tái)型，頂部半徑0.12 m，底部半徑0.1 m，高度0.1 m。材料模型設(shè)置采用軟件中crushable foam塑性模型，該模型用于分析可壓碎的泡沫，通常用作能量吸收結(jié)構(gòu)。網(wǎng)格采用六面體結(jié)構(gòu)化網(wǎng)格劃分，單元類型為C3D8R，網(wǎng)格尺寸為0.01，數(shù)量為5 740。泡沫鋁緩沖降載模型的建模及網(wǎng)格劃分如圖6所示。

圖6 泡沫鋁緩沖降載模型的建模及網(wǎng)格劃分Fig.6 Mesh division of foam aluminum Buffer Load Reduction Model

歐拉域分為兩部分，上側(cè)部分為空氣域，下側(cè)部分為水域，為了提高計(jì)算效率，其網(wǎng)格在水氣交界面進(jìn)行加密，加密處網(wǎng)格采用漸進(jìn)網(wǎng)格，靠近水氣交界面處網(wǎng)格最密，最大程度的保證撞水時(shí)刻的計(jì)算精度，最密處網(wǎng)格尺寸設(shè)置0.015。權(quán)衡計(jì)算精度與計(jì)算資源消耗，結(jié)構(gòu)體與歐拉域網(wǎng)格單元尺寸之比SR=1.5～2.0[18]。歐拉網(wǎng)格區(qū)域劃分和網(wǎng)格劃分如圖7所示。

圖7 歐拉域區(qū)域劃分及網(wǎng)格劃分Fig.7 Regionalization and mesh division of Euler Domain

初始時(shí)刻，航行體處于氣水交界面上方，其裝配示意圖如圖8所示。

圖8 入水模型裝配圖Fig.8 Assembly diagram of water entry model

歐拉域分水域和空氣域，流體力學(xué)材料模型以狀態(tài)方程的形式定義，狀態(tài)方程可以確定材料的體積力。線性Hugoniot(Us-Up)狀態(tài)方程可用于模擬由納維爾-斯托克斯方程控制的黏性層流。水域是用Mie-Gruneisen狀態(tài)方程的線性Hugoniot形式來(lái)模擬[17]。

Mie-Gruneisen狀態(tài)方程通常用于模擬高壓下的材料,是關(guān)于能量的線性方程，一般形式為

p-pH=Γρ(Em-EH)

(2)

式中：pH是Hugoniot的壓力；EH為Hugoniot的比能。

Γ是Grüneisen比率，其值為Γ=Γ0ρ0/ρ，其中：Γ0為材料常數(shù)，ρ0為參考密度。EH可以用PH表示

(3)

式中：η=1-ρ0/ρ為名義體積壓縮應(yīng)變。

從上式中消去Γ和EH，得到

(4)

在每個(gè)物質(zhì)點(diǎn)上需要同時(shí)求解狀態(tài)方程和能量方程，它們表示壓力和內(nèi)能的耦合方程。

一般形式的Hugoniot數(shù)據(jù)的擬合曲線為

(5)

式中：c0和s是沖擊波速度Us和質(zhì)點(diǎn)速度Up之間關(guān)系的參數(shù)。

Us=c0+sUp

(6)

基于上述假設(shè)，線性Us-UpHugoniot式子可以寫為

(7)

表2 水的材料參數(shù)Tab.2 Material parameters of water

空氣域選取理想氣體狀態(tài)方程，其狀態(tài)方程為

p+pA=ρR(θ-θZ)

(8)

式中：pA為周圍環(huán)境壓力，R為氣體常數(shù),θ為此刻溫度，θZ為所用溫度標(biāo)尺的絕對(duì)零度。這是一種理想化的氣體模型，可以近似地模擬適當(dāng)條件下的任意氣體?？諝庥虻某跏紬l件主要指定其初始密度ρ0，初始?jí)毫0，周圍環(huán)境壓力pA，絕熱質(zhì)數(shù)γ，定容比熱容cν，空氣的材料參數(shù)如表3所示。

表3 空氣材料參數(shù)Tab.3 Material parameters of air

歐拉域底面和四周固定，頂部施加大氣壓力，空氣域內(nèi)部設(shè)置101 325 Pa均勻壓力，水域設(shè)置隨深度變化的靜水壓力。設(shè)置重力加速度方向豎直向下。接觸條件設(shè)為通用接觸，物體之間設(shè)置摩擦因數(shù)為0.02。初始時(shí)刻航行體處在自由液面上方。

1.3 有效性驗(yàn)證

為驗(yàn)證CEL方法的有效性，采用數(shù)值仿真與文獻(xiàn)[19]中圓柱體入水實(shí)驗(yàn)對(duì)比的方法，圓柱體直徑42 mm，長(zhǎng)度208 mm，質(zhì)量為0.61 kg，以14.5 m/s的速度入水。實(shí)驗(yàn)中，采用phantom VEO 340S高速攝像機(jī)獲取入水圖像，使用TEMA軟件對(duì)實(shí)驗(yàn)視頻中模型上的標(biāo)記點(diǎn)進(jìn)行追蹤，得到速度、加速度變化曲線。數(shù)值方法中，對(duì)圓柱體采用rigid body約束，提取頭部中心點(diǎn)速度及加速度變化曲線，入水過(guò)程形態(tài)、速度曲線和加速度曲線如圖9和圖10所示。

圖9 實(shí)驗(yàn)結(jié)果和數(shù)值模擬形態(tài)對(duì)比(t=1.67 ms,t=4.58 ms,t=7.50 ms,t=10.40 ms,t=13.33 ms)Fig.9 Comparison of experimental results and numerical simulation morphology

圖10 數(shù)值模擬和實(shí)驗(yàn)結(jié)果運(yùn)動(dòng)參數(shù)對(duì)比Fig.10 Comparison of motion parameters between numerical simulation and experimental results

由數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果對(duì)比中可以看出，CEL方法模擬得到的結(jié)果在運(yùn)動(dòng)參數(shù)和空泡形態(tài)都與實(shí)驗(yàn)結(jié)果吻合的比較好，在入水過(guò)程中的加速度變化趨勢(shì)基本一致，且入水瞬間加速度峰值相差在10%以內(nèi)。值得注意的是，從入水空泡發(fā)展形態(tài)、速度和加速度的變化情況看，在航行體入水“早期”，即約10 ms時(shí)間內(nèi)，CEL方法模擬結(jié)果與實(shí)驗(yàn)符合程度較好，而在入水的“后期”，由于歐拉網(wǎng)格變形逐漸加大，累積誤差不斷增加，使空泡形態(tài)逐漸偏離，而且模擬中無(wú)法考慮空泡內(nèi)部氣體的運(yùn)動(dòng)和空化特性，導(dǎo)致入水“后期”數(shù)值誤差較大。但從本文關(guān)注的入水沖擊問(wèn)題和加速度峰值出現(xiàn)的時(shí)間看，沖擊峰值發(fā)生在航行體入水的“早期”，因此，可以認(rèn)為本文所采用的CEL方法是合適的。

2 計(jì)算結(jié)果與分析

以泡沫鋁作為緩沖材料，在航行體豎直入水過(guò)程中，依靠緩沖材料壓縮吸收能量，提取航行體頭部中心處加速度變化曲線，研究在不同相對(duì)密度的泡沫鋁作為緩沖材料對(duì)降載效果的影響，進(jìn)一步的，探究不同入水速度下的降載效果。

2.1 泡沫鋁相對(duì)密度對(duì)航行體沖擊性能的影響

泡沫鋁緩沖結(jié)構(gòu)隨航行體以100 m/s的速度入水。選取相對(duì)密度為8%，12%，14%，15%，17%，19%，23%的泡沫鋁進(jìn)行計(jì)算，對(duì)入水過(guò)程進(jìn)行數(shù)值分析，驗(yàn)證其降載效果。

圖11為不同相對(duì)密度μ泡沫鋁降載過(guò)程的加速度曲線。圖12為帶μ=17%的泡沫鋁入水過(guò)程圖，分別對(duì)應(yīng)0 ms，0.625 ms，1.25 ms時(shí)刻。泡沫鋁緩沖材料隨航行體入水時(shí)，緩沖材料先接觸到水面，入水阻力作用到緩沖材料上，受力壓縮，進(jìn)入短暫的彈性變形階段，經(jīng)由緩沖材料向航行體傳遞，緩沖材料向航行體傳遞的最大應(yīng)力為泡沫鋁材料的屈服應(yīng)力極限σs，超過(guò)這個(gè)應(yīng)力值時(shí)，泡沫鋁材料即進(jìn)入塑性變形階段，將沖擊能量吸收。此時(shí)向航行體傳遞的應(yīng)力即為泡沫鋁材料的平臺(tái)應(yīng)力σp，此時(shí)航行體的加速度曲線到達(dá)第一個(gè)峰值a1，圖 11中每條加速度曲線從0.2 ms開始到0.8 ms都有一個(gè)臺(tái)階，此階段即為泡沫鋁材料的塑性應(yīng)變階段，此時(shí)航行體受到緩沖材料所傳遞的力即為泡沫鋁材料的平臺(tái)應(yīng)力值。圖13為泡沫鋁前后受壓變形對(duì)比，由于泡沫鋁材料在此過(guò)程中，不是均勻壓潰，靠近水面一側(cè)壓潰程度大于緊貼航行體頭部一側(cè)，泡沫鋁中心受壓程度明顯大于四周，所以該階段的加速度曲線呈現(xiàn)抖動(dòng)。0.625 ms時(shí)刻，航行體的受力如圖12(b)所示，由于泡沫鋁材料在塑性變形階段，緩沖材料向航行體傳遞的應(yīng)力值幾乎不變，為σp，所以在此段時(shí)間內(nèi)加速度的變化量很小。隨后緩沖材料繼續(xù)壓縮，直到泡沫鋁材料塑性變形平臺(tái)化階段結(jié)束，進(jìn)入密實(shí)化階段，應(yīng)力才會(huì)繼續(xù)上升，同時(shí)，隨著航行體入水深度的增加，1.25 ms時(shí)刻，航行體直接入水發(fā)生接觸，部分入水阻力Fd直接加載在航行體上，此時(shí)刻的受力如圖12(c)所示，航行體的加速度值在剩余能量的作用下達(dá)到第二個(gè)峰值a2，隨后隨著航行體入水深度的增加，加速度值減小，逐漸趨于穩(wěn)定。

圖11 不同μ泡沫鋁降載過(guò)程的加速度曲線Fig.11 Acceleration curves of different μ foamed aluminum during load reduction

圖12 航行體入水過(guò)程Fig.12 Water entry process of vehicle

圖13 泡沫鋁壓縮變化對(duì)比Fig.13 Comparison of compression changes of foamed aluminum

表4列出了不同μ泡沫鋁降載過(guò)程中加速度曲線的最大值，圖14所示為表6中加速度峰值變化圖。

表4 不同相對(duì)密度的泡沫鋁降載后對(duì)應(yīng)的加速度最大值Tab.4 The maximum acceleration of foam aluminum with different relative densities after load reduction

圖14 不同泡沫鋁降載后的加速度最大值變化圖Fig.14 Maximum acceleration variation diagram after different foamed aluminum load reduction

從圖11中可以看出，泡沫鋁材料的μ越大，即泡沫鋁材料塑性階段平臺(tái)應(yīng)力越大，進(jìn)入塑性變形階段時(shí)航行體的加速度值越大，即a1越大。這是由于μ越大，σp越大，緩沖材料傳遞給航行體的應(yīng)力就會(huì)越大。而隨著μ的增大，a2降低，說(shuō)明隨著μ增大，泡沫鋁的吸能效果增加。當(dāng)泡沫鋁的μ較小時(shí)，a1

通過(guò)圖14加速度最大值的對(duì)比可以看出，入水過(guò)程中航行體的最大加速度值隨著泡沫鋁材料的μ值增大先降低后升高，在μ=17%時(shí)達(dá)到最低，也就是說(shuō)μ=17%時(shí)降載效果最好。這是由于μ偏小時(shí)，a1

將μ=17%的泡沫鋁降載過(guò)程和航行體直接入水的加速度曲線進(jìn)行對(duì)比，如圖15所示。對(duì)加速度時(shí)歷曲線進(jìn)行傅里葉變換，得到的頻譜曲線如圖15所示(為突出曲線的變化，對(duì)縱坐標(biāo)取lg對(duì)數(shù))。

圖15 帶μ=17%的泡沫鋁入水和直接入水的加速度對(duì)比Fig.15 Acceleration comparison of water entry with foamed aluminum μ=17% and direct water entry

入水后，泡沫鋁材料受到入水阻力的作用，先進(jìn)入短暫的彈性變形階段，0.2 ms時(shí)刻，進(jìn)入塑性變形階段，沖擊能量通過(guò)使泡沫鋁塑性變形而吸收，起到降載作用。0.75 ms時(shí)刻，泡沫鋁到達(dá)密實(shí)化階段，繼續(xù)產(chǎn)生很小的應(yīng)變就對(duì)應(yīng)的應(yīng)力值很大，此時(shí)泡沫鋁結(jié)構(gòu)向航行體傳遞的應(yīng)力也會(huì)迅速增加，且此時(shí)航行體頭部沾濕，水開始對(duì)航行體施加阻力，隨著航行體沾濕面積在速度方向上的投影面積增大，航行體在水中受到的阻力增大，到達(dá)1.25 ms時(shí)刻，航行體沾濕面積在速度方向上的投影面積達(dá)到最大，加速度達(dá)到峰值，為1 099 m/s2，之后隨著航行體速度的降低，航行體受到水的阻力逐漸減小，到達(dá)2.5 ms左右時(shí)，加速度區(qū)域穩(wěn)定，穩(wěn)定值約為400 m/s2。航行體直接入水的加速度曲線峰值為5 566 m/s2，而帶μ=17%的泡沫鋁緩沖結(jié)構(gòu)的航行體入水加速度峰值為1 099 m/s2，加速度的最大值降低80.3%。

圖16中，在0頻率時(shí)，直接入水要高于帶泡沫鋁入水，說(shuō)明泡沫鋁在入水過(guò)程吸收了一定的能量，直接入水曲線的峰值所在的頻率要高于到泡沫鋁入水。對(duì)比時(shí)歷曲線與頻率曲線，可以確定泡沫鋁在降載過(guò)程中起到了吸能作用，并且優(yōu)化載荷，降低沖擊。

圖16 帶μ=17%的泡沫鋁入水和直接入水的加速度頻譜Fig.16 Acceleration frequency spectrum of water entry with foamed aluminum μ=17% and direct water entry

2.2 入水速度對(duì)帶緩沖材料航行體沖擊性能的影響

當(dāng)航行體的入水速度不同時(shí)，入水時(shí)，緩沖材料受到入水沖擊力的大小也會(huì)產(chǎn)生差異，降載材料的最優(yōu)參數(shù)也會(huì)發(fā)生改變。帶泡沫鋁材料入水的航行體，加速度曲線中第一個(gè)加速度峰值和第二個(gè)加速度峰值基本相等時(shí)，降載效果達(dá)到最好。

2.2.1 50 m/s入水時(shí)的參數(shù)選擇

選取μ分別為6%、7%、8%、9%，10%。對(duì)入水過(guò)程進(jìn)行數(shù)值分析。

帶不同相對(duì)μ鋁材料以50 m/s的速度入水的加速度曲線對(duì)比如圖17所示。圖19為μ=8%的泡沫鋁材料入水與直接入水的加速度對(duì)比圖。對(duì)加速度進(jìn)行傅里葉變換，得到圖20中頻譜曲線。

圖17 不同μ泡沫鋁降載過(guò)程加速度曲線對(duì)比圖(50 m/s)Fig.17 Acceleration curves of different μ foamed aluminum during load reduction(50 m/s)

圖18 不同μ泡沫鋁降載過(guò)程加速度峰值變化曲線(50 m/s)Fig.18 Maximum acceleration variation diagram after different foamed aluminum load reduction(50 m/s)

圖19 帶μ=8%的泡沫鋁入水和直接入水的加速度對(duì)比Fig.19 Acceleration comparison of water entry with foamed aluminum μ=8% and direct water entry

從圖17可以看出，當(dāng)?shù)赼1=a2時(shí),加速度的最大值達(dá)到最低。圖18為50 m/s入水時(shí)不同μ泡沫鋁降載過(guò)程加速度峰值變化曲線，加速度的最大值隨著μ的增大先降低后升高，在μ=8%時(shí)達(dá)到最低。通過(guò)對(duì)比得出，當(dāng)μ=8%時(shí)，加速度最大值最低，為261 m/s2。入水速度為50 m/s時(shí)，μ=8%的泡沫鋁材料對(duì)入水沖擊載荷削峰效果最好。圖20中，直接入水在0頻率時(shí)峰值，以及其余峰值所在頻率，都高于帶泡沫鋁入水。

圖20 帶μ=8%的泡沫鋁入水和直接入水的加速度頻譜Fig.20 Acceleration frequency spectrum of water entry with foamed aluminum μ=8% and direct water entry

2.2.2 150 m/s入水時(shí)的參數(shù)選擇

帶不同μ泡沫鋁材料以150 m/s的速度入水的加速度曲線對(duì)比如圖21所示。圖23為帶μ=29%的泡沫鋁材料入水與直接入水的加速度對(duì)比圖。圖24中，對(duì)加速度進(jìn)行傅里葉變換，得到其中頻譜曲線。

圖21 不同μ泡沫鋁降載過(guò)程加速度曲線對(duì)比圖(150 m/s)Fig.21 Acceleration curves of different μ foamed aluminum during load reduction(150 m/s)

當(dāng)?shù)赼1=a2時(shí)從圖21加速度對(duì)比圖中可以看出，當(dāng)μ分別為21%、25%、29%時(shí)，μ越大，a1越大，而a2越小，μ=33%時(shí)，a1隨著泡沫鋁的μ增加而增大，而a2沒(méi)有明顯減小，μ=37%，航行體第二個(gè)加速度峰值反而較μ=33%時(shí)的大。圖22為150 m/s入水時(shí)不同μ泡沫鋁降載過(guò)程加速度峰值變化曲線，加速度的最大值隨著μ的增大先降低后升高，在μ=29%時(shí)達(dá)到最低，為2 853 m/s2。說(shuō)明當(dāng)入水速度為150 m/s時(shí)，μ=29%泡沫鋁材料對(duì)入水沖擊載荷削峰效果最好。

圖22 不同μ泡沫鋁降載過(guò)程加速度峰值變化曲線(150 m/s)Fig.22 Maximum acceleration variation diagram after different foamed aluminum load reduction(150 m/s)

圖23 帶μ=29%的泡沫鋁入水和直接入水的加速度對(duì)比Fig.23 Acceleration comparison of water entry with foamed aluminum μ=29% and direct water entry

圖24 帶μ=29%的泡沫鋁入水和直接入水的加速度頻譜Fig.24 Acceleration frequency spectrum of water entry with foamed aluminum μ=29% and direct water entry

直接入水加速度頻譜曲線與帶泡沫鋁入水相比，在0頻率時(shí)峰值近似，曲線峰值所在頻率更高。

不同入水速度下的最優(yōu)泡沫鋁材料的μ及其降載效果如表5所示。

表5 不同入水速度下的最優(yōu)μ及降載效果Tab.5 Optimum μ and load reduction effect under different water entry speeds

圖25為不同入水速度下最優(yōu)μ，隨著入水速度的增加，降載效果最好的泡沫鋁材料的μ也隨之增大，緩沖材料的降載效果隨之降低。由前一節(jié)的分析可知，航行體達(dá)到第一個(gè)加速度峰值a1時(shí)，為泡沫鋁材料進(jìn)入塑性變形平臺(tái)化階段的時(shí)候，a1的大小即為平臺(tái)應(yīng)力與泡沫鋁材料和航行體之間的接觸面積的乘積除以航行體的質(zhì)量，即a1=σp·A/m，其中A為泡沫鋁材料和航行體之間的接觸面積，m為航行體的質(zhì)量。泡沫鋁的μ越小，平臺(tái)應(yīng)力σp就越小,a1越小。a2則和航行體剩余動(dòng)能有關(guān)，剩余動(dòng)能越大，a2越大。

圖25 不同入水速度下最優(yōu)μFig.25 Optimal μ under different water entry speeds

泡沫鋁單位體積到達(dá)密實(shí)化階段時(shí)達(dá)到材料的最大吸能能力，此時(shí)吸收的最大能量由

(9)

得到。計(jì)算不同μ的泡沫鋁最大吸收的能量，得到圖22中最大吸能能力變化曲線。在μ<30%時(shí)，相對(duì)密度越大，泡沫鋁的吸能能力就越好，塑性變形吸收的能量就越多，航行體的剩余能量就越小，第二個(gè)加速度峰值a2也就越小，當(dāng)a1和a2達(dá)到平衡時(shí)，航行體加速度的最大值達(dá)到最低。

圖26 最大吸能能力變化曲線Fig.26 Variation curve of maximum energy absorption capacity

當(dāng)入水速度越小時(shí)，航行體的初始動(dòng)能也就較小，經(jīng)緩沖吸能后的剩余能量也就越小，此時(shí)第二個(gè)加速度峰值a2也越小，a1和a2的平衡點(diǎn)也就越低，所以入水速度越小，降載效果最好的泡沫鋁材料的μ也就越低。

隨著入水速度的增大，航行體初始動(dòng)能隨之增大，而緩沖材料的體積是一定的，雖然μ<30%時(shí)，隨著降載效果最好的泡沫鋁材料的相對(duì)密度的增加，其單位體積吸能能力有所增加，但其增幅遠(yuǎn)遠(yuǎn)低于航行體動(dòng)能的增加量，故隨著入水速度的增大，緩沖材料的緩沖效果降低。

3 結(jié) 論

本文采用了CEL方法，分析了泡沫鋁的力學(xué)性能，研究了帶緩沖材料航行體垂直入水時(shí)的動(dòng)力響應(yīng)問(wèn)題，通過(guò)對(duì)μ不同的泡沫鋁的緩沖效果進(jìn)行了對(duì)比，討論了對(duì)航行體在以一定速度入水時(shí)降載效果最好的泡沫鋁參數(shù)優(yōu)化問(wèn)題，找到了入水速度和最優(yōu)泡沫鋁參數(shù)之間的規(guī)律，同時(shí)驗(yàn)證了最優(yōu)參數(shù)泡沫鋁的降載效果，表明緩沖材料的性能直接影響到航行體緩沖裝置的降載效果。得出如下結(jié)論：

(1)帶緩沖材料的航行體入水時(shí)的受力過(guò)程存在兩個(gè)峰值，第一個(gè)峰值a1與σp的大小有關(guān)，隨著μ的增大而增大，第二個(gè)峰值a2與緩沖材料降載后的剩余動(dòng)能有關(guān)，隨著μ的增大而減小。當(dāng)所選泡沫鋁降載過(guò)程中，兩個(gè)峰值基本相等時(shí)，則此μ為最優(yōu)參數(shù)。

(2)泡沫鋁的μ不同時(shí)，吸能效果不同，其降載性能隨著μ的增加先升高再降低。當(dāng)μ較小時(shí)，航行體受力過(guò)程中，a1a2時(shí)，隨著μ的增加，整體降載效果反而降低。

(3)入水速度越大時(shí)，降載效果最好的泡沫鋁的μ越大。在緩沖材料的體積不能增加的情況下，入水速度越大，緩沖材料能起到的降載效果越低。

(4)緩沖材料能夠避免航行體觸水瞬間產(chǎn)生的瞬時(shí)巨大沖擊力，并能降低后續(xù)頭部入水階段由于入水阻力產(chǎn)生的加速度峰值。在入水速度分別為50 m/s，100 m/s和150 m/s時(shí)，加速度最大值可以分別降低90.8%，80.3%和69.7%，具有良好的緩沖降載效果。研究結(jié)果對(duì)今后航行體入水降載研究具有一定的借鑒作用，有助于指導(dǎo)降載結(jié)構(gòu)的設(shè)計(jì)。