《浅析拉格朗日有限元数值计算法(精品)》由会员分享,可在线阅读,更多相关《浅析拉格朗日有限元数值计算法(精品)(5页珍藏版)》请在金锄头文库上搜索。
1、成 绩 中 国 矿 业 大 学2011级硕士研究生课程考试试卷考试科目 岩土工程数值计算法 考试时间 2011.11.27 学生姓名 夏明 学 号 ZS11020068 所在院系 矿业工程学院 任课教师 徐志伟 中国矿业大学研究生院培养管理处印制对论文基于拉格朗日差分法的露天边坡稳定性研究中数值计算法浅析基于拉格朗日差分法的露天边坡稳定性研究一文以易门铜矿露天开采境界优化方案下的边坡为工程背景,通过岩体构造调查,质量分类,室内力学试验,力学参数的工程处理。建立铜厂露天开采边坡的三维地质模型,采用DIMINE数字矿山软件的耦合集成技术四面体网格化将地质模型转换为力学模型,应用基于拉格朗日法的有限
2、差分(FLAC3D)大变形方法对铜厂露天开采边坡的稳定性进行了数值模拟,分析了基于强度折减理论计算出的边坡安全系数以及基于莫尔库伦屈服准则的边坡开采后的位移、应力等的变化状况,得出了易门铜厂露天矿露天境界优化方案下的边坡的稳定性状况。岩质边坡稳定性评价的方法分主要有:极限平衡法、数值模拟计算、地质力学物理模拟试验和其它新方法。随着计算机技术和计算方法的发展,复杂的工程问题可以采用离散化的数值计算方法并借助计算机得到满足工程要求的数值解,数值模拟技术是现代工程学形成和发展的重要动力之一。通过计算模拟,可以模拟并得到模拟体内部的应力应变关系,再现其变形甚至破坏过程及其机制。在岩土工程数值分析中最常
3、用的数值方法有有限元法、离散元法、边界元法等。拉格朗日差分法(FLAC法)源于流体力学。它首先是Cundail在80年代提出来的,其基本原理类似于离散元法,但它却能像有限元那样适用于多种材料模式与边界条件的非规则区域的连续问题求解。在求解的过程中,FLAC又采用了离散元的动态松弛法,不需要求解大型联立方程,便于在微机上实现。另一方面,同以往的差分分析相比,FLAC在以下几个方面做了较大的改进和发展:它不但能处理一般的大变形问题,而且能模拟岩体沿某一弱面产生的滑移变形。一般有限单元法可以用来解决材料非线性问题,但对于大变形的几何非线性问题,有限单元法和边界元法都无能为力。拉格朗日法是分析非线性大
4、变形问题的数值方法,它依然遵循连续介质的假设,基于拖带坐标系的基本原理。用差分法或按时步显式迭代求解,不但可以解决几何非线性,也能解决材料非线性问题。下面是原文的部分描述:一、露天矿边坡三维地质模型的建立1.1计算几何模型范围及模拟方案设计 a模拟计算方案设计 决定边坡稳定性的因素很多,其中包括岩体强度,结构特征,水理性质和边坡几何参数等。根据现场工程地质调查,力学试验参数,对露天境界边坡进行修改,修改方案同第五章简化实体模型的建立。 b模型范围 模型范围需考虑两方面的内容:首先,必须能够覆盖铜厂露天矿边坡岩体的主要特征,即能容纳前面建立的实体模型;另外,为配合数值计算的需要,要考虑边坡岩体稳
5、定影响范围以及数值模拟软件。以南北向为Y方向,东西为x向,高程方向为z向。确定的数值计算范围参数见表7-1表7-1模型尺寸参数表c单元划分及计算几何模型 单元划分需要考虑三方面的问题。首先,能体现铜厂的地质特征,在岩性,断层和地表高度变化的位置需要细分单元块,以确保计算模型能够真实的反映岩体的几何形态;其次,保证数值计算中应力、位移等计算结果的连续性;第三,由于计算机计算能力和内存大小的限制,单元大小必须能保证计算程序能正常的运行。根据以上原则,对四面体单元单元块边长以及体积进行限制,四面体单元最小半径边比设为1. 414,落雪组四面体单元块最大体积设为500,火成胶结角砾岩组四面体单元块最大
6、体积设为10000 已经建立了边坡地质模型,根据地学实体模型,按照模型范围和单元划分规则,使用DIMINE软件四面体化功能,完成了从地质模型到数值计算模型的转换。网格的划分由先前地学软件DIMINE确定,网格的大小为245762格单元,42840个节点。 从以上两图中可以看出,通过转换程序的转换,得到了能完全反映地质结构以及岩性在空间上分布的精确三维力学模型,实现了三维地学模拟与数值模拟的耦合。从图7-4中可以看出,建立的计算模型能较准确地反映地形的微地貌形态以及能充分体现岩体、断层的尺寸和在空间中的位置,露天坑的坡度也在数值计算模型中也得以显现。1.2计算参数与约束条件a.屈服准则的选择 屈
7、服准则是判断材料进入塑性阶段的准则。岩土材料的屈服准则经过几十年的研究,提出的表达式不下几十种。以莫尔一库伦定律为基础的摩擦屈服准则在岩石力学与工程的实践中经受了考验,至今仍被广泛的应用。FLAC3D所采用的屈服准则之一为莫尔库伦与拉破坏准则结合的复合准则。b.计算力学参数 根据现场地质调查和岩石力学试验结果可确定岩体力学参数,由于FLAC3D中采用体积模量和剪切模量描述弹性模量和泊松比,所以根据式7-3 5式计算体积模量和剪切模量。其中:K体积模量 G剪切模量 E弹性模量 v泊松比考虑到计算机计算能力和内存大小的限制,在建立块段模型时将落雪组地层中的三种岩性(泥质白云岩,灰白色白云岩,青灰色
8、白云岩)合并成一种岩性,这样可以减少各种岩性间边界上的单元块次分,从而做到减少单元块数目。同时由第四章力学试验可得,落雪组地层的三种岩性力学参数很相似,所以在力学强度上将其合并也是可行的。计算中所需的力学参数如表7-2所示。表7-2 各类岩石物理力学参数c.边界约束模型边界约束采用位移(在FLAC3D中实质上是速度约束)约束的方式。底部所有节点取X、 Y、 Z三个方向的固定约束;南北边界为Y方向约束,东西边界进行X方向约束(即各面的法向方向的约束)。坡面为自由面。二、计算结果与分析按前述约束条件,在只考虑重力作用的情况下,进行本构模型为Mohr-Coulomb模型的弹塑性求解,直至系统达到平衡
9、。图7-7为数值计算过程中弹塑性求解阶段的系统不平衡力演化全过程曲线。图7一7系统不平衡力演化全过程曲线 设定体系最大不平衡力与典型内力比值下限为lxlo一5,迭代计算3683步后,系统达到近似平衡(见图7一7)。为便于分析,根据边坡的几何形态变化沿Y轴方分别在Y=2734650m、Y=2734750m、Y=2734850m和Y=2735000m切出4个典型剖面。下面就数值分析计算结果从安全系数、变形(位移)情况、应力情况等三个方面来分析边坡岩土体的力学响应特性,以及可能的内在变形破坏机理。浅析论文中的数值计算方法作者利用FLAC3D数值模拟软件对铜厂露天境界优化方案下的边坡稳定问题进行了数值
10、模拟分析,并对开挖后的边坡安全系数,位移、应力的变化状况进行了分析。这为寻找露天矿坑高边坡的稳定性与潜在破坏规律提供了捷径。数值模拟结果可作为生产现场必要的参考。但是,作者未考虑露天边坡的开挖以及台阶局部破坏同时顾及到计算规模和网格数目的庞大,没有模拟整个露天开采边坡,研究结果具有一定的局限性,在后期的研究中,可考虑更全面可行的模拟技术与方法;对建模耦合技术的理论研究不够深入,对其他方式的耦合方法未进行更深入的讨论和研究,其他耦合方法可能是后期研究的方向;对于四面体网格化还需进一步研究改进,使得网格更加符合数值模拟计算网格的需求。在这个学期的学期中,我深刻感受到徐老师渊博的专业知识,严谨踏实的治学态度,诲人不倦的育人精神,精益求精的工作作风和为科学事业忘我献身的崇高精神,也使我在学业上和思想上受益匪浅,这将是我终身学习的楷模。老师,您的上课内容非常生动,非常实用,很好的构造了一个知识体系,为我们的今后学习和工作有很大的帮助。希望老师在工作中更加顺利,生活更加美好!谢谢老师! 老师您辛苦了!