国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

高速破片穿透液艙的數(shù)值模擬研究

2019-04-03 01:17徐思博孟子飛劉文韜曹雪雁
振動與沖擊 2019年6期
關鍵詞:液艙彈片壁板

徐思博, 孟子飛, 劉文韜, 曹雪雁

(哈爾濱工程大學 船舶工程學院,哈爾濱 150001)

高速破片作為武器爆炸攻擊艦船的重要利器之一,其產(chǎn)生的高壓環(huán)境難以測量,初始速度通常在1 000 m/s左右,其形狀的隨機性以及強穿透性,使得如何防護破片毀傷成為研究的難點。如今國際公認的防護手段是在艙壁外加設防護液艙用來衰減破片的沖擊,研究水下近場爆炸高速破片在液艙中的衰減特性有極其重要的軍事戰(zhàn)略意義[1]。針對高速破片衰減特性,目前可行的方法主要有基于理想模型的數(shù)值驗證和一些實驗驗證。高速破片擊穿船舶液艙的整個過程大致分為沖擊、拉伸、空穴形成和潰滅以及穿透內(nèi)板四個階段[2]。沖擊階段,破片接觸液艙前壁并開始穿透,在艙壁后的液體內(nèi)產(chǎn)生擾動;拉伸階段,高速破片撞擊液體并產(chǎn)生沖擊波,之后在液體中前進;空穴的形成和潰滅階段,高速破片在液體中的速度被迅速衰減,形成空穴并演化直至潰滅;穿透內(nèi)板階段,高速破片抵達液艙內(nèi)壁并穿出。

穿透過程中有幾個關鍵性的難點,包括:材料的動態(tài)力學性能和力學響應,高速破片水中速度衰減問題,水中沖擊波形成的機理,空穴的形成和沖擊波的疊加等。20世紀初期, Worthington等[3]用高速攝像機記錄下圓球落水的動態(tài)過程中的飛濺和空泡現(xiàn)象,并以此為基礎開啟水彈道學的發(fā)展。我國在這項領域的研究開始于20世紀末期,僅做過少量的入水沖擊、入水空泡實驗[4-6]。陳先富[7]進行了不同形狀彈丸入水的試驗研究,記錄了空穴的產(chǎn)生潰滅以及脈動,但由于試驗工況較少,僅得出初始速度越高空穴增大越快這一結論。高速侵徹實驗危險性大,實驗條件要求高,同時隨著計算機技術的飛速發(fā)展,對侵徹問題應用數(shù)值仿真方式進行模擬研究成為了解決該問題的有效且經(jīng)濟的途徑之一。在21世紀,徐雙喜等[8]依照破片穿甲運動方程以及德·瑪爾模型,分析得出圓柱形破片穿透背水板的剩余速度公式,并采用Autodyn中的數(shù)值耦合算法模擬該過程,孔祥韶等[9]同時對雙發(fā)破片的相互耦合作用模式進行分析,其主要研究為穿透后的剩余速度,并未涉及空穴演化等現(xiàn)象。Varas等[10]開展了低速以及高速下的球型破片擊穿液艙的實驗,并用高速成像儀記錄下了高速破片形成的沖擊波和空穴,并運用ALE算法進行了數(shù)值仿真,但彈片衰減速度、水中壓力變化等仍有較大偏差。張婧等[11-13]則關注到多層防護結構及復合材料裝甲加持下破片對結構的毀傷問題,與此同時沈曉樂等[14]在試驗中對方形破片阻力系數(shù)進行修正,并通過試驗研究破片的侵徹能力以及剩余特性。近幾年,國內(nèi)針對破片形狀對侵徹的影響做了很多工作[15-17],但并未引進新的仿真方法。楊文山等[18-19]則應用SPH方法對侵徹過程進行仿真,這種無網(wǎng)格光滑粒子在處理破口界面處有極好的模擬效果。

本文根據(jù)Varas等開展的球型彈片擊穿液艙的實驗數(shù)據(jù),應用有限元軟件ABAQUS中的CEL算法對該過程進行數(shù)值仿真。試驗針對不同初始速度和裝水量下的壓力測點、空穴徑寬、彈片速度衰減與實驗數(shù)據(jù)進行對比,驗證了CEL在模擬高速侵徹過程的適用性問題。

1 CEL算法的基本方程

CEL方法同時具有拉格朗日方法和歐拉方法的優(yōu)勢,較好的處理大變形、碰撞及一些流固耦合問題[20]。CEL在處理耦合問題中,應用拉格朗日算法計算結構以保持其在響應上的精度,用歐拉算法計算流體介質(zhì)以解決其流動過程中的大變形問題。拉格朗日算法是一種離散化數(shù)值求解方法,其求解過程中,材料會依附在單元網(wǎng)格上,在網(wǎng)格的運動和變形的過程中伴隨著材料的流動,計算過程中,單元會發(fā)生變形但單元的初始質(zhì)量不變。具體控制方程為

(1)

(2)

(3)

式中:ρ,u,e,σ分別為拉格朗日流體介質(zhì)的密度、速度、比內(nèi)能和應力。

歐拉算法是一種固定網(wǎng)格而使材料在網(wǎng)格間運輸?shù)那蠼夥椒ǎ嬎氵^程中一個歐拉單元按用戶定義的材料分配,有著不同的物質(zhì),整個計算過程中空間網(wǎng)格將會固定,而材料網(wǎng)格將以拉格朗日單元的狀態(tài)變量實時傳輸?shù)娇臻g網(wǎng)格中,具體控制方程為

(4)

(5)

(6)

式中:ρ,u,e,σ,f分別為歐拉流體介質(zhì)的密度、速度、比內(nèi)能、應力和質(zhì)量力。

在耦合算法之前,這兩種算法是相互獨立的,同時在耦合計算時通過耦合面的控制,歐拉的材料流動也不會穿透拉格朗日單元。耦合作用則通過在界面交換壓力載荷和幾何約束來實現(xiàn),如圖1所示。

圖1 耦合界面處理示意圖Fig.1 Coupling interface processing

2 高速破片穿透液艙的數(shù)值分析

高速破片穿透液艙的數(shù)值模型使用有限元仿真進行建立運算分析,符合工程應用背景以及科學研究需求。

2.1 液艙結構

鑒于有限元分析精準度與網(wǎng)格尺寸密切相關,為提高運算精度以及計算效率,仿真分析將采用實驗整體模型的一半。根據(jù)Varas等實驗中液艙的尺寸結構,液艙用6063-T5鋁制成,長750 mm,寬150 mm,壁厚2.5 mm,液艙兩側用厚30 mm的PMMA有機玻璃封閉,在液艙內(nèi)部設有兩個壓力測點PTn與PTf,液艙外部前后壁板設有應力及撓度測點G1與G2,測點位置如圖4所示。破片為球形,直徑12.5 mm,質(zhì)量8 g,實驗中選取600 m/s和900 m/s了兩種初始速度。本次仿真針對該模型在壁厚方向(即彈片前進方向)設置3層網(wǎng)格單元,接近彈片侵徹位置附近皆為正六面體網(wǎng)格,網(wǎng)格尺寸為0.83 mm,最終液艙模型網(wǎng)格單元數(shù)為49 540個,如圖2和圖3所示。

仿真試驗采用J-C模型描述液艙的本構關系。J-C模型的關系式為

(7)

(8)

式中:θ為當前溫度;θmelt為熔化溫度;θtransition為轉化溫度。

在J-C失效模型中,斷裂應變定義式為

(9)

表1 結構材料表

圖2 液艙和破片數(shù)值模型Fig.2 Numerical modelling of fuel tanks and fragment

圖3 破片入射面網(wǎng)格圖Fig.3 Detail of the entry wall mesh

圖4 壓力測點位置示意圖Fig.4 Sketch of pressure transducers

2.2 歐拉域

鑒于高速沖擊問題帶來的大變形以及流體的不規(guī)則流動,仿真試驗將用歐拉網(wǎng)格建立水和空氣的模型。在彈片穿透過程中,結構的大變形同時會伴隨著流體的運動,因此歐拉域的設置應大于結構的區(qū)域以保證流體在流動過程中不會消失,保證能量守恒,所以歐拉域設置為190 mm寬。水和空氣的材料屬性通過體積分數(shù)方式離散進歐拉域內(nèi),最終歐拉網(wǎng)格單元數(shù)為4 061 300個,如圖5所示。

圖5 歐拉網(wǎng)格展示圖Fig.5 Mesh of the fluids

水的狀態(tài)方程采用Mie-Grüneisen狀態(tài)方程,定義式為

(10)

式中:ρ0為流體密度;S,c0,Γ0為流體常數(shù);η=1-ρ0/ρ??諝獠捎美硐霘怏w模型。水和空氣的具體材料數(shù)據(jù)見表2。

表2 水和空氣的材料參數(shù)表

3 數(shù)值計算結果分析

近場水下爆炸產(chǎn)生的破片速度通常在1 000 m/s左右,因此在本節(jié)中,將給出裝水量75%和60%,彈片初始速度從600~1 200 m/s速度梯度下的工況數(shù)值仿真結果并與Varas等的實驗數(shù)據(jù)進行對比,并分析裝水量75%時不同入射速度下前后壁板應力及撓度對比。

3.1 破片高速侵徹過程

由仿真的結果可知,破片侵徹過程分為四個階段,即沖擊、拉伸、空穴形成和潰滅以及穿透內(nèi)板,每個階段在數(shù)值仿真中得到了充足的展現(xiàn)。在彈片穿透的瞬間,結構的應力峰值達到257 MPa,并沿沖擊點向四周擴散圓形波,如圖6(a)所示。在沖擊階段,當彈片穿透壁板之后,彈片的動能迅速在流體內(nèi)形成高壓的環(huán)境,并向四周擴散沖擊波,如圖6(b)所示。在拉伸階段和空穴形成潰滅的階段,彈片穿透水箱在流體內(nèi)部留下一片空穴區(qū),流體在慣性的作用下不斷擴散使得空穴不斷變化,而空穴潰滅所需時間遠遠大于高速彈片穿透水箱所需時間,如圖6(c)所示。在穿透階段,由于流體內(nèi)沖擊波預先到達壁板,彈片在到達之前就在壁板處產(chǎn)生預應力,圖6(d)給出穿透階段液艙的應力云圖。

圖6 初始速度900 m/s裝水量75%的高速穿透過程Fig.6 HRAM phases in a tube impacted at 900 m/s and filled 75%

3.2 不同入射速度和裝載量對破片速度衰減的影響

圖7給出四組工況下速度衰減曲線與理論分析解和Varas等實驗值得對比情況。可以看出數(shù)值分析方法CEL的結果與理論分析結果非常接近,只是CEL算法在彈片穿透液艙時會有較大的減速。由圖7看出,CEL算法在彈片速度衰減趨勢上與理論分析和實驗有較好的擬合;由于彈片在液艙內(nèi)滯留時間極短,液體阻滯作用有限,同時液面距彈道位置較遠,液面效應較小,液艙內(nèi)裝水量對速度衰減趨勢并沒有影響;彈片在穿透入射面時速度大致衰減15%左右。

圖7 速度衰減時歷曲線對比(實驗值與理論值源自Varas等的研究)Fig.7 Comparison of velocity decay vs. time (data of experiment and analysis come from Varas et al)

3.3 不同液體裝載量對彈道空穴演化的影響

圖8和圖9給出由Varas等的實驗數(shù)據(jù)中特定時刻空穴大小與CEL同一時刻空穴大小對比??梢钥闯鯟EL仿真結果與實驗數(shù)據(jù)空穴大小、水面隆起的波高等十分接近。同時也發(fā)現(xiàn),入射面與射出面的箱體呈現(xiàn)向外隆起的變形,由于穿透瞬間與彈片接觸位置應力高達257 MPa,高于壁板的屈服應力使其瞬間形成穿透圓孔,而當彈片穿透之后,在內(nèi)部流體的作用下,壁板所受應力維持在100~200 MPa,使液艙壁板進入塑性變形階段發(fā)生屈曲卻不至于形成斷裂破損,而此時液體量的多少則是影響壁板屈曲變形的主要因素,如圖8和圖9所示,圖中框線為液面初始位置。

液體裝載量的變化決定了彈道距離自由液面的相對位置,并影響及導致彈道出現(xiàn)明顯的非對稱空穴現(xiàn)象,液體裝載量從75%減小到60%,即彈道與自由液面的距離從6倍彈片半徑減小到2.4倍彈片半徑,非對稱差異可由20%增加到48%,這種非對稱差異導致空穴內(nèi)壓不平衡從而對彈道穩(wěn)定性的影響很大。

圖8 初始速度600 m/s裝水量75%下空穴演化0.26 ms,0.46 ms,0.60 ms時刻實驗數(shù)據(jù)與CEL算法數(shù)據(jù)的對比Fig.8 Cavity evolution at 0.26 ms, 0.46 ms, 0.60 ms obtained from experiments and CEL simulation for a tube impacted at 600 m/s, filled 75%

圖9 初始速度600 m/s裝水量60%下空穴演化0.13 ms,0.24 ms,0.38 ms時刻實驗數(shù)據(jù)與CEL算法數(shù)據(jù)的對比Fig.9 Cavity evolution at 0.13 ms, 0.24 ms, 0.38 ms obtained from experiments and CEL simulation for a tube impacted at 600 m/s, filled 60%

3.4 不同入射速度對液艙隔板壓力和變形的影響

根據(jù)實驗中兩個壓力測點PTn,PTf與數(shù)值仿真結果進行對比,對比壓力時歷曲線如圖10和圖11所示。由于實驗與數(shù)值仿真中彈片速度存在偏差,因此反映到壓力時歷曲線中對應出現(xiàn)峰值的時刻不盡相同,同時測點壓力的峰值有所偏差,尤其是PTf,但兩次壓力峰值的脈寬十分接近,結果基本吻合;在Ptn測點,幾次仿真結果均出現(xiàn)沖擊波的二次加載,其峰值大致為第一次波峰的15%,兩次波峰間脈寬小于0.05 ms,根據(jù)Ptn測點空間位置,推測二次加載為液艙底面反射稀疏波疊加形成(若是由射出面反射形成則脈寬應在0.12 ms左右,假定沖擊波為球面波)。

圖10 裝水量75%初始速度900 m/s的實驗值與數(shù)值的壓力時歷曲線對比Fig.10 Experimental and numerical pressure time historyPtn and Ptf in a tube 75% filled impacted at 900 m/s

圖11 裝水量75%初始速度600 m/s的實驗值與數(shù)值的壓力時歷曲線對比Fig.11 Experimental and numerical pressure time historyPtn and Ptf in a tube 75% filled impacted at 600 m/s

結構在達到塑性變形之前,分別經(jīng)歷線彈性階段以及非線性彈性階段,前置壁板在不同速度下應力響應如圖12所示,不同入射速度下結構的應力響應僅因穿透時間不同而有一定的時域上的滯后性,而結構的應力響應只與材料本構方程有關,前置壁板在破片穿透過程中達到140 MPa左右應力,并在反射稀疏波抵達前置壁板前有一定的應力卸載,最終在稀疏波與水波的共同作用下達到動態(tài)屈服強度時應力達到最大值,此現(xiàn)象在后置壁板的應力響應中有更為明顯的表現(xiàn),如圖13所示。圖14給出在前后壁板測點位置最終撓度,隨著初始速度的增加,前后隔板的變形撓度呈現(xiàn)出不同特征,即前置隔板的變形近乎呈線性關系,但幅度平緩,而后置隔板的變形呈非線性增長趨勢。

圖12 裝水量75%不同速度下G1應力曲線Fig.12 Pressure time history G1 for a tube impacted at different velocity, filled 75%

圖13 裝水量75%不同速度下G2應力曲線Fig.13 Pressure time history G2 for a tube impacted at different velocity, filled 75%

圖14 裝水量75%壁板測點最終撓度Fig.14 Deflection in the end on the entry and exit wall of a tube impacted at different velocity, filled 75%

4 結 論

本文基于CEL算法研究了破片的高速侵徹問題,數(shù)值結果呈現(xiàn)出明顯的四個階段,即沖擊、拉伸、空穴形成和潰滅以及穿透內(nèi)板。文中首先以Varas等開展的高速彈片穿透鋁制液艙實驗為基礎模型,驗證了數(shù)值模擬結果的可靠性,進而模擬了不同破片速度、不同液艙裝載量時高速侵徹過程的空穴演化特征,得出了以下結論:

(1) 在當前液艙參數(shù)下,彈片在穿透前置隔板時速度大致衰減15%左右;當彈道遠離自由液面時,液艙裝載量對破片速度衰減規(guī)律影響較小。

(2) 近自由液面入射彈道會出現(xiàn)明顯的非對稱空穴演化特征,隨著與自由液面距離變小,非對稱性差異可由20%增加到48%,這對彈道的穩(wěn)定性影響很大。

(3) 隨著破片初始速度增加,前后隔板的變形撓度呈現(xiàn)出不同特征,在當前速度范圍內(nèi),前置隔板的變形隨著破片入射速度增加而略有增加,但幅度平緩,而后置隔板的變形隨入射速度增加呈非線性增長趨勢。

猜你喜歡
液艙彈片壁板
雙浮板液艙晃蕩特性的數(shù)值研究
模具內(nèi)組裝的接觸元件級進模設計
基于CFD的大型船舶液艙晃蕩研究
GDX1 包裝機四輪彈片清潔裝置
大將
某大型飛機復合材料壁板工藝仿真及驗證技術
液艙建模方法對深水半潛式平臺總體強度評估的影響*
矩形敞口水池壁板轉角處的精細化設計
淺析非洲的八音盒
機翼下壁板裂紋擴展分析