文档介绍：1我做地基中桩、土、承台相互作用。先算得重力作用下每个单元中心的六个应力,然后作为初始应力加到inp文件中,但我的inp文件是在CAE操作后得到的,所以part的单元结点编号以及各个module的特征在此inp文件中很明显,请问应该把initialconditions,type=stress这个keywords加在哪个位置合适呢?老是提示我关键词位置错误,那位大侠碰到过这个问题或者知道怎么解决,请告诉小弟啊,比较急。谢谢啦~~~~~~~~~insertbefore*STEPANSYS2【讨论】*initialconditions,type=stress,geostatic选项的功能Forpersonaluseonlyinstudyandresearch;mercialuse各位高人请指教:    在土工计算中,初始地应力是通过给出一点的竖向应力和侧压力系数来计算的,并假定竖向应力按静水压力分布。从ABAQUS帮助文件中知,ABAQUS中的*initialconditions,type=stress,geostatic选项可以实现这一功能,然后用*Geostatic求解程序来协调初始应力和外荷载之间的平衡,在这里初始应力主要是重力。这样求解以后,理想状态下位移应该为0。但我在应用中用*Geostatic求解后,位移总是很大,并且应力分布也不符合规律,不知为什么?    初始应力的关键字和语句行如下(平面应力)*Initialconditions,type=stress,geostatic  strip1,,8,-,2,:单元集,竖向应力值1,与竖向应力值1相应的高程,竖向应力值2,与竖向应力值2相应的高程,侧压力系数    重力荷载选项如下*Dload_PickedSet60,GRAV,,0.,-1.    请问以上语句操作和语句有什么错误?附件里有我做的一个简单模型。    多多指教!68948-()我看你的nodeset和elementset都用相同的名字,这样对操作不好的。同时只要对节点操作就可以,如果名字相同的话,就表示节点和单元都同时操作了。请检查你的边界条件。的确使用平面应变单元得出的结果不一样、但是在GEOSTATIC的回代计算中、我发现土的侧压力系数的设定影响到最后的变形极大。不知道seansheu大虾有没有办法把我附件中的例子、调整到GEOSTATIC完成后变形结果接近为0。73087-()你这个是用的SI单位制吗?若不是,而长度单位用mm吗?请予以说明。1 个人觉得你的网格有点小问题,四边形单元对网格形状敏感,而你用femap产生的网格质量并不高2 geostaticstep不用unsymm,也不用nlgeom,只要一个增量步即可完成,本身不存在非线性过程。3 我对你的数据文件稍微修改,我做得结果如图。单位是SI单位制我修改了你的模型将P3和p5合成了一层,inp文件如下。73547-()这个问题我想请教以下:应用yzf1111大侠的inp计算的结果:地表的沉降量是很小,可是地表和地底的竖向正应力和*initialconditions....里的值不符。不知什么原因?能否改为一致?非常感谢yzf1111师兄的热情帮助、谢谢你辛勤的工作、我在计算中的单位是SI。我发现你使用HyperMesh重新生成了一个正规四边形网格的模型。这对节理的考虑可能不确切。既然GhostFP阁下没有时间了、我再上传另外一个关于GeoSTATIC施加方法的INP(部分)、同样能达到初始应力和变形NearZero的水平。80956-()不同意楼上。nodeset和elementset当然可以用一样的名字。其实,初始应力的施加还不是那么简单,比如说有接触问题时,或者有的结构部分不需要加初始应力。好像ABAQUS有加接触应力,不知上面的兄弟有没有用过。对于不需要加初始应力的结构部分,如何考虑它的刚度影响,不知有没有人考虑过你的elementtype有問題,雖然土工我不太熟,不過土工計算應該用planestrainelementtype,應該用CPE8R而不是CPS8R,我改過後,得到的結果與你假設相符...變形圖但是在3D问题中,地压平衡这一步与土的弹性模量有很大关系,比如对天津这边来说,用5MPa,地压平衡基本是建立不起来的。但是若改成500MPa,就可以建立起来,不知道大家注意到没有,abqus中给的例子土的弹性模量都是在100MPa---1000MPa范围内的,所以这个问题还请大伙进一步讨论,Drcafe,你们既然用的是正版,可以向abaqus工程师咨询一下的。在下