0 引 言
近年来,受自然界生物结构的启发,仿生结构发展迅速。得益于优异的吸能能力,一批诸如蜂窝结构、蜘蛛网结构、螺旋结构等被应用于交通运输、航空航天、装甲防护等领域。
作为一种常见的多孔材料,蜂窝具有质量轻、强度高、刚度高、吸能效率高等良好的力学性能。目前,典型蜂窝结构主要有三角形蜂窝、四边形蜂窝、圆形蜂窝、负泊松比蜂窝等。Tan等
[1]在负泊松比蜂窝的基础上,结合分形设计,以正六边形和正三角形作为子结构,建立了两种分形负泊松比蜂窝结构RHH、RHT,在保持负泊松比特性的同时极大提高了结构的比吸能。于征磊等
[2]设计了简化的仿直蜂窝、仿斜蜂窝结构,研究了5种载荷工况下的结构承载特性。Yang等
[3]在三角形、四边形、六边形蜂窝等结构的基础上,提出了一种新型马蹄形铝蜂窝,有效提高了蜂窝结构的耐撞性能。Liang等
[4]以六边形蜂窝为基础结构,建立了7种不同尺寸参数的自相似分级蜂窝结构,试验结果表明:分级蜂窝结构具有变形模式稳定、多级吸能等优点,提出的参数化设计策略可为特定工程需求下结构的吸能过程调控提供参考。同样地,Zhang等
[5]基于蜘蛛网的生物几何六边形构造,在3种低阶分层结构的中心引入正六边形,建立了3种高阶分层结构,试验结果表明:分层结构破碎力波动较小,具有理想的变形模式和能量吸收能力。自然界中竹子具有良好的轴向力学性能,受竹子管壁微观结构的启发,Zhang等
[6]提出了一种具有结构层次的管状蜂窝结构。从实验、数值和理论3个方面对该结构的吸能性能进行了系统的研究,发现结构各项吸能指标受设计参数影响较大,平衡设计参数与结构吸能之间的关系可以获得较高吸能的结构。
自然界中,螺旋结构(Bouligand structure, BS)是一种奇特的轮廓,这类结构广泛存在于鱼类
[7,8]、虾类
[9,10]、螃蟹类
[11,12]、蝎子
[13,14]、甲虫
[15,16]等生物中,且具有高刚度、高强度、高延展性和高韧性的特点。巨骨舌鱼鳞片的胶原纤维片层中每层纤维取向是一致的,但相邻片层之间是错位的,构成了Bouligand结构,根据加载环境,这类结构允许层间片层重新定向,以适应不同方向的负载,有效地增强了鳞片的延展性和韧性
[17]。Chen等
[18]利用环保的硅基纳米纤维和生物聚合物构建了非连续的Bouligand结构纳米复合材料,与天然Bouligand结构材料相比,该类材料具有优异的层间耦合能力、较高的应力传递效率和抵抗弯曲裂纹扩展的能力。从螳螂虾螯棒微观Bouligand结构中获得灵感,Zhang等
[19]制备了不同层间夹角的仿生螺旋结构,并对其力学性能及承载时的损伤破坏机理进行了研究,结果表明:螺旋层状结构的承载能力和能量吸收特性与层间螺旋角密切有关,在一定范围内,随层间角的减小,螺旋结构的强度先增大后减小,而当螺旋角为12°时,吸能最大。Yang等
[20]提出了一种基于甲壳类螺旋结构断裂特性和增韧机理的理论建模方法,结合仿生设计,建立了高碳纤维增强聚合物交叉层合板,对其断裂性能进行研究,结果表明:螺旋纤维的取向对结构的断裂性能有重要影响,驱动倾斜扭曲裂纹的扩展将会消耗更多能量。同样地,结合甲壳类动物外骨骼观察到的独特螺旋形形态,Cheng等
[21]设计了3种仿生螺旋结构和一种基础结构,对4个结构的力学性能进行评估,发现仿生结构的抗弯刚度和强度、横向剪切模量等均优于常规基础结构,自然启发螺旋结构具有很大的应用潜力。
本文以螳螂虾虾螯微观螺旋结构为研究对象,选取旋转角度30°,建立了4种仿生螺旋层状结构:仿生六边形螺旋结构(FZ)、仿生圆形螺旋结构(YZ)、仿生四边形螺旋结构(SZ)、仿生负泊松比螺旋结构(FBZ)以及4种层状对比结构:六边形结构(F)、圆形结构(Y)、四边形结构(S)、负泊松比结构(FB)。以热塑性聚氨酯材料为基材,运用3D打印技术进行结构制备,结合试验和数值模拟对结构力学性能进行综合分析,旨在寻求高比吸能、耐撞性优良的轻量化结构设计方案。
1 仿生结构设计及材料分析
1.1 仿生螺旋结构设计
作为一种肉食性海洋甲壳类动物,螳螂虾利用其锤状螯可轻易击碎软体动物外壳进行捕食,速度接近23 m/s,且能承受上万次这样的冲击。相关研究表明
[9,10]:螳螂虾锤状螯棒结构由多相结合的取向晶体羟基磷灰石、无定形磷酸钙和碳酸盐的多相复合材料组成,其中由矿化几丁质纤维构成的螺旋(Bouligand)结构是该材料具有轻松抵御高能量冲击及阻碍裂纹扩展等性能的关键因素。如
图1(d)
[22]所示,进一步进行微观结构提取,发现螳螂虾螯微观结构是由螺旋形层状结构堆叠而成,每层间隔分布多个圆管形结构。
受BS启发,本文以六边形、圆形、四边形、负泊松比4种类型胞元作为层状基本结构,进行结构设计。如
图1(e)所示,以整体旋转六边形结构为例,将相邻层间角度设置为30°,层数为0~8层,角度从0°~210°。最终,建立了4种仿生螺旋结构FZ、YZ、SZ、FBZ以及4种对比结构F、Y、S、FB。如
图2所示,上述结构均设计为直径33 mm、高度41 mm的圆柱,其中各层高度为5 mm,4种胞元基本尺寸设置为4 mm,相邻层以及相邻胞元之间厚度为1 mm。
1.2 材料性能测试
以热塑性聚氨酯材料为基材,运用3D打印技术进行结构制备,打印设备型号为ZS500,每种模型均制备两个样件,共制备了16个压缩样件以及2个拉伸样件,打印工艺参数如下:切片厚度为0.15 mm,填充间距为0.2 mm,加工温度为95 ℃,激光功率为44 W,激光速度为4 m/s。
采用万能拉伸试验机DDL100对结构压缩性能和拉伸性能分别进行测试。试验设置为准静态压缩试验和准静态拉伸试验,试验速度分别设置为2、5 mm/min。为获取材料参数,对拉伸样件进行试验,通过固定样件一端,另一端以指定速度移动,直至拉伸样件断裂,具体拉伸曲线及试验设置如
图3所示。
具体的热塑性聚氨酯材料参数如下:密度为1.1 g/cm3,弹性模量为30 MPa,屈服强度8 MPa,泊松比为0.33。
1.3 有限元模型建立
本文利用有限元软件LSDYNA进行显式计算分析。为了复刻准静态压缩试验设置,建立
图4(a)所示的有限元分析模型。压缩样件底部被固定在刚性墙上,上端设置有匀速向下移动的刚性平板,压缩距离设置为12 mm。压缩样件与刚性平板之间设置为面面接触,以防止网格穿透,静摩擦因数与动摩擦因数均设置为0.3,以模拟真实压缩试验。选用四面体网格对结构进行划分,以SZ结构为例,选取0.8、1、1.2、1.4 mm 4种网格尺寸进行无关性分析,如
图4(b)所示,当网格尺寸由1.4 mm到0.8 mm,结构的比吸能开始下降,并在1 mm时趋于收敛,考虑到计算成本,最终确定网格尺寸为1 mm进行后续数值模拟。
1.4 分析评价指标
为了评价螺旋结构和对比结构的力学性能,采用两种分析指标对其进行分析,分别为吸能(Energy absorption,EA)、比吸能(Specific energy absorption,SEA)。
式中:F(x)为破碎力;x为瞬时破碎位移;d为有效破碎位移。
式中:M为结构的总质量。SEA值越高,能量吸收越好。
2 仿生螺旋结构性能分析
2.1 仿生螺旋结构力学特性分析
为了研究所建立结构样件的力学特性,对上述8组模型进行准静态压缩试验,得到每个模型样件的载荷力位移曲线。准静态压缩试验及有限元分析载荷力位移曲线如
图5所示。
如
图5所示,仿真曲线与试验曲线拟合较好。在4种类型的结构中,螺旋结构的载荷力位移曲线均明显高于对比结构,这表明螺旋结构具有较好的力学特性。如
图5(a)所示,FZ结构从压缩初始就表现出比F结构较大的载荷力。如
图5(b)所示,对于Y、YZ结构,试验曲线在达到峰值之后出现突然下降,这是由于结构压缩过程出现局部断裂,可能是结构打印过程存在缺陷,断裂区域位置成型质量较差,而仿真曲线没有考虑压缩过程中的断裂,仿真曲线与试验曲线表现出相似的趋势。如
图5(c)所示,S结构在达到峰值之后,载荷力开始下降,而SZ结构出现多个峰值,结合变形模式分析,这是由于SZ结构变形过程出现层间依次坍塌,导致载荷力位移曲线出现多个峰值,整体上SZ结构力学特性优于S结构。如
图5(d)所示,FB结构达到峰值后,出现平台期,载荷力平稳波动直至出现小幅上升,而FBZ的“平台期”出现斜率较小的力增长,随着压缩过程的进行,力继续增长,FBZ结构在压缩位移达到3 mm附近后,载荷力开始明显高于FB结构。整体来说,与对比结构相比,仿生螺旋结构具有较大的载荷力。
如
图6所示,分析了4种仿生螺旋结构及4种对比结构的吸能和比吸能情况,可以看出试验和仿真所得结果相差不大。如
图6(a)所示,F结构的吸能为6 652.14 mJ,比吸能为322.76 mJ/g,FZ结构的吸能为11 568.80 mJ,比吸能为531.65 mJ/g,与F结构相比,FZ结构比吸能提升64.72%,这表明螺旋结构的吸能特性得到了极大提升。如
图6(b)所示,Y结构的吸能为4 544.65 mJ,比吸能为240.84 mJ/g,YZ结构的吸能为6 348.38 mJ,比吸能为318.85 mJ/g,与Y结构相比,YZ结构比吸能提升32.39%。同样地,由
图6(c)(d)可知,SZ结构比吸能为135.33 mJ/g,FBZ结构比吸能为185.39 mJ/g,与S、FB结构相比,SZ、FBZ结构比吸能分别提升25.14%、55.84%,以上分析结果均证明了仿生螺旋结构设计的可行性。与对应对比结构相比,FZ结构比吸能提升最大,其次为FBZ、YZ、SZ,其中 FZ结构比吸能最大,为4种类型结构中的最优结构。
2.2 仿生螺旋结构变形分析
为了分析结构的压缩变形情况,本文记录了所有结构样件在准静态压缩试验过程中(
ɛ=0.1、0.2、0.3)的变形情况,并与有限元仿真的Von Mises云图进行了对照,仿生结构样件变形如
图7所示。
如
图7所示,F结构在中间位置出现相互正交的变形区域;FZ结构遵循整体屈曲变形,其压缩过程中也伴随着局部正交变形区域的出现。Y结构在
ɛ=0.2时在上部首先变形的区域出现断裂,之后在其对称下部区域产生变形,仿真结构变形与之类似,但未考虑断裂;YZ结构整体变形,随着压缩过程的进行,向一侧坍塌变形。S结构在变形初始阶段向右侧坍塌变形,随后结构上下对称变形,最终结构上下表面出现翘曲,仿真分析因结构底部区域固定,仅上部出现翘曲现象,但二者变形较为相似;SZ结构变形过程出现层间逐次变形,与试验相比,仿真变形层位置有所不同,这可能与结构实际打印成型质量有关,由于层间依次变形,导致载荷力-位移曲线出现多个峰值。FB结构变形与上述S结构较为相似,即向一侧坍塌,上下区域对称变形;FBZ结构整体屈曲变形。综上,与对比结构相比,得益于层间螺旋堆叠排列,仿生螺旋结构变形过程中避免了一些特殊变形的出现:Y、S、FB中的上下对称变形区域以及F结构的大面积正交变形区域,仿生螺旋结构变形较为均匀,有效增强了结构的承载能力。对于上述4种类型结构,在整个变形时间段内,仿真与试验变形基本吻合,验证了本文有限元模型的准确性。
通过上述试验结果发现,仿生螺旋结构具有较好的承载特性,为设计具有较高吸能和抗冲击特性的仿生结构提供了设计参考。
2.3 仿生螺旋结构变角度分析
为了对比仿生螺旋结构在不同层间角度下的性能,将每个仿生螺旋结构设置3个梯度(15°、30°、45°)进行仿真实验,压缩位移均设置为12 mm。
图8为4种仿生螺旋结构在3种层间螺旋角度下的载荷力位移曲线。如
图8(a)所示,3个结构在压缩初始阶段表现出相似的载荷力,随后FZ-45°结构表现出最大的载荷力。如
图8(b)所示,3个仿生圆形螺旋结构表现出相似的载荷力位移曲线。如
图8(c)所示,仿生四边形螺旋结构整体载荷力大小变化相似,然而由于层间螺旋角度的不同,3个结构峰值出现的时间不同。如
图8(d)所示,FBZ结构在较低层间螺旋角度下表现出较高的载荷力。
图9分析了上述结构的吸能情况。在3个旋转角度下,FZ、YZ、SZ、FBZ结构最佳吸能角度分别为45°、30°、45°和15°。整体来说,FZ结构吸能效果最好,其次为YZ、FBZ、SZ。
上述试验发现,仿生螺旋结构在不同层间角度下的吸能具有一定差异,未来会针对上述4种构型进行更多旋转角度梯度设计,旨在获得性能更为优异的结构。
3 结 论
(1)4种类型的仿生螺旋结构均表现出优于4种对比结构的吸能性能,这表明层间螺旋堆叠分布有效增强了结构的吸能特性,结构变形破坏将会吸收更多的能量。
(2)FZ结构性能最好,比吸能最大,可达531.65 mJ/g,其次为YZ、FBZ、SZ结构,且与对应对比结构相比,4种仿生螺旋结构比吸能分别提升64.72%、32.39%、55.84%、25.14%。
(3)4种仿生螺旋结构在不同层间螺旋角度下的吸能具有一定差异,FZ、YZ、SZ、FBZ结构最佳吸能角度分别为45°、30°、45°和15°。本文仿生螺旋结构设计思路可为建立轻质、高强、高比吸能仿生结构提供参考。