基于激光烧结纳米复材的高效吸能超材料

一、论文作者

• 袁上钦 西北工业大学无人系统技术研究院，西安 710072
• 王志豪 西北工业大学无人系统技术研究院，西安 710072

• 李  江  西北工业大学机电学院，西安 710072

※基金项目：国家自然科学基金 (联合基金U1930207，青年基金51905439）

本文选用激光烧结技术，采用纳米复合材料作为原料的，通过该技术制造具有拉胀复合晶格的机械超材料。该方法为设计和制造能量吸收可调控的系统提供了技术手段。通过对材料组分和结构的优化能够高效提升复材结构的能量吸收能力。拉胀型超材料的能量吸收能力与相对密度成指数关系，约为2.5^-3。经过对机械超材料进行合理的几何拓扑优化后，其具有较高的致密化强度（0.0195 MPa kg-1 m^-3）、能量吸收能力（6.29 MJ m^-3）和比能量吸收能力（20.42 J g^-1）。本研究提出的拉胀性超材料具有接近钛合金泡沫，优于塑料泡沫、铝合金泡沫以及其他增材制造成型结构的能量吸收能力。

关键词：曲屈机械超材料，增材制造，负泊松比晶格，能量吸收，选择性激光烧结

三、引言

近年来，研究者发现三维机械超材料具有天然材料不具备的机械性能，例如负泊松比，刚度和热膨胀性。尽管大多数已开发的超材料都可以用于工程项目，但由于增材制造使用的材料机械性能不足，以及制造过程中存在不可避免的结构缺陷，很难达到高效的吸收能量。超材料[1–5]是理想的人工材料，一般具有独特的结构，可以用于控制光[6]、声音[7,8]、机械应力[9,10]等，从而使其具有自然材料无法实现的特性，比如：光学性质[11]、声学性质[6]、机械性[12,13]和许多其他物理特性。研究者在金属晶体的变形过程中，首次观察到具有高度有序堆积的三维拉胀材料[14,15]。

迄今为止，机械超材料促进了原理上不相容的机械性能组合，例如刚度、阻尼能力和负泊松比[16–18]。这些属性均来自材料与结构的几何构型排列，以及材料组分的固有特性。但是，通过传统的材料开发（例如化学合成和自组装）和传统的制造方法（例如模制，加工和切割）很难实现材料配方的定制以及跨领域的三维结构设计。

三维机械超材料具有的的高能量吸收能力、可重复使用的机械减震能力、可控制的机械行为和可重构性等性能已被研究者预测[12,18,19]。Babaee等人设计了一组通过Bucklicrystals产生的机械超材料，并从理论上揭示了它们的屈曲引起的膨胀[12]。机械超材料与传统的能量吸收材料不同，后者是由空隙和组成材料的随机空间排列导致的固体泡沫组成的，它会导致局部结构异质性和不均匀性。金属或塑料的变形行为通常是不可控的，因此很难实现定向的拓扑设计。但是，实际加工的超材料很少具有较高的能量吸收能力。面向终端用途，关键在于超材料的设计和制造一体化，特别是材料配方和空间结构的理想组合。

3D打印或增材制造（AM）为合理设计和控制周期性变化的结构提供了机会，而结构的拓扑可以帮助实现其更好的性能，例如负刚度和负泊松比。例如，Fang和同事展示的轻型机械超材料具有可调节的负热膨胀[10]以及超高刚度[13]的性能，这要归功于先进的微立体光刻技术以及三维材料数字化功能。Frenzel等人还设计并制作了三维手征微结构，其扭曲自由度在数量上超过了Cauchy弹性。[9]

通过选择性激光烧结（SLS）进行的3D打印已被用于制造陶瓷、金属和聚合物复合材料，其尺度范围从微米到几米不等[18,20]。这种方法的本质是设计粉末的成分和微观结构，使其具有理想的流动能力并与激光源相互作用，最终冷凝成为块状复合材料[21,22]。激光烧结的拉胀超材料可以满足以下吸收能量的要求：（1）强韧；（2）结构中具有的柔性节点或连接点；（3）压缩后发生大变形；（4）制造大型结构。同时，与其他需要额外支撑结构、布局以及后处理十分复杂，与选择激光熔化（SLM），熔融沉积建模（FDM）的AM技术相比，SLS具有无支撑材料、易于后处理、成本低以及可以成型复杂结构的优点。

在航空制造领域，复合材料及结构的稳定性和可靠性尤其重要。同其他航空材料相比，复合材料及结构有更高的强度、更轻的重量、更好的减振以及降噪能力，这为其在航空领域的应用打下了良好基础。目前复合材料已经开始在客机上应用，空客A350的用量接近总质量的53%，波音787也超过了50%。为了满足航空应用的需求、获得适航批准的构件，需要使用复合材料制作部分零件，如：飞机前机身段、阻力板、机翼外翼、整流壁板，航空发动机的发动机的压气机叶片、盘、轴、机匣、传动杆等部件、机舱门锁拴轴和机翼肋条等，都可以使用复合材料制作。其中，负泊松比材料具有较高的抗冲击能力和减震性能，可以制成大型薄壁件的填充结构、机匣及机舱门等构件的内部支撑结构。材料的结构设计和组成体系需要高度耦合，同时增材制造复合材料为机械超材料在航空领域的应用提供了新的思路。因此，本研究开发了尼龙与纳米碳管的复合粉末，设计了不同孔数的体心立方结构（BCC）的拉胀晶格，通过SLS工艺制造负泊松比晶格结构，拓展增材制造机械超材料在航空领域的应用。

四、试验及方法

复合粉末制备方法：通过表面活性剂处理碳纳米管（CNT），通过溶液降温沉积法将CNT附着到尼龙（PA12）粉末表面，形成复合CNT/PA12粉末。PA12粉末（德国Krailling EOS GmbH，PA2200），直径为20–40nm或40–60nm。多壁CNT溶液长度大于5um（杭州Corker复合材料有限公司）。将水合胆酸钠（BoXtra，≥99％，Sigma-Aldrich）用作表面活性剂，以修饰去离子水中的多壁碳纳米管（MWCNT）。将CNT在90°C的加热的PA12粉末悬浮液中混合1小时，然后冷却至室温。吸附有CNT的粉末可能会沉淀出来，以进行进一步过滤和干燥，通过调节MWCNT的重量比来改变复合粉末的组成。

材料表征：为了表征复合粉末的表面形貌和微观结构，使用了场发射扫描电子显微镜（JSM-7600F，日本东京），电子束电压为2–5kV。为了使用光学显微镜（Olympus DP72，MA，0 2453，美国）捕获微观结构，需要对CNT/PA12样品进行良好的抛光以进行光学表征。

选择性激光烧结：激光烧结复合材料和结构均通过配备CO2激光器（λ=10.6µm，激光功率高达50W，德国EOSP395）。考虑了包括激光扫描速度、粉末层厚和阴影空间在内的参数，以优化烧结复合材料的机械性能。优化的工艺参数之后是EOS P395系统的商用PA12参数集。在将处理室冷却至室温之后处理完成。从机器中取出样品，并通过喷砂清理多余的粉末，并用压缩空气清理。拉伸和压缩样品分别按照ASTM标准（D638 14）和（D695 15）进行打印。抛光了15mm×20mm×25mm的立方试样，以检查其微观结构。

机械测试及评估：烧结复合材料的拉伸和压缩性能是使用Instron 3360负载框架（Instron，Norwood，MA）测量的。施加在样品上的拉伸速率和压缩速率分别为和 。每种材料UT的比拉伸韧性描述为等式

（1）

其中σ（ε）是与应变相对应的拉伸应力，而εd是断裂伸长率；ρs是激光烧结复合材料的密度。

在具有10 kN称重传感器的Instron 3343 Universal试验机上进行了拉胀晶格的压缩试验，并且以的相同位移速率固定了所施加的法向应变。样品的拉伸性能是平均值。每个晶格W的能量吸收能力由公式给出

（2）

其中εD是每个晶格在压缩时的致密化应变，并在具有最大吸收效率的应变下确定，该结构由方程式Emax给出

（3）

当达到Emax时，可以获得σD和εD。此后，如公式所述，通过比能量吸收Umax得出的单位质量能量吸收的最优值                             

（4）

其中Wmax是通过能量吸收能力压缩晶格时的最大能量吸收，而ρ是晶格密度。

五、结果与讨论

高性能纳米复合材料与理想的三维结构的匹配将使人们既可以利用拉胀型超材料的结构优势，又充分利用其能量吸收的潜力。本文对碳纳米管增强的纳米复合粉末进行激光烧结，以制备具有较高能量吸收能力的机械超材料，该材料超越了传统的金属随机泡沫或先前报道的增材制造复合材料。如图1所示，机械超材料的潜在优势需要通过增材制造工艺来实现。

由于CNT的网络使连续基体韧性增加，使用纳米复合材料制作的样品表现出了良好的拉伸韧性和比能量吸收能力。结构可产生屈曲主要的大形变和高度顺从效应，非常有利于能量的吸收。通常情况下，当尽可能多的结构单元避免与荷载方向一致，同时保证足够的结构失稳和力耗散时，结构的耗能效率最高。这些拉胀晶格具有各向同性，并具有高度的结构连通性，在静态压缩时表现出负泊松比和较大的几何非线性。

聚酰胺（PA12）是当前方法的良好实验原料，这是因为PA11和PA12占据激光烧结热塑性塑料90％的市场份额。激光烧结PA12具有约45MPa的拉伸强度和28％的断裂伸长率[23]。使用乳液沉积法，将平均粒径为60µm的原料粉末（PA12）附着表面活化处理的CNTs，以确保CNT在表面上均匀覆盖[24]。激光烧结工艺的系统粉体评价和工艺优化在我们之前的工作中得到了证明[21,23]。碳纳米管包覆的粉末可以有效吸收激光能量以促进其熔融凝结以及在聚合物基体内形成三维连续网络。但是，由于表面CNT膜的厚度与熔融复合材料的粘度之间存在正相关关系，覆膜加厚引起的熔体粘度增加可能会在短时间的激光照射下引起熔融障碍，从而阻止粉末的固结。所以，当涂层既能有效地增强能量吸收又能保持足够低的熔体粘度以保证复合材料熔体的快速熔融时，可以获得最佳的负载比[25]。

• a）用于选择性激光烧结（SLS）工艺的典型建筑平台的示意图。 
• b）具有均匀CNT覆膜的粉末。 
• c）烧结的CNT/PA12复合材料中的CNT网络。 
• d，e）不同放大倍数的激光烧结复合材料的微观结构。 
• f）具有BCC-6H结构的激光烧结复合超材料。
• g）BCC-6H晶格的代表性单元。 
• h）应力-应变曲线说明了BCC-6H晶格的能量吸收性能。机械变形的临界转变（例如平均应力，致密化应变和致密化应力）在应变-应力曲线中显示。

Figure1. Microstructure of CNT and energy absorption capacityof crystal lattice：a) Schematic illustration of a typical building platform forselective laser sintering (SLS) process. b) The feedstock powder with a layerof CNT coating. c) The networks of CNT s within sintered CNT/PA12 composites.d,e) The microstructures of laser-sintered composite with differentmagnifications. f) The laser-sintered composite metamaterials with BCC-6Hpattern. g) The representative cell unit of the BCC-6H lattice. h) Theenergy-absorption performance of the BCC-6H lattice illustrated by thestress–strain curve. The critical transitions of mechanical deformation such asaverage plateau stress, densification strain, and densification stress areindicated in the strain–stress curve.

为了获得最佳的激光烧结CNT/PA12复合材料的力学性能，在SLS系统中应用了多种材料组合。将直径范围为20–40或40–60nm的CNT作为原料涂覆在PA12粉末上作为进料材料。在图2a中，最佳的CNT/PA12复合材料可以明显改善机械性能，其中：杨氏模量E=1301±14.5 MPa，极限抗拉强度σu= 68±2.6 MPa，断裂伸长率εf= 33±4％，并且这些特性几乎与粉末的方向无关。

为了进一步研究各向同性的强化和增韧机理，对烧结复合材料及其变形成分进行了微观结构表征。观察结果表明，CNT网络的均匀分布可以阻止分离的聚合物相的移动，从而增强基质（图2b，c）。另一方面，变形复合材料的微观结构表明，外力在微观尺度上引起聚合物相的大变形。所得的CNT网络被密集地压实，并保留在这些分段的边界之间（图2d，e）。

这种精细压实的微观结构提供了一种适应应变的方法，从而防止了裂纹的产生和扩展，使基体坚实[26,27]。但是，CNTs（20–40 nm）的负载增加到1wt％会对聚合物基体的机械性能产生不利影响（图2a）。添加0.5wt％的CNT（40–60 nm）也会导致断裂伸长率从28％下降至24％。这主要是因为增加的涂层会提高复合材料的熔融粘度，然后阻止熔融粉末的固结，不充分的融合会引起巨大的微观缺陷或空洞，导致烧结复合材料的脆性破坏及其延展性的缺陷。因此，优化复合材料配方和有效控制烧结过程是在不牺牲其杨氏模量的情况下增强CNT纳米复合材料强度和韧性的关键步骤。

优化的CNT增强复合材料的机械性能（例如：韧性、比能吸收率和延展性）不仅优于纯聚合物，而且还超过了许多广泛使用的3D打印材料，包括塑料、复合材料和合金（图2f）。在SLS工艺中采用CNT增强粉末，对工程化所需的组成和微观结构起了关键作用，从而使断裂伸长率提高了33％，比吸收率达到20.68 J g-1。相比之下，通过FDM或直接添加的聚合物复合材料通常在聚合物相和微纤维之间存在界面缺陷和不可避免的空洞[28–32]。铝和钛合金通过SLM和电子束熔化（EBM）均表现出断裂伸长率降低的特征，这与快速熔化固结过程中形成的缺陷和残余应力有关[25]。

此外，通过激光烧结纳米复合材料制造可拉胀的超材料可能是一种有效的策略，可以实现将有益的机械性能植根于恒定的材料和结构。拉胀型超材料由一系列球形壳组成，在空间上排列以形成晶体结构。以体心立方模式组装的空心球形块由于高度的结构不稳定性而具有很高的可压缩性和变形性，可以通过更改壳体厚度、孔的数量和角度来进行拓扑优化。了解几何形状、机械性能和变形机制之间的关系至关重要，以便于计算出拉胀晶格允许的压应力和能量吸收能力。

为确定三维结构对能量吸收性能的影响，选择6孔（6H）和12孔（12H）构件来组装BCC模式，并调整壳厚度以更改绝对密度和相对密度。晶格分别为100-301 kg m-3和0.09-0.3（图3a）。该密度范围与气凝胶、氧化铝纳米晶格和其他超轻材料的密度范围相当[16,17,33]。压应力-应变曲线显示了准静态压缩时的变形和能量吸收（图3d，e）。复合材料的拉胀网格存在竞争失效机制，例如：材料断裂、节点屈服、连接单元的整体（欧拉）屈曲以及单个构件的局部（壳）屈曲，屈曲经常先于料断裂发生。破坏机制可能在变形过程中发生改变，这取决于在加载过程中在球形壳体和旋转节点中的应力状态[34–41]。

对于BCC-6H构型，晶格（ρ*> 20％）发生弹性变形以达到屈曲或塑性开始，然后由于整体屈曲和后屈服，其压缩应力在应力平稳区域内缓慢增加（图3d）。由于连接处周围壳体的广泛旋转，整个晶格承受较大的应力引起有效压缩应变。通过弹性力学的有限元分析，揭示了这种屈曲的诱导旋转现象和后屈曲效应。随着构件厚度的减小（ρ*<20％），破坏机制转变为由局部屈曲主导的破坏行为，应力-应变曲线的波动揭示了通过离散的层状裂缝引起的渐进形变（图3d）。

图2. CNT的微观结构以及力学性能

• a）重量百分比为0.5％和1％的PA12和CNT/PA12的拉伸应力-应变曲线。
• b，c）具有不同放大倍数的聚合物基质中CNT的微观结构。 
• d，e）具有不同放大倍数的CNT/PA12压缩或拉伸试样的微观结构。

• f）3D打印的轻质材料的比能量吸收率与断裂伸长率的关系图。[8,28,34–41]

Figure 2. The microstructure and mechanical properties of CNT：a) Tensile stress–strain curves of PA12 and CNT/PA12 with the weight percentages of 0.5% and 1%. Two sets of CNT s with the length range of 20–40 nm and 40–60 nm are selected. b,c) The original microstructures of CNT within the polymeric matrix with different magnifications. d,e) The microstructures of the compressed or elongated specimen of CNT/PA12 with different magnifications. f) The comparison diagram of specific energy absorption versus elongation at break for 3D-printed lightweight materials.[8,28,34–41]

图3. 不同类型晶格的结构以及力学性能

• a）BCC-6H和BCC-12H超材料的结构设计概念图（比例尺为10 mm）。
• b，c）拉胀BCC-6H和BCC-12H晶格相对于不同相对密度的结构变形。

• d，e）BCC-6H和BCC-12H晶格相对于不同相对密度的压缩应力-应变曲线。曲线下的有色 区域表示准静态压缩过程中结构吸收的能量。

Figure 3. The structure and mechanical properties of different crystal lattice：a) Illustration of the concept of structural designs of BCC-6H and BCC-12H metamaterials (Scale bar is 10 mm). b,c) The structural defor-mations of auxetic BCC-6H and BCC-12H lattices with respect to the different relative densities. d,e) The compressive stress–strain curves of BCC-6H and BCC-12H lattices with respect to the different relative densities. The colored-areas under curves indicate the energy absorbed by structures during the quasistatic compressions.

对于BCC-12H晶格受压直到其致密化，压应力与应变几乎呈线性关系。在相对密度从9.9％到33％的范围内，拉胀晶格具有显着的应变硬化效应。有限元（FE）模拟揭示了BCC-12H晶格的变形，整体屈曲起主导作用，引起了单个晶胞内的整体旋转以及分段旋转。

BCC-12H结构的竞争失效机制包括整体屈曲、节点屈服和晶胞局部断裂。每两个变形模式之间的过渡通常是由结构屈曲和材料破坏之间的激烈竞争驱动的。在严重变形的代表单元中给出了应力集中区和应变局部化区域，以预测引发每一种破坏机制所需的应力。

在图3c中，当网壳厚度接近临界状态时，一旦节点出现裂纹，破坏模式将从整体屈曲过渡到单元内的过度局部化断裂。简而言之，我们发现整体屈曲效应是引起拉胀的主要机制，主要受空心晶格行为的影响。在图3b，c中，具有较高密度的BCC-6H晶格与其他晶格相比，显示出负泊松比效应，这是十分关键的。它承受压缩后的最大体积减小，以此来消散机械能，从而显示出具有吸收能量的潜力。

古典泡沫理论允许使用Gibson和Ashby的对应关系来检测控制失效机理[8]。在图4a，b中，绘制了单位体积的能量吸收σp/Es和拉胀复合材料单位体积的峰值应力W/Es之间的关系图。这些图是针对特殊应用的能量吸收材料的选择指南，这些材料要能够控制应力和能量吸收能力。针对给定的类别，最优结果时在最大允许压应力σp下吸收最多的能量。每个晶格密度都有一个σp对应于能量曲线上的肩部。此后，最佳的能量吸收遵循各个W–σp曲线的包络线。对于泡沫结构，能量吸收能力可以表示为致密化应力σD和结构的相对密度的函数[8]。

（5）

其中A，Es，ρs，ρ和εD是经验参数，分别是杨氏模量，组成材料的密度，晶格的密度和晶格的致密化应变。在图4a，b中，包络线为线性函数，可以表示为

 （6）

η的值表示包络线的斜率，该包络线在每个σp=σD处对应不同相对密度的复合材料。在致密化应力的斜率η=1.01时，BCC-12H晶格遵循指数比例变化规律。这类似于在低应力下具有η≈1的以弯曲和屈曲为主的开孔弹性体随机泡沫[8]。然而，随着致密化应力增加几个数量级，η≈0.78的BCC-6H晶格表现出增量ΔW的下降趋势。

同时，多孔材料的 值随其相对密度 的变化以线性变化或以指数方式变化。能量吸收能力的最佳值Wmax随以下公式变化

 （7）

• a，b）拉胀型BCC-6H和BCC-12H晶格的能量吸收图，其中单位体积吸收的能量与外部压缩或冲击产生的峰值应力相对应。两个轴均通过实体模量标准化。刚好接触每条曲线的包络线就可以确定在特定应变率下泡沫的最佳选择。

• c）与弹性和塑料随机泡沫相比，高能网格（ρ*≥10％）的能量吸收能力与相对密度之间的关系。

Figure 4. Energy absorption performances of auxetic lattices. a,b) The energy absorption diagrams of auxetic BCC-6H and BCC-12H lattices with the energy absorbed per unit volume against peak stress generated by external compression or impact. Both axes are normalized by the solid modulus. The envelope which just touches each curve defines the optimum choice of foam at the specific strain-rate. c) The relationships of energy absorption capacity versus relative density for auxetic lattices (ρ*≥10％) as compared with the elastomeric and plastic stochastic foams.

图5.单位体积能量吸收与密度的阿什比图

这张图将迄今为止的拉胀超材料与其他最先进的金属和复合结构进行了比较[8,42–48]

Fig.5 Ashby map of energy absorption per unit volume versus density. This chart compares auxetic metamaterials against other most advanced metallic and composite structures so far.[8,42–48]

BCC-6H和BCC-12H的拉胀晶格比例因子γ为3和5/2。相比之下，橡胶和塑料的缩放系数只有2和3/2[8]。这些拉胀晶格具有高度各向同性，因此能量吸收能力与晶格密度不受复合材料晶格方向变化的影响，这一点通过沿不同方向压缩晶格可以证实。尤其是，BCC-6H晶格表现出显着的线性关系（γ=3），表明随着晶格密度略微增加，能量吸收会显著增加。与具有随机孔隙分布的泡沫相比，拉胀复合材料的每个晶格都经历了较大的体积变形和收缩，从而引起了普遍的体积收缩，有效地吸收了外部能量。实际上，拉胀结构的能量吸收能力也超过了蜂窝结构（平面内）、八角点阵结构和金字塔形等传统结构网格。这是因为受弯曲或拉伸应力支配的桁架格由于其局部支撑断裂或结构塌陷，导致无法同时变形以承受较大的压缩应变[13]。

这项工作解释了晶格的延展行为，并且探索了固有韧性纳米复合材料组成的非线性塑性变形结构。对于致密化应变，BCC-6H和BCC-12H的拉胀晶格分别在0.13–6.29和0.18–2.55 MJ m-3的范围内表现出优异的能量吸收能力。这些拉胀晶格的比能量吸收能力为1.6-20.42 J g-1，与镍、氧化铝纳米晶格和钛合金结构相当，并且超过了普通铝合金材料和其他增材制造结构，如图4所示[4]。相比之下，轻质钛合金泡沫具有出色的性能，然而它们的制造方法（例如EBM，SLM和激光熔覆）与目前用于轻型复合材料制造的方法（图5）相比价格昂贵且能耗高。迄今为止，微立体光刻法和直接打印技术仍受到尺寸的限制，难以实现大型结构的制作。总之，在诸如汽车和航空工业等领域中，激光烧结的拉胀复合材料的优点是能量吸收率高且经济效益好。

六、结论

1、本文通过激光烧结复合材料工艺制备了机械超材料，有效地将材料设计和结构选择结合，实现了可控制的形状转换并且验证了其出色的能量吸收能力。

2、激光烧结纳米复合超材料的能量吸收性能既取决于结构韧带的变形机理，又取决于组成材料的微观结构。

3、增材制造技术将高性能纳米复合材料与晶体的设计相结合，可以用于缓冲、防冲击、防爆以及包装等领域。

4、复合材料超材料由于其膨胀变形、可控制行为、高可预测性以及出色的能量吸收能力，在航空航天和军事等行业中具有巨大潜力。

• [1] MOHSENIZADEH M, GASBARRI F, MUNTHER M, et al.Additively-manufactured lightweight Metamaterials for energy absorption[J].Materials and Design, 2018, 139: 521–530.

• [2] CHEN Y, LI T, JIA Z, et al. 3D printed hierarchicalhoneycombs with shape integrity under large compressive deformations[J].Materials and Design, 2018, 137: 226–234.

• [3] MUELLER J, RANEY J R, SHEA K, et al. ArchitectedLattices with High Stiffness and Toughness via Multicore–Shell 3D Printing[J].Advanced Materials, 2018, 30(12): 1–8.

• [4] FAN S feng, ZHANG T, YU K, et al. Compressiveproperties and energy absorption characteristics of open-cell nickel foams[J].Transactions of Nonferrous Metals Society of China (English Edition), 2017,27(1): 117–124.

• [5] MASKERY I, ABOULKHAIR N T, AREMU A O, et al.Compressive failure modes and energy absorption in additively manufactureddouble gyroid lattices[J]. Additive Manufacturing, 2017, 16: 24–29.

• [6] KRUK S, KIVSHAR Y. Functional Meta-Optics and NanophotonicsGovern by Mie Resonances[J]. ACS Photonics, 2017, 4(11): 2638–2649.

• [7] LU M H, FENG L, CHEN Y F, et al. Phononic crystalsand acoustic metamaterials[J].Materials Today, 2009, 12(12): 34–42.

• [8] YUAN S, CHUA C K, ZHOU K, et al. 3D-Printed MechanicalMetamaterials with High Energy Absorption[J]. Advanced Materials Technologies,2019, 4(3): 1–9.

• [9] FRENZEL T, KADIC M, WEGENER M, et al.Three-dimensional mechanical metamaterials with a twist[J]. Science, 2017,358(6366): 1072–1074.

• [10] WANG Q, JACKSON J A, GE Q, et al. LightweightMechanical Metamaterials with Tunable Negative Thermal Expansion[J]. PhysicalReview Letters, 2016, 117(17): 1–6.

• [11] AKSELROD G M, ARGYROPOULOS C, HOANG T B, et al.Probing the mechanisms of large Purcell enhancement in plasmonicnanoantennas[J]. Nature Photonics, 2014, 8(11): 835–840.

• [12] BABAEE S, SHIM J, WEAVER J C, et al. 3D softmetamaterials with negative poisson’sratio[J]. Advanced Materials, 2013,25(36): 5044–5049.

• [13] ZHENG X, LEE H, WEISGRABER T H, et al. Ultralight,ultrastiff mechanical metamaterials[J]. Science, 2014, 344(6190): 1373–1377.

• [14] BAUGHMAN R H.Avoiding the shrink[J]. Nature, 2003, 425(6959): 667.

• [15] GREAVES G N, GREER A L, LAKES R S, et al. Poisson’sratio and modern materials[J]. Nature Materials, 2011, 10(11): 823–837.

• [16] MEZA L R, DAS S, GREER J R.,et al. Strong,lightweight, and recoverable three-dimensional ceramic nanolattices[J].Science, 2014, 345(6202): 1322–1326.

• [17] BAUER J, SCHROER A, SCHWAIGER R, et al. Approachingtheoretical strength in glassy carbon nanolattices[J]. Nature Materials, 2016,15(4): 438–443.

• [18] YUAN S, SHEN F, BAI J, et al. 3D soft auxeticlattice structures fabricated by selective laser sintering: TPU powderevaluation and process optimization[J]. Materials and Design, 2017, 120:317–327.

• [19] DUOSS E B, WEISGRABER T H, HEARON K, et al.Three-dimensional printing of elastomeric, cellular architectures with negativestiffness[J]. Advanced Functional Materials, 2014, 24(31): 4905–4913.

• [20] SHUAI C, GAO C, NIE Y, et al. Structure andproperties of nano-hydroxypatite scaffolds for bone tissue engineering with aselective laser sintering system[J]. Nanotechnology, 2011, 22(28).

• [21] SHEN F, YUAN S, CHUA C K, et al. Development ofprocess efficiency maps for selective laser sintering of polymeric compositepowders: Modeling and experimental testing[J]. Journal of Materials ProcessingTechnology, 2018, 254(July 2017): 52–59.

• [22] TIAN X, PENG G, YAN M, et al. Process prediction ofselective laser sintering based on heat transfer analysis for polyamidecomposite powders[J]. International Journal of Heat and Mass Transfer, 2018,120: 379–386.

• [23] YUAN S, BAI J, CHUA C K, et al. Material evaluationand process optimization of CNT-coated polymer powders for selective lasersintering[J]. Polymers, 2016, 8(10).

• [24] YUAN S, BAI J, CHUA C K, et al. Highly enhancedthermal conductivity of thermoplastic nanocomposites with a low mass fractionof MWCNTs by a facilitated latex approach[J]. Composites Part A: AppliedScience and Manufacturing, 2016, 90: 699–710.

• [25] YAN W, GE W, QIAN Y, et al. Multi-physics modelingof single/multiple-track defect mechanisms in electron beam selectivemelting[J]. Acta Materialia, 2017, 134: 324–333.

• [26] WAGNER H D, VAIA R A. Nanocomposites: Issues at theinterface[J]. Materials Today, 2004, 7(11): 38–42..

• [27] PAUL D R, ROBESON L M. Polymer nanotechnology:Nanocomposites[J]. Polymer, 2008, 49(15): 3187–3204.

• [28] COMPTON B G, LEWIS J A. 3D-printing of lightweightcellular composites[J]. Advanced Materials, 2014, 26(34): 5930–5935.

• [29] NING F, CONG W, HU Y, et al. Additive manufacturingof carbon fiber-reinforced plastic composites using fused deposition modeling:Effects of process parameters on tensile properties[J]. Journal of CompositeMaterials, 2017, 51(4): 451–462.

• [30] MOHAN N, SENTHIL P, VINODH S, et al. A review oncomposite materials and process parameters optimisation for the fuseddeposition modelling process[J]. Virtual and Physical Prototyping, 2017, 12(1):47–59.

• [31] LEWICKI J P, RODRIGUEZ J N, ZHU C, et al.3D-Printing of Meso-structurally Ordered Carbon Fiber/Polymer Composites withUnprecedented Orthotropic Physical Properties[J]. Scientific Reports, 2017,7(December 2016): 1–14.

• [32] RANEY J R, COMPTON B G, MUELLER J, et al. Rotational3D printing of damage-tolerant composites with programmable mechanics[J].Proceedings of the National Academy of Sciences of the United States ofAmerica, 2018, 115(6): 1198–1203.

• [33] ZHANG X, YAO J, LIU B, et al. Three-DimensionalHigh-Entropy Alloy-Polymer Composite Nanolattices That Overcome theStrength-Recoverability Trade-off[J]. Nano Letters, 2018, 18(7): 4247–4256.

• [34] STANEV L, DRENCHEV B, YOTOV A. CompressiveProperties and Energy Absorption Behaviour of AlSi10Mg Open-Cell Foam[J]. 2014,22(1): 44–53. .

• [35] PARTHASARATHY J, STARLY B, RAMAN S, et al.Mechanical evaluation of porous titanium (Ti6Al4V) structures with electronbeam melting (EBM)[J]. Journal of the Mechanical Behavior of BiomedicalMaterials, 2010, 3(3): 249–259.

• [36] GHITA O, JAMES E, DAVIES R, et al. High TemperatureLaser Sintering (HT-LS): An investigation into mechanical properties andshrinkage characteristics of Poly (Ether Ketone) (PEK) structures[J]. Materialsand Design, 2014, 61: 124–132.

• [37] GAIKWAD S, TATE J S, THEODOROPOULOU N, et al.Electrical and mechanical properties of PA11 blended with nanographeneplatelets using industrial twin-screw extruder for selective lasersintering[J]. Journal of Composite Materials, 2013, 47(23): 2973–2986.

• [38] BAI J, YUAN S, SHEN F, et al. Toughening ofpolyamide 11 with carbon nanotubes for additive manufacturing[J]. Virtual andPhysical Prototyping, 2017, 12(3): 235–240.

• [39] MARTIN J H, YAHATA B D, HUNDLEY J M, et al. 3Dprinting of high-strength aluminium alloys[J]. Nature, 2017, 549(7672):365–369.

• [40] RAFI H K, KARTHIK N V., GONG H, et al.Microstructures and mechanical properties of Ti6Al4V parts fabricated byselective laser melting and electron beam melting[J]. Journal of MaterialsEngineering and Performance, 2013, 22(12): 3872–3883.

• [41] MURR L E, GAYTAN S M, RAMIREZ D A, et al. MetalFabrication by Additive Manufacturing Using Laser and Electron Beam MeltingTechnologies[J]. Journal of Materials Science and Technology, 2012, 28(1):1–14.

• [42] MASKERY I, ABOULKHAIR N T, AREMU A O, et al. Amechanical property evaluation of graded density Al-Si10-Mg lattice structuresmanufactured by selective laser melting[J]. Materials Science and EngineeringA, 2016, 670: 264–274.

• [43] JUNG A, NATTER H, DIEBELS S, et al. Nanonickelcoated aluminum foam for enhanced impact energy absorption[J]. AdvancedEngineering Materials, 2011, 13(1–2): 23–28.

• [44] ALDOSHAN A, KHANNA S. Effect of relative density onthe dynamic compressive behavior of carbon nanotube reinforced aluminumfoam[J]. Materials Science and Engineering A, 2017, 689(January): 17–24.

• [45] CHOY S Y, SUN C N, LEONG K F, et al. Compressiveproperties of functionally graded lattice structures manufactured by selectivelaser melting[J]. Materials and Design, 2017, 131(May): 112–120.

• [46] XIE B, FAN Y Z, MU T Z, et al. Fabrication andenergy absorption properties of titanium foam with CaCl2 as a space holder[J].Materials Science and Engineering A, 2017, 708(July): 419–423.

• [47] ALVANDI-TABRIZI Y, RABIEI A. Use of Composite MetalFoam for Improving Absorption of Collision Forces[J]. Procedia MaterialsScience, 2014, 4(2004): 377–382.

• [48] QIAO J C, XI Z P, TANG H P, et al. Compressiveproperty and energy absorption of porous sintered fiber metals[J]. MaterialsTransactions, 2008, 49(12): 2919–2921.