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

    水下爆炸柱型裝藥與球形裝藥遠(yuǎn)場(chǎng)等效關(guān)系

    2017-10-13 02:03:25張弛宇劉榮忠楊永亮
    關(guān)鍵詞:遠(yuǎn)場(chǎng)長(zhǎng)徑裝藥

    張弛宇, 郭 銳, 劉榮忠, 陳 亮, 楊永亮

    ?

    水下爆炸柱型裝藥與球形裝藥遠(yuǎn)場(chǎng)等效關(guān)系

    張弛宇, 郭 銳, 劉榮忠, 陳 亮, 楊永亮

    (南京理工大學(xué)機(jī)械工程學(xué)院, 江蘇南京, 210094)

    柱型戰(zhàn)斗部和球形戰(zhàn)斗部作為2種最常規(guī)戰(zhàn)斗部形式, 其水下爆炸性能是其研究的重點(diǎn), 目前已存的水下爆炸經(jīng)驗(yàn)公式皆為球形裝藥。為研究柱型裝藥水下爆炸壓力場(chǎng)是否也存在類似球形場(chǎng)的經(jīng)驗(yàn)公式, 文中首先定性分析了作用方位和長(zhǎng)徑比對(duì)柱型裝藥水下爆炸壓力場(chǎng)的影響。通過(guò)對(duì)比分析柱型裝藥和等藥量球形裝藥壓力場(chǎng), 求得柱型裝藥與球形裝藥遠(yuǎn)場(chǎng)沖擊波關(guān)系式方程, 并結(jié)合球形裝藥經(jīng)驗(yàn)公式, 推導(dǎo)得到柱型裝藥壓力場(chǎng)近似經(jīng)驗(yàn)公式。仿真試驗(yàn)驗(yàn)證了上述公式的準(zhǔn)確性。文中的研究可為優(yōu)化水下戰(zhàn)斗部爆炸性能提供依據(jù)。

    水下爆炸;柱型裝藥;球形裝藥;經(jīng)驗(yàn)公式

    0 引言

    水下爆炸是水下武器及水下防護(hù)的基本課題,研究方向主要為炸藥的爆轟研究、沖擊波的產(chǎn)生和傳播以及氣泡脈動(dòng)的研究。球型裝藥作為戰(zhàn)斗部一種主要形式, 研究比較廣泛, 且目前已存的經(jīng)驗(yàn)公式皆為球形裝藥。而柱型裝藥由于其長(zhǎng)徑比變化多, 很難通過(guò)單一的研究判斷其壓力場(chǎng)規(guī)律。Sternberg[1]用2D軸對(duì)稱Lagrange方法計(jì)算了不同長(zhǎng)徑比的柱型裝藥水中爆炸壓力分布, 表明沖擊波參數(shù)受長(zhǎng)徑比控制, 且指出長(zhǎng)徑比接近1時(shí), 沖擊波可以近似為球形。Hammond[2]對(duì)柱型炸藥的研究表明, 雖然裝藥形狀對(duì)水中爆炸沖擊波有很大影響, 但超過(guò)一定距離時(shí), 柱型炸藥沖擊波場(chǎng)可以用等質(zhì)量球形炸藥近似, 與Sternberg研究相同。侯俊亮[3]通過(guò)數(shù)值仿真研究了不同形狀裝藥爆炸沖擊波場(chǎng)及對(duì)靶板作用效應(yīng)。李金河[4]等人用試驗(yàn)方法研究了柱形裝藥水下爆炸軸向和徑向沖擊波峰值變化, 表明其仍然符合相似律和指數(shù)衰減規(guī)律。但上述研究均未得出柱型裝藥沖擊波壓力場(chǎng)經(jīng)驗(yàn)公式。

    文中基于數(shù)值仿真, 通過(guò)對(duì)比、分析固定長(zhǎng)徑比柱型TNT裝藥與球形裝藥的沖擊波壓力場(chǎng), 得出一種具有普遍意義的柱型裝藥和球形裝藥遠(yuǎn)場(chǎng)沖擊波等效關(guān)系, 并結(jié)合球形裝藥經(jīng)驗(yàn)公式, 推導(dǎo)出柱型裝藥壓力場(chǎng)近似經(jīng)驗(yàn)公式, 為后續(xù)科學(xué)研究提供了理論依據(jù)。

    1 水下爆炸Colo經(jīng)驗(yàn)公式

    目前關(guān)于水下爆炸沖擊波最常用的公式多采用Colo總結(jié)的經(jīng)驗(yàn)公式[5], 具體如下

    (2)

    (3)

    式中: P為沖擊波壓力峰值壓力;P為沖擊波壓力;為衰減系數(shù);為爆距;0為藥包半徑;為TNT質(zhì)量。

    2 仿真模型和狀態(tài)方程

    2.1 炸藥的狀態(tài)方程

    炸藥采用JWL方程

    式中:為壓力;為相對(duì)體積;為炸藥內(nèi)能;,均為材料參數(shù);1,2和為常數(shù)。

    TNT狀態(tài)方程參數(shù)見(jiàn)表1。表中:為炸藥密度;為爆熱。

    表1 TNT狀態(tài)方程參數(shù)

    2.2 水的狀態(tài)方程

    水采用SHOCK狀態(tài)方程

    2.3 計(jì)算模型

    圖1為柱型裝藥有限元模型[6], 水域的大小為2 400×1 200, 裝藥為29 gTNT裝藥, 模型關(guān)于軸中心對(duì)稱, 水域和炸藥均采用歐拉網(wǎng)格。

    邊界條件設(shè)置為Flow-out。定義中心起爆方式, 炸藥中心點(diǎn)為坐標(biāo)原點(diǎn)(0,0), 炸藥軸向(水平方向軸)為0°, 炸藥徑向(豎直方向軸)為90°, 同一半徑上相鄰的監(jiān)測(cè)點(diǎn)間距22.5°, 相鄰半徑上同一方向角監(jiān)測(cè)點(diǎn)間距250 mm, 共設(shè)置20個(gè)監(jiān)測(cè)點(diǎn)。

    圖2為球形裝藥有限元模型[7], 采用1D楔形模型, 邊界條件設(shè)為全流出邊界條件, 楔形長(zhǎng)度為3 000 mm。

    3 仿真結(jié)果與分析

    3.1 柱型裝藥和球形裝藥等效關(guān)系式

    Sternberg[1]指出柱型裝藥沖擊波峰值壓力與爆距、作用方向及長(zhǎng)徑比()有關(guān)。設(shè)柱型裝藥沖擊波峰值壓力為(,,)。球形裝藥沖擊波峰值壓力為()。

    以=13/3的柱型TNT裝藥為代表, 根據(jù)AUTODYN仿真得出裝藥在不同爆距和方向角(,)下的沖擊波峰值壓力(MPa), 如表2所示。

    表2 柱型裝藥(L/D=13/3)不同爆距方向角時(shí)沖擊波峰值壓力

    表3為球形TNT裝藥在不同爆距的沖擊波峰值壓力大小。

    表3 球形裝藥不同爆距時(shí)的沖擊波峰值壓力

    對(duì)比表2和表3可知, 此布置下的柱形裝藥在徑向(90°方向)產(chǎn)生的沖擊波峰值壓力大于等質(zhì)量同位置的球形裝藥, 徑向峰值壓力最少提高了約10%, 在裝藥軸向(0°)同距離處沖擊波峰值壓力小于球形裝藥, 故推測(cè)爆炸遠(yuǎn)場(chǎng)存在一點(diǎn), 使得此處柱型裝藥和球形裝藥爆炸沖擊波峰值壓力相同, 在此可用柱型沖擊波完全代替球形[8]。

    圖3為1/4柱型裝藥(/=13/3)和球形裝藥爆炸云圖??梢钥闯? 球形裝藥沖擊波以球形傳播, 當(dāng)爆距較小時(shí), 柱型裝藥沖擊波云圖近似為矩形, 隨爆距增大, 趨于球形, 當(dāng)爆距達(dá)到一定距離時(shí), 繼續(xù)以球形傳播形式傳播[5]。這是由于水的可壓縮性很小, 當(dāng)爆距較小時(shí), 裝藥形狀對(duì)爆炸影響很大, 在爆距較遠(yuǎn)時(shí), 裝藥形狀幾乎不構(gòu)成影響[9]。

    表4 柱型裝藥(L/D=13/3)和球形裝藥沖擊波峰值壓力比值

    由圖4可知, 不同爆距下的沖擊波峰值壓力比值隨方向角的增大, 呈現(xiàn)出先分散后重合的特征, 且不同爆距的映射曲線上隨角度的遞增幾乎都交于一點(diǎn), 證明了在爆炸某處, 柱型裝藥在此處產(chǎn)生的沖擊波峰值壓力與球形裝藥完全相同, 驗(yàn)證了以上推測(cè)[10]。

    以相同的試驗(yàn)方法, 分別得到/=7/3時(shí)和/=19/3時(shí)柱型裝藥的爆炸云圖和壓力比值曲線, 見(jiàn)圖5。

    圖6和圖7分別為2種柱型裝藥與球形裝藥沖擊波峰值壓力比值曲線走勢(shì)圖。

    對(duì)比圖4~圖7可知, 不同長(zhǎng)徑比的柱型裝藥爆炸近場(chǎng)均近似以柱型傳播, 遠(yuǎn)場(chǎng)以球形傳播。且隨著長(zhǎng)徑比的增加, 軸向沖擊波峰值壓力減少, 徑向峰值壓力增加, 這是由于柱型裝藥幾何特點(diǎn)造成爆炸沖擊波的疊加, 當(dāng)長(zhǎng)徑比無(wú)窮大時(shí), 裝藥近似線性, 能量幾乎集中在徑向。

    以=13/3的柱型裝藥為例, 根據(jù)沖擊波峰值壓力比值走向, 取傅里葉函數(shù)表示該曲線, 設(shè)

    將表4不同點(diǎn)的比值系數(shù)帶入方程,θ為比值系數(shù)為1時(shí)對(duì)應(yīng)所得角度。通過(guò)擬合, 得到函數(shù)未知數(shù)(,,),見(jiàn)表5。

    表5 不同爆距下方程未知數(shù)擬合結(jié)果

    由表5數(shù)據(jù)知, 通過(guò)擬合得(0.25 m, 46.2o)、(0.5 m, 49.698o)、(0.75 m, 48.859o)、(1.0 m, 48.976o)處柱型裝藥產(chǎn)生的沖擊波峰值壓力與球形裝藥完全相同。

    為了得到柱型裝藥在不同(,)的與球形裝藥的關(guān)系, 進(jìn)一步分析表5數(shù)據(jù)。未知數(shù)穩(wěn)定在2.23附近, 取=2.23;也相對(duì)穩(wěn)定, 取為48o;在公式中影響因素不超過(guò)3%, 可忽略不計(jì)。因此, 沖擊波峰值比值走向方程

    (8)

    故, 在極坐標(biāo)下, 柱型(/=13/3)裝藥與球形裝藥沖擊波峰值關(guān)系式

    為了驗(yàn)證柱型裝藥與球形裝藥關(guān)系式的準(zhǔn)確性, 將柱型裝藥和球形裝藥擬合公式結(jié)果與球形裝藥結(jié)果進(jìn)行對(duì)比。

    表6為根據(jù)式(6)~式(9)計(jì)算所得不同爆距方向角(,)時(shí)的映射系數(shù)(,)值。

    表6 不同爆距方向角的映射值

    通過(guò)仿真驗(yàn)證, 由柱型裝藥和球形裝藥關(guān)系所得出的數(shù)據(jù)與AUTODYN所得出數(shù)據(jù)除(,)在(0.2 m, 0)處誤差超過(guò)5%, 其余點(diǎn)均小于5%,見(jiàn)表7, 原因在于: 1) 根據(jù)圖3映射系數(shù)走勢(shì), 方向角較小時(shí), 映射曲線分散大, 因此按照假設(shè)擬合產(chǎn)生的誤差也大; 2) 在大長(zhǎng)徑比情況下, 近距離的軸向和徑向所受長(zhǎng)徑比的影響較大。

    表7 柱形裝藥和球形裝藥等效關(guān)系式驗(yàn)證結(jié)果

    注:

    以相同的方法分析/=7/3和/=19/3時(shí)的柱型裝藥也滿足要求, 總的來(lái)說(shuō), 基于仿真得出柱型裝藥和球形裝藥沖擊波峰值關(guān)系式滿足要求, 證明了關(guān)系式的準(zhǔn)確性。

    3.2 柱型裝藥遠(yuǎn)場(chǎng)壓力場(chǎng)分析

    為得出柱型裝藥遠(yuǎn)場(chǎng)沖擊波壓力場(chǎng)經(jīng)驗(yàn)公式, 將柱型裝藥和球形裝藥等效關(guān)系式與Cole經(jīng)驗(yàn)公式對(duì)比, 見(jiàn)表8。表中,TNT為仿真數(shù)據(jù);P為經(jīng)驗(yàn)值計(jì)算數(shù)據(jù);為擬合值, 且

    因此

    (11)

    表8 仿真數(shù)據(jù)與經(jīng)驗(yàn)公式數(shù)據(jù)

    柱型裝藥水下爆炸遠(yuǎn)場(chǎng)壓力場(chǎng)分布

    4 結(jié)論

    為研究柱型裝藥水下爆炸壓力場(chǎng)是否也存在類似球形場(chǎng)的經(jīng)驗(yàn)公式, 文章定性分析了作用方位和長(zhǎng)徑比對(duì)柱型裝藥水下爆炸壓力場(chǎng)的影響, 并得出以下結(jié)論。

    1) 圓柱裝藥在極坐標(biāo)(,)下存在點(diǎn), 使得該點(diǎn)處柱型裝藥沖擊波峰值等于球形裝藥, 該點(diǎn)在一定長(zhǎng)徑比范圍內(nèi)穩(wěn)定于壓力場(chǎng)50°附近處。

    2) 基于數(shù)值仿真, 得出柱型裝藥與球形裝藥等效關(guān)系式, 結(jié)合Cole經(jīng)驗(yàn)公式推導(dǎo)得出了柱型裝藥壓力場(chǎng)分布公式均是較為準(zhǔn)確的。

    3) 隨著長(zhǎng)徑比的增加, 軸向沖擊波峰值壓力減少, 徑向峰值壓力增加。柱形裝藥在徑向(90°方向)產(chǎn)生的沖擊波峰值壓力大于等質(zhì)量同位置的球形裝藥, 徑向峰值壓力最少提高了約4%, 根據(jù)此特點(diǎn)可以柱型裝藥長(zhǎng)徑比優(yōu)化。

    4) 在一定范圍內(nèi), 隨著長(zhǎng)徑比的增加, 裝藥形狀對(duì)沖擊波的影響距離增加, 沖擊波傳播形式仍有近似柱型傳播變?yōu)榍蛐蝹鞑ァ?/p>

    文中所得柱型裝藥和球形裝藥遠(yuǎn)場(chǎng)爆炸壓力場(chǎng)等效關(guān)系和柱型裝藥遠(yuǎn)場(chǎng)壓力場(chǎng)計(jì)算公式, 已通過(guò)仿真初步驗(yàn)證其準(zhǔn)確性, 后續(xù)的研究將側(cè)重于通過(guò)試驗(yàn)進(jìn)一步驗(yàn)證理論的正確性。

    [1] Stemberg H M. Underwater Detonation of Pentolite Cylinders[J]. Physics of Fluids, 1995, 30(3): 761-769.

    [2] Hammond L. Underwater Shock Wave Characteristics of Cylindrical Charges[R]. Defence Science and Technology Organisation. AMRL, DSTO-GD-0029, 1995.

    [3] 侯俊亮, 蔣建偉, 門建兵, 等. 不同形狀裝藥爆炸沖擊波場(chǎng)及對(duì)靶板作用效應(yīng)的數(shù)值模擬[J]. 北京理工大學(xué)學(xué)報(bào), 2013, 33(6): 556-561.Hou Jun-liang, Jiang Jian-wei, Men Jian-bing, et al. Numerieal Simulation on Blast Wave Field and Deformation of Thin Plate Under Different-shape Charge Loading[J]. Transactions of Beijing Institute of Technology, 2013, 33(6): 556-561.

    [4] 李金河, 趙繼波, 池家春. 水中爆炸沖擊波傳播規(guī)律的實(shí)驗(yàn)研究[J]. 高能量密度物理, 2007, 22(4): 363-367.

    [5] Cole R H. Underwater Explosions[M]. New Jersey: Princeton University Press, 1948.

    [6] 賈憲振, 胡毅亭, 董明榮, 等. 基于ANSYS/LS-DYNA模擬水下爆炸沖擊波的等效質(zhì)量法[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2008, 28(3): 159-162.Jia Xian-Zhen, Hu Yi-ting, Dong Ming-rong, et al. Equivalent Mass Method of Underwater Explosionshock Wave Simulation Based on ANSYS/LS-DYNA[J]. Journal of Projectiles Rockets Missiles and Guidance, 2008, 28(3): 159-162.

    [7] 肖秋平, 陳網(wǎng)樺, 賈憲振, 等. 基于AUTODYN的水下爆炸沖擊波模擬研究[J]. 艦船科學(xué)技術(shù), 2009, 31(2): 38-43.Xiao Qiu-ping, Chen Wang-hua, Jia Xian-zhen, et al. Numerical Study of Underwater Explosion Shock Wave Based on Autodyn[J]. Ship Science and Technology, 2009, 31(2): 38-43.

    [8] 趙繼波, 譚多望, 李金河, 等. 柱形裝藥水中爆炸近場(chǎng)徑向壓力測(cè)試初探[J]. 高壓物理學(xué)報(bào), 2008, 22(3): 323-328.Zhao Ji-bo, Tan Duo-wang, Li Jin-he, et al. Primary Research on Side Pressure of Cylindrical TNT at Underwater Exposive Close-Field[J]. Chinese Journal of High Pressure Physics, 2008, 22(3): 323-328.

    [9] Li J, Rong J L. Experimental and Numerical Investigation of the Dynamic Response of Structuressubjected to Underwater Explosion[J]. European Journal of Mechanics B/Fluids, 2012, 32(3): 59-69.

    [10] Zhang A-man, Yang W S, Huang C, et al. Numerical Sim- ulation of Column Charge Underwater Explosion Based on SPH and BEM combination[J]. Computers & Fluids, 2013, 71: 169-178.

    (責(zé)任編輯: 楊力軍)

    Equivalent Relationship between Cylindrical Charge and Spherical Charge for Underwater Explosion

    ZHANG Chi-yuGUO RuiLIU Rong-zhongCHEN LiangYANG Yong-liang

    (School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China)

    Both cylindrical charge and spherical charge warheads are commonly used, so their underwater explosion performances attract more attention. However, the existing empirical formulas are only for spherical charge. To understand if the underwater explosion pressure field of cylindrical charge can be described by the empirical formulas for spherical charge, This paper qualitatively analyzes the influences of the orientation and the ratio of length to diameter on pressure field of cylinder charge underwater explosion, obtains a relationship expression of far field shock waves of cylinder charge and cycle charge by contrasting the pressure fields of cylindrical charge and same amount of spherical charge, and deduce an approximate empirical formula for cylinder charge by considering the empirical formulas for spherical charge. In addition, simulation is conducted to verify the accuracy of the deduced formula. This research may provide the basis for optimizing explosion performance of underwater warhead.

    underwater explosion; cylindrical charge; spherical charge; empirical formula

    10.11993/j.issn.1673-1948.2017.01.0013

    TJ630; TQ560.1

    A

    1673-1948(2017)01-0065-06

    2016-10-12;

    2016-11-04.

    高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金項(xiàng)目(20133219110019).

    張弛宇(1991-), 男, 在讀碩士, 主要研究方向?yàn)閺椝幘_化與智能化.

    猜你喜歡
    遠(yuǎn)場(chǎng)長(zhǎng)徑裝藥
    基于全三維動(dòng)網(wǎng)格技術(shù)的變長(zhǎng)徑比間隙環(huán)流的研究
    火炸藥學(xué)報(bào)(2022年3期)2022-07-04 07:31:00
    孔內(nèi)爆炸填塞效應(yīng)的數(shù)值模擬與分析
    玄武巖纖維長(zhǎng)徑比對(duì)混凝土力學(xué)性能的影響
    基于隨形冷卻的大長(zhǎng)徑比筆套注塑優(yōu)化
    某發(fā)射裝藥結(jié)構(gòu)改進(jìn)設(shè)計(jì)的新思路
    基于仿真與實(shí)測(cè)的列車遠(yuǎn)場(chǎng)氣動(dòng)噪聲分析
    某種陣列雷達(dá)發(fā)射通道遠(yuǎn)場(chǎng)校準(zhǔn)簡(jiǎn)易方法
    銦摻雜調(diào)控氧化鋅納米棒長(zhǎng)徑比
    深孔不耦合裝藥爆破技術(shù)卸壓效果驗(yàn)證
    成人三级黄色视频| 亚洲精品久久成人aⅴ小说| 午夜精品国产一区二区电影| 一级a爱片免费观看的视频| 色综合亚洲欧美另类图片| 性色av乱码一区二区三区2| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器 | 嫩草影院精品99| 电影成人av| 中文字幕色久视频| 女性被躁到高潮视频| 涩涩av久久男人的天堂| 性色av乱码一区二区三区2| 99久久精品国产亚洲精品| 我的亚洲天堂| 国内精品久久久久精免费| 日韩欧美免费精品| 国产91精品成人一区二区三区| 亚洲av成人不卡在线观看播放网| 99国产精品一区二区三区| 无限看片的www在线观看| av免费在线观看网站| 中文亚洲av片在线观看爽| 亚洲一区中文字幕在线| 国产亚洲精品久久久久久毛片| 老汉色∧v一级毛片| 成人三级黄色视频| АⅤ资源中文在线天堂| av在线天堂中文字幕| 欧美日韩亚洲综合一区二区三区_| videosex国产| 免费高清视频大片| 岛国视频午夜一区免费看| 亚洲国产精品999在线| 国产一区在线观看成人免费| 波多野结衣高清无吗| 国产在线观看jvid| 亚洲国产毛片av蜜桃av| 久热这里只有精品99| 亚洲成人国产一区在线观看| 欧美激情 高清一区二区三区| 深夜精品福利| 国产一区二区激情短视频| cao死你这个sao货| 人人妻人人澡人人看| 在线播放国产精品三级| 在线观看免费视频日本深夜| 中文字幕人妻丝袜一区二区| 国产成人av教育| 看黄色毛片网站| 伊人久久大香线蕉亚洲五| 欧美激情极品国产一区二区三区| 精品无人区乱码1区二区| 麻豆一二三区av精品| 老熟妇仑乱视频hdxx| 亚洲精品美女久久av网站| 如日韩欧美国产精品一区二区三区| 两个人看的免费小视频| 身体一侧抽搐| 久久久精品欧美日韩精品| x7x7x7水蜜桃| 欧美精品啪啪一区二区三区| 日韩欧美一区二区三区在线观看| 欧美乱妇无乱码| 亚洲在线自拍视频| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 日韩一卡2卡3卡4卡2021年| 99国产精品99久久久久| 免费看a级黄色片| videosex国产| 午夜免费激情av| 精品无人区乱码1区二区| 精品日产1卡2卡| 又黄又爽又免费观看的视频| 91麻豆精品激情在线观看国产| 嫁个100分男人电影在线观看| 欧美久久黑人一区二区| 999精品在线视频| 长腿黑丝高跟| av片东京热男人的天堂| 国产精品亚洲美女久久久| 搡老岳熟女国产| 最近最新中文字幕大全免费视频| 国产成人啪精品午夜网站| 少妇 在线观看| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 国产精品亚洲美女久久久| 一级a爱片免费观看的视频| 久久久久久久精品吃奶| 亚洲国产精品久久男人天堂| 亚洲成人久久性| 久久久久久国产a免费观看| 99riav亚洲国产免费| av片东京热男人的天堂| 不卡av一区二区三区| 男女下面进入的视频免费午夜 | or卡值多少钱| 熟女少妇亚洲综合色aaa.| 9热在线视频观看99| 99re在线观看精品视频| 国产精品久久久人人做人人爽| a级毛片在线看网站| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| 高清黄色对白视频在线免费看| √禁漫天堂资源中文www| 一边摸一边做爽爽视频免费| 国产精品爽爽va在线观看网站 | 久久伊人香网站| 久久香蕉精品热| 久久久水蜜桃国产精品网| 别揉我奶头~嗯~啊~动态视频| 我的亚洲天堂| 日本黄色视频三级网站网址| 久久精品国产亚洲av香蕉五月| 少妇裸体淫交视频免费看高清 | 又黄又爽又免费观看的视频| 久久人人精品亚洲av| 欧美黄色片欧美黄色片| 90打野战视频偷拍视频| 精品久久久久久久人妻蜜臀av | 麻豆久久精品国产亚洲av| 国产xxxxx性猛交| 亚洲精品粉嫩美女一区| 久久亚洲真实| 久久久久国产一级毛片高清牌| 国产精品一区二区免费欧美| 精品少妇一区二区三区视频日本电影| 男女做爰动态图高潮gif福利片 | 一区在线观看完整版| 校园春色视频在线观看| 免费观看精品视频网站| 亚洲午夜精品一区,二区,三区| 久久香蕉精品热| 精品久久久久久,| 久久精品国产亚洲av香蕉五月| 成人亚洲精品一区在线观看| 一级,二级,三级黄色视频| 性少妇av在线| 久久午夜综合久久蜜桃| 亚洲自偷自拍图片 自拍| 免费在线观看影片大全网站| 久久国产精品男人的天堂亚洲| 亚洲熟妇熟女久久| 欧美在线一区亚洲| 级片在线观看| 91成人精品电影| 亚洲国产精品成人综合色| 国产精品一区二区在线不卡| 久久人人97超碰香蕉20202| а√天堂www在线а√下载| 多毛熟女@视频| 999精品在线视频| 免费高清在线观看日韩| 天堂影院成人在线观看| 亚洲自拍偷在线| 一级片免费观看大全| 我的亚洲天堂| 久久亚洲精品不卡| 免费在线观看黄色视频的| 欧美乱色亚洲激情| АⅤ资源中文在线天堂| 欧美国产日韩亚洲一区| 午夜免费成人在线视频| 丝袜美腿诱惑在线| 99久久精品国产亚洲精品| 日韩大码丰满熟妇| 国产色视频综合| 国产激情欧美一区二区| 老司机在亚洲福利影院| 日韩欧美国产一区二区入口| 国产精品爽爽va在线观看网站 | 波多野结衣高清无吗| 亚洲电影在线观看av| 免费在线观看完整版高清| 久久久久久久久久久久大奶| 欧美丝袜亚洲另类 | 久久中文字幕一级| 欧美日韩福利视频一区二区| 99精品在免费线老司机午夜| 动漫黄色视频在线观看| 亚洲男人天堂网一区| 黄片大片在线免费观看| 欧美大码av| 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 欧美另类亚洲清纯唯美| av天堂在线播放| 午夜福利在线观看吧| 成人欧美大片| 亚洲激情在线av| 久久久精品国产亚洲av高清涩受| 亚洲三区欧美一区| 大型av网站在线播放| 中文字幕久久专区| 黄色女人牲交| aaaaa片日本免费| 在线观看66精品国产| 国产成人精品在线电影| 在线播放国产精品三级| 精品人妻1区二区| 90打野战视频偷拍视频| 亚洲欧美一区二区三区黑人| 搡老岳熟女国产| 90打野战视频偷拍视频| 好男人在线观看高清免费视频 | 亚洲精品国产色婷婷电影| 九色亚洲精品在线播放| 免费少妇av软件| 搡老妇女老女人老熟妇| 久久久国产成人精品二区| 免费一级毛片在线播放高清视频 | 国产高清有码在线观看视频 | 一个人免费在线观看的高清视频| 欧美丝袜亚洲另类 | 国产精品美女特级片免费视频播放器 | 国产一区二区三区综合在线观看| 人妻久久中文字幕网| 12—13女人毛片做爰片一| 满18在线观看网站| 国产亚洲精品一区二区www| 宅男免费午夜| 欧洲精品卡2卡3卡4卡5卡区| 一进一出抽搐gif免费好疼| videosex国产| 国产成人精品久久二区二区91| 变态另类丝袜制服| 国产熟女午夜一区二区三区| 91精品国产国语对白视频| 嫩草影视91久久| 国产精品av久久久久免费| 首页视频小说图片口味搜索| www日本在线高清视频| 国产精品久久久久久精品电影 | 久久久久久大精品| 一区在线观看完整版| 亚洲第一电影网av| 超碰成人久久| 99re在线观看精品视频| 国产成人精品久久二区二区91| 大码成人一级视频| 国产免费男女视频| 亚洲自拍偷在线| 最近最新中文字幕大全免费视频| 很黄的视频免费| 亚洲av片天天在线观看| 人人妻,人人澡人人爽秒播| 欧美日本中文国产一区发布| 久9热在线精品视频| 日韩一卡2卡3卡4卡2021年| 啪啪无遮挡十八禁网站| 狂野欧美激情性xxxx| 日日爽夜夜爽网站| 身体一侧抽搐| 深夜精品福利| 久久久精品国产亚洲av高清涩受| 久久精品91蜜桃| 男女之事视频高清在线观看| 亚洲av片天天在线观看| 日本精品一区二区三区蜜桃| 法律面前人人平等表现在哪些方面| 国产精品久久久久久亚洲av鲁大| 国产亚洲av嫩草精品影院| 亚洲一区中文字幕在线| 51午夜福利影视在线观看| 丁香六月欧美| 国产精品免费一区二区三区在线| www.精华液| 亚洲av成人不卡在线观看播放网| 黄频高清免费视频| 九色国产91popny在线| 伊人久久大香线蕉亚洲五| 亚洲av第一区精品v没综合| 高清在线国产一区| 亚洲三区欧美一区| 男女午夜视频在线观看| 欧美在线一区亚洲| 国产人伦9x9x在线观看| 免费在线观看视频国产中文字幕亚洲| 免费少妇av软件| 日本a在线网址| 国产精品久久电影中文字幕| 夜夜看夜夜爽夜夜摸| 欧美午夜高清在线| 中文字幕人妻熟女乱码| 咕卡用的链子| 日本五十路高清| av在线天堂中文字幕| 中文字幕色久视频| 国产精品99久久99久久久不卡| 久久狼人影院| 丝袜人妻中文字幕| 成熟少妇高潮喷水视频| 一级a爱视频在线免费观看| 男女床上黄色一级片免费看| 国内久久婷婷六月综合欲色啪| 无遮挡黄片免费观看| av福利片在线| 少妇粗大呻吟视频| 99国产精品一区二区蜜桃av| 国产精品一区二区免费欧美| 波多野结衣av一区二区av| 久久精品91无色码中文字幕| 久久影院123| 国产亚洲av高清不卡| 黄频高清免费视频| 超碰成人久久| 狠狠狠狠99中文字幕| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 成人手机av| 欧美日本视频| 久久精品国产亚洲av香蕉五月| 两性夫妻黄色片| 成人精品一区二区免费| 久久精品91蜜桃| 成人亚洲精品一区在线观看| 97超级碰碰碰精品色视频在线观看| 女性被躁到高潮视频| 国产三级在线视频| 久久精品91无色码中文字幕| 日韩大尺度精品在线看网址 | 日日干狠狠操夜夜爽| 国产精品一区二区在线不卡| 老汉色∧v一级毛片| 亚洲九九香蕉| 成人三级做爰电影| 激情视频va一区二区三区| 日本精品一区二区三区蜜桃| 成人国语在线视频| 精品少妇一区二区三区视频日本电影| 亚洲aⅴ乱码一区二区在线播放 | 欧美乱色亚洲激情| 久久精品国产亚洲av高清一级| 免费一级毛片在线播放高清视频 | 不卡av一区二区三区| 俄罗斯特黄特色一大片| 99久久99久久久精品蜜桃| 视频区欧美日本亚洲| 看免费av毛片| 欧美日本中文国产一区发布| 国产成人免费无遮挡视频| 十八禁网站免费在线| av天堂久久9| 国产精品久久久人人做人人爽| √禁漫天堂资源中文www| 日韩精品免费视频一区二区三区| 91麻豆av在线| 亚洲精品久久成人aⅴ小说| 精品久久久久久成人av| 男女做爰动态图高潮gif福利片 | 丁香欧美五月| 欧美日韩亚洲国产一区二区在线观看| 久热这里只有精品99| 韩国av一区二区三区四区| 美女高潮喷水抽搐中文字幕| 亚洲 欧美一区二区三区| 久久精品影院6| 精品一区二区三区av网在线观看| 啦啦啦 在线观看视频| 亚洲色图综合在线观看| 久久久久久大精品| 亚洲专区字幕在线| cao死你这个sao货| 亚洲五月色婷婷综合| 日本 欧美在线| 国产精品一区二区精品视频观看| 国产成人影院久久av| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 亚洲 欧美 日韩 在线 免费| 男人操女人黄网站| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 国产男靠女视频免费网站| aaaaa片日本免费| 亚洲全国av大片| 一区二区三区精品91| 国产亚洲av高清不卡| 亚洲人成电影观看| 亚洲国产欧美日韩在线播放| 久久久精品国产亚洲av高清涩受| 午夜福利18| 成年女人毛片免费观看观看9| 国产精品久久久av美女十八| 亚洲精品国产一区二区精华液| 日本 av在线| 欧美国产日韩亚洲一区| 午夜福利,免费看| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区| 欧美日本视频| 国产成人精品久久二区二区91| 欧美不卡视频在线免费观看 | 性色av乱码一区二区三区2| 欧美性长视频在线观看| 亚洲黑人精品在线| 夜夜看夜夜爽夜夜摸| 一本久久中文字幕| 成人特级黄色片久久久久久久| 黄片大片在线免费观看| 国产精品综合久久久久久久免费 | 亚洲精品在线美女| 日韩有码中文字幕| 日本五十路高清| 久久久国产成人精品二区| 国产xxxxx性猛交| √禁漫天堂资源中文www| 久久精品国产清高在天天线| 精品国内亚洲2022精品成人| 91成年电影在线观看| av电影中文网址| 一夜夜www| 免费无遮挡裸体视频| 麻豆av在线久日| 动漫黄色视频在线观看| 国产成人av教育| 999久久久精品免费观看国产| cao死你这个sao货| 亚洲av五月六月丁香网| 亚洲五月婷婷丁香| 午夜福利高清视频| 国产成人欧美在线观看| 一区福利在线观看| 久久草成人影院| 老熟妇仑乱视频hdxx| 免费在线观看黄色视频的| 伊人久久大香线蕉亚洲五| 禁无遮挡网站| 国产精品亚洲一级av第二区| tocl精华| 波多野结衣高清无吗| 国产单亲对白刺激| 免费av毛片视频| 人妻丰满熟妇av一区二区三区| 成人免费观看视频高清| 国产精品乱码一区二三区的特点 | 欧美一区二区精品小视频在线| 国产极品粉嫩免费观看在线| 成人三级黄色视频| 级片在线观看| 久久久久久国产a免费观看| 好看av亚洲va欧美ⅴa在| 51午夜福利影视在线观看| 一二三四在线观看免费中文在| 午夜福利视频1000在线观看 | 大码成人一级视频| 国产私拍福利视频在线观看| 一级毛片精品| 亚洲欧美一区二区三区黑人| 少妇 在线观看| 后天国语完整版免费观看| 日韩精品青青久久久久久| 可以免费在线观看a视频的电影网站| 日韩精品青青久久久久久| 69精品国产乱码久久久| 午夜福利18| 久久伊人香网站| 国产精品免费一区二区三区在线| 99国产精品99久久久久| 男女下面插进去视频免费观看| 女同久久另类99精品国产91| 一级作爱视频免费观看| 亚洲国产毛片av蜜桃av| 最新美女视频免费是黄的| 亚洲中文av在线| 一本大道久久a久久精品| 一级,二级,三级黄色视频| 亚洲三区欧美一区| 免费在线观看亚洲国产| 亚洲片人在线观看| 免费观看精品视频网站| 国产亚洲精品久久久久5区| 国产精品99久久99久久久不卡| 精品欧美国产一区二区三| 韩国精品一区二区三区| 色av中文字幕| 欧美黑人欧美精品刺激| 国产亚洲精品av在线| cao死你这个sao货| 免费在线观看影片大全网站| 欧美中文日本在线观看视频| 一进一出好大好爽视频| 99国产精品99久久久久| 色播亚洲综合网| 不卡一级毛片| 桃色一区二区三区在线观看| 一进一出好大好爽视频| 午夜日韩欧美国产| 久久亚洲真实| 日韩免费av在线播放| 国产成人欧美在线观看| 99热只有精品国产| 一本综合久久免费| 久久久久九九精品影院| 国产av在哪里看| 日本在线视频免费播放| 女生性感内裤真人,穿戴方法视频| 超碰成人久久| 午夜精品在线福利| 日韩成人在线观看一区二区三区| 精品第一国产精品| 黑丝袜美女国产一区| 亚洲av五月六月丁香网| 精品熟女少妇八av免费久了| 亚洲第一av免费看| 精品国产乱子伦一区二区三区| 叶爱在线成人免费视频播放| 亚洲人成网站在线播放欧美日韩| 性少妇av在线| 18禁国产床啪视频网站| 欧美乱码精品一区二区三区| 久久精品成人免费网站| 人人妻,人人澡人人爽秒播| 国内久久婷婷六月综合欲色啪| 国产三级黄色录像| 色老头精品视频在线观看| 好男人在线观看高清免费视频 | 日日干狠狠操夜夜爽| 久久久久久免费高清国产稀缺| 日本撒尿小便嘘嘘汇集6| 国产免费av片在线观看野外av| 真人做人爱边吃奶动态| 国产精品 国内视频| 精品电影一区二区在线| 日韩欧美在线二视频| av有码第一页| 在线a可以看的网站| 日韩 亚洲 欧美在线| 男女下面进入的视频免费午夜| 欧美最黄视频在线播放免费| 欧美成人性av电影在线观看| 国产高潮美女av| 一卡2卡三卡四卡精品乱码亚洲| АⅤ资源中文在线天堂| 日韩精品中文字幕看吧| 别揉我奶头~嗯~啊~动态视频| 午夜福利欧美成人| 天堂av国产一区二区熟女人妻| 日日干狠狠操夜夜爽| 成人av一区二区三区在线看| 精品99又大又爽又粗少妇毛片 | 久久久久国产精品人妻aⅴ院| 好男人在线观看高清免费视频| 婷婷色综合大香蕉| 国产欧美日韩精品一区二区| 毛片一级片免费看久久久久 | 一本一本综合久久| www日本黄色视频网| 欧美潮喷喷水| 欧美成人a在线观看| 男人舔奶头视频| 99久久成人亚洲精品观看| 国产精品精品国产色婷婷| 夜夜看夜夜爽夜夜摸| 久久久国产成人精品二区| 天堂动漫精品| 欧美日韩瑟瑟在线播放| 国产成人aa在线观看| 日韩欧美一区二区三区在线观看| 韩国av一区二区三区四区| 狂野欧美激情性xxxx在线观看| 成人高潮视频无遮挡免费网站| 久久久久九九精品影院| 日韩人妻高清精品专区| 国产单亲对白刺激| 色综合亚洲欧美另类图片| 啦啦啦韩国在线观看视频| av在线天堂中文字幕| 国产黄片美女视频| 日本一二三区视频观看| 日韩精品青青久久久久久| 99久国产av精品| 无遮挡黄片免费观看| 女的被弄到高潮叫床怎么办 | 非洲黑人性xxxx精品又粗又长| 嫩草影视91久久| 亚洲成a人片在线一区二区| 日本与韩国留学比较| 精品久久久久久成人av| 九色成人免费人妻av| 中国美白少妇内射xxxbb| 午夜福利在线观看吧| 久久精品人妻少妇| 校园春色视频在线观看| 国产黄a三级三级三级人| 亚洲欧美日韩东京热| 午夜免费激情av| 国产精品福利在线免费观看| 精品人妻一区二区三区麻豆 | 人人妻人人澡欧美一区二区| 又爽又黄a免费视频| or卡值多少钱| 国产白丝娇喘喷水9色精品| 1024手机看黄色片| 一级a爱片免费观看的视频| 午夜精品一区二区三区免费看| 天美传媒精品一区二区| 不卡一级毛片| 蜜桃亚洲精品一区二区三区| 久久久久久久久久黄片| 特级一级黄色大片| 国产亚洲av嫩草精品影院| 最好的美女福利视频网| 乱码一卡2卡4卡精品| 欧美一区二区亚洲| 亚洲电影在线观看av| 久久久国产成人精品二区| 天堂动漫精品| 成人永久免费在线观看视频| 亚洲最大成人手机在线| 午夜福利18|