改進(jìn)型FCC 晶格材料設(shè)計(jì)與吸能特性

郭 璐，劉志芳，李世強(qiáng)，吳桂英

（太原理工大學(xué)機(jī)械與運(yùn)載工程學(xué)院應(yīng)用力學(xué)研究所, 山西 太原 030024）

在航空航天、交通運(yùn)輸及軍事裝備等領(lǐng)域，蜂窩材料常被用于保護(hù)人員或結(jié)構(gòu)免受脈沖載荷的影響。晶格材料與蜂窩材料相似，具有超輕、高模量、高強(qiáng)度、延性好等特點(diǎn)，近年來得到了廣泛關(guān)注及研究[1-2]。與其他輕質(zhì)材料相比，晶格材料一般由精細(xì)、規(guī)則的單元結(jié)構(gòu)周期性排列而成，其力學(xué)性能高度依賴單元晶格的特性。光刻及增材制造技術(shù)可以實(shí)現(xiàn)對(duì)單元晶格結(jié)構(gòu)的小尺度幾何控制及精細(xì)制造[3]，從而獲得具有各種不同優(yōu)異力學(xué)性能的晶格材料。目前對(duì)晶格材料的研究大多集中于探索具有高比強(qiáng)度、高比剛度的晶格材料，并通過仿生、優(yōu)化等方法得到大量新型晶格結(jié)構(gòu)[4-6]。值得注意的是，晶格材料在能量吸收方面也具有巨大應(yīng)用潛力。

受金屬微觀體心立方（BCC）和面心立方（FCC）晶體結(jié)構(gòu)的啟發(fā)，本研究提出一種兼顧材料剛度和能量吸收性能的316L 不銹鋼FCC 晶格材料，并與常見的BCC 晶格材料進(jìn)行比較，利用 ABAQUS 有限元軟件對(duì)這兩種晶格材料進(jìn)行準(zhǔn)靜態(tài)和動(dòng)態(tài)壓縮數(shù)值模擬，分析兩種晶格材料在不同沖擊速度下的變形模式及應(yīng)力-應(yīng)變響應(yīng)，定量分析兩種晶格的能量吸收性能。

1 仿金屬材料晶體結(jié)構(gòu)的晶格

許多研究者對(duì)BCC 晶格材料（見圖1(a)）的彈塑性力學(xué)行為進(jìn)行了較全面的研究[14]，但對(duì)其能量吸收性能的探討較少?；诮饘傥⒂^FCC 晶體結(jié)構(gòu)，考慮到單元晶格中豎直桿的存在將使低密度周期性晶格材料極易發(fā)生彈性屈曲而導(dǎo)致晶格材料整體破壞，故與Pham 等[15]的研究不同，本研究設(shè)計(jì)了一種無豎直支撐桿件的FCC 晶格單元。與常見的Octet 結(jié)構(gòu)[16]相比，該FCC 晶格內(nèi)部未嵌套八面體單元晶格，在相同晶格大小及桿件尺寸下具有較低的相對(duì)密度。

圖1 兩種仿金屬晶體結(jié)構(gòu)的單元晶格Fig. 1 Unit cells of two lattices inspired by metal crystal structures

根據(jù)Deshpande 等[17]給出的Maxwell 準(zhǔn)則，BCC 晶格材料為彎曲主導(dǎo)型晶格，F(xiàn)CC 晶格材料為拉伸主導(dǎo)型晶格。

2 有限元模型及材料參數(shù)

2.1 有限元模型的建立

Cheng 等[18]的研究表明，當(dāng)立方體晶格材料在任意特征尺寸上陣列的單元晶格數(shù)大于或等于6 時(shí)，晶格材料的力學(xué)性能趨于穩(wěn)定。為此，將圖1 所示的兩種晶格分別周期性排列為6×6×6的立方體晶格材料，利用SolidWorks 建立三維有限元模型后，導(dǎo)入HyperMesh 進(jìn)行網(wǎng)格劃分及處理，最后利用ABAQUS 顯式算法對(duì)兩種晶格材料的單軸壓縮靜力學(xué)與動(dòng)力學(xué)響應(yīng)進(jìn)行數(shù)值模擬計(jì)算。圖2 顯示了有限元模型及相應(yīng)尺寸，晶格單元尺寸為5 mm×5 mm×5 mm。

圖2 有限元模型及模型尺寸Fig. 2 Finite element model and corresponding dimensions

為定量比較BCC 和FCC 晶格材料的能量吸收性能，取相同相對(duì)密度的BCC 和FCC 晶格材料進(jìn)行分析。假設(shè)晶格材料結(jié)點(diǎn)為完整球體，則單元晶格相對(duì)密度的理論模型可表示為

取相對(duì)密度為10.5%～10.6%的BCC 和FCC晶格材料進(jìn)行研究，即BCC 晶格材料的d/L=0.150，F(xiàn)CC 晶格材料的d/L=0.135。

Gümrük 等[14]通過數(shù)值模擬與實(shí)驗(yàn)對(duì)比發(fā)現(xiàn)，當(dāng)晶格材料的d/L＞0.087 時(shí)，細(xì)觀桿件的實(shí)體單元網(wǎng)格劃分較梁?jiǎn)卧獎(jiǎng)澐帜芨玫貜?fù)現(xiàn)實(shí)驗(yàn)結(jié)果，且能更好地描述晶格材料的后屈服行為。為此，單元晶格桿件采用實(shí)體單元C3D4 進(jìn)行劃分，通過網(wǎng)格敏感性分析，取網(wǎng)格尺寸為0.1 mm。如圖2 所示，在晶格材料頂部設(shè)置剛性板，通過剛性板的恒定速度加載來實(shí)現(xiàn)晶格的單軸壓縮。準(zhǔn)靜態(tài)壓縮時(shí)，剛性面板速度v取1 m/s，下面板保持固定；動(dòng)態(tài)壓縮時(shí)，剛性面板速度取10～100 m/s。

題目：圖書館門前建了兩個(gè)同樣大小的圓柱形花壇?；▔牡酌嬷睆綖?米，高為0.9米，往里裝泥土的高是0.7米，兩個(gè)花壇中共需要填土多少方？

圖3 相對(duì)密度理論模型與CAD 結(jié)果的比較Fig. 3 Comparison of the predictions of relative density between theoretical results and calculation with CAD models

2.2 材料參數(shù)

鑒于晶格材料制造工藝的特殊性，同時(shí)考慮到基體材料的選擇對(duì)晶格材料的變形模式等有較大影響[19]，金屬材料通常具有良好的韌性和較小的屈服應(yīng)變，為此采用選擇性激光熔化（selective laser melting，SLM）法制備的316L 不銹鋼金屬材料作為晶格基體材料。目前，已有大量關(guān)于SLM 制備316L 不銹鋼材料力學(xué)性能的研究。例如，Li 等[20]采用數(shù)值模擬與實(shí)驗(yàn)相結(jié)合的方法，研究了包括單軸拉伸、剪切等4 種不同工況下SLM 制備的316L 不銹鋼材料在準(zhǔn)靜態(tài)及高應(yīng)變率加載下的應(yīng)力-應(yīng)變響應(yīng)，并建立了材料的Johnson-Cook 模型

表1 SLM 制備316L 不銹鋼材料參數(shù)[20]Table 1 Material parameters of 316L stainless steel manufactured by SLM[20]

3 有限元方法驗(yàn)證

為評(píng)估ABAQUS 模擬晶格材料單軸壓縮響應(yīng)的準(zhǔn)確性，對(duì)模型中的相關(guān)能量進(jìn)行了驗(yàn)證，如圖4所示，其中：Et為總能，Ei為內(nèi)能，Ek為動(dòng)能，Ea為偽變形能。首先驗(yàn)證使用ABAQUS 顯式算法模擬材料準(zhǔn)靜態(tài)壓縮的合理性，如圖4(a)所示，當(dāng)上面板的速度v為1 m/s 時(shí)，BCC 和FCC 晶格材料主要通過晶格細(xì)觀桿件的彈塑性變形吸能，即內(nèi)能為材料吸能的主要形式，動(dòng)能的占比很小，且均小于內(nèi)能的5%，可忽略不計(jì)?？梢哉J(rèn)為，當(dāng)v=1 m/s 時(shí)，材料的壓縮過程是準(zhǔn)靜態(tài)的。

動(dòng)態(tài)加載下，應(yīng)保證系統(tǒng)的偽變形能與系統(tǒng)內(nèi)能之比小于5%[22]。以v=75 m/s 時(shí)的FCC 晶格為例，圖4(c)給出了壓縮過程中的能量平衡情況，其偽應(yīng)變能與系統(tǒng)內(nèi)能的比值遠(yuǎn)遠(yuǎn)小于極限值5%。可以認(rèn)為，模型是有效的。

圖4 相對(duì)密度為10.5% 的BCC 和FCC 晶格材料的能量驗(yàn)證Fig. 4 Energy verification of BCC and FCC lattice materials in case of relative density 10.5%

4 數(shù)值模擬結(jié)果及分析

4.1 晶格材料的初始彈性模量

研究表明，晶格材料的數(shù)值模擬計(jì)算能較精準(zhǔn)地預(yù)測(cè)材料的線彈性力學(xué)響應(yīng)[14,21]，故本研究采用數(shù)值模擬初始加載曲線的斜率定義材料剛度。計(jì)算得到BCC 和FCC 晶格材料的初始彈性模量分別為0.388 和2.576 GPa，即對(duì)相對(duì)密度為10.5%～10.6%的晶格材料而言，F(xiàn)CC 晶格材料的初始彈性模量可達(dá)BCC 晶格材料的6.64 倍，這是由于拉伸主導(dǎo)型晶格的比剛度普遍比彎曲主導(dǎo)型晶格高[17]。

4.2 晶格材料的準(zhǔn)靜態(tài)壓縮

4.2.1 晶格材料準(zhǔn)靜態(tài)壓縮變形模式

相對(duì)密度為10.5%～10.6%的BCC 和FCC 晶格材料在單軸壓縮載荷作用下的變形模式如圖5 所示。

圖5 BCC 和FCC 晶格材料的準(zhǔn)靜態(tài)壓縮變形模式Fig. 5 Quasi-static deformation modes for BCC and FCC lattice materials

從圖5 可以看出，BCC 晶格材料發(fā)生均勻變形，而FCC 晶格材料則發(fā)生45°剪切變形。分析認(rèn)為，F(xiàn)CC 晶格材料由拉伸主導(dǎo)型晶格材料周期排列而成，各細(xì)觀桿件主要承受軸力作用，小變形下易導(dǎo)致材料局部失穩(wěn)。

4.2.2 晶格材料準(zhǔn)靜態(tài)壓縮結(jié)果分析

對(duì)BCC 及FCC 晶格材料的準(zhǔn)靜態(tài)能量吸收能力進(jìn)行了分析，兩種晶格材料的準(zhǔn)靜態(tài)單軸壓縮應(yīng)力-應(yīng)變響應(yīng)如圖6 所示。與傳統(tǒng)蜂窩材料類似，定義材料的能量吸收效率f[2]

圖6 準(zhǔn)靜態(tài)壓縮晶格材料的應(yīng)力-應(yīng)變響應(yīng)Fig. 6 Stress-strain response curves of lattice materials under quasi-static compression

4.3 晶格材料的動(dòng)態(tài)壓縮

4.3.1 晶格材料的動(dòng)態(tài)壓縮變形模式

隨著壓縮速度的增大，晶格材料的變形模式會(huì)發(fā)生改變。圖7 顯示了BCC 和FCC 兩種晶格材料在v=75 m/s 的動(dòng)態(tài)壓縮下的變形模式，此時(shí)晶格材料的應(yīng)變率為2 500 s-1。在該應(yīng)變率下，BCC 和FCC 晶格材料靠近受壓側(cè)單元晶格產(chǎn)生較大的橫向形變，晶格材料整體發(fā)生逐層壓潰。與此同時(shí)，F(xiàn)CC 晶格材料在變形過程中還伴隨著剪切變形的發(fā)生，如圖7(b)所示。

圖7 v=75 m/s 時(shí)晶格材料的變形模式Fig. 7 Deformation modes of lattice materials at v =75 m/s

4.3.2 晶格材料動(dòng)態(tài)壓縮結(jié)果分析

晶格材料承受沖擊載荷作用時(shí)，隨著沖擊速度的提高，動(dòng)能占總能量的比例呈上升趨勢(shì)。以v=75 m/s為例，在FCC 晶格材料被壓實(shí)之前約12%的能量存儲(chǔ)為動(dòng)能。由圖8 可知，在1～100 m/s 速度范圍內(nèi)，F(xiàn)CC 晶格材料均可吸收更多的能量，該現(xiàn)象歸因于晶格材料內(nèi)能的增加而非動(dòng)能的增加。同時(shí)，晶格材料的塑性耗散能在內(nèi)能中的占比始終大于85%。

圖8 相同相對(duì)密度下不同晶格材料在1～100 m/s 速度范圍內(nèi)的能量吸收特性對(duì)比Fig. 8 Comparison of energy absorption for different lattice with the same relative density at various velocities of 1-100 m/s

不同沖擊速度下BCC 和FCC 晶格材料的動(dòng)態(tài)響應(yīng)如圖9 所示。顯然，當(dāng)沖擊速度為75 m/s 時(shí)，BCC 和FCC 晶格材料的應(yīng)力-應(yīng)變曲線均出現(xiàn)波動(dòng)，且FCC 晶格的應(yīng)力峰個(gè)數(shù)與晶格層數(shù)相同。上述數(shù)值模擬結(jié)果表明，當(dāng)壓縮應(yīng)變率大于或等于2 500 s-1時(shí)，SLM 制備的316L 不銹鋼BCC 和FCC 晶格材料的變形模式及平臺(tái)應(yīng)力均具有一定的應(yīng)變率敏感性。

圖9 不同沖擊速度下晶格材料的應(yīng)力-應(yīng)變響應(yīng)Fig. 9 Stress-strain curves of lattice materials subjected to impact at different velocities

高速?zèng)_擊下材料會(huì)產(chǎn)生應(yīng)力-應(yīng)變波動(dòng)，用承載能力波動(dòng)U（ULC）來定量分析壓縮應(yīng)力在平臺(tái)應(yīng)力周圍的波動(dòng)程度[23]

綜合上述兩種晶格在不同沖擊速度下的變形模式、動(dòng)能比例及壓縮力效率，定義相對(duì)低速（0～15 m/s）、中速（15～75 m/s）和高速（大于75 m/s）3 種沖擊速度模式。例如：低速模式下，BCC 晶格材料的動(dòng)能比例極小，如圖8(a)所示，其壓縮力效率在70%～80%之間；中速模式下，BCC 晶格材料的壓縮力效率在40%～70%之間，動(dòng)能所占比例較低；高速模式下，BCC 晶格材料的變形模式有明顯改變，如圖8(b)所示，晶格材料的動(dòng)能占比明顯增大，且壓縮力效率整體上小于40%。

此外，低速模式下FCC 晶格材料的壓縮力效率最高，可達(dá)90%。從低速模式向中速模式過渡時(shí)，壓縮力效率略有下降，在60%～80%區(qū)間，即隨著沖擊速度的提高，晶格材料應(yīng)力-應(yīng)變響應(yīng)中峰值應(yīng)力的出現(xiàn)導(dǎo)致壓縮力效率緩慢下降，如圖9 所示。當(dāng)沖擊速度大于75 m/s 時(shí)，晶格材料進(jìn)入高速模式，壓縮力效率有較大的下降。該現(xiàn)象可用承載能力波動(dòng)解釋，表現(xiàn)出與壓縮力效率相反的趨勢(shì)，如圖10(a)所示。對(duì)BCC 晶格材料而言， 100 m/s 高速?zèng)_擊時(shí)壓縮力效率的增大僅為一個(gè)數(shù)值波動(dòng)而不影響壓縮力效率的整體下降趨勢(shì)，如圖10(b)所示，且當(dāng)沖擊速度為150 和200 m/s 時(shí)，對(duì)應(yīng)的承載能力波動(dòng)分別為0.275 和0.327，壓縮力效率分別為30.534%和24.732%，即兩種晶格材料的壓縮力效率大體上隨沖擊速度的提高呈下降趨勢(shì)。

圖10 承載能力波動(dòng)及壓縮力效率與沖擊速度的關(guān)系Fig. 10 Trends of undulation of load-carrying capacity and CFE versus the impact velocity

4.3.3 晶格材料動(dòng)態(tài)壓縮的沖擊波模型

圖11 密實(shí)化應(yīng)變(a)、平臺(tái)應(yīng)力(b)以及塑性能量耗散(c)與沖擊速度的關(guān)系Fig. 11 Trends of onset strain of densification (a), plateau stress (b), and plastic energy dissipation (c) versus the impact velocity

5 結(jié)果討論

5.1 晶格材料的歸一化彈性模量

負(fù)泊松比材料具有較好的能量吸收性能，但比剛度較低[26-33]，而桁架晶格則能夠?qū)崿F(xiàn)較高的比強(qiáng)度和比剛度，其中材料的比剛度被定義為歸一化彈性模量與相對(duì)密度的比值。圖12 對(duì)比了BCC 和FCC 兩種晶格材料與常見負(fù)泊松比材料及桁架晶格的歸一化彈性模量-相對(duì)密度（E*/Es- ρ*/ ρs）關(guān)系。顯然，在相同相對(duì)密度下，F(xiàn)CC 晶格材料的歸一化彈性模量大于大多數(shù)材料，即較其他大多數(shù)材料而言，F(xiàn)CC 晶格材料具有更高的比剛度。

5.2 晶格材料的能量吸收性能評(píng)估

由4.2.2 節(jié)可知，相同密度的FCC 晶格材料比BCC 晶格材料具有更高的歸一化比吸能，其歸一化比吸能可達(dá)BCC 晶格材料的2.6 倍。由圖12 可知，拉伸主導(dǎo)型FCC 晶格比彎曲主導(dǎo)型BCC 晶格具有更高的比剛度，即材料能夠在較小的應(yīng)變下達(dá)到更高的強(qiáng)度，而且在準(zhǔn)靜態(tài)壓縮載荷作用下316L 不銹鋼FCC 晶格材料雖然發(fā)生局部失穩(wěn)，但是未發(fā)生破壞，且有一段較平滑的應(yīng)力平臺(tái)階段（見圖6(b)和圖9(b)），此外，相同相對(duì)密度的兩種晶格材料具有大致相同的密實(shí)化應(yīng)變，使得FCC晶格材料具有比BCC 晶格材料更高的比吸能。

圖12 BCC 及FCC 晶格材料與其他材料的歸一化彈性模量-相對(duì)密度關(guān)系比較Fig. 12 Comparison of normalized Young’s modulus versus relative density for BCC and FCC lattice materials with other materials

將BCC 晶格材料的吸能效率及壓縮力效率與FCC 晶格材料及傳統(tǒng)晶格材料[34-37]進(jìn)行比較，可以看出，相對(duì)密度為10.5%的FCC 晶格材料的壓縮力效率可達(dá)90%（見圖13(b)），明顯高于其他傳統(tǒng)晶格材料（能量吸收效率為59%），也高于其他相對(duì)密度為10%～16% 的傳統(tǒng)晶格材料。FCC 晶格材料的歸一化比吸能為0.228，高于其余幾種輕質(zhì)吸能材料（見圖13(c)）。

圖13 BCC 及FCC 晶格材料與其他材料的吸能效率(a)、壓縮力效率(b)和歸一化比吸能UM,n (c)比較Fig. 13 Comparison of energy absorption efficiency (a), compression force efficiency (b) and specific energy absorption UM,n (c) for BCC and FCC lattice materials with other materials

6 結(jié) 論

基于仿金屬微觀晶體結(jié)構(gòu)構(gòu)造了FCC 單元晶格，研究了BCC 和FCC 兩種晶格材料的單軸壓縮力學(xué)行為，選擇SLM 法制備的316L 不銹鋼材料為基體材料，對(duì)兩種晶格材料的抗沖擊性能進(jìn)行了數(shù)值模擬分析，定量分析了兩種晶格材料的能量吸收性能，給出了單軸動(dòng)態(tài)加載下晶格材料的壓縮平臺(tái)應(yīng)力及塑性能量耗散的半經(jīng)驗(yàn)公式，得到如下主要結(jié)論：

(1) 相同相對(duì)密度的FCC 晶格材料與常見的負(fù)泊松比材料及桁架晶格材料相比，具有較高的比剛度；

(2) 與BCC 晶格材料相比，相同相對(duì)密度的FCC 晶格材料在準(zhǔn)靜態(tài)壓縮載荷作用下的能量吸收性能更優(yōu)異，相對(duì)密度為10.5%～10.6%的FCC 晶格材料的比吸能是BCC 晶格材料的2.6 倍；

(3) 當(dāng)應(yīng)變率為2 500 s-1時(shí)，兩種晶格材料均表現(xiàn)出一定程度的應(yīng)變率敏感性，其變形模式及應(yīng)力-應(yīng)變響應(yīng)均有明顯改變，并且隨著沖擊速度的提高，晶格材料的密實(shí)化應(yīng)變及平臺(tái)應(yīng)力呈上升趨勢(shì)，故能量吸收效率增大，而壓縮力效率則由于應(yīng)力波動(dòng)的增加呈現(xiàn)下降趨勢(shì)；

(4) FCC 晶格材料既具有較高的比剛度，又具有較高的能量吸收效率、壓縮力效率及歸一化比吸能特性，有望用于夾芯梁及夾芯板等輕質(zhì)結(jié)構(gòu)中。