An experimental technique for medium strain-rate loading by a progressive cam
-
摘要: 研制了一种可以实现多次加载的凸轮递进式中应变率压缩实验系统。该实验装置采用伺服电机驱动蓄能飞轮转动,后蓄能飞轮带动加载凸轮压缩加载杆的方法,实现对试样中应变率的压缩;同时在一级压缩即将结束时步进电机迅速推动蓄能飞轮贴近加载凸轮,实现多级压缩。试样的动态压缩载荷通过两侧杆上粘贴的应变片所记录的应变信号得到;试样变形过程通过激光干涉测速系统测得的试样两侧杆端的运动速度信号得到。以纸蜂窝试样为例,基于研制的中应变率实验系统,并结合高速摄影图片,研究了厚度10 mm、直径14.5 mm的纸蜂窝试样在应变率3.5 s−1下的动态压缩力学性能,得到了单级压缩和两级压缩过程中纸蜂窝试样的应力-应变曲线和变形过程,并讨论了该实验系统的可靠性。此实验系统可以实现多级递进式中应变率加载;纸蜂窝试样在中等应变率下的峰值强度和平台应力对高应变率下的动态压缩实验数据和低应变率下的准静态实验数据进行了较好地衔接;试样的失效模式主要为准弹性变形后的外壁屈曲和面内剪切。Abstract: A medium strain rate compression experimental system based on a progressive cam was developed to realize multiple medium strain rate loading. The developed experimental system uses the servo motor to drive the energy storage flywheel to rotate at a certain speed, and when the clutch is started, the energy storage flywheel can drive the loading cam to rotate. The loading cam pushes the loading guide bar and the input bar to compress the sample. When the loading cam rotates one circle, a single medium strain-rate compression is completed. At the same time, when the first stage compression is about to end, the stepper motor rapidly pushes the energy storage flywheel close to the loading cam for the next compression, and the cycles repeat to achieve multiple medium strain rate compression. The load and deformation of the material during compression were measured by strain gauges and a velocity interferometer system for any reflector (VISAR), respectively. The strain gauges were affixed to the input bar and the support bar, respectively. The strain signals of the bars during compression were recorded by the strain gauges and the forces exerted on the sample were obtained based on these strain signals. Two fiber optic probes of the VISAR system were used to measure the velocities of the input bar and the support bar during compression. Based on the two velocity curves measured, the velocity difference curve between the two ends of the sample was obtained, and then the deformation of the sample was gained by integrating the velocity difference. The stress-strain curves were obtained from the load- and deformation-time curves. Taking the paper honeycomb sample as an example, the reliability of the developed medium strain rate experimental system was discussed based on the high-speed images. The dynamic compressive mechanical properties of the paper honeycomb samples with the thickness of 10 mm and the diameter of 14.5 mm at the strain rate of 3.5 s−1 were studied. The stress-strain curves and deformation processes of the paper honeycomb samples during single compression and double compression were obtained. The experimental system could realize multistage progressive medium strain rate loading. The peak strength and plateau stress of the paper honeycomb samples at medium strain rates well connect the dynamic compression results at high strain rates with the quasi-static compression results at low strain rates. The failure modes of the samples are mainly out-of-plane buckling and in-plane shear after quasi-elastic deformation.
-
Key words:
- medium strain rate /
- paper honeycomb /
- mechanical performance /
- compression experiment
-
薄壁结构作为吸能装置被广泛应用于吸能缓冲、抗冲击领域,近年来,折纸图案在管状薄壁结构上的应用越来越受到关注[1]。研究发现,如果在管件表面预制不同模式的折痕,可调控在轴向外力作用下结构的屈曲模式和能量吸收性能。在圆管表面引入波纹,可以产生较低的初始峰值力和更均匀的压缩过程[2]。Zhang等[3]通过数值模拟研究了两种具有金字塔形状的方形管,发现其能量吸收比传统方管分别提高了(15~33)%和(54~93)%。Zhao等[4]在折纸管中引入具有有限厚度和长度的铰链区域,并利用3D打印技术进行设计,这种设计充分展示了3D打印技术为增强折纸结构的承载能力和可控的变形模式提供了设计灵活性与制造效率。Ma等[5]提出了一种新的装置,该装置是根据一种特别设计的折纸图案,在方形管的表面预折而成。和传统方管相比,在该设计中,初始峰值力降低了20%以上,能量吸收增加了50%以上。Song等[6]应用了一种折纸图案做成方形管,使初始峰值力最小化。Zhou等[7]提出一种梯形折纸薄壁结构,其表面根据可展开的折纸图案预先折叠,并在方管表面引入一种梯形折痕以降低峰值力,将设计中最优的梯形折纸管与传统的等质量方管和八边形管比较,发现梯形折纸薄壁结构在吸能方面优于传统薄壁结构。周昳鸣等[8]针对这种折纹管进行了优化设计。郝文乾等[9]采用LS-DYNA软件研究了6种折纹管在轴向冲击载荷下的屈曲模式和能量吸收性能,发现折痕的引入有效地降低了初始峰值力,并且部分折纹管的压缩力效率和比总体效率高于方管的。Hou等[10]和Ma等[11]基于风筝形折纸模型,提出了一种截面为风筝形的管状结构,与传统方管相比,其峰值载荷降低了67.2%,能量吸收提升了36%。与传统的方形管相比,在保证能量吸收能力一定的情况下,具有初始折痕的薄壁管具有较低的初始峰值力和更均匀的压缩载荷。
本文中,将折纸艺术与3D打印技术相结合,设计并制备了面外方向具有初始折叠角度的弧形折纸薄壁管件单胞及其阵列结构,基于准静态轴向压缩实验,利用ABAQUS对弧形折痕薄壁管件在轴向压缩下的力学行为和吸能特性进行了数值模拟分析,探讨了结构的变形模式,分析了预折角的大小、单胞、面内阵列个数以及不同的冲击速度对其压缩行为的影响。由于压缩行程、质量是吸能结构设计与评价的主要参量,所以本文中将高度相等、质量近似相等(结构最大质量差为0.07 g,最小质量差为0.01 g)的方管及其阵列结构进行对比,发现弧形折痕的引入能够有效地降低初始压溃载荷峰值并减小冲击载荷作用下冲击力位移曲线的振荡幅值,可以作为一种变形可控的吸能结构。
1. 几何模型
本文中所研究的模型主要由一个矩形薄板沿着预定的折痕折叠而成,如图1(a)所示。其中,DE长度记为c,AD=DG=GK= h0,折痕倾角∠ABE记为α,折痕BE=CF=HL=IM=a,可得:
h0=asinα (1) 将基本单元沿折痕BE、EH、HL折叠成图1(b)所示的几何形状,其中,倾斜翼缘BCFE和水平面的夹角β、腹板DEHG和腹板FEHI的夹角ξ以及BE与水平面的夹角θ由以下几何关系进行计算:
sinβ=sinθsinα (2) cosξ2=cosαcosθ (3) 因此,几何模型可根据c、a、α、ξ来确定。为保证折纸管单胞的可周期阵列性,设计折痕时,令EF=2DE,OF=4DE,腹板DEHG和腹板FEHI的夹角ξ=60°。将2个沿折痕折叠成的单元对称连接,得到弧形折纸薄壁管件,如图1(c)所示。长、宽、高分别用l、w、h表示(图1中弧形实折线(CFIM)为向外的峰线,弧形虚折线为向内的谷线),其中:
l=2(2csinξ+asinαcosβ),w=OF=4c,h=asinα+2asinθ (4) 本文中所研究的模型均使用3D打印机打印制造,在建立3D打印模型时,连接面DEHG与其他竖直腹板的厚度一致。
2. 实验研究
本研究中所有样品均采用Volunex Mars Plus+ 3D打印机打印。使用的材料为PolyMaxTM PLA,成型后材料密度为1 200 kg/m3,泊松比为0.35,弹性模量为1.97 GPa。利用SHIMADZU AG-25TB/SFL-5AG材料试验机,参照拉伸测试规范ASTM D638-10测试了拉伸试样在准静态(应变率为10−4 s−1)拉伸时的力学参数,初始屈服强度为40 MPa,断裂应变约为0.3。拉伸试样的几何参数及真应力应变曲线如图2(a)所示。
3D打印试件力学性能与其打印工艺关系密切,对于熔融堆积法打印工艺,试件拉伸力学响应与打印角度相关性较大。本文研究的管状结构打印模式为沿着结构轴线向上铺层打印,为保证基体材料力学性能测试的可靠性,试件打印选取沿长度方向铺层打印,如图2(b)所示。
根据前文所述几何模型,利用Volunex Mars Plus+3D打印机制造表观尺寸(l×w×h)为80 mm×80 mm×60 mm、预折角α=60°、壁厚τ=1 mm、在面内方向阵列个数n为2×2的薄壁管状结构,在ABAQUS中建立了相应的有限元模型。实验中,采用SHIMADZU AG-25TB/SFL-5AG压缩试验机,压缩速度为1 mm/min;压头为钢制平板,与结构自由接触。结构在准静态压缩载荷下的实验和数值模拟对比结果如图3所示,两者吻合较好。数值模拟的模型尺寸与实验相同,有关单元类型、材料模型及接触关系的详细信息见第4节。由于数值模拟中未考虑摩擦效应的影响,结构的压缩力略低于实验值。
从图3可以看出,结构变形可定性地分为4个阶段。初始压溃阶段:当结构受到初始压缩载荷时,压缩力迅速上升至峰值载荷,约为7.5 kN。之后,结构进入预折角塑性旋转阶段:预折角开始发生塑性旋转变形,同时伴随有局部破坏,此时结构压缩力呈非线性下降趋势。随着预折角的进一步旋转,结构压缩力逐渐上升,当腹板两侧倾斜部分旋转至水平位置时,压缩力主要有腹板承担,在力-位移曲线上出现第2个峰值载荷,约为10 kN,这个过程为腹板塑性屈曲阶段。当腹板受力发生屈曲后,压缩力下降,随后整个结构逐步进入压缩密实化阶段。
3. 数值模拟
针对不同几何参数的弧形折纸模式薄壁结构开展准静态和动态压缩数值模拟。弧形折纸薄壁管单胞根据其预折角α的不同取3种,分别记为A50°、A60°、A70°,与之对比的方管记为SQU。由于常用的多孔夹芯抗冲击防护结构芯层厚度在15~40 mm之间[12-13],因此本节中所涉及的弧形折纸模式单胞结构(实验组)和方管单胞模型(对照组)的相关参数如表1所示。
表 1 模型参数Table 1. The parameters of the models模型 单胞质量/g h/mm w/mm l/mm τ/mm c/mm a/mm h0/mm ξ /(°) α/(°) θ/(°) β/(°) A50° 8.65 30 40 45.23 1.10 10 14.26 10.92 60 50 42 61 A60° 8.66 30 40 41.57 1.20 10 11.98 10.37 60 60 55 71 A70° 8.59 30 40 38.86 1.25 10 10.79 10.14 60 70 67 78 SQU 8.66 30 40 41.89 1.47 3.1 有限元模型及算法验证
结构单胞有限元模型如图4所示,折纸管采用S4R壳单元,上下板为刚体。对于2×2阵列模型,在建模时直接把所有的单胞合并为一个整体。经过网格敏感性和有限元算法验证后(图5~6),数值模拟中的模型网格尺寸确定为1 mm,上下板和折纸管采用自接触,忽略摩擦的影响,上板采用1 m/s的压缩速度来模拟准静态加载过程(总动能必须比压缩过程的总内能小得多,压缩力位移响应必须与所施加的速度无关[14-15])。
从图5可以看出:折纸管在①处产生初始压溃应力峰值,随后进入预折角塑性旋转阶段(①~②),该阶段称为第1平台段;预折角塑性旋转阶段结束后,应力值开始上升,此时主要承载部分为折纸管的竖直腹板,在折纸管出现第2个压溃应力峰值(③)之后,折纸管进入腹板塑性屈曲阶段(③~④)该阶段称为第2平台段;当竖直腹板完全压溃之后(④之后),整个模型进入密实化阶段。
3.2 参数分析
3.2.1 预折角及面内方向阵列个数的影响
本文中所涉及的应力应变曲线均为名义应力-名义应变曲线。名义应力σ为压板与薄壁管件接触面上的力与模型初始横截面面积之比,名义应变ε定义为模型沿压缩方向的整体压缩与原长度之比,即:
σ=flw,ε=Δlh (5) 式中:f为压缩载荷,Δl为压缩位移。
弧形折纸薄壁管和方管的单胞在准静态压缩作用下的应力-应变曲线如图7所示。由图7可知,随着预折角α的增大,折纸管在初始压溃应力峰值越大,折纸管的第一平台段应力应变曲线的变化趋势一致且应力值的大小基本保持恒定,在ε=0.62时,模型的主要承载部分由预折部分变为竖直的方管,折纸管产生第2个应力峰值,进入腹板塑性屈曲阶段。由于竖直管部分发生坍塌失效后,模型失去承载能力,折纸管的应力应变曲线会出现下降趋势。在此过程中,A50°单胞模型发生整体坍塌,模型承载能力急剧下降;A60°、A70°单胞模型产生渐进压溃模式,在压缩过程中,每个竖直板发生坍塌失效都会对应产生一个峰值应力,因此,A60°、A70°单胞模型的应力应变曲线会发生波动。方管的单胞模型在模拟中的初始压溃应力峰值最大。与只具有一个平台阶段的方管相比较,折纸管的第1平台段应力值与方管平台阶段的应力值相近,折纸管的第2平台段应力值高于方管平台阶段的应力值。
弧形折纸薄壁管和方管面内方向阵列个数n为2×2模型的应力-应变曲线如图8所示(方管阵列采用密排结构,相邻单胞方管间紧密接触,接触位置的壁厚为2倍的单壁厚(τ),在有限元模型中设置该处壳单元壁厚为2τ),腹板塑性屈曲阶段,A50°、A60°阵列模型在出现第2个峰值应力之后发生整体坍塌,应力急剧下降,A70°阵列模型竖直管发生整体失稳,第2平台段应力值变小,由于方管阵列模型中间连接部分腹板壁厚为侧边壁厚的两倍,在压缩过程中,当中间腹板发生整体坍塌时,应力有明显增高,屈曲失稳后失去承载能力,应力迅速下降。所有的弧形折纸薄壁管的第2平台段应力值均大于方管。图9为A60°单胞、面内方向阵列个数n为2×2、4×4的阵列模型在准静态作用下的应力-应变曲线,由图可知,2×2、4×4阵列模型的应力应变曲线基本相同,增加面内方向阵列个数对折纸管的应力应变关系没有明显的影响。由于阵列后胞元之间存在一定的横向约束,多胞阵列模型的第一平台段应力明显高于单胞模型,且随着单胞阵列个数的增加,折纸管的第2平台段应力值明显提高。
3.2.2 冲击速度的影响
图10为折纸管面内方向阵列个数n=2×2的阵列模型在冲击速度为10、20 m/s载荷作用下的应力-应变曲线。从图10可以看出,折纸管在冲击速度为10 m/s载荷作用下的初始压溃应力峰值均低于方管的,预折角的变化对折纸管的第1平台段应力-应变曲线变化趋势及其应力值影响较小。在腹板塑性屈曲阶段,所有折纸管在出现第2个峰值应力之后发生整体坍塌,应力急剧下降,但第2平台段应力值均大于方管的。当冲击速度为20 m/s时,A70°阵列模型在压缩过程中表现出逐层压溃模式,在ε=0.42时,折纸管上部预折部分首先达到密实化,产生第2个应力峰值,在ε=0.65时,折纸管预折部分全部密实化,产生第3个应力峰值,折纸管的第1、2平台段应力值均小于方管的。
3.2.3 弧形折纸薄壁管和方管的吸能效率
图11~12为方管和弧形折纸薄壁管的单胞、n=2×2阵列模型在准静态、冲击速度为10、20 m/s的比吸能对比。从图11~12可以看出:在准静态压缩作用下,对于单胞模型,A70°在压缩位移约22 mm之后,比吸能高于方管的,面内方向阵列个数n=2×2的阵列模型中,方管比吸能优于折纸管;在冲击速度为10、20 m/s的载荷作用下,n=2×2的阵列模型中,弧形折纸薄壁管中A70°模型具有较高的比吸能,但方管的比吸能均优于弧形折纸薄壁管。
对准静态,冲击速度为10、20 m/s的单胞、面内方向阵列个数n为2×2模型的数值模拟采用压缩力效率η和比总体效率δ作为薄壁管吸能特性的两种评价指标,压缩效率越大,模型的吸能效果越好。压缩力效率为平均载荷与初始压溃载荷峰值的比值,平均载荷是结构吸收的总能量与最后压缩量之比,比总体效率是单位长度下管的比吸能与初始压溃载荷峰值的比值,所有模型的初始压溃载荷峰值如表2所示。
表 2 不同模型的初始压溃载荷峰值Table 2. Initial peak crushing loads of different models冲击速度/(m·s−1) n 载荷/kN A50° A60° A70° SQU 1 1×1 0.23 0.65 1.72 6.59 2×2 1.72 3.91 8.47 32.46 10 2×2 4.92 9.06 13.22 42.84 20 2×2 10.47 19.77 23.72 46.14 图13~14为弧形折纸薄壁管和方管的压缩力效率-压缩位移曲线,从图中可以看出,在准静态压缩作用下,折纸管单胞及面内方向阵列个数n=2×2阵列模型的压缩力效率优于方管。在冲击速度为10 m/s的载荷作用下,A50°阵列模型具有较高的压缩力效率,且所有的弧形折纸薄壁管的压缩力效率均高于方管,如图13所示。在冲击速度为20 m/s的载荷作用下,当压缩位移d<19.73 mm时,方管阵列模型的压缩力效率高于折纸管,当压缩位移d>19.73 mm时,A50°折纸管阵列模型的压缩力效率高于方管,如图14所示。
图15~16为弧形折纸薄壁管和方管的比总体效率的比较。从图15可以看出,在准静态压缩及冲击速度为10 m/s的载荷作用下,折纸管单胞及面内方向阵列个数n=2×2阵列模型的比总体效率均优于方管。从图16可以看出,在冲击速度为20 m/s的载荷作用下,当压缩位移d<19.67 mm时,方管的阵列模型比总体效率高于折纸管,当压缩位移d>19.67 mm时,A50°折纸管阵列模型的比总体效率高于方管。
4. 结 论
利用3D打印技术制备了弧形折纸薄壁管件,基于准静态轴向压缩实验,并利用ABAQUS有限元软件对弧形折纸薄壁管件在准静态,冲击速度10、20 m/s下的压缩行为进行了数值模拟,探讨了结构的变形模式,分析了不同预折角以及单胞的阵列个数对其变形过程及承载能力的影响,并与高度相等、质量近似相等的方管作对比,得出以下结论。
(1)薄壁管件的变形过程可分为4个阶段,分别为初始压溃阶段、预折角塑性旋转阶段、腹板塑性屈曲阶段和完全压溃密实化阶段。
(2)设计的弧形折痕能够有效降低初始压溃阶段的载荷峰值,减小压缩载荷-位移曲线在预折角塑性旋转阶段、腹板塑性屈曲阶段的振荡幅值,使结构压缩过程可控、按照预定的弧形折痕发展。
(3)在准静态压缩下:对于单胞模型,A70°模型的比吸能优于方管的;对于多胞管件阵列模型,方管的比吸能均优于折纸管的;折纸管的压缩力效率和比总体效率均优于方管的,且A50°模型的压缩力效率和比总体效率最高。
(4)在本文涉及的动态冲击载荷下,分别通过比吸能、压缩力效率和比总体效率3个评价指标对比了不同结构的吸能特性:阵列方管的比吸能均优于阵列折纸管的;当冲击速度为10 m/s时,折纸管的压缩力效率和比总体效率均优于方管的,其中A50°模型的压缩力效率和比总体效率最高;当冲击速度为20 m/s时,仅有A50°模型的压缩力效率和比总体效率优于方管的。
-
-
[1] HUANG J, XU S, YI H, et al. Size effect on the compression breakage strengths of glass particles [J]. Powder Technology, 2014, 268: 86–94. DOI: 10.1016/j.powtec.2014.08.037. [2] 薛晓, 乔禹, 王鹏飞, 等. 碳纳米管纤维的动态拉伸力学性能研究 [J]. 实验力学, 2020, 35(5): 811–819. DOI: 10.7520/1001-4888-19-213.XUE X, QIAO Y, WANG P F, et al. Dynamic tensile mechanical properties of carbon nanotube fiber [J]. Journal of Experimental Mechanics, 2020, 35(5): 811–819. DOI: 10.7520/1001-4888-19-213. [3] MIAO C H, XU S L, SONG Y P, et al. Influence of stress state on dynamic breakage of quartz glass spheres subjected to lower velocity impacting [J]. Powder Technology, 2022, 397: 117081. DOI: 10.1016/j.powtec.2021.117081. [4] 袁良柱, 苗春贺, 单俊芳, 等. 冲击下混凝土试样应变率效应和惯性效应探讨 [J]. 爆炸与冲击, 2022, 42(1): 013101. DOI: 10.11883/bzycj-2021-0114.YUAN L Z, MIAO C H, SHAN J F, et al. On strain-rate and inertia effects of concrete samples under impact [J]. Explosion and Shock Waves, 2022, 42(1): 013101. DOI: 10.11883/bzycj-2021-0114. [5] 郑监, 卢芳云. 金属梁在预应力下的冲击响应 [J]. 爆炸与冲击, 2021, 41(3): 031401. DOI: 10.11883/bzycj-2020-0328.ZHENG J, LU F Y. On impact response of a prestressed metal beam [J]. Explosion and Shock Waves, 2021, 41(3): 031401. DOI: 10.11883/bzycj-2020-0328. [6] HUH H, LIM J H, PARK S H. High speed tensile test of steel sheets for the stress-strain curve at the intermediate strain rate [J]. International Journal of Automotive Technology, 2009, 10(2): 195–204. DOI: 10.1007/s12239-009-0023-3. [7] FROUSTEY C, LAMBERT M, CHARLES J L, et al. Design of an impact loading machine based on a flywheel device: application to the fatigue resistance of the high rate pre-straining sensitivity of aluminium alloys [J]. Experimental Mechanics, 2007, 47(6): 709–721. DOI: 10.1007/s11340-007-9082-4. [8] COLE B N, STURGES J L. The flying wedge: a method for high strain rate tensile testing. Part 2: characteristics of the device [J]. International Journal of Impact Engineering, 2003, 28(8): 891–908. DOI: 10.1016/S0734-743X(03)00003-4. [9] PETITEAU J C, OTHMAN R, GUÉGAN P, et al. A drop-bar setup for the compressive testing of rubber-like materials in the intermediate strain rate range [J]. Strain, 2014, 50(6): 555–562. DOI: 10.1111/str.12113. [10] PEROGAMVROS N, MITROPOULOS T, LAMPEAS G. Drop tower adaptation for medium strain rate tensile testing [J]. Experimental Mechanics, 2016, 56(3): 419–436. DOI: 10.1007/s11340-015-0112-3. [11] WHITTINGTON W R, OPPEDAL A L, FRANCIS D K, et al. A novel intermediate strain rate testing device: the serpentine transmitted bar [J]. International Journal of Impact Engineering, 2015, 81: 1–7. DOI: 10.1016/j.ijimpeng.2015.02.009. [12] ZHU D, RAJAN S D, MOBASHER B, et al. Modal analysis of a servo-hydraulic high speed machine and its application to dynamic tensile testing at an intermediate strain rate [J]. Experimental Mechanics, 2011, 51(8): 1347–1363. DOI: 10.1007/s11340-010-9443-2. [13] OTHMAN R, GUÉGAN P, CHALLITA G, et al. A modified servo-hydraulic machine for testing at intermediate strain rates [J]. International Journal of Impact Engineering, 2009, 36(3): 460–467. DOI: 10.1016/j.ijimpeng.2008.06.003. [14] LI Z, WANG T, JIANG Y, et al. Design-oriented crushing analysis of hexagonal honeycomb core under in-plane compression [J]. Composite Structures, 2018, 187: 429–438. DOI: 10.1016/j.compstruct.2017.12.066. [15] ZHOU H, XU P, XIE S C, et al. Mechanical performance and energy absorption properties of structures combining two Nomex honeycombs [J]. Composite Structures, 2018, 185: 524–536. DOI: 10.1016/j.compstruct.2017.11.059. [16] WANG Z G, LI Z D, XIONG W. Experimental investigation on bending behavior of honeycomb sandwich panel with ceramic tile face-sheet [J]. Composites Part B: Engineering, 2019, 164: 280–286. DOI: 10.1016/j.compositesb.2018.10.077. [17] WANG Z G, TIAN H Q, LU Z J, et al. High-speed axial impact of aluminum honeycomb: experiments and simulations [J]. Composites Part B: Engineering, 2014, 56: 1–8. DOI: 10.1016/j.compositesb.2013.07.013. [18] SIBEAUD J M, THAMIÉ L, PUILLET C. Hypervelocity impact on honeycomb target structures: experiments and modeling [J]. International Journal of Impact Engineering, 2008, 35(12): 1799–1807. DOI: 10.1016/j.ijimpeng.2008.07.037. [19] ZHOU X L, ZHOU H F, LI X Y, et al. Size effects on tensile and compressive strengths in metallic glass nanowires [J]. Journal of the Mechanics and Physics of Solids, 2015, 84: 130–144. DOI: 10.1016/j.jmps.2015.07.018. [20] ZHOU Q, MAYER R R. Characterization of aluminum honeycomb material failure in large deformation compression, shear, and tearing [J]. Journal of Engineering Materials and Technology, 2002, 124(4): 412–420. DOI: 10.1115/1.1491575. [21] HOU B, ZHAO H, PATTOFATTO S, et al. Inertia effects on the progressive crushing of aluminium honeycombs under impact loading [J]. International Journal of Solids and Structures, 2012, 49(19/20): 2754–2762. DOI: 10.1016/j.ijsolstr.2012.05.005. [22] WANG Z G, LI Z D, XIONG W. Numerical study on three-point bending behavior of honeycomb sandwich with ceramic tile [J]. Composites Part B: Engineering, 2019, 167: 63–70. DOI: 10.1016/j.compositesb.2018.11.108. [23] WANG D M. Impact behavior and energy absorption of paper honeycomb sandwich panels [J]. International Journal of Impact Engineering, 2009, 36(1): 110–114. DOI: 10.1016/j.ijimpeng.2008.03.002. [24] HEIMBS S, SCHMEER S, MIDDENDORF P, et al. Strain rate effects in phenolic composites and phenolic-impregnated honeycomb structures [J]. Composites Science and Technology, 2007, 67(13): 2827–2837. DOI: 10.1016/j.compscitech.2007.01.027. [25] SHAN J F, XU S L, ZHOU L J, et al. Dynamic fracture of aramid paper honeycomb subjected to impact loading [J]. Composite Structures, 2019, 223: 110962. DOI: 10.1016/j.compstruct.2019.110962. [26] 范成年, 张磊, 单俊芳, 等. 复合纸蜂窝结构的力学性能与吸能特性研究 [J]. 实验力学, 2021, 36(5): 627–637. DOI: 10.7520/1001-4888-21-001.FAN C N, ZHANG L, SHAN J F, et al. The mechanical performance and energy absorption properties of composite paper honeycomb structures [J]. Journal of Experimental Mechanics, 2021, 36(5): 627–637. DOI: 10.7520/1001-4888-21-001. [27] 徐松林, 单俊芳, 周李姜, 等. 对试样进行中等应变速率压缩的实验装置: CN109781555A [P]. 2019-05-21. 北京: 北京集佳知识产权代理有限公司, 2019. [28] RUAN D, LU G, CHEN F L, et al. Compressive behaviour of aluminium foams at low and medium strain rates [J]. Composite Structures, 2002, 57(1): 331–336. DOI: 10.1016/S0263-8223(02)00100-9. [29] XU S Q, RUAN D, BEYNON J H, et al. Experimental investigation of the dynamic behavior of aluminum foams [J]. Materials Science Forum, 2010, 654/655/656: 950–953. DOI: 10.4028/www.scientific.net/MSF.654-656.950. -