扰动冲击下弱胶结红砂岩的能量耗散与分形特征

张慧梅; 陈世官; 王磊; 程树范; 杨更社; 申艳军

doi:10.11779/CJGE202204004

扰动冲击下弱胶结红砂岩的能量耗散与分形特征 English Version

1.
西安科技大学力学系, 陕西西安 710054
2.
西安科技大学建筑与土木工程学院, 陕西西安 710054
3.
武汉大学土木建筑工程学院, 湖北武汉 430072
4.
西安科技大学地质与环境学院, 陕西西安 710054

基金项目:

国家自然科学基金项目 12172280

国家自然科学基金项目 42077274

国家自然科学基金项目 41907259

陕西省自然科学基金重点项目 2020JZ-53

详细信息

作者简介:
张慧梅（1968—），女，山西大同人，现任教授、博士生导师，主要从事岩土工程稳定性评价及岩石力学理论与应用。E-mail：zhanghuimei68@163.com

通讯作者:
陈世官, E-mail: 18992178070@163.com

中图分类号: TU458.3
计量
- 文章访问数: 248
- HTML全文浏览量: 24
- PDF下载量: 171
出版历程
- 收稿日期: 2021-06-09
- 网络出版日期: 2022-09-22
- 刊出日期: 2022-03-31

Energy dissipation and fractal characteristics of weakly cemented red sandstone under disturbance impact

1.
Department of Mechanics, Xi'an University of Science and Technology, Xi'an 710054, China
2.
College of Architecture and Civil Engineering, Xi'an University of Science and Technology, Xi'an 710054, China
3.
School of Civil and Architectural Engineering, Wuhan University, Wuhan 430072, China
4.
School of Geology and Environment, Xi'an University of Science and Technology, Xi'an 710054, China

摘要

摘要: 为探索动态扰动后西部矿区软岩夹层的能量耗散规律和破坏模式，利用分离式霍普金森压杆装置对弱胶结红砂岩进行动态冲击破坏试验，分析该类红砂岩在受到不同加载速率、不同次数扰动冲击以及是否扰动的条件下，试样在相同加载速率破坏性冲击过程中的能量耗散与分形特征。试验结果表明：在不同速率扰动冲击作用下，随着扰动冲击次数的增加反射能递增而透射能和耗散能呈减小趋势，其中较高速率扰动冲击下试样的反射能高于低速率扰动冲击，耗散能则相反，且低速率扰动冲击下试样的耗散能、能量耗散率和能量耗散密度高于较高速率扰动冲击，表明低速率扰动冲击下试样的能量利用率更高；在破坏性冲击试验中，随着扰动冲击次数的增加，低速率扰动后试样的破碎程度相较于未扰动与高速率扰动更为严重，对应分形维数D_b低速率扰动 > 未扰动 > 高速率扰动，表明分形维数与扰动冲击次数呈正相关，与扰动冲击速率呈负相关；在相同扰动冲击次数下，低速率扰动试样的D_b所对应累积耗散能和耗散能密度高于较高速率扰动试样，而对应的累积反射能则相反。
- 弱胶结砂岩 /
- 扰动冲击 /
- SHPB /
- 能量耗散 /
- 分形维数
Abstract: To explore the energy dissipation law and failure mode of soft rock interlayer in western mining areas of China after dynamic disturbance, the dynamic impact failure tests on weakly cemented red sandstone are carried out by using the separated Hopkinson compression bar device. Under the impact of this red sandstone under different loading rates, different times of disturbance and whether there is the disturbance or not, the energy dissipation and fractal characteristics of the samples during the same loading rate impact failure are analyzed. The experimental results show that under different disturbance impact rates, with the increase of disturbance impact times, the reflection energy increases, while the transmission energy and dissipation energy decrease. The reflection energy of the samples under the impact of high-speed disturbance is higher than that of the low-speed disturbance impact, while the dissipative energy is the opposite. Moreover, the dissipative energy of the samples under the impact of low-speed rate disturbance is opposite. The energy dissipation rate and energy dissipation density are higher than those of the high-speed disturbance impact, which indicates that the energy utilization rate of the samples is higher under the impact of low-speed disturbance. In the impact failure tests, with the increase of the number of disturbance impact, the fragmentation degree of the sample after the low-speed rate disturbance is more serious than that of the undisturbed and high-speed rate disturbance. The low-speed rate disturbance of fractal dimension D_b > undisturbed > high-speed rate disturbance shows that the fractal dimension is positively correlated with the number of disturbance shocks. The results show that the impact rate is negatively correlated with the disturbance. Under the same number of disturbance impact, the cumulative dissipation energy and energy density of D_b of the low-speed rate-disturbed samples are higher than those of the high-speed disturbed samples, while the cumulative reflection energy is opposite.
- weakly cemented sandstone /
- disturbance impact /
- SHPB /
- energy dissipation /
- fractal dimension

HTML全文

0. 引言

随着城市轨道交通的快速发展，列车荷载引起的环境振动问题受到越来越多的关注。不少学者对移动荷载作用下地基的动力响应开展了研究。Eason^[1]通过积分变换法求得了速度小于瑞丽波速的移动点荷载作用下均质弹性半空间位移、应力响应的数值解。Hung等^[2]研究了亚音速、跨音速及超音速3种不同速度范围内的移动荷载作用下黏弹性半空间的动力响应规律，并分析了移动荷载形状分布及加振频率的影响。De等^[3]采用传递反射矩阵法求得了成层弹性半空间在移动荷载的动力格林函数。自Biot^[4-5]提出饱和多孔介质的动力控制方程后，多数学者都基于该控制方程对饱和土体动力响应进行研究。Jin等^[6]利用梯形求积公式获得了振动荷载作用下半无限多孔饱和固体中的应力和孔隙水压力数值解。胡安峰等^[7]利用Fourier积分变换求解了下卧基岩饱和半平面在移动线荷载作用下的动力响应。Xu^[8]等研究了饱和成层土体上无限长Euler-Bernoulli梁在移动荷载作用下的动力响应。Cai等^[9]研究了移动矩形荷载作用下饱和土体的动力响应。

上述研究中，移动荷载大多作用在地基表面，而对地铁等地下构筑物来说，移动荷载作用在地表以下一定深度处。Pak^[10]采用势函数的方法求解了弹性半空间内部点源荷载作用下的动力响应。Senjunctichai等^[11]求得了二维均质多孔弹性半空间在内部荷载激励作用下的动力格林函数。陈胜立等^[12]通过Hankel变换变换方法，求解了埋置点源荷载的轴对称Lame问题。杜秦文等^[13]采用Hankel积分变换法，得到了埋置点源简谐荷载作用下Gibson土体的动力Green函数。上述文献中的埋置荷载多假定为固定位置动荷载。Metrikine等^[14]通过建立一个二维弹性层-梁耦合的简化模型，对隧道中列车运行引起的地表振动规律进行了研究。Yuan等 ^[15]在Metrikine等^[14]的基础上，求解了饱和地基-梁模型在移动荷载作用下的振动解析解。上述两个模型均为二维模型，忽略了垂直于荷载移动方向的地基振动。Forrest等^[16]通过建立筒中筒（PIP）模型，在柱坐标系下对移动点荷载作用下圆形衬砌隧道埋置于弹性全空间中的动力响应情况进行了求解，但未考虑地表边界的影响。本文建立了三维黏弹性半空间地基模型，并通过引入势函数及利用三维傅里叶积分变换及逆变换方法，在笛卡尔坐标系下求解了埋置移动点荷载作用下黏弹性半空间的动力响应积分形式解。最后，通过数值算例分析了埋置荷载移动速度，埋置深度等因素对地表振动传播分布及衰减规律的影响。

1. 控制方程及其通解

均质各向同性弹性体的控制方程以张量形式可表示为（不计体力）：

$\mu {u_{i, jj}} + (\lambda + \mu ) \cdot {u_{j, ji}} = \rho {\ddot u_i} 。$

(1)

式中： ${u_i}$ 为弹性体的位移， $i = x, y, z$ ； $\rho$ 为弹性介质的质量密度； $\lambda$ ， $\mu$ 为弹性体的Lame常数；（^..）为对时间 $t$ 的二阶导数。

根据Helmholtz分解定理，引入标量势函数 $\phi$ ，矢量势函数 ${\psi _j} = ({\psi _1}, {\text{ }}{\psi _2}, {\psi _3})$ 。对位移场作如下分解：

${u_i} = {\phi _{, i}} + {e_{ijk}}{\psi _{k, j}} 。$

(2)

将式（2）代入式（1）可得如下关于势函数的表达式：

${\nabla ^2}\phi - \frac{1}{{{c_{\rm{P}}}^2}} \cdot \frac{{{\partial ^2}\phi }}{{\partial {t^2}}} = 0 \text{, }\$

(3a)

${\nabla ^2}{{\mathit{\Psi}} _j} - \frac{1}{{{c_{\rm{S}}}^2}} \cdot \frac{{{\partial ^2}{{\mathit{\Psi}} _j}}}{{\partial {t^2}}} = 0 。$

(3b)

式中： ${c_{\rm{P}}}$ ， ${c_{\rm{S}}}$ 分别为弹性半空间体的压缩波与剪切波波速： ${c_{\rm{P}}} = \sqrt {\frac{{\lambda + 2\mu }}{\rho }}$ ， ${c_{\rm{S}}} = \sqrt {\frac{\mu }{\rho }}$ ； ${\nabla ^2}$ 为Laplace算子， ${\nabla ^2} = \frac{{{\partial ^2}}}{{\partial {x^2}}} + \frac{{{\partial ^2}}}{{\partial {y^2}}} + \frac{{{\partial ^2}}}{{\partial {z^2}}}$ 。

假定在直角坐标系中矢量势函数满足如下条件：

$\mu {\psi _{j, j}} = 0 。$

(3c)

定义对 $x$ ， $y$ ， $t$ 的三重Fourier变换及逆变换分别为

$\widehat{\overline{\overline{g}}}({k}_{1}, {k}_{2}, z, \omega )={\displaystyle {\int }_{-\infty }^{+\infty }{\displaystyle {\int }_{-\infty }^{+\infty }{\displaystyle {\int }_{-\infty }^{+\infty }g(x, y, z, t)}{\rm{e}}^{-{\rm{i}}{k}_{1}x}{\rm{e}}^{-{\rm{i}}{k}_{2}y}{\rm{e}}^{-i\omega t}{\rm{d}}x{\rm{d}}y{\rm{d}}t}}\text{, }\$

(4a)

$g(x, y, z, t)=\frac{1}{{\left(2{\rm{ \mathsf{ π}}} \right)}^{3}}{\displaystyle {\int }_{-\infty }^{+\infty }{\displaystyle {\int }_{-\infty }^{+\infty }{\displaystyle {\int }_{-\infty }^{+\infty }\widehat{\overline{\overline{g}}}({k}_{1}, {k}_{2}, z, \omega )}{\rm{e}}^{{\rm{i}}{k}_{1}x}{\rm{e}}^{{\rm{i}}{k}_{2}y}{\rm{e}}^{{\rm{i}}\omega t}{\rm{d}}{k}_{1}{\rm{d}}{k}_{2}{\rm{d}}\omega }}。$

(4b)

将式（3a）~（3c）进行式（4a）中的Fourier积分变换可得：

$\frac{{{\rm{d}}}^{2}\widehat{\overline{\overline{\phi }}}}{{\rm{d}}{z}^{2}}-{B}_{\rm{P}}{}^{2}\widehat{\overline{\overline{\phi }}}=0\text{ }\text{, }\$

(5a)

$\frac{{{{\rm{d}}^2}\widehat {\overline{\overline {{{\mathit{\Psi}} _j}}} }}}{{{\rm{d}}{z^2}}} - {B_{\rm{S}}}^2\widehat {\overline{\overline {{{\mathit{\Psi}} _j}}} } = 0 \text{, }\$

(5b)

$\frac{{{\rm{d}}\widehat {\overline{\overline {{{\mathit{\Psi}} _3}}} }}}{{{\rm{d}}z}} = - ({\rm{i}}{k_1}\widehat {\overline{\overline {{{\mathit{\Psi}} _1}}} } + {\rm{i}}{k_2}\widehat {\overline{\overline {{{\mathit{\Psi}} _2}}} }) 。$

(5c)

式中： ${B_{\rm{P}}}^2 = {k_1}^2 + {k_2}^2 - {k_{\rm{P}}}^2$ ； ${B_{\rm{S}}}^2 = {k_1}^2 + {k_2}^2 - {k_{\rm{S}}}^2$ ； ${k_{\rm{P}}} = \frac{\omega }{{{c_{\rm{P}}}}}$ ； ${k_{\rm{S}}} = \frac{\omega }{{{c_{\rm{S}}}}}$ 。其中 ${B_{\rm{p}}}$ ， ${B_{\rm{s}}}$ 的实部大于零。

由弹性介质本构方程及傅里叶变换公式可得波数-频率域内位移及应力响应分量表达式为

${\widehat {\overline{\overline u} }_x} = {\rm{i}}{k_1}\widehat {\overline{\overline \phi } } + {\rm{i}}{k_2}{\widehat {\overline{\overline \psi } }_3} - \frac{{\partial {{\widehat {\overline{\overline \psi } }}_2}}}{{\partial z}} \text{, }\$

(6a)

${\widehat {\overline{\overline u} }_y} = {\rm{i}}{k_2}\widehat {\overline{\overline \phi } } - {\rm{i}}{k_1}{\widehat {\overline{\overline \psi } }_3} + \frac{{\partial {{\widehat {\overline{\overline \psi } }}_1}}}{{\partial z}} \text{, }\$

(6b)

${\widehat {\overline{\overline u} }_z} = \frac{{\partial \widehat {\overline{\overline \phi } }}}{{\partial z}} + {\rm{i}}{k_1}{\widehat {\overline{\overline \psi } }_2} - {\rm{i}}{k_2}{\widehat {\overline{\overline \psi } }_1} \text{, }\$

(6c)

${\widehat{\overline{\overline{\sigma }}}}_{zz}=\lambda \left(\frac{{\partial }^{2}\widehat{\overline{\overline{\phi }}}}{\partial {z}^{2}}-{k}_{1}{}^{2}\widehat{\overline{\overline{\phi }}}-{k}_{2}{}^{2}\widehat{\overline{\overline{\phi }}}\right)+2\mu \left({\rm{i}}{k}_{1}\frac{\partial {\widehat{\overline{\overline{\psi }}}}_{2}}{\partial z}-{\rm{i}}{k}_{2}\frac{\partial {\widehat{\overline{\overline{\psi }}}}_{1}}{\partial z}+\frac{{\partial }^{2}\widehat{\overline{\overline{\phi }}}}{\partial {z}^{2}}\right)\text{, }\$

(6d)

${\widehat{\overline{\overline{\tau }}}}_{xz}=\mu \left(2{\rm{i}}{k}_{1}\frac{\partial \widehat{\overline{\overline{\phi }}}}{\partial z}+{k}_{1}{k}_{2}{\widehat{\overline{\overline{\psi }}}}_{1}-{k}_{1}{}^{2}{\widehat{\overline{\overline{\psi }}}}_{2}-\frac{{\partial }^{2}{\widehat{\overline{\overline{\psi }}}}_{2}}{\partial {z}^{2}}+{\rm{i}}{k}_{2}\frac{\partial {\widehat{\overline{\overline{\psi }}}}_{3}}{\partial z}\right)\text{, }\$

(6e)

${\widehat{\overline{\overline{\tau }}}}_{zy}=\mu \left(2{\rm{i}}{k}_{2}\frac{\partial \widehat{\overline{\overline{\phi }}}}{\partial z}+{k}_{2}{}^{2}{\widehat{\overline{\overline{\psi }}}}_{1}+\frac{{\partial }^{2}{\widehat{\overline{\overline{\psi }}}}_{1}}{\partial {z}^{2}}-{k}_{1}{k}_{2}{\widehat{\overline{\overline{\psi }}}}_{2}-{\rm{i}}{k}_{1}\frac{\partial {\widehat{\overline{\overline{\psi }}}}_{3}}{\partial z}\right)\text{ }。$

(6f)

2. 边值问题

给出了本文的理论研究模型示意图。沿 $x$ 轴正方向匀速移动的竖向点荷载作用在黏弹性半空间体表面以下深度 $z = h$ 处的平面上，移动荷载表达式为 $F = p{{\rm{e}}^{{\rm{i}}{\omega _0}t}}\delta (x - ct)\delta (y)$ ，其在变换域内荷载表达式为 $\widehat {\overline{\overline F} } = \frac{p}{{\sqrt {2{\rm{ \mathsf{ π}}} } }}\delta ({\omega _0} + \omega + c{k_1})$ 。其中p为荷载幅值，c为荷载移动速度， ${\omega _0}$ 为荷载自振频率， $\delta$ 为狄拉克函数。考虑荷载由 $x$ 轴的负无穷处向正无穷处移动，当 $t = 0$ 时，荷载恰好移动到坐标原点处。

图 1 埋置移动荷载作用于弹性半空间内部示意图

Figure 1. Diagram of elastic half-space with embedded moving loads

下载: 全尺寸图片幻灯片

根据波的辐射特性，将弹性半空间体分为上下两个区域：区域1 ${\text{(}}0 \leqslant z \leqslant {h_ - }{\text{)}}$ 与区域2 $({h_ + } \leqslant z \leqslant + \infty )$ ，如图 1所示。则由式（3a）、（3b）可解得变换域内势函数可表示为

当 $0 \leqslant z \leqslant {h_ - }$ 时，

$\left. \begin{array}{l} \widehat{\overline{\overline{{\mathit{\Phi}} }}}={a}_{1}{\rm{e}}^{{B}_{\rm{P}}z}+{a}_{2}{\rm{e}}^{-{B}_{\rm{P}}z}\text{, }\ \\ {\widehat{\overline{\overline{{\mathit{\Psi}} }}}}_{1}={a}_{3}{\rm{e}}^{{B}_{\rm{S}}z}+{a}_{4}{\rm{e}}^{-{B}_{\rm{S}}z}\text{, }\ \\ {\widehat{\overline{\overline{{\mathit{\Psi}} }}}}_{2}={a}_{5}{\rm{e}}^{{B}_{\rm{S}}z}+{a}_{6}{\rm{e}}^{-{B}_{\rm{S}}z}\text{, }\ \\ {\widehat{\overline{\overline{{\mathit{\Psi}} }}}}_{3}=-\frac{{\rm{ie}}^{{B}_{\rm{S}}z}}{{B}_{\rm{S}}}({a}_{3}{k}_{1}+{a}_{5}{k}_{2})+\frac{{\rm{ie}}^{-{B}_{\rm{S}}z}}{{B}_{\rm{S}}}({a}_{4}{k}_{1}+{a}_{6}{k}_{2})。 \end{array} \right\}$

(7a)

当 ${h_ + } \leqslant z \leqslant + \infty$ 时，

$\left. \begin{array}{l} \widehat{\overline{\overline{{\mathit{\Phi}} }}}={a}_{7}{\rm{e}}^{-{B}_{\rm{P}}z}\text{, }\ \\ {\widehat{\overline{\overline{{\mathit{\Psi}} }}}}_{1}={a}_{8}{\rm{e}}^{-{B}_{\rm{S}}z}\text{, }\ \\ {\widehat{\overline{\overline{{\mathit{\Psi}} }}}}_{2}={a}_{9}{\rm{e}}^{-{B}_{\rm{S}}z}\text{, }\ \\ {\widehat{\overline{\overline{{\mathit{\Psi}} }}}}_{3}=\frac{{\rm{ie}}^{-{B}_{\rm{S}}z}}{{B}_{\rm{S}}}({a}_{8}{k}_{1}+{a}_{9}{k}_{2})。 \end{array} \right\}$

(7b)

式中： ${a_j}$ 为未知常数， $j = 1\sim 9$ 。

由图 1可看出，半空间表面的边界条件及荷载作用平面处的连续条件可表示为

当 $z = 0$ 时，

$\left. \begin{array}{l} {\sigma }_{zz}(x, y, 0)=0\text{, }\ \\ {\tau }_{xz}(x, y, 0)=0\text{, }\ \\ {\tau }_{yz}(x, y, 0)=0。 \end{array} \right\}$

(8)

当 $z = h$ 时，

${\sigma _{zz}}(x, y, {h_ + }) - {\sigma _{zz}}(x, y, {h_ - }) = F \text{, }\$

(9a)

${\tau _{xz}}(x, y, {h_ + }) - {\tau _{xz}}(x, y, {h_ - }) = 0 \text{, }\$

(9b)

${\tau _{yz}}(x, y, {h_ + }) - {\tau _{yz}}(x, y, {h_ - }) = 0 \text{, }\$

(9c)

${u_{zz}}(x, y, {h_ + }) - {u_{zz}}(x, y, {h_ - }) = 0 \text{, }\$

(9d)

${u_{xz}}(x, y, {h_ + }) - {u_{xz}}(x, y, {h_ - }) = 0 \text{, }\$

(9e)

${u_{yz}}(x, y, {h_ + }) - {u_{yz}}(x, y, {h_ - }) = 0 。$

(9f)

观察式（7）~（9）可看到，边界条件个数与连续条件个数之和等于未知数个数，故可进行求解。

对边界条件与连续条件进行傅里叶变换后，将式（6a）~（6d）代入，可得到一组关于未知系数 ${a_j}$ 的线性方程组： ${A_{ij}}{a_j} = {f_i}$ ，其中矩阵 ${A_{ij}}$ 与 ${f_i}$ 的表达式见附录1。求得未知系数 ${a_j}$ 后，代入位移表达式, 即可得到变换域内的位移响应解 ${\text{\{ }}{\widehat {\overline{\overline u} }_x}, {\widehat {\overline{\overline u} }_y}, {\widehat {\overline{\overline u} }_z}{\text{\} }}$ 。由Fourier逆变换即可求得稳态响应的表达式：

$\begin{array}{l} \left\{ {{u_x}} \right., {u_y}, {u_z}\} = \frac{1}{{{{(2{\rm{ \mathsf{ π}}} )}^3}}}\int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{{\text{ + }}\infty } {\delta ({\omega _0} + \omega + c{k_1}){\text{\{ }}{{\widehat {\overline{\overline u} }}_x}, {{\widehat {\overline{\overline u} }}_y}, {{\widehat {\overline{\overline u} }}_z}{\text{\} }}} } } \cdot \\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {{\rm{e}}^{{\rm{i}}{k_1}x}}{{\rm{e}}^{{\rm{i}}{k_2}y}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}{k_1}{\rm{d}}{k_2}{\rm{d}}\omega 。 \end{array}$

(10)

3. 数值计算与对比验证

为考虑土体的黏滞性，引入如下黏滞阻尼模型： ${\lambda ^{\text{*}}}{\text{ = }}\lambda (1{\text{ + 2}}i\beta )$ ， ${\mu ^{\text{*}}}{\text{ = }}\mu (1{\text{ + 2}}i\beta )$ ，其中 $\beta$ 为土体介质的黏滞阻尼比。由于式（10）中的动力响应解为无穷积分形式，且其被积函数为一个复杂的振荡函数，故采用IFFT方法进行求解。

3.1 对比验证

当h=0时，本文解可退化为移动荷载作用在黏弹性半空间表面的动力响应解。Hung等^[2]研究了黏弹性半空间在不同类型的地表移动荷载作用下的动力响应。土体参数取值密度 $\rho = 2000{\text{ }}{\rm{kg}}/{{\rm{m}}^3}$ ，泊松比 $\nu = 0.25$ ，阻尼比 $\beta = 0.02$ ，横波波速 ${c_{\rm{s}}} = 100{\text{ }}{\rm{m/s}}$ ，纵波波速 ${c_{\rm{p}}} = 173.2{\text{ }}{\rm{m/s}}$ ，表面波波速 ${c_{\rm{R}}} = 92{\text{ }}{\rm{m/s}}$ 。定义观察点 $(x, y, z) = (0, 0, {z_0})$ （其中 ${z_0} = 1{\text{ }}{\rm{m}}$ ）处的无量纲化竖向位移 $V = 2{{\rm{ \mathsf{ π}}}}\mu {z_{\text{0}}}{u_z}/p$ ，无量纲化水平纵向位移 $W = 2{\rm{ \mathsf{ π}}} \mu {z_{\text{0}}}{u_x}/p$ 。为荷载移动速度 $c = 50{\text{ }}{\rm{m/s}}$ 时，本文位移响应退化解与Hung等^[2]结果的对比图，可见两者吻合得较好。

图 2 本文退化解与Hung等^[2]的对比图

Figure 2. Comparison between degenerated solutions and results of Hung et al^[2]

下载: 全尺寸图片幻灯片

3.2 数值算例分析

定义无量纲位移 $u_x^* = {{{u_x}\lambda {l_0}} / p}$ ， $u_z^* = {{{u_z}\lambda {l_0}} / p}$ ，其中 ${l_0}$ 为特征长度（ ${l_0} = 10{\text{ }}{\rm{m}}$ ）。土体Lame常数 $\lambda = 2 \times {10^7}{\text{ }}{\rm{N}}/{{\rm{m}}^2}$ ， $\mu = 2 \times {10^7}{\text{ }}{\rm{N}}/{{\rm{m}}^2}$ ，土体密度 $\rho =$ $2000{\text{ }}{\rm{kg}}/{{\rm{m}}^3}$ ，黏滞阻尼比 $\beta = 0.02$ ，移动荷载幅值 $p = 1{\text{ }}{\rm{N}}$ 。下面分析埋置移动点荷载作用下黏弹性半空间地表的振动情况。

（1）地表竖向位移响应分布

引入移动坐标系 ${x_t} = x - ct$ 。给出了不同速度的埋置移动荷载作用下，地表无量纲竖向位移 ${u_z}^*$ 在移动坐标系下的曲线分布图，其中 $y = 0$ ，移动荷载埋深分别取h为10，15，20 m。由可看到，地表竖向位移响应关于 ${x_t} = 0$ 基本呈轴对称分布，且荷载埋深越深，振动沿水平方向衰减越快，振动影响范围越小，地表振动幅值越小。此外，对比图 3（a），（b）可发现，荷载速度越大，地表最大竖向位移幅值也越大。

图 3 地表竖向位移

${u_z}^*$

Figure 3. Vertical displacements of surface

${u_z}^*$

下载: 全尺寸图片幻灯片

（2）地表纵向位移响应分布

在埋置移动荷载作用下，地表不仅会产生竖直方向的振动，也会产生水平方向的动力响应。为当荷载埋深 $h = 10{\text{ }}{\rm{m}}$ 时，地表 $x$ 方向的纵向水平位移响应沿移动坐标 ${x_t}$ 的分布图。由图 4可看到埋置移动点荷载作用下，地表纵向水平位移响应曲线关于移动坐标轴呈反对称分布，且垂直于荷载移动方向，距离荷载作用位置越远，即y值越大，纵向水平位移越小；对比图 4（a），（b）可看出，随荷载速度增大，纵向水平位移幅值降低。同时，对比图 3，4可看出，同一埋置荷载作用下，地表土体竖向位移响应幅值远大于纵向水平位移响应幅值，这与实测地面振动规律吻合^[17]。

图 4 地表纵向位移

$u_x^*$

Figure 4. Longitudinal displacements of surface

$u_x^*$

下载: 全尺寸图片幻灯片

（3）地表位移响应的频谱分析

由式（10）可得原点 $(0, 0, 0)$ 处的竖向位移响应表达式为

$\begin{array}{l} {u_z}(0, 0, 0, t) = \frac{1}{{{{{\text{(}}2{{{\rm{ \mathsf{ π}}} )}}}^2}}}\int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{{\text{ + }}\infty } {{{\widehat {\overline{\overline u} }}_z}} ({k_1}, {k_2}, 0, - ct - {\omega _0})} \cdot\\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {{\rm{e}}^{{\rm{i}}{k_1}(x - ct - {\omega _0})}}{{\rm{e}}^{{\rm{i}}{k_2}y}}{\rm{d}}{k_1}{\rm{d}}{k_2} 。 \end{array}$

(11)

则该点处的竖向振动频谱可表示为

$\begin{array}{l} {u_z}(f) = \int_{ - \infty }^{{\text{ + }}\infty } {{u_z}(0, 0, 0, t){{\rm{e}}^{ - {\rm{i}}2{\rm{ \mathsf{ π}}} ft}}} {\rm{d}}t\\ \ \ \ \ \ \ \ = \frac{1}{{2{\rm{ \mathsf{ π}}} }}\int_{ - \infty }^{ + \infty } {{{\widehat {\overline{\overline u} }}_z}\left( {\frac{{2{\rm{ \mathsf{ π}}} {f_0} - 2{\rm{ \mathsf{ π}}} f}}{c}, {k_2}, 0, 2{\rm{ \mathsf{ π}}} f} \right)} {\rm{d}}{k_2} 。 \end{array}$

(12)

同理可得原点处的水平位移频谱表达式 ${u_x}(f)$ 。

图 5给出了两种不同速度的常值移动荷载（h=10 m）作用下，地表原点处竖向位移与水平纵向位移响应的频谱图。可看到，当荷载速度较小时，位移频谱主要分布在频率较低的范围内，且竖向位移与水平位移几乎在同一频率达到峰值。当荷载速度变大时，位移频谱曲线的频率分布范围变大，位移幅值也增大，位移峰值所对应的频率也越大。

图 5 地表位移频谱图

Figure 5. Displacement spectra of surface

下载: 全尺寸图片幻灯片

4. 结论

采用傅里叶变换及Helmholtz分解的方法对埋置移动荷载作用下弹性半空间地基的动力响应进行理论求解，并通过数值计算，分析了荷载埋深、速度等因素对地表位移响应的影响，得到3点结论。

（1）移动荷载埋深越深，地表水平位移与竖向位移响应均越小；荷载移动速度越大，地表竖向位移响应幅值越大，纵向水平位移反而越小。

（2）垂直于荷载移动方向，距离荷载作用位置越远，地表水平位移响应越小。

（3）埋置荷载移动速度越大，位移响应频谱分布范围越大，位移幅值峰值对应的频率也越大。

图 1 红砂岩试样

Figure 1. Red sandstone sample

下载: 全尺寸图片幻灯片

图 2 试验方案流程图

Figure 2. Flow chart of experimental scheme

下载: 全尺寸图片幻灯片

图 3 动态应力平衡验证

Figure 3. Verification of dynamic stress equilibrium

下载: 全尺寸图片幻灯片

图 4 动态冲击应变时程曲线

Figure 4. Curves of strain and time under impact

下载: 全尺寸图片幻灯片

图 5 不同加载速率扰动冲击下入射能量时程曲线

Figure 5. Time-history curves of incident energy under disturbance impact with different loading rates

下载: 全尺寸图片幻灯片

图 6 扰动冲击下砂岩各能量随冲击次数的变化关系

Figure 6. Relationship between energy of sandstone and impact times under disturbance impact

下载: 全尺寸图片幻灯片

图 7 不同加载速率循环扰动后试样损伤形态

Figure 7. Damage morphology of samples after cyclic disturbance at different loading rates

下载: 全尺寸图片幻灯片

图 8 扰动冲击下耗散能密度与冲击次数的关系

Figure 8. Relationship between dissipative energy density and impact times under disturbance impact

下载: 全尺寸图片幻灯片

图 9 扰动与未扰动试样冲击破坏后的破碎形态

Figure 9. Fracture morphology of disturbed and undisturbed samples after impact failure

下载: 全尺寸图片幻灯片

图 10 扰动与未扰动试样冲击破坏后的破碎块度分布

Figure 10. Fragmentation distribution of disturbed and undisturbed samples after impact failure

下载: 全尺寸图片幻灯片

图 11 SHPB冲击试验中lg[M_R/M_T]–lgR曲线

Figure 11. Curve of lg[M_R/M_T]–lgR in SHPB impact tests

下载: 全尺寸图片幻灯片

图 12 分形维数与各累积能量关系曲线

Figure 12. Relationship between fractal dimension and cumulative energy

下载: 全尺寸图片幻灯片

表 1 红砂岩基本物理力学参数

Table 1 Basic physical and mechanical parameters of red sandstone

密度/(g·cm^-3)	纵波波速/(m·s^-1)	单轴抗压强度/MPa	弹性模量/GPa	孔隙度/%
1.874	1897	13.78	1.22	22.8

下载: 导出CSV

表 2 扰动作用下红砂岩冲击破碎块度的筛分结果

Table 2 Screening results of impact fragmentation of red sandstone under disturbing action

试件编号	筛分直径/mm										总质量/g	D_b
试件编号	0.080	0.160	0.315	0.630	1.250	2.500	5.000	10.000	15.000	30.000	总质量/g	D_b
A1-Z	2.18	3.56	1.52	1.20	0.80	3.44	6.48	2.53	35.82	39.18	96.71	2.365
A3-Z	2.78	4.41	1.33	1.20	0.72	2.65	7.46	3.33	33.96	37.35	95.19	2.403
A5-Z	5.49	10.54	3.23	3.35	2.00	6.95	13.39	13.18	9.37	34.07	101.57	2.592
C1-Z	7.84	5.34	1.38	1.35	0.70	3.35	13.89	24.45	36.46	0.00	94.76	2.586
C3-Z	8.68	12.39	3.50	3.43	2.15	6.68	12.68	15.64	30.57	0.00	95.72	2.644
C5-Z	9.40	17.37	4.56	3.35	1.80	7.91	15.73	16.52	16.6	0.00	93.24	2.673
DZ	3.66	10.41	2.95	2.28	1.19	4.11	4.87	8.50	34.84	23.01	95.82	2.537
注：试件编号说明，A和C分别表示较高速率扰动组和低速率扰动组，编号1，3，5表示需要进行扰动的次数，-Z表示每组所有扰动次数完成后的破坏性冲击试验，DZ表示对照组。

下载: 导出CSV

参考文献(26)

[1]	FAN L F, REN F, MA G W. Experimental study on viscoelastic behavior of sedimentary rock under dynamic loading[J]. Rock Mechanics and Rock Engineering, 2012, 45(3): 433–438. doi: 10.1007/s00603-011-0197-7
[2]	LI X B, WENG L. Numerical investigation on fracturing behaviors of deep-buried opening under dynamic disturbance[J]. Tunnelling and Underground Space Technology, 2016, 54: 61–72. doi: 10.1016/j.tust.2016.01.028
[3]	WENG L, HUANG L Q, TAHERI A, et al. Rockburst characteristics and numerical simulation based on a strain energy density index: a case study of a roadway in Linglong gold mine, China[J]. Tunnelling and Underground Space Technology, 2017, 69: 223–232. doi: 10.1016/j.tust.2017.05.011
[4]	谢和平. 分形-岩石力学导论[M]. 北京: 科学出版社, 1996. XIE He-ping. Fractal—Introduction to Rock Mechanics[M]. Beijing: Science Press, 1996. (in Chinese)
[5]	WENG L, WU Z J, LIU Q S, et al. Energy dissipation and dynamic fragmentation of dry and water-saturated siltstones under sub-zero temperatures[J]. Engineering Fracture Mechanics, 2019, 220: 106659. doi: 10.1016/j.engfracmech.2019.106659
[6]	WANG L, QIN Y, JIA H B, et al. Study on mechanical properties and energy dissipation of frozen sandstone under shock loading[J]. Shock and Vibration, 2020(4): 8893128.
[7]	赵忠虎, 谢和平. 岩石变形破坏过程中的能量传递和耗散研究[J]. 四川大学学报(工程科学版), 2008, 40(2): 26–31. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH200802006.htm ZHAO Zhong-hu, XIE He-ping. Energy transfer and energy dissipation in rock deformation and fracture[J]. Journal of Sichuan University (Engineering Science Edition), 2008, 40(2): 26–31. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH200802006.htm
[8]	黎立云, 徐志强, 谢和平, 等. 不同冲击速度下岩石破坏能量规律的实验研究[J]. 煤炭学报, 2011, 36(12): 2007– 2011. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201112009.htm LI Li-yun, XU Zhi-qiang, XIE He-ping, et al. Failure experimental study on energy laws of rock under differential dynamic impact velocities[J]. Journal of China Coal Society, 2011, 36(12): 2007–2011. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201112009.htm
[9]	谢和平, 鞠杨, 黎立云. 基于能量耗散与释放原理的岩石强度与整体破坏准则[J]. 岩石力学与工程学报, 2005, 24(17): 3003–3010. doi: 10.3321/j.issn:1000-6915.2005.17.001 XIE He-ping, JU Yang, LI Li-yun. Criteria for strength and structural failure of rocks based on energy dissipation and energy release principles[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(17): 3003–3010. (in Chinese) doi: 10.3321/j.issn:1000-6915.2005.17.001
[10]	LUNDBERG B, OKROUHLIK M. Efficiency of a percussive rock drilling process with consideration of wave energy radiation into the rock[J]. International Journal of Impact Engineering, 2006, 32(10): 1573–1583. doi: 10.1016/j.ijimpeng.2005.02.001
[11]	HONG L, ZHOU Z L, YIN T B, et al. Energy consumption in rock fragmentation at intermediate strain rate[J]. Journal of Central South University of Technology, 2009, 16(4): 677–682. doi: 10.1007/s11771-009-0112-5
[12]	李夕兵, 周子龙, 叶州元, 等. 岩石动静组合加载力学特性研究[J]. 岩石力学与工程学报, 2008, 27(7): 1387–1395. doi: 10.3321/j.issn:1000-6915.2008.07.011 LI Xi-bing, ZHOU Zi-long, YE Zhou-yuan, et al. Study of rock mechanical characteristics under coupled static and dynamic loads[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(7): 1387–1395. (in Chinese) doi: 10.3321/j.issn:1000-6915.2008.07.011
[13]	FENG J J, WANG E Y, SHEN R X, et al. Investigation on energy dissipation and its mechanism of coal under dynamic loads[J]. Geomechanics and Engineering, 2016, 11(5): 657–670. doi: 10.12989/gae.2016.11.5.657
[14]	MILLON O, RUIZ-RIPOLL M L, HOERTH T. Analysis of the behavior of sedimentary rocks under impact loading[J]. Rock Mechanics and Rock Engineering, 2016, 49(11): 4257–4272. doi: 10.1007/s00603-016-1010-4
[15]	GONG F Q, YE H, LUO Y. The effect of high loading rate on the behaviour and mechanical properties of coal-rock combined body[J]. Shock and Vibration, 2018(6): 1–9.
[16]	张文清, 石必明, 穆朝民. 冲击载荷作用下煤岩破碎与耗能规律实验研究[J]. 采矿与安全工程学报, 2016, 33(2): 375–380. https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201602030.htm ZHANG Wen-qing, SHI Bi-ming, MU Chao-min. Experimental research on failure and energy dissipation law of coal under impact load[J]. Journal of Mining & Safety Engineering, 2016, 33(2): 375–380. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201602030.htm
[17]	王利, 高谦. 基于损伤能量耗散的岩体块度分布预测[J]. 岩石力学与工程学报, 2007, 26(6): 1202–1211. doi: 10.3321/j.issn:1000-6915.2007.06.015 WANG Li, GAO Qian. Fragmentation distribution prediction of rock based on damage energy dissipation[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(6): 1202–1211. (in Chinese) doi: 10.3321/j.issn:1000-6915.2007.06.015
[18]	李成杰, 徐颖, 张宇婷, 等. 冲击荷载下裂隙类煤岩组合体能量演化与分形特征研究[J]. 岩石力学与工程学报, 2019, 38(11): 2231–2241. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201911008.htm LI Cheng-jie, XU Ying, ZHANG Yu-ting, et al. Study on energy evolution and fractal characteristics of cracked coal-rock-like combined body under impact loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(11): 2231–2241. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201911008.htm
[19]	李成杰, 徐颖, 叶洲元. 冲击荷载下类煤岩组合体能量耗散与破碎特性分析[J]. 岩土工程学报, 2020, 42(5): 981–988. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202005027.htm LI Cheng-jie, XU Ying, YE Zhou-yuan. Energy dissipation and crushing characteristics of coal-rock-like combined body under impact loading[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 981–988. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202005027.htm
[20]	戴兵, 罗鑫尧, 单启伟, 等. 循环冲击荷载下含孔洞岩石损伤特性与能量耗散分析[J]. 中国安全科学学报, 2020, 30(7): 69–77. https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK202007011.htm DAI Bing, LUO Xin-yao, SHAN Qi-wei, et al. Analysis on damage characteristics and energy dissipation of rock with a single hole under cyclic impact loads[J]. China Safety Science Journal, 2020, 30(7): 69–77. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK202007011.htm
[21]	邓树新, 王明洋, 李杰, 等. 冲击扰动下滑移型岩爆的模拟试验及机理探讨[J]. 岩土工程学报, 2020, 42(12): 2215–2221. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202012009.htm DENG Shu-xin, WANG Ming-yang, LI Jie, et al. Mechanism and simulation experiment of slip-type rock bursts triggered by impact disturbances[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(12): 2215–2221. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202012009.htm
[22]	谢和平, 鞠杨, 黎立云, 等. 岩体变形破坏过程的能量机制[J]. 岩石力学与工程学报, 2008, 27(9): 1729–1740. doi: 10.3321/j.issn:1000-6915.2008.09.001 XIE He-ping, JU Yang, LI Li-yun, et al. Energy mechanism of deformation and failure of rock masses[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(9): 1729–1740. (in Chinese) doi: 10.3321/j.issn:1000-6915.2008.09.001
[23]	谢和平, 彭瑞东, 鞠杨. 岩石变形破坏过程中的能量耗散分析[J]. 岩石力学与工程学报, 2004, 23(21): 3565–3570. doi: 10.3321/j.issn:1000-6915.2004.21.001 XIE He-ping, PENG Rui-dong, JU Yang. Energy dissipation of rock deformation and fracture[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(21): 3565–3570. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.21.001
[24]	马芹永, 高常辉. 冲击荷载下玄武岩纤维水泥土吸能及分形特征[J]. 岩土力学, 2018, 39(11): 3921–3928, 3968. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201811003.htm MA Qin-yong, GAO Chang-hui. Energy absorption and fractal characteristics of basalt fiber-reinforced cement-soil under impact loads[J]. Rock and Soil Mechanics, 2018, 39(11): 3921–3928, 3968. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201811003.htm
[25]	金解放, 吴越, 张睿, 等. 冲击速度和轴向静载对红砂岩破碎及能耗的影响[J]. 爆炸与冲击, 2020, 40(10): 42–55. https://www.cnki.com.cn/Article/CJFDTOTAL-BZCJ202010004.htm JIN Jie-fang, WU Yue, ZHANG Rui, et al. Effect of impact velocity and axial static stress on fragmentation and energy dissipation of red sandstone[J]. Explosion and Shock Waves, 2020, 40(10): 42–55. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BZCJ202010004.htm
[26]	平琦, 马芹永, 袁璞. 岩石试件SHPB劈裂拉伸试验中能量耗散分析[J]. 采矿与安全工程学报, 2013, 30(3): 401–407. https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201303017.htm PING Qi, MA Qin-yong, YUAN Pu. Energy dissipation analysis of stone specimens in SHPB tensile test[J]. Journal of Mining & Safety Engineering, 2013, 30(3): 401–407. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201303017.htm