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

    水下爆炸沖擊波在圓柱殼結構表面繞射衰減分布

    2022-07-12 12:41:58張迪洲何鎮(zhèn)宏何心怡
    水下無人系統(tǒng)學報 2022年3期
    關鍵詞:沖量藥量沖擊波

    張迪洲 , 何鎮(zhèn)宏 , 何心怡 , 李 營 , 盧 軍 , 陳 雙

    水下爆炸沖擊波在圓柱殼結構表面繞射衰減分布

    張迪洲1, 何鎮(zhèn)宏2, 何心怡1, 李 營2, 盧 軍1, 陳 雙1

    (1. 海軍研究院, 北京, 100161; 2. 北京理工大學 爆炸科學與技術國家重點實驗室, 北京, 100081)

    水下爆炸沖擊波作用于圓柱殼表面時, 除反射沖擊波外, 在背爆面還會形成繞射沖擊波。為了研究沖擊波壓力在圓柱殼表面的繞射分布特性, 采用ABAQUS 軟件的耦合歐拉-拉格朗日方法對水下爆炸沖擊波與圓柱殼結構的相互作用進行了數(shù)值仿真, 得到了圓柱殼周圍的壓力場分布情況, 分析了不同爆距和藥量對繞射沖擊波衰減的影響。研究結果表明: 沖擊波壓力峰值和沖量在迎爆面區(qū)域急速衰減, 隨著爆距的減小, 衰減到50%的角度變小, 衰減速度明顯變快, 而藥量對沖擊波壓力峰值和沖量衰減影響都很小。

    水下爆炸; 圓柱殼; 沖擊波; 繞射壓力衰減

    0 引言

    圓柱殼結構是艦艇及其他海洋結構物中廣泛使用的工程構件, 世界各國海軍密切關注水下爆炸對上述設施和結構的破壞問題, 其中, 沖擊波載荷是上述結構和設施安全的主要威脅之一, 其與圓柱殼體間的相互作用具有重大研究意義。

    在分析水下爆炸沖擊波波動效應時, 除了圓柱殼迎爆面的壓力分布和反射外, 也要考慮圓柱殼背爆面形成的繞射沖擊波。Huang [1]考慮了球面沖擊波與圓柱殼的相互作用, 計算了殼面幾個點的衍射壓力, 得到解析解。李裕春等[2]研究了水中平面沖擊波對剛性方形臺的繞射效應。姚熊亮等[3]采用ABAQUS 軟件, 仿真分析了雙層圓柱殼結構的抗沖擊性能。朱正洋[4]研究了近地空中爆炸沖擊波與柱形拱殼結構的相互作用, 分析了沖擊波在空氣介質中的傳播過程和對結構的繞射特性。邵宗戰(zhàn)[5]測量了水下爆炸作用下圓柱殼周圍的沖擊波壓力曲線。此外,姚熊亮等[6]還分析了不同板厚和爆距對沖擊波繞射圓柱殼的影響; 譚皓洋[7]采用ABAQUS的耦合歐拉-拉格朗日(coupled Eulerian-Lagrangian, CEL)方法, 仿真了水下爆炸沖擊波、氣泡脈動和氣泡射流全物理過程。但以上研究重點均圍繞水下爆炸圓柱殼載荷變化機理和結構響應規(guī)律, 對實尺度下沖擊波繞射規(guī)律和圖譜關注不夠, 限制了工程應用。

    文中以10 m 級圓柱殼為研究對象, 采用ABAQUS 有限元程序對水下爆炸沖擊波繞圓柱殼傳播的動態(tài)過程進行了數(shù)值仿真, 分析了沖擊波沿圓柱殼繞射的衰減過程, 得到了壓力峰值和沖量分布圖譜, 為水下武器毀傷評估的理論方法提供載荷輸入。

    1 數(shù)值仿真模型

    1.1 物理模型

    水下爆炸是一個復雜的流固耦合問題, 炸藥爆炸會在毫秒時間內產生高峰值壓力的沖擊波, 文中采用CEL 方法模擬沖擊波與圓柱殼的相互作用。計算模型如圖1 所示, 參考實際尺寸, 圓柱殼半徑為R=5 m, 藥包TNT 當量為100 kg, 藥包中心與殼體間的距離r=5 m。圓柱殼內部為空氣, 外部為水, 圓柱殼設為剛性結構并固定。為減小流場邊界影響, 流場大小為圓柱殼截面方向長50 m, 高20 m, 圓柱殼長度方向15 m, 流場最小網(wǎng)格和結構網(wǎng)格尺寸為0.1 m, 流場邊界采用無反射邊界條件。由于圓柱殼為對稱結構, 以1/2 模型進行計算,計圓柱殼迎爆面離藥包最近點角度為0°, 背爆面離藥包最遠點角度為180°, 為研究沖擊波沿圓柱殼的繞射, 在圓柱殼周圍每隔5°設置1 個壓力測點, 圖1 標記了幾個典型角度的測點。

    圖1 水下爆炸物理模型Fig. 1 Physical model of underwater explosion

    1.2 參數(shù)設置

    1.2.1 水介質狀態(tài)方程

    US-UP狀態(tài)方程被用來模擬不可壓的層狀流動, 并通過納維-斯托克斯(Navier-Stokes, N-S)方程控制運動, 能準確模擬水在沖擊波載荷作用下的運動過程。線性US-UPHugoniot 形式狀態(tài)方程為[8]

    式中:US為線性沖擊速度;UP為質點速度; 水中聲速C1=1 500 m/s; 材料常數(shù)S1=1.75, 水密度為1 000 kg/m3。

    1.2.2 空氣介質狀態(tài)方程

    理想氣體狀態(tài)方程為[8]

    式中: 空氣密度ρa=1.2 kg/m3; 氣體常數(shù)R=287;θ為當前溫度; θZ為絕對零度; 環(huán)境壓力pA設為0 .1 Mpa[8]。

    1.2.3 炸藥狀態(tài)方程

    水下爆炸產生的爆轟產物采用JWL(Jones-Wilkins-Lee)狀態(tài)方程

    式中:A,B,R1,R2,ω為試驗得到的材料常數(shù),A=3.7377×1011Pa ,B=3.7471×109Pa,R1=0.45,R2=0.9, ω=0.35; 炸藥單位質量內能Em=3.8×106J/kg,初 始炸藥密度ρ0=1 630 kg/m3; 爆轟產物密度為ρ。

    1.3 有效性驗證

    1.3.1 仿真結果與經驗公式對比

    為驗證數(shù)值仿真爆炸沖擊波載荷的準確性, 對自由場中的水下爆炸沖擊波進行計算。流場網(wǎng)格大小為0.1 m, TNT 當量為100 kg, 此時初始藥包半徑R0=0.245 m, 爆點設在藥包中心, 如圖2(a)所示。在距離爆心10R0、15R0和20R0處設置壓力測點P1,P2和P3。水下爆炸數(shù)值仿真計算的沖擊波壓力與Geersamp;Hunter 模型公式[9]對比曲線如圖2(b)~(d)所示, 從圖中可以看出, 峰值壓力數(shù)值仿真結果與經驗公式吻合良好。壓力曲線在峰值過后的多次壓力震蕩, Cole[10]認為是由于不同位置的炸藥顆粒產生的沖擊波疊加造成的, 而在數(shù)值仿真模型中的不同炸藥填充單元造成了峰值后的震蕩。

    圖2 爆炸沖擊波壓力數(shù)值仿真與經驗公式對比Fig. 2 Comparison of explosion shock wave pressures between numerical simulation and empirical formula

    1.3.2 仿真結果與試驗數(shù)據(jù)對比

    由水下爆炸沖擊波縮比理論[11]可知, 在相同爆距半徑比下, 水下爆炸沖擊波壓力值相同。根據(jù)文獻[12]中的試驗數(shù)據(jù), 試驗壓力傳感器位于TNT炸藥10 倍半徑處, 在相同爆距半徑比下, 水下爆炸沖擊波峰值仿真結果與文獻中的試驗結果如表1所示。從表中可以發(fā)現(xiàn), 仿真值與試驗值誤差較小, 驗證了仿真結果的正確性。

    表1 試驗與仿真沖擊波峰值壓力對比Table 1 Comparison of peak pressures of shock wave between test and simulation

    2 水下爆炸沖擊波繞射特性

    2.1 水下爆炸沖擊波繞射原理

    基于空氣動力學理論, 超聲速氣流沿外折微小角度的壁面流動時, 經過折角后流動空間擴大, 速度增大, 壓強、密度、溫度降低, 壓強變化的擾動沿折點向外傳播形成弱膨脹波, 與波前氣流方向形成夾角, 圓柱殼面可看作連續(xù)的折角變化, 在圓柱殼的每一處都形成弱膨脹波, 從而形成了傳播方向改變的繞射波[13]。

    2.2 水下爆炸沖擊波繞射過程

    數(shù)值仿真計算得到的水域壓力云圖如圖3 所示。炸藥起爆后, 沖擊波以球形向外傳播, 當?shù)竭_圓柱殼后, 形成反射, 如圖3(a)所示; 沖擊波在結構處反射后, 在結構處形成空化, 如圖3(b)所示; 沖擊波繼續(xù)沿著圓柱殼向前傳播, 到達圓柱殼頂部, 同時圓柱殼的反射波向炸藥方向傳播, 經過氣泡時產生截斷現(xiàn)象, 如圖3(c)所示; 沖擊波最后到達圓柱殼背爆面, 上下兩側沖擊波形成疊加效應, 如圖3(d)所示。

    2.3 不同角度下圓柱殼壓力時域特征

    不同角度下的圓柱殼壓力時域曲線如圖4 所示。迎爆面測點在沖擊波到達后出現(xiàn)第1 個峰值,隨后迅速衰減, 由于殼體的反射形成了第2 個波峰, 如圖4(a)所示; 背爆面測點在沖擊波到達后出現(xiàn)波峰, 在第1 個波峰衰減后仍然保持一定壓力,在180°處由于上下兩側沖擊波疊加, 壓力峰值疊加明顯大于170°和175°處, 同時也可以明顯看到,170°和175°處在第1 個峰值還未衰減完成時, 另一側沖擊波引起的2 次峰值疊加, 如圖4(b)所示。

    圖3 水域壓力云圖Fig. 3 Pressure contour of water

    圖4 壓力時域曲線Fig. 4 Time domain curves of pressure

    3 沖擊波繞射影響因素分析

    3.1 不同爆距下沖擊波繞射特性

    近距離爆炸下, 沖擊波繞射曲率與爆距密切相關。為研究爆距對沖擊波繞射衰減的影響, 在不改變藥量100 kg 的條件下, 改變爆距進行仿真計算, 設計爆距范圍為3~10 m。圓柱殼半徑R=5 m,爆 距為r, 定義無量綱相對爆距λ=r/R。

    3.1.1 壓力時域特征

    不同爆距下典型角度處壓力時域曲線如圖5所示, 可見在迎爆面沖擊波隨著爆距減小衰減較為明顯, 而背爆面沖擊波經過繞射后, 不同爆距下沖擊波峰值相差不大, 在150°處沖擊波峰值已經較為接近。

    3.1.2 壓力峰值空間分布

    在每個工況下, 將0°處的沖擊波壓力峰值記為Pmax, 測點處的沖擊波峰值記為Pm, 定義無量綱參數(shù)Pm/Pmax來衡量沖量的衰減程度, 如圖6 所示??梢钥吹經_擊波壓力峰值在迎爆面區(qū)域急速衰減, 而在背爆面衰減速度變慢。同時, 隨著爆距減小, 衰減速度也變快, 在λ=2時約68°處達到了50%的衰減, 在λ=0.6時約40°處便達到了50%的衰減。由于180°處上下兩側沖擊波的疊加, 沖擊波 峰值略微上升。

    3.1.3 沖量空間分布

    沖量是壓力場持續(xù)作用在目標物上壓力的時間積累量, 是造成結構毀傷的一個重要參數(shù), 對沖擊波壓力曲線的第1 個波峰時間段進行積分, 得到沖擊波在單位面積的沖量I。在每個工況下,將0°處的沖量記為Imax, 測點處的沖量記為Im, 定義無量綱參數(shù)Im/Imax, 來衡量沖量的衰減程度,如圖7 所示??梢钥吹經_量在迎爆面區(qū)域急速衰減, 而在背爆面衰減速度變慢。同時, 隨著爆距減小, 衰減速度也變快, 在λ=2時約70°處達到了50%的衰減, 在λ=0.6時約35°處便達到了50%的衰減。由于180°處上下兩側沖擊波的疊加, 沖量略微 上升。

    圖5 不同角度處典型爆距下沖擊波壓力時域曲線Fig. 5 Time domain curves of shock wave pressure with typical blast distance at different angles

    圖6 不同爆距下Pm/Pmax 隨角度變化曲線Fig. 6 Curves of Pm/Pmax varies with the angles in different blast distance

    圖7 不同爆距下Im/Imax 隨角度變化曲線Fig. 7 Curves of Im/I max varies with the angles in different blast distance

    3.2 不同藥量下沖擊波繞射特性

    沖擊波峰值與藥量相關, 為研究藥量對沖擊波繞射衰減的影響, 在不改變爆距5 m 的條件下, 通過改變藥量進行仿真計算, 設計范圍為25~400 kg。

    3.2.1 壓力時域特征

    不同藥量下典型角度處壓力時域曲線如圖8所示, 可見不論在迎爆面還是背爆面, 沖擊波壓力峰值都隨著藥量增大而顯著增大。

    3.2.2 壓力峰值空間分布

    無量綱參數(shù)Pm/Pmax在不同藥量下隨爆距的變化如圖9 所示, 可以看到在不同藥量下, 沖擊波壓力衰減趨勢相似。

    3.2.3 沖量空間分布

    無量綱參數(shù)Im/Imax在不同藥量下隨爆距的變化如圖10 所示, 可以看到在不同藥量下, 沖擊波壓力衰減趨勢相似。

    圖8 不同藥量下典型角度處沖擊波壓力時域曲線Fig. 8 Time domain curves of shock wave pressure with typical angles at different charge of TNT

    圖9 不同藥量下Pm/Pmax 隨角度變化曲線Fig. 9 Curves of Pm/Pmax varies with the angles in different charge of TNT

    圖10 不同藥量下Im/Imax 隨角度變化曲線Fig. 10 Curve of Im/Imax varies with the angles in different charge of TNT

    4 結論

    文中對水下爆炸沖擊波繞圓柱殼傳播的動態(tài)過程進行了數(shù)值仿真, 分析了沖擊波沿圓柱殼繞射的衰減過程, 得到了沖擊波壓力峰值、沖量在圓柱殼表面的分布圖譜, 探討了爆距和藥量對分布圖譜的影響, 為水下武器毀傷評估的理論方法提供載荷輸入。主要結論有:

    1) 水下爆炸沖擊波在傳播到圓柱殼表面時,會產生反射波和繞射波, 迎爆面處壓力迅速達到峰值隨后衰減, 反射波形成第2 個較小的峰值。

    2) 背爆面繞射沖擊波形成峰值, 衰減速度變慢, 在180°背爆面處由于兩側繞射波同時到達, 沖擊波疊加峰值變大。

    3) 沖擊波壓力峰值和沖量在迎爆面區(qū)域急速衰減, 隨著爆距的減小, 衰減到50%的角度變小,衰減速度明顯變快。而藥量對沖擊波壓力峰值和沖量衰減影響都很小。

    [1]H uang H. Scattering of Spherical Pressure Pulses by a Hard Cylinder[J]. Journal of the Acoustical Society of America,1975, 58: 310-317.

    [2]李 裕春, 高振儒, 馮彬, 等. 水中平面沖擊波對剛性方形臺的繞射效應研究[J]. 系統(tǒng)仿真學報, 2016, 28(4): 940-945.

    Li Yu-chun, Gao Zhen-ru, Feng Bin, et al. Study on Diffraction Effects of Planar Underwater Shock Interaction with Rigid Quadrilateral Structure[J]. Journal of System Simulation,2016, 28(4): 940-945.

    [3]姚 熊亮, 梁德利, 許維軍. 雙層殼結構抗沖擊性能仿真研究[J]. 哈爾濱工程大學學報, 2004, 25(3): 267-273.

    Yao Xiong-liang, Liang De-li, Xu Wei-jun. Research on the Simulation of Anti-shock Capability of Double Shell[J]. Journal of Harbin Engineering University, 2004, 25(3): 267-273.

    [4]朱 正洋. 近地空中爆炸沖擊波與柱形拱殼結構相互作用研究[J]. 建筑結構, 2011, 41(S1): 1486-1489.

    Zhu Zheng-yang. Research in the Interaction of Blast Wave with Semicircular Arched Structure Near the Ground Explosion[J]. Building Structure, 2011, 41(S1): 1486-1489.

    [5]邵 宗戰(zhàn). 水下爆炸作用下圓柱殼周圍壓力場分布[J]. 兵器裝備工程學報, 2017, 38(3): 38-41.

    Shao Zong-zhan. Pressure Field Distribution Around Cylindrical Shell Subjected to Underwater Explosion[J]. Journal of Sichuan Ordnance, 2017, 38(3): 38-41.

    [6]姚 熊亮, 屈子悅, 張樂, 等. 水下爆炸圓柱殼沖擊波繞射特性分析[J]. 哈爾濱工程大學學報, 2020, 41(1): 9-16.

    Yao Xiong-liang, Qu Zi-yue, Zhang Le, et al. Analysis of the Underwater Shock Wave Diffraction Characteristics of Cylindrical Shell[J]. Journal of Harbin Engineering University,2020, 41(1): 9-16.

    [7]譚 皓洋. 基于CEL方法的艦船近場水下爆炸全過程數(shù)值模擬研究[D]. 哈爾濱: 哈爾濱工程大學, 2017.

    [8]S imulia DCS. Abaqus 6.11 Analysis User ’s Manual[R].United States: Abaqus, 2011.

    [9]T homas G L. An Integrated Wave-effects Model for an Underwater Explosion Bubble[J]. The Journal of the Acoustical Society of America, 2002, 111: 1584-1601.

    [10]C ole R H. Underwater Explosion[M]. Princeton: New Jersey,1948.

    [11]張 效慈. 水下爆炸試驗相似準則[J]. 船舶力學, 2007, 11(1):108-118.

    Zhang Xiao-ci. Similarity Criteria for Experiment of Underwater Explosion[J]. Journal of Ship Mechanices, 2007, 11(1):108-118.

    [12]孫 遠翔, 田俊宏, 張之凡, 等. 含鋁炸藥近場水下爆炸沖擊波的實驗及數(shù)值模擬[J]. 振動與沖擊, 2020, 39(14): 171-178, 193.

    Sun Yuan-xiang, Tian Jun-hong, Zhang Zhi-fan, et al. Experiment and Numerical Simulation Study on the Near-field Underwater Explosion of Aluminized Explosive[J]. Journal of Vibration and Shock, 2020, 39(14): 171-178, 193.

    [13]屈 子悅. 水下爆炸作用下圓柱殼繞射特性及壓力分布特征研究[D]. 哈爾濱: 哈爾濱工程大學, 2019.

    Diffraction Attenuation Distribution of Underwater Explosion Shock Waves on the Surface of a Cylindrical Structure

    ZHANG Di-zhou1,HE Zhen-hong2,HE Xin-yi1,LI Ying2,LU Jun1,CHEN Shuang1

    (1. Naval Research Institute, Beijing 100161; 2. State Key Laboratory of Explosion Science and Technology, Beijing University of Technology, Beijing 100081)

    When an underwater explosion shock wave acts on the surface of a cylindrical shell, a shock wave forms on the backshock surface in addition to the reflected shock wave. To study the diffraction distribution features of shock wave pressure on the surface of a cylindrical shell, the coupled Eulerian-Lagrangian method in the ABAQUS soft ware is used to simulate the interactions of underwater explosion shock waves and cylindrical shell structures. The pressure field distribution around a cylindrical shell is obtained. The influence of different blast distances and charges of TNT on the diffraction effect of shock wave attenuation is analyzed. The results indicate that the peak pressure and impulse of shock waves decay rapidly in front of the blast surface. With a decrease in blast distance, the angle required for attenuation to 50% decreases and the attenuation speed increases significantly. However, the charge of TNT has little effect on the peak pressure and impulse attenuation of shock waves.

    underwater explosion; cylindrical shell; shock wave; diffraction

    TJ410; U661.44

    A

    2096-3920(2022)03-0371-07

    10.11993/j.issn.2096-3920.2022.03.013

    張迪洲, 何鎮(zhèn)宏, 何心怡, 等. 水下爆炸沖擊波在圓柱殼結構表面繞射衰減分布[J]. 水下無人系統(tǒng)學報, 2022,30(3): 371-377.

    2022-03-29;

    2022-05-08.

    張迪洲(1988-), 男, 博士, 助理研究員, 主要研究方向為裝備作戰(zhàn)效能評估技術.

    (責任編輯: 許 妍)

    猜你喜歡
    沖量藥量沖擊波
    連續(xù)3年銷量翻番,2022年欲沖量4000萬,福建這家動保企業(yè)正強勢崛起
    當代水產(2022年3期)2022-04-26 14:27:18
    碳化對飛灰穩(wěn)定化加藥量測評的影響
    兩倍藥量
    武漢沖擊波
    中國公路(2019年10期)2019-06-28 03:05:08
    能源物聯(lián)網(wǎng)沖擊波
    能源(2018年10期)2018-12-08 08:02:34
    醫(yī)生集團沖擊波
    小型水池抗爆炸沖擊極限藥量分析
    藥量適中
    小說月刊(2015年6期)2015-12-16 14:55:45
    超聲雙探頭聯(lián)合定位法在體外沖擊波碎石術中的應用
    變力沖量的教學研究
    物理通報(2011年8期)2011-01-24 07:39:38
    日本色播在线视频| 女人十人毛片免费观看3o分钟| 精品国内亚洲2022精品成人| 日韩av在线大香蕉| 成人欧美大片| 看片在线看免费视频| 校园春色视频在线观看| 欧美丝袜亚洲另类| 亚洲久久久久久中文字幕| 麻豆乱淫一区二区| 日本黄色片子视频| 在线免费观看不下载黄p国产| 在现免费观看毛片| 国产成人影院久久av| 久久综合国产亚洲精品| 男人舔女人下体高潮全视频| 欧美丝袜亚洲另类| 亚洲国产色片| 精品午夜福利在线看| 日韩视频在线欧美| 欧美极品一区二区三区四区| 在线国产一区二区在线| 久99久视频精品免费| 免费观看的影片在线观看| 亚洲av第一区精品v没综合| 青春草亚洲视频在线观看| 变态另类丝袜制服| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国模一区二区三区四区视频| 在线免费观看的www视频| 99国产极品粉嫩在线观看| 内射极品少妇av片p| 精品国内亚洲2022精品成人| 99久久成人亚洲精品观看| 亚洲美女搞黄在线观看| 成人鲁丝片一二三区免费| 国产精品嫩草影院av在线观看| 狠狠狠狠99中文字幕| 精品久久国产蜜桃| 国产一级毛片七仙女欲春2| 国产精品电影一区二区三区| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久精品电影小说 | 男人的好看免费观看在线视频| 国国产精品蜜臀av免费| 久久人人精品亚洲av| 2022亚洲国产成人精品| 日本免费一区二区三区高清不卡| 欧美精品一区二区大全| 小说图片视频综合网站| 麻豆久久精品国产亚洲av| av在线观看视频网站免费| 亚洲乱码一区二区免费版| 中文亚洲av片在线观看爽| 国产伦在线观看视频一区| av黄色大香蕉| 久久久久久久久久久丰满| 久久久久性生活片| 特大巨黑吊av在线直播| 日本三级黄在线观看| 免费在线观看成人毛片| 亚洲精品国产av成人精品| 1024手机看黄色片| 色播亚洲综合网| 最近2019中文字幕mv第一页| 国产真实乱freesex| 国产视频内射| 国产伦理片在线播放av一区 | 可以在线观看毛片的网站| 日日摸夜夜添夜夜添av毛片| 三级经典国产精品| 亚洲精品亚洲一区二区| 国产一区二区三区av在线 | 中文字幕人妻熟人妻熟丝袜美| 成人美女网站在线观看视频| 亚洲三级黄色毛片| 午夜视频国产福利| 亚洲经典国产精华液单| 三级国产精品欧美在线观看| 亚洲内射少妇av| 久久久久网色| 国内精品久久久久精免费| .国产精品久久| 欧美日韩综合久久久久久| 亚洲av第一区精品v没综合| 国产成人精品婷婷| 久久精品久久久久久噜噜老黄 | 激情 狠狠 欧美| 狂野欧美激情性xxxx在线观看| 国产精品日韩av在线免费观看| 欧美日韩乱码在线| 久久午夜亚洲精品久久| 国产精品伦人一区二区| 亚洲人成网站高清观看| 欧美+日韩+精品| 亚洲第一区二区三区不卡| 久久国产乱子免费精品| 床上黄色一级片| 亚洲婷婷狠狠爱综合网| 99riav亚洲国产免费| 1024手机看黄色片| 中文字幕免费在线视频6| 亚洲美女搞黄在线观看| 国产在视频线在精品| 久久精品国产亚洲av涩爱 | 欧美xxxx黑人xx丫x性爽| 中文字幕熟女人妻在线| 中出人妻视频一区二区| 在线免费观看不下载黄p国产| 婷婷色综合大香蕉| 国产成人一区二区在线| 久久久久久久久久久免费av| 欧美激情久久久久久爽电影| 亚洲经典国产精华液单| 日韩中字成人| 国产精品女同一区二区软件| 18禁在线无遮挡免费观看视频| 给我免费播放毛片高清在线观看| av卡一久久| 91久久精品国产一区二区成人| 日韩成人av中文字幕在线观看| 久久久久久国产a免费观看| 亚洲熟妇中文字幕五十中出| 爱豆传媒免费全集在线观看| 丰满的人妻完整版| 丰满人妻一区二区三区视频av| 久久久久久久久大av| 看非洲黑人一级黄片| 婷婷色综合大香蕉| 91麻豆精品激情在线观看国产| 日韩欧美三级三区| 听说在线观看完整版免费高清| 欧美日韩乱码在线| 99热这里只有是精品在线观看| 亚洲av成人精品一区久久| 国内精品一区二区在线观看| 不卡一级毛片| 久久人人精品亚洲av| 欧美bdsm另类| 97超视频在线观看视频| 久久人人爽人人爽人人片va| 亚洲精品粉嫩美女一区| 九九热线精品视视频播放| 乱人视频在线观看| 免费搜索国产男女视频| 一边摸一边抽搐一进一小说| 麻豆国产av国片精品| 又黄又爽又刺激的免费视频.| 特级一级黄色大片| 久久99蜜桃精品久久| 国产精品日韩av在线免费观看| 午夜免费激情av| 午夜精品一区二区三区免费看| 内地一区二区视频在线| 日韩一区二区视频免费看| 国产探花极品一区二区| 高清午夜精品一区二区三区 | 日本爱情动作片www.在线观看| 日韩av在线大香蕉| 高清毛片免费看| 99久久精品一区二区三区| 人人妻人人澡人人爽人人夜夜 | 国产精品一区二区三区四区久久| 亚洲国产精品成人综合色| 男人和女人高潮做爰伦理| 九九爱精品视频在线观看| 在线国产一区二区在线| 大型黄色视频在线免费观看| 亚洲欧美成人综合另类久久久 | 成人漫画全彩无遮挡| 久久午夜福利片| 国产成人aa在线观看| 蜜桃久久精品国产亚洲av| 国产精品野战在线观看| 国产精品.久久久| 日韩欧美精品v在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久中文看片网| 丝袜美腿在线中文| 中出人妻视频一区二区| 天堂av国产一区二区熟女人妻| 在线a可以看的网站| 免费观看a级毛片全部| 久久精品综合一区二区三区| 国产精品久久久久久久久免| 夜夜看夜夜爽夜夜摸| 一边亲一边摸免费视频| 欧美3d第一页| 成人毛片60女人毛片免费| or卡值多少钱| 国产成人a区在线观看| 在线观看av片永久免费下载| 精品一区二区三区视频在线| 欧美变态另类bdsm刘玥| 91午夜精品亚洲一区二区三区| 国产精品精品国产色婷婷| 校园春色视频在线观看| 日韩人妻高清精品专区| 欧美zozozo另类| 日韩大尺度精品在线看网址| 美女大奶头视频| 国产成人影院久久av| 国产欧美日韩精品一区二区| 午夜a级毛片| 麻豆成人午夜福利视频| 91精品一卡2卡3卡4卡| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久av不卡| 国产激情偷乱视频一区二区| 熟女电影av网| 国产成年人精品一区二区| 国模一区二区三区四区视频| 欧美色视频一区免费| 中出人妻视频一区二区| 赤兔流量卡办理| av在线蜜桃| 亚洲av男天堂| 亚洲国产精品成人综合色| 久久国内精品自在自线图片| 亚洲欧美中文字幕日韩二区| 国产一级毛片在线| 一本精品99久久精品77| 少妇高潮的动态图| 亚洲精品日韩av片在线观看| 免费大片18禁| 久久精品久久久久久久性| 国产av麻豆久久久久久久| 亚洲自拍偷在线| 久久午夜亚洲精品久久| 欧美日韩国产亚洲二区| 日韩三级伦理在线观看| 成熟少妇高潮喷水视频| 色哟哟哟哟哟哟| 丝袜美腿在线中文| 特级一级黄色大片| 国产一区二区三区在线臀色熟女| 尤物成人国产欧美一区二区三区| 免费看a级黄色片| 日本欧美国产在线视频| 久久久久免费精品人妻一区二区| www.色视频.com| 夜夜夜夜夜久久久久| 在线观看免费视频日本深夜| 国产视频首页在线观看| 高清毛片免费看| 亚洲真实伦在线观看| 亚洲婷婷狠狠爱综合网| 麻豆成人av视频| 亚洲熟妇中文字幕五十中出| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 青春草亚洲视频在线观看| 国产乱人偷精品视频| 中文字幕av成人在线电影| 两个人视频免费观看高清| 国产在线精品亚洲第一网站| 3wmmmm亚洲av在线观看| 国产黄色视频一区二区在线观看 | av国产免费在线观看| 午夜福利视频1000在线观看| 久久久久国产网址| 成人漫画全彩无遮挡| 老女人水多毛片| 国产真实伦视频高清在线观看| 嫩草影院新地址| 成人欧美大片| 久久精品影院6| 国产精品永久免费网站| 美女大奶头视频| 国产黄色视频一区二区在线观看 | 国产成人91sexporn| 欧美xxxx性猛交bbbb| 白带黄色成豆腐渣| 中国国产av一级| 国产探花在线观看一区二区| 久久久久久久久大av| 国产亚洲欧美98| 国产免费一级a男人的天堂| 狠狠狠狠99中文字幕| 欧美最黄视频在线播放免费| 亚洲av男天堂| 精品久久久久久久久久久久久| 国产三级中文精品| 亚洲人与动物交配视频| 精品熟女少妇av免费看| 日韩av在线大香蕉| 国产精品国产高清国产av| 精品久久久久久久久av| 中文字幕人妻熟人妻熟丝袜美| 2021天堂中文幕一二区在线观| 春色校园在线视频观看| 国产精品国产高清国产av| 久久精品综合一区二区三区| 国产成人精品一,二区 | 美女被艹到高潮喷水动态| 亚洲欧洲国产日韩| 91精品国产九色| 蜜臀久久99精品久久宅男| 国产亚洲精品久久久com| 久久久久国产网址| 亚洲av熟女| 伦精品一区二区三区| 亚洲18禁久久av| 国产乱人偷精品视频| 国产一级毛片七仙女欲春2| 一级毛片久久久久久久久女| 中文字幕久久专区| 国产精品一区二区三区四区久久| 又爽又黄a免费视频| 欧美日韩在线观看h| 国产片特级美女逼逼视频| 99久久九九国产精品国产免费| 淫秽高清视频在线观看| 99精品在免费线老司机午夜| 爱豆传媒免费全集在线观看| 国产av一区在线观看免费| 国产三级中文精品| 亚洲婷婷狠狠爱综合网| 热99re8久久精品国产| 99在线人妻在线中文字幕| 久久久色成人| 麻豆成人av视频| 欧美精品国产亚洲| 女的被弄到高潮叫床怎么办| 亚洲人成网站在线播放欧美日韩| 精品久久国产蜜桃| 直男gayav资源| 全区人妻精品视频| 综合色av麻豆| 日韩视频在线欧美| av国产免费在线观看| 午夜亚洲福利在线播放| 九九爱精品视频在线观看| 精品少妇黑人巨大在线播放 | 成人毛片60女人毛片免费| 成人漫画全彩无遮挡| 精品99又大又爽又粗少妇毛片| 少妇熟女aⅴ在线视频| 爱豆传媒免费全集在线观看| 亚洲精品影视一区二区三区av| 级片在线观看| 日韩 亚洲 欧美在线| 国产色婷婷99| 国产真实伦视频高清在线观看| 在线观看美女被高潮喷水网站| 高清在线视频一区二区三区 | 欧美高清成人免费视频www| 成人三级黄色视频| 春色校园在线视频观看| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 亚洲精品国产av成人精品| 尾随美女入室| 亚洲成人中文字幕在线播放| 国产精品久久视频播放| 日韩视频在线欧美| 亚洲精品成人久久久久久| 能在线免费看毛片的网站| 波多野结衣巨乳人妻| 精品人妻视频免费看| 六月丁香七月| 久久久欧美国产精品| 久久精品国产自在天天线| 国产老妇伦熟女老妇高清| 啦啦啦韩国在线观看视频| 精品人妻一区二区三区麻豆| 久久精品夜色国产| 亚洲无线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 狠狠狠狠99中文字幕| 亚洲天堂国产精品一区在线| 男的添女的下面高潮视频| 日本av手机在线免费观看| 国产在视频线在精品| 日韩一本色道免费dvd| 我要搜黄色片| 日韩大尺度精品在线看网址| 国产av在哪里看| 一本一本综合久久| 亚洲最大成人中文| 国产av不卡久久| 亚洲av不卡在线观看| 日本黄色视频三级网站网址| 亚洲最大成人av| 国内久久婷婷六月综合欲色啪| 男的添女的下面高潮视频| 国模一区二区三区四区视频| 三级男女做爰猛烈吃奶摸视频| 亚洲aⅴ乱码一区二区在线播放| 欧美激情国产日韩精品一区| 一个人免费在线观看电影| 91精品一卡2卡3卡4卡| 一级毛片电影观看 | 最近最新中文字幕大全电影3| 精品国产三级普通话版| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 国产91av在线免费观看| 草草在线视频免费看| 精品久久久久久成人av| 久久这里只有精品中国| 在线播放国产精品三级| 我的女老师完整版在线观看| 国产高清三级在线| 97在线视频观看| 日韩成人伦理影院| 日日干狠狠操夜夜爽| 天堂网av新在线| 国产av不卡久久| 看黄色毛片网站| 中文字幕熟女人妻在线| 国产伦一二天堂av在线观看| 欧美性猛交╳xxx乱大交人| 亚洲成av人片在线播放无| 国产视频首页在线观看| 国产伦理片在线播放av一区 | 亚洲三级黄色毛片| 久久久久久国产a免费观看| 欧美最黄视频在线播放免费| 国产精品日韩av在线免费观看| 一级毛片aaaaaa免费看小| 亚洲一级一片aⅴ在线观看| av黄色大香蕉| 国产精品三级大全| 亚洲在久久综合| 久久韩国三级中文字幕| 国产老妇女一区| 老司机影院成人| 日本免费一区二区三区高清不卡| 在线免费观看不下载黄p国产| 亚洲精品亚洲一区二区| 国产亚洲欧美98| 少妇丰满av| 高清在线视频一区二区三区 | 狂野欧美激情性xxxx在线观看| 人人妻人人看人人澡| 久久久久久久午夜电影| 看黄色毛片网站| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 日日摸夜夜添夜夜爱| 少妇高潮的动态图| 99久国产av精品国产电影| 亚洲中文字幕一区二区三区有码在线看| 国产人妻一区二区三区在| 成人毛片a级毛片在线播放| 99视频精品全部免费 在线| 国产黄片视频在线免费观看| 看十八女毛片水多多多| 一本久久精品| 亚洲精品乱码久久久久久按摩| 51国产日韩欧美| 国产真实伦视频高清在线观看| 国产v大片淫在线免费观看| 婷婷色综合大香蕉| av在线老鸭窝| or卡值多少钱| 国产白丝娇喘喷水9色精品| 天天躁日日操中文字幕| 日本成人三级电影网站| 99久久精品国产国产毛片| 亚洲一区二区三区色噜噜| 久久国内精品自在自线图片| 亚洲国产精品合色在线| 黄色视频,在线免费观看| 色综合亚洲欧美另类图片| 精品一区二区三区视频在线| 夫妻性生交免费视频一级片| 深爱激情五月婷婷| 国产精品电影一区二区三区| 美女大奶头视频| 久久久色成人| 美女脱内裤让男人舔精品视频 | 日日摸夜夜添夜夜爱| 久久久a久久爽久久v久久| 国产亚洲5aaaaa淫片| 国产老妇伦熟女老妇高清| 不卡视频在线观看欧美| 精品久久久久久久久av| 久久久久久久久久久丰满| 中文欧美无线码| 国产精品一二三区在线看| 两个人的视频大全免费| 成年女人永久免费观看视频| 日本三级黄在线观看| 久久99蜜桃精品久久| 简卡轻食公司| 日本免费一区二区三区高清不卡| 日韩高清综合在线| 美女脱内裤让男人舔精品视频 | 伦精品一区二区三区| 两个人视频免费观看高清| av天堂中文字幕网| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 老司机福利观看| 夫妻性生交免费视频一级片| 综合色丁香网| 女人被狂操c到高潮| 久久久久久久久久久丰满| 国产女主播在线喷水免费视频网站 | 亚洲国产精品成人久久小说 | 97热精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 激情 狠狠 欧美| 91精品一卡2卡3卡4卡| 床上黄色一级片| 国产精品电影一区二区三区| 深夜精品福利| 啦啦啦观看免费观看视频高清| 我的女老师完整版在线观看| 免费黄网站久久成人精品| 自拍偷自拍亚洲精品老妇| 最近2019中文字幕mv第一页| 国产一级毛片七仙女欲春2| 能在线免费看毛片的网站| 国产黄片美女视频| 亚洲精品影视一区二区三区av| 成人高潮视频无遮挡免费网站| 成人漫画全彩无遮挡| 日韩欧美在线乱码| 听说在线观看完整版免费高清| 久久这里有精品视频免费| 国产高清视频在线观看网站| 简卡轻食公司| 高清在线视频一区二区三区 | 女人被狂操c到高潮| 国产爱豆传媒在线观看| 国内精品宾馆在线| 精品一区二区三区人妻视频| 亚洲欧美清纯卡通| 身体一侧抽搐| 人妻少妇偷人精品九色| 两个人的视频大全免费| 久久国内精品自在自线图片| 美女脱内裤让男人舔精品视频 | 国产大屁股一区二区在线视频| 精品国产三级普通话版| 国产亚洲欧美98| 婷婷色综合大香蕉| 亚洲久久久久久中文字幕| 最新中文字幕久久久久| 老司机福利观看| 久久这里只有精品中国| 99热这里只有精品一区| 一区福利在线观看| 又爽又黄a免费视频| 久久亚洲精品不卡| 久久精品国产鲁丝片午夜精品| 精品国产三级普通话版| 色尼玛亚洲综合影院| 午夜a级毛片| 老女人水多毛片| 又粗又硬又长又爽又黄的视频 | 久久精品国产自在天天线| 久久久久久久亚洲中文字幕| 一级二级三级毛片免费看| 国产高清不卡午夜福利| 亚洲成人久久爱视频| 少妇被粗大猛烈的视频| 久久久久久久亚洲中文字幕| 18禁在线播放成人免费| 九九在线视频观看精品| 午夜激情欧美在线| 免费电影在线观看免费观看| 国产精品野战在线观看| 日韩av在线大香蕉| 此物有八面人人有两片| 91精品国产九色| 又粗又硬又长又爽又黄的视频 | 99在线人妻在线中文字幕| a级一级毛片免费在线观看| 91在线精品国自产拍蜜月| 人体艺术视频欧美日本| 狂野欧美白嫩少妇大欣赏| 成人国产麻豆网| 欧美高清成人免费视频www| 亚洲国产精品成人综合色| 成人毛片a级毛片在线播放| 亚洲欧美精品综合久久99| 国产精品人妻久久久影院| 亚洲成av人片在线播放无| 久久午夜福利片| 国产午夜精品久久久久久一区二区三区| 国产成人a区在线观看| 校园人妻丝袜中文字幕| 婷婷精品国产亚洲av| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月| 国产欧美日韩精品一区二区| 狠狠狠狠99中文字幕| 日本-黄色视频高清免费观看| 欧美丝袜亚洲另类| 欧美bdsm另类| 成人一区二区视频在线观看| 69av精品久久久久久| ponron亚洲| 亚洲天堂国产精品一区在线| 99九九线精品视频在线观看视频| 夫妻性生交免费视频一级片| 最近2019中文字幕mv第一页| 男女做爰动态图高潮gif福利片| 亚洲丝袜综合中文字幕| 免费搜索国产男女视频| 又粗又爽又猛毛片免费看| 高清毛片免费看| 一级毛片久久久久久久久女| 亚洲第一区二区三区不卡| 嘟嘟电影网在线观看| 亚洲国产欧洲综合997久久,| 人妻制服诱惑在线中文字幕| 亚洲国产精品成人久久小说 | 大又大粗又爽又黄少妇毛片口| а√天堂www在线а√下载| 身体一侧抽搐|