古瓦地下厂房弹塑性剖析.docVIP

古瓦地下厂房弹塑性剖析 　　摘要：基于三维非线性有限元法通过采用ADINA有限元软件建立弹塑性材料计算模型，对硕曲河古瓦水电站地下厂房（包括主副厂房、主变兼尾闸室、母线洞、交通洞、尾水连接洞以及尾水洞等结构）开挖施工进行数值模拟分析。研究了整个开挖过程中实施锚杆、喷混凝土等支护措施后围岩洞面的位移变形状态，应力状态，塑性区分布以及洞室群的整体稳定性。计算结果表明：洞室群整体变形不大，分布均匀连续，没有突变；围岩采取支护措施后应力基本都在许可范围之内；塑性区主要集中在母线洞与主厂房和主变室的交叉处，但总体积不大，与洞室开挖规模相比较小。地下厂房洞室群整体稳定，但是屈服区集中部位需要采取其他措施。 　　关键词：地下厂房；应力应变；塑性区；围岩稳定 　　中图分类号：TU751.2文献标志码：A文章编号2015）002006003 　　1厂区厂房布置及地质条件 　　古瓦水电站为硕曲河规划的“龙头”水库电站，具有年调节性能，工程规模中Ⅱ型。古瓦水电站采用混合式开发，利用落差约320 m，调节库容2641亿m3，装机容量222万kW，安装3台混流式机组。 　　[JP2]地下厂房 （包括副厂房和安装间）长6892 m，宽186 m，高4309 m，开挖高程3 0672 m，主机间尺寸4810 m×186 m×4309 m（长×宽×高），主机间内安装三台单机容量为74 MW的立轴混流式水轮发电机组，机组安装高程3 07580 m。主机间从下到上分四层布置，第一层为球阀层，层面高程3 07250 m；第二层为水轮机层，层面高程为3 07780 m；第三层为电缆层，层面高程3 08240 m；第四层为发电机层，层面高程3 08680 m；安装间长2082 m，宽186 m；副厂房长1302 m，[JP3]宽186m。主变室长5812 m，宽145 m，高167 m，底板高程3 08680 m；尾水闸门室开挖尺寸长5096 m，宽50 m，高763 m，底板高程3 08920 m。[JP]主变室开挖尺寸5812 m×145 m×167 m（长×宽×高），共分2层，第1层高程3 08680 m，高度850 m，第2层为电缆室，楼板高程3 09530 m，拱顶处净高550 m。[JP]主厂房与主变室之间设3条母线道，断面尺寸550 m×540 m（宽×高）。尾水洞连接洞尾端采用208%的反坡与尾水隧洞连接，尾水隧洞全长约10863 m，宽75 m，出口与硕曲河斜接；进厂交通洞自安装间左侧进厂，洞口设回车场，地坪高程3 09400 m。 　　2数值仿真分析 　　2.1采用ADINA 软件中的单元生死算法特点 　　采用有限元软件模拟岩土和地下空间的施工时，涉及到结构材料增加或减少的过程均需要频繁采用单元生死（Element Birth/Death）功能。如单元死（或Inactive）的功能常常用来模拟隧道中开挖的过程，而生（Active）则指开挖过程中围岩的支护过程，如初衬支护、二衬支护等。在单元生死的过程中，ADINA提供丰富的特征可更好处理各种复杂工程问题。ADINA定义的具有单元生死的重复单元中可以定义锚杆、土钉等单元，此功能主要用来模拟带有锚固、预应力锚固的支护结构。 　　2.2非线性求解过程 　　本文涉及的非线性问题属于材料非线性问题，进行非线性分析时采用弹塑性材料模型。弹塑性材料模型允许有能量耗散，塑性功是不可逆的。在程序中根据增量应变预测得到每一个积分点的中间增量步状态。对每步的第一个循环的预测是基于前一步的应变增量。在迭代过程中，从前一个迭代得到的应变，作为估算应变来进行刚度矩阵计算。若残余应力满足选择的容差，视为收敛将不再迭代。 　　3计算模型 　　3.1计算域及计算网格 　　本次古瓦水电站地下厂房洞室群开挖支护仿真分析中，模拟的范围为：上下游方向为300 m，为厂房上游边墙到主变室下游边墙距离5300 m的5倍；沿厂房轴线方向为280 m，相当于从安装间到副厂房8190 m总长度的4倍；底面取为2 93240 m高程，距厂房最低位置（高程3 06721 m）约135 m，主厂房拱顶高程3 11233 m，计算模型的顶部到地面，上部岩层厚度100～140 m。 　　3.2计算采用的坐标系 　　X轴：沿厂房横轴，顺水流向，从上游指向下游。Y轴：沿厂房纵轴，从副厂房水平指向安装间。Z轴：铅直向上，满足笛卡儿右手坐标关系。 　　3.3计算模型 　　本电站地下厂房开挖分析模型对主副厂房、主变兼尾闸室、母线洞、交通洞、尾水连接洞以及尾水洞等结构进行了模拟，支护分析模型模拟了锚杆、喷混凝土。计算中围岩采用四面体实体网格单元、锚杆采用三维杆单元、衬砌采用壳单元，整个计算域共划分为175 132个