張?jiān)浦?郭則慶,孫 帥

(南京理工大學(xué) 瞬態(tài)物理國家重點(diǎn)實(shí)驗(yàn)室， 南京 210094)

1 引言

射彈高速入水是一個(gè)包含空氣、水蒸氣、液態(tài)水等多相發(fā)展的復(fù)雜過程。當(dāng)射彈從空中高速進(jìn)入水中時(shí)，水被射彈排開形成空泡，空氣隨之填充空泡，與此同時(shí)空泡內(nèi)部表面由于壓力下降液體發(fā)生空化，形成了由空氣和水蒸氣混合的氣泡。射彈繼續(xù)運(yùn)動(dòng)后空泡在液面附近發(fā)生閉合，射彈周圍變?yōu)橥耆伤魵饨M成的超空泡。高速入水過程中空泡結(jié)構(gòu)發(fā)生快速變化，射彈結(jié)構(gòu)與空泡表面發(fā)生劇烈相互作用，射彈的入水穩(wěn)定性和阻力特性受到射彈結(jié)構(gòu)、入水角、攻角、初始擾動(dòng)等多個(gè)因素的影響，發(fā)展過程非常復(fù)雜。由于其復(fù)雜性和鮮明的軍事應(yīng)用價(jià)值，射彈高速入水研究受到了國內(nèi)外學(xué)者的廣泛關(guān)注。

在試驗(yàn)研究方面，Lee等[1]基于大量的試驗(yàn)數(shù)據(jù)、數(shù)值模擬結(jié)果和表面閉合理論，研究了剛性格球體高速垂直入水的空泡特性，提出了一個(gè)球體入水空泡閉合理論模型。Kulkarni[2]研究了超空泡彈丸的物理特性，建立了射彈與空泡壁接觸和射彈不與空泡壁接觸2個(gè)階段運(yùn)動(dòng)的數(shù)學(xué)模型。郭子濤[3]從實(shí)驗(yàn)和理論三方面研究不同彈頭的高速入水，建立彈丸入水空泡特性的解析模型。Song[4]對圓柱射彈以35°～90°、80～127.5 m/s進(jìn)行了試驗(yàn)，發(fā)現(xiàn)射彈的入水過程中舊空泡的塌陷與新空泡擴(kuò)張交替發(fā)生。李宜果[5]開展了圓柱射彈以45°～90°、150 m/s左右速度斜入水試驗(yàn)，提取極值位移與姿態(tài)角和空化數(shù)的關(guān)系并進(jìn)行曲面擬合，獲得噴濺環(huán)與空泡壁面極值位移的激勵(lì)關(guān)系。Chen[6]開展了射彈以20°～45°、150 m/s左右速度入水試驗(yàn)，研究發(fā)現(xiàn)射彈阻力系數(shù)主要受彈頭形狀影響。Truscott[7-8]主要利用高速攝像技術(shù)對不同形狀射彈以5°～15°入水角、380 m/s速度入水過程的跳彈現(xiàn)象進(jìn)行了研究，結(jié)果表明長徑比、空化器直徑等參數(shù)對入水彈道穩(wěn)定性有顯著影響。侯宇[9]在水靶道中試驗(yàn)研究了尾翼式超空泡射彈以8°入水角、650 m/s左右速度入水過程，得到了入水沖擊瞬間彈體姿態(tài)對噴濺形成、彈道軌跡變化和水下壓力波傳播的影響。

近年來，得益于計(jì)算機(jī)技術(shù)的進(jìn)步，射彈入水過程的數(shù)值研究取得了快速的發(fā)展。陳晨[10]考慮氣、汽、液三相流體的可壓縮性，對射彈以400 m/s垂直入水進(jìn)行數(shù)值模擬，分析了流體可壓縮性對激波、空泡演化、流場特性的影響。郝博[11]等采用數(shù)值方法研究了尾翼楔角對六尾翼槍彈的垂直入水阻力和運(yùn)動(dòng)穩(wěn)定性影響。郝亮[12]對圓柱形射彈在入水攻角-3°～3°的以1 500 m/s、45°入水角度的斜入水過程進(jìn)行仿真研究，分析得出射彈跨聲速入水過程中，液體可壓縮性對射彈受力特性和流場結(jié)構(gòu)有較大影響。王曉輝等[13-14]對尾翼式超空泡射彈帶初始攻角條件下的以1 000～1 500 m/s速度、60°入水角斜入水過程開展了數(shù)值研究，得到了攻角對尾翼射彈入水空泡特性和彈道穩(wěn)定性的影響。Akbari等[15-16]在Truscott試驗(yàn)研究[7-8]的基礎(chǔ)上，進(jìn)一步研究了帶階梯的圓柱錐射彈在10°～60°入水角范圍內(nèi)的入水穩(wěn)定性影響，結(jié)果表明階梯可提高射彈小角度入水的穩(wěn)定性。周可[17]對尖頭射彈以150 m/s、8°～25°斜入水進(jìn)行仿真研究，分析得到形成跳彈主要原因是射彈頭尾兩側(cè)受力不均。祁曉斌[18]對尾翼式超空泡射彈以825 m/s、5°～10°入水角度開展了仿真計(jì)算，結(jié)果發(fā)現(xiàn)入水初期的空泡不對稱性會(huì)導(dǎo)致彈體抬頭形成正攻角，引起跳彈和失穩(wěn)。

在跨介質(zhì)超空泡射彈設(shè)計(jì)時(shí)，為了兼顧空中和水下穩(wěn)定性，一般采用尾翼式結(jié)構(gòu)。在入水過程中尾翼不可避免會(huì)與水面發(fā)射碰撞沖擊，尾翼直徑對入水穩(wěn)定性和阻力特性都有較大的影響。因此，文中將研究不同直徑尾翼的射彈在10°～30°范圍內(nèi)的入水穩(wěn)定性，為射彈尾翼直徑的設(shè)計(jì)提供參考。

2 數(shù)值模擬方法

2.1 流體控制方程

流體體積(VOF)多相模型能夠解析混合相界面上的網(wǎng)格模擬幾種互不相溶流體的流動(dòng)。在高速小角度入水過程中，存在著液體、氣體、水蒸氣三相，VOF多相模型能準(zhǔn)確捕捉空泡的形狀和特征。流體的控制方程如下：

連續(xù)性方程，所有相的質(zhì)量守恒方程表示為：

(1)

式中：ui為流體微團(tuán)在x、y、z軸上的速度分量；ρm為混合相中各相的密度，其表達(dá)式為：

ρm=αlρl+αgρg+αvρv

(2)

式中：αl、αg、αv分別為水相、空氣相、水蒸氣相的體積分?jǐn)?shù)；ρl、ρg、ρv分別為水相、空氣相、水蒸氣相的密度。

動(dòng)量守恒方程為：

(3)

式中：uj為流體微團(tuán)在x，y，z軸上的速度分量；P為流體壓力；Fi為流體微團(tuán)在x，y，z軸上的受力；μm表示混合相動(dòng)力粘性系數(shù)。

考慮流體的可壓縮性，引入能量方程：

(4)

Tait方程用于忽略溫度影響可壓縮流體的狀態(tài)方程：

(5)

K=K0+nΔp

(6)

Δp=p-p0

(7)

式中：p0為參考壓力101 325 Pa；ρ0為參考壓力下的液體密度，取1 000 kg/m3；K0為參考壓力下的體積彈性模量，取2.2 GPa；n為密度指數(shù)，取7.15。

將空氣相和水蒸氣相視為理想氣體，狀態(tài)方程為：

p=ρgRT

(8)

式中：R為氣體常數(shù)；ρg為氣體密度。

2.2 空化模型

采用的空化模型為Schnerr-Sauer模型，由Rayleigh-Plesset空化模型簡化，忽略了空泡的加速度、黏性效應(yīng)和表面張力效應(yīng)的影響。其控制方程如下：

(9)

(10)

(11)

2.3 射彈運(yùn)動(dòng)方程

研究三維空間中射彈的入水過程，需要準(zhǔn)確定義射彈運(yùn)動(dòng)以及作用在射彈上的力和力矩。射彈的運(yùn)動(dòng)方程描述如下[15]：

(12)

(13)

(14)

式中：[X,Y,Z]和V=[VX,VY,VZ]分別是射彈的質(zhì)心位置和速度矢量;M和F分別是作用在射彈上的力矩和力;mp和I分別是射彈的質(zhì)量和慣性矩; [θ,φ,ψ]是射彈的旋轉(zhuǎn)角度，ω是射彈的角速度。

3 方法驗(yàn)證及模型設(shè)置

3.1 數(shù)值方法的驗(yàn)證

從兩方面對文中的數(shù)值模型進(jìn)行驗(yàn)證。首先是射彈入水產(chǎn)生的空泡與空泡經(jīng)驗(yàn)公式的對比。Savchenko等[19]通過大量試驗(yàn)總結(jié)了超空泡尺寸的經(jīng)驗(yàn)公式，具體表示如下：

(15)

圖1 模擬空泡與經(jīng)驗(yàn)公式曲線Fig.1 Comparison between calculated cavitation and empirical formula

此外，還計(jì)算了圓柱形射彈在380 m/s速度下小角度入水的空泡形態(tài)變化過程，并與文獻(xiàn)[7]中的實(shí)驗(yàn)結(jié)果進(jìn)行了對比，如圖2所示。數(shù)值計(jì)算的空泡變化、入水沖擊、水表面的擴(kuò)張、射彈尾拍和尾拍造成的空泡變化，可以看出試驗(yàn)拍攝的入水現(xiàn)象與計(jì)算結(jié)果基本一致。因此，使用的數(shù)值方法能以較高的精度模擬射彈的小角度入水問題。

圖2 空泡形態(tài)變化過程的計(jì)算結(jié)果與試驗(yàn)結(jié)果圖Fig.2 Comparison between calculation results and tests

3.2 計(jì)算模型

針對不同尾翼直徑的尾翼式超空泡射彈的小角度入水問題，進(jìn)行了入水角為10°、20°、30°的高速入水?dāng)?shù)值模擬研究。在Truscott試驗(yàn)[7-8]中可知，平頭空化器形狀的射彈入水性能較好，因此射彈空化器均采用平頭空化器形狀，如圖3所示。射彈空化器直徑為3 mm，尾翼直徑不同且為六尾翼結(jié)構(gòu)，整彈長徑比L/D=15。

射彈的質(zhì)心位置及轉(zhuǎn)動(dòng)慣量數(shù)據(jù)如表1所示。質(zhì)心位置以射彈彈底為基準(zhǔn)。射彈初始速度為1 000 m/s，初始自轉(zhuǎn)和章動(dòng)角速度均為零，3個(gè)射彈在均在指定角度和位置入水，總計(jì)9個(gè)工況。

接著是兩條五花龍睛，翩若驚鴻，宛若游龍。我照顧得無微不至。每天早起、回家，第一件事就是看它們是否活著。但最終還是死了，我有些難過，卻也大煞風(fēng)景地如釋重負(fù)。

圖3 不同尾翼直徑的射彈示意圖Fig.3 Projectiles with different tail diameters

表1 射彈結(jié)構(gòu)參數(shù)Table 1 Structural parameters of a projectile

3.3 邊界條件及求解設(shè)置

使用STAR-CCM+進(jìn)行數(shù)值模擬，計(jì)算方法基于有限體積法(FVM)，采用二階隱式格式來離散控制方程和VOF模型的擴(kuò)散對流項(xiàng)，利用SIMPLE算法的分離流求解器求解壓力和速度的耦合問題。為了清晰模擬液面濺射現(xiàn)象，需要對水面的網(wǎng)格局部加密。采用重疊網(wǎng)格和六自由度動(dòng)力學(xué)方法模擬射彈的入水過程。

計(jì)算域及邊界條件如圖4所示。其中θ為射彈的入水角，速度入口邊界速度設(shè)置為0，壓力出口邊界壓力與靜壓有關(guān)，采用靜水壓強(qiáng)的場函數(shù)進(jìn)行設(shè)置。圖5所示為射彈及重疊區(qū)域網(wǎng)格，對射彈空化器附近進(jìn)行局部加密，并在射彈外表面設(shè)置棱柱層以獲得高質(zhì)量的網(wǎng)格。

圖4 計(jì)算域示意圖(D=10 mm)Fig.4 Calculation domain diagram (D=10 mm)

圖5 計(jì)算網(wǎng)格圖Fig.5 Calculation grid diagram

3.4 網(wǎng)格無關(guān)性驗(yàn)證

對20 mm尾翼射彈入水進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證，建立3種不同數(shù)量的網(wǎng)格，數(shù)量分別為289萬(網(wǎng)格最小尺寸為0.5 mm)、538萬(網(wǎng)格最小尺寸為0.25 mm)、714萬(網(wǎng)格最小尺寸為0.125 mm)。入水角度為10°，如圖6所示為射彈入水0.5 ms時(shí)的速度曲線。由圖可知，最小網(wǎng)格尺寸減小到0.25 mm后，可以忽略因網(wǎng)格大小造成的計(jì)算誤差。綜合考慮計(jì)算精度和計(jì)算成本，采用的最小網(wǎng)格尺寸為0.25 mm(即case2對應(yīng)的網(wǎng)格規(guī)模)。整體網(wǎng)格數(shù)量因射彈結(jié)構(gòu)不同略有差異，但最小網(wǎng)格尺寸均控制在了0.25 mm。

圖6 不同網(wǎng)格數(shù)量的射彈入水速度曲線Fig.6 Calculation results of different grid numbers

4 結(jié)果與討論

4.1 入水空泡特性分析

圖7顯示了20 mm尾翼射彈在撞擊水面0.2 ms內(nèi)空泡的變化。射彈入水方向與水面成10°夾角，取射彈頭部與水面接觸時(shí)間為t=0。左側(cè)為水相體積分?jǐn)?shù)為10%的等值面圖，能清晰的看見空泡的輪廓，右側(cè)為射彈速度方向的后視角圖，可以看到射彈對水面沖擊造成的液體飛濺。t=0.05 ms時(shí)，頭部空化器沖擊水面形成空泡，空泡開始包裹射彈且處于擴(kuò)張階段。t=0.1 ms時(shí)，空泡已經(jīng)能包裹一半射彈，后部尾翼開始接觸水面，此時(shí)空泡尺寸還較小，尾翼刺穿空泡壁的深度較大，造成液面的二次飛濺。t=0.15 ms時(shí)，頭部空化器形成的空泡已經(jīng)包裹大部分射彈，尾翼刺穿空泡壁深度減小，兩者生成的空泡與頭部空化器形成的空泡未完全融合。空化器造成的水面濺射以及尾翼造成的二次飛濺加劇。t=0.2 ms時(shí)，空化器形成的空泡已經(jīng)能全部包裹住射彈，射彈只有頭部空化器處于沾濕狀態(tài)。

圖7 0.2 ms內(nèi)20 mm尾翼射彈入水示意圖Fig.7 Water entry diagram of 20 mm tail fin projectile in 0.2 ms

圖8為對應(yīng)的壓力云圖。從圖中可以看出：0.05 ms時(shí)，受到自由液面的影響，頭部空化器周圍壓力變化明顯且壓力場分布不對稱；前錐段下表面壓力有一定上升，這是因?yàn)樯鋸椷€沒有完全浸入水中，空泡發(fā)展尚不充分，前錐段下表面與水發(fā)生碰撞。

圖8 入水壓力云圖Fig.8 Cloud diagram of water inlet pressure

0.1 ms時(shí)空泡已包裹射彈主體，尾翼周圍壓力出現(xiàn)了較大變化，這是由于尾翼直徑較大，空泡尺寸不足以包裹住尾翼，尾翼與液體表面發(fā)生沖擊，如圖8(b)。0.15 ms時(shí)射彈繼續(xù)深入水面，由于空泡尺寸的變大，尾翼刺穿空泡壁的深度降低，如圖7所示，造成尾翼受到的壓力減小，表現(xiàn)為圖8(c)中尾翼周圍壓力場變化。0.2 ms時(shí)，空泡已充分發(fā)展并將射彈尾翼完全包裹，僅頭部受到了較大的壓力沖擊，并且此時(shí)頭部壓力場分布受水面影響減小而分布均勻。

圖9為尾翼直徑分別為15 mm和10 mm的射彈在0.1 ms和0.2 ms時(shí)的入水姿態(tài)水相等值面圖。與圖7對比，0.1 ms時(shí)，尾翼均開始觸碰到水面，由于尾翼大小的不同，10 mm尾翼幾乎沒有對水面造成形狀上的變化，15 mm尾翼也只是輕微改變了水面的形狀。結(jié)合圖8(b)和圖10中的0.1 ms時(shí)尾翼附近的壓力變化(圖中白色虛線圓圈標(biāo)注區(qū)域)，可以得出此時(shí)尾翼受到的壓力大小為20 mm尾翼>15 mm尾翼>10 mm尾翼。0.2 ms時(shí)，射彈均被空化器產(chǎn)生的空泡完全包裹，結(jié)合圖7還可看出尾翼激起的水面飛濺20 mm尾翼最大，15 mm尾翼次之，10 mm尾翼幾乎沒有。圖10壓力云圖中，0.2 ms時(shí)僅射彈的空化器周圍存在壓力場分布，并且與圖8的變化一致受水面影響減小而分布均勻。

圖9 15 mm(a)和10 mm(b)尾翼射彈 在0.1 ms和0.2 ms時(shí)的入水示意圖Fig.9 Water entry diagram of 15 mm (a) and 10 mm (b) tail fin projectiles at 0.1 ms and 0.2 ms

圖10 15 mm(a)和10 mm(b)尾翼射彈 在0.1 ms和0.2 ms時(shí)的壓力云圖Fig.10 Pressure nephogram of 15 mm (a) and 10 mm (b) tail projectiles at 0.1 ms and 0.2 ms

4.2 入水彈道特性分析

4.2.1阻力系數(shù)分析

(16)

式中：Fd為射彈受到的阻力；ρ為水的密度；v為射彈速度；r0為射彈半徑。

圖11為射彈在入水過程中受到的阻力系數(shù)曲線。由圖可以清晰的看到射彈在入水過程的阻力系數(shù)峰值以及其對應(yīng)的時(shí)刻。0.1 ms前尾翼沒有參與入水，3種射彈阻力系數(shù)變化一致，在頭部撞擊水面是出現(xiàn)第一次阻力系數(shù)峰值，約0.075。以圖11(a)為例，0.1～0.2 ms是尾翼參與入水的主要階段，20 mm尾翼射彈受到的阻力最大，阻力系數(shù)峰值為0.119，15 mm尾翼射彈的阻力系數(shù)峰值為0.080，10 mm尾翼射彈的阻力系數(shù)0.1 ms后基本不變。0.2 ms后，入水初始階段結(jié)束，阻力系數(shù)下降到了0.07左右，后續(xù)基本保持不變，這是因?yàn)榇藭r(shí)僅射彈頭部沾濕，沾濕面積達(dá)到最小，見圖7和圖9。但20 mm尾翼射彈在0.2 ms后阻力系數(shù)緩慢上升并在0.9 ms后繼續(xù)下降。這是因?yàn)樯鋸椀母┭鼋前l(fā)生變化導(dǎo)致尾翼刺穿上半部分空泡壁。圖12為0.8ms時(shí)射彈前側(cè)視圖的壓力云圖，可以看到，20 mm尾翼射彈的上半部分尾翼承受了較高的壓力，顯然已經(jīng)刺穿空泡壁，而10 mm和15 mm尾翼表面上壓力較低，說明其仍處于被空泡包裹狀態(tài)。同時(shí)射彈尾翼表面上不同的壓力變化，使得0.2 ms后20 mm尾翼射彈受到的阻力最大，10 mm尾翼射彈最小。符合圖12射彈在0.2 ms后阻力系數(shù)曲線的大小規(guī)律。

20°入水和30°入水的情況，如圖11(b)和圖11(c)，可以看出，同種射彈不同角度入水的阻力系數(shù)變化規(guī)律是一致的。對于20 mm尾翼射彈的阻力系數(shù)曲線，隨著入水角度的增加，其在0.2 ms內(nèi)的阻力系數(shù)峰值下降到了圖11(b)的0.106和圖11(c)的0.100，對應(yīng)于圖中的壓力云圖，其中白色虛線為自由液面的位置。0.2 ms后，另兩射彈并沒有因?yàn)槿胨堑脑黾佣写蟮淖兓?，阻力系?shù)穩(wěn)定維持在了0.07附近。

圖11 射彈入水阻力系數(shù)曲線Fig.11 Water entry resistance coefficient of a projectile

圖12 10°入水0.8 ms時(shí)射彈表面壓力云圖(前側(cè)視圖)Fig.12 Front side view of projectile surface pressure at 0.8 ms of 10° water entry

4.2.2速度分析

圖13(a)為射彈在10°入水隨時(shí)間變化的射彈速度曲線。由圖可以看出，不同的尾翼直徑對速度的影響較大。入水1 ms時(shí)，速度最高的是15 mm尾翼射彈，速度最低的則是 20 mm尾翼射彈。

本次計(jì)算中，射彈的轉(zhuǎn)動(dòng)慣量和質(zhì)心位置和質(zhì)量都是不同的，見表1。圖13(a)的放大部分，射彈在0.2 ms內(nèi)速度變化了2次。首先是0.1 ms內(nèi)的初期入水階段，此時(shí)有空化器和射彈前錐段與水接觸，3種射彈所受阻力基本相同，但由于質(zhì)量差異，20 mm尾翼射彈速度衰減最慢，10 mm尾翼射彈速度衰減最快。0.1 ms后射彈尾翼相繼入水，此時(shí)受尾翼影響，20 mm尾翼射彈所受阻力最大，導(dǎo)致其速度迅速衰減，變?yōu)樽畹停?5 mm尾翼射彈速度最快。

圖13(b)和圖13(c)為射彈20°入水和30°入水的曲線。隨著入水角度的增加，15 mm尾翼射彈的速度優(yōu)勢沒有變化。但是，20 mm與10 mm尾翼射彈的速度變化趨勢受角度影響。10°入水時(shí)，因?yàn)?0 mm尾翼射彈受到的阻力較大，10 mm尾翼射彈在0.2 ms后速度完全超過了前者；20°入水時(shí)，20 mm尾翼射彈的尾翼受到的阻力減小，速度曲線已經(jīng)緊貼10 mm尾翼射彈，僅在0.9 ms后，10 mm尾翼射彈的速度才稍微超過前者。30°入水的情況下，20 mm尾翼射彈的尾翼受到的阻力繼續(xù)減小，0.2 ms后10 mm尾翼射彈的速度曲線與前者完全重合。

4.2.3俯仰角分析

圖14(c)和圖14(d)分別為30°入水的俯仰角變化和俯仰力矩曲線。由圖可以看出，0.1 ms的初期入水，射彈空化器和前錐段使得射彈受到的俯仰力矩增加，射彈的俯仰角均因此減??；0.1 ms后尾翼開始接觸水面的影響，不同的尾翼規(guī)格使得俯仰運(yùn)動(dòng)出現(xiàn)了3種情況。首先是20 mm尾翼射彈，尾翼入水使其受到的相反的俯仰力矩最大，其俯仰角迅速增大。0.2 ms后尾翼入水的影響消失，俯仰力矩接近0，但由于大尾翼的影響，0.35 ms后射彈受到的俯仰力矩持續(xù)增加，這造成了射彈的俯仰角增加變緩且在0.9 ms后開始減??；15 mm尾翼射彈和10 mm尾翼射彈的情況類似，受到的俯仰力矩未能抵消初期入水的影響，射彈的俯仰角持續(xù)減小。其中15 mm尾翼射彈受到的影響更大，其俯仰角減小的速度更快。0.3 ms后兩者射彈受到的俯仰力矩均為0，直到最后射彈的俯仰運(yùn)動(dòng)不再發(fā)生改變。

圖13 射彈入水速度變化曲線Fig.13 Change of projectile entering water speed

圖14 射彈入水俯仰角及俯仰力矩曲線Fig.14 Pitching angle and moment of projectile entering water

如圖14(a)和圖14(b)，可以看出入水角度的變化不會(huì)影響射彈俯仰運(yùn)動(dòng)的運(yùn)動(dòng)趨勢，只會(huì)影響射彈俯仰運(yùn)動(dòng)的程度。隨著入水角度的增加，射彈俯仰運(yùn)動(dòng)程度變緩。從計(jì)算可以看出，15 mm尾翼射彈在入水的俯仰運(yùn)動(dòng)上具有明顯優(yōu)勢。

圖14(a)中20 mm尾翼射彈10°入水狀況，射彈俯仰角變化有往復(fù)運(yùn)動(dòng)的趨勢，將計(jì)算域擴(kuò)大，增加射彈運(yùn)動(dòng)的時(shí)間，得到圖15的俯仰力矩和俯仰角變化曲線?？梢钥吹?，射彈在初期入水影響后，俯仰力矩和俯仰角變化均呈現(xiàn)出振蕩變化規(guī)律，使得射彈的運(yùn)動(dòng)維持在一個(gè)相對穩(wěn)定的狀態(tài)。射彈形成尾拍后的最大俯仰角速度約150 rad/s，尾拍周期約為 2 ms，且4.5 ms俯仰角振蕩幅度逐漸擴(kuò)大的現(xiàn)象，與梁景奇[20]計(jì)算初速為1 000 m/s的尾翼射彈水中尾拍現(xiàn)象的俯仰角曲線趨勢一致，再次驗(yàn)證了計(jì)算的準(zhǔn)確性。

圖15 俯仰角與俯仰力矩曲線Fig.15 Angle and pitch moment

5 結(jié)論

對3種不同尾翼直徑的尾翼式超空泡射彈的小角度入水進(jìn)行了數(shù)值模擬，對空泡特性和彈道特性進(jìn)行分析，得出以下結(jié)論：

1) 0.1 ms前射彈入水狀態(tài)是一致的，0.1～0.2 ms之間，尾翼接觸到水面后，激起了水的飛濺并阻礙射彈的運(yùn)動(dòng)。0.2 ms之后，射彈被空化器產(chǎn)生的空泡完全包裹，此時(shí)受到的阻力大大減少。

2) 射彈速度受到射彈結(jié)構(gòu)參數(shù)以及尾翼與水的沖擊的影響，15 mm尾翼射彈存速能力最強(qiáng)。

3) 在不同入水角情況下，15 mm尾翼射彈的入水俯仰運(yùn)動(dòng)都最為平穩(wěn)。

4) 相同結(jié)構(gòu)射彈，隨著入水角度的增加，射彈的俯仰運(yùn)動(dòng)變緩。