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

    長方體密閉結(jié)構(gòu)內(nèi)爆炸沖擊波傳播與疊加分析模型

    2016-10-15 06:47:23楊亞東李向東王曉鳴
    兵工學(xué)報 2016年8期
    關(guān)鍵詞:爆源沖量測試點

    楊亞東,李向東,王曉鳴

    (1.重慶紅宇精密工業(yè)有限責(zé)任公司,重慶402760;2.南京理工大學(xué)智能彈藥技術(shù)國防重點學(xué)科實驗室,江蘇南京210094)

    長方體密閉結(jié)構(gòu)內(nèi)爆炸沖擊波傳播與疊加分析模型

    楊亞東1,2,李向東2,王曉鳴2

    (1.重慶紅宇精密工業(yè)有限責(zé)任公司,重慶402760;2.南京理工大學(xué)智能彈藥技術(shù)國防重點學(xué)科實驗室,江蘇南京210094)

    采用鏡像爆源構(gòu)建內(nèi)爆炸結(jié)構(gòu)壁面反射作用與非線性疊加處理多波耦合效應(yīng)相結(jié)合的方法,建立了長方體密閉結(jié)構(gòu)內(nèi)爆炸沖擊波傳播和疊加分析模型。該模型可以用于內(nèi)爆炸荷載的快速工程計算。算例分析表明,結(jié)構(gòu)壁面點上第一峰值超壓、脈沖寬度和比沖量的預(yù)測值與試驗結(jié)果的最大相對誤差分別為-10.5%、-19.2%和13.2%,模型預(yù)測內(nèi)爆炸沖擊波特征參量與試驗結(jié)果較為吻合,驗證了該分析模型的有效性。

    兵器科學(xué)與技術(shù);內(nèi)爆炸;密閉空間;爆炸荷載;沖擊波超壓;傳播;疊加

    0 引言

    發(fā)生于密閉或半密閉空間內(nèi)部的爆炸稱為內(nèi)爆炸。內(nèi)爆炸產(chǎn)生的高溫、高壓產(chǎn)物無法及時向外泄漏,沖擊波在結(jié)構(gòu)壁面上來回反射,各壁面反射波相互疊加使超壓峰值增大,作用時間加長,因此內(nèi)爆炸對結(jié)構(gòu)內(nèi)的人員和設(shè)備的毀傷作用更加巨大。由于內(nèi)爆炸彈藥能量利用率高,破壞效應(yīng)更強,使得內(nèi)爆炸相關(guān)研究成為國內(nèi)外兵器領(lǐng)域的研究熱點[1-6]。

    國內(nèi)外關(guān)于內(nèi)爆炸的研究主要集中在抗爆容器,其基本形狀多為球體、圓柱體或半球體[7],對于長方體結(jié)構(gòu)的研究較少。而精確制導(dǎo)武器攻擊地下掩體和彈藥庫、反艦導(dǎo)彈攻擊艦船、恐怖襲擊地鐵車站等都有一個共同的特點,爆炸都多發(fā)生于長方體結(jié)構(gòu)的內(nèi)部。由此可見,長方體結(jié)構(gòu)的內(nèi)爆炸研究更具有實際應(yīng)用價值。另外,由于內(nèi)爆炸沖擊波的耦合效應(yīng)顯著,對內(nèi)爆炸作用機(jī)理的認(rèn)識尚不完善,缺乏內(nèi)爆炸沖擊波傳播與疊加的理論計算模型。

    本文以沖擊波傳播理論為基礎(chǔ),用鏡像爆源結(jié)合非線性疊加原理描述內(nèi)爆炸的多波耦合效應(yīng),建立了長方體密閉結(jié)構(gòu)內(nèi)爆炸沖擊波傳播與疊加理論模型,并通過試驗對所建立的模型進(jìn)行了驗證。研究成果可以應(yīng)用于內(nèi)爆炸相關(guān)的工程快速計算。

    1 爆炸沖擊波的傳播與反射

    爆炸沖擊波入射峰值超壓的計算式[8]為

    沖擊波在壁面上的反射分為正反射、規(guī)則斜反射和馬赫反射3種類型。正反射是指入射角為0°的反射。正反射超壓的計算公式[8]為

    式中:p0為空氣初始壓力。

    規(guī)則斜反射是指入射角大于0°而小于規(guī)則反射極限角φ0c的反射??諝庵斜_擊波規(guī)則斜反射的超壓計算公式為

    由上式可知,對應(yīng)于給定的π和t0有兩個t2(t2-和t2+)和兩個Δpr值(Δpr-和Δpr+),實際上的反射對應(yīng)于較小的δ,即(4)式取負(fù)號[8-9]。

    φ0c可由(4)式在給定π和φ時只有一個實根的條件確定,即

    馬赫反射選用文獻(xiàn)[10]中的計算方法計算:

    沖擊波到達(dá)時間由(7)式計算[11]:

    式中:ca為聲速;rc為裝藥半徑。

    由于壁面反射作用在瞬間完成,未耦合前的單個反射沖擊波的到達(dá)時間和正壓持續(xù)時間僅與藥量和比例距離有關(guān),故反射沖擊波的到達(dá)時間和正壓持續(xù)時間的計算方法和上述方法相同。

    位置R處反射沖擊波超壓與時間的關(guān)系[8]為

    式中:a為控制衰減率的常數(shù)[12],

    由于負(fù)相峰值壓力遠(yuǎn)低于正相壓力,要精確測量和計算負(fù)相特性較為困難,故忽略負(fù)相的作用。聯(lián)立(1)式~(5)式,只要給出炸藥的裝藥量和爆心距離R,就可以得出R處的壓力時程曲線。

    沖擊波超壓的比沖量由(11)式計算:

    2 密閉空間內(nèi)爆炸沖擊波疊加原理

    2.1內(nèi)爆炸疊加原理及鏡像爆源法

    結(jié)構(gòu)壁面任意點處的內(nèi)爆炸荷載由爆源對測試點的直接作用和壁面反射沖擊波對測試點的作用共同組成。壁面反射沖擊波作用可看作爆源關(guān)于壁面對稱的鏡像爆源對測試點的作用,耦合爆源和所有鏡像爆源的作用即為測試點上總的內(nèi)爆炸荷載。爆源和鏡像爆源的耦合方法采用Needham提出的LAMB疊加原理[13],該原理按照質(zhì)量、動量和能量守恒定律建立,其表達(dá)如(12)式、(13)式和(14)式所示:

    將LAMB疊加原理應(yīng)用于密閉空間內(nèi)爆炸計算,其假設(shè)條件為:爆源產(chǎn)生的沖擊波和鏡像爆源產(chǎn)生的沖擊波作為各自的獨立波源分開計算,在傳播過程中相互不影響。

    圖1為長方體密閉結(jié)構(gòu)內(nèi)爆炸鏡像爆源分布及其沖擊波傳播示意圖。圖1中:S為爆源;Si為爆源關(guān)于墻i對稱的鏡像爆源;S24為鏡像爆源S2關(guān)于墻4對稱的二次鏡像爆源;hi為爆源到各壁面的距離,平行于頁面的墻5和墻6(圖1中未標(biāo)出),其處理方式和墻1~墻4的相同。壁面的一次反射以壁面1的反射為例,爆源S在壁面1上的A點反射后向P點入射等效為爆源關(guān)于壁面1對稱的一次鏡像爆源S1對P點的作用。壁面的二次反射以壁面4為例,鏡像爆源S2在壁面4上B點的反射沖擊波對P點的作用等效為鏡像爆源S2關(guān)于壁面4對稱的二次鏡像爆源S24對P點的作用。其余各面的一次反射、二次反射的處理方法和上述方法相同,三次反射以二次鏡像爆源進(jìn)行再鏡像處理。結(jié)構(gòu)壁面作用點的受力以反射壓力計算,去除自身壁面,共5個壁面起反射作用,因此每次反射共有5個鏡像爆源起作用。內(nèi)爆炸具有多峰值的特性,后續(xù)的內(nèi)爆炸荷載波峰衰減迅速[14]。本文計算僅考慮壁面的前3次反射對爆炸荷載的影響,考慮更多反射次數(shù)的意義不大。

    圖1 鏡像爆源分布及其沖擊波傳播示意圖Fig.1 Distribution of image burst and schematic diagram of shockwave propagation

    2.2內(nèi)爆炸荷載的計算方法

    以爆源直接作用于結(jié)構(gòu)壁面上點的荷載計算為例,計算示意圖如圖2所示,爆源S的位置坐標(biāo)為(xS,yS,zS),測試點的位置坐標(biāo)為(xP,yP,zP),坐標(biāo)原點設(shè)于受力壁面的左下角。爆炸荷載計算步驟如下:

    1)Dx、Dy和Dz的計算:

    2)爆源到測試點的入射距離R:

    3)爆源到測試點的比例距離Z:

    圖2 爆源直接作用于壁面的計算示意圖Fig.2 Schematic diagram of shockwave reflecting from the wall

    4)以(1)式和(17)式計算入射超壓。

    5)入射角φ的計算:

    6)當(dāng)入射角為0°的正反射時,以(2)式計算反射壓力;當(dāng)入射角不為0°時,以(5)式計算規(guī)則反射極限角,當(dāng)φ≤φ0c時,為規(guī)則斜反射,以(3)式計算反射壓力;當(dāng)φ>φ0c時,進(jìn)入馬赫反射區(qū),以(6)式計算反射壓力。

    7)以第6步所得的結(jié)果,結(jié)合(7)式~(10)式,可求得反射超壓的時程變化曲線。

    LAMB疊加方法結(jié)合鏡像爆源思想,可以確定壁面二次和三次反射沖擊波的起始位置,使長方體密閉結(jié)構(gòu)的二次鏡像爆源和三次鏡像爆源的確定有了理論依據(jù)。得到爆源、一次鏡像爆源、二次鏡像爆源和三次鏡像爆源等在作用點的爆炸荷載后,按(12)式~(14)式對它們進(jìn)行耦合疊加,所得結(jié)果即為作用點上的內(nèi)爆炸荷載,對耦合后的壓力時程曲線積分可得到比沖量時程曲線。

    3 試驗驗證與結(jié)果分析

    為了驗證分析模型的正確性,進(jìn)行了內(nèi)爆炸驗證試驗。試驗用鋼筋混凝土靶房,內(nèi)腔長480cm×寬480cm×高266cm,底板和側(cè)墻厚24cm,頂板厚15cm,墻體配筋率為0.5%,以C30混凝土澆筑。靶房開有寬50cm×高100cm的門洞,便于試驗人員出入,門洞配有寬60cm×高110cm×厚2cm的鋼門,鋼門安裝于靶房內(nèi)部,從內(nèi)往外關(guān),完全覆蓋住門洞以防止爆轟產(chǎn)物外泄。靶房整體形態(tài)如圖3(a)所示。裝藥位于密閉結(jié)構(gòu)中心位置,炸藥采用鑄裝梯恩梯炸藥,密度為1.6g/cm3.取后墻水平中線和結(jié)構(gòu)拐點f處的內(nèi)爆炸荷載計算結(jié)果與試驗結(jié)果進(jìn)行對比,后墻水平中線的測試點P1~P5的位置和拐角f如圖3(b)所示。各測試點預(yù)埋M61×2內(nèi)螺紋管,試驗時在螺紋管中安裝BA-YD205壓力傳感器,用NLG-202G數(shù)據(jù)采集系統(tǒng)采集數(shù)據(jù)。

    圖3 試驗靶房與壁面測試點分布示意圖Fig.3 Schematic diagram of reinforced concrete target room and distribution of test points

    圖4 720g梯恩梯裝藥在壁面中線的超壓峰值對比Fig.4 Comparison of shockwave overpressures of 720 g TNT explosive in the midline of wall

    3.1結(jié)構(gòu)壁面直線上的內(nèi)爆炸荷載結(jié)果分析

    720g裝藥作用下后墻壁面中線上的反射壓力如圖4所示。由圖4可以看出,壁面壓力關(guān)于中點左右對稱,呈W形分布。W形分布曲線的BC段受爆源直接作用的影響較大,各壁面的反射壓力對BC段的影響小于爆源的直接作用。隨著從中心P3點向兩端B、C點移動,作用點的入射距離R逐漸增大,反射壓力逐漸減小,使得BC段壓力呈凸形。當(dāng)達(dá)到B或C點后,爆源的影響不起主導(dǎo)作用,來自于其他壁面的反射沖擊波的耦合作用顯著加強,使得AB和CD段的壓力逐漸呈增長趨勢。BC段上理論與試驗的最大相對誤差在P4點產(chǎn)生,其最大相對誤差為6.7%.在W形分布的兩翼AB和CD段,理論與試驗的相對誤差隨著與鄰近側(cè)墻距離的減小而逐漸增大,其最大相對誤差為10.5%.

    圖5 壁面測試點壓力和比沖量時程曲線的理論與試驗值對比(720g梯恩梯炸藥)Fig.5 Comparison of experimental and theoretical pressures and impulse-time histories on the wall(720 g TNT)

    表1 測試點沖擊波超壓、脈沖寬度和比沖量(720g梯恩梯炸藥)Tab.1 Shockwave overpressures,pulse widths and impulses at different test points(720 g TNT)

    圖5為P1~P5點壓力和比沖量時程曲線的理論與試驗值對比。表1給出了峰值超壓和比沖量的對比情況。從圖5可以看出,處于W形分布曲線BC范圍內(nèi)的P3和P4點由爆源直接作用產(chǎn)生的波峰與壁面一次反射所產(chǎn)生的波峰彼此分開,即第一波峰由爆源所產(chǎn)生,壁面反射壓力在作用點上產(chǎn)生的波峰落后于爆源的波峰。P3和P4點第一峰值超壓的理論與試驗相對誤差為-6.0%和6.7%,脈沖寬度相對誤差分別為3.8%和-19.2%,比沖量相對誤差分別為8.6%和7.3%.由P1、P2和P5點的壓力時程曲線可以看出,當(dāng)處于AB和CD段時,鄰近壁面的一次反射波趕上爆源的波峰,使爆源和鄰近壁面的反射波相互疊加,第一波峰由單峰變?yōu)殇忼X狀的多峰,第一峰值由爆源和鄰近壁面的一次反射波共同作用,最大峰值由疊加后的齒狀多峰的最高值決定。P1、P2和P5點第一峰值超壓的理論與試驗相對誤差為-10.5%、-7.9%和5.8%,脈沖寬度相對誤差分別為12.1%、-12.4%和-17.4%,比沖量相對誤差分別為13.2%、9.7%和6.5%.由于理論沒有考慮負(fù)向超壓的影響,理論計算的比沖量時程曲線略高于試驗,比沖量的誤差略大于超壓峰值的誤差。由P1~P5的總體對比可以看出,理論和試驗所得的壓力和比沖量時程曲線吻合較好。

    3.2結(jié)構(gòu)拐角處的內(nèi)爆炸荷載結(jié)果分析

    結(jié)構(gòu)拐角處的內(nèi)爆炸增強效應(yīng)明顯,有必要對結(jié)構(gòu)上比較特殊的受力位置進(jìn)行單獨驗證。圖6為結(jié)構(gòu)拐角f在不同藥量作用下第一峰值超壓的理論與試驗對比。從圖6中可以看出,274.75g、322.5g、483.5g和636.5g梯恩梯作用下結(jié)構(gòu)拐角處的第一峰值超壓理論與試驗的相對誤差分別為9.6%、8.1%、10.3%和-2.4%,第一峰值超壓的理論值和試驗結(jié)果吻合良好。

    圖6 不同藥量作用下結(jié)構(gòu)拐角處第一峰值超壓的理論與試驗對比Fig.6 Comparison of shockwave overpressures at the corner for experiment and theoretical charge masses

    圖7 拐角處壓力和比沖量時程曲線的理論與試驗對比(483.5g梯恩梯)Fig.7 Comparison of experiment and theoretical pressures and impulse time-histories at the corner(483.5 g TNT)

    圖7為483.5g梯恩梯炸藥作用下結(jié)構(gòu)拐角f處壓力和比沖量時程曲線的理論與試驗值對比。第一峰值超壓、脈沖寬度和比沖量的理論預(yù)測值與試驗結(jié)果的最大相對誤差分別為10.3%、13.7%和11.4%,理論計算的壓力和比沖量時程曲線與試驗結(jié)果較為吻合。

    通過對長方體密閉結(jié)構(gòu)壁面直線和結(jié)構(gòu)拐角的內(nèi)爆炸荷載計算分析可以看出,本文模型有如下優(yōu)點:1)可以計算長方體結(jié)構(gòu)壁面任意位置處的內(nèi)爆炸荷載,能夠準(zhǔn)確表征內(nèi)爆炸的多峰值特性,能較好地描述內(nèi)爆炸的爆源沖擊波和壁面反射沖擊波之間的相互追趕疊加關(guān)系;2)考慮了入射角度和反射類型的關(guān)系;3)能夠處理結(jié)構(gòu)壁面上沖擊波的多次反射作用,能夠很好地計算內(nèi)爆炸在結(jié)構(gòu)棱邊、角點處的壓力增強效應(yīng)。

    4 結(jié)論

    1)以鏡像爆源結(jié)合非線性疊加原理構(gòu)建了長方體結(jié)構(gòu)內(nèi)爆炸沖擊波的傳播與疊加模型,計算了結(jié)構(gòu)壁面直線和拐角處的爆炸荷載和壓力分布規(guī)律。計算表明,該模型能夠準(zhǔn)確地表征內(nèi)爆炸的多峰值特性,且能較好地描述內(nèi)爆炸的爆源沖擊波和壁面反射沖擊波之間的相互追趕疊加關(guān)系。

    2)結(jié)構(gòu)壁面直線上和結(jié)構(gòu)拐角的驗證試驗表明,第一峰值超壓的理論與試驗最大相對誤差為-10.5%,脈沖寬度最大相對誤差為-19.2%,比沖量最大相對誤差為13.2%.

    3)結(jié)構(gòu)壁面內(nèi)爆炸荷載的理論和試驗結(jié)果的誤差在工程允許范圍之內(nèi),所建立的內(nèi)爆炸沖擊波傳播與疊加計算模型可以應(yīng)用于內(nèi)爆炸相關(guān)的工程快速計算。

    (References)

    [1] Anderson C E,Baker W E,Wauters D K,et al.Quasi-static pressure,duration,and impulse for explosions e.g.HE in structures[J].International Journal of Mechanical Sciences,1983,25(6):455-464.

    [2] Scott D H,Gregory D K.Development,characterization and testing of high blast thermataric compositions[C]∥The 31th International Pyrotechnics Seminer:Fort Collins,Colorado:IPSUSA Seminars Inc,2004.

    [3] 夏志成,許多,王靜,等.結(jié)構(gòu)內(nèi)爆炸荷載作用下鋼筋鋼纖混凝土抗爆墻設(shè)計探討[J].工程爆破,2008,14(2):8-11. XIA Zhi-cheng,XU Duo,WANG Jing,et al.Study of blast wall design under explosion load from its internal structure[J].Engineering Blasting,2008,14(2):8-11.(in Chinese)

    [4] 胡宏偉,宋浦,趙省向,等.有限空間內(nèi)部爆炸研究進(jìn)展[J].含能材料,2013,21(4):539-546. HU Hong-wei,SONG Pu,ZHAO Sheng-xiang,et al.Progress in explosion in confined space[J].Chinese Journal of Energetic materials,2013,21(4):539-546.(in Chinese)

    [5] 楊亞東,李向東,王曉鳴,等.密閉空間內(nèi)爆炸縮比相似模型研究[J].振動與沖擊,2014,33(2):133-139. YANG Ya-dong,LI Xiang-dong,WANG Xiao-ming,et al.Scale similarity mode of internal explosion in closed field[J].Journal of Vibration and Shock,2014,33(2):133-139.(in Chinese)

    [6] 楊亞東,李向東,王曉鳴,等.鋼筋混凝土結(jié)構(gòu)內(nèi)爆炸相似模型試驗研究[J].南京理工大學(xué)學(xué)報,2016,40(2):135-141. YANG Ya-dong,LI Xiang-dong,WANG Xiao-ming,et al.Experimental study on similarity mode of reinforced concrete structure under internal explosion[J].Journal of Nanjing University of Science and Technology,2016,40(2):135-141.(in Chinese)

    [7] 郭志昆,宋鋒良,劉峰,等.扁平箱形密閉結(jié)構(gòu)內(nèi)爆炸的模型試驗[J].解放軍理工大學(xué)學(xué)報,2008,9(4):345-350. GOU Zhi-kun,SONG Feng-liang,LIU Feng,et al.Experiment of closed flat box structure subjected to internal detonation[J].Journal of PLA University of Science and Technology,2008,9(4): 345-350.(in Chinese)

    [8] 北京工業(yè)學(xué)院八系《爆炸及其作用》編寫組.爆炸及其作用[M].北京:國防工業(yè)出版社,1979.“Explosion and Its Use”Editing Group,the Eight Department,Beijing Institute of Technology.Explosion and its use[M].Beijing:National Defense Industry Press,1979.(in Chinese)

    [9] 奧爾連科л Π.爆炸物理學(xué)[M].孫承緯,譯.北京:科學(xué)出版社,2011. OΡлeΗκο л Π.Explosion physics[M].SUN Cheng-wei,translated.Beijing:Science Press,2011.(in Chinese)

    [10] 易仰賢.空爆沖擊波馬赫反射近似計算[J].爆炸與沖擊,1983,3(2):44-49. YI Yang-xian.Approximate calculation of mach reflection of explosive shock waves in air[J].Explosion and Shock Waves,1983,3(2):44-49.(in Chinese)

    [11] Kinney G F,Graham K J.Explosive shocks in air[M].Berlin: Springer-Verlag,1985.

    [12] Wu C,Hao H.Modeling of simultaneous ground shock and airblast pressure on nearby structures from surface explosions[J]. International Journal of Impact Engineering,2005,31(6):699-717.

    [13] Needham C E.Blast waves[M].Berlin:Springer-Verlag,2010.

    [14] Baker W E.Explosion hazards and evaluation[M].New York: Elsevier Publishing Company,1983.

    An Analytical Model for Propagation and Superposition of Internal Explosion Shockwaves in Closed Cuboid Structure

    YANG Ya-dong1,2,LI Xiang-dong2,WANG Xiao-ming2
    (1.Chongqing Hongyu Precision Industrial Co.,Ltd,Chongqing 402760,China;2.Ministerial Key Laboratory of ZNDY,Nanjing University of Science and Technology,Nanjing 210094,Jiangsu,China)

    The method of using an image blast source to establish the reflection effect from the walls of structure is combined with the nonlinear superposition principle to deal with the multi-wave coupling effect.An analytical model of shockwave propagation and superposition in closed cuboid structure is established.The proposed model can be used for the rapid calculation of internal explosion loading in engineering projects.The results show that the maximum relative errors of theoretical and experimental first peak overpressures,pulse widths and specific impulses on the wall are less than-10.5%,-19.2% and 13.2%,respectively.The theoretical characteristic parameters of internal explosion are in good agreement with the experimental ones,which verifies the validity of this analytical model.

    ordnance science and technology;internal explosion;closed space;explosion loading;shockwave overpressure;propagation;superposition

    O383+.2

    A

    1000-1093(2016)08-1449-07

    10.3969/j.issn.1000-1093.2016.08.016

    2015-08-15

    國家自然科學(xué)基金項目(51278250)

    楊亞東(1979—),男,工程師,博士。E-mail:luxiya_2000_2000@tom.com;李向東(1969—),男,教授。E-mail:lixiangd@njust.edu.cn

    猜你喜歡
    爆源沖量測試點
    一種新型模擬電路故障字典測點選擇方法研究
    爆源和測點深度對水下爆炸沖擊波載荷的影響
    基于信息熵可信度的測試點選擇方法研究
    連續(xù)3年銷量翻番,2022年欲沖量4000萬,福建這家動保企業(yè)正強勢崛起
    邏輯內(nèi)建自測試雙重過濾測試點選取策略
    中型浮動沖擊平臺試驗爆源定位方法研究
    中國測試(2018年10期)2018-11-17 01:58:24
    不同攻角對實船爆炸試驗沖擊響應(yīng)影響研究
    某型艦抗沖擊試驗爆源定位方法研究*
    測試點的優(yōu)化選擇
    電子測試(2012年2期)2012-08-07 07:52:52
    變力沖量的教學(xué)研究
    物理通報(2011年8期)2011-01-24 07:39:38
    国语对白做爰xxxⅹ性视频网站| 1000部很黄的大片| 在线免费观看不下载黄p国产| 汤姆久久久久久久影院中文字幕 | 九九爱精品视频在线观看| 综合色丁香网| 婷婷色综合大香蕉| 国产一区二区三区综合在线观看 | 亚洲欧美日韩无卡精品| 男女那种视频在线观看| 插阴视频在线观看视频| 最近中文字幕2019免费版| 晚上一个人看的免费电影| 男女国产视频网站| 亚洲欧美日韩无卡精品| 中文资源天堂在线| 久久久久久久久久久免费av| 精品久久久久久久久亚洲| 床上黄色一级片| 国产一区亚洲一区在线观看| 亚洲精品日韩在线中文字幕| 国产午夜精品论理片| 国产在线一区二区三区精| 亚洲精品第二区| 精品午夜福利在线看| 亚洲精品日韩av片在线观看| 亚洲综合色惰| 春色校园在线视频观看| 午夜精品国产一区二区电影 | 日本一本二区三区精品| 欧美区成人在线视频| 又大又黄又爽视频免费| 精品久久久久久久末码| 日本三级黄在线观看| 久久亚洲国产成人精品v| 伦精品一区二区三区| 九九爱精品视频在线观看| 亚洲久久久久久中文字幕| 欧美精品国产亚洲| 日本黄大片高清| 天堂网av新在线| 久久人人爽人人片av| 一区二区三区乱码不卡18| 九九爱精品视频在线观看| 色5月婷婷丁香| 黄色配什么色好看| 在线免费观看不下载黄p国产| 国产真实伦视频高清在线观看| 国产精品综合久久久久久久免费| 搞女人的毛片| 日本黄色片子视频| videossex国产| 久久久久久久久久成人| 成人性生交大片免费视频hd| 国产视频首页在线观看| 亚洲熟妇中文字幕五十中出| 国内揄拍国产精品人妻在线| 国产av码专区亚洲av| 天堂av国产一区二区熟女人妻| 三级毛片av免费| 国产日韩欧美在线精品| 日韩av在线免费看完整版不卡| 777米奇影视久久| 美女内射精品一级片tv| 99九九线精品视频在线观看视频| 国产色爽女视频免费观看| 亚洲欧美一区二区三区国产| 嘟嘟电影网在线观看| 特级一级黄色大片| av播播在线观看一区| 校园人妻丝袜中文字幕| av免费观看日本| 黄色配什么色好看| 国产精品人妻久久久久久| 中文字幕久久专区| 99热这里只有精品一区| 亚洲欧美成人精品一区二区| 欧美区成人在线视频| 中文字幕人妻熟人妻熟丝袜美| 爱豆传媒免费全集在线观看| 男女啪啪激烈高潮av片| 丰满人妻一区二区三区视频av| 一个人观看的视频www高清免费观看| 欧美日韩亚洲高清精品| 国产成人午夜福利电影在线观看| 免费看光身美女| 久久久久精品性色| 青春草国产在线视频| 午夜福利成人在线免费观看| 欧美极品一区二区三区四区| 97精品久久久久久久久久精品| 久久久精品欧美日韩精品| 男女下面进入的视频免费午夜| 日韩av在线免费看完整版不卡| 精品国产露脸久久av麻豆 | a级毛色黄片| 国产精品美女特级片免费视频播放器| 深夜a级毛片| 久久久久久久久久人人人人人人| 尤物成人国产欧美一区二区三区| 97精品久久久久久久久久精品| 亚洲国产精品成人综合色| 欧美成人精品欧美一级黄| 精品一区二区三卡| 久久久久网色| 亚洲成人精品中文字幕电影| 天天躁夜夜躁狠狠久久av| 国产欧美另类精品又又久久亚洲欧美| 久久人人爽人人爽人人片va| 波野结衣二区三区在线| 亚洲色图av天堂| 三级男女做爰猛烈吃奶摸视频| 99热全是精品| 国内精品宾馆在线| 精品久久久久久久久久久久久| 网址你懂的国产日韩在线| 菩萨蛮人人尽说江南好唐韦庄| 搡老妇女老女人老熟妇| 婷婷色综合www| 天美传媒精品一区二区| 舔av片在线| 久久久久久久久久成人| 欧美成人午夜免费资源| 国产一级毛片在线| 亚洲电影在线观看av| 日日撸夜夜添| 毛片女人毛片| 亚洲精品456在线播放app| 男人狂女人下面高潮的视频| 国产片特级美女逼逼视频| 久久人人爽人人爽人人片va| 亚洲综合色惰| 亚洲国产欧美在线一区| 少妇的逼好多水| 男人和女人高潮做爰伦理| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 成年版毛片免费区| 欧美极品一区二区三区四区| 深爱激情五月婷婷| 97在线视频观看| 国内精品宾馆在线| 国产人妻一区二区三区在| 久久综合国产亚洲精品| 国产成人福利小说| 久久这里只有精品中国| 国内精品美女久久久久久| 一区二区三区高清视频在线| 搡老乐熟女国产| 久热久热在线精品观看| 国产探花在线观看一区二区| 亚洲三级黄色毛片| 99热全是精品| 国语对白做爰xxxⅹ性视频网站| 大片免费播放器 马上看| 舔av片在线| 欧美另类一区| 免费在线观看成人毛片| 51国产日韩欧美| 久久久久九九精品影院| 欧美一区二区亚洲| 夫妻午夜视频| 久久99热这里只频精品6学生| 国产三级在线视频| 在线免费观看的www视频| 精华霜和精华液先用哪个| 国产女主播在线喷水免费视频网站 | 欧美成人精品欧美一级黄| 三级男女做爰猛烈吃奶摸视频| 在线免费观看的www视频| 亚洲人与动物交配视频| 黄色配什么色好看| 欧美日韩在线观看h| av又黄又爽大尺度在线免费看| 久久国产乱子免费精品| 亚洲精品色激情综合| 亚洲欧洲日产国产| 狂野欧美激情性xxxx在线观看| 亚洲18禁久久av| xxx大片免费视频| 久久久久性生活片| av福利片在线观看| 自拍偷自拍亚洲精品老妇| 身体一侧抽搐| 免费无遮挡裸体视频| 日韩成人伦理影院| 久久久久久久亚洲中文字幕| 色综合色国产| 老司机影院毛片| 成人国产麻豆网| 99热网站在线观看| 亚洲av成人精品一二三区| 亚州av有码| 中文字幕av电影在线播放| av在线老鸭窝| 成人手机av| 人妻人人澡人人爽人人| 人妻 亚洲 视频| 亚洲欧美一区二区三区黑人 | 极品人妻少妇av视频| 成人漫画全彩无遮挡| 日本爱情动作片www.在线观看| av国产精品久久久久影院| 午夜激情久久久久久久| 一本久久精品| 少妇的逼水好多| 少妇人妻 视频| 国产精品久久久久久精品电影小说| 国产精品香港三级国产av潘金莲 | 人成视频在线观看免费观看| 国产在线视频一区二区| 男女下面插进去视频免费观看| 麻豆av在线久日| a级毛片在线看网站| 又大又黄又爽视频免费| 欧美精品一区二区大全| 午夜福利在线免费观看网站| 美女大奶头黄色视频| 啦啦啦中文免费视频观看日本| 久久人人97超碰香蕉20202| www.av在线官网国产| 婷婷色麻豆天堂久久| 久久这里只有精品19| 国产国语露脸激情在线看| 天堂8中文在线网| 丰满少妇做爰视频| 亚洲精品在线美女| 精品一品国产午夜福利视频| 男的添女的下面高潮视频| 一区二区三区精品91| 中文欧美无线码| 成人毛片a级毛片在线播放| 欧美日韩一级在线毛片| 久久热在线av| 欧美 亚洲 国产 日韩一| 国产日韩欧美视频二区| 亚洲精品美女久久久久99蜜臀 | 亚洲国产欧美网| 亚洲一区中文字幕在线| 蜜桃国产av成人99| 一级爰片在线观看| 国产一区二区三区综合在线观看| 亚洲第一区二区三区不卡| 国产精品 国内视频| 亚洲av.av天堂| 国产亚洲午夜精品一区二区久久| 色视频在线一区二区三区| 97在线视频观看| 国产激情久久老熟女| 一边亲一边摸免费视频| 国产成人免费无遮挡视频| 日韩 亚洲 欧美在线| 熟妇人妻不卡中文字幕| 热99久久久久精品小说推荐| 深夜精品福利| 亚洲,欧美,日韩| 一区二区日韩欧美中文字幕| 国产不卡av网站在线观看| 寂寞人妻少妇视频99o| 丝袜人妻中文字幕| 国产国语露脸激情在线看| 国产精品欧美亚洲77777| 久久影院123| av国产精品久久久久影院| 在线观看www视频免费| 国产日韩欧美亚洲二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美日韩亚洲高清精品| 免费黄频网站在线观看国产| 蜜桃在线观看..| 五月伊人婷婷丁香| 中文乱码字字幕精品一区二区三区| 在线免费观看不下载黄p国产| av在线app专区| 在线天堂最新版资源| av片东京热男人的天堂| 国产1区2区3区精品| 精品少妇久久久久久888优播| 久久午夜综合久久蜜桃| 激情视频va一区二区三区| h视频一区二区三区| 男女下面插进去视频免费观看| 国产一区二区在线观看av| 色播在线永久视频| 国产精品蜜桃在线观看| 精品视频人人做人人爽| 黑丝袜美女国产一区| 18在线观看网站| 国产成人91sexporn| 久久久精品94久久精品| 9热在线视频观看99| 久久久久久久亚洲中文字幕| av在线老鸭窝| 午夜福利视频在线观看免费| 午夜免费男女啪啪视频观看| videosex国产| 成人国语在线视频| 不卡av一区二区三区| 两个人免费观看高清视频| 久久久欧美国产精品| 欧美最新免费一区二区三区| 亚洲一区二区三区欧美精品| 免费播放大片免费观看视频在线观看| 国产极品天堂在线| 18禁裸乳无遮挡动漫免费视频| 777久久人妻少妇嫩草av网站| 亚洲成av片中文字幕在线观看 | 亚洲精品国产色婷婷电影| 亚洲av免费高清在线观看| 免费黄网站久久成人精品| 精品视频人人做人人爽| 免费女性裸体啪啪无遮挡网站| 欧美精品一区二区大全| 色婷婷av一区二区三区视频| 久久久久久久久免费视频了| 亚洲成人一二三区av| www.av在线官网国产| 国产成人免费观看mmmm| 国产在视频线精品| 韩国高清视频一区二区三区| 久久精品国产a三级三级三级| 国产免费又黄又爽又色| 亚洲精品美女久久av网站| 久久精品久久久久久噜噜老黄| 国语对白做爰xxxⅹ性视频网站| 男人添女人高潮全过程视频| 纵有疾风起免费观看全集完整版| 精品国产乱码久久久久久小说| 精品一区二区三区四区五区乱码 | 久久青草综合色| 亚洲国产欧美日韩在线播放| 国产成人午夜福利电影在线观看| 天天躁夜夜躁狠狠久久av| 91精品伊人久久大香线蕉| 亚洲,欧美精品.| 视频区图区小说| 自拍欧美九色日韩亚洲蝌蚪91| xxxhd国产人妻xxx| 老司机亚洲免费影院| 极品少妇高潮喷水抽搐| 国产精品一区二区在线观看99| 日日撸夜夜添| 女的被弄到高潮叫床怎么办| 男人操女人黄网站| 另类精品久久| 婷婷色综合大香蕉| 国产精品久久久久久久久免| 欧美精品一区二区免费开放| 黄色 视频免费看| 如日韩欧美国产精品一区二区三区| 国产精品99久久99久久久不卡 | 男女啪啪激烈高潮av片| 日本爱情动作片www.在线观看| 国产精品久久久久久精品电影小说| 女人被躁到高潮嗷嗷叫费观| 国产精品二区激情视频| 69精品国产乱码久久久| 久久精品国产综合久久久| 国产又爽黄色视频| 交换朋友夫妻互换小说| 校园人妻丝袜中文字幕| 国产精品蜜桃在线观看| 国产国语露脸激情在线看| av不卡在线播放| 黄片小视频在线播放| 午夜福利网站1000一区二区三区| 久久久久久久久久久免费av| 青草久久国产| 精品少妇内射三级| 久久久久精品人妻al黑| 日日摸夜夜添夜夜爱| 美女中出高潮动态图| 久久久精品国产亚洲av高清涩受| 我的亚洲天堂| 久久精品夜色国产| av在线app专区| 欧美精品一区二区大全| 午夜久久久在线观看| 久久精品久久久久久久性| kizo精华| 大香蕉久久成人网| 国产日韩欧美视频二区| 高清黄色对白视频在线免费看| www.精华液| 咕卡用的链子| 亚洲av在线观看美女高潮| 国产av码专区亚洲av| 亚洲国产av新网站| 久久久久久久大尺度免费视频| 一区二区三区精品91| 下体分泌物呈黄色| 精品少妇黑人巨大在线播放| 国产成人91sexporn| 啦啦啦中文免费视频观看日本| 久久99热这里只频精品6学生| 少妇熟女欧美另类| kizo精华| 亚洲美女黄色视频免费看| 久久午夜综合久久蜜桃| 又粗又硬又长又爽又黄的视频| 日韩制服丝袜自拍偷拍| 亚洲成国产人片在线观看| tube8黄色片| 亚洲综合色网址| 久久人人97超碰香蕉20202| 久久99精品国语久久久| 亚洲国产精品一区三区| 欧美日韩av久久| 亚洲久久久国产精品| 亚洲色图 男人天堂 中文字幕| 国产成人精品婷婷| av在线app专区| 亚洲av男天堂| 国产精品久久久av美女十八| 久久久久视频综合| 国产xxxxx性猛交| 欧美国产精品一级二级三级| 欧美精品高潮呻吟av久久| 久久精品国产自在天天线| av国产久精品久网站免费入址| 少妇人妻精品综合一区二区| 国产在线视频一区二区| 久久精品久久久久久噜噜老黄| 日本av手机在线免费观看| 一二三四在线观看免费中文在| 国产精品麻豆人妻色哟哟久久| 国产精品av久久久久免费| 久久久久久久久免费视频了| 韩国精品一区二区三区| 免费大片黄手机在线观看| 可以免费在线观看a视频的电影网站 | 少妇的丰满在线观看| 美女午夜性视频免费| 久久久精品94久久精品| 亚洲精华国产精华液的使用体验| 亚洲精品日韩在线中文字幕| 久久久久国产一级毛片高清牌| 人妻系列 视频| 午夜福利在线免费观看网站| 777米奇影视久久| 老汉色∧v一级毛片| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 国产又色又爽无遮挡免| 成人国语在线视频| 赤兔流量卡办理| 久久午夜综合久久蜜桃| 女人久久www免费人成看片| 国产在视频线精品| 视频在线观看一区二区三区| 丝袜在线中文字幕| 99久久中文字幕三级久久日本| 亚洲国产毛片av蜜桃av| 日本vs欧美在线观看视频| 欧美 日韩 精品 国产| 精品国产露脸久久av麻豆| 亚洲精品久久久久久婷婷小说| 91久久精品国产一区二区三区| 亚洲美女搞黄在线观看| 久久99蜜桃精品久久| 日韩一区二区三区影片| 日韩一本色道免费dvd| 我要看黄色一级片免费的| 90打野战视频偷拍视频| 欧美日韩亚洲国产一区二区在线观看 | av视频免费观看在线观看| 日韩av在线免费看完整版不卡| 高清不卡的av网站| 高清视频免费观看一区二区| 一边亲一边摸免费视频| 国产亚洲av片在线观看秒播厂| 亚洲国产欧美在线一区| 赤兔流量卡办理| 青青草视频在线视频观看| 一级爰片在线观看| 最近最新中文字幕免费大全7| 亚洲国产色片| 看免费成人av毛片| 大码成人一级视频| 亚洲美女搞黄在线观看| 最近最新中文字幕大全免费视频 | 一区二区三区激情视频| 国产有黄有色有爽视频| 色播在线永久视频| 大香蕉久久网| 亚洲综合色网址| 青春草国产在线视频| 久久人妻熟女aⅴ| 国产免费视频播放在线视频| 亚洲国产精品999| 男女边吃奶边做爰视频| 亚洲国产欧美日韩在线播放| 欧美日韩视频精品一区| 91成人精品电影| 又粗又硬又长又爽又黄的视频| 午夜福利视频在线观看免费| 国产成人免费观看mmmm| 久久久久久人人人人人| 亚洲av中文av极速乱| 高清不卡的av网站| 久久毛片免费看一区二区三区| 日韩三级伦理在线观看| 在线 av 中文字幕| 成人免费观看视频高清| 少妇猛男粗大的猛烈进出视频| 亚洲欧美成人精品一区二区| 亚洲一级一片aⅴ在线观看| 伊人久久国产一区二区| 男女边吃奶边做爰视频| 99re6热这里在线精品视频| 国产在视频线精品| 大香蕉久久成人网| 91成人精品电影| 国产精品一区二区在线不卡| 国产男女超爽视频在线观看| 侵犯人妻中文字幕一二三四区| 日本午夜av视频| 久久精品国产亚洲av天美| 看非洲黑人一级黄片| 亚洲精品在线美女| 色网站视频免费| 国产精品久久久久久av不卡| 久久av网站| 男女国产视频网站| 国产精品香港三级国产av潘金莲 | 天堂俺去俺来也www色官网| 美女大奶头黄色视频| 人妻少妇偷人精品九色| 国产 一区精品| 亚洲精品国产av蜜桃| 一边亲一边摸免费视频| 亚洲在久久综合| 人妻人人澡人人爽人人| 大香蕉久久成人网| 飞空精品影院首页| 国产极品粉嫩免费观看在线| 晚上一个人看的免费电影| 国产成人精品久久二区二区91 | 一区在线观看完整版| 欧美成人精品欧美一级黄| 亚洲国产精品一区二区三区在线| 色播在线永久视频| videosex国产| 欧美日韩国产mv在线观看视频| 最近最新中文字幕免费大全7| 国产无遮挡羞羞视频在线观看| 最新中文字幕久久久久| 国产精品香港三级国产av潘金莲 | 各种免费的搞黄视频| 亚洲成av片中文字幕在线观看 | 日韩精品免费视频一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜制服| 欧美国产精品va在线观看不卡| 亚洲男人天堂网一区| 国产精品 欧美亚洲| 免费看不卡的av| 国产成人精品婷婷| 十八禁高潮呻吟视频| 亚洲美女视频黄频| 国产野战对白在线观看| 黑人猛操日本美女一级片| 午夜福利在线免费观看网站| 三级国产精品片| 少妇的丰满在线观看| 亚洲国产日韩一区二区| 欧美精品亚洲一区二区| 亚洲国产av新网站| 欧美黄色片欧美黄色片| 老司机影院成人| 亚洲精品日韩在线中文字幕| 啦啦啦啦在线视频资源| 欧美亚洲 丝袜 人妻 在线| 少妇熟女欧美另类| 91午夜精品亚洲一区二区三区| 老司机亚洲免费影院| 中文字幕人妻丝袜制服| 亚洲综合精品二区| 国产有黄有色有爽视频| 亚洲综合精品二区| 大陆偷拍与自拍| 女性被躁到高潮视频| 成年美女黄网站色视频大全免费| 国产免费一区二区三区四区乱码| 亚洲av电影在线进入| 高清欧美精品videossex| 免费女性裸体啪啪无遮挡网站| 午夜91福利影院| 欧美日韩一区二区视频在线观看视频在线| 一本大道久久a久久精品| 成人毛片a级毛片在线播放| 老司机影院成人| 永久免费av网站大全| 亚洲精品在线美女| 亚洲国产av影院在线观看| 久久国产精品男人的天堂亚洲| 亚洲内射少妇av| 伦精品一区二区三区| 午夜免费鲁丝| 在线天堂中文资源库| 热re99久久国产66热| 国产成人精品久久二区二区91 | 中文欧美无线码| 青春草国产在线视频| 69精品国产乱码久久久| 精品亚洲乱码少妇综合久久| 日韩中文字幕视频在线看片| 欧美变态另类bdsm刘玥| 青春草亚洲视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产 精品1| av免费观看日本| 亚洲国产精品国产精品| 日本vs欧美在线观看视频|