第六篇 ：地应力平衡方法以及注意事项
地应力及其测量解析PPT教案
b)传统的岩石工程的开挖设计和施工的经验类比法有局限性。 小规模和地表的岩石工程，经验类比的方法往往是有效的。 大规模和深部的岩石工程，经验类比法有不足。
第3页/共74页
4
c) 是岩石工程数值分析方法计算分析的必要 前提条件。
第10页/共74页
4.1.3 地应力的成因
3）岩浆侵入
岩浆侵入、冷凝收 缩和成岩，均产生 应力场。
高温岩浆对周围岩 体施加均匀压力， 冷凝后收缩，并从 接触面逐渐向内部 发展。
产生局部应力场。
第11页/共74页
12
4.1.3 地应力的成因
4）岩体自重应力场
垂直应力： z H
侧压力： x y z
向接近水平的，而且最大水平 主应力一般第为5页/共垂74页直应力的1～2 倍，甚至更多。
6
4.1.3 地应力的成因
1）地幔热对流 硅镁质组成的地幔因高温，
上下对流、蠕动，深部地幔 上升到顶部时变成2股方向相 反的平流：与反向平流相遇， 转为下降流进入深部，形成 封闭的循环体系。 亚洲形成孟加拉湾延伸至贝 尔加湖的最低重力槽，是拉 伸的袋装区，西昌、攀枝花 到昆明的裂谷刚好处于这一 地区，形成以西藏为中心的 大对流环（上升），华北-山 西地堑有下降流。 产生很大的水平应力。
2)测量技术和计算理论 比 较成熟
3)精度相对较高的地应 力
4)目前使用最为广泛
是国际岩石力学学会推 荐使用的地应力测量方 法
第45页/共74页
46
c.测量要求与步骤
1)测量要求 钻孔的深 度必须超 过开挖 影响区， 才能测到 岩体内的 原始应力 ，否则测 出的是二 次应力。
地应力平衡总结
地应力平衡1、地应力平衡好坏评判标准1)地应力平衡后，位移云图中最大位移达到10-6量级或更低（接近于0）。
（主要判别条件）2)地应力平衡后，应力云图中应力有一定的数值。
（也就是应力不为0，但变形接近于0）2、进行地应力平衡的原因总的来说，如果不进行地应力平衡，而只施加重力，模型会在重力作用下产生变形，而实际工程中，我们施加荷载时，重力产生的而变形已经产生，实际上得到的是附加应力产生的变形。
1)我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
这就是所谓“地应力平衡”的含义、目的、作用。
2)地应力平衡中的外力和内力的问题。
地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
地应力测量方法及其需要注意的问题
地应力测量方法及其需要注意的问题地应力是指存在于地壳中的内应力。
主要由重力应力、构造应力、孔隙压力、热应力和残余应力等耦合而成,重力应力和构造应力是地应力的主要来源。
地应力测量是确定工程岩体力学属性、进行围岩稳定性分析、实现岩土工程开挖设计和决策科学化的前提。
地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用,近年来随着我国社会经济的持续快速发展,我国水电领域工程建设保持着较快增长势头,工程建设地点向江河源头、高山峡谷地带延伸,工程建设内容往往包含深埋长深隧道,大跨度、大尺度地下厂房等,在这种情况下,我国地应力测试事业也取得了长足的进步,各种试验手段、测试方法层出不穷,并取得一定的成果。
1地应力测量方法1.1 应力解除法应力解除法是以弹性理论为基础,它把一定范围内的岩体视为均质的、各向同性的完全弹性体。
这一测量方法的实质是在被测虚力场的岩体中选定测点,在测点位置安设测量元件,然后在所安装的测量元件周围掏槽或套孔,使安设有测量元件的岩石与周围岩体分离,也就是使这一部分岩石从被测应力场作用之下解脱出来。
此时,测点岩石将由于外力的消失而产生弹性恢复变形。
通过测量元件将这一变形记录下来,即可按弹性理论来确定被测应力场的3个主应力的大小、方向和倾角。
应力解除法测量地应力的方法有:孔底应变计、孔径应变计、孔壁应变计、空心包体应力计等方法,其中孔底应变计、孔径应变计只能测出二维应力,若用它测三维应力,则需要打交于一点互不平行的三个钻孔。
采用孔壁应变计和特殊制作的空心包体式孔壁的应力计只需要打1个钻孔就可测出三维应力。
1949年奥尔森(O.J.Olson)第一次将应力解除法用于岩石应力测试以来,套孔应力解除法发展为技术上比较成熟的一种原岩应力测量方法。
套孔应力解除法具有测量灵敏度高、测量结果可靠、可以在深孔中进行测量测点的三维应力状态(需要利用三孔交汇的方法)等特点。
因此,利用套孔应力解除法可以较为准确地测量矿山岩体的原岩应力。
第六篇 ：地应力平衡方法以及注意事项
第六篇：地应力平衡方法以及注意事项注意：只有采用弹塑性本构模型时需要地应力平衡，弹性本构不需要地应力平衡！第一部分地应力平衡方法简介地应力平衡主要有五种方法：（1）自动平衡：第一步创建分析步geostatic ，这种方法注意只能在第一步只能有土和重力的情况下能使用，有其他部件或者接触时计算不能收敛，效果是最好的，方便简单！（2）*initial conditions,type=stress,geostatic该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。
所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。
计算速度快,收敛性好。
缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。
高版本在CAE里也能操作用计算器算出每个分界面上的应力和坐标对应填入，也比较方便不需要修改关键字（3）*initial conditions,type=stress,geostatic,user导入ODB里的方法，也比较简单，高版本可在截面上操作，不需要修改关键字你放入ODB后，填入第一步不需要填名称就是1 ，增量步就是你第一步计算的最后一个增量步（4）*initial conditions,type=stress,input=FileName.csv(或inp)该方法中的文件FILENAME.INP获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。
但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。
地应力基本概念及测量方法
地应力基本概念及测量方法应力等因素导致岩体具有初始地应力（或简称地应力）是最具有特色的性质之一。
就岩体工程而言，如不考虑岩体地应力这一要素，就难以进行合理的分析和得出符合实际的结论。
岩体应力天然应力是指未经人为扰动的，主要是在重力场和构造应力场的综合作用下，有时也在岩体的物理、化学变化及岩浆侵入等的作用下所形成的应力状态，称为岩体天然应力或岩体初始应力，有时也称为地应力。
天然应力构成：岩体自重自重应力构造运动构造应力流体作用静水压力梯度，渗流应力其他（低温、地球化学作用）地壳岩体的天然应力状态与人类的工程活动关系极大，它不仅是决定区域稳定性的重要因素，而且往往对各类建筑物的设计和施工造成直接的影响。
比如，地下空间的开挖必然使围岩应力场和变形场重新分布并引起围岩损伤，严重时导致失稳、垮塌和破坏。
这都是由于在具有初始地应力场的岩体中进行开挖所致，因为这种开挖荷载通常是地下工程问题中的重要荷载。
由此可见，如何测定和评估岩体的地应力，如何合理模拟工程区域的初始地应力场以及正确和合理地计算工程问题中的开挖荷载，是岩石力学与工程问题中不可回避的重要问题。
已有的研究和工程实践表明，浅部地壳应力分布主要有如下的一些基本规律：地应力是一个具有相对稳定性的非稳定应力场，它是时间和空间的函数。
实测垂直应力基本等于上覆岩层的重量。
水平应力普遍大于垂直应力。
平均水平应力与垂直应力的比值随深度增加而减小，但在不同地区，变化的速度很不相同。
最大水平主应力和最小水平主应力也随深度呈线性增长关系。
最大水平主应力和最小水平主应力之值一般相差较大，显示出很强的方向性。
地应力的上述分布规律还会受到地形、地表剥蚀、风化、岩体结构特征、岩体力学性质、温度、地下水等因素的影响，特别是地形和断层的扰动影响最大。
高应力区实践表明，在高应力区，地表、地下工程施工期间所进行的岩体开挖工作，往往能在岩体内引起一系列与卸荷回弹和应力释放相联系的变形和破坏现象，其结果是不仅会恶化地基或边坡岩体的工程地质条件，而且作用的本身有时也会对建筑物造成直接的危害。
附录2-ADINA中关于地应力的四种处理方法
附录二：ADINA中关于地应力的处理方法方法一、导入初始应力方法此种方法的特点是针对线性材料；或非线性材料的线弹性阶段才能采用，如果非线性材料的部分单元进入塑性，则不能采用；非线性弹性材料、用户开发的材料则不能使用。
基本过程是将施加重力得到的应力从结果文件（Por）提取，通过前处理的Initial Condition 窗口中施加到模型上。
具体步骤如下：1.先对模型加重力载荷进行计算，计算完毕后进入后处理模块，用List Zone（Whole Model）列出所有节点的应力值；（注意：列出值时要选择smooth，否则没有结果，最后一步结果（Single Response）。
因为应力结果是单元的结果。
选择六个分量的顺序对于3Dsolid单元应该是stress11,stress22,stress33,stress12,stress13,stress23），对于2D solid，plate和shell单元初始应力轴的1-2轴必须位于平面内，3轴位于平面外，所以这些单元的六个分量的顺序为stress22,stress33,stress11,stress23,stress12,stress13。
如下图所示。
2.列出结果后点击Export将结果写到给定名字的TXT文件，之后将文件处理成ADINA的数据表格式，过程基本是用替换方式去掉所有的Node+空格；去掉多余文件说明；然后用Excel文件将连续空格转化称为制表符；然后输出新的TXT文件；见如下图示的过程说明：最后在文件第一行加上头信息或空行即可；3.重新打开原来的数据库，施加初始应力条件Model->Initial Condition，选择strain选项（这里没有Stress选项），然后用Import读入这个文件即可；4.在Control->Miscellaneous options中选择Input Strain use as Initial Stress cause deformation选项5.在所有的单元组中指定节点单元选择具有初始应变选项；6.现在可以计算。
abaqus地应力平衡的步骤
abaqus地应力平衡的步骤一、前期准备在进行地应力平衡分析之前，需要进行一系列的准备工作。
这包括收集地质勘察资料、了解场地工程地质条件、确定分析目标和评估范围等。
此外，还需要进行必要的模型简化与假设，以便在有限元分析中准确模拟实际地质结构和应力分布。
二、建立模型在abaqus中建立地应力平衡模型通常涉及以下步骤：创建模型空间、设置网格密度和类型、导入地质勘察数据（如岩石层分布、岩土性质等）以及为不同的材料设置相应的物理性质（如弹性模量、泊松比等）。
此外，还需根据实际工程需求和地质条件，设置边界条件和载荷，以模拟地质结构和应力分布。
三、设置材料属性材料属性的设置对于地应力平衡分析至关重要。
需要根据地质勘察资料和试验数据，为不同的岩土层设置准确的物理性质，包括密度、弹性模量、泊松比以及摩擦角和内聚力等。
这些参数将直接影响模型的计算结果和应力分布状态。
四、地应力平衡分析在进行地应力平衡分析时，需要采用合适的求解器和算法，以确保计算效率和准确性。
abaqus提供了多种求解器选择，如静力分析、动力分析以及流体分析等。
在地质工程领域，静力分析是最常用的求解器，用于模拟岩土体的稳态应力分布状态。
在设置求解器之后，需要设定合适的迭代步长和收敛准则，以便模型在迭代过程中逐渐达到平衡状态。
五、后处理与结果分析地应力平衡分析完成后，需要利用abaqus的后处理功能对结果进行详细分析。
这包括绘制应力分布云图、查看节理和断层的应力状态、评估潜在的破坏区域等。
通过对结果的深入分析，可以评估地质结构的稳定性和安全性，为进一步的地质工程设计和优化提供依据。
六、优化与迭代根据后处理阶段的结果分析，可能需要对模型进行优化和迭代。
这包括调整网格密度和类型、改进边界条件和载荷设置、调整材料参数等。
通过不断优化和迭代，可以提高模型的精度和可靠性，从而更准确地模拟实际地质结构和应力分布状态。
七、注意事项在地应力平衡分析过程中，需要注意以下几点：1.确保收集足够的地质勘察资料和试验数据，以便为模型提供准确的材料属性和边界条件。
地应力及地应力测量方法计划简介
地应力与地应力丈量方法简介地应力，又称原岩应力，也称岩体初始应力或绝对应力，是在漫长的地质年月里，因为地质结构运动等原由产生的。
在一准时间和必定地域内，地壳中的应力状态是各样发源应力的总和。
主要由重力应力、结构应力、孔隙压力、热应力和节余应力等耦合而成 ,重力应力和结构应力是地应力的主要根源。
地应力的形成主要与地球的各样动力运动过程相关，此中包含：板块界限受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非平均扩容等。
此外，温度不均、水压梯度、地表剥蚀或其余物理化学变化等也可惹起相应的应力场。
而重力作用和结构运动是惹起地应力的主要原由，此中尤以水平方向的结构运动对地应力的形成影响最大。
地应力丈量，就是确立拟开挖岩体及其四周地区的未受扰动的三维应力状态，这种丈量往常是经过多个点的量测来达成的。
地应力丈量是确立工程岩体力学属性、进行围岩稳固性剖析、实现岩土工程开挖设计和决议科学化的前提。
地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用，所以地应力研究是目前国际采矿界上的一个前沿性课题，近几十年来，世界上很多国家均展开了地应力的丈量及应用研究工作 ,获得了众多的成就。
跟着矿区开采现代化进度的不停提升和开采深度的不停增添，对矿区所处的地质条件和应力环境提出了更进一步的要求。
查明矿区深部煤炭资源的开采地质条件和应力环境，为深部矿井的设计、建设和生产供给更为精美靠谱的地质资料和数据，以便采纳有效技术手段和举措，防备和减少灾祸的发生，是实现矿井安全高效生产的重要保障。
地应力是惹起采矿工程围岩、支架变形和损坏、产生矿井动力现象的根本作使劲，在诸多的影响采矿工程稳固性要素中，地应力是最重要和最根本的要素之一。
正确的地应力资料是确立工程岩体力学属性，进行围岩稳固性剖析和计算，矿井动力现象地区展望，实现采矿决议和设计科学化的必需前提条件。
采矿规模的不停扩大和开采深度的纵深展开，地应力的影响越加严重，不考虑地应力的影响进行设计和施工常常造成露天边坡的失稳、地下巷道和采场的坍塌损坏、冲击地压等矿井动力现象的发生，以致矿井生产没法进行，并常常惹起严重的事故，造成人员伤亡和财富的重要损失。
地应力及井壁稳定机理讲课材料
H
h
来源 渤南井稳研究报告 王家祥 QHD32-6 研究报告 刘畅 刘建中 刘建中 刘建中 刘建中 根据破裂压力
2.1-2.2 2.32 1.96 2.23 2.10 2.10 2.10 2.10 2.25
1.8-2.2 2.92 2.19 2.35 2.59 2.97 2.59 2.74 2.10
v
H h
108°52′
278000
F2 中块3井区 F1
FA 中块4井区
F2A
20°50′ 20°49′
2305000 2304000 2303000
南
108°53′
280000
块
108°55′
282000 284000
0
0.5
1.0 km
最大水平主地应力方向N100E左右
108°54′ 108°56′
涠洲12－1－1井地层倾角测井数据统计表
井段（米） 1447 －1462 1462 －1486 1493 －1500 1560 －1567 1570 －1574 1585 －1675 1690 －1790 1947 －1950 1996 －2000 2285 －2295 2650 －2675 2805 －2890 一三臂井径 （in） 二四臂井径 （in） 井径差值（in ） 一号臂方位 （度） 12.3 12.3 12.3 12.3 12.3 12.3 12.2 14.5 16.5 17.5 14.8 14.2 14.5 13.8 14.2 14.1 14.1 15.2 14.8 12.3 12.3 12.3 12.3 12.3 2.3 1.5 1.9 1.8 1.8 2.9 2.6 2.2 4.3 5.2 2.5 1.9 230 －240 230 －250 240 －250 230 －235 230 －250 225 －235 220 －230 45 －55 50 －60 225 －235 220 －235 220 －240
地应力及其分布规律
地应力及其分布规律1 、地应力的基本概念地应力是存在于地层中的未受工程扰动的天然应力，也称岩体初始应力、绝对应力或原岩应力。
广义上也指地球体内的应力.它包括由地热﹑重力﹑地球自转速度变化及其他因素产生的应力。
地应力是各种岩石开挖工程变形和破坏的根本作用力;是确定工程岩体力学属性，进行围岩稳定性分析,实现开挖设计和决策科学化的必要前提条件。
此外地应力状态对地震预报、区域地壳稳定性评价、油田油井的稳定性、核废料储存、岩爆、煤和瓦斯突出的研究以及地球动力学的研究等也具有重要意义。
2、地应力的成因产生地应力的原因是十分复杂的，地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆浸入和地壳非均匀扩容等。
另外，温度不均、水压梯度、地表剥蚀或其它物理化学变化等也可引起相应的应力场。
其中，构造应力场和自重应力场为现今地应力场的主要组成部分。
当前的地应力状态主要由最近的一次构造运动所控制,但也与历史上的构造运动有关.由于亿万年来，地球经历了无数次大大小小的构造运动，各次构造运动的应力场也经过多次的叠加、牵引和改造,另外，地应力场还受到其他多种因素的影响,造成地应力状态的复杂性和多变性,地应力成因之一：地幔热对流（图１、图2）地应力成因之一：板块边界受压（图3)地应力成因之一:岩浆浸入(图4）３、地应力的影响因素地壳深层岩体地应力分布复杂多变,造成这种现象的根本原因在于地应力的多来源性和多因素影响,但主要还是由岩体自重、地质构造运动和剥蚀决定。
1)岩体自重的影响岩体应力的大小等于其上覆岩体自重,研究表明：在地球深部的岩体的地应力分布基本一致。
但在初始地应力的研究中人们发现，岩体初始应力场的形成因素众多，剥蚀作用难以合理考虑,在常规的反演分析中，通常只考虑岩体自重和地质构造运动2)地形地貌和剥蚀作用对地应力的影响地形地貌对地应力的影响是复杂的,剥蚀作用对地应力也有显著的影响，剥蚀前,岩体内存在一定数量的垂直应力和水平应力，剥蚀后，垂直应力降低较多，但有一部分来不及释放，仍保留一部分应力数量,而水平应力却释放很少，基本上保留为原来的应力数量，这就导致了岩体内部存在着比现有地层厚度所引起的自重应力还要大很多的应力数值。
6-1 地应力及其原理 (2)
水平方向的构造运动对地壳浅 以上表明： 层地应力的形成起控制作用。
平均水平应力/垂直主应力(%)
国家名称
中国 澳大利亚
加拿大 美国 挪威 瑞典 南非 前苏联 其他地区
<0.8 32 0 0 18 17 0 41 51 37.5
0.8~1.2 40 22 0 41 17 0 24 29 37.5
>1.2 28 78 100 41 66 100 35 20 25
动和变化中，由于重力和地壳中构造运动等因素的影
响，原岩在空间上存在着按一定规律分布的初始应力 场，称为原岩应力场(天然应力场）或地应力场。 --基本概念
3、人类的工程活动将会扰动原岩内地应力场的自然平
衡状态，使一定范围内的原岩应力状态发生改变，形
成新的应力分布状态，这种变化称为应力重分布，变 化后的应力称为次生应力(二次应力场）。
（3）水平应力普遍大于铅垂应力
在绝大多数地区均有两个主应力位于水平或接近水平 的平面内，其与水平面的夹角一般不大于30° 最大水平主应力普遍大于垂直应力，它们的比值一般 在0.5～5.5之间，在很多情况下比值大于2。 垂直应力在大多数情况下为最小主应力，在少数情况 下为中间主应力，只在个别情况下为最大主应力。
原岩应力实测方法因其所用传感器原理和结构的
不同而有许多种：孔径变形法、水压致裂法、孔
壁应变法、空心包体应变法。 前二种是二维应力测量，可以测量与钻孔轴线垂 直截面上的原岩应力分量； 后二种是三维应力测量，可以一次测得全部六个 应力分量；
测量方法：
1）测量手段的不同划分： 构造法、变形法、电磁法、地震法、放射性法。 2）测量原理的不同划分： 应力恢复法、应力解除法、应变恢复法、应变解除法、水压 致裂法、声发射法、X射线法、重力法 3）测量基本原理不同划分：
地应力平衡汇总
1、“地应力平衡”的含义、目的、作用我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
2、地应力如何平衡地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
abaqus的part模块绘图功能不是很强，因此常用AutoCAD绘出平面图后导入abaqus。
在abaqus6.10中，主导入要分为以下几个步骤：1，在AutoCAD中建好模型的平面图形，并且另存为dxf格式。
地应力平衡方法介绍
6.19 Enhancements to the geostatic procedureProducts: Abaqus/Standard Abaqus/CAEBenefits: The geostatic procedure for obtaining the initial equilibrium state has been enhanced so that you no longer have to specify initial stresses that are close to the equilibrium state to obtain a solution corresponding to the original configuration.Description: The geostatic procedure is normally used as the first step of a geotechnical analysis; in such cases gravity loads (and possibly other types of loads) are applied during this step. Ideally, the loads and initial stresses should exactly equilibrate and produce zero deformations. However, in previous releases of Abaqus the geostatic procedure did not enforce this condition. In complex problems it may be difficult to specify initial stresses and loads that equilibrate exactly. Consequently, the displacements corresponding to the equilibrium solution might be large unless a special procedure is used to enforce small displacements.The enhanced geostatic procedure allows you to obtain equilibrium in cases when the initial stress state is unknown or is known only approximately. Abaqus automatically computes the equilibrium corresponding to the initial loads and the initial configuration, allowing only small displacements within user-specified tolerances. The procedure is available with continuum and cohesive elements with pore pressure degrees of freedom and the corresponding stress/displacement elements. The elastic, porous elastic, Cam-clay plasticity, and Mohr-Culomb plasticity material models are supported. Although the list of supported materials includes materials that exhibit inelastic behavior, the procedure is intended to be used in analyses in which the material response is primarily elastic; that is, inelastic deformations are small.The new enhancements are available from the Incrementation tabbed page when you create or edit a geostatic step in Abaqus/CAE. You must select automatic incrementation to access the new controls. The default settings for increment size and maximum displacement change are shown in Figure 6–6.Figure 6–6 The Incrementation options for a geostatic step.Abaqus/CAE Usage:Step module:Create Step: General: Geostatic; IncrementationReferences:Abaqus Analysis User's Manual∙“Geostatic stress state,” Section 6.8.2Abaqus/CAE User's Manual∙“Configuring a geostatic stress field procedure” in “Configuring general analysis procedures,” Section 14.11.1Abaqus Keywords Reference Manual∙*GEOSTATICAbaqus Verification Manual∙“*GEOSTATIC, UTOL,” Section 5.1.9ABAQUS 地应力平衡2011-03-30 14:42:12关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!1.先施加重力荷载的作用,可以在cae中实现;2.在inp文件中的output request中写上*el prints,这样就会将施加重力荷载后的应力输出到*.dat文件中了;3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如ELEMENT T FOOT- S11 S22 S33 S12 NOTE1 1 -1.2598E+05 -1434. -3.1852E+04 892.72 1 -1.2249E+05 -6287. -3.2194E+04 1223.3 1 -1.1795E+05 -497.7 -2.9611E+04 1664.4 1 -1.1210E+05 -7240. -2.9834E+04 1992.5 1 -1.0485E+05 579.0 -2.6068E+04 2600.6 1 -9.5803E+04 -8272. -2.6019E+04 3031.7 1 -8.4709E+04 1915.-2.0698E+04 4083.8 1 -7.0634E+04 -9746. -2.0095E+04 4339.9 1 -5.1088E+04 5401.-1.1422E+04 8519.10 1 -2.4353E+04 -1.1150E+04 -8876. 1.2126E+0411 1 -1.2847E+05 268.1 -3.2050E+04 738.112 1 -1.2786E+05 -9868. -3.4433E+04 629.113 1 -1.2938E+05 -4224. -3.3402E+04 502.514 1 -1.3039E+05 -3458. -3.3461E+04 165.9 单独存为一个*.dat文件,4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance. ,即单元编号变为: instance名称.序号 ;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!5.接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),6.最后在inp文件的step之前写上*initial conditions,type=stress,input=文件名.csv即可!这种方法不需要用python,比较简单,希望能对大家有用!先说为什么要施加地应力：1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
岩石力学基本教程 教学PPT 第6章 地应力综述.资料
2．高地应力现象
（3）探洞和地下隧洞的洞壁产生剥离，岩体锤击为嘶哑声并有较大变形，在 中等强度以下的岩体中开挖探洞或隧洞，高地应力状况不会像岩爆那样剧烈， 洞壁岩体产生剥离现象，有时裂缝一直延伸到岩体浅层内部，锤击时有嘶哑 声。在软质岩体中洞体则产生较大的变形，位移显著，持续时间长、洞径明 显缩小。
（1）地应力是一个具有相对稳定性的非稳定应力场，它是时间和空间的函数。
➢ 地应力在绝大部分地区是以水平应力为主的三向不等压应力场。三个主应 力的大小和方向是随着时间和空间而变化的，因而它是一个非均匀的应力 场。
➢ 地应力在空间上的变化，从小范围来看，它在空间上的变化是比较明显的， 但就某个地区整体而言，其变化并不大。
相邻岩体的约束，不可能产生横向变形，即 x y 0 。而相邻岩体的阻 挡就相当于对单元体施加了侧向应力 x 及 y ，考虑广义虎克定律则有：
xE 1[x (yz)]0
（6.3a）
由此可得
yE 1[y (zx)]0 x y1 z1 H
（6.3b） （6.4）
式中，E为岩体的弹性模量， 为岩体的泊松比。令 (1) ，则有：
斯蒂芬森（O. Stephansson）等人根据实测结果给出了芬诺斯堪的亚古陆最大水平 主应力和最小水平主应力随深度变化的线性方程：
最大水平主应力
h ， m a x 6 .7 0 .0 4 4 4 H （ M P a ）
最小水平主应力
h ， m i n 0 .8 0 .0 3 2 9 H （ M P a ）
a.共性 ➢ h=（1/4-1/5）φ； ➢ 钻进过程差异卸荷回弹，破裂主要发生在一定高度的岩芯根部； ➢ 拉张和剪切复合机制； b.产生条件： ✓ 弹性高，储能条件好的岩性条件，如火成岩； ✓ 整体块状; ✓ 高地应力条件， max ≥ 30MPa。
地应力平衡方法
地应力平衡方法：第一步：建立模型，材料，分析步（GEOSTATIC）第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8SimWe第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON) (请严格按照这个格式，注意大小写的字母)仿第四步：提交J0B，完成后第五步：按以下步骤,[Module]Visualization-Report---Report Field Output---下拉菜单里面选择centriod，然后依次把s11、s22、s33、s12、s13、s23点选上，setup页面选择报告文件的名字---***.txt---Write中选择Field Output-------------ok!!!第六步：先打开excel，在excel中打开***.txt，选择“分隔符号”—选择“导入的起始行”---点下一步-----在出来的界面上点“完成”------保存为CSV逗号分隔文件***.csv（文件内容格式如下）单元号S11 S22 S33 S12 S13 S23 （请注意，在保存内容中没有这一行的）1 , . , . , . , . , . ,2 , . , . , . , . , . ,. , . , . , . , . , . ,. , . , . , . , . , . ,这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些这些数字这些数字。
在CSV文件中，数字本身是用逗号隔开的，不需要另加逗号。
第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上：*initial conditions,type=stress,input=***.csv完成仿第八步：重新提交JOB，OK。
地应力测试方法
地应力主要测试方法总结摘要：本文总结了目前使用较为广泛的26种地应力测试，并对这些方法的基本原理做了简要介绍。
这26种方法按照数据源途径可以分为5大类，分别为基于岩芯的方法、基于钻孔的方法、地质学方法、地球物理学方法以及基于地下空间的方法。
最后文章对这些方法进行了的优缺点和适用范围进行了分析对比。
蓄存在岩体内部未受到扰动的应力称之为地应力，地应力可以分为两类，原地应力和诱发应力，而原地应力主要来自五个方面:岩体自重、地质构造活动、万有引力、封闭应力和外部荷载。
地应力具有多来源性且受到多种因素的影响，因此地壳岩体地应力分布复杂多变。
从海姆假说认为“岩体中赋存的应力近似为静水压力状态，且等于上覆岩体自重”到金尼克假说认为“垂直应力等于上覆岩体自重，水平应力等于岩体泊松效应产生的应力”，人们对岩体应力的认识逐步提高，并利用实测数据否定了以上两种假说。
社会发展的需求直接催生了大量地应力测试和估算方法，而这些方法的发展又进一步促进了人类社会的基础设施建设、资源和能源开发。
随着人类对能源和矿产资源需求量的增加和开采强度的不断加大，浅部矿产资源日益减少，国内外矿山都相继进入深部资源开发状态，而深部开采中遇到的“三高”问题(高地应力、高地温、高水压) 将成为深部开采岩体力学研究中的焦点和难点问题。
准确确定深部开发空间区域的原地应力状态是解决以上难题的必要途径之一，这就需要进行地应力测试方法和技术的研究。
从地应力概念提出至今，各国科学家提出了数十种地应力测试方法，将其按照数据来源进行归类，大概可以分为五大类:基于岩芯的方法、基于钻孔的方法、地质学方法、地球物理方法( 或地震学方法)、基于地下空间的方法。
下面将对各种方法的测试原理和方法发展的脉络作一些简要介绍，表1包括了目前认可程度和使用范围较广的各种方法.表1 原地应力测试和估算方法汇总1 基于岩心的方法1.1 非弹性应变恢复法非弹性应变恢复法(ASR)是通过测量现场从井孔取得的定向岩芯与时间相关的应变松弛变形来反演原地应力场方向和量值的一种方法。
地应力平衡
ABAQUS地应力平衡：进行地应力平衡的原因陈述如下:我们建立的几何模型一般都和工程实际情况一致,例如边坡的几何模型与边坡实际尺寸相一致。
但是由于边坡的沉降和徐变作用，可以想像到，现在的边坡应该是由一个体积更大的原始边坡在很久以前由于受到重力作用和边界约束条件,逐渐形成了现今的边坡形态.但是对于那个原始的边坡形态，我们不得而知。
假如能准确知晓，我们就能够建立原始边坡的几何模型，接着对边坡施加重力和边界条件,受力后边坡形态应该和现在的边坡相一致，其内力就是初始应力场(地应力），这样就不用专门施加地应力了。
但现实情况是我们不能知晓原始边坡的形态.现在的边坡几何模型就是其实际形态，受力之后将会变成一个与现状不一致的边坡，这不符合现在的实际情况.如果我们计算出现今边坡的内力，并将其作为边坡的初始应力场，再去和外力平衡，这样我们建立的模型就和现实边坡情况相一致了。
对于涉及开挖、回填的动态岩土工程问题，地应力平衡是正确模拟施工过程的前提条件。
初始应力的加载必须满足地应力平衡,而地应力平衡就是为了使地基仅存在初始应力，而不存在初始应变。
当地基自重是产生地应力场的主要因素时,重力是外力，初始应力场是内力，将提取出的内力施加于模型后再施加重力,此时内力和外力平衡，该状态就是工程建设的初始状态。
在表面水平的情况下，ABAQUS中初始地应力场的平衡一般只和密度有关,土体的密度一样，平衡的效果就好，别的参数对地应力平衡的结果影响很小。
对于表面不平的情况,尽量通过inp文件导入初始应力的方法进行地应力平衡。
ABAQUS中进行地应力平衡的时间点的选择十分重要,地应力平衡是指在工程建设之前，地表的位移应为零, 而土体的应力却存在。
也就是说不管土体原来的样子如何（例如高山, 河流，丘陵，平原等), 进行地应力平衡的正确时间点应当是在我们对它做任何扰动之前。
具体采取的办法如下所述，我们对所建立的边坡几何模型施加和实际模型一致的重力和边界条件，得到变形以后边坡的内力，变形后边坡形状和原始边坡略有不同，其内力可近似作为现状边坡的内力,将其作为初始应力施加于现在的边坡中，接着施加外力(重力）来平衡初始应力，这样就建立了一个与现今边坡形态基本相同的边坡模型，这样之后的分析计算才是符合实际的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第六篇：地应力平衡方法以及注意事项
注意：只有采用弹塑性本构模型时需要地应力平衡，弹性本构不需要地应力平衡！
第一部分地应力平衡方法简介
地应力平衡主要有五种方法：
（1）自动平衡：第一步创建分析步geostatic ，这种方法注意只能在第一步只能有土和重力的情况下能使用，有其他部件或者接触时计算不能收敛，效果是最好的，方便简单！
（2）*initial conditions,type=stress,geostatic
该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。
所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。
计算速度快,收敛性好。
缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。
高版本在CAE里也能操作
用计算器算出每个分界面上的应力和坐标对应填入，也比较方便不需要修改关键字
（3）*initial conditions,type=stress,geostatic,user
导入ODB里的方法，也比较简单，高版本可在截面上操作，不需要修改关键字
你放入ODB后，填入第一步不需要填名称就是1 ，增量步就是你第一步计算的最后一个增量步
（4）*initial conditions,type=stress,input=FileName.csv(或inp)
该方法中的文件FILENAME.INP获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。
但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。
在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应,因此平衡后的效果不一定很理想,但无疑其适用性很强。
最麻烦，也是最容易出错的，不喜欢使用
（5）*initial conditions,type=stress,geostatic,user
该方法采用用户子程序SIGINI来定义初始应力场,可以定义其为应力分量为坐标、单元号、积分点号等变量的函数,要达到精确平衡需已知具体边界条件,在实际中应用较少。
注意：除了第一种方法geostatic 选的是自动，其他都是fixed,其他四种方法第一步除了土外可以有其他部件，效果一般很差，可以多循环几次！。