ABAQUS岩土隧道入门地应力平衡基本问题
abaqus 位移边界 地应力平衡
Abaqus 是一种广泛应用于工程领域的有限元分析软件，它可以用来模拟各种结构和材料在不同工况下的受力和变形情况。
其中，位移边界和地应力平衡是在地下工程和土木工程中常常需要考虑的重要问题。
一、Abaqus 软件简介Abaqus 是由达索系统公司开发的有限元分析软件，它可以模拟和分析各种结构和材料在受力和变形下的行为。
Abaqus 具有较为完善的功能和灵活的应用，广泛应用于航空航天、汽车制造、电子设备、地下工程和土木工程等领域。
二、位移边界的作用在有限元分析中，位移边界是指在结构或材料的某些边界上规定位移的边界条件。
通过位移边界的设置，可以模拟结构或材料在受力情况下的变形行为，并分析其受力性能和稳定性。
地下工程和土木工程中的许多问题都需要考虑位移边界的影响，如地基沉降、支撑结构的变形等。
1. 位移边界的设置方法在 Abaqus 软件中，可以通过定义边界条件的方式来设置位移边界。
用户可以选择固定边界、自由边界、受控边界等不同的边界条件，以模拟结构或材料在不同工况下的位移情况。
通过合理设置位移边界条件，可以准确模拟受力下的变形行为，并分析结构的稳定性和安全性。
2. 位移边界的影响位移边界的设置对结构或材料的受力和变形行为具有重要影响。
合理设置位移边界可以准确模拟受力下的变形行为，帮助工程师分析结构的稳定性和可靠性，指导工程设计和施工实践。
三、地应力平衡的问题地下工程和土木工程中，地应力平衡是一个重要的问题。
地下结构如地下室、隧道、地铁等的建设，往往需要考虑地应力平衡对结构的影响，以确保结构的稳定性和安全性。
1. 地应力平衡的影响地应力平衡对地下结构的影响主要体现在地应力的变化和传递上。
地下结构的施工和使用过程中，常常会引起地应力的变化，导致结构的变形和破坏。
合理分析和考虑地应力平衡对结构的影响，可以有效指导工程设计和施工实践，确保地下结构的稳定和安全。
2. Abaqus 软件在地应力平衡分析中的应用Abaqus 软件具有强大的地下工程分析功能，可以模拟地下结构在地应力作用下的受力和变形情况。
abaqus初始地应力法
Abaqus初始地应力法1. 什么是初始地应力法？初始地应力法是一种通过在有限元分析中引入初始地应力来模拟复杂的地下工程问题的方法。
它主要用于研究岩土体在开挖、填筑、注浆等施工过程中的变形和破坏行为，以及对结构物产生的影响。
2. 初始地应力的来源初始地应力主要来自以下几个方面：•地壳运动：包括板块运动、构造活动等，会引起地层的水平和垂直变形，进而产生初始地应力。
•地震活动：地震引起的强烈振动会导致岩土体发生变形，从而产生初始地应力。
•自重压力：岩土体受到自身重量作用时，会产生垂直于重力方向的初始地应力。
•水压作用：水对岩土体施加的压力也会导致初始地应力。
3. 初始地应力法在Abaqus中的实现步骤步骤1：建立模型首先，在Abaqus中建立模型。
可以选择二维平面或三维立体模型，根据实际情况选择适当的模型尺寸和几何形状。
步骤2：定义材料属性在Abaqus中定义岩土体的材料属性，包括弹性模量、泊松比、密度等。
根据实际情况选择合适的材料模型，如线弹性、Mohr-Coulomb等。
步骤3：划分单元网格将模型划分为有限元网格。
可以选择不同类型的单元，如三角形、四边形、六面体等。
根据问题的复杂程度和准确性要求，确定合适的网格密度。
步骤4：施加边界条件在Abaqus中施加边界条件，包括约束和加载。
约束条件用于限制结构或岩土体的运动自由度，加载条件用于施加外部荷载或位移。
步骤5：设置初始地应力在Abaqus中设置初始地应力。
可以通过以下方法进行设置：•直接指定初始地应力值；•根据已知地层压缩模量和厚度计算得到；•根据现场实测数据进行反演得到。
步骤6：求解并分析结果在Abaqus中进行数值求解，并对结果进行分析。
可以观察岩土体的变形、应力分布等情况，评估结构物的稳定性和安全性。
4. 初始地应力法的应用领域初始地应力法广泛应用于以下领域：•地下工程：如隧道、地铁、基坑开挖等。
•岩土工程：如边坡稳定性分析、地基承载力计算等。
abaqus地应力设置方法
abaqus地应力设置方法宝子，今天咱来唠唠Abaqus里地应力的设置方法哈。
在Abaqus里设置地应力呢，你得先明白啥是地应力。
简单说呢，就是地层里本来就存在的应力。
这就像咱们住在房子里，房子本身的结构是承受着一些压力的，地层也一样。
那咋在Abaqus里搞这个设置呢？一种常见的办法是通过初始应力场来设置。
你得在模型里定义这个初始的应力状态。
这就好比你给你的模型打个底，告诉它最开始的时候应力是啥样的。
比如说，你可以根据实际的地质情况，像地层的深度啊，岩石的类型啊这些，去估算出初始的应力值。
然后把这些值输入到Abaqus里相应的模块中。
还有哦，如果你有一些实测的数据，那就更好啦。
你可以直接把这些实测的地应力数据导入到Abaqus里。
这就像是给你的模型吃了一颗定心丸，让它按照真实的情况来模拟。
不过呢，导入数据的时候可不能马虎，得按照Abaqus要求的格式来哦。
另外呀，在设置地应力的时候，你得考虑模型的边界条件。
这就像是给你的模型设定一个活动范围一样。
如果边界条件没设置对，那地应力的模拟可能就会出岔子。
比如说，在一个简单的平面应变模型里，边界上的应力是怎么传递的，你得心里有数。
而且哦，不同的材料模型对于地应力的响应也不一样呢。
就像不同性格的人对同一件事的反应不同。
有些材料可能比较“强硬”，能承受较大的地应力，有些就比较“脆弱”。
所以你在选择材料模型的时候，也要把地应力的因素考虑进去。
宝子，设置Abaqus的地应力虽然有点小复杂，但只要你把这些关键的点都搞清楚了，就不会出大错啦。
多试试，多做做模拟，你就会越来越熟练的。
加油哦，希望你能在Abaqus里顺利设置好地应力，做出超棒的模拟结果！。
ABAQUS常见问题总结
ABAQUS注意事项1、建模前需要考虑的因素1)根据结构和荷载情况的特点，按照轴对称问题来建模。
2)对于大位移问题，应在step功能模块中把参数Nlgeom（几何非线性）设为ON。
3)根据问题的特点，选用合适的单元类型。
2、新建Part时，Approximate size (sketch) 的大小?Approximate size数值的大小，应根据模型的最大尺寸来确定：稍大于最大尺寸的2倍。
3、ABAQUS中平面应力、应变问题的截面属性为什么是实体而不是壳？那个壳设置是专门针对板壳单元的，用于板壳力学分析。
平面应力应变分析当然该选实体，因为是我们把三维实体分析简化成平面分析的。
3、ABAQUS有限元分析实例详解——石亦平注意事项1)P176页在单向压缩试验过程模拟时，试样冒设为解析刚体，建模时用一直线表示，而不是用一矩形表示，如果用矩形表示，其接触部位的尖角会造成错误的分析结果。
同时，代表试样冒的直线也必须绘制的足够长，因为压缩的过程中，试样的截面尺寸会不断增大，如果试样尺寸超出试样冒的尺寸，会造成接触分析的收敛问题。
4、选择主动面和从属面的几个原则1)Analytical rigid surfaces and rigid-element-based surfaces must always be the mastersurface.2)A node-based surface can act only as a slave surface and always uses node-to-surfacecontact.3)Slave surfaces must always be attached to deformable bodies or deformable bodiesdefined as rigid.4)Both surfaces in a contact pair cannot be rigid surfaces with the exception of deformablesurfaces defined as rigid5)当存在一个较小的面和一个较大的面时，一般将较小的面定义为从属面。
ABAQUS的平衡初始地应力的方法研究
２）岩土材料的特性与应力状态密切相关； ３）对于涉及开挖、填充的动态岩土工程问题，初始应力 场是正确模拟其施工过程的先决条件［５］。在岩土工程分析 当中，初始地应力场平衡一般都放在第一步，即对土体施加 重力。理想的 情 况 是 该 作 用 力 与 土 体 的 初 始 应 力 正 好 平 衡，使得土体的初始位移为 ０，使其对后续的分析不产生影
模型尺寸：１０ｍ×１０ｍ×１０ｍ，图 １为计算模型：３Ｄ实 体单元类 型。采 用 理 想 弹 塑 性 本 构 模 型，服 从 Ｍｏｈｒｃｏｕ ｌｏｍｂ屈服准则。
Ｘ方向约束 Ｕ１（Ｕ１ ＝０），Ｙ方 向 约 束 Ｕ２（Ｕ２ ＝０），Ｚ 方向约束 Ｕ３（Ｕ３ ＝０），重 力 加 速 度 作 用 在 Ｚ方 向 上，大 小为 １０。
导 入 法 进 行 结 果 对 比 研 究 。 首 先 ，通 过 一 个 简 易 模 型 对 两 种 平 衡 方 法 进 行 了 简 单 的 描 述 ，并 将 两 个 结 果 与 未 进 行 平 衡 的 结
果比较，说明了地应力平衡的重要性。然后根据一个工程实例通过上述两种方法进行计算，结果 发 现：在 土 体 模 型 简 单 的
第２０４２６卷０年第１２１１期 月
山 西 建 筑
ＳＨＡＮＸＩ ＡＲＣＨＩＴＥＣＴＵＲＥ
ＮＶｏｏｖｌ．．４ ６２Ｎ０ｏ２．０２１
·１·
·专家荐稿·
文章编号：１００９６８２５（２０２０）２１０００１０３
ＡＢＡＱＵＳ的 平 衡 初 始 地 应 力 的 方 法 研 究
Abaqus地应力平衡方法
：地应力平衡方法：第一步：建立模型，材料，分析步（GEOSTATIC）第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)(请严格按照这个格式，注意大小写的字母第四步：提交J0B，完成后第五步：按以下步骤,Roport---ReportField Output---选中S11,S22,S33,S12,S13,S23---Name:XX.INP---Write中选择Field Output-------------ok!!!第六步：用软件(excel就行)打开XX.INP，保存格式内容单元号S11S22S33S12S13S23（请注意，在保存内容中没有这一行的）,.,.,.,.,.,,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些数字。
第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上：*initial conditions,type=stress,input=xx.csv完成,第八步：重新提交JOB，OK在ABAQUS中对应变的部分理解1、E—总应变；Eij—应变分量2、EP---主应变；EPn----分为Minimum, intermediate, and maximum principal strains (EP1《EP2< EP3)3、NE----名义应变；NEP---主名义应变；4、LE----真应变（或对数应变）；LEij---真应变分量；LEP---主真应变；5、EE—弹性应变；6、IE---非弹性应变分量；7、PE---塑性应变分量；8、PEEQ---等效塑性应变---在塑性分析中若该值〉0，表示材料已经屈服；描述整个变形过程中塑性应变的累积结果；若单调加载则PEEQ=PEMAG ；9、PEMAG----塑性应变量（幅值Manitude）---描述变形过程中某一时刻的塑性应变，与加载历史无关；10、THE---热应变分量；。
基于ABAQUS岩土工程中地应力平衡的探讨
基于ABAQUS岩土工程中地应力平衡的探讨作者：马云峰来源：《科技与创新》2014年第08期摘要：平衡地应力是岩土工程计算分析中十分重要的部分，在岩土施工前保证地表位移近似为零，是数值模拟的首要前提。
利用ABAQUS对某基坑开挖前的土体进行地应力平衡，根据现场土体的物理性质指标，将地应力平衡前后土体的应力和位移进行对比。
进行地应力平衡后土体的位移数量级小于10-6，并且任意深度的竖向应力值等于该处深度与土体密度的乘积。
关键词：岩土工程；ABAQUS；地应力；土体位移中图分类号：TU431文献标识码：A 文章编号：2095－6835（2014）08－0061－02随着科学技术水平的提高，越来越多的有限元软件被应用于建筑、材料、机械、航天等领域。
这些有限元软件的出现对岩土工程来说是一个巨大的推动，其中，作为最先进的大型通用有限元分析软件之一的ABAQUS被广泛应用。
ABAQUS软件全面满足了静力许可，可以模拟多种不同模型的本构关系，结合有限元分析方法，可以使数值模拟不受几何形状的不规则、边界条件的多样性和材料部均匀性的限制，能够更加准确地分析应力和应变的过程，为现场实际操作提供有效的依据。
1本构模型土体是天然的、复杂的多孔材料，在受荷之后会表现出明显的压缩性、非线性、流变性、各向异性、剪切性等特征。
为了更好地描述土体的真实力学——变形特性，建立其应力应变和时间的关系，在各种试验和工程实践经验的基础上提出一种数学模型，即土体的本构模型。
世界各国学者提出过上百种不同的土体本构模型，其中，摩尔库仑模型、扩展D-P模型和修正剑桥模型的应用最广。
1.1摩尔库仑模型摩尔库仑的屈服准则假定作用在某一点的剪应力等于该点的抗剪强度时，该点发生破坏，剪切强度与作用在该面的正应力呈线性关系。
1.2扩展D-P模型扩展D-P模型允许材料各向同性硬化或软化，考虑了材料的剪胀性。
扩展D-P模型的屈服准则主要取决于屈服面在子午面的形状。
abaqus地应力平衡的步骤
abaqus地应力平衡的步骤一、前期准备在进行地应力平衡分析之前，需要进行一系列的准备工作。
这包括收集地质勘察资料、了解场地工程地质条件、确定分析目标和评估范围等。
此外，还需要进行必要的模型简化与假设，以便在有限元分析中准确模拟实际地质结构和应力分布。
二、建立模型在abaqus中建立地应力平衡模型通常涉及以下步骤：创建模型空间、设置网格密度和类型、导入地质勘察数据（如岩石层分布、岩土性质等）以及为不同的材料设置相应的物理性质（如弹性模量、泊松比等）。
此外，还需根据实际工程需求和地质条件，设置边界条件和载荷，以模拟地质结构和应力分布。
三、设置材料属性材料属性的设置对于地应力平衡分析至关重要。
需要根据地质勘察资料和试验数据，为不同的岩土层设置准确的物理性质，包括密度、弹性模量、泊松比以及摩擦角和内聚力等。
这些参数将直接影响模型的计算结果和应力分布状态。
四、地应力平衡分析在进行地应力平衡分析时，需要采用合适的求解器和算法，以确保计算效率和准确性。
abaqus提供了多种求解器选择，如静力分析、动力分析以及流体分析等。
在地质工程领域，静力分析是最常用的求解器，用于模拟岩土体的稳态应力分布状态。
在设置求解器之后，需要设定合适的迭代步长和收敛准则，以便模型在迭代过程中逐渐达到平衡状态。
五、后处理与结果分析地应力平衡分析完成后，需要利用abaqus的后处理功能对结果进行详细分析。
这包括绘制应力分布云图、查看节理和断层的应力状态、评估潜在的破坏区域等。
通过对结果的深入分析，可以评估地质结构的稳定性和安全性，为进一步的地质工程设计和优化提供依据。
六、优化与迭代根据后处理阶段的结果分析，可能需要对模型进行优化和迭代。
这包括调整网格密度和类型、改进边界条件和载荷设置、调整材料参数等。
通过不断优化和迭代，可以提高模型的精度和可靠性，从而更准确地模拟实际地质结构和应力分布状态。
七、注意事项在地应力平衡分析过程中，需要注意以下几点：1.确保收集足够的地质勘察资料和试验数据，以便为模型提供准确的材料属性和边界条件。
abaqus动力边界地应力平衡
abaqus动力边界地应力平衡1. 介绍在工程结构分析中，地应力是一个重要的影响因素，特别是对于地下工程和岩土工程来说。
abaqus是一款常用的有限元分析软件，可以用于分析各种工程结构的热、力学和动力学问题。
在abaqus中，通过设置动力边界条件来模拟地应力的作用，保证结构在地应力的作用下可以达到平衡状态。
2. 地应力的作用地应力是指在地下岩土中由地球自身重力和地球内部物质的重力所产生的应力。
它对地下工程结构的稳定性、变形和破坏具有重要的影响。
在地下工程中，结构体受到地应力的约束，必须在地应力的作用下达到平衡状态，因此需要进行地应力平衡的分析。
3. abaqus中动力边界条件的设置在abaqus中，可以通过设置动力边界条件来模拟地应力的作用。
首先需要建立模型，并对模型进行网格划分。
然后可以通过创建荷载来模拟地应力的作用，通常使用压力荷载来表示地应力的作用。
在abaqus中，可以选择不同的压力类型，比如静态压力、动态压力或者地震荷载等，根据实际情况选择合适的压力类型。
在设置动力边界条件时，需要考虑地应力的方向和大小，以及结构的几何形状和材料特性等因素。
4. 地应力平衡的分析方法在abaqus中，可以通过施加地应力荷载和设置动力边界条件来进行地应力平衡的分析。
首先需要对结构进行静力分析，确定结构的初始受力状态。
然后根据地应力的方向和大小，在结构表面施加相应的地应力荷载。
接着设置动力边界条件，对结构进行动力分析，通过迭代计算使结构达到地应力平衡状态。
最终得到结构在地应力作用下的平衡位移和应力分布情况。
5. 地应力平衡分析的应用地应力平衡分析在地下工程和岩土工程中具有广泛的应用。
比如在隧道、地铁、水坝、地下管道等工程中，地应力的影响必须得到合理的分析和处理。
通过使用abaqus软件进行地应力平衡分析，可以更准确地预测结构的受力状态，提高工程设计的安全性和可靠性。
6. 结论abaqus是一款强大的有限元分析软件，可以用于地应力平衡分析。
abaqus地应力平衡
初始地应力平衡应该满足：（1）竖向位移为零或者很小，一般至少10的负6次方，甚至更小（2）竖向应力s22要等于土体深度*密度*g关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!1.先施加重力荷载的作用,可以在cae中实现;2.在inp文件中的output request中写上*el prints,这样就会将施加重力荷载后的应力输出到*.dat文件中了;3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如ELEMENT T FOOT- S11 S22 S33 S12 NOTE1 1 +05 -1434. +042 1 +05 -6287. +04 1223.3 1 +05 +04 1664.4 1 +05 -7240. +04 1992.5 1 +05 +04 2600.6 1 +04 -8272. +04 3031.7 1 +04 1915. +04 4083.8 1 +04 -9746. +04 4339.9 1 +04 5401. +04 8519.10 1 +04 +04 -8876. +0411 1 +05 +0412 1 +05 -9868. +0413 1 +05 -4224. +0414 1 +05 -3458. +04单独存为一个*.dat文件,4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance. ,即单元编号变为: instance名称.序号 ;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!5.接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),6.最后在inp文件的step之前写上*initial conditions,type=stress,input=文件名.csv即可!这种方法不需要用python,比较简单,希望能对大家有用!先说为什么要施加地应力：1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
abaqus初始地应力法
abaqus初始地应力法【原创实用版】目录1.Abaqus 初始地应力法概述2.Abaqus 初始地应力法的应用3.Abaqus 初始地应力法的步骤4.Abaqus 初始地应力法的优缺点5.总结正文一、Abaqus 初始地应力法概述Abaqus 初始地应力法是一种在 Abaqus 软件中应用的求解地应力场的方法，主要通过指定主应力方向、应力值以及边界条件，来计算地应力场中的应力分布。
此方法在岩土工程、地质工程等领域具有广泛的应用。
二、Abaqus 初始地应力法的应用Abaqus 初始地应力法可以用于求解许多工程问题，例如：地下洞室的开挖、隧道的开挖与支护、边坡稳定性分析等。
通过该方法，工程师可以预测工程施工过程中可能产生的应力变化，从而采取相应的措施确保工程安全。
三、Abaqus 初始地应力法的步骤1.准备模型：首先需要创建一个三维模型，包括地基、岩土体等部分。
2.指定边界条件：根据实际情况，设置模型的边界条件，如自由表面、滑动面、固定面等。
3.指定初始地应力：在模型中指定主应力方向、应力值，作为计算的起始条件。
4.分网格：将模型划分为有限元网格，以进行数值计算。
5.求解：运用 Abaqus 软件求解地应力场，得到应力分布。
6.后处理：对计算结果进行可视化和分析，提取感兴趣的应力参数。
四、Abaqus 初始地应力法的优缺点优点：1.适用于复杂的地质条件和工程环境。
2.可以考虑多种边界条件，模拟实际情况。
3.能够得到应力分布的三维图像，便于分析和观察。
缺点：1.计算过程较为繁琐，需要一定的技术水平。
2.对模型的精度和网格划分要求较高，否则可能影响计算结果的准确性。
五、总结Abaqus 初始地应力法是一种实用的地应力场分析方法，适用于多种工程场景。
abaqus水平地应力平衡初始应力提取法
abaqus水平地应力平衡初始应力提取法文章标题：深入剖析abaqus水平地应力平衡初始应力提取法目录一、引言二、abaqus水平地应力平衡的意义三、abaqus水平地应力平衡初始应力提取法详解3.1 地应力平衡的基本原理3.2 初始应力提取法的流程和步骤3.3 初始应力提取结果的分析和应用四、我的理解与观点五、结语一、引言在地下工程设计与施工中，了解地下地应力分布对于合理设计和施工极为重要。
地应力平衡是地下工程中一个非常重要的概念，而abaqus 水平地应力平衡初始应力提取法则可以帮助我们更好地了解地下地应力情况。
本文将就abaqus水平地应力平衡初始应力提取法进行深入探讨，帮助读者更好地理解和应用这一方法。
二、abaqus水平地应力平衡的意义在地下工程中，地应力平衡是指地下岩土或者岩石在无外荷载或外荷载不变的情况下，达到了内部平衡状态。
了解地应力平衡对地下工程的设计、分析和施工具有重要意义。
abaqus水平地应力平衡初始应力提取法通过模拟地下地应力平衡状态，提取并分析地下初始应力分布情况，为地下工程设计和施工提供了重要参考依据。
三、abaqus水平地应力平衡初始应力提取法详解3.1 地应力平衡的基本原理地应力平衡是指地下岩土或者岩石在无外荷载或外荷载不变的情况下，达到了内部平衡状态。
地下岩土或者岩石的地应力状态与其应力历史、地层岩性以及地下水情况等密切相关。
了解地下地应力平衡状态对地下工程设计与施工具有重要意义。
3.2 初始应力提取法的流程和步骤1）模拟地下地应力平衡状态2）提取地下初始应力数据3）分析地下初始应力分布4）应用地下初始应力数据3.3 初始应力提取结果的分析和应用通过abaqus水平地应力平衡初始应力提取法，可以获取地下初始应力分布情况，为地下工程设计与施工提供重要参考。
地下初始应力数据不仅可以用于地下工程的稳定性分析，还可以用于地下工程的设计参数确定。
四、我的理解与观点作为地下工程设计与施工领域的从业者，我深切理解abaqus水平地应力平衡初始应力提取法对地下工程的重要意义。
基于ABAQUS岩土工程中地应力平衡的探讨
基于ABAQUS岩土工程中地应力平衡的探讨作者：马云峰来源：《科技与创新》2014年第08期摘要：平衡地应力是岩土工程计算分析中十分重要的部分，在岩土施工前保证地表位移近似为零，是数值模拟的首要前提。
利用ABAQUS对某基坑开挖前的土体进行地应力平衡，根据现场土体的物理性质指标，将地应力平衡前后土体的应力和位移进行对比。
进行地应力平衡后土体的位移数量级小于10-6，并且任意深度的竖向应力值等于该处深度与土体密度的乘积。
关键词：岩土工程；ABAQUS；地应力；土体位移中图分类号：TU431文献标识码：A 文章编号：2095－6835（2014）08－0061－02随着科学技术水平的提高，越来越多的有限元软件被应用于建筑、材料、机械、航天等领域。
这些有限元软件的出现对岩土工程来说是一个巨大的推动，其中，作为最先进的大型通用有限元分析软件之一的ABAQUS被广泛应用。
ABAQUS软件全面满足了静力许可，可以模拟多种不同模型的本构关系，结合有限元分析方法，可以使数值模拟不受几何形状的不规则、边界条件的多样性和材料部均匀性的限制，能够更加准确地分析应力和应变的过程，为现场实际操作提供有效的依据。
1本构模型土体是天然的、复杂的多孔材料，在受荷之后会表现出明显的压缩性、非线性、流变性、各向异性、剪切性等特征。
为了更好地描述土体的真实力学——变形特性，建立其应力应变和时间的关系，在各种试验和工程实践经验的基础上提出一种数学模型，即土体的本构模型。
世界各国学者提出过上百种不同的土体本构模型，其中，摩尔库仑模型、扩展D-P模型和修正剑桥模型的应用最广。
1.1摩尔库仑模型摩尔库仑的屈服准则假定作用在某一点的剪应力等于该点的抗剪强度时，该点发生破坏，剪切强度与作用在该面的正应力呈线性关系。
1.2扩展D-P模型扩展D-P模型允许材料各向同性硬化或软化，考虑了材料的剪胀性。
扩展D-P模型的屈服准则主要取决于屈服面在子午面的形状。
ABAQUS常见问题汇总 - 1.0版
ABAQUS常见问题汇总 - 1.0版作者：aba_aba2006.06.14在SIMWE论坛注册5个月以来，在ABAQUS版回帖860篇，知无不言，言无不尽。
在此整理一下，以方便大家查找。
下面绝大多数是我的回帖，有少部分是其他网友的回答，都注明了作者。
下面只摘录了帖子中的一些主要内容，有些地方可能上下文不太连贯，完整的讨论请大家根据相应链接去论坛上察看。
如果这些内容对你有帮助，希望你也能常上论坛来回答一下别人的问题。
众人拾柴，火焰才高。
目录点击小节标题，可以跳到相应的内容（有些WORD版本可能需要按住ctrl键）1. 论坛提问要诀 (4)2. 故障诊断基本方法 (4)3. 约束刚体位移 / NUMERICAL SINGULARITY (4)4．过约束（Overconstraint） (7)5．其他常见error和warning信息 (7)6．ABAQUS无法运行（安装、License、环境变量） (10)7．弹塑性分析 (14)7.1 塑性问题不收敛的常见现象 (14)7.2 接触问题和塑性材料不要用二阶单元 (14)7.3 不要在塑性材料上施加点载荷 (15)7.4 塑性材料参数 (15)7.5 屈服准则 (18)7.6 后处理 / 判断材料是否屈服 (20)8. 接触分析 (22)8.1 接触分析不收敛的常见现象和解决方法 (22)8.2 接触面上的网格密度 (22)8.3 接触面的法线方向 (23)8.4 过盈接触 (23)8.5 管土/桩土接触 (25)8.7 凹坑成型的接触问题 (34)8.8 刚体穿透 (37)8.9 接触力 (37)8.10 Explicit中的接触问题 (38)8.11 改变接触状态 (40)8.12 其他与接触有关的Error和Warning信息 (41)9. 各个量的单位 (42)9.1 各个量的单位要相互一致 (42)9.2 时间/增量步的含义 (44)10．导入 / 修改 / 求解inp文件 (45)10.1 把inp文件导入ABAQUS/CAE (45)10.2 运行inp文件 (46)10.3 修改inp文件 (47)11．用户子程序和FORTRAN (49)12．创建/修改part (50)13．在Assembly中定位 (51)14．划分网格 (52)15．选择单元类型 (54)16．截面属性和材料 (56)16.1 定义截面属性 (56)16.2 用OFFSET参数偏置shell的中性面 (57)16.3 材料属性 (57)17．载荷 / 速度场 (58)17.1 集中力 (58)17.2 线载荷 (59)17.3 面载荷 (60)17.4 Body Force / Gravity (60)17.5 随时间变化的载荷 / Amplitude (61)17.6 在实体单元上施加弯矩 / 旋转速度 (62)18．边界条件 (67)19．定义Set / 选取模型的局部 (67)20．弹簧 (68)21．约束（coupling, tie） (68)22．后处理 (69)22.1 显示局部坐标系上的结果 (69)22.2 绘制曲线（X–Y data） (69)22.3 field output / history output (70)22.4 显示结果 (71)22.5 输出shell积分点上的应力和应变 (73)23．多步骤分析 (74)23.1 重启动分析（Restart） (74)23.3 用*import 把计算结果传递到下一个过程中 (74)24．岩土分析 (75)24.1 初始地应力平衡 (75)24.2 基坑开挖 (77)24.3 钢筋混凝土（Rebar Layer和Embedded Element） (77)24.4 实体单元的截面力/弯矩/转角 (78)24.5 软土固结 (80)25．热分析 (81)26．动态分析 (83)26.1 Explicit分析 (83)26.2 分析结果与amplitude不一致 (84)26.3 固有频率 (84)26.4 *DYNAMIC分析 (85)27．屈曲分析 (86)28．裂纹 (87)29．连接单元（CONNECTOR） (88)30．惯性释放（INERTIA RELIEF） (88)31．非线性问题的算法 (89)32．ABAQUS和其它软件的比较 (89)33．使用ABAQUS帮助手册 (90)34．论坛的“全文搜索”功能 (90)35．关于找枪手 (91)1. 论坛提问要诀在论坛上的提问是否能够得到解答，在很大程度上取决于提问的方式。
abaqus地应力平衡
我也是新手这些资料你看看吧！都是在网上找的，看完应该能明白个大概！初始地应力平衡应该满足：（1）竖向位移为零或者很小，一般至少10的负6次方，甚至更小（2）竖向应力s22要等于土体深度*密度*g关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了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年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
初始地应力平衡方法
初始地应力平衡方法当模型并不是水平地面或地基成层分布时，用Initial conditions, type = stress, geostatic方法难以奏效，这里介绍一种可适用于复杂地形下的初始地应力平衡法。
1 CAE中建立仅有重力作用的模型，在提交Job前，输入Python命令产生没有Part信息的Input文件mdb.models['模型名字'].setValues(noPartsInputFile=ON)2 提交Job，求解完成后，进入CAE后处理部分1) Report——Report Field Output——下拉菜单中选择Centriod——勾选上S11、S22、S33、S12、S13和S23；2) Setup页面选择拟生成报告文件的名字***.CSV，其下方Write栏里仅选择Field Output。
3 用软件UltraEdit打开***.CSV文件进行编辑1) 删除CSV文件中不必要的内容，仅保留以下内容；单元号S11 S22 S33 S12 S13 S231 -29.3229E+03 -54.4569E+03 -29.3229E+03 23.1405E-12 339.116E-15 -53.1765E-122 -29.3229E+03 -54.4569E+03 -29.3229E+03 -14.2653E-12 -2.41752E-12 -39.8307E-12……n -29.3229E+03 -54.4569E+03 -29.3229E+03 -42.7399E-12 -4.38639E-12 -3.0521E-122) 利用UltraEdit的“列模式”，在各列中进行“插入/填充列”操作，添加逗号，并删除文件第一行和末尾的空行，保存文件。
1 , -29.3229E+03 , -54.4569E+03 , -29.3229E+03 , 23.1405E-12 , 339.116E-15 , -53.1765E-122 , -29.3229E+03 , -54.4569E+03 , -29.3229E+03 , -14.2653E-12 , -2.41752E-12 , -39.8307E-12……n , -29.3229E+03 , -54.4569E+03 , -29.3229E+03 , -42.7399E-12 , -4.38639E-12 , -3.0521E-12Tips:1) CSV文件中所保存的内容只有数字；2) 这里选用UltraEdit而不是Excel处理CSV文件是因为：Excel有行数限制，当n>65536时，将不能载入全部文件内容。
第六篇：地应力平衡方法以及注意事项
第六篇：地应力平衡方法以及注意事项注意：只有采用弹塑性本构模型时需要地应力平衡，弹性本构不需要地应力平衡！第一部分地应力平衡方法简介地应力平衡主要有五种方法：（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个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。
但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。
基于ABAQUS的无渗流偏压隧道开挖分析(包含地应力平衡分析)
基于ABAQUS的无渗流偏压隧道开挖分析（包含地应力平衡分析）1 隧道建模及地应力平衡1.1 工程概况本模型截取的一段隧道通过山体坡度30°左右一侧，隧道开挖及初衬断面为五心圆各项参数如表2.1所示。
隧道跨度13.36m，高11.71m，偏压部分覆盖层厚度约为20m。
隧道区地下水主要为第四系松散层孔隙水和基岩裂隙水，垂直隧道中轴线向东120米勘察孔水位为65米，向西120米勘察孔水位为10米。
施工过程中隧道单位正常涌水量3.43m3/d，中等富水。
具体围岩初衬参数将在后文建模时给出。
表2.1 隧道断面参数(单位：cm)r[1]R[1]r[4]R[4]h[1]h[2]h[3]a b h H B1.2 隧道建模及偏压判定结合隧道工程勘察报告和规范取围岩和初衬参数值如表2.2所示，隧道初衬厚度为30cm。
表2.2 围岩和初衬参数取值弹性模量/GPa泊松比密度/Kg/m3内摩擦角/°内聚力/MPaV级围岩 1.80.382200250.3初衬250.22400//文献指出，取模型边界为隧道开挖直径8倍时，其地震和静力计算对隧道的影响可以忽略不计。
本模型取隧道开挖直径的10倍，力求将边界影响降低到最小，同时兼顾计算工作量。
建模时先简化山体坡脚30°一侧为一顶角120°的圆锥，按照隧道实际位置确定隧道后，在山体上切下符合隧道开挖直径10倍的计算模型（如图2.2），这样会使模型上表面近似山体的曲面，较符合实际地形条件，模型自前向后偏压角度略有减小。
与直接将上表面简化为平面的模型相比，其计算结果和实际工程的符合性有较大提升，对隧道内部应力的分布也反映的更为准确。
本章将对无渗流条件下的静力偏压情况进行模拟，以便和后文中渗流作用下偏压隧道的应力场变化进行比较。
模型上隧道圈线和穿越隧道的斜线是为方便后续渗流计算所做分割面，不影响计算。
图2.2 模型切取示意图模型尺寸如上文所述，首先不考虑渗流场，仅考虑重力作用。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ABAQUS岩土隧道入门地应力平衡基本问题
一：Abaqus地应力平衡方法理解
Abaqus地应力平衡现常用分为两种方法：（6.10版以前那种笨拙修改csv文件和添加keywords自己计算每层土应力的方法，就真的很折腾，而且适应性还不好）
1.通过Geostatic中Automatic平衡，这种方法是自动平衡，通过设置最小位移精度，迭
代计算达到平衡的最小位移精度；算盘放小胖版主案例说法“依据小胖的经验，对于标准的隧道开挖，几何简单，采用1e-5的位移准则是可以的。
但如果比较复杂的模型，宝宝们也不要太吝啬，放宽到0.5 mm以下也是可以的。
毕竟我们玩的是大尺度模拟，半个毫米都不一定能测得出来。
”，而对于单元数量巨大的模型，本身计算一次就需要不短的时间，再通过迭代自动计算地应力平衡，这个时间。
；并且Automatic只有100个increment，如果100到了还没平衡好，虽然“可以在上次计算的应力基础上再平衡一次。
”但是这样下去如果遇见不收敛，就是何年何月才能算好地应力平衡。
So，再看第二种方法咯
2.通过导入自重变形结果odb文件，定义应力场来计算。
若用Geostatic分析步换成Fixed
平衡固定计算一次的自重变形结果odb；再通过Predefined field导入这个自重变形结果odb；具体小胖版主有实例截图，就是在Predefined field中initial分析步stress进行设置“从外部文件导入”这个自重变形结果，Geostatic分析步Fixed 下也只有一步step 和Increment，所有就应是1。
但是再计算七万别忘啦，再复制重命名或者新建一个job，要是覆盖了自重变形结果的odb，就白搭
若采用的static general分析步，控制初始增量步默认是Automatic，初始和最大控制一步计算就行，与geostatic的fixed的自重变形计算结果是一样滴，导入平衡方法也一样，结果可以比较一下，具体可以看看算盘坊小胖的地应力平衡初、中、高教程，这个只能算基础入门的概念理解
二：keywords语言基本单元生死法以及场变语言理解
（注：自用的6.14-1和新版，全部都可以通过UG窗口界面设置逻辑了，不用向上世纪那样自己编语言，当
然后期还是要学习一下编程）
▲原始方法地应力平衡语句：
*initial conditions, type=stress, geostatic
土层名,大力，竖向坐标，小力，竖向坐标，侧向系数
▲单元生死（可通过相互作用控制）：
*model change, add
XXXX
*model change, remove
XXXX
▲控制场变：
*field, variable=1
XXXX,2
三：文件存储路径问题
未汉化英文版是不能读取中文路径的，中文路径会有乱码
汉化版中路径文件夹命名不能带类似上面“.”的标点，比如odb文件这样就无法调用Job文件再命名可以用日期加“-”来隔断
调用odb文件需要在temp工作目录下。