Study on Anisotropic Crushing Behavior of the Functionally Gradient Aluminum Foam
-
摘要: 功能梯度泡沫金属因其密度连续变化，在轴向受压时可提供稳定增长的反馈载荷。然而当前研究多局限于其纵向压缩力学响应，考虑到实际应用中可能出现的横向冲击，基于低速冲击实验，考察梯度泡沫铝轴向和横向压缩力学响应的异同，并采用数字图像相关技术和数值模拟方法研究其宏细观压溃机制。结果表明：（1）在力学性能上，相比于纵向压缩加载的梯度泡沫铝，横向压缩加载下具有更高的抗压强度，而平台应力、致密化应变和能量吸收效果低于纵向压缩；（2）在失效变形模式上，纵向压缩变形模式为变形带渐进式压缩，而横向压缩变形模式的变形带则随机出现在试样的各个位置；（3）横向压缩下梯度泡沫铝致密化应变和比吸能的减小是由高孔隙率区的胞孔利用率降低导致的；（4）构建的弹性-塑性硬化-刚性模型能够较准确地描述梯度泡沫铝的纵向压缩力学行为。研究结果可为梯度泡沫金属在爆炸冲击结构防护工程中的设计提供理论参考。Abstract: Due to the continuous variation of material density in functionally gradient aluminum foams, a feedback load gradually increases within the compression process. Currently, most of the researches have been focused on the longitudinal compression mechanical response. However, considering the possible transverse impact in practical application, we investigated functionally gradient aluminum foams in terms of the axial and transverse compressive mechanical responses based on low-speed impact experiments, and also studied their macro and mesoscopic crushing mechanism via the digital image correlation and numerical simulation technology. The results were shown that: (1) Compared with that under longitudinal compression, the functional gradient aluminum under transverse compression has higher compressive strength and lower platform stress, densification strain and energy absorption. (2) In the aspect of failure deformation mode, the longitudinal compression deformation mode is progressive compression of deformation band, while the deformation band of transverse compression appears randomly at each position of the sample. (3) Under transverse compression, the densification strain and specific energy absorption of the functionally gradient aluminum foam get decreased by the reduction of cellular utilization in the higher porosity zone. (4) The established elastic-plastic hardening-locking (E-PH-L) constitutive model can well capture the longitude compression response of the functionally gradient aluminum foam. This study therefore can provide a reference for the theoretical design of functionally gradient aluminum foam in the protection engineering of explosive impact structures.
-
金属柱壳在内部爆炸加载条件下的动态变形、破碎特性有着重要的理论意义和应用价值。在军事领域，由于柱壳结构在战斗部中较为常见，其在爆炸载荷下的力学响应及终点效应一直是关注的热点[1-2]；在公共安全领域，为有效限制爆炸冲击波及爆轰产物的作用范围，对人员和设备实现有效保护，抗爆容器可用于紧急处理疑似爆炸物品[3-4]；在石油化工领域，压力容器经常在高温和高压下运行，在突发故障下内部易燃易爆介质可能发生爆炸，产生碎片引发多米诺事故也是不可避免的[5-6]。
国内外早期的研究主要集中在柱壳膨胀后期裂纹形成的物理机制、断裂应变分析和破片统计分布规律。Gurnery[7]根据能量平衡模型估算了柱壳外爆膨胀断裂后的破片速度。Taylor[8]基于拉伸断裂假设及弹塑性理论，提出柱壳内部的应力状态与断裂判据。Hoggatt等[9]观察到，随着爆轰压力的增加，柱壳的断裂模式由拉伸型向剪切型转变。Mott[10]和Grady[11]分别从统计学方法和能量守恒观点出发，系统地研究了柱壳破碎后形成大量破片的统计分布规律，并给出相应的破片分布模型。近些年随着计算机技术的进步，数值仿真成为广泛使用的一种研究方法。Hiroe等[12]基于实验与数值模拟结合的方法，研究了不同材料、缺陷结构、爆炸能量及起爆位置对柱壳变形及破碎行为的影响规律。Tanapornraweekit等[13]通过数值仿真研究了柱壳自然破碎后其破片初速度、飞散角和质量分布特征，并分析了不同材料性能对柱壳破碎性能的影响。Kong等[14]利用光滑粒子数值法计算了柱壳在内爆载荷下的动态断裂行为，破片分布及飞散速度与实验结果基本一致。由以上可见，系统研究不同条件下壳体破碎特性的报道并不多见[15]。针对这类强动载柱壳结构，国内外也未形成一种通用的破碎性评估准则。
本工作以圆柱钢壳为研究对象，分析不同装药条件对壳体膨胀断裂行为的影响，探讨破碎特性和装药与壳体质量比C/M(C为装药质量，M为壳体质量)的内在关系，得出一种新的适用于圆柱壳体的破碎性评估准则。
1. 实验
1.1 实验方法
试验材料为某典型高强合金钢，其主要化学成分包括0.6%C、1.27%Si、1.15%Mn、0.018%S、0.027%P。装药为B炸药，密度为1.67 g/cm3。所有圆柱钢壳在水井中爆破，且壳体周围预留足够空气区域供其膨胀断裂，按标准方法将破片回收，回收率均在97%以上。按质量分组并记录各组破片的质量和数量，其中质量小于0.04 g的破片只统计总质量。
1.2 破碎性表征
一般认为，壳体破片的质量分布符合修正Payman直线方程[16]
lgP=2−C0mM (1) 式中:m为规定的破片分组质量范围的下限; M为回收破片的总质量; P为质量大于m的破片质量百分数(0＜P≤100);直线的斜率C0称为修正Payman破碎参数，可用来表征壳体的破碎性能。一般来说，C0值越大，壳体破片的质量分布越合理，破碎性越好。
研究还发现，修正Payman破碎参数C0与壳体长度L成正比，即
Cu=C0/L (2) 式中：Cu即为单位长度修正Payman破碎参数，单位cm-1。Cu与壳体长度无关，只与壳体材料、装药类型及C/M比值有关。
2. 结果分析与讨论
2.1 壳体破碎性能评估准则
为研究壳体的破碎特性与C/M比值的关系，设计了如图 1所示的圆柱壳体，具体尺寸见表 1。为消除端部效应，药柱两端均突出壳体约20 mm。共进行5组试验，每组3发壳体。
表 1 壳体尺寸及装药条件Table 1. Dimensions of shell and chargeNo. D/mm d/mm L/mm C/M 1 30.64 20.64 70 0.18 2 41.99 30.13 70 0.23 3 40.13 30.13 70 0.28 4 50.87 39.52 70 0.33 5 49.52 39.52 70 0.38 破片回收后按质量分拣，并根据(1)式分析其质量分布规律，可得出修正Payman破碎参数C0与C/M值的关系。如图 2所示，随着C/M比值的增加，C0变大，壳体的破碎性能提高，该现象与壳体断裂及破片形成机制有关，与Reid等[17]的实验结果相符。
图 2 破碎参数C0与C/M值的关系Figure 2. Variation of Payman fragmentation parameter C0 with C/M ratio对每组C0取平均值并除以壳体长度，可得到单位长度修正Payman破碎参数Cu。如图 3所示，对Cu与C/M比值进行线性拟合，拟合直线方程可表示为
Cu=19.1+92.8C/M (3) 
图 3 破碎参数Cu与C/M比值的关系Figure 3. Variation of normalized Payman fragmentation parameter Cu with C/M ratio因此可认为，内爆炸加载条件下，任何圆柱钢壳与炸药组合一般都存在这类线性关系，即
Cu=a+bC/M (4) 式中：a和b被称为壳体材料与炸药组合的特征常数。对于给定的材料与炸药组合，a和b为定值。根据壳体破碎的相似性原则，利用(4)式可以评估在相同壳体材料及炸药条件下任何尺寸圆柱壳体的破碎参数，具有广泛的用途。由图 3还可以看出，当C/M比值为0.33和0.38时，破碎参数Cu值偏低。这是因为这两组壳体的直径较大，毛坯热锻后冷却特别慢，温度不均匀导致组织中铁素体含量增多，并呈块状形态，阻碍了裂纹扩展与贯通，最终造成壳体整体的破碎性降低。
2.2 端部效应对破碎性的影响
为研究端部效应对破碎性的影响，设计了如图 4所示的圆柱钢壳，其特点是起爆端药柱与壳体齐平。所有壳体分为两组：一组壁厚t固定，C/M比值变化，考察C/M比值对壳体破碎性的影响；另一组C/M比值保持不变，壁厚变化，考察壁厚对壳体破碎性的影响。壳体尺寸如表 2所示，得到的破碎性结果如图 5所示。可以看出, 当存在端部效应时，壳体破碎性随C/M比值的增大而提高，单位长度修正Payman破碎参数Cu与C/M比值同样可拟合成直线关系，其方程可表示为
Cu=4.7+93.1C/M (5) 表 2 考虑端部效应的壳体尺寸及装药条件Table 2. Dimensions of shell with end effect and chargeNo. D/mm d/mm L/mm t/mm C/M 1 30.64 20.64 60 5 0.18 2 40.13 30.13 60 5 0.28 3 49.52 39.52 60 5 0.38 4 58.84 48.84 60 5 0.48 5 32.11 24.11 60 4 0.28 6 48.16 36.16 60 6 0.28 7 56.19 42.19 60 7 0.28 
图 5 端部效应对壳体破碎参数Cu的影响Figure 5. Influence of end effect on the normalized Payman fragmentation parameter Cu图 5中虚线为(3)式拟合直线，表示壳体不存在端部效应时，其破碎参数Cu与C/M比值的变化关系。由此可以得出，无论C/M比值如何变化，端部效应将导致壳体破碎性能降低。这是因为该条件下装药起爆不能立刻达到最高压力，导致壳体起爆端不能及时断裂形成破片，整体破碎性能变小。需要注意的是，图 5中两条直线的斜率分别为92.8和93.1，近似相等，意味着不管C/M比值多大，由端部效应引起破碎参数Cu的降低量是个常数。
如果忽略两条直线斜率的微小差别，将(3)式与(5)式相减，即可得到如下关系式
Cu_end=Cu−34a (6) 式中：Cu_end和Cu分别为相同C/M比值条件下存在端部效应和无端部效应的破碎参数，a为材料与炸药组合的特性常数。
当C/M比值保持不变时，破碎参数Cu随壳体壁厚的变化趋势如图 6所示。可以看出，Cu与壁厚的关系基本为一条水平直线，具有很小的正斜率，表明端部效应与壁厚无关。
图 6 壁厚对壳体破碎参数Cu的影响Figure 6. Influence of wall thickness on the normalized Payman fragmentation parameter Cu2.3 无装药部分对破碎性的影响
为考察无装药部分对壳体破碎性的影响，设计了如图 7所示的圆柱钢壳。圆柱钢壳共分为3组：(1)壁厚t及无装药部分长度L′固定，C/M比值变化；(2)壁厚及C/M比值固定，无装药部分长度变化；(3) C/M比值及无装药部分长度固定，壁厚变化。壳体尺寸及装药条件见表 3。
表 3 考虑无装药部分的壳体尺寸及装药条件Table 3. Dimensions of shell with charge vacancy and chargeNo. D/mm d/mm L/mm L′/mm t/mm C/M 1 30.64 20.64 60 30 5 0.18 40.13 30.13 60 30 5 0.28 49.52 39.52 60 30 5 0.38 58.84 48.84 60 30 5 0.48 2 40.13 30.13 60 10 5 0.28 40.13 30.13 60 20 5 0.28 40.13 30.13 60 30 5 0.28 40.13 30.13 60 40 5 0.28 3 32.11 24.11 60 30 4 0.28 40.13 30.13 60 30 5 0.28 48.16 36.16 60 30 6 0.28 56.19 42.19 60 30 7 0.28 壳体装药部位及无装药部位的单位长度修正Payman破碎参数Cu与C/M比值的关系如图 8所示。可以看出，无装药部分的破碎参数Cu较小且变化不大，而壳体装药部分的破碎参数Cu与C/M比值之间依然存在线性关系，并且可用下式表示
Cu=4.8+91.7C/M (7) 
图 8 壳体破碎参数Cu与C/M比值的关系Figure 8. Variation of normalized Payman fragmentation parameter Cu with C/M ratio将(7)式与(5)式做对比发现，两者的截距及斜率几乎相等，因此可以近似认为壳体装药部分的破碎参数与C/M比值的关系和具有端部效应的壳体完全相同，即无装药部分的存在对装药部分的破碎性没有任何影响。对回收的破片仔细观察发现，壳体破碎时主要沿装药与无装药的连接处断裂，不产生既包括装药部分又包括无装药部分的破片，因此装药部分的破碎性与具有端部效应的壳体是等同的。
当装药部分的C/M比值和无装药部分长度L′保持不变时，壳体的破碎参数Cu随壁厚的变化关系如图 9所示。随着壳体壁厚的增加，装药部分及无装药部分的破碎性几乎保持不变，表明壳体壁厚改变时无装药部分对装药部分的破碎性无影响。对于无装药部分长度不同的壳体，其破碎性参数变化如图 10所示，可以看出，当无装药部分长度增加时，装药部分的破碎参数几乎不变。因此可得出，无装药部分长度不影响壳体装药部分的破碎性能。
图 9 壳体破碎参数Cu与壁厚的关系Figure 9. Variation of normalized Payman fragmentation parameter Cu with wall thickness
图 10 壳体破碎参数Cu与无装药部分长度的关系Figure 10. Variation of normalized Payman fragmentation parameter Cu with length of charge vacancy综合图 8、图 9和图 10还可以得出，无装药部分壳体的破碎参数事实上很小。当装药部分壳体向外膨胀时，无装药部分受影响而被“撕裂”成破片。由于两部分壳体向外膨胀的加速度相差悬殊，所以沿两部分的连接处断裂。因此可近似认为，装药部分壳体等效于存在端部效应的壳体。无装药部分的破碎性则主要取决于壳体材料本身性质，C/M比值即使有影响也不会很明显，而长度和壁厚的变化虽然会影响壳体质量和破片大小，但总破碎参数C0近似相等，因为C0是经过壳体质量(回收破片总质量)修正过的参数。
3. 结论
针对内爆炸载荷作用下金属壳体的动态碎裂问题，通过圆柱钢壳内爆加载试验研究了不同装药条件对破片质量分布的影响规律，得到如下主要结论。
(1) 对于不同的给定材料与炸药组合，壳体的单位长度修正Payman破碎参数Cu与C/M比值均存在线性关系，即Cu=a+bC/M，其中a、b为与壳体材料和炸药组合有关的特征常数；
(2) 无论壁厚和C/M比值如何变化，端部效应使壳体的破碎参数Cu降低至一个定值3a/4；
(3) 无装药部分对壳体的破碎性影响很小，当无装药部分长度或壁厚发生变化时，壳体的破碎参数Cu近似不变，且Cu与C/M比值的线性关系同样成立。
-
图 1 梯度泡沫铝和各向异性压缩测试试样
Figure 1. Functionally gradient aluminum foam specimens and anisotropic compression tests
图 2 均匀密度泡沫铝（a）、LC和TC（b）准静态压缩响应
Figure 2. Typical compressive response of uniform density aluminum foams (a)， LC and TC (b) under quasi-static compression
图 4 3种泡沫的压缩强度（a）和相对弹性模量（b）随相对密度的变化
Figure 4. Variation trends of strength (a) and relative elastic modulus (b) with relative density for three types of foams
图 5 3种泡沫的力学特性（a）和致密化应变（b）对比
Figure 5. Comparison of mechanical properties (a) and densification strain (b) for three types of aluminum foams
图 6 LC（a–d）和TC（e–f）试样准静态压缩典型力学行为
Figure 6. Typical quasi-static compression deformation behavior of LC (a–d) and TC (e–f) specimens
图 7 落锤冲击下LC试样变形行为（a–d）和应变场（e–h）
Figure 7. Deformation (a–d) and strain fields (e–h) of LC specimen under drop-weight impact
图 8 落锤冲击TC试样变形行为（a–c）和应变场（d–f）
Figure 8. Deformation (a–c) and strain fields (d–f) of TC specimen under drop-weight impact
图 10 数值模拟LC和TC的变形行为
Figure 10. Numerical simulation of deformation behavior for LC and TC samples
-
[1] 刘姗姗, 刘亚军, 张英杰, 等. 碳纤维-泡沫铝夹芯板低速冲击响应 [J]. 高压物理学报, 2020, 34(3): 034202.LIU S S, LIU Y J, ZHANG Y J, et al. Low-velocity impact response of carbon fiber-aluminum foam sandwich plate [J]. Chinese Journal of High Pressure Physics, 2020, 34(3): 034202. [2] 王永欢, 徐鹏, 范志强, 等. 球形孔开孔泡沫铝的力学特性及准静态压缩变形机制 [J]. 高压物理学报, 2019, 33(1): 014201.WANG Y H, XU P, FAN Z Q, et al. Mechanical characteristics and quasi-static compression deformation mechanism of open-cell aluminum foam with spherical cells [J]. Chinese Journal of High Pressure Physics, 2019, 33(1): 014201. [3] 张国栋, 李翔宇, 梁民族, 等. 内部爆炸载荷下泡沫铝夹心柱壳动态响应仿真研究 [J]. 高压物理学报, 2018, 32(2): 024101.ZHANG G D, LI X Y, LIANG M Z, et al. Dynamic response of sandwich cylinders cored with aluminum foam under internal blast loading [J]. Chinese Journal of High Pressure Physics, 2018, 32(2): 024101. [4] HE S, ZHANG Y, DAI G, et al. Preparation of density-graded aluminum foam [J]. Materials Science and Engineering: A, 2014, 618: 496–499. doi: 10.1016/j.msea.2014.08.087 [5] 杨旭东, 成莹, 郑远兴, 等. 梯度泡沫铝的制备及力学性能研究进展 [J]. 热加工工艺, 2019, 48(12): 16–20.YANG X D, CHENG Y, ZHENG Y X, et al. Research development of preparation and mechanical properties of graded aluminum foam [J]. Hot Working Technology, 2019, 48(12): 16–20. [6] HANGAI Y, MORITA T, UTSUNOMIYA T. Functionally graded aluminum foam consisting of dissimilar aluminum alloys fabricated by sintering and dissolution process [J]. Materials Science and Engineering A, 2017, 696: 544–551. doi: 10.1016/j.msea.2017.04.070 [7] XIA Y, WU C, LIU Z X, et al. Protective effect of graded density aluminum foam on RC slab under blast loading—an experimental study [J]. Construction and Building Materials, 2018, 111: 209–222. [8] 曹娇. 功能密度梯度泡沫铝及其填充结构的吸能性能研究[D]. 长沙: 湖南大学, 2015.CAO J. Study on crashworthiness of functionally graded density aluminum foam-filled tapered thin-walled structures [D]. Changsha: Hunan University, 2015. [9] ZHANG B, HU S, FAN Z. Anisotropic compressive behavior of functionally density graded aluminum foam prepared by controlled melt foaming process [J]. Materials, 2018, 11: 2470–2480. doi: 10.3390/ma11122470 [10] AMSTERDAM E, HOSSON J T M D, ONCK P R. On the plastic collapse stress of open-cell aluminum foam [J]. Scripta Materialia, 2008, 59(6): 653–656. doi: 10.1016/j.scriptamat.2008.05.025 [11] GIBSON L J, ASHBY M F. Cellular solids : structure and properties [M]. Cambridge, UK: Cambridge University Press, 1997. [12] MU Y, YAO G. Anisotropic compressive behavior of closed-cell Al–Si alloy foams [J]. Materials Science and Engineering: A, 2010, 527: 1117–1119. doi: 10.1016/j.msea.2009.09.045 [13] FAN Z, ZHANG B, GAO Y, et al. Deformation mechanisms of spherical cell porous aluminum under quasi-static compression [J]. Scripta Materialia, 2018, 142: 32–35. doi: 10.1016/j.scriptamat.2017.08.019 [14] AMSTERDAM E, HOSSON J T M, et al. Failure mechanisms of closed-cell aluminum foam under monotonic and cyclic loading [J]. Acta Materialia, 2006, 54: 4465–4472. doi: 10.1016/j.actamat.2006.05.033 -
下载:
下载:
点击查看大图
计量
- 文章访问数: 5199
- HTML全文浏览量: 1901
- PDF下载量: 34
