!Kegelstumpfschale mit Rillen
!aus dem Forum von schabim
!Schätzing, 5.1.2012


finish  $/clear
*abbr,schabim,/input,schabim,txt

!---Parameter----
radu=100e-3    !Radius unten
rado=50e-3     !Radius oben
hoch=100e-3    !Höhe
radri=4e-3     !Radius Rille
abst=2e-3      !Abstand zwischen zwei Rillen
pi=4*atan(1)

/prep7
et,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3


!---Außenkontur------
winkel=atan((radu-rado)/hoch)*180/pi
kante=sqrt((radu-rado)**2+hoch**2)
!----Anzahl der Rillen----
anz=kante/(abst+radri*2)-1
anz=nint(anz)

!---Lokales Koordinatensystem----
local,11,0,radu,0,0,winkel
wpcsys,,11

!---Generieren der Wellenkante----
i=1
k,i,0,0
k,i+1,0,abst
l,i,i+1  !L: i
k,i+2,0,abst+radri
k,i+3,-radri,abst+radri
k,i+4,0,abst+radri*2
larc,i+1,i+3,i+2,radri
larc,i+3,i+4,i+2,radri
!----Vervielfachen------
lgen,anz,all,,,0,abst+radri*2
!---Restlinie------
*get,kmax,kp,,num,max
k,kmax+1,0,kante
l,kmax,kmax+1


!------Rotation-------
wpcsys,,0
csys,0
!----Rotationsachse------
*get,kmax,kp,,num,max
k,kmax+1,0,0,0
k,kmax+2,0,hoch,0
*get,lanz,line,,num,max

*do,i,1,lanz
arotat,i,,,,,,kmax+1,kmax+2
*enddo

!------Vernetzung-----
nummrg,kp
mshkey,1
mshape,0,2d

amesh,all