地应力平衡总结
地应力平衡odb导入法详解
地应力平衡odb导入法详解在工程建设中，地应力平衡是一个非常重要的概念，它是指地下岩体内各点受到的应力相等，即地下岩体处于一种平衡状态。
地应力平衡对于岩土工程的稳定性、地下水的流动、地下开采等方面都有着重要的影响。
因此，准确地预测地应力平衡状态对于岩土工程的设计和实施具有非常重要的意义。
在岩土工程领域，常常需要对地下岩体进行数值模拟分析，以预测其稳定性和变形情况。
为了准确地预测地应力平衡状态，需要对岩体的物理性质、地质条件、地下水等因素进行全面的考虑，并采用适当的数值模拟方法。
其中，odb导入法是一种非常有效的数值模拟方法，可以很好地模拟地下岩体的应力平衡状态。
odb导入法是一种基于有限元分析的方法，它可以将实测的地下应力数据导入到有限元模型中，以模拟地下岩体的应力平衡状态。
具体来说，该方法的步骤如下：1. 对地下岩体进行野外实测，获取地下应力数据。
2. 建立岩体的有限元模型，包括岩体的几何形状、物理性质、边界条件等。
3. 将实测的地下应力数据导入到有限元模型中，作为岩体的边界条件。
4. 进行有限元分析，计算地下岩体的应力分布情况。
5. 根据计算结果，进行岩体的稳定性分析和变形预测等工作。
与传统的有限元分析方法相比，odb导入法具有以下优点：1. 可以很好地模拟地下应力平衡状态，提高分析结果的准确性。
2. 可以充分利用实测数据，避免了传统有限元分析中需要进行大量假设的情况。
3. 可以减少计算量，提高计算效率。
4. 可以更好地考虑地下水对岩体应力的影响。
在实际应用中，odb导入法已经被广泛应用于岩土工程领域。
例如，在地下隧道工程中，需要对隧道周围的地下岩体进行稳定性分析。
采用odb导入法，可以很好地模拟地下岩体的应力平衡状态，预测隧道周围岩体的变形情况，从而为隧道的设计和施工提供重要的依据。
总之，odb导入法是一种非常有效的数值模拟方法，可以很好地模拟地下岩体的应力平衡状态。
在岩土工程领域的应用前景非常广阔，可以为工程设计和实施提供重要的技术支持。
地应力平衡方法介绍
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年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
地应力平衡odb导入法详解
地应力平衡odb导入法详解随着科技的不断发展，计算机仿真技术在地质力学领域的应用越来越广泛。
在地应力平衡方面，有一种常用的方法是odb导入法。
本文将详细介绍该方法的原理、步骤和应用。
一、原理地应力平衡odb导入法是基于ABAQUS软件进行模拟的。
其基本原理是通过有限元分析，计算出地下岩层的应力分布情况，然后根据这些应力分布情况，结合地质条件和地下水情况，确定地下岩层的稳定状态，从而达到地应力平衡的目的。
二、步骤地应力平衡odb导入法的步骤如下：1、模型建立首先，需要根据实际工程情况，建立地下岩层的3D有限元模型。
模型中需要包括地质条件、岩层结构、地下水情况等因素，以及岩层内部的裂隙、孔洞等细节。
2、材料参数设定根据实际的岩层材料，需要设定相应的材料参数。
包括杨氏模量、泊松比、黏聚力、内摩擦角等。
3、加载条件设定根据实际工程情况，设定相应的加载条件。
包括荷载类型、荷载大小、荷载作用时间等。
4、求解根据上述设定，进行有限元计算，求解出地下岩层的应力分布情况。
5、结果分析根据计算结果，进行结果分析。
包括应力分布图、应力云图、位移云图等。
6、地应力平衡判断根据分析结果，结合实际的地质条件和地下水情况，判断地下岩层的稳定状态，从而达到地应力平衡的目的。
三、应用地应力平衡odb导入法在地质工程领域中有广泛的应用。
主要应用于以下几个方面：1、岩石工程稳定性分析地应力平衡odb导入法可以应用于岩石工程稳定性分析中，如隧道、坑道、地下室等工程。
通过模拟计算，可以得到地下岩层的应力分布情况，从而判断岩石工程的稳定性。
2、油气储层稳定性分析地应力平衡odb导入法可以应用于油气储层稳定性分析中。
通过模拟计算，可以得到地下岩层的应力分布情况，从而判断油气储层的稳定性。
3、地震灾害分析地应力平衡odb导入法可以应用于地震灾害分析中。
通过模拟计算，可以得到地下岩层的应力分布情况，从而判断地震灾害的危险程度。
4、地下水工程稳定性分析地应力平衡odb导入法可以应用于地下水工程稳定性分析中。
地应力平衡
地应力平衡地应力平衡，也称作“地质力场平衡”，是指地应力在一定区域内的均衡状态。
它是地球物理学中的一个主要原理，是地质结构以及地球物理效应的基础。
例如，地壳变形过程中的形变机理，地下水在地层内部流动的路径，以及各种地球物理现象是形成地应力平衡的基础。
地应力平衡是一个相对静态的概念。
它主要涉及到地壳构造、岩石结构和地表形态特征之间的关系，以及这些特征之间的力学相互作用。
它的核心思想是，地应力在一定的距离范围内，是一个相对稳定的力学状态。
当地应力强度发生变化时，它们会很快地调整自身的强度，以便保持平衡，而不会出现长期的不稳定状态。
地应力平衡也可以说是地震或地质活动的前兆，它是地应力能量发生变化时发生的一种重要物理反应。
准确地认识和分析地应力变化有助于提高预测地震的准确性和及时性，并有助于开发更有效的震防预报技术。
地应力平衡与地球物理学实验有着密切的联系。
地球物理实验通过量测岩石特性，如晶格参数、表面应力，以及其他物理属性来获取地应力的知识。
这些实验可以揭开地壳深处的秘密，有助于我们更好地了解地壳的变形过程，从而更好地控制地震的影响。
基于地应力平衡理论，地质学家和工程师创建了复杂的模型，用于评估地表物理现象的影响范围，以及地震发生的可能性。
地应力模型允许人们估算地震的影响范围、发生概率、以及可能造成的损失或伤害。
这些模型还能够反映地质构造的变化，以及不同地表学特征的影响，从而帮助地质学家对地质地质结构进行更为准确的分析。
此外，地应力平衡理论也被用于研究地质工程、矿物开采和矿山建设中的关键物理因素。
例如，开采过程中排放的气体会改变地应力，从而改变岩石的结构和组织，可能会导致岩体受到破坏。
在这种情况下，可以利用地应力模型来估算和模拟地应力的变化过程，帮助预测开采超过允许水平的可能性，从而保护地表面和地下环境。
综上所述，地应力平衡是一个十分重要的概念，它是地球构造变形和地表物理特征发生变化时出现的一种物理反应。
地应力平衡汇总
1、“地应力平衡”的含义、目的、作用我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
2、地应力如何平衡地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
abaqus的part模块绘图功能不是很强，因此常用AutoCAD绘出平面图后导入abaqus。
在abaqus6.10中，主导入要分为以下几个步骤：1，在AutoCAD中建好模型的平面图形，并且另存为dxf格式。
地应力及其分布规律
地应力及其分布规律1 、地应力的根本概念地应力是存在于地层中的未受工程扰动的天然应力，也称岩体初始应力、绝对应力或原岩应力。
广义上也指地球体内的应力。
它包括由地热﹑重力﹑地球自转速度变化及其他因素产生的应力。
地应力是各种岩石开挖工程变形和破坏的根本作用力；是确定工程岩体力学属性，进行围岩稳定性分析，实现开挖设计和决策科学化的必要前提条件。
此外地应力状态对地震预报、区域地壳稳定性评价、油田油井的稳定性、核废料储存、岩爆、煤和瓦斯突出的研究以及地球动力学的研究等也具有重要意义。
2、地应力的成因产生地应力的原因是十分复杂的，地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆浸入和地壳非均匀扩容等。
另外，温度不均、水压梯度、地表剥蚀或其它物理化学变化等也可引起相应的应力场。
其中，构造应力场和自重应力场为现今地应力场的主要组成局部。
当前的地应力状态主要由最近的一次构造运动所控制，但也与历史上的构造运动有关。
由于亿万年来，地球经历了无数次大大小小的构造运动，各次构造运动的应力场也经过屡次的叠加、牵引和改造，另外，地应力场还受到其他多种因素的影响，造成地应力状态的复杂性和多变性，地应力成因之一：地幔热对流〔图1、图2〕地应力成因之一：板块边界受压〔图3〕地应力成因之一：岩浆浸入〔图4〕3、地应力的影响因素地壳深层岩体地应力分布复杂多变，造成这种现象的根本原因在于地应力的多来源性和多因素影响，但主要还是由岩体自重、地质构造运动和剥蚀决定。
1)岩体自重的影响岩体应力的大小等于其上覆岩体自重，研究说明：在地球深部的岩体的地应力分布根本一致。
但在初始地应力的研究中人们发现，岩体初始应力场的形成因素众多，剥蚀作用难以合理考虑，在常规的反演分析中，通常只考虑岩体自重和地质构造运动2)地形地貌和剥蚀作用对地应力的影响地形地貌对地应力的影响是复杂的，剥蚀作用对地应力也有显著的影响，剥蚀前，岩体内存在一定数量的垂直应力和水平应力，剥蚀后，垂直应力降低较多，但有一局部来不及释放，仍保存一局部应力数量，而水平应力却释放很少，根本上保存为原来的应力数量，这就导致了岩体内部存在着比现有地层厚度所引起的自重应力还要大很多的应力数值。
地应力平衡
ABAQUS地应力平衡：进行地应力平衡的原因陈述如下:我们建立的几何模型一般都和工程实际情况一致,例如边坡的几何模型与边坡实际尺寸相一致。
但是由于边坡的沉降和徐变作用，可以想像到，现在的边坡应该是由一个体积更大的原始边坡在很久以前由于受到重力作用和边界约束条件,逐渐形成了现今的边坡形态.但是对于那个原始的边坡形态，我们不得而知。
假如能准确知晓，我们就能够建立原始边坡的几何模型，接着对边坡施加重力和边界条件,受力后边坡形态应该和现在的边坡相一致，其内力就是初始应力场(地应力），这样就不用专门施加地应力了。
但现实情况是我们不能知晓原始边坡的形态.现在的边坡几何模型就是其实际形态，受力之后将会变成一个与现状不一致的边坡，这不符合现在的实际情况.如果我们计算出现今边坡的内力，并将其作为边坡的初始应力场，再去和外力平衡，这样我们建立的模型就和现实边坡情况相一致了。
对于涉及开挖、回填的动态岩土工程问题，地应力平衡是正确模拟施工过程的前提条件。
初始应力的加载必须满足地应力平衡,而地应力平衡就是为了使地基仅存在初始应力，而不存在初始应变。
当地基自重是产生地应力场的主要因素时,重力是外力，初始应力场是内力，将提取出的内力施加于模型后再施加重力,此时内力和外力平衡，该状态就是工程建设的初始状态。
在表面水平的情况下，ABAQUS中初始地应力场的平衡一般只和密度有关,土体的密度一样，平衡的效果就好，别的参数对地应力平衡的结果影响很小。
对于表面不平的情况,尽量通过inp文件导入初始应力的方法进行地应力平衡。
ABAQUS中进行地应力平衡的时间点的选择十分重要,地应力平衡是指在工程建设之前，地表的位移应为零, 而土体的应力却存在。
也就是说不管土体原来的样子如何（例如高山, 河流，丘陵，平原等), 进行地应力平衡的正确时间点应当是在我们对它做任何扰动之前。
具体采取的办法如下所述，我们对所建立的边坡几何模型施加和实际模型一致的重力和边界条件，得到变形以后边坡的内力，变形后边坡形状和原始边坡略有不同，其内力可近似作为现状边坡的内力,将其作为初始应力施加于现在的边坡中，接着施加外力(重力）来平衡初始应力，这样就建立了一个与现今边坡形态基本相同的边坡模型，这样之后的分析计算才是符合实际的。
完整的地应力平衡方法
地应力平衡方法：第一步：建立模型，材料，分析步（GEOSTATIC）第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8 第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON) (请严格按照这个格式，注意大小写的字母)第四步：提交J0B，完成后第五步：按以下步骤,Roport---Report Field Output---选中S11,S22,S33,S12,S13,S23---Name:XX.INP---Write中选择Field Output-------------ok!!!第六步：用软件(推荐使用UltraEdit很好编辑的)打开XX.INP，保存格式内容单元号 S11 S22 S33 S12 S13 S23 （请注意，在保存内容中没有这一行的）1 , . , . , . , . , . , .2 , . , . , . , . , . , .. , . , . , . , . , . , .. , . , . , . , . , . , .这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些数字。
第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上：*initial conditions,type=stress,input=xx.inp 完成第八步：重新提交JOB，OK第九步：如果你还没有成功的话，那我只能说----------------我无语了。
：）ABAQUS的这项功能确实很不错。
：）这个功能让基坑开挖、隧道开挖等的初始应力，开挖后的残余应力很好的显示；也可以很好的模拟铁路设计中的工后沉降的概念，在地应力平衡后，加上荷载所得沉降即为工后沉降；也很好的模拟了桩土复合地基的问题，如果没有初始应力的模拟，使土对桩产生了挤压应力，从而通过设定摩擦系数就可以模拟了桩与土之间的摩擦力；除此之外，在进行挡土墙计算时也需要ABAQUS的这项功能，反正很多都用得着。
[测试成功]地应力平衡的一个简单例子(1)
![[测试成功]地应力平衡的一个简单例子(1)](https://img.taocdn.com/s3/m/d7fdd43e83c4bb4cf7ecd165.png)
地应力平衡方法熊志勇陈功奇第一部分地应力平衡方法简介地应力平衡有三种方法：（1）*initial conditions,type=stress,input=FileName.csv(或inp)该方法中的文件FILENAME.INP获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。
但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。
在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应,因此平衡后的效果不一定很理想,但无疑其适用性很强。
（2）*initial conditions,type=stress,geostatic该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。
所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。
计算速度快,收敛性好。
缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。
（3）*initial conditions,type=stress,geostatic,user该方法采用用户子程序SIGINI来定义初始应力场,可以定义其为应力分量为坐标、单元号、积分点号等变量的函数,要达到精确平衡需已知具体边界条件,在实际中应用较少。
第二部分 地应力平衡方法实例详解地应力平衡是岩土工程数值模拟分析的重要的内容，为了让师弟师妹们快点上手，我利用第一种方法做一个较简单的模型，希望对大家有用。
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年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
地面稳定性力学知识点总结归纳
地面稳定性力学知识点总结归纳
地面稳定性力学是研究地面工程中地基和岩土的稳定性问题的学科。
以下是地面稳定性力学的一些主要知识点的总结归纳：
1. 地基和岩土的力学性质：
- 地基和岩土的重要力学性质包括黏聚力、内摩擦角、孔隙比等。
- 地基和岩土的力学性质与其成分、颗粒大小和排列方式等因素有关。
2. 地基和岩土稳定性的判断：
- 地基和岩土的稳定性指其抵抗外力作用而不产生破坏或变形的能力。
- 进行地基和岩土稳定性判断时，需要考虑土层的强度、应力状态和外部荷载等因素。
3. 外部荷载的影响：
- 外部荷载是指施加到地基和岩土上的力，包括静力荷载和动力荷载。
- 外部荷载的大小和作用方式对地基和岩土的稳定性有重要影响。
4. 常见的地基和岩土稳定性问题：
- 地基沉降：地基沉降是指地基下沉或沉降引起的地面变形。
- 坡面稳定性：坡面稳定性是指坡面的抗滑性能。
- 基础稳定性：基础稳定性是指建筑物或结构物的基础抗倾覆或抗沉降能力。
- 地基裂缝：地基裂缝是指地基中出现的裂缝或断裂。
5. 地面加固与处理方法：
- 对于地基和岩土稳定性问题，可以采取不同的加固和处理方法，包括填土加固、挡土墙、地下支护等。
以上是地面稳定性力学的一些主要知识点的总结归纳。
地面稳定性力学的研究可以帮助我们理解地基和岩土的力学性质，并为地面工程中的土壤处理和加固提供指导。
地应力平衡方法
地应力平衡方法方法11）建立模型，材料，分析步（GEOSTATIC）。
2）施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8。
3）在JOB中提交分析。
4）按以下步骤,Report---Report Field Output---选中S11,S22,S33,S12,S13,S23---Name:cc.inp。
Write中只选择Field Output。
5）修改cc.inp，用excel，打开（分隔符，Tap键、空格键、逗号）6）删除都是1的那列。
在1,2,3,4等的前面加上（part instance）的name和小数点。
7）另存为，文件类型设置为CSV。
8）用文字编辑软件删除小数点后面的逗号。
9）最后变为soil-1.1,S11,S22,S33,S12,S13,S2310)另存为cc.dat11)在Edit keywords中材料属性后面加上*initial conditions,type=stress,input=cc.dat12）重新提交JOB，OK方法21）地表水平、土体材料在水平方向相同，可应用这种简单方法。
2）在Edit keywords中材料属性后面加上。
*initial conditions,type=stress,geostaticset-1,0.0,5,-392e3,-5,0.93）单元集名称、应力竖向分力第一个值、对应垂直坐标、应力竖向分力第二个值、对应垂直坐标、侧压力系数。
4）水平地应力通过竖向应力乘以侧压力系数得到。
补充6.10及6.11可以实现自动地应力平衡自动地应力平衡是新版本最为关注的新功能之一，因为它省去了计算自重应力以及生成相应初应力文件和导入的麻烦。
在地应力步中选择自动增量步就能使用自动地应力平衡功能，还能指定允许的位移变化容限。
不过自动地应力平衡功能仅支持有限的几种材料，D-P并不包含在内，而且对单元也有一定的要求。
虽然可以使用不支持的材料和单元，但可能自动地应力平衡不容易收敛或位移差值超过容限。
地应力平衡的一个简单例子
abaqus地应力平衡先说为什么要施加地应力：1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
这就是所谓“地应力平衡”的含义、目的、作用。
2. 地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
3. 地应力场的方向问题，有网友在论坛里问，既然重力是向下，为与重力平衡，那应力场的方向是不是向上呢，这同样是我开始接触abaqus 的疑问，相信很初学者也有这样的疑问，我的理解是内力是没有向上、向下或者向其它方向的概念的，内力只有拉力或压力或剪力之分，其方向也按是拉是压是顺时针或逆时针而分，内力往往都是成对出现，如地应力场中的应力以压应力为主，取一个微元，则压应力同时出现在向下和向上，你能说地应力就是向上，与重力反向吗？aba中初始地应力场平衡一般在表面水平的情况下仅仅和密度相关，密度一样的话平衡的结果很好，别的参数改变之后经过计算，差别很小。
ABAQUS应用于土石坝数值计算的技巧总结
ABAQUS应用于土石坝数值计算的技巧总结作者：dailinghuiABAQUS软件是国际著名的有限元通用软件，功能非常强大，但是在进行土石坝应力应变计算这一专业而又复杂的问题上，还存在一些局限性。
因此，必须对软件和计算问题有深入的了解，才能找到解决办法，突破局限性。
为了便于沟通和交流，本文对ABAQUS软件进行土石坝计算时的一些经验技巧加以总结和探讨。
共分为初始建模、网格划分、地应力平衡、无限元节点调整、增量步设置、计算过程控制、结果后处理几个部分进行总结。
一、初始建模1、ABAQUS CAE建模：ABAQUS软件的前处理功能比较强大，可以通过人性化的界面来建立几何模型。
但是广大工程人士最为了解和熟悉的画图软件当属AutoCAD。
ABAQUS软件有专门的接口，可以将格式为dxf的CAD文件导入进来。
在导入的过程中，可以采用以下步骤来进行：（1）首先将原始图形（坝体剖面图）导到Skectch模块；（2）在Part模块中建立Part（坝体）的时候，由于ABAQUS只允许外边框存在，而不能有内部线，因此，可以先把Skectch模块中的坝体剖面图全部添加进来，只保留坝体的外边框，将其余的线条都删去；（3）对Part进行Partition分割操作：把坝体剖面图全部添加进来，只保留分割坝体所需要用到的线条，把其他多余线条全部删去，包括坝体的外边框。
建第二个、第三个Part的时候，重复以上操作即可。
采用上述方法建模的好处是，整个建模过程用到的所有线条，都是基于原始图形（坝体剖面图）的，因此可以避免线条有交叉、结点不重合等情况的发生。
而且在Assembly模块里将各个Part 进行组装的时候，不但不需要再调整各个Part位置，还能保证各个Part之间接触面完全重合，为建立正确的接触关系打下基础。
2、其他前处理软件建模：ABAQUS软件的前处理功能是基于“自上而下”的思路设计的，即先建出模型的整体，然后根据不同部位的要求切出各个小块。
地应力平衡图解
*initial conditions,type=stress,input=0.csvmdb.models['Model-1'].setValues(noPartsInputFile=ON)（将网格划分精密一些，地应力平衡精度更高一些，2010-7-12注）注意：导入的inp、csv或其他文件里不得有空行，否则会出现element 0。
允许有空格。
--------------------------------------------------------------------方法一：桩土地应力平衡不容易通过，可先直接指定应力：*initial conditions,type=stress,geostaticSet-pile-soil,0,0,-510000,30,0.6,0.6将桩土摩擦系数设为0，容易平衡，输出含有应力s11、s22等的应力文件；导入应力文件，将摩擦系数提高为一个较小的数值（低于正常的摩擦系数），再计算，输出应力文件；再提高摩擦系数，导入应力文件，计算，重复上述操作，直到摩擦系数达到正常值。
上述如果不行，可以干脆将桩与的接触改为tie，2010-8-18的模型按照上面的步骤操作不行，后来tie就可以了，且收敛的速度很快。
地应力平衡方法：第一步：建立模型，材料，分析步（GEOSTATIC）第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8Y第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON) (请严格按照这个格式，注意大小写的字母) ，例如：mdb.models['Model-1'].setValues(noPartsInputFile=ON)第四步：提交J0B，完成后第五步：按以下步骤,Roport---Report Field Output---选中S11,S22,S33,S12,S13,S23---Name:XX.INP ---Write中选择Field Output-------------ok!!!第六步：用软件(推荐使用UltraEdit很好编辑的)打开XX.INP[注：先用UltraEdit打开，删除中间的英文，再用excel打开]，保存格式内容单元号 S11 S22 S33 S12 S13 S23 （请注意，在保存内容中没有这一行的）1 , . , . , . , . , . , .2 , . , . , . , . , . , .. , . , . , . , . , . , .. , . , . , . , . , . , .这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些数字。
地应力平衡汇总
1、“地应力平衡”的含义、目的、作用我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
2、地应力如何平衡地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
abaqus的part模块绘图功能不是很强，因此常用AutoCAD绘出平面图后导入abaqus。
在abaqus6.10中，主导入要分为以下几个步骤：1，在AutoCAD中建好模型的平面图形，并且另存为dxf格式。
地应力平衡方法 (自动保存的)3页
地应力平衡方法：. c9 Z A4 D6 k$ h 第一步：建立模型，材料，分析步（GEOSTATIC）( z4 A! K! _2 p' B% r( O6 o2 f; {第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8SimWe" F( ]& A3 q4 [8 [第三步：在命令行中输入mdb.models['Model-1'].setValues(noPartsInputFile=ON ) (请严格按照这个格式，注意大小写的字母)仿" T#R$ m/ W" p2 v/ ?第四步：提交J0B，完成后/ J. q5 _) p8 a"第五步：按以下步骤,[Module]Visualization-Report---Report FieldOutput---下拉菜单里面选择2 w5 n# Z2 U0 v( w centriod，然后依次把s11、s22、s33、s12、s13、s23点选上，setup页面选择报告文件的名字abaqus.csv---Write 中选择Field Output-------------ok!!!% F2 d: d8 u5 p- M) T第六步：, r% V# }& G5 d打开abaqus.csv，去掉英文部分，修改成下面形式（文件内容格式如下）5 U/ d; ]/ {( k" S) A单元号S11 S22 S33 S12 S13 S23 （请注意，在保存内容中没有这一行的）9 C; S4 e3 P" T6 d1 ~3 g0 U:R 1 , . , ., . , . , . , , ~: |1 e" o9R) m2 , . , . ,. , . , . , % w( y" ~0_7 d& r& S+ ~. , . , ., . , . , . ,* w& Y+ j$ F7 j2 l- l, F9 `! ^0 ]0A . , . , ., . , . , . , : x&O5 ?/ [ |# s, \& B这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些这些数字这些数字。
5.1.2地应力的分布规律
地形、地表剥蚀、风化、 岩体结构特征、岩体力学性质、 温度、地下水等因素
中国现代构造应力场图
（1）原岩应力场是相对稳定的非稳定场
地应力一般是以水平应力为主的三向不等压应力场。 主应力随空间和时间而变化，是非均匀应力场。 在空间上从小范围来看，应力空间分布不均。 整个区域而言，在长时间序列下相对稳定。
深度H/m
垂直应力σv/MPa
0
10
20
30
40
50
60
70
力γH,即：σv＝27H
在埋深小于1000m时，测量值与预测值可能差别 很大，有的甚至相差达到5倍。
深度H/m
平均水平应力与垂直应力的比值σh,av/σv
0
0.5 1.0 1.5 2.0 2.5
3.0 3.5
500 1000
k=1500/H+0.5
1500
2000
2500 3000
k=100/H+0.3
500 1000 1500 2000 2500 3000
各国垂直应力σv分布图
ห้องสมุดไป่ตู้
（2）实测垂直应力σv基本等于上覆岩层的重力
E.Hoek和E.T.Brown总结出的实测垂直应力随深 度H变化的规律，如下：
在25m～2700m深度内，垂直应力σv 随深度线性 增加，大致等于按平均重度γ=27kN/m3计算出来的重
平均水平应力与垂直应力的比值
（3）实测的水平应力普遍大于垂直应力
侧压力系数 λ＝σh/σv >1（构造应力起控制作用）
（4）平均水平应力σh与垂直应力σv的比值随 深度增加而减小
（5）最大、最小水平主应力随深度线性增加 σHmax＝6.7＋0.0444H (MPa) σHmin＝0.8＋0.0329H (MPa) （6）两个水平应力的关系一般相差较大
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
地应力平衡
1、地应力平衡好坏评判标准
1)地应力平衡后，位移云图中最大位移达到10-6量级或更低（接近于0）。
（主要判别条件）
2)地应力平衡后，应力云图中应力有一定的数值。
（也就是应力不为0，但变形接近于0）2、进行地应力平衡的原因
总的来说，如果不进行地应力平衡，而只施加重力，模型会在重力作用下产生变形，而实际工程中，我们施加荷载时，重力产生的而变形已经产生，实际上得到的是附加应力产生的变形。
1)我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型
和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。
真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
这就是所谓“地应力平衡”的含义、目的、作用。
2)地应力平衡中的外力和内力的问题。
地应力平衡中，显然，重力是外力，应力场是内力，
仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。
为的是内力和外力平衡。
) q0 F3 q6 H1 O#
3)地应力场的方向问题，有网友在论坛里问，既然重力是向下，为与重力平衡，那应力场
的方向是不是向上呢，这同样是我开始接触abaqus的疑问，相信很初学者也有这样的疑问，我的理解是内力是没有向上、向下或者向其它方向的概念的，内力只有拉力或压力或剪力之分，其方向也按是拉是压是顺时针或逆时针而分，内力往往都是成对出现，如地应力场中的应力以压应力为主，取一个微元，则压应力同时出现在向下和向上，你能说地应力就是向上，与重力反向吗？不怕各位笑话，以上几点在高手看来是很简单的问题，却是我经历了漫长而艰辛的摸索才得到的，今天也写给初学者，不要再走我的老弯路了。
aba中初始地应力场平衡一般在表面水平的情况下仅仅和密度相关，密度一样的话平衡
的结果很好，别的参数改变之后经过计算，差别很小。
表面不水平的情况则最好通过文件导入初应力的情况进行平衡。
一般来讲，表面不平的时候有很多因素造成误差很大。
3、地应力平衡的方法
地应力平衡有三种方法：
（1）*initial conditions,type=stress,input=FileName.csv(或inp)
该方法中的文件FILENAME.INP获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。
其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。
但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。
在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应,因此平衡后的效果不一定很理想,但无疑其适用性很强。
第一步：建立模型，材料，分析步（GEOSTATIC）
第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8Y 第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)(请严格按照这个格式，注意大小写的字母)
例如：mdb.models['Model-1'].setValues(noPartsInputFile=ON)
第四步：提交J0B，完成后
第五步：按以下步骤,Roport---Report Field Output---选中S11,S22,S33,S12,S13,S23---Name: XX.INP --- Write中选择Field Output-------------ok
第六步：用软件(推荐使用UltraEdit很好编辑的)打开XX.INP[注：先用UltraEdit打开，删除中间的英文，再用excel打开]，保存格式内容
单元号S11S22S33S12S13S23 （请注意，在保存内容中没有这一行的）1,.,.,.,.,.,.
2,.,.,.,.,.,.
.,.,.,.,.,.,.
.,.,.,.,.,.,.
这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2........这些数字。
例如：
}
另存为0.csv
当没有进行第三步操作时，前面的序号前应加上instance的名称。
见《ABAQUS有限元分析常见问题解答》P126。
第七步：在ABAQUS----Model---Edit keywords---Model-1(这就是你的Model名字)---在材料属性后面加上：*initial conditions,type=stress,input=0.csv
第八步：重新提交JOB，OK
第九步：如果你还没有成功的话，那我只能说----------------我无语了。
ABAQUS的这项功能确实很不错。
这个功能让基坑开挖、隧道开挖等的初始应力，开挖后的残余应力很好的显示；也可以很好的模拟铁路设计中的工后沉降的概念，在地应力平衡后，加上荷载所得沉降即为工后沉降；也很好的模拟了桩土复合地基的问题，如果没有初始应力的模拟，使土对桩产生了挤压应力，从而通过设定摩擦系数就可以模拟了桩与土之间的摩擦力；除此之外，在进行挡土墙计算时也需要ABAQUS的这项功能，反正很多都用得着。
（2）*initial conditions,type=stress,geostatic
该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。
所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。
计算速度快,收敛性好。
缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。
1）地表水平、土体材料在水平方向相同，可应用这种简单方法。
2）在Edit keywords中材料属性后面加上。
*initial conditions,type=stress,geostatic
instance名称.单元集名称,0.0,5,-392e3,-5,0.9
3）单元集名称、应力竖向分力第一个值、对应垂直坐标、应力竖向分力第二个值、对应垂直坐标、侧压力系数。
水平地应力通过竖向应力乘以侧压力系数得到。
（3）*initial conditions,type=stress,geostatic,user
该方法采用用户子程序SIGINI来定义初始应力场,可以定义其为应力分量为坐标、单元号、积分点号等变量的函数,要达到精确平衡需已知具体边界条件,在实际中应用较少。