作為快速加工和制造的重要技術(shù),3D打印可以借助計算機設(shè)計生產(chǎn)具有復(fù)雜幾何形狀的部件。特別地,使用熔融沉積成型 (FDM) 增材制造技術(shù)的連續(xù)纖維增強熱塑性復(fù)合材料 (CFRTPC) 受到廣泛關(guān)注和研究。本文對CFRTPCs進行了有限元分析,模擬其打印-冷卻-沉積過程,并獲得其力學(xué)性能。首先,發(fā)展增材制造過程的理論模型;其次,對3D打印CFRTPC孔洞含量的影響因素進行參數(shù)化分析,揭示了堆疊方式、走線長寬比、走線間寬度及高度對孔洞含量的影響規(guī)律;最后,采用宏細觀模擬方法,對3D打印的CFRTPCs進行了多尺度力學(xué)分析。本文對3D打印CFRTPCs的研究提供了模擬方法和理論依據(jù)。CFRTPCs的3D打印過程如圖1所示。熱塑性樹脂長絲進入擠出頭,并在噴嘴被加熱器加熱至熔融。同時,連續(xù)的纖維束從纖維供應(yīng)頭被傳送至噴嘴。在噴嘴中,連續(xù)的纖維束被熔融熱塑性聚合物滲透并涂覆。隨后,浸漬后的纖維束可以從噴嘴的下端被擠出。當(dāng)擠出的材料到達沉積平臺并迅速冷卻結(jié)晶后,連續(xù)的纖維由于前端的熱塑性樹脂凝固而連續(xù)地被拉出。按照設(shè)計軌跡,噴嘴可以沿著平面方向發(fā)生移動,這樣形成了3D復(fù)雜形狀的CFRTPCs構(gòu)件的第一層。圖1 CFRTPCs 3D打印過程示意圖
在單層完成后,噴嘴會沿著構(gòu)件厚度方向提升單個層厚度。擠出的材料到達零件表面,并會迅速固化并粘附到上一層。重復(fù)該過程,逐層完成構(gòu)件的沉積,直到完成其打印。通過上述3D打印方法,CFRTPC可以通過沉積熔融制成。上述3D打印成型過程包括3部分,分別為熱塑性樹脂浸潤纖維過程、溫度變化及熱傳導(dǎo)過程、熱塑性樹脂固化過程。圖2給出了多個典型階段的樹脂流動情況及對應(yīng)的橫截面3D打印珠粒形貌。圖 2 3D打印模擬中樹脂的流動情況及對應(yīng)的珠粒形貌:(a) 單珠粒打印;(b) 打印頭沿X軸正向轉(zhuǎn)彎;(c) 兩排珠粒的凝結(jié);(d) 第二層珠粒打?。?e) 兩層珠粒的凝結(jié);(f) 3D打印完整試件模擬
流體受到重力作用,同時它會由于固化產(chǎn)生體積收縮、密度增大和粘度增大,這些均導(dǎo)致其凝固之后的珠粒形貌與給定的橢圓形形貌不相同。珠粒的寬度在沉積-冷卻-凝固過程將變窄5.02%,其高度將提升1.59%,其整體的截面積將降低3.51%。隨后,打印頭沿著X軸正方向發(fā)生轉(zhuǎn)彎,并隨后轉(zhuǎn)向Z軸正方向繼續(xù)打印,新流入的高溫熔融樹脂將會令周圍的樹脂融化,并重新冷卻凝固,從而兩條并排的珠粒將發(fā)生粘接,其粘接的情況則與兩珠粒間距 相關(guān)。在一定程度內(nèi), 3D打印試件越致密,其力學(xué)性能將會越好。當(dāng)試件的第一層打印結(jié)束后,3D打印噴頭將沿著Y方向上升,并將第一層的終點作為第二層的起點繼續(xù)打印。上層的高溫熔融樹脂將融化掉下層已經(jīng)凝固的樹脂,兩層樹脂由于重新冷卻凝固而發(fā)生粘接,受到重力的影響,下層樹脂的形狀略有變化,其高度將略降低,并且致密性變好。最終CFRTPCs試件按照設(shè)計的路徑逐層完成3D打印。在打印過程中,可以看到在兩個珠粒之間存在樹脂未曾浸潤的區(qū)域,在本模型中,未曾浸潤區(qū)域的孔洞為12.23%,試件整體的纖維體積分數(shù)為26.33%。不同參數(shù)下得到的CFRTPCs試件的截面如圖3所示。發(fā)現(xiàn)堆疊方式及橢圓形珠粒的長寬比對孔洞含量的影響最為顯著。圖3 影響3D打印CFRTPCs孔洞含量的因素分析
為了能夠?qū)崿F(xiàn)對整體CFRTPCs試件的力學(xué)分析,需要建立隨機纖維代表性體積單元,如圖4所示?;谝阎腁BS熱塑性樹脂和T300碳纖維的力學(xué)性能,獲得樹脂與纖維的混合物的平均物理參量。圖4 碳纖維體積分數(shù)為30%的纖維隨機分布的代表性體積單元
對3D打印得到的CFRTPCs進行力學(xué)分析,需要將獲得的CFRTPCs試件的形貌輸出并導(dǎo)入固體分析軟件中。對試件的兩端設(shè)置邊界條件,令其一端固定,另一端沿著Z軸正方向以2mm/min的速度拉伸。本文描述的打印流程和材料參數(shù)與文獻中實驗中的參數(shù)相同,將兩種方法得到的應(yīng)力應(yīng)變曲線對比,結(jié)果如下:實驗中得到的3D打印的CFRTPCs的模量和強度分別為62.5GPa和986MPa,采用本文有限元方法獲得的試件模量與強度分別為64.72GPa和1001.17MPa,其對應(yīng)的誤差分別為3.55%和1.54%。結(jié)果表明本文提出的模擬手段可以對3D打印CFRTPCs的模量和強度給出很好的預(yù)測。將CFRTPCs中心位置的主承力區(qū)域放大,臨近失效時沿著纖維方向及垂直于纖維方向的損傷分布由圖5給出。試件的損傷首先發(fā)生在兩個珠粒之間的區(qū)域。這是由于在3D打印規(guī)則的橢圓形過程中,在孔洞缺陷附近將發(fā)生應(yīng)力集中,造成試件首先在該區(qū)域發(fā)生損傷。隨后,若CFRTPCs繼續(xù)受力,由于復(fù)合材料不同方向的強度特征,試件會沿著纖維方向開裂。圖5 宏觀CFRTPCs的損傷分布
該論文以Multi-scale analysis for 3D printed continuous fiber reinforced thermoplastic composites為題發(fā)表在《Composites Science and Technology》上。在清華大學(xué)姚學(xué)鋒教授的指導(dǎo)下,重慶大學(xué)弘深青年教師付宇彤為本文的第一作者。通訊作者清華大學(xué)姚學(xué)鋒教授目前主要從事先進復(fù)合材料結(jié)構(gòu)、航空橡膠密封結(jié)構(gòu)、飛機增升裝置、高通量設(shè)計、實驗力學(xué)與無損檢測等研究領(lǐng)域的科學(xué)研究工作。參與大型飛機C919、C929及支線客機ARJ21的重大科研攻關(guān)項目。兼任中國復(fù)合材料學(xué)會常務(wù)理事、中國力學(xué)學(xué)會實驗力學(xué)專業(yè)委員會委員、《Composite Structures》《Composites Part C》《 力學(xué)季刊 》《 現(xiàn)代機械 》等編委。