A theoretical model of PBXs’ tensile strength based on meso-structure parameters
-
摘要: 为给塑性黏结炸药(PBX)的力学强度设计提供支撑、探索材料细观特征量与材料强度之间的定量规律,应用微裂纹扩展区理论,将PBX炸药的单轴拉伸过程中力学响应特征的变化归结为扩展裂纹取向角度的增加,将扩展裂纹最大取向角与拉伸强度相关联,构建了基于材料细观特征量的拉伸强度理论模型,并采用不同温度的单轴拉伸实验验证了该理论模型的有效性。研究表明:该拉伸强度理论模型可以实现对PBX炸药拉伸强度与炸药微裂纹密度、颗粒/黏结剂界面性能以及颗粒/黏结剂体系的表观杨氏模量、泊松比等细观特征量之间关系的定量描述。Abstract: In practical applications, plastic bonded explosive (PBX) explosives are often used as load bearing structural components. Hence, the mechanical strength is an important parameter in PBXs’ design. It is of great interest how to adjust the microscopic characteristics of the material in the manufacturing process to obtain an PBX with the required strength. PBX consists of a large portion of energetic particles with a small portion of binder. Therefore, a large number of randomly distributed microcracks exist inside. Inspired by the successful application of microcrack extension mechanism in the PBXs’ mechanical response simulation in recent years, the theory related to microcrack extension was applied to the strength modeling study. The domain of microcrack growth (DMG) theory was applied to analyze PBXs’ uniaxial tension. The results shows that the increase in tensile load resulted in microcrack extension at different orientation angles, as well as a decrease in the crack interval. Meanwhile, recent studies on the cascading behavior of randomly distributed microcracks show that the crack cascading behavior with each other is controlled by the orientation angle and crack interval of adjacent cracks. The PBX tensile fracture comes from the instable extension of macroscopic cracks, and the macroscopic crack formation and extension both come from the cascading of microcracks. By treating the strength as the minimum load required for catastrophic crack cascade extension, there is a one-to-one correspondence between the strength and the maximum orientation angle of extended crack. Based on the DMG theory, this one-to-one correspondence is expressed as the theoretical model of tensile strength. This theoretical model establishes the connection between the tensile strength and the meso-structure parameters such as stress intensity factor, microcrack diameter and microcrack density. Comparison of tensile strength model predictions with wide temperature range (from –40 ℃ to 45 ℃) experimental data, indicated that this model is capable to describe the PBXs’ tensile strength. Suggestions for the design of strength enhancement of PBXs can be provided by this theoretical model.
-
Key words:
- PBX /
- tensile strength /
- theoretical model /
- meso damage /
- DMG theory
-
材料与结构的轻量化是科技发展的重要趋势,其中泡沫金属材料由于轻质化、孔隙率高和可变形量大等,被广泛应用于能量吸收和抗冲击减振的保护结构。为了提高泡沫材料的力学性能,研究人员提出了具有密度梯度的泡沫材料[1-3],以实现多目标优化设计。
由于材料制造技术的局限性,在实际应用中还不容易获得具有特定相对密度分布的梯度泡沫样本,也无法直接观察加载情况下胞元内部的变形情况,研究人员转而应用有限元方法研究梯度泡沫材料的力学响应[4-9]。研究结果显示,冲击载荷下梯度泡沫材料存在3种变形模态[10-12],且具有特定密度梯度的泡沫金属可以降低受保护结构的最大冲击应力,改善材料的能量吸收和抗冲击性能[13-19]。Zeng等[20]测试了密度梯度空心球金属泡沫的动态响应,将最硬层作为第1冲击层,将最弱层与受保护结构接触,可以减少传输力并获得更多的吸收能量。Ajdari等[21]研究了具有规则和不规则排列的功能梯度多胞结构的动态压溃,表明沿加载方向降低相对密度,可以增强蜂窝在压溃初期的能量吸收。Zhang等[22]建立了分层递变梯度蜂窝模型,建议把最弱的层放置在冲击端或输出端,最强的层被放置在中间阶段,以达到更高的能量吸收效率。Zheng等[23]从理论和数值计算两方面发现,均匀泡沫的能量吸收效果可能最好,密度梯度较大的正、中高梯度泡沫可以在撞击端保护物体,密度梯度较小的负梯度泡沫可以在支撑端保护物体,而大梯度的中低、中高、正梯度泡沫可以提高吸能效率。
目前,对梯度泡沫金属材料的研究主要集中在密度单向变化,为了进一步研究梯度泡沫金属材料密度在非单向变化下的变形及能量吸收机理,我们设计了一种密度径向连续排布的梯度泡沫材料,研究它在冲击载荷作用下的动态力学性能。首先,利用3D-Voronoi技术构建密度排布不同的层状、径向梯度泡沫结构;然后,导入有限元分析软件中生成梯度泡沫金属的冲击模型,观察他们在不同恒定冲击速度下的变形特点,比较他们的抗冲击性能和能量吸收能力,考察不同梯度泡沫材料在工程防护上的优劣;最后,分析密度梯度和相对密度参数对径向梯度泡沫材料力学性能的影响。
1. 数值模拟
1.1 三维随机Voronoi技术
3D-Voronoi模型是将N个成核点完全随机分布在体积为V0的立方体区域内生成,为了防止Voronoi模型产生过小胞元,立方体内任意两相邻核点之间的最小距离为:
δ=(1−k)d0=(1−k)√62(V0√2N)1/3 (1) 式中:无量纲变量k为不规则度,d0为正十四面体模型中任意两相邻成核点间的最小距离[24]。将成核点复制到周围区域内,编程生成Voronoi构型,并通过切割算法得到原先指定区域中的Voronoi结构作为最后的计算模型。
1.2 密度梯度
为了研究密度排布对泡沫金属力学性能的影响,基于密度沿着厚度方向上单向变化的层状梯度泡沫模型(见图1(a)),设计了一种径向梯度泡沫模型,即以正方形中心点出发密度沿着半径方向渐进变化,如图1(b)所示。
层状梯度泡沫的密度梯度γ可描述为:
ΔμΔL=γ¯μL (2) 由此,推演径向梯度泡沫密度梯度γ的定义为:
ΔμΔV=γ¯μV (3) 式中:γ为无量纲参数,
¯μ 为材料的平均相对密度,Δμ=μi+1−μi,μi为第i层的相对密度,L为沿梯度方向的总长度,V为泡沫模型的总体积。将泡沫模型平均分为N等份,ΔL、ΔV分别表示每个等份的长度和体积。当γ>0时,为径向正梯度泡沫,表示从冲击端至中心处密度逐渐增加,相反则为径向负梯度泡沫。曹国剑等[25]介绍了利用离心沉积法制备梯度泡沫铝的方法,可直接用来生产本文中提出的径向负梯度泡沫。
1.3 有限元模型
表1为层状正、负梯度与径向正、负梯度4种泡沫金属模型的材料参数。4种泡沫金属的冲击加载模型如图2所示,垂直于所示平面方向泡沫密度均匀,冲击加载沿梯度方向进行。模型尺寸为50 mm×50 mm×25 mm,不规则度为0.3,有限元模型的核点个数为900,胞壁厚度由给定的相对密度确定。在有限元收敛分析后,泡沫模型使用S3R和S4R壳单元进行混合网格划分,网格尺寸设置为0.2 mm。在左端固定一个刚性板,右端的刚性板以给定的恒定速度撞击泡沫模型。基体材料设为双线性应变强化模型,密度ρ=2 700 kg/m3,弹性模量E=69 GPa,泊松比ν=0.3,屈服应力σy=76 MPa,切线模量为0.69 GPa。泡沫模型定义为自接触,两端刚性板和泡沫模型之间定义为面面接触,设定动摩擦因数为0.2。
表 1 模型材料参数Table 1. Model material parameters泡沫模型 平均相对密度¯μ 密度梯度γ 层状正梯度 0.12 0.8 层状负梯度 0.12 −0.8 径向正梯度 0.12, 0.09 0.8, 0.4 径向负梯度 0.12, 0.09 −0.8, −0.4 2. 结果与讨论
2.1 变形特点和名义应力应变曲线
根据梯度泡沫在不同冲击速度下的准静态模态、过渡模态和冲击模态[26-27],将冲击速度设置为30、80和200 m/s。其中,层状正、负梯度泡沫的变形特点与以往文献描述的相似。图3为径向梯度泡沫在不同冲击速度下的应变云图。
当冲击速度较小 (v=30 m/s) 时,变形属于准静态模态。径向正梯度泡沫先在冲击端变形(低密度),随后支撑端(低密度)的胞元发生变形。这是由于,当速度相对较低时,泡沫内部处于均匀的应力状态,两端首先达到准静态屈服应力。正梯度泡沫中间区域的密度较高,且模型结构并不完全对称,在ε>0.5后出现沿支撑板两个切线方向滑移现象。径向负梯度泡沫因中间密度较小,变形先从中间区域处产生,两端在后期同样出现滑移现象。
当冲击速度中等(v=80 m/s)时,变形受准静态屈服应力和惯性效应两个因素的影响。径向正梯度泡沫的变形特点属于典型的过渡模态,而径向负梯度泡沫没有表现明显的过渡模态。负梯度泡沫先在冲击端发生变形,随后中间区域变形,并逐渐向支撑端传播,此现象与冲击模态类似。发现在冲击速度为60 m/s时,径向负梯度泡沫处于过渡模态。冲击端先发生少量变形,接着中间区域持续发生变形,待中间区域压实后,变形继续向支撑端传播,直到支撑端完全压实。从变形特点可以看出,径向负梯度泡沫比径向正梯度泡沫在更低速度下进入冲击模态。
当冲击速度较高(v=200 m/s)时,惯性作用成为了控制变形的主导因素,其影响效果远超过准静态屈服应力,泡沫进入冲击模态。径向正梯度泡沫在冲击端发生大量变形后,仅在支撑端发生少量变形;而径向负梯度泡沫只在冲击端发生变形,并逐渐变形至压实。
图4为不同冲击速度下4种泡沫冲击端与支撑端的名义应力应变曲线。随着冲击速度的增大,4种泡沫冲击端应力与压实应变逐渐增大,支撑端比冲击端越来越晚达到初始屈服应力,且应力平台区越来越长,应力值差异逐渐明显。层状正梯度泡沫的冲击端初始峰值应力较小,到达峰值应力后曲线逐渐上升,但支撑端曲线应力值高于层状负梯度泡沫,初始屈服应力较大。层状负梯度泡沫的冲击端初始峰值应力较大,到达峰值应力后曲线逐渐下降。层状梯度泡沫的数值模拟结果与其他文献描述相似,佐证了有限元模型的合理性。
与其他3种泡沫相比,径向正梯度泡沫由于两端密度小、相对较“软”,冲击端初始峰值应力和支撑端初始屈服应力均最低。径向正梯度泡沫中间区域密度大、强度高,冲击端曲线在应变0.2~0.6区间时先上升后下降。径向负梯度泡沫因两端密度大、中间区域密度小,冲击端曲线在应变0.2~0.6区间时先下降后上升,而支撑端曲线始终保持基本稳定。
图5为4种梯度泡沫冲击端与支撑端的最大应力和标准差。最大应力表示泡沫达到密实化应变前冲击端或支撑端所承受的最高载荷;标准差表示数据之间的离散程度,离散程度越小,曲线波动幅度越小。通过比较发现:冲击速度越高,4种泡沫冲击端的最大应力与波动幅度呈整体增加趋势;径向正梯度泡沫冲击端的最大应力在4种泡沫中最小;径向负梯度泡沫的两端应力波动幅度明显小于其他3种,支撑端应力尤其稳定。
2.2 能量吸收能力
能量吸收能力是判断泡沫材料抗冲击性能的一个典型指标。可以用材料压缩至某个变形量时所吸收的能量来衡量吸收能力的强弱,定义为:
W=∫ε0σ(ε)dε (4) 图6和表2为4种梯度泡沫能量吸收能力,可以看出,几种泡沫的能量吸收能力在变形过程中呈现交替变化的现象。
表 2 4种梯度泡沫的能量吸收能力Table 2. Energy absorption capacities of four graded foams泡沫模型 W/MPa v=30 m/s v=80 m/s v=200 m/s ¯ε=0.2 ¯ε=0.5 ¯ε=0.8 ¯ε=0.2 ¯ε=0.5 ¯ε=0.8 ¯ε=0.2 ¯ε=0.5 ¯ε=0.8 层状正梯度 0.55 1.98 4.22 0.82 2.73 5.36 2.36 7.18 13.51 层状负梯度 0.73 1.96 4.24 1.72 3.14 4.78 4.41 9.43 12.67 径向正梯度 0.62 2.00 3.97 1.00 3.06 4.91 2.84 8.34 12.63 径向负梯度 0.91 2.25 4.22 1.35 3.14 5.36 3.74 8.60 13.89 当速度为30 m/s时:在压缩初期,径向负梯度泡沫表现出优异的能量吸收能力;在压缩后期,径向负梯度泡沫与层状正、负梯度泡沫几乎相同,而径向正梯度泡沫能量吸收值相对较小。当速度提高到80和200 m/s时:在压缩初期,层状、径向负梯度泡沫的能量吸收能力大于层状、径向正梯度泡沫;在压缩后期,径向负梯度泡沫的能量吸收值仍较高,层状正梯度泡沫的能量吸收能力增长加快。
比较4种梯度泡沫的力学响应发现,径向正梯度泡沫不仅拥有层状正梯度泡沫冲击端应力较小、保护冲击端冲击物体的特点,而且继承了层状负梯度泡沫支撑端应力较小、保护支撑端受冲击物体的特点。径向负梯度泡沫的两端应力波动幅度较小,曲线变化趋势较为平缓,能够保证被保护物体受力稳定,同时它的能量吸收能力始终处于较高水平,在低速冲击或大变形工况下是最优的吸能结构。因此,径向梯度泡沫是缓冲减压吸能装置的良好选择。
2.3 径向密度梯度的影响
为了分析径向正、负梯度泡沫对密度梯度的敏感程度,图7~8分别给出了
|γ| =0.4、0.8时的应力应变曲线和能量吸收曲线。从图7可以看出:对于冲击端,当密度梯度改变时,泡沫初始峰值应力和密实化时的应力变化明显;对于支撑端,径向正梯度泡沫在两种密度梯度下的曲线出现分离,而负梯度泡沫曲线几乎相同。对于径向正梯度泡沫,γ越大,冲击端初始峰值应力与支撑端初始屈服应力越小;对于径向负梯度泡沫,当
|γ| 增加,冲击端初始峰值应力与支撑端初始屈服应力增大。从图8可以看出,泡沫的能量吸收能力对密度梯度大小不敏感、对梯度方向敏感。当密度梯度改变时,两种梯度下的泡沫所吸收的能量几乎相同;在达到完全密实化前,径向负梯度泡沫的能量吸收能力始终大于径向正梯度泡沫。由于试样具有相同的平均相对密度,不同梯度泡沫完全密实化后所吸收的总能量基本相同。
2.4 平均相对密度的影响
图9为不同平均相对密度下径向正、负梯度泡沫的名义应力应变曲线。可以发现,相对密度对金属泡沫的冲击响应具有显著影响,其影响程度远大于密度梯度。当相对密度
¯μ 由0.09上升到0.12时,泡沫金属的两端应力值都明显增大,而应力变化趋势并不发生改变。引入单位质量材料吸收的能量w和能量吸收效率ηw [28],分别为:
w=W¯μρ (5) ηw=wwf (6) 式中:wf是单位质量泡沫的总吸收能量。由式(5)~(6),不同密度下泡沫材料吸收的能量如图10所示。随着相对密度增大,径向梯度泡沫能量吸收增长趋势加快,能量吸收效率提高,材料的吸能效果更好;冲击速度越大,w和ηw越大,泡沫材料初始冲击阶段所吸收的能量越多。
当v=30 m/s、
¯ε =0.5时:¯μ =0.12的径向负梯度泡沫吸能最多,ηw为0.240;¯μ =0.09的径向正梯度泡沫吸能最少,ηw为0.225。当v=30 m/s、¯ε >0.5时,¯μ =0.12的径向正梯度泡沫吸能超过¯μ =0.09的径向负梯度泡沫。当v=80 m/s、
¯ε =0.5时,¯μ =0.12的径向负梯度泡沫ηw上升到0.404,¯μ =0.09的径向正梯度泡沫ηw上升到0.357。当v=80 m/s、¯ε >0.35时,¯μ =0.12的径向正梯度泡沫与¯μ =0.09的径向负梯度泡沫的吸能优势发生交换。当v=200 m/s时,两种密度下泡沫的能量吸收差异性逐渐消失,径向负梯度泡沫的能量吸收值始终大于径向正梯度泡沫。
3. 结 论
利用3D-Voronoi技术设计了层状正、负梯度与径向正、负梯度4种密度排布不同的梯度泡沫模型,通过有限元分析软件分别从变形特点、载荷响应历程及能量吸收机理比较了4种梯度泡沫金属在恒速冲击载荷下的动态响应结果,同时探讨了径向梯度泡沫中密度梯度、相对密度参数对力学性能的影响,得到如下结论。
(1)径向正梯度泡沫与层状正、负梯度泡沫相比冲击端与支撑端的应力值均处于较低水平,不仅保护冲击端冲击物体,而且保护支撑端受冲击物体;径向负梯度泡沫两端应力应变曲线平缓,波动幅度较小,能够保证被保护物体受力稳定。
(2)在不同冲击速度下,4种梯度泡沫金属的能量吸收能力交替变化,其中径向负梯度泡沫能量吸收能力始终处于较高水平,在低速冲击或大变形工况下是最优的吸能结构。
(3)能量吸收能力对密度梯度大小不敏感、对梯度方向敏感,在达到完全密实化前,径向负梯度泡沫的能量吸收能力始终大于径向正梯度泡沫;平均相对密度越大,径向正、负梯度泡沫两端应力越大、吸能效果越好。
-
表 1 与温度无关的PBX-3细观特征量
Table 1. Temperature independent mesoscopic characteristics of PBX-3
n/cm−3 a0/cm au/cm v 1372 0.006352 0.03703 0.3 表 2 与温度相关的PBX-3细观特征量
Table 2. Temperature dependent mesoscopic characteristics of PBX-3
T/℃ KⅠC/(MPa·cm1/2) KⅡC/(MPa·cm1/2) −40 0.1628 0.1911 −30 0.1502 0.1763 −25 0.1556 0.1827 −20 0.1520 0.1784 −10 0.1457 0.1711 25 0.1457 0.1711 45 0.1322 0.1552 -
[1] STEVENS R. A strength model and service envelope for PBX 9501: LA-UR-14-20696 [R]. Los Alamos: Los Alamos National Laboratory, 2014. [2] 肖磊, 刘杰, 郝嘎子, 等. 微纳米RDX颗粒级配对压装PBX性能影响 [J]. 含能材料, 2016, 24(12): 1193–1197. DOI: 10.11943/j.issn.1006-9941.2016.12.011.XIAO L, LIU J, HAO G Z, et al. Effects of nano-/micrometer RDX particle gradation on the property of PBX [J]. Chinese Journal of Energetic Materials, 2016, 24(12): 1193–1197. DOI: 10.11943/j.issn.1006-9941.2016.12.011. [3] 黄辉. 颗粒级配技术及其在含能材料中的应用 [J]. 含能材料, 2001, 9(4): 161–164. DOI: 10.3969/j.issn.1006-9941.2001.04.005.HUANG H. Particle grade technique and application on energetic materials [J]. Energetic Materials, 2001, 9(4): 161–164. DOI: 10.3969/j.issn.1006-9941.2001.04.005. [4] LV K Z, YANG K, ZHOU B, et al. The densification and mechanical behaviors of large-diameter polymer-bonded explosives processed by ultrasonic-assisted powder compaction [J]. Materials & Design, 2021, 207: 109872. DOI: 10.1016/j.matdes.2021.109872. [5] 黄辉, 王晓川. 偶联剂在HMX基浇注固化炸药中的作用 [J]. 含能材料, 2000, 8(1): 13–17. DOI: 10.3969/j.issn.1006-9941.2000.01.004.HUANG H, WANG X C. Behavior of coupling agent in HMX-based extrudable cast explosive [J]. Energetic Materials, 2000, 8(1): 13–17. DOI: 10.3969/j.issn.1006-9941.2000.01.004. [6] LI F, YE L, NIE F D, et al. Synthesis of boron-containing coupling agents and its effect on the interfacial bonding of fluoropolymer/TATB composite [J]. Journal of Applied Polymer Science, 2007, 105(2): 777–782. DOI: 10.1002/app.26117. [7] YANG X L, GONG F Y, ZHANG K, et al. Enhanced creep resistance and mechanical properties for CL-20 and FOX-7 based PBXs by crystal surface modification [J]. Propellants, Explosives, Pyrotechnics, 2021, 46(4): 572–578. DOI: 10.1002/prep.202000277. [8] 冯西桥, 余寿文. 准脆性材料细观损伤力学 [M]. 北京: 高等教育出版社, 2002: 20–41.FENG X Q, YU S W. Damage micromechanics of quasi-brittle materials [M]. Beijing: Higher Education Press, 2002: 20–41. [9] YANG Y Z, LI J L, CHEN J K. Experimental and numerical analysis of shear process of a high particle content bonding material [J]. Journal of Applied Mechanics, 2021, 88(6): 061003. DOI: 10.1115/1.4050277. [10] LIU C, THOMPSON D G. Crack initiation and growth in PBX 9502 high explosive subject to compression [J]. Journal of Applied Mechanics, 2014, 81(10): 101004. DOI: 10.1115/1.4028087. [11] LIU M, HUANG X C, WU Y Q, et al. Numerical simulations of the damage evolution for plastic-bonded explosives subjected to complex stress states [J]. Mechanics of Materials, 2019, 139: 103179. DOI: 10.1016/j.mechmat.2019.103179. [12] CURRAN D R, SEAMAN L, SHOCKEY D A. Dynamic failure of solids [J]. Physics Reports, 1987, 147(5/6): 253–388. DOI: 10.1016/0370-1573(87)90049-4. [13] BENNETT J G, HABERMAN K S, JOHNSON J N, et al. A constitutive model for the non-shock ignition and mechanical response of high explosives [J]. Journal of the Mechanics and Physics of Solids, 1998, 46(12): 2303–2322. DOI: 10.1016/s0022-5096(98)00011-8. [14] ZUO Q H, ADDESSIO F L, DIENES J K, et al. A rate-dependent damage model for brittle materials based on the dominant crack [J]. International Journal of Solids and Structures, 2006, 43(11/12): 3350–3380. DOI: 10.1016/j.ijsolstr.2005.06.083. [15] RAGASWAMY P, LEWIS M W, LIU C, et al. Modeling the mechanical response of PBX 9501: LA-UR-10-01416 [R]. Los Alamos: Los Alamos National Laboratory, 2010. [16] PIJAUDIER-CABOT G, BAŽANT Z P. Nonlocal damage theory [J]. Journal of Engineering Mechanics, 1987, 113(10): 1512–1533. DOI: 10.1061/(asce)0733-9399(1987)113:10(1512. [17] LIU C. Damage and fracture in high explosive materials: LA-UR-11-05618 [R]. Los Alamos: Los Alamos National Laboratory, 2011. [18] ZHANG S L, YANG W. Macrocrack extension by connecting statistically distributed microcracks [J]. International Journal of Fracture, 1998, 90(4): 341–353. DOI: 10.1023/A:1007440012284. 期刊类型引用(5)
1. 魏建辉,李旭,黄威,徐宏建,方搏鹏. 高速冲击载荷下梯度金属泡沫夹芯梁的动态响应与失效. 爆炸与冲击. 2023(05): 101-111 . 本站查看
2. 张健. 泡沫金属材料力学性能及能量分析. 产品可靠性报告. 2023(10): 68-70 . 百度学术
3. 陆飞,邓安仲,孙皓,段承军,牟泽龙. 能量耗散型抗冲击材料研究进展. 当代化工. 2022(03): 666-670 . 百度学术
4. 刘英志,何玉龙,王志华. 星形衍生多胞薄壁管的动态压缩. 应用力学学报. 2022(04): 758-765 . 百度学术
5. 贺璞,邓庆田,李新波. 层合多孔圆柱壳的轴向冲击吸能特性. 高压物理学报. 2022(04): 145-154 . 百度学术
其他类型引用(4)
-