!aus dem Forum von BauIng27
!Schätzing, 3.8.2016

Finish
/Clear,start
*abbr,bauing27,/input,bauing27,txt

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

!Radien
fak=7/8   !Faktor RadienReduzierung
r0=40/2			!Radius 0
r1=45/2			!Radius 1   
r2=50/2			!Radius 2   
r3=55/2			!Radius 3   
r4=r3
r0a=40/2*fak			
r1a=45/2*fak	   
r2a=50/2*fak	  
r3a=55/2*fak	
r4a=r3a
!Linienaufteilung
ele=30*h/250		!Linienaufteilung vertikal für Netzweite/Elementanzahl
rele=8 			!Linienaufteilung horizontal für Netzweite/Elementanzahl

!Höhen
h0=h			!Höhe 0 = 0
h1=(2/3)*h			!Höhe 1 = 2/3 Stabhöhe
h2=(1/3)*h			!Höhe 2 = 1/3 Stabhöhe
h3=h2/ele*2.4
h4=0

!------------------------------------------------------
/prep7       !Preprozessor starten
!seltol,0.00001
 
!---stark vereinfachte Geometrie----
ptxy,r0,h0,r0a,h0,r1a,h1,r1,h1
poly
ptxy,r1,h1,r1a,h1,r2a,h2,r2,h2
poly
ptxy,r2,h2,r2a,h2,r3a,h3,r3,h3
poly
ptxy,r3,h3,r3a,h3,r4a,h4,r4,h4
poly

ptxy,0,h0,r0a,h0,r1a,h1,0,h1
poly
ptxy,0,h1,r1a,h1,r2a,h2,0,h2
poly
ptxy,0,h2,r2a,h2,r3a,h3,0,h3
poly
ptxy,0,h3,r3a,h3,r4a,h4,0,h4
poly

aglue,all
numcmp,area  !8 Flächen

!------------------------------------------------------
!Elemente definieren:

!Beton
ET,1,SOLID186

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

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


! Volumen durch Rotation erzeugen
vrotat,all,,,,,,20,29,360  !32 Volumen

!Lösch überflüssige Volumen auf rechter Stabhälfte
vdele,13
vdele,21

esize,h/100
vmesh,all