完整的地应力平衡方法总结计划.docxVIP

地应力平衡方法： 第一步： 建立模型，材料，分析步（ GEOSTATIC） 第二步： 施加荷载， LOAD，选择施加重力 GRAVITY，在你想施加重力的方向输入数值 9.8 第三步： 在命令行中输入 mdb.models[模型名字 ].setValues(noPartsInputFile=ON) (请严格按照 这个格式，注意大小写的字母 ) 第四步： 提交 J0B，完成后 第五步： 按以下步骤 ,Roport---ReportField Output---选中 S11,S22,S33,S12,S13,S23--- Name: XX.INP ---Write 中选择 Field Output-------------ok!!! 第六步： 用软件 (推荐使用 UltraEdit 很好编辑的 )打开 XX.INP，保存格式内容 单元号 S11S22S33S12S13S23（请注意，在保存内容中没有这一行的） 1,.,.,.,.,.,. 2,.,.,.,.,.,. .,.,.,.,.,.,. 1 / 4 .,.,.,.,.,.,. 这个结果文件是最重要的，在所保存的文件中只有数值部分，没有英文字母，没有 上面那个 “单元号 ”这一行，而且单元号前面也没有什么 PART名字什么的，就是 1，2. .......这些数字。 第七步： 在 ABAQUS----Model---Edit keywords---Model-1(这就是你的 Model 名字 )---在 材料属性后面 加上： *initial conditions,type=stress,input=xx.inp 完成 第八步： 重新提交 JOB，OK 第九步： 如果你还没有成功的话，那我只能说 ----------------我无语了。：） ABAQUS的 这项功能确实很不错。： ）这个功能让基坑开挖、隧道开挖等的初始应力，开挖后的残余应力很好的显示；也可以很好的模拟铁路设计中的工后沉降的概念，在地应力平衡后，加上荷载所得沉降即为工后沉降；也很好的模拟了桩土复合地基的问题，如果没有初始应力的模拟，使土对桩产生了挤压应力，从而通过设定摩擦系数就可 以模拟了桩与土之间的摩擦力；除此之外，在进行挡土墙计算时也需要 ABAQUS 的这项功能，反正很多都用得着。 希望我的建议能帮助初学者（其实我也是初学者，呵呵）很快掌握这个技巧。初应力施加方法介绍（经过计算验证的） 2 / 4 一种方法： 1、模型建好 2、用 *INITIALCONDITIONS,TYPE=STRESS,GEOSTATIC命令施加计算得出的地应力，一般沿高度的线性分布，手工计算； 3、第一个 step 选用 geostatic，中用 *DLOAD 命令施加重力荷载，施加计算初始应力时用到的所有外载， geostatic 步会自动计算分析， 看你施加的初始应力与该步计算所得的应力是否平衡，如存有差异，会自动平衡该应力差，因此会产生沉降，以消除该差异。 如果沉降不大，就认为初应力施加是理想的。 4、继续下面的分析步。 第二种方法： 1、建立模型，一般采用线弹性材料，施加边界条件，不要把所有节点都约束，要跟正常计算一样，施加重力，如有 预加荷载，同样施加。 2、建立一个 static 分析步，即可获得模型中的应力分布，利用菜单 Report Field Output 项，在 Variable 对话框中的 应力分量（ 3－d）， 2－d 只有四个分量，存为 filenme.txt 文件 3、可用 excel 打开该文本文件， 2－d 要把不需要的应力值全部添上，为 0，然后存为 filenme.csv 文件，后改为 filenme.inp 文件， 如你的 excel 不能进行这种转换，可利用 fortran 程序转换。附后。 4、利用 initial conditions,type=stress,INPUT=filenme.inp 即可读入初始应力，abaqus自动将该应力分配给高斯积分点。 3 / 4 5、再利用该方法时，要修改 cae 形成的计算模型的 inp 文件，去掉所有instance,assembly，internal 等关键词，使 inp 文件等同于 我们手写的命令流格式。 6、计算时同样要建立 geostatic 步，否则在计算结果中可能未考虑初应力。 7、至于论坛中限制所有节点的方法，对重力荷载有效，如果在初始状态时就有预加荷载就不可行了。 4 / 4