陳開顏，陳 輝，魏海鵬，吉 康，裴金亮，程少華，王占瑩

（1.西北工業(yè)大學(xué)，西安 710072；2.北京宇航系統(tǒng)工程研究所，北京 100076）

0 引 言

氣囊裝置被廣泛應(yīng)用于重裝物資空投緩沖、航天器著陸回收、船舶下水和飛機(jī)水面迫降等領(lǐng)域[1-5]。回轉(zhuǎn)體安裝氣囊裝置，可實(shí)現(xiàn)回轉(zhuǎn)體的安全回收和重復(fù)使用，同時(shí)可避免回轉(zhuǎn)體的回落導(dǎo)致發(fā)射平臺(tái)受損[5]。國(guó)內(nèi)外專家對(duì)不同類型物體的落水和入水問題進(jìn)行了大量的研究，研究表明，數(shù)值仿真分析和落水試驗(yàn)是研究帶囊回轉(zhuǎn)體落水物理過程的有效手段。

在數(shù)值仿真分析方面，Shin和Chisum[6]采用CEL 方法對(duì)水下沖擊問題進(jìn)行了仿真分析；Hsu等[7]基于CEL 方法，采用三維氣泡動(dòng)力學(xué)模型對(duì)自由場(chǎng)氣泡載荷脈動(dòng)、射流過程等進(jìn)行了仿真分析；姚小虎等[8]基于CEL方法建立了水陸兩棲飛機(jī)機(jī)頭水上降落的流固耦合模型，分析了在不同入水角度下和垂直入水速度下的結(jié)構(gòu)響應(yīng)；王永虎等[9]基于CEL方法對(duì)失事客機(jī)（MH370）的水上迫降沖擊動(dòng)力學(xué)進(jìn)行了分析研究；王升等[10]基于CEL 方法對(duì)全雷入水沖擊問題進(jìn)行了建模與仿真分析，并結(jié)合入水試驗(yàn)，優(yōu)化了建模過程以及仿真方法。

在落水試驗(yàn)方面，Zhao 等[11]針對(duì)升角為30°的楔形體和典型船舷剖面結(jié)構(gòu)開展了入水沖擊試驗(yàn)；Panciroli等[12]針對(duì)彈性楔形體的板厚、底部升角、沖擊速度等開展了試驗(yàn)研究；何春濤[13]針對(duì)圓柱體開展了不同攻角、不同落水速度條件下的入水空泡試驗(yàn)研究；林賡[5]開展了低速小尺度帶囊圓柱體的原理試驗(yàn)研究，并基于CEL方法進(jìn)行了多因素的優(yōu)化仿真分析；Bardet[14]基于SAMCEF 有限元分析模型，設(shè)計(jì)了一種大尺度的落水試驗(yàn)系統(tǒng)，并結(jié)合落水試驗(yàn)，測(cè)試了一款大型氣囊回收系統(tǒng)（IRSS）的可行性。

本文基于CEL 方法和不同尺度落水試驗(yàn)對(duì)帶囊回轉(zhuǎn)體落水特征的影響因素、落水物理過程進(jìn)行辨識(shí)，并基于落水試驗(yàn)對(duì)不同囊-體連接結(jié)構(gòu)形式的落水特征進(jìn)行了分析，可為典型帶囊回轉(zhuǎn)體落水回收系統(tǒng)的設(shè)計(jì)提供參考。

1 CEL數(shù)值仿真方法

1.1 CEL算法

CEL 方法是Abaqus 軟件中分析流固耦合的關(guān)鍵技術(shù)，能夠解決有關(guān)大變形和材料破壞等諸多問題。利用Abaqus 軟件的動(dòng)力學(xué)分析模塊，建立落水載荷求解模型。CEL 方法是基于中心差分的動(dòng)態(tài)顯式求解算法，利用中心差分算法求解動(dòng)力學(xué)方程時(shí)，加速度和速度可以用位移的形式表示為[15]

1.2 氣-固耦合

氣囊是靠充入壓縮氣體起承載作用的，這是氣囊的基本特性。當(dāng)氣囊表面上的載荷發(fā)生變化時(shí)，氣囊腔內(nèi)的氣體壓力也隨之變化，從而實(shí)現(xiàn)氣囊反作用力和載荷的動(dòng)態(tài)平衡。

氣囊腔內(nèi)的氣體質(zhì)量一定，其容積表達(dá)式是腔內(nèi)充入氣體的壓力和溫度的函數(shù)，即=(p,θ,m)，其中p為氣體壓力，θ為氣體溫度，m為氣體質(zhì)量。根據(jù)虛功原理，我們可以得出

式中，δ∏*為氣囊所做的虛功，δ∏為沒有充入氣體時(shí)氣囊所作的虛功，V為實(shí)際氣體的容積，如果出現(xiàn)負(fù)值，說明氣腔內(nèi)容積增大，從而釋放能量，對(duì)外做功。對(duì)上式虛功等式進(jìn)行微分得

式中，pdδV為壓力載荷剛度，d/dp為氣體壓力容積。由于氣囊內(nèi)所有單元的壓力都相同，根據(jù)虛功原理，單個(gè)單元虛功可以用總虛功來表示：

根據(jù)充入氣囊內(nèi)氣體的總質(zhì)量等于氣體單元的總質(zhì)量，腔內(nèi)的氣體密度可以表示為

式中，θR和pR是參考密度ρR時(shí)的溫度和壓力，θA是在絕對(duì)零度時(shí)的溫度，pA是環(huán)境壓力。氣體的容積計(jì)算如下：

相應(yīng)的氣體壓力容積表達(dá)式為

其容積的變化可以表示為

2 落水試驗(yàn)系統(tǒng)及原理性試驗(yàn)

2.1 落水試驗(yàn)系統(tǒng)

本文建立的落水試驗(yàn)系統(tǒng)如圖1所示，包含回轉(zhuǎn)體、氣囊裝置、釋放裝置、地面支持設(shè)施（含水池、吊車等）和攝像系統(tǒng)等。

圖1 落水試驗(yàn)系統(tǒng)原理圖Fig.1 Sketch of falling experiment system

回轉(zhuǎn)體一般包含回轉(zhuǎn)結(jié)構(gòu)體、測(cè)量系統(tǒng)和充氣管路系統(tǒng)?；剞D(zhuǎn)結(jié)構(gòu)體是回轉(zhuǎn)體的結(jié)構(gòu)主體；測(cè)量系統(tǒng)用于測(cè)量落水過程中的參數(shù)（如氣囊囊內(nèi)壓力、氣囊連接帶拉力和慣性參數(shù)等）；充氣管路系統(tǒng)為氣囊裝置充氣并保證氣密。

氣囊裝置一般包含囊體、進(jìn)氣管、氣囊連接帶和囊-體連接結(jié)構(gòu)等。囊體為氣囊裝置的主結(jié)構(gòu)，用于儲(chǔ)存氣體，落水過程中為回轉(zhuǎn)體提供主要浮力；進(jìn)氣管為氣囊裝置的進(jìn)氣通道，充氣管路連接進(jìn)氣管為氣囊充氣；氣囊連接帶連接氣囊與連接結(jié)構(gòu)，落水過程中承受主要拉力載荷；囊體與回轉(zhuǎn)體之間通過連接結(jié)構(gòu)連接，典型囊-體連接結(jié)構(gòu)形式一般有點(diǎn)式（如圖2所示）和環(huán)式[14]（如圖3所示）。

圖2 典型點(diǎn)式連接結(jié)構(gòu)形式Fig.2 Typical point connection method

圖3 典型環(huán)式連接結(jié)構(gòu)形式Fig.3 Typical ring connection method

釋放裝置一般有電磁釋放和火工釋放兩種形式，本文采用電磁釋放裝置。地面支持設(shè)施包含水池或水槽、吊車和打撈回收設(shè)備等。攝像系統(tǒng)一般有常速攝像和高速攝像，用于記錄試驗(yàn)過程中的物理景象。

2.2 帶囊回轉(zhuǎn)體原理性落水試驗(yàn)研究

基于圖1 所示的落水試驗(yàn)系統(tǒng)原理圖，在實(shí)驗(yàn)室環(huán)境下利用1.2 m×1.2 m×2.4 m（長(zhǎng)×寬×高）的玻璃水箱開展直徑為50 mm、長(zhǎng)度為260 mm、重量約0.5 kg的帶囊回轉(zhuǎn)體落水原理性試驗(yàn)。并將高速攝像系統(tǒng)與水面保持水平，獲取了典型落水過程的物理景象。如圖4所示，典型帶囊回轉(zhuǎn)體落水過程可分為3個(gè)階段：

圖4 帶囊回轉(zhuǎn)體落水原理性試驗(yàn)物理景象Fig.4 Principle falling experiments of cylinder with airbags

（1）入水砰擊階段（A→B）：回轉(zhuǎn)體底部觸水至氣囊觸水。本階段主要是回轉(zhuǎn)體受到入水沖擊力作用，并在回轉(zhuǎn)體底部形成開口空泡。

（2）氣囊與回轉(zhuǎn)體瞬態(tài)作用階段（B→D）：氣囊接觸水面至氣囊相對(duì)回轉(zhuǎn)體靜止（完全浸入水中）。本階段回轉(zhuǎn)體繼續(xù)向下運(yùn)動(dòng)，回轉(zhuǎn)體表面形成空泡，逐漸閉合，氣囊與回轉(zhuǎn)體之間相互作用，并且氣囊繞回轉(zhuǎn)體壁面翻轉(zhuǎn)碰撞，受力特征較為復(fù)雜。

（3）氣囊與回轉(zhuǎn)體整體運(yùn)動(dòng)階段（D→E）：氣囊與回轉(zhuǎn)體一起運(yùn)動(dòng)。本階段回轉(zhuǎn)體表面空泡逐漸潰滅，氣囊與回轉(zhuǎn)體共同運(yùn)動(dòng)至最大入水深度，而后回轉(zhuǎn)體在氣囊浮力、自身浮力等綜合作用力作用下上浮。

通過小型帶囊回轉(zhuǎn)體的原理性試驗(yàn)辨識(shí)出帶囊回轉(zhuǎn)體的典型落水過程及物理景象，同時(shí)證明落水試驗(yàn)系統(tǒng)具有可行性。

3 不同影響因素下的落水特征仿真分析研究

對(duì)縮比回轉(zhuǎn)體模型在不同工況（不同落水速度、不同俯仰角度以及不同初始囊壓）下基于CEL 算法進(jìn)行數(shù)值仿真模擬。通過對(duì)計(jì)算結(jié)果進(jìn)行分析比較得到落水過程中落水深度、囊壓變化和氣囊總拉力受不同影響因素的變化規(guī)律。

3.1 物理量定義及初始邊界條件

文中對(duì)描述物理量均進(jìn)行了無量綱處理，其中：落水深度用h/L表示，h為回轉(zhuǎn)體浸入水面深度，L為回轉(zhuǎn)體軸向長(zhǎng)度；囊壓用P/P0表示，P為囊壓絕對(duì)值，P0為標(biāo)準(zhǔn)大氣壓；氣囊連接帶拉力用F/B表示，F(xiàn)為氣囊連接帶合力，B為回轉(zhuǎn)體全浸沒時(shí)的自身浮力；落水速度用V/gL表示，V為落水速度，g為重力加速度。俯仰角用兩氣囊連接點(diǎn)與回轉(zhuǎn)體軸線形成的平面與水平面之間的夾角表示，圖4 所示為90°俯仰角的入水姿態(tài)。

初始邊界條件：將帶囊回轉(zhuǎn)體整體設(shè)置于一個(gè)無限范圍的水域，整體流場(chǎng)模型寬度設(shè)置為回轉(zhuǎn)體直徑的10 倍，流場(chǎng)四周以及上下底面設(shè)置為剛性固定邊界。由于建立的流域足夠大，可認(rèn)為四周壁面反射的微小水波不會(huì)對(duì)整個(gè)落水過程產(chǎn)生影響。氣囊與回轉(zhuǎn)體位于自由液面的正上方。根據(jù)運(yùn)動(dòng)學(xué)原理V2= 2gh，賦予回轉(zhuǎn)體與氣囊初始下落速度以代替下落高度的效果，使回轉(zhuǎn)體與氣囊處于某一速度時(shí)，回轉(zhuǎn)體下表面恰好與水面相切。氣囊與回轉(zhuǎn)體在不受其他約束的條件下開始落水運(yùn)動(dòng)。

3.2 落水速度的影響分析

分別開展了俯仰角為90°、初始囊壓為P/P0=1.3、落水速度分別為V/gL=2.09和V/gL=1.74的數(shù)值仿真計(jì)算，結(jié)果見圖5～7。隨著落水速度的增大，直接影響帶囊回轉(zhuǎn)體落水時(shí)的初始動(dòng)能，繼而影響落水深度增加，囊壓峰值增大，氣囊總拉力增大。結(jié)合原理性落水試驗(yàn)的物理景象對(duì)仿真結(jié)果進(jìn)行分析，圖5～6標(biāo)示了典型落水物理特征點(diǎn)。

圖5 不同速度下落水深度仿真分析結(jié)果Fig.5 Simulation results of falling depth under different speeds

圖6 不同速度下囊壓變化仿真分析結(jié)果Fig.6 Simulation results of airbag pressure under different speeds

圖7 不同速度下氣囊總拉力仿真分析結(jié)果Fig.7 Simulation results of total pull on airbag under different speeds

3.3 俯仰角度的影響分析

分別開展了初始囊壓為P/P0=1.3、落水速度為V/gL=2.09、俯仰角分別為0°和30°的數(shù)值仿真計(jì)算，結(jié)果見圖8～10。隨著俯仰角的增大，帶囊回轉(zhuǎn)體的落水深度和氣囊總拉力變大。受氣囊接觸水時(shí)刻的速度影響，不同俯仰角下的最大囊壓不同，俯仰角為0°時(shí)的囊壓增值最大。

圖8 不同俯仰角度下落水深度仿真分析結(jié)果Fig.8 Simulation results of falling depth under different pitch angles

圖9 不同俯仰角度下囊壓變化仿真分析結(jié)果Fig.9 Simulation results of airbag pressure under different pitch angles

圖10 不同俯仰角度下氣囊總拉力仿真分析結(jié)果Fig.10 Simulation results of total pull on airbag under different pitch angles

3.4 初始囊壓的影響分析

分別開展了俯仰角為0°、落水速度為V/gL=2.09、初始囊壓分別為P/P0=1.1和P/P0=1.5的數(shù)值仿真計(jì)算，結(jié)果見圖11～13。隨著初始囊壓的增大，影響帶囊回轉(zhuǎn)體氣囊裝置入水后的囊壓峰值（絕對(duì)壓力）和落水后的囊壓峰值（絕對(duì)壓力）就越大。對(duì)落水深度和氣囊總拉力影響不明顯，僅高囊壓狀態(tài)入水深度稍小，氣囊總拉力差異較小。

圖11 不同初始囊壓下落水深度仿真分析結(jié)果Fig.11 Simulation results of falling depth under different initial pressures

圖12 不同初始囊壓下囊壓變化仿真分析結(jié)果Fig.12 Simulation results of airbag pressure under different initial pressures

圖13 不同初始囊壓下氣囊總拉力仿真分析結(jié)果Fig.13 Simulation results of total pull on airbag under different initial pressures

4 大尺度帶囊回轉(zhuǎn)體落水試驗(yàn)研究

基于上述仿真分析結(jié)果，對(duì)點(diǎn)式和環(huán)式連接結(jié)構(gòu)形式下的帶囊回轉(zhuǎn)體分別開展了兩次大尺度落水試驗(yàn)，對(duì)實(shí)際工程應(yīng)用的帶囊回轉(zhuǎn)體落水特征進(jìn)行分析，用于指導(dǎo)典型帶囊回轉(zhuǎn)體落水回收系統(tǒng)的設(shè)計(jì)及應(yīng)用。

4.1 落水深度和速度分析

測(cè)量系統(tǒng)測(cè)得各發(fā)次的落水深度和落水速度，如圖14所示。可見不同連接結(jié)構(gòu)形式下的帶囊回轉(zhuǎn)體落水深度和落水速度差別不大：在V/gL=1.45的相同落水速度下，點(diǎn)式最大落水深度h/L均值為1.3，環(huán)式最大落水深度h/L均值為1.1。

圖14 帶囊回轉(zhuǎn)體落水深度和落水速度Fig.14 Falling speed and depth of cylinder with airbags

4.2 囊壓分析

不同連接結(jié)構(gòu)形式下帶囊回轉(zhuǎn)體落水后的氣囊裝置囊壓變化差別較小，且趨勢(shì)相近，如圖15 所示。分別有兩個(gè)峰值，結(jié)合原理性落水試驗(yàn)物理過程分析可知：第一個(gè)峰值為氣囊撞水至完全浸沒后，囊壓升至第一個(gè)峰值；第二個(gè)峰值為氣囊與回轉(zhuǎn)體碰撞擠壓后，囊壓升至第二個(gè)峰值。

圖15 帶囊回轉(zhuǎn)體囊壓變化Fig.15 Airbag pressure of cylinder with airbags

4.3 連接帶受力分析

如圖16～18 所示，環(huán)式連接結(jié)構(gòu)與點(diǎn)式連接結(jié)構(gòu)相比，連接帶總拉力相差較小，但環(huán)式結(jié)構(gòu)連接帶拉力分布更均勻，單根帶拉力的幅值更小，主要原因是點(diǎn)式連接結(jié)構(gòu)存在附加彎矩，測(cè)力螺栓受力存在放大現(xiàn)象。因此，環(huán)式連接結(jié)構(gòu)能夠有效實(shí)現(xiàn)氣囊連接帶降載。

圖16 點(diǎn)式連接結(jié)構(gòu)各連接帶受力峰值情況Fig.16 Peak load on every belt with point connection method

圖17 環(huán)式連接結(jié)構(gòu)各連接帶受力峰值情況Fig.17 Peak load on every belt with ring connection method

圖18 氣囊連接帶合力圖Fig.18 Total load on connection belts

5 結(jié) 論

本文針對(duì)帶囊回轉(zhuǎn)體落水回收問題，通過建立帶囊回轉(zhuǎn)體落水回收試驗(yàn)系統(tǒng)，并基于落水仿真分析和不同尺度的帶囊回收體落水試驗(yàn)，得到以下結(jié)論：

（1）帶囊回轉(zhuǎn)體落水過程可分為三個(gè)階段：入水砰擊階段、氣囊與回轉(zhuǎn)體瞬態(tài)作用階段和氣囊與回轉(zhuǎn)體整體運(yùn)動(dòng)階段。

（2）基于CEL 數(shù)值仿真方法，建立了帶囊回轉(zhuǎn)體落水仿真分析模型，經(jīng)不同工況的仿真分析，辨識(shí)出落水過程的影響因素主要包括：落水速度、俯仰角度和初始囊壓。

（3）通過開展不同的囊-體連接形式下的多發(fā)次大尺度帶囊回轉(zhuǎn)體落水試驗(yàn)，有效驗(yàn)證了落水試驗(yàn)系統(tǒng)的可行性。

（4）在氣囊、連接結(jié)構(gòu)和回轉(zhuǎn)體組成的力學(xué)系統(tǒng)中，它們之間的作用力僅是系統(tǒng)的內(nèi)力，不影響在相同外力作用下的力學(xué)系統(tǒng)的宏觀運(yùn)動(dòng)，不同連接結(jié)構(gòu)形式對(duì)帶囊回轉(zhuǎn)體落水深度和囊壓變化影響較??；環(huán)式連接結(jié)構(gòu)的氣囊連接帶較點(diǎn)式連接結(jié)構(gòu)受力均勻，能夠?qū)崿F(xiàn)氣囊連接帶有效降載。