鐘 陽,王良明,李 巖,楊志偉

(南京理工大學(xué) 能源與動(dòng)力工程學(xué)院,江蘇 南京 210094)

20世紀(jì)90年代,美國陸軍裝備研發(fā)中心提出了低成本強(qiáng)力彈藥計(jì)劃。二維修正引信具有縱向橫向修正能力,精度高、成本低和對(duì)常規(guī)庫存彈藥兼容性好等優(yōu)點(diǎn),受到各國青睞,成為當(dāng)前的研究熱點(diǎn),其中采用固定鴨舵式的精確制導(dǎo)組件(precise guidance kit,PGK)最為經(jīng)典和成功。

氣動(dòng)特性技術(shù)是固定鴨舵二維修正彈的關(guān)鍵技術(shù)之一。常規(guī)彈藥改裝修正組件后,其氣動(dòng)特性必然發(fā)生改變,飛行過程中存在失穩(wěn)現(xiàn)象,氣動(dòng)特性是首先要考慮的問題[1]。吳萍采用實(shí)驗(yàn)手段對(duì)彈丸的氣動(dòng)特性展開研究[2];徐輝雯等人仿真分析了二維修正彈修正組件反旋與不旋時(shí)氣動(dòng)特性的差異[3];程杰通過對(duì)有攻角和有舵偏角2種狀態(tài)的分析,建立了基于2種狀態(tài)疊加的氣動(dòng)力工程模型[4];美國陸軍研究所的Sahu和Silton等人長期從事修正彈的氣動(dòng)特性數(shù)值計(jì)算,其數(shù)值計(jì)算結(jié)果與大量實(shí)驗(yàn)結(jié)果吻合較好,這也說明了CFD技術(shù)在獲取彈丸氣動(dòng)力方面的可靠性。隨著對(duì)修正彈氣動(dòng)特性研究的深入,其越來越細(xì)化,不僅對(duì)整彈氣動(dòng)參數(shù)進(jìn)行計(jì)算,近來,還針對(duì)關(guān)鍵性氣動(dòng)部件鴨舵展開氣動(dòng)特性研究。Sahu等人[5]采用數(shù)值仿真計(jì)算了鴨舵產(chǎn)生的滾轉(zhuǎn)力矩特性并將其應(yīng)用于彈丸的滾轉(zhuǎn)控制中;Silton等人[6-7]通過靜態(tài)掃描模擬技術(shù),數(shù)值計(jì)算了×型布局下,各片鴨舵產(chǎn)生的俯仰力矩隨攻角的變化歷程;彭程[8]采用CFD動(dòng)網(wǎng)格技術(shù)研究了帶鴨舵彈丸錐形運(yùn)動(dòng)和自轉(zhuǎn)運(yùn)動(dòng)耦合作用下的空氣動(dòng)力學(xué)特性,分析了鴨舵誘導(dǎo)渦系結(jié)構(gòu)的復(fù)雜發(fā)展情況以及周期性且有規(guī)律的氣動(dòng)力。含鴨舵的控制組件作為區(qū)別一般彈丸和二維修正彈的重要標(biāo)志,也是二維修正彈最關(guān)鍵的氣動(dòng)部件,差動(dòng)舵產(chǎn)生的軸向力矩決定著控制組件的轉(zhuǎn)動(dòng)情況,控制舵產(chǎn)生的法向力矩決定著彈丸飛行姿態(tài),而以上2種力矩均可由鴨舵法向力近似計(jì)算獲取。有必要建立一個(gè)較為精確的、滿足控制組件轉(zhuǎn)動(dòng)或者不轉(zhuǎn)動(dòng)的適合旋轉(zhuǎn)穩(wěn)定二維修正彈的鴨舵法向力模型。基于此,本文針對(duì)之前工作的不足,采用多元泰勒展開理論,并對(duì)參數(shù)進(jìn)行分析,建立關(guān)于彈丸攻角和舵偏角的二元鴨舵法向力計(jì)算模型。通過建立鴨舵坐標(biāo)系,并考慮彈丸運(yùn)動(dòng)和迎風(fēng)區(qū)背風(fēng)區(qū)的影響,將模型擴(kuò)展為適合于動(dòng)態(tài)控制組件下的動(dòng)態(tài)鴨舵法向力計(jì)算模型。為旋轉(zhuǎn)穩(wěn)定二維修正彈控制組件的滾轉(zhuǎn)空氣動(dòng)力矩計(jì)算和全彈氣動(dòng)力計(jì)算提供參考。

1 鴨舵法向力模型

旋轉(zhuǎn)穩(wěn)定二維修正彈是通過調(diào)整控制組件的轉(zhuǎn)動(dòng)情況,由一對(duì)反向差動(dòng)舵和一對(duì)同向控制舵產(chǎn)生期望的俯仰操縱力矩和偏航操縱力矩來控制彈丸姿態(tài),從而達(dá)到彈道修正的目的。鴨舵作為控制組件上的關(guān)鍵氣動(dòng)部件,有必要研究其氣動(dòng)特性。其中,法向力能夠產(chǎn)生使控制組件轉(zhuǎn)動(dòng)的軸向力矩,同時(shí)其產(chǎn)生的法向力矩也是彈丸俯仰操縱力矩和偏航操縱力矩的重要組成部分,因此,本節(jié)對(duì)鴨舵法向力展開研究,以期建立較為準(zhǔn)確的計(jì)算模型。

1.1 鴨舵垂直于攻角面時(shí)的法向力模型

對(duì)于0°滾轉(zhuǎn)角下的法向力模型,前人做過一些研究,文獻(xiàn)[9]將攻角和舵偏角作為合攻角,以合攻角為參數(shù),采用四階多項(xiàng)式來表示彈翼法向力系數(shù),考慮到攻角和舵偏角的耦合作用,模型中引入了2個(gè)相關(guān)系數(shù);文獻(xiàn)[4]將鴨舵法向力系數(shù)看作是0°攻角和0°舵偏角2種情況下的疊加,并基于該思想建立了工程氣動(dòng)模型。本節(jié)基于多元泰勒展開理論,考慮攻角和舵偏角的因素,并通過分析,對(duì)相關(guān)系數(shù)進(jìn)行取舍,建立鴨舵垂直于攻角面時(shí)的法向力系數(shù)模型。

氣動(dòng)力計(jì)算采用的坐標(biāo)系OAxAyAzA如圖1所示,原點(diǎn)OA在彈丸質(zhì)心處,OAxA軸與彈軸重合,沿彈尾指向彈頭為正;OAyA軸在攻角平面內(nèi)向上為正;OAzA軸滿足右手法則。

圖1 氣動(dòng)坐標(biāo)系示意圖

圖2為控制組件從彈尾指向彈頭方向的視圖。

圖2 控制組件后視圖及鴨舵坐標(biāo)系示意圖

鴨舵坐標(biāo)系是與鴨舵固連的坐標(biāo)系,定義鴨舵i(i=1,2,3,4)坐標(biāo)系Oxiyizi,Oxi軸與彈軸重合,由彈尾指向彈頭為正;Ozi軸過鴨舵形心并與彈軸相交,由翼根指向翼稍為正,交點(diǎn)O為原點(diǎn);Oyi軸滿足右手法則。Oyi軸繞Oxi軸正向轉(zhuǎn)到鴨舵i坐標(biāo)系中Ozi軸時(shí),轉(zhuǎn)過的角度為鴨舵i的滾轉(zhuǎn)角γi(i=1,2,3,4),控制組件的滾轉(zhuǎn)角與鴨舵1相等。

在給定空氣來流的馬赫數(shù)和雷諾數(shù)條件下,忽略舵片之間的相互干擾,靜態(tài)鴨舵法向力主要與彈丸攻角、舵偏角以及控制組件滾轉(zhuǎn)角有關(guān)。本節(jié)討論當(dāng)控制組件不滾轉(zhuǎn)且滾轉(zhuǎn)角為0°,即一對(duì)控制舵垂直于總攻角面時(shí),鴨舵對(duì)彈丸產(chǎn)生的法向力。此時(shí),操縱鴨舵的法向力可視為舵偏角和彈丸攻角的二元函數(shù),則鴨舵法向力系數(shù)為

CN=f(α,δ)

(1)

式中:α為彈丸攻角,δ為舵偏角。假設(shè)式(1)在點(diǎn)(0,0)的某鄰域內(nèi)具有n+1階連續(xù)偏導(dǎo)數(shù),采用泰勒展開:

(2)

式中:0<θ<1。考慮到非線性的影響,將式(2)保留到三階,且有彈丸攻角和舵偏角都為0°時(shí),鴨舵不產(chǎn)生法向力,即f(0,0)=0。重寫式(2),有:

(3)

式中:下標(biāo)α和δ均表示函數(shù)f對(duì)其求偏導(dǎo)?？紤]到:①當(dāng)彈丸攻角為0°時(shí),大小相等且方向相反的舵偏角產(chǎn)生的法向力也必然大小相等且方向相反,即式(3)有f(0,δ)=-f(0,-δ),得fδδ(0,0)=0;②同理,當(dāng)舵偏角為0°時(shí),大小相等且方向相反的彈丸攻角產(chǎn)生的法向力也必然大小相等且方向相反,即式(3)也有f(α,0)=-f(-α,0),得fαα(0,0)=0;③當(dāng)彈體攻角和舵偏角在保持大小不變的情況下同時(shí)反向,則法向力大小不變,方向反向,即式(3)有f(α,δ)=-f(-α,-δ),再綜合①和②可得fαδ(0,0)=0。

通過以上分析,并令:

(4)

將式(4)代入式(3)有:

(5)

式(5)為控制組件滾轉(zhuǎn)角為0°時(shí),控制鴨舵法向力系數(shù)關(guān)于彈丸總攻角和舵偏角的非線性模型。該模型中包含有攻角和舵偏角的線性、非線性項(xiàng)及交叉干擾項(xiàng)。各系數(shù)的獲取方法采用多元線性回歸,將模型(5)寫成一般多元線性回歸形式:

(6)

對(duì)于固定鴨舵二維修正彈,舵偏角可視為常量,令δ=δi,式(5)可寫為

CN(α,δ)|γF=0,δ=δi=k0+k1α+k2α2+k3α3

(7)

1.2 動(dòng)態(tài)任意滾轉(zhuǎn)角下的鴨舵法向力模型

上節(jié)建立了當(dāng)鴨舵垂直于攻角面時(shí),其關(guān)于攻角和舵偏角的法向力系數(shù)模型。當(dāng)鴨舵處于任意滾轉(zhuǎn)角時(shí),其法向力系數(shù)的獲取方法通常是先得到有效攻角,然后把問題轉(zhuǎn)化到0°滾轉(zhuǎn)角下解決。本文借鑒該思想,并考慮彈丸運(yùn)動(dòng)和迎風(fēng)區(qū)背風(fēng)區(qū)等因素,建立動(dòng)態(tài)任意滾轉(zhuǎn)角下的鴨舵法向力模型。首先對(duì)鴨舵坐標(biāo)系下的舵偏角和有效攻角進(jìn)行說明。

圖3為某鴨舵在其鴨舵坐標(biāo)系中的姿態(tài),當(dāng)該鴨舵沿Oxi軸正向運(yùn)動(dòng)時(shí),產(chǎn)生的空氣動(dòng)力在Oyi軸分量為負(fù),此時(shí)舵偏角δi為負(fù),反之則δi為正。

圖4中vri為鴨舵i壓心在鴨舵系下相對(duì)于空氣的速度矢量,αe,i為鴨舵i的有效攻角。法向力的大小和方向主要取決于空氣與鴨舵壓心的相對(duì)速度矢量。

圖3 鴨舵系下舵偏角示意圖

圖4 鴨舵i坐標(biāo)系及相關(guān)角度示意圖

圖5中rGC為彈丸質(zhì)心位置G到鴨舵壓心位置C的矢徑,rGCA為矢徑rGC在彈軸上的投影,rGCN為鴨舵壓心C到彈軸的距離。鴨舵壓心一般在鴨舵形心附近。當(dāng)鴨舵i滾轉(zhuǎn)角為γi時(shí),質(zhì)心到鴨舵i壓心矢徑為rGCi=(rGCAirGCNicosγirGCNisinγi)T。

圖5 彈丸質(zhì)心與鴨舵壓心位置關(guān)系示意圖

鴨舵i壓心速度在氣動(dòng)坐標(biāo)系OAxAyAzA下有矢量關(guān)系式:

(8)

OAxAyAzA坐標(biāo)系向鴨舵i坐標(biāo)系的轉(zhuǎn)換矩陣為

(9)

則鴨舵壓心速度矢量在鴨舵系中的映射為

(10)

式中:vri=(vxivyivzi)T。鴨舵i有效攻角αe,i滿足關(guān)系:

αe,i(γi,ωFx,ωy,ωz)=arctan(-vyi/vxi)

(11)

在有攻角的情況下,空氣繞流受到彈體的影響,存在迎風(fēng)壓縮區(qū)和背風(fēng)膨脹分離區(qū),不同區(qū)域?qū)喍嫣峁┓ㄏ蛄Φ男什煌?通常迎風(fēng)區(qū)效率高,背風(fēng)區(qū)效率低。文獻(xiàn)[9]數(shù)值計(jì)算了翼片在滾轉(zhuǎn)角為0°時(shí),即強(qiáng)度最高背風(fēng)區(qū)時(shí)翼片法向力效率損失隨攻角的變化曲線;滾轉(zhuǎn)角為180°時(shí),即強(qiáng)度最高迎風(fēng)區(qū)利用牛頓流理論構(gòu)造了壓縮增效函數(shù)。二維修正彈的鴨舵在滾轉(zhuǎn)過程中會(huì)周期性處在迎風(fēng)區(qū)和背風(fēng)區(qū),為了使得法向力模型能夠反映出迎風(fēng)區(qū)和背風(fēng)區(qū)的周期性影響,需要引入具有周期性的鴨舵效率模型。在鴨舵i坐標(biāo)系中,αe,i為有效攻角,則βe,i=arctan(vzi/vxi)可相應(yīng)稱之為鴨舵i的有效偏角。隨著控制組件周期性自轉(zhuǎn),有效偏角呈周期性變化,當(dāng)有效偏角為正時(shí)鴨舵處于迎風(fēng)區(qū),否則鴨舵處于背風(fēng)區(qū)。則可利用有效偏角來構(gòu)造鴨舵效率函數(shù):

(12)

式中:ke為鴨舵效率函數(shù),kde和kie分別為背風(fēng)減效因子和迎風(fēng)增效因子,其大小與βe相關(guān)。不妨將kde在βe=0附近進(jìn)行二階泰勒展開:

(13)

當(dāng)βe=0時(shí),易知鴨舵不處于迎風(fēng)區(qū)和背風(fēng)區(qū),鴨舵效率不提高也不減小,因此kde(0)=0。同理對(duì)減效因子kie進(jìn)行類似處理,于是式(12)可寫為

(14)

得到修正后的鴨舵坐標(biāo)系下的法向力系數(shù)模型:

CNM=ke(βe)CN(αe,δ)

(15)

通過轉(zhuǎn)換矩陣L-1(γ-π/2)和鴨舵系下的法向合力系數(shù)CNM=(0CNM0)T便可得到計(jì)算氣動(dòng)力坐標(biāo)系OAxAyAzA下的動(dòng)態(tài)鴨舵法向力系數(shù)模型:

(16)

2 數(shù)值仿真驗(yàn)證與分析

數(shù)值仿真采用有限體積法離散Navier-Stokes方程組,對(duì)流通量計(jì)算格式為SLAU2,粘性通量采用中心差分格式;湍流模型采用標(biāo)準(zhǔn)可壓Spalart-Allmaras(S-A)模型。幾何模型如圖1所示,鴨舵1和鴨舵3為反向差動(dòng)舵,舵偏角為4°,鴨舵2和鴨舵4為同向控制舵,舵偏角為8°。采用結(jié)構(gòu)網(wǎng)格對(duì)流場進(jìn)行劃分;遠(yuǎn)場邊界和近壁層網(wǎng)格設(shè)置參考文獻(xiàn)[10];二維修正彈的運(yùn)動(dòng)采用人工滑移邊界模擬。網(wǎng)格細(xì)節(jié)見圖6和圖7,網(wǎng)格總數(shù)約350萬。

圖6 流場網(wǎng)格剖面

圖7 滑移網(wǎng)格結(jié)構(gòu)

本文通過對(duì)比數(shù)值計(jì)算方法和飛行試驗(yàn)方法獲取的控制組件轉(zhuǎn)速來間接說明本文所采用的數(shù)值方法的可信度。本文研究對(duì)象在控制組件自由滾轉(zhuǎn)的條件下進(jìn)行了飛行試驗(yàn)。通過飛行過程中獲取的彈身轉(zhuǎn)速和控制組件相對(duì)彈身的轉(zhuǎn)速,處理出控制組件的絕對(duì)轉(zhuǎn)速ωFx。

圖8中,通過對(duì)比數(shù)值方法[11]處理出的控制組件平衡轉(zhuǎn)速和飛行試驗(yàn)獲得的轉(zhuǎn)速,其變化規(guī)律大致相同,各數(shù)值計(jì)算點(diǎn)和試驗(yàn)數(shù)據(jù)吻合較好,誤差小于15%?？刂平M件的滾轉(zhuǎn)空氣動(dòng)力矩主要是由鴨舵法向力產(chǎn)生的,說明通過本文數(shù)值方法得到的鴨舵法向力是合理的,計(jì)算結(jié)果具有參考價(jià)值。本文采用該數(shù)值方法對(duì)鴨舵法向力模型的合理有效性加以驗(yàn)證分析。

驗(yàn)證分為兩部分:第一部分,采用定常CFD方法仿真驗(yàn)證模型(5),能夠預(yù)測不同攻角和舵偏角下的鴨舵法向力;第二部分,采用非定常CFD方法仿真驗(yàn)證控制組件滾轉(zhuǎn)條件下動(dòng)態(tài)鴨舵法向力模型(16),能夠預(yù)測出每一片鴨舵法向力規(guī)律。

模型(5)中參數(shù)采用定常流場計(jì)算結(jié)果進(jìn)行擬合。通過網(wǎng)格旋轉(zhuǎn),舵偏角可以為±4°、±8°。定常流場工況:馬赫數(shù)2.0,壓力101 325 Pa,溫度288.15 K,密度1.225 kg/m3;攻角和舵偏角組合見圖9中數(shù)值計(jì)算點(diǎn)。為了擬合模型(5)中的系數(shù),從數(shù)值計(jì)算的單片鴨舵法向力系數(shù)中選取19個(gè)不同攻角和舵偏角組合,如表1所示。

圖8 控制組件轉(zhuǎn)速隨馬赫數(shù)的變化關(guān)系

圖9 不同舵偏角下法向力系數(shù)隨攻角的變化關(guān)系

編號(hào)δ/(°)α/(°)CN18-6-1.838 7×10-328-42.405 2×10-338-26.738 0×10-34801.113 8×10-25821.558 8×10-26842.004 5×10-27852.225 5×10-28862.442 5×10-29882.863 3×10-2108103.246 0×10-211812-1.149 4×10-2124-8-7.404 9×10-3134-6-3.153 4×10-3144-41.195 2×10-315405.584 3×10-316429.969 4×10-317441.437 6×10-218461.873 9×10-219482.303 3×10-2

利用表1中的參數(shù),通過式(6)獲得多元回歸方程組,進(jìn)而通過最小二乘擬合出模型(5)中各參數(shù)。在|α|≤12°,|δ|≤8°范圍內(nèi),鴨舵法向力系數(shù)擬合結(jié)果相對(duì)數(shù)值計(jì)算值誤差小于2.05%。擬合參數(shù)結(jié)果如表2所示。

表2 模型(5)參數(shù)擬合結(jié)果

從圖9中可以看出,采用4°和8°舵偏角擬合出的系數(shù),同樣適用于-4°和-8°,即在定馬赫數(shù)和雷諾數(shù)下,無需改變表2中的參數(shù),模型(5)計(jì)算出的鴨舵法向力系數(shù)與數(shù)值計(jì)算出的不同攻角、不同舵偏角下的結(jié)果吻合得很好,反映出本文針對(duì)彈丸攻角和舵偏角所建立的二元鴨舵法向力模型是合理有效的。由圖9中“δ=8°”曲線可以看出,當(dāng)攻角約10°時(shí)開始出現(xiàn)非線性現(xiàn)象。將舵偏角為8°時(shí)的鴨舵附近流場作剖面,對(duì)流場進(jìn)行觀察。由圖10可知,當(dāng)攻角為10°時(shí),空氣在舵面發(fā)生了分離,與仿真曲線表現(xiàn)出的非線性相一致。

圖10 舵偏角8°時(shí)不同攻角下空氣流動(dòng)情況

模型(16)中含有鴨舵效率模型(14),該模型中的相關(guān)參數(shù)采用定常流場計(jì)算結(jié)果擬合,工況為:鴨舵滾轉(zhuǎn)角0°和180°,舵偏角8°,彈丸攻角0°、2°、4°、6°、8°,其余條件與之前相同。當(dāng)滾轉(zhuǎn)角為180°時(shí),有效偏角βe與攻角相等;當(dāng)滾轉(zhuǎn)角為0°時(shí),有效偏角βe為攻角相反數(shù)。

為了得到控制組件滾轉(zhuǎn)時(shí)4個(gè)鴨舵對(duì)彈丸產(chǎn)生的法向力系數(shù),采用滑移網(wǎng)格來模擬轉(zhuǎn)動(dòng),非定常流場工況為:控制舵偏角大小為8°,差動(dòng)舵偏角大小為4°;馬赫數(shù)Ma=2.0;控制組件轉(zhuǎn)速ωFx=-100 rad/s;攻角α=2°,5°。雙時(shí)間步收斂精度為0.01,內(nèi)層最大迭代100步,每個(gè)時(shí)間步轉(zhuǎn)1°。

圖11 迎風(fēng)區(qū)背風(fēng)區(qū)鴨舵效率隨有效偏角關(guān)系曲線

圖12中,從原點(diǎn)指向曲線上任意一點(diǎn)所成的矢量為鴨舵法向合力系數(shù)矢量,當(dāng)控制組件以角速度ωFx=-100 rad/s滾轉(zhuǎn)時(shí),該矢量末端軌跡在面OAyAzA上形成周期性曲線,可稱為法向力矢端軌跡線。圖12中的軌跡線分別為CFD數(shù)值計(jì)算結(jié)果和模型(16)采用表2參數(shù)所計(jì)算出的結(jié)果。從圖12(a)和圖12(b)中CFD得到的鴨舵法向力矢端軌跡線可知,由于初始流場不精確,非定常數(shù)值計(jì)算結(jié)果初期不合理,但這種初始效應(yīng)隨著時(shí)間推進(jìn)逐漸消失,最后形成周期性的曲線;從圖12(a)和圖12(b)還可以看出,鴨舵1和鴨舵3的軌跡線幾乎一致,因?yàn)轼喍?和鴨舵3只有相位上的差別,它們轉(zhuǎn)到同一位置時(shí)的法向力一致;由于轉(zhuǎn)動(dòng),鴨舵2在整個(gè)過程中的有效攻角是減小的,鴨舵4反之,從圖12(a)可見鴨舵4的法向力系數(shù)曲線比鴨舵2整體大一圈;鴨舵2先進(jìn)入背風(fēng)區(qū)并產(chǎn)生OAzA軸負(fù)向的法向力分量,再進(jìn)入迎風(fēng)區(qū)并產(chǎn)生OAzA軸正向的法向力分量,因此正向分量絕對(duì)值比負(fù)向大,形成了非對(duì)稱形狀,鴨舵4則與之相反,這一點(diǎn)從圖12(b)中可看出。通過上述分析可知,非定常數(shù)值計(jì)算結(jié)果反映出的特性是合理的。

另外,由圖12(a)和圖12(b)攻角與舵偏角的關(guān)系大致來看,當(dāng)攻角大于舵偏角時(shí),軌跡線為雙環(huán)形,如圖12(b)鴨舵1。此時(shí),整個(gè)法向力矢端軌跡線幾乎都在OAyA軸正向,即無論控制組件如何轉(zhuǎn)動(dòng),鴨舵幾乎都不能產(chǎn)生OAyA軸負(fù)向的力;當(dāng)攻角與舵偏角相當(dāng)時(shí),軌跡為蘋果型,如圖12(a)鴨舵1和圖12(b)鴨舵2。此時(shí),鴨舵不僅能夠產(chǎn)生正向的力,也能夠產(chǎn)生較小負(fù)向的力;當(dāng)舵偏角大于攻角時(shí),軌跡為單環(huán)型,如圖12(a)鴨舵2。此時(shí),鴨舵能向各個(gè)方向提供較大的法向力。為了鴨舵能提供足夠的法向力,一種是增大鴨舵的氣動(dòng)面積,另一種是增加舵偏角。基于本文的分析,更傾向于在不發(fā)生分離的情況下盡量選擇大舵偏角,這樣更有可能為彈丸提供各個(gè)方向較大的法向力,提高鴨舵利用率。

圖12 不同攻角下法向力系數(shù)矢端軌跡線

對(duì)比圖12中數(shù)值和模型計(jì)算結(jié)果可以看出,動(dòng)態(tài)鴨舵法向力模型(16)在來流條件一定的情況下,僅采用一組系數(shù)便能預(yù)測出控制組件滾轉(zhuǎn)過程中不同攻角下每片鴨舵上的法向力系數(shù),預(yù)測結(jié)果與數(shù)值計(jì)算結(jié)果基本吻合,說明本文提出的模型是合理有效的。

3 結(jié)束語

本文針對(duì)旋轉(zhuǎn)穩(wěn)定二維彈道修正彈,基于多元泰勒展開理論,并對(duì)系數(shù)進(jìn)行分析取舍,得到鴨舵垂直于攻角面時(shí),其關(guān)于攻角和舵偏角的二元法向力模型。考慮到在實(shí)際飛行中,控制組件會(huì)發(fā)生滾轉(zhuǎn),且鴨舵會(huì)交替出現(xiàn)在迎風(fēng)區(qū)和背風(fēng)區(qū),通過建立鴨舵坐標(biāo)系和加入鴨舵效率模型,進(jìn)一步將該模型擴(kuò)展為動(dòng)態(tài)鴨舵法向力模型。最后通過數(shù)值仿真驗(yàn)證了該模型的合理性,為旋轉(zhuǎn)穩(wěn)定二維修正彈的設(shè)計(jì)提供參考。