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

    壓裝密實炸藥裝藥非沖擊點火反應傳播與烈度演化實驗研究進展*

    2020-04-09 11:14:30胡海波尚海林文尚剛
    爆炸與沖擊 2020年1期
    關鍵詞:柱殼烈度裝藥

    胡海波,傅 華,李 濤,尚海林,文尚剛

    (中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點實驗室,四川 綿陽 621999)

    1 概 述

    本文中關注的問題對象,是典型事故載荷作用條件下,以低孔隙率壓裝PBX 為代表的密實炸藥在特定裝藥約束結(jié)構(gòu)中的非沖擊點火反應傳播和演化行為。所謂非沖擊點火[1],泛指除高速碎片撞擊和射流侵徹之外,各類不經(jīng)由沖擊波直接引發(fā)爆轟的炸藥反應行為起始,其不存在固定演化模式和定常狀態(tài)特征,是典型事故反應研究中可觀測的顯性狀態(tài)突變起點。作為對事故反應基本認識進展的概要說明,本文暫避開炸藥非沖擊點火反應起源這個專門的問題,集中討論裝藥結(jié)構(gòu)中炸藥在以某種非沖擊方式點火后,以炸藥表面層流燃燒、產(chǎn)物沿縫隙流動及裝藥結(jié)構(gòu)變形響應耦合為特征機制的炸藥事故反應傳播、烈度演化的基本規(guī)律[2]。

    基于對非沖擊點火反應主導機制的認識,重新審視了傳統(tǒng)實驗設計理念,有針對性地改進診斷方法,通過實驗取得了更詳實展現(xiàn)對象過程本來面目的觀測結(jié)果。相關實驗中裝藥縫隙為殼體與炸藥的機加公差間隙,均未采用膠粘或封堵。作為代表案例,將介紹壓裝PBX 炸藥件之間的間隙中,燃燒產(chǎn)物對流傳播和烈度演化的物理行為特征[3];對毫米級薄壁柱殼弱約束下,分段藥柱在一端燃燒點火后的反應傳播、烈度演化,呈現(xiàn)整體爆燃的行為機理進行了解讀闡釋[4-5];對傳統(tǒng)DDT 實驗中20 mm 厚壁、大長度鋼質(zhì)柱殼約束下,壓裝PBX 分段藥柱爆炸級事故反應演化及爆燃轉(zhuǎn)爆轟的機理可行性,結(jié)合實測數(shù)據(jù)進行了澄清[6];細致跟蹤了20 mm 中等壁厚鋼球形殼約束下,密實炸藥球在從中心燃燒點火后,反應產(chǎn)物驅(qū)動動態(tài)裂紋傳播、多級分叉演化過程及晚期烈度的陡增表現(xiàn)[7];分析了一端開放的50 mm 大壁厚鋼柱殼在可自由滑動活塞質(zhì)量慣性約束下,從表面引燃的炸藥柱模擬事故反應呈現(xiàn)典型爆炸表現(xiàn)的成因[8];重點討論了在有限強度約束條件下,經(jīng)歷猛烈撞擊預破碎,或事故反應過程中高壓驅(qū)動下快速動態(tài)變形、破碎炸藥的反應行為特征[9],其可造成高超壓、高速碎片飛散的劇烈爆炸,但反應機理不能簡單歸類于因沖擊波引發(fā)的爆轟。試圖通過準確解讀這些實驗紛繁復雜的事故反應表象后所蘊含的物理規(guī)律,進一步澄清國內(nèi)外同行長期以來在炸藥事故反應行為機制解讀和模型描述中的基本概念誤區(qū)[10-14]。

    2 炸藥事故反應基本模式及其相互關聯(lián)

    對以壓裝PBX 為代表的密實炸藥,其作為含能材料可能發(fā)生的化學反應的最基本形式不外乎以下三種[15]:

    (1)分解反應,其涵蓋老化及溫升條件下炸藥基體中各部位炸藥晶體的持續(xù)緩慢分解,不存在反應波陣面,各部位炸藥反應速率取決于當?shù)販囟?,通常借助阿倫紐斯律表征;

    (2)燃燒反應,指熱傳導機制主導、以亞聲速向炸藥基體中傳播并轉(zhuǎn)化為氣態(tài)產(chǎn)物的反應陣面,故稱熱傳導燃燒、層流燃燒,高壓下的層流燃燒反應也被稱為爆燃,其燃速即反應陣面?zhèn)鞑ニ俣戎饕蕾嚥嚸嫔系姆磻a(chǎn)物壓力[16-17],而氣態(tài)反應產(chǎn)物壓力主要取決于其周邊結(jié)構(gòu)的約束狀態(tài);

    (3)爆轟,由高強度沖擊波引導、在炸藥基體中以超聲速傳播的化學反應,炸藥在跨越波陣面瞬間被壓縮并完成反應,高壓反應產(chǎn)物以亞聲速逆向膨脹飛散,爆轟定常傳播的速度和壓力與結(jié)構(gòu)約束原則上無關,僅取決于炸藥基體密度[18]。

    在安全性研究中,無波陣面形態(tài)的慢速分解反應被歸類為事故點火前反應,其生成的微量氣態(tài)產(chǎn)物以分散形態(tài)留存在炸藥基體微孔洞中,除引發(fā)溫升、體脹外尚未匯聚形成宏觀狀態(tài)分區(qū)效應。與分解反應劃界的關鍵狀態(tài)轉(zhuǎn)折點,即所謂事故點火被定義為炸藥特定部位因事故作用產(chǎn)生的極高溫度狀態(tài),宏觀上開始爆發(fā)性產(chǎn)生大量氣態(tài)產(chǎn)物[19-20],使反應進入有波陣面的燃燒,形成與炸藥基體狀態(tài)不同的氣態(tài)產(chǎn)物高溫區(qū)域的變化。從嚴格意義上講,由沖擊波作用引發(fā)的爆轟自然也可以歸入此定義,但其屬于極端事故刺激下,吉帕級高強度沖擊波作用、微秒級特征時間瞬時響應點火的特別案例;大多數(shù)事故條件下(如持續(xù)升溫至特定幅值后溫度突變爆發(fā)或劇烈撞擊條件下因局域化能量沉積生成極端高溫狀態(tài))的非沖擊點火,專指引發(fā)炸藥層流燃燒,但不直接形成爆轟的其他類別點火方式,其對應的響應也與微秒級時程特征的事故加載局域應力無綁定關系,除對極度敏感的始發(fā)藥外,這類事故初始階段反應所產(chǎn)生的應力波在各種主炸藥,尤其是典型PBX 炸藥中,不足以快速引發(fā)爆轟,事故反應呈現(xiàn)與結(jié)構(gòu)響應行為和炸藥變形、破壞相關的多階段反應接續(xù)演化特征,故被稱為非瞬時反應。對低孔隙率密實炸藥而言,在排除了燃燒陣面滲透進入炸藥基體中微介觀孔隙的經(jīng)典定義下的對流燃燒后,其非沖擊點火基本反應模式只可能是表面層流燃燒,具體涉及燃燒波陣面向炸藥基體中推進燃耗,以及氣態(tài)反應產(chǎn)物在密閉條件下的成壓和流動遷移。這是系統(tǒng)解讀低孔隙率密實炸藥非沖擊點火反應行為特性的基本源頭要素。在炸藥事故反應形成和發(fā)展過程中,上述三種基本反應模式間的相互轉(zhuǎn)化是所有研究者所關注和討論的基本內(nèi)容。

    需特別提醒注意的是,蘇聯(lián)、俄羅斯同行關于炸藥層流燃燒、對流燃燒、低速爆轟、爆轟四種基本反應形式的定義及其相互關系闡釋,尤其是其涉及對流燃燒、似爆轟的爆炸、低速爆轟、爆炸轉(zhuǎn)化等基本術(shù)語表述[21],源于相關學者早期對凝聚相多孔炸藥反應行為特征的現(xiàn)象學概括[22],若引鑒來概括所有種類炸藥反應行為,容易引起對這幾種反應型式注定會遞次轉(zhuǎn)化引發(fā)爆炸和爆轟的誤導,偏離原書作者的本義。長期以來,蘇聯(lián)、俄羅斯同行關于炸藥層流燃燒的研究非常系統(tǒng)、深入[16],但關聯(lián)至密實炸藥裝藥結(jié)構(gòu)中的事故反應演化機理解讀[23-24]仍不完全到位。

    此外,需要提醒注意應加以區(qū)別的是,在炸藥安全性研究領域,描述炸藥基本反應行為的術(shù)語,如燃燒、爆燃、爆轟[19-20],與彈藥事故反應烈度評價的術(shù)語,如燃燒、爆燃、爆炸、爆轟[25]的外在形式相同。前者針對作為含能材料的炸藥,用于描述特征反應事件在反應陣面上的瞬時進程狀態(tài),對應特定反應速率物理-數(shù)學模型;后者針對裝藥結(jié)構(gòu),是對作為彈藥的炸藥裝藥系統(tǒng)反應過程特征及整體后果的概括描述,側(cè)重表征后果輸出即毀傷效應分級。

    本文中的討論主體上側(cè)重從事故反應機理、反應演化行為跟蹤和對應過程物理-數(shù)學模型描述的角度來展開,此時使用的相關術(shù)語原則上屬于第一類定義;但在涉及討論事故反應烈度內(nèi)容時,自然應用第二類定義的烈度分級概念。在美國LANL、LLNL 和SNL 三大實驗室,對第二類定義中爆燃以上烈度的反應又被統(tǒng)稱為劇烈反應,包含前面分級中爆炸及爆轟,其中爆炸涵蓋從微量炸藥反應、不形成顯著超壓和危險碎片的有限烈度反應以外的各類相對劇烈反應,到定常爆轟,狀態(tài)跨度極大,通常還有系列細化分等,如輕度爆炸、劇烈爆炸等,但無嚴格的量化判別標準。

    從本質(zhì)上講,不同形態(tài)的爆炸只不過是炸藥層流燃燒反應,因炸藥基體在事故載荷下預破碎或在反應過程中動態(tài)破碎,形成網(wǎng)狀裂紋使燃燒面積劇增,在結(jié)構(gòu)和慣性約束下達到數(shù)百兆帕、甚至吉帕級水平的高壓,進入快速、高燃耗率反應狀態(tài),所呈現(xiàn)的綜合效應表現(xiàn)。這里所謂爆炸,并非炸藥事故反應的基本形式,其沒有嚴格的學術(shù)概念定義[26],也沒有相對固定的反應速率模型形式。在日常口語中,爆炸一詞還被用來泛指一切大量產(chǎn)氣、可形成結(jié)構(gòu)破壞、拋擲的劇烈過程,包含事故反應烈度分級中的爆炸和爆轟,以及槍彈發(fā)射時發(fā)射藥在膛內(nèi)的受控燃燒[18,26]。

    從非沖擊點火事故反應并非沖擊波引導的爆轟反應的共性看,無論是溫升烤燃,還是撞擊點火,裝藥結(jié)構(gòu)中事故反應的物理過程特征及反應行為演化描述的控制方程和物理-數(shù)學模型機理完全一致,基本落腳到炸藥燃燒區(qū)表面層流燃燒行為描述,以及高溫高壓產(chǎn)物在結(jié)構(gòu)間隙、炸藥基體裂紋中的輸運和熱傳導問題的模型描述上。其從根本上區(qū)別于描述爆轟反應點火、傳播和驅(qū)動效應的模型,無法從爆轟模型修正、派生,也不存在因具體事故類型和點火起源歷程特征而建立不同模型體系的需求。所不同的,僅僅是求解控制方程時,炸藥的初、邊界條件有所不同,如炸藥基體經(jīng)歷了特定歷程烤燃的溫度分布(炸藥基體力學性能、孔隙率會因此產(chǎn)生對應變化)或撞擊產(chǎn)生的斷裂等。在此大前提下,典型裝藥結(jié)構(gòu)中的事故反應傳播和演化,被納入與密實炸藥具體部位燃燒、斷裂發(fā)展及裝藥結(jié)構(gòu)變形響應的進程中來考慮。

    3 對炸藥非沖擊點火反應的基本認識

    在安全性研究起步關注的裝藥系統(tǒng)是否點火及如何點火的第一類基本問題外,本文中討論的基本視線將聚焦在炸藥非沖擊點火事件發(fā)生后,裝藥結(jié)構(gòu)約束以及炸藥基體質(zhì)量慣性約束狀態(tài)下,事故反應傳播及事故反應烈度演化問題上。其所涉及的幾方面的最基本物理行為機制,在早期的安全性實驗現(xiàn)象解讀和物理建模中普遍被忽略,潛意識地被借用于爆轟沖擊點火的概念和模型(認定反應主要發(fā)生在基體中,而不是亞音速層流燃燒波陣面上),或僅為方便數(shù)學處理而引入的假定前提(如忽略間隙和裂紋演化,忽略間隙中氣態(tài)區(qū)高溫與基體低溫的反差),而使得真實過程圖像得不到正確解讀和對應描述。因為漏掉了真實的主導機制和敏感因素,眾多學者的物理-數(shù)學建模嘗試思路與實際物理過程背道而馳,所依托的相關實驗校驗也多經(jīng)不起推敲,許多研究課題多年努力無功而返,甚至形成誤導[10-14]。在以層流燃燒為核心關注點的固體推進劑發(fā)動機工作狀態(tài)設計和安全性研究中,這些前提屬于常識[27-28],其在炸藥安全性研究領域曾被提出探討[29-31],但始終未得到深入研究、形成共識[32],未明確將其作為根本區(qū)別于爆轟行為的另外一類基本作用機制,嘗試建立獨立的模型框架體系。作為后續(xù)討論的基本出發(fā)點,也是化爆安全性問題研究物理-數(shù)學建模的基本立足點,本文中將幾個方面相互耦合并存的物理機制前提重新梳理加以表述[2],其適應于始發(fā)藥外絕大多數(shù)主炸藥類型的低孔隙率密實炸藥,尤其是PBX 主炸藥的安全性問題分析:

    (1)低幅值應力波,尤其是因燃燒引發(fā)或經(jīng)歷裝藥結(jié)構(gòu)調(diào)制到達炸藥的碰撞事故載荷前沿增長緩慢,無法引發(fā)密實炸藥基體中化學反應。即事故載荷無法激活炸藥基體中微孔洞作為基本起源熱點(典型PBX 主炸藥沖擊起爆熱點激活點火條件:~GPa、<μs)。碰撞類事故點火不可能由僅引發(fā)炸藥基體幾攝氏度、幾十攝氏度平均溫升的應力波本身引發(fā),而必須借助事故載荷在裝藥結(jié)構(gòu)中特定部位上非均勻能量沉積(結(jié)構(gòu)穿刺、縫隙擠壓,破碎、粉化及再壓實,達到數(shù)百甚至上千度瞬時溫升)而引發(fā)。

    (2)典型非沖擊點火類事故載荷所引發(fā)的燃燒反應通常只會在炸藥的表面進行。即便最初的點火發(fā)生在炸藥內(nèi)部,如慢烤過程中炸藥內(nèi)部相對高溫區(qū)中,因小幅度溫度漲落或其他種類炸藥成分夾雜造成的個別部位率先溫升失控點火爆發(fā),會因該點火起源點氣態(tài)產(chǎn)物大量生成,引發(fā)高壓使炸藥基體斷裂,而轉(zhuǎn)化為炸藥縫隙表面燃燒。此外,密實炸藥中初始孔隙率極低,高壓下還將進一步預壓縮變密實,通常在粉體炸藥或高孔隙率炸藥中典型的微介觀孔隙對流燃燒基本不會發(fā)生,更難形成主導機制。在密實炸藥安全性研究中,所簡稱的對流燃燒,嚴格講是指裝藥宏觀間隙和裂紋中反應產(chǎn)物的長程流動物質(zhì)流動及其后續(xù)燃燒、結(jié)構(gòu)響應行為,應始終牢記其與傳統(tǒng)對流燃燒定義的概念差別。對應不同壓力,炸藥表面層流燃燒有其對應尺度特征的波陣面結(jié)構(gòu),來維持能量釋放和向反應區(qū)的反饋。常壓下燃燒區(qū)典型層狀結(jié)構(gòu)為數(shù)毫米,高壓下降到亞毫米。為維系其準穩(wěn)態(tài)自持反應特征結(jié)構(gòu),微尺度縫隙中燃燒反應注定只能在高壓下進行。約束結(jié)構(gòu)破裂引發(fā)的急劇壓力下降,會使高壓下小尺度縫隙中燃燒波陣面結(jié)構(gòu)失穩(wěn)熄滅,這是爆燃類事故反應中部分反應的炸藥碎塊能夠被回收的原因。這類碎塊表面會留有燃燒反應痕跡,而基體中保持無反應狀態(tài)。

    (3)裝藥結(jié)構(gòu)中氣態(tài)產(chǎn)物的生成和流動是裝藥結(jié)構(gòu)事故響應的驅(qū)動力。從點火瞬間起,就開始有大量氣態(tài)高溫產(chǎn)物(>2 000 K)快速生成,其所占據(jù)的體積將遠大于反應前對應的固相部分體積。事故反應氣態(tài)高溫產(chǎn)物在密閉裝藥結(jié)構(gòu)中形成的壓力和流動去向,是主導任何一個裝藥系統(tǒng)事故反應演化的動力學因素。在化學反應無法借助應力波壓縮機制傳播的前提下,氣態(tài)高溫產(chǎn)物在裝藥結(jié)構(gòu)間隙、炸藥縫隙中的流動,即宏觀結(jié)構(gòu)中的對流,是燃燒反應快捷傳播的基本方式。炸藥基體燃耗層流燃速通常僅在1 mm/s 量級,但沿縫隙傳播的產(chǎn)物對流速度可達到超聲速,甚至超過爆速!

    (4)反應壓力驅(qū)動的炸藥基體變形、斷裂及燃燒表面積增長。氣態(tài)高溫產(chǎn)物進入結(jié)構(gòu)間隙或炸藥基體中裂紋,在縫隙壁面被引燃后的炸藥表面燃燒過程中,因反應產(chǎn)物無法從大長度的狹小縫隙快速流出導致?lián)砣?,引發(fā)高壓使燃速提升;同時,急劇增長的壓力將造成炸藥基體快速變形、破裂,不斷分叉形成枝狀裂紋,經(jīng)高溫產(chǎn)物對流火焰進入、引燃后變成新增燃燒反應表面。裂紋中燃燒產(chǎn)物高壓驅(qū)動的基體裂紋擴展速度極高,可達~103m/s!若受試炸藥脆性特征典型,在有限變形下就呈現(xiàn)斷裂、分叉、破碎,快速形成比表面積巨大的裂紋網(wǎng)絡,將加劇烈度增長正反饋趨勢。反之,若炸藥斷裂發(fā)展不顯著,則反應一定不劇烈。前一類行為多涉及低粘接劑組分的壓裝PBX 或熔鑄類炸藥,且相近配方比例的炸藥會因斷裂特性而呈現(xiàn)明顯差別;后一類反應烈度溫和的行為則是高粘接劑組分的大多數(shù)熱固型PBX 炸藥的典型表現(xiàn)。

    (5)炸藥層流燃燒速率特性影響。低黏結(jié)劑比分PBX 炸藥燃速對壓力相對較敏感,典型壓裝PBX 常壓下燃速不高于1 mm/s,即便到1 GPa 量級壓力段,燃速也僅能達到1 m/s 量級水平,遠低于數(shù)公里每秒的爆速。在毫秒級特征事故反應時段內(nèi),其對應的炸藥燃耗僅在微米至毫米量級。脆性炸藥的事故反應中通常會伴隨大規(guī)模裂紋的產(chǎn)生,使得參與燃燒的表面積增加,反應加速趨勢會隨縫隙的參與和反應縫隙中壓力的上升而加劇。在快速形成極度碎化裂紋網(wǎng)絡并在其中形成高壓的條件下,由其巨大的反應表面積折合至(對應爆轟反應截面的)裝藥初始截面積上的等效燃速會變得十分驚人。這是典型壓裝PBX 炸藥借助燃速極低的亞音速層流燃燒事故反應就能夠造成劇烈爆炸反應的直接原因。同時,取決于結(jié)構(gòu)強度和慣性約束能力,高壓反應持續(xù)時間可能足夠長,這可為表面燃燒的炸藥碎塊盡可能的份額地完成反應創(chuàng)造條件。相對而言,粘接劑比分較多的熱固炸藥(~10%)燃速對壓力不甚敏感,且炸藥斷裂隨變形的發(fā)展通常不嚴重,事故反應烈度增長通常會受到結(jié)構(gòu)解體強度限制,對應表現(xiàn)出較低的整體反應烈度。

    (6)事故反應裝藥結(jié)構(gòu)中溫度分布狀態(tài)特征。非沖擊點火事故反應過程中,應力波作用所引發(fā)的溫升可忽略,基體溫度基本保持起始前溫度狀態(tài),但炸藥燃燒波陣面及結(jié)構(gòu)縫隙、裂紋中氣態(tài)產(chǎn)物中的溫度極高,通常>2 000 K。這種溫度分區(qū)以炸藥表面熱傳導燃燒波陣面劃界,高溫產(chǎn)物氣體會因壓差沿縫隙流動,形成包覆炸藥碎塊燃燒的高溫氣體連通網(wǎng)絡。這里所謂熱傳導燃燒,僅指燃燒波陣面結(jié)構(gòu)內(nèi),燃燒反應高溫原則上無法通過熱傳導機制進入燃燒陣面前的炸藥固相基體中,所以,事故反應過程中,總是保持固冷、氣熱的兩極化溫度分區(qū)。在各類安全性實驗中回收的殘藥碎塊上,通??吹降谋砻鏌g、拋擲碰撞新斷面上炸藥無反應并存的特征,這正是這種溫度分區(qū)作用的直接表現(xiàn)。

    (7)裝藥結(jié)構(gòu)在事故反應壓力驅(qū)動下的變形響應。典型裝藥結(jié)構(gòu)強度有限,且存在連接弱環(huán),數(shù)十至數(shù)百兆帕、以毫秒級特征時間持續(xù)的內(nèi)部壓力,足以造成結(jié)構(gòu)從薄弱環(huán)節(jié)解體。這是非沖擊點火事故反應進程中,在燃速特性外,制約烈度增長的最基本負反饋因素。若此前正在燃燒的炸藥來不及完成反應,炸藥表面的燃燒會因陡然降壓而熄滅。但若炸藥燃速對壓力敏感、炸藥基體脆性典型而使得斷裂演化嚴重,在一定結(jié)構(gòu)強度,甚至僅在炸藥基體自身慣性約束下,事故反應烈度仍有可能在足夠短的時間內(nèi)躥升,使得結(jié)構(gòu)解體前,相當部分炸藥完成反應而形成劇烈爆炸。接下來的問題,便是這類劇烈反應能否在裝藥中及時形成沖擊波引發(fā)爆轟SDT。從典型事故反應毫秒級特征反應加速行為看,形成陡前沿、高強度沖擊波的可能性不大。在炸藥極度碎化、極高燃燒表面積條件下,對較大尺度裝藥,生成沖擊波轉(zhuǎn)爆轟的風險尚不能從理論上徹底排除。但高強度沖擊波形成之前的事故載荷作用和事故反應過程注定會與結(jié)構(gòu)變形耦合,結(jié)構(gòu)響應相關的特定對象在事故反應驅(qū)動下的壓縮效應評估,不能以最終是否可能轉(zhuǎn)爆轟作為唯一評判依據(jù),應具體考慮烈度演化過程中的壓力增長歷程與對象結(jié)構(gòu)變形的耦合。

    此外,從事故反應演化過程發(fā)展的整體時空尺度看,密實炸藥事故反應演化的主導機制雖明確,但因涉及炸藥裂紋復雜結(jié)構(gòu)中流固耦合,對應多物理、多過程反應演化走向的非線性特征極端典型。前述基本前提和影響因素會以不同形式相互關聯(lián),影響事故反應烈度演化走向和整體反應釋能。其中,受試高能炸藥的燃速特性是這種走向敏感性的基本影響因素,除其與原則上不會導致劇烈反應的鈍感炸藥基炸藥燃速的差別外,組分類似的不同高能炸藥配方之間因炸藥燃速性能差別,所對應的事故反應行為表現(xiàn)可能相去甚遠。其間會自動涉及到與力學性能和斷裂行為特性耦合,而炸藥變形和斷裂行為特性又顯著依賴于結(jié)構(gòu)約束形狀和圍壓狀態(tài),如此這般形成相互耦合,反應烈度演化與全系統(tǒng)狀態(tài)形成多重關聯(lián)。炸藥縫隙中反應驅(qū)動結(jié)構(gòu)變形的制約因素,在低壓段主要取決于來自于外部結(jié)構(gòu)約束的應力,但在高壓高應變率段,炸藥自身質(zhì)量的慣性約束,也會成為主導因素。事故反應不存在特征波陣面和固定傳播模式,事故載荷強弱和初始階段反應狀態(tài)并不完全決定其后續(xù)發(fā)展走向,對給定裝藥系統(tǒng)通常存在某種形式的事故反應特征行為。

    我們注意到,非沖擊點火反應事件通常呈分區(qū)異步形態(tài),在裝藥內(nèi)部各部位并行發(fā)展:反應較緩慢時,各部位狀態(tài)可通過應力波或縫隙對流作用相互感知、影響、協(xié)調(diào);在局域反應劇烈時,個別部位的快速壓力增長、變形破碎和產(chǎn)物流動或無法傳遞關聯(lián)至裝藥結(jié)構(gòu)遠區(qū)部位。對有一定尺度的裝藥結(jié)構(gòu),個別局部率先一步發(fā)展起來、有一定偶發(fā)性質(zhì)的瞬時高烈度局域反應,如大尺度裝藥慢烤過程中單點爆發(fā)引發(fā)的局部反應高壓,不一定會發(fā)展成全域性事件,卻可通過造成局部結(jié)構(gòu)破裂而主導整體事故反應烈度特征。

    基于上述基本認識,來反觀早期國內(nèi)外同行在安全性研究和反應行為解讀、建模中的系列假定,如借用SDT 過程點火增長模型調(diào)參數(shù),將炸藥反應與壓力、溫度綁定,假定炸藥基體內(nèi)部處處均存在反應、各部位按不同進程比例完成反應的流行套路;在建模背景分析中考慮炸藥基體的損傷、斷裂,在模型體系生成時,最終僅將斷裂歸入本構(gòu)模型的強度修正上考慮,整體上忽略炸藥在裝藥結(jié)構(gòu)系統(tǒng)中的宏觀斷裂,將炸藥基體作為無斷裂連續(xù)介質(zhì)處理來避開宏觀斷裂問題的業(yè)內(nèi)通行做法;在算法處理中,將殼體內(nèi)壁與炸藥外界面綁定,使炸藥形狀類似氣體一樣與殼體內(nèi)腔綁定的隱含操作;將僅適于疏松炸藥的對流燃燒概念直接引入密實炸藥事故反應進程模型,借助壓力閾值啟動所謂對流燃燒機制的生硬模型組合;經(jīng)典燃燒轉(zhuǎn)變爆轟DDT 實驗中,借用粉體炸藥DDT 實驗設計及概念,將密實炸藥與殼體間隙中對流火焰?zhèn)鞑ヒl(fā)的電離探針導通或光纖探針發(fā)光信號視為密實炸藥基體或固體推進劑中反應波陣面?zhèn)鞑ニ俣?,而不關注實驗中回收試件形態(tài),結(jié)合應變、速度時間歷程測試反復暴露出來的不自洽問題線索,去反思、求證這種假想的基體中波陣面是否真正存在的通行做法,迄今在國內(nèi)外同行中仍十分流行[10-14,29-30]。凡此種種,看似不經(jīng)意的小前提改變,卻都嚴重違背以科學為基礎的精密物理實驗的基本要義。這類基本認知的偏離,實質(zhì)性耽誤了炸藥安全性領域研究工作的深入。

    我們注意到,在固體推進劑事故機制早期研究中,上述物理機制圖像一開始就被作為問題研究的基本出發(fā)點[27-28]。在炸藥安全性研究領域,諸多研究團隊的思維方式,因受爆轟問題研究長期傳統(tǒng)方法影響,不斷試圖從已有爆轟、流體動力學模型框架做修正,來建立涵蓋炸藥安全性問題的單一理論體系框架。盡管在同行圈中曾存在不同看法和激烈爭議,但這類觀點卻長期沒有得到重視而被淹沒[31]。所幸的是,在近期炸藥化爆安全性研究進程中,以LANL 安全性研究團隊為代表的部分同行通過一系列的精密診斷的現(xiàn)象再發(fā)現(xiàn)[33-36],使化爆安全性研究的核心關注點重新回到燃燒裂紋這類主導機制和敏感因素的研究上,同時,對通過物理-數(shù)學建模實現(xiàn)事故反應行為模擬的前景有了更清醒認識。中物院化爆安全性研究團隊密切關注同行研究認識進展,在近些年所開展的系列實驗,能為前述作為基本前提闡述的事故反應主導機制和敏感因素認識提供系統(tǒng)的實驗佐證。

    4 實驗進展案例

    下面簡要介紹本團隊近年在化爆安全性實驗研究方面的相關進展結(jié)果和認識啟發(fā)。該系列分解模擬實驗主要圍繞事故反應演化行為機制觀測,為記錄非沖擊點火反應表面燃燒、間隙中產(chǎn)物對流和結(jié)構(gòu)響應行為表現(xiàn),有針對性地進行了實驗裝置和診斷項目設計。這些實驗中除個別采用直接撞擊模擬事故點火外,統(tǒng)一采取無預留空腔條件下電點火頭加1~2 g 黑火藥在裝藥的特定位置人為點火。類似點火系統(tǒng)能在緊鄰的炸藥端面上產(chǎn)生100 MPa 左右的高溫氣體點火條件。

    4.1 炸藥間隙中反應產(chǎn)物對流及表面燃燒反應行為特性[3]

    不同強度盒狀金屬結(jié)構(gòu)約束下,一端引燃后0.1~0.2 mm 的炸藥縫隙中燃燒產(chǎn)物對流過程無法以激波管中突躍波陣面方式傳播,沿100 mm 長的縫隙上早期階段壓力分布呈現(xiàn)緩坡狀逐漸上升趨勢,經(jīng)數(shù)毫秒持續(xù)結(jié)構(gòu)腔體內(nèi)壓力增長達到數(shù)兆帕,各測點間壓力漸趨平衡;隨起始端炸藥表面燃燒起始,縫隙中對流火焰?zhèn)鞑ニ俣葟膩喴羲僮兊匠羲?,包覆炸藥縫隙表面,此時縫隙中壓力呈現(xiàn)緩慢增長;其后,縫隙表面炸藥經(jīng)數(shù)百微秒延遲被對流高溫產(chǎn)物引燃,呈現(xiàn)二次火光爆發(fā);此后,壓力從先期階段的較慢增長后迅速躥升,且因狹縫變寬、聯(lián)通,沿炸藥縫隙長度方向各測點壓力基本持平,快速同步上升;壓力躥升造成約束結(jié)構(gòu)變形、解體,其后各測點壓力接近同步迅速下降(見圖1)。隨外部結(jié)構(gòu)約束強度增加,炸藥從一端燃燒點火到整體結(jié)構(gòu)解體的響應時間從~15 ms 降至~5 ms,最大壓力峰值呈現(xiàn)從~50 MPa到~300 MPa 數(shù)倍變化;全部反應過程中,參與反應的炸藥份額僅數(shù)百分點。該組實驗呈現(xiàn)了密實炸藥在密閉結(jié)構(gòu)內(nèi),先期燃燒反應高壓產(chǎn)物氣體沿炸藥結(jié)構(gòu)縫隙傳播和反應烈度增長的特征物理圖像,縫隙中產(chǎn)物對流、裂縫表面燃燒點火和反應烈度劇烈增長、結(jié)構(gòu)變形解體時序關系清楚,回收殘藥試樣也證明炸藥基體未參與反應。作為對比的透明有機玻璃管約束下炸藥與殼壁間隙中對流火焰?zhèn)鞑嶒灡砻?,間隙中火焰?zhèn)鞑ゲ⒉粚ㄋ幓w中假想的“爆燃波陣面”,除起始端部分炸藥反應、碎裂外,回收到的炸藥柱基本完整,因被火焰包覆時間偏短,大多只有淺色熏痕[37](見圖2)。

    圖 1 從一端引燃的中強約束下炸藥裂縫中氣態(tài)產(chǎn)物流動、表面燃燒演化過程Fig. 1 Convective gaseous products flow along explosive slot and surface combust evolution in middle strength confinement

    圖 2 炸藥與有機玻璃柱殼間的縫隙中的產(chǎn)物流動情況及實驗后的殘藥狀態(tài)Fig. 2 Convective products flow along seam between HE charges and PMMA tube wall,and explosive residual

    4.2 薄壁金屬柱殼裝藥的爆燃行為[4-5]

    裝藥一端引燃后,其結(jié)構(gòu)響應變形時間大幅晚于應力波傳播時間;借助產(chǎn)物間隙對流,引發(fā)分段炸藥柱表面燃燒;柱殼膨脹過程中,燃燒的炸藥柱仍處在原位基本保持初始直徑,炸藥柱端面因燃燒產(chǎn)氣相互小幅度分離,部分炸藥柱局部破裂,反應產(chǎn)物進入引發(fā)燃燒開裂;柱殼軸向長度隨反應進程先伸長后縮短;殼體徑向膨脹平均速度僅1~10 m/s 量級,加速度趨近于零,說明反應過程中內(nèi)腔反應壓力始終與殼體約束強度大致保持平衡,即炸藥表面燃燒維持在30~50 MPa 準等壓狀態(tài),直至殼體破裂卸壓使反應終止;從破裂柱殼溢出的反應產(chǎn)物膨脹速度略低于1 km/s;大量破碎炸藥散落在附近,部分炸藥柱接近完整殘留,超壓等效反應份額僅百分點水平,呈現(xiàn)輕度爆燃特征(見圖3)。這組實驗確認了縫隙對流、炸藥表面燃燒過程中,炸藥基體中并無所謂反應波陣面?zhèn)鞑ィ徊煌贒DT 管中粉體炸藥行為,密實炸藥柱直徑并不會因內(nèi)部反應而伴隨柱殼膨脹同步變形,而是在表面燃燒過程中輕微減小。

    圖 3 薄壁柱殼約束下炸藥柱的爆燃反應行為演化Fig. 3 Deflagration evolution in thin wall tube confinement

    早期針對弱約束金屬柱殼中壓裝PBX 的DDT 實驗[38]顯示,大多數(shù)藥柱外形相對完整殘留,部分藥柱表面有燒蝕,反應并未進入炸藥基體內(nèi)部(見圖4)。這一關鍵細節(jié)是數(shù)值模擬建模和校驗的最基本證據(jù),絕不能被無視、忽略[13]。應反思壁面上電離探針或光纖探針觀測到的信號是否對應假想的炸藥基體中的反應陣面,并仔細結(jié)合柱殼變形、膨脹及壓力測試數(shù)據(jù)的時序、幅值,來解讀炸藥事故反應傳播及烈度增長的真實過程物理圖像。

    圖 4 殼壁帶測試開孔的薄壁柱殼約束炸藥爆燃反應后的典型殘藥狀態(tài)Fig. 4 Typical explosive residual after experiment in thin wall tube confinement with diagnostic drills in tube wall

    4.3 厚壁金屬柱殼裝藥模擬事故反應表現(xiàn)[6]

    實驗記錄顯示在點火頭動作后,炸藥非沖擊點火反應發(fā)展至厚壁柱殼斷裂的多階段響應行為共持續(xù)7~8 ms(而爆轟~50 μs):在起點段燃燒壓力驅(qū)動下,分段炸藥柱表面經(jīng)1 ms 量級特征時間的間隙對流火焰?zhèn)鞑c火和感應點火,進入上、下游同時進行的藥柱表面包覆燃燒階段;炸藥柱表面被引燃后,在壓力增長達到裝置設計所選取的柱殼環(huán)向箍縮屈服應力前,持續(xù)燃燒了數(shù)毫秒;在壓力超過柱殼屈服極限后,柱殼開始膨脹變形,但各測點膨脹起始順序呈現(xiàn)出與點火端起上、下游無關的隨機狀態(tài);對兩種HMX 基壓裝PBX 殼體膨脹直至柱殼斷裂,分別經(jīng)~200 μs 和~300 μs 才被緩慢加速到~150 m/s 和~200 m/s,無爆轟驅(qū)動的沖擊波前沿起跳特征,且上、下游各測點斜率,即柱殼加速度相互間很接近、且基本保持常數(shù),對應每微秒的速度增量僅~0.5 m/s、~1 m/s,扣除柱殼環(huán)向應力后按加速度折算的壓力為~0.6 GPa 和~1.2 GPa,可認定柱殼斷裂前,炸藥事故反應在對應幅值的準定常壓力下進行,而在爆轟標定實驗中,厚壁柱殼被瞬間加速至~400 m/s,再被爆轟產(chǎn)物驅(qū)動緩慢加速,爆壓~39 GPa;超壓等效反應份額30%~40%,殘藥極少,斷口上有未反應炸藥涂覆痕跡,整體反應烈度屬中輕度爆炸。采用20、460 mm 兩種藥柱長度的實驗測量結(jié)果顯示,反應烈度與藥柱長度無明顯關聯(lián)。早期的系列厚壁DDT 管實驗中[38-40],借助壁面探針反應前沿速度測量結(jié)果或個別壓力計數(shù)據(jù),聲稱觀測到低速爆轟、DDT 轉(zhuǎn)化行為、測量出受試炸藥在DDT 轉(zhuǎn)化特征距離的研究結(jié)論所依托的基本證據(jù)需要重新解讀、認識。

    圖 5 厚壁柱殼約束大長度DDT 管炸藥反應歷程特征和柱殼斷裂狀態(tài)Fig. 5 Non-shock initiation reaction evolution in thick wall confinement experiment and tube wall fracture

    4.4 厚壁球殼約束裝藥中心非沖擊點火實驗觀測[7]

    圖 6 點火頭從帶窗厚壁約束壓裝炸藥中心點火后炸藥球外表裂紋演化早期過程圖像Fig. 6 The crack evolution at the early half after center ignition (20 mm steel wall with PMMA window)

    圖 7 點火頭從帶窗20 mm 厚壁約束壓裝炸藥中心點火后殼體速度及約束結(jié)構(gòu)中壓力增長歷程Fig. 7 The pressure and velocity profiles inside confinement of 20 mm steel wall with PMMA window

    圖 8 點火頭從帶窗厚壁約束壓裝炸藥中心點火后,約束結(jié)構(gòu)中炸藥裂紋及反應演化進程示意圖Fig. 8 The crack evolution and reaction growth after center combustion ignition (20 mm steel wall with PMMA window)

    4.5 自由活塞柱強慣性約束藥柱模擬事故反應表現(xiàn)[8]

    一端開口、帶自由活塞的50 mm 厚壁柱殼設計的結(jié)構(gòu)連接強度為零,但活塞柱總質(zhì)量數(shù)倍于柱殼厚度對應的質(zhì)量。受試藥柱一端表面被引燃后,腔體內(nèi)的自由狀態(tài)活塞在數(shù)百微秒內(nèi)被100 MPa 的產(chǎn)物壓力驅(qū)動移動數(shù)毫米,反應壓力快速躥升至2~3 GPa 水平,驅(qū)動鋼質(zhì)柱殼壁膨脹斷裂;從柱殼變形和活塞運動加速度陡變特征評估,反應晚期藥柱可能產(chǎn)生了碎裂,引發(fā)可燃表面積增長;實驗后,裝置近區(qū)未找到殘藥,但厚壁柱殼斷裂面上未反應炸藥涂覆痕跡明顯;采用相對低感的壓裝PBX 炸藥也呈現(xiàn)出爆炸事件特征。

    4.6 相對弱約束下事故撞擊預碎化炸藥事故反應表現(xiàn)

    在數(shù)毫米厚的鋼殼約束栓塞試驗[41]、厚炸藥樣品Steven 試驗[9]中,炸藥反應超壓等效反應百分比隨撞擊速度增加變大,可超過等效全額爆炸當量的一半,甚至接近、超過滿當量。測得的典型點火延遲時間為數(shù)百微秒,反應持續(xù)時間在百微秒水平,反應壓力通常在數(shù)百兆帕。薄殼體栓塞實驗及Susan 試驗中,炸藥在反應點火前的預破碎是高烈度表現(xiàn)的直接誘因(見圖9)。對厚炸藥樣品Steven 試驗,一旦點火便呈現(xiàn)高烈度響應的表現(xiàn),則更可能與事故反應驅(qū)動的炸藥基體斷裂相關;相比之下,薄炸藥樣品Steven 試驗的放能份額隨加載速度略有增長,但總體份額始終較低(見圖10)。更大尺度的實驗中觀測到壓力數(shù)百MPa、持續(xù)時間近10 ms 的反應事件,超壓等效反應比例>2/3,則可能與基體嚴重破碎和持續(xù)反應過程中大多數(shù)碎裂表面點火燃燒相關。在這個過程中,結(jié)構(gòu)約束未能夠扮演反應烈度限幅主導角色。在Susan 試驗研究中,觀測到炸藥的點火發(fā)生在脆性藥柱在碰撞后大幅變形、破碎后狀態(tài)下[42]。

    圖 9 典型壓裝炸藥在碰撞類事故載荷作用下的變形、破碎特征[42-43]Fig. 9 Character of deformation and fragmentation of typical pressed explosive under impact loading[42-43]

    圖 10 不同厚度樣品炸藥在Steven 試驗中的反應放能特征差別[9]Fig. 10 The difference of energy release behavior of explosive in Steven test with varied thickness[9]

    5 討 論

    前述相關實驗中取得的主要認識進展進一步歸納如下:

    (1)金屬結(jié)構(gòu)約束下的密實炸藥裝藥非沖擊點火反應,本質(zhì)上是限制在炸藥表面上進行的亞音速層流燃燒,其反應速度亦稱燃耗、基體回退速度主要取決于反應產(chǎn)物壓力;此時,高壓產(chǎn)物氣體在裝藥結(jié)構(gòu)或炸藥基體縫隙中的對流(氣態(tài)產(chǎn)物快速遷移、輸運)扮演了最關鍵角色,早期反應部位的高溫產(chǎn)物氣體在宏觀結(jié)構(gòu)中的對流,是引發(fā)其他部位炸藥燃燒的基本機制途徑,上游產(chǎn)物熱焰通過縫隙到達新的縫隙表面后,經(jīng)特定弛豫時間,可使該部位炸藥表面升溫燃燒,當然也可能因熱傳導降溫熄滅;結(jié)構(gòu)間隙、炸藥縫隙中的這種反應傳播也被簡稱為對流燃燒,其涉及宏觀結(jié)構(gòu)中反應產(chǎn)物長程流動和炸藥表面燃燒兩個階段,應注意其不同于早期同行針對多孔或粉體炸藥定義的微孔隙中的對流燃燒,密實炸藥的可透過性原則上可忽略不計;

    (2)典型裝藥結(jié)構(gòu)中的非沖擊點火事故反應相對于事故刺激作用過程而言,是一種非瞬時反應,嚴格意義上,其是一個涉及結(jié)構(gòu)響應、斷裂與炸藥化學反應耦合的多物理、多階段、多部位的反應過程;事故點火后的烈度增長需經(jīng)歷涉及初始燃燒反應通過縫隙或炸藥裂紋傳播、燃面擴展的低壓燃燒特征階段。對典型的密實炸藥,通常會因結(jié)構(gòu)的變形和過早解體而使事故反應烈度被限制在僅少量炸藥反應弱驅(qū)動輸出的爆燃水平,其相應的毀傷效應極其有限;非沖擊點火反應演化過程中,劇烈反應的生成注定與非鈍感炸藥的強約束狀態(tài),以及炸藥基體在點火前預破碎或反應過程驅(qū)動的炸藥基體動態(tài)破碎有關聯(lián);

    (3)炸藥表面亞音速燃燒形成的應力波因幅值和前沿寬度限制,無法直接激發(fā)炸藥內(nèi)部微介觀熱點形成點火,所以,假定密實炸藥基體內(nèi)存在分布式點火、按比例緩慢反應增長的相關模型缺乏源頭概念立足依據(jù);對典型主炸藥而言,表面燃燒引發(fā)的應力波的壓力增長因燃速限制,很難在炸藥基體中形成強沖擊波和急劇體積壓縮,引發(fā)類似粉體炸藥中的DDT。通過近期開展的系列厚壁、大長度DDT 管實驗的詳實證據(jù),糾正了早期實驗中因基本現(xiàn)象解讀錯誤,對相關PBX 形成DDT 轉(zhuǎn)化的誤判;

    (4)因事故碰撞嚴重預破壞或因非沖擊點火反應進程中高壓驅(qū)動引發(fā)的炸藥破碎,可在相對較弱的約束條件下表現(xiàn)出極為劇烈的爆炸表現(xiàn)。這種超大比表面積的反應形態(tài),從遠區(qū)超壓看非常類似爆轟,甚至可能給出高于對應炸藥完全爆轟的超壓當量,碎片速度也較高,但不能據(jù)此認定事故反應過程轉(zhuǎn)入了爆轟。在這種反應進程中,劇烈變形的裝藥結(jié)構(gòu)中,反應區(qū)原位壓力可能遠低于爆轟壓力,但反應持續(xù)時間極長,多為毫秒級特征,使得相當部分細碎炸藥來得及完成反應,形成大量氣態(tài)產(chǎn)物驅(qū)動殼體飛散和高水平超壓。對這種區(qū)別于一維厚壁構(gòu)型DDT 過程機制、被部分同行稱為延期爆轟,即XDT 的實驗現(xiàn)象,值得從真實過程演化細節(jié)上重新再審視,落實同行對其發(fā)生爆轟的判定是基于外在效應似然,還是符合對密實炸藥爆轟定義的實證數(shù)據(jù)。

    6 結(jié) 論

    通過對炸藥非沖擊點火事故反應現(xiàn)象實驗診斷結(jié)果的分析解讀,明確了非沖擊點火反應主導機制是炸藥表面層流燃燒與氣態(tài)燃燒產(chǎn)物的縫隙對流,而非炸藥基體中的反應在主導。從本質(zhì)上講,炸藥非沖擊點火事故反應屬于層流燃燒,但發(fā)生在極端復雜的結(jié)構(gòu)中,基體中反應裂紋演化、產(chǎn)物對流和結(jié)構(gòu)變形、破壞,使反應烈度轉(zhuǎn)變變得極度復雜。這一基本認識是炸藥事故反應行為解讀、實驗設計和數(shù)值模擬物理-數(shù)學模型建模、數(shù)值模擬校驗的基本出發(fā)點。

    猜你喜歡
    柱殼烈度裝藥
    火炸藥學報(2022年3期)2022-07-04 07:31:00
    孔內(nèi)爆炸填塞效應的數(shù)值模擬與分析
    高烈度區(qū)域深基坑基坑支護設計
    電磁驅(qū)動薄壁鋁合金柱殼的動態(tài)膨脹變形特征
    某發(fā)射裝藥結(jié)構(gòu)改進設計的新思路
    高烈度地震區(qū)非規(guī)則多跨長聯(lián)連續(xù)梁抗震分析
    偏心柱殼自由振動的級數(shù)變換求解方法
    1/3含口蓋復合材料柱殼后屈曲性能
    深孔不耦合裝藥爆破技術(shù)卸壓效果驗證
    變截面加筋柱殼結(jié)構(gòu)振動特性分析
    国产有黄有色有爽视频| 成人手机av| 国产成人午夜福利电影在线观看| 亚洲国产精品成人久久小说| 免费观看无遮挡的男女| 日韩制服骚丝袜av| 大陆偷拍与自拍| 青春草视频在线免费观看| 国产亚洲一区二区精品| 国产精品久久久久久av不卡| 日本猛色少妇xxxxx猛交久久| 国产精品国产三级专区第一集| 成年女人在线观看亚洲视频| 高清在线视频一区二区三区| 最近手机中文字幕大全| 女人久久www免费人成看片| 老司机亚洲免费影院| 一级毛片 在线播放| 99re6热这里在线精品视频| 精品少妇久久久久久888优播| 狠狠婷婷综合久久久久久88av| 香蕉精品网在线| 爱豆传媒免费全集在线观看| 免费看光身美女| 久热久热在线精品观看| 久久国产亚洲av麻豆专区| 免费高清在线观看视频在线观看| 少妇的逼水好多| 亚洲精品成人av观看孕妇| 一级片'在线观看视频| 满18在线观看网站| 中文字幕亚洲精品专区| 搡老乐熟女国产| av网站免费在线观看视频| 亚洲成人一二三区av| 亚洲国产最新在线播放| 免费观看a级毛片全部| 永久网站在线| 久久久欧美国产精品| 日韩视频在线欧美| 国产一级毛片在线| 国产 精品1| 王馨瑶露胸无遮挡在线观看| 伦理电影免费视频| 国产永久视频网站| 国产在线免费精品| 久久99一区二区三区| 男的添女的下面高潮视频| 男女免费视频国产| 亚洲精品自拍成人| av福利片在线| 91aial.com中文字幕在线观看| 久久热精品热| 你懂的网址亚洲精品在线观看| 成人免费观看视频高清| 欧美激情极品国产一区二区三区 | 国产免费一级a男人的天堂| 欧美 日韩 精品 国产| 欧美日韩国产mv在线观看视频| 性色avwww在线观看| av电影中文网址| 亚洲国产精品成人久久小说| 亚洲婷婷狠狠爱综合网| 一本久久精品| 免费日韩欧美在线观看| 亚洲精品国产av成人精品| 久久97久久精品| 免费观看性生交大片5| 国产成人精品婷婷| 十分钟在线观看高清视频www| 少妇人妻 视频| 亚洲国产av新网站| 我的女老师完整版在线观看| www.色视频.com| 成人黄色视频免费在线看| 26uuu在线亚洲综合色| 中文欧美无线码| 日本黄大片高清| 久久女婷五月综合色啪小说| 你懂的网址亚洲精品在线观看| 午夜福利在线观看免费完整高清在| 九九在线视频观看精品| 国产一区二区三区综合在线观看 | 国产在视频线精品| 在线观看免费高清a一片| av卡一久久| 高清欧美精品videossex| 亚洲婷婷狠狠爱综合网| 亚洲精品第二区| 国产精品国产三级国产专区5o| 欧美 亚洲 国产 日韩一| 能在线免费看毛片的网站| 亚洲综合精品二区| 午夜老司机福利剧场| 久久97久久精品| 亚洲欧美中文字幕日韩二区| 在线免费观看不下载黄p国产| 99久久精品国产国产毛片| 国产在视频线精品| 午夜激情av网站| 下体分泌物呈黄色| 亚洲欧美中文字幕日韩二区| 亚洲四区av| 青春草国产在线视频| 国产精品久久久久久精品古装| 女的被弄到高潮叫床怎么办| 亚洲精品第二区| 成人二区视频| 久久精品国产a三级三级三级| 欧美成人精品欧美一级黄| 久久久久久久大尺度免费视频| 国产精品久久久久久久久免| 内地一区二区视频在线| 特大巨黑吊av在线直播| 天堂俺去俺来也www色官网| 三级国产精品片| 日本黄色片子视频| 久久午夜综合久久蜜桃| 黄色怎么调成土黄色| 啦啦啦视频在线资源免费观看| 亚洲高清免费不卡视频| 国产精品一区二区在线不卡| 久久青草综合色| 日本猛色少妇xxxxx猛交久久| 少妇被粗大猛烈的视频| 一区二区av电影网| 国产精品秋霞免费鲁丝片| www.色视频.com| 高清av免费在线| 色视频在线一区二区三区| 黄片播放在线免费| 高清av免费在线| av国产精品久久久久影院| 久久久久网色| 涩涩av久久男人的天堂| 国产在线视频一区二区| 精品国产一区二区三区久久久樱花| 国产免费现黄频在线看| 寂寞人妻少妇视频99o| 超色免费av| 亚洲怡红院男人天堂| 看非洲黑人一级黄片| 午夜激情久久久久久久| 一级毛片电影观看| 成人漫画全彩无遮挡| 久久韩国三级中文字幕| 中文乱码字字幕精品一区二区三区| 九色亚洲精品在线播放| 91精品国产国语对白视频| 亚州av有码| 男女边摸边吃奶| .国产精品久久| 精品一品国产午夜福利视频| 国产无遮挡羞羞视频在线观看| 高清欧美精品videossex| 嫩草影院入口| 久久ye,这里只有精品| 国产成人91sexporn| 男的添女的下面高潮视频| 九九在线视频观看精品| 97超碰精品成人国产| 在线观看一区二区三区激情| 国产精品一二三区在线看| 毛片一级片免费看久久久久| 中文字幕最新亚洲高清| 伦理电影免费视频| a级毛片免费高清观看在线播放| 18禁在线无遮挡免费观看视频| 黄色视频在线播放观看不卡| 久久精品久久久久久久性| 国产深夜福利视频在线观看| 久久久a久久爽久久v久久| 精品一区二区三卡| 91精品一卡2卡3卡4卡| 少妇精品久久久久久久| 精品视频人人做人人爽| 国产伦理片在线播放av一区| 人人妻人人澡人人爽人人夜夜| 中文字幕人妻丝袜制服| 亚洲精品,欧美精品| 免费少妇av软件| 老女人水多毛片| 内地一区二区视频在线| 久久精品国产鲁丝片午夜精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产欧美亚洲国产| 人妻夜夜爽99麻豆av| 高清黄色对白视频在线免费看| 18在线观看网站| 日本午夜av视频| 亚洲伊人久久精品综合| 在线观看三级黄色| 91久久精品国产一区二区三区| 亚洲av电影在线观看一区二区三区| xxxhd国产人妻xxx| 国产免费一区二区三区四区乱码| 高清欧美精品videossex| 日本黄色片子视频| 极品少妇高潮喷水抽搐| av在线app专区| 秋霞伦理黄片| 又黄又爽又刺激的免费视频.| 日韩视频在线欧美| 午夜视频国产福利| 黄色配什么色好看| 精品少妇黑人巨大在线播放| 男女边摸边吃奶| 久久久久视频综合| 搡老乐熟女国产| 嫩草影院入口| 亚洲四区av| 欧美 日韩 精品 国产| 亚洲无线观看免费| 大香蕉久久网| h视频一区二区三区| 精品久久久精品久久久| 亚洲高清免费不卡视频| 春色校园在线视频观看| 午夜激情福利司机影院| 久久国产精品男人的天堂亚洲 | 中文字幕制服av| 午夜福利视频在线观看免费| 久久精品国产鲁丝片午夜精品| 亚洲丝袜综合中文字幕| 国产精品不卡视频一区二区| 久久久精品免费免费高清| 99久久人妻综合| 美女中出高潮动态图| 国产毛片在线视频| 国产免费现黄频在线看| 纵有疾风起免费观看全集完整版| 日韩人妻高清精品专区| 成人手机av| av免费在线看不卡| 国产视频首页在线观看| 中文字幕最新亚洲高清| 午夜免费鲁丝| 我要看黄色一级片免费的| 在线免费观看不下载黄p国产| 内地一区二区视频在线| 久热这里只有精品99| 国产免费福利视频在线观看| 最近最新中文字幕免费大全7| 久久综合国产亚洲精品| 日本爱情动作片www.在线观看| 午夜激情福利司机影院| 久久婷婷青草| 免费高清在线观看日韩| 久久99精品国语久久久| 久久久国产一区二区| 少妇猛男粗大的猛烈进出视频| 国内精品宾馆在线| 久久久久久久大尺度免费视频| 伊人久久国产一区二区| 亚洲在久久综合| 亚洲精品aⅴ在线观看| 日韩不卡一区二区三区视频在线| 精品一区二区免费观看| 丰满迷人的少妇在线观看| 99久久精品国产国产毛片| 在线观看www视频免费| av在线老鸭窝| 伊人久久精品亚洲午夜| 在线看a的网站| 少妇人妻 视频| 大片电影免费在线观看免费| 综合色丁香网| 亚洲国产精品成人久久小说| 国产爽快片一区二区三区| 国产精品国产av在线观看| 最近中文字幕2019免费版| 黄片无遮挡物在线观看| 18禁裸乳无遮挡动漫免费视频| 一级二级三级毛片免费看| 国内精品宾馆在线| 纯流量卡能插随身wifi吗| h视频一区二区三区| 国产av精品麻豆| 69精品国产乱码久久久| 啦啦啦在线观看免费高清www| av国产精品久久久久影院| 亚洲国产色片| 亚洲少妇的诱惑av| 国产不卡av网站在线观看| 久久久久久久久久久免费av| 成人亚洲欧美一区二区av| 熟女电影av网| 3wmmmm亚洲av在线观看| 欧美另类一区| av电影中文网址| 国产一级毛片在线| 久久av网站| 亚洲国产精品一区二区三区在线| 五月伊人婷婷丁香| 中国国产av一级| 亚洲欧美一区二区三区国产| 欧美精品高潮呻吟av久久| 成人亚洲欧美一区二区av| 少妇被粗大猛烈的视频| 少妇的逼好多水| 黄色毛片三级朝国网站| 蜜桃久久精品国产亚洲av| 高清视频免费观看一区二区| 亚洲国产精品一区二区三区在线| 视频区图区小说| 狂野欧美激情性xxxx在线观看| 一本一本综合久久| 中文字幕免费在线视频6| 国产精品人妻久久久久久| 国产精品人妻久久久影院| 日韩强制内射视频| 亚洲av免费高清在线观看| 国产乱人偷精品视频| 看免费成人av毛片| 91aial.com中文字幕在线观看| 大片电影免费在线观看免费| 精品国产国语对白av| 亚洲人成网站在线播| av有码第一页| 免费观看a级毛片全部| 欧美日韩视频高清一区二区三区二| 久久精品国产亚洲网站| 亚洲av成人精品一区久久| 岛国毛片在线播放| 国产午夜精品久久久久久一区二区三区| 国产精品免费大片| 99久久精品国产国产毛片| 亚洲欧美色中文字幕在线| 男的添女的下面高潮视频| 夫妻午夜视频| 欧美亚洲日本最大视频资源| 中文字幕人妻丝袜制服| 制服人妻中文乱码| 一边摸一边做爽爽视频免费| 国产成人一区二区在线| 久久99蜜桃精品久久| 久久99热这里只频精品6学生| 欧美日韩亚洲高清精品| 国产欧美亚洲国产| 又大又黄又爽视频免费| 99久久综合免费| 亚洲精品一区蜜桃| 免费高清在线观看日韩| 人人妻人人澡人人看| av国产久精品久网站免费入址| 伊人久久国产一区二区| 高清午夜精品一区二区三区| 亚洲欧洲日产国产| 欧美3d第一页| 欧美日本中文国产一区发布| 国产黄色免费在线视频| 中国三级夫妇交换| a级毛片黄视频| 一级毛片aaaaaa免费看小| av有码第一页| 午夜福利网站1000一区二区三区| 精品久久久久久久久亚洲| 搡女人真爽免费视频火全软件| 亚洲成人手机| 亚洲av在线观看美女高潮| 成人亚洲欧美一区二区av| 女人精品久久久久毛片| 在线免费观看不下载黄p国产| 欧美人与善性xxx| 久久久久精品久久久久真实原创| 欧美精品亚洲一区二区| 国产精品99久久久久久久久| 国产 一区精品| 日日摸夜夜添夜夜爱| 中文字幕免费在线视频6| av有码第一页| 精品久久久久久电影网| 亚洲成人手机| 91成人精品电影| 在线天堂最新版资源| 日本vs欧美在线观看视频| 免费观看的影片在线观看| 又黄又爽又刺激的免费视频.| 久久久久久久久久人人人人人人| 久久久久久久久大av| 午夜福利影视在线免费观看| 丁香六月天网| 人妻夜夜爽99麻豆av| 国产男女内射视频| 亚洲五月色婷婷综合| 欧美人与善性xxx| 日韩不卡一区二区三区视频在线| 好男人视频免费观看在线| 精品人妻一区二区三区麻豆| 你懂的网址亚洲精品在线观看| 国产成人精品婷婷| 久久精品夜色国产| 熟女av电影| 国产熟女午夜一区二区三区 | 欧美一级a爱片免费观看看| 国产免费福利视频在线观看| 亚洲精品色激情综合| 婷婷色综合www| 高清在线视频一区二区三区| 久久人妻熟女aⅴ| 国产高清国产精品国产三级| 成人二区视频| 女性生殖器流出的白浆| 国产成人免费观看mmmm| 亚洲四区av| 日韩欧美一区视频在线观看| 亚洲美女黄色视频免费看| 欧美性感艳星| 中文字幕制服av| 中文字幕人妻丝袜制服| 久久久国产一区二区| 亚洲欧美日韩另类电影网站| 欧美+日韩+精品| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| tube8黄色片| 天天躁夜夜躁狠狠久久av| 中国国产av一级| 婷婷色av中文字幕| 看十八女毛片水多多多| 亚洲国产精品专区欧美| 免费观看无遮挡的男女| 欧美成人午夜免费资源| 日本免费在线观看一区| 亚洲av福利一区| 啦啦啦啦在线视频资源| 欧美最新免费一区二区三区| 日韩三级伦理在线观看| 亚洲精品成人av观看孕妇| 久久精品人人爽人人爽视色| 黑人高潮一二区| 一边摸一边做爽爽视频免费| 伦理电影大哥的女人| 久久久久视频综合| 99热6这里只有精品| 亚洲第一区二区三区不卡| 在线观看国产h片| 十八禁高潮呻吟视频| 成人亚洲精品一区在线观看| 成人国产av品久久久| 午夜久久久在线观看| 国产精品蜜桃在线观看| 免费观看在线日韩| 黑人猛操日本美女一级片| 亚洲国产欧美日韩在线播放| 看十八女毛片水多多多| 国产欧美另类精品又又久久亚洲欧美| 美女国产高潮福利片在线看| 91精品一卡2卡3卡4卡| 99热国产这里只有精品6| 国产精品久久久久成人av| 国产精品不卡视频一区二区| 日日撸夜夜添| 秋霞在线观看毛片| 大陆偷拍与自拍| 亚洲欧洲日产国产| 亚洲成人一二三区av| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产精品专区欧美| 亚洲国产色片| 久久久久国产网址| 在线天堂最新版资源| 高清av免费在线| 亚洲人成网站在线播| 欧美精品人与动牲交sv欧美| 一级二级三级毛片免费看| 人妻 亚洲 视频| 亚州av有码| 飞空精品影院首页| 22中文网久久字幕| 中文字幕av电影在线播放| 国产色婷婷99| av免费观看日本| 国产精品久久久久久精品电影小说| 亚洲精品美女久久av网站| 亚洲一级一片aⅴ在线观看| 香蕉精品网在线| 全区人妻精品视频| 久久久久人妻精品一区果冻| 久久精品国产鲁丝片午夜精品| 婷婷色麻豆天堂久久| 性色avwww在线观看| 亚洲国产精品999| 国产欧美另类精品又又久久亚洲欧美| 三级国产精品片| 肉色欧美久久久久久久蜜桃| 9色porny在线观看| 亚洲国产精品国产精品| 国产黄色免费在线视频| 亚洲国产最新在线播放| 久久亚洲国产成人精品v| 亚洲国产精品专区欧美| 老熟女久久久| 亚洲色图 男人天堂 中文字幕 | 日韩伦理黄色片| 亚洲欧美一区二区三区国产| 免费高清在线观看日韩| 黑人欧美特级aaaaaa片| 精品久久久噜噜| 精品国产一区二区久久| 少妇 在线观看| 国产av码专区亚洲av| 免费观看无遮挡的男女| 国产精品99久久久久久久久| 国产成人精品婷婷| 97在线人人人人妻| av在线观看视频网站免费| 国内精品宾馆在线| 亚洲中文av在线| 久久ye,这里只有精品| 免费少妇av软件| www.色视频.com| 如何舔出高潮| 蜜臀久久99精品久久宅男| 亚洲国产最新在线播放| 蜜桃在线观看..| 亚洲精品乱久久久久久| 中文字幕免费在线视频6| 91国产中文字幕| 中文乱码字字幕精品一区二区三区| 亚洲国产日韩一区二区| 丰满乱子伦码专区| 新久久久久国产一级毛片| 一区在线观看完整版| 秋霞在线观看毛片| 久久人妻熟女aⅴ| 国产精品久久久久久av不卡| 久久亚洲国产成人精品v| 国产成人免费观看mmmm| 国产一级毛片在线| 亚洲激情五月婷婷啪啪| 热99国产精品久久久久久7| 国产白丝娇喘喷水9色精品| 久久ye,这里只有精品| 成人综合一区亚洲| 美女国产视频在线观看| 亚洲成人一二三区av| 日韩av不卡免费在线播放| 99九九在线精品视频| 高清视频免费观看一区二区| 丰满乱子伦码专区| 极品少妇高潮喷水抽搐| 五月伊人婷婷丁香| 久久99热6这里只有精品| 国产永久视频网站| 美女cb高潮喷水在线观看| 久久精品人人爽人人爽视色| 久久精品久久精品一区二区三区| 人妻 亚洲 视频| 又粗又硬又长又爽又黄的视频| 纯流量卡能插随身wifi吗| 街头女战士在线观看网站| 国产免费一区二区三区四区乱码| 久久狼人影院| 日韩三级伦理在线观看| 亚洲不卡免费看| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 大话2 男鬼变身卡| 大香蕉97超碰在线| 国产精品久久久久久久电影| 免费观看无遮挡的男女| 岛国毛片在线播放| 免费大片18禁| 国产不卡av网站在线观看| 精品久久久久久久久av| 国产高清三级在线| 99热网站在线观看| 女的被弄到高潮叫床怎么办| 如日韩欧美国产精品一区二区三区 | 日日撸夜夜添| 久久精品国产a三级三级三级| 国产熟女午夜一区二区三区 | av.在线天堂| av电影中文网址| 国产av码专区亚洲av| 男男h啪啪无遮挡| 老司机影院毛片| 久久 成人 亚洲| 精品酒店卫生间| 精品国产国语对白av| 三级国产精品欧美在线观看| 99热这里只有是精品在线观看| 久久精品久久久久久噜噜老黄| 日日啪夜夜爽| 少妇的逼好多水| 亚洲国产精品999| 一二三四中文在线观看免费高清| 免费av不卡在线播放| 亚洲精品国产av成人精品| 精品人妻偷拍中文字幕| freevideosex欧美| 欧美日韩av久久| 日本av免费视频播放| 亚洲av在线观看美女高潮| 日韩中文字幕视频在线看片| 一区二区日韩欧美中文字幕 | 菩萨蛮人人尽说江南好唐韦庄| 免费少妇av软件| 青春草亚洲视频在线观看| 岛国毛片在线播放| 日韩欧美精品免费久久| 亚洲av欧美aⅴ国产| 婷婷成人精品国产| 色婷婷av一区二区三区视频| 免费人成在线观看视频色| 黄色怎么调成土黄色| 欧美日韩综合久久久久久| 日韩亚洲欧美综合| 2022亚洲国产成人精品| 特大巨黑吊av在线直播| 成人国产av品久久久| 国产成人aa在线观看| 欧美日韩精品成人综合77777| 久久久久久久亚洲中文字幕|