fini
/clear,nostart
/title,composite without contact
/filname, composite 



/PREP7  
ET,1,Solid186

!KEYOPT,1,3,0
KEYOPT,1,2,1

!Material steel
MP,EX,1,210000		
MP,EY,1,210000
MP,EZ,1,210000
MP,PRXY,1,0.3
MP,PRYZ,1,0.3
MP,PRXZ,1,0.3
MP,GXY,1,80769
MP,GYZ,1,80769
MP,GXZ,1,80769
MP,DENS,1,7850.0E-9

TB,BKIN,1,1             ! DEFINE NON-LINEAR MATERIAL PROPERTY FOR STEEL
TBTEMP,0
TBDATA,1,300,0			
tbplot,bkin
/axlab,x,Dehnung
/axlab,y,Spannung

!Material wood
MP,EX,2,1340,		
MP,EY,2,1340,	
MP,EZ,2,16150,	
MP,PRXY,2,0.35,
MP,PRYZ,2,0.04,
MP,PRXZ,2,0.04,
MP,GXY,2,84,	
MP,GYZ,2,840,	
MP,GXZ,2,840,	 
MP,DENS,2,5.82e-10,	

!tubes
cyl4,0,0,119,0,135,360,2500		!wood
cyl4,0,0,136,0,137.5,360,2500		!steel
	

TYPE,   1   		!Wood
MAT,       2
REAL,   
ESYS,       0   
SECNUM, 
!*  
ESIZE,30,0, 		!Elementgröße
CM,_Y,VOLU  
VSEL, , , ,       1 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
TYPE,   1   		!steel
MAT,       1
REAL,   
ESYS,       0   
SECNUM, 
!*  
CM,_Y,VOLU  
VSEL, , , ,       2 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2 

/eshape,1,1
eplot

!Folgend lineare Eigenwertanalyse nach http://www.mece.ualberta.ca/tutorials/ansys/IT/Buckling/Buckling.html
/sol
antype,0			!statische Berechnung
PSTRES,1		

!NLGEOM,1
NSUBST,20,1000,1
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
LNSRCH,1
NEQIT,1000 
 
NSUBST,100,1000,10  
OUTRES,ERASE
OUTRES,ALL,-20

!Lagerung Fußpunkt
DA,7,ux,
DA,7,uy,
DA,7,uz,
DA,1,ux,
DA,1,uy,
DA,1,uz,

!Lagerung Stützenkopf (nur Formholz)
da,2,ux
da,2,uy

solve

!Einslast anteilig auf die Querschnitte
SFA,8,1,PRES,0.006523		!Belastung NUR auf Stahlrohr
SFA,2,1,PRES,0.064619		!Belastung NUR auf Holzrohr

solve
Finish

/SOLU   
ANTYPE,1			!Eigen Buckling
BUCOPT,LANB,1,0,0,CENTER
solve
finish

/SOLU
EXPASS,1
MXPAND,1,0,0,0,0.001,  
solve  

/POST1  
SET,LAST
PLDISP,0
/VIEW, 1 ,1,1,1 
/ANG, 1 
/REP,FAST