!Kegelstumpfschale mit einer Längsrille
!aus dem Forum von schabim
!Schätzing, 7.1.2012


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

!---Parameter----
radu=100e-3    !Radius unten
rado=50e-3     !Radius oben
hoch=100e-3    !Höhe
radri=6e-3     !Radius Rille
pi=4*atan(1)

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

!----Kegelstumpfschale-----
cone,radu,rado,0,hoch,-90, 90
cone,radu,rado,0,hoch, 90,270


!---Außenkontur------
winkel=atan((radu-rado)/hoch)*180/pi
kante=sqrt((radu-rado)**2+hoch**2)

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

!---Generieren der Rille als Zylinder----
cylind,,radri,0,kante*1.1

vsbv,1,3

vdele,all
/pnum,area,1
aplot

wpcsys,,0
csys,0
asel,s,loc,z,0
asel,a,loc,z,hoch
adele,all

asel,all
asel,s,loc,x,0
adel,all
asel,all

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

amesh,all