!-------------------------------------------------------
!	Modelagem sistema rotativo
!	JHSP - 175976
!-------------------------------------------------------

/CLEAR
/UNITS,SI

!---Parâmetros
PI = acos(-1)

d_e = 12e-3 	!Diametro do eixo
L_e = 600e-3	!Comprimento do eixo

d_d = 95e-3		!Diametro do disco
L_d	= 47e-3		!Comprimento do disco

ro = 7850		!Densidade	                                        
Ee = 200e9		!Módulo de Elasticidade				
nu = 0.3		!Coeficiente de Poisson

alpha = 0		!Amortecimento Proporcional Massa
beta = 1e-4		!Amortecimento Proporcional Rigidez

!---Cálculos Iniciais
m_d = 0.25*PI*(d_d**2-d_e**2)*L_d*ro
Id_d = m_d*(3*((d_d/2)**2-(d_e/2)**2)+L_d**2)/12
Ip_d = 0.5*m_d*((d_d/2)**2-(d_e/2)**2)


/PREP7

!---Elementos utilizados

ET, 1, BEAM188 		!Elemento de viga com dois nós
KEYOPT, 1, 1, 0		!Desativa Warping (6gdl's apenas)
KEYOPT, 1, 3, 0 	!Função de Forma Linear

ET, 2, MASS21		!Elemento de massa estrutural (1 nó)
KEYOPT, 2, 1, 0
KEYOPT, 2, 2, 0
KEYOPT, 2, 3, 0 

!Constantes do disco
R, 2, m_d, m_d, m_d, Id_d, Id_d, Ip_d

!---Materiais
MP, EX, 1, P_E
MP, DENS, 1, P_RO
MP, NUXY, 1, P_nu

!---Seções do eixo
sectype, 1, BEAM, csolid
secdata, d_e/2, 8, 2

!---Nós
N, 01,,, 0.000
N, 02,,, L_e/2
N, 03,,, L_e

!---Criação dos elementos de viga
TYPE, 1
MAT, 1

SECNUM, 1
E, 1,2
E, 2,3

/VIEW, 1, 1, 1, 1
/ESHAPE, 1.0
EPLOT

!---Criação dos discos
TYPE, 2
REAL, 2
E, 2

!---Criação dos mancais
D, 1, all, 0
D, 3, all, 0

FINISH