!Versuch: Schaukeln einer halben Zylinderschale
!Schätzing, 21.1.2012

finish  $/clear
*abbr,schaukeln,/input,schaukeln,inp

rad=60e-3        !Radius
winkel=15        !Start-Winkel
dick=rad/12
!dick=rad/24

!------------------------------
/prep7

!-------------------------
et,1,plane183
et,2,targe169
et,3,conta172,,,,2,3
keyopt,3,10,2
keyopt,3,15,3

r,3,,,0.0001

!mp,mu,1,0.001
mp,ex,1,210e9
mp,nuxy,1,0.3
mp,dens,1,7800

!----Geometrie-----
rectng,-rad,rad,-rad,-rad-dick

local,11,0,,,,winkel
wpcsys,,11
pcirc,rad-dick,rad,180,270
pcirc,rad-dick,rad,270,360
aglue,2,3
wpcsys,,0
csys,0

!--------Vernetzung--------
type,1
mat,1
esize,dick
mshkey,1
amesh,all

!---------------------------
lsel,s,,,5
lsel,a,,,13
nsll,s,1
type,3
real,3
esurf
!-----------------------
lsel,s,,,3
nsll,s,1
type,2
esurf

!-----Randbedingungen--------
lsel,s,line,,1
nsll,s,1
d,all,uy,0
d,all,ux,0

allsel
!-------Last-----
acel,,9.81

!------Lösung---------
/solu
antype,trans
nlgeom,on
kbc,1
outres,all,all
timint,on
deltim,1e-4
time,1
time,.1
!autots,on
allsel
solve

!------------Auswertung------
/post1
/dscale,1,1
/focus,1,0,0,0
/dist,1,0.1
set,first
pldi,2

andata,0.01,,2,,,1,,1