- 积分
- 0
- 注册时间
- 2009-4-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-8-29 16:12:50
|
显示全部楼层
来自 上海
本帖最后由 fifasong 于 2011-8-29 16:19 编辑
我跟楼主的做法部分类似。第一步重力计算,将桩的弹性模量调小,大致与土体一个数量级(或者用几层土的平均值);第二步应力平衡;第三步加载;
初始应力文件读取可用pyhon编个简单的程序,如下:
#export ELEMENT CENTROID STRESS COMPONENT
from odbAccess import *
from abaqusC**tants import *
odb=openOdb(path='xx-grav.odb')
myfile=open('xx-grav.txt','w')
partinstance=odb.rootAssembly.instances['SOIL-1']
step='Step-1'
frame=odb.steps[step].frames[-1]
stressField=frame.fieldOutputs['S']
field=stressField.getSubset(region=partinstance,position=CENTROID)
fieldValues = field.values
for v in fieldValues:
print >>myfile,"SOIL-1.%s,%20.6e,%20.6e,%20.6e,%20.6e,%20.6e,%20.6e" %(v.elementLabel,v.data[0],v.data[1],v.data[2],v.data[3],v.data[4],v.data[5])
myfile.close()
odb.close()
若有其他instance,将上面代码复制一下,把SOIL-1改为相应的instance即可。
|
|