/CLEAR, START
/PREP7  

!*****************************************
!* nastavení materiálových vlastností    *
!*****************************************

ET,1,SOLID186	!* solid prvek pro modelování
ET,2,COMBIN39	!* prvek nelineární pružina
ET,3,MASS21		!* prvek hmotnost
ET,4,BEAM189	!* nosníkový prvek
ET,5,COMBIN14	!* lineární pružina a nelin. tlumič.

!* hřídel
MP,EX,1,2.1e11 		!* younguv modul pruznosti v tahu
MP,PRXY,1,0.3		!* poissonovo cislo
MP,DENS,1,7800		!* hustota oceli
MP,DAMP,1,.1e11		!* tlumení pružiny (konst)


!* vlastnosti pruziny
K1=1e2			!* Konst. tuhosti 1
m1sc=1e-3		!* měřítko pro graf tuhosti
R,1,-1.5*m1sc,-3.375*K1,-0.9*m1sc,-0.725*K1,-0.5*m1sc,-0.125*K1
RMORE,0.0,0.0,.5*m1sc,.125*K1,.9*m1sc,.725*K1 
RMORE,1.5*m1sc,3.375*K1,2*m1sc,8*K1,3*m1sc,27*K1   
RMORE,4*m1sc,64*K1,5*m1sc,125*K1,6*m1sc,216*K1 
RMORE,7*m1sc,343*K1,8*m1sc,512*K1

!* vlastnosti hmotnosti
R,2,.02,.02,.02, , , ,

!* vlastnosti inerciální soustavy
DOMEGA,1.5,0,0, 	!* úhlové zrychlení soustavy

!* vlastnosti beamu
SECTYPE,   1, BEAM, CSOLID, , 0 
SECOFFSET, CENT 
SECDATA,0.005,0,0,0,0,0,0,0,0,0 

m_alfa=3 			!* koeficient alfa materiálového tlumení (hřídel)
m_beta=.001			!* koeficient beta materiálového tlumení (hřídel)


!*****************************************
!* model hřídele - rozumná frekvence 2800*
!*****************************************

!* model
K,1,0,0,0
K,2,.3,0,0

K,3,0,0.002,0
K,4,0,0,0.002
K,5,0,-0.002,0
K,6,0,0,-0.002

K,7,.3,.002,0
K,8,.3,0,.002
K,9,.3,-.002,0
K,10,.3,0,-.002

L,1,3
L,1,4
L,1,5
L,1,6
L,1,2
L,2,7
L,2,8
L,2,9
L,2,10

KWPLAN,-1,1,3,4
CYL4,0,0,0.03
VOFFST,1,0.30,

N,1,0,0
N,2,0,0
N,3,.3,0
N,4,.3,0

!* meshovani pruziny
Type,2
E,1,2
E,3,4

!* meshovani  valce

ESIZE,0.01,0
MSHAPE,1,3D 
MSHKEY,0

MAT,1             	!*HRIDEL
TYPE,1           
VMESH,1

D,1,ALL
D,3,ALL
LNSRCH,ON

!* svázání posuvů čel hřídele - pokud chcete použít tak odkomentovat
!* ASEL,S, , ,1 
!* NSLA,S,1
!* CP,1,UX,ALL
!* ALLSEL,ALL

!* ASEL,S, , ,2 
!* NSLA,S,1
!* CP,2,UX,ALL
!* ALLSEL,ALL


!*************************************************
!* geometrie pro jednostranně vetknutý nosník    *
!*************************************************

/PREP7
K,1,,,, 
K,2,1,,,
LSTR,1,2    !* vytvoří přímou line mezi dvěma kpointy.
TYPE,4   
MAT,1
REAL,1   
ESYS,0   
SECNUM,1 
LMESH,1
/ESHAPE,1.0 
/EFACET,1 
/REPLOT

!* podmínky vetknutí
DK,1,ALL,
 

!*************************************************
!* výpočet modální analýzy                       *
!*************************************************

/SOL
ANTYPE,MODAL 
MSAVE,0   
MODOPT,LANB,14  
EQSLV,SPAR  
MXPAND,14, , ,1 
LUMPM,1 
PSTRES,0
MODOPT,LANB,14,0,0, ,OFF
  
/STATUS,SOLU
SOLVE   
FINISH 
 
/POST1  
SET,LIST

!*************************************************
!* Výpočet přechodového jevu           			 *
!*************************************************

/SOL
!*  
ANTYPE,4
!*  
TRNOPT,FULL 
LUMPM,1 
!*  
DELTIM,0.05,0,0 
OUTRES,ERASE
OUTRES,ALL,ALL  
TIME,5  
!*  
OUTPR,NSOL,ALL, 
/STATUS,SOLU
SOLVE   
FINISH  


!*************************************************
!* vykreslení time hist. postprocessor  		 *
!*************************************************

/POST26 
FILE,'diplomka2','rst','.'  
/UI,COLL,1  
NUMVAR,200  
SOLU,191,NCMIT  
STORE,MERGE 
FILLDATA,191,,,,1,1 
REALVAR,191,191 
!*  
NSOL,2,2,U,Y, UY_2  
STORE,MERGE 
!*  
NSOL,3,2,U,Z, UZ_3  
STORE,MERGE 
XVAR,1  
PLVAR,2,3,  
FINISH  

!*************************************************
!* harmonická analýza superponovaná na modální	 *
!*************************************************
 
/SOL
ANTYPE,HARMIC
HROPT,full 			!* plná harmonická analýza (další možnosti: reduc - redukovaná, msup - superpozice k modální)
HROUT,OFF			!* výstup uvádět jako reálné a komplexní složky přetvoření
LUMP,0				!* nepoužívat redukci
KBC,1               !* STEP CHANGE FORCE
HARFRQ,0,1000       !* Provádět by se měla analýza do 1000 Hz
NSUBST,200		    !* počet Substepů na LoadStep
OUTPR,NSOL,ALL,
OUTRES,,1
SOLVE

!*************************************************
!* přidání nelineární pružiny  			  		 *
!*************************************************

N,6,1,				!* Přidat keypoint na konci
TYPE,2   			!* Type 2 = combin 39
MAT,1
REAL,1   
ESYS,0   
SECNUM,1 
TSHAP,LINE  		!* liniový element
E,6,2,				!* Vložit element mezi nody 2 a 6
D,6,,0,,,,UX,UY,UZ	!* zamezit posuvům v nodu 6

