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

    基于有限元-離散元結(jié)合方法的Steven 實(shí)驗(yàn)三維數(shù)值模擬

    2020-10-31 07:57:36黃彬彬劉倉(cāng)理
    含能材料 2020年10期
    關(guān)鍵詞:平頭彈頭炸藥

    黃彬彬,傅 華,喻 寅,劉倉(cāng)理

    (中國(guó)工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點(diǎn)實(shí)驗(yàn)室,四川 綿陽 621900)

    1 引言

    炸藥在跌落、翻倒等低速撞擊事故中易經(jīng)受低幅值(數(shù)十到數(shù)百兆帕)、長(zhǎng)脈沖(數(shù)百微秒)的非沖擊載荷,可能造成炸藥發(fā)生意外點(diǎn)火,引發(fā)安全性事故。但炸藥非沖擊點(diǎn)火的實(shí)驗(yàn)精細(xì)觀測(cè)存在多方面的困難與限制,不易獲得點(diǎn)火演化的細(xì)節(jié)與機(jī)制。因此,開展非沖擊點(diǎn)火的數(shù)值模擬研究,對(duì)深入揭示炸藥在低速撞擊下的響應(yīng)規(guī)律具有重要意義。

    Chidester 等[1-3]設(shè) 計(jì) 了Steven 實(shí) 驗(yàn) 用 于 研 究 帶殼體約束的炸藥經(jīng)受低速撞擊下的非沖擊點(diǎn)火問題,建立了撞擊點(diǎn)火模型并采用實(shí)驗(yàn)獲得的臨界撞擊速度對(duì)模型參數(shù)進(jìn)行標(biāo)定。此后,眾多國(guó)內(nèi)外學(xué)者開展了相應(yīng)的數(shù)值模擬工作,分析了不同形狀彈頭加載下的點(diǎn)火速度閾值和點(diǎn)火位置等[4-7]。目前相關(guān)的數(shù)值模擬工作主要采用了有限元方法與離散元方法。其中有限元方法(Finite Element Method,F(xiàn)EM)具有高精度、高效率等優(yōu)點(diǎn),是目前最成熟且應(yīng)用最廣泛的力學(xué)響應(yīng)計(jì)算方法,但是依賴網(wǎng)格的計(jì)算方法普遍存在不適合模擬大變形和損傷斷裂現(xiàn)象的缺點(diǎn),雖然基于微裂紋界面的摩擦生熱的點(diǎn)火模型(Visco-SCRAM 模型)在炸藥低速撞擊點(diǎn)火數(shù)值模擬中也有較多應(yīng)用[8-11],但在處理撞擊加載下炸藥局域點(diǎn)火前的壓剪斷裂與斷裂之后的裂紋摩擦等非連續(xù)性問題時(shí)還是面臨較大的困難。離散元方法(Discrete Element Method,DEM)的基本思想是把不連續(xù)體或連續(xù)體離散為具有一定物理意義的獨(dú)立“微元”或“粒子”,粒子間通過虛彈簧保持連接,粒子間的相互作用力通過彈簧的拉伸和壓縮來實(shí)現(xiàn),因此在處理斷裂和破碎問題時(shí)具有天然優(yōu)勢(shì),目前也被廣泛應(yīng)用于炸藥力學(xué)響應(yīng)機(jī)理的模擬研究[12-18],但是由于離散元不擅長(zhǎng)描述材料的塑性變形,同時(shí)出于計(jì)算尺度和效率的考量,在應(yīng)用于炸藥材料模擬時(shí),通常模擬中未考慮炸藥外部的金屬殼體對(duì)其非沖擊點(diǎn)火行為的影響。

    本研究中建立了有限元-離散元結(jié)合方法,利用有限元模擬裝藥結(jié)構(gòu)響應(yīng),利用離散元描述炸藥局部斷裂損傷,針對(duì)脆性炸藥PBX-2 的Steven 實(shí)驗(yàn)開展了三維數(shù)值模擬,獲得了裝藥受載應(yīng)力歷史、裂紋分布、溫升曲線及溫度分布。以炸藥裂紋摩擦以及自反應(yīng)放熱為主導(dǎo)點(diǎn)火機(jī)制,分析了由炸藥材料非均質(zhì)特性引起的的概率點(diǎn)火現(xiàn)象,并討論了鈍頭彈與平頭彈撞擊下點(diǎn)火位置的分布特性。

    2 計(jì)算方法及模型

    2.1 FEM-DEM 結(jié)合方法

    FEM-DEM 結(jié)合方法的核心問題在于有限元單元與離散元粒子間接觸作用力的計(jì)算。由于離散元單元本身具有一定尺寸,區(qū)別于有限元節(jié)點(diǎn)僅僅是一個(gè)點(diǎn),因此離散元粒子-有限元面單元的接觸可以根據(jù)其接觸位置的不同分為三類,即:粒子-表面(Paticles-to-Facets,簡(jiǎn)稱PTF)、粒子-棱邊(Paticles-to-Edges,簡(jiǎn)稱PTE)、粒子-頂點(diǎn)(Paticles-to-Vertices,簡(jiǎn)稱PTV),如圖1所示。

    三種接觸方式下的離散元單元所受到的作用力為;

    式中,κ為自定義接觸系數(shù);δn為侵徹距離,m;n為單位法向量;E*為等效楊氏模量,MPa;r*為等效半徑,m;k法向剛度,N·m-1。

    圖1 離散元單元與有限元面單元接觸模型Fig.1 Contact model between a discrete-element particle and a finite-element facet unit

    2.2 DEM 斷裂模型

    PBX-2 各組分的質(zhì)量百分比為HMX/TATB/黏結(jié)劑+鈍感劑=87/7/6,屬于典型的脆性炸藥,因此采用彈脆性本構(gòu)對(duì)其進(jìn)行描述。假設(shè)材料為各向同性,建立的離散元模型如圖2 所示,圖2a 為球形離散元結(jié)構(gòu)示意圖,圖2b 為相鄰單元間作用力示意圖。

    圖2 離散元模型示意圖Fig.2 Schematic diagram of discrete-element model

    當(dāng)兩個(gè)處于連接狀態(tài)的離散元單元間法向作用力fn與切向作用力fs滿足關(guān)系式(5)中關(guān)系時(shí),則判定其連接斷裂,兩個(gè)單元由連接狀態(tài)變?yōu)榻佑|狀態(tài):

    式中,fn0和fs0分別為定義的最大法向與切向作用力,N。

    為驗(yàn)證該計(jì)算方法以及所選取模型的有效性,對(duì)脆性炸藥PBX-2 的分離式霍普金森桿實(shí)驗(yàn)(SHPB)進(jìn)行了數(shù)值模擬。模型包含炸藥試樣、子彈、整形器、入射桿和透射桿,炸藥試樣尺寸為φ12 mm×6 mm,子彈尺寸為φ20 mm×300 mm,入射桿尺寸為φ20 mm×2000 mm,透射桿尺寸為φ20 mm×1000 mm。炸藥試樣采用離散元?jiǎng)澐志W(wǎng)格,半徑為125 μm,其余部分采用有限元?jiǎng)澐志W(wǎng)格,網(wǎng)格尺寸為2 mm,如圖3 所示。分別于A、B、C三點(diǎn)處記錄入射波、反射波和透射波信息。

    “民營(yíng)企業(yè)反映突出的痛點(diǎn)、難點(diǎn)、堵點(diǎn),給出了6方面35條政策,讓我們真正從政策中增強(qiáng)獲得感。”記者在采訪謝清森時(shí),他用這句話表達(dá)自己內(nèi)心的最大感受。

    圖3 SHPB 實(shí)驗(yàn)計(jì)算模型Fig.3 Calculated model of SHPB test

    計(jì)算模型中采用彈脆性本構(gòu),彈性段彈性模量E=5.5 GPa,平均斷裂應(yīng)變?chǔ)?0.01。模擬結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比如圖4 所示。從圖4a中可以看出,計(jì)算得到的三波(入射波、透射波、反射波)波形與實(shí)驗(yàn)結(jié)果符合較好,表明該計(jì)算方法能較好地處理FEM-DEM 接觸界面間的物理量傳遞。圖4b中應(yīng)力-應(yīng)變本構(gòu)關(guān)系模擬結(jié)果與實(shí)驗(yàn)結(jié)果大致吻合,呈現(xiàn)明顯的彈脆性特征,表明所采用的斷裂模型適用于模擬脆性炸藥材料的力學(xué)響應(yīng)。

    圖4 SHPB 實(shí)驗(yàn)?zāi)M結(jié)果Fig.4 Simulation results of SHPB experiment

    2.3 溫升點(diǎn)火模型

    考慮到相比炸藥塑性功產(chǎn)生的熱量,顆粒間的摩擦溫升更容易產(chǎn)生較高的局部溫升導(dǎo)致點(diǎn)火,因此主要計(jì)算炸藥斷裂后的摩擦生熱以及炸藥自反應(yīng)放熱作為熱源。

    兩個(gè)接觸的離散元粒子i與j之間的摩擦力可表示為:

    式中,eij為耗散能量,J;Δt為時(shí)間步長(zhǎng),s。

    由于加載過程為瞬態(tài)過程,一個(gè)時(shí)間步長(zhǎng)內(nèi)單元間相互作用可視為絕熱過程,則粒子i的摩擦溫升為:

    采用Arrhenius 方程[17]描述炸藥自反應(yīng)放熱

    式中,S為源項(xiàng),J?s-1;ρ為炸藥密度,g?cm-3;Q為反應(yīng)熱,J?kg-1;Z為指前因子,s-1;E為活化能,J?mol-1;R為普適氣體常數(shù),8.314 J?mol-1?K-1;T為溫度,K;α 為已反應(yīng)炸藥分?jǐn)?shù)。初始溫度設(shè)為298 K,炸藥反應(yīng)動(dòng)力學(xué)參數(shù)見表1[18]。

    表1 PBX-2 炸藥反應(yīng)動(dòng)力學(xué)參數(shù)[18]Table 1 Kinetic parameters of explosive reaction for PBX-2[18]

    當(dāng)某個(gè)炸藥單元的溫度超過某個(gè)臨界溫度,根據(jù)Arrhenius 方程,將大量放熱導(dǎo)致溫度近乎直線上升超過3000 K,這種情況下我們可以判定炸藥發(fā)生了點(diǎn)火。因此,可以將溫度出現(xiàn)快速上升的拐點(diǎn)作為點(diǎn)火臨界點(diǎn),其對(duì)應(yīng)的溫度為點(diǎn)火臨界溫度。

    2.4 Steven 實(shí)驗(yàn)計(jì)算模型

    模型包含炸藥試樣、蓋板、約束環(huán)、樣品盒和彈丸。炸藥試樣尺寸為Φ98 mm×13 mm,壓板尺寸為Φ120 mm×3.5 mm,樣品盒底部高度19 mm,彈丸主體部分半徑為31 mm。炸藥試樣采用離散元?jiǎng)澐志W(wǎng)格(半徑400 μm),其余部分采用有限元?jiǎng)澐志W(wǎng)格(網(wǎng)格尺寸4 mm),建立的FEM-DEM 模型如圖5 所示。

    圖5 不同形狀彈頭Steven 實(shí)驗(yàn)計(jì)算模型Fig.5 Calculated model of Steven impact test for different projectile heads

    模型中所用樣品為PBX-2 炸藥,彈丸、蓋板和樣品盒的材料為Q235 鋼,約束環(huán)材料為聚四氟乙烯,模型中各材料性能參數(shù)見表2[18]。

    實(shí)際加載中,炸藥材料的孔洞坍縮、炸藥-粘結(jié)劑脫粘等行為會(huì)進(jìn)一步導(dǎo)致響應(yīng)局域化,因此模型中需要考慮炸藥材料的非均勻性。為模擬炸藥材料的非均勻性,選取隨機(jī)分布的占比為1%的離散元單元,設(shè)置其最大法向與切向作用力為其余單元的一半。

    表2 模型中各材料性能參數(shù)Table 2 Performance parameters for various materials in the model

    3 結(jié)果與討論

    3.1 力學(xué)響應(yīng)與損傷裂紋

    采用有限元/離散元結(jié)合方法得到鈍頭彈在撞擊速度為45 m·s-1時(shí)炸藥樣品底部中心處壓力的模擬結(jié)果,將其與文獻(xiàn)[15]實(shí)測(cè)曲線進(jìn)行對(duì)比,如圖6 所示。由圖6 可以看出,試樣底部中心處的壓力峰值約為0.10 GPa。彈丸撞擊后壓力迅速升高,在約170 μs 處達(dá)到峰值,之后隨彈丸反彈壓力持續(xù)下降,進(jìn)入卸載階段,卸載階段持續(xù)380 μs 左右。數(shù)值模擬的壓力峰值及壓力脈沖寬度都與實(shí)驗(yàn)測(cè)試曲線符合較好,表明該方法能較準(zhǔn)確地模擬Steven 實(shí)驗(yàn)中炸藥的受力過程。

    圖6 炸藥底部中心實(shí)測(cè)和計(jì)算壓力曲線(鈍頭彈彈速45 m·s-1)Fig.6 Experimental and simulated pressures at the bottom center of explosives(the velocity of blunt projectile is 45 m·s-1)

    圖7 和圖8 為不同彈丸撞擊下炸藥試樣內(nèi)部產(chǎn)生的力學(xué)響應(yīng)及損傷裂紋分布情況。在鈍頭彈實(shí)驗(yàn)中,裂紋由彈頭軸線處沿徑向發(fā)散,形成多段貫穿裂紋(圖7a 和圖7b);在平頭彈實(shí)驗(yàn)中,裂紋沿彈頭圓周附近呈環(huán)狀分布(圖7c 和圖7d)。對(duì)比圖8 中的應(yīng)力分布情況可知,裝藥試樣受到鈍頭彈撞擊后,撞擊中心的炸藥受到彈丸擠壓(圖7a),該區(qū)域的炸藥向遠(yuǎn)離撞擊中心的方向運(yùn)動(dòng),應(yīng)力最大值始終集中于軸線附近(圖7b);裝藥試樣受到平頭彈撞擊后,在彈頭邊緣與蓋板接觸處發(fā)生強(qiáng)烈的剪切作用(圖7c),從而形成繞圓周分布的損傷裂紋,而彈頭下方的炸藥在軸向被壓實(shí)(圖7d)的同時(shí)也會(huì)產(chǎn)生較分散的損傷裂紋。

    3.2 點(diǎn)火閾值與點(diǎn)火位置

    圖7 不同形狀彈頭撞擊下PBX-2 炸藥的裂紋分布Fig7 Crack distribution of PBX-2 under the impact of different projectile heads

    圖8 不同形狀彈頭撞擊下PBX-2 炸藥的應(yīng)力分布Fig.8 Stress distribution of PBX-2 under the impact of different projectile heads

    圖9 為鈍頭彈與平頭彈以不同速度撞擊時(shí)炸藥內(nèi)部最高溫度單元的溫度-時(shí)間曲線。由圖9 可以看出,當(dāng)鈍頭彈與平頭彈速度低于40 m·s-1時(shí),炸藥溫度上升比較緩慢,當(dāng)達(dá)到一定溫度后便不再升高,未發(fā)生點(diǎn)火。當(dāng)鈍頭彈與平頭彈的速度分別大于45 m·s-1和50 m·s-1時(shí),炸藥溫度以更快的速率上升,達(dá)到點(diǎn)火溫度后溫度急劇上升發(fā)生點(diǎn)火,且點(diǎn)火時(shí)間隨彈速的增加而減小。分析可知,撞擊速度較低時(shí),炸藥裂紋出現(xiàn)的時(shí)間較晚,摩擦溫升比較緩慢;隨著撞擊速度增加,裂紋間摩擦加劇,自反應(yīng)放熱加快,形成局部熱點(diǎn)導(dǎo)致點(diǎn)火。

    圖9 不同速度彈丸撞擊下的溫升曲線Fig.9 Temperature rise curve under the impact of different projectile velocities

    圖10 為炸藥達(dá)到點(diǎn)火臨界點(diǎn)時(shí)的溫度分布。鈍頭彈加載下,達(dá)到點(diǎn)火溫度的炸藥單元集中于彈頭軸線下方附近;平頭彈加載下,達(dá)到點(diǎn)火溫度的炸藥單元主要位于彈頭圓周下方,呈環(huán)狀分布。結(jié)合圖7 中的裂紋分布與圖8 中的應(yīng)力分布可知,鈍頭彈加載下,軸線下方炸藥單元所承受的壓力更大,發(fā)生斷裂后單元間的摩擦作用更強(qiáng),導(dǎo)致更高的摩擦溫升,加速自反應(yīng)放熱,從而更快到達(dá)點(diǎn)火溫度;平頭彈加載下,彈頭邊緣與蓋板接觸處發(fā)生強(qiáng)烈的剪切作用,發(fā)生斷裂后此處的炸藥單元間相對(duì)位移速度更大,導(dǎo)致更高的摩擦溫升,加速自反應(yīng)放熱,從而更快到達(dá)點(diǎn)火溫度。

    圖10 不同形狀彈頭撞擊下炸藥的溫度分布圖Fig.10 Temperature distribution of explosives under the impact of different projectile heads

    表3 鈍頭彈不同撞擊速度下的點(diǎn)火情況Table 3 Ignition situation of blunt projectile under the impact of different velocities

    重復(fù)隨機(jī)生成炸藥非均勻模型,對(duì)每速度撞擊下的點(diǎn)火情況各進(jìn)行25 次模擬,模擬結(jié)果見表3。借鑒落錘實(shí)驗(yàn)炸藥感度判定標(biāo)準(zhǔn)[16](一定落高下重復(fù)實(shí)驗(yàn)點(diǎn)火概率超過50%),可得數(shù)值模擬中使炸藥點(diǎn)火的速度閾值分別約為43 m·s-1(鈍頭彈)和47 m·s-1(平頭彈),與文獻(xiàn)中結(jié)果(鈍頭彈43 m·s-1,平頭彈47 m·s-1)對(duì)比大致吻合,表明本文中建立的計(jì)算模型可以較好地模擬Steven 實(shí)驗(yàn)。

    圖11 為表3 中點(diǎn)火概率100%時(shí)的25 次模擬中點(diǎn)火位置的平面分布統(tǒng)計(jì):鈍頭彈撞擊后,點(diǎn)火位置集中分布于距離軸心4~11 mm 處,即點(diǎn)火點(diǎn)圍繞撞擊軸產(chǎn)生;對(duì)于平頭彈,點(diǎn)火位置集中分布于距離軸心27~32 mm 處,點(diǎn)火點(diǎn)呈環(huán)狀分布,與彈頭直徑大致相當(dāng)。圖11 中的點(diǎn)火位置分布與文獻(xiàn)[5]中實(shí)驗(yàn)觀測(cè)到的結(jié)果大致相符。

    圖11 不同形狀彈頭撞擊下炸藥點(diǎn)火位置平面分布統(tǒng)計(jì)Fig.11 Plane distribution statistics of explosive ignition positions under the impact of different projectile heads

    4 結(jié)論

    本研究采用FEM-DEM 結(jié)合方法對(duì)Steven 實(shí)驗(yàn)中脆性炸藥PBX-2 的低速撞擊破壞和點(diǎn)火行為進(jìn)行了研究。

    (1)建立了基于FEM-DEM 結(jié)合方法的三維非沖擊點(diǎn)火計(jì)算模型并編寫了Fortran 計(jì)算程序,通過對(duì)PBX-2 炸藥的SHPB 實(shí)驗(yàn)?zāi)M,驗(yàn)證了程序處理FEM-DEM 界面物理量傳遞的可靠性,以及所采用的DEM 斷裂模型模擬脆性炸藥斷裂破壞過程的合理性。

    (2)考慮了炸藥斷裂破碎、裂紋摩擦升溫與炸藥熱分解,模擬了Steven 實(shí)驗(yàn)中在鈍頭彈與平頭彈加載下PBX-2 炸藥的受載應(yīng)力歷史、裂紋分布、溫升曲線等結(jié)果。鈍頭彈實(shí)驗(yàn)中,裂紋由彈頭軸線處沿徑向發(fā)散,形成多段貫穿裂紋,中心底部形成明顯凹陷;在平頭彈實(shí)驗(yàn)中,裂紋沿彈頭圓周附近呈環(huán)狀分布。

    (3)基于考慮了非均質(zhì)特性的足夠數(shù)量PBX 炸藥樣本的模擬,給出了考慮隨機(jī)分散特性的Steven 實(shí)驗(yàn)中的炸藥點(diǎn)火閾值速度與點(diǎn)火位置的統(tǒng)計(jì)結(jié)果。模擬獲得的鈍頭彈與平頭彈撞擊下引起點(diǎn)火的閾值速度分別為43 m·s-1和47 m·s-1。在鈍頭彈加載下點(diǎn)火位置集中分布在撞擊中軸線附近;在平頭彈加載下形成點(diǎn)火環(huán)狀帶,且環(huán)狀直徑與彈頭直徑相當(dāng)。本研究從炸藥裂紋摩擦以及自反應(yīng)放熱的角度解釋脆性炸藥的非沖擊點(diǎn)火現(xiàn)象,模擬結(jié)果與實(shí)驗(yàn)報(bào)道定性相符。

    猜你喜歡
    平頭彈頭炸藥
    “炸藥”驚魂
    熱身管下小口徑槍彈彈頭殼材料對(duì)其膛內(nèi)運(yùn)動(dòng)的影響規(guī)律
    議論火炸藥數(shù)字化制造
    “平頭哥”的悲劇之路
    雜文月刊(2020年1期)2020-02-20 14:28:46
    動(dòng)物奇異行為觀察報(bào)道(三)
    導(dǎo)彈突防概率解析與仿真計(jì)算結(jié)果差異分析
    一字之差
    遼河(2014年5期)2014-05-23 10:45:33
    Al粉對(duì)炸藥爆炸加速能力的影響
    DNAN基熔鑄復(fù)合炸藥的爆轟性能
    碳納米管填充PTFE的制備及平頭型小沖桿測(cè)試研究
    国产成年人精品一区二区| e午夜精品久久久久久久| 男插女下体视频免费在线播放| 最新在线观看一区二区三区| 欧美高清成人免费视频www| а√天堂www在线а√下载| 一本一本综合久久| 少妇人妻一区二区三区视频| 色吧在线观看| 青草久久国产| 亚洲 欧美 日韩 在线 免费| 身体一侧抽搐| 欧美日韩精品网址| 国产精品av久久久久免费| 亚洲熟女毛片儿| 亚洲专区国产一区二区| 国产乱人视频| 他把我摸到了高潮在线观看| 人妻丰满熟妇av一区二区三区| 国内精品一区二区在线观看| 少妇丰满av| 我要搜黄色片| 久久热在线av| 免费在线观看日本一区| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久av网站| 老鸭窝网址在线观看| e午夜精品久久久久久久| 亚洲国产看品久久| 观看美女的网站| 国产淫片久久久久久久久 | 精品久久久久久,| 淫秽高清视频在线观看| 亚洲av五月六月丁香网| 巨乳人妻的诱惑在线观看| 久久中文字幕人妻熟女| 久久久久国产一级毛片高清牌| 亚洲国产日韩欧美精品在线观看 | 国产视频内射| 亚洲自偷自拍图片 自拍| 精品无人区乱码1区二区| 久久国产精品人妻蜜桃| 1000部很黄的大片| 精品一区二区三区av网在线观看| 曰老女人黄片| 他把我摸到了高潮在线观看| 国产一区二区三区视频了| 巨乳人妻的诱惑在线观看| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 国产高清videossex| 黄色日韩在线| 他把我摸到了高潮在线观看| 后天国语完整版免费观看| 一本精品99久久精品77| 亚洲中文字幕一区二区三区有码在线看 | 欧美一级a爱片免费观看看| 两性午夜刺激爽爽歪歪视频在线观看| a级毛片a级免费在线| 国产熟女xx| 丝袜人妻中文字幕| 男人和女人高潮做爰伦理| 国产午夜福利久久久久久| 免费在线观看日本一区| 99热精品在线国产| 日本一本二区三区精品| 老司机福利观看| 日本在线视频免费播放| 在线播放国产精品三级| 特级一级黄色大片| 国产三级中文精品| 久久久久久久久中文| 制服人妻中文乱码| 99热精品在线国产| 香蕉久久夜色| 不卡一级毛片| 人妻夜夜爽99麻豆av| 99国产精品一区二区三区| 日本熟妇午夜| 亚洲欧美日韩高清在线视频| 亚洲国产色片| 日本黄色视频三级网站网址| 午夜久久久久精精品| 在线观看日韩欧美| 1024手机看黄色片| 亚洲国产欧美一区二区综合| 欧美成人免费av一区二区三区| 美女黄网站色视频| 久久精品91蜜桃| 日本在线视频免费播放| 成在线人永久免费视频| 神马国产精品三级电影在线观看| 日韩中文字幕欧美一区二区| 亚洲精品美女久久久久99蜜臀| 国产高清视频在线播放一区| 精品不卡国产一区二区三区| 国语自产精品视频在线第100页| 精品一区二区三区四区五区乱码| 99在线人妻在线中文字幕| 日本与韩国留学比较| 精品电影一区二区在线| 国产成+人综合+亚洲专区| 国产精品,欧美在线| 一个人免费在线观看的高清视频| 两个人视频免费观看高清| 校园春色视频在线观看| 国产一级毛片七仙女欲春2| 国产av不卡久久| 午夜精品久久久久久毛片777| www.www免费av| 级片在线观看| а√天堂www在线а√下载| 日韩人妻高清精品专区| 舔av片在线| 色视频www国产| 久久久久国产精品人妻aⅴ院| 久久亚洲精品不卡| 亚洲人成电影免费在线| 免费看光身美女| 午夜a级毛片| 手机成人av网站| 99久国产av精品| 精品久久久久久久久久久久久| 男女视频在线观看网站免费| 人妻久久中文字幕网| 三级男女做爰猛烈吃奶摸视频| 国产精品免费一区二区三区在线| 国产淫片久久久久久久久 | 日本精品一区二区三区蜜桃| 琪琪午夜伦伦电影理论片6080| 国产97色在线日韩免费| 99久国产av精品| 最近最新中文字幕大全电影3| 欧美高清成人免费视频www| 精品国产乱子伦一区二区三区| 亚洲精品美女久久久久99蜜臀| 99热只有精品国产| 母亲3免费完整高清在线观看| 国产亚洲欧美98| 国产激情偷乱视频一区二区| 欧美日韩综合久久久久久 | 麻豆一二三区av精品| 久久久精品欧美日韩精品| 变态另类丝袜制服| 国产精品日韩av在线免费观看| 亚洲人与动物交配视频| 久久久久久久久久黄片| 欧美性猛交黑人性爽| 成人午夜高清在线视频| 黄色视频,在线免费观看| 色综合站精品国产| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 搞女人的毛片| 色哟哟哟哟哟哟| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人av| 精品午夜福利视频在线观看一区| 日韩中文字幕欧美一区二区| 亚洲av日韩精品久久久久久密| 精品久久蜜臀av无| www.www免费av| 精品乱码久久久久久99久播| 无人区码免费观看不卡| 日韩欧美免费精品| 国产精品 欧美亚洲| 亚洲 欧美 日韩 在线 免费| 观看美女的网站| 级片在线观看| 国产三级在线视频| 搡老岳熟女国产| 国产蜜桃级精品一区二区三区| 成人欧美大片| 9191精品国产免费久久| 欧美绝顶高潮抽搐喷水| 狂野欧美激情性xxxx| 免费观看人在逋| 久久亚洲精品不卡| 成人亚洲精品av一区二区| 在线播放国产精品三级| 色综合欧美亚洲国产小说| 91av网一区二区| 欧美午夜高清在线| 国产成人精品无人区| 日本撒尿小便嘘嘘汇集6| 国产激情偷乱视频一区二区| 法律面前人人平等表现在哪些方面| 日韩av在线大香蕉| 免费观看精品视频网站| 国产极品精品免费视频能看的| 深夜精品福利| 国产精品久久视频播放| 国产综合懂色| 成年女人看的毛片在线观看| 国产高清激情床上av| 午夜福利在线观看免费完整高清在 | 巨乳人妻的诱惑在线观看| 在线观看美女被高潮喷水网站 | 他把我摸到了高潮在线观看| 床上黄色一级片| 欧美成狂野欧美在线观看| 十八禁人妻一区二区| 免费在线观看成人毛片| 日本撒尿小便嘘嘘汇集6| 欧美中文综合在线视频| 黑人欧美特级aaaaaa片| 成人国产综合亚洲| 欧美日本视频| xxx96com| 九色国产91popny在线| 亚洲男人的天堂狠狠| 成人18禁在线播放| 亚洲精品乱码久久久v下载方式 | 一卡2卡三卡四卡精品乱码亚洲| netflix在线观看网站| 亚洲欧美日韩卡通动漫| 天天躁日日操中文字幕| 一本一本综合久久| 亚洲人与动物交配视频| 久久久国产成人免费| 国产亚洲精品久久久com| 最近最新中文字幕大全电影3| 免费在线观看影片大全网站| 亚洲精华国产精华精| 18禁黄网站禁片午夜丰满| 99久久99久久久精品蜜桃| 久久久久久九九精品二区国产| 久久久久久大精品| 美女高潮的动态| 中文资源天堂在线| 欧美+亚洲+日韩+国产| 日韩 欧美 亚洲 中文字幕| 美女大奶头视频| 一a级毛片在线观看| 国产91精品成人一区二区三区| 18禁国产床啪视频网站| 国产成人精品久久二区二区91| 美女cb高潮喷水在线观看 | 免费无遮挡裸体视频| 精品电影一区二区在线| 国产精品免费一区二区三区在线| 国产午夜精品久久久久久| 国产一区二区三区视频了| 国产精品香港三级国产av潘金莲| 亚洲欧美一区二区三区黑人| 蜜桃久久精品国产亚洲av| 亚洲人成网站高清观看| 亚洲美女黄片视频| 欧美日韩综合久久久久久 | 国产成人av教育| 青草久久国产| 精品午夜福利视频在线观看一区| 精品一区二区三区av网在线观看| 免费看光身美女| 亚洲 国产 在线| 成在线人永久免费视频| 亚洲第一欧美日韩一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成网站在线播放欧美日韩| 无遮挡黄片免费观看| 午夜福利高清视频| 亚洲成人精品中文字幕电影| 中文字幕人成人乱码亚洲影| 18禁国产床啪视频网站| 国产精品久久久久久人妻精品电影| 欧美中文综合在线视频| 免费无遮挡裸体视频| 国产精品一区二区三区四区免费观看 | 99视频精品全部免费 在线 | 亚洲无线在线观看| 免费看十八禁软件| 级片在线观看| 国产精品久久久久久久电影 | 久久国产精品影院| 男人舔女人的私密视频| 国产精品久久久人人做人人爽| 国产精品一区二区三区四区久久| 国产精品一区二区免费欧美| 美女被艹到高潮喷水动态| 又紧又爽又黄一区二区| 白带黄色成豆腐渣| 激情在线观看视频在线高清| 精品久久久久久,| 一个人免费在线观看的高清视频| 国产美女午夜福利| 高潮久久久久久久久久久不卡| 九色国产91popny在线| 日韩精品中文字幕看吧| 亚洲精品久久国产高清桃花| 免费看日本二区| 日本一本二区三区精品| 在线a可以看的网站| 欧美日本亚洲视频在线播放| 哪里可以看免费的av片| 日韩欧美精品v在线| 成人特级av手机在线观看| 日韩三级视频一区二区三区| www国产在线视频色| 国产成人福利小说| 美女高潮喷水抽搐中文字幕| 19禁男女啪啪无遮挡网站| 亚洲无线观看免费| 国产成人福利小说| 日韩欧美三级三区| 亚洲五月天丁香| 免费看光身美女| 午夜亚洲福利在线播放| 十八禁人妻一区二区| 欧美成人免费av一区二区三区| 中文字幕精品亚洲无线码一区| 国产成+人综合+亚洲专区| 他把我摸到了高潮在线观看| 国产精品亚洲av一区麻豆| 国产精品亚洲美女久久久| 真人一进一出gif抽搐免费| 欧美黄色淫秽网站| 精品福利观看| 午夜久久久久精精品| 日本黄色视频三级网站网址| 丁香六月欧美| 日本五十路高清| 成人永久免费在线观看视频| 免费看日本二区| 老熟妇乱子伦视频在线观看| 黑人欧美特级aaaaaa片| 毛片女人毛片| 国产1区2区3区精品| 村上凉子中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看| 日本在线视频免费播放| 色av中文字幕| 国产午夜福利久久久久久| 综合色av麻豆| 国产精品 欧美亚洲| 悠悠久久av| 久久人妻av系列| 精品人妻1区二区| 又黄又爽又免费观看的视频| 精品国产美女av久久久久小说| 老熟妇仑乱视频hdxx| 俄罗斯特黄特色一大片| 热99在线观看视频| 国产精华一区二区三区| 国产成+人综合+亚洲专区| 亚洲天堂国产精品一区在线| 999精品在线视频| 欧美成人性av电影在线观看| 免费观看精品视频网站| 亚洲人与动物交配视频| 国产毛片a区久久久久| 亚洲精品久久国产高清桃花| 国产激情偷乱视频一区二区| 亚洲美女视频黄频| 欧美日韩综合久久久久久 | 美女免费视频网站| 一本综合久久免费| 国产v大片淫在线免费观看| 亚洲精品国产精品久久久不卡| av天堂在线播放| 久久中文字幕一级| 我的老师免费观看完整版| 日本免费a在线| 国产精品乱码一区二三区的特点| bbb黄色大片| 久久精品国产清高在天天线| 中文字幕高清在线视频| 国产亚洲精品一区二区www| 熟女人妻精品中文字幕| 可以在线观看毛片的网站| 国产人伦9x9x在线观看| 女人被狂操c到高潮| 日韩有码中文字幕| 久久性视频一级片| a级毛片在线看网站| 香蕉丝袜av| 精品国产三级普通话版| 91九色精品人成在线观看| 欧美乱码精品一区二区三区| 欧美中文综合在线视频| 免费在线观看亚洲国产| 一进一出好大好爽视频| 深夜精品福利| 人妻丰满熟妇av一区二区三区| 国产精品野战在线观看| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 日本免费一区二区三区高清不卡| 国内精品久久久久精免费| 久久国产精品人妻蜜桃| 1024手机看黄色片| 88av欧美| 变态另类成人亚洲欧美熟女| 无限看片的www在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲成人久久性| 国产成人啪精品午夜网站| 午夜日韩欧美国产| 免费观看精品视频网站| 亚洲第一电影网av| aaaaa片日本免费| 国产高清激情床上av| 日本成人三级电影网站| 国产激情久久老熟女| 在线看三级毛片| 真实男女啪啪啪动态图| 观看美女的网站| 少妇裸体淫交视频免费看高清| 亚洲自偷自拍图片 自拍| 亚洲成av人片在线播放无| 久久国产精品影院| 色视频www国产| 国产一区二区激情短视频| 免费看美女性在线毛片视频| 99久久99久久久精品蜜桃| 两个人视频免费观看高清| 欧美在线一区亚洲| 99re在线观看精品视频| 免费观看精品视频网站| 在线播放国产精品三级| 深夜精品福利| 欧美+亚洲+日韩+国产| 99久久精品国产亚洲精品| 亚洲av电影不卡..在线观看| 国产精品亚洲av一区麻豆| 99热6这里只有精品| 久99久视频精品免费| 好男人在线观看高清免费视频| 久久精品91无色码中文字幕| 在线免费观看的www视频| 成人高潮视频无遮挡免费网站| 国内久久婷婷六月综合欲色啪| 国产欧美日韩一区二区精品| 桃色一区二区三区在线观看| 嫁个100分男人电影在线观看| www日本黄色视频网| 叶爱在线成人免费视频播放| netflix在线观看网站| 波多野结衣高清无吗| 国产成人av激情在线播放| 亚洲欧美精品综合久久99| 精品国产美女av久久久久小说| 九色国产91popny在线| 亚洲国产精品久久男人天堂| av女优亚洲男人天堂| 欧美另类亚洲清纯唯美| 青春草国产在线视频| 国产女主播在线喷水免费视频网站 | 亚洲美女视频黄频| 丰满少妇做爰视频| 天天躁日日操中文字幕| 久久久a久久爽久久v久久| 色综合色国产| 精品久久国产蜜桃| 男人狂女人下面高潮的视频| 日韩成人伦理影院| 国产成人91sexporn| 中文字幕精品亚洲无线码一区| kizo精华| 美女黄网站色视频| 麻豆乱淫一区二区| 欧美激情久久久久久爽电影| 少妇的逼好多水| 亚洲精品456在线播放app| 国产精品不卡视频一区二区| 久久久久精品久久久久真实原创| 久久久国产成人免费| 亚洲在久久综合| 99久久精品热视频| 最近2019中文字幕mv第一页| 久久这里有精品视频免费| 少妇丰满av| 99久久精品国产国产毛片| 国产午夜精品一二区理论片| 亚洲美女视频黄频| 久久久久网色| 亚洲国产精品合色在线| 日本猛色少妇xxxxx猛交久久| 久久久久久久国产电影| 国产爱豆传媒在线观看| 丝袜美腿在线中文| 级片在线观看| 欧美日韩在线观看h| 亚洲四区av| 一个人看视频在线观看www免费| 久久精品国产自在天天线| 国产精品熟女久久久久浪| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美三级三区| 国产精品蜜桃在线观看| 永久网站在线| 中国国产av一级| 日本三级黄在线观看| 成年女人看的毛片在线观看| 九九热线精品视视频播放| 女的被弄到高潮叫床怎么办| 免费观看在线日韩| 久久精品人妻少妇| 久久久久久大精品| 国产综合懂色| av播播在线观看一区| 国产人妻一区二区三区在| 久久久欧美国产精品| 国产高清视频在线观看网站| 欧美3d第一页| 欧美区成人在线视频| 长腿黑丝高跟| 欧美高清性xxxxhd video| 国产精品综合久久久久久久免费| 久久午夜福利片| 搞女人的毛片| 亚洲av免费在线观看| 少妇人妻精品综合一区二区| 黄色日韩在线| 精品久久国产蜜桃| 天堂影院成人在线观看| av在线蜜桃| 狂野欧美激情性xxxx在线观看| 亚洲欧美一区二区三区国产| 色5月婷婷丁香| 午夜a级毛片| 国产亚洲精品久久久com| 成人午夜精彩视频在线观看| 国产精品福利在线免费观看| 性插视频无遮挡在线免费观看| 国产又色又爽无遮挡免| 亚洲av.av天堂| 成人二区视频| 精品一区二区三区人妻视频| 午夜免费激情av| 亚洲人成网站在线观看播放| 色哟哟·www| 最近最新中文字幕大全电影3| av国产免费在线观看| 国产乱来视频区| 国产高潮美女av| 一本久久精品| 色综合色国产| 女人十人毛片免费观看3o分钟| 日韩在线高清观看一区二区三区| 高清av免费在线| 国产淫语在线视频| 免费电影在线观看免费观看| 久久久久久久国产电影| 五月伊人婷婷丁香| 午夜福利成人在线免费观看| .国产精品久久| 嫩草影院精品99| 成人午夜精彩视频在线观看| 精品久久久久久久久av| 亚洲人成网站高清观看| 国产精品久久久久久精品电影| videos熟女内射| 精品久久久久久电影网 | 我的老师免费观看完整版| 亚洲婷婷狠狠爱综合网| 一个人看视频在线观看www免费| 青青草视频在线视频观看| 伦精品一区二区三区| 国产精品三级大全| 天堂影院成人在线观看| 国产精品久久久久久精品电影| 欧美成人午夜免费资源| 午夜日本视频在线| 搡老妇女老女人老熟妇| 中文资源天堂在线| 欧美性猛交黑人性爽| 成人美女网站在线观看视频| 波多野结衣高清无吗| 国产老妇伦熟女老妇高清| 欧美三级亚洲精品| 午夜老司机福利剧场| 老司机影院毛片| 18禁裸乳无遮挡免费网站照片| h日本视频在线播放| 午夜福利在线观看吧| 精品国产露脸久久av麻豆 | 中文精品一卡2卡3卡4更新| 国产91av在线免费观看| 汤姆久久久久久久影院中文字幕 | 亚洲自偷自拍三级| 亚洲在线观看片| 亚洲熟妇中文字幕五十中出| 国产淫语在线视频| 久久精品综合一区二区三区| 在线播放国产精品三级| 亚洲欧美精品专区久久| 国产黄片视频在线免费观看| 大香蕉97超碰在线| 日韩国内少妇激情av| 一级毛片电影观看 | 免费一级毛片在线播放高清视频| 国产av不卡久久| 亚洲国产精品久久男人天堂| 少妇熟女aⅴ在线视频| 色播亚洲综合网| 国产高清不卡午夜福利| 国产人妻一区二区三区在| 成人三级黄色视频| 观看免费一级毛片| 激情 狠狠 欧美| av天堂中文字幕网| 简卡轻食公司| 22中文网久久字幕| 亚洲av电影不卡..在线观看| 黑人高潮一二区| 小蜜桃在线观看免费完整版高清| 久久精品久久精品一区二区三区| 在线a可以看的网站| 国产亚洲av嫩草精品影院| 亚洲av日韩在线播放| 麻豆av噜噜一区二区三区| 国产伦精品一区二区三区视频9| 亚洲国产精品久久男人天堂| 有码 亚洲区| 一级毛片aaaaaa免费看小|