CN111753401A - 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 - Google Patents
基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 Download PDFInfo
- Publication number
- CN111753401A CN111753401A CN202010468866.0A CN202010468866A CN111753401A CN 111753401 A CN111753401 A CN 111753401A CN 202010468866 A CN202010468866 A CN 202010468866A CN 111753401 A CN111753401 A CN 111753401A
- Authority
- CN
- China
- Prior art keywords
- slope
- blasting
- vibration
- rock
- cantilever beam
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005422 blasting Methods 0.000 title claims abstract description 98
- 239000011435 rock Substances 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000005474 detonation Methods 0.000 claims description 7
- 238000012417 linear regression Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 238000009412 basement excavation Methods 0.000 abstract description 13
- 238000010276 construction Methods 0.000 abstract description 4
- 239000002344 surface layer Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000004141 dimensional analysis Methods 0.000 description 2
- 230000009172 bursting Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域
本发明涉及爆破振动技术领域,特别是涉及一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法。
背景技术
目前,在露天矿山开采、大型水电开挖、岩质深基坑开挖等岩质高边坡工程作业中,仍然以岩体的爆破开挖施工方法为主。其中,岩质高边坡爆破开挖诱发边坡的振动速度是边坡围岩损伤控制的关键指标,可指导实际的爆破工程。
现有常用的爆破振动预测公式为萨道夫斯基预测公式,在平整地形条件下预测地面的爆破振动质点速度具有较高的精度,其具体形式为:
式中K为考虑地质影响及地形影响的系数;α为与地质条件有关的爆破地震波衰减系数;Q为最大单响药量;R为测点与爆源的距离。
由于岩质高边坡地形的爆破振动的传播受到了底部约束的影响,传统的萨道夫斯基预测公式对于岩质高边坡地形的爆破振动速度预测的准确性不高,随着大型水电开挖、岩质深基坑开挖等岩质高边坡工程越来越多,传统萨道夫斯基预测公式已经不能满足实际工程对于爆破振速精确性和安全性的要求。
发明内容
本发明的目的在于至少解决现有技术中存在的技术问题之一,提供一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,能够提高岩质高边坡地形的爆破坡面振动速度预测的精确度。
根据本发明的实施例,提供一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,包括以下步骤:
S1.将岩质高边坡简化为悬臂梁力学模型,建立振动速度预测公式
式中k为考虑地质影响、边坡底部约束及地形影响的系数,β1为与地质条件有关的爆破地震波衰减系数,β2为边坡底部约束影响因子,Q为最大单响药量,R为测点与爆源的距离,h为测点与爆源的高差,H为岩质高边坡总高度;
S2.进行现场爆破试验,确定系数k、β1和β2,得到岩质高边坡爆破振动预测公式的具体形式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S2包括以下步骤:
S210.在边坡坡面进行钻孔装药,记录最大单响药量Q,沿边坡坡面布置爆破振动测试仪,记录爆破振动过程不同高程差位置的震动波形,获得爆破振动速度峰值V,同时测量并记录测点与爆源距离R、测点与爆源高差h、岩质高边坡总高度H;
S220.根据所述步骤S210获得的数据开展二元线性回归分析,确定系数k、β1、β2,以得到岩质高边坡爆破振动预测公式具体形式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S220包括以下步骤:
S223.根据所述线性方程的斜率和截距,求得系数k、β1、β2,将所述k、β1、β2代入振动速度预测公式,得到岩质高边坡爆破坡面的振动速度预测公式的具体形式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S1包括以下步骤:
S101.建立悬臂梁的受迫振动方程;
S102.结合量纲分析方法得到振动速度预测公式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S101中悬臂梁的受迫振动方程为
式中f(t)为等效爆破荷载,q(H)sin(pt)为等效梁受围岩的作用荷载,只考虑第一阶振型,则得到梁上任意一点的最大速度峰值与整个梁上最大速度峰值的比值为sin(πh/2H),作为表征边坡各点速度与整个边坡最大速度峰值比值特性的无量纲量。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S102中,具体步骤为:爆破振动受到测点与爆源的高差h、岩质高边坡总高度H、测点爆心距R、地表岩体质点振动位移μ、地表岩体质点振动加速度a、岩体的固有振动频率f、最大单响药量Q、岩体密度ρ、振动波传播速度c、以及爆轰时间t的影响,由量纲分析得到
V=Φ(Q,R,H,h,f,c,a,t,μ,ρ),
式中取质量Q,速度c,长度R为基本量纲,则有
不同无量纲量πi的乘积和乘方仍为无量纲量,得到
岩质高边坡的振动速度表示为
有益效果:岩质高边坡爆破开挖过程中,岩质高边坡的表层围岩底部受到约束作用,具有典型的悬臂梁特性,因此将岩质高边坡简化为悬臂梁力学模型,据此建立振动速度预测公式,并且通过现场爆破试验记录和测量相关参数,代入振动速度预测公式中获得系数k、β1和β2的值,得到岩质高边坡爆破振动预测公式的具体形式,本基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,将影响边坡振动速度的边坡底部约束考虑在内,提高了岩质高边坡地形的爆破坡面振动速度预测的精确度,可以广泛应用于岩质高边坡的爆破开挖工程中,用以指导爆破施工。
附图说明
下面结合附图对本发明作进一步说明:
图1是本发明实施例现场爆破试验的结构示意图;
图2是本发明实施例爆破振动速度V与高程差h的关系图。
具体实施方式
本部分将详细描述本发明的具体实施例,本发明之较佳实施例在附图中示出,附图的作用在于用图形补充说明书文字部分的描述,使人能够直观地、形象地理解本发明的每个技术特征和整体技术方案,但其不能理解为对本发明保护范围的限制。
在本发明的描述中,需要理解的是,涉及到方位描述,例如上、下、前、后、左、右等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,若干的含义是一个或者多个,多个的含义是两个以上,大于、小于、超过等理解为不包括本数,以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
本发明的描述中,除非另有明确的限定,设置、安装、连接等词语应做广义理解,所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。
本发明实施例提供一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,包括以下步骤:
S1.将岩质高边坡简化为悬臂梁力学模型,建立振动速度预测公式
式中k为考虑地质影响、边坡底部约束及地形影响的系数,β1为与地质条件有关的爆破地震波衰减系数,β2为边坡底部约束影响因子,Q为最大单响药量,R为测点与爆源的距离,h为测点与爆源的高差,H为岩质高边坡总高度。其中,步骤S1包括步骤S101和步骤S102。
S101.通过简化的悬臂梁力学模型建立悬臂梁的受迫振动方程:
式中f(t)为等效爆破荷载,q(H)sin(pt)为等效梁受围岩的作用荷载,通常对于岩质高边坡这类刚度较大的岩土体结构,只考虑第一阶振型,则求解得到梁上任意一点的最大速度峰值与整个梁上最大速度峰值的比值为sin(πh/2H),以此作为表征边坡各点速度与整个边坡最大速度峰值比值特性的无量纲量。
S102.采用量纲分析方法得到振动速度预测公式:
结合上述动力学分析结果可知岩质高边坡的爆破振动响应跟测点与爆源的高差h、岩质高边坡总高度H、以及两者比值的正弦函数相关。
此外,爆破振动还主要受到地形地貌、地质条件、测点爆心距R、地表岩体质点振动位移μ、地表岩体质点振动加速度a、岩体的固有振动频率f、最大单响药量Q、岩体密度ρ、振动波传播速度c、以及爆轰时间t等因素影响,由量纲分析得到:
V=Φ(Q,R,H,h,f,c,a,t,μ,ρ),
式中取质量Q,速度c,长度R为基本量纲(根据π定理取独立变量为(Q,R,c),则有7个π分量,以πi代表无量纲量),则有:
不同无量纲量πi的乘积和乘方仍为无量纲量,得到:
岩质高边坡的振动速度表示为:
S2.进行现场爆破试验,确定未知系数k、β1和β2,得到岩质高边坡爆破振动预测公式的具体形式,具体地,步骤S2包括步骤S210和步骤S220。
S210.参照图1,在边坡坡面进行钻孔装药,记录最大单响药量Q,沿边坡坡面布置爆破振动测试仪,记录爆破振动过程不同高程差位置的震动波形,获得爆破振动速度峰值V,同时测量并记录测点与爆源距离R、测点与爆源高差h、岩质高边坡总高度H,获得至少3组测试数据;
S220.根据步骤S210获得的数据开展二元线性回归分析,确定系数k、β1、β2,以得到岩质高边坡爆破振动预测公式具体形式。其中,步骤S220还包括以下步骤:
S223.根据线性方程的斜率和截距,采用自然常数e为底的指数函数,求得系数k、β1、β2的值,将k、β1、β2的值代入振动速度预测公式,得到岩质高边坡爆破坡面的振动速度预测公式的具体形式。
本实施例的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,主要针对传统的萨道夫斯基预测公式在岩质高边坡爆破坡面振动速度预测上的局限性,岩质高边坡爆破开挖过程中,岩质高边坡的表层围岩底部受到约束作用,具有典型的悬臂梁特性,因此将岩质高边坡简化为悬臂梁力学模型,据此建立振动速度预测公式,并且通过现场爆破试验记录和测量相关参数,代入振动速度预测公式中获得系数k、β1和β2的值,得到岩质高边坡爆破振动预测公式的具体形式,以此修正传统的萨道夫斯基预测公式。
在本实施例中,k=106.10,β1=1.55,β2=0.18,从而得到岩质高边坡爆破振动预测公式具体形式将现场实测值、萨道夫斯基预测值和使用本实施例的爆破振动速度预测公式得到的预测值(即修正公式预测值)绘制于同一表格和折线图中,参照表一和图2。
表一
高程差/m | 实测值m/s | 萨道夫斯基公式预测值m/s | 修正公式预测值m/s |
1 | 4.13 | 4.28 | 3.91 |
2 | 4.01 | 3.96 | 4.63 |
3 | 4.65 | 3.83 | 4.71 |
4 | 4.95 | 3.50 | 4.61 |
5 | 4.59 | 3.23 | 4.41 |
6 | 4.29 | 2.82 | 3.97 |
7 | 3.56 | 2.58 | 3.67 |
8 | 3.23 | 2.37 | 3.40 |
由表一和图2可知,使用本发明实施例的公式的预测精度明显高于目前广泛应用的萨道夫斯基公式,能够实现对岩质高边坡爆破坡面振动速度的精确预测。
本实施例的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,将影响边坡振动速度的边坡底部约束考虑在内,提高了岩质高边坡地形的爆破坡面振动速度预测的精确度,可以广泛应用于岩质高边坡的爆破开挖工程中,用以指导爆破施工,例如根据相关规范(爆破安全规程GB6722-2017)中对结构物的安全振动速度的规定确定V,反推爆破开挖区的最大单响药量Q的限制,具体关系式为:
上面结合附图对本发明实施例作了详细说明,但是本发明不限于上述实施例,在所述技术领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (6)
2.根据权利要求1所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,所述步骤S2包括以下步骤:
S210.在边坡坡面进行钻孔装药,记录最大单响药量Q,沿边坡坡面布置爆破振动测试仪,记录爆破振动过程不同高程差位置的震动波形,获得爆破振动速度峰值V,同时测量并记录测点与爆源距离R、测点与爆源高差h、岩质高边坡总高度H;
S220.根据所述步骤S210获得的数据开展二元线性回归分析,确定系数k、β1、β2,以得到岩质高边坡爆破振动预测公式具体形式。
4.根据权利要求1所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,所述步骤S1包括以下步骤:
S101.建立悬臂梁的受迫振动方程;
S102.结合量纲分析方法得到振动速度预测公式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010468866.0A CN111753401B (zh) | 2020-05-28 | 2020-05-28 | 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010468866.0A CN111753401B (zh) | 2020-05-28 | 2020-05-28 | 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111753401A true CN111753401A (zh) | 2020-10-09 |
CN111753401B CN111753401B (zh) | 2024-04-09 |
Family
ID=72673590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010468866.0A Active CN111753401B (zh) | 2020-05-28 | 2020-05-28 | 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111753401B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112611805A (zh) * | 2020-12-11 | 2021-04-06 | 武汉理工大学 | 一种基于衰减系数的评价围岩松动圈范围的方法 |
CN112632759A (zh) * | 2020-12-14 | 2021-04-09 | 武汉科技大学 | 一种毫秒延时起爆网路理论单响药量的概率计算方法 |
CN114723863A (zh) * | 2022-01-28 | 2022-07-08 | 北京理工大学 | 露天台阶爆破的振动速度云图预测方法、系统及电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104390537A (zh) * | 2014-09-23 | 2015-03-04 | 武汉大学 | 一种基于爆破振动测试的边坡预裂爆破开挖损伤控制方法 |
CN104654948A (zh) * | 2015-01-16 | 2015-05-27 | 武汉科技大学 | 一种边坡与近接隧道协同爆破施工的爆破控制方法 |
-
2020
- 2020-05-28 CN CN202010468866.0A patent/CN111753401B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104390537A (zh) * | 2014-09-23 | 2015-03-04 | 武汉大学 | 一种基于爆破振动测试的边坡预裂爆破开挖损伤控制方法 |
CN104654948A (zh) * | 2015-01-16 | 2015-05-27 | 武汉科技大学 | 一种边坡与近接隧道协同爆破施工的爆破控制方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112611805A (zh) * | 2020-12-11 | 2021-04-06 | 武汉理工大学 | 一种基于衰减系数的评价围岩松动圈范围的方法 |
CN112611805B (zh) * | 2020-12-11 | 2024-04-05 | 武汉理工大学 | 一种基于衰减系数的评价围岩松动圈范围的方法 |
CN112632759A (zh) * | 2020-12-14 | 2021-04-09 | 武汉科技大学 | 一种毫秒延时起爆网路理论单响药量的概率计算方法 |
CN114723863A (zh) * | 2022-01-28 | 2022-07-08 | 北京理工大学 | 露天台阶爆破的振动速度云图预测方法、系统及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111753401B (zh) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111753401A (zh) | 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 | |
CN103150871B (zh) | 利用地下水位与位移实时监测的滑坡预测方法 | |
CN107505043A (zh) | 一种地下水封洞库围岩爆破损伤程度评价方法 | |
CN111814372B (zh) | 一种联合数值计算和块度筛分的爆破块度控制方法 | |
CN105527652A (zh) | 一种岩石脆性的测井方法和装置 | |
CN111751872A (zh) | 地下洞室边墙中心区域爆破振动速度预测系统及方法 | |
CN115660487A (zh) | 一种隧道光面爆破综合质量评价方法 | |
CN111784025B (zh) | 基于简支梁模型的竖井边墙爆破振动速度预测系统和方法 | |
CN107506556B (zh) | 一种确定新鲜完整岩块声波纵波速度值的方法 | |
CN113569487B (zh) | 一种基于bp神经网络预测台阶爆破抛掷效果的方法 | |
Guo et al. | Numerical simulation of surface vibration propagation in tunnel blasting | |
CN110823028B (zh) | 一种优化控制露天矿山台阶爆破爆堆宽度的方法 | |
Zhang et al. | Study on elevation effect of blast wave propagation in high side wall of deep underground powerhouse | |
CN110887419B (zh) | 台阶状边坡受爆破振动影响的监测方法、存储介质和系统 | |
CN112985657A (zh) | 一种高强预应力锚杆支护应力评估方法及相关设备 | |
CN112064617A (zh) | 一种土石混合料地基质量检测方法 | |
Kennedy et al. | Initial Evaluation of the Free-Field Response of the Large Blast Load Generator | |
CN112923879B (zh) | 一种利用声波速度评价洞室围岩松弛厚度及等级的方法 | |
CN116307045B (zh) | 隧道爆破下建筑物振速计算方法、系统、设备及介质 | |
CN111062087A (zh) | 地下工程中基于位移差/梯度的锚固件长度设计方法 | |
CN117784244B (zh) | 基于变纵波速度的细粒混积岩孔隙压力预测方法及系统 | |
CN116793870A (zh) | 一种岩石爆破裂隙区半径概率计算分析方法 | |
US11920472B2 (en) | Reasonable millisecond time control method for excavation blasting of tunnel | |
RU2102697C1 (ru) | Способ определения относительной энергии взрывчатых веществ, переданной грунтовому массиву при взрыве | |
CN115112771A (zh) | 一种测定围岩爆破开挖损伤区的装置及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |