! Command File mode of Beam - Kreisvollprofil

/title, Beam Test

! Geometrie-Parameter

L = 300
R = 5					! Radius Querschnitt
N = 12					! Querschnitt unterteilt gradweise
T = 20					! Querschnitt unterteilt von innen nach außen
a = L/20				! Anzahl Nodes in Line = a+1
Fres = 10				! resultierende Gesamtlast aus Streckenlast, Last muss zu Beginn kleiner sein als max aufnehmbare Last -> es wird nur nach oben hin iteriert / gesteigert 	
F = Fres/(a+1)				! Last pro Node
grenzwert = 23.5			! maximale aufnehmbare Spannung von S235

/Prep7					!Enter Preprocessor

! Keypoints
K,1,0,0,0
K,2,L,0,0

! Lines
L,1,2

! Element Type und Section
ET,1,beam189
SECTYPE,1,BEAM,CSOLID
secdata,R,N,T


! Real Constants			! entfällt für Beam188

! Material Properties
	! elastic
MP,EX,1,21000				! E-Modul = 21000 kN/cm^2
MP,EX,1,0.3				! Querdehnzahl = 0.3
	! multilinear elastic



! Mesh Size
LESIZE,ALL, , ,a, ,1, , ,1,

! Line Mesh
!REAL,1
LMESH,ALL				! Mesh Line 1 to 1 in increments of 1

FINISH					! Finishes pre-processing



/solu					! Enters Solution Mode

! Analysis Type
ANTYPE,0				! Static Analysis

! Displacement				! Auflager
DK,1, , , ,0,UX,UY,UZ, , , ,
DK,2, , , ,0,UY,UZ, , , , ,

! Loads	als Schleife auf Nodes aufgeteilt -> Streckenlast

!*DO,i,1,L/a+1,1			! Do-Schleife für Last auf Nodes	
!	F,i,FZ,-F		! Streckenlast auf Node 1 bis a aufteilen
!*Enddo

! Lasten variieren

SOLU
*DO,i,1,6,1
    varLast= F*i                              	! bsp variable Last 'varLast' via *DO index i
!    SFA,,,,varLast                              ! aktuelle Last aufbringen
	*DO,j,1,2*a+2,1				! Do-Schleife für Last auf Nodes	
		F,j,FZ,-varLast			! Streckenlast auf Node 1 bis a aufteilen
	*Enddo
    SOLVE
    /POST1
    SET,LAST                                    ! aktuelles Ergebniss
    !*GET,SMAX,NODE,15,S,X                            ! schreibe aktuelle Spannung [o.a. Grenzwert] in User-Variable 'varSpng'
    !*GET,SMAX,ELEM,15,SDSG
    *IF,SMAX,GE,grenzWert,EXIT         	! Falls gewuenschte Grenzbeanspr. erreicht beende Lasterhoehung
    /SOLU
!    ANTYPE,,REST            			! gegebenenfalls, nichtlineare Analyse
*ENDDO

SOLVE   
 
! Volumenkörper darstellen und Spannung in x-Richtung anzeigen
/SHRINK,0   
/ESHAPE,1.0 
/EFACET,1   
/RATIO,1,1,1
/CFORMAT,32,0   
/REPLOT 
FINISH  
/POST1    
/EFACET,1   
PLNSOL, S,X, 0,1.0  


! Ansicht
/VIEW,1,,-1 
/ANG,1  
/REP,FAST