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

    初始分離條件對航彈與載機分離安全性影響的數(shù)值模擬研究

    2016-04-18 07:40:35雷娟棉牛健平王鎖柱周奇北京理工大學宇航學院北京100081
    兵工學報 2016年2期
    關鍵詞:數(shù)值模擬

    雷娟棉,牛健平,王鎖柱,周奇(北京理工大學宇航學院,北京100081)

    ?

    初始分離條件對航彈與載機分離安全性影響的數(shù)值模擬研究

    雷娟棉,牛健平,王鎖柱,周奇
    (北京理工大學宇航學院,北京100081)

    摘要:為研究投放條件對航彈與載機分離安全性的影響,采用非定常計算流體力學數(shù)值模擬方法和動網(wǎng)格技術,同時耦合求解六自由度彈道方程,對航彈與載機的分離過程進行模擬。給出載機在不同飛行馬赫數(shù)、攻角、側滑角、飛行高度及航彈在不同初始下拋速度、角速度條件下,航彈從載機投放后的分離軌跡和姿態(tài)變化規(guī)律,研究了這些因素對分離安全性的影響。研究結果表明:初始分離過程中載機對航彈有很強的氣動干擾,對航彈的氣動特性、分離軌跡及彈體姿態(tài)影響很大;隨著分離馬赫數(shù)、投放攻角增大,載機對航彈的氣動干擾增強,航彈的分離安全性變差;對于掛載于左側機翼下的航彈,一定的負向側滑角有利于彈體與載機安全分離;飛行高度越高,越有利于航彈與載機安全分離;一定的初始下拋速度和適當?shù)南聮伋跏冀撬俣扔欣诎踩蛛x。

    關鍵詞:兵器科學與技術;機彈分離;數(shù)值模擬;動網(wǎng)格;氣動特性

    0 引言

    隨著國際形勢的發(fā)展變化,機載武器在軍事上的使用越來越廣泛,其外形也越來越復雜,對其性能的要求也越來越高。機載武器的裝載方式由傳統(tǒng)的外掛式向內埋式轉變,投放的速度范圍也在增大。這些因素都增加了載機與武器周圍流場的復雜性,加之載機與投放武器之間的氣動干擾,使得機載武器的分離軌跡更加難以預測。如果不以安全的方式投放,武器可能與載機發(fā)生碰撞,危及載機的安全,也不能確保機載武器精確命中目標。因此,研究機載武器初始分離條件對機彈分離安全性的影響,給出安全投放分離包絡,預測武器投放后的飛行彈道與姿態(tài),具有十分重要的意義。

    保證機載武器與載機之間的安全分離是新型武器研制過程中的關鍵點之一。為了保證機載武器在各種條件下能夠與載機安全、可靠地分離,國內外學者已開展了大量的機彈分離安全性投放飛行試驗、地面試驗及計算流體力學(CFD)數(shù)值模擬研究。飛行試驗成本高昂,也有較高的危險性,所以目前的研究以地面試驗和CFD數(shù)值模擬為主。在解決多體分離問題時,國內主要采用捕獲軌跡系統(tǒng)(CTS)試驗方法,如趙忠良等[1]利用CTS方法,研究了子母彈在高超聲速條件下的分離問題,反映了母彈頭部激波對子彈的干擾特性。雷娟棉等[2]進行了布撒器和子彈間的氣動干擾CTS試驗研究,得出了各擾動參數(shù)隨子彈位置的變化規(guī)律和氣動干擾區(qū)的范圍。此外地面試驗方法還有風洞中的自由釋放試驗[3-4]等。相比于飛行試驗和地面風洞試驗,CFD數(shù)值模擬成本低,效率高,并且能方便地獲得各種數(shù)據(jù),為設計和改型提供依據(jù);且隨著技術的發(fā)展,數(shù)值模擬結果的可信度在不斷提高,逐漸成為研究機彈分離問題的重要手段。采用CFD數(shù)值方法對機彈分離問題進行模擬,主要有以下幾類方法:采用多塊嵌套網(wǎng)格來解決機彈分離問題是一種常用的方法。唐志共等[5]采用非結構重疊網(wǎng)格并求解歐拉方程,模擬了外掛投放問題,計算結果與試驗數(shù)據(jù)吻合較好。張玉東等[6]采用分區(qū)拼接網(wǎng)格數(shù)值模擬方法,聯(lián)立氣動方程和彈道方程,研究了子母彈分離過程中子彈從母彈中分離并穿越激波的過程。李亭鶴等[7]改進了重疊網(wǎng)格技術中的割補法,編制了通用程序并給出了三維算例,達到了較好的效果。采用非結構網(wǎng)格或混合網(wǎng)格求解動態(tài)邊界問題是另一種廣泛使用的方法。雷娟棉等[8]采用分區(qū)非結構網(wǎng)格模擬了某型火箭彈的拋撒分離過程。張軍等[9]采用自適應非結構網(wǎng)格計算了二維多體干擾流場。此外,還有學者使用笛卡爾網(wǎng)格來研究機彈分離問題,如桑為民等[10]應用笛卡爾非結構切割網(wǎng)格,數(shù)值模擬了機翼帶外掛物投放的問題。其他采用數(shù)值模擬方法的例子如文獻[11-17]。

    本文采用的基于彈簧法和局部網(wǎng)格重構的非結構動網(wǎng)格技術,對復雜外形問題具有很好的適應性,并能在計算域大變形的情況下仍保證較好的網(wǎng)格質量,有利于提高計算精度。本文研究了各種分離條件對航彈與載機分離安全性的影響,對各種工況進行了比較全面的數(shù)值仿真,得到了具有一定可信度的結果,為空投試驗時初始分離條件的確定提供了支撐,也可作為同類產(chǎn)品的工程設計與實際使用時的參考。

    1 數(shù)值方法

    本文以非定常三維Navier-Stokes方程為基礎,采用Spalart-Allmaras湍流模型,借助非結構動網(wǎng)格及其局部重構技術,耦合求解六自由度彈道方程,在不同條件下對航彈與載機的初始分離過程進行了數(shù)值模擬。

    1.1控制方程

    在航彈與載機分離的過程中,彈體表面作為流場域的內邊界一直在運動。而動網(wǎng)格技術可用于解決由于流場邊界運動造成流場域隨時間變化的流場計算問題。本文通過基于非結構網(wǎng)格的彈簧法網(wǎng)格變形和局部網(wǎng)格重構,解決含動邊界的機彈分離非定常流場數(shù)值模擬問題。在航彈與載機的分離過程中,考慮由于航彈運動引起的網(wǎng)格運動。將Navier-Stokes方程用于基于運動網(wǎng)格的流場域求解時,基本的控制方程組可表示為

    式中:?Ω(t)代表控制體積Ω(t)的邊界面;dV表示對體積的積分;ρ、u和e分別為控制體中流體密度、速度和單位體積的總能;D為網(wǎng)格運動速度,對于不運動網(wǎng)格D =0;n為?Ω(t)的外法向單位向量;P為流體靜壓張量;τ為粘性應力張量;q為由熱傳導產(chǎn)生的熱流量。

    (1)式即是笛卡爾坐標系下,邊界以任意速度運動的控制體流動控制方程組。下面將其寫為計算流體力學慣用的Navier-Stokes方程形式:

    定義守恒變量Q為

    邊界運動速度D為

    則笛卡爾坐標系下,邊界運動的控制體流動控制方程組(1)式可寫為

    式中:分量E、F、G及Ev、Fv、Gv的表達式分別為

    式中:流體速度u、υ、w及控制體邊界運動速度ug、υg、wg均為笛卡爾坐標系下的速度分量;p為流體壓力;ρ為流體密度;e為單位體積總能。

    1.2湍流模型

    選擇標準的Spalart-Allmaras湍流模型。Spalart-Allmaras模型為單方程模型,求解變量是~ν,代表除了近壁(粘性影響)區(qū)域以外的湍流運動粘性系數(shù)。~ν的輸運方程為

    式中:Gν是湍流粘性產(chǎn)生項;Yν是由于壁面阻擋與粘性阻尼引起的湍流粘性的減少;σν~和Cb2是常數(shù);μ是動力粘性系數(shù)。

    1.3網(wǎng)格和邊界條件

    采用非結構網(wǎng)格對初始流場域進行離散。在機彈分離非定常流場模擬過程中,運動的彈體表面邊界位移較小的情況下,采用彈簧法網(wǎng)格變形調整網(wǎng)格節(jié)點的位置。當邊界的位移與網(wǎng)格尺寸相比較大時,彈簧法可能導致網(wǎng)格扭曲嚴重,質量變差。因此,需要將那些扭曲率過大或超出預先指定的網(wǎng)格尺寸范圍的網(wǎng)格單元與周圍的網(wǎng)格局部聚合,把聚合后的區(qū)域當作一個新的空間域,然后在這個空間域內重新生成新的網(wǎng)格單元。

    在數(shù)值模擬機彈分離過程時,計算域外邊界設置為壓力遠場邊界條件,彈體及載機表面為無滑移絕熱壁面條件。

    1.4六自由度運動軌跡計算

    利用非定常氣動力耦合六自由度運動方程計算機彈分離過程中航彈的運動軌跡。在載機慣性坐標系下,根據(jù)牛頓定律給出剛體平動方程:

    式中:m是航彈質量;vb、χ分別是航彈質心在慣性坐標系下的速度矢量和位置矢量;F是外掛物所受合力,包括氣動力、重力等。

    在彈體坐標系下,給出剛體轉動方程。外掛物的動量矩定義為

    式中:I是彈體的慣性張量。

    根據(jù)動量矩定理

    將(8)式代入(9)式得

    1.5數(shù)值格式

    對Navier-Stokes方程的空間離散采用Roe-FDS通量格式,梯度插值采用基于節(jié)點的格林-高斯方法,流動項和修正的湍流粘性項采用2階迎風格式離散。

    2 數(shù)值驗證

    二維翼型受迫振蕩繞流場是典型的非定常流場,本節(jié)通過對NACA0012翼型的受迫振蕩非定常流場進行數(shù)值模擬,并將數(shù)值模擬結果與風洞試驗數(shù)據(jù)[18]進行對比,以此來檢驗本文計算機彈分離問題時用到的非定常流動計算方法的有效性。

    翼型攻角隨時間周期性變化的公式為

    式中:α0為初始攻角;αm為振幅;k為減縮頻率,k = ωc/2υ∞,ω為角頻率,c為弦長;t為無量綱時間,振蕩運動軸心定義為χm/ c.各參數(shù)取值如表1所示。

    表1 NACA0012翼型振蕩參數(shù)表Tab.1 Oscillation parameters of NACA 0012 airfoil

    使用與前述機彈分離問題相同的動網(wǎng)格方法和數(shù)值方法,對振蕩翼型流動進行模擬。所用網(wǎng)格和得到的法向力系數(shù)滯回曲線如圖1所示。

    圖1 網(wǎng)格和法向力系數(shù)計算結果Fig.1 Mesh and numerical result of normal force coefficient

    由圖1可見,計算得到的法向力系數(shù)滯回曲線與實驗值較為符合,也說明了計算采用的數(shù)值方法是可信的。

    3 計算外形和條件

    計算外形如圖2所示,航彈掛載于載機左側機翼內側掛架。在掛載狀態(tài)時,航彈前翼、尾翼呈“×”形,本文主要研究航彈與載機分離初期的氣動干擾特性,模擬分離后1 s的物理時間內的分離情況。

    圖2 機彈分離計算模型Fig.2 Computation model

    彈體在投放慣性坐標系Oχyz中,χ、y、z分別表示彈體質心坐標在3個方向上的分量;θχ、θy、θz分別表示彈體滾轉角、偏航角和俯仰角。載機機翼有1°的安裝角,因此彈體相對投放慣性坐標系有-1°的初始俯仰角,而偏航和滾轉角為0°.投放慣性坐標系的原點取在載機頭部頂點,彈體初始質心位置坐標為(8.322,0.079,2.572)。

    在不同初始分離條件下,對航彈投放初期與載機的分離過程進行數(shù)值模擬。計算條件:

    1)載機來流馬赫數(shù)Ma分別為0.6、0.8、1.1;

    2)載機攻角α分別為3°、5°、7°;

    3)載機側滑角β分別為-4°、-3°、3°、4°;

    4)載機飛行高度h分別為500 m、3 000 m、12 000 m;

    5)航彈初始下拋速度υb0分別為0 m/ s、5 m/ s、10 m/ s;

    6)航彈初始下拋角速度ω0分別為0°/ s、15°/ s、30°/ s.

    4 初始分離條件對分離過程影響的數(shù)值模擬結果

    4.1來流馬赫數(shù)對非定常分離的影響

    圖3和圖4分別為來流馬赫數(shù)Ma分別為0.6、0.8、1.1,來流攻角α=3°,初始下拋速度υb0=0 m/ s,初始下拋角速度ω0= 0°/ s條件下,航彈的分離軌跡及彈體姿態(tài)角隨分離時間變化的曲線圖。

    圖3 不同來流馬赫數(shù)下機彈分離軌跡Fig.3 Trajectories of bomb at different Mach numbers

    圖4 不同來流馬赫數(shù)下彈體姿態(tài)角隨時間變化曲線Fig.4 Bomb attitude angle variation over time at different Mach numbers

    從圖3和圖4可以看出,隨著來流馬赫數(shù)的增大,分離過程中彈體與載機間的氣動干擾增強,其中當來流馬赫數(shù)Ma = 1.1時,干擾最強烈,導致彈體位移和姿態(tài)角的變化,此時彈體處于跨聲速流場,受到的阻力以及載機側向洗流的影響都比較大,因此χ方向和側向位移都比較大,尤其是側向位移,隨分離時間的增加彈體迅速遠離載機機身向翼梢方向運動。在整個分離過程中,彈體始終受到抬頭力矩,彈體的俯仰角是隨著分離時間的增加而增大的,導致了升力的不斷增大。

    4.2來流攻角對非定常分離的影響

    圖5和圖6分別為來流馬赫數(shù)Ma = 0.6,來流攻角α分別為3°、5°、7°,初始下拋速度υb0=0 m/ s,初始下拋角速度ω0= 0°/ s條件下,航彈的分離軌跡及彈體姿態(tài)角隨分離時間變化的曲線圖。

    圖5 不同來流攻角下機彈分離軌跡Fig.5 Trajectories of bomb at different angles of attack

    從圖5、圖6可以看出,隨著來流攻角的增大,分離過程中彈體與載機之間的氣動干擾增強,對彈體質心位置與姿態(tài)變化的影響增大。來流攻角越大,彈體俯仰角增大越快,隨之升力增大,下落變慢。彈體在χ方向位移量隨著來流攻角的增大變化不大,但在z方向,攻角越大,彈體受到側向洗流的作用越明顯,彈體向翼梢方向運動越快。

    4.3來流側滑角對非定常分離的影響

    圖7和圖8分別為來流馬赫數(shù)Ma = 0.6,載機側滑角β分別為-4°、-3°、3°、4°,攻角α=3°,航彈初始下拋速度υb0=0 m/ s,初始下拋角速度ω0=0°/ s 時,航彈的分離軌跡及彈體姿態(tài)角隨分離時間變化的曲線圖。

    由圖7和圖8可以看出,不同側滑角下,彈體在χ方向的位移量變化不大;在y方向上,負側滑角比正向側滑角下落得快,并且負側滑角絕對值越大,下落得越快,然而正向側滑角則相反,側滑角絕對值越大,下落得越慢;在z方向上,負向側滑角下,彈體向靠近機身的方向運動,而正側滑角下,彈體向靠近翼梢的方向運動,并且運動趨勢會大很多。

    圖6 不同來流攻角下彈體姿態(tài)角隨時間變化曲線Fig.6 Bomb attitude angle variation over time at different angles of attack

    總而言之,負側滑角下,彈體受到的氣動力與氣動力矩都較為平緩,這是由于彈體受到載機側向洗流的影響,會受到側向力作用而向著翼梢方向運動,正側滑角加劇了這一作用,使彈體向著翼梢方向運動得更快;而負側滑角則抑制了側向洗流的作用,使彈體向豎直方向較快地脫離載機洗流的影響,因此一定的負側滑角有利于彈體與載機安全分離。

    4.4載機飛行高度對非定常分離的影響

    圖9和圖10分別為載機飛行高度h分別為500 m、3 000 m、12 000 m,攻角α= 5°,來流馬赫數(shù)Ma =0.8,初始下拋速度υb0=0 m/ s,初始下拋角速度ω0=0°/ s時,航彈的分離軌跡及彈體姿態(tài)角隨分離時間變化的曲線圖。

    圖7 不同來流側滑角下機彈分離軌跡Fig.7 Trajectories of bomb at different angles of sideslip

    由圖9和圖10可以看出,隨著飛行高度的增加,空氣密度降低,分離過程中彈體所受氣動力減小,載機對彈體的影響也減弱。在12 000 m高空,彈體所受的氣動力和力矩都相對較小,并且變化不劇烈,而在500 m低空,彈體所受氣動力及氣動力矩變化幅度較大,對分離軌跡和姿態(tài)有較大的影響。

    由于高空環(huán)境氣動力和力矩作用的減小,彈體在分離過程中所受氣動干擾及姿態(tài)角變化都隨飛行高度的增加而減小,阻力減小使得χ方向的位移減小;y方向重力的作用遠大于升力的作用,使得彈體下落得更快;側向洗流影響減小也使z方向的位移變化不大??傊?在較高的飛行高度下,機彈分離過程中載機與彈體間干擾作用較小,彈體能較快地脫離載機流場,有利于分離。

    4.5初始下拋速度對非定常分離的影響

    圖11和圖12分別為來流馬赫數(shù)Ma = 0.6,載機攻角α=5°,航彈初始下拋速度υb0分別為0 m/ s、5 m/ s、10 m/ s,初始下拋角速度ω0= 0°/ s時,航彈的分離軌跡及彈體姿態(tài)角隨分離時間變化的曲線圖。

    由圖11和圖12可以看出,隨著初始下拋速度的增加,機彈分離過程中機彈間的氣動干擾并沒有減弱,對彈體氣動特性影響仍然較大,但由于初始速度較大,彈體能迅速通過載機的下洗流場,載機對彈體影響時間縮短,使彈體質心位置與姿態(tài)角還未發(fā)生較大變化時與載機安全分離。彈體質心在χ方向,隨初始下拋速度增大的影響并不大;y方向上,隨著初始下拋速度的增加,相同時間內彈體在豎直方向上的位移越大,迅速遠離載機;z方向上,彈體初始下拋速度越大,位移量越小。因此,一定的初始下拋速度有利于機彈安全分離。

    圖8 不同來流側滑角下彈體姿態(tài)角隨時間變化曲線Fig.8 Bomb attitude angle variation over time at different angles of sideslip

    4.6初始下拋角速度對非定常分離的影響

    圖13和圖14分別為來流馬赫數(shù)Ma = 0.6,載機攻角α=5°,航彈初始下拋速度υb0=0 m/ s,初始下拋角速度ω0分別為0°/ s、15°/ s、30°/ s時,航彈的分離軌跡及彈體姿態(tài)角隨分離時間變化的曲線圖。

    圖9 不同飛行高度下機彈分離軌跡Fig.9 Trajectories of bomb at different flight altitudes

    圖10 不同飛行高度下彈體姿態(tài)角隨時間變化曲線Fig.10 Bomb attitude angle variation over time at different flight altitudes

    由圖13和圖14可以看出,加初始下拋角速度后,彈體迅速低頭,出現(xiàn)了負攻角,彈體所受升力為負,由于該彈是不穩(wěn)定的,彈體受到低頭的俯仰力矩,導致負攻角增大,所受負升力變大,彈體下落也越快。同時,隨著初始下拋角速度的增大,彈體負攻角增大,因此受到的阻力也隨之迅速增大,彈體質心在χ方向的位移也就越大。由于載機側向洗流對彈體的作用并沒有較大變化,因此彈體所受側向力和彈體質心在z方向的位移沒有明顯變化,彈體仍向翼梢方向運動。

    圖11 不同初始下拋速度的機彈分離軌跡Fig.11 Trajectories of bomb at different initial separation speeds

    但是過大的初始下拋角速度也使得彈體的負攻角迅速增大,并可能導致彈體的翻轉,因此適當?shù)某跏枷聮伣撬俣扔欣诎踩蛛x。

    5 航彈與載機分離安全性分析

    對本文所計算的各個工況,在機彈分離過程中彈體質心在χ方向的位移并不大,沒有威脅到載機的安全,所以本節(jié)主要關注其在y方向和z方向的位移,以及彈體在Oyz平面內的軌跡,以判斷機彈之間的分離是否安全。

    圖15~圖17分別為來流馬赫數(shù)Ma分別為0.6、0.8、1.1條件下,載機在不同攻角下投放航彈時,航彈的分離軌跡。

    由圖15可見,馬赫數(shù)Ma = 0.6時,各攻角下航彈都能與載機安全分離,但隨著攻角的增大,彈體的側向位移增大,彈體下落得更慢。由圖16可見,在Ma =0.8,α=5°時,彈體所受側向力和升力都比較大,彈體幾乎水平向翼梢運動,并有向上靠近機翼的趨勢,對載機十分危險。繼續(xù)增大攻角或來流馬赫數(shù),彈體將直接撞上掛架,見圖16(c)和圖17(b).

    圖12 不同初始下拋速度的彈體姿態(tài)角隨時間變化曲線Fig.12 Bomb attitude angle variation over time at different initial separation speeds

    通過數(shù)值模擬得到了航彈與載機安全分離的投放包絡,表2和表3是根據(jù)數(shù)值模擬結果進行整理后給出的航彈與載機分離情況表。

    表2 不同馬赫數(shù)、攻角下分離安全情況表Tab.2 Separation safety situation at different Mach numbers and angles of attack

    圖13 不同初始下拋角速度的機彈分離軌跡Fig.13 Trajectories of bomb at different initial angular velocities

    表3 不同馬赫數(shù)、攻角、側滑角下分離安全情況表Tab.3 Separation safety situation at different Mach numbers,angles of attack and angles of sideslip

    6 結論

    從各種初始分離條件下,航彈與載機分離過程的數(shù)值模擬與分析結果,可得出以下結論:

    1)隨著來流馬赫數(shù)的增加,分離過程中載機對航彈的氣動干擾增強,航彈向下與載機的分離變慢,但沿機翼翼展方向的運動變得更快。

    2)隨著載機攻角增大,分離過程中航彈受到載機的氣動干擾增強,對彈體質心位置與姿態(tài)變化影響增大。升力和側向力的增大導致航彈向下與載機的分離減慢,向側向有顯著運動,投放時載機的攻角越大對安全分離越不利。

    3)當載機的側滑角為正時,位于左側機翼下的航彈向機翼翼梢方向的運動加快,向下分離運動減慢,不利于安全分離;當載機側滑角為負時,有利于左側機翼下的航彈的安全分離。

    圖14 不同初始下拋角速度的彈體姿態(tài)角隨時間變化曲線Fig.14 Bomb attitude angle variation over time at different initial angular velocities

    圖15 Ma =0.6與不同攻角時彈體的分離軌跡圖Fig.15 Trajectories of bomb at Ma =0.6 and different angles of attack

    4)飛行高度增加,分離過程中彈體所受的氣動力減小,受到載機的干擾作用較小,彈體能較快脫離載機所處流場,有利于安全分離。

    圖16 Ma =0.8與不同攻角時彈體的分離軌跡圖Fig.16 Trajectories of bomb at Ma =0.8 and different angles of attack

    圖17 Ma =1.1與不同攻角時彈體的分離軌跡圖Fig.17 Trajectories of bomb at Ma =1.1 and different angles of attack

    5)初始下拋速度增大,可使彈體迅速通過載機的下洗流場,有利于航彈與載機的安全分離;初始下拋角速度可改變分離過程中彈體姿態(tài)和所受氣動力,適當?shù)某跏枷聮伣撬俣扔欣诎踩蛛x。

    參考文獻(References)

    [1]趙忠良,龍堯松,余立,等.高超聲速風洞子母彈分離干擾測力試驗技術[J].流體力學實驗與測量,2004,18(3):32-35.ZHAO Zhong-liang, LONG Yao-song, YU Li, et al.Hypersonic wind tunnel interference force test of sub-munition reentry separation [J].Experiments and Measurements in Fluid Mechanics, 2004,18(3):32-35.(in Chinese)

    [2]雷娟棉,吳甲生,肖雅彬.布撒器-子彈氣動干擾風洞實驗研究[J].兵工學報,2005,26(4):535-539.LEI Juan-mian, WU Jia-sheng, XIAO Ya-bin.Aerodynamic interference wind tunnel experiment investigation for dispenser and submunition [J].Acta Armamentarii, 2005,26(4):535-539.(in Chinese)

    [3]達興亞,周為群,趙忠良,等.窄條翼導彈模型搖滾運動動力學特性研究[J].空氣動力學學報, 2014, 32(2):154-158.DA Xing-ya, ZHOU Wei-qun, ZHAO Zhong-liang, et al.Research on the kinetics of rock motion for a missile model with strake wings [ J].Acta Aerodynamica Sinica, 2014, 32 (2): 154-158.(in Chinese)

    [4]孫新利,蔡星會,王少龍,等.子母彈靜態(tài)開艙拋射實驗[J].兵工學報,2002,23(2):258-260.SUN Xin-li, CAI Xing-hui, WANG Shao-long, et al.Projectile tests from a bomb dispenser [J].Acta Armamentarii, 2002,23(2):258-260.(in Chinese)

    [5]唐志共,李彬,鄭鳴,等.飛行器外掛投放數(shù)值模擬[J].空氣動力學學報,2009,27(5):592-596.TANG Zhi-gong, LI Bin, ZHENG Ming, et al.Store separation simulation using overset unstructured grid [J].Acta Aerodynamica Sinica, 2009,27(5):592-596.(in Chinese)

    [6]張玉東,紀楚群.子母彈分離過程的數(shù)值模擬方法[J].空氣動力學學報,2003,21(1):47-52.ZHANG Yu-dong, JI Chu-qun.The numerical simulation of submunition separation processes from dispenser [J].Acta Aerodynamica Sinica, 2003,21(1):47-52.(in Chinese)

    [7]李亭鶴,閻超,李躍軍.重疊網(wǎng)格技術中割補法的研究與改進[J].北京航空航天大學學報,2005,31(4):402-406.LI Ting-he, YAN Chao, LI Yue-jun.Investigation and enhancement of cut-paste algorithm in overlapping grid [J].Journal of Beijing University of Aeronautics and Astronautics, 2005,31(4): 402-406.(in Chinese)

    [8]雷娟棉,苗瑞生,居賢銘.戰(zhàn)術火箭子母戰(zhàn)斗部第一次拋撒分離多體干擾流場數(shù)值模擬[J].北京理工大學學報,2004, 24(9):766-769,785.LEI Juan-mian, MIAO Rui-sheng, JU Xian-ming.Numerical investigation of multi-body flow fields in the first separation of tactical rocket's submunitions [J].Transactions of Beijing Institute of Technology, 2004,24(9):766-769,785.(in Chinese)

    [9]張軍,譚俊杰,耿繼輝.含運動邊界二維多體干擾流場的數(shù)值模擬[J].空氣動力學學報,2003,21(4):449-453.ZHANG Jun, TAN Jun-jie, GENG Ji-hui.Numerical simulation of 2D multi-bodies with moving boundaries [J].Acta Aerodynamica Sinica, 2003,21(4):449-453.(in Chinese)

    [10]桑為民,李鳳蔚,鄂秦.應用笛卡爾非結構切割網(wǎng)格進行外掛物投放的數(shù)值模擬[J].計算物理,2002,19(3):268-272.SANG Wei-min, LI Feng-wei, E Qin.Numerical simulation of the store separation using unstructured Cartesian cutted grid [J].Chinese Journal of Computational Physics, 2002,19(3):268-272.(in Chinese)

    [11]陶如意,王浩,趙潤祥,等.超音速子母彈分離激波干擾特性研究[J].兵工學報, 2011, 32(10):1206-1211.TAO Ru-yi, WANG Hao, ZHAO Run-xiang,et al Research on shock/ shock wave disturbance characteristics in separation of supersonic cluster munition[J].Acta Armamentarii, 2011, 32(10): 1206-1211.(in Chinese)

    [12]Kim D H, Choi J H, Kwon O J.Detached eddy simulation of weapons bay flows and store separation[J].Computers & Fluids, 2015, 121:1-10.

    [13]Berglind T, Tysell L.Numerical investigation of the impact of maneuver on store separation trajectories[C]∥28th AIAA Applied Aerodynanics, Conference.Chicago, Illinois: AIAA, 2010: 4241.

    [14]Purdon M.F-35 pre-flight store separation analyses: innovative techniques for affordability[C]∥47th AIAA Aerospace Sciences Meeting.Oriando, Florida: AIAA, 2009: 102.

    [15]Flora T J, Reeder M F, Lofthouse A, et al.Dynamic store release of ice models from a cavity into mach 2.9 flow[J].Journal of Aircraft, 2014, 51(6):1927-1941.

    [16]Baker W, Keen S, Morgret C, et al.Validation of weapon separation predictions using F/ A-22 flight test results[C]∥USAF Developmental Test and Evaluation Summit.Woodland Hills, California , US: Air Force, 2004.

    [17]Johnson R A, Stanek M J, Grove J E.Store separation trajectory deviations due to unsteady weapons bay aerodynamics[C]∥46th AIAA Aerospace Sciences Meeting and Exhibit.Reno, NV, US: AIAA,2008.

    [18]Landon R H.NACA 0012 oscillating and transient pitching,compendium of unsteady aerodynamic measurement, AGARD report 702[R].Batina, JT, US: AGARD, 1982.

    Numerical Simulation about the Effect of Initial Separation Condition on Safety of Aerial Bomb Separated from an Aircraft

    LEI Juan-mian, NIU Jian-ping, WANG Suo-zhu, ZHOU Qi
    (School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China)

    Abstract:In order to investigate the effect of initial release conditions on the safety of aerial bomb separated from aircraft, the separation process of aerial bomb from aircraft is simulated by using unsteady CFD numerical simulation method, dynamic mesh technique and 6-DOF trajectory equations.The trajectory and attitude variation of bomb after release are given under different conditions of Mach number, angle of attack, angle of sideslip, flight altitude, initial separation velocity and initial angular velocity.The results show that the aircraft has a strong aerodynamic interference on the bomb when separating, and has a great influence on the bomb's aerodynamic characteristics, trajectory and attitude.As the Mach number and angle of attack increase, the interference of aircraft on the bomb enhances, and the safety gets worse.For the bomb mounted under the left wing, a certain negative angle of sideslip is favorable to safe separation.The higher the flight altitude is, the safer the separation is.The certain initial separation speed and proper angular velocity are benefit for separation safety.

    Key words:ordnance science and technology; weapon separation; numerical simulation; dynamic mesh; aerodynamic characteristics

    作者簡介:雷娟棉(1968—),女,教授,博士生導師。E-mail: leijm@ bit.edu.cn;王鎖柱(1984—),男,博士后。E-mail: wangsuozhu@163.com

    收稿日期:2015-05-08

    DOI:10.3969/ j.issn.1000-1093.2016.02.023

    中圖分類號:V211.3

    文獻標志碼:A

    文章編號:1000-1093(2016)02-0357-10

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學中的應用
    科教導刊(2016年28期)2016-12-12 06:22:00
    基于氣象信息及風場信息的風機輪轂處風速預測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預處理和加工研究報告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    午夜免费成人在线视频| 久99久视频精品免费| 夜夜夜夜夜久久久久| 观看免费一级毛片| 成人特级黄色片久久久久久久| avwww免费| 久久中文字幕一级| 97超级碰碰碰精品色视频在线观看| 啦啦啦观看免费观看视频高清| 国产黄片美女视频| 九九久久精品国产亚洲av麻豆 | a级毛片在线看网站| www国产在线视频色| 亚洲中文日韩欧美视频| 神马国产精品三级电影在线观看| 一个人免费在线观看的高清视频| 久99久视频精品免费| 亚洲欧美日韩高清在线视频| 麻豆成人av在线观看| 在线视频色国产色| 国产高清有码在线观看视频| 校园春色视频在线观看| 欧美日韩精品网址| 亚洲激情在线av| 99久久无色码亚洲精品果冻| 精品久久久久久久久久久久久| 国内精品一区二区在线观看| 草草在线视频免费看| 性色av乱码一区二区三区2| 国产日本99.免费观看| 99久久精品一区二区三区| 全区人妻精品视频| 欧美午夜高清在线| 香蕉国产在线看| 两个人看的免费小视频| 1000部很黄的大片| 久久香蕉国产精品| 99久久国产精品久久久| 亚洲一区二区三区不卡视频| 亚洲美女黄片视频| 精品国产乱码久久久久久男人| 免费在线观看亚洲国产| 日韩国内少妇激情av| 亚洲精品一区av在线观看| 午夜福利免费观看在线| 亚洲一区高清亚洲精品| 国产精品国产高清国产av| 搡老岳熟女国产| 香蕉丝袜av| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 国产真人三级小视频在线观看| 欧美日韩乱码在线| 法律面前人人平等表现在哪些方面| 亚洲国产精品合色在线| 欧美午夜高清在线| 别揉我奶头~嗯~啊~动态视频| 精品国产乱子伦一区二区三区| 国内精品美女久久久久久| 97超视频在线观看视频| 日韩有码中文字幕| 一级毛片精品| 999久久久国产精品视频| 亚洲精品粉嫩美女一区| 桃色一区二区三区在线观看| 国产亚洲精品av在线| 久久久久亚洲av毛片大全| 九九热线精品视视频播放| 可以在线观看的亚洲视频| 国产精品电影一区二区三区| 无遮挡黄片免费观看| 国产伦人伦偷精品视频| 欧美3d第一页| 国产黄片美女视频| АⅤ资源中文在线天堂| 在线看三级毛片| netflix在线观看网站| 神马国产精品三级电影在线观看| 在线观看美女被高潮喷水网站 | 精品国产乱码久久久久久男人| 性欧美人与动物交配| 欧美乱色亚洲激情| 国产精品久久久久久精品电影| 欧美成人一区二区免费高清观看 | 国产一区二区在线av高清观看| 午夜精品久久久久久毛片777| 脱女人内裤的视频| 2021天堂中文幕一二区在线观| 国内毛片毛片毛片毛片毛片| www.精华液| 变态另类成人亚洲欧美熟女| 九色成人免费人妻av| or卡值多少钱| 成人三级做爰电影| av天堂中文字幕网| 国产成人福利小说| 老司机午夜福利在线观看视频| 日本熟妇午夜| 欧美黄色片欧美黄色片| 美女 人体艺术 gogo| 亚洲av片天天在线观看| 国产又色又爽无遮挡免费看| 午夜免费激情av| 天堂网av新在线| 老熟妇仑乱视频hdxx| or卡值多少钱| 欧美国产日韩亚洲一区| 一进一出抽搐gif免费好疼| 亚洲精品国产精品久久久不卡| 亚洲人与动物交配视频| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 日本黄色片子视频| 亚洲成人免费电影在线观看| 听说在线观看完整版免费高清| 最新在线观看一区二区三区| 免费无遮挡裸体视频| 国产蜜桃级精品一区二区三区| www.www免费av| 男女那种视频在线观看| 久久久久久人人人人人| 99国产精品一区二区蜜桃av| 国产探花在线观看一区二区| 一进一出抽搐gif免费好疼| 亚洲狠狠婷婷综合久久图片| 色综合欧美亚洲国产小说| 网址你懂的国产日韩在线| 久久天堂一区二区三区四区| 99在线视频只有这里精品首页| 久久久久久久久中文| 天天躁狠狠躁夜夜躁狠狠躁| 日韩高清综合在线| 9191精品国产免费久久| 国产精品一及| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 99精品欧美一区二区三区四区| 久99久视频精品免费| 国产伦在线观看视频一区| 巨乳人妻的诱惑在线观看| 色精品久久人妻99蜜桃| 看片在线看免费视频| 婷婷精品国产亚洲av在线| 黄频高清免费视频| 精品一区二区三区av网在线观看| 亚洲欧美日韩无卡精品| 午夜免费激情av| 99久久精品一区二区三区| 亚洲乱码一区二区免费版| 日日干狠狠操夜夜爽| 99久久成人亚洲精品观看| 国产一级毛片七仙女欲春2| 久久精品综合一区二区三区| 人妻久久中文字幕网| 亚洲人成网站在线播放欧美日韩| 欧美一区二区国产精品久久精品| 欧美中文综合在线视频| 91字幕亚洲| 最近最新中文字幕大全免费视频| 淫妇啪啪啪对白视频| 欧美又色又爽又黄视频| 日本 欧美在线| 精品国产乱子伦一区二区三区| 日韩欧美国产一区二区入口| 嫩草影院入口| 少妇的逼水好多| 97人妻精品一区二区三区麻豆| 午夜福利在线观看吧| 亚洲在线自拍视频| 99re在线观看精品视频| 亚洲五月婷婷丁香| 一本综合久久免费| 88av欧美| 亚洲avbb在线观看| 亚洲 国产 在线| 69av精品久久久久久| 国产在线精品亚洲第一网站| 欧美在线黄色| 亚洲人成电影免费在线| 网址你懂的国产日韩在线| 老汉色av国产亚洲站长工具| 精品99又大又爽又粗少妇毛片 | 国产成人欧美在线观看| 女同久久另类99精品国产91| 在线观看美女被高潮喷水网站 | 午夜精品在线福利| 亚洲国产精品合色在线| 99在线人妻在线中文字幕| 精品久久蜜臀av无| 日韩人妻高清精品专区| 高清毛片免费观看视频网站| 亚洲欧美日韩东京热| 丰满人妻一区二区三区视频av | 欧美中文综合在线视频| 亚洲真实伦在线观看| 欧美zozozo另类| 亚洲在线自拍视频| 热99在线观看视频| 欧美黄色淫秽网站| 麻豆久久精品国产亚洲av| 91av网站免费观看| 国产精品野战在线观看| 久久中文看片网| 国内精品久久久久久久电影| 久久精品91蜜桃| 国产一区二区三区在线臀色熟女| 黑人欧美特级aaaaaa片| 国产一区二区在线观看日韩 | 色视频www国产| 午夜两性在线视频| 亚洲aⅴ乱码一区二区在线播放| 一级毛片女人18水好多| 久久这里只有精品中国| 亚洲人与动物交配视频| 黑人巨大精品欧美一区二区mp4| 国产一区二区在线av高清观看| 一卡2卡三卡四卡精品乱码亚洲| 三级国产精品欧美在线观看 | 免费av不卡在线播放| 国产精品爽爽va在线观看网站| 色尼玛亚洲综合影院| 国产精品影院久久| 久久午夜亚洲精品久久| 一个人看的www免费观看视频| 国产午夜精品久久久久久| 亚洲av电影不卡..在线观看| av中文乱码字幕在线| 久久精品91无色码中文字幕| 在线国产一区二区在线| 久久久久久久久久黄片| 美女扒开内裤让男人捅视频| 欧美日本亚洲视频在线播放| 精品熟女少妇八av免费久了| 国产亚洲精品久久久久久毛片| 亚洲狠狠婷婷综合久久图片| 日日夜夜操网爽| 在线观看一区二区三区| 国产视频内射| 日本免费一区二区三区高清不卡| 嫁个100分男人电影在线观看| 亚洲专区字幕在线| 哪里可以看免费的av片| 成人18禁在线播放| 国产又色又爽无遮挡免费看| 欧美3d第一页| 欧美又色又爽又黄视频| 精品熟女少妇八av免费久了| 亚洲成人免费电影在线观看| a级毛片在线看网站| 亚洲国产日韩欧美精品在线观看 | 免费观看的影片在线观看| 国产伦在线观看视频一区| 天堂动漫精品| 中文字幕最新亚洲高清| 亚洲欧美日韩无卡精品| 男人的好看免费观看在线视频| 熟妇人妻久久中文字幕3abv| 国产精品爽爽va在线观看网站| 最新美女视频免费是黄的| 色在线成人网| 国产真实乱freesex| h日本视频在线播放| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧洲精品一区二区精品久久久| 村上凉子中文字幕在线| 国产美女午夜福利| 亚洲一区二区三区色噜噜| 国产精品久久久久久人妻精品电影| 小蜜桃在线观看免费完整版高清| 一二三四社区在线视频社区8| 午夜福利免费观看在线| 天堂影院成人在线观看| 久久久国产欧美日韩av| x7x7x7水蜜桃| 大型黄色视频在线免费观看| 精品欧美国产一区二区三| 婷婷亚洲欧美| 最近最新中文字幕大全免费视频| 国产精品99久久久久久久久| 露出奶头的视频| 淫妇啪啪啪对白视频| 日本 av在线| 欧美丝袜亚洲另类 | 亚洲人成伊人成综合网2020| 老熟妇仑乱视频hdxx| 一本久久中文字幕| 黑人巨大精品欧美一区二区mp4| 日本a在线网址| 一级毛片高清免费大全| 亚洲中文av在线| av黄色大香蕉| 精品国产美女av久久久久小说| 少妇的丰满在线观看| 日本免费一区二区三区高清不卡| 久久久国产成人免费| 99视频精品全部免费 在线 | 国内少妇人妻偷人精品xxx网站 | 啦啦啦免费观看视频1| 久久久久性生活片| 男人舔女人的私密视频| 成年人黄色毛片网站| av天堂中文字幕网| 久久亚洲真实| 欧美日韩精品网址| 一个人观看的视频www高清免费观看 | 亚洲国产精品久久男人天堂| 男人舔女人下体高潮全视频| 一卡2卡三卡四卡精品乱码亚洲| 国产野战对白在线观看| 99久久精品国产亚洲精品| 国产高清视频在线播放一区| 国产真实乱freesex| 一区二区三区激情视频| 十八禁人妻一区二区| 亚洲午夜精品一区,二区,三区| 香蕉丝袜av| 国产亚洲精品综合一区在线观看| 亚洲无线在线观看| 真人做人爱边吃奶动态| 天堂√8在线中文| 久久午夜综合久久蜜桃| 18禁国产床啪视频网站| 久久热在线av| 88av欧美| 啦啦啦观看免费观看视频高清| 91麻豆精品激情在线观看国产| 欧美中文综合在线视频| 99久久精品热视频| 国产成人啪精品午夜网站| 一级a爱片免费观看的视频| 夜夜夜夜夜久久久久| 嫩草影视91久久| 欧美日韩福利视频一区二区| 大型黄色视频在线免费观看| 中文字幕精品亚洲无线码一区| 亚洲国产高清在线一区二区三| 在线观看午夜福利视频| 小说图片视频综合网站| 男人的好看免费观看在线视频| 一二三四在线观看免费中文在| 国产成人精品无人区| 无限看片的www在线观看| 极品教师在线免费播放| 国产伦在线观看视频一区| 久久久精品欧美日韩精品| 亚洲av中文字字幕乱码综合| 女生性感内裤真人,穿戴方法视频| 高潮久久久久久久久久久不卡| 1024香蕉在线观看| xxx96com| www.999成人在线观看| 国产人伦9x9x在线观看| 亚洲av免费在线观看| 真实男女啪啪啪动态图| 深夜精品福利| 欧美最黄视频在线播放免费| 久久久久久久午夜电影| 亚洲中文日韩欧美视频| 久久中文字幕人妻熟女| 老司机深夜福利视频在线观看| 成人特级av手机在线观看| 亚洲 欧美一区二区三区| 长腿黑丝高跟| 免费一级毛片在线播放高清视频| 综合色av麻豆| 亚洲国产日韩欧美精品在线观看 | 国产69精品久久久久777片 | 日本一本二区三区精品| 免费高清视频大片| 麻豆一二三区av精品| 国产亚洲欧美在线一区二区| 小蜜桃在线观看免费完整版高清| 99国产综合亚洲精品| 亚洲无线观看免费| 观看美女的网站| 特级一级黄色大片| 欧美激情在线99| 欧美不卡视频在线免费观看| a级毛片在线看网站| 亚洲专区国产一区二区| 嫁个100分男人电影在线观看| 成人特级黄色片久久久久久久| 男插女下体视频免费在线播放| 亚洲电影在线观看av| 美女cb高潮喷水在线观看 | e午夜精品久久久久久久| 最好的美女福利视频网| 精品日产1卡2卡| 亚洲国产精品sss在线观看| 久久草成人影院| 国产一区在线观看成人免费| 亚洲国产中文字幕在线视频| 久久精品国产清高在天天线| 精品国产乱子伦一区二区三区| 国产av不卡久久| 窝窝影院91人妻| 久久久久久久午夜电影| 国产亚洲精品av在线| 精品国产乱码久久久久久男人| 99国产精品一区二区三区| 91在线精品国自产拍蜜月 | 久久精品夜夜夜夜夜久久蜜豆| 国产成人啪精品午夜网站| 色综合亚洲欧美另类图片| 99久久成人亚洲精品观看| 色吧在线观看| 看免费av毛片| 天天添夜夜摸| 日韩精品中文字幕看吧| 99riav亚洲国产免费| 成人特级黄色片久久久久久久| 久久午夜亚洲精品久久| www.自偷自拍.com| 国产一区二区在线观看日韩 | 国产精品久久电影中文字幕| 色吧在线观看| 日日摸夜夜添夜夜添小说| a级毛片a级免费在线| 国产伦在线观看视频一区| 欧美成人性av电影在线观看| 国产精品一区二区精品视频观看| a在线观看视频网站| 亚洲电影在线观看av| 麻豆国产av国片精品| 国产高清视频在线播放一区| 桃色一区二区三区在线观看| avwww免费| 久久午夜综合久久蜜桃| 国产三级黄色录像| 老司机深夜福利视频在线观看| 欧美在线一区亚洲| 久久国产精品影院| 亚洲五月婷婷丁香| 又粗又爽又猛毛片免费看| 最新在线观看一区二区三区| 精品无人区乱码1区二区| 在线观看66精品国产| 亚洲五月婷婷丁香| 国产三级中文精品| 夜夜夜夜夜久久久久| 国产一区二区在线观看日韩 | 中文字幕高清在线视频| 91av网站免费观看| 一级黄色大片毛片| 欧美日韩一级在线毛片| 天天躁狠狠躁夜夜躁狠狠躁| 99久久成人亚洲精品观看| 中文字幕人妻丝袜一区二区| 国产成人精品久久二区二区免费| 精品免费久久久久久久清纯| 九九在线视频观看精品| 国产真实乱freesex| 精品人妻1区二区| 欧美成人免费av一区二区三区| 国产高清激情床上av| 亚洲av成人精品一区久久| 亚洲av日韩精品久久久久久密| 亚洲专区字幕在线| www.www免费av| 日韩欧美 国产精品| 精品无人区乱码1区二区| av国产免费在线观看| 国产精品 欧美亚洲| 国产高清视频在线观看网站| 黑人操中国人逼视频| 69av精品久久久久久| 夜夜看夜夜爽夜夜摸| 一进一出抽搐动态| 午夜福利欧美成人| 国内久久婷婷六月综合欲色啪| 亚洲精品美女久久久久99蜜臀| 成人精品一区二区免费| 成人午夜高清在线视频| 国产极品精品免费视频能看的| 午夜激情福利司机影院| 18禁黄网站禁片免费观看直播| 久久久久精品国产欧美久久久| 国产精品一区二区三区四区久久| 香蕉久久夜色| 又大又爽又粗| 91字幕亚洲| 99久久综合精品五月天人人| 亚洲av免费在线观看| 亚洲国产欧美网| 男女之事视频高清在线观看| 女同久久另类99精品国产91| 中文字幕精品亚洲无线码一区| 国产又色又爽无遮挡免费看| 精品国产乱子伦一区二区三区| 极品教师在线免费播放| 黑人巨大精品欧美一区二区mp4| 国产成人aa在线观看| 精华霜和精华液先用哪个| 欧美大码av| 九九热线精品视视频播放| 国产精品国产高清国产av| 国产亚洲av嫩草精品影院| 精品人妻1区二区| 国产精品 欧美亚洲| 国产午夜精品久久久久久| 不卡av一区二区三区| 午夜免费激情av| 99国产极品粉嫩在线观看| 国产久久久一区二区三区| 亚洲熟妇熟女久久| 不卡av一区二区三区| 国产亚洲精品久久久com| www.自偷自拍.com| 欧美国产日韩亚洲一区| 欧美一区二区国产精品久久精品| 亚洲色图av天堂| 亚洲色图 男人天堂 中文字幕| 日本五十路高清| 伦理电影免费视频| 欧美性猛交黑人性爽| 欧美日韩精品网址| 亚洲人与动物交配视频| 草草在线视频免费看| 日日摸夜夜添夜夜添小说| 亚洲av熟女| 悠悠久久av| 特大巨黑吊av在线直播| 国产av不卡久久| 欧美精品啪啪一区二区三区| 好男人在线观看高清免费视频| 亚洲国产精品999在线| 国产综合懂色| 免费看美女性在线毛片视频| 国产亚洲精品久久久久久毛片| 国内精品久久久久久久电影| 国产欧美日韩一区二区三| 看黄色毛片网站| 夜夜爽天天搞| 欧美+亚洲+日韩+国产| 日本在线视频免费播放| 国产精品爽爽va在线观看网站| 一二三四在线观看免费中文在| 热99re8久久精品国产| 男女做爰动态图高潮gif福利片| 亚洲 欧美一区二区三区| 欧美zozozo另类| 老熟妇仑乱视频hdxx| 国产精品久久电影中文字幕| 成人特级av手机在线观看| 一个人看的www免费观看视频| bbb黄色大片| 黄色女人牲交| 亚洲国产日韩欧美精品在线观看 | 久久久久久久久免费视频了| 在线观看美女被高潮喷水网站 | 脱女人内裤的视频| 成人午夜高清在线视频| 日韩免费av在线播放| 久久伊人香网站| 精品久久久久久久毛片微露脸| 国产视频一区二区在线看| 成年免费大片在线观看| 88av欧美| 亚洲国产精品合色在线| 国产成人aa在线观看| 国产伦精品一区二区三区视频9 | 国内精品一区二区在线观看| 亚洲人成电影免费在线| 九九久久精品国产亚洲av麻豆 | 变态另类丝袜制服| 国产黄片美女视频| aaaaa片日本免费| 日韩欧美三级三区| 757午夜福利合集在线观看| 国产真实乱freesex| 热99在线观看视频| 亚洲国产欧美网| 日本一二三区视频观看| 搡老岳熟女国产| 国产欧美日韩精品一区二区| 国产精品 国内视频| 波多野结衣高清无吗| 国产人伦9x9x在线观看| 黄频高清免费视频| 欧美一区二区精品小视频在线| 成人18禁在线播放| 国内毛片毛片毛片毛片毛片| 色哟哟哟哟哟哟| 蜜桃久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 动漫黄色视频在线观看| 午夜福利成人在线免费观看| 90打野战视频偷拍视频| 成年女人看的毛片在线观看| 人妻久久中文字幕网| 90打野战视频偷拍视频| 亚洲 欧美 日韩 在线 免费| 午夜福利成人在线免费观看| 露出奶头的视频| 亚洲精品在线美女| 国产激情久久老熟女| 欧美性猛交╳xxx乱大交人| 日本黄大片高清| 日韩欧美在线二视频| 久久中文字幕一级| 国产高清视频在线播放一区| 久久久久精品国产欧美久久久| 制服丝袜大香蕉在线| 久久久久久人人人人人| 国产人伦9x9x在线观看| 91在线精品国自产拍蜜月 | 国产精品精品国产色婷婷| 欧美日韩黄片免| 国产精品 欧美亚洲| 精品国产乱码久久久久久男人| 老熟妇仑乱视频hdxx| 国产又黄又爽又无遮挡在线| 国产成人精品久久二区二区91| 久久久久亚洲av毛片大全| 怎么达到女性高潮|