羅書舟, 陳 超, 伍乾坤, 胡德安
(湖南大學(xué) 特種裝備先進(jìn)設(shè)計技術(shù)與仿真教育部重點(diǎn)實驗室,長沙 410082)
近幾十年來,隨著中國經(jīng)濟(jì)的飛速發(fā)展,汽車工業(yè)特別是家用轎車的發(fā)展越來越快。汽車的大量使用給人們帶來了便利,但是也帶來了如能源過度消耗、環(huán)境污染等問題。目前,節(jié)能減排、低碳生活成為全球汽車工業(yè)的研究熱點(diǎn),汽車輕量化是實現(xiàn)汽車低油耗和低排放的重要途徑,而使用輕質(zhì)高效能的材料可有效的實現(xiàn)汽車輕量化[1],如碳纖維復(fù)合材料在汽車工業(yè)中的應(yīng)用。在工程實際中,復(fù)合材料結(jié)構(gòu)需離型制造,再進(jìn)行裝配,所以復(fù)合材料間的連接性能至關(guān)重要[2]。與傳統(tǒng)的機(jī)械連接相比,膠接具有結(jié)構(gòu)輕、減震性能好等優(yōu)點(diǎn),越來越廣泛地應(yīng)用于復(fù)合材料結(jié)構(gòu)件間的連接。因此,對復(fù)合材料膠接接頭的研究越來越引起學(xué)者和工程師的重視。
從20世紀(jì)30年代開始,國內(nèi)外的學(xué)者就對膠接接頭進(jìn)行了理論研究[3-9]。近年來,陳煊等[10]探討了板厚、溫度和試樣端部的加載速度對單搭接膠接接頭強(qiáng)度的影響。鄭小玲等[11]探討了膠瘤對單搭接膠接接頭強(qiáng)度的影響。Kadioglu等[12]實驗探討了單搭接接頭在擺錘沖擊下的力學(xué)性能。梁祖典等針對T300碳纖維層合板單搭接膠接接頭進(jìn)行了實驗研究和數(shù)值模擬,得到了膠接連接的失效模式和極限載荷與膠接長度和膠接厚度密切相關(guān)的結(jié)論,并分析了膠接連接在拉伸載荷下的破壞起始和演化過程。Stazi等[13]通過實驗和仿真研究了不同濕度、水分、紫外線條件下的不同膠的拉伸性能,得到了外界條件對單搭接接頭拉伸性能的影響。Wu等[14]通過三點(diǎn)彎曲和沖擊實驗,探討了膠接長度和寬度對單搭接接頭的力學(xué)性能的影響。Shi等[15]通過實驗和仿真研究了復(fù)合材料在低速沖擊時的面內(nèi)損傷和層間損傷,分析了復(fù)合材料中裂紋的演化過程。開展復(fù)合材料沖擊實驗和仿真的研究越來越多,然而,對于復(fù)合材料膠接接頭在不同膠參數(shù)下的失效模式、破壞形貌和能量吸收缺乏對比研究。
本文針對HTS40/977-2碳纖維層合板單搭接接頭在三種不同膠膠接下承受沖擊載荷的特性進(jìn)行研究。首先對碳纖維層合板低速沖擊進(jìn)行數(shù)值模擬,通過與實驗數(shù)據(jù)對比,驗證了材料模型參數(shù)的合理性。其次,建立了采用不同膠的碳纖維層合板單搭接接頭的有限元模型,基于Hashin準(zhǔn)則的本構(gòu)關(guān)系預(yù)測層合板的面內(nèi)損傷,應(yīng)用內(nèi)聚力單元分析層合板的分層及不同膠膠接的失效。最后,基于數(shù)值模擬結(jié)果分析不同膠膠接接頭的失效模式和能量吸收特性。
論文采用商業(yè)軟件ABAQUS/Explicit進(jìn)行數(shù)值模擬,有限元模型考慮面內(nèi)損傷(基體和纖維的損傷)、層間分離和膠層損傷。復(fù)合材料面內(nèi)損傷通過Hashin準(zhǔn)則來預(yù)測,層間分離和膠層損傷引入內(nèi)聚力單元進(jìn)行分析。
Hashin準(zhǔn)則考慮如下四種模式的損傷失效[16]。
纖維損傷
(1)
(2)
基體損傷
(3)
(4)
圖1 纖維拉伸雙線性損傷演化模型
本文所研究的碳纖維層合板HTS40/977-25是用單層板層層堆疊后加壓固化而成的,兩相鄰單層板之間的界面在沖擊載荷作用下會產(chǎn)生損傷并擴(kuò)展。論文引入內(nèi)聚力(Cohesive)單元研究相鄰單層板之間的分層損傷和擴(kuò)展,同時層合板單搭接膠層也采用內(nèi)聚力單元進(jìn)行模擬。
內(nèi)聚力單元示意圖如圖2所示,由連接相鄰單層板的上下兩個表面組成,這兩個表面沿厚度方向的運(yùn)動表征內(nèi)聚力單元的張開和閉合,沿著垂直于厚度方向的運(yùn)動表征內(nèi)聚力單元的橫向剪切行為。內(nèi)聚力單元包括三個應(yīng)力,即厚度方向的正應(yīng)力σn和垂直于厚度方向的兩個切應(yīng)力σs和σt,δn、δs和δt分別為三個應(yīng)力所對應(yīng)的位移。
圖2 內(nèi)聚力單元示意圖
通常情況下內(nèi)聚力單元的厚度很小,忽略法向變形幾乎不影響切向的變形。內(nèi)聚力單元線彈性狀態(tài)下本構(gòu)關(guān)系的表達(dá)式為
(5)
本文內(nèi)聚力單元采用雙線性本構(gòu)模型,如圖3所示。當(dāng)δ=δ0時,表示材料開始產(chǎn)生失效;當(dāng)δ≥δmax時,表示材料已經(jīng)完全失效;當(dāng)δ0≤δ≤δmax時,表示分層失效開始擴(kuò)展。
圖3 內(nèi)聚力單元的雙線性本構(gòu)模型曲線
內(nèi)聚力單元的損傷主要由損傷的起始和擴(kuò)展這兩個階段構(gòu)成。損傷起始判據(jù)采用平方應(yīng)力準(zhǔn)則,當(dāng)對應(yīng)的界面應(yīng)力和臨界界面強(qiáng)度比的平方達(dá)到1時,表示損傷的起始。損傷擴(kuò)展采用考慮了不同裂紋模式應(yīng)變能釋放率的B-K準(zhǔn)則來判定。式(6)和式(7)分別給出了平方應(yīng)力準(zhǔn)則和B-K準(zhǔn)則的具體形式。平方應(yīng)力準(zhǔn)則為
(6)
(7)
式(7)分別代表材料在受拉和受剪作用下的損傷破壞。B-K準(zhǔn)則[17]為
(8)
基于前面所述的面內(nèi)損傷和層間損傷模型,參照文獻(xiàn)[18]中的實驗,在ABAQUS中建立碳纖維復(fù)合材料層合板HTS40/977-2的沖擊模型。沖擊實驗滿足ASTM D7136/D7136M-07標(biāo)準(zhǔn)[19]。層合板厚度為2 mm,鋪層方式為[0/90]2s,實驗試樣由兩塊中空鋼板組成的夾具夾持,外露直徑為75 mm的圓形區(qū)域。沖頭為半球型,直徑為15 mm,整個實驗過程中沖頭跌落高度為0.75 mm,沖頭的質(zhì)量分別為1 kg、1.5 kg和2 kg,沖擊能量分別為7.35 J、11.03 J和14.7 J。
論文建立碳纖維復(fù)合材料層合板HTS40/977-2承受沖擊載荷的有限元模型,如圖4所示,碳纖維層合板由8層單層板組成,每層厚度為0.25 mm,采用減縮積分8節(jié)點(diǎn)連續(xù)殼單元SC8R進(jìn)行離散,其材料參數(shù)如表1所示。碳纖維層合板層間加入有厚度的內(nèi)聚力單元,每層厚度為0.000 75 mm,內(nèi)聚力單元的單元類型為COH3D8,材料屬性見表2。連續(xù)殼單元和內(nèi)聚力單元的單元刪除損傷參量極限值設(shè)為0.99。層合板邊界設(shè)置為6個自由度完全固定,沖頭定義為剛體,限制住除z方向之外的所有自由度,沖擊速度為3.83 m/s,通過改變沖頭的質(zhì)量改變沖擊能量。
圖4 層合板沖擊模型
參數(shù)E1/GPaE2/GPaE3/GPav12v13數(shù)值15310.310.30.30.3參數(shù)v23G12/GPaG13/GPaG23/GPa數(shù)值0.410.310.30.3
表2 內(nèi)聚力單元材料參數(shù)
仿真結(jié)果和試驗結(jié)果的對比見圖5,圖中力-時間曲線是從剛性沖頭上提取得到的。從圖中可以看出,沖頭在達(dá)到最大反作用力值后,其反作用力慢慢減小至0,該階段為回彈過程,與實驗現(xiàn)象一致。圖5(a)中,在沖頭能量為7.35 J時,實驗得到的反作用力峰值3 383 N與模擬得到的峰值3 277 N較吻合。圖5(b)中,在沖頭能量為11.03 J時,實驗得到的反作用力峰值3 765 N,略高于模擬得到的峰值3 388 N。圖5(c)中顯示了在沖頭能量為14.7 J下模擬得到的反作用力峰值4 358 N,該峰值略低于實驗值4 605 N。而論文[18]中模擬得到的最大反作用力分別為3 116 N、3 338 N、3 917 N,在不同沖擊能量下本文模擬精度均高于文獻(xiàn)[18]。以上分析表明,在不同的沖擊能量下,本文仿真得到的力-時間曲線和實驗測量的力-時間曲線基本吻合,驗證了碳纖維復(fù)合材料層合板HTS40/977-2材料參數(shù)的有效性。
(a) 沖擊能量7.35 J
(b) 沖擊能量11.03 J
(c) 沖擊能量14.7 J
為研究不同膠對單搭接接頭力學(xué)性能的影響,結(jié)合復(fù)合材料在汽車結(jié)構(gòu)上的應(yīng)用,設(shè)計出的單搭接接頭模型如圖6所示,該模型由兩塊HTS40/977-2層合板膠接而成,層合板鋪層方式為[0/90]4s,每塊層合板厚度為4 mm,長寬為115 mm×30 mm,膠層厚度為0.2 mm。試件具體尺寸如圖6所示。
碳纖維復(fù)合材料由16層單層材料組成,對單層材料采用8節(jié)點(diǎn)一階減縮積分單元C3D8R進(jìn)行離散。所采用的膠分別為Aradldite AV138、Araldite 2015和Sikaforce 7752,膠的韌度依次增大,膠層參數(shù)見表3。為了研究不同沖擊能量下不同膠的動態(tài)特性,選擇3 J和4 J的沖擊能量進(jìn)行分析,其中沖頭沖擊速度為1.5 m/s,通過改變沖頭的質(zhì)量改變沖擊能量。碳纖維復(fù)合材料層合板HTS40/977-2采用表1和2所示的經(jīng)過驗證的材料參數(shù)進(jìn)行模擬。
有限元模型一共有108 544個連續(xù)殼單元、1 104個實體單元和103 350個內(nèi)聚力單元。碳纖維層合板每層由連續(xù)殼單元離散,層間通過無厚度內(nèi)聚力單元連接,膠層采用厚度為0.2 mm的三維內(nèi)聚力單元模擬。邊界條件約束如圖7所示。模型中沖頭與試樣之間定義接觸,計算方法采用動力顯式格式進(jìn)行分析。
圖6 單搭接膠接接頭模型示意圖
本文數(shù)值模擬結(jié)果顯示單搭接膠接接頭主要失效形式有三種:層合板面內(nèi)失效、層合板層間分離和膠層失效。不同膠的膠接接頭在沖擊載荷下呈現(xiàn)出不同的破壞形貌。
表3 不同膠的材料參數(shù)[20]
圖7 有限元模型
通過數(shù)值模擬發(fā)現(xiàn),當(dāng)膠層為Aradldite AV138和沖擊能量為3 J時,如圖8(a)所示,膠層在沖頭的沖擊下,部分單元損傷系數(shù)達(dá)到極限值,膠層表面發(fā)生了破壞。上層合板邊界約束一端的膠層最先失效,然后向另一端擴(kuò)展,與Delale等的實驗現(xiàn)象吻合。當(dāng)沖擊能量為4 J時,膠層在沖頭的作用下,全部發(fā)生破壞,如圖8(b)所示。
當(dāng)膠層為Aradldite 2015和沖頭能量為3J時,模擬的破壞形貌如圖8(c)所示,膠層破壞也從左側(cè)端部開始,隨著沖擊的進(jìn)行,最終膠層破壞長度為10 mm,膠層長度共30 mm。當(dāng)沖擊能量增大到4 J時,膠層破壞長度增大到14 mm,如圖8(d)所示。
當(dāng)膠層為Sikaforce 7752時,破壞形貌如圖8(e)和8(f)所示。沖擊能量為3 J時,由于Sikaforce 7752膠層韌度大,損傷未到極限值,未出現(xiàn)膠層分離現(xiàn)象。當(dāng)沖擊能量為4 J時,也未出現(xiàn)脫膠。
對比三種膠膠接接頭在兩種沖擊能量下膠層的破壞形貌可以發(fā)現(xiàn),層合板受到沿厚度方向的沖擊時,由于上下層合板結(jié)構(gòu)的不對稱導(dǎo)致膠層受到的應(yīng)力不對稱,上層合板一端的膠層上表面應(yīng)力值高于膠層下表面應(yīng)力值,因此膠層失效從上層合板一端的膠層上表面開始。另外,Aradldite AV138、Aradldite 2015和Sikaforce 7752三種膠的韌度依次增加,故在沖頭的沖擊載荷作用下,Aradldite AV138、Aradldite 2015接頭膠層失效,且Aradldite AV138膠層的失效程度大于Aradldite 2015膠層,而沖擊能量未達(dá)到Sikaforce 7752膠層的失效判據(jù),故Sikaforce 7752膠層未發(fā)生脫膠現(xiàn)象。
(a) Aradldite AV138膠層-3 J
(b) Aradldite AV138膠層-4 J
(c) Aradldite 2015膠層-3 J
(d) Aradldite 2015膠層-4 J
(e) Sikaforce 7752膠層-3 J
(f) Sikaforce 7752膠層-4 J
圖8 膠接接頭的破壞形貌
Fig.8 Failure morphology of adhesive bonded single-lap joint
為了進(jìn)一步分析層合板的失效模式,圖9給出了單搭接膠接接頭的剖視圖。由于上層合板出現(xiàn)的層間分離現(xiàn)象更顯著,因此,對圖9剖視圖中的上層合板進(jìn)行局部放大,得到了三種膠接接頭在不同能量沖擊載荷下的破壞形貌??偟膩砜?,在3 J和4 J的沖擊能量下,層合板都出現(xiàn)了層間損傷,且損傷區(qū)域隨著能量的增加而增加,主要層間損傷出現(xiàn)在上層合板沖頭接觸區(qū)域附近。
圖9 剖視示意圖
當(dāng)膠層為Aradldite AV138時,在3 J和4 J的沖擊能量下,層合板在沖頭接觸附近區(qū)域出現(xiàn)較嚴(yán)重的層間失效,如圖10所示。同時,由于膠層破壞嚴(yán)重,上下層合板位移較大,貼近膠層一側(cè)的上層合板也出現(xiàn)了部分層間失效,并且在沖擊能量為4 J時,下層合板也出現(xiàn)了局部區(qū)域?qū)娱g失效。
當(dāng)膠層為Aradldite 2015時,層合板靠近沖頭附近的層間單元出現(xiàn)層間分離現(xiàn)象,如圖11(a)和11(b)所示,并且貼近膠層的一側(cè)的層合板也存在部分層間失效。與Aradldite AV138膠接比較,由于膠層失效更少,上下層合板位移較小,貼近膠層一側(cè)的層間失效區(qū)域更少。
當(dāng)膠層為Sikaforce 7752時,靠近沖頭附近的上層合板發(fā)生層間失效,沖擊能越大,失效單元越多。而膠層未失效,靠近膠層的層合板也未出現(xiàn)失效,如圖12(a)和12(b)所示。
(a) Aradldite AV138-3 J
(b) Aradldite AV138-4 J
Fig.10 Interlaminar failure of the upper laminate-Aradldite AV138
(a) Aradldite 2015-3 J
(b) Aradldite 2015-4 J
Fig.11 Interlaminar failure of the upper laminate-Aradldite 2015
(a) Sikaforce 7752-3 J
(b) Sikaforce 7752-4 J
Fig.12 Interlaminar failure of the upper laminate-Sikaforce 7752
比圖8、10、11和12三種不同膠的破壞形貌與層間失效可以發(fā)現(xiàn),Aradldite AV138膠的韌度是三種膠中最低的,其膠接接頭膠層全部失效,同時上下層合板層間均發(fā)生層間分離,并且上層合板層間分離區(qū)域面積遠(yuǎn)比下層合板大。對于Aradldit 2015膠接接頭,在沖頭沖擊下,部分膠層單元發(fā)生膠層失效,上層合板部分層間分離,下層合板幾乎沒有發(fā)生層間分離現(xiàn)象。Sikaforce 7752是三種膠中韌度最高的,在沖擊過程中,膠層未發(fā)生破壞,上層合板部分區(qū)域出現(xiàn)層間分離。以上分析可以看出,不同膠膠接接頭會明顯影響層合板的損傷形貌,在相同沖擊能量下,Aradldite AV138膠接層合板損傷最嚴(yán)重,其次是Aradldite 2015膠接層合板,Sikaforce 7752膠接層合板損傷程度最小。
在3 J和4 J能量沖擊下接觸力-位移曲線如圖13所示,從圖13可知,在沖擊過程中,沖頭都會回彈。沖頭的能量越大,沖頭的最大位移也增大。
膠接接頭吸收能量-時間曲線可以通過以下公式計算得到
(9)
式中:Ea(t)表示t時刻吸收的能量,vi表示沖頭開始時
(a) 沖擊能量3 J
(b) 沖擊能量4 J
的速度。vt表示t時刻沖頭沖擊的速度,δ(t)表示t時刻沖頭的位移。
吸收能量-時間曲線如圖14所示,從圖14可知,能量吸收-時間曲線在開始階段隨時間上升,達(dá)到峰值后下降并維持在一個恒定的水平。峰值代表相應(yīng)的沖頭能量,這一規(guī)律表明本文所建立的仿真模型能量守恒,最后恒定值代表試件最終吸收的能量。
(a) 沖擊能量3 J
(b) 沖擊能量4 J
三種膠膠接的試件吸收的能量見表4。從圖表中可以看出,膠接試件損傷越嚴(yán)重,其吸收的能量越多,而膠的韌性越大,接頭抵抗破壞的能力則越強(qiáng)。因此,膠的性能與膠接試件的力學(xué)性能密切相關(guān),選擇合適的膠在復(fù)合材料連接中至關(guān)重要。
表4 接頭吸收能量數(shù)據(jù)
論文研究了不同膠膠接接頭的低速沖擊性能,對不同膠膠接的HTS40/977-2復(fù)合材料層合板單搭接膠接接頭低速沖擊問題進(jìn)行了數(shù)值模擬。在低速沖擊時,膠的屬性會明顯影響接頭的失效模式和能量吸收。在相同沖擊能量條件下,膠層為Aradldite AV138和Aradldite 2015時,接頭失效模式為膠層和層合板失效。膠層為Aradldite AV138時,膠層完全失效;膠層為Aradldite 2015時,膠層部分失效;膠層為Sikaforce 7752時,失效模式僅為層合板失效。在本文三種膠中,膠的韌度越低,接頭損傷越嚴(yán)重,能量吸收越多;膠的韌度越高,接頭抵抗損傷的能力越強(qiáng),其能量吸收越少。本文研究可為復(fù)合材料膠接連接提供依據(jù)和參考。