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

    爆轟驅(qū)動下軸向預(yù)制破片飛散特性研究

    2023-10-14 01:00:56李正川王傳婷
    彈道學(xué)報 2023年3期
    關(guān)鍵詞:戰(zhàn)斗部內(nèi)層破片

    李正川,何 源,郭 磊,王傳婷,何 勇,馬 躍

    (南京理工大學(xué) 機械工程學(xué)院,江蘇 南京 210094)

    傳統(tǒng)預(yù)制破片戰(zhàn)斗部是將破片布置在裝藥周圍,炸藥爆炸后驅(qū)動周向破片高速飛出,利用破片的高速動能實現(xiàn)對目標結(jié)構(gòu)的毀傷,然而破片周向排列方式將導(dǎo)致戰(zhàn)斗部在單一方向上有效破片數(shù)量少,打擊效率降低。為更大效率地發(fā)揮戰(zhàn)斗部對固定方向上的毀傷威力,有學(xué)者提出了定向戰(zhàn)斗部[1]的概念,軸向預(yù)制破片戰(zhàn)斗部是其中的一種,它將預(yù)制破片排列于裝藥端面,引爆裝藥后,預(yù)制破片會沿裝藥軸線方向向前運動,具有定向毀傷攔截的特點,并且可以改變炸藥端面的曲率以控制破片場的有效毀傷范圍。目前,以軸向預(yù)制破片形式設(shè)計的攔截彈已應(yīng)用于武裝直升機等機載主動防御系統(tǒng)中,如美國“直升機主動RPG防護”系統(tǒng)[2]。

    為了研究軸向戰(zhàn)斗部多層預(yù)制破片層間速度大小、飛散角分布的影響因素,國內(nèi)外學(xué)者針對軸向預(yù)制破片戰(zhàn)斗部破片飛散情況展開了系統(tǒng)的研究。JONES[3]基于能量守恒定律和Gurney假設(shè)推導(dǎo)了多層破片速度的理論計算方法,為爆轟驅(qū)動多層預(yù)制破片的相關(guān)研究提供了理論支撐。文獻[4-5]針對單層軸向預(yù)制破片戰(zhàn)斗部進行了研究,基于Gurney公式,得到了軸向單層預(yù)制破片速度和飛散角的計算公式。李明星[6]利用理論和數(shù)值方法對單層球形、方形和圓柱形軸向預(yù)制破片的飛散特性進行研究,結(jié)果表明球形破片速度和飛散角均大于其他形狀破片。MA等[7]針對軸向單層預(yù)制破片在不同襯層下的爆轟驅(qū)動過程進行了試驗和數(shù)值研究,并與無襯層形式進行對比,結(jié)果表明,凱夫拉/環(huán)氧樹脂襯層使得破片的初速度低于無襯層形式下的破片速度,并且球形破片在爆轟驅(qū)動方面優(yōu)于圓柱形破片。文獻[8-9]基于改進后的Gurney公式得到了多層球形破片在爆轟驅(qū)動后的破片平均速度,結(jié)合Taylor公式得到了破片的飛散范圍角,并通過實驗驗證了理論結(jié)果的合理性。文獻[10-11]研究了弧形裝藥結(jié)構(gòu)下多層球形破片的飛散情況,得到了帶弧形圓柱裝藥下,從內(nèi)往外球形破片的平均初速逐層增加,以及隨著球形破片直徑減小最內(nèi)層破片初速增大的規(guī)律。

    目前,有關(guān)多層軸向預(yù)制破片的研究通常將金屬破片層看作一個金屬板[12-14],分析金屬板的初速得到破片的近似速度,并未涉及到破片層間的相互作用關(guān)系。本文針對多層預(yù)制破片在炸藥爆轟驅(qū)動后的層間速度關(guān)系問題,基于Gurney假設(shè)和動量守恒定律得到了軸向預(yù)制破片戰(zhàn)斗部多層破片初速的近似計算模型,并進行了破片爆轟驅(qū)動過程動力學(xué)仿真,分析了不同破片形狀對于多層破片初速的影響,以及破片層數(shù)對于破片初速的影響,為提高戰(zhàn)斗部的毀傷威力設(shè)計提供一定的參考作用。

    1 破片層初速理論計算模型

    1.1 速度計算模型

    本文采用圓柱形裝藥結(jié)構(gòu),不同形狀的預(yù)制破片排布在炸藥端部的隔板上,并用樹脂黏結(jié)劑固化黏結(jié)在一起。為建立不同形狀的多層破片初速理論計算公式,作如下假設(shè):

    ①破片前置且按軸向分層排列,球形破片和圓柱形破片均采用密排方式,破片間隙采用樹脂黏結(jié)劑填充,忽略其對破片飛散特性的影響,只考慮質(zhì)量。炸藥采用中心點起爆方式。

    ②炸藥瞬時爆炸,釋放的能量全部用于破片和爆轟產(chǎn)物的飛散,爆轟產(chǎn)物不通過預(yù)制破片的間隙向外溢出。

    ③爆轟產(chǎn)物從起爆端面到球形破片內(nèi)切面、圓柱破片底面、方形破片底面之間的速度呈線性分布,并且與第一層破片接觸的爆轟產(chǎn)物的速度與該層破片速度相同。

    在炸藥爆炸驅(qū)動多層破片飛散過程中,分別用mi,vi(i=1,2,…,n)表示從靠近炸藥到遠離炸藥的第i層破片的質(zhì)量和速度,由于相鄰破片層之間的排列方式相同,假設(shè)相鄰破片層之間的速度比值相同,即vi+1/vi=k。

    圖1 理論模型示意圖Fig.1 Schematic diagram of the theoretical model

    由動量守恒得到預(yù)制破片與爆轟產(chǎn)物的運動方程:

    (1)

    式中:mc為炸藥爆炸前的質(zhì)量,mc/2為炸藥爆炸后爆轟產(chǎn)物的虛擬質(zhì)量,m為破片層中樹脂黏結(jié)劑、隔板和殼體的總質(zhì)量,S為爆轟產(chǎn)物作用于端面第一層預(yù)制破片的有效面積,p為某瞬間爆轟產(chǎn)物作用于預(yù)制破片層的壓強。同時,假設(shè)炸藥半徑為r0,炸藥密度為ρ0,r為爆轟產(chǎn)物的膨脹后與裝藥軸線的距離。

    (2)

    設(shè)S=λS0,λ為軸向預(yù)制破片數(shù)量和r的因變量。當預(yù)制破片結(jié)構(gòu)尺寸和多層破片排布方式確定后,定有常數(shù)0<ε≤1,使得:

    (3)

    將式(2)、式(3)代入式(1),并整理可得到:

    (4)

    式中:r≥r0。所以最內(nèi)層破片初速為

    (5)

    對于球形破片,有:

    (6)

    對于方形破片,有:

    (7)

    對于圓柱破片,有:

    (8)

    vi=k(i-1)v1=k(i-1)aψ(n1)v0

    式中:k與破片層數(shù)i、最內(nèi)層破片數(shù)n1有關(guān)。針對于不同形狀多層破片的k值,可通過數(shù)值仿真的結(jié)果來確定。

    1.2 破片飛散角計算模型

    圖2 飛散角計算示意圖Fig.2 Schematic diagram of dispersion angle calculation

    2 仿真分析

    2.1 戰(zhàn)斗部有限元模型

    本文使用非線性有限元計算程序LS-DYNA和網(wǎng)格劃分軟件True-Grid分別對球形、方形、圓柱形3種破片進行爆轟驅(qū)動數(shù)值模擬分析,戰(zhàn)斗部結(jié)構(gòu)包括炸藥、空氣、隔板、殼和預(yù)制破片。戰(zhàn)斗部數(shù)據(jù)如表1所示。戰(zhàn)斗部殼體厚度為1 mm,各工況下殼體質(zhì)量如表2所示。

    表1 戰(zhàn)斗部數(shù)據(jù)Table 1 Warhead data

    表2 各工況殼體質(zhì)量Table 2 Shell quality under different conditions

    其中,每層球形破片和圓柱破片均為密排方式,樹脂黏結(jié)劑采用環(huán)氧樹脂。以球形破片為例,戰(zhàn)斗部模型如圖3所示,在各工況下,將球形、方形和圓柱形多層破片的每層破片按圖4的方式分為不同列,并記為第1,2,…,i列。

    圖3 兩層球破片戰(zhàn)斗部示意圖Fig.3 Schematic diagram of two-layer ball-fragment warhead

    圖4 破片分列示意圖Fig.4 Schematic diagram of fragmentary distribution

    為節(jié)約仿真計算時間,提高計算效率,考慮到戰(zhàn)斗部的對稱關(guān)系,采用1/4戰(zhàn)斗部模型進行仿真分析。破片、殼體和隔板采用Lagrange網(wǎng)格,炸藥、空氣和樹脂黏結(jié)劑采用Euler網(wǎng)格。針對流體網(wǎng)格和固體網(wǎng)格在數(shù)值計算中的優(yōu)缺點,本文采用ALE多材料、流體/結(jié)構(gòu)耦合算法進行計算。對于Lagrange區(qū)域重疊歐拉網(wǎng)格區(qū)域的問題,本文仿真中采用*CONSTRAINED_LAGRANGE_IN_SOLID關(guān)鍵字來實現(xiàn)流體-結(jié)構(gòu)間的相互作用,拉格朗日網(wǎng)格與ALE(流體)網(wǎng)格在空間上重疊,每個網(wǎng)格都用獨立的節(jié)點id來定義,LS-DYNA搜索拉格朗日網(wǎng)格和ALE網(wǎng)格之間的空間交點,當網(wǎng)格重疊時,發(fā)生相互作用。破片與破片、隔板與破片以及破片和殼體間設(shè)置為自動單面接觸,炸藥采用單點起爆方式,起爆點設(shè)置在炸藥底面圓心處。

    2.2 材料模型和參數(shù)

    空氣材料采用MAT_NULL模型和線性多項式狀態(tài)方程(EOS_LINEAR_POLYNOMIAL)。炸藥采用高能炸藥燃燒材料模型(HIGH_EXPLOSIVE_BURN)和JWL狀態(tài)方程,材料參數(shù)如表3所示。表中,pCJ為爆轟壓力;A,B,R1,R2,ω為炸藥材料參數(shù)。

    表3 炸藥材料參數(shù)Table 3 Material parameters of explosive

    殼體和隔板材料模型均采用Johnson-Cook模型和Gruneisen狀態(tài)方程,材料參數(shù)如表4所示。表中,A為材料初始屈服應(yīng)力,B為應(yīng)變硬化系數(shù),n為應(yīng)變硬化指數(shù),C為應(yīng)變率系數(shù),m為熱軟化指數(shù)。

    表4 殼體和隔板材料參數(shù)Table 4 Material parameters of shell and baffle

    破片材料采用MAT_PLASTIC_KINEMATIC模型,并添加沙漏(HOURGLASS)控制,材料參數(shù)如表5所示。表中,E為楊氏模量,ν為泊松比,σ0為屈服應(yīng)力。

    表5 破片材料參數(shù)Table 5 Material parameters of explosive

    樹脂材料采用MAT_ELASTIC_PLASTIC_HYDRO模型和GRUNEISEN狀態(tài)方程,材料參數(shù)如表6所示。表中,G為剪切模量,σ為屈服應(yīng)力。

    表6 環(huán)氧樹脂材料參數(shù)Table 6 Material parameters of epoxy resin

    3 分析與討論

    3.1 破片飛散仿真分析

    本文針對球形、方形和圓柱形破片的爆轟驅(qū)動數(shù)值仿真分別設(shè)置了兩層和三層破片的工況,共6種工況。由于各工況下爆轟產(chǎn)物與破片間的作用效果相似,故本文以球形破片戰(zhàn)斗部為例,爆轟產(chǎn)物與破片間的作用示意圖如圖5所示??梢钥闯銎破c爆轟產(chǎn)物之間存在相互作用關(guān)系,破片在爆轟產(chǎn)物的作用下向前運動。

    圖5 爆轟產(chǎn)物與破片相互作用示意圖Fig.5 Schematic diagram of interaction between detonation product and fragment

    不同形狀破片情況下爆轟波的傳播情況相似,各工況下多層破片爆轟驅(qū)動飛散過程如圖6所示。從圖6可以看出爆轟波在t=90 μs左右時到達最內(nèi)層破片的下表面,戰(zhàn)斗部軸線周圍破片首先受到爆轟波的推動作用開始向前運動,t=180 μs時,破片已在空氣中運動。炸藥爆炸后產(chǎn)生的爆轟波以及爆轟產(chǎn)物首先作用到隔板上,由于在極短的時間內(nèi)隔板處爆轟產(chǎn)物壓力急劇增大,最高可達61.39 GPa,且由于隔板薄以及隔板本身材料的原因,隔板迅速變形以致破壞失效,然后爆轟產(chǎn)物再作用到第一層破片下表面,破片開始向前飛散并與第二層破片發(fā)生碰撞,相同時刻下,外層破片飛行距離較遠。

    圖6 各工況下多層破片飛散示意圖(單位:GPa)Fig.6 Schematic diagram of multi-layer fragmentation dispersion under various working conditions(unit:GPa)

    圖7為破片碰撞速度變化示意圖。可以看出,第一層破片在與第二層破片發(fā)生碰撞后,第一層破片速度會發(fā)生衰減,同樣第二層破片在和第三層破片發(fā)生碰撞后也會發(fā)生衰減,然后在爆轟產(chǎn)物的作用下繼續(xù)加速,最后各層破片速度達到穩(wěn)定的狀態(tài)。同時,第一層破片在與第二層破片、第二層破片與第三層破片發(fā)生碰撞前,第二層破片和第三層破片已經(jīng)受到爆轟波的加速作用。

    圖7 破片碰撞速度變化示意圖Fig.7 Variation of fragments collision velocity

    圖8為不同工況同一時刻下,多層破片飛散速度云圖??梢钥闯?炸藥爆炸后,破片會沿著戰(zhàn)斗部軸向和徑向飛散,由于爆轟波先到達戰(zhàn)斗部軸線附近的破片,故相同時刻下破片的最大速度出現(xiàn)在破片層中心位置。從各工況的速度云圖可以看出,最大動能均出現(xiàn)在破片層的中心位置,也即紅色區(qū)域,從內(nèi)往外破片動能減少,兩層球破片和三層球破片的第一層邊緣破片最大平均速度分別為977.6 m/s和749.3 m/s,超過了藍色區(qū)域破片的速度,并且從圖7(a)中可以明顯看到,第一層邊緣破片已超過第二層藍色區(qū)域破片,其余圖中也可看出邊緣破片超過非邊緣破片的趨勢,這是由于爆轟產(chǎn)物在破片層邊緣位置發(fā)生繞流,并且邊緣位置破片約束較少,進一步加速了邊緣位置的破片,其中球形破片受到繞流加速作用的效果最為明顯。

    圖8 多層破片飛散示意圖(單位:m/s)Fig.8 Schematic diagram of multi-layer fragment dispersion(unit:m/s)

    3.2 數(shù)值仿真與理論計算結(jié)果分析

    圖9為破片速度穩(wěn)定后,各層各列破片速度平均值以及每層破片速度平均值??梢钥闯龈鞴r下的數(shù)值仿真結(jié)果中,各層破片速度分布合理,沿半徑方向同列破片速度梯度越來越大,兩層破片情況下,各列破片最大速度差值:球破片為42.9 m/s,方破片為30.78 m/s,圓柱破片為12 m/s;三層破片情況下,第一層和第二層破片間,各列破片最大速度差值:球破片為48.2 m/s,方破片為19.2 m/s,圓柱破片為8.3 m/s,;第二層和第三層破片間,各列破片最大速度差值:球破片為57.2 m/s,方破片為28.2 m/s,圓柱破片為10.1 m/s,顯然球破片的速度梯度相對方形和圓柱形破片更大,同時,對于不同類型破片,以上差值分別取自不同的列位置,其中只有球形破片取自破片層邊緣位置,說明在爆轟產(chǎn)物繞流作用、邊緣位置破片束縛較少以及球形破片受到的阻力相對于方形和圓柱形較小的多方影響因素下,同列位置內(nèi)層破片速度更高于外層破片。三層球形破片和方形破片的工況仿真結(jié)果中,存在明顯的破片層邊緣處破片受到繞流加速的作用后速度變大的情況,與前文分析一致。

    圖9 各列破片速度Fig.9 Fragment velocity of each column

    表7~表9為不同工況下破片速度計算值和仿真值??梢钥闯?對應(yīng)層數(shù)破片工況下,球形破片初速平均值大于圓柱破片初速平均值,對應(yīng)層數(shù)圓柱破片初速平均值大于方破片初速平均值,其原因為:①球形破片相對于其他兩種類型破片的平均迎風(fēng)面積小,受到的空氣阻力小,導(dǎo)致在相同飛散距離內(nèi)破片速度衰減較弱;②在相同破片材料情況下,單枚破片質(zhì)量均相同,而不同形狀破片的空間占有率不同,所以相同層數(shù)破片數(shù)目不同,最終戰(zhàn)斗部裝填比不同,導(dǎo)致球形破片一開始初速也會高于其他破片;③邊緣位置破片在爆轟產(chǎn)物繞流加速后導(dǎo)致整體速度平均值上升。兩層破片情況下,圓柱破片速度計算值與仿真值最為接近,誤差為2.2%,三層破片情況下,球破片速度計算值與仿真值最為接近,誤差為1.2%;隨著破片層數(shù)以及破片數(shù)量的增加,最內(nèi)層破片初速逐漸減小;相鄰層之間破片平均初速的比值幾乎相等。

    表7 球破片初速理論計算值和仿真值Table 7 Theoretical calculation and simulation values of initial velocity of ball fragments

    表8 圓柱破片初速理論計算值和仿真值Table 8 Theoretical calculation and simulation values of initial velocity of cylinder fragments

    表9 方破片初速理論計算值和仿真值Table 9 Theoretical calculation and simulation values of initial velocity of square fragments

    3.3 破片飛散角討論

    仿真結(jié)果結(jié)合簡化平均飛散角計算模型得到的破片飛散角如圖10和圖11所示。

    從圖中可以看出,兩層破片情況下,方破片內(nèi)層破片飛散角度小于外層破片,同層破片從內(nèi)往外破片飛散角度先增大后減小;對于球破片,整體上內(nèi)層破片飛散角小于外層破片的規(guī)律,而在破片層邊緣位置,則存在內(nèi)層破片飛散角大于外層破片的情況,這也說明了內(nèi)層邊緣破片受到爆轟產(chǎn)物繞流加速作用后速度高于外層破片;圓柱破片內(nèi)層和外層破片飛散角較為一致,均在17°~22°范圍內(nèi),不同層破片間飛散角差值較小,同層破片的每列破片飛散角分布較為均勻。三層破片情況下,方破片飛散角度從內(nèi)層破片到外層破片逐漸增大,對于同層破片,沿半徑方向破片飛散角呈先增大后減小的趨勢;對于球破片而言,最外層破片飛散角均小于同列內(nèi)層破片,而沿半徑方向破片飛散角逐漸增大;同列圓柱破片從最內(nèi)層到最外層破片飛散角呈逐層增加的趨勢,但整體上內(nèi)層和外層破片角度分布在19°~24°范圍內(nèi),同層每列破片的飛散角分布較為均勻。

    整體而言,對于不同類型破片的飛散角存在一定的差異,圓柱破片飛散角較為均勻,球破片沿半徑方向飛散角逐漸增大,而方破片的邊緣位置破片飛散角為最小值,本文認為是邊緣位置破片數(shù)量較少導(dǎo)致數(shù)據(jù)存在一定的偏差,再者考慮到破片形狀的差異、迎風(fēng)面積的影響導(dǎo)致各類型破片的飛散程度存在差異。

    4 結(jié)束語

    本文基于Gurney假設(shè)和動量守恒,建立了軸向預(yù)制破片不同層間的初速近似計算簡化模型,首先通過理論計算,得到多層破片的第一層破片初速,然后利用LS-DYNA軟件進行數(shù)值仿真,得到多層破片各層破片速度平均值和層間速度比值,最后結(jié)合層間速度比值得到各層破片的平均速度理論計算值。研究結(jié)果表明:

    ①沿半徑方向破片速度梯度先減小后變大,相同層數(shù)時,球形破片速度梯度最大;隨著破片層數(shù)增加,破片速度梯度增加;相同形狀破片的平均速度從最內(nèi)層到最外層逐漸增加。

    ②對于三種破片形狀,球形破片的平均速度最大,并且層間速度差值最大,兩層破片排布與三層破片排布時層間速度差值相接近;球破片、圓柱破片和方破片的理論模型計算值與仿真值間的誤差最低分別為1.2%,2.2%,8.8%,再次說明速度計算模型可以較好反映出圓柱裝藥下不同形狀軸向多層預(yù)制破片初速情況。

    ③各列圓柱破片飛散角分布在相同范圍內(nèi),球破片沿半徑方向飛散角有增大趨勢;在相同層數(shù)情況下,方破片的飛散角度小于圓柱破片和球破片,破片飛散角的分布規(guī)律和破片形狀密切相關(guān)。綜合破片平均速度和飛散角可以認為球形破片飛散性能較好。

    猜你喜歡
    戰(zhàn)斗部內(nèi)層破片
    ◆ 裝飾板材
    ◆ 裝飾板材
    裝飾板材
    ◆ 裝飾板材
    破片群作用下復(fù)合材料層合板近場動力學(xué)損傷模擬*
    爆炸與沖擊(2022年3期)2022-04-11 03:14:58
    彈藥動態(tài)加載下破片測試方法
    中國測試(2018年10期)2018-11-17 01:58:40
    半預(yù)制破片戰(zhàn)斗部破片威力分析
    隨動定向戰(zhàn)斗部的抗振動沖擊控制研究
    幾種結(jié)構(gòu)的三段離散桿戰(zhàn)斗部數(shù)值模擬
    前向戰(zhàn)斗部破片散布均勻性研究*
    国产伦精品一区二区三区视频9 | 美女高潮喷水抽搐中文字幕| 啪啪无遮挡十八禁网站| e午夜精品久久久久久久| 亚洲av电影不卡..在线观看| 国产成人欧美在线观看| 精品99又大又爽又粗少妇毛片 | 久久久精品大字幕| 男女那种视频在线观看| 色精品久久人妻99蜜桃| 又大又爽又粗| 中国美女看黄片| 麻豆成人午夜福利视频| 午夜精品久久久久久毛片777| 无人区码免费观看不卡| 在线观看日韩欧美| 美女 人体艺术 gogo| 欧美中文日本在线观看视频| 亚洲精品乱码久久久v下载方式 | 一本综合久久免费| 婷婷精品国产亚洲av| 亚洲国产精品999在线| 亚洲精品粉嫩美女一区| 亚洲成人久久爱视频| 两个人看的免费小视频| 色哟哟哟哟哟哟| 午夜免费观看网址| 男人的好看免费观看在线视频| 国产亚洲欧美98| 亚洲中文日韩欧美视频| 久久久久精品国产欧美久久久| 亚洲av片天天在线观看| 国产男靠女视频免费网站| 国产一区二区在线av高清观看| 久久九九热精品免费| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品合色在线| 草草在线视频免费看| 男女之事视频高清在线观看| 久久久国产成人免费| 亚洲在线观看片| 国产亚洲av高清不卡| 国产真人三级小视频在线观看| 免费大片18禁| 色综合欧美亚洲国产小说| 人人妻人人看人人澡| 国产高清视频在线播放一区| 中文字幕熟女人妻在线| 久久久国产成人精品二区| 国产亚洲精品av在线| 欧美+亚洲+日韩+国产| 久久香蕉国产精品| 最新在线观看一区二区三区| 天天一区二区日本电影三级| av国产免费在线观看| 色综合站精品国产| 亚洲一区高清亚洲精品| 又黄又爽又免费观看的视频| 欧美午夜高清在线| 久久午夜亚洲精品久久| 国产伦精品一区二区三区四那| av在线蜜桃| 亚洲精品色激情综合| 在线观看午夜福利视频| 午夜精品一区二区三区免费看| 亚洲中文字幕日韩| 国产黄a三级三级三级人| 后天国语完整版免费观看| 国内精品久久久久精免费| 午夜福利在线观看免费完整高清在 | 伦理电影免费视频| 色哟哟哟哟哟哟| 久久精品影院6| 久久久色成人| 国产成人精品久久二区二区免费| 少妇人妻一区二区三区视频| 免费搜索国产男女视频| 久久热在线av| 国产精品影院久久| 欧美成人性av电影在线观看| av黄色大香蕉| 12—13女人毛片做爰片一| 变态另类丝袜制服| 日本黄色片子视频| 国产精品99久久99久久久不卡| 婷婷六月久久综合丁香| 后天国语完整版免费观看| 窝窝影院91人妻| 久久久久久九九精品二区国产| 久久久久久久午夜电影| 久久久久国产一级毛片高清牌| 免费看光身美女| 亚洲av片天天在线观看| 精品国产超薄肉色丝袜足j| 国产私拍福利视频在线观看| 午夜福利视频1000在线观看| 成人精品一区二区免费| av国产免费在线观看| 免费搜索国产男女视频| 亚洲在线自拍视频| 波多野结衣高清无吗| 舔av片在线| 色视频www国产| 熟妇人妻久久中文字幕3abv| 久久婷婷人人爽人人干人人爱| 午夜福利在线观看吧| 观看免费一级毛片| 午夜福利在线在线| 欧美黑人巨大hd| 国内久久婷婷六月综合欲色啪| 在线a可以看的网站| 91麻豆精品激情在线观看国产| 给我免费播放毛片高清在线观看| 9191精品国产免费久久| 久久婷婷人人爽人人干人人爱| 黄片大片在线免费观看| 亚洲av片天天在线观看| 九九久久精品国产亚洲av麻豆 | 99国产精品一区二区三区| 中文字幕人成人乱码亚洲影| 毛片女人毛片| 香蕉久久夜色| 香蕉av资源在线| 国产亚洲av嫩草精品影院| 免费观看精品视频网站| 午夜影院日韩av| 日本a在线网址| 法律面前人人平等表现在哪些方面| 亚洲人成网站在线播放欧美日韩| 久久久水蜜桃国产精品网| 美女 人体艺术 gogo| 一个人看视频在线观看www免费 | 麻豆国产av国片精品| 麻豆一二三区av精品| 无遮挡黄片免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品乱码久久久v下载方式 | 中文字幕熟女人妻在线| 日本在线视频免费播放| 一进一出抽搐gif免费好疼| 90打野战视频偷拍视频| 蜜桃久久精品国产亚洲av| 夜夜看夜夜爽夜夜摸| 久久精品综合一区二区三区| svipshipincom国产片| 欧美日韩黄片免| 久久久国产精品麻豆| 免费看日本二区| 欧美日韩精品网址| 变态另类丝袜制服| 午夜福利成人在线免费观看| 身体一侧抽搐| 一夜夜www| 中文字幕最新亚洲高清| 亚洲精品在线观看二区| 亚洲精品久久国产高清桃花| 欧美高清成人免费视频www| 国产真人三级小视频在线观看| 99riav亚洲国产免费| 99久久精品一区二区三区| 免费看a级黄色片| 天天躁日日操中文字幕| 真人做人爱边吃奶动态| 国产不卡一卡二| 午夜福利视频1000在线观看| 亚洲国产欧美网| 搡老岳熟女国产| 欧美xxxx黑人xx丫x性爽| 亚洲在线观看片| 国产av在哪里看| 日本 欧美在线| 亚洲 欧美一区二区三区| 色在线成人网| 男插女下体视频免费在线播放| 丁香六月欧美| 久久99热这里只有精品18| 1024手机看黄色片| 成人高潮视频无遮挡免费网站| 久久中文看片网| 女警被强在线播放| 国产精品99久久99久久久不卡| АⅤ资源中文在线天堂| 狂野欧美激情性xxxx| 变态另类成人亚洲欧美熟女| 90打野战视频偷拍视频| 九色国产91popny在线| 成人国产一区最新在线观看| 91久久精品国产一区二区成人 | 欧美日本视频| 美女黄网站色视频| 免费看日本二区| 亚洲国产高清在线一区二区三| 午夜福利免费观看在线| 久久人妻av系列| 此物有八面人人有两片| 国产成人影院久久av| 亚洲av中文字字幕乱码综合| 色吧在线观看| www国产在线视频色| 国产午夜福利久久久久久| 国产精品女同一区二区软件 | 亚洲av成人av| 日本三级黄在线观看| 国内精品美女久久久久久| 99国产综合亚洲精品| 高清在线国产一区| 日韩欧美精品v在线| 青草久久国产| 又黄又爽又免费观看的视频| 国产一区二区三区在线臀色熟女| xxx96com| 午夜福利欧美成人| 99热6这里只有精品| 成年版毛片免费区| 日韩欧美免费精品| 亚洲 欧美一区二区三区| 五月伊人婷婷丁香| 国产美女午夜福利| 精品国产超薄肉色丝袜足j| 制服丝袜大香蕉在线| 亚洲人成网站高清观看| 男人舔女人下体高潮全视频| 久久久国产成人免费| 国产精品国产高清国产av| 变态另类丝袜制服| 真实男女啪啪啪动态图| 搡老岳熟女国产| 99在线视频只有这里精品首页| 日韩国内少妇激情av| 天堂动漫精品| 日本三级黄在线观看| 日韩欧美国产一区二区入口| 国产亚洲av嫩草精品影院| 国产免费男女视频| 久久精品综合一区二区三区| 国产成人av激情在线播放| 97人妻精品一区二区三区麻豆| 精品一区二区三区视频在线 | 久久亚洲精品不卡| 成年免费大片在线观看| 男人舔女人的私密视频| 美女cb高潮喷水在线观看 | 麻豆av在线久日| 在线观看舔阴道视频| 深夜精品福利| 久久伊人香网站| 亚洲av成人不卡在线观看播放网| 国产精华一区二区三区| 757午夜福利合集在线观看| 十八禁网站免费在线| 99久国产av精品| 999久久久国产精品视频| 中国美女看黄片| 亚洲一区二区三区不卡视频| 九九热线精品视视频播放| 丰满人妻熟妇乱又伦精品不卡| 免费看a级黄色片| 亚洲成av人片免费观看| 久久国产精品影院| 免费搜索国产男女视频| 天天一区二区日本电影三级| 偷拍熟女少妇极品色| 亚洲在线观看片| 最新中文字幕久久久久 | 日本黄大片高清| 欧美成人性av电影在线观看| 午夜日韩欧美国产| 亚洲精品美女久久av网站| 国产亚洲av高清不卡| 老熟妇仑乱视频hdxx| 天天躁日日操中文字幕| 757午夜福利合集在线观看| 久久久国产精品麻豆| xxx96com| 宅男免费午夜| 国产成人欧美在线观看| 国产男靠女视频免费网站| 给我免费播放毛片高清在线观看| 亚洲午夜精品一区,二区,三区| 中文字幕熟女人妻在线| 成人精品一区二区免费| 老汉色av国产亚洲站长工具| 观看免费一级毛片| 亚洲激情在线av| 中出人妻视频一区二区| 欧美一区二区精品小视频在线| a级毛片在线看网站| 韩国av一区二区三区四区| 99在线人妻在线中文字幕| 成人18禁在线播放| 久久亚洲精品不卡| 国产激情久久老熟女| av天堂在线播放| 香蕉国产在线看| 亚洲精品一卡2卡三卡4卡5卡| 偷拍熟女少妇极品色| 免费一级毛片在线播放高清视频| 丰满人妻一区二区三区视频av | av黄色大香蕉| 久久这里只有精品中国| www.精华液| 色在线成人网| 欧美丝袜亚洲另类 | 国产成人一区二区三区免费视频网站| 一区二区三区激情视频| 在线观看午夜福利视频| 在线a可以看的网站| 国产成+人综合+亚洲专区| 日韩三级视频一区二区三区| 成人一区二区视频在线观看| 欧美黑人欧美精品刺激| 一本精品99久久精品77| 国产蜜桃级精品一区二区三区| 欧美成人免费av一区二区三区| 国产男靠女视频免费网站| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 两性夫妻黄色片| 亚洲精品一卡2卡三卡4卡5卡| 狠狠狠狠99中文字幕| 成人一区二区视频在线观看| 床上黄色一级片| 制服人妻中文乱码| 色噜噜av男人的天堂激情| 成人鲁丝片一二三区免费| 波多野结衣高清无吗| 别揉我奶头~嗯~啊~动态视频| 国产亚洲av高清不卡| 天天躁日日操中文字幕| av福利片在线观看| 欧美成人性av电影在线观看| 欧美激情久久久久久爽电影| 亚洲国产精品成人综合色| 国产精品亚洲美女久久久| 噜噜噜噜噜久久久久久91| 国产欧美日韩精品一区二区| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| 国产精品一区二区免费欧美| 亚洲18禁久久av| 国产高清videossex| 日韩 欧美 亚洲 中文字幕| 狂野欧美激情性xxxx| 村上凉子中文字幕在线| 最新在线观看一区二区三区| 美女大奶头视频| xxxwww97欧美| 国产三级在线视频| 性欧美人与动物交配| 国产精品久久久久久亚洲av鲁大| 三级毛片av免费| 亚洲 国产 在线| 亚洲成人久久爱视频| 久久久久久人人人人人| 久99久视频精品免费| 国产高清视频在线播放一区| 黄色片一级片一级黄色片| 12—13女人毛片做爰片一| 欧美色视频一区免费| 999久久久精品免费观看国产| 日韩高清综合在线| av天堂中文字幕网| 男人的好看免费观看在线视频| 一二三四在线观看免费中文在| 国产成人精品久久二区二区91| 欧美乱色亚洲激情| 噜噜噜噜噜久久久久久91| 精品久久久久久久久久免费视频| 精华霜和精华液先用哪个| 日本精品一区二区三区蜜桃| 欧美午夜高清在线| 特大巨黑吊av在线直播| 国产黄片美女视频| 嫩草影院精品99| 色综合欧美亚洲国产小说| 动漫黄色视频在线观看| 18禁美女被吸乳视频| 啦啦啦观看免费观看视频高清| 他把我摸到了高潮在线观看| 亚洲av电影不卡..在线观看| 免费看a级黄色片| 久久这里只有精品19| 色哟哟哟哟哟哟| 午夜免费激情av| 国产精品一区二区三区四区免费观看 | 日本a在线网址| 天堂影院成人在线观看| 欧美午夜高清在线| 午夜精品一区二区三区免费看| 亚洲精品美女久久久久99蜜臀| 啦啦啦观看免费观看视频高清| 九色成人免费人妻av| 国产激情欧美一区二区| 欧美性猛交黑人性爽| 国产亚洲精品综合一区在线观看| 欧美日韩乱码在线| h日本视频在线播放| 亚洲国产欧美一区二区综合| 亚洲国产高清在线一区二区三| 香蕉国产在线看| 一区二区三区高清视频在线| 国产日本99.免费观看| 久久人人精品亚洲av| 99久国产av精品| 精品国产三级普通话版| 最近最新中文字幕大全电影3| 精品乱码久久久久久99久播| 亚洲av成人精品一区久久| 亚洲熟妇熟女久久| 久久久久久久久中文| 亚洲狠狠婷婷综合久久图片| 成人高潮视频无遮挡免费网站| 免费搜索国产男女视频| 久久精品91蜜桃| 久久久色成人| 大型黄色视频在线免费观看| 又爽又黄无遮挡网站| 搡老岳熟女国产| 国产激情久久老熟女| 午夜免费观看网址| 成人欧美大片| 亚洲性夜色夜夜综合| 国产成+人综合+亚洲专区| 成人午夜高清在线视频| 国产亚洲精品久久久久久毛片| 免费看十八禁软件| 国产高清视频在线播放一区| 色视频www国产| 美女午夜性视频免费| 丰满人妻熟妇乱又伦精品不卡| 99热只有精品国产| 国产成人精品久久二区二区免费| 一夜夜www| 国产精品一区二区精品视频观看| tocl精华| 女人高潮潮喷娇喘18禁视频| 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 国产欧美日韩精品一区二区| 午夜亚洲福利在线播放| 日本成人三级电影网站| 午夜福利视频1000在线观看| 999精品在线视频| 可以在线观看的亚洲视频| 久久久久久九九精品二区国产| 亚洲成人久久爱视频| 日韩欧美国产一区二区入口| 桃色一区二区三区在线观看| 国产又色又爽无遮挡免费看| 高清在线国产一区| 亚洲精品美女久久久久99蜜臀| 国产激情偷乱视频一区二区| 欧美性猛交黑人性爽| 亚洲国产日韩欧美精品在线观看 | 国产三级在线视频| 国内精品久久久久久久电影| av国产免费在线观看| 法律面前人人平等表现在哪些方面| 草草在线视频免费看| 久久热在线av| 欧美日韩中文字幕国产精品一区二区三区| 91九色精品人成在线观看| 亚洲性夜色夜夜综合| 不卡av一区二区三区| 天天躁日日操中文字幕| 人妻夜夜爽99麻豆av| 欧美日韩一级在线毛片| 久久久久亚洲av毛片大全| 久久精品国产99精品国产亚洲性色| 在线观看一区二区三区| 国产人伦9x9x在线观看| 亚洲18禁久久av| 天堂av国产一区二区熟女人妻| av在线天堂中文字幕| 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片| 国产高潮美女av| 亚洲成人久久爱视频| 国产野战对白在线观看| 国产v大片淫在线免费观看| 久久99热这里只有精品18| 亚洲欧美日韩高清在线视频| 亚洲av日韩精品久久久久久密| 欧美乱色亚洲激情| 欧美成人一区二区免费高清观看 | 99re在线观看精品视频| 黄色丝袜av网址大全| 波多野结衣高清作品| 在线国产一区二区在线| 欧美黑人巨大hd| 亚洲成人久久性| 桃色一区二区三区在线观看| 欧美午夜高清在线| 国产亚洲精品综合一区在线观看| 色精品久久人妻99蜜桃| 91字幕亚洲| 亚洲avbb在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲一区二区三区不卡视频| 俺也久久电影网| 欧美又色又爽又黄视频| 亚洲男人的天堂狠狠| 亚洲精品美女久久av网站| 免费在线观看视频国产中文字幕亚洲| 丝袜人妻中文字幕| 精品人妻1区二区| 国产精品久久电影中文字幕| 亚洲精品国产精品久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看| 毛片女人毛片| 亚洲成人精品中文字幕电影| 国产在线精品亚洲第一网站| 制服丝袜大香蕉在线| 全区人妻精品视频| 最新在线观看一区二区三区| 黄片小视频在线播放| 人人妻人人看人人澡| 国产精品 国内视频| 99久久久亚洲精品蜜臀av| 亚洲一区高清亚洲精品| 岛国视频午夜一区免费看| 国产亚洲精品一区二区www| 最近最新免费中文字幕在线| 日韩欧美国产一区二区入口| 精品久久久久久久久久久久久| 国产精品九九99| 免费无遮挡裸体视频| 久99久视频精品免费| 男女下面进入的视频免费午夜| 成人一区二区视频在线观看| 91久久精品国产一区二区成人 | 97人妻精品一区二区三区麻豆| 男女那种视频在线观看| 制服丝袜大香蕉在线| 动漫黄色视频在线观看| 少妇的丰满在线观看| 亚洲熟妇中文字幕五十中出| 亚洲av成人一区二区三| 国产成人影院久久av| 2021天堂中文幕一二区在线观| 高清在线国产一区| 极品教师在线免费播放| 久久精品91蜜桃| 老司机午夜福利在线观看视频| 中文字幕人妻丝袜一区二区| 国产免费男女视频| 1024手机看黄色片| 男女下面进入的视频免费午夜| 国产av麻豆久久久久久久| 99久久无色码亚洲精品果冻| 日韩免费av在线播放| 亚洲avbb在线观看| 欧美性猛交黑人性爽| 日本成人三级电影网站| 色视频www国产| 99精品欧美一区二区三区四区| 国产成+人综合+亚洲专区| 香蕉国产在线看| 国产一区在线观看成人免费| 99精品在免费线老司机午夜| 天天躁日日操中文字幕| av天堂在线播放| 国产在线精品亚洲第一网站| а√天堂www在线а√下载| 欧美日韩黄片免| 成人亚洲精品av一区二区| 国产精品日韩av在线免费观看| 亚洲欧美精品综合一区二区三区| 国产成人一区二区三区免费视频网站| 91麻豆精品激情在线观看国产| 国产精品亚洲美女久久久| 久久久色成人| 久久亚洲精品不卡| 亚洲国产日韩欧美精品在线观看 | 日韩大尺度精品在线看网址| 免费人成视频x8x8入口观看| 国产精品久久久久久亚洲av鲁大| 九色国产91popny在线| 91在线精品国自产拍蜜月 | 久久久久久国产a免费观看| 一二三四在线观看免费中文在| 国产成人一区二区三区免费视频网站| 18禁国产床啪视频网站| 国产成+人综合+亚洲专区| 国产淫片久久久久久久久 | 一级黄色大片毛片| 国产精品亚洲av一区麻豆| 亚洲熟妇熟女久久| av天堂中文字幕网| 久久精品国产99精品国产亚洲性色| 亚洲熟妇熟女久久| 久久国产精品影院| 少妇的丰满在线观看| 国产精品九九99| 国产精品亚洲一级av第二区| 精品久久蜜臀av无| 国产一区二区在线av高清观看| 欧美日韩综合久久久久久 | 男人舔奶头视频| 亚洲欧洲精品一区二区精品久久久| 少妇裸体淫交视频免费看高清| 久久久久久人人人人人| 级片在线观看| 狠狠狠狠99中文字幕| 亚洲av成人不卡在线观看播放网| 国产在线精品亚洲第一网站| 色播亚洲综合网| 日本与韩国留学比较| 99热精品在线国产| 成人无遮挡网站| 18禁黄网站禁片免费观看直播| 亚洲中文日韩欧美视频| 久久九九热精品免费| 亚洲美女视频黄频|