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,104

!Horizontale Linien

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

!Vertikale Linien

l,105,106
l,106,107
l,107,108

l,112,108
l,111,109

!Linien für Kerbfläche

l,104,109
l,109,111
l,111,110
l,110,104

!Flächen

!!Drei Hauptflächen
al,1,4,9,5
al,2,5,10,6
al,3,6,11,7

!!Fläche Kerbe
al,13,14,15,16

!!Fläche unterster Stabteil (Kerbfläche mit eingeschlossen)
al,7,16,8,12

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

!Elemente definieren:

!Beton
ET,1,SOLID186

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

!Vernetzen des Betonvolumens
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,
lesize,9,,,ele/2.4,
lesize,10,,,ele/2.4,
lesize,11,,,ele/2.4,

lesize,16,,,1			!1 für ein Element
lesize,13,,,1			!1 für ein Element
lesize,12,,,1			!1 für ein Element

!Linienunterteilung horizontale Linien
lesize,4,,,rele,
lesize,5,,,rele,
lesize,6,,,rele,
lesize,7,,,rele,
lesize,8,,,rele,

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

!Lösch überflüssige Volumen auf rechter Stabhälfte
vdele,14,,,1
vdele,15,,,1
vdele,9,,,1
vdele,10,,,1

!Überlappende Volumen linke Stabhälfte
vovlap,5,3
vovlap,20,18
vovlap,4,10
vovlap,19,5

!löscht überflüssige Volumen auf linker Stabhälfte
vdele,14,,,1
vdele,20,,,1
vdele,10,,,1
vdele,15,,,1

!neue Linienunterteilung linke untere Stabhälfte
lesize,50,,,ele/2.4-1
lesize,51,,,ele/2.4-1
lesize,52,,,ele/2.4-1
lesize,83,,,ele/2.4-1

vmesh,all