大坝溢流坝段坝踵处应力大于坝趾处？请水电专业的人帮忙检查命令流~

suking_love

我分析的是溢流坝段坝体内的应力分布。
查看等效应力时，发现坝踵处的应力总大于坝趾处，可是一般都是坝趾处较大~
不好意思，我还不能上传图片。

现附上命令流，请热心的人帮忙检查检查~
我考虑的外荷载有，重力、扬压力、静水压力

运行完之后，需要进行荷载步的相减才能查看结果~

PS：荷载步的相减是为了消除重力引起的初始位移！

fini
/cle
/units,si
h=26.5             !定义相关变量
h1=31.94
h2=19.22
h3=11
gm=1000
/filename,lisuqin
/prep7             !定义单元材料及材料属性
et,1,plane42       !外围混凝土单元
mp,ex,1,2.55e10
mp,nuxy,1,0.167    !假定泊松比
mp,dense,1,2400
keyopt,1,3,0
keyopt,1,5,0
keyopt,1,6,0
et,2,plane42       !浆砌石单元
mp,ex,2,1.15e10
mp,nuxy,2,0.16     !假定泊松比比混凝土小
mp,dense,2,2300
keyopt,2,3,0
keyopt,2,5,0
keyopt,2,6,0
et,3,plane42        !基岩单元
mp,ex,3,4e10
mp,nuxy,3,0.167
mp,dense,3,2660
keyopt,3,3,0
keyopt,3,5,0
keyopt,3,6,0
!定义关键点
!坝剖面
K,1,0,0
K,2,5,0
K,3,6,1
K,4,27.6,1
K,5,27.6,9.3
!材料分线
K,10,6,1.5
K,11,27.1,1.5
K,12,27.1,6.2
K,13,17.58198,6.2
K,14,3.38198,25
K,15,0.6,25
K,16,0.6,8
K,17,6,8
K,18,8.51815,18.2
K,19,0.6,18.2
K,20,0,25
K,21,27.6,6.2
!溢流面
!圆弧段
K,6,27.1,9.3
k,7,16.5011,10.16
!中间直线段
K,8,7.87419,21.66615
!顶部曲线段
K,9,0,25.98761
K,30,1.16549,26.5
K,31,5.16549,24.64157
!灌浆廊道
!圆心
K,22,2.5,3
K,23,2.5,5.25
K,24,2,5.25
K,25,5.5,5.25
K,26,5,5.25
K,27,5,3
K,28,3.75,5.25
k,29,22.1811,7.32
circle,28,1.25,,,180
circle,28,1.75,,,180
larc,6,7,29
larc,30,31,14,6.259
larc,9,30,15,1.32
nummrg,all
!地基
K,32,-0.6799,0.97128
K,37,-40.4299,0.97128
K,34,-40.4299,-52.02872
K,35,80.6,-52.02872
K,38,80.6,1
!定义面
a,15,14,18,19
a,19,18,13,12,11,10,17,16
a,9,20,15,14,18,13,12,21,5,6,7,8,31,30
a,20,15,19,16,17,10,11,12,21,4,3,2,1
a,22,27,26,25,36,24,23
asba,4,5
a,22,27,26,25,36,24,23
a,23,33,26,27,22
asba,4,5
a,37,34,35,38,4,3,2,1,32
wprota,,-90
wpoff,6,,8
asbw,2
wpoff,11.51815,,-1.8
asbw,8
wpoff,9.58185,,-4.7
asbw,6
wpoff,-21.1,,-0.5
asbw,10
wpoff,-6,,7
asbw,11
wpcsys,-1

esize,1                  !@浆砌石划分网格
mshape,0,2d
mshkey,1
allsel
asel,s,,,5
amap,5,18,19,16,39
amesh,5
allsel
asel,s,,,2
amap,2,39,17,40,13
amesh,2
allsel
asel,s,,,9
amap,9,40,10,11,12
amesh,9
esize,1                  !外围混凝土划分网格
mshape,0,2d
mshkey,1
asel,s,,,8
amesh,8
asel,s,,,12
amap,12,42,43,41,4
amesh,12
asel,s,,,6
amap,6,43,1,2,3
amesh,6
asel,s,,,10
amap,10,44,16,15,20
amesh,10
esize,1
mshape,0,2d
mshkey,0
asel,s,,,13
amesh,13
esize,2                !基岩划分网格
mshape,0,2d
mshkey,1
asel,s,,,4
amap,4,37,34,35,38
amesh,4
esize,0.1               !现浇拱划分网格
mshape,0,2d
mshkey,1
asel,s,,,7
amap,7,24,23,26,25
amesh,7
esize,0.5                 !内部C10划分网格
mshape,0,2d
mshkey,1
asel,s,,,1
amesh,1
esize,1                 !溢流面划分网格
mshape,0,2d
mshkey,0
asel,s,,,3
amesh,3


nsel,all
nsel,s,loc,x,-40.4299
nsel,r,loc,y,-52.02872,0.97128
d,all,ux,0
nsel,all
nsel,s,loc,x,80.6
nsel,r,loc,y,-52.02872,0.97128
d,all,ux,0
nsel,all
nsel,s,loc,x,-40.4299,80.6
nsel,r,loc,y,-52.02872
d,all,all,0
allsel


finish
/solu
save
time,1
allsel
acel,,9.8
allsel,all
iswrite,on
solve


/solu
time,2
isfiled,read,lisuqin,ist,,2
allsel

nsel,all
nsel,s,loc,x,0
sfgrad,pres,,y,0,-gm*9.8
sf,all,pres,gm*9.8*(h1+1)

nsel,all
nsel,s,loc,x,27.6
nsel,r,loc,y,1,9.3
sfgrad,pres,,y,1,-gm*9.8
sf,all,pres,gm*9.8*h2

nsel,all
nsel,s,loc,y,9.3
sfgrad,pres,,x,27.6,1
sf,all,pres,gm*9.8*(h2-9.3)

nsel,all
nsel,s,loc,x,11.7611,16.5011
nsel,r,loc,y,10.16,19.22
sfgrad,pres,,y,16.48,-gm*9.8
sf,all,pres,

nsel,all
nsel,s,loc,x,16.5011,22.1811
nsel,r,loc,y,7.32,10.16
sfgrad,pres,,y,10.16,-gm*9.8
sf,all,pres,gm*9.8*(h2-10.16)

nsel,all
nsel,s,loc,x,22.1811,27.1
nsel,r,loc,y,7.32,9.3
sfgrad,pres,,y,9.3,-gm*9.8
sf,all,pres,gm*9.8*(h2-9.3)

esel,all
esel,s,,,411,416,1
sfgrad,pres,,x,0,-21053
sfe,all,3,pres,0,313012
esel,all
sfe,411,4,pres,0,0.24*gm*9.8*(h1-h2)+9.8*gm*h2-1295.13*2
esel,all
sfe,405,4,pres,0,0.24*gm*9.8*(h1-h2)+9.8*gm*h2-1295.13*3
esel,all
esel,s,,,381,404,1
sfgrad,pres,,x,0,-1295.13
sfe,all,2,pres,0,0.24*gm*9.8*(h1-h2)+9.8*gm*h2-1295.13*4


allsel
solve