• <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| 午夜爱爱视频在线播放| 性色avwww在线观看| 能在线免费看毛片的网站| 最近手机中文字幕大全| 亚洲va在线va天堂va国产| 亚洲高清免费不卡视频| 国产精品久久久久久av不卡| 国产高清有码在线观看视频| 69人妻影院| 偷拍熟女少妇极品色| 乱人视频在线观看| 亚洲av一区综合| 黄色视频,在线免费观看| 国产一区二区激情短视频| 久久久久久久久久黄片| 亚洲av免费在线观看| 亚洲18禁久久av| 精品久久久久久久久av| 国产一区二区在线av高清观看| 综合色av麻豆| 欧美性猛交╳xxx乱大交人| 偷拍熟女少妇极品色| 欧美激情在线99| 性插视频无遮挡在线免费观看| 国产精品久久久久久av不卡| 一级毛片我不卡| 欧美成人一区二区免费高清观看| 男女视频在线观看网站免费| 欧美在线一区亚洲| 亚洲精品久久国产高清桃花| 春色校园在线视频观看| 国产片特级美女逼逼视频| 久久久久网色| 麻豆精品久久久久久蜜桃| 麻豆乱淫一区二区| 内地一区二区视频在线| 观看免费一级毛片| 别揉我奶头 嗯啊视频| 色尼玛亚洲综合影院| 日韩欧美国产在线观看| 久久精品国产鲁丝片午夜精品| 狂野欧美白嫩少妇大欣赏| 欧美成人精品欧美一级黄| 国内久久婷婷六月综合欲色啪| 日韩欧美在线乱码| 成人欧美大片| 久久久久久久久中文| 伊人久久精品亚洲午夜| 在线天堂最新版资源| 国产精品综合久久久久久久免费| 日产精品乱码卡一卡2卡三| 久久精品国产自在天天线| 色综合亚洲欧美另类图片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 级片在线观看| 99热这里只有是精品50| 久久久久久久久中文| 两性午夜刺激爽爽歪歪视频在线观看| 搞女人的毛片| av天堂中文字幕网| 精品熟女少妇av免费看| 亚洲电影在线观看av| 国内精品美女久久久久久| 色哟哟哟哟哟哟| 亚洲av二区三区四区| 国产精品久久久久久久电影| 国产精品一区二区在线观看99 | 日韩中字成人| avwww免费| 精品一区二区免费观看| 长腿黑丝高跟| 亚洲国产欧美人成| 一夜夜www| 国产成人影院久久av| 久久精品夜色国产| 特大巨黑吊av在线直播| 九九爱精品视频在线观看| 亚洲人成网站高清观看| 日日摸夜夜添夜夜添av毛片| 少妇人妻精品综合一区二区 | 午夜视频国产福利| 人妻制服诱惑在线中文字幕| 人体艺术视频欧美日本| 国产免费一级a男人的天堂| 日本五十路高清| 日韩精品有码人妻一区| 五月玫瑰六月丁香| 丰满乱子伦码专区| 两个人视频免费观看高清| 欧美色视频一区免费| 又爽又黄无遮挡网站| 国产久久久一区二区三区| 国产日韩欧美在线精品| 亚洲性久久影院| 免费观看在线日韩| 我的老师免费观看完整版| 蜜臀久久99精品久久宅男| 色综合亚洲欧美另类图片| 丰满人妻一区二区三区视频av| 变态另类丝袜制服| 青青草视频在线视频观看| 禁无遮挡网站| 亚洲成a人片在线一区二区| 高清在线视频一区二区三区 | 亚洲人与动物交配视频| 一级av片app| 日韩大尺度精品在线看网址| 中文精品一卡2卡3卡4更新| 亚洲国产欧洲综合997久久,| 亚洲无线在线观看| 一个人观看的视频www高清免费观看| 精品欧美国产一区二区三| 欧美成人精品欧美一级黄| 日韩视频在线欧美| 久久久久国产网址| 国产极品精品免费视频能看的| 日本与韩国留学比较| 欧美精品国产亚洲| 日韩av在线大香蕉| 国产精品一区www在线观看| 中文字幕制服av| 99热精品在线国产| av在线亚洲专区| 在线观看免费视频日本深夜| 国产极品天堂在线| av免费在线看不卡| 乱人视频在线观看| 美女cb高潮喷水在线观看| 蜜桃亚洲精品一区二区三区| 亚洲三级黄色毛片| 久久久久国产网址| 91午夜精品亚洲一区二区三区| 色综合亚洲欧美另类图片| 日本免费a在线| 亚洲最大成人手机在线| 亚洲精品456在线播放app| 成人综合一区亚洲| 婷婷亚洲欧美| 国产高潮美女av| www日本黄色视频网| 中文精品一卡2卡3卡4更新| 亚洲五月天丁香| 乱系列少妇在线播放| 亚洲av电影不卡..在线观看| 久久精品国产亚洲网站| 久久精品国产鲁丝片午夜精品| 中文资源天堂在线| 人妻夜夜爽99麻豆av| 日韩欧美一区二区三区在线观看| 久久久久久久午夜电影| 18禁在线无遮挡免费观看视频| 99热6这里只有精品| 国产色爽女视频免费观看| 天天一区二区日本电影三级| 久久久久久九九精品二区国产| 国产极品天堂在线| 中文精品一卡2卡3卡4更新| 日韩av在线大香蕉| 麻豆av噜噜一区二区三区| 国产成人精品一,二区 | 中文精品一卡2卡3卡4更新| 色播亚洲综合网| a级毛色黄片| 免费观看精品视频网站| av福利片在线观看| 久久99精品国语久久久| 99热只有精品国产| 综合色丁香网| 在线免费十八禁| 亚洲国产精品国产精品| 国产精品人妻久久久久久| 日本黄色视频三级网站网址| 色噜噜av男人的天堂激情| 成年女人永久免费观看视频| 毛片一级片免费看久久久久| 青春草国产在线视频 | 国产午夜福利久久久久久| 国产伦精品一区二区三区视频9| 成人漫画全彩无遮挡| 在现免费观看毛片| 色尼玛亚洲综合影院| 婷婷亚洲欧美| 亚洲真实伦在线观看| 哪里可以看免费的av片| 在线观看午夜福利视频| 色视频www国产| 国产色婷婷99| 日本成人三级电影网站| 听说在线观看完整版免费高清| 一夜夜www| 久久久国产成人精品二区| 真实男女啪啪啪动态图| 伦理电影大哥的女人| 蜜臀久久99精品久久宅男| 亚洲18禁久久av| 久久久久九九精品影院| 特大巨黑吊av在线直播| 在线播放无遮挡| 边亲边吃奶的免费视频| 黑人高潮一二区| av卡一久久| 亚洲av电影不卡..在线观看| 国产精品不卡视频一区二区| 身体一侧抽搐| 在线免费观看的www视频| 免费搜索国产男女视频| 性欧美人与动物交配| 少妇的逼水好多| 1000部很黄的大片| 日日撸夜夜添| 一区二区三区四区激情视频 | 国产成人精品久久久久久| 亚洲欧洲国产日韩| 我要看日韩黄色一级片| 一进一出抽搐动态| 啦啦啦韩国在线观看视频| 日本黄大片高清| 国产男人的电影天堂91| 亚洲四区av| 日本一本二区三区精品| 变态另类成人亚洲欧美熟女| 亚洲欧洲国产日韩| 99国产精品一区二区蜜桃av| 国产一级毛片在线| 亚洲欧美成人精品一区二区| 欧美精品一区二区大全| 国产一区二区激情短视频| 久久久色成人| 国产成人精品一,二区 | 免费观看的影片在线观看| av在线老鸭窝| 夫妻性生交免费视频一级片| 久久精品夜夜夜夜夜久久蜜豆| 91久久精品国产一区二区成人| 国产精品蜜桃在线观看 | 搡女人真爽免费视频火全软件| 岛国在线免费视频观看| 干丝袜人妻中文字幕| 亚洲无线观看免费| 蜜臀久久99精品久久宅男| 91久久精品国产一区二区三区| 精品无人区乱码1区二区| 免费一级毛片在线播放高清视频| 两个人的视频大全免费| 亚洲av成人av| 国产精品一区www在线观看| 国产av麻豆久久久久久久| av福利片在线观看| 久久久久久九九精品二区国产| 深夜精品福利| 日韩强制内射视频| 日韩 亚洲 欧美在线| 国产av麻豆久久久久久久| 一进一出抽搐gif免费好疼| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一区二区激情短视频| 日韩 亚洲 欧美在线| 免费观看精品视频网站| 一区二区三区高清视频在线| 国产成人福利小说| 亚洲国产色片| 丰满乱子伦码专区| 高清毛片免费观看视频网站| 国产日韩欧美在线精品| 永久网站在线| 床上黄色一级片| 久久人人爽人人片av| 国产伦理片在线播放av一区 | 欧美一区二区国产精品久久精品| 99在线视频只有这里精品首页| 国产高清三级在线| 你懂的网址亚洲精品在线观看 | 国产精品免费一区二区三区在线| 高清毛片免费观看视频网站| 成人综合一区亚洲| 国产乱人视频| 舔av片在线| 国产精品久久久久久久久免| 免费无遮挡裸体视频| 日日撸夜夜添| 国产精品嫩草影院av在线观看| 99精品在免费线老司机午夜| 亚洲激情五月婷婷啪啪| 丰满乱子伦码专区| 美女高潮的动态| 精品少妇黑人巨大在线播放 | 18禁黄网站禁片免费观看直播| 69人妻影院| 别揉我奶头 嗯啊视频| 亚洲欧美精品自产自拍| 丰满乱子伦码专区| 成熟少妇高潮喷水视频| 三级男女做爰猛烈吃奶摸视频| 爱豆传媒免费全集在线观看| 禁无遮挡网站| 久久久a久久爽久久v久久| av在线蜜桃| 亚洲国产欧美人成| 色视频www国产| 少妇高潮的动态图| 久久精品国产亚洲av香蕉五月| 国产亚洲av嫩草精品影院| 精品久久久噜噜| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 如何舔出高潮| 免费av不卡在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女高潮的动态| 黄色一级大片看看| 免费av毛片视频| 好男人在线观看高清免费视频| 赤兔流量卡办理| 免费看a级黄色片| 亚洲av熟女| 久久久久久久亚洲中文字幕| 两个人的视频大全免费| 丰满乱子伦码专区| 极品教师在线视频| .国产精品久久| 爱豆传媒免费全集在线观看| 久久久色成人| 亚洲av中文字字幕乱码综合| 久久人人精品亚洲av| 嫩草影院入口| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 国产黄色小视频在线观看| 别揉我奶头 嗯啊视频| 国产午夜福利久久久久久| a级毛片免费高清观看在线播放| 干丝袜人妻中文字幕| 女同久久另类99精品国产91| 一边亲一边摸免费视频| 在线免费观看不下载黄p国产| 国产精品国产高清国产av| 国产精品一及| 亚洲欧美中文字幕日韩二区| 国产色婷婷99| 91精品一卡2卡3卡4卡| 中国美白少妇内射xxxbb| 天堂网av新在线| 国产色爽女视频免费观看| 久久久久久久久久久丰满| 久99久视频精品免费| 深爱激情五月婷婷| 在现免费观看毛片| 国产探花在线观看一区二区| 欧美又色又爽又黄视频| 亚洲av熟女| 最近最新中文字幕大全电影3| 国产三级中文精品| av视频在线观看入口| 国产色婷婷99| 黄片无遮挡物在线观看| 久久午夜福利片| 久久久久久久久中文| 精品国产三级普通话版| 亚洲欧美中文字幕日韩二区| 亚洲性久久影院| 丰满的人妻完整版| 最近2019中文字幕mv第一页| 国产精品.久久久| 国产亚洲91精品色在线| 欧美丝袜亚洲另类| 人妻制服诱惑在线中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| 在线免费观看的www视频| 国产69精品久久久久777片| 人体艺术视频欧美日本| a级毛片免费高清观看在线播放| 国产v大片淫在线免费观看| 天堂中文最新版在线下载 | 精品一区二区免费观看| 熟女人妻精品中文字幕| 国产毛片a区久久久久| 欧美性猛交黑人性爽| 国产v大片淫在线免费观看| 国产一级毛片七仙女欲春2| 一夜夜www| 亚洲三级黄色毛片| 色噜噜av男人的天堂激情| 麻豆国产av国片精品| 搞女人的毛片| 中文精品一卡2卡3卡4更新| 国产91av在线免费观看| 亚洲精品亚洲一区二区| 久久午夜福利片| 国产伦精品一区二区三区四那| 亚洲美女视频黄频| 国产精品一区www在线观看| 国产黄片美女视频| 日韩强制内射视频| 午夜久久久久精精品| 国产av在哪里看| 男人的好看免费观看在线视频| 18+在线观看网站| 99热只有精品国产| 国产一级毛片在线| 边亲边吃奶的免费视频| 久久6这里有精品| 99久久精品一区二区三区| 久久精品久久久久久噜噜老黄 | 老司机影院成人| 日韩,欧美,国产一区二区三区 | 久久精品综合一区二区三区| 国产精品久久电影中文字幕| 波多野结衣高清无吗| 免费观看的影片在线观看| 免费无遮挡裸体视频| 国产精品爽爽va在线观看网站| 国产亚洲av片在线观看秒播厂 | 禁无遮挡网站| 国产成人freesex在线| 国产男人的电影天堂91| 中文字幕制服av| 欧美+亚洲+日韩+国产| 免费观看的影片在线观看| 亚洲在久久综合| 国产国拍精品亚洲av在线观看| 特大巨黑吊av在线直播| 国产爱豆传媒在线观看| 亚洲婷婷狠狠爱综合网| 国产伦在线观看视频一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 97热精品久久久久久| 欧美激情国产日韩精品一区| 国产av不卡久久| 99热这里只有精品一区| 久久久午夜欧美精品| 给我免费播放毛片高清在线观看| 欧美人与善性xxx| 男女边吃奶边做爰视频| 亚洲国产精品国产精品| 国产 一区精品| 国产片特级美女逼逼视频| 成年版毛片免费区| 嫩草影院精品99| 成人综合一区亚洲| 国产伦精品一区二区三区视频9| a级毛色黄片| 国产成人精品久久久久久| 天美传媒精品一区二区| 日韩精品青青久久久久久| 日韩av在线大香蕉| 日日摸夜夜添夜夜添av毛片| 亚洲av成人av| 欧美一级a爱片免费观看看| 99久久精品热视频| 嘟嘟电影网在线观看| 2022亚洲国产成人精品| 国产片特级美女逼逼视频| 国产亚洲精品久久久久久毛片| 日本av手机在线免费观看| 看黄色毛片网站| 久久精品影院6| 国产综合懂色| 精品国产三级普通话版| 97热精品久久久久久| 午夜福利成人在线免费观看| 欧美一级a爱片免费观看看| 欧美bdsm另类| 99热这里只有精品一区| 99久久精品热视频| 亚洲精品国产av成人精品| 伦理电影大哥的女人| 美女xxoo啪啪120秒动态图| 国产午夜精品论理片| 欧美精品国产亚洲| 男女那种视频在线观看| 精品一区二区免费观看| 性欧美人与动物交配| 麻豆成人av视频| eeuss影院久久| 国产一区二区在线观看日韩| 免费黄网站久久成人精品| 高清日韩中文字幕在线| 国产在线精品亚洲第一网站| 亚洲av免费在线观看| 日韩视频在线欧美| 床上黄色一级片| 中文字幕久久专区| 长腿黑丝高跟| 欧美高清成人免费视频www| 五月玫瑰六月丁香| 成人一区二区视频在线观看| 男女边吃奶边做爰视频| 午夜视频国产福利| 国产黄a三级三级三级人| 亚洲丝袜综合中文字幕| 亚洲熟妇中文字幕五十中出| 国产色婷婷99| 亚洲欧美成人综合另类久久久 | 国产视频内射| 亚洲真实伦在线观看| 麻豆精品久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 国产麻豆成人av免费视频| 我的老师免费观看完整版| 亚洲四区av| 日本黄色视频三级网站网址| 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av| videossex国产| 婷婷精品国产亚洲av| 亚洲精品影视一区二区三区av| 色综合色国产| 美女内射精品一级片tv| 国语自产精品视频在线第100页| 综合色丁香网| 色综合色国产| 狠狠狠狠99中文字幕| 国产毛片a区久久久久| 亚洲高清免费不卡视频| 成人特级黄色片久久久久久久| 春色校园在线视频观看| 永久网站在线| 能在线免费观看的黄片| 蜜臀久久99精品久久宅男| 啦啦啦啦在线视频资源| 身体一侧抽搐| 国产 一区精品| 免费在线观看成人毛片| 人人妻人人澡欧美一区二区| 韩国av在线不卡| av在线观看视频网站免费| 久久人人爽人人片av| 亚洲精品久久久久久婷婷小说 | 久久精品国产亚洲av香蕉五月| 少妇的逼好多水| 国内精品久久久久精免费| 日韩欧美三级三区| 国产精品麻豆人妻色哟哟久久 | 最近最新中文字幕大全电影3| 亚洲国产欧美在线一区| 精品人妻一区二区三区麻豆| 日韩三级伦理在线观看| 免费人成在线观看视频色| av在线天堂中文字幕| 欧美激情久久久久久爽电影| 免费搜索国产男女视频| 网址你懂的国产日韩在线| 六月丁香七月| 久久久久免费精品人妻一区二区| av福利片在线观看| 国产精品女同一区二区软件| 国产av麻豆久久久久久久| 国产黄a三级三级三级人| 成人欧美大片| 亚洲精品粉嫩美女一区| 国产精品野战在线观看| 中文精品一卡2卡3卡4更新| 又爽又黄a免费视频| 国产成人91sexporn| 又爽又黄a免费视频| 国产精品一及| 国产色爽女视频免费观看| 国产亚洲91精品色在线| 国产色婷婷99| av免费在线看不卡| 直男gayav资源| 日韩制服骚丝袜av| 亚洲国产精品国产精品| 久久久久久久久久黄片| 成人综合一区亚洲| 淫秽高清视频在线观看| 久久久久免费精品人妻一区二区| 久久综合国产亚洲精品| 亚洲国产精品合色在线| 69人妻影院| 久久久久久久久久久免费av| 插阴视频在线观看视频| 欧美三级亚洲精品| 1024手机看黄色片| 男的添女的下面高潮视频| 乱系列少妇在线播放| 人妻夜夜爽99麻豆av| 亚洲人与动物交配视频| 99国产极品粉嫩在线观看| 色视频www国产| 热99在线观看视频| 国产精品久久电影中文字幕| 久久精品国产亚洲av香蕉五月| 又黄又爽又刺激的免费视频.| 免费不卡的大黄色大毛片视频在线观看 | 欧美日韩精品成人综合77777| 性欧美人与动物交配| 精华霜和精华液先用哪个| 人妻久久中文字幕网| 搡女人真爽免费视频火全软件| 日本三级黄在线观看| 久久久久久久久久黄片| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 99热精品在线国产| 尾随美女入室| 亚洲av.av天堂| 精品久久久久久久久久免费视频| 精品久久国产蜜桃| 久久亚洲国产成人精品v| 97热精品久久久久久| 精品人妻偷拍中文字幕| 国产单亲对白刺激| 日日撸夜夜添| 欧美潮喷喷水| 亚洲精品粉嫩美女一区| 99热这里只有精品一区| 小说图片视频综合网站| 久久九九热精品免费| 亚洲av中文av极速乱| 性欧美人与动物交配|