!Ermittlung der Eigenfrequenzen
!einer Glocke
finish  $/clear
!schaetzing  25.12.03
*abbr,glocke,/inp,glocke.inp
!-------------------
rado=0.5        !Radius Oben
wadi=25e-3      !Wand-Dicke
radu=0.8        !Radius unten
hoch=1.6        !Höhe
faktor=0.5
rado=rado*faktor
radu=radu*faktor
wadi=wadi*faktor
hoch=hoch*faktor
!---------------------
/prep7
et,1,solid185
et,2,mesh200,6   !4 Knoten
!---Material1-----Stahl--------
mp,ex,1,210e9
mp,nuxy,1,0.3
mp,dens,1,7850
!----Geometrie----
/color,outl,blac
cyl4,,,rado-wadi,0,rado,90
k,11,radu,rado-hoch
k,12,radu-wadi/2,rado-hoch
bspline,12,4,,,,,1,-1,,0,1  !L5
bspline,11,1,,,,,1,-1,,0,1  !L6
l,4,1                       !L4
l,11,12                     !L7
al,5,4,6,7  !A2
aglue,all
type,2
amesh,1,2
!-----Rotation--------
k,13,0,rado
k,14,0,rado-hoch
type,1  $mat,1
esize,,12
vrotat,1,2,,,,,13,14
!----Hilfselemente-löschen-
aclear,1,2
eplot
!-----Mittelelemente-raus--------------
nsel,s,loc,x,0
nsel,r,loc,z,0
esln
nsle
!-----Randbedingung-----------
d,all,all
allsel
!-----------------
/solu
antype,modal
modopt,lanb,10,10,1000
mxpand,10
solve
!--------------
/post1
set,first
/triad,off
/view,,1,-1,1
pldi,0

!nsel,s,loc,x,0,radu !Schnitt
esln
pldi,0
anmode

/eof
!------------------
set,first
pldi,0
anmode
!------------------
set,next
pldi,0
anmode