趙博博，劉榮忠，郭 銳，張 迪，袁 軍，陳 亮

(1.南京理工大學(xué) 機(jī)械工程學(xué)院，南京 210094；2.中國(guó)人民解放軍73917部隊(duì)，南京 290014；3.中國(guó)人民解放軍陸軍軍官學(xué)院，合肥 230000)

?

扭曲尾翼彈箭的馬格努斯數(shù)值研究

趙博博1, 2，劉榮忠1，郭 銳1，張 迪3，袁 軍1，陳 亮1

(1.南京理工大學(xué) 機(jī)械工程學(xué)院，南京 210094；2.中國(guó)人民解放軍73917部隊(duì)，南京 290014；3.中國(guó)人民解放軍陸軍軍官學(xué)院，合肥 230000)

為提升掠飛攻頂彈箭較高轉(zhuǎn)速下的飛行穩(wěn)定性，運(yùn)用數(shù)值計(jì)算方法研究了彈體-扭曲尾翼組合體在飛行過(guò)程中的馬格努斯效應(yīng)氣動(dòng)機(jī)理，并應(yīng)用標(biāo)準(zhǔn)尾翼彈(BFM)模型的實(shí)驗(yàn)數(shù)據(jù)對(duì)數(shù)值方法進(jìn)行了驗(yàn)證。分別研究了帶有平板尾翼和扭曲尾翼彈體模型的馬格努斯力和力矩隨攻角的變化規(guī)律，并針對(duì)彈體彈翼組合體產(chǎn)生馬格努斯效應(yīng)的機(jī)理深入分析。結(jié)果表明，扭曲尾翼可有效改善翼面的壓力分布，并降低彈體對(duì)翼面馬格努斯效應(yīng)的干擾，在大攻角時(shí)其表現(xiàn)更勝一籌；彈體所受馬格努斯力較大，主要集中在受到渦對(duì)稱畸變的尾錐部；尾翼主要由于彈體干擾以及幾何外形的影響馬格努斯力集中在尾部，兩者產(chǎn)生的馬格努斯力矩?cái)?shù)值相差不大，但方向相反。

流體力學(xué)；扭曲尾翼；旋轉(zhuǎn)尾翼彈丸；馬格努斯效應(yīng)；數(shù)值仿真

0 引言

掠飛攻頂末敏彈是一種采用非直瞄方式向裝甲目標(biāo)方位發(fā)射的新型末敏彈，應(yīng)用旋轉(zhuǎn)的彈體作為掃描平臺(tái)，彈丸一邊飛行一邊高速旋轉(zhuǎn)，敏感器對(duì)前側(cè)方進(jìn)行探測(cè)掃描，掃描區(qū)域?yàn)橹螀^(qū)，一旦探測(cè)到目標(biāo)，即按一定的決策規(guī)則起爆MEFP戰(zhàn)斗部，所形成的多枚彈丸從裝甲目標(biāo)上方擊毀目標(biāo)[1-2]。旋轉(zhuǎn)的彈體不僅可有效降低由質(zhì)量偏心、尾翼加工誤差、推力偏心等因素的引起的彈丸散布，還作為敏感探測(cè)器的掃描平臺(tái)[3]。普通尾翼彈轉(zhuǎn)速僅僅為30 r/s以內(nèi)，掠飛攻頂末敏彈則要求轉(zhuǎn)速為其2～3倍[4]。但隨著轉(zhuǎn)速的增加，旋轉(zhuǎn)產(chǎn)生的馬格努斯效應(yīng)不再可忽略，將引起彈丸攻角增加，甚至破壞彈丸靜穩(wěn)定性，導(dǎo)致飛行失穩(wěn)[5]，扭曲尾翼作為一種新型高轉(zhuǎn)速尾翼彈結(jié)構(gòu)，具有增加轉(zhuǎn)速、降低阻力的優(yōu)勢(shì)，對(duì)其開(kāi)展研究也未見(jiàn)于報(bào)道。由此可見(jiàn)，研究馬格努斯力與力矩對(duì)較高轉(zhuǎn)速要求的掠飛末敏彈具有重要意義。

旋轉(zhuǎn)彈馬格努斯效應(yīng)是由彈體表面附面層內(nèi)復(fù)雜三維流動(dòng)產(chǎn)生，引起彈體相反的兩面壓力不對(duì)稱[6]，其不僅與彈體轉(zhuǎn)速和攻角相關(guān)，還與彈丸質(zhì)心、彈體外形、馬赫數(shù)、雷諾數(shù)等有關(guān)；尾翼穩(wěn)定彈在零轉(zhuǎn)速時(shí)，即存在馬格努斯效應(yīng)，這是由尾翼的差動(dòng)配置與彈體之間相互影響產(chǎn)生[7]；過(guò)去對(duì)馬格努斯效應(yīng)的研究主要通過(guò)實(shí)驗(yàn)的方法，文獻(xiàn)[8-9]通過(guò)風(fēng)洞實(shí)驗(yàn)方法研究了標(biāo)準(zhǔn)尾翼模型(BFM)在超音速、跨音速下的馬格努斯力與力矩，指出尾翼彈馬格努斯力具有非線性的特性，在跨音速時(shí)受尾翼影響較大；文獻(xiàn)[10]則研究了旋轉(zhuǎn)彈在不同馬赫數(shù)與轉(zhuǎn)速下的馬格努斯效應(yīng)，但風(fēng)洞實(shí)驗(yàn)方法對(duì)設(shè)備要求非常高、周期長(zhǎng)、花費(fèi)昂貴；數(shù)值仿真在附面層的畸變，彈體受力分析上具有直觀、清晰的優(yōu)勢(shì)，但由于馬格努斯力相對(duì)較小，數(shù)值仿真較難實(shí)現(xiàn)[11-12]；文獻(xiàn)[13]對(duì)BFM模型進(jìn)行數(shù)值仿真，并與美國(guó)阿諾德實(shí)驗(yàn)中心風(fēng)洞數(shù)據(jù)進(jìn)行對(duì)比，成功實(shí)現(xiàn)了平均雷諾下的N-S方程對(duì)大攻角條件下尾翼彈進(jìn)行數(shù)值模擬；文獻(xiàn)[14]通過(guò)數(shù)值和實(shí)驗(yàn)相對(duì)比的方法，驗(yàn)證了數(shù)值方法對(duì)旋轉(zhuǎn)彈模擬的精確度；文獻(xiàn)[15-16]也對(duì)馬格努斯效應(yīng)開(kāi)展了一定研究。

將掠飛攻頂末敏彈簡(jiǎn)化為帶有尾桿結(jié)構(gòu)的扭曲尾翼彈丸，由于尾桿直徑遠(yuǎn)小于彈體直徑，導(dǎo)致風(fēng)洞實(shí)驗(yàn)特別是縮比實(shí)驗(yàn)受到尾桿直徑的限制而難以實(shí)現(xiàn)。所以，對(duì)此類結(jié)構(gòu)的尾翼彈實(shí)驗(yàn)較少。本文通過(guò)數(shù)值仿真的方法，采用滑移網(wǎng)格技術(shù)，先對(duì)BFM模型數(shù)值模擬，與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，進(jìn)行算法驗(yàn)證，然后計(jì)算了扭曲尾翼彈丸的馬格努斯力與力矩，并深入分析，對(duì)新型尾翼結(jié)構(gòu)彈丸設(shè)計(jì)具有較大參考意義。

1 數(shù)值方法

數(shù)值計(jì)算方法具有清晰、直觀的優(yōu)勢(shì)，能從機(jī)理上解釋馬格努斯效應(yīng)產(chǎn)生的原因。

1.1 掠飛末敏彈模型

模型為帶有尾桿結(jié)構(gòu)的扭曲尾翼彈丸，如圖1所示。尾翼彈丸總長(zhǎng)為6.5D，頭部長(zhǎng)為1.5D，頭部為圓弧形母線，半徑為3.2D，圓柱部長(zhǎng)2.4D，尾錐錐角為45°，尾桿直徑為0.4D、長(zhǎng)度為2.05D，尾翼翼展為2.72D，翼弦為0.28D，厚度為0.024D。彈體轉(zhuǎn)動(dòng)方向?yàn)閺膹椢蚕驈楊^方向看為順時(shí)針旋轉(zhuǎn)。

圖1 “十字形”布局扭曲尾翼彈丸結(jié)構(gòu)圖Fig.1 The cruciform layout twist fin projectile model

圖2中，扭曲尾翼為一種沿翼展均勻扭曲的尾翼，在翼根處安裝攻角較小，沿翼展安裝攻角均勻增大。引入扭曲率K(y)和平均攻角εm表示尾翼的幾何特征，定義如下：扭曲率K(y)表征沿翼展延伸方向尾翼面的扭曲程度，將安裝角δ(y)對(duì)翼展y求偏導(dǎo)，得到扭曲率為

K(y)=?δ(y)/?y

(2)

平均攻角εm的正弦值等于扭曲尾翼在yoz平面上的投影面積與尾翼側(cè)表面積的比值，是扭曲尾翼與斜置平板尾翼相對(duì)比的重要參考量，即

εm=arcsin(S1/Sw)

(3)

式中S1為尾翼面在垂直于彈軸平面上的投影面積；Sw為翼面的側(cè)表面積。

圖2 扭曲尾翼幾何模型Fig.2 The geometric model of twisted fin

采用笛卡爾右手坐標(biāo)系，坐標(biāo)原點(diǎn)在彈尾中心，x軸正向從彈頭部指向彈尾部，y軸指向上方，z軸垂直于xoy平面，指向外方向?？芍?，當(dāng)攻角向上，尾翼彈丸順時(shí)針旋轉(zhuǎn)時(shí)，根據(jù)馬格努斯力的經(jīng)典解釋，旋轉(zhuǎn)彈的馬格努斯力指向z軸正向。

取兩個(gè)模型計(jì)算：第一個(gè)模型為特殊的扭曲尾翼，即扭曲率為0的平板尾翼，由平板尾翼偏置13°，記為F1；第二個(gè)模型為翼根處攻角小，翼梢處攻角大，沿翼展方向扭曲率為0.11°/mm的扭曲尾翼，扭曲尾翼的平均攻角仍為13°，記為F2。

1.2 網(wǎng)格劃分及邊界條件

滑移網(wǎng)格技術(shù)是一種易實(shí)現(xiàn)的動(dòng)網(wǎng)格技術(shù)，具有計(jì)算速度快、占用內(nèi)存少、精度高等優(yōu)點(diǎn)。將流場(chǎng)劃分為內(nèi)外2個(gè)具有交接面的區(qū)域，滑移網(wǎng)格技術(shù)通過(guò)設(shè)置內(nèi)外2個(gè)區(qū)域的不同轉(zhuǎn)速，實(shí)現(xiàn)不同域內(nèi)網(wǎng)格的繞軸轉(zhuǎn)動(dòng)。馬格努斯效應(yīng)主要是由邊界層的復(fù)雜三維流動(dòng)引起。所以，網(wǎng)格的質(zhì)量是影響求解結(jié)果的關(guān)鍵因素。如圖3(a)所示，將流場(chǎng)劃分為內(nèi)外兩個(gè)區(qū)域，采用圓柱形流場(chǎng)區(qū)域，以防止激波溢出。為減少網(wǎng)格數(shù)目，采取內(nèi)部流場(chǎng)區(qū)域網(wǎng)格較外部區(qū)域密，交接面上節(jié)點(diǎn)不用對(duì)齊。為了得到準(zhǔn)確的翼面氣動(dòng)力特性、以及邊界層的復(fù)雜流動(dòng)，使用O型拓?fù)浣Y(jié)構(gòu)對(duì)彈體周圍加密。如圖3(b)、(c)所示，壁面第一層網(wǎng)格的厚度采用0.002 mm，以滿足y+<1，隨著網(wǎng)格向外延伸，控制網(wǎng)格厚度增長(zhǎng)率小于1.3，以保證網(wǎng)格的均勻增長(zhǎng)。

(a)總體網(wǎng)格結(jié)構(gòu) (b)彈頭部網(wǎng)格結(jié)構(gòu) (c)尾翼部網(wǎng)格結(jié)構(gòu)

湍流模型采用具有更高可信度和精度的剪切壓力運(yùn)輸(SST)k-ω模型，該模型結(jié)合了k-ε模型與k-ω模型的優(yōu)點(diǎn)，將混合功能和雙模型疊加在一起，此模型避免了應(yīng)用壁面函數(shù)時(shí)對(duì)網(wǎng)格的y+要求較低的近似方法，而是直接對(duì)附面層內(nèi)進(jìn)行模擬，對(duì)網(wǎng)格要求很高。

采用有限體積法對(duì)空間進(jìn)行離散，采取基于密度的隱式耦合算法，時(shí)間項(xiàng)采用穩(wěn)態(tài)方法來(lái)計(jì)算模型在給定狀態(tài)下的氣動(dòng)特性，采用迎風(fēng)型矢通量分裂格式(AUSM+)計(jì)算對(duì)流和流體壓力，AUSM+ 格式兼有Roe格式的間斷高分辨率和VanLeer格式的計(jì)算效率，而且克服了二者的缺點(diǎn)。

遠(yuǎn)場(chǎng)入口邊界采用壓力遠(yuǎn)場(chǎng)，設(shè)置來(lái)流值，遠(yuǎn)場(chǎng)周向邊界采用自由流條件，外部區(qū)域不轉(zhuǎn)動(dòng)，內(nèi)部區(qū)域采用網(wǎng)格移動(dòng)條件，設(shè)置轉(zhuǎn)速。彈體采用壁面邊界條件，邊界采用無(wú)滑移邊界條件，壁面采用運(yùn)動(dòng)方式，運(yùn)動(dòng)為旋轉(zhuǎn)方法，旋轉(zhuǎn)軸為彈軸，壁面與內(nèi)部運(yùn)動(dòng)區(qū)相關(guān)聯(lián)。

2 網(wǎng)格收斂性及算法驗(yàn)證

采用標(biāo)準(zhǔn)尾翼彈模型(BFM)進(jìn)行算法驗(yàn)證[9]，幾十年來(lái)此模型在國(guó)外被作為參考模型，已經(jīng)積累了大量風(fēng)洞實(shí)驗(yàn)數(shù)據(jù)，是較成熟的研究彈型，其尺寸為模型尺寸為彈身直徑D=45 mm，彈長(zhǎng)為10D，錐形頭部角度20°，翼展3D，弦長(zhǎng)為D，翼寬0.08D，質(zhì)心距離彈頭部6.1D，尾翼偏置角度為3°，求解不同條件下的馬格努斯力及力矩，進(jìn)行算法的有效性驗(yàn)證。計(jì)算中，分別采取130萬(wàn)、280萬(wàn)、430萬(wàn)的網(wǎng)格數(shù)目進(jìn)行求解，驗(yàn)證網(wǎng)格的收斂性。結(jié)果表明，在網(wǎng)格數(shù)目為280萬(wàn)時(shí)，滿足精度要求。

圖4 BFM模型仿真結(jié)果對(duì)比Fig.4 The simulation comparison with experimental result

將仿真數(shù)據(jù)與風(fēng)洞實(shí)驗(yàn)數(shù)據(jù)相比較[10]。圖4為在2種條件下，BFM馬格努斯力數(shù)值仿真值與實(shí)驗(yàn)值對(duì)比結(jié)果?？芍?，數(shù)值仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)相吻合。其中，對(duì)于Ma=1.1、ω*=0.025時(shí)，吻合較好。對(duì)于Ma=2.5、轉(zhuǎn)速較低為ω*=0.025時(shí)，誤差較大，最大誤差為25%，但基本規(guī)律相吻合?？傻肅FD方法在求解小攻角，高轉(zhuǎn)速的情況下有較高的精確度，驗(yàn)證了算法的正確性。

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

計(jì)算條件：Ma=1.5，總溫T0=300 K，攻角α=0°～12°，無(wú)量綱轉(zhuǎn)速ω*=ωr/v∞=0.122(其中ω為彈丸轉(zhuǎn)速，r為彈體最大半徑，v∞為來(lái)流速度)，采用氣體為理想氣體。分別計(jì)算模型F1，F(xiàn)2，模型F1尾翼為平板尾翼，模型F2尾翼為相同平均偏置角下的扭曲尾翼。

3.1 整體結(jié)果分析

圖5為模型尾翼、彈體及組合體馬格努斯力系數(shù)隨攻角的變化曲線。對(duì)于模型F1與F2，由于彈體外形及飛行條件完全一樣，且在尾翼后端受尾翼干擾的彈體區(qū)域較小。所以，認(rèn)為2種模型彈體的馬格努斯效應(yīng)完全一致。模型彈體受到z軸正向的馬格努斯力，隨著攻角增加，當(dāng)攻角大于一定值后變化不大；模型F1、F2尾翼上馬格努斯力比彈體小得多，與彈體反向指向z軸負(fù)向，且F1的尾翼隨著攻角增大馬格努斯力出現(xiàn)換向，F(xiàn)2的尾翼上的力隨攻角一直增加；組合體的馬格努斯力為小攻角時(shí)模型F1大于F2 ，較大攻角時(shí)相反，組合體總體小于單獨(dú)彈體的馬格努斯力。

圖6為模型尾翼、彈體及組合體的馬格努斯力矩變化曲線。由圖6可知，尾翼的馬格努斯力矩和彈體數(shù)值差別不大，但與彈體呈反向；在小攻角時(shí)，模型F2的馬格努斯力矩略大于模型F1，在大攻角時(shí)，遠(yuǎn)小于F1，這說(shuō)明扭曲尾翼能有效改善模型的馬格努斯力矩。這是由于F2尾翼上的馬格努斯力矩隨攻角先增大、后減小，但減小幅度不大，始終與彈體保持反向；F1尾翼的馬格努斯力矩在較大攻角換向，提供與彈體馬格努斯力矩方向一致的力，導(dǎo)致組合體力矩增加。

圖5 模型F1、F2的馬格努斯力系數(shù)的變化曲線Fig.5 Change of Magnus force coefficient of model F1,F2

圖6 模型F1、F2的馬格努斯力矩變化曲線Fig.6 Change of Magnus moment coefficient of model F1,F2

3.2 尾翼馬格努斯效應(yīng)分析

已知平板尾翼的馬格努斯力隨著攻角的增加呈現(xiàn)反向的趨勢(shì)，而扭曲尾翼馬格努斯力隨著攻角的增大方向不變。針對(duì)尾翼上馬格努斯效應(yīng)產(chǎn)生機(jī)理單獨(dú)分析。

圖7為翼面中點(diǎn)處的赤道面與上下翼面的交線上的壓力差沿y軸的分布，D為彈體直徑。圖9中，左邊為尾翼3，右端為尾翼1。由圖7、圖8知，對(duì)于F1的尾翼1翼根處零攻角時(shí)產(chǎn)生z軸正向的力，尾翼3翼根處產(chǎn)生z軸負(fù)向力，有攻角時(shí)打破平衡，尾翼3翼根上壓強(qiáng)更高，提供z軸正向力；同樣，在模型F1的尾翼1、3的翼梢處，由于攻角引起的壓力變化，使尾翼3壓強(qiáng)高受到z軸正向力，兩者相互抵消，導(dǎo)致模型F1的馬格努斯力隨著攻角的變化先增大，后反向。對(duì)于模型F2的尾翼3，由于扭曲翼面上實(shí)際攻角均勻，平衡轉(zhuǎn)速下幾乎為零，在有攻角δ時(shí)，彈體對(duì)尾翼3翼面的干擾帶來(lái)的壓力變化并不影響翼面壓力分布；而模型F2的尾翼1在翼根處受到尾桿的干擾，攻角越大產(chǎn)生的渦越強(qiáng)烈，對(duì)尾翼1翼根處產(chǎn)生越大的z軸負(fù)向力。所以，模型F2的尾翼馬格努斯力隨著攻角增加而增加。可見(jiàn)，尾翼的馬格努斯力產(chǎn)生的原因主要是因?yàn)閺楏w對(duì)尾翼的干擾，有攻角時(shí)，扭曲尾翼同樣改善壓力隨翼展的分布。

圖7 δ=9°時(shí)尾翼1、3翼面z方向壓強(qiáng)分布Fig.7 z direction component of pressure at wing 1,3

尾翼上馬格努斯力矩一部分由尾翼1、3的馬格努斯力產(chǎn)生，另一部分源于差動(dòng)配置的尾翼2、4翼面上壓強(qiáng)在x軸的分量。

圖8為翼面中點(diǎn)處赤道面與尾翼上下翼面的交線上的壓力差及其在x軸(對(duì)安裝角取正弦)上投影沿z軸的分布。右邊為尾翼4，左邊為尾翼2，D為彈體直徑。尾翼2、4的壓強(qiáng)主要源于安裝角、彈體攻角、以及轉(zhuǎn)速引起的附加攻角共同作用。模型F1為平板尾翼，相比較于模型F2的扭曲尾翼，則翼根處安裝角較大，翼梢處安裝角小。所以，模型F1尾翼4在翼根處壓力大于模型F2，在翼梢處小于模型F2；模型F1的尾翼2在翼根處提供較大的向下的安裝角。所以，在翼根處壓力較小，在翼梢處的壓強(qiáng)較大。對(duì)比可知，尾翼4上由于彈體攻角與安裝角同向，導(dǎo)致尾翼4壓力分布最不均勻，扭曲尾翼同樣改善翼面上的壓力分布。

圖8 δ=9°時(shí)尾翼2、4翼面壓強(qiáng)及翼面x方向壓強(qiáng)分布Fig.8 Pressure distribution and x direction component on the surface of the wing 2,4

由于圖8中尾翼2、4壓力差在x軸上分量方向相反，所以產(chǎn)生y軸正向的力矩。模型F2的尾翼2在翼梢處壓力雖然不大，但此處尾翼安裝角較大，仍存在較大的x軸正向分力，尾翼4同樣存在較大x軸負(fù)向分力。所以，模型F2尾翼的馬格努斯力矩大于模型F1，主要是由于扭曲尾翼翼梢處安裝角較大而導(dǎo)致；模型F2尾翼的馬格努斯力矩與彈體相反，達(dá)到降低組合體的馬格努斯力矩的效果。

3.3 彈體馬格努斯效應(yīng)分析

彈體的馬格努斯效應(yīng)同樣是影響彈丸馬格努斯效應(yīng)的重要因素。單獨(dú)彈體的馬格努斯效應(yīng)主要由于空氣的粘性作用產(chǎn)生，邊界層位移厚度的非對(duì)稱畸變、邊界層轉(zhuǎn)捩非對(duì)稱畸變、徑向壓力梯度的非對(duì)稱畸變、周向切應(yīng)力的非對(duì)稱畸變和體渦的非對(duì)稱畸變是馬格努斯效應(yīng)產(chǎn)生的重要原因[6]。

下面對(duì)主要產(chǎn)生彈體馬格努斯效應(yīng)的主要因素進(jìn)行分析。

3.3.1 彈體流場(chǎng)曲線

彈體的馬格努斯效應(yīng)同樣是影響彈丸馬格努斯效應(yīng)的重要因素。圖9為彈丸周圍流場(chǎng)的流線圖，在小攻角時(shí)，流體沿頭部、圓柱部的壁面流動(dòng)，由層流逐漸轉(zhuǎn)化為湍流，在尾錐部形成明顯的渦，此時(shí)渦仍較對(duì)稱；而大攻角時(shí)，流體在背風(fēng)面渦會(huì)脫離附面層，特別是尾錐部的渦明顯呈現(xiàn)非對(duì)稱?？梢?jiàn)，在不同攻角時(shí)，導(dǎo)致馬格努斯效應(yīng)的機(jī)理不同。

3.3.2 彈體圓柱部

圖10為Ma=1.5旋轉(zhuǎn)彈丸和無(wú)旋彈丸在圓柱部(x=-3.5D)赤道面彈體周圍速度等值線，彈體部不旋轉(zhuǎn)時(shí)，邊界層相對(duì)于彈體攻角平面對(duì)稱，且下部等值線較密，上部稀疏；當(dāng)彈體旋轉(zhuǎn)后左側(cè)等值線變薄，速度增加，由于超音速下速度大壓強(qiáng)小，所以左側(cè)壓強(qiáng)減；右側(cè)等值線變稀疏，說(shuō)明速度減小，壓強(qiáng)增大，在彈體圓柱部產(chǎn)生z軸正向的力。這是因?yàn)樵诠ソ堑臈l件下，彈體受到由下向上的橫流，由于氣體本身具有粘性，當(dāng)彈體旋轉(zhuǎn)時(shí)帶動(dòng)周圍流場(chǎng)沿彈體旋轉(zhuǎn)產(chǎn)生環(huán)流。由于彈丸左側(cè)環(huán)流速度與橫流速度相同，所以氣流速度增大，且彈體附面層變??； 右側(cè)環(huán)流速度與橫流速度相反，流速降低，則附面層變厚，壓強(qiáng)增大，即產(chǎn)生z軸正向的側(cè)向力。

(a)δ=2° (b)δ=9°

(a)ω*=0 (b)ω*=0.122

圖11為Ma=1.5、δ=9°時(shí)，無(wú)旋彈丸和旋轉(zhuǎn)彈丸在圓柱部(x=-4.5D)與赤道面交線的壓力分布曲線。在彈體不旋轉(zhuǎn)時(shí)兩側(cè)壓力對(duì)稱，且迎風(fēng)面壓力遠(yuǎn)大于背風(fēng)面壓力，此時(shí)彈丸受到強(qiáng)烈的正升力；旋轉(zhuǎn)彈丸周向壓力分布不再對(duì)稱，主要表現(xiàn)為在彈丸在背風(fēng)面右側(cè)的壓力大于左側(cè)的壓力，使彈丸受到z軸正向的壓力。這主要由于邊界層位移厚度的非對(duì)稱畸變、邊界層轉(zhuǎn)捩非對(duì)稱畸變引起；在迎風(fēng)面左側(cè)壓力大于右側(cè)，彈丸受到z軸負(fù)向的壓力，主要由于周向切應(yīng)力的非對(duì)稱畸變以及氣體離心力的作用，但彈丸總體受到z軸正向的壓力。

3.3.3 尾錐部馬格努斯效應(yīng)

圖12為流場(chǎng)在尾錐部(x=-2.3D)的赤道面內(nèi)流線及密度云圖?？芍鲌?chǎng)在尾錐部流動(dòng)非常復(fù)雜，但主要受到體渦的非對(duì)稱畸變及邊界層的畸變的影響。圖12(a)為旋轉(zhuǎn)彈體在攻角為零時(shí)，尾錐部的流線及密度云圖。此時(shí)，由于沒(méi)有攻角，彈體附面層內(nèi)有均勻環(huán)流，帶動(dòng)周圍流場(chǎng)均勻旋轉(zhuǎn)，由于流場(chǎng)在尾錐部受到膨脹波的影響，速度增加密度減小，使密度呈現(xiàn)梯度變化，密度也是均勻軸對(duì)稱，此時(shí)沒(méi)有側(cè)向力。圖12(b)為彈丸在攻角為9°不旋轉(zhuǎn)時(shí)的密度云圖。此時(shí)，迎風(fēng)面密度高，背風(fēng)面密度低，云圖不在是軸對(duì)稱，而是相對(duì)于攻角平面對(duì)稱，仍沒(méi)有側(cè)向力，流線不再是均勻的環(huán)流，呈現(xiàn)面對(duì)稱的流場(chǎng)，同時(shí)在彈體背風(fēng)面產(chǎn)生兩個(gè)對(duì)稱的渦，這是因?yàn)榱骶€脫體，在背風(fēng)面與附面層分離。圖12(c)為旋轉(zhuǎn)彈丸在攻角為9°時(shí)的密度云圖。此時(shí)，旋轉(zhuǎn)彈體背風(fēng)面內(nèi)的流動(dòng)分離產(chǎn)生的渦已經(jīng)嚴(yán)重非對(duì)稱。其中，左側(cè)的渦由于旋轉(zhuǎn)向彈體靠近，已依附在彈體上，而右側(cè)的渦順時(shí)針?lè)较蛞苿?dòng)，遠(yuǎn)離彈體。由于渦的中心壓強(qiáng)較低，外圍壓強(qiáng)高，所以左側(cè)渦依附在彈體上形成低壓區(qū)，右側(cè)的渦低壓區(qū)遠(yuǎn)離彈體，彈體表面為高壓區(qū)，產(chǎn)生z軸正向力。

圖11 彈體在圓柱部(x=-4.5D)處周向壓力分布Fig.11 Computed Circumferential surface pressure distribution δ=9°, x=-4.5D

(a)δ=0°,ω*=0.122 (b)δ=9°,ω*=0 (c)δ=9°,ω*=0.122

4 結(jié)論

(1)尾翼上的馬格努斯力主要由彈體存在攻角時(shí)，彈體對(duì)氣流的干擾，流過(guò)尾翼的上下翼片所引起，導(dǎo)致上下翼片壓力改變。平板尾翼翼面上翼根、翼梢受力不均衡，更容易受壓力改變的影響，馬格努斯力變化劇烈；扭曲尾翼能有效改善翼面上的實(shí)際攻角，從而改善壓力分布，增強(qiáng)了翼面的抗干擾性，馬格努斯力變化規(guī)律，在大攻角時(shí)表現(xiàn)更勝一籌。

(2)彈體整體受到z軸正向的馬格努斯力，圓柱部由于主要為層流的影響、尾桿處雖然為湍流，但直徑較小，馬格努斯效力并不大；尾錐處有較強(qiáng)的渦畸變產(chǎn)生，此處馬格努斯力最大，但由于此處靠近彈丸質(zhì)心，彈體整體的馬格努斯力矩并不大，這為掠飛攻頂末敏彈氣動(dòng)優(yōu)化提供依據(jù)。尾翼處遠(yuǎn)離質(zhì)心，受到z軸負(fù)向的馬格努斯力，雖然作用在尾翼上馬格努斯力并不大，但由于力臂長(zhǎng)，所以尾翼馬格努斯力矩和彈體上數(shù)值相似，方向相反。

(3)通過(guò)數(shù)值仿真的方法，對(duì)彈體彈翼組合體的馬格努斯效應(yīng)產(chǎn)生機(jī)理，指出了彈體、彈翼的馬格努斯效應(yīng)主要影響因素并進(jìn)行分析，這對(duì)降低掠飛攻頂末敏彈的馬格努斯效應(yīng)及彈丸氣動(dòng)優(yōu)化具有重要意義。

[1] 楊紹卿. 靈巧彈藥工程[M]. 北京: 國(guó)防工業(yè)出版社, 2010.

[2] Megson, Thomas H.Aircraft structures for engineering stud-ents[M]. Access Online via Elsevier, 2012.

[3] 胡志鵬, 劉榮忠, 郭銳. 兩種典型尾翼形狀對(duì)無(wú)傘末敏彈氣動(dòng)特性的影響[J]. 南京理工大學(xué)學(xué)報(bào), 2012, 36(5).

[4] 趙博博, 劉榮忠, 郭銳. 扭曲尾翼飛行器氣動(dòng)特性數(shù)值研究[J].國(guó)防科技大學(xué)學(xué)報(bào),2014(3).

[5] Oh S Y, Kim S C, Lee D K, et al. Magnus and spin-damping measurements of a spinning projectile using design of experiments[J]. Journal of Spacecraft and Rockets, 2010, 47(6): 974-980.

[6] Schlichting H. Boundary-layer theory[M]. McGraw-hill Bo-ok Company, 1979.

[7] Morote J, Liao G. Prediction of nonlinear rolling and magnus coefficients of cruciform-finned missiles[J]. Journal of Aircraft, 2010, 47(4): 1413-1425.

[8] Arnan Seginer, Izhak Rosenwasser. Magnus effects on spinning transonic missiles[R]. AIAA 83-246.

[9] Jenke L M. Experimental roll-damping magnus, and static stability characteristics of two slender missile configurations at high angles of attack (0 to 90 deg)and Mach numbers 0.2 through 2.5[R]. AEDC-TR 76-58.

[10] Sturek W B, Dwyer H A, Kayser L D. Computations of magnus effects for a yawed, spinning body of revolution[J]. AIAA Journal, 1978, 16(7): 687-692.

[11] Sahu J. Numerical computations of dynamic derivatives of a finned projectile using a time accurate CFD method[R]. AIAA 2007.

[12] James Despirito. CFD prediction of magnus effect in subsonic to supersonic flight[J]. Army Research Laboratory, ARL-TR-4929, 2009.

[13] Bhagwandin V, Sahu J. Numerical prediction of dynamic stability derivatives for finned projectiles[C]//Ballistics 2013: 27th International Symposium on Ballistics. DEStech Publications, Inc, 2013: 186.

[14] Klatt D, Hruschka R, Leopold F. Numerical and experimental investigation of the magnus effect in supersonic Flows[R]. AIAA 2012-3230.

[15] 薛幫猛，楊永. 旋轉(zhuǎn)彈丸馬格努斯力數(shù)值計(jì)算[J]. 彈箭與制導(dǎo)學(xué)報(bào)，2005，25(2): 85-87.

[16] 王智杰，陳偉芳，李浩. 旋轉(zhuǎn)彈丸空氣動(dòng)力特性數(shù)值解法[J]. 國(guó)防科技大學(xué)學(xué)報(bào), 2003, 25(4): 15-19.

(編輯：崔賢彬)

Numerical prediction of the Magnus effect for twist fin swept flight projectile

ZHAO Bo-bo1, 2，LIU Rong-zhong1，GUO Rui1，ZHANG Di3，YUAN Jun1，CHEN Liang1

(1.School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China; 2.73917 Troops of the Chinese People's Liberation Army, Nanjing 210094, China; 3.The Chinese People's Liberation Army Military Academy, Hefei 230000, China)

In order to improve flight stability of swept flight assault roof (SFAR)projectile under high rotation speed, Magnus effect pneumatic mechanism of projectile bodies-twist fin combination in flight was studied by using numerical method and the numerical method was validated by using experimental data of BFM model.On the basis of using standard tail experimental data to validate the numerical method, the variation of Magnus effect with the changing rule of attack angle of flat tail and twist tail projectile model was studied. In view of the projectile wing assembly, Magnus effect mechanism was analyzed thoroughly. The results show that twist fin can effectively improve the pressure distribution of wing surface and reduce Magnus effect interference of projectile bodies to wing surface, especially for big attack angle; larger projectile body Magnus force is mainly focused on coccygeal vertebra which is distorted of the vortex symmetry closed to centroid ; Due to projectile interference and geometric profile effect ,empennage lateral force is mainly concentrated on tail.The values of torgue are not much different from each other, while the directions are quite the contrary.

fluid mechanics；twist fin；rotate finned projectiles；Magnus effect；numerical simulation

2014-05-14；

：2014-08-11。

國(guó)家自然科學(xué)基金(11372136)；國(guó)家部委資助項(xiàng)目。

趙博博(1989—)，男，博士生，研究方向?yàn)閺椉傮w設(shè)計(jì)及氣動(dòng)優(yōu)化。E-mail:zhao-bo-bo@163.com

V211

A

1006-2793(2015)04-0465-07

10.7673/j.issn.1006-2793.2015.04.003