岩土工程数值分析方法
岩土勘察数据处理与分析方法
岩土勘察数据处理与分析方法岩土勘察是土木工程中一项至关重要的工作,它的目的是获取有关地质、土壤和岩石等岩土体的详细信息。
这些信息将直接影响到土木工程设计和施工的效果。
然而,采集到的原始勘察数据并不直接适用于工程设计,因此需要进行数据处理与分析,以便得出准确、可靠的结果。
本文将介绍一些常用的岩土勘察数据处理与分析方法。
第一步,在进行数据处理与分析之前,需要对原始数据进行质量控制。
这包括对数据的可靠性进行评估,如检查测试设备的准确性、观测者的技术能力等。
同时,还需要进行数据的整理,删除异常值和错误数据。
只有经过质量控制的数据才能保证分析结果的准确性和可靠性。
第二步,数据处理与分析常用的方法之一是数据统计。
通过对采集到的数据进行统计分析,可以得到数据的分布规律、集中趋势等信息。
例如,可以计算岩土体的平均值、中位数、标准差等指标,揭示岩土体的基本特征。
此外,还可以通过绘制频率分布图、箱线图等图表,直观地展示数据的分布情况。
第三步,数据处理与分析常用的方法之二是插值方法。
在进行土层分析和地下水分析时,常常需要根据已知样点的数据推算出未知位置的数据。
这时,插值方法可以发挥重要作用。
常用的插值方法有反距离权重法、克里金法等。
这些方法不仅可以揭示未知位置的数据分布情况,还可以为工程设计提供有效的参考。
第四步,数据处理与分析常用的方法之三是地质统计方法。
在进行土质分析时,需要对岩土体的成分和特性进行分析。
地质统计方法包括聚类分析、因子分析、主成分分析等。
这些方法可以识别出不同土质类型的主要控制因素,为工程设计提供有力依据。
第五步,数据处理与分析常用的方法之四是地质力学方法。
岩土体的力学特性对土木工程的承载能力和稳定性至关重要。
地质力学方法包括岩土体参数的反演、试验与模拟分析等。
通过这些方法,可以确定岩土体的强度参数、变形特性等,为工程设计提供可靠的依据。
综上所述,岩土勘察数据处理与分析是土木工程设计不可或缺的一环。
数值分析在岩土工程中的应用课件
数值分析可以模拟土壤的力学行为和地基 工程的地震、渗流等复杂因素,提高地基 工程的可靠性和安全性。
地下水工程与水环境
数值模拟与预测
数值分析可以模拟地下水的流动和污染物 扩散等过程,为水资源开发和环境保护提 供支持。
数值分析可以通过对历史数据的模拟和分 析,预测岩土工程的未来趋势和可能遇到 的问题,为决策提供科学依据。
数值分析在岩土工程中的应用成果已经得到了广泛认可和应用,为工程实践提供了 更为有效的技术支持和解决方案。
数值分析的未来发展趋势
随着计算机技术和数值计算方法的不断发展,数值分析在岩土工程中的 应用将更加广泛和深入。
未来,数值分析将进一步考虑更多的物理和化学效应,建立更为精细和 全面的模型,提高模拟的准确性和精度。
数值分析的重要性
01
随着岩土工程规模的扩大和复杂 性的增加,传统实验方法已经难 以满足实际工程的需求。
02
数值分析可以弥补实验方法的不 足,通过计算机模拟实验过程, 提高分析的精度和效率。
岩土工程中的数值分析应用
岩石力学与岩石稳定性分析
土力学与地基工程
数值分析可以模拟岩石的力学行为和稳定 性,为矿山、隧道等岩石工程的设计提供 依据。
土壤水盐分运移分析 地下水资源保护及利用方案设计
其他工程实例
桩基承载力及沉降分析
边坡支护结构选型及设 计
01
02
03
地基处理方案优化设计
04
工程材料力学性能及参 数确定
05
结论与展望
数值分析在岩土工程中的应用总结
数值分析在岩土工程中得到了广泛应用,为工程设计和施工提供了重要的支持和辅 助。
数值分析能够模拟和预测岩土工程中的各种复杂现象和问题,为工程实践提供了更 为精确和可靠的技术手段。
2021年浅谈各种数据分析方法在岩土工程中的运用比较
浅谈各种数据分析方法在岩土工程中的运用比较1. 1 边界元法的基本原理边界元法,顾名思义,就是只在边界上剖分单元,把边界分方程转变为线性代数方程,从而得出各边界单元处特定的边界值,然后再利用把边界值同域内数值联系起来的解析公式,这样就可以将计算区域内的任何一点的函数值求算出来。
按照边界元法的求算途径可分为两种类型: 直接法和间接法,前者利用具有明确物理意义的变量来建立边界积分方程; 后者却是利用不很明确的变量,一般,边界被加上了虚拟力和虚拟位移,这种虚拟力和虚拟位移是按照一定规律分布的,作为基本数,建立离散化的方程,待求出这些变量后再计算边界域内的位移和应力。
1. 2 边界元法在岩土工程中的应用在深基坑工程支护方式中,土钉墙越来越多得应用于工程中。
一般方法中,人们将主动区视为刚形体,用安全系数数值的大小来判断土钉墙的安全可靠性。
在实际的工作中,土钉墙主动区并非刚性体,它会随着开挖深度等因素的不同而发生大小不同的变形,所以,极限平衡法的应用有一定的局限性。
另外,有限元法也有应用在分析土钉墙稳定性的案例,用位移或应力值的大小来反映它的受力机制和状况。
虽然,在一定程度上,与极限平衡法相比,有限元法要相对合理一些,但是它和在工程中广为应用的安全系数直接联系起来较为困难,加上在计算过程中所需土体参数常规实验不易确定等因素,在工程中的应用也受到了一定程度的限制。
通过分析深基坑工程中土钉墙力学性状和对实际应用的分析,对提出的拥有土钉墙安全系数意义上的弹塑性边界元法的合理性给予了合理性证明。
工程实践证明,边界元法是较为理想的数值分析方法。
2 非线性数值分析在岩土工程中的应用下面本文将以 FLAC 为例对有限元非线性分析法在岩土工程中的应用进行简要分析。
FLAC,即快速拉格郎日差分分析,是力学计算的数值分析方法之一,主要应用于岩体力学中分析计算一般岩石的应力和应变。
它的基本原理类似于离散元的机构,但它使用于多种材料模式与边界条件的非规则区域的连续问题求解。
岩土工程的数值分析
有限元法的发展:20世纪60年代初,有限 元法在岩土工程中得到应用。由于它能够较 容易地处理分析域的复杂形状及边界条件、 材料的物理非线性和几何非线性,所以有限 元的应用和发展非常快。它在土体渗流、固 结、稳定和变形分析等各个领域得到广泛应 用。它被应用于浅基础、桩基础及各类深基 础、挡土墙、堤坝、基坑和隧道等各类岩土 工程问题的分析。有限元法不仅用于分析静 力问题、动力问题,还用于分析上部结构、 基础和地基的相互作用等问题。 在有限元 分析中可以采用总应力分析法,也可以采用 有效力分析法。
边界元法的发展: 20世纪20年代,边界元法的理论基 础就已初步形成。到20世纪60年代中期,边界元法在工 程技术问题中的应用逐渐斩露头角。尤其是自70年代末 直至现在,边界元法发展迅速。边界元法的思路是把所 要求解的微分方程转化成相应的边界积分方程,然后采 用边界积分方程的数值法求得原问题的数值解。这种方 法的特点是通过边界上的量来确定区域内部的未知量, 它与有限元法相比,具有信息准备工作少等优点。边界 元法采用类似于有限元法的离散技术来离散边界。离散 化所引起的误差仅来源于边界,提高了计算精度。边界 元法在分析边坡稳定性、地下水渗流等方面取得了良好 的效果。边界元法与有限元法相比,具有降低所求问题 的维数,计算量和计算时间相对减少,计算精度相对高 的特点,应该说比有限元法具有更强的生命力。但事实 上边界元法并没有得到像有限元法那样广泛的应用。其 中原因有,一是边界元法对变系数或非线性问题的适用 性不如有限元法;其次,边界元方程没有较好的前后置 处理技术,计算结果表示也不很直观和形象。
岩土数值分析 1-1
第1章 土的本构模型岩土工程的数值分析,离不开岩土材料的本构关系。
岩土工程数值分析的精度在很大程度上取决于所采用的本构模型的实用和合理性。
本构关系广义上是指自然界中作用于该作用产生的效用两者之间的关系。
如电学中的电压与电流的关系,力学中的力与变形之间的关系,热学中的温差与热流之间的关系,力学中水利梯度与渗流之间的关系。
土的本构关系主要是指土的力学本构关系,即土的应力-应变关系。
描述土的本构关系的数学表达式就成土的本构方程,或称本构模型。
一般将本构模型分为以下几类:弹性模型、弹塑性模型、粘弹塑性模型、内蕴时间性模型和损伤模型等。
要研究材料的塑性本构关系和塑性极限荷载,首先要建立材料产生屈服于破坏的条件与准则。
本章将介绍各种经典的和近年来提出的各种适于岩土类材料的屈服于破坏准则,并对于这些准则从理论与实践方法上进行评价。
这些准则包括区瑞斯卡屈服准则、米塞斯准则、摩尔-库仑屈服准则、辛克维兹-潘德准则、松港元提出的SMP 屈服准则及俞茂宏提出的双剪应力屈服与破坏准则。
在实际工程中,土体常常既有很复杂的应力应变特性,如非线性、弹性、塑性、粘性及剪胀性、应变硬化(软化)和各项异性等,同时应力路径、应力历史即土的状态、组成、结构和温度等均对其有不同程度的影响。
因此,为了反映土体真实的力学性状,必须建立较为复杂的本构模型。
而在实际工程应用中,在满足一定的精度条件下,又要求本构模型简单实用。
因此,目前的本构理论研究呈现两种倾向:一种是建立用于解决实际工程问题的本构模型;另一种则是比较精细的模型研究,目的在于进一步解释土体应力应变特性的内在规律。
1.1应力变力分析1.1.1应力张量在直角坐标系中,一点的应力状态可以用该点的九个应力分量表示⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=z zy zx yz y yz xz xy x ij στττστττσσ (1-1)式中,x σ、y σ、z σ为应力法向的分量;xy τ、xz τ、yx τ、yz τ、zx τ、zy τ为应力的切向分量。
岩土工程的数值方法 ppt课件
3.快速性
FLAC3D运行在所有的Window平台上,命令模式 可以用于标准输出窗口。FLAC3D提供内置的图元图形, 和快速、高分辨的图形能力来加速建模过程。求解参 数可以由用户指定,让用户可以在最大程度上控制持 续时间、长度和模型运行的效率。通过FLAC3D的强有 力的内建语言——FISH,用户可以进行附加的控制和 专用定制。
岩土工程的数值方法
第六章 岩土工程的数值方法
FLAC3D的应用介绍
一、关于FLAC3D
FLAC3D 是针对工程进行力学计算的三维 有限差分程序,它可以模拟由土体、岩体或其 它材料组成的结构体的三维表现行为,包括这 些结构体在满足屈服极限后发生塑性流动的行 为。
由于FLAC3D 的应用对象是地质岩土工程 , 它可以模拟各种不同的地质体材料及其各种不 同的行为模式,因而,广泛应用于斜坡稳定性、 基础工程、矿山开采、隧道工程等等岩土工程 中的数值模拟分析和工程设计。
算两种方式,但是,两种方式不能同时激活, 因为时步相差太大。蠕变和动态计算在动态/蠕 变分析中固定在一起。
(4)用户定义的基本模式: 用户定义的基 本模式现在可以用C++编写,并编译为DLL (动态连接库)文件,它们在FLAC3D 模拟需 要的时候加载上。DLL模式的运行与内置模式 的速度一样。用Visual C++ 6.0版本的编译器来 编译DLL文件。对所有FLAC3D C++ 模块的源 文件都提供给用户,相同的DLL 也被FLAC和即 将到来的3DEC 和UDEC所采用。这样,一个 用户定义的单个模块可以被好几个Itasca程序 利用。
4.检验性
FLAC3D已经应用了七年,应用于工程、咨询和 大学教育和研究之中。FLAC3D目前的使用许可遍布42 个国家,超过500用户,使它成为世界上应用最广泛 的岩土工程分析的三维数值建模工具。
数值分析在岩土工程中的应用
精选ppt
7
1、基本 概 念
数值分析的主要求解 方 法
变分法是讨论泛函的极值问题,对上述差分法及有限元法都可起推导 基本公式的作用,而这方法本身,也是数值方法中最古老的方法。
岩土工程中的基本方程包括土体平衡(或运动)方程、物理(或本构)方程、 几何方程、有效应力原理、孔隙流体(水)平衡方程、连续方程等。总控制 方程(即Biot动力固结方程)由这些基本方程组合而成。
总控制方程的推导基于以下假定: (1) 土体是完全饱和的横观各向同性弹性体。 (2) 土体的变形是微小的。 (3) 土颗粒和孔隙水不可压缩。 (4) 孔隙水相对于土骨架的渗流运动服从Darcy定律,其惯性力可不计。 (5) 应力应变的正负号法则与弹性力学相反。
离散单元法与其他数值方法(如有限单元法、边界单元法等)耦合更能发 挥各自方法的优点。例如,用边界单元法考虑远场应力的影响以模拟弹性的 性质,用有限单元法作为中间过渡考虑塑性变形,再用离散单元法考虑近场 不连续变形的情况,从而极大地扩展了数值方法的解题范围。
精选ppt
13
2、几种常见数值分析方法的主要 特 点
正确的剖分
不正确的剖分
精选ppt
10
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method (FEM)
研究表明,对于大多数岩土工程问题,无论是进行总应力分析还是进 行有效应力分析,均可归结为对Biot动力固结方程的求解。因此,可将Biot 动力固结方程作为岩土工程问题的总控制方程。
图中数值方法列出最常用的五种:差分法、有限无法、边界元法、变 分法和加权余量法。
关于岩土工程的数值计算方法的综述
关于岩土工程的数值计算方法的综述学院:资源与土木工程学院专业:岩土工程学号:姓名:数值计算方法其主要有有限单元法、有限差分法、边界元法、离散元法和流形元法等。
有限单元法:有限单元法发展非常迅速,至今已经成为求解复杂工程问题的有力工具,并在岩土工程领域广泛的采用,主要的分析软件ANSYS。
有限单元法的最基本的元素是单元和节点,基本计算步骤的第一步为离散化,问题域的连续体被离散为单元与节点的组合,连续体内部分的应力及位移通过节点传递,每个单元可以具有不同的物理特征,这样,便可以得到在物理意义上与原来的连续体相近似的模型。
第二步为单元分析,一般以位移法为基本方法,建立单元的刚度矩阵。
第三步由单元的刚度矩阵集合成总体刚度矩阵,并由此建立系统的整体方程组。
第四步进入计算模型的边界条件,求解方程组,求得节点位移。
第五步求出各单元的应变、应力及主应力。
有限差分法:有限差分法在岩土工程中是应用非常广泛的方法,在数值计算模拟上有很大的贡献,主要的应用软件为FLAC3D。
基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解。
然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
边界单元法:边界单元法在岩土工程领域也有很大优势,主要的应用软件是二维边界元法软件THBEM2和三维边界元法软件THBEM3,它们在复杂工程问题的线弹性应力分析以及弹性力学辅助教学等方面的应用有很大优势。
积分法统称为边界单元法,有直接法和间接法两类,它们都是利用了简单奇异问题的解析解,并可近似满足每个边界单元的应力和位移边界条件。
该法仅仅限定和离散问题的边界,可把问题的重点转移到边界上,可以有效地使已知条件降维,从而减小方程组的规模,大大提高计算效率。
岩土工程数值计算方法
岩土工程数值计算方法
岩土工程数值计算方法牛不牛?那绝对超厉害!咱先说说这步骤哈。
首先得收集岩土工程的各种数据,就像大厨准备食材一样,一点都不能马虎。
然后建立数学模型,这就好比给房子搭框架,得结实。
接着进行计算求解,这过程就像赛车冲刺,紧张又刺激。
注意事项可不少呢!数据得准确呀,要是数据错了,那不就像在沙漠里找大海,瞎忙活嘛!模型选择也得合适,不然就像穿小鞋走路,难受得很。
再说说安全性和稳定性。
这可太重要啦!要是不稳定,那不是像在摇摇欲坠的桥上走,提心吊胆嘛!所以在计算过程中一定要确保结果的可靠性,不然出了问题可不得了。
应用场景那可多了去了。
比如在建筑工程中,可以预测地基的沉降,这就像给大楼安了个保险。
在隧道工程中,能分析围岩的稳定性,就像给隧道穿上了铠甲。
优势也很明显啊,省时省力还精准,比起传统方法,那简直是鸟枪换炮。
举个实际案例,有个大型建筑项目,用了岩土工程数值计算方法,提前预测了各种问题,及时调整方案,最后顺利完工。
这效果,杠杠的!
岩土工程数值计算方法就是这么厉害,能解决实际问题,让工程更安全、更高效。
咱就该大胆地用起来,让它为我们的工程建设助力。
基于岩土工程数值方法应用
大坝稳定性数值模拟
01
02
03
大坝应力分析
通过数值方法模拟大坝在 不同工况下的应力分布, 评估大坝的整体和局部稳 定性。
大坝变形分析
分析大坝在不同水位下的 变形情况,预测大坝可能 出现的裂缝和滑动等失稳 现象。
大坝抗震性能评估
利用数值方法模拟地震作 用下大坝的响应,评估大 坝的抗震性能和安全性。
边坡稳定性数值模拟
详细描述
岩土工程涉及岩石、土壤和地下水等自然资源的工程应用,主要研究岩石和土的 物理性质、力学行为和工程应用。它具有地质性、材料性、力学性等特点,需要 综合考虑地质学、材料科学、力学等多个学科的知识。
岩土工程的重要性
总结词
岩土工程在国家建设、经济发展、环境保护等方面具有重要 意义,是土木工程领域的重要分支。
详细描述
离散单元法在岩土工程中常用于分析岩石崩塌、滑坡等问题。该方法能够真实地模拟岩土材料的破裂和块体运动 ,但计算效率较低,需要处理大量的离散单元。
无单元法
总结词
无单元法是一种不依赖于网格的数值分析方 法,能够避免网格生成和扭曲问题。
详细描述
无单元法在岩土工程中常用于分析动态破坏 、冲击波传播等问题。该方法具有较高的计 算效率和鲁棒性,但需要处理无网格节点的 插值和积分问题。
边坡应力分析
通过数值方法分析边坡在不同工况下的应力分布 ,评估边坡的整体和局部稳定性。
边坡变形分析
分析边坡在不同荷载下的变形情况,预测边坡可 能出现的滑坡和崩塌等失稳现象。
边坡加固方案优化
根据数值分析结果,优化边坡加固方案,提高边 坡的稳定性。
地下水渗流数值模拟
地下水流动分析
通过数值方法模拟地下水 在土体中的流动规律,预 测地下水位的变化趋势。
岩土工程数值分析
土的本构模型
有效应力(也称应力强度、或广义剪应力)
3J2
1 2
(1 2 )2 ( 2 3)2 ( 3 1)2
在简单拉伸时,应力强度还原为简单拉应力
1
1 2 0 故
有效剪应力(也称剪应力强度)
T
J2
1 6
(1 2 )2 ( 2 3)2 ( 3 1)2
土的本构模型
应变不变量
II12''
x y z 1 2 2
x
y
y z
z
x
2 xy
V
2 yz
2 zx
1 2
23
31
I
' 3
1 2 3
偏应变不变量
JJ12''
6(
2 xy
2 yz
2 zx
)
1 6
(1 2 )2 ( 2 3 )2 ( 3 1)2
2 3
(12
2 2
2 3
)
J
3
SxSySz
2 xy
yz zx
S
x
2 yz
S
y
2 zx
S
z
2 xy
S1S 2 S3
弹塑性本构关系中,J 2 J 3 反映切应力大小及方向。
破坏比:
R
(1 3 ) f (1 3 )ult
岩土工程中的数值方法
中 图分 类 号
T 42 U 5
文献标识码
A
1 引 言
元 法 ( F M)非 连 续 变 形 分 析 ( D 1流 形 Co g 先 将有 限单 元法 用 于 土 坝 的 应 力 RE 、 D A、 l h首 u
方 法 ( n od Me o )无 单 元 法 ( sl s 和 变 形 分析 . 内沈 珠 江 1 7 Mai l t d 、 f h Mehe s 国 9 7年 首 先 将 有
工程 中 , 土体一般是 非连续 、 均匀 、 岩 非 各 形 分 析 方 法 两 大 类 。 连 续 变 形 分 析 法 主 要 以 及 材 料 的非 均 匀性 和各 向 异 性 .还 可 以 向异 性 的 介 质 。 存 于 自重 应 力 场 、 造 应 有 : 限 差 分 法 、 限单 元 法 、 限单 元 法 有 效 地 模 拟 材 料 的 非 线 性 应 力 应 变 关 系 . 赋 构 有 有 无
题 , 数 值 方 法 成 为 岩 土 力 学 f 题 的 重 要 研 究 手段 。 绍 了岩 土 工 程 中九 种 数 值 方 法 的特 点 . 论 了 当前 岩 土数 使 - - I 介 讨 值 方 法存 在 的 问题 和 发 展 前 景 。
关键 词 岩 土 工 程 数值方法 特点 发 展 前 景
有 限元 法 为 了推 广有 限单 元 法 .使 其
了解 , 同时 往 往 耗 费 巨大 , 且受 地形 、 质 、 容 易 编 制 程 序 , 2 地 从 0世 纪 4 0年 代 以 来 , 至 能 够 处 理 简 单 的非 连 续性 问 题 。人 们 提 出 施 工 条 件 限制 。 验 结 果 不 具 代 表 性 , 难 今 仍 被 广 泛 应 用 。4 试 很 0年 代 后 期 . 工 中 的 了 各 种 能 反 映非 连 续 性 质 的 简 单 力 学 元 件 土 推 广 到 其 它 工程 。 而 岩 土 力学 的 理 论研 究 渗 流 及 固 结 等 实 际 问 题 被 差 分 法 成 功 解 或 特 殊 界 面 单元 用 于 模 拟 单 一 非 连 续 界 面
岩土工程中的数值方法综述
对 于岩 土工程数值 分析 方法 按照其 特点 可 以分 为连续 变 形 应用 问题 。 无 网格伽辽金 法 的前 处 理只 要节 点 位 置信 息 , 不 用 网格 信 分析和非 连续变形分析两大类方法 , 主要包括 : 1 ) 有 限差分法 ; 2 ) 息, 比有限元法简单 。 边界元法 ; 3 ) 无网格方 法 ; 4 ) 耦合 方法 等。 由于数值 分析 方法 能 无 网格伽 辽金法计算的结果是光 滑连续 的, 不必再 进行应 力 够在较短 时间内完成大量的计算分 析工作 , 计算 成本低 、 效 率高 , 光顺化等 后处理。 因此在 大型水 利 、 土木 、 矿山、 隧道 、 交通 等工 程 的设计 分 析 中起 但是 , 它在 以下几 个方 面还存 在着 需要 完善 的地 方 , 极大 地 到 了重要 的作用 。随着数值分析方法 的不 断发展 , 也为 岩土工 程
比较 了工程应用 中有 限元 法的劣势及无 网格法 的优势 , 为岩土工程 的数值分析 与模 拟提供参考。 关键 词 : 岩 土工程 , 数值方 法 , 计算 , 优 势
中图分 类号 : T U 4 5 2 文献标识码 : A
1 岩 土 工程 中常用 的数值 分 析方 法
的函数 系列 , 适用 于分析 各类 具有 高梯度 、 奇 异性 等特殊 性 质 的
与其他数值分析 方法 , 比如有 限元法 的联 合使 用 , 充分 发挥 各 自 法的优点并进行耦合 提高。 这些数值方法 中 , 以有限元法为代 表的连续 变形 分析方法 在 的优势 , 提高 了计 算效率和计算精度 。
1 . 2 无 网格 法 无 网格方法是在数值 计算中不需要 生成 网格 , 而是按 照一些 任 意分布的坐标点构 造插值函数离散 控制方程 , 就可 方便 的模 拟
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
E p 0
有限元法的基本方程
单元位移函数
(i,1) | ( x1 , y1 ) | (u1 , v1 )
( j,2) | ( x2 , y2 ) | (u2 , v2 )
假设任意点的位移:
(k ,3) | ( x3 , y3 ) | (u3 , v3 )
第7章 岩土工程数值分析方法
有限元法
边界元法
有限差分法
离散单元法
概述
1 有限元法
基本思路:将复杂的结构看成由有限个仅在结点处
联结的整体,首先对每一个单元分析其特性,建立 相关物理量之间的相互联系。然后,依据单元之间 的联系再将各单元组装成整体,从而获得整体特性 方程,应用方程相应的解法,即可完成整个问题的 分析
A Sa
Wi dxdy
T A
We Wi
最小势能原理
定义1:外力从位移状态退回到无位移的初始状态时所作的功 称为外力势能
T T E* ( F ldxdy F p b S lds) A Sa
定义2:形变势能和外力势能的和称为总势能
1 E p T dxdy ( FbT ldxdy FST lds ) A Sa 2 A 最小势能原理:实际发生的位移总能使对应l的势能一阶变分
具有一定厚度的单元可按四边形等参元处理
k B DBdV B DB J dd
e T 1 1 T V 1 1
当厚度很小时按等厚度或无厚度节理单元处理
s K s n 0
0 us Kn u n
多节理岩体的模拟
等效连续体
利用节理单元来模拟密集分布或随机分布的节理与裂隙是不 适宜的,会给离散化与计算带来诸多困难与麻烦 从总体上考虑节理裂隙对岩体的影响,将岩体视为等效的正 交异性、各向异性或各向同性体,即等效连续体 将岩体中的节理裂隙当成存在于岩体材料内的一种损伤,因 如损伤力学的原理建立损伤模型来考虑节理裂隙的影响
N i x Bi 0 N i y 0 N i y N i x
B2 Bm
N i N i x 1 N J N i i y
或: u
[ N ] v
N1 [N ] 0 0 N1 N2 0 0 N2 N3 0 0 N3
T
u1
v1 u2 v2 u3 v3
插值函数(形函数)
N1 [N ] 0
0 N1
N2 0
0 N2
N3 0
0 N3
K
1 3 L Kn 4
2K n 0 Kn 0 Kn 0 2K n
(对 2K s 0 2K s 0 2K s 0 2K n 0 2K n 0 Kn 2K s 0 Ks 0
称)
2K n 0 Kn 2K s 0
2K n
变厚度节理单元
* T T * T T * T An A An
F F k
* T e * T * T e * T e e
T e B D B tdA A n
单元刚度矩阵:
k B DBtdA B DBt
e T T An
总体刚度矩阵
由于虚位移
的任意性,等式两边与其相乘的矩阵相等,则:
*
F k 设结构体剖分成n个单元,根据虚功有:
e e
n * T e n
F k
* T e e i 1 i 1
K U P
N1 x B 0 N 1 y
0 N1 y N1 x
N 2 x 0 N 2 y
0 N 2 y N 2 x
N 3 x 0 N 3 y
0 N 3 B1 y N 3 x
B2
分析过程:
结构离散化 集成总体特性 确定单元位移模式 接方程求未知量 单元特性分析
概述
工程界比较流行、被广泛使用的大型有限元软件:
MSC(航空航天领域,是目前规模最大的有限元分析 系统);Marc(非线性分析软件);Adina(可进行结构、 流体、热的耦合计算,具有隐式和显式两种时间积 分算法,非线性运算功能强大);ANSYS(可进行 结构、流体、热、电磁场的计算,是有限元分析的 通用软件包)
0 us K n un
y
n
s s
n
x
节理单元本构模型
Fn
Fs
kn tan
Fs max
ቤተ መጻሕፍቲ ባይዱ
ks c n tan
0
un
0
us
节理单元刚度矩阵(假定位移沿单元长度线性变化)
2K s 0 Ks tl 0 e k 6 Ks 0 2 K s 0
i
空单元
i 1
v
施工过程模拟
空单元
第1步开挖
h H
空单元
第2步开挖
v H
空单元
第3步开挖
节理及不连续面的模拟
平面问题节理单元-Goodman
Goodman单元是无厚度4结点单元
结点传递切向力与法向力
节理应力-应变关系
s K s n 0
2K n 0 Kn 0 Kn 0 2K n
(对 2K s 0 2K s 0 2K s 0 2K n 0 2K n 0 Kn 2K s 0 Ks 0
称)
2K n 0 Kn 2K s 0
2K n
考虑嵌入的节理单元模拟(考虑转动)
初始地应力场
自然状态的岩体处于一定的初始地应力状态,在结构荷载作 用下,岩体内的应力为荷载产生的应力与初始地应力之和 由于初始地应力的存在,开挖将导致部分岩体卸荷,通常采 用沿开挖面作用着与地应力等价的“开挖释放荷载”
释放荷载
施工建造过程的模拟
开挖释放荷载
i 1
h
Phi
Pvi
总体刚度矩阵
K 总体刚度矩阵 U 结点位移列阵 P 荷载列阵
等参元分析
平面任意四边形单元
4
3
0
1
2
结点位移矩阵:
u1
v1 u2 v2 um vm
T
插值函数:
N1 [N ] 0
0 N1
N2 0
0
Nm 0
N2
0 Nm
s K s n 0 M 0 0
2K s 0 Ks 0 tl e k 4 Ks 0 2 K s 0
0 Kn 0
0 u s 0 u n K
ai bi x ci y N i ( x, y) ( x1 y2 x2 y3 x3 y1 ) ( x2 y1 x3 y2 x1 y3 )
形函数特点:
N ( x , y ) 1, i 1,2,3 i i i ,i j N i ( x j , yi ) 0, i 1,2,3
岩土工程涉及无限域或半无限域,但处理问题时只能对有限 域进行离散化 模型范围可取结构体轮廓尺寸的3~4倍
边界效应
以小变形理论为基础的有限元法中,力与变形的影响范围是 无限域,因此,设定有限域,并假定模型边界的位移为零或 为受力边界就会带来误差,靠近边界越近误差越大,靠近边 界越远误差越小
初始地应力场与释放荷载
1 N1 ( , ) (1 )(1 ) 4 1 N 2 ( , ) (1 )(1 ) 4
1 N 3 ( , ) (1 )(1 ) 4 1 N 4 ( , ) (1 )(1 ) 4
几何矩阵:
B B1
单元应力矩阵:
S DB DB1
B2
B3 S1 S2
S3
Si DBi
单元刚度矩阵
岩土体或结构体发生虚位移,单元结点的虚位移为 应变为 * ,则根据虚功原理有:
,相应的虚
*
N F tdA N P tds tdA
优点:可分析几何形状及受荷条件复杂、非均质的
各种实际结构;可在计算中模拟各种复杂的材料本 构关系、边界条件等;前后处理技术先进。
有限单元法的理论基础
虚位移原理:受给定外力的变形体处于平衡状态的充
要条件是,对一切虚位移,外力所作总虚功恒等于内 力总虚功
We FbT ldxdy FST lds
i 1,2,3, , m
单元刚度矩阵:
k B DBtdA
e T An
1 1 m 1
B DBt J dd
T m 1 T i
1 1
(B D B i t J ) lpWlW p
模型范围与边界效应
模型范围
u ( x, y ) a1 a2 x a3 y v( x, y ) a4 a5 x a6 y
单元位移函数: u ( x, y )