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

    帶隔板裝藥爆轟波驅(qū)動(dòng)藥型罩?jǐn)?shù)值模擬研究*

    2019-12-28 06:23:58張建軍侯云輝劉燕寧
    關(guān)鍵詞:藥型罩隔板裝藥

    潘 建,張建軍,侯云輝,劉燕寧

    (西安現(xiàn)代控制技術(shù)研究所, 西安 710065)

    0 引言

    在常規(guī)武器系統(tǒng)的設(shè)計(jì)中,聚能裝藥戰(zhàn)斗部扮演著越來越重要的角色,提高炸藥能量利用率是戰(zhàn)斗部設(shè)計(jì)者追求的目標(biāo)之一。通常而言,采用環(huán)形或多點(diǎn)起爆、高低爆速的夾層裝藥和帶隔板裝藥等結(jié)構(gòu)可有效提高炸藥能量利用率。其中,帶隔板裝藥不僅結(jié)構(gòu)簡(jiǎn)單而且可以在減少裝藥量以及不改變裝藥長(zhǎng)徑比的條件下,提高其對(duì)藥型罩的壓垮做功能力。

    國(guó)內(nèi)外學(xué)者針對(duì)帶隔板裝藥驅(qū)動(dòng)藥型罩的作用過程開展了深入的研究。張先鋒等人[1]對(duì)帶隔板裝藥EFP成型進(jìn)行了相關(guān)研究,數(shù)值模擬和實(shí)驗(yàn)結(jié)果均表明,相對(duì)于無隔板裝藥,帶隔板裝藥所形成EFP具有較高的速度和較大的長(zhǎng)徑比。張彥國(guó)等[2]利用JWL和Lee-Tarver模型對(duì)帶大隔板裝藥爆轟波的傳播過程和EFP成型分別進(jìn)行了數(shù)值計(jì)算,計(jì)算得出帶大隔板裝藥結(jié)構(gòu)可以提高EFP的速度。門建兵等[3]分別利用3種敏感探測(cè)結(jié)構(gòu)(平板天線、波導(dǎo)管中饋和波導(dǎo)管邊饋)對(duì)EFP的成型進(jìn)行了數(shù)值模擬研究,結(jié)果表明中饋波導(dǎo)管結(jié)構(gòu)通過進(jìn)一步優(yōu)化可以得到滿意的EFP形狀。朱傳勝等[4]對(duì)帶隔板中空裝藥的EFP成型進(jìn)行了數(shù)值研究,最終優(yōu)化設(shè)計(jì)了一種弧錐結(jié)合罩結(jié)構(gòu),試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果吻合較好。

    文中的研究在前人工作基礎(chǔ)上,利用AUTODYN有限元計(jì)算程序中基于炸藥沖擊起爆的Lee-Tarver模型對(duì)三種帶隔板裝藥(SCJ、JPC和EFP)爆轟波驅(qū)動(dòng)藥型罩進(jìn)行數(shù)值模擬計(jì)算。

    1 帶隔板裝藥爆轟波傳播過程分析

    爆轟波繞過隔板后,在軸線位置處相互作用后發(fā)生正規(guī)反射或馬赫反射(超壓爆轟),反射壓力將超過CJ壓力,以pentolite(50/50 TNT/pentaerythritol tetranitrate)為例[5]。爆轟波在帶隔板裝藥的傳播過程中,隨入射角的增大,發(fā)生正規(guī)反射和馬赫反射后反射壓力與CJ壓力的比值如圖1所示。計(jì)算結(jié)果表明:當(dāng)發(fā)生正規(guī)反射時(shí),反射波壓力與CJ壓力的比值基本恒定,達(dá)到2.4;當(dāng)剛發(fā)生馬赫反射時(shí),馬赫桿上的壓力與CJ壓力比值最大為3.8,最小為2.7;隨著入射角再次增大,馬赫桿上的壓力值迅速降低,最終達(dá)到CJ值。根據(jù)文獻(xiàn)[6]知,帶隔板裝藥結(jié)構(gòu)爆轟波相互作用發(fā)生馬赫反射后,馬赫桿的高度值不斷變化。

    圖1 反射波壓力隨入射角變化的計(jì)算結(jié)果

    通過計(jì)算可以得出:在相同裝藥結(jié)構(gòu)中,適當(dāng)增加隔板直徑不僅可以有效改變爆轟波傳播方向,提高爆轟波對(duì)藥型罩的壓垮能力,而且可以有效增大馬赫桿的高度。為下一步帶隔板裝藥驅(qū)動(dòng)藥型罩的數(shù)值計(jì)算提供理論依據(jù)。

    2 數(shù)值模擬

    2.1 有限元計(jì)算方案

    根據(jù)對(duì)帶隔板裝藥爆轟波驅(qū)動(dòng)藥型罩過程的初步分析[7],在保證隔板外側(cè)爆轟波可以正常傳播的條件下,隔板直徑尺寸越大,調(diào)整爆轟波形狀的能力越強(qiáng),在裝藥結(jié)構(gòu)中爆轟波驅(qū)動(dòng)藥型罩的能力越強(qiáng)。

    帶隔板裝藥爆轟波驅(qū)動(dòng)藥型罩是一種多物質(zhì)相互作用的大變形運(yùn)動(dòng)過程,文中的研究采用多物質(zhì)耦合的Euler算法分析帶隔板裝藥爆轟波驅(qū)動(dòng)藥型罩過程。

    2.2 材料模型

    爆轟產(chǎn)物的JWL狀態(tài)方程可以精確描述爆轟波未到達(dá)隔板以及在隔板側(cè)面時(shí)爆轟波膨脹傳播過程。

    (1)

    式中:P為爆轟產(chǎn)物的壓力;V為爆轟產(chǎn)物相對(duì)比容;E為爆轟產(chǎn)物的比內(nèi)能;A、B、R1、R2、w為待擬合參數(shù)。

    藥型罩材料為高導(dǎo)無氧銅(OFHC),采用Johnson-Cook模型,狀態(tài)方程為Gruneisen方程。隔板材料為環(huán)氧樹脂(BORON EPXY),狀態(tài)方程為Puff狀態(tài)方程。藥型罩、隔板材料和炸藥參數(shù)主要來自文獻(xiàn)[8-9]及AUTODYN中內(nèi)置的參數(shù),其中TNT、B炸藥和8701炸藥的爆轟參數(shù)分別為:TNT炸藥的ρ=1.63 g/cm3,D=6.93 km/s,PCJ=21.0 GPa;B炸藥的ρ=1.717 g/cm3,D=7.980 km/s,PCJ=29.5 GPa;8701炸藥的ρ=1.70 g/cm3,D=8.425 km/s,PCJ=30.0 GPa。

    3 數(shù)值仿真結(jié)果與分析

    3.1 隔板結(jié)構(gòu)對(duì)高速射流(SCJ)數(shù)值仿真計(jì)算

    SCJ裝藥結(jié)構(gòu)參數(shù)來自于文獻(xiàn)[10]中的參數(shù)。采用中心起爆和Comp B炸藥,其狀態(tài)方程參數(shù)如2.2節(jié)所示,裝藥直徑為86.3 mm,長(zhǎng)度為152.4 mm。采用單錐形藥型罩,其錐角角度為42°,壁厚為2.69 mm,罩口直徑為86.3 mm。

    為了探討隔板結(jié)構(gòu)對(duì)高速射流的速度和長(zhǎng)度的影響,首先開展了有無隔板結(jié)構(gòu)壓垮藥型罩的數(shù)值模擬,高速射流裝藥結(jié)構(gòu)及射流成形結(jié)果如圖2所示。與不帶隔板裝藥相比,帶隔板裝藥結(jié)構(gòu)尺寸相同,圓臺(tái)形隔板最大直徑為60.42 mm,厚度為25 mm,采用酚醛樹脂材料。

    圖2 有無隔板裝藥結(jié)構(gòu)及射流成形結(jié)果

    從數(shù)值計(jì)算結(jié)果發(fā)現(xiàn):

    1)50 μs時(shí)帶隔板裝藥射流頭部速度為7 943.1 m/s,射流長(zhǎng)度為294.5 mm;不帶隔板裝藥射流頭部速度為6 788.2 m/s,射流長(zhǎng)度 253.5 mm,文獻(xiàn)[10]中試驗(yàn)結(jié)果為7 010 m/s,兩者誤差為3.3%;

    2)SCJ聚能裝藥中內(nèi)置隔板后,雖然裝藥量減少111.3 g,隔板體積占整體的10.7%,但是隔板的存在優(yōu)化了爆轟波形,最終使得射流速度提高14.5%,射流長(zhǎng)度提高13.9%,從而提高了射流的侵徹威力。

    在裝藥種類、長(zhǎng)徑比和藥型罩材料不變的情況下,影響藥型罩成型質(zhì)量因素主要有隔板直徑與裝藥直徑的比值Z、隔板厚度L、藥型罩錐角α、壁厚δ和隔板距起爆位置距離等參數(shù),文中研究的帶隔板裝藥主要針對(duì)小長(zhǎng)徑比的聚能裝藥(比如末敏彈等),該結(jié)構(gòu)中各個(gè)部件比較緊湊,隔板距起爆位置的距離通常保證裝藥可以正常起爆以及爆轟波可以正常傳播即可。因此通過改變Z、L、α和δ使其相互配合。根據(jù)正交優(yōu)化設(shè)計(jì)方法[11],將隔板直徑與裝藥直徑的比值Z、隔板厚度L、藥型罩錐角α和壁厚δ作為本次優(yōu)化設(shè)計(jì)的4個(gè)因素,各因素水平取值如表1所示。

    表1 正交優(yōu)化設(shè)計(jì)各因素水平取值表

    通過數(shù)值仿真計(jì)算得到高速射流的頭部速度和長(zhǎng)度,將其作為本次優(yōu)化設(shè)計(jì)的兩個(gè)評(píng)定指標(biāo)。用L16(45)正交表對(duì)4個(gè)因素所對(duì)應(yīng)的16個(gè)水平進(jìn)行合理安排,得到表2中的16個(gè)數(shù)值模擬計(jì)算方案。

    表2 正交優(yōu)化設(shè)計(jì)仿真計(jì)算方案

    按照表2所示的SCJ裝藥仿真計(jì)算方案,分別建立對(duì)應(yīng)的有限元仿真計(jì)算模型,依次計(jì)算得到聚能射流頭部速度、長(zhǎng)度和比動(dòng)能等模擬結(jié)果,如表3所示。

    通過計(jì)算得出:射流頭部速度隨著Z減小而減小,隨著α和δ減小而增大;射流長(zhǎng)度隨δ減小而增大;比動(dòng)能隨著Z減小而減小,隨δ減小而增大。

    根據(jù)極差分析法知:Z對(duì)頭部速度v和比動(dòng)能E都是第一重要目標(biāo),Z取0.9時(shí)v和E都存在極大值;L對(duì)v和E存在極大值時(shí)取19 mm;α和δ對(duì)v都是第一重要目標(biāo)影響較大,分別取42°和2.69 mm。所以,最終優(yōu)化的方案為編號(hào)3,此時(shí)射流頭部速度為9 112.5 m/s,射流長(zhǎng)度為310.5 mm。圖3為編號(hào)3在不同時(shí)刻所形成SCJ的成形外形,反映了SCJ的成形過程。

    表3 50 μs時(shí)高速射流頭部速度、長(zhǎng)度和比動(dòng)能

    圖3 SCJ成形過程數(shù)值模擬結(jié)果

    3.2 隔板結(jié)構(gòu)對(duì)桿式射流(JPC)數(shù)值仿真計(jì)算

    JPC裝藥中隔板結(jié)構(gòu)以倒錐形為主,其結(jié)構(gòu)參數(shù)來自于文獻(xiàn)[12]的參數(shù)。采用中心起爆和8701炸藥,裝藥直徑為110 mm,長(zhǎng)度為125 mm。采用壁厚和高度不變的偏心亞半球形藥型罩,選用紫銅材料,曲率半徑為90 mm,偏心距為39 mm,壁厚2.2 mm。

    為了探討隔板結(jié)構(gòu)對(duì)桿式射流的速度和長(zhǎng)度的影響,開展了有無隔板結(jié)構(gòu)壓垮藥型罩的數(shù)值模擬,桿式射流裝藥結(jié)構(gòu)及射流成形結(jié)果如圖4所示。與無隔板裝藥相比,帶隔板裝藥結(jié)構(gòu)尺寸相同,隔板直徑為100 mm,其圓柱段厚度和錐角分別為4.4 mm和110°,選用酚醛樹脂材料。

    從數(shù)值計(jì)算結(jié)果發(fā)現(xiàn):

    1)帶隔板裝藥射流頭部速度為8 693.6 m/s,60 μs時(shí)射流長(zhǎng)度為341 mm;不帶隔板裝藥射流頭部速度為7 324.3 m/s,60 μs時(shí)射流長(zhǎng)度為286.5 mm,計(jì)算結(jié)果與文獻(xiàn)[12]中結(jié)果吻合,從而驗(yàn)證了數(shù)值模型的正確性。

    圖4 有無隔板K裝藥結(jié)構(gòu)及射流成形結(jié)果

    2)JPC聚能裝藥中內(nèi)置隔板后,雖然裝藥量減少214.2 g,隔板體積占據(jù)整體的18.7%,但是隔板的存在優(yōu)化了爆轟波形,在裝藥中產(chǎn)生超壓爆轟,最終使得射流速度提高18.6%,從而可大大提高侵徹威力。

    在裝藥種類、長(zhǎng)徑比和藥型罩材料不變的情況下,影響射流質(zhì)量的主要因素是Z、偏心亞半球型藥型罩的圓心距軸線距離s、藥型罩曲率半徑r和壁厚δ等參數(shù)。根據(jù)正交優(yōu)化設(shè)計(jì)方法,將Z、s、r和δ作為本次優(yōu)化設(shè)計(jì)的4個(gè)因素,各因素水平取值如表4所示。

    表4 正交優(yōu)化設(shè)計(jì)各因素水平取值表

    通過數(shù)值仿真計(jì)算得到桿式射流的頭部速度、長(zhǎng)度和比動(dòng)能,將其作為本次優(yōu)化設(shè)計(jì)的3個(gè)評(píng)定指標(biāo)。用L16(45)正交表對(duì)4個(gè)因素所對(duì)應(yīng)的16個(gè)水平進(jìn)行合理安排,得到表5中的16個(gè)數(shù)值模擬計(jì)算方案。

    按照表5所示的K裝藥仿真計(jì)算方案,分別建立對(duì)應(yīng)的有限元仿真計(jì)算模型,依次計(jì)算得到桿式射流頭部速度、長(zhǎng)度和比動(dòng)能等模擬結(jié)果,如表6所示。

    通過計(jì)算得出:JPC裝藥頭部速度和長(zhǎng)度均隨Z、s減小而減小,隨r增大而減小,隨δ減小而增大;比動(dòng)能隨s減小而增大,其他因素在水平二和水平三存在極大值。

    根據(jù)極差分析法知:Z為頭部速度和長(zhǎng)度的第一重要指標(biāo),所以Z取為0.91;s和r為頭部速度和長(zhǎng)度的第二重要指標(biāo),所以s和r分別取39 mm和90 mm;δ對(duì)比動(dòng)能影響最大應(yīng)該取2.0 mm,但對(duì)于頭部速度和長(zhǎng)度來講應(yīng)該取1.9 mm。

    表5 正交優(yōu)化設(shè)計(jì)仿真計(jì)算方案

    表6 50 μs時(shí)桿式射流頭部速度、長(zhǎng)度和比動(dòng)能

    因此,δ應(yīng)該在1.9 mm和2.0 mm中取舍,而方案Z=0.91,s=39 mm,r=90 mm并不在表4所給的方案中,需要重新建立數(shù)值模型編號(hào)為17和18,計(jì)算結(jié)果如表7所示。

    表7 優(yōu)化方案仿真計(jì)算結(jié)果

    通過優(yōu)化方案計(jì)算結(jié)果得出,方案18的各項(xiàng)計(jì)算指標(biāo)均優(yōu)于表6中的16個(gè)方案,其頭部速度略低于方案17,但是比動(dòng)能卻高出11.6%,因此,該結(jié)構(gòu)戰(zhàn)斗部的優(yōu)化方案為:Z=0.91,s=39 mm,r=90 mm,δ=2.0 mm。圖5為優(yōu)化后方案在不同時(shí)刻所形成JPC的成形外形,反映了SCJ的成形過程。

    圖5 SCJ成形過程數(shù)值模擬結(jié)果

    3.3 隔板結(jié)構(gòu)對(duì)爆炸成型彈丸(EFP)數(shù)值仿真計(jì)算

    帶隔板EFP裝藥,其結(jié)構(gòu)參數(shù)來自于文獻(xiàn)[13]的參數(shù)。采用中心起爆和8701炸藥,裝藥直徑為110 mm,長(zhǎng)度為88 mm。采用弧錐結(jié)合型藥型罩,選用紫銅材料,其外半徑為75 mm,壁厚7.5 mm,罩口徑104 mm。

    為了探討隔板結(jié)構(gòu)對(duì)EFP的速度和長(zhǎng)度的影響,開展了有無隔板結(jié)構(gòu)壓垮藥型罩的數(shù)值模擬,裝藥結(jié)構(gòu)及EFP成形結(jié)果如圖6所示。與無隔板裝藥相比,帶隔板裝藥結(jié)構(gòu)尺寸相同,隔板直徑為66 mm,錐角90°,圓柱段長(zhǎng)度為3 mm,采用酚醛樹脂材料。

    圖6 有無隔板EFP裝藥結(jié)構(gòu)及射流成形結(jié)果

    從數(shù)值計(jì)算結(jié)果發(fā)現(xiàn):

    1)帶隔板裝藥EFP速度為1 336.9 m/s,長(zhǎng)徑比為3.7;不帶隔板裝藥EFP速度為1 229.0 m/s,長(zhǎng)徑比為1.4,計(jì)算結(jié)果與文獻(xiàn)[13]中結(jié)果吻合,從而驗(yàn)證了數(shù)值模型的正確性。

    2)EFP裝藥中內(nèi)置隔板后,雖然裝藥量減少了82.0 g,隔板體積占據(jù)整個(gè)體積的6.7%,但是隔板的存在優(yōu)化了爆轟波波形,在裝藥中產(chǎn)生超壓爆轟,最終使得EFP的速度提高了8.8%,其長(zhǎng)徑比提高164.3%,從而大大提高了飛行穩(wěn)定性。

    在裝藥種類、長(zhǎng)徑比和藥型罩材料不變的情況下,影響EFP質(zhì)量的主要因素是Z、隔板圓柱段厚度L、弧錐型藥型罩曲率半徑r和δ等參數(shù)。根據(jù)正交優(yōu)化設(shè)計(jì)方法,將Z、L、r、δ作為本次優(yōu)化設(shè)計(jì)的4個(gè)因素,各因素水平取值如表8所示。

    表8 正交優(yōu)化設(shè)計(jì)各因素水平取值表

    通過數(shù)值仿真計(jì)算得到EFP速度和長(zhǎng)度,將其作為本次優(yōu)化設(shè)計(jì)的兩個(gè)評(píng)定指標(biāo)。用L16(45)正交表對(duì)4個(gè)因素所對(duì)應(yīng)的16個(gè)水平進(jìn)行合理安排,得到表9中的16個(gè)數(shù)值模擬計(jì)算方案。

    表9 正交優(yōu)化設(shè)計(jì)仿真計(jì)算方案

    按照表9所示的EFP裝藥仿真計(jì)算方案,分別建立對(duì)應(yīng)的有限元仿真計(jì)算模型,依次計(jì)算得到EFP頭部速度、長(zhǎng)度和比動(dòng)能等模擬結(jié)果,如表10所示。

    通過計(jì)算得出:EFP的長(zhǎng)徑比、頭部速度和長(zhǎng)度均隨Z減小而減小,隨δ減小而增大;其他因素在變化中存在極值。

    根據(jù)極差分析法知:Z和δ為頭部速度和長(zhǎng)度的第一重要指標(biāo),所以Z和δ分別取為0.7和6.0 mm;L和r為頭部速度和長(zhǎng)度的第二重要指標(biāo),且存在極值,故取1 mm和75.0 mm。計(jì)算結(jié)果如表10中的方案4所示,但EFP在成形和飛行過程中由于速度過高而發(fā)生斷裂,所以最終優(yōu)化方案為編號(hào)10,Z=0.5,δ=6.5 mm,L=1 mm,r=75.0 mm。圖7為優(yōu)化后方案在不同時(shí)刻所形成EFP的成形外形,反映了EFP的成形過程。

    表10 帶隔板EFP速度、長(zhǎng)度、比動(dòng)能和長(zhǎng)徑比

    圖7 EFP成形過程數(shù)值模擬結(jié)果

    4 結(jié)論

    對(duì)帶隔板裝藥爆轟波驅(qū)動(dòng)藥型罩開展數(shù)值仿真模擬計(jì)算得出,在裝藥種類、長(zhǎng)徑比和藥型罩材料不變的情況下,通過匹配隔板和藥型罩結(jié)構(gòu)提高藥型罩成型質(zhì)量。主要結(jié)論如下:

    1)裝藥結(jié)構(gòu)中內(nèi)置隔板可改變爆轟波傳播路徑,爆轟波在軸線處相互作用產(chǎn)生超壓爆轟,可明顯提高炸藥對(duì)藥型罩的驅(qū)動(dòng)能力;

    2)通過隔板與藥型罩結(jié)構(gòu)之間的正交優(yōu)化得到兩者最佳匹配關(guān)系,優(yōu)化后SCJ裝藥的頭部速度提高25.1%、長(zhǎng)度提高11%,優(yōu)化后JPC裝藥的頭部速度提高18.6%,優(yōu)化后EFP裝藥的頭部速度提高8.8%,長(zhǎng)徑比提高164.3%。

    猜你喜歡
    藥型罩隔板裝藥
    《混合炸藥及裝藥??沸蜓?/a>
    銅基非晶合金雙層藥型罩射流形成及侵徹性能
    孔內(nèi)爆炸填塞效應(yīng)的數(shù)值模擬與分析
    基于有限元分析的隔板工藝支撐技術(shù)研究
    壓力容器隔板的一種設(shè)計(jì)方法
    某發(fā)射裝藥結(jié)構(gòu)改進(jìn)設(shè)計(jì)的新思路
    藥型罩材料對(duì)三層串聯(lián)EFP成型影響研究*
    橫隔板參數(shù)對(duì)裝配式箱梁寬橋受力的影響
    新型復(fù)合藥型罩設(shè)計(jì)研究
    雙頻帶隔板極化器
    国产亚洲av嫩草精品影院| 成人av在线播放网站| 又爽又黄无遮挡网站| 夜夜躁狠狠躁天天躁| 亚洲av五月六月丁香网| 午夜a级毛片| 一个人免费在线观看的高清视频| 真人一进一出gif抽搐免费| 身体一侧抽搐| АⅤ资源中文在线天堂| 久久精品91无色码中文字幕| 亚洲av电影不卡..在线观看| 免费看光身美女| 婷婷精品国产亚洲av在线| 十八禁网站免费在线| 国模一区二区三区四区视频 | 女同久久另类99精品国产91| 国产精品av视频在线免费观看| 国产精品女同一区二区软件 | 国产av在哪里看| 麻豆国产97在线/欧美| 青草久久国产| 免费观看精品视频网站| 日韩欧美 国产精品| 美女黄网站色视频| 1024手机看黄色片| 99久久精品一区二区三区| 国产午夜福利久久久久久| 亚洲激情在线av| 久久亚洲真实| 国产综合懂色| 亚洲国产高清在线一区二区三| 男女做爰动态图高潮gif福利片| 日日干狠狠操夜夜爽| 色视频www国产| 国产欧美日韩精品一区二区| 日韩欧美国产一区二区入口| 黄频高清免费视频| 免费看日本二区| 香蕉久久夜色| 久久国产乱子伦精品免费另类| 一个人免费在线观看的高清视频| www国产在线视频色| 国产激情偷乱视频一区二区| 午夜久久久久精精品| 日本五十路高清| 人人妻人人澡欧美一区二区| 国产乱人视频| 国产男靠女视频免费网站| 成人国产一区最新在线观看| 亚洲专区国产一区二区| 亚洲欧美日韩高清在线视频| 婷婷丁香在线五月| 久久精品国产99精品国产亚洲性色| 一边摸一边抽搐一进一小说| 亚洲aⅴ乱码一区二区在线播放| 淫秽高清视频在线观看| 亚洲av成人av| 国产精品一区二区三区四区久久| 日本 av在线| 久久久国产欧美日韩av| 神马国产精品三级电影在线观看| 变态另类丝袜制服| 美女 人体艺术 gogo| 一区二区三区国产精品乱码| 最好的美女福利视频网| 69av精品久久久久久| 色综合欧美亚洲国产小说| 欧美绝顶高潮抽搐喷水| av国产免费在线观看| 国产成人精品无人区| 蜜桃久久精品国产亚洲av| 婷婷精品国产亚洲av| 久久精品aⅴ一区二区三区四区| 中亚洲国语对白在线视频| 国产成人啪精品午夜网站| 一二三四在线观看免费中文在| 草草在线视频免费看| 久久亚洲真实| 国产又黄又爽又无遮挡在线| 两个人的视频大全免费| 精品国产乱码久久久久久男人| e午夜精品久久久久久久| 黑人操中国人逼视频| 日本与韩国留学比较| 精品国产乱码久久久久久男人| 老司机午夜福利在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 一级a爱片免费观看的视频| 男女之事视频高清在线观看| 国产高清videossex| 国产亚洲av高清不卡| 久久天堂一区二区三区四区| 在线观看免费午夜福利视频| 久久久久免费精品人妻一区二区| 亚洲片人在线观看| 久久久色成人| 亚洲国产日韩欧美精品在线观看 | 国产蜜桃级精品一区二区三区| 亚洲五月天丁香| 伦理电影免费视频| 我的老师免费观看完整版| 变态另类丝袜制服| 婷婷丁香在线五月| 黄色日韩在线| 成年女人永久免费观看视频| 两人在一起打扑克的视频| 亚洲专区字幕在线| 国产精品久久久久久精品电影| 这个男人来自地球电影免费观看| 国产激情欧美一区二区| 欧美+亚洲+日韩+国产| 国产精品自产拍在线观看55亚洲| 99热6这里只有精品| 久久久久免费精品人妻一区二区| 丰满的人妻完整版| 国产精华一区二区三区| 欧美高清成人免费视频www| 国产成人系列免费观看| 久久久久九九精品影院| 欧美日本视频| 国产蜜桃级精品一区二区三区| 国产亚洲精品久久久com| 啦啦啦观看免费观看视频高清| 日韩免费av在线播放| 噜噜噜噜噜久久久久久91| 国产97色在线日韩免费| 一个人看的www免费观看视频| 国内精品一区二区在线观看| 精品人妻1区二区| 一夜夜www| 校园春色视频在线观看| 午夜两性在线视频| 国产一区二区三区在线臀色熟女| 精品久久蜜臀av无| 动漫黄色视频在线观看| 噜噜噜噜噜久久久久久91| 精品久久蜜臀av无| 国产精品久久久久久久电影 | 别揉我奶头~嗯~啊~动态视频| 真人做人爱边吃奶动态| 亚洲国产精品sss在线观看| 母亲3免费完整高清在线观看| 国产综合懂色| 精品久久久久久,| 男女视频在线观看网站免费| 男女那种视频在线观看| 午夜影院日韩av| 91久久精品国产一区二区成人 | 99久久99久久久精品蜜桃| www日本黄色视频网| 欧美xxxx黑人xx丫x性爽| 亚洲在线观看片| 欧美+亚洲+日韩+国产| 国产高清有码在线观看视频| 最新美女视频免费是黄的| 日韩欧美 国产精品| 国产精品女同一区二区软件 | 国产av不卡久久| 免费在线观看亚洲国产| 在线观看免费视频日本深夜| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 国产精品综合久久久久久久免费| 欧美性猛交╳xxx乱大交人| 精品久久久久久成人av| 好男人电影高清在线观看| 欧美日韩乱码在线| 精品久久久久久久久久免费视频| 丁香欧美五月| 国产淫片久久久久久久久 | 亚洲一区高清亚洲精品| 午夜福利在线观看免费完整高清在 | 国产激情久久老熟女| 欧美三级亚洲精品| 51午夜福利影视在线观看| 亚洲欧洲精品一区二区精品久久久| 成人av一区二区三区在线看| cao死你这个sao货| 麻豆一二三区av精品| 18禁黄网站禁片免费观看直播| 亚洲精品久久国产高清桃花| 亚洲一区二区三区不卡视频| 亚洲18禁久久av| 亚洲一区高清亚洲精品| 国产高潮美女av| а√天堂www在线а√下载| 99久久成人亚洲精品观看| 欧美一区二区精品小视频在线| 欧美日韩乱码在线| 99久久精品一区二区三区| 亚洲国产看品久久| 99久久国产精品久久久| 欧美三级亚洲精品| 亚洲黑人精品在线| 9191精品国产免费久久| 美女高潮喷水抽搐中文字幕| 亚洲欧美精品综合久久99| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 少妇的逼水好多| 精品国产三级普通话版| 欧美丝袜亚洲另类 | 操出白浆在线播放| 少妇熟女aⅴ在线视频| 99精品欧美一区二区三区四区| 国产69精品久久久久777片 | 成年女人毛片免费观看观看9| 午夜两性在线视频| 亚洲国产精品成人综合色| 怎么达到女性高潮| 丰满的人妻完整版| 成人三级做爰电影| 97碰自拍视频| 青草久久国产| 欧美激情久久久久久爽电影| 一个人免费在线观看的高清视频| 亚洲欧洲精品一区二区精品久久久| 婷婷精品国产亚洲av在线| 18禁裸乳无遮挡免费网站照片| 91九色精品人成在线观看| 丰满人妻一区二区三区视频av | 欧美国产日韩亚洲一区| 久久国产乱子伦精品免费另类| 久久九九热精品免费| 天天一区二区日本电影三级| 嫩草影院精品99| 午夜成年电影在线免费观看| 色在线成人网| 三级毛片av免费| 婷婷精品国产亚洲av| 一个人免费在线观看的高清视频| 两个人看的免费小视频| 啦啦啦免费观看视频1| 成人av在线播放网站| 网址你懂的国产日韩在线| 国产亚洲欧美在线一区二区| 99久久成人亚洲精品观看| 人人妻人人澡欧美一区二区| 国产亚洲av嫩草精品影院| 亚洲最大成人中文| 视频区欧美日本亚洲| 日本在线视频免费播放| 国产成人系列免费观看| 99久久精品国产亚洲精品| 亚洲欧美精品综合一区二区三区| 国产高清有码在线观看视频| 欧美日韩国产亚洲二区| 成年女人永久免费观看视频| 欧美大码av| 丁香欧美五月| 国产亚洲欧美98| 精品国产超薄肉色丝袜足j| 最近最新中文字幕大全免费视频| 香蕉久久夜色| 中文字幕熟女人妻在线| 啦啦啦韩国在线观看视频| 国内揄拍国产精品人妻在线| 我的老师免费观看完整版| 免费高清视频大片| 天堂影院成人在线观看| 久久久国产欧美日韩av| 无人区码免费观看不卡| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 9191精品国产免费久久| 国产一区二区激情短视频| 成人18禁在线播放| 一a级毛片在线观看| 欧美乱妇无乱码| 亚洲在线自拍视频| 国产av麻豆久久久久久久| 日本免费a在线| 久久这里只有精品19| 中文字幕久久专区| 成人国产综合亚洲| 国产在线精品亚洲第一网站| 在线免费观看不下载黄p国产 | 欧美在线一区亚洲| or卡值多少钱| 18禁观看日本| 中出人妻视频一区二区| 国产成人av教育| 九色国产91popny在线| 国产人伦9x9x在线观看| 91老司机精品| 午夜久久久久精精品| 欧美另类亚洲清纯唯美| 亚洲av中文字字幕乱码综合| 天堂网av新在线| 制服丝袜大香蕉在线| 韩国av一区二区三区四区| 午夜日韩欧美国产| 夜夜看夜夜爽夜夜摸| 国产99白浆流出| 日本黄色视频三级网站网址| 好男人电影高清在线观看| 男女做爰动态图高潮gif福利片| 嫩草影院精品99| 久久精品aⅴ一区二区三区四区| 两性午夜刺激爽爽歪歪视频在线观看| 久久中文字幕一级| 夜夜爽天天搞| 欧美乱妇无乱码| 欧美日韩一级在线毛片| 久久久久精品国产欧美久久久| 国模一区二区三区四区视频 | 亚洲在线自拍视频| 午夜成年电影在线免费观看| 国产单亲对白刺激| 最新在线观看一区二区三区| 久久久久性生活片| 国产淫片久久久久久久久 | 成人三级做爰电影| 日韩欧美在线二视频| 精品一区二区三区av网在线观看| 国产免费男女视频| 亚洲国产欧美网| 亚洲片人在线观看| 国产精品av视频在线免费观看| 久久精品国产亚洲av香蕉五月| 亚洲av美国av| 国产精品免费一区二区三区在线| av欧美777| 丰满人妻一区二区三区视频av | 亚洲欧美日韩高清在线视频| 啦啦啦观看免费观看视频高清| 成年女人毛片免费观看观看9| 一级毛片女人18水好多| 757午夜福利合集在线观看| 夜夜爽天天搞| 1024手机看黄色片| 国产精品日韩av在线免费观看| 人人妻,人人澡人人爽秒播| 又黄又粗又硬又大视频| 亚洲av片天天在线观看| 女警被强在线播放| 波多野结衣巨乳人妻| 变态另类成人亚洲欧美熟女| 麻豆成人午夜福利视频| 好男人电影高清在线观看| 日日摸夜夜添夜夜添小说| 国产亚洲欧美98| 2021天堂中文幕一二区在线观| 亚洲色图 男人天堂 中文字幕| 岛国在线免费视频观看| 国产麻豆成人av免费视频| 久久久色成人| 可以在线观看的亚洲视频| 99国产精品一区二区蜜桃av| 欧美色视频一区免费| 国产高清有码在线观看视频| 久久亚洲真实| 国产精品香港三级国产av潘金莲| 国产精品久久久久久久电影 | 欧美一区二区国产精品久久精品| 视频区欧美日本亚洲| h日本视频在线播放| 免费在线观看影片大全网站| 69av精品久久久久久| 国产单亲对白刺激| 亚洲第一电影网av| 女同久久另类99精品国产91| av国产免费在线观看| 午夜福利欧美成人| 一本综合久久免费| 亚洲成人中文字幕在线播放| 亚洲国产精品999在线| 婷婷六月久久综合丁香| 草草在线视频免费看| 欧美性猛交╳xxx乱大交人| 操出白浆在线播放| 伦理电影免费视频| 动漫黄色视频在线观看| 国产精品一区二区三区四区免费观看 | 岛国视频午夜一区免费看| 日本五十路高清| 观看美女的网站| 在线播放国产精品三级| www日本在线高清视频| 亚洲自拍偷在线| 日本a在线网址| av片东京热男人的天堂| 一二三四在线观看免费中文在| 少妇丰满av| 国产精品自产拍在线观看55亚洲| 亚洲,欧美精品.| 三级毛片av免费| 免费看日本二区| 国产精品亚洲av一区麻豆| 久久久国产成人精品二区| 婷婷精品国产亚洲av在线| 亚洲电影在线观看av| 黄色成人免费大全| 午夜免费成人在线视频| 亚洲乱码一区二区免费版| 国产在线精品亚洲第一网站| 999久久久国产精品视频| 亚洲一区二区三区色噜噜| 香蕉丝袜av| 亚洲欧美日韩高清在线视频| 美女免费视频网站| 成人无遮挡网站| 国产伦在线观看视频一区| 19禁男女啪啪无遮挡网站| 校园春色视频在线观看| 国产高潮美女av| 性色av乱码一区二区三区2| 欧美午夜高清在线| 噜噜噜噜噜久久久久久91| 成人av在线播放网站| 国产精品永久免费网站| 欧美成人免费av一区二区三区| 欧美+亚洲+日韩+国产| 久久婷婷人人爽人人干人人爱| cao死你这个sao货| 欧美一级毛片孕妇| 国产视频一区二区在线看| cao死你这个sao货| 免费一级毛片在线播放高清视频| 欧美黄色片欧美黄色片| 国产单亲对白刺激| 此物有八面人人有两片| 国产黄色小视频在线观看| 两性夫妻黄色片| 日日摸夜夜添夜夜添小说| 久久久国产欧美日韩av| 亚洲国产欧美人成| aaaaa片日本免费| 一个人看的www免费观看视频| av在线天堂中文字幕| 蜜桃久久精品国产亚洲av| 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频| 女警被强在线播放| 国内精品久久久久精免费| 很黄的视频免费| 中文字幕最新亚洲高清| 成人特级黄色片久久久久久久| 国产高清视频在线观看网站| 桃红色精品国产亚洲av| 中文资源天堂在线| 女生性感内裤真人,穿戴方法视频| 欧美黄色淫秽网站| 精品人妻1区二区| 亚洲黑人精品在线| 欧美高清成人免费视频www| 女人被狂操c到高潮| 曰老女人黄片| АⅤ资源中文在线天堂| 成人三级做爰电影| 国产av在哪里看| 国产精华一区二区三区| 免费看光身美女| 亚洲精品美女久久av网站| 欧美国产日韩亚洲一区| 美女大奶头视频| 热99re8久久精品国产| 亚洲电影在线观看av| 香蕉丝袜av| 色综合婷婷激情| 免费高清视频大片| 天堂影院成人在线观看| 国产又色又爽无遮挡免费看| x7x7x7水蜜桃| 三级国产精品欧美在线观看 | 99热这里只有是精品50| 国产蜜桃级精品一区二区三区| 精品久久蜜臀av无| 精品人妻1区二区| 午夜激情欧美在线| 亚洲专区中文字幕在线| 1024手机看黄色片| 精品国产亚洲在线| 这个男人来自地球电影免费观看| 色综合婷婷激情| 听说在线观看完整版免费高清| 亚洲va日本ⅴa欧美va伊人久久| 又大又爽又粗| 好男人电影高清在线观看| 国产v大片淫在线免费观看| 1024香蕉在线观看| 亚洲av片天天在线观看| 国产精品 国内视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人一区二区三区免费视频网站| 搡老岳熟女国产| 九色国产91popny在线| 90打野战视频偷拍视频| 一个人看的www免费观看视频| 亚洲aⅴ乱码一区二区在线播放| 欧美乱色亚洲激情| or卡值多少钱| 免费一级毛片在线播放高清视频| 亚洲中文字幕一区二区三区有码在线看 | 成人18禁在线播放| 麻豆久久精品国产亚洲av| 日本 欧美在线| 亚洲成av人片在线播放无| 久久香蕉国产精品| 脱女人内裤的视频| 国产av不卡久久| 午夜成年电影在线免费观看| 精品国产三级普通话版| e午夜精品久久久久久久| 99久久久亚洲精品蜜臀av| 欧美日韩乱码在线| 丁香六月欧美| 亚洲国产欧洲综合997久久,| 精品久久蜜臀av无| 男女午夜视频在线观看| 俄罗斯特黄特色一大片| 最新中文字幕久久久久 | 免费看日本二区| 亚洲精品在线观看二区| 日韩有码中文字幕| 国产爱豆传媒在线观看| 国产精品美女特级片免费视频播放器 | 国产私拍福利视频在线观看| 一本精品99久久精品77| 宅男免费午夜| 亚洲在线自拍视频| 天堂√8在线中文| 99久久精品一区二区三区| 成年女人永久免费观看视频| 中文在线观看免费www的网站| 久久精品国产清高在天天线| 国产 一区 欧美 日韩| 他把我摸到了高潮在线观看| 日韩欧美在线二视频| 亚洲片人在线观看| 精品久久久久久成人av| 精品国产三级普通话版| 欧美国产日韩亚洲一区| 色吧在线观看| 观看免费一级毛片| 在线免费观看的www视频| 国产综合懂色| 国产av不卡久久| 小蜜桃在线观看免费完整版高清| 色综合站精品国产| 国产伦精品一区二区三区四那| 99在线人妻在线中文字幕| 99久久精品一区二区三区| 88av欧美| www.精华液| 日韩成人在线观看一区二区三区| 午夜精品久久久久久毛片777| 国内精品一区二区在线观看| 午夜精品一区二区三区免费看| 老汉色∧v一级毛片| 又大又爽又粗| 亚洲中文字幕日韩| a在线观看视频网站| 又黄又粗又硬又大视频| 一夜夜www| 天天一区二区日本电影三级| 亚洲欧美日韩卡通动漫| 网址你懂的国产日韩在线| 国产精品女同一区二区软件 | 窝窝影院91人妻| 深夜精品福利| 观看免费一级毛片| 日韩欧美三级三区| 最近最新中文字幕大全电影3| 99热只有精品国产| 国内精品久久久久久久电影| 国产99白浆流出| 亚洲欧美日韩高清专用| 两个人的视频大全免费| 在线免费观看不下载黄p国产 | 波多野结衣巨乳人妻| 久久久精品欧美日韩精品| 中文资源天堂在线| 午夜亚洲福利在线播放| 少妇熟女aⅴ在线视频| 老鸭窝网址在线观看| 国产精品99久久99久久久不卡| 免费看十八禁软件| 日韩国内少妇激情av| 欧美zozozo另类| 亚洲中文字幕日韩| 国产麻豆成人av免费视频| 好男人电影高清在线观看| 国产视频内射| 国产真人三级小视频在线观看| 首页视频小说图片口味搜索| 日韩欧美 国产精品| x7x7x7水蜜桃| 国内精品美女久久久久久| 亚洲色图 男人天堂 中文字幕| 一二三四在线观看免费中文在| 精品欧美国产一区二区三| 亚洲熟女毛片儿| 嫁个100分男人电影在线观看| 999精品在线视频| 中文字幕精品亚洲无线码一区| 亚洲成人精品中文字幕电影| 一卡2卡三卡四卡精品乱码亚洲| www.999成人在线观看| 岛国在线观看网站| 很黄的视频免费| 十八禁网站免费在线| 国产精品1区2区在线观看.| 欧美色欧美亚洲另类二区| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品合色在线| 亚洲人成伊人成综合网2020| 天堂动漫精品| 亚洲av五月六月丁香网| 亚洲成a人片在线一区二区| 亚洲精品美女久久av网站| 在线a可以看的网站| 老司机午夜福利在线观看视频| 国产一区二区在线av高清观看|