! Versuch Wasserdruck

fini
/clear

/vup,1,z
/VIEW,1,1,1,1
/TRIAD,LTOP

/PREP7

/graphics,full

L=5.0		! [m]
B=2.0
H=2.0
T=1.0

!Erdbeschleunigung
grav=9.81

!Dichte Wasser
rho=1025.0

mp,dens,1,7850
mp,ex,1,2E11
mp,nuxy,1,0.3

et,1,shell63


mat,1
type,1
r,1,0.004

k,1,0,0,0
k,2,L,0,0
k,3,0,B,0
k,4,L,B,0
k,5,0,0,H
k,6,L,0,H
k,7,0,B,H
k,8,L,B,H

a,1,2,4,3
a,1,5,6,2
a,2,4,8,6
a,3,4,8,7
a,1,3,7,5
a,6,8,7,5

alls
aglue,all

esize,0.1
amesh,all

!Unterseite
nsel,s,loc,z,0
esln,s,1
*get,emax1,elem,,count
*do,i,1,emax1
*get,nelem,elem,,num,min
*get,zloc1,elem,nelem,cent,z
druck=rho*grav*(T-zloc1)
sfe,nelem,,pres,,druck
esel,u,,,nelem
*enddo

alls

!vorne
nsel,s,loc,x,L
nsel,r,loc,z,0,T
esln,s,1
*get,emax2,elem,,count
*do,i,1,emax2
*get,nelem,elem,,num,min
*get,zloc2,elem,nelem,cent,z
druck=rho*grav*(T-zloc2)
sfe,nelem,,pres,,druck
esel,u,,,nelem
*enddo

alls

!hinten
nsel,s,loc,x,0
nsel,r,loc,z,0,T
esln,s,1
*get,emax3,elem,,count
*do,i,1,emax3
*get,nelem,elem,,num,min
*get,zloc3,elem,nelem,cent,z
druck=rho*grav*(T-zloc3)
sfe,nelem,,pres,,druck
esel,u,,,nelem
*enddo

alls

!links
nsel,s,loc,y,B
nsel,r,loc,z,0,T
esln,s,1
*get,emax4,elem,,count
*do,i,1,emax4
*get,nelem,elem,,num,min
*get,zloc4,elem,nelem,cent,z
druck=rho*grav*(T-zloc4)
sfe,nelem,,pres,,druck
esel,u,,,nelem
*enddo

alls

!rechts
nsel,s,loc,y,0
nsel,r,loc,z,0,T
esln,s,1
*get,emax5,elem,,count
*do,i,1,emax5
*get,nelem,elem,,num,min
*get,zloc5,elem,nelem,cent,z
druck=rho*grav*(T-zloc5)
sfe,nelem,,pres,,druck
esel,u,,,nelem
*enddo

alls

/psf,pres,norm,2
/pnum,esys,1

!Randbedingungen

nsel,s,loc,x,L/2
nsel,r,loc,y,B/2
nsel,r,loc,z,0

d,all,ux
d,all,uy

alls
nsel,s,loc,z,H
nsel,r,loc,x,1,2
d,all,uz

alls
nsel,s,loc,z,H
nsel,r,loc,x,3,4
d,all,uz

alls

/solu
solve

/post1
set,last

! DIAGRAMM
! VERSCHIEBUNG IN Z AUF BASIS

path,uz,2
ppath,1,node(0,B/2,0)
ppath,2,node(L,B/2,0)
pdef,uz,u,z
plpa,uz


! DIAGRAMM
! VERSCHIEBUNG IN Z AUF GESAMTHÖHE

path,uz,2
ppath,1,node(0,B/2,H)
ppath,2,node(L,B/2,H)
pdef,uz,u,z
plpa,uz