论坛汇总第一部分

好像就是这么做的，当然水工计算相对说要粗一些，不知说的对不对。不吝赐教！ 四、请教一个初始应力问题，多谢！

我是新手，加了初始应力后，可是没有反应， 不知道什么地方出了问题 请各位大哥分析，谢谢了 ** INTERACTION PROPERTIES **

*Surface Interaction, name=IntProp-1 1., *Friction 0., **

*INITIAL CONDITIONS,TYPE=STRESS bulk, 7800., 7800. soil, 1600., 1600.

** BOUNDARY CONDITIONS **

五、请问如何解决初始地应力问题

请教诸位大侠，在基坑分析当中，如何解决初始地应力的问题？

1，能否先计算自然土层（开挖前）的每个单元节点的竖向位移（重力作用下），然后再在以后的开挖回填荷载步中，将重力作用下的位移场减去

2，或者能否先计算重力作用下每个单元节点的应力状态，再把这个应力场作为初始应力场加在模型上，然后再在这个被加上初始应力场的模型下进行开挖计算

这些方法都是在ansys或是adina里面可以用的，可是我不知道如何在abaqus里面如何实现。请教各位，恳请指点一二，谢谢了

abaqus的documentation里面的standard/manual中的analysis目录下面的coupled porous fluid flow中有关于地应力场如何处理的说明，abaqus中线建立一个初始平衡场，有应力没有位移，用*geostatic *initial conditions, type=stress, geostatic

六、我们在处理土木工程的时候，一般都要遇到初始地应力的处理问题，基本都要使用如下的关键词： *INITIAL CONDITIONS TYPE=STRESS, GEOSTATIC 在手册中队其数据行的解释是如下： Data lines for TYPE=STRESS, GEOSTATIC: First line:

1. Element number or element set label.

2. First value of vertical component of (effective) stress. 3. Vertical coordinate corresponding to the above value. 4. Second value of vertical component of (effective) stress. 5. Vertical coordinate corresponding to the above value.

6. First coefficient of lateral stress. This coefficient defines the X-direction stress components.

7. Second coefficient of lateral stress. This coefficient defines theY -direction stress component in three-dimensional cases and the thickness or hoop direction component in plane or axisymmetric cases. If this value is omitted, it is assumed to be the same as the first lateral stress coefficient given in the previous field.

Repeat this data line as often as necessary to define an initial geostatic stress state in various elements or element sets. 对其中我有几点认识不清的地方，首先我对(effective) stress的理解就是初始地应力，在数据行中2为Y方向的垂直地应力，

- 9 -

论坛汇总第一部分

3为其对应的坐标，但是对4种的解释地应力的垂直分量的第二个值如何理解？7侧应力的第二系数，如何理解，在ABAQUS中Y方向不就是垂直方向吗？如何把垂直方向按照是Z方向来理解，则7中的解释可以理解，也就是两个水平应力的计算中的系数，但是按照手册中是如何理解的呢？我现在没有弄明白，希望有朋友可以帮助我一下，谢谢。

我看ABAQUS的帮助里，初始地应力好像是通过给定两点的高程，及相应的z方向应力值、水平方向应力系数来计算单元集的初始地应力，可以试着*initial conditions,type=stress（即不用geostatic参数）每个单元所有的应力分量值。 Data lines for TYPE=PRESSURE STRESS: First line:

1. Node set or node number. 2. Equivalent pressure stress, .P

楼上的朋友你说的是不是这个啊？如果是的话，也许你的解释有点问题，因为这个参数的数据行只有两个，一个是节点集或节点号，另一个是等效压应力。您可以进一步解释一下你的理解吗？谢谢。 七、【讨论】由初始条件地应力引起的接触问题如何处理?

考虑重力场，需要在初始条件中加地应力,然后Step中施加Gravity与此相平衡，虽然地表产生的沉降不大，但不可避免总有一些沉降。我在地表放一基础，有Assembly将地基和基础方在一起。由于Grvaity产生的沉降，将两者分离，引起计算不能进行。如何解决这个问题？请高手帮助。谢谢

“由于Grvaity产生的沉降，将两者分离，引起计算不能进行。 ” 1如果地基和基础在一个part里，不会产生分离。 2如果是两个part，应该设置接触条件。 3可以把基础单元先remove掉，geosattic后再激活，这样模拟施工顺序要好一点。 4把inp文件传上来看看嘛

八、【原创】哈哈，终于建立地应力平衡了，用CAM-CLAY的平衡

摸索了好几天，都快崩溃了。最后终于解决。之前也参考了我们这网上的一些文章。但总是解决不了。最后我发现原来主要问题在这里:只要你的模型的密度×加速度值×你的模型的重力作用方向的尺度应该与你的*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC里设置的应力相对应，想想：如果不对应模型就不平衡，然后就得去寻求平衡的街！！！我试了一下，改变侧压力系数并不会对初应变有太大的影响。 这是我的模型，CONC是用来以后用的

239202-jk1.inp (0.0k) 这是变形图

- 10 -

论坛汇总第一部分

九、本人用ABAQUS软件对某地铁工程的施工过程进行了三维弹塑性数值模拟分析，现将有关内容显示在附件中。

340189-ABAQUS.doc (0.0k) 看来大家还是很关注这个工程事例，现把其他几个文件给出，想必对大家关于ABAQUS的实际应用会有很多启发和帮助吧！希望斑竹能够给予鼓励，如果希望还会有更精彩的工程应用出现，我等待。。。。。 342333-other.rar (0.0k) 你好，我刚运行了你的文件，老是出错，提示：ERROR: SPECIFIED STANDARD_MEMORY VALUE OF 256 MB IS TOO SMALL TO RUN THE

ANALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE VALUE IS 436 MB. LOOK AT MEMORY ESTIMATES SECTION OF .dat FILE FOR FURTHER INFORMATION.

我的电脑配置是512内存，p4处理器，应该不会有问题吧？

控制面板／管理工具／计算机管理，进入后，在计算机管理上点右键，编辑属性／高级／性能／性能选项，更改虚拟内存，把最大和最小都改的大一些 侧压系数是通过泊松比来体现的！

请教各位：我运行了n次gremms大侠的inp文件，错误提示：ERROR: SPECIFIED STANDARD_MEMORY VALUE OF 256 MB IS TOO SMALL TO RUN THE

ANALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE VALUE IS 269 MB. LOOK AT MEMORY ESTIMATES SECTION OF .dat FILE FOR FURTHER INFORMATION. 我已经改了虚拟内存，而且我的内存是512mb,不知道是什么原因 急盼各位大侠解答，急！

- 11 -

论坛汇总第一部分

你需要修改ABAQUS_V6.ENV中的内存设置,在WIN2000中ABAQUS支持2GB的物理内存，而在WINXP中支持3GB的物理内存，你的ABAQUS_V6.ENV中的默认只有256MB，你改为大于或等于512MB就可以了。goog luck！ 你在安装的\下面找到该文件后打开,在文件的开始就有,一看便知 现在把inp文件传上。

341818-BB8.inp (0.0k) 十、【寻找】弹塑性模型中*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC *ELASTIC,TYPE=ISOTROPIC 1.0E7, 0.30 *DRUCKER PRAGER 36, 1.0, 36

*DRUCKER PRAGER HARDENING 5E4 *DENSITY 2.26E3

参数取值如上,为何初始地应力不平衡,而把 *DRUCKER PRAGER HARDENING 5E4

的参数调大,8E5,地应力变平衡,为什么?地应力平衡与DRUCKER PRAGER HARDENING有关吗? 弹性模型也演算过, *ELASTIC,TYPE=ISOTROPIC 1.0E7, 0.30平衡,位移为1E-6以上,为什么?

GeoStatic是将模型恢复到初始的真实应力状态的一个平衡过程。你的模型在弹性阶段和塑性阶段对应的应该是同一个初始状态。你是想把达到塑性破坏的状态作为initial condition吗?

斑竹大人,你好,我的计算模型分别采用线弹性材料模型和弹塑性材料(DRUCKER),所加的初始应力相同 *INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC rock,0,0, -1.909E6,-86.2,0.5用弹性材料,初始位移很小(1E-6),地应力平衡,结果合理用弹塑性材料,初始位移较大(1E-1),地应力不平衡,结果不合理接着我就遇到了上面的问题,不知道我说明白了没有? 我的意思是要在GEOstatic内保证初始平衡,这个阶段是不用考虑材料的塑性啊,刚性啊,流动性啊的。因为你使用的是TYPE=STRESS

对呀,GEOstatic内保证初始平衡,这个阶段是不用考虑材料的塑性啊,刚性啊,流动性啊的。因为使用的是TYPE=STRESS但是后来的计算需要先前定义的模型材料参数材料属性分弹性和弹塑性同一模型,可以用不同材料参数,而其他部分一样当用弹性材料时,地应力很容易平衡,但是我改成弹塑性材料,地应力就不平衡了,是这个意思,不是同一过程的不同阶段而是两个分别的过程在GEOstatic内保证初始平衡我也并没有考虑材料的属性,而是前面定义的材料属性

我当然知道你已经是不同材料的不同计算case。但是平衡是使用的材料参数不就是一个density或指定应力吗?

谢谢斑竹！问题已解决，在用GEOSTATIC平衡时，我认真阅读了ＡＢＡＱＵＳ的例子，发现在地应力平衡与DRUCKER PRAGER HARDENING参数有一定的观系，当DRUCKER PRAGER HARDENING较大时，与弹性材料平衡一致，当DRUCKER PRAGER HARDENING较小时 ，需要添加新的边界条件，强制使其初始位移为０． 十一、求助:initial conditions的一个问题

不好意思，我看了以前的关于这方面的讨论，但是，还是没有解决自己的问题，实在没有办法，拿出来占大家的时间，谢谢！我做的类似桩基，是土体中加入桶形的桩基础。现在对土进行了iinitial conditions，现在土的初始位移基本为0了，然后我对桶同样用initial conditions，想使桶的初始位移也接近于0，但是，发现这个initial conditions对于桶不起作用，换了很多的参数，题目也不报错，好多天了，没有解决，希望得到大家的帮助，谢谢，谢谢！

315031-Job-1.inp (0.0k) *Material, name=SOIL *Density 1600.,

*Drucker Prager

- 12 -