D0I:10.13374/i.issm1001053x.2003.06.005 第25卷第6期 北京科技大学学报 Vol.25 No.6 2003年12月 Journal of University of Science and Technology Beijing Dee.2003 爆破震动对高速路边坡影响的数值模拟 林士炎)李长洪)乔兰)刘宝许 1)河北省政府办公厅建设环保处,石家庄0500532)北京科技大学土木与环境工程学院,北京100083 摘要以京秦高速公路某段路堑边坡爆破开挖为例,采用数值模拟的方法分析了爆破震 动波在边坡介质中的应力场、位移场、速度场等的分布规律,多方面研究了燥破震动效应对 山区高等级公路岩质边坡稳定性的影响,现场振动测试结果验证了数值模拟计算的可靠性, 关键词路堑开挖:爆破震动:岩质边坡:数值模拟:稳定评价 分类号TU451:TD235 边坡稳定性问题是许多山区高速公路工程 挖,地表岩层基本全部棵露,开挖前旧路边坡经 所面临的重要技术问题,来自路堑边坡开挖爆破 常有松动的岩块滑落,坡面岩体节理裂隙发育, 所产生的震动效应,作为一种影响边坡稳定的 多为垂向节理,边坡剖面示意图见图1. 外在因素,得到日益广泛的关注.在不稳定或临 界稳定的边坡上进行边坡施工,爆破震动效应显 图例 得尤为重要,有必要对爆破震动效应进行综合研 Ars]片麻岩 究,将爆破震动对边坡稳定性的影响降到最低。 k5+960钻孔及 京秦高速公路青龙连接线所经区域地处燕 最终边坡 20编号标高 山山脉东段,属山岭重丘区,穿越的地形、地貌及 Ys2花岗岩 地质条件环境复杂,沿线多数路段需采用爆破的 方法进行开挖施工四.工程完工后将出现深路堑 图1边坡剖面示意图 Fig.1 Schematic diagram of the slope section plane 边坡,台阶为具有三个自由面的半路堑开挖体, 上部预留边坡高度最高可达45m,山顶部位有古 1.2岩体力学参数 长城遗址需保护,路基则为高台阶临河边坡 岩体力学参数根据边坡岩样物理力学性质 实验得到,研究路段岩性较为均一,主要是混合 1边坡地质环境特征 岩,多组岩样力学性质基本相似,在模拟计算过 1.1地层岩性 程中采用的物理力学参数取其均值见表1, 采用爆破工程施工路段出露的地层多以燕 山期侵入岩为主,代表岩性为花岗岩、闪长岩:其 表1岩体物理力学参数 Table 1 Parameters of physical-mechanical properties of 次为太古界变质岩,代表岩性为片麻岩、混合岩 the rock 及碎裂状大理岩等:再次为第四系松散堆积物 p/(g.cm)V/(m.s)o/MPa d/MPa E/GPa 花岗岩多为强风化似斑状粗粒结构,全路段穿越 2.89 6520.9 228.159.80 93.49 地形和地质构造条件复杂,有包括褶皱、大断裂 G/GPa C/MPa) 和各种节理裂隙在内的多种地质构造.边坡岩性 118.88 0.283 33.17 43 以碎裂状混合岩和闪长玢岩为主,垂直节理较为 发育,岩石硬度较大,该路段主要沿旧路进行开 2爆破震动数值模拟 2.1爆破数值分析模型的选择 收稿日期2003-03-17 林士炎男,40岁,项士 *国家自然科学基金资助项目(No.50074002) 边坡开挖数值模拟分析采用了DYNA-2D程
第 卷 第 期 年 月 北 京 科 技 大 学 学 报 心】 爆破震动对高速路边坡影响的数值模拟 林士 炎 ‘, 李长 洪 , 乔 兰 ” 刘 宝 许 ” 河 北省政府 办公 厅 建设 环 保处 , 石 家庄 北京 科技 大 学土 木 与环 境工 程 学 院 , 北京 摘 要 以京秦 高速 公 路 某 段 路 堑 边坡 爆破 开 挖 为 例 , 采用 数 值模拟 的方 法 分 析 了爆 破 震 动 波在 边 坡介质 中的应 力场 、 位 移 场 、 速 度场 等 的分 布 规 律 , 多方 面研 究 了爆 破 震 动 效应 对 山区 高等 级 公 路岩质 边 坡 稳 定 性 的影 响 现 场 振 动测 试 结果 验 证 了数 值 模拟 计 算 的可 靠 性 关键词 路 堑 开 挖 爆 破 震 动 岩质 边 坡 数值模拟 稳 定评 价 分 类号 边 坡 稳 定 性 问题 是 许 多 山 区 高速 公 路 工 程 所 面 临 的重 要技术 问题 , 来 自路 堑边坡 开 挖爆破 所 产 生 的震 动 效 应 ‘ , 作 为 一 种 影 响边 坡 稳 定 的 外在 因素 , 得 到 日益 广 泛 的关 注 在 不 稳 定 或 临 界稳 定 的边 坡 上进 行 边坡 施 工 , 爆破 震 动 效应 显 得 尤 为重要 , 有 必 要对 爆破 震动 效应进行 综合研 究 , 将 爆 破 震 动 对 边 坡 稳 定性 的影 响 降 到最 低 京 秦 高速 公 路 青 龙 连 接 线 所 经 区 域 地 处 燕 山 山脉 东段 , 属 山岭 重丘 区 , 穿越 的地 形 、 地 貌及 地 质条件环 境 复杂 , 沿 线 多数路 段 需采用 爆破 的 方 法 进 行 开挖 施 工 ‘ , 工 程 完 工 后 将 出现 深 路 堑 边 坡 , 台阶 为 具 有 三 个 自由面 的半 路 堑 开 挖 体 , 上 部 预 留边 坡 高度 最 高可达 , 山顶 部位 有 古 长 城 遗 址 需保 护 , 路 基 则 为 高 台阶 临河 边 坡 挖 , 地 表 岩 层 基 本 全 部裸 露 , 开 挖 前 旧 路 边 坡 经 常 有松 动 的岩 块 滑 落 , 坡 面 岩 体 节 理 裂 隙发育 , 多 为垂 向节 理 , 边 坡 剖 面 示 意 图见 图 图例 还至到片 麻岩 最终边坡 钻 孔 及 医疆编黔号标 高 匝习花 岗岩 图 边 坡剖面 示 意 图 边 坡地 质环 境 特 征 地 层 岩 性 采 用 爆 破 工 程 施 工 路 段 出 露 的地 层 多 以燕 山期侵 入岩 为主 , 代 表 岩性 为花 岗岩 、 闪长 岩 其 次 为太 古 界变质 岩 , 代 表 岩 性 为 片 麻 岩 、 混 合 岩 及 碎 裂状 大 理 岩 等 再 次 为第 四系 松 散堆 积 物 花 岗岩 多为强风 化似 斑状 粗 粒 结构 全 路 段 穿越 地 形 和 地 质 构造 条件 复 杂 , 有 包 括 褶 皱 、 大 断裂 和 各种 节 理裂 隙在 内的多种 地 质 构造 边 坡 岩 性 以碎 裂状 混 合 岩和 闪长纷 岩 为主 , 垂 直节 理较 为 发育 , 岩 石 硬 度较 大 该 路 段 主 要 沿 旧 路 进 行 开 收稿 日期 刁 一 林 士 炎 男 , 岁 , 硕 士 国家 自然 科 学基金 资助 项 目 岩 体 力学 参数 岩 体 力 学 参 数 根 据 边 坡 岩 样 物 理 力 学 性 质 实验 得 到 , 研 究 路 段 岩 性 较 为均 一 , 主 要 是 混 合 岩 , 多组 岩 样 力 学性 质 基 本 相 似 , 在模拟 计 算过 程 中采用 的物 理 力学 参 数 取 其 均 值 见表 表 岩 体物 理 力学参数 一 功 · 刁 · 一 , 口刀吐 曰 价 刀 爆 破震 动 数值模 拟 爆 破 数 值 分析模 型 的选 择 边 坡 开挖 数 值 模 拟 分 析 采 用 了 一 程 DOI :10.13374/j .issn1001-053x.2003.06.005
·508· 北京科技大学学报 2003年第6期 序,其状态方程为非线性结构材料的动力响应模 (1)应力在炮孔周围产生巨大的应力集中, 型,对于模型的处理采用有限元的方法.结合本 并出现了一定范围的拉应力区,拉应力值最大达 工程实际,所研究边坡为岩质边坡,其材料模型 到302MPa,远远大于工作边坡岩石的抗拉强度, 选取弹塑性模型. 这正是岩石爆破所要达到的效果,而在边坡岩体 2.2数值模型的建立 内部,由于岩体内节理裂隙的存在,应力波迅速 数值模拟计算范围包括了整个边坡,通过研 衰减,压应力和拉应力值都很小,未对边坡内部 究边坡中爆破地震波所产生的应力场、速度场、 的岩体造成实质性破坏, 加速度场以及位移场的综合分布情况,分析爆破 (2)坡脚位置质点震动的最大速度,水平方 地震波对整个边坡稳定性的影响程度.除了分析 向为20cm/s,垂直方向为44cms,大于边坡的其 在整个平面上的计算结果外,还对一些典型的位 他位置,因此坡脚是边坡受震动最容易破坏的区 置进行了具体分析,这些位置包括:(1)坡底面 域之一. 点,距离爆破点水平距离2.5m;(2)爆破台阶的 (3)通过分析波形图可以发现,坡底的振幅 坡脚点,距离爆破点4.6m;(3)坡中点;(4)坡顶 最大值出现在波形的第一个振幅上,然后衰减: 面上与坡面走向平行的方向上三个点.边坡单元 而坡顶的最大振幅出现在第二个振幅以后,然后 剖分见图2. 衰减,说明由于坡顶的放大效应造成了最大振幅 的相对滞后 (4)坡顶由于夹制作用小,振动时间持续长, 频率低:而坡底由于围岩的夹制作用,尽管振幅 大,但是频率高,因此,结合最大振幅和频率特 性,坡顶和坡底同样存在破坏危险的区域. (S)各点的位移值均比较小,最大为1.6mm, 因此不是引起破坏的主要因素, 3 模拟计算值与实测值比较 图2边坡计算模型单元划分示意图 数值模拟分析表明,坡底和坡顶处的速度和 Fig.2 Meshes of the slope-computing model 加速度是爆破敏感参量,所以在对爆破现场的实 为了研究爆破震动对已形成边坡的影响及 际监测中,对其进行了重点分析,爆破测试所用 爆破地震波在边坡岩体内的传播规律,进行了水 仪器为美国SAULS公司的NCSC5000型测震仪, 平和垂直坡面两个方向剖面的边坡模拟计算,计 测试参量为测点三个互相垂直方向上的最大速 算模型采用的边坡尺寸及爆破参数见表2.炸药 度分量及最大加速度、频率和位移,计算结果与 参数:p=1gcm',爆速D=3600m/s,C-J压力Pc-, 现场及经验数值相比较,模拟计算数值偏大,一 =1.03GPa,采用JWL状态方程. 般认为爆破产生的震动速度小于22cms时,爆 破震动对边坡不会产生影响,而计算结果值却要 表2模拟边坡尺寸及爆破参数 大许多(见图3). Table 2 The size and blasting parameters of the simulating 对这一差异,可以从以下两个方面来解释: slope (1)量级的差别.模拟计算时间精确度很高,计算 H/m W/m A/m ho/m 中所用的爆破响应时间仅为微秒级,而现场实测 6.0 1.15 1.45 1.15 的数据如爆破震动速度,爆破震动的响应时间最 φ/mm (ha+h)/m o/kg q/(kg-m) 32 5.6 3.45 0.35 小为毫秒级,大多是在秒级上,时间量级不一致: (2)测震距离的影响.由于采用有限元划分网格 23模拟计算结果分析 的局限,所建模型的大小一般取为一到几个抵抗 通过对边坡体内应力、速度、加速度和位移 线的大小,即1~7m:进行实测时,测震仪距起爆 场数值模拟结果的综合分析,可以得到以下主要 点相对较远,所测数据衰减量更大.另外,计算采 计算结果: 用的响应时间很短,计算得出的应该是冲击波的
北 京 科 技 大 学 学 报 年 第 期 序 , 其 状 态 方程 为非线 性 结构 材料 的动 力 响应 模 型 , 对 于 模 型 的处 理 采 用 有 限元 的方 法 结 合 本 工 程 实 际 , 所 研 究 边 坡 为 岩 质 边 坡 , 其 材 料 模 型 选 取 弹 塑 性 模 型 数值模 型 的建 立 数 值 模拟 计 算 范 围包 括 了整 个 边 坡 , 通 过研 究 边 坡 中爆 破 地 震 波 所 产 生 的应 力场 、 速 度 场 、 加速 度 场 以及 位 移 场 的综 合 分 布 情 况 , 分 析 爆 破 地 震 波 对 整 个 边坡稳 定性 的影 响程 度 除 了分 析 在整 个平 面 上 的计 算结果 外 , 还 对 一 些 典型 的位 置 进 行 了具 体 分 析 这 些 位 置 包 括 坡 底 面 点 , 距 离爆 破 点水 平 距 离 爆 破 台 阶 的 坡 脚 点 , 距 离爆 破 点 坡 中 点 坡 顶 面 上 与坡 面 走 向平 行 的方 向上 三 个 点 , 边 坡 单 元 剖分 见 图 千「 飞开 于 宁 于咔一 开于 于任 门勺广, ‘ 二 」仁二勺 必 」 刁灌刁刃 扁升干 、 砚声一」卜于 「 示 呈 厂勺「甲丫尸 丫州产咚一 卜 」 燕 热跳 瀚黯 嘟 抹 称挪遭 滤律幸溉 滋 粼日麟 尔牡陈 二 姗圈圈拼霍母栩 牙眯牛拼本下不要莽翔卜瘫三孚 鬓鬓鑫蓦乏 , 一 一 一 一 、 阳 攀 蓄卉 人幼因叮 宝 卜 ‘ 「攀 专 定 限二洲 止〔 庵 百 誉 , 愁 , 一】 厂戛万「厂 拿 圣 仁 」 旧 卜布 厂丁二口口 奎 「工二口 主 黯撇落井十一 了一 十 日到以 洲 、 一信一 不 厂飞二更 奎 厂 应 力 在 炮 孔 周 围产 生 巨大 的应 力集 中 , 并 出现 了一 定 范 围 的拉 应 力 区 , 拉应 力值 最 大达 到 , 远远 大 于 工 作边 坡 岩石 的抗 拉 强度 , 这 正 是岩 石 爆破 所 要 达 到 的效 果 而 在边 坡 岩 体 内部 , 由于 岩 体 内节 理 裂 隙 的存 牲 , 应 力波 迅速 衰 减 , 压 应 力 和 拉 应 力 值 都 很 小 , 未对 边 坡 内部 的岩 体 造 成 实质 性 破 坏 坡 脚 位 置 质 点震 动 的最 大 速 度 , 水 平方 向为 耐 , 垂 直 方 向为 而 , 大 于 边 坡 的其 他位 置 , 因此 坡脚 是边坡 受 震 动 最 容 易破 坏 的区 域 之 一 通 过 分 析 波 形 图可 以发 现 , 坡 底 的振 幅 最 大值 出现 在波 形 的第 一 个 振 幅 上 , 然 后 衰减 而 坡 顶 的最 大振 幅 出现 在 第 二 个 振 幅 以后 , 然 后 衰 减 , 说 明 由于 坡 顶 的放 大 效应 造 成 了最 大振 幅 的相 对 滞 后 坡 顶 由于 夹 制 作 用 小 , 振动 时 间持续 长 , 频 率 低 而 坡 底 由于 围岩 的夹 制 作用 , 尽 管 振 幅 大 , 但 是 频 率 高 因此 , 结 合 最 大振 幅 和 频 率特 性 , 坡 顶 和 坡底 同样 存 在 破 坏 危 险 的 区 域 各 点 的位 移 值均 比较 小 , 最 大 为 , 因此 不 是 引起 破 坏 的主 要 因 素 图 边 坡 计 算模型 单 元 划 分 示 意 图 一 为 了研 究 爆 破 震 动 对 已 形 成 边 坡 的 影 响 及 爆破 地 震 波在 边坡 岩 体 内的传 播规 律 , 进 行 了水 平 和 垂 直 坡 面 两 个 方 向剖面 的边坡 模 拟 计 算 , 计 算模 型采 用 的边 坡 尺 寸及 爆 破 参 数 见表 炸 药 参 数 夕 留 ,, 爆 速 , 一 压 力 一 , 采 用 状 态 方 程 · 表 模 拟边坡 尺 寸 及 爆 破参数 山 牙 。 沪 叮 一 , 模 拟 计 算 结 果 分 析 通过 对 边坡 体 内应 力 、 速 度 、 加速 度 和位 移 场 数 值 模 拟 结 果 的综 合 分 析 , 可 以得 到 以下 主 要 计 算 结 果 模 拟 计 算值 与 实测 值 比较 数 值 模 拟 分 析 表 明 , 坡 底 和 坡 顶 处 的速 度 和 加速度 是爆破敏感 参 量 , 所 以在 对爆破现场 的实 际 监 测 中 , 对 其 进 行 了重 点 分 析 爆 破 测试 所 用 仪 器 为 美 国 公 司 的 型 测 震 仪 , 测 试 参 量 为 测 点 三 个 互 相 垂 直 方 向上 的最 大 速 度 分 量 及 最 大 加 速 度 、 频 率 和 位 移 计 算 结 果 与 现 场 及 经 验 数值‘ 相 比较 , 模 拟 计 算 数值偏 大 , 一 般 认 为爆破 产 生 的震 动 速 度 小 于 耐 时 ‘ , 爆 破 震 动 对 边坡 不会 产 生 影 响 , 而 计 算 结 果值 却 要 大许 多 见 图 对 这 一 差 异 , 可 以从 以下 两 个 方 面 来 解 释 量 级 的差 别 模 拟 计 算 时 间精确 度 很 高 , 计 算 中所 用 的爆破 响应 时 间仅 为微秒 级 , 而 现 场 实测 的数据 如 爆破 震 动 速度 , 爆破 震 动 的响应 时 间最 小 为毫 秒 级 , 大 多是在 秒 级 上 , 时 间量 级 不 一致 测 震 距 离 的影 响 由于 采用 有 限元 划 分 网格 的局 限 , 所 建 模 型 的大 小一 般 取 为一 到几 个 抵 抗 线 的大 小 , 即 一 进 行 实 测 时 , 测 震 仪 距 起 爆 点相对 较 远 , 所 测 数 据 衰 减 量 更 大 另外 , 计算 采 用 的 响应 时 间很 短 , 计 算得 出的应 该 是 冲 击 波 的
Vol.25 No.6 林士炎等:爆破震动对高速路边坡影响的数值模拟 ·509* 200 (a) 12.5-25cm/s,而实际监测的震速则为21.12cm/s, 150 v=275140261 可以说数值模拟计算的结果与实测结果很接近, 邑100 从而印证了数值模拟计算的可行性, 50 4结论 4 (1)爆炸开始后极短的时间内,在炮孔内壁 0 500 1000 及邻近范围内产生了较大的应力集中,其值远大 t/μs 200 于岩石强度,造成岩石破碎,爆破产生的应力波 (b) 在岩体内部的传播过程中,应力极值逐步衰减. 150F v=7×1021o7s (2)坡脚位置质点震动的最大速度大于边坡 怎100 的其他位置,是边坡最易破坏的区域之一.坡顶 位置的震动持续时间长,频率低,也是边坡容易 50F 产生破坏的区域.距离边坡越近,破坏的可能性 越大.在实际爆破施工过程中,应重点对坡脚处 250 500 750 1000 和坡顶位置进行爆破现场监测,保证其质点震速 t/us 在允许的范围内, 图3坡脚及一定距离处质点震速与时间关系 Fig.3 Relations of particle vibration velocity to time in the 参考文献 slope foot and certain distance 1周晓杰,黄海,王振中,等.爆炸震动效应对公路边 效应:而实测时得到的则是应力波的效应,甚至 坡的影响[)西部探矿工程,2001(⑤):73 是应力波经过衰减以后的效应,由于波的传播在 2王德志,谭文辉,乔兰,高速路连接线路堑开挖中爆 时间和空间上的衰减很明显,经过衰减以后的波 破震动监测与评价).工程爆破,2002(1):10 3王在泉,陆文星。高边坡爆破开挖振动传播规律及 的效应会有很大的改变,因此导致了模拟计算的 质量控制[).爆破,1994,11(3):2 数值偏大.基于以上考虑,为了使计算更符合实 4徐顺香,临长高速公路路堑开挖中的控制爆破技术 际,对计算结果进行了折减,通过对大量计算结 .爆破,2002.19(1):6 果的分析与比较以及参考相关资料,确定折减 5陈寿如,宋光明。确定露天边坡爆破度动临界振速 系数为0.25较为合适. 的-种方法[.矿业研究与开发,2001,21(4):9 在已形成的边坡上,t=0.38ms时的质点震动 6中国力学学会工程燥破专业委员会,爆破T程M 速度最大,为50-100cm/s,经折减后,速度值为 北京:治金工业出版社,1992 Numerical Simulation on the Influence of Blasting Vibration on the Freeway Slope LIN Shiyan",LI Changhong",QIAO Lan",LIU Baoxu 1)Construction and Environment Protection Department,Government of Hebei Province,Shijiazhuang 050053,China 2)Civil and Environment Engineering School,University of Science and Technology,Beijing.Beijing 100083.China ABSTRACT Taking the cutting excavation of Beijing-Qinhuangdao freeway slope as an example,the dynamic features and the distribution regularity of blasting vibration waves such as stress field,displacement field and vel- ocity field were analyzed and various aspects of the influence of blasting vibration on the mountainous freeway-cut- ting slope were studied by use of the numerical simulation method.The comparison of simulation results with in- situ measured results shows that they are approximate. KEY WORDS cutting excavation;blasting vibration;rocky slope;numerical simulation;stability assessment
林 士 炎等 爆 破 震 动 对 高速 路 边 坡 影 响 的 数值 模拟 一 一 耐 , 而 实 际监 测 的震 速 则 为 耐 , 可 以说数 值模 拟 计算 的结 果 与实测 结果 很接近 , 从 而 印证 了数 值 模拟 计 算 的可行 性 哪 八“ 甲澎日。, 一 一‘ 一 山 一一一一一 协 图 坡脚 及 一 定距 离处 质 点 艇 速 与时间关 系 · 效应 而 实测 时得 到 的则 是应 力 波 的效应 , 甚至 是应 力波经 过 衰减 以后 的效应 由于波 的传播 在 时 间和 空 间上 的衰减 很 明显 , 经 过 衰减 以后 的波 的效应会 有 很 大 的改变 因此 导 致 了模拟 计 算 的 数值偏 大 基 于 以上 考 虑 , 为 了使 计 算 更 符 合 实 际 , 对 计 算 结 果进 行 了折 减 , 通 过 对 大 量 计 算 结 果 的分 析与 比较 以及 参考 相 关 资料 〔, , ,, 确 定折 减 系数 为 较 为 合适 在 己 形 成 的边 坡 上 , 时 的质 点震 动 速 度 最 大 , 为 一 而 , 经 折 减 后 , 速 度 值 为 结 论 爆 炸 开 始 后 极 短 的 时 间 内 , 在 炮 孔 内壁 及邻 近 范 围 内产 生 了较 大 的应 力集 中 , 其值 远 大 于 岩 石 强度 , 造 成 岩 石 破 碎 爆 破 产 生 的应 力 波 在 岩 体 内部 的传 播 过 程 中 , 应 力 极 值 逐 步 衰减 坡 脚 位 置 质 点震 动 的最 大 速 度 大 于 边 坡 的其 他 位 置 , 是 边 坡 最 易破 坏 的 区 域 之 一 坡 顶 位 置 的震 动 持 续 时 间长 , 频 率低 , 也 是 边坡 容 易 产 生 破 坏 的 区 域 距 离边 坡越 近 , 破 坏 的可 能性 越 大 在 实 际爆 破 施 工 过 程 中 , 应 重 点对 坡脚 处 和 坡 顶位 置进 行 爆 破 现 场 监 测 , 保 证 其 质 点震速 在 允 许 的范 围 内 参 考 文 献 周 晓杰 , 黄 海 , 王 振 中 , 等 爆 炸 震 动 效应 对 公路 边 坡 的影 响 西 部探矿 工 程 , 王 德 志 , 谭 文 辉 , 乔 兰 高速 路连 接 线 路 堑 开 挖 中爆 破 震 动 监 测 与 评 价 工 程爆 破 , 王 在泉 , 陆 文 星 高边坡 爆破 开 挖 振动 传播规 律及 质 量 控制 爆 破 , , 徐顺 香 临长 高速 公路路 堑 开 挖 中 的控 制 爆破 技术 爆破 , , 陈寿 如 , 宋 光 明 确 定 露 天 边 坡爆 破 震 动 临界振速 的 一 种 方 法 【 矿 业 研 究 与 开 发 , , 中国力学 学会 工 程 爆 破 专 业 委 员会 爆破 工 程 北 京 冶 金 工 业 出版 社 , , “ 沙 口 刀 , 例 刀 , 如 城 , 石山 , , 既 , , “ 一 , 叮 , 一 材 叩