中圖分類號:TB332 DOI: 10.16579/j. issn.1001.9669.2025.07.011
0 引言
復合材料具有強度高、密度低、疲勞性能好等優(yōu)點,在航空航天、新能源領域得到廣泛應用[1-3]。復合材料層合板由單向碳纖維與樹脂基體通過層合熱壓而成,在使用過程中可能會受到冰雹、碎片等撞擊。這會導致復合材料層合板的剩余強度和抗疲勞性能顯著下降,削減了材料的使用性能[4]。因此,研究復合材料層合板疲勞壽命成為國內(nèi)外研究的熱點之一[5]。復合材料層合板疲勞試驗周期長,耗費大量的物力,其最直接、有效的方法就是通過建立疲勞壽命預測模型對疲勞性能進行評價。
復合材料層合板的疲勞損傷機制十分復雜[6-7]在疲勞過程中,各種損傷相互影響,其材料性能逐漸退化。復合材料層合板的疲勞壽命預測模型分為三大類: ① 根據(jù)S-N曲線來預測特定載荷下的疲勞壽命; ② 根據(jù)復合材料層合板疲勞漸進損傷模型來預測疲勞壽命; ③ 根據(jù)剩余強度或剩余模量建立疲勞壽命預測模型。FAWAZ等將S-N疲勞壽命曲線變換為對數(shù)形式的表達式。SHEN等[0]研究了準靜態(tài)和疲勞載荷作用下正交鋪設的玻璃纖維復合材料層合板損傷演化規(guī)律,分析不同鋪設順序對材料損傷演化的影響。ATTIA等[\"研究了沖擊后復合材料在疲勞載荷作用下的損傷擴展問題,利用應變能釋放率來預測沖擊損傷后復合材料的疲勞壽命。徐穎等針對沖擊后復合材料層合板,研究了沖擊損傷層合板的壓-壓疲勞壽命預測方法。劉儉輝等[13修正了應力場強法,建立了疲勞壽命預測模型,并驗證了模型的可靠性。
復合材料層合板疲勞損傷過程十分復雜,得到準確的疲勞壽命通常需要做大量的試驗。為了降低試驗成本,提高預測模型的精度,考慮沖擊損傷層合板的壓縮剩余強度,選用開口等效法[1417]將損傷區(qū)域等效為圓形開孔,圓的直徑取決于最大損傷尺寸,并認為該區(qū)域不再具有承載能力。此外,子層屈曲法18]將沖擊損傷看作多個矩形分層,并認為壓縮破壞是各個子層發(fā)生屈曲失效的過程,該方法在工程實際中有一定的局限性。軟化夾雜法[19-21]將沖擊損傷等效為形狀規(guī)則的軟化夾雜,然后用應力準則、應變準則或其他準則判定層合板的失效。損傷累積法[22-25]利用動態(tài)有限元模擬層合板的沖擊損傷,將對應退化后的剛度作為層合板的初始損傷,再用損傷累積法模擬層合板的壓縮破壞過程,并計算剩余壓縮強度。
借助低速沖擊試驗分析損傷區(qū)域尺寸,并進行等效孔徑計算,通過壓縮試驗得到壓縮剩余強度,進而求解外載荷應力,通過15J沖擊能量下的壓-壓疲勞試驗獲得其疲勞壽命?;谝陨辖Y果,擬合等效損傷系數(shù),提出適用于低速沖擊損傷狀態(tài)下的T300/69層合板疲勞壽命預測模型,為工程應用提供理論依據(jù)。本文只考慮沖擊能量對疲勞壽命的影響,未考慮鋪層方式、鋪層角度、沖擊大小及方向對疲勞性能的影響。
1低速沖擊后壓縮試驗
1.1 T300/69試驗件
T300/69 碳纖維復合材料層合板采用熱壓罐工藝制作而成,單層板厚度為 0.15mm ,材料鋪層順序為[0/45/-45/90]3s ,共24層。使用線切割技術將復合材料層合板切割成長為 75mm 、寬為 35mm 的試驗件。T300/69 層合板彈性性能參數(shù)和強度性能參數(shù)如表1所示。
1.2基于沖擊試驗獲取凹陷損傷數(shù)據(jù)
如圖1所示,沖擊試驗在落錘試驗臺上完成,通過調(diào)整落錘高度來獲得沖擊能量,沖頭的形狀可以根據(jù)實際的工況進行調(diào)整。
圓形凹陷,層間出現(xiàn)少量分層;當沖擊能量為15J時,損傷區(qū)域的凹陷程度更明顯,其部位產(chǎn)生纖維翹起等現(xiàn)象(圖2)。
由試驗可以得到不同沖擊能量下層合板損傷形貌尺寸,如表2所示。沖擊產(chǎn)生的凹坑長度和寬度隨沖擊能量增大而增加。相比于凹陷寬度,損傷區(qū)域凹陷長度對于沖擊能量的變化更加敏感。觀察層合板的損傷區(qū)域發(fā)現(xiàn),其主要損傷模式有纖維壓縮破壞、基體擠壓及基體裂紋擴展導致的分層。
1.3基于低速沖擊后壓縮試驗獲取壓縮剩余強度
按照ASTMD7137[26標準,壓縮試驗在WDW-300微機控制電子萬能試驗機上進行。壓縮試驗采用載荷控制的加載方式,壓縮試驗加載速率為 2mm/min 試件夾持情況如圖3所示。
復合材料層合板邊緣沖擊的損傷形貌與沖擊能量有顯著關系。當沖擊能量為5J時,損傷區(qū)域產(chǎn)生橢
損傷區(qū)域隨沖擊能量的增大而增大;當沖擊能量為5J時,壓縮失效試件基本觀察不到因壓縮失效而產(chǎn)生的裂紋[圖4(a);當沖擊能量為15J時,壓縮失效試件表面產(chǎn)生裂紋[圖4(d)],在沖擊點附近已經(jīng)造成了剪切損傷。由試件壓縮側面損傷圖可以看出,邊緣沖擊后層合板壓縮失效主要是以分層損傷為主,產(chǎn)生的分層對稱分布在層合板的兩側。分層損傷是由邊緣低能量沖擊造成試件的損傷起始,使得層合板在壓縮載荷的作用下不斷擴展,最終失效,
不同沖擊能量后,層合板最大壓縮載荷及壓縮剩余強度如表3所示。邊緣沖擊對層合板壓縮強度的影響較大;隨著沖擊能量的增加,壓縮剩余強度隨之降低。同時也做了無損傷壓縮試驗,無損傷的壓縮強度為 333.25MPa 。
2沖擊損傷層合板的失效理論分析
開孔等效法認為,沖擊損傷層合板壓縮失效模式和開孔層合板相似,都是由孔(沖擊損傷區(qū))邊緣應力集中導致纖維屈曲,從而使層合板發(fā)生壓縮破壞2的。
2.1 平均應力判據(jù)
平均應力破壞準則是將距離孔邊一定距離上的平均應力作為標準來表征開孔層合板的破壞。假設在 y 軸上從 r 到 r+a0 這段距離的平均應力 σy 等于或者大于無缺口層合板的強度 σ0 ,就認為開孔層合板失效。
式中, σ?N∞ 為遠場應力,即壓縮剩余強度; KT∞ 為開孔層合板孔邊應力集中系數(shù); r 為孔徑; a0 為特征長度。
2.2 等效孔徑與特征長度的確定
運用開口等效法將沖擊造成的損傷處理為一圓 形開孔,根據(jù)表3將損傷區(qū)域等效為半圓孔,凹陷長度 等效為開孔孔徑。
張彥等[28]認為,含孔復合材料層合板特征長度 a0 與層合板的材料體系有關。
通過層合板沖擊后的靜壓縮強度和平均失效準則便可以推導出不同孔徑層合板的特征長度。損傷寬度擴展到層合板寬度的 50% 時,試件基本失去承載能力,經(jīng)過少數(shù)加載循環(huán)后就會失效。試驗表明,選擇孔徑,到 50% 板寬距離的平均應力會更適用[29]。如表4所示,沖擊能量越大,等效孔徑越大,特征長度越小。
3壓-壓疲勞試驗及壽命預測模型
3.1基于壓-壓疲勞試驗獲取疲勞壽命
復合材料層合板低速沖擊后的壓-壓疲勞試驗在試驗機(圖5)進行,試驗均在室溫下進行,疲勞加載采用載荷控制的加載方式,載荷波形為等幅正弦波。
T300/69 復合材料層合板沖擊后壓-壓疲勞損傷如圖6所示。由圖6可以看出,應力水平越大,其疲勞損傷越嚴重,層間分層越明顯。在壓-壓疲勞載荷作用下,試件邊緣沖擊造成的分層由沖擊點附近沿寬度方向向另一側擴展,這個過程十分緩慢。隨著循環(huán)次數(shù)的增加,直至臨近破壞時,整個試件才會發(fā)生貫穿式分層和子層屈曲現(xiàn)象。
由表5可知,即使是同一沖擊能量、同一應力水平,疲勞壽命也存在分散性。應力水平越高,循環(huán)次數(shù)越少。從邊緣沖擊試驗到壓-壓疲勞試驗,層合板局部損傷的復雜多樣性也對疲勞壽命產(chǎn)生了影響。試驗發(fā)現(xiàn),15J能量沖擊后層合板在 80% 應力水平下壽命相差較大。
3.2 數(shù)值分析
建立層合板沖擊后的疲勞壽命模型:
式中, σy(x,0) 與 σs 分別為開孔層合板的外載荷應力和無損傷層合板的靜壓縮強度。為了更加真實地反映出應力-壽命情況,根據(jù)等應力等損傷原則對開孔層合板的疲勞壽命模型進行修正,等效損傷系數(shù)為α ,式(2)變?yōu)?/p>
根據(jù)無損傷層合板疲勞壽命模型的指數(shù)函數(shù)規(guī)律,假設損傷層合板的對數(shù)壽命與應力水平之間也存在指數(shù)函數(shù)規(guī)律,用試驗數(shù)據(jù)進行擬合,得到了低速沖擊后層合板的疲勞壽命:
1)計算孔徑應力集中系數(shù) KT∞ ,即
2)計算外載荷應力 σy ,即
等效孔徑 r 和特征長度 a0 由表4確定,將5J、15J能量沖擊后的層合板的疲勞壽命代入式(4)得到 α 值,如表6所示。
將15J沖擊能量在不同應力水平下獲得的等效損傷系數(shù)求和再取平均值,得到的疲勞模型等效損傷系數(shù) α=0.4970 ,則疲勞壽命模型為
4試驗驗證與分析
4.1 試驗驗證
如圖7所示,將15J能量沖擊后復合材料層合板所獲得的試驗結果與預測結果進行對比。結果表明,誤差范圍較小,具體誤差分析見第4.2節(jié)。平均應力水平為 60% 時,相對誤差最大為 6.541% ;平均應力水平為 70% 時,相對誤差最大為 5.896% 。
如圖8所示,將低速沖擊損傷層合板的疲勞壽命試驗值與預測值進行了對比。低速沖擊損傷層合板的疲勞壽命預測結果說明,本文所采用的疲勞壽命預測方法是合理的。
4.2 誤差分析
誤差來源可以歸納為以下幾點: ① 復合材料層合板的疲勞壽命存在一定的分散性,疲勞性能易受到?jīng)_擊載荷、環(huán)境溫度和濕度的影響。邊緣沖擊后的壓-壓疲勞試驗,層合板局部損傷的復雜多樣性也對疲勞壽命產(chǎn)生了影響。 ② 試驗發(fā)現(xiàn),層合板在 80% 應力水平下的壽命相差較大。一方面,因為層合板對分層損傷較為敏感,15J能量的低速沖擊造成的分層損傷嚴重影響了復合材料層合板的疲勞性能;另一方面,邊緣沖擊能量越高,層合板損傷區(qū)域越大,層合板承受的壓縮載荷變?nèi)酰瑧惺沟脫p傷區(qū)域的應力變大,從而影響層合板的疲勞壽命。 ③ 碳纖維復合材料本身為各向異性材料,存在一定的分散性。 ④ 將沖擊損傷區(qū)域等效為圓形孔,雖然兩者的力學特性相似但并不完全相同,導致等效過程存在差異。
5結論
針對低速邊緣沖擊下的 T300/69 復合材料層合板疲勞壽命開展一些研究。考慮層合板壓縮剩余強度對疲勞壽命的影響,基于開口等效法建立了疲勞壽命預測模型;通過預測值與試驗值,驗證了該方法的合理性和可靠性,得出以下結論:
1)考慮低速沖擊后復合材料層合板復雜的損傷機制,將邊緣沖擊損傷層合板等效為開孔層合板,通過測量凹陷損傷尺寸,確定開孔孔徑及特征長度,其計算過程簡潔高效。
2)考慮壓縮剩余強度對疲勞壽命的影響,結合平均應力準則計算外載荷應力,通過試驗數(shù)據(jù)擬合疲勞壽命模型,得到損傷參數(shù),從而獲得 T300/69 復合材料層合板的疲勞壽命預測模型。
3)將疲勞壽命預測與疲勞壽命試驗相比,驗證該方法可以準確地預測邊緣低速沖擊下的層合板疲勞壽命,且預測精度較高,預測誤差在 10% 以內(nèi)。
參考文獻(References)
[1] 葉林,賀鵬飛,李文曉.含z向纖維增強的連續(xù)纖維增強樹脂基 復合材料沖擊性能研究進展[J].高科技纖維與應用,2021,46 (3):44-52. YE Lin,HE Pengfei,LI Wenxiao. A review on impact property of continuous fiber reinforced resin matrix compositeswith z -direction fiber reinforced[J].Hi-Tech Fiber and Application,2021,46(3): 44-52. (In Chinese)
[2]陳縣輝,張彥飛,杜瑞奎,等.沖擊方向對層合板低速沖擊響應的 模擬研究[J].工程塑料應用,2014,42(5):89-92. CHEN Xianhui, ZHANG Yanfei,DU Ruikui,et al. Simulation of effects of impact orientations on low velocity impact responses of composite laminates[J].Engineering Plastics Application,2014,42 (5):89-92.(In Chinese)
[3]RHEAD A T,MARCHANT D,BUTLER R. Compressive strength of composite laminates following free edge impact[J]. Composites Part A:Applied Science and Manufacturing,2010,41(9) : 1056- 1065.
[4] 賈耀雄,敖清陽,張文正,等.碳纖維復合材料層壓板低速沖擊損 傷性能分析[J].兵器材料科學與工程,2022,45(5):170-174. JIA Yaoxiong,AO Qingyang,ZHANG Wenzheng,et al. Analysis on low-velocity impact damage properties of carbon fiber compositelaminates[J].Ordnance Material Science and Engineering, 2022,45(5):170-174.(In Chinese)
[5]俞鳴明,朱雪莉,劉雪強,等.低速多次沖擊下碳纖維/環(huán)氧樹脂 基復合材料層合板失效機制及剩余強度評估[J].復合材料學 報,2023,40(9):5359-5370. YU Mingming,ZHU Xueli,LIU Xueqiang,et al. Failure mechanism and assessment of residual strength of carbon fiber/epoxy resin matrix composite laminates under multiple impacts at low velocities[J].Acta Materiae Compositae Sinica,2023,40(9): 5359-5370.(In Chinese)
[6]MAHMOUDI A,MOHAMMADI B,HOSSEINI-TOUDESHKY H T. Damage behaviour of laminated composites during fatigue loading[J].Fatigueamp;Fracture of EngineeringMaterialsamp; Structures, 2020,43(4):698-710.
[7]ZUO P Y,SRINIVASAN D V, VASSILOPOULOS A P. Review of hybrid composites fatigue[J].Composite Structures,2021,274: 114358.
[8]DEGRIECK J, VAN PAEPEGEM W. Fatigue damage modeling of fibre-reinforced composite materials: review[J].Applied Mechanics Reviews,2001,54(4):279-299.
[9]FAWAZ Z,ELLYIN F. Fatigue failure model for fibre-reinforced materials under general loading conditions[J]. Journal of Composite Materials,1994,28(15):1432-1451.
[10]SHEN HJ,YAO W X,QI W X,et al. Experimental investigation on damage evolution in cros-ply laminates subjected to quasi-static and fatigue loading[J]. Composites Part B: Engineering,2017, 120:10-26.
[11]ATTIAO,KINLOCHAJ,MATTHEWSFL.Thepredictionof fatigue damage growth in impact-damaged composite skin/stringer structures.part I: theoretical modelling studies [J].Composites Science and Technology,2003,63(10):1463-1472.
[12]]徐穎,溫衛(wèi)東,崔海濤.復合材料層合板沖擊后壓-壓疲勞壽命預 測方法[J].復合材料學報,2007,24(2):159-167. XU Ying,WEN Weidong,CUI Haitao. Prediction method for fatiguelifeof laminated compositesafterimpactunder compression-compression loading[J]. Acta Materiae Compositae Sinica,2007,24(2):159-167.(In Chinese)
[13]劉儉輝,薛文卓,魏泰.含沖擊損傷復合材料層合板疲勞壽命預 測[J].西安理工大學學報,2022,38(2):164-170. LIU Jianhui,XUE Wenzhuo,WEI Tai. Fatigue life prediction of composite laminates with impact damage[J].Journal of Xi'an UniversityofTechnology,2022,38(2):164-170.(In Chinese)
[14]CHEN P H,SHEN Z,WANG JY.A new method for compression after impact strength prediction of composite laminates[J]. Journal of Composite Materials,2002,36(5):589-610.
[15]SOUTIS C,CURTIS P T. Prediction of the post-impact compressive strength of CFRP laminated composites[J]. Composites Science and Technology,1996,56(6):677-684.
[16]SOUTIS C,SMITH F C,MATTHEWS F L. Predicting the compressive engineering performance of carbon fibre-reinforced plastics [J].Composites Part A:Applied Science and Manufacturing, 2000,31(6) :531-536.
[17]童谷生,孫良新,劉英衛(wèi),等.復合材料層壓板低能量沖擊后剩余 抗壓強度的工程估算[J].機械工程材料,2004,28(3):19-21. TONG Gusheng,SUN Liangxin,LIU Yingwei,et al. A simplified method of evaluating compression after impact strength for composite laminate[J].Materials for Mechanical Engineering,2004,28 (3):19-21. (In Chinese)
[18]KUTLU Z,CHANG F K. Composite panels containing multiple through-the-width delaminations and subjected to compression. Part I:analysis[J]. Composite Structures,1995,31(4):273-296.
[19]JOHNSON W S,MASTERS JE,O'BRIEN TK,et al. Residual strength of impacted composites: analysis and tests[J]. Journal of Composites Technology and Research,1994,16(3):244-255.
[20]XIONG Y,POON C,STRAZNICKY P V,et al.A prediction method for the compressive strength of impact damaged composite laminates[J]. Composite Structures,1995,30(4):357-367.
[21]燕瑛,曾東.復合材料層板低速沖擊剩余強度的研究[J].航空學 報,2003,24(2):137-139. YAN Ying,ZENG Dong. Study on the post-impact compressive strength of composite laminates[J]. Acta Aeronautica et Astronautica Sinica,2003,24(2):137-139.(In Chinese)
[22]程小全,正能.復合材料層合板低速沖擊后壓縮的損傷累積模 型[J].應用數(shù)學和力學,2005(5):569-576. CHENG Xiaoquan,LI Zhengneng. Damage progressive model of comperssion composite laminates after lowvelocity impact[J]. Applied Mathematics and Mechanics,2005(5):569-576.(In Chinese)
[23]崔海坡,溫衛(wèi)東,崔海濤.復合材料層合板沖擊損傷及剩余強度 研究進展[J].材料科學與工程學報,2005,23(3):466-472. CUIHaipo,WEN Weidong,CUI Haitao.Advancesof study on lowvelocity impact damaged and residual strength of composite laminates[J].Journal ofMaterialsScienceand Engineering,2005, 23(3):466-472.(In Chinese)
[24]崔海坡,溫衛(wèi)東,崔海濤.復合材料層合板沖擊損傷及剩余強度 分析方法[J].固體力學學報,2006,27(3):237-242. CUI Haipo,WENWeidong,CUI Haitao.Research on low velocity impact damageand residual strength ofcomposite laminates[J]. ActaMechanica Solida Sinica,2006,27(3):237-242.(In Chinese)
[25]崔海坡,溫衛(wèi)東,崔海濤.層合復合材料板的低速沖擊損傷及剩 余壓縮強度研究[J].機械科學與技術,2006,25(9):1013-1017. CUI Haipo,WENWeidong,CUI Haitao.Research on lowvelocity impact damage and residual compressive strength of laminated composites[J].Mechanical Science and Technology,2006,25(9): 1013-1017.(InChinese)
[26]UDA N,ONO K,KUNOO K.Compression fatigue failure of CFRPlaminateswith impactdamage[J].Composites Science and Technology,2009,69(14):2308-2314.
[27]梁小林.復合材料層合板沖擊后的疲勞壽命研究[D].南京:南 京航空航天大學,2014:36-38. LIANG Xiaolin.The study on fatigue life of composite laminates after impact[D]. Nanjing:NanjingUniversity of Aeronautics and Astronautics,2012:36-38.(In Chinese)
[28]張彥,朱平,來新民,等.低速沖擊作用下碳纖維復合材料鋪層板 的損傷分析[J].復合材料學報,2006,23(2):150-157. ZHANG Yan,ZHU Ping,LAI Xinmin,etal.Low-velocity impact damageanalysisin carbon fiber reinforced composite laminated plates[J]. Acta Materiae Compositae Sinica,2006,23(2):150- 157.(In Chinese)
[29]朱煒垚.含低速沖擊損傷復合材料層板剩余強度及疲勞性能研 究[D].南京:南京航空航天大學,2012:99-101. ZHUWeiyao.Research on residual strength and fatigue performance of composite laminates with low-velocity impact damage [D].Nanjing:Nanjing University of Aeronautics and Astronautics, 2012:99-101.(In Chinese)
Abstract:Afterlow-velocityimpactat theedge,delamination and matrix extrusion occurinside thecomposite laminates, which willhaveaseriousimpactonthesafeuseand lifeofthecompositelaminates.Therefore,itisofpractical engieering significancetoestablishafatiguelifepredictionmodelforlow-velocityimpactattheedge.Thedentdamagesize,copressive residualstrengthandfatiguelifeofthefatiguelifepredictionmodelwereobtainedbylow-spedimpacttest,compressiontest andcompresson-compressionfatigue test.Basedontheaverage stress failure criterion,theimpactdamageareaof the laminated platewasequivalenttothecorrespondingaperturebycombiningtheopeningequivalentmethod,andtheequivalent damage coeffcient of diferent impact energy was proposed.A fatigue life prediction model considering the compressive residual strengthofimpactdamagedlaminates wasestablished,andthe predictionresults werecompared with thetestresults. The resultsshow that the fatigue life prediction accuracy of the model is high,the eror is controlled within 10% ,andthe model has good prediction ability.
Keywords:Compositelaminate;Fatigue lifepredictionmodel;Opening equivalentmethod;Low-velocityedge impact
Correspondingauthor:LIUJianhui,E-mail: liujh@lut.edu.cn
Fund:National Natural Science Foundation of China (51605212)
Received:2023-11-23 Revised: 2024-01-15