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

    艙室密閉空間中爆炸載荷燃燒增強(qiáng)效應(yīng)試驗(yàn)研究

    2020-02-18 03:21:02孔祥韶況正鄭成吳衛(wèi)國(guó)
    兵工學(xué)報(bào) 2020年1期
    關(guān)鍵詞:氦氣藥量靜態(tài)

    孔祥韶,況正,鄭成,吳衛(wèi)國(guó)

    (1.武漢理工大學(xué) 綠色智能江海直達(dá)船舶與郵輪游艇研究中心,湖北 武漢 430063;2.武漢理工大學(xué) 交通學(xué)院 船舶、海洋與結(jié)構(gòu)工程系,湖北 武漢 430063)

    0 引言

    當(dāng)爆炸發(fā)生在艦船艙室密閉空間內(nèi)部時(shí),約束環(huán)境對(duì)爆炸載荷的影響甚為明顯,將造成結(jié)構(gòu)的嚴(yán)重毀傷[1]。梯恩梯(TNT)為典型的負(fù)氧型炸藥,其爆炸過(guò)程中爆轟產(chǎn)物燃燒釋放的能量占總能量比重的近40%[2],且作用時(shí)間尺度為毫秒量級(jí),一般在結(jié)構(gòu)的動(dòng)態(tài)響應(yīng)時(shí)間范圍內(nèi)。與敞開(kāi)環(huán)境不同,封閉空間內(nèi)部爆炸的特征主要表現(xiàn)為沖擊波效應(yīng)顯著增強(qiáng)、熱效應(yīng)明顯、爆轟產(chǎn)物膨脹做功的有效作用時(shí)間更長(zhǎng)等[3],使得爆炸載荷的燃燒增強(qiáng)效應(yīng)對(duì)結(jié)構(gòu)響應(yīng)和破壞的影響更加顯著。

    針對(duì)封閉空間內(nèi)爆炸產(chǎn)物燃燒效應(yīng)問(wèn)題,國(guó)內(nèi)外開(kāi)展了試驗(yàn)及數(shù)值計(jì)算方面的研究工作。Ornellas[4]發(fā)現(xiàn)了爆炸產(chǎn)物的燃燒反應(yīng)對(duì)約束空間TNT爆炸能量和壓力的影響。Ferguson等[5]和Kuhl等[6-7]建立了描述TNT在封閉空間內(nèi)爆炸后發(fā)生燃燒反應(yīng)的氣動(dòng)模型,同時(shí)開(kāi)展了沖擊分散炸藥(SDF)在不同尺寸密閉容器內(nèi)部爆炸的試驗(yàn)[8],探索了艙室容積和幾何形狀及艙內(nèi)氣體環(huán)境對(duì)燃燒完全性的影響,并在試驗(yàn)研究的基礎(chǔ)上,通過(guò)引入詳細(xì)的化學(xué)反應(yīng)動(dòng)力學(xué)過(guò)程和湍流混合模型,建立了可較好地模擬封閉空間內(nèi)炸藥爆炸過(guò)程的數(shù)值計(jì)算方法[9-10]。金朋剛等[11-12]的研究發(fā)現(xiàn)氧氣濃度對(duì)封閉空間內(nèi)TNT爆炸產(chǎn)物燃燒的影響顯著。鐘巍等[13-14]通過(guò)開(kāi)展試驗(yàn)研究、理論分析和數(shù)值計(jì)算的方法,提出了化學(xué)反應(yīng)動(dòng)力學(xué)影響的約束空間TNT爆炸的準(zhǔn)靜態(tài)壓力計(jì)算方法。

    在準(zhǔn)靜態(tài)壓力的影響因素方面,王等旺等[15]通過(guò)試驗(yàn)和理論分析得到了封閉空間TNT爆炸準(zhǔn)靜態(tài)壓力和裝藥量與體積比的擬合關(guān)系。張玉磊等[16-17]基于爆炸相似律與理想氣體狀態(tài)方程分析建立了沖擊波超壓及準(zhǔn)靜態(tài)壓力計(jì)算模型,并開(kāi)展了溫壓炸藥和TNT的內(nèi)爆炸試驗(yàn),發(fā)現(xiàn)溫壓炸藥內(nèi)爆沖擊波壓力峰值和準(zhǔn)靜態(tài)壓力峰值比TNT高,藥量—容積比是準(zhǔn)靜態(tài)壓力的主要影響因素。Feldgun等[18]基于熱力學(xué)理論,開(kāi)發(fā)了計(jì)算限制空間爆炸壓力的熱力學(xué)模型,研究了準(zhǔn)靜態(tài)壓力對(duì)艙內(nèi)氣體比熱比和內(nèi)能的敏感性,結(jié)果表明艙內(nèi)氣體的比熱比對(duì)準(zhǔn)靜態(tài)壓力的影響大于內(nèi)能。在考慮了艙內(nèi)氣體比熱比和爆炸釋放總能量的變化后,準(zhǔn)靜態(tài)壓力取決于炸藥質(zhì)量與容器體積比。

    目前,國(guó)內(nèi)外研究人員在密閉空間TNT內(nèi)爆的沖擊波壓力特性、后燃燒效應(yīng)及準(zhǔn)靜態(tài)壓力的影響因素方面已經(jīng)開(kāi)展了系列的研究工作。然而,在密閉空間內(nèi)爆炸載荷的燃燒增強(qiáng)效應(yīng)對(duì)結(jié)構(gòu)響應(yīng)影響方面的研究開(kāi)展得較少。本文開(kāi)展空氣和氦氣環(huán)境中艙室密閉空間內(nèi)的TNT爆炸載荷以及結(jié)構(gòu)響應(yīng)的對(duì)比試驗(yàn),分析艙內(nèi)爆炸載荷的燃燒增強(qiáng)效應(yīng)對(duì)結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的影響。

    1 TNT封閉空間內(nèi)爆炸試驗(yàn)

    1.1 試驗(yàn)裝置

    為開(kāi)展對(duì)比試驗(yàn),設(shè)計(jì)了如圖1所示的爆炸筒試驗(yàn)裝置。爆炸筒由高強(qiáng)度鋼焊接成型,內(nèi)部空間尺寸為900 mm×400 mm×400 mm,左右兩側(cè)端部開(kāi)螺栓孔,通過(guò)壓板及螺栓固定端部的試件。筒身焊接縱橫加強(qiáng)筋保證結(jié)構(gòu)強(qiáng)度,下部通過(guò)固定工裝與地基連接,確保裝置的穩(wěn)固。在筒身不同部位設(shè)置溫度傳感器及壓力傳感器安裝基座,同時(shí)在筒身對(duì)角處設(shè)置閥門(mén),實(shí)現(xiàn)艙內(nèi)氣體的替換。傳感器布置位置如圖2所示,其中P1~P8和SP1~SP3為壓力傳感器,T1、T2、T4、T5、T7、T8為溫度傳感器。

    圖1 爆炸試驗(yàn)裝置示意圖

    圖2 溫度及壓力傳感器布置圖

    本試驗(yàn)中采用的TNT由粉末壓制成型,密度為1.55 g/cm3.壓力傳感器P1、P2、P3、P4、P5、P6、P7、P8采用美國(guó)PCB公司生產(chǎn)的102B型壓力傳感器,壓力傳感器SP1和SP3采用美國(guó)PCB公司生產(chǎn)的113B22型壓力傳感器,兩種傳感器的量程均為34.5 MPa;溫度傳感器為美國(guó)NANMAC公司生產(chǎn)的C2-7-K-L型熱電偶,采用法國(guó)Texense公司生產(chǎn)的THNF-A型熱電偶放大器,量程為0~1 250 ℃.采集系統(tǒng)為德國(guó)HBM公司生產(chǎn)的Genesis 7T型高速數(shù)據(jù)采集系統(tǒng),采樣頻率為1×106Hz,在每個(gè)測(cè)試系列的開(kāi)始和結(jié)束時(shí)所有的測(cè)點(diǎn)通道都進(jìn)行校準(zhǔn)。

    試驗(yàn)中采用德國(guó)GOM公司生產(chǎn)的數(shù)字圖像相關(guān)(DIC)方法測(cè)量系統(tǒng)對(duì)內(nèi)爆載荷下試件的變形歷程進(jìn)行測(cè)量,試驗(yàn)布置如圖3所示。根據(jù)測(cè)試范圍(ROI)的大小(400 mm×400 mm)決定分辨率設(shè)置為960×960,幀率為1×104幀/s,每個(gè)像素長(zhǎng)度約0.4 mm.DIC方法的匹配精度為5.4%像素長(zhǎng)度,本試驗(yàn)中其變形測(cè)量精度為0.02 mm.

    圖3 3D-DIC布置形式

    1.2 試驗(yàn)工況

    空氣工況下,關(guān)閉氦氣進(jìn)氣閥和出氣閥。將TNT懸掛于試驗(yàn)裝置(見(jiàn)圖1)內(nèi)部的中心位置,在試件鋼板表面用白色油漆噴涂均勻,并用黑色油漆噴涂散斑達(dá)到DIC方法的測(cè)量要求。TNT懸掛安裝后采用螺栓和壓板將試件固定在爆炸筒端部,形成封閉環(huán)境,隨后引爆TNT測(cè)量數(shù)據(jù)。

    氦氣工況下,待試件安裝之后,打開(kāi)氦氣進(jìn)氣閥與空氣出氣閥,并在出氣口用氧氣濃度檢測(cè)儀檢測(cè)氧氣濃度,至氧氣濃度下降到2.2%時(shí)停止充氣,關(guān)閉出氣閥與進(jìn)氣閥,引爆TNT并測(cè)量數(shù)據(jù)。

    本次試驗(yàn)設(shè)置了從7.50~30.00 g 5種TNT藥量、空氣與氦氣艙內(nèi)環(huán)境下共10種試驗(yàn)工況。試驗(yàn)記錄工況編號(hào)具體含義見(jiàn)圖4,具體工況設(shè)置如表1所示。

    圖4 試驗(yàn)記錄工況編號(hào)含義

    2 試驗(yàn)結(jié)果與分析

    2.1 各工況下艙內(nèi)爆炸載荷對(duì)比分析

    2.1.1 試驗(yàn)結(jié)果與分析

    由于爆炸筒的端部試件在爆炸過(guò)程中會(huì)產(chǎn)生變形,端部螺栓孔附近會(huì)產(chǎn)生縫隙,造成少量的泄壓,因此壓力歷程曲線在后續(xù)準(zhǔn)靜態(tài)壓力階段出現(xiàn)了線性緩降的過(guò)程。為更好地反映艙內(nèi)準(zhǔn)靜態(tài)壓力情況,本文對(duì)試驗(yàn)測(cè)得的壓力歷程曲線以8 ms為區(qū)間進(jìn)行平均處理,得到圖5中所示的準(zhǔn)靜態(tài)壓力曲線。

    表1 工況設(shè)置

    圖5 15.00 g TNT藥量工況下P8測(cè)點(diǎn)的壓力歷程曲線

    由于數(shù)據(jù)較多,選取15.00 g TNT工況下空氣和氦氣環(huán)境中各測(cè)點(diǎn)的壓力峰值和準(zhǔn)靜態(tài)壓力數(shù)據(jù)列于表2,同時(shí)選取該工況P8測(cè)點(diǎn)的沖擊波壓力和準(zhǔn)靜態(tài)壓力歷程曲線如圖5所示。通過(guò)分析數(shù)據(jù)發(fā)現(xiàn),不同位置測(cè)點(diǎn)的沖擊波壓力峰值有一定差異,而準(zhǔn)靜態(tài)壓力則相差不大,表明準(zhǔn)靜態(tài)壓力在封閉空間內(nèi)具有均勻壓力載荷的特性。7.50 g、11.25 g、15.00 g、22.50 g、30.00 g TNT爆炸工況下空氣和氦氣環(huán)境中各測(cè)點(diǎn)的準(zhǔn)靜態(tài)壓力平均值數(shù)據(jù)如表3所示。通過(guò)對(duì)比發(fā)現(xiàn),由于氦氣對(duì)爆炸產(chǎn)物燃燒效應(yīng)的抑制,各藥量工況下其準(zhǔn)靜態(tài)壓力較空氣環(huán)境中降低幅度在38.81%~46.85%之間。

    從以上5組試驗(yàn)數(shù)據(jù)的對(duì)比中可以發(fā)現(xiàn),當(dāng)TNT爆炸發(fā)生在封閉空間內(nèi)時(shí),爆炸產(chǎn)物的燃燒效應(yīng)對(duì)反射沖擊波和準(zhǔn)靜態(tài)壓力具有明顯的增強(qiáng)作用??諝夤r下,壓力峰值出現(xiàn)在初始沖擊波到達(dá)之后的極短時(shí)間內(nèi),由于氧氣充足,且空間限制使得爆炸產(chǎn)物與空氣充分混合,燃燒釋放的能量使內(nèi)部壓力增加,進(jìn)而與艙壁反射的沖擊波疊加,使反射沖擊波的峰值壓力和后續(xù)的準(zhǔn)靜態(tài)壓力明顯提高。而在氦氣環(huán)境中,壓力峰值為初始沖擊波峰值,由于艙內(nèi)氧氣濃度僅為2.2%,其惰性氣體的性質(zhì)抑制了TNT爆炸產(chǎn)物的進(jìn)一步燃燒,導(dǎo)致TNT爆炸產(chǎn)物的能量不能完全釋放。試驗(yàn)現(xiàn)象表明,在氦氣工況下,TNT爆炸后筒壁內(nèi)有大量未燃燒的碳粉,如圖6所示,說(shuō)明TNT爆炸產(chǎn)物基本未發(fā)生燃燒效應(yīng)。而在空氣工況下,TNT爆炸后,筒壁內(nèi)并未發(fā)現(xiàn)殘留碳粉,如圖7所示。

    表2 15.00 g TNT藥量工況下各測(cè)點(diǎn)的壓力峰值及準(zhǔn)靜態(tài)壓力

    表3 不同工況下準(zhǔn)靜態(tài)壓力平均值對(duì)比

    圖6 氦氣工況下TNT爆炸后艙內(nèi)情況

    圖7 空氣工況下TNT爆炸后艙內(nèi)情況

    此外,通過(guò)對(duì)比各測(cè)點(diǎn)記錄的空氣和氦氣環(huán)境中沖擊波壓力時(shí)程曲線可以發(fā)現(xiàn),TNT爆炸產(chǎn)物燃燒造成的反射沖擊波增強(qiáng)發(fā)生在初始沖擊波到達(dá)之后很短的時(shí)間內(nèi)(1 ms),且持續(xù)時(shí)間在3 ms左右。這一時(shí)間段內(nèi)的爆炸產(chǎn)物燃燒釋放能量直接導(dǎo)致了艙內(nèi)準(zhǔn)靜態(tài)壓力的大幅提高。

    2.1.2 準(zhǔn)靜態(tài)壓力理論計(jì)算

    Feldgun等[18]給出了密閉空間中爆炸載荷準(zhǔn)靜態(tài)壓力的計(jì)算公式為

    (1)

    式中:p0=101.3 kPa為大氣壓力;γ、γ0分別為爆炸后混合氣體和爆炸前艙內(nèi)氣體的絕熱指數(shù);ρE為炸藥密度;W為炸藥質(zhì)量;V為封閉空間的體積;Q為單位質(zhì)量炸藥釋放的能量。當(dāng)W/V<0.387 kg/m3時(shí),由于爆炸后混合氣體成分難以確定,為簡(jiǎn)化計(jì)算,認(rèn)為空氣工況下γ=γ0=1.4,氦氣工況下γ=γ0=1.621,氦氣工況下γ0由表4得到。

    氦氣工況下,爆炸前封閉艙室內(nèi)氣體主要為大部分氦氣與少量空氣的混合氣體,氧氣濃度檢測(cè)儀在氦氣工況下測(cè)得艙內(nèi)氧氣濃度為2.2%,因此艙內(nèi)空氣體積占比為10.48%,氦氣的體積占比為89.52%,爆炸前艙內(nèi)混合氣體的參數(shù)如表4所示。

    表4 氦氣工況下爆炸前艙內(nèi)混合氣體參數(shù)

    將試驗(yàn)測(cè)量得到的準(zhǔn)靜態(tài)壓力和其他具體數(shù)值代入(1)式,并計(jì)入雷管的爆轟能量(1 g TNT當(dāng)量)可得到各工況下導(dǎo)致準(zhǔn)靜態(tài)壓力升高的等效能量Qe,如表5所示。

    表5 不同工況下導(dǎo)致準(zhǔn)靜態(tài)壓力升高的等效能量

    5種不同藥量TNT在空氣和氦氣環(huán)境封閉空間內(nèi)釋放的能量如表5所示。從數(shù)據(jù)對(duì)比來(lái)看,由于氦氣對(duì)爆炸產(chǎn)物燃燒效應(yīng)的抑制,較空氣工況下等效能量下降的比例在63.4%~68.1%之間。

    在氦氣工況下若不考慮燃燒效應(yīng),TNT的爆轟能量為4 200 kJ/kg;在空氣工況下,若TNT完全燃燒,則單位質(zhì)量TNT釋放的能量為14 827.5 kJ/kg.分析表5數(shù)據(jù)可知,7.50 g TNT氦氣工況的等效能量略大于初始爆轟能量,原因可能是小部分TNT爆炸產(chǎn)物在低濃度氧氣中發(fā)生反應(yīng),使等效能量偏大??諝夂秃夤r下,隨著TNT質(zhì)量的升高,用于提高艙內(nèi)壓力的TNT等效能量逐漸降低,且均低于TNT的理論釋放能量;這說(shuō)明一部分爆轟能量用于結(jié)構(gòu)變形及傳遞至外部耗散,且隨著TNT質(zhì)量的增大,用于結(jié)構(gòu)變形及耗散的能量比例增大。對(duì)比氦氣及空氣工況,不同TNT質(zhì)量下氦氣相較空氣工況等效能量下降的比例變化不大,約為65.5%,與理論情況的71.67%接近。說(shuō)明相同TNT質(zhì)量下,空氣與氦氣工況用于升高氣體壓力的能量比例基本不變。因此在理論計(jì)算各工況準(zhǔn)靜態(tài)壓力時(shí),使用表5中的等效能量進(jìn)行計(jì)算。

    艙室密閉空間內(nèi)中爆炸的能量釋放大致分為兩部分:初始爆轟階段和燃燒增強(qiáng)階段。在初始爆轟階段,能量釋放時(shí)間為微秒量級(jí),產(chǎn)生大量的高溫氣體產(chǎn)物,形成初始沖擊波;燃燒增強(qiáng)階段在初始爆轟后的1 ms時(shí)間內(nèi),高溫爆轟產(chǎn)物擴(kuò)散,爆轟產(chǎn)物之間、產(chǎn)物與空氣中的氧氣相互反應(yīng),放熱升溫導(dǎo)致密閉空間壓力繼續(xù)上升,兩階段釋能最終形成密閉空間內(nèi)的高溫高壓環(huán)境。

    在氦氣工況下,不考慮燃燒增強(qiáng)效應(yīng),艙內(nèi)準(zhǔn)靜態(tài)壓力由爆轟產(chǎn)物在密閉艙室中產(chǎn)生的氣體壓力升高pg和初始爆轟能量加熱氣體導(dǎo)致的壓力升高pt所組成。根據(jù)理想氣體狀態(tài)方程:

    pg=p0V0(m/V),

    (2)

    式中:V0為T(mén)NT爆容;m為T(mén)NT質(zhì)量。

    假定爆炸反應(yīng)釋放的能量全部用來(lái)加熱氣體,那么氣體的溫度升高為

    (3)

    式中:ΔT為密閉空間內(nèi)氣體溫度升高;QV為T(mén)NT爆熱;mg為爆炸后密閉空間內(nèi)混合氣體質(zhì)量;cV為爆炸后密閉空間內(nèi)混合氣體定容比熱。

    由理想氣體狀態(tài)方程可得

    (4)

    故理想狀態(tài)下準(zhǔn)靜態(tài)壓力pqs的計(jì)算公式為

    (5)

    式中:V0=800 L/kg;QV=Qe;R=8.314(kPa·L)/(mol·K).

    TNT爆炸的方程式為C7H5O6N3=2.5H2O+3.5CO+3.5C+1.5N2,即1 mol TNT爆炸產(chǎn)生2.5 mol水蒸氣、3.5 mol一氧化碳和1.5 mol氮?dú)?,?5 ℃下的定容比熱容分別為1.400 kJ/(kg·K)、0.743 kJ/(kg·K)、0.741 kJ/(kg·K)。將上述參數(shù)代入(5)式,并代入表5中的Qe,考慮雷管的爆轟能量(1 g TNT當(dāng)量),計(jì)算各藥量下的準(zhǔn)靜態(tài)壓力,如表6所示,繪制準(zhǔn)靜態(tài)壓力試驗(yàn)值與計(jì)算值隨藥量體積比的變化關(guān)系,如圖8所示。

    表6 氦氣工況下各藥量的準(zhǔn)靜態(tài)壓力計(jì)算值

    圖8 氦氣工況下準(zhǔn)靜態(tài)壓力計(jì)算值與試驗(yàn)值隨m/V的變化關(guān)系

    在空氣工況下,由于藥量體積比小于0.371 3 kg/m3,TNT爆轟產(chǎn)物與空氣中氧氣發(fā)生如下反應(yīng):

    假定試驗(yàn)中TNT均完全反應(yīng),則1 mol TNT爆炸燃燒消耗5.25 mol氧氣,產(chǎn)生2.5 mol水蒸氣、7 mol二氧化碳、1.5 mol氮?dú)?。氧氣與二氧化碳在25 ℃下的定容比熱容分別為0.657 kJ/(kg·K)和0.638 kJ/(kg·K)。由于反應(yīng)中有氧氣的消耗,將空氣成分假設(shè)為由79%氮?dú)馀c21%氧氣組成,氮?dú)馀c氧氣在25 ℃下的密度分別為1.25 kg/m3和1.43 kg/m3.將上述參數(shù)代入(5)式,并代入表5中的Qe,考慮雷管的爆轟能量(1 g TNT當(dāng)量),不考慮雷管的燃燒能量,計(jì)算各藥量下的準(zhǔn)靜態(tài)壓力,如表7所示,繪制準(zhǔn)靜態(tài)壓力試驗(yàn)值與計(jì)算值隨藥量體積比的變化關(guān)系,如圖9所示。

    通過(guò)對(duì)比表6、表7和圖8、圖9中的數(shù)據(jù)發(fā)現(xiàn),在氦氣和空氣工況中理論計(jì)算值均低于試驗(yàn)值,且隨著藥量體積比的增大,試驗(yàn)值與理論計(jì)算值的差距增大。出現(xiàn)這一現(xiàn)象的主要原因是:使用(1)式計(jì)算Qe時(shí),假設(shè)爆炸后混和氣體的絕熱系數(shù)與爆炸前艙內(nèi)氣體的絕熱系數(shù)相同,但實(shí)際情況中,氦氣工況下TNT爆炸產(chǎn)物的絕熱系數(shù)均低于氦氣的絕熱系數(shù),空氣工況下TNT爆炸產(chǎn)物二氧化碳的絕熱系數(shù)低于空氣的絕熱系數(shù)。因此,爆炸后艙內(nèi)混合氣體的絕熱指數(shù)比爆炸前偏小,即γ<γ0,從而計(jì)算得到的Qe偏小,進(jìn)而計(jì)算得到的準(zhǔn)靜態(tài)壓力偏小;同時(shí),隨著TNT質(zhì)量的增加,爆炸產(chǎn)物增加,Qe的偏差會(huì)增大,因此計(jì)算得到的準(zhǔn)靜態(tài)壓力的偏差也會(huì)增大。

    表7 空氣工況下各藥量的準(zhǔn)靜態(tài)壓力計(jì)算值

    圖9 空氣工況下準(zhǔn)靜態(tài)壓力計(jì)算值與試驗(yàn)值隨m/V的變化關(guān)系

    2.2 溫度測(cè)試結(jié)果與分析

    2.2.1 試驗(yàn)結(jié)果分析

    艙內(nèi)爆炸及燃燒過(guò)程產(chǎn)生的溫度場(chǎng)不規(guī)則,并且溫度傳感器伸入筒內(nèi)105 mm進(jìn)行測(cè)量,不同溫度測(cè)點(diǎn)距離爆炸中心的距離不同(T1、T4、T7測(cè)點(diǎn)離爆炸中心的距離為249 mm,T2和T5測(cè)點(diǎn)離爆炸中心的距離為415 mm),因此各溫度測(cè)點(diǎn)的溫度數(shù)據(jù)在同一工況下有一定差異。通過(guò)對(duì)距爆炸中心距離相同測(cè)點(diǎn)的溫度峰值平均值來(lái)表征各工況的溫度數(shù)據(jù)。

    空氣環(huán)境中不同TNT藥量工況下的T7測(cè)點(diǎn)溫度變化如圖10所示;氦氣環(huán)境中不同TNT藥量工況下T1測(cè)點(diǎn)的溫度變化如圖11所示。通過(guò)對(duì)比發(fā)現(xiàn),空氣和氦氣工況下,艙內(nèi)溫度隨TNT炸藥質(zhì)量的增大而升高,TNT質(zhì)量相同時(shí)氦氣工況的艙內(nèi)溫度較空氣工況的溫度大幅度降低。對(duì)距離爆炸中心249 mm的T1、T4和T7測(cè)點(diǎn)數(shù)據(jù)取平均值、對(duì)距離爆炸中心415 mm的T2、T5測(cè)點(diǎn)數(shù)據(jù)取平均值列于表8.通過(guò)對(duì)比距爆炸中心不同距離的溫度測(cè)點(diǎn)數(shù)據(jù)發(fā)現(xiàn):對(duì)于較小藥量的氦氣及空氣環(huán)境下爆炸,距離爆炸中心越近,溫度越高;對(duì)于TNT質(zhì)量為30.00 g、艙內(nèi)環(huán)境為空氣的工況,爆炸后產(chǎn)生明顯燃燒效應(yīng),導(dǎo)致爆炸中心向外輻射的大范圍區(qū)域均為高溫燃燒區(qū)域,因此距離爆炸中心不同距離的測(cè)點(diǎn)溫度差異不明顯。進(jìn)一步處理得到氦氣工況相對(duì)空氣工況溫度峰值下降的比例,如表9所示,由于氦氣對(duì)TNT爆炸產(chǎn)物燃燒過(guò)程的抑制,爆炸后艙內(nèi)溫度大幅度降低,5種藥量工況下溫度降幅在57.5%~76.4%之間,進(jìn)一步說(shuō)明爆炸產(chǎn)物的燃燒是艙內(nèi)氣體高溫的主要因素。

    圖10 不同TNT藥量在空氣工況下T7測(cè)點(diǎn)的溫度歷程曲線

    圖11 不同TNT藥量在氦氣工況下T1測(cè)點(diǎn)的溫度歷程曲線

    表8 不同工況下不同測(cè)點(diǎn)位置的溫度峰值平均值

    表9 氦氣工況相對(duì)空氣工況的溫度峰值下降比例

    2.2.2 溫度升高理論計(jì)算

    將表5中Qe和表6、表7中的mg、cV數(shù)據(jù)代入(3)式,計(jì)算得到各工況下艙內(nèi)溫度的升高,如表10所示。

    表10 各工況下溫度升高計(jì)算值

    很明顯,理論計(jì)算的溫度升高明顯大于試驗(yàn)測(cè)量值,主要是因?yàn)闊犭娕嫉捻憫?yīng)時(shí)間為20 ms,測(cè)得的溫度具有一定的滯后性,無(wú)法反映艙內(nèi)溫度的瞬態(tài)變化歷程。

    2.3 試件變形結(jié)果及分析

    2.3.1 試件變形試驗(yàn)結(jié)果

    本文試驗(yàn)中測(cè)量了爆炸筒端部安裝的金屬板試件的變形響應(yīng),可用以分析艙內(nèi)爆炸燃燒增強(qiáng)效應(yīng)對(duì)結(jié)構(gòu)響應(yīng)的影響程度。通過(guò)分析DIC的測(cè)量結(jié)果,提取試驗(yàn)試件中心點(diǎn)的變形響應(yīng)歷程曲線。其中MB1.5-5-11.25(Air)和MB3-1-22.5(Air)工況的DIC數(shù)據(jù)在試驗(yàn)過(guò)程中未成功采集,在試驗(yàn)結(jié)束后采用激光掃描儀測(cè)量該工況下試件的最終變形,各工況下試件中心點(diǎn)的變形歷程曲線如圖12~圖16所示。在封閉空間內(nèi)爆載荷作用下,各試件的變形歷程相似,開(kāi)始響應(yīng)之后1.3 ms內(nèi)即達(dá)到初始響應(yīng)峰值,在初始響應(yīng)峰值之后發(fā)生回彈,并在一定范圍內(nèi)振蕩。隨著藥量的增加,回彈量值和振蕩范圍逐漸減小。相對(duì)比來(lái)看,氦氣工況下達(dá)到初始響應(yīng)峰值的時(shí)間與空氣工況下的差異很小,然而初始響應(yīng)峰值和振蕩平衡位置大大降低,且初始峰值之后的回彈量增加。從圖16可以看出,在30 g TNT藥量?jī)?nèi)爆工況下,空氣工況下較大的準(zhǔn)靜態(tài)壓力使得試件初始響應(yīng)峰值和振蕩平衡位置處的差別較小,且出現(xiàn)振蕩平衡位置略微增加的趨勢(shì)。

    圖12 7.5 g TNT藥量在空氣與氦氣工況下試件中心的變形歷程曲線

    圖13 11.25 g TNT藥量在氦氣工況下試件中心的變形歷程曲線

    圖14 15 g TNT藥量在空氣與氦氣工況下試件中心的變形歷程曲線

    圖15 22.5 g TNT藥量在氦氣工況下試件中心的變形歷程曲線

    圖16 30 g TNT藥量在空氣與氦氣工況下試件中心的變形歷程曲線

    將變形歷程曲線的振蕩階段取平均值得到試件中心的最終變形,結(jié)合三維掃描得到的最終變形數(shù)據(jù),各工況中試驗(yàn)試件中心最終變形如表11所示。對(duì)比表11中數(shù)據(jù)發(fā)現(xiàn),不同TNT藥量下氦氣工況試件最終變形均遠(yuǎn)小于空氣工況。5種工況下試件最終變形的降幅分別達(dá)到了48.9%、38.1%、36.9%、19.1%和25.0%.

    2.3.2 試件變形理論分析

    文獻(xiàn)[19-20]給出了在均布爆炸載荷作用下金屬薄板中點(diǎn)變形- 厚度比的計(jì)算公式,即

    (6)

    式中:δ為板的中點(diǎn)變形;h和L分別為板的厚度和邊長(zhǎng);I為爆炸載荷沖量;ρ和σy分別為材料的密度和屈服應(yīng)力,其中ρ=7.83 g/cm3,1.5 mm板的屈服應(yīng)力為314.6 MPa,3 mm板的屈服應(yīng)力307 MPa.

    表11 各工況下試件中心的最終變形

    根據(jù)(6)式,可得到本試驗(yàn)中內(nèi)爆載荷的等效沖量:

    Ieq=4.17×δhL(ρσy)1/2.

    (7)

    金屬薄板在爆炸載荷作用下的響應(yīng)具有飽和特性,即當(dāng)金屬板結(jié)構(gòu)受到持續(xù)時(shí)間較長(zhǎng)的爆炸沖擊波載荷作用時(shí)發(fā)生彈塑性動(dòng)態(tài)響應(yīng),繼而產(chǎn)生較大的塑性變形;當(dāng)板的變形撓度遠(yuǎn)大于其板厚時(shí),由于四周邊界的約束作用,在板的中面內(nèi)會(huì)產(chǎn)生較大的膜力效應(yīng),并使得板結(jié)構(gòu)的受載能力進(jìn)一步增強(qiáng),使之變形在達(dá)到一定值之后就不會(huì)再繼續(xù)增加;此時(shí)金屬板的動(dòng)態(tài)響應(yīng)達(dá)到“飽和狀態(tài)”,載荷的后續(xù)階段并不會(huì)對(duì)板結(jié)構(gòu)的變形響應(yīng)產(chǎn)生影響[21-22],并且在內(nèi)爆載荷作用下,板結(jié)構(gòu)達(dá)到最大變形響應(yīng)的時(shí)間僅與板長(zhǎng)度、材料密度、屈服強(qiáng)度等因素有關(guān)[23]。根據(jù)試驗(yàn)結(jié)果可以發(fā)現(xiàn),在不同工況下板試件達(dá)到最大變形的時(shí)間(飽和響應(yīng)時(shí)間)基本保持不變。將各工況中試驗(yàn)試件的詳細(xì)參數(shù)代入(7)式,可得到作用于試驗(yàn)試件的等效沖量,進(jìn)而可得到在板的飽和響應(yīng)時(shí)間范圍內(nèi)作用于板的等效平均壓力pa,如表12所示,表明在試件的飽和響應(yīng)時(shí)間內(nèi)爆炸產(chǎn)物的燃燒效應(yīng)對(duì)試件的動(dòng)態(tài)響應(yīng)起到了明顯的增強(qiáng)作用。

    表12 不同工況下內(nèi)爆載荷的等效沖量和壓力

    從艙內(nèi)爆炸載荷的試驗(yàn)數(shù)據(jù)分析可以發(fā)現(xiàn),對(duì)爆炸產(chǎn)物燃燒效應(yīng)的抑制可有效地降低反射沖擊波和準(zhǔn)靜態(tài)壓力,然而,由于結(jié)構(gòu)響應(yīng)的“飽和特性”,能否有效減小金屬板在內(nèi)爆載荷作用下的響應(yīng)則取決于能否“及時(shí)”抑制艙內(nèi)爆炸載荷。對(duì)比表11與表12中7.50 g、15.00 g和30.00 g TNT藥量工況可知,氦氣工況相對(duì)于空氣工況,板中心點(diǎn)位移的下降比例與等效壓力下降比例接近,因此造成板最終變形的主要原因是飽和響應(yīng)時(shí)間內(nèi)的等效壓力。試驗(yàn)中板的飽和響應(yīng)時(shí)間約為開(kāi)始變形后的1.1 ms,板開(kāi)始變形的時(shí)刻即為載荷初始沖擊波達(dá)到的時(shí)刻。分析各工況沖擊波載荷歷程曲線,以11.25 g TNT藥量工況為例,繪制爆炸后5 ms內(nèi)空氣及氦氣工況的載荷歷程曲線如圖17所示,在空氣和氦氣工況下,初始沖擊波的峰值大小相差不大,而隨著爆炸產(chǎn)物的燃燒,空氣工況下在初始沖擊波到達(dá)后1 ms內(nèi)的反射沖擊波顯著增強(qiáng)。這種增強(qiáng)效應(yīng)發(fā)生在試件的飽和響應(yīng)時(shí)間內(nèi),載荷的增強(qiáng)與試件的響應(yīng)過(guò)程疊加,造成了空氣工況下試件的初始響應(yīng)峰值和后續(xù)振蕩平衡位置較氦氣工況下顯著提高。因此,在結(jié)構(gòu)內(nèi)爆響應(yīng)的評(píng)估中不能忽略爆炸載荷燃燒增強(qiáng)效應(yīng)的影響。

    圖17 11.25 g TNT藥量在爆炸工況下P5測(cè)點(diǎn)的壓力歷程曲線

    3 結(jié)論

    本文開(kāi)展了TNT炸藥爆炸產(chǎn)物的燃燒釋放能量對(duì)封閉空間中爆炸載荷影響的試驗(yàn)研究,測(cè)量得到了艙內(nèi)環(huán)境氣體環(huán)境為氦氣和空氣工況下5種不同TNT藥量爆炸載荷的壓力時(shí)程和準(zhǔn)靜態(tài)壓力、艙內(nèi)氣體溫度變化和試驗(yàn)試件的動(dòng)態(tài)響應(yīng)歷程和最終變形。通過(guò)試驗(yàn)數(shù)據(jù)和試驗(yàn)現(xiàn)象分析,得出了以下結(jié)論:

    1)TNT炸藥爆炸產(chǎn)物的燃燒對(duì)封閉空間內(nèi)反射沖擊波和準(zhǔn)靜態(tài)壓力的增強(qiáng)效應(yīng)明顯,且該增強(qiáng)過(guò)程在初始沖擊波到達(dá)后的極短時(shí)間內(nèi)發(fā)生。通過(guò)采用在封閉內(nèi)注入氦氣的方法抑制爆炸產(chǎn)物的進(jìn)一步燃燒釋放能量之后,TNT藥量為7.50 g、11.25 g、15.00 g、22.50 g、30.00 g時(shí),最大準(zhǔn)靜態(tài)壓力相較空氣工況分別下降了38.81%、46.85%、43.87%、45.05%、42.90%,溫度峰值相較空氣工況分別下降了70.73%、76.35%、66.45%、57.53%、62.92%,降低了艙內(nèi)爆炸載荷和艙內(nèi)氣體環(huán)境溫度。

    2)較空氣工況下試驗(yàn)試件的響應(yīng),氦氣工況下受到5種不同TNT藥量的爆炸載荷時(shí),試件的初始響應(yīng)峰值和最終變形均大幅減小,最終變形降幅分別達(dá)到了48.9%、38.1%、36.9%、19.1%和25.0%,表明爆炸產(chǎn)物燃燒的能量釋放率較高,能量釋放時(shí)間與金屬板試件的響應(yīng)時(shí)間在同一量級(jí)。在結(jié)構(gòu)內(nèi)爆響應(yīng)分析和安全評(píng)估中需考慮爆炸載荷燃燒增強(qiáng)效應(yīng)的影響。

    猜你喜歡
    氦氣藥量靜態(tài)
    神奇的氦氣
    碳化對(duì)飛灰穩(wěn)定化加藥量測(cè)評(píng)的影響
    靜態(tài)隨機(jī)存儲(chǔ)器在軌自檢算法
    兩倍藥量
    跟氣球上天
    廉政瞭望(2020年17期)2020-11-17 07:37:32
    坦桑尼亞發(fā)現(xiàn)巨型氦氣礦
    低溫與特氣(2018年1期)2018-04-16 13:19:36
    飛走的氦氣球能飛多高?
    小型水池抗爆炸沖擊極限藥量分析
    藥量適中
    機(jī)床靜態(tài)及動(dòng)態(tài)分析
    夜夜躁狠狠躁天天躁| 日韩熟女老妇一区二区性免费视频| 久久中文看片网| 欧美午夜高清在线| 黑丝袜美女国产一区| 大香蕉久久成人网| 丰满人妻熟妇乱又伦精品不卡| 91在线观看av| 亚洲av成人不卡在线观看播放网| 女性生殖器流出的白浆| 国产精品美女特级片免费视频播放器 | 建设人人有责人人尽责人人享有的| 国产真人三级小视频在线观看| 国产一区二区三区在线臀色熟女 | 久久青草综合色| 黑人巨大精品欧美一区二区mp4| 亚洲美女黄片视频| 亚洲av第一区精品v没综合| 国产欧美日韩一区二区精品| 夜夜爽天天搞| 999久久久精品免费观看国产| 国产乱人伦免费视频| 精品福利观看| av线在线观看网站| 亚洲成人国产一区在线观看| 亚洲av电影在线进入| 一个人免费在线观看的高清视频| 极品教师在线免费播放| 捣出白浆h1v1| 日韩欧美一区视频在线观看| av线在线观看网站| av在线播放免费不卡| 久久久久精品人妻al黑| 不卡一级毛片| 久热爱精品视频在线9| 91精品国产国语对白视频| 桃红色精品国产亚洲av| 18禁裸乳无遮挡免费网站照片 | 日本撒尿小便嘘嘘汇集6| 久久国产亚洲av麻豆专区| 成人av一区二区三区在线看| 久久久国产精品麻豆| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 精品亚洲成a人片在线观看| 久久久国产精品麻豆| 日本一区二区免费在线视频| 王馨瑶露胸无遮挡在线观看| 日本欧美视频一区| 高清黄色对白视频在线免费看| 少妇的丰满在线观看| 色播在线永久视频| 在线视频色国产色| 亚洲精品中文字幕一二三四区| 中文字幕色久视频| 欧美日韩精品网址| 80岁老熟妇乱子伦牲交| 久久精品亚洲av国产电影网| 亚洲精品av麻豆狂野| 国产精品综合久久久久久久免费 | 超碰97精品在线观看| 亚洲精品一二三| 国产精品久久视频播放| 久久影院123| 日韩中文字幕欧美一区二区| 欧美激情极品国产一区二区三区| 日本撒尿小便嘘嘘汇集6| 韩国av一区二区三区四区| 日韩欧美国产一区二区入口| 精品国产乱码久久久久久男人| 亚洲成a人片在线一区二区| 欧美乱码精品一区二区三区| 女人高潮潮喷娇喘18禁视频| 在线免费观看的www视频| 欧美午夜高清在线| 最新美女视频免费是黄的| 多毛熟女@视频| 精品午夜福利视频在线观看一区| 首页视频小说图片口味搜索| 欧美精品人与动牲交sv欧美| 午夜精品久久久久久毛片777| 麻豆av在线久日| 天天操日日干夜夜撸| 99re在线观看精品视频| 亚洲精品久久成人aⅴ小说| 黑人巨大精品欧美一区二区蜜桃| 叶爱在线成人免费视频播放| 国产精品欧美亚洲77777| 精品福利观看| 亚洲熟妇熟女久久| 两人在一起打扑克的视频| 精品久久久精品久久久| 老熟女久久久| 亚洲av熟女| 欧美中文综合在线视频| 不卡一级毛片| 精品国产乱码久久久久久男人| 777米奇影视久久| 757午夜福利合集在线观看| 一a级毛片在线观看| 老鸭窝网址在线观看| 伦理电影免费视频| 久久香蕉精品热| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到| 国产99久久九九免费精品| 午夜精品国产一区二区电影| 亚洲av第一区精品v没综合| 美女扒开内裤让男人捅视频| 妹子高潮喷水视频| 国产精品自产拍在线观看55亚洲 | 亚洲av成人不卡在线观看播放网| 最近最新中文字幕大全免费视频| 国产高清videossex| 一区二区三区激情视频| 国产精品一区二区在线观看99| 两人在一起打扑克的视频| 亚洲va日本ⅴa欧美va伊人久久| 国产主播在线观看一区二区| 亚洲人成电影免费在线| 国产成+人综合+亚洲专区| 一级a爱片免费观看的视频| cao死你这个sao货| 50天的宝宝边吃奶边哭怎么回事| 女人高潮潮喷娇喘18禁视频| 国产成人影院久久av| 亚洲在线自拍视频| 91成人精品电影| 久久久久久久久久久久大奶| 久久青草综合色| 99国产极品粉嫩在线观看| 亚洲七黄色美女视频| 午夜久久久在线观看| 国产亚洲av高清不卡| 欧美日本中文国产一区发布| 久久久久视频综合| 丰满饥渴人妻一区二区三| 久久久久精品人妻al黑| 自拍欧美九色日韩亚洲蝌蚪91| 一二三四在线观看免费中文在| 欧美激情 高清一区二区三区| 成人黄色视频免费在线看| 真人做人爱边吃奶动态| 久久精品aⅴ一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 国产精品免费一区二区三区在线 | 一边摸一边做爽爽视频免费| 亚洲精品美女久久av网站| 日韩一卡2卡3卡4卡2021年| 欧美精品啪啪一区二区三区| 亚洲熟妇熟女久久| 国产成人啪精品午夜网站| 一区二区三区激情视频| а√天堂www在线а√下载 | 久久人妻福利社区极品人妻图片| 久久影院123| 韩国精品一区二区三区| a级毛片在线看网站| 精品一区二区三区av网在线观看| 啦啦啦免费观看视频1| 嫁个100分男人电影在线观看| 久久 成人 亚洲| 色在线成人网| 久久精品91无色码中文字幕| 一区二区三区国产精品乱码| av视频免费观看在线观看| 国产熟女午夜一区二区三区| 又黄又爽又免费观看的视频| 19禁男女啪啪无遮挡网站| 女人精品久久久久毛片| 老汉色∧v一级毛片| 日韩欧美免费精品| 亚洲在线自拍视频| 久久精品熟女亚洲av麻豆精品| 黄片大片在线免费观看| 一本大道久久a久久精品| 少妇被粗大的猛进出69影院| 国产亚洲欧美98| 欧美激情久久久久久爽电影 | 精品国产美女av久久久久小说| 男人操女人黄网站| 美女国产高潮福利片在线看| videosex国产| 久久人人爽av亚洲精品天堂| 天天影视国产精品| 午夜精品在线福利| 欧美日韩精品网址| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区 | 日韩有码中文字幕| 在线观看66精品国产| 国产午夜精品久久久久久| 中文字幕色久视频| 嫩草影视91久久| av超薄肉色丝袜交足视频| 国产99久久九九免费精品| 巨乳人妻的诱惑在线观看| 亚洲五月婷婷丁香| 亚洲国产精品合色在线| 亚洲成av片中文字幕在线观看| 18禁裸乳无遮挡动漫免费视频| 天堂俺去俺来也www色官网| 亚洲人成77777在线视频| 欧洲精品卡2卡3卡4卡5卡区| 新久久久久国产一级毛片| 午夜福利在线观看吧| 久久 成人 亚洲| 一区福利在线观看| 久久精品成人免费网站| 久久午夜综合久久蜜桃| 亚洲av片天天在线观看| 不卡av一区二区三区| 在线观看一区二区三区激情| 三级毛片av免费| 激情在线观看视频在线高清 | 亚洲专区字幕在线| 日韩三级视频一区二区三区| 国产aⅴ精品一区二区三区波| 国产野战对白在线观看| 99精品久久久久人妻精品| www日本在线高清视频| 欧美黑人欧美精品刺激| 午夜久久久在线观看| 成人国产一区最新在线观看| 视频在线观看一区二区三区| 亚洲专区字幕在线| 午夜福利免费观看在线| 久久ye,这里只有精品| 老熟妇乱子伦视频在线观看| 国产精品.久久久| www.精华液| 999久久久国产精品视频| 电影成人av| 国产成人精品久久二区二区免费| 精品人妻熟女毛片av久久网站| 自拍欧美九色日韩亚洲蝌蚪91| 老熟妇仑乱视频hdxx| 国产成人精品久久二区二区91| 久久 成人 亚洲| 另类亚洲欧美激情| 成熟少妇高潮喷水视频| 精品久久久久久,| 欧美人与性动交α欧美精品济南到| 手机成人av网站| 国产三级黄色录像| 黑人欧美特级aaaaaa片| 一进一出好大好爽视频| 成熟少妇高潮喷水视频| 亚洲国产精品sss在线观看 | 美女扒开内裤让男人捅视频| 精品久久久久久,| 在线观看日韩欧美| 不卡一级毛片| 嫁个100分男人电影在线观看| av不卡在线播放| 啦啦啦 在线观看视频| 国产精品久久电影中文字幕 | 后天国语完整版免费观看| 手机成人av网站| 成在线人永久免费视频| 91精品国产国语对白视频| 亚洲成人国产一区在线观看| 国产精品永久免费网站| 巨乳人妻的诱惑在线观看| 国产成人一区二区三区免费视频网站| 国产无遮挡羞羞视频在线观看| 久久精品国产亚洲av高清一级| 99国产精品99久久久久| www.999成人在线观看| www.精华液| 亚洲成人免费av在线播放| 夜夜夜夜夜久久久久| 热re99久久精品国产66热6| 狠狠狠狠99中文字幕| 亚洲精品av麻豆狂野| 国产不卡av网站在线观看| 久久久久久亚洲精品国产蜜桃av| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 不卡一级毛片| 成人手机av| 老司机在亚洲福利影院| 人成视频在线观看免费观看| 不卡av一区二区三区| 久久精品人人爽人人爽视色| 中文字幕av电影在线播放| 国产主播在线观看一区二区| 亚洲精品国产一区二区精华液| 男女午夜视频在线观看| 久久久精品国产亚洲av高清涩受| 三上悠亚av全集在线观看| 一区二区三区精品91| 曰老女人黄片| 久久精品国产亚洲av高清一级| 亚洲国产欧美网| 亚洲成人国产一区在线观看| 国产人伦9x9x在线观看| 精品国产美女av久久久久小说| 精品午夜福利视频在线观看一区| 欧美精品高潮呻吟av久久| 精品一区二区三区四区五区乱码| 韩国精品一区二区三区| 无人区码免费观看不卡| 久久人人97超碰香蕉20202| a在线观看视频网站| 欧美日韩瑟瑟在线播放| 亚洲av日韩在线播放| 国产成人精品久久二区二区免费| 久久精品人人爽人人爽视色| 成人精品一区二区免费| 久久久久久久午夜电影 | 日本黄色视频三级网站网址 | 成人国产一区最新在线观看| 国产精品美女特级片免费视频播放器 | 色尼玛亚洲综合影院| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品亚洲av一区麻豆| av线在线观看网站| 日韩一卡2卡3卡4卡2021年| 国产欧美日韩一区二区三区在线| 99香蕉大伊视频| 精品国产超薄肉色丝袜足j| 久久久久久久久久久久大奶| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av第一区精品v没综合| 女人爽到高潮嗷嗷叫在线视频| 久久精品人人爽人人爽视色| 亚洲精品av麻豆狂野| 十八禁网站免费在线| 国产精品一区二区在线不卡| 亚洲视频免费观看视频| 成人18禁高潮啪啪吃奶动态图| 久久亚洲精品不卡| 法律面前人人平等表现在哪些方面| 波多野结衣av一区二区av| 日本撒尿小便嘘嘘汇集6| 黑人巨大精品欧美一区二区蜜桃| 极品人妻少妇av视频| 欧美+亚洲+日韩+国产| 久9热在线精品视频| 热re99久久精品国产66热6| 欧美日韩成人在线一区二区| 黄网站色视频无遮挡免费观看| 免费人成视频x8x8入口观看| 国产乱人伦免费视频| 天天操日日干夜夜撸| 亚洲国产精品合色在线| 久久亚洲真实| 国产亚洲av高清不卡| 久久精品91无色码中文字幕| 一边摸一边做爽爽视频免费| 成人手机av| 精品一区二区三区视频在线观看免费 | 老司机影院毛片| 免费在线观看日本一区| 岛国毛片在线播放| 色综合婷婷激情| 露出奶头的视频| 真人做人爱边吃奶动态| 91成年电影在线观看| 美女福利国产在线| 国产三级黄色录像| 午夜91福利影院| 国产免费现黄频在线看| 亚洲av成人一区二区三| 黑丝袜美女国产一区| 99久久99久久久精品蜜桃| 国产精品一区二区在线观看99| 色综合婷婷激情| 男女床上黄色一级片免费看| 中文字幕人妻丝袜制服| 99国产精品一区二区三区| 精品免费久久久久久久清纯 | 久久精品成人免费网站| 丰满饥渴人妻一区二区三| 满18在线观看网站| 欧美最黄视频在线播放免费 | 久久久久国产精品人妻aⅴ院 | 人妻久久中文字幕网| 国产成人免费无遮挡视频| 国产男靠女视频免费网站| videosex国产| 亚洲熟女精品中文字幕| 看免费av毛片| 国产无遮挡羞羞视频在线观看| 日本精品一区二区三区蜜桃| 欧美日韩福利视频一区二区| xxx96com| 又黄又爽又免费观看的视频| 91成年电影在线观看| 免费在线观看日本一区| 精品乱码久久久久久99久播| 久久久久精品国产欧美久久久| 人妻 亚洲 视频| 正在播放国产对白刺激| 99国产精品免费福利视频| 我的亚洲天堂| 国产精品九九99| 午夜日韩欧美国产| 欧美激情高清一区二区三区| 色婷婷久久久亚洲欧美| 亚洲av欧美aⅴ国产| 久热爱精品视频在线9| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区视频在线观看免费 | 国产成人精品无人区| 亚洲国产精品合色在线| 这个男人来自地球电影免费观看| 婷婷成人精品国产| 久久这里只有精品19| av中文乱码字幕在线| av片东京热男人的天堂| 久久狼人影院| 两个人看的免费小视频| 岛国毛片在线播放| 国产在视频线精品| 亚洲一码二码三码区别大吗| 变态另类成人亚洲欧美熟女 | av超薄肉色丝袜交足视频| videos熟女内射| av片东京热男人的天堂| 美女午夜性视频免费| 夫妻午夜视频| 99热只有精品国产| 在线观看免费高清a一片| 国产av精品麻豆| 久久精品亚洲熟妇少妇任你| 亚洲一区二区三区欧美精品| 亚洲精品在线美女| 18禁黄网站禁片午夜丰满| 国产亚洲精品第一综合不卡| 午夜久久久在线观看| 国产91精品成人一区二区三区| 欧美久久黑人一区二区| av有码第一页| 免费日韩欧美在线观看| 欧美亚洲 丝袜 人妻 在线| 飞空精品影院首页| 色婷婷av一区二区三区视频| 欧美在线黄色| 中文字幕最新亚洲高清| 国产精品.久久久| 午夜福利影视在线免费观看| 国产精品 国内视频| 精品国产亚洲在线| 一进一出抽搐动态| 国产又爽黄色视频| av天堂久久9| 成年人免费黄色播放视频| 久久精品国产综合久久久| 在线看a的网站| 一二三四在线观看免费中文在| 中文欧美无线码| 丝瓜视频免费看黄片| 久久草成人影院| 久久久精品免费免费高清| 我的亚洲天堂| 叶爱在线成人免费视频播放| 亚洲精品一卡2卡三卡4卡5卡| 如日韩欧美国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 色综合欧美亚洲国产小说| 国产不卡av网站在线观看| 国精品久久久久久国模美| 黑丝袜美女国产一区| 9热在线视频观看99| 欧美乱妇无乱码| 人人妻,人人澡人人爽秒播| 精品久久蜜臀av无| 91麻豆精品激情在线观看国产 | 乱人伦中国视频| 成熟少妇高潮喷水视频| 亚洲久久久国产精品| 一区二区三区国产精品乱码| 亚洲欧美一区二区三区黑人| 天堂俺去俺来也www色官网| 无人区码免费观看不卡| avwww免费| 人人妻人人添人人爽欧美一区卜| 国产精品自产拍在线观看55亚洲 | 欧美日韩视频精品一区| 色在线成人网| 久久影院123| 成人亚洲精品一区在线观看| 免费日韩欧美在线观看| 如日韩欧美国产精品一区二区三区| 婷婷精品国产亚洲av在线 | 成年人免费黄色播放视频| 美国免费a级毛片| 免费人成视频x8x8入口观看| 女人被躁到高潮嗷嗷叫费观| 777米奇影视久久| 免费在线观看完整版高清| 欧美色视频一区免费| 午夜免费成人在线视频| 午夜成年电影在线免费观看| 捣出白浆h1v1| 免费日韩欧美在线观看| 国产av又大| 久久国产乱子伦精品免费另类| √禁漫天堂资源中文www| 亚洲久久久国产精品| 日韩免费高清中文字幕av| 亚洲一区二区三区欧美精品| 日韩欧美国产一区二区入口| 国产精品 欧美亚洲| 岛国毛片在线播放| 一边摸一边抽搐一进一出视频| 男人操女人黄网站| 亚洲第一av免费看| 精品国产乱子伦一区二区三区| 露出奶头的视频| 水蜜桃什么品种好| bbb黄色大片| 在线观看免费视频网站a站| 50天的宝宝边吃奶边哭怎么回事| 狠狠婷婷综合久久久久久88av| 午夜两性在线视频| e午夜精品久久久久久久| 热99久久久久精品小说推荐| 精品卡一卡二卡四卡免费| 亚洲精华国产精华精| 老熟妇仑乱视频hdxx| 熟女少妇亚洲综合色aaa.| 日韩三级视频一区二区三区| 一边摸一边抽搐一进一出视频| 久久久国产成人免费| 国产真人三级小视频在线观看| 午夜免费鲁丝| 久久人妻熟女aⅴ| 亚洲中文字幕日韩| av片东京热男人的天堂| 国产精品免费一区二区三区在线 | 亚洲精品粉嫩美女一区| 激情在线观看视频在线高清 | 99re在线观看精品视频| 窝窝影院91人妻| 人成视频在线观看免费观看| 亚洲专区中文字幕在线| 黑人猛操日本美女一级片| 国产精品九九99| 午夜激情av网站| 国产成人一区二区三区免费视频网站| 国产精品久久久人人做人人爽| 中文字幕高清在线视频| 久久精品91无色码中文字幕| 婷婷精品国产亚洲av在线 | 老司机影院毛片| 黑人巨大精品欧美一区二区mp4| 日韩一卡2卡3卡4卡2021年| 国产91精品成人一区二区三区| 精品久久久精品久久久| 久久九九热精品免费| 国产精品久久久av美女十八| 50天的宝宝边吃奶边哭怎么回事| 久久99一区二区三区| 午夜精品国产一区二区电影| 色尼玛亚洲综合影院| 亚洲性夜色夜夜综合| 丰满饥渴人妻一区二区三| 中文字幕色久视频| 啦啦啦 在线观看视频| 97人妻天天添夜夜摸| 久久这里只有精品19| a级片在线免费高清观看视频| 咕卡用的链子| 成年人黄色毛片网站| 亚洲国产毛片av蜜桃av| 国精品久久久久久国模美| a级毛片在线看网站| 午夜免费成人在线视频| 精品久久久久久久久久免费视频 | 午夜精品在线福利| 国产又色又爽无遮挡免费看| 看黄色毛片网站| 大型av网站在线播放| 亚洲欧美日韩另类电影网站| 亚洲熟女精品中文字幕| 日韩欧美国产一区二区入口| 极品人妻少妇av视频| 波多野结衣av一区二区av| 法律面前人人平等表现在哪些方面| 国产一卡二卡三卡精品| 亚洲精品久久午夜乱码| 女性生殖器流出的白浆| 一本大道久久a久久精品| 一级片免费观看大全| 建设人人有责人人尽责人人享有的| 日韩有码中文字幕| tube8黄色片| 亚洲精品国产色婷婷电影| 日本精品一区二区三区蜜桃| 色综合婷婷激情| 高清黄色对白视频在线免费看| 中文亚洲av片在线观看爽 | 18禁观看日本| 91麻豆av在线| 亚洲av美国av| 成人精品一区二区免费| 久久精品国产亚洲av高清一级| 国产亚洲av高清不卡| 免费在线观看黄色视频的| 无遮挡黄片免费观看| 伦理电影免费视频| av有码第一页| 日本五十路高清| 午夜久久久在线观看| 99在线人妻在线中文字幕 | 一级作爱视频免费观看| 一进一出好大好爽视频| 欧美性长视频在线观看| 一级作爱视频免费观看| 1024香蕉在线观看| 欧美激情久久久久久爽电影 | 久99久视频精品免费|