UNIVERSIDAD NACIONAL DE MISIONES

CUADERNO de CÁTEDRA

MECÁNICA RACIONAL

Sergio Edgardo Katogui

Jorge Luis López
Rubén Darío Ferreyra
Víctor Javier Sommer

Año 2007
 EDITORIAL UNIVERSITARIA DE MISIONES

San Luis 1870

Posadas - Misiones – Tel-Fax: (03752) 428601
Correos electrónicos:
edunam-admini@arnet.com.ar
edunam-direccion@arnet.com.ar
edunam-produccion@arnet.com.ar
edunam-ventas@arnet.com.ar

Colección: Cuadernos de Cátedra

Coordinación de la edición: Claudio Oscar Zalazar

Armado de interiores: Amelia E. Morgenstern

Corrección: Gabriela Domínguez - Gisel Kabut

ISBN 978-950-579-074-6
Impreso en Argentina
©Editorial Universitaria

Katogui, Sergio
Mecánica racional. - 1a ed. - Posadas: EDUNaM - Editorial Universitaria
de la Univ. Nacional de Misiones, 2007.
150 p.; 30x21 cm.
ISBN 978-950-579-074-6
1. Mecánica. 2. Cinemática. 3. Dinámica. I. Título
CDD 531.112

Fecha de catalogación: 27/07/2007.

LOS AUTORES

KATOGUI, Sergio Edgardo

Ingeniero Electromecánico. Facultad de Ingeniería. UNaM.
Profesor Adjunto Regular.

LÓPEZ, Jorge Luis

Ingeniero. Facultad de Ingeniería. UNNE.
Jefe de Trabajos Prácticos Regular.

FERREYRA, Rubén darío

Ingeniero. Facultad de Ingeniería. UNaM.
Ayudante de Primera.

SOMMER, Víctor Javier

Estudiante Facultad de Ingeniería. UNaM.
Ayudante deSegunda.
 Índice

U nidad 1:Álgebra Tensorial..................................................................................9

Transformación de coordenadas ............................................................................10
Propiedades de los tensores ...................................................................................11
Operaciones con tensores.......................................................................................13
Tensores simétricos y antisimétricos .....................................................................13
Diagonalización de un tensor.................................................................................13
Diagonalización de una matriz simétrica real........................................................14
Forma cuadrática asociada.....................................................................................15
Operadores lineales................................................................................................16
Tensor de inercia....................................................................................................16
El tensor de inercia y la Transformación de coordenadas......................................17
Tensor de inercia como operador lineal..................................................................18

U nidad 2:Cinem ática de la Partícula...................................................................19

Vector de Posición, Velocidad y Aceleración .......................................................19
Coordenadas cartesianas....................................................................................... 20
Coordenadas cilíndricas........................................................................................ 21
Coordenadas esféricas ...........................................................................................23
Coordenadas intrínsecas ........................................................................................26
Sistemas Referenciales.......................................................................................... 28
Sistemas de coordenadas generalizadas................................................................ 28
Radio Vector o Vector de Posición ...................................................................... 28
Base de coordenadas generalizadas ...................................................................... 29
Coordenadas Cilíndricas....................................................................................... 29
Coordenadas Esféricas.......................................................................................... 31
Coordenadas intrínsecas ....................................................................................... 32
Transformación de Coordenadas .......................................................................... 33
Cinem ática Relativa
Sistemas de Referencia en Movimiento ............................................................... 34
Sistemas Rígidos y Deformables .......................................................................... 34
Sistemas en Traslación (únicamente) ................................................................... 35
Sistemas en Rotación............................................................................................ 36
Sistemas en Traslación y Rotación Simultánea .................................................... 40

U nidad 3: Ciném atica delCuerpo Rígido

Condición de rigidez............................................................................................. 41
Condición cinemática de velocidades................................................................... 42
Condición cinemática de aceleraciones ................................................................ 43
Movimientos del Sólido en el Espacio ................................................................. 43
Translación pura ................................................................................................... 44
Rotación pura........................................................................................................ 45
Velocidad en la rotación ....................................................................................... 45
Vector velocidad angular...................................................................................... 46
Aceleración en la rotación .................................................................................... 47
Campo de velocidades .......................................................................................... 47
Grados de libertad................................................................................................. 48
Movimiento polar del sólido................................................................................. 48
 Ángulo de Euler................................................................................................... 49
Ángulos de precesión ........................................................................................... 49
Ángulo de nutación............................................................................................... 50
Ángulo de rotación propia .................................................................................... 50
Movimiento rototraslatorio................................................................................... 51
Teorema de Euler – Chassles.................................................................................51
Polo de reducción.................................................................................................. 51
Invariante vectorial............................................................................................... 52
Invariante escalar .................................................................................................. 52
Descomposición impropia .................................................................................... 54
Descomposición propia ........................................................................................ 54
Eje central del movimiento................................................................................... 55
Ecuación del eje central ........................................................................................ 56
Movimiento helicoidal tangente............................................................................ 57
Axoides................................................................................................................. 57
Movimiento paralelo............................................................................................. 57
Movimiento polar ................................................................................................. 58
Movimiento Plano ................................................................................................ 58
Centro instantáneo de rotación ............................................................................. 59
Trayectorias polares.............................................................................................. 60
Propiedades del Centro instantáneo de rotación ................................................... 60
Ecuación de la base............................................................................................... 61
Ecuación de la rodante.......................................................................................... 62
Centro instantáneo de velocidad........................................................................... 63
Estado de aceleración en el movimiento plano..................................................... 63
Polo de aceleraciones .............................................................................................64
Coordenadas del polo de aceleraciones................................................................. 64

U nidad 4: D inám ica de la Partícula

Principio de Inercia (1º Ley de Newton) ...............................................................67
Principio de masa (2º Ley de Newton)..................................................................67
Principio de Acción y Reacción (3º Ley de Newton)........................................... 67
Principio de independencia de acción de fuerzas.................................................. 67
Movimiento unidimensional de la partícula ......................................................... 68
Fuerza dependiente del tiempo ............................................................................. 69
Fuerza dependiente de la velocidad.......................................................................70
Fuerza dependiente de la Posición........................................................................ 70
Definición de energía potencial.............................................................................70
Curvas de potencial............................................................................................... 71
Estados de Equilibrio............................................................................................ 72
Movimiento Curvilíneo de la Partícula................................................................. 74
Cantidad de movimiento....................................................................................... 75
Momento cinético ................................................................................................. 76
Teorema del Momento cinético .............................................................................76
Vínculos o Ligaduras............................................................................................ 77
Reacciones de vínculos......................................................................................... 78
Vínculos Rugosos.................................................................................................. 79
Movimiento sobre una curva lisa............................................................................79
Movimiento sobre una superficie lisa......................................................................80
 Dinámica Relativa de la Partícula......................................................................... 80
Principio de D´Alembert....................................................................................... 81

U nidad 5: D inám ica de los Sistem as de Partículas

Masa del sistema................................................................................................... 83
Cantidad de movimiento....................................................................................... 84
Conservación de la cantidad de movimiento ........................................................ 85
Momento cinético ..................................................................................................86
Conservación del momento cinético..................................................................... 86
Trabajo y Energía ................................................................................................. 87
D inám ica delSólido Rígido
Masa del sólido rígido .......................................................................................... 89
Cantidad de movimiento........................................................................................ 89
Momento cinético ..................................................................................................90
Momento y productos de inercia .......................................................................... 91
Conservación del momento cinético......................................................................93
Diagonalización del tensor de Inercia................................................................... 94
Momento cinético................................................................................................. 95
Teorema de los ejes paralelos................................................................................ 95
D inám ica relativa delsólido
Cantidad de movimiento........................................................................................ 96
Momento cinético .................................................................................................. 97
Reacciones dinámicas de un sólido en rotación ................................................... 98
Equilibrio dinámico .............................................................................................100
Movimiento de un sólido con un punto fijo ........................................................100
Movimiento de un sólido por Inercia ..................................................................101

U nidad 6: Sistem as vibratorios de un grado de libertad

Fuerza recuperadora elástica................................................................................105
Modelación dinámica...........................................................................................105
Deformación de sistemas elásticos ..................................................................... 106
Combinación de constantes elásticas.................................................................. 109
Estudio de Vibraciones Mecánicas Armónicas .................................................. 110
Vibraciones armónicas libres no amortiguadas .................................................. 111
Oscilador Amónico Torsional Libre................................................................... 112
Oscilaciones con influencia del peso................................................................. 112
Métodos energéticos ........................................................................................... 113
Vibraciones Armónicas amortiguadas................................................................ 114
Vibraciones Subamortiguadas ............................................................................ 115
Vibraciones Sobreamortiguadas .........................................................................116
Vibración Críticamente Amortiguada................................................................. 116
Factor o grado de amortiguamiento.....................................................................117
Decremento Logarítmico.................................................................................... 117
Vibraciones Armónicas Forzadas ....................................................................... 118
Vibraciones Forzadas sin Amortiguamiento ...................................................... 119
Vibraciones Forzadas con Amortiguamiento ..................................................... 120
Representación fasorial....................................................................................... 124
Vibración armónica Libre................................................................................... 124
Vibración forzada no amortiguada ..................................................................... 124
Vibración Forzada con Amortiguamiento .......................................................... 125
 Desbalance rotatorio ........................................................................................... 126
Movimiento del soporte...................................................................................... 127
Aislamiento Vibratorio ....................................................................................... 130

U nidad 7: O scilaciones de dos o m ás G rados de Libertad

Modos normales de vibración............................................................................. 133
Valores y vectores propios ................................................................................. 137
Vibración armónica forzada ............................................................................... 138
Amortiguador de vibraciones ............................................................................. 139

U nidad 8: M ecánica analítica

Principio de los desplazamientos virtuales ......................................................... 141
Trabajos virtuales................................................................................................ 141
Ecuación general de la dinámica......................................................................... 142
Coordenadas generalizadas................................................................................. 142
Fuerzas Generalizadas ........................................................................................ 143
Ecuaciones de Lagrange ..................................................................................... 144
Dinámica impulsiva ............................................................................................ 145
Índice Alfabético..................................................................................................147
Bibliografía......................................................................................................... 149
Unidad 1: Álgebra Tensorial

Introducción

Se estudiará en esta unidad el álgebra de las funciones vectoriales lineales, o

tensores, instrumentos matemáticos que serán utilizados para el estudio de la dinámica de
los sólidos.
Para introducir el concepto de tensor se partirá de un ejemplo: se sabe que cuando un
sólido gira alrededor de un eje de simetría, el vector momento angular L tiene la misma
dirección que el vector velocidad angular ω, y la relación entre ambos está dada por:
L=Iω
donde I representa el momento de inercia del sólido respecto del eje de rotación.
Por ejemplo, consideremos el caso de un sólido
ω cilíndrico que rota sobre su eje.
En este caso, I = Momento de inercia del cuerpo
L respecto del eje de rotación.
En el dibujo se destaca la colinealidad entre el
momento angular y la velocidad angular.
La colinealidad entre los vectores L y ω no se
verifica en general, sino que se trata de un caso
particular, con las condiciones mencionadas.

A continuación se analiza un caso general.

En la figura de la derecha, se muestra un sólido no L
regular que rota sobre un eje cualquiera, y se representa un
caso general en que no existe colinealidad entre los
vectores L y ω. ω
Sabemos, sin embargo, la íntima relación entre estas
magnitudes y, por otra parte, se puede verificar que a cada
vector ω le corresponde un vector L perfectamente
definido.
Podemos pensar entonces en una operación lineal tal
que realizada sobre ω da por resultado L. Esto puede
representarse mediante la siguiente ecuación operacional:
L = [I] ω
[I] representa el operador, que “transforma” a ω en L.
Se verá que [I] se representa por medio de una matriz de 3x3, y cada componente de
esa matriz depende de la forma en que está distribuida la masa del sólido, y del sistema
referencial en que se expresa.
De la ecuación operacional L = [I] ω podemos “despejar” [I] = L /ω ω. Esta operación
(el cociente entre dos vectores) no es factible de realizar y la indicamos al solo efecto de
introducir el siguiente razonamiento: el cociente entre dos cantidades puede ser de tipo
diferente y “más complicado” que el de estas; por ejemplo, el cociente de dos enteros
puede ser un número racional. De modo análogo, el cociente entre vectores no es un
vector. No es extraño entonces que [I] sea una magnitud de un nuevo tipo: un tensor de
segundo orden. Este tensor se denomina tensor de inercia.

9
 Otro ejemplo de tensor es el tensor producto tensorial de dos vectores. Dados los
vectores a:[a1,a2,a3], y b:[b1,b2,b3], se define el producto tensorial a⊗b de la siguiente
manera:
 a1b1 a1b2 a1b3 
a ⊗ b = a 2 b1 a 2 b2 a 2 b3 
 a3b1 a3 b2 a3b3 

Aquí se ha representado el tensor producto tensorial en forma matricial. No debe

confundirse la representación del tensor con el tensor mismo, se verá que, en distintos
sistemas referenciales, el tensor adoptará distintas representaciones.

Representación con subíndices

Si bien se trabajará con la representación matricial, los tensores también se

pueden representar utilizando subíndices, así, el tensor producto tensorial se puede
representar como aibj=[a1b1,a1b2,.....,a3b2,a3b3].
Para indicar en forma compacta sumatorias de términos, se adopta la siguiente
convención: “Cuando en una expresión monomia, figuran dos índices repetidos, se
entenderá que se trata de una suma en la que los índices repetidos van sumados de 1 a la
dimensión en que se está trabajando”.
Por ejemplo:
n
aibi = ∑ ai bi = a1b1+a2b2+.....+anbn
i =1
n
aibjci= ∑ ai b j ci =a1bjc1+a2bjc2+...+anbjcn
i =1

Transformación de coordenadas

Las componentes aibj son las componentes del tensor en el sistema referencial en el
cual se representan los vectores a y b. En un sistema referencial distinto, ligado al primero
por la matriz de transformación o de cambio de base [M], las componentes de tensor serán
distintas, es decir, la representación del tensor será distinta.

Con convención de subíndices la transformación de un vector se escribe:

a’m = hmiai
donde hmi son los elementos de la matriz de transformación.

Los vectores a y b se transforman el nuevo sistema referencial utilizando la matriz de

transformación:
a’ = [M]a
b’ = [M]b
Se entiende que en las operaciones indicadas los vectores se escriben como columnas,
para realizar la aplicación de la matriz de transformación a los mismos.
En la nueva base, se realiza el producto tensorial:

10
  a1b1 a1b2 a1b3 
a ⊗ b = a 2 b1 a 2 b2 a 2 b3 
 a3b1 a3 b2 a3b3 

Los elementos del tensor producto tensorial [a⊗b] se transforman como sigue:
[a’⊗b’] = hmiaihmjbj
Utilizando el álgebra matricial puede transformarse el tensor de un sistema referencial
a otro utilizando la matriz de cambio de base [M]:
a’⊗b’=[M][a⊗b][M]T
donde [M]T indica la transposición de la matriz de cambio de base, equivalente a la
inversión cuando las matrices son ortonormales.

Generalización del concepto de tensor

Los tensores son entes matemáticos que pueden representar propiedades o

características de cuerpos o materiales. Además del tensor de inercia, podemos mencionar
el tensor de tensiones, cuyos elementos son tensiones normales y tensiones tangenciales, y
que representa el estado tensional de un sólido.
Ante una transformación de coordenadas puede cambiar la representación del tensor,
el cual conservará propiedades que no cambian cuando se cambia el sistema referencial. Se
puede generalizar el concepto de tensor, diciendo que tmn elementos constituyen un tensor
si las componentes del mismo en dos sistemas referenciales distintos están ligados por:
t’mn = hmihnjtij
Cuya notación en forma matricial es:
[t’] =[M][t][M]T

Orden de los tensores

Los tensores hasta aquí presentados a modo de ejemplo son tensores de segundo
orden y su representación se realiza mediante matrices. Los escalares son tensores de
orden 0, en tanto que los vectores son tensores de orden 1.
El número de componentes de un tensor se determina elevando la dimensión del
espacio de representación del tensor al orden del tensor. Llamando n a la dimensión del
espacio, y p al orden del tensor, resulta:

nº de componentes = np

En esta asignatura se utilizaran tensores de orden 2, específicamente, el tensor de

inercia, para representar características de un sólido ligadas a la forma en que está
distribuida su masa.

Propiedades de los tensores:

Un tensor posee propiedades independientes del sistema de representación,

denominadas invariantes.
Ejemplos de invariantes:
1. Dado el tensor de orden 0 o escalar, el mismo no cambia en distintos sistemas
referenciales.
2. Dado el tensor de primer orden o vector, a = [a1,a2,a3]

11
 El módulo de a es un invariante ante cambios de coordenadas.
Dados los vectores a y b
El producto escalar de los mismos es un invariante.
El ángulo entre los mismos es un invariante.
En general, las operaciones entre invariantes de tensores da nuevos invariantes.
3. Para un tensor de segundo orden distinguimos los siguientes invariantes:
Dado el tensor [t], y su representación en dos sistemas referenciales:

t11 t12 t13  t '11 t '12 t '13 

t t 23  t ' t ' 23 
[t]=  21 t 22 ; [t’]=  21 t ' 22
t 31 t 32 t 33  t '31 t '32 t '33 
Podemos distinguir los siguientes invariantes.
1. La suma de los componentes de la diagonal t11+t22+t33 = t’11+t’22+t’33=k
2. Dado un tensor de segundo orden, del mismo se puede deducir un tensor de primer
orden o vector: v = [(t23-t32),(t31-t13),(t12-t21)]. Es evidente que en el caso de tensores
simétricos, este vector será nulo.
3. A un tensor de segundo orden le corresponde un escalar invariante dado por su
determinante, k = Det[T]

Ejemplos de invariantes de tensores de segundo orden

1. Tensor producto tensorial
 a1b1 a1b2 a1b3 
Dados los vectores a y b, su producto tensorial es: a⊗b= a b a 2 b2 a 2 b3 
 2 1
 a3 b1 a3b2 a3b3 
Los invariantes son:
a. Invariante escalar = a1b1 + a2b2+a3b3 = a•b (producto escalar entre a y b)
b. El vector v = [(a2b3-a3b2),( a3b1-a1b3), ( a1b2-a2b1)] es el producto vectorial a x b.
c. El tensor derivado:
Dado el campo vectorial B(x,y,z)=[B1(x1,x2,x3),B2(x1,x2,x3),B3(x1,x2,x3)], siendo las
Bi(x,y,z) continuas y derivables, las nueve derivadas parciales constituyen un tensor,
denominado tensor derivado:

 ∂B1 ∂B1 ∂B1 

 
 ∂x1 ∂x 2 ∂x3 
∂ ( B1 , B2 , B3 )  ∂B2 ∂B2 ∂B2 
∂ ( x1 , B2 , B3 )
=  ∂x ∂x 2 ∂x3 
 ∂B1 ∂B3 ∂B3 
 3
 ∂x1 ∂x 2 ∂x3 

En este tensor los invariantes son:

a. Invariante escalar:
La suma de los elementos de la diagonal constituyen la divergencia de B:
∂B ∂B ∂B
∇•B = 1 + 2 + 3
∂x1 ∂x 2 ∂x3

12
 b. Invariante vectorial:
 ∂B2 ∂B3   ∂B3 ∂B1   ∂B1 ∂B2 
El vector  − ,  − ,  −   es el rotor de B.
 3∂x ∂x 2   ∂x1 ∂x 3   ∂x 2 ∂x1 
Operaciones con tensores:

1. Adición y sustracción: esta operación está definida para tensores del mismo orden,
dados dos tensores tmn y umn, su suma es otro tensor smn cuyas componentes son las sumas
de las componentes del mismo subíndice: smn = tmn + umn. En forma matricial: [s] = [t]+[u].
2. Multiplicación por un escalar: se define la multiplicación de un tensor por un
escalar, al nuevo tensor que resulta de multiplicar cada componente del tensor por dicho
escalar: smn = k tmn, o [s]=k [t]
3. Producto tensorial: el producto de un tensor de orden p por otro de orden q, da por
resultado un tensor de orden p+q, cuyas componentes son el producto de las componentes
del primero por las componentes del segundo.

Tensores simétricos y antisimétricos:

Un tensor de segundo orden es simétrico si tij = tji, y es antisimétrico si tij = -tji. Las
propiedades de simetría y antisimetría son independientes del sistema referencial. Todo
tensor puede expresarse como la suma de un tensor simétrico y un tensor antisimétrico:
tij= ½(tij+tji)+ 1/2(tij-tji)
En forma matricial: [T] = ½ ([T]+[T]T)+1/2([T]-[T]T).

Diagonalización de un tensor

En esta sección se trabajará con la diagonalización de un tensor, es decir, dado un

tensor se buscará la manera de encontrar un sistema referencial en el cual su
representación sea diagonalizada. Esta operación implica la determinación de los vectores
y valores característicos (autovectores y autovalores) de la matriz que representa al tensor.

Se considera una transformación lineal de la forma a’m=hmi.ai, que en forma matricial

se indica [v’]T = [H][v]T donde [H] es la matriz de transformación.
Se buscarán ahora vectores tales que al aplicarles el operador [H], se transforman en
vectores colineales. Significa esto que el vector transformado es múltiplo escalar del vector
dado, v’ = λv

[H][v]T = λ[v]T = λ [1][v]T [1] : Matriz identidad

que puede reescribirse:

[H][v]T - λ [1][v]T= ([H] - λ [1] )[v]T = 0

Desarrollando los productos, queda el siguiente sistema de ecuaciones:

 (h11 − λ )v1 + h12 v 2 + h13 v3 = 0


h21v1 + (h22 − λ )v 2 + h23 v3 = 0
 h v + h v + ( h − λ )v = 0
 31 1 32 2 33 3

13
 Siendo este sistema homogéneo tendrá soluciones no triviales si el determinante de
los coeficientes es nulo, o sea:

(h11 − λ ) h12 h13

h21 (h22 − λ ) h23 =0
h31 h32 (h33 − λ )

Al resolver el determinante, queda un polinomio de tercer grado:

λ3 + Aλ2 + Cλ +D = 0
Esta ecuación se denomina ecuación característica y sus raíces, λ1, λ2, λ3 se
denominan valores característicos, autovalores o eigenvalues.

Reemplazando en el sistema de ecuaciones los autovalores (de a uno por vez), se

obtienen los vectores v1 (v11,v12,v13), v2 (v21,v22,v23), v3 (v31,v32,v33) denominados vectores
característicos, vectores propios, autovectores o eigenvectores.

Se verifican los siguientes corolarios y teoremas:

a. Los autovalores de una matriz de componentes real y simétrica son reales.

b. Un autovector no puede corresponder a dos autovalores distintos.
c. Si v1, v2, v3 son autovectores que corresponden a valores propios distintos, entonces
son linealmente independientes.
d. Si vi, vj son autovectores que corresponden a los autovalores λi, λj de una matriz
simétrica real, entonces son ortogonales.

Diagonalización de una matriz simétrica real.

Sean λ1, λ2, λ3 autovalores de [A], matriz real y simétrica que representa a un tensor en
un sistema referencial ortogonal y v1,v2,v3, los autovectores correspondientes a los λ1, λ2,
λ3.
Dividiendo cada autovector por su módulo, y ordenándolos en forma adecuada, se
forma con ellos una base ortogonal directa. Se forma la matriz [M], tal que sus filas sean
los autovectores divididos por su módulo:

 v11 v12 v13 

[M]= v 21 v 22 v 23  [M]: Matriz Modal

v31 v32 v33 
Esta matriz [M] es la matriz de transformación que liga el referencial original con un
sistema referencial cuya base son los autovectores.
Premultiplicando [A] por [M], y postmultiplicandolo por [M]T, resulta una matriz
diagonal, siendo los elementos de la diagonal los autovalores λ1, λ2, λ3.
λ1 0 0
 λ2 0 
[M][A][M]T=[D] =  0
 0 0 λ3 
Esta matriz diagonal es la representación del tensor en el referencial cuya base son los
autovectores.
14
Casos particulares

Los tensores simétricos pueden ponerse en forma diagonalizada mediante una

transformación ortogonal apropiada. Los elementos de la diagonal son únicos, aunque
pueden no aparecer en el mismo orden en la diagonal de la matriz. Los ejes del sistema
referencial en el cual el tensor tiene representación diagonalizada se denominan ejes
principales. Si los autovalores son distintos estos ejes son únicos y ortogonales.
Como casos particulares pueden aparecer autovalores iguales. Si aparecen dos autovalo-
res iguales se dice que el autovalor es doblemente degenerado, si los tres autovalores
son iguales, el mismo es triplemente degenerado.
En el caso de triple degeneración el tensor tiene la misma representación el todo
sistema referencial, y cualquier eje es eje principal. Como ejemplo, el tensor de inercia de
una esfera de densidad uniforme es diagonalizado en todo sistema referencial cuyo origen
es el centro de la esfera y cualquier eje que pase por su centro es un eje principal de
inercia.
En el caso de doble degeneración, el eje principal asociado al autovalor distinto es
único, y perpendicular a un plano. Cualquier eje que esté contenido en este plano y corte al
eje principal único es eje principal.

Forma cuadrática asociada

Se llama Forma cuadrática a una expresión de segundo grado en “n” variables. Por
ejemplo, para n = 3, Q(x1,x2,x3) = a11x12+a12x1x2+a13x1x3+….+a21x2x1+a33x32.
La forma cuadrática en xi puede ser expresada en forma matricial:

Q[x] = [x][A][x]T

 a11 a12 a13 

[x] = [x1, x2, x3], A = a 21 a 22 a 23 
 a31 a32 a33 

Si la matriz [A] es la representación de un tensor, Q[x]= [x][A][x]T es la forma

cuadrática asociada al tensor.
La cuadrática asociada a un tensor es un invariante ante cambios de sistemas
referenciales. En particular, para un sistema referencial cuya base asociada son los
autovectores la representación del tensor será diagonal, y la forma cuadrática se reduce a:

Q[x] = λ1x12+ λ2x22+ λ3x32

Definiendo la siguiente operación:

[x][A][x]T=1

y considerando que [A] es un tensor simétrico de segundo orden, se obtiene:

a11x12+ 2a12x1x2 + a 22x22+2a13x1x + 2a23x2x3+a33x32 = 1

15
 Se trata de una ecuación de segundo grado que define una superficie de segundo
orden en el espacio tridimensional. Como no tiene términos lineales se trata de un
elipsoide o un hiperboloide de una o dos hojas. Esta superficie se denomina superficie
característica o cuádrica característica, y permite representar un tensor en forma gráfica.
En el sistema referencial para el cual el tensor es diagonalizado, la forma cuadrática
se reduce a:

λ1x12 + λ2x22 + λ3x32 = 1

En el caso del tensor de inercia, dado que los autovalores son momentos de inercia,
siempre positivos, la cuádrica es siempre un elipsoide, que en el caso de doble
degeneración es un elipsoide de revolución, y en el caso de triple degeneración es una
esfera.

Operadores lineales

Se considera la siguiente operación: c = [A] b

Donde [A] es una matriz, c y b son vectores. En esta operación la matriz puede ser
considerada como una transformación lineal, o bien como un operador lineal, que al vector
b le hace corresponder el vector c.

Anexo

El tensor de Inercia

Se estudiará un ejemplo de tensor: el tensor de Inercia. Este tensor de segundo orden

resume las características de un sólido, relativas a la cantidad de masa que posee, y a la
forma en que ésta está distribuida. Su importancia para esta asignatura radica en la
utilización del mismo en el estudio de la dinámica de rotación de los sólidos.
Los componentes del tensor de inercia de un sólido referido a un sistema referencial
“xyz”, son los momentos de inercia del sólido respecto de los ejes, elementos que se
encuentran en la diagonal, y los productos de inercia.

 I xx
I xz  I xy
Tensor de inercia: [I] =  
 I yx
I yz  I yy
 I zx
I zy I zz 

Donde: Ixx, Iyy, Izz son los momentos de inercia, Ixy = Iyx, Ixz=Izx, Iyz=Izy son los
productos de inercia.
Estos momentos y productos de inercia, para distribuciones continuas de masa, se
calculan de acuerdo a las siguientes expresiones:

∫∫∫ (y )
+ z 2 ρ [x, y, z ]dV
2
Ixx =
V

∫∫∫ (x + z2 )ρ [x, y, z ]dV

2
Iyy = V

∫∫∫ (x )
+ y 2 ρ [x, y, z ]dV
2
Izz = V

Ixy=Iyx = − ∫∫∫V xyρ [x, y, z ]dV

16
 Ixz =Izx= −
∫∫∫ xzρ [x, y, z ]dV
V

Iyz =Izy= − ∫∫∫V yzρ [x, y, z ]dV

En estas definiciones para los momentos y productos de inercia, se tiene.

ρ [x, y, z ] = densidad de masa V = volumen ocupado por el sólido

x, y, z son las coordenadas en el sistema “xyz” del elemento de volumen “dV”

ρ [x, y, z ] dV es la masa (dm) del elemento de volumen “dV”

Para sistemas formados por mi masas puntuales, las definiciones para los momentos de
inercia son las siguientes:

Ixx= ∑(y 2
i + z i2 )mi
Iyy= ∑ (x i
2
+ z i2 )mi
Izz= ∑ (x 2
i + yi2 )mi

Ixy=Iyx = − ∑ ( xi y i )mi

Ixz =Izx= − ∑ ( xi z i )mi

Iyz =Izy= − ∑ ( y i z i )mi

Donde xi,yi,zi son las coordenadas en el sistema “xyz” de la masa puntual mi.

En base a las definiciones de los elementos que componen el tensor de inercia, se

puede afirmar que el tensor de inercia es simétrico, y sus coeficientes son números reales.

El tensor de inercia y la transformación de coordenadas

Sea el tensor de inercia [I], referido a un sistema referencial “xyz”. Si se considera un

sistema referencial x’y’z’, rotado respecto del primero, el tensor de inercia [I] referido a
este sistema referencial es el mismo, pero su representación cambiará. Significa esto que
las componentes del tensor tendrán valores distintos. La relación entre ambas
representaciones del tensor se obtienen mediante las ecuaciones de transformación:
[I]’ = [R][I][R]T
Donde [R] es la matriz modal o matriz de cambio de base, que relaciona
representaciones en xyz con representaciones en x’y’z’.

Para cada tensor de inercia existe por lo menos un sistema referencial en el cual la
representación de la matriz es diagonalizada. Para determinar estas direcciones puede
utilizar el procedimiento descrito en la sección correspondiente, las direcciones que
coinciden con las de ese sistema referencial se denominan ejes principales de inercia.
Ejemplos: para un sólido cilíndrico de densidad de masa uniforme, el tensor de inercia
referido a cualquier sistema referencial de ejes ortogonales tal que su origen esté ubicado

17
en el centro de gravedad, y tenga un eje coincidente con el eje del cilindro, es diagonal
izado.

El tensor de inercia como operador lineal

Para un sólido en rotación, el momento angular del mismo puede obtenerse aplicando
como un operador lineal el tensor de inercia al vector velocidad angular.
L =[ I] ω
Si ω tiene la dirección de un eje principal de inercia, es un autovector del operador,
es decir, para este caso L será colineal con ω
Según lo visto en Física 1, la variación del momento angular de un cuerpo, expresada
como su derivada respecto del tiempo, es igual al torque o momento externo que actúa
sobre el cuerpo.

18
Unidad 2: Cinemática de la Partícula
Conceptos generales:

La cinemática estudia el movimiento de los cuerpos con prescindencia de las causas

que lo producen. Describe la posición, la velocidad, la aceleración y la trayectoria que
puede tener un cuerpo, sin considerar las fuerzas que lo ocasionan.
En esta unidad se describirá el movimiento de partículas puntuales, sin dimensiones ni
estructura interna. Partícula o punto es una abstracción que permite simplificar el estudio
del movimiento y desarrollar sus conceptos fundamentales. En determinadas circunstancias
y condiciones algunos sistemas mecánicos se pueden abordar razonable y suficientemente,
considerándolos como “puntuales”.
Para estudiar el movimiento resulta indispensable adoptar sistemas de referencia
convenientes. Así, se tienen los Sistemas de Referencia “Fijos” o “Absolutos” y sistemas
de Referencia móviles o “Relativos”.

Vector de Posición, Velocidad y Aceleración

El Vector de posición de la partícula P, respecto del sistema fijo “O” es r(t), función
vectorial que depende del tiempo. El Vector de posición tendrá siempre su origen
coincidente con el origen del sistema fijo, en tanto que el otro extremo coincidirá o seguirá
la posición de la partícula.
La trayectoria de un móvil es la curva alabeada c descripta paramétricamente por la
función r(t).
Existen muchos tipos de movimientos diferentes, incluso sobre una misma trayectoria.
Estos se distinguen por el ritmo con que cambia el vector de posición al transcurrir el
tiempo.
Se ve en la figura la partícula en dos instantes diferentes, “t” y poco tiempo después,
“t + ∆t”. Durante ese intervalo de tiempo, la posición pasa de ser descripta por el vector
r[t], al vector r[t+∆t]
En el dibujo se han omitido los ejes del sistema referencial.

El cambio de posición que ha

experimentado la partícula, está dado por el
∆r vector desplazamiento ∆r:
∆r = r[t+∆t] – r[t]
r[t] Generalmente, este vector no coin-
c
cide con la trayectoria, salvo el caso par-
r[t+∆t] ticular del movimiento rectilíneo.
O
La velocidad, o ritmo con que cambia o
se modifica el vector posición, es la derivada
del vector posición respecto del tiempo, límite del cociente incremental ∆r/∆t cuando el
incremento “∆t” tiende a cero:
r [t + ∆t ] − r [t ] dr [t ]
v = Lím =
∆t → 0 ∆t dt
En el límite, cuando el intervalo de tiempo se hace cada vez cada vez más pequeño, el
vector ∆r se aproxima a la trayectoria y el vector velocidad queda tangente a la misma.

19
 Este vector velocidad dependiente del tiempo se indicará de la siguiente manera:
dr(t ) • r( t + ∆t ) − r(t )
v (t ) = = r( t ) = lim
dt ∆t

Siendo r(t) el vector de posición referido a un sistema fijo, la velocidad v(t) definida
es la velocidad absoluta, y es independiente del origen “O” considerado.

La aceleración se define como:

dv (t ) d 2 r(t ) •• v ( t + ∆t ) − v ( t )
a(t ) = = = r( t ) = lim
dt dt 2
∆t
Podemos entender a la aceleración como la causa de los cambios en la velocidad. Más
adelante se verá que el cambio en el módulo del vector velocidad se deben a una
componente de la aceleración denominada aceleración tangencial, en tanto que los cambios
en la dirección se deben a la aceleración normal.

Interesa analizar r(t), v(t) y a(t) para los sistemas de coordenadas más usuales o
prácticos: cartesianas, cilíndricas, esféricas, intrínsecas. Para cada uno de estos sistemas de
coordenadas se define una base en la cual se expresan los vectores. Para más detalles
acerca de los sistemas referenciales y los versores que conforman sus bases, consultar el
Anexo 1, y la bibliografía recomendada.

En coordenadas cartesianas:

Se toma un sistema de coordenadas xyz, de base ortonormal {i,j,k}, con origen en el

punto
fijo “O”. El vector de posición r(t) tomará la forma:
r(t) = x(t) i + y(t) j + z(t) k
La velocidad es:
• dz • • •
dx dy
v(t)= r(t) = i+ j + k = x(t) i + y(t) j+ z(t) k
dt dt dt

La aceleración resulta:
• •• d 2x d 2y d 2z •• •• ••
a(t)= v(t) = r(t) = i+ j+ k = x(t) i + y(t) j+ z(t) k
dt 2 dt 2 dt 2
Nota: se considera el referencial de coordenadas cartesianas como fijo, es decir, que
los versores i,j,k no varían en el tiempo y tienen derivada nula.

20
En coordenadas cilíndricas:

Se utilizan las coordenadas (ρ, φ, z) y la base móvil ortonormal {eρ,eφ,k}.

El vector de posición es:
r(t)= ρ(t) eρ+z(t)k
z
k
z eφ
ρ
eρ

r(t) z(t)k

φ y
ρ(t)eρ
ρ
x

La velocidad se obtiene derivando r(t) respecto del tiempo. Al realizar esta derivada,
se tiene en cuenta que el vector de posición está referido a una base móvil, es decir, que los
versores que la forman cambian su dirección en el tiempo.
• • •
dr(t)
= ρ eρ + ρ ρ
e
v(t)= + z k
dt •
e
En la ecuación anterior, ρ representa la derivada del versor eρ respecto del
tiempo.
Siendo:
eρ = Cos (φ) i + Sen(φ) j, y eφ = -Sen (φ) i + Cos(φ) j
• • •
e ρ = φ (-Sen (φ) i + Cos(φ) j) = φ e
φ

Es decir, que la derivada del versor eρ respecto del tiempo tiene la dirección de eφ.
Utilizando esta relación se reformula v(t):
• • •
v(t) = ρ eρ + ρ φ eφ + z k

En la descripción de la velocidad en coordenadas cilíndricas aparecen tres

componentes:
•
ρ , componente de la velocidad en la dirección de e , será distinto de cero cuando el
ρ
movimiento sea tal que ρ varíe en el tiempo.
•
ρ φ es la componente de la velocidad en la dirección de eφ, y

21
 •
z la componente de la velocidad en la dirección de k.

Se remarca que, en un caso general de movimiento en el espacio, el vector velocidad

v(t) siempre tiene la dirección de la tangente a la trayectoria. {eρ,eφ,k} forman un triedro
ortogonal, cuya orientación cambia “acompañando” al móvil.

Un caso particular es el movimiento circular, si el sistema de coordenadas se elige de

manera tal que su origen coincide con el centro de la circunferencia descrita por el móvil, y
el eje z tal que resulte perpendicular al plano en que se realice el movimiento, resultarán
• • •
z = 0 y ρ = constante = radio de la circunferencia; ρ = z =0; φ = ω (velocidad
angular);
ρ ω = velocidad tangencial (v ). La expresión de la velocidad será:
T
ρ
v(t) = ω eφ = vT eφ
vT es la conocida velocidad tangencial del movimiento circular.

•• • • • • ••
•• • •
dv(t) ρ ρ e ρ ρ φ ρ φ ρ φ e φ
La aceleración es: a(t) = = eρ + + eφ + eφ + +z k
dt

• • •
será φ = - φ ( Cos (φ) i + Sen(φ) j) = - φ eρ
e
Siendo eφ = -Sen (φ) i + Cos(φ) j,
Agrupando términos:
•• • •• • • ••
a = ( ρ - ρ φ 2) eρ + ( ρ φ +2 ρ φ ) eφ + z k

Se analizarán a continuación los distintos componentes de la expresión de la

aceleración en coordenadas cilíndricas:

Componentes en eρ:
••
ρ : Puede interpretarse como una aceleración radial.
•
ρ φ 2: Es la causante del cambio de dirección de la componente de velocidad en e . En
φ
el caso de un movimiento circular, es la aceleración centrípeta.

Componentes en eφ:
••
ρ φ : Es originada por la variación de la velocidad angular con que rota el radio ρ . Si
el movimiento es circular, es la aceleración tangencial.

• •
2 ρ φ : Es la aceleración complementaria (ac), aparece únicamente cuando tanto el
radio ρ como el ángulo φ varían en el tiempo. Proviene de dos contribuciones simultáneas
• •
de igual valor ρ φ pero de diferente origen:

22
 • •
1- Debido al cambio de dirección de ρ . Suponiendo que el movimiento posee ρ ≠0, y
• •
de módulo constante. Si φ ≠0, la dirección del vector ρ cambia contínuamente, como se
muestra en la figura:

En la figura, se aprecia que:

• • • • •
∆ρ
∆ ρ =∆φ ρ 1 =∆φ ρ 2 =∆φ ρ
•
• Formando el cociente ∆ ρ = ∆φ •
ρ2 ∆φ ρ
∆t ∆t
•
ρ1
y en el límite cuando ∆t→0:
• •
(½)ac= ρ φ ac = Aceleración complementaria
• •
2- Debido al desplazamiento sobre el radio: Suponiendo φ =constante, si ρ ≠0, el
valor de ρ se modifica, lo que producirá un cambio en la componente de la velocidad en la
dirección de eφ
•
∆v = v2 – v1 = φ (ρ2-ρ1)
v2 ∆v
• ρ 2 − ρ1
Formando el cociente = φ
v1 ∆t ∆t

• •
ρ1 ρ2
En el límite cuando ∆t→0 b( 1/2 )ac = φ . ρ

En coordenadas esféricas:
Se utilizan las coordenadas generalizadas esféricas ρ,θ,λ y la base ortonormal
{er,eθ,eλ}
El vector de posición se expresa en la base de coordenadas esféricas como:
r = r er

23
 z

er

eλ

eθ
θ θ r.er

y
λ λ

Para las coordenadas esféricas se cumplen las siguientes relaciones entre los versores
de la base:
er = Sen(θ)Cos(λ) i + Sen(θ )Sen(λ) j + Cos (θ)k

∂e r = Cos(θ)Cos(λ) i + Cos(θ )Sen(λ) j - Sen (θ)k = eθ

∂θ
∂e r
= - Sen(θ)Sen(λ) i + Sen(θ )Cos(λ) j = Sen(θ)eλ
∂λ
El vector velocidad en coordenadas esféricas resulta:
• • •
dr
v = = r er + r θ eθ + r λ Sen(θ)eλ
dt
En coordenadas esféricas, las curvas a r y λ constantes se denominan meridianos; las
curvas a r y θ constantes son los paralelos.
En la expresión de la velocidad, aparecen tres componentes:
•
r : es la componente en la dirección del radio r. Será siempre nula en los movimientos
con r constantes.
•
r θ : componente en la dirección de eθ, tangente a un meridiano.
•
r λ : componente en la dirección de eλ, tangente a un paralelo.

24
 z

er

eλ
r
θ
eθ
y

La aceleración resulta:
dv • ••
•• • • • •
a= dt = ( r - r – r Sen (θ) ) er + (r + 2 r θ - r Sen(θ)Cos(θ) λ 2) eθ +
θ 2 2 λ 2 θ
•• • • • •
+(r Sen(θ) λ +2 r Sen(θ) λ +2 r Cos(θ) θ λ )eλ
Las componentes en er:
••
r : Aceleración en la dirección radial, originada en el cambio de módulo de la
componente de la velocidad en er.
•
r θ 2: Es una aceleración centrípeta debida a la rotación sobre el meridiano.
•
r Sen (θ) λ 2: Proyección sobre er de la aceleración centrípeta debido a la rotación
2

sobre el paralelo.

• •
r Sen(θ) λ 2 r Sen(θ)Cos(θ) λ 2

Paralelo θ
•
r Sen2(θ) λ 2
•
r Sen2(θ) λ 2
θ

Componentes en eθ

25
 ••
r θ : aceleración tangencial debida a rotación en meridiano.
• •
2 r θ : Aceleración complementaria.
•
r Sen(θ)Cos(θ) λ 2: Proyección sobre la tangente al meridiano de la aceleración centrípeta,
debida a la rotación sobre el paralelo.

Componentes en eλ:
••
r Sen(θ) λ : Aceleración tangencial debido a rotación sobre el paralelo.
• •
2 r Sen(θ) λ : Aceleración complementaria, proyección sobre el plano paralelo de la
•
velocidad radial ( r Sen(θ)).
• •
2 r Cos(θ) θ λ : Aceleración complementaria en paralelo debido a la proyección sobre
•
el plano paralelo de la velocidad tangencial en el meridiano (r Cos(θ) θ ).

••
r Sen(θ) λ

r Sen(θ)
λ
θ
••
rSen(θ) λ

En coordenadas intrínsecas

En las coordenadas intrínsecas o de trayectoria, la descripción del movimiento de la

partícula se realiza en términos de componentes que son tangentes o perpendiculares a la
trayectoria. Esto resulta particularmente útil cuando se conoce la curva que sigue la
partícula, por ejemplo un automóvil que recorre un camino sinuoso.
El movimiento de una partícula sobre una curva en el espacio puede describirse en
función de la coordenada generalizada q1=s, la longitud de arco medida desde un punto de
referencia.
El vector de posición se expresa en función de esta coordenada generalizada:
r(s) = x(s) i + y(s) j + z(s) k

26
 En esta expresión, “s” es depende del tiempo (a medida que transcurre el tiempo, el
móvil se desplaza sobre su trayectoria, y la longitud del arco recorrido aumenta).
Para determinar la velocidad se aplica la regla de la cadena:
dr(s) dr(s) ds dr(s) • dx(s) dy(s) dz(s)
s • • •
dt ds d t ds ds s ds s ds s
v(t) = = = = i+ j+ k=
dx(s) dy(s) dz(s)
•
=s ( ds i + ds j + ds k)
dx(s) d y(s) d z(s)
Siendo ( ds i+j ds + ds k) = T, el versor tangente (T es tangente a

la curva)
El vector velocidad resulta:
•
v(t) = s T .

Esta última ecuación refleja el hecho de que el vector velocidad es tangente a la

•
trayectoria, por lo tanto, tiene solamente componente en el versor tangente T. s es la
rapidez o celeridad.

La aceleración
dv(s) se obtienedderivando
T la velocidad respecto del tiempo.
•• • •
a(t) = dt =sT+ s ds s
dT N
Siendo ds = ρ , donde ρ es el radio de curvatura (no debe ser confundido con

el radio ρ de coordenadas cilíndricas), la aceleración en coordenadas

intrínsecas se puede expresar
2 como:•
dv(s) s
••
a(t) = dt =sT+ ρ N

En coordenadas intrínsecas, aparecen dos componentes de la aceleración: la

••
aceleración tangencial ( s ), responsable del cambio en el valor de la velocidad, y la
•2
s
aceleración normal( ρ ) responsable del cambio de dirección de la velocidad.

En el dibujo de la r
T r r
derecha se muestra un aN aT
r
movimiento curvilíneo N
que tiene lugar en un
plano. Se muestran tres
instantes diferentes del
r
movimiento y en cada V

27
uno de ellos el radio de curvatura y su correspondiente círculo osculador. También se
muestran el versor normal y tangencial en el primer círculo oscilador. En el segundo
círculo oscilador solo se muestra la velocidad y, por último, en el tercer círculo se indican las
aceleraciones normal y tangencial.
Si bien la descripción de la velocidad y aceleración en función de la coordenada
intrínseca “s” es sencilla, se debe tener en cuenta que los versores T y N cambian de
dirección a medida que la partícula evoluciona en su movimiento, además, el radio de
curvatura ρ se mide desde un centro de rotación cuya posición se modifica con-
tínuamente.

Anexo 1

Sistemas Referenciales
1. Sistemas de coordenadas generalizadas

El número de coordenadas necesarias para describir un sistema es igual al número de

grados de libertad que el mismo posee. Así, para determinar la posición de una partícula en
el plano son necesarias dos coordenadas. Si bien por defecto se piensa en distancias
medidas sobre ejes perpendiculares (Sistema de ejes cartesianos), este par de coordenadas
puede tener cualquier dimensionalidad.
Se denomina sistema de coordenadas generalizadas a cualquier sistema de
coordenadas que permita describir un sistema físico.
Un sistema que posee “n” grados de libertad, puede ser expresado por q1, q2, ..., qn
coordenadas generalizadas. Casos particulares son las coordenadas cartesianas, las
coordenadas cilíndricas y las coordenadas esféricas.
Cartesianas Cilíndricas Esféricas
q1 x ρ r
q2 y φ λ
q3 z z θ

Estas coordenadas generalizadas definen un espacio al que se denomina espacio de

representación. En el caso de las coordenadas cartesianas, el espacio de representación es
coincidente con el espacio físico o real.
Para los demás sistemas de coordenadas generalizadas, la relación entre los espacios se
establece mediante relaciones, por ejemplo, para el espacio de tres dimensiones:

x = x(q1, q2,q3)
y = y(q1, q2,q3) [1]
z = z(q1, q2,q3)

Radio vector o Vector de Posición:

La ubicación de un punto en el espacio cartesiano viene dado por el radio vector

r = xi + yj + zk
donde r es una combinación lineal de los versores bases.

Utilizando las relaciones [1], el radio vector en coordenadas generalizadas queda:

r = x(q1, q2,q3)i + y(q1, q2,q3)j + z(q1, q2,q3)k [2]

28
Base de coordenadas generalizadas

Partiendo del vector posición en coordenadas generalizadas, ecuación [2], derivando

este vector posición respecto de las coordenadas generalizadas, se obtiene otro vector con
los cuales se puede formar una base:
∂r ∂r ∂r
b1 = , b2 = , b3 = [3]
∂q1 ∂q 2 ∂q 3
Estos vectores no necesariamente son ortogonales, pero siempre cumplen la condición
de no coplanaridad, lo que asegura que puedan formar una base para un espacio
tridimensional.
Dividiendo por el módulo de cada vector, se obtiene una base de versores:
∂r ∂r ∂r
∂q 1 ∂q 2 ∂q 3
e1 = ; e2= ; e3= [4]
∂r ∂r ∂r
∂q 1 ∂q 2 ∂q 3
z
Coordenadas Cilíndricas

En este caso en el
espacio de configuración se
adoptan como coordenadas
r
y dos distancias y un ángulo:
z q1 = ρ ; q2 = φ ; q3 =
z
φ ρ
La interpretación de estas
x variables se muestra en la
figura de la izquierda.

Las tres coordenadas independientes que definen la posición del punto son la proyección
de r sobre el plano horizontal, el ángulo que forma esta proyección con el eje “x”, y la
altura “z”.

Las relaciones entre coordenadas cilíndricas y cartesianas son:

 x = ρ cos φ

 y = ρ sen φ
z = z
 [5]
ρ = x 2 + y 2

  y
φ = ArcTang  
 x
z = z

Utilizando la base de coordenadas cartesianas, y las relaciones entre las coordenadas
cartesianas y cilíndricas, Los versores base del sistema referencial de coordenadas
29
cilíndricas se definen utilizando las ecuaciones [4]. Para ello, se escribe la ecuación del
vector posición utilizando las relaciones [5]

r ( ρ , φ , z )= ρ cos(φ ) i + ρ sen(φ ) j + z k
∂r
∂ρ
eρ= = cos (φ) i + sen (φ)j [6]
∂r
∂ρ

∂r
∂φ
eφ= = -sen (φ) i + cos (φ)j [7]
∂r
∂φ

∂r
∂k
ek= =k [8]
∂r
∂k

Se verifica que estos versores son perpendiculares entre si, y en la secuencia indicada
definen una base directa.

Las componentes del vector posición en la base de coordenadas cilíndricas se obtienen

realizando los productos escalares:
r ( ρ ,φ , z ) • e = ρ
ρ

r ( ρ , φ , z )• eφ = 0
r ( ρ , φ , z )• e = z
k

El vector posición en la base de coordenadas cilíndricas resulta:

r ( ρ ,φ , z ) = ρ e + z k
ρ [9]

Partiendo de las ecuaciones [6],[7] y [8], puede escribirse la matriz de cambio de base
de coordenadas cartesianas a coordenadas cilíndricas. Para ello, se escriben los versores
base de coordenadas cilíndricas como filas:
 cos(φ ) sen(φ ) 0
 
[M]cart-> cil= − sen(φ ) cos(φ ) 0 [10]
 0 0 1
Aplicando esta matriz a vectores escritos en la base cartesiana, resulta el vector escrito
en la base de coordenadas cilíndricas.
Si se compara esta matriz de cambio de base, se comprueba es igual a la matriz de
rotación alrededor del eje z.

30
Coordenadas Esféricas

En las denominadas coordenadas esféricas, las tres coordenadas generalizadas toman

la forma:
q1 = r ; q2 = λ ; q3 = θ
En forma similar a lo hecho con las coordenadas cilíndricas, caracterizamos
las coordenadas esféricas en base a sus relaciones con un referencial cartesiano:
Se observa que “r” es la distancia desde el punto hasta el origen de coordenadas
(módulo del vector de posición), λ es el ángulo medido desde el eje “x” hasta la
proyección del vector de posición sobre el plano “xy” y θ es el ángulo entre el eje “z” y el
vector de posición.

r
θ
y

Las relaciones que ligan las coordenadas esféricas con las cartesianas son:
 x = r sen(θ ) cos(λ )

 y = r sen(θ ) sen(λ )
 z = r cos(θ )
 [11]

r = x 2 + y 2 + z 2

  x2 + y2 
  
θ = ArcTang  
 z
 

λ = ArcTang  y 
  x
En base a las relaciones [11], el vector de posición tiene la forma:
r[r,λ,θ] = r sen(θ) cos(λ) i + r sen(θ)sen(λ) j + r cos(θ)k [12]
Utilizando en mecanismo de derivación parcial respecto de las coordenadas
generalizadas, se obtiene la terna de versores base de coordenadas esféricas:
∂r
er= ∂r = sen(θ) cos(λ) i + sen(θ)sen (λ)j + cos(θ)k [13]
∂r
∂r
∂r
∂θ
eθ= ∂r = cos(θ)cos (λ) i +cos(θ) sen (λ)j – sen(θ)k [14]
∂θ
31
 ∂r
eλ= ∂λ = -sen(λ)i + cos(λ) j [15]
∂r
∂λ
Definidos los versores de la base de coordenadas esféricas, se determinan las
componentes del vector de posición en esa base:
r ( r, λ ,θ ) • e = r r

r ( r , λ , θ )• eλ = 0
r ( r , λ , θ )• eθ = 0
Quedando entonces, el vector de posición en coordenadas esféricas:
r ( r, λ ,θ ) = r e [16]
r
De la misma manera que en coordenadas cilíndricas, partiendo de las ecuaciones
[13],[14]y [15], se puede formar la matriz de cambio de base de cartesianas a esféricas.

Coordenadas intrínsecas

Sea un movimiento sobre una curva en el espacio. La posición de una partícula sobre
esa curva se establece mediante una coordenada generalizada, la cual puede ser la longitud
medida desde un punto arbitrario de la curva, que se toma como origen.
q1 = s

Partiendo de la ecuación del vector de posición r[t] en coordenadas cartesianas:

r[t] = x[t] i + y[t] j +z[t]k

la longitud de arco (s) se determina mediante la integral:

( x ') + ( y ') + ( z ')
t

s= ∫
2 2 2
dt
0 = s[t]
Siendo esta función siempre creciente, en teoría es factible determinar su inversa:
t = t[s]
Reemplazando en la ecuación del vector de posición en función del tiempo, resulta:
r[s] = x[s] i + y[s] j +z[s]k

Se definen los versores T (Tangente), N (Normal) y B (Binormal):

T = dr[s ]
ds

dT [s ]
ds dT [s ]
dT [s ]
N= ds = k= curvatura de flexión
ds
B=T×N

Estos versores forman el triedro intrínseco, y pueden utilizarse como una base para
expresar los vectores posición, velocidad y aceleración.
El recíproco de la curvatura de flexión es el “radio de curvatura” ρ:
1
ρ=
k

32
Anexo 2
Transformación de Coordenadas

En el Anexo 1 se detallan aspectos de los sistemas referenciales. Conocidos los

componentes de un vector en un sistema de coordenadas, se buscará la manera de conocer
las componentes del mismo en otro sistema.
La relación entre sistemas de coordenadas se dan a través de las matrices de
transformación. Se verá de qué manera de definen estas matrices de transformación entre
sistemas referenciales ortogonales.

Sean dos sistemas ortogonales “xyz”, “x’y’z’” a las cuales se asocian las bases ijk,
i’j’k’.
Los versores de la base i’j’k’ pueden escribirse como combinación lineal de los versores
dados por i j k:
i’ = a11i + a12 j + a13 k
j’ = a21i + a22 j + a 23 k [A2-1]
k’ = a31i+ a32 j + a33 k

donde aij son los cosenos directores de los versores i’,j’,k’ en el sistema ijk.
Dado el vector A, de componentes Ax,Ay,Aj en el sistema “xyz”, sus componentes en
el sistema “x’y’z’” se determinan mediante los productos escalares:
Ax’= A•i’ = A• (a11i + a12 j + a13 k)
Ay’= A•j’ = A• (a21i + a22 j + a 23 k) [A2-2]
Az’= A•k’ = A• (a31i+ a32 j + a33 k)

Utilizando notación matricial, las relaciones listadas en [A-2] resultan:

A’ = [M]A [A2-3]

donde [M] , matriz de transformación o cambio de base:

 a 11 a 12 a 13 
 a 23 
[M] = a 21 a 22 [A2-4]
a 31 a 32 a 33 

En la ecuación [A2-3], los vectores A’ y A se escriben como vectores columnas.

Siendo los sistemas xyz y x’y’z’ ortogonales, la matriz [M] es una matriz ortonormal,
y su matriz inversa es equivalente a la matriz transpuesta. Premultiplicando ambos
miembros de la ecuación [A2-3] por la inversa de [M]:

[M]-1 A’ = [M]-1[M] A = [1] A = A

Reemplazando por la transpuesta, resulta la relación:

A = [M]TA’ [A2-5]

33
Unidad 2: Cinemática Relativa

Este Capítulo también es conocido como "Cinemática del Movimiento Compuesto

de la Partícula" o, también, "Cinemática en Sistemas de Coordenadas Móviles".
Tiene importancia como herramienta para facilitar el análisis de movimientos
complejos, que provienen de la combinación o efecto simultáneo de distintos movimientos
sencillos, o bien conocidos.
Se trabaja con sistemas de referencia simultáneos uno fijo o "absoluto" y al menos
uno móvil o "relativo".

Sistemas de Referencia en Movimiento

Dado un sistema de referencias fijo Oxyz de centro O diremos que el sistema O'x'y'z'
está en movimiento respecto a aquel cuando ocurra alguno de los siguientes casos:
z z´

k´ y´
j´
i´
R
x´
k

j y
i
x
( ( (
dR(t ) di dj dk
≠ 0, ó ≠ 0, ó ≠ 0, o bien, ≠0
dt dt dt dt

Sistemas Rígidos y Deformables

Trabajaremos en nuestro estudio con sistemas de referencia ortogonales y rígidos, es

decir, que conservan la mutua perpendicularidad de sus ejes, durante todo el movimiento
(es decir, para todo t).
(( (( ((
La condición de RIGIDEZ es i . j = j .k = k .i = 0 en tal caso
Los sistemas que no cumplen la misma se dicen DEFORMABLES

34
1) Movimiento con Sistemas en Traslación (únicamente)
z´
z P
r
k´ rM
y´
i´ j´
x´
r r
R (t ) rF

j y
i
x

- El movimiento de P puede describirse tanto desde Oxyz, cómo desde O'x'y'z'.

- Oxyz es un SISTEMA FIJO (F)
- O'x'y'z' es un SISTEMA MÓVIL (M) en TRASLACIÓN.
Un sistema móvil evoluciona en Traslación, respecto de otro, fijo cuando su origen O'
r
se mueve según una ley cualquiera R(t ) pero sus ejes se mantienen paralelos a sí mismos.
(Los ejes de O'x'y'z' se dibujaron paralelos a los Oxyz únicamente como criterio
simplificativo. Pueden tener otras direcciones).

r r r
rF = R + rM (1)

r
Aquí rF : Vector posición referido al Sistema Fijo
r
rM : Vector posición referido al Sistema Móvil
r
R: Vector posición del Origen del Sistema Móvil respecto del Sistema
Fijo

Derivando (1), obtenemos la VELOCIDAD

r r r
r dr dR dr
vF = = + (2)
dt F dt dt M

En esta ecuación:
dR
: velocidad de arrastre debida al movimiento del centro O'
dt

35
 r
dr r
= v M : velocidad relativa
dt M
r
dr r
= v F : velocidad absoluta
dt F

Derivando (2), obtenemos la aceleración

r r r
dv F d 2 rF d 2 R d 2 rM
aF = = = 2 + (3)
dt dt dt dt 2
donde:
r
a F : aceleración absoluta
r
r d 2 rM
aF = : aceleración relativa
dt 2
r
d 2R
: aceleración de arrastre, por movimiento del origen del Sistema Móvil o,
dt 2
respecto del Sistema Fijo.

2) Movimiento con Sistemas en Rotación (únicamente)

En la figura se muestran dos sistemas, el Oxyz que consideramos como fijo, y el
O´x´y´z´, animados de un movimiento de rotación con velocidad angular ω. Los
orígenes O y O´ permanecen siempre coincidentes. La posición de la partícula P se
describe mediante la del vector de posición R(t), el cual puede referirse a ambos
sistemas:

z
z´

R(t)
ω

k´
k y´
j´
y
i j
i´
x

x´

36
Vector de posición referido al sistema móvil
r ( ( (
R (t ) M = x' i '+ y ' j '+ z ' k ' (4)
Vector de posición referido al sistema fijo
r ( ( (
R (t ) F = xi + yj + zk
Se reescribe la ecuación (4) indicando que las componentes x´y´z´ dependen del tiempo
r ( ( (
R(t ) = x' (t )i (t )'+ y ' (t ) j ' (t ) + z ' (t )k ' (t ) (4')
r r
Nótese que R(t ) = R (t ) , constituyen el mismo vector, es decir, es válida la igualdad
F M

vectorial, aunque estén expresados en bases diferentes.

r
Evaluamos la derivada de R(t )
r ( ( (
d R (t ) • ( • ( • ( di (t ) dj (t ) dk (t )
= x' (t )i '+ y ' (t ) j '+ z (t )k '+ x' + y' + z'
dt dt dt dt

(5) (6)
Al evaluar las derivadas se ha tomado en cuenta que los versores de la base móvil cambian
en el tiempo y tienen por lo tanto derivada no nula.
( ( (
El término (5), al mantener i ' , j ' , k ' fijos, o sin variar, constituye una velocidad
• • •
r
dR (t )
( x' , y ' , z ) en el sistema móvil O'x'y'z'. Lo llamaremos, por lo tanto M.
dt
Un observador solidario a la base móvil vería a la partícula P animada de esta velocidad.
En cuanto a (6) vale lo siguiente: Por condición de rigidez del sistema O'x'y'z'
( (
( ( d ( ( di ' ( ( dj '
i '. j ' = 0 (i '. j ' ) = . j '+i '. = 0
dt dt dt
( (
( ( d ( ( dj ' ( ( dk '
j '.k ' = 0 ⇒ ( j '.k ' ) = .k '+ j '. =0
dt dt dt
( (
( ( d ( ( dk ' ( ( di '
k '.i ' = 0 (k '.i ' ) = .i '+ k '. =0
dt dt dt

( (
di ' ( ( dj '
Llamaremos wz ' (t ) = . j ' = −i '.
dt dt
( (
dj ' ( ( dk '
wx ' (t ) = .k ' = − j '.
dt dt
( (
dk ' ( ( di '
w y ' (t ) = .i ' = −k '.
dt dt
Con wx ' (t ), w y ' (t ), wz ' (t ) funciones escales de t.

37
 Y por otra parte, recordando que la derivada de un vector de magnitud constante es
ortogonal al vector:
(
dj ' (
. j '= 0
dt
(
dk ' (
.k ' = 0
dt
( ( (
d i ' dj ' dk ' ( ( (
Ahora expresamos , y en la base i ' , j ' , k ' quedando
dt dt dt

( ( ( (
dj '  dj ' (  (  dj ' (  (  dj ' (  (
=  .i ' .i ' + . j ' . j ' + .k ' .k '
dt  dt   dt   dt 
( ( ( (
d i '  di ' (  (  d i ' (  (  d i ' (  (
=  .i ' .i ' + . j ' . j ' + .k ' .k '
dt  dt   dt   dt 
( ( ( (
dk '  dk ' (  (  dk ' (  (  dk ' (  (
= .i ' .i ' + . j ' . j ' + .k ' .k '
dt  dt   dt   dt 
Sabiendo que los siguientes productos escalares son iguales a cero
( ( (
 di ' (  (  dj ' (  (  dk ' (  (
 .i ' .i ' = 0  . j ' . j ' = 0  .k ' .k ' = 0
 dt   dt   dt 

Volviendo a la ecuación (6) y definiendo el vector ϖ = (ω x ; ω y ; ω z ) se tendrá

( ( (
di (t ) d j (t ) dk (t )
x' + y' + z' =
dt dt dt
( ( ( ( ( (
x' [ωz ′(t ). j ′ − ωy ′(t ).k ′] + y ′[− ωz ′(t ).i′ + ωx′(t )k ′] + z ′[ωy ′(t )i′ − ωx′(t ) j ′] =
(
[ωy ′(t ).z′ − ωz ′(t ). y′]i(′ + [ωz ′(t ).x′ − ωx′(t ).z′](j ′ + [ωx′(t ). y′ − ωy′(t ).x′]k ′ =
i′ j′ k′
ω x ′ ω y ′ ω z ′ = ϖ × RM
x′ y′ z′
Finalmente queda
r r
 dR   dR  r
  =   + ϖ × RM (7)
 dt  F  dt  M
Admitiendo la existencia de un vector ϖ = (ω x′;ω y ′;ω z ′) que denominaremos
VECTOR ROTACIÓN, es posible definir un operador δ aplicable a cualquier vector
expresado en sendos sistemas de referencia fijo y móvil.

38
  d ( )  d ( )
δ →  =  + (ϖ × ( ) )M
 dt  F  dt  M
Analizando la VELOCIDAD y ACELERACIÓN en un sistema que rota respecto de
otro fijo se tendrá:
z
z´

rF = rM
ω

k´ k y´
j´
y
i
j
i´
x

En el sistema móvil O x´y´z´es

r ( ( r
rM = x′.i ′ + y ′. j ′ + z ′.k ′
y como vectorialmente
r r
rM = rF
r ( ( r
rF = x′.i ′ + y ′. j ′ + z ′.k ′
Aplicando el oprerador δ a esta última expresión se tendrá
r r
 d (r )   d (r )  r
  =  + (ϖ × (r ) )M
 dt  F  dt  M
Resultando en definitiva la expresión de la velocidad
(vr )F = (vr )M + (ϖ × ( rr ) )M
donde la velocidad referida al sistema fijo o ABSOLUTO es igual a la velocidad del
sistema móvil o RELATIVO más la velocidad de arrastre por efecto de la rotación del
sistema móvil dado por ϖ = (ω x′;ω y ′;ω z ′)
La aceleración se obtiene aplicando el operador δ a la velocidad absoluta quedando

r r r
 d (v )   d ((v )M + (ϖ × (r ) )M )  r r
  =  + (ϖ × ((v )M + (ϖ × (r ) )M ) )M
 dt  F  dt M
39
 r r r
 d (v )   d (v )   d (ϖ × (r ) )M   d r  r
  =  +  + ϖ × (r ) + ϖ × (ϖ × rM )
 dt  F  dt  M  dt M  dt  M

Resultando en definitiva

r r
 d (v )   d (v )  • r
 d r  r
  =  + ϖ × rM + 2.ϖ × (r ) + ϖ × (ϖ × rM )
 dt  F  dt  M  dt  M
Donde la aceleración absoluta referida al sistema fijo es igual a la aceleración relativa
al sistema móvil más la aceleración de arrastre debida a la variación de ϖ (aceleración de
rotación) más la aceleración de arrastre complementaria más la aceleración de arrastre
centrípeta.

3) Movimiento con Sistemas en Traslación y Rotación Simultánea

De la ecuación (1)
r r r
rF = R + rM
r
r r r dR
(v )F = (v )M + (ϖ × (r ) )M +
dt
Los dos últimos términos del segundo término de la ecuación anterior representan a la
velocidad de ARRASTRE de TRANSPORTE. Posee un término debido a la ROTACIÓN
r
DEL SISTEMA MÓVIL (ϖ × (r ) )My otro debido al movimiento del ORIGEN MÓVIL
r
dR
O´ respecto del sistema fijo .
dt
Para obtener la aceleración se debe derivar la anterior respecto del tiempo obteniendo

r
d 2 (R ) • r
r  d 2 (rr ) 
 d (v )  r r
(a )F =  =  2 
 = + ϖ × rM + ϖ × (ϖ × rM ) + 2.(ϖ × v )M
 dt  F  dt  F
2
dt
Donde los términos son conocidos.
r
d 2 (R )
Se debe hacer notar que el término se incorpora a la descripción del
dt 2
movimiento como ARRASTRE o TRANSPORTE originado en el movimiento del origen
del sistema móvil respecto del sistema fijo. El resto de los términos son las aceleraciones
de arrastre inducidas por la rotación del sistema móvil.
IMPORTANTE: En la expresiones (1) ; (10) y (11) las igualdades vectoriales
imponen, a los efectos de la resolución de problemas prácticos, conocer necesariamente la
transformación lineal [A] que permita relacionar las bases (i,j,k) con (i´,j´,k´) de los
sistemas fijo y móvil.
40
Unidad 3: Cinemática del Cuerpo Rígido
Introducción

En el estudio de la cinemática del sólido se toman en cuenta las dimensiones de los

cuerpos. Se verá que el movimiento general de un cuerpo siempre se puede describir como
la combinación de movimientos sencillos de traslación y rotación. Se analizarán casos
particulares del movimiento de sólidos, tales como el movimiento plano, el movimiento
polar o con un punto fijo.

Definición

Cuerpo rígido: se define como cuerpo sólido rígido aquel elemento constituido por
partículas que mantienen sus distancias relativas constantes.
Cuerpo deformable: es aquel en el cual las distancias relativas entre partículas pueden
variar.
Estudiaremos la cinemática del sólido rígido, o sea las posibilidades y movimientos
que puede presentar en el espacio.
El concepto de sólido rígido es una idealización; en la realidad física, el mismo no
existe. Nos encontramos en ella con sólidos reales: deformables.
En general, las deformaciones de los sólidos reales son pequeñas comparadas con sus
dimensiones y podremos aplicar los conceptos y definiciones relativas al sólido rígido
ideal.
En los casos en que las deformaciones sean importantes respecto de la geometría del
sólido, al mismo no podrán aplicarse los conceptos referidos al sólido rígido.

Condiciones Cinemáticas del Sólido Rígido

Partiendo de la definición dada anteriormente, estableceremos relaciones que ligan las

distintas variables cinemática del sólido rígido.

Condición analítica de rigidez

La condición de rigidez la podremos expresar por:
c = cte

P1 c = cte distancia entre dos puntos

C P1 y P2 dos puntos del rígido

P2

Referenciando al sólido rígido respeto de un punto de referencia "O" quedará:

41
 P2 − P1 = cte
( P2 − P1 ) es el vector diferencia entre los vectores posición de los puntos P y P . Su
1 2
módulo es la distancia entre los mismos.

Condición cinemática de velocidades

P1
r
( r1
k
O r
c
( r
( j r2
i
P2

( P2 − P1 ) = cte
( P2 − P1 ).( P2 − P1 ) = cte = ( P2 − P1 ) 2
Si el sólido está en movimiento las posiciones de los puntos P1 y P2 son funciones del
tiempo entonces:
d  dP dP 
( P2 − P1 ) 2 = 0 2.( P2 − P1 ). 2 − 1  = 0
dt ⇒  dt dt 

dP2
v2 =
dt
dP
v1 = 1
2( P2 − P1 ).(v2 − v1 ) = 0 donde dt

r r
( P2 − P1 ).v2 = ( P2 − P1 ).v1 (1)

Dividiendo por el módulo de

( P2 − P1 )

r r
( P2 − P1 ) r ( P2 − P1 ) r e(21 .v2 = e(21.v1
.v2 = .v1
P2 − P1 P2 − P1
queda

42
 Denominada condición cinemática para las velocidades. En definitiva la proyección de
las velocidades de dos puntos sobre la dirección de la línea que une dichos puntos es
constante.
r
e(21 .v1
P1
r
v2
r
v1
r
P2 e(21 .v2

Condición cinemática de aceleraciones

Partiendo de la condición cinemática de velocidades y realizando las derivadas se tiene

r r
 dP 2 dP 2  r dv 2 dv1  dP 2 dP 2  r
 − .v2 + ( P2 − P1 ). = ( P2 − P1 ). + − .v1
 dt dt  dt dt  dt dt 
r r r r r r
v21 .v2 + ( P2 − P1 ).a2 = v21 .v1 + ( P2 − P1 ).a1
r r r r r
v21 .(v2 − v1 ) + ( P2 − P1 ).a2 = ( P2 − P1 ).a1

r 2
v21 r r
+ e(21.a2 = e(21.a1 r
P2 − P1 v21 : velocidad relativa
r
a2 : aceleración de P
2
r
a1 : aceleración de P
1

r 2
r r v21 condición cinemática de aceleraciones
a 2 .e(21 = a1 .e(21 −
P2 − P1
( P2 − P1 ) es
O sea que la proyección de la aceleración en el punto P2 sobre la dirección
r
a
igual a la proyección de 1 sobre la misma dirección menos el cuadrado escalar de la
velocidad relativa dividido por la distancia que separa los puntos.

Movimientos del Sólido en el Espacio

El sólido rígido no vinculado en el espacio puede realizar movimientos o sea tener
distintas ubicaciones y orientaciones en el espacio a través del tiempo. Estos movimientos
pueden ser puros o compuestos.

43
Movimientos puros del sólido en el espacio

Translación pura

Este tipo de movimiento se da cuando la recta que une dos puntos cualquiera conserva
su módulo y dirección en las distintas posiciones que ocupa el sólido a través del tiempo.

P´1
Trayectoria del punto
P1
P"1
r
c

P1
r
P´2 c

r
c
P"2

P2 Trayectoria del punto

P2

En t = t se tiene P1 , P2
en t = t´ se tiene P´1 , P´2
Para t = t" se tiene P"1 , P"2

Características cinemática del movimiento de traslación del sólido

ƒ Todos los puntos del rígido tienen idénticas velocidades instantáneas:
r lim ( P´1 − P1 ) r lim ( P´2 − P2 )
v1 = v2 =
∆t → 0 ∆t ∆t → 0 ∆t (1)
P2 = P1 + c P´2 = P´1 +c (2)
Reemplazando (2) en (1) resulta:
r r
v1 = v2
ƒ Los vectores desplazamiento c y velocidad de traslación v son vectores libres.
ƒ Toda línea del sólido se traslada manteniéndose paralela a si misma.
Supongamos tener las bielas P1 - O1 y P2 - O2 que pivotean en O1 y O2; la placa que
une P1 - P2 pasa a ocupar la posición P´1 - P´2 después de la traslación el vector c se
traslada paralelo a si mismo. La velocidad de todos los puntos es idéntica.

44
 P1 P2
C

P´1 P´2

O1 O2

ƒ Las trayectorias de todos los puntos son idénticas y superponibles.

Rotación pura

Definición: cuando un sólido se mueve de tal manera que permanecen fijos dos de sus
puntos se dice que está en rotación.
En realidad permanecen fijos todos los puntos que se encuentran sobre la recta
definida por esos dos puntos.
Si O1 y O2 permanecen fijos y O3 se encuentra alineado con ellos se tendrá:
Vp
O3 − O1 = α .O2 − O1
derivando la anterior quedará P
(V3 − V1 ) = α .(V2 − V1 ) O3
Pero
V2 = V1 = 0 O2
Entonces se deduce que
V3 = 0
Como O3 es genérico todos los puntos alineados O1
O − O2 permanecerán fijos.
en 1

Velocidad en la rotación
El estado de velocidad del sólido se estudia partiendo de la condición de rigidez. Para
un punto cualquiera P fuera del eje:
( P − O1 ) 2 = cte
( P − O2 ) 2 = cte
Considerando que la posición del punto es función del tiempo se deriva quedando:
2.( P − O1 ).(V p − V1 ) = 2.( P − O1 ).V p = 0 pues V1 = 0
2.( P − O2 ).(V p − V2 ) = 2.( P − O2 ).V p = 0
de lo anterior surge que vp es perpendicular a c y a
( P − O2 ) , luego es perpendicular al

plano
P, O O
1 2 .

La posición del sólido queda determinado por tres puntos. Luego el movimiento del
sólido queda determinado por el movimiento de un plano del mismo.

45
 OO
La recta 1 2 permanece fija en ese plano, luego el movimiento del sólido está dado
por el movimiento del punto. La velocidad del punto P1 y en general la de todos los puntos,
es perpendicular al eje.

En un tiempo ∆t, el plano barre un ángulo ∆ϕ. ∆ϕ

P1 P2
Se define la velocidad angular media como
∆ϕ
ωm =
∆t O2
y la velocidad angular instantánea como
∆ϕ
ω = lim
∆t →0 ∆t
Un punto cualquiera del sólido por ejemplo Pi
que esta a una distancia ri del eje en el
O1
intervalo ∆t describe un arco AB de trayectoria
circular de radio ri.
AB = ∆l = ri .∆ϕ
La velocidad del punto Pi ri
B
dl ∆ϕ dϕ A
Vi = = ri .lim = r. = ri .ω Vi
dt ∆t →0 ∆t dt
O2
La velocidad de un punto es igual al
producto de la velocidad angular por el
radio del punto respecto del eje de rotación.

O1
Vector velocidad angular
Es conveniente otorgarle al vector ω (velocidad angular) carácter vectorial.
1). El módulo es:
dϕ
ϖ =
dt
2) Dirección: la del eje de rotación
3) Sentido: según la regla de la mono derecha
4) Punto de aplicación: cualquier punto del eje de rotación (vector deslizantes)
Vpi = ϖ .(Pi − O1 ).sen α
Nos queda
Vpi = ϖ × (Pi − O1 )
O sea que la velocidad de un punto es el momento del vector rotación respecto de ese
punto (por eso se le otorga carácter vectorial).
El movimiento se denomina uniforme cuando:
dϕ
ω= = cte
dt
46
 ϖ(t)

P1 (x1 , y1 , z1 )

Así definido el movimiento polar, el sólido posee 3 (tres) grados de libertad ya que le
hemos restringido el movimiento a un punto "O" (tres grados de libertad).

Orientación del sólido en el espacio en un movimiento polar

Para orientar en el espacio un sólido que posee un punto fijo; adoptaremos dos
sistemas referenciales: uno fijo (absoluto); y un sistema solidario con el sólido (en
movimiento); ambos ortogonales. Estos sistemas estarán relacionados mediante una
transformación de rotación expresada por una matriz [ A] , ortogonal.
Esta rotación puede descomponerse en tres rotaciones independientes de tal manera
que el producto de las matrices que representan esas rotaciones de como resultado la
matriz:
[A] = [A3 ][. A2 ][. A1]
Ángulos de Euler
Definiremos tres rotaciones independientes que compuestas darán como resultado una
rotación total que nos dará la relación entre un sistema fijo y un sistema solidario con el
sólido (que fija la orientación del sólido).

1º rotación: Ángulo de precesión

Supongamos un sistema inicial (i,j,k,) y realizamos una rotación sobre el eje k. El eje
k' será el mismo k original, los versores i, j rotarán sobre un ángulo de precesión ϕ.
La matriz de transformación es:
 cosϕ sen ϕ 0
[Aϕ ] = − sen ϕ cosϕ 0

 0 0 1
z

j´ j
i´
k =k´
ϕ
y
i j ϕ

x i

49
 La rotación se provoca en el plano x-y. El eje de rotación es el eje z.

2º rotación: Ángulo de nutación

El sistema rota sobre el eje i' (línea de nodo) un ángulo de nutación θ. El sistema pasa
de i',j',k' a i'',j'',k''. El eje i'≡i'' o sea se conserva la dirección i'.
La matriz de transformación es:
k´´ z
1 0 0 
[Aθ ] = 0 cosϕ senθ 
 j´´
0 − senθ cosθ 
θ
θ
j´
ψ
j
ψ
i = i´

3º rotación: Ángulo de rotación propia

El sistema rota sobre el eje k'',la rotación se denomina rotación propia del ángulo ψ.
El sistema pasa de i'',j'',k'' a i1, j1,k1.
La matriz de rotación es: k´´ =
k´´
 cosψ senψ 0
[Aψ ] = − senψ cosψ 0

ψ
θ
 0 0 1  θ
ψ
j

La transformación total resulta: θ ψ

[A] = [Aψ ][. Aθ ][. Aϕ ]
Donde [ A] es una matriz ortogonal, o sea que su determinante es 1 y su inversa es igual
a su transpuesta.
det[ A] = 1
[A]−1 = [A]T
δ ij = 1
i=j
aik .aij = δ ij
δ ij = 0
i≠j

Los tres ángulos ϕ, θ y ψ son en general funciones del tiempo.

ϕ = ϕ(t)
50
θ = θ(t)
ψ = ψ(t)

Movimiento rototraslatorio

Teorema de Euler - Chassles: ϕ - Rotación

El movimiento más general
de un sólido es el de una
traslación seguida de una
C - Traslación
rotación.
Siempre es posible descom-
poner el movimiento de un sóli-
do en una traslación de un punto
más una rotación respecto de un
eje que pasa por dicho punto más
una rotación respecto de un eje
que pasa por dicho punto. V2 ϖ1
Sistema de traslación y rotaciones
Supongamos tener un rígido,
cuyo movimiento resulta de la super- A1
ϖ2 V1
posición de varias traslaciones y A2
rotaciones; podemos considerar que
el estado de movimiento del mismo Ai
queda descripto en cada instante ϖi
por un sistema de vectores.

ϖi = vectores velocidad de rotación

en un eje que pasa por Ai Vi
Vi = vectores velocidad de traslación.

Polo de reducción.
Según el teorema de Euler-Chasles es posible describir el movimiento general anterior,
mediante la traslación de un punto
"O" y una rotación respecto de un V2 ϖ1
eje que pasa por "O".
Al punto "O" que elegiremos
arbitrariamente le denominamos A1
polo de reducción del sistema P V1
ϖ2
rototraslatorio. A2
Si consideramos la velocidad de
un punto cualquiera P. El mismo
estará sujeto a todas las velocidades Ai ϖi
de traslación y rotación. Entonces: O

Vi

51
VP = ∑Vi + ∑ϖ × (P − Ai )
VP = VO + ϖ × (P − O )
Como
(P − Ai ) = (O − Ai ) + (P − O )
Reemplazando queda
VP = ∑Vi + ∑ϖ i × [(O − Ai ) + (P − O )]
VP = ∑Vi + ∑ϖ i × (O − Ai ) + ∑ϖ i ×(P − O )
VP = VO + ϖ × (P − O )
donde ϖ = ∑ϖ i VP2
Es decir hemos expresado la velocidad VP1
de un punto de un sólido en función de la P1
P2
velocidad de un polo "O" y una rotación
respecto del mismo.

Para dos puntos cualquiera P1 y P2 V0

VP1 = VO + ϖ × P1 − O
(
VP 2 = VO + ϖ × (P2 − O )
) ϖ0

Proyectando sobre las recta que une ambos puntos

e1−2 ⋅ VP1 = e1−2 ⋅ VO + e1−2 ⋅ [ϖ × (P1 − O )]
e1−2 ⋅ VP1 = e1−2 ⋅ VO + e1−2 ⋅ [ϖ × (P2 − O ) + (P1 −P2 )]
e1−2 ⋅ VP1 = e1−2 ⋅ VO + e1−2 ⋅ ϖ × (P2 − O ) + e1−2 ⋅ ϖ × (P1 −P2 )
e1−2 ⋅ VP1 = e1−2 ⋅ VO + e1−2 ⋅ ϖ × (P2 − O ) + 0
e1−2 ⋅ VP1 = e1−2 ⋅ VP 2
Cumple con la condición cinemática para
las velocidades.
O sea que nuestro sistema queda reducido a:

ϖ = ∑ϖ i V0
V0 = ∑Vi + ∑ϖ i × ( Ai − O ) ϖ0

Invariante vectorial
Si elegimos un polo arbitrario O´, allí también se tendrá que
ϖ ′ = ∑ϖ i . Luego
ϖ = ϖ ′ es un invariante vectorial para los distintos polos elegidos.

Invariante escalar
Suponiendo haber reducido el sistema al polo "O" y conociendo ϖ y V0.
52
 En O1 se conoce ϖ y V0. Respecto de "O" se tendrá:
VP = VO + ϖ × (P − O )
Tomando "O1" como polo de reducción ϖ V0´
O1
se tendrá:
VP = VO '+ϖ × (P − O1 )
Sabiendo que Vp = VP´ y además
(P − O ) = (P − O1 ) + (O1 − O )
V0 P
igualando:
VO '+ϖ × (P − O1 ) = VO + ϖ × (P − O ) ϖ
O
VO '+ϖ × (P − O1 ) = VO + ϖ × (P − O1 ) + ϖ × (O1 − O )

′
V0 = VO + ϖ × (O1 −O )

Trasladamos el vector traslación; para ello se debe adicionar el momento del vector
rotación por la distancia que separa los dos polos.
r
Proyectando sobre ω .
′ ϖ ϖ ϖ
V0 . = VO . + [ϖ × (O1 −O )].
ϖ ϖ ϖ
Como el segundo término del segundo miembro es cero queda

′ ϖ ϖ
V0 . = VO .
ϖ ϖ
= INVARIANTE ESCALAR

La proyección de las velocidades de traslación sobre la dirección del vector rotación es

constante para cualquier polo elegido, en el mismo instante de tiempos: a esta proyección
se lo denomina invariante escalar.
Los invariantes vectorial y escalar, varían en el tiempo.
Los invariantes I.V. e I.E. definen estados instantáneos del movimiento.

ω .P Vo ω
Vo ω´´ Vo´´
O ´ O´
O ´
´
´
[t]
[t´´]
[t´]

r r r r
En cada punto, el estado de movimiento queda definido por
( v o , w) ,
( v o ' , w' ) ,
r r
(vo ' ' , w' ' ) y el espacio recorrido por el punto P en dt.
Vp.dt = Vo.dt + [ϖ × (P − O )].dt

53
 r r
La proyección de
vo sobre w es constante para los distintos polos, para un mismo
instante de tiempo.

Descomposición

Llamamos descomponer el movimiento de un sólido a elegir un polo de reducción y

r vr
establecer los invariantes; o sea definir el movimiento en términos de w y o.
Podemos elegir infinitos polos de reducción.
Descomposición impropia
r r
v
Se da cuando o no se mantiene paralelo a w en los distintos instantes de tiempo para
el polo "O" elegido.
Descomposición propia
r r
Se da cuando para el polo elegido,
vo se mantiene paralelo a w en los distintos
instantes de tiempo.

Ejemplo: supongamos un cilindro que rueda sobre un plano que posee una traslación
perpendicular la rotación del cilindro.
Descomposición propia: se logra tomando como polo "O" la generatriz de contacto en el
instante (t). En el instante siguiente la generatriz de contacto mantiene la condición de
r r
paralelismo entre
vo y w .
V

z
α
Vo //
y ω
Vo
x
O

ω
Descomposición impropia
r r
Si elegimos el polo "O" en otro punto
vo ' no será paralelo a w y tenemos
descomposición impropia:
V

z
α V
Vo // ω
y O Vo
´
x ω x (0´-
O 0)

Ejemplos de composición de rotaciones

Analizaremos los siguientes casos particulares.

54
Cuerpo rígido sometido a dos rotaciones iguales y de sentido contrario
La velocidad de un punto genérico
Vp = ϖ 1 × (P − O2 ) − ϖ 1 × (P − O1 )
Vp = ϖ 1 × (O1 − O2 ) = cte ∀ P
(O1 − O2 ) no depende del punto. P

El movimiento del sólido es una traslación ϖ1

ϖ2 O1
O2
ϖ1
ϖ1

ϖ1
Dos rotaciones concurrentes
La velocidad del punto P es:
Vp = ϖ 1 × (P − O ) + ϖ 2 × (P − O ) •P

Vp = (ϖ 1 + ϖ 2 ) × (P − O ) ϖ
Vp = ϖ R × (P − O )
r r
Donde
w 1 w
y 2 no están fijos en el
espacio. Se trata de un movimiento
polar. Se puede describir el estado ϖR
instantáneo del movimiento mediante
un único vector ωR aplicado en el ϖ
punto fijo "O".
ϖ1

Eje central del movimiento

r r
vo y el invariante vectorial w
Para un polo "O" arbitrario elegido tendremos un .
r vr
Buscaremos ahora el lugar geométrico donde se cumple que sea w // o o sea que de
una traslación paralela a la rotación.
Vc = α .ϖ
ϖ
r Vc = V 0 .
v
El módulo de c será
ϖ , invariante escalar.
Podemos determinar el punto "C" de su recta de acción multiplicando vectorialmente
r
r ( )
por w a la expresión V c = V o + ϖ × C − O y considerando que c se supone paralelo
v
r
a w o sea su producto escalar cero.

55
Vc × ϖ = V o × ϖ + ϖ × (C − O ) × ϖ = 0
desarrollando ϖ × (C − O ) × ϖ = 0
(C − O ).(ϖ .ϖ ) − ϖ .[(C − O ).ϖ ] = 0 Vc
En particular tomando un punto tal que
r
(C − O ) ⊥ w ϖ
nos queda:
V o × ϖ + (C − O ).ω 2 = 0 C ϖ

(C − O ) = V o ×2ϖ O
ω Vo
nos da la posición del punto "C" del eje central respecto del
polo arbitrario "O".
r
w
w
La recta definida por el punto "C" y la dirección de se denomina eje central del
movimiento y es el lugar de los puntos que tomados como polos de reducción dan lugar a
una traslación paralela a la rotación.
Los puntos del eje central son de velocidad mínima y el movimiento se denomina
helicoidal tangente.
ϖ ϖ
V o. V c.
ϖ ϖ
Vc

V0

Ecuación del eje central

r
Sea (xc yc zc) las coordenadas de "C" referido a un polo arbitrario "O" y w el
invariante vectorial; la ecuación del eje central será:
P − C = α .ϖ Ecuación vectorial

P(xyz)  x − xC = α .ω X

ω  y − yC = α .ω Y
 z − z = α .ω
 C Z Ecuación paramétrica

C x − xC y − yC z − zC
= =
ωX ωY ωZ Ecuación del eje central
r r
Si conocemos en el polo arbitrario "O" a
vo y w . Tendremos que
r r r ωr
Vc = Vmin = Vo . = Iescalar
ω .
El eje central puede o no estar dentro del sólido dado.

56
Movimiento helicoidal tangente

Observamos que un movimiento rototraslatorio se puede reducir a un movimiento de

r r
rotación w v
alrededor de un eje y una traslación o paralela a dicho eje, de magnitud igual
al invariante escalar. Dicho eje se denomina eje central y el movimiento helicoidal.

ϖ
•
C
VC

V
r
Si w se mantiene constante en el tiempo cada punto describe arcos de hélices
r
cilíndricas y su velocidad proviene de la descomposición de una rotación w y una
traslación con la misma dirección.

Axoides

El eje central va tomando distintas posiciones en el espacio, y describe una superficie

reglada que se llama axoide.
De acuerdo a cual sea la terna de referencia que se tome se tendrá:
Axoide fijo: superficie generada por el eje central visto desde una terna fija.
Axoide móvil: superficie generada por el eje central, vista desde un sistema solidario al
sólido considerado.
Ejemplos
Movimiento paralelo

Cilindro que rueda sin resbalar sobre un plano:

z
y

z x

y
ϖ
x 1 2 3

- Visto desde el sistema fijo el eje central se ubica siempre en el plano x-y. El axoide fijo
es un plano.
- Visto desde el sistema de movimiento el eje central se ubica en un cilindro.
- El movimiento en este caso, puede ser descripto como el rodar de un cilindro sobre el
plano.

57
Movimiento polar
Dos conos con generatriz de contacto
El sólido rota respecto z
r
ω
de z con rotación 2 y y
respecto de z con rotación z y x
r
ω1 manteniendo fijo el x
polo "O".
El eje central coincide en dirección con la
composición de las velocidades angulares. Los Eje ϖf
axoides que se generan son: central
- Visto desde un sistema fijo un cono con eje
z. ϖ1
- Visto desde un sistema móvil, solidario con ϖ2
el sólido, el axoide es un cono con eje z.
- Los dos axoides son conos circulares que
ruedan sin resbalar. Su generatriz de contacto O
es el eje central.
Como "O" permanece fijo, en este caso
r
vc = 0

Movimiento Plano

Introducción
Cuando un sistema se mueve de tal manera que todos sus puntos tienen velocidades
paralelas a un plano fijo, el sistema tiene un movimiento "plano". Es suficiente considerar
solo el movimiento de una sección plana, paralela al plano fijo. Es decir, se reduce al
estudio de la cinemática de una figura plana.
Desde otro punto de vista se puede considerar el movimiento como el rodar de dos
axoides cilíndricos con invariante escalar nulo.
La intersección del axoide fijo con un plano dará una curva denominada base, y la del
axoide móvil, rodante. Luego el movimiento plano puede ser considerado como el rodar
de la "ruleta" sobre la "base"
Eje central

Axoide
móvil ϖ
Centro
instantáneo
de rotación
Ruleta

C Base

58
Centro instantáneo de rotación
Ri Vi
El movimiento más general de un sólido es el
rototraslatorio que se reduce a una rotación ω aplicada en un
punto O y una traslación vo. Si el movimiento es plano;
r
ωr debe ser perpendicular al plano π o y vo paralelo al
mismo. Entonces un punto cualquiera Ri, tendrá velocidad: ϖ
• Vo
V i = V 0 + ϖ × (Ri − O ) (Ec. 1)

Siempre existirá un punto "C" donde:

V c = V 0 + ϖ × (C − O ) = 0

O sea un punto de velocidad nula.

Como "C" tiene velocidad nula.

Tomando como polo de reducción "C":

V i = V 0 + ϖ × (Ri − C )

Igualando (1) y (2)

V0 + ϖ × (Ri − O ) = ϖ × (Ri − C )
V0 + ϖ × [(Ri − O ) − (Ri − C )] = 0
V0 + ϖ × (C − O ) = 0

multiplicando por ωr

ϖ × V0 + ϖ × ϖ × (Ri − O ) = ϖ × ϖ × (Ri − C ) = 0
ϖ × V0 = −ϖ × ϖ × (C − O )

(C − O ) = ϖ ×2V0 (3)
ϖ

(C − O ) da la ubicación del centro de rotación "C".

Las coordenadas respecto del polo "O"
(
ωr = ωr .k
( (
(C − O) = X c .i + Yc . j
( (
Vo = Vo . cosα .i + Vo .sen α . j
resolviendo (3)

59
 Vo
Xc = .sen α
ω
Vo
Yc = . cosα
ω
coordenadas del centro instantáneo de rotación

Propiedades del Centro instantáneo de rotación

El CIR tiene velocidad nula. El mismo está definido Vi

para un instante determinado.
En otro instante la ubicación del CIR cambia:
a) En cada instante las normales de las
trayectorias que describen los puntos de
C.I.R.
la trayectoria de la figura móvil pasan
por el CIR.
b) La velocidad de un punto cualquiera es
normal a la recta que lo une con el CIR. Trayectoria

Trayectorias polares
Supongamos dos planos superpuestos, uno de ellos sostiene la figura móvil, el otro
fijo. Si marcamos sobre los dos planos, las ubicaciones del CIR (o polo de velocidades),
obtendremos dos curvas: una sobre el plano fijo llamado base; y otra sobre el plano móvil
llamada ruleta.
O sea que se puede reproducir el movimiento, mediante el rodar sin resbalar de la
ruleta sobre la "base".

Ejemplo
Consideremos la barra AB que se mueve de tal manera que el punto A se desliza sobre
el eje x y el punto B sobre el eje y.
AB=l
B C

C´

Trazando por A y por B las normales a las trayectorias, en su intersección obtenemos

el CIR. Si trazamos para cada posición obtendremos la curva "base", en este caso una
circunferencia de radio "l"= (AB).
Las distintas posiciones del punto CIR respecto de un sistema solidario con la barra,
describe la curva "ruleta" en este caso es una circunferencia de diámetro "l".

60
Deducción de la ecuación de la base

Y y Q
x

j i
a(t) ϕ
O
b(t)
X

a(t) y b(t) son las coordenadas del origen O del sistema xy solidario a la figura móvil.
a y b son funciones del tiempo
ϕ: parámetro

r Larvelocidad
r
de un punto Q respecto del sistema fijo
VQ = VO + ωx(Q − O)
r
V
Si en particular tomamos el centro instantáneo de rotación "C" c
= 0 por definición
de centro instantáneo de rotación.
r r r
VQ = Vo + ωx(Q _ O) = 0

r dQ ( db (
Vo = .I + .J
dt dt
. dϕ
r .( ϕ =
ω = ω .k dt
( (
(C − O) = ( X c − a ).I + (Yc − b).J
donde Xc y Yc son las coordenadas del punto "C" referido al sistema fijo.
O sea
 I J K
r r  dϕ 
Vc = Vc +  0 =0
dt 
0
( X − a ) (Y − b)
 c c 0 

61
Sus proyecciones son:
da
Yc = b + dt
dϕ
dt
da dϕ db
− .(Yc − b) = 0
dt dt
X c = a − dt
db dϕ dϕ
+ .( X c − a ) = 0
dt dt dt

Pero

da db
dt = da dt = db
dϕ dϕ dϕ dϕ
dt ; dt
Nos queda
da
Yc = b +
dϕ
db
Xc = a −
dϕ (4) Ecuaciones paramétricas de la base

Deducción de la ecuación de la rodante

La coordenadas vectorial de la rodante respecto del sistema móvil es
(C − O )

(C − O) = (C − O' ) − (O − O')
Llamando Xc Yc a las coordenadas de C en el sistema móvil.

Tendremos:
( ( ( ( ( (
X c .i + yc j = X c .I + Yc .J − aI − bJ

Proyectando sobre los ejes móviles para lo cual multiplicamos escalarmente ambos
miembros de la igualdad anterior:
( ( (
i : (cosϕ .I + sen ϕ .J )
( (
j : (− sen ϕ + cosϕ .J )

nos queda:
X c = ( X c − a ). cosϕ + (Yc − b).sen ϕ
Yc = −( X c − a ).sen ϕ + (Yc − b). cosϕ
reemplazando (4)

62
 db da
Xc = − . cos ϕ + .sen α
dϕ dϕ
db da
Yc = .sen ϕ + . cosϕ
dϕ dϕ Ecuaciones paramétricas de la ruleta.

Velocidad de alternación del centro instantáneo de velocidad

Hemos definido el centro instantáneo de velocidad o polo de velocidades, a aquel

punto que en cada instante se cumple.
Vi = ωx( Ri − C ) Es decir, el plano móvil gira en un movimiento de rotación
pura respecto de "C".
Este punto singular está definido para cada instante en un punto del plano fijo.
La posición de este punto varía de instante en instante.
Las distintas posiciones de dicho punto referida al sistema fijo constituyen la base.
Podemos decir que el Centro instantáneo recorre la base.
Simultáneamente, referido al sistema solidario con el plano móvil, el CIR recorre la
ruleta.
La velocidad con que el polo de velocidades recorre la base, se denomina "velocidad
de alternación del polo de velocidades".
Vc : (Vxc ,V yc )
Las componentes de esta velocidad, teniendo presente las ecuaciones de base
da
Xc = b +
dϕ
db
Yc = a −
dϕ
Nos queda:
.
• db d  da 
VX c = X c = +  
dt dt  dϕ 
.
• da d  db 
VYc = Y c = −  
dt dt  dϕ 

Estado de aceleración en el movimiento plano

La velocidad de un punto cualquiera es:
dRi
= ωx( Ri − C )
dt Donde "C" es el centro instantáneo de rotación

Derivando respecto del tiempo

d 2 Ri dω  dRi dc 
= x( Ri − C ) + ωx − 
dt dt  dt dt 

63
o bien
d 2 Ri dω
= x( Ri − C ) + ωx[ωx( Ri − C )] − ωxVc
dt dt
d 2 Ri dω
2
= x( Ri − C ) − ω 2 .( Ri − C ) − ωxVc
dt dt

Donde:
Vc : velocidad de alternación del polo.
d 2 Ri
dt 2 : aceleración total del punto Ri
dω
dt : aceleración angular
C: coordenada del polo de velocidades

Polo de aceleraciones
La expresión de aceleración para un punto cualquiera

d 2 Ri dω
2
= x( Ri − C ) − ω 2 .( Ri − C ) − ωxVc
dt dt (1)

Buscamos un punto de aceleración nula

d C o dω
2

2
= x(Co − C ) − ω 2 .(Co − C ) − ωxVc = 0
dt dt (2)

restando (2) (1)

d 2 Ri dω
2
= x( Ri − Co ) − ω 2 .( Ri − Co )
dt dt
que es la derivada con respecto al tiempo de:

d 2 Ri dω
= x( Ri − Co )
dt 2 dt
respecto de las aceleraciones todo acontece como si la figura móvil tuviese un movimiento
de rotación pura en torno a Co.
Este punto Co se denomina polo de las aceleraciones.

Coordenadas del polo de aceleraciones

Por definición de polo de aceleraciones
d 2 C o dω
2
= x(Co − C ) − ω 2 .(Co − C ) − ωxVc = 0
3 dt dt

64
 dω dω ( ( (
= .k Vc = V X c i + VY c j
dt dt
( ( (
ωr = ω .k (Co − C ) = X o i + Yo j

reemplazando
( ( ( ( ( (
i j k   i j k
 dω  ( (  
 − ω ( X oi + Yo j ) −  0 ω = 0
2
0 0 0
 dt 
V X VX c 0 
X o Yo 0   c

Obtenemos:
( dω
i (−Yo. − ω 2 X o + VYc .ω ) = 0
dt
( dω
j ( X o. − ω 2Yo + V X c .ω ) = 0
dt

dω
=0
cuya solución da los valores de Xo y Yo. En particular si la aceleración dt

VYc
X o=
ω
VX c
Y o= −
ω

65
Unidad 4: Dinámica de la Partícula

Principio de Inercia (1º Ley de Newton)

“Todo cuerpo continúa en su estado de reposo o de movimiento rectilíneo y uniforme,

a menos que sea obligado a cambiar de estado por fuerzas que actúan sobre él”
Es equivalente a expresar que “es posible adoptar por lo menos un sistema de
referencia, para el cual toda partícula aislada continúa en su estado de reposo o de
movimiento rectilíneo y uniforme”.
Recuerde se denomina partícula a todo cuerpo cuyas dimensiones y orientación
resultan despreciables en la descripción de su movimiento.

A los sistemas de referencia en reposo o que se desplazan con movimiento rectilíneo y

uniforme (es decir, no acelerado) se los denomina Sistemas de referencia inerciales.
Una partícula aislada es aquella que no interactúa con otras o sea, no está afectada por
las acciones (“fuerzas”) que estas otras pudieran ejercer sobre ella.

Principio de masa (2º Ley de Newton)

“La derivada de la cantidad de movimiento respecto del tiempo es proporcional a la

fuerza que actúa sobre el cuerpo y tiene su dirección y sentido”.
→ →
Siendo p = m v la “cantidad de movimiento” o “momento lineal”

→ →
→ dp dv →
F= =m =m a
dt dt

Principio de Acción y Reacción (3ª Ley de Newton)

“A toda acción ejercida sobre una partícula se opone una reacción igual y de
sentido contrario que ejerce la propia partícula”.
Se refiere a la interacción entre dos partículas y en otras palabras señala que toda
vez que una partícula ejerce una fuerza sobre otra, ésta ejerce sobre la primera una fuerza
de igual intensidad y dirección pero de sentido contrario.

Principio de independencia de acción de fuerzas

“Todo cuerpo, bajo la acción conjunta de dos fuerzas, describe la diagonal del
paralelogramo en el mismo tiempo que emplearía en describir los lados del mismo, bajo la
acción de cada una de las fuerzas”.

Si dos o más fuerzas actúan sobre una partícula, la aceleración resultante es igual a
la suma de las aceleraciones que adquiriría si cada una de aquellas actuara aisladamente.

67
 Este principio expresa el de superposición de
efectos.

F1

F1 + F2

F2

Movimiento unidimensional de la partícula

En este apartado se estudia el movimiento de una partícula de masa m a lo largo de una

recta. Al quedar de tal modo definida la dirección, es habitual prescindir del tratamiento
vectorial de los distintos elementos que intervienen. El tratamiento vectorial es de todos
modos, general, y es estrictamente indispensable, en el estudio del movimiento curvilíneo,
es decir, en dos o en tres dimensiones.

Según el Principio de Masa:

dp d d  dx  d2x
= (m v) =  m  =m 2 =F =m a
dt dt dt  dt  dt
Expresiones válidas para un partícula cuya masa se acepta como constante.
dp
Este principio, expresado como = F = m a constituye la Expresión diferencial del
dt
Momento Lineal, cuya integración conduce a:
p2 t2 t2
∫ p1
dp = ∫ F dt
t1
⇒ p2 – p1 = ∫ F dt
t1

Que es la expresión integral del Teorema del momento lineal.

t2
∫ t1
F dt se denomina impulso de F, o impulso suministrado por la fuerza F.

d2x
La expresión F= m es la “Ecuación fundamental de la Dinámica”, y es una ecuación
dt 2
diferencial ordinaria cuya solución es x = x(t), y en la que F puede ser función de t (fuerza
variable en el tiempo), x (Fuerza que depende de la posición) o de v (Fuerza que depende
de la velocidad).

dv dv
A partir de F = m , se puede expresar Fv = mv
dt dt
d 1
⇒ Fv = ( mv2)
dt 2
Que es la “Expresión diferencial del Teorema de la Energía”. En esta última expresión
1
T= mv2 se define como Energía Cinética de la partícula.
2
La forma integral del Teorema de la Energía es:
t2 d  1
2
t 2 dT t2
∫t1 dt  2 mv dt = ∫t1 dt dt ⇒ T2 – T1 = ∫t1 F v dt
En la expresión recuadrada, los términos que intervienen representan:
68
 t2
∫ t1
F v dt : Trabajo realizado por la fuerza F entre t1 y t2.
F v : Potencia suministrada por la fuerza F.
1 1 x2
Así, T2 – T1= mv 22 − mv12 = ∫ F dx
2 2 x1

d2x
La Ecuación General de la Dinámica, m = F permite abordar los dos grandes grupos de
dt 2
problemas del estudio del movimiento, a saber:
A) Conocidas las fuerzas o acciones de otras partículas o cuerpos sobre una
partícula dada, obtener la descripción completa del movimiento, identificando para todo
instante, la posición, la velocidad y la aceleración resultantes (Proceso de integración).
B) Conocida la “ley del movimiento”, o sea, x(t), determinar (ademas de la
velocidad y la aceleración) la fuerza resultante o acción que la determina (Proceso de
derivación).

Como se ha dicho, la fuerza F puede ser función de una o más variables, por
•
ejemplo: F(x, x ,t) o de una combinación de ellas, pudiendo inclusive existir relaciones
funcionales que las liguen entre sí.
Ejemplos:
9 La fuerza de rozamiento es función de la velocidad.
9 La fuerza ejercida sobre un cuerpo cargado por un campo electrostático
dependiente del tiempo.
9 La fuerza ejercida por el campo gravitatorio sobre una masa y la fuerza de
recuperación elástica son funciones de la posición.
9 La fuerza de resistencia del aire, en general, es función de la posición y de
la velocidad.

d2x •
Finalmente, la solución de la ecuación m 2
= F(x, x ,t) requiere los conocimientos de
dt
•
condiciones iniciales o de referencia, es decir, conocer los valores que asumen x y x para
un determinado y conocido instante t0.

para t = t 0

Es decir: x = x(t 0 ) = x 0
• •
x = x ( t 0 ) = v 0

Caso de Fuerza dependiente del tiempo, únicamente

d2x • d2x d 2 x dv
La ecuación m 2 =F(x, x ,t) se expresa como m 2 =F(t), y siendo = ,
dt dt dt 2 dt
dv
Será m =F(t) , EDO de primer grado (requiere el conocimiento de una condición
dt
inicial, en este caso, v(0). Esta EDO se puede resolver por separación de variables seguida
de un proceso de integración, o utilizando transformada de Laplace.

69
 Por separación de variables:
dv v t
m =F(t) ⇒ m ∫ dv = ∫ F(t)dt
dt v0 t0

t 1 t
mv – mv0 = ∫ F(t)dt ⇒ v(t) = v0 +   ∫ F(t)dt
t0
m t0

dx t
Conocido v(t), y siendo v(t) = ⇒ dx = v(t)dt ⇒ (x – x0)= ∫ v(t)dt , y reemplazando
dt t0

1 t t
v(t) será x(t) =x0 +v0(t – t0)+   ∫  ∫ F(t)dt dt
 m  t 0  t 0 

Caso de Fuerza dependiente de la velocidad, únicamente

En este caso, es F = F(v), y se supone conocida esta relación funcional.

La ecuación general se expresa como:
dv dv dt
m =F(v), y separando variables: =
dt F(v) m
v dv 1 t 1
⇒ ∫ =   ∫ dt =   (t - t 0 )
v 0 F(v)
 m  t0 m
v dv
Siendo F(v) conocida, ∫ =ϕ(v). En principio, siempre es posible obtener ϕ(v) t de
v 0 F(v)

t - t0
ésta, despejar v. Entonces, resultará v(t) = ψ( v0, ), y conocidos v0, t0,m es v=ψ(t).
m
Conocida la relación funcional entre la velocidad y el tiempo, se puede obtener la ley del
dx t
movimiento. Para esto, se parte de =v ⇒ x – x0 = ∫ ψ(t)dt
dt t0
t
x(t) = x0 + ∫
t0
ψ(t)dt

Caso de Fuerza dependiente de la Posición, únicamente

dv
Para este caso, la ecuación general queda expresada como: m =F(x)
dt
dv
Multiplicando ambos miembros por dx: m dx=F(x)dx ⇒ mv dv = F(x)dx.
dt
Integrando entre las posiciones x0 y x, y para v0 y v:
x v
∫x0
F(x)dx = ∫ mv dv
v0
[I]
donde x0 , v0 son la coordenada y velocidad en un punto de referencia P0
x y v son la coordenada y velocidad en un punto genérico P.
x
∫x0
F(x)dx es el trabajo de la fuerza F(x) al pasar de x0 a x,

Definición de Energía Potencial: Es el trabajo necesario para llevar una partícula desde
un punto genérico hasta el punto que se considera como de referencia.
0 x
Esto se expresa: V(x) = ∫ F(x)dx = - ∫ F(x)dx
x 0

70
Utilizando el potencial, la expresión [I] resulta:
x 1 1
∫x 0 F(x)dx = 2 mv − 2 mv 0 = V(x0) – V(x)
2 2
[ II ]
Si la fuerza que actúa sobre la partícula realiza trabajo positivo, ello representará en una
disminución de la función V(x) (Energía Potencial).
Siendo (1/2)mv2 = T(x), la Energía Cinética, la expresión [ II ] se puede escribir:

T(x) + V(x) = T(x0)+V(x0) = Constante = E [ III ]

donde E es la Energía mecánica total del sistema.

[ III ] expresa la ley (teorema) de la conservación de la energía mecánica total, y significa

que si la fuerza depende solo de la posición, la totalidad del trabajo se produce a expensas
de la energía potencial y da lugar al incremento de la energía cinética.

Un sistema de fuerzas que satisface [ III ] se dice es conservativo, y es condición que la

fuerza dependa solo de la posición.

Cuando aparecen fuerzas que no dependen solo de la posición, aparecen otras formas de
energía, y se dice que intervienen fuerzas “disipativas”, y el sistema será no conservativo.

2
De T(x) + V(x) = E ⇒ (1/2)mv2 + V(x) = E ⇒ v = (E - V(x))
m

dx m x dx t

2 ∫x 0 E - V(x) ∫t 0
y siendo v = ⇒ = dt =(t-t0)
dt
La ley de movimiento podría obtenerse resolviendo la integral y despejando x(t).

x
Teniendo presente que V(x) = - ∫
x0
F dx resulta:
dV(x)
− =F(x)
dx
Es decir que la energía potencial es una función escalar cuya derivada, cambiada de signo,
dV(x)
da la fuerza. De aquí surge que todo punto cuya coordenada cumpla que =0, al
dx
hacer F(x) =0, constituye un punto de equilibrio.

No debe perderse de vista el carácter escalar del potencial, y el carácter vectorial de la

dV(x)
fuerza. La última expresión, es la componente en x (coordenada a la que hemos
dx
restringido el movimiento en este estudio) del Gradiente de V(x). (Se sabe que el
gradiente asigna a un campo escalar un vector cuya dirección es la dirección de máxima
variación, y cuyo módulo es el valor de la máxima derivada).

Curvas de Potencial

Para un movimiento unidimensional, es posible graficar V(x) en función de x.

71
 V(x)
5
E5

3 4
E3

1 6
E1 2
V(x)
E0
0

x5 x3 x1 x0 x2 x4 x6
x

m x dx t
La ecuación ∫ = ∫ dt =(t-t0) expresa que, para una energía dada E, la
2 x 0 E - V(x) t 0
partícula sólo puede evolucionar en un intervalo en que V(x)≤E. El análisis de la curva, de
tal modo, conduce a las siguientes conclusiones:

9 Si la energía inicial es E = E0, la partícula sólo podrá permanecer en x0. Allí estará
en un punto de equilibrio: la tangente geométrica a la curva es horizontal, y ello
dV(x)
indica que = 0.
dx x = x 0
9 Los puntos 4,5 y 6 corresponden a las coordenadas x4,x5 y x6 para las que
igualmente se verifica equilibrio de fuerzas.
9 Si las condiciones iniciales condujeran a un nivel de energía E = E1, el movimiento
de la partícula quedaría confinado entre x1 y x2.
9 Como E = V(x) + T(x) se tendría, si se adopta como nivel de referencia V(x0)=0,
que:
¾ en x0, E1 = T(x0) : La energía es enteramente cinética y además, T(x)
alcanza su valor máximo.
¾ en los puntos 1 y 2, E1 = V(x1) = V(x2): La energía es enteramente
potencial, y esta es máxima para el intervalo [x1,x2].
¾ La velocidad disminuye al acercarse a x1 o a x2, puntos de retorno. El
movimiento cesa y se invierte su sentido al alcanzarlos. El tramo [x1,x2]
suele denominarse “valle de potencial”.

9 Para que la partícula pueda alcanzar la posición x4, la energía inicial debería
alcanzar a E =E2. En x 4, la partícula podrá o bien regresar hacia el valle de x 0,
o escapar hacia el valle de x6.

Estados de Equilibrio

De particular interés resulta el análisis de las situaciones de equilibrio, el carácter

del mismo, y las condiciones de su existencia, para fuerzas actuantes dependientes de la

72
posición. Se ha dicho que los estados de equilibrio se presentan donde las fuerzas aplica-
das sobre la partícula son nulas, o sea F(x) = 0.
dV(x)
Siendo − = F(x), las condiciones necesarias y suficientes para que la
dx
dV(x)
partícula se encuentre en un punto de equilibrio, es que dado V(x); = 0.
dx
Los puntos (1), (2) y (3) cumplen con la
dV(x)
V(x) condición = 0, y se denominan
dx
puntos de equilibrio.
(1) (3)
dV(x)
Desarrollando en series de
dx
potencias:
dV(x)  dV(x)   d 2 V(x) 
(2) =  +   x +
 dx  0  dx  0
2
dx
1  d 3 V(x)  2
+   x +....
x 2  dx 3  0
Despreciando los términos de orden superior y considerando el origen de coordenadas en
dV(x)
el punto donde = 0 queda:
dx
dV(x)  d 2 V(x) 
- = -  2
 x = F
dx  dx 
dV(x)  d 2 V(x) 
En el punto x0 en que se verifica = 0,  2
 puede adoptar un valor >0,<0 o
dx  dx 
igual a 0. Ello conduce a los siguientes casos:
1. Equilibrio Estable:
 d 2 V(x) 
Si  2
 >0, la ecuación general del movimiento es:
 dx  0
d2x  d 2 V(x)  d2x 1  d 2 V(x) 
m 2 = F = -  2
 x ⇒ +   x = 0
dt  dx 0 dt 2 m  dx 2
0
1  d 2 V(x) 
llamando ω 02 =  2
 , se puede rescribir la ecuación como:
m  dx  0
••
x + ω 02 x = 0 cuya solución es del tipo x(t) = A Cos( ω 0 t + θ)

El movimiento así descrito es armónico y ocurre en torno del punto de equilibrio; A y θ

dependen de las condiciones iniciales.
 d 2 V(x)  dV(x)
La condición  2
 >0, en un punto donde = 0, implica que el punto en
 dx  0 dx
cuestión es un mínimo relativo. Puede expresarse entonces que si la función potencial
posee un mínimo en el punto de equilibrio, entonces este es estable. Implica que si una
partícula es alejada del punto de equilibrio, tiende a retornar.

73
 2. Equilibrio Inestable
 d 2 V(x) 
Si  2
 < 0, la ecuación general del movimiento es:
 dx  0
d2x  d 2 V(x)  d2x 1  d 2 V(x) 
m 2 = F =   x ⇒ -   x=0
dt 2
 dx  0 dt 2 m  dx 2  0
1  d 2 V(x) 
llamando ω 02 =  2
 , se puede reescribir la ecuación como:
m  dx  0
••
x − ω 02 x = 0
La solución a esta ecuación diferencial se puede expresar como:
x(t) = C1e ω0 t + C 2 e ω0 t
donde C1 y C2 dependen de las condiciones iniciales. En esta expresión, los coeficientes de
las exponenciales son positivos, lo que se interpreta de la siguiente manera: apartada la
partícula de su posición de equilibrio, la partícula tiende a alejarse indefinidamente.
Las condiciones de equilibrio para un punto de equilibrio inestable son:
 dV
 dx = 0
 2 condición de máximo en la curva de potencial.
d V < 0
 dx 2
3. Equilibrio indiferente.

 d 2 V(x)  d2x  d 2 V(x) 

Si  2

 = 0 ⇒ m = F =   x = 0
 dx  0 dt 2 2
 dx  0

d2x
⇒ m 2 =0 (La fuerza actuante es nula)
dt
Esta última ecuación diferencial tiene por solución x(t) = C1t+C2, donde las
constantes dependen de las condiciones iniciales. Esta ecuación corresponde a un
movimiento a velocidad constante, lo que resulta lógico pues la fuerza que se ejerce
sobre la partícula es nula.
Las condiciones de equilibrio indiferente son entonces:
 dV
 dx = 0
 2 En la curva de potencial, indican potencial constante.
 d V
=0
 dx 2

Movimiento Curvilíneo de la Partícula

Consiste en el estudio del movimiento en dos o tres dimensiones. Es imprescindible el

tratamiento vectorial de los distintos elementos intervinientes.
Según el Principio de masa:
→ →
ma = F [ 1]
Por el principio de independencia, o de Superposición de efectos:

74
 → → → →
F = ∑ Fi , donde F i son las fuerzas actuantes sobre la partícula y F , la Resultante,
→
o fuerza equivalente a la acción combinada del conjunto F i .
Para un tratamiento en coordenadas cartesianas, la ecuación [ 1 ] es equivalente al
tratamiento simultáneo de:
•• •• ••
m x = Fx m y = Fy mz = Fz [2]
→ •• •• ••
puesto que m a = m(x i + y j + z k)= Fx i + Fy j + Fz k

En coordenadas intrínsecas:
→ dv v2
ma =m T+m N = FT T + FNN FB= 0
dt ρ
→ → → → →
La fuerza F es en general, F = F ( r , v , t )
La integración del sistema de 3 ecuaciones diferenciales de 2º orden como [ 2 ],
→
genera 6 constantes arbitrarias, y plantea el conocimiento de condiciones iniciales r 0 y
→
v 0 . La resolución es laboriosa, aunque se siguen los mismos lineamientos desarrollados
para el caso unidimensional. Particularmente si ocurre la independencia dimensional dada
por:
•
Fx = Fx( x ,x,t)
•
Fy = Fy( y ,y,t)
•
Fz = Fz( z ,z,t)
Las tres ecuaciones resultan independientes entre sí y el tratamiento en x, y o z es
similar al ya estudiado.
A continuación, se presenta un resumen de las expresiones vectoriales, válidas para dos
y tres dimensiones, de los teoremas enunciados:

Cantidad de Movimiento o Momento Lineal:

→ →
p=mv
Principio de Masa:
→
d p d  → →
= m v= F
dt dt  
→
Impulso de F (Teorema del Momento Lineal, forma integral):
→ → → → t2 →
p 2 - p1 = m v 2 - m v1 = ∫
t1
Fdt
Teorema de la Energía (forma integral)
2 2
 →   →  b→ →
T2 – T1 = ( /2)m  v 2  -(1/2)m  v1  = ∫ F • dr
1

    a

Teorema de la Energía (Forma diferencial)

dT d  1  → →
=  mv 2  = F • v
dt dt  2 

75
Momento Angular o Momento Cinético

Es el momento del vector Momento lineal de la partícula, respecto de un punto 0.

→ →
Sea p = m v el momento angular de la partícula mv
→
P de masa m y v la velocidad respecto del punto
P
“O” se define
→ →
→
L0
=r × p c
omo el momento angular es entonces el vector
→ →
perpendicular al plano determinado por r y v . En el
caso de un movimiento en el espacio, este vector
cambiará de módulo y dirección a medida que la O
partícula se desplace.
Un caso particular se da cuando la partícula se mueve en un plano y el punto “O” está
contenido en el mismo, en este caso, el vector momento angular será siempre
perpendicular al plano.

Derivando el momento angular respecto del tiempo:

→
 → → → →
d L0 d → → d r dv
=m ( r× v )=m  × v+ r ×
dt dt  dt dt 
 
→
dr →
como × v = 0 por ser paralelos
dt
→ →
d L0 → dv → →
⇒ = m r× =m ( r× F)
dt dt
→
d L0 → →
⇒ = m ( r × F) Forma diferencial del Teorema del momento cinético:
dt
La variación en el tiempo del momento angular respecto de un punto “0” es igual al
momento de la fuerza actuante sobre la partícula respecto de “0”.
→ → t2 → →
La forma integral de este teorema es: L 2 - L1 = ∫ ( r × F)dt
t1

76
Unidad 4: Vínculos o Ligaduras.

Hasta aquí se ha analizado el movimiento libre de la partícula. A la acción de las

fuerzas exteriores, ella responderá sin condicionamientos.
Cuando aparecen restricciones al movimiento, aquella libertad queda condicionada.
Estas condiciones se materializan a través de vínculos o ligaduras.
Se llama vínculo a cualquier sistema material capaz de impedir o condicionar el
libre movimiento de una partícula o cuerpo.
Los vínculos pueden ser simples, dobles, etc, según eliminen uno, dos, etc., grados
de libertad.
La eliminación de grados de libertad constituye restricciones que se imponen al
movimiento. En consecuencia, tres ligaduras simples, convenientemente dispuestas,
eliminan los tres grados de libertad de una partícula, produciendo su inmovilización.
Las restricciones al movimiento suelen presentarse cuando la partícula se halla
condicionada a desplazarse sobre una curva o sobre una superficie. Constituyen casos de
movimientos muy usuales, y de las más habituales aplicaciones prácticas. Representan el
caso del MOVIMIENTO RESTRINGIDO.
Sobre una curva, el movimiento de la partícula se halla restringido en dos grados
de libertad, pudiendo la partícula moverse únicamente según la ley R(s), siendo q = s la
coordenada que determina en cualquier instante la posición.
Sobre una superficie la posición de la partícula queda determinada por R = R(q1,q2),
o sea, función de dos coordenadas. El movimiento tiene dos grados de libertad, en
consecuencia, restringir el movimiento a una superficie representa imponer una condición
de vínculo.
Desde el punto de vista matemático, los vínculos se representan por relaciones de
vínculos.

Ejemplo 1: Una partícula está obligada a deslizarse sobre un alambre cuyo eje describe una
curva plana y = f(x), ∀x / a ≤ x ≤ b. Para P, el alambre constituye un vínculo mecánico,
el que matemáticamente se expresa por la relación de vínculo:
y - f(x)= 0, a ≤ x ≤ b.

Ejemplo 2: Una partícula puede moverse en el interior de una esfera de radio r. Al poder
ocupar cualquier posición dentro del volumen limitado por la superficie, incluido su
contorno, la relación de vínculo resulta:
x2 + y2 + z2 ≤ r2, x ≤ r, y ≤ r, z ≤ r

Ejemplo 3: Una partícula está sujeta al extremo de una varilla indeformable de longitud
L, y fija en su extremo opuesto P0. La relación de vínculo para el sistema de referen-
cia adoptado es:
(x-x0)2 + (y-y0)2+(z-z0)2=L2
La partícula evolucionará siempre sobre una superficie esférica de radio L.

Observaciones:
1. La relación de vínculos expresa las posibilidades que tiene el movimiento, pero no
caracteriza físicamente al vínculo. Así, una partícula sujeta a una cuerda inextensible de
longitud l, constituye un vínculo físicamente diferente de los indicados en los ejemplos 2 y
2. La relación de vínculo es, sin embargo, igual a la del ejemplo 2.
3. La expresión de la relación de vínculo depende del sistema de referencia que se adopte.

77
 Las consideraciones precedentes permiten la siguiente agrupación: los vínculos
mecánicos pueden ser completos o bilaterales, o incompletos o unilaterales.

El vínculo es completo cuando la partícula está

obligada a permanecer en todo momento en P
contacto con él.
Matemáticamente, la relación de vínculo ds
constituye una ecuación o igualdad. En este tipo
de vínculo para cada desplazamiento ds es
siempre posible que ocurra el opuesto -ds. P
-ds

El vínculo es incompleto cuando solo

impone una condición de frontera entre dos
semiespacios y la partícula queda condicionada
únicamente a no trasponerla, pudiendo evolucionar en cualquier punto del semiespacio
habilitado. El contacto con la frontera es posible, pero no forzoso. Matemáticamente, la
relación de vínculo constituye una inecuación o desigualdad. En estos casos, para algunos
desplazamientos posibles ds2, no son posibles sus opuestos -ds2.

ds
P
P
ds2
-ds

Los vínculos podrán ser rugosos o lisos, según introduzcan o no fuerzas disipativas
de fricción o rozamiento.

También es posible considerar vínculos de naturaleza y acción dinámica variables o

constantes en el tiempo.

Reacciones de vínculos

La restricción al movimiento que introduce un vínculo, constituye una acción dinámica

actuante sobre la partícula. Es decir, se materializa mediante fuerzas conocidas como
Reacciones de vínculo. Así, los vínculos son responsables de la aparición de fuerzas
reactivas, en tanto que a las restantes acciones actuantes sobre la partícula, se denominan
fuerzas activas.
Son características de las reacciones de vínculo:
™ dependen de las fuerzas activas actuantes.
™ Si las fuerzas se anulan, las reacciones de vínculo se anulan.
™ Las reacciones de vínculo no pueden, de por sí, producir movimiento.

Diagrama de cuerpo libre: esquema en que se representan todas las fuerzas actuantes:
activas por interacción de la partícula con otros sistemas, y reactivas introducidas por los
vínculos.

78
Vínculos Rugosos. Rozamiento.

Se ha dicho que los vínculos rugosos son los que introducen fuerzas disipativas de

rozamiento. Estas acciones se oponen al deslizamiento relativo entre partículas y

vínculos, y se ejercen según la dirección tangencial.

Para describir las características del rozamiento de primera especie, es conveniente

considerar un cuerpo (que se puede asimilar a una partícula) apoyado sobre una superficie
rugosa que actúa como vínculo, y se puede señalar que:
♦ La magnitud de la fuerza de rozamiento es independiente de la extensión de la
superficie de contacto.
♦ Depende de la naturaleza de las superficies en contacto.
♦ Es proporcional a la acción de las fuerzas activas normales a la superficie o vínculo.
♦ Equilibra cualquier componente tangencial de la fuerza exterior, es decir, impide
desplazamientos, mientras no se supere el valor límite fN, en que f es el coeficiente de
rozamiento estático (siempre menor que la unidad).

El movimiento en presencia de vínculos

En presencia de vínculos, la ecuación general del movimiento se expresa por:
Fa + Fvin = ma
Fa es la fuerza activa, Fvin es la resultante de todas las reacciones de vínculos. Cada
T N T
una de estas reacciones admite una descomposición (Fvin)i = ( Fvin )i +( Fvin )i, donde ( Fvin )i
N
es la reacción de vínculo tangencial y ( Fvin )i la reacción de vínculo normal.
T
Cuando el vínculo es liso, ( Fvin )i = 0, es decir, la reacción es normal a la superficie
o curva en que se materializa el vínculo.
T
En vínculos rugosos ( Fvin )i es debida a las fuerzas de fricción o rozamiento.

Movimiento de una partícula sobre una curva lisa.

Para una curva lisa y fija se puede considerar el movimiento en coordenadas
intrínsecas r = r(s). La reacción de vínculo actúa en el plano normal a la trayectoria (tiene
componentes según T (versor Tangente) y N (versor normal)). Entonces, si la fuerza activa
es R=RtT + RnN + RbB (B=versor binormal)
Aplicando Fa + Fvin = ma,
Tendremos las siguientes ecuaciones en componentes:
 ••
F
 t = m s
 •
 s2
R n + Fvin = m
N

 ρ
R + F = 0
b
 b vin


79
Movimiento sobre una superficie lisa

Una superficie en el espacio está expresada por : f(x,y,z,t) = 0 [1].

∂f
La superficie será fija cuando: = 0. Si es lisa, debe ser (Fvin)= (FvinN) = λ ∇f ,
∂t
donde λ es un escalar y ∇f es el gradiende de f.
La ecuación fundamental es R + ( FvinN
)= R + λ∇f = ma , (R es la fuerza activa) que
se puede expresar en componentes:
 d2x ∂f
m 2 = R x + λ
 dt ∂x
 d 2 y ∂f
m 2 = R y + λ [2]
 dt ∂y
 d2z ∂f
m 2 = R z + λ
 dt ∂z
Mediante las ecuaciones [1] y [2] se puede resolver el movimiento si se conocen las
condiciones iniciales (r0,v0).

Existirá equilibrio toda vez que la resultante de las fuerzas activas sea nula, o bien
normal a la superficie (caso de vínculo bilateral). En caso de vínculo incompleto, además
es necesario que R tenga un sentido único: el que posibilite presión de la partícula sobre la
superficie.

Dinámica Relativa de la Partícula

En Cinemática hablamos de Sistemas de Referencia Fijos y Móviles. En Dinámica,

a su vez, distinguimos Sistemas Inerciales y No Inerciales. Entre los sistemas móviles
basta una rotación de los mismos para que estén en acelerados y, en consecuencia,
constituyen sistemas no inerciales.
En Cinemática Relativa, hemos establecido la relación entre la aceleración de la
partícula referida a un sistema fijo cuando interviene además, al menos un sistema
móvil (en el caso más general, animado de movimiento de rototraslación).

•• •
a F = a M + R + ω× r M + 2 ω× v M + ω × (ω × r M )
Por la ecuación fundamental de la dinámica: F = ma considerando a un sistema fijo
como Inercial y al sistema en rototraslación como no inercial:
•• •
F = ma F = ma M + m R + m ω× r M + 2m ω× v M + mω × (ω × r M )
En esta expresión, se define a:
F M = ma M (masa x aceleración relativa)

•• •
Fuerza inercial de arrastre: -( m R + m ω× r M + 2m ω× v M + mω × (ω × r M ) = F A )
= (masa x aceleración de arrastre).

80
 Esta fuerza inercial de arrastre constituye una "fuerza correctiva que permite ajustar
la ecuación fundamental cuando está referida a un sistema no inercial, o sea:
F M = ma M = F F − F A
La fuerza inercial de arrastre se compone de una de arrastre por traslación y otra de
arrastre por rotación:
F A = F A ( Traslación ) + F A(Rotación )
••
F A ( Traslación ) = - m R

•
F A(Rotación ) = - ( m ω× r M + 2m ω× v M + mω × (ω × r M ) )
Entre los términos de la fuerza de arrastre por rotación, se distinguen:

2m ω× v M = Fuerza de inercia de Coriolis.

•
m ω× r M = Fuerza de inercia debida a la aceleración de la rotación.

mω × (ω × r M ) = Fuerza de inercia Centrífuga.

Principio de D´Alembert

Se trata de un Principio que permite plantear los problemas de dinámica en forma similar a
problemas de estática, es decir mediante ecuaciones de equilibrio.
En el apartado anterior, se ha introducido el concepto de fuerza inercial, la cual surge de la
multiplicación de la masa de la partícula por una aceleración.
En base a esta idea la ecuación de Newton
F = m a escrita como F – m a = 0
Considerando ahora la cantidad - m a como una fuerza inercial FIN, la ecuación anterior
puede pensarse como una ecuación de equilibrio
F + FIN = 0
Esta última ecuación permite tratar un problema de dinámica como uno de estática
planteando la ecuación de movimiento como la indicada mas arriba.
Esto constituye la base del principio de D´Alembert para un sistema. Cabe resaltar que las
ecuaciones basadas en este principio y las correspondientes a la segunda ley de Newton
son enteramente equivalentes.

81
Unidad 5: Dinámica de los Sistemas de Partículas . Parte I

Introducción
Llamaremos sistema de partículas a aquel m
conjunto de puntos materiales ligados entre si a
través de fuerzas internas.
r mi m
Cada partícula posee una masa mi y respecto de m ri
un sistema referencial se las ubica por sus r
vectores posición. r1

Fuerzas Actuantes:

m Sobre las partículas se encuentran

mj r
ri F jii actuando fuerzas. Estas pueden ser
F21 r exteriores o internas:
r
F12i
Fiji
1-1) Fuerzas Exteriores F e : ( )
r r M mi Provienen de la acción de sistemas o
F3e m rCt r partículas no pertenecientes al sistema
m ri
considerado.
( )
1-2) Fuerzas Internas F i : Son
debidas a la acción mutua entre las
partículas. Las fuerzas internas se
presentan “de a pares”. Por el principio
de acción y reacción Fij = − F ji ; o sea
que la interacción entre partículas tiene
para cada una igual intensidad y recta
de acción pero sentido contrario.

Masa del sistema:

La masa total del sistema es:
n
M = ∑ mi mi: masa de cada partícula, o sea que queda definida como la suma total
i =1
de las masas de las partículas.

Centro de masa:
Se define como momento de primer orden al producto de la masa por la distancia a un
punto.
d e = ri ∗ mi
para todo el sistema:
De = ∑ d e = ∑ ri ∗ mi i=1....n
Podemos encontrar un punto tal que:
De = ∑ ri ∗ mi = M ∗ rg

rei =
∑r ∗m i i
=
∑r ∗m
i i

M ri ∑m i

83
 Se denomina centro de masa, y es un punto en el que se puede suponer concentrada toda
la masa, respecto del momento de primer orden o momento estático.
Las coordenadas cartesianas del centro de masa son:

Xg =
∑x ∗m i i

∑m i

Yg =
∑ y ∗m
i i

∑m i

Zg =
∑z ∗m
i i

∑m i

La ubicación del centro de masa es independiente del sistema referencial elegido.

Cantidad de Movimiento:

r
Si cada partícula mi se mueve con velocidad vi ; la cantidad de movimiento de cada
partícula es:
r r
p = mi ∗ v i
Para todo el sistema es:
r r r
P = ∑ p i = ∑ mi ∗ v i i=1,....n
Sabemos que:
r
rG ∗ M = ∑ ri ∗ mi
La posición de la partícula es función del tiempo; luego, derivando:
r r r r s
r& ∗ M = vG ∗ M = ∑ r&i ∗ mi = ∑ vi ∗ mi = P
Nos queda que:

r r
m P = M ∗ vG
R r
r vri
W r R vi
m ri i
Rr Si el sistema de partículas se encuentra en
v0 un movimiento roto-traslatorio reducido al
polo “O” la velocidad de la partícula mi

r r
vi = vo + W × (Ri − O )

84
 r r
P = ∑ mi ∗ v i
r r
P = ∑ mi ∗ vo + ∑ mi ∗ w × (Ri − O )
r r
P = vo ∗ ∑ mi + ∑ mi ∗ w × (Ri − O )
La cantidad de movimiento: r r
r r
P = vo ∗ M + W × M ∗ rG

r r r
∑m i ∗ w × (Ri − O ) = W × ∑ mi ∗ (Ri − O ) = W × M ∗ rG
Si el movimiento es una traslación pura:
r r r
W = 0 ⇒ P = M ∗ vo
Si el polo de reducción se toma en el centro de masa, ”O” coincide con “G” y:
r r r
rG = 0 ⇒ P = M ∗ vG
Y si es una rotación pura de “O” que no coincide con “G”:
r r r
P = M ∗ W × rG
Conservación de la cantidad de movimiento:
Aplicando el principio de inercia a cada partícula:
r r r
∂p i r& r
= pi = mi ∗ v&i = Fi e + Fiji
∂t
Sumando para todo el sistema:
r r& r r r
∑ p& i = P = ∑ mi ∗ v&i = ∑ Fi e + ∑∑ Fiji
i j

Pero:
ri
∑∑ F
i j
ij =0 Por el principio de acción y reacción
re r
∑F i= Fe Fuerza exterior total
Nos queda:
r r&
Fe = P
“La suma de todas las fuerzas exteriores es igual a la derivada de la cantidad de
movimiento del sistema respecto del tiempo”.
r& r
P = Fe = 0
Si r
P = cte

Esta ecuación Fext = M dvG/dt, indica que dado un sistema de partículas, el centro de
masas de dicho sistema se mueve como una partícula de masa igual a la masa total del
sistema, sometida a una fuerza igual a la fuerza exterior total que actúa sobre el sistema.

85
 Si la suma de las fuerzas exteriores es nula, la cantidad de movimiento del sistema no
varía. Desde el punto de vista del centro de masas, esta continuará en su estado de reposo o
movimiento rectilíneo uniforme.
Ejemplo: Dos masas que colisionan entre sí, en ausencia de fuerzas externas. Sea el
choque elástico o plástico, la cantidad de movimiento se conserva porque no actúan fuerzas
exteriores al sistema. En cuanto al centro de masas, esta posee movimiento rectilíneo
uniforme antes de la colisión, y conserva su estado después. Si el choque es plástico, las
dos masas unidas continúan en la trayectoria del centro de masa.

Momento Cinético:
El momento cinético de una partícula respecto de
un punto “O” es: r
r r vi
Li = (Ri − O ) × mi ∗ vi m
r
Para un sistema de partículas: ri
R
r r r
L = ∑ Li = ∑ [(Ri − O ) × mi ∗ vi ]
Si el sistema se halla sometido a un movimiento
roto-traslatorio con polo de reducción en “O”
r r
vi = vo + W × (Ri − O )

El momento cinético total respecto de “O”

r r
L = ∑ [(Ri − O )× mi ∗ (vo + W × (Ri − O ))]
r r
L = ∑ (Ri − O ) × mi ∗ vo + mi ∗ (Ri − O ) × W × (Ri − O )
r r
L = ∑ (Ri − O ) × mi ∗ vo + ∑ mi ∗ (Ri − O ) × W × (Ri − O )
En esta ecuación tenemos:
r r r
∑ (Ri − O ) × mi ∗ vo = M ∗ rG × vo
Nos queda:
r r r
L = M ∗ rG × vo + ∑ mi ∗ (Ri − O ) × W × (Ri − O )
Llamaremos Momento Cinético de Rotación Pura (L0) a:
r
L0 = ∑ mi ∗ (Ri − O ) × W × (Ri − O )
Entonces:
r r r r
L = M ∗ rG × vo + L0
Conservación del momento cinético
r r r
∂p r r
Vimos que para cada partícula: i = p& i = mi ∗ v&i = Fi e + Fiji
∂t
Tomando momento respecto de “O”
r r
(R i − O ) × m i ∗ ar i = (R i − O ) × F i e + (R i − O ) × F iji
La expresión del momento cinético:
r r r r r
L = M ∗ rG × vo + L0 = ∑ (Ri − O ) × mi ∗ vi
Derivando:

86
 ∂L0 r r  ∂R ∂O  r r
+ M ∗ r&G × vo + M ∗ rG × v&o = ∑  i +  × mi ∗ vi + ∑ (Ri − O ) × mi ∗ a i
∂t  ∂t ∂t 
Llamando
r r
M e = ∑ (Ri − O ) × Fi e
r r
M e = ∑ (Ri − O ) × Fiji = 0 Por principio de acción y
reacción
Nos queda:
r
re r& r r r r r ∂L0
M = M ∗ rG × vo + M ∗ rG × a o + vo × p +
∂t
Si el polo de reducción se elige en “O” = “Centro de masas”, el vector posición del
centro de masas y sus derivadas serán nulas, y como la cantidad del movimiento puede
expresarse en función de la velocidad del centro de masas, el tercer término de la ecuación
se anula. Resulta entonces:
r r
r e ∂L0 re ∂L0 r
M = ⇒M =0⇒ = 0 ⇒ L0 = Cte.
∂t ∂t
El momento angular de un sistema de partículas se mantiene constante si el momento
exterior que actúa sobre el mismo es nulo.

Por otra parte, si el momento exterior no es nulo, como la derivada del momento
angular respecto del tiempo es igual al momento exterior, se interpreta que la variación del
momento angular coincide en dirección con el momento exterior.
Ejemplo: Considere el movimiento de un trompo que rota a gran velocidad sobre su eje.
Si bien es una aproximación, se puede asumir que el momento angular coincide en
dirección con el vector rotación del trompo. El momento exterior que actúa sobre el
trompo es el debido al peso del mismo. El momento angular variará de manera tal que dL
coincida en dirección con el torque. Se origina, por lo tanto, el movimiento del precesión,
rotación alrededor de un eje vertical. La dirección del vector rotación de precesión depende
de la dirección del vector rotación propia del trompo.

Trabajo y Energía:
r
Si Fi e son las fuerzas actuantes sobre las partículas mi; en un desplazamiento elemental:
r r r r
∂Ti = Fi ∗ ∂ri Como ∂ri = vi ∗ ∂t
r r
∂Ti = Fi ∗ vi ∗ ∂t
r r
∂T = ∑ Fi ∗ vi ∗ ∂t
Para un sistema roto-traslatorio
r r
vi = vo + W × (Ri − O )
r r r
∂T = ∑ Fi ∗ v o ∗ ∂t + ∑ Fi ∗ (W × (Ri − O )) ∗ ∂t
r r r
( )
∂T = F e ∗ vo + M e ∗ W ∗ ∂t
Dividiendo todo por “dt” obtenemos la potencia

r r r
PW =
∂T
∂t
(
= F e ∗ vo + M e ∗ W)
Definimos a la energía cinética de un sistema como la suma de las energías cinéticas
de las partículas que componen el sistema.

87
 1 r
Ti = ∗ mi ∗ vi2
2
Para todo el sistema
1 r
T = ∑ ∗ mi ∗ vi2
2
En un movimiento roto-traslatorio
r r
vi = vo + W × (Ri − O )
1 r
T = ∗ ∑ mi ∗ [vo + W × (Ri − O )]
2

2
1
[ r r
]
T = ∗ ∑ mi ∗ vo2 + ∑ mi ∗ 2 ∗ vo ∗ (W × (Ri − O )) + ∑ mi ∗ (W × (Ri − O ))
2
2

1 r r 1
T = ∗ M ∗ vo2 + vo ∗ ∑ mi ∗ (W × (Ri − O )) + ∗ ∑ mi ∗ (W × (Ri − O ))
2

2 2
1 2 3
1 r2
1) T1 = ∗ M ∗ vo Energía cinética de traslación
2
r 1 r
1
( ) ( r
3) T3 = ∗ ∑ mi ∗ W × (Ri − O ) = ∗ W 2 ∗ ∑ mi ∗ eWV × (Ri − O )
2
2

2
) 2

( rV
)
eW × (Ri − O ) : Cuadrado de la distancia
2

( r
)
I W = ∑ mi ∗ eWV × (Ri − O ) : Momento de inercia respecto del eje ew
2

1 r
Nos Queda T3 = ∗ W 2 ∗ I W
2
r r
2) T2 = vo ∗ ∑ mi ∗ (W × (Ri − O )) = vo ∗ W × ∑ mi ∗ (Ri − O )
r r
T2 = M ∗ (v o ∗ W × rG ) : Energía cinética mixta
Esta cantidad se anula cuando el centro de referencia se encuentra en reposo, dado que para
r
ese caso vo = 0 o cuando coincide con el centro de masas, por ser rG
En estos casos la expresión de la energía toma la forma simplificada
1 r 1
T = M ∗ vG2 + ∗ W 2 ∗ Iw
2 2

88
Dinámica del Sólido Rígido
Introducción:
Consideramos al sólido rígido como aquel cuerpo constituido por un sistema de
partículas que cumplen con la condición de rigidez detenida, es decir, mantienen sus
distancias relativas constantes.

Masa del sólido rígido:

Para un sistema de partículas, la masa total es:
M = ∑ mi
i
Para el sólido rígido:
M = ∫ ρ ∗ dv ρ: Densidad
v
Se entiende que el proceso de integración se extiende a todo el volumen ocupado por el

sólido, es decir, se involucra una integral triple en un marco referencial adecuado.

si ρ = Cte.: M = ρ ∗ ∫ dv = ρ ∗ V
v

Centro de masa:
El momento de primer orden es:
De = ∫ mi ∗r i
v
Y el centro de masa: r m
r rG
M ∗ rG = De = ∫ mi ∗r i r
ri
v

r ∫ ρ ∗r ∗dv i

rG = v

∫ ρ ∗ dv
v

Cantidad de Movimiento:
Si el sólido se encuentra en movimiento rototraslatorio reducido al polo “O” la
velocidad de cada punto será:
r r
vi = vo + W × (Ri − O )
La cantidad de movimiento total:
r r r
P = ∫ ρ ∗ dv ∗ vi = ∫ ρ ∗ dv ∗ vo + ∫ ρ ∗ dv ∗ (W × (Ri − O ))
v v v
r r
P = vo ∗ ∫ ρ ∗ dv + W × ∫ ρ ∗ dv ∗ (Ri − O )
v v

r r r ∫ ρ ∗ dv ∗ (R
i − O ) = De
P = vo ∗ M + W × M ∗ rG (1) donde: v
r
De = M ∗ rG

89
Conservación de la cantidad de movimiento:
Derivando (1):
r r
∂P r ∂W r r r
= ao ∗ M + × M ∗ rG + W × M ∗ r&G
∂t ∂t r r
W vi
Sabemos que la variación de la cantidad de r
ri
movimiento respecto al tiempo es:
r r
v0 r R
∂P r e r e rG
= Fi = F
∂t
Entonces nos queda:
r r r& r r r
F e = a o ∗ M + W × M ∗ rG + W × M ∗ vG

Momento cinético:
Para un sistema de partículas teníamos, en un
movimiento rototraslatorio con polo de
reducción en “O”:
r r
r r L = ∑ ( Ri − 0 ) × m i ∗ v i
W vi i

r Para un sólido:
r r r
r
v0
ri L = ∫ ρ ∗ (r × v ∗ dv )
r R r r
v
rG v = v ( x, y , z )
r r
r = r ( x, y , z )
r r r r
Pero: v = vo + W × r
Reemplazando, el movimiento cinético nos
queda:
r r r r
r r
(
L = ∫ ρ ∗ r × vo ∗ dv + ∫ ρ ∗ r × W × r ∗ dv )
v v
r r
L1 L2
r
La primera integral ( L1 ) es:
r r r r r r r
L1 = ∫ ρ ∗ r × vo ∗ dv = −vo × ∫ ρ ∗ r ∗ dv = −vo × M ∗ rG
v v
r r r r
L1 = −vo × M ∗ rG Se anula si vo = 0 o cuando la reducción
se hace al centro de masa.
r
La segunda integral ( L2 ) es:
r r r r
( ) [
r r r r r r
L2 = ∫ ρ ∗ r × W × r ∗ dv = ∫ ρ ∗ (r ∗ r ) × W − r ∗ W ∗ r ∗ dv ( ) ]
v v
Expresando en coordenadas cartesianas:
r
[( ) (
r r r
) r r r
(
L2 = ∫ ρ ∗ x 2 + y 2 + z 2 × W x i + W y j + W z k − (x ⋅ W x + y ⋅ W y + z ⋅ W z ) ∗ xi + yj + zk ∗ dv )]
v
r
Efectuando los productos, operando y poniendo: L 2 : ( L2 x , L 2 y , L2 z )
Teniendo presente que se trata de una igualdad vectorial.
La componente:

90
 r
( )
L x = W x ∗ ∫ ρ ∗ y 2 + z 2 ∗ dv − W y ∗ ∫ ρ ∗ ( yx ) ∗ dv − W z ∗ ∫ ρ ∗ ( zx ) ∗ dv
v v v

Momento y productos de inercia:

Definiremos el momento de inercia respecto de un eje “e” como:
Para un sistema de partículas:
I e = ∑ mi ∗ d i2 Donde
i
di : Distancia de la partícula al eje
considerado
- Para un sólido:
I e = ∫ ρ ∗ d 2 ∗ dv
v
di
mi Definiremos como producto de
inercia respecto de los ejes “e” y
“h”
I eh = −∑ mi ∗ d ie ∗ d ih
i
Para un sólido:

I eh = − ∫ ρ ∗ d e ∗ d h ∗ dv
v
Los productos de inercia son:

 I xy = − ∫ ρ ∗ x ∗ y ∗ dv
 v
 I = − ρ ∗ x ∗ z ∗ dv
 xz ∫v mi
 d ih d ie
 I yz = − ρ ∗ y ∗ z ∗ dv
 ∫v
En las fórmulas de cálculo de los
productos de inercia se trabaja con
coordenadas, es decir, pueden tener
signo positivo o negativo. Esto hace
que el resultado de la integral de
volumen puede resultar nulo. De
hecho, esto ocurre cuando los ejes son
ejes principales de inercia.

r
La expresión del momento cinético L2 nos queda:
r
 L x = I xx ⋅ W x + I xy ⋅ W y + I xz ⋅ W z
r
 L y = I yx ⋅ W x + I yy ⋅ W y + I yz ⋅ W z
r
 L z = I zx ⋅ W x + I zy ⋅ W y + I zz ⋅ W z

91
 Estas expresiones escalares podemos poner en forma matricial:

[ ]
r T
L2 = [I ] ⋅ [W ] Siendo:
T
r
L 2 =
r
L 2x +
r
L[ ] (
2y +
r
L 2z )
[W ] = (Wx ,W y ,Wz )
Y la matriz [I ] es la representación en el sistema referencial adoptado del tensor de
inercia.
I xx I xy I xz
[I ] = I yx I yy I yz
I zx I zy I zz
r r r
Finalmente la expresión del momento cinético es: L = L1 + L2
r r r
L = M ⋅ rG × vo + [I ] ⋅ [W ]
T

Carácter Tensorial de [I ]
Supongamos la matriz [ A] que transforma del sistema referencial x, y, z ortogonal al x’,
y’, z’
también ortogonal.
[A] es ortogonal. Luego [A]−1 = [A]T
Tenemos que:
r T
[ ]
L2 = [I ] ⋅ [W ]
T
(1)

Transformando:
r
[] r T
L ' = [ A] ⋅ [L ] ⇒ [L ] = [ A] ⋅ L '
T T T
[]
r
[ ] r T
[ ]
r T
W ' = [A] ⋅ W ⇒ W = [ A] ⋅ W '
T
[ ] r T
[ ]
Reemplazando en (1):
r T
[] r T
[A]T ⋅ L' = [I ] ⋅ [A]T ⋅ W ' [ ]
Multiplicando por [ A]

[A] ⋅ [A]T ⋅ [L '] [ ]

r r T
= [A] ⋅ [I ] ⋅ [ A] ⋅ W '
T T

[U ] [T ]
Nos queda
[]
r T
L ' = [I '] ⋅ [W ']
T
Expresión en x’, y’, z’
Donde
[I '] = [A] ⋅ [I ] ⋅ [A]T
Vemos que [I ] se transforma según la ley de transformaciones de tensores de segundo
orden, se concluye que [I ] es la representación de un tensor.

Operador Lineal
r r r r r
La expresión de L = L1 + L2 = M ⋅ rG × vo + [I ] ⋅ [W ]
T

r
[]
Tiene carácter tensorial, es decir que L es un invariante de los distintos sistemas
referenciales.
r
[ ] r r
L1 : Proviene del producto vectorial rG × vo y como tal tiene carácter invariante.
r
[ ]
L1 : Proviene de la aplicación del operador lineal [I ] (tensor de segundo orden)
sobre el tensor de primer orden [W ] ; luego, es un invariante.
T

92
 r r r
L = L1 + L2 : Es un invariante por ser la suma de dos invariantes.

Conservación Del Momento Cinético

r r r r r
El momento cinético es L = L1 + L2 = M ⋅ rG × vo + [I ] ⋅ [W ] .
T

Derivando con respecto al tiempo:

r
∂L
∂t
r r r r
= M ⋅ r&G × vo + M ⋅ rG × a o + [I ] ⋅ W&[ ]
T

r
∂L r
Pero sabemos que =Me
∂t
re r& r
G o
r r
M = M ⋅ r × v + M ⋅ r × a + [I ] ⋅ W&
G o [ ]
T

Energía Cinética
La energía cinética para un sistema de partículas
1 r r 1
T = ∗ M ∗ vo2 + vo ∗ ∑ mi ∗ (W × (Ri − O )) + ∗ ∑ mi ∗ (W × (Ri − O ))
2

2 2
Para el caso de un sólido nos queda:
r r
1 r
( r
)1
T = ∗ M ∗ vo2 + M ∗ v o ∗ W × rG + ∗ I w ∗ W 2
2 2
Si el polo de reducción coincide con el centro de masa:
1 r 1
T = ∗ M ∗ vG2 + ∗ I w ∗ W 2 = Tt + TR
2 2
Tt: Energía de traslación
1
TR: Energía de rotación ∗ I w ∗ W 2 = TR
2
Expresando en función de P y L
r
r
[] []
r
r
P = M ∗v o
r r2 r r r r
1
( )
TR = ∗ ∫ ρ ∗ dv ∗ W × r = ∫ ρ ∗ dv ∗ W × r ∗ W × r
2 v
( )( )
v
Desarrollando en coordenadas cartesianas y operando obtenemos:
1
TR = [W ] ∗ [I ] ∗ [W ] Donde [I ] : Tensor de inercia
T

2
[W ] ∗ [I ] ∗ [W ]T Cuadrática de energía asociada
al tensor de inercia
Finalmente:
1 r r r r r 1 r rT
T = ∗ P × vo + P ∗ W × rG + ∗ W L
2 2
[ ][ ]
Cuadráticas asociadas al tensor de Inercia
11-1) Energía cinética de rotación
eWV
1
Vimos que TR = [W ] ∗ [I ] ∗ [W ] es invariante.
T

2
11-2) Momento de inercia respecto de un eje:
1 1
TR = [W ] ∗ [I ] ∗ [W ] = ∗ I w ∗ W 2
T

2 2
IW: Momento de inercia

93
respecto del eje de dirección e wv
1/
[W ] ∗ [I ] ∗ [W ]T = 1/ ∗ I w ∗ W 2
2/ 2/
[W ] ∗ [I ] ∗ [W ] = I
T

w
W W
[ew ] ∗ [I ] ∗ [ew ]T = I w
O sea que el momento de inercia respecto de un eje es la cuadrática asociada al tensor
con respecto al versor de dirección de ese eje.
ewv = (α , β , γ )
I w = I xx ∗ α 2 + I yy ∗ β 2 + I zz ∗ γ 2 + 2 ∗ α ∗ β ∗ I xy + 2 ∗ α ∗ γ ∗ I xz + 2 ∗ β ∗ γ ∗ I yz

Radio de Giro
Se define el radio de giro de un sólido, respecto de un eje cuya dirección es e wv como:
Iw
rk =
M
Al momento de inercia lo podemos expresar:
I w = M ∗ r 2k
Si se divide en forma cuadrática [ew ] ∗ [I ] ∗ [ew ] = I w por I w
T

[X ] ∗ [I ] ∗ [X ] = 1 Donde [X ] = [e ] = [e ]
v v
T w w

Iw M ∗ rk
En forma desarrollada:
I w = I xx ∗ X 2 + I yy ∗ Y 2 + I zz ∗ Z 2 + 2 ∗ X ∗ Y ∗ I xy + 2 ∗ X ∗ Z ∗ I xz + 2 ∗ Y ∗ Z ∗ I yz
Es la ecuación de un elipsoide denominado elipsoide de inercia donde (x, y, z) son las
componentes del vector:
r 1 1
r= ∗ env = ∗ env
Iw M ∗ rk
Es decir que todo semieje del elipsoide de inercia es proporcional a la inversa del radio
de giro.

Diagonalización del tensor de Inercia

El tensor [I ] se transforma mediante una rotación ortogonal
[I '] = [A] ∗ [I ] ∗ [A]T
Existe una transformación [M ] tal que:
[I '] = [M ] ∗ [I ] ∗ [M ]T = [I d ]
Donde:
I1 0 0
[I d ] =
0 I2 0
0 0 I3
I1, I2, I3 son los autovalores del tensor de inercia, denominados momentos de inercia
principales.
Por tratarse el tensor de inercia de un tensor real y simétrico los autovalores son reales.

94
 Los vectores propios asociados V1, V2, V3 establecen direcciones denominadas ejes
principales de inercia.
[I ] ∗ [U ]T = I ∗ [v]T
Expresiones en la base principal
Para los ejes principales de inercia, las variables que intervienen tienen expresiones más
simples.

Momento cinético
Lx I1 0 0 Wx
[L] = [I 0 ] ∗ [W ] ⇒ L y = 0
T T
I2 0 ∗ Wy
Lz 0 0 I 3 Wz
Nos queda:

Lx = I 1 ∗ Wx ; Ly = I 2 ∗Wy ;
Lz = I 3 ∗ Wz

O sea que en la base principal el momento

cinético tiene la dirección del eje de rotación.

Las cuadráticas asociadas al tensor también adquieren una forma simple.

Momento según un eje.

I e = [ee ] ∗ [I d ] ∗ [ee ] = I 1 ∗ α 2 + I 2 ∗ β 2 + I 3 ∗ γ 2
T

Energía cinética de rotación.

2 ∗ Tr = [W ] ∗ [I 0 ] ∗ [W ] = I 1 ∗ W x + I 2 ∗ W y + I 3 ∗ W z
T 2 2 2

Elipsoide de inercia.
1 = [x ] ∗ [I 0 ] ∗ [x ] = I 1 ∗ X 2 + I 2 ∗ Y 2 + I 3 ∗ Z 2
T

Teorema de los ejes paralelos

Sea un sólido e [IG] su tensor de inercia referido a una terna x, y, z cuyo origen coincide
con el centro de masa G del sólido.
Si se desea referir el tensor a los ejes x´,y´,z´, paralelos a la terna x, y, z, se aplica el
teorema de los ejes paralelos. De acuerdo a este teorema se debe adicionar a [IG] un tensor
de inercia correspondiente a la masa total del sólido concentrada en el centro de masas,
referido a x´, y´, z´,. Es decir que : [Ix´, y´,z´] = [IG] + [ITRASLACION] donde: el tensor de
inercia referido a la terna x´,y´,z´ es igual a la suma de el tensor de inercia baricéntrico mas
el tensor de inercia de la masa concentrada en el centro de masa referido a x´,y´,z´.
Utilizando el producto tensorial esta operación de traslación se puede realizar mediante la
expresión
[Ix´, y´,z´] = [IG] + M (rG2 [ID]- rG ⊗ rG )
Donde:
M es la masa total del sólido
rG el vector posición del centro de masa respecto al sistema x´,y´,z´.
ID la matriz identidad.
95
Dinámica relativa del sólido
Introducción
Las causas del movimiento del sólido en el espacio son las fuerzas exteriores y sus
momentos. Estas cantidades se relacionan con las derivadas respecto del tiempo del
momento lineal o cantidad de movimiento y el momento angular o momento cinético del
sólido. La cantidad de movimiento puede expresarse en función de la masa total del
sistema y la velocidad del centro de masa, en tanto que el momento angular depende del
tensor de inercia del sólido. Si el tensor de inercia está referido a ejes fijos, sus
componentes no permanecen constantes al cambiar el sólido de posición, por lo que el trato
de las ecuaciones que ligan la variación del momento angular con el momento de las
fuerzas exteriores puede ser complicado. Esta dificultad puede eludirse introduciendo
sistemas referenciales relativos, fijos al cuerpo en estudio. Luego, la derivada de
magnitudes expresadas en este referencial móvil pueden obtenerse aplicando el operador
de derivación para sistemas móviles en rototraslación.

En lo que sigue, trataremos el movimiento del sólido, con sus distintas variables
dinámicas, referidos a sistemas de ejes en movimiento respecto de una terna fija.

Sólido referido a un sistema rototraslatorio

Supongamos un sólido referido a un sistema
que se encuentra en un estado de rototraslación
respecto de una terna fija. El sólido se encuentra
fijo respecto del sistema móvil.
r
vo
Cantidad de Movimiento
Habíamos visto que la expresión de la cantidad
de movimiento referida al sistema x, y, z, era:
r r
P = vro ∗ M + W × M ∗ rrG
m
Derivando respecto del sistema fijo, recordando
que la derivación relativa de un vector se hacía
mediante:
d d r
= + (W × )
dt F dt m
Para el vector cantidad de movimiento:
r r
r r
+ W × P  =
dP dP
=
dt F dt m  m
r r& r r r r r r r r
v& o ∗ M + W × M ∗ rG + W × M ∗ r&G + W × v o ∗ M + W × W × M ∗ rG
r
r&G = 0 , pues el sólido se encuentra fijo respecto de x, y, z. Sabemos que las fuerzas
aplicadas respecto de un sistema fijo son:
r
re re dP
∑ F = F F = dt
F
Nos queda:
r r r& r r r r r r
Fe = a o ∗ M + W × M ∗ rG + W × v o ∗ M + W × W × M ∗ rG
F

96
Se verá que esta ecuación se simplifica en forma notable cuando el sólido rota alrededor de

un eje fijo, y se toma el punto “O” sobre el mismo. En ese caso, será:

r r r r r r
Fe = W × M ∗ rG + W × W × M ∗ rG
F
r r
por ser a o = 0 y v o = 0

Momento Cinético

La expresión de momento cinético en el sistema en movimiento es:

r r r
L = M ⋅ rG × vo + [I ] ⋅ [W ]
T
m
Se considera que el tensor de inercia ha sido referido al sistema móvil solidario al
sólido, y por lo tanto sus componentes permanecen constantes al cambiar el sólido de
posición.
Derivando respecto del sistema fijo:
r r
r r
+ W × L  =

dL dL
=
dt dt  m
F m

r r r r r&
M ∗ r&G × v o + M ∗ rG × a o + [I ]⋅ W [ ] + Wr × M ∗ rr × vr + Wr × [I ]⋅ [Wr ]
T
G o
T

Teniendo presente la relación entre la derivada del momento angular y el momento de

las fuerzas externas:
r
dL
=Me
dt F
F

Nos queda:
Me
F
r r r r
[ ] + Wr × M ∗ rr × vr + Wr × [I ]⋅ [Wr ]
r&
= M ∗ r&G × v o + M ∗ rG × a o + [I ] ⋅ W
T
G o
T

e
Rotación de un Sólido con eje fijo
En este caso el sólido mantiene su eje
de rotación fijo en el espacio. r
Entonces: rG
r
vo = 0
r
ao = 0
Tendremos en este caso:
r r& r r r r e
Fe = W × M ∗ rG + W × W × M ∗ rG
F

M e
F
[]r& T r
[ ]
r T
= [I ] ⋅ W + W × [I ] ⋅ W

97
Reacciones dinámicas de un sólido en rotación

Supongamos tener un sólido rotando sobre un eje fijo en el espacio. Este eje lo
consideramos vinculado al sistema fijo a través de cojinetes. Como ejemplo, puede
considerarse el rotor de una turbina hidráulica o el impulsor de una bomba de agua. Se
trata de sólidos montados sobre un eje, el cual es sustentado por cojinetes que permiten la
rotación.
El sólido en su rotación puede generar fuerzas inerciales que presionan sobre los
vínculos. Estas presiones dinámicas (debidas al movimiento) no existen si el sólido se
encuentra estático, en cuyo caso solo habrá reacciones debidas al peso u otras fuerzas
estáticas.
Las presiones dinámicas aparecen cuando el sólido se encuentra en rotación, y así como
las cargas estáticas originan reacciones en los soportes, estas cargas dinámicas harán
aparecer reacciones, que se denominan reacciones dinámicas. La magnitud de estas
reacciones estarán directamente ligadas a la cantidad de masa del sólido, a la forma en que
la masa está distribuida y a las condiciones cinemáticas del movimiento.

Las condiciones dinámicas de este sólido establecen:

r r& r r r r
Fe = W × M ∗ rG + W × W × M ∗ rG
F

M e
F
[]r& T r r T
= [I ] ⋅ W + W × [I ] ⋅ W [ ]
En ambas ecuaciones aparece un término que depende de la aceleración angular, su
existencia se comprende con facilidad puesto que para acelerar una masa, es necesaria una
fuerza o un momento. Estos términos se anulan en el caso de que el sólido rote a velocidad
constante. El otro término involucrado, depende solamente de la velocidad de rotación, de
la cantidad de masa y la forma en que está distribuida, resumidas en el vector posición del
centro de masas y el tensor de inercia.
Las acciones exteriores están formadas por acciones activas y acciones reactivas, es
decir:
r r r
F e = Fae + FRe
r r r
M e = M ae + M Re
El equilibrio del sólido es:
r r r r& r r r r
F e = Fae + FRe = W × M ∗ rG + W × W × M ∗ rG
r r r r& T r
[]
M e = M ae + M Re = [I ] ⋅ W + W × [I ] ⋅ W
r T
[ ]
Que nos provee las ecuaciones para obtener el
valor de las reacciones del sólido.
Por ejemplo si:
r
W e = (0 0 W z ) = cte
r
W e = cte
I xx I xy I xz 0
r r r
M e = W × [I ] ⋅ W [ ] T r
= W × I yx I yy I yz ⋅ 0
I zx I zy I zz Wz
r
= W × (I xz ⋅ Wz I yz ⋅ Wz I zz ⋅ Wz )

98
 i j k
r r r
Me = 0 0 Wz = − I yz ⋅ W z2 ⋅ i + I xz ⋅ W z2 ⋅ j
I xz ⋅ W z I yz ⋅ W z I zz ⋅ W z

Planteamos:
r r r r
M ae + M Re = − I yz ⋅ W z2 ⋅ i + I xz ⋅ W z2 ⋅ j
Debemos tener en cuenta que los vectores de la izquierda de la igualdad están en el
sistema fijo, mientras que los de la derecha están en el sistema en movimiento, lo cual

Fae

M ae

involucra una transformación de coordenadas.

r r r r r r
Para nuestro ejemplo, de acuerdo a los sistemas adoptados i j k ≡ I J K
Luego la transformación es una identidad. Esta ecuación vectorial implica:
r re
M axe + M Rx = − I yz ⋅ W z2
r re (1)
M aye + M Ry = I xz ⋅ W z2
r r r r
Las otras ecuaciones nos proveen: F e
= W × W × M ∗ rG
r
Si las coordenadas del centro de masa son: rG = xG yG xG
i j k
r r r r
W × M ∗ rG = 0 0 W z = − yG ⋅ W z ⋅ i + xG ⋅ W z ⋅ j
xG yG zG
i j k
r r r r r
W × W × M ∗ rG = 0 0 W z = xG ⋅ W z2 ⋅ i − y G ⋅ W z2 ⋅ j
− yG ⋅ W xG ⋅ W z 0
r r r
Nos queda: F e = xG ⋅ W z2 ⋅ i − y G ⋅ W z2 ⋅ j
Que nos da las ecuaciones:
r r
Faxe + FRxe = xG ⋅ W z2
r r (2)
Faye + FRye = y G ⋅ W z2
Que resuelven junto con (1) el estado de equilibrio.
r r
Vemos en las ecuaciones (1) y (2) que para que sea F e = 0 y M e = 0 , deben ser
99
 I xy = I xz = I yz = 0
xG = y G = z G = 0

Para que se cumplan estas condiciones, el eje de rotación debe coincidir con un eje
principal de inercia y debe contener al centro de masas del sólido.
Los momentos centrífugos no nulos y el desplazamiento del centro de masa dan origen
a reacciones dinámicas. Se puede concluir que estas magnitudes “miden” el grado de
desbalanceo que posee el sólido.

Equilibrio dinámico

Dado un sólido en desequilibrio dinámico, se busca equilibrarlo de tal manera que

sea posible eliminar las reacciones dinámicas que, en general, dan origen a la aparición de
vibración en el sistema. Por ejemplo, las ruedas de un automóvil deben “balancearse” luego
de un cierto uso, puesto que el desgaste asimétrico
da origen a reacciones dinámicas. El balanceo se
realiza agregando masas en el caso de las ruedas,
también se puede remover o quitar masa. El ob-
jetivo de este agregado o remoción de masas, es
llevar al eje de rotación a ser nuevamente un eje
principal de inercia conteniendo el centro de masas.
A continuación se plantea el problema del
balanceo o equilibrio dinámico de un sólido en
rotación.
Dado el sólido con:
I xz ≠ 0
I yz ≠ 0
X G , YG , Z G ≠ 0
Se buscará anular los productos de inercia que
tienen un subíndice coincidente con el eje de rotación, y los componentes de vector
posición del centro de masas de subíndice distintos del eje de rotación.
Se ubican las masas m1 y m2 de tal manera que:
I ' yz = I yz + m1 ⋅ Y1 ⋅ Z 1 + m2 ⋅ Y2 ⋅ Z 2 = 0
I ' xz = I xz + m1 ⋅ X 1 ⋅ Z 1 + m2 ⋅ X 2 ⋅ Z 2 = 0
Y 'G = M ⋅ YG + m1 ⋅ Y1 + m2 ⋅ Y2 = 0
X 'G = M ⋅ X G + m1 ⋅ X 1 + m2 ⋅ X 2 = 0
Tenemos ocho incógnitas (m1, m2, X1, X2, Y1, Y2, Z1, Z2), y solo son cuatro las
ecuaciones. Esto se soluciona fijando cuatro y despejando las otras. En general se fija
posiciones a veces limitadas por la geometría y del sólido a “balancear”. Por ejemplo, en
las ruedas del automóvil, los contrapesos pueden ubicarse únicamente en la periferia de la
llanta.

Movimiento de un sólido con un punto fijo

Ecuación de Euler
Vamos a considerar ahora el movimiento de un sólido que mantiene un punto fijo; es
decir un movimiento polar.

100
 Supondremos un sistema fijo y un
sólido vinculado a dicho sistema a
través de un punto, y otro móvil y
solidario con el sólido y en las
direcciones principales de inercia.
El sólido se encuentra vinculado
en su centro de gravedad:
r r
W = W(t ) o sea que la velocidad
angular varía en módulo y dirección.
La ecuación que regula el
movimiento es:
r
M e = [I ] ⋅ W[ ] + Wr × [I ]⋅ [Wr ]
r& T T

I1 0 0
donde: [I ] = 0 I2 0
0 0 I3
Haciendo las operaciones obtendremos:
r r
W × [I ] ⋅ W [ ] = (I
T r
⋅ Wz ⋅ W y − I 2 ⋅ W y ⋅ Wz ) ⋅ i +
3
r
(I 1 ⋅ W x ⋅ W z − I 3 ⋅ W z ⋅ W x ) ⋅ j +
r
(I 2 ⋅ W y ⋅ W x − I 2 ⋅ Wx ⋅ W y ) + k
r r r r
(
M e = M xe , M ye , M ze )
r
M xe = I 1 ⋅ W& x + (I 3 − I 2 ) ⋅ W y ⋅ W z
r r
Me M ye = I 2 ⋅ W& y + (I 1 − I 3 ) ⋅ W x ⋅ W z Ecuaciones de Euler
r
M ze = I 3 ⋅ W& z + (I 2 − I 1 ) ⋅ W y ⋅ W x

Se encuentra expresado en los versores en movimiento, si queremos expresarlos en el

sistema fijo debemos aplicar la transformación correspondiente.

Movimiento de un Sólido por Inercia

Analizaremos el movimiento del sólido
que no está sometido a fuerzas y momentos.
El sólido está vinculado a su centro de masa
y lo mantiene fijo a través de un apoyo. La
resultante del peso se equilibra con la
reacción de apoyo.
r r
Fae + FRe = 0
Establecemos un sistema referencial
solidario con el sólido que coincide con los
ejes principales de inercia.
I1 0 0
[I ] =
0 I2 0
0 0 I3
Como no hay pares exteriores aplicados:
101
 r r
r dL
Me =0 =0 L = cte
dt F

Vamos a suponer un sólido giroscópico, es decir un sólido donde:

I1 = I 2 = I
Las Ecuaciones de Euler quedan:
0 = I ⋅ W& x + (I 3 − I ) ⋅ W y ⋅ W z
0 = I ⋅ W& y − (I 3 − I ) ⋅ W x ⋅ W z
0 = I 3 ⋅ W& z
De la tercera se deduce que:
I 3 ⋅ W& z = 0 W& z = 0 W z = cte
O sea que la rotación sobre el eje “z” se realiza en ausencia de momento exterior y bajo
r
condición L = cte con W z = cte
Las otras dos ecuaciones dan:
0 = W& x + 3
(I − I ) ⋅ W ⋅ W
y z
I
(I − I )
0 = W& y − 3 ⋅ Wx ⋅ Wz
I
W& x + k ⋅ W y = 0 (I − I )
Se puede poner: (1) Donde: k = 3 ⋅ Wz
W& y − k ⋅ W x = 0 I
Derivando la primera: W&& + k ⋅ W& = 0
x y

Introduciendo en la segunda: W x + k 2 ⋅ W& x = 0

&&
Cuya solución es: W x = A ⋅ cos(k ⋅ t ) )
Luego derivando y llevando a (1): W y = A ⋅ sen(P ⋅ t )
El vector velocidad
angular:
r r r r
W = Wx ⋅ i + W y ⋅ j + Wz ⋅ k
(
Llamando: A 2 = W x2 + W y2 )
r
Por tanto el módulo de W es:
r r
W 2 = W z2 + A 2
El vector velocidad angular
r
W efectúa un movimiento de
presesión alrededor del eje “z”
con pulsación:
(I − I )
k= 3 ⋅ Wz
I
Para el caso particular en que
r r
la condición principal es tal que L tiene la dirección del eje de rotación. I 3 ⋅ W z = L
W x = W y = 0 A = 0 ; no hay presesión, el sólido mantiene la dirección del eje “z”
constante.

102
 r
re
rT r
[ ]
r
M 0e = [I o ] ⋅ [α ] + W × [I o ] ⋅ W
T

r
∑ = M ⋅ acm
F
r r r
a cm = W 2 ⋅ rcm + α ⋅ rcm
El tensor de inercia respecto a un sistema
de ejes cuyo origen se encuentra en un
punto que no coincide con el centro de
masas, se puede calcular con la expresión
tensorial
[I o ] = [I G ] + M ⋅ [[r0G ]2 ⋅ [I D ] + r0G ⊗ r0G ]

103
Unidad 6: Sistemas vibratorios de un grado de libertad

Conceptos generales
El estudio de los sistemas vibratorios trata sobre el comportamiento oscilatorio de los
cuerpos, es decir, las oscilaciones alrededor de una posición de equilibrio.
Se estudiará el comportamiento dinámico de masas sometidas a fuerzas elásticas. El
análisis se efectúa sobre masas puntuales, es decir, sistemas materiales de dimensiones
despreciables. Posteriormente se tratarán sistemas con dimensiones considerables (cuerpos
sólidos) en los cuáles esas deben ser tenidas en cuenta a efectos del estudio de su
comportamiento dinámico.
Los sistemas vibratorios pueden clasificarse de acuerdo a distintos puntos de vista.
Pueden ser sistemas lineales o no lineales. Un sistema es lineal cuando la composición
de las respuestas a las fuerzas actuantes consideradas por separados es igual a la respuesta a
la composición de dichas fuerzas. Los sistemas reales son generalmente no lineales pero en
el caso de pequeñas oscilaciones pueden ser aproximados mediante sistemas lineales.
Los sistemas pueden ser de uno o varios grados de libertad, de acuerdo al número de
coordenadas independientes necesarias para su descripción.
En cuanto a la vibraciones, estas pueden ser libres o forzadas. Las vibraciones libres se
producen cuando se aparta al sistema de su posición de equilibrio y luego este oscila bajo
la acción de fuerzas propias del sistema. En el caso de vibraciones forzadas, el sistema
vibra sometido a la acción de fuerzas forzadoras externas.

Fuerza recuperadora elástica:

Es aquella dependiente de la posición (por lo tanto, admiten un potencial) y que se
materializa mediante elementos elásticos tales como barras, resortes, flejes, etc.
Entendemos por elástico a cualquier sistema susceptible de deformarse bajo
determinada solicitación y que recupera totalmente su forma o configuración original al
suprimirse la misma.
La fuerza recuperadora elástica es la respuesta del sistema a la interacción que lo
deforma.
Los cuerpos deformados almacenan energía potencial, denominada energía potencial
elástica. Esta energía es la necesaria para producir la deformación, y es devuelta al
recuperarse la forma original.

Modelación dinámica:
Las interacciones dinámicas entre fuerzas y masas se
estudian mediante modelos teóricos o abstracciones que,
suficientemente analizadas y próximas a la realidad,
representan satisfactoriamente a los sistemas reales.
Constituyen esquemas simplificados que permiten
progresivamente abordar el análisis de sistemas crecientes en
grado de complejidad. En este capítulo serán utilizados L
modelos físicos, representaciones esquemáticas del sistema
real, y modelos matemáticos, es decir, representaciones
analíticas mediante una ecuación diferencial o un sistema de
ecuaciones diferenciales establecidas en base a leyes físicas,
como la segunda ley de Newton. Estos modelos matemáticos
δ
describen el comportamiento dinámico del modelo físico en
cuestión.

105
Deformación de sistemas elásticos típicos y fuerzas recuperadoras elásticas

6.2.1. Se analiza en primer lugar el comportamiento de una barra elástica sometida a

tracción, de longitud l y sección transversal A.
Se asume el material cumple la Ley de Hooke, que relaciona la deformación resultante
(respuesta) con la fuerza actuante (solicitud o interacción).

δ = Deformación
F = Solicitación
L = Longitudid original (sin deformación) de la barra
A = Área sección transversal
E = Módulo de elasticidad (módulo de Young)

La relación entre la deformación (δ) y la solicitación (F) está dada por:

FL
δ= [6-1]
AE
L 1
Haciendo =
AE K

y eligiendo la coordenada x = δ, se puede escribir:

F = Kx [6-2]
relación que expresa que la fuerza es proporcional a la deformación o elongación. El
coeficiente K depende de las características geométricas (L, A) y elásticas (E) del
material que constituye la barra.
Desde el punto de vista energético, el trabajo (W) necesario para deformar la barra es
almacenado como Energía potencial elástica o de deformación:
x 1
W = ∫ − Kxdx = - Kx 2 = ∫x Kxdx =V(0)-V(x)
0
0 2
1
Si V(0) = 0 ⇒ V(x) = Kx 2 [6-3]
2
Viga elástica simplemente apoyada, sometida a una carga concentrada en el centro

Sean:
E = Módulo de elasticidad del material de la viga. F
J = Momento de inercia de la sección transversal
de la viga. x
L = Luz de la viga.
F = Carga o solicitación.
x = deformación (deflexión) de la
viga en el punto central
L/2 L/2
La flecha o deformación “x” en el centro de
acuerdo con la teoría de la resistencia de materiales vale:
FL3
x=
48EJ
48EJ
y haciendo =K [6-4]
L3

106
 48EJ
F(x) = Kx = x [6-5]
L3
24EJ
V(x) = (½)Kx2 = 3
x2
[6-6]
L
Barra elástica sometida a torsión (se aplica un momento torsor M en su extremo).

Sean:
G= Módulo de elasticidad transversal del material que
constituye la barra.
JP = Momento de Inercia Polar de la sección transversal de la
barra.
L = Longitud de la barra.

Al aplicar el momento torsor M la sección transversal extrema

experimenta una rotación relativa (deformación angular) φ que
responde a: M
ML M
φ= =
GJ P Kφ
GJ P
con Kφ = [6-7]
L
GJ
Así: M = Kφ φ = P φ [6-8]
L
GJ P 2
V(φ) = (1/2) Kφφ2 =(1/2) φ [6-9]
L

Modelación

A efectos de lograr una mejor interpretación conceptual del comportamiento de

sistemas dinámicos, los casos (sistemas) típicos examinados son susceptibles de ser
reducidos a modelos como los que se indican:

Barra sometida a tracción:

Se sustituye la barra elástica por un resorte axial cuya

constante elástica es K.
K

107
Barra sometida a flexión (viga con carga concentrada en el centro)

Si la carga está aplicada

en el centro de la viga, F
esta puede ser
reemplazada por un x
resorte cuya constante K
vale, según [6-4]:

48EJ
K= L/2 L/2
L3

Barra sometida a torsión.

Una barra cilíndrica sometida a torsión puede asimilarse a un resorte

torsional, cuya constante es, de acuerdo
con [6-7]:

GJ P
Kφ =
L
Esta constante relaciona el momento torsor
aplicado con el ángulo de deformación, sus
unidades son distintas de las constantes
anteriores.
M

Movimiento de una partícula en las inmediaciones de un punto de equilibrio

Según se ha visto en el capítulo de dinámica de la partícula, en presencia de fuerzas de
dependen de la posición únicamente, se define la función de energía potencial V(x).
Esta función V(x) puede desarrollarse en series de Taylor:
∞
1
V(x) = ∑ V (n) ( x 0 )( x - x 0 ) n [6-10]
n =0 n!
V(x) = V(x0)+V(1)(x0)(x-x0)+(1/2)V(2)(x0)(x-x0)2+....

Derivando respecto de “x”:

dV(x) dV(x 0 ) d 2 V(x 0 ) 1 d 3 V(x 0 )
= + 2
(x-x0)+ 3
(x-x0)2+…
dx dx dx 2 dx

Sabemos que si F = F(x), admite una función potencial tal que se verifica:
d2x dV(x)
m 2 = F(x) = -
dt dx
dV(x 0 )
y en un punto x = x0 donde se verifica equilibrio, se cumple F(x0) = - =0
dx
Si se adopta como origen del sistema referencial el punto donde se verifica equilibrio,
queda:

108
 •• d 2 V(x 0 ) 1 d 2 V(x 0 ) 2
m x = F(x) = - x- x - ...
dx 2 2 dx 2
••
m x = F(x) = -Kx – ax2 – b x3 - … [6-11]

dV(x 0 ) 1 d 2 V(x 0 )
donde: K= , a= , .....
dx 2 dx 2

Si no se desprecian los términos de orden superior, la solución de la ecuación

diferencial [6-11] proporciona la ecuación de movimiento del Oscilador Anarmónico o no
lineal.
Si, en cambio, se toma solamente el término lineal, despreciando los de orden superior,
queda la ecuación diferencial:
••
m x = -Kx [6-12]

Esta ecuación diferencial de segundo orden es la ecuación del movimiento del

Oscilador armónico o lineal. Debe quedar claro que el oscilador armónico es una
aproximación del caso real, inarmónico. Esta aproximación es buena cuando las
deformaciones son pequeñas.
La ecuación [6-12] expresa la dinámica de un sistema material constituido por una
partícula de masa m sobre la que actúa una fuerza recuperadora elástica (proveniente de un
“resorte” idealizado de constante elástica “K”).
Visto de otro modo, si se desprecian los términos de orden superior en el desarrollo en
serie de V(x), y adoptando x0 = 0, resulta:
V(x) = (1/2)Kx2 [6-13]
d 2 V(x 0 )
siendo K =
dx 2

K= Constante de proporcionalidad entre la fuerza “F” y la

m
deformación “x”. El modelo físico para este oscilador armónico
o lineal es el mostrado en la la figura. Si bien la masa se puede
dibujar como un bloque o cuadro, se interpretará se comporta
x
como una partícula.

Combinación de constantes elásticas

Un sistema mecánico puede estar constituido por

una serie de elementos elásticos (resortes) actuando K1
sobre la o las masas que lo constituyen. K2
Este conjunto de resortes que oponen resistencia a su
deformación, pueden ser reducidos o sustituidos por uno
solo de tal manera que el efecto dinámico producido
resulta equivalente. δ
Se distinguen dos modos de conexión de los resortes
componentes:

109
6.5.1. Los resortes experimentan la misma deformación, al aplicar la (o las) fuerza activa
sobre los mismos.

La masa dispone de un solo grado de libertad, y el único movimiento posible es la

traslación en dirección vertical.
“δ” es la deformación producida por la fuerza “F”.
La fuerza total “F” es absorbida por los resortes
de constantes elásticas “K1” y “K2”. La deformación δK2 δK1
que sufren ambos resortes es la misma, e igual a “δ”.

El diagrama de cuerpo libre para la masa “m” se

muestra en el dibujo de la izquierda.
La ecuación de equilibrio correspondiente es: F
F = δK1 + δK2
F = δ(K1 + K2) = δ KT
Llamando KT = K1 + K2, Constante total o equivalente.
Generalizando, para “n” resortes que sufren igual deformación:
KT = ∑ K i , i=1,2,...n [6-14]

Los resortes presentan deformaciones diferentes bajo la acción de una misma

solicitación

Según la disposición indicada, y por condición de

equilibrio se debe cumplir: K1
F1 = F2 = F (Las fuerzas recuperadoras de
cada resortes son iguales a la solicitación)
La deformación total del sistema será la suma de las δ1
deformaciones de cada resorte: K2
δT = δ1 + δ2
Haciendo δT = FT/KT = F1/K1 + F2/K2
= =F(1/K1+1/K2)
Simplificando: δ2
1/KT = (1/K1+1/K2)
En general, para “n” resortes:
F
1/KT = Σ 1/Ki [6-15]

Se puede realizar una analogía con las resistencias en un circuito eléctrico. Las fuerzas
serían los equivalentes a las tensiones, y las deformaciones equivalentes a corrientes
eléctricas. Los resortes que tienen la misma deformación (el primer caso expuesto) pueden
asimilarse a resistencias conectadas en serie (la misma corriente), en tanto que los resortes
sometidos a la misma fuerza equivalen a resistencias conectadas en paralelo (la misma
tensión aplicada).

Estudio de Vibraciones Mecánicas Armónicas

En lo que sigue, se considerará el Oscilador armónico o lineal. Se estudiarán

vibraciones de pequeña amplitud, que producen desplazamientos muy pequeños en torno

110
de una posición de equilibrio estable. Se analizarán en primer lugar, las vibraciones
armónicas libres, es decir las oscilaciones que se producen en ausencia de fuerzas externas.
Luego, serán estudiadas las vibraciones forzadas.

Vibraciones armónicas libres no amortiguadas.

Las vibraciones libres no amortiguadas son las que se producen en ausencia de fuerzas
disipativas, por lo que unas vez iniciadas, las oscilaciones se mantienen en forma
indefinida.
La ecuación del movimiento es:
•• F(x)
m x = F(x)= - Kx [6-16]
F(x) representa la fuerza recuperadora elástica. m
de [6-16]:
••
m x + Kx =0 x
k
llamando a ω o = ( ω o es la frecuencia
2

m
propia o natural del sistema, o frecuencia circular Posición de
de oscilación armónica. También se suele equilibrio
denominar pulsación).

••
x + ω o . x = 0 [6-17]
2
Asi
La [6-17] es una ecuación diferencial lineal homogénea de 2° orden.
Resolver esta ecuación implica hallar la ley del movimiento x(t), es decir la relación
entre la deformación “x” y el tiempo “t”.
Para hallar la solución se sigue el método clásico para resolver ecuaciones
diferenciales. Se propone como solución a x = ept ; donde p es una
constante.
• ••
Sus derivadas resultan: x = p . ept ; x = p2. ept
Sustituyendo en (1) queda: p2 + ω o = 0 − ωo
2 2
p=
p2 = ± i ω o
Entonces es válido x(t) = C1.ei ωo t + C2.e-i ωo t, que mediante un cambio de constante de
integración:
A A
C1 = .e iθ C 2 = .e −iθ
2 2

A i (ω ot +θ ) A −i (ω ot +θ )
X= .e + .e [6-18]
2 2
Utilizando las relaciones de De Moivre, queda:

X(t) = A cos (ωo t + θ) [6-19]

A, θ, son constantes a determinar una base a condiciones iniciales.

Así, si para t = to = 0 las condiciones iniciales son

•
x = xo; v = x = vo
111
 1/ 2
vo 2  v 2 
Será xo = A . cos θ xo +2
=A 2
→ A =  xo 2 + o 
ωo2  ωo2 
 
vo v
vo = - A ω sen θ y ω o .tgθ = → tgθ = − o
xo x oω o

De tal modo, la ley de movimiento [6-19] corresponde a un movimiento armónico de

amplitud A y desfasaje θ.
El ángulo de fase θ depende de las condiciones iniciales. En particular si vo = 0
entonces θ = 0; si xo= 0 y v0 ≠ 0, resulta θ = 90°.
2π
El movimiento es periódico de periodo: τ = [6-20]
ωo
1 ωo
La frecuencia del mismo es f = = . Así: ω o = 2π . f [6-21]
τ 2π
Nótese que la frecuencia natural es independiente de las condiciones iniciales, es decir,
que ante distintas configuraciones para las condiciones iniciales cambiarán la amplitud y el
ángulo de fase de la oscilación, pero este tendrá siempre la misma frecuencia igual a la
frecuencia natural.

Oscilador Amónico Torsional Libre

En este caso θ, es la coordenada del movimiento. La ecuación del movimiento es:
M=-kθ Aplicando la Ley de Newton
•• ••
M=I θ y I θ =-kθ donde I es el momento de inercia.

•• k
θ + .θ = 0 [6-22]
I
θ(t) = θo. cos (ωo t + β) Ley del movimiento

M: par aplicado
••
θ = α : aceleración angular
I: momento de inercia de la masa, respecto del polo O

Oscilaciones con influencia del peso

Un par de consecuencias importantes surgen de considerar el
modelo Masa-Resorte con influencia del peso, por ejemplo si el K
movimiento se realiza en dirección vertical.
En equilibrio estático, el resorte estará deformado por la
acción del peso. Luego resulta:
P = k . δest [6-23] (la fuerza del peso es igual a la δEst
recuperadora elástica). P x
Se introduce un sistema referencial cuyo origen coincide con la
posición de equilibrio, entonces el planteo de la 2° Ley de
Newton resulta:
••
m x = - k (δest + x) + P
112
 ••
m x = - k δest - K x + P y por [6-23]:
••
mx = -kx

Como conclusión
P mg
1) De [6-23] k = = relación que proporciona un modo expeditivo para
δ est δ est
determinar K (Habitualmente P es conocido y δest, se puede medir).

2) Teniendo en cuenta que la fuerza P predeforma al resorte en δest, ella no influye en las
características del movimiento.

Métodos energéticos

Para sistemas conservativos (fuerzas dependientes de la posición) vale:

E = T + V = constante T = energía cinética

V = energía potencial

Todo sistema de fuerzas recuperadoras elásticas del tipo

1 2
F=-kk admite V ( x) = kx (potencial elástico)
2
dV( x )
Tal que: − = F ( x) = − kx
dx

En un sistema mecánico con vibraciones libres, es decir, sin disipaciones debidas a

amortiguamientos, la energía es en parte cinética y en parte potencial.
La ecuación de movimiento surge de:

dE d
E = T + V = constante o = (T + V ) = 0
dt dt

Para un sencillo sistema masa-resorte

d  1 • 
2
dE d
m x + kx 2  = 0
1
= (T + V ) =
dt dt dt  2 2 
 
••
Que conduce a m x + k x =0

Si el sistema es más complejo, corresponde evaluar apropiadamente T y V.

dE
= 0 conduce a la ecuación de movimiento.
dt

113
Vibraciones Armónicas amortiguadas

Modelo dinámico
Al esquema se incorpora un elemento, el amortiguador, en el cual se origina una fuerza
dependiente de la velocidad, actuando de tal modo de oponerse a todo instante al sentido
de esta última. Esta fuerza se llama fuerza amortiguadora. En el caso más general es
función de una potencia de la velocidad
Fb = - k vn [6-25]

Para el estudio de oscilaciones armónicas, se adoptará n=1, que representa

razonablemente situaciones reales. La constante del amortiguador se identificará como “b”,
siendo sus unidades [Newton]/[m/s].
En el modelo dinámico, el amortiguador se simboliza como un cilindro-pistón.
Para obtener la ecuación de movimiento, se plantea el diagrama de cuerpo libre, sobre
la masa actúan dos fuerzas: la debida al resorte y la debida al amortiguador. Aplicando la
ley de
Newton, queda:
•• • m
mx =-kx–b x
kx •
•• •
mx + b x +kx=0 [6-27] b x

•• b • k k
que se puede escribir x+ x+ x = 0
m m b

k b
y llamando ω o 2 = y γ=
m 2m
la ecuación anterior queda:

•• •
x + 2γ x + ω o 2 .x = 0 [6-28]
Ecuación diferencial homogénea de las vibraciones libres amortiguadas.

Proponiendo como solución x(t) = ept con derivadas

• ••
x = p ept , x = p2 ept
Sustituyendo en [6-27]:
ept (m p2 + b p + k) = 0

Y para que ello se cumpla para todo t, debe ser :

2
mp + bp + k = 0 [6-29]
Esta ecuación es llamada ECUACIÓN CARACTERÍSTICA.
Las raíces de esta ecuación son:
2
b  b  k
p1, 2 =− ±   − [6-30]
2m  2m  m

O bien
p1, 2 = −γ ± γ 2 − ωo
2
[6-31]

114
 El análisis de [6-30] y [6-31] conduce a tres situaciones físicas diferenciadas, según el
radicando ( γ 2 − ω o ) negativo, positivo o nulo. Si el radicando es negativo, el caso se dice
2

es subamortiguado, y el movimiento puede describirse como oscilaciones amortiguadas, es

decir, que se producen oscilaciones cuya amplitud decae en el tiempo. Este hecho se debe a
la naturaleza disipativa de las fuerzas amortiguadoras. Cuando el radicando γ 2 − ωo es
2

positivo, se dice que el caso es sobreamortiguado, las fuerzas amortiguadoras son de mayor
importancia y no llegan a producirse oscilaciones. Cuando γ 2 − ωo es nulo, es caso se
2

denomina crítico, y marca el límite entre los casos anteriores.

Vibraciones Subamortiguadas
En este caso
(γ 2
)
− ω o 2 〈 0 , o también (γ2 < ωo2) [6-32]

Puede interpretarse como que el término en que interviene el amortiguamiento es de

“menor significación” en relación con ωo; la influencia del amortiguamiento tiende a
desaparecer cuando γ<<ωo y el movimiento se aproxima a las vibraciones libres no
amortiguadas.
Se define ω1 = γ 2 − ω o
2
[6-33]
cantidad que se denomina pseudopulsación.
Entonces:

γ 2 − ω o 2 = − 1. γ 2 − ωo 2 = iω1 donde i es la unidad imaginaria.

Las raíces pueden escribirse entonces como p1 = - γ + i ω1; p2 = - γ - i ω1

y la solución puede escribirse como la combinación lineal de las soluciones particulares:

x(t ) = C1 .e p1t + C 2 .e p2t = C1 .e (−γ +iω1 )t + C 2 .e (−γ −iω2 )t

Haciendo:
A A
C1 = e iθ = .(cos θ + isenθ )
2 2
A −iθ A
C 2 = e = .(cos θ − isenθ )
2 2
La ley del movimiento x(t) queda:

x(t ) = A.e −γt . cos(ω1t + θ ) [6-34]

Donde ω1 es la pseudopulsación o frecuencia circular de oscilaciones amortiguadas.

Las vibraciones subamortiguadas son efectivamente oscilatorias y de amplitud
decreciente en el tiempo. Su período de oscilación es:
1 2π o
τ1 = = [6-35]
f1 ω1

115
 Como se ve en la figura, las
oscilaciones son tales que su
amplitud va decayendo. El
Ae-γt
caso expuesto corresponde a
condiciones iniciales x(0)>0 y
•
x >0. La fuerza del
amortiguador disipa la energía
disponible, y el movimiento
cesa una vez que se disipa
toda la energía que
inicialmente se le comunica al
sistema.

Vibraciones Sobreamortiguadas (movimiento no oscilatorio)

En este caso:
(γ
− ω o 2 〉0
2
) [6-36]
Cuando se cumple esta condición, el término de amortiguamiento es importante en
relación a ωo, cuanto mayor sea γ frente a ωo la influencia del amortiguamiento será más
fuerte.
p1 = - γ + γ 12 − ωo 2
p2 = - γ - γ 12 − ωo 2

La solución es:
x (t ) = C1 .e − p1t + C 2 .e − p2t [6-37]
que representa un movimiento exponencialmente decreciente del tiempo, y aperiódico, es
decir que no se producen oscilaciones.

Vibración Críticamente Amortiguada

Ocurre cuando se cumple:
(γ 2
− ωo2 ) =0 o sea, γ 2 = ω o 2 [6-38]

Y la ecuación [6-30] conduce a una raíz doble p1 = p2 = - γ

− p1t − p2t −γ t −γ t
x (t ) = C1.e + C 2 .e = (C1 + C 2 ).e = C.e
O también, según se puede demostrar:
x (t ) = e −γ t .( A1 + A2 .t )
Con las constantes A1 y A2 dependientes de las condiciones iniciales:
to = 0
x = xo A1 = x o
•
x = vo A 2 = vo + γ x o

116
 La forma de la ley de X (t)
movimiento depende de los
valores de xo, vo.
Para un determinado xo se V(0) > 0
representan gráficamente distintas
evoluciones de la deformación, X0 V(0) = 0
variando el valor de la velocidad
inicial.
Tener velocidad inicial V(0)
significa que se impulsa a la masa t
en la misma dirección que la
deformación X(0), por lo que el V(0) <0
valor de la deformación primero
aumentará para luego decaer bajo la acción de la fuerza elástica del resorte. Velocidad
inicial negativa corresponde al caso en que la impulsa a la masa en dirección contraria a la
deformación, la cantidad de energía disponible (la suma de la almacenada en el resorte y la
cinética que se le comunica a la masa) puede ser suficiente para que la masa supere la
posición de equilibrio estático.

Factor o grado de amortiguamiento

De la condición de amortiguamiento crítico:
2
 b  k
γ 2
=  = ωo =
2
 2m  m

Llamando bc: coeficiente de amortiguamiento crítico.

bc = 2mω o = 2 mk [6-39]

El Grado de Amortiguamiento suele entonces expresarse como la razón:

b  b γ 
ζ=  ζ = = ⇒ γ = ζ .ω o  [6-40]
bc  2mω o ω o 

b < bc (ζ < 1) vibración subamortiguada

Si b = bc (ζ = 1) vibración crítica
b > bc (ζ > 1) vibración sobreamortiguada

Para Vibraciones subamortiguadas la ecuación [6-34] resulta:

x (t ) = A.e −ζω ot . cos(ω o (1 − ζ 2 )1 / 2 t + θ ) [6-41]

Donde (
ω1 = ω o 2 − γ 2 )1/ 2
(
= ω o 2 − ζ 2ω o 2 )
1/ 2
(
= ωo2 1 − ζ 2 )
1/ 2
[6-42]

Decremento Logarítmico
Cuando el factor de amortiguación ζ es menor que uno, se producen oscilaciones
amortiguadas, es decir, oscilaciones cuya amplitud decae en el tiempo. La razón entre dos
amplitudes máximas consecutivas es:
117
 2π
−γ 2π
x n +1 A.e . cos[ω 1 (t + τ ) + θ ] e .e
−γ (t +τ ) −γt ω1 −γ
ω1
= = =e
xn A.e −γ (t ) . cos(ω 1t + θ ) e −γt

2π
x n +1 − .γ
= e ω1 [6-42]
xn
2π
− .γ
ω1
Entonces: x n +1 = x n .e = x n .e −γτ1 [6-43]
2π
.γ
xn ω
y por lo tanto =e 1
x n +1

xn 2π
ln = .γ = δ (constante)= decremento logarítmico [6-44]
x n +1 ω1
El logaritmo de la razón entre una amplitud máxima y la inmediata siguiente es una
constante denominada decremento logarítmico δ.

2π 2π .τ 2π .ζ
δ = .γ = = [6-45]
ω1 ω o (1 − ζ ) 2 1/ 2
1−ζ 2

Cuando ζ→0 δ≅2πξ [6-46]

Si se puede medir dos amplitudes máximas consecutivas, en una vibración

amortiguada (se logra conocer δ) es posible obtener b, (despejando ζ en función de δ). Esto
indica una forma experimental de medir el coeficiente de amortiguamiento.

En cuanto a la Energía Potencial Elástica acumulada en sucesivas amplitudes máximas.

2  2π 
− 2 γ
E n+1 1 / 2.k .x n+12  x n+1   ω1 
= = = e [6-46]
En 1 / 2.k .x n 2 x 
 n 

La energía decrece al doble de rapidez que lo que decrece la amplitud (comparar con [6-
43]).

E n+1 = E n .e −2γτ1 [6-47]

Vibraciones Armónicas Forzadas

Ocurren vibraciones forzadas cuando la masa se encuentra sometida a la acción de una
fuerza forzadora o excitación. El análisis de diferencia según el cómo es la fuerza
forzadora.
Las soluciones más sencillas son las que se obtienen para fuerzas forzadoras
armónicas, es decir, cuando responden a la ley F(t) = F0 sen ωt. Para este tipo de
excitación, la deformación resultante será también armónica y de la misma frecuencia, pero
en general estará desfasada.

118
 Las fuerzas forzadoras periódicas pero no sinusoidales pueden descomponerse en
armónicas siguiendo el desarrollo en series trigonométricas de Fourier, luego la respuesta
será la composición de las respuestas a las excitaciones armónicas consideradas por
separado. Este procedimiento es aceptable si el comportamiento del sistema es lineal.

Vibraciones Forzadas Sin Amortiguamiento

Modelo Dinámico:
F(x) = -K.x
Sobre la partícula de masa m actúa
m la fuerza perturbadora exterior armónica:
F(t) = Fo . sen ω t
+x
Donde ω es la frecuencia circular de la fuerza perturbadora
(es frecuentemente producida por desbalances en maquinarias rotatorias).
En ausencia de amortiguación, la ecuación del movimiento resulta:
d 2x
m. 2 = −kx + Fo .senωt [6-48]
dt

k Fo
Llamando ω o 2 = y Po = (dimensiones de aceleración)
m m

Reemplazando en la ecuación queda

••
x + ω o . x = Po. sen ωt
2
[6-49]

Ecuación diferencial no homogénea de segundo orden de las vibraciones forzadas, sin

amortiguamiento.
La solución general, según la teoría de las ecuaciones diferenciales es:

x(t) = xh(t) + xp(t)

donde:
x(t) = ecuación del movimiento, da la posición en función del tiempo.
xh(t): solución de la ecuación diferencial homogénea asociada
xp(t): solución particular.

Según lo visto anteriormente: xh(t) = A cos (ωot + θ)

Para la solución particular, xp(t), se prueba con:
xp(t) = a senωt
derivando respecto del tiempo:
•
x p = a ω cos ωt
••
x p = - a ω2 sen ωt

Reemplazando en la ecuación [6-49], y para que se cumpla para todo t, debe ser:
Po
a= 2
[6-50]
ωo − ω 2

119
 ω
Llamaremos λ = (razón de frecuencias)
ωo
Po 1 Po Fo
a= . [6-51] = = δ est [6-52]
ωo 2
(1 − λ )
2
ωo 2
m
k
m

δest: deformación que acusa el resorte si Fo se aplica en forma estática, o también la

deformación que produce por acción del máximo valor que puede asumir F(t).

Entonces la ley de movimiento resulta:

x(t) = xh + xp = A. cos (ωo t +θ) + a . sen ωt [6-53]

Siempre que se cumpla el valor de a, según [6-51]

Puede escribirse:
a 1
=. [6-54]
a δ est 1 − λ2
δ est
Siendo δest una constante interesa
analizar la relación funcional entre a
ω
λ= y λ.
ωN
(para λ > 1, a < 0, pero se
acostumbra rebatir la gráfica sobre el
semiplano en que a > 0),

En la gráfica se observa que cuando

la relación entre la frecuencia
forzadora y la frecuencia natural es
igual a 1, la amplitud de la oscilación tiende a infinito, situación conocida como
ω
Resonancia: estado para el que λ = = 1 o sea ω = ω b , conduce a amplitud que tiende
ωb
a a →∝

Vibraciones Forzadas con Amortiguamiento

Modelo dinámico: En el esquema, se muestra un sistema de un grado de libertad,

cuyos elementos son la masa (m), el resorte (k) y el
F0Sen(ωt) amortiguador (b). También se muestran las fuerzas que
m
actúan sobre la masa: La fuerza forzadora F0senωt , la
kx • •
b x fuerza del resorte (kx) y la fuerza del amortiguador (b x ).
Considerando un desplazamiento “x” desde la posición de
equilibrio estático, positivo dirigido hacia arriba, se
k plantea la ecuación de movimiento:
b •• •
m x = - k x - b x + Fo senωt [6-55]

120
 k b Fo b
ωo2 = γ = Po = ζ = 〈1
m 2m m bc

•• •
x + 2 γ x + ωo2 x = Po senωt [6-56]
Solución general
x(t) = xh(t) + xp(t)
xh(t) = A e-γt cos (ω1 t +θ)

Para la solución particular xp, se prueba:

xp = a . sen (ω t - β) [6-57]
•
x p = a ω cos (ωt - β)
••
x p = - a ω2 sen (ωt - β)

Haciendo: ωt - β = ∅ o bien ωt = ∅ + β

Probando en [6-56]

a (- ω2 + ωo2) sen ∅ + 2 γ ω a cos ∅ = Po (cos β sen ∅ + sen β cos ∅)

Lo cual, para que se cumpla, requiere que:

a (ωo2 - ω2) = Po cos β [6-58]

2 γ ω a = Po sen β [6-59]

Que permiten obtener los valores de la amplitud a y el ángulo de fase β de [6-57]

Po 2γω
a= [6-59] tgβ = [6-60]
(ω o 2 − ω 2 ) 2 + 4γ 2ω 2 ωo2 − ω 2

Algunas conclusiones:
1. Para distintas relaciones ω/ωo se obtienen amplitudes de oscilaciones forzadas
diferentes.
2. Si ω = 0, o si ω << ωo, λ→0.
a
→ 1 y la amplitud de las oscilaciones forzadas tiende a δest.
δ est
3. Si ω >> ωo la fuerza excitatriz F(t) = Fo. sen ωt cambia de signo muy rápidamente, no
a
dando tiempo a la masa a reaccionar. →0 y a→0
δ est
4. Si ω ≅ ωo a →∝, situación de resonancia.

ω
El ángulo de desfasaje β puede graficarse en función de λ = . La siguiente gráfica
ωo
corresponde a un factor de amortiguación muy pequeño.

121
 β
180º

λ=1
λ

∀λ < 1 (ω < ωo) → β = 0 fuerza y desplazamiento tienen igual sentido.

∀λ > 1 (ω > ωo) → β = π fuerza y desplazamiento tienen sentido opuestos.

Aclaración: las vibraciones forzadas sin amortiguamiento comprenden un caso ideal, útil a
efectos de presentar el problema de excitaciones exteriores. Para amortiguación muy
débiles (ζ << 1).
De tal modo, la solución completa es de:

x(t) = A e-γt cos (ω1 t + θ) + a sen ((ω t - β) [6-61]

(1) (2)

El primer sumando (1) se amortigua relativamente pronto. Llamemos ts: tiempo de

establecimiento, de tal modo que, si por ejemplo consideramos que las vibraciones
transitorias se desprecian cuando su valor es el 1 % de la amplitud de oscilaciones
forzadas.

Ae −γt = 0,01.a
1 100 A
ts = . ln
γ a

Las vibraciones transitorias (1) dependen de las condiciones iniciales. Las vibraciones
forzadas (2) de las constantes físicas (k, b, m) y de la fuerza perturbadora (Po, ω).
Resulta de alto interés analítico las relaciones funcionales [6-59]y [6-60]
A tal fin se definen:
ω b
λ= h=
ωo ωo

Po Fo
= = δ est (def. estática)
ωo2 k
La [6-59] queda:
a 1
= [6-62]
δ est (1 − λ2 ) 2 + 4h 2 λ2

y de la [6-60]
2 hλ
tgβ = [6-63]
1 − λ2

122
 En la figura se muestran
a gráficas de la relación
δ est a
adimensional para
δ est
distintos valores de h.

Se observa que la máxima

amplitud de la oscilación
1 se produce a un valor
ligeramente distinto de
λ=1.
ω
λ =1 λ=
ω0

Si en [6-62] llamamos ε = λ2, la amplitud resulta máxima cuando el denominador es un

mínimo.
Asi
f(ε) = (1 - λ2)2 + 4 h2 λ2 = (1 - ε)2 + 4 h2 ε

tiene un mínimo f’(ε) = 0 = - 2 ε + 4 h2 → ε = λ2 = 1 – 2 h2

Llamando resonancia al valor λ que conduce a amax:

λres = (1 – 2 h2)1/2

Por consiguiente, en presencia de amortiguamiento, las vibraciones forzadas (ω) conducen

a resonancia, para valores λ un poco inferiores a 1 (es decir, para ω < ωo).
Algunas situaciones particulares:
a
a) Si λ << 1, (o sea, ω << ωo), resulta → 1 y a ≅ δest, en tanto β ≅ 0.
δ est
b) Si λ es muy grande (ω >> ωo) el valor de a se toma muy pequeño y β → 180°.
c) Si λ ≅ 1 se produce resonancia.

Para la relación [6-63], se grafica el ángulo de desfaseje para distintos valores de factor de
amortiguación:

ω
ω0

123
Conclusiones
1. La amplitud de las oscilaciones forzadas no dependen de las condiciones iniciales del
movimiento.
2. En presencia de resistencia (medio viscosos) las vibraciones forzadas no se amortiguan.
3. La frecuencia ω de las oscilaciones forzadas no depende de las constantes físicas del
sistema (como ωo y ω1, por ejemplo) y está determinada y coincide con la de la fuerza
perturbadora.
4. Aún para Fo pequeñas, se pueden producir amplitudes forzadas muy elevadas.
5. Para ω>>ωo se puede minimizar las amplitudes, incluso para valores elevados de Fo.

Representación fasorial

La representación fasorial de las ecuaciones de movimiento aportan elementos

analíticos adicionales. Se recuerda que un fasor es un vector rotante que se utiliza para
representar magnitudes que son funciones sinusoidales del tiempo. Se utilizan para
representar tensiones y corrientes alternas, y aquí serán utilizadas para representar las
fuerzas forzadoras, fuerzas elásticas, de amortiguación y de inercia. Estas últimas son
armónicas de la misma frecuencia que la fuerza forzadora, pues dependen de la
deformación y de sus derivadas ( la velocidad y la aceleración).
El módulo de cada fasor se hace igual a la amplitud de la magnitud que representa.
Así, si la deformación sigue la ley x(t) = X seno(ωt), el fasor tendrá módulo igual a X, y se
entiende es un vector que rota con velocidad angular ω. Se recuerda que e valor
instantáneo de la deformación se obtiene como la proyección sobre el eje vertical del fasor.

Vibración Armónica Libre

•• •
mx +k x =0

•
x = - A ωo sen (ωo t + θ)
••
x = - A ωo2 cos (ωo t + θ) ⇒ - ωo2 x

Sea x el vector desplazamiento.

•
La fase del vector velocidad x está adelantada en 90° resp. de x
••
La fase del vector aceleración x está adelantada 180° resp. de x

Vibración Forzada no Amortiguada

•• •
m x + k x = Fo sen ωt

Caso (I); ω < ωo, β = 0

Fo y x no tienen desfasaje
Caso (II), ω > ωo, β = 180º

124
Fo y x están desfasados en 180º

Vibración Forzada con Amortiguamiento

•• •
m x + 2γ x + k x = Fo sen ωt

Interesa visualizar tres casos:

Cuando: 1) ω << ωo
2) ω ≅ ωo
3) ω >> ωo
1) ω << ωo Mω X
2

bωX
Fo: fuerza excitatriz
kX: fuerza de recuperación elástica
b m X: fuerza de amortiguamiento F0
m ω2 X: fuerza de inercia β
kX

1. Las fuerzas de inercia, así como las de

amortiguamiento son pequeñas.
2. El ángulo β es muy pequeño
3. Fo tiende a equilibrar a la fuerza elástica.
4. La excitación exterior tiende a aparecer como aplicada estáticamente.
2
Mω X
2) ω = ωo bωX
En este caso, el ángulo de fase entre la
deformación y la fuerza forzadora es de 90º.
La fuerza de inercia está equilibrada por la
fuerza del resorte, y la fuerza excitatriz se F0
β=90º
emplea para vencer a la fuerza de kX
amortiguación.

3) ω >> ωo
En este caso, al ángulo de fase
entre la deformación y la fuerza 2
Mω X bωX
forzadora supera los 90º, y si la
frecuencia tiende a un valor muy alto,
el ángulo tiende a 180º, quedando la
deformación en oposición de fase con F0
β>90º
la fuerza forzadora. Si se mantiene kX
constante la amplitud de la fuerza
forzadora, al aumentar la frecuencia se
verá que las amplitudes de la deformación y la fuerza en el amortiguador se harán
pequeñas, y para altas frecuencias, la deformación será prácticamente nula.

125
Desbalance Rotatorio

El desbalance de masas en máquinas rotatorias es el origen de fuerzas excitadoras. Se

considera a continuación un sistema masa-resorte, restringido a moverse en el sentido
vertical y sometido a una fuerza forzadora originada por el desbalance de una máquina
rotatoria.
Sean:
M = masa total del sistema
m = masa excéntrica
e= excentricidad
ω= velocidad angular de rotación
x = desplazamiento de la masa no rotante (M –m)
x + e Sen ωt = desplazamiento vertical de m.

La ecuación diferencial del movimiento es la

siguiente:
•• d2 •
(M – m) x +m (x+e Senωt)= - kx - b x [6-64]
dt 2

Desarrollando la derivada, y simplificando la

ecuación del movimiento queda:
•• •
M x + b x + kx = meω2Sen(ωt)
Haciendo: F0 = meω2
queda la ecuación:
•• •
M x + b x + kx = F0Sen(ωt) [6-65]
Que resulta idéntica a la ecuación diferencial correspondiente al caso de oscilación
forzada por una fuerza sinusoidal. En el proceso de resolución de esta ecuación diferencial,
la amplitud de la fuerza F0 puede tomarse como constante.
La respuesta permanente se puede escribir entonces como:
x(t) = X Sen(ω t - φ) [6-66]
Siendo la amplitud X :
meω 2
X= [6-67]
(k − Mω 2 ) + (bω ) 2
bω
y el ángulo de fase φ = ArcTang [6-68]
k − Mω 2
Utilizando las siguientes definiciones:
ϖn = k = frecuencia natural de oscilación no amortiguada.
masa
bc= 2 k × masa = 2 masaϖn = amortiguamiento crítico
ζ = b/bc = Factor de amortiguamiento

puede escribirse la ecuación de la amplitud en forma no dimensional:

126
 2
 ω 
 
MX
=  ωn  [6-69]
me   ω 2   ω  2
1 −    + 2ς 
  ω n    ω n 
 
el desfasaje quedará:
 ω 
2ς  
 ω n 
tan(φ) = [6-70]
2
 ω 
1 −  
 ωn 
En la siguiente gráfica, se muestra la amplitud “X”y el ángulo de fase “φ” versus la
frecuencia “ω”, utilizándose el factor de amortiguación como parámetro.

Movimiento del soporte

En muchos
casos la excitación
del sistema no es x
debida a una m
fuerza que actúa
sobre la masa,
sino por el
movimiento del
soporte. En esta
sección, se
y • • k(x-y)
considerará como b  x − y 
excitación el  
movimiento del
127
soporte y se verá de qué manera plantear la ecuación diferencial que rige el movimiento.
Sea “x” la coordenada que indica la posición de la masa, medida desde la posición del
•
equilibrio estático. La coordenada “y” indica el movimiento del soporte, y es la velocidad,
••
y la aceleración . En el gráfico de la derecha se indican las fuerzas no balanceadas (el
peso no se considera puesto que es compensada por la fuerza del resorte kδest). La
deformación del resorte es (x-y), y la velocidad relativa entre las piezas del amortiguador
• •
es  x − y  .
 
Se plantea la ecuación diferencial que rigen el movimiento de la masa:
•• • •
m x = - k (x-y) – b  x − y  [6-71]
 
En esta ecuación diferencial, intervienen las variables “x” e “y”, y sus derivadas.
Haciendo la sustitución:
z = x –y
la ecuación del movimiento se transforma en:
•• • ••
m z + b z +k z = - m y [6-72]
Para el caso en que el movimiento del soporte sea sinusoidal, y = Y sen(ωt), siendo
“Y” la amplitud del movimiento del soporte y “ω” su frecuencia angular.
Derivando “y” dos veces, y reemplazando la ecuación diferencial queda:
•• •
m z + b z +k z = mω2Y sen (ωt) [6-73]

La forma de esta ecuación diferencial es similar a la correspondiente al desbalance

rotatorio, con “mω2Y” reemplazando a “meω2”. La solución de la ecuación diferencial
puede entonces escribirse como:
z = Z sen(ωt - φ)
mYω 2
Con Z= [6-74]
2 2
(k − mω ) + (cω )

Haciendo x = z + y se obtendrá el movimiento absoluto “x” de la masa. Por lo

tanto:
x(t) = z(t) + y(t) = Y sen (ωt) + Z sen(ωt - φ )
Siendo x(t) la suma de dos funciones sinusoidales de frecuencia angular “ω”, será
también una función sinusoidal de la misma frecuencia:
x(t) = X sen(ωt - ψ )

donde “X” es la amplitud del movimiento

Y
de la masa, y “ψ” el desfasaje entre el
movimiento de la masa y el movimiento
del soporte.
X

Para obtener el valor de “X” y “ψ”, Z

se recurre a la representación fasorial:
128
Estos fasores pueden ser representados mediante el uso de números complejos:
y = Y ejωt
z= Z ej(ωt - φ)
x = X e j(ωt - ψ)

x= Y ejωt + Z ej(ωt - φ)

Utilizando la expresión deducida [6-74] para “Z”, se llega a que:

k 2 + (ωc) 2
X= Y [6-75]
(k − mω 2 ) 2 + (cω ) 2
mcω 3
tan(ψ) = [6-76]
k (k − mω 2 ) + (ωc) 2
En la gráfica siguiente, se muestra el grafo de relación de amplitudes X/Y, y el ángulo de
desfasaje “ψ” versus la relación ω/ωn , siendo el factor de amortiguación el parámetro.

[6-75] y [6-76]

En el grafo se observa que:

• La relación X/Y vale 1 (lo que implica que X = Y, es decir que la amplitud del
movimiento de la masa es igual a la amplitud del movimiento del soporte),
cualquiera sea el factor de amortiguación para ω/ωn = 2 .
• Para valores de ω/ωn > 2 , X/Y es menor que la unidad, lo que significa que
la oscilación de la masa será atenuada por la combinación del resorte y el
amortiguador. Debe notarse que la atenuación será más efectiva cuanto menor sea
el factor de amortiguación para frecuencia ω > ωn.

129
Aislamiento Vibratorio

Entendemos por aislamiento vibratorio la atenuación de los efectos de las vibraciones

originadas en máquinas y motores. Puede tratarse de aislar fuerzas perturbadoras
producidas por desbalances rotatorios o movimientos del soporte. Se verá que ambos
problemas pueden tratarse de manera similar. La solución de este problema consistirá en el
montaje de la máquina sobre resortes, o combinación de resortes y amortiguadores.
A continuación se analizará el aislamiento de la fuerza perturbadora originada por el
desbalance rotatorio. En la figura se muestra el diagrama de fasores, indicándose la fuerza
transmitida a través del resorte y el amortiguador como FT.

F0Sen(ωt)
mω2X
bωX
FT

kX
F X

2
 bω 
FT = (kX ) + (bωX ) = kX 1 + 
2 2

 k 

X = Amplitud del movimiento originado por una fuerza forzadora F0 Sen (ω t)

Reemplazando X, se puede escribir la relación (FT/F0) como:
2
 bω 
2
 ω 
1+   1 +  2ς 
FT  k   ωn 
= = [6-77]
F0 2
 mω 2   bω 
2 2
  ω 2   ω  2
1 −  + 1 −    + 2ς
 k   k    ω n    ω n 


Esta relación es idéntica a la relación X/Y hallada para el caso del movimiento del
soporte, por lo que se concluye que el aislamiento de una máquina de las perturbaciones
del soporte es idéntico al problema de aislar fuerzas perturbadoras. Las razones (FT/F0) y
(X/Y) se conocen como Transmisibilidad de fuerza o desplazamiento.
Del análisis del gráfico que corresponde a la relación X/Y se deduce que:

• Para lograr una transmisibilidad menor que 1 (de esta manera se atenuaría la fuerza
transmitida) los resortes sustentadores deben ser tales que la frecuencia natural
debe ser baja comparada con la frecuencia de la perturbación.
• En la región (ω/ωn) > 2 , a menor factor de amortiguación menor es la
transmisibilidad. Como contrapartida, en la región ω/ωn < 2 esta relación se
invierte.

130
Respuesta a fuerzas periódicas

En el caso de fuerzas forzadoras periódicas no sinusoidales, aplicando el análisis de

Fourier pueden descomponerse en una serie de componentes armónicos (suma de senos y
cosenos) de distintas frecuencias.
Una fuerza periódica se caracteriza por su período, el cual se identifica con el símbolo
T, y se cumple que F(t)=Ff(t+T), es decir que el valor de la fuerza para un instante t se
repite T segundos después.
Si F(t) es una función de período T, puede ser representada por medio de una serie de
Fourier:
∞ ∞
F(t) = a0+ ∑a
1
n cos ω n t + ∑ bn senω n t
1
[6-78]

Los coeficientes a0, an y bn se determinan mediante las fórmulas de Euler.

Las frecuencias de las armónicas son múltiplos de la frecuencia de la armónica
fundamental, la cual es igual a la frecuencia de la fuerza forzadora, es decir, ωn=nω1.
La serie de senos y cosenos pueden ser combinados en una serie de solo senos, y la
fuerza forzadora se escribe entonces:
∞
F(t) = ∑ F sen(nω t − φ
1
n 1 n ) [6-79]

Deben tomarse los términos necesarios para una correcta aproximación. Considerando
que los elementos del sistema se comportan en forma lineal, es factible hallar la respuesta a
cada armónico por separado, luego componer estas respuestas para determinar la respuesta
del sistema a la fuerza excitadora periódica.

131
Unidad 7: Oscilaciones Mecánicas de Dos o Más Grados de Libertad

Introducción

En el capítulo anterior se ha visto que el número de grados de libertad de un sistema es

igual al número de coordenadas independientes necesarias para describir su movimiento.
Para introducir conceptos relacionados con estos sistemas de varios grados de libertad, se
iniciará con el estudio de sistemas de dos grados de libertad.
Como primera característica de estos sistemas, los sistemas de dos grados de libertad
poseen dos frecuencias naturales. En casos particulares (que corresponden a ciertas formas
de iniciar el movimiento, es decir una combinación particular de condiciones iniciales) la
vibración libre del sistema ocurre a una de estas frecuencias naturales, existiendo una
relación definida entre las amplitudes de las dos coordenadas y la configuración
correspondiente es un modo normal de vibración. El sistema de dos grados de libertad
tendrá entonces dos modos normales de vibración, cada uno de los cuales corresponde a
una de las frecuencias naturales. Cuando el movimiento se origine bajo condiciones
generales, las vibraciones libres de las masas ocurrirán como la superposición de los
modos normales de vibración (esto quiere decir que el movimiento de cada masa será la
suma o superposición de dos movimientos armónicos, cada uno de los cuales tendrá una
frecuencia igual a la frecuencia natural).
Si el sistema está sometido por una fuerza forzadora, la vibración ocurrirá a la
frecuencia de la fuerza forzadora y las amplitudes tenderán a un máximo cuando esa
frecuencia coincida con las frecuencias naturales.

Modos normales de Vibración

En esta sección se determinarán las configuraciones de los modos normales de

vibración. En primer lugar, se modelizará un sistema de dos grados de libertad. Tomando
como base el sistema mostrado en la figura x.x, donde se han despreciado las
amortiguaciones, las ecuaciones diferenciales correspondientes son:
••
m x 1 = k (x 1 − x 2 ) − kx 1 [7-1]
••
m x 2 = k(x 1 − x 2 ) − kx 2

Un modo normal de vibración es aquel en el cual cada masa experimenta un

movimiento armónico de la misma frecuencia y fase u oposición de fase, es decir, pasan
simultáneamente por la posición de equilibrio.
Para tal movimiento, y utilizando números complejos, las ecuaciones de las
deformaciones pueden escribirse de la siguiente manera:
x1 = A1e jωt [7-2]
x 2 = A 2 e jω t
En estas ecuaciones, A1, A2 son las amplitudes de las armónicas, y ω la frecuencia a la
que suceden.
Sustituyendo en las ecuaciones diferenciales, queda el siguiente sistema de ecuaciones
homogéneo:
(2k - ω 2 mA1 − kA 2 = 0 [7-3]
− kA1 + (2k - 2ω 2 m)A2 = 0

133
 Estas ecuaciones tienen soluciones A1, A2 distintas de la trivial, si el determinante
siguiente es nulo:
( 2k - ω 2 m) -k
=0 [7-4]
-k (2k - 2ω 2 m)
Haciendo ω2=λ, al resolver el determinante se llega a la ecuación característica,
polinomio de segundo grado:
2
 k 3 k 
λ − 3 λ +   = 0
2
[7-5]
 m 2m
Las raíces de esta ecuación son los cuadrados de las frecuencias naturales del sistema.
De la resolución de la ecuación característica:
k
λ1= λ1 = 0,634 [7-6]
m

k
λ 2= λ2 = 2,366 [7-8]
m

Sustituyendo estas frecuencias naturales en las ecuaciones, se calculan las relaciones

entre las amplitudes, (A1/A2)(i), donde “i” indica la frecuencia para la cual se calcula dicha
relación. Para ω1:
(1)
 A1  k
  = =0.731 [7-9]
 A2  (2k - ω12 m)
( 2)
 A1  k
  = =-2.73 [7-10]
 A2  (2k - ω 22 m)
Nota: Los valores numéricos calculados corresponden al ejemplo que se está
desarrollando.
Estas razones de amplitud son las formas modales que corresponden a cada modo
normal. La razón de amplitud para ω1, 0.731 es positivo, significando esto que para la
frecuencia ω1 las dos masas oscilan en fase, en tanto que para ω2 el valor –2.73 indica que
las dos masas oscilan en contrafase.
Estos modos normales de vibración se producen para condiciones iniciales particulares.
Si el movimiento se inicia con condiciones iniciales diferentes de las de los modos
normales, las oscilaciones contienen a los modos normales simultáneamente
(superposición de armónicas de distintas frecuencias). Para determinar las ecuaciones
correspondientes, será necesaria la resolución de un sistema de ecuaciones formadas por
las dos ecuaciones de movimiento y las relaciones de amplitudes. Así, si el movimiento se
genera apartando las masas de su posición de equilibrio mediante las deformaciones x1(0)
• •
y x2(0), y las velocidades iniciales x 1 (0) y x 2 (0) , se reemplazan estos valores en las
ecuaciones de las deformaciones y sus derivadas:
x1(t)= A1(1) sen(ω1t+θ1) + A1(2) sen(ω2t+ϕ2) [7-11]
x2(t)= A2(1) sen(ω1t+θ1) + A2(2) sen(ω2t+ϕ2)
A estas cuatro ecuaciones se agregan las relaciones de amplitudes. De esta manera, se
genera un sistema de 6 ecuaciones con 6 incógnitas, cuya resolución permite calcular las
amplitudes (A1(1), A1(2), A2(1), A2(2)) y los ángulos de fase (θ1, ϕ2).

134
Modelización de un sistema de varios grados de libertad

En este apartado se describe el proceso de modelización mediante el cual se obtiene el

sistema de ecuaciones diferenciales para un
sistema de varios grados de libertad. Como
primer caso, se considera un sistema de dos
k1 grados de libertad sin amortiguación, tal
k1x1
como se muestra en la figura.
x1 m1
k2(x1-x2) Las masas m1 y m2 están sujetas por los
k2 resortes de constantes k1, k2 y k3. Respecto de
k2(x1-x2) las posiciones de equilibrio estático, la masa
x2 m2 m1 está separada la cantidad x1, y la masa m2
k3x2 la cantidad x2. Si se hace x1 > x2, el resorte k1
estará elongado, en tanto que los resortes k2 y
k3
k3 estarán comprimidos.
Considerando que los pesos de las masas
están compensados por la deformación
estática, no son incluidas en los diagramas de
cuerpo libre. Para la masa m1, la ecuación de Newton queda:
m1 &x&1 =-k1x1-k2(x1-x2)
En tanto que para la masa m2 será:
m2 &x&2 = k2(x1-x2)-k3x2
Estas ecuaciones diferenciales no son independientes entre sí, puesto que las
coordenadas x1 y x2 aparecen en ambas. Agrupadas forman un sistema de ecuaciones
diferenciales, cuya resolución dará como resultado las deformaciones x1 y x2 como función
del tiempo. El sistema de ecuaciones resultante es el siguiente:
m1 &x&1 + (k1 + k 2 ) x1 − k 2 x2 = 0

m2 &x&2 − k 2 x1 + (k 2 + k 3 ) x2 = 0
Haciendo k1 + k 2 = k11 , − k 2 = k12 , k 2 + k 3 = k 22 , el sistema de ecuaciones puede
reescribirse como sigue:
m1 &x&1 + k11 x1 + k12 x2 = 0

m2 &x&2 + k12 x1 + k 22 x2 = 0
Forma que permitirá escribir en forma matricial el sistema de ecuaciones diferenciales.
En el caso de que se tomen en cuenta las amortiguaciones, se cambia el modelo físico
introduciendo los elementos amortiguadores. En el dibujo se han agregado los amortigua-
dores c1, c2, c3.
Las velocidades de las masas son x&1 y x& 2 , se considera que la velocidad de la masa m1
es mayor que la velocidad de la masa m2 ( x&1 > x& 2 ).
Las fuerzas originadas en los amortiguadores están representadas en los diagramas de
cuerpo libre.

En base a la segunda ley de Newton, se escriben las ecuaciones de movimiento de las

dos masas, formando un sistema de ecuaciones diferenciales. Para la masa m1:
m1 &x&1 = − k1 x1 − k 2 ( x1 − x 2 ) − c1 x&1 − c2 ( x&1 − x& 2 )
Que se reescribe:

135
 k1 c1 x&1
c1
k1x1
x1 m1
x&1 ( x&1 - x& 2 )c2 k2(x1-x2)
k2 c2
( x&1 - x& 2 )c2 k2(x1-x2)
x2 m2
c3 x& 2 k3x2
x& 2
k3 c3

m1&x&1 + (k1 + k 2 ) x1 − k 2 x2 + (c1 + c2 ) x&1 − c2 x& 2 = 0

Para la masa m2:
m2 &x&2 = k 2 ( x1 − x2 ) − k 3 x2 + c 2 ( x&1 − x& 2 ) − c3 x& 2
m2 &x&2 − k 2 x1 + (k 2 + k 3 ) x2 − c2 x&1 + (c2 + c3 ) x& 2 = 0
Se concluye con el sistema de ecuaciones diferenciales:
m1 &x&1 + (k1 + k 2 ) x1 − k 2 x2 + (c1 + c2 ) x&1 − c2 x& 2 = 0

m2 &x&2 − k 2 x1 + (k 2 + k 3 ) x2 − c2 x&1 + (c2 + c3 ) x& 2 = 0
Nótese que en este caso, al igual que en el anterior, si se eliminan el resorte k2 y
el amortiguador c2 desaparece el acoplamiento y resultan dos ecuaciones diferenciales
independientes.

Acoplamiento de coordenadas.
Las ecuaciones diferenciales que describen sistemas de varios grados de libertad se
llaman acopladas porque las distintas coordenadas aparecen en cada ecuación.
Para un caso general de dos GL, las ecuaciones tienen la forma:
•• ••
m11 x 1 +m12 x 2 +k11x1+k12x2=0 [7-12]
•• ••
m21 x 1 +m22 x 2 +k21x1+k22x2=0

Las ecuaciones diferenciales pueden escribirse en forma matricial:

 m11 m 12   x1   k11 k 12   x1  0
••

m   ••  +    =  [7-13]
 21 m 22   x 2  k 21 k 22   x 2  0
 m11 m 12   k11 k 12 
Donde [M]=   es la matriz de masas, [k]=   es la matriz de
m 21 m 22  k 21 k 22 
rigidez. Si la matriz de masas es no diagonal, existe acoplamiento dinámico o de masa, en
tanto que si la matriz de rigidez es diagonal, existe acoplamiento estático o de rigidez. La
matriz de masas está multiplicando al vector de aceleraciones, en tanto que la matriz de
rigidez al vector de deformaciones.
136
 Los acoplamientos dependen del sistema de coordenadas elegido. Un sistema de
coordenadas en el cual no aparece acoplamiento alguno se denomina de coordenadas
principales o coordenadas normales.
Si el sistema presenta amortiguaciones, aparecerá la matriz de amortiguamiento
c 11 c12 
[c]=  
c 21 c 22 
Las ecuaciones diferenciales de un sistema de dos grados de libertad con amortiguación
se escribirán de la siguiente manera:
 m11 m 12   x1  c 11 c12   x1   k11 k 12   x1  0
•• •

m   ••  +    • +    =   [7-14]
 21 m 22   x 2  c 21 c 22   x 2  k 21 k 22   x 2  0

Valores y vectores propios

En esta sección se generalizarán los conceptos de modos normales y frecuencias
naturales para sistemas de más de dos grados de libertad. Para calcular estas cantidades, se
utilizarán los autovalores y autovectores.
Para la vibración libre de un sistema de varios grados de libertad sin amortiguamiento
se aplica la ecuación matricial [xx], que se puede escribir en forma abreviada:
[M] x  +[k] {x} = {0}
••
[7-15]
 
 
Premultiplicando por la matriz inversa de la matriz de masas, queda:
[M]-1[M]=[I] [I]=Matriz identidad
-1
[M] [k]=[A] [A]=Matriz de sistema o matriz dinámica
Utilizando estas definiciones, se puede reescribir la ecuación matricial:
[I] •x• +[A] {x} = {0} [7-16]
 
Asumiendo que los desplazamientos son movimientos armónicos, la amplitud de cada
aceleración será igual a la amplitud de la deformación multiplicada por el cuadrado de la
frecuencia:
••
X=λ X
donde λ = ω 2
La ecuación matricial queda:
[A − λ [I] ] {x} = {0} [7-17]
Igualando a cero el determinante de la matriz se obtiene la ecuación característica del
sistema:
A − λ[ I] =0 [7-18]
Las raíces de la ecuación característica son los autovalores o valores propios, y las
frecuencias naturales del sistema se calculan mediante la relación:
λ =ω2 [7-19]
Sustituyendo los autovalores en la ecuación matricial [7-17], se calculan los
autovectores. Estos autovectores se caracterizan porque en los mismos se encuentran las
formas modales, las que se determinan mediante las relaciones entre los componentes del
autovector. Por ejemplo, para un sistema con dos grados de libertad, para la primer
frecuencia natural se tendrá el autovector v1: [v11,v12]. La forma modal de la primer
frecuencia natural (A1/A2)(1) será igual a la relación (v11/v12).

137
 A cada frecuencia natural le corresponde un vector propio, el cual da la configuración
de modo normal que corresponde a dicha frecuencia.

Vibración armónica forzada

En los sistemas de varios grados de libertad, al ser excitado el sistema con una fuerza
sinuosidad las masas oscilan con la misma frecuencia que la fuerza excitadora. Si se varía
la frecuencia de la fuerza forzadora, se verificará que el sistema entra en resonancia cuando
la frecuencia coincide con una de las frecuencias naturales del sistema, es decir que un
sistema de varios grados de libertad tiene tantas frecuencias de resonancia como grados
de libertad.
Se analizará un sistema de dos grados de libertad, tal como se muestra en la figura. La
fuerza forzadora F1senωt actúa sobra la masa m1. El sistema de ecuaciones diferenciales
que corresponde a este sistema es:
m1 0   x1   k11 k12   x1   F1 
••

 0 m   ••  + k   =   sen(ωt)
k 22   x   0 
[7-20]
 2   x   21
 2  2
Considerando que la solución es:
 x1   X 1 
  =   sen(ωt)
 x2   X 2 

Es decir, las oscilaciones de cada masa tendrán la misma frecuencia que la fuerza
forzadora.
Sustituyendo el vector de deformaciones y su derivada segunda en la ecuación queda:

m1 0   − ω 2 X 1   k11 k12   X 1   F1 

 0 m   2  sen(ωt)+ k    sen(ωt)=   sen(ωt) [7-21]
 2  − ω X 2   21 k 22   X 2  0
Cancelando sen(ωt) y reordenando:
(k11 − m1ω 2 ) k12   X 1   F1 

(k 22 − m2ω )  X 2   0 
2  
= [7-22]
 k 21
Resolviendo este sistema de ecuaciones, se obtienen las amplitudes de las
deformaciones, las cuales pueden escribirse en función de las frecuencias naturales ω1, ω2
del sistema:
(k 22 − m2ω 2 ) F0
m1m2 (ω12 − ω 2 )(ω 22 − ω 2 )
X1 = [7-23]

− k12 F0
m1m2 (ω1 − ω 2 )(ω 22 − ω 2 )
X2= 2
[7-24]

En estas ecuaciones de las amplitudes, puede observarse que al variar la frecuencia de la

fuerza excitadora, el denominador se anulará dos veces, cada vez que la frecuencia
coincida con las frecuencias naturales del sistema. Estas son las condiciones de resonancia
en el sistema ideal sin amortiguaciones, las amplitudes de las deformaciones tenderían a un
valor muy grande.

138
En la gráfica se han representado las amplitudes X1 (línea continua) y X2 (línea a trazos) en
forma adimensional: se observan las dos frecuencias de resonancia.

Amortiguador de vibraciones

En el caso analizado anteriormente, la fuerza forzadora actúa sobre una de las masas, y
(
en la ecuación de la amplitud de esta masa, en el numerador aparece el factor k 22 − m2ω 2 )
(Ec [7-23]). Este factor se anula para una combinación adecuada de k 22 , m2 y ω . Al ocurrir
esto, la amplitud de la masa m1 será nula, es decir que la masa sobre la cual actúa la fuerza
forzadora no oscilará. Esto se debe a que para esta combinación de valores, la deformación
de la masa será tal que la fuerza que el resorte de constante k 2 ejerce sobre la masa m1 es
de igual magnitud y opuesta a la fuerza forzadora.
Este hecho puede ser utilizado para eliminar las vibraciones de una masa mediante el
agregado de una segunda masa vinculada a la primera por medio de un resorte,
(
seleccionando el valor de la masa masa m2 y el resorte k 2 de manera tal que k 22 − m2ω 2 )
se anule.

139
Unidad 8: Mecánica analítica

Principio de los desplazamientos virtuales

Introducción
En este capítulo se trabajará con los trabajos virtuales, principio que permitirá
establecer las condiciones de equilibrio de cualquier sistema mecánico.
La particularidad de este principio es que en la aplicación de este método la acción de
las ligaduras no se calcula introduciendo reacciones desconocidas, sino considerando los
desplazamientos que pueden ser comunicados a los puntos del sistema, si este se saca de la
posición ocupada. Estos desplazamientos son los denominados desplazamientos virtuales.
Los desplazamientos virtuales de los puntos de un sistema deben satisfacer las
siguientes condiciones:

1. Deben ser infinitamente pequeños, puesto que si fueran finitos el sistema pasará a otra
posición donde las condiciones de equilibrio pueden ser diferentes.
2. Deben ser tales que respeten las ligaduras impuestas al sistema.

Se adoptará como definición para los desplazamientos virtuales la siguiente:

Desplazamiento virtual de un sistema es un conjunto de desplazamientos infinitamente
pequeños que son compatibles, en un momento dado, con todas las ligaduras impuestas a
este sistema.
El desplazamiento virtual de un punto de un sistema será representado por un vector
∆s, dirigido en la dirección del desplazamiento.
Las ligaduras impuestas a un sistema determina un número finito de desplazamientos
virtuales independientes, número que coincide con el número de grados de libertad del
sistema.
Los demás desplazamientos virtuales pueden obtenerse como la suma de los
desplazamientos virtuales independientes.

Trabajos virtuales
Se denomina trabajo virtual al trabajo elemental que puede realizar una fuerza aplicada
a un punto material durante un desplazamiento que coincida con el desplazamiento virtual
de dicho punto. Sobre todo sistema sujeto a ligaduras actúan fuerzas activas y reactivas. El
trabajo virtual de la fuerza activa FA será ∆WA=FA•∆s y el trabajo ∆WR=∆FR•∆s.
Si se considera que las ligaduras son ideales, la suma de los trabajos elementales
debido a las fuerzas reactivas será nula, es decir que:
Σ∆WR = 0
Si se elige un punto arbitrario Bk del sistema, si se llama FAk a la suma de todas las
fuerzas activas que actúan sobre el punto, y Nk a la suma de todas las fuerzas reactivas tan-
to internas como externas será:
FAk + Nk= 0
FAk= -Nk
Por lo tanto, para cualquier desplazamiento virtual del punto Bk la suma de los trabajos
virtuales debido a fuerzas activas y reactivas debe ser nula. Si por otra parte se consideran
ideales las ligaduras impuestas, el trabajo virtual de las fuerzas reactivas será nula, por lo
que se puede establecer la siguiente ecuación:
∆WAk=0
Considerando todos los puntos del sistema:
Σ∆WAk=0
141
 Condición que deben satisfacer los sistemas mecánicos con ligaduras ideales que se
encuentran en equilibrio. Por lo tanto, la condición necesaria y suficiente para que un
sistema mecánico con ligaduras ideales se encuentre en equilibrio es que la suma de los
trabajos virtuales sea nula. Analíticamente se expresa esto de la siguiente manera:
Σ ∆WAk=Σ FAk•∆sk= Σ (FAkx∆x+FAky∆y+FAkz∆z) =0
Donde ∆x, ∆y, ∆z son las componentes de ∆s.

El principio de los trabajos virtuales se puede utilizar para resolver sistemas estáticos.
En general, en estos problemas, dada una configuración de cargas que actúan sobre el
sistema, se desea conocer el valor de las reacciones de los vínculos o ligaduras. El
procedimiento a seguir se puede dividir en los siguientes pasos:
a. Representar las fuerzas activas que actúan sobre el sistema.
b. Asignar al sistema un desplazamiento virtual independiente.
c. Calcular los trabajos elementales de todas las fuerzas activas para el
desplazamiento dado.
d. Repetir el procedimiento para todos los grados de libertad del
sistema.

Ecuación general de la dinámica

Tal como se ha visto, el principio de los trabajos virtuales proporciona un método
general para resolver los problemas de la Estática. El Principio de D’Alembert permite
aplicar los métodos de la Estática para resolver problemas dinámicos. Por lo tanto,
aplicando ambos métodos simultáneamente, es posible formular un método general para
resolver problemas de Dinámica.
Dado un sistema de puntos materiales al cual se han impuesto ligaduras ideales. A las
fuerzas activas y reacciones de ligaduras, se agregan las fuerzas inerciales Fink= - mkak,
luego, de acuerdo con el Principio de D’Alembert, el sistema estará en equilibrio.
Aplicando a estas fuerzas el principio de los trabajos virtuales:
Σ ∆WAk+Σ∆WRk +Σ∆Wink=0
Pero siendo ideales las condiciones de ligadura, Σ∆Wink=0, por lo tanto será:
Σ ∆WAk+Σ∆WRk =0
Esta ecuación es la ecuación general de la dinámica y significa que en el movimiento
de un sistema sujeto a ligaduras ideales, en cada instante la suma de los trabajos
elementales de todas las fuerzas activas aplicadas y de todas las fuerzas de inercia, para
cualquier desplazamiento virtual, es nulo.

Condiciones de Equilibrio y ecuaciones del movimiento del sistema en coordenadas

generalizadas

Coordenadas generalizadas y velocidades generalizadas

Dado un sistema con ligaduras, el número de coordenadas independientes que
determinan la posición de un sistema es igual al número de grados de libertad del mismo.
Si bien en el desarrollo de los capítulos de cinemática y dinámica se han trabajado con
distintos tipos de coordenadas (cartesianas, cilíndricas, esféricas ) donde las coordenadas
son en general distancias medidas sobre rectas o ángulos, se pueden elegir parámetros de
cualquier dimensión y significado geométrico, como ser longitudes de arcos, áreas,
volúmenes, etc.
Los parámetros de cualquier dimensionalidad, independientes entre sí, en número igual
al número de grados de libertad del sistema y que determinan la posición del mismo, se
denominan las coordenadas generalizadas del sistema. Estas coordenadas generalizadas
142
serán identificadas con la letra q, y los desplazamientos virtuales de cada coordenada, con
∆q.
Como ejemplo, una partícula que puede moverse libremente en el espacio tiene tres
grados de libertad. Las coordenadas cartesianas x,y,z permiten localizar está partícula, al
igual que las coordenadas esféricas r,θ,λ, las cilíndricas r,z,λ o las coordenadas
generalizadas q1,q2,q3.
Las coordenadas cartesianas pueden expresarse como funciones de las coordenadas
generalizadas, por ejemplo para una partícula con tres grados de libertad:
x=x(q1,q2,q3)
y=y(q1,q2,q3)
z=z(q1,q2,q3)

Para un sistema descrito por las coordenadas generalizadas q1,…,qn se denomina

velocidad generalizada asociada a la coordenada qk a su derivada respecto del tiempo,
•
qk . La dimensionalidad de estas velocidades generalizadas depende de las dimensiones
de la coordenada generalizada, por ejemplo, en coordenadas esféricas, para la coordenada
θ, su velocidad generalizada será medida en radianes/seg.

Fuerzas Generalizadas

Si sobre un sistema de n partículas cuya posición está determinada por coordenadas

generalizadas y sobre ellas actúan fuerzas F1, F2,…,Fn. Si se aplica a la coordenada q1 el
desplazamiento virtual ∆q1 mientras que las demás coordenadas permanecen invariables.
Luego, cada vector posición rk tendrá una variación ∆rk, que puede escribirse en
función de las coordenadas generalizadas de la siguiente manera:
Siendo rk=rk(q1,q2,....,qs)
Considerando que únicamente varía la coordenada q1, aplicando el diferencial:
∂r
(∆rk)1= k ∆q1
∂q1
El subíndice 1 de ∆rk indica que se trata de la variación en el vector posición rk debido
al desplazamiento virtual ∆q1.
El trabajo virtual debido a este desplazamiento virtual se calcula como:
 ∂r 
∆W1= ∑ Fk • ∆rk = ∑ Fk •  k ∆q1
 ∂q1 
Donde el punto representa el producto escalar. Si se extrae ∆q1 como factor común
queda:
  ∂r  
∆W1=  ∑ Fk •  k  ∆q1
  ∂q1  
 ∂r 
El término ∑F •  k  se denomina fuerza generalizada correspondiente a la
 ∂q1 
k

coordenada q1.
El razonamiento hasta aquí planteado puede extenderse para todas las coordenadas, de
esa manera se definen las fuerzas generalizadas correspondientes a cada coordenada.
 ∂r 
Qi= ∑ Fk •  k 
 ∂qi 

143
Condiciones de equilibrio en Coordenadas generalizadas
La condición de equilibrio expresada en coordenadas generalizadas es la siguiente:
∑ Qi ∆qi = 0
Como los desplazamientos ∆qi son independientes entre sí, esta ecuación se cumple
únicamente si las fuerzas generalizadas Qi son todas nulas.
Por tal motivo, se puede establecer como condición de equilibrio para un sistema
mecánico que todas las fuerzas generalizadas que corresponden a las coordenadas elegidas
sean iguales a cero.

Ecuaciones de Lagrange
Para resolver un sistema en coordenadas generalizadas, un camino posible es plantear
las ecuaciones de Newton en coordenadas cartesianas para luego realizar una
transformación de coordenadas. Otra forma es establecer directamente las ecuaciones del
movimiento en coordenadas generalizadas. Este método se debe a Lagrange, y las
ecuaciones diferenciales del movimiento en coordenadas generalizadas se denominan
Ecuaciones de Lagrange.
Para establecer las ecuaciones del movimiento, se parte de la ecuación:
Σ ∆Wk +Σ∆Wink=0
Asumiendo que no todas las ligaduras son ideales, el primer término puede incluir
también fuerzas de rozamiento.
En coordenadas generalizadas:
Σ ∆Wk = ∑ Qi ∆qi
De la misma manera que se generaron las fuerzas generalizadas para las fuerzas ac-
∂r
tivas, se introducen las fuerzas de inercia generalizada Qi in = ∑ Fkin • k .
∂qi
Con estas definiciones, queda:
+Σ∆Wink = ∑ Qiin ∆qi
Reemplazando la ecuación general de la dinámica:
Σ ∆Wk +Σ∆Wink= ∑ Qi ∆qi + ∑ Qiin ∆qi = 0
Ecuación que puede reordenarse de la siguiente manera:
( )
∑ Qi + Qiin ∆qi = 0
( )
Siendo las ∆qi independientes, para que se cumpla esta ecuación debe ser Qi + Qiin = 0 .
Desarrollando esta ecuación para cada coordenada generalizada se llegará a un sistema
de ecuaciones que pueden aplicarse directamente para resolver problemas de dinámica. El
planteo de Lagrange se basa en expresar las fuerzas de inercia generalizadas en función de
la energía cinética del sistema.
∂r
En la ecuación Qiin = ∑ Fkin • k , las Fkin pueden escribirse en base a las aceleraciones
∂qi
ak:
dv
Fkin = − mk a k = − mk k
dt
reemplazando en la ecuación de las fuerzas inerciales generalizadas, queda:
dv ∂r
Qiin = ∑ − mk k • k
dt ∂qi

144
 Resta ahora dejar toda la ecuación en función de la velocidad vk. En primer lugar se
puede escribir:
dv k ∂rk d  ∂r  d  ∂r 
• =  vk • k  - vk •  k 
dt ∂qi dt  ∂qi  dt  ∂qi 
dqi •  ∂r 
Además = qi y  k  = v k
dt  ∂t 
Entonces, haciendo:
d  ∂rk  d  ∂rk  ∂vk
 =  =
dt  ∂q1  dq1  ∂t  ∂q1
 ∂r 
La derivada parcial  k  es el límite del cociente incremental ∆rk/∆q1, por lo tanto,
 ∂q1 
aplicando la regla de L’Hospital queda:
 ∂rk   ∂ rk   ∂vk 
•

  =  •  =  • 
 ∂q1   ∂ q1   ∂ q1 
Utilizando estas relaciones:
 
dv k ∂rk d  ∂T  ∂T
• = -
dt ∂qi dq1  ∂ q•  ∂q1
 1
1
donde T = ∑ mk vk2 es la energía cinética del sistema.
2
Finalmente planteando esta ecuación para cada coordenada generalizada:

Estas son las ecuaciones diferenciales del movimiento en coordenadas generaliza-

das o ecuaciones de Lagrange.

Dinámica impulsiva

Se denomina movimiento impulsivo el originado por fuerzas de gran magnitud que

actúan sobre tiempos muy cortos, siendo la integral de la fuerza en función del tiempo
finita, como ocurre en los choques, explosiones o el disparo de proyectiles.
Sea una fuerza que actúa durante el lapso de
tiempo ∆t, con magnitud 1∆t
F(t) El impulso de la fuerza está dado por:
t + ∆t
I = ∫t 0 F (t )dt
1 0
F=
∆t Si se considera que el impulso es unitario cuando:
t + ∆t
I = ∫t 0 F (t )dt =1
0
t0
Cuando el impulso es unitario y la duración es
t muy pequeña, es decir, se está en el límite cuando ∆t
tiende a cero, la fuerza se denomina función delta de
∆t
Dirac, con la siguiente definición:
 1
δ(t-t0) =  ∆t si t 0 < t < t 0 + ∆t
 0 ∀ t < t0 o t > t 0

145
 De acuerdo con la definición, una fuerza que sigue la función delta de Dirac se debe
entender como de magnitud infinita, durante un lapso de tiempo infinitesimal. En la
realidad, la duración de una fuerza es impulsiva no es infinitesimal, pero se puede tomar
esto como una aproximación aceptable. Esta aproximación introduce una simplificación
notable en la resolución de problemas, ya que se puede considerar que la partícula
permanece inmóvil durante la aplicación del impulso.
Cuando actúan fuerzas de tipo impulsivo o percusiones, el análisis del movimiento se
puede hacer siguiendo los siguientes puntos:
1. Actuando simultáneamente fuerzas “ordinarias” e impulsivas, el efecto de las
primeras se desprecia.
2. Las partículas se consideran inmóviles durante el lapso de actuación de las
fuerzas impulsivas.
De acuerdo al teorema del momento angular, el impulso de las fuerzas es igual al
cambio de la cantidad de movimiento. En base a los principios enunciados se desprecian
las fuerzas no impulsivas, quedando la expresión del momento angular:
t0 + ∆t t + ∆t
∫t0 F (t )dt = ∫ 0
t0
dp

I = p1 − p 0
Es decir, que el impulso de las percusiones es igual al incremento en la cantidad de
movimiento.
En cuanto al momento cinético, y considerando que durante la acción de las
percusiones las partículas se consideran inmóviles, se llega a que la suma de los momentos
de las fuerzas impulsivas es igual al incremento del momento cinético del sistema.
t0 + ∆t t + ∆t t + ∆t
∫t0 M (t )dt = ∫ 0
t0
dL + ∫ 0
t0
v 0 × p dt

Llamando N a la suma de los momentos de las fuerzas impulsivas, queda:

N = L1-L0

146
 Índice Alfabético
A
Aceleración en la rotación.......................................... 47 Dinámica de la Partícula. ........................................... 67
Aislamiento vibratorio ............................................. 130 Dinámica de los Sistemas de
Álgebra Tensorial......................................................... 9 Partículas ................................................................... 83
Amortiguador de vibraciones ................................... 139 Dinámica del Sólido Rígido. ...................................... 89
Ángulo de nutación .................................................... 50 Dinámica impulsiva ................................................. 145
Ángulo de precesión................................................... 49 Dinámica Relativa de la Partícula .............................. 80
Ángulo de rotación propia.......................................... 50 Dinámica relativa del sólido....................................... 96
Ángulos de Euler........................................................ 49
Axoides ...................................................................... 57
E
B Ecuación de la base .................................................... 61
Base de coordenadas generalizadas............................ 29 Ecuación de la rodante ............................................... 62
Ecuación del eje central ............................................. 56
C Ecuación general de la dinámica .............................. 142
Campo de velocidades ............................................... 47 Ecuaciones de Lagrange........................................... 144
Cantidad de movimiento .......................... 75, 89, 90, 96 Eje central del movimiento ........................................ 55
Centro instantáneo de rotación............................. 59, 60 Equilibrio dinámico.................................................. 100
Centro instantáneo de velocidad ................................ 63 Estado de aceleración en el movimiento plano........... 63
Cinemática de la Partícula.......................................... 19 Estados de Equilibrio ................................................. 72
Estudio de Vibraciones Mecánicas Armónicas ........ 110
Cinématica del Cuerpo Rígido. .................................. 41
Cinemática Relativa ................................................... 34 F
Combinación de constantes elásticas ....................... 109
Condición cinemática de Factor o grado de amortiguamiento.......................... 117
aceleraciones .............................................................. 43 Forma cuadrática asociada ......................................... 15
Condición cinemática de Fuerza dependiente de la Posición ............................. 70
velocidades................................................................. 42 Fuerza dependiente de la velocidad............................ 70
Condición de rigidez .................................................. 41 Fuerza dependiente del tiempo................................... 69
Conservación de la cantidad de Fuerza recuperadora elástica .................................... 105
movimiento ................................................................ 85 Fuerzas Generalizadas.............................................. 143
Conservación del momento G
cinético................................................................. 86, 93 Grados de libertad ...................................................... 48
Coordenadas cartesianas ............................................ 20 I
Coordenadas cilíndricas ....................................... 21, 29 Invariante escalar ....................................................... 52
Coordenadas del polo de Invariante vectorial..................................................... 52
aceleraciones .............................................................. 63 M
Coordenadas esféricas.......................................... 23, 31 Masa del sistema ........................................................ 83
Coordenadas generalizadas ...................................... 142 Masa del sólido rígido................................................ 89
Coordenadas intrínsecas....................................... 26, 32 Mecánica analítica.................................................... 141
Curvas de potencial.................................................... 71
Métodos energéticos................................................. 113
D Modelación dinámica ............................................... 105
Modos normales de vibración. ................................. 133
Decremento logarítmico........................................... 117 Momento cinético .............................. 76, 86, 90, 95, 97
Definición de energía potencial.................................. 70 Momento y productos de inercia ................................ 91
Deformación de sistemas elásticos........................... 106 Movimiento Curvilíneo de la
Desbalance rotatorio ................................................ 126 Partícula ..................................................................... 74
Descomposición impropia.......................................... 54 Movimiento de un sólido con un
Descomposición propia.............................................. 54 punto fijo.................................................................. 100
Diagonalización de un tensor ..................................... 13 Movimiento de un Sólido por Inercia....................... 101
Diagonalización de una matriz Movimiento del soporte ........................................... 127
simétrica real.............................................................. 14 Movimiento helicoidal tangente................................. 57
Diagonalización del tensor de Movimiento paralelo .................................................. 57
inercia ........................................................................ 94 Movimiento Plano...................................................... 58
Movimiento polar....................................................... 58

147
Movimiento polar del sólido ...................................... 48 T
Movimiento rototraslatorio ........................................ 51
Movimiento sobre una curva lisa ............................... 79 Tensor de inercia.................................................. 16, 17
Movimiento sobre una superficie lisa ........................ 80 Tensores simétricos y antisimétricos.......................... 13
Movimiento unidimensional de la Teorema de Euler - Chassles...................................... 51
partícula ..................................................................... 68 Teorema de los ejes paralelos..................................... 95
Movimientos del Sólido en el Teorema del momento cinético .................................. 76
Espacio....................................................................... 43 Trabajo y Energía....................................................... 87
Trabajos virtuales ..................................................... 141
O Transformación de coordenadas........................... 10, 17
Operaciones con tensores:.......................................... 13 Transformación de Coordenadas................................ 33
Operadores lineales .................................................... 16 Translación pura......................................................... 44
Trayectorias polares ................................................... 60
Oscilaciones con influencia del peso........................ 112
Oscilaciones de Dos o Más Grados de Libertad....... 133 V
Oscilador Amónico Torsional Libre......................... 112
Valores y vectores propios ....................................... 137
P Vector de Posición, Velocidad y
Polo de aceleraciones ................................................. 64 Aceleración ................................................................ 19
Polo de reducción....................................................... 51 Vector velocidad angular ........................................... 46
Principio de Acción y Reacción ................................. 67 Velocidad en la rotación............................................. 45
Principio de D´Alembert ............................................ 81 Vibración armónica forzada ..................................... 138
Principio de independencia de Vibración armónica Libre ........................................ 124
acción de fuerzas........................................................ 67 Vibración Críticamente Amortiguada ...................... 116
Principio de Inercia (1º Ley de Vibración Forzada con
Newton) ..................................................................... 67 Amortiguamiento ..................................................... 125
Principio de los desplazamientos Vibración forzada no amortiguada ........................... 124
virtuales.................................................................... 141 Vibraciones Armónicas amortiguadas...................... 114
Principio de masa ....................................................... 67 Vibraciones Armonicas Forzadas............................. 118
Propiedades de los tensores........................................ 11 Vibraciones armónicas libres no amortiguadas ........ 111
Vibraciones Forzadas con
R Amortiguamiento ..................................................... 120
Radio vector o Vector de Posición ............................. 28 Vibraciones Forzadas Sin
Reacciones de vínculos .............................................. 78 Amortiguamiento ..................................................... 119
Reacciones dinámicas de un sólido en rotación ......... 98 Vibraciones Sobreamortiguadas............................... 116
Representación fasorial ............................................ 124 Vibraciones Subamortiguadas.................................. 115
Rotación pura ............................................................. 45 Vínculos o Ligaduras. ................................................ 77
Vínculos Rugosos ...................................................... 79
S
Sistemas de coordenadas
generalizadas.............................................................. 28
Sistemas de Referencia en
Movimiento................................................................ 34
Sistemas en Rotacion ................................................. 36
Sistemas en Traslación (únicamente) ......................... 35
Sistemas en Traslación y Rotacion Simultanea.......... 40
Sistemas Referenciales............................................... 28
Sistemas Rígidos y Deformables................................ 34
Sistemas vibratorios de un grado
de libertad. ............................................................... 105

148
BIBLIOGRAFÍA BÁSICA

Boresi, A. y Schmidt, R. (2001). Ingenieria Mecánica, Dinámica. Thomson Learning.

Pytel, A. y Kiusalaas, J. (1999). Ingeniería Mecanica, Dinámica. Thomson Learning.

Golstein, H. (1979). Macánica Clásica. Editorial Aguilar.

Leon, J. (1979). Mecánica. Editorial Limusa.

Symon, K. (1977). Mecánica. Editorial Aguilar.

Targ, S.M. (1986). Curso breve de Mecánica teórica. Editorial Mir.

Longhini, P.(1951). Lecciones de Mecánica. Tomo 2. Editorial El Ateneo.

Hertig, R. (1978). Mecánica Teórica. Editorial El Ateneo.

Beer, F. (1993). Mecánica Vectorial para Ingenieros: Dinámica. Editorial McGraw-Hill.

Spiegel, M.(1991). Análisis Vectorial. Editorial McGraw-Hill.

McLean, N. (1980). Mecánica para Ingenieros. Editorial McGraw-Hill.

Timoshenko, S. , Young, D. (1971). Dinámica superior. Editorial Urmo.

Mersherski, I. (1974). Problemas de mecánica teórica. Editorial Mir.

Thomson, W. (1992). Teoría de las vibraciones. Prentice Hall.

BIBLIOGRAFÍA COMPLEMENTARIA

Macaulay, D., Ardley, N. (1999). Como funcionan las cosas. Ed. Atlántida.

Roederer, J. (2002). Mecánica Elemental, Eudeba.

Alonso M. Finn, E. s/f . Física. Volumen I. Fondo Educativo Americano.

Lea, S. Burke, J. (1999). Física, la naturaleza de las cosas. (1999). Tomson Editores.

Serway, R. y Jewett, J. (2005). Física para ciencias e ingenierías. Volumen I. Tomson

Editores.

Sears, F. Freedman, R., Young, H. (2005). Física Universitaria con fisica moderna.
Pearson Educación.

Tipler, P. (1994). Fisica. (Tomo1). Ed Reverte.

149
Halliday, D. Resnick, R. Kaen, K. (2003). Física. Volumen I. Compañía Editorial
Continental.

Bueche, F. (1998). Física general. Mc Graw Hill.

Tipens, P. (1999). Física básica. Mc Graw Hill.

Hewitt, P. (2004). Física conceptual. Pearson Educación.

150