FLAC数值模拟分析学习计划.docxVIP

FLAC数值模拟分析学习计划 FLAC数值模拟分析学习计划 PAGE PAGE5 FLAC数值模拟分析学习计划 PAGE 工程大要 研究此段为中条山地道 K9+450～K10+560段，此处地道最大埋深约 540m，主要由太古界涑水 群表壳岩组合解州片麻岩（ Hgn）地层组成，组成中条山地道分水岭北侧的主体；此段组成地道的 围岩岩性主要为变粒岩、花岗片麻岩等。该套地层岩性复杂，组合无规律。岩层产状整体倾向南东， ｏ ｏ 倾角一般在50～70间变化。在AK9+900～AK10+000段为地域性断层影响段，此断层为破碎岩石 组成，将为基岩裂隙水下渗供应通道， 地道开挖必然引起涌水，同时此段围岩稍差，施工时易引起 3D 在此处，详尽运用FLAC进行模拟的区段均取洞身 YK10+100～YK10+180段。该区段为V级围 岩地域，埋深为505～512m，为断层，周边太古界涑水群花岗片麻岩、黑云斜长片麻岩、岩石破碎。地道断面为SVc型，如图2-2所示。 图2-2SVc型地道断面图 地道衬砌按新奥法原理设计，采用SVc型复合式衬砌，该衬砌适用于地道洞身 V级断层影响带 及柔弱破碎围岩段的初期支护及衬砌，超前支护各环采用 42×4mm注浆小导管超前预加固围岩， 长，环向间距35cm，搭接长度，斜插角10ｏ～15ｏ，每环37根；初衬以喷、锚、网为主要支护手段：钢拱架为I20a型钢，纵向间距75cm，每榀钢拱架之间采用φ22钢筋连接，环形间距 ；锚杆采用D25中空注浆锚杆，长，间距75cm（纵）×100cm（环），与钢拱架交叉部署；喷C25早强混凝土26cm。二次衬砌和仰拱均为C30钢筋混凝土结构，厚50cm。 数值计算模型 依照中条山地道工程的本质状况，为提高计算速度，在保证计算精度的前提下，取桩号 3D YK10+100～YK10+160段采用大型有限差分软件 FLAC进行建模解析。 对于全断面法、预留中心土法、台阶法，由于整个地道模型左右对称，为减少计算量，可取隧 道模型的一半计算，地道的计算模型 I如图2-3所示。 图2-3计算模型 I表示图 x、y、z各方向的长度分别为 60m、60m和140m。模型的表面为自由界线，底部为固定界线 条件，四周为法向位移拘束界线条件。模型共有 10912个节点，9360个实体单元。 对于单侧壁导坑法，由于整个地道模型左右不对称， 可取地道模型的一半计算， 地道的计算模 型II如图2-4所示。 x、y、z各方向的长度分别为 120m、60m和140m。模型的表面为自由界线，底部为固定界线 条件，四周为法向位移拘束界线条件。模型共有 25637个节点，23220个实体单元。 图2-4计算模型 II表示图 3D FLAC供应了丰富的单元库供用户选择， 这样能够对各种资料进行模拟。计算中V级围岩均采 用摩尔—库伦资料，初期支护系统中的喷层混凝土视为线弹性体，二次衬砌混凝土及仰拱均采用-本源网络，仅供个人学习参照 3D FLAC软件内置的壳单元（shell）进行模拟，系统锚杆和锁脚锚杆均分别采用锚索单元（cable）和桩单元（pile）进行模拟，超前导管采用桩单元（pile）进行模拟。各资料参数见表2-1。 数值计算说明 需要说明的是，由于研究该段埋深很大，最大处达到540m。此处假设地道埋深为500m，假设地应力场仅考虑自重应力场，则依照等效方法在模型顶部竖直方向应施加的垂直向应力为 z=2650××。地应力平衡后，模型底部垂直向应力 z=2650×× ，地道中心垂直向应力 z=2650××现依照上述方法进行平衡，得 到的竖向应力云图如图 2-5及图2-6所示。 图2-5 模型Ⅰ平衡地应力后垂直向应力分布云图（单位： Pa） 图2-6模型Ⅱ平衡地应力后垂直向应力分布云图（单位： Pa） 地道开挖进尺为2m。纵向每个网格长为 2m。纵向有30个网格，即不论哪一种工法，开挖到30 步时，地道模型即贯穿。 表2-1资料参数汇总表 资料 密度 弹性模量 泊松比 黏聚力（MPa） 内摩擦角 （kg/m3） （GPa） (°) V级围岩 2650 5 38 初衬 —— —— 二次衬砌 2500 30 —— —— 系统锚杆 2479 —— —— —— 锁脚锚杆 2549 —— —— —— 超前导管 2549 —— —— —— 对于初期支护中钢拱架和钢筋网的支护作用采用等效方法计算，立刻钢拱架和钢筋网的弹性模量折算给混凝土，其计算方法为： E E0 AgEg（2-1） Ac 式中： E——折算后的混凝土弹性模量（ Pa）； E0——原混凝土的弹性模量（ Pa）； Ag——初期支护钢拱架的截面积（ m2）； Eg——初期支护钢拱架的弹性模量（ Pa）； Ac——混凝土的截面积（m2）。 对于系统锚杆及锁脚锚杆的参数采用，亦同样