3D打印梯度Gyroid结构的动态冲击响应

厉雪; 肖李军; 宋卫东

doi:10.11858/gywlxb.20210701

3D打印梯度Gyroid结构的动态冲击响应

doi: 10.11858/gywlxb.20210701

北京理工大学爆炸科学与技术国家重点实验室，北京　100081

基金项目: 国家自然科学基金（11972092，12002049）

详细信息

作者简介:
厉　雪（1996－），女，硕士研究生，主要从事材料与结构的冲击动力学研究. E-mail：15152103981@163.com

通讯作者:
宋卫东（1975－），男，博士，教授，主要从事材料与结构的冲击动力学研究. E-mail：swdgh@bit.edu.cn

中图分类号: O347.1
计量
- 文章访问数: 4317
- HTML全文浏览量: 1921
- PDF下载量: 77
出版历程
- 收稿日期: 2021-01-04
- 修回日期: 2021-01-19

Dynamic Behavior of 3D Printed Graded Gyroid Structures under Impact Loading

State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 利用ANSYS/LS-DYNA对均匀及梯度Gyroid结构进行准静态与动态压缩数值模拟，分析其应力分布、变形模式、承载能力以及吸能特性。对3D打印的316L不锈钢试样实施了单轴拉伸实验，获取了相应的材料参数，建立了Gyroid结构有限元模型，进而对其动态力学响应进行了数值仿真。结果表明：均匀结构呈现出较均匀的变形模式，梯度结构为低密度端向高密度端传播的逐层变形模式；两种结构均呈现明显的应变率敏感性，且负梯度结构的应变率敏感性最明显；在相同的加载速度下，负梯度结构的吸能效率最高，且具有最低的支撑端应力，是最佳的防护结构。研究结果可为冲击载荷下防护结构的设计选型提供参考。
- 最小三维周期曲面 /
- Gyroid结构 /
- 3D打印 /
- 抗冲击性能 /
- 变形模式 /
- 吸能特性
Abstract: The quasi-static and dynamic compression properties of uniform and graded Gyroid structures were studied by using the commercial software ANSYS/LS-DYNA. The stress distributions, deformation mode, loading bearing and energy absorption abilities were analyzed accordingly. Some numerical material parameters of SLM (selective laser melting) printed 316L stainless steel was obtained by tensile tests. Finite element models of Gyroid structure were established, and the dynamic mechanical response was numerically simulated. The results indicate that the uniform structure performs relatively uniform deformation patterns, while the graded specimen shows deformation propagation from the low density end to the high density end. All the studied structures show obvious strain rate sensitivities, which is the most apparent in the negative structures. Moreover, the negative gradient structure absorbs the most energy and possesses the smallest supporting stresses compared with the other structures under the similar loading velocity, which denotes that the negative arrangement is the perfect protect structures. The results can provide guidance for the design of the protect structures under impact loading.
- triplyperiodicminimalsurface /
- Gyroid structures /
- 3D printing /
- impact resistance /
- deformation mode /
- energy absorption abilities

HTML全文

高性能轻质复合材料越来越多地应用于工程材料领域，尤其是贝壳珍珠层的应用最广泛。贝壳珍珠层的最大优点是高强度和高抗破裂韧性，这些优越性能得益于其独特的多层次微结构。珍珠层由体积分数为95%的文石碳酸钙和5%的有机质组成，两者交叉堆叠成有序的层状结构^[1]，即浆砌层合结构，其中有机质的桥连作用以及诱导裂纹偏转作用强化了材料的韧性^[2-3]。根据贝壳这种浆砌层合结构，可制备具有硬材料基体和黏结层的复合材料。其多层次微结构如图1^[4]所示。

图 1 珍珠层微结构^[4]

Figure 1. Microstructure of nacre^[4]

下载: 全尺寸图片幻灯片

鉴于贝壳珍珠层这种浆砌层合结构的诸多优点，随之出现了很多浆砌结构材料制造技术，其中3D打印技术作为新兴的制造技术，被越来越多地用在工业生产和人们的生活中，如航空航天中的机翼、生物材料中的牙齿等，同时也使仿生浆砌层合结构复合材料的制备更方便。马骁勇等^[5]研究发现，三维打印贝壳仿生结构断裂发生在类橡胶软材料内部，因此浆砌层合结构内部层间的分层与脱黏对整体结构力学性能有着重要影响。而浆砌层合结构的软硬交替界面的力学性能可通过分析复合材料的层间断裂韧性来研究，为此本工作着眼于3D打印复合材料层间断裂韧性。

近些年来，针对复合材料层间断裂韧性的有限元计算模型发展迅速。Xu等^[6]的内聚力模型发展较完善，且具有自身的理论原理，广泛应用于复合材料断裂与分层破坏研究中。其中，许多学者将内聚力模型引入标准复合材料层间扩展试验模型中，如Hosseini等^[7]基于内聚力模型研究了编织玻璃环氧基复合材料的断裂韧性，Hua等^[8]借助内聚力模型研究了玻璃钢铝层合板的非对称双悬臂梁，宗要武^[9]基于内聚力模型开展了钢纤维水泥基材料的界面性能分析。众多学者还将内聚力模型引入冲击载荷下复合材料分层损伤的有限元模拟中，较好地反映了复合材料结构分层损伤破坏的全过程。如Alfaro等^[10]基于内聚力模型研究了三维层合板的断裂和分层，Liu等^[11]基于内聚力模型研究了层合板动态双悬臂梁，赵丽滨等^[12]总结了纤维增强复合材料层合板分层行为研究进展以及内聚力方法。在上述模型基础上，本研究在内聚力模型中考虑了3D打印浆砌层合结构的裂纹萌生和扩展行为，并分析了模型参数的影响。

基于此，针对3D打印浆砌层合结构复合材料的I型和II型断裂韧性进行有限元模拟，通过有限元数值模拟与模型试验的对比分析，验证有限元模拟的正确性，随后考察初始裂纹长度、断裂韧性、起始界面刚度、界面强度、黏结层厚度以及净距对3D打印浆砌层合结构复合材料层间力学性能的影响。

1. 层间断裂韧性的有限元模拟

1.1 内聚力本构模型及损伤判据

采用双线性拉伸分离准则来模拟黏结层接触面的黏性力学行为，其本构关系如图2所示。

图 2 双线性拉伸分离准则

Figure 2. Bilinear stretching separation criteria

下载: 全尺寸图片幻灯片

用内聚力单元来模拟材料的损伤破坏过程，分为两个阶段：损伤起始和损伤演化^[13]。

针对损伤起始，采用二次应力准则

${\left( {\frac{{\langle {\sigma _n}\rangle }}{{{N_{\max }}}}} \right)^2} + {\left( {\frac{{{\sigma _s}}}{{{S_{\max }}}}} \right)^2} + {\left( {\frac{{{\sigma _t}}}{{{T_{\max }}}}} \right)^2} = 1$

(1)

式中：n是单元法向，对应于I型断裂；s、t是单元的两个切线方向，对应于II型和III型断裂； ${\sigma _n}$ 、 ${\sigma _s}$ 、 ${\sigma _t}$ 分别为法向和两个切向的应力， ${N_{\max }}$ 、 ${S_{\max }}$ 、 ${T_{\max }}$ 分别为法向和两个切向的峰值强度。

针对损伤演化，采用基于能量的损伤演化BK准则

${G\rm_C} = {G\rm_{IC}} + ({G\rm_{IIC}} - {G\rm_{IC}}){\left( {\frac{{{G\rm_{shear}}}}{{{G\rm_T}}}} \right)^{\eta} }$

(2)

${G\rm_{shear}} = {G\rm_{II}} + {G\rm_{III}}$

(3)

${G\rm_T} = {G\rm_I} + {G\rm_{shear}}$

(4)

式中： ${G\rm_{IC}}$ 、 ${G\rm_{IIC}}$ 分别为法向和II型切向的断裂韧性值， ${G\rm_I}$ 、 ${G\rm_{II}}$ 、 ${G\rm_{III}}$ 分别为法向和两个切向的能量释放率，η为材料常数。

损伤起始后会有应力的变化，需要引入刚度弱化系数D

$D = \frac{{\delta _m^{\rm f}\left( {\delta _m^{\max } - \delta _m^0} \right)}}{{\delta _m^{\max }\left( {\delta _m^{\rm f} - \delta _m^0} \right)}}$

(5)

$\delta _m^0 = \sqrt {{{\langle {\delta _n}\rangle }^2} + {\delta _s^2} + {\delta _t^2}}$

(6)

$\delta _m^{\rm f} = \frac{{2{G\rm_C}}}{{\sqrt {{{\langle {\sigma _n}\rangle }^2} + {\sigma _s^2} + {\sigma _t^2}} }}$

(7)

式中： $\delta _m^0$ 为破坏起始时单元节点的张开量， $\delta _m^{\rm f}$ 为破坏演化时节点的张开量， $\delta _m^{\max }$ 为在荷载历程中节点的最大张开量。

1.2 I型和II型断裂韧性有限元数值模型

采用ABAQUS进行I型材料断裂韧性的有限元数值模拟，如图3所示。数值模拟选用三维建模，长240 mm，宽37 mm，初始裂纹长32 mm。采用动态模拟准静态，进行复合材料裂纹张开过程模拟，边界条件：保持试样右端固定，在距试样左端14 mm处绑定一个铁片，施加竖向位移荷载，在铁片控制点处实现材料的张开。黏结层采用内聚力单元和扫略网格描述，且布置局部边界种子，尺寸为0.5 mm。

图 3 I型有限元模型及示意图

Figure 3. Finite element model and schematic of model-I

下载: 全尺寸图片幻灯片

采用ABAQUS进行II型材料断裂韧性的有限元数值模拟，如图4所示。数值模拟采用三维模型，长140 mm，宽25 mm，净距80 mm，初始裂纹长20 mm。采用动态模拟准静态，进行复合材料裂纹错开过程模拟，试样上表面中心线与参考点绑定，其边界条件：保持两个支座固定，通过对参考点施加竖向位移荷载实现材料的错开。黏结层采用内聚力单元和扫略网格描述，且布置局部边界种子，尺寸为0.5 mm。考虑到支座处和初始裂纹处上下表面的接触，采用了通用接触。

图 4 II型有限元模型及示意图

Figure 4. Finite element model and schematic of model-II

下载: 全尺寸图片幻灯片

通过对实验数据结果的分析，确定了内聚力模型层间断裂韧性和节点张开量。由图2的本构关系可知，层间断裂韧性为图形与坐标轴围成的面积，进而确定界面强度，初始斜率为起始界面刚度。具体参数如表1所示，其中T为黏结层厚度，K_nn和K_ss为I型和II型模型的刚度参数。

表 1 内聚力参数

Table 1. Cohesive parameters

T/mm	N_max/ MPa	S_max/ MPa	G_I/(N∙mm⁻¹)	G_II/(N∙mm⁻¹)	K_nn/(N∙mm⁻³)	K_ss/(N∙mm⁻³)
0.16	0.35	2.0	0.40	0.73	15 625	31 250
0.20	0.50	2.5	0.55	0.65	12 500	25 000

下载: 导出CSV

| 显示表格

1.3 数值模拟与实验结果对比验证

为了验证有限元计算结果的正确性，进行了I型和II型断裂韧性的模型试验。试样由Objet 260 Connex 3光固化三维打印机制备。通过SolidWorks软件设计模型并导入三维打印机程序中，然后利用双喷头3D打印机将模型打印成最终试样。试样尺寸与有限元模型尺寸一致，试样上下两层硬材料采用不透明灰色树脂材料VeroGray，通过准静态实验，得到其弹性模量和泊松比分别为1 250 MPa和0.3；中间黏结层采用类橡胶高分子材料TangoBlack，通过准静态实验，得到其弹性模量和泊松比分别为0.4 MPa和0.33。采用ASTM D5528-01^[14]标准进行I型断裂韧性试验，如图5(a)所示，通过试样拉伸实现裂纹张开，得到荷载-位移曲线；采用NASA/TM-2010-216838^[15]标准进行II型断裂韧性试验，如图5(b)所示，通过试样弯曲实现层间裂纹错开，得到荷载-位移曲线。

图 5 断裂韧性测试试验照片

Figure 5. Picture of fracture toughness test

下载: 全尺寸图片幻灯片

DCB和ENF试验与模拟结果的对比分别如图6和图7所示，其中A点和B点分别代表黏结层局部损伤起始点和损伤演化点^[16-17]。由图6和图7可知，不同黏结层厚度的数值模拟结果与试验测量所得荷载-位移曲线具有相同的变化趋势。损伤起始点和损伤演化点的数值模拟结果与试验实测值吻合较好，从而验证了数值建模过程的可行性及模拟结果的可靠性。因此，本工作采用该建模方法进行后续数值模拟研究。

图 6 不同黏结层厚度的I型对比

Figure 6. Comparison chart of different bonding layer thicknesses in model-I

下载: 全尺寸图片幻灯片

图 7 不同黏结层厚度的II型对比

Figure 7. Comparison chart of different bonding layer thicknesses in model-II

下载: 全尺寸图片幻灯片

对于I型试样，以黏结层厚度为0.20 mm的工况为例进行分析，如图6所示。由图6可知，荷载-位移曲线的数值模拟结果与实测值具有相同的变化趋势，且能较好地预测局部损伤起始点和损伤演化点，但因存在损伤致使试样实际变形大于数值模拟结果，所以数值模拟所得拉伸力比实测值略小。

对于II型试样，以黏结层厚度为0.16 mm的工况为例进行分析，如图7所示。可见，荷载-位移曲线的数值模拟结果与试验模型具有相同的变化趋势，且能较好地预测损伤起始点。损伤演化点的荷载变小，这是由于试验过程中试样产生微弱错动，一定程度上延缓了裂纹的损伤，增大了试样的弯曲力。

2. 数值模拟结果分析

2.1 I型数值模拟结果分析

针对中间黏结层厚度为0.20 mm的I型模拟结果进行初始裂纹长度、断裂韧性、起始界面刚度与界面强度的研究分析，并讨论黏结层厚度的影响。

2.1.1 初始裂纹长度的影响

不同初始裂纹长度C的荷载-位移曲线如图8所示。由图8可知，I型试样荷载峰值所对应的位移随初始裂纹长度的增大而增大，拉伸力峰值略有增大。由修正梁理论^[14]公式推断可知：当断裂韧性不变时，失效裂纹长度增大，在拉伸力峰值增幅不大的情况下，荷载峰值所对应的位移会增大。

图 8 初始裂纹长度对I型荷载-位移曲线的影响

Figure 8. Influence of initial crack length on model-I load-displacement curve

下载: 全尺寸图片幻灯片

此外，在局部损伤起始点之前，荷载增长速率随初始裂纹长度的增大而减小。如图9所示，在损伤演化点之前，达到损伤演化点的裂纹尖端开口位移(CTOD)相同，裂纹长度增大182%，即由22 mm增大到62 mm；拉伸位移增大160%，即由20 mm增大到52 mm，表明试样的裂纹拉伸张开呈非线性，这是由试样弯曲变形和尖端黏结层阻碍共同导致。在损伤演化后，载荷下降速度随初始裂纹长度的增大而加快，表明初始裂纹长度越小，材料的韧性越强，材料不会瞬时失去承载力。

图 9 初始裂纹长度对裂纹尖端开口位移曲线的影响

Figure 9. Influence of initial crack length on crack tip opening displacement curves

下载: 全尺寸图片幻灯片

2.1.2 断裂韧性的影响

不同断裂韧性的荷载-位移曲线如图10所示。由图10可知，拉伸力峰值随断裂韧性的增强而增大。这是由于断裂韧性增强，层间抵抗破坏能力增大，因此需要更大的拉伸力峰值使材料破坏。保持相同的断裂韧性增加值0.1 N/mm，则拉伸力峰值的增加幅度为4 N，说明断裂韧性对I型试样层间破坏具有显著影响。

图 10 断裂韧性对I型荷载-位移曲线的影响

Figure 10. Influence of fracture toughness on model-I load-displacement curve

下载: 全尺寸图片幻灯片

2.1.3 起始界面刚度的影响

不同起始界面刚度的荷载-位移曲线如图11所示，K_nn = 12 500 N/mm³为模拟试样黏结层厚度为0.20 mm时的刚度参数值。由图11可知，拉伸力峰值随起始界面刚度的增大而略有下降。这是由于起始界面刚度的改变仅影响破坏起始时的单元节点张开量，不会影响整体材料的破坏能。因此刚度参数值在一定范围内变化时，对荷载-位移曲线的影响并不显著。

图 11 起始界面刚度对I型荷载-位移曲线的影响

Figure 11. Influence of interface stiffness on model-I load-displacement curve

下载: 全尺寸图片幻灯片

2.1.4 界面强度的影响

不同界面强度的荷载-位移曲线如图12所示。由图12可知，拉伸力峰值随界面强度的增加无明显变化，而荷载峰值所对应的位移减小。这是由于界面强度越大，损伤演化时的节点错开量越小，导致试样破坏时张开位移变小，荷载峰值所对应的位移减小。

图 12 界面强度对I型荷载-位移曲线的影响

Figure 12. Influence of interface strength on model-I load-displacement curve

下载: 全尺寸图片幻灯片

2.1.5 黏结层厚度的影响

不同黏结层厚度的荷载-位移曲线如图6所示。由图6可知，拉伸力峰值随着黏结层厚度的增加而增加，且达到拉伸力峰值对应的位移增大。保持相同的黏结层厚度增加值0.04 mm，则拉伸力峰值增量减小，由4.72 N降为1.51 N。当黏结层厚度增加到一定程度后，拉伸力峰值将不会产生明显的增长。

2.2 II型数值模拟结果分析

针对中间黏结层厚度为0.16 mm的II型模拟结果进行初始裂纹长度、断裂韧性、起始界面刚度、界面强度与净距的研究分析，并讨论黏结层厚度的影响。

2.2.1 初始裂纹长度的影响

不同初始裂纹长度C的荷载-位移曲线如图13所示。由图13可知，II型试样的弯曲力峰值随初始裂纹长度的增大而减小。初始裂纹长度越小，起始刚度越大，达到弯曲力峰值后，黏结层产生破坏并引起弯曲力下降，且下降幅度较大。这是由于黏结层剪切力增大，层间损伤抵抗力增强，导致弯曲力峰值增大。当试样处于损伤演化后期，瞬时黏性层破坏长度较大，导致荷载骤然下降，曲线出现剧烈抖动。同样表明初始裂纹长度越小，材料的韧性越强，越不容易破坏，但材料破坏后承载力降幅很大。

图 13 初始裂纹长度对II型荷载-位移曲线的影响

Figure 13. Influence of initial crack length on model-II load-displacement curve

下载: 全尺寸图片幻灯片

2.2.2 断裂韧性的影响

不同断裂韧性的荷载-位移曲线如图14所示。由图14可知，弯曲力峰值随断裂韧性的增加而增大。这是由于黏结层抵抗破坏的能力随断裂韧性的增强而增强，试样越不容易发生破坏。

图 14 断裂韧性对II型荷载-位移曲线的影响

Figure 14. Influence of fracture toughness on model-II load-displacement curve

下载: 全尺寸图片幻灯片

此外，弯曲力峰值与断裂韧性之间的关系曲线如图15所示。由图15可知，保持相同的断裂韧性增加值0.05 N/mm，弯曲力峰值增长率变化依次为120、40、30，表明断裂韧性的增强对弯曲力峰值的影响越来越小。

图 15 断裂韧性对应的荷载峰值变化曲线

Figure 15. Loading peak curve corresponding to fracture toughness

下载: 全尺寸图片幻灯片

2.2.3 起始界面刚度的影响

不同起始界面刚度的荷载-位移曲线如图16所示。由图16可知，荷载-位移曲线随起始界面刚度的增大无显著变化。这是由于起始界面刚度仅影响试样破坏起始时的单元节点错开量，不会影响整体试样的破坏能。

图 16 起始界面刚度对II型荷载-位移曲线的影响

Figure 16. Influence of interface stiffness on model-II load-displacement curve

下载: 全尺寸图片幻灯片

2.2.4 界面强度的影响

不同界面强度的荷载-位移曲线如图17所示。由图17可知，弯曲力峰值随界面强度的增加而增大。这是由于保持断裂韧性不变，破坏演化阶段的节点错开量随起始界面刚度的增大而减小，达到单元破坏所需的应力增大，因此所需的外荷载增大。

图 17 界面强度对II型荷载-位移曲线的影响

Figure 17. Influence of interface strength on model-II load-displacement curve

下载: 全尺寸图片幻灯片

图17中S_max = 2.0 MPa为模拟II型试样黏结层的应力参数值。可见，当S_max大于2.0 MPa时，节点错开量减小，中间黏结层在达到损伤演化点前发生破坏，导致曲线在弯曲力峰值点前出现抖动；当S_max小于2.0 MPa时，节点位移错开量增加，中间黏结层能够抵挡损伤破坏，曲线平滑。

2.2.5 净距的影响

不同净距L下的荷载-位移曲线如图18所示。由图18可知，弯曲力峰值随净距的增大而减小，荷载峰值对应的位移增大。

图 18 净距对II型荷载-位移曲线的影响

Figure 18. Influence of clear distance on model-II load-displacement curve

下载: 全尺寸图片幻灯片

净距的改变影响了初始裂纹长度和试样跨距两部分。弯曲力峰值随初始裂纹长度的增大而减小(见图13)。随着跨距增大，试样整体抗弯刚度降低，位移加载相同距离时所对应的弯曲力减小，因此试样破坏时所对应的位移增大。因此初始裂纹长度和试样跨距共同导致了弯曲力峰值的减小以及荷载峰值对应位移的增大。

2.2.6 黏结层厚度的影响

不同黏结层厚度条件下的荷载-位移曲线如图7所示。由图7可知，弯曲力峰值随黏结层厚度的增大而减小，弯曲力峰值对应的位移在15.5 mm左右。保持相同的黏结层厚度增加值0.04 mm，则弯曲力峰值减小量降低，由17.93 N降为2.74 N。因此当黏结层厚度增大到一定程度后，弯曲力峰值同样不会产生明显的增长。

3. 结　论

通过有限元模拟对3D打印浆砌层合结构复合材料的I型和II型断裂韧性进行了分析，研究了初始裂纹长度、断裂韧性、起始界面刚度、界面强度、黏结层厚度以及净距等参数对3D打印浆砌层合结构复合材料层间力学性能的影响，得到以下结论。

(1) 对于3D打印浆砌层合结构复合材料I型模型，减小初始裂纹长度、增加断裂韧性值以及增大黏结层厚度都能够增强材料韧性，使层间具有较高的承载能力。其中损伤演化后荷载下降速度随初始裂纹长度的减小而平稳缓慢下降，起始界面刚度和界面强度的改变对荷载峰值无显著变化，界面强度的增大会引起载荷峰值对应位移的减小。

(2) 对于3D打印浆砌层合结构复合材料II型模型，减小初始裂纹长度、增加断裂韧性、增强界面强度、减小黏结层厚度和减小净距都能够增强材料韧性，使层间具有较高的承载能力。其中II型断裂韧性的增大对荷载峰值的影响越来越小；净距的减小会引起载荷峰值对应位移的减小，引起材料瞬间破坏；起始界面刚度对试样层间性能无显著影响。

图 1 t与Gyroid结构相对密度的关系

Figure 1. Relationship of t and relative density of the Gyroid structures

下载: 全尺寸图片幻灯片

图 2 均匀及梯度Gyroid结构

Figure 2. Uniform and graded Gyroid structures

下载: 全尺寸图片幻灯片

图 3 SLM打印的316L不锈钢试样的工程应力-应变曲线

Figure 3. Engineering stress-strain curve of SLM printed 316L stainless steel specimens

下载: 全尺寸图片幻灯片

图 4 SLM打印的均匀TPMS结构冲击实验与数值模拟对比

Figure 4. Comparison of the experiment and numerical simulation of the SLM printed TPMS structures under impact loading

下载: 全尺寸图片幻灯片

图 5 均匀和梯度Gyroid结构在不同加载速度下的冲击端应力-应变曲线

Figure 5. Impact end stress-strain curves of the uniform and graded Gyroid structures under different loading velocities

下载: 全尺寸图片幻灯片

图 6 均匀结构（a）、正梯度结构（b）、负梯度结构（c）的准静态能量吸收效率曲线以及不同加载速度下各结构的密实应变（d）

Figure 6. Efficient energy curves of uniform structures (a), positive gradient structures (b) and negative gradient structures (c) under quasi-static loading, and the densification strain of different structures under different loading velocities (d)

下载: 全尺寸图片幻灯片

图 7 不同结构冲击端和支撑端的屈服应力对比

Figure 7. Stress comparison between the impact end and support end of different structures

下载: 全尺寸图片幻灯片

图 8 均匀和梯度Gyroid结构在30 m/s加载速度下的变形模式：（a）均匀结构，（b）正梯度结构，（c）负梯度结构

Figure 8. Deformation evolution of the uniform and graded Gyroid structures under v = 30 m/s: (a) uniform structures, (b) positive gradient structures, (c) negative gradient structures

下载: 全尺寸图片幻灯片

图 9 冲击波模式下均匀结构（a）和负梯度结构（b）的应力-应变曲线

Figure 9. Stress-strain curves of the uniform (a) and negative gradient (b) structures under shock wave

下载: 全尺寸图片幻灯片

图 10 冲击波模式下均匀结构（a）和负梯度结构（b）的变形模式

Figure 10. Deformation evolution of the uniform (a) and negative gradient (b) structures under shock wave modes

下载: 全尺寸图片幻灯片

图 11 不同加载速度下各结构的比吸能对比

Figure 11. Comparison of specific energy absorption between different structures under different loading velocities

下载: 全尺寸图片幻灯片

表 1 不同加载条件下各结构冲击端屈服应力 ${\sigma }_{\mathrm{y}}$ 、平台应力 ${\sigma }_{\mathrm{p}}$ 和应力放大因子 ${\sigma _{{\rm{DIF}}}}$

Table 1. Impact end yield stress ${\sigma }_{\mathrm{y}}$ , plateau stress ${\sigma }_{\mathrm{p}}$ , and stress increased factors ${\sigma _{{\rm{DIF}}}}$ of different structures under different loading velocities

Structures	v/(m·s⁻¹)	$\sigma{\rm{_y}}$ /MPa	$\sigma{\rm{_p}}$ /MPa	$\sigma{\rm{_{DIF}}}$
Uniform structures	1	26.80	39.86	1.00
	10	50.20	61.56	1.87
	30	74.00	66.91	2.76
	50	96.64	70.02	3.61
Positive gradient structures	1	14.81	44.87	1.00
	10	21.66	61.54	1.46
	30	30.60	77.04	2.07
	50	68.68	78.70	4.64
Negative gradient structures	1	15.69	40.86	1.00
	10	48.53	75.24	3.09
	30	70.61	76.76	4.50
	50	139.68	81.14	8.90

下载: 导出CSV

参考文献(22)

[1]	MONTAZERIAN H, DAVOODI E, ASADI-EYDIVAND M, et al. Porous scaffold internal architecture design based on minimal surfaces: a compromise between permeability and elastic properties [J]. Materials & Design, 2017, 126: 98–114. doi: 10.1016/j.matdes.2017.04.009
[2]	SHEN M H, QIN W, XING B H, et al. Mechanical properties of 3D printed ceramic cellular materials with triply periodic minimal surface architectures [J]. Journal of the European Ceramic Society, 2021, 41(2): 1481–1489. doi: 10.1016/j.jeurceramsoc.2020.09.062
[3]	Al-KETAN O, ROWSHAN R, AL-RUB R K A. Topology-mechanical property relationship of 3D printed strut, skeletal, and sheet based periodic metallic cellular materials [J]. Additive Manufacturing, 2018, 19: 167–183. doi: 10.1016/j.addma.2017.12.006
[4]	MASKERY I, ABOULKHAIR N T, AREMU A O, et al. Compressive failure modes and energy absorption in additively manufactured double gyroid lattices [J]. Additive Manufacturing, 2017, 16: 24–29. doi: 10.1016/j.addma.2017.04.003
[5]	BONATTI C, MOHR D. Mechanical performance of additively-manufactured anisotropic and isotropic smooth shell-lattice materials: simulations & experiments [J]. Journal of the Mechanics and Physics of Solids, 2019, 122: 1–26. doi: 10.1016/j.jmps.2018.08.022
[6]	CHEN Z Y, XIE Y M, WU X, et al. On hybrid cellular materials based on triply periodic minimal surfaces with extreme mechanical properties [J]. Materials & Design, 2019, 183: 108109. doi: 10.1016/j.matdes.2019.108109
[7]	MASKERY I, ASHCROFT I A. The deformation and elastic anisotropy of a new gyroid-based honeycomb made by laser sintering [J]. Additive Manufacturing, 2020, 36: 101548. doi: 10.1016/j.addma.2020.101548
[8]	YÁNEZ A, CUADRADO A, MARTEL O, et al. Gyroid porous titanium structures: a versatile solution to be used as scaffolds in bone defect reconstruction [J]. Materials & Design, 2018, 140: 21–29. doi: 10.1016/j.matdes.2017.11.050
[9]	HARUN W S W, KAMARIAH M S I N, MUHAMAD N, et al. A review of powder additive manufacturing processes for metallic biomaterials [J]. Powder Technology, 2018, 327: 128–151. doi: 10.1016/j.powtec.2017.12.058
[10]	BAI L, GONG C, CHEN X H, et al. Mechanical properties and energy absorption capabilities of functionally graded lattice structures: experiments and simulations [J]. International Journal of Mechanical Sciences, 2020, 182: 105735. doi: 10.1016/j.ijmecsci.2020.105735
[11]	NOVAK N, KRSTULOVIĆ-OPARA L, REN Z R, et al. Compression and shear behaviour of graded chiral auxetic structures [J]. Mechanics of Materials, 2020, 148: 103524. doi: 10.1016/j.mechmat.2020.103524
[12]	XIAO L J, SONG W D, XU X. Experimental study on the collapse behavior of graded Ti-6Al-4V micro-lattice structures printed by selective laser melting under high speed impact [J]. Thin-Walled Structures, 2020, 155: 106970. doi: 10.1016/j.tws.2020.106970
[13]	刘宇, 郝琪, 田钰楠, 等. 负泊松比梯度蜂窝结构研究 [J]. 湖北汽车工业学院学报, 2020, 34(3): 64–68, 74. doi: 10.3969/j.issn.1008-5483.2020.03.013 LIU Y, HAO Q, TIAN Y N, et al. Study on structure of negative poisson’s ratio gradient honeycomb [J]. Journal of Hubei University of Automotive Technology, 2020, 34(3): 64–68, 74. doi: 10.3969/j.issn.1008-5483.2020.03.013
[14]	张权, 高松林, 杜志鹏, 等. 星形梯度负泊松比蜂窝结构面内冲击动态响应 [J]. 武汉理工大学学报(交通科学与工程版), 2020, 44(5): 886–891. doi: 10.3963/j.issn.2095-3844.2020.05.023 ZHANG Q, GAO S L, DU Z P, et al. In-plane impact dynamic response of star-shaped gradient negative poisson's ratio honeycomb structure [J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2020, 44(5): 886–891. doi: 10.3963/j.issn.2095-3844.2020.05.023
[15]	XIAO L J, SONG W D. Additively-manufactured functionally graded Ti-6Al-4V lattice structures with high strength under static and dynamic loading: experiments [J]. International Journal of Impact Engineering, 2018, 111: 255–272. doi: 10.1016/j.ijimpeng.2017.09.018
[16]	ZHAO M, ZHANG D Z, LIU F, et al. Mechanical and energy absorption characteristics of additively manufactured functionally graded sheet lattice structures with minimal surfaces [J]. International Journal of Mechanical Sciences, 2020, 167: 105262. doi: 10.1016/j.ijmecsci.2019.105262
[17]	MASKERY I, STURM L, AREMU A O, et al. Insights into the mechanical properties of several triply periodic minimal surface lattice structures made by polymer additive manufacturing [J]. Polymer, 2018, 152: 62–71. doi: 10.1016/j.polymer.2017.11.049
[18]	ZHANG L, FEIH S, DAYNES S, et al. Energy absorption characteristics of metallic triply periodic minimal surface sheet structures under compressive loading [J]. Additive Manufacturing, 2018, 23: 505–515. doi: 10.1016/j.addma.2018.08.007
[19]	MCKOWN S, SHEN Y, BROOKES W K, et al. The quasi-static and blast loading response of lattice structures [J]. International Journal of Impact Engineering, 2008, 35(8): 795–810. doi: 10.1016/j.ijimpeng.2007.10.005
[20]	冯根柱, 于博丽, 李世强, 等. 多层级夹芯结构的变形与能量吸收 [J]. 高压物理学报, 2019, 33(5): 055902. doi: 10.11858/gywlxb.20180707 FENG G Z, YU B L, LI S Q, et al. Deformation and energy absorption of multi-hierarchical sandwich structures [J]. Chinese Journal of High Pressure Physics, 2019, 33(5): 055902. doi: 10.11858/gywlxb.20180707
[21]	XIAO L J, LI S, SONG W D, et al. Process-induced geometric defect sensitivity of Ti-6Al-4V lattice structures with different mesoscopic topologies fabricated by electron beam melting [J]. Materials Science and Engineering: A, 2020, 778: 139092. doi: 10.1016/j.msea.2020.139092
[22]	ZHENG Z J, LIU Y D, YU J L, et al. Dynamic crushing of cellular materials: continuum-based wave models for the transitional and shock modes [J]. International Journal of Impact Engineering, 2012, 42: 66–79. doi: 10.1016/j.ijimpeng.2011.09.009

施引文献

资源附件(0)

访问统计

点击查看大图

图(11) / 表(1)

计量

文章访问数: 4317
HTML全文浏览量: 1921
PDF下载量: 77

1. 层间断裂韧性的有限元模拟
1.1 内聚力本构模型及损伤判据
1.2 I型和II型断裂韧性有限元数值模型
1.3 数值模拟与实验结果对比验证
2. 数值模拟结果分析
2.1 I型数值模拟结果分析
2.2 II型数值模拟结果分析
3. 结　论

3D打印梯度Gyroid结构的动态冲击响应

doi: 10.11858/gywlxb.20210701

作者简介: 厉 雪（1996－），女，硕士研究生，主要从事材料与结构的冲击动力学研究. E-mail：15152103981@163.com

通讯作者: 宋卫东（1975－），男，博士，教授，主要从事材料与结构的冲击动力学研究. E-mail：swdgh@bit.edu.cn

计量

出版历程

Dynamic Behavior of 3D Printed Graded Gyroid Structures under Impact Loading

1. 层间断裂韧性的有限元模拟

1.1 内聚力本构模型及损伤判据

1.2 I型和II型断裂韧性有限元数值模型

1.3 数值模拟与实验结果对比验证

2. 数值模拟结果分析

2.1 I型数值模拟结果分析

2.1.1 初始裂纹长度的影响

2.1.2 断裂韧性的影响

2.1.3 起始界面刚度的影响

2.1.4 界面强度的影响

2.1.5 黏结层厚度的影响

2.2 II型数值模拟结果分析

2.2.1 初始裂纹长度的影响

2.2.2 断裂韧性的影响

2.2.3 起始界面刚度的影响

2.2.4 界面强度的影响

2.2.5 净距的影响

2.2.6 黏结层厚度的影响

3. 结 论

计量

出版历程

目录

1. 层间断裂韧性的有限元模拟

1.1 内聚力本构模型及损伤判据

1.2 I型和II型断裂韧性有限元数值模型

1.3 数值模拟与实验结果对比验证

2. 数值模拟结果分析

2.1 I型数值模拟结果分析

2.2 II型数值模拟结果分析

3. 结 论

作者简介:
厉　雪（1996－），女，硕士研究生，主要从事材料与结构的冲击动力学研究. E-mail：15152103981@163.com

通讯作者:
宋卫东（1975－），男，博士，教授，主要从事材料与结构的冲击动力学研究. E-mail：swdgh@bit.edu.cn

3. 结　论

3. 结　论