! Command File mode of Beam

finish   $/clear
*abbr,macnite4,/input,macnite4,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 	
fres=20
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

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

tbplot,biso,1
/wait,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   

!----Randbedingungen----
nsel,s,loc,z,0
nsel,r,loc,y,0
d,all,ux,0
d,all,uy,0
d,all,uz,0

nsel,s,loc,z,l
nsel,r,loc,y,0
d,all,ux,0
d,all,uy,0

!------------Last--------------
i=5
i=6.8

nsel,s,loc,y,h
*get,n_anz,node,,count
varlast=Fres*i/n_anz
f,all,fy,-varlast

!----Lösung---------
/solu
allsel
NLGEOM,1
NSUBST,20,1000,1
OUTRES,ALL,ALL  
AUTOTS,on
LNSRCH,1
NEQIT,1000 
SOLVE            			
!------------Auswertung--------   
/POST1
set,last  !letzter Schritt
set,prev  !vorletzter Schritt
/VIEW,1,1   
/dscale,,1
plns,s,z

/eof

set,first
plns,s,z   $set,next