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

    沖擊波和破片聯(lián)合作用下I型夾層板毀傷仿真

    2015-02-21 02:17:03段新峰程遠(yuǎn)勝張攀劉均李勇
    中國(guó)艦船研究 2015年6期
    關(guān)鍵詞:變形

    段新峰,程遠(yuǎn)勝,張攀,劉均,李勇

    華中科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074

    沖擊波和破片聯(lián)合作用下I型夾層板毀傷仿真

    段新峰,程遠(yuǎn)勝,張攀,劉均,李勇

    華中科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074

    基于非線性有限元軟件LS-DYNA,通過在TNT炸藥底部布置預(yù)制破片模擬戰(zhàn)斗部爆炸產(chǎn)生的沖擊波與破片聯(lián)合作用載荷,計(jì)算3種TNT炸藥當(dāng)量下I型夾層板的毀傷響應(yīng),分析沖擊波單獨(dú)作用及沖擊波與破片聯(lián)合作用下I型夾層板失效模式的差異,研究夾層板芯層配置以及上、下面板厚度配置對(duì)其失效模式的影響,并與等效實(shí)體板的抗毀傷性能進(jìn)行對(duì)比。同時(shí),從吸能的角度分析不同載荷工況下I型夾層板的吸能特性。數(shù)值仿真結(jié)果表明:在沖擊波與破片聯(lián)合作用下,結(jié)構(gòu)的毀傷程度遠(yuǎn)大于沖擊波單獨(dú)作用時(shí);當(dāng)載荷強(qiáng)度較小時(shí),I型夾層板的抗毀傷性能優(yōu)于等效實(shí)體板;載荷強(qiáng)度、載荷類型(沖擊波單獨(dú)作用或沖擊波與破片聯(lián)合作用)及上、下面板厚度配置對(duì)I型夾層板的失效模式有較大影響;從吸能特性來看,在沖擊波單獨(dú)作用下,上面板和芯層是主要的吸能構(gòu)件,而在沖擊波與破片聯(lián)合作用下,上面板和下面板是主要的吸能構(gòu)件。

    夾層結(jié)構(gòu);沖擊波;破片;LS-DYNA;失效模式

    0 引 言

    隨著各種精確制導(dǎo)武器的快速發(fā)展,艦船結(jié)構(gòu)遭受接觸爆炸及近距非接觸爆炸的概率顯著提高,生命力受到嚴(yán)重威脅。在這種近距條件下,常規(guī)武器戰(zhàn)斗部爆炸產(chǎn)生的沖擊波和大量高速破片幾乎同時(shí)作用于船體結(jié)構(gòu),因此在評(píng)估防護(hù)結(jié)構(gòu)的效能時(shí)必須考慮這兩種破壞效應(yīng)的耦合作用[1]。可見,研究艦船結(jié)構(gòu)在沖擊波與破片聯(lián)合作用下的變形及破壞模式具有重要意義。

    國(guó)內(nèi)外學(xué)者針對(duì)沖擊波與破片的聯(lián)合毀傷效應(yīng)開展了一系列試驗(yàn)及數(shù)值研究。在試驗(yàn)方面,侯海量等[2]采用帶殼裝藥和裸裝藥分別進(jìn)行艙內(nèi)爆炸模型試驗(yàn),分析了爆炸沖擊波和高速破片聯(lián)合作用對(duì)艙室結(jié)構(gòu)破壞模式的影響,結(jié)果表明,對(duì)于彈孔密集作用區(qū),彈孔間的邊界會(huì)由于應(yīng)力集中而撕開形成相互連通的大破口;張成亮等[3]采用粘貼了預(yù)制破片和未粘貼預(yù)制破片的TNT裝藥分別進(jìn)行鋼—玻璃鋼—鋼夾層結(jié)構(gòu)空中近場(chǎng)爆炸模型試驗(yàn),分析了結(jié)構(gòu)變形破壞模式及沖擊波與高速破片的聯(lián)合毀傷機(jī)制。

    目前,對(duì)于沖擊波和破片聯(lián)合作用問題的數(shù)值研究,主要采取3種方式:

    1)通過預(yù)制孔或預(yù)制溝槽模擬破片的毀傷作用,然后采用流固耦合方法或者直接施加壓力時(shí)程曲線計(jì)算結(jié)構(gòu)在后續(xù)沖擊波作用下的響應(yīng)[4-6]。這種方法僅適用于破片先于沖擊波到達(dá)目標(biāo)結(jié)構(gòu),并且破片對(duì)結(jié)構(gòu)的實(shí)際毀傷與預(yù)先打孔等方式存在著一定差距的情況。

    2)通過ConWep或壓力時(shí)程曲線實(shí)現(xiàn)沖擊波加載,采用Lagrange接觸算法模擬破片的侵徹作用[7-9]。這種方法一般很難用來求解結(jié)構(gòu)在近爆時(shí)的響應(yīng),因?yàn)榻瑫r(shí)沖擊波與結(jié)構(gòu)間有著較強(qiáng)的流固耦合效應(yīng),沖擊波載荷難以用常規(guī)的壓力時(shí)程曲線準(zhǔn)確表達(dá)。

    3)分別通過流固耦合方法和Lagrange接觸算法實(shí)現(xiàn)沖擊波載荷和破片載荷的加載[10-11]。目前,采用此種方法進(jìn)行的研究主要是針對(duì)單個(gè)破片和簡(jiǎn)單結(jié)構(gòu),較少涉及破片群和復(fù)雜結(jié)構(gòu),并且一般會(huì)忽略破片的加速過程,直接對(duì)破片施加初始速度。

    夾層結(jié)構(gòu)具有比強(qiáng)度高、抗沖擊性能好、可設(shè)計(jì)性強(qiáng)、隔熱、隔聲等優(yōu)良性能,在艦船結(jié)構(gòu)中有著廣闊的應(yīng)用前景。為探討沖擊波和破片群對(duì)鋼質(zhì)夾層結(jié)構(gòu)的聯(lián)合毀傷效應(yīng),本文將采用非線性有限元軟件LS-DYNA模擬I型夾層板在近場(chǎng)爆炸沖擊波和破片群聯(lián)合作用下的響應(yīng),分析炸藥當(dāng)量、芯層配置以及上、下面板厚度配置對(duì)其變形及破壞模式的影響,比較I型夾層板與等效實(shí)體板的抗毀傷性能,并分析I型夾層板的吸能特性。

    1 仿真方法和計(jì)算模型

    1.1 耦合及接觸算法

    在本文的計(jì)算中,炸藥和空氣采用ALE網(wǎng)格描述,預(yù)制破片、I型夾層板及等效實(shí)體板結(jié)構(gòu)采用Lagrange網(wǎng)格描述。為了模擬破片在沖擊波作用下的加速過程,以及夾層板在沖擊波與破片聯(lián)合作用下的響應(yīng),通過罰函數(shù)耦合算法定義破片、夾層板結(jié)構(gòu)與炸藥和空氣之間的流固耦合,使用侵蝕接觸定義破片與夾層板之間以及破片和夾層板自身間的接觸。

    LS-DYNA程序中罰函數(shù)耦合算法的基本原理是:如果在一個(gè)時(shí)間步內(nèi)有ALE流體材料穿透Lagrange結(jié)構(gòu)表面,則對(duì)相應(yīng)的結(jié)構(gòu)和流體分別施加耦合力以限制穿透并實(shí)現(xiàn)結(jié)構(gòu)與流體間力學(xué)參量的傳遞,其中,耦合力的大小與穿透深度成正比[12]。

    侵蝕接觸是侵徹模擬時(shí)常用的一種接觸類型。接觸通常用來處理Lagrange結(jié)構(gòu)之間的相互作用,其中最常用的接觸算法是罰函數(shù)法,其原理與上述罰函數(shù)耦合算法基本一致,主要區(qū)別在于一個(gè)用于Lagrange結(jié)構(gòu)間的接觸碰撞分析,另一個(gè)用于Lagrange結(jié)構(gòu)與ALE流體間的流固耦合分析[12]。本文使用的侵蝕接觸基于罰函數(shù)法,與其他接觸類型不同的是,侵蝕接觸可在表面單元失效刪除后繼續(xù)考慮剩余單元的接觸,可用來模擬破片侵徹問題。

    1.2 計(jì)算模型

    為分析I型夾層板在近距爆炸沖擊波與破片群聯(lián)合作用下的毀傷響應(yīng),本文選取了如圖1所示的計(jì)算模型。由于由15個(gè)胞元組成的I型夾層板結(jié)構(gòu)與5個(gè)胞元的類似,故圖中僅給出5個(gè)胞元的示意圖。I型夾層板結(jié)構(gòu)由上面板、下面板及芯層組成,材料均為304不銹鋼,特征參數(shù)有:夾層板長(zhǎng)度a、夾層板寬度b、上面板厚度tf、下面板厚度tb、胞元寬度bc、芯層高度Hc及芯層壁板厚度tc。本文所有I型夾層板的長(zhǎng)度a=400 mm,寬度b=400 mm,其余參數(shù)如表1所示,總的胞元個(gè)數(shù)為b/bc,邊界均為四周剛性固定。所有工況中,所用炸藥均為圓柱形TNT,裝藥位置為夾層板正中心,保持爆距(即炸藥底端面與夾層板上表面的距離)R=50 mm不變,均采用炸藥頂端中心起爆方式。通過在炸藥底端面布滿長(zhǎng)方體預(yù)制破片來模擬爆炸沖擊波和破片聯(lián)合作用載荷。預(yù)制破片采用穿甲能力較強(qiáng)的鎢合金材料,單個(gè)破片尺寸為6 mm×6 mm×3 mm,質(zhì)量約為1.92 g,破片數(shù)量由炸藥半徑?jīng)Q定(布滿炸藥底端面所需的最少破片),其平面布置方式如圖2所示。

    圖1 計(jì)算模型剖面示意圖(5個(gè)胞元)Fig.1 The cross section of calculation model(5 cells)

    圖2 不同炸藥半徑時(shí)破片平面布置示意圖Fig.2 The schematic layout of fragments at different radius of TNT explosive

    考慮到炸藥、破片及夾層板結(jié)構(gòu)(或等效實(shí)體板結(jié)構(gòu))的對(duì)稱性,為減少計(jì)算規(guī)模,進(jìn)行數(shù)值計(jì)算時(shí)只建立1/4模型,并在對(duì)稱面上設(shè)置對(duì)稱邊界條件。另外,由于近場(chǎng)爆炸所產(chǎn)生的沖擊波載荷具有局部特性,為節(jié)約計(jì)算時(shí)間,借鑒文獻(xiàn)[13]的做法,將空氣域的尺寸取為100 mm×100 mm× 260 mm,并在2個(gè)對(duì)稱面上施加對(duì)稱邊界條件,其余4個(gè)外表面施加無反射邊界條件。為了更好地反映結(jié)構(gòu)在爆炸沖擊波與破片聯(lián)合作用下的變形及損傷情況,夾層板、等效實(shí)體板及預(yù)制破片均采用三維實(shí)體單元Solid 164建模,并采用Lagrange算法,炸藥和空氣采用ALE算法。夾層板、等效實(shí)體板及預(yù)制破片基本采用1 mm的大小進(jìn)行網(wǎng)格劃分,但對(duì)于2 mm厚的上、下面板和芯層,調(diào)整其厚度方向份數(shù)為3份,以便更好地模擬其彎曲和侵徹效果且不會(huì)導(dǎo)致因網(wǎng)格尺寸過小而限制求解時(shí)間步長(zhǎng)??諝庥虿捎? mm的大小進(jìn)行均勻網(wǎng)格劃 分 ,并 通 過 正確 設(shè) 置*INITIAL_VOL?UME_FRACTION_GEOMETRY關(guān)鍵字完成圓柱形炸藥的填充。以上、下面板厚度均為3 mm的I型夾層板為例,其整體計(jì)算模型如圖3所示。

    表1 I型夾層板計(jì)算工況及仿真結(jié)果Tab.1 Computational conditions and numerical results of I-core sandwich panels

    圖3 I型夾層板計(jì)算模型示意圖(1/4模型)Fig.3 Schematic of calculation model for I-core sandwich panel(one quarter of model)

    1.3 材料模型

    TNT炸藥采用 LS-DYNA程序中的*MAT_HIGH_EXPLOSIVE_BURN模型及JWL狀態(tài)方程描述,空氣采用*MAT_NULL材料模型及*EOS_LINEAR_POLYNOMIAL狀態(tài)方程描述,具體參數(shù)見文獻(xiàn)[14]。

    夾層板及實(shí)體板結(jié)構(gòu)材料(304不銹鋼)采用爆炸沖擊問題中常用的Johnson-Cook材料模型及等效塑性應(yīng)變失效準(zhǔn)則描述,具體參數(shù)見文獻(xiàn)[15],失效應(yīng)變?nèi)?.42。

    預(yù)制破片采用預(yù)制破片戰(zhàn)斗部常用的鎢合金材料。為保守估計(jì)預(yù)制鎢合金破片的毀傷性能并簡(jiǎn)化計(jì)算,本文忽略了破片在加速過程以及侵徹靶板過程中的變形或損傷,采用剛體材料模型來描述,密度為17 800 kg/m3,彈性模量為357 GPa,泊松比為0.303。

    1.4 計(jì)算工況

    本文考慮表1所示的計(jì)算工況。表中:工況I-1~I(xiàn)-10和I-15~I(xiàn)-22中I型夾層板的結(jié)構(gòu)參數(shù)僅上、下面板厚度配置可能不同,后文將以“tf-tb”來表示這些工況下上、下面板厚度分別為tf和tb的I型夾層板;工況I-11~I(xiàn)-14中的夾層板具有不同的芯層配置,由15個(gè)胞元構(gòu)成;M為夾層板質(zhì)量,本文所有夾層板結(jié)構(gòu)質(zhì)量均為9.354 kg;TNT當(dāng)量W=157,74,35 g對(duì)應(yīng)的圓柱形TNT尺寸分別為(d=70 mm,L=25 mm),(d=48 mm,L=25 mm)和(d=48 mm,L=12 mm);n=0表示裸裝藥工況(沖擊波單獨(dú)作用),n>0表示聯(lián)合作用工況(沖擊波與破片聯(lián)合作用);破口表示穿孔破壞,損傷表示表層有單元失效刪除但沒有形成貫穿破口;下標(biāo)f,b,c分別表示上面板、下面板及芯層。根據(jù)I型夾層板的結(jié)構(gòu)尺寸,為保證面內(nèi)尺寸不變,由質(zhì)量相等原則可計(jì)算出等效實(shí)體板的厚度為7.4 mm。為對(duì)比I型夾層板與等效實(shí)體板的抗毀傷能力,對(duì)等效實(shí)體板在相同載荷下的毀傷進(jìn)行了計(jì)算,具體計(jì)算工況如表2所示。

    2 仿真方法驗(yàn)證

    在用仿真方法研究參數(shù)影響規(guī)律之前,一般需要根據(jù)試驗(yàn)數(shù)據(jù)對(duì)所采用的方法進(jìn)行有效性驗(yàn)證。但由于目前關(guān)于結(jié)構(gòu)在沖擊波與破片聯(lián)合作用下的試驗(yàn)數(shù)據(jù)較少,因此本文僅對(duì)ALE流固耦合算法計(jì)算裸裝藥爆炸問題、Lagrange接觸算法計(jì)算侵徹問題進(jìn)行驗(yàn)證。

    2.1 裸裝藥爆炸下實(shí)體板響應(yīng)仿真驗(yàn)證

    Zhang等[15]針對(duì)三角形波紋金屬夾層板及等效實(shí)體板在空中近場(chǎng)爆炸載荷下的響應(yīng)進(jìn)行了試驗(yàn)研究。本文選取 Zhang等給出的典型等效實(shí)體板試驗(yàn)工況進(jìn)行仿真計(jì)算,以驗(yàn)證前述ALE流固耦合算法計(jì)算裸裝藥近場(chǎng)爆炸問題的準(zhǔn)確性。

    試驗(yàn)用等效實(shí)體板的有效尺寸為300 mm×288 mm,板厚為3.753 mm,邊界為四邊固支,材料為304不銹鋼。選取的載荷工況為:TNT當(dāng)量55 g,爆距50 mm。根據(jù)結(jié)構(gòu)及載荷的對(duì)稱性,只建立1/4有限元模型,并在對(duì)稱面上施加對(duì)稱邊界條件。其中,空氣域的面內(nèi)尺寸70 mm×70 mm是根據(jù)文獻(xiàn)[13]的建議并綜合計(jì)算精度和計(jì)算時(shí)間而選取的。實(shí)體板結(jié)構(gòu)和空氣域均采用三維實(shí)體單元Solid 164離散,其中實(shí)體板采用1 mm的網(wǎng)格,空氣域采用1.5 mm的網(wǎng)格。通過體積初始化方法填充炸藥,并使用罰函數(shù)耦合算法定義平板與流體(炸藥和空氣)之間的耦合。另外,為了控制大變形問題所用單點(diǎn)積分單元帶來的沙漏模態(tài),必須選取合適的沙漏控制方法及沙漏系數(shù)。仿真中所用材料模型及材料參數(shù)詳見第1.3節(jié)。

    表2 等效實(shí)體板計(jì)算工況及仿真結(jié)果Tab.2 Computational conditions and numerical results of equivalent solid plates

    圖4給出了試驗(yàn)與仿真得到的實(shí)體板中心橫截面輪廓示意圖。從中可以看出,仿真計(jì)算得到的實(shí)體板變形模式與試驗(yàn)結(jié)果具有較好的一致性,均為整體變形疊加局部大變形。計(jì)算得到實(shí)體板中心點(diǎn)的最終塑性變形為28.9 mm,與試驗(yàn)值間的誤差為13%(此工況下中心點(diǎn)變形試驗(yàn)值為33.24 mm)。由定性和定量結(jié)果可以看出,本文采取的ALE流固耦合算法可以較準(zhǔn)確地預(yù)測(cè)結(jié)構(gòu)在近距爆炸沖擊波下的響應(yīng)。

    圖4 實(shí)體板中心橫截面輪廓示意圖Fig.4 Schematic of center cross section profile of solid plate

    2.2 實(shí)體板侵徹仿真驗(yàn)證

    梅志遠(yuǎn)等[16]對(duì)破片侵徹945鋼靶板進(jìn)行了試驗(yàn)研究。為方便起見,本文僅選取典型試驗(yàn)工況進(jìn)行仿真計(jì)算,以驗(yàn)證前述Lagrange接觸算法計(jì)算侵徹問題的可行性。

    試驗(yàn)用靶板尺寸為80 cm×80 cm,板厚4 mm,材料為945鋼。破片為邊長(zhǎng)7.5 mm的立方體,材料為45鋼,質(zhì)量3.2 g,入射速度為1 066 m/s。同樣,根據(jù)正面垂直侵徹時(shí)結(jié)構(gòu)和載荷的對(duì)稱性,只建立1/4有限元模型,并在對(duì)稱面上施加對(duì)稱邊界條件。破片和靶板均采用三維實(shí)體單元Solid 164離散,網(wǎng)格大小為1 mm。采用侵蝕接觸考慮破片與靶板間的相互作用,并設(shè)置接觸剛度縮放因子為0.5,以得到合適的接觸剛度。仿真中,破片和靶板均采用Johnson-Cook材料模型,具體材料參數(shù)詳見文獻(xiàn)[16]。

    圖5所示為仿真得到的破片速度時(shí)程曲線。由圖可知,破片的剩余速度為490.4 m/s,與試驗(yàn)值的誤差為9%(試驗(yàn)測(cè)得的剩余速度為449.5 m/s)。可見,本文采用的Lagrange接觸算法可以較準(zhǔn)確地計(jì)算破片對(duì)靶板的侵徹效應(yīng)。

    圖5 破片采用不同材料模型時(shí)破片速度時(shí)程曲線Fig.5 Velocity history curves of fragment with different material models for fragment

    圖6(a)給出了靶板的破壞模式。從中可明顯看出,破片穿透靶板后,靶板有一塊被沖出,為典型的沖塞破壞模式。同時(shí)在此過程中,破片本身也出現(xiàn)了一定程度的變形和侵蝕。為了簡(jiǎn)便起見,也可將破片當(dāng)成剛體處理。圖5和圖6(b)分別給出了破片采用剛體模型時(shí)破片的速度時(shí)程曲線以及靶板的破壞模式。通過比較可知,破片采用剛體模型時(shí),其穿透靶板后剩余速度較大,但靶板的破壞模式基本相同??梢?,將破片當(dāng)成剛體處理是一種較保守的方式,本文在沖擊波與破片聯(lián)合作用分析時(shí)也將采取這種簡(jiǎn)化方式(即將預(yù)制鎢合金破片當(dāng)成剛體)。

    圖6 破片采用不同材料模型時(shí)靶板破壞模式Fig.6 Failure modes of target plate with different material models for fragment

    3 計(jì)算結(jié)果分析

    在分別驗(yàn)證ALE流固耦合算法計(jì)算裸裝藥爆炸問題和Lagrange接觸算法計(jì)算侵徹問題有效性的基礎(chǔ)上,本文采用第1.1節(jié)中所述方法計(jì)算了I型夾層板及等效實(shí)體板在沖擊波載荷及沖擊波與破片聯(lián)合作用載荷下的毀傷響應(yīng)。

    3.1 I型夾層板動(dòng)態(tài)響應(yīng)過程

    沖擊波單獨(dú)作用與沖擊波和破片聯(lián)合作用下夾層板的響應(yīng)過程有些差異,故本節(jié)選取典型工況I-9和I-10進(jìn)行說明。

    圖7給出了裸裝藥爆炸下沖擊波的傳播及I型夾層板的響應(yīng)過程。初始時(shí)刻,空氣具有初始?jí)毫?.1 MPa(圖7(a))。炸藥起爆后,在極短的時(shí)間內(nèi)轉(zhuǎn)變?yōu)楦邷?、高壓的爆轟產(chǎn)物,并急劇膨脹,強(qiáng)烈擠壓周圍空氣,在空氣中形成爆炸空氣沖擊波。炸藥起爆后約15 μs(圖7(b)),沖擊波傳播到夾層板上表面并與其發(fā)生相互作用,使上面板獲得動(dòng)量。隨著時(shí)間的推進(jìn),上面板中心區(qū)域出現(xiàn)凹陷變形,并擠壓芯層壁板,芯層發(fā)生塑性變形(圖7(c)~圖7(d))。大約200 μs后,由于壓力較低,耦合作用很弱,夾層板結(jié)構(gòu)主要在其自身慣性力的作用下繼續(xù)變形、振蕩。

    圖7 沖擊波單獨(dú)作用下I型夾層板響應(yīng)過程(工況I-10)Fig.7 The response process of I-core sandwich panel subjected to blast loading(case I-10)

    圖8給出了底部布有預(yù)制破片的炸藥爆炸后破片的運(yùn)動(dòng)及I型夾層板的毀傷過程。初始時(shí)刻,預(yù)制破片緊貼炸藥底端面,并與夾層板上表面間有一定的距離(圖8(a))。炸藥起爆后,破片在爆轟產(chǎn)物的作用下開始加速,約70 μs時(shí),位于中心部位的破片開始接觸上面板(圖8(b))。因速度較高,中心部位破片逐漸穿透上面板,并且上面板形成部分二次碎片(圖8(c))。由于靠近外圍的破片獲得的速度較低并且撞擊上面板時(shí)入射角較大,因此沒能穿透上面板而是發(fā)生跳飛現(xiàn)象。穿過上面板的破片以及上面板二次碎片繼續(xù)向下運(yùn)動(dòng),對(duì)下面板進(jìn)行作用并部分穿透下面板(圖8(d))。整個(gè)過程中,沖擊波除了加速破片之外,也對(duì)夾層板結(jié)構(gòu)產(chǎn)生作用,加劇了其毀傷程度。

    圖8 沖擊波和破片聯(lián)合作用下I型夾層板響應(yīng)過程(工況I-9)Fig.8 The response process of I-core sandwich panel subjected to combined loading(case I-9)

    3.2 I型夾層板失效模式

    圖9和圖10分別給出了沖擊波單獨(dú)作用及沖擊波與破片聯(lián)合作用下上、下面板厚度均為3 mm的I型夾層板的失效模式。通過對(duì)比兩圖可以明顯看出,夾層板在沖擊波與破片聯(lián)合作用下的毀傷程度均大于同等當(dāng)量沖擊波單獨(dú)作用下夾層板的毀傷,并且不同當(dāng)量下夾層板的失效模式也有差異。

    圖9 不同TNT當(dāng)量沖擊波單獨(dú)作用下I型夾層板的失效模式(tf=3,tb=3)Fig.9 The failure modes of I-core sandwich panel at different TNT charges subjected to blast loading(tf=3,tb=3)

    圖10 不同TNT當(dāng)量沖擊波與破片聯(lián)合作用下I型夾層板的失效模式(tf=3,tb=3)Fig.10 The failure modes of I-core sandwich panel at different TNT charges subjected to combined loading(tf=3,tb=3)

    沖擊波單獨(dú)作用下,當(dāng)TNT當(dāng)量W=35,74 g時(shí),夾層板上面板僅發(fā)生塑性大變形,沒有形成破口;當(dāng)TNT當(dāng)量W=157 g時(shí),夾層板上面板沿著與芯層連接的部位撕開且中心部位在沖擊波的作用下形成了二次碎片(圖11),下面板在沖擊波及上面板二次碎片的共同作用下產(chǎn)生塑性大變形及局部表層損傷。

    圖11 157 g當(dāng)量沖擊波單獨(dú)作用下夾層板上面板失效模式(工況I-8)Fig.11 Failure mode of front plate of the I-core sandwich panel subjected to blast loading with 157 g TNT(case I-8)

    在沖擊波與破片的聯(lián)合作用下,當(dāng)TNT當(dāng)量W=35 g時(shí),上面板出現(xiàn)中心大破口、小穿孔及局部表層損傷等破壞模式并形成二次碎片(圖12),下面板在沖擊波、預(yù)制破片及上面板二次碎片的共同作用下產(chǎn)生塑性大變形及局部表層損傷;當(dāng)TNT當(dāng)量W=74,157 g時(shí),夾層板上、下面板均出現(xiàn)中心大破口、局部表層損傷及局部塑性大變形,并且在W=157 g時(shí)上、下面板的損傷范圍更大(中心兩側(cè)的胞元也有較大損傷),這主要是因?yàn)樵诖朔N當(dāng)量下炸藥半徑更大、破片數(shù)量更多。

    圖12 35 g當(dāng)量沖擊波與破片聯(lián)合作用下夾層板上面板失效模式(工況I-15)Fig.12 Failure mode of front plate of the I-core sandwich panel subjected to combined loading with 35 g TNT(case I-15)

    為研究芯層配置對(duì)I型夾層板失效模式的影響,在原有5個(gè)胞元結(jié)構(gòu)(工況I-9和I-10)的基礎(chǔ)上,設(shè)計(jì)了2種含有15個(gè)胞元的等重夾層板結(jié)構(gòu)(一種是芯層高度Hc=70 mm不變,根據(jù)質(zhì)量相等計(jì)算得到芯層壁板厚度tc=0.571 mm;另一種是芯層高度減為一半,即Hc=35 mm,同樣根據(jù)質(zhì)量相等計(jì)算得到芯層壁板厚度tc=1.143 mm),并計(jì)算了其在74 g當(dāng)量沖擊波單獨(dú)作用及沖擊波與破片聯(lián)合作用下的毀傷響應(yīng)(工況I-11~I(xiàn)-14)。

    圖13和圖14分別給出了TNT當(dāng)量W=74 g時(shí)沖擊波單獨(dú)作用及沖擊波與破片聯(lián)合作用下3種不同芯層配置I型夾層板的失效模式。沖擊波單獨(dú)作用下,在工況I-10和I-12中,夾層板上面板均沒有破壞而是產(chǎn)生了塑性大變形,并且工況I-12中夾層板芯層的屈曲變形由于剛度的減小而更加明顯(圖13(b));工況I-14中,夾層板由于芯層剛度較大,在上面板與中心芯材連接部位附近產(chǎn)生了局部破壞(圖13(c))。另外,由于芯層高度減小,I型夾層板整體抗彎剛度也大幅減小,下面板中心點(diǎn)變形增大,與工況I-10相比,I-14中夾層板下面板中心點(diǎn)的最大變形增大了58.4%(變形數(shù)據(jù)見表1)。

    圖13 74 g當(dāng)量沖擊波單獨(dú)作用下不同芯層配置I型夾層板的失效模式Fig.13 The failure modes of I-core sandwich panel with different core configuration subjected to blast loading with 74 g TNT

    在沖擊波與破片聯(lián)合作用下,3種芯層配置I型夾層板的上、下面板均產(chǎn)生了中心破口(圖14)。在工況I-11和I-13中,夾層板由于芯層比較密集,故增大了破片撞擊芯層的概率,提高了芯層的參與程度,并減小了中心破口大小。

    圖14 74 g當(dāng)量沖擊波與破片聯(lián)合作用下不同芯層配置I型夾層板的失效模式Fig.14 The failure modes of I-core sandwich panel with different core configuration subjected to combined loading with 74 g TNT

    3.3 上、下面板厚度配置對(duì)I型夾層板失效模式影響分析

    為了分析上、下面板厚度配置對(duì)I型夾層板失效模式的影響,本文對(duì)3種較典型厚度配置的夾層板進(jìn)行了計(jì)算(工況I-1~I(xiàn)-10及I-15~I(xiàn)-22)。其中,上、下面板厚度均為3 mm(簡(jiǎn)記為“3-3”)的I型夾層板毀傷結(jié)果已在第3.2節(jié)中給出,上、下面板厚度分別為4,2 mm(簡(jiǎn)記為“4-2”)以及2,4 mm(簡(jiǎn)記為“2-4”)的I型夾層板毀傷結(jié)果如圖15~圖18所示。各工況下夾層板上、下面板的變形及塑性應(yīng)變結(jié)果見表1。

    圖15 不同TNT當(dāng)量沖擊波單獨(dú)作用下I型夾層板的失效模式(tf=4,tb=2)Fig.15 The failure modes of I-core sandwich panel at different TNT charges subjected to blast loading(tf=4,tb=2)

    圖16 不同TNT當(dāng)量沖擊波與破片聯(lián)合作用下I型夾層板的失效模式(tf=4,tb=2)Fig.16 The failure modes of I-core sandwich panel at different TNT charges subjected to combined loading(tf=4,tb=2)

    圖17 不同TNT當(dāng)量沖擊波單獨(dú)作用下I型夾層板的失效模式(tf=2,tb=4)Fig.17 The failure modes of I-core sandwich panel at different TNT charges subjected to blast loading(tf=2,tb=4)

    圖18 不同TNT當(dāng)量沖擊波與破片聯(lián)合作用下I型夾層板的失效模式(tf=2,tb=4)Fig.18 The failure modes of I-core sandwich panel at different TNT charges subjected to combined loading(tf=2,tb=4)

    對(duì)比同一尺寸夾層板在相同當(dāng)量下下面板的變形或破壞程度可以明顯看出,與沖擊波單獨(dú)作用相比,在沖擊波與破片聯(lián)合作用下,夾層板下面板的毀傷更為嚴(yán)重。其中,TNT當(dāng)量W=35 g時(shí),與沖擊波單獨(dú)作用相比,在沖擊波與破片聯(lián)合作用下,“4-2”,“3-3”及“2-4”夾層板下面板中心點(diǎn)的最大變形分別增大了389.3%,352.4%和375.6%。

    在沖擊波單獨(dú)作用下,當(dāng)TNT當(dāng)量W=35,74 g時(shí),3種厚度配置夾層板的下面板均沒有損傷,只發(fā)生塑性變形,并且“2-4”夾層板下面板中心點(diǎn)的變形最?。ū?)。這主要是由于本文分析的I型夾層板芯層高度較大,在上面板未破壞時(shí),上面板中心區(qū)域在變形過程中不會(huì)接觸下面板,下面板的變形主要由中心附近I型芯層傳遞的力造成,因此,當(dāng)下面板較厚時(shí)中心點(diǎn)的變形較小。需要指出的是,如果夾層板的芯層高度、芯層間距或芯層拓?fù)湫褪桨l(fā)生變化,以致上面板中心區(qū)域在變形過程中接觸到了下面板,那么這一結(jié)論就有可能不再成立[15]。

    在沖擊波單獨(dú)作用下,當(dāng)TNT當(dāng)量W=157 g時(shí),“2-4”夾層板下面板發(fā)生局部損傷并且與芯層的連接處被撕裂(圖17(c)),“3-3”夾層板下面板發(fā)生局部損傷(圖9(c)),“4-2”夾層板下面板沒有損傷,只有塑性變形(圖15(c))。造成這種差異的主要原因是當(dāng)上面板較薄時(shí),在沖擊波的作用下,上面板容易發(fā)生材料破碎而形成較多的二次碎片(圖19),二次碎片對(duì)下面板有著較強(qiáng)的毀傷作用;當(dāng)上面板較厚時(shí),在沖擊波的作用下,上面板主要發(fā)生撕裂破壞,形成的二次碎片較少(圖20),因此此時(shí)下面板局部損傷較小或沒有局部損傷。

    圖19 157 g當(dāng)量沖擊波單獨(dú)作用下夾層板上面板失效模式(工況I-18)Fig.19 Failure mode of front plate of the I-core sandwich panel subjected to blast loading with 157 g TNT(case I-18)

    圖20 157 g當(dāng)量沖擊波單獨(dú)作用下夾層板上面板失效模式(工況I-2)Fig.20 Failure mode of front plate of the I-core sandwich panel subjected to blast loading with 157 g TNT(case I-2)

    在沖擊波與破片聯(lián)合作用下,當(dāng)TNT當(dāng)量W=74,157 g時(shí),3種厚度配置夾層板的上、下面板均產(chǎn)生了中心破口。在沖擊波與破片聯(lián)合作用下,當(dāng)TNT當(dāng)量W=35 g時(shí),3種厚度配置夾層板下面板均沒有產(chǎn)生中心破口,但“3-3”和“2-4”夾層板下面板均有一定的表層損傷,而“4-2”夾層板下面板則沒有局部損傷,只發(fā)生了塑性大變形。產(chǎn)生這種差異的主要原因是當(dāng)上面板較薄時(shí),在沖擊波與破片聯(lián)合作用下,上面板破碎較嚴(yán)重,形成的二次碎片較小(圖21),并且預(yù)制破片穿透上面板后的剩余速度較大,對(duì)下面板造成的毀傷更為嚴(yán)重;而當(dāng)上面板較厚時(shí),在沖擊波與破片聯(lián)合作用下,上面板形成了較大的二次碎片(即沖塞塊),預(yù)制破片穿過上面板后剩余速度較小,并且大塊二次碎片與下面板大面積接觸,基本阻擋了預(yù)制破片與下面板的直接小面積接觸(圖22),因此,此時(shí)下面板的局部損傷較小或沒有局部損傷,但塑性變形較大。

    圖21 35 g當(dāng)量沖擊波與破片聯(lián)合作用下夾層板上面板失效模式(工況I-21)Fig.21 Failure mode of front plate of the I-core sandwich panel subjected to combined loading with 35 g TNT(case I-21)

    圖22 35 g當(dāng)量沖擊波與破片聯(lián)合作用下夾層板上面板失效模式(工況I-5)Fig.22 Failure mode of front plate of the I-core sandwich panel subjected to combined loading with 35 g TNT(case I-5)

    3.4 I型夾層板與等效實(shí)體板抗毀傷性能對(duì)比

    為了對(duì)比I型夾層板與等效實(shí)體板的抗毀傷性能,計(jì)算了等效實(shí)體板在相同載荷工況下的毀傷情況,如圖23和圖24所示。1.4節(jié)中的表2列出了各工況下等效實(shí)體板的變形及等效塑性應(yīng)變結(jié)果。結(jié)合圖表可以發(fā)現(xiàn),3種當(dāng)量下沖擊波單獨(dú)作用時(shí),等效實(shí)體板均只產(chǎn)生了塑性大變形,并且塑性應(yīng)變及中心點(diǎn)的變形是隨當(dāng)量的增加而增大。當(dāng)沖擊波與破片聯(lián)合作用時(shí),實(shí)體板在35 g當(dāng)量下僅產(chǎn)生塑性大變形及局部表層損傷,沒有形成破口,在74 g和157 g當(dāng)量下,實(shí)體板均產(chǎn)生中心破口,并且在157 g當(dāng)量時(shí)中心破口更大。

    圖23 不同TNT當(dāng)量沖擊波單獨(dú)作用下等效實(shí)體板的失效模式Fig.23 The failure modes of equivalent solid plate at different TNT charges subjected to blast loading

    圖24 不同TNT當(dāng)量沖擊波與破片聯(lián)合作用下等效實(shí)體板的失效模式Fig.24 The failure modes of equivalent solid plate at different TNT charges subjected to combined loading

    仔細(xì)對(duì)比分析夾層板及等效實(shí)體板的計(jì)算結(jié)果(表1和表2)可以發(fā)現(xiàn),沖擊波單獨(dú)作用時(shí),在35 g和74 g當(dāng)量下,3種厚度配置夾層板的下面板及等效實(shí)體板均只產(chǎn)生塑性變形,并且3種夾層板下面板的變形均小于實(shí)體板,其中“2-4”夾層板的變形最小。若以實(shí)體板中心點(diǎn)的變形為基準(zhǔn),則35 g當(dāng)量下,“4-2”,“3-3”及“2-4”夾層板下面板中心點(diǎn)的最大變形分別減小了71.7%,73.3%和77.7%;在74 g當(dāng)量下,“4-2”,“3-3”及“2-4”夾層板下面板中心點(diǎn)的最大變形分別減小了67.1%,72.4%和79.5%。在157 g當(dāng)量下,“3-3”和“2-4”夾層板的下面板均有損傷,其抗毀傷性能不如等效實(shí)體板(因?yàn)樵诖溯d荷下等效實(shí)體板只有塑性變形沒有損傷),而“4-2”夾層板的下面板只有塑性變形,并且其變形比實(shí)體板小55.0%,因此在此種載荷下,“4-2”夾層板的抗毀傷性能要優(yōu)于等效實(shí)體板。實(shí)際上,若TNT當(dāng)量繼續(xù)增大,夾層板和等效實(shí)體板可能會(huì)出現(xiàn)僅有一個(gè)發(fā)生破口的臨界狀態(tài)。Zhang等[15]通過試驗(yàn)研究發(fā)現(xiàn),當(dāng)沖擊波載荷較強(qiáng)時(shí),三角形波紋夾層板會(huì)發(fā)生破口,但實(shí)體板則僅產(chǎn)生塑性變形,此時(shí)夾層板的抗毀傷能力不如實(shí)體板。

    沖擊波與破片聯(lián)合作用時(shí),在74 g和157 g當(dāng)量下,3種夾層板的下面板及等效實(shí)體板均產(chǎn)生了中心破口,不易比較抗毀傷性能的優(yōu)劣;在35 g當(dāng)量下,等效實(shí)體板產(chǎn)生局部損傷,“3-3”和“2-4”夾層板的下面板也有損傷但其變形分別比等效實(shí)體板小20.1%和29.7%,“4-2”夾層板的下面板沒有損傷,只有塑性變形,因此在此種載荷下,3種夾層板的抗毀傷性能均優(yōu)于實(shí)體板,并且“4-2”夾層板的性能最好。

    3.5 I型夾層板吸能特性分析

    圖25給出了不同當(dāng)量沖擊波單獨(dú)作用下3種不同上、下面板厚度配置的I型夾層板及等效實(shí)體板的總吸能量對(duì)比。由圖可以明顯看出,在同種當(dāng)量下,3種I型夾層板的吸能量均高于實(shí)體板,體現(xiàn)了沖擊波單獨(dú)作用下I型夾層板在能量吸收效率上的優(yōu)勢(shì)。若以等效實(shí)體板吸能量為基準(zhǔn),則在35 g當(dāng)量下,“4-2”,“3-3”及“2-4”夾層板的吸能量分別提高了95.9%,158.7%和281.7%;在74 g當(dāng)量下,“4-2”,“3-3”及“2-4”夾層板的吸能量分別提高了86.7%,146.1%和211.5%;在157 g當(dāng)量下,“4-2”,“3-3”及“2-4”夾層板的吸能量分別提高了65.6%,44.6%和22.2%。在大當(dāng)量(157 g)下,I型夾層板的吸能優(yōu)勢(shì)有所減小,這是因?yàn)殡S著載荷強(qiáng)度的提高,夾層板的上面板因厚度較薄容易過早發(fā)生破壞,導(dǎo)致部分材料沒有充分變形,從而減少了吸能。

    圖25 沖擊波單獨(dú)作用下夾層板與實(shí)體板總吸能Fig.25 Total energy absorbed by sandwich and equivalent solid plate under blast loading

    不同當(dāng)量沖擊波單獨(dú)作用下夾層板各部件(上面板、芯層和下面板)的吸能特性如圖26所示。可以明顯看出,上面板和芯層是主要的吸能構(gòu)件,如“3-3”夾層板在74 g當(dāng)量下上面板、芯層及下面板的吸能比例分別為77.4%,21.0%和1.6%。在157 g當(dāng)量時(shí),由于載荷強(qiáng)度過大,夾層板的上面板過早發(fā)生了破壞,使得沖擊波及上面板的二次碎片直接作用于下面板,從而導(dǎo)致下面板變形增大、吸能比例上升。

    圖26 沖擊波單獨(dú)作用下夾層板各部件吸能Fig.26 Energy absorbed by different parts of sandwich panel under blast loading

    圖27給出了不同當(dāng)量沖擊波與破片聯(lián)合作用下3種不同上、下面板厚度配置的I型夾層板及等效實(shí)體板的總吸能量對(duì)比。從中可以看出,此時(shí)I型夾層板與等效實(shí)體板的總吸能量差別較小。在沖擊波與破片聯(lián)合作用下,I型夾層板各部件的吸能特性如圖28所示。此時(shí),上面板和下面板是主要的吸能構(gòu)件,并且“4-2”及“2-4”夾層板上面板的吸能比例分別為最高和最低。

    圖29給出了3種不同芯層配置I型夾層板的吸能特性。從中可以看出,芯層加密后,芯層的吸能比例提高了(比如工況I-10,I-12及I-14中夾層板芯層的吸能比例分別為21.0%,38.8%和46.9%),這與前面失效模式的分析結(jié)果一致。

    圖27 沖擊波與破片聯(lián)合作用下夾層板與實(shí)體板吸能Fig.27 Total energy absorbed by sandwich and equivalent solid plate under combined loading

    圖28 沖擊波與破片聯(lián)合作用下夾層板各部件吸能Fig.28 Energy absorbed by different parts of sandwich panel under combined loading

    圖29 不同芯層配置夾層板各部件吸能Fig.29 Energy absorbed by different parts of sandwich panel with different core configuration

    4 結(jié) 論

    本文通過在炸藥底部布置預(yù)制破片來模擬戰(zhàn)斗部爆炸產(chǎn)生的沖擊波與破片聯(lián)合作用載荷,利用動(dòng)力分析有限元軟件LS-DYNA計(jì)算了3種TNT炸藥當(dāng)量下I型夾層板的毀傷響應(yīng),分析了沖擊波單獨(dú)作用及沖擊波與破片聯(lián)合作用下夾層板失效模式的差異,研究了夾層板芯層配置以及上、下面板厚度配置對(duì)其失效模式的影響,并與等效實(shí)體板的抗毀傷性能進(jìn)行了對(duì)比。同時(shí),還從吸能的角度分析了不同載荷工況下I型夾層板的吸能特性,得出如下結(jié)論:

    1)與沖擊波單獨(dú)作用相比,同當(dāng)量沖擊波與破片的聯(lián)合作用對(duì)實(shí)體板和I型夾層板結(jié)構(gòu)造成的毀傷更為嚴(yán)重。當(dāng)TNT當(dāng)量W=35 g時(shí),與沖擊波單獨(dú)作用相比,沖擊波與破片聯(lián)合作用下“4-2”,“3-3”及“2-4”夾層板下面板中心點(diǎn)的最大變形分別增大了389.3%,352.4%和375.6%;當(dāng)TNT當(dāng)量W=74,157 g時(shí),夾層板下面板在沖擊波單獨(dú)作用下沒有產(chǎn)生中心破口,而在沖擊波與破片聯(lián)合作用下則產(chǎn)生了較大的中心破口。

    2)載荷強(qiáng)度與I型夾層板結(jié)構(gòu)的抗毀傷性能密切相關(guān)。當(dāng)載荷強(qiáng)度較小(TNT當(dāng)量較?。r(shí),I型夾層板的抗毀傷性能優(yōu)于等效實(shí)體板(當(dāng)TNT當(dāng)量W=35 g時(shí),在沖擊波的單獨(dú)作用下,“4-2”,“3-3”及“2-4”夾層板下面板中心點(diǎn)的最大變形分別比等效實(shí)體板小71.7%,73.3%和77.7%;而在沖擊波與破片聯(lián)合作用下,“4-2”,“3-3”及“2-4”夾層板下面板中心點(diǎn)的最大變形分別比等效實(shí)體板小8.3%,20.1%和29.7%);而當(dāng)載荷強(qiáng)度較大(TNT當(dāng)量較大)時(shí),I型夾層板及等效實(shí)體板均會(huì)在中心區(qū)域產(chǎn)生破口失效。

    3)載荷強(qiáng)度、載荷類型(沖擊波單獨(dú)作用或沖擊波與破片聯(lián)合作用)及上、下面板厚度配置對(duì)I型夾層板的失效模式有較大影響。在沖擊波單獨(dú)作用下,I型夾層板的上面板在厚度較薄時(shí)容易出現(xiàn)大范圍的材料破碎而形成較多的二次碎片,從而影響I型夾層板的抗毀傷性能;而在厚度較厚時(shí)則會(huì)使得上面板主要出現(xiàn)撕裂破壞,能在一定程度上改善I型夾層板的抗毀傷性能。在沖擊波與破片聯(lián)合作用下,當(dāng)上面板的厚度較薄時(shí),破片容易穿透上面板而直接對(duì)下面板造成損傷;而當(dāng)上面板的厚度較厚、TNT當(dāng)量較小時(shí),在破片的作用下,上面板會(huì)出現(xiàn)沖塞破壞,下面板在大塊二次碎片和破片的共同作用下會(huì)發(fā)生較大的塑性變形,但可能不會(huì)出現(xiàn)損傷現(xiàn)象。

    4)從吸能特性分析結(jié)果來看,在沖擊波單獨(dú)作用下,上面板和芯層是主要的吸能構(gòu)件,并且隨著載荷強(qiáng)度的增大,上面板容易出現(xiàn)破壞,導(dǎo)致沖擊波和上面板的二次碎片直接作用于下面板,從而使下面板變形增大、吸能比例上升。在沖擊波與破片聯(lián)合作用下,上面板和下面板是主要的吸能構(gòu)件,并且加密芯材能夠增加破片撞擊芯層的概率,在一定程度上提高了芯層的吸能比例。

    [1] 孔祥韶.爆炸載荷及復(fù)合多層防護(hù)結(jié)構(gòu)響應(yīng)特性研究[D].武漢:武漢理工大學(xué),2013.

    [2] 侯海量,朱錫,李偉,等.爆炸沖擊波和高速破片聯(lián)合作用下艙室結(jié)構(gòu)破壞模式試驗(yàn)研究[C]//中國(guó)鋼結(jié)構(gòu)協(xié)會(huì)海洋鋼結(jié)構(gòu)分會(huì)2010年學(xué)術(shù)會(huì)議暨第六屆理事會(huì)第三次會(huì)議論文集,2010.

    [3] 張成亮,朱錫,侯海量,等.爆炸沖擊波與高速破片對(duì)夾層結(jié)構(gòu)的聯(lián)合毀傷效應(yīng)試驗(yàn)研究[J].振動(dòng)與沖擊,2014,33(15):184-188. ZHANG Chengliang,ZHU Xi,HOU Hailiang,et al. Tests for combined damage effect of blast waves and high-velocity fragments on composite sandwich plates[J].Journal of Vibration and Shock,2014,33(15):184-188.

    [4] 金乾坤.破片和沖擊波毀傷圓柱靶的數(shù)值仿真[J].兵工學(xué)報(bào),2006,27(2):215-218. JIN Qiankun.Simulation of cylindrical shell damage by fragments and shock waves[J].Acta Armamentarii,2006,27(2):215-218.

    [5] RAKV?G K G,UNDERWOOD N J,SCHLEYER G K,et al.Transient pressure loading of clamped metal?lic plates with preformed holes[J].International Jour?nal of Impact Engineering,2013,53:44-55.

    [6] 姚志敏,雷灝,尉廣軍,等.破片和沖擊波復(fù)合作用下靶板毀傷仿真[J].火力與指揮控制,2014,39(7):152-154. YAO Zhimin,LEI Hao,WEI Guangjun,et al.Dam?age simulation of target metal plate by fragments and shock waves[J].Fire Control&Command Control,2014,39(7):152-154.

    [7] 郭淼,袁俊明,劉玉存,等.破片和沖擊波對(duì)相控陣?yán)走_(dá)天線罩的復(fù)合毀傷研究[J].爆破,2014,31(1):114-117,149. GUO Miao,YUAN Junming,LIU Yucun,et al.Simu?lation of composite damage for phased array radome by fragments and shock wave[J].Blasting,2014,31(1):114-117,149.

    [8] 劉剛,李向東,張媛.破片和沖擊波對(duì)直升機(jī)旋翼聯(lián)合毀傷仿真研究[J].計(jì)算機(jī)仿真,2013,30(6):68-71,181. LIU Gang,LI Xiangdong,ZHANG Yuan.Numerical simulation of combined damage of fragments and shock wave on helicopter rotor[J].Computer Simulation,2013,30(6):68-71,181.

    [9] 楊曙光.常規(guī)武器爆炸破片與沖擊波復(fù)合破壞效應(yīng)數(shù)值模擬研究[D].蘭州:蘭州大學(xué),2008.

    [10] 呂勇,石全,錢芳,等.不同時(shí)序破片和沖擊波在對(duì)天線的復(fù)合毀傷分析[J].火力與指揮控制,2014,39(3):128-132,139. LV Yong,SHI Quan,QIAN Fang,et al.Analysis of multiplex damages for fragment and shock wave in dif?ferent time sequence to antenna[J].Fire Control& Command Control,2014,39(3):128-132,139.

    [11] HATCH-AGUILAR T,NAJJAR F,SZYMANSKI E. Computational hydrocode study of target damage due to fragment-blast impact[C]//26th International Sym?posium of Ballistics.Miami,US,2011.

    [12] HALLQUIST J O.LS-DYNA theory manual[M].Cal?ifornia:Livermore Software Technology Corporation,2006.

    [13] CHUNG KIM YUEN S,LANGDON G S,NURICK G N,et al.Response of V-shape plates to localised blast load:experiments and numerical simulation[J].Inter?national Journal of Impact Engineering,2012,46:97-109.

    [14] CHENG D S,HUNG C W,PI S J.Numerical simula?tion of near-field explosion[J].Journal of Applied Science and Engineering,2013,16(1):61-67.

    [15] ZHANG P,CHENG Y S,LIU J,et al.Experimental and numerical investigations on laser-welded corru?gated-core sandwich panels subjected to air blast loading[J].Marine Structures,2015,40:225-246.

    [16] 梅志遠(yuǎn),朱錫.利用MSC/DYTRAN程序仿真分析導(dǎo)彈戰(zhàn)斗部立方體破片的侵徹威力[J].海軍工程大學(xué)學(xué)報(bào),2002,14(2):39-42. MEI Zhiyuan,ZHU Xi.The resistance ability of steel plate for penetration by cubic fragment through MSC/ Dytran[J].Journal of Naval University of Engineer?ing,2002,14(2):39-42.

    [責(zé)任編輯:盧圣芳]

    Numerical analysis of the damage on I-core sandwich panels subjected to combined blast and fragment loading

    DUAN Xinfeng,CHENG Yuansheng,ZHANG Pan,LIU Jun,LI Yong
    School of Naval Architecture and Ocean Engineering,Huazhong University of Science and Technology,Wuhan 430074,China

    In this paper,the combined blast and fragment loads due to close-in warhead explosion is simu?lated by placing prefabricated fragments at the bottom of bare TNT explosive,and the corresponding dam?age on I-core sandwich panels are evaluated using the software LS-DYNA.The differences of failure modes of sandwich panels subjected to only blast and combined blast and fragment loads are analyzed,and the effects of core configuration and thickness configuration between the front and back plates on the fail?ure modes are studied.The anti-damage performance of sandwich panels is then compared with that of the equivalent solid plate.In addition,the energy absorption characteristics of I-core sandwich panels under different loads are also analyzed.Numerical results show that the overall damage of both the sandwich pan?els and solid plates subjected to combined blast and fragment loads is more severe than that caused by blast only.The performance of I-core sandwich panels is superior to equivalent solid plates when the load?ing intensity is small.It is also found out that the loading intensity,loading type(blast loading or combined blast and fragment loading)and the thickness configuration between the front and back plates yield signifi?cant influences on the failure modes of the I-core sandwich panel.From the results of energy analysis,it can be concluded that the front face plate and the core are the main energy-absorbing parts under blast loading,while the front and back face plates become the main energy-absorbing parts when combined blast and fragment loading is considered.

    sandwich structure;blast;fragment;LS-DYNA;failure mode

    U661.43

    A

    10.3969/j.issn.1673-3185.2015.06.008

    http://www.cnki.net/kcms/detail/42.1755.TJ.20151110.1026.022.html期刊網(wǎng)址:www.ship-research.com

    段新峰,程遠(yuǎn)勝,張攀,等.沖擊波和破片聯(lián)合作用下I型夾層板毀傷仿真[J].中國(guó)艦船研究,2015,10(6):45-59. DUAN Xinfeng,CHENG Yuansheng,ZHANG Pan,et al.Numerical analysis of the damage on I-core sandwich panels subjected to combined blast and fragment loading[J].Chinese Journal of Ship Research,2015,10(6):45-59.

    2015-03-17 < class="emphasis_bold"> 網(wǎng)絡(luò)出版時(shí)間:

    時(shí)間:2015-11-10 10:26

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51209099)

    段新峰,男,1988年生,碩士生。研究方向:爆炸沖擊動(dòng)力學(xué)。E-mail:xfduan@hust.edu.cn程遠(yuǎn)勝(通信作者),男,1962年生,博士,教授,博士生導(dǎo)師。研究方向:結(jié)構(gòu)沖擊動(dòng)力學(xué)與防護(hù)設(shè)計(jì),結(jié)構(gòu)分析與優(yōu)化,結(jié)構(gòu)振動(dòng)與噪聲控制。E-mail:yscheng@hust.edu.cn

    猜你喜歡
    變形
    變形記
    談詩(shī)的變形
    柯西不等式的變形及應(yīng)用
    “變形記”教你變形
    不會(huì)變形的云
    “我”的變形計(jì)
    會(huì)變形的折紙
    童話世界(2018年14期)2018-05-29 00:48:08
    變形巧算
    例談拼圖與整式變形
    會(huì)變形的餅
    免费在线观看影片大全网站| 亚洲欧美色中文字幕在线| 韩国精品一区二区三区| 十八禁网站免费在线| 亚洲国产欧美一区二区综合| 美女扒开内裤让男人捅视频| 欧美日韩av久久| 无遮挡黄片免费观看| 岛国毛片在线播放| 黄片大片在线免费观看| 淫妇啪啪啪对白视频 | 男女无遮挡免费网站观看| 亚洲av日韩精品久久久久久密| 国产黄频视频在线观看| 欧美精品av麻豆av| 久久久精品区二区三区| 99国产精品免费福利视频| 这个男人来自地球电影免费观看| 亚洲av电影在线观看一区二区三区| 黄色视频在线播放观看不卡| 久热这里只有精品99| 满18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利乱码中文字幕| 最近最新中文字幕大全免费视频| 欧美日韩亚洲国产一区二区在线观看 | cao死你这个sao货| 91麻豆av在线| a 毛片基地| 国产xxxxx性猛交| 香蕉丝袜av| 国产一区二区在线观看av| 窝窝影院91人妻| 999久久久精品免费观看国产| 午夜福利一区二区在线看| 满18在线观看网站| 国产高清videossex| 老司机在亚洲福利影院| 考比视频在线观看| 国产色视频综合| 成人影院久久| 国产精品国产三级国产专区5o| 两个人免费观看高清视频| 天天操日日干夜夜撸| 色综合欧美亚洲国产小说| 欧美 亚洲 国产 日韩一| 日韩欧美国产一区二区入口| 在线天堂中文资源库| 国产一区二区三区av在线| 制服人妻中文乱码| 亚洲国产欧美日韩在线播放| 午夜福利免费观看在线| 在线av久久热| 成年人午夜在线观看视频| 巨乳人妻的诱惑在线观看| 肉色欧美久久久久久久蜜桃| 2018国产大陆天天弄谢| 丰满饥渴人妻一区二区三| 亚洲成人国产一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻人人添人人爽欧美一区卜| 91字幕亚洲| av有码第一页| 成人国产av品久久久| 国产老妇伦熟女老妇高清| 午夜福利视频在线观看免费| 69av精品久久久久久 | 丝袜脚勾引网站| 日韩精品免费视频一区二区三区| 亚洲精品日韩在线中文字幕| 亚洲精品久久午夜乱码| av在线app专区| 99国产综合亚洲精品| 99久久精品国产亚洲精品| 丰满少妇做爰视频| 日韩精品免费视频一区二区三区| av超薄肉色丝袜交足视频| 多毛熟女@视频| 亚洲欧美精品综合一区二区三区| 国产又色又爽无遮挡免| 国产精品成人在线| 亚洲av片天天在线观看| av在线老鸭窝| 男女边摸边吃奶| 欧美大码av| 久久久久久免费高清国产稀缺| 亚洲国产日韩一区二区| 国产成人av教育| 久久精品国产综合久久久| 一个人免费看片子| 免费看十八禁软件| 窝窝影院91人妻| 99久久人妻综合| 国内毛片毛片毛片毛片毛片| 欧美另类一区| 天堂俺去俺来也www色官网| 美女高潮喷水抽搐中文字幕| 国产麻豆69| 人妻人人澡人人爽人人| 久久久精品免费免费高清| 久久久久国内视频| 少妇猛男粗大的猛烈进出视频| 日韩制服丝袜自拍偷拍| 蜜桃在线观看..| 叶爱在线成人免费视频播放| 高清黄色对白视频在线免费看| 又大又爽又粗| 亚洲一卡2卡3卡4卡5卡精品中文| 各种免费的搞黄视频| 搡老乐熟女国产| 国产免费视频播放在线视频| 亚洲午夜精品一区,二区,三区| 悠悠久久av| 亚洲国产看品久久| 亚洲中文av在线| a级毛片在线看网站| 精品久久久久久电影网| √禁漫天堂资源中文www| 久久 成人 亚洲| 在线 av 中文字幕| 日本五十路高清| 捣出白浆h1v1| 97精品久久久久久久久久精品| 在线观看www视频免费| 午夜日韩欧美国产| 日韩免费高清中文字幕av| a 毛片基地| 91字幕亚洲| 日韩熟女老妇一区二区性免费视频| 国产熟女午夜一区二区三区| 法律面前人人平等表现在哪些方面 | 免费在线观看日本一区| 妹子高潮喷水视频| 在线精品无人区一区二区三| 在线观看免费日韩欧美大片| 高清在线国产一区| 欧美日韩av久久| 在线观看免费视频网站a站| 狠狠狠狠99中文字幕| 国产成人av教育| 真人做人爱边吃奶动态| 午夜两性在线视频| 1024香蕉在线观看| 交换朋友夫妻互换小说| 国产精品亚洲av一区麻豆| 99九九在线精品视频| 在线亚洲精品国产二区图片欧美| 久久国产精品男人的天堂亚洲| 色94色欧美一区二区| 中文字幕人妻熟女乱码| 午夜激情久久久久久久| 在线观看舔阴道视频| 女人高潮潮喷娇喘18禁视频| 老司机影院毛片| 亚洲av男天堂| 国产精品久久久久成人av| videos熟女内射| 久久中文字幕一级| 精品欧美一区二区三区在线| 色婷婷久久久亚洲欧美| 国产又爽黄色视频| 国产一区二区三区在线臀色熟女 | 亚洲 欧美一区二区三区| 久久久久久免费高清国产稀缺| 亚洲久久久国产精品| 制服诱惑二区| 黑人猛操日本美女一级片| 欧美国产精品va在线观看不卡| 精品人妻1区二区| 国产一区二区三区在线臀色熟女 | 国产精品 欧美亚洲| 丝袜美足系列| 亚洲精品国产区一区二| 一二三四社区在线视频社区8| 两性夫妻黄色片| 亚洲精品国产色婷婷电影| 亚洲欧美一区二区三区久久| 动漫黄色视频在线观看| 搡老岳熟女国产| 99精国产麻豆久久婷婷| 欧美成狂野欧美在线观看| 91大片在线观看| 欧美日韩av久久| 80岁老熟妇乱子伦牲交| 一级片免费观看大全| 各种免费的搞黄视频| 少妇猛男粗大的猛烈进出视频| 嫩草影视91久久| 19禁男女啪啪无遮挡网站| 青草久久国产| 国产男人的电影天堂91| 大片免费播放器 马上看| 国产欧美日韩一区二区三区在线| 亚洲精品一二三| 嫁个100分男人电影在线观看| 黄片播放在线免费| 纯流量卡能插随身wifi吗| 黑丝袜美女国产一区| 2018国产大陆天天弄谢| 亚洲熟女毛片儿| 中文字幕人妻丝袜一区二区| 精品熟女少妇八av免费久了| 在线观看免费高清a一片| 成年av动漫网址| 十八禁网站免费在线| 欧美精品一区二区免费开放| 高清在线国产一区| 色老头精品视频在线观看| 少妇精品久久久久久久| 国产精品.久久久| 久久香蕉激情| 99久久99久久久精品蜜桃| 91成人精品电影| 亚洲精品自拍成人| 国产高清videossex| 亚洲欧美精品综合一区二区三区| 精品乱码久久久久久99久播| 久久久国产成人免费| 亚洲精品国产av成人精品| 男人添女人高潮全过程视频| 欧美在线一区亚洲| av片东京热男人的天堂| 欧美精品啪啪一区二区三区 | 亚洲国产日韩一区二区| 亚洲精华国产精华精| 两人在一起打扑克的视频| 久久影院123| 又紧又爽又黄一区二区| 久久久欧美国产精品| 亚洲精品久久久久久婷婷小说| 亚洲精品乱久久久久久| 真人做人爱边吃奶动态| 不卡一级毛片| 欧美变态另类bdsm刘玥| 欧美日韩成人在线一区二区| 亚洲精品久久午夜乱码| 丰满饥渴人妻一区二区三| 国产精品麻豆人妻色哟哟久久| 欧美日韩亚洲综合一区二区三区_| 成人国语在线视频| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线| 国产一区二区三区综合在线观看| 啪啪无遮挡十八禁网站| 欧美少妇被猛烈插入视频| 精品少妇内射三级| 婷婷色av中文字幕| 捣出白浆h1v1| 99精品久久久久人妻精品| 亚洲欧洲精品一区二区精品久久久| 免费在线观看视频国产中文字幕亚洲 | 咕卡用的链子| 免费日韩欧美在线观看| 国产视频一区二区在线看| 亚洲三区欧美一区| 人人妻人人爽人人添夜夜欢视频| 精品国产乱码久久久久久小说| 亚洲自偷自拍图片 自拍| 热99久久久久精品小说推荐| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 热re99久久精品国产66热6| 国产精品国产av在线观看| 国产福利在线免费观看视频| 亚洲成国产人片在线观看| 精品卡一卡二卡四卡免费| 91成人精品电影| 91av网站免费观看| 欧美黄色淫秽网站| 不卡av一区二区三区| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 男女高潮啪啪啪动态图| 黑人巨大精品欧美一区二区mp4| 精品乱码久久久久久99久播| 多毛熟女@视频| 午夜福利视频在线观看免费| 日本vs欧美在线观看视频| 国产一级毛片在线| 欧美成狂野欧美在线观看| 午夜激情av网站| 日韩欧美免费精品| 亚洲av成人一区二区三| 成人亚洲精品一区在线观看| 久久热在线av| 午夜影院在线不卡| 这个男人来自地球电影免费观看| 久久久久久久精品精品| 亚洲成人免费av在线播放| 亚洲国产欧美一区二区综合| 黄色a级毛片大全视频| videos熟女内射| 国产成人精品无人区| 国产片内射在线| 亚洲欧美激情在线| 日本黄色日本黄色录像| 一本—道久久a久久精品蜜桃钙片| 夜夜骑夜夜射夜夜干| 成年av动漫网址| 香蕉丝袜av| 欧美另类亚洲清纯唯美| 国产日韩欧美视频二区| 国产欧美日韩精品亚洲av| 免费在线观看黄色视频的| 国产亚洲精品一区二区www | av片东京热男人的天堂| 欧美性长视频在线观看| 国产真人三级小视频在线观看| 国产成人影院久久av| 9色porny在线观看| 国产欧美日韩一区二区三区在线| 精品免费久久久久久久清纯 | 另类亚洲欧美激情| 国产成人影院久久av| 高潮久久久久久久久久久不卡| 成人国产一区最新在线观看| 天天添夜夜摸| 国产亚洲欧美精品永久| 建设人人有责人人尽责人人享有的| 亚洲中文av在线| 大片电影免费在线观看免费| 水蜜桃什么品种好| 免费在线观看完整版高清| 午夜精品久久久久久毛片777| 美女福利国产在线| 欧美av亚洲av综合av国产av| 日韩人妻精品一区2区三区| 9色porny在线观看| 新久久久久国产一级毛片| 亚洲欧美一区二区三区黑人| 亚洲精品国产色婷婷电影| 少妇粗大呻吟视频| 成在线人永久免费视频| 夫妻午夜视频| 老汉色∧v一级毛片| 51午夜福利影视在线观看| 中文字幕最新亚洲高清| 日韩制服丝袜自拍偷拍| 肉色欧美久久久久久久蜜桃| 黑人猛操日本美女一级片| 建设人人有责人人尽责人人享有的| 国产精品一区二区在线观看99| √禁漫天堂资源中文www| 悠悠久久av| 国产一区二区三区在线臀色熟女 | 久久热在线av| 好男人电影高清在线观看| 91av网站免费观看| 国产片内射在线| 午夜免费成人在线视频| 纵有疾风起免费观看全集完整版| 精品一区二区三区av网在线观看 | 欧美激情久久久久久爽电影 | 婷婷色av中文字幕| 黄色视频,在线免费观看| 婷婷丁香在线五月| 777久久人妻少妇嫩草av网站| 欧美激情极品国产一区二区三区| 亚洲伊人久久精品综合| 丝袜脚勾引网站| 啦啦啦 在线观看视频| 久久久久国产一级毛片高清牌| 一区二区三区激情视频| 国产成人精品在线电影| 亚洲avbb在线观看| 丰满饥渴人妻一区二区三| 大香蕉久久网| 91麻豆av在线| 动漫黄色视频在线观看| 精品久久久久久久毛片微露脸 | av不卡在线播放| 精品人妻在线不人妻| 中文欧美无线码| 精品国产国语对白av| 91麻豆精品激情在线观看国产 | 中文字幕色久视频| 亚洲中文字幕日韩| 国产免费现黄频在线看| 夫妻午夜视频| 黄色 视频免费看| 国产伦人伦偷精品视频| 黑人欧美特级aaaaaa片| 啦啦啦啦在线视频资源| 91精品三级在线观看| 亚洲欧美成人综合另类久久久| 9色porny在线观看| 老司机靠b影院| 精品国产乱码久久久久久男人| 国产精品欧美亚洲77777| 久久久久久人人人人人| 肉色欧美久久久久久久蜜桃| 高潮久久久久久久久久久不卡| 一个人免费看片子| 别揉我奶头~嗯~啊~动态视频 | 国产三级黄色录像| 国产男女内射视频| 午夜精品国产一区二区电影| 电影成人av| 国产免费视频播放在线视频| 青春草视频在线免费观看| 午夜精品久久久久久毛片777| 国产精品久久久久成人av| 91精品三级在线观看| 可以免费在线观看a视频的电影网站| 动漫黄色视频在线观看| 狠狠婷婷综合久久久久久88av| 一级毛片女人18水好多| 免费少妇av软件| 精品国产乱码久久久久久小说| 9热在线视频观看99| 亚洲av电影在线进入| 久久影院123| 国产真人三级小视频在线观看| 男人操女人黄网站| 97人妻天天添夜夜摸| 五月开心婷婷网| 天堂8中文在线网| 国产在线观看jvid| 亚洲美女黄色视频免费看| 欧美激情久久久久久爽电影 | 777久久人妻少妇嫩草av网站| 亚洲av日韩精品久久久久久密| 亚洲av电影在线观看一区二区三区| 在线观看人妻少妇| 国产熟女午夜一区二区三区| 日韩 亚洲 欧美在线| 中文字幕另类日韩欧美亚洲嫩草| 国精品久久久久久国模美| 18禁国产床啪视频网站| 日韩大片免费观看网站| 悠悠久久av| 1024视频免费在线观看| 亚洲专区字幕在线| 精品少妇黑人巨大在线播放| 国产激情久久老熟女| tube8黄色片| 精品熟女少妇八av免费久了| 亚洲专区国产一区二区| 国产有黄有色有爽视频| 午夜福利视频精品| 中文字幕人妻丝袜一区二区| 大片电影免费在线观看免费| 嫩草影视91久久| 国产av又大| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲 | 精品人妻熟女毛片av久久网站| 中文字幕色久视频| 亚洲精品中文字幕在线视频| 亚洲av日韩在线播放| 国精品久久久久久国模美| 日韩欧美一区二区三区在线观看 | 一区二区三区激情视频| 国产精品久久久av美女十八| 国产熟女午夜一区二区三区| 亚洲精品一卡2卡三卡4卡5卡 | 免费在线观看黄色视频的| 一级a爱视频在线免费观看| 在线看a的网站| 久久 成人 亚洲| 亚洲,欧美精品.| av网站在线播放免费| 亚洲av成人不卡在线观看播放网 | 久久国产精品影院| 欧美日韩亚洲综合一区二区三区_| 黄色毛片三级朝国网站| 欧美老熟妇乱子伦牲交| 亚洲国产中文字幕在线视频| 午夜精品国产一区二区电影| 麻豆av在线久日| 搡老熟女国产l中国老女人| 中文字幕人妻丝袜一区二区| 国产区一区二久久| 飞空精品影院首页| 免费人妻精品一区二区三区视频| 曰老女人黄片| 国产有黄有色有爽视频| 午夜日韩欧美国产| 捣出白浆h1v1| 黑人欧美特级aaaaaa片| 国产不卡av网站在线观看| 精品欧美一区二区三区在线| 亚洲 国产 在线| 国产欧美日韩一区二区精品| 不卡一级毛片| 桃红色精品国产亚洲av| 久久久久国产精品人妻一区二区| 精品国产乱码久久久久久男人| 青青草视频在线视频观看| avwww免费| 日韩中文字幕视频在线看片| 乱人伦中国视频| 国产欧美日韩一区二区精品| 亚洲五月色婷婷综合| 免费观看av网站的网址| 一级毛片女人18水好多| 欧美午夜高清在线| 人妻久久中文字幕网| 国产日韩欧美亚洲二区| 最黄视频免费看| 又紧又爽又黄一区二区| 纵有疾风起免费观看全集完整版| 精品久久久精品久久久| 成年人免费黄色播放视频| av天堂久久9| 少妇精品久久久久久久| 日韩大片免费观看网站| 天天躁夜夜躁狠狠躁躁| 波多野结衣一区麻豆| 黄色 视频免费看| 91国产中文字幕| 999精品在线视频| 一级片免费观看大全| 亚洲av成人一区二区三| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲综合色网址| 国产精品影院久久| 美女主播在线视频| 曰老女人黄片| 另类亚洲欧美激情| 香蕉国产在线看| 99九九在线精品视频| 久久久精品免费免费高清| 高清av免费在线| 国产真人三级小视频在线观看| 777米奇影视久久| av在线老鸭窝| 少妇的丰满在线观看| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 黑人巨大精品欧美一区二区蜜桃| 波多野结衣一区麻豆| 日韩中文字幕欧美一区二区| 电影成人av| 手机成人av网站| 在线看a的网站| 免费在线观看视频国产中文字幕亚洲 | 免费看十八禁软件| 亚洲欧美激情在线| 欧美日韩成人在线一区二区| 动漫黄色视频在线观看| 如日韩欧美国产精品一区二区三区| 婷婷成人精品国产| 精品熟女少妇八av免费久了| 999久久久国产精品视频| 一区二区av电影网| 男女午夜视频在线观看| av又黄又爽大尺度在线免费看| 男女无遮挡免费网站观看| 伦理电影免费视频| 亚洲,欧美精品.| 国产av国产精品国产| 天天添夜夜摸| 国产精品熟女久久久久浪| 国产精品 国内视频| 日韩中文字幕欧美一区二区| 一区在线观看完整版| 成人国产一区最新在线观看| 成人国产av品久久久| 嫩草影视91久久| 99热网站在线观看| 热99久久久久精品小说推荐| 国产亚洲精品一区二区www | 亚洲精品自拍成人| 国产野战对白在线观看| 精品国产乱子伦一区二区三区 | 黄色 视频免费看| 高清av免费在线| 国产免费现黄频在线看| 午夜两性在线视频| 亚洲第一欧美日韩一区二区三区 | 肉色欧美久久久久久久蜜桃| 亚洲国产欧美日韩在线播放| 午夜成年电影在线免费观看| 黑丝袜美女国产一区| 动漫黄色视频在线观看| 人妻久久中文字幕网| 岛国毛片在线播放| 精品少妇一区二区三区视频日本电影| 欧美黄色淫秽网站| 在线天堂中文资源库| av免费在线观看网站| 国产成人免费观看mmmm| 精品人妻熟女毛片av久久网站| 欧美日韩av久久| 99久久精品国产亚洲精品| 国产精品二区激情视频| 国产精品久久久久久精品电影小说| 午夜老司机福利片| 午夜福利,免费看| 亚洲一区中文字幕在线| 国产极品粉嫩免费观看在线| 91麻豆精品激情在线观看国产 | 亚洲色图 男人天堂 中文字幕| 日韩精品免费视频一区二区三区| 亚洲av成人一区二区三| 香蕉国产在线看| 久久亚洲精品不卡| 人妻人人澡人人爽人人| 午夜精品久久久久久毛片777| h视频一区二区三区| 王馨瑶露胸无遮挡在线观看| 日本av免费视频播放| 精品一品国产午夜福利视频| 王馨瑶露胸无遮挡在线观看| 日本av免费视频播放| 一区二区日韩欧美中文字幕| 久久精品国产综合久久久| 亚洲天堂av无毛| 在线永久观看黄色视频| 我要看黄色一级片免费的| 国产欧美日韩综合在线一区二区| 18禁国产床啪视频网站| 不卡av一区二区三区| 久久这里只有精品19|