Method of theoretically calculating spherical stress wave field in linear-hardening materials under impact load

Wang Zhi-Huan; Jia Lei-Ming; He Zeng; Tian Zhou

doi:10.7498/aps.71.20210954

Abstract

Authors and contacts

FullText HTML

1. 引　言
2. 弹塑性球面应力波场的理论计算方法
2.1 物理模型
2.2 载荷卸载速率对应力波传播的影响
2.3 应力波参量求解
3. 结果与讨论
3.1 恒定载荷作用下的弹塑性应力波
3.2 指数衰减载荷作用下的弹塑性应力波
3.3 理论方法适用性及误差分析
4. 结　论

References

Cited By

Abstract

Based on the linear hardening plastic constitutive model, the theoretical solution of the elastic-plastic spherical stress wave field under impact load is established. Firstly, the influence of the impact load unloading rate on the propagation of the spherical stress wave is analyzed, and three different types of propagation images are obtained. On this basis, the method of theoretically calculating the spherical wave equation in elastic stage, plastic loading stage and unloading stage is established separately, and the calculation scheme of particle displacement, particle velocity, stress and strain is given. Compared with the existing theoretical methods, this method takes into account the different propagation patterns of the stress waves under different unloading rates, and shows how to calculate the stress wave parameters in unloading stage, which is more applicable. This method is used to calculate the elastic-plastic spherical stress wave field under constant shock load and exponential attenuation shock load. The calculated results are in good agreement with those from the existing theoretical method and numerical simulation results in the elastic stage and also in the plastic loading stage. In the unloading stage, the existing theoretical method is no longer applicable, while the results obtained in this paper are in good agreement with the numerical simulation results, which verifies the correctness of the theoretical method.

Keywords:

PACS:

83.60.Uv	(Wave propagation, fracture, and crack healing)
83.10.Gr	(Constitutive relations)
81.40.Lm	(Deformation, plasticity, and creep)

Authors and contacts

1.
Department of Engineering Physics, Tsinghua University, Beijing 100084, China
2.
Northwest Institute of Nuclear Technology, Xi’an 710024, China

Corresponding author: Wang Zhi-Huan, wangzhihuan@nint.ac.cn

FullText HTML

1. 引　言

无限大固体介质中球面应力波的传播问题一直被关注和研究. 基于弹性^[1-9]和黏弹性^[10-15]本构关系, 大量理论研究已经开展, 分析了球面应力波传播特征和规律. 当载荷较高时, 介质发生塑性变形, 需要利用弹塑性本构模型对介质进行描述, 因此, 开展弹塑性介质中球面应力波传播理论研究具有重要的意义.

针对弹塑性球面应力波, 多位学者已经开展了部分理论研究. Yang^[16]利用特征线理论研究了线性硬化塑性材料在线性衰减的冲击载荷作用下球面应力波的传播情况, 分析了塑性剪切模量的影响, 但对卸载区的求解借助了数值方法. Morland^[17]通过假设跨越塑性波波阵面的应力间断的具体表达形式, 限定了外加载荷的形式, 在此条件下给出了弹塑性应力波参量的解析表达式. Milne等^[18,19]研究了恒定冲击载荷、连续卸载冲击载荷和渐增冲击载荷作用下弹塑性球面应力波的传播情况. Rapoport等^[20]对理想弹塑性介质中球面弹塑性波的传播进行了研究, 针对球形空腔壁面施加恒定速度边界条件和恒定压力边界条件两种情况, 给出了弹塑性波参量的解析表达式. Santos等^[21]考虑了黏塑性Perzyna型材料中空腔在恒定载荷作用下的动态响应, 给出了介质中应力的计算公式并与有限元模拟结果进行了对比. 肖建华等^[22]对点爆炸震源产生的地震子波进行了求解. Chen等^[23]研究了爆炸载荷加载下的弹塑性球面应力波问题, 通过求解波动方程, 给出了弹性区和塑性区应力波参量的解析解, 得到了弹塑性边界上的应力时间历程.

目前针对弹塑性介质中球面应力波的理论研究, 都没有考虑载荷卸载速率不同情况下应力波的传播情况. 事实上, 外加载荷卸载速率会影响弹塑性球面应力波传播情况, 载荷卸载速率不同情况下应力波传播过程也存在区别. 利用应力波理论对不同载荷卸载速率情况下的弹塑性球面应力波传播过程进行分析, 可以加深对应力波传播衰减机理的认识, 便于进一步分析弹塑性球面应力波的传播规律. 此外, 现有理论研究缺少对介质进入塑性后卸载区域的计算方法, 也需要进一步补充研究.

本文基于线性硬化塑性本构模型, 结合球面波波动方程, 分析了不同卸载速率的冲击载荷加载下弹塑性球面应力波的传播过程, 并通过解析求解波动方程, 给出了应力波参量的精确解. 通过与现有理论方法、数值模拟结果进行对比, 验证了本文理论模型和方法的正确性, 并对理论方法的适用性和计算过程中引入假设造成的误差进行了分析.

2. 弹塑性球面应力波场的理论计算方法

2.1 物理模型

设无限大的固体介质中有一个半径为r₀的球形空腔, 冲击载荷均匀作用于空腔内壁, t时刻空腔内壁所受压力为p(t). 其中p(t)在t = 0时刻从0突跃至峰值P₀, 随后在t > 0区间内逐渐减小.

在球对称应力应变状态下, 固体材料容变律可以写为

${\sigma _r} + 2{\sigma _\theta } = \left( {3\lambda + 2\mu } \right)\left( {{\varepsilon _r} + 2{\varepsilon _\theta }} \right),$

(1)

其中 ${\sigma _r}$ 和 ${\sigma _\theta }$ 分别为径向应力和环向应力, ${\varepsilon _r}$ 和 ${\varepsilon _\theta }$ 分别为径向应变和环向应变, λ和μ为拉梅常数. 畸变律采用线性硬化塑性模型:

${\sigma _\theta } - {\sigma _r} = 2\mu \left( {{\varepsilon _\theta } - {\varepsilon _r}} \right),{\text{ elastic}},$

(2)

$\begin{split} {\sigma _\theta } - {\sigma _r} =\;& 2{G_{\text{p}}}\left( {{\varepsilon _\theta } - {\varepsilon _r} - {\theta _0}} \right) {+} 2\mu {\theta _0}, \\ & {\text{ plastic, loading}},\end{split}$

(3)

$\begin{split} {\sigma _\theta } - {\sigma _r} =\;& 2\mu \left( {{\varepsilon _\theta } - {\varepsilon _r}} \right) + 2\left( {{G_{\text{p}}} - \mu } \right)\left( {{\theta _{\text{m}}} - {\theta _0}} \right), \\ & {\text{ unloading}}, \end{split}$

(4)

此处忽略了反向塑性加载. 式中 ${\theta _0}$ 为到达初始屈服强度时的偏应变; $Y = 2\mu {\theta _0}$ 为屈服强度; ${\theta _{\text{m}}}$ 为最大偏应变; G_p为塑性剪切模量. 畸变律曲线如图1所示.

$图 1 畸变律曲线\r\nFig. 1. Curve of the distortion law.$

图 1 畸变律曲线

Fig. 1. Curve of the distortion law.

下载: 全尺寸图片幻灯片

在小变形条件下, 球面波传播的波动方程为

$\frac{{\partial {\sigma _r}}}{{\partial r}} + \frac{{2\left( {{\sigma _r} - {\sigma _\theta }} \right)}}{r} = \rho \frac{{{\partial ^2}u}}{{\partial {t^2}}},$

(5)

其中, r和t分别为空间径向坐标和时间坐标, u为介质位移, ρ为介质密度. 几何方程为

${\varepsilon _r} = \frac{{\partial u}}{{\partial r}},$

(6)

${\varepsilon _\theta } = {u}/{r} .$

(7)

模型的初始条件和边界条件分别为

${\left. u \right|_{t = 0}} = 0,$

(8)

${\left. {{\sigma _r}} \right|_{r = {r_0}}} = - P\left( t \right) .$

(9)

2.2 载荷卸载速率对应力波传播的影响

对上述模型进行解析求解的难点在于, 介质本构方程中畸变律在不同阶段分别满足方程(2)—(4), 只有确定了介质处于哪一阶段, 才能代入对应的畸变律方程进行求解.

根据介质本构关系, 弹性球面应力波和塑性球面应力波的传播速度均为常数. 由于 ${G_{\text{p}}} < \mu$ , 所以塑性波波速小于弹性波波速, 在弹塑性球面波传播过程中, 强间断弹性波波阵面和强间断塑性波波阵面在介质中依次向前传播. 尽管载荷是连续卸载的形式, 但外加载荷对应于边界处的径向应力, 而介质的加卸载条件考察的是偏应力, 因此介质中塑性波波阵面后仍可能有一段塑性加载过程, 此加载过程的长短与外加载荷卸载速率相关.

假设介质中任意位置处完成塑性加载过程进入卸载后, 不再进入塑性加载阶段, 则介质中的弹塑性球面应力波传播过程可以分为3种情况, 在r-t平面上表示如图2所示.

$图 2 应力波传播示意图　(a) 载荷卸载较快; (b) 载荷卸载较慢; (c) 载荷卸载速率介于上述二者之间\r\nFig. 2. Schematic diagrams of the stress wave propagation (a) Rapid unloading; (b) slow unloading; (c) middle unloading.$

图 2 应力波传播示意图　(a) 载荷卸载较快; (b) 载荷卸载较慢; (c) 载荷卸载速率介于上述二者之间

Fig. 2. Schematic diagrams of the stress wave propagation (a) Rapid unloading; (b) slow unloading; (c) middle unloading.

下载: 全尺寸图片幻灯片

当外加载荷卸载较快时, 介质在塑性波波阵面后直接进入卸载过程, 没有塑性加载阶段, 在r-t平面上的应力波传播过程如图2(a)所示. 图中的粗实线表示强间断, 细实线表示弱间断. OE为强间断弹性波波阵面, 其斜率为弹性波波速的倒数, 在OE前方为未扰动区. OP为强间断塑性波波阵面, 其斜率为塑性波波速的倒数, 跨越OP后立即进入卸载阶段. 当强间断塑性波波阵面传播至P点时, 偏应力衰减到屈服强度, 塑性波完全衰减, 形成以弹性波波速向前后分别传播的两个弱间断PC和PD, PD在边界发生反射, 成为弱间断DF并向前传播. 根据上述分析, 区域EOr为未扰动区, 区域EOPP₁为弹性区, 在塑性波波阵面OP上发生塑性加载, 区域P₁POt为卸载区.

当外加载荷卸载较慢时, 介质在塑性波波阵面后继续发生塑性加载行为, 且塑性加载区较大, 加卸载边界不会与强间断塑性波波阵面相交, 在r-t平面上的应力波传播过程如图2(b)所示. 初始的应力波传播情况与图2(a)相同, 但跨越OP后介质继续发生塑性加载. 当强间断塑性波波阵面传播至P点时, 偏应力衰减到屈服强度, 塑性波继续以弱间断的形式向前传播. 因为几何扩散效应会使满足屈服条件的塑性应力波随传播而衰减, 从而转化为弹性波, 所以弱间断塑性波的传播速度小于强间断塑性波^[24], 且其传播速度不再是常数. 弱间断塑性波波阵面沿PD传播至D点后完全衰减为弹性波. 边界在C点处开始进入卸载, 并向前传播弱间断CG, CD为加卸载边界.

当外加载荷卸载速率介于上述两种情况之间时, 介质在塑性波波阵面后继续发生塑性加载行为, 但塑性加载区较小, 加卸载边界与强间断塑性波波阵面相交, 在r-t平面上的应力波传播过程如图2(c)所示. 此时加卸载边界CD与塑性波波阵面相交于点D. 跨越OD段后, 介质继续进行塑性加载, 跨越DP段后介质直接进入卸载阶段.

由此确定了不同情况下r-t平面上弹性区、塑性加载区和卸载区的分布情况, 就可以在各区域将对应的本构方程代入波动方程进行解析求解, 得到应力波传播过程中参量的解.

2.3 应力波参量求解

在不同区域求解对应的应力波波动方程, 可以得到应力波参量的解. 弹性区波动方程为

$\frac{{{\partial ^2}u}}{{\partial {r^2}}} + \frac{2}{r}\frac{{\partial u}}{{\partial r}} - \frac{{2u}}{{{r^2}}} = \frac{1}{{{c^2}}}\frac{{{\partial ^2}u}}{{\partial {t^2}}},$

(10)

其中c为弹性波波速:

$c = \sqrt {\frac{{\lambda + 2\mu }}{\rho }},$

(11)

引入位移势函数ϕ, 满足:

$u = \frac{{\partial \phi }}{{\partial r}},$

(12)

则波动方程(10)可化为

$\frac{{{\partial ^2}\left( {r\phi } \right)}}{{\partial {r^2}}} = \frac{1}{{{c^2}}}\frac{{{\partial ^2}\left( {r\phi } \right)}}{{\partial {t^2}}},$

(13)

由于弹性前驱波中只存在右行波, 上式的通解形式为

$\phi \left( {r,t} \right) = \frac{1}{r}f\left( {t - \frac{{r - {r_0}}}{c}} \right),$

(14)

结合边界条件可以得到f的表达式, 从而计算得到介质位移及其他应力波参量.

塑性加载区波动方程为

$\frac{{{\partial ^2}u}}{{\partial {r^2}}} + \frac{2}{r}\frac{{\partial u}}{{\partial r}} - \frac{{2u}}{{{r^2}}} = \frac{1}{{c_{\text{p}}^2}}\frac{{{\partial ^2}u}}{{\partial {t^2}}} + \frac{3}{r}\frac{{{c^2} - c_{\text{p}}^2}}{{c_{\text{p}}^2}}{\theta _0},$

(15)

其中c_p为塑性波波速:

${c_{\text{p}}} = \sqrt {\frac{{3\lambda + 2\mu + 4{G_{\text{p}}}}}{{3\rho }}},$

(16)

利用位移势函数可得:

$\frac{{{\partial ^2}\left( {r\phi } \right)}}{{\partial {r^2}}} = \frac{1}{{c_{\text{p}}^2}}\frac{{{\partial ^2}\left( {r\phi } \right)}}{{\partial {t^2}}} + \frac{{{c^2} - c_{\text{p}}^2}}{{c_{\text{p}}^2}}3{\theta _0}r\ln \frac{r}{{{r_0}}},$

(17)

上式通解形式为

$\begin{split} \phi \left( {r,t} \right) =\;& \frac{1}{r}{f_1}\left( {t - \frac{{r - {r_0}}}{{{c_{\text{p}}}}}} \right) + \frac{1}{r}{f_2}\left( {t + \frac{{r - {r_0}}}{{{c_{\text{p}}}}}} \right) \\&+ \frac{{\left( {{c^2} - c_{\text{p}}^2} \right){r^2}{\theta _0}}}{{6c_{\text{p}}^2}}\left( {3\ln \frac{r}{{{r_0}}} - \frac{5}{2}} \right), \\[-15pt]\end{split}$

(18)

结合边界条件可以得到f₁和f₂的解.

卸载区波动方程为

$\begin{split} &\dfrac{{{\partial ^2}u}}{{\partial {r^2}}} + \dfrac{2}{r}\dfrac{{\partial u}}{{\partial r}} - \dfrac{{2u}}{{{r^2}}}\\ =\;& \dfrac{1}{{{c^2}}}\dfrac{{{\partial ^2}u}}{{\partial {t^2}}} + \dfrac{{c_{\text{p}}^2 - {c^2}}}{{{c^2}}}\left[ {\dfrac{{{\text{d}}{\theta _{m} }}}{{{\text{d}}r}} + \dfrac{{3\left( {{\theta _{m} } - {\theta _0}} \right)}}{r}} \right], \end{split}$

(19)

将位移u进行分解^[16]:

$u\left( {r,t} \right) = U\left( {r,t} \right) + g\left( r \right),$

(20)

其中U满足弹性区波动方程, 可以利用求解弹性区方程的方法进行求解, 但此时需同时考虑左行波和右行波. g满足方程:

$\frac{{{\text{d}}g}}{{{\text{d}}r}} - \frac{g}{r} = \frac{{c_{\text{p}}^2 - {c^2}}}{{{c^2}}}\left( {{\theta _{\text{m}}} - {\theta _0}} \right),$

(21)

结合边界条件易得g的表达式.

利用上述方法, 在应力波传播的r-t图中不同区域分别求解对应的方程, 区域交界处采用位移连续条件, 就可以得到全场位移的解. 将位移对时间求导即为质点速度. 利用几何关系(6)和(7)以及本构关系(1)—(4), 即可得到应变和应力的解. 需要注意的是, 为便于求解, 本文引入了部分假设, 主要包括将图2(b)中弱间断塑性波波阵面PD和加卸载边界CD近似为直线, 忽略特征线DH和HK, 以及将图2(c)中加卸载边界CD近似为直线, 假设CD与左行特征线平行等.

3. 结果与讨论

本节利用弹塑性球面应力波传播过程分析方法和应力波参量解析计算方法, 计算弹塑性应力波参量在传播过程中的变化情况, 并与现有理论方法以及动力学软件AUTODYN计算结果进行对比, 验证本文方法的合理性. 最后对本文方法的适用性及计算过程中引入假设造成的误差进行了分析.

3.1 恒定载荷作用下的弹塑性应力波

Rapoport等^[20]给出了恒定冲击载荷作用下理想弹塑性材料中应力波参量的解析计算方法. 在本文模型中, 取塑性剪切模量G_p = 0, 外加载荷p(t) = P₀为恒定值, 也可以得到此情况下应力波参量的解. 同时, 利用动力学软件ANSYS AUTODYN ALE求解器对相同的模型进行数值模拟. 相关参数取值为: 密度ρ = 2630 kg·m^–3, 弹性模量E = 50 GPa, 泊松比ν = 0.16, 屈服强度Y = 40 MPa, 空腔半径r₀ = 0.2 m, 载荷峰值P₀ = 0.2 GPa.

采用dr = 5.6 × 10^–4 m和dr = 2.8 × 10^–4 m两种网格方案进行数值模拟, 图3给出了两种网格方案条件下计算得到空腔壁面处的位移时间历程和径向应力的空间分布结果, 两种方案条件下的计算结果基本吻合, 验证了网格收敛性, 后续的模拟均采用网格dr = 2.8 × 10^–4 m进行计算.

$图 3 网格无关性　(a) 空腔壁面位移时间历程; (b) 径向应力空间分布\r\nFig. 3. Gird independence: (a) Time history of the cavity surface displacement; (b) radial stress distribution.$

图 3 网格无关性　(a) 空腔壁面位移时间历程; (b) 径向应力空间分布

Fig. 3. Gird independence: (a) Time history of the cavity surface displacement; (b) radial stress distribution.

下载: 全尺寸图片幻灯片

图4给出了利用本文理论方法、文献[20]中理论方法和数值模拟方法得到的不同位置处径向应力时间历程和不同时刻径向应力的空间分布. 由于载荷没有卸载, 此情况下应力波传播过程为图2(b)所示的情况. Rapoport等^[20]没有考虑弱间断塑性应力波的传播和后续卸载过程, 因此其方法在图4(b), 4(c), 4(e)和4(f)中不再适用. 本文中的方法考虑了应力波传播的整个过程, 但在求解波动方程过程中引入了部分假设, 因此得到的结果与Rapoport等^[20]的理论方法以及数值模拟得到的结果有些偏差. 但本文的方法适用范围更广, 对于Rapoport等^[20]的理论方法不适用的区域, 本文方法也可以求解. 在r和t较小的区域, 计算结果不受本文假设近似的影响, 理论方法得到的结果与文献[20]中的理论方法结果完全相同. 在其他区域, 和数值模拟结果能够基本吻合, 验证了理论方法的合理性.

$图 4 径向应力时间历程　(a) r = 0.5 m; (b) r = 0.8 m; (c) r = 1.4 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.3 ms (恒定冲击载荷)\r\nFig. 4. Time history of the radial stress: (a) r = 0.5 m; (b) r = 0.8 m; (c) r = 1.4 m and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.3 ms (constant impact loading).$

图 4 径向应力时间历程　(a) r = 0.5 m; (b) r = 0.8 m; (c) r = 1.4 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.3 ms (恒定冲击载荷)

Fig. 4. Time history of the radial stress: (a) r = 0.5 m; (b) r = 0.8 m; (c) r = 1.4 m and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.3 ms (constant impact loading).

下载: 全尺寸图片幻灯片

3.2 指数衰减载荷作用下的弹塑性应力波

取塑性剪切模量G_p = 12.5 MPa, 外加载荷p(t) = P₀exp(–t/T₀)为指数衰减形式, 其中P₀为载荷峰值, T₀为载荷衰减指数. T₀的大小会影响载荷卸载速率, 进而影响应力波传播过程.

给定T₀, 首先需要判断应力波传播过程属于图2中的哪种情况. 首先按照图2(b)所示情况进行计算并得到C点的时间坐标t_C, 若t_C = 0, 则应力波传播中不包含塑性加载区, 为图2(a)所示的情况. 若t_C > 0, 再判断加卸载边界CD是否与强间断塑性波波阵面OP相交, 若不相交, 则为图2(b)所示的情况, 否则为图2(c)所示的情况.

按照上述方法计算, 载荷衰减指数T₀分别取0.01 ms, 0.2 ms和0.03 ms时, 对应的应力波传播过程分别为图2(a)—2(c)所示的情况. 图5—图7分别给出了衰减指数T₀的3个取值情况下, 本文理论方法和数值模拟方法得到的不同位置处径向应力时间历程和不同时刻径向应力的空间分布.

$图 5 径向应力时间历程　(a) r = 0.3 m; (b) r = 0.5 m; (c) r = 0.6 m和径向应力空间分布 (d) t = 0.05 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T0 = 0.01 ms)\r\nFig. 5. Time history of the radial stress: (a) r = 0.3 m; (b) r = 0.5 m; (c) r = 0.6 m and the radial stress distribution: (d) t = 0.05 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T0 = 0.01 ms).$

图 5 径向应力时间历程　(a) r = 0.3 m; (b) r = 0.5 m; (c) r = 0.6 m和径向应力空间分布 (d) t = 0.05 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.01 ms)

Fig. 5. Time history of the radial stress: (a) r = 0.3 m; (b) r = 0.5 m; (c) r = 0.6 m and the radial stress distribution: (d) t = 0.05 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.01 ms).

下载: 全尺寸图片幻灯片

$图 7 径向应力时间历程　(a) r = 0.4 m; (b) r = 0.6 m; (c) r = 0.7 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T0 = 0.03 ms)\r\nFig. 7. Time history of the radial stress: (a) r = 0.4 m; (b) r = 0.6 m; (c) r = 0.7 m; and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T0 = 0.03 ms).$

图 7 径向应力时间历程　(a) r = 0.4 m; (b) r = 0.6 m; (c) r = 0.7 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.03 ms)

Fig. 7. Time history of the radial stress: (a) r = 0.4 m; (b) r = 0.6 m; (c) r = 0.7 m; and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.03 ms).

下载: 全尺寸图片幻灯片

从上述结果可以看出, 3种情况下本文理论方法计算得到的结果与数值模拟结果能够基本吻合, 验证了应力波传播过程分类和参数理论计算方法的合理性. 在冲击载荷作用下, 介质中首先产生弹塑性双波结构, 随着应力波的传播, 塑性波逐渐衰减, 最终强间断完全消失, 成为弹性波向前传播. 由于几何扩散效应, 在传播过程中弹性波峰值也逐渐衰减. 由于发生了塑性变形, 在应力波传播过后, 塑性区仍会存在残余应力, 残余应力不会衰减, 而是永久留存于介质中. 当载荷衰减指数增大时, 介质中传播的应力波脉宽以及残余应力均会增大.

3.3 理论方法适用性及误差分析

本文的理论计算方法是在小变形假设下得到的, 因此载荷峰值P₀不能过高, 以保证介质满足小变形条件. 当载荷峰值P₀很高时, 介质发生大变形, 控制方程形式将发生改变, 本文方法不再适用.

在对应力波传播过程分析中, 本文假设介质中任意位置处完成塑性加载过程进入卸载后, 不再进入塑性加载阶段, 因此要求冲击载荷P(t)在峰值后为连续卸载形式, 不能出现二次加载. 此外, 对于介质本构模型, 本文忽略了反向塑性加载阶段, 实际问题中介质可能出现反向塑性加载, 但只存在于边界r = r₀附近, 因此当出现反向塑性加载现象时, 本文理论方法的计算结果在边界r = r₀附近不再适用, 但对其他区域仍然适用.

在满足上述物理模型假设条件下, 本文结果与数值模拟结果仍存在一定偏差, 这是由于求解过程中引入的假设造成的. 本文引入假设导致应力波峰值后的计算结果与实际有一定的偏差, 因此定义径向应力误差:

$e = \bigg({{\displaystyle\int_{{x_1}}^{{x_2}} {\left| {{\sigma _r} - {\sigma _{r,{\rm{N}}}}} \right|{\text{d}}x} }}\bigg)\bigg/\bigg({{\displaystyle\int_{{x_1}}^{{x_2}} {\left| {{\sigma _{r,{\rm{N}}}}} \right|{\text{d}}x} }}\bigg),$

(22)

其中 ${\sigma _r}$ 为本文计算得到的径向应力, ${\sigma _{r, {\text{N}}}}$ 为数值模拟得到的径向应力, 对于时间历程取x = t, 空间分布则取x = r, 则e可以反映本文计算结果与数值模拟结果在[x₁, x₂]区间内的平均相对偏差大小. 对图4—图7中的所有结果, 按 (22) 式计算, 得到的误差如表1所示.

表 1 图4—图7中各曲线相对偏差

Table 1. Relative errors of the curves in Figs. 4-7

图4—图7	e/%
图4—图7	(a)	(b)	(c)	(d)	(e)	(f)
图4 恒定冲击载荷	6.1	6.7	6.9	3.3	6.2	11.2
图5 T₀ = 0.01 ms	4.2	8.2	8.8	6.9	8.3	6.4
图6 T₀ = 0.2 ms	6.5	11.2	20.4	2.7	3.1	13.4
图7 T₀ = 0.03 ms	6.6	3.4	4.4	3.8	12.6	7.1

下载: 导出CSV

| 显示表格

从表1可以看出, 大多数计算结果的相对偏差均小于10%. 图4和图6中存在部分结果相对偏差较大, 这两个算例均为载荷卸载较慢的情况, 对应的应力波传播情况为图2(b), 由于直线化了弱间断塑性波波阵面并忽略了卸载区弱间断传播的特征线DH和HK, 导致应力波峰值后的计算结果与实际有所偏差. 图7(e)结果相对偏差较大, 这是由于在其对应的应力波传播情况中, 假设了加卸载边界CD与左行特征线平行, 导致边界r = r₀附近的计算结果与实际有所偏差. r = r₀附近的计算结果与实际有所偏差.

$图 6 径向应力时间历程　(a) r = 0.5 m; (b) r = 0.75 m; (c) r = 1.25 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T0 = 0.2 ms)\r\nFig. 6. Time history of the radial stress: (a) r = 0.5 m; (b) r = 0.75 m; (c) r = 1.25 m and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T0 = 0.2 ms).$

图 6 径向应力时间历程　(a) r = 0.5 m; (b) r = 0.75 m; (c) r = 1.25 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.2 ms)

Fig. 6. Time history of the radial stress: (a) r = 0.5 m; (b) r = 0.75 m; (c) r = 1.25 m and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.2 ms).

下载: 全尺寸图片幻灯片

4. 结　论

在球形空腔内部冲击载荷作用下, 介质中产生弹塑性球面应力波并向前传播. 本文从理论分析的角度研究了弹塑性球面应力波的传播过程, 并建立了应力波参量的解析求解方法.

外加冲击载荷的卸载速率对弹塑性球面应力波传播过程有很大影响. 在r-t平面上, 载荷卸载速率会影响塑性加载区的大小. 当载荷卸载较快时, r-t平面上不存在塑性加载区. 随着载荷卸载速率的降低, r-t平面上开始出现塑性加载区, 但塑性加载区较小时, 加卸载边界会与强间断塑性波波阵面相交. 载荷卸载速率进一步降低时, r-t平面上塑性加载区面积扩大, 加卸载边界不再与强间断塑性波波阵面相交, 因此在强间断塑性波完全衰减后, 还会以弱间断的形式继续传播.

在上述应力波传播分析结果基础上, 利用解析方法对应力波参量进行求解. 在不同载荷加载条件下, 利用本文方法得到应力波参量的解并与Rapoport等^[20]理论方法、数值模拟得到的结果进行对比. 结果显示, 本文方法与现有理论方法、数值模拟方法得到的结果基本吻合, 能够反映应力波传播过程中的主要变化规律, 且绝大多数结果相对数值模拟结果偏差均小于10%. 现有理论方法没有考虑载荷卸载速率对应力波传播的影响, 也没有给出卸载区域的求解方法, 本文理论方法考虑了上述问题, 相对现有方法具有更广泛的适用范围.

References

[1]	Wlodarczyk E, Zielenkiewicz M 2009 J. Theor. App. Mech. 47 127
[2]	Wlodarczyk E, Zielenkiewicz M 2009 J. Theor. App. Mech. 47 761
[3]	Wlodarczyk E, Zielenkiewicz M 2009 Shock Waves 18 465 Google Scholar
[4]	Wlodarczyk E, Zielenkiewicz M 2011 J. Theor. App. Mech. 49 457
[5]	Biswas S 2021 Indian J. Phys. 95 705 Google Scholar
[6]	李孝兰 2000 爆炸与冲击 20 186 Google Scholar Li X L 2000 Explo. Shock Waves 20 186 Google Scholar
[7]	李孝兰 2000 爆炸与冲击 20 283 Google Scholar Li X L 2000 Explo. Shock Waves 20 283 Google Scholar
[8]	肖建华 2001 石油地球物理勘探 36 160 Google Scholar Xiao J H 2001 Oil Geophys. Prospect. 36 160 Google Scholar
[9]	肖建华 2004 石油地球物理勘探 39 249 Google Scholar Xiao J H 2004 Oil Geophys. Prospect. 39 249 Google Scholar
[10]	Xu C J, Chen Q Z, Zhou J, Cai Y Q 2015 J. Civ. Eng. 19 2035 Google Scholar
[11]	卢强, 王占江, 王礼立, 赖华伟, 杨黎明 2013 爆炸与冲击 33 463 Google Scholar Lu Q, Wang Z J, Wang L L, Lai H W, Yang L M 2013 Explo. Shock Waves 33 463 Google Scholar
[12]	卢强, 王占江 2015 物理学报 64 108301 Google Scholar Lu Q, Wang Z J 2015 Acta Phys. Sin. 64 108301 Google Scholar
[13]	Lu Q, Wang Z J 2016 J. Sound Vib. 371 183 Google Scholar
[14]	邓德全, 李兆权 1992 第三届全国岩石动力学学术会议论文集, 桂林, 1992, 第152页 Deng D Q, Li Z Q 1992 Proceedings of the 3rd National Conference on Rock Dynamics Guilin, 1992 p152 (in Chinese)
[15]	赖华伟, 王占江, 杨黎明, 王礼立 2013 爆炸与冲击 33 1 Google Scholar Lai H W, Wang Z J, Yang L M, Wang L L 2013 Explo. Shock Waves 33 1 Google Scholar
[16]	Yang C Y 1970 Int. J. Solids Struct. 6 757 Google Scholar
[17]	Morland L W 1971 J. Mech. Phys. Solids 19 295 Google Scholar
[18]	Milne P C, Morland L W, Yeung W 1988 J. Mech. Phys. Solids 36 15 Google Scholar
[19]	Milne P C, Morland L W 1988 J. Mech. Phys. Solids 36 215 Google Scholar
[20]	Rapoport L, Katzir Z, Rubin M B 2011 Wave Motion 48 441 Google Scholar
[21]	Santos T, Brezolin A, Rossi R, Rodriguez-Martinez J A 2020 Acta Mech. 231 2381 Google Scholar
[22]	肖建华, 孙文涛 1997 石油地球物理勘探 32 809 Xiao J H, Sun W T 1997 Oil Geophys. Prospect. 32 809
[23]	Chen S, Wu J, Zhang Z 2017 J. Eng. Mech. 143 04017034 Google Scholar
[24]	王礼立 2005 应力波基础 (第二版) (北京: 国防工业出版社) 第236页 Wang L L 2005 Foundation of Stress Waves (2nd Ed.) (Beijing: National Defense Industry Press) p236 (in Chinese)

Cited By

期刊类型引用(0)

其他类型引用(2)

图 1 畸变律曲线

Figure 1. Curve of the distortion law.

DownLoad: Full-Size Img PowerPoint

图 2 应力波传播示意图　(a) 载荷卸载较快; (b) 载荷卸载较慢; (c) 载荷卸载速率介于上述二者之间

Figure 2. Schematic diagrams of the stress wave propagation (a) Rapid unloading; (b) slow unloading; (c) middle unloading.

DownLoad: Full-Size Img PowerPoint

图 3 网格无关性　(a) 空腔壁面位移时间历程; (b) 径向应力空间分布

Figure 3. Gird independence: (a) Time history of the cavity surface displacement; (b) radial stress distribution.

DownLoad: Full-Size Img PowerPoint

图 4 径向应力时间历程　(a) r = 0.5 m; (b) r = 0.8 m; (c) r = 1.4 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.3 ms (恒定冲击载荷)

Figure 4. Time history of the radial stress: (a) r = 0.5 m; (b) r = 0.8 m; (c) r = 1.4 m and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.3 ms (constant impact loading).

DownLoad: Full-Size Img PowerPoint

图 5 径向应力时间历程　(a) r = 0.3 m; (b) r = 0.5 m; (c) r = 0.6 m和径向应力空间分布 (d) t = 0.05 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.01 ms)

Figure 5. Time history of the radial stress: (a) r = 0.3 m; (b) r = 0.5 m; (c) r = 0.6 m and the radial stress distribution: (d) t = 0.05 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.01 ms).

DownLoad: Full-Size Img PowerPoint

图 7 径向应力时间历程　(a) r = 0.4 m; (b) r = 0.6 m; (c) r = 0.7 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.03 ms)

Figure 7. Time history of the radial stress: (a) r = 0.4 m; (b) r = 0.6 m; (c) r = 0.7 m; and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.03 ms).

DownLoad: Full-Size Img PowerPoint

图 6 径向应力时间历程　(a) r = 0.5 m; (b) r = 0.75 m; (c) r = 1.25 m和径向应力空间分布 (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.2 ms)

Figure 6. Time history of the radial stress: (a) r = 0.5 m; (b) r = 0.75 m; (c) r = 1.25 m and the radial stress distribution: (d) t = 0.1 ms; (e) t = 0.2 ms; (f) t = 0.6 ms (T₀ = 0.2 ms).

DownLoad: Full-Size Img PowerPoint

表 1 图4—图7中各曲线相对偏差

Table 1. Relative errors of the curves in Figs. 4-7

图4—图7 e/%
(a) (b) (c) (d) (e) (f)

图4 恒定冲击载荷 6.1 6.7 6.9 3.3 6.2 11.2
图5 T₀ = 0.01 ms 4.2 8.2 8.8 6.9 8.3 6.4
图6 T₀ = 0.2 ms 6.5 11.2 20.4 2.7 3.1 13.4
图7 T₀ = 0.03 ms 6.6 3.4 4.4 3.8 12.6 7.1

DownLoad: CSV

[1]	Wlodarczyk E, Zielenkiewicz M 2009 J. Theor. App. Mech. 47 127
[2]	Wlodarczyk E, Zielenkiewicz M 2009 J. Theor. App. Mech. 47 761
[3]	Wlodarczyk E, Zielenkiewicz M 2009 Shock Waves 18 465 Google Scholar
[4]	Wlodarczyk E, Zielenkiewicz M 2011 J. Theor. App. Mech. 49 457
[5]	Biswas S 2021 Indian J. Phys. 95 705 Google Scholar
[6]	李孝兰 2000 爆炸与冲击 20 186 Google Scholar Li X L 2000 Explo. Shock Waves 20 186 Google Scholar
[7]	李孝兰 2000 爆炸与冲击 20 283 Google Scholar Li X L 2000 Explo. Shock Waves 20 283 Google Scholar
[8]	肖建华 2001 石油地球物理勘探 36 160 Google Scholar Xiao J H 2001 Oil Geophys. Prospect. 36 160 Google Scholar
[9]	肖建华 2004 石油地球物理勘探 39 249 Google Scholar Xiao J H 2004 Oil Geophys. Prospect. 39 249 Google Scholar
[10]	Xu C J, Chen Q Z, Zhou J, Cai Y Q 2015 J. Civ. Eng. 19 2035 Google Scholar
[11]	卢强, 王占江, 王礼立, 赖华伟, 杨黎明 2013 爆炸与冲击 33 463 Google Scholar Lu Q, Wang Z J, Wang L L, Lai H W, Yang L M 2013 Explo. Shock Waves 33 463 Google Scholar
[12]	卢强, 王占江 2015 物理学报 64 108301 Google Scholar Lu Q, Wang Z J 2015 Acta Phys. Sin. 64 108301 Google Scholar
[13]	Lu Q, Wang Z J 2016 J. Sound Vib. 371 183 Google Scholar
[14]	邓德全, 李兆权 1992 第三届全国岩石动力学学术会议论文集, 桂林, 1992, 第152页 Deng D Q, Li Z Q 1992 Proceedings of the 3rd National Conference on Rock Dynamics Guilin, 1992 p152 (in Chinese)
[15]	赖华伟, 王占江, 杨黎明, 王礼立 2013 爆炸与冲击 33 1 Google Scholar Lai H W, Wang Z J, Yang L M, Wang L L 2013 Explo. Shock Waves 33 1 Google Scholar
[16]	Yang C Y 1970 Int. J. Solids Struct. 6 757 Google Scholar
[17]	Morland L W 1971 J. Mech. Phys. Solids 19 295 Google Scholar
[18]	Milne P C, Morland L W, Yeung W 1988 J. Mech. Phys. Solids 36 15 Google Scholar
[19]	Milne P C, Morland L W 1988 J. Mech. Phys. Solids 36 215 Google Scholar
[20]	Rapoport L, Katzir Z, Rubin M B 2011 Wave Motion 48 441 Google Scholar
[21]	Santos T, Brezolin A, Rossi R, Rodriguez-Martinez J A 2020 Acta Mech. 231 2381 Google Scholar
[22]	肖建华, 孙文涛 1997 石油地球物理勘探 32 809 Xiao J H, Sun W T 1997 Oil Geophys. Prospect. 32 809
[23]	Chen S, Wu J, Zhang Z 2017 J. Eng. Mech. 143 04017034 Google Scholar
[24]	王礼立 2005 应力波基础 (第二版) (北京: 国防工业出版社) 第236页 Wang L L 2005 Foundation of Stress Waves (2nd Ed.) (Beijing: National Defense Industry Press) p236 (in Chinese)

[1]	Lei Yu, Zhao Dan-Ning, Qiao Hai-Hua. Nonlinear characteristics of variations of Earth’s rotation rate. Acta Physica Sinica, 2024, 73(19): 199101. doi: 10.7498/aps.73.20240815
[2]	Huang Kun, Wang Teng-Fei, Yao Ji. Nonlinear plate theory of single-layered MoS₂ with thermal effect. Acta Physica Sinica, 2021, 70(13): 136201. doi: 10.7498/aps.70.20210160
[3]	Wang Ru-Jia, Wu Shi-Ping, Chen Wei. Propagation of thermoviscoelastic wave in inhomogeneous alloy melt with varying temperature. Acta Physica Sinica, 2019, 68(4): 048101. doi: 10.7498/aps.68.20181923
[4]	Song Xu, Lu Yong-Jun, Shi Ming-Liang, Zhao Xiang, Wang Feng-Hui. Effects of plastic deformation in current collector on lithium diffusion and stress in bilayer lithium-ion battery electrode. Acta Physica Sinica, 2018, 67(14): 140201. doi: 10.7498/aps.67.20180148
[5]	Yin Jian-Wei, Pan Hao, Wu Zi-Hui, Hao Peng-Cheng, Hu Xiao-Mian. A growth study of the Richtmyer-Meshkov flow in the elastoplastic solids under explosive loading. Acta Physica Sinica, 2017, 66(7): 074701. doi: 10.7498/aps.66.074701
[6]	Wang Hong-Ming, Zhu Yi, Li Gui-Rong, Zheng Rui. Plasticity and microstructure of AZ31 magnesium alloy under coupling action of high pulsed magnetic field and external stress. Acta Physica Sinica, 2016, 65(14): 146101. doi: 10.7498/aps.65.146101
[7]	Lu Qiang, Wang Zhan-Jiang. Characteristics of spherical stress wave propagation in the standard linear solid material. Acta Physica Sinica, 2015, 64(10): 108301. doi: 10.7498/aps.64.108301
[8]	Li Jing-Tian, Wang Jian-Lu, Zhang Bang-Qiang, Rong Xi-Ming, Ning Xi-Jing. A statistical model to predict the steady-state creep rate. Acta Physica Sinica, 2014, 63(2): 028101. doi: 10.7498/aps.63.028101
[9]	Zhang Jian-Wen, Ren Yong-Hua, Wu Run-Heng, Feng Tao. The global attractor of nonlinear thermoelastic coupled Sine-Gordon system. Acta Physica Sinica, 2012, 61(11): 110404. doi: 10.7498/aps.61.110404
[10]	Li Xue-Mei, Yu Yu-Ying, Zhang Lin, Li Ying-Hua, Ye Su-Hua, Weng Ji-Dong. Elastic-plastic response of shocked 100 LiF and its window correction at 1550 nm wavelength. Acta Physica Sinica, 2012, 61(15): 156202. doi: 10.7498/aps.61.156202
[11]	Zhang Jian-Wen, Li Jin-Feng, Wu Run-Heng. Global attractor of strongly damped nonlinearthermoelastic coupled rod system. Acta Physica Sinica, 2011, 60(7): 070205. doi: 10.7498/aps.60.070205
[12]	Zhang Feng-Guo, Zhou Hong-Qiang, Zhang Guang-Cai, Hong Tao. Inertial and elastic-plastic effect on spallationdamage of ductile metals. Acta Physica Sinica, 2011, 60(7): 074601. doi: 10.7498/aps.60.074601
[13]	Cen Min, Zhang Yue-Guang, Chen Wei-Lan, Gu Pei-Fu. Influences of deposition rate and oxygen partial pressure on residual stress of HfO2 films. Acta Physica Sinica, 2009, 58(10): 7025-7029. doi: 10.7498/aps.58.7025
[14]	Sun Hong-Xiang, Xu Bai-Qiang, Wang Ji-Jun, Xu Gui-Dong, Xu Chen-Guang, Wang Feng. Numerical simulation of laser-generated Rayleigh wave by finite element method on viscoelastic materials. Acta Physica Sinica, 2009, 58(9): 6344-6350. doi: 10.7498/aps.58.6344
[15]	He An-Min, Shao Jian-Li, Qin Cheng-Sen, Wang Pei. Molecular dynamics study on the plastic behavior of monocrystalline copper under shock loading and unloading. Acta Physica Sinica, 2009, 58(8): 5667-5672. doi: 10.7498/aps.58.5667
[16]	Chen Bin, Peng Xiang-He, Fan Jing-Hong, Sun Shi-Tao, Luo Ji. A thermo-elastoplastic constitutive equation including phase transformation and its applications. Acta Physica Sinica, 2009, 58(13): 29-S34. doi: 10.7498/aps.58.29
[17]	Zhang Liang-Ying, Cao Li, Jin Guo-Xiang. Stochastic resonance of amplitude modulated wave in a linear model of single-mode laser. Acta Physica Sinica, 2006, 55(12): 6238-6242. doi: 10.7498/aps.55.6238
[18]	Bian Bao-Min, Yang Ling, Zhang Ping, Ji Yun-Jing, Li Zhen-Hua, Ni Xiao-Wu. General self-simulating motion mode of spherical strong shock waves in ideal gas. Acta Physica Sinica, 2006, 55(8): 4181-4187. doi: 10.7498/aps.55.4181
[19]	DING HONG-YU, WU BAI-SHI, YIN DA-JUN, LI JIN-QUAN. RATE EQUATION THEORY OF DEGENERATE FOUR WAVE MIXING. Acta Physica Sinica, 1988, 37(3): 408-415. doi: 10.7498/aps.37.408
[20]	LI FU-BIN. PLANE AND SPHERICAL WAVES SOLUTIONS OF THE FIELD EQUATION FOR A COMPLEX SCALAR FIELD WITH A CUBIC NONLINEAR TERM IN FOUR DIMENSIONAL MINKOWSKI SPACE. Acta Physica Sinica, 1988, 37(9): 1461-1470. doi: 10.7498/aps.37.1461

期刊类型引用(0)
其他类型引用(2)

Catalog

Vol. 71, Issue 1 - 2022-01-05

Metrics

Abstract views: 6800
PDF Downloads: 117
Cited By: 2

姓名
邮箱
手机号码
标题
留言内容
验证码

Search

Article

留言板