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

    膛口初始流場對火藥燃氣射流的影響

    2017-05-07 06:11:31李子杰
    含能材料 2017年4期
    關鍵詞:炮口火藥彈丸

    李子杰, 王 浩

    (南京理工大學能源與動力工程學院, 江蘇 南京 210094)

    1 引 言

    對膛口流場的研究涉及氣體動力學、氣動聲學等多個學科,因此膛口流場的機理研究是一個具有挑戰(zhàn)性的課題。在膛內火藥燃氣推動作用下,加速運動的彈丸不斷壓縮彈前空氣,形成初始流場[1-2]。高溫高壓火藥燃氣在彈丸出炮口后迅速噴出,受到彈丸的阻礙并與初始流場融合,在膛口處形成具有復雜波系的膛口流場[3-4]。大口徑高超速火炮膛口流場的有害擾動極其強烈,但由于炮口焰、炮口煙、電磁干擾等的存在,使得瞬態(tài)變化的膛口流場的激波結構很難測量和計算[5]。

    人們對膛口流場形成的發(fā)展機理、膛口沖擊波場的分布規(guī)律、膛口氣流對彈丸的后效作用等方面都有了較為清晰的認識[6]。Love[7]等通過實驗和理論方法對軸對稱自由射流進行了研究,分析了射流馬赫數(shù)、噴管形狀等因素對軸線上的膛口氣流參數(shù)分布的影響。李鴻志[8]等分析了膛口沖擊波與射流的特征及形成機理,提出了膛口沖擊波的變能量、動球心的球形物理模型,并總結出膛口沖擊波遠場傳播規(guī)律; 姜孝海[9]等用基于ALE方程及動網(wǎng)格技術對44 mm口徑彈丸由膛內到膛外直到飛離初始流場的整個過程進行數(shù)值模擬。江坤[10]利用結構網(wǎng)格對122 mm車載火炮的膛口流場進行數(shù)值模擬,對初速為713 m·s-1的彈丸發(fā)射過程中的流動現(xiàn)象進行了分析,得到了較為可信的結果。代淑蘭[11]采用三維非定?;瘜W反應流控制方程組對帶制退器的7.62 mm口徑、735 m·s-1初速的彈丸膛口燃燒流場進行了數(shù)值模擬,清晰的描述了膛口流場的發(fā)展過程結構和性質及彈丸與流場的相互影響。郭澤慶[12]等以7.62 mm彈道槍為研究對象,基于ALE方程及分區(qū)結構貼體網(wǎng)格對膛口初始流場及后效期過程進行了試驗分析。Jiang[13]基于二維結構網(wǎng)格分別用TVD格式和耗散控制格式離散Euler方程,對軸對稱膛口流場進行了數(shù)值模擬,考慮了圓柱形彈丸對流場的影響,并分析了彈丸與膛口氣流間的相互作用。張煥好[14]基于三維Euler方程,結合Roe格式與結構化動網(wǎng)格技術對含有制退器的20 mm槍炮膛口流場進行了數(shù)值模擬。

    前人所研究的重點多為燃氣流場的特性分析,且多以155 mm口徑以下,1000 m·s-1以下初速的槍炮為研究對象。而對于300 mm大口徑1730 m·s-1超高速彈丸的初始流場對火藥燃氣射流結構的影響分析未見研究。本研究采用有限體積法并結合Realizablek-ε湍流模型,使用結構動網(wǎng)格技術,建立兩種二維軸對稱數(shù)值仿真模型,研究超高速彈丸含初始流場和簡化的不含初始流場情況下的火藥燃氣流場結構特性,得到了該發(fā)射條件下初始流場對火藥燃氣流場分布的影響,為大口徑高超速火炮的設計與試驗提供理論指導,為深入研究奠定良好基礎。

    2 數(shù)值方法

    2.1 控制方程

    當不考慮外加熱和徹體力的影響時,二維軸對稱可壓縮非定常N-S方程組為[15-16]:

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    其中應力項為:

    (10)

    系數(shù)σ決定流動的類型,當σ=1時,為二維軸對稱流動模型,當σ=0時,為二維平面流動模型[11]。

    壓力由理想氣體方程給出,即:

    (11)

    2.2 湍流模型

    本研究所采用的湍流模型為Realizablek-ε模型[17]。該模型比起標準k-ε模型有兩個主要的不同點: ①Realizablek-ε模型為湍流粘性增加了一個公式; ②為耗散率增加了新的傳輸方程。

    引入Boussinesq的線性渦粘假設,雷諾應力表達式為:

    (12)

    其中:

    (13)

    (14)

    其中Reynolds剪切應力的Schwarz不等式:

    (15)

    3 邊界條件與計算模型

    3.1 網(wǎng)格劃分

    本研究涉及到兩種計算模型,第一種計算域不考慮初始流場的影響,為彈丸從炮口處開始運動,如圖1a所示。將整個計算域分為三個區(qū)域: 彈后Ⅰ區(qū)、彈前Ⅱ區(qū)、及膛口周圍流場Ⅲ區(qū)。膛口周圍流場區(qū)域長13 m,寬3 m。最小網(wǎng)格尺寸為2 mm,位于膛口附近。計算域邊界處網(wǎng)格相對稀疏。均采用結構化網(wǎng)格進行劃分,共10萬個網(wǎng)格,如圖1b所示。

    第二種計算域為彈丸自膛底運動出炮口整個過程,炮管總長29.6 m,管壁厚度為50 mm,如圖1c所示。膛口周圍流場區(qū)域與第一種模型相同,長13 m,寬3 m。最小網(wǎng)格尺寸為2 mm,位于膛口附近。計算域邊界處網(wǎng)格相對稀疏。共28萬個網(wǎng)格,如圖1d所示。

    3.2 邊界條件

    不考慮初始流場的計算模型中,設定膛口處為壓力入口邊界條件,使用自編的UDF程序對入口壓力進行設定,膛口周圍流場邊界為壓力出口邊界條件。設定身管、炮口為固壁邊界條件,彈丸運動由六自由度程序控制,設定對稱軸為軸邊界。壁面均假定為絕熱,壁面上溫度為臨近壁面網(wǎng)格點溫度。身管外流場區(qū)域初始壓力為1.01×105MPa,初始溫度為300 K。

    a. sketch model without precursor flow field

    b. mesh model without precursor flow field

    c. sketch model with precursor flow field

    d. mesh model with precursor flow field

    圖1計算模型示意圖

    Fig.1Schematicc diagrams of calculation models

    考慮初始流場的計算模型中,彈丸在膛內運動階段由內彈道程序控制,出膛口后由六自由度程序控制,設定膛底為壓力入口邊界條件,設定身管、炮口為固壁邊界條件,彈丸為移動固壁邊界條件。

    3.3 內彈道與后效期

    以300 mm大口徑超高速平衡炮為研究對象,彈丸質量為160 kg,總裝藥量為240 kg,平衡體重2000 kg,圖2a、圖2b分別給出了內彈道過程的彈底壓力和速度隨時間的變化規(guī)律。彈丸炮口初度為1730.4 m·s-1,炮口壓力為46.34 MPa,運動時間為37.6 ms。

    針對本算例,后效期炮口壓力隨時間變化公式如下[18-20]:p=46.3×106e-24.886t

    彈丸后效期的炮口壓力隨時間呈指數(shù)分布規(guī)律,其壓力分布如圖2c所示。

    a. pressure-time curve of projectile base for interior ballistic

    b. velocity-time curve for interior ballistic

    c. pressure-time curve in after-effect period

    圖2內彈道計算參數(shù)分布曲線

    Fig.2Interior ballistic calculation parameter distribution curves

    3.4 網(wǎng)格無關性驗證

    對于網(wǎng)格收斂性測試,本研究采用網(wǎng)格數(shù)為六十萬的無初始流場模型作為對比,以彈底壓力為參考值,如圖3所示,當兩種不同網(wǎng)格數(shù)量下彈底壓力相差很小時,可以認為計算結果與網(wǎng)格無關,由圖3可以看出,兩者誤差小于1%,由此可以認為本計算結果與網(wǎng)格無關。

    圖3作為網(wǎng)格無關性驗證的彈底壓力-時間曲線

    Fig.3Pressure-time curves of projectile base for grid convergence tests

    4 結果分析

    4.1 初始流場對燃氣射流的影響分析

    模擬結果基于對有無初始流場的對比上,圖4a~圖4d顯示了高超音速彈丸初始流場的形成過程。彈丸在膛內推動彈前空氣,形成了一系列壓縮波,壓縮波不斷向炮口方向傳播。隨著彈丸的加速運動,彈前空氣被擠出炮口形成初始流場。當彈前空氣被壓縮流出膛口后,氣體迅速膨脹,壓力降低,速度增大,并向周圍傳播。彈丸的圓柱部運動出炮口時,彈前氣體全部被壓出,被壓縮的空氣向前傳播使得后方壓力降低,周圍空氣流向該區(qū)域。

    a. t=34.3 ms

    c. t=35.4 ms

    b. t=34.6 ms

    d. t=37.6 ms

    圖4超高速彈丸初始流場速度等值線圖

    Fig.4Velocity contour of the hyper-velocity projectile precursor flow field

    隨著彈丸運動出炮口,彈后高溫高壓火藥氣體迅速噴出,形成了復雜而強烈的火藥氣體射流結構,圖5a~圖5d和圖6a~圖6d分別給出了無初始流場和有初始流場情況下的火藥燃氣射流波系結構。

    高溫高壓火藥氣體首先從彈尾周圍逸出,向側前方迅速發(fā)展傳播,由圖5a~圖5d所示,當無初始流場擾動時,火藥燃氣流場波陣面近似為球形,火藥氣體不能追趕并超越彈丸,火藥燃氣剛噴出時速度約為2169 m·s-1。彈丸側面形成了明顯的起始于彈丸頭部脫離彈丸側面向外延伸的剪切層。因為炮口噴射出的氣體速度大于彈丸的運動速度,彈丸底部形成了清晰的彈底激波。彈底激波在膛口沖擊波的作用下進一步加強并阻礙馬赫盤的生成,隨著彈底激波的作用越來越弱,馬赫盤逐漸增大,瓶狀激波不斷擴大,膛口沖擊波不斷向遠處傳播并衰減。炮口形成了由冠狀沖擊波、彈底激波、馬赫盤等構成的完整波系。

    a. t=0.134 ms

    b. t=0.754 ms

    c. t=1.410 ms

    d. t=3.973 ms

    圖5無初始流場條件下燃氣射流速度等值線圖

    Fig.5Velocity contour of the combustion gas jet flow without precursor flow field

    a. t=0.134 ms

    b. t=0.754 ms

    c. t=1.410 ms

    d. t=3.973 ms

    圖6有初始流場條件下燃氣射流速度等值線圖

    Fig.6Velocity contour of the combustion gas jet flow with precursor flow field

    由圖6a~圖6d可以看出有初始流場的火藥燃氣流場結構與無初始流場影響下的燃氣流場結構明顯不同?;鹚幦細馔苿优诳谥車驯怀跏剂鲌鰯_動過的空氣,形成了炮口沖擊波。在炮口初始流場的影響下,火藥燃氣追趕并超越彈丸,在彈體周圍形成了復雜的火藥燃氣團,火藥燃氣剛噴出膛口時速度約為2359 m·s-1,比無初始流場影響下的速度提高約200 m·s-1。隨著彈丸的運動,彈丸逐漸擺脫了火藥氣體對其的包圍。隨著彈丸的運動,彈底激波的作用越來越弱,馬赫盤逐漸增大,瓶狀激波不斷擴大,膛口沖擊波不斷向遠處傳播并衰減,同樣形成了由冠狀沖擊波、彈底激波、馬赫盤等構成的完整波系。

    在初始流場的干擾下,火藥燃氣流場的各參數(shù)也明顯不同。圖7給出了計算模型中監(jiān)測點的分布狀況, 其坐標分布見表1。為了研究無初始流場干擾下和有初始流場干擾下的火藥燃氣流場參數(shù)變化,在彈前膛口周圍處監(jiān)測了四個點的壓力和溫度變化,表2給出了監(jiān)測點在有初始流場和無初始流場情況下的最大壓力和最大溫度。

    圖7監(jiān)測點分布情況

    Fig.7The distribution of monitoring points

    表1監(jiān)測點的坐標分布

    Table1The distribution of monitoring points

    monitorpoint1234coordinates(x/mm,y/mm)(1000,150)(2000,150)(2000,500)(3000,200)

    圖8a~圖8d和圖9a~圖9d給出了有無初始流場情況下的監(jiān)測點壓力和溫度的對比曲線。由圖8可以看出,四個監(jiān)測點的最大壓力隨膛口距離的增加而減小。對于監(jiān)測點1,有初始流場擾動情況下的初始壓力比無初始流場擾動情況下的初始壓力低,而監(jiān)測點2、3正好相反。這是因為監(jiān)測點1位于膨脹波影響的區(qū)域內,膨脹波的影響使得此處的壓力低于環(huán)境壓力。而監(jiān)測點2、3位于壓縮波影響的區(qū)域內,壓縮波使得此處壓力高于環(huán)境壓力。而此時監(jiān)測點4位于初始流場沒有影響到的區(qū)域內,所以有無初始流場擾動情況下的初始壓力相同。值得注意的是,有初始流場擾動情況下的最大壓力低于無初始流場擾動情況下的最大壓力,這是因為初始流場推動了膛口前方空氣,形成了負壓,因此當火藥燃氣射流傳播至此時,有初始流場擾動的最大壓力低于無初始流場擾動下的最大壓力。

    a. monitor point 1

    b. monitor point 2

    c. monitor point 3

    d. monitor point 4

    圖8有無初始流場的監(jiān)測點壓力對比

    Fig.8Comparison of the pressure of monitor points with and without precursor flow field

    表2有無初始流場的最大壓力和溫度對比

    Table2The maximum pressure and temperature comparison with and without precursor flow field

    monitorpointcoordinates(x/mm,y/mm)maximumpressure/MPawithoutprecursorflowfieldwithprecursorflowfieldmaximumtemperature/Kwithoutprecursorflowfieldwithprecursorflowfield1(1000,150)1.7000.7392143.923475.662(2000,150)1.0320.5582120.752950.513(2000,500)1.3670.5651141.202152.864(3000,200)0.8650.6161246.511747.75

    a. monitor point 1

    c. monitor point 3

    b. monitor point 2

    d. monitor point 4

    圖9有無初始流場的監(jiān)測點溫度對比

    Fig.9Comparison of the temperature of monitor points with and without precursor flow field

    圖10為初始流場的溫度等值線圖。由圖10可以看出,經初始流場擾動過的彈前空氣的溫度分布有了明顯的變化,無初始流場下的膛外溫度為室溫300 K,而有初始流場的膛外溫度一般在1000~1500 K左右。因此,經初始流場擾動過的膛口前方溫度高于沒有初始流場擾動過的溫度。

    4.2 初始流場對彈丸運動的影響

    由于初始流場的存在,使得彈丸出炮口前的流場區(qū)域已經受到復雜的干擾,火藥燃氣對彈丸底部的作用也受到了影響。圖11a和圖11b分別給出了無初始流場和有初始流場情況下的彈底壓力和馬赫數(shù)分布情況。

    由圖11可以看出,無初始流場時后效期彈底最大壓力為40.3 MPa,有初始流場情況下為38.9 MPa。無初始流場情況下,彈底馬赫數(shù)最高達到4.13,在初始流場的影響下,彈底馬赫數(shù)出現(xiàn)波動,最大馬赫數(shù)為4.09。

    圖10初始流場溫度等值線圖

    Fig.10Temperature contour of the precursor flow field

    初始流場波系強度較弱,與火藥燃氣流場相差一個數(shù)量級,但初始流場的存在影響了火藥燃氣對彈丸的作用,使得彈底壓力減小,馬赫數(shù)減小。

    a. projectile-base pressure

    b. Mach number

    圖11有無初始流場彈底壓力和馬赫數(shù)分布

    Fig.11The projectile-base pressure andMach number disstribution with and without precursor flow field

    5 結 論

    針對300 mm大口徑1730 m·s-1超高速平衡炮,建立兩種數(shù)學模型,探討了無初始流場和有初始流場下火藥燃氣流場的異同,以及初始流場的存在對彈丸運動的影響。本研究編寫內彈道程序控制彈丸在膛內的運動,編寫后效期程序對炮口壓力進行設定、編寫六自由度程序控制彈丸在膛外的運動。

    (1) 當沒有初始流場影響時,火藥燃氣射流呈球狀,燃氣不能追趕上彈丸。而在初始流場的影響下,火藥燃氣可以追趕上彈丸,但由于初速高,并不能超越并包圍彈丸。

    (2) 初始流場的存在使得火藥燃氣剛噴出膛口時的速度增高。

    (3) 初始流場的干擾使得燃氣射流的最大壓力降低近一倍,溫度增加1000 K以上,彈底壓力和彈丸馬赫數(shù)稍有減小。

    參考文獻:

    [1] 余文杰, 余永剛. 環(huán)狀噴口對底排尾部流場影響的數(shù)值模擬[J]. 含能材料, 2014, 22(3): 318-324.

    YU Wen-jie, YU Yong-gang. Numerical simulation on effects of annular jet in base flow of base bleed[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2014, 22(3): 318-324.

    [2] WANG Jiang-feng, WU Yi-zhao. Parallel computation of 3-d hypersonic flows on unstructured hybrid meshes[J].TransactionsofNanjingUniversityofAeronautics&Astronautics, 2015, 22(3): 300-305.

    [3] 張輝, 譚俊杰, 催東明. 高速彈丸穿越不同帶膛口裝置的流場數(shù)值模擬[J]. 火炮發(fā)射與控制學報, 2007, (2): 48-51.

    ZHANG Hui, TAN Jun-jie, CUI Dong-ming. Numerical simulation method of flow field with muzzle attachments [J].JiournalofGunLaunch&Control, 2007, (2): 48-51.

    [4] Rehman H, Chung H, Joung T, et al. CFD analysis of sound pressure in tank gun muzzle silencer[J].JCentSouthUnivTechnol, 2011, (18): 2015-2020.

    [5] 馬大為, 樂貴高, 余世海. 膨脹不足超音速自由射流的數(shù)值模擬[J]. 計算物理學報, 1994, 11(3): 362-366.

    MA Da-wei, LE Gui-gao, YU Shi-hai. Numerical simulation of under expanded supersonic free-jet flow[J].ChineseJournalofComputationalPhysics, 1994, 11(3): 362-366.

    [6] 王仕松, 鄭堅, 賈長治, 等. 帶制退器的膛口流場數(shù)值模擬[J].火力與指揮控制學報, 2011, 36(2): 148-151.

    WANG Shi-song,ZHENG Jian,JIA Chang-zhi, et al. Numerical simulation of muzzle blast flow field with muzzle brake[J].FireControl&CommandControl, 2011, 36(2): 148-151.

    [7] Love E S. Initial inclination of the mixing boundary separating an exhausting supersonic jet from a supersonic ambient stream[J].UgeskriftforLaeger, 1999, 161(36): 5038-5039.

    [8] 李鴻志,劉曉利. 膛口變能量沖擊波特性分析—彈道初始參量的影響[J]. 兵工學報, 1993 (3): 17-21.

    LI Hong-zhi, LIU Xiao-li. An analysis on the characteristics of Variable energy blast wave at the muzzle-the influence of ballistics initial parameters[J].ActaArmamentarii, 1993(3): 17-21.

    [9] JIANG Xian-hai, CHEN Zhi-hua, FAN Bao-chun, et al. Numerical simulation of blast flow fields induced by a high-speed projectile[J].ShockWaves, 2008 (18): 205-212.

    [10] 江坤, 王浩. 基于動網(wǎng)格技術的膛口流場數(shù)值模擬[J].火炮發(fā)射與控制學報, 2010(3): 1-4.

    JIANG Kun, WANG Hao. Numerical simulation of muzzle flow field based on dynamic meshing technique[J].JournalofGunLaunch&Control, 2010(3): 1-4.

    [11] 代淑蘭, 許厚謙, 肖忠良. 帶制退器的膛口燃燒流場并行數(shù)值模擬[J].彈道學報, 2009, 21(4): 84-87.

    DAI Shu-lan, XU Hou-qian, SIAO Zhong-liang. Numerical simulation of gun muzzle flow field including movable boundary[J].JournalofBallistics, 2009, 21(4): 84-87.

    [12] 郭則慶, 王楊, 姜孝海,等. 膛口初始流場對火藥燃氣流場影響的數(shù)值研究[J]. 兵工學報, 2012, 33(6): 663-668.

    GUO Ze-qing, WANG Yang, JIANG Xiao-hai, et al. Numerical study on effects of precursor flow on muzzle propellant flow field[J].ActaArmamentarii, 2012, 33(6): 663-668.

    [13] 蘇曉鵬, 錢林方, 戴勁松.帶炮口裝置時某火炮膛口流場數(shù)值仿真[J].計算機仿真學報, 2009, 26(9): 15-18.

    SU Xiao-peng, QIAN Lin-fang, DAI Jin-song. Muzzle flow field simulation of gun with a muzzle attachment [J].ComputerSimulation, 2009, 26(9): 15-18.

    [14] 張煥好, 陳志華, 姜孝海, 等. 高速彈丸穿越不同制退器時的膛口流場波系結構研究[J].兵工學報, 2012, 33(5): 623-629.

    ZHANG Huan-hao, CHEN Zhi-hua, JIANG Xiao-hai, et al. Investigation on the blast wave structures of a high-speed projectile flying through different muzzle brakes[J].ActaArmamentarii, 2012, 33(5): 563-567.

    [15] 劉昕, 姜宗林, 王春, 等. 欠膨脹超聲速流不穩(wěn)定性機理的數(shù)值研究[J].力學學報, 2008, 40(5): 577-584.

    LIU Xin, JIANG Zong-lin, WANG Chun, et al. Simulation study of supersonic flow instability mechanism[J].ChineseJournalofTheoreticalandAppliedMechanics, 2008, 40(5): 577-584.

    [16] ZHANG Huan-hao, CHEN Zhi-hua, JIANG Xiao-hai, et al. Investigations on the exterior flow field and the efficiency of the muzzle brake[J].JournalofMechanicalScienceandTechnology, 2013, 27(2): 95-101.

    [17] 郭則慶, 姜孝海, 王楊. 后效期火藥燃氣加速彈丸的數(shù)值研究[J]. 高壓物理學報, 2012, 26(5): 564-570.

    GUO Ze-qing, JIANG Xiao-hai, WANG Yang. Numerical investigations on acceleration of projectile in post-effect period[J].ChineseJournalofHighPressurePhysics, 2012, 26(5): 564-570.

    [18] 高樹滋, 趙潤祥, 馬大為. 火炮計算流體力學[M]. 北京: 兵器工業(yè)出版社, 1995: 165-172.

    [19] 陶如意, 孫繼兵, 黃明, 等. 高低壓室平衡炮內彈道數(shù)值模擬及試驗研究[J]. 南京理工大學學報, 2006, 30(4): 478-485.

    TAO Ru-yi, SUN Ji-bing, HUANG Ming, et al. Interior ballistics numerical Ssimulation and experiment of balance gun with low and high pressure chamber[J].JournalofNanjingUniversityofScienceandTechnology, 2006, 30(4): 478-485.

    [20] 張博孜, 王浩, 王珊珊, 等. 大口徑平衡炮內彈道一維兩相流建模與仿真[J]. 彈道學報, 2014, 26(1): 17-21.

    猜你喜歡
    炮口火藥彈丸
    長生不老藥和火藥
    炮口制退器研究現(xiàn)狀及發(fā)展趨勢
    超高速撞擊下球形彈丸破碎特性仿真研究
    神奇的火藥
    大灰狼畫報(2022年2期)2022-05-23 12:57:10
    神秘的『彈丸』
    火藥的來歷
    “火藥弟弟”
    炮口結構參數(shù)對炮口性能的影響分析
    彈丸對預開孔混凝土靶體侵徹的實驗研究
    圓柱形彈丸高速撞擊薄板的碎片云特性數(shù)值模擬
    黄片大片在线免费观看| 一本大道久久a久久精品| 在线观看www视频免费| 青春草视频在线免费观看| 成人亚洲精品一区在线观看| 99热国产这里只有精品6| 欧美久久黑人一区二区| 性色av一级| 国产免费av片在线观看野外av| 99精国产麻豆久久婷婷| 日韩一卡2卡3卡4卡2021年| 国产99久久九九免费精品| 亚洲全国av大片| 国产一卡二卡三卡精品| 欧美日本中文国产一区发布| 亚洲人成电影免费在线| 99国产精品99久久久久| 国产精品影院久久| 99精国产麻豆久久婷婷| av又黄又爽大尺度在线免费看| 精品国产国语对白av| 三上悠亚av全集在线观看| 精品乱码久久久久久99久播| 免费看十八禁软件| 亚洲国产精品一区三区| 国产精品久久久久久精品古装| 蜜桃在线观看..| 亚洲国产欧美一区二区综合| 正在播放国产对白刺激| 欧美成人午夜精品| 人人妻人人爽人人添夜夜欢视频| 亚洲天堂av无毛| 精品第一国产精品| 亚洲精品中文字幕在线视频| 亚洲男人天堂网一区| 大型av网站在线播放| 女人精品久久久久毛片| 久久久久网色| 国产精品久久久人人做人人爽| 永久免费av网站大全| 色婷婷av一区二区三区视频| 咕卡用的链子| 精品亚洲乱码少妇综合久久| 男女之事视频高清在线观看| 亚洲免费av在线视频| 精品卡一卡二卡四卡免费| 亚洲专区中文字幕在线| 美女脱内裤让男人舔精品视频| 99精品久久久久人妻精品| 高潮久久久久久久久久久不卡| 午夜两性在线视频| 人人妻人人爽人人添夜夜欢视频| 国产国语露脸激情在线看| 国精品久久久久久国模美| 美女视频免费永久观看网站| 久久精品人人爽人人爽视色| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美一区二区三区久久| 免费少妇av软件| 欧美xxⅹ黑人| 精品一区在线观看国产| 亚洲精品国产区一区二| 亚洲av日韩精品久久久久久密| 在线观看免费午夜福利视频| 丰满少妇做爰视频| 欧美97在线视频| 亚洲国产精品一区二区三区在线| 丰满人妻熟妇乱又伦精品不卡| 水蜜桃什么品种好| 欧美黑人欧美精品刺激| 久9热在线精品视频| 亚洲色图综合在线观看| 人人妻,人人澡人人爽秒播| 国产成人精品在线电影| 大码成人一级视频| 亚洲精品第二区| 午夜福利乱码中文字幕| 亚洲精华国产精华精| 高清在线国产一区| 老熟妇仑乱视频hdxx| 岛国毛片在线播放| 久久久精品国产亚洲av高清涩受| 老司机深夜福利视频在线观看 | 日韩大码丰满熟妇| 精品国内亚洲2022精品成人 | 日日爽夜夜爽网站| 涩涩av久久男人的天堂| 国产有黄有色有爽视频| 欧美国产精品一级二级三级| 亚洲精品成人av观看孕妇| 国产精品国产av在线观看| 电影成人av| 国产成人免费观看mmmm| 日日爽夜夜爽网站| 在线观看一区二区三区激情| 大香蕉久久网| 69精品国产乱码久久久| 一本久久精品| 啦啦啦免费观看视频1| 高清黄色对白视频在线免费看| 国产成人啪精品午夜网站| 久热爱精品视频在线9| 91麻豆精品激情在线观看国产 | 亚洲少妇的诱惑av| 高清黄色对白视频在线免费看| 国产精品 欧美亚洲| 又紧又爽又黄一区二区| 岛国毛片在线播放| 国产亚洲欧美精品永久| 国产精品久久久久久精品古装| 亚洲中文av在线| 老司机午夜十八禁免费视频| 一级黄色大片毛片| 日韩熟女老妇一区二区性免费视频| 女警被强在线播放| 国产在视频线精品| 热99久久久久精品小说推荐| 亚洲天堂av无毛| 亚洲av片天天在线观看| 在线永久观看黄色视频| 99久久国产精品久久久| 日韩 亚洲 欧美在线| 99精品欧美一区二区三区四区| 欧美日韩视频精品一区| 国产精品久久久久成人av| 精品亚洲成国产av| 亚洲av国产av综合av卡| 一级毛片精品| 久久亚洲精品不卡| 美女主播在线视频| 五月天丁香电影| 另类亚洲欧美激情| 亚洲欧美成人综合另类久久久| 无限看片的www在线观看| 日韩精品免费视频一区二区三区| 欧美成狂野欧美在线观看| 久久 成人 亚洲| 自线自在国产av| 免费观看人在逋| 免费av中文字幕在线| 丰满饥渴人妻一区二区三| www.av在线官网国产| 最新在线观看一区二区三区| 少妇 在线观看| 久久中文看片网| av福利片在线| 一本大道久久a久久精品| 精品人妻在线不人妻| 国产成人系列免费观看| 视频在线观看一区二区三区| 欧美乱码精品一区二区三区| 久久狼人影院| 国产欧美日韩综合在线一区二区| 欧美日韩精品网址| 国产成人精品在线电影| 午夜久久久在线观看| 韩国高清视频一区二区三区| videosex国产| 男女高潮啪啪啪动态图| 精品亚洲成a人片在线观看| 满18在线观看网站| 窝窝影院91人妻| 精品亚洲成国产av| 亚洲熟女精品中文字幕| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲| 最黄视频免费看| 一区二区三区乱码不卡18| 精品一品国产午夜福利视频| 精品久久久久久电影网| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 在线av久久热| 超碰97精品在线观看| 午夜福利在线免费观看网站| 大片电影免费在线观看免费| 亚洲专区中文字幕在线| 十八禁人妻一区二区| 国产成人免费无遮挡视频| 国产一区有黄有色的免费视频| 久久精品成人免费网站| 99国产精品99久久久久| 亚洲成人国产一区在线观看| 欧美一级毛片孕妇| 狠狠精品人妻久久久久久综合| 精品免费久久久久久久清纯 | 91麻豆av在线| 不卡一级毛片| 99国产精品一区二区蜜桃av | 搡老乐熟女国产| 黑人巨大精品欧美一区二区蜜桃| 大香蕉久久成人网| 国产成+人综合+亚洲专区| 国产一卡二卡三卡精品| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美日韩另类电影网站| 国内毛片毛片毛片毛片毛片| 日日爽夜夜爽网站| 久久精品aⅴ一区二区三区四区| 欧美日韩视频精品一区| 亚洲欧美精品综合一区二区三区| 亚洲性夜色夜夜综合| 正在播放国产对白刺激| 巨乳人妻的诱惑在线观看| 国产色视频综合| 欧美乱码精品一区二区三区| 在线天堂中文资源库| 大码成人一级视频| 日韩欧美国产一区二区入口| 国产黄频视频在线观看| 中文字幕精品免费在线观看视频| 五月开心婷婷网| 欧美另类亚洲清纯唯美| videos熟女内射| 青草久久国产| 国产精品免费视频内射| 亚洲精品中文字幕一二三四区 | 纯流量卡能插随身wifi吗| 精品免费久久久久久久清纯 | 亚洲欧洲精品一区二区精品久久久| 久久精品国产a三级三级三级| 永久免费av网站大全| 国产成人免费无遮挡视频| avwww免费| 亚洲人成电影免费在线| 国产精品一二三区在线看| 中文字幕人妻丝袜制服| 91精品国产国语对白视频| 日本欧美视频一区| 精品久久久久久电影网| 午夜视频精品福利| 中文欧美无线码| 老司机亚洲免费影院| 99re6热这里在线精品视频| 国产精品亚洲av一区麻豆| 亚洲三区欧美一区| 天天操日日干夜夜撸| 国产成人a∨麻豆精品| 亚洲视频免费观看视频| 久久久国产一区二区| 午夜视频精品福利| 日韩大码丰满熟妇| 国产精品熟女久久久久浪| 亚洲精品国产精品久久久不卡| 熟女少妇亚洲综合色aaa.| 亚洲精品国产区一区二| 日韩大码丰满熟妇| 精品福利观看| 精品国产一区二区三区久久久樱花| 国产亚洲午夜精品一区二区久久| 丁香六月天网| 久久精品久久久久久噜噜老黄| 久久国产精品影院| 99国产综合亚洲精品| 黑人欧美特级aaaaaa片| 美女午夜性视频免费| 精品熟女少妇八av免费久了| 老司机在亚洲福利影院| 精品久久久久久久毛片微露脸 | 巨乳人妻的诱惑在线观看| 在线观看免费午夜福利视频| 欧美日韩黄片免| 国产成人精品无人区| 亚洲av成人一区二区三| 国产精品九九99| 国产片内射在线| 午夜影院在线不卡| 亚洲伊人久久精品综合| 动漫黄色视频在线观看| 伊人亚洲综合成人网| 欧美日韩国产mv在线观看视频| 久久av网站| 亚洲avbb在线观看| 丝袜喷水一区| 国产成人免费观看mmmm| 99久久人妻综合| 久久九九热精品免费| 五月天丁香电影| 黑丝袜美女国产一区| 久久久久久久久免费视频了| 黄片小视频在线播放| 91麻豆av在线| 少妇精品久久久久久久| 亚洲国产av影院在线观看| 性色av一级| 母亲3免费完整高清在线观看| netflix在线观看网站| 国产精品二区激情视频| 久久ye,这里只有精品| 亚洲欧美色中文字幕在线| 精品国产乱码久久久久久小说| 亚洲三区欧美一区| 51午夜福利影视在线观看| 国产成人免费无遮挡视频| 午夜福利,免费看| 高清av免费在线| 国产在线一区二区三区精| 中文欧美无线码| 97精品久久久久久久久久精品| 中文字幕人妻丝袜制服| 日韩免费高清中文字幕av| www.av在线官网国产| 欧美久久黑人一区二区| 人妻一区二区av| 国产亚洲午夜精品一区二区久久| av欧美777| 最新在线观看一区二区三区| 男女午夜视频在线观看| 最新在线观看一区二区三区| 日本91视频免费播放| 欧美国产精品一级二级三级| 亚洲国产看品久久| 一级片'在线观看视频| 午夜免费成人在线视频| 99久久精品国产亚洲精品| 亚洲精品国产区一区二| 亚洲五月婷婷丁香| 国产亚洲一区二区精品| 在线观看免费日韩欧美大片| 亚洲精品在线美女| 国产一区二区三区av在线| 国产日韩欧美视频二区| 免费人妻精品一区二区三区视频| 亚洲精品第二区| 国产亚洲精品久久久久5区| 99久久99久久久精品蜜桃| 国产在线免费精品| 午夜免费鲁丝| 亚洲成av片中文字幕在线观看| svipshipincom国产片| 国产日韩欧美在线精品| 欧美一级毛片孕妇| 青春草视频在线免费观看| 精品久久蜜臀av无| 亚洲欧美色中文字幕在线| 久久精品久久久久久噜噜老黄| 热re99久久精品国产66热6| 免费看十八禁软件| 手机成人av网站| 满18在线观看网站| 手机成人av网站| 在线看a的网站| 99九九在线精品视频| 午夜福利免费观看在线| 精品少妇内射三级| 五月开心婷婷网| 久久99热这里只频精品6学生| 午夜免费观看性视频| 各种免费的搞黄视频| 91成人精品电影| 国产成人免费观看mmmm| 欧美av亚洲av综合av国产av| 婷婷丁香在线五月| 日韩视频在线欧美| av网站在线播放免费| 亚洲成人免费av在线播放| 夜夜夜夜夜久久久久| 国产日韩欧美亚洲二区| 可以免费在线观看a视频的电影网站| 国产亚洲av片在线观看秒播厂| 国产麻豆69| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 飞空精品影院首页| 最黄视频免费看| 婷婷成人精品国产| 午夜精品久久久久久毛片777| 999久久久精品免费观看国产| 久久九九热精品免费| 热re99久久精品国产66热6| av超薄肉色丝袜交足视频| 欧美国产精品va在线观看不卡| 国产男女超爽视频在线观看| 精品国产乱子伦一区二区三区 | 精品高清国产在线一区| 成人亚洲精品一区在线观看| 国产在视频线精品| 精品免费久久久久久久清纯 | 久久精品国产综合久久久| 精品久久久久久久毛片微露脸 | 黄色视频不卡| 99国产精品99久久久久| 国内毛片毛片毛片毛片毛片| 亚洲一区中文字幕在线| 人成视频在线观看免费观看| 日韩 亚洲 欧美在线| 视频在线观看一区二区三区| 婷婷色av中文字幕| 三级毛片av免费| 好男人电影高清在线观看| 久久这里只有精品19| 国产深夜福利视频在线观看| 日本vs欧美在线观看视频| 在线观看人妻少妇| 亚洲国产欧美在线一区| 午夜免费观看性视频| 久久久久精品国产欧美久久久 | 2018国产大陆天天弄谢| 日韩精品免费视频一区二区三区| 中国国产av一级| 日本精品一区二区三区蜜桃| 啦啦啦中文免费视频观看日本| 五月开心婷婷网| 国产成人a∨麻豆精品| 日韩视频一区二区在线观看| 五月天丁香电影| 亚洲国产欧美一区二区综合| 久久久久久免费高清国产稀缺| 69精品国产乱码久久久| 99国产极品粉嫩在线观看| 亚洲av男天堂| 亚洲精品久久成人aⅴ小说| 夫妻午夜视频| 久久久精品免费免费高清| 男女国产视频网站| 日本黄色日本黄色录像| 色播在线永久视频| 少妇裸体淫交视频免费看高清 | 国产淫语在线视频| 高清av免费在线| 国产又色又爽无遮挡免| 人人妻人人添人人爽欧美一区卜| 90打野战视频偷拍视频| 欧美精品啪啪一区二区三区 | 久久精品亚洲熟妇少妇任你| 亚洲一区二区三区欧美精品| 12—13女人毛片做爰片一| 在线观看免费午夜福利视频| 免费黄频网站在线观看国产| 多毛熟女@视频| 久久av网站| 亚洲精品av麻豆狂野| 色94色欧美一区二区| 高清欧美精品videossex| av网站免费在线观看视频| 一级黄色大片毛片| 亚洲人成电影免费在线| av在线老鸭窝| 久久ye,这里只有精品| 捣出白浆h1v1| 亚洲精品国产av成人精品| 啪啪无遮挡十八禁网站| 国产精品久久久久久精品古装| 女警被强在线播放| 日本欧美视频一区| 久久国产精品人妻蜜桃| 国产一区二区激情短视频 | 18在线观看网站| 免费在线观看影片大全网站| 十分钟在线观看高清视频www| av欧美777| 黑人欧美特级aaaaaa片| 国产一区二区三区在线臀色熟女 | 国产欧美日韩一区二区精品| 久久精品国产a三级三级三级| 亚洲av成人不卡在线观看播放网 | 国产av一区二区精品久久| 国产精品一区二区在线不卡| www.av在线官网国产| av电影中文网址| 久久久久久久久久久久大奶| 丝袜在线中文字幕| 成人国语在线视频| 80岁老熟妇乱子伦牲交| 亚洲精品av麻豆狂野| 中文字幕人妻丝袜制服| 国产日韩欧美视频二区| 国产麻豆69| 精品视频人人做人人爽| 久久久久久久精品精品| 俄罗斯特黄特色一大片| 黄色视频在线播放观看不卡| 别揉我奶头~嗯~啊~动态视频 | av天堂久久9| av有码第一页| 亚洲一区中文字幕在线| 久久精品亚洲av国产电影网| 十八禁高潮呻吟视频| 大片电影免费在线观看免费| www.999成人在线观看| 搡老熟女国产l中国老女人| 母亲3免费完整高清在线观看| 国产成人欧美| 欧美黄色淫秽网站| 亚洲成人国产一区在线观看| 午夜日韩欧美国产| 交换朋友夫妻互换小说| 婷婷色av中文字幕| 一本色道久久久久久精品综合| 一区二区日韩欧美中文字幕| 中文字幕av电影在线播放| 男女之事视频高清在线观看| 啦啦啦啦在线视频资源| 首页视频小说图片口味搜索| 欧美激情 高清一区二区三区| 精品少妇内射三级| 亚洲精华国产精华精| 三上悠亚av全集在线观看| 久久久久国产一级毛片高清牌| 亚洲,欧美精品.| 精品卡一卡二卡四卡免费| 动漫黄色视频在线观看| 99久久精品国产亚洲精品| 80岁老熟妇乱子伦牲交| 91精品伊人久久大香线蕉| 男女下面插进去视频免费观看| 久久久久国产精品人妻一区二区| 99久久国产精品久久久| 日韩欧美一区视频在线观看| 青青草视频在线视频观看| 免费高清在线观看日韩| 成年人午夜在线观看视频| 欧美日韩视频精品一区| 国产日韩欧美在线精品| 日本av手机在线免费观看| 少妇裸体淫交视频免费看高清 | 亚洲黑人精品在线| 老汉色∧v一级毛片| 中文字幕色久视频| 人妻一区二区av| 国产在视频线精品| 久久性视频一级片| 欧美精品高潮呻吟av久久| 国产精品一区二区免费欧美 | 欧美日韩精品网址| 一进一出抽搐动态| 精品一品国产午夜福利视频| 国产1区2区3区精品| 啪啪无遮挡十八禁网站| 91字幕亚洲| 夫妻午夜视频| 精品久久久久久久毛片微露脸 | 高清黄色对白视频在线免费看| 美女视频免费永久观看网站| 中文精品一卡2卡3卡4更新| 亚洲男人天堂网一区| 无限看片的www在线观看| 久久精品亚洲av国产电影网| 国产成人系列免费观看| xxxhd国产人妻xxx| 永久免费av网站大全| 国产欧美日韩一区二区三区在线| √禁漫天堂资源中文www| 大香蕉久久网| 亚洲专区国产一区二区| 色婷婷av一区二区三区视频| 精品乱码久久久久久99久播| 免费在线观看视频国产中文字幕亚洲 | 99久久综合免费| 亚洲av日韩在线播放| 老汉色∧v一级毛片| 人妻久久中文字幕网| 日本黄色日本黄色录像| 9热在线视频观看99| 黄色片一级片一级黄色片| 午夜福利一区二区在线看| netflix在线观看网站| 超色免费av| 欧美av亚洲av综合av国产av| 啦啦啦中文免费视频观看日本| 精品福利永久在线观看| 老熟女久久久| 久久国产精品人妻蜜桃| 狠狠狠狠99中文字幕| 纯流量卡能插随身wifi吗| 夫妻午夜视频| 最近中文字幕2019免费版| 在线观看www视频免费| 婷婷丁香在线五月| 性高湖久久久久久久久免费观看| 国产精品1区2区在线观看. | 精品免费久久久久久久清纯 | 韩国高清视频一区二区三区| 中文字幕色久视频| 岛国毛片在线播放| 97精品久久久久久久久久精品| 亚洲av男天堂| 精品人妻熟女毛片av久久网站| 欧美性长视频在线观看| 久久久国产成人免费| 久久久精品免费免费高清| 欧美乱码精品一区二区三区| 后天国语完整版免费观看| 久久精品熟女亚洲av麻豆精品| 久久 成人 亚洲| 久久久久国内视频| 国产精品麻豆人妻色哟哟久久| 国产av国产精品国产| 老鸭窝网址在线观看| 极品少妇高潮喷水抽搐| 高潮久久久久久久久久久不卡| 国产精品久久久久久人妻精品电影 | 天天影视国产精品| 真人做人爱边吃奶动态| 黄片大片在线免费观看| 亚洲国产看品久久| netflix在线观看网站| tocl精华| 少妇被粗大的猛进出69影院| 99九九在线精品视频| 成人黄色视频免费在线看| 日韩精品免费视频一区二区三区| 精品卡一卡二卡四卡免费| 如日韩欧美国产精品一区二区三区| 成年动漫av网址| 咕卡用的链子| 午夜精品久久久久久毛片777| 他把我摸到了高潮在线观看 | 欧美在线黄色| 色播在线永久视频| 国产一区二区三区av在线| 两个人免费观看高清视频| 又黄又粗又硬又大视频| 国产国语露脸激情在线看| 美女高潮喷水抽搐中文字幕|