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

    基于SPH的橢圓雙極線性聚能藥包控制爆破數(shù)值模擬研究

    2022-07-20 01:13:08徐世祥黃宗輝
    煤炭科學(xué)技術(shù) 2022年5期
    關(guān)鍵詞:藥型罩藥包雙極

    吳 波,韋 漢,徐世祥,許 杰,黃宗輝

    (1.廣西大學(xué) 土木建筑工程學(xué)院,廣西 南寧 530004;2.教育部工程防災(zāi)與結(jié)構(gòu)安全重點(diǎn)實(shí)驗(yàn)室,廣西 南寧 530004;3.廣西防災(zāi)減災(zāi)與工程安全重點(diǎn)實(shí)驗(yàn)室,廣西 南寧 530004)

    0 引 言

    聚能裝藥廣泛應(yīng)用于軍事領(lǐng)域,比如爆炸成型彈丸、破甲彈設(shè)計(jì)等,此外聚能裝藥也被廣泛應(yīng)用于工業(yè)生產(chǎn),尤其是在實(shí)現(xiàn)巖石的定向斷裂控制爆破技術(shù)中。線性聚能裝藥是通過擠壓聚能藥型罩形成聚能射流的,從而達(dá)到定向斷裂控制爆破的效果。近年來,人們對線性聚能裝藥爆破技術(shù)進(jìn)行了廣泛的研究,由于爆破的瞬時性,實(shí)驗(yàn)和理論難以捕捉爆炸的整個過程,而數(shù)值模擬可以再現(xiàn)爆炸的過程,使之成為續(xù)實(shí)驗(yàn)和理論之后的一種研究手段。MOLINARI[1]采用有限元來研究聚能裝藥中銅射流的形成、破碎和侵徹過程。王成等[2]通過ALE研究不同起爆方式對線性聚能射流形成的影響進(jìn)行分析。張先鋒等[3]也通過ALE對3種不同典型線性裝藥的聚能射流形成及侵徹過程進(jìn)行了數(shù)值模擬。李必紅[4]從理論、實(shí)驗(yàn)與數(shù)值模擬研究橢圓雙極線性聚能藥包的爆破機(jī)理,提出了有效裝藥邊界函數(shù)并通過理論推導(dǎo)確定。

    以上實(shí)現(xiàn)的數(shù)值分析都是基于網(wǎng)格的拉格朗日方法、歐拉方法或者兩者的結(jié)合,然而實(shí)現(xiàn)的拉格朗日方法由于大變形的存在往往使得計(jì)算終止,而歐拉方法或ALE方法對物質(zhì)邊界的捕捉存在一定的困難。由于無網(wǎng)格法具有處理大變形和破壞問題的特點(diǎn),使得其成為處理極端變形條件下有限元無法比擬的優(yōu)勢,尤其是無網(wǎng)格SPH方法在線性聚能裝藥中的應(yīng)用。SPH方法最早于1977年由LUCY等[5]提出并將其天體物理學(xué),隨后廣泛應(yīng)用于各個領(lǐng)域的研究,SWEGLE等[6]首次將其應(yīng)用于水下爆炸數(shù)值研究,實(shí)現(xiàn)SPH在爆炸領(lǐng)域的應(yīng)用,LIU等[7-8]等基于SPH對無藥型罩聚能裝藥爆轟過程及聚能射流形成過程進(jìn)行研究。YANG等[9]基于二維模型對線性聚能射流進(jìn)行SPH模擬研究,并對比分析不同數(shù)值分析方法對聚能射流的影響。李磊等[10-11]通過SPH實(shí)現(xiàn)聚能裝藥的三維數(shù)值模擬,并將SPH方法與有限元耦合實(shí)現(xiàn)聚能射流對鋼板的侵徹效果。楊建輝[12]通過SNSYS/LS-DYNA對切縫聚能管參數(shù)進(jìn)行優(yōu)化,分析不同參數(shù)下對巖石應(yīng)力狀態(tài)的影響。陳帥志等[13]基于LS-DYNA分析切縫藥包與普通藥包爆炸在巖石的傳播特性。岳中文等[14]結(jié)合數(shù)值模擬及試驗(yàn)研究炮孔間距對切縫藥包爆破引起的裂紋擴(kuò)展的影響。張盛等[15]將橢圓雙極線性聚能藥包在巷道中的應(yīng)用,并結(jié)合現(xiàn)場試驗(yàn)對參數(shù)進(jìn)行優(yōu)化研究。

    綜上所述,目前對聚能裝藥形成數(shù)值分析主要基于有限元分析,而有限元模型無法動態(tài)捕捉爆生氣體跟聚能罩的動態(tài)演化規(guī)律,或?qū)Ь勰懿鄣木匦嗡幇蚯锌p藥包進(jìn)行聚能分析,而基于無網(wǎng)格實(shí)現(xiàn)橢圓雙極線性聚能裝藥結(jié)構(gòu)的研究較少。鑒于SPH方法在處理爆炸這種極端變形等復(fù)雜條件下具有有限元無法比擬的優(yōu)勢。筆者基于SPH方法研究聚能裝藥的形成,并與已有試驗(yàn)及模擬做對比,驗(yàn)證本文實(shí)現(xiàn)的SPH方法的有效性;然后基于SPH研究橢圓雙極線性聚能藥包爆破機(jī)理并分析外殼及其厚度、藥型罩及其厚度對橢圓雙極線性聚能藥包實(shí)現(xiàn)聚能射流及其對巖石定向斷裂控制爆破技術(shù)的影響。

    1 SPH控制方程理論

    SPH方法的核心思想就是通過構(gòu)造一個核函數(shù)對粒子的核估計(jì)來計(jì)算梯度的相關(guān)量,不需要求解偏微分方程,而是將偏微分控制方程轉(zhuǎn)化為積分形式,偏微分方程的轉(zhuǎn)化包括函數(shù)核近似和粒子近似兩步,函數(shù)核〈f(x)〉[11]近似如下:

    (1)

    式中,W(x-x′,h)為光滑函數(shù)或核函數(shù);Ω為積分域;x為空間坐標(biāo);x′為x緊支域內(nèi)某一點(diǎn)的坐標(biāo);h為光滑長度。

    對于連續(xù)形式的核近似式,函數(shù)f(x)的粒子近似可對粒子i緊支域內(nèi)所有粒子j進(jìn)行加權(quán)求和[16],即

    (2)

    式中:mj,ρj分別為粒子i的臨近粒子j的質(zhì)量和密度。

    光滑函數(shù)需要滿足歸一化、δ函數(shù)性質(zhì)和緊支性3個基本條件,常見的光滑函數(shù)有3次B樣條函數(shù)、5次B樣條函數(shù)、高斯核函數(shù)、指數(shù)核函數(shù)等,本文采用的是最常用的3次B樣條核函數(shù)[11],形式如下:

    (3)

    在粒子近似過程中,如果采用固定光滑長度,粒子集中時將影響計(jì)算效率,同時對于具有強(qiáng)度的材料可能造成張力不穩(wěn)定的現(xiàn)象,粒子散射時會出現(xiàn)粒子緊支域內(nèi)粒子數(shù)減少造成數(shù)值不穩(wěn)定,因此本文采用變光滑長度法,其光滑長度的時間積分[17]格式如下:

    (4)

    式中:d為空間維數(shù),本文取d=2;v為粒子速度;h(t)為粒子光滑長度;div為輔助函數(shù)。

    2 模型驗(yàn)證

    2.1 試驗(yàn)描述問題

    通過選取與文獻(xiàn)[18]相同的幾何參數(shù)模型,模型參數(shù)分別為藥型罩頂角W=80,厚度D=1.35 mm,裝藥長度L=25.4 mm,裝藥高度H=28.45 mm,模型如圖1所示。通過LS-DYNA刪除Lagrange單元信息再添加SPH相關(guān)屬性,實(shí)現(xiàn)的SPH二維分析模型如圖2所示。試驗(yàn)?zāi)P蛯?shí)現(xiàn)的SPH粒子間距為0.008~0.010 cm,炸藥和藥型罩的最小最大光滑長度分別為初始光滑長度的0.2倍和2倍,由于模型對稱,計(jì)算時取1/2模型進(jìn)行計(jì)算。

    圖1 數(shù)值模型Fig.1 Numerical model

    圖2 SPH離散模型Fig.2 SPH discrete model

    炸藥采用B炸藥進(jìn)行分析模擬,并采用MAT_HIGH_EXPLOSIVE_BURN模型描述材料本構(gòu),狀態(tài)方程為JWL:

    (5)

    式中:p為爆轟壓力;E為比內(nèi)能;V為相對比容;A、B、R1、R2、ω為狀態(tài)方程系數(shù),炸藥材料參數(shù)及狀態(tài)方程參數(shù)見表1[18]。

    表1 B炸藥材料參數(shù)Table 1 B explosive material parameters

    金屬藥型罩同樣采用文獻(xiàn)[18]的材料1006號鋼,本構(gòu)模型為MAT_JOHNSON_COOK模型,狀態(tài)方程為EOS_GRUNEISEN:

    (6)

    表2 1006號鋼材料參數(shù)Table 2 Material parameters of No.1006 steel

    2.2 模擬結(jié)果分析

    圖3為聚能射流形成過程中不同時刻射流形態(tài)及速度分布。由圖3可知,基于LS-DYNA實(shí)現(xiàn)的SPH線性聚能射流與文獻(xiàn)[18]基于高斯光滑核函數(shù)的自編程序?qū)崿F(xiàn)的聚能射流在不同時刻基本一致,且都在t=4 μs 時爆轟基本完成,射流初步形成,如圖3a所示。

    圖3 聚能射流形成不同時刻速度云圖Fig.3 Cloud chart of velocity at different time of jet formation

    圖4為t=6.5μs時聚能射流的速度云圖,此刻射流頭部平均最大速度達(dá)到3 507 m/s,與文獻(xiàn)[19]電磁線圈測得射流頭部最大平均速度值3 520 m/s基本一致。文獻(xiàn)[18]實(shí)現(xiàn)的SPH方法獲取的射流頭部峰值速度為3 294 m/s,文獻(xiàn)[19]通過X光技術(shù)測得的射流頭部峰值速度為3 300~3 500 m/s,本文實(shí)現(xiàn)的SPH聚能射流形成穩(wěn)定之后射流頭部平均速度為3 180 m/s,說明本文實(shí)現(xiàn)的SPH方法和文獻(xiàn)[18]及文獻(xiàn)[19]基本一致,表明本文基于DYNA實(shí)現(xiàn)的SPH方法的有效性。

    圖4 t=6.5 μs時刻聚能射流速度云圖Fig.4 Cloud chart of jet velocity at t=6.5 μs

    3 橢圓雙極線性聚能藥包爆破數(shù)值分析

    基于前文實(shí)現(xiàn)的SPH方法,可知SPH可以處理極端變形復(fù)雜條件下的破壞問題,為此對橢圓雙極線性聚能藥包爆破機(jī)理及其影響因素(外殼及其厚度、藥型罩及其厚度)進(jìn)行分析。

    3.1 橢圓雙極線性聚能藥包形成

    3.1.1 模型的建立

    圖5為橢圓雙極線性聚能藥包四分之一模型,幾何參數(shù)如圖5a所示,長短軸比為11∶15,離散模型如圖5b所示,為考慮薄外殼與藥型罩具有足夠的臨近粒子數(shù),減小外殼與藥型罩的粒子間距,考慮模型尺寸較小,將炸藥粒子間距設(shè)置為外殼及藥型罩相同的粒子間距0.008~0.010 cm。模型材料參數(shù)取值與上節(jié)材料參數(shù)一致,外殼選取與藥型罩相同的材料及本構(gòu)方程和狀態(tài)方程。整個模型采用變光滑長度的計(jì)算方法,炸藥、外殼和藥型罩的最小最大光滑長度分別為初始光滑長度的0.2倍和2倍,計(jì)算時采用對稱邊界處理。

    圖5 四分之一聚能藥包模型Fig.5 Quarter model of shaped-charge

    3.1.2 聚能藥包作用機(jī)理分析

    圖6為聚能射流形成不同時刻速度分布云圖。t=0.8 μs時爆轟波先傳到藥型罩錐角頂端,開始擠壓藥型罩形成射流,隨后到達(dá)90°非聚能方向的外殼,開始沖擊作用于90°非聚能方向的外殼,使之開始發(fā)生壓縮變形與向外運(yùn)動。t=4.0 μs時射流初步形成,射流頭部平均最大速度達(dá)到3 224 m/s,外殼不斷發(fā)生塑性拉伸及向外運(yùn)動,同時外殼與藥型罩連接部位受拉力作用不斷拉伸縮頸。在t=14.0 μs時連接部位首先發(fā)生斷裂,爆生氣體從該位置溢出,此時射流不斷被拉伸,并與存在速度梯度的杵體部分逐步分離,由于射流不再受杵體、爆生氣體和外殼的影響,射流頭部平均速度基本保持不變。但最大速度并未發(fā)生在射流頭部,而是在藥型罩與外殼連接處的炸藥粒子,由于射流形成穩(wěn)定后不再受其他因素的影響,而接觸部位炸藥粒子在脫離藥型罩與外殼的約束作用之后受到內(nèi)部相鄰炸藥粒子的直接加速作用。t=20.0 μs時,脫離外殼與藥型罩的炸藥粒子不斷加速追趕射流頭部,此時炸藥粒子最大速度達(dá)到3 772 m/s,而射流頭部平均最大速度為2 985 m/s,同時在90°非聚能方向的外殼不斷發(fā)生縮頸甚至斷裂。

    圖6 聚能射流形成不同時刻速度分布云圖Fig.6 Cloud chart of velocity distribution at different time of jet formation

    圖7為不同位置粒子的時間歷程曲線,A′為射流頭部粒子編號為302614,B′為藥型罩與外殼接觸部位內(nèi)側(cè)的炸藥粒子編號為292551,C′為90°非聚能方向外殼粒子編號為302831,D′為90°非聚能方向外殼內(nèi)側(cè)炸藥粒子編號為290904。由圖可知,聚能射流在t=4 μs時射流基本形成,與圖6b表述相吻合。t=14.0 μs時爆生氣體溢出藥型罩與外殼連接位置,最終在t=15.0 μs后粒子以3 440 m/s的速度穩(wěn)定向前運(yùn)動。90°非聚能方向外殼粒子與外殼內(nèi)部炸藥粒子的速度增加的總體趨勢一致,但由于炸藥粒子受到爆轟波的往復(fù)作用,速度波動幅度相對外殼粒子要大。此外,射流頭部速度發(fā)展較快,其他粒子速度存在明顯的時間滯后現(xiàn)象。

    圖7 不同位置處粒子速度隨時間變化曲線Fig.7 Curve of particle velocity with time at different positions

    3.2 聚能藥包參數(shù)優(yōu)化

    3.2.1 問題描述

    不同幾何參數(shù)的橢圓雙極線性聚能藥包四分之一模型如圖8所示,在保證長短軸比、聚能槽角度和裝藥量一致的情況下,各個模型采用相同的材料參數(shù)和粒子間距(同上節(jié))。粒子的最大、最小光滑長度分別為0.2倍和2倍初始光滑長度,為保證初始時刻每個粒子的支持域內(nèi)具有相同的質(zhì)量,每個粒子初始臨近粒子數(shù)設(shè)置為600。

    圖8 不同幾何參數(shù)的聚能藥包四分之一模型Fig.8 Quarter models of shaped-charge cartridge with different geometric parameters

    3.2.2 數(shù)值模擬結(jié)果及分析

    圖9為不同幾何參數(shù)聚能藥包在t=14.0 μs時速度云圖,由圖可知不同幾何參數(shù)的聚能藥包實(shí)現(xiàn)的定向效果不同。圖9a為2 mm厚的藥型罩1 mm厚外殼的聚能藥包速度分布圖,由圖可知除射流頭部速度比較大之外,外殼向外運(yùn)動的速度也較大,整體速度大于杵體速度,主要是由于外殼比較薄,爆生氣體對外殼做的功大部分轉(zhuǎn)化為外殼動能。對比圖9b無外殼情況可知,圖9a形成的射流速度較大且藥型罩大部分的粒子形成射流,而圖9b藥型罩大部分粒子用于形成杵體,表明外殼對聚能射流的形成產(chǎn)生較大影響。圖9c為1 mm厚藥型罩和2 mm厚外殼的聚能藥包速度分布圖,對比圖6c的2 mm厚藥型罩情況可知,兩者在同一時刻形成的射流形態(tài)相似,但前者形成的射流頭部速度更大,射流頭部平均最大速度達(dá)到3 670 m/s,更有利于巖石的侵徹作用。圖9d由于無藥型罩而無法形成射流,圖中只給出爆生氣體從外殼切縫處溢出形態(tài),該形態(tài)與申濤等[20]基于ALE算法模擬切縫藥包爆破相似,由圖可見粒子散射比較明顯,實(shí)現(xiàn)的定向效果不顯著。圖9e為1 mm厚藥型罩和1 mm厚外殼的聚能藥包速度分布圖,對比圖6c可知同樣前者形成的射流更明顯,速度更大,表明藥型罩和外殼厚度同時減小時有利于聚能射流的形成。圖9f為無外殼無藥型罩藥包的爆生氣體速度分布圖,由圖可見,并未形成明顯射流,可見外殼與藥型罩對聚能射流的形成起決定作用。

    圖10為不同幾何參數(shù)聚能藥包射流頭部速度隨時間變化曲線,通過比較圖5a、圖8a和圖8b所對應(yīng)的曲線可知,隨著外殼厚度不斷減小,最終射流頭部速度不斷減小,但三者的速度變化趨勢基本一致,再次表明外殼對橢圓雙極線性聚能藥包射流頭部的形成有一定的影響,外殼越厚更有利于穩(wěn)定爆轟作用,從而有利于射流的形成,射流速度越大;比較圖5a、圖8c所對應(yīng)的曲線,可知當(dāng)藥型罩為1 mm,外殼2 mm時形成射流頭部速度最大,最大值為3 700 m/s,而圖8d由于無藥型罩而無法形成聚能射流,對應(yīng)圖10的曲線只不過是聚能方向外殼邊緣處粒子的速度時間歷程曲線,表明藥型罩對聚能射流的形成起決定性作用;對比圖5a和圖8e所對應(yīng)的曲線圖可知,隨著藥型罩和外殼厚度一致減小,聚能射流頭部速度增大,但對比圖9f可知當(dāng)兩者減小為0時不形成射流。對比圖8d和圖8f所對應(yīng)的曲線,由于兩者均沒有藥型罩,只能借助炸藥粒子分析,由圖可見,兩者速度曲線重合,即最終穩(wěn)定速度完全一致,表明當(dāng)無藥型罩時并沒有明顯的聚能效應(yīng)。

    圖9 不同幾何參數(shù)聚能藥包在t=14.0 μs時速度云圖Fig.9 Velocity nephogram of shaped charge with different geometric parameters at t=14.0 μs

    圖10 不同幾何參數(shù)聚能藥包射流頭部速度時間歷程曲線Fig.10 Time history curve of jet head velocity of shaped-charge with different geometric parameters

    4 結(jié) 論

    1)橢圓雙極線性聚能藥包通過壓垮藥型罩形成聚能射流,聚能射流形成的同時爆生氣體不斷撞擊外殼向非聚能方向運(yùn)動,變形到達(dá)一定程度之后外殼與藥型罩連接部位斷開,爆生氣體從斷開部位溢出,溢出的爆生氣體速度大于射流頭部速度。

    2)隨著外殼厚度不斷增大,射流頭部速度也在增大,且不同外殼厚度的聚能藥包射流頭部速度隨時間變化趨勢趨于一致,表明外殼對橢圓雙極線性聚能藥包射流頭部速度有一定影響,外殼越厚越利于穩(wěn)定爆轟,從而利于射流的形成,射流速度越大。

    3)當(dāng)外殼厚度一定時,藥型罩厚度不斷減小,爆生氣體用于壓垮藥型罩做的功越少,轉(zhuǎn)換為聚能射流動能越多,聚能射流頭部速度越大,但當(dāng)藥型罩厚度為0時未能形成明顯的聚能效應(yīng)。

    猜你喜歡
    藥型罩藥包雙極
    基于雙極化解耦超表面的線轉(zhuǎn)圓極化反射陣列天線設(shè)計(jì)
    無線電工程(2024年9期)2024-10-24 00:00:00
    雙極直覺模糊超圖*
    銅基非晶合金雙層藥型罩射流形成及侵徹性能
    柱狀藥包爆破漏斗效率函數(shù)及其等效爆破作用
    更正
    含能材料(2020年5期)2020-01-14 00:19:10
    藥型罩材料對三層串聯(lián)EFP成型影響研究*
    水下爆炸中水面效應(yīng)以及藥包形狀對沖擊波的影響
    中國測試(2018年10期)2018-11-17 01:58:26
    新型復(fù)合藥型罩設(shè)計(jì)研究
    藥型罩切分方式對射流形成影響的數(shù)值模擬
    鼻內(nèi)鏡下雙極電凝治療嚴(yán)重鼻出血的療效
    操美女的视频在线观看| 亚洲 国产 在线| 人人澡人人妻人| 久久香蕉激情| 黄频高清免费视频| 美女高潮到喷水免费观看| 成人三级黄色视频| 九色亚洲精品在线播放| 国产高清国产精品国产三级| bbb黄色大片| 国产熟女午夜一区二区三区| 99国产综合亚洲精品| av欧美777| ponron亚洲| 国产精品98久久久久久宅男小说| 真人做人爱边吃奶动态| 国产高清videossex| 日本wwww免费看| 一区二区三区精品91| 人成视频在线观看免费观看| 午夜福利免费观看在线| 国产一区二区激情短视频| 宅男免费午夜| 日本五十路高清| 老熟妇仑乱视频hdxx| 国产精品1区2区在线观看.| 午夜91福利影院| 在线播放国产精品三级| 亚洲精品中文字幕一二三四区| 日韩精品免费视频一区二区三区| 波多野结衣高清无吗| 欧美激情久久久久久爽电影 | 午夜福利免费观看在线| 国产蜜桃级精品一区二区三区| 国产成人欧美| 国产精品久久久久成人av| 日韩大尺度精品在线看网址 | 久久精品亚洲av国产电影网| 国产成人精品无人区| 亚洲,欧美精品.| 自线自在国产av| 日韩精品中文字幕看吧| 精品卡一卡二卡四卡免费| 9色porny在线观看| 桃色一区二区三区在线观看| av网站免费在线观看视频| 亚洲中文av在线| 天堂动漫精品| 亚洲片人在线观看| 大型黄色视频在线免费观看| 国产视频一区二区在线看| 高清欧美精品videossex| xxxhd国产人妻xxx| 国产成人啪精品午夜网站| 国产色视频综合| 女性生殖器流出的白浆| 长腿黑丝高跟| 一进一出好大好爽视频| 日韩成人在线观看一区二区三区| 亚洲国产欧美日韩在线播放| 亚洲精品成人av观看孕妇| 亚洲一区二区三区不卡视频| 人人妻人人添人人爽欧美一区卜| 9色porny在线观看| 免费在线观看视频国产中文字幕亚洲| 午夜91福利影院| 超碰成人久久| 日韩精品中文字幕看吧| 黄色成人免费大全| 亚洲国产欧美日韩在线播放| 久久久久久人人人人人| 18禁观看日本| 最好的美女福利视频网| 80岁老熟妇乱子伦牲交| 男女床上黄色一级片免费看| 久久久国产精品麻豆| 黄色视频,在线免费观看| 久久人妻熟女aⅴ| 夜夜爽天天搞| 欧美成狂野欧美在线观看| 极品人妻少妇av视频| 久久久久国产一级毛片高清牌| 正在播放国产对白刺激| 人妻丰满熟妇av一区二区三区| 看黄色毛片网站| 日韩国内少妇激情av| 国产精品99久久99久久久不卡| 国产在线精品亚洲第一网站| 国产极品粉嫩免费观看在线| 精品欧美一区二区三区在线| 高潮久久久久久久久久久不卡| 热99国产精品久久久久久7| 欧美亚洲日本最大视频资源| 欧美在线黄色| 国产午夜精品久久久久久| 99国产精品一区二区蜜桃av| 欧美日韩一级在线毛片| 高清av免费在线| 男男h啪啪无遮挡| 一边摸一边做爽爽视频免费| 国产精品久久久av美女十八| 日韩三级视频一区二区三区| 一个人观看的视频www高清免费观看 | 久久 成人 亚洲| 大型黄色视频在线免费观看| 丝袜美腿诱惑在线| 久久久久国产一级毛片高清牌| 亚洲av五月六月丁香网| 亚洲精品一二三| 波多野结衣一区麻豆| 黄色毛片三级朝国网站| 久久人人97超碰香蕉20202| 国产精品av久久久久免费| www国产在线视频色| 日韩精品中文字幕看吧| 精品国产美女av久久久久小说| 国产精华一区二区三区| 窝窝影院91人妻| 乱人伦中国视频| 男女之事视频高清在线观看| 高清在线国产一区| 激情视频va一区二区三区| 日本黄色日本黄色录像| 性少妇av在线| 亚洲国产精品合色在线| 亚洲激情在线av| 国产色视频综合| 在线观看免费高清a一片| 日本五十路高清| 在线看a的网站| 国产高清国产精品国产三级| 精品电影一区二区在线| av在线天堂中文字幕 | 伊人久久大香线蕉亚洲五| 日韩视频一区二区在线观看| 高清av免费在线| 国产亚洲精品久久久久久毛片| 亚洲久久久国产精品| 午夜老司机福利片| 老司机深夜福利视频在线观看| 男人舔女人下体高潮全视频| 如日韩欧美国产精品一区二区三区| 可以在线观看毛片的网站| 新久久久久国产一级毛片| 亚洲av熟女| 久久国产精品人妻蜜桃| 12—13女人毛片做爰片一| 天天添夜夜摸| 国内久久婷婷六月综合欲色啪| 不卡一级毛片| 亚洲九九香蕉| 在线播放国产精品三级| 亚洲国产欧美日韩在线播放| 性色av乱码一区二区三区2| 色尼玛亚洲综合影院| 热re99久久精品国产66热6| 久久久久久大精品| 色综合欧美亚洲国产小说| 又紧又爽又黄一区二区| 亚洲av成人av| 亚洲第一欧美日韩一区二区三区| 成人特级黄色片久久久久久久| 女性生殖器流出的白浆| 精品欧美一区二区三区在线| 国产无遮挡羞羞视频在线观看| 91成年电影在线观看| 精品午夜福利视频在线观看一区| 久久香蕉精品热| videosex国产| www.精华液| 悠悠久久av| 少妇裸体淫交视频免费看高清 | 村上凉子中文字幕在线| 日本免费a在线| 亚洲成国产人片在线观看| 亚洲欧美日韩另类电影网站| 亚洲中文av在线| 精品日产1卡2卡| 亚洲成人免费电影在线观看| 国产精品自产拍在线观看55亚洲| 多毛熟女@视频| 悠悠久久av| 日本a在线网址| 黑人欧美特级aaaaaa片| 乱人伦中国视频| av超薄肉色丝袜交足视频| 一级毛片女人18水好多| 亚洲成人免费av在线播放| 夜夜爽天天搞| 热re99久久国产66热| 国产真人三级小视频在线观看| 黄片大片在线免费观看| 91在线观看av| 18禁国产床啪视频网站| 久久久精品国产亚洲av高清涩受| 一级黄色大片毛片| 90打野战视频偷拍视频| 午夜成年电影在线免费观看| 自线自在国产av| 亚洲五月天丁香| 国产精品野战在线观看 | 国产成人av激情在线播放| 午夜影院日韩av| a在线观看视频网站| 亚洲午夜理论影院| 精品人妻在线不人妻| 国产精品美女特级片免费视频播放器 | 国产精品一区二区精品视频观看| 亚洲av日韩精品久久久久久密| 欧洲精品卡2卡3卡4卡5卡区| 国产精品综合久久久久久久免费 | 黄片大片在线免费观看| 91在线观看av| 国产91精品成人一区二区三区| aaaaa片日本免费| 在线观看免费视频日本深夜| 免费观看人在逋| 淫妇啪啪啪对白视频| 国产在线观看jvid| 波多野结衣一区麻豆| 欧美精品啪啪一区二区三区| 国产精华一区二区三区| 一本大道久久a久久精品| 亚洲,欧美精品.| 国产单亲对白刺激| 黄网站色视频无遮挡免费观看| 俄罗斯特黄特色一大片| 国产在线精品亚洲第一网站| 日韩 欧美 亚洲 中文字幕| 免费在线观看黄色视频的| 九色亚洲精品在线播放| 精品久久久久久久毛片微露脸| 国产高清国产精品国产三级| 亚洲第一av免费看| 久99久视频精品免费| 97人妻天天添夜夜摸| 日本黄色日本黄色录像| 午夜久久久在线观看| 国产深夜福利视频在线观看| 男人的好看免费观看在线视频 | 在线观看午夜福利视频| 国产伦一二天堂av在线观看| 天堂动漫精品| 夫妻午夜视频| 久久国产乱子伦精品免费另类| 亚洲五月天丁香| 韩国av一区二区三区四区| 久久中文看片网| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 在线观看午夜福利视频| 又黄又粗又硬又大视频| 成年人免费黄色播放视频| 成人三级做爰电影| 18禁观看日本| 成年人黄色毛片网站| 国产精品乱码一区二三区的特点 | 国产精品秋霞免费鲁丝片| av免费在线观看网站| 大陆偷拍与自拍| 欧美久久黑人一区二区| 日本欧美视频一区| 精品熟女少妇八av免费久了| 麻豆成人av在线观看| 日韩视频一区二区在线观看| 脱女人内裤的视频| 两个人免费观看高清视频| 免费在线观看亚洲国产| 国产99白浆流出| 欧美在线一区亚洲| 黑人操中国人逼视频| 国产精品二区激情视频| 久久伊人香网站| 精品久久久久久,| 别揉我奶头~嗯~啊~动态视频| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 高清在线国产一区| 亚洲伊人色综图| 免费看十八禁软件| 欧美乱码精品一区二区三区| 天堂动漫精品| 国产亚洲欧美精品永久| 男人舔女人的私密视频| 久久精品成人免费网站| 女人爽到高潮嗷嗷叫在线视频| 欧美黄色片欧美黄色片| 久久精品国产99精品国产亚洲性色 | 午夜精品久久久久久毛片777| 免费日韩欧美在线观看| 国产精品九九99| 精品久久久久久成人av| 自线自在国产av| 久久精品成人免费网站| 一本综合久久免费| 亚洲熟妇中文字幕五十中出 | 国产xxxxx性猛交| 国产色视频综合| 欧美精品啪啪一区二区三区| 90打野战视频偷拍视频| 国产日韩一区二区三区精品不卡| 国产成人系列免费观看| 午夜精品久久久久久毛片777| 亚洲色图av天堂| 国产成人精品在线电影| 最近最新中文字幕大全电影3 | 91麻豆av在线| 久久久国产欧美日韩av| 国产精品秋霞免费鲁丝片| 老司机靠b影院| 亚洲人成电影观看| 成年版毛片免费区| 中文字幕人妻丝袜制服| 亚洲av片天天在线观看| 国产黄a三级三级三级人| 国产aⅴ精品一区二区三区波| 久久久国产欧美日韩av| 精品一区二区三区四区五区乱码| 悠悠久久av| 成人手机av| 在线观看免费视频日本深夜| 中国美女看黄片| 免费在线观看视频国产中文字幕亚洲| 亚洲九九香蕉| 国产亚洲精品综合一区在线观看 | 久久草成人影院| 亚洲成人久久性| 国产成人精品无人区| 两个人看的免费小视频| 国产又色又爽无遮挡免费看| 自拍欧美九色日韩亚洲蝌蚪91| 日韩欧美国产一区二区入口| 精品一区二区三区四区五区乱码| 国产深夜福利视频在线观看| 国产麻豆69| 大型黄色视频在线免费观看| 黄色片一级片一级黄色片| 欧美日韩亚洲高清精品| 男人操女人黄网站| 久久精品91蜜桃| 丝袜人妻中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 69av精品久久久久久| 又黄又粗又硬又大视频| 久久中文字幕人妻熟女| 91九色精品人成在线观看| 90打野战视频偷拍视频| 99国产精品99久久久久| 成年人免费黄色播放视频| 高潮久久久久久久久久久不卡| 精品乱码久久久久久99久播| 9色porny在线观看| 欧美黄色片欧美黄色片| 丰满迷人的少妇在线观看| 亚洲精品一二三| 亚洲欧洲精品一区二区精品久久久| 日本 av在线| 国产亚洲欧美精品永久| 国产不卡一卡二| 色综合欧美亚洲国产小说| 久久国产精品男人的天堂亚洲| 性欧美人与动物交配| 女人高潮潮喷娇喘18禁视频| 国产97色在线日韩免费| 亚洲自拍偷在线| 视频区图区小说| 久久人人97超碰香蕉20202| 18禁黄网站禁片午夜丰满| 午夜精品国产一区二区电影| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 精品久久蜜臀av无| 一级a爱片免费观看的视频| 可以免费在线观看a视频的电影网站| 欧美日本亚洲视频在线播放| 精品久久蜜臀av无| 国产一区在线观看成人免费| 亚洲avbb在线观看| 男女做爰动态图高潮gif福利片 | 午夜日韩欧美国产| 亚洲av电影在线进入| 精品国产国语对白av| av在线天堂中文字幕 | 久久久久久久精品吃奶| 曰老女人黄片| 国产成人影院久久av| 亚洲人成网站在线播放欧美日韩| 久久午夜亚洲精品久久| 成人av一区二区三区在线看| 精品国产一区二区久久| 日韩精品青青久久久久久| 性少妇av在线| 巨乳人妻的诱惑在线观看| 久久久国产一区二区| 在线观看免费视频网站a站| √禁漫天堂资源中文www| av视频免费观看在线观看| 中文字幕最新亚洲高清| 国产成年人精品一区二区 | 日本免费a在线| 免费日韩欧美在线观看| 欧美日韩乱码在线| 欧美性长视频在线观看| 日日干狠狠操夜夜爽| 国产精品国产高清国产av| 亚洲aⅴ乱码一区二区在线播放 | 久久这里只有精品19| 亚洲少妇的诱惑av| 国产单亲对白刺激| 精品一品国产午夜福利视频| 久久国产精品人妻蜜桃| 国产一区二区三区综合在线观看| 国产亚洲精品综合一区在线观看 | 侵犯人妻中文字幕一二三四区| 在线看a的网站| 亚洲av五月六月丁香网| 18禁裸乳无遮挡免费网站照片 | 脱女人内裤的视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品永久免费网站| 在线免费观看的www视频| 午夜免费鲁丝| 日韩视频一区二区在线观看| 757午夜福利合集在线观看| 亚洲国产欧美一区二区综合| 亚洲成人免费电影在线观看| 国产熟女午夜一区二区三区| 亚洲精品一区av在线观看| 精品久久蜜臀av无| 国产av一区二区精品久久| 国产日韩一区二区三区精品不卡| 在线天堂中文资源库| 成人精品一区二区免费| 午夜亚洲福利在线播放| 97人妻天天添夜夜摸| 黄片大片在线免费观看| ponron亚洲| av片东京热男人的天堂| 欧美在线黄色| 久久国产精品人妻蜜桃| 亚洲精品国产精品久久久不卡| 精品国产国语对白av| 亚洲七黄色美女视频| 欧美精品亚洲一区二区| 午夜精品在线福利| 久久欧美精品欧美久久欧美| 91麻豆av在线| 欧美激情 高清一区二区三区| 亚洲精品久久午夜乱码| 午夜福利免费观看在线| 国产精品乱码一区二三区的特点 | 97碰自拍视频| 中国美女看黄片| 妹子高潮喷水视频| 日本一区二区免费在线视频| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 在线天堂中文资源库| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av| 精品国产超薄肉色丝袜足j| 99在线视频只有这里精品首页| 99久久人妻综合| 国产精品亚洲av一区麻豆| 久久亚洲精品不卡| 看免费av毛片| 国内毛片毛片毛片毛片毛片| av免费在线观看网站| 午夜精品在线福利| 最好的美女福利视频网| 成人18禁在线播放| 亚洲中文av在线| 老司机深夜福利视频在线观看| 久久精品91无色码中文字幕| 亚洲国产欧美日韩在线播放| 亚洲免费av在线视频| 黄色女人牲交| 久久国产亚洲av麻豆专区| 男女之事视频高清在线观看| 亚洲人成网站在线播放欧美日韩| 午夜老司机福利片| 国产男靠女视频免费网站| 校园春色视频在线观看| 亚洲成人国产一区在线观看| xxxhd国产人妻xxx| 一夜夜www| 夜夜看夜夜爽夜夜摸 | 黄网站色视频无遮挡免费观看| 女警被强在线播放| 日本欧美视频一区| 我的亚洲天堂| 亚洲九九香蕉| 亚洲熟女毛片儿| 精品国产一区二区久久| 国产高清激情床上av| 日韩大尺度精品在线看网址 | 中文字幕色久视频| 美女高潮喷水抽搐中文字幕| 韩国av一区二区三区四区| 中亚洲国语对白在线视频| 老熟妇仑乱视频hdxx| 欧美一级毛片孕妇| 亚洲精品美女久久久久99蜜臀| 成人精品一区二区免费| 一级片免费观看大全| 一进一出抽搐gif免费好疼 | 国产欧美日韩精品亚洲av| 黑人操中国人逼视频| 中亚洲国语对白在线视频| 国产精品综合久久久久久久免费 | 久久人人爽av亚洲精品天堂| 99精国产麻豆久久婷婷| 国产精品亚洲一级av第二区| 午夜老司机福利片| av视频免费观看在线观看| 麻豆一二三区av精品| 亚洲国产精品一区二区三区在线| 超碰97精品在线观看| 性欧美人与动物交配| 日韩欧美三级三区| 女人爽到高潮嗷嗷叫在线视频| 国产黄色免费在线视频| av中文乱码字幕在线| 国产精品永久免费网站| 露出奶头的视频| av免费在线观看网站| 国产极品粉嫩免费观看在线| 一区二区日韩欧美中文字幕| 国产成人免费无遮挡视频| 免费人成视频x8x8入口观看| 国产成人欧美| 欧美日韩亚洲国产一区二区在线观看| 成年人免费黄色播放视频| 久久久久国产精品人妻aⅴ院| 麻豆国产av国片精品| 亚洲国产精品sss在线观看 | 男女做爰动态图高潮gif福利片 | 精品福利观看| 国产极品粉嫩免费观看在线| 一区二区日韩欧美中文字幕| 国产成人精品久久二区二区免费| 色精品久久人妻99蜜桃| 嫩草影视91久久| 国内久久婷婷六月综合欲色啪| 91成年电影在线观看| 久久国产精品影院| e午夜精品久久久久久久| 国产精品免费一区二区三区在线| 午夜福利影视在线免费观看| 成年人黄色毛片网站| 亚洲第一av免费看| 可以在线观看毛片的网站| 交换朋友夫妻互换小说| 伊人久久大香线蕉亚洲五| 高清在线国产一区| 叶爱在线成人免费视频播放| 999久久久精品免费观看国产| 狂野欧美激情性xxxx| 日本黄色日本黄色录像| 老司机在亚洲福利影院| 精品一区二区三区四区五区乱码| 国产91精品成人一区二区三区| 午夜福利影视在线免费观看| 如日韩欧美国产精品一区二区三区| 搡老岳熟女国产| 亚洲色图综合在线观看| 精品人妻1区二区| 91大片在线观看| 嫩草影院精品99| 日韩中文字幕欧美一区二区| 首页视频小说图片口味搜索| 欧美日本中文国产一区发布| 国产精品99久久99久久久不卡| 精品一区二区三区四区五区乱码| a级片在线免费高清观看视频| 久久精品亚洲精品国产色婷小说| 一区在线观看完整版| 国产又色又爽无遮挡免费看| 中文字幕高清在线视频| 久久久久九九精品影院| 日本欧美视频一区| 真人做人爱边吃奶动态| 国产成人影院久久av| 成人三级黄色视频| 国产精品野战在线观看 | 嫩草影视91久久| 最新美女视频免费是黄的| 久久人妻熟女aⅴ| ponron亚洲| 亚洲欧洲精品一区二区精品久久久| 香蕉丝袜av| 免费在线观看影片大全网站| 午夜成年电影在线免费观看| 99国产综合亚洲精品| 久久婷婷成人综合色麻豆| 欧美日韩亚洲国产一区二区在线观看| 90打野战视频偷拍视频| 丰满饥渴人妻一区二区三| 国产亚洲欧美精品永久| 免费高清在线观看日韩| 久久久久久免费高清国产稀缺| 99久久综合精品五月天人人| 五月开心婷婷网| 久久香蕉精品热| 色综合欧美亚洲国产小说| 亚洲av电影在线进入| 亚洲精品国产区一区二| 欧美中文综合在线视频| 国产精品av久久久久免费| 12—13女人毛片做爰片一| 国产片内射在线| 亚洲七黄色美女视频| 久久欧美精品欧美久久欧美| 我的亚洲天堂| 一级,二级,三级黄色视频|