• <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)證
    免费看美女性在线毛片视频| 在线免费观看的www视频| 99热只有精品国产| 欧美在线一区亚洲| 美女 人体艺术 gogo| 欧美黄色淫秽网站| 真实男女啪啪啪动态图| 国产成人av激情在线播放| 欧美高清成人免费视频www| 亚洲国产精品合色在线| 男人舔奶头视频| 亚洲中文av在线| 91久久精品国产一区二区成人 | 国产欧美日韩精品亚洲av| 国产精品av久久久久免费| 一个人观看的视频www高清免费观看 | 一个人看的www免费观看视频| 网址你懂的国产日韩在线| 国产免费男女视频| 天堂av国产一区二区熟女人妻| 免费看光身美女| 色老头精品视频在线观看| 熟女电影av网| 久久亚洲精品不卡| 日本免费a在线| 热99re8久久精品国产| 亚洲人成网站高清观看| 免费无遮挡裸体视频| 无遮挡黄片免费观看| 久久久久久大精品| 可以在线观看的亚洲视频| 亚洲人成网站高清观看| 悠悠久久av| 国产精品,欧美在线| 听说在线观看完整版免费高清| 国产淫片久久久久久久久 | 久久精品国产亚洲av香蕉五月| 国产成年人精品一区二区| 国内毛片毛片毛片毛片毛片| 久久九九热精品免费| 两人在一起打扑克的视频| 国产伦精品一区二区三区四那| 欧美日韩瑟瑟在线播放| 亚洲真实伦在线观看| 久久精品国产清高在天天线| 国产精品久久久久久亚洲av鲁大| 国产1区2区3区精品| 国产精品久久久久久久电影 | 精品国产乱码久久久久久男人| 欧美成人一区二区免费高清观看 | 免费看a级黄色片| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区激情短视频| 午夜影院日韩av| 日本 av在线| 日韩欧美精品v在线| 久久精品国产亚洲av香蕉五月| 极品教师在线免费播放| 久久久久国内视频| 久久热在线av| 国产精品99久久99久久久不卡| 国产三级中文精品| 亚洲精品粉嫩美女一区| 黄频高清免费视频| 国产69精品久久久久777片 | 久久人人精品亚洲av| 午夜精品久久久久久毛片777| 国产一区二区激情短视频| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲美女久久久| 国产爱豆传媒在线观看| 老汉色av国产亚洲站长工具| 黄片大片在线免费观看| 中文字幕最新亚洲高清| av福利片在线观看| 午夜a级毛片| 国产精品av久久久久免费| 一级作爱视频免费观看| 久久热在线av| 亚洲电影在线观看av| 色综合亚洲欧美另类图片| av在线天堂中文字幕| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 变态另类成人亚洲欧美熟女| 午夜a级毛片| 日本免费a在线| 热99re8久久精品国产| 精品一区二区三区av网在线观看| 全区人妻精品视频| 97超视频在线观看视频| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 巨乳人妻的诱惑在线观看| 亚洲七黄色美女视频| 老司机午夜福利在线观看视频| 在线观看舔阴道视频| 午夜影院日韩av| 日韩欧美免费精品| 亚洲欧洲精品一区二区精品久久久| 一级黄色大片毛片| 色综合亚洲欧美另类图片| 高清毛片免费观看视频网站| 日日夜夜操网爽| 久久久久久九九精品二区国产| 亚洲中文av在线| 免费在线观看亚洲国产| 欧美不卡视频在线免费观看| 国产美女午夜福利| 国产亚洲精品一区二区www| 亚洲自拍偷在线| 日本a在线网址| 亚洲第一欧美日韩一区二区三区| 又紧又爽又黄一区二区| 久久久水蜜桃国产精品网| 在线观看免费午夜福利视频| 久久精品亚洲精品国产色婷小说| 岛国视频午夜一区免费看| av女优亚洲男人天堂 | 国产伦精品一区二区三区视频9 | av片东京热男人的天堂| 欧美不卡视频在线免费观看| 亚洲av免费在线观看| 日本一二三区视频观看| 少妇的逼水好多| 嫩草影视91久久| 国产黄色小视频在线观看| 精品久久久久久成人av| 国产视频一区二区在线看| 免费高清视频大片| 亚洲欧美激情综合另类| 免费无遮挡裸体视频| 美女黄网站色视频| 97超级碰碰碰精品色视频在线观看| 国产精品99久久久久久久久| 久久精品aⅴ一区二区三区四区| 狠狠狠狠99中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 动漫黄色视频在线观看| 俄罗斯特黄特色一大片| 国产av一区在线观看免费| 日本成人三级电影网站| 久久久久久久午夜电影| av片东京热男人的天堂| www日本在线高清视频| 国内揄拍国产精品人妻在线| 99精品久久久久人妻精品| 黑人操中国人逼视频| 五月伊人婷婷丁香| 国产欧美日韩一区二区精品| 两性午夜刺激爽爽歪歪视频在线观看| 99精品欧美一区二区三区四区| 熟女人妻精品中文字幕| 欧美激情久久久久久爽电影| 国产高清视频在线观看网站| 两人在一起打扑克的视频| 手机成人av网站| 免费电影在线观看免费观看| av片东京热男人的天堂| 一区二区三区激情视频| 欧美性猛交╳xxx乱大交人| 成人三级做爰电影| 亚洲成人免费电影在线观看| a级毛片a级免费在线| 在线看三级毛片| 色老头精品视频在线观看| 免费电影在线观看免费观看| 久久久国产成人精品二区| 成年版毛片免费区| 一二三四社区在线视频社区8| 国产成人一区二区三区免费视频网站| 国产精品1区2区在线观看.| 亚洲av免费在线观看| 中文字幕精品亚洲无线码一区| 久久久久免费精品人妻一区二区| 1024香蕉在线观看| 在线看三级毛片| 老司机在亚洲福利影院| 欧美一区二区国产精品久久精品| 高清在线国产一区| 久久精品人妻少妇| 99国产综合亚洲精品| 久久久久国产一级毛片高清牌| 一级黄色大片毛片| 极品教师在线免费播放| 亚洲欧洲精品一区二区精品久久久| 国产精品综合久久久久久久免费| 亚洲狠狠婷婷综合久久图片| 国产aⅴ精品一区二区三区波| 国产高清视频在线播放一区| 在线观看美女被高潮喷水网站 | 91久久精品国产一区二区成人 | 国产精品99久久久久久久久| 亚洲片人在线观看| 午夜激情福利司机影院| 日本熟妇午夜| www.999成人在线观看| tocl精华| 中文资源天堂在线| 日本精品一区二区三区蜜桃| 美女 人体艺术 gogo| 99国产精品99久久久久| 99re在线观看精品视频| 国产亚洲精品久久久com| 99热这里只有是精品50| 成年版毛片免费区| 国产亚洲欧美在线一区二区| 精品久久蜜臀av无| 欧美激情久久久久久爽电影| 亚洲精品色激情综合| 好男人在线观看高清免费视频| 91老司机精品| 日本与韩国留学比较| 一夜夜www| 少妇的丰满在线观看| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 麻豆成人午夜福利视频| 国产精品99久久久久久久久| 久久中文看片网| 久久久国产成人免费| 黑人欧美特级aaaaaa片| 国产一区在线观看成人免费| 给我免费播放毛片高清在线观看| 免费一级毛片在线播放高清视频| 亚洲片人在线观看| 操出白浆在线播放| www日本在线高清视频| АⅤ资源中文在线天堂| 国产亚洲欧美98| 成人高潮视频无遮挡免费网站| 一a级毛片在线观看| 69av精品久久久久久| 伦理电影免费视频| 97人妻精品一区二区三区麻豆| 免费电影在线观看免费观看| 在线观看免费午夜福利视频| 亚洲人成网站在线播放欧美日韩| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 国产欧美日韩一区二区三| 免费大片18禁| 国产精品野战在线观看| 午夜a级毛片| 国产麻豆成人av免费视频| netflix在线观看网站| 国产精品久久久av美女十八| 午夜激情福利司机影院| 国产一区二区在线观看日韩 | 在线看三级毛片| 特级一级黄色大片| 女生性感内裤真人,穿戴方法视频| 精品无人区乱码1区二区| 国产亚洲精品av在线| 亚洲av免费在线观看| 很黄的视频免费| 岛国在线观看网站| 精品久久久久久久毛片微露脸| 精品无人区乱码1区二区| 亚洲最大成人中文| 久久99热这里只有精品18| 欧美中文日本在线观看视频| 国产99白浆流出| 九九热线精品视视频播放| 中文字幕最新亚洲高清| 久久欧美精品欧美久久欧美| 欧美xxxx黑人xx丫x性爽| 亚洲国产日韩欧美精品在线观看 | 亚洲无线在线观看| 欧美中文综合在线视频| 一个人看的www免费观看视频| 国产亚洲欧美在线一区二区| 精品久久久久久久久久久久久| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交╳xxx乱大交人| 午夜成年电影在线免费观看| 一进一出抽搐gif免费好疼| 俄罗斯特黄特色一大片| 久久人妻av系列| 999久久久国产精品视频| 一个人免费在线观看电影 | 国产av不卡久久| 国产真实乱freesex| 在线观看舔阴道视频| 男女床上黄色一级片免费看| 人妻久久中文字幕网| 亚洲电影在线观看av| 两个人视频免费观看高清| 真实男女啪啪啪动态图| 午夜影院日韩av| 成年女人永久免费观看视频| 精品无人区乱码1区二区| 曰老女人黄片| a级毛片a级免费在线| 精品一区二区三区av网在线观看| 在线视频色国产色| 两性夫妻黄色片| 在线国产一区二区在线| 亚洲,欧美精品.| 亚洲人成伊人成综合网2020| 婷婷亚洲欧美| 成人特级av手机在线观看| 在线观看午夜福利视频| 亚洲午夜理论影院| 亚洲av免费在线观看| 中文字幕精品亚洲无线码一区| 精品一区二区三区视频在线观看免费| 两个人视频免费观看高清| 久久久精品大字幕| 国产私拍福利视频在线观看| 少妇的丰满在线观看| 久久香蕉精品热| 亚洲国产精品成人综合色| 怎么达到女性高潮| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类 | 在线永久观看黄色视频| 久99久视频精品免费| 国产精品国产高清国产av| 亚洲精品美女久久久久99蜜臀| 非洲黑人性xxxx精品又粗又长| 动漫黄色视频在线观看| 夜夜爽天天搞| 99久久精品一区二区三区| 97超视频在线观看视频| 成人亚洲精品av一区二区| 一区二区三区高清视频在线| 后天国语完整版免费观看| 岛国在线观看网站| 嫩草影院精品99| 精品国产乱子伦一区二区三区| 91av网一区二区| 亚洲色图 男人天堂 中文字幕| 丁香六月欧美| 一进一出好大好爽视频| 国产日本99.免费观看| 人妻夜夜爽99麻豆av| 亚洲av熟女| 天堂网av新在线| 五月伊人婷婷丁香| 在线看三级毛片| 亚洲 国产 在线| a级毛片在线看网站| 国产成人精品久久二区二区免费| 少妇丰满av| 国产精品美女特级片免费视频播放器 | 日本黄大片高清| 久久久久久久午夜电影| 天堂av国产一区二区熟女人妻| 在线永久观看黄色视频| aaaaa片日本免费| 欧美乱码精品一区二区三区| 亚洲,欧美精品.| 在线免费观看不下载黄p国产 | 免费在线观看视频国产中文字幕亚洲| 亚洲人成网站在线播放欧美日韩| 免费搜索国产男女视频| 又黄又粗又硬又大视频| 中文字幕熟女人妻在线| 精品久久久久久久久久免费视频| 亚洲自拍偷在线| 欧美日韩一级在线毛片| av天堂中文字幕网| 听说在线观看完整版免费高清| 国产aⅴ精品一区二区三区波| 国产精品综合久久久久久久免费| 日韩av在线大香蕉| 又黄又粗又硬又大视频| 亚洲精品美女久久久久99蜜臀| 日本三级黄在线观看| 18禁观看日本| 琪琪午夜伦伦电影理论片6080| 午夜激情欧美在线| 99国产极品粉嫩在线观看| 午夜免费激情av| 国产视频一区二区在线看| а√天堂www在线а√下载| 久久久久精品国产欧美久久久| 久久99热这里只有精品18| 成年免费大片在线观看| 国产精品综合久久久久久久免费| 国产成人一区二区三区免费视频网站| 人妻夜夜爽99麻豆av| 亚洲狠狠婷婷综合久久图片| 欧美zozozo另类| 九色国产91popny在线| 五月伊人婷婷丁香| 亚洲第一电影网av| 国产精品,欧美在线| 久久国产精品人妻蜜桃| 日韩三级视频一区二区三区| 美女被艹到高潮喷水动态| 国产午夜精品久久久久久| 欧美中文综合在线视频| 999久久久国产精品视频| 国产成人精品无人区| 一级毛片精品| 男女床上黄色一级片免费看| 两个人看的免费小视频| 国产av不卡久久| 一级毛片高清免费大全| 舔av片在线| 中文字幕最新亚洲高清| 99精品欧美一区二区三区四区| 亚洲avbb在线观看| 久久亚洲精品不卡| 亚洲精品456在线播放app | 精品免费久久久久久久清纯| 久久精品亚洲精品国产色婷小说| 深夜精品福利| 久久午夜综合久久蜜桃| 久久亚洲精品不卡| 久久人妻av系列| 国产一区在线观看成人免费| 一卡2卡三卡四卡精品乱码亚洲| 久久精品人妻少妇| 99久久无色码亚洲精品果冻| 国产精品98久久久久久宅男小说| cao死你这个sao货| 熟女少妇亚洲综合色aaa.| 日韩欧美国产一区二区入口| 成年人黄色毛片网站| 99久久久亚洲精品蜜臀av| 亚洲国产欧美一区二区综合| 日本成人三级电影网站| 久久天躁狠狠躁夜夜2o2o| 欧美日本视频| а√天堂www在线а√下载| 国产成人精品久久二区二区免费| 亚洲精品在线观看二区| a级毛片在线看网站| 国内精品一区二区在线观看| 变态另类丝袜制服| 99久久99久久久精品蜜桃| 久久婷婷人人爽人人干人人爱| 色视频www国产| 国产一区在线观看成人免费| 国产日本99.免费观看| 一区二区三区高清视频在线| 亚洲成a人片在线一区二区| 久久久久久人人人人人| 欧美成狂野欧美在线观看| 国产精品综合久久久久久久免费| 成人av在线播放网站| 国产视频一区二区在线看| a在线观看视频网站| 亚洲专区国产一区二区| 两性夫妻黄色片| 亚洲国产高清在线一区二区三| 丝袜人妻中文字幕| 在线观看免费视频日本深夜| 成人无遮挡网站| 亚洲国产高清在线一区二区三| 又大又爽又粗| av福利片在线观看| 亚洲人成电影免费在线| 国产午夜福利久久久久久| 国产成人一区二区三区免费视频网站| 每晚都被弄得嗷嗷叫到高潮| 桃色一区二区三区在线观看| 免费看a级黄色片| 午夜激情欧美在线| 久久久国产欧美日韩av| 国产精品女同一区二区软件 | av天堂在线播放| 一进一出抽搐动态| 美女午夜性视频免费| 首页视频小说图片口味搜索| 精品国产三级普通话版| 美女扒开内裤让男人捅视频| 草草在线视频免费看| 成人18禁在线播放| 久久久水蜜桃国产精品网| 99热精品在线国产| 日韩欧美免费精品| 成人国产综合亚洲| 亚洲精品中文字幕一二三四区| 一二三四在线观看免费中文在| 啪啪无遮挡十八禁网站| 欧美日韩黄片免| 成人av在线播放网站| 色哟哟哟哟哟哟| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线观看二区| www.精华液| 最新中文字幕久久久久 | 日韩大尺度精品在线看网址| 哪里可以看免费的av片| 操出白浆在线播放| 美女午夜性视频免费| 人妻久久中文字幕网| 欧美黑人欧美精品刺激| 人妻丰满熟妇av一区二区三区| 国产私拍福利视频在线观看| 国产欧美日韩精品一区二区| 亚洲18禁久久av| 欧美高清成人免费视频www| or卡值多少钱| 国产精品国产高清国产av| 亚洲精华国产精华精| 国产精品一区二区精品视频观看| 精品熟女少妇八av免费久了| 此物有八面人人有两片| 国产1区2区3区精品| 1024香蕉在线观看| 国产亚洲精品av在线| 久久天躁狠狠躁夜夜2o2o| 淫妇啪啪啪对白视频| 哪里可以看免费的av片| 一级毛片高清免费大全| 亚洲av电影在线进入| 特级一级黄色大片| 毛片女人毛片| 久久久国产欧美日韩av| 观看免费一级毛片| 久久精品国产清高在天天线| 中文字幕高清在线视频| 91久久精品国产一区二区成人 | 成人一区二区视频在线观看| 淫秽高清视频在线观看| 麻豆成人av在线观看| 在线观看美女被高潮喷水网站 | 97超视频在线观看视频| 亚洲人成网站高清观看| 亚洲精品粉嫩美女一区| 老鸭窝网址在线观看| 日韩人妻高清精品专区| 法律面前人人平等表现在哪些方面| 五月伊人婷婷丁香| 国产欧美日韩精品一区二区| 老鸭窝网址在线观看| 国产精品久久久久久久电影 | 国产亚洲欧美98| 欧美午夜高清在线| 岛国在线免费视频观看| 老熟妇仑乱视频hdxx| 色av中文字幕| 精品国产美女av久久久久小说| 中文字幕av在线有码专区| 国产成人福利小说| 国产伦在线观看视频一区| 午夜免费激情av| 免费无遮挡裸体视频| 黄色片一级片一级黄色片| 久久99热这里只有精品18| 久久久久久久久免费视频了| 国产黄色小视频在线观看| 黄色 视频免费看| 无人区码免费观看不卡| 欧美午夜高清在线| 国产午夜精品论理片| 午夜亚洲福利在线播放| 色哟哟哟哟哟哟| 一卡2卡三卡四卡精品乱码亚洲| 久久亚洲精品不卡| 狂野欧美白嫩少妇大欣赏| 久久久国产欧美日韩av| 一边摸一边抽搐一进一小说| 又大又爽又粗| 免费看光身美女| 久久精品国产亚洲av香蕉五月| 麻豆av在线久日| 九九在线视频观看精品| 99国产综合亚洲精品| 757午夜福利合集在线观看| 色精品久久人妻99蜜桃| 法律面前人人平等表现在哪些方面| 国产又色又爽无遮挡免费看| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美一区二区三区在线观看| 色视频www国产| 中文字幕久久专区| 曰老女人黄片| 精品无人区乱码1区二区| 又粗又爽又猛毛片免费看| 久久伊人香网站| 国产又色又爽无遮挡免费看| 国内久久婷婷六月综合欲色啪| a级毛片a级免费在线| 两人在一起打扑克的视频| 国产高清视频在线观看网站| 日日摸夜夜添夜夜添小说| 国产成年人精品一区二区| 淫妇啪啪啪对白视频| 久久中文看片网| 国产精品野战在线观看| 亚洲最大成人中文| 男女视频在线观看网站免费| 亚洲中文字幕日韩| 成熟少妇高潮喷水视频| 美女黄网站色视频| 国产亚洲精品av在线| 最近最新中文字幕大全免费视频| 国产精品99久久久久久久久| 视频区欧美日本亚洲| 男人和女人高潮做爰伦理| 大型黄色视频在线免费观看| 欧美3d第一页| 99久久久亚洲精品蜜臀av| a在线观看视频网站| 一二三四在线观看免费中文在| 欧美3d第一页| 国产激情欧美一区二区| 丝袜人妻中文字幕| www.自偷自拍.com| 亚洲欧美精品综合久久99| 深夜精品福利| 91九色精品人成在线观看| 欧美成人性av电影在线观看| 999久久久精品免费观看国产| 亚洲国产精品成人综合色| 欧美日本亚洲视频在线播放| 最新在线观看一区二区三区| 91字幕亚洲|