• <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è)試研究
    看黄色毛片网站| 国产精品国产av在线观看| 热99国产精品久久久久久7| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧洲精品一区二区精品久久久| 中文字幕人妻熟女乱码| 女人久久www免费人成看片| 国产不卡av网站在线观看| 下体分泌物呈黄色| 99精品在免费线老司机午夜| 精品亚洲成国产av| 99久久国产精品久久久| 身体一侧抽搐| 久久午夜亚洲精品久久| 黄网站色视频无遮挡免费观看| 国产有黄有色有爽视频| 怎么达到女性高潮| 欧美午夜高清在线| 色精品久久人妻99蜜桃| 一边摸一边抽搐一进一出视频| av有码第一页| 国产成人免费无遮挡视频| 女人久久www免费人成看片| 中文字幕av电影在线播放| 国产精品久久久人人做人人爽| 大片电影免费在线观看免费| 日本欧美视频一区| 人妻久久中文字幕网| 国产精品香港三级国产av潘金莲| 久久草成人影院| 美女午夜性视频免费| 国产亚洲欧美精品永久| 国产真人三级小视频在线观看| 久久久久国产一级毛片高清牌| 免费在线观看黄色视频的| 久久香蕉国产精品| 男人舔女人的私密视频| 国产在视频线精品| 夫妻午夜视频| 午夜老司机福利片| 丝瓜视频免费看黄片| 99久久精品国产亚洲精品| 国产又爽黄色视频| 免费av中文字幕在线| 最近最新中文字幕大全免费视频| 亚洲五月天丁香| 中文欧美无线码| 午夜免费成人在线视频| 女人精品久久久久毛片| 欧美在线黄色| 777米奇影视久久| 午夜福利一区二区在线看| 成人国产一区最新在线观看| 亚洲国产欧美网| 中文字幕色久视频| 久久天躁狠狠躁夜夜2o2o| 亚洲第一青青草原| 12—13女人毛片做爰片一| 很黄的视频免费| 国产又色又爽无遮挡免费看| 丰满迷人的少妇在线观看| 日韩熟女老妇一区二区性免费视频| 美国免费a级毛片| 日韩人妻精品一区2区三区| 国产99白浆流出| 极品人妻少妇av视频| 日韩精品免费视频一区二区三区| 国产片内射在线| 黄色视频不卡| 无限看片的www在线观看| 国产精品久久久久久人妻精品电影| 97人妻天天添夜夜摸| tocl精华| 国产野战对白在线观看| 欧美中文综合在线视频| 天堂动漫精品| 好看av亚洲va欧美ⅴa在| 最近最新免费中文字幕在线| 美女视频免费永久观看网站| 人人妻人人澡人人看| 丁香欧美五月| 久久国产精品人妻蜜桃| videosex国产| 欧美日韩瑟瑟在线播放| 午夜精品在线福利| 色在线成人网| 亚洲熟妇熟女久久| 高清视频免费观看一区二区| 最新在线观看一区二区三区| 嫩草影视91久久| 国产国语露脸激情在线看| 国产精品1区2区在线观看. | 丰满饥渴人妻一区二区三| 99精品久久久久人妻精品| tube8黄色片| xxx96com| 日韩欧美免费精品| 一边摸一边做爽爽视频免费| 热99re8久久精品国产| 午夜影院日韩av| 色尼玛亚洲综合影院| 捣出白浆h1v1| 日韩免费高清中文字幕av| 老鸭窝网址在线观看| 黄色 视频免费看| 在线av久久热| 亚洲成人国产一区在线观看| 色综合婷婷激情| 国产黄色免费在线视频| 亚洲伊人色综图| 欧美精品av麻豆av| 国产一区二区三区在线臀色熟女 | 国产精品 欧美亚洲| 成年人免费黄色播放视频| 91在线观看av| 久久精品成人免费网站| 丝袜在线中文字幕| 亚洲综合色网址| 在线天堂中文资源库| 超碰成人久久| 操出白浆在线播放| 久久久精品免费免费高清| 黄色毛片三级朝国网站| 久久久久久久精品吃奶| 中文字幕制服av| www.熟女人妻精品国产| 久久中文字幕一级| 热re99久久精品国产66热6| 中出人妻视频一区二区| 精品人妻熟女毛片av久久网站| 色老头精品视频在线观看| 视频区图区小说| 亚洲视频免费观看视频| 又黄又粗又硬又大视频| 国产熟女午夜一区二区三区| 一区福利在线观看| 丝袜美足系列| 国产亚洲欧美98| 欧美大码av| 99久久综合精品五月天人人| 女警被强在线播放| 欧美亚洲 丝袜 人妻 在线| 高清黄色对白视频在线免费看| 国产伦人伦偷精品视频| 色在线成人网| 韩国av一区二区三区四区| 女警被强在线播放| 中文字幕人妻丝袜一区二区| 黄色成人免费大全| 亚洲综合色网址| 欧美日韩乱码在线| 人人妻,人人澡人人爽秒播| 中文字幕av电影在线播放| 两性夫妻黄色片| 香蕉丝袜av| 久久久国产精品麻豆| 日本黄色日本黄色录像| 黄片小视频在线播放| 在线观看免费视频网站a站| 色在线成人网| 久久久精品免费免费高清| 国产av精品麻豆| 欧美+亚洲+日韩+国产| 国精品久久久久久国模美| 国产av又大| 一本一本久久a久久精品综合妖精| 亚洲专区国产一区二区| 夫妻午夜视频| 一级,二级,三级黄色视频| 精品国产一区二区三区四区第35| 欧美乱色亚洲激情| 精品久久久久久电影网| 国产成人精品在线电影| 成人国产一区最新在线观看| 国产1区2区3区精品| av福利片在线| 19禁男女啪啪无遮挡网站| 飞空精品影院首页| 亚洲精品久久午夜乱码| tocl精华| 一二三四社区在线视频社区8| 亚洲精品国产区一区二| 欧美大码av| 欧美精品一区二区免费开放| 99热只有精品国产| 久久久精品区二区三区| 国产aⅴ精品一区二区三区波| 国产精品久久视频播放| av视频免费观看在线观看| 777久久人妻少妇嫩草av网站| 法律面前人人平等表现在哪些方面| 一级片免费观看大全| 一区二区日韩欧美中文字幕| 国产精品自产拍在线观看55亚洲 | av免费在线观看网站| 国产高清视频在线播放一区| 热99久久久久精品小说推荐| 亚洲精品中文字幕一二三四区| 欧美精品高潮呻吟av久久| 亚洲自偷自拍图片 自拍| 久久久久视频综合| 男女免费视频国产| av在线播放免费不卡| 亚洲色图av天堂| av天堂久久9| 国产免费男女视频| 欧美一级毛片孕妇| 日韩人妻精品一区2区三区| 久久天堂一区二区三区四区| 丰满饥渴人妻一区二区三| 一级毛片女人18水好多| 国产精品美女特级片免费视频播放器 | 免费观看a级毛片全部| 亚洲av成人av| 日韩免费av在线播放| 国产亚洲精品一区二区www | 九色亚洲精品在线播放| 成人18禁高潮啪啪吃奶动态图| 亚洲色图 男人天堂 中文字幕| 在线天堂中文资源库| 男女午夜视频在线观看| www日本在线高清视频| 伦理电影免费视频| tocl精华| 狠狠婷婷综合久久久久久88av| 搡老岳熟女国产| 欧美精品亚洲一区二区| 亚洲精品国产色婷婷电影| 国产熟女午夜一区二区三区| 天天添夜夜摸| 国产xxxxx性猛交| 亚洲专区国产一区二区| 一级毛片女人18水好多| 日本撒尿小便嘘嘘汇集6| 亚洲av欧美aⅴ国产| 日本精品一区二区三区蜜桃| 亚洲国产欧美一区二区综合| 丝袜人妻中文字幕| 久久久久视频综合| 99热国产这里只有精品6| 国产亚洲精品久久久久5区| 777米奇影视久久| 天天操日日干夜夜撸| 亚洲男人天堂网一区| 精品人妻1区二区| 精品亚洲成国产av| 悠悠久久av| av一本久久久久| 18在线观看网站| 久久精品人人爽人人爽视色| av欧美777| 999久久久国产精品视频| 9热在线视频观看99| 亚洲色图综合在线观看| 中出人妻视频一区二区| 欧美成狂野欧美在线观看| 国产精品av久久久久免费| 亚洲精品一二三| 欧美av亚洲av综合av国产av| 国产一卡二卡三卡精品| 国产精品九九99| 男女之事视频高清在线观看| 丝袜人妻中文字幕| 亚洲人成伊人成综合网2020| svipshipincom国产片| 在线av久久热| 18禁黄网站禁片午夜丰满| 成年人黄色毛片网站| 国产国语露脸激情在线看| 18禁裸乳无遮挡免费网站照片 | xxxhd国产人妻xxx| 法律面前人人平等表现在哪些方面| 国产麻豆69| 亚洲五月婷婷丁香| 69av精品久久久久久| 日韩熟女老妇一区二区性免费视频| 久久亚洲真实| 在线观看日韩欧美| 亚洲综合色网址| 高潮久久久久久久久久久不卡| 久久精品国产综合久久久| 久久国产精品影院| www.自偷自拍.com| 老司机影院毛片| 黄色片一级片一级黄色片| 在线观看免费日韩欧美大片| a级毛片在线看网站| 国产精品亚洲av一区麻豆| 夜夜爽天天搞| 国产亚洲精品久久久久5区| 久久久久久久午夜电影 | 国产精品98久久久久久宅男小说| av电影中文网址| 午夜日韩欧美国产| 最新美女视频免费是黄的| 欧美另类亚洲清纯唯美| 亚洲精品乱久久久久久| 国产成人一区二区三区免费视频网站| 91国产中文字幕| 麻豆av在线久日| 人成视频在线观看免费观看| 免费观看人在逋| 老司机影院毛片| av有码第一页| www.999成人在线观看| 久久狼人影院| 在线观看一区二区三区激情| 国产深夜福利视频在线观看| 妹子高潮喷水视频| 欧美成人午夜精品| 手机成人av网站| 超碰成人久久| 一本综合久久免费| netflix在线观看网站| 丁香六月欧美| 免费观看精品视频网站| 人人澡人人妻人| 国产日韩一区二区三区精品不卡| 高清在线国产一区| 伦理电影免费视频| 国产亚洲av高清不卡| 国产精品二区激情视频| bbb黄色大片| 好看av亚洲va欧美ⅴa在| 两性午夜刺激爽爽歪歪视频在线观看 | 人人澡人人妻人| 亚洲精品美女久久av网站| 黄色女人牲交| 亚洲精华国产精华精| 一级片免费观看大全| 交换朋友夫妻互换小说| 精品免费久久久久久久清纯 | 国产片内射在线| 伦理电影免费视频| 国产片内射在线| 最新美女视频免费是黄的| 国产在线精品亚洲第一网站| 国产精品一区二区免费欧美| 国产人伦9x9x在线观看| 日韩中文字幕欧美一区二区| 亚洲国产欧美日韩在线播放| 每晚都被弄得嗷嗷叫到高潮| 三上悠亚av全集在线观看| 一进一出抽搐gif免费好疼 | 黑人操中国人逼视频| 日韩免费高清中文字幕av| 99国产精品一区二区三区| 欧美激情久久久久久爽电影 | 成年动漫av网址| 午夜福利一区二区在线看| 69精品国产乱码久久久| www日本在线高清视频| 亚洲久久久国产精品| 老司机午夜福利在线观看视频| 精品国产国语对白av| 国产免费男女视频| 两性午夜刺激爽爽歪歪视频在线观看 | 性少妇av在线| 最新在线观看一区二区三区| 色综合欧美亚洲国产小说| 国产精品欧美亚洲77777| www.熟女人妻精品国产| 精品久久久久久久久久免费视频 | 午夜精品久久久久久毛片777| av电影中文网址| 80岁老熟妇乱子伦牲交| 午夜福利,免费看| 五月开心婷婷网| 国产精品av久久久久免费| 大型av网站在线播放| 中亚洲国语对白在线视频| 久久精品亚洲av国产电影网| 欧美日韩一级在线毛片| 国产亚洲欧美98| www.精华液| 一级毛片女人18水好多| 午夜91福利影院| 女性被躁到高潮视频| 黑人巨大精品欧美一区二区mp4| 成人特级黄色片久久久久久久| 18禁观看日本| 亚洲国产精品sss在线观看 | 人人澡人人妻人| netflix在线观看网站| 成年人午夜在线观看视频| 亚洲中文字幕日韩| 日韩欧美国产一区二区入口| 巨乳人妻的诱惑在线观看| 午夜免费观看网址| 成人手机av| 丁香六月欧美| 亚洲 国产 在线| 美女高潮到喷水免费观看| 99久久国产精品久久久| 成人精品一区二区免费| 老熟女久久久| 国产有黄有色有爽视频| 免费在线观看亚洲国产| 正在播放国产对白刺激| 在线国产一区二区在线| 又紧又爽又黄一区二区| av天堂在线播放| 宅男免费午夜| 香蕉丝袜av| 成人黄色视频免费在线看| 成年女人毛片免费观看观看9 | 9191精品国产免费久久| 久热爱精品视频在线9| 国产精品影院久久| 一进一出好大好爽视频| 精品国产一区二区三区久久久樱花| 中国美女看黄片| 国产在线精品亚洲第一网站| 女人被躁到高潮嗷嗷叫费观| 国产三级黄色录像| 精品国产一区二区久久| 在线观看舔阴道视频| 精品久久久久久,| 51午夜福利影视在线观看| 欧美中文综合在线视频| 亚洲熟女毛片儿| 成人免费观看视频高清| 国产成人影院久久av| 成人三级做爰电影| 成人精品一区二区免费| 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 国产免费现黄频在线看| 青草久久国产| 亚洲熟妇熟女久久| 校园春色视频在线观看| 日韩欧美国产一区二区入口| 中亚洲国语对白在线视频| 老司机深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全电影3 | 精品国产超薄肉色丝袜足j| 露出奶头的视频| 99久久国产精品久久久| 久久影院123| 在线播放国产精品三级| 中文亚洲av片在线观看爽 | 精品久久久久久电影网| 欧美乱色亚洲激情| 久热这里只有精品99| 日韩欧美一区二区三区在线观看 | 久热爱精品视频在线9| 久久久久国产精品人妻aⅴ院 | 三上悠亚av全集在线观看| 国产亚洲一区二区精品| 又紧又爽又黄一区二区| 一二三四在线观看免费中文在| 老司机深夜福利视频在线观看| 99久久99久久久精品蜜桃| 久久 成人 亚洲| 一级毛片女人18水好多| 黄色怎么调成土黄色| 国产在线精品亚洲第一网站| 国产日韩一区二区三区精品不卡| av一本久久久久| 国产人伦9x9x在线观看| www.熟女人妻精品国产| 丰满人妻熟妇乱又伦精品不卡| 日日夜夜操网爽| 热re99久久精品国产66热6| 黄色女人牲交| 久久久久久久精品吃奶| 亚洲精品乱久久久久久| 99久久99久久久精品蜜桃| 中文字幕制服av| 中文字幕人妻丝袜一区二区| 亚洲av电影在线进入| 久久精品亚洲av国产电影网| av一本久久久久| 欧美日韩精品网址| 国产在线观看jvid| 国产有黄有色有爽视频| 亚洲欧美激情综合另类| 看片在线看免费视频| 国产高清视频在线播放一区| 欧美黑人欧美精品刺激| 亚洲五月婷婷丁香| av天堂久久9| 国产精品亚洲一级av第二区| 亚洲精品久久成人aⅴ小说| 交换朋友夫妻互换小说| 大片电影免费在线观看免费| 欧美国产精品va在线观看不卡| 91成人精品电影| 99国产综合亚洲精品| 两人在一起打扑克的视频| 亚洲国产欧美网| 十八禁高潮呻吟视频| 老司机影院毛片| 夜夜夜夜夜久久久久| 黑人操中国人逼视频| 涩涩av久久男人的天堂| 叶爱在线成人免费视频播放| 国产成人免费观看mmmm| 午夜福利影视在线免费观看| 久久国产乱子伦精品免费另类| 国产有黄有色有爽视频| 久久精品国产99精品国产亚洲性色 | 一个人免费在线观看的高清视频| 精品熟女少妇八av免费久了| 久久国产亚洲av麻豆专区| 欧美+亚洲+日韩+国产| 免费高清在线观看日韩| 欧美丝袜亚洲另类 | 国产人伦9x9x在线观看| 成人特级黄色片久久久久久久| 高清视频免费观看一区二区| 麻豆国产av国片精品| 看片在线看免费视频| 新久久久久国产一级毛片| 中文字幕人妻丝袜制服| 婷婷丁香在线五月| 老司机在亚洲福利影院| 另类亚洲欧美激情| 黄片大片在线免费观看| 人人澡人人妻人| 亚洲 国产 在线| 91av网站免费观看| 欧美av亚洲av综合av国产av| 高清黄色对白视频在线免费看| 国产精品 欧美亚洲| 丰满饥渴人妻一区二区三| 日韩三级视频一区二区三区| 精品国产一区二区久久| 欧美 亚洲 国产 日韩一| 国产av一区二区精品久久| 成年人黄色毛片网站| 中文字幕高清在线视频| 搡老岳熟女国产| 国产成人影院久久av| 亚洲色图综合在线观看| 涩涩av久久男人的天堂| 亚洲aⅴ乱码一区二区在线播放 | 色综合婷婷激情| 国产一区在线观看成人免费| 欧美亚洲 丝袜 人妻 在线| 国精品久久久久久国模美| 男女午夜视频在线观看| 亚洲国产欧美一区二区综合| 亚洲精品久久午夜乱码| 新久久久久国产一级毛片| 欧美精品人与动牲交sv欧美| 亚洲在线自拍视频| 女人高潮潮喷娇喘18禁视频| 高清黄色对白视频在线免费看| 久久人妻av系列| 极品少妇高潮喷水抽搐| 国产麻豆69| 亚洲专区中文字幕在线| 久久九九热精品免费| 九色亚洲精品在线播放| 性色av乱码一区二区三区2| 国产激情欧美一区二区| 久久中文看片网| 日韩欧美在线二视频 | 日韩有码中文字幕| 免费日韩欧美在线观看| 国产欧美日韩精品亚洲av| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜精品| 久久中文字幕一级| xxx96com| 99re在线观看精品视频| 在线观看日韩欧美| 精品福利观看| 十八禁网站免费在线| av网站在线播放免费| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看| 国产99久久九九免费精品| 日韩 欧美 亚洲 中文字幕| 在线观看免费午夜福利视频| 下体分泌物呈黄色| 女人久久www免费人成看片| 欧洲精品卡2卡3卡4卡5卡区| 韩国精品一区二区三区| 这个男人来自地球电影免费观看| 91九色精品人成在线观看| 亚洲久久久国产精品| 亚洲情色 制服丝袜| 精品国产美女av久久久久小说| 多毛熟女@视频| 精品国产一区二区三区四区第35| 久久久国产成人精品二区 | 激情视频va一区二区三区| a级毛片在线看网站| 老司机在亚洲福利影院| 亚洲av欧美aⅴ国产| 又黄又爽又免费观看的视频| 90打野战视频偷拍视频| 淫妇啪啪啪对白视频| 欧美黑人精品巨大| 欧美黄色淫秽网站| 国产高清视频在线播放一区| 视频在线观看一区二区三区| 高清欧美精品videossex| 久久香蕉国产精品| 久久国产精品人妻蜜桃| 精品视频人人做人人爽| 两人在一起打扑克的视频| 女人被躁到高潮嗷嗷叫费观| 国产黄色免费在线视频| 女性被躁到高潮视频| 久久中文字幕人妻熟女| 啦啦啦视频在线资源免费观看| 宅男免费午夜| 日韩欧美在线二视频 | 两性夫妻黄色片| 亚洲专区国产一区二区| 高清黄色对白视频在线免费看|