基于GeoStudio下尾礦坝地应力场的数值计算.doc

基于GeoStudio下尾矿坝地应力场的数值计算尾矿坝应力场410.20 文献标志码:A 0引言 岩土工程的数值分析方法比较多，主要有：有限元法、边界元法、离散元法及拉格朗日元法、DDA法等。其中，有限元法是目前已广应用于岩土工程与结构分析的有力工具。该法是把一个实际的结构物或连续体用一种由多个彼此相联系的单元体所组成的近似等价物理模型来代替，通过结构及连续体力学的基本原理及单元的物理特性建立起表征力和位移关系的方程组，解方程组求其基本未知物理量，并由此求得各单元的应力、应变以及其他辅助量值。尾矿坝的数值模拟计算 （）计算模型 计算模型的几何尺寸及边界条件的确定不但涉及到计算量的问题（即单元大小及数量），而且影响到计算的精度及整个坝体的力学分析，因此，合理确定计算剖面、剖面的尺寸及边界条件非常重要。结合现场实际情况，确定出模型的几何尺寸。 边界条件的确定也有很多种，根据以往所做类似工程的有限元数值分析实例，再结合现场实际情况而定。具体如图，其中U=0，V0，底部边界为水平、垂直均约束U=0，V=0。  （a）计算模型 （b）实际剖面 图（）材料的本构模型与参数 材料模型选择非线性材料模型能更好地模拟土工问题。一般弹性非线性模型是把弹性模量和泊松比看作随应力状态变化而改变的变量，即和是应力状态的函数。对土体而言应用较多的模型为邓肯（Duncan）和张（Chang）提出了的双曲线模型。这种模型已被广泛应用于土坝和地基等土工问题的有限元分析，在很多情况下可以得到满意的结果。 表 材料参数一览表 天然重度　 静力 水平 有效应力 弹性模量 泊松比 Kv （10-4cm/s） Kh （10-4cm/s） γ （KN/m3） cˊ （KPa） φˊ （°） E (MPa) μ R　 ①1 2.5 2.61 17.5 14.4 21.6 55.0 0.3 0.62 尾粉土①2 2.5 2.61 19.2 14.4 21.6 55.0 0.3 0.60 尾粉土①3 2.5 2.61 19.5 14.4 23.2 55.0 0.3 0.57 尾粉土①4 2.5 2.61 19.5 16.0 24.8 55.0 0.3 0.58 尾粉质黏土①4-2 1.5 1.5 18.5 12.0 20.0 45.0 0.35 0.58 采矿弃渣 （未碾压） 12.0 12.0 20.5 20.8 28.0 150.0 0.18 0.70 初期坝体 12.0 12.0 21.0 15.0 32.0 150.0 0.18 0.70 粉质黏土② 0.1 0.1 18.5 28.0 25.0 25.0 0.2 0.70 泥质粉砂岩③ 0.1 0.1 26.0 100.0 38.0 2000.0 0.15 0.70 材料的物理力学参数是模拟计算的基础资料。根据该尾矿的物理力学实验测试结果，参照国内外有关尾矿和堆石坝等相关资料，综合分析后，确定出各材料的计算参数，见表.1所示。数值模拟计算结果分析计算模型单元网格划分如图2。 （a）+1924m标高尾矿坝的计算网格（现状） （b）+1934m标高的尾矿坝计算网格（最终标高） 图2 计算模型单元网格划分 在计算中考虑到流固耦合，即地下水的作用，为此计算了两种情况，一是洪水情况（干滩面为0m，即规范中的最小长度），二是正常情况（按照10m）分别进行计算。计算结果分别如图～图。 图16 加高后（+1934m）位移矢量图 图17 加高后（+1934m）x方向应力云图（单位KPa） 图18 加高后（+1934m）y方向应力云图（单位KPa） 图19 加高后（+1934m）剪应力云图（单位KPa） 从计算结果（图）中可以看出： ①尾矿坝体的应力、应变等分布规律与其他资料获取尾矿坝的应力应变等规律基本相同。 ②坝体，0.179m，，最高点（即坝顶处）达到m左右，沉降。 ③坝体，。水平方向的变形位移不大垂直方向的变形位移较大沉降GeoStudio软件在岩土工程中的应用 56 1 2 济南大学学报(自然科学版) 第22卷 U=0，V=0 U=0 V≠0 +1924m U=0，V=0 U=0 V≠0 +1934m 尾粉土 1:5.0 基底（基岩） 初期坝 +1934m +1924m +1850m 尾粉质粘土