地应力平衡方法

方法1

1）建立模型，材料，分析步（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,S23

10)另存为cc.dat

11)在Edit keywords中材料属性后面加上

*initial conditions,type=stress,input=cc.dat

12）重新提交JOB，OK

方法2

1）地表水平、土体材料在水平方向相同，可应用这种简单方法。

2）在Edit keywords中材料属性后面加上。

*initial conditions,type=stress,geostatic

set-1,0.0,5,-392e3,-5,0.9

3）单元集名称、应力竖向分力第一个值、对应垂直坐标、应力竖向分力第二个值、对应垂直坐标、侧压力系数。

4）水平地应力通过竖向应力乘以侧压力系数得到。

补充

6.10及6.11可以实现自动地应力平衡

自动地应力平衡是新版本最为关注的新功能之一，因为它省去了计算自重应力以及生成相应初应力文件和导入的麻烦。在地应力步中选择自动增量步就能使用自动地应力平衡功能，还能指定允许的位移变化容限。不过自动地应力平衡功能仅支持有限的几种材料，D-P并不包含在内，而且对单元也有一定的要求。虽然可以使用不支持的材料和单元，但可能自动地应力平衡不容易收敛或位移差值超过容限。虽然可以用塑性模型，但帮助文件中说应该用在主要为弹性的情况下。我认为材料限制应该不算太大问题，D-P仍可以使用，即使不收敛只要做一些调整比如减小容限等应该一样可以得到收敛的平衡状态。