不同速度下负泊松比弧形结构的变形行为

罗伟洪; 何婉青; 吴文军; 李世强; 王志勇

doi:10.11883/bzycj/2022-0520

不同速度下负泊松比弧形结构的变形行为

doi: 10.11883/bzycj/2022-0520

太原理工大学机械与运载工程学院应用力学研究所，山西太原 030024

基金项目: 国家自然科学基金（12272257）；山西省自然科学基金面上项目（202203021211169）

详细信息

作者简介:
罗伟洪（1998－　），男，硕士研究生，754314343@qq.com

通讯作者:
王志勇（1982－　），男，博士，副教授，wangzhiyong@tyut.edu.cn

中图分类号: O347
计量
- 文章访问数: 709
- HTML全文浏览量: 184
- PDF下载量: 178
- 被引次数: 7
出版历程
- 收稿日期: 2022-11-16
- 修回日期: 2023-05-12
- 网络出版日期: 2023-05-31
- 刊出日期: 2023-11-17

Deformation behavior of curved structures with negative Poisson’s ratio under diverse loading velocities

Institute of Applied Mechanics, College of Mechanical and Vehicle Engineering, Taiyuan University of Technology, Taiyuan 030024, Shanxi, China

摘要

摘要: 高孔隙率的负泊松比蜂窝结构在能量吸收的过程中往往伴随剧烈的应力波动和显著的峰值应力，极易造成蜂窝结构的局部损坏，影响能量的连续吸收。为了减少局部破坏的产生，基于传统内凹六边形蜂窝胞元设计了一种反对称的弧形胞元，并通过不同的阵列方向，得到了2种新型反对称负泊松比弧形蜂窝结构。采用准静态压缩试验和有限元模拟的方法，探究了速度梯度对新型反对称弧形蜂窝结构模型的整体变形模式，不同层水平应变分布，变形机理，以及抗冲击性的影响。研究结果表明：不同于传统负泊松比蜂窝模型中出现大量的局部密实化区域，新型反对称负泊松比弧形蜂窝结构中局部密实带明显减少，结构中多层胞元组成的变形区域同时参与变形，整体表现出十分稳定的变形模式。这与最大水平应变的提高以及新型蜂窝结构抗冲击性的增强密切相关，特别是在中速模式下，新型反对称弧形蜂窝模型抗冲击性明显增强，冲击载荷效率达到78%，远高于传统蜂窝模型43%的冲击载荷效率；此外，反对称弧形蜂窝结构胞元还带动了相邻胞元之间的胞壁发生向上弯曲来抵抗弯矩。在低速模式下，2种新型反对称弧形蜂窝模型的最大水平应变分别提高了100%、36%；在中速模式下，2种模型均提高了39%。
- 蜂窝结构 /
- 胞元设计 /
- 变形 /
- 负泊松比效应 /
- 速度效应
Abstract: High-porosity structures with negative Poisson’s ratio often experience severe stress fluctuations and significant peak stresses during energy absorption, which can easily cause local damage to the honeycomb structure and affect continuous energy absorption. In order to reduce the occurrence of local damage, an anti-symmetric arc-shaped cell element is designed based on the traditional negative Poisson’s ratio honeycomb cell element, and two new anti-symmetric negative Poisson’s ratio arc-shaped honeycomb structures are obtained through different array directions. Through 0.0025 m/s (quasi-static) compression test and 10 m/s (low velocity), 50 m/s (medium velocity) and 100 m/s (high velocity) finite element simulation, the effect of velocity gradient on the overall deformation pattern, horizontal strain distribution of different layers, deformation mechanism, and impact resistance of the new anti-symmetric arc-shaped honeycomb structure model are revealed. The research results show that unlike the large number of local densification areas that appear in traditional negative Poisson’s ratio honeycomb models, the local densification bands in the new anti-symmetric negative Poisson’s ratio arc-shaped honeycomb structure are significantly reduced. The deformation areas composed of multiple layers of cells in the structure participate in deformation at the same time, showing a very stable deformation pattern as a whole. This is closely related to the increase in maximum horizontal strain and the enhancement of impact resistance of the new honeycomb structure. Especially under the medium-speed loading, the impact resistance of the new anti-symmetric arc-shaped honeycomb model is significantly enhanced, and the impact load efficiency reaches 78%, which is much higher than the 43% impact load efficiency of the traditional honeycomb model; in addition, the anti-symmetric arc-shaped honeycomb structure cells also drive the cell walls between adjacent cells to bend upwards to resist bending moments, further increasing the maximum horizontal strain. Under low-speed loading, the maximum horizontal strain of the two types of new anti-symmetric arc-shaped honeycomb models increases by 100% and 36%, respectively. Under medium-speed loading, it increases by 39% for both types.
- honeycomb structure /
- cell design /
- deformation behavior /
- negative Poisson’s ratio effect /
- loading velocity effect

HTML全文

Zr基非晶合金既具有高强度、高硬度、极大的弹性极限等优异的力学性能，又具有高能量密度和快速释能特性，是一种能综合利用化学能和动能以提高战斗部毁伤效应的新型功能材料。

目前，针对高速冲击下Zr基非晶合金破片反应机制的研究尚处于探索阶段，是前沿研究领域之一。Luo等^[1]利用弹道枪和准密封箱实验，研究了W/Zr基非晶合金复合材料破片的冲击释能特性，并推导了材料的冲击诱发化学反应阈值。Wang等^[2]研究了Zr₅₅Cu₃₀Ni₅Al₁₀非晶合金破片的冲击诱发化学反应过程，通过能谱仪（energy disperse spectroscopy，EDS）分析了Zr₅₅Cu₃₀Ni₅Al₁₀非晶合金破片冲击释能反应后形成碎片的微区元素，发现材料中存在大量氧元素，推测发生了Zr与氧气的燃烧反应。Wei等^[3]利用EDS、三轴晶高分辨X射线衍射（TAXRD）等手段，分析了ZrHf非晶合金破片冲击靶板的反应过程及反应产物，认为ZrO₂为燃烧反应的主产物，但随着碎片厚度增加，氧气向碎片内部扩散能力减弱，碎片内部也存在ZrO、ZrO₃等反应产物。

随着研究的进一步深入，学者们认识到了冲击破碎对Zr基非晶合金破片反应特性有重要影响，Wei等^[4]提出了基于高速摄影图像的像素偏移算法，以统计破片碎片的尺寸分布。Ji等^[5]认为Zr基非晶合金破片冲击破碎后，其碎片分布规律符合幂次律，并收集了冲击实验中的碎片进行验证分析。结果表明，幂次律可以较好地解释Zr基非晶合金碎片分布规律；20 μm是材料碎片是否发生燃烧反应的尺寸阈值，碎片尺寸小于20 μm时，该碎片参与燃烧反应，而当碎片尺寸大于20 μm时，该碎片不参与燃烧反应。上述研究成果初步阐释了Zr基非晶合金破片的冲击破碎反应过程，但对冲击破碎反应机制的研究不够深入，所提出的理论模型未考虑破片冲击升温、碎片尺寸分布等对碎片反应阈值的影响，计算结果与实验结果相差较远。

本文中，通过Zr基非晶合金破片准密封箱冲击超压测试，测量材料的冲击释能规律，分析破片碎片的粒径分布、成分组成；利用冲击升温和碎片燃烧理论分析、验证材料的冲击破碎反应机制。

1. 实　验

1.1 实验方案

1.1.1 准密封箱冲击超压测试

Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片的准密封箱冲击超压测试布置如图1所示， $\varnothing$ 8 mm×10 mm的圆柱形破片装入尼龙弹托中，弹托刻有凹槽，保证出枪口后弹托、破片分离；破片及尼龙弹托由14.5 mm弹道枪发射，通过调节装药量，可以改变破片初速；2块间距为1 m的铝箔断通靶纸与六通道计时仪相连，用以测量破片飞行速度，忽略破片飞行过程中的速度降，测得的平均飞行速度记为破片撞击速度；靶纸置于距枪口3 m处，避免火药气体对测量结果的影响。

图 1 准密封箱冲击超压实验布置

Figure 1. Experimental layout of quasi-sealed chamber shock overpressure test

下载: 全尺寸图片幻灯片

圆柱体准密封箱容积为35.2 L，箱体前端有圆形开口并可由靶板封闭，箱体底端固定有30 mm厚硬质钢板，破片贯穿箱体靶板后撞击硬质钢板并激发化学反应，该装置可模拟预制破片与毁伤目标间的作用过程。破片的冲击破碎反应导致箱体内空气温度、压力升高，箱体内的超压信号可通过箱壁上不同位置上的2个瞬时超压传感器（分别距箱体前端20和40 cm，传感器轴线相互垂直）测量，其触发阈值设置为0.02 MPa。由于超压信号测试过程中箱体靶板已被破片贯穿，故称为准密封箱冲击超压实验。密封箱侧面有圆形窗口，以防弹玻璃封闭，可由高速摄影机观察破片的冲击破碎反应过程，高速摄影频率设置为10 000 s⁻¹。

实际上，破片在撞击并贯穿靶板过程中，有部分材料破碎且无法进入箱体，这部分质量损失对破片的超压测试有一定影响，为减小这部分影响，在实验中测试了破片在500～1600 m/s速度范围内冲击0.5 mm厚Q235钢靶板的靶后超压信号。

1.1.2 粒度测试

实验过程中收集箱体内的破片碎片，并利用2种方法进行粒度测试，分析碎片尺寸分布情况：

称重法，利用1 000、640、320、160、80、40、20和10 μm孔径的样品筛依次筛取碎片试样，通过精密电子天平称取不同粒度下碎片质量；

激光衍射法，利用Beckman Coulter LS13320激光粒度分析仪测量收集试样的粒度分布，仪器的测量范围为0.02～2 000 μm，但当试样粒度分布区间较大时，仪器测量结果精度较低，因此仅用于分析称重法筛取的粒径0.2～640 μm范围碎片的碎片尺寸分布。

实际上，实验中无法完全收集全部碎片，一般地，假设未收集到的碎片与收集碎片粒径分布相同，以收集到的破片碎片试样的粒径分布代表破片全部碎片的粒径分布情况。

1.1.3 X射线衍射（XRD）分析

利用Empyrean X射线衍射仪对碎片进行XRD分析，确定反应产物中存在的物相。

1.2 实验结果

1.2.1 冲击诱发化学反应过程

图2为破片飞行过程中的高速摄影图像，破片用圆圈标出，破片后方为弹托碎片，由于弹托凹槽的设计，破片及弹托出炮口后可以有效分离，弹托与破片分离后弹道性能变差，阻力增大，导致其在破片后方飞行，不会随破片进入密封箱，保证了实验结果的可靠性。

图 2 破片飞行过程高速摄影图像

Figure 2. High-speed photograph of fragment flight

下载: 全尺寸图片幻灯片

图3为高速摄影拍摄的破片分别以552、1105和1485 m/s速度撞击箱体的图像，令拍摄到的破片撞击靶板瞬间时刻为t=0 ms。破片撞击靶板时可以看到明显的火光，破片撞击并贯穿靶板的过程中小部分破碎并激发化学反应，导致了靶板外侧火光的产生（图3(a) t=0 ms）。未破碎的破片主体及部分破片碎片继续飞行，并撞击箱体底部的硬质钢板（图3(a)中t=0.5 ms）。猛烈冲击造成破片进一步破碎并形成碎片云，高温碎片云与空气中的氧气发生剧烈化学反应（图3(a)中t=1.0 ms）。破片化学反应导致箱体内空气温度、压力升高，由于箱体系统内外压力不平衡，箱体内碎片伴随高压气体由靶板弹孔处喷出，形成火舌，直至箱体内外压力相同（图3 (b)中t=10.0 ms）。

图 3 破片冲击破碎反应过程高速摄影图像

Figure 3. High-speed photographs of process for reaction of fragments

下载: 全尺寸图片幻灯片

撞击速度对破片的冲击破碎反应影响很大，由图3 可以看出，随着破片撞击速度的提升，破片反应的火光范围变大，箱体内碎片向外喷射的距离增加，碎片数目增多，如图3(a)中t=10.0 ms基本观察不到有碎片喷射的现象，而图3(c)中t=10.0 ms燃烧碎片的喷射现象明显且剧烈。随着破片撞击速度的提升，破片撞击箱底硬质钢板后破碎程度增加，碎片数量增多，碎片尺寸减小，其与空气接触的比表面积增加，破片的反应程度增大，导致火光范围增大，箱体内瞬时超压更高，箱体内碎片喷射距离更远、数目更多。

1.2.2 冲击超压测试结果

图4为破片撞击速度1105 m/s时，箱体内瞬时超压传感器的实测信号及准静态压力信号。图4(a)和(b)实线分别为2个不同位置传感器的实测信号，由于传感器位置不同，传感器受到的初始冲击等干扰程度不同，2条曲线的噪音幅值差异较大；实测曲线经降噪处理后得到的准静态压力信号在图4(a)和(b)由虚线对应标出，并同时列于图4(c)对比，两条曲线差异极小，表明对于该实验装置，传感器位置对测量结果的影响较小，符合实验要求。

图 4 实测压力曲线与准静态压力曲线

Figure 4. Actual pressure and quasi-static pressure curves

下载: 全尺寸图片幻灯片

图5为实验测得的准静态压力-时间历程，可以看到，曲线基本分为2部分，急速上升段和缓慢卸载段，破片冲击靶板后，激发剧烈的氧化反应，箱体内部压力快速升高，导致曲线上升段斜率较大；当化学反应导致的压力上升速率小于高压气体通过弹孔排出导致的压力下降速率时，曲线开始向下倾斜；箱体内高压气体仅能通过破片贯穿靶板产生的弹孔排出，故而曲线下降段较为平缓。

图 5 准静态压力-时间曲线

Figure 5. Quasi-static pressure-time curves

下载: 全尺寸图片幻灯片

表1为破片在不同速度下的超压峰值数据，m表示破片质量，超压峰值 ${{\Delta }}p$ 表示材料反应过程中超压时程曲线的最大值。

表 1 超压峰值实验数据

Table 1. Experimental data of overpressure peak data

试验编号	v/(m·s⁻¹)	m/g	${{\Delta }}p$ /MPa
1	552	3.43	0.019
2	723	3.4	0.033
3	834	3.42	0.042
4	1105	3.49	0.058
5	1152	3.42	0.069
6	1242	3.36	0.105
7	1406	3.41	0.124
8	1485	3.47	0.151

下载: 导出CSV

| 显示表格

1.2.3 碎片粒度分布规律

Zr基非晶合金受冲击破碎后，碎片着火温度、反应程度均与碎片尺寸相关，建立精确的碎片粒度分布模型对分析破片释能效应具有重要意义。微米级碎片对燃烧反应贡献度较大，为保证测试精度，利用Beckman Coulter LS13320激光粒度分析仪测试0.2～640 μm范围内碎片粒度分布的统计规律。图6为不同撞击速度条件下，Zr基非晶合金破片碎片粒度分析结果。随着撞击速度提升，小粒径碎片体积分数增大；累积体积分数曲线前半部分斜率明显大于曲线后半部分，其拐点位于20 μm。实验仅测试了收集碎片试样质量占破片总质量分数70%以上的样本，以使测试结果更接近实际结果，减小系统误差。

图 6 碎片粒度分布测试结果

Figure 6. Test results of debris size distribution

下载: 全尺寸图片幻灯片

图6(b)数据点为利用称重法得到的碎片样本粒度分布结果，受实验条件限制，称重法的粒度梯度较大，该结果用以检验、校核激光粒度分析仪测试结果。对比2种方法测试结果表明，激光粒度分析仪测试结果与称重法测试结果误差较小，激光粒度分析仪测试可作为进一步理论分析的实验依据。

大量实验数据的统计结果表明，脆性材料动态破碎后，碎片分布符合幂次律，常用的有Rosin-Rammler分布，其质量分布规律为^[6-7]：

$M\left( {{d_{\text{p}}}} \right) = m{\left( {{{{d_{\text{p}}}} \mathord{\left/ {\vphantom {{{d_{\text{p}}}} s}} \right. } s}} \right)^{{k_{\text{s}}}}}$

(1)

式中：M为粒径小于d_p的碎片的累积质量，m为破片质量，k_s为与材料性质相关的形状参数， $s = {\left[ {{{\sqrt {24} {K_{{\text{{\text{Ⅰ}}C}}}}}/ {(\rho c\dot \varepsilon )}}} \right]^{{2/ 3}}}$ 为碎片的特征尺寸，K_ⅠC为材料的断裂韧性，c为材料声速， $\dot \varepsilon$ 为破片的应变率。

通过实验数据拟合，即可得到Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金碎片粒度的Rosin-Rammler分布：

$\lg \left( {\frac{M}{m}} \right) = 0.75\lg \left( {\frac{{{d_{\text{p}}}}}{s}} \right)$

(2)

小尺寸碎片对Zr基非晶合金的燃烧反应贡献较大^[5]，由图6(b)，碎片累积体积分数曲线在20 μm左右有一明显拐点，式(1)无法体现出这种变化。因此，采用多项式形式幂次律描述碎片的尺寸分布：

$\lg \left( {\frac{M}{m}} \right) = {a_{\text{s}}}\lg \left( {\frac{{{d_{\text{p}}}}}{s}} \right) + {b_{\text{s}}}$

(3)

式中：a_s、b_s为待定参数。依据粒度测试结果拟合得到的分段多项式为：

$\lg \left( {\dfrac{M}{m}} \right) = \left\{ \begin{array}{ll} 4.48\lg \left( {\dfrac{{{d_{\text{p}}}}}{s}} \right) + 12.16 & \lg \left( {\dfrac{{{d_{\text{p}}}}}{s}} \right) {\text{≤}} - 3.18 \\ 0.72\lg \left( {\dfrac{{{d_{\text{p}}}}}{s}} \right) + 0.22 & \lg \left( {\dfrac{{{d_{\text{p}}}}}{s}} \right) {\text{＞}} - 3.18 \end{array} \right.$

(4)

图7为得到的碎片粒度分布模型，可以看出，lg(d_p/s)=−3.18约为实验数据拐点，lg(d_p/s)>−3.18 （d_p >20 μm）时，Rosin-Rammler分布可较好地描述Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片碎片的分布特征，但当lg(d_p/s)<−3.18时，实验数据lg(M/m)-lg(d_p/s)曲线斜率增大，Rosin-Rammler分布与实验数据误差较大；相较于Rosin-Rammler分布，式(4)所示分段模型更好地描述了小尺寸碎片粒度分布规律。

图 7 碎片粒度分布模型

Figure 7. Models of debris size distribution

下载: 全尺寸图片幻灯片

1.2.4 XRD分析结果

图8为不同粒径尺度碎片的XRD分析结果，在反应产物分析图谱中可以清晰地观察到ZrO₂、单质Al的特征峰，其他成分特征峰不明显，表明：Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片的冲击破碎反应以Zr的燃烧反应为主，而Al参与燃烧反应较少，推测与碎片温度未达到Al着火点有关。在不同粒径尺度碎片试样中均观察到ZrO₂的特征峰，粒径尺度越小的碎片其特征峰越明显，表明粒径尺度小的碎片ZrO₂结晶程度更高。

图 8 碎片XRD分析结果

Figure 8. Results of XRD analysis of debris

下载: 全尺寸图片幻灯片

2. 理论分析

2.1 能量释放特性

Zr基非晶合金冲击释能主要源自于Zr碎片的燃烧，其反应式为：

${\text{Zr}} + {{\text{O}}_2} \to {\text{Zr}}{{\text{O}}_2},\;\;\;\;{{ \Delta }}{H_{{\text{Zr}}}} = 1.2 \times {10^4}{\text{ kJ}}/{\text{kg}}$

(5)

式中： ${{\Delta }}{H_{{\text{Zr}}}}$ 为反应焓。由于化学反应开始到准密封箱内压力达到超压峰值时间极短，因此，在分析时该段时间内气体泄漏可忽略不计^{[5, 8]}；经计算分析，材料燃烧反应消耗氧气体积占箱体容积极小，也可忽略不计，则箱体内能增量与超压峰值的关系为^{[2, 5]}：

${{\Delta }}p = \frac{{{\gamma _{\text{s}}} - 1}}{{{V_{\text{C}}}}}{{\Delta }}{E_{\text{C}}}$

(6)

式中：V_C为箱体容积， ${\gamma _{\text{s}}}$ =1.4为空气的比热系数，箱体内能增量 $\Delta {E_{\text{C}}}$ 由氧化反应释放化学能 $\Delta Q$ 和部分破片动能E_k组成：

${{\Delta }}{E_{\text{C}}} = {{\Delta }}Q + {E_{\text{k}}}$

(7)

由文献[2]的实验结论，E_k约为破片动能的30%。破片完全反应释放能量为Q_t，则破片的反应程度可由下式计算：

$y = \frac{{{{\Delta }}Q}}{{{Q_{\text{t}}}}}$

(8)

将实验结果代入式(6)～(8)即可计算得到对应撞击速度下的能量及反应程度等，计算结果如表2所示。

表 2 实验数据及计算结果

Table 2. Results of calculation and experimental data

实验	v/(m·s⁻¹)	${{\Delta }}{E_{\text{C}}}$ /kJ	E_k/kJ	${{\Delta }}Q$ /kJ	Q_c/kJ	y/%
1	552	1.67	0.16	1.52	29.32	5.17
2	723	2.90	0.27	2.64	29.06	9.08
3	834	3.70	0.36	3.34	29.23	11.42
4	1 105	5.10	0.64	4.47	29.83	14.97
5	1 152	6.07	0.68	5.39	29.23	18.44
6	1 242	9.24	0.78	8.46	28.72	29.47
7	1 406	10.91	1.01	9.90	29.15	33.97
8	1 485	13.29	1.15	12.24	29.66	40.93

下载: 导出CSV

| 显示表格

可以看到，在500 ～1500 m/s，撞击速度提升导致了材料冲击压力增大，进而诱发了不同程度的化学反应，材料的反应程度与撞击速度、冲击压力正相关，最大冲击释能反应程度达到40.93%。通常，含氧化剂的MESMs如Al/Ni系材料反应程度较高，而需与空气混合反应的负氧平衡类MESMs反应程度较低，Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金材料属于负氧平衡类MESMs，因此在中低速撞击条件下，其反应程度小于50 %。

如图9所示，Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金的反应程度表现出两种特征的速度相关性：低速时，反应程度随速度提升而缓慢增加，高速时，反应程度随速度提升而急剧增加。这也表明，Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片适用于交汇速度高的工况，例如作为防空反导破片使用时，其毁伤效应将更为显著。

图 9 反应程度-撞击速度曲线

Figure 9. Curve of extent of reaction vs. impacting velocity

下载: 全尺寸图片幻灯片

2.2 冲击破碎反应机制

综合破片冲击破碎反应实验结果及反应过程分析，建立冲击破碎反应理论模型以实现Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片反应程度的快速计算。

2.2.1 冲击升温

由破片-靶板撞击的一维正冲击波假设，破片的冲击压力为^[9]：

${p_{\text{H}}} = {\rho _{\text{p}}}{D_{\text{p}}}{u_{\text{p}}}$

(9)

式中：ρ_p、D_p、u_p分别为破片密度（本文中ρ_p=6.76 g/cm³）、冲击波速、材料粒子速度。一般情况下，可假设固态金属材料的冲击波速-粒子速度为线性关系，由破片-靶板界面连续性^[10]：

${u_{\text{p}}} = v - \frac{{ - \left( {{\rho _{\text{t}}}{c_{\text{t}}} + {\rho _{\text{p}}}{c_{\text{p}}} + 2{\rho _{\text{p}}}{\lambda _{\text{p}}}v} \right) \pm {\varDelta ^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}}}{{2\left( {{\rho _{\text{t}}}{\lambda _{\text{t}}} - {\rho _{\text{p}}}{\lambda _{\text{p}}}} \right)}}$

(10)

式中：v为破片撞击速度，ρ_t为靶板密度（本文中ρ_t=7.85 g/cm³），c_p=3 868 m/s、c_t=4 570 m/s、λ_p=3.97、λ_t=1.49^[11]分别为破片、靶板材料的Hugoniot参数， $\varDelta$ 的表达式为^[10]：

$\varDelta = {\left( {{\rho _{\text{t}}}{c_{\text{t}}} + {\rho _{\text{p}}}{c_{\text{p}}} + 2{\rho _{\text{p}}}{\lambda _{\text{p}}}v} \right)^2} + 4{\rho _{\text{p}}}\left( {{\rho _{\text{t}}}{\lambda _{\text{t}}} - {\rho _{\text{p}}}{\lambda _{\text{p}}}} \right)\left( {{c_{\text{p}}}v + {\lambda _{\text{p}}}{v^2}} \right)$

(11)

根据固态材料绝热线和等熵线的关系，假定材料等容比热容c_V=0.32 J·g⁻¹·K^−1[11]为常数，得到冲击温度 ${T_{\text{H}}}$ 的表达式为^[9]：

${T_{\text{H}}} = {T_0}\exp \left( { - \int_{{V_0}}^{{V_{\text{H}}}} {\frac{V}{\gamma }{\text{d}}V} } \right) + \frac{V}{\gamma }\frac{{{p_{\text{H}}} - {p_{\text{S}}}}}{{{c_{{V}}}}}$

(12)

式中：T₀为常态温度，V为比容，下标0、H分别表示初始状态和冲击绝热状态，γ为Grüneisen系数，p_S为与冲击绝热线同一起始状态出发的等熵线，其解析式为^[9]：

${p_{\text{S}}}\left( \eta \right) = {\rho _{{\text{p}}0}}c_{\text{p}}^2\exp \left( {{\gamma _0}\eta } \right)\int_0^\eta {\frac{{\left( {1 + {\lambda _{\text{p}}^{} - {\gamma _0}} } \right)\eta '}}{{{{\left( {1 - \lambda _{\text{p}}^{}\eta '} \right)}^3}}}} \exp \left( { - {\gamma _0}\eta '} \right){\text{d}}\eta '$

(13)

式中： $\eta = {\mu \mathord{\left/ {\vphantom {\mu {\left( {1 + \mu } \right) = }}} \right. } {\left( {1 + \mu } \right) = }}1 - {V \mathord{\left/ {\vphantom {V {{V_0}}}} \right. } {{V_0}}}$ 为压缩度，ρ_p0为破片初始密度，γ₀为初始Grüneisen系数。

将材料参数代入式(9)～(13)即可分别求得破片冲击压力、冲击温度随破片撞击速度的变化情况，如图10所示。

图 10 Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金的冲击压力、温度与撞击速度的关系

Figure 10. Relationships among shock pressure, shock temperature and impact velocity of Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆ bulk metallic glass

下载: 全尺寸图片幻灯片

2.2.2 碎片燃烧

Zr基非晶合金破片受撞击后破碎，形成大量碎片并与空气中的氧气发生燃烧反应。假设碎片云的燃烧反应为大量单一碎片燃烧反应的总和，依据单一碎片燃烧反应规律及破片碎片尺寸分布规律即可求得Zr基非晶合金破片的冲击破碎反应规律。金属碎片自持燃烧的临界温度即为着火温度，其判据的数学表达式为^[12]：

${Q_{\text{e}}} = {Q_{\text{r}}} + {Q_{\text{c}}}$

(14)

式中：Q_e为燃烧反应产热率，Q_r和Q_c分别为金属碎片通过热辐射和热对流效应对环境的放热率。

金属碎片通过热辐射对周围环境的散热率为^[13]：

${Q_{\text{r}}} = {\varepsilon _{\text{r}}}{\sigma _{{\text{SB}}}}\left( {{\text{π}}{d_{\text{p}}^2}} \right)\left( {{T_{\text{p}}^4} - {T_{\text{e}}^4}} \right)$

(15)

式中：d_p为碎片直径， ${\varepsilon _{\text{r}}}$ 为辐射系数， ${\sigma _{{\text{SB}}}}$ 为Stefan-Boltzman常数，T_p为碎片温度，T_e为环境温度。

对微米、纳米级碎片，其与环境的热交换机制由Knudsen数K控制，K的表达式为^[13]：

$K = \frac{{4{k_{\text{g}}}}}{{5{r_{\text{p}}}{p_{\text{g}}}}}\sqrt {\frac{{{m_{\text{g}}}{T_{\text{e}}}}}{{2{k_{\text{B}}}}}}$

(16)

式中：k_g为气体热导率，r_p为碎片半径，p_g为气体压力，m_g为气体平均分子质量，k_B为Boltzman常数。当金属碎片足够小（K>10）时，热对流机制为自由分子机制，热交换率取决于金属碎片、空气分子间碰撞概率^[13-15]，此时热对流效应放热率Q_FM为：

${Q_{{\text{FM}}}} = {\text{π}}{r_{\text{p}}^2}\frac{{{p_{\text{g}}}}}{2} {\left( {\frac{{8{k_{\text{B}}}{T_{\text{e}}}}}{{{\text{π }} {m_{\text{g}}}}}} \right)^{{1 / 2}}}\frac{{\gamma _{{\text{ab}}} + 1}}{{\gamma _{{\text{ab}}} - 1}}\left( {\frac{{{T_{\text{p}}} - {T_{\text{e}}}}}{{{T_{\text{e}}}}}} \right)$

(17)

式中： $\gamma _{{\text{ab}}}^{}$ 为气体平均绝热指数。当金属碎片足够大（K<0.01）时，热对流机制为连续体机制^[13-14]，此时热对流效应放热率Q_D为：

${Q_{\text{D}}} = 4{\text{π}} {r_{\text{p}}}\bar k_{\text{g}}\left( {{T_{\text{p}}} - {T_{\text{e}}}} \right)$

(18)

式中： $\bar k_{\text{g}}^{}$ 为平均热导率。大多数情况下，K介于0.01和10之间，此时热对流机制为过渡机制，其热交换率可定义为^[13]：

$\frac{1}{{{Q_{\text{c}}}}} = \frac{1}{{{Q_{\text{D}}}}} + \frac{1}{{{Q_{{\text{FM}}}}}}$

(19)

通过式(16)计算金属碎片的Knudsen数，根据对应的热交换机制，即可得到金属碎片对环境的放热率。

纳米、微米级金属碎片的燃烧主要由氧通过氧化层的固态扩散机制控制^[16]，假设破片碎片为球形，氧气在氧化层中的固态扩散为等温稳态扩散且摩尔体积不变，燃烧反应产热率可表示为^[17]：

${Q_{\text{e}}} = {{\Delta }}H {m_{\text{p}}} \dfrac{{3{{\left( {{r_0} - \xi } \right)}^2}}}{{{r_{0}^3}}} \dfrac{{{\varOmega _0}\exp \left( { - \dfrac{{{E_{\text{a}}}\left( {{r_{\text{p}}}} \right)}}{{RT}}} \right)}}{{\left( {1 - {\dfrac{\xi }{{{r_0}}}} } \right)\xi }}$

(20)

式中：m_p为碎片质量，r₀为碎片初始半径，ξ为氧化层厚度，Ω₀为常数，R为普适气体常数，E_a为燃烧反应表观活化能，对于微米乃至纳米级别的金属碎片，E_a为与粒径r_p相关的变量。对于Zr碎片，其E_a（单位：kJ/mol）与d_p（单位：μm）的关系为^[16]：

${E_{\text{a}}}\left( {{\text{Zr}}} \right) = 40\,080\, \lg \left( {{d_{\text{p}}}} \right) + 296\,200\left( { \pm 30\,200} \right)$

(21)

由式(14)～(21)，通过数值法即可得到Zr碎片自持燃烧的着火温度与粒径的关系，如图11所示。

图 11 Zr碎片着火温度的尺寸效应

Figure 11. Size effect of ignition temperature for Zr fragments

下载: 全尺寸图片幻灯片

由反应物XRD分析结论，实验中，由于Al的着火温度较高，Al参与反应极少，Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金的化学反应主要来为Zr的燃烧。计算得到着火温度为冲击温度时所对应的Zr碎片粒径d_p1，将d_p1代入材料碎片尺寸分布规律（式(5)），则Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合破片的反应程度可估算为：

$y = \left\{ \begin{array}{ll} \dfrac{{{m_1}m{{\left( {{{{d_{{\text{p1}}}}} / s}} \right)}^{4.48}} \times {{10}^{12.16}} \times {{\Delta }}{H_{{\text{Zr}}}}}}{{m{{\Delta }}{H_{{\text{BMG}}}}}} &{\text{lg}}\left( {{{{d_{{\text{p1}}}}} / s}} \right){\text{≤}} - 3.18 \\ \dfrac{{{m_1}m{{\left( {{{{d_{{\text{p1}}}}}/ s}} \right)}^{0.72}} \times {{10}^{0.22}} \times {{\Delta }}{H_{{\text{Zr}}}}}}{{m{{\Delta }}{H_{{\text{BMG}}}}}} & {\text{lg}}\left( {{{{d_{{\text{p1}}}}} /s}} \right) {\text{＞}} - 3.18 \\ \end{array} \right.$

(22)

式中：m₁为Zr在Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金中的质量分数， ${{\Delta }}{H_{{\text{BMG}}}}$ 为Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金的氧化焓。

图12为Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片的反应程度-撞击速度关系和反应程度-冲击压力关系。可以看出，随着撞击速度、冲击压力上升，材料的反应程度也逐步提高，理论计算结果与实验结果较为吻合，表明冲击破碎燃烧模型可较好地阐释Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片冲击破碎反应的物理过程和反应机制。

图 12 反应深度与冲击速度、冲击压力关系的实验与理论计算结果对比

Figure 12. Comparison between experimental results and theoretical results of relationships among reaction extent, impact velocity and shock pressure

下载: 全尺寸图片幻灯片

3. 结　论

(1) 材料冲击释能反应规律为：随着撞击速度提高，Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片的反应程度逐渐提高，其反应程度表现出2种特征的速度相关性，低速时，反应程度随速度提升而缓慢增加，高速时，反应程度随速度提升而急剧增加。

(2) 材料冲击破碎反应机制为：破片撞击目标产生冲击波，冲击波在破片内的传播使材料处于高温状态；破片-靶板界面材料处于高压状态，局部压力大于破片强度，从断裂力学角度考虑，破片局部累积的弹性应变能大于其断裂表面能，导致破片变形、破碎；达到着火条件的碎片与空气中的氧气发生燃烧反应（主要燃烧产物为ZrO₂），释放能量。破片撞击速度为控制碎片温度、碎片粒径分布、破片反应程度的关键控制变量。

(3) 建立了Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片冲击破碎反应理论模型，该模型在500～1 600 m/s速度范围内能较好地解释冲击作用下Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金破片的反应规律，理论计算与实验结果吻合度较高。

图 1 RH、CRH-1和CRH-2三种蜂窝模型设计思路

Figure 1. RH, CRH-1 and CRH-2 Three types of cellular model design ideas

下载: 全尺寸图片幻灯片

图 2 3D打印模型尺寸

Figure 2. 3D printed model size

下载: 全尺寸图片幻灯片

图 3 准静态下3种模型应力-应变曲线

Figure 3. Quasi-static stress-strain curves for three models

下载: 全尺寸图片幻灯片

图 4 准静态下3种结构的变形模式

Figure 4. Three kinds of model deformation under quasi-static loading

下载: 全尺寸图片幻灯片

图 5 模型加载边界条件

Figure 5. Loading boundary conditions

下载: 全尺寸图片幻灯片

图 6 有限元模型

Figure 6. Finite element models

下载: 全尺寸图片幻灯片

图 7 RH模型的变形模式

Figure 7. Deformation modes of the RH model

下载: 全尺寸图片幻灯片

图 8 名义应力-应变曲线

Figure 8. Nominal stress-strain curves

下载: 全尺寸图片幻灯片

图 9 不同速度下RH模型每一层的水平应变分布

Figure 9. Horizontal strain distribution of each layer of the RH model at different velocities

下载: 全尺寸图片幻灯片

图 10 不同速度下CRH-1模型每一层的水平应变分布

Figure 10. Horizontal strain distribution of each layer of the CRH-1 model at different velocities

下载: 全尺寸图片幻灯片

图 11 不同速度下CRH-2模型每一层的水平应变分布

Figure 11. Horizontal strain distribution of each layer of the CRH-2 model at different velocities

下载: 全尺寸图片幻灯片

图 12 低速模式

Figure 12. Low-velocity mode

下载: 全尺寸图片幻灯片

图 13 高速模式

Figure 13. High-velocity mode

下载: 全尺寸图片幻灯片

图 14 RH、CRH-1和CRH-2模型在不同速度下的应变-应力曲线

Figure 14. Strain-stress curves of RH, CRH-1 and CRH-2 models under different velocities

下载: 全尺寸图片幻灯片

表 1 材料属性

Table 1. Material properties

材料密度/（kg·m⁻³）弹性模量/GPa 泊松比初始屈服强度/MPa 断裂应变

PolyMax^TMPLA 1180 1.97 0.35 40 0.3

下载: 导出CSV

表 2 3种模型在10 m/s 速度下的变形模式

Table 2. Deformation patterns of three types of models under 10 m/s velocity

ε RH CRH-1 CRH-2

0.1

0.3

0.5

0.7

下载: 导出CSV

表 3 3种模型在50 m/s 速度下的变形模式

Table 3. Deformation patterns of three types of models under 50 m/s velocity

ε RH CRH-1 CRH-2

0.1

0.3

0.5

0.7

下载: 导出CSV

表 4 3种模型在100 m/s速度下的变形模式

Table 4. Deformation patterns of three types of models under 100 m/s velocity

ε RH CRH-1 CRH-2

0.1

0.3

0.5

0.7

下载: 导出CSV

参考文献(23)

[1]	YAO Y Z, SHEN Y X, ZHU L Q, et al. Preliminary study and bioinformatics analysis on the potential role of CagQ in type IV secretion system of H. pylori [J]. Microbial Pathogenesis, 2018, 116: 1–7. DOI: 10.1016/j.micpath.2017.12.076.
[2]	KUMAR P, KUCHEROV L, RYVKIN M. Fracture toughness of self-similar hierarchical material [J]. International Journal of Solids and Structures, 2020, 203: 210–223. DOI: 10.1016/j.ijsolstr.2020.07.011.
[3]	WANG Z G. Recent advances in novel metallic honeycomb structure [J]. Composites Part B: Engineering, 2019, 166: 731–741. DOI: 10.1016/j.compositesb.2019.02.011.
[4]	SPADONI A, RUZZENE M, SCARPA F. Dynamic response of chiral truss-core assemblies [J]. Journal of Intelligent Material Systems and Structures, 2006, 17(11): 941–952. DOI: 10.1177/1045389x06060219.
[5]	QI C, JIANG F, YANG S, et al. Dynamic crushing response of novel re-entrant circular auxetic honeycombs: numerical simulation and theoretical analysis [J]. Aerospace Science and Technology, 2022, 124: 107548. DOI: 10.1016/j.ast.2022.107548.
[6]	LI S Q, YU B L, KARAGIOZOVA D, et al. Experimental, numerical, and theoretical studies of the response of short cylindrical stainless steel tubes under lateral air blast loading [J]. International Journal of Impact Engineering, 2019, 124: 48–60. DOI: 10.1016/j.ijimpeng.2018.10.004.
[7]	ALMGREN R F. An isotropic three-dimensional structure with Poisson’s ratio=−1 [J]. Journal of Elasticity, 1985, 15(4): 427–430. DOI: 10.1007/BF00042531.
[8]	WILT J K, YANG C, GU G X. Accelerating auxetic metamaterial design with deep learning [J]. Advanced Engineering Materials, 2020, 22(5): 1901266. DOI: 10.1002/adem.201901266.
[9]	韩会龙, 张新春. 星形节点周期性蜂窝结构的面内动力学响应特性研究 [J]. 振动与冲击, 2017, 36(23): 223–231. DOI: 10.13465/j.cnki.jvs.2017.23.033. HAN H L, ZHANG X C. In-plane dynamic impact response characteristics of periodic 4-point star-shaped honeycomb structures [J]. Journal of Vibration and Shock, 2017, 36(23): 223–231. DOI: 10.13465/j.cnki.jvs.2017.23.033.
[10]	XIAO D B, DONG Z C, LI Y, et al. Compression behavior of the graded metallic auxetic reentrant honeycomb: experiment and finite element analysis [J]. Materials Science and Engineering: A, 2019, 758: 163–171. DOI: 10.1016/j.msea.2019.04.116.
[11]	ZHANG J J, LU G X. Dynamic tensile behaviour of re-entrant honeycombs [J]. International Journal of Impact Engineering, 2020, 139: 103497. DOI: 10.1016/j.ijimpeng.2019.103497.
[12]	HU L L, ZHOU M Z, DENG H. Dynamic crushing response of auxetic honeycombs under large deformation: theoretical analysis and numerical simulation [J]. Thin-Walled Structures, 2018, 131: 373–384. DOI: 10.1016/j.tws.2018.04.020.
[13]	JIN X C, WANG Z H, NING J G, et al. Dynamic response of sandwich structures with graded auxetic honeycomb cores under blast loading [J]. Composites Part B: Engineering, 2016, 106: 206–217. DOI: 10.1016/j.compositesb.2016.09.037.
[14]	FENG J W, FU J Z, YAO X H, et al. Triply periodic minimal surface (TPMS) porous structures: from multi-scale design, precise additive manufacturing to multidisciplinary applications [J]. International Journal of Extreme Manufacturing, 2022, 4(2): 022001. DOI: 10.1088/2631-7990/ac5be6.
[15]	DONG Z C, LI Y, ZHAO T, et al. Experimental and numerical studies on the compressive mechanical properties of the metallic auxetic reentrant honeycomb [J]. Materials & Design, 2019, 182: 108036. DOI: 10.1016/j.matdes.2019.108036.
[16]	韩会龙, 张新春, 王鹏. 负泊松比蜂窝材料的动力学响应及能量吸收特性 [J]. 爆炸与冲击, 2019, 39(1): 013103. DOI: 10.11883/bzycj-2017-0281. HAN H L, ZHANG X C, WANG P. Dynamic responses and energy absorption properties of honeycombs with negative Poisson’s ratio [J]. Explosion and Shock Waves, 2019, 39(1): 013103. DOI: 10.11883/bzycj-2017-0281.
[17]	GUO Y G, ZHANG J, CHEN L M, et al. Deformation behaviors and energy absorption of auxetic lattice cylindrical structures under axial crushing load [J]. Aerospace Science and Technology, 2020, 98: 105662. DOI: 10.1016/j.ast.2019.105662.
[18]	XIAO D B, KANG X, LI Y, et al. Insight into the negative Poisson’s ratio effect of metallic auxetic reentrant honeycomb under dynamic compression [J]. Materials Science and Engineering: A, 2019, 763: 138151. DOI: 10.1016/j.msea.2019.138151.
[19]	LIU W Y, WANG N L, LUO T, et al. In-plane dynamic crushing of re-entrant auxetic cellular structure [J]. Materials & Design, 2016, 100: 84–91. DOI: 10.1016/j.matdes.2016.03.086.
[20]	KOOISTRA G W, DESHPANDE V S, WADLEY H N G. Compressive behavior of age hardenable tetrahedral lattice truss structures made from aluminium [J]. Acta Materialia, 2004, 52(14): 4229–4237. DOI: 10.1016/j.actamat.2004.05.039.
[21]	GIBSON L J, ASHBY M F. Cellular solids: structure and properties [M]. 2nd ed. Cambridge: Cambridge University Press, 1997: 1-13. DOI: 10.1017/CBO9781139878326.
[22]	QIU X M, ZHANG J, YU T X. Collapse of periodic planar lattices under uniaxial compression, part Ⅱ: dynamic crushing based on finite element simulation [J]. International Journal of Impact Engineering, 2009, 36(10/11): 1231–1241. DOI: 10.1016/j.ijimpeng.2009.05.010.
[23]	SUN D Q, ZHANG W H, ZHAO Y C, et al. In-plane crushing and energy absorption performance of multi-layer regularly arranged circular honeycombs [J]. Composite Structures, 2013, 96: 726–735. DOI: 10.1016/j.compstruct.2012.10.008.

施引文献

期刊类型引用(3)

1.	毛光辉，王成，王万里，徐文龙. 变形可控内凹蜂窝结构抗冲击性能研究. 兵工学报. 2025(03): 116-130 . 百度学术
2.	石南南，张伟晨，李振宝，王利辉，刘晗，张雷，夏阳. 泡沫填充四韧带反手性结构和内凹结构的面内压缩性能. 复合材料学报. 2024(06): 2935-2946 . 百度学术
3.	李娜，刘述尊，张新春，张英杰，齐文睿. 新型星-菱形负泊松比蜂窝结构的动态力学特性. 复合材料学报. 2024(09): 4956-4967 . 百度学术

其他类型引用(4)

资源附件(0)

访问统计

图(14) / 表(4)

计量

文章访问数: 709
HTML全文浏览量: 184
PDF下载量: 178
被引次数: 7

1. 实　验
1.1 实验方案
1.2 实验结果
2. 理论分析
2.1 能量释放特性
2.2 冲击破碎反应机制
3. 结　论

材料	密度/（kg·m⁻³）	弹性模量/GPa	泊松比	初始屈服强度/MPa	断裂应变
PolyMax^TMPLA	1180	1.97	0.35	40	0.3

ε	RH	CRH-1	CRH-2
0.1
0.3
0.5
0.7

ε	RH	CRH-1	CRH-2
0.1
0.3
0.5
0.7

ε	RH	CRH-1	CRH-2
0.1
0.3
0.5
0.7

不同速度下负泊松比弧形结构的变形行为

doi: 10.11883/bzycj/2022-0520

作者简介: 罗伟洪（1998－ ），男，硕士研究生，754314343@qq.com

通讯作者: 王志勇（1982－ ），男，博士，副教授，wangzhiyong@tyut.edu.cn

计量

出版历程

Deformation behavior of curved structures with negative Poisson’s ratio under diverse loading velocities

1. 实 验

1.1 实验方案

1.1.1 准密封箱冲击超压测试

1.1.2 粒度测试

1.1.3 X射线衍射（XRD）分析

1.2 实验结果

1.2.1 冲击诱发化学反应过程

1.2.2 冲击超压测试结果

1.2.3 碎片粒度分布规律

1.2.4 XRD分析结果

2. 理论分析

2.1 能量释放特性

2.2 冲击破碎反应机制

2.2.1 冲击升温

2.2.2 碎片燃烧

3. 结 论

期刊类型引用(3)

其他类型引用(4)

计量

出版历程

目录

1. 实 验

1.1 实验方案

1.2 实验结果

2. 理论分析

2.1 能量释放特性

2.2 冲击破碎反应机制

3. 结 论

作者简介:
罗伟洪（1998－　），男，硕士研究生，754314343@qq.com

通讯作者:
王志勇（1982－　），男，博士，副教授，wangzhiyong@tyut.edu.cn

1. 实　验

3. 结　论

1. 实　验

3. 结　论