ABAQUS 如何施加残余应力

abaqus自带有地应力平衡功能，用GEOStatic就能完成，但是比较多的做法是利用初始应力方法实现的，网上有很多具体做法，如：

1.先施加重力荷载的作用,可以在cae中实现2.在inp文件中的output

request中写上*el

prints,这样就会将施加重力荷载后的应力输出到*.dat文件中了

3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如ELEMENT

T

FOOT-

S11

S22

S33

S12

单独存为一个*.dat文件,4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance.

,即单元编号变为:

instance名称.序号

注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!5.接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),6.最后在inp文件的step之前写上*initial

conditions,type=stress,input=文件名.csv

有人提出如下问题：

在这里我有两个问题，麻烦了解的告诉我一下：

1.关于第二步分析步选项，选择Static，General和Geostatic，两者都可以计算出重力下的应力，具体区别在哪？或者说选Static不对的原因？

2.我发现缺少第三步照样可以算出来，这是为什么？换句话第三步的目的和必要性是什么？

我的回答如下：

地应力是只有应力没有位移的，采用静力计算会产生位移，利用地应力平衡只需要应力不差生位移。

第三步骤，主要是输出的问题，有没有都没有关系

是没有差别的，两者应该是一样的。Geostatic

那个本来就是地应力平衡的，可能是后续的分析步有限制，所以才用了初始应力的方法吧

abaqus6.12以及以上版本可以在load里面有个预定义场，选择在initial步时设置initial

stress。当然在具体实施中可能会遇到一些无法施加的问题，那样就要针对具体问题进行相应改动。若abaqus版本较低，则在inp文件中增加*initial

condition,type=stress,input=XXX.csv命令行，具体实施参见帮助文档或者去simwe论坛查看。

混凝土密度是2.5E-9的话和重力加速度输入9800的单位不一样（话说2.5E-9是不是太小了点了，你这是咋出来的数据呀），所以你把重力加速度减小能算。虽然ABAQUS里的数值都不用输入单位，但你输入的所有数值的单位都要统一，建议都按国际单位换算后输入，这样不容易出问题。

abaqus 中，如果用“重力GRAV“”分布载荷定义时，abaqus/standard基于总孔压进行分析；若模型重力通过body force来实现，则abaqus/standard采用的是超孔压。

分析中我们关心载荷引起的超孔压分布及消散，因此大多数问题的初始超孔压为0，无需定义；

而若采用总孔压进行分析，则需通过*Initial conditions,type=pore pressure来定义初始总孔压的分布，如静水压力分布。

---