• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    結(jié)構(gòu)誤差對(duì)旋轉(zhuǎn)穩(wěn)定彈丸氣動(dòng)特性影響的數(shù)值模擬

    2014-04-30 07:25:08陳東陽(yáng)LaithAbbas芮筱亭王國(guó)平
    關(guān)鍵詞:彈體湍流氣動(dòng)

    陳東陽(yáng),Laith K.Abbas,芮筱亭,王國(guó)平

    (南京理工大學(xué) 發(fā)射動(dòng)力學(xué)研究所,江蘇 南京 210094)

    結(jié)構(gòu)誤差對(duì)旋轉(zhuǎn)穩(wěn)定彈丸氣動(dòng)特性影響的數(shù)值模擬

    陳東陽(yáng),Laith K.Abbas,芮筱亭,王國(guó)平

    (南京理工大學(xué) 發(fā)射動(dòng)力學(xué)研究所,江蘇 南京 210094)

    基于剪切應(yīng)力傳輸(SST)k-ω湍流模型,對(duì)考慮結(jié)構(gòu)誤差情況下的彈箭模型進(jìn)行數(shù)值模擬計(jì)算。首先采用CFD和工程經(jīng)驗(yàn)公式相結(jié)合的方法,得到了不同馬赫數(shù)、不同轉(zhuǎn)速情況下無(wú)結(jié)構(gòu)誤差模型的多種氣動(dòng)特性參數(shù)。通過(guò)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比,阻力、法向力、俯仰力矩、壓心位置誤差在10%以內(nèi),驗(yàn)證了該方法的可行性和準(zhǔn)確性。然后,建立了考慮結(jié)構(gòu)誤差即質(zhì)量分布不對(duì)稱、彈體不同軸和無(wú)結(jié)構(gòu)誤差模型,并分別進(jìn)行了氣動(dòng)特性模擬計(jì)算。結(jié)果表明,質(zhì)量偏心對(duì)滾轉(zhuǎn)阻尼力矩系數(shù)、馬格努斯力矩系數(shù)影響很大;彈體不同軸對(duì)法向力系數(shù)導(dǎo)數(shù)、俯仰力矩系數(shù)導(dǎo)數(shù)、馬格努斯力矩系數(shù)導(dǎo)數(shù)和壓心位置都產(chǎn)生很大影響。結(jié)構(gòu)誤差使模型的多種氣動(dòng)特性參數(shù)產(chǎn)生很大變化,將影響其飛行彈道和穩(wěn)定性。

    結(jié)構(gòu)誤差;湍流模型;計(jì)算流體力學(xué);旋轉(zhuǎn)

    0 引 言

    在彈箭的生產(chǎn)制造過(guò)程中,由于設(shè)備、工藝以及技術(shù)水平等諸多因素的影響,彈箭結(jié)構(gòu)總會(huì)和理論設(shè)計(jì)的有一定的偏差。對(duì)于彈箭,主要的結(jié)構(gòu)誤差有彈體質(zhì)量分布不均勻?qū)е碌馁|(zhì)量偏心、由于安裝誤差導(dǎo)致的彈體不同軸、翼面安裝角誤差等。這些誤差對(duì)高速飛行彈箭的氣動(dòng)特性和飛行穩(wěn)定性都會(huì)帶來(lái)影響。因此,對(duì)彈箭氣動(dòng)參數(shù)的準(zhǔn)確獲取和對(duì)結(jié)構(gòu)誤差可能帶來(lái)的氣動(dòng)特性影響進(jìn)行精確仿真分析,十分必要。隨著計(jì)算機(jī)的發(fā)展和數(shù)值模擬軟件的計(jì)算能力的提高與完善,采用數(shù)值計(jì)算方法獲取飛行器氣動(dòng)數(shù)據(jù)成為可能。相對(duì)于風(fēng)洞實(shí)驗(yàn)和工程估算,CFD仿真有其特有的優(yōu)勢(shì)。CFD仿真可以很好的預(yù)測(cè)復(fù)雜幾何體的氣動(dòng)特性參數(shù)和流體流動(dòng)現(xiàn)象[1-5],可以模擬風(fēng)洞實(shí)驗(yàn)無(wú)法模擬的條件,并且可以減少實(shí)驗(yàn)次數(shù),大大節(jié)省設(shè)計(jì)成本。表1[6]給出了三種方法預(yù)測(cè)氣動(dòng)力的誤差統(tǒng)計(jì)結(jié)果。其中飛行試驗(yàn)預(yù)測(cè)氣動(dòng)力的結(jié)果要好于風(fēng)洞實(shí)驗(yàn)和工程程序。文獻(xiàn)[7]基于CFD技術(shù),比較了RNGk-ε、k-εRealizable、DES等湍流模型的計(jì)算結(jié)果,證明了CFD程序有很高的計(jì)算精度,甚至有時(shí)可以代替飛行試驗(yàn),為飛行器彈道計(jì)算、飛行控制提供準(zhǔn)確氣動(dòng)特性數(shù)據(jù)。

    表1 空氣動(dòng)力系數(shù)的預(yù)估誤差(%)[6]Table 1 Expected aerodynamic coefficient percent errors

    本文采用基于SST湍流模型的高精度CFD技術(shù),和高質(zhì)量的網(wǎng)格生成方法對(duì)M910旋轉(zhuǎn)彈丸模型進(jìn)行數(shù)值計(jì)算,獲得模型的各種氣動(dòng)參數(shù)并與實(shí)驗(yàn)數(shù)據(jù)比較,來(lái)驗(yàn)證該數(shù)值方法的可行性與準(zhǔn)確性。然后通過(guò)對(duì)帶有結(jié)構(gòu)誤差的兩種模型采用同樣的數(shù)值方法進(jìn)行數(shù)值模擬計(jì)算,精確獲得其各種氣動(dòng)參數(shù),然后分析結(jié)構(gòu)誤差對(duì)旋轉(zhuǎn)彈丸模型氣動(dòng)特性,飛行彈道和穩(wěn)定性可能帶來(lái)的影響。網(wǎng)格劃分軟件為ICEMCFD,流體解算器為Fluent。

    1 幾何模型

    選用M910旋轉(zhuǎn)穩(wěn)定彈丸為假設(shè)的模型,幾何參數(shù)、飛行速度、旋轉(zhuǎn)速度等參數(shù)詳見文獻(xiàn)[7]。本文建立了三種幾何模型即無(wú)結(jié)構(gòu)誤差模型、質(zhì)量偏心模型、不同軸模型。如圖1~圖3所示。假設(shè)質(zhì)心相對(duì)無(wú)結(jié)構(gòu)誤差模型的質(zhì)心向右平移了1mm,如圖2所示。如圖3所示,彈體軸線相對(duì)無(wú)結(jié)構(gòu)誤差模型做了不同軸假設(shè),使彈頭處向下偏轉(zhuǎn)了2°。

    圖1 無(wú)結(jié)構(gòu)誤差的模型Fig.1 Model without structure errors

    圖2 有質(zhì)量偏心誤差的模型Fig.2 Model with mass asymmetry errors

    圖3 有同軸誤差的模型Fig.3 Model with coaxial error

    2 計(jì)算方法

    2.1 計(jì)算流體力學(xué)控制方程和湍流模型

    控制方程通用形式[8]:

    式中:φ為通用變量,可代表速度、溫度等求解變量;Γ為廣義擴(kuò)散系數(shù);S為廣義源項(xiàng)。

    式(1)中各項(xiàng)依次為瞬態(tài)項(xiàng)、對(duì)流項(xiàng)、擴(kuò)散項(xiàng)和源項(xiàng)。表2給出了三個(gè)符號(hào)與各特定方程的對(duì)應(yīng)關(guān)系。

    表2 通用控制方程中各符號(hào)的具體形式Table 2 Specific expressions of three parameters in general control equation

    對(duì)于湍流,如果直接求解高度非線性的三維瞬態(tài)N-S方程,需要采用對(duì)計(jì)算機(jī)內(nèi)存和速度要求很高的直接模擬方法,目前還不可能在實(shí)際工程中采用此方法。工程上對(duì)瞬態(tài)的N-S方程做時(shí)間平均處理,得到雷諾平均的NS方程(RANS):

    其中:

    本文選用渦粘模型(EVM)中的SSTk-ω湍流模型[9],它結(jié)合了在邊界層內(nèi)能很好模擬低雷諾數(shù)流動(dòng)的標(biāo)準(zhǔn)k-ω湍流模型,在邊界層外能很好模擬完全湍流流動(dòng)的k-ε湍流模型的優(yōu)點(diǎn)。其間通過(guò)一個(gè)混合函數(shù)來(lái)過(guò)渡,屬于兩方程渦粘性模式。k和ω的輸運(yùn)方程如下所示:

    混合函數(shù)

    其中渦粘系數(shù)定義為:

    其中Ω是渦量的絕對(duì)值。

    混合函數(shù)F2定義為:

    SST湍流模型中常數(shù)通過(guò)式(9)混合:

    其中集合(φ1)代表標(biāo)準(zhǔn)的k-ω湍流模型中的常數(shù),集合(φ2)代表標(biāo)準(zhǔn)的k-ε湍流模型中的常數(shù)。

    2.2 邊界條件與數(shù)值格式

    遠(yuǎn)場(chǎng)邊界按黎曼不變化無(wú)反射條件來(lái)處理,在物面上采用無(wú)滑移條件來(lái)處理,對(duì)計(jì)算模型設(shè)置壁面旋轉(zhuǎn)。采用基于有限體積法的AUSM格式進(jìn)行空間離散,AUSM兼有Roe格式的間斷高分辨率和Van Leer格式的計(jì)算效率高的優(yōu)點(diǎn)。在對(duì)流項(xiàng),湍流粘性項(xiàng)采用二階迎風(fēng)格式。同時(shí),采用隱式時(shí)間離散格式,對(duì)于超聲速粘性流動(dòng),物體近壁面處的流場(chǎng)會(huì)產(chǎn)生急劇的變化,因此在近壁面處和激波處需要加密網(wǎng)格。

    2.3 計(jì)算網(wǎng)格

    流體與壁面間相互作用,很多因變量具有較大的梯度,而且粘度對(duì)傳輸過(guò)程有很大的影響。k-ω模型(包括SST模型)采用的是自動(dòng)壁面函數(shù),即可以在近壁處,將壁面函數(shù)自動(dòng)調(diào)整為低雷諾數(shù)壁面方程。一般k-ω模型要求近壁處有y+<5[10]。為了準(zhǔn)確模擬滾轉(zhuǎn)阻尼等系數(shù)。在邊界層內(nèi)對(duì)網(wǎng)格進(jìn)行加密,保證y+≤0.5,以保證邊界層內(nèi)有十層以上的網(wǎng)格。本文采用多Block生成流場(chǎng)拓?fù)浣Y(jié)構(gòu)、外O Block生成彈體邊界層的方法,生成了高質(zhì)量的結(jié)構(gòu)化網(wǎng)格。

    為了減小計(jì)算量,將每套網(wǎng)格數(shù)量都?jí)嚎s在50萬(wàn)左右,僅為文獻(xiàn)[7]網(wǎng)格數(shù)量的1/3~1/2。但通過(guò)提高網(wǎng)格質(zhì)量,保證了計(jì)算精度。由于來(lái)流速度范圍較大,流場(chǎng)均分別采用兩套計(jì)算網(wǎng)格。部分流場(chǎng)拓?fù)浣Y(jié)構(gòu)如圖4、圖5所示,部分流場(chǎng)計(jì)算網(wǎng)格如圖6、圖7所示。

    圖4 不同軸模型流場(chǎng)拓?fù)鋱D(Ma>1.3)Fig.4 Fluid field topology graph of the model with coaxial error(Ma>1.3)

    圖5 不同軸模型流場(chǎng)拓?fù)鋱D(Ma<1.3)Fig.5 Fluid field topology graph of the model with coaxial error(Ma<1.3)

    圖6 無(wú)結(jié)構(gòu)誤差模型流場(chǎng)計(jì)算網(wǎng)格(Ma>1.3)Fig.6 Flow field grid of the model without structure errors

    圖7 不同軸模型流場(chǎng)網(wǎng)格剖面圖Fig.7 Pr of ile of flow field grid of the model with coaxial error

    3 計(jì)算公式

    CFD求解氣動(dòng)特性參數(shù)的原理是,先對(duì)流場(chǎng)進(jìn)行求解,得到流場(chǎng)每一個(gè)網(wǎng)格上的速度、密度、壓力等參數(shù),然后求出氣動(dòng)系數(shù)。

    動(dòng)穩(wěn)定性導(dǎo)數(shù)簡(jiǎn)稱動(dòng)導(dǎo)數(shù),是飛行器彈道計(jì)算和飛行控制所不可缺少的原始?xì)鈩?dòng)參數(shù)。對(duì)于外形簡(jiǎn)單,小攻角的彈箭可采用CFD與經(jīng)驗(yàn)公式結(jié)合的方法[11-12],基于SST湍流模型,較好地計(jì)算了旋轉(zhuǎn)彈的滾轉(zhuǎn)阻尼等系數(shù)。

    類似于風(fēng)洞自由滾轉(zhuǎn)技術(shù)原理公式[11],假設(shè)彈繞x軸以定角速度ωx做勻角速度旋轉(zhuǎn),則有:

    其中,Cl0為彈箭無(wú)滾轉(zhuǎn)時(shí)的滾轉(zhuǎn)力矩。Cl為氣動(dòng)滾轉(zhuǎn)力矩,以保證彈體做勻角速度旋轉(zhuǎn),由定常的CFD方法求出。Jx表示模型在滾轉(zhuǎn)方向上的慣性矩。由上式推導(dǎo)出滾轉(zhuǎn)阻尼力矩系數(shù):

    同理:俯仰阻尼力矩系數(shù):

    馬格努斯力矩系數(shù):

    對(duì)于本文的光彈體計(jì)算模型,Cl0、Cm0、Cn0的值都為零。

    馬格努斯力系數(shù)導(dǎo)數(shù)公式:

    法向力系數(shù)導(dǎo)數(shù)公式:

    馬格努斯力矩系數(shù)導(dǎo)數(shù)公式:

    俯仰力矩系數(shù)導(dǎo)數(shù)公式:

    4 結(jié)果與分析

    對(duì)三種模型分別進(jìn)行了數(shù)值模擬計(jì)算,馬赫數(shù)分別為0.6、0.9、1.2、1.4、2、2.5、3.5、4.5,攻角為3°。圖8~圖11給出了Ma為2,攻角為3°時(shí),流場(chǎng)速度云圖、壓力云圖、密度云圖和繞彈體的流線圖。從圖8~圖10可以看出頭部和尾部形成了清晰的斜激波在彈肩和彈底部區(qū)域形成了低壓區(qū)。

    由于存在攻角,模型周圍的壓力和密度在彈軸兩側(cè)不對(duì)稱。彈體下方壓力值和密度值高于彈體上方。從圖11可以看到繞彈體周圍速度流線圖,圖中顯示出在彈底部形成許多渦旋。這些激波和渦旋形成的波阻和渦阻給彈體增加了阻力。仿真結(jié)果符合空氣動(dòng)力學(xué)規(guī)律。

    圖8 速度云圖Fig.8 Velocity contour

    圖9 壓力云圖Fig.9 Pressure contour

    圖10 密度云圖Fig.10 Density contour

    圖12~圖14為馬赫數(shù)2、攻角為3°時(shí),x=0.05m處的橫截面上速度矢量云圖,壓力等值線云圖和縱截面上的速度矢量云圖。從圖12、圖13看出,邊界層內(nèi)速度矢量分布和橫截面壓力等值線分布不對(duì)稱,這是有攻角的彈體高速旋轉(zhuǎn)造成的,也是馬格努斯力產(chǎn)生的原因。圖14中也可以看到,貼近壁面處的速度矢量是彈體旋轉(zhuǎn)的方向,然后漸變?yōu)檩S向方向,說(shuō)明模擬精度較高。

    圖11 流線圖Fig.11 Streamline contour

    圖12 x=0.05m處橫截面速度矢量圖Fig.12 Section of velocity vector contour at x=0.05m

    圖13 x=0.05m處橫截面壓力等值線圖Fig.13 Section of pressure contour at x=0.05m

    圖14 縱截面速度矢量圖Fig.14 Longitudinal section of velocity vector contour

    圖15為CFD結(jié)果與PRODAS軟件計(jì)算結(jié)果以及試驗(yàn)結(jié)果的比較,其中模型實(shí)驗(yàn)數(shù)據(jù)來(lái)自文獻(xiàn)[7],從圖15可以看出基于SST湍流模型的CFD數(shù)值模擬能力要比PRODAS軟件模擬結(jié)果更好。

    阻力系數(shù)、法向力系數(shù)導(dǎo)數(shù)、俯仰力矩系數(shù)導(dǎo)數(shù)誤差均在15%以內(nèi)。滾轉(zhuǎn)阻尼系數(shù)誤差在20%以內(nèi),誤差都在工程計(jì)算誤差范圍內(nèi),但是對(duì)于馬格努斯力矩系數(shù)導(dǎo)數(shù),在超聲速部分模擬較好,在亞聲速和跨聲速部分模擬誤差比較大。兩方程渦粘模型不能滿足模擬要求。

    此時(shí),對(duì)于馬格努斯力矩系數(shù)導(dǎo)數(shù)的計(jì)算,湍流模型需要采用大渦模擬或者分離渦模型模擬效果較好[5,7]。這些非定常的計(jì)算方法所需的計(jì)算機(jī)資源過(guò)大、時(shí)間太長(zhǎng),對(duì)計(jì)算機(jī)的硬件性能和內(nèi)存資源都有很大的需求。從圖15看出,阻力系數(shù)在跨聲速處為極值,符合氣動(dòng)特性規(guī)律。由于平衡狀態(tài)下模型滾轉(zhuǎn)阻尼力矩與給定平衡的力矩方向是相反的,所以滾轉(zhuǎn)阻尼系數(shù)值為負(fù)值。壓心位置計(jì)算時(shí)選用彈徑為參考長(zhǎng)度。壓心位置隨馬赫數(shù)的變化趨勢(shì)如圖15(f)所示,和實(shí)驗(yàn)數(shù)據(jù)對(duì)比,誤差較小。

    圖16為無(wú)結(jié)構(gòu)誤差模型、質(zhì)量偏心模型、不同軸模型的氣動(dòng)特性參數(shù)計(jì)算結(jié)果比較圖。從圖16(a)看出,質(zhì)量偏心對(duì)計(jì)算模型的阻力幾乎沒有影響,說(shuō)明阻力的大小只和模型外形有關(guān),與質(zhì)心偏移幾乎無(wú)關(guān)。不同軸模型在亞聲速和跨聲速部分,阻力系數(shù)比無(wú)結(jié)構(gòu)誤差模型大,在超聲速部分阻力系數(shù)比無(wú)結(jié)構(gòu)誤差模型小。最大誤差達(dá)到15%,發(fā)生在亞聲速處。這是由于模型不同軸時(shí),導(dǎo)致模型外形改變,摩阻、渦阻、波阻都有所改變所造成的。

    圖15 氣動(dòng)特性隨馬赫數(shù)的變化Fig.15 Effect of Mach numbers on aerodynamic characteristics

    圖16 三種模型氣動(dòng)特性的比較Fig.16 Comparision of aerodynamic characteristics of three models

    圖16(b)為三種模型的法向力系數(shù)導(dǎo)數(shù)隨馬赫數(shù)變化趨勢(shì)。對(duì)于質(zhì)量偏心模型,導(dǎo)致旋轉(zhuǎn)軸線偏移,在亞聲速和跨聲速部分,計(jì)算結(jié)果和無(wú)結(jié)構(gòu)誤差模型結(jié)果一致;在超聲速時(shí),使法向力系數(shù)導(dǎo)數(shù)相對(duì)無(wú)結(jié)構(gòu)誤差模型的計(jì)算結(jié)果產(chǎn)生一些波動(dòng)。由于不同軸模型的彈頭處攻角相對(duì)理論模型變小,從而導(dǎo)致法向力系數(shù)導(dǎo)數(shù)比無(wú)結(jié)構(gòu)誤差模型計(jì)算結(jié)果小,并且隨著馬赫數(shù)的增大,法向力系數(shù)導(dǎo)數(shù)下降趨勢(shì)變大,最大誤差達(dá)到80%。質(zhì)量偏心導(dǎo)致的旋轉(zhuǎn)軸偏移和不同軸導(dǎo)致的外形改變對(duì)彈的法向力有一定影響,尤其是外形的改變對(duì)法向力影響很大。

    從圖16(c)看出,不同軸模型對(duì)俯仰力矩系數(shù)導(dǎo)數(shù)影響很大,最大誤差達(dá)到91%,這是由于不同軸模型對(duì)法向力系數(shù)影響很大造成的。從圖16(d)可以看出,質(zhì)量偏心導(dǎo)致的旋轉(zhuǎn)軸偏移,從而導(dǎo)致滾轉(zhuǎn)阻尼系數(shù)變化很大,最大誤差已超過(guò)100%。而不同軸導(dǎo)致的外形微小改變對(duì)滾轉(zhuǎn)阻尼系數(shù)幾乎沒有影響。馬格努斯力矩的產(chǎn)生和攻角、旋轉(zhuǎn)有關(guān)。圖16(e)可以看出,不同軸導(dǎo)致的彈頭部攻角改變,和質(zhì)量偏心導(dǎo)致的旋轉(zhuǎn)軸偏移,都會(huì)對(duì)馬格努斯力矩系數(shù)導(dǎo)數(shù)產(chǎn)生了較大影響,尤其是質(zhì)量偏心的情況,對(duì)馬格努斯力矩影響更大,使馬格努斯力矩系數(shù)導(dǎo)數(shù)絕對(duì)值整體變大,最大誤差超過(guò)100%。從圖16(f)可以看出,主要是不同軸導(dǎo)致的彈體外形變化對(duì)壓心位置影響很大,最大誤差為40%。

    旋轉(zhuǎn)穩(wěn)定彈飛行穩(wěn)定性必須滿足陀螺穩(wěn)定條件和動(dòng)態(tài)穩(wěn)定條件[14-15]。其中陀螺穩(wěn)定因子與彈的幾何尺寸、飛行速度、自轉(zhuǎn)速度、極轉(zhuǎn)動(dòng)慣量、赤道轉(zhuǎn)動(dòng)慣量、俯仰力矩系數(shù)導(dǎo)數(shù)、馬格努斯力系數(shù)導(dǎo)數(shù)有關(guān)。動(dòng)態(tài)穩(wěn)定性條件與陀螺穩(wěn)定因子、馬格努斯力矩系數(shù)導(dǎo)數(shù)、俯仰阻尼系數(shù)導(dǎo)數(shù)、阻力系數(shù)、法向力系數(shù)導(dǎo)數(shù)、彈道傾角、重力加速度有關(guān)。從本文計(jì)算結(jié)果發(fā)現(xiàn),三種模型的馬格努斯力矩系數(shù)導(dǎo)數(shù),俯仰力矩系數(shù)導(dǎo)數(shù),法向力系數(shù)導(dǎo)數(shù)差異很大,所以結(jié)構(gòu)誤差對(duì)旋轉(zhuǎn)彈的飛行穩(wěn)定性和飛行彈道將造成很大影響。

    5 結(jié) 論

    通過(guò)計(jì)算結(jié)果比較分析,采用本文生成的低數(shù)量、高質(zhì)量結(jié)構(gòu)化網(wǎng)格,基于SST湍流模型邊界層內(nèi)外的優(yōu)點(diǎn),保證y+≤0.5,并與工程經(jīng)驗(yàn)公式相結(jié)合的方法,可以較快的并保證精度的計(jì)算出各種復(fù)雜的氣動(dòng)特性參數(shù)。兩種結(jié)構(gòu)誤差對(duì)氣動(dòng)特性參數(shù)都有影響,主要表現(xiàn)在質(zhì)量偏心對(duì)滾轉(zhuǎn)阻尼力矩系數(shù),馬格努斯力矩系數(shù)影響很大;彈體不同軸對(duì)法向力系數(shù)導(dǎo)數(shù),俯仰力矩系數(shù)導(dǎo)數(shù),馬格努斯力矩系數(shù)導(dǎo)數(shù)和壓心位置都產(chǎn)生很大影響。由于對(duì)這些氣動(dòng)特性參數(shù)產(chǎn)生的很大影響,將對(duì)彈箭飛行穩(wěn)定性和飛行彈道產(chǎn)生很大影響。

    參考文獻(xiàn):

    [1]WU J,GU Z Q,ZHONG Z H.The application of SST turbulence model in the aerodynamic simulation of the automobile[J].Automotive Engineering,2003,25(4):326-329.(in Chinese)

    吳軍,谷正氣,鐘志華.SST湍流模型在汽車?yán)@流仿真中的應(yīng)用[J].汽車工程,2003,25(4):326-329.

    [2]MENTER F R.Two-equation eddy-viscosity models for engineering applications[J].AIAA Journal,1994,32:1598-1605.

    [3]KLEB W.Computational aeroheating predictions for X-34[J].Journal of Spacecraft and Rockets,1999,36(2):179-188.

    [4]SAHU J,EDGE H L,HEAVEY K R,et al.Computational fluid dynamics modeling of multi-body missile aerodynamic interference[R].U.S.Army Research Laboratory,ARL-TR-1765,Aberdeen Proving Ground,MD,1998.

    [5]DESPIRITO J,PLOSTINS P.CFD prediction of M910 projectile aerodynamics:unsteady wake effect on magnus moment[R].AIAA 2007-6580.

    [6]Arrow Tech Associates.PRODAS version 3 technical manual[M].South Burlington,VT,2002.

    [7]DESPIRITO J,HEAVEY K R.CFD computation of magnus moment and roll damping moment of a spinning projectile[R].AIAA 2004-2713,2004.

    [8]ZHANG S S.Fluid dynamics and application[M].Wuhan:Huazhong University of Science and Technology Press,2011:5-10.(in Chinese)

    張師帥.計(jì)算流體動(dòng)力學(xué)及其應(yīng)用[M].武漢:華中科技大學(xué)出版社,2011:5-10.

    [9]MENTER F R.Zonal two equationk-ωturbulence models for aerodynamic flows[R].AIAA 93-2906,1993.

    [10]ZHANG C C,REN L Q,WANG J.Simulation on flow control for drag reduction of revolution body using bionic dimpled surface[J].ACTA Armamentarii,2009,30(8):1066-1072.(in Chinese)

    張春成,任露泉,王晶,等.旋成體仿生凹坑表面流場(chǎng)控制減阻仿真分析[J].兵工學(xué)報(bào),2009,30(8):1066-1072.

    [11]SUN Z W,CHENG Z Y,BAI J Q.A high efficient method for computing dynamic derivatives of aircraft based on quasi-steady CFD method[J].Flight Dynamics,2010,28(2):28-30.(in Chinese)

    孫智偉,程澤蔭,白俊強(qiáng),等.基于準(zhǔn)定常的飛行器動(dòng)導(dǎo)數(shù)的高效計(jì)算方法[J].飛行力學(xué),2010,28(2):28-30.

    [12]DENG F.Investigations of grid fins aerodynamic shape design and the roll characteristics of wing-body configuration[D].[Ph D.Dissertation].Nanjing:Nanjing University of Science&Technology,2011.(in Chinese)

    鄧帆.柵格翼氣動(dòng)外形設(shè)計(jì)及其翼身組合體滾轉(zhuǎn)特性的研究[D].[博士學(xué)位論文].南京:南京理工大學(xué),2011.

    [13]DESPIRITO J,PLOSTINS P.CFD Prediction of M910 projectile aerodynamics:unsteady wake effect on magnus moment[R].AIAA 2007-6580,2007.

    [14]XU A Y.Study of two-dimension trajectory correction projectile aerodynamic characteristic[D].[Master Thesis].Nanjing:Nanjing University of Science&Technology,2010.(in Chinese)

    許安勇.二維彈道修正彈的氣動(dòng)特性研究[D].[碩士學(xué)位論文].南京:南京理工大學(xué),2010.

    [15]XU M Y.Rocket exterior ballistics[M].Defense Industry Press,1980:224-226.(in Chinese)

    徐明友.火箭外彈道學(xué)[M].國(guó)防工業(yè)出版社,1980:224-226.

    Numerical simulation of a spinning stabilized projectile aerodynamic characteristics effected by structure errors

    CHEN Dongyang,LAITH K.Abbas,RUI Xiaoting,WANG Guoping
    (Institute of Lunch Dynamics,Nanjing University of Science&Technology,Nanjing 210094,China)

    Based on the(Shear Stress Transport)SSTk-ωturbulent model,the effets of vehicle′s model with structure errors are investigated using computational fluid dynamic(CFD)method.Firstly,the aerodynamic coefficients of the model without structure errors are calculated in the conditions of different Mach numbers and spinning speeds,and the results are compared with that of the experimental data.CFD calculations predicted the drag,normal force,pitching moment,and normal force center of pressure very well,and their errors are within 10% of the experimental data which verified the accuracy and feasibility of the CFD output.Then,the models without structure errors and the models with structure errors including projectile′s mass asymmetry and coaxial error are established,and their aerodynamic characteristics are numerically calculated by CFD.The results shown that the model of mass asymmetry have a great influence on its roll damping coefficient and magnus moment coefficient derivative,and the model with coaxial error have a great influence on its normal force coefficient derivative,pitching moment coefficient derivative,magnus moment coefficient derivative and center of pressure location.Projectiles′aerodynamic characteristics coefficients may have a prodigious change due to structure errors,and the flight stability of the model may be effected by these coefficients′prodigious change.The flight trajectory and stability also will be affected accordingly.

    structure errors;turbulent model;CFD;spinning

    V448.15+3

    Adoi:10.7638/kqdlxxb-2013.0059

    0258-1825(2014)05-0705-07

    2013-05-31;

    2013-08-07

    新世紀(jì)優(yōu)秀人才支持計(jì)劃資助(NCET-10-0075)

    陳東陽(yáng)(1988-)男,漢族,江蘇連云港人,碩士研究生,研究方向?yàn)榱鞴恬詈希瑲鈩?dòng)熱彈性仿真.E-mail:cdy_1988@sina.cn

    Laith K.Abbas(1965-),男,伊拉克人,教授,博導(dǎo),研究方向?yàn)闅鈩?dòng)彈性、結(jié)構(gòu)力學(xué)、多體系統(tǒng)動(dòng)力學(xué).E-mail:laithabbas@yahoo.com

    陳東陽(yáng),Laith K.Abbas,芮筱亭,等.結(jié)構(gòu)誤差對(duì)旋轉(zhuǎn)穩(wěn)定彈丸氣動(dòng)特性影響的數(shù)值模擬[J].空氣動(dòng)力學(xué)學(xué)報(bào),2014,32(5):705-711.

    10.7638/kqdlxxb-2013.0059. CHEN D Y,Laith K Abbas,RUI X T,et al.Numerical simulation of a spinning stabilized projectile aerodynamic characteristics effected by structure errors[J].ACTA Aerodynamica Sinica,2014,32(5):705-711.

    猜你喜歡
    彈體湍流氣動(dòng)
    中寰氣動(dòng)執(zhí)行機(jī)構(gòu)
    尾錐角對(duì)彈體斜侵徹過(guò)程中姿態(tài)的影響研究
    橢圓截面彈體斜侵徹金屬靶體彈道研究*
    爆炸與沖擊(2022年2期)2022-03-17 07:28:44
    基于NACA0030的波紋狀翼型氣動(dòng)特性探索
    重氣瞬時(shí)泄漏擴(kuò)散的湍流模型驗(yàn)證
    基于反饋線性化的RLV氣動(dòng)控制一體化設(shè)計(jì)
    STOPAQ粘彈體技術(shù)在管道施工中的應(yīng)用
    上海煤氣(2018年6期)2018-03-07 01:03:22
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    KJH101-127型氣動(dòng)司控道岔的改造
    亚洲av不卡在线观看| 欧美成狂野欧美在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产伦人伦偷精品视频| 亚洲国产精品成人综合色| 真人一进一出gif抽搐免费| 亚洲av熟女| 日本免费a在线| 禁无遮挡网站| 国产成人av教育| xxxwww97欧美| 久久久久性生活片| xxxwww97欧美| 一级毛片女人18水好多| 免费人成在线观看视频色| 搡老熟女国产l中国老女人| 91在线观看av| 在线天堂最新版资源| 久久欧美精品欧美久久欧美| 婷婷精品国产亚洲av在线| 少妇丰满av| 国产高清videossex| 国产精华一区二区三区| avwww免费| 亚洲精品粉嫩美女一区| 欧美一区二区国产精品久久精品| 男女做爰动态图高潮gif福利片| 久久久久精品国产欧美久久久| 亚洲av免费高清在线观看| www.999成人在线观看| av国产免费在线观看| 成年版毛片免费区| 免费在线观看日本一区| 色综合欧美亚洲国产小说| 99久久无色码亚洲精品果冻| 亚洲aⅴ乱码一区二区在线播放| 怎么达到女性高潮| 少妇的丰满在线观看| 日韩欧美精品免费久久 | 老汉色av国产亚洲站长工具| 大型黄色视频在线免费观看| 婷婷亚洲欧美| 一区二区三区国产精品乱码| 国产又黄又爽又无遮挡在线| 18禁在线播放成人免费| 国产av在哪里看| 一进一出抽搐gif免费好疼| 国产综合懂色| 在线观看美女被高潮喷水网站 | 亚洲国产精品999在线| 女生性感内裤真人,穿戴方法视频| 俺也久久电影网| 91麻豆av在线| 国产乱人伦免费视频| 在线视频色国产色| 亚洲欧美日韩卡通动漫| 亚洲国产中文字幕在线视频| 欧美最新免费一区二区三区 | 丰满人妻一区二区三区视频av | 国产免费男女视频| 男人舔女人下体高潮全视频| 国产精品久久久久久人妻精品电影| 最新在线观看一区二区三区| 熟妇人妻久久中文字幕3abv| 亚洲精品成人久久久久久| 亚洲性夜色夜夜综合| 一级作爱视频免费观看| a级一级毛片免费在线观看| 国产午夜福利久久久久久| 久久久成人免费电影| 日韩av在线大香蕉| 高清在线国产一区| 亚洲 欧美 日韩 在线 免费| 国产精品一区二区三区四区久久| 亚洲人与动物交配视频| 亚洲午夜理论影院| 色在线成人网| 亚洲欧美日韩卡通动漫| 天堂影院成人在线观看| 午夜福利在线观看免费完整高清在 | 亚洲欧美精品综合久久99| 少妇的逼水好多| 悠悠久久av| 久久久成人免费电影| 国产精品 国内视频| 精品久久久久久成人av| 日本成人三级电影网站| 国产国拍精品亚洲av在线观看 | 一个人看的www免费观看视频| 18禁黄网站禁片免费观看直播| 欧美性感艳星| 蜜桃亚洲精品一区二区三区| 99久久九九国产精品国产免费| 国产一区在线观看成人免费| 嫩草影院入口| 国产高清三级在线| 日韩中文字幕欧美一区二区| 国产高清视频在线观看网站| 国产美女午夜福利| 国产午夜精品论理片| a级一级毛片免费在线观看| 蜜桃亚洲精品一区二区三区| 国产中年淑女户外野战色| 久久性视频一级片| 极品教师在线免费播放| 国产91精品成人一区二区三区| 久久久久久久久大av| 九九久久精品国产亚洲av麻豆| 黄色女人牲交| 男女下面进入的视频免费午夜| eeuss影院久久| av国产免费在线观看| 久久香蕉精品热| 99热精品在线国产| 高潮久久久久久久久久久不卡| 国产高清videossex| 精品久久久久久久久久久久久| 国产v大片淫在线免费观看| 亚洲一区高清亚洲精品| 欧美黄色淫秽网站| 久久草成人影院| 成人无遮挡网站| 男人舔女人下体高潮全视频| 日本免费一区二区三区高清不卡| 久久精品综合一区二区三区| 一级作爱视频免费观看| 内射极品少妇av片p| 美女大奶头视频| 搡老妇女老女人老熟妇| 免费看十八禁软件| 丁香欧美五月| 国产午夜福利久久久久久| 日本三级黄在线观看| 有码 亚洲区| 国产欧美日韩一区二区三| 亚洲精品日韩av片在线观看 | 免费电影在线观看免费观看| 国产午夜精品久久久久久一区二区三区 | 两个人看的免费小视频| 欧美极品一区二区三区四区| 中文在线观看免费www的网站| 亚洲欧美日韩东京热| 一级a爱片免费观看的视频| 中文字幕人妻丝袜一区二区| 日韩欧美三级三区| 又粗又爽又猛毛片免费看| 真人一进一出gif抽搐免费| 欧美色视频一区免费| 国产精品久久久久久久久免 | 国产精品三级大全| 成人无遮挡网站| 特大巨黑吊av在线直播| 国产色婷婷99| 亚洲va日本ⅴa欧美va伊人久久| 中文亚洲av片在线观看爽| 久久99热这里只有精品18| 香蕉av资源在线| 麻豆国产av国片精品| 日韩欧美在线乱码| 一进一出好大好爽视频| 老司机午夜十八禁免费视频| 国产老妇女一区| 久久精品影院6| 成人鲁丝片一二三区免费| 网址你懂的国产日韩在线| 国产精品免费一区二区三区在线| 亚洲人成伊人成综合网2020| 波野结衣二区三区在线 | 特大巨黑吊av在线直播| 午夜福利在线观看免费完整高清在 | 午夜影院日韩av| 欧美成人性av电影在线观看| 观看免费一级毛片| 亚洲自拍偷在线| 最近最新免费中文字幕在线| 精品国内亚洲2022精品成人| 亚洲欧美日韩无卡精品| 在线观看日韩欧美| 欧美日韩精品网址| 女生性感内裤真人,穿戴方法视频| 免费在线观看成人毛片| 999久久久精品免费观看国产| 一二三四社区在线视频社区8| 成人鲁丝片一二三区免费| 九九在线视频观看精品| 亚洲av第一区精品v没综合| 两个人看的免费小视频| 男人舔女人下体高潮全视频| bbb黄色大片| 一个人看视频在线观看www免费 | 少妇的逼水好多| 男女下面进入的视频免费午夜| 好男人在线观看高清免费视频| 欧美乱色亚洲激情| 男女床上黄色一级片免费看| 亚洲av不卡在线观看| 日本黄大片高清| 偷拍熟女少妇极品色| 日本三级黄在线观看| 中文字幕av成人在线电影| 91在线精品国自产拍蜜月 | 日韩欧美免费精品| 中文字幕久久专区| 亚洲精品日韩av片在线观看 | 免费av不卡在线播放| 少妇的逼好多水| 亚洲精品456在线播放app | 婷婷亚洲欧美| 欧美性感艳星| 国产高清三级在线| 热99在线观看视频| 别揉我奶头~嗯~啊~动态视频| 一个人看视频在线观看www免费 | 又黄又粗又硬又大视频| 午夜免费观看网址| 国产v大片淫在线免费观看| 日本黄色视频三级网站网址| 国产一区二区亚洲精品在线观看| 国产久久久一区二区三区| 国产精品综合久久久久久久免费| 一区二区三区免费毛片| 小蜜桃在线观看免费完整版高清| 亚洲18禁久久av| 成人鲁丝片一二三区免费| 又紧又爽又黄一区二区| 国内毛片毛片毛片毛片毛片| 淫妇啪啪啪对白视频| av中文乱码字幕在线| 久久亚洲真实| 久久久久免费精品人妻一区二区| 亚洲熟妇中文字幕五十中出| 国产精品久久电影中文字幕| 日本与韩国留学比较| a在线观看视频网站| 午夜日韩欧美国产| 19禁男女啪啪无遮挡网站| 99久久综合精品五月天人人| 精品电影一区二区在线| 午夜福利视频1000在线观看| 激情在线观看视频在线高清| 欧美性感艳星| 中文字幕高清在线视频| 99精品欧美一区二区三区四区| 亚洲内射少妇av| 国产精品久久久久久久久免 | 日韩免费av在线播放| 国产精品久久视频播放| 亚洲va日本ⅴa欧美va伊人久久| 国产av不卡久久| 国产精品 欧美亚洲| 嫩草影视91久久| 亚洲黑人精品在线| 亚洲va日本ⅴa欧美va伊人久久| 19禁男女啪啪无遮挡网站| svipshipincom国产片| 久久久久亚洲av毛片大全| 高清毛片免费观看视频网站| 少妇高潮的动态图| 日韩欧美免费精品| 亚洲成av人片在线播放无| 香蕉av资源在线| 狂野欧美激情性xxxx| 成人特级黄色片久久久久久久| 舔av片在线| 在线看三级毛片| 99久久九九国产精品国产免费| 十八禁网站免费在线| 青草久久国产| 久久久久国产精品人妻aⅴ院| 国产一区二区三区在线臀色熟女| 啦啦啦观看免费观看视频高清| 精品久久久久久,| 亚洲国产色片| 亚洲一区二区三区色噜噜| 亚洲国产精品成人综合色| 中文字幕av在线有码专区| 日韩欧美一区二区三区在线观看| 91麻豆av在线| 亚洲精品日韩av片在线观看 | 国产一区二区三区视频了| 亚洲欧美日韩卡通动漫| netflix在线观看网站| 88av欧美| 国产高清三级在线| 在线免费观看的www视频| 老司机福利观看| 最好的美女福利视频网| 一级作爱视频免费观看| 我的老师免费观看完整版| 18美女黄网站色大片免费观看| 在线播放无遮挡| 五月玫瑰六月丁香| 久久久精品大字幕| 99精品欧美一区二区三区四区| 色尼玛亚洲综合影院| 亚洲真实伦在线观看| 毛片女人毛片| h日本视频在线播放| 欧美性感艳星| 亚洲一区二区三区不卡视频| 亚洲av电影在线进入| 色av中文字幕| 18禁黄网站禁片免费观看直播| 免费在线观看成人毛片| 国产精品影院久久| 精品熟女少妇八av免费久了| 99视频精品全部免费 在线| 免费人成视频x8x8入口观看| 欧美不卡视频在线免费观看| 露出奶头的视频| 天堂影院成人在线观看| 高清日韩中文字幕在线| 国产又黄又爽又无遮挡在线| 国产高清激情床上av| 欧美中文日本在线观看视频| 日本五十路高清| 成人性生交大片免费视频hd| 精品国产超薄肉色丝袜足j| 夜夜爽天天搞| 国产精品99久久久久久久久| 国产精品爽爽va在线观看网站| 午夜精品一区二区三区免费看| 亚洲欧美日韩无卡精品| 性欧美人与动物交配| 大型黄色视频在线免费观看| 亚洲成人久久爱视频| 日韩中文字幕欧美一区二区| 午夜老司机福利剧场| 老司机深夜福利视频在线观看| 久久久久久国产a免费观看| 国产在线精品亚洲第一网站| 中文字幕人妻丝袜一区二区| 国产又黄又爽又无遮挡在线| 高清日韩中文字幕在线| 两个人看的免费小视频| 久久精品国产清高在天天线| 精品久久久久久久久久久久久| 一级a爱片免费观看的视频| 热99re8久久精品国产| 精品国产亚洲在线| 十八禁人妻一区二区| 免费无遮挡裸体视频| 日韩欧美免费精品| 亚洲无线观看免费| 免费看美女性在线毛片视频| 中文资源天堂在线| 亚洲狠狠婷婷综合久久图片| 18+在线观看网站| 欧美黑人欧美精品刺激| 尤物成人国产欧美一区二区三区| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 亚洲欧美日韩无卡精品| 免费看日本二区| 久久国产乱子伦精品免费另类| 亚洲人成伊人成综合网2020| 嫩草影视91久久| 狂野欧美白嫩少妇大欣赏| 夜夜看夜夜爽夜夜摸| 亚洲精华国产精华精| 手机成人av网站| 日本a在线网址| 中文字幕av在线有码专区| 免费高清视频大片| 一级毛片高清免费大全| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 国产精品免费一区二区三区在线| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 亚洲男人的天堂狠狠| 精品国产三级普通话版| 高清在线国产一区| 色吧在线观看| 最近最新中文字幕大全电影3| 国产欧美日韩一区二区三| 内地一区二区视频在线| 欧美bdsm另类| 日韩欧美在线二视频| 美女大奶头视频| 桃色一区二区三区在线观看| 日韩欧美精品v在线| 在线看三级毛片| 欧美又色又爽又黄视频| e午夜精品久久久久久久| 人人妻,人人澡人人爽秒播| 欧美日韩一级在线毛片| 婷婷精品国产亚洲av在线| 国产美女午夜福利| 神马国产精品三级电影在线观看| e午夜精品久久久久久久| 日韩成人在线观看一区二区三区| 国产v大片淫在线免费观看| 国产一区在线观看成人免费| 亚洲国产精品久久男人天堂| 香蕉av资源在线| 国产综合懂色| 麻豆国产97在线/欧美| 天堂网av新在线| 亚洲av美国av| 脱女人内裤的视频| 有码 亚洲区| 亚洲国产欧美网| 色吧在线观看| 日日摸夜夜添夜夜添小说| 亚洲熟妇熟女久久| 在线国产一区二区在线| 在线观看美女被高潮喷水网站 | 亚洲中文字幕日韩| 麻豆国产av国片精品| 日韩人妻高清精品专区| 黄色片一级片一级黄色片| 人人妻,人人澡人人爽秒播| 97超级碰碰碰精品色视频在线观看| 无遮挡黄片免费观看| 国产中年淑女户外野战色| 村上凉子中文字幕在线| 成人亚洲精品av一区二区| eeuss影院久久| 亚洲第一电影网av| 丰满人妻一区二区三区视频av | 国产欧美日韩一区二区三| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产亚洲精品综合一区在线观看| 国产精品爽爽va在线观看网站| 国产亚洲精品av在线| 51午夜福利影视在线观看| 日韩国内少妇激情av| 丝袜美腿在线中文| 在线观看免费视频日本深夜| 亚洲欧美日韩高清专用| 久久久久久人人人人人| 97超级碰碰碰精品色视频在线观看| 美女大奶头视频| 亚洲久久久久久中文字幕| 色吧在线观看| 国产av麻豆久久久久久久| 久久久国产精品麻豆| 国产熟女xx| 亚洲av第一区精品v没综合| 老熟妇乱子伦视频在线观看| 男插女下体视频免费在线播放| 12—13女人毛片做爰片一| 无人区码免费观看不卡| 日韩欧美精品v在线| 色av中文字幕| 88av欧美| 午夜日韩欧美国产| 欧美成人性av电影在线观看| 亚洲国产精品合色在线| 久久久久久久午夜电影| 欧美高清成人免费视频www| 黄色视频,在线免费观看| 国产伦在线观看视频一区| 亚洲av熟女| 在线a可以看的网站| 一边摸一边抽搐一进一小说| 99视频精品全部免费 在线| 亚洲国产色片| 99国产极品粉嫩在线观看| 精品日产1卡2卡| 色综合亚洲欧美另类图片| 国产 一区 欧美 日韩| 深夜精品福利| 久久久久久久久中文| 深爱激情五月婷婷| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久视频播放| 亚洲国产高清在线一区二区三| 国产精华一区二区三区| 日本五十路高清| 亚洲精品456在线播放app | 久久人人精品亚洲av| 99久国产av精品| 丰满的人妻完整版| 窝窝影院91人妻| 久久久久久久午夜电影| 伊人久久精品亚洲午夜| 成人无遮挡网站| 夜夜看夜夜爽夜夜摸| av女优亚洲男人天堂| 日本五十路高清| 亚洲欧美精品综合久久99| 国产精品,欧美在线| 国产成人a区在线观看| 亚洲av成人av| 欧美日本视频| 午夜老司机福利剧场| 三级国产精品欧美在线观看| 欧美中文日本在线观看视频| 欧美午夜高清在线| 啦啦啦韩国在线观看视频| 中文字幕熟女人妻在线| 久久久久免费精品人妻一区二区| 欧美黑人巨大hd| 身体一侧抽搐| 在线免费观看不下载黄p国产 | 1000部很黄的大片| 久久久久性生活片| 国产不卡一卡二| 一本精品99久久精品77| 国产精品三级大全| 成人高潮视频无遮挡免费网站| 午夜福利视频1000在线观看| ponron亚洲| 免费在线观看亚洲国产| 国产一区在线观看成人免费| 精品人妻一区二区三区麻豆 | 国产单亲对白刺激| 国产高清激情床上av| 亚洲欧美日韩高清在线视频| 亚洲av日韩精品久久久久久密| 少妇熟女aⅴ在线视频| 高清毛片免费观看视频网站| 91字幕亚洲| 国产单亲对白刺激| 免费电影在线观看免费观看| 亚洲av中文字字幕乱码综合| 久久亚洲真实| 我的老师免费观看完整版| 免费看日本二区| 国产欧美日韩一区二区三| 免费人成在线观看视频色| 观看美女的网站| 熟女电影av网| 热99re8久久精品国产| 搡女人真爽免费视频火全软件 | 亚洲美女黄片视频| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久久黄片| 看黄色毛片网站| 久久天躁狠狠躁夜夜2o2o| 午夜福利在线观看免费完整高清在 | 国产精华一区二区三区| 嫩草影视91久久| 校园春色视频在线观看| 五月玫瑰六月丁香| 天天一区二区日本电影三级| 亚洲精品久久国产高清桃花| 男人舔奶头视频| 熟女少妇亚洲综合色aaa.| 亚洲18禁久久av| 搞女人的毛片| 欧美成人一区二区免费高清观看| 午夜福利高清视频| 亚洲男人的天堂狠狠| 一二三四社区在线视频社区8| 日本熟妇午夜| 欧美乱妇无乱码| 在线视频色国产色| 亚洲在线自拍视频| 中文字幕熟女人妻在线| 精品不卡国产一区二区三区| 国产亚洲欧美98| 一进一出好大好爽视频| 亚洲一区二区三区不卡视频| www.www免费av| 亚洲精品一区av在线观看| 高清在线国产一区| 亚洲av日韩精品久久久久久密| 亚洲专区国产一区二区| 国产真实乱freesex| 色噜噜av男人的天堂激情| 看片在线看免费视频| 欧美zozozo另类| 久久精品人妻少妇| 特大巨黑吊av在线直播| 色综合站精品国产| 国产乱人伦免费视频| 国产美女午夜福利| 久久久久久久久大av| 日韩欧美在线乱码| 日本一二三区视频观看| 99国产极品粉嫩在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 麻豆国产97在线/欧美| 欧美精品啪啪一区二区三区| 午夜精品久久久久久毛片777| 中文字幕人成人乱码亚洲影| 一区二区三区国产精品乱码| 国产精品 国内视频| 人人妻人人澡欧美一区二区| 在线播放无遮挡| 夜夜爽天天搞| 免费看a级黄色片| 国产高清三级在线| 香蕉丝袜av| 色哟哟哟哟哟哟| 日韩欧美在线二视频| 亚洲国产欧美人成| av天堂中文字幕网| 免费在线观看成人毛片| 日本精品一区二区三区蜜桃| 欧美乱色亚洲激情| 午夜福利视频1000在线观看| 国产一区二区在线观看日韩 | 色吧在线观看| 日韩国内少妇激情av| 欧美色视频一区免费| 久久精品国产亚洲av香蕉五月| netflix在线观看网站| 男插女下体视频免费在线播放| 五月伊人婷婷丁香| 看片在线看免费视频| 黄片大片在线免费观看| 无遮挡黄片免费观看| av专区在线播放| 欧美激情在线99| 国产av在哪里看| 久久国产精品人妻蜜桃| 91在线精品国自产拍蜜月 | 色av中文字幕| 成人特级av手机在线观看| 在线a可以看的网站| 亚洲第一欧美日韩一区二区三区|