15

112Equation Chapter 2 Section 12.- ANALISIS DE MECANISMOS

2.1 INTRODUCCION

El análisis cinemático de mecanismos se refiere al estudio del movimiento de los mismos

(desplazamientos, velocidades, aceleraciones, y sobreaceleraciones de puntos y miembros
constituyentes), sin atender para nada la causa que produce tales movimientos. El análisis
cinemático en sí mismo tiene poca utilidad. Más bien debe ser considerado como un medio para
llegar a resolver el análisis dinámico y por lo tanto el análisis de esfuerzos, así como para
comprobar los resultados obtenidos en la síntesis de mecanismos. El análisis cinemático de
mecanismos puede acometerse por diversos métodos:

2.1.1 MÉTODOS GRÁFICOS

Los métodos gráficos han sido los más antiguos en su aparición y los tradicionalmente empleados.
Pese a su relativa facilidad estos métodos eran sumamente tediosos si se necesitaba determinar
todo el análisis cinemático para todo el ciclo de movimiento del mecanismo, pues cada nueva
posición requiere que el trazo de un sistema completamente nuevo de diagramas vectoriales y
nada de lo hecho anteriormente sirve para esta nueva posición.

2.1.2 MÉTODOS NUMÉRICOS

Los métodos numéricos tienen la ventaja de ser completamente generalistas, ya que no precisan
del desarrollo de un conjunto de ecuaciones para cada mecanismo particular. Por tanto, no se
requiere acceder al programa informático fuente, para cambiar las ecuaciones, cada vez que se
desea analizar un nuevo mecanismos. En base de este método están elaborados los programas
de simulación multicuerpo, SimWise 4D, Working Model 2D, 3D, 4D, Nastran, Cosmos Motion,
Adams, Universal Mechanism, etc.

Figura 2.1. Mecanismo realizado en Msc Adams1

2.1.3 METODOS ANALÍTICOS

1
https://www.google.com/software-simulation-rmb-dynamic-multibody-Adams-Flexible-Body
 16

Los métodos analíticos de empleo mucho más moderno (datan de la década del 50) se apoyan,
tanto en la puesta a su servicio de la matemática como en el uso de la moderna herramienta que
es el computador personal PC, lo cual ha evolucionado la práctica de la ingeniería. Estos métodos
están alcanzando un gran desarrollo y a niveles prácticos desbancaron a los métodos gráficos,
frente a los que presintan indudables ventajas:

 Permiten el análisis en ciclo completo

 Dan soluciones rápidas
 Permiten planteamientos generales
 Su precisión es definitivamente mayor que los métodos gráficos

El método analítico más utilizado es el Método de Raven o de Algebra Compleja

 LOS NUMEROS COMPLEJOS COMO VECTORES

Hay muchos métodos de represintar vectores. Estos se pueden definir en coordenadas polares,
por su magnitud y su ángulo, o en coordenadas cartesianas, mediante las componentes x y y.
Un vector de posición cualquiera puede ser represintado mediante vectores unitarios o mediante
notación de números complejos.

Formato en coordenadas polares Formato en coordenadas cartesianas

R @< θ ^ R sen ( θ ) ^j
R cos ( θ ) i+
R e iθ R cos ( θ ) + R sen ( θ ) i

En este caso la componente en la dirección X se denomina parte real, y la componente en la

dirección Y, recibe termino de parte imaginaria. Este número imaginario se usa en el análisis
como un operador y no como un valor. En la figura 2.2 se muestra el plano complejo en el que el
eje real represinta la dirección de la componente X del vector en el plano, y el eje imaginario
represinta la dirección del componte Y del mismo vector

Figura 2.2. Plano complejo2

Advierta en la figura 2.3. que cada multiplicación por el operador i resulta en una rotación de 90 o en
sintido contraria a las manecillas del reloj del vector.

2
Fuente propia
 17

Figura 2.3. Operador rotación3

Comprobación: Dado el vector

e i .45°
Al multiplica por el operador i se obtiene :

i∙ e i .45 °=i ( cos ( 45 )+ isin ( 45 ) )=i ( √22 +i √22 )=i √22 − √22
Ahora rotando el vector 90 0 en sintido antihorario, se observa que coinciden:

−√ 2 2
e i .(45° +90 °)=cos(135)+isin (135)= +i. √
2 2

La ventaja de utilizar esta notación de números complejos para represintar vectores en el plano
proviene de la identidad de Euler

e i θ=cos ( θ )+i sen (θ)

No hay función más fácil de derivar o integrar, ya que tal función es su propia derivada

d iθ
e → iexp(iθ )
dθ

Se puede utilizar el Mathematica como herramienta para derivar

D [ exp [ Iθ ] ,θ ] ⅈ ⅇⅈθ

Se utiliza esta notación de números complejos para los vectores, con el fin de desarrollar y deducir
las ecuaciones para la posición, velocidad y aceleración de eslabonamientos

2.1.4 WOLFRAM MATHEMATICA

3
Fuente propia
 18

En el estudio a llevarse a cabo, es imprescindible generar los módelos matemáticos por medio de
formulas que expresan las variables cinemáticas, las cuales necesariamente deben ser graficadas,
como por ejemplo: trayectorias, velocidades, aceleraciones o torques. Gracias a la computadora es
posible generar estos gráficos con el auxilo de muchos programas de análisis númerico, tales
como Wolfram Mathematica (Figura 2.4), MathCAD, MATLAB, MAPLE, TK Solver, etc. Los análisis
llevados a cabo en este texto se basarán en la aplicación de Wolfram Mathematica el cual puede
ser descargado gracias a la RED CEDIA, de la página:
https://www.cedia.org.ec/respositorios/wolfram-mathematica

Figura 2.4. Logo Wolfram Mathematica4

2.2 MECANISMO MANIVELA - CORREDERA

4
https://www.wolfram.com/mathematica/
 19

Figura 2.5. Mecanismos biela manivela a) Motor de combustión interna, b) Tamizador de aridos gruesos5, c)
Aplastador de latas6, d) Aplastador de papas
El mecanismo manivela corrdera se trata de un mecanismo capaz de transformar el movimiento
circular de una manivela en un movimiento lineal alternativo o viceversa. Dicho sistema esta
formado por un elemento giratorio denominado manivela que va conectado con una barra rígida
llamada biela, de tal forma que al girar la manivela la corredera se ve obligada a retroceder y
avanzar, produciendo un movimiento alternativo. Es sin duda uno de los mecanismos más
construidos en el mundo, es la base de los motores de combustión interna (figura 2.5 a), bombas
de desplazamiento positivo, compresores, troqueladoras, (figura 2.6) y diferente tipo de
maquinaria.

5
www.proetisa.com
6
http://www.regalosfrikis.com/aplastador-de-latas/
 20

Figura 2.6, a) Troqueladora b) Plegadoras7

Para efectuar el análisis, se establece una ecuación vectorial en base de los vectores de posición:

2.2.1 ANALISIS DEL MOVIMIENTO

Partiendo del siguiente circuito vectorial, dado por la cadena cinemática (Figura 2.7)

Figura 2.7. Circuito vectorial8

r 2+ ⃗
⃗ r 3=⃗
r1 (2.1)

Y utilizando la notación en álgebra compleja donde el eje x es el eje real y el eje y es el imaginario:

r 2 e iθ 2 +r 3 e iθ 3 =r 1 (2.2)

La identidad de Euler se obtiene por expansión en serie de Taylor de e i θ y se genera la siguiente

equivalencia:
e i θ=cos ( θ )+i sin(θ) (2.3)

Reemplazando en la expresión (2.1) se tiene:

r 1 ( cos 0 0+i sin 00 ) =r 2 ( cos ( θ2 ) +i sin ( θ2 ) ) + r 3 (cos ( θ3 ) +i sin(θ3 ))

Igualando la parte real y la parte imaginaria se obtiene el siguiente sistema de ecuaciones no

lineales.
r 1=r 2 cos ( θ 2 ) +r 3 cos ⁡( θ 3) (2.4)

0=r 2 sin ( θ2 ) +r 3 sin ⁡(θ3 ) (2.5)

7
http://www.catalogometalurgico.com/empresas/view/558

8
Fuente propia
 21

θ2 es la variable independiente que varia de 0 a 360o y la solución de las ecuaciones es:

−r 2 (2.6)
θ3 =sin
−1
( r3
sin ( θ 2) )
r 1=r 2 cos ( θ 2 ) + √ r 32−r 22 sin 2 (θ2) (2.7)

2.2.2 GRÁFICOS EN WOLFRAM MATHEMATICA

ClearAll[r2,r3,x,v,a,θ]
r2=150;
r3=300;
x=r2*Cos[θ Degree]+r3*Sqrt[1-(r2/r3*Sin[θ Degree])^2]
θ3=ArcSin[(-r2*Sin[θ Degree])/r3 ]
Plot[x,{θ,0,360}]
Plot[θ3 180/π,{θ,0,360}]

Figura 2.8. Desplazamiento de la corredera y Rotación de la biela

450 30

400 20

350 10

300
50 100 150 200 250 300 350
250 10

200 20

150
30
50 100 150 200 250 300 350

En el grafico (Figura 2.8) del desplazamiento de la corredera se puede determinar el valor de la

carrera efectuando la diferencia entre 450−150=300 que como se ve es igual a 2 r 2. Se puede
hacer uso de la manipulación interactiva que ofrece el programa Wolfram Mathematica (Figura 2.9)

r2=150
Manipulate[Plot[r2 Cos[θ]+r3 Sqrt[1-(r2/r3 Sin[θ])^2],{θ,0,2*π}],{r3,150,300}]
 22

r3

350

300

250

200

150

100

50

1 2 3 4 5 6

Figura 2.9. Desplazamiento de la corredera con la sintencia Manipulate

2.2.3 ANALISIS SIMPLIFICADO DE DESPLAZAMIENTO

Debido a que este mecanismos es uno de los más estudiados y útiles que existen, se ha buscado
maneras de simplificar el tratamiento matemático con el fin de obtener expresiones más simples
que son de mucha ayuda para el análisis vibratorio de motores, ya sean mono o multicilindricos.En
primer lugar se establece un ligero cambio en la nomenclatura de los eslabones:

Figura 2.10. Circuito vectorial9

Es importante indicar que esta simplificación es valida únicamente para eslabonamientos con
relación
r 1
≤
l 3

La expresión ha simplificar es entonces la siguiente:

9
Fuente propia
 23

2
r
√
x ( θ2 )=r cos ( θ 2 ) +l 1− ( l
sin(θ 2) ) (2.8)

Donde se desarrolla el radical mediante la teoría del binomio

n an n n−1 n ( n−1 ) n −2 2 n ( n−1 ) ∙ ( n−2 ) n−3 3 (2.9)

(a+ b) = + ∙ a ∙ b+ ∙a ∙b + ∙ a ∙b + ...
0! 1! 2! 3!

1 −1 (2.10)
r 2 1
1
2 2
∙( )
2 2 2 2

[ (
1−
l
∙sin (θ2) ) ] = 11 + 12 ∙[−( rl ∙ sin (θ )) ]
2
2 +
2
r
[[
∙ − ∙ sin ( θ2 )
l ]] +…

Puesto que el valor r /l es menor que uno puedo omitirse todos los términos menos los dos
primeros, truncándose la serie
2 (2.11)
r
r 2 1
l (
∙sin ( θ 2)
r2 )
[ (
1−
l
∙sin (θ2) )]2
=1−
2
=1− 2 ∙ sin(θ2 )2
2 ∙l

De la identidad trigonométrica

cos ( 2∙ θ2 )=cos ⁡(θ 2)2−sin (θ2)2 =1−2 ∙sin ⁡(θ2 )2 (2.12)

Y reemplazando 2.12 en 2.10

1−cos ⁡( 2θ2 ) (2.13)
sin ⁡(θ2)2 =
2

2 1 (2.14)
r r2
[ ( 1− ∙sin (θ2)
l )]2
=1− 2 ∙(1−cos ⁡(2θ 2))
4 ∙l

Se obtiene el desplazamiento del pistón

(2.15)
r2 r2 r
x ( θ 2 ) =r cos ( θ2 ) +l 1−
[ 4l 2 ( 1−cos ⁡(2 θ 2 ) ) ]
=l−
4l [
+r cos ( θ2 ) + cos ⁡(2 θ2 )
4l ]
Puesto que θ2=ω t , se logra la expresión del desplazamiento en función del tiempo y sus
respectivas derivadas:

r2 r (2.16)
x ( t )=l− +r (cos ( ωt ) + cos ( 2ωt ) )
4l 4l

r (2.17)
v ( t )=−rω(sin ( ωt ) + sin ( 2 ωt ) )
2l
 24

r r (2.18)
(
a ( t )=−rα sin ( ωt ) +
2l )
sin (2 ωt ) −r ω 2(cos ( ωt ) + cos ( 2 ωt ))
l

2.2.3.1 PROGRAMACIÓN EN MATHEMATICA FUNCIONES SIMPLIFICADAS

Derivación simbólica

x=l-r^2/(4*l)+r*(Cos[ω*t]+r/(4*l)*Cos[2*ω*t]);
x'[t]
x''[t]

Graficación

radio=150;
long=300;
freq=100;
ω2=freq*2*π/60;
Plot[long-radio^2/(4*long)+radio*(Cos[ω2*t]+radio/(4*long)*Cos[2*ω2*t]),{t,0,1},AxesLabel-
>{"t","desplazamiento"} ,PlotLabel->{"Desplazamiento"}]
Plot[radio (-ω2 Sin[t ω2]-(radio ω2 Sin[2 t ω2])/(2 long)),{t,0,1},AxesLabel->{"t","velocidad"} , PlotLabel-
>{"Velocidad"}]
Plot[radio (-ω2 Cos[t ω2]-(radio ω2 Cos[2 t ω2])/long),{t,0,1},AxesLabel->{"t","aceleración"}, PlotLabel-
>{"Aceleración"} ]

Los siguientes son los gráficos de desplazamiento,velocidad y aceleración

 25

D e s p la z a m ie n t o V e lo c id a d
d e s p la z a m ie n to v e lo c id a d

450
1500
400
1000

350 500

300 t
0 .2 0 .4 0 .6 0 .8 1 .0
250 500

200 1000

1500
t
0 .2 0 .4 0 .6 0 .8 1 .0

A c e le ra c ió n
a c e le r a c ió n

1000

500

t
0 .2 0 .4 0 .6 0 .8 1 .0
500

1000

1500

2000

2500

Figura 2.11. Desplazamiento, Velocidad y Aceleración

2.2.4 ANALISIS DE VELOCIDAD

Derivando la ecuación vectorial :

(2.19)
r˙1=⃗
⃗ r˙2+ ⃗
r˙3

r˙1 ∙ ei ∙θ= r˙2 ∙ e i∙ θ 2+ r 2 ∙ i∙ ω2 ∙ ei ∙θ 2 + r˙3 ∙ ei ∙θ 3 +r 3 ∙i ∙ ω3 ∙ e i ∙θ 3 (2.20)

r˙1 ∙ ( cos ( 0 ) +i∙ sin ( 0 ) ) =r 2 ∙ i∙ ω2 ∙ ( cos ( θ 2) + i∙ sin ( θ2 ) ) +r 3 ∙ i∙ ω3 ∙ ( cos ( θ3 ) +i ∙sin ( θ 3 ) ) (2.21)

Igualando los términos reales e imaginarios:

r˙1=−r 2 ∙ ω2 sin ( θ2 ) −r 3 ∙ ω 3 ∙ sin ( θ3 ) (2.22)

0=r 2 ∙ ω 2 ∙ cos ( θ2 ) + r 3 ∙ ω3 ∙ cos ( θ3 ) (2.23)

Se llega finalmente un sistema lineal de dos ecuaciones con dos incógnitas, las incógnitas son la
velocidad angular 3 y la velocidad del pistón r˙1=v (θ2 )
 26

−r 2 cos ⁡(θ2 ) (2.24)

ω 3 ( θ 2) =ω 2
r 3 cos ⁡(θ3 ( θ2 ))

sin ⁡(θ3 ( θ2 )−θ 2) (2.25)

v ( θ 2) =r 2 ω2
cos ⁡(θ 3 ( θ 2) )

2.2.4.1 GRAFICAS DE VELOCIDAD EN MATHEMATICA

ClearAll[r2,r3,x,v,a,θ]
r2=150;
r3=300;

Desplazamiento;
x=r2*Cos[θ]+r3*Sqrt[1-(r2/r3*Sin[θ])^2];
θ3=ArcSin[(-r2*Sin[θ])/r3 ];

Velocidad;
ω2=100;
ω2=%*2*π/60;
ω3=-(r2/r3)*ω2*Cos[θ]/Cos[θ3];
v=r2*ω2*Sin[θ3-θ]/Cos[θ3];

Plot[x,{θ,0,6},PlotLabel->"Desplazamiento Pistón"]
Plot[θ3,{θ,0,6},PlotLabel->"Oscilación Angular Biela"]
Plot[ω3,{θ,0,6},PlotLabel->"Velocidad Angular Biela"]
Plot[v,{θ,0,6},PlotLabel->"Velocidad Pistón"]

Del siguiente gráfico, Fig.2.12, se puede sacar como conclusión que las velocidades del pistón es
cero en los puntos muertos 0 y 180 grados, debido a que existe un cambio de dirección, que la
velocidad del pistón no es uniforme y que es máxima en 75 o (1.309 rad) y 285o (4.974 rad),
dependiendo de las proporciones del mecanismo

V e lo c id a d P is t ó n

1500

1000

500

1 2 3 4 5 6
500

1000

1500

Figura 2.12. Velocidad Pistón

 27

Del siguiente gráfico se puede sacar como conclusión que la velocidad angular de la biela es cero
cuando el ángulo de la manivela es de 90 y 270 grados, y que es máxima en 0 y 180 grados

V e lo c id a d A n g u la r B ie la

1 2 3 4 5 6

Figura 2.13- Velocidad Angular Biela

2.2.5 ANALISIS DE LA ACELERACION

Derivando el circuito vectorial nuevamente:

r¨1=⃗
⃗ r ¨3+ ⃗
r¨2 (2.26)

Lo que significa derivar la velocidad 2.27 obteniendose:

r˙1 ∙ ei ∙θ=r 2 ∙ i∙ ω2 ∙ ei ∙θ +r 3 ∙i ∙ ω3 ∙e i ∙θ
2 3 (2.27)

r̈ 1 ∙ ei ∙θ=r 2 ∙ i∙ ( α 2 ∙ e i ∙θ + ω22 ∙i ∙e i ∙θ ) +r 3 ∙ i ∙ ( α 3 ∙ e i∙ θ + ω32 ∙i ∙e i ∙θ )

2 2 3 3 (2.28)

r̈ 1 ( cos ( 0 ) +i sin ( 0 ) )=r 2 i ( α 2 +ω 22 i ) ( cos ( θ2 ) +isin ( θ 2) ) +r 3 i ( α 3 +ω 32 ∙ i ) ( cos ( θ3 ) +i sin ( θ3 ) ) (2.29)

Resolviendo el sistema de ecuaciones lineales se obtienen las aceleraciones

r 2 ∙ ω 22 ∙ sin ( θ2 ) +r 3 ∙ω 3 ( θ 2 )2 ∙ sin(θ3 (θ 2)) (2.30)

α 3 ( θ 2 )=
r 3 ∙ cos(θ 3 ( θ 2) )

ap ( θ 2 )=−r 2 ∙ ω22 ∙ cos ( θ 2 )−r 3 ∙ α 3 ( θ2 ) ∙ sen ( θ3 ( θ2 ) ) −r 3 ∙ ω3 ( θ2 )2 ∙cos (θ3 ( θ2 ) ) (2.31)

Donde ap ( θ 2 ) es la aceleración del pistón r̈ 1

 28

2.2.5.1 GRAFICAS DE ACELERACIÓN CON MATHEMATICA

ClearAll[r2,r3,x,v,α3,ap,θ]
r2=150;
Manipulate[x=r2*Cos[θ]+r3*Sqrt[1-(r2/r3*Sin[θ])^2];
θ3=ArcSin[(-r2*Sin[θ])/r3];
ω3=-(r2/r3)*ω2*Cos[θ]/Cos[θ3];
v=r2*ω2*Sin[θ3-θ]/Cos[θ3[θ]];
α3=(r2*(ω2^2)*Sin[θ]+r3*((ω3)^2)*Sin[θ3])/(r3*(Cos[θ3]));
ap=-r2*(ω2^2)*Cos[θ]-r3*α3*Sin[θ3]-r3*((ω3)^2)*Cos[θ3];
Plot[ap,{θ,0,6},PlotLabel->"Aceleración Pistón"],{r3,140,600},{ω2,10,100}]

Figura 2.14. Aceleraciòn del Pistón

La aceleración angular es la siguiente. Las aceleraciones angulares máximas ocurren en 90 y 270

grados.
 29

Figura 2.15. Aceleraciòn angular de la Biela

Observando la suavidad de las curvas se concluye que mientras mayor sea la relación radio de
manivela a longitud de biela (r /l) peor es el comportamiento del mecanismo. Un funcionamiento
recomendable desde el punto de vista dinámico se obtiene con relaciónes mayores a r /l=0.3

2.2.6 ACELERACIONES DE LOS CENTROS DE GRAVEDAD

Para calcular las aceleraciones de los centros de gravedad se dibujan vectores de posición que se
extienden desde un punto fijo a cualquier punto que nos interese en este caso los centros de
gravedad :

Figura 2.16. Diagrama vectorial para el cálculo de aceleraciones en los centros de gravedad

Posición del centro de gravedad del eslabón 3 (biela)

r CG 3=r 2 ∙ ei ∙θ + r´3 ∙ ei ∙θ
2 3 (2.32)

Derivando la expresión indicada se obtiene la velocidad absoluta del centro de gravedad del
eslabón 3. El vector velocidad del centro de gravedad de 3

vCG 3=( r 2 ∙ ω2 ) ∙ i∙ ei ∙θ +( r´3 ∙ ω 3) ∙i ∙ ei ∙θ

2 3 (2.33)

Donde las expresiones (r ∙ ω) son los valores absolutos de la velocidad tangencial y las
expresiones (i∙ e ¿¿ i∙ θ)¿indican que su dirección es perpendicular al eslabón. El vector
aceleración absoluta del centro de gravedad de 3, se obtiene de igual manera derivando la
expresión anterior
i ∙θ 2 i ∙θ 2 i∙ θ3 i ∙θ 3 (2.34)
a CG3=( r 2 α 2 ) i e − ( r 2 ω 22) e + ( ŕ 3, α 3 ) i e −( ŕ 3 ω 32) e

Aceleración tangencial Aceleración normal

Donde las expresiones (r 22 ) son las aceleraciones normales y las expresiones ( r  ) son las
aceleraciones tangenciales, la expresión i ei2 indica que es perpendicular al radio y la expresión
 30

- ei2 indica que la aceleración normal esta dirigida hacia el centro de rotación. Finalmente las
componentes reales e imaginarias del vector aceleración absoluta

a CG3 x ( θ2 ) =−r 2 ∙ ω 22 ∙ cos ( θ2 ) −r 3 p ∙ ω3 ( θ2 )2 ∙ cos ( θ3 ( θ2 ) ) −r 3 p ∙ α 3 ( θ2 ) ∙sin ( θ3 ( θ 2 ) ) (2.35)

a CG3 y ( θ2 )=−r 2 ∙ ω22 ∙sin ( θ 2) −r 3 p ∙ω 3 ( θ2 )2 ∙ sin ( θ3 ( θ2 ) ) −r p 3 ∙ α 3 ( θ2 ) ∙cos ( θ3 ( θ2 ) ) (2.36)

2.2.6.1 PROGRAMACIÓN EN MATHEMATICA DE LA ACELERACIÓN ABSOLUTA DEL

CENTRO DE GRAVEDAD

ClearAll[r2,r3,x,v,α3,ap,θ]
r2=150;r3=450;
ω2=100*2*π/60;
Manipulate[x=r2*Cos[θ]+r3*Sqrt[1-(r2/r3*Sin[θ])^2];
θ3=ArcSin[(-r2*Sin[θ])/r3];
ω3=-(r2/r3)*ω2*Cos[θ]/Cos[θ3];
v=r2*ω2*Sin[θ3-θ]/Cos[θ3];α3=(r2*(ω2^2)*Sin[θ]+r3*((ω3)^2)*Sin[θ3])/(r3*(Cos[θ3]));
ap=-r2*(ω2^2)*Cos[θ]-r3*α3*Sin[θ3]-r3*((ω3)^2)*Cos[θ3];
aG3x=-r2*(ω2^2)*Cos[θ]-f*r3*α3*Sin[θ3]-f*r3*((ω3)^2)*Cos[θ3];
aG3y=-r2*(ω2^2)*Sin[θ]-f*r3*α3*Cos[θ3]+f*r3*((ω3)^2)*Sin[θ3];
ParametricPlot[{aG3x,aG3y},{θ,0,6.3}],{f,0,1}]

Figura 2.17. Aceleraciones en el centro de gravedad

2.2.7 ANIMACIÓN DEL MECANISMO MANIVELA CORREDERA

El software Wolfram Mathematica ofrece la posibilidad de simular el movimiento del mecanismo

manivela corredera en forma bastante sencilla:
 31

r2=70;
r3=140;

x[θ2_]=r2*Cos[θ2]+r3*Sqrt[1-((r2/r3)*Sin[θ2])^2];

R12={0,0};
R23[θ2_]:={r2 Cos[θ2],r2 Sin[θ2]};
R14[θ2_]={x[θ2],0};

L2[θ2_]:=Graphics[{Red,Thickness[0.01],Line[{R12,R23[θ2]}]}];
L3[θ2_]:=Graphics[{Blue,Thickness[0.01],Line[{R23[θ2],R14[θ2]}]}];

Animate[Column[{Show[{L2[θ2],L3[θ2]},PlotRange->{{-50,150},{-50,150}}],Plot[x[θ],{θ,0,2π},PlotRange->{-
200,200},Epilog->{PointSize[0.05],Point[{θ2,x[θ2]}]},AxesLabel->{Subscript[θ,in],Subscript[x,out]}]}],{θ2,0,2 π}]

Figura 2.17. Animación resultante

2.2.8 TAREAS PARA EL ESTUDIANTE

1. Usando la ecuación generada abajo, hallar una expresión que tome en cuenta el tercer
termino de la serie y mediante identidades trigonométricas halle sus armónicas respectivas
y la formula completa, luego graficar en MATHEMATICA

1 −1
r 2 1
1
2
1 1 r 2
∙( )
2 2 r 2 2

[ (
1− ∙ sen(θ 2)
l )] 2
[(
= + ∙ − ∙ sen ( θ 2 )
1 2 l )] +
2 [(
∙ − ∙ sen ( θ2 )
l ) ] +…
 32

2. Realizar el análisis del mecanismo yugo escocés o biela manivela de longitud infinita a
mano.

3. Realizar el análisis del mecanismo biela manivela excéntrico y todas sus gráficas en
MATHEMATICA, efectuar la animación respectiva