地应力平衡图解(2)

以上3个步骤可以代替以下命令： *el print s,

手工输入地应力平衡

*initial conditions,type=stress,geostatic Set-pile,0,30,-500000,10,0.2,0.2 set-soil,0,30,-540000,0,0.6,0.6

对于很厚的土层必须先进行地应力平衡，否则有些单元会出现过大的扭曲变形，最后运算不了，出现ERROR: FIXED TIME INCREMENT IS TOO LARGE（即使在step中打开非线性开关也不行），必须先手工输入地应力平衡初始条件，见上面。在inp中加入 *el print s,

将各单元的应力导出（*.csv），再*initial conditions,type=stress,input=*.csv导入初试应力（程序自动将取排在后面的初始应力为准，自动把第一个初始应力除去如下图，*initial conditions,type=stress,geostatic和下面一行自动被除去），这样迭代平衡效果非常好。

对于采用c3d8等非缩减积分单元，会出现如下图的应力文件：

用下面的方法就OK：（不要选择积分点，而是质心）