Apuntes sobre Máquinas y Mecanismos
Apuntes sobre Máquinas y Mecanismos
Diciembre 2023
Índice general
1. Introducción 1
1.1. Sobre Máquinas y Mecanismos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2. Algunas definiciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3. Cadena cinemática y mecanismo . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.4. Grados de libertad y pares cinemáticos . . . . . . . . . . . . . . . . . . . . . . . . 3
1.5. Representación esquemática . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6. Restricciones redundantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.7. Inversiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.8. Degeneraciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.9. Montajes y puntos singulares de los mecanismos . . . . . . . . . . . . . . . . . . 9
1.10. Consideraciones sobre la equivalencia de mecanismos . . . . . . . . . . . . . . . . 10
1.11. Ángulo de transmisión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.12. Ley de Grashof . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
i
ÍNDICE GENERAL ii
3.3. Cálculo de la velocidad de un tercer punto conocidas las velocidades de otros dos
de la misma barra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4. Cálculo gráfico de velocidades en un mecanismo de biela manivela . . . . . . . . 61
3.5. Cálculo de velocidades de un mecanismo de retorno rápido . . . . . . . . . . . . . 61
3.6. Cálculo gráfico de velocidades en un actuador lineal . . . . . . . . . . . . . . . . 66
3.7. Cálculo gráfico de velocidades en un mecanismo de leva-seguidor . . . . . . . . . 69
3.8. Centro Instantáneo de Rotacion (CIR) . . . . . . . . . . . . . . . . . . . . . . . . 71
3.9. Teorema de los tres centros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
3.10. Centros instantáneos de rotación en diferentes mecanismos. . . . . . . . . . . . . 72
3.11. Análisis de velocidad de un cuatro barras usando Centros Instantáneos de Rotación 74
5. Dinámica 93
5.1. Conceptos básicos de mecánica . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
5.1.1. Centro de masas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
5.1.2. Cantidad de movimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
5.1.3. Momento cinético del sistema . . . . . . . . . . . . . . . . . . . . . . . . . 94
5.1.4. Ecuaciones de movimiento del sólido rı́gido . . . . . . . . . . . . . . . . . 95
5.2. Teorema de las potencias virtuales . . . . . . . . . . . . . . . . . . . . . . . . . . 96
5.3. Análisis dinámico inverso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.3.1. Análisis dinámico inverso de mecanismos de 1 DOF y con fuerzas de inercia
despreciables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.3.2. Amplificación mecánica . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
5.3.3. Método de superposición . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
5.3.4. Método matricial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
5.3.5. Método de las potencias virtuales . . . . . . . . . . . . . . . . . . . . . . . 106
5.4. Análisis dinámico Directo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
5.4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
5.4.2. Método de Quinn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
5.4.3. Inercia y par reducidos a un eje . . . . . . . . . . . . . . . . . . . . . . . . 109
ÍNDICE GENERAL iii
[Link] 181
10.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 181
10.2. Tensiones cı́clicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
10.3. El experimento de la viga rotatoria . . . . . . . . . . . . . . . . . . . . . . . . . . 184
10.4. Curva S-N . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 184
10.4.1. Fatiga de bajo número de ciclos . . . . . . . . . . . . . . . . . . . . . . . . 186
10.4.2. Fatiga de alto número de ciclos . . . . . . . . . . . . . . . . . . . . . . . . 187
10.4.3. Ejemplo de cálculo de tensiones y ciclos para el ensayo de viga rotatoria . 187
10.5. Factores reductores del lı́mite a fatiga . . . . . . . . . . . . . . . . . . . . . . . . 188
10.5.1. Factor de gradiente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
10.5.2. Factor de acabado de la superficie . . . . . . . . . . . . . . . . . . . . . . 189
10.5.3. Factor de temperatura . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
10.5.4. Factor de corrosión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 191
10.5.5. Factor de fiabilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
10.5.6. Ejemplo de uso de factores reductores de Sn . . . . . . . . . . . . . . . . . 192
10.6. Influencia de la concentración de tensiones . . . . . . . . . . . . . . . . . . . . . . 193
10.6.1. Ejemplo de fatiga con concentradores de tensiones . . . . . . . . . . . . . 194
10.7. Fatiga con Torsor y Cortante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 195
10.7.1. Ejemplo de cálculo a fatiga de una viga sometida a torsor . . . . . . . . . 196
ÍNDICE GENERAL v
Introducción
1
CAPÍTULO 1. INTRODUCCIÓN 2
manera que se restringe la capacidad de movimiento que tendrı́an ambas barras por separado.
En estos apuntes se estudiará un modelo de máquina formada por sólidos indeformables
unidos por unos determinados pares cinemáticos. Con este modelo se pretende responder a
preguntas tales como qué dimensiones debe tener el mecanismo para que cumplan su misión o
qué potencia y par deben impulsarlos para obtener unas determinadas posiciones, velocidades o
aceleraciones.
Si una barra en 2D tiene 3DOF, dos barras tendrán 6. Pero si están unidas con una bisagra,
esos 6 DOF se reducen. Ahora, para definir totalmente la posición de ese conjunto solo es
necesario determinar una de esas barras (3 datos) y añadir el dato del ángulo entre esas barras,
ya que la bisagra impone un punto en común de ambas piezas, el eje de esa bisagra.
Esa bisagra es un par cinemático que ha reducido el número de grados de libertad del
conjunto de las dos barras. Los pares cinemáticos restringen los grados de libertad de las barras
del mecanismo. Hay muchas maneras de conectar diferentes barras, y por tanto habrá muchas
maneras de modelar esas conexiones. A los pares cinemáticos que dejan libres i grados de libertad
se les llama pi . De este modo una bisagra será un par p1 al que también llamaremos par
de revolución, puesto que solo permite la variación del ángulo relativos entre las barras que
conecta.
En los mecanismos hay un barra que se considera fija. Esta barra no tiene capacidad de
movimiento y por tanto ha perdido sus 6DOF (o sus 3DOF si estamos hablando de máquinas
de movimiento plano).
Para calcular los grados de libertad de un mecanismo tridimensional se usa la expresión:
donde n es el numero de DOF del mecanismo y b es el número de barras, contando la barra fija.
En la figura 1.1 se representan algunos pares cinemáticos junto al número de grados de
libertad que permiten. En los siguientes apartados se comentan los más usuales.
Par de revolución
Este par cinemático, ya comentado, modela las conocidas bisagras y algunos tipos de ro-
damientos o cojinetes. Un par de revolución solo deja libre el giro entre las dos barras que
conecta, por tanto restringe 5DOF en mecanismos tridimensionales y 2 en los planos. La figura
1.1 muestra este par cinemático, junto al giro que permite indicados con una flecha.
Par prismático
Un par prismático deja libre el desplazamiento en una dirección, impidiendo tanto el giro
como el desplazamiento perpendicular a la dirección del par prismático. Deja libre 1 DOF y
restringe 5DOF en mecanismos tridimensionales y 2 en los 2D.
Par esférico
Con él se modelan, por ejemplo, las elementos que se encuentran en las barras de la dirección
de un coche o en los extremos de los actuadores en un simulador de vuelo. Este es un par
tridimensional y permite los tres giros entre dos barras, restringiendo los tres desplazamientos
relativos.
CAPÍTULO 1. INTRODUCCIÓN 5
Par helicoidal
Esta es el par que conecta un tornillo y una tuerca. Aunque la tuerca y el tornillo tienen dos
movimientos, el de giro y el desplazamiento en la dirección del eje del tornillo, en realidad es un
par con solo un DOF. Esto es ası́ porque el giro y el desplazamiento no son cualesquiera, sino que
ambos están ligados. El giro de una vuelta de la tuerca obliga a un determinado desplazamiento
del tornillo llamado paso. Este es un par tridimensional.
Par cilı́ndrico
Este par modeları́a la conexión que hay entre un eje y su alojamiento, cuando no hay una
chaveta ni arandelas de entre ellos y se permite el desplazamiento axial. A diferencia del par
anterior, ahora sı́ que hay 2 DOF, estos están asociados al desplazamiento axial y al giro, porque
ambos desplazamientos se pueden efectuar de manera independiente. Es decir, puede haber giro
sin desplazamiento, desplazamiento sin giro o una combinación de ambos.
Par de leva
Este par modela el contacto de dos superficies no conformes en las que se permite el des-
lizamiento. Hay 2DOF asociados al desplazamiento y al giro ya que ambos desplazamientos se
pueden efectuar de manera independiente.
circulo concéntrico o bien un triángulo indica que ese par de rotación une una barra con la barra
fija. Esa barra fija es otra modelización que representa a un sólido que no sufrirá movimiento
independientemente del valor de las fuerzas que actúen sobre él.
1.7. Inversiones
La inversión es la operación según la cual se cambia la barra considerada fija el mecanismo.
El movimiento relativo de las barras no se ve afectado por esta operación, aunque el movimiento
absoluto sı́, ya que este se refiere a la barra considerada fija, y esta ha cambiado con la inversión.
CAPÍTULO 1. INTRODUCCIÓN 8
1.8. Degeneraciones
La degeneración es la operación mediante la cual se cambia de manera drástica la longitud
de alguna barra. Como resultado los pares cinemáticos pueden cambiar de un tipo a otro. En la
figura 1.5 se representa un ejemplo de esto. En la figura 1.5a) se representa un mecanismo de
cuatro barras convencional en donde las barras tienen dimensiones finitas. Pero si la longitud
CAPÍTULO 1. INTRODUCCIÓN 9
O O O
de la barra 4 se alarga mucho, el punto de unión de la barra 3 y 4 pasa a describir arcos de
cada vez mayor radio hasta que, en el lı́mite, la trayectoria de este punto pasa a ser una linea
O
vertical. En esta circunstancia la barra 4 y sus dos pares de rotación podrı́an ser sustituidos
por un par prismático, tal y como se muestra en la figura 1.5c). La cinemática de este cuatro
barras degenerado pasa a ser la del mecanismo biela-manivela.
O
CD. Cuando alguno de los semicı́rculos de a pase a contactar con el arco ED el mecanismo
equivalente será otro. Y si por ejemplo, la curva EC no fuera un arco de cı́rculo y tuviera
diferentes centros de curvatura, cuando el contacto se efectuara en un punto con otro centro de
curvatura, el mecanismo equivalente también serı́a otro.
Para usar la equivalencia de mecanismos en un cálculo cinemático completo, hay que tener
cuidado en asegurarse de que la igualdad cinemática incluya las aceleraciones. Es decir, dos
mecanismos que son iguales en cuanto a sus velocidades pueden no serlo para sus aceleraciones.
Por ejemplo, considérese la figura 1.9. En ella se representa un biela manivela que mueve un
patı́n sobre una guı́a recta a) y otra curva b). Incluso si en la situación de estudio los mecanismos
están en la misma posición, con las mismas dimensiones, y también con la misma dirección del
patı́n, los mecanismos no serán equivalentes, ya que las aceleraciones son diferentes. Las barras
OO
a b y los puntos A y B, si tendrán la misma velocidad en ese instante en los mecanismos 1.9
O
a) y 1.9 b), pero se puede observar que la barra c en el caso a) no tiene velocidad angular,
mientras que en el caso b) sı́. El mecanismo en b) en cambio, sı́ será equivalente al c) tanto a
O
nivel de velocidades como de aceleraciones, incluida la barra c , que tiene la misma velocidad
angular en ambos casos. En general, será con esta equivalencia como se tratarán en estos apuntes
las guı́as curvas.
que existe entre dos vectores. El primero es la velocidad absoluta del punto de la barra de salida
que está unido al acoplador. El segundo vector es la velocidad relativa del punto del acoplador
unido a la barra de salida, respecto a la barra de entrada (es decir, la velocidad de ese punto
considerando la barra de entrada como fija).
Recuérdese, que se denomina velocidad absoluta a la velocidad que mide un observador unido
a la barra considerada fija. Una velocidad relativa, en cambio, es la velocidad que registra un
observador móvil. Lo usual en este libro será que se utilicen letras mayúsculas para nombrar
vectores registrados por el observador fijo, y las minúsculas para los de los observadores móviles.
También se usará la negrita para denotar las magnitudes vectoriales en vez de añadir encima
una pequeña flecha a la letra en cuestión.
En la figura 1.10 se muestran dos ejemplos de ángulos de transmisión, notado con la letra
griega µ.
En el mecanismo de biela manivela de la figura 1.10 a), el movimiento es generado por el
CAPÍTULO 1. INTRODUCCIÓN 13
O
pistón marcado como 2 , por lo que este es considerado como la barra de entrada. Éste se
O O
transmite al cigüeñal, que es la barra de salida 4 mediante la barra 3 . En el mecanismo de
cuatro barras de la figura 1.10 b), el movimiento es generado por la barra OA (etiquetada como
O O
2 ), y se transmite mediante el acoplador 3 a la barra de salida QB 4 . O
O O
Para ambos esquemas, el punto B pertenece a la barra 4 y el C a la barra 3 . La velocidad
absoluta del punto B será perpendicular al segmento QB (representado por VB ), mientras que
la velocidad relativa del punto C respecto a la barra de entrada, será perpendicular al segmento
AC (representada por el vC ).
El ángulo de transmisión µ será el ángulo entre los vectores VB y vC . Para una transmisión
óptima del movimiento este ángulo debe estar próximo a π/2, debiéndose evitar que µ rebase
los lı́mites del intervalo [π/4, 3π/4].
Cuando el ángulo de transmisión es 0 o π, el mecanismo se dice que está en punto muerto. En
esta situación, cualquier esfuerzo en la barra de entrada no producirı́a movimiento en la barra
de salida.
l+s≤p+q (1.3)
Para ilustrar esta condición se utilizará la figura 1.11. En ella se muestran dos mecanismos
de cuatro barras. En el mecanismo en 1.11 a) las dimensiones de las barras son:
En este mecanismo 420 + 100 ≤ 300 + 400 y como se puede comprobar, en las articulaciones
en O y A se producen giros completos.
CAPÍTULO 1. INTRODUCCIÓN 14
Figura 1.11: Ley de Grashof: a) Mecanismo con giros completos en dos rótulas. b)Mecanismo sin giros
completos.
El mecanismo en 1.11 b) tiene las siguientes dimensiones para sus barras dimensiones de las
barras son:
En este mecanismo 440 + 100 6≤ 200 + 200 y en él ninguna articulación puede realizar giros
completos.
Capı́tulo 2
2.1. Introducción
Un vector en dos dimensiones no es más que un segmento orientado dibujado en un papel,
una pizarra o una pantalla de ordenador. Este entidad nos puede servir para representar posi-
ciones, velocidades, aceleraciones, fuerzas... Este segmento también se suele mostrar con un par
de números, que serı́an dos cantidades en dos direcciones diferentes y perpendiculares que usual-
mente son llamadas Eje X y Eje Y. Ası́, para un cierto vector r, para el que estas cantidades
son rX y rY este vector se representa según
!
rX
r=
rY
Las ecuaciones que ligarı́an la longitud del vector r con este par de números, ası́ como el
ángulo θ con alguna referencia (usualmente el eje X positivo) serı́an
q
ρ = rX 2 + r2
Y
rY
θ = tan−1
rX
El hecho de que solo sean necesarios dos cantidades para describir un vector en 2D permite
utilizar otra conocida entidad, el número complejo. De esta manera el complejo r pasa a ser:
r = Re(r) + i Im(r)
donde
Re(r) = ρ cos θ
Im(r) = ρ sin θ
r = ρ eiθ (2.1)
15
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 16
Para hacer más legible la ecuación, se denominará ω a θ̇ y α a θ̈. Debe tenerse en cuenta
ω = θ̇ y α = θ̈ no tienen por qué corresponder con alguna velocidad y aceleración angular. Esto
solo será ası́ si el vector con ese θ gira al unı́sono con alguna barra.
Con esta nueva notación las ecuaciones anteriores pasarán a ser
ṙ = ρ̇ eiθ + ρω ieiθ
(2.2)
r̈ = (ρ̈ − ρω 2 ) eiθ + (ρα + 2ρ̇ω) ieiθ
Si el vector rA une un punto fijo y otro punto móvil A, entonces la velocidad de ese punto será:
VA = ṙA . Análogamente la aceleración de A será: AA = r̈A . También se podrı́a utilizar una
cadena de n vectores para alcanzar el punto A. Entonces, si
n V = Pn ṙ
X A i=1 i
rA = ri ⇒
A = Pn r̈
i=1 A i=1 i
Usando los complejos todos los vectores están referenciados a unos únicos ejes, y por tanto a
un único observador, el observador fijo. Por esta razón las derivadas de vectores que representen
posición siempre serán velocidades absolutas, y si se deriva dos veces, representarán aceleraciones
absolutas.
Para resolver la cinemática de un mecanismo se debe tener en cuenta que los mecanismos
están formados por cadenas cinemáticas. Éstas son conjuntos de barras unidas mediante pares
cinemáticos. Los complejos conecatarı́an los puntos representativos de esas barras formando
bucles o lazos. Cada uno de estos bucles se puede escribir como una ecuación compleja. Estos
lazos de vectores se usarán para resolver los tres problemas cinemáticos de un mecanismo, que
son el cálculo de posición, velocidad y aceleración.
El problema de posición es el primero a resolver a la hora de estudiar la cinemática de
un mecanismo. Consiste en la obtención de los ángulos de las diferentes barras respecto a una
determinada referencia, ası́ como el cálculo de los puntos interesantes del mecanismo respecto
a un determinado origen. En el problema de posición la geometrı́a de las diferentes barras del
mecanismo es conocida, incluyendo la posición de los diferentes pares cinemáticos que unen las
barras. Para resolver el problema de posición deberá proporcionarse el mismo número de datos
iniciales que grados de libertad posea el sistema. Un mecanismo de 1 DoF necesitará un dato
inicial, otro de 2 DoF, requerirá una pareja de datos iniciales, y ası́ sucesivamente. Con todos
estos datos se usarán ecuaciones de lazo complejas para determinar los ángulos y posiciones
relevantes del mecanismo. Para obtener una solución se deberán tener la mitad de ecuaciones
complejas que el número de incógnitas, ya que de cada ecuación compleja se obtienen otras dos,
correspondientes a la parte real e imaginaria.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 17
rD = rA + ri + rj − rk
Figura 2.2: Mecanismo de cuatro barras: a) Esquema. b) Vectores propuestos para su resolución.
barras, lo que incluye la posición de los pares cinemáticos y de los puntos de interés de cada
barra. También debe ser un dato la posición velocidad y aceleración de alguna barra. A esta
barra la llamaremos barra de entrada o barra conductora. También se debe especificar cual es
el montaje realizado del mecanismo, ya que puede haber más de uno, como se comentó en el
primer capı́tulo. El esquema de un mecanismo de cuatro barras se muestra en la figura 2.2. El
doble cı́rculo simboliza un par de rotación ente dos barras en el que una de ellas es la barra fija.
moderno. Realizando un análisis numérico se deben obtener dos soluciones que corresponden a
los dos posibles montajes del mecanismo de cuatro barras. Si, por ejemplo, las longitudes de las
√
barras fueran ρ1 = 0.36 m, ρ2 = 0.1 2 m, ρ3 = 0.4 m y ρ4 = 0.24 m, para θ2 = 120◦ y θ1 = 0, se
obtendrı́a las dos siguientes parejas de soluciones:
Por tanto
ρ2 sin θ2
θ5 = arctan
ρ2 cos θ2 − ρ1
p
ρ5 = (ρ2 cos θ2 − ρ1 )2 + (ρ2 sin θ2 )2
El otro posible montaje deriva de las múltiples soluciones que tiene la función arc cos(x).
Concretamente como arc cos(x) = α y como arc cos(x) = −α.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 20
Figura 2.4: Acoplador con un tercer punto de estudio: a) Esquema. b) Vectores propuestos para el cálculo
de la posición velocidad y aceleración del punto C.
Una vez calculado θ4 , sea con un montaje u otro, θ3 puede obtenerse con
ρ4 sinθ4 − ρ2 sin θ2
θ3 = arcsin
ρ3
r6 = r2 + r7 + r8 (2.5)
El análisis de velocidad consistirá en, conocida la posición de todas las barras del mecanismo,
y dada una velocidad angular en la barra de entrada, el cálculo de la velocidad angular de todas
las demás barras y la velocidad de los puntos de interés. Todas estas velocidades estarán referidas
a la barra fija y por tanto serán velocidades absolutas, las que permiten calcular los complejos.
Para esto derivaremos la ecuación (2.3), que queda
Para resolver el problema de velocidad se puede derivar el lazo de control, o también se puede
sustituir ri por sus valores segun la expresión (2.2), suprimiendo los elementos que proceda. En
este caso, como los complejos unen la posición de los pares de rotación en cada barra, y la
distancia ente ellos no cambia, esto lleva a que todos los sumandos con ρ̇i sean igual a cero.
Además θ̇1 también se anulará al ser el giro de la barra fija. De este modo
ρ̇2 eiθ2 + ρ2 θ̇2 ieiθ2 + ρ̇3 eiθ3 + ρ3 θ̇3 ieiθ3 = ρ̇1 eiθ1 + ρ1
θ̇1 ieiθ1 + ρ̇4 eiθ4 + ρ4 θ̇4 ieiθ4
que usando (2.1) puede dividirse en dos ecuaciones escalares que corresponderı́an a la parte real
y la imaginaria de (2.6).
Este sistema de ecuaciones, con ω3 y ω4 como incógnitas, es lineal y fácilmente resoluble por
la regla de Cramer
Si se quisiera conocer la velocidad de un punto de alguna de las barras se utilizará una cadena
de vectores que uniese el origen en O con el punto en cuestión. Si se desea calcular la velocidad
del punto C del acoplador definido anteriormente usarı́amos la expresión (2.5) que al derivarse
da
Usando 2.2 y eliminando todos los términos con ρ̇i esta ecuación quedarı́a
ρ3 iθ3
ṙ6 = ω2 ρ2 ieiθ2 + ω7 ρ7 ieiθ7 + ω8 ρ8 ieiθ8 = ω2 ρ2 ieiθ2 + ω3 ie − ω3 ρ8 eiθ3
2
donde se ha tenido en cuenta que:
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 22
ρ8 eiθ8 = ρ8 ieθ3 .
En la ecuación (2.7) todo a la derecha del signo igual es conocido, de esta forma sustituyendo
los valores de ρ2 , ρ7 , ρ8 , ω2 , ω3 , θ2 y θ3 se obtendrı́a la velocidad de C.
de ella, eliminando todos los términos con ρ̇i y ρ̈i , ya que ninguno de los complejos utilizados
en la expresión (2.8) cambia su módulo, queda entonces
− ω22 ρ2 eiθ2 + α2 ρ2 ieiθ2 − ω32 ρ3 eiθ3 + α3 ρ3 ieiθ3 = −ω42 ρ4 eiθ4 + α4 ρ4 ieiθ4 (2.9)
|{z} |{z} |{z}
dato ? ?
y de la parte imaginaria:
y usando (2.2) y eliminando todos los sumandos con ρ̇i y ρ̈i esta ecuación quedarı́a
ρ8 eiθ8 = ρ8 ieθ3 .
entonces
ρ3 iθ7 ρ3
AC = −ω22 ρ2 eiθ2 + α2 ρ2 ieiθ2 − ω72 e + α7 ieiθ7 − ω32 ρ8 ieiθ3 − α3 ρ8 eiθ3
2 2
Ası́ la aceleración de C, que notaremos como AC , quedarı́a como
( )
−ω22 ρ2 cos θ2 − α2 ρ2 sin θ2 − ω32 ρ23 cos θ3 − α2 ρ23 sin θ3 + ω32 ρ8 sin θ3 − α3 ρ8 cos θ3
AC =
−ω22 ρ2 sin θ2 + α2 ρ2 cos θ2 − ω32 ρ23 sin θ3 + α2 ρ23 cos θ3 − ω32 ρ8 cos θ3 − α3 ρ8 sin θ3
En la ecuación anterior todo a la derecha del signo igual es conocido, de esta forma la
aceleración AC se obtendrı́a sustituyendo los valores de ρ2 , ρ3 , ρ8 , θ2 , θ3 , ω2 , ω3 , α2 y α3 .
1 function ou = cinem4b ( i n )
2 r5 = sqrt ( ( i n . r 2 * cos ( i n . t 2 )= i n . r 1 )ˆ2+( i n . r 2 * sin ( i n . t 2 ) ) ˆ 2 ) ;
3 t5 = atan2 ( ( i n . r 2 * sin ( i n . t 2 ) ) , ( i n . r 2 * cos ( i n . t 2 )= i n . r 1 ) ) ;
4 ou . t 4 = t5 =acos ( ( i n . r 4 ˆ2+ r 5 ˆ2= i n . r 3 ˆ 2 ) / ( 2 * i n . r 4 * r 5 ) ) ;
5 ou . t 3 = acos ( ( i n . r 4 * cos ( ou . t 4 )= r 5 * cos ( t 5 ) ) / i n . r 3 ) ;
6
7 mv11 = = i n . r 3 * sin ( ou . t 3 ) ;
8 mv12 = i n . r 4 * sin ( ou . t 4 ) ;
9 mv21 = = i n . r 3 * cos ( ou . t 3 ) ;
10 mv22 = i n . r 4 * cos ( ou . t 4 ) ;
11 tiv1 = i n . r 2 * i n . w2 * sin ( i n . t 2 ) ;
12 tiv2 = i n . r 2 * i n . w2 * cos ( i n . t 2 ) ;
13 ou . w3 = det ( [ t i v 1 , mv12 ; t i v 2 , mv22 ] ) / det ( [ mv11 , mv12 ; mv21 , mv22 ] ) ;
14 ou . w4 = det ( [ mv11 , t i v 1 ; mv21 , t i v 2 ] ) / det ( [ mv11 , mv12 ; mv21 , mv22 ] ) ;
15
16 ma11 = i n . r 3 * sin ( ou . t 3 ) ;
17 ma12 = = i n . r 4 * sin ( ou . t 4 ) ;
18 ma21 = = i n . r 3 * cos ( ou . t 3 ) ;
19 ma22 = i n . r 4 * cos ( ou . t 4 ) ;
20 t i a 1 = = i n . r 2 * i n . w2ˆ2 * cos ( i n . t 2 )= i n . r 2 * i n . a2 * sin ( i n . t 2 ) . . .
21 = i n . r 3 * ou . w3ˆ2 * cos ( ou . t 3 )+ i n . r 4 * ou . w4ˆ2 * cos ( ou . t 4 ) ;
22 t i a 2 = = i n . r 2 * i n . w2ˆ2 * sin ( i n . t 2 )+ i n . r 2 * i n . a2 * cos ( i n . t 2 ) . . .
23 = i n . r 3 * ou . w3ˆ2 * sin ( ou . t 3 )+ i n . r 4 * ou . w4ˆ2 * sin ( ou . t 4 ) ;
24 ou . a3 = det ( [ t i a 1 , ma12 ; t i a 2 , ma22 ] ) / det ( [ ma11 , ma12 ; ma21 , ma22 ] ) ;
25 ou . a4 = det ( [ ma11 , t i a 1 ; ma21 , t i a 2 ] ) / det ( [ ma11 , ma12 ; ma21 , ma22 ] ) ;
26 end
1 >> Cinem4bar
2 output =
3 s t r u c t with f i e l d s :
4 t4 : 1.7678
5 t3 : 0.2861
6 w3 : 0 . 1 1 3 9
7 w4 : 0 . 5 7 5 0
8 a3 : 0 . 2 5 2 1
9 a4 : 0 . 4 2 8 0
Si se llama a la función de cálculo con diferentes datos iniciales se pueden obtener las gráficas
con las propiedades cinemáticas del acoplador y la barra de salida. Aquı́ hay que tener en cuenta
que la velocidad y la aceleración de entrada están relacionadas, concretamente ω̇2 = α2 . Por
esto, quizás lo más útil sea definir un vector de tiempos y asignar a cada tiempo un θ2 , un ω2 y
una α2 .
Por ejemplo si el ángulo de entrada θ2 toma la forma:
2π 1
θ2 = + t + t2
3 2
entonces, derivando la expresión anterior se obtiene la velocidad angular y la aceleración angular
de entrada.
θ̇2 = ω2 = 1 + t
θ̈2 = α2 = 1
11 f o r i =1:1000
12 tiempo ( i ) = t ;
13 input . t 2 = 2 * pi/3+ t+t ˆ 2 / 2 ;
14 input . w2 = 1+t ;
15 output = cinem4b ( input ) ;
16 omega3 ( i ) = output . w3 ;
17 omega4 ( i ) = output . w4 ;
18 t = t+t m p o f i n / 1 0 0 0 ;
19 end
20
21 hold on ;
22 plot ( tiempo , omega3 ) ;
23 plot ( tiempo , omega4 ) ;
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 26
El resultado se muestra en las figuras 2.5. En ella se puede apreciar como conforme pasa el
tiempo los módulos de ω3 y ω4 van creciendo. Esto es debido a la aceleración de la barra de
entrada α2 . A menudo los mecanismos trabajan a velocidad de entrada constante, en este caso
las gráficas de las velocidades de salida toman forman cı́clicas como las que se muestran en la
figura 2.6.
Figura 2.5: Velocidades angulares ω3 y ω4 del mecanismo de 4 barras para unos datos de entrada
ω2 = 1 rad/s y α2 = 1 rad/s2 .
Figura 2.6: velocidades angulares ω3 y ω4 del mecanismo de 4 barras para unos datos de entrada ω2 =
1 rad/s y α2 = 0.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 27
Figura 2.7: Cinemática de un cuatro barras:a) Componentes de los complejos usados. b) multiplicación
de un complejo por i, y sus coordenadas.
En la figura 2.9b) se puede apreciar que las coordenadas del complejo girado son las mismas
salvo que se intercambian las partes reales e imaginarias y se le cambia el signo a una de ellas.
De este modo ( ) ( )
iθ2 −0.0707 iθ2 −0.1225
r2 = ρ2 e = ⇒ i r2 = ρ2 ie =
0.1225 −0.0707
Sabiendo esto se pueden escribir directamente las coordenadas de los conjuntos ρj eiθj y ρj ieiθj .
En el mecanismo de cuatro barras que nos ocupa, se sustituirán ṙj = eiθj + ωj ρj ieiθj , ya que
ρ̇j
en este problema todos los vectores tienen módulo constante. Con ello se llega a
Ahora es sencillo sustituir los complejos por sus coordenadas, con lo que se llega a
( ) ( ) ( )
−0.1129 −0.2354 −0.1225
−ω3 + ω4 = ω2
0.3837 −0.0470 −0.0707
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 28
que se puede reescribir como el siguiente sistema del que se obtendrán las incógnitas ω3 y ω4 ,
teniendo en cuenta el dato ω2 = 1 rad/s.
" #( ) ( ) (
0.1129 −0.2354 ω3 −0.1225 ω3 = 0.1138 rad/s
= ⇒
−0.3837 −0.0470 ω4 −0.707 ω4 = 0.5750 rad/s
En ella se sustituirá r̈j teniendo en cuenta que ρ̇j = ρ̈j = 0 para todos los vectores de la ecuación.
De este modo r̈j = −ωj2 ρj eiθj + αj ρj ieiθj , obteniéndose
−ρ2 ω22 eiθ2 + ρ2 α2 ieiθ2 − ρ3 ω32 eiθ3 + ρ3 α3 ieiθ3 = −ρ4 ω42 eiθ4 + ρ4 α4 ieiθ4
que si se reordena para tener los términos multiplicados por incógnitas en un miembro y los
términos conocidos en el otro, lleva a
En esta ecuación se reemplazarán ρ2 eiθ2 , ρ2 ieiθ2 , ρ3 eiθ3 , ρ3 ieiθ3 , ρ4 eiθ4 y ρ4 ieiθ4 por sus coor-
denadas para obtener el sistema
( ) ( ) ( ) ( ) ( ) ( )
−0.1129 −0.2354 2 −0.0707 −0.1225 2 0.3837 2 −0.0470
α3 −α4 = ω2 −α2 +ω3 −ω4
0.3837 −0.0470 0.1225 −0.0707 0.1129 0.2354
donde todo en el lado izquierdo es conocido ya que ω2 y α2 son datos, y ω3 y ω4 se han calculado
previamente en el análisis de velocidad. Sustituyendo queda
" #( ) ( ) (
−0.1129 0.2354 α3 0.0723 α3 = 0.2521 rad/s2
= ⇒
0.3837 0.0470 α4 0.1168 α4 = 0.4281 rad/s2
Figura 2.8: Diferentes pares cinemáticos con su equivalencia a cadenas de pares de rotación.
un total de 5 DOF = 2 × 3 − 1 × 1 y las tres barras y los dos pares de rotación equivalentes
tendrı́an 5 DOF = 3 × 3 − 2 × 2.
En la 2.8 c) se representa otro par cinemático P2, en este caso un contacto con deslizamiento.
O
La cadena con pares de rotación equivalente añade otra barra, la e y dos pares de rotación en
los dos centros de curvatura de las superficies en el punto de contacto. De este modo también
se mantiene el número de grados de libertad.
En la figura 2.9 a) se muestran dos barras unidas mediante un par prismático. Es habitual
esta configuración en la que hay sendas rótulas en los extremos de las barras de manera que
el eje de la guı́a movimiento pasa por ambos puntos cada uno perteneciente a cada una de las
O O
barras. En este ejemplo A ∈ a y B ∈ b . En este caso se añade un único vector que representa a
ambas barras, siendo el eje de la guı́a y el complejo añadido colineales. Este vector ri añadido,
a diferencia del caso anterior, no tiene su módulo constante y por tanto ρi 6= 0. Si por ejemplo
esta configuración modelase un cilindro hidráulico en que se inyecto aun fluido de manera que el
velocidad relativa de la camisa del pistón y el émbolo tiene un valor determinado v mm/s, este
valor serı́a v = ρ̇i .
En la figura 2.9 b) la guı́a del par prismático no es colineal con las rótulas en los puntos A y
B. Esta configuración se intenta evitar en la práctica ingenieril puesto que añade tensiones en los
componentes mecánicos encargados de la realización material de esta restricción (rodamientos
lineales, patinetes, guı́as), aunque en ocasiones pueda ser imposible evitar esta configuración.
Independientemente de las razones por las que el sistema tenga esta geometrı́a es necesario una
correcta modelización cinemática. En este caso para salvar esta falta de colinealidad además del
vector ri , análogo al caso anterior, se añade un vector rj perpendicular en todo momento a ri tal
y como se muestra en la figura 2.9 b). Este complejo rj tendrá módulo constante y la variación
de su argumento será igual a la ri , de manera que θ̇i = θ̇j . A efectos cinemáticos, puesto que
O O
el par prismático obliga a una traslación entre las barras a y b da igual si el complejo ri es
colineal con la guı́a o si el complejo rj se coloca en el punto A o B.
En la figura 2.9 c) hay un par de pin ranura P2, donde el pin está en el punto C. Para el
modelado de esta configuración debe haber un vector que llegue hasta este punto C. En este
O
caso, como a es la barra con la ranura, el par de complejos perpendiculares ri y rj deberán
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 30
O
partir de A ∈ a , aunque también será indiferente el orden de estos complejos. En el caso de que
la linea de la ranura pase por A, no se añadirı́a el complejo rj .
En la figura 2.9 d) aparece un par cinemático P2 de contacto con deslizamiento entre una
O O
superficie curva b y una plana a (si ambas superficies fueran planas se tendrı́a uno de los
casos anteriores). Para ello se añadirá el complejo ri de módulo variable, similar al añadido en
los casos anteriores. El complejo ri será paralelo a la superficie plana. Perpendicular a este se
O
añade el complejo rj que enlaza el centro de curvatura del sólido b con el complejo ri . Este
complejo tendrá módulo constante e igual al radio de curvatura, al menos hasta que cambie la
curvatura de la superficie de contacto en el contacto.
donde ωb = θ̇b , ωc = θ̇c , Rb = |O1 C| y Rc = |O2 C|. En el caso de que los centros de curvatura
del mecanismo no estén fijados habrá que descontar el giro de esta barra en la expresión (2.11),
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 31
rk
rj
ri
Figura 2.10: Equivalencia de esquemas entre un mecanismo con un contacto curva con plano y otro con
pares prismáticos y de rotación.
donde Rb = |O1 C| y ωb = θ̇k . Al igual que en el caso anterior, si esta barra asociada a rj girase,
se le descontarı́a su giro dejando la ecuación (2.13) como
Figura 2.11: Esquema de un biela manivela y los vectores para resolver su cinemática.
El complejo r1 cambia de módulo pero su ángulo es fijo e igual a π/2 y por tanto ω1 = 0.
Los complejos r2 y r3 cambian de ángulo pero no de módulo, lo que lleva a ρ̇2 = ρ̇3 = 0.
Por lo que se obtendrá una ecuación compleja y dos incógnitas donde se sustituye ṙi por la
expresión (2.2)
1+i0 −0.1−i 0.4 −0.2+i 0.2
z}|{ z }| { z }| {
iθ1 iθ1 iθ3 iθ3 iθ2
ṙ1 + ṙ3 = ṙ2 ⇒ ρ̇1 e +ρ1H
ωH
1 ie +@
ρ̇@
3e + ω3 ρ3 ie =@
ρ̇@
2e + ω2 ρ2 ieiθ2
|{z} |{z} |{z}
? ? 10rad/s
El comando de Matlab para resolver este sistema serı́a algo tan sencillo como:
1 >> inv ( [ 1 , = 0.1; 0 , = 0 . 4 ] ) * 1 0 * [ = 0 . 2 ; 0 . 2 ]
2 ans =
3 = 2.5000
4 = 5.0000
La derivación dos veces del lazo de vectores se corresponde con la sustitución de r̈j por la
expresión (2.2). El dato inicial en este análisis es α2 = 0.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 33
ρ@
+(@¨3 − ρ3 ω32 ) eiθ3 + (ρ3 α3 +X
2ρX
˙3 ωX
3 ) ie
iθ3
|{z}
?
¨2 −
ρ@
=(@ ρ2 ω22 ) eiθ2 ρ2X
+ (X X˙ ωX) ieiθ2
2 + 2ρX
αX 2 2
de la que se obtiene
( ) ( ) ( ) ( )
1 −0.1 0.2 −0.4 ρ̈ = −25.625 m/s2
1
ρ̈1 + α3 = −102 + 52 ⇒
0 −0.4 0.2 0.1 α3 = −43.75 rad/s2
aceleración se propuso ún método de solución analı́tica directa. Lo más eficiente para tratar el
problema de posición en la actualidad, es la resolución mediante programas de CAD o incluso
su resolución numérica con programas de cálculo matemático, en ves de recurrir a la engorrosa
trigonometrı́a. Por esta razón en los sucesivo se obviará el problema de posición con geometrı́as
sencillas ya que, como se ha comprobado con el cuatro barras, no hay pérdida de generalidad
para los problemas de velocidad y aceleración.
El retorno rápido que se estudiará con los vectores mostrados en 2.13 tiene las siguientes
√ √ √
dimensiones y datos iniciales: ρ1 = 0.3 m, ρ2 = 0.1 5 m, ρ3 = 0.2 5 m, ρ4 = 0.3 5 m, θ1 =
π/2 rad, θ2 = tan−1 0.5, θ3 = θ4 = tan−1 2, ω2 = 1 rad/s y α2 = 0.
En este punto se pueden aclarar algunos conceptos y advertir de algunos errores a los que
quizás pudo llevar el análisis del cuatro barras. No hay una asociación biunı́voca entre barras y
complejos, que lo pareciese para el 4 barras solo es debido a que todos los pares eran de rotación y
cada barra solo tenı́a dos de ellos. Pero los complejos y las barras son cosas diferentes, las barras
son modelizaciones de los elementos reales conectados por bisagras, guı́as, etc.. Los complejos
son segmentos que unen dos puntos que pueden pertenecer o no a la misma barra. Si por la
elección de puntos de ese segmento, el complejo y la barra giran al unı́sono, entonces la derivada
temporal del argumento coincidirá con la velocidad angular de la barra, pero en cualquier otro
caso no. E igual sucede con la aceleración angular.
r1 + r2 = r3 (2.15)
Debe tenerse en cuenta que, a diferencia del cuatro barras, ahora no todos los términos con
ρ̇i se anulan, ya que ρ̇3 6= 0. Esto es debido a que el vector r3 une el punto Q con un punto de
la deslizadera que no es fijo. El desarrollo de (2.16) serı́a entonces
Figura 2.13: Dimensiones y vectores utilizados para la resolución del mecanismo de retorno rápido.
Para obtener las dos ecuaciones escalares esta vez se recurrirá directamente al método rápido
utilizado con el mecanismo de cuatro barras. Con él se llega a un sistema lineal de dos ecuaciones
con dos incógnitas, que son ρ̇3 y ω3 . En este caso aparece el término ρ̇3 eiθ3 , para escribirlo se
tendrá en cuenta que eiθ3 es un vector unitario cuyas componentes son:
( )
iθ3 1 1
e =√
5 2
√
donde 5 es el módulo del complejo 1 + i 2 que es paralelo a r3 .
La incógnita ω3 es la velocidad angular de la barra de salida, y ρ̇3 es la velocidad a la que el
O O
carro c viaja por la guı́a de b .
( ) ( ) ( )
0.3
ρ̇3 = √
−0.1 ρ̇3 1 −0.4 5
= 0.1342 m/s
1 =√ + ω3 ⇒
0.2 5 2 0.2 ω = 0, 4 rad/s
3
r̈2 = r̈3
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 36
Si se requiriese conocer la aceleración del punto C, extremo e la barra 3, esta será la derivada
segunda de r4 . Debe recordarse que, a diferencia de lo que sucedı́a con r3 , ahora sı́ ρ̇4 = ρ̈4 = 0,
y por tanto
Figura 2.15: Esquema de un mecanismo de actuador lineal: a) Esquema cinemático, b) Vectores para
resolver su cinemática.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 38
Análisis de velocidad
siendo el dato inicial ρ̇2 = 0.03 m/s y las incógnitas ω2 y ω3 . Reordenando la expresión anterior,
esta queda
ρ̇2 eiθ2 = −ω2 ρ2 ieiθ2 + ω3 ρ3 ieiθ3
de la que se obtiene el sistema de ecuaciones lineal del que hallar las incógnitas.
( ) ( ) ( ) √
0.03 2 −0.2 −0.1 ω = 0.9 5 = 0.0503 rad/s
2 40
√ = −ω2 + ω3 ⇒ √
5 1 0.4 −0.2 ω = − 3 5 = −0.1677 rad/s
3 40
Para el cálculo de aceleraciones se derivará dos veces la expresión del lazo de complejos,
teniéndose como dato inicial que ρ̈2 = 0, por lo que
r̈3 = r̈2
(ρ̈2 − ρ2 ω22 ) eiθ2 + (ρ2 α2 +2ρ̇2 ω2 ) ieiθ2 = (ρ̈3 − ρ3 ω32 ) eiθ3 + (ρ3 α3 +
2ρ̇ iθ3
3 ω3 ) ie
| {z } | {z }
α2 ? α3 ?
Figura 2.16: Esquema de mecanismo de leva seguidor: a) Esquema, b) Vectores para resolver su ci-
nemática.
y en general en donde se necesita que una barra siga una determinada función de movimiento
periódica.
En un contacto no conforme entre piezas hay una tangente común en el punto de contacto
(punto C). Esto implica que los vectores que unirı́an los centros de curvatura de la superficie
de contacto con el punto de contacto serı́an colineales y perpendiculares a la tangente de las
superficies en el contacto. En este caso el contacto es entre una superficie con curvatura finita y
otra con curvatura infinita, es decir plana. En este caso no es posible trazar ese vector que una
los centros de curvatura y se utilizará la construcción de la figura 2.16. En ella se utilizan los
O
vectores perpendiculares r3 y r4 . El complejo r3 une el centro de curvatura de arco AB ∈ a con el
O
extremo de r4 . El complejo r4 es paralelo a la recta de contacto DE ∈ b . LA perpendicularidad
O O
ente r3 y r4 es una condición del contacto entre los cuerpos a y b . Ya que de otro mkodo
habrı́a interferencia ente los cuerpos. Debe recordarse que las ecuaciones planteadas serán válidas
mientras el contacto se produzca en el arco AB y en la linea DE.
r2 + r4 = r1 + r3 (2.17)
El complejo r1 es fijo.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 40
Por tanto,
En estos datos hay que tener en cuenta que ρ̇3 no es la velocidad de deslizamiento entre los
O O
sólidos a y b en el punto de contacto C. Su valor serı́a la diferencias de las velocidades del
O O
punto C ∈ a y un punto C 0 ∈ b coincidente con C en el instante de estudio. Es decir:
h ρ4 i h ρ4 i
VC − VC0 = ω2 ρ2 ieiθ2 + ieiθ4 − ω3 ρ3 ieiθ3 − ieiθ4
| {z 2 } | {z 2 }
VC V C0
El análisis de aceleración realizará derivando dos veces (2.17) con las mismas consideraciones
que para la velocidad.
r̈2 + r̈4 = r̈3
se desarrolla como
−ω22 ρ2 eiθ2 − ω32 ρ4 eiθ4 + α3 ρ4 ieiθ4 = ρ̈3 eiθ3 − ω32 ρ2 eiθ3 + α3 ρ3 ieiθ3 + 2ρ̇3 ω3 ieiθ3
y reordenando, se obtiene
−ω22 ρ2 eiθ2 − ω32 ρ4 eiθ4 + ω32 ρ3 eiθ3 − 2ρ̇3 ω3 ieiθ3 = ρ̈3 eiθ3 + α3 ρ3 ieiθ3 − α3 ρ4 ieiθ4
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 41
E1:10
0 0.1 0.2 m
r7
6 r6
D
r5
�2 5
P 2 r1
O 4 r2
3 C
A, B r3 r4
Esta ecuación compleja implica un sistema de dos ecuaciones en las que las dos incógnitas son
ρ̇3 y ω3 . El sistema que se obtiene de esta ecuación es
( ) ( ) ( )
0.1 0 1 ρ̇ = 1 m
3 4 s
2.5 = ω3 + ρ̇3 ⇒
0.2 0.3 0 ω3 = 5 rad 3 s
Al calcular ω3 ahora ya sı́ se conoce el valor de la velocidad de entrada del segundo mecanismo,
el de biela manivela. Para resolverlo se planteará
O
La velocidad del punto D∈ 6 serı́a
√ ( √2 ) ( 1 )
− 2 −3 m
VD = ρ̇7 eiθ7 = 2√
2
= 1
3 − 2 3
s
Con la velocidad angular de la barra 3 se puede calcular la aceleración del punto D del
mecanismo de biela manivela. Para ello se parte del bucle r1 + r6 + r7 = r4 + r5 , que derivado
dos veces da r̈7 = r̈4 + r̈5 . Desarrollando se llega a
En esta ecuación debe tenerse en cuenta que α3 = α4 . Con esta consideración quedarı́an dos
incógnitas; ρ̈7 y α5 .
( √ ) ( ) ( ) ( ) ( )
2 2 2
2√ 5 −0.6 25 0 5 0.4 −0.2
ρ̈7 = 2 − + 2 + α5
− 22 3 0 36 0.6 3 −0.2 0.4
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 43
" √ #( ) ( )
2
0.2 ρ̈7 − 95 ρ̈ = −0.0655 m
2√ 7 s2
= ⇒
− 2
2
0.4 α5 − 35
36
α = −2.5463 rad
5 s2
La aceleración de D será
( √ ) ( )
2
−0.04629 m
AD = ρ̈7 eiθ7 = −0.0655 2√ =
− 2
2
0.04629 s2
actúa sobre la barra que se mueve solidariamente con el vector r2 , la barra de entrada del primer
mecanismo de cuatro barras. Este primer mecanismo de cuatro barras, formado por los vectores
r1 , r2 , r3 y r4 , ha sido sintetizado mediante el método de generación de funciones para otorgar
al punto extremo E un movimiento con velocidad relativamente constante cuando αF = Θ̈F = 0.
La barra de salida del primer mecanismo de cuatro barras es la barra etiquetada como a . Los O
vectores r4 y r6 están unidos rı́gidamente a la barra a . O
El segundo mecanismo de cuatro barras ha sido sintetizado utilizando puntos de precisión
para darle una trayectoria recta horizontal al punto E en la fase de apoyo en tierra de mecanismo
de la pata. Será estudiado mediante el bucle de vectores que forman los vectores r5 , r6 , r7 y r8 .
O
La barra fija del mecanismo es r5 . La barra de entrada es la a y su movimiento es el del vector
O
r6 . El acoplador es la b que gira al unı́sono con el vector r7 . Para sintetizar el mecanismo se
consideró que el ángulo ΘF permanece constante. Esto es ası́ porque durante la fase de apoyo
en tierra ese ángulo, que conduce el tercer mecanismo de cuatro barras, permanecerá estático.
O O
Es decir, durante la fase de apoyo a tierra las barras b y c junto con las barras asociadas a
los vectores r11 y r12 se moverán como un único sólido rı́gido.
El tercer mecanismo de cuatro barras se estudiará con los vectores r10 , r11 , r12 y r13 . Ha
sido sintetizado mediante el uso de puntos de precisión, pero en este caso la trayectoria para
el punto P es vertical. En este caso la barra fija de ese mecanismo es b , que también es el O
acoplador del segundo mecanismo de cuatro barras. La barra de entrada es la que está asociada
al vector r13 . El motor para el movimiento vertical conduce el ángulo ΘF que varı́a el ángulo
entre los vectores r10 y r13 .
11 q1 = 3.0802; q5 = 5.8016;
12 al = 0.9203; be = 2.2703;
13 ga = 0.9606; ep = 0.7101;
14
Figura 2.18: Mecanismos de cuatro barras del movimiento frontal, vertical y lateral de la pata del robot:
a) Mecanismos de cuatro barras para el movimiento frontal y lateral, b) mecanismo de movimiento lateral.
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 46
19 e5 = r 5 * cos ( q5)+ r 8 * cos ( q8)+ r 9 * cos ( q7=be)+ r 1 1 * cos ( q11)+ r 1 4 * cos ( q12+ep)==xE ;
20 e6 = r 5 * sin ( q5)+ r 8 * sin ( q8)+ r 9 * sin ( q7=be)+ r 1 1 * sin ( q11)+ r 1 4 * sin ( q12+ep)==yE ;
21
22 ecs = [ e1 , e2 , e3 , e4 , e5 , e6 ] ;
23 var = [ q6 , q7 , q8 , q11 , q12 , q13 ] ;
24 ini = [ 2 ; 2; 3; 5; 4; 5 ] ;
25 ss = vpasolve ( ecs , var , i n i ) ;
26 q6 = s s . q6 ;
27 q7 = s s . q7 ;
28 q8 = s s . q8 ;
29 q11 = s s . q11 ;
30 q12 = s s . q12 ;
31 q13 = s s . q13 ;
32
33 QV = eval ( q7=ga=q13 ) ;
34 i f QV <0
35 QV = 2 * pi+QV;
36 end
37
r1 + r4 = r2 = r3
CAPÍTULO 2. CINEMÁTICA MEDIANTE NÚMEROS COMPLEJOS 47
que utilizando la ecuación de Euler se podrá dividir en una parte real y otra imaginaria, que son
y del tercero
ρ cos θ + ρ cos θ + ρ cos θ = ρ cos θ
10 10 11 11 12 12 13 13
r10 + r11 + r12 = r13 ⇒ (2.19)
ρ sin θ + ρ sin θ + ρ sin θ = ρ sin θ
10 10 11 11 12 12 13 13
rE = r5 + r8 + r9 + r11 + r14
Pero también, observando los vectores, se pueden establecer una serie de relaciones ente los
ángulos
θ6 = θ4 + α θ9 = θ7 − β θ10 = θ7 − γ
θ13 = θ7 − γ − ΘF θ14 = θ12 + ε
En este momento se puede plantear un sistema de seis ecuaciones con 6 incógnitas. Las
ecuaciones son las etiquetadas como (2.18), (2.19) y (2.20), donde se sustituirán: θ9 = θ7 − β,
θ10 = θ7 − γ, y θ14 = θ12 + ε. Como resultado solo se mantendrán como incógnitas: θ6 , θ7 , θ8 ,
θ11 , θ12 y θ13 . El sistema quedará como
ρ6 cos θ6 = ρ5 cos θ5 + ρ8 cos θ8 + ρ7 cos θ7
θ6
θ7
ρ6 sin θ6 = ρ5 sin θ5 + ρ8 sin θ8 + ρ7 sin θ7
ρ13 cos θ13 = ρ10 cos(θ7 − γ) + ρ11 cos θ11 + ρ12 cos θ12 θ8
⇒
ρ13 sin θ13 = ρ10 sin(θ7 − γ) + ρ11 sin θ11 + ρ12 sin θ12
θ11
xE = ρ5 cos θ5 + ρ8 cos θ8 + ρ9 cos(θ7 − β) + ρ11 cos θ11 + ρ14 cos(θ12 + ε) θ12
yE = ρ5 sin θ5 + ρ8 sin θ8 + ρ9 sin(θ7 − β) + ρ11 sin θ11 + ρ14 sin(θ12 + ε) θ13
Una vez resuelto el sistema, el ángulo de movimiento vertical se obtiene inmediatamente con
la diferencia
ΘF = θ7 − γ − θ13
calcular la TCD con las expresiones anteriores, pero para ofrecer una alternativa, se resolverá
numéricamente mediante Matlab. En este caso el código será
1 function [ xE , yE ] = TCD 2D(QF,QV)
2 syms q2 q3 q4 q7 q8 q11 q12 q13 xE yE ;
3
10 q1 = 3.0802; q5 = 5.8016;
11 al = 0.9203; be = 2.2703;
12 ga = 0.9606; ep = 0.7101;
13
14 q2 = QF;
15
34 xE = r 5 * cos ( q5)+ r 8 * cos ( q8)+ r 9 * cos ( q7=be)+ r 1 1 * cos ( q11)+ r 1 4 * cos ( q12+ep ) ;
35 yE = r 5 * sin ( q5)+ r 8 * sin ( q8)+ r 9 * sin ( q7=be)+ r 1 1 * sin ( q11)+ r 1 4 * sin ( q12+ep ) ;
36 end
El cálculo gráfico ha sido muy utilizado desde siempre en la ingenierı́a, sobre todo cuando no
existı́an las calculadoras ni los ordenadores. Se realizaba con papel, lápiz, escuadra y cartabón.
Era un medio rápido de hacer cálculos estructurales o cinemáticos y dinámicos, aunque tenı́a el
problema de una precisión deficiente. A mediados del siglo XX con los primeros ordenadores,
el cálculo gráfico se vio relegado por los análisis analı́ticos. Pero en las últimas décadas, con la
aparición de potentes programas de CAD (Computer Aid Design), el cálculo gráfico resuelve su
gran problema de precisión manteniendo su versatilidad e inmediatez a la hora de hacer análisis
en la mecánica [7].
En este capı́tulo se expondrán las herramientas que permitirán el cálculo gráfico de las
velocidades de un mecanismo. Para ello se dibujarán diagramas con las velocidades, en la que
las velocidades de los puntos del mecanismo están representadas por segmentos orientados (los
vectores). A la figura que forman los extremos de los vectores velocidad de los diferentes puntos de
una barra, trazados desde un determinado origen, lo denominaremos cinema de velocidades,
siendo el polo del cinema ese origen. En estos apuntes usualmente usaremos la letra Q para
notar al polo del cinema de velocidades.
Para la construcción de los cinemas de velocidades se utilizarán expresiones vectoriales que
ligarán las velocidades de puntos de las barras conectadas mediante pares cinemáticos. La opera-
tiva será similar a la realizada con complejos, es decir, se plantearán ecuaciones de lazo a través
de los puntos del mecanismo. La diferencia estriba en que esta vez la resolución será gráfica, en
vez de con un sistema de ecuaciones e incógnitas.
Las expresiones que se usarán en el cálculo gráfico serán vectoriales. Aquı́ estos apuntes
difieren de lo que normalmente se puede encontrarse en la bibliografı́a de Teorı́a de Máquinas
como, por ejemplo, [8], [5] y [6]. Lo usual es construir estas expresiones usando calculo vec-
torial mediante el método de las velocidades relativas e introduciendo diferentes observadores
con movimiento relativo entre ellos. De esta forma se explica este método en su formulación
tridimensional, para luego particularizarse para los mecanismos planos. Pero, puesto que el 99 %
de los mecanismos que se diseñan y el 99,99 % de los fabricados son planos, en estos apuntes se
llegará a las expresiones directamente para el caso plano a partir de las expresiones obtenidas
con el análisis mediante complejos. También se utilizará un solo tipo de observador intermedio
que es el ligado a una barra con un par prismático que la une a otra. Si se requiriese analizar
50
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 51
un mecanismo tridimensional, se puede demostrar que las expresiones obtenidas son igualmente
válidas, encontrándose la demostración en bibliografı́a como la referenciada en este párrafo.
Al igual que se hizo con cálculo cinemático mediante complejos, el método se desarrollará
usando ejemplos. Estos ejemplos tendrán geometrı́as sencillas, con distancias fáciles de calcular
y con ángulos para las barras con valores tales como: π/2, pi/6, arctan(1/4), y similares. Se
hará ası́ por una cuestión didáctica, para que ası́ se aprecien mejor las relaciones geométricas
de paralelismo y perpendicularidad. Pero estas geometrı́as no implican pérdidas de generalidad.
Para los métodos gráficos da igual unas distancias o ángulos que otros, y es con estos programas
modernos de CAD como realmente adquieren efectividad y potencia los métodos gráficos, espe-
cialmente con los programas de CAD paramétrico, tales como SolidWorks, Catia, CREO CAD,
etc.
El complejo eiθ tiene módulo unidad y dirección θ, por lo que ei(θ+pi/2) es un complejo
también de módulo 1 y girado π/2 respecto al anterior.
En la figura 3.1 también se muestran los vectores unitarios eiθ y ei(θ+π/2) , ası́ como las dos
componentes perpendiculares del complejo ṙ, que son ρ̇ eiθ y θ̇ρ ei(θ+π/2) . En la figura 3.1 b) se
muestra, por tanto, la ecuación de lazo (3.1).
En este momento se abren dos posibles casos, que se muestran en la figuras 3.2 y 3.3:
Esta ecuación anterior se puede escribir gráficamente con el lazo de tres vectores que se
muestra en la figura 3.2 b) con origen en Q. Esta expresión indica que las velocidad de un punto
B perteneciente al mismo sólido que otro punto A es igual a la velocidad de A más el vector
ωa × r, que es perpendicular a r y con módulo igual al producto de la velocidad angular del
sólido por la distancia entre A y B.
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 53
Figura 3.3: Relación entre velocidades de dos puntos coincidentes pero pertenecientes a diferentes sólidos.
Caso b). Dos puntos coincidente pertenecientes a dos barras unidas por un par
prismático.
Este caso se muestra en la figura 3.3, y en él se establece la relación entre las velocidades
VA y VB de dos puntos coincidentes momentáneamente en el espacio, pero pertenecientes a dos
O O
sólidos diferentes. El punto A ∈ a y el punto B ∈ b . Estas dos barras están unidas mediante
un par prismático, mediante uno de pin-ranura, u otro que obligue al desplazamiento en linea
recta entre los puntos A y B. Esta dirección de movimiento relativo está marcada con la linea
de trazos en la figura 3.3a) y tiene la dirección θ. Como el segmento r tiene longitud 0 (ρ = 0)
ya que coinciden los puntos A y B en el instante de estudio, se tendrá que la expresión 3.1 pasa
a
VB = VA + ρ̇ eiθ + θ̇ρe
A
i(θ+π/2)
⇒ VB = VA + vB |a (3.3)
La expresión 3.3 conforma el lazo de los tres vectores que se muestra en la figura 3.3 b)
con origen en Q. Esta expresión indica que las velocidad de un punto B coincidente con A es
igual a la velocidad de A más vB |a , donde vB |a es un vector con la dirección del movimiento
O O
relativo entre a y b y módulo la variación de la distancia entre A y B. Este vector vB |a es la
O O
velocidad con la que el observador A ∈ a , considerando a la barra a como fija, ve moverse al
O
punto B ∈ b . Es decir, puede entenderse como la velocidad de B pero habiendo realizado una
O
inversión en el mecanismo de manera que la barra a pasa a ser la barra fija.
De hecho, reservaremos esta notación con minúscula vP |j para notar la velocidad relativa
O
con la que un observador en la barra j ve moverse al punto P. Si se sabe la barra a la cual nos
referimos, la notación podrá simplificarse de manera que vP |j = vP , bastando con el uso de la
minúscula para marcar que es una velocidad relativa a una barra con un par prismático.
O O O
barras será e para la barra de entrada, a para el acoplador, s para la barra de salida y f O
para la barra fija.
Al igual que en de la resolución mediante complejos, el paso previo obliga a que se haya
realizado el cálculo de posición. Para ello, en este ejemplo se conocen las longitudes de la barra
fija, la de entrada, la del acoplador y la de salida, que son: lf , le , la y ls . En estos apuntes usaremos
la notación li para notar una longitud geométrica que no cambia durante el funcionamiento de
la máquina.
Otro dato necesario será el ángulo inicial de alguna barra, a esta barra la denominaremos
barra de entrada y al ángulo en cuestión θe . Las incógnitas del problema de posición serı́an θa y
θs que son los ángulos de la barra acoplador y la barra de salida. En este ejemplo en concreto,
dada su sencillez geométrica, es inmediato. Si se utilizase un ordenador con un programa de
CAD, la resolución del problema de solución serı́a igualmente sencilla para cualquier geometrı́a,
y bastarı́a con dibujar el mecanismo para tenerlo resuelto, ofreciendo el programa todas las
posiciones y ángulos que se requiriesen.
En estos apuntes, puesto que se dibujan los esquemas y cinemas, es necesario incluir una
escala para saber las longitudes originales del mecanismo que se está estudiando, ası́ como las
cantidades reales de metros por segundo de las velocidades de la resolución gráfica del problema
de velocidad. Y dependiendo de las velocidades, será habitual que las escalas para los esquemas
y los cinemas de velocidad y aceleración no sean las mismas. Pero si se estuviera resolviendo la
cinemática en un ordenador con software CAD, esto ya no serı́a necesario gracias a la facilidad
para hacer zoom en las diferentes partes del dibujo y los cinemas sin pérdida de precisión.
En la figura 3.4 se muestra el mecanismo a resolver en la posición de estudio y con el montaje
correcto. Los datos iniciales toman los valores:
√
lf = 0.4 m le = 0, 2 m la = 0.2 2 m
√
ls = 0.2 5 m θf = 0 θe = π/2 rad
y del análisis de posición se obtendrı́a los ángulos restantes, cuyo cálculo es obvio para esta
sencilla geometrı́a. En este caso esos ángulos son:
π
θa = rad
4
π 1
θs = + tan−1 rad
2 2
Como dato inicial para resolver este problema de velocidades se deberá conocer la velocidad
de la barra de entrada ωe = 1 rad/s. Dada ωe , el cálculo gráfico de las velocidades consiste en
obtener la velocidad angular de la barra de salida s . O
Las ecuaciones vectoriales que permiten este cálculo derivan de la expresión ??. Vectorial-
mente esta resolución consiste en el cálculo de la velocidad del punto B, de la que se desconoce
tanto módulo como argumento. Relacionando a los puntos A y B como pertenecientes ambos al
O
acoplador a mediante la ecuación 3.2, se puede escribir:
VB = VA + ωa × rAB
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 55
Figura 3.4: Esquema de cuatro barras y vectores usados para resolver su cinemática.
donde la velocidad VA se puede calcular teniendo considerando que tanto O como A son puntos
O
que pertenecen a la barra de entrada e y por tanto
VA = H
VHO + ωe × rOA
teniendo en cuenta que O, es un punto fijo del chasis del mecanismo y por tanto su velocidad es
nula. Por otro lado, desde la barra de salida también se puede escribir
VB = H
VHP + ωs × rPB
P + ωs × rPB = VH
VH
H H O + ωe × rOA + ωa × rAB
| {z } | {z } | {z }
módulo? móduloX módulo?
direcciónX direcciónX direcciónX
La ecuación compleja anterior tiene dos componentes escalares, y solo dos incógnitas que
son ωa y ωs . Los vectores que aparecen en la ecuación se puede proyectar en unos ejes para ası́
poder realizar los productos vectoriales. Es decir, se escribirı́a
i j k i j k i j k
0 0 ωs = 0 0 ωe + 0 0 ωa
ls cos θs ls sin θs 0 le cos θe le sin θe 0 la cos θa la sin θa 0
El desarrollo de esta ecuación vectorial da dos ecuaciones escalares significativas donde ωe
es la velocidad angular de la barra de entrada, que es el dato de partida. Por tanto, las únicas
dos incógnitas serán ωa y ωs . Por tanto, se obtiene un sistema lineal de 2 ecuaciones con 2
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 56
con las incógnitas ωa y ωs y donde ωe es dato. Este sistema coincide con el obtenido en la sección
2.3
Pero la ecuación compleja que de la que se obtiene el sistema también se puede resolver
gráficamente. La resolución gráfica de esta ecuación compleja se muestra en la figura 3.5 en
cuatro pasos.
PASO 1:
El primer paso consiste en dibujar la velocidad del punto A. Esta es muy fácil de obtener
O
puesto que A pertenece a la barra de entrada. Ası́, puesto que O y A pertenecen a a , se puede
usar la expresión 3.2, que toma la forma
VA = H
VHO + ωe × rOA (3.4)
O O
donde VO = 0 puesto que O además de a a también pertenece a la barra fija f . De (3.4) se
infiere que la velocidad de A es un vector perpendicular a rOA , siendo su módulo ωe le = 0.2 m/s.
En el cinema de velocidades, el extremo de la velocidad del punto A se representará con esa misma
letra pero en minúscula (a). Ası́ el segmento Qa se dibujará con la longitud de 0.2 m. Si el papel
donde se dibuja el cinema no lo permitiese, este vector se dibujarı́a escalado. Afortunadamente,
en un programa de CAD, ese problema con las escalas se puede obviar, ya que la precisión se
mantiene independientemente del tamaño del dibujo.
PASO 2:
VB = VA + ωa × rAB (3.5)
PASO 3:
O
Volviendo a utilizar la ecuación 3.2, esta vez para la barra de salida s , se obtiene
VB = H
VHP + ωs × rPB
Figura 3.5: Pasos para resolución gráfica del problema de velocidad en un mecanismo de 4 barras.
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 58
PASO 4:
Pero si el extremo del vector VB debe estar en algún punto de la recta perpendicular a rAB
y también en alguno de la perpendicular a rPB , solo podrán cumplirse ambas condiciones si la
velocidad VB tiene su extremo en la intersección de ambas, tal y como se muestra en el PASO
4 de la figura 3.5. El módulo de la velocidad de B se obtiene de la medición en el gráfico, ya sea
en el papel o en el ordenador. Conocido el módulo de la velocidad de B la velocidad angular de
la barra de salida se obtiene de dividirlo entre la distancia PB.
0.1491
|VB | = 0.1491 m/s ⇒ ωs = √ = 0.333 rad/s
0, 2 5
0.0943
|ωa × rAB | = 0.0943 m/s ⇒ ωa = √ = 0.333 rad/s
0.2 2
Debe tenerse en cuenta que con las divisiones anteriores se obtienen los módulos de las
velocidades angulares. En el ejemplo se obtiene que las velocidades angulares del acoplador y
de la barra de salida son iguales, pero el sentido es inverso, siendo horario para el acoplador
y antihorario para la barra de salida. Estos sentidos son fáciles de determinar simplemente
observando las direcciones de los vectores. Ası́ para determinar el sentido de giro del acoplador
advertiremos que rAB apunta arriba a la derecha mientras ωa × rAB apunta hacia abajo a la
derecha, tal y como se muestra en el PASO 4 de la figura 3.5. Análogamente comprobarı́amos
el sentido de la velocidad angular de la barra de salida mirando las direcciones de rPB y de
ωs × rPB .
Una vez que se obtiene el cinema de velocidades, este no deja de ser un simple triángulo.
Ası́ que si no se dispone de un ordenador con un programa de CAD y se quiere obtener más
precisión que con lápiz y regla, siempre se puede resolver trigonométricamente el triángulo del
cinema de velocidades.
VC = VA + ωa × rAC = VB + ωa × rBC
|{z} |{z} | {z } |{z} | {z }
módulo? móduloX módulo? móduloX módulo?
dirección? direcciónX direcciónX direcciónX direcciónX
donde los dos últimos miembros representan un sistema de dos ecuaciones con dos incógnitas.
El proceso se muestra con la figura 3.6.
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 59
Figura 3.6: Pasos para determinar la velocidad de un tercer punto de una barra.
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 60
Figura 3.7: Evolución del acoplador y determinación del sentido mediante el cinema.
PASO 1:
Se parte de las velocidades de dos puntos, en este caso VA y VB . También debe ser conocida
la posición del punto en cuestión del cual se quiere conocer la velocidad. En este caso se pretende
averiguar VC , estando C situado donde se indica en la figura 3.6.
PASO 2:
VC = VA + ωa × rAC
PASO 3:
VC = VB + ωa × rAC
PASO 4:
lo hace hacia la derecha y arriba, lo que muestra el sentido de giro horario del acoplador cuando
la barra de entrada gira en sentido antihorario.
VB = VA + ωb × rAB
|{z} |{z} | {z }
módulo? móduloX módulo?
direcciónX direcciónX direcciónX
Paso 1:
Paso 2:
Usando la expresión 3.2 se obtiene una recta donde ubicar el final del vector VA , ya que
VB = VA + ωb × rAB
Paso 3:
Por otro lado el pistón tiene movimiento horizontal, y por tanto, la VB también será hori-
zontal. De este modo se obtiene otra recta donde ubicar el final del vector VB
Paso 4:
Figura 3.10: a) Evolución del mecanismo de retorno rápido. b) Evolución de la inversión del mecanismo
O
anterior fijando c
VB = VA + vB |c
VB = H
VHP + ω b × rPB (3.6)
e igualando las expresiones anteriores se llega a la ecuación vectorial (que representa dos ecu-
ciones escalares)
VA + vB |c = ωb × rPB (3.7)
|{z} |{z} | {z }
móduloX módulo? módulo?
direcciónX direcciónX direcciónX
de donde se obtienen las incógnitas, que son los módulos de vB |c y ωb × rP B . Ası́, una vez que
se ha medido el módulo de ωb × rPB , el módulo de ωb se obtiene dividiendo por rP B .
La resolución gráfica se muestra en los siguientes pasos de la figura 3.9
Paso 1:
En este primer paso se dibuja la velocidad VA aplicando la expresión (3.2), a los puntos A
√
y O. esta velocidad tiene módulo 20 5 m/s y dirección perpendicular a rOA .
VA = H
VHO + ωa × rOA
Paso 2:
VB = H
VHP + ωb × rPB
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 66
Paso 3:
En este paso, a diferencia de los anteriores, se usará la expresión (3.3) que liga las velocidades
de dos puntos coincidentes en el instante de cálculo pero pertenecientes a diferentes barras. En
O O
concreto, los puntos serán el A ∈ c y el B ∈ b .
VB = VA + vB |c
donde vB |c es, como ya se ha comentado, la velocidad con la que el punto A ve moverse al punto
O
B, pero considerando a la barra c como fija. Es decir realizando la inversión pertinente.
La velocidad vB |c , como ya hemos visto, tendrá la dirección del par prismático, por lo que
la velocidad VB tendrá su fin en la lı́nea de puntos paralela al par prismático en el punto de
estudio, tal y como se ha dibujado en la figura 3.9 Paso 3.
Paso 4:
O
Conocida la velocidad de B, la velocidad angular de la barra b serı́a
VB
ωb =
rP B
VB = VA + vB |c
VB = H
VHP + ωc × rPB
donde
VA + vAC = ωc × rPC
|{z} |{z} | {z }
módulo? móduloX módulo?
direcciónX direcciónX direcciónX
Paso 1:
En este primer paso se dibuja la velocidad dato del problema, que es la velocidad relativa
O O
vB |a = 30 mm/s. Esta es paralela a la linea de movimiento de b respecto de a en la situación
de estudio. Dada la baja velocidad del mecanismo los cinemas de velocidades se podrán dibujarse
a escala 1:1.
Paso 2:
O
Posteriormente se traza la dirección que tendrá la velocidad de B ∈ c , que es perpendicular
a rPC ya que
VC = H VHP + ωc × rPC
Paso 3:
VB = vB |a + VA
Paso 4:
En este paso se determina VB como el vector con origen en Q y fin en las dos rectas trazadas
O
en los pasos 2 y 3. Si además se quisiera determinar la velocidad de un punto D ∈ c se calcuları́a
conforme a lo explicado en la sección 3.3. Según este método , la velocidad vD tendrı́a su extremo
donde se cortan las perpendiculares a los segmentos AD y P D, tal y como se muestra en el paso
4 de la figura 3.11.
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 69
VB = H
VHP + ωb × rPB
VB = VA + vB |a
O
donde VA es conocida por ser A ∈ a perteneciente a la barra de entrada. La dirección de vB |a ,
O
que es la velocidad con la que A∈ a observa el movimiento de B, también estás determinada,
al ser paralela a la parte plana del seguidor. Por esto
ωb × rPB = VA + vB |a
| {z } |{z} | {z }
módulo? móduloX módulo?
direcciónX direcciónX direcciónX
Paso 1:
Con la expresión 3.3 se obtiene una recta en cuyo final debe estar VB ya que
VB = VA + vB |a
Paso 3:
También, con la expresión 3.2 se obtiene una segunda recta para determinar VB , debido a
O
que se sabe la velocidad nula del punto P ∈ b , al ser la articulación con la barra fija de la barra
Ob.
VB = H
VHP + ωb × rPB
Paso 4:
EJEMPLO 1
En un mecanismo de cuatro barras es sencillo colocar los CIR: I21 , I32 , I43 e I41 , por lo que
quedan por determinar I31 e I42 . Del teorema de los tres centros se sabe que I31 debe de estar
en la misma recta que I21 e I32 , también debe estar en la recta que une I43 e I41 , por lo que la
posición de I31 se encintrará en la interescción de las dos rectas anteriores.
Análogamente, I42 se encontrará la la intersección de las rectas I32 I43 e I21 I41
EJEMPLO 2
En este mecanismo de biela-manivela las posiciones de I21 , I32 e I43 son obvias. Entre las
O O
barras 4 y 1 existe un par prismático de dirección vertical. Para entender donde estarı́a
el CIR I41 se puede recurrir al mecanismo de cuatro barras de cuya degeneración deviene el
biela manivela. En este el par prismático del piston es sustituido por una barra horizontal e
CAPÍTULO 3. CÁLCULO GRÁFICO DE VELOCIDADES 73
infinitamente larga con un extremo en I43 y el otro en el infinito en dirección horizontal. Allı́
estarı́a el I41 . Esta información del I41 es suficiente para establecer que mediante el toerema
de Kennedy, el CIR I31 debe de estar en algun punto de la recta horizontal que pasa por I43 .
También al conocerse los emplazamientos de I21 e I32 , se sabe que I31 debe estar en la recta que
los une. Ası́ el I31 estará en la interesección de las dos rectas anteriores.
Para determinar I42 se actuará de forma similar, ya que, por estar I41 muy, muy alejado,
y en dirección horizontal, la recta horizontal que pasa por I21 también servirá para usando el
teoremo de los tres centros, interseccionarla con la que contiene a I43 e I32 para ası́ determinar
I42 .
EJEMPLO 3
O O
En este mecanismo de retorno rápido el par prismático que une a 3 y 4 ahora no es
vertical, está inclinado, y pr tanto el i43 habrá que buscarlo también en el infinito, pero ahora
no en dirección horizontal si no en la dirección perpendicular al par prismático. tal y como se
muestra en la figura 3.14. Por tanto, las rectas que se usarán en el teorema de Kennedy tendrán
O
dirección perpendicular a 3 y pasarán por los puntos I42 e I31 . Ası́, I32 se encontrará en la
O
interescción de la recta que pasa por I4 2 y es perpendicular a la barra 3 , y la recta qeu pasa
por I2 1 e I31 . Análogamente, I41 está en la interesección de I42 I21 con la recta perpendicular a
O3 que pasa por I31 .
EJEMPLO 4
PASO 1
PASO 2
Se calculan los CIR del mecanismo como se indicó en la sección 3.10, concretamente será
necesario localizar I21 , I41 e I42 .
PASO 3
O
Se calcula la velocidad de I42 ∈ 2 respecto a la barra fija.
VI42 = ω2 × O1 I42
PASO 4
O
Se calcula la velocidad de I42 ∈ 4 respecto a la barra fija.
VI42 = ω4 × O2 I42
y por tanto
|O1 I42 |
ω4 = ω2
|O2 I42 |
Capı́tulo 4
En este capı́tulo se expondrán las herramientas que permitirán el cálculo gráfico de las
aceleraciones de un mecanismo. Para ello se dibujarán diagramas con las aceleraciones, en la que
éstas están representadas por segmentos orientados (los vectores). A este conjunto de vectores
que representan las aceleraciones lo llamaremos cinemas de aceleraciones.
Como paso previo al análisis de aceleraciones se habrá de realizar el análisis de posición y
velocidad. Estos serán imprescindibles para poder ejecutar el análisis de aceleraciones.
77
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 78
Figura 4.2: Vectores que relacionan la aceleración de dos puntos pertenecientes a una misma barra.
De la fórmula de Euler se desprende que el complejo eiθ tiene módulo unidad y dirección θ,
por lo que ei(θ+pi/2) es un complejo también de módulo 1 y girado π/2 respecto al anterior. Por
tanto la expresión anterior pasa a
O
Este caso se muestra en la figura 4.2. Considérese dos puntos A, B ∈ a , diferentes, unidos
por rAB , con |rAB | = ρAB > 0, siendo ρAB constante. Esto implicará que ρ̇AB = ρ̈AB = 0 y la
ecuación (4.1) pasará a tomar la forma
donde ωa = θ̇ y αa = θ̈.
El término −ρAB ωa2 eiθ tiene de módulo ωa2 ρAB y misma dirección pero sentido contrario a
rAB lo que hace que se pueda escribir como el siguiente doble producto vectorial
O
Además, al estar el complejo rAB rı́gidamente unido al solido a al que pertenecen A y B,
θ̈AB será igual a la aceleración angular de la barra a la que llamaremos αa .
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 79
Figura 4.3: Vectores que relacionan la aceleración de dos puntos pertenecientes a distintas barras pero
coincidentes.
Caso b). Dos puntos pertenecientes a distintas barras pero coincidentes en el espacio
Este caso se muestra en la figura 4.3, y en él se establece la relación entre las aceleraciones
AA y AB de dos puntos coincidentes momentáneamente en el espacio, pero pertenecientes a
O O
dos sólidos diferentes: A ∈ a y B ∈ b . Las dos barras en cuestión están unidas mediante un
par prismático, mediante uno de pin-ranura, u otro que obligue a que el desplazamiento entre
los puntos A y B tengan una determinada dirección rectilı́nea. Esta dirección de movimiento
relativo está marcada con la linea de trazos en la figura 4.2 b) y tiene la dirección θ. Como el
segmento rAB tiene longitud 0 ya que coinciden los puntos A y B en el instante de estudio, se
tendrá que la expresión 4.1 pasa a
2 iθ
AB = AA + (ρ̈ − Z
ρθ̇Z) e + (ρθ̈ + 2ρ̇θ̇) ei(θ+pi/2)
SS
O O
El término ρ̈ eiθ es la aceleración de B ∈ b considerando fija a la barra a . Este término tiene
O O
dirección θ, como muestra el estar multiplicado por eiθ y que a y b están unidas por el par
O
prismático de dirección θ. A ρ̈ eiθ se le llamará aB |a , la aceleración relativa de B ∈ b respecto
O
a A ∈ a . Como ya se hizo para la velocidad relativa, usaremos la letra minúscula aP |j para
remarcar que es la aceleración del punto P, pero que se desde el sólido j. Esto es, la aceleración
que tendrı́a P si se hiciera una inversión de manera que la barra j pasase a ser la fija.
El otro término que no se cancela es 2ρ̇θ̇ ei(θ+pi/2) , que es perpendicular a la dirección del par
prismático y coincidente con el productor vectorial 2ωa × vB |a . Por tanto la ecuación anterior
pasa a escribirse como
AB = AA + aB |a + 2ωa × vB |a (4.3)
AA = H
AHO + ωe × (ωe × rOA ) + αe × rOA
AB = H
AHP + ωs × (ωs × rPB ) + αs × rPB
e igualando las expresiones anteriores se tiene un sistema vectorial plano (dos ecuaciones) donde
se tienen dos incógnitas, que son los dos módulos marcados con el sı́mbolo de interrogación.
Paso 1
En este primer paso de la figura 4.4 se dispone de las soluciones gráficas al problema de
posición y velocidad.
Paso 2
O
Como dato inicial se debe conocer la aceleración angular de la barra e . Conocida este dato
inicial, la aceleración del punto A se obtendrá con la expresión 4.2:
AA = H
AHO + ωe × (ωe × rOA ) + αe × rOA
donde todo es conocido en la expresión anterior, por lo que la AA está totalmente determinada.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 81
Paso 3
El siguiente paso es conocer la velocidad del punto B, para lo cual se usará de nuevo la
expresión 4.3. En este caso
en donde el termino ωa ×(ωa ×rAB ) es totalmente conocido, ya que de él se sabe que su dirección
es la rAB , sentido contrario a rAB y módulo ωa2 ρAB . En este paso se dibuja este término.
Paso 4
Pero todavı́a se desconoce el valor de la aceleración angular del acoplador αa por lo que del
termino αa ×rAB solo se conoce la dirección, que es lo que representa la linea de trazos dibujada
en este paso.
Paso 5
AB = H
AHP + ω s × (ω s × rPB ) + αs × rPB
Paso 6
Pero del término αs × rPB solo se conoce la dirección, hecho que se refleja en el dibujo de
la linea de puntos que se coloca en este paso.
Paso 7
PASO 1
Como dato inicial se tiene la aceleración de los puntos A y B de la barra triangular con
vértices en los puntos ABC.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 83
Figura 4.5: Pasos para el cálculo gráfico de aceleraciones conocida la de otros dos puntos de una barra.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 84
PASO 2
AB = AA + ω × (ω × rAB ) + α × rAB
Por tanto, se puede trazar una lı́nea paralela a rAB y otra perpendicular de manera que ese
par de lineas arranque en el punto final del vector AA y termine en el punto final de AA . Ası́,
ese bucle de vectores reproducirı́a la ecuación anterior.
Conocido ω × ω × rAB y α × rAB se obtiene fácilmente tanto la aceleración α como la
velocidad angular de la barra ω, que serán
s
|ω × ω × rAB |
ω=
rAB
|α × rAB |
α=
rAB
PASO 3
La expresión de la aceleración de C es
AC = AA + ω × (ω × rAC ) + α × rAC
donde se puede trazar el termino ω × (ω × rAC ) dado que es paralelo a rAB y se conocen los
módulos de los vectores de la expresión.
PASO 4
4.3.2. Método 2
El triángulo que forman los extremos de las aceleraciones en el cinema es semejante al que
forman los puntos de la barra. Este triángulo este escalado y también esté girado un cierto
ángulo, que a diferencia del cinema de velocidades, no tiene por qué ser pi/2.
La demostración se obtiene de las expresiones que ligan la aceleración del punto C con la de
los puntos A y B.
AC = AA + ω × (ω × rAC ) + α × rAC
AC = AB + ω × (ω × rBC ) + α × rBC
AC − AA = ac = ω × (ω × rAC ) + α × rAC
AC − AB = bc = ω × (ω × rBC ) + α × rBC
AA − AB = ba = ω × (ω × rBA ) + α × rBA
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 85
4.3.3. Método 3
Un tercer método paa el cálculo de la aceleración del punto C conocidas las de A y B se
muestra en la figura 4.6.
PASO 1
Como dato inicial se tiene la aceleración de los puntos A y B de la barra triangular con
vértices en los puntos ABC.
PASO 2
AC = AA + ω × (ω × rAC ) + α × rAC
PASO 3
Del término α×rAC se conoce su dirección, que es perpendicular a rAB , pero es una incógnita
su módulo, por tanto solo podemos trazar una lı́nea de puntos de manera que en algún punto
de ella termine la aceleración AA .
PASO 4
PASO 5
En este paso se dibuja otra lı́nea de trazos donde también deberá terminar AC ya que
tampoco se conoce la dirección de α × rBC .
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 86
Figura 4.6: Pasos para el cálculo gráfico de aceleraciones conocida la de otros dos puntos de una barra.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 87
PASO 6
AB = AA + aB |c + 2ωc × vB |c
AA = H
AHO + ωa × (ωa × rOA ) + αa ×
hhhrh
OA
h
AB = H
AHP + ωb × (ωb × rPB ) + αb × rPB
e igualando las expresiones anteriores se tiene un sistema vectorial plano (dos ecuaciones) donde
se tienen dos incógnitas, que son los dos módulos marcados en
PASO 1
Figura 4.7: Pasos para el cálculo gráfico de aceleraciones de un mecanismo de retorno rápido.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 89
PASO 2
O
La aceleración de A ∈ a es conocida e igual a AA = ωa × (ωa × rOA ), por lo que se dibujará
un segmento de longitud ωa2 rOA = 12 ×0.2236 m/s2 , dirección paralela a rOA y sentido contrario
al de rOA .
PASO 3
PASO 4
Se desconoce el módulo de aB |c , pero se sabe su dirección, que por ser una aceleración relativa
O O
es paralela a la velocidad relativa y por tanto paralela a la guı́a que une las barras a y s . De
este modo se dibuja una linea de puntos con dirección paralela a la guı́a. La aceleración de B
terminará en algún punto de esa linea, teniendo su origen en el cinema de aceleraciones.
PASO 5
PASO 6
PASO 7
la posición velocidad y aceleración relativa con la que se mueven las dos piezas principales que
O O
constituyen el actuador, que en la figura están etiquetadas como a y b .
En el PASO 1 se muestra el cinema de velocidades. En la sección 3.6 se resolvió la velocidad
del mecanismo para una velocidad de entrada de vAB = 30 mm/s. Como resultado se obtuvo
que VA = 22.5 mm/s y VB = VC = 37.5 mm/s, lo que lleva a ωa = ωb = VA /rOA = 22.5/447.2 =
0.05 rad/s y ωc = 37.5/223.6 = 0.168 rad/s.
Y para el cálculo de aceleraciones el dato inicial será que la velocidad relativa del actuador
es constante y por tanto aAB = 0.
O O O
Para resolver el problema se considerarán tres puntos A ∈ a , B ∈ b y C ∈ c . Ası́, usando la
expresión 4.3 que relaciona las aceleraciones de dos puntos de diferentes barras pero coincidentes
en el instante de cálculo, y teniendo en cuenta el dato inicial de v̇AB = 0 se tendrá
AB = AA + X X + 2ω × vAB
aAB
X
siendo la aceleración de A
AA = H
AHO + ωa × (ωa × rOA ) + αa × rOA
AC = H
AHP + ωc × (ωc × rPC ) + αc × rPC
Con las expresiones anteriores se podrá montar la siguiente ecuación, donde solo queda por
conocer las dos incognitas marcadas.
PASO 1
PASO 2
PASO 3
PASO 4
Del término αa × rOA solo se conoce la dirección. Por ello solo se puede trazar la lı́nea de
puntos.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 91
Figura 4.8: Pasos para el cálculo gráfico de aceleraciones de un mecanismo de retorno rápido.
CAPÍTULO 4. CÁLCULO GRÁFICO DE ACELERACIONES 92
PASO 5
PASO 6
Del término αc × rPC solo se conoce la dirección. Por ello solo se puede trazar la lı́nea de
puntos.
PASO 7
Dinámica
Hasta ahora se han analizado los mecanismos solo en lo referente a su movimiento, sin
importar la causa que hacı́a que las barras girasen y/o se desplazasen con un determinada
velocidad o aceleración. Pero las máquinas solo se mueven si en el algún momento se les imprimió
fuerzas o pares. Además, solo mantendrán su movimiento si se continúa aplicando fuerzas, salvo
en el hipotético caso de que no exista rozamiento. En este capı́tulo se verá la relación entre
el movimiento de las máquinas y las fuerzas externas que lo generan. A esta relación se la
denominará análisis dinámico, y puede ser de dos tipos: directo e inverso.
El paso previo al análisis dinámico, en cualquiera de sus tipos, es el análisis cinemático. Éste
es imprescindible ya que relaciona las velocidades de las diferentes barras del mecanismo.
Para un sólido rı́gido en el que se puede considerar la masa como continua, el sumatorio pasa
a ser un integral, siendo entonces la expresión a usar
Z
1
rG = ρ r dV
m V
93
CAPÍTULO 5. DINÁMICA 94
y por tanto la cantidad de movimiento del sistema o sólido es igual a la cantidad de movimiento
que tendrı́a una única particula situada en G, moviéndose como este punto y que tuviera una
masa igual a la total del sistema, es decir:
p = m VG
T = Ḣ
donde T es la suma de los momentos aplicados sobre el cuerpo más los momentos de las fuerzas
exteriores, considerando que los momentos se toman respecto al centro de gravedad o un punto
fijo.
El momento cinético de un sólido rı́gido se escribe:
HG = IG ω
CAPÍTULO 5. DINÁMICA 95
donde IG es el tensor de inercia del cuerpo y ω es la velocidad angular. Esta ecuación tensorial se
simplifica en el caso de movimiento plano, es decir con ω perpendicular al plano de movimiento.
Ixx Ixy Ixz 0
Ixz
HG = Iyx Iyy Iyz 0 = Iyz ω
Izx Izy Izz G ω Izz G
HG = IG ω
donde IG es el momento de inercia respecto de un eje perpendicular al plano que pasa por el
centro de gravedad. Las condiciones anteriores aseguran que el momento es constante con el
tiempo y, por tanto, la derivada temporal de HG se puede escribir como
ḢG = IG ω̇
F = mAG
(5.1)
TG = IG ω̇
La ecuación 5.1 también vale para un punto fijo O, quedando entonces como
F = mAG
(5.2)
TO = IO ω̇
Las ecuaciones de 5.1 pueden reescribirse según el principio de D’Alambert, que quedarı́an
como
F + FIG = 0
(5.3)
TG + TGI = 0
donde
FIG = −mAG
TGI = −IG ω̇
Lo usual es que las fuerzas de inercia se reduzcan al punto G, por esto si este es el punto
O
elegido se podrá omitir la letra G, de este modo para una determinada barra a , en vez de TGaI y
FGa a a O
I , en adelante se escribirá T I y F I . Para una determinada barra a las ecuaciones anteriores
quedarán como
Fj + FIj = 0
I
TGj + TGj =0
Las fuerzas y pares con el superı́ndice I se refieren a los pares y fuerzas de inercia. Por otro
lado, respecto a las fuerzas y pares que aparecen sobre la barra, se distinguirá dos tipos. El
primero es el de las fuerzas externas aplicadas, marcadas con el superı́ndice A, y el segundo
O
aparece por la existencia de pares cinemáticos que unen a la barra j con otras barras, que
llevará el superı́ndice R. Por tanto las ecuaciones anteriores se pueden reescribir como
b
X
FA
j + FR I
jk + Fj = 0
k=1
Xb
A R I
TGj + TGjk + TGj =0
k=1
donde FR O
kj es la fuerza que aparece en el sólido k por la restricción que mantenı́a con el sólido
Oj.
Se puede multiplicar la ecuación de balance de fuerzas por un determinado conjunto de
desplazamientos y giros δr y δθ, a los que solo le exigirá respetar las restricciones del mecanismo.
Esta condición implicarı́a que, por ejemplo, las barras no cambiasen sus longitudes, ya que están
modelados como solidos inextensibles. Como otro ejemplo, el desplazamiento relativo de dos
barras ligados con una restricción prismática, deberá ser en la dirección de movimiento de ese
O
par. Ası́, para una barra j se podrı́a escribir
b
X b
X
FA
j · δrA,j + FR I A
jk · δrj + Fj · δrj + TGj δθj +
R
TGjk I
δθj + TGj δθj (5.5)
k=1 k=1
Las fuerzas de ligadura entre las barras del mecanismo responden a la ley de acción y reacción
de Newton y por tanto FR R R
kj = −Fjk . El trabajo Fjk · δrj será el mismo pero de signo contrario
que el FRkj · δrk . Esto permite la cancelación de todos los trabajos de las fuerzas de ligadura en la
expresión que se obtendrı́a de la suma, para todas las barras del mecanismo, de la expresión 5.5.
La suma de todos las expresiones 5.5 pero particularizadas para todas las barras del mecanismo,
quedarı́a como:
X b Xb X b
A I
TjA + TGj
I
Fj · δrA,j + Fj · δrj + · δθj = 0
j=1 j=1 j=1
Los desplazamientos y giros utilizados anteriormente no tienen por que ser pequeños, y en
este caso los trabajos vendrán con la integración temporal de las fuerzas y pares anteriores. Pero
CAPÍTULO 5. DINÁMICA 97
se puede imponer δrj → 0 y δθj → 0. En este caso se divide por el tiempo transcurrido en ese
desplazamiento, que también será δt → 0. El resultado es que la ecuación de trabajos anterior
pasa a la siguiente ecuación de potencias
b
X b
X b
X
0
FIj · vj0 +
0
FA
j · vA,j + TjA + TGj
I
· ωj = 0
j=1 j=1 j=1
donde vj0 y ωj0 representan a un conjunto de velocidades lineales y angulares, que no tienen por
que ser las reales, y a las que solo se les exige compatibilidad con las restricciones del mecanismo.
La prima se coloca para diferenciarlas de las velocidades reales. Por otro lado vA,j0 es la velocidad
del punto de aplicación de la fuerza en la barra j. Evidentemente, las velocidades del mecanismo
en un problema real podrı́an ser utilizadas en la expresión del principio de las potencias virtuales.
El principio de las potencia virtuales enuncia:
La suma de las potencias generadas por todas las fuerzas y pares, tanto aplica-
dos como de inercia, para un sistema de velocidades virtuales que cumpla con las
restricciones del problema, es nula.
nf np
X X
0= Fj · Vj + Ti ωi
j=1 i=1
! !
VCX −0.0067
0 = (0, −Mp g) + Tm ωe = (0, −1000) + 0.2 Tm
VCY 0.0133
donde el signo menos del producto Tr ωs indica que los sentidos de Tr y ωs están definidos como
opuestos.
La amplificación mecánica se define como el cociente entre el par entrante y el par saliente
del mecanismo. En este caso será, por tanto:
Tr ωe
Ap = =
Tm ωs
Por tanto, cuanto menor sea la velocidad de salida frente a la de entrada mayor será la
amplificación mecánica. Si la amplificación mecánica es mayor que uno en una determinada
posición, el mecanismo actúa como una reductora. Si Ap < 1 entonces se comportarı́a como una
multiplicadora.
El comportamiento del mecanismo de la figura 5.2 en cuanto a su aspecto amplificador de
par se puede evaluar gráficamente gracias al CIR Iaf , que por definición es el punto de las
O O
barras a y f que en el instante de estudio tienen la misma velocidad.
O O
Consideremos dos puntos C ∈ a y D ∈ f , entonces
! !
0 0
VC = = VD =
OC PD
ωe PD
Ap = =
ωs OC
De la expresión anterior se puede inferir el crecimiento de la amplificación dinámica cuan-
do, OC → 0, es decir, cuando Iaf se aproxima al punto O. Como podemos ver, esta circunstacia
O O
se produce cuando las barras a y e tienden a alinearse, lo que lleva al mecanismo a un punto
singular.
Figura 5.3: Mecanismo de cuatro barras con todas sus barras con masa y los vectores para la resolución
de su dinámica.
inciso porque la superposición no puede aplicarse a la posición, que no es problema lineal, y por
tanto no será válido usar este método para dos posiciones del mecanismo diferentes.
Una descomposición viable consiste en obtener problemas parciales que se obtienen como
resultado de descargar las barras de manera que se puedan calcular las reacciones fácilmente
para cada uno de los problemas. En la figura 5.3 se muestra un sencillo mecanismo cuya dinámica
se resolverá mediante el método de superposición. Concretamente se requiere conocer el par
motor necesario para que el mecanismo la barra de entrada se mueva con velocidad constante
ω2 = 1 rad/s. Las barras del mecanismo tienen una densidad lineal de masa de λ = 10 kg/m.
La resolución de este problema de dinámica inversa se requiere el cálculo de la cinemática.
El problema de posición es obvio. Las dimensiones y ángulos de los vectores son
√
ρ2 = 0.2 m ρ3 = 0.2 2 m ρ4 = 0.4 m
√
ρG2 = 0.1 m ρG3 = 0.1 2 m ρG4 = 0.2 m
θ2 = θG2 = π/2 rad θ3 = θG3 = π/4 rad θ4 = θG4 = π/2 rad
α3 = −0.741 rad/s2
! ! ! ! ! (
0 −0.3 0 −0.3 −0.2
α3 + α4 = −12 + 0.6672 + 12 ⇒
−0.3 −0.2 0.3 0 0.3 α4 = 1.111 rad/s2
Con todas las fuerzas y pares aplicados, ahora es cuando se comienza a utilizar el método de
superposición. Usándolo, el problema original se divide en otros dos a los que se llamará A y B.
En el A se ha descargado el acoplador, por lo que al ser éste una barra biarticulada y sin cargas,
las reacciones son colineales con la linea que une los pares de rotación, es decir son paralelas
O
a r3 . En el B es la barra 4 la que se encuentra descargada y es en ella donde las reacciones
son paralelas a r4 . Esta descomposición se muestra en la figura 5.4. Los pares motores que
equilibran los mecanismos A y B serán Tm A y T B , siendo el par motor requerido del problema
m
A
original Tm = Tm + Tm .B
CAPÍTULO 5. DINÁMICA 102
Figura 5.4: Descomposición del mecanismo original según el método de superposición con las flechas en
gris indicando dónde existen fuerzas y pares de inercia.
En este primer subproblema el acoplador se encuentra descargado por lo que las reacciones
tienen la dirección de la lı́nea que une las uniones de la barra. Teniendo esto en cuenta, se puede
calcular la reacción R mostrada en la figura 5.5. Para ello se planteará que el sumatorio de
O
momentos en la barra 4 respecto del punto O2 debe ser cero. Este toma la forma genérica de
X
MO2 |4 = 0 ⇒ r4 × R + rG4 × FI4 + TI4 = 0 (5.6)
pero si se tiene en cuenta que las reacciones articulares R y R̂ son paralelas a la barra 3 , es O
decir horizontales, entonces se puede escribir
i j k i j k 0
−0.2 0.3 0 + −0.1 0.15 0 + 0 = 0
este signo negativo de R se deriva que el acoplador no se encuentra trabajando a tracción como
se deducirı́a de las fuerzas dibujadas en 5.5 si hubieran resultado positivas, sino que trabaja a
compresión. Con este valor de R se puede plantear el sumatorio de momentos respecto a O1 en
O
la barra 2 para obtener TA m.
X
MO1 |2 = 0 ⇒ TA I
m + r4 × R̂ + rG2 × F2 = 0
Nótese que en la primera columna, tercera fila del primer determinante de la expresión anterior
aparece el valor 0.517, y no -0.517, ya que en este caso se utiliza R̂.
O
Planteando de forma genérica el balance de momentos en la barra 3 respecto al punto P,
se tiene
X
MP |3 = 0 ⇒ (−r3 ) × RQ + (−rG3 ) × FI3 + TI3 = 0
En principio la reacción R tiene dos componentes desconocidas, y del sistema anterior solo es
significativa la tercera ecuación por lo que se tienen dos ecuaciones con una sola incógnita. La
O
ecuación que falta se obtiene de considerar que en este caso es la barra 4 la que se encuentra
O
descargada. Esto implica que la reacción RQ es paralela a la barra 4 y por tanto es conocida
su dirección, tal y como se muestra en la figura 5.6. Esto proporciona una nueva ecuación con
la que resolver el problema, ya que
Q
RX 2 Q
= − ⇒ RX = −0.0667 RYQ
RYQ 3
CAPÍTULO 5. DINÁMICA 104
quedando la reacción RQ
−0.031
RQ = 0.046 N
0
Conocida RQ se puede determinar RP mediante los sumatorios de fuerzas horizontales y verti-
cales. Del sumatorio de fuerzas horizontales se tiene
P I Q
RX + F3,X + RX =0
P P
RX + 1.4 − 0.031 = 0 ⇒ RX = −1.369N
y del de verticales
RYP + F3,Y
I
+ RYQ = 0
RYP + 0.1 + 0.046 · 0.3 = 0 ⇒ RYP = −0.146N
B se obtendrá con el sumatorio de momentos respecto O en la barra 2
El par Tm 2 O
X
MO1 |2 = 0 ⇒ TB P
m + r2 × R̂ = 0
Por lo que el par motor Tm necesario para que el mecanismo se mueva con ω2 = 1 rad/s constante
será Tm = Tm A + T B = 0.12 − 0.205 = 0.085 N m.
m
Figura 5.7: Fuerzas, pares y reacciones en las barras utilizados para la resolución de la dinámica de un
cuatro barras por el método matricial.
tres barras móviles se pueden obtener tres ecuaciones correspondientes al sumatorio de fuerzas,
verticales, horizontales y sumatorio de momentos, resulta que se puede plantear un sistema de
nueve ecuaciones con nueve incógnitas.
O
Para la barra 2 del sumatorio de momentos respecto a O1 se obtiene
X
MO1 |2 = 0 ⇒ TB P I
m + r2 × R̂ + rG2 × F2 = 0
O
Para la barra 3 del sumatorio de momentos respecto a P se obtiene
X
MP |3 = 0 ⇒ TI3 + (−r3 ) × RQ + (−rG3 ) × FI3 = 0
O
siendo las ecuaciones de 4
O2 Q
RX + 0.240 + R̂X =0
RYO2 + 0.941 + R̂YQ = 0
−0.105 − 0.2 R̂YQ − 0.3 R̂X
Q
=0
3 r3
r4
r2
Tm 2
�� Figura 5.8: Análisis dinámico di-
O recto de un mecanismo de biela-
manivela.
tantos sistemas como grados de libertad tenga el mecanismo y también se podrán plantear el
mismo número de ecuaciones de potencias virtuales linealmente independientes.
Aunque el método de las potencias virtuales es rápido y directo, tiene el problema de no
permitir calcular las reacciones en las uniones de las barras. Para obtenerlas será necesario
romper las uniones y plantear las ecuaciones de igualdad de momentos y fuerzas en la barra en
las barras.
O
De esta manera se obtendrı́a la velocidad angular de la barra de entrada 2 para las diferentes
posiciones del mecanismo. Con esta evolución de ω2 para diferentes ω2 se podrı́a obtener también
O
la aceleración angular de 2 . Lo que no se sabrı́a es el tiempo invertido en llegar a cada una de
estas posiciones. Para ello se tendrı́a en cuenta que ω2 = dθ 2
dt
Z θ2
dt 1 1
dt = dθ2 = dθ2 ⇒ t = dθ2
dθ2 ω2 θ=0 ω2
Por ejemplo, para el caso del biela manivela, la integral serı́a
s
Z θ2
m4 ρ22 sin2 (θ3 − θ2 )
t= dθ2
θ=0 2Tm ∆θ2 sin2 θ3
Este teorema se concretiza en la definición de unos complejos i 6= f (ωe ) que toman la forma
1
Eci mi v 2 + 1 IGi ω 2
i = Pnb = Pnb2 1 Gi 2 2 1 i 2
j=1 Ecj j=1 2 mj vGj + 2 IGj ωj
lo que significa que los factores εi se pueden determinar para cada posición del mecanismo con
una velocidad de entrada del mecanismo arbitraria.
Puede resultar contradictorio que en la definición matemática de i aparezcan las velocidades
angulares, pero que después se diga que el valor de i no depende de ellas. Esto es ası́ ya que al ser
el problema de velocidades lineal todas las velocidades podrı́an escribirse como proporcionales a
una cierta velocidad ωe . De esta manera wj = kj ωe y VGj = kj0 ωe , siendo kj y kj0 constantes que
solo dependerı́an de la posición. Ası́ sustituyendo en la ecuación anterior, estas ωe2 se canceları́an
en el denominador y numerador, quedando una expresión de i independiente de las velocidades
angulares.
CAPÍTULO 5. DINÁMICA 109
El trabajo realizado por las fuerzas exteriores (∆W ) desde la posición inicial hasta otra
cualquiera se invierte en incrementar la energı́a cinética del mecanismo:
Ecf in = Ecini + ∆W
y puesto que ahora se conoce la fracción de ese trabajo que incrementa la energı́a cinética de la
barra en cuestión, a partir de ella se podrı́a obtener fácilmente la velocidad de esa barra.
Por ejemplo, para el caso del mecanismo de biela manivela anterior, pero en el que ahora la
O
barra 2 ya tiene una I2O no despreciable. Por el momento es desconocida la velocidad angular.
A este valor desconocido se le llamará ω20 . Puesto que el punto O es un punto fijo entonces se
puede escribir la energı́a cinética de esa barra. También se podrı́a escribir la energı́a cinética de
O 0
4 poniendo V4 en función de ω como
2
0 1 1 sin2 (θ3 − θ2 )
Ec2 = I2O ω202 0
Ec4 = m4 ρ22 ω202
2 2 sin2 θ3
quedando 2 como
1 02
2 I2O ω2 I2O
2 = 2 = 2
1 02 1 2 02 sin (θ3 −θ2 ) (θ3 −θ2 )
I2O + m4 ρ22 sinsin
2 I2O ω2 + 2 m4 ρ2 ω2 sin2 θ3 2θ
3
3 r3
r4
r2
2 Figura 5.9: Inercia y Par reducidos
Tm
�� a un eje en un mecanismo de biela-
O manivela en el actúa un par Tm y una
fuerza F.
igualando las expresiones anteriores se obtiene el valor de la inercia reducida a la barra O, que
2
es
sin2 (θ3 − θ2 )
IR = I2O + m4 ρ22
sin2 θ3
donde se ve la independencia con ω2 y la dependencia con θ2
Por otro lado, los mecanismo, en general, están sometidos a fuerzas y pares aplicados en
distintas barras del mecanismo. La potencia añadida o substraı́da del mecanismo por todas
estas acciones será:
Xn
P = Fi · vi + Ti ωi
i=1
de este modo también se puede definir un par asociado a un eje de reducción (TR ) de manera
que la potencia desarrollada por ese par reducido sea la misma que la del conjunto de fuerzas y
pares actuantes. De este modo
n Pn
X
i=1 Fi· vi + Ti ωi
Fi · vi + Ti ωi = TR ωR ⇒ TR =
ωR
i=1
Por ejemplo, si en el mecanismo de biela manivela se añade una fuerza hacia abajo de valor
F aplicada en el émbolo, tal y como se muestra en la figura 5.9, entonces el par reducido al eje
O
que pasa por O y gira con 2 serı́a
sin(θ3 − θ2 )
P = Tm ω2 − F ρ̇4 = Tm ω2 − F ρ2 ω2 = TR ω2
sin θ3
donde el signo negativo en el término de potencia de la F viene porque la fuerza va hacia abajo
y ρ̇4 está definido positivo hacia arriba. Se puede despejar TR y, cancelando las ω2 , queda
sin(θ3 − θ2 )
TR = Tm − F ρ2
sin θ3
Lo habitual es que el eje de reducción de la inercia y el par sea el mismo. En este caso la
inercia y el par reducidos a un eje permiten analizar la dinámica del mecanismo como si fuese un
único eje de inercia variable IR (θ) sometido al par TR (θR ). Ası́, el trabajo desarrollado por TR
CAPÍTULO 5. DINÁMICA 111
entre las posiciones θRini y θ f in está relacionado con la velocidad del eje de reducción mediante
R
la expresión
Z θf in
R 1 f in 2 f in 1 ini 2 ini
TR dθR = IR (θR )ωr (θR ) − IR (θR )ωr (θR )
θR =θRini 2 2
Para calcular la aceleración del eje de reducción se utilizará la ecuación del momento cinético
d
TR = (IR ωR )
dt
pero en general IR no es constante, y por tanto
dIR dIR 2
TR = ωR + IR αR = ω + IR αR
dt dθR R
de donde se puede despejar αR .
Capı́tulo 6
6.1. Problema
O
Calcule el par motor Tm que es necesario aplicar a la barra 2 para mover el mecanismo
de la figura 6.1. Se usará el método de las potencias virtuales. Las fuerzas de inercia no se
despreciarán en esta ocasión. La velocidad de entrada del mecanismo es constante y de valor
O
ω2 = 1 rad/s. El carro 4 tiene su centro de gravedad en G, su masa es M = 4 kg y su inercia
es IG = 0.01 kgm2 . La gravedad actúa en dirección Y negativa.
112
CAPÍTULO 6. PROBLEMAS RESUELTOS DE DINÁMICA 113
Dentro del análisis cinemático de velocidades, será necesario conocer la velocidad de los
puntos con masa, la velocidad de los puntos de aplicación de fuerzas, la velocidad angular de las
barras con momento de inercia no despreciable y la velocidad angular de las barras donde estén
aplicados pares externos. En este caso será necesario conocer VG y ω4 (ω2 es el dato inicial).
Para ello se usará el bucle r1 + r2 = r3 . Derivándolo se obtiene
y análogamente
ρ3 eiθ3 = 0.2 + i 0.2 ⇒ ρ3 ieiθ3 = −0.2 + i 0.2
En el caso del sumando ρ̇3 eiθ3 se prescindirá de ρ̇3 . Las componentes del vector unitario eiθ3
se obtendrán dividiendo ρ3 eiθ3 por ρ3 .
0.2 + i 0.2 1 1
eiθ3 = √ = √ + i√
0.2 2 2 2
Ahora se pueden escribir las partes reales e imaginarias de la ecuación ω2 ρ2 ieiθ2 = ρ̇3 eiθ3 +
ω3 ρ3 ieiθ3 , y ası́ obtener el sistema a resolver
! √ ! !
0.1 1/ 2 −0.2
1 = ρ̇3 √ + ω3
0.2 1/ 2 0.2
r̈1 + r̈2 = r̈3 ⇒ −ω22 ρ2 eiθ2 = ρ̈3 eiθ3 − ω32 ρ3 eiθ3 + α3 ρ3 ieiθ3 + 2ρ̇3 ω3 ieiθ3
ecuación cuyas componentes reales e imaginarias se pueden escribir de forma análoga a como se
ha procedido con las de velocidad.
! √ ! ! ! √ !
2 −0.2 1/ 2 1 0.2 −0.2 0.3 1 −1/ 2
1 = ρ̈3 √ − 2 + α3 + 2√ √
0.1 1/ 2 4 0.2 0.2 2 4 1/ 2
O
Las fuerza y par de inercia en 4 serán:
!
I 0.55
F = −M AG = N ; T I = −IG α3 = −0.0037 N m
−0.75
En este momento ya están calculados todos los elementos requeridos. Ahora se puede plantear
la ecuación de las potencias virtuales. En ella el par motor T m es la única incógnita.
FI · VG + M g · VG + T m ω2 + T I ω3 = 0
! !
0.15 0.15
0.55 −0.75 + 0 −39.2 + T m · 1 − 0.0037 · 0.25 = 0
−0.25 −0.25
CAPÍTULO 6. PROBLEMAS RESUELTOS DE DINÁMICA 115
6.2. Problema
Calcule la fuerza F necesaria que debe suministrar el actuador de la figura 6.2 para que la
O O
barra 3 gire con una velocidad angular constante de ω3 = 1 rad/s. La barra 3 tiene una masa
M = 100 kg y una inercia respecto del centro de gravedad G de IG = 100 kgm2 . La masa e
inercia de las demás barras se considerará despreciable. La gravedad actúa en dirección vertical
negativa.
Restando la segunda ecuación (parte imaginaria) a la primera (parte real) multiplicada por
√
2 se despeja ω4 = −1 rad/s. Y sustituyendo ω4 en la primera ecuación se obtiene ρ̇2 = 5 m/s.
También es necesaria VG = ṙ4 + ṙ3 .
! ! !
iθ4 iθ3 −1 −3 −2
VG = ω4 ρ4 ie + ω3 ρ3 ie = −1 +1 = m/s
−3 1 4
Las velocidades y aceleraciones necesarias para plantear la ecuación del Teorema de las
Potencias Virtuales también se pueden obtener gráficamente. Los cinemas de velocidad y ace-
O
leración se muestran en la figura 6.3. Para la resolución gráfica se usarán los puntos A ∈ 2 y
O
B ∈ 3 . Las ecuaciones que se usarán para calcular las velocidades serán:
VA = ω3 × rOA
VB = VA + vAB
VB = ω4 × rPB
VG = VB + ω3 × rBG
AA = ω3 × (ω3 × rOA )
AB = AA + aAB + 2ω3 × vAB
AB = ω4 × (ω4 × rPB ) + α4 × rPB
AG = AB + ω3 × (ω3 × rBG )
Una vez obtenidas las velocidades y aceleraciones, ya sea de forma gráfica o mediante com-
plejos, ya se puede plantear la ecuación de las potencias virtuales para el cálculo de T m . Se
usará como campo de velocidades virtuales las velocidades reales calculadas anteriormente. La
diferencia respecto al primer ejemplo es que en esta caso la fuerza la proporciona el actuador. En
CAPÍTULO 6. PROBLEMAS RESUELTOS DE DINÁMICA 117
O O
este sistema aparecen dos fuerzas, una sobre 2 y otra igual y contraria en 3 . Estas fuerzas son
colineales y actúan en la linea de puntos de la figura 6.3. Supongamos que esas fuerzas actúan
O O
en los puntos A ∈ 2 y B ∈ 3 . En la ecuación de las potencias virtuales aparecerá por tanto
O O
donde vAB es la velocidad con que el punto A ∈ 2 ve moverse al B ∈ 3 . Ahora se debe tomar
en consideración la colinealidad de la fuerza F con vAB . Esto implica que el producto escalar
F · vAB se convierta en el producto de los módulos de F y vAB . Pero la velocidad con la que se
despliega el actuador es la variación por unidad de tiempo del módulo de r2 , que por definición
es ρ̇2 , y por tanto
... + F · vAB + ... = ... + F ρ̇2 + ...
0 = −M AG · VG + M g · VG + F ρ̇2
que sustituyendo da
! !
−2 −2 √
0 = −100 5 5 + 0 −100 · 9.8 +F 5
4 4
1000 + 3920
F = √ = 2200 N
5
CAPÍTULO 6. PROBLEMAS RESUELTOS DE DINÁMICA 119
6.3. Problema
1. Calcule la velocidad y aceleración que hay que suministrar en el actuador de la figura 6.4
O
para que la barra 3 gire con una velocidad angular de ω3 = 1 rad/s y una aceleración
angular de α3 = 1 rad/s2 .
2. Calcule la fuerza F necesaria que debe suministrar el actuador para obtener las condiciones
O
cinemáticas del punto anterior. La barra 3 tiene una masa M = 100 kg y una inercia
respecto del centro de gravedad G de IG = 100 kg m2 . La masa e inercia de las demás
barras se considerará despreciable.
El bucle de vectores que se utilizará será r1 + r4 + r2 = r3 , que derivado queda ṙ4 + ṙ2 = ṙ3 .
Esta ecuación pasa a:
Estos complejos se pueden reescribir fácilmente observando la cuadrı́cula de la figura 6.5. Te-
niendo en cuenta que ω3 = ω2 = 1
! ! ! ! (
−1 0 −2 0 ω4 = −2 rad/s
ω4 + ρ̇2 +1 =1 ⇒
−3 1 0 3 ρ̇2 = −3 m/s
Para las aceleraciones se derivará de nuevo el bucle anterior, obteniéndose r̈4 + r̈2 = r̈3 . Del que
se obtiene:
−ω42 ρ4 eiθ4 + α4 ρ4 ieiθ4 + ρ̈2 eiθ2 − ω22 ρ2 eiθ2 + α2 ρ2 ieiθ2 + 2ρ̇2 ω2 ieiθ2 = −ω32 ρ3 eiθ3 + α3 ρ3 ieiθ3
6.4. Problema
El mecanismo de la figura denomina yugo inglés y una de sus aplicaciones es la de herramienta
O
para machacar. En el mecanismo de la figura la barra 3 tiene una masa de M3 = 60 kg y una
O
IG3 = 3 kg m2 . Respecto a la barra 4 , su centro de masas está en G4 siendo M4 = 4 kg y
IG4 = 0.02 kg m2 . En el punto A se encuentra un par de rotación que une a las barras 2 y 4 . O O
O
La barra impulsora es la 2 que tiene una masa de M2 = 40 kg y un centro de inercia respecto
al punto O de IO = 10 kg m2 .
2. Calcule el par T2 para que el sistema venza una fuerza en B de F = 10000 N cuando la
barra de entrada gira a ω2 = 10 rad/s constante.
Al punto B se llega mediante la suma de los vectores ṙB = ṙ1 + ṙ5 + ṙ6 , donde ṙ5 y ṙ6 son
vectores cuyo módulo y argumento no cambia, por lo que
!
−1
ṙ@
VB = ṙ1 + @ ṙ@
5 +@ 6 =2 m/s
0
Para el cálculo de las aceleraciones, se derivará dos veces el bucle inicial, obteniéndose
r̈1 + r̈3 = r̈2 ⇒ ρ̈1 eiθ1 + ρ̈3 eiθ3 = −ω22 ρ2 eiθ2 + α2 ρ2 ieiθ2
y sustituyendo
!
−2
T1 · 10 + (−10 − 4(0.2 + 0.4 ))0.5 · 10 + (−60) 39.9 0 = 0 ⇒ T1 = −472.6N m
| {z }
2
0
ρ2
siendo la aceleración AG3 = r̈1 = (40, 0)T . Por lo que el nuevo PTV será
y sustituyendo
! !
−2 −2
T2 · 10 + (−60) 40, 0 + 10000, 0 = 0 ⇒ T2 = 1520 N m
0 0
CAPÍTULO 6. PROBLEMAS RESUELTOS DE DINÁMICA 124
6.5. Problema
Un mecanismo de cuatro barras se encuentra girando con θ2 = 1 rad s en el instante inicial
mostrado en la figura 6.7. A partir de este momento y durante una vuelta se aplica un par
Tm = 10 N m. La barra tiene una inercia respecto al punto fijo O de I2,O = 5 kgm2 y el acoplador
tiene una masa de m3 = 20 kg y una inercia respecto de su centro de gravedad de I3,G = 12 kgm2 .
O
¿Qué velocidad angular tendrá 2 tras esta actuación de Tm ?
E1:100
0 1 2m r3
G G
3 r3/2
4
r4
2 r2
Tm
r1
O
Por tanto la energı́a cinética final de la barra 2 será
f in 1
Ec,2 = 2 · Ecf in = 9.98 J = 5 · (ω2f in )2
2
O
La velocidad angular de 2 tras la actuación del par Tm será entonces
r
f in 9.98 · 2 rad
ω2 = = 1.998
5 s
CAPÍTULO 6. PROBLEMAS RESUELTOS DE DINÁMICA 126
6.6. Problema
1. Calcule la velocidad del punto G del mecanismo de la figura. La velocidad angular de
entrada es ω2 = 1 rad/s. Los contactos son sin rozamiento.
E1:10
0 0.1 0.2 m
3 C
B G
2
A
Tm
O
E1:10
0 0.1 0.2 m rD
D r6
3 C r4
r3 r2
B G
2 r5
A
Tm r1
O
donde desaparece ṙ5 ya que este vector no cambia ni de dirección ni de módulo. Si se tiene en
cuenta que θ4 = θ2 + π2 , por lo que ω2 = ω4 , entonces este bucle se desarrolla como
La velocidad VG será
! ! !
1 1 −0.4 1 0.1 −0.3 m
VG = ṙ1 − ṙ3 = −ρ̇1 − ω3 ρ3 ieiθ3 = − (−2) ⇒ VG =
2 2 0 2 −0.3 −0.3 s
se les colocará una prima para indicar que nos son las energı́as cinéticas que tienen las barras
en este segundo apartado.
0 1 2
Ec2 2 IO,2 ω2
2 = =
Ec0 1 2 1 2 1 2
2 IO,2 ω2 + 2 IG,3 ω3 + 2 m3 VG
1
2 0.015 · 12
2 = 1 2 1 2 1 2
⇒ 2 = 0.0194
2 0.015 · 1 + 2 0.01 · (−2) + 2 4 · (2 · 0.3 )
O
A la velocidad angular de 2 en este segundo apartado se le llamará ω̂ 2 . La energı́a cinética
O
real de la barra 2 para este segundo apartado será
1
Ec2 = IO,2 ω̂22 = 2 Ec = 2 · ∆θ2 · Tm
2
expresión de la que se puede despejar ω̂2
s r
22 · ∆θ2 · Tm 2 · 2.067 × 10−3 · 0.1 · 2 rad
ω̂2 = = ⇒ ω̂2 = 0.7184
IO,2 0.015 s
! !
ω̂2 V̂G −0.3 −0.2155 m
= 0.7184 = ⇒ V̂G = 0.7184 · VG = 0.7184 ⇒ V̂G =
ω2 VG −0.3 −0.2155 s
Capı́tulo 7
Sı́ntesis de Mecanismos
129
CAPÍTULO 7. SÍNTESIS DE MECANISMOS 130
barra de salida se conceta al émbolo, que en su movimiento le dará soporte a la inspiración del
paciente cuando suba, y le ayudará en en la expiración mientras baja.
La función de volumen de aire en función del tiempo se le suele llamar aleta de tiburón. Esta
función se encuentra en la literatura médica, y tiene una enorme bastante variabilidad según
las personas, aunque está usualmente establecido que la inspiración toma aproximadamente el
doble de tiempo que la expiración.
La función se muestra en la figura 7.1, y se ha ajustado a ciertas figuras geométricas simples
para un tratamiento matemático más sencillo. Tampoco es crucial un ajuste escrupuloso a una
determinada gráfica, ya que la variabilidad de estas funciones entre las distintas personas es
bastante elevada, por esto es más importante que se adapte a la forma general, manteniendo,
eso sı́, la distribución de tiempos entre la ida y la vuelta del mecanismo, es decir entre la
inspiración y la expiración. También es importante que los caudales máximos de inspiración y
expiración se den en las partes adecuadas del ciclo. Estos flujos son la variación con el tiempo
del volumen, es decir, la derivada de la función 7.1.
Figura 7.2: Esquema de respirador: a) con un mecanismo de retorno rápido. b) con un mecanismo de
cuatro barras.
espacio imponiendo precisamente que los tiempos de ida y vuelta estén en proporción de 2:1.
Se considerá que la barra de salida generará un arco cuyos puntos iniciales y finales B inf y
B sup estarán en un segmento vertical. Esto añade como variable el ángulo θ1 que forman
las dos articulaciones con la barra fija del mecanismo (O1 y O2 ).
La distancia entre los puntos B inf y B sup será igual a 1. Posteriormente a la sı́ntesis, el me-
canismo se escalará para llegar a la carrera deseada en función del volumen de respiración
y la superficie del émbolo del fuelle.
El ángulo que forman los segmentos O1 Asup y O1 Ainf será igual 2π/3, esto establecerá la
proporción 2:1 entre expiración e inspiración.
en los puntos extremos superior de la carrera los segmentos O1 Asup y Asup B sup serán
colineales y superpuestos parcialmente
en los puntos extremo inferior de la carrera los segmentos O1 Ainf y Ainf B inf serán coli-
neales, pero solo tendrán en común el punto Ainf .
CAPÍTULO 7. SÍNTESIS DE MECANISMOS 132
Figura 7.3: Esquema para el cálculo de las variables del mecanismo de cuatro barras.
O
Para obtener la longitud de la barra fija 1 , ρ1 = |O1 O2 | se necesitarán previamente el ángulo
α, que se determina teniendo en cuenta que el triángulo O2 B 0 B es isósceles, y también γ que se
obtendrá con el teorema del seno sobre el triángulo O1 B sup B inf .
1
α = arc cos
2ρ4
√ (7.2)
1 ρ3 − ρ2 3(ρ3 − ρ2 )
= ⇒ γ = α − arcsin
sin(π/3) sin(α − γ) 2
Figura 7.4: Vectores para el cálculo del mecanismo para una posición θ2 de la barra de entrada.
π
δ= −α−β
2
8π 19π (7.6)
θ2ini = + δ ⇒ θ2ini = −α−β
3 6
Una vez que el mecanismo está totalmente determinado se calculará la posición del émbolo
O
6 del fuelle a partir del ángulo θ2 . Para ello se utilizarán los complejos rj de la figura 7.5.
Para conocer la posición de B primero se calculará la longitud del segmento O2 A, que es
q
|O2 A| = ρ21 + ρ22 + 2ρ1 ρ1 cos(θ2 + φ) (7.7)
a partir de ella se obtienen los ángulos O\1 O2 A y AO2 B mediante el teorema del coseno, del que
\
se obtendrá el ángulo θ4 de la figura 7.5
A partir de la posición de B se puede calcular fácilmente la posición del émbolo (el punto B )
O
conocidas la longitud de la barra 5 y la situación del segmento DC, que es el eje del émbolo.
En principio se podrı́an incluir la distancia O2 B y ρ5 = |r5 | en el análisis de optimización, pero
es probable que afecten poco. Si se añaden, además habrá que añadir alguna guı́a que garantice
el movimiento vertical del punto C.
Si se dejan como parámetros ρ2 y ρ4 para el análisis de optimización, lo lógico es establecer
la longitud del segmento O2 B como la media entre ρ4 y ρ4 cos(α/2). Con esta construcción, si
consideramos el origen de los ejes en O2 , el punto medio de la carrera del émbolo será (D, −ρ5 ).
La posición vertical del punto YC (θ2 ) para un determinado ángulo θ2 se obtendrá a partir del
ángulo que falta, θ5 . A partir de él la distancia YC (θ2 ) = |DC| que representa la posición del
émbolo, será:
ρ4 cos θ4 − |O2 D|
θ5 = arc cos
ρ5 (7.9)
YC (θ2 ) = |O2 D| = ρ4 sin θ4 − ρ5 sin θ5
72
X
S(ρ2 , ρ4 ) = (YC (nπ/36) − F obj (nπ/36)2 (7.10)
n=1
Figura 7.5: Vectores para el cálculo del mecanismo para una posición θ2 de la barra de entrada.
Capı́tulo 8
8.1. Introducción
En los capı́tulos precedentes se calcularon las posiciones velocidades, aceleraciones, fuerzas
y pares que actúan sobre las barras de mecanismos con una determinada geometrı́a. Con la
sı́ntesis se determinan las dimensiones de las barras para que cumplan determinadas funciones.
Esto afecta a la posición de las articulaciones y de los puntos de aplicación de fuerzas y pares
tales como pares motores y resistentes, fuerzas de inercia, gravedad, etc. Pero para que la barra
cumpla su función también es necesario especificar espesores y materiales. Si la elección no es
correcta, cuando se apliquen las cargas, las deformaciones o las tensiones pueden ser excesivas
y dar lugar a un mal funcionamiento. Por tanto, un análisis de las cargas es fundamental en el
diseño de máquinas.
La selección adecuada de un elemento de máuina, o su diseño, consiste en calcular los es-
fuerzzos y deformaciones que se esperan del servicio de la máquina y, luego, elegir los tamaños,
materiales o espesores adecuados de manera que no se excedan las tensiones y deformaciones
crı́ticas. El primer paso para el cálculo de las tensiones y deformaciones es la determinación
exacta de las cargas.
136
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 137
4. Carga cı́clica: Una carga cı́clica varı́a, incluso en su signo, de manera periódica con una
frecuencia caracterı́stica.
1. Carga concentrada: La carga se aplica en una porción muy pequeña del área total del
elemento. Un ejemplo serı́a el contacto entre una rail y la rueda de un tren.
3. Carga volumétrica: La carga volumétrica o por unidad de volumen es la carga que tiene
una pieza en cada uno de sus puntos con masa, ya sean externos o internos. Un ejemplo
de estas cargas serı́an el peso o las fuerzas de inercia.
1. Carga normal (N): La carga pasa a través del centroide 1 de la sección. La carga puede ser
de tracción, que se considera positica, o de compresión, que se toma como negativa. Si la
carga no pasa por el centroide, tendrá un equivalente a una carga normal más otras cargas
flexionantes. Se considera que la carga normal es positiva si es de tracción y negativa para
la compresión.
2. Carga cortante (V): Esta carga es colineal al esfuerzo cortante y lo provoca. Si esta carga
no es colineal al cortante, sino simplemente paralela, será equivalente a una carga cortante
más una carga de torsión. El criterio de signos para los esfuerzos cortantes será el siguiente:
Una fuerza cortante es positiva si la dirección de la fuerza y la dirección normal de la sección
son ambas positivas o ambas negativas.
3. Carga flexionante (M): Es la que produce esfuerzo flector. Aunque su distinción respecto
a la anterior puede ser difusa, está relacionada con las zonas de aplicación y geometrı́a de
la pieza. La carga en el extremo libre de una viga esbelta empotrada produce sobre todo
esfuerzo flector hasta el punto que el esfuerzo cortante puede ser despreciado en el análisis
de tensiones, este serı́a un ejemplo de carga flexionante. Se considerará que un flector es
positivo si la zona de tensiones de tracción que origina ese flector está en la parta positiva
del eje. Ası́, para un flecotr MZ , si se han dibujados los ejes de la sección con el eje Y
hacia arriba, este será positivo si las tensiones de tracción están en la parta de arriba. En
la figura 8.1a) se muestran flectores positivos y en 8.1a) negativos.
4. Carga de torsión (T): Esta carga retuerce el elemento según su dirección principal. El
criterio de signos para los torsores será el de la regla de la mano derecha.
1
El centroide de de un área, o centro de gravedad de un área, se refiere al punto C(xC , yC ) que define su centro
R R
geométrico. Matemáticamente, el centroide cumple A (x − xC )dA = 0 y A (y − yC )dA = 0
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 138
5. Carga combinada: Es la carga que produce más de uno de los efectos de las cargas ante-
riores. Una carga normal pero que no pase por el centroide.
f G (rP ) = ρ(rP )g
AP = AG + α × rP + ω × (ω × rP )
f I (rP ) = −ρ AG − ρ α × rP − ρ ω × (ω × rP )
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 139
Figura 8.3: Fuerzas por unidad de longitud en una barra esbelta para el cálculo de tensiones internas
En la figura 8.2 se muestra una barra esbelta de longitud l, masa m y densidad uniforme que
gira con velocidad angular ω y aceleración angular α. El punto central coincide con el centro de
gravedad y su aceleración es AG . En este ejemplo, en vez de las fuerzas por unidad de volumen,
se calcularán las fuerzas por unidad de longitud. Para ello se usará la densidad lineal λ = m/l. El
diferencial de masa estará designado por el vector rP . También se supone que actúa la gravedad.
Para este caso las fuerzas por unidad de longitud f dl (rP ) serán:
f dl (rP ) = λ g − λ AG − λ α × rP − λ ω × (ω × rP ) (8.1)
En el lado derecho de la figura 8.2 se muestran la fuerzas de inercia que se han utilizado en el
análisis dinámico de los mecanismos. Pero si se requiere calcular las tensiones dentro de la barra,
estas fuerzas puntuales deben sustituirse por la distribución de fuerzas por unidad de longitud
que se muestran en la figura 8.3, donde en la parte inferior izquierda se muestra la fuerza por
unidad de longitud que corresponderı́a al primer termino del segundo miembro de la ecuación
(8.1). Arriba a la izquierda se muestra la distribución constante que corresponde al segundo
término. Abajo a la derecha se muestra la distribución que corresponde al tercer término. Esta
distribución es lineal y depende de la distancia del elemento de masa al punto G. Por último,
arriba a la derecha se muestra la distribución que representa el último termino de la ecuación
(8.1). El valor de ésta también depende linealmente de la distancia a G.
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 140
2. Si la viga está empotrada en un extemo y libre en el otro se dice que es una viga en
voladizo.
3. Una viga suspendida tiene uno o ambos extemos extendiéndose libremente más allá de sus
apoyos.
d
x x
W·d W
W·d
y
W
Mz
W T Vy
MZ x
W·d
Vz N
VY x x
z My
-W
Figura 8.5: Diagrama de cuerpo libre en una viga en voladizo carga da en su extremo libre con una
fuerza W.
Figura 8.6: Diagramas de cuerpo libre en una viga biapoyada cargada en su centro y una viga suspendida
cargada en sus extremos.
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 142
Figura 8.7: Diagrama de cuerpo libre en una eje biapoyado con dos discos cargados.
distancia entre discos es d3 = 0.14 m. Los radios de los discos son r1 = 0.1 m y r2 = 0.15 m. Con
estos datos se calcularán las reacciones, esfuerzos y el valor de la fuerza F2 .
P
Del sumatorio de momentos respecto al eje X MX = 0 se obtiene
X
MX = 0 ⇒ −F1 r1 + F2 r2 = 0 ⇒ F2 = 80 N
Una vez calculada F2 , se obtienen las reacciones en los dos extremos con los sumatorios de
momentos y fuerzas.
P
MY |O = 0 ⇒ F2 (d1 + d3 ) − R2Z (d1 + d3 + d2 ) = 0 ⇒ R2Z = 56.47 N
P
MZ |O = 0 ⇒ F1 (d1 ) − R2Y (d1 + d3 + d2 ) = 0 ⇒ R2Y = 35.29 N
P
FY = 0 ⇒ −R1Y + F1 − R2Y = 0 ⇒ R1Y = 84.76 N
P
FZ = 0 ⇒ R1Z − F2 + R2z = 0 ⇒ R1Z = 23.53 N
Con las reacciones es fácil obtener los esfuerzos. Estos se muestran en la figura 8.7.
Con todos los esfuerzos obtenidos, quedarı́a designar la sección más desfavorable en el pro-
blema. De la figura 8.7 se infiere que las esfuerzos tiene sus valores más altos en valor absoluto en
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 143
la parte central. De hecho entre los discos hay esfuerzo un torsor importante. Lso esfuerzos flec-
tores tiene sus valores más latos cerca de los discos. De mabos hechos se infiere que las secciones
problemáticas serán aquellas cercanas a los discos pero entre ellos, concretamente las secciones A
y B marcadas en la figura. Puesto que ambas secciones tiene el mismo flector, quedarı́a dilucidar
cual de ellas tiene el momento flector más grande. El valor de flector total será la composición
de las componentes Y y Z, que será
p
MA = 2.3532 + 8.4762 = 8.796 N m
p
MB = 3.5242 + 5.6472 = 6.656 N m
1. La normal y las componentes de las tensiones tangeciales que actúan en lados opuestos del
elemento diferencial deben ser iguales en magnitud pero opuestos en dirección.
2. Para que se cumpla el equilibrio de momentos, es necesario que las fuerzas tangenciales
sean simétricas, es decir τij = τji .
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 144
Figura 8.9: Pieza en situación de tensión plana y elemento diferencial de volumen en vista tridimensional
y en vista plana.
Con estos condicionantes los valores diferentes de las tensiones que actúan en el cubito se
reducen de 9 a 6. Para el tratamiento de esta nueva entidad, por tanto no podrá usarse un
vector, y se utilizará el tensor de dimensión 3 × 3 simétrico σ.
σxx τxy τxz
σ = τxy σyy τyz
0 0 σ3
Figura 8.10: Fragmento de eje cargado y elemento diferencial de volumen en vista tridimensional y en
vista plana.
σφ A = τxy A sin φ cos φ + τyx A cos φ sin φ + σx A cos φ cos φ + σy sin φ sin φ
σ +σ
no principales, siendo este ( x 2 y , 0). El radio será σ1 −σ
2 . Este radio también se puede obtener
2
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 146
Figura 8.12: Cı́rculo de Mohr para tensión plana. Criterios de signos para las tensiones tangenciales y
normales.
Además vemos que el argumento del seno y del coseno es 2θ. Por esto los giros en el cı́rculo
de Mohr para pasar de las coordenadas (σ, τ ) de una dirección a otra, serán el doble del giro
que se da en la orientación del material. El criterio de signos utilizado será el propuesto en este
libro [9] que aparece en la bibliografı́a [9]. Para la σ, serán consideradas positivas las de tracción
y negativas las de compresión. Para las tensiones tangenciales τ , las positivas en la cara x serán
hacia abajo, y las negativas en la cara y serán hacia arriba.
Para ilustrar el manejo del circulo de Mohr se ha incluido la figura 8.13. Supongase un
estado un punto del cuerpo donde existe el estado tensional de la figura 8.13a) donde se tiene
una σx = 30M P a, una σy = 10M P a y una τxy = 15M P a. El cı́rculo de Mohr tendrı́a su
√
centro en ( 30+10
2 , 0), siendo su radio 102 + 152 = 18 M P a. En la figura 8.13b) se mostrarı́an
las tensiomes normales y tangenciales para las dos orientaciones x e y en el circulo de Mohr del
cubo diferencial.
Se pueden encontrar las direcciones principales fácilmente a través del Cı́rculo de Mohr. El
proceso se explica con la figura 8.14. Supóngase que se tiene un punto con el estado tensional
de la figura 8.13. Las direcciones principales son aquellas sin τ . Por tanto se deberá pasar del
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 147
Figura 8.14: Cálculo de las direcciones principales mediante el circulo de Mohr. El giro necesario para
ir alas direcciones principales en el cı́rculo de Mohr será 2θ = tan−1 10
15
antihorario. Por tanto el giro en
el sólido será igual a θ = 2 tan−1 10 .
1 15
Figura 8.15: Cálculo de las orientación de mayor tensión tangencial. El giro en el cı́rculo de Mohr para
llegar a la situación de máxima τ será 2θ = tan−1 10
15 en sentido horario. Por tanto el giro en el sólido
1 −1 10
será igual a θ = 2 tan 15 .
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 148
punto inicial con σ = 15 y τ = 10 de la cara x y etiquetado con el punto gris oscuro, al punto
gris claro. El giro en el cı́rculo de Mohr será 2θ = tan− 1 10
15
. Por tanto el giro en el sólido será
1 15
igual a θ = 2 10 . Puesto que el radio del circulo de Mohr es 18 M P a, las tensiones principales
serán σ1 = 38 M pa y σ2 = 2 M pa.
También se puede encontrar la dirección de máxima tensión tangencial mediante el cı́rculo
de Mohr. El proceso se muestra en la figura 8.15.
Lı́mite elástico o tensión de plastificación (Punto B). Este punto está muy cercano al
anterior y marca la tensión máxima que el material puede soportar sin sufrir una deforma-
ción permanente. Si el material, en su funcionamiento, recibe una tensión σ 0 < σf , cuando
esta carga es retirada el material vuelve a la situación de partida tanto en tensiones como
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 149
Tensión de rotura o tensión última (Punto U). En este punto se llega al mayor nivel
de tensiones que puede sufrir el material, esta es σu .
Figura 8.17: Diagramas σ- para un acero. a) Puntos significativos. b) Proceso de carga y descarga
elástico. c) Proceso de carga y descarga a una tensión superior al lı́mite elástico. d) Proceso de carga y
descarga en un material que ya ha plastificado.
Figura 8.18: Circulos de Mohr del material durante el ensayo de tracción para σ < σ 0 < σ 00 .
Material E [MPa]
Bronce 120 000
Cobre 110 000
Hierro colado < 175 000
Hierro forjado 190 000
Acero 210 000
Magnesio 45 000
Titanio 107 000
Nı́quel 22 000
Aleaciones de Mg 42 000
Aleaciones de Al 70 000
Latón 110 000
Zafiro 420 000
Diamante sintetizado 491 000
Grafeno 1 000 000
Polietileno, Nylon 1 400
εtrans.
ν=−
εlong.
Para un material isótropo elástico perfectamente incompresible, este es igual a ν = 0.5. La
mayor parte de los materiales prácticos en la ingenierı́a rondan entre ν ∈ [0, 0.5]. En la tabla
8.3 se relata ν para una serie de materiales.
Material G [MPa]
Acero 85 000
Aluminio 26 300
Latón 39 200
Fundición gris (4.5 % C) 41 000
Bronce 41 000
Cobre 42 500
Hierro colado < 65 000
Hierro forjado 73 000
1
dU = σz εz dx dy dz
2
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 154
Por tanto, si el elemento de máquina es sometido a una tensión normal uniaxial, la energı́a
de deformación será
Z
σ
U= dv
V 2
Si el material se comporta de manera elástica lineal, puede utilizarse la ley de Hooke. La
energı́a entonces pasa a ser
σ2
Z
U= dv (8.2)
V 2E
1
dU = τ dx dy γ dz
2
Al integrar a todo el volumen, la energı́a de deformación que almacena la pieza es
Z
τγ
U= dv
V 2
τ2
Z
U= dv (8.3)
V 2G
El hecho de que la tensión sea uniforme en el sección implica que la forma de la sección no
influye en la resistencia o la rigidez ante este esfuerzo. Esto es totalmente cierto para el caso
de tracción, pero para la compresión solo es válido si las longitudes y/o las cargas no son tan
elevadas como para que aparezca pandeo, en cuyo caso sı́ que es muy importante la inercia de
la sección.
En realidad, la uniformidad de las tensiones solo es cierta si la sección se encuentra alejada
lo suficientemente del punto de aplicación del esfuerzo axil. Si este es aplicado de forma puntual,
la tensión no será constante en un área del orden del tamaño de la sección. Al ser la tensión
provocada por el axil uniforme en la sección, se podrı́a decir que este esfuerzo es el que mejor
aprovecha la sección, y por tanto, si es posible, se intentarán diseñar las barras para que trabajen
preferiblemente con este esfuerzo.
En una barra de longitud l sometida a axil, asumiendo la expresión (8.2), la energı́a de
deformación elástica es
P2
Z
U= 2
dv
v 2E A
pero si dv = A dx y siendo A constante, entonces
P 2V P 2l
U= =
2EA2 2EA
Figura 8.22: Barra sometida a un flector M. Las tensiones normales que provoca son lineales con la
coordenada perpendicular a la dirección del flector y para un flector My toman el valor σx = − M z 0
Izz y .
tensiones normales nulas. A esta zona se le llama superficie o lı́nea neutra. Para, por ejemplo,
My las tensiones normales serán lineales con la coordenada z de la sección. El valor de estas
tensiones es
Mz 0
σx = − y
Izz
donde la inercia Izz = A y 0 dx0 dy 0 .
R
Teniendo en cuenta la expresión anterior, para minimizar las tensiones se deberá maximizar
el valor de Izz . Además la tensión toma sus valores máximos en los extremos. Por esto los perfiles
adecuados para aprovechar el material cuando el esfuerzo flector es el dominante serı́an de la
forma que se muestra en la figura 8.23.
z 0
Para un momento flector, usando la expresión (8.2) y σx = − M
Izz y , la energı́a de una barra
es
M 2 y2
Z
U= dv
v 2 E I2
que teniendo en cuenta que dv = dA dx, entonces
l
M2
Z
U= dx
0 2E I
Por ejemplo, para una sección rectangular de dimensiones h × b, donde b(y 0 ) = b ∀y y donde el
1
valor de la inercia es Izz = 12 b h3 , entonces
V h2 V
→ τxy = τxy,max = 8yIzz = 32 Ay
Z h/2 2
y=0
Vy 0 0 Vy h 2
τxy = y Ab dy = −y
Izz Ab y 2 Izz 4 y = y
max → τxy = 0
Como vemos de los ejemplos anteriores, al contrario que el flector, el cortante provoca ten-
siones tangenciales τxy que toman valor cero en los extremos y son máximas en el centro de la
sección. Esto implica que la seccción que es buena para trabajar a cortante, no lo es para el
flector. Pero las tensiones que provoca el cortante son significativamente menores que aquellas
provocadas el flector y el torsor, por esta razón el perfil de la viga no se elige pensando en este
esfuerzo.
En cuento a la energı́a, si se aplica la expresión (8.3) a un perfil rectangular con altura h y
ancho b, esta es
3V 2 l
U=
5Gbh
Para una sección circular Ixx = J = 12 A r2 y por tanto, la tensión tangencial máxima, que
aparecerá en la superficie exterior, será
16T
τxθ,max =
πd3
Los perfiles más adecuados para trabajar a torsor son los perfiles huecos, ya que la tensión
tangencial circula de forma continua en toda la sección. Pero esto será ası́, y esto es importante,
siempre que la sección no se vea interrumpida, ya que si el perfirl hueco tiene una discontinuidad,
el circuito de tensiones tangenciales no podrá establecerse y las tensiones tomarán la disposición
mostrada en la figura 8.26b). Si ambas distribuciones, la cerrada y la abierta deben equilibrar un
mismo torsor T , entonces es fácil de inferir que las tensiones τxθ en el perfil de la figura 8.26a)
serán mucho menores que aquellas que aprareerı́an en el 8.26b).
A partir de la ley de Hooke y de la expresión (8.3), la energı́a elástica de una barra sometida
a torsor es:
l
T 2 r2 T2
Z Z
U= dv ⇒ U = dx
v 2GJ 2 0 2GJ
En la tabla 8.5 se recogen los valores de las tensiones normales y tangenciales máximas para
las secciones más comunmente usadas y en la tabla 8.6 la energı́a de deformación en barras para
cuatro tipos de cargas.
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 160
M V T
σmax = τmax = k τmax =
W A WP
W k WP
π d3 4 π d3
32 3 16
π d3 (1 − n4 ) 4(1 + n + n2 ) π d3 (1 − n4 )
32 3(1 + n2 ) 16
n=d0 /d
b h2 3 b h2
6 2 3 + 1.8n
n=b/h ≤1
B H 3 − b h3 V
τmax = 2 bm hm emin
6H Aalma
emin =mı́n[e1 , e2 ]
bm =(b+B)/2
hm =(h+H)/2
Cuadro 8.5: Tensiones máximas producidas por los esfuerzos flector, cortante y torsor para diferentes
secciones
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 161
P 2l P2
Rl
Axial P, E, A U= 2EA U= 0 2EA dx
M 2l M2
Rl
Flexionante M, E, I, U= 2EI U= 0 2EI dx
T 2l T2
Rl
Torsión T, G, J U= 2GJ U= 0 2GJ dx
3V 2 l 3V 2
Rl
Cortante V, G, A U= 5GA U= 0 5GA dx
(Sec. rectangular)
Solución
Con los dato de partida es fácil obtener el diagrama de esfuerzos flectores y cortantes. Este
es el que se muestra en la figura 8.28. De él se obtiene que las secciones más desfavorables son
aquellas adyacentes a las ruedas. De entre estas, aquellas secciones entre el apoyo del vagón y la
rueda además del flector de valor M = 40 kN m, también poseen constantes de v = 80 kN , por
tanto serán estas las más desfavorables. En esta sección se seleccionarán los puntos A y B por
ser estos donde se maximizan el flector y el cortante respectivamente.
En el punto A la tensión normal máxima debida al flector, usando la tabla 8.5 toma el valor
32M 32 · 40 × 103
σA = 3
= = 50.9 M P a
πd π · 0.23
siendo la tensión tangencial τA = 0. En el punto B se encuentra el máximo de la tensión
tangencial. Usando la misma tabla, esta es
4V 4 · 80 × 103
τB = = 2 = 3.4 M P a
3A 3 · π 0.24
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 162
Solución
Los esfuerzos torsor y cortante son constantes en el eje con las cargas de la figura, pero
en cambio, el flector es máximo en el extremo empotrado por tanto será la sección pegada
al empotramiento la más desfavorable. El valor de estos esfuerzos es en esa sección es M =
1000 · 0.05 = 50 N m, V = 1000 N y T = 1000 · 0.15 = 150 N m. En esta sección se estudiarán los
puntos A, B, y D situados donde muestra la figura 8.30.
En la sección de estudio, las tensiones máximas debidas al flector cortante y torsor son:
M 32M 32 · 50
σmax = 3
= = 63.7 M P a
πd π · 0.023
V 4 V 16 · 1000
τmax = = = 4.2 M P a
3 πd2 /4 3 · π · 0.022
T 16T 16 · 150
τmax = 3
= = 95.5 M P a
πd π · 0.023
En los puntos de estudio las tensiones son:
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 163
Como resultado:
σ A = 63.7 M P a τ A = −95.5 M P a
σ B = −63.7 M P a τ B = 95.5 M P a
σC = 0 τ C = −4.2 − 95.5 = −99.7 M P a
En la figura 8.31 se muestran las tensiones normales y tangenciales resultantes en los puntos
A, B y C. Las caras coloreadas son las caras externas y libres de tensiones.
de un agujero en una placa infinita es tres veces mayor a la tensión de la placa lejos del agu-
jero. Pero no todos los problemas de discontinuidades en materiales han podido ser resueltos
analı́ticamente. Estos problemas son tales como engrosamientos en ejes, ranuras, muescas, etc.
Un concentrador de tensiones es una parte de la geometrı́a de la pieza que altera la distribu-
ción de tensiones amplificándola. Se define el factor de concentración de tensiones Kt como
el factor que relaciona el esfuerzo máximo real en la discontinuidad con el esfuerzo promedio sin
la discontinuidad. Esto es
8.10.1. Gráficas de Kt
Los concentradores de tensiones que se suele apartar en la bibliografı́a más común en los
elementos de máquinas, como por ejemplo en [14] y [13], se suelen limitar a dos tipos de geo-
metrı́as, la placa plana y el eje macizo. Para otras geometrı́as pueden consultarse las referencias
[15] y [16]. Estas no son las únicas geoemtrı́as posibles, pero sı́ que son las más utilizadas. De
todas maneras, con la extensión del uso de los programas de cálculo numérico mediante elemen-
tos finitos (FEM), el uso de los concentradores de tensiones ha perdido relevancia. De todas
maneras su conocimiento sigue siendo pertinente ya que ayuda a interiorizar el fenómeno, ası́
como también puede servir de comprobación o número previo mediante FEM.
En la figuras 8.32 se muestra el concentrador de tensiones para un placa plana con tres tipos
de geometrı́as para los esfuerzos axil y flexionante. En la figura 8.33 se añaden los concentradores
de tensiones para una barra redonda y maciza con dos tipos de geometrı́a y para los esfuerzos
axil, flexionante y torsor.
De estas gráficas se puede inferir varias caracterı́sticas de la concentración de tensiones:
2. Está muy afectado por la geometrı́a. Puede apreciarse como a medida que disminuye el
radio asociado a la discontinuidad se incrementa la concentración de las tensiones.
Figura 8.32: Factores de concentración de tensiones para una placa con agujero, entalla y filete, para
los esfuerzos axil y flxionante. [17]
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 166
Figura 8.33: Factores de concentración de tensiones para una barra redonda y maciza con filete y ranura,
para los esfuerzos axil, flxionante y torsor. [17]
CAPÍTULO 8. ESFUERZOS Y TENSIONES EN BARRAS DE MECANISMOS 167
los resultados de las gráficas. Pero esta reorganización tensional evita la formación de grieta,
proporcionándole cierta reserva de resistencia. Por esto si la aplicación no requiere una excesiva
tolerancia dimensional y una gran rigidez, los concentradores de tensiones puede no ser tomados
totalmente en cuenta. Sin embargo esto no es posible en materiales frágiles donde el lı́mite
de plastificación coincide o se encuentra muy próximo al de rotura, siendo en este caso muy
necesario un correcto estudio de las concentraciones de tensiones.
Solución
La sección con el flector más importante es aquella donde actúa la fuerza F, pero en ese
punto no hay concentración de tensiones. Esta se encuentra en las secciones C y D. LA sección
C se encuentra más cercana al máximo de tensión y además el eje tiene ahı́ menos diámetro,
por lo que esta será la sección más desfavorable. Para encontrar el momento en C MC , primero
se calculará el valor de la reacción en FA y a partir de ella, el momento MC . Estos aparecen en
la figura 8.35.
b 0.325
FA = F = 6675 = 3556 N
a+b 0.25 + 0.035 + 0.1 + 0.225
MC = FA · AC = 3556 · 0.25 = 889.0 Nm
La tensión sin tener en cuenta el concentrador de tensiones la notaremos como una prima.
Ası́, la tensión máxima en la sección C sin el concentrador de tensiones será
0 32M 32 · 889.0
σmax = 3
= = 335 MPa
πd π 0.033
Para usar la tabla del concentrador de tensiones se calcularán los cocientes r/d = 3.2
30 = 0.1
38
y D/d = 30 = 1.27. Con estos datos y la gráfica del concentrador se obtiene kt = 1.95, tal y
como se muestra en la figura 8.35. El valor de la tensión máxima será, por tanto de
Teniendo en cuenta que la tensión de fluencia es de Sf = 600 MPa, se puede ver que se
sobrepasarı́a la tensión de fluencia y se tendrı́a plastificación en el eje.
Figura 8.35: Solución al ejemplo de uso de concentrador de tensiones en un eje biapoyado. Momentos
sobre el eje y consulta del valor del concentrador de tensiones para r/d = 3.2/30 = 0.1 y D/d = 38/30 =
1.27.
1. Obténgase una expresión para la energı́a de deformación que incluya las cargas actuantes
y una fuerza ficticia Q que actúa en el punto y dirección deseados
Solución
Puesto que la carga se aplica en la mitad de la vifga, el problema es simétrico. Debido a esto
puede calcularse la energı́a de media viga, obteniéndose la energı́a total multiplicando por 2. Las
cargas que se aplican son un flector M = P x/2 y un cortante constante V = P/2. Asumiendo
que la sección es rectangular y usando la tabla 8.6, la energı́a se expresa como
Z d/2 Z d/2 Z d/2 2 2 Z d/2
M2 3V 2 P x 3P 2
U =2 dx + 2 dx ⇒ U = 2 dx + 2 dx
0 2EI 0 5GA 0 4EI 0 10GA
Usando el Teorema de Castigliano y considerando que P, E, I, G y A no son funciones de x se
obtiene la flecha, que es
Z d/2 Z d/2
∂U P x2 3P
−δP = = dx + dx
∂P 0 8EI 0 10GA
Para una sección rectangular de sección A = b × h, I = bh3 /12, y por tanto
P d3 3P d
δP = − 3
−
|16Ebh
{z } |20Gbh
{z }
flector cortante
Para el acero G/E = 0.383, y en una viga, lo usual es que d/h sea al menos 10. Particularizando
para estos valores se obtiene que la parte de la flecha δP debido al flector es 162 veces superior
a la de cortante. Por esta razón, en la mayorı́a de las aplicaciones, a la hora del diseño, la flecha
del cortante suele despreciarse.
y
d
b
x
Q P Figura 8.37: Ejemplo 2, uso del teo-
rema de Castigliano. Fuerzas y reac-
Qb ciones aplicadas en una viga en vola-
Qd+P(d-b)
dizo con una carga P.
Solución
En este caso se generará una fuerza ficticia, puesto que ninguna fuerza actúa en el extremo
libre. El flector será
Qx para 0 < x < b
M=
Qx + P (x − b) para b < x < d
b Z d
M2 M2
Z
U= dx + dx
0 2EI b 2EI
Z b Z d
1 1
U= Q2 x2 dx + [Qx + P (x − b)]2 dx
2EI 0 2EI b
Z b Z d
1 2 2 1 2 2
Q x + P 2 (x − b)2 + 2QxP (x − b) dx
U= Q x dx +
2EI 0 2EI b
Fallo estático
9.1. Introducción
En este capı́tulo se estudiará el fallo de un elemento de máquina debido a la acción de las
cargas. Es evidente que las piezas se diseñen para que no fallen, pero conocer este estado último,
y cuándo y cómo se alcanza permitirá diseñar para que trabajen correctamente.
171
CAPÍTULO 9. FALLO ESTÁTICO 172
E: Impacto económico
En la tabla 9.1 se proporcionan los calores para Nsx para varias condiciones de A, B y C.
Para usar esta tabla se estima cada caracterı́stica para una determinada aplicación como muy
buena (mb), buena (B), regular (r) o pobre (p). En la tabla 9.2 se dan los valores de nsy para
las condiciones D y E. En esta tabla se entra con una estimación de cada caracterı́stica que se
valora como muy seria (ms), seria (S) o no seria (ns). Con ambas tablas se construye el factor
de seguridad.
Al factor de seguridad se le ha llegado a denominar factor de ignorancia porque en cierto
modo compensa la ignorancia de ciertas variables del problema. Por esto el método Pugsley
debe tomarse como meramente orientativo. Este método es además excesivamente conservador.
De hecho, en la normativa de ingenierı́a los factores de seguridad suele ser bastante mayores a
los que proporciona este método.
σ1 = Sf t
σ3 = Sf c
Figura 9.1: a) Zonas admisibles (en color sólido) y no admisibles (zona rallada) junto con ejemplos de
cı́rculos de Mohr de estados válidos (en lı́nea contı́nua), y fallidos (en lı́nea discontı́nua), utilizando el
criterio de la tensión máxima para plastificación en un estado de tensión bidimensional con σ3 = 0: a)
En el espacio de las tensiones normales y tangenciales. b) Para las tensiones principales σ1 y σ2 .
σ1 − σ2 sf
τmax = ≤ ⇒ σ1 − σ2 ≤ sf
2 2
En la figura 9.2 se representan las zonas admisibles (en color sólido) y no admisibles (zona
rallada) utilizando el criterio de la tensión tangencial máxima para plastificación en un estado
de tensión bidimensional con σ 6= 0, σ2 6= 0, σ3 = 0. A la izquierda en el espacio de las tensiones
normales y tangenciales, y a la derecha en el espacio de las tensiones principales σ1 y σ2 .
En el caso de la presión hidroestática, este criterio no pronosticarı́a plasticación, por lo que
en principio serı́a válido, y de hecho se utiliza en bastantes materiales para diferentes estados.
El problema para los metales reside en que la tensión de plastificación que pronostica no es la
correcta cuando se somete una probeta a un ensayo de cortante puro, (análogo al de tracción pero
aplicando un esfuerzo cortante en vez de uno axial). En un ensayo de cortante puro, el cı́rculo
de Mohr está centrado en el origen y por tanto σmax = τmax , siendo las tensiones principales
σ1 = −σ2 y con σ3 = 0. En este caso el criterio de tensión tangencial máxima pronosticarı́a
que la plastificación ocurrirı́a para una tensión tangencial τ = sf /2, valor que proviene de
σ1 = −σ2 = τ . Pero en realidad la plastificación sucede para τ = sf /1.73. Por esta razón se
buscaron nuevos criterios que pronosticaran el fallo por fluencia con mayor exactitud.
CAPÍTULO 9. FALLO ESTÁTICO 175
Figura 9.2: a) Zonas admisibles (en color sólido) y no admisibles (zona rallada) junto con ejemplos de
cı́rculos de Mohr de estados válidos (en lı́nea contı́nua), y fallidos (en lı́nea discontı́nua), utilizando el
criterio de la tensión tangencial máxima para plastificación en un estado de tensión bidimensional con
σ3 = 0: a) En el espacio de las tensiones normales y tangenciales. b) Para las tensiones principales σ1 y
σ2 .
1 ν
u= (σ12 + σ22 + σ32 ) − (σ1 σ2 + σ2 σ3 + σ1 σ3 ) (9.1)
2E E
De los resultados experimentales se desprendı́a que los materiales dúctiles sometidos a pre-
sión hidroestática tenı́an tensiones de fluencia muy superiores a las aparecidas en el ensayo de
tracción. De este hecho se podı́a intuir que la fluencia no estaba tan influenciada por la energı́a
asociada a los valores absolutos de las tensiones como por la asociada a la deformación angular
o distorisión del elemento tensionado. Por ello se dividió la energı́a total de deformación en dos
tipos, una energı́a de deformación de cambio de volumen, a la que se denominó uv y otra energı́a
de distorsión del elemento, pero sin cambio de volumen, notada como ud . De este modo
u = uv + ud
La energı́a de cambio de volumen está asociada a una tensión hidroestática que será la media
de las tensiones principales. De este modo se define
σ1 + σ2 + σ3
σm =
3
siendo uv la energı́a de deformación asociada al trabajo de esa tensión. Introduciendo σmed en
la expresión (9.1)
CAPÍTULO 9. FALLO ESTÁTICO 176
3σm2
uv = (1 − 2ν) (9.2)
2E
por lo que uv se obtendrá de la diferencia de las energı́as de las expresiones (9.1) y (9.2). Esta
expresión, es
1 + ν (σ1 − σ2 )2 + (σ2 − σ3 )2 + (σ1 − σ3 )2
ud = u − uv = (9.3)
3E 2
De este modo, la teorı́a de la energı́a de distorsión pronostica la falla por fluencia cuando
la energı́a de distorsión de una unidad de volumen sea igual a la energı́a de distorsión de una
misma cantidad de volumen de material sometida a la tensión que produce la fluencia en el
ensayo de tracción.
Durante el ensayo de tracción se tiene que σ1 = σ 0 y σ2 = σ3 = 0. Por lo tanto la energı́a de
distorsión de la unidad de volumen en el ensayo de tracción es
1 + ν 02
us = σ (9.4)
3E
Igualando las expresiones (9.3) y (9.4) se puede obtener el valor de σ 0 asociada a la energı́a
de distorsión, a esta tensión se la denominará tensión equivalente de Von Mises (σeq = σ 0 ),
que se define como
12
(σ1 − σ2 )2 + (σ2 − σ3 )2 + (σ1 − σ3 )2
σeq = (9.5)
2
de este modo el criterio de plastificación de la energı́a de distorsión, o criterio de Von Mises
anticipa la fluencia del material siempre que
σeq ≥ sf (9.6)
De la expresión 9.5 se desprende que cuando el elemento está sometido a una expresión
hidroestática de manera que σ1 = σ2 = σ3 , la tensión equivalente de Von Mises es nula, por lo
que cumplirı́a el criterio.
Para el caso de cortante puro, en este caso σ1 = −σ2 = τ , con σ3 = 0, usando la expresión
(9.5), se obtiene
2 4τ 2 + τ 2 + τ 2 τ
σeq = ⇒ σeq = √
2 3
Por lo que usando (9.6) el criterio pronosticaba que la plastificación ocurrirı́a para τ = √1 sf
3
tal y como habı́an avanzado los resultados experimentales.
En la superficie exterior sin carga de barras sometidos a los esfuerzos axil, cortante, flector y
torsor, que además es donde se dan las tensiones máximas, se tiene el estado tensional mostrado
en la figura 9.3a).
CAPÍTULO 9. FALLO ESTÁTICO 177
Figura 9.3: Circulo de Mohr en los puntos de la superficie de una barra sometida a esfuerzos axiles,
cortantes, flectores y torsores.
A partir de estos datos es sencillo obtener las tensiones principales σ1 y σ2 . Estas se obtienen
a partir del centro del cı́rculo de Mohr sumándole y restándole el valor del radio del cı́rculo. Este
se obtiene fácilmente de la figura, siendo
q
σ σ 2
r
σ 2 σ 1 = 2 + + τ2
R= + τ2 ⇒ q 2
2 σ = σ − σ 2
+ τ2
2 2 2
√
sf ≤ σ 2 + 3τ 2 (9.7)
9.4.1. Solución
El par entrante de Te = 150 Nm se convierte en un par de salida de
Z2 40
Ts = Te = 150 = 316 N m
Z1 19
En la figura 9.5 se muestra la fuerza de contacto ası́ como los radios base de la transmisión Rb1
y Rb2 y su relación con la distancia entre ejes. Donde dada la relación de transmisiónn = 1 : 2.105
se cumple Rb2 = 2.105 Rb1 .
Conocida la fuerza de contacto FC , el flector y el torsor en el punto medio entre apoyos son
Fc 5378 Fc
Mz 0 = e = 0.08 · = 215 N m Vy 0 = = 2689 N
2 2 2
Además del flector, también actúa el torsor Ts = 316 Nm y el cortante Vy0 = FC /2 = 2689 N.
En las direcciones del flector y el torsor se hace referencia a unos ejes x’y’z’. El eje y’ x’ coincide
CAPÍTULO 9. FALLO ESTÁTICO 179
con el eje 2, el eje y’ tiene la dirección de la fuerza de contacto Fc que viene determinada por el
ángulo de presión α y el eje z’ es perpendicular al anterior, como se muestra en la figura 9.5.
O
La sección más desfavorable estará en el eje 2 ya que flectores y cortantes son iguales para
O
ambos ejes, pero el torsor es más del doble mayor en el 2 . Dentro de ese eje, la sección con
mayores tensiones será la adyacente al piñon, pero por el lado derecho, ya que es la que sufre el
torsor.
En los puntos A y B se dan las tensiones normales máximas debidas al flector,y también las
tangenciales máximas debidas al torsor. En los puntos A y B la tensión debida al cortante es
nula. En el punto C se da el máximo valor de la tensión tangencial debida al cortante y debida
al torsor.
Además hay que tener en cuenta la concentración de tensiones. Los valores para entrar en
las tablas son dr = 302
= 0.067 y D 34
d 30 = 1.13. Con ellos se obtienen el concentrador para flector
KtF = 1.85, y para el torsor KtT = 1.4. Si no se disponen de tablas para el esfuerzo cortante, se
puede tomar el concentrador de torsor también para el cortante.
4V 4 2689
τV,max = KtT = 1.4 = 5.3 MPa
3A 3 π0.0152
32M 32 · 215
σ[Link] = KtF 3
= 1.85 = 150 MPa
πd π0.033
16T 16 · 316
τT,max = KtT 3 = 1.4 = 83 MPa
πd π0.033
Las tensiones equivalentes de Von Mises serán
p p
σeq,A = σeq , B = σ 2 + 3τ 2 = 1502 + 3 · 832 = 207 MPa
p
σeq,C = 3(5.3 + 83)2 = 153 MPa
CAPÍTULO 9. FALLO ESTÁTICO 180
495
n= = 2.38
207
Un coeficiente de seguridad mayor que uno quiere decir que la pieza no plastificarı́a. Este
coeficiente de seguridad es para la carga aplicada estáticamente. Si la es aplicada un grasn número
de veces este co3eficiente de seguridad se verá reducido. Concretamente, para este problema,
puesto que los ejes de la reductora giran, eso implica que la carga de flexión se para un mismo
punto del eje, pasa de compresión a tracción en cada vuelta. Esto generará problemas que se
verán en el siguiente capı́tulo.
Capı́tulo 10
Fatiga
10.1. Introducción
A finales del siglo XIX la ingenierı́a habı́a avanzado lo suficiente para poder diseñar máquinas
que trabajaban con fuerzas y velocidades cada vez mayores. También habı́a se habı́an desarro-
llado las herramientas para calcular las tensiones que sufrı́an los elementos de las máquinas
para poder ası́ compararlas con los de los ensayos del material. Se pensó que bastaba con no
superar las tensiones de plastificaban obtenidas experimentalmente. Pronto de comprobó que
esta prevención no serı́a suficiente.
En los primeros ferrocarriles se observó que los ejes se rompı́an con tensiones claramente
inferiores a las tensiones máximas obtenidas en ensayos de tracción. En esos primeros casos ya
se constató que las roturas se originaban en esquinas abruptas y ya desde esa fecha se comprobó
la utilidad del diseño con geometrı́as más fluidas y radios de acuerdo amplios. Las primeras
investigaciones exhaustivas sobre este fenómeno fueron realizadas por el ingeniero alemán August
Wöhler [19]. Mediante la realización de numerosos ensayos de laboratorio llegó a la conclusión
de que eran las cargas fluctuantes las que disminuı́an la vida de los ejes. Observó que la rotura
prematura dependı́a tanto del número de ciclos como de su amplitud.
Aunque en un primer momento se achacaron los fallos a la vibración, pronto se desestimó
esta teorı́a para introducir el fenómeno de la fatiga. Se ha avanzado mucho en el estudio del
antes extraño fenómeno que ocasionaba la rotura inesperada de los elementos de las máquinas,
pero todavı́a queda mucho camino por recorrer. Un siglo después de los trabajos de Wöhler los
primeros barcos fabricados mediante soldadura se partı́an abruptamente de forma aparentemente
inexplicable. En este caso fue la fatiga de las uniones soldadas las causantes de la catástrofe
[20]. Este ejemplo pone de manifiesto como la aparición de un nuevo procedimiento de unión,
como en este caso fue la soldadura, generó problemas con la fatiga desconocidos hasta la fecha.
Por tanto quedan todavı́a quedan muchas cuestiones abiertas en el estudio de la fatiga. Este
fenómeno está ı́ntimamente ligado a las caracterı́sticas de los materiales, y aparecen nuevos de
ellos continuamente. Las fibras de carbono es un ejemplo de gran importancia en la industria
aeronáutica moderna. El estudio de la fatiga en ellas es un campo todavı́a abierto en mucho de
sus aspectos.
La fatiga aparece tanto a grandes como a pequeñas escalas. Puede ser la causante de la
181
CAPÍTULO 10. FATIGA 182
rotura prematura tanto de un pequeño rodamiento como de una barco completo. En ambos casos
concluye con el colapso por unas cargas en principio asumibles. Una caracterı́stica importante y
peligrosa es que el fallo por fatiga se asemeja a una rotura frágil, aun tratándose de materiales
dúctiles.
Otro aspecto a tener en cuenta es que el diseño a fatiga no se realiza en el peor de los
escenarios posibles, como por ejemplo con la superficie en muy mal estado, la vida a fatiga
tendrı́a valores tan bajos que llevarı́a a un diseño excesivamente conservador, o incluso imposible.
Por ello este enfoque se abandona para la fatiga. El peor de los casos se deshecha para pasar a
diseñar en un estado mejor, pero siempre con el cuidado de evitar escenarios desfavorables.
Como se puede apreciar, este es un fenómeno complejo. Esencialmente consiste en la gene-
ración de pequeñas grietas que progresan rápidamente cuando adquieren un tamaño crı́tico. La
fatiga es algo que un ingeniero siempre debe tener en cuenta en presencia de esfuerzos cı́clicos de
cierta entidad. Las grietas usualmente comienzan en el exterior de los ejes y barras, que son los
lugares donde se encuentran las tensiones máximas generados por los esfuerzos flector y torsor.
También en esta zona suelen encontrarse discontinuidades en la geometrı́a con concentración
de tensiones. Teniendo en cuenta los fenómenos de la generación y la propagación de la grieta,
se considera que la vida del componente es la suma de los tiempos necesarios de aparición y
desarrollo de las grietas que derivan en rotura.
Por tanto, la vida del componente se maximizará siguiendo unas pautas:
3. Maximizado del tiempo de propagación. Las grietas se propagan más rápidamente por las
zonas entre granos, por lo que ciertas estructuras metalográficas son más propensas a una
rápida propagación de las grieta que otras.
Resulta obvio que las operaciones para aumentar la vida a fatiga tienen costes. Pero también
la minimización de las tensiones exige mayor cantidad de material, lo que aumentará el peso de
la máquina, factor crı́tico en la industria aeronáutica. Por todo esto, como suele ser habitual
en la ingenierı́a, el diseño a fatiga es un proceso con complicados compromisos entre el peso, la
duración y el coste de los componentes.
CAPÍTULO 10. FATIGA 183
σmax + σmin
σm =
2
El rango de tensiones σr es la diferencia entre la tensión máxima y la mı́nima, esto es
σr = σmax − σmin
σmax − σmin
σa =
2
La relación de esfuerzo RS es la relación de las tensiones mı́nima y máxima:
σmin
Rs =
σmax
Y la relación de amplitud resulta del cociente entre la amplitud de esfuerzo y el esfuerzo
medio, y por tanto
σa σmax − σmin 1 − Rs
Aa = = =
σm σmax + σmin 1 + Rs
Ciertos tipos de tensiones cı́clicas tienen una denominación especial, estos son algunos ejem-
plos:
ferrosos presentan una zona en la que por debajo de una determinada carga no ocurren fallos
independientemente del número de ciclos. La carga que admite un número de ciclos infinitos
en el ensayo de vida rotatoria se le denomina lı́mite de fatiga que se escribirá como Sn0 . Con
la prima se denota que es para el ensayo de viga rotatoria. Al lı́mite a fatiga en una situación
general se le retirará la prima para ser simplemente Sn . Este lı́mite a fatiga suele encontrarse
alrededor del 50 % de la tensión última del material para esfuerzos de flexión. La figura 10.4
muestra los lı́mites a fatiga en función de la tensión última (Su ) para diferentes aceros forjados.
Las aleaciones metálicas no ferrosas, como las de aluminio, no presentan esta zona, ası́ que
las cargas, por bajas que sean, siempre fatigan el material y por tanto nunca tendrán vida
infinita. En la figura 10.5 se muestra el comportamiento a fatiga de las aleaciones de aluminio.
No obstante estos materiales presentan una inclinación menor en las curvas S-N para N > 106
que para 103 < N < 106 . En la tabla 10.1 se muestran las tensiones para N = 108 lo que en
estos materiales suele considerarse como vida infinita.
La fatiga no es un fenómeno exclusivo de metales. En la figura 10.6 se representa el com-
portamiento a fatiga de diferentes polı́meros. Pero dado que los polı́meros tienen una gran
variabilidad en sus caracterı́sticas mecánicas, esta tabla debe tomarse a efecto indicativo de su
CAPÍTULO 10. FATIGA 186
Material N Tensión
Aleaciones de Mg 108 Sn0 = 0.35 Su
Aleaciones de Cu 108 0.25 Su < Sn0 < 0.5 Su
Aleaciones de Ni 108 0.35 Su < Sn0 < 0.5 Su
Titanio 108 0.45 Su < Sn0 < 0.65 Su
Aleaciones de Al 5 × 108 Sn0 = 0.45 Su (Su < 330 M P a)
Sn0 = 130 M P a(Su ≥ 330 M P a)
Cuadro 10.1: Lı́mite a fatiga aproximado para aleaciones metálicas no ferrosas [21]
industrial o aeronáutica. Por esto, la fatiga que se contemplará en lo sucesivo será únicamente
la alto número de ciclos.
El primer paso para la resolución del ejemplo consiste en la reconstrucción de la curva S-N.
Esta se muestra en la figura 10.8. Esta tarea es muy sencilla ya que se conocen las tensiones
para 103 y 106 ciclos.
La pendiente de la recta entre 103 < N < 106 será:
log(0.9 Su ) − log(0.5 Su )
m=
6−3
log(1080) − log(600) 3.03 − 2.78
m= 6 3
= = 0.085
log(10 ) − log(10 ) 3
Con ella y sabiendo que el punto 0.5 Su para N = 106 está incluida en ella el cálculo de la
CAPÍTULO 10. FATIGA 188
Para el cálculo del número de ciclos N que se esperan soportar para σ = 750 M P a, la
ecuación será:
log(750) − log(600)
m = 0.085 =
6 − log(N )
log(750) − log(600)
log(N ) = 6 − = 4.86
0.085
N = 104.86 = 72400 ciclos
Sn = Sn0 · CG · CS · CR · ...
CAPÍTULO 10. FATIGA 189
CG diámetro d
1 d < 0.30 o d ≤ 8 mm
0.869 · d−0.112 0.30 < d ≤ 100
1.189 · d−0.112 8 mm < d ≤ 250 mm
Cuadro 10.2: Factor de gradiente CG para diferentes ejes en función de su diámetro [13]
encuentran en el exterior sufren muy bajas temperaturas cuando la aeronave vuela a grandes
alturas, lo que conlleva la fragilización del material. Llegado a este punto el diseñador puede, o
bien utilizar los datos corregidos para esas temperaturas, o bien usar el factor de temperatura
CT .
Probabiliad de supervivencia ( %) CR
50 1.00
90 0.90
95 0.87
99 0.82
99.9 0.75
99.99 0.70
Figura 10.14: Ejemplo de uso de factores reductores de Sn . Probeta rectangular sometida a axil fluc-
tuante.
La carga última se calcula multiplicando la tensión última Su por el área de la sección, que
es de 30 × 30 mm.
y frena el desarrollo de la grieta. Por tanto, el valor de Kt deberá ser moderado para su uso en
fatiga.
En lugar del factor de concentración de esfuerzos convencional Kt , se usará el factor de
concentración de esfuerzos Kf , definido como
S˜n
Kf =
Sn
donde S˜n es el lı́mite a fatiga para el espécimen con muesca y Sn el lı́mite para el espécimen sin
muesca.
A diferencia del Kt que solo dependı́a de la geometrı́a y del tipo de esfuerzo, Kf sı́ que depende
del material. La relación entre Kt y Kf se establece a través del factor de sensibilidad a la
entalla qn , y toma la forma
Kf = 1 + (Kt − 1) qn ; qn ∈ [0, 1]
Pu = 0.032 · Su = 673 kN
Pf P
σf = 427, 5 M P a = Kt = 2.5 ⇒ Pf = 154 KN
A 0.032
Para el cálculo del lı́mite a fatiga para axil alternante se utilizarán los factores de reducción
de Sn0 de gradiente y el de superficie con valores CG = 0.8 y CS = 0.74. De este modo el lı́mite
a fatiga será:
748
Sn = Sn0 · CG · CS = 0.8 · 0.74 = 221 M P a
2
La sensibilidad a la entalla aparece en la figura 10.15. Para una tensión última de Su =
748 M P a y r = 2.5mm el valor es qn ' 0.86.
Sn
Pn = A = 87 kN
kf
Como es habitual en fatiga, la manera en que se abordó la fatiga con tensiones tangencial fue
experimental. Para ello se realizaron ensayos para especı́menes sometidos a tensiones normales
y tangenciales, obteniéndose los lı́mites a fatiga para ambas situaciones, Sn y Tn . Los resultados
se muestran en la tabla 10.4.
Material Tn /Sn
Acero 0.60
Aleaciones de Al 0.55
Cu y aleaciones de Cu 0.56
Aleaciones de Mg 0.54
Titanio 0.48
Fundición de Fe 0.90
Fundición de Al y Mg 0.85
Cuadro 10.4: Cociente entre los lı́mites a fatiga τn /Sn para diferentes materiales [25]
Para estos esfuerzos el inicio de la grieta también se produce por la acumulación de disloca-
ciones en las bandas de máxima tensión tangencial, como sucede con el axil cı́clico o la flexión
rotatoria. Esto llevó a pensar en utilizar el criterio de Tresca de la tensión tangencial máxima
para incluir las tensiones tangenciales. Pero a excepción del titanio, el cociente entre ambos
√
resultó acercarse más al predicho por el criterio de Von Mises (0.577 = 1/ 3).
También es remarcable el resultado de Tn /Sn para las aleaciones fundidas. Pero se puede
considerar que las escamas de grafito de la fundición gris actúan como muescas iniciales para las
microfisuras. Estas están aleatoriamente orientadas, lo que lleva a que la tensión máxima en los
extremos de las escamas tendrá un carácter equivalente para tensiones normales y tangenciales,
y por tanto implicarı́a Tn = Sn .
En realidad, las razones fı́sicas por las que aplicar el criterio de plastificación de Von Mi-
ses en la fatiga son bastantes cuestionables, sobre todo teniendo en cuenta la importancia del
deslizamiento cı́clico de en la nucleación de las grietas. Pero los datos experimentales son más
poderosos que las razones microestructurales. Por esta razón se usa la tensión equivalente de
Von Mises para tratar barras sometidas a esfuerzo que provocan tanto tensiones normales como
tangenciales.
Al no disponerse de más datos del material, visitando la página Matweb [28] y utilizando
su herramienta de búsqueda con S355J2, aparecen algunos aceros. Se elige EN S355J2G4 High
Manganese, Structural, Hot Rolled, Quality Steel y aparece la tabla 10.5
Cuadro 10.5: Propiedades del acero EN S355J2G4 High Manganese, Structural, Hot Rolled, Quality
Steel según la página web [28]
Su 490 × 106
Sn = CG CS = 0.85 · 0.78 · = 162 M P a
2 2
Los torsores admisibles son aquellos que producen una tensión equivalente tal que
FS σeq ≤ Sn
CAPÍTULO 10. FATIGA 198
Estos esfuerzos producen tensión tangencial de valor τ = 16T /(πd3 ). Usando el criterio de
Von Mises la tensión equivalente sera:
√ 16 T
q
2
σeq = Kt @ σ@ + 3τ 2 = Kt 3 3
πd
Pero para el lı́mite a fatiga se deberá usar el concentrador Kf = 1 + (kt − 1)qn . De las tablas
se obtiene que Kt = 1.54, siendo la sensibilidad a la entalla q = 0.76 para Su = 890 M P a, radio
de entalla r = 2 mm y esfuerzo de tensión. De aquı́ se tiene que Kf = 1.4.
El valor máximo de T se obtiene igualando
√ 16 T
Sn = FS · Kf 3 3
πd
resultando
n1 n2
+ ≥1
N1 N2
Y en general para nC bloques de cara
nC
X ni
≥1 (10.1)
Ni
i=1
Si NT es el número total de ciclos para llegar al fallo, se pueden definir unos ratios αi que
corresponderı́an a la proporción en la que aparece la carga i respecto al total de cargas aplicadas.
Es decir αi = ni /NT . Entonces
ni = αi NT (10.2)
nC
X αi 1
≥ (10.3)
Ni NT
i=1
α1 α2 α3 α4 1
+ + + =
N1 N2 N3 N4 NT
y sustituyendo, se obtiene el número total de ciclos que soportará la barra. Obsérvese como las
cargas por debajo del lı́mite de fatiga no contribuirı́an al daño.
σa σm
+ =1
Sn Su
σa σm
+ =1
Sn Sf
σa + σm = Sf
Para definir el coeficiente de seguridad de una determinada carga cı́clica con sus correspon-
dientes σm y σa se utilizará otro punto que posea la misma relación de amplitud e incluido en la
lı́nea del criterio. Para el punto A = (σmA , σ A , ) ese punto serı́a el Ag = (n σ A , n σ A , ) mostrado
a g m g a
en la figura 10.21 y donde ng es el coeficiente de seguridad para el criterio de Goodman. Puesto
que el punto AG cumple la ecuación (10.4) entonces
ng σa ng σm
+ =1 (10.4)
Sn Su
y por tanto
1
ng = σa σm
Sn + Su
Sf Sf
nf = o nf =
σa + σm σmax
Como se ha comentado, el criterio de Goodman modificado añade una comprobación de plas-
tificación, al criterio de Goodman. El coeficiente de seguridad, por tanto, será n = mı́n(ng , nf ).
En la zona rallada en horizontal en la figura 10.21 se cumplirá que ng < nf y en la rallada en
horizontal nf < ng .
CAPÍTULO 10. FATIGA 203
2. Calcule el coeficiente a fatiga y la vida en numero de ciclos si el acero pasa a ser un AISI
1015 de Su = 386.1 M P a y Sf = 284.4 M P a
Figura 10.23: Ejemplo de cálculo a fatiga para barras sometida a esfuerzos cı́clicos con tensiones medias
diferentes de cero.
Su
S n = CG CS = 315 M P a
2
En este problema el flector fluctúa y el torsor es constante. Las tensiones normales provocadas
por el flector serán las que aportarán la parte alternante. La tensión tangencial provocada por
el flector dará la componente media de la excitación. El valor de las tensiones provocadas por
flector y torsor sin tener en cuenta el concentrador de tensiones son:
32 M 32 · 259.5
σ= 3
= = 97.9 M P a
πd π 0.033
16 T 16 · 300
τ= = = 56.6 M P a
π d3 π 0.033
CAPÍTULO 10. FATIGA 204
Para el cálculo del coeficiente de seguridad según el criterio de Goodman modificado primero
se calculará el coefciente de seguridad frente a la plastificción.
p p
σeq,max = σ 2 + 3 τ 2 = (1.85 · 97.9)2 + 3(1.4 · 56.6)2 = 227 M P a
siendo nf
Sf 495
nf = = = 2.2
σeq,max 227
Puesto que en este problema se presentan tanto tensiones normales como tangenciales deberá
que usarse tensiones equivalentes. Para la comprobación de plastificación no hay dudas de que
se debe usar Von Mises, pero para la comprobación de nG no está tan claro. Hay diferentes
propuestas dependiendo del autor sobre cómo calcular esas tensiones equivalentes. En el cociente
σm
Su se está comparando la carga media (no la alterna) con una tensión última. En este caso
muchos autores piensan que no deberı́a usarse una tensión equivalente de plastificación sino una
para rotura, como podrı́a ser la de tensión máxima.
La tensión equivalente con un criterio de rotura será la tensión máxima, es decir la tension
principal σ1 (recuérdese que se toma σ1 > σ2 > σ3 ). Esta será notada como σeq∗ y para un estado
1 1
ng = σeq,a ∗
σeq,m = 1.70×97.9 56.6 = 1.68
+ 315 + 1000
Sn Su
Puesto que ng < nf el coeficiente dee seguridad será n = ng = 1.68.
1 1
ng = σeq,a ∗
σeq,m = 1.54×97.9 56.6 = 0.81 < 1
+ 139.3 + 386.1
Sn Su
Por lo que el criterio de Goodman no se cumple y la pieza no tendrı́a vida infinita. Para
estimar el número de ciclos deberemos extrapolar el punto en el espacio de las tensiones medias
y alternas al eje de las tensiones alternas. Es decir, se deberá calcular una tensión alterna
modificada σ ca que con la misma vida que la pieza sometida a σm y σa . Con esta tensión alterna
pura ya sı́ se podrı́a entrar en la gráfica S-N para calcular la vida del espécimen.
Para el cálculo de esta tensión σ
ca se prolongará una lı́nea desde el punto Su en el eje de σm
hasta la intersección con en el eje de σa , tal y como se muestra en la figura 10.25
Para ello se tendrá en cuenta la proporcionalidad que se puede establecer por la linea de
puntos en la parte izquierda de la figura 10.25.
σa σ
ca σa
= ⇒ σ
ca =
Su − σm Su 1 − σSmu
que para el presente ejemplo se convierte en
Figura 10.25: Cálculo del número de ciclos para esfuerzos cı́clicos con σm 6= 0.
σa ) − log(Sn )
log(c log(0.9 Su ) − log(Sn )
6
=
log(10 ) − log(N ) log(106 ) − log(103 )
2.25 − 2.14 2.54 − 2.14
=
6 − log(N ) 6−3
⇒ N = 150 000 ciclos
Capı́tulo 11
11.1. Problema
El sistema con husillo de la figura 11.1 está realizado con un acero de tensión de plastificación
de sf = 330 M P a, tensión última su = 450 M P a y lı́mite a fatiga (sin concentradores de de
tensión) de sn = 180 M P a. El rozamiento entre superficies de acero es µ = 0.2. El sistema esta
diseñado para levantar una compuerta unida al actuador con una barra a 45◦ . El husillo tiene
un filete trapezoidal de diámetros mı́nimo y máximo de 20 y 25 mm respectivamente siendo
paso p = 5 mm. El ángulo de entre lados del trapecio es de β = 30◦ . En los puntos A y B los
concentradores de tensión para cortante y torsor son ktV = ktT = 1.4 y de ktA = ktM = 1.6 para
axil y flector. La sensibilidad a la entalla es de q = 0.8. El motor tiene un par de Tm = 6 N m y
la reducción en la transmisión de engranajes es r = 2 (i = 1 : 2), con ángulo de presión α = 30◦ .
En el punto C hay una articulación que permite el giro.
1. Calcule el máximo peso de compuerta que puede alzarse con el sistema. Si se incrementa
el paso un 25 % respecto al máximo para que la puerta no caiga sola, calcule el nuevo par
motor y el freno necesario para frenar la compuerta.
207
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 208
∅max + ∅min 20 + 25
Rm = = = 11.25 mm
2·2 4
p
α = tan−1 = 4.05◦
2π · Rm
µ
µ0 = = = 0.207
cos β/2
φ = tan−1 µ0 = 11.70◦
Thusillo = Tm · r = Tm /i = 6 · 2 = 12 N m
Con el paso igual a p0 = 1.25 · pr = 18.28 mm el nuevo ángulo de la hélice del husillo sera:
p0
α0 = tan−1 = 14.50◦
2π · Rm
y el nuevo par será, por tanto
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 209
0
Thusillo = Mc · g · Rm · tan(α0 + φ) = 21 N m
0 0
Tm = Thusillo /2 = 10.5 N m
y el par de frenada en el husillo que serı́a necesario proporcionar serı́a
Tf,m = Tf,husillo /r = 1 N m
En el lado izquierdo de la figura 11.2 se muestran las reacciones con las flechas en gris y las
fuerzas externas con las flechas en blanco. Los esfuerzos son las flechas en negro dibujadas en
una sección genérica en medio de la barra.
El punto más desfavorable de la viga es el B. En él aparece un flector M B = 0.1 · Mc g =
√
392 N m. Las secciones también sufren axil N B = 400 · 9.8/ 2 = 2772 N y cortante V = 2772 N ,
de manera que su suma sea igual a la reacción RB , VB + NB = RB = Mc g.
El cortante produce su máximo en la linea neutra, a diferencia del flector que lo proporciona
en los puntos más alejados de ella, por esto y porque no es un perfil de alma delgada, no se
computarán su tensiones. Ası́ que para el diseño se contará con el flector y el axil.
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 210
Si los diámetros exterior e interior del tuvo son de y di , entonces tensiones máximas provo-
cadas por el flector M B y el axil N B en un tubo redondo hueco son:
F 32 M 1
σmax = KtF 4
π d3e di
1− de
N 4N
σmax = KtN
π d2e − π d2i
Y puesto que di = de −2 e, donde e = 3×10−3 m, para obtener el diámetro del tubo necesario
se resolverı́a la ecuación
F N
sf = σmax + σmax
32 · 392 · 1.6 1 4 · 2772 · 1.6
330 × 106 = 4 +
π d3e π d2e − π(de − 0.003)2
de −0.003
1− de
La solución a esta ecuación es de = 43.8 mm. Este diámetro de tubo no está disponible en
los almacenes de acero ası́ que se usarı́a el más cercano disponible, que por ejemplo podrı́a ser
de = 50 mm. Si este fuese un problema de examen y no tuviéramos acceso a un ordenador,
la solución se encontrarı́a probando, sabiendo que tampoco se tendrı́an que realizar muchas
iteraciones porque éstas serı́an con los diámetros disponibles. De este modo se comprobarı́a que
para de = 40 mm ⇒ σmax B > sf y para de = 50 mm ⇒ σmax B < sf
sf sf d3
nf = √ =q
σ2 + τ 2 32 M F 2
16 T T 2
π kt +3 π kt
Y por tanto
1
ng =
1.48 · 32 M 16 T
+
π d3A sn π d3A su
Sustituyendo y despejando se obtiene el otro dA .
d3A 1/3
2= ⇒ dA = 2 · 2.22 × 10−6 = 0.0164 m
1.48 · 32 · 24.9 16 · 12
+
π · 180 × 106 π · 450 × 106
El diámetro elegido para garantizar el coeficiente de seguridad de 2 será el mayor de ellos
dA = 16.4 mm.
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 212
11.2. Problema
El conjunto de la figura 11.4 está realizado con un acero de tensión de plastificación de
sf = 500 M P a, tensión última su = 650 M P a y lı́mite a fatiga (sin concentradores de de tensión)
O
de sn = 227.5 M P a. En el punto A el concentrador de tensión para flector es KtF = 1.4,
O
y para torsor KtT = 1.1. En el punto B hay una masa de 40 kg. El sistema gira con un
movimiento armónico entre los valores ±π/12 respecto a la referencia horizontal. La gravedad
actúa verticalmente. La sensibilidad a la entalla es q = 0.8 tanto para esfuerzos de flector como
de torsor.
2. Calcule la máxima frecuencia a la que el sistema puede oscilar sin que falle por fatiga.
32 · M 32 · 40 · 9.8 · 0.3
σ = KtF F 3
= 1.4 = 62.1 × 106 Pa
π·d π · 0.033
16 · T 16 · 40 · 9.8 · 0.3
τ = KtT 3
= 1.1 = 24.4 × 106 Pa
π·d π · 0.033
En esta cuenta se ha considerado que el ángulo π/12 es pequeño por lo que el coseno de este
ángulo se tomará como 1. El coeficiente de seguridad será
sf 500 × 106
n= √ =√ ⇒ n = 6.66
σ2 + 3 τ 2 62.12 + 3 · 24.42 × 106
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 213
Figura 11.5: Alzado de sistema oscilante del Ejemplo 3 con las fuerzas de inercia y el peso.
π
AI = −0.3 · ω 2 sin ωt
12
π I
F I = 40 · 0.3 · ω 2 sin ωt ⇒ Fmax = ω2 π
12
Puesto que existe un problema de fatiga, habrá que calcular los concentradores de tensión
kfF y kfT . Estos son
32 · ω 2
π · 0.3
σa = 1.32 = 0.469 × 106 ω 2 Pa
π · 0.033
16 · ω 2
π · 0.3
τa = 1.08 = 0.192 × 106 ω 2 Pa
π · 0.033 p
p
σeq,a = σa2 + 3 τa2 = ω 2 0.4692 + 3 · 0.1922 × 106 Pa = 0.575 × 106 ω 2 Pa
∗ 32 · M 32 · 40 · 9.8 · 0.3
σm = = = 44.4 × 106 Pa
π · d3 π · 0.033
∗ 16 · T 16 · 40 · 9.8 · 0.3
τm = 3
= = 22.2 × 106 Pa
π·d q π · 0.033
∗
σm ∗
∗
σeq,m = 2 + ( σ2m )2 + (τm
∗ )2 = 58.7 × 106 Pa
1 1
ng = σeq,a ∗
σeq,m = 0.575×106 ω 2 58.7×106
= 1 ⇒ ω = 19.0 rad/s ⇒ f = 3.02 Hz
sn + su 227.5×106
+ 650×106
Respecto a nf
11.3. Problema
O
En el sistema de la figura, el actuador 3 provoca una flecha de valor δ = 0.07 m en el
O
extremo de la viga empotrada de sección circular 2 . El actuador genera el movimiento mediante
un tornillo de diámetro exterior 24 mm, diámetro interior 20 mm, paso 4 mm y ángulo entre
flacos de 30o . El tornillo va unido al motor con una reductora de ı́ndice de reducción i = 10.11.
Los concentradores de tensiones en el empotramiento de la viga son de KtF = 1.6 y KtT = 1.4.
El momento de inercia del área es I = π4 r4 (W = Ir ). Datos: sf = 500 M P a, su = 650 M P a,
sn = 230 M P a (sin concentradores de tensión), E = 210 GP a, q = 0.8, µ = 0.2.
1. Calcule el diámetro máximo de la viga para que no haya plastificación, con un coeficiente
de seguridad de 2.
2. Calcule el par del motor necesario para generar la flecha mencionada con una barra de 60
mm de diámetro. Calcule el paso de tornillo máximo para que el sistema mantenga esa
flecha sin aporte de par del motor. Calcule el par motor para ese nuevo paso.
3. Si se tiene una viga de 80mm de diámetro y el actuador flecta la viga hasta δ = 0.07 m y
la relaja alternativamente a una frecuencia relativamente baja. Calcule la flecha máxima
en la viga para que este proceso se pueda realizar indefinidamente sin que la viga rompa
por fatiga.
1.5m
4mm
30º
3
1.5m
� 20mm
3m 2 24mm
Fv L3 Fv L3 Fv 1 4L3
δ= = π 4 = 3 (11.1)
3EI 3E 4 r r r 3Eπ
Para calcular el diámetro máximo de la viga se limitarán las tensiones máximas a la mitad de
la tensión de plastificación. Téngase en cuenta que cuanto mayor sea el diámetro, menos flecha
aceptará la viga sin plastificar. Por otro lado, al estar el actuador inclinado a 45o , aparece una
Fh
fuerza horizontal Fh = Fv . Esta fuerza generará axil, y éste, tensiones normales iguales a πr 2.
Pero estas tensiones serán mucho menores que las provocadas por el flector, por lo que al elegir
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 216
el diámetro de la viga, redondearemos hacia abajo y comprobaremos la sección con todos los
esfuerzos.
sf 32M 4L Fv
σmax = = KtF = KtF (11.2)
2 πd3 π r3
donde el momento en el empotramiento es M = Fv L, siendo L la longitud de la viga (L = 3 m).
Si de la expresión 11.2 se despeja Fr3v y se sustituye en 11.1, esta ecación toma la forma
1 4L3
sf π
δ=
2KtF 4L r 3Eπ
donde todo es conocido salvo el radio de la sección de la viga (r ). Por tanto
sf L2 500 × 106 · 32
r= = = 0.0319 m
δ2kKtF 3E 0.07 · 2 · 1.6 · 3 · 2.1 × 1011
Con este resultado se elige una sección de 60 mm de diámetro y se comprueban las tensiones
en la sección con el flector y el axil. No hace falta usar tensiones equivalentes de Von Mises
porque no hay tensiones tangenciales en el punto de comprobación, el extremos superior de la
sección, donde se suman las tensiones de flector y axil (en el extremo inferior se restarı́an). Las
tensiones provocadas por el cortante no se tendrán en cuenta ya que estas actúan en el centro
de la sección, donde solo hay tensión normal, que es muy inferior a la tensión provocada por el
flector.
Para calcular los esfuerzos primero se obtendrá Fv de la ecuación 11.1
De este modo
4Fv L N 4 · 1039 · 3 1039 500
σ = KtF 3
+ KtF 2 = 1.6 3
+ 1.6 2
= 236 M P a < MPa
πr πr π0.03 π0.03 2
∗
Thus = Fhus rm tan(2ϕ) = 1469 · 0.011 · tan(2 · 11.69◦ ) = 6.99 Nm
1 4.6914 × 104 4 · 3
∗ 1 4M N
σm = + 2 = + 1 δ = 1.4047 × 109 · δ
2 πr3 πr 2 0.042 π 0.04
Para la tensión alternante sı́ que tendrán que añadirse el concentrador de tensión, pero en este
caso se disminuirá teniendo en cuenta la sensibilidad a la entalla, por lo que se usará el valor
Por lo que la flecha máxima será la mı́nima de las dos flechas calculadas anteriormente
11.4. Problema
La culata del compresor que se representa en la figura 11.6 está fijada por 6 tornillos de
diámetro medio dm = 10 mm, α = 2.5◦ , β = 30◦ , µ = 0.8. Se propone para la junta un material
de módulo de Young E J = 30 × 109 Pa. El resto del material es un acero de E = 210 × 109 Pa.
1. Determine el par de apriete en los tornillos para que la junta quede sellada cuando el com-
presor esté en funcionamiento de manera que la relajación en la fuerza en la junta sea igual
a la fuerza que quede remanente en la junta cuando el sistema entre en funcionamiento.
Calcular la fuerza qe debe soportar el tornillo en esta situación.
2. Determine el nuevo par de apriete si se desea que queden 5000 N de fuerza en la junta
cuando el compresor esté funcionando. Calcule la nueva fuerza que está soportando el
tornillo.
18.9 mm
10
3
12 10
p = 100 bar
A = 80 cm2
8.5
11.5
Se denominará Ptotal a la fuerza total en la culata que deberá ser soportada por los 6 tornillos.
Por tanto, la fuerza en un tornillo será P = Ptotal /6.
80000
Ptotal = p · A = 100 × 105 · 80 × 10−4 = 80 kN ⇒ P = = 13.33 kN
6
Los incrementos de fuerza en el tornillo y en la placa que se relaja, es decir en la junta, son:
kt0 kr
∆Ft = P ∆Fpr = P
kt0 + kpr kt0 + kpr
Para evaluar estos incrementos es necesario calcular las rigideces kt0 y kr . Para ello se hará
uso de la gráfica 11.7. En ella se muestran las porciones de bridas y de junta que se computan
para obtener la rigidez de las distintas placas.
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 220
32.2
30.4
18.9
12
10
30º
3
10
24.7
31.3
La rigidez kt0 engloba dos rigideces en serie, la propia del tornillo y la de las dos placas que
se comprimen, que son las bridas de acero. De este modo
1 2 1
0 = +
kt kpc kt
donde kp c es la rigidez de una de las placas que se comprimen, que son las de acero, y kt es
la rigidez del tornillo. La rigidez del tornillo es la de una barra de acero de longitud 23 mm y
sección circular de diámetro, el diámetro medio del tornillo (dm ). Por tanto
2 2
EA E · π d4m 2.1 × 1011 · π 0.01
kt = = = 4
= 7.17 × 108 N/m
L L 0.023
La rigidez de la placa de acero será la de un cilindro hueco de diámetro interior 12 mm y
exterior 24.7 mm, que es el diámetro medio de la zona en gris más oscura de la brida de acero,
tal y como se muestra en la figura 11.7
2 d2i 2 0.0122
E · π( d4e − 4 ) 2.1 × 1011 · π( 0.0247 − 4 )
kpc = = 4
= 3.34 × 109 N/m
L 0.023
Por tanto
1 2 1
0 = 9
+ ⇒ kt0 = 5.02 × 108 N/m
kt 3.34 × 10 7.17 × 108
Queda por calcular la rigidez de la junta, que es la única placa que se relaja (kpr ). Esta es
la de un cilindro de sección hueca con diámetros 12 y 31.3 mm, siendo la altura 3 mm.
2 0.0122
30 × 109 · π( 0.0313 + 4 )
kpr = 4
= 6.56 × 109 N/m
0.003
Con los resultados anteriores ya se puede calcular el la variación de fuerza de compresión en
la junta por el hecho de la entrada en funcionamiento del compresor (∆Fpr )
6.56 × 109
∆Fpr = 13333 = 12385 N
6.56 × 109 + 5.02 × 108
Si tras entrar en servicio, la junta mantiene la misma fuerza que la que perdió por entrar en
servicio respecto a la pretensión inicial, entonces la fuerza que queda en la junta es Fj = ∆Fpr =
CAPÍTULO 11. PROBLEMAS RESUELTOS DE ELEMENTOS DE MÁQUINAS 221
12385 N. La fuerza de pretensión F , con la que se calculará el par de apriete en el tornillo será
F = 2 · ∆Fpr = 2 · 12385 = 24770 N. Por tanto, el par de apriete será
2 re3 − ri3
Ta = F Rm tan(α + φ) + µc 2
3 re − ri2
donde
µ 0.2
φ = arctan = arctan = arctan(0.207) = 11.7◦
cos β cos 15◦
18.9
re = = 9.45 mm
2
12
ri = = 6 mm
2
F = 24770 N
y por tanto
2 0.009453 − 0.0063
◦ ◦
Ta = 24770 0.01 · tan(2.5 + 11.7 ) + 0.2 · = 101.5 Nm
3 0.009452 − 0.0062
La fuerza Ft que soportará el tornillo será la carga de servicio P más la fuerza remanente en
la junta (que es la misma que ∆Fpr ), tal y como se muestra en la figura 11.8a)
5000 Fj 5000
5000
k't kpr
0 0
22.94 22.95 22.96 22.97 22.98 22.99 23.00 22.94 22.95 22.96 22.97 22.98 22.99 23.00
L[mm] L[mm]
a) b)
Figura 11.8: Fuerzas y desplazamientos en la unión atornillada: a) fuerzas para el apartado 1. b) fuerzas
para el apartado 2.
Solución al apartado 2
La nueva situación se representa en la figura 11.8b). En este nuevo caso P, sigue siendo la
misma, al igual que las rigideces de las placas y el tornillo, por lo que ∆Ft y ∆Fpr son los mismos
del apartado anterior. Si ahora hay que asegurar una fuerza en la junta de Fj = 5000 N. entonces
la nueva fuerza de pretensión F será
2 0.009453 − 0.0063
Ta0 = 17385 0.01 · tan(2.5◦ + 11.7◦ ) + 0.2 · = 71.2 Nm
3 0.009452 − 0.0062
Y la nueva fuerza en el tornillo será
[2] T. McGeer, “Passive dynamic walking,” International Journal of Robotics Research, pp. 62
– 82, 1990.
[3] C. Gosselin and J. Angeles, “Singularity analysis of closed-loop kinematic chains,” Robotics
and Automation, IEEE Transactions on, vol. 6, pp. 281 – 290, 07 1990.
[4] B. Karan and M. Vukobratović, “Calibration and accuracy of manipulation robot models,”
Mechanism and Machine Theory, vol. 29, pp. 479–500, 04 1994.
[7] Jorge Angeles, Damiano Pasini, Fundamentals of Geometry Construction. The Math Behind
the CAD. Springer International Publishing, 2020.
[10] A. M. Howatson, P. G. Lund, J. D. Todd , Engineering Tables and Data. Springer, 978-94-
010-9314-9.
[11] J. F. Schackelford, Introducción a la ciencia de los materiales para ingenieros, 6ª ed. Pearson
Prentice Hall, 2008.
[13] B.J. Hamrock, B. J. Jacobson, S.R. Schmid, Elementos de máquinas. McGraw-Hill, 1999.
[14] J. E. Shigley, L.D. Mitchell, Mechanical Engineering Design. McGraw-Hill, New York, 1983.
223
BIBLIOGRAFÍA 224
[16] W. C. Young, Roark’s Formulas for Stress and Strain. McGraw-Hill, New York, 1989.
[17] J. A. Collins, Failure of Materials in Mechanical Design. Wiley, New York, 1981.
[19] R.A. Smith, S. Hillman, “A brief historical overview of the fatigue of railway axles.,” Procee-
dings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit
218, vol. 4, pp. 267 – 277, 2004.
[20] C. Tassava, “Weak seams: Controversy over welding theory and practice in american ship-
yards, 1938 1946.,” History, and Technology, vol. 19, no. 2, pp. 87 – 108, 2003.
[21] R.C. Juvinall, K. M. Marshek, Fundamentals of Mchine component design. Wiley, New
York, 1991.
[23] C. Lipson, R.C. Juvinall, Handbook of Stress and Strength. Macmillan, New York, 1963.
[29] A. Palmgren, “The fatigue life of ball-bearings,” Z. VDI, vol. 68, pp. 339 – 341, 1924. [en
alemán].