! Command File mode of Beam - Kreisvollprofil

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


/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
Fres=2e-3				! 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
grenzwert = 0.5e-2			! maximale aufnehmbare Spannung von S235

/Prep7					!Enter Preprocessor

ET,1,beam189
SECTYPE,1,BEAM,CSOLID
secdata,R,N,T
MP,EX,1,21000				! E-Modul = 21000 kN/cm^2
MP,EX,1,0.3				! Querdehnzahl = 0.3

!------Geometrie-----------------
K,1,0,0,0
K,2,L,0,0
L,1,2
LESIZE,ALL, , ,a, ,1, , ,1,
LMESH,ALL

!-------Randbedingungen----------------
DK,1, , , ,0,UX,UY,UZ, , , ,
DK,2, , , ,0,UY,UZ, , , , ,

!-----------Lösung-------------
/solu
*do,i,1,6
    varLast= F*i                              	! bsp variable Last 'varLast' via *DO index i
	*DO,j,1,2*a+1,1				! Do-Schleife für Last auf Nodes	
		F,j,FZ,-varLast			! Streckenlast auf Node 1 bis a aufteilen
	*Enddo
SOLVE
!*enddo
!-------------Auswertung---------------------
/POST1
/dscale,1,1
/eshape,1
/view,,,-1
plns,s,x

*get,smax,plns,,max    

*if,smax,ge,grenzwert,exit

/solu
*enddo

!---------Endauswertung-------
/post1
smax=smax
Fmax=F*i