h=180				!Ӹ߶
fini_days=30			!ɴӽܵ
ht=h/fini_days			!ÿɵĴӸ߶
h1=100				!ˮλ߶
h2=80				!ˮλ߶
h3=100				!岻ͬϵķֽߣΪ100m
gm=1000				!gmΪˮ
fia1=90-atan(0.75)*180/3.1415926	!бн
fia2=90-fia1			!תǶ
/prep7
et,1,plane42			!֮
et,2,solid65			!Ԫ
et,3,solid45			!ʯԪ	
keyopt,2,5,2			!solid65ԪĽڵӦ

!1Ų
mp,ex,1,2.85e10			!100m»Ĳ 
MP,PRXY,1,0.167 
TB,CONC,1,1,9, 
TBDATA,,0.3,1,1.96e6,22e6 	!ſдϵΪ0.3ǿfc=1.96e6pa,ѹǿfs=22e6pa,
mp,dens,1,2400

!2Ų
mp,ex,2,2.6e10			!100mϻĲ 
MP,PRXY,2,0.167 
TB,CONC,2,1,9, 
TBDATA,,0.3,1,1.2e6,17.5e6 	!ſдϵΪ0.3ǿfc=1.2e6pa,ѹǿfs=17.5e6pa,
mp,dens,2,2400

!3Ų
mp,ex,3,2.9e10
mp,PRXY,3,0.3
mp,dens,3,2600

!ģ
!ؼ
k,1
k,2,0.9*0.75*h+0.1*h
k,3,0.1*h,0.9*h
k,4,0,0.9*h
!
l,1,2
l,2,3
l,3,4
l,4,1
!
al,1,2,3,4
RECTNG,0,0.1*h,0.9*h,h				!Ӷ
RECTNG,-1.5*h,0,-2*h,0				!ӻ
RECTNG,0,0.9*0.75*h+0.1*h,-2*h,0		!ӻ·
RECTNG,0.9*0.75*h+0.1*h,(2+0.9*0.75)*h,-2*h,0	!ӻ
RECTNG,-1.5*h,(2+0.9*0.75)*h,-2*h,h		!ܾ
AOVLAP,all					!浥Ԫճ
/pnum,line,1					!棬ߵĺ뿪
/pnum,area,1
nummrg,all					!ϲظԪ
numcmp,all
aplot

!3,Ķ
lsel,s,,,3,5,2
lesize,all,,,5
lsel,s,,,12,13,1
lesize,all,,,2
amesh,3

!1񣬰²
lsel,s,,,2,4,2
lesize,all,,,18
lsel,s,,,1
lesize,all,,,5
amesh,1
eplot

!2񣬰·
lsel,s,,,11
lesize,all,,,5
lsel,s,,,9,10,1
lesize,all,,,8,4
amesh,2

!4
lsel,s,,,14
lesize,all,,,5,4
lsel,s,,,7
lesize,all,,,8,4
lsel,s,,,6
lesize,all,,,5,0.25
amesh,4
eplot

!5
lsel,s,,,15,16,1
lesize,all,,,8,4
lsel,s,,,8
lesize,all,,,8,0.25
amesh,5
eplot

!4ߺ13ߣ7
lsel,s,,,4,13,9
LCCAT,all
lsel,s,,,19
lesize,all,,,5,4
lsel,s,,,20,21,1
lesize,all,,,20
amesh,7
eplot

!2ߺ12ߣ6
lsel,s,,,2,12,10
lccat,all
lsel,s,,,17
lesize,all,,,8,4
lsel,s,,,18,22,4
lesize,all,,,20
amesh,6
eplot

!ɾǰӵԪ
allsel
lsel,r,lcca
ldele,all

!ɰ嵥Ԫsolid65Ԫ1Ųϣıh3ϰĲ
EXTOPT,ESIZE,8,0,   			!ķ     
type,2
mat,1
VEXT,1,3,2,,,-0.75*h			!1ź3һģ
/view,1,1,1,1         
vplot
eplot

!ײԪԪsolid45ŵԪ3Ų
EXTOPT,ESIZE,8,0,        
type,3
mat,3
VEXT,2,,,,,-0.75*h
VEXT,4,5,1,,,-0.75*h
eplot

!ɴʯԪ
allsel
EXTOPT,ESIZE,5,4,        
EXTOPT,ACLEAR,1				!ɾ浥Ԫ
type,3
mat,3
VEXT,1,7,1,,,h		
/pnum,mat,1
eplot

!öԳƷԪģ
local,11,0,,,-0.75*h			!ֲϵ
csys,11					!ֲϵ
/PSYMB,CS,1				!ʾͬľֲϵ
DSYS,11					!ʾֲϵ
NSYM,Z,20000,all
ENSYM,30000, ,20000,all
nummrg,all			!ϲظڵ͵Ԫ
numcmp,all
eplot

!ı
esel,s,mat,,1			!ѡϱΪ1ĵԪ
eplot
nsle,s				!ѡԪϵĽڵ
nplot
nsel,r,loc,y,h3+1,h+1		!ѡh3ϵĵԪԱת
nsel,r,loc,z,-0.75*h+0.1,0.75*h-0.1
nplot
esln,s
eplot
MPCHG,2,all,			!ıΪ2Ų
/REPLOT 
allsel
eplot
save,dam_model,db		!ģļ
finish

/solu
!ʩӱ߽
csys,0
dsys,0
nsel,s,loc,x,(2+0.9*0.75)*h
nsel,a,loc,x,-1.5*h
nplot
d,all,ux
allsel

nsel,s,loc,z,h
nsel,a,loc,z,-(1+1.5)*h
nplot
d,all,uz
allsel

nsel,s,loc,y,-2*h
nplot
d,all,uy
allsel
gplot

!ʩ
ACEL,0,9.8,0,

!ѡΰ浥Ԫʩˮƽˮѹ
ESEL,S,TYPE,,2
eplot
nsel,s,loc,x,0
nsel,r,loc,z,-1.5*h+0.1,-0.1
nsel,r,loc,y,0.1,h1-5		!(h1-0.1)ѡĵԪ,
				!ʹʩӵֲسָ
esln,s
nplot
eplot
/PSF,PRES,NORM,2,0,1		!趨ʾѹķ
sfgrad,pres,0,y,0,-gm		!غݶΪ-gm
sfe,all,2,pres,,gm*h1
/replot

!ΰ浥Ԫʩˮƽˮѹ
allsel
ESEL,S,TYPE,,2
nsle,s
eplot
local,12,0,0.9*0.75*h+0.1*h,,,fia2
csys,12
/PSYMB,CS,1
DSYS,12

nsel,s,loc,y,0.1,h2/sin(0.75)-25	!(h2/sin(0.75)-1)ѡĵԪ,
					!ʹʩӵֲسָ
nsel,r,loc,z,-1.5*h+0.1,-0.1
nsel,u,loc,x,-1000,-2
esln,s
nplot
eplot
/PSF,PRES,NORM,2,0,1		!趨ʾѹķ
sfgrad,pres,0,y,0,-gm		!غݶΪ-gm
sfe,all,4,pres,,gm*h2
/replot

!ʩˮλh1͸ѹ
!ͬʱʩˮλh2ĸ
!ѡһְӵ׵Ԫ
csys,0
dsys,0
/psymb,cs,1
allsel
nsel,s,loc,y,0
nsel,r,loc,z,-135+0.1,0.1
esln,s
esel,r,type,,2
eplot
SFCUM,PRES,ADD			!غǵӵ
sfe,all,5,pres,,gm*h2		!ˮλĸ
p0=gm*h1/(0.9*0.75*h+0.1*h)
sfgrad,pres,0,x,0,-p0
sfe,all,5,pres,,gm*h1
/replot
!ѡڶְӵ׵Ԫ
allsel
nsel,s,loc,y,0
nsel,r,loc,z,-270+0.1,-135-0.1
esln,s
esel,r,type,,2
eplot
p0=0
sfgrad,pres,0,x,0,-p0		!бΪ0
sfe,all,3,pres,,gm*h2		!ˮλĸ
p0=gm*h1/(0.9*0.75*h+0.1*h)
sfgrad,pres,0,x,0,-p0
sfe,all,3,pres,,gm*h1
SFCUM,PRES,repl			!غ

!ӵ׵غǷȷ
nsel,s,loc,y,0
esln,s
esel,r,type,,2
/replot
save,dam_static_solu,db
allsel
outres,all,all
autots,on
nsubst,20
solve
save,dam_static_rst,db

/post1
esel,s,type,,2

set,1,last			!ȡһغɲӲĽ
PLNSOL,U,X
PLNSOL,U,Y
PLNSOL,EPTO,1,0,1		!һӦ
PLNSOL,S,1,0,1			!һӦ
PLNSOL,S,3,0,1			!Ӧ
/DEVICE,VECTOR,1
PLCRACK,0,0			!ʾλ
!쿴ڸֵϵĿ
ETABLE,11,NMISC,53		!ʾѵѷſպϣѹ
ETABLE,22,NMISC,60
ETABLE,33,NMISC,67
ETABLE,44,NMISC,74
ETABLE,55,NMISC,81
ETABLE,66,NMISC,88
ETABLE,77,NMISC,95
ETABLE,88,NMISC,102
PLETAB,33,AVG   
PLETAB,44,AVG   
finish








