修改灰坝命令流

finish
/clear
/filname,姚孟电厂
/title,灰坝
/units,si
!定义单元类型和材料属性
/prep7
!坝截面拖拉成实体
et,1,solid95
!坝体截面
et,2,plane82
!1号材料粉煤灰
mp,ex,1,3.43e10
mp,prxy,1,0.33
mp,dens,1,1240
!建模
k,1
k,2,280,0
k,3,105,50
k,4,100,50
k,5,82.5,55
k,6,76,55
k,7,58,46
k,8,114,46
k,9,84.925,52.5
k,10,49.75,62.55
k,11,43.25,62.55
k,12,18.125,52.5
!创建线面形成一级子坝截面
l,1,2
l,2,3
l,3,4
l,4,1
al,1,2,3,4
!创建线面形成二级子坝截面
l,7,8
l,8,5
l,5,6
l,6,7
al,5,6,7,8
!创建线面形成三级子坝截面
l,12,9
l,9,10
l,10,11
l,11,12
al,9,10,11,12
rectng,0,84.925,0,52.5
rectng,0,114,0,46
aplot
!对断面进行布尔操作
aovlap,all
/pnum,line,1
/pnum,area,1
nummrg,all
numcmp,all
allsel
aadd,all
nummrg,all
numcmp,all
aplot
!对面进行网格划分
allsel
type,2
mat,1
asel,s,,,1
mshkey,0
mshape,0,2d
amesh,all
eplot
!拉伸成坝体单元，采用solid95单元和1号材料
!纵向拉伸80份
allsel
extopt,esize,80,0,
type,1
mat,1
vext,all,,,,,400
/view,1,1,1,1
/angle,1
vplot
eplot
!施加边界条件
/solu
csys,0
dsys,0
allsel
nsel,s,loc,y,0
nplot
d,all,all
allsel
gplot
!施加外载。外载包括，重力，净水压力，由于下部积压的是沉淀粉煤灰，不考虑扬压力（渗透压力和上举力）
allsel
acel,0,9.8,0
!选择挡水坝面，施加水平静水压力
allsel
asel,s,loc,x,0
esla,s,,,all
!设定将显示压力的方向
/psf,pres,norm,2,0,1
!给定荷载梯度
sfgrad,pres,0,y,0,-2200
sfe,all,4,pres,,2200*60
/replot
save,dam_static_solu,db
!静力求解
allsel
outres,all,all
autots,on
nsubst,20
solve
save,dam_static_solu,db
finish

不错点资料哈，谢谢分享