高台阶抛掷爆破速度规律的数值模拟

黄永辉; 刘殿书; 李胜林; 李祥龙; 王家磊

doi:10.11883/1001-1455(2014)04-0495-06

高台阶抛掷爆破速度规律的数值模拟

doi: 10.11883/1001-1455(2014)04-0495-06

1.
昆明理工大学电力工程学院，云南昆明 650500
2.
中国矿业大学(北京) 力学与建筑工程学院，北京 100083
3.
昆明理工大学国土资源工程学院，云南昆明 650093

基金项目: 国家自然科学基金项目(51304087);云南省级人培项目(KKSY201404056, KKSA201121083)

详细信息

作者简介:
黄永辉(1981—), 男, 博士, 讲师

中图分类号: O389;TU313
计量
- 文章访问数: 3938
- HTML全文浏览量: 626
- PDF下载量: 676
- 被引次数: 0
出版历程
- 收稿日期: 2012-11-22
- 修回日期: 2014-06-05
- 刊出日期: 2014-07-25

Numerical simulation on pin-point blasting of sloping surface

1.
Faculty of Electic Power Engineering, Kunming University of Science and Technology, Kunming 650500, Yunnan, China
2.
School of Mechannics and Civil Engineering, China University ofMining and TechnologyBeijing), Beijing 100083, China
3.
Faculty of Land Resources Engineering, Kunming University of Science and Technology, Kunming 650093, Yunnan, China

Funds: Supported bythe National Natural Science Foundation of China (51304087)

More Information

Corresponding author: Huang Yong-hui, 8176309@qq.com

摘要

摘要: 为了获得高台阶抛掷爆破作用下岩石的抛掷速度变化规律，结合现场实验，采用理论分析、高速摄影技术、数值计算方法开展研究。研究结果表明，炮孔内炸药起爆后，坡面岩石抛掷初速度达到最大值的时间在93~105 ms之间，最大抛掷初速度在18~28 m/s之间；坡面岩石的最大抛掷初速度在延炮孔内传爆方向呈增长、稳定、下降的趋势；岩石抛掷运动过程的最后阶段呈自由落体运动，个别岩块的运动速度存在增减现象，主要是由于岩体破碎后岩块间存在的相互碰撞作用；高速摄影实验结果验证了数值计算结果的正确性、RHT材料本构模型及参数在高台阶抛掷爆破数值计算过程中的可用性。
- 爆炸力学 /
- 速度场 /
- RHT材料本构模型 /
- 高台阶抛掷爆破 /
- 高速摄影
Abstract: In order to obtain the casting speed of rock in pin-point blasting of sloping surface, a research was carried out by the field test, theoretical analysis, the high-speed photography and numerical simulation. The results of research show that maximum casting speed of rock was in the range of 18-28 m/s which was reached in 93-105 ms after detonation; along the direction of detonation propagating, maximum casting speed of rock increased to a platform firstly, then declined with the same motion resistance; in the end of pin-point blasting, the rock fall as free-faller, and the casting speed of rock vibrated due to the collision of rock; the process of pin-point blasting of sloping surface was simulated with the RHT constitutive model and parameters assorted. The result of numerical simulation was proved by the measurement of high-speed photography
- mechanics of explosion /
- velocity field /
- RHT material constitutive model /
- pin-point blasting of sloping surface /
- high-speed photography

HTML全文

高台阶抛掷爆破技术主要应用在大型露天矿生产中, 能大幅度降低剥离成本, 提高露天矿的生产水平。对高台阶抛掷爆破技术有了大量研究, 在爆堆形态预测方面应用较成熟的理论为弹道理论模型和WEIBULL模型^[1-2], 然而高台阶抛掷爆破技术成功与否的关键在于岩石抛掷过程中的最大起始速度。由于爆破过程的复杂性, 对爆破现象进行详细的观测比较困难, 利用数值模拟技术, 可以再现爆破过程。软件AUTODYN在求解碰撞冲击、爆炸、冲压成型等瞬态问题方面有独到的优势, 在相关领域已经得到了广泛的应用^[3-4]。本文中, 选用软件AUTODYN, 算法上采用光滑粒子流体动力学方法(smoothed particle hydrodynamics, SPH)和有限元法(finite element method, FEM), 对高台阶抛掷爆破进行模拟, 同时通过现场高速摄影设备及其技术方法, 分析爆炸作用下的台阶坡面岩石速度场分布特征, 并对数值计算过程中所采用的RHT(Riedel-Hiermaier-Thoma)材料本构模型的实用性及其参数的可靠性进行验证。

1. 计算模型及材料参数

1.1 模型构建

以哈尔乌素露天矿高台阶抛掷爆破工程为背景, 构建二分之一等比例的三维高台阶抛掷爆破模型。模型中高台阶坡面及地表面部分定义为自由边界, 其他与岩层接触部位定义为无反射边界; 模型的主体部分采用粒子单元进行网格划分, 考虑粒子单元无法定义边界条件, 在需要定义无反射边界条件的区域构建有限元单元, 采用SPH-FEM耦合的方法进行边界约束, 再对有限元边界定义无反射边界条件。台阶高度35 m, 坡面延伸方向长20 m, 坡面倾角65°, 坡面背侧方向上坡顶及坡底长度分别为17.00和33.29 m; 炮孔布置方面, 模型中构建一个半炮孔, 炮孔采用与台阶坡面平行的倾斜炮孔, 直径310 mm, 孔距11 m, 炮孔至坡面的水平距离为9 m; 炮孔填塞长度7 m; 起爆方式采用孔底起爆。在与对称面交叉部位的台阶坡面上, 由上向下依次在距离台阶顶端7、22、30 m处选定3个高斯记录点。有限元及其他部分布置情况、采用软件AUTODYN构建的二分之一三维模型见图 1。

图 1 高台阶抛掷爆破的三维模型

Figure 1. Three dimensional model of high bench cast blasting

下载: 全尺寸图片幻灯片

1.2 本构模型及参数

RHT材料本构模型综合考虑了在材料破坏过程中所具有的应变硬化、应变率敏感性、压力依赖性和压缩损伤软化等特性, 考虑了压静水和拉静水区应变速率敏感性的差异性, 引入了偏应力张量第三不变量对破坏面形状的影响, 采用了不同的动力放大系数。另外, RHT模型还引入了拉伸损伤, 拉伸和压缩损伤均取决于等效塑性应变, 与材料塑性体积变化无关。RHT模型可应用到弹体贯穿、侵彻、爆炸等数值模拟研究中。该模型引入了最大失效面、弹性极限面和残余失效面3个控制破坏面^[5-7]。失效方程为:

$\sigma_{\mathrm{eq}}^{*}(p, \theta, {\dot \varepsilon })=Y_{\mathrm{TXC}}^{*}(p) R_{3}(\theta) F_{\mathrm{rate}}({\dot \varepsilon })$

(1)

式中: 表示压缩子午线上的等效应力强度, 为应变率强化因子, 分别为:

$Y_{\mathrm{TXC}}^{*}(p)=A\left(p^{*}-p_{\mathrm{spall}}^{*} F_{\mathrm{rate}}({\dot \varepsilon })\right)^{N}$

(2)

${F_{{\rm{rate }}}}(\dot \varepsilon ) = \left\{ {\begin{array}{*{20}{l}} {{{\left( {\frac{{\dot \varepsilon }}{{{{\dot \varepsilon }_0}}}} \right)}^\alpha }}&{p > \frac{{{f_c}}}{3},\quad {{\dot \varepsilon }_0} = 3 \times {{10}^{ - 5}}{{\rm{s}}^{ - 1}}}\\ {{{\left( {\frac{{\dot \varepsilon }}{{{{\dot \varepsilon }_0}}}} \right)}^\alpha }}&{p > \frac{{{f_{\rm{t}}}}}{3},\quad {{\dot \varepsilon }_0} = 3 \times {{10}^{ - 6}}{{\rm{s}}^{ - 1}}} \end{array}} \right.$

(3)

R₃(θ)是相似角θ以及偏平面上失效曲线在拉、压子午线处的偏应力的比Q₂的函数, 可以描述失效强度与偏应力张量第二不变量J₂、第三不变量J₃及净水压力有关的实验结果:

$R_{3}(\theta)=\frac{2\left(1-Q_{2}^{2}\right) \cos \theta+\left(2 Q_{2}-1\right)\left(4\left(1-Q_{2}^{2}\right) \cos ^{2} \theta+5 Q_{2}^{2}-4 Q_{2}\right)^{\frac{1}{2}}}{4\left(1-Q_{2}^{2}\right) \cos ^{2} \theta+\left(1-Q_{2}^{2}\right)^{2}}$

(4)

$\begin{array}{r} \theta=\frac{1}{3} \arccos \left(\frac{3 \sqrt{3} J_{3}}{2 J_{2}^{\frac{1}{2}}}\right) \quad 0 \leqslant \theta \leqslant \frac{\pi}{3} \\ Q_{2}=\frac{r_{1}}{r_{\mathrm{c}}}=Q_{0}+B_{Q} P^{*} \quad 0.51 \leqslant Q_{2} \leqslant 1.00 \end{array}$

式中:量纲一层裂强度为层裂强度; r_t是拉子午线处的偏应力, r_c是压子午线处的偏应力, B_Q是压力影响系数, A、N、α、δ和Q₀是材料常数。

计算中采用的岩石RHT本构模型失效强度参数分别为: 其他参数采用文献[5]中的低强度混凝土参数。

炸药采用铵油炸药, 结合JWL状态方程, 模拟炸药爆炸过程中的压力与体积的关系:

$p=A\left(1-\frac{\omega \eta}{R_{1}}\right) \mathrm{e}^{\frac{R_{1}}{\eta}}+B\left(1-\frac{\omega \eta}{R_{2}}\right) \mathrm{e}^{\frac{R_{1}}{\eta}}+\omega \rho e$

(5)

式中:A、B、R₁、R₂、ω为实验拟合参数; e为比内能; η=ρ/ρ₀, ρ₀为初始密度。采用的炸药材料参数见文献[8]。

2. 数值计算结果

运用软件AUTODYN, 按上述模型尺寸及本构参数构建模型并进行计算, 获得了高台阶坡面体中岩石在抛掷过程中的形态和速度变化规律。在模型的对称面部位进行切片处理, 提取坡面上设置的3个高斯点的速度时程曲线, 如图 2所示; 提取6个时间点的岩石抛掷形态及速度云图, 如图 3所示。

图 2 高斯点速度

Figure 2. The velocity curves of Gauss points

下载: 全尺寸图片幻灯片

图 3 起爆后不同时间岩石运动形态

Figure 3. Variety of damage and deformation of rock when different sequence

下载: 全尺寸图片幻灯片

对计算结果进行分析研究后得到如下结论:

(1) 起爆点设置在孔底部位, 起爆方式采用反向起爆, 整个炮孔内的炸药由起爆点开始起爆至爆轰波传递完成耗时约10 ms, 和装药长度与炸药爆轰速度的比保持一致。

(2) 起爆后, 台阶坡面岩石由静止到产生明显位移时历时约50 ms, 大约在100 ms时抛掷速度达到最大, 最大初速度在22~30 m/s之间, 3个高斯点的最大初速度为22.5、29.6、24.1 m/s。

(3) 台阶坡面岩石运动初期, 各个测点的运动速度差别不大, 随着时间的推移, 在速度增大过程中, 台阶中部岩石移动速度逐渐高出上下两端, 致使岩体整体运动状态呈现中间鼓起、两端略低的形态。

(4) 台阶坡面岩石最大水平抛掷距离约60 m, 整个抛掷过程历时约3.95 s。

3. 高速摄影实验

3.1 实验

哈尔乌素露天矿爆破爆区范围80 m×670 m, 高台阶抛掷爆破现场实验所用的高速摄影设备为Mega Speed Corp公司MS系列CCD和CMOS高速摄影机。为保证人员和仪器的安全, 高速摄影机架设点布置在东观礼台附近(见图 4(a)), 采用手持式GPS获取测点和摄影设备布设点的三维坐标, 确定观测距离为1 073 m, 摄影主光轴倾角和斜角分别为7.01°和17.37°。地面标志点作为高速摄影实验的主要参照物, 为保证在爆破过程中不产生大地位移, 在高速摄影设备一侧距爆区边界50 m处, 在地面树立旗子, 旗子大小120 cm×180 cm(见图 4(b))。

图 4 摄影机架设及地面标志点

Figure 4. Erect high-speed cameras and ground mark

下载: 全尺寸图片幻灯片

3.2 结果

通过现场高速摄影实验, 采用图像分析软件AVI Player, 通过Grid Setings功能在图像上设置网格, 图像分析中采用框标坐标系, 以左右两边的中点连线为x轴, 上下两边的中点连线为y轴, 以两联系交点为原点。读取观测点的坐标, 每个观测点及其坐标测试结果见表 1。

表 1 观测点坐标

Table 1. Imagecoordinates of observation points

t/ms	l=7 m		l=22 m		l=30 m
t/ms	x/m	y/m	x/m	y/m	x/m	y/m
42	37.65	17.57	37.82	13.40	39.36	10.92
95	37.95	17.85	38.15	13.87	39.65	11.25
203	38.25	18.09	38.37	14.05	39.88	11.46
310	38.60	18.37	38.80	14.20	40.01	11.67
428	38.92	18.60	39.12	14.50	40.50	11.91
541	39.21	18.75	39.38	14.73	40.75	12.20
655	39.39	19.05	39.63	15.05	41.17	12.35
778	39.90	19.37	39.95	15.31	41.60	12.55
889	40.55	19.65	40.38	15.43	41.79	12.75
1 009	40.90	19.87	40.80	15.62	42.25	12.93

下载: 导出CSV

| 显示表格

运用高速摄影分析原理及方法^[9], 对高速摄影录像分析后得出如下结论:

(1) 根据高速摄影过程, 将台阶坡面岩石抛掷运动过程分为静止、加速运动和自由落体运动3个阶段。

(2) 根据测点、地面标志点及拍摄点的空间坐标及记录时间的关系, 获得了台阶坡面上3个记录点的位移及速度变化规律, 如图 5所示。台阶坡面岩石抛掷过程中, 由静止到产生明显移动所需时间在33~42 ms之间, 达到最大抛掷速度的时间在93~105 ms之间, 最大初速度在18~28 m/s之间。

图 5 高速摄影速度曲线

Figure 5. The velocity curves from high-speed photography

下载: 全尺寸图片幻灯片

(3) 岩石抛掷速度达到最大值后, 呈现带有起伏波动的下降趋势, 这主要由岩体破碎后岩块间的相互碰撞作用引起的。

4. 结论

(1) 构建了高台阶抛掷爆破的三维数值计算模型, 运用软件AUTODYN对岩石抛掷过程进行了数值模拟, 模拟结果与高速摄影实验结果非常接近, 成功获得了高台阶抛掷爆破岩石抛掷速度变化规律。

(2) 采用高速摄影技术, 对哈尔乌素露天矿高台阶抛掷爆破工程进行现场高速摄影实验, 得到了爆炸瞬间完整的岩石鼓包过程, 将台阶坡面岩石抛掷运动过程分为静止、加速运动和自由落体运动3个阶段。岩石抛掷过程中, 台阶中部岩石移动速度逐渐高出上下两端, 使岩体整体运动状态呈现中间鼓起、两端略低的形态。

(3) 通过对标记点的研究分析, 获得了岩石鼓包过程中的最大抛掷速度, 与数值计算结果一致, 验证了数值计算结果的正确性。说明RHT本构模型在高台阶抛掷爆破数值模拟中具有实用性。

图 1 高台阶抛掷爆破的三维模型

Figure 1. Three dimensional model of high bench cast blasting

下载: 全尺寸图片幻灯片

图 2 高斯点速度

Figure 2. The velocity curves of Gauss points

下载: 全尺寸图片幻灯片

图 3 起爆后不同时间岩石运动形态

Figure 3. Variety of damage and deformation of rock when different sequence

下载: 全尺寸图片幻灯片

图 4 摄影机架设及地面标志点

Figure 4. Erect high-speed cameras and ground mark

下载: 全尺寸图片幻灯片

图 5 高速摄影速度曲线

Figure 5. The velocity curves from high-speed photography

下载: 全尺寸图片幻灯片

表 1 观测点坐标

Table 1. Imagecoordinates of observation points

t/ms	l=7 m		l=22 m		l=30 m
t/ms	x/m	y/m	x/m	y/m	x/m	y/m
42	37.65	17.57	37.82	13.40	39.36	10.92
95	37.95	17.85	38.15	13.87	39.65	11.25
203	38.25	18.09	38.37	14.05	39.88	11.46
310	38.60	18.37	38.80	14.20	40.01	11.67
428	38.92	18.60	39.12	14.50	40.50	11.91
541	39.21	18.75	39.38	14.73	40.75	12.20
655	39.39	19.05	39.63	15.05	41.17	12.35
778	39.90	19.37	39.95	15.31	41.60	12.55
889	40.55	19.65	40.38	15.43	41.79	12.75
1 009	40.90	19.87	40.80	15.62	42.25	12.93

下载: 导出CSV

参考文献(9)

[1]	李祥龙, 何丽华, 栾龙发, 等.露天煤矿高台阶抛掷爆破爆堆形态模拟[J].煤炭学报, 2011, 36(9): 1457-1462. http://www.cqvip.com/QK/96550X/201109/39487592.html Li Xiang-long, He Li-hua, Ruan Long-fa, et al. Simulation model for muckpile shape of high bench cast blasting in surface coal mine[J]. Journal of China Coal Society, 2011, 36(9): 1457-1462. http://www.cqvip.com/QK/96550X/201109/39487592.html
[2]	尤浩生.台阶爆破的计算机模拟[D].北京: 北京科技大学, 1997.
[3]	AUTODYN Theory Manual[Z]. Horsham, UK: Century Dynamics Ltd, 2003.
[4]	杨秀敏.爆炸冲击现象数值模拟[M].合肥: 中国科技大学出版社, 2010.
[5]	Hansson H, Skoglund P. Simulation of concrete penetration in 2D and 3D with the RHT material model[R]. Tumba, Sweden: Swedish Defence Research Agency, 2002.
[6]	Lu Yong, Wang Zhong-qi, Chong K. A comparative study of buried structure in soil subjected to blast loads using 2D and 3D numerical simulations[J]. Soil Dynamics and Earthquake Engineering, 2005, 25(4): 275-288. doi: 10.1016/j.soildyn.2005.02.007
[7]	Brara A, Klepaczko J R. Fracture energy of concrete at high loading rates in tension[J]. International Journal of Impact Engineering, 2007, 34(3): 424-435. doi: 10.1016/j.ijimpeng.2005.10.004
[8]	佟强.中深孔掏槽爆破与矸石块度控制模拟试验研究[D].北京: 中国矿业大学(北京), 2010.
[9]	于起峰.摄像测量学原理与应用研究[M].北京: 科学出版社, 2009.

施引文献

资源附件(0)

访问统计

图(5) / 表(1)

计量

文章访问数: 3938
HTML全文浏览量: 626
PDF下载量: 676
被引次数: 0

1. 计算模型及材料参数
1.1 模型构建
1.2 本构模型及参数
2. 数值计算结果
3. 高速摄影实验
3.1 实验
3.2 结果
4. 结论

高台阶抛掷爆破速度规律的数值模拟

doi: 10.11883/1001-1455(2014)04-0495-06

作者简介: 黄永辉(1981—), 男, 博士, 讲师

计量

出版历程