層合多孔結(jié)構(gòu)的裂紋擴展規(guī)律研究

關(guān)鍵詞：層合多孔結(jié)構(gòu)；三點彎曲實驗；裂紋擴展；聚乳酸材料

0 引言

近年來，隨著輕質(zhì)結(jié)構(gòu)的迅速發(fā)展，多孔結(jié)構(gòu)因具有質(zhì)量輕、滲透性好、比強度高等優(yōu)點，從一系列輕質(zhì)結(jié)構(gòu)中脫穎而出，廣泛應(yīng)用于工程與醫(yī)療領(lǐng)域［1-2］。

迄今為止，現(xiàn)有的研究主要集中于多孔結(jié)構(gòu)的單層模型，大多由相同材料制成的周期性單元組成，許多學(xué)者已經(jīng)使用實驗仿真［3］和理論［4-6］方法探究了其裂紋擴展模式。張欣茹等［7-8］基于有限元仿真和實驗研究了拓?fù)湫螤?、單元方向、孔隙率、預(yù)制裂紋參數(shù)和相對密度對平面多孔結(jié)構(gòu)裂紋擴展的影響。GIBSON等［9］提出了一些經(jīng)驗公式來描述多孔結(jié)構(gòu)的斷裂韌性與其孔隙率、孔洞形狀和分布等結(jié)構(gòu)參數(shù)之間的關(guān)系。也有學(xué)者探究梯度多孔結(jié)構(gòu)的裂紋擴展模式［10-11］。SHU等［12］提出了在均勻蜂窩結(jié)構(gòu)中改變支板厚度的概念，通過非均質(zhì)設(shè)計提高均勻蜂窩結(jié)構(gòu)的韌性。對多孔復(fù)合材料也有較多的研究［13］。LIU等［14］以Si3N4晶須為界面層制備了仿生ZrB2-SiC層狀陶瓷，研究其抗彎強度和斷裂韌性。也有對復(fù)合材料層合板的研究［15-16］。趙麗濱等［17］從試驗研究、理論分析和數(shù)值模擬3個方面介紹了單向復(fù)合材料層間斷裂韌性測試方法和原理，對多向?qū)雍习宸謱訑U展行為進(jìn)行了試驗研究。牟浩蕾等［18］對纖維增強復(fù)合材料層合試驗樣件進(jìn)行拉伸、壓縮剪切試驗，分析其破壞模式。有學(xué)者將層合板理論引入多孔結(jié)構(gòu)中，建立層合多孔結(jié)構(gòu)。馬南芳［19］采用層合結(jié)構(gòu)的方式提高蜂窩結(jié)構(gòu)的能量吸收性能。賀璞等［20］研究了層合多孔圓柱殼在軸向壓縮載荷下的承載力性能和變形破壞模式。鮮有對多孔層合結(jié)構(gòu)裂紋擴展的研究。本文在多孔結(jié)構(gòu)的基礎(chǔ)上，以層合板理論為依托，建立層合多孔結(jié)構(gòu)模型，探究了層間效應(yīng)對多孔結(jié)構(gòu)裂紋擴展規(guī)律的影響。

本文基于復(fù)合材料力學(xué)層合板理論，設(shè)計了層合多孔結(jié)構(gòu)板，通過合理調(diào)整每層多孔結(jié)構(gòu)板的孔的構(gòu)型和鋪設(shè)方向獲取等效力學(xué)性能。利用SolidWorks軟件進(jìn)行模型設(shè)計，導(dǎo)入仿真軟件進(jìn)行模擬計算，并通過3D打印技術(shù)制備了實驗?zāi)Ｐ停ㄟ^實驗和仿真研究結(jié)果分析驗證設(shè)計和方法的合理性。為功能多孔結(jié)構(gòu)的設(shè)計和制備提供了一種新的思路和方法。

1 層合多孔結(jié)構(gòu)

采用復(fù)合材料層合板設(shè)計方法，以聚乳酸（PolylacticAcid，PLA）材料為基礎(chǔ)材料，設(shè)計了層合多孔結(jié)構(gòu)板。多孔結(jié)構(gòu)板的基本尺寸為140 mm×70 mm×10 mm（長×寬×厚），單層結(jié)構(gòu)厚度為10 mm，雙層結(jié)構(gòu)每層厚度為5 mm。本文引入了正六邊形和內(nèi)凹形2種基本構(gòu)型，并對構(gòu)型進(jìn)行旋轉(zhuǎn)操作，圖1為2種基本構(gòu)型示意圖，其夾角α 分別為120°和60°，基本構(gòu)型旋轉(zhuǎn)角度為β，其中，正六邊形旋轉(zhuǎn)角度分別為15°和30°，內(nèi)凹形則包括了15°、30°、45°和90° 4種不同的旋轉(zhuǎn)角度，基本構(gòu)型的幾何參數(shù)如表1所示。在雙層模型中，對所選取構(gòu)型進(jìn)行分層組合，模型及其簡稱如表2所示。為了分析裂紋的擴展規(guī)律，在下邊緣中部預(yù)置一個長11 mm、寬2 mm的三角形貫穿裂紋，如圖2所示。

2 實驗測試與有限元仿真

2. 1 實驗測試

采用3D打印機熔融絲進(jìn)行模型制備，打印速度為30 mm/s，如圖3所示。通過三點彎曲測試層合多孔結(jié)構(gòu)裂紋擴展規(guī)律，實驗裝置如圖4所示（以NAX-1C-0為例），使用微機控制電子萬能試驗機進(jìn)行測試，加載速度為1 mm/min，如圖5所示。

2. 2 有限元仿真

使用Abaqus有限元仿真軟件對平面多孔結(jié)構(gòu)板進(jìn)行了準(zhǔn)靜態(tài)三點彎曲模擬，并選用聚乳酸作為材料，材料參數(shù)如表3所示。實驗裝置底端兩支撐間距為120 mm。頂端壓頭及底端兩支撐與平面多孔結(jié)構(gòu)板之間的接觸形式均為摩擦接觸，摩擦因數(shù)為0. 3。采用四面體網(wǎng)格對平面多孔結(jié)構(gòu)板進(jìn)行劃分，模型劃分網(wǎng)格單元尺寸為2 mm，網(wǎng)格劃分如圖6所示。設(shè)置底端兩支撐面固定端約束，并對頂端壓頭施加分步位移載荷直至50 mm。

3 結(jié)果分析

3. 1 破壞形式

六邊形模型的實驗結(jié)果和有限元仿真結(jié)果如表4所示。由表4可知，在六邊形模型中，LBX-1C-15裂紋擴展方向與豎直方向的夾角大約為15°，LBX-2C-0-15胞元15°鋪層的裂紋擴展方向相對于LBX-1C-15的裂紋擴展方向向豎直方向偏轉(zhuǎn)；LBX-1C-30的裂紋擴展方向大致為豎直方向，整個裂紋擴展路徑呈“S”形，LBX-2C-0-30胞元30°鋪層的裂紋擴展方向為豎直方向，相較于LBX-1C-30的裂紋擴展方向更靠近豎直方向。其他六邊形模型的裂紋擴展方向都大致為豎直方向。由于失穩(wěn)現(xiàn)象，LBX-2C-0-15上方萌生了新裂紋，且每層新裂紋的擴展方向與預(yù)置裂紋擴展方向一致，其他六邊形模型均未出現(xiàn)此現(xiàn)象。

單層模型和雙層反對稱模型都在壓頭附近出現(xiàn)輕微致密化現(xiàn)象，LBX-2C-0-15和LBX-2C-0-30在壓頭與支座附近均出現(xiàn)明顯致密化現(xiàn)象；且由于失穩(wěn)現(xiàn)象的出現(xiàn)，在LBX-2C-0-15中，壓頭附近的致密化現(xiàn)象只出現(xiàn)在胞元15°鋪層中，0°鋪層產(chǎn)生了較大變形但只在支座附近出現(xiàn)輕微致密化現(xiàn)象；在LBX-2C-0-30中，胞元30°鋪層的致密化現(xiàn)象較0°鋪層更為嚴(yán)重；LBX-2C-0-30在壓頭下方出現(xiàn)對稱“V”形致密化區(qū)域，其他六邊形模型均未出現(xiàn)此現(xiàn)象。

內(nèi)凹形模型的實驗結(jié)果和有限元仿真結(jié)果如表5所示。由表5 可知，在內(nèi)凹形模型中，NAX-1C-15、NAX-1C-30裂紋擴展方向與豎直方向的夾角同基本構(gòu)型轉(zhuǎn)動的角度一致，分別為15°和30°；NAX-1C-45的裂紋擴展方向與豎直方向的夾角大約為20°。NAX-2C-0-15、NAX-2C-0-30和NAX-2C-0-45的非0°鋪層裂紋擴展方向相對于對應(yīng)單層模型的裂紋擴展方向向豎直方向偏轉(zhuǎn)。其余模型的裂紋擴展方向都為豎直方向。由于失穩(wěn)現(xiàn)象，NAX-2C-0-15上方萌生了新裂紋，且每層新裂紋的擴展方向與預(yù)置裂紋擴展方向一致，其他模型均未出現(xiàn)此現(xiàn)象。

在內(nèi)凹形單層模型中，NAX-1C-90產(chǎn)生了較大變形，在壓頭兩側(cè)與支座附近出現(xiàn)了大面積致密化區(qū)域；NAX-1C-30和NAX-1C-45產(chǎn)生了輕微變形，在靠近壓頭處出現(xiàn)了輕微致密化現(xiàn)象。在內(nèi)凹形雙層模型中，由于失穩(wěn)現(xiàn)象，NAX-2C-0-15和NAX-2C-0-90的0°鋪層出現(xiàn)了較大變形，但未出現(xiàn)致密化區(qū)域，非0°鋪層則在壓頭附近產(chǎn)生了嚴(yán)重的致密化現(xiàn)象。

在所有模型中，反對稱結(jié)構(gòu)裂紋擴展路徑大致為豎直方向，其余模型的非0°層相對于單層模型都發(fā)生了向豎直方向的偏轉(zhuǎn)。在雙層模型中，0°層總是先出現(xiàn)裂紋，由于層間黏結(jié)力的存在，非0°層受到0°層裂紋擴展力的影響，從而使其裂紋擴展方向向0°層裂紋擴展的方向發(fā)生偏轉(zhuǎn)。發(fā)生失穩(wěn)現(xiàn)象的模型，均是朝0°層偏轉(zhuǎn)，導(dǎo)致非0°層產(chǎn)生嚴(yán)重致密化現(xiàn)象。這與雙層模型兩層的強度差有關(guān)，強度低的一側(cè)會發(fā)生更大的變形，從而導(dǎo)致模型向強度高的一側(cè)偏轉(zhuǎn)。

3. 2 承載力和變形

選取壓頭下壓距離和下壓過程中的峰值力為判斷指標(biāo)，壓頭的下壓距離越大，則模型韌性越強；壓頭下壓過程中的峰值力越大，則模型的承載力越強；對于壓頭下壓至30 mm前產(chǎn)生嚴(yán)重失穩(wěn)的模型，選取發(fā)生失穩(wěn)時壓頭的下壓距離為判斷指標(biāo)，發(fā)生失穩(wěn)時壓頭的下壓距離越大，則結(jié)構(gòu)穩(wěn)定性越好。

在六邊形模型中，所有單層模型和LBX-2C-15-15均在壓頭下壓至30 mm之前發(fā)生了脆性斷裂，導(dǎo)致無法繼續(xù)加載。六邊形模型壓頭下壓距離和下壓過程中的峰值力如表6所示，載荷-位移曲線如圖7、圖8所示。

由表6、圖7和圖8可知，六邊形單層模型中，0°模型承載力和韌性最弱；30°模型承載力和韌性最強，相較于單層0°和15°，承載力分別提升了12. 2% 和0. 9%，韌性分別提升了33. 1%和20. 6%。六邊形雙層模型中，LBX-2C-0-15相較于單層0°和15°模型，承載力分別提升了21. 9%和9. 6%，韌性分別提升了236. 3%和204. 9%；LBX-2C-15-15的承載力和韌性最弱，承載力相較于單層15°模型下降了18. 3%，韌性相較于單層15°模型提升了9. 7%；LBX-2C-0-30 承載力和韌性最強，承載力相較于單層0°、30°，雙層15°/-15°和0°/15°模型分別提升了29. 2%、15. 2%、42. 2%和6%，韌性相較于單層0°、30°和雙層15°/-15°模型分別提升了236. 3%、204. 9%和178%。

在內(nèi)凹形模型中，NAX-2C-0-15 和NAX-2C-0-90在壓頭下壓至30 mm前發(fā)生了嚴(yán)重失穩(wěn)導(dǎo)致不能加載，單層90°模型和NAX-2C-45-45壓頭下壓至30 mm實驗停止，其余模型均在壓頭下壓至30 mm之前發(fā)生了脆性斷裂導(dǎo)致無法繼續(xù)加載。內(nèi)凹形模型壓頭的下壓距離和下壓過程中的峰值力如表7所示，載荷-位移曲線如圖9、圖10所示。

由表7、圖9可知，內(nèi)凹形單層模型中，0°模型的承載力和韌性最弱；90°模型的承載力和韌性遠(yuǎn)大于其他模型，相較于單層0°、15°、30°和45°，承載力分別提升了99. 8%、98. 8%、44. 7% 和64. 5%，韌性分別提升了523. 7%、274. 1%、72. 8%和96. 7%。在內(nèi)凹形雙層模型中，承載力最弱的模型為NAX-2C-15-15，韌性最弱的模型為NAX-2C-30-30；承載力最強的模型為NAX-2C-0-15，韌性最強的模型為NAX-2C-45-45。

由表7和圖10可知，內(nèi)凹形雙層模型中，NAX-2C-0-15相較于單層0°和15°，承載力分別提升了59. 2%和58. 4%，韌性分別提升了365. 7% 和179. 3%；NAX-2C-0-30 相較于單層0°，承載力提升了34. 2%，韌性提升了120. 6%；NAX-2C-0-30 和NAX-2C-30-30 相較于單層30°，承載力分別降低了2. 8% 和17. 3%，韌性分別降低了38. 9% 和53. 1%；NAX-2C-0-45相較于單層0°和45°，承載力分別提升了38. 9%和14. 4%，韌性分別提升了222% 和1. 5%；NAX-2C-45-45 相較于單層45°，承載力提升了26. 9%，韌性提升了96. 7%。NAX-2C-0-90相較于單層0°，承載力提升了48%，韌性提升了383. 6%；相較于單層90°，承載力降低了25. 9%，韌性降低了22. 5%。

3. 3 誤差分析

圖11～圖14所示分別為六邊形單層、六邊形雙層、六邊形0°/15°單雙層對比和六邊形0°/30°單雙層對比模型擬合載荷-位移曲線。本文有限元仿真僅考慮了彈塑性的本構(gòu)關(guān)系，未考慮材料和結(jié)構(gòu)的斷裂屬性和損傷規(guī)律。由圖11～圖14可知，彈性階段的實驗和仿真結(jié)果較吻合，且規(guī)律類似，彈性階段之后的仿真結(jié)果相較于實驗結(jié)果偏大（圖11～14中，T為實驗結(jié)果，S為仿真結(jié)果）。

4 結(jié)論

基于六邊形和內(nèi)凹形兩種基本構(gòu)型研究了單元方向、角度組合多孔結(jié)構(gòu)的破壞形式、承載力和變形的影響，經(jīng)實驗測試與有限元分析得到以下主要結(jié)論：

1）由于層間黏結(jié)力的存在，雙層模型的裂紋擴展路徑會向裂紋擴展抑制能力較差層的裂紋擴展路徑偏轉(zhuǎn)。

2）在雙層模型中，強度低的一側(cè)會發(fā)生更大的變形，從而導(dǎo)致模型向強度高的一側(cè)偏轉(zhuǎn)，發(fā)生失穩(wěn)現(xiàn)象，進(jìn)而會導(dǎo)致兩側(cè)模型發(fā)生不同程度的致密化現(xiàn)象。

3）六邊形單層模型中，胞元30°鋪層模型承載力和韌性最強；六邊形雙層模型中，0°與30°層合模型的承載力和韌性最強；內(nèi)凹形單層模型中，胞元90°鋪層模型的承載力和韌性最強；內(nèi)凹形雙層模型中，承載力最強的模型為0°與15°層合模型，韌性最強的模型為45°與-45°層合模型。

4）在一些角度組合下，雙層模型承載力和韌性相較于對應(yīng)的單層模型有明顯提升，六邊形0°與30°層合模型承載力和韌性提升幅度均最大，內(nèi)凹形0°與15°層合模型承載力提升幅度最大，內(nèi)凹形0°與90°層合模型韌性提升幅度最大。