地应力平衡方法以及注意事项

 注意：只有采用弹塑性本构模型时需要地应力平衡，弹性本构不需要地应力平衡！
第一部分  地应力平衡方法简介

地应力平衡主要有五种方法：

（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,其他四种方法第一步除了土外可以有其他部件，效果一般很差，可以多循环几次！