2. 中国工程物理研究院化工材料研究所,四川 绵阳 621999;
3. 西南科技大学,四川 绵阳 621010;
4. 中国工程物理研究院研究生部,四川 绵阳 621999
2. Institute of Chemical Materials, CAEP, Mianyang 621999, China;
3. Southwest University of Science and Technology, Mianyang 621010, China;
4. Graduate School of China Academy of Engineering Physics, Mianyang 621999, China
炸药部件在武器系统中除需提供能量外,还要承受装配、服役环境变化等产生的力学载荷,美国Los Alamos实验室的研究[1]表明,破坏的炸药部件不仅影响武器系统的战场表现,更有可能诱发意外爆炸。目前炸药工程中常采用单轴强度准则(Uniaxial-strength,认为某点的最大拉/压应力超过其单轴拉伸/压缩强度时材料破坏,又称最大主应力准则或第一强度准则[2])评估炸药部件的承载能力,实践中炸药部件的意外破坏事件表明了单轴强度准则的局限性,研究高聚物粘结炸药(Polymer Bonded Explosive,PBX)复杂应力下的强度特性及其适用准则对于保障武器的安全性和可靠性具有现实意义。
强度问题命题简单,求解复杂,自15世纪达芬奇和伽利略研究铁丝和石料的单轴拉伸强度[2]以来,人们迄今已针对不同材料提出了“上百个模型或准则”[3],遗憾的是还未有满足所有材料各种应力状态的准则出现。俞茂宏[4]根据强度准则描述中主剪应力数目的多寡将众多的准则分为单剪、双剪和三剪三个系列,该方法在业界获得了广泛的认同。炸药材料的强度问题系关炸药部件的结构完整性和服役安全性,是破坏力学研究和结构强度评估的基础,问题的重要性引起了业内科研人员的关注,但由于多轴加载手段的缺失,目前的研究大都只能基于单轴加载手段进行。Li等[5]采用单剪Mohr-Coulomb准则[6]分析了三种炸药材料的单轴动态破坏试验,结果表明仅有一种材料的分析结果满足工程需要;傅华[7]、李俊玲[8]基于霍普金森杆研究了单轴应力下PBX破坏强度等参量的率相关性;蓝林钢[9]、赵玉刚[10]分别就单轴动态拉伸时PBX的力学响应和强度特性进行了研究。上述工作有助于认识PBX单轴应力下的强度特性,而对于复杂应力下的强度问题的指导意义不大,事实上,目前尚未系统开展过PBX的多轴力学性能研究,人们对复杂应力下PBX的强度认识还较初浅,不能回答炸药部件装配预紧等复杂应力环境下的强度特性和破坏行为,笔者基于端部约束压溃试验获得了三轴压缩应力下3种典型强度准则描述PBX强度特性的适应性[11],而实际工程中拉伸应力主导破坏的风险远高于压缩应力,这是因为PBX单轴拉伸强度常不足单轴压缩强度的三分之一。本研究基于单轴加载技术,依据端部约束拉溃试验,并在先前研究[11]的基础上采用试验测试结合数值模拟的方法,分析Uniaxial-strength[2]、Mohr-Coulomb[6]和Twin-shear[3]三种典型强度准则描述某TATB基PBX三轴拉伸应力下强度特性的适应性。
2 三轴拉伸应力破坏试验 2.1 试验设计研究对象为某TATB基PBX,试验温度(20±2) ℃,试验测得该批次试件的杨氏模量为7.37 GPa,单轴拉伸强度为8.83 MPa,单轴压缩强度27.61 MPa,泊松比0.365,图 1为20 ℃下三发PBX样品单轴拉伸时的应力应变曲线。
研究中基于单轴材料试验机通过调节试样边界条件实现复杂应力下样品的破坏,试验设计遵循两个基本原则:一是试件形状和边界设置尽可能简单,以便于测试;二是试件破坏时的应力尽可能复杂。为此,选取不锈钢(1Cr18Ni9Ti,杨氏模量206 GPa,泊松比0.3)为工装材料,其模量为TATB基PBX的27.95倍,边界约束时可视为刚性约束;不同约束方式的有限元计算表明,端部约束条件下的拉伸破坏试验应力场既满足试验要求(详见3.1),又便于实现;进一步的探索试验表明,试件同平头底座间由于粘接力局限不能造成试件的破坏,经改进,最终确定了以Φ20 mm×20 mm圆柱体试件为测试对象,粘接头开槽(深度3 mm直径Φ20 mm)粘接(图 2);试验中采用粘接性能优良且与该PBX相容性良好的自制配方粘接胶,粘接后样品端部及嵌入的3 mm柱面段均可视为边界约束。
试验在INSTRON5582型材料试验机上进行,拉溃试验过程中5发样品的载荷-时间曲线如图 3所示,轴向破坏载荷数值见表 1,由表 1可知轴向破坏载荷的均值为841.60 N;试件的破坏面形态如图 4(因仅三发样品的破坏面较清晰,故未全部列出),最先破坏的位置发生在试件与工装的粘接处,整个破坏面呈现较为规则的圆锥状,依据试验样品正视图的测试结果,破坏面的水平夹角约14°~15°。
试件尺寸为Φ20 mm×20 mm,边界约束条件为距试件端部3 mm以内的区域为刚性约束。如图 1所示,该材料拉伸应力下的线弹性特征明显,研究中采用线弹性本构、轴对称PLANE182单元对临界轴向载荷下(841.60 N)样品内部的应力场进行计算,获得的三向主应力分布如图 5,有限元模型中相邻节点间距1 mm。
由图 5可知,最大第一、第二主应力出现在试件与工装的粘接处,最大数值分别为7.554 MPa和2.977 MPa(数值为正表示拉应力);最大第三主应力出现在试件端部与工装接触的区域,数值在1.08~1.38 MPa之间,最先破坏区域单元的三向主应力如表 2,该单元的节点编号如图 5c右上,最先破坏位置的节点编号为56,它所处的应力状态为三轴拉伸应力;另外结合图 4和图 5可以发现,破坏面的水平夹角与三向主应力分布存在关联,研究TATB基PBX复杂应力下的强度问题需要兼顾三向主应力的综合作用。
现有研究[11]获得了三轴压缩应力下该TATB基PBX的强度特性规律:Uniaxial-strength准则(式(1))的理论预测值低于试验值最多,安全阈度最高;Mohr-Coulomb准则(式(2))可用于表征表观强度下限,Twin-shear准则(式(3))可用于描述表观强度上限,式中σt、σc及α分别表示单轴拉伸强度、单轴压缩强度及二者之比,σt和σc的单位为MPa,均取正值,α为无量纲数。本研究基于端部约束拉溃试验数据,考察上述三个强度准则描述三轴拉伸应力破坏的适应性,描述误差修正为公式(4)中的算法。
$ {\sigma _1} = {\sigma _{\rm{t}}}\;\;\;{\rm{or}}\;\;\;\;-{\sigma _3} = {\sigma _{\rm{c}}} $ | (1) |
$ {\sigma _1}-\frac{{{\sigma _{\rm{t}}}}}{{{\sigma _{\rm{c}}}}}{\sigma _3} = {\sigma _{\rm{t}}} $ | (2) |
$ \left\{ \begin{array}{l} {\sigma _1}-\frac{\alpha }{2}\left( {{\sigma _2} + {\sigma _3}} \right) = {\sigma _{\rm{t}}}, \;\;\;\;\left( {{\sigma _2} \le \frac{{{\sigma _1} + \alpha {\sigma _3}}}{{1 + \alpha }}} \right)\\ \frac{1}{2}\left( {{\sigma _1} + {\sigma _2}} \right)-\alpha {\sigma _3} = {\sigma _{\rm{t}}}, \;\;\;\left( {{\sigma _2} \ge \frac{{{\sigma _1} + \alpha {\sigma _3}}}{{1{\rm{ + }}\alpha }}} \right)\\ \alpha {\rm{ = }}\frac{{{\sigma _{\rm{t}}}}}{{{\sigma _{\rm{c}}}}} \end{array} \right. $ | (3) |
$ {\rm{Representation\;\;error}} = \frac{{{\sigma _{1{\rm{theory}}}}-{\sigma _{1{\rm{test}}}}}}{{{\sigma _{1{\rm{test}}}}}} \times 100\% $ | (4) |
依据表 2的数据,采用Uniaxial-strength准则(式(1))描述时,σ1theory为单轴拉伸强度,8.83 MPa,理论预测的描述误差较试验值偏大,为16.90%。
$ {\sigma _{{\rm{1theory}}}} = {\sigma _{\rm{t}}} = 8.83 $ |
$ \begin{array}{l} {\rm{Representation\;\;error}} = \frac{{8.83-7.5536}}{{7.5536}} \times 100\% \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; = 16.90\% \end{array} $ |
采用Mohr-Coulomb准则(式(2))描述时,σ1theory为9.0353 MPa,描述误差中理论值较试验值偏大,为19.62%。
$ \alpha = \frac{{{\sigma _{\rm{t}}}}}{{{\sigma _{\rm{c}}}}} = \frac{{8.83}}{{27.61}} = 0.319812 $ |
$ \begin{array}{l} {\sigma _{{\rm{1theory}}}} = {\sigma _{\rm{t}}} + \frac{{{\sigma _{\rm{t}}}}}{{{\sigma _{\rm{c}}}}}{\sigma _3}\\ \;\;\;\;\;\;\;\;\; = 8.83 + 0.319812 \times 0.6419\\ \;\;\;\;\;\;\;\;\; = 9.0353 \end{array} $ |
$ \begin{array}{l} {\rm{Representation\;\;error}} = \frac{{9.0353-7.5536}}{{7.5536}} \times 100\% \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; = 19.62\% \end{array} $ |
采用Twin-shear准则(式(3))描述时,根据该准则的使用判据,应使用第一个描述式,此时的σ1theory为9.4087 MPa,理论值比试验值偏大24.56%。
$ \begin{array}{l} \frac{{{\sigma _1} + \alpha {\sigma _3}}}{{1 + \alpha }} = \frac{{7.5536 + 0.319812 \times 0.6419}}{{1 + 0.319812}}\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\; = 5.8788 \ge {\sigma _2} = 2.9774 \end{array} $ |
$ \begin{array}{l} {\sigma _{{\rm{1theory}}}} = {\sigma _{\rm{t}}} + \frac{\alpha }{2}\left( {{\sigma _2} + {\sigma _3}} \right)\\ \;\;\;\;\;\;\;\;\; = 8.83 + \frac{{0.319812}}{2}\left( {2.9774 + 0.6419} \right)\\ \;\;\;\;\;\;\;\;\; = 9.4087 \end{array} $ |
$ \begin{array}{l} {\rm{Representation\;\;error}} = \frac{{9.4087-7.5536}}{{7.5536}} \times 100\% \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; = 24.56\% \end{array} $ |
由上可知,TATB基PBX三轴拉伸应力下的强度特性较三轴压缩应力下的强度特性[11]存在明显差异:炸药工程中现用的Uniaxial-strength准则在描述三轴压缩应力下的破坏时,描述误差的理论值较试验值偏小43.77%(式(4)),体现出了数个强度准则中最好的安全阈度,而在描述三轴拉伸应力破坏时,描述误差的理论值却较试验值偏大16.90%,表明最为保守的Uniaxial-strength准则也不能有效地预测三轴拉伸应力下炸药材料的许用强度上限;Twin-shear准则和Mohr-Coulomb准则可分别用于描述三轴压缩应力下TATB基PBX的表观强度上下限,但在三轴拉伸应力下各自描述的理论值均大于试验值,最多偏大24.56%,三种强度准则的描述精度从高到低的排序为Uniaxial-strength准则、Mohr-Coulomb准则和Twin-shear准则。
试验中三轴拉伸主应力的比例为11.768:4.638:1,此时的最大主应力仅为单轴拉伸强度的0.831倍(σ1test=0.831σt),GB50010-2010[12]规定混凝土结构设计中三轴拉伸应力状态下的强度不能超过0.9倍单轴拉伸强度σt,可见三轴拉伸强度较单轴拉伸强度低在工程材料中具有普适性,目前炸药材料尚没有类似的标准,建立该标准还需要开展更多的工作。
另外,强度理论的研究表明[2-5],Twin-shear准则和单剪Mohr-Coulomb准则分别处于强度体系的上下限,材料的强度特性不仅受应力状态影响,还与静水压力(σm,又称平均应力,等于三个主应力的平均值)有关,因此,本研究中Mohr-Coulomb准则的表观强度低于试验测试结果的可能关键原因是静水压力的影响,建立该TATB基PBX高描述精度的强度准则除了要考虑应力状态外,还须考虑这一影响。
4 结论(1) 现炸药工程中使用的Uniaxial-strength准则的安全阈度相对最高,优于Mohr-Coulomb准则和Twin-shear准则,但仍不能安全预测三轴拉伸应力下的TATB基PBX破坏,试验测试结果较理论预测值偏大约16.90%。
(2) 描述精度方面,三种强度准则的描述精度由高到低的排序(误差)为:Uniaxial-strength准则(16.90%)>Mohr-Coulomb准则(19.62%)>Twin-shear准则(24.56%),误差均为试验值偏大理论值。
(3) 三轴拉伸应力(主应力比例11.768:4.638:1)破坏时的最大拉应力σ1为0.831σt,炸药部件的结构设计需考虑这种关系,建立TATB基PBX描述精度高的强度准则须同时考虑应力状态和静水压力的影响。
[1] |
Scott I Jackson, Larry G. Runaway reaction due to gas-dynamic choking in solid explosive containing a single crack[J].
Proceedings of the Combustion Institute, 2009, 32(2): 2307-2313. DOI:10.1016/j.proci.2008.05.089 |
[2] |
Mao-hong Y U. Advances in strength theories for materials under complex stress state in the 20th Century[J].
Appl Mech Rev, 2002, 55(3): 169-218. DOI:10.1115/1.1472455 |
[3] |
俞茂宏, YoshimineM, 强洪夫, 等. 强度理论的发展和展望[J].
工程力学, 2004, 21(6): 1-20. YU Mao-hong, Yoshimine M, QIANG Hong-fu, et al. Advances and prospects for strength theory[J]. Engineering Mechanics, 2004, 21(6): 1-20. |
[4] |
俞茂宏. 强度理论百年总结[J].
力学进展, 2004, 34(4): 529-560. YU Mao-hong. Advances in strength theories for materials under complex stress state in the 20th century[J]. Advances in Mechanics, 2004, 34(4): 529-560. DOI:10.6052/1000-0992-2004-4-J2002-054 |
[5] |
Jun-ling LI, Fang-yun LI. Effects of temperature and strain rate on the dynamic responses of three polymer-bonded explosives[J].
Journal of Strain Analysis, 2011, 47(2): 104-112. |
[6] |
俞茂宏, 昝月稳, 范文, 等. 20世纪岩石强度理论的发展——纪念Mohr-Coulomb强度理论100周年[J].
岩石力学与工程学报, 2000, 19(5): 545-550. YU Mao-hong, ZAN Yue-wen, FAN Wen, et al. Advances in strength theory of rock in 20 century 100 years in memory of the Mohr-Coulomb strength theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2000, 19(5): 545-550. |
[7] |
傅华, 李俊玲, 谭多望. PBX炸药本构关系的实验研究[J].
爆炸与冲击, 2012, 32(3): 231-236. FU Hua, LI Jun-ling, TAN Duo-wang. Experimental study on constitutive relations for plastic-bonded explosives[J]. Explosion and Shock Wave, 2012, 32(3): 231-236. DOI:10.11883/1001-1455(2012)03-0231-06 |
[8] |
李俊玲, 卢芳云, 傅华, 等. 某PBX炸药的动态力学性能研究[J].
高压物理学报, 2011, 25(2): 159-164. LI Jun-ling, LU Fang-yun, FU Hua, et al. Research on the dynamic behavior of a PBX explosive[J]. Chinese Journal of High Pressure Physics, 2011, 25(2): 159-164. DOI:10.11858/gywlxb.2011.02.012 |
[9] |
蓝林钢, 温茂萍, 李明, 等. PBX一维动态拉伸性能实验研究[J].
含能材料, 2011, 19(6): 684-688. LAN Lin-gang, WEN Mao-ping, LI Ming, et al. One dimension dynamic tensile properties of PBX[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2011, 19(6): 684-688. |
[10] |
赵玉刚, 傅华, 李俊玲, 等. 三种PBX炸药的动态拉伸力学性能[J].
含能材料, 2011, 19(6): 194-199. ZHAO Yu-gang, FU Hua, LI Jun-ling, et al. Dynamic tensile mechanical propert ies of three types of PBX[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2011, 19(6): 194-199. |
[11] |
唐维, 李明, 温茂萍, 等. 四种强度准则在高聚物粘结炸药强度分析中的适应性[J].
固体力学学报, 2013, 34(6): 56-61. TANG Wei, LI Ming, WEN Mao-ping, et al. Adaptability of four strength criterions in polymer bonded explosives strength analysis[J]. Chinese Journal of Solid Mechanics, 2013, 34(6): 56-61. |
[12] |
混凝土结构设计规范国标编写组. 混凝土结构设计规范[S]. 中华人民共和国国家标准GB 50010-2010, 北京: 中国建筑工业出版社, 2010.
|
Strength property of TATB based polymer bonded explosive(PBX) under triaxial tension stress state was studied by an end-restraint cylinder tension collapse method. The stress field under critical load was quantitatively analyzed through finite element simulation according to the tested boundary conditions. The adaptability of three typical criterions(uniaxial-strength,Mohr-Coulomb,twin-shear) in strength analysis was illustrated using principal stresses of initial failure position.