米 鵬，梁 森，張義霞

(1.青島理工大學(xué) 機(jī)械工程學(xué)院，青島 266033;2.新南威爾士大學(xué) 工程及信息技術(shù)學(xué)院國(guó)防軍事學(xué)院校區(qū)，堪培拉 2600)

嵌入式共固化纖維增強(qiáng)樹(shù)脂基復(fù)合材料阻尼結(jié)構(gòu)是將三種不同性質(zhì)的材料(如碳纖維、樹(shù)脂和粘彈性阻尼材料)，通過(guò)物理或化學(xué)的方法復(fù)合而成的一種多相固體。從它的組成與結(jié)構(gòu)上分析，其中有一相在層內(nèi)基本上是連續(xù)的被稱為基體相(如樹(shù)脂)，而另一相是分散而被基體所包容的被稱為增強(qiáng)相(如碳纖維)，還有一相是各向同性的粘彈性阻尼材料。其基體相、增強(qiáng)相和粘彈性阻尼材料在性能上起協(xié)調(diào)作用，實(shí)現(xiàn)大幅度地提高復(fù)合材料構(gòu)件阻尼的目的，得到單一材料難以比擬的綜合力學(xué)性能［1］。由于這種事先阻尼處理結(jié)構(gòu)是將阻尼材料嵌在結(jié)構(gòu)內(nèi)部的，因此其阻尼材料具有不易脫落、抗老化等優(yōu)點(diǎn)，在高速列車、航空、航天等高科技領(lǐng)域有著廣闊的應(yīng)用前景［1－3］。這種新型的復(fù)合材料阻尼結(jié)構(gòu)一經(jīng)提出，很快成為研究的熱點(diǎn)［2－5］。目前對(duì)此的相關(guān)研究主要集中在粘彈性阻尼材料薄膜的制作工藝、共固化工藝以及在實(shí)驗(yàn)基礎(chǔ)上的嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)動(dòng)力學(xué)特性［4－5］，針對(duì)嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)沖擊性能的研究相對(duì)較少。而沖擊破壞是復(fù)合材料經(jīng)常出現(xiàn)的破壞形式之一，因此有必要對(duì)嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)的抗沖擊性能進(jìn)行深入研究。

本文采用LS-DYNA顯式動(dòng)力學(xué)分析軟件，在相關(guān)實(shí)驗(yàn)的基礎(chǔ)上對(duì)嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)的低速?zèng)_擊性能進(jìn)行數(shù)值模擬，并將數(shù)值模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較，探討新型復(fù)合材料阻尼板對(duì)沖擊能量和沖擊峰值載荷的吸收能力。

1 數(shù)值模擬過(guò)程

LS-DYNA顯式動(dòng)力學(xué)分析程序，能夠模擬各種復(fù)雜非線性問(wèn)題，適合求解各種二維、三維非線性結(jié)構(gòu)的高速碰撞、爆炸和金屬成型等非線性動(dòng)力學(xué)沖擊問(wèn)題。LS-DYNA分析過(guò)程主要包括:前處理、求解和后處理［6］。

1.1 前處理

整個(gè)幾何模型包括沖擊錘與被沖擊試件兩部分。本文前處理在ANSYS/LS-DYNA中完成，主要包括:選擇單元類型和算法、定義實(shí)常數(shù)、定義材料屬性、構(gòu)建實(shí)體模型、有限元網(wǎng)格劃分、創(chuàng)建PART和定義接觸類型等。

本文網(wǎng)格單元全部采用8節(jié)點(diǎn)六面體實(shí)體模型SOLID164單元，采用默認(rèn)Lagrangian算法［7］，該單元無(wú)需定義實(shí)常數(shù)。

為節(jié)省計(jì)算時(shí)間，在模型中采用剛體材料來(lái)模擬實(shí)驗(yàn)中的鋼質(zhì)沖擊錘，用關(guān)鍵字*MAT_RIGID定義［8］。由于建模過(guò)程中只取沖擊錘底部的沖擊部分，未包含傳感器以及相關(guān)固定裝置，因此對(duì)該部分的密度進(jìn)行了相應(yīng)的修正，以確保模型質(zhì)量和實(shí)際質(zhì)量相等，其參數(shù)見(jiàn)表1。

表1 沖擊錘材料模型參數(shù)Tab.1 Material model parameters of the impact hammer

嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)中的復(fù)合材料層采用Chang-Chang失效模型，通過(guò)關(guān)鍵字*MAT_COMPOSITE_DAMAGE定義。該模型由Chang等［9］提出。它包含三個(gè)失效準(zhǔn)則:基體開(kāi)裂失效、壓縮失效和纖維斷裂失效。確定上述失效準(zhǔn)則除需要三個(gè)方向的彈性模量、剪切模量和泊松比九個(gè)參數(shù)外，還需要確定其他五個(gè)參數(shù)，依次為:主方向拉伸強(qiáng)度XT，橫向拉伸強(qiáng)度YT，面內(nèi)剪切強(qiáng)度SC，橫向壓縮強(qiáng)度YC和非線性剪切應(yīng)力參數(shù)ALPH。其XT，YT，SC，YC大小是通過(guò)材料的強(qiáng)度測(cè)試獲得，ALPH則是由材料的剪切應(yīng)變來(lái)確定［10－11］，具體值見(jiàn)表 2。

表2 樹(shù)脂基纖維增強(qiáng)型復(fù)合材料材料模型參數(shù)Tab.2 Material model parameters of fiber reinforced resin matrix composite

嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)中的粘彈性阻尼材料采用不可壓縮Mooney-Rivlin橡膠模型來(lái)模擬。此模型需要輸入系數(shù)A，B以及泊松比υ，以確定其應(yīng)變能密度函數(shù):

式中:

I1，I2，I3分別為變形張量不變量［10－11］。粘彈性阻尼材料模型具體參數(shù)值見(jiàn)表3。

表3 粘彈性阻尼材料模型參數(shù)Tab.3 Viscoelastic material model parameters

根據(jù)實(shí)際幾何尺寸來(lái)建立實(shí)體模型，其中沖擊錘是由一個(gè)半徑和高度均為10 mm的圓柱體，和一個(gè)半徑為10 mm的半球體組成。被沖擊的圓形試件共2個(gè)，試件1為樹(shù)脂基纖維增強(qiáng)型復(fù)合材料，試件2為嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)，它們的幾何參數(shù)、鋪層順序、鋪層角度具體見(jiàn)表4，其結(jié)構(gòu)見(jiàn)圖1。實(shí)驗(yàn)時(shí)試件約束邊界圓周的直徑為81 mm，在建模過(guò)程中以試件邊界約束尺寸作為試件的實(shí)際尺寸，分別建立兩個(gè)試件低速?zèng)_擊的實(shí)體模型。由于模型幾何結(jié)構(gòu)的對(duì)稱性，因此采用1/2模型進(jìn)行計(jì)算，從而減少計(jì)算量。

表4 試件的鋪層順序和幾何參數(shù)Tab.4 Specimen stacking sequence and geometric parameters

為獲得較理想的模擬結(jié)果，應(yīng)盡量將實(shí)體模型劃分為較為規(guī)則的六面體單元。在模擬過(guò)程中發(fā)現(xiàn)單元尺寸過(guò)大或過(guò)小均會(huì)影響模擬結(jié)果:當(dāng)單元尺寸過(guò)大時(shí)，則單元失效時(shí)吸收能量過(guò)多;反之則吸收能量太少［8，12］。經(jīng)過(guò)反復(fù)比較，本文模擬中單元的尺寸為1.5 mm左右時(shí)，能獲得最佳模擬結(jié)果。沖擊錘底部半球體單元尺寸與試件單元尺寸相近。劃分網(wǎng)格后的模型見(jiàn)圖2。

圖1 試件結(jié)構(gòu)示意圖Fig.1 Schematic configuration of specimens

圖2 1/2有限元模型Fig.2 1/2 Finite element model

在ANSYS/LS-DYNA前處理器中，可根據(jù)每個(gè)實(shí)體材料模型的不同自動(dòng)生成相應(yīng)的PART，兩個(gè)PART之間需要定義兩者的接觸類型。沖擊錘與復(fù)合材料之間以及沖擊錘與粘彈性阻尼材料之間的接觸均設(shè)置為面與面之間侵蝕接觸，該接觸形式是通過(guò)關(guān)鍵字*CONTACT_ERODING_SURFACE_TO_SURFACE定義的，其中包括四個(gè)默認(rèn)的選項(xiàng)外，還要定義一個(gè)Option Card A選項(xiàng)。由于沖擊錘的材料和網(wǎng)格大小與試件的材料與網(wǎng)格大小均有差別，因此采用基于段(Segment)的接觸算法，將該選項(xiàng)的參數(shù)SOFT設(shè)為2;不定義載荷曲線，參數(shù) LCIDAB設(shè)為0;參數(shù) SOFSCL取為0.1，BSORT 設(shè)為15，其他參數(shù)設(shè)置采用默認(rèn)數(shù)值［7－8］。由于在相關(guān)的沖擊實(shí)驗(yàn)中，嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)沒(méi)有出現(xiàn)復(fù)合材料層與層之間或復(fù)合材料與粘彈性材料之間的脫層現(xiàn)象，因此數(shù)值模擬中不考慮脫層現(xiàn)象的影響，復(fù)合材料層間以及復(fù)合材料層與粘彈性阻尼材料層間均采用共節(jié)點(diǎn)方式連接［12］。

1.2 求解及后處理

根據(jù)試驗(yàn)情況，約束試件圓周上節(jié)點(diǎn)的所有自由度，在試件和沖擊錘的剖面上施加對(duì)稱約束，并給沖擊錘施加了沿軸線方向指向試件的初速度。在試驗(yàn)中，試件1和試件2的沖擊初速度分別為2.165 m/s和1.3615 m/s，模擬中初速度大小與實(shí)驗(yàn)相同。

整個(gè)沖擊過(guò)程求解時(shí)間設(shè)置為0.4×10－1s，輸出結(jié)果文件的時(shí)間間隔為0.1×10－3s。最后，前處理器輸出關(guān)鍵字文件(k文件)。對(duì)無(wú)法在前處理器中修改的文件中部分關(guān)鍵字進(jìn)行手動(dòng)修改，將修改完畢后的關(guān)鍵字文件導(dǎo)入LS－DYNA求解器進(jìn)行求解后，將輸出的d3plot結(jié)果文件導(dǎo)入后處理器，獲得每個(gè)試件沖擊過(guò)程中的相關(guān)數(shù)據(jù)。

2 結(jié)果分析

2.1 能量吸收分析

將求解器輸出的d3plot文件輸入LS-PREPOST后處理器，獲取沖擊錘的動(dòng)能變化曲線，速度變化曲線，同文獻(xiàn)［13］中實(shí)驗(yàn)獲得的曲線及數(shù)據(jù)相比較。

試件1和試件2的數(shù)值模擬曲線分別如圖3和圖4所示。

試件的動(dòng)能、速度變化量的模擬值與實(shí)驗(yàn)結(jié)果的比較見(jiàn)表5。

表5 動(dòng)能、速度變化量的模擬值與實(shí)驗(yàn)結(jié)果的比較Tab.5 Comparison simulated kinetic energy and velocity variation with experiment data

從圖3、圖4以及表5可以看出:在本文所模擬的試件1和試件2的沖擊過(guò)程中，沖擊錘的速度變化曲線和動(dòng)能變化曲線與實(shí)驗(yàn)測(cè)試曲線吻合較好，動(dòng)能變化量、速度變化量的模擬值也與實(shí)驗(yàn)結(jié)果較為接近。

圖3 試件1的模擬曲線Fig.3 Simulation curve of specimen 1

圖4 試件2的模擬曲線Fig.4 Simulation curve of specimen 2

在模擬的沖擊過(guò)程中，由于不存在其他形式的能量損失，沖擊錘損失的動(dòng)能等于試件吸收的能量。通過(guò)表5可知，試件1在沖擊過(guò)程中吸收的能量為5.96 J，試件2沖擊過(guò)程中吸收的能量為7.00 J。兩者相比，試件2在沖擊過(guò)程中吸收的能量比試件1增加了17.5%。

2.2 沖擊載荷峰值分析

在LS-PREPOST后處理器中查看沖擊錘受到的沖擊載荷的峰值，然后將其與文獻(xiàn)［13］實(shí)驗(yàn)中的沖擊載荷峰值相比較，結(jié)果如表6所示。

表6 模擬的沖擊力峰值與實(shí)驗(yàn)結(jié)果比較Tab.6 Comparison between simulation peak force and experiment data

從表6可看出:模擬沖擊力峰值與實(shí)驗(yàn)的測(cè)試值非常接近，本文數(shù)值模擬較準(zhǔn)確的反映了沖擊過(guò)程中沖擊力峰值的大小。

將表6中試件1和試件2的沖擊力峰值的模擬值進(jìn)行比較可知:在模擬的沖擊過(guò)程中，試件2受到的沖擊力峰值比試件1增加了約140%。

2.3 破壞過(guò)程分析

在數(shù)值模擬中，當(dāng)有限元模型中某些單元的應(yīng)力或應(yīng)變達(dá)到預(yù)先設(shè)定的失效準(zhǔn)則時(shí)，這些單元將會(huì)被從模型中刪除，以模擬材料的破壞或裂縫的產(chǎn)生過(guò)程。圖5、圖6分別將試件1和試件2的模擬破壞過(guò)程與實(shí)際破壞過(guò)程進(jìn)行比較。

由圖5、圖6可知，在數(shù)值模擬中，當(dāng)沖擊錘接觸試件時(shí)，在相互垂直的直徑方向上的單元，由內(nèi)而外依次出現(xiàn)變形。隨著沖擊錘不斷向下運(yùn)動(dòng)，單元的變形不斷增大，一旦達(dá)到預(yù)先設(shè)定的失效準(zhǔn)則，相應(yīng)的單元被刪除，從而出現(xiàn)裂縫。隨著沖擊過(guò)程的進(jìn)行，裂縫由內(nèi)而外不斷擴(kuò)展，直至沖擊過(guò)程結(jié)束。最終，被沖擊破壞的試件在中心處產(chǎn)生“十”狀的裂縫，這與實(shí)際的破壞情形一致。

比較圖5和圖6可知:在沖擊過(guò)程中，試件2的剛度小于試件1，其破壞部分的面積和產(chǎn)生的裂紋長(zhǎng)度遠(yuǎn)小于試件1;試件1呈脆性破壞，而試件2呈現(xiàn)塑性破壞的特點(diǎn)。上述特點(diǎn)與文獻(xiàn)［13］的結(jié)論一致。

3 結(jié)論

本文采用LS-DYNA顯式動(dòng)力學(xué)軟件建立了嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)的有限元分析模型，對(duì)這種新型阻尼結(jié)構(gòu)的低速?zèng)_擊性能進(jìn)行數(shù)值模擬，并將數(shù)值模擬中的能量吸收、沖擊力峰值以及破壞形式與實(shí)驗(yàn)結(jié)果進(jìn)行了比較，結(jié)論如下:

圖5 試件1模擬破壞過(guò)程與實(shí)際破壞過(guò)程的比較Fig.5 Comparison between failure simulation and actual failure of specimen 1

圖6 試件2模擬破壞過(guò)程與實(shí)際破壞過(guò)程的比較Fig.6 Comparison between failure simulation and actual failure of specimen 2

(1)本文所采用的數(shù)值模擬方法能夠比較準(zhǔn)確的反映嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)在低速?zèng)_擊條件下的破壞過(guò)程。

(2)嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)在沖擊破壞過(guò)程中吸收的能量高于傳統(tǒng)的纖維增強(qiáng)樹(shù)脂基復(fù)合材料。

(3)與普通纖維增強(qiáng)型樹(shù)脂基復(fù)合材料相比，嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)在沖擊破壞過(guò)程中的作用力峰值遠(yuǎn)遠(yuǎn)高于普通的復(fù)合材料。

(4)在低速?zèng)_擊條件下，普通樹(shù)脂基纖維增強(qiáng)型復(fù)合材料會(huì)產(chǎn)生脆性破壞，而嵌入式共固化復(fù)合材料結(jié)構(gòu)的破壞呈現(xiàn)塑性材料的特點(diǎn)。

［1］梁 森，梁 磊，米 鵬.嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)的新進(jìn)展［J］.應(yīng)用力學(xué)學(xué)報(bào)，2010，27(4):767 －771.

［2］張少輝，陳花玲.國(guó)外纖維增強(qiáng)樹(shù)脂基復(fù)合材料阻尼研究綜述［J］.航空材料學(xué)報(bào)，2002，22(1):58 －62.

［3］Gibson R F，F(xiàn)inegan I C.Recentresearch on enhancement of damping in polymer composites［J］.Composite Structures，1999，44(2):89 －98.

［4］李 烜，梁 森，常園園，等.嵌入式共固化復(fù)合材料阻尼結(jié)構(gòu)阻尼性能的實(shí)驗(yàn)研究［J］.科學(xué)技術(shù)與工程，2010，10(6):1510－1513.

［5］張少輝，陳花玲.共固化復(fù)合材料粘彈阻尼結(jié)構(gòu)的損耗因子研究［J］.航空材料學(xué)報(bào)，2005，25(1):53－57.

［6］時(shí)黨勇，李裕春，張勝民.基于ANSYS/LS-DYNA8.1進(jìn)行顯示動(dòng)力學(xué)分析［M］.北京:清華大學(xué)出版社，2005:119－130.

［7］Loikkanen M，Praveen G，Powell D.Simulation of ballistic impact on composite panels［C］.10thInternational LS-DYNA User Conference，Detroit:Livermore Software Technology Corporation，2008.

［8］Bhushan S，Thatte G S，Chandekar A D K，et al.Studies on behavior of carbon and fiber glass epoxy composite laminate under low velocity impact loading using LS-DYNA［C］.10thInternational LS-DYNA User Conference，Detroit:Livermore Software Technology Corporation，2008.

［9］Chang F K，Chang K Y.A progressive damage model for laminated composites containing stress concentration［J］.Journal of Composite Materials，1988，19(4):834 －855.

［10］Hallquist J O. LS-DYNA theoretical manual［M］.Livermore:Livermore Software Technology Corporation，1998:227－245.

［11］Hallquist J O.LS-DYNA keyword user's manual［M］.Livermore:Livermore Software Technology Corporation，1998:843－846.

［12］Palazotto A N，Herup E J，Gummadi L N B.Finite element analysis of low-velocity impact on composite sandwich plates［J］.Composite Structures，2000，49(2):209 －227.

［13］Liang S，Zhang Y X.Impact performance of embedded and co-cured composite damping panels［C］.16th International Conference on Composite Structures ICCS 16 FEUP，Porto，2011.