《岩土力学数值方法》结课论文
- 格式:doc
- 大小:39.50 KB
- 文档页数:4
岩土工程中数值流形方法的应用及研究【摘要】岩土工程中数值流形方法是一种新兴的计算方法,能够更准确地模拟岩土工程中的复杂问题。
本文首先概述了数值流形方法在岩土工程中的应用及研究意义,介绍了研究背景和目的。
然后详细解释了数值流形方法的基本原理,并给出了岩土工程中数值流形方法的具体应用实例。
探讨了数值流形方法在岩土工程中的优势和局限性,以及与传统方法的比较。
展望了岩土工程中数值流形方法的发展趋势,探讨了应用前景和未来研究方向。
数值流形方法有望在岩土工程领域取得更广泛的应用,为工程实践提供更准确、可靠的数值模拟结果,推动岩土工程技术的进步和发展。
【关键词】岩土工程、数值流形方法、应用、研究、基本原理、实例、优势、局限性、比较、发展趋势、应用前景、总结、展望、未来研究方向1. 引言1.1 岩土工程中数值流形方法的应用及研究概述数值流形方法是一种基于数学流形概念的数值计算技术,通过对数据集进行降维和重构,从高维数据中提取关键信息,实现对复杂系统进行建模和分析。
在岩土工程中,数值流形方法可以用于建立土体的模型,预测其稳定性和变形情况。
相比传统的数值模拟方法,数值流形方法具有更高的精度和效率。
本文将介绍数值流形方法在岩土工程中的具体应用实例,并分析其优势和局限性。
我们将比较数值流形方法与传统方法的差异,探讨岩土工程中数值流形方法的发展趋势。
通过对岩土工程中数值流形方法的应用及研究进行概述,我们可以更好地了解其在该领域的潜力和前景,为未来的研究提供指导和启示。
1.2 数值流形方法在岩土工程中的意义在岩土工程领域,数值流形方法的应用具有重要意义。
数值流形方法能够帮助工程师更准确地分析岩土体的结构和性质,从而预测土体的变形和破坏情况。
通过对土体的数值建模和仿真,可以更好地理解岩土工程中复杂的地质问题,为工程设计和施工提供更可靠的依据。
数值流形方法能够帮助工程师更快速地进行参数优化和工程优化,提高工程效率和质量。
通过数值流形方法的应用,可以有效地降低设计成本和风险,同时提高工程的安全性和可靠性。
计算岩土力学现有软件及其应用情况1.简介岩体是一种具有不连续性、非均质性、各向异性和非线性的天然地质体[1],由于岩体介质的复杂性, 分析边界和工况的多变, 用经典力学求解复杂的岩体工程问题的解析解已经很难, 甚至不可能。
随着计算机科学的进步,数值计算方法也不断发展,并且在岩土工程领域得到了广泛的应用。
例如,20世纪60年和70年代开始出现用于岩土工程稳定性计算的数值计算方法,主要是有限元法;20世纪80年代有限元、边界元及其混合模型得到广泛应用[2];20世纪90年代以来,离散元法、DDA法等方法层出不穷。
随着数值计算方法的百花齐放,计算岩土力学也从无到有发展起来,并且根据不同的计算方法形成了众多的计算软件,本文主要介绍较为常见且使用较广的几款计算软件。
数值计算方法大致分为连续变形分析方法和非连续变形分析方法两大类[3]。
连续变形分析方法主要有:有限单元法(FEM)、有限差分法(FDM)、边界单元法(BEM)等,其中以有限单元法应用最为广泛;非连续变形分析方法主要有:离散元法(DEM)、非连续变形分析(DDA)、流形元法(MM)、无单元法(EFM)等。
[4]2.有限单元法有限元分析中最基本的思想就是单元离散,即将求解区域剖分为若干单元,把一个连续的介质换成为一个离散的结构物,然后就各单元进行分析,最后集成求解整体位移。
通用的有限元软件主要有ABAQUS、ANSYS、COMSOL Multiphysics、ADINA、ALGOR、LS-DYNA、Nastran等,岩土工程专用有限元软件主要有PLAXIS、Midas GTS 、CRISP、Geoslope、GeoStudio、Rocsciences、Z-Soil和SoilVision等,本文仅简单介绍一些常见的有限元软件。
2.1ABAQUSABAQUS是一款由达索SIMULIA公司进行开发、维护及售后的有限元分析软件,在岩土工程数值分析方面,ABAQUS具有以下优点[5]:1)拥有能够真实反映岩土体实反映土体性状的本构模型,如土体的屈服特性、剪胀特性等。
岩土工程专业论文范文精选3篇我校对毕业设计质量十分重视,但近年来出现毕业设计质量下降的现象,为保证毕业设计的教学质量,在对有关院校毕业设计指导经验进行调研的基础上,结合我校岩土工程专业毕业设计的现状,分析影响毕业设计质量的主要原因,并针对如何提高毕业设计质量问题进行探讨。
一、影响岩土工程专业毕业设计(论文)质量的主要原因1.就业与考研的影响。
毕业设计(论文)通常安排在大四的最后一个学期进行。
大四学生在做毕业设计(论文)的同时,还面临着找工作、考公务员、研究生复试等诸多问题。
所以说学生本身主观上非常想投入大量的时间和全部的精力进行毕业设计(论文),力求呈现高质量的设计成果。
但现实情况迫使学生不得不压缩毕业设计的时间,从而影响毕业设计(论文)的完成质量。
2.指导教师自身实践能力的影响。
目前高校教师的主要来源是高校毕业生,这部分教师所占比例较大。
尽管他们具有较高的学历,但均是从校门到校门,没有企业工作的实践经验,缺乏实际动手能力,指导学生进行毕业设计(论文)的能力不足,从而造成学生的毕业设计成果质量下降。
3.毕业设计成果与实际需求脱节。
本科生毕业设计(论文)存在为设计而设计的倾向,毕业设计(论文)选题脱离实际,通常只是为了方便学生进行设计,从而对工程条件进行不恰当的简化。
这样即便学生在某一方面进行详细的设计,但设计方法单一,考虑实际问题不全面,不利于培养学生全面分析、解决实际问题的能力。
另外毕业设计成果缺乏展示平台,没有进行社会转化的机会,更产生不了行业价值、社会价值和经济价值,从而在一定程度上影响学生和指导教师的积极性。
二、提高岩土毕业设计(论文)质量的措施1.建立健全毕业设计(论文)监管机制。
学院作为毕业设计质量(论文)监管的二级单位,对学生毕业设计(论文)进行统一治理。
资环学院按照我校的本科生毕业设计(论文)治理办法,认真制定本学院本科生毕业设计(论文)治理细则,将监管工作落实到每一个环节,严把质量关。
岩土工程数值方法摘要:逐渐发展起来的一些岩土分析手段与数学理论,如信息量法、层次分析法、随机模拟法、无网络法、数值流形法、离散元法、分形理论、可靠度分析、人工神经元网络和智能岩石力学等,已经呈现出综合应用的趋势,对于岩体力学研究而言,岩石破坏过程的渐进性、岩体内部初始损伤的存在及块体之间的不连续特征是必须考虑的因素,因此建立在连续介质力学基础上的传统有限单元法具有明显的局限性。
各种新方法的涌现从不同方面推动了岩石力学数值计算方法的进步。
关键词:岩土数值模拟有限元法无网络伽辽金法扩展有限元法数值流形法离散元法Abstract: gradually developed some geotechnical analysis method and mathematical theory, such as information method, the analytic hierarchy process (ahp), random simulation method, the numerical manifold method, no network, discrete element method, fractal theory, reliability analysis, artificial neural network and intelligent rock mechanics etc, has presented a comprehensive application trend, for research in rock mechanics, rock failure process of rock mass progressive, the existence of the internal initial damage and block the discontinuous characteristics between is must consider factors so based on continuum mechanics on the basis of the traditional finite element method has obvious limitation. All kinds of the emerging of the new method from different aspects promote the rock mechanics numerical calculation method of progress.Keywords: geotechnical numerical simulation finite element method without network petro-galerkin method was expanded numerical manifold method finite element method of discrete element method中图分类号:O241 文献标识码:A文章编号:岩土数值模拟是否正确,其解决问题的重要基础仍然是地质工作,“地质体运动真实行为的理解比精确计算更为重要”。
岩土试验力学课程论文第一篇:岩土试验力学课程论文岩土试验力学课程论文题目:岩土试验力学发展现状和前景专业:岩土工程一、岩土力学试验1.岩土力学试验概况要很好的解决岩土工程问题、防灾、治灾,必须首先进行勘察与测试、试验与分析,并利用土力学、岩石力学、基础工程、工程地质学等的理论与方法,对各类工程进行系统研究。
因此,岩土力学试验是岩土工程规划设计、防灾的前期工程,也是地基与基础设计,治理地质灾害的不可缺少的重要环节。
2.岩土力学试验目的(1)了解岩石本身的物理和力学性质;(2)岩体质量分级、工程地质条件与问题评价;(3)边坡、地基和隧道围岩变形及稳定性分析,地质灾害防治工程方案论证等;(4)为岩土工程设计与施工提供参数和依据;(5)揭示岩土的变形规律和强度特征及破裂机理,建立其数学力学模型,进行岩土工程结构的力学分析。
3.岩土力学试验内容(1)岩石物理性质试验含水率、颗粒密度、块体密度;(2)岩石水理性质试验吸水性、渗透性、膨胀性、耐崩解性和冻融性。
(3)岩石力学性质试验单轴压缩强度和变形试验、三轴压缩强度和变形试验、抗拉强度试验、直剪强度试验和点荷载强度。
二、岩土试验力学概况岩土试验力学是土木工程岩土专业的一个分支,它是一门十分重要的技术基础课。
它主要包括学习岩土实验力学的基本理论,知道岩土的物理力学性质、强度变形计算、稳定性分析、挡土墙及基坑围护的设计与计算、地基承载力等岩土力学基本理论与方法。
结合有关交通土建、建筑工程、土木工程的理论和施工知识,分析和解决岩体工程及地基基础问题。
三、岩土试验力学的发展现状1.计算方面由于岩土材料比较特殊,那么在研究岩土试验力学方面就会比较复杂。
岩土体本身就是一个复杂的系统,具有不确定性,不规则性和不明确性。
目前,我国的岩土试验力学工作者倾向于采用理想数学模型和力学模型建立和描述岩土的各类特性,结果往往不是很理想,甚至出现很大的偏差。
那么,为解决这一现状,为突破创新,新的方法和技术是必不可少的。
岩土工程中的数值分析方法在研究生教学中的改革探讨发表时间:2020-10-20T11:19:51.620Z 来源:《教学与研究》2020年7月第18期作者:吴月秀姜清辉位伟卢海峰[导读] 随着经济的飞速发展,城市地铁、穿湖隧道等岩土工程朝着“高、深、大”的方向发展吴月秀姜清辉位伟卢海峰武汉大学土木建筑工程学院湖北武汉 430072摘要:随着经济的飞速发展,城市地铁、穿湖隧道等岩土工程朝着“高、深、大”的方向发展,如何应用数值分析方法来模拟岩土工程设计或施工、为实际工程提供指导是研究的热点和重点,也是培养行业内专业人才的一个重要方向。
然而由于岩土工程中的数值分析方法涉及的计算方法繁多、理论知识抽象难懂、编程要求高等原因,导致教学困难。
本文基于多年的实践教学,从岩土体特性的复杂性及数值分析方法的局限性两方面分析了该课程的特点,阐述了各方法在使用时的优缺点,拓宽了学生的专业视野;总结了当前教学中存在的问题如教学内容单一、教学方法过于抽象、教学与实际工程脱节、学生条件和诉求差异性大等,最后从调动学生学习兴趣、培养学生编程习惯、组织学生小组讨论、进行工程实例教学及改变考核方式几方面提出了相应的教学改革方案。
1、引言根据岩土工程的定义可知,其主要研究对象是岩土体,通过利用、整治或改造岩土体性质来解决工程问题,按照工作内容可分为岩土工程勘察、岩土工程设计、岩土工程施工、岩土工程检测及岩土工程管理,各个环节紧密联系,不可分割;涉及的工程有地基工程、边坡工程、隧道工程、铁路工程、矿山工程、环境工程、城市地铁工程、水利水电工程、地下工程、军事工程、近海工程等。
在《2019-2025年中国岩土工程行业市场全景调研与竞争格局分析报告》[1]一书中,详细阐述了岩土工程行业在国民经济中的重要性、行业发展规划、行业发展对社会发展的影响,分析了行业主要技术人才的现状。
在岩土工程实施过程中常使用的方法有:监测、数值模拟、室内试验、模型试验等。
岩土工程中常用的数值模拟方法综述-岩土工程论文-土木建筑论文——文章均为WORD文档,下载后可直接编辑使用亦可打印——岩土工程原理论文第八篇:岩土工程中常用的数值模拟方法综述摘要:从数值模拟实验的基本原理和典型实验案例介绍岩土工程问题分析中常用的数值模拟实验方法.通过对传统分析方法与数值分析方法的对比分析, 总结数值模拟实验的特点及适用范围.参数的选择和对比是目前数值模拟实验在实际应用中遇到的主要问题, 针对上述问题提出了多种数值模拟方法综合应用的解决思路.数值模拟实验在与地质力学模型实验及现场工程结合验证方面有广泛的应用前景, 是解决岩土工程问题的有效工具.关键词:数值模拟实验; 岩土工程; 应用实例; 参数对比;Application and Prospect of Numerical Simulation Experiment in Geotechnical EngineeringPENG Yanyan LIU Yuhang WANG Tianzuo DU Wei ZHENG ZhibinSchool of Civil Engineering, Shaoxing University Center of Rock Mechanics and Geohazards, Shaoxing UniversityAbstract:This paper introduces the numerical simulation experiment methods commonly used in the analysis of geotechnical engineering problems.Expanded from the basic principles of numerical simulation experiments and typical experimental cases, and through the comparative analysis of traditional analysis methods and numerical analysis methods, the characteristics and application scope of numerical simulation experiments are summarized.The selection and comparison ofparameters are the main problems encountered in the current practical application of numerical simulation experiments.This paper proposes a solution to the comprehensive application of varied numerical simulation methods for the above problem.Numerical simulation experiments, which have broad application prospects in combination with geomechanical model experiments and field engineering verification, are effective tools for solving geotechnical engineering problems..岩土是一般材料, 也是一种地质结构体, 它具有非连续、非均质、非线性的特性及复杂的加卸载条件和边界条件, 这使得岩土工程问题通常无法简单求解.因此数值模拟方法成了解决某些岩土工程问题的有效工具之一.从20世纪50年始, 人们就利用数值模拟方法对岩土工程问题进行大量研究, 经过60多年的发展, 针对岩土工程问题的数值模拟方法逐渐成熟, 成为岩土工程学科一个重要的研究方向, 使复杂岩土工程问题的设计发生了根本性的变化.岩土工程数值模拟不仅了传统线弹性力学实验, 而且也在岩石工程非线性实验中显示出极大的优势[1].岩土工程的数值模拟实验是对岩土工程活动和自然环境变化过程中岩体及工程结构的力学行为进行数值模拟的一种手段.在进行岩土工程数值模拟的时候, 需要对岩体进行分类:一是基于连续岩体力学的数值模拟, 二是基于不连续岩体力学的数值模拟.有限元法、边界元法、有限差分法把岩体看作连续介质进行模拟, 而离散元法、不连续变形分析法、数值流形法把岩体看作不连续介质进行模拟, 本文会对上述方法和应用进行总结介绍.1 数值模拟实验在岩土工程中的优势岩体本身具有非均质、非连续、非线性及复杂的加卸载条件和边界条件的特点, 加之岩体所处环境也比较复杂, 岩体工程开挖前就受地应力、地下水、周围温度等耦合作用, 所以很难建立完善的地质力学模型[2].1.1 传统研究方法的局限性由于岩体的不连续性与非均质性, 实验室岩石试样具有明显的尺寸效应, 并且通过实验得到的结论往往与实际工程相差甚大[3].通过现场原位实验的方式取得数据, 也非常艰难.岩体力学原位实验一般耗资巨大, 且受地形、地质、施工条件限制, 得到的实验结果不具代表性, 难以推广到其他工程[2].1.2 数值模拟在岩土工程中的优越性数值模拟在岩土工程中的适用范围非常广, 并且节约资金,它不仅能模拟岩体复杂结构特性, 还能研究岩土工程活动对周围环境的影响, 并对工程灾害进行预报.通过对现场原位实验的实测与反分析, 可以获得节理岩体的等效力学模型, 逐步成本较高的原位实验, 加快工程进度, 且可应用于各种地形、地质与施工条件, 推广实验结果的应用范围与使用条件[3].因此用数值模拟的方法来解决岩土工程问题是行之有效的.2 岩土工程中常用的数值模拟方法2.1 有限元法2.1.1 有限元法原理有限元法是利用变分的原理去求解数学物理问题的一种数值模拟方法.有限元法最早由布理克(W.Blake) 在1966年引入岩土工程领域, 用来解决岩土工程问题.有限元法基于最小总势能原理通过解方程组的方法来求解, 是目前岩土工程领域中应用最广泛的数值模拟方法.有限元法是用多个彼此相联系的单元体所组成的近似等价物理模型来代替实际的结构或者连续物体, 通过结构及连续体力学的基本原理及单元的物理特性建立起表征力和位移的关系去建立方程组, 解方程求其基本未知物理量, 并由此求得各个单元的应力、应变及其他辅助值[4].2.1.2 有限元法的应用有限元法由线性发展到非线性和大变形问题的应用(二维发展到三维) , 目前还可考虑流变、温度与应力场耦合, 损伤、渗流、断裂以及波动和动力效应[5].刘庭金等[6]利用有限元分析矿山、地铁等地下工程由于洞室开挖引起的围岩卸载过程中, 洞室孔壁围岩附近发生的损伤演化和应力场调整全过程进行分析.郑颖人等[7]对有限元强度折减法的计算精度和影响因素进行了详细分析, 包括屈服准则、流动法则、有限元模型本身以及计算参数对安全系数计算精度的影响, 并给出了提高计算精度的具体措施.应用于岩质边坡的稳定分析, 得到了岩质边坡的滑动面和安全系数, 开创了求节理岩质边坡滑动面与稳定安全系数的先例.2.2 边界元法2.2.1 边界元法原理边界元法同有限差分法和有限元法一样, 都是一种用来解决边值问题的数值分析方法, 它可以用来解决弹性力学、塑性力学, 以及热传导、地下水力学等方面的问题, 发展历史悠久, 但直到20世纪60年代后期在计算机技术得到发展时, 边界元法才成为实际可行的一种数值模拟方法[8].由于岩土工程问题的复杂性, 边界元法在1976年才引入到岩土工程中.边界元法是通过求解边界积分的方法来求解边值问题, 在边界元上划分单元, 求边界积分方程的解, 进而求出区域内任意点的场变量, 所以边界元法也称边界积分方程法.边界元法又分以互等功原理为基础的直接法, 和以叠加原理建立起来的间接法.2.2.2 边界元法的应用边界元法和有限元法比起来, 可以用降维的方法来简化计算(三维问题二维化, 二维问题一维化) , 不但计算起来方便, 而且计算精度高, 但是面对非连续、非线性介质问题边界元法则比较难适应.目前边界元法主要在地下工程开挖、土体结构相互作用及地下水流动过程的一般应力和变形分析有着应用.虽然边界元法的适用范围有限, 但是和其他数值方法的联合使用能充分发挥其优越性, 为解决岩土工程问题开辟了新的途径.例如, 在线弹性区域或无限域、半无限域可以采用边界元法, 在非线性的区域采用有限元法, 发挥两种算法各自的优势, 使计算效率及精度得到提高, 对工程实际应用有很大的帮助[9].马天寿等[10]用边界元法对页岩地层井眼坍塌问题进行了分析并得出弹性模量各向异性、水平地应力差异和钻井液密度等对井壁应力分布影响较大, 而泊松比各向异性的影响较小的结论, 如图1.图1 边界元方法解出的井周应力分布图[10]2.3 有限差分法2.3.1 有限差分法原理有限差分方法是以最小势能原理, 通过解方程的方式进行求解.这种方法是一种最古老的求解方程组的数值方法, 在计算机出现以前一般的手摇计算器也可求解.20世纪80年代末由美国ITASCA公司开发的FLAC程序广泛采用差分方法进行求解, 并且在岩土工程数值计算中得到了广泛应用[5].2.3.2 FLAC3D的应用鉴于有限差分法单独在岩土工程中的应用并不多, 一般都是基于有限差分法的FLAC程序进行岩土工程计算, 所以这里讲的是FLAC程序的应用.FLAC程序可用来模拟地质材料的大变形、失稳、动力、流变、支护、建造及开挖等问题, 同时还可以模拟渗流场和温度场对岩土工程的影响.李为腾等[11]解决了FLAC3D中CABLE单元无法实现锚杆(索) 破断失效的问题, 并采用Fish语言编程, 将修正模型嵌入到FLAC3D 主程序中, 实现锚杆破断失效的单元化, 如图2.图2 FLAC模拟修正影响系数随地应力变化曲线[11]2.4 离散元法2.4.1 离散元法原理离散元法是Cundall在1971年所提出来的, 后经Voegele等人的发展, 成为一种新的数值模拟方法[12].离散元法是以牛顿运动定律的显示求解的数值方法, 离散元法也要将区域划分为单元, 但是单元因受节理、劈理等不连续面的控制, 在以后的运动过程中, 单元节点可以分离, 即一个单元与其邻近单元可以接触, 也可以分开.单元之间相互作用的力可以根据力和位移的关系求出, 而个别单元的运动则完全可以根据单元所受的不平衡力和不平衡力矩的大小按牛顿运动定律确定[9].2.4.2 离散元法的应用离散单元法可分为动态松弛法和静态松弛法两种.目前常用的大多是动态松弛法.动态松弛法用解决动力学问题的方法去求解非线性静力学问题, 用显式中心差分法近似对运动方程进行积分计算, 并假设块体在运动时将动能转化成热能耗散掉, 把人工黏性阻尼引入计算, 使系统达到平衡状态、运动趋于稳定[13].常晓林等[14]用离散元的方法模拟了岩石介质从小变形到大变形再到破坏的全过程.李晓柱等[15]用离散元的方法分析了堆石坝现场碾压实验, 验证离散元数值模拟方法应用于堆石坝碾压特性研究的可行性, 更直观地从细观角度解释堆石体碾压过程中宏观参数(如干密度等) 的变化规律, 为大坝I区堆石料选取科学合理的碾压施工参数, 为堆石体碾压特性研究提供新的途径.王贵君[16]针对国内某高速公路隧道工程, 应用离散单元法对节理裂隙岩体中不同埋深无支护暗挖隧洞的稳定性及其机理进行了数值模拟.2.5 不连续变形分析法(DDA法)2.5.1 不连续变形分析法原理不连续变形分析方法一般简称DDA法(Discontinuous Deformation Analysis) , 由石根华首创, 基于岩体介质非连续性发展起来的, 以模拟复杂加载条件下离散块体系统不连续大变形的力学行为为目的的一种数值方法[1], 该方法以最小势能原理为基础, 通过解方程的方式求解问题, 适用于发生大变形, 岩体发生非连续破坏的情况.石根华在1988至1989年期间开发了二维DDA程序, 用来进行前处理、正分析、反分析和后处理的计算, 程序包括DDACUT, DDA FOR-WARD, DDA BACKWARD和DDAGRAPH, 4个部分[1].2.5.2 不连续变形分析法的应用DDA法在岩土工程当中主要应用于两个方面:一是在爆破工程研究中的应用, 二是在地下岩体开挖工程中的应用.张丽娟等[17]采用DDA方法对台阶爆破抛掷的全过程进行了模拟分析, 给出了人工边界、爆破荷载的处理算法.邬爱清等[18]利用DDA的方法分析了块体稳定性验证及其在岩质边坡性分析中的应用, 如图3.吴建宏等[19]将DDA法应用在岩石边坡失稳数值仿真当中.图3 DDA计算后块体系统变形与破坏[18]2.6 数值流形法2.6.1 数值流形法原理在1995年, 石根华提出数值流形法(NMM) , 数值流形法是利用现代数学中流形的有限覆盖技术建立起来的一种新的数值计算方法, 将有限元、不连续变形分析(DDA) 和解析方法统一到一种计算方法中, 它吸收了有限元、DDA和解析法各自的优点, 通过分片光滑的覆盖函数, 对连续和非连续问题统一了计算格式, 是一种十分适合于岩土工程分析的数值方法[13].2.6.2 数值流形法的应用数值流形法可同时处理非连续问题与连续问题.但由于网格的连接与单元划分的限制, 数值流形法在开裂计算方面仍存在一定的困难.目前研究主要集中在连续与非连续问题的求解和裂纹扩展的模拟[20].李树忱等[21]充分利用数值流形方法中两套网格的特点, 采用围线积分法来计算应力强度因子和最大周向应力确定裂纹的扩展角, 模拟裂纹的扩展过程.钱莹[22]等利用三角形元素来构成数学网格, 利用流变法来分析爆破振动对边坡动态稳定性的影响, 如图4.图4 流形元模型网格[22]3 数值模拟实验存在的问题及发展尽管数值模拟方法取得了很多研究成果, 但从当前数值模拟在岩土工程中的应用来看, 主要存以下几个问题:(1) 参数选取的问题.无论哪种数值模拟方法都必须准确选定岩石或岩体的物理力学参数.但由于岩体本身与所处环境的复杂性, 确定这些参数并非易事[1], 因此得到的数据和工程实际也有出入.(2) 计算机储存量不足和计算速度受限的问题[23].由于越来越多大型工程的兴建, 当前岩土工程的数值模拟迫切需要发展并行计算方法, 用多核或多联计算机将原有的程序做重构, 以提高计算速度.(3) 岩体进入破坏模式后力学机制发生转变, 数值模拟中的相关破坏准则难以准确描述破坏后的岩体力学行为[24].岩土工程中的数值模拟方法有很多, 每种模拟方法都有其优点, 但也各局限性.要想用单一的方法解决岩土工程问题是不现实的, 为了更好地发挥每一种方法的长处, 多种数值模拟方法综合应用成了近年来岩土工程数值模拟实验的新趋势.如在围岩影响区采用非连续分析法, 在围岩影响区之外的原岩应力区用连续介质分析法.但是要解决两种不同算法在围岩影响区和原岩应力区交界面的位移和应力的协调问题.遗憾的是, 目前专注于上述问题的学者有限, 尚未有突破性进展.4 总结本文介绍了数值模拟实验在岩土工程中的研究进展, 列举了先进的数值模拟实验方法, 并分别介绍了相关原理及在岩土工程中的应用, 最后指出了目前数值模拟实验存在的问题及其发展方向.总而言之, 随着岩土工程规模日益扩大, 对岩土工程建设的科研设计水平和建设精度的要求越来越高, 用数值模拟去解决岩土工程问题是一种有效手段.数值模拟实验对掌握岩土工程围岩应力变化规律和变形破坏规律有重要意义, 能够对实际工程的支护设计提供理论依据.参考文献[1]唐广慧, 刘发祥, 唐升贵, 等.当前岩石力学数值计算方法应用探讨[J].西部探矿工程, 2007 (12) :25-29.[2]李宁, 辛有良.岩石力学数值方法的作用与地位浅析[J].陕西水力发电, 1997 (2) :19-22.[3]李宁, G SWOBODA.当前岩石力学数值方法的几点思考[J].岩石力学与工程学报, 1997 (5) :104-107.[4]朱立仁, 黄玉凯, 焦向东.有限元法在宝日希勒露天煤矿边坡稳定性分析中的应用[J].露天采矿技术, 2011 (5) :7-8+11.[5]佘诗刚, 董陇军.从文献统计分析看中国岩石力学进展[J].岩石力学与工程学报, 2013, 32 (3) :442-4 .[6]刘庭金, 朱合华, 唐春安.围岩卸载损伤演化及应力场调整有限元分析[J].地下空间, 2002 (4) :310-313+319-371.[7]郑颖人, 赵尚毅.有限元强度折减法在土坡与岩坡中的应用[J].岩石力学与工程学报, 2004 (19) :3381-3388.[8]王泳嘉.边界元法在岩石力学中的应用[J].岩石力学与工程学报, 1986 (2) :205-222.[9]鄢建华, 汤雷.水工地下工程围岩稳定性分析方法现状与发展[J].岩土力学, 2003 (S2) :681-686.[10]马天寿, 陈平.应用边界元法分析页岩地层井眼坍塌问题[J].中南大学学报(自然科学版) , 2016, 47 (3) :839-849.[11]李为腾, 杨宁, 李廷春, 等.FLAC3D中锚杆破断失效的实现及应用[J].岩石力学与工程学报, 2016, 35 (4) :753-767.[12]王泳嘉.离散元法及其在岩石力学中的应用[J].金属矿山, 1986 (8) :13-17+5.[13]孔德森, 栾茂田.岩土力学数值分析方法研究[J].岩土工程技术, 2005 (5) :249-253.[14]常晓林, 胡超, 马刚, 等.模拟岩体失效全过程的连续-非连续变形体离散元方法及应用[J].岩石力学与工程学报, 2011, 30 (10) :2004-2011.[15]李晓柱, 刘洋, 吴顺川.堆石坝现场碾压试验与离散元数值分析[J].岩石力学与工程学报, 2013, 32 (S2) :3123-3133[16]王贵君.节理裂隙岩体中不同埋深无支护暗挖隧洞稳定性的离散元法数值分析[J].岩石力学与工程学报, 2004 (7) :1154-1157.[17]张丽娟, 查珑珑, 郭剑鸿.非连续变形分析(DDA) 法及在岩石力学中的应用[J].山西建筑, 2009, 35 (8) :3-5.[18]邬爱清, 丁秀丽, 卢波, 等.DDA方法块体稳定性验证及其在岩质边坡稳定性分析中的应用[J].岩石力学与工程学报, 2008 (4) :6 -672.[19]吴建宏, 大西有三, 石根华, 等.三维非连续变形分析(3D DDA) 理论及其在岩石边坡失稳数值仿真中的应用[J].岩石力学与工程学报, 2003 (6) :937-942.[20]王书法, 朱维申, 李术才, 等.加锚岩体变形分析的数值流形方法[J].岩石力学与工程学报, 2002 (8) :1120-1123.[21]李树忱, 程玉民.数值流形方法及其在岩石力学中的应用[J].力学进展, 2004 (4) :446-454.[22]钱莹, 杨军.流形元法分析爆破振动对边坡动态稳定性影响[J].有色金属(矿山部分) , 2007 (5) :35-38.[23]朱维申, 赵成龙, 周浩, 等.当前岩石力学研究中若干关键问题的思考与认识[J].岩石力学与工程学报, 2015, 34 (4) : 9-658.[24]彭岩岩.层状岩体深部巷道变形破坏物理模拟及红外探测研究[D]. :中国矿业大学, 2015.。
岩土数值极限分析方法的发展与应用一、本文概述随着科学技术的不断进步和工程实践的日益深化,岩土工程的数值极限分析方法在工程安全评估、优化设计以及风险控制等方面发挥着越来越重要的作用。
本文旨在全面概述岩土数值极限分析方法的发展历程、现状以及未来趋势,并深入探讨其在各类岩土工程中的应用。
本文将首先回顾岩土数值极限分析方法的起源与发展,梳理其从早期的简单理论模型到现代复杂数值分析技术的演变过程。
接着,文章将重点介绍当前主流的数值极限分析方法,包括有限元法、有限差分法、离散元法等,并分析它们各自的优缺点和适用范围。
本文还将探讨岩土数值极限分析方法在岩土工程中的应用案例,如边坡稳定性分析、隧道开挖模拟、地下工程安全评估等,以展示其在实际工程中的重要作用。
本文将展望岩土数值极限分析方法的未来发展趋势,包括技术创新、方法优化、多学科交叉融合等方面,以期为相关领域的研究人员和实践工作者提供有益的参考和启示。
通过本文的阐述,希望能够推动岩土数值极限分析方法在岩土工程领域的进一步发展与应用。
二、岩土数值极限分析方法的发展历程岩土数值极限分析方法的发展历程可以追溯到20世纪中期,随着计算机技术的飞速发展和数值计算方法的不断创新,岩土数值极限分析逐渐成为一种重要的研究手段。
其发展过程大致可以分为以下几个阶段:初期探索阶段:在20世纪50至60年代,研究者开始尝试运用数值方法对岩土体的极限状态进行分析。
当时主要采用有限元法等基本的数值计算方法,对岩土体的应力、应变和位移等进行了初步的探索。
这一阶段的研究虽然较为基础,但为后续的发展奠定了坚实的基础。
方法发展阶段:随着计算机技术的不断进步和数值计算方法的日益成熟,岩土数值极限分析方法在20世纪70至80年代得到了快速发展。
研究者开始尝试运用更加复杂和精确的数值方法,如离散元法、边界元法、有限差分法等,对岩土体的力学特性、破坏模式和极限承载能力等进行了深入的研究。
这些方法的出现极大地丰富了岩土数值极限分析的手段,提高了分析的准确性和可靠性。
岩土工程数值方法发表时间:2018-07-25T16:28:04.553Z 来源:《建筑学研究前沿》2018年第7期作者:林杰[导读] 逐渐发展起来的一些岩土分析手段与数学理论,如信息量法、层次分析法、随机模拟法、无网络法林杰广东金山河建筑基础工程有限公司 523000 摘要:逐渐发展起来的一些岩土分析手段与数学理论,如信息量法、层次分析法、随机模拟法、无网络法、数值流形法、离散元法、分形理论、可靠度分析、人工神经元网络和智能岩石力学等,已经呈现出综合应用的趋势,对于岩体力学研究而言,岩石破坏过程的渐进性、岩体内部初始损伤的存在及块体之间的不连续特征是必须考虑的因素,因此建立在连续介质力学基础上的传统有限单元法具有明显的局限性。
各种新方法的涌现从不同方面推动了岩石力学数值计算方法的进步。
文网 http://www.x 关键词:岩土工程;数值分析;方法在解决不同的岩土工程问题时,基本方程中的运动微分方程、有效力原理、连续方程和几何方程的表达式基本相同;但方程却与前几种有着较大的差异" 当涉及到具体的岩土工程问题时,我们可以根据具体的边界条件以及相关的初始条件来解决上述问题,所采用的方法一般为数值分析法" 通过以上论述我们可以发现,岩土工程一般都需借助数值分析法进行解决。
下面我们了解这门科学。
一、数值分析在岩上工程中的地位众所周知,岩土是自然的产物,具有较强的区域性,初始应力也较难预测,在分析岩土工程时,首先需要掌握的就是工程的地质条件,此外还要掌握岩上的工程性质,再者需要掌握的就是力学中一些最为基木的概念,并且能够在此基础上利用公式以及数值分析方法来解决问题、在计算时,需要能够做到因地制宜,要能够对具体问题做出具体的分析,然后将得出的结论运用到工程建设之中;在实际的岩土分析过程中,数值分析所得出的结果对于工程师的综合判断而言是极其重要的。
在实际的操作之中,在对岩土工程的对象进行分析时要做好岩土材料特性的分析工作,还要注意结合岩土工程的初始条件以及边际条件进行综合的分析;就目前的岩土分析发展状况而言,还不能得出具体的解析解,只能适用于定向分析,岩土工程在进行设计时要将重点要放在概念设计上,并且参照岩土工程师的判断结果,岩土数值分析所得出的结果也是岩土分析的关键参数。
岩土力学数值方法的应用及发展前景班级:硕研14级4班姓名:吴纪达学号:140431055 专业方向:结构工程1 引言随着计算机技术进入各个学科领域,近年来岩土力学数值分析方法有了迅速的发展,越来越多的生产、设计和科研部门采用数值方法对岩土工程进行分析,出现了一系列的数值方法,这些数值方法可大致分为连续变形分析方法和非连续变形分析方法两大类。
本文简要讨论了目前已较为成熟应用于岩土工程中的多种数值方法,介绍了最新尝试应用于岩土工程中几种新的数值方法,并在此基础上对今后岩土工程数值方法的发展作了一些展望。
分析了采用连续介质力学分析岩土工程问题的关键,并讨论分析了岩土本构理论发展现状,提出对岩土本构理论发展方向的思考,最后对数值分析在岩土工程分析中的地位作了分析。
分析表明,岩土工程数值分析结果是岩土工程师在岩土工程分析过程中进行综合判断的重要依据之一;采用连续介质力学模型求解岩土工程问题的关键是如何建立岩土的工程实用本构方程;建立多个工程实用本构方程结合积累大量工程经验才能促使数值方法在岩土工程中由用于定性分析转变到定量分析。
2 岩土工程中当前较为常用的几种数值方法2. 1 有限单元法有限元法自20世纪70年代提出发展至今,已经相当成熟,是目前最广泛使用的数值方法之一。
有限元分析可以处理复杂的边界条件以及材料的非均匀性和各向异性,还可以有效地模拟材料的非线性应力应变关系,得到围岩的整体应力场、位移场,而得出围岩的可能破坏部位,可对围岩的力学特性有全面了解。
但如何从有限元分析的结果得出合适的衡量围岩稳定程度的定量指标还需进一步研究。
有限单元法也遇到了本身不可克服的困难,如在前处理困难,计算应力应变解答不连续,难以进行任意路径开裂计算等。
2. 2 界面单元有限元法为了推广有限单元法,使其能够处理简单的非连续问题,人们提出了各种能够反映非连续变形性质的简单力学模型和特殊界面单元,并将其应用于单一性非连续界面力学行为的模拟中。
界面单元有限元法将岩土介质视为准连续介质,仍以连续分析为主,但可以对个别具有控制作用的宏观非连续面的变形与破坏力学行为进行重点分析,在工程中得到了广泛应用。
同时,该方法也存在诸多的局限性,如只能对原生的非连续界面进行计算,对次生的非连续界面无法进行处理、界面单元数目不能设置太多、界面弹簧刚度的选取较为困难等。
因此,该方法在处理复杂的非连续变形问题时仍显得无能为力。
2. 3 离散元法(DEM)、块体理论及刚体有限元法这是分析多节理岩体的另一类重要方法。
在这类方法中通常把岩体当成是由节理切割的岩块组成的结构体。
离散元法、块体理论、刚体有限元法是三种有代表性的不连续块体方法。
离散单元法是Cundall于1971年提出来的,该方法与在时域中进行的其他显示计算相似例如与解抛物线型偏微分方程的显示差分格式相似。
离散单元法也像有限单元法那样将区域划分成单元。
但是单元因受节理等不连续面控制在以后的运动过程中单元节点可以分离即一个单元与其邻近单元可以接触也可以分开适合节理岩体的应力分析。
单元之间相互作用的力可以根据力和位移的关系求出而个别单元的运动则完全根据单元所受的不平衡力和不平衡力矩的大小按牛顿运动定律确定。
其基本思想是岩块之间的相互作用同时受表征位移—力的物理方程和反应力—加速度(位移、速度)的运动方程支配通过迭代求解显示岩体的动态破坏过程。
块体理论和刚体有限元法则属于静态分析的方法。
仅考虑在恒定荷载作用下系统的平衡及失稳条件。
块体理论实质上是一种几何学的方法,根据岩体中实际存在的节理倾角及其方位,利用块体之间的相互作用条件找出具有移动可能性的岩体及其位置,故也称关键块理论。
该理论是一种通过判别和描述洞室围岩最危险岩石块体运动来确定岩体稳定性的分析方法,是由石根华提出并在美与Goodman 合作完善起来的。
认为岩体完全被交错节理所切割,形成各种形状的空间块体,而在这种岩体中开挖的洞室,从方向和位置不利且以开挖临空面为一个或若干面的块体的运动开始的,这个块体称为关键块体。
对这种理论的一个新发展是边边接触块体理论,它除了与关键块体理论相同的假定之外还假定块体之间接触可以沿节理全长或部分长度分布。
刚体有限元法则把岩块视为刚体块单元,块体之间以具有相应刚度的弹簧(阻尼器)连结形成总的系统,建立相应的基本方程。
用刚体有限元法( RBFEM) 给出的应力精度不会低于甚至高于位移精度,优于传统的以位移法为基础的有限元数值解。
此外RBFEM 大大地减小了总体的半带宽及体积,大幅度降低了计算量。
这些优点有利于结构分析,特别是非线性问题的数值解。
这类模型如有刚体-弹簧模型,刚性有限元,分块刚体位移-界面应力元法,块体-弹簧模型。
在这类模型中,界面特性均假设服从Coulomb摩擦定律,考虑体系的静力约束条件比较充分,对于连续状态的应力分析往往可给出较高的计算精度,对于临界状态能够用于估算出极限荷载,并可有效地用于少量块体界面间的摩擦接触分析。
但这类方法过分强调岩土体结构面作用,而对于岩体的结构体的变形未给予足够重视。
虽然可以用于分析原生界面的破坏分析,但不能模拟实际岩体的破坏发展过程,对于破坏发展导致的次生界面的非连续变形行为以及块体失稳的运动过程也无法进行模拟。
以上三种不连续块体方法在模拟节理岩体方面,克服了连续体力学的局限性,在模拟岩体运动趋势及失稳条件方面有其优点,然而也存在一些缺点,如:离散元法中关于块体之间阻尼系数、运算的时间步长等参数的确定带有极大的随意性和盲目性,至今没有作为确定这些参数可遵循的原则;三种方法中关于块体单元的划分带有随意性,对最终成果的影响尚有待探讨。
2. 4 耦合方法除了上述各种数值方法之外,耦合的方法也得到了发展与应用。
如有限元与边界元的耦合、有限元与离散元的耦合、有限元与界面元的耦合、离散元与边界元耦合等,这些耦合方法可分别发挥各种方法的优点并进行耦合提高,将非连续体的模拟与无限介质辐射阻尼的模拟统一到一个模型中,为地下结构和岩质边坡的抗震稳定性分析提供了全新的手段。
3 岩土工程中最新应用的几种数值方法3 . 1 界面元法在岩体工程中,无论是变形或是强度或是稳定性,往往受控于软弱结构面的变形和受力特性,因此软弱结构面部位是岩体工程问题分析的热点和难点,也是重点,倘若能建立一类以界面变形的受力分析为重点的数值离散模型无疑将赢得工程界的重用和青睐。
我国学者卓家寿等人前人的基础上提出了界面元离散模型和界面元法。
首先在模型中以界面元取代弹簧元,这充分体现建模中将单元变形累积于界面的思想,且可用不同类型元作(如弹性元件、塑性元件、粘性元件、开裂元件和接触元件等)反映单元的弹性、塑性和粘性变形,描述结构的开裂过程或接触界面上错位、滑移和张开度;其次,在略去高阶小量的假定基础上根据平衡律、协调律和本构律,导出非均匀、各向异性体中相邻单元交界上任意一点的界面应力公式;进而,引入自然坐标导出了任意块体形状界面元的有关列式,给出求解块体内部应力的途径和公式。
在以上工作基础上进而完成各类非线性、动力、稳定、梁板、问题和耦合问题的界面元法研究及其工程应用,从而使界面元法不仅能求解有限元所能计算的问题,而且可以弥补有限元在求解不连续介质问题和开裂等问题的不足,显示了界面元法的强劲生命力,形成了独具一格的新的数值分析体系。
界面元法具有的优点:适用性强;精度高;计算便捷;相容性强。
3. 2 非连续变形分析(DDA)方法非连续变形分析方法(DDA)是继离散单元法之后,于20世纪80年代末发展起来的一种更新的用于模拟散体系统力学响应的数值分析方法。
与有限单元法相似,非连续变形分析方法用隐式方程表达,所不同的是它引入了运动方程,用最小势能原理把块体之间的接触问题和块体本身的变形问题统一到矩阵的求解中,理论严密,精度较高,而且把静力和动力、正分析和反分析统一起来,不仅可以计算破坏前的小位移,也可以计算破坏后的大变形,对滑动、崩塌、爆炸和贯入等问题也十分有效,是一种不同于DEM的新的数值计算方法。
由于DDA方法在运动约束方面做得比较充分,理论上比较严密,且静力和动力分析采用了统一的数值计算格式。
因此,在结构、岩体和土体的非连续大变形力学过程模拟方面发挥了较大潜力。
3. 3 流形方法在以上讨论的方法中,有限元等方法将岩体理想化为完全连续介质,离散单元法等又将岩体抽象为完全非连续介质,这其实是走向了问题的两个极端。
将连续变形与非连续变形统一起来的方法更适合实际岩体变形分析。
于是有了一种更新颖的数值方法,这种新颖的数值方法是流形方法,它引起了广大学者的兴趣。
流形方法由石根华于1992 年提出,它是应用流形覆盖技术建立的一种把有限元法、非连续变形分析和解析方法包含在内的全新的统一计算方法。
这种新的方法由数学覆盖和物理网格组成。
数学覆盖由占整个材料体的许多有限重叠的数学覆盖组成,而物理网格包括材料体的边界、裂缝、块体和不同材料区域的交接面。
该方法的优点是具有相对完善的非连续变形处理功能,可以在统一的数学理论框架下同时处理连续与非连续问题,较之有限元法更适合于开裂模拟。
3. 4 无单元法对于岩体这种介于连续介质与完全非连续介质之间的材料,沟通连续变形与完全离散体的非连续变形的桥梁是岩体中的地质非连续面的扩展、贯通等的破坏行为,它是联接岩体变形过程中初始阶段与终了阶段的一条纽带,在整个岩体变形过程中,这些非连续面的变形、扩展以及失稳对整个工程岩体结构的稳定是至关重要的。
而现有的基于网格划分的方法均是基于单元离散思想的,在分析裂纹扩展时,都面临两方面的问题:网格重新划分的问题;裂纹扩展结果受到网格划分形式严重影响的问题。
无单元类方法因具有无须单元网格划分、前后处理简单、较传统有限单元法更适合断裂问题的计算分析等优点而受到学术界的广泛关注。
无单元类方法在进行裂纹扩展模拟时不再存在传统有限单元法的重新剖分网格的困难,而仅仅在裂尖局部区域内布置节点,大大简化了前处理过程,在剖分策略上,无单元法较传统的有限单元法更适于断裂问题的计算分析。
对于以无单元伽辽金法(EFGM)为代表的无单元类方法,现有研究成果大多集中在裂纹扩展的模拟方面,在将无单元方法应用于岩土力学的数值计算方面仍存在很多困难,如非连续材料插值函数的构造、摩擦接触问题和多体相互作用问题的处理等还需作进一步的研究。
4 岩土工程数值分析方法发展展望随着岩土工程领域的不断扩展与延伸,数值分析方法得到了长足的发展,分析模式不断改进,分析精度不断提高,,以有限元为代表的连续性数值分析方法在实际岩土工程中应用较为广泛,从学科发展来看,非连续变形分析将具有广阔的发展空间。
但岩土工程材料是一种复杂的地质材料,具有高度非连续性、非均匀性以及各向异性等地质特点,在力学性质上表现出强烈的非线性、非弹性等力学行为。