外部爆炸载荷下表面粗糙度对45钢柱壳剪切带行为的影响

杨智程; 刘龙飞; 刘炼煌; 殷鹏志; 吴志强

doi:10.11858/gywlxb.20220506

外部爆炸载荷下表面粗糙度对45钢柱壳剪切带行为的影响

doi: 10.11858/gywlxb.20220506

湖南科技大学材料科学与工程学院, 湖南湘潭　411201

基金项目: 国家自然科学基金（11772127，51704112）

详细信息

作者简介:
杨智程（1996－），男，硕士，主要从事冲击动力学研究. E-mail：444356771@qq.com

通讯作者:
刘龙飞（1975－），男，博士，教授，主要从事冲击动力学研究. E-mail：lfliu1@hnust.cn

中图分类号: O347; O521.2
计量
- 文章访问数: 298
- HTML全文浏览量: 160
- PDF下载量: 26
出版历程
- 收稿日期: 2022-01-27
- 修回日期: 2022-03-02
- 网络出版日期: 2022-07-19
- 刊出日期: 2022-07-28

Effect of Surface Roughness on Shear Band Behavior of 45 Steel Cylindrical Shell under External Explosion Load

College of Materials Science and Engineering, Hunan University of Science and Technology, Xiangtan 411201, Hunan, China

摘要

摘要: 通过厚壁圆筒爆轰坍塌实验和有限元数值模拟，研究了45钢柱壳内壁表面粗糙度对其剪切带行为的影响。实验结果表明：外爆加载下，金属柱壳内壁表面粗糙度显著改变试样剪切带的成核位置和数量；当柱壳内壁表面粗糙度增大时，剪切带的数量和长度均增加，部分剪切带的扩展速度也增大；而当柱壳内壁表面的峰谷单元的平均宽度减小时，其剪切带的成核点数量增加，相邻剪切带之间的相互作用增强。有限元模拟和分析结果表明，柱壳中的最大剪切应力产生于柱壳内壁表面峰谷单元的谷底两侧，表面粗糙度的增大和峰谷单元平均宽度的减小均会提升柱壳内表面的剪切应力，促进柱壳中剪切带的成核和扩展，进而导致剪切带成核数量增多，主剪切带发展加快，剪切带屏蔽效应增强。
- 45钢 /
- 有限元法 /
- 剪切带 /
- 表面粗糙度 /
- 成核
Abstract: The effects of different inner wall surface roughness of 45 steel cylindrical shell on its shear band behavior were studied based on detonation collapse experiments of thick walled cylinder and finite element numerical simulation. The experimental results show that the surface roughness of the inner wall of the metal cylindrical shell significantly changes the nucleation position and the quantity of the shear bands under external explosion load. When the surface roughness of the inner wall of the cylindrical shell increases, the quantity and length of the shear bands increase, and the propagation speed of some shear bands increases as well. When the average width of the valleys of the inner wall decreases, the interaction between adjacent shear bands enhances with the increase of the nucleation points in the shear band. The finite element simulation and the theoretical analysis show that the maximum shear stress in the cylindrical shell generates on both sides of the valley bottom of the inner wall in the cylindrical shell. The shear stress on the inner surface of the cylindrical shell is enhanced with the increase of surface roughness, or with the decrease of average width of valleys. As the shear stress augments, the nucleation and the propagation of shear bands in the cylindrical shell are promoted, the nucleation quantity of the shear bands increase, the development of the main shear bands accelerate, and the shielding effect of the shear bands enhance.
- 45 steel /
- finite element method /
- shear band /
- surface roughness /
- nucleation

HTML全文

轴对称柱壳结构是战斗部件的典型代表，其在爆炸载荷下的变形和破坏问题一直是武器研究的热点之一^[1-10]。目前，对金属柱壳变形剪切的影响研究主要集中在不同材料的自身特性、显微组织和缺陷结构上^[11-17]，而表面加工介观状态对柱壳剪切带影响的研究相对较少。Nesterenko等^[18]首次通过厚壁圆筒（thick-walled-cylinder，TWC）装置研究了外爆下柱壳内剪切带的演化。而后，Nesterenko等^[19]和Meyers等^[20]先后对不锈钢和钛进行了TWC外爆实验，发现剪切带产生于柱壳内壁，并由内向外发展。黄西成^[21]经过理论分析得出，在外爆加载下柱壳的最大应力产生于柱壳的内表面。Yang等^[22]通过开展内表面预制缺口的TWC实验，发现内表面缺口能诱发剪切带优先形成，并对周围剪切带的形成产生屏蔽作用，表明内表面缺陷影响剪切带的形成和发展。任国武等^[23]通过钛合金柱壳内爆膨胀断裂实验研究发现，柱壳内壁率先诱发剪切带，剪切带由内向外贯穿，最终导致试样断裂。在爆炸载荷下，剪切带率先在柱壳内壁形核，柱壳内壁的表面加工介观状态会对剪切带产生重要影响，其中，内壁表面粗糙度如何影响剪切带行为以及影响规律如何仍有待探索。

随着计算技术的发展，有限元模拟在研究材料及结构的高速变形及断裂机理方面具有较明显的优势。与爆炸实验相比，有限元模拟有助于获得爆炸过程的详细数据和图像、理解其内在物理机制，常被用于爆炸冲击下材料及结构的变形、损伤和断裂等问题的研究^[24-33]。黄西成^[21]在外爆304不锈钢柱壳的数值模拟中发现，材料的初始微缺陷会诱发损伤，其随机分布会影响损伤演化，导致柱壳的非对称断裂。刘明涛等^[34-35]通过概率型材料本构实现了材料的不均匀性，由此模拟了金属柱壳在爆轰加载下的剪切断裂过程。Lovinger等^[36]在数值模拟中，将柱壳内壁正弦波边界作为表面缺陷，研究结果表明，只有当缺陷较大时，主剪切带才会影响周围的细小剪切带。以上研究主要通过调整材料参数来改变柱壳内表面的应力状态，并没有改变柱壳内表面材料的物理结构，柱壳内表面的粗糙度对其剪切行为影响的物理机制仍不清楚。45钢是一种优质的碳素结构钢，因调质处理后具有良好的综合力学性能而被广泛用于机械制造、交通运输和国防工业等领域，其在冲击载荷下的变形和断裂是人们非常关注的问题^{[12, 37-38]}。本研究通过制备具有不同表面粗糙度的45钢柱壳，采用TWC技术研究在爆炸载荷下45钢柱壳内表面粗糙度对其剪切带数量和轨迹的影响，并采用ANSYS/Workbench有限元软件对爆轰实验过程进行模拟，为理解柱壳在高应变率下剪切变形行为的物理图像和相关工程设计提供有价值的参考。

1. TWC实验及有限元模拟

1.1 样品及实验方法

本实验材料选择商业45钢棒材，其化学成分（质量分数）如表1所示。

表 1 45钢的化学成分

Table 1. Chemical composition of 45 steel %

C	Si	Mn	S	P	Ni	Cr	Cu	Fe
0.44	0.23	0.56	0.012	0.014	0.04	0.06	0.07	Rest

下载: 导出CSV

| 显示表格

将45钢棒材进行调质处理后，加工成内径为15 mm、外径为21 mm、长67 mm的圆管，将加工好的45钢圆管进行退火（550 ℃保温2 h），然后通过电火花线切割机采用不同的电参数对内表面进行切割，采用Marsurf M300C粗糙度仪测量内壁的粗糙度，结果如图1所示。其中：R_z表示轮廓中多个峰高和谷深的平均值之和，其值越大，表面轮廓中峰谷纵深越大；R_sm表示相邻轮廓峰谷单元的平均宽度，其值越小，表面轮廓中峰谷数量越多。图1(a)对应粗糙度较小的样品，R_z=5.8 μm，R_sm=455 μm，此时曲线整体较平整，峰高较小，峰谷间距较大；随着R_z逐渐增大到20.9 μm，R_sm依次减小到179 μm，曲线密集波动，峰高增大，峰谷间距缩小，如图1(e)所示。从图1(a)～图1(e)可以观察到内表面轮廓曲线大致呈现周期性变化，表面粗糙度越大的样品，曲线的峰高和谷深越大，峰谷数量越多，轮廓的峰谷单元平均宽度减小。

图 1 45钢样品的粗糙度轮廓线

Figure 1. Roughness contour line of 45 steel sample

下载: 全尺寸图片幻灯片

观察具有不同内表面粗糙度的45钢圆管横截面，如图2所示，可以看出：随着表面粗糙度的增大，表面轮廓峰谷的纵深差值不断增大，峰谷的宽度不断减小，峰谷的纵深差值和间距与图1中的轮廓曲线基本一致。

图 2 45钢圆管的表面形貌

Figure 2. Surface morphology of 45 steel cylindrical shells

下载: 全尺寸图片幻灯片

为了研究内表面粗糙度对45钢柱壳剪切带行为的影响，采用TWC实验技术对加工好的样品开展外爆加载实验。在圆管样品内外分别套上铜管作为制动管和驱动管，为减少应力波在管壁之间界面上反射，在圆管之间填充环氧树脂。外壳材料选用PVC管，在PVC管和驱动管及顶盖之间均匀填充炸药，炸药选取改性铵油炸药，炸药的密度为0.78 g/cm³，爆速为3.2 km/s。实验装置如图3所示，采用一端顶点起爆方式起爆，顶盖的斜坡设计是为了减少爆轰波对圆管轴向的加载，驱动管外炸药产生爆轰波，从而驱使铜管和45钢圆管向内坍塌发生塑性变形，诱发剪切带形成。实验装置的尺寸如表2所示。

表 2 实验装置的尺寸

Table 2. Size of experimental device

Component	Internal diameter/mm	Outer diameter/mm	Length/mm	Thickness/mm
Copper stopper	11	15	67	2
Specimen	15	21	67	3
Copper driver	21	25	67	2
Explosive	25	73		24

下载: 导出CSV

| 显示表格

图 3 TWC实验装置

Figure 3. Thick-walled cylinder experimental device

下载: 全尺寸图片幻灯片

为研究不同表面粗糙度对多重剪切带的自组织行为的影响，对爆炸后的样品进行回收。为了比较不同样品在相同应变下的剪切带行为，计算其等效应变ε_eff，即

${\varepsilon _{{\rm{eff}}}} = \frac{2}{{\sqrt 3 }}\ln {\frac{{{r_0}}}{{{r_{\text{f}}}}}}$

(1)

式中：r₀和 r_f分别为样品参考点的初始半径和最终半径。

对实验后的样品进行回收并切片，机械抛光后，用4%硝酸酒精溶液腐蚀，然后进行金相观察。通过测量实验后切片样品的内径，由式(1)计算等效应变，比较相同应变下不同样品的剪切带行为。

1.2 有限元数值模拟

在外爆加载下金属柱壳的坍塌过程中，由于径向爆轰速度比轴向高一个数量级，因此可以将滑移爆轰过程中圆管横截面近似为平面应变模型^{[12, 39]}，建立二维平面应变模型，见图4。模型中各部件材料尺寸均与实验一致，如表2所示。采用二维平面应变自适应网格，网格尺寸为30 μm，可以满足计算所需的精度要求，并且自适应网格能够较好地避免高应变率下材料变形计算过程中由网格破坏造成的计算不收敛导致的计算终止。部件之间的接触采用面-面无摩擦接触。

图 4 表面粗糙度的二维模型

Figure 4. Two dimensional model of surface roughness

下载: 全尺寸图片幻灯片

在实际的表面粗糙度加工中，峰谷之间的纵向高度差和峰谷数量是难以精确控制的。因此，利用ANSYS/Workbench模拟45钢柱壳时，根据实际情况中柱壳的内表面粗糙度轮廓，采用模型对其进行简化，通过模型函数中变量的变化来表征表面粗糙度中相关参数的差异，具体如下

$R = {R_0} + {R_\textit{z}}\cos \left( {\frac{{\text{π }}}{2}nt} \right)$

(2)

式中：R₀为柱壳初始内径；t=360°，为柱壳角度位置参量；n为峰谷数。2πR₀/n即为表征表面粗糙度的峰谷单元的平均宽度。

用式(2)表征内表面的粗糙度变化时，通过改变R_z和n的大小模拟不同内表面粗糙度，表3列出了模拟中选取的表面粗糙度参数。

表 3 表面粗糙度模拟参数

Table 3. Simulation parameters of surface roughness

No.	R_z/μm	R_sm/μm
1	10	471, 235, 157
2	50	471, 235, 157

下载: 导出CSV

| 显示表格

选用Johnson-Cook本构模型^[40]作为材料的本构关系，其表达式如下

${\sigma _{{\rm{eff}}}} = \left( {A + B\varepsilon _{\text{p}}^n} \right)\left( {1 + C\ln {{\dot \varepsilon }^*}} \right)\left( {1 - {T^{*m}}} \right)$

(3)

式中：σ_eff为等效应力；ε_p为等效塑性应变； ${\dot \varepsilon ^*}$ = $\dot \varepsilon$ / ${\dot \varepsilon _0}$ 为无量纲应变率， $\dot \varepsilon$ 为塑性应变率， ${\dot \varepsilon _0}$ 为参考应变率；T ^*=(T−T_r)/(T_m−T_r)为无量纲温度，T_m为材料熔点，T_r为环境参考温度，T为材料环境温度；A、B、C、m和n均为材料参数。

与之对应，采用Johnson-Cook失效模型，表达式为

${\varepsilon _{\rm{f}}} = \left[ {{D_1} + {D_2}\exp ({D_3}{\sigma ^*})} \right]\left( {1 + {D_4}\ln {{\dot \varepsilon }^*}} \right){(1+D_5}{T^*)}$

(4)

式中：ε_f为失效应变；D₁、D₂、D₃、D₄、D₅为材料参数；σ^*≈p/σ_eff，p为压力。

表4和表5列出了45钢和无氧铜的Johnson-Cook本构模型和失效模型参数。

表 4 材料的Johnson-Cook本构模型参数^[41]

Table 4. Johnson-Cook constitutive parameters of materials^[41]

Material	A/MPa	B/MPa	C	n	m
45 steel	507	320	0.064	0.28	1.06
Oxygen free copper	90	292	0.025	0.31	1.09

下载: 导出CSV

| 显示表格

表 5 材料的Johnson-Cook失效模型参数

Table 5. Johnson-Cook failure parameters of materials

Material	D₁	D₂	D₃	D₄	D₅
45 steel	0.10	0.76	1.57	0.005	−0.84
Oxygen free copper	0.45	3.44	−2.12	0.002	1.09

下载: 导出CSV

| 显示表格

模拟中炸药选用TNT炸药，采用JWL状态方程描述。表6列出了JWL状态方程参数A_TNT、B_TNT、R₁、R₂、ω，此外，还列出了密度ρ、爆速D、爆压p_CJ和格尼能E。

表 6 TNT的JWL状态方程参数^[42]

Table 6. Parameters of JWL equation of state of TNT^[42]

A_TNT/GPa	B_TNT/GPa	R₁	R₂	ω	E/(GJ·m⁻³)	p_CJ/GPa	ρ/(g·cm⁻³)	D/(km·s⁻¹)
371.2	3.231	4.51	0.95	0.3	7.0	21	1.63	6.93

下载: 导出CSV

| 显示表格

2. 实验及模拟结果

2.1 实验结果分析

图5为爆炸后等效应变ε_eff分别为0.73±0.01和1.22±0.01时样品的宏观图像。图5(a)表明，当柱壳的塑性变形较小时，内表面有少许“刺状”出现，预示着剪切带短小，将其称为剪切带成核阶段。随着应变增大，柱壳内表面变形明显，如图5(b)所示，内表面产生大量明显的“倒钩”，同时剪切带逐渐向外扩展，开始出现少量裂纹，称其为剪切带扩展阶段。通过金相观察，记录样品截面上所有剪切带的形貌、数量及位置，绘制成剪切带分布图，如图6所示。

图 5 不同应变下柱壳截面的宏观图像

Figure 5. Macro picture of cylinder cross section under different strains

下载: 全尺寸图片幻灯片

图 6 不同表面粗糙度样品的剪切带分布（ε_eff=0.73±0.01）

Figure 6. Distribution of shear band of samples with different roughnesses (ε_eff=0.73±0.01)

下载: 全尺寸图片幻灯片

图6(a)～图6(e)显示了等效应变为0.73±0.01，柱壳内表面R_z分别为5.8、9.7、13.8、17.3、20.9 μm时样品的剪切带分布，可以看出，剪切带数量分别为36、43、46、51、57，即剪切带数量随着粗糙度的增大而增加。结合图1、图2和图6，可以发现：粗糙度明显影响外爆载荷下柱壳中剪切带的成核位置；当R_sm减小时，表面轮廓总峰谷数量增加，在爆炸载荷作用下，柱壳的最大剪切应力产生于内壁，而内壁表面的峰谷会促使应力集中，使其成为剪切带潜在的成核点，从而诱发更多的剪切带形成。

图7(a)～图7(e)显示了等效应变为1.22±0.01，柱壳内表面R_z分别为5.8、9.7、13.8、17.3、20.9 μm时样品的剪切带分布，与之对应的剪切带数量分别为66、74、83、88、93。可以看出，该阶段的剪切带数量较其成核阶段（ε_eff=0.73±0.01）多，部分剪切带发展较快，其长度明显比周围剪切带长度大得多，长剪切带之间夹着多条短小的剪切带。当表面粗糙度增大时，柱壳中长剪切带的数量增加，短剪切带变得更短，剪切带长度的差异增大，屏蔽效应增强，并且部分剪切带出现单旋现象。分析该实验现象可知：当R_sm减小时，样品内表面轮廓的峰谷总数增多，致使剪切带的成核数量增多；柱壳内壁表面R_z的增大和R_sm的减小都会增大表面应力的集中程度，诱使某些剪切带率先成核并扩展，增强对周围剪切带的屏蔽效应，压制其周围剪切带发展，使其一直处于短小形态。因此，形成了剪切带长短差异增大和部分单旋的现象。

图 7 不同表面粗糙度样品的剪切带分布（ε_eff=1.22±0.01）

Figure 7. Distribution of shear band in samples with different roughnesses (ε_eff =1.22±0.01)

下载: 全尺寸图片幻灯片

2.2 模拟结果分析

为了进一步理解外爆加载下表面粗糙度对45钢柱壳剪切带行为影响的物理机制，对该过程进行了有限元模拟，以期提供有价值的参考。图8为R_z=10.0 μm，R_sm分别为471、235和157 μm时45钢柱壳在ε_eff =0.66时的等效塑性应变云图。图8(a)中，剪切带处在成核期，少量剪切带发展；图8(b)和图8(c)中，剪切带开始发展，随着R_sm的减小，剪切带数量增多。整体来看，图8中的剪切带都处于早期阶段，没有完全发展，剪切带只出现在样品的内壁，与图6所示的实验结果一致。

图 8 R_z=10.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳的等效塑性应变云图（ε_eff=0.66）

Figure 8. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm) of surface roughnesses at R_z=10.0 μm (ε_eff=0.66)

下载: 全尺寸图片幻灯片

图9为R_z=50.0 μm，R_sm分别为471、235、157 μm，ε_eff=0.66时45钢柱壳的等效塑性应变云图。当R_sm减小时，图9(a)、图9(b)和图9(c)中剪切带数量依次增加，并且图9(c)中长剪切带周围开始出现短小剪切带。另外，图9(a)中剪切带大多数呈现“两两相交”之势，而图9(b)和图9(c)中剪切带却没有出现该现象，部分剪切带朝着同一方向发展。随着R_sm不断地减小，剪切带的数量不断增加，同时短小的剪切带也增多，说明屏蔽效应增强。由图8和图9可知， R_z增大、R_sm减小，会促进剪切带的成核和扩展，剪切带数量增加，剪切带扩展的屏蔽效应增强。这与图6(c)、图6(d)和图6(e)所示的实验结果一致，部分长剪切带夹着短小剪切带，大部分短小剪切带方向一致，个别短小剪切带呈相交趋势。

图 9 R_z=50.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳等效塑性应变云图（ε_eff=0.66）

Figure 9. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm ) of surface roughnesses at R_z=50.0 μm (ε_eff=0.66)

下载: 全尺寸图片幻灯片

图10为ε_eff =1.15，R_z=10.0 μm，R_sm分别为471、235、157 μm时45钢柱壳的等效塑性应变云图。此时，剪切带处于扩展阶段，大部分接近柱壳外壁，少数贯穿柱壳壁。图10(a)中，剪切带整体呈蛛网状，相互交错；图10(b)中，开始出现少量剪切带方向一致的单旋现象；图10(c)中，单旋现象较图10(b)中更加明显。

图 10 R_z=10.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳的等效塑性应变云图（ε_eff=1.15）

Figure 10. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm ) of surface roughnesses at R_z=10.0 μm (ε_eff=1.15)

下载: 全尺寸图片幻灯片

图11为ε_eff=1.15，R_z=50.0 μm，R_sm分别为471、235、157 μm时45钢柱壳的等效塑性应变云图。3个样品中大部分剪切带都已贯穿柱壳壁。图11(a)中的剪切带“蛛网”比图10(a)明显，R_z的增大促进了剪切带的发展。图11(c)比图11(b)中的单旋现象更明显，这是因为随着R_sm的减小，剪切带之间的屏蔽效应增强，某个最大剪切应力方向（如45°方向）的剪切带形成后，屏蔽了另外一个最大剪切应力方向（如135°方向）的剪切带扩展，剪切带更易发展成单旋现象。综合图7、图10和图11可以看出，随着R_sm减小，剪切带成核数量增加。当R_sm较大时，成核点的间距较大，每个成核阶段的剪切带都能吸收能量进而发展演化，如模拟结果图10(a)、图11(a)和实验结果图7(a)所示，剪切带长度的差异不明显。随着R_sm减小，图10(b)、图10(c)、图11(b)、图11(c)与图7(b)～图7(e)的现象较为一致，剪切带成核数量增加导致剪切带间距减小，剪切带之间的相互影响增强，部分成核剪切带率先成核生长，对相邻的剪切带产生屏蔽作用，抑制了一些生长缓慢的剪切带，增强了形成单旋剪切带的趋势。

图 11 R_z=50.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳的等效塑性应变云图（ε_eff =1.15）

Figure 11. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm) of surface roughnesses at R_z=50.0 μm (ε_eff=1.15)

下载: 全尺寸图片幻灯片

图12(a)和图12(b)为具有相同R_sm、不同R_z的柱壳内表面的剪切应力分布云图。可以看出，最大剪切应力出现在峰谷单元中靠近谷底两侧45°和135°方向，且粗糙度大的样品，最大剪切应力大。具有相同R_sm、不同R_z的单元几何结构如图12(c)表示，可以将柱壳内表面粗糙度形成的谷底看作表面为半椭圆的裂纹，R_z的变化改变了半椭圆裂纹的切口深度。采用线弹性力学中的断裂力学方法，对切口尖端附近区域的应力场进行定性分析，可以得出该区域内任意一点剪切应力分量τ的表达式

图 12 表面粗糙度对柱壳内表面剪切应力的影响（R_sm=157 μm）

Figure 12. Effect of surface roughness on internal surface shear stress of cylindrical shell (R_sm=157 μm)

下载: 全尺寸图片幻灯片

$\tau {\text{ = }}\frac{{\sigma \sqrt {{\text{π}}{R_{{\textit{z}}}}} }}{{\sqrt {{{2\text{π} }}\gamma } }}{\text{cos}}\,\frac{\theta }{{\text{2}}}\,{\text{sin}}\,\frac{\theta }{{\text{2}}}\,{\text{cos}}\,\frac{{{\text{3}}\theta }}{{\text{2}}}$

(5)

式中：σ为距离尖端足够远处的均匀拉伸应力，θ为τ与径向的夹角。

图12(d)为45°方向柱壳内表面无量纲最大剪切应力随R_z的变化曲线。随着R_z的增加，无量纲最大剪切应力增大，更易诱发剪切带的成核与扩展。

图13(a)和图13(b)为具有相同R_z、不同R_sm柱壳的剪切应力分布云图。从图13(a)和图13(b)中可以看出，最大剪切应力同样出现在峰谷单元中靠近谷底两侧的45°和135°方向，且R_sm小的样品，最大剪切应力大。具有不同R_sm、相同R_z的单元几何结构如图13(c)所示，可以将表面粗糙度形成的谷底看作表面为半椭圆的裂纹，R_sm的变化则改变了半椭圆裂纹的切口宽度。同样地，基于断裂力学方法，切口尖端附近区域内任意一点的剪切应力分量可用式(5)描述。此时，在相同的R_z下，随着R_sm的减小，式(5)中的θ增大。为了考察R_sm对金属柱壳中剪切带行为的影响，基于式(5)获得了R_z=50.0 μm时无量纲剪切应力τ/σ随θ的变化曲线，如图13(d)所示。随着θ的增大，无量纲剪切应力增大，即无量纲剪切应力随着R_sm的减小而增大，较小的R_sm更易诱发剪切带的成核与扩展。

图 13 表面粗糙度峰谷单元平均宽度对柱壳内表面剪切应力的影响（R_z=50.0 μm）

Figure 13. Effect of surface roughness peak valley element average width on internal surface shear stress of cylindrical shell (R_z=50.0 μm)

下载: 全尺寸图片幻灯片

图12和图13表明，在相同的外爆加载条件下，表面粗糙度的增大以及峰谷单元平均宽度的减小都将导致柱壳内表面剪切应力增大，共同促使柱壳中剪切带的成核与扩展，即柱壳剪切带行为受其内表面粗糙度的影响显著，实验、模拟与理论分析结果一致。

3. 结　论

利用线切割加工技术，获得了具有不同表面粗糙度的45钢圆管，采用外爆加载下的TWC实验技术，结合实验观察、统计分析和有限元模拟，得出如下结论：

(1) 柱壳内表面表征粗糙度的峰谷纵深（R_z）的增大促进了剪切带的成核与扩展，R_z由5.8 μm增大到20.9 μm时，剪切带的数量和最大长度随R_z的增大而增大，剪切带长度差异随R_z的增大而加速增大；

(2) 柱壳内表面表征粗糙度的相邻轮廓峰谷单元的平均宽度（R_sm）减小同样促进剪切带的成核与扩展，R_sm由455 μm减小到179 μm时，剪切带的数量随R_sm减小而增加，剪切带间的相互作用增强；

(3) 柱壳最大剪切应力产生于表面粗糙度峰谷单元中谷底两侧的45°和135°方向，R_z由5 μm增大到50 μm以及R_sm由471 μm减小到157 μm的过程中柱壳内表面的剪切应力均增大，共同促使柱壳中剪切带的成核与扩展。

图 1 45钢样品的粗糙度轮廓线

Figure 1. Roughness contour line of 45 steel sample

下载: 全尺寸图片幻灯片

图 2 45钢圆管的表面形貌

Figure 2. Surface morphology of 45 steel cylindrical shells

下载: 全尺寸图片幻灯片

图 3 TWC实验装置

Figure 3. Thick-walled cylinder experimental device

下载: 全尺寸图片幻灯片

图 4 表面粗糙度的二维模型

Figure 4. Two dimensional model of surface roughness

下载: 全尺寸图片幻灯片

图 5 不同应变下柱壳截面的宏观图像

Figure 5. Macro picture of cylinder cross section under different strains

下载: 全尺寸图片幻灯片

图 6 不同表面粗糙度样品的剪切带分布（ε_eff=0.73±0.01）

Figure 6. Distribution of shear band of samples with different roughnesses (ε_eff=0.73±0.01)

下载: 全尺寸图片幻灯片

图 7 不同表面粗糙度样品的剪切带分布（ε_eff=1.22±0.01）

Figure 7. Distribution of shear band in samples with different roughnesses (ε_eff =1.22±0.01)

下载: 全尺寸图片幻灯片

图 8 R_z=10.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳的等效塑性应变云图（ε_eff=0.66）

Figure 8. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm) of surface roughnesses at R_z=10.0 μm (ε_eff=0.66)

下载: 全尺寸图片幻灯片

图 9 R_z=50.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳等效塑性应变云图（ε_eff=0.66）

Figure 9. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm ) of surface roughnesses at R_z=50.0 μm (ε_eff=0.66)

下载: 全尺寸图片幻灯片

图 10 R_z=10.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳的等效塑性应变云图（ε_eff=1.15）

Figure 10. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm ) of surface roughnesses at R_z=10.0 μm (ε_eff=1.15)

下载: 全尺寸图片幻灯片

图 11 R_z=50.0 μm时具有不同表面粗糙度峰谷单元平均宽度的45钢柱壳的等效塑性应变云图（ε_eff =1.15）

Figure 11. Equivalent plastic strain nephogram of 45 steel cylindrical shell with different average widths of peak valley element (R_sm) of surface roughnesses at R_z=50.0 μm (ε_eff=1.15)

下载: 全尺寸图片幻灯片

图 12 表面粗糙度对柱壳内表面剪切应力的影响（R_sm=157 μm）

Figure 12. Effect of surface roughness on internal surface shear stress of cylindrical shell (R_sm=157 μm)

下载: 全尺寸图片幻灯片

图 13 表面粗糙度峰谷单元平均宽度对柱壳内表面剪切应力的影响（R_z=50.0 μm）

Figure 13. Effect of surface roughness peak valley element average width on internal surface shear stress of cylindrical shell (R_z=50.0 μm)

下载: 全尺寸图片幻灯片

表 1 45钢的化学成分

Table 1. Chemical composition of 45 steel %

C Si Mn S P Ni Cr Cu Fe

0.44 0.23 0.56 0.012 0.014 0.04 0.06 0.07 Rest

下载: 导出CSV

表 2 实验装置的尺寸

Table 2. Size of experimental device

Component Internal diameter/mm Outer diameter/mm Length/mm Thickness/mm

Copper stopper 11 15 67 2
Specimen 15 21 67 3
Copper driver 21 25 67 2
Explosive 25 73 24

下载: 导出CSV

表 3 表面粗糙度模拟参数

Table 3. Simulation parameters of surface roughness

No. R_z/μm R_sm/μm

1 10 471, 235, 157
2 50 471, 235, 157

下载: 导出CSV

表 4 材料的Johnson-Cook本构模型参数^[41]

Table 4. Johnson-Cook constitutive parameters of materials^[41]

Material A/MPa B/MPa C n m

45 steel 507 320 0.064 0.28 1.06
Oxygen free copper 90 292 0.025 0.31 1.09

下载: 导出CSV

表 5 材料的Johnson-Cook失效模型参数

Table 5. Johnson-Cook failure parameters of materials

Material D₁ D₂ D₃ D₄ D₅

45 steel 0.10 0.76 1.57 0.005 −0.84
Oxygen free copper 0.45 3.44 −2.12 0.002 1.09

下载: 导出CSV

表 6 TNT的JWL状态方程参数^[42]

Table 6. Parameters of JWL equation of state of TNT^[42]

A_TNT/GPa B_TNT/GPa R₁ R₂ ω E/(GJ·m⁻³) p_CJ/GPa ρ/(g·cm⁻³) D/(km·s⁻¹)

371.2 3.231 4.51 0.95 0.3 7.0 21 1.63 6.93

下载: 导出CSV

参考文献(42)

[1]	刘明涛, 汤铁钢. 爆炸加载下金属壳体膨胀断裂过程中的关键物理问题 [J]. 爆炸与冲击, 2021, 41(1): 011402. doi: 10.11883/bzycj-2020-0351 LIU M T, TANG T G. Key physical problems in the expanding fracture of explosively driven metallic shells [J]. Explosion and Shock Waves, 2021, 41(1): 011402. doi: 10.11883/bzycj-2020-0351
[2]	朱建士, 陈裕泽. 核武器研制中的力学问题 [J]. 力学与实践, 2002, 24(1): 67–71. doi: 10.3969/j.issn.1000-0879.2002.01.024 ZHU J S, CHEN Y Z. The mechanics of nuclear weapons development [J]. Mechanics in Engineering, 2002, 24(1): 67–71. doi: 10.3969/j.issn.1000-0879.2002.01.024
[3]	MOTT N F. A theory of the fragmentation of shells and bombs [M]//GRADY D. Fragmentation of Rings and Shells: The Legacy of N. F. Mott. Berlin: Springer Publishing, 2006: 243−294.
[4]	任国武, 温上捷, 张茹, 等. 约束层对金属柱壳膨胀变形影响的数值模拟 [J]. 爆炸与冲击, 2017, 37(6): 946–951. doi: 10.11883/1001-1455(2017)06-0946-06 REN G W, WEN S J, ZHANG R, et al. Numerical simulation of influence of constrained layer on expanding deformation of metal cylindrical shell [J]. Explosion and Shock Waves, 2017, 37(6): 946–951. doi: 10.11883/1001-1455(2017)06-0946-06
[5]	杨云川, 朱建军, 郑宇, 等. 战斗部壳体爆炸破片体/线分形维数研究 [J]. 兵工学报, 2018, 39(8): 1499–1506. doi: 10.3969/j.issn.1000-1093.2018.08.006 YANG Y C, ZHU J J, ZHENG Y, et al. Research on the volume and line fractal dimensions of fragments from the explosion of warhead shell [J]. Acta Armamentarii, 2018, 39(8): 1499–1506. doi: 10.3969/j.issn.1000-1093.2018.08.006
[6]	汤铁钢, 胡海波, 李庆忠, 等. 外部爆轰加载过程中金属圆管断裂实验研究 [J]. 爆炸与冲击, 2002, 22(4): 333–337. doi: 10.3321/j.issn:1001-1455.2002.04.008 TANG T G, HU H B, LI Q Z, et al. Studies on the fracture of steel cycinder under external explosive loading [J]. Explosion and Shock Waves, 2002, 22(4): 333–337. doi: 10.3321/j.issn:1001-1455.2002.04.008
[7]	刘龙飞, 周强. 表面粗糙度对6061铝合金薄壁管冲击膨胀断裂性能的影响 [J]. 爆炸与冲击, 2018, 38(4): 749–758. doi: 10.11883/bzycj-2016-0389 LIU L F, ZHOU Q. Effect of surface roughness on impact expansion fracture of 6061 aluminum alloy thin-walled cylindrical tube [J]. Explosion and Shock Waves, 2018, 38(4): 749–758. doi: 10.11883/bzycj-2016-0389
[8]	胡海波, 汤铁钢, 胡八一, 等. 金属柱壳在爆炸加载断裂中的单旋现象 [J]. 爆炸与冲击, 2004, 24(2): 97–107. doi: 10.3321/j.issn:1001-1455.2004.02.001 HU H B, TANG T G, HU B Y, et al. An study of uniform shear bands orientation selection tendency on explosively loaded cylindrical shells [J]. Explosion and Shock Waves, 2004, 24(2): 97–107. doi: 10.3321/j.issn:1001-1455.2004.02.001
[9]	MEYERS M A, NESTERENKO V F, LASALVIA J C, et al. Shear localization in dynamic deformation of materials: microstructural evolution and self-organization [J]. Materials Science and Engineering: A, 2001, 317(1/2): 204–225. doi: 10.1016/S0921-5093(01)01160-1
[10]	罗渝松, 李伟兵, 陈志闯, 等. 内爆加载下金属柱壳的冻结回收方法 [J]. 爆炸与冲击, 2020, 40(10): 104101. doi: 10.11883/bzycj-2020-0041 LUO Y S, LI W B, CHEN Z C, et al. A freezing recovery method for metallic cylinder shells under internal explosive loading [J]. Explosion and Shock Waves, 2020, 40(10): 104101. doi: 10.11883/bzycj-2020-0041
[11]	吴文苍, 董新龙, 庞振, 等. TA2钛合金开口柱壳外爆碎片分布研究 [J]. 力学学报, 2021, 53(6): 1795–1806. doi: 10.6052/0459-1879-21-017 WU W C, DONG X L, PANG Z, et al. Study on fragments distribution of explosively driven cylinders for TA2 titanium alloy [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(6): 1795–1806. doi: 10.6052/0459-1879-21-017
[12]	汤铁钢, 李庆忠, 孙学林, 等. 45钢柱壳膨胀断裂的应变率效应 [J]. 爆炸与冲击, 2006, 26(2): 129–133. doi: 10.11883/1001-1455(2006)02-0129-05 TANG T G, LI Q Z, SUN X L, et al. Strain-rate effects of expanding fracture of 45 steel cylinder shells driven by detonation [J]. Explosion and Shock Waves, 2006, 26(2): 129–133. doi: 10.11883/1001-1455(2006)02-0129-05
[13]	禹富有, 董新龙, 俞鑫炉, 等. 不同填塞装药下金属柱壳断裂特性的实验研究 [J]. 兵工学报, 2019, 40(7): 1418–1424. doi: 10.3969/j.issn.1000-1093.2019.07.011 YU F Y, DONG X L, YU X L, et al. Fracture characteristics of metal cylinder shells with different charges [J]. Acta Armamentarii, 2019, 40(7): 1418–1424. doi: 10.3969/j.issn.1000-1093.2019.07.011
[14]	NESTERENKO V F, LAZARIDI A N, PERSHIN S A. Localization of deformation in copper by explosive compression of hollow cylinders [J]. Fizika Goreniyai Vzryva, 1989, 25(4): 154–155.
[15]	CHEN Y J, MEYERS M A, NESTERENKO V F. Spontaneous and forced shear localization in high-strain-rate deformation of tantalum [J]. Materials Science and Engineering: A, 1999, 268(1/2): 70–82. doi: 10.1016/S0921-5093(99)00110-0
[16]	LOVINGER Z, RITTEL D, ROSENBERG Z. An experimental study on spontaneous adiabatic shear band formation in electro-magnetically collapsing cylinders [J]. Journal of the Mechanics and Physics of Solids, 2015, 79: 134–156. doi: 10.1016/j.jmps.2015.04.007
[17]	杨涛, 刘龙飞, 杨智程, 等. 表面粗糙度对TC4钛合金柱壳剪切带形成的影响 [J]. 力学学报, 2021, 53(3): 813–822. doi: 10.6052/0459-1879-20-433 YANG T, LIU L F, YANG Z C, et al. Effect of surface roughness on the formation of shear band in Ti-6Al-4V alloy cylindrical shell [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(3): 813–822. doi: 10.6052/0459-1879-20-433
[18]	NESTERENKO V F, BONDAR M P. Investigation of deformation localization by the “thick-walled cylinder” method [J]. DYMAT Journal, 1994, 1(3): 245–251.
[19]	NESTERENKO V F, XUE Q, MEYERS M A. Self-organization of shear bands in stainless steel: grain size effect [J]. Journal of Physics (Ⅳ), 2000, 10: 269–274.
[20]	MEYERS M A, XUE Q, NESTERENKO V F. Evolution in the patterning of adiabatic shear bands [J]. AIP Conference Proceedings, 2002, 620(1): 567–570. doi: 10.1063/1.1483602
[21]	黄西成. 内爆与外爆加载下壳体的力学状态及破坏模式分析 [D]. 北京: 中国工程物理研究院, 2010. HUANG X C. Analysis of mechanical states and failure modes of shells subjected to implosive and explosive loadings [D]. Beijing: China Academy of Engineering Physics, 2010.
[22]	YANG Y, LI X M, CHEN S W, et al. Effects of pre-notches on the self-organization behaviors of shear bands in aluminum alloy [J]. Materials Science and Engineering: A, 2010, 527(20): 5084–5091. doi: 10.1016/j.msea.2010.04.079
[23]	任国武, 郭昭亮, 汤铁钢, 等. 高应变率加载下金属柱壳断裂的实验研究 [J]. 兵工学报, 2016, 37(1): 77–82. doi: 10.3969/j.issn.1000-1093.2016.01.012 REN G W, GUO Z L, TANG T G, et al. Experimental research on fracture of metal case under loading at high strain rate [J]. Acta Armamentarii, 2016, 37(1): 77–82. doi: 10.3969/j.issn.1000-1093.2016.01.012
[24]	汤铁钢, 谷岩, 李庆忠, 等. 爆轰加载下金属柱壳膨胀破裂过程研究 [J]. 爆炸与冲击, 2003, 23(6): 529–533. doi: 10.3321/j.issn:1001-1455.2003.06.008 TANG T G, GU Y, LI Q Z, et al. Expanding fracture of steel cylinder shell by detonation driving [J]. Explosion and Shock Waves, 2003, 23(6): 529–533. doi: 10.3321/j.issn:1001-1455.2003.06.008
[25]	张世文, 金山, 刘仓理. 含缺陷厚壁圆管爆轰膨胀断裂的数值模拟 [J]. 应用力学学报, 2010, 27(3): 622–625. ZHANG S W, JIN S, LIU C L. Simulation of the dynamic expanding process of thick-walled cylinder with defects [J]. Chinese Journal of Applied Mechanics, 2010, 27(3): 622–625.
[26]	任国武, 郭昭亮, 张世文, 等. 金属柱壳膨胀断裂的实验与数值模拟 [J]. 爆炸与冲击, 2015, 35(6): 895–900. doi: 10.11883/1001-1455(2015)06-0895-06 REN G W, GUO Z L, ZHANG S W, et al. Experiment and numerical simulation on expansion deformation and fracture of cylindrical shell [J]. Explosion and Shock Waves, 2015, 35(6): 895–900. doi: 10.11883/1001-1455(2015)06-0895-06
[27]	宋鹏飞, 董新龙, 周刚毅, 等. 帽型试样绝热剪切演化实验及数值模拟 [J]. 宁波大学学报(理工版), 2018, 31(2): 50–54. doi: 10.3969/j.issn.1001-5132.2018.02.010 SONG P F, DONG X L, ZHOU G Y, et al. Experimental and numerical analysis of adiabatic shear evolution in hat-shaped specimen [J]. Journal of Ningbo University (Natural Science & Engineering Edition), 2018, 31(2): 50–54. doi: 10.3969/j.issn.1001-5132.2018.02.010
[28]	REMINGTON T P, OWEN J M, NAKAMURA A M, et al. Numerical simulations of laboratory-scale, hypervelocity-impact experiments for asteroid-deflection code validation [J]. Earth and Space Science, 2020, 7(4): e2018EA000474. doi: 10.1029/2018EA000474
[29]	滕凌虹, 曹伟伟, 朱波, 等. ABAQUS在模拟弹丸高低速冲击金属和复合材料靶板方面的应用及研究进展 [J]. 材料导报, 2021, 35(11): 11145–11153. doi: 10.11896/cldb.19110218 TENG L H, CAO W W, ZHU B, et al. Application and development of ABAQUS in simulating high and low velocity impact metallic and composite targets of projectiles [J]. Materials Reports, 2021, 35(11): 11145–11153. doi: 10.11896/cldb.19110218
[30]	吴思思, 董新龙, 俞鑫炉. 45钢柱壳爆炸膨胀断裂的SPH模拟分析 [J]. 爆炸与冲击, 2021, 41(10): 103101. doi: 10.11883/bzycj-2021-0172 WU S S, DONG X L, YU X L. An investigating on explosive expanding fracture of 45 steel cylinders by SPH method [J]. Explosion and Shock Waves, 2021, 41(10): 103101. doi: 10.11883/bzycj-2021-0172
[31]	谢富佩, 徐绯, 曾卓, 等. 复合圆柱壳冲击压缩数值模拟及稳定性研究 [J]. 爆炸与冲击, 2021, 41(11): 112201. doi: 10.11883/bzycj-2020-0431 XIE F P, XU F, ZENG Z, et al. Numerical simulation on stability of composite cylindrical shell under impact compression [J]. Explosion and Shock Waves, 2021, 41(11): 112201. doi: 10.11883/bzycj-2020-0431
[32]	彭克锋, 崔世堂, 潘昊, 等. 冲击载荷作用下柱壳链中的弹性波传播简化模型及其解析解 [J]. 爆炸与冲击, 2021, 41(1): 011403. doi: 10.11883/bzycj-2020-0246 PENG K F, CUI S T, PAN H, et al. Simplified model of elastic wave propagation in cylindrical shell chain under impact load and its analytical solution [J]. Explosion and Shock Waves, 2021, 41(1): 011403. doi: 10.11883/bzycj-2020-0246
[33]	舒旗, 董新龙, 俞鑫炉. 基于Hopkinson压杆的M型试样动态拉伸实验方法研究 [J]. 爆炸与冲击, 2020, 40(8): 084101. doi: 10.11883/bzycj-2019-0433 SHU Q, DONG X L, YU X L. A dynamic tensile method for M-shaped specimen loaded by Hopkinson pressure bar [J]. Explosion and Shock Waves, 2020, 40(8): 084101. doi: 10.11883/bzycj-2019-0433
[34]	刘明涛. 剪切带演化模型及其在模拟柱壳内外爆剪切失稳中的应用 [D]. 合肥: 中国科学技术大学, 2014. LIU M T. A multi-stage probabilistic model for shear band and its application in expanding fracture of cylinder and thick-walled cylinder [D]. Hefei: University of Science and Technology of China, 2014.
[35]	LIU M T, REN G W, FAN C, et al. Experimental and numerical studies on the expanding fracture behavior of an explosively driven 1045 steel cylinder [J]. International Journal of Impact Engineering, 2017, 109: 240–252. doi: 10.1016/j.ijimpeng.2017.07.008
[36]	LOVINGER Z, RITTEL D, ROSENBERG Z. Modeling spontaneous adiabatic shear band formation in electro-magnetically collapsing thick-walled cylinders [J]. Mechanics of Materials, 2018, 116: 130–145. doi: 10.1016/j.mechmat.2017.01.010
[37]	胡昌明, 贺红亮, 胡时胜. 45号钢的动态力学性能研究 [J]. 爆炸与冲击, 2003, 23(2): 188–192. doi: 10.3321/j.issn:1001-1455.2003.02.017 HU C M, HE H L, HU S S. A study on dynamic mechancial behaviors of 45 steel [J]. Explosion and Shock Waves, 2003, 23(2): 188–192. doi: 10.3321/j.issn:1001-1455.2003.02.017
[38]	宫能平, 李贤. 45^#钢动态断裂韧性测试的试验研究 [J]. 安徽理工大学学报(自然科学版), 2007, 27(4): 65–68. GONG N P, LI X. Experimental study of dynamic fracture toughness of 45^# steel [J]. Journal of Anhui University of Science and Technology (Natural Science), 2007, 27(4): 65–68.
[39]	俞鑫炉, 董新龙, 潘顺吉. 不同爆炸载荷下TA2钛合金圆管膨胀破坏过程 [J]. 爆炸与冲击, 2018, 38(1): 148–154. doi: 10.11883/bzycj-2017-0014 YU X L, DONG X L, PAN S J. Fracture behaviors of explosively driven TA2 alloy cylinders under different loadings [J]. Explosion and Shock Waves, 2018, 38(1): 148–154. doi: 10.11883/bzycj-2017-0014
[40]	JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [J]. Engineering Fracture Mechanics, 1983, 21: 541–548.
[41]	汤祁. 面向45^#钢切削残余应力仿真的JC本构模型参数和刀: 屑摩擦系数的确定 [D]. 武汉: 华中科技大学, 2015. TANG Q. Identification of JC constitutive model parameters and tool-chip friction coefficient for 45^# steel cutting residual stress simulation [D]. Wuhan: Huazhong University of Science and Technology, 2015.
[42]	宋浦, 杨凯, 梁安定, 等. 国内外TNT炸药的JWL状态方程及其能量释放差异分析 [J]. 火炸药学报, 2013, 36(2): 42–45. doi: 10.3969/j.issn.1007-7812.2013.02.010 SONG P, YANG K, LIANG A D, et al. Difference analysis on JWL-EOS and energy release of different TNT charge [J]. Chinese Journal of Explosives & Propellants, 2013, 36(2): 42–45. doi: 10.3969/j.issn.1007-7812.2013.02.010