UNIVERSIDAD MAYOR REAL Y PONTIFICIA DE

SAN FRANCISCO XAVIER DE CHUQUISACA

FACULTAD DE TECNOLOGÍA

MATERIA: MEC 334

TÍTULO DE LA PRÁCTICA: DEFENSA EJERCICIO 4-17
UNIVERSITARIOS (AS):

o BARRIOS DURAN ARON PEDRO

o CALVIMONTES MUÑOZ LUIS VALENTIN
o LOPEZ ZULETA HECTOR PABLO

FECHA DE ENTREGA: 16/05/19

SUCRE-BOLIVIA
2019
INTRODUCCIÓN:

En principio, el análisis de un sistema de “n” grados de libertad no es muy diferente del de los grados de libertad;
se pueden emplear vías y métodos análogos para hacer sus análisis y hallar sus soluciones. Sin embargo, debido
al gran número de posibilidades que existen, la labor de encontrar las frecuencias naturales y evaluar las
respuestas, generalmente requiere mucho más trabajo matemático.

LA VIBRACIÓN TORSIONAL:

Ocurre debido al desequilibrio en sistemas giratorios, tales como desalineamiento de un eje de rotación o de
un acoplador débil que permita los movimientos pequeño-indeseados a lo largo del eje de la rotación. Las piezas
se diseñan para hacer girar con una velocidad constante o, se requieren a veces para acelerar o para retrasar.
Pocas vibraciones precipitadas o al azar que a la pieza de rotación experimenta mientras que es en
funcionamiento, el más largo su vida. Muchos componentes torsionales se diseñan con los materiales que
pueden soportar daño torsional de largo plazo, también conocidos como fatiga torsional.

MATRICES:

El uso de matrices en el análisis vibracional no solo simplifica el trabajo, sino que también ayuda a comprender
el procedimiento utilizado en la solución. Esto es particularmente cierto en sistemas de varios grados de
libertad.
Las ecuaciones diferenciales de movimiento de un sistema de n masas puede expresarse en general así:

O se pueden expresar simplemente como:

Donde (M) se llama la matriz inercia y (K) la matriz rigidez.

ECUACIÓN DE MOVIMIENTO:

El movimiento de un sistema vibratorio de “n” grados de libertar se representa por “n” ecuaciones diferenciales
de movimiento, las cuales se pueden obtener por la ley del movimiento de Newton, por la ecuación de Lagrange
o por el método de los coeficientes de influencia. Como las ecuaciones de movimiento no son eternamente
independientes, se necesita la evaluación completa de determinantes de orden “n” para obtener la solución
simultanea de estas ecuaciones. La evaluación de tales determinantes producirá todas las frecuencias naturales
del sistema.
Otros métodos utilizados son: el método Stodola, el método de Holzer y la iteración matricial. Estos son medios
numéricos más directos, que han sido empleados frecuentemente en el análisis de sistemas vibratorios de
múltiples grados de libertad.
ECUACIÓN DE EJERCICIO:

LA REGLA DE CRAMER:
Proporciona la solución de sistemas de ecuaciones lineales compatibles determinados (con una única solución)
mediante el cálculo de determinantes. Se trata de un método muy rápido para resolver sistemas, sobre todo,
para sistemas de dimensión 2x2 y 3x3. Para dimensiones mayores, los determinantes son bastante más
engorrosos.
Recordad que un sistema de ecuaciones puede escribirse en forma matricial como

donde
A es la matriz de coeficientes del sistema,
X es la matriz con las incógnitas,
B es la matriz con los términos independientes de las ecuaciones.
Para poder aplicar Cramer, la matriz A tiene que ser cuadrada y regular (determinante distinto de 0).
La regla de Cramer establece que la incógnita 𝑋𝐾 de la solución del sistema, cuyos coeficientes están en la
columna k de A, es

Donde 𝐴𝐾 es como la matriz A pero cambiando su columna número k por la columna de términos
independientes, B.

Paso 1.- Indicar el sentido del desplazamiento para cada cuerpo (movimiento en mismo sentido).
Paso 2.- Determinar que desplazamiento es mayor 𝜃2 > 𝜃1
Paso 3.- Analizar
𝐾1 : Como uno de sus lados esta empotrado, el otro lado siempre tendrá un torque en sentido opuesto al
movimiento y su valor es 𝑘1 𝜃1 .
𝐾2 : Como 𝜃2 > 𝜃1 , 𝐾2 tendrá torque para 𝐽1 𝑦 𝐽2 .
Para 𝐽1 es el mismo sentido de movimiento, y su valor es 𝐾2 ( 𝜃2 − 𝜃1 )
Para 𝐽2 en sentido contrario al movimiento, y su valor es 𝐾2 ( 𝜃2 − 𝜃1 )
Obtenemos dos ecuaciones para nuestro ejercicion:
EJERCICIO PROPUESTO:

Calcular las frecuencias naturales del sistema torsional que se muestra en la figura 4-17

Figura (4-17): Grafica:

Sumatoria de Momentos en cada masa:

𝜮𝑴𝟏 = 𝑱𝟏 ∗ 𝜽̈𝟏

𝐽1 ∗ 𝜃̈1 + (𝑘2 + 𝑘1 ) ∗ 𝜃1 − 𝑘2 ∗ 𝜃2 = 0

𝜮𝑴𝟐 = 𝑱𝟐 ∗ 𝜽̈𝟐

𝐽2 ∗ 𝜃̈2 − 𝑘2 ∗ 𝜃1 + (𝑘3 + 𝑘2 ) ∗ 𝜃2 − 𝑘3 ∗ 𝜃3 = 0

𝜮𝑴𝟑 = 𝑱𝟑 ∗ 𝜽̈𝟑

𝐽3 ∗ 𝜃̈3 − 𝑘3 ∗ 𝜃2 + 𝑘3 ∗ 𝜃3 = 0

Sistema de Ecuaciones

𝐽1 ∗ 𝜃̈1 + (𝑘2 + 𝑘1 ) ∗ 𝜃1 − 𝑘2 ∗ 𝜃2 = 0

𝐽2 ∗ 𝜃̈2 − 𝑘2 ∗ 𝜃1 + (𝑘3 + 𝑘2 ) ∗ 𝜃2 − 𝑘3 ∗ 𝜃3 = 0

𝐽3 ∗ 𝜃̈3 − 𝑘3 ∗ 𝜃2 + 𝑘3 ∗ 𝜃3 = 0
Movimiento Armonico Simple

𝜃1 = 𝐴 sin(𝑤𝑛 ∗ 𝑡 + ∅)
𝜃2 = 𝐵 sin(𝑤𝑛 ∗ 𝑡 + ∅)
𝜃3 = 𝐶 sin(𝑤𝑛 ∗ 𝑡 + ∅)
Derivando
 𝜃1̈ = −𝐴 𝑤𝑛 2 ∗ sin(𝑤𝑛 ∗ 𝑡 + ∅)

𝜃2̈ = −𝐵 𝑤𝑛 2 ∗ sin(𝑤𝑛 ∗ 𝑡 + ∅)

𝜃3̈ = −𝐶 𝑤𝑛 2 ∗ sin(𝑤𝑛 ∗ 𝑡 + ∅)

Reemplazando, simplificando las funciones y ordenando los valores de 𝜃 𝑦 𝜃̈ se obtiene el sistema de ecuaciones

−𝐽1 ∗ 𝐴𝑤𝑛 2 + 𝐴(𝑘2 + 𝑘1 ) − 𝐵𝑘2 = 0

−𝐽2 ∗ 𝐵𝑤𝑛 2 − 𝐴𝑘2 + 𝐵(𝑘3 + 𝑘2 ) − 𝐶𝑘3 = 0

−𝐽3 ∗ 𝐶𝑤𝑛 2 − 𝐵𝑘3 + 𝐶𝑘3 = 0

Ordenamos los terminos e introducimos a una matriz

(𝑘2 + 𝑘1 − 𝐽1 𝑊𝑛2 ) −𝑘2 0

[ −𝑘2 (𝑘3 + 𝑘2 − 𝐽2 𝑊𝑛2 ) −𝑘3 ]
2
(𝑘3 − 𝐽3 𝑊𝑛 )
0 −𝑘3

Reemplazando k1 k2 k3 y J1 J2 J3 por k 2k 3k y J 2J 3J respectivamente

(5𝐾 − 𝐽𝑊𝑛2 ) −2𝑘 0

[ −2𝑘 (3𝑘 − 2𝐽𝑊𝑛2 ) −𝐾 ]
2
(𝐾 − 3𝐽𝑊𝑛 )
0 −1

Resolviendo la Matriz

(5𝐾 − 𝐽𝑊𝑛2 ) −2𝑘 0

[ −2𝑘 (3𝑘 − 2𝐽𝑊𝑛2 ) −𝐾 ]
2
(𝐾 − 3𝐽𝑊𝑛 )
0 −1
(5𝐾 − 𝐽𝑊𝑛2 ) −2𝑘 0
[ ]
−2𝑘 (3𝑘 − 2𝐽𝑊𝑛2 ) −𝐾
Regla de Cramer:

(5𝐾 − 𝐽𝑊𝑛2 )(3𝑘 − 2𝐽𝑊𝑛2 )(𝐾 − 3𝐽𝑊𝑛2 )

15𝐾 2 − 10𝐾𝐽𝑊𝑛2 − 3𝐾2𝐽𝑊𝑛2 + 2𝐽𝑊𝑛4 (𝐾 − 3𝐽𝑊𝑛2 )

15𝐾 3 − 10𝐾 2 𝐽𝑊𝑛2 − 3𝐾 2 𝐽𝑊𝑛2 + 2𝐾𝐽𝑊𝑛4 - 45𝐾𝐽𝑊𝑛2+30𝐾𝐽𝑊𝑛4+9𝐾𝐽𝑊𝑛4-6𝐽3 𝑊𝑛6

15𝐾 3 − 58𝐾 2 𝐽𝑊𝑛2 + 39𝐾𝐽𝑊𝑛4 + 2𝐾𝐽𝑊𝑛2 - 6𝐽𝑊𝑛6-[(5𝐾 − 𝐽𝑊𝑛2 ) × (−𝐾)(−𝐾) + (𝐾 − 3𝐽𝑊𝑛2 )(−2𝐾)(−2𝐾) ]

15𝐾 3 − 58𝐾 2 𝐽𝑊𝑛2 + 39𝐾𝐽𝑊𝑛4 + 2𝐾𝐽𝑊𝑛2 - 6𝐽𝑊𝑛6-[(5𝐾 − 𝐽𝑊𝑛2 ) × (𝐾 2 ) + (𝐾 − 3𝐽𝑊𝑛2 )(4𝐾 2 ) ]

15𝐾 3 − 58𝐾 2 𝐽𝑊𝑛2 + 39𝐾𝐽𝑊𝑛4 + 2𝐾𝐽𝑊𝑛2 - 6𝐽𝑊𝑛6-[5𝐾 3 − 𝐾 2 𝐽𝑊𝑛2 + 4𝐾 3 − 12𝐾 2 𝐽𝑊𝑛2 ]

15𝐾 3 − 58𝐾 2 𝐽𝑊𝑛2 + 39𝐾𝐽𝑊𝑛4 + 2𝐾𝐽𝑊𝑛2 - 6𝐽𝑊𝑛6-5𝐾 3 ∓ 𝐽𝑊𝑛2 − 4𝐾 3 + 12𝐾 2 𝐽𝑊𝑛2

6𝐾 3 − 45𝐾 2 𝐽𝑊𝑛2 + 41𝐾𝐽2 𝑊𝑛4 - 6𝐽3 𝑊𝑛6 = 0

Reemplazando J y k por sus valores (J=1, k=1)

6 − 45𝑊𝑛2 + 41𝑊𝑛4- 6𝑊𝑛6 = 0

Los resultado de la ecuacion son

𝑘 𝑘
𝜔1 = 0.39 √ 𝜔2 = −0.39 √
𝑗 𝑗

𝑘 𝑘
𝜔3 = 1.1 √ 𝜔4 = −1.1 √
𝑗 𝑗

𝑘 𝑘
𝜔5 = 2.35 √ 𝜔6 = −2.35 √
𝑗 𝑗

Se desprecian los valores negativos y solo se toma como solocuiones los valores positivos.

CODIGO DE MATLAB:
%Extremos Fijos
clc,clear;
format short
%Sistemas Masa-Resorte de Extremos Fijos
display('Introducir las constantes de las Ecuaciones:');
A=input('Introducir el valor del Resorte 1= ');
B=input('Introducir el valor de la Masa 1= ');
C=input('Introducir el valor del Resorte 2= ');
D=input('Introducir el valor de la Masa 2= ');
E=input('Introducir el valor del Resorte 3= ');
F=input('Introducir el valor de la Masa 3= ');
syms W
matriz=[((A+C)-(B*W^2)), -C, 0;-C, (E+C-(D*W^2)), -E;0, -E, (E-(F*W^2))]
DET=det(matriz)
a=input('Introducir el valor de la constante de W^6=');
b=input('Introducir el valor de la constante de W^4=');
c=input('Introducir el valor de la constante de W^2=');
d=input('Introducir el valor de la constante=');
%Ecuacion de tercer grado
wn1=sqrt(((-b/(3*a))-((2^(1/3).*(-b^2+3*a*c))/(3*a*(-2*b^3.+9*a*b*c-27*a^2.*d+sqrt(4*(-b^2.+3*a*c)^3.+(-
2*b^3.+9*a*b*c-27*d*a^2.)^2.))^(1/3)))+((((-2*b^3.+9*a*b*c-27*a^2.*d+sqrt(4*(-b^2.+3*a*c)^3.+(-
2*b^3.+9*a*b*c-27*d*a^2.)^2.))^(1/3))/(3*2^(1/3).*a)))));
Wn1=real(wn1);
wn2=sqrt(((-b/(3*a))+(((1-(sqrt(3))*i)*(-b^2.+3*a*c))/((a*3*2^(2/3))*(-2*b^3.+9*a*b*c-27*a^2.*d+sqrt(4*(-
b^2.+3*a*c)^3.+(-2*b^3.+9*a*b*c-27*d*a^2.)^2.))^(1/3)))-(((1+sqrt(3)*i)*(-2*b^3.+9*a*b*c-
27*a^2.*d+sqrt(4*(-b^2.+3*a*c)^3.+(-2*b^3.+9*a*b*c-27*d*a^2.)^2))^(1/3))/(a*6*2^(1/3)))));
Wn2=real(wn2);
wn3=sqrt(((-b/(3*a))+(((1+(sqrt(3))*i)*(-b^2.+3*a*c))/((a*3*2^(2/3))*(-2*b^3.+9*a*b*c-27*a^2.*d+sqrt(4*(-
b^2.+3*a*c)^3.+(-2*b^3.+9*a*b*c-27*d*a^2.)^2.))^(1/3)))-(((1-sqrt(3)*i)*(-2*b^3.+9*a*b*c-
27*a^2.*d+sqrt(4*(-b^2.+3*a*c)^3.+(-2*b^3.+9*a*b*c-27*d*a^2.)^2))^(1/3))/(a*6*2^(1/3)))));
Wn3=real(wn3);
%Valores de las Frecuencias
disp('Las Frecuencias Naturales son:')
Wn1
Wn2
Wn3
%Grafica
w = -2.5:0.0001:2.5;
y = -6*w.^6 + 41*w.^4 - 45*w.^2 + 6;
z = zeros(size(w));
Z = zeros(size(y));
plot( w,y )
hold on
plot ( w,z, 'r' )
plot ( z,y, 'r' )

SIMULACION EN INVENTOR:
Pasos a seguir:
1. Asignar material
2. Hacer la restricción
3. Introducir las cargas
4. Malla
(La malla es para restringirla geometría de operaciones de trabajo)
Los resultados que aparecen en rojo representan tensión o deformaciones altas, o un factor de seguridad bajo.

RESULTADOS
TENSION DE VON MISES
Es una magnitud física proporcional a la energía de distorsión, se usa en el contexto de las teorías de fallo
como indicador de buen diseño para materiales dúctiles.

PRIMERA TENSIÓN PRINCIPAL

Permite indicar el valor de la tensión que es normal al plano en el que la tensión de corte es cero. Permite
verificar la tensión de elasticidad máxima inducida en la pieza debido a las condiciones de carga.
TERCERA TENSIÓN PRINCIPAL
Permite la actuación normal al plano en el que la tensión de corte es cero. Permite verificar la tensión de
compresión máxima inducida en la pieza debido a las condiciones de carga.

DESPLAZAMIENTO
Los resultados muestran la forma deformada del modelo después de aplicar la solución. Permite indicar la
magnitud de la deformación con respecto a la forma original.
COEFICIENTE DE SEGURIDAD
Permite indicar las áreas del modelo que es probable que fallen bajo la carga.

BIBLIOGRAFIA:
Vibraciones Mecánicas (Balakumar Balachandran)
Seto, William - Vibraciones Mecánicas