力学学报, 2021, 53(6): 1795-1806 DOI: 10.6052/0459-1879-21-017

生物、工程及交叉力学

TA2钛合金开口柱壳外爆碎片分布研究1)

吴文苍, 董新龙,2), 庞振, 周风华

*宁波大学冲击与安全工程教育部重点实验室, 浙江宁波 315211

宁波大学机械工程与力学学院, 浙江宁波 315211

STUDY ON FRAGMENTS DISTRIBUTION OF EXPLOSIVELY DRIVEN CYLINDERS FOR TA2 TITANIUM ALLOY1)

Wu Wencang, Dong Xinlong,2), Pang Zhen, Zhou Fenghua

*Key Laboratory of Impact and Safety Engineering (Ningbo University), Ministry of Education, Ningbo 315211, Zhejiang, China

Faculty of Mechanical Engineering and Mechanics, Ningbo University, Ningbo 315211, Zhejiang, China

通讯作者: 2)董新龙, 教授, 主要研究方向: 冲击动力学. E-mail:dongxinlong@nbu.edu.cn

收稿日期: 2020-01-5   接受日期: 2021-03-25   网络出版日期: 2021-06-18

基金资助: 1)国家自然科学基金资助项目.  11672143
国家自然科学基金资助项目.  11932018

Received: 2020-01-5   Accepted: 2021-03-25   Online: 2021-06-18

作者简介 About authors

摘要

金属柱壳爆炸膨胀断裂机制及其对碎片分布、特征尺寸的影响是应用物理、力学、兵器工程等领域共同关心的重要课题, 但目前除数值模拟外, 考虑断裂机制的简单二维碎裂模型尚未出现.开展TA2钛合金开口柱壳在不同装药条件下的碎裂实验研究, 通过对软回收碎片的统计及微观分析, 探讨金属柱壳外爆断裂模式及二维碎片分布规律, 结果显示:(1) TA2钛合金柱壳在实验爆压(7 $\sim$ 25 GPa)下宏观断口均为剪切断裂模式, 但机制不同, 在较高爆压下柱壳剪切断裂由多重绝热剪切带破坏控制, 在较低压力下为剪切破坏;(2) 与一维拉伸碎裂相比, 柱壳爆炸碎裂不充分, 碎片质量更符合$\beta=1$ (或更接近1)的指数分布; 爆炸碎裂越充分, 碎片越小并趋于均匀, $\beta$趋于较小的值, 趋向Mott和Linfoot提出的泊松统计分布形式;(3) Rayleigh分布可以较好描述柱壳碎片的宽度分布规律, 不同爆压下柱壳碎片宽度归一化尺寸分布具有相似性, 呈现"量子化"特性, 即存在最小的特征尺寸; (4) TA2柱壳碎片特征尺寸远大于G-K剪切断裂公式预测的尺寸, G-K剪切式描述的是多重绝热剪切带间距.本研究为金属柱壳碎片特征、分布规律及其模型分析提供了重要参考.

关键词: 金属柱壳 ; 碎裂 ; 断裂模式 ; 绝热剪切 ; 特征尺寸

Abstract

Understanding the characteristics of fracture mode, fragment distribution and its controlling factor for explosively driven metal cylinders are important topics in applied physics, mechanics, weapon engineering and other fields. However, except for numerical simulation, a simple two-dimensional fragmentation model considering fracture mechanism has not put forward. In this paper, the fragmentation of TA2 titanium cylinders with varying charge were carried out experimentally. Through the analysis of macroscopic fracture and microscopic metallographic for the recovered fragments, the fracture mode and its fragmentation distribution of TA2 metal are discussed. The results show that:(1) The fragmentation of TA2 titanium alloy cylinders is all shear fracture modes under different detonation pressure of 7 $\sim$ 25 GPa, but the mechanism is different, which the fracture controlled by multiple adiabatic shear bands at higher explosion pressure; (2) different from the one-dimensional tensile fragmentation, the fragment mass distribution is modeled as an exponential distribution with $\beta = 1$ due to the insufficient fragmentation of metal cylinders. The $\beta $ can tend to be smaller when the fragmentation of explosively driven cylinders is more sufficient, which would be closer to the Poisson statistical distribution proposed by Mott and Linfoot; (3) The width distribution of cylinder fragments can be well described by the Rayleigh distribution. The normalized width distribution of fragments under different explosion pressures is similar and presents the characteristic of "quantization", which means a minimum characteristic size exists; (4) The characteristic sizes of TA2 shell fragments are much larger than that predicted by the G-K formula, as it is used to describe the distances between multiple adiabatic shear bands. Many factors affect the explosion of cylindrical shell, and the formation of each fragment is random, but the statistical regularity shows that the explosion process and fragment distribution are still predictable overall. The research provides an important reference for fragmentation characteristics, distribution, and model analysis of metal cylinders.

Keywords: metal cylinder ; fragmentation ; fracture mode ; adiabatic shear ; characteristic size

PDF (10885KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

吴文苍, 董新龙, 庞振, 周风华. TA2钛合金开口柱壳外爆碎片分布研究1). 力学学报, 2021, 53(6): 1795-1806 DOI:10.6052/0459-1879-21-017

Wu Wencang, Dong Xinlong, Pang Zhen, Zhou Fenghua. STUDY ON FRAGMENTS DISTRIBUTION OF EXPLOSIVELY DRIVEN CYLINDERS FOR TA2 TITANIUM ALLOY1). Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(6): 1795-1806 DOI:10.6052/0459-1879-21-017

引言

金属材料在爆炸、冲击等高速变形中发生的碎裂现象, 一直是冲击动力学与材料领域的研究热点[1-26].最早在二战时期, Gurney[1], Taylor[2], Mott[3-4]等针对炮弹、手榴弹等壳体装药武器在爆炸载荷下的碎裂问题开展分析研究.其中, Mott[3-4]在研究炸弹壳体的动态碎裂时, 将二维柱壳碎裂简化为金属环的拉伸断裂问题, 并假设断裂的发生具有随机性与瞬时性, 提出断裂卸载稀疏波控制碎片特征尺寸的一维统计分析方法, Mott卸载波思想为后来众多研究奠定了学术基础. 随后, 一些学者从不同理论出发对金属、岩石、玻璃等材料的碎裂现象进行解释, 各自提出了碎片尺寸分布的计算公式, 并且做了大量实验加以验证[5-7].尤其, Kipp等[8-9]在Mott卸载波的传播距离控制碎片尺度的思想基础上, 引入内聚断裂模型, 给出了Mott模型中的特征时间, 通过Mott波传播距离来描述碎裂过程中产生碎片的平均尺寸, 其公式形式简单, 不需要做关于碎裂过程的随机假设, 便于应用.陈磊等[10-11]的一维数值模拟工作表明:Kipp和Grady[8-9]的韧性碎片尺寸公式(G-K公式)能较好地预测碎片平均尺寸的下限, 其中碎片个数和应变率2/3次幂关系的趋势与Grady[12]、Altynova等[13], Zhang等[14]学者的大量一维韧性碎裂实验结果基本一致.Mott卸载波传播距离控制碎片尺度的思想在一维应力理论框架上很好地描述了固体在冲击拉伸作用下的碎裂过程.

事实上, 从军事应用需求角度来说, 二维碎裂问题比一维应力膨胀环碎裂问题更具有工程意义. 但目前除数值模拟外, 简单的包含裂纹传播过程的二维碎裂理论模型尚未出现.对于在二维尺度上均匀膨胀的固体, Grady[12]曾定性地将二维碎裂问题简单地分解为两个一维情况的正交, 将二维应力状态下的卸载波简单处理为两个方向上相互独立的一维Mott波, 得到了二维均匀拉伸载荷下的平板碎片的面积尺度, 但碎片尺度公式显著高估了双轴拉伸碎片面积. 在实验分析方面, Grady和Hightower[15]通过爆炸驱动4410钢厚壁圆筒实验, 发现对于柱壳爆炸膨胀碎裂, 不是单纯的拉伸破坏, 而出现了拉伸$+$剪切的断裂模式.Grady等[16]建议可将断裂能定义为拉伸与剪切两部分的贡献, 并假设金属柱壳的碎片尺度与一维拉伸尺寸相同, 但未给出两者的比例及进一步的分析. Hoggart和Recht[17]实验观察到低爆轰压力下, 拉伸断裂为主要破坏模式; 而在高爆轰压力下, 碎片中更多的是剪切失效, 他们认为爆轰压力控制了断裂模式的转变, 剪切断裂占断口比例与爆炸压力有关.显然, 对于柱壳的碎裂过程及碎片分布, 其影响因素要比一维应力复杂, 而且还需要考虑断裂模式的影响.

本文对爆炸作用下TA2钛合金开口柱壳破坏模式及碎裂分布开展研究, 探究不同爆炸压力(应变率)下, 韧性金属的碎裂特征及其对碎片分布的影响, 比较G-K理论预测公式对于柱壳剪切断裂模式的适用性, 探讨剪切断裂能的意义, 以期为相关工程分析提供参考.

1 实验设计及碎片回收

1.1 试样材料及实验方法

实验柱壳采用工业纯钛TA2, 其密度为}4.51 g/cm$^{3}$, 成分组成比为Ti占99.6%, C占0.021%, Fe占0.06%, N占0.025%, H占0.004%, O占0.11%; 退火处理后原始组织为$\alpha $等轴晶.设计加工成内径19 mm、 壁厚3 mm、 长度60 mm的开口柱壳试样.

柱壳外爆实验原理如图1(a), 柱壳试样内填充不同的炸药, 一端设置传爆药及雷管, 如图1(b)所示. 实验中, 通过对柱壳内部填充不同密度、不同混合比炸药的方式, 来改变实验爆炸压力, 获得不同压力特征的载荷.

图1

图1   柱壳外爆实验原理及设计

Fig.1   Experimental assembly and site of cylindrical shell


爆炸实验在一个直径为800 mm的容器中完成, 实验中, 将柱壳密封并置于充水的塑料袋中, 对碎片进行软回收, 减小碎片的二次损伤, 如图1(c)所示.

1.2 炸药种类及爆轰特性

实验中, 在柱壳内部实心填充不同密度、不同混合比炸药, 以获得不同的爆炸压力.炸药分别为:密度为0.86 g/cm$^{3}$的黑索金粉、密度为1.48 g/cm$^{3}$的黑索金$+$添加剂的药柱以及密度为1.65 g/cm$^{3}$的B炸药药柱.采用探针法分别测试了三种炸药的爆速$D$, 虽然炸药爆炸过程相当复杂[27], 但此处可以忽略不同炸药的化学性质及爆炸产物的影响, 采用简化公式[28]对炸药的爆炸C-J压力$P_C\text{-J}$进行估算, 如表1所示.

表1   不同装药条件下的爆速及爆压

Table 1  Detonation velocity and detonation pressure under different charge conditions

新窗口打开| 下载CSV


2 实验结果及分析

2.1 碎片特征

2.1.1 宏观断裂特征

图2所示为不同爆压实验下回收得到的典型碎片, 主要碎片呈长条状, 长度与原始柱壳长度几乎相同, 少量在轴向上发生断裂.随爆炸压力增大, 回收到的碎片数量增加, 宽度更趋于均匀.

图2

图2   不同爆压下的回收碎片

Fig.2   Recovered fragments under different detonation pressures


对实验回收的碎片宏观观察显示: 柱壳均为剪切断裂破坏, 碎片横截面存在五种典型的特征, 如图3所示, 未见拉伸及拉剪混合断裂模式.绝大部分碎片呈图3所示的a, b, c, d及e 五种特征, 其中, d, e碎片都是轴向长度较短的碎片, 基本是a, b, c 三种主要条状碎片交界处剥落产生的小碎片. 值得注意的是:实验中对每块碎片逐一观测, 没有发现碎片中有与a形状相反的内壁短、外壁长的倒梯形状碎片, 这可能与柱壳断裂起始及Mott卸载波传播相关.

图3

图3   典型碎片截面示意图

Fig.3   Typical cross-section of fragments


进一步观察发现: (1) 在较低爆压(7 GPa)下, 回收碎片主要呈长条状, 长度基本与原始柱壳的长度相同, 碎片的断口方向与径向近似成45$^\circ$或135$^\circ$, 碎片截面主要呈a, b, c 三种类型;个别宽度较大的碎片内、外表面存在较多未断开的轴向裂纹, 这是由于裂纹沿轴向并未发展完全, 最终没有形成独立存在的新碎片, 如图4(a)所示. (2) 当爆压为22 GPa时, 除长条状碎片, 沿轴向断裂的小碎片数量增多. 另外, 碎片内表面存在沿环向分布台阶及沿轴向的微裂纹, 如图4(b)所示. (3) 当爆压为25 GPa时, 碎片尺寸较为均匀, 整体形状更加规则, 断口平滑光亮, 呈发蓝、发紫颜色, 但外表面存在一些沿轴向发展的变形集中带或微裂纹如图4(c)所示.特别值得注意的是, 碎片截面的剪切断口方向均为同一旋向, 即图3中碎片b形式.胡海波等[29]曾在多种金属爆炸碎片中观察到类似的单旋破坏现象, 但到目前为止其机制并不清楚.

图4

图4   不同爆压下的典型碎片

Fig.4   Typical fragments under different detonation pressures


2.1.2 微观断裂机制

对三种爆压下典型碎片横截面进行微观金相分析, 如图5所示. 微观金相显示: 在较低爆压(7 GPa)下, 碎片截面较为平整, 虽然在内、外表面局部存在裂纹, 前端可见绝热剪切局域化带特征, 但碎裂断口不是由绝热剪切带破坏导致的, 如图5(a)所示. 在22 GPa爆压下, 断口为与径向成45$^{\circ}$的剪切断裂, 但与7 GPa爆炸压力作用下不同, 在试样内表面分布着一组沿剪切方向(与径向45$^\circ$ 与135$^\circ$)交叉、平行的多重绝热剪切带, 多重剪切带间相对滑移, 一些裂纹已沿剪切带扩展, 在碎片内表面形成裂纹或台阶, 与宏观碎片内表面所见的沟槽及微裂纹相对应.断口是裂纹沿剪切带扩展至外表面形成的绝热剪切断裂模式, 如图5(b). 在25 GPa爆压下, 典型碎片截面成近似平行四边形, 内表面存在一组与径向成45$^\circ$、平行的多重绝热带及沿剪切带发展的裂纹, 但趋于单旋的绝热剪切断裂模式, 如图5 (c)所示.

图5

图5   不同爆压下的碎片截面金相

Fig.5   Metallography of cross section of debris


对不同爆炸压力下的柱壳碎片截面的宏观及微观分析可见:虽然TA2钛合金柱壳宏观断口均呈现剪切型断裂, 但微观机制不同.在较高爆压(22, 25 GPa)下, 为绝热剪切断裂机制控制; 而在较低压力(7 GPa)下, 为非绝热剪切断裂. 另外, 对于绝热剪切带控制的断裂, 试样内表面呈现的台阶分布与多重绝热剪切带滑移位置相对应, 剪切带之间的相对错动, 在碎片内表面形成沿轴向分布的沟槽及微裂纹.

2.2 碎片统计方法及分布规律

2.2.1 统计方法

回收碎片特征显示: 不同爆炸压力下的碎片主要呈长条状, 与试样长度相近.小碎片主要是由于局部剥落, 或不同截面剪切裂纹发展交错造成断裂. 因此, 柱壳爆炸膨胀断裂过程可近似假设为定常爆轰作用下的平面应变问题, 采用统计碎片宽度方法来表征碎片的特征尺寸. 由于剪切裂纹沿轴向呈台阶分布, 不同截面尺寸存在较大差别, 直接多点测量碎片平均宽度误差较大. 另外, 截面存在图3所示的多种形状, 难以统一定义宽度.沿碎片轴向多点厚度测量发现: 同一试验回收的碎片沿轴向厚度较为均匀, 为此, 采用称重及测量碎片平均厚度、长度来统计碎片. 进一步假设材料变形体积不变, 计算碎片的名义宽度$w$

$w_{i} =m_{i}/(\rho l_{i} h_{i})$

式中, $w_{i}$为碎片宽度, $m_{i}$为碎片质量, $\rho $为材料密度, $l_{i}$为碎片长度, $h_{i}$为平均厚度, 下标$i$代表第$i$块碎片.

2.2.2 碎片质量分布规律

统计时, 采用天平对每块碎片进行称重, 并同时测量碎片的长度及上中下3个位置厚度.表2为不同压力下典型试样碎片平均厚度及碎片的回收率$\eta$. 可见随爆炸压力增大, 碎片厚度减小.

表2   不同爆炸压力下回收碎片统计

Table 2  Statistics of the recovery fragments

新窗口打开| 下载CSV


图6给出了不同爆炸压力下回收碎片数量分布, 可以看到:在较低爆压下的碎片质量分布区间范围更大.

图6

图6   碎片质量分布

Fig.6   Mass distribution of fragmentation


由于柱壳断裂起始位置存在随机性, Mott和Linfoot[30]在对炸弹碎片特性分析中, 采用统计方法建立统一尺度比拟模型, 以便评估参数变化对弹头碎片特性影响, 并提出了碎片质量的泊松统计分布关系.Grady和Kipp[31]考虑到炸药、装药特征和外壳金属广泛性, 提出采用具有更广适应性的形式

$N(>m)=N_{0}e^{-(m/\mu)^{\beta }}$

其中, $N$ ($>m$)为质量大于$m$的碎片数, $N_{0}$为碎片总数, $m$为碎片质量, $\mu $为碎片质量参数, $\beta $为分布形状参数.这一关系事实上包含了Mott模型的特例, $\beta=1/2$即对应于Mott和Linfoot[30]提出的泊松统计分布形式.图6中同时给出了按式(2)中$\beta =1$时拟合的曲线, 可见取$\beta =1$与实验碎片质量分布趋势符合较好, 但对于大碎片的分布描述并不够准确. 这主要有两个原因:一方面是由于较多的小质量碎片未完全被回收; 另一方面, 与一维应力拉伸断裂不同, 二维柱壳断裂发展不充分, 如图4图5所示, 大尺寸碎片表面存在系列轴向裂纹, 如果充分碎裂, 碎片宽度应更小.

在碎片质量分布实验研究方面, Mock和Holt[32]曾对工业用纯铁和两种不同热处理的厚壁HF-1钢管(径厚比为3 $\sim$ 4)爆炸碎裂碎片分布进行了统计, 指出$\beta=1/2$的指数模型可以很好描述碎片质量分布规律. 另外, 在一维应力碎裂中, 很多学者也都采用$\beta =1/2$的分布.而Grady与Kipp[31]针对金属球壳动态破碎研究中, 认为采用$\beta =1$的模型能更好描述了碎片的分布, 并且明确认为$\mu$为特征尺寸碎片的质量.Odintsov[33]对Mott形式的碎片质量分布规律做过更详细的研究, 也认识到不同大小碎片段分布的断裂机制可能不同, 规律也可能不同.

事实上, 碎片质量分布的规律性与Mott卸载波控制相关. 分析认为, 对于一定的装药系统及爆炸碎裂过程, 存在随机的不充分碎裂, 导致实际碎片尺寸在统计上的不均匀性, 使分布变宽.按式(2)统计拟合时, $\beta $值趋于变大.如在相同柱壳壁厚下, 内径越大, 爆炸碎裂越充分, 碎片越小并趋于均匀, $\beta$的值应会变得较小. 与柱壳碎片相比, 一维应力动态膨胀环碎片分布更接近$\beta =1/2$.

虽然断裂存在随机性, 但在给定的加载率下, 碎片分布的规律性与Mott卸载波相关, 二维碎片分布应具有特征尺寸, 下面进一步就柱壳碎裂碎片尺寸分布特性进行分析.

2.2.3 碎片宽度分布及特征尺寸

采用碎片的环向宽度来表征碎片的尺寸, 图7给出按式(1)计算的TA2柱壳不同爆压下的碎片宽度分布, 可见尺寸分布存在相似性, 大多碎片尺寸集中带靠近较小尺寸区域, 爆炸压力越低碎片尺寸分布宽度越大. 采用Weibull积累分布函数拟合

$N\left( {>s} \right)=N_{0} \exp \left[ {-\left( {\frac{s-s_{\min} }{s_{0} }} \right)^{n}} \right]$

图7

图7   碎片宽度分布及概率拟合

Fig.7   Fragment width distribution and probability fitting


其中, $N$ ($>s$)表示碎片尺寸大于$s$的碎片个数, $N_{0}$为碎片总个数, $s_{\min}$为有效最小碎片尺寸, $n$为形状参数, $s_{0}$为标度参数, $s>s_{\min}$.周风华等[34]、郑宇轩等[35]对一维应力韧性金属圆环碎片长度尺寸研究认为:$n=2$时的Weibull函数, 即Rayleigh分布可以较好描述碎片宽度尺寸分布.对于平面应变柱壳碎片宽度分布拟合比较发现, Rayleigh同样可以较好地来表征柱壳碎片的尺寸分布, 如图7所示, 即满足

$n\left( s \right)= {\frac{2\left( {s-s_{\min} } \right)}{s_{0}^{2} }} \exp \left[ {-\left( {\frac{s-s_{\min} }{s_{0} }} \right)^{2}} \right]$

其中$s>s_{\min}$, 从图7可见: 虽然碎片数量较少, 但不同爆炸压力下平面应变柱壳碎片宽度分布仍与Rayleigh分布趋势吻合较好, 并且随爆炸压力增大, 碎片数量增加, 分布曲线趋势符合更好.

不同爆压下碎片分布规律结果显示: (a)Rayleigh分布形式可以描述柱壳碎片尺寸分布规律, 并随着加载率的提高, Rayleigh概率密度分布趋向变窄; (b) 随着爆压(加载率)的增大, 碎片特征尺寸呈减小的趋势. 另外, 爆压为25 GPa时的回收得到的碎片数量虽然增加, 但最小碎片尺寸、特征尺寸都较爆压为22 GPa的略大.这与爆压为25 GPa时裂纹大多为沿单旋绝热剪切带发展相关, 由于断裂面单旋发展, 这样交叉的裂纹少, 使得剥落的小碎片数量减小.

按Rayleigh分布可计算得到不同爆炸压力下的碎片特征尺寸

$s_ave =s_{\min} +\sqrt \pi s_{0} /2$

计算得到的爆压7 GPa, 22 GPa, 25 GPa下碎片的特征尺寸分别为5.67 mm, 3.87 mm, 4.59 mm.可见: 爆炸压力(应变率)增大, 碎片特征尺寸随应变率增大出现先减小, 随后又增大的现象. 按Mott卸载波传播距离控制碎片特征尺寸思想, 理论上, 应该随应变率提高, 碎片尺寸减小.但爆炸压力25 GPa时碎片特征宽度反而比22 GPa增大.这与25 GPa时断裂扩展模式改变有关, 在25 GPa时, 断裂裂纹呈单旋方向, 碎片截面为平行四边形, 使裂纹交叉产生的剥落小碎片少, 实际碎片分布跨度相对要小. 事实上, 从控制TA2钛合金在较高爆压下碎裂绝热剪切断裂机制看, 较高爆压25 GPa时的多重绝热剪切带间距还是比22 GPa时减小(如图4图5剪切带间距).

2.2.4 尺寸归一化

对不同爆压下碎片的特征尺寸对每块碎片宽度分布归一化处理, 即

$\frac{N(>s)}{N_{0} }=\exp \left[-(\frac{\bar{{s}}-\bar{{s}}_{\min} }{\bar{{s}}_{0} } )^{2} \right]$

式中取$\overline s_{\min} =s_{\min} /s_ave $, $\bar{{s}}_{0} =s_{0}/s_ave $. 得到TA2柱壳在不同爆压下的归一化累积分布$N$ ($>s$)/$N_{0}$和归一化尺寸$s$/$s_ave$分布, 如图8所示.可以看出: (1) 数据基本都处在同一条分布区域带内, 说明碎片归一化分布具有相似性, 可以采用同一条归一化Rayleigh曲线来描述; (2)在不同爆压下碎片尺寸分布均出现"台阶状"分布特征, 表明碎片的分布并不是连续的, 而是集中在某些特殊宽度附近, 即所谓的量子化特性.而且在高压25 GPa (高应变率)时, 由于单旋断裂, 碎片特征更能反映Mott波控制规律, 故量子化特征更明显. 表明: 虽然柱壳碎裂的裂纹起始存在随机性, 但碎片宽度具有自组织特征, 碎片宽度为某一"特征尺寸"的倍数, 多源损伤破坏碎裂尺寸是受材料内在特性控制的.

图8

图8   不同爆压下归一化尺寸累积分布

Fig.8   Cumulative distribution of normalized sizes under different detonation pressures


近似地, 对图8不同爆压下碎片归一化尺寸"台阶状"数据分布集中带取平均值$\overline s_{i} $, 计算其相邻台阶间距并取最小值$\{\min(\overline s_{i+1} -\overline {s}_{i})$, $i=1, 2, \cdots , n\}$, 计算得到7 GPa, 22 GPa和25 GPa下"台阶状"碎片尺寸最小间距分别为1.71 mm, 0.33 mm和0.42 mm.对于碎片的自组织"特征尺寸", Mott提出是由拉伸断裂过程产生的卸载波(Mott波)传播距离控制, Grady-Kipp基于材料拉伸断裂软化演化特征时间的一维应力韧性碎裂模型抓住了韧性碎裂过程中的主要矛盾, 周风华等[10-11]的数值模拟工作表明Grady-Kipp模型所预测的韧性碎片特征尺寸公式较好地预测碎片的平均尺寸的下限, 与Grady[12], Altynova等[13], Zhang等[14]学者的一维韧性碎裂实验结果也一致, 下文对剪切特征尺寸预测及影响因素进行讨论.

3 关于二维碎裂特征尺寸

3.1 剪切碎裂G-K公式

对于柱壳拉伸碎裂破坏, Mott[4]提出了卸载波控制碎片特征尺寸的分析思想, 给出了一维应力膨胀环碎片的特征尺寸分析方法.Grady等[36]、Kipp等[8]进一步扩展了Mott思想, 采用内聚断裂模型描述断裂能量耗散和应力释放率, 给出了Mott模型中的特征时间, 无需引入碎裂起始随机假设, 直接得到了Mott波传播距离, 确定碎片的特征尺寸

$s=\left( {\frac{12W}{\rho {{{\dot{{\varepsilon }}}}}^{2}}} \right)^{1/3}$

式中, $\rho $为金属密度, $\dot{{\varepsilon }}$为膨胀应变率, $W$为材料拉伸断裂能.

Grady等[16, 37]进一步又对拉伸碎裂模型进行修改, 来解决二维剪切局域化自组织问题, 认为在绝热剪切情况下, 式(7)仍适用, 只要将$W$变为绝热剪切带断裂时单位面积耗散的塑性功, 则绝热剪切带的自组织间距为

$s=\left( {\frac{12W}{\rho {{{{\dot{{\gamma }}}}}}^{2}}} \right)^{1/3}$

其中, ${{{{\dot{{\gamma }}}}}}$为剪切应变率.

进一步, 如假定剪切断裂时的应力释放由热软化引起, 假设有以下形式

$\tau =\tau_{Y} [1+\alpha (T_{o} -T)]$

式中, $\tau_{Y}$为流动应力, $\alpha $为热软化系数, $T-T_{0}$为绝热温升.简单取$\alpha =1/\Delta T$, 即线性热软化, 考虑剪切带内绝热温升, 得到剪切带形成的最小局域化时间, 导出多重剪切带间距$s_{0}$为

$s_{o} =2(\frac{9c^{2}\rho \chi }{\tau_{Y} \alpha^{2}{{\dot{{\gamma }}}}^{3}} )^{1/4}$

其中, $c$为材料比热, $\chi $为热扩散系数.由式(8)、式(10)得到了在绝热剪切局域化过程中的能量耗散

$W_{0} =\frac{2\rho c}{\alpha }\left( {\frac{9\rho ^{3}c^{2}\chi^{3}}{\tau_{Y}^{3} \alpha^{2}{{\dot{{\gamma }}}}}} \right)^{1/4}$

对于金属圆筒的碎片的特征尺寸, Grady认为同样可由界面断裂能$W$表征, 形式与一维应力碎裂情况下推导所得式(7)相同.只是当断裂面的形成由拉伸断裂和剪切断裂两种失效模式共同主导时, 断裂能需要考虑拉伸和剪切两种破坏在断裂面中所占比例, 以适当的方法加权计算.拉伸失效断裂能可采用金属断裂韧性$K_c$估算

$W=\frac{K_c^{2} }{2E}$

而剪切失效断裂能可采用式(11)的能量耗散计算.

3.2 特征尺寸及断裂能

3.2.1 碎片特征尺寸

实验显示: TA2钛合金柱壳外爆碎裂均为剪切断裂, 未见拉伸及拉伸$+$剪切混合断裂模式. 因此, 采用材料剪切失效断裂能$W$, 按G-K剪切式(8)分析碎片特征尺寸.设柱壳爆炸膨胀变形过程近似为平面应变状态, 其应变近似取为

$\varepsilon_{\theta \theta } =\left( {R(t)-R_{0} } \right)R_{0}$

式中, $R(t)$为$t$时刻的柱壳半径, $R_{0}$为初始柱壳半径.

开口柱壳断裂近似作为平面应变膨胀破坏, 则剪切应变$\gamma $为

$\gamma =\frac{2}{\sqrt 3 }\varepsilon_{\theta \theta }$

柱壳剪切应变率$\dot{\gamma}$近似为

$\dot{\gamma}=\frac{d\gamma }{dt}=\frac{2V_0}{\sqrt 3 R_{0} }$

其中, $V_{0}$为柱壳断裂时刻碎片速度.

Gurney[1]基于炸药能量平衡, 提出了一个估算壳体破片最终速度公式.Lloyd[38]考虑柱壳为有限长度的影响, 两端稀疏作用降低了壳体内炸药爆轰产物膨胀做功能力, 会使壳体破片最终速度$V_0$降低, 对Gurney速度公式进行修正

$V=\sqrt {\frac{2E}{(M/C+0.5)(1+R_{0} /L)}}$

式中, $M$/$C$是柱壳与炸药质量之比, $E$为格尼能通过炸药密度及爆速估算[39-40], $R_{0}$/$L$是柱壳的内径与长度比.Goto等[41]和Lambert等[42]应用该公式预测柱壳断裂时速度与测试结果均符合较好.

假设柱壳膨胀断裂后, 碎片的厚度不再变化, 则可采用回收碎片的平均厚度来估算柱壳断裂半径

$R^{\ast }=\frac{\Sigma \left( {l_{i} w_{i} h_{i} } \right)}{2\pi L_{0} h\cdot \eta }$

式中, $L_{0}$为柱壳高度, $h$为壁厚, $l_{i}$, $w_{i}$, $h_{i}$为同一碎片的长度、宽度、厚度, $\eta $为碎片质量回收率.

采用上述方法对柱壳断裂时速度及半径估算虽然与实际存在偏差, 但在同一条件下的实验中应具有一定的可比性.表3分别给出了不同炸药作用下按式(16)、式(17)、式(15)计算柱壳的膨胀速度$V_0$、断裂半径$R^*$及其应变率$\dot{\gamma}$的结果, 则按式(14)计算的断裂应变分别为0.37, 0.43, 0.45, 即随着爆炸压力的提高、应变率增大, 断裂应变增大.

表3   不同炸药作用下柱壳断裂参数估算

Table 3  Estimation of parameters of cylindrical shell under different explosives

新窗口打开| 下载CSV


对TA2纯钛的动态剪切性能及变形过程的绝热温升进行分析, 近似取流动应力$\tau_{Y}$为其屈服强度$Y=500$ MPa, 材料比热$c=500$ J/(kg$\cdot$K), 材料热软化系数$\alpha =1/\Delta T=2.5\times 10^{-3}$ K$^{-1}$[43].不考虑压力对热传导系数的影响, 则$\chi =\kappa /(\rho c)$, TA2钛等压热传导系数$\kappa =17$ W/(m$\cdot$K).图9给出了实验回收碎片统计分布分析得到的碎片特征尺寸及采用G-K剪切式(8)预测的碎片特征尺寸随应变率变化的比较, 其中, 剪切断裂能按式(11)计算. 可见:G-K式(8)预测的碎片特征尺寸与实际碎片相差很大.G-K公式预测的碎片特征尺寸(0.95 mm, 0.75 mm和0.72 mm)远小于回收碎片特征尺寸(5.67 mm, 3.87 mm和4.59 mm), 而更接近上节7 GPa, 22 GPa和25 GPa下归一化分析时近似得到的"量子化特征间距" 1.710 mm, 0.332 mm和0.417 mm.

图9

图9   碎片尺寸随应变率

Fig.9   Variation of fragment size with strain rate


实际上, G-K式(8)建议采用的剪切断裂能是基于绝热剪切带形成分析得到的最小能量释放率, 而实验回收碎片微观分析显示, 不同爆压下柱壳的碎片断裂模式并不相同, 在7 GPa压力下碎片断口虽为剪切断裂, 但并不是绝热剪切带破坏.图9中同时还给出22 GPa和25 GPa爆压下典型碎片微观分析测量得到的多重绝热剪切带的平均间距, 可见: G-K剪切式预测的碎片特征尺寸与多重绝热剪切带间距随应变率变化更为接近.

对于实验碎片特征尺寸与G-K预测值的差别, 主要包括以下几个方面原因: (1)G-K剪切断裂公式(8)建议的剪切断裂能是绝热剪切破坏的耗能, 即使对高爆压下柱壳绝热剪切碎裂, 其断裂并也不充分, 碎片中还包含多条平行、交叉绝热剪切带, 实际中并不是每条剪切带最终都形成裂纹或断口, 高压下量子化的"特征尺寸"反映的是多重绝热剪切带的间距, 碎片尺寸应为多重绝热剪切带特征间距的倍数; (2)许多裂纹沿柱壳轴向没有得到充分发展断开, 尤其对于二维问题, 不同截面断裂起始的随机性, 导致实际和预测模型差别较膨胀环碎裂相差更大; (3)对于低爆压, 由于不是绝热剪切断裂破坏, 采用绝热剪切断裂能分析预测, 造成两者的差别.

可见:Grady-Kipp剪切式(8)、式(11)是对多重剪切带间距的分析, 其结果与实际碎片尺寸相差两个数量级, Grady-Kipp给出的剪切断裂公式与实验多重绝热剪切带间距符合较好.

3.2.2 绝热剪切断裂能

Grady-Kipp的碎片特征尺寸分析公式依赖于断裂能释放率, 因此, 实际分析中明确断裂模式及其断裂能极为关键. 对于绝热剪切演化, 周刚毅等[44]曾采用平板帽型剪切试样对相同TA2钛的绝热剪切破坏演化特性开展实验研究, 测试了从剪切变形到绝热剪切破坏过程的等效应力-应变全曲线及相应剪切区的绝热温升, 如图10所示, 阴影面积$W$实际上即为应变率为3100 s$^{-1}$时的绝热剪切破坏的耗能, 近似为6 kN/m.

图10

图10   TA2绝热剪切破坏演化的热-力学响应[44]

Fig.10   Thermal-mechanical response of A TA2 adiabatic shear failure[44]


同时, 可以利用实验回收碎片的特征尺寸按G-K式(8)反演剪切断裂能, 按回收碎片中多重绝热剪切带平均间距按式(11)反演绝热剪切破坏断裂能.图11分别给出了按柱壳碎片尺寸、多重绝热剪切带平均间距反演的断裂能与实验测试得到的绝热剪切带破坏过程的断裂能的比较.可见: 绝热剪切断裂能明显低于剪切断裂能, 而G-K公式导出的绝热剪切带形成的最小耗能比帽型剪切实验得到的断裂能略小, 两者比较相近. 低爆压下的剪切断裂时, 材料剪切断裂能可能远大于绝热剪切断裂能.

图11

图11   TA2剪切断裂能分析

Fig.11   Estimation and comparison of TA2 shear fracture energy


可见:二维柱壳碎片的自相似特征尺寸不仅依赖于材料的应变率、断裂模式(断裂能不同), 另外还与柱壳断裂的充分程度相关.

4 结论

对TA2钛合金开口柱壳在不同装药条件下的碎裂特性开展实验研究, 通过对回收碎片的统计及典型碎片的宏微观分析, 探讨碎裂断裂模式及碎裂特性, 结果显示:

(1) TA2钛合金柱壳在实验爆压下宏观断口均为剪切断裂模式, 但不同压力下, 断裂机制不同. 在较低压力下(7 GPa), 为非绝热剪切带破坏, 而随压力增大至22 GPa, 25 GPa柱壳发生多重绝热剪切断裂, 并且在25 GPa时, 趋于单旋绝热剪切破坏现象.

(2) 随爆压增大(加载率)增大, 柱壳碎片质量、尺寸减小, 但25 GPa时的碎片特征宽度反而较22 GPa时略有增大, 这是由于25 GPa时, 多重绝热剪切带及裂纹趋于单旋发展, 碎裂时裂纹交叉少, 崩落产生的小碎片数量较少造成的, 而25 GPa时碎片内表面处多重绝热剪切带间距, 较22 GPa还是减小的.

(3) 柱壳爆炸膨胀碎裂的碎片质量分布符合$\beta =1$的指数分布, 分析认为对于一定的装药系统及爆炸碎裂过程, 由于存在随机的不充分碎裂, 使分布变宽, $\beta $值趋于变大. 如在相同柱壳壁厚下, 内径越大, 爆炸碎裂越充分, 碎片越小并趋于均匀, $\beta $值应会变得较小, 更趋于Mott和Linfoot[30]提出的泊松统计分布形式.

(4) Weibull分布中的Rayleigh形式可以较好描述柱壳碎片的宽度分布规律, 并且随着加载率的提高, Rayleigh概率密度分布逐渐变窄.不同应变率下的碎片尺寸累积分布归一化结果显示, 碎片分布具有相似性, 呈现"量子化"特性, 即碎片尺寸为特征尺寸的倍数.

(5) TA2柱壳碎片特征尺寸远大于G-K剪切断裂公式预测的特征尺寸, G-K剪切式与高爆炸压力下的多重绝热剪切带间距相近;G-K剪切断裂公式采用的是绝热剪切带形成的最小耗能分析, 不适用于低压下的非绝热剪切断裂情况.

参考文献

Gurney RW.

The initial velocities of fragments from bombs, shells and grenades

Ballistic Research Laboratories, 1943, 405: 1-22

[本文引用: 3]

Taylor GI. The Fragmentation of Tubular Bombs. London: Cambridge University Press, 1963: 387-390

[本文引用: 1]

Mott NF.

A Theory of the Fragmentation of Shells and Bombs//Grady D ed, Fragmentation of Rings and Shells

Shock Wave and High Pressure Phenomena, Berlin, Heidelberg: Springer, 2006: 243-294

[本文引用: 2]

Mott NF.

Fragmentation of shell cases

Proceedings of the Royal Society of London, 1947, 189(1018): 300-308

[本文引用: 3]

Grady DE.

Local inertial effects in dynamic fragmentation

Journal of Applied Physics, 1982, 53(1): 322-325

DOI      URL     [本文引用: 1]

Grady DE, Kipp ME.

Mechanisms of dynamic fragmentation: Factors governing fragment size

Mechanics of Materials, 1985, 4(3-4): 311-320

DOI      URL    

熊迅, 王珠, 郑宇轩 .

石英玻璃杆Taylor撞击实验的数值模拟

力学学报, 2019, 51(4): 1082-1090

DOI      [本文引用: 1]

采用离散元方法模拟石英玻璃杆Taylor撞击问题,再现了其破坏过程:在撞击端,杆以压缩失效波的形式破坏;在自由端,出现了密集的拉伸层裂破坏. 分析表明:层裂是失效波阵面应力快速下降引起的追赶卸载波,与弹性压缩前驱波在自由端反射引起的迎面卸载波相互作用的结果;随着撞击速度的增大,撞击端失效波造成的压缩破坏区域损伤程度增大,反射端层裂破坏损伤区域减小. 进一步对失效波阵面的结构变化及其波速问题进行了研究,发现失效区域随着扩张变成一段裂纹逐渐由密到稀的区段,将此区段分为高损伤区和低损伤区,研究发现由稀疏微裂纹组成的低损伤区的前端面传播速度和弹性前驱波速基本相同,为固定值;而高损伤区前端面的裂纹密度随着传播距离的增加变稀,直至过渡为低损伤区,其传播具有显著的速度衰减、端面模糊直至停止的过程. 高损伤前端面的平均速度随着撞击速度的增大而增大,并逐渐趋近于弹性波速. 最后与已有实验做了对比,发现实验中高速摄影观察到的玻璃中"失效波"阵面实际上是高损伤前端面,而稀疏的低损伤微裂纹很难捕捉.

(Xiong Xun, Wang Zhu, Zheng Yuxuan, et al.

Numerical simulations of taylor impact experiments of quartz glass bars

Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(4): 1082-1090 (in Chinese))

[本文引用: 1]

Kipp ME, Grady DE.

Dynamic fracture growth and interaction in one dimension

Journal of the Mechanics & Physics of Solids, 1985, 33(4): 399-415

[本文引用: 3]

Grady DE, Olsen ML.

A statistics and energy based theory of dynamic fragmentation

International Journal of Impact Engineering, 2003, 29(1-10): 293-306

DOI      URL     [本文引用: 2]

陈磊, 周风华, 汤铁钢 .

韧性金属圆环高速膨胀碎裂过程的有限元模拟

力学学报, 2011, 43(5): 861-869

DOI      [本文引用: 2]

用数值方法模拟了韧性金属圆环的自由膨胀碎裂过程. 采用Johnson-Cook热黏塑性本构模型描述材料的动态变形和热软化特性, 采用包含内聚力失稳断裂准则的Johnson-Cook型损伤断裂模型描述材料的破坏和分离过程, 采用结合单元消去技术的ABAQUS/Explicit程序进行分析. 在特定膨胀速度下对多个圆环进行碎裂数值实验, 获得碎片样本集合. 研究证实, Grady-Kipp基于塑性卸载波传播机制的韧性碎片尺寸公式可以较好地预测碎片的平均长度. 模拟再现了圆环碎裂过程中塑性卸载波的传播,揭示了Grady-Kipp公式合理性的物理基础.

(Chen Lei, Zhou Fenghua, Tang Tiegang, et al.

Finite element simulations of the high velocity expansion and fragmentation of ductile metallic rings

Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(5): 861-869 (in Chinese))

[本文引用: 2]

郑宇轩, 胡时胜, 周风华.

韧性材料的高应变率拉伸碎裂过程及材料参数影响

固体力学学报, 2012, 33(4): 358-369

[本文引用: 2]

(Zheng Yuxuan, Hu Shisheng, Zhou Fenghua.

High strain rate tensile fragmentation process of ductile materials and the effects of material parameters

Acta Mechanica Solida Sinica, 2012, 33(4): 358-369 (in Chinese))

[本文引用: 2]

Grady D.

Fragmentation of Rings and Shells

Berlin, Heidelberg: Springer, 2006: 153-195

[本文引用: 3]

Altynova M, Hu X, Daehn GS.

Increased ductility in high velocity electromagnetic ring expansion

Metallurgical and Materials Transactions A, 1996, 27(7): 1837-1844

DOI      URL     [本文引用: 2]

Zhang H, Ravi-Chandar K.

Dynamic fragmentation of ductile materials

Journal of Physics D Applied Physics, 2009, 42(21): 1-16

[本文引用: 2]

Grady DE, Hightower MM. Natural Fragmentation of Exploding Cylinders//Meyers M A, Murr L E, Staudhammer K P, eds. Shock Wave and High-Strain Rate Phenomena in Materials. San Diego, CA (USA), 1990.New York: Marcel Dekker, Inc., 1992: 713-721

[本文引用: 1]

Grady DE, Kipp ME.

The growth of unstable thermoplastic shear with application to steady-wave shock compression in solids

Journal of the Mechanics and Physics of Solids, 1987, 35: 95-119

DOI      URL     [本文引用: 2]

Hoggatt CR, Recht RF.

Fracture behavior of tubular bombs

Journal of Applied Physics, 1968, 39(3): 1856-1862

DOI      URL     [本文引用: 1]

高光发, 雷天刚, 戴兰宏.

第三届全国爆炸与冲击动力学青年学者学术研讨会报告综述

力学学报, 2020, 52(4): 1211-1219

(Gao Guangfa, Lei Tiangang, Dai Lanhong.

Review of the third national symposium on explosion and impact dynamics for young scholars

Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(4): 1211-1219 (in Chinese))

李涛, 胡海波, 尚海林 .

强约束球形装药反应裂纹传播和反应烈度表征实验

爆炸与冲击, 2020, 40(1): 18-25

(Li Tao, Hu Haibo, Shang Hailin, et al.

Propagation of reactive cracks and characterization of reaction violence in spherical charge under strong confinement

Explosion and Shock Waves, 2020, 40(1): 18-25 (in Chinese))

卢秋虹, 王宁, 范诚 .

壁厚对HR2钢柱壳爆轰加载下膨胀断裂行为的影响

材料研究学报, 2020(4): 241-246

(Lu Qiuhong, Wang Ning, Fan Cheng, et al.

Effect of shell thickness on expanding fracture behavior of HR2 steel cylinders under explosive loading

Chinese Journal of Materials Research, 2020(4): 241-246 (in Chinese))

杨云川, 朱建军, 郑宇 .

战斗部壳体爆炸破片体/线分形维数研究

兵工学报, 2018, 8: 1499-1506

(Yang Yunchuan, Zhu Jianjun, Zheng Yu, et al.

Research on the volume and line fractal dimensions of fragments from the explosion of warhead shell

Acta Armamentarii, 2018, 8: 1499-1506 (in Chinese))

罗渝松, 李伟兵, 陈志闯 .

内爆加载下金属柱壳的冻结回收方法

爆炸与冲击, 2020, 40(10): 87-96

(Luo Yusong, Li Weibing, Chen Zhichuang, et al.

A freezing recovery method for metallic cylinder shells under internal explosive loading

Explosion and Shock Waves, 2020, 40(10): 87-96 (in Chinese))

李茂, 侯海量, 朱锡 .

模拟破片杀伤战斗部空爆冲击波与高速破片群联合作用的等效试验方法

振动与冲击, 2020(1): 184-190

(Li Mao, Hou Hailiang, Zhu Xi, et al.

Equivalent test method to simulate combined damage action of air blast shock wave and high speed fragment group of fragment killing warhead

Journal of Vibration and Shock, 2020(1): 184-190 (in Chinese))

Dany F.

Estimating the metal acceleration ability of high explosives

Defence Technology, 2020, 1: 225-231

Abrosimov NA, Igumnov LA, Novosel'Tseva NA.

Numerical analysis of the effect of strain rate on the dynamic strength of cylindrical metal-plastic shells under explosive loading

Journal of Applied Mechanics and Technical Physics, 2020, 61(2): 267-276

DOI      URL    

Guo ZW, Huang GY, Liu CM, et al.

Velocity axial distribution of fragments from non-cylindrical symmetry explosive-filled casing

International Journal of Impact Engineering, 2018, 118: 1-10

DOI      URL     [本文引用: 1]

任会兰, 储著鑫, 栗建桥 .

B炸药爆炸过程中电磁辐射研究

力学学报, 2020, 52(4): 1199-1210

DOI      [本文引用: 1]

炸药爆炸过程中产生明显的电磁辐射干扰在很多领域引起了重视, 爆炸产生电磁辐射的机理和理论模型都不成熟, 实验仍然是研究这一现象的重要手段. 本文设计了不同质量的B炸药爆炸电磁辐射实验, 测量了爆炸整个过程中不同时刻的电磁辐射信号. 通过小波分析方法对实测信号进行处理, 得到真实信号和信号的时频特点, 小波分析得到真实信号的频谱主要分布在0$\sim$50 kHz. 利用自主开发的EXPLOSION-3D软件对实验工况进行数值模拟, 得到爆炸过程不同时刻流场的特点. 对比实验结果和数值模拟结果得出: 第一个脉冲信号是由B炸药爆轰产生的高温高压等离子体直接产生的电磁脉冲; 第二个脉冲信号是由地面反射的冲击波追上前沿冲击波, 空气冲击波阵面处形成的等离子体产生的电磁脉冲; 第三个脉冲信号是冲击波撞击测量线圈引起的无效信号. 其中, 第一个电磁脉冲幅值与当量的1/3次方基本呈线性关系, 其到达时间对炸药药量不敏感; 第二个电磁脉冲出现的时间与炸药的当量成指数关系, 药量越大, 出现时间越晚. 通过对第二个脉冲信号的分析, 得出了冲击波反射形成电磁脉冲信号时的爆炸波流场特征, 为后续相关理论机理研究提供了验证数据.

(Ren Huilan, Chu Zhuxin, Li Jianqiao, et al.

Research on electromagnetic radiation during the explosion progress of composition B explosives

Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(4): 1199-1210 (in Chinese))

[本文引用: 1]

л.п. 奥尔连科, 爆炸物理学(第3版上). 孙承纬, 译. 北京: 科学出版社, 2011: 170-174

[本文引用: 1]

(Oly. Orlenko, Explosion Physics. 3 (First). Sun Chengwei, Translate. Beijing: Science Press, 2011: 170-174(in Chinese))

[本文引用: 1]

胡海波, 汤铁钢, 胡八一 .

金属柱壳在爆炸加载断裂中的单旋现象

爆炸与冲击, 2004(2): 97-107

[本文引用: 1]

(Hu Haibo, Tang Tiegang, Hu Bayi, et al.

An study of uniform shear bands orientation selection tendency on explosively loaded cylindrical shells

Explosion and Shock Waves, 2004(2): 97-107 (in Chinese))

[本文引用: 1]

Mott NF, Linfoot EH.

A Theory of Fragmentation//Grady D ed, Fragmentation of Rings and Shells

Shock Wave and High Pressure Phenomena., Berlin, Heidelberg: Springer, 2006: 207-226

[本文引用: 3]

Grady DE, Kipp ME.

Fragmentation properties of metals

International Journal of Impact Engineering, 1997, 20(1-5): 293-308

DOI      URL     [本文引用: 2]

Mock W, Holt WH.

Fragmentation behavior of armco iron and HF-1 steel explosive filled cylinders

Journal of Applied Physics, 1983, 54: 2344-2351

DOI      URL     [本文引用: 1]

Odintsov VA.

Hyper exponential spectra of exponential fracture

Mechanics of Solids, 1992, 27(5): 42-48

[本文引用: 1]

Zhou F, Molinari JF, Ramesh KT.

Characteristic fragment size distributions in dynamic fragmentation

Applied Physics Letters, 2006, 88(26): 1210-192

[本文引用: 1]

郑宇轩, 陈磊, 胡时胜 .

韧性材料冲击拉伸碎裂中的碎片尺寸分布规律

力学学报, 2013, 45(4): 580-587

DOI      [本文引用: 1]

利用有限元方法模拟韧性金属圆环高速膨胀过程中的碎裂过程, 获得不同初始膨胀速度下碎片的样本集合. 通过对碎片的尺寸进行统计分析发现:(1)无论初始膨胀速度如何, 碎片的归一化尺寸分布具有相似性, 可以用一个具有初始阈值的Weibull分布描述, 近似地, 这个分布还可以简化为Rayleigh分布;(2)碎片尺寸的累积分布曲线呈现阶梯特性, 表现出较明显的"量子化"特性.在上述发现基础上, 建立一个Monte-Carlo模型:碎裂点来自于颈缩点, 颈缩之间的间距满足某种连续的Weibull分布, 而碎片的尺寸为随机的若干个颈缩间距之和.概率模拟表明:除非早期的颈缩间距分布很宽, 否则选择的离散性必然导致碎片尺寸分布呈现某种量子化特性.采用L04工业纯铝和无氧铜试件进行了爆炸膨胀碎裂实验, 回收得到的碎片尺寸分布结果与理论分析基本一致.

(Zheng Yuxuan, Chen Lei, Hu Shisheng, et al.

Characteristics of fragment size distribution of ductile materials fragentized under high strainrate tension

Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(4): 580-587 (in Chinese))

[本文引用: 1]

Grady DE, Benson DA.

Fragmentation of metal rings by electromagnetic loading

Experimental Mechanics, 1983, 23(4): 393-400

DOI      URL     [本文引用: 1]

Grady DE, Kipp ME, Benson DA. Mechanical Properties of Materials at High Rates of Strain, London, Institute of Physics, 1984: 315

[本文引用: 1]

Lloyd R.

Conventional warhead systems physics and engineering design

American Institute of Aeronautics and Astronautics, Inc, 1998: 25-27

[本文引用: 1]

Sun G, Wang X, Gao W, et al.

Expansion fracture behavior of metallic cylindrical shell caused by explosive detonation

IOP Conference Series: Earth and Environmental Science, 2019, 267(4): 1-8

[本文引用: 1]

沈飞, 王辉, 袁建飞 .

炸药格尼系数的一种简易估算法

火炸药学报, 2013, 36(6): 36-38

[本文引用: 1]

(Shen Fei, Wang Hui, Yuan Jianfei, et al.

A simple estimation method for gurney coefficient of explosive

Chinese Journal of Explosives & Propellants, 2013, 36(6): 36-38 (in Chinese))

[本文引用: 1]

Goto DM, Becker R, Orzechowski TJ, et al.

Investigation of the fracture and fragmentation of explosively driven rings and cylinders

International Journal of Impact Engineering, 2008, 35(12): 1547-1556

DOI      URL     [本文引用: 1]

Lambert DE, Weiderhold J, Hopson MV, et al.

Controlled loading fragmentation: Experiments and continuum damage modeling

AFRL-RW-EG-TR-2010-106, USA: Air Force Research Laboratory, 2010: 16-23

[本文引用: 1]

周刚毅.

TA2钛合金绝热剪切破坏特性及应力状态、晶粒度影响. [博士论文]

宁波: 宁波大学, 2018: 48-83

[本文引用: 1]

(Zhou Gangyi.

Adiabatic shearing behavior of TA2 titanium alloy and its influence of stress state, grain size. [PhD Thesis]

Ningbo: Ningbo University, 2018: 48-83 (in Chinese))

[本文引用: 1]

周刚毅, 董新龙, 付应乾 .

不同加载状态下TA2钛合金绝热剪切破坏响应特性

力学学报, 2016, 48(6): 1353-1361

DOI      [本文引用: 3]

一般认为绝热剪切现象在宏观上表现为材料动态本构失稳,即热软化大于应变硬化.本文采用帽型受迫剪切试样研究TA2钛合金的动态力学特性和本构失稳过程.首先对剪切区加载应力状态进行理论和数值分析,通过合理设计帽型试样,剪切区变形可近似按剪切状态处理;结合二维数字图像相关法(two-dimensional digitalimage correlation,DIC-2D)直接测试试样剪切区应变演化,给出帽型受迫剪切实验的等效应力-应变响应曲线.进一步,利用Hopkinson压杆对TA2钛合金开展动态压缩及帽型剪切对比试验研究,比较压缩、剪切试验得到的等效应力-应变曲线,采用“冻结”试样方法分析试样中绝热剪切局域化演化过程,探讨不同加载状态下TA2钛合金的绝热剪切破坏现象及其动态力学响应特性.实验结果表明,在塑性变形初始阶段,动态压缩及剪切加载下的等效应力-应变曲线符合较好,但随塑性损伤发展及绝热剪切带形成,两者出现分离,表明损伤及绝热剪切演化过程与应力状态相关.剪切试样实验得到的本构“软化”特性能够反映绝热剪切带起始、破坏演化过程的力学响应特性,而在动态压缩实验中,即使试样中已出现双锥形的绝热剪切带及局部裂纹分布,其表观等效应力-应变曲线并不出现软化特征,动态压缩实验无法得到关于绝热剪切起始、发展以及破坏的本构软化响应特性.

(Zhou Gangyi, Dong Xinlong, Fu Yingqian, et al.

An experimental study on adiabatic shear behavior of TA2 titanium alloy subject to different loading condition

Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(6): 1353-1361 (in Chinese))

[本文引用: 3]

/