Análisis Dinámico Por Métodos Numéricos

Capı́tulo 5
Análisis dinámico por

métodos numéricos
5.1. Introducción
Por medio del análisis dinámico tratamos de encontrar las ecuaciones di-
ferenciales que gobiernan el movimiento de un mecanismo como consecuen-
cia de su velocidad inicial o de las fuerzas exteriores. La Mecánica Clásica
proporciona diversos caminos para obtener dichas ecuaciones: aislando ca-
da elemento y aplicándoles las ecuaciones de Newton-Euler, mediante las
ecuaciones de Lagrange o empleando el Teorema de los Trabajos Virtuales.
Cualquiera que sea el camino seguido, el resultado es un conjunto de tantas
ecuaciones diferenciales ordinarias como grados de libertad.
El empleo de las coordenadas naturales no introduce ninguna variación
conceptual respecto de los procedimientos ya conocidos. Simplemente apor-
ta un soporte matemático adecuado para dar a los procedimientos un en-
foque general y sistemático. Gracias a las coordenadas naturales podremos
desarrollar las ecuaciones del movimiento de cualquier sistema mecánico,
sin importar su número de elementos, pares y grados de libertad. Como
veremos, para obtener las ecuaciones diferenciales del movimiento necesi-
taremos recurrir a las resoluciones de los problemas cinemáticos, vistos en
el capı́tulo anterior. Asimismo, necesitaremos presentar algunos conceptos
sobre la integración numérica de ecuaciones diferenciales.
La resolución del análisis dinámico nos proporciona la trayectoria de los
grados de libertad z(t). Para determinar las reacciones en los pares recurri-
remos al método de los multiplicadores de Lagrange, un método empleado
tradicionalmente en la minimización de funciones sujetas a restricciones,
151
152 Capı́tulo 5. Análisis dinámico por métodos numéricos
xe
ye P
e
rp j
y0 r
i
{0} ri
x0
Figura 5.1
que en este caso proporcionará la forma de asociar las fuerzas interiores a

las ecuaciones de restricción.
5.2. Fuerzas de inercia de un elemento

Partimos de un sólido definido mediante dos puntos i y j. Además
del sistema de referencia fijo {0}, definimos un sistema de referencia {e},
cuyo origen se encuentra en el punto i y cuyo eje xe atraviesa el punto
j. Denominamos Ae a la matriz de cambio de base entre los sistemas de
referencia {e} y {0}.
Tratemos de escribir las coordenadas de una partı́cula material de masa
dm situada en el punto P , en función de las coordenadas de los puntos i y
j. Para ello, llamando e rP a las coordenadas del punto P en el sistema de
referencia {e}, podemos escribir
� � � � �e �
e xi 1 (xj − xi ) −(yj − yi ) xP
r = ri + Ae r P = + ey (5.1)
yi Lij (yj − yi ) (xj − xi ) P
Mediante una simple reordenación de términos, esta expresión se puede

escribir en forma de producto de matriz por vector,
 

 xi 
1
�
Lij − e xP ey e x −e y �   
P P P yi
r= = TP qe (5.2)
Lij −e yP Lij − e xP e yP e xP   xj 

 
yj
donde TP es una matriz de tamaño 2 × 4, constante en el tiempo, pero

variable de una partı́cula a otra, y qe es el vector de coordenadas naturales
del elemento.
5.2. Fuerzas de inercia de un elemento 153
Para calcular la aceleración del punto P , basta con derivar dos veces la
ecuación anterior con respecto al tiempo, obteniendo
r̈ = TP q̈e (5.3)
También podemos calcular el desplazamiento virtual del punto P em-

pleando el operador virtual δ sobre la ecuación 5.2,
δr = TP δqe (5.4)
La fuerza de la inercia de la partı́cula P es sencillamente −r̈dm. El

trabajo virtual δWP producido en un desplazamiento virtual δr del punto P
se calculará mediante el producto escalar de la fuerza por el desplazamiento,
obteniendo
δWP = −δrT r̈ dm = −δqTe TTP TP q̈e dm (5.5)
Para calcular el trabajo virtual de todo el elemento, deberemos sumar

los trabajos virtuales de todas las partı́culas del volumen V del sólido,
mediante la integral
� �
e
δWiner = δWP = − δqTe TTP TP q̈e dm (5.6)
V V
En esta ecuación los términos δqe y q̈e no varı́an de una partı́cula a

otra, pues contienen las coordenadas de los puntos i y j, por lo que podemos
sacarlos fuera de la integral,
��
e T T
δWiner = −δqe TP TP dm q̈e (5.7)
V
El paréntesis representa una matriz simétrica de tamaño 4 × 4, que

multiplica al vector q̈e . El resultado de dicho producto es un vector de 4 × 1
que expresa las fuerzas de inercia del elemento mediante las coordenadas
naturales. Por el sentido fı́sico del paréntesis, que es el de una masa, recibe
el nombre de matriz de masas del elemento, que denominaremos Me, y
cuya expresión es
�
Me = TTP TP dm (5.8)
V
Con la ayuda de la ecuación 5.2, que contiene la expresión explı́cita de

la matriz TP , podemos expandir la ecuación anterior y obtener la expresión
c Alejo Avello. Tecnun – Universidad de Navarra

�
ye
1 G 2 xe
Figura 5.2
explı́cita de la matriz de masas del elemento. Operando y realizando algunas

simplificaciones que omitimos, se obtiene finalmente:
 
M + a − 2bx 0 bx − a −by
 M + a − 2bx by bx − a 
 
Me =   (5.9)
 a 0 
sim. a
con
Ii M e xG M e yG
a= b x = by = (5.10)
L2ij Lij Lij
�
M= dm (5.11)
V
�
Ii = (e x2P + e yP2 ) dm (5.12)
V
� �
e 1 e e 1 e
xG = xP dm yG = yP dm (5.13)
M V M V
siendo Lij la longitud entre los puntos i y j, M la masa del elemento, e xG ,

e y las coordenadas del centro de gravedad en el sistema de referencia {e},
G
e Ii el momento del inercia del elemento con respecto del punto i.
Finalmente, podemos escribir abreviadamente el trabajo virtual de las
fuerzas de inercia de un elemento como
e
δWiner = −δqTe Me q̈e (5.14)
Ejemplo 5.1 Calcular la matriz de masas Me del elemento de la figura

5.2, consistente en una barra de longitud L y masa M uniformemente dis-
tribuida, modelizada mediante dos puntos 1 y 2 situados en los extremos de
la barra.
En primer lugar, situamos el sistema de referencia {e} en el punto 1 y
5.3. Fuerzas de inercia de un mecanismo 155
dibujamos el eje xe apuntando hacia el punto 2. En este caso:
e L
xG =
2
e
yG = 0
M L2
I1 =
3
Sustituyendo estas ecuaciones en la expresión general de la matriz de
masas 5.9, se obtiene:
 
2 0 1 0
M 0 2 0 1
Me =
6 1 0 2 0
0 1 0 2
5.3. Fuerzas de inercia de un mecanismo

Un mecanismo está compuesto por múltiples elementos unidos entre
sı́ mediante pares cinemáticos. Para hallar el trabajo virtual de todas las
fuerzas de inercia de un mecanismo durante un desplazamiento virtual ar-
bitrario, podemos sumar le trabajo virtual de cada elemento, expresado en
la ecuación 5.14,
� �
e
δWiner = δWiner =− δqTe Me q̈e (5.15)
e e
Esta ecuación contiene un sumatorio extendido a cada elemento del

mecanismo en el que aparecen los vectores qe de coordenadas naturales de
cada elemento. Mediante una simple reordenación algebraica de los térmi-
nos, podemos transformar dicha ecuación, eliminando el sumatorio, para
reescribir el trabajo virtual de las fuerzas de inercia en función del vector
q de coordenadas naturales de todo el mecanismo,
δWiner = −δqT M q̈ (5.16)
El proceso de reordenación anterior se conoce comúnmente como en-

samblado de la matriz de masas, y consiste en generar la matriz de masas
del mecanismo a partir de las matrices de masas de los elementos que lo
componen. Una justificación completamente general de este proceso con-
lleva cierta complejidad, por lo que se ha optado por explicarlo mediante
ejemplos sencillos.

�
M,L 2
M,L
1
3
Figura 5.3
Comencemos con el mecanismo de la figura 5.3, compuesto por dos

barras articuladas de longitud L y masa M uniformemente distribuida. El
mecanismo se modeliza mediante los puntos 1, 2 y 3. Al elemento 1 le
corresponden los puntos 1 y 2, mientras que al elemento 2 le corresponden
los puntos 2 y 3.
El trabajo virtual de las fuerzas de inercia del mecanismo es la suma de
los trabajos virtuales de los dos elementos,
1 2
δWiner = δWiner + δWiner =
  
2 0 1 0  ẍ1 
� 0   
M� 2 0 1 ÿ1 +
= δx1 δy1 δx2 δy2 1
6 0 2 0 
 ẍ 
 2 
0 1 0 2 ÿ2
  
2 0 1 0  ẍ2 
� 0   
M� 2 0 1 ÿ2
+ δx2 δy2 δx3 δy3 1 (5.17)
6 0 2 0 
 ẍ 
 3 
0 1 0 2 ÿ3
Podemos reescribir la ecuación empleando el vector ri para referirnos

a las coordenadas xi , yi de los puntos 1 y 2, y empleando las submatrices
Meij para referirnos a la partición 2 × 2 de la matriz de masas del elemento
e correspondiente a la fila del punto i y la columna del punto j.
� ��
� T T � M111 M112 r̈1
δWiner = δr1 δr2 +
M121 M122 r̈2
� ��
� T T � M222 M223 r̈2
+ δr2 δr3 2 2 =
M32 M33 r̈3
= δrT1 (M111 r̈1 + M112 r̈2 )+
+ δrT2 (M121 r̈1 + M122 r̈2 + M222 r̈2 + M223 r̈3 )+
+ δrT3 (M232 r̈2 + M233 r̈3 ) (5.18)
2M
1 4 2 3
M M M
1 2 3
A B C
Figura 5.4
Finalmente, podemos reordenar todos los términos y escribir una ecua-

ción en la que aparezca q como único vector de coordenadas naturales,
 
1 1
 M r̈ + M r̈ 
� T T T� 1 11 1 12
1
2
2 2

δWiner = δr1 δr2 δr3 M21 r̈1 + M22 r̈2 + M22 r̈2 + M23 r̈3 =

M2 r̈ + M2 r̈ 

32 2 33 3
 1  
� T T T � M11 M112 0 r̈1 
= δr1 δr2 δr3 M121 M122 + M222 M223  r̈2 =
 
0 M232 M233 r̈3
 
2 0 1 0 0 0
 0 2 0 1 0 0  
 
� T T T � M 1 0 4 0 1 0 r̈1 
= δr1 δr2 δr3   r̈2 =
6 0 1 0 4 0 1 r̈3 

 0 0 1 0 2 0
0 0 0 1 0 2
= δqT Mq̈ (5.19)
Como se puede comprobar, la matriz de masas del mecanismo se ha

obtenido distribuyendo las submatrices de cada punto en la matriz global
del mecanismo. Cuando dos submatrices de elementos diferentes se super-
ponen, por corresponderles la misma posición en la matriz de masas global,
ambas submatrices se suman. El proceso puede generalizarse sin dificultad.
Para calcular la matriz de masas de un mecanismo, se obtienen en primer
lugar las matrices de masas de cada elemento, seguidamente se identifican
las submatrices de cada punto y, finalmente, se ensamblan las submatrices
en la matriz global, sumando aquéllas en las que se produzca superposición.
Ejemplo 5.2 Calcular la matriz de masas del mecanismo de la figura 5.4,

sabiendo que los cuatro elementos son barras de masa uniformemente dis-
tribuida.

�
Puesto que los cuatro elementos tienen masa distribuida y están modeli-
zados con puntos en sus extremos, se puede aplicar la expresión del ejemplo
5.1, obteniendo:
A 1 B 2
   
2 0 1 0 2 0 1 0
 A B
M 0 2 0 1

M 0 2 0 1
M1 = M2 =
6 1 0 2 0
1 6 1 0 2 0
2
0 1 0 2 0 1 0 2
C 3 1 3
   
2 0 1 0 2 0 1 0
C 1
M 0 2 0 1 M 0 2 0 1
M3 = M4 =
6 1 0 2 0
3 3 1 0 2 0
3
0 1 0 2 0 1 0 2
expresión en la que se indican el punto a que corresponden las filas y co-

lumnas.
El ensamblado de la matriz de masas requiere distribuir estas 16 sub-
matrices en la matriz de masas global. Omitiendo las submatrices que co-
rresponden a puntos fijos, resulta
 
M111 + M411 0 M413
M= 0 M222 0 
4 3 4
M31 0 M33 + M33
Y, operando, se obtiene finalmente
1 2 3
 
6 0 0 0 2 0
0 1
 6 0 0 0 2
M
0 0 2 0 0 0 2
M=
6
0 0 0 2 0 0
2 0 0 0 6 0
3
0 2 0 0 0 6
Ejemplo 5.3 Calcular la matriz de masas del mecanismo biela-manivela

de la figura 5.5, compuesto por tres elementos: la manivela (elemento 1)
tiene forma circular de radio R y masa M uniformemente distribuida; la
biela (elemento 2), es una barra de masa M uniformemente distribuida; la
deslizadera (elemento 3) tiene masa M .
1
2 M
R 2 M
A
3
1
2M
Figura 5.5
El vector de coordenadas naturales q está compuesto por las cuatro coor-

denadas de los puntos 1 y 2, es decir, x1 , y1 , x2 , y2 . Puesto que el meca-
nismo tiene un único grado de libertad, serán precisas tres ecuaciones de
restricción: a) distancia constante entre los puntos A y 1; b) distancia cons-
tante entre los puntos 1 y 2; c) coordenada y2 = 0.
Una opción alternativa consiste tomar un vector q que sólo incluya las
coordenadas x1 , y1 , x2 , aprovechando que y2 = 0, por lo que no constituye
una incógnita. En este caso, será suficiente con imponer las dos primeras
ecuaciones de restricción. Para obtener la matriz de masas adoptaremos es-
ta elección de coordenadas, por lo que la matriz de masas M tendrá tamaño
3 × 3. Las matrices de masas de los tres elementos son las siguientes:
A 1 1 2
   
3 0 −1 0 2 0 1 0
 0 3 A 1
0 −1  M 0 2 0 1
M1 = M 
 −1 0
 M2 =  
1 0  6 1 0 2 0
1 2
0 −1 0 1 0 1 0 2
� 2 �
1 0
M3 = M 2
0 1
La matriz M1 se obtiene particularizando para el elemento 1 la ex-

presión general de la matriz de masas de la ecuación 5.9, con masa 2M ,
longitud Lij = R, momento de inercia Ii = M R2 , y posición del centro de
gravedad e xG = e yG = 0. La matriz M3 es trivial, ya que el elemento 3
puede considerarse una masa puntual y el trabajo virtual de las fuerzas de
inercia se escribe directamente como
� ��
3
� � 1 0 ẍ2
δWiner = − δx2 δy2
0 1 ÿ2
Para ensamblar la matriz de masas del mecanismo, debemos distribuir

�
dr F
xe
ye P
e
rp j
y0 r
i
{0} ri
x0
Figura 5.6
las submatrices correspondiente a las filas y columnas del vector q. Operan-

do, se obtiene
x1 y1 x2
8 0 1 x1
M 
M= 0 8 0 y1
6 1 0 8 x2
5.4. Fuerzas exteriores puntuales

Consideramos una fuerza puntual F actuando sobre una partı́cula ma-
terial P del elemento e, como se muestra en la figura 5.6. El trabajo virtual
producido por F en un desplazamiento virtual arbitrario δr se obtiene me-
diante el producto escalar de los dos vectores,
δWF = δrT F (5.20)
Anteriormente, en la ecuación 5.4, hemos calculado el valor del des-

lazamiento δr en función de las coordenadas naturales del elemento qe .
Sustituyendo,
δWF = δqTe TTP F = δqTe QF (5.21)
donde el vector QF se denomina fuerza generalizada, y de acuerdo con la

ecuación 5.2, toma el valor
 
Lij − e xP −e yP � �
1 
ey
P Lij − e xP 
 Fx
QF = TTP F = (5.22)
Lij  e xP ey
P
 Fy
− e yP ex
P
5.5. Fuerzas introducidas por los actuadores 161
s
2
1
2
s 1
2 j
1
F F
2 2
1 1
t t
Figura 5.7 Figura 5.8
Un caso particular es la fuerza de la gravedad, en que la fuerza F se

encuentra aplicada en el centro de gravedad del elemento, y toma el valor
� �
0
F= (5.23)
−M g
También es frecuente el caso en que existan muelles o amortiguado-

res actuando en uno más de los elementos del mecanismo. Esta situación
se puede tratar incluyendo fuerzas iguales y de signo opuesto, cada una
actuando sobre uno de los extremos de los muelles o amortiguadores.
5.5. Fuerzas introducidas por los actuadores

La función de los actuadores es proporcionar fuerzas exteriores para
gobernar el movimiento del mecanismo. Pueden ser de dos tipos: lineales o
angulares. Los actuadores lineales introducen fuerzas F iguales y de signo
opuesto entre sus dos puntos de anclaje, de forma análoga a los muelles
(figura 5.7); a diferencia de éstos, que son elementos internos, los actuado-
res producen una fuerza generada por una fuente de energı́a externa, que
produce una ley de control pre-establecida. Los actuadores angulares fun-
cionan de forma análoga, pero en lugar de fuerzas introducen momentos
τ , iguales y de signos opuestos entre dos elementos del mecanismo (figura
5.8).
Para calcular el trabajo virtual producido por F o τ , es recomendable
definir coordenadas relativas lineales (s) o angulares (ϕ), según correspon-
da. En el caso del actuador lineal, el trabajo virtual se calcula como el
producto del módulo de la fuerza F por el desplazamiento virtual de la

�
2
2
1
1
3
M1g M2g
Figura 5.9
distancia s,
δWF = F δs (5.24)
Por su parte, el trabajo virtual del actuador angular se obtiene mediante

el producto del momento τ por el desplazamiento virtual del ángulo ϕ,
δWτ = τ δϕ (5.25)
5.6. Ensamblado del vector de fuerzas generaliza-

das
En un mecanismo pueden actuar numerosas fuerzas exteriores de los ti-
pos descritos en los apartados anteriores. El trabajo virtual de todas ellas se
calculará sumando los trabajos virtuales de cada una de las fuerzas. Análo-
gamente a lo visto en el apartado 5.3 para las fuerzas de inercia, la suma
de las fuerzas generalizadas es equivalente al ensamblado de los vectores
particulares de cada fuerza generalizada para constituir un único vector
Q. Llamando nF al número de fuerzas generalizadas exteriores, podemos
escribir la expresión
�
δWF ext = δWF = δqT Q (5.26)
nF
Para comprender mejor el proceso de ensamblado del vector Q, recurri-

mos al ejemplo de la figura 5.9, compuesto por dos barras rı́gidas de masa
M1 y M2 uniformemente distribuida, y articuladas en el punto 2. Pues-
to que los pesos actúan en los centros de gravedad de cada elemento, que
5.7. Ecuaciones del movimiento en coordenadas independientes 163
coinciden con los centros geométricos de las barras, se puede escribir

   
1/2 0 � � 
 0 

    1
T  0 1/2  0 −M 1 g/2
Q1 = TP F =  = (5.27)
1/2 0  −M1 g 
 0 

  2
0 1/2 −M1 g/2
   
1/2 0 � � 
 0 

    2
0 1/2 0 −M 2 g/2
Q2 = TP F = 
T 
 1/2 0  −M2 g = (5.28)

 0 

  3
0 1/2 −M2 g/2
Ensamblando los dos vectores Q1 y Q2 en un único vector Q de fuerzas

generalizadas del mecanismo, se obtiene:
 

 0 


 
 1

 −M1 g/2 

 
0
Q= 2 (5.29)

 −(M 1 + M2 )g/2 


 


 0 

  3
−M2 g/2
Este proceso puede generalizarse fácilmente para cualquier tipo de fuer-

za exterior. El concepto es siempre el mismo: primero se forma el vector
particular de cada fuerza generalizada, y después se ensamblan en el vector
global, encajándolos en las posiciones correspondientes a las coordenadas
afectadas. Cuando dos o más términos coincidan en una misma posición,
se suman.
5.7. Ecuaciones del movimiento en coordenadas

independientes
Del teorema de los trabajos virtuales sabemos que
δWiner + δWext = 0 (5.30)
En los apartados anteriores hemos calculado las expresiones de δWiner

y δWext en función del vector q. Sustituyendo las expresiones 5.16 y 5.26,
se obtiene
−δqT Mq̈ + δqT Q = 0 (5.31)

�
Cambiando de signo y sacando factor común,
δqT (Mq̈ − Q) = 0 (5.32)
No podemos concluir que el paréntesis sea nulo, debido a que las compo-
nentes de q son variables dependientes, relacionados entre sı́ por las ecuacio-
nes de restricción. Para justificar esta afirmación, consideremos un ejemplo
de dos dimensiones en que las ecuaciones anteriores se particularizan arbi-
trariamente con los siguientes valores:
��
� � 2 1 ẍ 1
δx δy − =0 (5.33)
1 2 ÿ 1
y supongamos que las dos variables x, y están relacionadas mediante la
ecuación
y = 2x (5.34)
Precisamente porque esta relación existe, la ecuación 5.33 puede anu-

larse sin necesidad de que el paréntesis sea nulo. Concretamente, cuando
el paréntesis toma el valor {−2 1}T o cualquier múltiplo de él, se verifi-
cará que
� �
� � −2
δx δy = −2δx + δy = 0 (5.35)
1
ecuación que siempre es nula en virtud de la relación 5.34.
A la vista de que no podemos eliminar los desplazamientos virtuales por
no ser independientes, trataremos de reducir el número de variables hasta
que su número coincida con el número de grados de libertad del mecanismo.
Llamaremos z al vector de coordenadas independientes, de dimensión g ×
1, obtenido extrayendo tantas componentes del vector q como grados de
libertad.
Conocido el vector z podemos obtener el vector q resolviendo el pro-
blema de posición. Por tanto, podemos pensar en q como una función del
vector z, escribiendo
q = f(z) (5.36)
Derivando respecto al tiempo,

 ∂f ∂f1

∂z1
1
· · · ∂z g
∂f  . .. .. 
q̇ = ż = 
 .
. . 
.  ż = R ż (5.37)
∂z ∂fn
∂z1 · · · ∂f
∂zg
n
donde R es una matriz de tamaño n × g, y representa la matriz jacobiana

de la función vectorial f respecto de z. La ecuación anterior demuestra que,
para una posición dada, la relación existente entre q̇ y ż es lineal.
Sabemos que los desplazamientos virtuales satisfacen las mismas ecua-
ciones que las velocidades. Entonces, podemos escribir una ecuación análoga
a la 5.37,
δq = R δz (5.38)
Sustituyendo en 5.32, se obtiene
δzT RT (Mq̈ − Q) = 0 (5.39)
Pero en esta ecuación los desplazamientos virtuales δz son arbitrarios e

independientes, por lo que debe cumplirse
RT (Mq̈ − Q) = 0 (5.40)
o bien,
RT Mq̈ = RT Q (5.41)
Ahora debemos sustituir el vector q̈ por z̈, y para ello derivamos la

ecuación 5.37,
d(Rż)
q̈ = = Rz̈ + Ṙż (5.42)
dt
Sustituyendo este valor en la ecuación 5.41, se obtiene finalmente
RT MRz̈ = RT Q − RT MṘż (5.43)
Esta expresión representa un sistema de ecuaciones diferenciales ordi-

narias de segundo orden, que rigen el movimiento del mecanismo. Las di-
mensiones de los distintos vectores y matrices se pueden apreciar de forma
más intuitiva en la figura 5.10.
.. ..
g RT z= RT Q- RT Rz
n M R M
Figura 5.10

�
Ya hemos visto anteriormente cómo calcular la matriz de masas M y

el vector de fuerzas generalizadas Q. Ahora veremos cómo calcular en la
práctica la matriz R y el vector Ṙż.
Para calcular la matriz R recurrimos a la ecuación 5.37, que reprodu-
cimos a continuación,
q̇ = Rż (5.44)
De entre todos los vectores ż posibles, escojamos arbitrariamente aquél

que tiene todas sus componentes nulas, excepto la primera, que toma valor
unidad, es decir, ż = {1 0 · · · 0}T . Con esta elección particular, podemos
realizar el producto Rż, indicando cada una de las g columnas de la matriz
R mediante un superı́ndice,
 

 1 

0 
� 1 � 1
q̇ = R R2 · · · Rg .. = R (5.45)

 . 


  
0
Podemos interpretar este resultado diciendo que R1 , primera columna

de la matriz R, es el vector q̇ correspondiente a las velocidades de los grados
de libertad ż = {1 0 · · · 0}T . El proceso de calcular el vector q̇ para un vector
ż dado constituye el problema cinemático de velocidad, cuya resolución ya
se ha abordado.
El razonamiento anterior es fácilmente generalizable: para calcular Ri ,
columna i de la matriz R, deberemos resolver el problema de velocidad con
un vector ż tal que żi = 1 y żj = 0, ∀j �= i. Repitiendo g veces este proceso,
se obtendrán todas las columnas de R.
Para calcular el producto Ṙż recurrimos a la ecuación 5.42, que esta-
blece
q̈ = Rz̈ + Ṙż (5.46)
En el caso particular en que z̈ = 0, se observa que el producto Ṙż

coincide exactamente con el vector de aceleraciones q̈. Por tanto, el pro-
ducto Ṙż representa fı́sicamente la aceleración del mecanismo debida a la
velocidad, que es suma de los términos centrı́fugos y de Coriolis.
En términos prácticos, para calcular Ṙż necesitamos hallar la acelera-
ción q̈ correspondiente a z̈ = 0. Como se recordará, al proceso de calcular
q̈ conocido el vector z̈ lo denominamos problema cinemático de aceleración,
y su resolución ha sido ya estudiada. Por tanto, el problema de calcular el
vector Ṙż se reduce a resolver un problema de aceleración.
1
L= 2 L= 2
.
j
45º 2
A
Figura 5.11
Ejemplo 5.4 Calcular la aceleración instantánea del mecanismo biela-manivela

de la figura 5.11, suponiendo que se encuentra en un plano horizontal
y que se deja evolucionar libremente con velocidad angular instantánea
ϕ̇ = 2 rad/s. Las dos barras tienen masa m uniformemente distribuida.
Empleando como coordenadas naturales el vector qT = {x1 y1 x2 ϕ},
las ecuaciones de restricción y su correspondiente matriz jacobiana (parti-
cularizada para la posición instantánea), serán:
 
 (x1 − xA )2 + (y1 − yA )2 − 2 
Φ = (x2 −√x1 )2 + (y2 − y1 )2 − 2
 
x1 − 2 cos ϕ
Y la matriz jacobiana,
 
2(x1 − xA ) 2(y1 − yA ) 0 0
Φq =  −2(x2 − x1 ) −2(y2 − y1 ) 2(x2 − x1 ) √ 0 
1 0 0 2 sen ϕ
que particularizada para la posición indicada, resulta

 
2 2 0 0
Φq =  −2 2 2 0
1 0 0 1
Para calcular la matriz R debemos realizar sucesivos análisis de veloci-

dades dando valor unitario a cada grado de libertad y cero a los demás. Co-
mo en este caso tenemos un solo grado de libertad, la matriz R tendrá una
única columna, que obtendremos mediante un análisis de velocidad con va-
lor ϕ̇ = 1 rad/s. Conviene resaltar que éste no es el valor real de la velocidad
(el enunciado dice que la velocidad es ϕ̇ = 2 rad/s) sino un valor ficticio
que empleamos con la única finalidad de calcular la matriz R. Para ello,
emplearemos la ecuación
Φq q̇∗ = −Φt = 0

�
en la que empleamos el superı́ndice “*” para indicar que se trata de velo-

cidades ficticias. Al sistema de ecuaciones anterior le que añadiremos la
ecuación ϕ̇∗ = 1 como una fila adicional para introducir el valor conocido
de la velocidad. Con ello,
  ∗   
2 2 0 0   ẋ1 
 0 
 −2 2 2 0   ẏ1∗   0 
 
 1 0 0 1   ẋ∗2  =  0 

 ∗    
0 0 0 1 ϕ̇ 1
La solución de este sistema nos proporciona el valor de la matriz R:

 ∗  

 ẋ   −1 
 ∗1      
ẏ1 1
R= =

 ẋ∗   −2 
 2∗      
ϕ̇ 1
A continuación podemos calcular las velocidades reales del mecanismo

mediante la expresión q̇ = Rż, resultando
     

 ẋ1 
 
 −1   
 −2 
      
ẏ1 1 2
q̇ = = 2 =

 ẋ   −2   −4 
 2      
  
ϕ̇ 1 2
� ��
R ż
Seguidamente, calculamos el término Ṙż, que se obtiene mediante un

análisis de aceleraciones con z̈ = 0, es decir con ϕ̈ = 0. Resaltamos de
nuevo que dicha aceleración es ficticia, ya que el valor de ϕ̈ es descono-
cido (precisamente es lo que se pide en el enunciado). Para el análisis de
aceleraciones empleamos la ecuación:
Φq q̈∗ = −Φ̇q q̇
a la que añadiremos la ecuación ϕ̈∗ = 0 como fila adicional, para introducir

el valor conocido de la aceleración. Desarrollando las ecuaciones, obtene-
mos
  ∗   2   
2 2 0 0   ẍ1 
 
 2ẋ1 + 2ẏ12 
 
 16 
 −2   ∗  2 2    
 2 2 0  ÿ1 2(
√ ẋ2 − ẋ1 ) + 2(ẏ2 − ẏ1 ) 16
 1 = − = −
0 0 1  ẍ∗   2 cos ϕ ϕ̇2   4 
 2∗  
 
 
  
0 0 0 1 ϕ̈ 0 0
Y resolviendo estas ecuaciones, obtenemos

 ∗  

 ẍ1  
 −4 

 ∗  
ÿ1 −4
Ṙż = =
 ẍ∗2 
   −8 
 ∗   

ϕ̈ 0
La matriz de masas M se obtiene ensamblando las matrices de masas

de los dos elementos, como se ha mostrado en el apartado 5.3. Según se ha
visto, las matrices de cada elemento son:
A 1 1 2
   
2 0 1 0 2 0 1 0
A 1
M
 0 2 0 1 M0 2 0 1
M1 =  M2 =
6 1 0 2 0 6 1 0 2 0
1 2
0 1 0 2 0 1 0 2
Y la matriz de masas del mecanismo resulta:
x 1 y1 x 2 ϕ
 
4 0 1 0 x1
M
 0 4 0 0 
 y1
M=
6 1 0 2 0  x2
0 0 0 0 ϕ
El vector Q de fuerzas generalizadas es nulo, ya que no actúan fuer-

zas exteriores (el peso no actúa porque el mecanismo se halla en un plano
horizontal).
Sustituyendo los términos anteriores en la ecuación diferencial del mo-
vimiento 5.43 particularizada para el instante actual, tenemos:
  
4 0 1 0   −1 

   10
� � M 0 4 0 0 
 1
RT MR = −1 1 −2 1 = M
6  1 0 2 0  −2 
 3
 
0 0 0 0 1
  
4 0 1 0   −4 

  
� � M 0 4 0 0 
 −4
RT MṘż = −1 1 −2 1 = −8M
6  1 0 2 0  −8 

 
0 0 0 0 0
RT Q = 0

�
y, finalmente, se obtiene:
10
M ϕ̈ = −8M
3
de donde se deduce la solución
12
ϕ̈ = −
5
5.8. Integración numérica de las ecuaciones del

movimiento
Hemos visto que las ecuaciones del movimiento se describen matemáti-
camente mediante un sistema de ecuaciones diferenciales ordinarias, no li-
neales y de segundo orden. En este apartado revisaremos brevemente al-
gunos conceptos básicos para resolver numéricamente las ecuaciones dife-
renciales. Puesto que las fórmulas habituales de integración numérica están
pensadas para sistemas de primer orden de la forma ẏ = f(y, t), será preciso
reducir previamente el orden de nuestro sistema de ecuaciones, pasando de
segundo a primer orden.
5.8.1. Reducción del orden de las ecuaciones diferenciales

del movimiento
Todo sistema de ecuaciones diferenciales de segundo orden puede ser
transformado en otro de primer orden duplicando el número de ecuaciones
y de incógnitas. En nuestro caso, necesitamos reducir el orden del sistema
de g ecuaciones diferenciales dado por la expresión 5.43, que reproducimos
a continuación,
RT MRz̈ = RT Q − RT MṘż (5.47)
A estas g ecuaciones, podemos añadirles otras g nuevas ecuaciones tri-

viales para obtener el siguiente sistema ampliado,
� T �� T �
R MR 0 z̈ R Q − RT MṘż
= (5.48)
0 I ż ż
Ahora, definimos el vector y, de tamaño 2g × 1

� � � �
y1 ż
y= = (5.49)
y2 z
5.8. Integración numérica de las ecuaciones del movimiento 171
en el que y1 equivale a ż, e y2 equivale z. Teniendo en cuenta que las

matrices R y M dependen exclusivamente de z, y que los vectores Q y
Ṙż dependen de z y ż, podemos concluir que cada término de la ecuación
5.49 depende de y, puesto que dicho vector ya incluye a z y ż. Con estas
consideraciones, y despejando de la ecuación anterior, se obtiene,
� � � �−1 � �
ẏ1 RT MR 0 RT Q − RT MṘż
ẏ = = = f(y, t) (5.50)
ẏ2 0 I y1
que representa un sistema de 2g ecuaciones diferenciales de primer orden.

A este sistema podemos aplicarle cualquiera de los métodos de integración
numérica descritos en la bibliografı́a, e incluidos en programas de cálculo
matemático como Maple o Matlab.
5.8.2. Métodos de integración explı́citos e implı́citos

Los métodos de integración numérica tratan de hallar una solución dis-
creta aproximada del sistema de ecuaciones diferenciales, cuya solución
exacta es continua. Dichos métodos se clasifican en explı́citos e implı́ci-
tos. Para entender la diferencia entre ambos podemos recurrir al método
de integración de Euler, el más sencillo de todos ellos. Consideremos la
ecuación de una variable
ẏ = f (y, t) (5.51)
Llamamos yn al valor de la función en el instante tn , es decir, yn = y(tn ),

y denominamos h a la diferencia tn+1 −tn . En el método de Euler explı́cito, a
partir el valor de yn aproximamos la solución en yn+1 empleando la derivada
en tn ,
yn+1 = yn + f (yn , tn ) h (5.52)
cuya interpretación geométrica se muestra en la figura 5.12.

El método de Euler implı́cito aproxima la solución en yn+1 utilizando
la derivada en tn+1 ,
yn+1 = yn + f (yn+1 , tn+1 ) h (5.53)
con la interpretación geométrica de la figura 5.13. El método de Euler

implı́cito requiere resolver esta última ecuación iterativamente, ya que el
valor de yn+1 aparece a ambos lados de la ecuación sin que, en general, sea
posible despejarla. Normalmente, los métodos explı́citos son más fáciles de

�
y yn+1 y
yn+1
yn h yn
h
t t
tn tn+1 tn tn+1
programar pero adolecen de falta de estabilidad. Por ello, es frecuente tener

que recurrir a los métodos implı́citos.
Matlab dispone de diversas funciones para resolver sistemas de ecuacio-
nes diferenciales ordinarias: ode23, ode45, ode113, ode15s, ode23s, ode23t,
ode23tb. Cada una de ellas implementa distintas fórmulas de integración,
y están más indicadas para ciertos tipos de problemas. Las funciones ode45
y ode113 suelen ser una buena elección para muchos problemas mecánicos
de tamaños pequeño y mediano.
5.9. Cambio de coordenadas independientes

En el apartado 5.7 hemos definido el vector z como una extracción de g
componentes de q. Si se ha realizado una buena elección, las componentes
de z constituyen los grados de libertad del mecanismo, por lo que definen
su posición. Sin embargo, a medida que el movimiento del mecanismo hace
cambiar su posición, los grados de libertad pueden volverse inservibles.
Por ejemplo, en la figura 5.14 la elección de y1 como grado de libertad
es adecuada inicialmente, pero cuando el cuadrilátero alcanza la posición
de la figura 5.15, la coordenada y1 no sirve como grado de libertad. Al
encontrarse la barra en posición vertical ẏ1 debe ser siempre nulo, por
lo que no proporciona ninguna información sobre el movimiento real del
sistema.
Cuando las coordenadas independientes se vuelven inservibles es preciso
cambiarlas. En el ejemplo anterior, cuando la barra se encuentra vertical
podemos sustituir y1 por la coordenada x1 . En el caso más general, la
selección de un conjunto de grados de libertad adecuado no siempre es
tarea sencilla, especialmente en los mecanismos tridimensionales.
El cambio de coordenadas independientes se puede evitar casi siempre
escogiendo grados de libertad válidos para todo el rango del movimiento,
5.10. Ecuaciones del movimiento en coordenadas dependientes 173
y1 x1
y1
j
j q j q
lo que habitualmente se consigue mediante coordenadas relativas de ángulo

o de distancia.
Conviene señalar que el mero hecho de escoger coordenadas relativas
como grados de libertad no garantiza su validez en todo el rango de movi-
miento. En ciertos problemas, la necesidad de cambiar los grados de libertad
es inherente el propio problema, por lo que es inevitable. Por ejemplo, en
un cuadrilátero articulado con movimiento de doble balancı́n, los topes del
movimiento obligan al cambio de coordenadas independientes. En la figura
5.16 el ángulo ϕ es válido como grado de libertad, pero al alcanzar la posi-
ción de la figura 5.17 se vuelve inservible. Lo contrario le ocurre al ángulo
θ.
5.10. Ecuaciones del movimiento en coordenadas

dependientes
Según hemos visto, la aplicación del teorema de los trabajos virtuales
con coordenadas naturales conduce a la expresión
δqT (Mq̈ − Q) = 0 (5.54)
Sabemos que de esta ecuación no podemos eliminar los desplazamien-

tos virtuales δq, por no ser independientes. Una primera alternativa, recién

�
vista, consiste en convertir las coordenadas en independientes. Otra al-

ternativa –que veremos en este apartado– consiste en emplear la técnica
matemática de los multiplicadores de Lagrange, estudiada en Cálculo para
la minimización de funciones sujetas a restricciones.
Partimos de las m ecuaciones de restricción que representábamos me-
diante la expresión
Φ(q) = 0 (5.55)
Consideramos un desplazamiento virtual arbitrario δq, compatible con

los enlaces, que provoca una variación Φ en las ecuaciones de restricción,
∂Φ
δΦ = δq = Φq δq = 0 (5.56)
∂q
Añadamos ahora a la ecuación 5.54 un término adicional de valor nulo,

resultado de multiplicar escalarmente el vector δΦ por un vector λ arbitra-
rio, de tamaño m × 1, denominado vector de multiplicadores de Lagrange,
δqT (Mq̈ − Q) + δΦT λ = 0 (5.57)
Sustituyendo en ella la ecuación 5.56, se obtiene
δqT (Mq̈ − Q + ΦTq λ) = 0 (5.58)
En esta ecuación, el vector δq sigue sin poder ser eliminada por el mis-
mo motivo que antes, pues sus componentes son dependientes. Pero algo
ha cambiado: ahora tenemos un vector λ arbitrario cuyas componentes
podemos escoger según convenga. Pues bien, de los infinitos vectores posi-
bles escogeremos, precisamente, aquél que obligue al paréntesis a ser nulo,
cualquiera que sea el vector δq. Con ello, conseguimos eliminar el vector
δq aunque sus componentes sean dependientes, a costa de optar por una
elección concreta del vector λ.
Expliquemos este proceso mediante el ejemplo de dos dimensiones visto
anteriormente. Consideramos el sistema descrito por la ecuación
��
� � 2 1 ẍ 1
δx δy − =0 (5.59)
1 2 ÿ 1
sujeto a la restricción
φ = y − 2x = 0 (5.60)
5.10. Ecuaciones del movimiento en coordenadas dependientes 175
La matriz jacobiana es,

� �
∂φ ∂φ � �
φq = = −2 1 (5.61)
∂x ∂y
Añadiendo el producto δΦT λ a la ecuación 5.59, obtenemos

��
� � 2 1 ẍ 1 −2
δx δy − + λ =0 (5.62)
1 2 ÿ 1 1
Operando,
δx(2ẍ + ÿ − 1 − 2λ) + δy(ẍ + 2ÿ − 1 + λ) = 0 (5.63)
Debido a que δx y δy están relacionados, no podemos afirmar que los dos

paréntesis sean nulos para cualquier desplazamiento virtual. Ahora bien,
podemos escoger un valor de λ arbitrario y, por conveniencia, tomamos
el valor que anula el segundo paréntesis, es decir, λ=1 − ẍ − 2ÿ. Con esta
elección, el primer paréntesis tiene que ser forzosamente nulo, ya que el des-
plazamiento virtual δx es arbitrario. Es decir, al introducir el multiplicador
de Lagrange podemos eliminar los desplazamientos virtuales dependientes,
siempre y cuando la elección de λ sea la apropiada.
Consecuentemente, podemos eliminar el término δq de la ecuación 5.58,
obteniendo
Mq̈ − Q + ΦTq λ = 0 (5.64)
que representa un sistema de n ecuaciones con n + m incógnitas, que son

las coordenadas dependientes q y los multiplicadores de Lagrange λ. Pa-
ra complementar este sistema, podemos recurrir al análisis cinemático de
aceleración que, como se vio, es
Φq q̈ = −Φ̇q q̇ − Φ̇t (5.65)
El conjunto de las dos ecuaciones 5.64 y 5.65 constituyen el siguiente

sistema de n + m ecuaciones
� ��
M ΦTq q̈ Q
= (5.66)
Φq 0 λ −Φ̇q q̇ − Φ̇t

�
Ry
Rx
A A
L L 1
45º 45º A
M M
45º
Mj L
..
Mg 2
5.11. Interpretación fı́sica de los multiplicadores

de Lagrange
Cada multiplicador de Lagrange se encuentra asociado a una –y sólo
una– ecuación de restricción. Se puede demostrar con relativa facilidad que
el producto ΦTq λ de la ecuación 5.64 tiene dimensiones de fuerza generaliza-
da y que fı́sicamente representa la fuerza necesaria para que las restricciones
se cumplan. Quizá la forma más sencilla de abordar una demostración no
rigurosa sea mediante un ejemplo.
Consideremos la barra rı́gida mostrada en la figura 5.18, de longitud L
y con una masa puntual M en el extremo. La barra se encuentra en reposo,
inclinada 45o en un plano vertical, y se desea calcular la aceleración inicial
del extremo y las reacciones en la articulación. Las ecuaciones de la mecáni-
ca clásica fácilmente nos permiten calcular estos términos. Basándonos en
la figura y tomando momentos en el punto A, obtenemos
√
g 2
ϕ̈ = − (5.67)
2L √
2 g
ẍ = ÿ = ϕ̈L =− (5.68)
2 2
Realizando los sumatorios de fuerzas horizontales y verticales, se obtiene
mg
Rx = − (5.69)
2
mg
Ry = (5.70)
2
Resolvamos ahora el mismo problema empleando coordenadas natura-
les. Para ello, empleamos la modelización mostrada en la figura 5.19, en la
que se utilizan los puntos 1 y 2. Por tanto, nuestro vector de coordenadas
5.11. Interpretación fı́sica de los multiplicadores de Lagrange 177
dependientes es
� �
qT = x 1 y1 x 2 y2 (5.71)
Puesto que el mecanismo tiene un solo grado de libertad, deben existir

n − m = 3 ecuaciones de restricción, que son
 
 (x2 − x1 )2 + (y2 − y1 )2 − L2 
Φ = x1 =0 (5.72)
 
y1
La correspondiente matriz jacobiana es

 
−2(x2 − x1 ) −2(y2 − y1 ) 2(x2 − x1 ) 2(y2 − y1 )
Φq =  1 0 0 0 = (5.73)
0 1 0 0
 √ √ √ √ 
− 2L 2L 2L − 2L
= 1 0 0 0 
0 1 0 0
Y la matriz de masas del elemento,

 
0 0 0 0
 0 0 0 0 
M=
 0
 (5.74)
0 M 0 
0 0 0 M
Finalmente, podemos escribir el sistema completo ensamblado,

 √    
0 0 0 0 −√ 2L 1 0  ẍ1   0 
 0    
 0 0 0 √2L 0 1
 ÿ 

  1 




  0 



   
 0 0 M 0 √2L 0 
0   ẍ2   0 
  
 
 √0 √0 √0 M
√ − 2L 0 0  ÿ2 = −mg (5.75)
  
 
 

 − 2L 2L 2L − 2L 0 0 0  λ1   0 
   
 

 1 0 0 0 0 0 0  λ
 
 2

 

 
 0 



0 1 0 0 0 0 0 λ 3 0
Pese al tamaño del sistema, el gran número de términos nulos hace que

�
la resolución a mano sea una tarea sencilla, de la que se obtiene la solución:
ẍ1 = ÿ1 = 0 (5.76)

g
ẍ2 = ÿ2 = − (5.77)
2
mg
λ1 = √ (5.78)
2 2L
mg
λ2 = (5.79)
2
mg
λ3 = − (5.80)
2
La fuerza generalizada correspondiente a las ecuaciones de restricción
segunda y tercera resulta ser
   mg 

 1 − 2 

  mg   

T 0 0
Q2 = −Φq2 λ2 = − = (5.81)

 0 2 
 0 
   
0 0
   

 0  0 
  � mg �  mg 
  
T 1 2
Q3 = −Φq3 λ3 = − − = (5.82)

 0 2  0 
   
 

0 0
resultados que coinciden con los obtenidos por el método clásico para Rx y
Ry en las ecuaciones 5.69 y 5.70.
5.12. Cálculo de las reacciones

Cada ecuación de restricción lleva asociado su multiplicador de Lagran-
ge y su correspondiente interpretación fı́sica. En este apartado nos pre-
guntamos por el sentido fı́sico de cada una de las restricciones, lo que nos
permitirá hallar las reacciones en cada par cinemático.
5.12.1. Reacciones de sólido rı́gido

Los multiplicadores de Lagrange asociados a las restricciones de sólido
rı́gido no tienen interés práctico inmediato. Representan las fuerzas internas
entre partı́culas para mantener constantes sus posiciones relativas a lo largo
del movimiento. No representan reacciones en los pares.
5.12. Cálculo de las reacciones 179
1 l1 l1 2
1 2
l2 l2
5.12.2. Reacciones en la articulación

Para hallar las reacciones de una articulación se precisa no compartir los
puntos entre elementos adyacentes, como en la figura 5.20, de forma que la
ecuación de restricción del par aparezca explı́citamente. Las dos ecuaciones
del par son
Φ 1 = x1 − x 2 = 0 (5.83)
Φ2 = y1 − y2 = 0 (5.84)
Las filas y columnas de la matriz jacobiana asociadas a estas dos res-
tricciones son
� �
Φq1 = 1 0 −1 0 (5.85)
� �
Φq2 = 0 1 0 −1 (5.86)
Llamando λ1 y λ2 a los correspondientes multiplicadores de Lagrange
obtenemos las fuerzas generalizadas Q1 y Q2
   

 1  
 −λ1  
   
T 0 0
Q1 = −Φq1 λ1 = − λ1 = (5.87)

 −1  
 λ1  
   
0 0
   

 0 
 
 0 

   
T 1 −λ 2
Q2 = −Φq2 λ2 = − λ2 = (5.88)

 0  
 0  
   
−1 λ2
El sentido fı́sico de dicho resultado se muestra en la figura 5.21. Las dos
primeras filas de cada vector representan las componentes de la reacción
sobre el punto 1, y las dos segundas filas representan las componentes de
la reacción sobre el punto 2. Por el principio de acción y reacción, ambas
reacciones deben ser iguales y de signo contrario.
De las dos ecuaciones anteriores comprobamos que los dos multiplica-
dores de Lagrange λ1 y λ2 tienen un claro sentido fı́sico: λ1 representa el
módulo de la reacción Rx , y λ2 el módulo de la reacción Ry .

�
L23
1 lL23
3
90º 1
3 3
1 L12 90º
2
2 L13 2 lL23
Figura 5.22 Figura 5.23 Figura 5.24
5.12.3. Reacciones en la deslizadera articulada

Como se vio, una deslizadera articulada se modeliza mediante tres pun-
tos, como se muestra en la figura 5.22. Su ecuación de restricción debe
garantizar el alineamiento entre los puntos 1, 2 y 3, lo que se consigue
mediante el producto vectorial entre los vectores r12 y r13 , cuya expresión
es
Φ = (x2 − x1 )(y3 − y1 ) − (y2 − y1 )(x3 − x1 ) = 0 (5.89)
La fila correspondiente en la matriz jacobiana es,
� 1 2 3 � (5.90)
Φq = −(y3 − y2 ) (x3 − x2 ) (y3 − y1 ) −(x3 − x1 ) −(y2 − y1 ) (x2 − x1 )
Y la fuerza generalizada,
 

 (y3 − y2 )   1

 −(x3 − x2 )  

 

 
T −(y 3 − y 1 )
Q = −Φq λ = λ 2 (5.91)

 (x3 − x1 )  

 

 (y2 − y1 )  

  3
−(x2 − x1 )
Observemos que las reacciones aparecen como fuerzas generalizadas dis-

tribuidas en los puntos 1, 2 y 3, como se muestra en la figura 5.23. Fácil-
mente se puede comprobar que la suma de las fuerzas de los puntos 2 y
3, coincide con la fuerza del punto del 1 cambiada de signo, como postula
el principio de acción y reacción. La acción de estas dos fuerzas equivale a
una sola fuerza aplicada en la deslizadera, como se ilustra en la figura 5.24.
5.12. Cálculo de las reacciones 181
1 lL34 1
2 4 N
1 4 lL34 2 4
a
2 lL12
a lL12
3 N
3 a 3
Figura 5.25 Figura 5.26 Figura 5.27
5.12.4. Reacciones en la deslizadera rı́gida

La deslizadera rı́gida se modeliza mediante cuatro puntos, tres de ellos
alineados a lo largo del eje, como en la figura 5.25. Es preciso imponer dos
ecuaciones de restricción: la primera, de producto vectorial nulo entre los
vectores r23 y r24 obliga a que los puntos 2, 3 y 4 permanezcan alineados;
la segunda ecuación, de producto escalar constante entre los vectores r12 y
r34 , obliga a que los segmentos 12 y 34 formen ángulo α constante. Omiti-
remos la primera de las ecuaciones, por ser análoga a la recién vista para
la deslizadera articulada. La segunda ecuación establece
Φ = (x2 − x1 )(x4 − x3 ) + (y2 − y1 )(y4 − y3 ) − c = 0 (5.92)
Su fila en la matriz jacobiana es,
� 1 2 3 4 � (5.93)
Φ q = x 3 − x 4 y3 − y4 x 4 − x 3 y4 − y 3 x 1 − x 2 y1 − y2 x 2 − x 1 y 2 − y1
Y la fuerza generalizada,
 
 x3 − x4 

  1

 y3 − y4  


 


 x 4 − x 3 


 
 2
T y 4 − y 3
Q = −Φq λ = −λ (5.94)

 x1 − x2  

 
 3

 y 1 − y 2 


 


 x 2 − x 1 

  4
y2 − y1
Podemos ver la representación de la fuerzas generalizadas en la figura
5.26. Fácilmente comprobamos que la resultante neta en cada elemento es
un momento, que se muestra en la figura 5.27. Además, comprobamos que
los momentos en cada elemento tienen el mismo módulo
N = λL12 L34 sen α (5.95)

�
pero signos contrarios, como establece el principio de acción y reacción.

Por tanto, y como cabı́a esperar, la fuerza generalizada Q asociada a la
restricción de ángulo constante es un momento.
5.12.5. Reacciones de coordenada relativa

Los multiplicadores de Lagrange asociados a las restricciones de coorde-
nada relativa no tienen significado fı́sico, por lo que tampoco tienen interés
práctico.
5.12.6. Esfuerzos motores

Los esfuerzos motores son los que se aplican en los grados de libertad
para producir el movimiento deseado. En el problema dinámico directo, a
partir del valor de los esfuerzos motores τ (t) calculamos la trayectoria del
mecanismo z(t). En el problema dinámico inverso hacemos exactamente
lo contrario: a partir de la trayectoria del mecanismo z(t) calculamos los
esfuerzos motores.
Existen varias formas de resolver este problema. Una de ellas consiste
en contemplar el movimiento conocido del mecanismo bajo la perspectiva
de las ecuaciones de restricción. Si el mecanismo se mueve según un la ley
zc (t), podremos escribir g ecuaciones de restricción de la forma
Φ = z(t) − zc (t) = 0 (5.96)
Estas g ecuaciones tienen como finalidad imponer a cada grado de liber-

tad la condición de que su movimiento coincida con el de la función zc (t).
Entonces, el sentido fı́sico de los g multiplicadores de Lagrange asociados a
estas restricciones es, precisamente, el de esfuerzos motores. La deducción
es lógica: el multiplicador asociado con cada restricción indica la fuerza ge-
neralizada necesaria para que se cumpla dicha restricción. Puesto que las
restricciones anteriores imponen un movimiento conocido, la fuerza corres-
pondiente es la que garantiza exactamente dicho movimiento, es decir, el
esfuerzo motor.
Un simple ejemplo puede servir para aclarar este concepto. Conside-
remos el cuadrilátero articulado de la figura 5.28, cuyo movimiento viene
definido por ϕc (t). Deseamos obtener el par motor τ (t) capaz de causar
dicho movimiento. Para ello, introducimos la nueva ecuación de restricción
Φ = ϕ(t) − ϕc (t) = 0 (5.97)

5.13. Problema dinámico inverso 183
1 2
t(t)
jc(t) B
A
Figura 5.28
La fuerza generalizada es, en este caso

   
 0
   
 0 


  1  
0   0 
 
Q = −ΦTq λ = −λ 0 = 0 (5.98)

  2  
0
  
 0 

   
 

1 ϕ −λ
lo que indica que el par motor τ es, precisamente, el multiplicador de La-

grange λ cambiado de signo.
5.13. Problema dinámico inverso

Para resolver el problema dinámico inverso, suponemos conocido el mo-
vimiento de los grados de libertad y deseamos calcular las reacciones en
algunas o en todas las articulaciones. Según hemos visto en los apartados
anteriores, el primer paso para determinar las reacciones consiste de esta-
blecer ecuaciones de restricción de par. Por ejemplo, en una articulación,
debemos imponer las restricciones x2 − x1 = 0 e y2 − y1 = 0, para que los
puntos 1 y 2 coincidan.
A partir del movimiento conocido z(t) y ż(t), resolviendo el problema
cinemático de posición podemos obtener q(t) y q̇(t). A partir de ellos,
podemos recurrir a la ecuación 5.66, que establece
� ��
M ΦTq q̈ Q
= (5.99)
Φq 0 λ −Φ̇q q̇ − Φ̇t
En este caso, podemos contemplar la expresión anterior como un sistema

lineal de ecuaciones del que podemos calcular el vector de las aceleraciones
q̈ y el vector λ de multiplicadores de Lagrange. A partir de λ, el cálculo de
las reacciones en cada par se realiza de acuerdo a lo descrito en el apartado
5.12.

�

Análisis Dinámico Por Métodos Numéricos

Cargado por

Copyright:

Formatos disponibles

Análisis Dinámico Por Métodos Numéricos

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Análisis Dinámico Por Métodos Numéricos

Cargado por

Copyright:

Formatos disponibles

Capı́tulo 5

Análisis dinámico por

que en este caso proporcionará la forma de asociar las fuerzas interiores a

5.2. Fuerzas de inercia de un elemento

Mediante una simple reordenación de términos, esta expresión se puede

donde TP es una matriz de tamaño 2 × 4, constante en el tiempo, pero

También podemos calcular el desplazamiento virtual del punto P em-

La fuerza de la inercia de la partı́cula P es sencillamente −r̈dm. El

δWP = −δrT r̈ dm = −δqTe TTP TP q̈e dm (5.5)

Para calcular el trabajo virtual de todo el elemento, deberemos sumar

En esta ecuación los términos δqe y q̈e no varı́an de una partı́cula a

El paréntesis representa una matriz simétrica de tamaño 4 × 4, que

Con la ayuda de la ecuación 5.2, que contiene la expresión explı́cita de

c Alejo Avello. Tecnun – Universidad de Navarra

explı́cita de la matriz de masas del elemento. Operando y realizando algunas

siendo Lij la longitud entre los puntos i y j, M la masa del elemento, e xG ,

Ejemplo 5.1 Calcular la matriz de masas Me del elemento de la ﬁgura

dibujamos el eje xe apuntando hacia el punto 2. En este caso:

5.3. Fuerzas de inercia de un mecanismo

Esta ecuación contiene un sumatorio extendido a cada elemento del

δWiner = −δqT M q̈ (5.16)

El proceso de reordenación anterior se conoce comúnmente como en-

c Alejo Avello. Tecnun – Universidad de Navarra

Comencemos con el mecanismo de la ﬁgura 5.3, compuesto por dos

Podemos reescribir la ecuación empleando el vector ri para referirnos

Finalmente, podemos reordenar todos los términos y escribir una ecua-

Como se puede comprobar, la matriz de masas del mecanismo se ha

Ejemplo 5.2 Calcular la matriz de masas del mecanismo de la ﬁgura 5.4,

c Alejo Avello. Tecnun – Universidad de Navarra

expresión en la que se indican el punto a que corresponden las ﬁlas y co-

Y, operando, se obtiene ﬁnalmente

Ejemplo 5.3 Calcular la matriz de masas del mecanismo biela-manivela

El vector de coordenadas naturales q está compuesto por las cuatro coor-

La matriz M1 se obtiene particularizando para el elemento 1 la ex-

Para ensamblar la matriz de masas del mecanismo, debemos distribuir

c Alejo Avello. Tecnun – Universidad de Navarra

las submatrices correspondiente a las ﬁlas y columnas del vector q. Operan-

5.4. Fuerzas exteriores puntuales

δWF = δrT F (5.20)

Anteriormente, en la ecuación 5.4, hemos calculado el valor del des-

δWF = δqTe TTP F = δqTe QF (5.21)

donde el vector QF se denomina fuerza generalizada, y de acuerdo con la

Un caso particular es la fuerza de la gravedad, en que la fuerza F se

También es frecuente el caso en que existan muelles o amortiguado-

5.5. Fuerzas introducidas por los actuadores

c Alejo Avello. Tecnun – Universidad de Navarra

Por su parte, el trabajo virtual del actuador angular se obtiene mediante

5.6. Ensamblado del vector de fuerzas generaliza-

Para comprender mejor el proceso de ensamblado del vector Q, recurri-

coinciden con los centros geométricos de las barras, se puede escribir

Ensamblando los dos vectores Q1 y Q2 en un único vector Q de fuerzas

Este proceso puede generalizarse fácilmente para cualquier tipo de fuer-

5.7. Ecuaciones del movimiento en coordenadas

δWiner + δWext = 0 (5.30)

En los apartados anteriores hemos calculado las expresiones de δWiner

−δqT Mq̈ + δqT Q = 0 (5.31)

c Alejo Avello. Tecnun – Universidad de Navarra

Cambiando de signo y sacando factor común,

δqT (Mq̈ − Q) = 0 (5.32)