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

    水下氣動(dòng)發(fā)射裝置發(fā)射性能參數(shù)靈敏度分析

    2015-06-05 09:00:42牛清勇李天勻朱翔朱顯明
    關(guān)鍵詞:發(fā)射管發(fā)射裝置氣瓶

    牛清勇,李天勻,朱翔,朱顯明

    (1.華中科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074;2.華中科技大學(xué)船舶和海洋水動(dòng)力湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢430074;3.中國艦船研究設(shè)計(jì)中心船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室,湖北武漢430064)

    作為一種重要的水聲對抗器材,某型水下氣動(dòng)發(fā)射裝置,利用高壓氣瓶、發(fā)射閥和發(fā)射管等組件將武器以一定速度發(fā)射出管。建立描述該裝置發(fā)射過程的數(shù)學(xué)模型并對內(nèi)彈道性能進(jìn)行參數(shù)影響規(guī)律和靈敏度分析,對該發(fā)射系統(tǒng)的優(yōu)化設(shè)計(jì)和發(fā)射試驗(yàn)的參數(shù)估計(jì)等均有重要意義。由于保密原因該型發(fā)射裝置的可收集資料較少。因其與魚雷或水雷發(fā)射裝置在部分構(gòu)件上原理相似,研究該發(fā)射裝置時(shí)可參考其他發(fā)射系統(tǒng)的分析方法。針對魚雷或水雷的發(fā)射系統(tǒng),已經(jīng)在內(nèi)彈道數(shù)學(xué)建模和內(nèi)流場仿真分析方面開展有許多研究[1-5]。在建立了內(nèi)彈道模型的基礎(chǔ)上,內(nèi)彈道性能的參數(shù)影響規(guī)律的研究也隨之開展起來。對渦輪泵式魚雷發(fā)射裝置[6-7]、水雷發(fā)射裝置[8]和液壓平衡式魚雷發(fā)射裝置[9]內(nèi)彈道性能的參數(shù)影響規(guī)律的研究表明,發(fā)射氣瓶初壓、發(fā)射時(shí)潛艇深度和速度以及發(fā)射閥的開啟規(guī)律等對內(nèi)彈道性能均有影響。發(fā)射管內(nèi)武器的移動(dòng)速度、發(fā)射管內(nèi)膛壓既隨發(fā)射氣瓶初始壓強(qiáng)的增加而增大,又隨發(fā)射閥開啟面積的增大而增大。發(fā)射深度和艇速等對內(nèi)彈道性能也有一定的影響。已有文獻(xiàn)在分析模型參數(shù)對內(nèi)彈道性能的影響時(shí),只討論了參數(shù)本身的變化對內(nèi)彈道性能的影響規(guī)律,并未針對不同參數(shù)之間的影響程度即靈敏度進(jìn)行比較。

    本文以某型水下氣動(dòng)發(fā)射裝置為對象,在建立了描述發(fā)射過程數(shù)學(xué)模型的基礎(chǔ)上,利用局部靈敏度分析方法對內(nèi)彈道性能的主要參數(shù)進(jìn)行了影響規(guī)律分析,并對不同參數(shù)之間的靈敏度進(jìn)行了比較,確定了武器出管速度和最大膛壓的主要影響參數(shù)。

    1 發(fā)射系統(tǒng)數(shù)學(xué)模型

    圖1為某型水下氣動(dòng)發(fā)射裝置的結(jié)構(gòu)簡圖。發(fā)射閥開啟后,高壓空氣被引入發(fā)射管,武器被推動(dòng)并沖破端蓋處金屬密封片以一定速度出管。利用熱力學(xué)和動(dòng)力學(xué)知識(shí),可建立描述該發(fā)射過程的數(shù)學(xué)模型。

    圖1 某水下氣動(dòng)發(fā)射裝置結(jié)構(gòu)簡圖Fig.1 Sketch of an underwater pneumatic launcher

    1.1 發(fā)射氣瓶模型

    由理想氣體狀態(tài)方程和絕熱過程狀態(tài)方程可推出發(fā)射氣瓶的溫度和壓強(qiáng)的變化規(guī)律:

    式中:TB0、TB分別為初始時(shí)刻(定義為發(fā)射閥開啟時(shí)刻)、發(fā)射過程中發(fā)射氣瓶內(nèi)空氣的溫度,K;pB0、pB分別為初始時(shí)刻、發(fā)射過程中發(fā)射氣瓶內(nèi)空氣的壓強(qiáng),Pa;mB0、mB分別為初始時(shí)刻、發(fā)射過程中發(fā)射氣瓶內(nèi)空氣的質(zhì)量,mg;κ為絕熱指數(shù),κ=1.4。

    不計(jì)空氣泄露,經(jīng)發(fā)射閥流入發(fā)射管的空氣質(zhì)量為mi,質(zhì)量流量為

    1.2 發(fā)射閥模型

    式(3)表示的質(zhì)量流量又可寫成如下形式[10]:

    式中:φ為發(fā)射閥的流量系數(shù),經(jīng)驗(yàn)值取0.6~0.65;Ai為發(fā)射閥的瞬時(shí)開啟面積,m2;ρi為發(fā)射閥出口截面處的空氣密度,kg/m3;ui為發(fā)射閥出口截面處的氣體流速,m/s。由噴管理論可知,發(fā)射閥出口處流速和空氣密度分別為

    式中:R為理想氣體常數(shù),取287 J/(kg·K);pC為發(fā)射過程中發(fā)射管的內(nèi)壓,Pa;β=[2/(κ +1)]κ/(κ-1)為臨界壓力比;Ti為發(fā)射閥出口截面處的溫度:

    式中:pi為發(fā)射閥出口截面處的壓強(qiáng),Pa。

    1.3 發(fā)射閥開啟規(guī)律

    發(fā)射閥的開啟規(guī)律對于內(nèi)彈道控制起著重要的作用,其開啟面積通??梢员硎境呻S時(shí)間變化的連續(xù)函數(shù)或階段函數(shù)。由于下文中的模型試驗(yàn)中發(fā)射閥開啟到最大面積后一直保持直至武器出管后才關(guān)閉,為便于與試驗(yàn)結(jié)果比較,假定發(fā)射閥為理想的快開閥模型,其開啟面積隨時(shí)間的變化規(guī)律如圖2所示。其中t0為發(fā)射閥開啟動(dòng)作的響應(yīng)時(shí)間,AV0為最大開啟面積。

    圖2 發(fā)射閥開啟規(guī)律Fig.2 Open law of the launch valve

    1.4 發(fā)射管模型

    因武器的保護(hù)套管與發(fā)射管壁的間隙很小,又由于武器兩端的端蓋處有密封圈,因此密封片被沖破后外部水沒有立即進(jìn)入發(fā)射管且發(fā)射管內(nèi)氣體也無溢出。發(fā)射前后,發(fā)射管內(nèi)空氣體積為連續(xù)變化過程,忽略其與管壁間的熱交換,發(fā)射管內(nèi)空氣的熱力學(xué)變化過程可用熱力學(xué)第一定律描述:

    式中:UC、UC0分別為發(fā)射過程中和初始時(shí)刻發(fā)射管內(nèi)空氣的內(nèi)能,J;Hi為從發(fā)射閥流入發(fā)射管的空氣的熱焓,J;W為推動(dòng)武器運(yùn)動(dòng)的力所做的功,J。

    發(fā)射管內(nèi)空氣的溫度TC和壓強(qiáng)pC滿足:

    式中:mC為發(fā)射管內(nèi)空氣的質(zhì)量,kg;cV為空氣的定容比熱,J/(kg·K);VC為發(fā)射管內(nèi)的充氣容積,m3。

    通過發(fā)射閥注入發(fā)射管的空氣,其熱焓注入率為dHi/dt=cpTBdmi/dt。其中,cp為空氣的定壓比熱,J/(kg·K)??紤]到空氣經(jīng)發(fā)射閥節(jié)流后其滯止溫度有所降低,則節(jié)流后的空氣溫度滿足經(jīng)驗(yàn)關(guān)系:

    因此,經(jīng)由發(fā)射閥注入發(fā)射管的空氣熱焓注入率為

    發(fā)射管內(nèi)空氣的質(zhì)量變化率為

    式中:S為武器的橫截面面積,m2;vg為武器的運(yùn)動(dòng)速度,m/s;Fcr為推動(dòng)武器開始運(yùn)動(dòng)所需克服的初始阻力,主要來源于密封片的沖斷阻力。與魚雷發(fā)射裝置不同,密封片的存在除了保護(hù)管內(nèi)武器外,還使武器在運(yùn)動(dòng)前就已受到較大的推力作用。

    1.5 武器在發(fā)射管內(nèi)的運(yùn)動(dòng)

    武器的運(yùn)動(dòng)規(guī)律滿足:

    發(fā)射開始后,作用于武器底面的推力使武器克服初始阻力開始運(yùn)動(dòng)。此過程發(fā)射管內(nèi)充氣容積的變化率為

    式中:mg為武器質(zhì)量,kg;x為武器在發(fā)射管內(nèi)的位移;Rg為運(yùn)動(dòng)開始后武器所受的阻力,N,可表達(dá)為

    式中:ρw為水的密度,kg/m3;CD為武器水下運(yùn)動(dòng)的阻力系數(shù),s2/(m·kg);μ為武器與發(fā)射管管壁間的摩擦系數(shù);G為武器的重量,N;Fg為武器所受浮力,N;ph為發(fā)射管外的流體壓強(qiáng),Pa。

    運(yùn)動(dòng)過程中,武器所受推力的做功功率為

    利用MATLAB編寫程序求解方程(1)~(17)所示的發(fā)射系統(tǒng)數(shù)學(xué)模型,可對發(fā)射過程內(nèi)彈道結(jié)果進(jìn)行仿真分析。

    2 模型的有效性驗(yàn)證

    利用圖3所示的發(fā)射試驗(yàn)裝置進(jìn)行了模型發(fā)射試驗(yàn),圖4為某次發(fā)射過程中武器水下運(yùn)動(dòng)的截圖。為驗(yàn)證理論模型的有效性,試驗(yàn)中測量了武器在發(fā)射管內(nèi)的運(yùn)動(dòng)時(shí)間和發(fā)射后氣瓶氣體壓強(qiáng)。

    模型試驗(yàn)和仿真計(jì)算所用參數(shù)相同,主要包括:發(fā)射氣瓶容積0.04 m3,發(fā)射管內(nèi)初始空氣體積0.000 5 m3??諝獾亩▔罕葻?cp為 1 006 J/(kg·K),定容比熱cV為719 J/(kg·K)。武器質(zhì)量為10 kg,武器發(fā)射管內(nèi)行程為0.8 m,由試驗(yàn)確定的武器運(yùn)動(dòng)初始阻力Fcr為9 000 N。發(fā)射閥開啟動(dòng)作的響應(yīng)時(shí)間為50 ms,最大開啟截面為20 mm2。武器水下運(yùn)動(dòng)阻力系數(shù)CD取0.8 s2/(m·kg)。

    對發(fā)射氣瓶初始壓強(qiáng)為2.50 MPa的發(fā)射工況,利用數(shù)學(xué)模型仿真計(jì)算得到的武器管內(nèi)運(yùn)動(dòng)時(shí)間為0.138 s,發(fā)射后氣瓶壓強(qiáng)為 2.29 MPa。由水下模型發(fā)射試驗(yàn)測到的武器出管時(shí)間為0.133 s,發(fā)射后氣瓶壓強(qiáng)為2.25 MPa。數(shù)學(xué)模型結(jié)果與試驗(yàn)結(jié)果基本吻合,驗(yàn)證了仿真分析方法的有效性。

    圖3 發(fā)射試驗(yàn)裝置Fig.3 Experimental facilities of the launch

    圖4 發(fā)射過程截圖Fig.4 Screenshot of the launching process

    3 內(nèi)彈道參數(shù)分析

    3.1 局部靈敏度分析方法

    靈敏度分析用于研究一個(gè)系統(tǒng)的狀態(tài)或輸出變化對系統(tǒng)參數(shù)變化的敏感程度,其中局部靈敏度分析方法是應(yīng)用較廣泛的一種方法。使用該方法時(shí),一般假定系統(tǒng) y 是 n 個(gè)參數(shù) x1,x2,..,xn的函數(shù)。分析某個(gè)參數(shù)xi的靈敏度時(shí),通常只改變該參數(shù)的值,同時(shí)保持其余參數(shù)不變,由此計(jì)算函數(shù)值的改變量以衡量參數(shù)的影響程度。某個(gè)參數(shù)xi的局部靈敏度可用如下的偏導(dǎo)數(shù)表示[11]:

    式中:Si定義為函數(shù)在設(shè)計(jì)點(diǎn)處關(guān)于參數(shù)xi的靈敏度,的值越大,表示參數(shù)xi越靈敏。

    對于不同的參數(shù),由式(18)計(jì)算得到的結(jié)果單位不同,為便于比較,將其無量綱化:

    式中:yi是函數(shù)在設(shè)計(jì)點(diǎn)處的函數(shù)值。某些系統(tǒng)的解析形式難于得到,為進(jìn)行局部靈敏度分析,可將參數(shù)改變一個(gè)微小幅度(如0.5%)采用有限差分法計(jì)算靈敏度。

    3.2 局部靈敏度結(jié)果

    水下氣動(dòng)發(fā)射系統(tǒng)的內(nèi)彈道設(shè)計(jì)目標(biāo)是在發(fā)射膛壓允許的條件下使武器以設(shè)定的規(guī)律出管,因此可取發(fā)射過程中最高膛壓pC,max和武器的出管速度vg,exit作為內(nèi)彈道性能的目標(biāo)函數(shù),分析主要參數(shù)武器運(yùn)動(dòng)初始阻力Fcr、發(fā)射閥最大開啟面積AV0、發(fā)射氣瓶初壓pB0和體積VB等的影響規(guī)律和靈敏度。以水下試驗(yàn)的參數(shù)值為基準(zhǔn),將待分析參數(shù)分別減小和增大5%、10%和15%并保持其他參數(shù)不變,計(jì)算最高膛壓和武器出管速度的變化規(guī)律及參數(shù)靈敏度結(jié)果,分別如圖5、6所示。

    由圖5(a)可知,發(fā)射管最高膛壓與各參數(shù)間近似成線性關(guān)系。由圖5(b)可看出,武器運(yùn)動(dòng)初始阻力的靈敏度值最大。隨著武器運(yùn)動(dòng)初始阻力的增加,發(fā)射管最高膛壓值也增加。最高膛壓與初始阻力這種關(guān)系與發(fā)射過程內(nèi)彈道變化有關(guān):當(dāng)武器沖破密封片開始運(yùn)動(dòng)后,發(fā)射管內(nèi)空氣體積快速增加引起管內(nèi)膛壓由峰值處迅速下降,因此密封片被沖破瞬時(shí)恰好對應(yīng)著武器克服初始阻力和膛壓達(dá)到最高值的時(shí)刻。

    由圖6(a)可以看出,武器出管速度與各參數(shù)間也近似成線性關(guān)系,且隨參數(shù)的增加而增大。由圖6(b)可知,各參數(shù)中發(fā)射閥最大開啟面積和氣瓶初壓的靈敏度值最大且變化規(guī)律接近,武器運(yùn)動(dòng)初始阻力的靈敏度值較小,氣瓶體積的靈敏度值最小。發(fā)射閥最大開啟面積和發(fā)射氣瓶初壓,與發(fā)射過程中空氣推力對武器的做功大小直接相關(guān),因而對武器出管速度的影響最為明顯。

    圖5 最高膛壓的參數(shù)影響規(guī)律Fig.5 Effect law of parameters on the maximum chamber pressure of the launch tube

    圖6 出管速度的參數(shù)影響規(guī)律Fig.6 Effect law of parameters on the exit velocity of the weapon

    通過以上分析可知,為滿足發(fā)射管最高膛壓的允許條件,調(diào)整運(yùn)動(dòng)初始阻力(主要是調(diào)整密封金屬片的厚度)效果最顯著。為保證武器的出管速度,調(diào)整發(fā)射氣瓶初壓和發(fā)射閥最大開啟面積效果最明顯,調(diào)整運(yùn)動(dòng)初始阻力效果次之。

    4 結(jié)論

    對某型水下氣動(dòng)發(fā)射系統(tǒng)進(jìn)行了內(nèi)彈道數(shù)學(xué)建模和參數(shù)靈敏度分析。得到如下結(jié)論:

    1)通過逐個(gè)建立發(fā)射系統(tǒng)各環(huán)節(jié)數(shù)學(xué)模型得到了發(fā)射系統(tǒng)的整體模型,發(fā)射過程武器管內(nèi)運(yùn)動(dòng)時(shí)間和發(fā)射氣瓶剩余壓強(qiáng)的理論和試驗(yàn)結(jié)果吻合,驗(yàn)證了該型發(fā)射裝置內(nèi)彈道模型的正確性。

    2)通過局部靈敏度分析可知:初始阻力對最高膛壓的影響最顯著,阻力越大最高膛壓越大。發(fā)射閥最大開啟面積、氣瓶初壓和氣瓶體積對最高膛壓的影響較小。發(fā)射閥最大開啟面積和發(fā)射氣瓶初壓對出管速度的影響最明顯,出管速度隨這2個(gè)參數(shù)的增加而增加,運(yùn)動(dòng)初始阻力對出管速度的影響較小,氣瓶體積的影響最小。

    3)單參數(shù)變化的局部靈敏度分析方法,未能考慮多參數(shù)同時(shí)變化的影響,更深入的參數(shù)影響分析可采用全局靈敏度方法。

    4)內(nèi)彈道性能主要影響參數(shù)的確定和影響規(guī)律的研究,可為該型發(fā)射裝置內(nèi)彈道的優(yōu)化設(shè)計(jì)提供參考。

    [1]練永慶,王樹宗.魚雷發(fā)射裝置設(shè)計(jì)原理[M].北京:國防工業(yè)出版社,2012:169-280.

    [2]田兵,王樹宗,練永慶.液壓蓄能式魚雷發(fā)射裝置內(nèi)彈道建模與仿真[J].魚雷技術(shù),2011,19(1):68-71.

    TIAN Bing,WANG Shuzong,LIAN Yongqing.Inner trajectory modeling and simulation of hydraulic energy accumulated torpedo launcher[J].Torpedo Technology,2011,19(1):68-71.

    [3]練永慶,田兵,王樹宗,等.空氣渦輪泵發(fā)射系統(tǒng)發(fā)射過程仿真分析[J].兵工學(xué)報(bào),2011,32(2):155-162.

    LIAN Yongqing,TIAN Bing,WANG Shuzong,et al.Simulation and analysis of the launching process of air-turbine pump launch system[J].Acta Armamentarii,2011,32(2):155-162.

    [4]宋保維,黃海.氣動(dòng)不平衡式魚雷發(fā)射器內(nèi)彈道流場仿真研究[J].魚雷技術(shù),2008,16(2):46-49.

    SONG Baowei,HUANG Hai.Simulation of inner trajectory flow field for pneumatic imbalance torpedo launcher[J].Torpedo Technology,2008,16(2):46-49.

    [5]劉樂華,張宇文,袁緒龍.水下大深度垂直發(fā)射內(nèi)流場的數(shù)值研究[J].水動(dòng)力學(xué)研究與進(jìn)展,2005,20(1):90-94.

    LIU Lehua,ZHANG Yuwen,YUAN Xulong.Numerical simulation of the interior flow for launch under deep water[J].Journal of Hydrodynamics,2005,20(1):90-94.

    [6]練永慶,徐勤超,王樹宗.空氣渦輪機(jī)帶負(fù)載啟動(dòng)過程仿真[J].機(jī)械工程學(xué)報(bào),2011,47(8):191-198.

    LIAN Yongqing,XU Qinchao,WANG Shuzong.Simulation of startup process of air turbine with load[J].Journal of Mechanical Engineering,2011,47(8):191-198.

    [7]吳朝暉,王樹宗,馬世杰.旋轉(zhuǎn)泵式魚雷發(fā)射系統(tǒng)建模與仿真[J].海軍工程大學(xué)學(xué)報(bào),2000,12(5):54-58.

    WU Zhaohui,WANG Shuzong,MA Shijie.Building model and simulation of submarine rotary pump torpedo launch system[J].Journal of Naval University of Engineering,2000,12(5):54-58.

    [8]張永,逄洪照,黃波,等.氣動(dòng)不平衡發(fā)射裝置發(fā)射水雷內(nèi)彈道特性[J].艦船科學(xué)技術(shù),2011,33(12):93-96.

    ZHANG Yong,PANG Hongzhao,HUANG Bo,et al.Research on interior trajectory characteristics of the seamine launched by the aerodynamic non-balance torpedo launcher[J].Ship Science and Technology,2011,33(12):93-96.

    [9]陳宜輝,王樹宗,練永慶,等.水壓平衡式水下發(fā)射裝置的可調(diào)節(jié)發(fā)射閥研究[J].機(jī)械科學(xué)與技術(shù),2006,25(2):210-211.

    CHEN Yihui,WANG Shuzong,LIAN Yongqing,et al.Research on the adjustable discharge valve of underwater hydraulic balanced launcher[J].Mechanical Science and Technology,2006,25(2):210-211.

    [10]喬汝椿.渦輪泵式潛艇魚雷發(fā)射裝置設(shè)計(jì)分析[J].魚雷技術(shù),2004,12(2):43-46.

    QIAO Ruchun.A design analysis of pneumatic turbine pump-type torpedo launcher[J].Torpedo Technology,2004,12(2):43-46.

    [11]束龍倉,王茂枚,劉瑞國,等.地下水?dāng)?shù)值模擬中的參數(shù)靈敏度分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2007,35(5):491-495.

    SHU Longcang,WANG Maomei,LIU Ruiguo,et al.Sensitivity analysis of parameters in numerical simulation of groundwater[J].Journal of Hohai University:Natural Sciences,2007,35(5):491-495.

    猜你喜歡
    發(fā)射管發(fā)射裝置氣瓶
    UUV 變截面自航發(fā)射管線型優(yōu)化仿真分析
    某發(fā)射裝置隱形設(shè)計(jì)研究與分析
    航空火箭發(fā)射器平行度測試方法研究
    科技視界(2020年15期)2020-08-04 16:20:35
    新形勢下艦載垂直發(fā)射裝置發(fā)展趨勢
    火箭炮方形發(fā)射管流場仿真分析
    《壓縮氣體氣瓶充裝規(guī)定》解讀
    Worthington公司宣布推出全回收復(fù)合材料氣瓶
    復(fù)雜戰(zhàn)場環(huán)境導(dǎo)彈發(fā)射裝置隱身防護(hù)技術(shù)研究
    大同專項(xiàng)檢查氣瓶充裝站要求嚴(yán)格落實(shí)操作規(guī)章
    同心筒垂直發(fā)射裝置底部導(dǎo)流分析
    国产成人一区二区在线| 国产欧美日韩精品一区二区| av在线老鸭窝| 色综合色国产| 日韩av在线大香蕉| 亚洲va在线va天堂va国产| 91精品国产九色| 国产单亲对白刺激| 黄色一级大片看看| 亚洲精品日韩av片在线观看| av在线天堂中文字幕| 国产亚洲精品久久久com| 亚洲精品国产av蜜桃| 一级二级三级毛片免费看| 欧美xxxx性猛交bbbb| 亚洲欧美日韩东京热| 男女下面进入的视频免费午夜| 六月丁香七月| 免费人成在线观看视频色| 久久久成人免费电影| 亚洲成人久久爱视频| 在线免费观看的www视频| 国产色婷婷99| 国产视频首页在线观看| 99re6热这里在线精品视频| 一个人免费在线观看电影| 亚洲人成网站在线观看播放| 国产麻豆成人av免费视频| 亚洲高清免费不卡视频| 亚洲av二区三区四区| 免费观看无遮挡的男女| 大陆偷拍与自拍| 久久草成人影院| 床上黄色一级片| 日本色播在线视频| 久久精品久久久久久久性| 日韩精品有码人妻一区| 中文字幕亚洲精品专区| 国产成人免费观看mmmm| 中文字幕人妻熟人妻熟丝袜美| 久久久午夜欧美精品| 波多野结衣巨乳人妻| 一级毛片电影观看| 午夜免费男女啪啪视频观看| 国产黄色小视频在线观看| 午夜福利在线观看免费完整高清在| 听说在线观看完整版免费高清| 久久久久国产网址| 亚洲四区av| 午夜日本视频在线| 成年av动漫网址| 国产黄a三级三级三级人| 国产午夜精品论理片| 欧美成人a在线观看| 国产爱豆传媒在线观看| 午夜精品一区二区三区免费看| 男女视频在线观看网站免费| 中文天堂在线官网| 欧美高清性xxxxhd video| 久久久精品欧美日韩精品| 婷婷六月久久综合丁香| 亚洲最大成人中文| 久久99热6这里只有精品| 日日干狠狠操夜夜爽| 高清日韩中文字幕在线| 午夜福利高清视频| 成年版毛片免费区| 久久久久久久国产电影| 久久久久九九精品影院| 亚洲精品色激情综合| 日韩人妻高清精品专区| 国产一区二区三区综合在线观看 | 天堂中文最新版在线下载 | 色哟哟·www| 91午夜精品亚洲一区二区三区| 国产在线男女| 欧美成人午夜免费资源| av女优亚洲男人天堂| 嫩草影院新地址| 亚洲精品色激情综合| 免费观看a级毛片全部| 国产高潮美女av| 美女被艹到高潮喷水动态| 国产精品久久久久久久久免| 午夜福利视频1000在线观看| 国产精品女同一区二区软件| 男女下面进入的视频免费午夜| 国产女主播在线喷水免费视频网站 | 亚洲四区av| 精品久久国产蜜桃| 亚洲国产av新网站| 成人毛片a级毛片在线播放| 中文在线观看免费www的网站| 性色avwww在线观看| 综合色丁香网| 久久精品夜夜夜夜夜久久蜜豆| 最近中文字幕2019免费版| 国产精品不卡视频一区二区| 欧美人与善性xxx| 精品一区二区三区人妻视频| 在线a可以看的网站| a级一级毛片免费在线观看| 天堂影院成人在线观看| 成年人午夜在线观看视频 | 国产精品久久视频播放| 国产亚洲最大av| 热99在线观看视频| 国产爱豆传媒在线观看| 免费观看av网站的网址| 伦理电影大哥的女人| 国产女主播在线喷水免费视频网站 | 深爱激情五月婷婷| 国产淫片久久久久久久久| 18禁动态无遮挡网站| 久久国内精品自在自线图片| 一本久久精品| 色吧在线观看| 国产三级在线视频| 久久99精品国语久久久| 亚洲av日韩在线播放| 国内精品一区二区在线观看| 午夜精品一区二区三区免费看| 又爽又黄a免费视频| 高清毛片免费看| 18禁动态无遮挡网站| 国产91av在线免费观看| 亚洲图色成人| 国产综合懂色| 国产成人a区在线观看| 青春草视频在线免费观看| 九九爱精品视频在线观看| 超碰97精品在线观看| 国产中年淑女户外野战色| 久久精品夜色国产| 51国产日韩欧美| 中文乱码字字幕精品一区二区三区 | 一本久久精品| 久久韩国三级中文字幕| 天天躁日日操中文字幕| 国产黄a三级三级三级人| 最近最新中文字幕免费大全7| 久久久久久久久大av| 久久99蜜桃精品久久| 亚洲欧美日韩无卡精品| 婷婷色麻豆天堂久久| 色综合色国产| 看黄色毛片网站| 免费人成在线观看视频色| 欧美97在线视频| 老女人水多毛片| 成人二区视频| 人妻系列 视频| 午夜老司机福利剧场| 美女xxoo啪啪120秒动态图| 亚洲美女视频黄频| 免费高清在线观看视频在线观看| 国产高清不卡午夜福利| 中国美白少妇内射xxxbb| 亚洲久久久久久中文字幕| 麻豆精品久久久久久蜜桃| 国产美女午夜福利| av女优亚洲男人天堂| 亚洲av免费高清在线观看| 美女内射精品一级片tv| 美女内射精品一级片tv| 在线观看免费高清a一片| av在线天堂中文字幕| 亚洲欧洲日产国产| 国产成年人精品一区二区| 久久精品久久久久久噜噜老黄| 麻豆久久精品国产亚洲av| 夫妻性生交免费视频一级片| 国产成人a区在线观看| 亚洲国产av新网站| 2022亚洲国产成人精品| 国产成人一区二区在线| 国产精品一二三区在线看| 婷婷色综合www| 免费大片18禁| 久久精品综合一区二区三区| 色哟哟·www| 久久这里只有精品中国| 人人妻人人澡欧美一区二区| 国产爱豆传媒在线观看| 夜夜看夜夜爽夜夜摸| 久久鲁丝午夜福利片| 国产免费又黄又爽又色| 国产免费又黄又爽又色| av线在线观看网站| 中文字幕免费在线视频6| 男人狂女人下面高潮的视频| 乱系列少妇在线播放| 亚洲av免费高清在线观看| 亚洲国产色片| 在线免费十八禁| 精品人妻熟女av久视频| 狂野欧美白嫩少妇大欣赏| 男插女下体视频免费在线播放| 少妇猛男粗大的猛烈进出视频 | 亚洲精品久久久久久婷婷小说| freevideosex欧美| 婷婷六月久久综合丁香| 爱豆传媒免费全集在线观看| 国产成人福利小说| 听说在线观看完整版免费高清| 国产免费视频播放在线视频 | 成人亚洲精品一区在线观看 | av女优亚洲男人天堂| 国产麻豆成人av免费视频| 特大巨黑吊av在线直播| 又黄又爽又刺激的免费视频.| 一级毛片aaaaaa免费看小| 欧美人与善性xxx| 精品不卡国产一区二区三区| av播播在线观看一区| 日韩三级伦理在线观看| 久久综合国产亚洲精品| 深爱激情五月婷婷| 精品久久久精品久久久| 亚洲精品第二区| 高清在线视频一区二区三区| 免费av观看视频| 中文资源天堂在线| 美女脱内裤让男人舔精品视频| 国产亚洲5aaaaa淫片| 婷婷色综合大香蕉| 床上黄色一级片| 日本免费a在线| 最近手机中文字幕大全| 狂野欧美白嫩少妇大欣赏| 欧美日韩亚洲高清精品| 一级爰片在线观看| 在线播放无遮挡| 免费av不卡在线播放| 亚洲欧美一区二区三区黑人 | 美女内射精品一级片tv| 两个人视频免费观看高清| videossex国产| 一区二区三区乱码不卡18| 韩国av在线不卡| 少妇熟女欧美另类| 久久久久精品性色| 久久国内精品自在自线图片| 一级黄片播放器| 成人一区二区视频在线观看| 国产爱豆传媒在线观看| 99久久精品热视频| 国产精品人妻久久久影院| 国产视频首页在线观看| 在线观看免费高清a一片| 国模一区二区三区四区视频| 99热6这里只有精品| 亚洲精品国产av蜜桃| 国产综合懂色| 亚洲最大成人中文| 精品不卡国产一区二区三区| 日韩人妻高清精品专区| 久久精品夜夜夜夜夜久久蜜豆| 国产成人福利小说| 在线观看免费高清a一片| 国产亚洲av嫩草精品影院| 亚洲av电影不卡..在线观看| 国产激情偷乱视频一区二区| av在线亚洲专区| 欧美区成人在线视频| 免费高清在线观看视频在线观看| 久久99热这里只有精品18| 亚洲av电影在线观看一区二区三区 | 成年版毛片免费区| 亚洲第一区二区三区不卡| 汤姆久久久久久久影院中文字幕 | 亚州av有码| 精品一区在线观看国产| 日日摸夜夜添夜夜爱| 大香蕉97超碰在线| 99热全是精品| 精品亚洲乱码少妇综合久久| 91精品一卡2卡3卡4卡| 欧美人与善性xxx| 边亲边吃奶的免费视频| 国产不卡一卡二| 精品一区在线观看国产| 亚洲精品国产av成人精品| 日本黄色片子视频| 亚洲三级黄色毛片| 一个人看的www免费观看视频| 国产极品天堂在线| 七月丁香在线播放| 免费大片18禁| 亚洲精品国产av成人精品| 亚洲精品日本国产第一区| 国产精品蜜桃在线观看| 亚洲在线自拍视频| 国产精品99久久久久久久久| 精品酒店卫生间| 精品不卡国产一区二区三区| av免费观看日本| 国产亚洲最大av| 有码 亚洲区| 一级毛片久久久久久久久女| 深爱激情五月婷婷| 国产成人a∨麻豆精品| 亚洲va在线va天堂va国产| 麻豆av噜噜一区二区三区| 色尼玛亚洲综合影院| 久久久a久久爽久久v久久| 中文在线观看免费www的网站| 国产久久久一区二区三区| 国产色爽女视频免费观看| 久久久久九九精品影院| 春色校园在线视频观看| 男人爽女人下面视频在线观看| 国产亚洲一区二区精品| 亚洲精品中文字幕在线视频 | 春色校园在线视频观看| 亚洲成人久久爱视频| 午夜福利成人在线免费观看| 婷婷色av中文字幕| 久久久久久久久中文| 91久久精品电影网| 国产极品天堂在线| 亚洲综合色惰| 丰满人妻一区二区三区视频av| 97在线视频观看| 天堂√8在线中文| 国产精品女同一区二区软件| 青春草亚洲视频在线观看| 国产精品女同一区二区软件| 欧美激情在线99| 亚洲精品影视一区二区三区av| 欧美xxxx性猛交bbbb| 国产精品蜜桃在线观看| 精品久久久噜噜| 午夜视频国产福利| 成人毛片60女人毛片免费| av播播在线观看一区| 国产精品国产三级国产av玫瑰| 大香蕉97超碰在线| 日韩电影二区| 日韩,欧美,国产一区二区三区| 91久久精品电影网| 亚洲精品,欧美精品| 国产成人91sexporn| xxx大片免费视频| 中国国产av一级| 国产精品久久久久久av不卡| 成年人午夜在线观看视频 | 国产精品久久久久久av不卡| 久久精品夜夜夜夜夜久久蜜豆| 七月丁香在线播放| 偷拍熟女少妇极品色| 熟女人妻精品中文字幕| 中文在线观看免费www的网站| 国产又色又爽无遮挡免| 男女边摸边吃奶| 日本猛色少妇xxxxx猛交久久| 午夜精品在线福利| 九九爱精品视频在线观看| 熟妇人妻久久中文字幕3abv| 免费观看a级毛片全部| 免费大片18禁| 久久久久性生活片| kizo精华| 久久99蜜桃精品久久| 日韩强制内射视频| 老司机影院成人| 国产精品麻豆人妻色哟哟久久 | 美女xxoo啪啪120秒动态图| 国产黄色视频一区二区在线观看| 黄色欧美视频在线观看| 天堂中文最新版在线下载 | 国产成人freesex在线| 婷婷色综合大香蕉| 亚洲高清免费不卡视频| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 国产极品天堂在线| 深爱激情五月婷婷| 久久久久性生活片| 欧美激情国产日韩精品一区| 日韩av不卡免费在线播放| 夫妻午夜视频| 国产成人aa在线观看| 成人一区二区视频在线观看| 日本-黄色视频高清免费观看| 国产成人免费观看mmmm| 少妇熟女欧美另类| 日本免费在线观看一区| 亚洲av一区综合| 亚洲人成网站在线播| 看十八女毛片水多多多| 又大又黄又爽视频免费| 精品久久久精品久久久| 一级爰片在线观看| 我的老师免费观看完整版| 丰满人妻一区二区三区视频av| 午夜福利视频1000在线观看| 久久精品国产鲁丝片午夜精品| 大香蕉久久网| 亚洲精品视频女| 舔av片在线| 欧美成人精品欧美一级黄| 亚洲,欧美,日韩| 日韩不卡一区二区三区视频在线| 一个人观看的视频www高清免费观看| 久久精品国产鲁丝片午夜精品| 国产精品1区2区在线观看.| 欧美日韩精品成人综合77777| 深爱激情五月婷婷| 女的被弄到高潮叫床怎么办| 久久久久国产网址| 国产精品av视频在线免费观看| 国产高清国产精品国产三级 | 国产成人午夜福利电影在线观看| 永久免费av网站大全| 丰满乱子伦码专区| 亚洲成人一二三区av| 欧美不卡视频在线免费观看| 91av网一区二区| 特级一级黄色大片| 成人毛片60女人毛片免费| 美女国产视频在线观看| 亚洲成人久久爱视频| 亚洲,欧美,日韩| 日韩制服骚丝袜av| 亚洲国产精品成人综合色| 国产成人a区在线观看| 男人舔奶头视频| 97在线视频观看| 日韩av在线大香蕉| 亚洲人成网站高清观看| 91av网一区二区| 日本黄色片子视频| 日韩av在线免费看完整版不卡| 亚洲国产最新在线播放| 亚洲在久久综合| or卡值多少钱| 秋霞伦理黄片| 99re6热这里在线精品视频| 亚洲无线观看免费| 亚洲,欧美,日韩| 欧美97在线视频| 国内精品一区二区在线观看| 久久国产乱子免费精品| 久久久久精品性色| 精品亚洲乱码少妇综合久久| 好男人视频免费观看在线| 国语对白做爰xxxⅹ性视频网站| 青春草国产在线视频| 久久久久久久大尺度免费视频| 男人舔奶头视频| 能在线免费观看的黄片| 天天一区二区日本电影三级| 亚洲综合色惰| 国产精品国产三级国产专区5o| 大陆偷拍与自拍| 嫩草影院精品99| 国产一区有黄有色的免费视频 | 国产视频首页在线观看| 免费看日本二区| 99视频精品全部免费 在线| 我要看日韩黄色一级片| 久久精品国产鲁丝片午夜精品| 男女边摸边吃奶| 好男人视频免费观看在线| 免费大片黄手机在线观看| 色综合站精品国产| 如何舔出高潮| 日韩欧美 国产精品| 免费看美女性在线毛片视频| 男女那种视频在线观看| 只有这里有精品99| 性色avwww在线观看| 伦精品一区二区三区| 一级二级三级毛片免费看| av免费观看日本| 一本久久精品| 人人妻人人看人人澡| 一级毛片 在线播放| 在线免费十八禁| 不卡视频在线观看欧美| 一级片'在线观看视频| 国产91av在线免费观看| 人妻夜夜爽99麻豆av| 一区二区三区四区激情视频| 国产色婷婷99| 国内揄拍国产精品人妻在线| 午夜老司机福利剧场| 国产伦理片在线播放av一区| 色综合亚洲欧美另类图片| 国产乱来视频区| av在线蜜桃| 欧美日韩国产mv在线观看视频 | 久久久精品94久久精品| 伊人久久国产一区二区| 亚洲精品亚洲一区二区| 一二三四中文在线观看免费高清| 可以在线观看毛片的网站| 99久久精品热视频| 国产精品伦人一区二区| 边亲边吃奶的免费视频| 丰满少妇做爰视频| 日韩国内少妇激情av| 亚洲自偷自拍三级| 成年av动漫网址| 国产三级在线视频| 国产精品一区www在线观看| 精品不卡国产一区二区三区| 校园人妻丝袜中文字幕| 国产 一区 欧美 日韩| eeuss影院久久| 三级国产精品片| 国产一级毛片七仙女欲春2| 性插视频无遮挡在线免费观看| 欧美高清成人免费视频www| 亚洲综合色惰| 干丝袜人妻中文字幕| 在线观看av片永久免费下载| 久久久久九九精品影院| 免费看光身美女| 国产国拍精品亚洲av在线观看| 丝瓜视频免费看黄片| 大话2 男鬼变身卡| 能在线免费看毛片的网站| 天天一区二区日本电影三级| 国产成人福利小说| 18禁动态无遮挡网站| 精品一区二区三区人妻视频| 欧美成人一区二区免费高清观看| 91狼人影院| 日韩强制内射视频| 啦啦啦韩国在线观看视频| 国产成人91sexporn| 久久99热6这里只有精品| 国产有黄有色有爽视频| 女人十人毛片免费观看3o分钟| 亚洲无线观看免费| or卡值多少钱| 夜夜爽夜夜爽视频| 99热网站在线观看| 欧美极品一区二区三区四区| 久久99精品国语久久久| 2021少妇久久久久久久久久久| av天堂中文字幕网| 一二三四中文在线观看免费高清| 熟妇人妻久久中文字幕3abv| 青春草亚洲视频在线观看| 波多野结衣巨乳人妻| 日韩大片免费观看网站| www.色视频.com| 亚洲av不卡在线观看| 日本色播在线视频| 超碰97精品在线观看| freevideosex欧美| 51国产日韩欧美| 大话2 男鬼变身卡| 日韩欧美精品免费久久| 久久久久久久久久久免费av| 国精品久久久久久国模美| 日韩,欧美,国产一区二区三区| 97在线视频观看| a级一级毛片免费在线观看| 亚洲精品色激情综合| 欧美成人午夜免费资源| 亚洲熟女精品中文字幕| 啦啦啦啦在线视频资源| 亚洲精品亚洲一区二区| 禁无遮挡网站| 免费观看av网站的网址| 欧美xxxx黑人xx丫x性爽| h日本视频在线播放| 六月丁香七月| 国产黄片美女视频| 大片免费播放器 马上看| 少妇熟女aⅴ在线视频| 观看美女的网站| 国产成人一区二区在线| 国产综合精华液| 内射极品少妇av片p| 校园人妻丝袜中文字幕| 亚洲在久久综合| 日日干狠狠操夜夜爽| 少妇熟女aⅴ在线视频| 久久97久久精品| 日韩国内少妇激情av| 日日干狠狠操夜夜爽| 少妇熟女aⅴ在线视频| 久久97久久精品| 网址你懂的国产日韩在线| 美女内射精品一级片tv| av福利片在线观看| 18禁裸乳无遮挡免费网站照片| 精品国产一区二区三区久久久樱花 | 九色成人免费人妻av| 中文资源天堂在线| 欧美成人a在线观看| 人体艺术视频欧美日本| 春色校园在线视频观看| 永久网站在线| 婷婷色综合www| 夜夜看夜夜爽夜夜摸| 午夜精品在线福利| 91久久精品电影网| 在线免费观看不下载黄p国产| 亚洲精品日韩在线中文字幕| 免费在线观看成人毛片| 街头女战士在线观看网站| 91精品国产九色| 精华霜和精华液先用哪个| 亚洲人成网站在线播| www.色视频.com| 男女那种视频在线观看| 欧美精品国产亚洲| 精品久久久久久久人妻蜜臀av| 乱人视频在线观看|