Abaqus速成技巧 - 图文(2)

论坛汇总第一部分

地应力问题

一、在桩土接触分析中，地应力步中，老是提示我this step can't continue by the select increment .我只把*geostatic改成*static就能够计算到30多个increment.为什么会出现这种问题？请各位大虾赐教。 你的geostatic使用不恰当,因为在有接触面的情况下,要注意力的方向。

楼上大哥，我还是不太明白，我的inp文件是由hypermesh生成的，既然方向不对，应该用*static也不能计算阿，请说的明白一点，非常感谢，我现在把我的inp文件传上，望赐教。

268204-1.inp (0.0k)

二、不平地面初始地应力的处理方法

245572-geo3.ps 命令流文件

245574-geo3.inp (0.0k) 我看ABAQUS的帮助里，初始地应力好像是通过给定两点的高程，及相应的z方向应力值、水平方向应力系数来计算单元集的初始地应力，这样对于地面不是水平的怎么办？能不能通过外部文件给单元赋初始地应力,如何操作 可以！！通过inp文件

还有，如果按重力计算出来的应力作为初始应力，请问侧向压力系数怎么体现？急！ 请看本人通过文件进行初始应力场平衡的工程事例（也是刚刚才发出的）吧！ 我有这样一个想法：

我们在一些情况下有了地应力，但是，我们对地表的初始变形不满意，那么我们先在第一步计算出地应力合初始变形，然后我们将地应力写入为一个文件，在这个文件里我们得到单元积分点上保持平衡的应力（或者其他，我还在想）然后我们通过USER SUBROUTINE 或者INPUT读入，那样我们不就得到了应力却没有如何初始变形吗！网上也有这样的帖子，但是过程不是很明了。

我看了帮助，INPUT好像不行，因为它输入的是单元号和几个应力，然后ABAQUS将这几个应力平均的施加到单元的面上。 USER SUBROUTINE 好像可以它输入的是单元号，积分点号，然后就是应力分量。而我们可以在第一步就输出这几个量，只需稍微操作一下，加一些逗号就可以。只是还有一个问题，USER SUBROUTINE 一次只能返回一个单元号一个积分点的应力分量，但是我们需要一批单元的应力分量，暂时还不知道如何处理！！！不知有无大侠知道如何处理！！！ 纠正一点，不是INPUT，是指忽略GEOSTATIC,而直接用*initial conditions,type=stress输入应力分量。 很好的想法，我也在想通过子程序施加，但是ABAQUS的例子太简单，看完 不知所云，：）

我仔细看了帮助，现在面对的难题是：SIGINI里输入的是一个ELEMENT的SRESS，注意是一个ELEMENT而不是一个ELEMENT的一个积分点的STRESS，但是我找了一下，好像找不到一个ELEMENT的STRESS输出，不知兄台有无办法输出！！！如果能找到，就应该能实现目标！

努力值得赞赏，但是你的物质的弹模过高，土的弹模一般达不到GPa的量级，而且如果采用弹塑性本构后引入非线性，收敛情况将会恶化，呵呵

假定你计算的是个土质边坡，假设土的粘结力为零，内摩擦角为20，那么此时边坡的安全系数为1.2。和想和原作者一起探讨边坡稳定性分析的问题 请教zhangpeiwen 兄： 你的INP文件里 *GEOSTATIC 1E-6,1E-6

这个1E-6是什么参数啊？我在HELP里看到*GEOSTATIC里是没有参数的。 是初始时间增量，和最终时间，一些例题里采用这样的参数，应该没有问题．