08 Dinámica de Lagrange
08 Dinámica de Lagrange
08 Dinámica de Lagrange
Cinemtica Sistemas de Partculas Ecuaciones Cardinales de la Dinmica Dinmica del Movimiento Plano Campo de Fuerzas Oscilador Lineal Unidimensional Funcin Lagrangiana
Este cuadernillo:
08 Dinmica de Lagrange
2011
Diego E. Garca
142
Dinmica de Lagrange
Presentacin
El presente texto Teoremas de la Dinmica, corresponde a una seleccin de contenidos tomados del amplio campo de la Mecnica Clsica. Est destinado a los estudiantes de diversas ramas de Ingeniera, que hayan completado los cursos de Anlisis Matemtico y de Fsica. Parte de ellos, corresponden a las clases Mecnica Analtica (Ingeniera Civil), en el Departamento de Fsica de la F. de C.E.F y N. de la UNC, a mi cargo entre 2001 y 2010. El trabajo lo desarroll a partir de una Mecnica de nivel intermedio, necesariamente limitado a un curso breve, por las restricciones de tiempo en la formacin de grado en Ingeniera. En este cuadernillo se analizan la configuracin de un sistema de partculas, los conceptos de desplazamientos reales y virtuales, coordenadas generalizadas y se presentan, sin demostracin, las Ecuaciones de Lagrange, expresadas a partir de la Funcin Lagrangiana, todo ello a efectos de servir de introduccin para la resolucin de algunos problemas por el mtodo de Lagrange. En cuadernillo separado se desarrolla este tema con mayor amplitud.
8 Dinmica de Lagrange
Contenidos de este cuadernillo Determinacin de la posicin de un sistema de partculas, parmetros de configuracin con vnculos fijos, (sistema esclernomo) Configuracin de un sistema de partculas en el caso de vnculos mviles, (sistema renomo) Desplazamiento virtual y desplazamiento real de las partculas de un sistema Grados de libertad, sistemas holnomos y anholnomos Coordenadas generalizadas Algunos ejemplos de coordenadas generalizadas Funcin Lagrangiana, ecuaciones de Lagrange Resolucin del pndulo simple mediante las ecuaciones de Lagrange Algunas aplicaciones de la funcin lagrangiana 1 Pndulo simple con punto de suspensin oscilante 2 - Cadena deslizante 3 Cilindro rodante en plano inclinado 4 - Partcula en un plano vertical rotante Referencias Bibliogrficas
Pgina
143 144 145 146 147 147 149 150 152 152 159 161 163 168
Derechos reservados. Ley 11723 ISBN 978-987-05-4041-0 Impreso en la ciudad de Crdoba, Argentina 2 edicin, marzo 2008 Reimpresiones en 2009, 2010, 2011
Teoremas de la dinmica.
143
Dinmica de Lagrange
8 Dinmica de Lagrange
Determinacin de la posicin de un sistema de partculas, parmetros de configuracin con vnculos fijos, (sistema esclernomo)
Un sistema con vnculos fijos, se denomina esclernomo. (del griego sklers, duro rgido y de nomos, ley). Imaginemos una partcula que se puede mover libremente en el espacio, en ese caso, su posicin queda determinada cuando se conocen, en cada instante, sus 3 coordenadas. Si se tienen N partculas, habr que fijar 3N coordenadas para determinar la posicin del sistema. El trmino configuracin del sistema se refiere a las posiciones que ocupan las partculas en un determinado instante: cuando se conocen las ubicaciones de las mismas, decimos que el sistema est configurado. Ahora bien: puede ocurrir que algunos puntos del sistema estn vinculados entre s, por ejemplo: 2 masas puntuales libres en el plano, suman 4 grados de libertad y hacen falta 4 coordenadas para configurar el sistema. Pero supongamos que las dos partculas estn vinculadas entre s por un segmento inextensible de longitud l : en este caso bastar fijar las 2 coordenadas de una de ellas ms el ngulo que forma el segmento con la horizontal, en total 3 coordenadas independientes entre s, para tener determinado el sistema. Ciertamente, esta restriccin de vnculo, se podr establecer mediante una determinada ecuacin que debern satisfacer las coordenadas. Esa ecuacin es:
x2
x1
y2
y1
l2
Se llaman parmetros de configuracin , o coordenadas de Lagrange, a las coordenadas que se usan, o que se pueden usar, para configurar un sistema. Decimos, que se pueden usar, porque, en realidad, la posicin de las dos masas acopladas por una barra, queda definida tambin usando los 4 parmetros x1 x2 x3 x4 , aunque si la definimos as, estaramos usando un parmetro superabundante: Llamaremos por lo tanto parmetros de Lagrange a todas las coordenadas de configuracin, ya sean estas superabundantes o ya sean independientes. Designaremos con qk a estos parmetros, con
k 1, 2,3,...n .
De acuerdo con el ejemplo que se ha mostrado, resulta entonces, que, si los parmetros de configuracin deben satisfacer una ecuacin de vnculo, eso hace que el sistema se configure con 3, en vez de con 4 parmetros. Hablando en general, si llamamos n a la cantidad de parmetros de Lagrange y m a la cantidad de ecuaciones que representan las restricciones de vnculo, podemos expresar:
Cantidad de parametros independientes n m Si llamamos P O al vector posicin de una partcula i del sistema, podemos expresar i
finalmente la configuracin de un sistema de partculas mediante la siguiente funcin vectorial:
Pi O
tiempo.
Pi qk
con
k 1, 2,3,... n m
(1)
en donde los qk son parmetros independientes entre s. Cada uno de estos parmetros es funcin del Tambin podemos expresar, en forma escalar, las coordenadas de posicin xi , yi , zi de una partcula i del sistema, en funcin de los n
m parmetros independientes:
xi yi xi qk yi qk
(1) con
k 1, 2,3,... n m
zi
zi qk
Si N es la cantidad de partculas del sistema tendremos 3 N ecuaciones como las (1) que nos permitirn conocer las coordenadas de cada una de las N partculas en funcin de los parmetros qk independientes. Entonces podemos llamar a las (1) o (1), ecuaciones de transformacin, porque nos permiten pasar de los parmetros de configuracin, a las coordenadas de posicin de las partculas. Tambin podramos designarlas como ecuaciones de configuracin.
Diego E. Garca
144
Dinmica de Lagrange
Debemos observar que en las (1) y (1), no aparece el tiempo en forma explcita: ello significa que si bien las coordenadas xi , yi , zi dependen del tiempo, lo hacen a travs a travs de las qk .
Esta no dependencia en forma explcita del tiempo, se debe a que hemos considerado que los vnculos permanecen fijos en el tiempo, por ejemplo, si un punto se mueve sobre una curva en el espacio, esa curva no cambia su forma con el tiempo, entonces, si dejamos constante un valor del parmetro ( por ejemplo el recorrido sobre la curva), la posicin de la partcula no cambia porque la curva no se mueve.
Veamos el siguiente ejemplo, figura 2: Una partcula P se mueve en un plano sujeto a un vnculo tal que su distancia al punto O es constante e igual a l , figura 1
y l
f x
q1
q2
y y la ecuacin de
ligadura es
2 1
2 2
0.
q1 , en cuyo caso q2 es superabundante. Tambin podemos tomar como parmetro independiente al ngulo en cuyo caso q1 y q2 pasan a ser
Podemos tomar como parmetro independiente a superabundantes. Si consideramos los 3 parmetros q1 , q2, , n
ecuaciones de ligadura ( m 2 ), de manera que la cantidad de parmetros independientes siempre ser 3 2 1. Las ecuaciones de transformacin que dan las coordenadas de la partcula en funcin del parmetro son:
x l cos y lsen
donde es una cierta funcin del tiempo. Como vemos, el tiempo no aparece explcitamente en la expresiones que dan las coordenadas de P , cuando los vnculos son fijos.
o i
o1
O1 O
sen t i (2)
, ya que
j
y
Figura 2
f l i m
La posicin queda determinada con una nica coordenada la posicin de O1 est predeterminada por la (2). La coordenada vectorial de posicin de la partcula P ser:
P O
o bien:
sen t l cos
i lsen j (3)
(3)
Como vemos, se trata de un sistema renomo y en la (3) aparece explcitamente el tiempo, cosa que no ocurre cuando los vnculos no se mueven.
Teoremas de la dinmica.
145
Dinmica de Lagrange
Pi O
con la siguiente ecuacin de ligadura:
q1 i q2 j , con q1
x y q2
q1 sen t
2 q2 l 2
0 (4)
Pi q1 , t . En esta funcin, el
tiempo aparece en forma explcita, por tratarse de un sistema renomo. Entonces, podemos expresar el desplazamiento real de la partcula i , como el diferencial de la funcin vectorial P , que se obtiene como i la suma de los 2 diferenciales parciales con respecto a cada variable: Pi Pi Vector desplazamiento real dPi dq1 dt (5) de la partcula i , con una q1 t nica coordenada q1 Si se tratase de un sistema de partculas con h coordenadas independientes, el desplazamiento real de una partcula, de un sistema renomo se expresara como:
k h
dPi
k 1
Pi dqk qk
Pi dt t
(6)
Si ahora, en cambio, queremos expresar el desplazamiento virtual, de una partcula del pndulo de la figura 3, desaparece el segundo trmino de la (5), porque suponemos que la ligadura no se mueve. Entonces, se tendr Vector desplazamiento Pi Pi q1 (7) virtual de la partcula i , con q1 una nica coordenada q1 Si se tratase de un sistema de partculas con h coordenadas independientes, el desplazamiento virtual de una partcula, de un sistema renomo se expresara como:
k h
Pi
k 1
Pi qk qk
(8)
Observamos que si se trata de un sistema esclernomo, o de vnculos fijos, las (5) y (6), se transforman en las (7) y (8). En las (7) y (8) se han indicado tanto la variacin de q1 como el desplazamiento de la partcula, con la letra en lugar de usar la letra d de diferencial de una funcin. El uso de
Pi para el
Diego E. Garca
146
Dinmica de Lagrange
desplazamiento, es para indicar que se trata de un desplazamiento virtual, que, como se dijo, se realiza imaginando que los vnculos se inmovilizan mientras ocurre la variacin de las coordenadas qk . El uso de
qk es para poner de manifiesto que dicha coordenada experimenta una variacin, mientras el
t1 t t2 .
Llamamos isocrnica a esa variacin
expresada que, mientras dura la variacin de las qk , el vnculo queda transitoriamente inmovilizado en la
t1 t t2 etc. Podemos establecer el siguiente resumen con respecto a las caractersticas salientes del desplazamiento virtual de una partcula i del sistema:
posicin que ocupaba en un determinado instante t
Suponiendo el vnculo fijo, debe ser compatible con la condicin de vnculo, es decir moverse segn la tangente a la trayectoria que le impone dicha condicin. Si los vnculos son mviles, se consideran inmovilizados mientras dura la variacin de las qk . En el caso que los vnculos sean fijos, no se establece diferencias entre el vector desplazamiento real dPi y el vector desplazamiento virtual Pi de cada partcula. Tal sera el caso que el vnculo fuese realmente fijo en la figura 3. Como no lo es , para tener dPi el vector desplazamiento de arrastre de Pi debido al movimiento del vnculo. No necesariamente debe ser un valor diferencial: puede ser un valor finito, siempre cuando los desplazamientos de las partculas del sistema, sean compatibles con los vnculos. El smbolo usado para indicar un desplazamiento virtual, cuando ste es una magnitud diferencial y participa de las mismas propiedades usuales del smbolo d usado como diferencial. Un desplazamiento virtual puede ser un desplazamiento de un punto, como tal, pero tambin puede corresponder a la variacin de un ngulo. O sea que puede medirse en metros o en radianes.
Grados de libertad
Se llaman grados de libertad a cada uno de los desplazamientos virtuales, independientes entre s, que pueden tener las partculas de un sistema. Supongamos un sistema de N partculas materiales, tal que su configuracin est dada por n parmetros de Lagrange. q j con j 1, 2,3......n . El sistema est sometido a m ecuaciones de restriccin de vnculo, o tambin puede no tener vnculos. Esta definicin de grados de libertad es vlida independientemente de la forma que tengan las ecuaciones de vnculo, es decir, ya sean estas ecuaciones finitas o diferenciales, integrables o no. Entonces, si llamamos h a la cantidad de grados de libertad de un sistema, sta queda dada por: h n m (9)
m ecuaciones finitas, podemos obtener m parmetros del sistema en funcin de los otros n m
parmetros que hemos tomado como independientes. En ese caso, siempre ser posible considerar el sistema como si fuese un sistema libre, sin vnculos, pero que, en vez de tener 3N grados de libertad, tuviese h 3 N m grados de libertad.
En un sistema holnomo, la cantidad de parmetros o coordenadas necesarios para configurar el sistema, coincide con la cantidad de grados del libertad del mismo.
Teoremas de la dinmica.
147
Dinmica de Lagrange
Decimos que un sistema es no holnomo o anholnomo si una, o ms, de las m relaciones de vnculo, est expresada como una ecuacin diferencial no integrable. Si la ecuacin diferencial que expresa el vnculo fuese integrable, dicha expresin se transformara en finita y el sistema pasara a ser holnomo. Es importante hacer notar que, si el sistema es anholnomo, no es posible encontrar los m parmetros de configuracin en funcin de los restantes. En este tipo de sistemas, la cantidad de grados de libertad sigue siendo h , pero es posible demostrar que ya no resultan suficientes h parmetros para configurar el sistema, sino que para ello son necesarios h de ecuaciones de vnculo que no son integrables. (ver nota 1)
Nota 1: puede consultarse al respecto, el captulo XIII del texto Mecnica Analtica del Ing. Nilo Penazzi, Ctedra de Mecnica, Ingeniera Aeronutica de la Facultad de Ingeniera de la UNC, 1998.
Coordenadas generalizadas
A los parmetros qk independientes entre s, que determinan la configuracin instantnea de un sistema y ya mencionados al hablar de la configuracin de un sistema, es usual y resulta cmodo, darles el nombre de coordenadas generalizadas. En los sistemas holnomos, como ya se dijo, la cantidad de coordenadas generalizadas coincide con la cantidad de grados de libertad h del sistema. En consecuencia podemos expresar las ecuaciones de transformacin que dan las coordenadas cartesianas de las N partculas, en funcin de las h coordenadas generalizadas qk , de la siguiente forma:
xi xi q1 , q2 ,... qh , t yi yi q1 , q2 ,... qh , t
i 1, 2,3,... N
y
k 1, 2,3,.........h
zi zi q1 , q2 ,... qh , t
En las ecuaciones (1) aparece el tiempo en forma explcita, porque se considera el caso general de vnculos mviles.
Algunos ejemplos de coordenadas generalizadas a) Partcula que se mueve sobre vnculo con forma de elipse
Pm
x a
Figura 1 La partcula tiene un solo grado de libertad. Tomamos como coordenada generalizada el ngulo , en forma tal que (figura 1):
y x
b) Cilindro que rueda sobre un plano inclinado
b sen a cos
Ecuaciones de transformacin
Si la rodadura se hace sin deslizamiento se tiene una sola coordenada generalizada, que puede ser la coordenada x ( ver en la figura 2). El ngulo se relaciona con x mediante la relacin:
Diego E. Garca
148
Dinmica de Lagrange
Figura 2 En cambio si hay deslizamiento en el punto de contacto, el sistema tiene 2 grados de libertad y las coordenada generalizadas son x y .
O1
a x
m
Y
Figura 3
Se tiene una barra puede girar en un plano vertical alrededor de la articulacin O1 , como se muestra en la figura 3. La masa m , puede deslizar sobre la barra y a su vez, est sujeta por un resorte a la articulacin O1 . El sistema est formado por la barra y la masa m y al resorte, lo consideramos sin masa. La ligadura del sistema con el exterior es la articulacin O1 , la que, a su vez, puede, o no, deslizarse a lo largo del eje X . a es una determinada longitud del resorte, a partir de la cual medimos las elongaciones x de la masa m. En la tabla siguiente, se analizan diversas alternativas que pueden presentarse
Situacin
Ecuaciones de Transformacin
Vnculo
XP
X P (x , )
O1 Fijo,
variacin de ley
conocida: Distancia O
t
siguen
Y P Y P (x , ) X P X P (x , t )
Y P Y P (x , t ) X P X P (x , t )
Esclernomo
Renomo
O1 y
O1
x
t
Coordenadas generalizadas
YP
Y P (x , t )
Renomo
Ecuaciones de Transformacin
Vnculo Renomo
O1 se puede mover sobre colisas deslizante y sigue ley conocida: Slo la distancia O
t O1 sigue
O1 , x
,
XP YP
XP YP
X P (x , O Y P (x , O
Y P (x , , t ) X P (x , , t )
O1 , t ) O1 , t )
Renomo
Teoremas de la dinmica.
149
Dinmica de Lagrange
x, O
O1
XP YP
X P (x , , O Y P (x , , O
O1 ) O1 )
Esclernomo
d) Pndulo doble
En la figura 4 se muestra un sistema formado por dos varillas articuladas que se conoce como pndulo doble. Las coordenadas generalizadas son 1 y 2 . Como se observa en la figura, 1 puede variar independientemente de
2
y adems
1.
!
1
2
Figura 4
x1
Nudo
x2 x2
Polea Figura 5
T
i 1
1 m i Pi 2
(10)
donde N es la cantidad de partculas del sistema. A su vez, el vector posicin Pi , es una funcin de las coordenadas generalizadas qk y del tiempo. Entonces, de acuerdo con la (10), la energa cintica de un sistema de partculas, puede expresarse como una funcin de las coordenadas generalizadas qk del sistema, de sus derivadas con respecto al tiempo, qk y del tiempo:
f qk , qk , t
(11)
Diego E. Garca
150
Dinmica de Lagrange
Supongamos ahora, que la totalidad de las fuerzas que obran sobre el sistema son de carcter conservativo: En este caso y de acuerdo con la definicin de fuerzas conservativas, es posible definir una cierta funcin energa potencial, correspondiente a dichas fuerzas. Sabemos, por otra parte, que en un campo de fuerzas, la energa potencial depende de la posicin de las partculas del sistema y del tiempo, pero no depende de la velocidad de las mismas. Ello significa que la energa potencial, ser una funcin de las coordenadas generalizadas y del tiempo, pero no de las derivadas con respecto al tiempo de dichas coordenadas generalizadas. Entonces, si V es la energa potencial del sistema, podemos expresar la misma como funcin de las coordenadas qk y del tiempo, de la siguiente forma: V f qk , t (12) De acuerdo con las consideraciones precedentes, en aquellos sistemas conservativos que admiten energa potencial V , vamos a definir una cierta funcin que llamaremos Funcin Lagrangiana, de la siguiente manera: T V (13) Funcin Lagrangiana Las ecuaciones que se expresarn a continuacin, se conocen como Ecuaciones de lagrange y fueron desarrolladas por Louis de Lagrange, en Mcanique Analytique, Paris, 1788. En ellas interviene la Funcin Lagrangiana y su importancia consiste en que permiten la resolucin de problemas, sin intervencin de las ecuaciones cardinales de la dinmica. En ese sentido, constituyen un mtodo autnomo de la dinmica. Las ecuaciones de Lagrange, que presentamos sin demostracin, se expresan de la siguiente manera:
d dt
de Lagrange:
qk
qk
(14)
1, 2,..., h
Podemos hacer el siguiente resumen de los pasos a seguir para plantear las ecuaciones
Determinar el nmero de grados de libertad del sistema y elegir las coordenadas generalizadas. Expresar la energa cintica T del sistema qk y qk y del tiempo (si se tienen vnculos mviles). Calcular la energa potencial V en funcin de
V .
Calcular las correspondientes derivadas del Lagrangiano con respecto a introducidas en las ecuaciones de Lagrange:
d dt
Es decir:
qk d dt
d dt d dt
qk
1, 2,..., h
q1
q2
q1
q2
0
0
(14)
........................
qh
qh
Teoremas de la dinmica.
151
Dinmica de Lagrange
l
cota cero
m m
Figura 1 La variacin de energa potencial cuando el peso mg experimenta un desplazamiento dxi , es igual, por definicin, al trabajo de producto escalar:
dV
mg dxi
dV
mgdx
x , resulta:
V
Pero:
mgx
x l 1 cos
(1)
l cos
mgl 1 cos
1 mv 2 2
donde v es la velocidad de la partcula. El mdulo de la misma lo calculamos multiplicando el largo del pndulo por su velocidad angular que vale:
d . Entonces: dt v l 1 m l 2
T
2
(2)
V
(3)
1 2 ml 2
d dt
O bien:
mgl 1 cos
De acuerdo con las (14) y como se tiene un solo grado de libertad, la nica ecuacin de Lagrange que se tiene en este caso es:
q1
q1
0
d dt
(4)
Seguidamente, efectuaremos las correspondientes derivadas del Lagrangiano, para poder reemplazar en la (4). Observamos que la (3) es funcin de
q1
y de q 1
explcita, porque se trata de un sistema con ligadura fija. Para derivar la (3) con respecto a debemos considerar a como constante y cuando derivemos con respecto a constante. Derivando (3):
q1
, supondremos
Diego E. Garca
152
Dinmica de Lagrange
ml 2
(5)
Adems, pare obtener el primer trmino de la (4), falta an derivar la (5) con respecto al tiempo, para lo cual deberemos tener presente que la velocidad angular es funcin del tiempo:
d dt
d ml 2 dt
d dt
ml 2
k
(6)
ml 2
Para oscilaciones pequeas, resulta que: Entonces, la (7) queda:
mgl sen
sen
(7)
l
O bien:
g l
A cos t
en donde:
2
0 (8)
Bsen t
La (8) es la conocida ecuacin diferencial de un movimiento armnico simple, cuya solucin es:
g ; l
2 T
g l
2 g l
l g
Ejemplos de aplicacin de la funcin Lagrangiana 1 Pndulo simple con punto de suspensin oscilante
Hallar, por el mtodo de Lagrange, la ecuacin diferencial del movimiento de un pndulo cuyo punto de suspensin oscila segn la ley x1 a sen t , suponiendo que la amplitud de las oscilaciones es pequea.
Resolucin: Grados de libertad El sistema tendra dos grados de libertad, si la posicin de la articulacin O1 no estuviese determinada. Pero como la posicin del punto o1 est definida por la funcin x1 (t ) que es un dato, tenemos slo un grado de libertad. Tomaremos el ngulo Energa cintica Las coordenadas de posicin de la masa son: (ver en figura 1) como coordenada generalizada.
Teoremas de la dinmica.
153
Dinmica de Lagrange
x1
o1
x
lcos
v
x1 l sen
Figura 1
x y
x1 l cos
l sen
O1
yj
v
Derivamos con respecto al tiempo para obtener las componentes del vector velocidad: (ver en figura 2)
x y
Entonces, podemos expresar:
x1
l cos
l sen
Figura 2
xi
v2
x2
y2
1 2
mv 2 , resulta:
T T T
1 2 1 2 1 2
x1
l cos l2 l2
2
l sen 2x 1l cos
m x 12 m x 12
cos2
l2
sen 2
2x 1l cos
y (aparte
De la expresin precedente surge que la energa cintica es funcin de x1 , x1 , de las constantes m y l Energa Potencial )
Tomamos el potencial cero de referencia en el eje x, y teniendo en cuenta la definicin de potencial resulta: (ver en figura 3) dV = -m g dy V=-mgy
y y
mg
Figura 3 Entonces:
mgl cos
Diego E. Garca
154
Dinmica de Lagrange
L
L L
m x1 m x1
l2 l2
a continuacin efectuaremos las derivadas que corresponden para poder escribir el binomio de Lagrange:
1 2
m 2l 2
m l2
d L dt
Tngase en cuenta que:
derivamos ahora con respecto al tiempo para tener el 1 trmino del binomio:
m l2
x1l cos
x1l sen
(1)
d cos dt
sen
Para hallar el segundo trmino del binomio de Lagrange, para lo cual derivamos el lagrangiano con respecto a :
mlx1 sen
mgl sen
(2)
d L dt
resulta:
m l2
x1l cos l2
0 0 0
x1l cos
2
x1 cos l
si x1
g sen l
(3)
a sen t , ser:
x1 x1 x1
reemplazando (4) en (3), queda:
a sen t a cos t a
2
(4)
sen t
d2 dt 2
a l
sen t cos
g sen l
Teoremas de la dinmica.
155
Dinmica de Lagrange
Esta ecuacin diferencial permitira encontrar la funcin solucin de la misma caso ms general en que la amplitud de las oscilaciones no sea pequea.
f (t ) , en el
Se trata de una ecuacin diferencial ordinaria de segundo orden que no es lineal. Teniendo en cuenta las condiciones impuestas en el enunciado de pequeos valores de , podemos considerar que:
cos sen
la cual permitir encontrar
1 0
d2 dt 2
Los valores de g; a; l; y diferencial constar de dos partes:
g l
a l
sen t (5)
Una parte, ser la solucin de la homognea asociada La otra parte, ser la solucin particular. La ecuacin homognea asociada es:
d 2 dt 2
g l
0 (6)
. La solucin de
esta ecuacin es la misma del oscilador armnico simple (ver la (4) de ese tema), donde en vez de la elongacin
x , aparece el ngulo
. Entonces es:
hom
C 1cos
g t l
C 2sen
g t (7) l
donde
, C2
0 0
A cos C2 . C1
Mcos t
(7)
donde A
C 12
2 C2
arc tan
N sen t
(8)
Si reemplazamos la funcin
part
a l
2 0
0,
2 2
Diego E. Garca
156
Dinmica de Lagrange
A cos
Nsen t
(9)
0
, lo
que ya no es una funcin armnica simple, el valor de un determinada elongacin mxima, no es igual a la que le sigue. En forma cualitativa, la curva
Figura 4 La funcin
f t
0,
es anloga a la funcin x t
Si fuese
frecuencia natural del pndulo, la funcin (8) ya no es solucin de la ecuacin diferencial (5). Debe buscarse una solucin de la forma (ver ):
part
M cos t
Nsen t t
(10)
Puede consultarse el texto clsico de Anlisis Matemtico N. Piskunov, Clculo diferencial e Integral, captulo Ecuaciones Diferenciales, artculo ecuaciones diferenciales de las oscilaciones mecnicas.
a 2l
0 o bien, siendo
a 0 2l
(11)
a 0 t cos 0t 2l
(12)
La elongacin
A cos
Como vemos, cuando
0
a 0 t cos 0t 2l
(13)
f t
(segundo
trmino de la (13), es anloga a la funcin x part del oscilador forzado sin amortiguamiento, ver el tema oscilador forzado sin amortiguamiento caso a), segundo trmino de la expresin (9). Ver nota 2
Teoremas de la dinmica.
157
Dinmica de Lagrange
La grfica de la funcin representada por la solucin particular, en el caso de resonancia, se muestra en la figura 5 siguiente y es anloga a la del oscilador forzado sin amortiguamiento en resonancia, tema oscilador forzado sin amortiguamiento caso a), figura 4 de ese tema.
part
Figura 5 No obstante, debe recordarse que este problema ha sido resuelto con la aproximacin de 1 y sen 0 y que tanto la ecuacin diferencial (5), como sus soluciones (9) y (13), son vlidas para pequeas amplitudes de oscilacin. Entonces la figura (5) ser representativa mientras que las amplitudes de la oscilacin se ajusten a esta hiptesis. Cuando las amplitudes sean mayores, la ecuacin diferencial que describe el comportamiento del sistema, ser de la forma
cos
d2 dt 2
Nota 2
g sen l
a l
cos sen t .
Podemos plantear cierta analoga entre la ecuacin diferencial (8) de oscilador forzado sin amortiguamiento, captulo oscilador lineal, y la ecuacin diferencial (5), en donde Adems, de este problema sera equivalente al frecuencia natural del resorte.
a l
sen t
sera anlogo al
0
Q0 cos t m
de la (8).
Es interesante destacar que mediante la formulacin de Lagrange hemos hallado la ecuacin diferencial que permite resolver el sistema, por un camino prescindente de las ecuaciones cardinales de la dinmica; si hubiremos usado estas ecuaciones, habramos llegado a la misma ecuacin diferencial. Resolucin usando la ley de Newton en un sistema no inercial Elegimos un sistema de referencia x1 e y1 que tiene movimiento armnico simple de traslacin, dado por = a sen t, se trata de un sistema llamado no inercial por que todos los puntos del mismo estn sometidos a una aceleracin de arrastre que vale:
d2 dt 2
vale:
sen t
(14)
En estas condiciones, cualquier masa de ese sistema est sometida a fuerzas de inercia que
d2 dt 2
Las fuerzas que actan sobre la masa en este sistema no inercial se muestran en la figura 6 y son:
y1
O1
T
x1
m
mg
Figura 6
F in
Diego E. Garca
158
Dinmica de Lagrange
T cos90 mg sen
a
Fin cos
ma
x resulta que a
mg
Adems:
d2 dt 2
mx1
(15)
x1 x1
Reemplazamos en (15)
l sen l ; x1 l y x1 l (16)
d2 y x1 dt 2
g a 2 sen t l o bien dividiendo ambos miembros por l , llegamos finalmente a la misma ecuacin diferencia (5) que
habamos obtenido por el mtodo de lagrange:
g l
a l
sen t
Teoremas de la dinmica.
159
Dinmica de Lagrange
2 - Cadena deslizante Una cadena cuya densidad lineal es k y apoyada sobre un plano horizontal liso como se muestra en la figura, se suelta, desde una posicin inicial en reposo. En el instante inicial, la coordenada del extremo inferior de la cadena es y y 0 .
Se pide encontrar la coordenada y en funcin del tiempo Resolucin Se trata de un sistema con un solo grado de libertad y adoptamos la coordenada y como coordenada generalizada. Energa cintica La energa cintica del sistema se expresa como
1 2
mv 2
(1)
donde m es la masa del sistema, o sea, en este caso la masa de la cadena y la misma para todos los puntos de la cadena. Es:
v es la velocidad, que es
m l
kl
(2) (3)
v
Reemplazamos (2) y (3) en (1):
dy dt
v
1 2
y
(4)
k l y2
Energa potencial Tomamos como referencia el plano horizontal, entonces la energa potencial de la cadena coincide con la energa potencial del tramo y:
mg
y 2
k yg
y 2
1 2
k g y 2 (5)
L T V
Reemplazamos los valores de T y V de (4) y (5):
L
L y
1 2
k l y2
1 2
k g y 2 (6)
L y
kl y
d dt
L y
kl y
(7)
Ahora, buscamos el segundo trmino del binomio de Lagrange, para lo cual derivamos el lagrangiano con respecto a y
k gy
(8)
d L dt y
L y
Diego E. Garca
160
Dinmica de Lagrange
y
cuya solucin es:
g y l
(9)
y
en donde:
A cos h ( t )
B sen h ( t ) (10)
g l
y0
de donde surge que:
La segunda condicin inicial, es la de velocidad cero, y 0 Para encontrar la constante B , derivamos la (10):
y
Valuamos la (12) en t
(13)
y 0 cosh
g t l
Teoremas de la dinmica.
161
Dinmica de Lagrange
3 Cilindro rodante en plano inclinado Se tiene un cilindro de radio R y masa m que rueda sin resbalar sobre un plano inclinado de ngulo y altura h , como se muestra en la figura siguiente. El cilindro parte del reposo.
X
R
h
Se pide encontrar la ecuacin diferencial del movimiento, mediante el mtodo de lagrange Parmetros de configuracin, coordenada generalizada La coordenada de posicin x del centro de la esfera y el ngulo , que ha girado el cilindro, desde su posicin inicial hastala posicin x , son parmetros de configuracin. Como no hay resbalamiento no son independiente una del otro, sino que estn relacionadas con la siguiente condicin de vnculo (ver en la figura 1)
x x0 x
G
x0
Rcos
h- xsen
Figura 1 El sistema tiene un solo grado de libertad. Elegimos como coordenada generalizada, a la coordenada x . Funcin Lagrangiana Se trata de un sistema conservativo, porque las fuerzas actuantes son gravitatorias y porque no hay disipacin de energa, debido a que, por hiptesis, el cilindro no desliza sobre el plano. En consecuencia, podemos usar la funcin lagrangiana para plantear la ecuacin de Lagrange correspondiente a la coordenada x . Encontraremos, en primer trmino, la energa cintica del cilindro rodante, usamos el teorema de Knig, con polo en el centro G
1 2
m x2
1 2
IG
El 1 trmino es la energa cintica de traslacin y el 2 es la energa cintica de rotacin; I G es el momento de inercia del cilindro con respecto al eje que pasa por G. En lo que sigue suprimiremos por simplicidad el subndice G. La velocidad del centro de masas G es x
x . R
Si reemplazamos
Diego E. Garca
162
Dinmica de Lagrange
1 2 mx 2
1 x I 2 R
I x2 R2 2
(1)
Para determinar la energa potencial del cilindro rodante, tomamos como referencia el plano horizontal. La altura del punto de contacto sobre la base, vale: h x sen , ver en la figura 1. Para encontrar la altura del centro G , habr que sumarle Rcos :
hG
La energa potencial es:
h x sen
R cos
V
o bien:
mghG
R cos
(2)
mg h
x sen
L
I x2 R2 2
T V
(3)
mg h
x sen
R cos
Efectuamos las derivadas parciales correspondientes para formar las ecuaciones de Lagrange:
L x L x
derivando ahora con respecto al tiempo:
I x R2 (4)
mg sen
d L dt x
I x R2 d L dt x L x
0
(5)
0
x mg sen I m R2
I x R2
mg sen
IG
mR 2 ; entonces 2
mg sen mR 2 m 2R2
x 2 g sen 3
g sen 1 1 2
Teoremas de la dinmica.
163
Dinmica de Lagrange
4 - Partcula en un plano vertical rotante Una partcula de masa m se mantiene sobre una placa sin masa que rota con velocidad angular constante alrededor del eje Z y no hay rozamiento entre la placa y la partcula. x , y , z son los ejes del cuerpo que rotan con la placa y X ,Y , Z son los ejes fijos. Las condiciones iniciales son, para t= 0
x0 , z
z0 ; x 0
0, z
Se pide encontrar las ecuaciones paramtricas de la trayectoria, referidas al sistema del cuerpo, planteando: a) Las leyes de newton, en el sistema no inercial x , y , z b) La funcin lagrangiana y las ecuaciones de lagrange c) Encontrar la fuerza de reaccin N que el plano ejerce sobre la partcula
mg
Y
t
X
a) Planteo de las Leyes de Newton en el sistema no inercial x y z Un punto cualquiera del plano x , z , tiene una aceleracin centrpeta que vale:
ac
z m F in
xi
mg
Figura 1
En consecuencia, la masa m estar sometida en ese plano rotante, a una fuerza de inercia de arrastre de sentido contrario a la aceleracin centrpeta ac y que vale (ver en la figura 1):
Fin Fin
x x
m ac m
2
(1)
Esta es la nica fuerza actuante en la direccin del eje x, por lo tanto la ecuacin diferencial de Newton en esta direccin es:
Diego E. Garca
164
Dinmica de Lagrange
Fin
d 2x (2) dt 2 x
2
De acuerdo con (1) y (2), la ecuacin diferencial del movimiento, resulta entonces:
mx
x 0 (3)
Se trata de una ecuacin de 2 orden lineal, homognea y cuya solucin, para las condiciones iniciales establecidas, da una de las ecuaciones paramtricas de la trayectoria, correspondiente a la coordenada x : (ver en la figura 2)
x 0 cos h( t )
x0
Figura 2
La pendiente de esta curva es la velocidad en x de la partcula; inicialmente es cero y es mayor a medida que se aleja del eje de rotacin. Con respecto a la ecuacin paramtrica correspondiente a la coordenada z : La nica fuerza actuante en la direccin de z es mg , de manera que la ecuacin de Newton en esta direccin es:
mg
d 2z m 2 dt
La partcula, como no poda ser de otro modo al no haber rozamiento, tiene, en la direccin una aceleracin g . Entonces:
z,
z
b) Planteo por el Mtodo de Lagrange
z0
1 2 g t (4) 2
Coordenadas generalizadas En principio, el sistema tiene tres grados de libertad, que corresponderan a las coordenadas x , z y . t , que es dato del problema. Pero el valor de ya est determinado previamente por la funcin Entonces, tenemos 2 grados de libertad y adoptaremos como coordenadas generalizadas a las coordenadas cartesianas x , z de la partcula en el plano rotante. Energa cintica El vector velocidad de la partcula, referido o visto desde el sistema fijo , pero expresado en funcin de sus proyecciones sobre los ejes en rotacin es:
xi
xj
zk (5)
En la figura (3) se ha dibujado el plano rotante x , y y all se han indicado las componentes v x y
vx
xi
m
Figura 3
Teoremas de la dinmica.
165
Dinmica de Lagrange
v2
La energa cintica de la partcula es:
x2
( x )2
z 2 (6)
T
2
1 mv 2 2 1 2 mz 2 1 m ( x )2 (7) 2
1 mx 2 2
Energa potencial, funcin lagrangiana, ecuaciones Planteamos V = mgz (8) L=TV Reemplazamos en la funcin lagrangiana los valores de (7) y (8):
1 mx 2 2 d dt L x d dt L z
1 2 mz 2
L x m L z mx
2
1 m ( x )2 2
(9) (10) (11) (12)
mgz
mz mg
Las ecuaciones de Lagrange , referidas a cada una de las coordenadas generalizadas x , z son:
d L dt x d L dt z
L x L z
0 0
Reemplazando (9) (10) (11) y (12) en las ecuaciones de Lagrange resultan las siguientes ecuaciones diferenciales, para cada una de las coordenadas x , z :
mx mz
m mg
x 0
x z g
x 0
(13)
La solucin de la segunda de las (13), para las condiciones iniciales dadas, es inmediata:
z z vz
g z0 gt 1 2 gt 2
(14)
Diego E. Garca
166
Dinmica de Lagrange
x
Para t
A sen h ( t )
B cos h ( t ) (15)
Determinaremos seguidamente las constantes, de acuerdo con las condiciones iniciales. 0 es x x 0 , entonces:
x0
senh 0
B cos h 0
x0
x
Pata t
A cos h( t )
B sen h( t ) (16)
0 , es x
x
La velocidad en
x 0 cos h( t ) (17)
x de la partcula es:
x 0 sen h ( t ) (18)
ma rel
interesa la proyeccin, en la direccin del eje
F i ma arr
ma cor (1)
y , de la ecuacin (1):
Fi ma arr
y y
ma rel
y
ma cor (2)
Con respecto a las fuerzas que intervienen en la (2), podemos decir que:
ma rel
y
0 (3)
ma arr
y
0 (4)
in x
porque la nica fuerza de inercia de arrastre que acta en la partcula, es la fuerza de inercia centrfuga, dada por la (1) de la parte (a) de este problema, F
La nica fuerza de interaccin que acta sobre la partcula en la direccin normal al plano, es
N j (5)
ma coriol
2 x j (6)
(7) :
Nj
m2 x j
Si expresamos en forma escalar y reemplazamos el valor de x de la (7) por su igual, dado por la (18) del tramo (b), resulta la siguiente expresin para el mdulo de la reaccin N
2m
x0 sen h( t )
(8)
En la figura 4 se observan la fuerza de Coriolis y la reaccin normal que actan sobre la partcula.
Teoremas de la dinmica.
167
Dinmica de Lagrange
Plano rotante
F iner cor
Figura 4 La (8) nos permite expresar las siguientes observaciones: La reaccin normal es creciente a medida que la partcula se aleja del eje de rotacin El valor de la reaccin normal coincide, en valor absoluto, con el de la fuerza de Coriolis. El sentido de la fuerza de Coriolis es tal, que hace que la partcula presione, empuje, contra el plano rotante y a su vez, ste ejerce una reaccin N sobre la partcula. En la direccin normal al plano, la fuerza resultante sobre la partcula es nula y corresponde a la suma de los vectores de igual mdulo y sentido opuesto F iner cor y N . En cambio, sobre el plano y en su direccin normal slo acta una fuerza - N : esta fuerza produce un momento con respecto al eje de rotacin en sentido contrario al sentido de giro del plano rotante. Esto significa que si pretendemos que la velocidad angular se mantenga constante, tal como indica el enunciado del problema ser necesario aplicar un momento exterior al sistema plano rotante masa, para que el plano no pierda velocidad angular. Significa tambin que debe entregarse una potencia al sistema para que la velocidad rotacin se mantenga constante. Esta potencia suministrada se traducir, a su vez, en una energa cintica creciente de la masa.
Diego E. Garca
168
Dinmica de Lagrange
Referencias bibliogrficas
Bedford Fowler, Mecnica para Ingeniera, Dinmica, Addison Wesley Iberoamericana, 1996 y siguientes. Boresi Arthur, Schmidt Richard, Ingeniera Mecnica:Dinmica, Universidad de Wyoming, editora Thomson Learning, Mxico 2002 Argello Luis R, Mecnica, editorial Answer Just in Time, Mxico, 2000 Meriam L. - Kraige L.G., Mecnica para Ingenieros, Dinmica, editorial Revert, S.A. 1998 Marsicano, Fnix R. Mecnica, UBA, Facultad e Ingeniera, 1985 y siguientes Wells, Dare A.. Teora y Problemas de Dinmica de Lagrange Mxico, Mx : McGraw-Hill, 1972 Steidel R.F. Jr, Introduccin al estudio de las vibraciones mecnicas, editorial Cecsa, 1981. Spiegel, Murray R.. Teora y Problemas de Mecnica Terica, Mxico, McGraw-Hill,1979-1980 Synge-Griffith, Principios de Mecnica,.Mc. Graw Hill Book 1965 y siguientes. Symon Keith, Mecnica, Editorial Aguilar, 1979 y siguientes. Beer, F Jhonston, E Mecnica Vectorial para Ingenieros. Mxico, McGraw-Hill, 1998-2003 Longhini Pedro, Mecnica Racional, Palumbo Editores, Buenos Aires 1960 Marsicano Fnix, Mecnica, Facultad de Ingeniera de la UBA, 1980. Hertig Ricardo, Mecnica Terica, el Ateneo editorial, 1976. Targ S. Curso Breve de Mecnica terica, editorial Mir 1971 Yavorski Detlaff, Manual de Fsica, editorial Mir, Mosc, 1975 y siguientes
La digitalizacin y edicin de los originales se efectu durante 2004, 2005, 2006 y 2007. Reimpresiones en 2008, 2009 y 2010. Crdoba, Argentina.
This document was created with Win2PDF available at http://www.win2pdf.com. The unregistered version of Win2PDF is for evaluation or non-commercial use only. This page will not be added after purchasing Win2PDF.