★地应力平衡图解

*initialconditions,type=stress,input=0.csvmdb.models['Model-1'].setValues(noPartsInputFile=ON)（将网格划分精密一些，地应力平衡精度更高一些，2010-7-12注）注意：导入的inp、csv或其他文件里不得有空行，否则会出现element0。允许有空格。--------------------------------------------------------------------方法一：桩土地应力平衡不容易通过，可先直接指定应力：*initialconditions,type=stress,geostaticSet-pile-soil,0,0,-510000,30,0.6,0.6将桩土摩擦系数设为0，容易平衡，输出含有应力s11、s22等的应力文件；导入应力文件，将摩擦系数提高为一个较小的数值（低于正常的摩擦系数），再计算，输出应力文件；再提高摩擦系数，导入应力文件，计算，重复上述操作，直到摩擦系数达到正常值。上述如果不行，可以干脆将桩与的接触改为tie，2010-8-18的模型按照上面的步骤操作不行，后来tie就可以了，且收敛的速度很快。桩土之间用摩擦容易出现“Theratioofdeformationspeedtowavespeedexceeds1.0000inatleastoneelement.”的错误。而用tie，对于动力学更容易收敛，一般不会出现上述错误提示（2010-8-19的锤击桩模型将土分为内外两层，内层网格密度与桩划分的密）。方法二：桩土地应力平衡方法：（建议此法）1、按照实际工程要求建立桩土模型model-1；2、模型树上的model-1上面copy，复制成另外一个模型model-2；3、把mode-2里的所有与土无关的东西全部删除，包括part、assembly、set、interaction、load等；4、在load模块的BC里约束原来土与桩接触的面，约束住土不要进入桩孔内（竖向不要约束）；5、注意：土不能受到重力方向的约束，包括土四周的边界，让土自由下沉。6、不用导入初始应力，提交model-2的job，产生应力文件-2，编辑成ABAQUS要求的格式；6、将应力文件-2导入model-1，提交计算model-1的job，产生应力文件-1，编辑成ABAQUS要求的格式，重新导入model-1，提交，平衡完毕；7、删除辅助用的model-2及其产生的文件。该法通过率高，精度高（2010-8-19用锤击桩模型试算，用此法平衡达到10e-4，而用其它方法平衡的精度达10e-3）。问题：2010-8-19的模型（锤击桩内外土tie的模型）用该法没有平衡成功，还是用方法一平衡成功。地应力平衡方法：

第一步：建立模型，材料，分析步（GEOSTATIC）

第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8Y

第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

(请严格按照这个格式，注意大小写的字母)，例如：mdb.models['Model-1'].setValues(noPartsInputFile=ON)第四步：提交J0B，完成后

第五步：按以下步骤,Roport---ReportFieldOutput---选中S11,S22,S33,S12,S13,S23---Name:XX.INP

---Write中选择FieldOutput-------------ok!!!第六步：用软件(推荐使用UltraEdit很好编辑的)打开XX.INP[注：先用UltraEdit打开，删除中间的英文，再用excel打开]，保存格式内容

单元号

S11

S22

S33

S12

S13

S23（请注意，在保存内容中没有这一行的）

1

,

.

,

.

,

.

,

.

,

.

,

.

2

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有

上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2.

.......这些数字。例如：*elprints,*initialconditions,type=stress,geostaticSet-all,0,0,-510000,30,0.6,0.6*initialconditions,type=stress,geostaticSet-pile,0,0,-250000,10,0.2,0.2Set-soil,0,0,-255000,15,0.6,0.6一个part分层土地应力平衡SHAPEPart-1是个部件，要在assembly组装成实体instance，instance的名字一般是在部件的名字后面加“-1”，例如部件part-1形成实体part-1-1。两个part分层土地应力平衡SHAPE另外一种方法1：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

(请严格按照这个格式，注意大小写的字母)2：提交J0B，完成后3：按以下步骤,Roport---ReportFieldOutput---选中S11,S22,S33,S12,S13,S23---Name:XX.INP---Write中选择FieldOutput-------------ok!!!以上3个步骤可以代替以下命令：*elprints,手工输入地应力平衡SHAPE*initialconditions,type=stress,geostaticSet-pile,0,30,-500000,10,0.2,0.2set-soil,0,30,-540000,0,0.6,0.6对于很厚的土层必须先进行地应力平衡，否则有些单元会出现过大的扭曲变形，最后运算不了，出现ERROR:FIXEDTIMEINCREMENTISTOOLARGE（即使在step中打开非线性开关也不行），必须先手工输入地应力平衡初始条件，见上面。在inp中加入*elprints,将各单元的应力导出（*.csv），再*initialconditions,type=stress,input=*.csv导入初试应力（程序自动将取排在后面的初始应力为准，自动把第一个初始应力除去如下图，*initialconditions,type=stress,geostatic和下面一行自动被除去），这样迭代平衡效果非常好。对于采用c3d8等非缩减积分单元，会出现如下图的应力文件：用下面的方法就OK：（不要选择积分点，而是质心）由于老板项目多，自己论文忙，好长一段时间都没来SIMWE了，一上来就收到不少朋友的消息，问我的地应力平衡方法，现在我把整个过程详细说明一下，希望对大家有所帮助。

地应力平衡方法：

第一步：建立模型，材料，分析步（GEOSTATIC）

第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8

第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

(请严格按照

这个格式，注意大小写的字母)

第四步：提交J0B，完成后

第五步：按以下步骤,Roport---ReportFieldOutput---选中S11,S22,S33,S12,S13,S23---Name:XX.INP

---Write中选择FieldOutput-------------ok!!!

第六步：用软件(推荐使用UltraEdit很好编辑的)打开XX.INP，保存格式内容

单元号

S11

S22

S33

S12

S13

S23（请注意，在保存内容中没有这一行的）

1

,

.

,

.

,

.

,

.

,

.

,

.

2

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有

上面那个“单元号”这一行，而且单元号前面也没有什么PART名字什么的，就是1，2.

.......这些数字。

第七步：在ABAQUS----Model---Editkeywords---Model-1(这就是你的Model名字)---在材料属性后面

加上：*initialconditions,type=stress,input=xx.inp

完成

第八步：重新提交JOB，OK

第九步：如果你还没有成功的话，那我只能说----------------我无语了。：）

ABAQUS的这项功能确实很不错。：）这个功能让基坑开挖、隧道开挖等的初始应力，开挖后的残余应力很好的显示；也可以很好的模拟铁路设计中的工后沉降的概念，在地应力平衡后，加上荷载所得沉降即为工后沉降；也很好的模拟了桩土复合地基的问题，如果没有初始应力的模拟，使土对桩产生了挤压应力，从而通过设定摩擦系数就可以模拟了桩与土之间的摩擦力；除此之外，在进行挡土墙计算时也需要ABAQUS的这项功能，反正很多都用得着。

希望我的建议能帮助初学者（其实我也是初学者，呵呵）很快掌握这个技巧。原文地址：\o"有关ABAQUS数据的后处理"有关ABAQUS数据的后处理作者：\o"举举"举举三种方法：1．

XYDATA

通过Tools-XYdata-edit,thencopytoexcelormatlab,etc.2．

通过Plugins-Tools-excelutilities直接在Excel中输出数据和图形3．

通过Report-XY，将结果输出至abaqus.rpt，此文件后缀名可改为.txt，通过MATLAB读取。桩土共同作用对于新手来说一直是个难题，本人结合3年做桩土分析的经验，将桩土分析中遇到的问题做个总结：

主要思想是在计算土体自重时将桩杀死，然后约束土体中桩孔的横向自由度，建立step1，通过静力分析步计算土体的自重应力，然后将土体的自重应力编辑后，通过initialconditions，type=stress，input=自重应力文件（可以是txt或是inp文件）导入到原模型中，然后添加step2分析步，在step2分析步中，激活桩，同时建立接触，施加桩的自重及其他荷载，即可完成桩土的共同作用分析。具体操作过程如下：

1.首先建立模型1，如model1，在该模型里建立桩和土的几何模型，并赋予材料属性，将所有的桩设置为一个集合，并命名（如z1），建立step1，使用静力分析步，注意，在editkeyword中在step1定义的后面添加两行语句，第一行为*modelchange，remove，第二行为z1（即前面定义的桩的集合），然后在荷载里定义土体的重力加速度，约束土体的位移，同时约束土体中桩孔的横向位移（如，z为重力加速度方向，则约束x和y的自由度），计算分析，得到土体的自重应力，编辑自重应力文件，命名为XXX.txt。

2.进入到模型1中，将模型1另存为模型2，此时，模型2和模型1是相同的，在模型2中，添加step2，编辑editkeyword，在边界条件定义前，添加*initialconditions，type=stress，input=XXX.txt，然后在step2定义后面添加两行语句，第一行为*modelchange，add第二行为z1（桩的集合名），然后在step2中建立桩和土的接触，施加桩的自重，然后在step中杀死土体中桩孔的横向自由度（对应step1中桩孔的横向自由度约束），后面的荷载可以在step2中施加，也可以建立新的step施加。

3.补充：建立接触的分析步（本例中为step2，初始步要设置的小的，本人一般设置为0.01，这样有利于接触的顺利建立）；当模型较大时，地应力文件的编辑工作量很大，本人找朋友编了个小程序，可以将地应力文件自动生成abaqus要求的格式，该程序将在附件中给出。附件太大，无法上传，悲剧[local]1[/local]\o"汇总"汇总,

\o"经验"经验查看3D实体截面的内力，如截面合力、截面总弯矩，或者截面各方向的分力、分弯矩，按照aba-aba大侠编写的《ABAQUS常见问题汇总》，可以设置cutsurface来实现。但这种方法较难，而且对每一个需要查看内力的截面都要单独设置，但要查看内力的截面较多时，就更麻烦。通过odb文件可以较为简单的获取截面内力。具体操作步骤如下：

a)

打开odb文件，在窗口中只保留要获取内力的部件，如某根柱、梁等；若保留整个模型，则获取的截面内力是整个模型在该截面的内力；