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

    炸藥顆粒床沖擊加熱多尺度數(shù)值計(jì)算

    2013-02-23 06:44:54張新明吳艷青黃風(fēng)雷
    兵工學(xué)報(bào) 2013年2期
    關(guān)鍵詞:構(gòu)形內(nèi)能活塞

    張新明,吳艷青,黃風(fēng)雷

    (北京理工大學(xué) 爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100081)

    0 引言

    塑料黏結(jié)炸藥(PBXs)可看作是一種特殊的顆粒復(fù)合材料,其主要組成部分包括HMX 含能顆粒(質(zhì)量分?jǐn)?shù)約占95%)、黏結(jié)劑、增塑劑及鈍感劑等(質(zhì)量分?jǐn)?shù)約占5%).PBXs 在受到?jīng)_擊加載時(shí),HMX 炸藥顆粒為主要的承載對(duì)象,因此,在進(jìn)行PBXs 沖擊起爆研究時(shí),為了簡(jiǎn)化模型,有時(shí)忽略黏結(jié)劑、增塑劑等的作用,將PBXs 看成由HMX 炸藥顆粒構(gòu)成的顆粒床進(jìn)行研究。另外,對(duì)于無空隙的均質(zhì)含能材料,由于化學(xué)構(gòu)成的不穩(wěn)定性、老化、損傷等因素會(huì)出現(xiàn)孔隙,研究這類材料的性能,也經(jīng)常將它們看成顆粒炸藥床進(jìn)行研究[1]。

    燃燒轉(zhuǎn)爆轟(DDT)的研究已有幾十年的歷史[2-10],也是顆粒含能材料受沖擊加載時(shí)常常產(chǎn)生的一個(gè)現(xiàn)象。文獻(xiàn)[3 -8]均表明顆粒含能材料的燃燒轉(zhuǎn)爆轟是顆粒材料的壓實(shí)(Material Compaction)和熱量釋放(Heat Release)的相互作用,而多數(shù)情況下,燃燒轉(zhuǎn)爆轟都與致密波(Compaction Wave)作用密切相關(guān)[11]。所謂致密波,即為使顆粒床孔隙率發(fā)生變化的行進(jìn)波[11],Sandusky 等[12]最早在試驗(yàn)中觀察到穩(wěn)態(tài)致密波。在致密區(qū)內(nèi),由于顆粒與顆粒之間的摩擦、應(yīng)力集中、塑性變形等原因會(huì)導(dǎo)致材料的局部溫度大幅升高,形成所謂的熱點(diǎn)(Hot-spot),在整體平均溫度很低,甚至還沒達(dá)到熔點(diǎn)的情況下,就有可能使材料局部熔化,發(fā)生相變,激起劇烈的化學(xué)反應(yīng),并最終導(dǎo)致爆燃甚至是爆轟。自上世紀(jì)80年代始,國(guó)內(nèi)外很多學(xué)者開展了這方面的實(shí)驗(yàn)、理論和數(shù)值研究[13-19]。

    但宏觀尺度的整體平均模型[13-17]不能有效地反映出顆粒尺度的響應(yīng),而基于單個(gè)顆粒的統(tǒng)計(jì)模型[18-19]由于計(jì)算量太大,不適合工程計(jì)算。近年來,Gonthier 等[20]根據(jù)Hertz 接觸力學(xué)建立了顆粒材料沖擊加熱的多尺度模型,該模型能夠有效地反映出顆粒尺度的響應(yīng),且適合大量工程計(jì)算。但是,Gonthier[20]沒有考慮顆??蓧嚎s性和固液相變對(duì)壓縮過程的影響,并且顆粒構(gòu)形應(yīng)力函數(shù)由準(zhǔn)靜態(tài)壓縮試驗(yàn)測(cè)得,工作量大。本文選取Tait 狀態(tài)方程,由顆粒構(gòu)形應(yīng)力的定義式來確定該函數(shù),避免進(jìn)行大量試驗(yàn),從而減小工作量;同時(shí)引入化學(xué)反應(yīng)進(jìn)程變量,并考慮炸藥顆粒的可壓縮性以及固液相變的影響,對(duì)活塞沖擊顆粒床進(jìn)行多尺度數(shù)值計(jì)算。

    1 計(jì)算模型

    Baer 等[7]的兩相流模型,通常稱為BN 模型,成為國(guó)內(nèi)外學(xué)者研究致密波及含能材料燃燒轉(zhuǎn)爆轟基礎(chǔ),本文在BN 模型的基礎(chǔ)上進(jìn)行簡(jiǎn)化,并結(jié)合顆粒尺度熱點(diǎn)模型和阿倫尼烏斯反應(yīng)速率方程,對(duì)活塞沖擊顆粒床進(jìn)行數(shù)值計(jì)算。

    1.1 一維單相流模型

    根據(jù)Baer 等[13]和Lovey 等[21]的工作,氣相對(duì)致密波的影響很小,本文忽略氣相作用,一維單相流模型如下:

    式中:t 為時(shí)間;x 為位置;ρs為固相密度;φs為固相體積分?jǐn)?shù);vs為顆粒速度;ps為壓強(qiáng);E 為單位體積炸藥的總能量,其中E =es+B,es為單位體積炸藥的內(nèi)能,B 為單位體積炸藥的壓縮潛能;β 為顆粒構(gòu)形應(yīng)力函數(shù),表示顆粒床受到擠壓后,顆粒之間的相互作用力;常數(shù)μc為壓縮粘性。(1)~(3)式分別為質(zhì)量、動(dòng)量和能量守恒方程;(4)式表示由于固相壓力與顆粒構(gòu)形應(yīng)力存在壓差而引起的固相體積分?jǐn)?shù)的變化。

    1.2 顆粒尺度的熱點(diǎn)模型

    圖1為Gonthier[20]根據(jù)Hertz 接觸力學(xué)建立的熱點(diǎn)模型。圖中,rc為接觸球半徑,r0為接觸球半徑的最大值。當(dāng)顆粒床受到外載荷時(shí),顆粒與顆粒之間由于摩擦、應(yīng)力集中、塑性變形等因素使得顆粒接觸面附近溫度升高,形成所謂的“熱點(diǎn)”,在整體平均溫度很低的情況下有可能發(fā)生燃燒轉(zhuǎn)爆轟。

    圖1 顆粒與顆粒間接觸的幾何模型Fig.1 Illustration of the intergranular contact geometry used for the thermal energy localization strategy

    為簡(jiǎn)化模型,假設(shè)HMX 顆粒不在發(fā)生破碎和聚并,顆粒半徑為R,單位體積內(nèi)的顆粒數(shù)n,顆粒數(shù)守恒方程為

    本文忽略顆粒與顆粒之間的摩擦作用,單位體積內(nèi)塑性功功率等于內(nèi)能的變化率,

    式中:nc為單位體積內(nèi)的接觸球個(gè)數(shù);Vc=4/3r3c為單個(gè)接觸球的體積;σY=3.0Y 為塑形應(yīng)力,Y 為HMX的屈服強(qiáng)度;des/dt 為內(nèi)能變化率。

    阿倫尼烏斯反應(yīng)速率方程[20]

    式中:λ(x,r,t)(0≤λ≤1)為沿接觸球徑向(r 方向,見圖1)變化的化學(xué)反應(yīng)過程變量,λ =0 表示無化學(xué)反應(yīng)發(fā)生,λ=1 表示完全反應(yīng);z 為阿倫尼烏斯因子;Tα為活化溫度。

    內(nèi)能變化率

    式中:q 為傳導(dǎo)熱通量,q = -k ?T(x,r,t)/?r,k 為熱傳導(dǎo)率;ρs0為固相初始密度。(8)式右端第1 項(xiàng)為由熱傳導(dǎo)引起的內(nèi)能變化率,第2 項(xiàng)為由化學(xué)反應(yīng)引起的內(nèi)能變化率,第3 項(xiàng)為由顆粒材料孔隙率變化引起的內(nèi)能變化率Sφ以及由顆粒密度變化引起的內(nèi)能變化率Sρ.

    HMX 在常溫常壓下的熔點(diǎn)T0m為552 K,根據(jù)Kraut-Kennedy[22]關(guān)系,熔點(diǎn)Tm= T0m(1 +αΔV/V0),其中,常數(shù)α =1.53,V 為比容,V0為初始比容;在ps≈5 GPa 的條件下,ΔV≈0.1,熔點(diǎn)升高到約636 K,而本文條件下ps?5 GPa,因此假設(shè)Tm=T0m,顆粒床發(fā)生等溫相變:

    式中:cvs為定容比熱;q0m為熔化潛熱;0≤χ≤1 為液相質(zhì)量分?jǐn)?shù)。

    1.3 本構(gòu)關(guān)系

    為使方程組封閉,還需要給出固相狀態(tài)方程和顆粒構(gòu)形應(yīng)力函數(shù)β.本文對(duì)Powers 等[23]給出的狀態(tài)方程進(jìn)行修改,引入化學(xué)反應(yīng)進(jìn)程變量λs,Helmholtz 自由能的表達(dá)式為:

    式中:Ts0為固相初始溫度;ps0為固相初始?jí)簭?qiáng);pg0為氣相初始?jí)簭?qiáng);φs0為固相初始體積分?jǐn)?shù)。

    狀態(tài)方程

    顆粒構(gòu)形應(yīng)力函數(shù)β 的表達(dá)式如下:

    Tait 狀態(tài)方程中的參數(shù)如表1所示。

    表1 Tait 狀態(tài)方程中的參數(shù)Tab.1 Parameters in Tait EOS for solid HMX

    2 計(jì)算模型的無量綱化

    本文以HMX 顆粒為例,對(duì)其致密波傳播特性進(jìn)行分析,表2為HMX 的物性參數(shù)及初始條件。

    表2 HMX 物理性質(zhì)參數(shù)及初始條件Tab.2 HMX material properties and initial conditions

    假設(shè)致密波以波速D 向右傳播,利用伽利略坐標(biāo)變換公式:相對(duì)位置ξ=x-D·t 和相對(duì)速度v=vs-D 以及定常假設(shè),?/?t ≡0.定義如下無量綱量:

    將計(jì)算模型寫成無量綱形式

    3 計(jì)算結(jié)果

    為確定壓縮過程的終態(tài),令β =ps,即顆粒床達(dá)到壓力平衡。圖2為顆粒構(gòu)形應(yīng)力函數(shù)β 與固相體積分?jǐn)?shù)φs之間的關(guān)系。由圖可知,β 為體積分?jǐn)?shù)φs的單調(diào)遞增函數(shù)。圖3為不同致密波速D,固相壓力Ps與體積分?jǐn)?shù)φs之間的關(guān)系。由圖可知,隨著致密波速度的增加,終態(tài)體積分?jǐn)?shù)φs逐漸增加,當(dāng)D >500 m/s 時(shí),終態(tài)體積分?jǐn)?shù)φs≈1,即顆粒床處于完全壓實(shí)狀態(tài)。

    圖2 顆粒構(gòu)形應(yīng)力β 與固相體積分?jǐn)?shù)φs之間的關(guān)系Fig.2 Intergranular stress β vs.solid volume fraction φs

    圖3 固相壓力ps與固相體積分?jǐn)?shù)φs之間的關(guān)系Fig.3 Solid pressure ps vs.solid volume fraction φs

    3.1 惰性顆粒藥床

    圖4 初始體積分?jǐn)?shù)φs0 =0.76,活塞速度up =100 m/s 時(shí)模型物理量在ξ 方向上的分布Fig.4 Variation of model parameters in ξ direction for φs0 =0.76 and up =100 m/s

    圖4為活塞速度up=100 m/s 撞擊顆粒藥床,產(chǎn)生亞音速致密波(D=418.7 m/s <C0),固相體積分?jǐn)?shù)φs、壓力ps、顆粒間的應(yīng)力β、整體平均溫度Ts和顆粒密度ρs的分布。由圖4可知,各物理量由初始狀態(tài)光滑地過渡到平衡狀態(tài),終態(tài)值φs=0.995,ps=67.73 MPa,Ts=304.08 K,ρs=1 906.62 kg/m3.由圖4(b)可知,在初始階段顆粒床處于平衡狀態(tài),即Ps0=β(φs0),隨著活塞向右運(yùn)動(dòng),壓力增加,顆粒間構(gòu)形應(yīng)力增加,當(dāng)ps=β(φs)時(shí),顆粒床再次達(dá)到壓力平衡狀態(tài)。

    圖5為up=100 m/s 撞擊顆粒床時(shí),顆粒半徑R、接觸球半徑rc和顆粒溫度T 的分布。由圖5(a)可知,壓縮過程中顆粒半徑R 幾乎保持不變,接觸球半徑rc由初始值rc0≈1.859 8 μm 光滑地過渡到穩(wěn)態(tài)值rc≈4.849 4 μm;由圖5(b)可知,由于顆粒之間的擠壓,發(fā)生塑性變形,在顆粒接觸面附近,顆粒溫度明顯升高,最大值T≈973 K,超過了HMX 的熔點(diǎn)T0m=552 K,并可以觀察到相變發(fā)生。而由圖4(c)可知,相同的加載條件下,整體平均溫度Ts≈304 K,遠(yuǎn)低于顆粒尺度溫度的最大值,即整體平均模型不能有效地反應(yīng)顆粒尺度的響應(yīng)。

    圖5 顆粒半徑和顆粒尺度溫度的分布(up =100 m/s)Fig.5 Variation of grain radius in ξ direction and grainscale temperature in ξ and radial directions at up =100 m/s

    3.2 考慮化學(xué)反應(yīng)

    圖6為考慮固相可壓縮性和相變的情況下,化學(xué)反應(yīng)進(jìn)程變量λ、顆粒溫度T、整體平均反應(yīng)進(jìn)程變量λs和整體平均溫度Ts的分布。由(19)式可知,由于考慮了固相可壓縮性,接觸球半徑rc會(huì)增大,這就意味著內(nèi)能將會(huì)分布在顆粒接觸面附近更大的區(qū)域上,從而使顆粒溫度降低;另外由于相變會(huì)消耗掉一部分能量,也將使顆粒溫度降低,因此在相同的初始條件下要激發(fā)化學(xué)反應(yīng)必須增大活塞速度。當(dāng)活塞速度up=79.0 m/s 時(shí),致密波速度D=348.75 m/s,由圖6(a)和圖6(b)可知,在顆粒接觸面附近化學(xué)反應(yīng)進(jìn)程變量λ =1,即完全反應(yīng),相應(yīng)地溫度出現(xiàn)了突跳,顆粒溫度的最大值為1 638 K,并且可以觀察到相變的發(fā)生。由圖6(c)和圖6(d)可知,整體平均化學(xué)反應(yīng)進(jìn)程變量λ≈7.5 ×10-7,整體平均溫度Ts≈302.5 K.結(jié)果再次表明整體平均模型不能有效地反應(yīng)顆粒尺度的響應(yīng)。

    3.3 參數(shù)研究

    圖7為活塞速度up、初始體積分?jǐn)?shù)φs0、顆粒半徑R 等模型參數(shù)對(duì)顆粒床壓縮過程的影響。由圖7(a)可知,隨著活塞速度up的增加,固相體積分?jǐn)?shù)φs逐漸增加,up~φs曲線的斜率逐漸減小,表明顆粒床更難被壓實(shí);當(dāng)活塞速度up>100 m/s 時(shí),固相體積分?jǐn)?shù)φs≈1,即顆粒藥床被完全壓實(shí)。圖7(b)為活塞速度up與致密波速度D 之間的關(guān)系,由圖可知,隨著活塞速度up的增加,致密波速度D 逐漸增加,且兩者之間近似成線性關(guān)系。由圖7(c)可知,隨著初始體積分?jǐn)?shù)φs0的增加,要想激起化學(xué)反應(yīng),所需壓力逐漸增加;當(dāng)φs0>0.90 時(shí),所需壓力急劇增加,也即致密的顆粒床更難發(fā)生化學(xué)反應(yīng)。由圖7(d)可知,隨著顆粒半徑R 的增加,激起化學(xué)反應(yīng)所需的活塞速度逐漸減小,即粗顆粒床更容易發(fā)生化學(xué)反應(yīng),這一結(jié)果與已有的試驗(yàn)和數(shù)值模擬結(jié)果[25-28]相一致。

    4 結(jié)論

    本文對(duì)活塞沖擊HMX 顆粒床進(jìn)行多尺度數(shù)值模擬。結(jié)果表明,對(duì)于亞音速致密波,固相體積分?jǐn)?shù)、壓力、密度和整體平均溫度光滑地從初始值過渡到平衡狀態(tài),致密波速度隨著活塞速度的增加而逐漸增加,顆粒床的致密程度隨著致密波速度的增加而增加,直到藥床被完全壓實(shí)(φs≈1).顆粒尺度的溫度分布表明,在相同的加載條件下(up=100 m/s),顆粒尺度的溫度明顯高于整體平均溫度,并超過了熔點(diǎn),同時(shí)可以觀察到固液相變的發(fā)生,結(jié)果表明整體平均模型不能有效地反映出材料顆粒尺度的響應(yīng)。在考慮化學(xué)反應(yīng)的情況下,由于顆粒材料可壓縮性和固液相變的影響,在相同的初始條件下,必須增大活塞速度才能激起化學(xué)反應(yīng)。計(jì)算結(jié)果還表明,隨著固相初始體積分?jǐn)?shù)的增大,發(fā)生化學(xué)反應(yīng)所需的外界載荷越大,即致密床更難發(fā)生化學(xué)反應(yīng);隨著顆粒半徑的增大,發(fā)生化學(xué)反應(yīng)所需的活塞速度越小,即粗顆粒床更容易發(fā)生化學(xué)反應(yīng)。

    圖6 化學(xué)反應(yīng)進(jìn)程變量和溫度的分布Fig.6 Predicted chemical reaction progress variable and temperature through the compaction zone

    圖7 模型參數(shù)研究Fig.7 Parametric sensitivity analysis results

    References)

    [1] Greenaway M W.Measurement of intergranular stress and porosity during dynamic compaction of porous beds of HMX[J].Journal of Applied Physics,2005,97(9):093521 -093521-7.

    [2] Gonthier K A,Powers J M.A high resolution numerical method for a two-phase model of deflagration-to-detonation transition[J].Journal of Computational Physics,2000,163(2):376 -433.

    [3] Cochran M T,Powers J M.Computation of compaction in compressible granular material[J].Mechanics Research Communications,2008,35(1 -2):96 -103.

    [4] Bernecker R R,Price D.Studies in the transition from deflagration to detonation in granular explosives-Ⅰ.Experimental arrangement and behavior of explosives which fail to exhibit detonation[J].Combust and Flame,1974,22(1):111 -118.

    [5] Bernecker R R,Price D.Studies in the transition from deflagration to detonation in granular explosives-Ⅱ.Transitional characteristics and mechanics observed in 91/9 RDX/wax[J].Combust and Flame,1974,22(1):119 -129.

    [6] Butler P B,Krier H.Analysis of deflagration to detonation transition in high-energy solid propellant[J].Combust.Flame,1986,63(1 -2):31 -48.

    [7] Bear M R,Nunziato J W.A two-phase mixture theory for the deflagration to detonation transition (DDT)in reactive granular materials[J].International Journal of Multiphase Flow,1986,12(6):861 -889.

    [8] Xu S J,Stewart D S.Deflagration-to-detonation transition in porous energetic materials:a comparative model study[J].Journal of Engineering Mechanics,1997,31(2 -3):143 -172.

    [9] Kapila A K,Menikoff R,Bdzil J B,et al.Two-phase modeling of deflagration-to-detonation transition in granular materials:reduced equations[J].Physics of Fluids,2001,13(10):3002 -3024.

    [10] Butler P B,Krier H.Modeling deflagration to detonation transition in granular explosives pentaerythritol tetranitrate[J].Journal of Applied Physics,2008,104(4):043519-(043519-14).

    [11] Powers J M,Stewart D S,Krier H.Analysis of steady compaction waves in porous materials[J].Journal of Applied Mechanics,1989,56(1):15 -24.

    [12] Sandusky H W,Liddiard T P.Dynamic compaction of porous beds[R].US:Naval Surface Warfare Center Silver Spring MD,1985.

    [13] Baer M R.Numerical studies of dynamic compaction of inert and energetic granular materials[J].Journal of Applied Mechanics,1988,55:36 -43.

    [14] Gonthier K A,Menikoff R,Son S E,et al.Molding compactioninduced energy dissipation of granular HMX[C]∥11th Symposium (International)on Detonation.1998:153 -161.

    [15] 孫錦山,朱建士,賈祥瑞.顆粒材料中致密波結(jié)構(gòu)研究[J].力學(xué)學(xué)報(bào),1999,31(4):423 -433.SUN Jin-shan,ZHU Jian-shi,JIA Xiang-rui.An analysis of compaction wave in granular material[J].Journal of Theoretical and Applied Mechanics,1999,31(4):423 -433.(in Chinese)

    [16] Low C A,Greenaway M W.Compaction process in granular beds composed of different particle sizes[J].Journal of Applied Physics,2005,98(12):123519-(123519 -12).

    [17] Low C A,Longbottom A W.Effect of particle distribution on the compaction behavior of granular beds[J].Physics of Fluids,2006,18(6):066101-(066101 -18).

    [18] Menikoff R,Kober E.Compaction waves in granular HMX[J].Shock Compression of Condensed Matter,1999,505:397 -400.

    [19] Bardenhagen S G,Brackbill J U.Dynamic stress bridging in granular material[J].Journal of Applied Physics,1998,83(11):5732 -5740.

    [20] Gonthier K A.Modeling and analysis of reactive compaction for granular energetic solids[J].Combustion Science and Technology,2003,175:1679 -1709.

    [21] Lovey A,Ben-Dor G,Sorek S,et al.Jump conditions across strong compaction waves in gas saturated rigid porous media[J].Journal of Shock Waves,1993,(3):105 -111.

    [22] Menikoff R,Sewell T D.Constituent properties of HMX needed for mesoscale simulations[J].Combustion Theory and Modeling,2002,(6):103 -125.

    [23] Powers J M.Two-phase viscous modeling of compaction of granular materials[J].Physics of Fluids,2004,16(8):2975 -2990.

    [24] Elban W L,Chiarito M A.Quasi-static compaction study of coarse HMX explosive[J].Power Technology,1986,46(2 -3):181 -193.

    [25] Moulard H,Kury J W,Declos A.The effect of RDX particle size on the shock sensitivity of cast PBX formulations[C]∥8th Symposium (International)on Detonation.1985:902 -913.

    [26] Moulard H.Particular Aspect of the explosive particle size effect on shock sensitivity of cast PBX formulations[C]∥9th Symposium (International)on Detonation.1989:18 -24.

    [27] Simpson R L,Helm F H,Crawford P C,et al.Particle size effects in the initiation of explosives containing reactive and nonreactive continuous phases[C]∥9th Symposium (International)on Detonation.1989:25 -38.

    [28] Gustavsen R L,Sheffield S A,Alcon R R.Low pressure shock initiation of porous HMX for two grain size distribution and two densities[J].American Institute of Physics,1996,370:851 -854.

    猜你喜歡
    構(gòu)形內(nèi)能活塞
    “內(nèi)能”“內(nèi)能的利用”綜合測(cè)試題
    一種活塞擠壓式室內(nèi)通風(fēng)設(shè)備的研發(fā)與應(yīng)用
    雙星跟飛立體成像的構(gòu)形保持控制
    通有構(gòu)形的特征多項(xiàng)式
    “內(nèi)能”“內(nèi)能的利用”綜合測(cè)試題
    “內(nèi)能和內(nèi)能的利用”易錯(cuò)點(diǎn)剖析
    “內(nèi)能”“內(nèi)能的利用”綜合測(cè)試題
    對(duì)一個(gè)幾何構(gòu)形的探究
    KS Kolbenschmidt公司的新型鋼活塞
    低噪聲活塞——降低活塞銷的拍擊噪聲
    免费观看的影片在线观看| 一个人观看的视频www高清免费观看 | 亚洲国产色片| 视频区欧美日本亚洲| 国产精品av视频在线免费观看| 五月伊人婷婷丁香| 国产精品久久久人人做人人爽| 精品国产美女av久久久久小说| 黄片小视频在线播放| 久久精品影院6| 国产精品av久久久久免费| 亚洲aⅴ乱码一区二区在线播放| 一个人观看的视频www高清免费观看 | 精品欧美国产一区二区三| 国产欧美日韩一区二区精品| 亚洲五月婷婷丁香| 69av精品久久久久久| 欧美日韩乱码在线| а√天堂www在线а√下载| 亚洲精品中文字幕一二三四区| 国产高清有码在线观看视频| 香蕉丝袜av| 欧美av亚洲av综合av国产av| 少妇熟女aⅴ在线视频| 精品99又大又爽又粗少妇毛片 | 精品免费久久久久久久清纯| 免费电影在线观看免费观看| 热99在线观看视频| 婷婷六月久久综合丁香| 久久亚洲精品不卡| 国内揄拍国产精品人妻在线| 男人舔女人下体高潮全视频| 999久久久国产精品视频| 黄色成人免费大全| 国产亚洲精品综合一区在线观看| 天堂网av新在线| 日韩欧美在线二视频| 97超级碰碰碰精品色视频在线观看| 男人和女人高潮做爰伦理| 免费搜索国产男女视频| 我要搜黄色片| 久久久久久久午夜电影| 精品乱码久久久久久99久播| av欧美777| 在线视频色国产色| www.自偷自拍.com| 国产99白浆流出| 免费av不卡在线播放| 国产精品99久久久久久久久| 中文字幕久久专区| 级片在线观看| 午夜成年电影在线免费观看| 91麻豆av在线| 男人的好看免费观看在线视频| 天天添夜夜摸| 亚洲国产欧美人成| 国产91精品成人一区二区三区| 国产真人三级小视频在线观看| 一区二区三区激情视频| 成人国产综合亚洲| 久久久久久久久中文| 久久人妻av系列| 欧美一级毛片孕妇| 最新在线观看一区二区三区| 国产美女午夜福利| 毛片女人毛片| 欧美xxxx黑人xx丫x性爽| 91av网站免费观看| 亚洲中文av在线| 久久99热这里只有精品18| 成人三级黄色视频| 国产三级中文精品| 身体一侧抽搐| 亚洲欧美精品综合一区二区三区| 9191精品国产免费久久| 男女那种视频在线观看| 欧美成人免费av一区二区三区| 热99re8久久精品国产| 欧美成狂野欧美在线观看| 日本 欧美在线| 舔av片在线| 波多野结衣高清作品| 五月伊人婷婷丁香| 午夜福利在线在线| 亚洲电影在线观看av| www日本在线高清视频| 小蜜桃在线观看免费完整版高清| 午夜久久久久精精品| 欧美成人免费av一区二区三区| 国产美女午夜福利| 久9热在线精品视频| 桃红色精品国产亚洲av| 日本免费a在线| 亚洲狠狠婷婷综合久久图片| 久久精品国产清高在天天线| 亚洲欧美日韩无卡精品| 99国产极品粉嫩在线观看| 日本黄大片高清| 国产视频一区二区在线看| 757午夜福利合集在线观看| 欧美日韩国产亚洲二区| 一级毛片高清免费大全| 精品一区二区三区av网在线观看| 日韩欧美国产在线观看| 成熟少妇高潮喷水视频| 午夜精品一区二区三区免费看| 亚洲精品乱码久久久v下载方式 | 很黄的视频免费| 久久九九热精品免费| 免费搜索国产男女视频| 搡老妇女老女人老熟妇| 久久久久久久午夜电影| 免费高清视频大片| 日韩欧美一区二区三区在线观看| 欧美av亚洲av综合av国产av| 国产精品久久久久久亚洲av鲁大| 亚洲电影在线观看av| 老司机在亚洲福利影院| 天天一区二区日本电影三级| 国产探花在线观看一区二区| 首页视频小说图片口味搜索| 亚洲aⅴ乱码一区二区在线播放| 搡老岳熟女国产| 女人高潮潮喷娇喘18禁视频| 一区二区三区高清视频在线| 亚洲电影在线观看av| 国产精品亚洲一级av第二区| 色老头精品视频在线观看| 老司机深夜福利视频在线观看| 波多野结衣巨乳人妻| 丁香六月欧美| 亚洲成av人片在线播放无| 一本综合久久免费| 国产精品99久久久久久久久| 国产伦在线观看视频一区| 欧美高清成人免费视频www| 精品久久久久久久人妻蜜臀av| 国产免费男女视频| 日本成人三级电影网站| 三级国产精品欧美在线观看 | 俺也久久电影网| 99久国产av精品| 观看免费一级毛片| 国产午夜精品论理片| 国产免费男女视频| 最好的美女福利视频网| 色综合站精品国产| 高清毛片免费观看视频网站| 丰满人妻熟妇乱又伦精品不卡| 久久久久国产精品人妻aⅴ院| 此物有八面人人有两片| 日韩三级视频一区二区三区| 男人和女人高潮做爰伦理| 色综合欧美亚洲国产小说| 欧美高清成人免费视频www| 一进一出抽搐gif免费好疼| 日韩欧美三级三区| 一本久久中文字幕| 久久伊人香网站| 岛国在线观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品女同一区二区软件 | 亚洲人成电影免费在线| 久久精品国产99精品国产亚洲性色| 亚洲在线观看片| 动漫黄色视频在线观看| 变态另类丝袜制服| 性色avwww在线观看| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 日韩欧美一区二区三区在线观看| 欧美激情在线99| 日本黄色片子视频| 久久久精品大字幕| 欧美日韩综合久久久久久 | 国产午夜福利久久久久久| 亚洲精品在线观看二区| 国产精品久久久人人做人人爽| 男女之事视频高清在线观看| 国产欧美日韩一区二区精品| 国产 一区 欧美 日韩| 久久国产精品人妻蜜桃| 免费看光身美女| 久久精品人妻少妇| 国产精品亚洲一级av第二区| 亚洲精品在线观看二区| 久久久久久国产a免费观看| 午夜精品在线福利| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| 无遮挡黄片免费观看| 色综合站精品国产| 亚洲人成电影免费在线| 人妻丰满熟妇av一区二区三区| 国产亚洲精品久久久久久毛片| 午夜福利在线在线| 丰满人妻一区二区三区视频av | 亚洲美女黄片视频| 欧美黄色片欧美黄色片| 神马国产精品三级电影在线观看| 国产成人系列免费观看| 国产欧美日韩一区二区三| 国产午夜精品久久久久久| 国产黄色小视频在线观看| 成人无遮挡网站| 精品一区二区三区视频在线 | 成人av在线播放网站| 国产麻豆成人av免费视频| 国产精品一区二区免费欧美| 99国产精品一区二区三区| 男人舔奶头视频| 国产高清视频在线观看网站| 母亲3免费完整高清在线观看| 欧美丝袜亚洲另类 | 啦啦啦观看免费观看视频高清| 成人永久免费在线观看视频| 在线观看免费视频日本深夜| 婷婷精品国产亚洲av在线| 99久久精品热视频| 一个人免费在线观看的高清视频| 成年人黄色毛片网站| 亚洲第一欧美日韩一区二区三区| 国产伦精品一区二区三区四那| 99视频精品全部免费 在线 | 好男人在线观看高清免费视频| 美女高潮的动态| 亚洲无线观看免费| 国产精品影院久久| 国产午夜精品久久久久久| 久久久久九九精品影院| 级片在线观看| 在线观看免费视频日本深夜| 久久精品aⅴ一区二区三区四区| or卡值多少钱| 国内精品久久久久久久电影| 99国产精品99久久久久| 国产三级中文精品| 国产精品av久久久久免费| 夜夜夜夜夜久久久久| 99国产精品一区二区三区| 国内揄拍国产精品人妻在线| 久99久视频精品免费| 国产精品98久久久久久宅男小说| 偷拍熟女少妇极品色| 久久热在线av| 久久婷婷人人爽人人干人人爱| 国内精品久久久久久久电影| 国产爱豆传媒在线观看| www.www免费av| 欧美乱码精品一区二区三区| 亚洲精品美女久久av网站| 色av中文字幕| 亚洲国产精品sss在线观看| 麻豆国产av国片精品| 老熟妇乱子伦视频在线观看| 国产一区二区三区在线臀色熟女| 国产激情欧美一区二区| 国产成年人精品一区二区| 免费大片18禁| 丁香六月欧美| www.999成人在线观看| 色综合欧美亚洲国产小说| 成人国产一区最新在线观看| 亚洲五月婷婷丁香| 久久国产精品影院| 99久久无色码亚洲精品果冻| 欧美精品啪啪一区二区三区| 91av网站免费观看| 国产亚洲精品av在线| 亚洲av免费在线观看| 少妇的丰满在线观看| 国产爱豆传媒在线观看| 午夜免费激情av| 岛国在线观看网站| 51午夜福利影视在线观看| 亚洲 欧美一区二区三区| 国产精品久久久久久精品电影| 国产单亲对白刺激| 99国产极品粉嫩在线观看| 成年人黄色毛片网站| 网址你懂的国产日韩在线| 淫秽高清视频在线观看| 一级毛片女人18水好多| 免费av毛片视频| 日本黄色片子视频| 国产男靠女视频免费网站| 亚洲精品久久国产高清桃花| 国产精品一区二区三区四区免费观看 | 亚洲精品美女久久av网站| 在线播放国产精品三级| 国产又色又爽无遮挡免费看| 成年女人毛片免费观看观看9| 999久久久国产精品视频| 午夜福利视频1000在线观看| 国产亚洲精品久久久久久毛片| 国内毛片毛片毛片毛片毛片| www.精华液| 老熟妇仑乱视频hdxx| 日本免费一区二区三区高清不卡| xxx96com| 999久久久国产精品视频| 国产精品香港三级国产av潘金莲| 免费看十八禁软件| 国产伦在线观看视频一区| 国产精品,欧美在线| 夜夜夜夜夜久久久久| 精品国产亚洲在线| 国产av不卡久久| 可以在线观看的亚洲视频| 国产精品99久久久久久久久| 丁香欧美五月| 丝袜人妻中文字幕| 在线国产一区二区在线| 国产又色又爽无遮挡免费看| 狂野欧美白嫩少妇大欣赏| 国产精品98久久久久久宅男小说| 亚洲狠狠婷婷综合久久图片| 色综合欧美亚洲国产小说| 国产精品一区二区三区四区免费观看 | 老鸭窝网址在线观看| 人妻丰满熟妇av一区二区三区| 免费看a级黄色片| 又大又爽又粗| h日本视频在线播放| 久久中文字幕一级| 国产免费男女视频| 欧美日韩乱码在线| 天天一区二区日本电影三级| 亚洲七黄色美女视频| 久久精品亚洲精品国产色婷小说| 在线观看美女被高潮喷水网站 | 热99re8久久精品国产| 黄色 视频免费看| 两性夫妻黄色片| 国产精品久久久久久久电影 | 国产av一区在线观看免费| 五月伊人婷婷丁香| 成熟少妇高潮喷水视频| 国内揄拍国产精品人妻在线| 久久伊人香网站| avwww免费| 狂野欧美激情性xxxx| 神马国产精品三级电影在线观看| 99久久精品国产亚洲精品| 成人国产一区最新在线观看| 亚洲人成电影免费在线| www.自偷自拍.com| e午夜精品久久久久久久| 亚洲精品在线观看二区| 最新在线观看一区二区三区| 美女大奶头视频| 首页视频小说图片口味搜索| 巨乳人妻的诱惑在线观看| 久久久国产成人精品二区| 最好的美女福利视频网| 狂野欧美激情性xxxx| 在线视频色国产色| 欧美日韩亚洲国产一区二区在线观看| 国产高清激情床上av| 亚洲国产中文字幕在线视频| 色噜噜av男人的天堂激情| 欧美极品一区二区三区四区| 男女做爰动态图高潮gif福利片| 特大巨黑吊av在线直播| 热99在线观看视频| 天堂影院成人在线观看| 88av欧美| 国产av麻豆久久久久久久| 麻豆国产97在线/欧美| 免费在线观看日本一区| 色综合站精品国产| 欧美日韩国产亚洲二区| 岛国在线免费视频观看| 白带黄色成豆腐渣| 日本 欧美在线| 一二三四在线观看免费中文在| 国产精品 国内视频| 国产精品98久久久久久宅男小说| 十八禁人妻一区二区| 国产精品久久久久久久电影 | 天堂网av新在线| 国产高清视频在线观看网站| 国产激情偷乱视频一区二区| 欧美性猛交╳xxx乱大交人| 久久久久久久久中文| 日韩欧美一区二区三区在线观看| 国产成人av教育| 国产熟女xx| 亚洲欧美日韩无卡精品| www.www免费av| 亚洲午夜理论影院| 免费电影在线观看免费观看| x7x7x7水蜜桃| 成人特级黄色片久久久久久久| 欧美性猛交╳xxx乱大交人| 性欧美人与动物交配| 特级一级黄色大片| 国产三级黄色录像| 亚洲一区高清亚洲精品| 又黄又粗又硬又大视频| 成人午夜高清在线视频| 亚洲专区国产一区二区| 久久午夜亚洲精品久久| 99久久久亚洲精品蜜臀av| 亚洲人成电影免费在线| 最好的美女福利视频网| 成人亚洲精品av一区二区| 国产主播在线观看一区二区| 男女那种视频在线观看| 女生性感内裤真人,穿戴方法视频| 成人精品一区二区免费| 在线观看美女被高潮喷水网站 | 婷婷丁香在线五月| 国产欧美日韩一区二区三| 国产在线精品亚洲第一网站| 琪琪午夜伦伦电影理论片6080| 国产1区2区3区精品| 嫩草影院精品99| 超碰成人久久| 欧美黑人巨大hd| 日本成人三级电影网站| 精品一区二区三区视频在线 | 噜噜噜噜噜久久久久久91| 国产伦在线观看视频一区| 欧美性猛交╳xxx乱大交人| 少妇的逼水好多| 久久中文字幕一级| 精品福利观看| 一本一本综合久久| 婷婷精品国产亚洲av在线| 欧美中文日本在线观看视频| 国产免费av片在线观看野外av| 两个人看的免费小视频| 久久久久久久午夜电影| 熟妇人妻久久中文字幕3abv| 夜夜爽天天搞| 色哟哟哟哟哟哟| 国产真实乱freesex| 午夜福利在线在线| 欧美日韩亚洲国产一区二区在线观看| 丰满人妻一区二区三区视频av | 欧美绝顶高潮抽搐喷水| 亚洲黑人精品在线| 全区人妻精品视频| 欧美日韩福利视频一区二区| 国产野战对白在线观看| 后天国语完整版免费观看| 国产免费男女视频| 国产午夜精品论理片| 国产精品一及| 99久久综合精品五月天人人| 老汉色∧v一级毛片| 91在线精品国自产拍蜜月 | 最新在线观看一区二区三区| ponron亚洲| 久久九九热精品免费| 白带黄色成豆腐渣| 日本a在线网址| 黑人欧美特级aaaaaa片| 国产精品av视频在线免费观看| www国产在线视频色| 亚洲精品一区av在线观看| 亚洲成人久久爱视频| 精品国产亚洲在线| 亚洲精华国产精华精| 91老司机精品| 99国产精品一区二区三区| 亚洲熟妇熟女久久| 国产精品电影一区二区三区| 久久久久久国产a免费观看| 看片在线看免费视频| 久久精品影院6| 脱女人内裤的视频| 国产av不卡久久| 91老司机精品| www.999成人在线观看| 国内精品久久久久精免费| 床上黄色一级片| 一级毛片女人18水好多| 黄片小视频在线播放| www.熟女人妻精品国产| 婷婷六月久久综合丁香| 美女扒开内裤让男人捅视频| 国产又色又爽无遮挡免费看| 国产精品永久免费网站| 国产精品免费一区二区三区在线| 亚洲天堂国产精品一区在线| 国产精品,欧美在线| 啦啦啦观看免费观看视频高清| 男女床上黄色一级片免费看| 中文字幕精品亚洲无线码一区| 国产三级在线视频| 中出人妻视频一区二区| 国产一区在线观看成人免费| 国产97色在线日韩免费| 久久欧美精品欧美久久欧美| 色视频www国产| 午夜福利18| 国产日本99.免费观看| 午夜福利在线在线| 操出白浆在线播放| 欧美三级亚洲精品| 中国美女看黄片| 黄色日韩在线| 91在线精品国自产拍蜜月 | 91老司机精品| 欧美午夜高清在线| 夜夜躁狠狠躁天天躁| 极品教师在线免费播放| 黑人巨大精品欧美一区二区mp4| 可以在线观看的亚洲视频| 国产一区在线观看成人免费| 中文字幕久久专区| 亚洲精品中文字幕一二三四区| 精品国产亚洲在线| 欧美中文日本在线观看视频| 88av欧美| 国产免费男女视频| 亚洲成av人片免费观看| 亚洲欧美日韩无卡精品| 精品一区二区三区av网在线观看| 观看美女的网站| 精品国产亚洲在线| 国产成人欧美在线观看| 亚洲18禁久久av| 搞女人的毛片| 午夜福利视频1000在线观看| 无遮挡黄片免费观看| 两个人看的免费小视频| 欧美日韩乱码在线| 色哟哟哟哟哟哟| 精品久久蜜臀av无| 久久精品影院6| 国内精品一区二区在线观看| 久久久色成人| 国产一区二区在线观看日韩 | 免费看光身美女| 精品欧美国产一区二区三| 少妇人妻一区二区三区视频| 欧美一级a爱片免费观看看| 性色av乱码一区二区三区2| 国产一级毛片七仙女欲春2| 久久精品国产99精品国产亚洲性色| 岛国视频午夜一区免费看| 特级一级黄色大片| 亚洲欧美日韩高清专用| 少妇丰满av| 啦啦啦韩国在线观看视频| 搡老妇女老女人老熟妇| 免费在线观看影片大全网站| 成熟少妇高潮喷水视频| 琪琪午夜伦伦电影理论片6080| h日本视频在线播放| 午夜激情福利司机影院| 国产精品一及| 亚洲专区字幕在线| 1024手机看黄色片| www日本黄色视频网| 可以在线观看毛片的网站| 一个人观看的视频www高清免费观看 | 亚洲18禁久久av| 欧美精品啪啪一区二区三区| 亚洲人成电影免费在线| 色噜噜av男人的天堂激情| 天堂av国产一区二区熟女人妻| 色精品久久人妻99蜜桃| ponron亚洲| 欧美丝袜亚洲另类 | 成人av在线播放网站| 婷婷精品国产亚洲av在线| 最新美女视频免费是黄的| 日本免费a在线| 亚洲成人久久性| 亚洲国产精品久久男人天堂| 舔av片在线| 成年女人看的毛片在线观看| 色av中文字幕| 淫秽高清视频在线观看| 成人特级黄色片久久久久久久| 精品久久蜜臀av无| 成年女人永久免费观看视频| 女同久久另类99精品国产91| 又爽又黄无遮挡网站| 五月玫瑰六月丁香| 少妇熟女aⅴ在线视频| 狂野欧美白嫩少妇大欣赏| 一进一出好大好爽视频| 搞女人的毛片| 日韩有码中文字幕| 久久天堂一区二区三区四区| 搞女人的毛片| 一进一出抽搐动态| 18禁国产床啪视频网站| 亚洲在线观看片| 夜夜爽天天搞| 人人妻人人澡欧美一区二区| 国内精品久久久久精免费| 母亲3免费完整高清在线观看| av中文乱码字幕在线| 国内精品久久久久精免费| 天天躁日日操中文字幕| 亚洲美女黄片视频| 国产精品电影一区二区三区| 国模一区二区三区四区视频 | 日本一本二区三区精品| 欧美zozozo另类| 99国产极品粉嫩在线观看| 久久久精品大字幕| 国产精品 国内视频| 村上凉子中文字幕在线| 黑人欧美特级aaaaaa片| 视频区欧美日本亚洲| 欧美激情久久久久久爽电影| 男插女下体视频免费在线播放| 国产精品影院久久|