Mecánica Clásica 1 Prof. Cayetano Di Bartolo Departamento de Fı́sica Universidad Simón Bolı́var Esta guı́a está basada en los manuscritos que elaboré para los cursos de Mecánica que dicté en la Universidad Simón Bolı́var. La guı́a todavı́a requiere de modificaciones y correcciones, y es mi esperanza que en algún momento se convierta en un libro. Si el lector desea hacerme alguna observación puede escribirme a la dirección [email protected] AGRADECIMIENTOS El libro se está realizando con la magnı́fica colaboración de mi esposa Jacqueline Geille, quién contribuye en todos los aspectos de su elaboración. También agradezco al Profesor Lorenzo Leal, de la Universidad Central de Venezuela, que muy amablemente me facilitó sus notas para el curso de Mecánica. Ultima actualización: Julio de 2004 2 Dinámica de un sistema de partı́culas En este capı́tulo trataremos con sistemas de partı́culas. Se definirán “cantidades extensivas” que permitirán describir el movimiento global del sistema. Se definirán los momentos lineal y angular del sistema y sus energı́as cinética y potencial. Se estudiarán los sistemas con masa variable y se tratará el problema de dos cuerpos. 2.1 Centro de masa y momentum. Consideremos un sistema formado por N partı́culas. Llamaremos mα a la masa de la partı́cula α-ésima y rα a su vector posición en algún referencial S. Supondremos que existen fuerzas externas e internas al sistema. Llamaremos Fαe a la fuerza externa total sobre la partı́cula α-ésima y Fα,β a la fuerza que sobre dicha partı́cula aplica la β-ésima. m1 mN O rα mα S (referencial) Para determinar la evolución del sistema, i.e. el movimiento de cada partı́cula, debemos resolver las ecuaciones de movimiento de todas ellas. Si el referencial S es inercial el conjunto de ecuaciones es N e Fα,β , α ∈ {1, ..., N } (2.1) mα r̈α = Fα + β=1 donde Fα,α = 0. Se trata de 3N ecuaciones diferenciales de segundo orden acopladas. En muchas ocasiones este problema es difı́cil de resolver si no imposible. Sin embargo todavı́a cabe la posibilidad de estudiar la traslación global del sistema y para ello es útil el concepto de centro de masa del sistema. El centro de masa del sistema se define como un punto cuyo vector posición es rcm ≡ 1 mα rα M α con M≡ α 25 mα (2.2) (2.3) 26 Mecánica Clásica 1 C. Di Bartolo (Julio de 2004) la masa total del sistema. Nótese que si todas las masas son iguales el centro de masa coincide con el centro geométrico del sistema. Una cantidad ı́ntimamente relacionada con el centro de masa es el momentum total del sistema. Llamaremos pα = mα ṙα al momentum (o momento lineal) de la partı́cula α-ésima. Se define el momentum total del sistema en el referencial S como pα = mα ṙα (2.4) P ≡ α α y es inmediato que P = M ṙcm . (2.5) Para obtener la ecuación de movimiento que satisface el centro de masa en el referencial S (supuesto inercial) comencemos sumando las ecuaciones de movimiento (2.1) mα r̈α = α Fαe + α Fα,β . β Como Fα,β = −Fβ,α ⇒ α Fα,β = 0 β se llega a la ecuación mα r̈α = F Ext , (2.6) α donde F Ext ≡ Fαe (2.7) α es la fuerza externa neta que actúa sobre el sistema. De derivar respecto al tiempo las ecuaciones (2.4) y (2.2) se obtiene finalmente que la ecuación (2.6) es equivalente a M r̈cm = Ṗ = F Ext . (2.8) Esta ecuación es la ecuación de movimiento para el centro de masa y nos dice que el centro de masa se mueve igual de como lo harı́a una partı́cula con masa igual a la total del sistema y bajo una fuerza igual a la fuerza externa neta que actúa sobre el sistema. Al resolver la ecuación se obtiene la trayectoria del centro de masa que proporciona una idea del desplazamiento global del sistema. De (2.8) se obtiene la ley de conservación del momento lineal (válida en referenciales inerciales): F Ext = 0 ⇔ P = cte ⇔ ṙcm = cte . (2.9) C. Di Bartolo 27 Dinámica de un sistema de partı́culas En particular el centro de masa de un sistema aislado se mueve con velocidad constante respecto a cualquier referencial inercial. 2.2 Sistemas con masa variable. Llamamos sistemas de masa variable a sistemas de partı́culas con una distribución de masa no discreta cuya masa total varı́a en el tiempo. En esta sección escribiremos la ecuación de movimiento que satisface el centro de masa de tales sistemas. Consideremos un sistema de partı́culas con una distribución continua de masa. Etiquetaremos los puntos del sistema con uno o varios parámetros reales designados colectivamente con la letra λ y llamaremos r(t, λ) a sus vectores posición. El elemento de masa en el punto r(t, λ) del sistema lo escribiremos como ρdλ r(t, λ) dm = ρ(λ)dλ donde dλ puede ser el elemento de longitud, área o volumen y ρ(λ) es la densidad de masa correspondiente. Esto significa que la masa total del sistema es ρ(λ)dλ (2.10) M= Λ donde Λ es la región de integración o conjunto de valores de λ. Para este sistema la ecuación (2.6) toma la forma ∂2 F = r(t, λ)ρ(λ)dλ 2 Λ ∂t (2.11) con F la fuerza neta sobre el sistema. El centro de masa y el momento lineal del sistema satisfacen ∂r d ρdλ . (2.12) M ṙcm = r(t, λ)ρdλ y P = dt Λ Λ ∂t Si Λ = Λ(t) y las velocidades son funciones continuas del tiempo se cumple que d ∂r F = P con P = ρdλ = M ṙcm . dt Λ ∂t (2.13) Pero si Λ es función del tiempo como en los sistemas de masa variable la ecuación anterior no es válida. ¿Cómo obtener una ecuación de movimiento para el centro de masa de tales sistemas? Veremos que algunos problemas de masa variable pueden pensarse como sistemas 28 Mecánica Clásica 1 C. Di Bartolo (Julio de 2004) que coliden o interaccionan. Comenzaremos estudiando cómo se relacionan los cambios de momentum de dos sistemas que chocan. Sea un sistema continuo de partı́culas SA ∪ SB compuesto de dos subsistemas que no varı́an su masa, supondremos que las velocidades de todas las partı́culas son funciones continuas del tiempo. Sean PA y PB sus momentos y F la fuerza neta sobre el sistema compuesto. Entonces (2.13) conduce a F = SA ∪ SB SA SB d (PA + PB ) . dt (2.14) Si integramos la ecuación anterior en el intervalo de tiempo [t, t + ∆t] y llamamos ∆P a los cambios en los momentos se obtiene que t+∆t F · dt = ∆PA + ∆PB = mA ∆vA + mB ∆vB (2.15) t donde hemos llamado mA , mB , vA y vB a las masas y velocidades de los centros de masa de los subsistemas SA y SB . Algunos problemas de masa variable pueden entenderse como el problema de subsistemas que acabamos de estudiar. Consideremos un sistema S(t) de masa variable M (t) y centro de masa con velocidad V (t). Supondremos que SB es la porción de este sistema que se agrega o se pierde en el intervalo [t, t + ∆t]; llamemos m = mB a la masa de esta porción y ∆v = ∆vB al cambio de velocidad de su centro de masa en el intervalo. Entonces el cambio de masa del sistema S(t) en el intervalo [t, t + ∆t] es |M (t + ∆t) − M (t)| = m ≡ mB . (2.16) Supondremos que cuando SB forma parte del sistema S la velocidad de su centro de masa coincide con la del resto del sistema,i.e., vB = vA = vA∪B . Si se trata de un sistema que gana masa podemos tomar que S(t) = SA M (t) = mA S(t + ∆t) = SA ∪ SB ∆V = vA∪B (t + ∆t) − vA (t) = ∆vA . Si se trata de un sistema que pierde masa tomamos S(t) = SA ∪ SB M (t) = mA + mB S(t + ∆t) = SA ∆V = vA (t + ∆t) − vA∪B (t) = ∆vA . Si ahora suponemos que el intervalo de tiempo ∆t es muy pequeño (infinitesimal) al igual que ∆V y ∆M mientras que ∆v es finito, entonces la ecuación (2.15) conduce para cualquiera C. Di Bartolo 29 Dinámica de un sistema de partı́culas de los dos casos considerados a la ecuación dM dV ∆v + F =M dt dt (2.17) donde M , V y F son respectivamente la masa del sistema, su velocidad y la fuerza que actúa sobre el mismo; todas son cantidades al tiempo t. Por otro lado ∆v es el cambio en la velocidad (entre t y t + dt) de la masa que deja o se agrega al sistema. Ejemplo 2.1 (Un cohete). Un cohete con masa inicial M0 y velocidad inicial V0 (respecto a un sistema de referencia inercial) se mueve quemando combustible y expulsando los gases resultantes a una velocidad constante µ respecto al cohete. Hallaremos la velocidad del cohete en función de su masa M (M < M0 ). Partiremos de la ecuación (2.17). Como no hay fuerzas externas al sistema cohete+gases colocamos F = 0. Llamaremos M (t) a la masa del cohete, V (t) a su velocidad y v a la velocidad del combustible expelido por el cohete respecto al referencial inercial. El cambio de v entre los tiempos t y t + dt es ∆v = (V + µ) − V = µ, luego dM dV 1 dM d dV − µ ⇒ 0=M − µ = M [V − µln(M )] 0=M dt dt dt M dt dt ⇒ V − µ ln(M ) = cte = V0 − µ ln(M0 ) por lo cual V = V0 + µ ln 2.3 M M0 Momento angular y torque. En esta sección definiremos el momento angular y el torque de un sistema de partı́culas. Estas cantidades están relacionadas con la rotación del sistema como veremos mas adelante. Dada una partı́cula de masa mα definimos su momento angular en el referencial S y respecto a un punto Q como d LαQ |S ≡ mα rα,Q × (2.18) rα,Q , dt S mα rα rα,Q Q O S rQ 30 Mecánica Clásica 1 C. Di Bartolo (Julio de 2004) donde rα,Q = rα − rQ es el vector posición de la partı́cula respecto al punto Q. El punto Q no necesariamente está fijo en S. Nótese que dos referenciales con distinto origen pero que no rotan entre sı́ miden los mismos momentos angulares. Para un sistema de N partı́culas se define el momento angular del sistema respecto al punto Q como d rα,Q . LαQ |S = mα rα,Q × LQ |S ≡ dt S α=1 α=1 N N (2.19) Dejaremos como un pequeño ejercicio para el lector demostrar que si se cambia de referencial se cumple que N LQ |S = LQ |S + (2.20) mα rα,Q × WS |S × rα,Q , α=1 y si se cambia de punto de referencia entonces LQ |S = LP |S + M ( rcm,P × ṙP,Q |S + rP,Q × ṙcm,Q |S ) . (2.21) En particular, si tomamos Q como el origen y P como el centro de masa se obtiene LO |S = Lcm |S + M rcm × ṙcm |S . (2.22) A la cantidad Lcm |S cuando S es inercial se le suele llamar momento angular interno. La derivada temporal del momento angular está relacionada con el torque. El torque de un sistema de partı́culas respecto al punto Q se define por NQ ≡ rα,Q × Fα = α α rα,Q × Fαe + α rα,Q × Fα,β . El último término de esta expresión se anula si suponemos que la fuerza interna entre cualquier par de partı́culas es paralela a la lı́nea de unión entre las dos partı́culas, rα,β × Fα,β = 0 . (2.23) β rα,β mα (2.24) Esta suposición se conoce con el nombre de “versión restringida de la tercera ley de Newton”. Veamos la demostración: α,β 1 (rα,Q × Fα,β + rβ,Q × Fβ,α ) 2 α,β 1 1 (rα,Q − rβ,Q ) × Fα,β = rα,β × Fα,β = 0 . = 2 α,β 2 α,β rα,Q × Fα,β = Fα,β mβ Dinámica de un sistema de partı́culas 31 En consecuencia sólo contribuyen al torque las fuerzas externas rα,Q × Fαe . NQ = (2.25) C. Di Bartolo α A continuación veamos cómo se relacionan el torque y la derivada temporal del momento angular. Calculemos la derivada de LQ en un referencial inercial S. d mα rα,Q × r̈α,Q |S = rα,Q × (Fα − mα r̈Q |S ) LQ = dt S α α luego d LQ = NQ − M rcm,Q × r̈Q |S . dt S (2.26) Esta ecuación juega un papel muy importante en la dinámica de rotación del sistema, principalmente en sistemas rı́gidos (i.e. sistemas en los cuales no varı́an las distancias relativas de las partı́culas que los componen). Como caso particular esta ecuación conduce a coincide con el centro de masa o d (2.27) LQ = NQ si Q dt está fijo en algún referencial inercial S Volveremos sobre esta ecuación cuando estudiemos la rotación de un cuerpo rı́gido. Por último para los puntos Q que se mencionan en esta última ecuación se cumple la ley de conservación del momento angular, NQ = 0 2.4 ⇔ LQ = cte . (2.28) Energı́a en un sistema de partı́culas. En esta sección definiremos para un sistema de partı́culas sus energı́as cinética, potencial y mecánica total. Energı́a cinética: La energı́a cinética de la partı́cula α-ésima se define como 1 Tα ≡ mα |ṙα |2 2 (2.29) y en un referencial inercial su variación es igual al trabajo realizado por la fuerza total que actúa sobre la partı́cula. Si llamamos Γα a la trayectoria que sigue la partı́cula en el lapso [t0 , tf ] el trabajo total sobre la partı́cula es 32 Mecánica Clásica 1 WΓα (Fα ) ≡ C. Di Bartolo (Julio de 2004) Γα Fα · drα = mα r̈α · drα = Γα tf d 1 mα |ṙα |2 · dt dt 2 t0 luego WΓα (Fα ) = Tα . (2.30) La energı́a cinética del sistema se define como la suma de las contribuciones individuales T ≡ 1 α 2 mα |ṙα |2 (2.31) y debido a (2.30) cumple que su variación es igual al trabajo neto realizado por todas las fuerzas internas y externas que actúan sobre el sistema, esto es ∆T = W Ext +W Int = α Γα Fαe · drα + α,β Γα Fα,β · drα . (2.32) Por último dados tres referenciales S, S y S dejamos al lector la demostración de las siguientes dos relaciones: T |S = T |S + WS |S · Lcm |S + 1 mα |WS |S × rα |2 2 α si O = O = c.m. (2.33) y 1 T |S = T |S + M |vcm |2 2 si O = c.m. y WS |S = 0 (2.34) donde vcm es la velocidad del centro de masa en el referencial S. Al término M |vcm |2 /2 se le llama energı́a cinética de traslación del sistema en el referencial S y a T |S energı́a cinética de rotación del sistema en S. Energı́a potencial. Una fuerza F (r) es conservativa si y solo si existe una función potencial V (r) que satisfaga F (r) = −∇V (r) (2.35) donde ∇ es el operador gradiente respecto a la coordenada r. Si llamamos ·Γ y Γ· a los puntos inicial y final de una curva Γ entonces el trabajo realizado por la fuerza a lo largo de dicha curva es C. Di Bartolo 33 Dinámica de un sistema de partı́culas WΓ (F ) ≡ Γ F (r) · dr = V (r ·Γ) − V (r Γ·) = −∆V . (2.36) Las fuerzas que actúan sobre cada partı́cula del sistema pueden ser conservativas y no conservativas. Llamaremos Vα (rα ) a la energı́a potencial asociada a la fuerza conservativa, neta y externa que actúa sobre la partı́cula α-ésima. El trabajo realizado por las fuerzas externas es entonces W Ext = α Γα Fαe · drα = −∆V Ext + W ExtNC , (2.37) donde Γα es la trayectoria recorrida por la partı́cula α-ésima, V Ext es el potencial externo neto dado por Vα (rα ) , (2.38) V Ext ≡ α y W ExtNC es el trabajo realizado por las fuerzas externas no conservativas. A continuación trataremos el problema de calcular el trabajo realizado por las fuerzas internas. En primer lugar si suponemos que el espacio es homogéneo debe cumplirse que las fuerzas internas son solamente función de las posiciones relativas de las partı́culas, i.e., Fα,β (rα , rβ ) = Fα,β (rα − rβ ). En segundo lugar supondremos que las fuerzas internas son conservativas (2.39) Fα,β (rα,β ) = −∇αβ Vαβ (rα,β ) , donde Vαβ es el potencial asociado a la fuerza Fα,β y ∇αβ es el gradiente en las coordenadas rα,β . Como Fα,β = −Fβ,α entonces −∇αβ Vαβ (rα,β ) = ∇βα Vβα (rβ,α ) = −∇αβ Vβα (rβ,α ) luego Vαβ (rα,β ) = Vβα (rβ,α ) . (2.40) También supondremos que se satisface la versión restringida de la tercera ley de Newton, luego Fα,β (rα,β ) = Fα,β (rα,β )ûα,β . mα (2.41) rα,β Si expresamos en coordenadas esféricas el gradiente en (2.39) y lo comparamos con la expresión anterior se obtiene que Vαβ (rα,β ) = Vαβ (rαβ ) con rαβ = |rα,β | . (2.42) Esto significa que Fα,β = Fα,β (rαβ )ûα,β es una fuerza central. mβ ûα,β 34 Mecánica Clásica 1 C. Di Bartolo (Julio de 2004) Calculemos el trabajo realizado por las fuerzas internas cuando el sistema evoluciona desde el instante t0 al instante tf . Int W0→f luego 1 tf = · drα = Fα,β · ṙα dt = [Fα,β · ṙα − Fα,β · ṙβ ] dt 2 t0 αβ Γα t0 α αβ 1 tf 1 tf = Fα,β (rα,β ) · ṙα,β dt = − ∇αβ Vαβ (rα,β ) · ṙα,β dt 2 t0 αβ 2 t0 αβ FαInt tf Int = −∆V Int W0→f (2.43) con V Int el potencial interno neto definido como 1 Vαβ (rαβ ) = Vαβ (rαβ ) . V Int ≡ 2 αβ α β<α (2.44) Note que hay una única contribución a la energı́a potencial por cada par de partı́culas. Al introducir (2.43) y (2.37) en (2.32) se obtiene que con ∆E = W ExtNC (2.45) E ≡ T + V Ext + V Int (2.46) la energı́a del sistema. Nótese que si todas las fuerzas son conservativas entonces la energı́a se conserva, i.e., ∆E = 0. 2.5 Problema de dos cuerpos. Un problema de gran importancia es el del movimiento de dos partı́culas que interactúan entre sı́. Este problema se simplifica considerablemente si lo descomponemos en el movimiento del centro de masa y el movimiento relativo de una partı́cula respecto a la otra. Hagamos el cambio de variables (r1 , r2 ) → (rcm , r) (2.47) r2 M2 c.m. O rcm r1 r M1 con M1 r1 + M2 r2 M1 + M 2 r = r1,2 = r1 − r2 . rcm = (2.48) (2.49) C. Di Bartolo Dinámica de un sistema de partı́culas 35 Es inmediato que µ r M1 µ r2 = rcm − r M2 M 1 M2 , µ ≡ M1 + M 2 r1 = rcm + (2.50) (2.51) (2.52) donde µ se conoce como la masa reducida del sistema. Las ecuaciones de movimiento de las dos partı́culas son M1 r̈1 = F1e + F1,2 M2 r̈2 = F2e + F2,1 . (2.53a) (2.53b) A su vez estas ecuaciones conducen a ecuaciones de movimiento para las variables rcm y r, M r̈cm = F1e + F2e e F1 F2e µr̈ = F1,2 + µ − M1 M2 (2.54a) (2.54b) donde M = M1 + M2 es la masa total del sistema. Los conjuntos de ecuaciones (2.53) y (2.54) son equivalentes. La ecuación diferencial (2.54a) proporciona la evolución del centro de masa del sistema y la ecuación (2.54b) nos dice como evoluciona el vector r visto desde un sistema de referencia inercial. Esta última ecuación admite otra lectura; nos dice cómo es la trayectoria de M1 vista por un observador no inercial con origen en M2 pero sin rotación respecto a los sistemas inerciales. Nótese que si F1e Fe = 2 M1 M2 (2.55) el problema se convierte en el de dos cuerpos que no interaccionan. Esto puede suceder si, por ejemplo, no hay fuerzas externas o si M1 y M2 sienten sólo la fuerza gravitatoria de partı́culas muy alejadas comparadas con r. En este caso se cumple que M r̈cm = F Ext µr̈ = F1,2 . (2.56) (2.57) Por último evaluemos en las nuevas variables la energı́a cinética y el momentum angular del sistema. Llamaremos S a un sistema de referencia inercial cualquiera y S al referencial solidario con el centro de masa pero tal que WS |S = 0. Se tiene entonces que 1 1 1 T |S = M1 (ṙ1, cm )2 + M2 (ṙ2, cm )2 = µ|ṙ|2 2 2 2 (2.58) 36 Mecánica Clásica 1 C. Di Bartolo (Julio de 2004) 1 1 1 T |S = T |S + M |ṙcm |2 = µ|ṙ|2 + M |ṙcm |2 . 2 2 2 (2.59) La energı́a cinética aparece como la suma de las energı́as cinéticas de dos partı́culas ficticias, una de masa M en el centro de masa y otra de masa µ con vector posición r. De manera similar para el momento angular también se encuentra que se puede escribir como la suma de los momentos angulares de las dos partı́culas ficticias, esto es Lcm = µ r × ṙ (2.60) L0 = Lcm + M rcm × ṙcm = µ r × ṙ + M rcm × ṙcm . (2.61) Finalmente de (2.27) se obtiene la ecuación que obedece la derivada del momento angular d d Ext . Lcm = (µ r × ṙ) = Ncm dt dt S 2.6 (2.62) Resortes ideales. En esta sección estudiaremos parcialmente el problema de dos partı́culas interactuando por medio de un resorte ideal. Un estudio completo tendrá que esperar al capı́tulo sobre fuerzas centrales. Un resorte ideal de constante elástica K y longitud natural l es un dispositivo sin masa que aplica en sus extremos una fuerza restitutiva proporcional a la diferencia de su longitud r = |r1 − r2 | y su longitud natural, i.e., F2,resorte = +K(r − l) ûr F1,resorte = −K(r − l) ûr , (2.63) (2.64) r 2 K r2 1 ûr r1 donde r (2.65) con r ≡ r1,2 = r1 − r2 . r Si lo deseamos podemos tomar al resorte como parte de una de las partı́culas o cómo un tipo de fuerza de interacción entre ellas y escribir simplemente ûr = F1,2 = −F2,1 = −K(r − l) ûr . (2.66) Debe tenerse en cuenta que esta fuerza siempre es paralela a la lı́nea de unión entre las dos partı́culas,“los resortes ideales no se doblan”. C. Di Bartolo Dinámica de un sistema de partı́culas 37 A continuación demostremos que la energı́a potencial asociada a un resorte es 1 V = K(r − l)2 , 2 (2.67) F1,2 = −∇r V (r) . (2.68) esto significa que Para demostrarlo usamos la expresión del gradiente en esféricas, ∂ µ̂θ ∂ µ̂ϕ ∂ 1 −∇V (r) = − µ̂r + + K(r − l)2 = −K(r − l) ûr . ∂r r ∂θ r senθ ∂ϕ 2 Si este sistema está aislado su energı́a se conserva y de acuerdo a (2.46) y (2.59) viene dada por 1 1 1 E = µ|ṙ|2 + M |ṙcm |2 + K(r − l)2 , (2.69) 2 2 2 con µ la masa reducida y M la masa total del sistema. También se conserva su momentum por lo cual su centro de masa se mueve con velocidad constante, P = M Vcm = constante . (2.70) De acuerdo a (2.62) otra cantidad que se conserva es el momento angular en el referencial centro de masa, (2.71) Lcm = µ r × ṙ = constante . Las consecuencias de esta última ecuación las estudiaremos en el capı́tulo de fuerzas centrales, sin embargo mencionaremos que una de ellas es que las dos partı́culas se mueven en un plano constante. La forma de tratar en general con la ecuación de movimiento del vector r también será abordada en el capı́tulo de fuerzas centrales. Nosotros trataremos ahora sólo con el caso particular de las dos partı́culas moviéndose sobre una lı́nea recta que pasa por el centro de masa y tiene dirección ûr = û =constante. En este caso la ecuación 2.57 conduce a µr̈ + K(r − l) = 0 . (2.72) Esta ecuación se conoce como la ecuación de movimiento de un oscilador armónico unidimensional y su solución general es r − l = A cos(w t + θ0 ) (2.73) donde A y θ0 son constantes de integración y w= K . µ El movimiento relativo de las partı́culas es oscilante con un perı́odo τ = 2π/w. (2.74)
© Copyright 2024