纖維增強復合材料三明治板破片穿甲數(shù)值仿真

徐豫新， 戴文喜， 王樹山, 趙曉旭

(1.北京理工大學 爆炸科學與技術(shù)國家重點實驗室，北京 100081；2. 中國艦船研究設(shè)計中心，武漢 430064)

夾層復合材料由性質(zhì)不同的表面材料與芯材組合而成[1]，不僅能克服單一材料缺陷，且具有高硬度、高強度、高韌性、低密度、低成本等特性，廣泛應用于航空航天、兵器工業(yè)、交通運輸及工程防護等領(lǐng)域。對三明治板作為夾層復合材料抗沖擊載荷已開展實驗研究，若干種三明治板的吸能機制及防護性能亦有報道[2-8]。但該研究集中于金屬泡沫或周期性點陣三明治板在爆炸載荷或中低速(0～500 m/s)[9]彈體沖擊下局部及整體凹陷變形及損傷規(guī)律[2-8，10-13]，加載、材料響應時間均在毫秒量級，研究內(nèi)容屬于低應變率的結(jié)構(gòu)力學行為問題。破片高速(500～2 000 m/s)[9]沖擊下，加載及材料響應時間為微秒量級，彈靶界面瞬間形成高溫、高壓、高應變率區(qū)域，慣性效應及應變率效應在局部顯著，輕質(zhì)金屬泡沫夾層材料及周期性點陣結(jié)構(gòu)的相對強度及剛度下降，吸能比率隨撞擊速度的提高不斷降低[7]，雖可通過合理的材料尺寸配比后與前后面板復合作用，提高結(jié)構(gòu)整體防護效果，但因結(jié)構(gòu)較單板復雜，設(shè)計難度增加，較難用于艦體防護艙壁。纖維增強復合材料比強度、比模量均高于金屬材料，且輕質(zhì)、吸能高，與鋼板有機組合的纖維增強復合材料三明治板可用于艦船內(nèi)艙體防護結(jié)構(gòu)。

對由鋼板、纖維增強復合材料板及鋼板疊合而成的纖維增強復合材料三明治板，在結(jié)構(gòu)設(shè)計中需據(jù)防護位置的重要程度、纖維增強復合材料夾層板厚度實現(xiàn)全船結(jié)構(gòu)優(yōu)化。因此出現(xiàn)在前、后鋼板厚度確定條件下，夾層板厚度與不同質(zhì)量破片貫徹結(jié)構(gòu)所需最小速度是否具有相關(guān)性問題。由于纖維增強復合材料夾層板具有的各向異性力學特性會增加分析其在破片高速撞擊下?lián)p傷模式及自身力學性能對結(jié)構(gòu)整體防護性能影響難度。故多以實驗研究[14]為主。國外基于有限元的數(shù)值仿真[15-16]亦不夠系統(tǒng)。而國內(nèi)研究尚未見報道。

本文采用有限差分分析程序AutoDyn進行破片模擬彈丸(FSP)對纖維增強復合材料三明治板高速穿甲數(shù)值仿真研究，通過材料模型及參數(shù)合理選擇，獲得FSP彈丸[17-19]對纖維增強復合材料三明治板的彈道極限。將數(shù)值仿真結(jié)果與實驗結(jié)果[14]對比，驗證數(shù)值仿真的可信度，并進一步研究破片臨界貫穿條件下纖維增強復合材料三明治板各組成部分的吸能比率與結(jié)構(gòu)尺寸相關(guān)性。建立破片貫徹纖維增強復合材料三明治板彈道極限速度計算方法。

1 數(shù)值仿真

1.1 仿真控制方程

AutoDyn程序以質(zhì)量守恒、動量守恒、能量守恒為基礎(chǔ)進行計算。數(shù)值仿真遵循的控制方程為：

(1) 質(zhì)量守恒

(1)

(2) 動量守恒

(2)

(3) 能量守恒

(3)

式中：ρ為材料密度；ui為速度；fi為單位質(zhì)量外力；σij為應力張量；e為總和比能，即比動能與比內(nèi)能之和：

(4)

1.2 仿真模型

1.2.1 幾何模型及離散化

據(jù)實驗研究[14]采用FSP彈丸及靶體結(jié)構(gòu)尺寸，見圖1。選cm-μs-g-Mbar單位制與1/2面對稱建模方法建立數(shù)值仿真模型。模型中靶體長(寬為120 mm×120 mm，結(jié)構(gòu)各層厚度由具體(16種)工況確定，破片與靶體破壞區(qū)域分別用0.5 mm、1.0 mm網(wǎng)格尺寸。通過TrueGrid軟件建立數(shù)值仿真幾何模型離散化后導入AutoDyn程序中，靶體四周施加固定約束，破片施加初始速度后進行仿真計算。數(shù)值仿真模型見圖2。

圖1 破片模擬彈丸(FSP)結(jié)構(gòu)尺寸(單位：mm)

圖2 數(shù)值仿真模型

1.2.2 材料模型及參數(shù)

文獻[14]中破片采用35CrMnSi鋼，質(zhì)量10 g，抗拉強度1 620 MPa(實測)，屈服強度1 275 MPa(實測)；三明治板面板、背板材料均為Q235A鋼，其屈服極限305 MPa(實測)，拉伸極限420 MPa(實測)，延伸率33%(實測)；夾層材料為芳綸、玻璃纖維板兩類。

由于破片材料35CrMnSi鋼為各向同性材料，在彈性段、塑性段應變隨應力的提高呈線性增加。因此，據(jù)AutoDyn材料手冊，用Von Mises模型描述材料強度，用Linear模型描述材料狀態(tài)變化，用Plastic Strain值判定材料失效與侵蝕[20]，并據(jù)文獻[21]研究結(jié)果獲得模型參數(shù)設(shè)置，見表1。對前、后面板材料Q235A鋼，用Johnson-Cook[22]模型描述本構(gòu)關(guān)系， 用Mie-Grüneisen[23]模型描述狀態(tài)變化，據(jù)文獻[24]及實測強度所得參數(shù)見表2。復合材料本構(gòu)行為的合理描述為材料學、力學領(lǐng)域研究熱點之一[25-28]。通常采用兩種方法進行處理：①將復合材料作為均質(zhì)各向異性材料；②由細觀角度考慮復合材料非均勻性。前者已取得滿意結(jié)果；后者雖與實際接近，但因大量細觀參數(shù)缺失，仍難以用于宏觀斷裂定量分析。借鑒文獻[16]研究結(jié)果，用Puff模型描述纖維材料狀態(tài)，用Von Mises模型描述纖維材料強度。據(jù)文獻[14-16]中實測抗拉強度及復合材料參數(shù)設(shè)置，獲得kevlar-129纖維材料及E-Glass纖維材料模型參數(shù)，見表3。

表1 35CrMnSi鋼的材料模型及參數(shù)

表2 Q235A鋼材料模型及參數(shù)

1.2.3 接觸算法及參數(shù)設(shè)置

文獻[14]中前、后板及芯板緊密貼合通過板四角螺栓予以固定，板與板間無縫隙。數(shù)值仿真中在三層板和與破片間設(shè)置Largrange/Lagrange耦合接觸，結(jié)構(gòu)引入小間隙確定子區(qū)域是否相互作用。該間隙對每個相互作用面定義一個作用的發(fā)現(xiàn)區(qū)。一旦節(jié)點進入發(fā)現(xiàn)區(qū)將被推出。但若間隙尺寸小于最小相互作用面尺寸的1/10或大于其1/2，Autodyn將不進行計算。因此，可通過網(wǎng)格共節(jié)點方式將三層板進行界面耦合。在仿真中間隙尺寸(Gap Size)設(shè)為0.002 cm，工作單元尺寸(Work Unit Size)設(shè)為平均(Average)模式，接觸安全因子(Safety Factor)設(shè)為0.2。

表3 纖維增強復合板材料模型及參數(shù)

1.3 數(shù)值仿真結(jié)果及實驗校驗

對16種實驗工況[14]進行破片數(shù)值仿真，獲得典型穿甲過程見圖3。用兩射彈彈道極限法[29]，獲得FSP彈丸貫穿不同結(jié)構(gòu)三明治板彈道極限，見表4。由表4看出，16種實驗工況數(shù)值仿真與實驗結(jié)果誤差均在10%以內(nèi)，表明采用網(wǎng)格劃分、接觸設(shè)置及材料模型參數(shù)模擬所得結(jié)果具有可信度，亦為用數(shù)值仿真的進一步研究奠定基礎(chǔ)。

圖3 破片對纖維增強復合材料三明治板穿甲仿真

2 數(shù)值仿真結(jié)果分析

2.1 臨界貫穿條件下能量轉(zhuǎn)換

彈道極限指彈體以給定著角貫穿一定厚度指定靶板的最小著速，為50%穿透概率條件下的特征著速。通常總將彈道極限視為彈體對靶體的臨界貫穿速度。破片以臨界貫穿速度侵徹靶體過程中初始動能轉(zhuǎn)化為崩濺塞塊的動能、彈體靶體的熱能、靶體各部分及破片的內(nèi)能等。據(jù)文獻[7]對多孔金屬夾層板吸能特性的研究方法，設(shè)RF，RM，RR分別為前板、夾層板、后板的吸能比率(即內(nèi)能增加量與彈體初始動能之比)，則：

表4 FSP彈丸貫穿每種結(jié)構(gòu)三明治板彈道極限

(5)

據(jù)數(shù)值仿真結(jié)果，通過式(5)獲得前、后板厚度相同條件下(前后、面板分別厚4 mm，6 mm)，破片以彈道極限速度侵徹時，三明治板各部分吸能比率見圖4。由圖4看出，不同厚度夾層板的吸能比率(C值)基本恒定，通過線性回歸獲得芳綸及玻璃纖維板的吸能比率分別為10.41%，2.68%，見圖5，芳綸纖維夾層板的吸能比率為玻璃纖維夾層板的3.88倍。

圖4 三明治板各部分吸能率

通過仿真獲得前、后板厚度不同時夾層纖維板的吸能比率見表5。由表5看出，前、后板厚度不同，夾層板吸能比率亦不同，前、后板厚度和與夾層板厚度之比越大，夾層板吸能比率越低，即夾層板吸收能量越少，前后板吸收能量越多。

表5 前、后板不同厚度夾層纖維板吸能比率

圖6 三明治板整體吸能比率隨夾層板厚度的變化

2.2 最小臨界速度計算方法

由于纖維增強復合材料三明治板的吸能性能隨中間夾層厚度的增加呈二次函數(shù)規(guī)律遞增[14]，與單一結(jié)構(gòu)纖維增強復合材料抗彈體貫穿性質(zhì)[30]類似。仿真結(jié)果見圖8。并由圖8擬和獲得破片以彈道極限速度侵徹時夾層板內(nèi)能增加量隨夾層板厚度函數(shù)關(guān)系式為：

(6)

式中：T為夾層板厚度， mm；A，B,C為系數(shù)，前面板厚度4 mm，后面板厚度8 mm時，A，B,C的取值與擬和曲線相關(guān)系數(shù)R2見表6。

表6 系數(shù)取值與擬和曲線相關(guān)系數(shù)

由式(5)、(6)可得破片貫徹纖維增強復合材料三明治板所需最小速度計算式為：

(7)

式中：v50為破片對纖維增強復合材料三明治板彈道極限速度， m/s；m為破片質(zhì)量， kg；前面板厚度4 mm，后面板厚度6 mm時，A,B,C值見表6，Rm值由圖5獲得。

2.3 計算方法驗證

針對前面板厚度4 mm，后面板厚度6 mm條件，用式(7)計算得：

(1) 試驗工況夾層板厚度(9～17 mm間)的彈道極限速度見表7，并與試驗結(jié)果對比；

(2) 20 mm夾層板厚度三明治板的彈道極限速度見表8，并與數(shù)值仿真值對比。

表7 10.0 g FSP對試驗工況三明治板的彈道極限

表8 10.0 g FSP彈丸對夾層板厚度20 mm三明治板彈道極限

由表7、表8看出，式(7)所得計算結(jié)果與試驗值、仿真值誤差均在10%以內(nèi)，說明：① 上述分析方法具有一定可靠性；② 本文給出的A,B,C,Rm系數(shù)適用于前面板厚度4 mm，后面板厚度6 mm的纖維增強復合材料三明治板彈道極限計算。

3 結(jié) 論

通過合理選擇材料模型與參數(shù)，采用有限差分程序AutoDyn可實現(xiàn)破片對纖維增強復合材料三明治板穿甲的數(shù)值仿真結(jié)果具有可信度，對數(shù)值仿真結(jié)果分析，結(jié)論如下：

(1) 破片以彈道極限速度穿甲過程中，前、后板厚度確定時，不同厚度夾層板的吸能比率恒定；

(2) 破片以彈道極限速度侵徹時，纖維增強復合材料三明治板的夾層板內(nèi)能隨厚度增加呈二次函數(shù)遞增；

(3) 對彈體結(jié)構(gòu)及前、后板厚度確定的纖維增強復合材料三明治板，可據(jù)不同厚度夾層板的吸能比率及纖維增強復合材料三明治板的夾層板內(nèi)能隨厚度增加規(guī)律獲得破片對三明治板的彈道極限計算式，且具一定可靠性。

參 考 文 獻

[1]李順林，王興業(yè). 復合材料結(jié)構(gòu)設(shè)計基礎(chǔ)[M]. 武漢：武漢工業(yè)大學出版社，1993.

[2]倪長也, 金 峰, 盧天健, 等. 3種點陣金屬三明治板的抗侵徹性能模擬分析[J]. 力學學報，2010,42(6):1125-1136.

NI Chang-ye, JIN Feng, LU Tian-jian, et al. Penetration and perforation performance of three pyramidal lattice-cored sandwich plates numerical simulations[J]. Chinese Journal of Theoretical and Applied Mechanics, 2010,42(6):1125-1136.

[3]倪長也, 金 峰, 盧天健. 超輕金屬點陣三明治板結(jié)構(gòu)抗侵徹性能分析[J]. 兵工學報，2009,30 (S2):94-96.

NI Chang-ye, JIN Feng，LU Tian-jian. Numerical simulation of the ballistic performance of composite pyramidal lattice truss structure[J]. Acta Armamentarii，2009,30 (S2):94-96.

[4]Phoenix S L,Porwal P K. A new membrane model for the ballistic impact response and V50performance of multi-ply fibrous systems[J]. International Journal of Solids and Structures,2003,40(24):6723-6765.

[5]Hoo Fatt M S, Sirivolu D. A wave propagation model for the high velocity impact response of a composite sandwich panel[J]. International Journal of Impact Engineering，2010,37(2):117-130.

[6]Dean J, S-Fallah A, Brown P M,et al. Energy absorption during projectile perforation of lightweight sandwich panels with metallic fibre cores[J]. Composite Structures，2011，93(3):1089-1095.

[7]趙桂平，盧天健. 多孔金屬夾層板在沖擊載荷作用下的動態(tài)響應[J]. 力學學報,2008,40(2):194-206.

ZHAO Gui-ping, LU Tian-jian. Dynamic response of cellular metallic sandwich plates under impact loading[J]. Chinese Journal of Theoretical and Applied Mechanics, 2008,40(2):194-206.

[8]石少卿，劉仁輝，汪 敏. 鋼板-泡沫鋁-鋼板新型復合結(jié)構(gòu)降低爆炸沖擊波性能研究[J]. 振動與沖擊，2008,27(4): 143-146.

SHI Shao-qing，LIU Ren-hui，WANG Min. Shock wave reduction behavior of a new compound structure composed of a foam aluminum layer between two steel plates[J]. Journal of Vibration and Shock, 2008,27(4):143-146.

[9]王曉強，朱 錫，梅志遠. 纖維增強復合材料抗侵徹研究綜述[J]. 玻璃鋼/復合材料,2008,34(5):47-55.

WANG Xiao-qiang, ZHU Xi, MEI Zhi-yuan. The development of fiber-reinforced composites under ballistic impact[J]. Fiber Reinforced Plastics/Composites，2008，34(5):47-55.

[10]Shen J H, Lu G X, Wang Z H,et al. Experiments on curved sandwich panels under blast loading[J]. International Journal of Impact Engineering，2010,37(9):960-970.

[11]Gustin J, Joneson A, Mahinfalah M, et al. Low velocity impact of combination Kevlar/carbon fiber sandwich composites[J]. Composite Structures,2005,69(4):396-406.

[12]董永香，馮順山，李學林. 爆炸波在硬-軟-硬三明治介質(zhì)中傳播特性的數(shù)值分析[J]. 彈道學報,2007,19(1):59-63.

DONG Yong-xiang，F(xiàn)ENG Shun-shan，Li Xue-lin. Numerical analysis of propagation characteristics of explosive wave in the hard-soft-hard sandwich media[J]. Journal of Ballistics,2007,19(1):59-63.

[13]Jackson M,Shukla A. Performance of sandwich composites subjected to sequential impact and air blast loading[J].Composites: Part B,2011,42(2):155-166.

[14]徐豫新，王樹山，嚴文康，等.纖維增強復合材料三明治板的破片穿甲實驗[J]. 復合材料學報, 2012,29(3)：72-77.

XU Yu-xin, WANG Shu-shan, YAN Wen-kang, et al. Armor-piercing experiment on fragment against sandwich plate with fiber reinforced composite cores[J]. Acta Materiae Compositae Sinica, 2012,29(3):72-77.

[15]Mahfuz H,Zhu Y H,Haque A,et al. Investigation of high-velocity impact on integral armor using finite element method[J]. International Journal of Impact Engineering,2000,24(2):203-217.

[16]Preece D S, Berg V S, Risenmay M A. Bullet impact on steel and kevlar(/steel armor-experimental data and hydrocode modeling with eulerian and lagrangian method[C]. 22ed Int.Symp.on Ballistics，Vancouver，BC Canada，14th～18th Nov. 2005:971-979.

[17]Wambua P, Vangrimde B, Lomov S, et al.The response of natural fibre composites to ballistic impact by fragment simulating projectiles[J]. Composite Structures,2007, 77(2):232-240.

[18]Zukas J R. High velocity impact dynamics[M].John Wiley and Sons, Inc. 1990.

[19]Iremonger M J, Went A C. Ballistic impact of fibre composite armours by fragment-simulating projectiles[J]. Composites Part A ,1996,27(7):575-581.

[20]Autodyn_Chapter9_Material_Models [EB/OL]. http://wenku.baidu.com/view/381b136527d3240c8447ef6a.html. 2012.

[21]王 琳，王富恥，王 魯，等. 空心彈體侵徹金屬靶板的數(shù)值模擬和實驗研究[J]. 兵器材料科學與工程，2001，24(6):13-17.

WANG Lin, WANG Fu-chi, WANG Lu, et al. Numerical simulation and impact test study of hollow projectiles impacting steel plates[J]. Ordnance Material Science and Engineering,2001,24(6):13-17.

[22]Johnson G R, Cook W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperature[C]. Proceeding of the 7th International Symposium on Ballistics. Hague, Netherlands: International Ballistics Committee, 1983:541-547.

[23]Steinberg D J, Cochran S G, Guinan M W. A constitutive model for metals applicable at high strain rate[J]. J.Appl.Phys, 1980,51(3):1498-1503.

[24]陳 剛，陳小偉，陳忠富,等. A3鋼鈍頭彈撞擊45鋼板破壞模式的數(shù)值分析[J]. 爆炸與沖擊，2007,27(5): 390-397.

CHEN Gang, CHEN Xiao-wei, CHEN Zhong-fu, et al. Simulation of A3 steel blunt projectiles impacting 45 steel plates[J]. Explosion and Shock Waves, 2007,27(5): 390-397.

[25]Grujicic M, Sun Y P, Koudela K L. The effect of covalent functionalization of carbon nanotube reinforcements on the atomic-level mechanical properties of polyvinylestereposxy [J]. Applied Surface Science,2007, 253:3009-3021.

[26]Grujicic M, Bell W C, Thompson L L,et al. Ballistic-protection performance of carbon-nanotube-doped poly-vinyl-ester-epoxy matrix composite armor reinforces with E-glass fiber mats[J]. Materials Science&Engineering, 2008,479(1-2):10-22.

[27]Grujicic M, Glomski P S, He T, et al.Material modeling and ballistic-resistance analysis of armor-grade composites reinforced with high-performance fibers[J]. Journal of Materials Engineering and Performance , 2009,18(9):1169-1182.

[28]Gama B A,John W, Gillespie J. Finite element model of impact, damage evolution and penetration of thick-section composites[J]. International Journal of Impact Engineering, 2011,38(4):181-197.

[29]隋樹元，王樹山. 終點效應學[M]. 北京：國防工業(yè)出版社,2000.

[30]馬曉青，韓 峰. 高速碰撞動力學[M]. 北京：國防工業(yè)出版社,1998.

[31]梅志遠，朱 錫，劉燕紅，等. 纖維增強復合材料層合板彈道沖擊研究進展[J].力學進展,2003,33(3):375-385.

MEI Zhi-yuan, ZHU Xi, LIU Yan-hong, et al. The developments of fibre-composite laminates under ballistic impact[J]. Advances In Mechanics,2003,33(3):375-385.