预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共30页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

施工爆破地振动波传播规律的数值模拟分析 1工程概况 龙头山隧道左右线进口最小净距23m,洞身左右线最大净距51m,出口最小净距20.8m。 龙头山隧道采用光面爆破开挖,根据不同围岩类别采用不同的炮眼布置和不同的装药量。 Ⅳ、Ⅴ类围岩采用台阶法开挖,掏槽眼3.5m,其它眼深度3.3m,预计进尺3.0m。采用非电毫秒雷管及2号岩石硝铵炸药。Ⅳ、Ⅴ类围岩采用四臂台车钻孔,孔径Φ43mm,周边眼采用Φ25mm药卷间隔装药。Ⅳ、Ⅴ类围岩爆破孔布置如图1所示,掏槽眼布置图如图2所示。(一)部掏槽眼采用连续装药结构,布设10个孔,平均每孔装药量2.34kg,小计装药量总合23.4kg。起爆雷管段别为1、3。 图1Ⅳ、Ⅴ类围岩开挖炮眼布置图 图2掏槽眼布置图 掏槽眼段采用2号岩石硝铵炸药,其爆速为3000m/s,密度为1g/cm3。计算输入参数如表1所示。 完整岩石的力学参数表1 物性参数 围岩弹性模量E (GPa)泊松比 (ν)内摩擦角φ()粘聚力c (KPa)重度γ (KN/m³)纵波速度cρ(m/s)Ⅱ类围岩1.00.3424.582116.3850~1450Ⅲ类围岩1.30.339150202500~3200Ⅳ类围岩6.00.339700233400~3800Ⅴ类围岩200.24501500253900~4600 2计算模型的建立 根据龙头山地质剖面图,截取一定范围纵向160m,横向近500m,高度模拟一定坡度建立如图3所示计算模型图。模型共划分节点,单元95158个。其中对已经施工完成的结构用板单元模拟,定义成弹性材料属性。围岩采用实体单元建模,并且定义材料属性为摩尔-库仑。边界采用曲面弹簧单元模拟。图4为隧道网格图,图5为隧道内部结构图,图6为油库位置图,图7为模型网格图。计算采用MIDAS-GTS有限元分析软件进行。 对于爆破振动引起的冲击荷载,采用时程函数来模拟,并转化成作用到孔壁上的孔壁压力,图8给出了计算过程中施加在爆破面上的面压力。它是一个时程函数。其荷载衰减形式如图9所示。 模型模拟情况为:沿隧道开挖方向取Ⅱ类围岩60m,Ⅲ类围岩40m,Ⅳ类围岩40m,Ⅴ类围岩20m。其中左右线二衬施作完成10m,拆除临时支撑但尚未施作二衬段30m,临时支撑尚未拆除段30m,30m仅开挖了两侧壁上台阶导洞段30m,即爆破面距离临时支撑最近位置为30m。计算模拟掏槽爆破点在左线中导洞开挖里程ZK5+940位置。计算掏槽眼一次最大装药量23.4kg爆破,分析比较对左右线已完成结构、已开挖洞室初期支护及临时支撑以及距离隧道300多米远处油库的影响。评价指标主要选取各测点速度是否满足爆破安全规程。 采用有限元法分析爆破震动影响的一项关键工作是建立爆破加载模型,包括确定爆破激振力的大小、作用位置和方向、峰值时刻和持续时间等方面的内容。 本报告根据计算和爆破的实际情况,在不失一般性的条件下作了以下的假设: 爆破荷载以压力形式的均布荷载作用在隧道壁上,方向垂直于洞壁。根据计算情况,输入爆破荷载曲线为脉冲形式。如图9所示。 为了解爆破振动波在岩体中的传播规律,取计算时间为2s。 Ⅲ类围岩 Ⅱ类围岩 Ⅳ类围岩 Ⅴ类围岩 图3爆破计算模型 Ⅱ-Ⅱ Ⅲ-Ⅲ 右线未拆除临时中隔壁段 Ⅰ-Ⅰ 右线拆除临时中隔壁段 右线施作二衬段 图4隧道网格图 临时中隔壁 图5隧道内部结构图 油库 图6油库位置图 图7模型网格图 图8掏槽眼孔壁面上面压力示意图 3边界条件的定义 3.1特征值分析 为了进行特征值分析通过弹性边界来定义支座的边界条件。计算通过曲面弹簧定义弹性边界,弹性系数根据道路设计规范的地基反力系数计算。 竖直地基反力系数:(KN/m³) 水平地基反力系数:(KN/m³) 注: 模型各方向截面积如表2所示: 模型各方向截面积表2 围岩 面积(m2) Ⅱ类围岩 Ⅲ类围岩 Ⅳ类围岩 Ⅴ类围岩 Ax8352609763863241.94638309230921546Ay66853.9--64501.4Az3125420836208310418 计算竖直及水平地基反力系数如表3所示: 竖直及水平地基反力系数表3 围岩 K(KN/m³) Ⅱ类围岩 Ⅲ类围岩 Ⅳ类围岩 Ⅴ类围岩 Kx45713.8979266872.0236303326.31841303766.52756996.1679186262.73661398135.70741721060.445Ky20955.45818--424777.2394Kz27869.5581542180.10513194677.4083841551.2116 特征值分析得到第1,2振型周期值:1.416644s,0.841708s。 3.2时程分析 对于动力分析建立一般的边界条件会由