Finish
/Clear,start

!Parameter definieren:
*set,h,980			!Stabhöhe

!Radien
*set,r0,40/2			!Radius 0
*set,r1,45/2			!Radius 1   
*set,r2,50/2			!Radius 2   
*set,r3,55/2			!Radius 3   

!Höhen
*set,h0,0			!Höhe 0 = 0
*set,h1,(2/3)*h			!Höhe 1 = 2/3 Stabhöhe
*set,h2,(1/3)*h			!Höhe 2 = 1/3 Stabhöhe

!Linienaufteilung
*set,ele,30*h/250		!Linienaufteilung vertikal für Netzweite/Elementanzahl
*set,rele,8 			!Linienaufteilung horizontal für Netzweite/Elementanzahl

!Radius/Höhe für Kerbe
*set,r4,(r3)-((r3)/8)		!Radius 4 = Radius 3 - eine Elementbreite
*set,h3,(h0)+((h2)/(ele/2.4))	!Höhe 3 = Höhe 0 + eine Elementhöhe

!------------------------------------------------------

/prep7       !Preprozessor starten

seltol,0.00001 

!Knoten erzeugen

k,101,-r0,h,0       
k,102,-r1,h1,0
k,103,-r2,h2,0
k,104,-r3,h0,0

k,105,0,h,0 
k,106,0,h1,0 
k,107,0,h2,0 
k,108,0,h0,0


k,109,-r4,h0,0
k,110,-r3,h3,0
k,111,-r4,h3,0

k,112,0,h3,0

spline,101,102,103,110,104

!Horizontale Linien

l,105,101
l,106,102
l,107,103
l,112,111
l,111,110
l,108,109
l,109,104

!Vertikale Linien

l,105,106
l,106,107
l,107,112
l,112,108
l,111,109

!Flächen

!!Zwei oberen Hauptflächen aus r0 bis r2 und h bis h2
al,1,5,12,6
al,2,6,13,7

!!Fläche bis Höhe Kerbe von h3 bis h2
al,3,7,14,8,9

!!Fläche Kerbe
al,4,9,16,11

!!Fläche neben Kerbe
al,16,8,15,10
!-----------------------------------------------------

!Elemente definieren:

!Beton
ET,1,SOLID186

!------------------------------------------------------

!Vernetzen
Type,1         !Element festlegen
mat,1          !Material festlegen

!!Linienunterteilung für Netzweite

! Vertikale Linien
lesize,1,,,ele/2.4,
lesize,2,,,ele/2.4,
lesize,3,,,(ele/2.4)-1,    !-1 wegen Kerbe = einer Elementhöhe = eine Unterteilung
lesize,12,,,ele/2.4,
lesize,13,,,ele/2.4,
lesize,14,,,(ele/2.4)-1,    !-1 wegen Kerbe = einer Elementhöhe = eine Unterteilung

! Horizontale Linien
lesize,5,,,rele,
lesize,6,,,rele,
lesize,7,,,rele,
lesize,8,,,rele,
lesize,10,,,rele,

! Verkleben der Flächen
aglue,all

! Volumen durch Rotation erzeugen
vrotat,all,,,,,,105,106,360

!Löscht die Volumen aus Kerbflächen der linken Stabhälfte um Kerbe zu erzeugen
vdele,19,,,0
vdele,4,,,0

vmesh,all