! Command File mode of Beam - Kreisvollprofil

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

/title, Solid45 mit Laststeigerung

!------------Geometrie-Parameter

L = 300					! Länge Träger
H = 12					! Höhe Querschnitt
B = 6					! Breite Querschnitt
a = 3					! Element-Größe
Fres = 5				! resultierende Gesamtlast aus Streckenlast, Last muss zu Beginn kleiner sein als max aufnehmbare Last -> es wird nur nach oben hin iteriert / gesteigert 	
F = Fres/(L/a+1)/3				! Last pro Node
grenzwert = 23.5			! maximale aufnehmbare Spannung von S235

/Prep7					!Enter Preprocessor

block,0,B,0,H,0,L

!------------Element Type und Section
ET,1,Solid45


!------------Real Constants			! entfällt für Beam188

!------------Material Properties
	! elastic
MP,EX,1,21000				! E-Modul = 21000 kN/cm^2
MP,PRXY,1,0.3				! Querdehnzahl = 0.3
	! multilinear elastic 
TB,BISO,1,1,2,  
TBTEMP,0
TBDATA,,23.5,2.1,,,,



!------------Meshing
ESIZE,a,0
CM,_Y,VOLU  
VSEL, , , ,       1 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
MSHAPE,0,3d 
MSHKEY,1
VMESH,_Y1   

/solu					! Enters Solution Mode

!------------Analysis Type
ANTYPE,0				! Static Analysis

!------------Displacement				! Auflager
D,16, , , ,0,UX,UY,UZ, , , ,
D,17, , , ,0,UX,UY,UZ, , , ,
D,18, , , ,0,UX,UY,UZ, , , ,
D,2, , , ,0,UX,UY, , , , ,
D,8, , , ,0,UX,UY, , , , ,
D,12, , , ,0,UX,UY, , , , ,

!------------Laststeigerung

SOLU
*DO,i,1,5,1
    varLast= F*i+0.28				! bsp variable Last 'varLast' via *DO index i +0.99 um schneller an Ergebnis zu kommen
	F,16,FY,-varLast
	F,17,FY,-varLast
	F,18,FY,-varLast
	F,2,FY,-varLast
	F,8,FY,-varLast
	F,12,FY,-varLast	        
!    SFA,,,,varLast                             ! aktuelle Last aufbringen
	*DO,j,31,327,1				! Do-Schleife für Last auf Nodes	
		F,j,FY,-varLast			! Streckenlast auf Node 1 bis a aufteilen
	*Enddo

NLGEOM,1
NSUBST,20,1000,1
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
LNSRCH,1
NEQIT,1000 
    SOLVE            			
!------------Auswertung   
 
/POST1
/VIEW,1,1   
/ANG,1  
/REP,FAST 

plns,s,z
*get,smax,plns,,max
*get,smin,plns,,min    
*if,smax,ge,grenzwert,exit
*if,smin,le,-grenzwert,exit  

/solu
*enddo

!---------Endauswertung-------
/post1
smax=smax
smin=smin
Faktor=i
qmax=F*i*300/L