Mecanica Racional
Mecanica Racional
Mecanica Racional
Virginia Azurmendi
MECÁNICA
RACIONAL
MECÁNICA RACIONAL
edUTecNe
Buenos Aires, 2020
Ercoli, Liberto y Arzumendi, Virginia .
Mecánica racional / Liberto Ercoli ; Virginia Arzumendi. - 1a ed . - Ciudad Autónoma de
Buenos Aires : edUTecNe, 2020.
Libro digital, PDF
Reservados todos los derechos. No se permite la reproducción total o parcial de esta obra, ni su incorporación a un
sistema informático, ni su transmisión en cualquier forma o por cualquier medio (electrónico, mecánico, fotocopia, grabación u
otros) sin autorización previa y por escrito de los titulares del copyright. La infracción de dichos derechos puede constituir un
delito contra la propiedad intelectual.
El Agradecimiento
a nuestras familias, por su generosa cesión de
tiempo que pudo haber sido compartido.
Contenido
Capítulo 1: Cinemática del Punto
Video Introductorio
Presentación del libro por los autores, Dr. Ing. berto Ercoli e Ing. Virginia Azurmendi
http://www.youtube.com/embed/CzXyQSxWWOc
http://www.edutecne.utn.edu.ar/mec_racional/mec_racion-anexos.html
http://www.edutecne.utn.edu.ar/mec_racional/mec_racion-cv.html
PRÓLOGO
Liberto Ercoli
Virginia Azurmendi
10
Capítulo 1- Cinemática del Punto
La Cinemática consiste en estudiar el movimiento sin hacer referencia a las fuerzas que
lo originan analizando las sucesivas posiciones del punto en el espacio en función del
tiempo.
El movimiento de los puntos puede describirse especificando las coordenadas lineales
y angulares y sus derivadas respecto a marcos de referencia que pueden ser
considerados “fijos” (análisis del movimiento absoluto), o “móviles” (análisis del
movimiento relativo).
Los sistemas de coordenadas se eligen arbitrariamente en el marco de referencia
adoptado en función de la geometría del problema en cuestión. Las coordenadas más
corrientes son: cartesianas, polares, cilíndricas, esféricas, intrínsecas y generalizadas.
Como convención, se acuerda que en el caso de dibujar solamente dos ejes, el tercero
sale de la hoja hacia el lector.
11
Mecánica Racional
La base está constituida por el versor radial ê que tiene sentido positivo alejándose
G
del polo en la dirección del vector posición e (con origen en el polo y extremo en el
punto), y por el versor transversal Tˆ que es perpendicular al ê con sentido positivo en
la dirección creciente del ángulo T .
Como se verá más adelante, todos los sistemas coordenados son relacionables entre sí.
Trayectoria: es el lugar geométrico de las posiciones ocupadas por un punto móvil con
el transcurso del tiempo.
12
Capítulo 1- Cinemática del Punto
l = l (t) es la ecuación del movimiento en forma escalar. Para el instante inicial (t = 0), el
móvil ocupará en general una posición dada por el arco lo (espacio inicial) respecto al
origen de la posición. l1 es una posición de alejamiento máximo desde el punto tomado
como referencia.
En t2 el móvil está en la posición donde se comenzaron a medir los tiempos (l2 = l0). En
t3 el móvil pasa por el punto tomado como referencia para medir las posiciones.
La función l = l (t) es continua por cuanto un punto no puede ocupar más de una
posición para un instante dado. Esta forma de dar la posición se usa cuando se conoce
la trayectoria.
G
r (t ) x (t )iˆ y (t ) ˆj z (t )kˆ
Las proyecciones del móvil sobre cada eje coordenado constituyen los movimientos
proyectados:
G G G
rx x (t ) i ; ry y (t ) j ; rz z (t ) k
13
Mecánica Racional
G
³t 0
t t
dl dx 2 dy 2 dz 2 { dr ; l lo dx 2 dy 2 dz 2
en: t t1 o l1 l t1
en: t t 2 o l2 l t1 ' t
'l l2 l1 l t 1 't l t 1
14
Capítulo 1- Cinemática del Punto
'l
A la relación se la denomina expresión escalar de la velocidad media
't
'l
Vm
't
Luego, esta velocidad es una magnitud escalar y es la rapidez con que se recorren los
espacios en el tiempo. Cuando esta expresión se mantiene constante cualquiera sea t y
' t, el movimiento es uniforme. De igual forma:
G
'r G
Vm expresión vectorial de la velocidad media
't
Para intervalos de tiempo muy pequeños (' t o 0)
'l dl
lím V expresión escalar de la velocidad instantánea.
't o 0 't dt
G G
' r dr G
y lím V expresión vectorial de la velocidad instantánea. Veamos
't o 0 ' t dt
si ambas expresiones se relacionan entre sí. Siendo r una función del tiempo a través
de l, se tiene:
G G
r r > l (t )@
G G G
dr d r dl dr
V
dt dl dt dl
G G
pero en el límite d r es tangente a la curva y dl es el módulo de d r (y un vector
sobre su módulo es un versor):
G
En el límite, es d r { dl
G
dr
Luego, l versor tangente a la
d1
trayectoria
G
y V V l (1.1)
15
Mecánica Racional
El versor normal n que se elegirá está dirigido de P a O' y su recta de acción coincide
con la normal principal a la curva en P.
El plano que contiene a ly n se llama "plano osculador", el cual se mantiene fijo si
el movimiento es plano, por lo que es en este tipo de movimiento donde estas
coordenadas tienen su mayor aplicación.
16
Capítulo 1- Cinemática del Punto
b l n
G
Sea r ( t ) x (t ) i y (t ) j z (t ) k
Luego
G G
G. dx dy dz G dr
V r i j k , r
dt dt dt dt
G
V Vx i Vy j Vz k
ó G G G G
V Vx Vy Vz
2 2 2
dl § dx · § dy · § dz ·
dl dx dy dz
2 2 2
y V ¨ ¸ ¨ ¸ ¨ ¸ (1.1’)
dt © dt ¹ © dt ¹ © dt ¹
dT
donde T y T sen T i cosT j
dt
17
Mecánica Racional
G de dz
V eˆ e T Tˆ kˆ (1.2)
dt dt
ó
G
V Ve eˆ VT Tˆ Vz kˆ
§ de · 2 2 § dz ·
2
El valor del módulo será: ¨ ¸ e T ¨ ¸
(1.2’)
© dt ¹ © dt ¹
Cuando 't o 0:
'T dT
lím Z (1.3)
't o 0 't dt
18
Capítulo 1- Cinemática del Punto
G 1G G
'S r 'r
2
1 1
S rV sen\ rVT
2 2
Por lo tanto:
1 2
S r T (1.4’)
2
19
Mecánica Racional
G
'V G
A la relación am
't
se la denomina expresión vectorial de la aceleración media.
Tomando límite:
G G G
'V dV G d 2r
lím a (1.5)
't o 0 't dt dt 2
G dl
V V l l
dt
G
dV G dV d l
luego: a l V
dt dt dt
dV d 2l G dV ˆ
at ; at l
dt dt 2 dt
20
Capítulo 1- Cinemática del Punto
dl
Para calcular l nótese que cuando la partícula se mueve a lo largo del arco dl en
dt
el tiempo dt, l conserva su magnitud unitaria cambiando sin embargo su dirección, de
modo que se vuelva l c .
Si se grafica a continuación el cambio (variación) dl de l :
Aquí dl se extiende entre dos puntos que están sobre un arco infinitesimal dT de
radio l 1
Por lo tanto, dl tiene un módulo dl l dT y su dirección se define por n
.Consecuentemente:
dl . .
dl dT n ó l T n
dt
.
dT l V
y como: dl U dT o T
dt U U
.
V
por lo tanto lˆ nˆ
U
G
Por lo tanto, la expresión final del vector a cuando se lo refiere a la terna intrínseca es
G dV V 2
a l n (1.6)
dt U
G dV ˆ G V2
Con at l y an nˆ
dt U
G
donde an es la aceleración normal o centrípeta.
21
Mecánica Racional
Gráficamente:
Conclusiones:
G
a) El vector a no tiene componente según la binormal y por lo tanto está contenido
en el plano osculador en el punto considerado.
G
b) La at tendrá siempre la dirección de la velocidad y define la variación del módulo de
la misma. G
G
c) La at puede tener el mismo sentido o contrario que la V de donde resultan los
movimientos acelerados y desacelerados respectivamente.
G
d) La an está siempre orientada hacia el centro de curvatura y define el cambio de
dirección del vector velocidad a lo largo de la trayectoria.
22
Capítulo 1- Cinemática del Punto
Detalle en
P
G
r P O x (t ) i y (t ) j z (t ) k
G dr (t )
V Vx i Vy j Vz k y ?
dt
G G 2 2 2
G dV dr 2 d x d y d z
a 2
i 2
j 2
k
dt dt 2 dt dt dt
G
luego: a a xi a y j az k (1.7)
y su módulo a a x2 a 2y a z2
23
Mecánica Racional
G
Sea r (t ) e e z k
G G
dr de de dz
V e e k
dt dt dt dt
de dT dT
dt dt
senT i cosT j dt
T
con T A e
G de dT dz
Luego: V e e T k
dt dt dt
Y derivando:
G
G dV d 2e de de de dT d 2T dT dT d 2 z
a
e T e 2 T e k
dt dt 2 dt dt dt dt dt dt dt dt 2
Pero:
dT dT dT
dt
d
dt
senT iˆ cos T ˆj dt
cos T iˆ senT ˆj
dt
eˆ
Luego:
de dT de dT d 2T § dT ·
2
G d 2e d 2z
a e T T e 2 T e ¨ ¸ e 2 k
dt 2 dt dt dt dt dt © dt ¹ dt
24
Capítulo 1- Cinemática del Punto
Y agrupando:
G ª d 2 e § dT · 2 º ª de dT d 2T º ˆ d 2 z ˆ
a « 2 e¨ ¸ » «
ˆ
e 2 e »T 2 k (1.8)
«¬ dt © dt ¹ »¼ ¬ dt dt dt 2 ¼ dt
G G G G
ó a ae aT a z
G G
donde: ae es la aceleración radial y aT la transversal
G
az es la aceleración en el sentido de la cota y coincide en coordenadas
cartesianas o cilíndricas.
Este gráfico muestra las componentes de un mismo vector aceleración en los distintos
sistemas de coordenadas. Obsérvese que las componentes varían de un sistema a otro,
G
pero el vector a es único.
25
Mecánica Racional
'Z
A la relación: Jm
't
se la denomina aceleración
angular media en el lapso
't.
'Z dZ d 2T
lím J (1.9)
't o 0 't dt dt 2
G. G.
Sean S 1 y S las velocidades areolares en los instantes (t + 't) y t respectivamente. El
vector:
.
G. G. G.
'S S 1 S
define la aceleración areolar media.
't 't
y:
G. G. G
'S dS d 2S G
..
lím S (1.10)
't o0 't dt dt 2
G. 1G G
Conociendo que S r V y derivando:
2
..G
1 G G 1G G ..G
S V V r a o 1G G
2 2 S r a (1.10’)
2
26
Capítulo 1- Cinemática del Punto
(1.12)
V3
U G G (1.13)
V a
Con esta expresión se puede hallar el radio de curvatura de una curva dada con
independencia del movimiento. En efecto, suponiendo un movimiento según una de
las coordenadas y conociendo la curva, se tiene la otra coordenada; luego se hallan
G G
V y a y se aplica la última expresión.
El mismo procedimiento puede emplearse para pasar de coordenadas cartesianas a
polares o de polares a intrínsecas.
Para una mejor comprensión de las expresiones de Darboux se sugiere realizar la
gráfica de la trayectoria y en el punto dibujar los vectores posición, velocidad y
aceleración con sus respectivas componentes.
27
Mecánica Racional
1.6. Diagramas:
Es común en la mecánica representar gráficamente algunos de
los parámetros estudiados. Los diagramas usados más comúnmente son:
dl
Como V , la pendiente de la tangente a la curva l = l(t) en cualquier instante
dt
representa el módulo de la velocidad instantánea: V = tg M
b) Diagrama Velocidad-Tiempo:
aquí es V = V (t)
dl
pero V
dt
? dl = V dt
t
l l1 ³ V dt
t1
28
Capítulo 1- Cinemática del Punto
c) Curva Hodógrafa:
Sea una trayectoria cualquiera sobre la que se dibujan los vectores velocidad;
Es decir que la velocidad del punto P’ que recorre la hodógrafa es la aceleración del
punto P que describe la trayectoria.
d) Diagrama de aceleración-tiempo:
aquí es at = at (t); y dv = at dt
t
V V1 ³t1
at dt
29
Mecánica Racional
1 ª ciclos º
f « »
en 1 seg. Æ T ¬ seg . ¼ o [Hertz]
30
Capítulo 1- Cinemática del Punto
Para la aceleración:
G ª .. . 2º
ª . . ..
º ˆ
a « r r T » eˆ « 2 r T r T » T
¬ ¼ ¬ ¼
31
Mecánica Racional
G ..
a r Z 2 e r T T
Luego:
G
a J r l Z 2 r n (1.15)
at J r ; an Z 2 r
..
T J o dZ J dt y Z Zo J (t t o )
t2
si to 0 o es: Z J t Zo y T J Zo t T o
2
1 Z
y la frecuencia: f
T 2S
@KZ@Z0 Æ
T Zt T o (1.16)
32
Capítulo 1- Cinemática del Punto
Para P’ es:
O = centro de oscilación
x = elongación o posición
r = semiamplitud
2r = amplitud
2S
Z = pulsación o frec. circular =
T
T = fase o ángulo de fase
y: xo = r cos To
.
luego: V x r Z sen (Z t T o ) (1.18)
..
a x r . Z 2 . cos Z t T o (1.19)
33
Mecánica Racional
Puede demostrarse que este movimiento es periódico, para lo cual basta con evaluar x,
V y a, a partir de un instante cualquiera tomando tiempos iguales al período del
movimiento:
x r cos T o
° o
t = 0 ® Vo Z r sen T o
°
¯ ao Z 2 . r . cosT o
xt
r cos 2S T
o
rcosT o xo
2S °°
t=T= ®Vt rZ sen2S T o rZ senT o Vo
Z °
°¯ at rZ cos2S T o rZ cosT o
2 2
a0
2S
Así, el período de este movimiento es : T
Z
1) Elongación nula (x = 0)
x
r cos Z t T o
? cos Z t + T o 0
S
? Z t To 2. n 1 n 0,1,2....
2
el instante de nulidad, será:
* 2n 1 S / 2 T o
t
Z
** n S To
t
Z
34
Capítulo 1- Cinemática del Punto
a
2
Z r cos Z t T o
? cos Z t T o 0 ? Z t T o 2. n 1 S 2
* 2 n 1 S 2 T o
t
Z
4) Aceleración máxima:
.
a amáx Z 2 ro a 0
. n S To
Z r sen Zt T o ? t
3 **
a
Z
Dado que el eje de abscisas representa la fase del movimiento, el diagrama (V-T) no
constituye la gráfica (V-t), aunque T = T (t). Para ello, debe reemplazarse T por t en las
escalas, luego:
35
Mecánica Racional
T Z t To
Para t = 0 (origen de tiempos) será T = To por lo tanto, puede ubicarse el origen de los
tiempos para el valor To en los gráficos anteriores. Con respecto a la escala, se está
condicionado por aquella que se adoptó para los T. Veamos por ejemplo: si se hubiese
2 S rad T seg
tomado Escala T correspondería tomar Escala t .
10 cm 10 cm
En otras palabras, al segmento que representa T 2 S rad , le corresponde en el
2S
gráfico de tiempos, el mismo segmento pero representando T seg.
Z
que es la ecuación diferencial del MOA (segundo orden, primer grado, homogénea)
2 2
O Z 0
de donde: O1,2 riZ
Luego, existen dos valores de O que hacen que la función propuesta sea solución de la
(1.20). La solución general será una combinación lineal de ambas:
A + B = C1 ; A - B = - i C2
C1 C2 C1 C2
con lo que se obtiene: A i ; B i
2 2 2 2
36
Capítulo 1- Cinemática del Punto
Resulta importante pasar esta función a una forma equivalente; sean dos nuevas
constantes X y \ tal que:
C1 X cos \
C2 X sen \
y reemplazando en (1.21):
? x (t ) X cos (Z t \ ) (1.22)
con
§ C ·
2 2 2
X C1 C2 , \ arc tg ¨¨ ¸¸
© C1 ¹
Con esta última forma de la expresión de x(t) se observa la coincidencia entre las
soluciones aportadas por el razonamiento gráfico-cinemático (1.17) y el analítico
{| !
!}+
!
$
$+!{|`~!
*
!{Z0).
Estas constantes deben determinarse de las condiciones iniciales del problema:
x ( o) x o
°
t 0®.
°̄ x ( o ) V
o
.
? x ( t ) X Z sen Z t \
y xo X cos\ , Vo X Z sen\
de donde:
§ Vo · 2
X x o2 ¨ ¸
©Z ¹
§ Vo ·
\ arc tg ¨ ¸
© Z xo ¹
Ejemplo de aplicación:
Sea el movimiento del punto P que describe la circunferencia de
radio r = 0.27 m, con una velocidad angular de 20 rad/seg. Encontrar el movimiento de
P’ sobre la dirección x si en t = 0 era xo = 0.1m, vo = 5m / seg.
x(t ) X cosZ t \
.
x Z X senZt \
37
Mecánica Racional
xo 0.1m X cos\
1
Vo 5m / s 20 X sen\
seg
Y dividiendo miembro a miembro
50 20tg\
tg\ 5 / 2
\ 1.19rad .
Ordenando:
0.1m X cos\
5/ 20m X sen\
2 2
X 0.1 (5 / 20) 0.27 m
Respuesta:
x(t) = 0.27 cos (20 t - 1.19) m
x1 r1 cos T 1
x2 r2 cos T 2
38
Capítulo 1- Cinemática del Punto
T1 Zt T o1
Donde:
T2 Zt T o 2
? T1 T 2 T o1 T o 2 D cte.
x x1 x 2 r cos T
r cos Zt T o
Para resolver el problema se debe calcular su semiamplitud r y su ángulo de fase inicial
To, ya que Z es la misma por girar todos los vectores con igual velocidad angular.
Para ello se proyectan los vectores posición sobre ambos ejes en t = 0
2 2 ry
r rx r y To arctg
rx
y el movimiento resultante será:
x
r cos Zt T o
Como se observa, el movimiento resultante es otro MOA con la misma pulsación que
los que lo generan, semiamplitud r y ángulo de fase inicial To.
2S
El período es T
Z
39
Mecánica Racional
x1 r1 cos T 1 r1 cos Z1t T o1
x2 r2 cos T 2
r2 cos Z 2 t T o 2
Z1t o T o1 Z2 t o T o2 To
T o 2 T o1
y ? to
Z1 Z 2
G G
Siendo to el tiempo necesario para que r2 alcance a r1 . Con este to en las anteriores
se calcula To. En estas condiciones se puede simplificar el problema ya que,
comenzando a contar los tiempos a partir de to, los dos movimientos componentes
tendrían la misma fase inicial:
x1
r1 cos Z1t T o ½°
¾ x x1 x 2
x2
r2 cos Z 2 t T o °¿
El problema podría simplificarse aún más girando el eje x un ángulo To, con lo que se
anularía la fase inicial, en tal caso:
x1c r1 cos Z1 t
x 2c r2 cos Z 2 t
40
Capítulo 1- Cinemática del Punto
p Z1
Número Racional = T2/T1
q Z2
Luego para que el movimiento resultante sea periódico la relación de pulsaciones tiene
que dar un número racional (cociente entre números enteros).
El período resultante T será el tiempo mínimo, único y común para ambos
movimientos, es decir, el mínimo común múltiplo.
T = p T1 = q T2
Ejemplo:
Sean los movimientos: x1 4 cos 4S t S 3
x2
6 cos 3S t S 6
41
Mecánica Racional
Z1 4S 4
? Número Racional
Z2 3S 3
Z1 p 4
Como: p 4 y q 3
Z2 q 3
2S 1 2S 2 1 2
Así : T1 ; T2 y T 4 3 2 s
Z1 2 Z2 3 2 3
Gráficamente:
42
Capítulo 1- Cinemática del Punto
Teniendo en cuenta que las pulsaciones son iguales, la diferencia de fase entre ambos
movimientos será constante, como se vio en el caso (a.1).
T 2 T1 T O 2 T O1 cte.
43
Mecánica Racional
y
cos Zt cos T o 2 sen Zt sen T o 2
r2
y x
sen T o1 sen T o 2
r2 r1
sen T o1 T o 2
y x
cos T o1 cos T o 2
r2 r1
sen Zt
sen T o1 T o 2
elevando al cuadrado y sumando:
2 2
x y x y 2
2
2 2 cos T o1 T o 2 sen T o1 T o 2 0
r1 r2 r1 r2
44
Capítulo 1- Cinemática del Punto
Casos particulares
a) T o1 T o 2 0 rG1A rG2
? cos T o1 T o 2 1
sen T o1 T o 2 0
Fig a)
luego:
2 2
x y 2 xy
2
2 0 que es la ecuación de una elipse degenerada
r1 r2 r1 r2
§ x y ·2
¨¨ ¸¸ 0 dos rectas coincidentes
© r1 r2 ¹
r2
de donde y x
r1
b) T o1 T o 2 S 2 rG1 { rG2
aquí: cos T o1 T o 2 0
sen T o1 T o 2 1
2 2
x y
? 2
2
1
r1 r2
45
Mecánica Racional
G G cosS 1
c) T o1 T o 2 S r1A r2 ®
¯sen S 0
2 2
x y 2x y
2
2
0
r1 r2 r1 r2
2
§x y· r2
ó ¨ ¸ 0 ? y=- x
© r1 r2 ¹ r1
3
d) T o1 T o 2 S (r1 { r2 pero en sentidos opuestos)
2
3
° cos S 0
° 2
®
° 3
°¯ sen 2 S 1
e) Para fases intermedias entre las analizadas más arriba, los movimientos resultantes
responderán a trayectorias elípticas más o menos achatadas e inclinadas según se
acerquen a la recta o a la elipse.
46
Capítulo 1- Cinemática del Punto
En la siguiente figura, se observa que para cada diferencia de fase resultará una
trayectoria distinta, dando lugar a una familia de curvas que reciben el nombre de
“figuras de LISSAJOUS”.
47
Mecánica Racional
x r1 cos Z1t T o1
y r2 cos Z2 t T o 2
p Z1
Número racional
q Z2
Siendo en tal caso el período resultante el mínimo común múltiplo de dos períodos
componentes:
T = p T1 = q T2
Para cada relación de pulsaciones habrá una trayectoria distinta, dando lugar a figuras
de Lissajous como las mostradas, donde se ha tomado r1 2
48
Capítulo 1- Cinemática del Punto
*
+!
k !
+
!+$$
*
+
!
!%
!`
+!
k
Estos movimientos pueden ser convergentes o divergentes según que hacia el centro C
concurran las aceleraciones o sus prolongaciones.
G
Si el movimiento es referido al centro C, tomado como polo, el vector a sólo tendrá
componente radial.
G
a a ae e e T 2
En el movimiento central, la
aceleración areolar será siempre
nula, por cuanto:
G
1G G
S ea
2 G G
Como los vectores e y a son
colineales, resulta:
G
G G G
0 S vector cte.
dS
S
dt
49
Mecánica Racional
G 1G G G
Siendo S e V , el vector V y el punto C determinan un plano invariable en el
2
que se desarrolla la trayectoria.
Por tanto, todos los movimientos centrales tienen trayectorias planas que se barren
con
G G
S = constante y S 0. La recíproca también se cumple.
e e(e, T )
°
Datos ® G 1 2
°̄ S 2 e T cte. (*)
e e T
2
y sea: a e a e (t ) (**)
2S O
T con O 2S Constante de Binet
e2 e2
es e e > T ( t )@
de dT
luego: e
dT dt
O de 2 de d (1 / e )
Así: e Oe O
e
2
dT dT dT
2 2 2
d ª d 1/ e º d 1/ e O d 1/ e
Ahora: e « O » O 2 T 2
dt ¬ dT ¼
2
dT e dT
reemplazando en (**):
50
Capítulo 1- Cinemática del Punto
O ª d (1 / e ) 1º
2 2 2 2 2
O d 1/ e O
ae e 2« » (1.23)
e «¬ dT e »¼
2 2 4 2
e dT e
§ d1 / e ·
2
O2 ª
2 § d1 / e ·
2
1º
V O ¨
2 2
¸ 2 O «¨ ¸ 2»
© dT ¹ e ¬«© dT ¹ e ¼» (1.24)
Movimiento de Kepler:
Uno de los movimientos centrales más interesantes es el llamado kepleriano que
cumplen los planetas alrededor del sol.
El problema consiste en determinar la trayectoria de un punto material atraído por un
centro fijo en razón inversa al cuadrado de la distancia, es decir:
G k
ae eˆ (1.25)
e2
51
Mecánica Racional
T To 0
G G
V Vo
Datos t = 0 M Mo
G G
e eo
G k
a eo e
eo2
1º
O2 ª d 2 1 / e k
« »
e 2 «¬ dT 2 e» e2
¼
d 2 1/ e 1 k
ó (1.26)
dT 2 e O2
1
[ (T )
e
1 k
[ AcosT E
(T ) e O2
O2 / k
ó e
AO2
(1.27)
1 cos T E
k
52
Capítulo 1- Cinemática del Punto
Esta es la ecuación de una cónica en coordenadas polares referida a un foco y cuyo eje
focal forma un ángulo E con el eje polar.
Donde:
O2 b2
P parámetro de la cónica.
k a
2
AO c
H excentricidad de la cónica
k a
Recordando que:
Æ Circunferencia
Æ Elipse
Si
Æ Parábola
Æ Hipérbola
Cálculo de O:
1
O 2S 2 e 2T eeNT e V cosM
2 VT
pero siendo S = cte., su valor se mantendrá para toda posición del movimiento. En
particular para t = 0, se tendrá:
O eo Vo cos M o
Cálculo de A y E:
Se necesitan dos ecuaciones; la (1.27) es una. Para la otra, se
expresa (1.24) en función de las incógnitas:
1 k
se tiene: A cos T E 2
e O
53
Mecánica Racional
2
1 k 2 kA
A cos T E 4 2 cos T E
2 2
luego: 2
e O O
d 1/ e
y también: A sen T E
dT
reemplazando en (1.24):
k2
V2 O2 A 2 2 Ak cosT E (1.28)
O2
O2 / k
° eo {|
° AO2
® 1 cosE
° k
°¯Vo2 O2 A 2 k 2 / O2 2 Ak cos E {|
1 k
De (1.27’) A cos E 2 (1.27”)
eo O
en (1.28’)
1/ 2
2 k2 § O2 · k ½
°Vo 2 2k ¨¨ 1¸¸ 2 °
° O © keo ¹O °
A ® ¾
° O 2
°
° °
¯ ¿
§ 1 / e k / O2 ·
o
Luego: E arccos ¨ ¸
© A ¹
Velocidades orbitales:
O2 / k
e (1.27. bis)
1 H cosT
54
Capítulo 1- Cinemática del Punto
Si se elige la velocidad V0 +
& !
k%
! !
@ K !
anterior expresión resulta:
O2
e = cte
k
Y como es O eo Vo cos M o +
0 = 0 resulta O eo Vo , de donde
k
Vo2
eo
por debajo de esta velocidad inicial el vehículo no entrará en órbita, sino que caerá
describiendo una especie de espiral. Entre este valor de la velocidad y el
correspondiente a la parábola, !
`
!+ @ !
k
(1.27.bis) se tiene
O2
e = cte
2k
Y procediendo de la misma manera que antes, resulta
2k
Vo2
eo
Para este valor de V0 la trayectoria será abierta y el vehículo no regresará, motivo por
!
!!
$
#!
+!
`
+%!
para la cual un foco actúa como centro de atracción y el otro como centro de
repulsión, por ello para una aceleración del tipo
G k
a 2 e
e
sólo puede usarse la rama que hace que el movimiento sea de atracción.
Ejemplo de Aplicación:
55
Mecánica Racional
Solución:
S
2
O / 2 29075 Km / s
AO2
b) H Debemos determinar A. Del sistema (1.27’ – 1.28’), y teniendo en cuenta
k
que para este sistema coordenado es E = 0, resulta:
1 k 1 39,88 x10 4 Km 3 / s 2
A 2
eo O 6978Km 58150 2 Km 4 / s 2
H
0,253 x 10
4
58150 Km / s
2
2 2
0,215
4 3 2
39 ,88 x 10 Km / s
H 1 elipse
c) V 2 O2 A 2 k 2 / O2 2 Ak cosT
39,88 x 104 2
2
Vapg 58150 0,253 x 10
2 4 2
581502
2 0,253 x 104 39,88 x 104 cos S
2,1644 47,03 20,18 29,019 Km2 / s2 .
56
Capítulo 2 – Cinemática de los Sistemas de Puntos
Cuando un conjunto de partículas se mueven ligadas por ciertas relaciones se dice que
constituyen un sistema. Éste puede ser discreto o continuo según esté constituido por
una cantidad finita o infinita de puntos. Cuando se trata de determinar la posición de
todos los puntos de un sistema, en apariencia se necesitarían tantos parámetros como
puntos materiales tenga el sistema por tres Pi (xi, yi, zi), pero como los puntos
materiales están ligados por ciertas relaciones entre sus coordenadas (o parámetros),
dicho número será mucho menor.
2.1. Definiciones:
Ejemplo: sean tres puntos que forman un sistema material en el cual las distancias
entre ellos deberán permanecer constantes. Se tiene:
57
Mecánica Racional
1º) x1 x2 x3 0
d 1 2 cte.
2º) d 2 3 cte.
d o1 cte.
y1 , z1
°
m = 6 ® y2 , z2 (2.1)
°y , z
¯ 3 3
z 2 z1 2 y 2 y1 2 d122
°°
n = 3 ® z 3 z 2 y 3 y 2
2 2
d 232 (2.2)
° 2
¯° z1 y1 d o1
2 2
Esto significa que en (2.1) pueden tomarse a voluntad tres coordenadas como
libres, pero las otras tres quedan dependientes de éstas. Es decir, si se da una variación
(movimiento del sistema) a las libres, las dependientes se moverán obligatoriamente
de determinada manera compatible con las condiciones de vínculo (2.2) impuestas.
a) m = qj = T1 ,T 2 ,T 3 ( j 1,2,3)
b) n = 0
c) k = 3 - 0 = 3 gl
58
Capítulo 2 – Cinemática de los Sistemas de Puntos
Se observa que los grados de libertad son independientes del sistema de coordenadas
utilizado.
Cuando las restricciones pueden expresarse en formas algebraicas sencillas del tipo de
las anteriores, se denominan holónomas, y en ellas las ecuaciones de ligadura o de
vínculo pueden expresarse en función de las coordenadas solamente, o de las
coordenadas y el tiempo si los marcos de referencia son móviles; en éstas, el número
de coordenadas generalizadas es igual al número de grados de libertad. Existe un
grupo de problemas para los que las restricciones resultan ser ecuaciones diferenciales
no integrables que se denominan no holónomas, los cuales no serán tratados en este
curso.
En los sistemas móviles siempre resulta m > n dado que los gl representan los distintos
grados de movilidad del sistema por ser el número de parámetros libres que pueden
variarse en forma arbitraria. Se los denomina sistemas móviles o mecanismos. Cuando
m = n, el sistema material queda inmóvil por cuanto no habrá parámetros libres y las
coordenadas del sistema resultarán de resolver el sistema de ecuaciones formado por
las condiciones de vínculo. Son los llamados sistemas isostáticos. Existen sistemas
materiales donde m < n, en este caso tampoco habrá parámetros libres y el sistema
seguirá siendo inmóvil. Son los sistemas hiperestáticos. Desde el punto de vista
cinemático los dos últimos no tienen aplicación. Este curso se orientará hacia el
estudio de los primeros.
rG
j
G 2
ri
G
rij2 d ij2 cte. (2.3)
59
Mecánica Racional
rG1 rG2 2 2
d12 cte.
rG1 rG3 2 2
d13 cte.
rG2 rG3 2 2
d 23 cte.
x x 2 y y 2 z z 2 d 2
° 1 2 1 2 1 2 12
°
n = 3 ® x1 x 3 y1 y 3 z1 z3
2 2 2 2
d13
°
°¯ x 2 x 3 y 2 y 3 z2 z3
2 2 2 2
d 23
x1 , y1 , z1
°
Conjunto de coordenadas m = 9 ® x2 , y2 , z2
°x , y , z
¯ 3 3 3
rG rG
2
j i d ij2 cte.
60
Capítulo 2 – Cinemática de los Sistemas de Puntos
ordenando:
G G
rG rG V rG rG V
j i j j i i
G G
Dividiendo m.a m. por d ij r j ri
G G
V j rˆij Vi rˆij (2.4)
G G
r j ri
con: rˆi j G G
r j ri
a.1) Traslación
°a) Estados simples ®
° ¯a.2) Rotación
Movimientos de los ° b.1) Suma de Traslaciones
Sistemas Materiales ® °
Rígidos ° ° b.2) Suma de Rotaciones
b) Estados compuestos ®
°
° b.3) Suma de Traslaciones +
° °
¯ ¯ Suma de Rotaciones
61
Mecánica Racional
G G G
rj ri rij = vector constante (2.5)
Es decir que un segmento definido por dos puntos de un cuerpo permanece paralelo a
sí mismo durante el movimiento de traslación.
1) Las trayectorias de todos los puntos del sistema son congruentes, es decir, son
idénticas pero se dan en distintos lugares del espacio, son superponibles:
de (2.3)
G G G G
r j ri r jc ri c. = vector cte.
G G G G
? ri c ri rjc rj
G G
ó ' ri ' rj
si G G
't o 0 'ri { dri
G G
'rj { drj
G G y
dri drj
62
Capítulo 2 – Cinemática de los Sistemas de Puntos
Luego, si todos los puntos tienen la misma velocidad, ésta estará representada por un
G
vector libre llamado vector traslación W (t) cuyo módulo, dirección y sentido es el de
la velocidad de cualquier punto del sistema. Cuando esta propiedad se cumple para un
solo instante, se dice que la traslación es instantánea.
Un ejemplo característico de este último caso es el mecanismo biela-manivela. En la
figura siguiente se observa que en la posición OBA el movimiento no es de traslación
ya que los vectores velocidad de los puntos A y B son distintos, mientras que en la
posición OB’A’ hay traslación instantánea porque los vectores velocidad de A’y B’ son
iguales.
Se dice que un cuerpo rígido se encuentra en este movimiento cuando dos de sus
puntos permanecen fijos durante el mismo respecto a un observador situado en 0,
perteneciente al marco de referencia. Sean esos los puntos 01 y 02 del gráfico:
63
Mecánica Racional
G G G
Vo1 Vo 2 0
64
Capítulo 2 – Cinemática de los Sistemas de Puntos
G G G G
V r1 Vo1 r1 0
G G G G
V r Vo r 0
G G G
de donde se deduce que V A r y r1 como se quería demostrar.
Esta propiedad unida a la condición geométrica de rigidez dice que en el movimiento
de rotación cada punto del sistema rígido realiza movimientos circulares con centro en
el eje de rotación y en planos normales al mismo. Es decir que si Z es la velocidad
angular de dicho movimiento, el módulo de la velocidad de un punto cualquiera es:
G
V Z r1
así
G G G
V Z r
G
V Z r senM Z r1
G
Obsérvese que mientras Z es un vector axial (puede desplazarse sobre su recta de
G
acción), V resulta ser aplicado (propio de cada punto material del sistema).
Ahora no sólo deberá estudiarse la variación de la velocidad de un punto en el tiempo
G G G
(a ) sino también la variación de Z en el tiempo (J ) . Supóngase para ello un vector
rotación en un eje definido por un versor Z :
65
Mecánica Racional
G G G
V Zr
G G G
G dV dZ G G dr G G G G
a r Z J r Z V
dt dt dt
G G G G G G
a J r Z Z r (2.8)
Se acostumbra denominar:
G G G
J r aceleración tangencial = at
G G G G
Z Z r aceleración normal = an
G G G
luego: a at a n
Ejemplo de aplicación:
El cono de la figura gira con Z = 3t + 2 [rad/seg] alrededor del eje k . Sabiendo que en
t = 0, era T = 0, determinar: a) velocidad y aceleración de P en t = 3 seg; b) ¿Cuántas
vueltas habrá girado hasta ese instante?
J 3
Solución: a) T Zot t 2 2t t 2
2 2
66
Capítulo 2 – Cinemática de los Sistemas de Puntos
G
r (t ) 3(cosT iˆ senT ˆj ) 12kˆ
G
V (t )
G
G G
>
Z r 3t 2kˆ 3 cosT iˆ senT ˆj 12kˆ @
V (t ) >
33t 2 cosT ˆj senT iˆ @
T (3) 19,5rad
G
V (3) > @
311 0,79 ˆj 0,60iˆ 19,9iˆ 26 ˆj>cm / s @
G G G G G G
a J r Z Z r
G
a (t )
G
> @ >
3kˆ 3 cosT i senT ˆj 12kˆ 3t 2 kˆ 33t 2 cosT ˆj senT iˆ @
a (t ) 9 cos T j sen T i 3 3t 2 cos T i sen T j
2
G
a (3)
7,11 ˆj 5,4iˆ 363 0,79iˆ 0,60 ˆj
>
292,17iˆ 210,69 ˆj cm / s 2 @
b) T ( 3) 19,5 rad 1117,26º 3,1 vueltas .
67
Mecánica Racional
G G
'r1 W 1 't
G G
'r2 W 2 't
#
G G
'rn W n 't
G
Nótese que W , que es la velocidad resultante, es otra traslación. Por lo tanto, el
movimiento resultante de varias traslaciones es una traslación. Tal vez sirva al lector
para la visualización del tema, el imaginar un punto sobre una hoja de papel que se
traslada sobre una mesa en traslación con respecto a una plataforma que a su vez se
traslada.
68
Capítulo 2 – Cinemática de los Sistemas de Puntos
G G G G
En tales condiciones, P tendrá las velocidades V1 y V2 originadas por Z1 y Z2 , según
las siguientes expresiones:
G G G
V1 Z1 r
G G G
V2 Z2 r
sumando miembro a miembro
G G G G G
V1 V2 Z1 Z 2 r
G G G
? V Zr (2.10)
Luego el movimiento resultante es una rotación pero con un solo punto fijo que es el
G
0. Por lo tanto, el vector rotación resultante Z no es fijo y por ello da lugar a la
denominada rotación instantánea. Al punto fijo O se lo llama polo y al movimiento se
lo denomina polar.
G G G
Generalizando la demostración con Z1 ,Z 2 , " ,Z n rotaciones concurrentes se deduce
que el cuerpo sometido a ellas posee un movimiento de rotación instantánea sobre un
eje que pasa por el polo, siendo en cada instante la rotación el vector resultante de las
rotaciones dadas.
n
G
Z ¦ZGi (2.11)
i 1
69
Mecánica Racional
aquí será: G G G
V1 Z1 r1
G G G
V2 Z2 r2
sumando miembro
G G aG miembro
G G G
V1 V2 Z1 r1 Z2 r2
G G G G G
V Z1 r1 Z2 r2
b.2.3.Par de Rotaciones:
G G
Si en el sistema rígido anteriormente dibujado las rotaciones Z1 y Z2 a que está
sometido el disco son además de paralelas, de la misma intensidad y de sentidos
opuestos, se tendrá:
70
Capítulo 2 – Cinemática de los Sistemas de Puntos
G G G
V1 Z1 r1
G G G
V2 Z2 r2
____________________________
G G G G G
V Z1 r1 Z 2 r2
G G G G
Z1 r1 Z1 r2
G G G G G
Z1 >r1 r2 @ Z1 r12 vectorcte
M1 Z1 't
M 2 Z1 't M1
G
Por lo tanto, V = vector constante
en cada posición. Es decir, todos los
puntos del disco tienen la misma
velocidad, puesto que si se hubiera
tomado otro punto para el análisis,
por ejemplo el P’, el resultado
hubiese sido el mismo. Luego, el
movimiento resultante es una
traslación curvilínea.
71
Mecánica Racional
G G G G
V W Z r (2.12)
72
Capítulo 2 – Cinemática de los Sistemas de Puntos
G
W S
G G G G
Vp W Z rp
G G G G
En este caso, es posible encontrar un punto C para el cual la suma W Z rc 0
En la siguiente figura se observa que cualquier punto ubicado sobre la recta oc posee
dos vectores velocidad colineales de sentidos opuestos: uno debido a la traslación y
G
otro a la rotación. Habrá entonces un punto C a una distancia rc desde O para el cual
ambos vectores sean iguales y de sentidos opuestos.
Todos los puntos de la normal al plano que pasa por C tendrán velocidad nula. Luego el
movimiento resultante será una rotación alrededor de ese eje de velocidades nulas, el
cual, por otra parte cambiará de posición con el tiempo debido a que es “arrastrado”
por la traslación.
Por lo tanto, se trata de una rotación instantánea. Al punto C se lo denomina centro de
rotación instantánea o polo de velocidades `
! $
!
!+!
&+
+
lo llama eje instantáneo de rotación.
G G G G G G
Z W Z Z rc 0
73
Mecánica Racional
G G G GG G G G G
Z W Z Zrc rc ZZ 0
G G G G
y tomando Z A rc Z rc 0
G G
y : Z Z Z2
G G
G Z W
resulta: rc (2.13)
Z2
Notar que la velocidad absoluta del punto del cuerpo que coincida con el centro
instantáneo de rotación es nula en dicho instante, pero su aceleración no lo es. Luego,
el centro instantáneo de rotación considerado como punto del cuerpo, no puede ser
tomado como centro instantáneo de aceleración nula.
G G
Volviendo ahora al caso general b.3 en que Z y W forman un ángulo cualquiera, el
mismo puede ser analizado por superposición de los dos casos anteriores:
74
Capítulo 2 – Cinemática de los Sistemas de Puntos
G G
Supóngase en el punto O una rotación Z y una traslación W que forman entre sí un
ángulo [.
G G
Componiendo Z con la proyección del vector traslación normal a ella W n se obtendrá
G
una rotación instantánea en el centro instantáneo C. Y siendo W P , la proyección del
G
vector traslación en la dirección de Z un vector libre, puede ser trasladado al punto C,
resultando un caso de composición de una rotación instantánea con una traslación
G G
paralela ( Z W ).
Luego, se desarrolla un movimiento helicoidal instantáneo, cuyo eje que pasa por C
cambia de posición con el tiempo. Este movimiento es el más general que puede tener
un sistema rígido.
La velocidad de un punto P cualquiera del cuerpo será:
G G G G
V W Z r
G G
G Z W
y la posición del centro instantáneo: rc
Z2
75
Mecánica Racional
G G G
V1 Z1 r1
G G G
V2 Z2 r2
#
G G G
Vn Zn rn
G nG n
G G
?VP ¦ i
V
i 1
¦ Z r
i 1
i i
G G
Las rotaciones Z1 y Z1 en 0 y R1 respectivamente, forman un par de rotaciones y por
lo tanto dan lugar a una traslación perpendicular al plano que ellas determinan,
ocurriendo lo mismo con el resto de los pares así formados. Siendo las traslaciones
vectores libres, se pueden llevar al punto 0, donde se tendrán las traslaciones
G G G
W 1 , W 2 ," ,W n que admitirán una resultante
n
G
W ¦WGi
i 1
76
Capítulo 2 – Cinemática de los Sistemas de Puntos
77
Mecánica Racional
Ejemplo demostrativo:
G
Sea un disco de radio r que gira con Z2 respecto del bastidor 0 0’ 0” el cual a su vez
G
gira con Z1 alrededor del eje 0 0”.
78
Capítulo 2 – Cinemática de los Sistemas de Puntos
G
La velocidad V1 del punto P1 del disco será:
G G G G G
V1 Z1 r1 Z 2 roc1
79
Mecánica Racional
Derivando con respecto al tiempo la (2.12) o la (2.14) puede ser encontrada la ley de
distribución de aceleraciones en el movimiento rototraslatorio.
Partiendo de (2.12):
G
G G dW G G G G
a V J r Z V (forma propia)
dt
tomando (2.14)
G
G dV1 G G G G G G
a2 J r2 r1 Z V2 V1
dt
o
G G G G G G G G G
a2 >
a1 J r2 r1 Z Z r2 r1 @ (2.15)
80
Capítulo 2 – Cinemática de los Sistemas de Puntos
Luego, se estudia el
movimiento del sistema
rígido analizando el de la
terna móvil solidariamente
unida al mismo.
Los problemas fundamentales que se plantean en el estudio del movimiento son los
siguientes:
1) Conocer la configuración (o posición) del sistema rígido en cada instante.
2) Conocer el estado de velocidad del sistema rígido, lo que implica poder calcular los
vectores velocidad de todos sus puntos.
3) Conocer el estado de aceleración.
81
Mecánica Racional
El proceso se inicia haciendo girar el sistema original (0, i, j , k) un ángulo \ sobre el
eje fijo k en sentido contrario al horario, obteniéndose el sistema (01 , iˆc, ˆj c, kˆ). Este
movimiento recibe el nombre de precesión y \ ángulo de precesión.
82
Capítulo 2 – Cinemática de los Sistemas de Puntos
Si se tienen estas funciones, se conoce para cada instante t la posición del cuerpo.
83
Mecánica Racional
y01 = l sen
z01 = h
2) Estado de velocidad
Se ha dicho que conocer el estado de velocidad de un cuerpo implica conocer los
vectores velocidad de todos sus puntos. Para ello, se tomará el origen de la terna móvil
como centro de reducción del movimiento, suponiendo conocida la velocidad de dicho
G G
punto V01 y el vector rotación Z en el mismo 01.
Sea la terna (01 , iˆ1 , ˆj1 , kˆ1 ) solidariamente unida al cuerpo. Para un punto genérico P del
cuerpo su vector posición con respecto a la terna en movimiento será:
84
Capítulo 2 – Cinemática de los Sistemas de Puntos
G
r x1 i1 y1 j1 z1 k1 (2.16)
con x1 , y1 , y z1 constantes
G G G
pero con respecto a la terna absoluta: r l l01
y derivando con respecto al tiempo:
G G G
dr dl dl01 G G
V V01
dt dt dt
Se analizará el significado de las derivadas de los versores con respecto al tiempo. Para
ello, comparando la expresión (2.17) con la forma impropia (2.14)
G G G G
V V01 Z r
pero es:
G G G
Z r Z x1 i1 y1 j1 z1 k1
G G G
x1 Z i1 y1 Z j1 z1 Z k1
di1 G
° Z i1
dt
°
° dj1 G
® Z j1 (2.18)
° dt
° dk1 G
° Z k1
¯ dt
85
Mecánica Racional
G G G G
V2 V1 Z r12 ½°
G G G G ¾ (2.19)
V3 V1 Z r13 °¿
Las ecuaciones (2.19) son vectoriales y de las mismas surgirán seis ecuaciones
G G G
escalares con doce parámetros: nueve componentes de las V1 , V2 y V3 y tres
G
componentes de Z .
Por lo tanto, doce parámetros menos seis ecuaciones de vínculo, resultan de nuevo
seis parámetros libres, por lo que el estado de velocidad de un cuerpo quedará
determinado cuando se conozcan seis parámetros de velocidad de tres de sus puntos
no alineados.
Las seis ecuaciones escalares obtenidas de las (2.19) serán:
86
Capítulo 2 – Cinemática de los Sistemas de Puntos
V2 x V1x Z y z2 z1 Z z y 2 y1
V2 y V1 y Z z x2 x1 Z x z2 z1
V2 z V1z Z x y 2 y1 Z y x 2 x1
V3 x V1x Z y z3 z1 Z z y 3 y1 (2.19’)
V3 y V1 y Z z x 3 x1 Z x z3 z1
V3z V1z Z x y 3 y1 Z y x 3 x1
V1x , V1 y , V1z
°
° V2 x , V2 y , V2 z
y los doce parámetros ®
° V3 x , V3 y , V3z
° Z x , Z y , Zz
¯
de los cuales seis deberán ser dados para poder determinar el estado de velocidad.
Un cuerpo rígido se mueve con respecto del marco de referencia absoluto. Los
datos geométricos y cinemáticos aportados están referenciados a la terna móvil fija al
cuerpo, y son:
\x = 2 rad/s
v1x = - 3 cm/s }
Resolución:
Haciendo uso de las ecuaciones
(2.19) y (2.19´), se obtiene:
v3x = v1x + 4 \y
v3y = v1y – 8 – 7 \z
v3z = v1z + 7 \y
v2x = v1x – 5 \z
v2y = v1y – 7 \z
87
Mecánica Racional
v2y = v1z + 10 + 7 \y
\y = 2 rad/s; \z = -1 rad/s; v1y = 6 cm/s; v2y = 13 cm/s; v1z = -22 cm/s y v3z = -8 cm/s
n
G
Z ¦ZGi (2.20)
i 1
Esto expresa que los vectores velocidad de un sistema material rígido proyectados en
un determinado instante sobre la dirección del vector rotación son una constante que
recibe el nombre de “invariante escalar P“.
88
Capítulo 2 – Cinemática de los Sistemas de Puntos
G
G G Z
P V01 Z V01 (2.21)
Z
G G
c) Si P r V01 debe ser Z V01 porque al proyectar, lo hace con su verdadero
valor Movimiento helicoidal permanente.
P z0 ½
d) Si G ¾ Movimiento helicoidal instantáneo
P z V01 ¿
En este último caso, tal como se vio en el apartado b.3), la velocidad o traslación
forma un ángulo G distinto de 0º ó 90º y teniendo en cuenta que la componente paralela
G
a Z de la V de cualquier punto es constante, al pasar de un punto al otro la
G
velocidad varía sólo por su componente perpendicular a Z . Existe entonces una recta
G
paralela a Z en cuyos puntos se anula la componente de la velocidad perpendicular a
dicha dirección y en tal caso, esa velocidad toma su valor mínimo. El lugar geométrico
de los puntos de velocidad mínima recibe el nombre de eje central del movimiento o
eje instantáneo del movimiento helicoidal. Dicha velocidad mínima representa el
vector traslación del movimiento helicoidal instantáneo:
G
W P Z (2.22)
G G G G G G G
Z V01 Z Z r 01 E
PZ Zˆ 0
y resolviendo:
G G G G G G G G G
Z V01 Z Z r0 E r0 E
Z Z 0
1
1
G G
tomando r01 E A Z Z2
GG
es Zr01 E 0
G G
G Z V01
finalmente; r01 E (2.23)
Z2
Una vez determinado E que es un punto del eje central, éste puede ser tomado
como centro de reducción para usar la forma propia de la ley de distribución de
velocidades:
G G G G
Vi W Z rEi
90
Capítulo 2 – Cinemática de los Sistemas de Puntos
Resolución:
91
Mecánica Racional
G G G
Vo Z2 rpo Z 2 i1 l j1 h k1
G
Vo Z2 l k1 h j1
Cualquier punto de la chapa que pasare
(hipotéticamente) por 0 tendrá siempre la
G
velocidad Vo calculada antes. Al colocar los
G G
vectores Z y Vo en 0 hemos reducido el
sistema a ese punto y a partir de él podemos
calcular la velocidad de cualquier punto de la
chapa.
G Z1 Z 2 l Z1 Z 2 l
W Z2 i1 Z1 k1
Z 2 i1 Z1 k1
Z12 Z 22
2
Z 22 Z12
4)
G G
G Z Vo Z 22 hkˆ1 Z 22 l ˆj1 Z1 Z 2 hiˆ1
r0 E
Z2 Z 22 Z12
92
Capítulo 2 – Cinemática de los Sistemas de Puntos
G
1º) Z Z2 i1 Z1 k1
G
2º) V p
G
Z1 r0 p
G
Z1 kˆ1 l ˆj1 hkˆ1 Z1 l iˆ1
G Z1 Z 2 l
3º) P V p Z
Z12 Z 22
G Z1 Z 2 l Z1 Z 2 l
W P Z Z2 i1 Z1 k1 2 Z 2 i1 Z1 k 1
Z1 Z 2
2 2
Z12 Z 22
4º)
G G
G Z Vp Z12 l ˆ
rPE j
Z2 Z12 Z22 1
G
rPE medido ahora desde P.
3) Estado de aceleración:
Tomando los mismos puntos del cuerpo que en el apartado anterior y trabajando
análogamente con la forma impropia de la ley de distribución de aceleraciones:
G G G G G G G
a2 a1 J r21 Z Z r12
®G G G G G G G (2.24)
¯ a3 a1 J r31 Z Z r13
G
Siendo: J J x i J y j J z k
93
Mecánica Racional
94
Capítulo 2 – Cinemática de los Sistemas de Puntos
G G G
dr dl dl01 dx1 dy1 dz1 di dj dk
i1 j1 k1 x1 1 y1 1 z1 1 (2.25’)
dt dt dt dt dt dt dt dt dt
G G
pero siendo l y l01 vectores de posición con respecto a la terna absoluta, sus
derivadas
G G temporales darán las velocidades de P y 01 con respecto al marco absoluto:
V y V01 .
Con respecto a los tres últimos sumandos del lado derecho de la igualdad, pueden
aplicarse las fórmulas de Poisson, obteniéndose:
G G G G G G
x1 : i1 y1 : j1 z1 : k1 : x1 i1 y1 j1 z1 k1 : r
Por lo tanto y teniendo en cuenta que los tres primeros sumandos representan la
velocidad de P como si la terna móvil estuviese quieta:
G G G G G
V Vrel V01 : r (2.26)
G
donde: V = velocidad absoluta de P
G G
dr
Vrel = velocidad relativa de P =
dt rel
donde el subíndice “rel” indica que se deriva el vector posición como si los ejes de
la terna de arrastre permanecieran inmóviles.
G G G
V01 : r = sería la velocidad de P como si el cuerpo estuviese solidariamente
unidoGa la terna móvil y fuese arrastrado por ella (velocidad de arrastre); así, rotaría
con : y 01 sería el centro de reducción del movimiento.
Luego:
G G G
V Vrel Varr (2.26’)
95
Mecánica Racional
G G § dx dy dz ·
= a rel : ¨ 1 i1 1 j1 1 k1 ¸
© dt dt dt ¹
G G G
= a rel : Vrel
G G G G G G G G G G G G G
= : r : Vrel : r : r : Vrel : : r
Reemplazando en (21)
G G G G G G G G
G
a
G
a 01 : r : : r a rel 2: Vrel (2.28)
G
donde: a aceleración absoluta de P
G
a rel aceleración relativa de P
G G G G G
G
a 01 : r : : r es la forma impropia de la ley de distribución de
aceleraciones en un sistema rígido (tal como si éste fuese arrastrado por la terna
móvil) y se denomina aceleración de arrastre.
G G
2: Vrel es la aceleración complementaria o de Coriolis. Aparece por la rotación
de los ejes de la terna móvil y representa la diferencia en aceleración de P como si
fuera medida a partir de unos ejes (0,i,j,k) no giratorios y de otros (01, i1, j1, k1)
giratorios, ambos con origen en 01. Se anula si no hay rotación o bien si no hay
movimiento relativo y también en los movimientos helicoidales permanentes donde
G G
: Vrel .
G G G G (2.28’)
Así resulta: a a rel a arr a com
96
Capítulo 2 – Cinemática de los Sistemas de Puntos
Resolución:
Se adoptará como marco de referencia
relativo el bastidor 00´P. En él, se
ubicará la terna móvil o de arrastre con
G G
origen en P, es decir, rotando con : Z
. Luego,
G G G
VrelR Z2 rPR
G G G G
VarrR VP Z1 rPR
G G G
donde V P Z1 rOP
G G G G
arelR Z2 Z2 rPR
G G G G G G G
aarrR aP J 1 rPR Z1 Z1 rPR
donde
G G G G
aP Z1 Z1 rOP
G
G d Z1 G
J1 0
dt
G G G
acomR 2Z1 VrelR
Las velocidades de todos los puntos que se encuentran sobre un radio, por ejemplo
P0, son proporcionales a sus distancias a 0, por lo que si se conoce una de esas
97
Mecánica Racional
G
velocidades, las demás surgen por proporcionalidad. La velocidad instantánea V y la
G
aceleración a de un punto cualquiera como el P del cuerpo vienen dadas por:
G G G
V Zr
G G G G G G
a Z r Z Z r
como en la rotación alrededor de un eje fijo pero con la única diferencia que si el
G G G
eje de rotación es fijo entonces J Z está dirigido según el eje fijo Z y representa la
G
variación del módulo de Z por unidad de tiempo; mientras que cuando el eje no está
G G
fijo en el cuerpo o en el espacio, el vector J ya no sólo reflejará la variación del Z
G G
sino también de la dirección de Z y no estará dirigido según el eje Z .
En general, en el caso de un cuerpo que gire alrededor de un punto fijo, el eje
instantáneo variará de posición tanto en el espacio como en el cuerpo.
Cuando el eje se mueve en el espacio genera un cono espacial o herpolodio y cuando
el eje se mueve respecto al cuerpo genera un cono relativo al cuerpo llamado cono del
cuerpo o polodio. Estos conos son tangentes a lo largo del eje instantáneo de rotación
G
Z
y el movimiento del cuerpo se puede describir como la rodadura del cono del cuerpo
sobre el espacial.
98
Capítulo 2 – Cinemática de los Sistemas de Puntos
En este caso, la poloide y la herpoloide son círculos y los conos del espacio y del
cuerpo son conos circulares rectos:
Caso 1: Caso 2:
Se dice que un sistema rígido está en movimiento plano cuando las velocidades
de cualquiera de sus puntos son paralelas a un plano fijo. Teniendo en cuenta esta
definición y la condición de rigidez, todos los puntos ubicados sobre una normal a
dicho plano tendrán igual velocidad, por lo que basta para estudiar el movimiento el
análisis de un único plano.
99
Mecánica Racional
Los sistemas rígidos que se mueven con estas características suelen ser
denominados “chapas rígidas” y en consecuencia se estudia el movimiento de una
chapa sobre el plano que la contiene.
En el movimiento plano se utilizarán dos ejes ortogonales como sistema de
referencia contenidos en el plano del movimiento y el tercero perpendicular a dicho
plano.
x x 2 y y 2 d 2
° 1 2 1 2 12
°° x x y y
2 2 2
d13
1 3 1 3
n=3 ® (2.29)
° x2 x3 y2 y3
2 2 2
d23
°
°¯ z1 z2 z3 0
x1 x2 2 y1 y2 2 2
d12 n 1
x1 y1
con 4 parámetros ® m=4
¯ x2 y 2
de donde: k = m - n = 3 gl
100
Capítulo 2 – Cinemática de los Sistemas de Puntos
1) Configuración:
Dado que en el plano el planteo se reduce notablemente, para hallar la
configuración es conveniente utilizar los cosenos directores de la terna de fija al
cuerpo. En efecto, los tres parámetros libres serían en este caso las 2 coordenadas del
origen de la terna móvil 01 , x 01 , y 01 y uno de los cosenos directores:
Aquí:
i1 c11 i c12 j
j1 c21 i c22 j
con
c11 cos D 11
c12 sen D 11
c21 sen D 11
c22 cos D 11
Con lo que se observa que con un único argumento se conocen los cuatro cosenos
directores.
2) Estado de velocidad:
Considerando los vectores velocidad de dos puntos cualquiera de la chapa:
G
°V1 V1x i1 V1 y j1
®G
°̄V2 V2 x i1 V2 y j1
Aplicando la condición
cinemática de rigidez entre ellos:
G G G G
G r1 r2 G r1 r2
V1 G G V2 G G (2.30)
r1 r2 r1 r2
101
Mecánica Racional
se obtiene una ecuación escalar con cuatro parámetros, de donde resultan tres
parámetros libres de velocidad. Dándole valores a éstos, se puede calcular el cuarto
mediante la ecuación (2.30).
G
Si el vector rotación Z no es conocido, se hace necesario determinarlo, aunque ya
se sabe que su dirección deberá estar sobre una perpendicular al plano del
movimiento.
Para lograrlo, se aplica la forma impropia de la ley de distribución de velocidades
entre los dos puntos tomando uno de ellos como centro de reducción:
G G G G G
V2 V1 Z r2 r1 (2.31)
V2 x V1x Z y 2 y1
(2.31’)
V2 y V1 y Z x 2 x1
Con las que se resuelve el problema del estado de velocidades en el plano. Por
supuesto que si se conociese el vector velocidad del origen 01, éste punto podría
tomarse como centro de reducción para determinar la velocidad de cualquier otro
punto Pi de la chapa mediante:
G G G G
Vi V01 Z ro1i
102
Capítulo 2 – Cinemática de los Sistemas de Puntos
3) Estado de Aceleración:
G
El planteo es totalmente análogo, pero debe tenerse en cuenta ahora el vector J .Si
G G G
se conociesen a 01 ,Z y J , la aceleración de cualquier punto quedaría expresada por:
G G G G G G G
a a 01 J r Z Z r
G
con: J J k1
G
Z Z k1
Cuando los vectores no son conocidos, hay que determinarlos a partir de ciertos
parámetros libres de la aceleración y cuyo número podría determinarse aplicando:
G G G G G G G
a2 a1 J r12 Z Z r12 (2.32)
°a
2x a1x J y 2 y1 Z 2 x 2 x1
®
a1 y J x 2 x1 Z 2 y 2 y1
°̄a 2 y (2.32’)
V Ax 1 cm/s
G
aA 3 i1 j1 cm/s2
a By 0
103
Mecánica Racional
G G
Una vez determinados los estados de velocidad y aceleración calcular V01 y a 01 y
hallar el centro instantáneo de rotación.
Solución:
G
a) Para determinar el estado de velocidad, se tendrá un parámetro por Z y cuatro
por las velocidades de A y B: cinco en total. Tomando A como centro de reducción se
obtienen dos ecuaciones escalares que las relacionan y por lo tanto tres parámetros
libres.
G G G G
VB V A Z rAB
V BX V AX Z y B y A
®
¯V BY V AY Z x B x A
reemplazando
4 = 1 Z (4 0)
3
Z ; V AY 8.
4
2 = VAy + Z (0 8)
G
VA i1 8 j1 (cm/s)
G
VB 4 i1 2 j1 (cm/s)
G 3
Z kˆ1 (rad/s)
4
°a
BX a AX J y B y A Z 2 x B x A
®
°̄a BY a AY J x B x A Z 2 y B y A
reemplazando
104
Capítulo 2 – Cinemática de los Sistemas de Puntos
§3·
2
°a BX 3 J 4 ¨ ¸ (8)
° ©4¹
® 2
° § 3· 5 (rad/s2)
°O 1 J ( 8) ¨ ¸ 4J
¯ ©4¹ 32
5 9 65
a BX 3 (cm / s 2 )
8 2 8
y por lo tanto
G
aA 3iˆ1 ˆj1 (cm / s 2 )
G 65 ˆ
aB i1 (cm / s 2 )
8
G 5ˆ
J k1 (rad / s 2 )
32
G G G
c) V01 V A Z rA01
G
Z rA 01
3
kˆ1 8iˆ1
4
6 ˆj1
G
V01 iˆ1 8 ˆj1 6 ˆj1 iˆ1 2 ˆj1 (cm/s)
G G G G G G G
a 01 a A J rA01 Z Z rA01
G G
J rA01
5 ˆ
32
k1 8iˆ1 5ˆ
4
j1
G G G 3 9ˆ
Z Z rA01 kˆ1 6 ˆj1 i1
4 2
G 13 ˆ 9 ˆ
a01 i1 j1 (cm/s2)
2 4
G G G G G G
Conociendo V01 , a 01 , Z , J pueden ser determinados los vectores Vi , ai para
G G
cualquier punto del plano, aplicando las leyes de distribución de V , a .
105
Mecánica Racional
G G G G
° Vi V01 Z ri
®G G G G G G G
°̄ai a 01 J ri Z Z ri
§ 3·
°°0 1 ¨ ¸ yc
©4¹ § 8 4·
® x x , yc ¨ , ¸ (cm)
°0 § 3· © 3 3¹
°¯ 2 ¨ ¸ x c
© 4¹
106
Capítulo 2 – Cinemática de los Sistemas de Puntos
G G
G Z V01
rc (2.33)
Z2
Ahora, reemplazando en esta última expresión los vectores referidos a las ternas
absoluta y móvil respectivamente, se obtienen las ecuaciones de la base y la ruleta.
a) Curva base:
Los vectores referidos a la terna absoluta son los siguientes:
G
rc xc x01 i yc y01 j
G
Z Z k
G dx 01 dy
V01 i 01 j
dt dt
107
Mecánica Racional
Z kˆ § dx01 ˆ dy 01 ˆ · 1 dx01 ˆ 1 dy 01 ˆ
xc x01 iˆ y c y 01 ˆj ¨ i j¸ j i
Z 2 © dt dt ¹ Z dt Z dt
1 dy 01
° xc x 01
Z dt
? ® (2.34)
° y 1 dx 01
y 01
¯ c Z dt
que son las coordenadas del polo con respecto a la terna absoluta o curva base.
Ahora teniendo en cuenta que:
dM
Z
dt
y reemplazando en la (2.34):
dy01
°° xc x01
dM
® (2.35)
° dx
°¯ yc y 01 01
dM
Estas últimas también son las ecuaciones paramétricas de la curva base, y permiten
observar que dicha curva es independiente del tiempo. En otras palabras, no importa
cuál sea la velocidad de la chapa al describir su movimiento, la trayectoria de C será la
misma.
b) Curva ruleta:
Referidos a la terna móvil, los vectores de la (2.33) son:
G
rc x c1 i1 y c1 j1
G
Z Z k 1
G
para obtener V01 expresado en el plano móvil, se proyectan sus componentes:
G ª§ dx 01 dy 01 · º ª§ dx dy · º
V01 «¨© dt i dt j ¸ i1 » i1 «¨ 01 i 01
¹ ¼ ¬© dt
j ¸ j1 » j1
¹ ¼
¬ dt
y siendo:
108
Capítulo 2 – Cinemática de los Sistemas de Puntos
iˆ iˆ1 cos M
ˆj iˆ1 cos M 270D sen M
iˆ ˆj1 cos M 90D sen M
ˆj ˆj1 cos M
resulta:
G § dx dy · § dx dy ·
V01 ¨ 01 cos M 01 sen M ¸ i1 ¨ 01 sen M 01 cos M ¸ j1
© dt dt ¹ © dt dt ¹
O, vectorialmente
G G ª cos M senM º
V c RV01 con R « senM
¬ cos M »¼
Reemplazando en (2.33):
ª§ dx dy ·
Z k1 «¨ 01 cos M 01 sen M ¸ i1
¬© dt dt ¹
x c1 i1 y c1 j1
Z2
§ dx 01 dy 01 · º
¨ sen M cos M ¸ j1 »
© dt dt ¹ ¼
Z2
resultando:
1 § dx01 dy01 ·
° x c1 ¨ sen M cos M ¸
° Z © dt dt ¹
®
° 1 § dx01 dy ·
° yc1 ¨ cos M 01 sen M ¸
¯ Z © dt dt ¹
Estas ecuaciones dan las coordenadas del polo en los ejes móviles en forma
paramétrica (curva ruleta). Como antes, puede eliminarse el tiempo obteniéndose:
dx01 dy
°° x c1 sen M 01 cos M
dM dM
® (2.36)
° y dx01 dy
°¯ c1 cos M 01 sen M
dM dM
x 01 x 01 (M ) , y 01 y 01 (M )
109
Mecánica Racional
Solución:
Para determinar las coordenadas del polo C en el plano absoluto se aplican las
ecuaciones (2.35);
° x01 A cos 180q 90q M A cos 90q M A sen M
®
°̄ y01 0
reemplazando:
xc A sen M
®
¯ yc A cos M
° xc1 A cos M sen M
®
°̄ yc1 A cos 2 M
eliminando el parámetro M:
yc1
xc1 A cos M 1 cos 2 M ; cos 2 M
A
110
Capítulo 2 – Cinemática de los Sistemas de Puntos
yc1 y
? xc1 A 1 c1 Ayc1 yc21
A A
x a 2 y b
2
r2
En este caso es a = 0
x 2 y 2 2 y b b2 r2
y tomando b = 1/2 = r
x2 y 2 2 b y A y
Por lo tanto:
BASE : circunferencia de centro 0 y radio l
RULETA : circunferencia desplazada del origen 01 sobre el eje j1 de radio l /2.
111
Mecánica Racional
Propuesta para el alumno: Dibujar los ejes de referencia absolutos en una hoja papel.
Dibujar los ejes de referencia móviles y la barra sobre otra hoja, en lo posible
transparente. Superponerlos partiendo de 0 coincidente con 01 y comenzar el
deslizamiento del móvil sobre el fijo manteniendo los extremos de la barra en contacto
con los ejes fijos. Marcar con la punta de un lápiz en varias posiciones la posición de C
en los dos sistemas simultáneamente. Separar y observar las trayectorias del polo.
112
Capítulo 3- Cinética del Punto Material
Las ideas actuales acerca del movimiento de los cuerpos se remontan a Galileo y
Newton. Antes de ellos, se creía en las ideas de Aristóteles, quien sostenía que el
estado natural de un cuerpo era estar en reposo y que éste sólo se movía si era
empujado por una fuerza o un impulso. De ello se deducía que un cuerpo pesado debía
caer más rápido que uno liviano, porque sufría una atracción mayor hacia la Tierra.
Galileo demostró que las anteriores ideas estaban equivocadas dejando caer esferas
de distintos pesos a lo largo de un plano inclinado. Demostró que cada cuerpo
aumentaba su velocidad al mismo ritmo, independientemente de su peso. Por
supuesto que una esfera de plomo caerá más rápida que una pluma, pero ello se debe
únicamente a que la pluma es frenada por la resistencia de aire. En los experimentos
de Galileo, cuando un cuerpo caía rodando, siempre actuaba sobre él la misma fuerza
(su peso) y el efecto que se producía consistía en acelerarlo en forma constante. Esto
demostraba que el efecto real de una fuerza era el de cambiar la velocidad del cuerpo,
en vez de simplemente ponerlo en movimiento, como se pensaba antes. Ello también
significaba que siempre que sobre un cuerpo no actuara ninguna fuerza, éste se
mantendría moviéndose en una línea recta con la misma velocidad.
Los experimentos de Galileo sirvieron de base a Newton para la obtención de sus leyes
del movimiento.
En el presente Capítulo se estudiarán las relaciones entre el movimiento de los puntos
materiales y las causas que los provocan. Recuérdese que para la Mecánica las
partículas o puntos materiales, son entes desprovistos de dimensiones (punto
geométrico en el sentido del tamaño), aunque se admite la abstracción de que pueden
poseer masa.
113
Mecánica Racional
F1 F2 Fn
! m
a1 a2 an
G G
en todos los casos el cociente será el mismo y como F y a son vectores, ambos
tendrán igual dirección y sentido. La magnitud modular m recibe el nombre de masa
de la partícula y proporciona una medida cuantitativa de la resistencia que opone la
misma a que su vector velocidad sea cambiada, propiedad que se denomina inercia.
Si la masa de la partícula es m, esta primera ley queda expresada matemáticamente
como: G G
F ma (3.1)
n G
G
¦ Fi m a (3.1’)
i 1
El primer término, comprende todas las fuerzas actuantes sobre la partícula, es decir la
suma de las fuerzas activas y las reactivas.
114
Capítulo 3- Cinética del Punto Material
Considérese
G aGmodo de ejemplo una partícula de masa m sometida a la acción de dos
fuerzas F1 y F2 :
El diagrama de cuerpo libre representa gráficamente a todas las fuerzas que actúan
n G G G
sobre la partícula ¦
Fi F1 F2 , mientras que el diagrama cinético toma en cuenta
i 1
G
gráficamente el vector m a .
Como se observa en los diagramas, la partícula se acelera en la dirección de
G G G
¦ Fi F1 F2
G
de manera que el módulo de a resulta:
G
a
¦Fi
G m G
Se comprueba que si en (3.1’) es ¦ Fi 0, la aceleración también es nula, en cuyo
caso la partícula permanecerá en reposo o en movimiento rectilíneo uniforme.
G
V constante vectorial (MRU)
G
a 0
G G
V 0 (Reposo)
m1 m2
F G (3.2)
r2
115
Mecánica Racional
La formulación de Mach:
Debido a las críticas que sufrieron los postulados de Newton, Mach reformuló los
fundamentos de la Mecánica estudiando las interacciones entre dos y luego tres o más
puntos aislados, arribando a las mismas conclusiones presentadas por Newton en su
teoría, aunque con mayor rigurosidad físico-matemática.
116
Capítulo 3- Cinética del Punto Material
x = X – Vo´ t
y=Y (3.3´)
t = t´
Derivando a (3.3) dos veces con respecto al tiempo (y teniendo en cuenta que Vo’ y D
son constantes), resulta:
G G
d 2r d 2 R
dt 2 dt 2
o, escalarmente
d 2x d 2X d 2 y d 2Y
;
dt 2 dt 2 dt 2 dt 2
Terna en Rotación:
Se tomará para el análisis el caso común que se presenta con el uso de la Tierra como
marco de referencia para el estudio de problemas de ingeniería. Los sistemas de
coordenadas solidarios a la Tierra se denominan “ternas del lugar”.
El eje k c está dirigido desde el centro de la tierra hacia afuera, el eje j c es tangente a
un paralelo y el eje i c tangente a un meridiano. Se considera inercial a un sistema
solidario con el Sol, suponiendo que el centro de éste no posee aceleración.
G
Posicionando a P con el vector r c respecto de 0’, se estudiará el movimiento relativo
de P respecto del sistema móvil con la Tierra (0c, i c, j c, k c) . Pero la aceleración
G G G G
absoluta de P respecto del Sol estará dada por a a rel a arr a comp
117
Mecánica Racional
O2 ª d 2 1 e 1 º
Según Binet: ae 2« » (3.4)
e ¬ dT 2 e¼
O 2 s
H = excentricidad de la órbita.
p = parámetro de la cónica =
b 2 a.
y como es:
1 1 H cosT (3.5)
e p
entonces:
d21 e H
cosT (3.6)
dT 2 p
O21
Reemplazando (3.5) y (3.6) en (3.4) resulta a e
p e2
4S 2 a 3 cm
resulta ae 2 #1 # 0.001g
e T
4
seg 2
118
Capítulo 3- Cinética del Punto Material
a R Z 2 r cos M
M latitud
2S
Z
86164 seg (tiempo de un dia sideral )
RT r 6378 Km
cm
a R 2,5 0.0025 g
seg 2
Por lo que la aceleración de arrastre oscila en el orden de los 3,5 cm/s2 (0,0035 g).
G
Dada la pequeñez de Z y tomando Vrel = 60 m/seg, por ejemplo, resulta:
G G
Por lo tanto, a a rel sólo alcanza milésimos del valor de la aceleración gravitatoria
G
terrestre g 0,005 g
Este hecho permite concluir que las aceleraciones producidas por el movimiento
terrestre pueden despreciarse en la mayoría de los cálculos de ingeniería, por lo que
los problemas dinámicos relacionados con movimientos sobre o cerca de la superficie
terrestre pueden resolverse usando como marco de referencia a la Tierra y un sistema
de coordenadas fijo a él, que se supone inercial.
Cuando se estudian los movimientos de cohetes y satélites r !! RT es necesario
considerar el sistema inercial de referencia como fijo a las estrellas para obtener
mayor precisión en los cálculos.
¦ Fx m ax ¦ Fy m ay ¦ Fz m az (3.7)
119
Mecánica Racional
G
¦F ¦ F ¦ F ¦ F 2 2 2
Con: x y z
G
y: a a a x2 a y2 a z2
¦ Fe m ae , ¦ FT m aT , ¦ Fz m az (3.7’)
¦F n m an , ¦Ft m at
V2 l
donde: an V T U T 2 y at V
U
G
Al término m a suele denominárselo reacción de inercia (igual y de sentido puesto
a la fuerza directriz).
Esta nueva forma (que no debe interpretarse como un simple pasaje de términos)
establece un equilibrio dinámico entre fuerzas actuantes y la resistencia al cambio de
velocidad de la partícula. Gráficamente:
120
Capítulo 3- Cinética del Punto Material
G G
y llamando reacción de incercia o fuerza ficticia a: ma FI , queda:
G G G G
¦F ¦F A R FI 0
121
Mecánica Racional
G G G G
¦ r
F 0 ? Fp 0
en un punto material libre no hay fuerza perdida, toda la fuerza activa produce
movimiento:
G G G G
Fp ¦F A FI 0
?
G G G
¦ FA FI ma
G
es decir que las FA por sí solas constituyen la fuerza directriz del movimiento.
c) Como conclusión de los casos límites, puede observarse que en todos los casos
intermedios parte de las fuerzas activas dan lugar al movimiento mientras que el resto
se pierde en los vínculos.
3. 4. Fuerzas naturales
La Mecánica establece como hipótesis que las fuerzas actuantes sobre un punto
material dependen en el caso más general de la posición de la partícula, su velocidad y
el tiempo. Así:
G G G G G§ dx dy dz ·
F F r ,V , t F ¨ x , y , z , , , , t ¸
© dt dt dt ¹
122
Capítulo 3- Cinética del Punto Material
Estas fuerzas, por su naturaleza,G dan lugar a los campos de fuerza (que son el espacio
donde está definida la función F ) donde para cada tres valores x, y, z corresponde un
único vector fuerza.
En estos campos vectoriales, más que la fuerza definida en cada punto del campo,
interesa conocer la intensidad del mismo en cada punto, que sería la fuerza que
corresponde a un punto material de masa unitaria colocado en el campo, es decir, la
fuerza por unidad de masa. G
En los campos de fuerza que se presentan en la naturaleza se verifica que si H es la
intensidad y m la masa, la fuerzaG que actúa
G en cada punto es
F m H
Un campo de fuerzas, entonces, queda definido por su intensidad, dado que la fuerza
depende de la masa que se coloque en el campo. Un ejemplo de tales campos es el
G
gravitatorio, en el cual la intensidad es el vector aceleración de la gravedad g , ya que
el peso depende de la masa que se coloque en cada punto:
G G
P m g
G
Un campo se dice que es uniforme cuando su vector intensidad H es constante para
toda la región del espacio en la que está definido. Cuando se trabaja en una porción
reducida del espacio, el campo gravitacional se toma como uniforme.
123
Mecánica Racional
G G G
H dr 0 (3.10)
ó
iˆ ˆj kˆ
G
Hx Hy Hz 0
dx dy dz
Por lo tanto
Hy dz Hz dy iˆ
Hz dx Hx dz ˆj
G
Hx dy Hy dx kˆ 0
dx dy dz (3.10’)
Hx Hy Hz
n
k P o P1 v + P 2 v 2 " P n v n ¦P v
i 0
i
i
y donde:
124
Capítulo 3- Cinética del Punto Material
G rG2 G G G G
W1 2 F ³G dr F r2 r1 Fd12 cos D
r1
G G G
siendo D el ángulo que subtienden F y r2 r1 ; en el caso que sea nulo, la expresión
del trabajo realizado por la fuerza constante colineal con el desplazamiento toma la
conocida expresión
W1-2 = F d1-2
125
Mecánica Racional
G G
b) Fuerza dependiente del tiempo: F F (t )
G G
dW F (t ) dr
G
para integrar, se debe conocer dr en función del tiempo:
G G
dr V (t ) dt
luego:
G G
dW F (t ) V (t ) dt
dW Fx Vx Fy V y Fz Vz dt f (t ) dt
³
t2
? W1 2 f (t ) dt
t1
G G G
c) Fuerza dependiente de la posición: F F (r )
dW Fx dx Fy dy Fz dz
³ ³ ³
x2 y2 z2
W1 2 Fx dx Fy dy Fz dz
x1 y1 z1
G wu wu wu
F u ( x , y , z ) i j k
wx wy wz
126
Capítulo 3- Cinética del Punto Material
iˆ ˆj kˆ
G w w w G
rot F 0
wx wy wz
Fx Fy Fz
w Fx w Fy w Fz w Fx w Fz w Fy (3.12)
; ;
wy wx wx wz wy wz
Resulta importante analizar si para todos los campos irrotacionales las fuerzas derivan
de una función escalar uniforme (simplemente valuada).
El trabajo aquí será:
§w u w u w u ·
dW ¨ i j k ¸ dx i dy j dz k
©w x wy wz ¹
wu wu wu
dx dy dz du
wx dy wz
u2
W1 2 ³
u1
du u 2 u1 (3.13)
G G G
³
l
F dr ³³s
ˆ
n rot F ds
G
en el cual la integral de línea representa el trabajo de circulación de la fuerza F y la
integral de superficie es el flujo del rotor del campo a través de la superficie limitada
por l.
G
Si se tiene un campo irrotacional, es decir rot F = 0, resultará:
127
Mecánica Racional
G G
³ dr
F
l
0 ó W12 u 2 u1 0 si P1 { P2
El trabajo a lo largo de una línea cerrada es cero, lo que será válido si se cumplen las
condiciones del Teorema de Stockes, es decir, que el campo sea simplemente conexo.
G
El trabajo de las fuerzas F del campo a
lo largo de un camino que vaya de P1 a
P2 por la línea 11 será W12 = u2 - u1 ya
que esta región es de simple conexión.
Sin embargo, si se va a lo largo de la
línea 12 (espacio de múltiple conexión) el
trabajo W sería diferente, ya que se lo
podría considerar como la suma de los
siguientes trabajos:
W1cc 2cc * módulo del campo (trabajo de circulación en una línea cerrada que
limita la zona no definida)
128
Capítulo 3- Cinética del Punto Material
W12 u 2 u1 *
Esta última expresión establece que el trabajo entre dos puntos para un camino
múltiplemente conexo es igual al trabajo para un camino simplemente conexo más
tantas veces el módulo del campo como vueltas desarrolla el camino alrededor de la
zona no definida.
Como conclusión de este análisis, se observa que las ecuaciones (3.12) resultan ser una
condición necesaria pero no suficiente para que un campo sea conservativo. Estas
aseguran que el campo es irrotacional, pero será conservativo si además la función
potencial u = u(x,y,z) es uniforme, es decir, está definida en toda la región. Esto
garantizará que el trabajo W a lo largo de una línea cerrada sea nulo.
El peso de una partícula y la fuerza de un resorte elástico son dos ejemplos de fuerzas
conservativas que se encuentran a menudo en Mecánica. En el primer caso, será:
G G
F mg j ; dr dx i dy j dz k
G
r2 G G y2
W1 2 ³
G
r1
mg dr m g ³ dy
y1
Para el resorte, tomando x=0 en la posición del mismo sin deformar (x1 y x2 coinciden
con su elongación), es:
129
Mecánica Racional
W1 2
1 2 K x 22 x12
El trabajo sólo depende de las longitudes inicial x1 y final x2 del resorte.
d) Fuerzas de fricción:
El trabajo de una fuerza de fricción ejercida sobre un objeto móvil por un medio fijo
depende de la trayectoria (cuanto más larga sea ésta, mayor será el trabajo). Por
consiguiente, las fuerzas de fricción no son conservativas. Generalmente el trabajo se
disipa en forma de calor.
3.5.2. Potencia:
Este concepto asocia el trabajo W al tiempo empleado en producirlo.
Potm = W / 't
G
La potencia media en 't = t1 –t0 expresa el valor medio del trabajo efectuado por F en
un intervalo de tiempo.
Si 't o 0,
G G
Pot = dW/dt = F v = Fx vx + Fy vy + Fz vz Potencia instantánea
130
Capítulo 3- Cinética del Punto Material
G
G G dV
¦F ma m
dt
G
y multiplicando miembro a miembro por dr :
G G
G G dr G G
¦ F drG
dV G
m dr m dV m V dV
dt dt
G
Introduciendo m V dentro del diferencial:
§ V2·
dW d ¨m ¸
© 2 ¹
m V 2
El término del miembro derecho recibe el nombre de energía cinética (e).
2
Luego:
dW de (3.14)
Esta ecuación describe el principio del trabajo y la energía (o Teorema de las fuerzas
vivas). Integrando entre dos instantes t1 y t2:
³ ³
P2 V2
dW de
P1 V1
ó
1 1
W1 2 mV22 mV12 (3.14’)
2 2
Aquí W1-2 representa el trabajo total realizado por todas las fuerzas que actúan sobre
la partícula cuando la misma se mueve desde el punto P1 al P2. Los términos de la
derecha son cantidades positivas por cuanto no dependen de la dirección de la
velocidad. En otras palabras: “la variación de la energía cinética de un punto material
en un intervalo de tiempo t2 - t1, iguala al trabajo realizado por las fuerzas en ese
mismo intervalo”.
G
En el caso particularmente importante en el cual la fuerza F es posicional y además
es conservativa (es gradiente de una función escalar u), utilizando coordenadas
cartesianas se tendrá:
donde u(x,y,z) es una función escalar cualquiera que se denomina función potencial.
Ahora bien, observando la expresión anterior es lógico inferir que si las fuerzas del
campo realizan un trabajo positivo generando una energía cinética, lo tendrán que
hacer a expensas de un consumo de trabajo que las mismas están capacitadas para
realizar.
131
Mecánica Racional
Así, puede definirse en el campo conservativo una cierta función energía potencial (p)
como energía capaz de realizar trabajo:
dW = - dp (3.15)
ya que a un trabajo positivo de las fuerzas del campo le corresponderá una cantidad
igual pero de signo contrario de p.
La función energía potencial pG será igual y de signo contrario a la función escalar
u(x,y,z) generadora del campo F .
La relación de p con la energía cinética es
de = - dp (3.16)
Si se pasa de un punto P1 con la velocidad V1, a un punto P2 con la velocidad V2, siendo
p1 y p2 los respectivos valores de la función energía potencial, de la integración de
(3.16) se tiene:
³ ³
V2 p2
de dp
V1 p1
e2 e1 p2 p1
? e2 p2 e1 p1 cte. E (3.17)
132
Capítulo 3- Cinética del Punto Material
G G G G
Luego: F dt mdV
d mV dQ
G G G
A F dt dI se lo denomina impulso elemental y a d (mV ) cantidad de movimiento
G
elemental dQ .
G G
Así: dI dQ (3.18)
Si ocurre que no hay fuerzas que actúen sobre el punto P, o bien la resultante de ellas
es nula, de (3.18) se tiene:
G G G
dQ 0Q vector constante
G
lo que expresa que en este caso se conserva el vector cantidad de movimiento Q .
G N G
Para un conjunto de N puntos materiales, es: K( o) ¦ rGi mi Vi
1
133
Mecánica Racional
dV
luego: mg PV m
dt
m dV
dt dV
m g PV P
g V
m
§ P · m
Integrando por sustitución: ¨ u g V ¸, dV du
© m ¹ P
1 § P ·
t ln ¨ g V ¸ C1
P /m © m ¹
x xo 0
recurriendo a las condiciones iniciales: t 0®
¯V Vo 0
m m
0 ln g C1 C1 ln g
P P
m § P · m
? t ln ¨ g V ¸ ln g
P © m ¹ P
§ P ·
m g ¨ t
V 1 e m ¸¸ (3.20)
P ¨© ¹
dx
V
dt
? dx
mg
P 1 e P / m t
dt ;
integrando
134
Capítulo 3- Cinética del Punto Material
mg m g e P / m t
x t C2
P P P /m
pero en t = 0 es x = xo = 0, de donde:
m2
C2 g y:
P2
m m2 g P / m t m2 g
x gt e 2 (3.21)
P P 2
P
de manera que para conocer el tiempo de caída tc basta con reemplazar x por H en la
(3.21).
mg
Vf en x o f (3.23)
P
135
Mecánica Racional
K = constante de proporcionalidad =
GM
G
r x i y j z k
r2 x2 y2 z2
G
r
r
r
G
La expresión de F en coordenadas cartesianas, será:
G ª º
mK« x y z
F 2 i j k»
x y 2 z 2 «¬ x 2 y 2 z 2 x2 y2 z2 2 2
x y z 2 »
¼
es decir:
G
F
mK
3 xi yj zk
x2 y2 z2 2
donde es:
mK x
Fx
x 2
y2 z2
3
2
mK y
Fy
x 2
y2 z2
3
2
mKz
Fz 3
x2 y2 z2 2
de donde puede verificarse que tal campo es irrotacional, puesto que resulta de (3.12):
w Fx w Fy w Fy w Fz w Fx w Fz
; ;
wy wx wz wy wz wx
136
Capítulo 3- Cinética del Punto Material
mK mK
?dp dr p C
r2 r
y tomando p o 0 si r o f , resulta: C = 0
mK mK G
luego: p y será : u F u
r r
dx dy
Kx kKy
x 2
y z2 2
3
2
x 2
y2 z2 3
2
dx dy
?
x y
Integrando:
ln x ln y ln C ln x ln C y
? y = C1 x
haz de rectas por el origen contenidas en el plano coordenado (x,y). De igual forma
pueden conocerse las líneas de fuerza en los otros planos coordenados.
r2 x2 y2 z2 m2 K 2 / p2
Para cada valor p0 de p existe una superficie equipotencial. Estas son una familia de
esferas concéntricas en 0 y por lo tanto son perpendiculares a las líneas de fuerza.
137
Mecánica Racional
e1 p1 e2 p2
mV12 m K mV22 m K
de donde:
2 r1 2 r2
y por lo tanto:
V12 K K
V2 1,41
2 r1 r2
En movimiento, se plantea:
d2y
¦F m
dt 2
138
Capítulo 3- Cinética del Punto Material
d2y
? Ky m
dt 2
d2y K
y 0 (3.24)
dt 2 m
Que es la ecuación diferencial del movimiento oscilatorio armónico, cuya solución es:
Ecuación característica:
K K
O2 0 O1, 2 ri
m m
K K
? y (t ) C1 cos t C 2 sen t
m m
y (t ) Y sen Zn t M
y (t ) Y cosZ n t T
y(t 0) y o
° (3.25)
® dy
°̄ dt (t 0) Vo
K K K K
y (t ) C1 sen t C2 cos t
m m m m
K m
? y (0) C 2 Vo C 2 Vo
m K
K m K
y (t ) y o cos t Vo sen t (3.26)
m K m
139
Mecánica Racional
yo Zn
tg M
Vo
Luego: y (t ) Y sen M cos Znt Y cosM sen Znt Y sen Znt M
y (t ) Y sen Zn t M (3.28)
La expresión (3.28) permite asegurar que un ciclo del movimiento tiene lugar cuando
ZnT varía 2 S rad. Es decir que el período T, verifica:
2S
Z n T 2S T
Zn
Como la frecuencia es fn = 1/T, resulta que la frecuencia natural del sistema viene dada
por:
1
fn K/m (ciclos/seg o Hertz)
2S
En conclusión: el movimiento del sistema es armónico (la relación entre y y t está dada
por funciones trigonométricas), con frecuencia natural:
Zn
fn
2S
140
Capítulo 3- Cinética del Punto Material
El análisis de sistemas con amortiguamientos resulta muy complicado; sin embargo,
existen “modelos” de amortiguamiento ideal que se adaptan bien a ciertos casos
particulares. Uno de ellos es el que permite el tratamiento matemático más simple y se
basa en la hipótesis de que la fuerza del amortiguador es proporcional a la velocidad.
Tal como se vio en 3.4 b), a esta fuerza generada por el amortiguador se la suele llamar
dy
de resistencia viscosa F = - c. , donde c recibe el nombre de coeficiente de
dt
amortiguamiento viscoso.
d2y dy
m 2 c Ky
dt dt
d2y d y
? m 2
c K y 0 (3.29)
dt dt
d 2 y c dy K
y 0
dt 2 m dt m
c K
O2 O 0 donde:
m m
141
Mecánica Racional
§ c 2 ·
¨ r §¨ c ·¸ 4 K ¸
¨ m ©m¹ m¸
O1, 2 © ¹
2
c § c ·2 K
O1,2 r ¨ ¸ (3.30)
2m © 2m ¹ m
2
§ c · K
a) Si ¨ ¸ ! lo que ocurre es que c tiene valor elevado y hay un gran
© 2m ¹ m
amortiguamiento. En este caso O1 y O2 resultan ser reales y por lo tanto la solución
general de (3.29) viene dada por:
c § c ·2 K
Llamando: D y E = ¨ ¸
2m © 2m ¹ m
O1,2 D r E
y (t ) C1 e (D E )t C 2 e (D E )t (3.32)
2 K
como es D ! 0 y E ! 0 y E D2 D 2 ? E D
m
resulta D E ! 0 y D +E >0
Luego, los exponentes de las funciones exponenciales son negativos, es decir que la
solución y(t) viene dada por la suma de dos exponenciales decrecientes. Así, en este
caso, el sistema no oscila. Al crecer t la masa tiende a la posición de equilibrio estático
sin oscilar. Este movimiento se llama sobreamortiguado, y sus gráficas para distintas
condiciones iniciales son:
142
Capítulo 3- Cinética del Punto Material
§ c ·2 K
b) Si es: ¨ ¸ , la solución general de (3.29) resulta:
© 2m ¹ m
c
C1 C 2 t e
t
2m
y (t ) (3.33)
2
§c · K
En este caso, c se designa con cc; despejando de ¨ c ¸ , resulta:
© 2m ¹ m
cc 2 mk 2mZ n
C1 C sen M
con: ®
¯C2 C cos M
? y (t ) C e Dt sen Et M (3.34)
143
Mecánica Racional
2
K § c ·
La frecuencia circular natural reducida E ¨ ¸ , es mayor cuanto más
m © 2m ¹
pequeño es c, y por lo tanto las oscilaciones son más rápidas.
Donde:
T = 2S
E es ahora el pseudoperíodo.
K
Nótese que cuando c o 0 E o y resulta el movimiento de vibraciones libres
m
sin amortiguamiento.
Efectuando el cociente:
144
Capítulo 3- Cinética del Punto Material
c
t c
y1 e 2m 1 T
c
e 2m
y2 (t T )
e 2m 1
luego:
y1 c
ln T
y2 2m
donde:
2m y1 c 1 y
c ln ó ln 1
T y2 cc T Zn y2
dy
y (o ) yo ; (o) Vo (3.36)
dt
145
Mecánica Racional
a) F(y,t) = Fo y. t
b) F(y,t) = Fo y2. sen Zt (3.37)
2
c) F(y,t) = Fo y
En el caso a) se tiene una ecuación diferencial con coeficientes variables, en b) la
ecuación es no lineal con coeficientes variables y en c) es no lineal.
Si la ecuación diferencial no es lineal, la solución del problema se complica y aún para
los casos más simples de la función F(y) es necesario recurrir a métodos aproximados.
Se analizará el caso en que la fuerza externa F aplicada al sistema es una función
armónica de t:
F(t) = Fo cos Zf t (3.38)
d2y
Sea: m K y Fo cosZ f t o
dt 2
d2y K Fo
2
y cosZ f t (3.39)
dt m m
yp A cos Z f t B sen Z f t
siendo yp A Z f sen Z f t B Z f cosZ f t
yp A Z 2f cosZ f t B Z 2f sen Z f t
Reemplazando en (3.39):
A Z 2f cosZ f t B Z 2f sen Z f t
k
m
>
A cos Z f t B sen Z f t @ Fo
m
cos Z f t
Igualando coeficientes:
ª Kº ª K º Fo
cosZ f t « AZ 2f A » senZ f t « B BZ 2f » cosZ f t
¬ m¼ ¬ m ¼ m
°A K Fo Fo ª 1
º
A Z 2f A « »m
° m m m «¬ K m Z 2f »¼
? ®
° K
°B B Z 2f 0B 0
¯ m
146
Capítulo 3- Cinética del Punto Material
§ 1 ·
¨ ¸
yp Fo ¨ 2 ¸ cos Z f t
© K m Z f ¹
ª 1 º
?y (t ) yc y p C1 cosZ n t C 2 senZ n t Fo « 2 »
cosZ f t
¬« K mZ f ¼»
Por lo tanto, la solución general del problema de condiciones iniciales (3.39 – 3.36)
d2y K Fo y (o) yo
y cos Z f t ®
dt 2 m m ¯ y(o) Vo
resulta ser:
§ Fo · § F · V
y (t ) ¨¨ yo ¸ cos Zn t ¨ o ¸ cos Z f t o sen Zn t
K Z 2f m ¸¹ ¨ K Z 2 m¸ Zn
© © f ¹
Vo § F · § F ·
yo cosZnt ¨
sen Znt ¨ o ¸ ¨ o ¸
2 ¸ cos Znt ¨ 2 ¸ cos Z f t
y (t ) (3.40)
Zn © K Z m
f ¹ © K Z m
f ¹
Debe notarse que los tres primeros sumandos tienen la frecuencia natural Zn del
sistema, mientras que el último la frecuencia de la fuerza externa Zf.
§ Fo · V Fo
y (t ) ¨ yo ¸ cos Zn t o sen Zn t cos Z f t
¨
© K Z 2 ¸
m
f ¹
Z n 2
K Z f m
los dos primeros sumandos pueden reemplazarse por
C sen Zn t M
entonces es:
147
Mecánica Racional
Fo
y (t ) C sen Zn t M cos Z f t (3.41)
K Z 2f m
Por lo tanto la respuesta del sistema viene dada por la superposición de dos
movimientos oscilatorios armónicos en la misma dirección, cuyo análisis se realizó en
el Apartado 1.7.3 a). Una de ellas posee la frecuencia natural Zn del sistema y la otra la
frecuencia forzada de la entrada Zf, como se observa en la figura:
El Fenómeno de Resonancia.
Una situación muy importante para analizar desde el punto de vista de la ingeniería, la
Z Zf .
constituye el caso cuando n
K
Zn2
Recordando que m y reemplazando en (3.41):
Fo
y (t ) C sen Zn t M cos Z f t .
ª §Z ·2 º
f
K«1 ¨ ¸ »
«¬ © Zn ¹ »¼
Z f o Zn
Cuando resulta
Fo
of
ª §Z ·2 º
f
K«1 ¨ ¸ »
«¬ © Zn ¹ »¼
148
Capítulo 3- Cinética del Punto Material
1
Se denomina factor de multiplicación a M
ª §Zf ·
2
º
«1 ¨¨ ¸
¸ »
«¬ © Z n ¹ »¼
y su gráfica en función de Zf es:
d2y Fo
2
Zn2 y cos Zn t
dt m
Dada (3.40):
Vo Fo 1
y (t ) y o cosZ n t senZ n t (cosZ f t cosZ n t ) (3.42)
Zn K ª §Z ·
2
º
«1 ¨¨ f ¸
¸ »
«¬ © Z n ¹ »¼
Vo Fo t senZ t
y o cosZ n t senZ n t
f
y (t )
Zn K
(2) §¨¨ Z1 ·¸¸
Zf
Zn
© n ¹
149
Mecánica Racional
Vo Fo (3.43)
y (t ) yo cosZn t sen Zn t t sen Zn t
Zn 2mZn
Se observa que el movimiento de la masa aumenta sin límite a medida que el tiempo
transcurre.
§ Fo ·
¨ t sen Zn t ¸
© 2m Z n ¹
Ejemplo de Aplicación:
El motor a explosión monocilíndrico de la figura está montado sobre un bloque de
cimentación que está apoyado en resortes. Describir la vibración del estado
permanente del sistema si el bloque y el motor tienen una masa total de m = 80 Kg y el
motor cuando está funcionando crea una fuerza de F = 50 sen 2t [N], donde t se mide
en seg. Suponer que el sistema vibra solamente en la dirección vertical, con el
desplazamiento positivo medido hacia abajo, y que la rigidez total de los resortes
puede representarse como K = 2000 N/m.
a) Expresar la ecuación de movimiento correspondiente al estado permanente de la
vibración forzada; y b) Calcular la velocidad rotacional del motor que producirá la
resonancia del sistema.
150
Capítulo 3- Cinética del Punto Material
K Fo Fo / K
y y sen Z f t y p 2
sen Z f t
m m §Z ·
1 ¨ f ¸
© Zn ¹
K 1
? Zn 5
m s
50
a) y p 2000 sen 2t 0,0297 sen 2t [m]
1 25
2
b) Z f 5s 1
d2y dy
m c K y Fo cos Z f t (3.44)
dt 2 dt
dy
(o) Vo (3.46)
dt
151
Mecánica Racional
d 2 y c dy K
y 0
dt 2 m dt m
Solución complementaria:
yc C1 e O1t C 2 e O2t
Las distintas formas de la solución de esta ecuación han sido halladas en vibraciones
libres de un sistema con amortiguamiento.
yp A cos Z f t B sen Z f t
Por consiguiente:
y p AZ f senZ f t BZ f cosZ f t
y p AZ 2f cosZ f t BZ 2f senZ f t
Reemplazando en (3.44):
AZ 2f cosZ f t BZ 2f senZ f t
c
m
> @
AZ f senZ f t BZ f cosZ f t
K
m
>AcosZ f t B senZ f t @ Fo
m
cosZ f t
c K Fo
°° AZ f BZ f m A m
2
m
®
° BZ 2 c AZ K B 0
°¯ f
m
f
m
§ K· §Z f c · Fo
° A¨ Z f ¸ B¨¨ ¸
2
° © m¹ ¸
© m ¹ m
®
° § c · §K 2 ·
° A¨© m Z f ¸ B¨ Z f ¸
¹ ©m ¹
0
¯
De donde:
° A K Z 2f m BcZ f Fo
(3.47)
®
°̄ A cZ f B K mZ f
2
0
152
Capítulo 3- Cinética del Punto Material
A Z f c B K mZ 2f 0
se obtiene
§ Z 2f c 2 · Fo Z f c
B¨ k Z 2
m ¸
¨ K Z2m f
¸ K Z 2f m
© f ¹
Tomando denominador común y como 2 K Z 2f m Z 4f m 2 = K Z 2f m entonces
2
§ Z 2f c 2 K 2 2 KZ 2f m Z 4f m 2 · Fo Z f c
B¨ ¸
¨ K Z 2f m ¸ K Z 2f m
© ¹
Zfc
B Fo
K Z m 2
f
2
Z 2f c 2
B K Z 2f m
Fo Z f c K Z 2f m
A
Zf c >
Z f c K Z 2f m Z 2f c 2
2
@
K Z 2f m
A Fo
K Z m 2
f
2
Z 2f c 2
Se supone que K Z 2f m Z 2f c 2 z 0
2
Fo Fo
B=0 y A ? yp cos Z f t
K Z 2f m K Z 2f m
153
Mecánica Racional
Esta es la solución del problema (3.44), (3.45), (3.46) con c = 0, que lógicamente
coincide con la expresión obtenida al estudiar vibraciones forzadas de un sistema no
amortiguado.
Para el caso en que c z 0, se tiene que la solución particular yp de (3.44) viene dada
por:
K Z 2f m Zfc
y (t ) Fo cosZ f t Fo senZ f t (3.48)
K Z m 2
f
2
Z 2f c 2 K Z m 2
f
2
Z 2f c 2
Fo §K 2 · cZ f
¨ Zf ¸ Fo
m ©m ¹ m2
y (t ) 2
cosZ f t 2
senZ f t
§Z f c · §Z f c ·
2 2
§K 2 · §K 2 ·
¨ Z f ¸ ¨¨ ¸
¸ ¨ Z f ¸ ¨¨ ¸
¸
©m ¹ © m ¹ ©m ¹ © m ¹
y (t )
C cos Z f t I (3.49)
donde:
Fo K
C (3.50)
2
ª §Z f ·
2
º § c Zf ·
2
«1 ¨¨ ¸
¸ » ¨¨ 2 ¸
¸
«¬ © Z n ¹ »¼ © cc Z n ¹
§ ª §Z · 2 º·
c Zf
arc tg ¨ 2 «1 ¨ f ¸ »¸¸ (3.51)
¨ cc Z n «¬ © Zn ¹ »¼¹
©
154
Capítulo 3- Cinética del Punto Material
155
Mecánica Racional
Ejemplo de aplicación:
Un motor eléctrico obligado a desplazarse verticalmente gira a 1470 rpm con un
tornillo prisionero de masa m = 20 g situado a 10 cm del eje de rotación. La masa del
conjunto motor-estructura es M = 50 Kg y está montado sobre un elastómero de
K = 152.000 N/m y c = 2500 Kg/s.
Se desea encontrar:
a) Calcular c / cc y Z f / Zn y visualizar la zona de trabajo en el gráfico del factor de
amplificación M vs Z f /Z n
b) Máximo desplazamiento en el estado estable para un ciclo del movimiento
G
> @
c) Diferencia de fase entre la entrada Fo t y la salida > y t @ . Verificar en el gráfico de
vs Z f /Z n el resultado obtenido
d) Expresión analítica del movimiento estable.
e) Frecuencia natural y seudoperíodo del movimiento.
Solución:
152.000
Zn 55,13 1/ s
a) 50
Zf S 1470
2,79
Zn 30 55,13
cc 2 KM 2 M Zn 5.513 Kg/s
c 2500
cc 0,45
5513
b) Es
y (t ) C cos Z f t I
Fo K
y máx C
2
ª §Z f ·
2
º § c Zf ·
2
«1 ¨¨ ¸
¸ » ¨¨ 2 ¸
¸
«¬ © Z n ¹ »¼ © cc Z n ¹
47,37 1
? y má x 43 P m
152.000
>1 2,79 @
2 2
2 0,45 2,79
2
156
Capítulo 3- Cinética del Punto Material
Zf
2c
cc Zn 20,452,79
c) arctg arctg 20 D 0,35rad
§Z ·
2
1 2,79 2
1 ¨ f ¸
© Z n ¹
e)
2 2
K § c · 152.00 § 2500 ·
E ¨ ¸ ¨ ¸ 49,14 1/s
M © 2M ¹ 50 © 2 50 ¹
T = 2S = 0,13 s
E
Se ha visto que el estado permanente o estable del sistema está descrito por la
ecuación (3.49), en la cual C –dada por (3.50)- representa la amplitud de la oscilación y
I(3.51) el ángulo de fase, es decir, el ángulo de retraso existente entre el movimiento
(respuesta) y la fuerza exterior (entrada).
y = C cos (Zf t - I)
Estos vectores pueden ser representados por la proyección sobre la dirección dada por
y, de los vectores de amplitudes: K C, c Zf C, m Zf2 C y F0:
157
Mecánica Racional
La Ley de Newton, o la ecuación (3. 44), que es lo mismo, requiere que la suma
vectorial de las cuatro fuerzas sea nula en todo instante (polígono funicular cerrado):
6 FV = 0 K C – m Zf2 C - F0 cosI = 0
6 FH = 0 c Zf C - F0 senI = 0
Caso a) Para vibraciones muy lentas (Zf# 0), las FA yFI resultan despreciables, siendo F0
= K C para I = 0.
Aumentando la frecuencia, el vector FA crece y la FI lo hace aún más rápido. El ángulo
de fase no puede seguir siendo nulo, puesto que F0 deberá tener una componente
horizontal hacia la izquierda para equilibrar a c Zf C.
Para este caso, las FI y FA son pequeñas. Es decir, que la FE se emplea casi totalmente
en equilibrar a la FR .
158
Capítulo 3- Cinética del Punto Material
Ky KC cosZ f t I
c y cZ f C senZ f t I
159
Mecánica Racional
y por lo tanto:
FT K C 2 cZ f C 2
C K 2 cZ f
2
FT
TRF
Fo Fo
§ Zf ·2
1 ¨2 c c Zn ¸¹
© c
TRF (3.53)
2
ª §Z ·2 º § c Z f ·2
f
«1 ¨ Z ¸ » ¨ 2 c Zn ¸¹
¬ © n¹ ¼ © c
Zf
El efecto de sobre TRF se muestra la siguiente gráfica:
Zn
160
Capítulo 3- Cinética del Punto Material
Zf
Puede observarse que para una buena aislación debe ser ! 1,41 necesitándose por
Zn
lo tanto un bajo valor de Zn, lo que implica baja rigidez, es decir, un montaje altamente
flexible.
Esto no siempre es aceptable en la práctica donde usualmente es necesaria una cierta
rigidez mínima para satisfacer condiciones de operación. Asimismo, el
amortiguamiento juega un rol importante cuando el sistema debe trabajar cercano a la
Zf
resonancia, pero cumple un efecto adverso para valores Zn ! 1,41 . Es por este
motivo que los montajes antivibratorios se fabrican mayormente con relaciones de
c/cc| 0,15. Asimismo, los estudios de confort en automotores indican que un buen
sistema de suspensión debe poseer valores de fn| 1 Hz y de c/cc| 0,25 (Aparicio
Izquierdo, Teoría de los vehículos automóviles, Madrid, 1995). Estos bajos valores de
amortiguamiento y frecuencia natural contribuyen a lograr una baja transmisibilidad.
Ejemplo de aplicación:
161
Mecánica Racional
Solución:
En el catálogo comercial suministrado por la empresa fabricante de amortiguadores
antivibratoriosVibrachoc (http://www.vibrachoc.es/productos.html) se observa que
para máquinas rotativas girando por encima de 400 rpm se cuenta con amortiguadores
de fn entre 3 y 4 Hz, debiendo usarse de acuerdo al peso del conjunto (3000 Kg) 6
elementos V1136-25 que soportan 420 a 620 Kg cada uno).
Se observa que de acuerdo a estos valores resulta ff / fn = 2 y siendo para este tipo de
elementos c / cc# 0,15, del gráfico de TR resulta TR# 0,4.
m x c y x K y x
o:
m x c x Kx c y K y
c AZ f cosZ f t K AsenZ f t
y tomando:
CZ f Bcos E
K B sen E
resulta:
m x c x K x
B Acos Z f t E (3.44’)
con:
K 2 cZ f 2
B
K
E arctg
cZ f
Por lo tanto, teniendo en cuenta la similitud entre las expresiones (3.44) y (3.44’), se
tendrá:
162
Capítulo 3- Cinética del Punto Material
A K 2 cZ f 2
x(t ) K
2
cosZ f t E I
§ c Zf
>1 Z Zn @
2 2
¨¨ 2
·
¸
¸
© cc Z n
f
¹
§ c Zf ·2
1 ¨2 c Zn ¸¹
C1 © c
TRM (3.54)
A ª §Z 2 º2 2
· § Zf ·
«1 ¨ f Z ¸ » ¨ 2 c c Z ¸
¬ © n¹ ¼ © c n¹
Así, las transmisibilidades de fuerza (3.53) y de movimiento (3.54) resultan ser iguales:
TRF = TRM
Ejemplo de aplicación:
La vibración del piso donde se desea instalar un espectrógrafo de masas delicado es
armónica simple a una frecuencia en el rango 15-60 Hz. El equipo, que debe ser aislado
del movimiento del piso, está sujeto a una pequeña plataforma montada sobre tres
resortes iguales, cada uno soportando la misma carga. El movimiento ocurre en la
dirección vertical. La masa combinada del equipo y la plataforma es de 40 Kg y el
coeficiente de amortiguamiento viscoso de la suspensión es c/cc = 0,2.
Encontrar el valor máximo de la constante de resorte, si la amplitud de la vibración
transmitida tiene que ser menor del 10% de la vibración del piso sobre el rango de
frecuencias dado.
Solución:
Aplicando la ecuación (3.53), es:
2 º2
ª
«¬1 Z
f Z
n »¼ 0,4 Z f Z n
2
> @ ª
2º
100 «1 0,4 Z f Zn »
¬ ¼
es decir:
Z f Zn 4 17,84 Z f Zn 2 99 0
de donde: Z f Zn 4,72
y dado que Zn K m y m = 40 Kg
163
Mecánica Racional
KE = 15.935 N/m
K = KE / 3 = 5300 N/m
La amplitud de la vibración transmitida del piso al equipo será así menor del 10% a
frecuencias superiores a 15 Hz.
164
Capítulo 4- Cinética de los Sistemas Materiales
Los problemas que aquí se plantean son idénticos en esencia a los que aparecen en la
cinética de una partícula. Para llegar a plantear las ecuaciones que resuelven estos
problemas se comenzará por encontrar las expresiones -para el caso de sistemas- de
dos conceptos fundamentales: trabajo y energía cinética, así como también será
necesaria la formulación de algunos teoremas.
4.1. Trabajo elemental de las fuerzas que actúan sobre un sistema material
G G
pero como dri Vi dt
G G
resulta dWi Fi Vi dt
y el trabajo elemental de todas las fuerzas que actúan sobre el sistema será:
N G G
dW ¦ F V dt
i 1
i i
(4.1)
165
Mecánica Racional
N G G
¦F
i 1
int i 0
¦ ¦ G
Fi
G
Fact i ¦ G
Freact i ¦ G
FE i (4.2)
N G G
donde: ¦F
i 1
exti = resultante de las fuerzas exteriores = FE
G G
¦ F >V @
N
G G
dW i 01 Z ri dt
i 1
G G
N N G
G G
¦ Fi V01 dt ¦ Fi Z ri dt
i 1 i 1
(4.4)
G
En el primer sumando, V01 puede sacarse como factor común fuera de la sumatoria,
puesto que al ser la velocidad del centro de reducción, se considera que todos los
puntos del cuerpo se trasladan con ella. Teniendo en cuenta (4.2), se obtiene:
N G G G N G G G
¦ i 01
1
F V dt V01 dt ¦ Fi 1
FE V01 dt
En cuanto al producto mixto del segundo sumando, puede, sin alterarse el orden
cíclico de los vectores, escribirse de la siguiente manera:
N G G G G N G G G N G G
¦ i
F Z ri Z ¦ ri Fi Z ¦ ri FE
i
1 1 1
N
y
¦rG FG
1
i Ei
G
M E ( 01)
166
Capítulo 4- Cinética de los Sistemas Materiales
G G
ei 1 / 2mi Vi 2 1 / 2mi Vi Vi
1 N
e ¦
2i1
mi Vi 2 (4.6)
Reemplazando en (4.6)
N N
e 1/ 2 ¦i 1
G
G 2
mi Voi Vri 1/ 2 ¦
1
G G
mi Vo2 Vr2 2 Vo1 Vri
1 i
167
Mecánica Racional
N N N
1 / 2 Vo2
1 ¦
1
mi 1 / 2 ¦
1
mi Vr2
i
G
Vo1 ¦
1
G
mi Vri
N N
e 1 / 2 m Vo2
1
¦ 1
1 / 2 mi Vr2
i
G
Vo1 ¦
1
G
mi Vri (4.7)
Esta expresión da la energía cinética (que tendrá el mismo valor cualquiera sea el
punto 01 elegido) del sistema material y como se observa está compuesta de tres
sumandos, cada uno de ellos con un importante significado:
El primer sumando, e 1 / 2 m Vo21 , recibe el nombre de energía cinética de arrastre o
de traslación y es la que tendría el sistema en el supuesto que toda la masa estuviera
concentrada en el centro de reducción, siendo generada por la velocidad de éste
último.
N
El segundo sumando e2 ¦1 / 2m V
1
i
2
ri
se denomina energía cinética relativa y está
valor depende del centro de reducción. Esta e3 puede anularse si se toma como centro
G G
de reducción a un punto fijo del sistema (si lo hubiera), así resultaría V01 0 , por
lo tanto e3 = 0, siendo en este caso e1 también nula.
Pero el caso más importante de anulación de e3 es cuando se toma como centro de
reducción el baricentro “G” del sistema, lo que da lugar al teorema de König: “La
energía cinética de un sistema material cualquiera es en cada instante igual a la
energía cinética que corresponde al baricentro supuesto que en él está concentrada
toda la masa, más la energía cinética que le corresponde al sistema en su movimiento
relativo al baricentro”.
Para demostrarlo, es conveniente recordar que se define como centro de masa de un
sistema de N puntos materiales al punto cuyo vector posición está dado por:
G
rG
¦ mi rGi xG
¦ mi xi
; yG
¦ mi yi
; zG
¦ mi zi
¦ mi m m m
G
En efecto entonces, si 01 es el centro de reducción, luego rG es el vector posición de G
respecto de 01.
G G
Si 01 { G rG 0
N
Por otro lado: m rGG
¦
1
G
mi ri
168
Capítulo 4- Cinética de los Sistemas Materiales
G N G
derivando: m drG
dt ¦ 1
mi
dri
dt
G
G G drG G G G
Siendo rG 0 será 0 y por lo tanto ¦ mi Vr 0
dt i
G G G G
Vi Vo1 Z ri
N N
e 1/ 2 ¦
1
G
mi Vi2 1/ 2 ¦m VG
1
i o1
G G 2
Z ri
N N N
1/ 2 ¦
1
mi Vo2
1
1/ 2 ¦ 1
G G 2
mi Z ri 1 / 2 ¦m 2 VG
1
i o1
G G
Z ri
N N
1 / 2 m Vo2
1
1/ 2 Z 2
¦ 1
G 2 G G
mi Z ri Vo1 Z ¦1
G
mi ri
169
Mecánica Racional
Z rGi 2 1 ri sen D 2 d i2
G 2
En consecuencia mi Z ri mi d i2 es el momento de inercia de la masa mi
respecto al eje Ẑ . Así
N
G
¦ m Zˆ r
2
i i I ZZ 01
1
G
resulta ser el momento de inercia del sistema material rígido respecto al eje Z
pasante por 01.
G G
En el tercer sumando, se tiene ¦ mi ri rG y por lo tanto, la expresión general de la
energía cinética para un cuerpo rígido es:
1 1 G G G (4.8)
e mVo21 I ZZo Z 2 mVo1 Z rG
2 2 1
En un cuerpo rígido siempre es posible anular e3. Una forma es usar como centro de
reducción a un punto fijo, lo que será posible si el cuerpo está en movimiento de
rotación o de rotación instantánea, en este caso es:
1 G G
e I ZZo Z 2
2 1
V o1
0
Pero también se puede anular e3 tomando el baricentro como centro de reducción; en
G G
este caso rG 0 e3 0 y:
1 1
e mVG2 I ZZG Z 2
2 2
170
Capítulo 4- Cinética de los Sistemas Materiales
G G
donde Vo1 es la velocidad de una de las partículas del sistema y Vri es la de todas las
otras respecto de ella. Así:
N N
G G
G
¦ G
¦ G N
Si se tomara como centro de reducción al baricentro G del sistema, se tendría (01 { G):
N
G
Q
G
m VG ¦1
G
mi Vri( G )
Esta última expresión implica que la cantidad de movimiento total del sistema
material es la que tendría su baricentro en el supuesto de que toda la masa estuviese
concentrada en él.
Para los sistemas rígidos la (4.9) constituye la expresión más directa para obtener la
cantidad de movimiento.
Siendo para
G una partícula: G
G
K( o1)i ri mi Vi
Se definirá como momento cinético de un sistema de partículas respecto de un punto
01 a la suma de los momentos cinéticos de cada partícula:
N
G
K( o1) ¦ 1
G G
ri mi Vi
N N
G
K( 01) ¦
1
G G
ri mi Vo1 ¦ 1
G G
ri mi Vri
N N
G
V01 ¦ 1
G
mi ri ¦1
G G
ri mi Vri
171
Mecánica Racional
N
?
G
K( o1)
G G
rG m Vo1 ¦ G G
ri mi Vri (4.10)
G
1 G
donde rG es la posición del baricentro respecto al centro de momentos 01 y Vo1 es la
velocidad de ese punto.
El primer sumando del término de la derecha en (4.10) expresa que una parte del
momento cinético respecto del punto 01 sería el que tendría toda la masa como si ésta
estuviese concentrada en el punto G y con la velocidad de 01. Recibe el nombre de
momento cinético de arrastre u orbital.
El segundo sumando es debido a las velocidades relativas a 01 y se denomina momento
cinético relativo o intrínseco.
Si se ubicara la terna sobre el baricentro G, tomándoselo como centro de momentos,
G G
se tendría rG 0 y:
N
G
K( G ) ¦ 1
G G
ri mi Vri (4.10’)
que es el momento cinético relativo al baricentro. Por lo tanto, respecto del baricentro
el momento cinético total y el relativo son iguales, no existiendo momento cinético
orbital. Los movimientos de las partículas con respecto al baricentro son los que
generan la existencia del momento cinético relativo al mismo. Se verá ahora qué forma
adopta la (4.10) para el caso específico de un cuerpo rígido. En este caso se tomará
como centro de reducción el origen de la terna en un punto 01 del mismo, donde ahora
G G
estarán aplicados los vectores característicos V01 yZ , siendo:
G G G G G G
Vi V01 Vri V01 Z ri
luego, reemplazando:
N N
G
K( o1) ¦ 1
G G
ri mi V01 ¦
1
G
> G G
mi ri Z ri @
El primer sumando resulta
G
G
rG m V01
Para resolver el segundo sumando se aplica la fórmula de Gibbs:
172
Capítulo 4- Cinética de los Sistemas Materiales
G G G G G G G G G
a b c a c b a b c
obteniéndose
G G G N
G G G G G G
K (01) rG m V01 ¦ mi ª¬ ri ri Z ri Z ri º¼ (4.11)
1
con:
G G
ri ri ri2 xic 2 yic 2 zic 2
por lo que:
G G
ri2 Z xc i c Z cy j c Z zc k c
ri 2 Z
o, trabajando en forma matricial:
ª1 0 0º ª Z xc º
» « »
G
ri2 Z «
xic 2 yic 2 zic 2 « 0 1 0» « Z cy»
G
ri2 1 Z (4.12)
«¬ 0 0 1¼» «¬ Z c »
z ¼
Donde 1 es la matriz unidad.
Y el producto:
rGi ZG rGi xic ic yic j c zic kc xic Z xc yic Zic zic Z zc
Donde X i es una matriz cuadrada de 3 x 3 originada por el producto binario del vector
G
posición ri por sí mismo.
Luego, introduciendo (4.12) y (4.13) en (4.11), resulta:
G
K o1
G G N
G
rG mVo1 ¦ mi ri 21 X i Z
1
(4.14)
173
Mecánica Racional
d i2 yic 2 z ic 2
?mi d i2 I xcxic
174
Capítulo 4- Cinética de los Sistemas Materiales
que representa el tensor de inercia del cuerpo rígido respecto de la terna de origen 01.
En el Apéndice I -Tensores Cartesianos- se repasan conceptos físicos y de Álgebra
Lineal concernientes a estas magnitudes tensoriales.
Notar que I o1 está referido a la terna con respecto a la cual se expresa la velocidad
G
angular Z .
De esta forma, la expresión (4.14) queda:
G G G G
K o1 rG m Vo1 I o1 Z (4.16)
Según qué punto se tome como centro de reducción, se tendrán diversos casos:
a) Si se toma un punto fijo 01 del cuerpo, será:
G G G G
Vo1 0 K o1 I o1 Z (4.16a)
G G G
ri rG Ui (4.18)
G G
N G N G G G N G G G G G
K o1 ¦ G i i ¦ U i mi Vi
r
1
m V
1
rG ¦ mi Vi K G rG Q K G
1
luego: G G G
G
K o1 K G rG Q (4.19)
175
Mecánica Racional
ªi c j c k cº ª I xx I xy I xz º ªZ x º
G « » « »« »
K o1 m « xG yG zG » « I yy I yz » «Z y »
«¬Vo1x Vo1y Vo1z »¼ «¬ Sim I zz »¼ «¬Z z »¼
Ejemplo de aplicación
G
Sea la rueda B de la figura sometida a la rotación Z1 de dirección i3 que rueda sin
deslizar sobre el plano i1 , i2 . Se considera que la masa de
la varilla A01 es
despreciable. Los parámetros mecánicos y geométricos del problema son:
a = 20 cm
r = h = 10 cm
T(t) = 20 S t
G
Z1 20S i3c 1/s
g
U = densidad = 8
cm 3
e = espesor = 1 cm
Nota: Despreciar e para el
cálculo de los momentos de
inercia
176
Capítulo 4- Cinética de los Sistemas Materiales
Se desea determinar:
a) Cantidad de movimiento del disco.
b) Energía cinética del disco tomando como centro de reducción: b.1) al punto 01, b.2)
al punto P, b.3) al punto 0.
Solución:
Por comodidad, se referirán los vectores respecto a la terna fija a la varilla A01 , ^iic`.
G G § Kgm ·
Q m VG 31,58 iˆ2c ¨ ¸
© s ¹
b.1)
G G
e1 1 / 2mVo21 , VG = V01 ?
2
1 2 m
e1 2,513kg 4S 2 198,5 Joule.
2 s
G
Z 20 S iˆ3c 40 S iˆ1c (1/ s ) ? Z 2 20S 40 S
2 2
19.740 1/ s 2
1
I11c m R 2 0,01256kg m 2
2
1
c
I 22 m R 2 0.00628kg m 2
4
1
c
I 33 m R 2 0,00628kg m 2
4
I12c I13c I 23c 0.
177
Mecánica Racional
ª 0,01256 0 0 º
« »
? I o1^i c` « 0.00628 0 »
i
«¬ Sim. 0,00628»¼
Z2
tg M 2M 63,42º
Z1
I ZZo1 I 3cc3cc 2
c31 I11c c32
2
c c33
I 22 2
c
I 33
c31
cos iˆ3cc ; iˆ1c cos 900 63, 420 0,89431
c32 cos 90º 0.
c33 cos 63,42º 0,44744.
? I ZZo1 0,011302kg m 2
1
? e2 0,011302 19740 111,5 Joule
2
además:
G G
e3
0 porque 01 { G rG 0
luego
e = e1 + e2 + e3 # 310 Joule
b.2)
1
e1 mV p2
2
G G G G
Vp
Z 2 2 Riˆ3c 800S iˆ2c cm / s siendo VP
G G
VP ' Z rP ' P
2
1 2 m
e1 2,513kg 8S 793,7 Joule.
2 seg 2
178
Capítulo 4- Cinética de los Sistemas Materiales
1
e2 I ZZp Z 2
2
I ZZp I ZZG md 2
0,011302 2,5130,0894
2
I ZZp 0,03138Kgm 2
1
e2 0,03138 19740 # 310 Joule.
2G
e3
G G
> @
mVp Z rG 2,513kg 8S iˆ2 20S iˆ3c 40S iˆ1c 0,1iˆ3c = -793,7 Joule
Así,
e e1 e2 e3 # 310 Joule
b.3) Se deja para el alumno demostrar que si se toma como centro de reducción el
punto 0, se obtiene: e1 = 198,40 J; e2 = 905,00 J y e3 = - 793,40 J, resultando e # 310 J.
G
Verificar que la distancia entre los ejes paralelos a Z que pasan por 0 y 01 resulta igual
a 0,18 m.
G G G G
c) K o K G rG Q
G G § m·
Q m VG 2,513 4S iˆ2c 31,58 ¨ kg ¸ iˆ2c
© s¹
G G
rG Q
a iˆ1c R iˆ3c 31,58 iˆ2c
2
m
6,316 iˆ3c 3,158 iˆ1c ( Kg )
s
G § m2 ·
? K 0 6,316iˆ3c 3,158iˆ1c 1,578iˆ1c 0,394iˆ3c ¨ Kg ¸
© s ¹
G § m2 ·
K 0 4, 73 iˆ1c 6, 71 iˆ3c ¨ Kg ¸
© s ¹
A continuación se enunciará una serie de teoremas que son válidos para la cinética de
los sistemas en general, aunque se particularizarán aquí para los cuerpos rígidos.
179
Mecánica Racional
G
G G dVi G G G
FI i FEi mi con FEi = FAi FRi
G G dt
donde FI i y FE i son las resultantes de las fuerzas interiores y exteriores que actúan
sobre la partícula.
N N N G
¦ 1
G
FE i ¦1
G
FI i ¦
1
mi
d Vi
dt
N G G
para un cuerpo rígido, es ¦F
1
Ii 0 (por cuanto éstas actúan de a pares, con el mismo
N
módulo, colineales y opuestas); luego ¦ 1
G
FE i
G
FE . También es:
N G N G
¦ 1
mi
d Vi
dt
d
dt ¦
1
G
mi Vi
dQ
dt
Luego: G
G dQ
FE (4.20) Ecuación de Newton
dt
G G G G G
Como Q m VG , se observa que si FE 0 resulta VG constante, lo que indica que el
movimiento del baricentro no puede variar si la resultante de las fuerzas exteriores es
nula.
Si se ha tomado una terna de referencia en movimiento para referir los vectores se
tendrá:
G
Q Qxiˆ1 Qy ˆj1 Qz kˆ1
y G
dQ d Qx d Qy d Qz d i d j1 d k1
i1 j1 k1 Qx 1 Qy Qz
dt dt dt dt dt dt dt
G
Los tres primeros sumandos expresan la variación de Q respecto de la terna móvil
G
ªd Qº
como si ésta estuviese fija y se denotará: « » . Los tres siguientes sumandos
¬dt¼ rel
d i1 G
pueden expresarse teniendo en cuenta las fórmulas de Poisson: : i1 , etc.,
dt
180
Capítulo 4- Cinética de los Sistemas Materiales
G G G
como : Q , siendo : la velocidad angular de la terna móvil (que si es solidaria al
G
cuerpo coincidirá con la Z del mismo).
Por lo tanto:
G
G ª dQ º G G
FE « » :Q (4.20’)
¬ d t ¼ rel
G N
G G
Recordando que: K o1 ¦r m V
1
i i i
G
donde ri es el vector posición de la partícula genérica i, con origen en el centro de reducción 01
G N G N G
¦ G
¦
d K o1 d ri G d Vi
Luego: mi Vi ri mi
dt dt dt
1 1
G G G
d ri
Pero: Vi Vo1 (velocidad relativa a la terna móvil)
dt
G
Notar que aquí Vo1 es la velocidad absoluta de 01 porque a diferencia de los apartados
precedentes, en este caso no se lo ha utilizado como centro de reducción para aplicar
la forma impropia de la ley de distribución de velocidades.
Reemplazando
G N N N G
¦ G G
¦ G G
¦
d K o1 G d Vi
Vi mi Vi Vo1 mi Vi ri mi
dt dt
1 1 1
Donde el primer producto es nulo por cuanto se trata de vectores colineales.
G N G
G G
¦
d K o1 G d Vi
V o1 Q ri mi (4.21)
dt dt
1
Como se vió,
G
G d Vi G G G G G
ri mi
dt
ri Fi ri FEi FIi
G
G G G G G d Vi (4.22)
? ¦ ri FEi ¦ ri FIi ¦ ri mi
G
G
G dt
M E o1 M I o1 0
181
Mecánica Racional
En (4.22), MEo1 es el momento de todas las fuerzas exteriores respecto del centro
G
momentos 01; y M I o1 es el momento de las fuerzas interiores respecto de 01, que
por actuar éstas de a pares y contrarias, en un rígido resulta nulo.
Reemplazando (4.22) en (4.21)
G
G d K o1 G G (4.23) Ecuación de Euler
M E o1 Vo1 Q
dt
G
G G G G d K o1
Si el punto 01 es fijo u 01 { G Vo1 Q 0 y M E o1 (4.23bis)
dt
G
Si la terna respecto de la cual está referido K o1 no es inercial, es decir que rota con
G
: , se tendrá:
G
K o1 K x iˆ1 K y ˆj1 K z kˆ1
y
G
d K o1 dKx ˆ dKy ˆ dKz ˆ d iˆ d ˆj d kˆ
i1 j1 k1 K x 1 K y 1 K z 1
dt dt dt dt dt dt dt
que se transforma en:
G G
d K o1 ª d K o1 º G G
« » : K o1
dt ¬« d t ¼» rel
G
G ª d K G º G G
?M EG « » : K G
«¬ d t »¼ rel
Como
G antes, si la terna elegida para expresar los vectores está adherida al cuerpo, será
G
:= Z
Las ecuaciones (4.20’) y (4.23’) constituyen las denominadas ecuaciones cardinales de
la cinética o del movimiento de los cuerpos rígidos. En el espacio, de estas dos
ecuaciones vectoriales se obtienen seis escalares que relacionan entre sí a los
parámetros intervinientes en la mayoría de los problemas de la cinética. Esto permite:
o bien determinar las seis coordenadas que fijan su posición en función del tiempo,
cuando se conoce en cada instante el sistema de fuerzas que actúan sobre el cuerpo; o
bien encontrar el sistema de fuerzas actuantes si lo que se conoce es el movimiento.
d Wi d ei
y extendiendo a todo el sistema:
dW de (4.24)
Esta expresión dice que el trabajo elemental de todas las fuerzas que actúan sobre el
sistema es igual a la variación de la energía cinética del mismo. Debe notarse que en el
trabajo W intervienen fuerzas interiores y exteriores teniendo presente que si bien
G G
¦ FIi 0 , no siempre es ¦WFI = 0, resultando nulo sólo para el caso de los
i
sistemas rígidos. Por este motivo, la aplicación de la Gecuación (4.24) sólo se hace
posible cuando se conoce el comportamiento de las FI , tal como es el caso de los
rígidos o de los sólidos elásticos. Conviene recordar que las expresiones (4.5) y (4.8)
nos dan la forma final de la (4.24) para un sistema rígido:
§G G G G·
dW © FE Vo1 M E Z ¹ dt
o1
ª1 1 G G G º
de d « mVo21 I ZZo1 Z 2 mVo1 Z rG »
¬2 2 ¼
Igualando estas expresiones resulta una ecuación diferencial de aplicación directa.
4.5.4 Potencia
Utilizando la definición de Potencia Instantánea del Apartado 3.5.2, surge de la
expresión anterior del trabajo que para un cuerpo rígido, es
dW G G G G
P FE Vo1 M E o1 Z
dt
En el caso que no exista traslación, resulta:
dW G G
P M E o 1 Z
dt
Esta última expresión es de gran aplicación en los bancos de prueba de motores, en los
cuales a determinada velocidad de rotación se mide el par motor para calcular la
Potencia entregada.
183
Mecánica Racional
G
d Si
¦ mi d t vector constante
184
Capítulo 4- Cinética de los Sistemas Materiales
G G G
Si FE 0 Q vector constante
G G G G G
Si MFE o1 0 para Vo1 0 K o1 vector constante
G G
K G I G Z vector constante
G G
K o1 I o1Z vector constante
185
Mecánica Racional
a) La terna móvil estará fija al cuerpo y el punto fijo 0 será el origen del sistema de
coordenadas por lo que serán nulos los momentos de las reacciones en 0; también uno
de sus ejes coincidirá con el de rotación, luego:
G
Z Z k 1
b) El centro de gravedad del cuerpo se ubicará dentro de uno de los planos
coordenados, eliminándose así una coordenada (en el caso de la figura: yg). En otras
palabras, la proyección del punto G sobre el plano i, j determina la dirección del
eje i1 .
Es conveniente referir todos los vectores al sistema móvil, puesto que conociendo M =
M(t) es muy sencillo pasar al sistema fijo por transformación de coordenadas. Por otra
parte, la distribución espacial de la masa del cuerpo permanecerá constante respecto
de este sistema fijo a él, resultando constante el tensor de inercia.
Se supondrá que el cuerpo está sujeto a un sistema de fuerzas activas cuya resultante
es:
N
¦ i 1
G
FAi
G
FA ¦ Fx i1 ¦ Fy j1 ¦ Fz k1
G G G G G
y se hallarán las reacciones giratorias R Xo ;RYo ;RZ o ;R X1 ; R Y1 aplicando las ecuaciones
cardinales de la dinámica. Estas reacciones se denominan estáticas, en el caso que el
cuerpo se analice en estado de reposo y dinámicas cuando se generan como
consecuencia del movimiento (desbalanceo, desalineación, etc). A la composición
entre estáticas y dinámicas se las denomina reacciones totales.
Para ello, por lo expresado anteriormente, se tomarán momentos respecto del punto
fijo 0. Así:
G
M FE o ¦ M x RY 1 h iˆ1 ¦ M y
R X1 h ˆj1 ¦ Mz kˆ1
G G G G
K o K G rG Q
186
Capítulo 4- Cinética de los Sistemas Materiales
G G G G G G
con K G I GZ o K o I o Z rG mVo
Como se dijo, el tensor de inercia se referirá a los ejes móviles. Así resultará, de (4.16.
a):
ª* * I xz º ª 0 º
G «* * I » « 0 » I Ziˆ I Z ˆj I Z kˆ
K o « yz » « » xz 1 yz 1 zz 1
«¬* * I zz »¼ «¬Z »¼
donde los asteriscos indican las componentes del tensor que no hacen falta ser
calculadas por no intervenir en el producto matricial (debido a las componentes nulas
G
de Z ). Por lo tanto:
G
ª d K o º
« » I xz Z iˆ1 I yz Z ˆj1 I zz Z kˆ1
¬« d t ¼» rel
y:
G G
Z K o I yz Z 2 iˆ1 I xz Z 2 ˆj1
Habiendo encontrado todos los términos que componen las ecuaciones cardinales, se
está en condiciones de proseguir con su aplicación.
De la primera de ellas:
G G G G G G
FE FA FR Q > Q @ rel ZG Q Expresión de Newton
¦ Fx R X o R X1 m X G Z 2
°
°
®¦ Fy RY o RY1 m X G Z (4.26.a)
°
°¦ Fz RZ o 0
¯
La última de las tres igualdades expresa que las reacciones sobre k no están afectadas
por el movimiento de acuerdo con las condiciones que se han elegido.
De la segunda de las ecuaciones cardinales (ecuación de Euler):
G
G ª d K o º G G
M o « » Z K o
¬« d t ¼» rel
187
Mecánica Racional
¦ Mx RY1 h I xz Z I yz Z 2
°°
®¦ My R X1 h I yz Z I xz Z (4.26.b)
2
°
°¯¦ Mz I zz Z
Estas seis últimas expresiones (4.26 a y b) permiten analizar el caso general en el cual
G
la velocidad angular Z varía con el tiempo (por ejemplo, en el arranque de una
turbina, motor eléctrico, etc).
Algunos casos particulares son:
G
a) Caso en que Z cte
Se tendrá:
¦ Fx R X o R X1 m xG Z 2
°°
®¦ Fy RYo RY1 0
°
°¯¦ Fz RZ o 0
¦ Mx RY1 h Iyz Z 2
°°
®¦ My R X1 h I xz Z
2
°
°¯¦ Mz 0
Esta última igualdad a cero es una condición necesaria y no una ecuación, puesto que
si existiera momento en k, habría aceleración angular J, y ?Z Z (t ) z cte
Se nota que a pesar de considerar Z = cte, algunas reacciones están influenciadas por
el movimiento.
b) Si la única fuerza activa fuese el peso (siempre con Z = cte), se tendría:
G
FA ¦ ¦
m g k Fx Fy 0
¦ Fz mg
Mx Mz 0
My m g xG
y las ecuaciones cardinales:
RX o RX 1 m X G Z 2
°
° RYo RY1 0
°
° m g RZo 0
®
° RY1 h I yz Z 2
°
°m g X G RX1 h I xz Z
2
°
¯ 0 0
G
c) Cuando el centro de masas del cuerpo G está sobre k, siendo mg la única fuerza
actuante, se tiene:
188
Capítulo 4- Cinética de los Sistemas Materiales
R R 0
° Xo X1
° RYo RY1 0
°
® m g R Z o 0
°
° RY1 h I yz Z
2
°
¯ R X1 h I xz Z
2
I yz Z 2
RY1 RYo RY1 RZo mg
h
Nótese que a pesar de no existir fuerzas activas en las direcciones x1 e y1, sí existen
reacciones de vínculo en esas mismas direcciones. Por lo tanto, éstas surgen como
consecuencia del movimiento y no de la carga. RZo no depende del movimiento ya que
sólo equilibra al peso.
Por otra parte, siendo:
R X o R X1 y RYo RY1
surge claramente que las reacciones dinámicas generan cuplas.
d) Sistema de fuerzas activas cuya resultante (incluido el peso) pasa por el punto 0
se tendrá:
Mx = My = Mz = 0
y por lo tanto:
¦ Fx R X o R X1 0
¦ Fy R Yo RY1 0
¦ Fz RZ o 0
189
Mecánica Racional
RY1 h I yz Z 2
R X 1 h I xz Z 2
G
Se observa que el vínculo en 0 absorbe a la resultante FA .
Las reacciones RX1 y RY1 dependen exclusivamente del movimiento (puesto que si Z =
0, éstas no existen) y por ello reciben
G el nombre de reacciones dinámicas. Así mismo,
RXo, RYo, RZo dependen de las FA y del movimiento denominándose reacciones
totales. La diferencia vectorial entre las reacciones totales y las dinámicas, da como
resultado las reacciones estáticas. Es sencillo observar que para que no se produzcan
esfuerzos dinámicos en el sistema material rígido, deberán anularse Iyz e Ixz para lo cual
el eje de rotación debe coincidir con un eje principal de inercia.
Así por ejemplo, si la chapa triangular de la figura está en reposo no ejerce ningún
empuje lateral sobre los soportes, debido a que su centro de gravedad está sobre k ;
sólo existirá la reacción estática RZo. Se dice que la chapa se halla estáticamente
equilibrada.
Sin embargo, y dado que z1 no es eje principal de inercia, es Iyz z 0 y las reacciones
G
RYo y RY1 serán proporcionales a Z2 cuando la chapa comience a girar con Z Z k1 .
Estas reacciones giran con la chapa originando desgaste en los cojinetes y vibraciones
en las estructuras que los soportan. Un mecanismo está dinámicamente equilibrado si
las fuerzas antes y durante el movimiento son las mismas.
En el caso de la chapa en forma de paralelogramo de la figura de la izquierda, se
observa que si no hay movimiento la misma está en equilibrio estático (con reacciones
en 0 vertical y sin reacción en 01).
190
Capítulo 4 – Cinética de los Sistemas Materiales
PROBLEMA EJEMPLO
Se aplicarán las ecuaciones de Newton (4.20’) y de Euler (4.23’), tomando como terna
G
^`
de referencia la iˆicc solidaria a la chapa, es decir, rotando con Z alrededor de ^ i3cc`
G
G G ªd Qº G G
¦ FE « » Z Q
¬ d t ¼rel
G
G ª d K ( A) º G G G G
¦ M FE( A ) « » Z K ( A) V A Q
«¬ d t »¼ rel
191
Mecánica Racional
Solución:
Ecuación de Newton:
G G G G K
Q mVG 0 ¦ FE 0
³Z t dt
t
M t
0
G G
W mg mg (cos M iˆ1cc sen M )iˆ2cc
G G G G
¦ FE W RA RB
m g (cos I iˆ1cc sen I iˆ2cc) RA1 iˆ1cc
G
cc cc cc cc
RA2 i2 RA3 i3 RB1 i1 RB2 i2 0
ˆ ˆ ˆ ˆ
Luego:
m g cos M R R 0
° A1 B1
® m g sen M R A R B 0
2 2
°
¯ R A3 0
Ecuación de Euler:
G G
Ahora K A I A Z , y de tablas:
192
Capítulo 4 – Cinética de los Sistemas Materiales
I11
m 2
12
a ; I 22
12
m 2
a b 2 ; I 33
m 2
12
b , luego es:
ª I 11 0 0 º
^`
I G iˆi «0
« I 22 0 »»
«¬ 0 0 I 33 »¼
c22 cos 0º 1
c33 cos J 11
luego:
0
I11c c112 I11 c122 I 22 c132 I 33
c
I12 0
c
I13 c11 c31 I11 c13 c33 I 33
c
I 21 0
c
I 22 I 22
c
I 23 0
c
I 31 c31 c11 I11 c13 c33 I 33 c
I13
c
I 32 0
193
Mecánica Racional
2 2
c
I 33 c31 I11 c33 I 33
ª a 2 b2 º
cc
I11 c m«
I11 »
¬ 4 ¼
cc
I12 c m00 0
I12
0
cc
I13 c m 1 / 2
I13 a 2 b2 0 c
I13
cc
I 21 0
§ a 2 b2 ·
cc
I 22 c m¨
I 22 ¸
© 4 ¹
cc
I 23 0
cc
I 31 c 0
I 31 c
I 31
cc
I 32 0
cc
I 33 c
I 33
Así resulta el tensor de inercia en A referido a los ejes ^ iicc` , y aplicando Euler:
ª I "11 0 I "13 º ª0 º
G « 0 I"
K A
« 22 0 »» ««0 »» I 13cc Z iˆ1cc I 33
cc Z iˆ3cc
«¬ I "31 0 I "33 »¼ «¬Z »¼
Recordar que al trabajar numéricamente los momentos de inercia centrífugos deben
cambiarse de signo al introducirlos en el tensor de inercia.
G
ª d K ( A) º
« » I13cc Z iˆ1cc I 33
cc Z iˆ3cc
«¬ d t »¼ rel
G G
Z K A I13cc Z 2 iˆ2cc
G G G
V A Q 0
194
Capítulo 4 – Cinética de los Sistemas Materiales
G G G G G
1/ 2 a 2 b 2
1/2
¦ M FE A rG mg rB RB iˆ3cc mg cos I iˆ1cc sen I i2cc
a 2 b2
1/2
iˆ3cc RB1 iˆ1cc RB2 iˆ2cc
1/2 ª§ m g sen I · § m g cos I · º
a 2
b 2 «¨ RB2 ¸ iˆ1cc ¨ RB1 ¸ iˆ2cc»
¬© 2 ¹ © 2 ¹ ¼
e igualando en componentes:
§ m g sen M ·
cc Z
° I13 ¨
© 2
R B2 ¸ a 2 b 2
¹
1/ 2
°
° § m g cos M ·
cc Z 2 ¨
® I13
© 2
R B1 ¸ a 2 b 2
¹
1/ 2
°
° I 33
cc Z 0
°
¯
De la última Z 0 Z CTE M Z t , lo que es lógico puesto que no hay ni
G
momento de las FE según i3 ni una cupla aplicada en esa dirección.
Así, estas últimas expresiones quedan:
m g sen M
° R B2 0
° 2
® m g cos M cc Z 2
I13
° R B1
°
¯
2
a 2 b 2 1/ 2
y por lo tanto:
m g sen M
° R B2 2
°
® cc Z 2
I13 m g cos M
° RB
°
¯
1
a 2 b2
1/ 2
2
Retomando las ecuaciones cartesianas de Newton podemos calcular las reacciones en
A:
m g cos M cc Z 2
I13
°RA
° 1
®
2
a 2 b 2 1/ 2
° m g sen M
° R A2
¯ 2
195
Mecánica Racional
Para M = 0 se tiene
mg
R B1 cte
2
R B2 0
mg
R A1 cte
2
R A2 0
Para M = S/2
R B1 cte
mg
R B2
2
R A1 cte
mg
R A2
2
Puede observarse claramente que mientras las reacciones estáticas generadas en este
caso por el peso permanecen inalteradas con el movimiento, la cupla dinámica rota
con la chapa.
Se vio que cuando el cuerpo tiene un punto fijo 0, las ecuaciones cardinales de la
cinética toman la forma:
G
G ªd Qº G G
FE « » : Q
¬ d t ¼rel
G
G ª d K o º G G
M E o « » : K o
«¬ d t »¼ rel
196
Capítulo 4 – Cinética de los Sistemas Materiales
en las cuales el subíndice rel indica que las derivadas deben efectuarse como si la terna
G
en rotación 0 i1 , j1 , k1 permaneciera inmóvil, y : es la rotación absoluta de esta
última terna respecto de la 0 i, j , k , que cuando los ejes i1 , j1 , k1 son solidarios al
G
cuerpo es la rotación absoluta de éste último Z .
Se verá como se aplican en este caso las ecuaciones cardinales.
Para el estudio que se propone encarar se tomará una terna de referencia móvil
unida solidariamente al cuerpo, la cual, para mayor simplicidad coincidirá con los ejes
principales de inercia del mismo con origen en el punto fijo 0. En la Figura se ilustra el
elipsoide de inercia de tal cuerpo (Ver Apéndice I, Tensores Cartesianos – A1.4).
ª I xx 0 0º
I o ^iˆ1 , ˆj1 ,kˆ1 ` «0
« I yy 0 »»
«¬ 0 0 I zz »¼
Además, el cuerpo tendrá tres grados de libertad, es decir que su posición quedará
expresada por los tres ángulos de Euler (\, T, M) y la velocidad de cada punto se
conocerá por el vector rotación:
G G
:{Z Z x i1 Z y j1 Z z k1
G
Luego el vector K o será
197
Mecánica Racional
G G
K o I o Z I xx Z x iˆ1 I yy Z y ˆj1 I zz Z z kˆ1
y por lo tanto
G
ª d K o º
« » I xx Z x iˆ1 I yy Z y ˆj1 I zz Z z kˆ1
¬« d t ¼» rel
G G
El producto Z K o resultará
i1 I zzZ yZ z I yyZ yZ z j1 I xxZ xZ z I zzZ zZ x k1 I yyZ yZ x I xxZ xZ y
de manera que reemplazando en la ecuación de momentos
M I xx Z x Z y Z z I zz I yy
° Ex (0)
°
®M E y ( 0 ) I yy Z y Z x Z z I xx I zz (4.27)
°
°¯M E z ( 0 ) I zz Z z Z x Z y I yy I xx
G G
Estas son las llamadas ecuaciones de Euler. En el caso más general FE y ME serán
G
funciones de t, de r y de la velocidad de los puntos en que actúan.
G
Las componentes de Z podrán expresarse también en función de los ángulos de
Euler:
G G G G
Z \ M T
198
Capítulo 4 – Cinética de los Sistemas Materiales
luego:
Estas tres ecuaciones junto con las (4.27) son seis ecuaciones diferenciales de primer
orden en Z x ,Z y ,Z z ,\ ,T ,M que permiten calcular los ángulos de Euler en función del
tiempo, que es el objetivo que se persigue (aunque su integración analítica sólo es
posible en algunos casos particulares).
Las ecuaciones de Newton (4.20’) resultan útiles en este problema para determinar las
reacciones en el punto fijo 0; se tendrá:
G G G G
Q m VG m Z rG
199
Mecánica Racional
G
y FE Rox ¦ Fx i1 Roy ¦ Fy j1 Roz ¦ Fz k1
donde Rox, Roy, Roz son las reacciones (incógnitas) en 0 y ¦ Fx , ¦ Fy , ¦ Fz , son las
fuerzas activas (datos).
Las ecuaciones de Euler (4.27) no son en general integrables, pero resuelven algunos
casos particulares, de los cuales los dos más importantes para el ingeniero son:
El anillo exterior posee una rótula que le permite el giro sobre el eje k .El anillo
interior está sujeto con una rótula al exterior que le permite el giro en i. El giróscopo
está vinculado al anillo interior con una unión que le permite girar sobre su eje de
simetría. Los dos aros se encuentran en el plano S 1 i, k cuyo versor normal es ĵ ; y
el giróscopo en el plano coordenado iˆ, ˆj de versor normal k̂ .
Sacando al giróscopo de esta posición, se lo puede “seguir” con los ángulos de Euler. El
primer movimiento será el ángulo de precesión \ alrededor de k (único eje fijo), lo
G
que hará que los aros se trasladen al plano S 2 i c, k , ; mientras el giróscopo continúa
en el plano (i c, j c),
200
Capítulo 4 – Cinética de los Sistemas Materiales
El segundo movimiento será el ángulo de nutación T, por el cual el aro interior gira
alrededor de iˆc pasando a encontrarse sobre el plano S 3 iˆc,kˆ1 ; mientras que el
giróscopo se ubica en el plano (iˆc, ˆj cc) normal a k̂1 .
201
Mecánica Racional
G
Z \ k T i1 M k1
o, dado que k está en el plano j1 , k1 y es k cos T k1 sen T j1 ; resultan
G
Z T i1 \ sen T j1 \ cos T M k1 (4.28)
G
: T i1 \ sen T j1 \ cos T k1 (4.29)
Las expresiones (4.28) y (4.29) dan los vectores rotación expresados en función de los
ángulos de Euler y sus derivadas y están referidos a la terna móvil.
Zx T
Zy \ sen T
Zz \ cos T M (4.30)
Se convendrá ahora en llamar Izz al momento de inercia principal, respecto del eje de
rotación propia y se llamará I a los momentos de inercia respecto de los otros dos ejes,
los cuales serán iguales entre sí.
donde:
ªI 0 0 º ªT º
G G « » « »
K o I GZ « I 0 » «\ sen T »
«¬ sim I zz ¼ ¬\ cos T M »¼
» «
derivando:
202
Capítulo 4 – Cinética de los Sistemas Materiales
G
ª d K o º
« »
ITiˆ1 I \ senT \ TcosT ˆj1 M \cosT \ TsenT I zz kˆ1 (4.31)
¬« d t ¼» rel
y:
ªi j1 k1 º
G G «1 »
: K o «T \ sen T \ cos T »
« »
¬I T I \ sen T I zz M \ cos T ¼
Mx
I T \ 2 senT cosT I zz \ M \ cosT senT
My I\ senT 2 I\TcosT I TM \ cosT
zz (4.33)
Mz M \cosT \TsenT I zz
M
° x
> I zz M \ cos T I \ cos T @ \ sen T
® My 0
°
¯ Mz 0
o,
G
M o > I zz Z z I \ cos T @ \ sen T i1
203
Mecánica Racional
Lo que significa que al aplicar un momento respecto al punto 0 según una dirección
normal al plano que determinan la rotación propia y la precesión o sea k y k1 se
logra poner al giróscopo en precesión estable.
Un caso muy común en ingeniería es el denominado precesión estable normal en el
que T 90º , de donde sen T = 1 y cos T = 0, por lo que
Mx I zz M \
G
ó M o I zz M \ i1
vectorialmente G
G G
M o I zz \ M
I zz \ ˆj1 M kˆ1
I zz \ M i1
204
Capítulo 4 – Cinética de los Sistemas Materiales
G
En el caso de un rotor industrial girando con spin M , la introducción de una precesión
normal produce la aparición del momento que tenderá a volcar el rotor generando
reacciones en contrafase en los cojinetes.
b) Trompo o giróscopo pesado. El cuerpo -que es simétrico con respecto a un eje (es
decir, cuyo elipsoide de inercia es de revolución: cono, cilindro, pirámide regular,
disco, etc.)- posee un punto fijo (0) y además la única fuerza actuante es el peso, que
está aplicada en G que pertenece al eje de revolución (0 z G ) . Es el caso denominado
de Lagrange-Poisson.
205
Mecánica Racional
G
El M o referido a la terna móvil para este caso es:
G
M o l k1 P k l k1 P sen T j1 cos T k1
G
? M o P l sen T i1
Mx P l sen T
es decir: ®
¯ My Mz 0
P l sen T
I T \ 2 sen T cos T I zz \ M \ cos T sen T ½
°°
0 I \ sen T 2 I \ T cos T I T M \ cos T
¾ (4.34)
°
0 I zz M \ cos T \ T sen T °¿
d d
0 I zz
dt
M \ cos T I zz Z
dt z
206
Capítulo 4 – Cinética de los Sistemas Materiales
de donde resulta
M \ cos T Zz cte.
es decir que en este movimiento la velocidad angular total del trompo en la dirección
k1 permanece constante.
Se analizará en qué condiciones se producirá una precesión estable; para esto es
T cte T o
Así resulta que el trompo tiene una precesión alrededor del eje k con un ángulo
constante de inclinación To.
De la primera de las (4.34) prosigue que para esta condición:
I zz I cos T o \ 2 I zz M \ P l 0
> M 4 P l I zz I cos T o @
2 2 1/ 2
I zz M r I zz
\1,2 (4.35)
2 I zz I cos T o
como puede observarse, para un mismo ángulo de nutación To, hay en general dos
valores para la rapidez de precesión \ . Estas se denominan precesión “rápida” y
“lenta” según se use el signo + ó - en la expresión (4.35) respectivamente.
Para cualquier movimiento “posible” del trompo, \ debe ser un número real, lo que
implica que en la (4.35) debe ocurrir siempre la condición:
M 4 P l I zz I cos T o t 0
2 2
I zz
207
Mecánica Racional
G G
K o I o Z # I zz Z z k1
que representa
G la velocidad del extremo
del vector K o . Por lo tanto, el
G
momento del peso M o trata de
G
desplazar el extremo de K o en la
G G
dirección de M o y como K o debe
K
coincidir con Z , el eje de rotación inicia
un movimiento de precesión, el cual da
lugar al principio de aplicación de los giróscopos.
I zz M \ Pl
Pl
\
I zz M
208
Capítulo 4 – Cinética de los Sistemas Materiales
Cuando la única fuerza externa que actúa sobre un cuerpo es causada por la
gravitación, el movimiento general del cuerpo se conoce como movimiento libre de
par motriz. Es característico de los planetas, satélites artificiales y proyectiles
(despreciando los efectos de fricción con el aire). Se supondrá que el cuerpo tiene su
masa distribuida axialmente simétrica. El origen de los ejes coordenados i1 , j1 , k1 se
localiza en el centro de masa G. Se nota que I xx I yy I .
Arrojando el cuerpo hacia el aire, la suma de momentos respecto al baricentro es cero,
lo que significa -ecuación (4.23 bis)- que el momento cinético del cuerpo es constante:
G G G
M G 0 KG cte
G
Tomando k { K G y al eje j1 que quede en el plano formado por los ejes k1 y k. El
G
ángulo de nutación T será el formado por k y k1 , y como K G está en este plano, o
sea en el j1 , k1 , se tendrá:
G
K G KG sen T j1 KG cos T k1 . (4.36)
ª I 0 0 º ªZ x º
G G «0 I 0 » «Z »
K G I GZ « »« y »
«¬0 0 I zz »¼ «¬Z z »¼
209
Mecánica Racional
°Z 0
° x
° KG
®Z y sen T
° I
° KG
°Z z I zz
cosT
¯
o:
G KG K
Z sen T j1 G cos T k1
I I zz
T 0
KG
\ sen T sen T
I
KG
\ cos T M cos T
I zz
°T cte
°
° KG
®\ (4.38)
° I
° I I zz
°¯M I I zz
KG cos T
G G
Así, para un movimiento como éste, se observa que el ángulo T (entre K G y M )
G G
permanece constante durante el movimiento. Además K G , \ y M se conservan
durante todo el tiempo con lo que éste resulta en una precesión estable.
Combinando en (4.38) se tiene la siguiente relación entre el spin y la precesión:
I I zz
M \ cos T (4.39)
I zz
El cuerpo precesa alrededor del eje k (fijo) mientras que gira alrededor del k1
(móvil).
210
Capítulo 4 – Cinética de los Sistemas Materiales
El cono espacial (que define la precesión) está fijo, ya que la precesión lo está,
mientras que el cono del cuerpo gira alrededor de la superficie exterior del cono
espacial sin deslizar.
Como el espín es una función de los momentos de inercia I e Izz del cuerpo puede
ocurrir que Izz > I y el espín resulta negativo, siendo la precesión positiva.
211
Capítulo 5- Dinámica Analítica
5. DINÁMICA ANALÍTICA
Desde que Newton enunció formalmente como axiomas las “leyes del movimiento”
(Principia Mathematica, Londres, 1686) la Mecánica fue expandiéndose como ciencia a
través de numerosos aportes tanto desde la Física como desde la Matemática.
El método de la Dinámica Analítica o de Lagrange presentado un siglo después
(Mécanique Analytique, París, 1788), constituye un punto de vista completamente
diferente al de sus predecesores. Más que seguir el movimiento de cada masa
individual de un sistema, Lagrange introdujo coordenadas generalizadas (qj) en igual
cantidad que los grados de libertad del sistema (k). Suele decirse que el esfuerzo
colectivo de organización de la Mecánica del Siglo XVIII culminó con Lagrange, es decir,
con la casi perfección formal de la Mecánica Clásica. A partir de aquí se entra en el
terreno didáctico.
El método es aplicable a un gran conjunto de problemas de partículas y de cuerpos
rígidos que abarca desde los más sencillos hasta los más complejos. Las ventajas de
este procedimiento sobre los métodos convencionales, por las razones que se indican
a continuación, son de gran importancia.
En gran parte, el método de Lagrange reduce todo el campo de la estática, de la
dinámica de partículas y de la dinámica de los cuerpos rígidos a un solo procedimiento
que implica las mismas etapas básicas, independientemente del número de masas
consideradas, del tipo de coordenadas empleadas, del número de restricciones o
vínculos sobre el sistema y de que éstos y el marco de referencia estén o no en
movimiento. Por lo tanto, los métodos especiales se reemplazan por un método
general único.
Pueden utilizarse una gran variedad de coordenadas generalizadas. Es decir, que las
ecuaciones de Lagrange son válidas para cualquier sistema de coordenadas (inerciales
o una combinación de inerciales y no inerciales) que resulten adecuadas para
representar la configuración del sistema. Estas conducen directamente a las
ecuaciones del movimiento en cualquier sistema de coordenadas que se escoja,
aunque en ningún caso a la integración de las mismas. No se requiere introducir
primero el método vectorial formal y luego hacer la traducción a las coordenadas
deseadas.
Las fuerzas de restricción, para los casos de restricciones holónomas lisas, quedan
eliminadas automáticamente y no aparecen en las ecuaciones de Lagrange. Como se
definió en el Capítulo 2, en las restricciones holónomas las ecuaciones de ligadura o de
vínculo pueden expresarse en función de las coordenadas solamente, o de las
coordenadas y el tiempo y en ellas el número de coordenadas generalizadas es igual al
número de grados de libertad. La eliminación de las fuerzas de restricción por métodos
convencionales puede presentar dificultades enormes.
El método de Lagrange está basado en gran parte en magnitudes escalares: energía
cinética (e), energía potencial (p), trabajo virtual elemental (GW), y en muchos casos, la
213
Mecánica Racional
potencia (P). Todas ellas pueden expresarse generalmente sin ninguna dificultad en
cualquier sistema de coordenadas adecuado. Desde luego, la naturaleza vectorial de la
fuerza, la velocidad, la aceleración, etc, debe tenerse en cuenta en el estudio de los
problemas de dinámica. Sin embargo, las ecuaciones de Lagrange, basadas en las
anteriores magnitudes escalares tienen en cuenta en forma completa y automática
estas magnitudes vectoriales sin necesidad de recurrir a métodos vectoriales formales.
Por complejo que sea un sistema, los términos de la ecuación del movimiento de
Lagrange consisten en las componentes adecuadas de la fuerza y la aceleración
expresadas en el sistema de coordenadas elegido.
Las ideas básicas empleadas en la deducción de las ecuaciones de Lagrange son
sencillas. La aplicación de las ecuaciones de Lagrange a los problemas prácticos resulta
sencilla aún cuando se trate de problemas complicados. Excepto en los problemas muy
elementales, el procedimiento es en general más sencillo y requiere invertir menos
tiempo que otros métodos.
Sobre la base de los estudios de Lagrange, en 1853 el matemático irlandés William R.
Hamilton introdujo el concepto de las “funciones fuerza” (energía potencial con signo
negativo). Sus ecuaciones canónicas transforman el sistema de k ecuaciones
diferenciales de segundo orden de Lagrange para el caso de sistemas holónomos
conservativos, en un sistema de 2k ecuaciones de primer orden, con la ventaja de
incluir una única función, la hamiltoniana (H).
Por cuestiones de tiempo y espacio, este Capítulo pretende presentar en forma
resumida los principios básicos de la dinámica de Lagrange y ofrecer al alumno
suficiente entrenamiento en las técnicas –físicas y matemáticas- necesarias para
aplicar sus ecuaciones. Debido a que esta teoría abarca a la totalidad de la Mecánica,
no es posible incluir los aspectos más sofisticados de la misma, alentando al estudiante
con inquietudes a consultar la excelente y profusa bibliografía especializada. En cada
apartado se incluyen varios problemas ejemplos para facilitar la fijación de conceptos.
Asimismo, dado que la aplicación de las ecuaciones de Lagrange conduce en forma
directa a las ecuaciones diferenciales del movimiento, gran parte de ellas de muy
complicada integración, se deja para los alumnos la resolución de las mismas,
incentivando la utilización de la herramienta computacional mediante diferentes
software de ingeniería.
Sea una partícula P de masa m sometida a una fuerza externa aplicada, resultante de
las cargas activas y reactivas, de componentes Fx, Fy, Fz. Multiplicando cada una de las
componentes cartesianas escalares de la ecuación del movimiento de Newton para m
(3.1) por el desplazamiento virtual respectivo, se tiene:
Fx Gx mxGx
Fy Gy myGy
Fz Gz mzGz
214
Capítulo 5- Dinámica Analítica
x1 = l1 sen q1
y1 = l1 cos q1
x2 = l1 sen q1 + l2 sen (q1 + q2)
y2 = l1 cos q1 + l2 cos (q1 + q2)
Así,
x1 = x1 (q1)
Ecuaciones de las
y1 = y1 (q1)
restricciones para
x2 = x2 (q1, q2)
dos grados de
y2 = y2 (q1, q2)
libertad
En general, se tendrá:
Estas son las ecuaciones de las restricciones para el sistema de k grados de libertad del
cual forma parte m. No se estudiarán en este curso coordenadas rectangulares que
sean funciones del tiempo además de las qj, como así tampoco ejes coordenados y
restricciones móviles.
215
Mecánica Racional
wx 1
^@
wq1
`@ wy 1
wq1
wz 1
@
wq1
wx wy wz § wx wy wz ·¸
Fx Fy Fz Gq1 m¨¨ x y z Gq1 (5.2)
wq1 wq1 wq1 © wq1 wq1 wq1 ¸¹
a) Miembro izquierdo:
dW F dr Fds cos F , dr
Considerando
+!
$
%
* $
!
!&
necesariamente sobre la trayectoria dada, a este desplazamiento se lo llama “virtual”:
{^`|!“trabajo virtual elemental”
!
+!
*
+
216
Capítulo 5- Dinámica Analítica
b) Miembro derecho:
Partiendo de la identidad
d § wx · wx d § wx ·
¨ x ¸ x x ¨¨ ¸
dt ¨© wq1 ¸¹ wq1 dt © wq1 ¸¹
es:
wx d § wx · d § wx ·
x ¨¨ x ¸¸ x ¨¨ ¸
wq1 dt © wq1 ¹ dt © wq1 ¸¹
dx wx wx
x q1 q 2 ...
dt wq1 wq 2
!
!!+$$
!+&¥^¦¥&j es función solamente de qj y el
último sumando se anula por cuanto q 2 z f( q1 ), ya que son independientes.
wx wx
? (procedimiento de “eliminación de puntos”)
wq1 wq1
217
Mecánica Racional
wx d § wx · § wx ·
x ¨ x ¸ ¨ x ¸
wq1 dt ¨© wq1 ¸¹ ¨© wq1 ¸¹
§ wx · w § x 2 · § wx · w § x 2 ·
¨¨ x ¸¸ ¨ ¸ y ¨¨ x ¸¸ ¨ ¸
© wq1 ¹ wq1 ¨© 2 ¸¹ © wq1 ¹ wq1 ¨© 2 ¸¹
Por lo tanto
wx d ª w § x 2 · º w § x 2 ·
x « ¨ ¸ » ¨ ¸ , y por similitud
wq1 dt « wq1 ¨© 2 ¸¹ » wq1 ¨© 2 ¸¹
¬ ¼
wy d ª w § y ·º w § y 2 ·
2
y « ¨ ¸» ¨ ¸,
wq1 dt ¬ wq1 ¨© 2 ¸¹¼ wq1 ¨© 2 ¸¹
wz d ª w § z 2 ·º w § z 2 ·
z « ¨ ¸¸» ¨¨ ¸¸
wq1 dt ¬ wq1 ¨© 2 ¹¼ wq1 © 2 ¹
1 d ª w 2 º 1 w 2 ½
GWq1 Q1Gq1
m®
« x y z »
2 2
x y 2 z 2 ¾Gq1
¯ 2 dt ¬ wq1 ¼ 2 wq1 ¿
pero como : e
1
2
m x 2 y 2 z 2 resulta
d § we · we
Q1 ¨ ¸
dt ©¨ wq1 ¸¹ wq1
d §¨ we ·¸ we
Qj ; j 1,2,..., k (5.5)
dt ¨© wq j ¸¹ wq j
que son las ecuaciones de Lagrange para un punto material sometido a cualquier clase
de fuerzas.
Estas k ecuaciones de segundo orden para cada grado de libertad j, darán la ecuación
de movimiento del punto material que corresponde a cada grado de libertad.
218
Capítulo 5- Dinámica Analítica
? e
1 2 1
2
mv
2
m x 2 y 2 z 2
y las fuerzas generalizadas Fx, Fy, Fz { Fq j
219
Mecánica Racional
Así : e
1
2
m x 2 y 2
pero y 2bxx e
1
2
m x 2 4b 2 x 2 x 2
1 2
2
mx 1 4b 2 x 2
Aplicando la ecuación (5.5)
we we d §¨ we ·
wq j wx
mx 1 4b 2 x 2 ; dt ¨© wq j
¸
¸
mx 1 4b 2 x 2 8mb 2 x 2 x
¹
we we
4mx 2 xb 2
wq j wx
Por lo tanto;
mx 1 4b 2 x 2 4mxxb 2 Fx Qj j 1
Para hallar Fx se supone que el eje y es vertical y que la única fuerza activa es la
atracción de la gravedad. Imprimiendo a la esf
+!
$ ¢^
$ $#
+ ¢`+ !
!
$%, como
se observa en la Figura!
%
+!
*
#
¨ @ - $`
ª« !
k !
k `@%^^
Fx 2mgbx
mx 1 4b 2 x 2 4mx 2 xb 2 2mgbx
G
A los fines ilustrativos, se considera una fuerza f aplicada a la esfera mediante un hilo
en el plano x-`` ! !
k con la horizontal se mantiene constante
(los cosenos directores c11y c12 son constantes).
G
f
f cos M iˆ senM ˆj
cos M @11
senM @12
Para un desplazamiento pequeño, se tiene según
(5.4):
GW mgGy c12 fGy c11 fGx
y como Gy 2bxGx, se tiene
GW - 2mgbx c11 f 2bxc12 f Gx y
Fx 2mgbx c11 f 2bxc12 f 2bxc12 f mg c11 f
K
Capítulo 5- Dinámica Analítica
mx 1 4b 2 x 2 4mx 2 xb 2 2bxc12 f mg c11 f
G
v rrˆ rTTˆ
e
1
2
m r 2 r 2T 2
Aplicando la ecuación (5.5):
we
mr
we wr
wq j we
mr 2T
wT
d § we ·
· ¨ ¸ mr
d §¨ we ¸ dt © wr ¹
dt ¨© wq j ¸ d § we ·
¹ ¨ ¸ mr T 2mrTr
2
dt © wT ¹
we
we mrT 2
wr
wq j we
0
wT
mr mrT 2 Fr
?® 2
¯mr T 2mrrT FT
221
Mecánica Racional
de donde: Fr k r r0 mg cos T
y similarmente: FT mgrsenT .
wx wy wz
b) Por aplicación directa de la ecuación (5.2). Es Q j Fx Fy Fz
wq j wq j wq j
Las componentes x e y de las fuerzas sobre m, son:
Fx k r r0 senT
Fy mg k (r r0 ) cos T
x rsenT
Pero como:
y r cos T
wx
Resulta para :
wq j
wx wy wz
senT ; cos T ; 0
wr wr wr
wx wy wz
r cos T ; rsenT ; 0
wT wT wT
222
Capítulo 5- Dinámica Analítica
!
&£^`Z
y la de la barra:
eB
1
2m V A 2 1 I AT 2 2mV A T rG
2 2 A
eB
1
2
1 1
2 3 ¬
2mVA 2 2m 2a T 2 2mVAiˆ Tkˆ a ª senT iˆ cos T ˆj º
2
¼
mVA 2
4 2 2
3
ma T
¬
2mVAiˆ ª Ta senT ˆj cos T iˆ º
¼
4 2 2
mVA 2 ma T 2amTVA cos T
3
1 2 § 4 ·
e mx m¨ x 2 a 2T 2 2aT.x cos T ¸
2 © 3 ¹
we
mx 2mx 2maT cos T
we wx
wq j we 8
m a 2T 2amx cos T
wT 3
223
Mecánica Racional
d §¨ we ·
¸
dt ¨© wq j ¸
¹
d § we ·
¨ ¸ 3mx 2maT cos T 2amT senT
2
dt © wx ¹
d § we · 8
¨ ¸ ma 2T 2axm cos T 2amxTsenT
dt © wT ¹ 3
we we we
0; 2amT.xsenT
wq j wx wT
Si
F u
wu ˆ wu ˆ wu ˆ
F i j k
wx wy wz
dW F dr du
y como du - dp,
wp wp ˆ wp ˆ
es : F iˆ j k p
wx wy wz
Se ha visto en (5.2) que la fuerza generalizada para la variación de la coordenada qj, es:
224
Capítulo 5- Dinámica Analítica
wx wy wz
Qj Fx Fy Fz
wq j wq j wq j
Puede entonces escribirse
§ wp wx wp wy wp wz · wp
Qj ¨ ¸¸
¨ wx wq wy wq wz wq wq j
© j j j ¹
en (5.5) :
d § we · we wp
¨ ¸¸
dt ¨© wq j ¹ wq j wq j
O sea:
d § we · w
¨ ¸¸ e p 0 (5.5´)
dt ¨© wq j ¹ wq j
L=e–p
d §¨ wL · wL
¸ 0 j=1,2,…,k (5.6)
dt ¨© wq j ¸ wq
¹ j
Que son las ecuaciones de Lagrange para un punto material bajo la acción de
fuerzas conservativas.
Solución:
Sabiendo que:
1
2
m r 2 r 2T 2
e
la energía potencial acumulada en el resorte será:
1
k r r0
2
pr
2
y la energía potencial gravitatoria respecto a un nivel de referencia coincidente
con el eje x:
p g mgr cos T
225
Mecánica Racional
1
? k r r0 mgr cos T
pr p g
2
p
2
Así , L
1
2
1
2
m r 2 r 2T 2 k r r0 mgr cos T
2
wL d wL
mr · mr
wL wr d §¨ wL ¸ dt wr
wq j wL dt ¨© wq j ¸ d wL
mr 2T ¹ 2mrrT mr 2T
wT dt wT
wL
wL mrT 2 k r r0 mg cos T
wr
wq j wL
mgrsenT
wT
m.r 2
T 2S
k1.S1 k 2 .S 2
2 2
1
e .m.(T.r ) 2
2
226
Capítulo 5- Dinámica Analítica
En tanto que la energía potencial estará compuesta por la gravitacional más la de los
resortes:
pG m.g .(h r.T )
1
p r1 .k1 .(h L1 s1 .T ) 2
2
1
pr 2 .k 2 .(h L2 s 2 .T ) 2
2
wL
m.T.r 2
wT
d § wL · 2
¨ ¸ m.T .r
dt © wT ¹
wL
[m.g .r k1.s1 (h L1.s1 T .s1 ) k2 .s2 .(h L2 T .s2 )]
wT
(k .s k .s )
2 2
T 1 1 2 2 2 T 0
m.r
227
Mecánica Racional
k1 .s1 k 2 .s2
2 2
Z 2
m.r 2
por lo tanto:
k1 .s1 k 2 .s2
2 2
Z
m.r 2
2.S 2.S
T
Z k1 .s1 k 2 .s 2
2 2
m.r 2
o sea:
m.r 2
T 2S
k1.S1 k 2 .S 2
2 2
n ° d § we · wei ½° n
¦ ® ¨¨
i 1 °̄ dt wq
i ¸
¸ wq ¾° ¦Q ji j 1,2,....., k
© j ¹ j ¿ i 1
n
pero ¦e
i 1
i es la energía cinética de todo el sistema, y
n
¦Q
i 1
ji Q j1 Q j 2 ! Q jn es la fuerza generalizada resultante que elimina
n G
las fuerzas interiores del sistema, puesto que ¦ Fint i
i 1
0
d § we · we
Entonces: ¨ ¸¸ Qj j 1, 2,..., k (5.5.a)
dt ¨© wq j ¹ wq j
que son las ecuaciones de Lagrange para un sistema general de puntos materiales,
conservativo o no.
d § wL · wL
y: ¨ ¸¸ 0 j 1, 2,..., k (5.6.a)
dt ¨© wq j ¹ wq j
228
Capítulo 5- Dinámica Analítica
Teniendo en cuenta que el sistema tiene dos grados de libertad (¡@), en el esquema
simplificado de la figura se adoptan como coordenadas generalizadas q1 @^ y q2 @
1
La energía cinética del carro es: ec .m1 .x 2
2
1 1 G G G
.m2 .Vo I 0I 2 m2Vo I roG
2
y la energía cinética de la varilla: ev
2 2
229
Mecánica Racional
ev
1
2
.m2 .x 2
11
23
A
>
m2 A 2I 2 m2 xiˆ Ikˆ senIiˆ cos Iˆj
2
@
1 1 A
ev .m2 .x 2 m2 A 2I 2 m2 xI cos I
2 6 2
m1 2 m2 § 2 1 2 2 ·
L .x ¨ x A I AxI cos I gA cos I ¸
2 2 © 3 ¹
wL w ª m1 2 m2 § 2 1 2 2 ·º
« .x ¨ x A I AxI cos I gA cos I ¸»
wx wx ¬ 2 2 © 3 ¹¼
wL
m1 m2 x m2 AI cos I
wx 2
d wL
m1 m2 x m2 AI cos I I 2 senI
dt wx 2
y
wL
0
wx
d wL m2 § 2 2 ·
¨ A I Ax cos I AxIsenI ¸
dt wI 2 ©3 ¹
wL
wI
m2
2
A xI.senI g .senI
K
Capítulo 5- Dinámica Analítica
2) Una partícula P de masa m se encuentra sobre una mesa lisa horizontal, sujeta a una
cuerda larga e inextensible de longitud l que pasa por un agujero liso en su superficie.
El otro extremo de la cuerda soporta una partícula Q de masa km, que cuelga
libremente. Se impulsa a la partícula P a partir del reposo formando ángulos rectos con
la cuerda, con una velocidad 8ag cuando está a una distancia a del agujero.
Encontrar las ecuaciones del movimiento para el sistema de dos partículas.
+@¡$
e
1
2
1
m r 2 r 2T 2 kmr 2
2
1 1
? L mk 1r 2 mr 2T 2 kmgr
2 2
wL wL wL
m(k 1)r mr 2T
wq j wr wT
d §¨ wL ·
¸ ${¡¢| r mr 2T 2mrrT
dt ¨© wq j ¸
¹
wL wL wL
kmg mrT 2 ; 0
wq j wr wT
231
Mecánica Racional
.I 0 .Z 2 M .Vo m1 .Vo V1 m2 .Vo V2
1 1 1 2 1 2
e r et
2
e
2 2 2 2
1 1 1 1
e r et .I 0 .Z 2 M .Vo m1 .Vo V1 m2 .Vo V1
2 2 2
e
2 2 2 2
y1
Z
r
y las velocidades:
232
Capítulo 5- Dinámica Analítica
Vo y
V1 y1
2
1 § y1 · 1 1 1
.I 0 .¨ ¸ M . y 2 m1 . y y1 m2 . y y1
2 2
e
2 © r ¹ 2 2 2
p p g pe
Energía potencial gravitatoria:
y1 S .r y2 A
y2 C y1
@l –r
pg M .g. y m1.g. y y1 m2 .g. y C y1
ª!
!
k
@– p, será:
ª 1 § y1 · 2 1 1 1 º
L « .I .¨ ¸ M . y 2 m1 . y y1 m2 . y y1 »
2 2
¬« 2 © r ¹ 2 2 2 ¼»
ª1 º
« .k . y y0 M .g. y m1 .g. y y1 m2 .g. y C y1 »
2
¬2 ¼
L L^y1 , y, y1 , y `
Habrá tantas ecuaciones como grados de libertad tenga el sistema. En este caso se
tienen dos grados de libertad.
233
Mecánica Racional
wL
M . y m1. y y1 m2. y y1
wL wy
wq j wL I
. y1 m1. y y1 m2 y y1
wy1 r2
d § wL ·
¨ ¸ My m1 y1 m2 (
y y
y1 )
d § wL · dt © wy ¹
¨ ¸
dt ¨ wq ¸ d § wL · I
© j ¹ ¨ ¸
y1 m1 (
y
y1 ) m2 (
y
y1 )
dt © wy1 ¹ r2
wL
k . y y 0 M .g m1.g m 2 .g
wL wy
wq j wL
m1.g m 2 .g
wy1
d §¨ wL · § wL
¸¨
·
¸ 0
dt ¨© wq j ¸ ¨ wq
¹ © j
¸
¹
Resulta
My m1 (
y
y1 ) m2 (
y
y1 ) [k .( y y0 ) Mg m1 g m2 g ] 0
§ I0 ·
¨ 2 .
y1 m1 (
y
y1 ) m2 ( y1 ) ¸ m1.g m2 .g 0
y
©r ¹
Reordenando:
y. M m1 m2
y1 (m1 m2 ) k . y y0 g . M m1 m2 0
§ I ·
y1. ¨ 2 m1 m2 ¸ g . m1 m2 0
y (m1 m2 )
©r ¹
234
Capítulo 5- Dinámica Analítica
d §¨ we ·¸ w
e p @ ¡
dt ¨© wq j ¸¹ wq j
Introduciendo nuevas variables \j relacionadas con las qj mediante
we
\j (5.7)
wq j
las ecuaciones (5.5´) toman la forma:
d\ j w (e p ) (5.5”)
dt wq j
que constituyen un sistema de k ecuaciones lineales que vinculan los k parámetros \j y
permiten calcular las q j en función de las pj, de las qj y de t, después de reemplazar en
las ecuaciones de Lagrange las q j por los valores obtenidos.
A estas nuevas variables \j suelen denominárselas “ímpetus” del sistema, y están
directamente vinculadas al concepto de “cantidad de movimiento” cuando se trabaja
con las ecuaciones de Newton.
Si además se define ahora una función H dada por la ecuación
k k
H ¦\
j 1
j q j (e p ) ¦\
j 1
j q j L
wH d\ j
\ j (5.9)
wq j dt
Las ecuaciones (5.8) y (5.9) forman 2k ecuaciones diferenciales de primer orden que
reemplazan a las ecuaciones de Lagrange que, como se ha visto, son de segundo
orden. Se denominan Ecuaciones de Hamilton, siendo H la función de Hamilton o
“función hamiltoniana”.
235
Mecánica Racional
Como medio para analizar la mayor parte de los problemas aplicados, este método no
es tan conveniente como el de Lagrange. Su utilidad crece en problemas que se tratan
por técnicas de transformación, perturbación, mecánica estadística y cuántica.
\1 \2 \3 1 \1 \ 2
2
\ 2 2
\1 \ 2 \ 3 m( 2 2 32 ) mgz
m m m 2 m m m
1
2m
\ 1 2 \ 2 2 \ 3 2 mgz
las ecuaciones (5.8) en este caso, dan:
wH \1 wH \2 wH \3
x; y ; z
w\ 1 m w\ 2 m w\ 3 m
y las (5.9)
wH wH wH wH wH wH
0 \ 1 ; 0 \ 2 ; mg \ 3
wq1 wx wq 2 wy wq3 wz
!
¡@
®
$!n.
Tomando la derivada de las (5.8) con respecto al tiempo y eliminando \j de las (5.9) se
obtienen las ecuaciones usuales del movimiento. Es decir,
d wH \ 1 d wH \ 2 d wH \ 3
( ) x ; ( ) y ; ( ) z
dt w\ 1 m dt w\ 2 m dt w\ 3 m
236
Capítulo 5- Dinámica Analítica
mx 0
my 0
mz mg
Estas últimas tres ecuaciones son las relaciones halladas por aplicación directa de las
ecuaciones de Newton primero y de Lagrange después. Por lo tanto, es evidente que
las ecuaciones de Hamilton en este problema no presentan ninguna ventaja.
q1 @^&2 @`
2
H ¦\ j q j L \ 1 m
\ 2
m 2 ©m
m¨¨ 12 22
m
¸ k x2 y2
¸ 2
j 1 ¹
1
2m
\ 21 \ 2 2 k x 2 y 2
1
2
Luego, de (5.8):
wH \1 wH \2
x; y
w\ 1 m w\ 2 m
y de (5.9):
237
Mecánica Racional
wH wH
kx \ 1
wq1 wx
wH wH
ky \ 2
wq 2 wy
Ahora:
d § wH · \ 1
¨ ¸ x \ 1 mx
dt ¨© w\ 1 ¸¹ m
d § wH · \ 2
¨ ¸¸ y \ 2 my
dt ¨© w\ 2 ¹ m
y reemplazando en (5.9)
k
mx kx x x
m
k
my ky y y
m
Estas últimas son las ecuaciones diferenciales de dos movimientos oscilatorios
armónicos en x e y. El movimiento del punto P resulta entonces de la composición de
k
dos movimientos oscilatorios armónicos ortogonales de pulsación Z , dando
m
lugar a la correspondiente figura de Lissajous.
238
Apéndice 1
APÉNDICE 1
G G G G G G
T x y T x T y x;yHV
G G
y T O x O T x (O número real)
T’ = R T RT
ó:
t ijc cim c jn t mn (A.1)
Donde cim y c jn resultan ser los cosenos directores entre los ejes girados y los
originales. La expresión A1 está escrita en la notación indicial de Einstein, donde se
entiende que los subíndices que se repiten en un mismo miembro suman,
permaneciendo fijos los que no varían; así por ejemplo:
c
t13 c11 c31 t11 c11 c32 t12 c11 c33 t13
c12 c31 t 21 c12 c32 t 22 c12 c33 t 23
c13 c31 t 31 c13 c32 t 32 c13 c33 t 33
239
Mecánica Racional
gráficamente:
que cuando se lo refería al sistema {01 , iˆ1 ; iˆ2 ; iˆ3 } tenía como expresión
Tensor de inercia.
Tomando ahora la expresión (4.15) se tiene que el tensor de inercia de un cuerpo
referido a un sistema coordenado con origen en 01 es:
240
Apéndice 1
ª I I xy I xz º
« xx »
I o1 « I yx I yy I yz »
« »
¬ I zx I zy I zz ¼
Es muy fácil demostrar que este tensor de inercia es simétrico, por cuanto:
N N
I xy ¦ 1
mi xic yic y I yx ¦ 1
mi yic xic
Luego, en la aplicación de la expresión (4.16) del momento cinético, para el caso del
tensor de inercia respecto a ejes rotados, será:
a) Para un elemento con subíndices iguales (momentos de inercia o axiales)
I cyy c
I 22 2
c21 I11 c21 c22 I12 c21 c23 I13 ......... c23
2
I 33
I zy
c I 32
c c31 c21 I11 c31 c22 I12 c31 c23 I13 ........ c33 c23 I 33
Notar que los momentos de inercia axiales son directamente elementos del tensor,
pero los momentos centrífugos (o productos de inercia) son los elementos del tensor
con signo cambiado.
I xG xG md x2 I yG yG m( y´G z´G )
2 2
I xcxc
(A.2)
y
I xcyc I xG yG m x´G y 'G
donde:
{0, x’, y’, z’} son ejes paralelos a los ejes baricéntricos (xG, yG, zG)
dx, dy, dz son las distancias entre los ejes xG, yG, zG y los x’, y’, z’, respectivamente.
x´G, y’G, z´G, son las distancias de m a los planos coordenados denotados por sus
versores normales i, j y k.
241
Mecánica Racional
N
m ¦ 1
mi masa del sistema
Así también:
I yG yG md y2 I yG yG m( x´G z´G )
2 2
I ycyc
I zG zG md z2 I yG yG m( x´G y´G )
2 2
I zczc
I xczc I xG zG m x'G z 'G
I yczc I yG z G m y ' G z ' G
242
Apéndice 1
También se supondrán conocidos todos los momentos de segundo orden del cuerpo
con respecto al sistema o base ^ii `. Esto es, se conoce el tensor de inercia
ª I 11 I 12 I 13 º
I o «
« I 22 I 23 »»
«¬ Sim I 33 »¼
Tomando la recta como eje i1c de un nuevo sistema de coordenadas ^iic` y calculando
el momento de inercia del cuerpo respecto de este nuevo eje,
2
c
I11 c11 c11 I11 c11 c12 I12 c11 c13 I13 c12 c11 I12 c12 I 22 "
2 2 2
= c11 I11 c12 I 22 c13 I 33 2 c11 c12 I12 2 c11 c13 I13 c12 c13 I 23
(A.3)
Con esta expresión se obtiene el momento de inercia del cuerpo respecto a un eje
r i1c cualquiera en función de los momentos de inercia respecto a las direcciones
coordenadas originales y de los productos de inercia (o momentos centrífugos)
respecto a los pares de planos coordenados.
Con el fin de establecer de una manera gráfica la ley de variación de I’11, se toma sobre
1
la recta un punto R situado a una distancia rR de 0. Se tiene que
I '11
2 1 x x1
r y teniendo en cuenta que: cos 11 = c11 = 1
c
I11 rR rR
y también que
x2 x3
c12 c13
rR rR
entonces:
243
Mecánica Racional
2 2 2
1 cc x1 x2 I 22
I11 cc I 33
cc x3 (A.5)
x2 y2 z2
2 2 1
a2 b c
244
Apéndice 1
1 1 1
a1 ; a2 ; a3
I11cc 1/ 2
I 22
cc
1/ 2
I 33
cc
1/ 2
Siendo nulos los momentos centrífugos para estas direcciones de los ejes, el tensor de
inercia toma la forma
ª I11
cc 0 0 º
« »
I o « cc
I 22 0 »
«¬ Sim cc »¼
I 33
En cuanto a la forma del elipsoide, si el cuerpo es por ejemplo una esfera o un cubo, el
elipsoide será una esfera. A medida que el cuerpo se alarga o se acorta, también lo
hace el elipsoide, acompañando la forma del cuerpo.
A modo de ejemplo, analícense las expresiones de los momentos de inercia para el
caso de un bloque rectangular de lados a, b, h, tomadas del apartado siguiente (A1.4):
I xx
m
12
b2 h2
I yy
m
12
a 2 h2
I zz
m
12
a 2 b2
I xy I xz I yz 0
245
Mecánica Racional
En la figura siguiente se nota que si aumenta h, I33 (Izz) permanece constante, pero I11
(Ixx) e I22 (Iyy) aumentan, por lo que disminuyen los semiejes en x e y del elipsoide.
Gráficamente:
Asimismo, al alejarse del baricentro el semieje disminuye (por Steiner), pues aumenta
el momento de inercia.
Esfera
2
I xx I yy I zz mR 2
5
I xy I xz I yz 0
246
Apéndice 1
______________________________________________________________________
Hemisferio
I xx I yy 0.259 mR 2
2
I zz mR 2
5
I xy I xz I yz 0
______________________________________________________________________
Cilindro
1
I xx I yy m ( 3R 2 h 2 )
12
1
I zz mR 2
2
I xy I xz I yz 0
_____________________________________________________________________
Bloque rectangular
I xx
m
12
b2 h2
I yy
m
12
a 2 h2
I zz
m
12
a 2 b2
I xy I xz I yz 0
247
Mecánica Racional
______________________________________________________________________
m 2
I xx b
12
m 2
I yy a
12
I zz
m 2
12
a b2
I xy I xz I yz 0
______________________________________________________________________
Barra delgada
1
I xx I yy ml 2
12
1 2
I x cx c I y cy c ml
3
I zz 0
I xy I xz I yz 0
______________________________________________________________________
1
I xx I yy mR 2
4
1
I zz mR 2
2
I xy I xz I yz 0
248
Apéndice 1
______________________________________________________________________
Anillo delgado
1
I xx I yy mR 2
2
I zz mR 2
I xy I xz I yz 0
______________________________________________________________________
Cono
I xx I yy
3
80
m 4R2 h2
3
I zz mR 2
10
I xy I xz I yz 0
______________________________________________________________________
Semicilindro
m 2
I xx 0.0699 mR 2 h
12
I yy 0.320 mR 2
I zz
1
12
m 3R 2 h 2
I xy I xz I yz 0
249
Apéndice 2
APÉNDICE 2
251
Mecánica Racional
Calcular
Resolución
G n
G G G
Z ¦Z
i 1
i Z1 Z 3 Z1 Kˆ Z 3 ( senTJˆ cos TKˆ )
252
Apéndice 2
G
Para T 30q , se tiene: Z 15 Jˆ 34 Kˆ rad/s
G
Tomando al punto B para el cálculo del invariante escalar (2.21), será: P VB Zˆ
G G
Es, VB Z1 0,2( Iˆ Jˆ ) 1,6( Iˆ Jˆ ) m/s
y Z 15 2 34 2 37,16 rad/s
G
Z
Zˆ 0,4 Jˆ 0,91Kˆ
Z
Con este valor del invariante escalar, surge que el movimiento es del tipo helicoidal
instantáneo, lo cual se corresponde con las dos rotaciones alabeadas que afectan al
disco.
2) Velocidad del punto A en la periferia del disco aplicando el método del movimiento
absoluto, tomando como centro de reducción al punto O y al punto B.
G G G G
VA VB Z rBA
Del apartado anterior, es:
G G
VB 1,6( Iˆ Jˆ ) m/s y Z 15 Jˆ 34 Kˆ rad/s
G
rBA 0,2( sen30 Jˆ cos 30 Kˆ ) 0,1( cos 30 Jˆ sen30 Kˆ ) 0,0134 Jˆ 0,22 Kˆ m
253
Mecánica Racional
G G
Así: Z rBA 2,86 Iˆ
G
y, resulta finalmente: VA 1,26 Iˆ 1,6 Jˆ m/s
G G G G
V0 Z 3 rB 0 (15 J 25,98 Kˆ ) 0,2( Iˆ Jˆ ) 5,196 Iˆ 5,196 Jˆ 3Kˆ
G G
Z r0 A (15 Jˆ 34 Kˆ ) (0,2 Iˆ 0,213Jˆ 0,22 Kˆ ) 4,24 Iˆ 6,8 Jˆ 3Kˆ
G
Y por lo tanto, VA 1,26 Iˆ 1,6 Jˆ m/s
G
Por lo que, finalmente: aA 12,8 Iˆ 110,61Jˆ 44,3Kˆ m/s2
5) Velocidad del punto A, aplicando el método del movimiento relativo.
254
Apéndice 2
G G G
VrelA Z 3 rBA (15 Jˆ 25,98Kˆ ) (0,0134 Jˆ 0,22 Kˆ ) 2,96 Iˆ m/s
G G G G
VarrA V B Z1 rBA 1,6( Iˆ Jˆ ) 8 Kˆ (0,0134 Jˆ 0,22 Kˆ ) 1,704 Iˆ 1,6 Jˆ m/s
G
Resultando, V A 1,26 Iˆ 1,6 Jˆ m/s
Por lo tanto:
G G G G G G G
a arrA a B J 1 rBA Z1 Z1 rBA 12,8( Iˆ Jˆ ) 8 Kˆ 0,107( Iˆ) 12,8 Iˆ 13,65 Jˆ
G G G G
a relA Z 3 Z 3 rBA (15 Jˆ 25,98Kˆ ) 2,96 Iˆ 76,9 Jˆ 44,4 Kˆ m/s2
G G G
a comA 2Z 1 VrelA 2 8 Kˆ 2,96 Iˆ 47,36 Jˆ m/s2
G
Resultando: aA 12,8 Iˆ 110,61Jˆ 44,3Kˆ m/s2
255
Mecánica Racional
1 1
e e1 e 2 mV B2 I ZZ B Z 2
2 2
G G G
Ya que es e3 =0 por cuanto rG rB 0
G
Siendo VB 1,6( Iˆ Jˆ ) m/s, resulta VB2 = 5,12 (m/s)2 y por lo tanto
e1 = 0,5 x 2 x 5,12 = 5,12 J
Según los datos, los momentos de inercia y centrífugos del conjunto rotante con
respecto a los ejes x, y, z fijos a él, resultan:
256
Apéndice 2
ª0,02 0 0 º
IB « 0 0,01 0 »»
«
«¬ 0 0 0,02»¼
G
En esta base, el vector Z toma la forma:
G
Z 8(cos Tˆj senTkˆ) 30 ˆj 36,93 ˆj 4kˆ
Por lo que el ángulo entre el vector rotación y el eje coordenado y, resulta:
4
tgE 0,108 `@¹
36,93
Para el cálculo del momento de inercia I ZZ B del conjunto, se asignará a la dirección del
G
vector Z la notación y”, por lo que será I ZZ B I 2"2" . Utilizando la notación indicial
(A.1), se tendrá:
257
Mecánica Racional
Donde:
c 2"1 cos 90q 0
c 2" 2 cos 6,18q 0,994
c 2"3 cos 276,18q 0,108
Y por lo tanto:
I 2" 2" 0,994 2 .0,01 0,108 2 .0,02 0,0101 Kg m2
Luego,
e2 0,5.0,0101.(36,96 2 4 2 ) 6,97 J
Resultando, finalmente:
e e1 e 2 12,1 J
9) Momento cinético del conjunto respecto del punto B tomando como sistema
coordenado el {B, x, y z}.
G G
Según la expresión (4.16.b), es K G I G Z .
Por lo tanto:
ª0,02 0 0 ºª 0 º ª 0 º
G G «
KB I B Z « 0 0,01 0 »» ««36,93»» «0,37 »
« »
«¬ 0 0 0,02»¼ «¬ 4 »¼ «¬0,08»¼
o:
G
K B 0,37 ˆj 0,08kˆ Kg m2/s
G
10) Verificación del valor obtenido para K B en el inciso anterior tomando como
$
! ³µ } ¶· &
\1, paralelo al que se observa en la
figura con origen en 0.
258
Apéndice 2
Denotando X = 1’, Y = 2’ y Z = 3’, los ángulos entre los ejes coordenados X, Y, Z y los x,
y, z, son:
Por lo que los momentos de segundo orden del conjunto respecto de los ejes X,Y,Z,
resultan:
I 1'1' I 11 0,02 Kg m2
Por lo tanto:
ª0,02 0 0 ºª 0 º ª 0 º
G G «
KB I B Z « 0 0,0175 0,00433»» ««15 »» « 0,12 »
« »
«¬ 0 0,00433 0,0125 »¼ «¬34»¼ «¬0,36»¼
259
Mecánica Racional
G
o: K B 0,12 Jˆ 0,36 Kˆ Kg m2/s
Verificación:
G
reemplazando en la expresión para K B :
G
K B 0,12(0,5 ˆj 0,866kˆ) 0,36(0,866 ˆj 0,5kˆ) 0,37 ˆj 0,08kˆ Kg m2/s
12) Cálculo de las reacciones estáticas en el punto 0’de unión entre el bastidor y el eje
Z.
Utilizando las ecuaciones de la Estática:
G G
¦ F FE 0
G G
¦ M FE B 0
260
Apéndice 2
Rx = 0
Ry = 0
Rz – mg = 0 , de donde Rz = mg
De las segundas:
X) mg .0, 2 M 0' x 0
Y) mg .0, 2 M 0' y 0
Z) M 0' z 0
G G
Z1 Q 25,6( Iˆ Jˆ ) N
y en consecuencia:
RX = -25,6 N
RY = -25,6 N
RZ = 0
Estas reacciones en el plano del bastidor representan el esfuerzo que debe ejercer el
G
vínculo para soportar la fuerza centrífuga del conjunto generada por la rotación Z1 .
Teniendo en cuenta ahora que las reacciones presentes en 0’ -vínculo de soporte entre
el bastidor y el eje Z- ejecutan momento respecto al punto B a través de la estructura,
se tiene:
G
rB 0 ' R x Iˆ 0,2 R x Kˆ 0,1R x Jˆ
G
rB 0 ' R y Jˆ 0,2 R y Kˆ 0,1R y Iˆ
G G
rB 0 ' R z Kˆ 0
G
donde: rB 0 ' 0,2 Iˆ 0,2 Jˆ 0,1Kˆ m, y por lo tanto, aplicando la ecuación de Euler en
componentes:
MBx – 0,1 Ry = -0,96 N m
MBy + 0,1 Rx = 0
MBz - 0,2 Rx + 0,2 Ry = 0
de donde resulta:
MBx = 1,6 N m ; MBy = -2,56 N m ; MBz = 0 N m
263
BIBLIOGRAFÍA
LIBROS
LIBROS DE LECTURA
ENLACES DE INTERÉS
APÉNDICE 1
APÉNDICE 2
El presente trabajo sobre Mecánica Racional cons tuye un aporte que es fruto de muchos años de
experiencia en la materia y que ha sido concebido como una solución para el alumno en lo que respecta
al seguimiento de las clases con un único texto.
La obra está dirigida a los cursantes de la asignatura mecánica Racional del tercer nivel de la carrera de
grado de Ingeniería Mecánica de la Universidad Tecnológica Nacional, aunque puede ser usada como
base por alumnos de otras ingenierías, quienes debieran adaptarla a los contenidos de sus programas.
Liberto Ercoli es Profesor en Docencia Superior, Ingeniero Mecánico y Doctor en Ingeniería. Es Decano y
Profesor Ordinario Titular de Mecánica Racional en el Departamento de Ingeniería Mecánica de la
Facultad Regional Bahía Blanca, Universidad Tecnológica Nacional, donde desde 1979 se desempeña
como docente inves gador. Ha realizado numerosos cursos de perfeccionamiento en el país y en el
exterior, especializándose en acús ca, vibraciones mecánicas y diagnós co de maquinarias.Fue
miembro ac vo del Vibra on Ins tute y de la American Academy of Mechanics, siéndolo actualmente de
la Acous cal Society of America. En 1987 se hizo acreedor a uno de los premios "Dr. B. Houssay"
otorgados por el CONICET.
Ha dirigido proyectos de inves gación, par cipado como jurado en numerosos concursos docentes y
como evaluador de tesis de posgrado y de proyectos en diferentes universidades, además de revistar
como referí de diversas publicaciones cien ficas.
Posee más de sesenta publicaciones en revistas internacionales con referato, como: Journal of Sound
and Vibra on, Nuclear Engineering and Design, Applied Acous cs, Ocean Engineering, Acus ca, Journal
of Mechanical Engineering Educa on, Journal of the Acous cal Society of America, J. of Building
Acous cs, etc. Sus trabajos han par cipado en numerosos congresos y reuniones cien ficas de la
especialidad.
Virginia Azurmendi es Ingeniera Civil. Profesora Titular en los Departamentos de Ciencias Básicas, de
Ingeniería Mecánica e Ing. Eléctrica de la Facultad Regional Bahía Blanca, Universidad Tecnológica
Nacional y es Profesora Asociada del Departamento de Ingeniería de la Universidad Nacional del Sur.
Ha sido jurado de gran can dad de concursos docentes, es miembro de comisiones evaluadoras de
carreras académicas y es docente inves gadora de la UTN.
Tiene cargos direc vos de ges ón universitaria en la Universidad Nacional del Sur y en la Facultad
Regional Bahía Blanca de la Universidad Tecnológica Nacional siendo en la actualidad Directora del
Departamentos Ciencias Básicas.
ISBN 978-987-1869-23-3
9 789871 896233