• <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一区二区精品久久| 亚洲,欧美,日韩| 免费av不卡在线播放| 人人妻人人爽人人添夜夜欢视频 | 秋霞伦理黄片| 97超视频在线观看视频| 久久精品熟女亚洲av麻豆精品| 久久99热这里只频精品6学生| av网站免费在线观看视频| 亚洲综合色惰| 亚洲精品第二区| 国内精品宾馆在线| 精品国产一区二区久久| 亚洲国产欧美日韩在线播放 | 伦理电影免费视频| 亚洲欧美精品专区久久| 国产乱人偷精品视频| 亚洲欧美日韩卡通动漫| 久久久亚洲精品成人影院| 人妻制服诱惑在线中文字幕| 免费人成在线观看视频色| 99热全是精品| 久久久精品免费免费高清| 亚洲国产日韩一区二区| 国产一区二区三区综合在线观看 | 亚洲人成网站在线观看播放| 两个人的视频大全免费| 亚洲国产欧美在线一区| 男男h啪啪无遮挡| 人妻系列 视频| 99热这里只有是精品50| 插逼视频在线观看| www.av在线官网国产| 亚洲精品国产av成人精品| 美女cb高潮喷水在线观看| av免费观看日本| h视频一区二区三区| 大码成人一级视频| 免费大片黄手机在线观看| 亚洲av日韩在线播放| 亚洲欧洲国产日韩| 在线观看免费高清a一片| 亚洲av综合色区一区| 少妇被粗大的猛进出69影院 | 精品人妻熟女av久视频| 免费黄色在线免费观看| 亚洲图色成人| 国产在线免费精品| 99久久精品一区二区三区| 欧美3d第一页| 中文字幕免费在线视频6| 中文在线观看免费www的网站| 国产乱人偷精品视频| h视频一区二区三区| 久久这里有精品视频免费| a级片在线免费高清观看视频| 日本黄色片子视频| 高清毛片免费看| 日韩精品免费视频一区二区三区 | 婷婷色av中文字幕| 亚洲中文av在线| av福利片在线| 91在线精品国自产拍蜜月| 在线观看国产h片| 精品卡一卡二卡四卡免费| 自线自在国产av| 日韩熟女老妇一区二区性免费视频| 精品国产乱码久久久久久小说| 天美传媒精品一区二区| 十八禁高潮呻吟视频 | 亚州av有码| 亚洲精品国产成人久久av| 日韩 亚洲 欧美在线| 精品视频人人做人人爽| 日日啪夜夜撸| 丰满迷人的少妇在线观看| 国产视频内射| 欧美日本中文国产一区发布| 久久久久久伊人网av| 爱豆传媒免费全集在线观看| 免费高清在线观看视频在线观看| 国产一区二区三区av在线| 精品国产一区二区三区久久久樱花| 日韩强制内射视频| 少妇 在线观看| 各种免费的搞黄视频| 波野结衣二区三区在线| 少妇高潮的动态图| 日韩熟女老妇一区二区性免费视频| 成人综合一区亚洲| 午夜福利影视在线免费观看| 国产男人的电影天堂91| 日韩av在线免费看完整版不卡| 人妻一区二区av| 国产一区有黄有色的免费视频| 另类精品久久| 国产在视频线精品| 国产欧美另类精品又又久久亚洲欧美| 国产精品伦人一区二区| 高清不卡的av网站| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美精品一区二区大全| 国产黄片美女视频| 久久女婷五月综合色啪小说| 亚洲国产精品成人久久小说| 久久人人爽av亚洲精品天堂| 一个人看视频在线观看www免费| 视频中文字幕在线观看| av一本久久久久| √禁漫天堂资源中文www| 久久久久久伊人网av| 七月丁香在线播放| 国产成人精品无人区| www.av在线官网国产| 国产欧美日韩精品一区二区| 国产视频内射| 能在线免费看毛片的网站| 免费观看在线日韩| 亚洲精华国产精华液的使用体验| 又黄又爽又刺激的免费视频.| 国产深夜福利视频在线观看| 国产精品久久久久久精品电影小说| 丰满饥渴人妻一区二区三| 天堂俺去俺来也www色官网| 国产精品久久久久久精品古装| 国产免费又黄又爽又色| 日韩一区二区视频免费看| 80岁老熟妇乱子伦牲交| 色视频www国产| 中文字幕av电影在线播放| 国模一区二区三区四区视频| 国国产精品蜜臀av免费| 精品午夜福利在线看| 亚洲精品乱码久久久久久按摩| 国产一区二区在线观看日韩| 日本色播在线视频| 99热国产这里只有精品6| 有码 亚洲区| tube8黄色片| 亚洲婷婷狠狠爱综合网| 久久久久久伊人网av| 97在线人人人人妻| 久热久热在线精品观看| 亚洲av免费高清在线观看| 午夜久久久在线观看| 免费高清在线观看视频在线观看| 亚洲成人av在线免费| 女性被躁到高潮视频| 少妇裸体淫交视频免费看高清| 日韩三级伦理在线观看| 精品一区在线观看国产| 久久精品国产亚洲网站| 国产中年淑女户外野战色| 伊人久久国产一区二区| 国产视频内射| 青青草视频在线视频观看| 最近的中文字幕免费完整| 涩涩av久久男人的天堂| 国产午夜精品一二区理论片| 免费播放大片免费观看视频在线观看| 亚洲精品国产色婷婷电影| 一级av片app| 欧美变态另类bdsm刘玥| 如日韩欧美国产精品一区二区三区 | 高清av免费在线| 亚洲欧美清纯卡通| 亚洲av在线观看美女高潮| 欧美日韩一区二区视频在线观看视频在线| 国产毛片在线视频| 日韩一区二区视频免费看| 在线播放无遮挡| 深夜a级毛片| 在线观看美女被高潮喷水网站| 国产亚洲精品久久久com| 精品亚洲成国产av| 三级国产精品欧美在线观看| 精品久久久久久久久av| 国产视频首页在线观看| 久久精品国产a三级三级三级| 国产探花极品一区二区| 亚洲伊人久久精品综合| 王馨瑶露胸无遮挡在线观看| 一级毛片aaaaaa免费看小| .国产精品久久| 99热国产这里只有精品6| 18禁在线无遮挡免费观看视频| 丝袜在线中文字幕| 精品国产露脸久久av麻豆| 午夜av观看不卡| 欧美日韩国产mv在线观看视频| 国产精品久久久久久久久免| 国产精品99久久久久久久久| 26uuu在线亚洲综合色| 9色porny在线观看| 丝袜喷水一区| 成年人免费黄色播放视频 | av免费在线看不卡| 日本av手机在线免费观看| 亚洲不卡免费看| 色94色欧美一区二区| 男女边吃奶边做爰视频| 亚洲欧美一区二区三区国产| 久久久久久久国产电影| 日韩av不卡免费在线播放| 亚洲av国产av综合av卡| 欧美精品一区二区免费开放| 日韩成人av中文字幕在线观看| 免费观看性生交大片5| 亚洲精品一区蜜桃| 夫妻性生交免费视频一级片| 亚洲精品国产成人久久av| 精品少妇内射三级| 国产免费一区二区三区四区乱码| 99re6热这里在线精品视频| 日本-黄色视频高清免费观看| 亚洲激情五月婷婷啪啪| 少妇被粗大猛烈的视频| 高清不卡的av网站| 国产av一区二区精品久久| 亚洲精品自拍成人| 国产免费一区二区三区四区乱码| 久久6这里有精品| 国产黄片视频在线免费观看| 国产精品久久久久久av不卡| 国语对白做爰xxxⅹ性视频网站| 国产乱人偷精品视频| 天天操日日干夜夜撸| av天堂久久9| 国产一级毛片在线| 国产精品一区二区三区四区免费观看| 国产黄频视频在线观看| 国产精品熟女久久久久浪| 男人添女人高潮全过程视频| 看非洲黑人一级黄片| av国产精品久久久久影院| 成人综合一区亚洲| 伦精品一区二区三区| 国产在线男女| 欧美bdsm另类| 国产免费一区二区三区四区乱码| a级片在线免费高清观看视频| 亚洲激情五月婷婷啪啪| 免费久久久久久久精品成人欧美视频 | 日韩三级伦理在线观看| 久久久久久人妻| 最近中文字幕2019免费版| 免费观看a级毛片全部| 啦啦啦中文免费视频观看日本| 国产精品嫩草影院av在线观看| 免费看光身美女| 精品国产乱码久久久久久小说| 欧美少妇被猛烈插入视频| 国产黄片视频在线免费观看| 精品一区在线观看国产| 婷婷色综合大香蕉| 精品亚洲成a人片在线观看| 免费看不卡的av| 国产精品一区二区在线观看99| 色94色欧美一区二区| 乱码一卡2卡4卡精品| 亚洲人成网站在线播| 国产日韩欧美亚洲二区| 另类亚洲欧美激情| 久久99一区二区三区| 汤姆久久久久久久影院中文字幕| 日本色播在线视频| 免费人妻精品一区二区三区视频| 亚洲人成电影观看| 国产成人欧美在线观看 | bbb黄色大片| 亚洲av电影在线观看一区二区三区| 老司机福利观看| 成人影院久久| 一区二区三区乱码不卡18| 国产精品99久久99久久久不卡| av在线app专区| 国产伦人伦偷精品视频| 亚洲欧洲精品一区二区精品久久久| 亚洲中文字幕日韩| 国产免费福利视频在线观看| 成人影院久久| 999精品在线视频| 中国美女看黄片| 久久青草综合色| 考比视频在线观看| 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区 | 亚洲欧美日韩高清在线视频 | 80岁老熟妇乱子伦牲交| 中文字幕色久视频| 99热全是精品| 永久免费av网站大全| 久久久精品94久久精品| 国产日韩欧美亚洲二区| 精品人妻1区二区| 我要看黄色一级片免费的| 99精国产麻豆久久婷婷| 美女高潮喷水抽搐中文字幕| 久久亚洲精品不卡| 两个人免费观看高清视频| 搡老岳熟女国产| 巨乳人妻的诱惑在线观看| 精品亚洲成a人片在线观看| 在线观看人妻少妇| 另类亚洲欧美激情| 日韩三级视频一区二区三区| 80岁老熟妇乱子伦牲交| 成人免费观看视频高清| 视频区欧美日本亚洲| 99国产综合亚洲精品| 电影成人av| 极品人妻少妇av视频| 日本黄色日本黄色录像| 欧美性长视频在线观看| 国产一区二区 视频在线| 美女福利国产在线| av电影中文网址| 午夜激情av网站| 午夜精品国产一区二区电影| 欧美av亚洲av综合av国产av| 国产亚洲精品第一综合不卡| 黄片播放在线免费| 国产亚洲精品一区二区www | 五月天丁香电影| 岛国在线观看网站| 人人妻人人添人人爽欧美一区卜| 99国产精品99久久久久| 一本大道久久a久久精品| 亚洲精品成人av观看孕妇| av线在线观看网站| 一区二区三区激情视频| 天堂俺去俺来也www色官网| 狠狠精品人妻久久久久久综合| 老司机靠b影院| 伊人久久大香线蕉亚洲五| 国产深夜福利视频在线观看| 亚洲全国av大片| 婷婷色av中文字幕| 国产精品99久久99久久久不卡| 精品久久蜜臀av无| 操出白浆在线播放| 777米奇影视久久| 亚洲精品自拍成人| 母亲3免费完整高清在线观看| 51午夜福利影视在线观看| 国产亚洲精品久久久久5区| 99精国产麻豆久久婷婷| 国产精品 国内视频| 日韩大码丰满熟妇| 欧美日韩视频精品一区| 99国产精品一区二区蜜桃av | 亚洲综合色网址| 国产欧美日韩一区二区精品| 久久香蕉激情| 一本—道久久a久久精品蜜桃钙片| 国产片内射在线| 国产av国产精品国产| 黄色怎么调成土黄色| 王馨瑶露胸无遮挡在线观看| 日韩 亚洲 欧美在线| 另类精品久久| www.自偷自拍.com| 亚洲国产毛片av蜜桃av| 精品亚洲成国产av| 国产又色又爽无遮挡免| 国产欧美日韩精品亚洲av| 欧美日韩亚洲综合一区二区三区_| 大型av网站在线播放| 亚洲精品中文字幕一二三四区 | 久久久精品94久久精品| 正在播放国产对白刺激| 免费一级毛片在线播放高清视频 | 国产色视频综合| 亚洲黑人精品在线| 在线观看免费视频网站a站| 少妇被粗大的猛进出69影院| 国产在线免费精品| 欧美午夜高清在线| 男人舔女人的私密视频| 妹子高潮喷水视频| 国产免费一区二区三区四区乱码| 丝袜人妻中文字幕| 欧美激情高清一区二区三区| 亚洲,欧美精品.| 亚洲熟女毛片儿| 夜夜骑夜夜射夜夜干| 欧美国产精品一级二级三级| 免费人妻精品一区二区三区视频| 国产欧美日韩一区二区三区在线| 国产免费视频播放在线视频| 久久久精品94久久精品| 日韩有码中文字幕| 免费看十八禁软件| 五月开心婷婷网| 无遮挡黄片免费观看| 精品卡一卡二卡四卡免费| 亚洲欧洲日产国产| 欧美日韩亚洲国产一区二区在线观看 | 欧美人与性动交α欧美精品济南到| 正在播放国产对白刺激| 亚洲专区国产一区二区| 久久精品国产亚洲av香蕉五月 | xxxhd国产人妻xxx| 天堂俺去俺来也www色官网| 久久亚洲精品不卡| 69av精品久久久久久 | 天天添夜夜摸| 91国产中文字幕| 亚洲欧洲日产国产| 在线 av 中文字幕| 乱人伦中国视频| 深夜精品福利| 中国美女看黄片| 亚洲伊人久久精品综合| 国产精品免费大片| 欧美性长视频在线观看| 丝袜喷水一区| 日本a在线网址| 如日韩欧美国产精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 精品福利观看| 日韩电影二区| 中文字幕人妻丝袜制服| 两人在一起打扑克的视频| 国产精品亚洲av一区麻豆| 国产成人免费无遮挡视频| 一边摸一边做爽爽视频免费| 蜜桃在线观看..| 99精国产麻豆久久婷婷| 久久久精品国产亚洲av高清涩受| 日日夜夜操网爽| 天堂俺去俺来也www色官网| 人妻一区二区av| tocl精华| 精品欧美一区二区三区在线| 少妇猛男粗大的猛烈进出视频| 欧美日韩亚洲高清精品| 亚洲九九香蕉| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区久久| 女人被躁到高潮嗷嗷叫费观| 国产一区二区三区综合在线观看| av福利片在线| 在线观看免费高清a一片| 丝袜美足系列| 嫩草影视91久久| 国产精品一区二区在线不卡| 国产色视频综合| 人妻久久中文字幕网| 色播在线永久视频| 国产精品.久久久| 秋霞在线观看毛片| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| tube8黄色片| 一区二区av电影网| 不卡一级毛片| 纵有疾风起免费观看全集完整版| 日日爽夜夜爽网站| 精品亚洲成a人片在线观看| 欧美精品一区二区免费开放| 黄色视频在线播放观看不卡| 久久人妻福利社区极品人妻图片| 日日夜夜操网爽| 美女高潮喷水抽搐中文字幕| 午夜激情久久久久久久| 国产在线一区二区三区精| 一区二区三区精品91| 美女国产高潮福利片在线看| 国产精品自产拍在线观看55亚洲 | 岛国在线观看网站| 欧美av亚洲av综合av国产av| 国产伦人伦偷精品视频| 欧美黄色淫秽网站| 在线观看免费高清a一片| 18禁观看日本| 国产男女超爽视频在线观看| 桃花免费在线播放| 欧美精品啪啪一区二区三区 | 国产av一区二区精品久久| 午夜福利在线观看吧| 亚洲情色 制服丝袜| 午夜精品国产一区二区电影| 大香蕉久久网| 操美女的视频在线观看| 丝袜脚勾引网站| 另类亚洲欧美激情| 欧美日韩一级在线毛片| 大片免费播放器 马上看| 久久国产精品男人的天堂亚洲| 亚洲人成电影观看| 精品少妇久久久久久888优播| 每晚都被弄得嗷嗷叫到高潮| 免费高清在线观看日韩| 精品久久久久久电影网| 欧美人与性动交α欧美软件| 人人妻人人澡人人爽人人夜夜| 久久久久久久精品精品| 国产精品久久久人人做人人爽| 欧美黄色淫秽网站| 天天躁夜夜躁狠狠躁躁| 国产精品.久久久| 黑人欧美特级aaaaaa片| www.自偷自拍.com| 久久99一区二区三区| 久久久久久久大尺度免费视频| 成人国产一区最新在线观看| 亚洲一码二码三码区别大吗| 久久久久久久久免费视频了| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕一级| 男女国产视频网站| 成年动漫av网址| 亚洲成人国产一区在线观看| 国产日韩欧美在线精品| 久久精品国产a三级三级三级| 亚洲熟女毛片儿| 欧美日韩精品网址| 丝袜美足系列| 日日摸夜夜添夜夜添小说| 久久久国产一区二区| 精品熟女少妇八av免费久了| 午夜91福利影院| 人人妻,人人澡人人爽秒播| 久久人妻福利社区极品人妻图片| 亚洲avbb在线观看| 亚洲精品成人av观看孕妇| 国产欧美日韩精品亚洲av| 99国产精品一区二区蜜桃av | 男女无遮挡免费网站观看| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久| 亚洲精品第二区| 女人久久www免费人成看片| 色综合欧美亚洲国产小说| 国产精品久久久久久精品电影小说| 男女无遮挡免费网站观看| 一个人免费在线观看的高清视频 | 国产成人影院久久av| 2018国产大陆天天弄谢| 不卡av一区二区三区| 国产男人的电影天堂91| 中亚洲国语对白在线视频| 99国产精品免费福利视频| netflix在线观看网站| 亚洲 欧美一区二区三区| 在线观看舔阴道视频| 亚洲精品成人av观看孕妇| 国产成人啪精品午夜网站| 久久这里只有精品19| 亚洲中文日韩欧美视频| 午夜福利,免费看| 欧美日韩视频精品一区| 久久国产亚洲av麻豆专区| 悠悠久久av| 欧美黑人欧美精品刺激| 1024视频免费在线观看| 老司机午夜福利在线观看视频 | 91麻豆av在线| 看免费av毛片| 一本综合久久免费| 国产日韩欧美亚洲二区| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 五月天丁香电影| 最近中文字幕2019免费版| 91麻豆av在线| 在线 av 中文字幕| 在线观看www视频免费| av欧美777| 亚洲国产欧美一区二区综合| 欧美日韩福利视频一区二区| 精品少妇内射三级| 精品人妻一区二区三区麻豆| 久久国产精品大桥未久av| 女人精品久久久久毛片| 久热爱精品视频在线9| videosex国产| av片东京热男人的天堂| 欧美日韩精品网址| 精品乱码久久久久久99久播| 国产一区二区三区av在线| 99九九在线精品视频| 黄片播放在线免费| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| 高清av免费在线| 美女视频免费永久观看网站| 丝袜美腿诱惑在线| 正在播放国产对白刺激| 日本五十路高清| 亚洲一区二区三区欧美精品| 美女福利国产在线| 亚洲国产毛片av蜜桃av| 亚洲五月婷婷丁香| 搡老熟女国产l中国老女人| 1024香蕉在线观看| 桃花免费在线播放| 久久久精品区二区三区| 国产精品九九99| 亚洲国产欧美在线一区| 亚洲精品一区蜜桃| 汤姆久久久久久久影院中文字幕| 亚洲av片天天在线观看| kizo精华| 女人精品久久久久毛片| 亚洲精品中文字幕在线视频| 最近最新中文字幕大全免费视频| 窝窝影院91人妻|