生成地应力错误求助，图片和命令流上传在附件

acufo

生成地应力错误求助，图片和命令流上传在附件.土体中设置了一条沟，土和沟的材料参数与单元模型是一样的，求解后只有沟的一段有沉降，找不到错误的原因。命令流在附件可以查看

acufo

知道了，2个不同的面要一起拉伸才行

acufo

修改后的命令流：finish               !隐示建模不移除自重位移大变形/clear
/filname,cant35
/prep7
k,1,0,0,0
k,2,5,5,0
k,3,10,5,0
k,4,15,0,0
k,5,15,10
k,6,10,10
k,7,5,10
k,8,0,10
a,1,4,5,8
l,7,2
l,6,3
l,2,3
asbl,1,all
kgen,2,1,,,,,-20,,,
l,1,9
lesize,1,,,15
lesize,4,,,10
lesize,2,,,10
lesize,10,,,5
lesize,8,,,5
lesize,9,,,5
lesize,5,,,5
lesize,6,,,5
lesize,7,,,5
et,1,mesh200
et,2,mesh200
et,3,solid185
keyopt,1,1,6
keyopt,2,1,10
mp,dens,1,2200
mp,ex,1,0.5e9
mp,nuxy,1,0.3
mp,dens,2,2200
mp,ex,2,0.5e9
mp,nuxy,2,0.3
asel,all
asel,s,area,,2
aatt,1,,1
mshape,0
mshkey,2
amesh,all
asel,all
asel,s,area,,3
aatt,2,,1
mshape,0
mshkey,2
amesh,all
asel,all
asel,s,area,,2  
asel,a,area,,3  !2个面一起拉伸
type,3                 !拉伸成体的单元类型         
mat,1                 !拉伸成体的材料参数
extopt,esize,20        !拉伸成体路径方向的网格数目
vdrag,all,,,,,,3   !拉伸路径
nsel,all               
nsel,s,loc,x,0
nsel,a,loc,x,15
d,all,ux             !x约束
nsel,all
nsel,s,loc,y,0
d,all,uy          !y约束
nsel,all
nsel,s,loc,z,0
nsel,a,loc,z,-20
d,all,uz            !z约束
allsel,all
acel,,9.8
/solu
solve
/post1
prnsol,s,comp       !节点应力分量
plnsol,s,x
plnsol,u,y