Apunte UChile - Sistemas Newtonianos
Apunte UChile - Sistemas Newtonianos
Apunte UChile - Sistemas Newtonianos
Departamento de Física
Facultad de Ciencias Físicas y Matemáticas
Universidad de Chile
Marzo 2010
Indice
I Información General 9
1.1. Presentación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.2. Programa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3. Material Docente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3.1. Material teórico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.2. Material complementario . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.3. Guía de prácticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.4. Guía de ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.5. Advertencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.4. Asistencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.5. Informes de Prácticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.6. Gráficos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.7. Sala Galileo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
II Material Docente 17
1
Información General Sistemas Newtonianos 2
4C.1.Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 214
4C.2.Guía Práctica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 214
8
Parte I
Información General
9
Información General Sistemas Newtonianos 10
1.1. Presentación
La física es una disciplina inherentemente vinculada a la experimentación y a la cuantifica-
ción de sus resultados. Para esto diseñar sistemas controlados es de suma importancia. Por
otra parte, la caracterización y síntesis de los fenómenos físicos mediante las matemáticas
formales también ha sido de crucial importancia para poner a prueba el potencial predictivo
de las teorías y modelos. En la asignatura de Sistemas Newtonianos (SN) se busca un acer-
camiento integrado de ambos aspectos, donde coexisten la experimentación, el formalismo
y el tratamiento numérico de sistemas mecánicos.
El esquema de trabajo semanal de SN consiste en una sesión de cátedra, una clase auxiliar
y una sesión práctica de tres horas de duración.
Las cátedras y clases auxiliares se realizarán en salas estándares. La sesión práctica se
realizará en en la Sala Galileo (segundo piso, lado oriente del Departamento de Física), bajo
la dirección del profesor de la cátedra y el apoyo de dos profesores auxiliares. Esta sala
cuenta con recursos multimediales para la presentación de las materias, la realización de
experimentos y de cálculo numérico.
Por lo general, el desarrollo de cada unidad contempla una introducción por parte del profesor
de la cátedra, quien presentará el soporte teórico del tema a tratar. Además supervisará la
realización de una práctica experimental y/o numérica que culmina con la elaboración de
un informe. En la realización de las prácticas, cada sección será estructurada en equipos de
trabajo de tres estudiantes, los cuales serán definidos previamente por el profesor.
1.2. Programa
Semana Unidad Contenidos
1 0: Introducción, Matlab Uso de Matlab.
Uso de Matlab, solución de problemas alge-
2 1: Métodos Numéricos
braicos, solución de ecuaciones diferenciales.
Mediciones, promedio, error, desviación están-
3 2: Métodos Experimentales
dar, histogramas. Uso del sensor de fuerzas.
Produto cruz, torque y momento angular.
4 4A: Estática de Sólidos Centro de masas. Leyes de la estática. Apli-
cación de las leyes de estática a sólidos.
Momento de Inercia. Teoremas de Steiner
5 4B: Dinámica Plana de Sólidos
Energía. Uso de webcam para medir tiempos
6 4C: Dinámica Plana de Sólidos Rotación de Sólidos. Momento angular.
Rotacion con traslación, movimiento en plano
7 4D: Dinámica Plana de Sólidos
inclinado.
Introducción a tipos de oscilaciones. Defini-
8 5A: Oscilaciones
ción y medición de periodo, amplitud y fase.
Oscilación amortiguada y caída de un cuerpo
9 5B: Oscilaciones
con roce.
10 5C: Oscilaciones Resonancia
Introducción a medios continuos, ondas pro-
11 6A: Ondas
pagativas.
12 6B: Ondas Modos normales. Resonancia.
Presión. Experimento de presión colisional.
13 7A: Fluidos
Identificación microscópica de la presión.
14 7B: Fluidos Leyes de Pascal y Arquímides.
15 Examen de Laboratorio Examen individual de técnicas de laboratorio.
Material teórico
Material complementario
Guía de prácticas
Guía de ejercicios
Este material complementa el Material Teórico con lecturas más avanzadas, material espe-
cífico de algunos temas o guías de problemas propuestos.
Esta guía describe las actividades prácticas (laboratorios y/o cálculos numéricos) que se van
a realizar en la clase.
La lectura de este material también será evaluada en el Control de Lectura al inicio de cada
sesión práctica.
Es una guía de problemas asociada a cada unidad, indicando el nivel y tipo de problemas
que serán evaluados en el Ejercicio Semanal. Este ejercicio será tomado al inicio de las clase
auxiliar que se realice inmediatamente después a la cátedra respectiva.
1.3.5. Advertencia
Es importante recalcar que NO toda la materia del curso se encuenta en estos apuntes.
Los alumnos deberán complementar estos apuntes con las lecturas recomendas (libros de
texto como el Tipler o Serway) donde la materia se expone en muchos casos con mayor
profundidad y detalle.
1.4. Asistencia
La asistencia a las sesiones prácticas es obligatoria. Una inasistencia mayor al 30 % es
causal de reprobación del curso.
El inicio de las sesiones está programada para las 08:30 hrs en la mañana, y 14:30 en
la tarde. El acceso al recinto se cierra a esta hora, permitiéndose la entrada 20 minutos
después. Posterior a ello no se permitirá el acceso al recinto de clases.
Descripción Se describe en algún detalle los pasos y protocolos seguidos y las elecciones
de parámetros o valores tomados.
Criterio de evaluación: Una correcta exposición le permitiría reproducir el experimento
a cualquier persona.
Resultados, análisis y discusión Se presenta los datos obtenidos y los graficos respecti-
vos. Se realiza además un análisis respecto a los posibles errores y la consistencia con
la teoría. Se plantaen posibles caminos para corregir las falencias, se refutan o corrigen
supuestos, etc.
Criterio de evaluación: Una correcta presentación de resultados indica los valores de
las medidas y sus desviaciones estándar o errores. Los gráficos deben indicar los ejes
y unidades y deben estar en las escalas adecuadas (ver sección Guía sobre Gráficos).
Por último, un buen análisis y discusión de los resultados permitiría comprender si se
han cumplido los objetivos de la experiencia, si los resultados son consistente y si hay
alguna dificultad propia a la actividad.
1.6. Gráficos
Los gráficos, para que éstos sean útiles para representar información física, deben tener cierta
estructura mínima.
Los gráficos deben indicar claramente las magnitudes que se están graficando (las etiquetas
de los ejes) y las unidades en las que están medidas. Además los ejes deben estar en una
escala que permita leer e interpretar correctamente los datos.
Para ilustrar lo anterior, a continuación se presenta un gráfico con toda la información
relevante:
60
50
40
v [m/s]
30
BUENO
20
10
0
10 12 14 16 18 20 22
t [s]
donde se indica que se grafica la velocidad (en metros por segundo) en función del tiempo
(en segundos). Los puntos indican las mediciones hechas y la línea es simplemente una guía
visual para leer mejor el gráfico.
A contiunuación se presentan algunos gráficos con algunas deficiencias que se deben evitar:
60
50
40
Se indica que se grafica la velocidad pero
30
no se señalan las unidades. Tampoco se
v
10
0
10 12 14 16 18 20 22
1060
1000
1050
800
1040
v [m/s]
v [m/s]
600
MALO 1030
400 BUENO
1020
200 1010
0 1000
10 12 14 16 18 20 10 12 14 16 18 20
t [s] t [s]
H G F E
MESA PROFESOR
I J
A B C D
Material Docente
17
Unidad 1: Métodos Numéricos
1.1. Introducción
En la descripción cuantitativa de los fenómenos físicos, no siempre es posible resolver aná-
liticamente el modelo que lo describe e, incluso en algunos casos en que esto es posible, la
solución analítica no siempre es fácil de interpretar y visualizar.
Como una herramienta complementaria a las técnicas analíticas, es usual en la física utilizar
las herramientas numéricas, que con la ayuda de un computador permiten resolver, analizar
y graficar diversos modelos físicos.
En esta unidad se ilustrará cómo se utiliza el computador, mediante el programa Matlab
con estos fines.
Ejemplo:
Consideremos el ejemplo de una persona parada sobre
una tarima. La persona está sujeta a dos resortes de
largo natural L, cada uno de los cuales está amarrado
a la misma altura de la persona en direcciones opuestas
y que inicialmente no están elongados y en posición h
horizontal.
En cierto momento la personal se suelta de la tarima y
cae debido a la fuerza de gravedad. Se desea encontrar
la profundidad hasta la que cae.
18
Métodos Numéricos Sistemas Newtonianos 19
Una vez que la persona se lanza, adquiere energía potencial gravitatoria, elástica y cinética.
En la elongación máxima la energía cinética se anula y su energía es:
p
E = −mgh + k( h2 + L2 − L)2
Si esta expresión se expandiera daría lugar a una ecuación de tercer grado lo que no permite
realizar un análisis simple de la solución.
Sin embargo, si se dan valores numéricos es posible obtener el valor de h pedido. Por ejemplo,
si m = 80kg, L = 10m, k = 10N/m y g = 9,8 m/s2 , en Matlab se hace:
Como la ecuación es no lineal, ésta puede tener muchas soluciones. Para escoger la que uno
desea debe darse una adivinanza inicial. Matlab va a buscar la solución más cercana a la
adivinanza inicial, de manera que es conveniente usar una buena adivinanza.
También se podría preguntar como depende h de la masa m, manteniendo los otros valores
fijos. Para eso, lo mejor es hacer un gráfico de la forma:
>> % Se puso punto y coma (;) para que no despliegue los resultados:
>> % las masas y hasol son arreglos y pueden ocupar mucha pantalla
Tarea: Ejecute este programa y analice si los resultados obtenidos son razonables.
Si t ≈ t0 ⇒ x(t) ≈ x(t0 )
Luego, si se quiere representar una función x(t) en un intervalo [TA , TB ] lo primero que se
hace es discretizar el tiempo usando un espaciado pequeño ∆t, de manera que
ti = TA + i × ∆t ; con i = 0, 1, 2, . . .
Luego, la función x(t) se representa de manera discreta, dando los valores de ésta sólo en
los instantes discretos
xi = x(ti )
t0 t1 ti−1 ti ti+1 t
∆t
Entonces, cuando se quiera resolver la ecuación de Newton en el computador, en vez de
buscar la función completa x(t), vamos a buscar los valores de xi = x(ti ) para un ∆t
dado. Es importante destacar que mientras más pequeño sea el valor de ∆t, más fiel será la
representación de la función; pero de todos modos no es necesario exagerar pues las funciones
son continuas.
Hay que notar que esta representación discreta de las funciones es lo que usualmente se hace
cuando se miden magnitudes físicas que varían en el tiempo o el espacio. Así, por ejemplo, la
temperatura en Santiago se mide una vez cada minuto y no de manera continua. Lo mismo
con la posición de un objeto móvil que al filmarlo con una cámara se tiene un muestreo cada
cierta fracción de segundos solamente.
t0 t1 ti−1 ti ti+1 t
∆t
t0 t1 ti−1 ti ti+1 t
∆t
Promediando las dos expresiones, de la derivada hacia adelante y hacia atrás, se obtiene la
llamada derivada centrada
1 xi+1 − xi xi − xi−1
ẋ(ti ) ≈ +
2 ∆t ∆t
xi+1 − xi−1
≈ (1.3)
2∆t
x(ti + ∆t) − x(ti − ∆t)
≈ (1.4)
2∆t
t0 t1 ti−1 ti ti+1 t
2∆t
De las tres expresiones para la derivada, si ∆t está fijo, se puede demostrar que la derivada
centrada es más precisa que las otras, pero a veces por comodidad o facilidad de cálculo se
usarán las otras dos. Esto se verá con detalle en el curso Métodos Numéricos.
A su vez,
x(t + ∆t) − x(t)
ẋ(t + ∆t/2) =
∆t
x(t) − x(t − ∆t)
ẋ(t − ∆t/2) =
∆t
Reemplazando en la expresión anterior se obtiene
x(t+∆t)−x(t)
∆t − x(t)−x(t−∆t)
∆t
ẍ(t) =
∆t
x(t + ∆t) − 2x(t) + x(t − ∆t)
=
∆t2
Luego, evaluando en un punto de la discretización
xi+1 − 2xi + xi−1
ẍ(ti ) = (1.5)
∆t2
que se llama derivada centrada de segundo orden.
Tarea: Muestre que si x(t) = αt + β (α y β constantes), es decir, una línea recta, entonces
la expresión aproximada de la derivada se anula tal como la expresión exacta.
En resumen, se tienen la siguientes expresiones aproximadas para la derivadas
xi+1 − xi
ẋ(ti ) ≈ (1.6)
∆t
xi − xi−1
ẋ(ti ) ≈ (1.7)
∆t
xi+1 − xi−1
ẋ(ti ) ≈ (1.8)
2∆t
xi+1 − 2xi + xi−1
ẍ(ti ) ≈ (1.9)
∆t2
mẍ = −kx
Esta expresión indica que si se conoce la posición en un instante (xi ), en otro previo (xi−1 )
y el valor de la aceleración en ese instante (−kxi /m), entonces es posible calcular la posición
en un instante posterior. Dado lo anterior surge la idea de realizar una iteración pues con el
nuevo valor obtenido ahora será posible calcular el siguiente y así para adelante.
t0 t1 ti−1 ti ti+1 t
Como para cada paso de la iteración se necesitan dos valores previos de la posición, surge
el problema de cúales deben ser los primeros valores. Como se sabe, en la física newtoniana
se debe indicar la posición y velocidad inicial del cuerpo para poder resolver el movimiento.
Es decir, se dan x0 y v0 . Recordando que la velocidad es la primera derivada de la posición
y usando la expresión (1.6) se tiene
v0 = ẋ(0)
x1 − x0
=
∆t
de la que se despeja
x1 = x0 + v0 ∆t (1.11)
Se tienen así todos los pasos del llamado algoritmo de Verlet (L. Verlet, Physical Review
1967)
• Dados: x0 , v0
• Se calcula: x1 = x0 + v0 ∆t
• Se itera desde i = 1 hasta el tiempo final:
k
xi+1 = 2xi − xi−1 − xi ∆t2
m
>> v0=2.0
>>
>> % Se define el arreglo de tiempos
>> t=0:dt:Tfin;
>> % Se define el arreglo de posiciones
>> x=zeros(1,length(t));
>> % Se ponen las condiciones iniciales
>> x(1) = x0;
>> x(2) = x0+v0*dt;
>> % Se itera usando el algoritmo de Verlet
>> for i=2:length(x)-1
>> x(i+1) = 2*x(i)-x(i-1) - k*x(i)*dt^2/m;
>> end
>> plot(t,x)
Para el caso de una fuerza cualquiera F (x) es resultado es análogo y el paso de iteración es
Si el sistema tiene fuerzas que dependen de la velocidad, por ejemplo fuerzas de roce viscosas
o turbulentas, F (x, ẋ), se debe escribir la derivada de manera discreta. Lo más práctico en
este caso es usar la derivada hacia atrás. En efecto, en ese caso se puede hacer
xi − xi−1
F (x, ẋ) ≈ F xi ,
∆t
que permite hacer el cálculo numérico pues es explícita, es decir, el valor de xi+1 sólo depende
de valores previos (que ya son conocidos).
x*
t0 t1 t2 t3 t4 t5 t6 t7 t8 t
Una vez que se sabe en qué intervalo ocurre la intersección, se debe asignar un tiempo.
Dado que la función se conoce sólo con una precisión ∆t, la respuesta también tendrá esta
indeterminación. Así, las siguientes respuestas son igualmente válidas.
t∗ ≈ ti (1.15)
∗
t ≈ ti+1 (1.16)
∗
t ≈ (ti + ti+1 )/2 (1.17)
Una mejor estimación se puede obtener haciendo una interpolación lineal, pero esto se dejará
para más adelante.
¿Cómo graficar dos o más curvas con una sola instrucción plot?
Utilice plot(x1,y1,’*’,x2,y2,’-’,x3,y3,’k–’), donde los símbolos entre comillas
son opcionales.
z=x./y: define el arreglo z formado por el cuociente de los elementos de x con los de
y.
Pregunta 3: Si x(t) = sin(t) exp(−0,1t), escriba las líneas de Matlab que permitirían gra-
ficar esa función en el intervalo t ∈ [0 : 10] usando un espaciado ∆t = 0,1.
Pregunta 5: Si x(t) = cos(t) exp(−2t), escriba las líneas de Matlab que permitirían graficar
esa función en el intervalo t ∈ [0 : 4] usando un espaciado ∆t = 0,1.
Pregunta 7: Suponga que tiene dos arreglos de N elementos, T y X. T (i) indica el tiempo
i-ésimo y X(i) la posición de una partícula en dicho instante. Escriba un programa
que encuentre el instante aproximado en que la partícula cruza en origen.
ẍ + cẋ + kx = 0
Use el método de Verlet visto en clases para resolver estas ecuaciones acopladas. Con-
sidere los siguientes valores de las constantes: G = 1 y M = 1. Además considere como
Ejercicio 7: Se desea modelar un billar donde las bolas se mueven en un plano con roce
viscoso sobre la superficie y rebotan contra las paredes de manera elástica (es decir, el
ángulo de entrada es igual al ángulo de salida en el rebote).
Para simplificar el ejercicio, se considerá sólo una pared (en vez de 4 que tiene el billar).
Esta pared es horizontal y está en y = 0.
Ejercicio 8: Las nuevas micros del TranSantiago dispondrán de GPS, aparato que les entre-
ga la posición de la micro (x, y) cada cierto intervalo de tiempo ∆t. Se desea incorporar
a las micros un mecanismo de control que, usando los datos del GPS, permita medir
la velocidad ~v y la aceleración ~a en cada instante de manera que suene una alarma si
|~v | > Vc ó |~a| > Ac .
Considere por simplicidad que el movimiento es puramente bidimensional (es decir,
Santiago es plano).
Como el sistema tiene roce, si el péndulo se suelta del reposo desde un ángulo inicial
φ0 , los ángulos máximos que alcance (indicados por un círculo en la figura) serán cada
vez menores.
Se busca resolver numéricamente la dinámica del sistema para obtener cómo van dis-
minuyendo estos ángulos máximos. Para eso:
Ejercicio 10: Se desea calcular el alcance de un proyectil que se lanza, desde el nivel del
suelo, en un ángulo α respecto a la horizontal con velocidad V0 . Sobre el proyectil
actúa el roce con el aire, que le ejerce una fuerza F~ = −γ~v .
1. Escriba la parte del programa Matlab que permita calcular la trayectoria (x(t) e
y(t)) del proyectil. No olvide poner las condiciones iniciales.
2. Escriba la parte del programa Matlab que entregue el alcance del proyectil.
2.1. Introducción
La física es una ciencia natural, y por ende el rol de la experimentación es clave. Toda ciencia
natural se basa en el método científico, el cual puede resumirse brevemente en tres pasos: 1)
Se realizan observaciones de un fenómeno natural. 2) Mediante un proceso de razonamiento
se postulan hipótesis que explican los fenómenos observados. Con estas hipótesis es posible
predecir otros fenómenos. 3) Se realizan experimentos para verificar las predicciones. En el
caso de que éstas no se verifiquen, los nuevos resultados experimentales sirven de observación
(se vuelve al paso 1) para poder plantear nuevas hipótesis.
El proceso de medida es también esencial en la Ingeniería, pues en muchas aplicaciones o pro-
cesos es necesario controlar o monitorear un grupo de cantidades físicas, ya sean mecánicas,
térmicas o eléctricas, como también cantidades químicas o biológicas.
En esta unidad se presentarán métodos experimentales básicos, en particular la medición de
cantidades físicas como fuerza, posición, velocidad y aceleración. Se presentarán los conceptos
básicos de un análisis estadístico y se discutirán los tipos genéricos de errores de medición.
Se introducen luego algunas instrucciones de Matlab que le permitirán realizar un análisis
estadístico básico de un conjunto de medidas. Finalmente se dan detalles del uso de la tarjeta
de adquisición de datos que se usará en este curso. En el curso Métodos Experimentales se
extregarán herramientas experimentales más avanzadas.
35
Métodos Experimentales Sistemas Newtonianos 36
Figura 1: Sensor de fuerza marca Vernier, modelo Dual-Range sensor. La fuerza puede ser
aplicada en ambos sentidos en la dirección del eje que se muestra en la foto.
tiempo respectivamente. Una manera simple de medir fuerzas es usando masas conocidas.
Por ejemplo, si se quiere aplicar tensión a una cuerda podemos colocar astutamente una
masa m a un extremo la cuerda. Si el sistema es estático, entonces las leyes de Newton nos
dicen que la tensión, una fuerza, es igual a mg. También se puede usar un Dinamómetro,
el cual mediante una escala y una aguja muestra la medida de fuerza tal como una balanza
casera.
Otra manera de medir una fuerza, la cual aplicaremos en este curso, es mediante un sensor de
fuerza. El tipo de sensor que usaremos, llamado strain gage en inglés, consta de una pequeña
y delgada lámina metálica la cual es deformada al aplicar una fuerza. Su deformación consiste
básicamente en un cambio de espesor de la película, lo que induce un cambio de resistencia
eléctrica de ella. El sensor entrega entonces un voltaje proporcional al cambio de resistencia,
el que a su vez es proporcional a su espesor. El cambio de espesor es a su vez proporcional
a la fuerza aplicada (Ley de Hooke, o elasticidad lineal). Entonces, finalmente, el voltaje
medido es proporcional a la fuerza aplicada.
En la sección Material Docente de U-Cursos podrá encontrar el archivo dfs-bta.pdf donde se
dan los detalles técnicos del sensor de fuerza. El sensor tiene dos rangos de medición, ±10
N y ±50 N, a escojer según la aplicación. La salida de voltaje del sensor varía en el rango 0
a 5 V (en estas expresiones V corresponde a la unidad de Volt y N a Newton). La relación
entre fuerza y voltaje, es decir la calibración, es la siguiente:
F = A · U + B,
donde F es la fuerza y U el voltaje medido. En el rango ±10N del sensor, se tiene por defecto
los siguientes parámetros de calibración: A = −4,9 N/V y B = 12,25 N. Para el rango ±50
N, se tiene A = −24,5 N/V y B = 61,25 N. Ambas situaciones implican que para F = 0,
entonces el voltaje medido será U = 2,5 V. El sensor tiene una sensibilidad de 0,01 N y 0,05
N en los rangos ±10 N y ±50 N respectivamente. Atención que el sensor es suficientemente
sensible como para detectar el peso del gancho que usa. Para una configuración dada (hori-
zontal o vertical) se puede imponer F=0 para el voltaje medido que puede ser ligeramente
diferente a 2,5 V. En todo caso, cada sensor tiene constantes de calibración ligeramente
diferentes a las nominales. Estrictamente uno debería calibrar independientemente cada uno
de ellos. Sin embargo, en este curso usaremos las constantes de calibración indicadas por el
fabricante.
La medida de voltaje que entrega el sensor de fuerza será realizada con una tarjeta de
conversión análoga/digital la cual se contecta a un computador mediante un puerto USB.
Esta tarjeta tiene varios canales (vias) de medición y daremos más detalles de su uso más
adelante (sección 2.8).
Instrumentos simples para medir posición y tiempo son una regla y un cronómetro. Pero si
el movimiento es muy rápido, estos elementos en general no son suficientemente precisos.
En este curso usaremos una cámara web, a partir de la cual obtendremos una secuencia
de imágenes, ya sea un video o un conjunto de fotos. Las imágenes deberán ser calibradas,
por ejemplo tomando una foto de una regla y obteniendo una conversion entre pixeles y
centímetros. Con esta cámara se medirá entonces la distancia recorrida por un objeto entre
dos cuadros consecutivos, con lo cual se puede definir una velocidad. En otras aplicaciones
se medirá el ángulo recorrido por un objeto entre dos cuadros consecutivos, obteniéndose
una medida de velocidad angular.
Pregunta 1: ¿A partir de un conjunto de imágenes consecutivas, cómo determinaría usted
la aceleración de un objeto? ¿Cuántas imágenes consecutivas necesitaría?
30 1400
N = 100 N = 10000
25 1200
1000
20
Ocurrencia
Ocurrencia
800
15
600
10
400
5 200
0 0
6 8 10 12 14 6 8 10 12 14
C C
De los datos de la figura 26 con N = 100 medidas se obtiene hCi = 10,05 y σ = 1,12, y con
los datos con N = 10000, hCi = 10,01 y σ = 1,00. Se observa que las diferencias obtenidas
no son muy significativas, pero evidentemente mientras más datos o medidas se tenga de
una cantidad C, mejor será su estimación. En particular, la distribución es evidentemente
más simétrica entorno al valor medio para el caso con N = 10000 medidas.
Para terminar daremos una interpretación del parámetro σ. En el caso de una distribución
Un ejemplo de error sistemático es la medición de una longitud utilizando una regla plástica
que ha sido deformada, por ejemplo estirada en una dirección. Entonces, al medir una
longitud con esta regla, se obtendrá un valor sistemáticamente menor al valor real. Otro
ejemplo es la medición del ángulo de un plano inclinado. Supongamos que este plano esta
puesto sobre una mesa, con un cierto ángulo α respecto al plano de la mesa. Ya sea mediante
el uso de un transportador o de una imagen obtenida con una cámara se puede determinar
el ángulo α. Pero si se trata de usar este ángulo para comprobar alguna predicción realizada
con el plano inclinado, es muy probable que encontremos un error. La razón de esto es
que el ángulo a utilizar en nuestro cálculo es el ángulo con respecto a la horizontal que
es perpendicular a la dirección de gravedad. Si la mesa no se encuentra bien nivelada con
respecto a esta horizontal, es decir hay un ángulo β 6= 0 entre la horizontal y el plano de la
mesa, entonces nuestra medida presentará un error sistemático. La forma de eliminar este
error es simplemente nivelando bien la mesa (tratando de obtener β = 0), o bien, en lugar
de usar el plano de la mesa como referencia se puede usar una plomada, la cual definirá la
dirección vertical de la gravedad.
Los errores aleatorios se pueden presentar por un problema de medición, como por la pre-
cisión de un instrumento, o bien porque la cantidad misma a medir tiene una variación
intrínseca. Por ejemplo, si se quiere medir una barra metálica de longitud L = 12,46 cm con
una regla común, lo más probable es que ciertas veces midamos 12,4 cm, otras veces 12,5, y
tal vez interpolando la graduación milimétrica de la regla, algunas veces 12,45 cm. Si hace-
mos un gran número de medidas, lo más probable es que al promediarlas obtengamos 12,45,
lo que es una muy buena estimación del valor real. Otro ejemplo es la determinación de la
altura de los niños de un curso de primero básico. La altura de éstos varía intrinsecamente,
no por un problema de medición si no porque simplemente es imposible que todos midan
lo mismo. Una representación estadística nos dirá cual es el valor medio, supongamos 1,14
m, con una desviación estándar de 12 cm (0,12 m). Pero si tomamos un niño para medir su
altura, este valor no será una buena estimación. Lo mismo con el ejemplo de la medida de
la longitud L de la barra metálica.
Error absoluto: se asocia como error absoluto de una medida el valor de la desviación
estándar σ. La cantidad C queda representada entonces por
C = hCi ± ∆C,
con ∆C = σ. En términos vagos, esto quiere decir que C toma valores típicamente
(un 68 % de las veces) entre hCi − ∆C y hCi + ∆C.
Suma: c = a + b
p
c = hci ± ∆c = (hai + hbi) ± (∆a)2 + (∆b)2 ,
Resta: c = a − b
p
c = hci ± ∆c = (hai − hbi) ± (∆a)2 + (∆b)2 ,
Multiplicación: c = a · b
s 2 2
∆a ∆b
c = hci ± ∆c = (hai · hbi) ± (hai · hbi) + ,
hai hbi
División: c = a/b
s 2 2
hai hai ∆a ∆b
c = hci ± ∆c = ± + ,
hbi hbi hai hbi
Se debe notar que al sumar o restar valores, el error absoluto siempre aumenta. Por ejemplo:
(2,0 ± 0,3) − (1,9 ± 0,4) = 0,1 ± 0,5
Tarea: Explique por qué este resulltado es razonable.
En otras situaciones se tiene que aplicar una función f (x) a una cantidad representada por
a = hai ± ∆a. Para ello se aplica la siguiente fórmula:
Función f (x):
df
f (a) = f (hai ± ∆a) = f (hai) ± · ∆a,
dx x=hai
es muy buena, por motivos obvios. Se habla entonces que “en Chile hay aprox 15 millones de
personas”, es decir 15 × 106 habitantes. Nunca se dice “en Chile hay 15.123.876 habitantes”.
La primera frase incluye una cantidad con dos cifras significativas, la segunda, con siete!
¿Cuál es correcta? ¿Cuál es la precisión? Probablemente la respuesta correcta no es ninguna
de las dos, pero se acerca bastante más a la primera.
Definiremos las cifras significativas de una cantidad como el número total de cifras que tiene
el número respecto a la posición de la cifra “importante” que tiene su error, lo que equivale de
cierta manera a su precisión. Esto quedará mas claro dando algunos ejemplos. Supongamos
que se mide una longitud con una regla y se concluye que el valor es de 14,4 ± 0,1 cm. La
cifra importante del error es el “1”, o sea los ceros no cuentan. En este caso el número 14,4
tiene tres cifras significativas. Si por algún motivo, la medida hubiese dado 14,4 ± 1 cm,
entonces lo correcto es decir que tiene dos cifras siginificativas y se deja el resultado como
14 ± 1 cm. Si en vez de 14, 4 es 14, 5, entonces el resultado ese número se aproxima a 15,
siendo el resultado final 15 ± 1 cm.
Un error clásico, que lamentablemente se “acarrea” hasta muy adelante en las carreras
científicas y de ingeniería, es que al presentar un resultado, por el ejemplo el promedio y
la desviación estándar de una cantidad, un alumno suele entregar C = 12,23453 ± 0,01236.
¿Se tiene la precisión para entregar todas esas cifras? ¿Cuáles son las cifras significativas de
este número? Lo cierto es que no se puede presentar un resultado así. Lo correcto es darse
cuenta que la cifra importante del error es el “1” en la centésima, por lo tanto el resultado
correcto es C = 12,24 ± 0,01, y el resultado tiene cuatro cifras significativas.
Cuando se realizan operaciones artiméticas con varias cantidades físicas, entonces las cifras
significativas estarán dadas por el error final, dadas por las fórmulas de la sección precedente.
Por ejemplo, al multiplicar 14, 4 ± 0, 1 por 5, 3 ± 0, 1, para calcular un área por ejemplo, se
obtiene 76,32... ± 1,53..., donde los “...” representan más digitos. Lo correcto en este caso es
dejar el resultado como 76 ± 2, siendo 76,3 ± 1,5 también aceptable.
Pregunta 4: Indique el número de cifras significativas y escriba la forma correcta de pre-
sentar las siguientes cantidades: a = 23,34 ± 0,19; b = 234,56 ± 12,21; a + b; a − b; a · b;a/b;
a · b2 ;
>> d = load(’datos.txt’);
Esta instrucción crea una variable llamada d, que tiene una estructura de vector de dimen-
siones 100 × 1. Una advertencia, para que Matlab lea bien los datos éstos deben tener como
separador decimal un punto y no una coma (datos del tipo 3.45 los entiende como un sólo nú-
mero pero 3, 45 no, de hecho los entiende como dos datos diferentes). A veces Windows está
configurado para grabar los números con comas, lo cual puede provocar errores en Matlab.
Esto se corrige cambiando la opción en Configuración regional de Windows.
A continuación usted puede calcular el valor medio y la desviación estándar de este conjunto
de datos usando las instrucciones
>> p = mean(d);
>> s = std(d);
>> hist(d);
Lo cual creará una figura con el histograma obtenido con 10 intervalos (llamados bins en
inglés). Si se quiere cambiar el número de intervalos a 25, se debe hacer
>> hist(d,25);
Supongamos ahora que se quiere encontrar el porcentaje de datos que se encuentran dentro
del rango
La primera línea crea una variable K que tiene como valores los índices (posiciones) del
vector de datos d que cumplen con la condición de estar en el rango [p − s, p + s].
Otras datos que pueden ser útiles son el máximo y mínimo de un conjunto de números. Estos
se obtienen con
>> max(d)
>> min(d)
Pregunta 5: Usando los datos del archivo datos.txt, ¿cuál es el porcentaje de medidas fuera
del rango [hCi − 3σ, hCi + 3σ]? ¿Coincide bien esta medida con el resultado que obtuvo para
la pregunta 2?
La tarjeta tiene una velocidad máxima de adquisición de 10000 datos por segundo,
pero esta velocidad es regulable desde el programa que se usa con la tarjeta. Normal-
mente una velocidad de 100 a 1000 datos por segundo debería ser suficiente según la
aplicación.
Hay ocho canales (vias) de entrada de señal análoga, las cuales se pueden usar en modo
simple (single mode en inglés) o en modo diferencial. En modo simple cada entrada
de voltaje se mide con respecto a una referencia común (tierra), idealmente igual a
0. En este caso se pueden usar efectivamente los ocho canales, es decir una señal
independiente en cada vía de entrada. En el modo diferencial se mide la diferencia
entre dos canales vecinos (canal ai0 con ai1, analog input 0 con el canal analog input
1, etc). En este caso se pueden usar hasta un máximo de cuatro pares de canales. Con
el sensor de fuerza usaremos la tarjeta en modo simple (RSE).
Existen dos salidas análogas de voltaje, con una velocidad de salida de 150 puntos por
segundo. Se pueden usar estas salidas como fuentes de voltaje, en un rango de ±5 V.
De hecho, usaremos una de estas salidas para alimentar al sensor de fuerza con un
voltaje constante de +5 V.
Pregunta 2: Describa las diferencias entre los errores de medición sistemáticos y los alea-
torios.
Pregunta 3: Explique cómo usando una regla y una cámara web se pueden medir velocidad
y aceleración de un cuerpo. Suponga que usted conoce el tiempo δt entre fotos sucesivas.
Pregunta 5: Explique cómo usando una regla y un resorte de rigidez y largo natural co-
nocidos se puede medir una fuerza.
P2. [1 pt] Se da un valor como 35,562 ± 0,2. Reescríbalo con el número adecuado de cifras
significativas. Si el valor se diera como 35,562 ± 0,7 % ¿Cómo debiera escribirse?
P3. [4 pt] Se realiza una serie de medidas de temperatura T con una termocupla a través
de una medida de voltaje, que denotaremos U . Se usa además un segundo sensor de
temperatura que entrega una medida directa en grados Celsius (◦ C). La conversión
entre U y T para la termocupla está dada por una ecuación a determinar.
El cuadro 1 tabula estas medidas. A partir de estos datos grafique la curva temperatura
versus voltaje en el rectángulo cuadriculado que se encuentra en el reverso de esta hoja.
Suguiera y justifique la ecuación que describe la conversión entre estas dos cantidades.
Medida N◦ U (Volts) T (◦ C)
1 0.18 -77.7
2 1.22 -66.3
3 1.88 -57.4
4 2.83 -45.2
5 4.19 -24.3
6 4.87 -12.9
7 6.07 10.5
8 7.23 34.7
9 8.08 52.8
10 9.19 78.4
11 9.75 93.5
Ejercicio 3: Se realiza una serie de medidas de tensión T en un hilo, cada medida con
sus respectivo error ∆T . Todos los datos se anotan en una tabla. Encuentre el valor
medio de T y el error absoluto asociado. Con los valores de tensión, realice además un
histograma con 5 intervalos (bines).
Ejercicio 4: Un sensor de fuerza tiene una relción de calibración dada por F = A·U 2 +B·U ,
donde U es el voltaje medido y F la fuerza correspondiente. El sensor puede ser utili-
zado en dos rangos de fuerza, con un módulo máximo de 5 N y 25 N respectivamente.
Describa un método que permita calibrar este sensor, es decir que permita encontrar
los valores las constantes A y B en cada rango de medida. Se sabe que estas constantes
son ambas positivas. ¿Qué diferencia este sensor del modelo que se usa en este curso?
Nota: Por ahora sólo preocúpese de los valores “promedios" de A y B, sin considerar
sus errores absolutos ∆A y ∆B. El proceso para obtener estos errores será descrito
más adelante.
Ejercicio 5: Un estudiante decide hacer una medida de la velocidad media con el cual
camina desde la estación de metro a su universidad. Para ello decide recorrer el mismo
camino durante cinco días (de lunes a viernes) y utiliza un cronómetro, con precisión
de 0,1 s, para registrar el tiempo recorrido Tr . Los resultados obtenidos se muestran
en la siguiente tabla.
Medida N◦ Tr (s)
1 605.5
2 623.3
3 598.6
4 655.3
5 640.9
Ejercicio 6: Determine si falta alguna información relevante o si hay algún error evidente
para todas las siguientes expresiones que involucran alguna cantidad física:
Considere las siguientes cantidades: a = 4,3 ± 0,1, b = 10,12 ± 0,04 y c = −6,08 ± 0,03.
Evalúe las siguientes expresiones y tenga cuidado con las cifras significativas de su
resultado final.
a/b − c
9,81 · exp(−b · c/10)
a + 3,2 · b − c
a/b + c
5,5a exp(−bc)
a ln(2b) + c
a + 7b − c
abc
a − 3 sin(bc)
Ejercicio 9: Considere los siguientes valores de velocidad entregados por una persona que
no sabe de cifras significativas ni errores. Corríjalos para que sean consistentes
3.1. Introducción
Durante las primeras aproximaciones a la física han sido expuestos a las ideas de Galileo y
Newton sobre la naturaleza del movimiento. En estas nociones introductorias se ha usado el
concepto de partícula puntual. Entendemos por partícula puntual un objeto sin dimensiones
espaciales, dotado de una posición y masa como únicos atributos físicos. Dicho concepto
es una idealización, solo aplicable, cuando en los fenómenos que estamos estudiando la
extensión espacial de los objetos puede ser despreciada. De este modo, Newton fue capaz
de explicar el movimiento de la traslación de la tierra en su conjunto, sin preocuparse de la
otros movimientos, como la rotación. Del mismo modo, ud. ha trabajado con conceptos como
la velocidad de un automovil, la aceleración de un paracaidista, etc. Todos los ejemplos son
sistemas complejos, cada cual con una estructura interna capaz de cambiar. En este capítulo
construiremos las herramientas básicas que permiten la comprensión del movimiento de los
objetos reales.
Para poder plantear el estudio de un sistema complejo, usando las herramientas que nos
entrega la mecánica Newtoniana, debemos establecer un modelo del mismo. El hacer un
modelo, implica dejar de lado una serie de fenómenos asociados con el sistema, hacer apro-
ximaciones, y por último, plantear representaciones simplificadas del sistema. Todo esto es
en principio de gran dificultad, y requiere una gran intuición, pero por sobre todo hace
necesaria la constante comprobación experimental. No es posible establecer modelos físicos
que abarquen la totalidad de los fenómenos asociados a un sistema. La validez de nuestros
modelos es limitada a ciertos aspectos y ciertos rangos de valores.
En un intento de clasificación muy general, los tipos de sistemas que podemos contemplar
son
i.- disgregados, tales como galaxias de estrellas, sistemas granulares, gases, etc.;
ii.- líquidos, donde las moléculas constituyentes mantienen cohesión, pero permiten que
dos moléculas en contacto puedan, luego de algún tiempo, estar muy distantes entre
sí.
iii.- medios elásticos, donde el conjunto de moléculas vecinas se mantiene pero sus y dis-
tancia de separación puede variar moderadamente ante deformaciones.
51
Sistemas Extendidos Sistemas Newtonianos 52
iv.- sólidos indeformables, donde las distancias entre moléculas vecinas es invariable en el
tiempo.
Figura 4: Choque de dos galaxias interretado a partir de datos obtenidos por el telecopio espacial
Hubble.
III.- La fuerza que un agente externo ejerce sobre el cuerpo es igual en magnitud, pero de
sentido opuesto, a la fuerza que el cuerpo ejerce sobre el agente.
El Principio de Superposición se refiere a que las fuerzas son aditivas en el sentido vectorial.
(a)
(b)
+ + + + + + + + + + (c)
(d)
N
X
~ = mi
R ~ri .
M
i=1
~ como el
Tan solo al mirar la definición podemos interpretar geométricamente el punto R
promedio de las posiciones de las partículas, ponderando cada una relativamente con la
contribución que hacen a la masa total. Esta interpretación en términos de promedio, permite
utilizar la intuición para identificar el CM de sistemas simples, por ejemplo el CM de un
círculo esta en su centro. En coordenadas cartesianas, si denotamos
~ = (Rx , Ry , Rz ) ,
R (3.1)
~ri = (xi , yi , zi ) , (3.2)
entonces,
N
1 X
Rx = mi xi , (3.3)
M
i=1
XN
1
Ry = mi yi , (3.4)
M
i=1
N
1 X
Rz = mi zi . (3.5)
M
i=1
111111111111111111
000000000000000000
M, L
0 xi= i d
del sistema es la suma de la contribición de cada una. Así, si la celda i-ésima es de masa
mi y su coordenada con respecto al nivel cero de energía potencial es yi , entonces la energía
potencial total Ug es
N
X PN
i=1 mi yi
Ug = mi gyi = M g = M gYCM ,
M
i=1
donde claramente YCM representa la coordenada Y del centro de masas del sistema.
Como vimos anteriormente, el centro de masas de una barra uniforme se ubica en su punto
medio. Nos preguntamos por la ubicación del centro de masas de una ’L’formada por dos
barras idénticas, unidas perpendicularmente en uno de sus extremos. Este cálculo es bastante
simple si se recurre al siguiente teorema:
~ ~
~ = MA RA + MB RB
R
MA + M B
~ A = (P
Pero, MA R ~ B . Con ello,
i∈A mi~
ri ), análogamente para R
~ = MA R
MR ~ A + MB R
~B ,
el resultado buscado.
M × (L/2) + M × 0 L M × 0 + M × (L/2) L
Rx = = ; Ry = = .
M +M 4 M +M 4
Este resultado representa el punto medio entre los puntos medios de cada barra,
como se ilustra en la figura de más abajo.
y
L 1
0
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
L/2 0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
x
11111111111
00000000000
O L/2 L
simplemente el primero desplazado hacia la derecha por a/2: CM2 = a2 + X2 , Y2 , mientras
que el tercero es el primero desplazado hacia arriba por b/2: CM3 = X2 , 2b + Y2 . Por último,
el cuarto esta invertido y desplazado hacia (a/2, b/2), CM4 = a2 − X2 , 2b − Y2 .
El CM del todo es el CM de las partes (propiedad mágica). Como cada triángulo tiene un
cuarto de la masa tenemos:
M M
4 CM1 + 4 CM2 +M4 CM3 +
M
4 CM4
CM = (3.6)
M
1
= (CM1 + CM2 + CM3 + CM4 ) (3.7)
4
usando la parte anterior, tenemos,
1
CM = (a + X, b + Y ) (3.8)
4
pero sabemos que el CM es (X, Y ), de modo que:
1
(X, Y ) = (a + X, b + Y ) (3.9)
4
Obteniendo:
1
(X, Y ) = (a, b) . (3.10)
3
Entonces
P~ = M V
~ (3.11)
Hemos denotado por V ~ la velocidad del centro de masas de un sistema. Su tasa de variación
por unidad de tiempo corresponde a la aceleración del centro de masas y la denotamos por
~a. En algunos textos se usa la notación ~aCM . Entonces,
~
dV
~a = .
dt
El uso de las leyes 2da y 3ra de Newton sobre cada una de las N componentes del sistema
permite encontrar una ecuación para el movimiento del centro de masas del sistema. Esta
ecuación resulta independiente de las fuerzas internas que ligan al sistema y del eventual
cambio de geometría que este pueda experimentar.
Consideremos un sistema formado por N partículas interactuando entre ellas. El sistema
no necesariamente es un sólido; puede tener cualquier constitución. Al analizar la i-ésima
partícula, observamos que sobre ella puede actuar resultante externa que denotamos F~i .
También, sobre la misma partícula i-ésima interactúan las restantes N − 1 partículas. Por
simplicidad, asumiremos que las fuerzas de interacción actúan entre pares de partículas. Es
fácil generalizar las conclusiones que obtendremos, si consideramos fuerzas más complejas,
pero la discusión se vería seriamente entorpecida4 . Hecho este comentario, denotemos f~j/i
la fuerza que ejerce la componente j sobre la i, naturalmente f~i/i = 0, entonces la fuerza
4
Es posible que el estudiante haya sido educado en la creencia de que las fuerzas elementales son siempre
entre pares. Por ejemplo, las leyes de gravitación Universal y de interacción electrostática son de este tipo. Es
necesario aclarar que la generalización de este principio al conjunto de todas las interacciones fundamentales
carece de un sentido no preciso, pues las leyes fundamentales de la naturaleza no se escriben en el lenguaje
entregado por las leyes de Newton. Sin embargo, y este es el principal punto de este apartado, incluso si fuese
cierto, el supuesto de interacciones de a pares no se reflejaría necesariamente en las interacciones de nuestras
partículas. Las celdas en las que dividimos los sistemas son complejas y estan compuestas de muchos átomos.
P ~
neta de todas las componentes sobre la i-ésima es N j=1 fj/i . La ecuación de movimiento de
Newton para la partícula i es:
XN
d~pi
= F~i + f~j/i . (3.12)
dt
j=1
Aquí observamos
PN d~
pi ~
dP
que la suma de derivadas es igual a la derivada de la suma, por lo cual i=1 dt = dt ;
que la suma de todos los pares de fuerzas internas es nulo, en virtud al principio de
acción y reacción (f~i/j = −f~j/i ).
Con lo anterior
dP~ d(M V ~)
~˙ = M~a .
F~ = = = MV (3.14)
dt dt
Esta ecuación resume el comportamiento del centro de masas de un sistema cualquiera
cuando sobre este actúa una fuerza externa neta F~ . Este es un resultado totalmente general,
independiente de si el sistema es sólido, líquido, gaseoso, granular, plástico, amorfo, etc. El
movimiento del centro de masas está determinado exclusivamente por las fuerzas externas.
Notemos que este resultado es lo que a priori nos permitía el uso de las Leyes de Newton
para sistemas extendidos. A lo largo de todo el primer semestre cuando hablabamos del
moviento de un cuerpo en respuesta a fuerzas, nos referiamos implicitamente al movimiento
del centro de masas del cuerpo.
j
i
fj/i
ω
ρi ωρi
1
0
0
1
Figura 8: Un cuerpo plano rotando en torno a un eje fijo; a la derecha se ilustra una vista desde
arriba.
Apéndices
1 Fx
Fx = (m + m)aCM ⇒ aCM = .
2 m
Por otro lado, al considerar sólo la partícula B, su aceleración es simplemente aB = Fx /m.
Con ello notamos que
1
aCM = aB ,
2
vale decir, el centro de masas acelera a una tasa igual a la mitad de lo que ocurre con B.
a CM aB
111
000 11
00
000
111 00
11 Fx
xA xB
Lo anterior es consistente con el siguiente resultado para la localización del centro de masas:
1
Rx = (xA + xB )
2
Si xA no es alterada por la fuerza externa, entonces
1 1
Ṙx = ẋB ⇒ R̈x = ẍB
2 |{z} 2 |{z}
aCM aB
2 πrd
1111111111111
0000000000000
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
r 0000000000000
1111111111111
r
0000000000000
1111111111111
0000000000000
1111111111111 r
0000000000000
1111111111111
2π
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
r=d/2:d:1
d franja extendida
Pregunta 6: Considere un lápiz de grafito e indique un eje que lo atraviese con respecto al
cual su momento de inercia sea máximo, y un eje con respecto a cuál sea mínimo.
Pregunta 7: Considere el sistema conformado por granos de arena apilados sobre una mesa.
IdentiÞque fuerzas internas en el sistema y fuerzas externas actuando sobre el sistema.
Pregunta 8: Considere un sistema formado por dos partículas distantes en 1 cm. Indique
un eje con respecto al cual su momento de inercia es nulo.
Pregunta 11: Un canasto con manzanas posa sobre una mesa. En el contexto del sistema
formado por el canasto y las manzanas, indique dos fuerzas ex- ternas y dos fuerzas
internas.
Pregunta 12: Describa dos ejemplos de cuerpos cuyos centros de masa se ubican en el
exterior de ellos. Justifique su respuesta.
1. (2pt.) Determine a que distancia del vértice (donde se dobla la barra) se encuentra
el centro de masa del objeto.
2. (2pt.) Verifique que el resultado tiene sentido si α = 0, π/2 y π
3. (2pt.) Un frabricante entrega barras de L = 1m con un ángulo α no muy preciso
alrededor de π/2. En el control de calidad se puede medir la posición del centro de
masas y de ahí obtener la calidad del ángulo. ¿Con qué precisión se debe medir la
posición del centro de masa para poder determinar si ángulo está entre π/2 − 0,01
y π/2 + 0,01?
En la respuesta indique simplemente si se debe medir con una precisión de 1m,
10cm, 1 cm, 1mm ó 0,1mm. Justifique.
4A.1. Introducción
En este capítulo, continuaremos el desarrallo de las herramientas necesarias para la descrip-
ción de la mecánica de los sistemas extendidos. Como vimos en el capítulo anterior, además
del movimiento del CM los objetos pueden moverse de diversas otras formas. Incluso si un
cuerpo es rígido, este puede rotar sin modificar la posición de su CM. Estudiaremos la leyes
de la estática, que nos permitiran evaluar cuando un cuerpo permanecera en reposo. Las
leyes de la estática son uno de los conocimientos empíricas más antiguos de nuestra civili-
zación. Este a permitido la construcción de monumentos formidables que aún, con todos los
recursos tecnológicos disponibles, nos sorprenden. El descubrimiento básico fué que era posi-
ble amplificar, mediante formas adecuadas, el efecto de una fuerza. De hecho, se le atribuye
a Arquímedes haber dicho “Dame un punto de apoyo y levantaré el Mundo”.
Para que el centro de masas no se mueva se exige que la fuerza neta (suma de las
fuerzas externas) sea nula.
Para que el objeto no rote se exige que el torque neto (suma de los externos) sea nulo.
66
Estática Sistemas Newtonianos 67
La fuerza necesaria para cerrar una puerta es menor si se hace lejos de las bisagras.
Por el contrario mientras la fuerza se hace más cerca del borde fijo de la puerta, más
cuesta cerrarla.
brazo
fza
(−)
11 (+)
00
00
11
00
11
A 1
0
0
1
11
00
00
11
00
11
B 1
0
0
1
11
00
00
11
00
11
(−) 1
0
0
1
0
1
00
11 0
1 00
11 0
1 00
11
Como vimos anteriormente, el uso de palancas para levantar cuerpos es un hecho empírico
en el cual, mediante una combinación de fuerza, puntos de apoyo y lugares donde se aplica la
fuerza, podemos levantar cuerpos que mediante levantamiento directo sería imposible. Este
hecho es explicado mediante las ecuaciones de Newton, para lo cual requeriremos de una
nueva construcción vectorial que llamaremos producto cruz.
Sean A ~yB ~ dos vectores, entonces el elemento
~ =A
C ~×B
~,
es un vector cuya
~ y B;
dirección es perpendicular a ambos, A ~
~ y B;
tamaño AB sin(θAB ), con θAB el ángulo entre los vectores A ~ y
B
A
~×A
1. A ~ = 0, para todo A;
~
~×B
2. A ~ = −B
~ × A;
~
~ × (λB)
3. A ~ = λ(A
~ × B);
~ con λ un escalar;
~τ = ~r × F~ ,
donde ~r es el vector qu une el orige O con el punto P . Físicamente, este producto representa
la habilidad de la fuerza F~ de inducir un giro del cuerpo en torno al origen O. Este origen
es totalmente arbitrario; sin embargo, una vez es escogido, ha de mantenerse durante el
desarrollo de las ecuaciones en cuestión.
P F P
F
r
r
O P F
111111111
000000000
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
111111111
000000000
Σ ri x w i RCM x W
donde hemos identificado R~ CM , el vector centro de masas del sólido con respecto a O. Este
resultado nos permitirá una notable simplicidad en el estudio de sólidos.
Para que el sistema esté en equilibrio estático es necesario que se cumplan las siguientes
condiciones:
III.- ~v = ~0;
~ P = ~0.
IV.- ω
Donde hemos indicado explícitamente el punto P con respecto al cual se determina los tor-
ques y la velocidad angular.
OBSERVACIONES:
1. La suma de fuerzas es vectorial, de modo que muchas veces resulta conveniente ex-
presarla en términos de componentes ortogonales, usualmente según ejes ’x’, ’y’, ’z’
escogidos adecuadamente.
2. Los torques son calculados con respecto a un punto P arbitrario. Recomendamos tomar
torque con respecto a un punto donde se simplifique de mejor forma el cálculo.
3. En este curso nos centraremos en problemas en los cuales todas las fuerzas están
contenidas en un plano. Si el punto P está en el mismo plano, entonces los torques
serán perpendiculares a éste. Con ello, basta con especificar si el torque
Uno puede convenir el signo (+) del torque para cualquiera de estos casos, quedando
los torques en sentido opuesto con signo (–).
θA θB
A M,L
SOLUCION:
• Las fuerzas sobre (tablón+carga): tensiones T~A y T~B , peso carga m~g y peso tablón M~g .
θA y θB
+
x
T
B
T x
A
A
L/2
Mg mg
• De lo anterior tenemos:
vale decir, si los ∠’s son iguales, la carga debe ubicarse simétricamente en el tablón.
Pregunta 3: Una viga de largo L y masa M esta firmemente empotrada a una pared vertical
sobre la superficie de la Tierra. Dibujar la viga, un sistema de referencia, y listar las
componentes de la reacción de la pared sobre la viga en dicho sistema.
Pregunta 5: Una viga de densidad uniforme tiene masa M y largo L. Si la viga se apoya
en un punto fijo al piso a una distancia L/4 de un extremo; ¿Qué fuerza es necesario
aplicar y donde para que la viga este horizontal en equilibrio estático?
Pregunta 7: Una lámina sólida de densidad constante se encuentra en reposo sobre una
mesa. Si la forma de la lámina es la indicada en la figura:
Pregunta 8: ¿Cuáles son las condiciones para que un cuerpo esté en equilibrio?
P2. [5 puntos] Una barra de masa M y largo L puede girar libremente en torno a
su punto de apoyo O en una mesa horizontal en un laboratorio terrestre. La masa se
mantiene en equilibrio estático con una masa m y una cuerda como indica la figura.
Encuentre el ángulo α de equilibrio si m/M = 0,5.
Ejercicio 2: P1. [1 punto] Se quiere soltar una tuerca (nut) muy apretada y para ello
dispone de una llave francesa (o punta-corona, wrench) y una barra (rod). Ordene
las siguientes configuraciones desde la más eficiente a la menos eficiente para soltar la
tuerca (menor a mayor fuerza necesaria)
P2.[5 puntos] Una barra de largo total 3m y masa 10kg se dobla en un ángulo π/4rad
a un metro de uno de sus extremos y se cuelga de una cuerda como muestra la figura.
Encuentre el valor de la tensión T y la masa M si el sistema se encuentra en equilibrio
estático.
Ejercicio 3: P1.[1 punto] Dos masas idénticas están unidas por una barra rígida sin masa.
Si se aplica una fuerza F durante un pequeño intervalo de tiempo ∆t en dos posiciones
distintas como se indica en la figura. ¿En cuál caso el centro de masa del sistema
adquiere una mayor rapidez?
en caso (a);
en caso (b);
es igual;
depende del momento de inercia.
P2. [5 puntos] Un vaso cilíndrico homogéneo (abierto por arriba) de radio basal a y
altura b posa sobre un plano inclinado en un ángulo β c/r a la horizontal y no resbala
debido a un pequeño tope (de altura despreciable) fijo sobre el plano. Suponga que
el vaso se puede modelar como compuesto de una lámina sin espesor y de densidad
superficial (masa por unidad de superficie) constante.
Ejercicio 4: P1. [1 punto] Una masa de 1kg se ata (con una cuerda ideal sin masa) al
extremo de una barra de 1m de longitud. Cuál es la masa de la barra si el sistema se
encuentra en equilibrio al aplicar una fuerza de soporte en la marca de 0,25m como
indica la figura?
P2. [5 puntos] Una esfera homogénea de radio R y masa M se corta en dos mitades.
El sistema se dispone con las mitades cara a cara y la superficie de corte vertical.
A fin de evitar que las semi-esferas se separen, se dispone una cuerda ideal con dos
masas m idénticas como indica la figura. No hay roce entre la cuerda y las semi-esferas.
Determinar las masas mínimas a atar en los extremos de la cuerda de modo que las
semi-esferas no se separen.
Wneto = Kf − Ki (4B.1)
Consideremos que sobre una partícula solo actúa una fuerza conservativa F~c . Entonces,
combinando (4B.1) y (4B.2) obtenemos que ∆K = Wc = −∆U , de donde:
∆(K + U ) = 0 (4B.3)
78
Energía de Rotación Sistemas Newtonianos 79
Ejemplo: Un objeto se suelta desde una altura H, ¿A qué velocidad impacta contra el
suelo?
Solución: Despreciando el roce con el aire, podemos aplicar la conservación de energía
mecánica entre el instante inicial (objeto se suelta desde el reposo) y final (objeto impacta
el suelo):
1 1
mv02 + mgH = mvf2 + mg0
2 2
de donde obtenemos directamente
1
vf = (2gH) 2
Consideremos el caso simple en que la única fuerza externa es el peso. Entonces, para la
partícula n-ésima (n de un total de N ) se puede escribir la ecuación (4B.4)
X
Wnk (F~nc ) = ∆(Kn + Ugn ) (4B.5)
k
donde Wnk (F~nc ) es el trabajo que le realiza la partícula k a la partícula n (de la fuerza
interna), Kn es la energía cinética de la partícula n y Ugn es la energía potencial gravitatoria
de la partícula n. Esta expresión se ve compleja, pero si se suma para todos los valores de
n, se obtiene que se cancelan los trabajos de las fuerzas internas, obteniéndose
!
X X
∆ Kn + Ugn (4B.6)
n n
Si el cuerpo esta rotando con respecto a un eje fijo, todas la partículas tienen igual velocidad
angular, ω , y la velocidad (lineal) de la partícula n-esima es ρn ω, donde ρn es el radio de
la orbita que describe la partícula en torno al eje de rotación. Entonces, la energía cinética
del sistema resulta ser:
X X 1
1X 1
2
K= Kn = m n vn = (mn ρ2n )ω 2 = Iω 2 (4B.7)
2 2 2
La cantidad I se denomina Momento de Inercia, y es una propiedad del cuerpo que juega
un papel
P de gran importancia en la dinámica de las rotaciones. Consideremos su definición
I = (mn ρ2n ). Es claro que I crece con la masa total del cuerpo, pero depende críticamente
de su forma y del eje de rotación. A igual masa, dos cuerpos pueden tener muy diferentes
momentos de inercia en virtud de cuan cercana o lejana este la masa del eje de rotación.
La energía potencial gravitatoria de cada partícula del sólido es Ugn = gmn yn . La la expre-
sión (4B.6) aparece la suma de estas energías que, de acuerdo a la definición del centro de
masa, es X X
Ugn = gmn yn = M YCM g (4B.8)
n n
donde YCM es la altura del centro de masa. Luego, la energía potencial gravitatoria de un
sólido es igual a la de una partícula puntual de masa igual a la del sólido y ubicada a la
altura del centro de masa del mismo.
4B.2.3. Resumen
∆(K + Ug ) = 0 (4B.9)
con K = Iω 2 /2 y Ug = M YCM g.
Figura 11: Momento de Inercia de una lámina, una esfera llena y un cascarón esférico
como la suma de dos momentos de inercia: uno de la barra al lado derecho del centro y la
otra al izquierdo, ambas de masa M/2 y largo L/2.
1 M L2 1
Ic = Iizq + Ider = 2Io = 2[ ] = M L2
3 2 2 12
Propuesto: Calcule el momento de inercia (I) de una barra en torno a un punto ubicado a
1
4 de uno de sus extremos.
Ejemplo: un sistema formado por una barra de largo L y masa M distribuída en forma
uniforme, unida a una esfera pequeña de masa 2M en uno de sus extremos. En un cierto
instante la barra esta dispuesta en forma horizontal, pero luego rotará en torno a uno de
sus extremos hasta quedar en posición vertical. Pongamos nuestro sistema de referencia en
el punto de rotación.
~ ~
~ = 2M Resf era + M Rbarra = − 2M L + M L/2 ß̂ = − 5 Lß̂
Ri
3M 3M 6
Notar que hemos empleado el hecho del que el centro de masa de una barra homogénea está
en su centro (Unidad 3). El signo − y el vector unitario ß̂ en la posición del CM aparecen
por la posición específica del sistema en este instante. Más importante, sabemos que el CM
se ubica a 65 del largo desde el extremo donde no está la esfera.
¿Cuál es la posición del CM cuando el sistema está vertical? No necesitamos calcular la
posición del CM nuevamente. Por simple inspección podemos escribir: R~f = − 65 Læ̂
Cuando un cuerpo tiene algún tipo de simetría, es relativamente sencillo conocer el momento
de inercia c/r a un eje que pasa por el centro de masa. El Teorema de Steiner nos permite
encontrar fácilmente el momento de inercia c/r a cualquier eje paralelo al anterior, desplazado
una distancia R. En su demostración empleamos la geometría y notación de la figura adjunta:
P ~ + r~i , luego
Por definición IO = mi x2 , y ICM = Σmi r2 . Expresamos el vector x~i = R
i i
~ + r~i )2 = R
x~i 2 = (R ~ 2 + 2R~ · r~i + r~i 2
X
IO = mi [R~ 2 + 2R~ · r~i + r~i 2 ]
luego X
IO = ~ 2 + r~i 2 ] = M R2 + ICM
mi [R
Ejemplo: Consideremos una barra delgada en forma de L, con dos partes iguales cada una
de masa M y largo L. Calculemos el momento de inercia respecto su extremo B. Recuerde
que el momento de inercia de una barra de largo L y masa M con respecto a su CM es
1 2
12 M L .
Aunque se trata de un cuerpo sólido, haremos el cálculo descomponiendo la L en sus dos
partes (vertical y horizontal): IB = IV + IH . En cada caso aplicamos Steiner:
1 1 1
IH = ICM + M (L/2)2 = M L2 + M L2 = M L2
12 4 3
1 5 4
IV = ICM + M D 2 = M L2 + M L2 = M L2
12 4 3
luego
5
IB = M L2
3
Podemos resolver este problema fácilmente considerando la energía mecánica del sistema
E = K + Ug . Para Ug tomemos como referencia el nivel donde la barra esta horizontal.
Como la rotación ocurre en torno a un extremo, I0 = 31 M L2 . Suponiendo que la barra
es uniforme, su CM está en su centro geométrico. Si la barra parte del reposo se tiene
Ei = Ki + Ugi = 0 + 0 y a energía final es
1 11 1
Ef = I0 ωf2 + M gyf = M L2 ωf2 + M g(− L) = 0
2 23 2
expresión de la que se obtiene la velocidad angular final y la velocidad final del centro de
masa
ωf = (3g/L)1/2
1 1
vCM f = Lωf = (3gL)1/2
2 2
Figura 14: Dos cuerpos de igual masa pero distinta forma, ¿cuál cae más rápido?
I⊥ = Ixx0 + Iyy0 .
Ixx’
I yy’
??
xx’
1111111111
0000000000
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
?? (A) yy’
(B)
Antes de demostrar este teorema examinemos sus ventajas. Consideremos primero un aro
de masa M y radio R (recuadro (A) de más arriba). Sabemos que el momento de inercia
con respecto a un eje que pasa axialmente por su centro es
I⊥ = M R2 .
Si consideramos dos ejes perpendiculares, en el plano del aro, que se cruzan mutuamente
en el centro, entonces los momentos de inercia con respecto a cada uno de estos ejes será
el mismo (por simetría). Llamemos a cada uno de ellos Io . Aplicando el teorema anterior
tenemos
M R2 = I⊥ = Io + Io ⇒ Io = M R2 /2 .
2
MR
2
MR /2
ML 2/3 ML 2/3
M, L
M, L
Usando el teorema de los ejes perpendiculares podemos calcular el momento de inercia con
de la placa rectangular con respecto a un eje perpendicular a su plano, que pasa por una de
sus esquinas. Escribimos (ver recuadro (B)) de la figura de la página anterior),
M a2 M b2 M 2
I⊥ = Ixx0 + Iyy0 = + → I⊥ = (a + b2 ) .
3 3 3
000000000000
111111111111
x’
000000000000
111111111111
000000000000
111111111111
y’
z 111111111111
000000000000
000000000000
111111111111
000000000000
111111111111
000000000000
111111111111
000000000000
111111111111
000000000000
111111111111
y
000000000000
111111111111
000000000000
111111111111
i
x y
000000000000
111111111111
0000000000000
1111111111111
ρ
i 0000000000000
1111111111111
y’ x’
y 0000000000000
1111111111111
0000000000000
1111111111111
x (xi,y i) 0000000000000
1111111111111
i
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
xi
0000000000000
1111111111111
0000000000000
1111111111111
x y
En el estado A, K = 0 y Ug = M ga/2.
En el estado B: K = (1/2)Iω 2 y Ug = −M g/2.
El momento de inercia en este caso resulta I = 2M a2 /3, con lo cual,
ω 2 = 3g/a .
Lectura Suplementaria
Fisica. Serway Tomo 1. McGraw Hill
Secciones 10.4 y 10.5
Pregunta 2: En relación a la sesión práctica de esta semana, ¿qué rol desempeña el mo-
mento de inercia en las ecuaciones utilizadas para el estudio de la ‘T’ ?
Ejercicio 3: Una masa 2M se ata a una cuerda que esta enrollada en el borde de un disco
de radio R y masa M . El disco puede girar sin roce respecto a un eje que pasa por su
centro. La masa posa sobre un plano inclinado si roce entre ambos.
La masa se suelta del reposo y comienza a bajar sobre el plano inclinado por acción de
la gravedad. La cuerda se va desenrollando del disco, haciéndolo girar de manera que
en todo momento vmasa = Rω, donde vmasa es la rapidez con que desciende la masa y
ω la velocidad angular del disco.
m~r¨ = F~ (4C.1)
~ = m~r × ~v
L (4C.2)
pero ~r˙ = ~v de lo que resulta que el primer término es m~v ×~v , que es nulo por las propiedades
del producto cruz. Además, el segundo término puede ser escrito como ~r × m~a que, usando
la ley de Newton queda
L~˙ = ~r × F~ (4C.4)
Esta última es precisamente la definición de torque ~τ = ~r × F~ , obteniéndose la llamada
ecuación de torque para una partícula.
~˙ = ~τ
L (4C.5)
94
Torque y Momento Angular Sistemas Newtonianos 95
4C.1.1. Ejemplo
Para ver cómo se usa esta ecuación consideremos un péndulo simple de masa m que cuelga
de una cuerda ideal de largo R.
T
m
mg
~ =A
C ~×B
~
es un vector cuya
~ y B;
dirección es perpendicular a ambos, A ~
~ y B;
tamaño AB sin(θAB ), con θAB el ángulo entre los vectores A ~ y
La partícula está a una distancia R del origen y el módulo de la velocidad es v = Rθ̇ por
ser un movimiento circunferencial. Además, como la velocidad es perpendicular al vector
posición se tiene sin θ = 1. Por último la dirección de L~ es perpendicular a la posición y
velocidad, es decir, sale del plano del papel hacia afuera. Si definimos el vector k̂ como el
que es perpendicular al plano se tiene
~ = mR2 θ̇k̂
L (4C.6)
~˙ = mR2 θ̈k̂
L (4C.7)
2
= mR αk̂ (4C.8)
Propuesto:
Considere que el mismo péndulo experimenta además una fuerza de roce viscoso F~v = −γ~v .
Determine la ecuación de movimiento usando el método del torque.
En lo que sigue vamos a considerar que el sólido tiene un movimiento plano. Es decir, la
velocidad del sólido está en un plano (movimiento bidimensional) y la rotación ocurre con
respecto a un eje perpendicular al plano de movimiento. Además, en esta unidad conside-
raremos sólo el movimiento de rotación en torno a un punto fijo. En la próxima unidad se
verá el caso general.
V
ω
Figura 15: Ejemplos de movimientos planos. En el ejemplo de la izquierda el sólido se traslada
y el eje de rotación es perpendicular a la traslación. En el ejemplo de la derecha el movimiento
es sólo de rotación en torno al punto fijo.
En el caso de tener un sólido rígido que gira en torno a un punto fijo O (tal como los descritos
en la Unidad 3 y Unidad 4B), se puede descomponer el sólido en N partículas individuales
y luego se hace tender N → ∞. Como es un sólido rígido, la distancia de cada una de estas
partícula al punto fijo es constante ρi , describiendo un movimiento circular de ese radio en
torno al punto fijo.
111111111111
000000000000
O
ρi vi
ω
De esta forma, el momento angular total del sólido que rota en torno a un punto fijo es
X
~O =
L mi~ri × ~vi (4C.16)
i
X
= mi ρ2i ω k̂ (4C.17)
i
!
X
= mi ρ2i ω k̂ (4C.18)
i
~ O = IO ω
L ~ (4C.19)
donde se ha puesto el subíndice “O” para indicar explícitamente que se mide el momento
angular respecto al punto fijo. Al pasar de la tercera a la cuarta línea se identificó el momento
de inercia respecto al punto fijo O. Además, se definió la velocidad angular vectorial ω ~ = ω k̂
como el vector que tiene la maginitud ω = θ̇ y cuya dirección está dada por el eje de giro y
sentido por la regla de la mano derecha.
Se obtiene entonece que el momento angular total de un sólido es proporcional a su momento
de inercia, que es una propiedad intrínseca del cuerpo, y a la velocidad angular que mide
el estado de rotación en cada instante. Esta relación es análoga a la del mometum lineal
p~ = m~v , donde m es una propiedad intrínseca del cuerpo y ~v mide el estado de traslación
en cada instante.
~˙ = d
L (IO ω
~) (4C.23)
dt
~˙
= IO ω (4C.24)
= IO θ̈k̂ (4C.25)
donde se ha usado que el momento de inercia es una propiedad del sólido y si éste es rígido,
entonces es constante en el tiempo.
Esta ecuación es análoga a la ecuación de Newton m~a = F~total donde IO juega el rol de la
masa asociada al movimiento de rotación. La aceleración angular (cambio en la velocidad
angular) es producida por los torques sobre el cuerpo. Dado un torque fijo, un cuerpo de
mayor momento de inercia tendrá una menor aceleración angular. Es decir, el momento
de inercia indica la dificultad para cambiar (acelerar o frenar) el estado de rotación de un
cuerpo.
Propuestos:
Se tienen dos discos de masa M y radio R, uno de ellos con la masa distribuida uniformemente
y el otro con la masa sólo en la circunferencia exterior. Si se aplica el mismo torque τ sobre
ambos discos, ¿cuál acelererá más?
Para los mismos discos anteriores, si ambos están girando con velocidad angular ω. ¿A cuál
se le debe aplicar un mayor torque para que ambos se frenen en el mismo tiempo?
El torque total tiene una componente interna y otra externa. Calculemos primero la com-
ponente interna: XX
int
~τtotal = ~ri × f~ik (4C.31)
i k
donde para pasar de la primera a la segunda línea se usó el principio de acción y reacción,
f~ki = −f~ik . Como (4C.31) y (4C.33) son válidas, el torque interno se puede escribir también
como el promedio de las dos expresiones
!
1 XX XX
int
~τtotal = ~ri × f~ik − ~rk × f~ik (4C.34)
2
i k i k
1 XX
= (~ri − ~rk ) × f~ik (4C.35)
2
i k
El vector ~rik ≡ ~ri − ~rk es paralelo a la línea que une los centros de las partículas. Por otro
lado, se sabe que las fuerzas fundamentales de la naturaleza (electromagnetismo, gravitación,
fuerzas atómicas, fuerzas nucleares,...) cumplen con la propiedad que la dirección de la
fuerza es paralela a la línea que une los centros de las partículas (por ejemplo, la fuerza de
gravitación universal es Fik = −Gm1 m2 r̂ik /rik 2 ). Debido a esa propiedad, los productos cruz
donde se deja explícita la indicación que el brazo de los torques se mide respecto al punto
fijo O.
4C.5. Resumen
En resumen, la ecuación de torque para un sólido rígido respecto a un punto fijo O se escribe
de las siguientes maneras
~˙ = ~τOext
IO ω (4C.38)
IO θ̈k̂ = ~τOext (4C.39)
4C.6. Ejemplos
Se llama péndulo físico al caso de un sólido rígido de forma arbitraria que puede girar
libremente respecto a un punto fijo bajo la acción de la gravedad.
θ
G
mg
En la figura se representa el péndulo físico que está sujeto del punto fijo O. El centro de
masa del cuerpo está indicado por una G. Las fuerzas externas que actúan sobre el cuerpo
son su peso m~g y la fuerza de sujeción en el punto fijo R.~ Esta fuerza, también llamada
de reacción es la que impide que el cuerpo caiga (por eso tiene una componente vertical),
como también impide que el punto O se mueva hacia los lados (por eso también tiene una
componente horizontal). De esta forma, la fuerza de reacción tiene una magnitud y dirección
en principio desconocidas que se deben determinar de las ecuaciones de movimiento.
La fuerza de reacción actúa en el punto O. Luego, su torque con respecto a O es nulo pues
el brazo es nulo. Aquí nuevamente se ve la utilidad del método de torques porque permite
describir el movimiento eliminando algunas fuerzas desconocidas.
El torque de las fuerzas externas se reduce entonces al torque del peso que, como se vio en
la Unidad 4A, actúa sobre el centro de masa del cuerpo.
donde RG es la distancia del centro de masa al punto fijo y el signo se obtuvo de la regla de
la mano derecha.
La ecuación de movimiento es entonces
ecuación que es muy parecida a la del péndulo simple, pero ahora depende de la forma del
cuerpo a través del momento de inercia.
Propuesto:
Encuentre la ecuación de movimiento de un péndulo físico formado por una barra de largo
L y masa M , que cuelga de un extremo.
¿La aceleración angular es mayor o menor que la un péndulo simple? ¿Por qué?
Considere una polea de masa M , radio R y momento de inercia I que está sujeta por su
centro en O. De la polea cuelga una cuerda ideal en cuyo extremo está sujeta una masa m.
Se busca saber cómo gira la polea.
111111111111
000000000000
ω
M, R, I
O
Sobre la polea actúan tres fuerzas externas: su peso M~g , la reacción del soporte en O R ~ y
~
la tensión de la cuerda T . Las dos primeras fuerzas no ejercen torque respecto al punto fijo
O pues su brazo es nulo. La tensión se aplica en el extremo del círculo, a una distancia R
del centro. Como el vector que va del punto O al punto de aplicación es perpendicular a la
tensión, el brazo es simplemente R. Por último el sentido del torque es, de acuerdo a la regla
de la mano derecha, un vector que sale del plano hacia afuera (según k̂). Luego
~τOext = RT k̂ (4C.45)
Se debe notar que si la cuerda hubiera estado sujeta en el borde derecho del círculo, el torque
habría sido −RT k̂.
Luego, la ecuación de movimiento de la polea es
IO α = RT (4C.46)
con α = θ̈.
El movimiento de la masa m se determina de la ecuación de Newton. Tomando el eje y
vertical hacia arriba se tiene
ma = T − mg (4C.47)
Se tiene, además, la relación cinemática que relaciona lo que baja la masa con lo que gira la
polea
a = −Rα (4C.48)
donde el sigo “-” aparece porque cuando la polea gira en su sentido positivo, la masa baja (es
decir, se mueve en sentido negativo). ¿Cómo se escribiría la relación cinemática si la masa
colgara del lado derecho?
Reemplazando (4C.48) en (4C.47) se obtiene
T = mg − mRα (4C.49)
a = −Rα (4C.51)
g
= − (4C.52)
1 + IO /M R2
lo que indica que acelera más lento que g. Nuevamente, esto es producto de la inercia de la
polea.
Propuesto:
Determine la aceleración angular de 111111111111
000000000000
una polea que tiene dos radios don-
de se enrollan las cuerdas y de la que ω
R1 M, I
cuelgan dos masas. Indique la condi-
R2
ción crítica para que la polea gire en O
uno u otro sentido.
m1
. m2
Tipler y Mosca, Física para la Ciencia y la Ingeniería, Secciones 9.3, 9.4, 9.5, 10.1
(opcionales 10.2 y 10.3)
Pregunta 3: Para el movimiento de un sólido rígido con punto fijo, indique cuál es la
relación análoga a F = m · a
1. Determine la magnitud de todas las fuerzas externas que actúan sobre la esfera.
2. Calcule cuánto tiempo tarda en detenerse la esfera debido al roce.
45 45
Ejercicio 3: Una rueda de masa M , radio R y momento de inercia I está en contacto con
una superficie horizontal rugosa, caracterizada por un coefienete de roce estático µe y
dinámico µd . El centro de la rueda está unido a una cuerda ideal cuyo otro extremo
está unido a una pared fija a la misma altura del centro de la rueda.
Ejercicio 4: Considere el sistema de la figura, compuesto por dos masas unidas por una
barra rígida de largo 2L y masa M . La barra puede girar libremente en torno al
punto medio de ellaen O. En el extremo P de la barra se encuentra una masa m y
en el extremo Q una masa 2m. Todo el sistema se mueve bajo la acción de la gravedad.
θ̈ = f (θ, θ̇)
4D.1. Introducción
En una primera caracterización del movimiento de cuerpos sólidos, en relación a sus contactos
con elementos externos, podemos identificar dos casos de interés particular:
Ciertamente estos no son los únicos casos concevibles. Podemos agregar movimientos com-
puestos de rotación y traslación (como un boomerang en el aire), y de rotación con resbala-
miento (como el neumatico de un vehículo al partir resbalando).
Si pensamos específicamente en una rueda, en la figura de más abajo ilustramos una rueda
dentada que gira, sin deslizar, sobre una superficie igualmente dentada. Cuanto los dientes P
y Q entran en contacto no hay deslizamiento mutuo, de modo que la velocidad instantánea
relativa entre ellos es nula. Puesto que el engranaje inferior está en reposo, entonces Q
necesariamente está en reposo instantáneo. Esta idea se extiende a la rueda que no resbala
de la derecha, donde observamos que el punto Q de la rueda en contacto con el plano está
en reposo instantáneo.
to o
mien ient
Q sbala lam
sin r
e
Q esba
sin r
P P
Entre P y Q la velocidad relativa es casi nula Entre P y Q la velocidad relativa es nula
109
Rodadura Sistemas Newtonianos 110
τo = Io α .
R
R
T T’ Mg
En la figura de más arriba se ilustran dos casos de rotación en torno a un eje fijo: el de la
máquina de Atwood (izquierda) y el de una rueda que puede girar (a modo de péndulo) en
torno al eje P. Ambas se rigen por la misma ecuación del movimiento τf ijo = If ijo α, donde
el eje fijo se ubica en O y P, respectivamente.
De forma muy general, la ecuación de torques (4D.1) es válida cuando ellos se evalúan
con respecto a un eje instantáneamente en reposo (eje instantáneo de rotación), siendo el
momento de inercia evaluado con respecto a ese eje. Por lo tanto podemos aplicarla a una
rueda que rota sin resbalar, donde el punto de contacto está instantáneamente en reposo.
Así,
~C
dL
~τC = = IC α
~. (4D.2)
dt
Estudiemos el caso de una rueda de masa M y radio R rodando sobre un plano inclinado con
roce. No hay resbalamiento (gracias al roce), de modo que el punto de la rueda en contacto
con el piso está instantáneamente en reposo. Denotamos ese punto por ’C’. Denominando
óbjetoá la rueda, entonces las fuerzas externas actuando sobre ella son:
mg
N
f C
x
Podemos escribir la ecuación para el movimiento traslacional del centro de masas (CM),
F~ = M~a, donde ~a representa la aceleración del centro de masas de la rueda. Para este
ejemplo específico,
M~g + f~ + N~ = M~a ,
M g sin β − f = M ax (4D.3)
−M g cos β + N = 0 → N = M g cos β , (4D.4)
Ahora, si aplicamos la ecuación de torques con respecto al punto de contacto, es fácil verificar
τC (M~g ) + τC (f~) + τC (N
~ ) = IC α → M gR sin β = IC α
pues los torques de la fuerza normal y de la fuerza de roce se anulan pues sus brazos son
nulos.
En el caso de una rueda que rota sin resbalar sobre una superficie, el desplazamiento de
su centro se relaciona de forma muy simple con la rotación angular que experimenta. Para
fijar ideas, en la figura de más abajo se ilustra una rueda de radio R en contacto con una
superficie rectilínea. Al desplazar el centro de la rueda en δx, esta rota angularmente en δθ.
La huella impresa sobre la superficie coincide con el arco Rδθ. Así,
δx = Rδθ
R δθ
δ s=R δθ
v = Rω .
a = R ω̇ = R α
Si queremos obtener ax , de la Ec. (4D.7) se tiene que α = M gR sin β/IC , que sustituida en la Ec.
(4D.8) nos da para ax
M R2
ax = g sin β .
IC
Recordar que IC es el momento de inercia con respecto a un eje que pasa por el punto de contacto.
Si denotamos por ICM al momento de inercia que pasa por el centro de masas de la rueda, entonces
el teorema de Steiner nos permite afirmar que IC = M R2 + ICM . Con lo anterior,
g sin β
ax = .
1 + ICM /M R2
Un par de alcances destacables sobre este resultado
Notar que la aceleración del centro de masas es constante, de modo que su velocidad
asociada se relaciona con el desplazamiento mediante
v 2 − vo2 = 2ax ∆x .
Notese que si ICM = 0, se obtiene el resultado conocido para un cuerpo resbalando sin
roce: ax = g sin β. ¿Es esto razonable?.
La ecuación (4D.5) para el roce f sugire una aceleración límite (¿máxima o mínima?)
que garantice no resbalamiento. Se deja propuesto determinar el ángulo β máximo que
garantice que la rueda no resbala.
Una esfera maciza tiene un momento de inercia 2M R2 /5 con respecto a un eje que pasa
por su centro. Verifique que en tal caso
1111
0000
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
H
ESTADO A EVOLUCION A−>B
ESTADO B
KB = KA + M gH (4D.9)
∆x
H
β
De esta última relación inferimos por simple inspección que la aceleración lineal del centro
de la rueda es
g sin β
ax = ,
1 + ICM /M R2
coincidente con la obtenida mediante la aplicación directa de las leyes de Newton. Era un
resultado exigible, dado que la relación energía-trabajo se obtiene de las mismas leyes.
~τCM = ICM α
~.
Fi
ri
P i r = R + ρi
Fijo ρ i
R i
V CM
ω,α
V2 ~ .
~a = − R̂ + α
~ ×R (4D.10)
R
~ +ρ
6. Sustituyendo ~ri = R ~i , y reagrupando términos,
~1 × F~1 + ρ
ρ ~2 × F~2 + · · · + R
~ × (F~1 + F~2 + · · · ) = IP α
~
7. Claramente ρ ~1 × F~1 +~
ρ2 × F~2 +· · · = ~τCM , el torque neto con respecto a CM. Además de
lo anterior, reconocemos F~1 + F~2 +· · · = M~a, con ~a la aceleración del CM. Sustituyendo
en (4D.11) obtenemos
~τCM + R ~ × (M~a) = IP α~.
~ × R̂ = 0, se tiene
8. Al reemplazar ~a de la Ec. (4D.10), usando que R
~ × (M α
~τCM + R ~ = IP α
~ × R) ~.
~ × (M α
9. Se pide que verifique la identidad R ~ = M R2 α
~ × R) ~ , con la cual
~τCM = (IP − M R2 ) α
~,
| {z }
ICM
vale decir,
~τCM = ICM α
~.
4D.1.6. Apéndice
Una aplicación MatLab interesante es dibujar la trayectoria de los puntos de una rueda
cuando este rueda sin resbalar.
x=X+rho*cos(theta)
y=R−rho*sin(theta)
X = Rθ , Y =R.
Esta construcción garantiza que para θ = 0, el centro de la rueda se ubica en (0, R). Las
coordenadas (x, y) de un punto P en la rueda está dada simplemente por
x = X + ρ cos θ (4D.12)
y = R − ρ sin θ . (4D.13)
con 0 ≤ h ≤ 1. Juegue con la siguiente rutina MatLab, modificándola a gusto, para visualizar
lo que resulta. Al hacer h = 0 se está identificando el centro de la rueda, cuya trayectoria
debiera ser recta. Al hacer h = 1 se identifica un punto en la periferia de la rueda. Cualquier
valor intermedio identificará puntos al interior de la rueda.
Referencias:
La rodadura perfecta de un cuerpo rígido se presentan en
Pregunta 3: Escriba -sin resolver- la ecuaciones dinámicas para el disco (masa M y mo-
mento de inercia I) que se muestra en la figura. El cable esta enrollado en un carrete de
radio r (r < R) y se mantiene paralelo al plano inclincado. El disco rueda sin resbalar.
Pregunta 4: Dibuje el DCL para el disco (radio R y masa M) que se muestra en la figura.
Pregunta 5: Dibuje el DCL y escriba Ðsin resolverÐ la ecuaciones dinámicas para el ci-
lindro (masa M y momento de inercia I) que se muestra en la figura. El cable tira al
cilindro desde su centro.
Ejercicio 2: Un disco de radio R y masa M rueda sin resbalar sobre una superficie hori-
zontal rugosa, tirada hacia la derecha por una cuerda ideal que se mantiene paralela
al plano. La tensión de la cuerda es T (constante). La cuerda se va desenrollando sin
resbalar de un carrete de radio r concéntrico al disco (ver figura).
Haga un gráfico de la fuerza de roce (incluyendo signos y magnitud) como función del
radio r.
Ejercicio 3: Un disco de radio R y masa M puede rodar sin resbalar sobre una superficie
horizontal rugosa. El centro del disco esta enganchado a un resorte de constante elástica
k. Inicialmente, el resorte se encuentra en su largo natural y súbitamente el centro de
masa del disco adquiere una velocidad v0 tal que el resorte se comienza a comprimir.
Calcule la comprensión máxima del resorte y compare con el valor que
obtendría si en vez del disco se tratara de una partícula de masa M sin
roce con la superficie.
Indicación: Resuelva este problema mediante trabajo/energía
Ejercicio 4: Una partícula de masa m puede deslizar sin roce sobre un plano inclinado
(ángulo α con respecto a la horizontal). La partícula esta unida al extremo de una
cuerda ideal que se mantiene paralela al plano inclinado. El otro extremo de la cuerda
se enrolla en el borde de un disco de masa M y radio R. La cuerda se va desenrollando
sin resbalar del borde del disco.
Si la partícula se suelta del reposo, determine su velocidad una vez que ha
descendido una distancia D medida a lo largo del disco.
Indicación: Resuelva este problema mediante trabajo/energía
Ejercicio 5: Dos cilindros de radio R y masa M pueden rodar sin resbalar por un plano
inclinado en un ángulo α. Los cilindros están unidos por su centro mediante una cuerda
ideal que se mantiene tensa. El cilindro que está delante tiene su masa distribuida
uniformemente en su volumen, mientras que el cilindro que está detrás tiene su masa
concentrada en el perímetro de radio R.
Si la posición del cuerpo cambia en el tiempo, entonces podemos expresar un pequeño cambio
en la posición como
~ = ds ẑ ≈ dr ẑ
dφ
r r
122
Oscilaciones Sistemas Newtonianos 123
donde hemos utilizado la regla de la mano derecha para expresar el ángulo de forma vectorial;
ds es la longitud del arco de circumferencia entre ambas posiciones, considerado positivo si
el ángulo aumenta. Naturalmente dφ es adimensional (cuociente entre cantidades con las
mismas unidades), para mayor claridad se define el radián como una dimensión angular.
Si el movimiento circular es uniforme, entonces por definición la tasa de cambio del ángulo
θ en el tiempo es constante y se denomina velocidad angular:
~
dφ
ω
~ = = constante, luego
dt
φ(t) = φ0 + ωt (5A.1)
Notamos que las coordenadas x(t) e y(t) del vector posición oscilan entre ±R y las podemos
escribir como
x(t) = R cos(φ0 + ωt) ; y(t) = R sin(φ0 + ωt).
R representa entonces la amplitud del movimiento en x o en y. Al ángulo φ(t) lo llamamos
fase, y al ángulo inicial φ0 lo llamamos constante de fase. El movimiento es periódico, repi-
tiéndose cuando la fase cambia en 2π. El período del movimiento lo obtenemos de imponer
ωt = 2π:
2π
T = (5A.2)
ω
En el caso de movimiento circular, ω representa la velocidad angular, que podemos relacionar
con la velocidad tangencial ~v = ω
~ × ~r. Sin embargo, es muy común poder expresar movi-
mientos periódicos a través de funciones trigonométricas. En estos casos el ángulo φ puede
no representar la posición de un cuerpo en una circumferencia, pero siempre representa la
fase del movimiento, como veremos a continuación.
donde ~x0 es la posición de equilibrio, que corresponde al largo natural del resorte en el
caso de no haber otras fuerzas presentes, y k es una constante de proporcionalidad positiva
(constante del resorte, unidades N/m).
Eligiendo el origen del sistema de coordenadas en el largo natural del resorte y orientando
el eje x a lo largo del resorte podemos escribir la segunda ley de Newton como
F = −kx = ma
o bien
k
ẍ + x=0 (5A.3)
m
donde los puntos sobre la x representan derivadas c/r al tiempo. Para encontrar las soluciones
a esta ecuación, suponemos soluciones del tipo exponencial, es decir x = aebt . Al reemplazar
en al ecuación anterior obtenemos una ecuación algebraica para la constante b:
k
b2 + =0
m
cuyas
p soluciones son imaginarias. Para simplificar la notación, definimos la constante ω0 =
k/m de modo que podemos expresar
b± = ±iω0
Vemos que efectivamente funciones exponenciales son soluciones de la ecuación (5A.3). Sin
embargo, estas son exponenciales de exponente imaginario:
donde a± son constantes y en la última igualdad hemos utilizado la expresión para la función
exponencial compleja. La solución general de la ecuación (5A.3) es una combinación lineal de
las dos soluciones x± . Esta solución incluye dos constantes de integración (a± ) que dependen
de las condiciones iniciales del problema. Esta solución general se puede expresar de una
variedad de formas, siempre con dos constantes de integración, dejamos al lector verificar
las siguientes expresiones:
Donde la última expresión es la que utilizaremos frecuentemente. En este caso las constantes
de integración son la amplitud A y la constante de fase θ0 . La cantidad ω0 = 2π/T la
llamamos frecuencia angular; T es el período del movimiento. Las unidades de la frecuencia
angular son [rad/s]. Notamos que ω0 tiene unidades de frecuencia, pero NO corresponde a
una velocidad angular como en el caso de movimiento circunferencial uniforme, ya que no
hay ningún ángulo en la definición del problema.
En resumen, el movimiento armónico simple se puede caracterizar a través de
Esta es una ecuación muy sencilla, pero al no ser polinomial puede tener un comportamiento
no intuitivo para los alumnos. Les recomendamos graficar la función (5A.8) jugando con los
parámetros para explorar el efecto que tienen las variaciones de cada parámetro en la función
x(t).
Para encontrar una solución para la posición definida por el ángulo θ(t) lo más sencillo es
utilizar la ecuación de torque c/r al punto de apoyo del hilo ya que la tensión no realiza
Esta ecuación es aparentemente muy similar a la ec. (5A.3) sin embargo el que aparezca la
función sin θ en vez del ángulo θ puede introducir diferencias sustanciales en sus soluciones.
Afortunadamente para el caso en que el péndulo oscila con ángulos pequeños, podemos
aproximar la función sin θ ≈ θ si θ << 1 rad. Notamos también que las constantes g y R
son positivas, luego podemos definir la constante
g
ω02 = (5A.13)
R
obteniendo la siguiente ecuación para describir el movimiento del péndulo
θ̈ + ω02 θ = 0 (5A.14)
donde es evidente ahora la similitud con la ecuación (5A.3) que describe el movimiento de
una masa pegada a un resorte ideal. Podemos adivinar inmediatamente la solución de esta
ecuación:
θ(t) = θmax cos(ω0 t + θ0 ) (5A.15)
donde la amplitud A = θmax es un ángulo. Es muy importante distinguir en este caso
la constante ω0 = Rg que representa la frecuencia angular del movimiento de la velocidad
angular ω = θ̇ que para empezar no es constante y tiene un significado físico distinto. El
período del movimiento es s
2π R
T = = 2π . (5A.16)
ω0 g
Un péndulo físico es cualquier péndulo real que no pueda ser bien aproximado por una
masa puntual atada a un hilo ideal. Es decir, cualquier caso en que el momento de inercia
del sistema c/r al punto de apoyo difiera de mR2 . Para distinguirlo del péndulo simple,
utilizaremos M para la masa, d para la posición del centro de masa c/r al punto de apoyo,
I para el momento de inercia del sistema c/r al punto de apoyo. El torque del peso c/r al
punto de apoyo es ahora
τg = −dM g sin θ = Iα = I θ̈
En general notamos que el período de un péndulo físico es mayor que el período de un péndulo
simple con la misma masa concentrada en su centro de masa y menor que el período de un
péndulo simple de la misma longitud.
Notamos que la solución general para cualquier amplitud del movimiento de un péndulo
(simple o físico) no representa un movimiento armónico simple (M.A.S.) que se reproduce
bien por la ecuación (5A.8). Sin embargo el movimiento si es periódico y tiene período y
frecuencia angular bien definidos. Su solución general, sin embargo es matemáticamente
mucho más compleja.
Lectura Suplementaria
El capítulo 13 “Movimiento Armónico Simple"de los apuntes de Massmann son ligeramente
avanzados para este curso. Contiene abundantes ejercicios, algunos de ellos resueltos. Las
secciones 13.1 y 13.2 son apropiadas como material docente de esta semana.
El capítulo 15 del libro de Serway y capítulo 14 del Tipler contienen la materia necesaria
para la unidad 5 (A, B, y C).
t[s]
Pregunta 2: Una masa de m=0.3 kg se ata a un resorte de constante natural k=60 N/m.
Desde la posición de equilibrio se estira el resorte en 0.2 m y se suelta:
Pregunta 3: Determinar el largo que debe tener un péndulo simple sobre la superficie
terrestre de modo que su período sea de 60 s para pequeñas oscilaciones en torno a la
vertical local.
Un hamster se coloca en una rueda circular que puede girar libremente (sin roce) en torno a
un pivote que pasa por su centro P. Para comodidad del hamster se coloca una plataforma
horizontal (bajo el pivote). El sistema se libera con el hamster inicialmente en reposo en
uno de los extremos de la plataforma. El pobre animal, ante la inminencia de su caída
comienza a correr hacia el centro de la plataforma. Tras un poco de práctica el ingenioso
roedor descubre una forma de movimiento que mantiene la plataforma (y la rueda que la
sostiene) estacionaria y horizontal. Determine dicha forma de movimiento.
5A.5.2. Solución
Dado que el hamster debe comenzar a moverse desde el reposo, necesariamente la plataforma
ejerce una fuerza que le permita acelerar. A su vez, el hamster ejerce la misma fuerza sobre
la plataforma. Dicha fuerza es naturalmente −ma, donde m es la masa del hamster y a su
aceleración. La posición del hamster, tomando como origen el centro de la rueda es (x, δ),
donde δ es la distancia entre el pivote y el centro de la plataforma horizontal. La fuerza
neta en dicho punto es (−ma, mg), de modo que el torque neto es: ~r × F~ = m(xg + δa)ẑ.
Igualando a cero el torque obtenemos a = − gδ x, i.e. la ecuación de un oscilador armónico.
La solución debe cumplir las condiciones iniciales x(0) = −L/2 (donde L es el largo de la
plataforma horizontal) y v(0) = 0. De este modo obtenemos:
r
L g
x(t) = − cos t . (5A.18)
2 δ
1. Determine la ecuación de movimiento para las desviaciones del sistema desde su posi-
ción de equilibrio vertical.
1. (1 punto) Un objeto descanza en equilibrio estático desde un resorte colgado del te-
cho del laboratorio. Al bajar el objeto una pequeña distancia, la suma de su energía
potencial elástica y gravitacional:
a) Permanece constante,
b) Aumenta,
c) Disminuye.
Ejercicio 2:
a) Igual,
b) Mayor,
c) Menor.
133
Osc. Amortiguadas Sistemas Newtonianos 134
átomos, ellos se llevarán energía y el sistema a los ojos del obervador macroscópico parecera
perder energía. En la medida en que nos interesemos sólo por el movimiento macroscópico
podemos hablar de una fuerza neta que hace el aire sobre el globo. En general podemos
representar la fuerza del aire sobre el globo de la forma
F~ = −f (|~v |)v̂ ,
con −v̂ representando un vector unitario en sentido opuesto al desplazamiento. Con esta
expresión se asume muy poco sobre los mecanismos que llevan a la forma de f . La forma
general de f (v) depende del regimen de velocidades, densidad del aire, propiedades termo-
dinámicas, geometría del cuerpo, entre muchas otras. En la práctica el verdadero sentido de
la ecuación anterior se completa al medir f mediante experimentos. Esta forma de teoría se
conoce como fenomenológica. La forma anterior posee una gran ventaja, que el estudiante
puede verificar con gran facilidad, para f > 0 (pero, de otro modo, completamente arbita-
rio) la energía mecánica del sistema es siempre decreciente. Es decir, esta forma de roce es
disipativa.
En general el comportamiento de la fuerza de roce viscoso debe ser determinada en un
experimento6 . Aunque es un comportamiento complicado y no universal, se reconocen, sin
embargo, dos regímenes de velocidades donde f (v) toma formas particularmente simples:
f (v) → fs (v) ∝ v ,
f (v) → fr (v) ∝ v 2 ,
11
00 11
00 dx/dt
00
11 00
11
00
11
00
11 00
11
00
11
(t=0) F=− b dx/dt
6
En la actualidad existen diversas formas de determinar la fuerza usando una simulación computacional.
vx (t) = A e−t/τ .
Aquí, A representa una constante que está restringida por la condición inicial, vale decir,
información sobre la velocidad en t = 0. Recordamos que en t = 0, vx = vo , entonces A = vo .
De esta forma,
vx (t) = vo e−t/τ .
Notar que la velocidad se atenua exponencialmente (1/e ≈ 0.368), lo que se ilustra en la
tabla siguiente:
t vx
0 1.000 vo
τ 0.368 vo
2τ 0.135 vo
3τ 0.050 vo
4τ 0.020 vo
5τ 0.007 vo
El resultado anterior para vx (t) puede ser utilizado para obtener x(t). Escribimos esta vez
dx
vx (t) = vo e−t/τ ⇒ = vo e−t/τ .
dt
Esta vez nos preguntamos por aquella función que al derivar resulta una exponencial en t.
La respuesta nuevamente es una exponencial. Planteamos
con C una constante que ha de ser determinada por la condición inicial. Si exigimos que en
t = 0 la posición del globo coincide con el origen (x = 0), entonces C = vo /τ , con lo cual
x(t) = vo τ (1 − e−t/τ ) .
Un gráfico de esta función se ilustra en la figura siguiente, y se observa que para un tiempo
muy grande, x → vo τ . El globo no sobrepasará esa distancia. Cuan lejos queda tal punto
dependerá de τ = m/b; mientras más chico sea b (la fricción), más lejos llegará el globo. Lo
mismo ocurre si es más masivo.
1.0
Velocidad [ vo ]
0.8
0.6 v(t)
0.4
0.2
0.0
0 1 2 3 4 5 6
alcance
1.0
Posicion [ vo τ ]
0.8
0.6
x(t)
0.4
0.2
0.0
0 1 2 3 4 5 6
Tiempo [ τ ]
La caída vertical por gravedad cuando actúa el roce viscoso se puede representar mediante
la ecuación ÿ = g − (1/τ )ẏ. Si el objeto parte del reposo en y = 0, entonces y(t) queda dado
por
2 t −t/τ
y(t) = gτ +e −1 .
τ
La velocidad terminal Vt es aquella que adquiere el cuerpo cuando deja de acelerar. En este
caso es cuando el peso (mg) equipara la fuerza por roce (bVt ). Así,
mg
Vt = = gτ .
b
FStokes = 6πRη v
donde v es la rapidez del fluído. Nótese que la fuerza es proporcional al radio de la esfera, y
no a su área transversal.
1 1 − e−βgt
vy (t) =
β 1 + e−βgt
La cantidad r
1 c
= ,
β m
representa la velocidad terminal (Vt ), de modo que podemos re-
escribir convenientemente
1 − e−gt/Vt
vy (t) = Vt .
1 + e−gt/Vt
O x(t)
0.8 -1
τ=2 s; Ω=1 s
0.6
0.4
nodos
0.2
x [m]
-0.2
-0.4
-t/2τ
-0.6
x=Ae cos(Ωt)
-0.8
-1
0 1 2 3 4 5 6 7 8 9 10
Tiempo [ s ]
La descomposición vectorial conveniente es, en este caso, mediante las coordenadas polares.
Para un movimiento circunferencial de radio L, la velocidad y aceleración quedan expresadas
por
~v = Lφ̇ φ̂ ~a = Lφ̈ φ̂ − Lφ̇2 r̂ .
L
(!bL! )
! (T)
m !
(mg) r
Una persona lanza verticalmente un proyectil hacia arriba, con rapidez v0 . Despreciando la
resistencia del aire, el tiempo que demora en subir hasta su altura máxima es igual al tiempo
que demora en caer desde ese punto hasta su posición original (ambos son iguales a v0 /g).
Considerando la resistencia del aire: £cual de los dos tramos demora más?.
5B.8.2. Solución
Consideremos una altura dada. El proyectil pasa por dicha altura dos veces una de ida y otra
de vuelta. Como en ambas circunstancias la energía potencial es la misma, la diferencia en
energía cinética deben igualar a la diferencia entre energías totales. Dicha energía es siempre
mayor en el tramo de subida. Esto es fácil de ver, pues en el intervalo entre que sube y baja
el roce viscoso solo pudo restar energía. La energía cinética es menor en el tramo de bajada.
Esto es cierto para cada altura en el recorrido total, se sigue que demora más en bajar que
en subir. En este problema no hemos hecho supuestos sobre la forma del roce viscoso, solo
que disipa energía.
FS = 6πRη v
Pregunta 4: Indique las precauciones que hay que tener durante el manejo del globo du-
rante la primera práctica.
Pregunta 8: La velocidad terminal de un objeto al caer en un medio viscoso está dada por
vt = mg/b. ¿A qué corresponde b y qué tipo de fuerza de roce se trata?
1. Estime la fuerza de arrastre del agua sobre el nadador cuando se mueve hori-
zontalmente a 2 m/s. Suponga Cd = 0, 7 y el área transversal del orden de 0,1
m2 .
2. Estime la rapidez con que debe mover las manos, con respecto al agua detenida,
a fin de compensar la fuerza de arrastre obtenida en la primera parte. Haga
una buena estimación de la sección transversal de las manos, y use en este caso
Cd =1,0.
Para velocidades grandes domina la primera de ellas (Fa ), mientras que para pequeñas
velocidades prima la de Stokes. Para R= 1 cm, determine la rapidez a la cual estas fuer-
zas son comparables en el caso del aire (η ∼0,017 mPa s=0,017 milipascal×segundo),
y agua (η ∼1 mPa s). Estime en cada caso la rapidez mínima para que el movimiento
resulte amortiguado dominantemente por viscosidad.
Ejercicio 6: ¿En qué unidades queda representado el cuociente η/ρ? A esta cantidad se le
denomina viscosidad cinemática y es simbolizada por ν. Calcule ν del agua y el aire.
5C.1. Introducción
En las unidades anteriores se estudiaron oscilaciones armónicas mecánicas. Se consideró
también el efecto del roce viscoso (de algún tipo) sobre estas oscilaciones, lo cual se llama un
oscilador amortiguado. Este oscilador puede ser representado gráficamente como se muestra
en la figura 1, el cual consta de tres elementos: una masa (inercia), una fuerza restitutiva
(resorte) y el roce (amortiguador). Es importante notar que no se ha especificado aún el
origen del roce en este sistema. Puede ser un roce superficial con el plano en la cual se apoya
la masa, o puede ser roce viscoso con el aire o algún líquido lubricante. Notemos también
que el largo natural del resorte simplemente nos dice que en equilibrio, es decir cuando la
suma de fuerzas es nula, la masa se enuentra naturalmente a una distancia lo de la pared,
lo que define x = 0. En el caso de un péndulo formado por un globo, la fuerza restitutiva
está dada por la gravedad, y el roce está dominado por la fuerza de arrastre de Rayleigh.
En esta unidad estudiaremos el caso de un oscilador amortiguado forzado. Esto significa que
a la imagen anterior agregaremos un forzaje con una dependencia explicita en el tiempo. La
idea es que si uno fuerza el sistema a una frecuencia bien particular, el movimiento oscila-
torio será amplificado. Esta frecuencia es denominada frecuencia de resonancia o frecuencia
natural de vibración. Un primer ejemplo es un columpio. Se trata de un péndulo, en el cual
se realiza un forzaje periódico en el tiempo moviendo las piernas de cierta manera hacia
adelante y hacia atrás. Aquí, se deben mover las piernas sincronizadamente con el péndulo,
si uno lo hace más rapido, o más lento, el columpio no funciona.
Por simplicidad, supondremos en esta unidad que el amortiguamiento está dado por un roce
viscoso lineal. Esto ayuda en el tratamiento matemático y numérico del problema, aunque
puede ser poco realista en algunas situaciones. El énfasis estará dado en el efecto del forzaje
más que en considerar bien cuál es el mecanismo dominante de roce en estos sistemas.
145
Osc. Forzadas Sistemas Newtonianos 146
k, lo
M
0 x(t)
Figura 18: Esquema de un oscilador amortiguado compuesto por un resorte, una masa y
algún tipo de roce. El forzaje externo puede ser incluido como una fuerza F~ (t) que actúa
sobre la masa M .
Figura 19: Derrumbe del Hotel Continental ocurrido durante el terremoto de 1985 en cuidad
de México. Este edificio entró en resonancia con la frecuencia de forzaje del temblor, y sus
pisos superiores colapsaron. Otros edificios contiguos, más altos o más bajos, no sufrieron
los mismos daños que este hotel.
en la figura 1, es el más simple que puede considerarse. Más adelante, ustedes irán agregando
grados de complejidad a los sistemas bajo estudio. Aquí mencionamos algunos ejemplos
muy generales de sistemas oscilantes forzados que presentan las caracteristicas básicas que
discutiremos a continuación. En principio cualquier sistema que se encuentra en una posición
de equilibrio estable realizará oscilaciones en torno a dicho punto en caso de ser perturbado.
Es interesante notar que, de este modo, el comportamiento dictado por el ejemplo en la
figura 1 es un comportamiento universal. Por esto queremos decir que muchos sistemas
desarrollan oscilaciones análogas a las déscritas en la figura 1, independientemente de su
naturaleza (sistemas mecánicos, electrónicos, electromagnéticos, geológicos, atmosféricos,
etc). El simple modelo de la figura 1 sirve para caracterizar diversos sistemas.
Los instrumentos musicales son buenos ejemplos de osciladores forzados. Sin embargo,
pertenecen a una clase más compleja pues son medios contínuos y por lo tanto pre-
sentan un conjunto de frecuencias de resonancia (notas musicales de una cuerda de
guitarra por ejemplo). Se estudiarán los modos de resonancia de una cuerda un poco
más adelante en el semestre.
Los edificios tienen frecuencias naturales de vibración, por muy compleja que sea su
estructura. Si un terremoto ocurre a una frecuencia de resonancia del edificio, este
puede llegar a romperse. Por supuesto que técnicas actuales reducen esta posibilidad.
Un ejemplo de un colapso real se muestra en la figura 2.
Existe la leyenda de que durante el siglo XIX un ejército Francés produjo el colapso
de un puente al cruzarlo debido al forzaje resonante que produjieron los soldados al
marchar. Se dice que en la actualidad de advierte a los soldados de no marchar sobre
un puente.
con 2
2 1
Ω = ωo2 − . (5C.4)
2τ
Nótese que esta solución es válida si Ω2 > 0. Cuando el roce es pequeño, τ es grande y por
lo tanto Ω ≈ ωo . Las constantes A y φo se determinan de las condiciones iniciales x(0) y
ẋ(0).
En el caso de un oscilador amortiguado forzado, con una fuerza F (t) = Fo sin(ωt), la ecuación
de Newton resulta
M ẍ = −kx − bẋ + Fo sin(ωt). (5C.5)
0.05 0.1
(a) (b)
0.04 0.08
0.03 0.06
0.02 0.04
posición (m)
posición (m)
0.01 0.02
0 0
!0.01 !0.02
!0.02 !0.04
!0.03 !0.06
!0.04 !0.08
!0.05 !0.1
0 10 20 30 40 50 0 10 20 30 40 50
Tiempo (s) Tiempo (s)
Figura 20: Ejemplos de soluciones x(t) obtenidas con M = 0,778 kg, τ = 4,72 s, k = 78
N/m, Fo = 0,25 N, x(0) = 0,05 m, ẋ(0) =p0. En (a), ω = 7 rad/s y en (b) ω = 9,8 rad/s. La
frecuencia natural de resonancia es ωo = k/M = 10 rad/s.
con
ω
tan δ = . (5C.8)
τ (ωo2 − ω2)
Al final de esta guía se indican los pasos a seguir para demostrar que esta es efectivamente
la solución de la ecuación (5C.6). Las constantes A y φo se determinan de las condiciones
iniciales x(0) y ẋ(0). El primer término de esta solución corresponde a un transiente, es
decir después de un cierto tiempo este término decae a cero, tal como en la solución de un
oscilador amortiguado sin forzaje externo. Es el segundo término, proporcional a Fo , el que
da la solución estacionaria, es decir aquella que perdurará en el tiempo. De esta manera se
puede interpretar δ como una diferencia de fase entre la solución estacionaria y el forzaje.
Cuando b → 0, se tiene τ → ∞, por lo que e−t/2τ → 1. Además, en este límite, la amplitud
de la parte estacionaria B(ω) → Fo /(M (ωo2 − ω 2 )), y por lo tanto presenta una divergencia
en ω = ωo .
Ejemplos de la solución x(t) para dos valores de la frecuencia de forzamiento se muestran
en la figura 3, lejos (a) y cerca (b) de la frecuencia de resonancia natural ωo . Es importante
0.16 3.5
(a) (b)
0.14 3
0.12
2.5
0.1
2
" (rad)
B (m)
0.08
1.5
0.06
1
0.04
0.5
0.02
0
5 10 15 5 10 15
! (rad/s) ! (rad/s)
Figura 21: (a) Curva de resonancia de la amplitud B como función de ω. (b) Diferencia de
fase δ como función de ω. Ambas curvas fueron obtenidas con M = 0,778
p kg, τ = 4,72 s,
k = 78 N/m, Fo = 0,25 N. La frecuencia natural de resonancia es ωo = k/M = 10 rad/s.
notar de que a pesar de que el decaimiento τ = 4,72 s, en este caso (por las condiciones
iniciales) el transiente tiene un efecto importante sobre la solución para tiempos más largos,
hasta aproximadamente 40 s.
Si nos fijamos únicamente en la solucion estacionaria (e(−t/2τ ) → 0), ésta tiene una amplitud
dada por
Fo /M
B=q 2 . (5C.9)
(ωo2 − ω 2 )2 + ωτ
Ejemplos de la amplitud B y el desfase δ como funciones de ω se muestran en la figura 4. Se
observa que la función B(ω) tiene un máximo. Imponiendo dB(ω)/dω = 0 se obtiene que el
máximo ocurre cuando 2
2 2 2 1
ω = ωr = ωo − 2 . (5C.10)
2τ
De nuevo, cuando el roce es pequeño, τ es grande y por lo tanto el sistema es resonante, es
decir su respuesta es máxima, cuando ω ≈ ωo .
La interpretación de δ es simple. Para ω ωr , δ ≈ 0, es decir siendo la amplitud de
oscilación pequeña, está en fase con el forzaje. Por el contrario, para ω ωr , también la
respuesta de la amplitud de oscilación es pequeña pero el desfase es δ ≈ 180◦ . Finalmente,
para ω = ωr , se tiene δ = 90◦ .
Figura 22: Distintas posibilidades para el incremento en energía mecánica debido al forza-
miento. La esfera representa la partícula, de masa M que se mueve con velocidad v y es
forzada mediantela fuerza F . En los casos en que P > 0, la rapidez de la masa se incrementa
debido a la fuerza.
independientemente, la potencia neta entregada será nula. En el caso en que existe una
sincronización entre F y v, i.e. que oscilan a la misma frecuencia, la tasa neta de entrega de
enegŠa en cada ciclo será positiva y por lo tanto la energía cinética se incrementa en cada
oscilación, es decir tenemos una resonancia. Como la velocida debe estar sincronizada con
la fuerza, entre la posición y la fuerza hay un desfase de π/2.
Al incluir disipación la esencia del argumento se mantiene, pero ahora hay una perdida
adicional de energŠa debido al roce viscoso. Como esta perdida se incrementa con la rapidez,
Proce = −mv 2 /τ , siempre negativo, la potencia entregada por el forzamiento externo, será
eventualmente compensada porla disipación. Esto nos permite entender porque se llega a un
regimen estacionario, donde la partícula oscila con amplitud finita, incluso bajo condiciones
de resonancia. Este argumento permite predecir el valor de la amplitud final. Si la partícula
oscila con amplitud A, la potencia entregada por la fuerza externa es Pext = ωF A sin2 (ωt),
mientras que la potencia disipada es Proce = −mA2 ω 2 sin2 (ωt)/τ . De este modo, la condición
de compensación es simplemente: A = (F/M )τ /ω, consistente con el resultado analítico.
donde ahora τ = (M + m)/b y ωo2 = k/(M + m). Suponemos ahora que el carro pequeño
tiene algún tipo de motor que permite asegurar que y(t) = yo sin(ωt) para todo t, entonces
ÿ = −yo ω 2 sin(ωt). Con esto, la ecuación de Newton resulta
1 myo ω 2
ẍ + ẋ + ωo2 x = sin(ωt), (5C.15)
τ M +m
que básicamente igual a la ecuación (5C.6) salvo que la constante que va delante de sin(ωt),
llamada Fo , ya no es independiente de ω. La solución analítica (5C.7) no cambia, sólo
debe cambiarse Fo /M por myo ω 2 /(M + m). De igual manera, la amplitud de la solución
estacionaria B(ω) cambia sólo en el reemplazo Fo /M → myo ω 2 /(M + m), es decir si cambia
como función de ω. Explícitimante tenemos
myo ω 2 /(M + m)
B=q 2 . (5C.16)
(ωo2 − ω 2 )2 + ωτ
y(t)
0' m
k, lo
M
0 x(t)
Figura 23: Esquema de un oscilador amortiguado compuesto por un resorte, una masa y
algún tipo de roce. El forzaje externo sobre la masa M se obtiene por un movimiento
oscilatorio impuesto a la masa m, es decir con y(t) = yo sin(ωt).
0.16
(a) (b)
!1
0.14 10
0.12
0.1
B (m)
B (m)
0.08
!2
10
0.06
0.04
0.02
!3
10
5 10 15 5 10 15
! (rad/s) ! (rad/s)
Figura 24: Curvas de resonancias de la amplitud B como función de ω, en escala lineal (a) y
en escala semilogaritmica (b). La curva segmentada es la misma presentada en la figura 4a.
La curva contínua se obtuvo con la pecuación (5C.16) con M = 0,735 kg, m = 43 g, yo = 5,8
cm, τ = 4,72 s, k = 78 N/m, ωo = k/(M + m) = 10 rad/s. Los datos (◦) corresponden a
medidas experimentales.
(3) Se escribe
B1 sin(ωt) + B2 cos(ωt) = B sin(ωt − δ),
y se encuentran B y δ.
Fo /M
x(t) = Ae−t/2τ cos(Ωt + φo ) + q sin(ωt − δ),
ω 2
(ωo2 − ω 2 )2 + τ
Recordamos que las oscilaciones armómicas (amortiguadas o no) ocurren en sistemas que
son pertubados cuando están cerca de un punto de equilibrio. Así, hay una fuerza que tiende
a llevar al sistema de vuelta al equilibrio, pero la inercia (masa) hace que primero le tome
un tiempo llegar al punto de equilibrio pero además que se pase de largo, haciendo que lleve
a cabo un movimiento oscilatorio.
Las ondas aparecen por un mecanismo similar pero a diferencia de la las oscilaciones, las
ondas involucran la oscilación de muchos grados de libertad (en principio un número infinito).
Estas oscilaciones, al estar acopladas se sincronizan de muchas maneras dando origen al
comportamiento ondulatorio. Por ejemplo, ya no es una partícula puntual o un sólido ríguido
156
Ondas Propagativas Sistemas Newtonianos 157
que cuelga de un péndulo, sino que es toda una cuerda, por ejemplo, la que puede ser
deformada (sacada del equilibrio) punto a punto.
Así una cuerda tensa entre sus extremos tiende a estar recta. Si es deformada, cambiándole
la forma, tratará de volver a su forma original, generando un movimiento que llamaremos
ondulatorio.
Cuando uno tiene una cuerda tensa, una membrana de tambor o la superficie del agua se
observa que dependiendo de cómo se deforme el sistema se presentan los siguientes fenóme-
nos:
Oscilación colectiva En una cuerda de guitarra, se observa que toda la cuerda oscila
colectivamente de arriba a abajo con una amplitud que depende de la posición. A este
tipo de movimientos se le llama modos normales y serán vistos en la Unidad 6B.
Vamos a ver que los dos tipos de fenómenos se describen de manera unificada como fenó-
menos ondulatorios.
Vista en perspectiva
Tal como se ve en la figura, las varillas pueden girar en torno a su centro, pero como están
soldadas al hilo, eso tuerce el hilo, lo que provoca un torque a las varillas vecinas.
Rotulemos con i = 1, . . . , N las varillas y sea θi el ángulo que forma la varilla i-ésima con la
horizontal. Calculemos el torque sobre esta varilla. Si la varrila que le se sigue (la i + 1) tiene
el mismo ángulo que ésta, el hilo no está torcido y el torque es nulo; si el la varilla forma un
ángulo mayor con la horizonal (θi+1 > θi ) entonces el hilo está torcido y le hace un torque a
la varilla i-ésima que tiende a hacer que θi crezca para que los dos ángulos sean iguales, es
decir un torque positivo (noten que por acción y reacción el torque que i le hará a i + 1 es
negativo); finalmente si θi+1 < θi el torque será negativo. Considerando que las diferencias
de ángulos no son muy grandes, se puede considerar entonces como una buena aproximación
que el torque sobre i debido a i + 1 será proporcional a la diferencia de ángulos
donde T es una constante que habrá que medir y que depende de las propiedades del hilo.
De igual manera se puede calcular el torque sobre la varilla i debido a la anterior (la i − 1),
resultando τi debidoa i−1 = T (θi−1 − θi )
Sumando los dos resultados se tiene el torque total sobre la varilla i
Por otro lado, sabemos que el movimiento de las varillas está dada por la ley de Newton
para sólidos rígidos con un punto fijo (donde están soldadas al hilo)
d2 θi
I = τi (6A.1)
dt2
d2 θi
I 2 = T (θi+1 − 2θi + θi−1 ) (6A.2)
dt
d2 θ d2 θ
I 2
= T ∆2 2 (6A.5)
dt dx
que se puede escribir como
d2 θ 2
2d θ
= c (6A.6)
dt2 dx2
donde se ha definido por comodidad
p
c= T ∆2 /I (6A.7)
Esta ecuación se llama ecuación de ondas y se debe entender que es para el ángulo θ que
depende tanto de x como de t. Es decir, θ(x, t).
La ecuación (6A.5) se puede entender como una ecuación de Newton, donde a la izquierda
está la inercia y a la derecha las fuerzas. Si la colección de varillas tiene una gran variación
de ángulos (d2 θ/dx2 grande) entonces habrá una fuerza mayor que provocará una mayor
aceleración en cada punto.
6A.2.2. La cuerda
Un segundo ejemplo de ondas en una dimensión se puede obtener del análisis del movimiento
de una cuerda que está atada en sus extremos y que se mantiene tensa, con tensión T . La
cuerda tiene una masa M y largo L, de la que se obtiene una densidad de masa lineal
ρ = M/L. Este ejemplo es ampliamente analizada en los textos, en los que pueden buscar
explicaciones alternativas.
y
x
La cuerda es levemente extensible y se puede deformar verticalmente (se dice que se deforma
transversalemente a la dirección en la que está estirada).
Para describir la dinámica de la cuerda, estudiemos lo que pasa en una vecindad de un punto
x de la cuerda. Llamaremos y(x, t) la deformación vertical de la cuerda en el punto x en el
instante t, de igual manera como llamabamos θ(x, t) a la torsión de la varilla en un punto x
en el instante t.
Consideremos un trozo de cuerda que está entre x − ∆ y x + ∆. 7 . Al hacer el DCL de ese
trozo, las fuerzas que aparecen son las tensiones a ambos lados. Si llamamos θ1 al ángulo
que forma la tangente a la cuerda en x − ∆ y θ2 al ángulo que forma la tangente a la cuerda
en x + ∆, la fuerza total sobre ese trozo es
T
θ2
θ1
x−∆ x x+ ∆
La figura muestra una cuerda muy deformada, pero vamos a considerar el caso en que la
deformación es pequeña. Eso se caracteriza porque vamos a suponer que los ángulos que
forma la tangente a la cuerda con la horizontal son siempre pequeños (θ 1). Si no se
cumple esta hipótesis debemos hacer un análisis mucho más complejo.
Si los ángulos son pequeños, entonces, la fuerza se simplifica pues cos θ ≈ 1 y sin θ ≈ θ ≈
tan θ,
Por otro lado, la ley de Newton, dice que la fuerza es masa por aceleración. La masa de ese
trozo de cuerda es la densidad por el largo que es 2∆. Luego la masa es m = 2∆ρ. Como el
movimiento es puramente vertical, la aceleración es
d2 y
~a = ŷ (6A.17)
dt2
Combinando todo, se encuentra que la ecuación para el trozo de cuerda es
d2 y dy dy
2∆ρ 2 ŷ = T (x + ∆) − (x − ∆) ŷ (6A.18)
dt dx dx
Calculemos las derivadas temporales y espaciales de u. Las derivadas espaciales son directas
y se calculan como
du
= f 0 (x − ct) (6A.25)
dx
d2 u
= f 00 (x − ct) (6A.26)
dx2
donde f 0 y f 00 son la primera y segunda derivada de f respecto a su argumento. Para calcular
las derivadas de u respecto al tiempo hay que usar la regla de la cadena pues se debe derivar
el argumento (x − ct) respecto a t. El resultado es
du
= (−c)f 0 (x − ct) (6A.27)
dt
d2 u
= (−c)2 f 00 (x − ct) = c2 f 00 (x − ct) (6A.28)
dt2
Notamos que si reemplazamos estas derivadas en la ecuación de onda, se produce una can-
celación a ambos lados para cualquier f .
Luego, hemos encontrando una solución de la ecuación de ondas.
u(x, t) = f (x − ct) (6A.29)
De manera análoga se puede mostrar (queda de tarea) que si g(x) es una función cualquiera
entonces
u(x, t) = g(x + ct) (6A.30)
también es solución.
Por último, una combinación cualquiera
u(x, t) = f (x − ct) + g(x + ct) (6A.31)
es solución (verificar!).
El análisis de estas soluciones y su discusión se verá en la clase y en la práctica de laboratorio.
6A.4.1. Problema
1. Determine que el tiempo total que tarda un pulso en propagarse desde el extremo
inferior al superior. (Indicación: Recuerde que la tensión de la cuerda cambia según la
altura).
2. ¿Como cambia su respuesta si una masa M se cuelga desde el extremo inferior de la
cuerda?.
6A.4.2. Solución
.
La velocidad depun pulso depende de la tensión y de la densidad de masa, a través de
la relación c = T /ρ. En este problema la tensión depende de la posición debido a que la
cuerda debe sostener su propio peso. En el punto y medido desde el punto inferior, la tensión
debe sostener un peso ρgy. De este modo tenemos que la velocidad satisface la relación:
v 2 = gy. La velocidad satisface la relaciŮn de un movimiento uniformemente acelerado:
vf2 − vi2 = 2a∆. Con vi = 0, ∆ = y y a = g/2. El tiempo es entonces entregado por la
relación:
1
y(t) = y0 + vi t + at2 = L (6A.32)
2
q
obteniendo la solución: t = 2 Lg La clave para este problema es identiÞcar la dependencia de
la tensión, con respecto a la posición. Una vez hecho eso, el resultado es directo si logramos
identiÞcar las ecuaciones resultantes con las de un movimiento acelerado uniformemente.
Esta identiÞcación no deja de ser extraña pues el valor de la aceleración es igual a la mitad
de g ¡pero hacia arriba!. Esto se debe a mĞs arriba más peso debe soportar la cuerda debido
a su propio peso
En el caso en que una masa cuelga del extremo inferior de la cuerda, el único cambio
corresponde a la tensión T = ρgy + M g. De este modo obtenemos: v 2 = gy + M/ρg.
Haciendo nuevamente la analogŠa con el movimiento uniforme acelerado, obtenemos:
s r
M m
t=2 1+ −1 , (6A.33)
ρg M
donde c es constante y dado por la densidad deqla cuerda y la tensión generada por la masa
externa. En el otro caso límite tenemos, t ∼ 2 Lg . Trivialmente lo calculado en la primera
parte.
d2 u
= c2 ddxu2
2
dt2
d2 u
= c2 ddxu2
2
dt2
(a) Frecuencia
(b) Velocidad de propagación de la onda
(c) Longitud de onda
(d) Período
P3. [4 pts] Un pulso se mueve en dirección x en un sistema de varillas acopladas por torsión
τ , todas las varillas de igual largo L = 0,2 m, masa m = 0,3 kg, separación ∆ = 0,01
m, y de diámetro muy pequeño. El pulso está descrito por
θ(x, t) = A e−(ax+bt) ,
2
168
FI1A2 - SISTEMAS NEWTONIANOS GUIA TEORICA
Departamento de Fı́sica Unidad 6B: Ondas Estacionarias
Facultad de Ciencias Fı́sicas y Matemáticas Profs: H. Arellano, D. Mardones, N. Mujica
Universidad de Chile Semestre 2008-1
1. Introducción
Todos los cuerpos exhiben algún grado de flexibilidad en cuanto pueden experimentar pequeñas
deformaciones, sean estas de tipo longitudinal (a lo largo del cuerpo) o transversal (normales al
cuerpo). Estas perturbaciones, inicialmente forzadas por un agente externo, pueden viajar a través
del medio dando lugar a ondas y pulsos.
Un caso simple, pero muy relevante, es una cuerda tensa dispuesta en forma horizontal. Supon-
dremos el eje x alineado con la cuerda. En este caso la deformación vertical y(x,t) corresponde a
los pequeños cambios de posición vertical de las partı́cula que forman la cuerda. Las partı́culas se
refieren a un elemento infinitesimal de la cuerda entre x y x + dx. Al aplicar la segunda ley de
Newton dicho elemento infinitesimal de la cuerda, obtenemos la ecuación de onda:
∂2y 1 ∂2y
= (1)
∂x2 c2 ∂t2
p
donde c = τ /ρ, con τ tensión de la cuerda y ρ densidad lineal de masa. Notar que c depende
exclusivamente de las propiedades del medio y no de las condiciones iniciales o amplitud de las
deformaciones.
En la clase anterior se demostró que la función
satisface la ecuación de onda. Tal como se discutió en la unidad 6A, las funciones f y g representan
la forma de una onda viajera o pulso que se desplaza a la derecha y a la izquierda, respectivamente,
con una rapidez c (velocidad de fase).
En esta unidad, estudiaremos el caso en que f y g son funciones sinusoidales, las cuales al sumarse,
y bajo condiciones apropiadas, describen ondas estacionarias. Estas ondas parecen no moverse ni a
la derecha o izquierda en el tiempo, pero son el resultado de la suma de ondas que individualmente
si lo hacen (f y g).
2. Ondas Armónicas
1
Otoño 2008 Sistemas Newtonianos – Oscilaciones 2
Notar que en x = nλ la forma sinusoidal se reproduce nuevamente. Por esta razón λ se denomina
longitud de onda.
Cuando t > 0 la onda comienza a avanzar con velocidad de fase c, y supongamos que lo hace hacia
la derecha. Usando (3) en (2) con g = 0 obtenemos entonces:
donde hemos definido T = λ/c. Interpretemos ahora este nuevo parámetro. Por simplicidad consi-
deremos x = 0, de manera que (4) se reduce a y(0, t) = A sin(2πt/T ) cuyo gráfico se muestra en
la figura 2. Claramente la partı́cula en x = 0 (y cualquier otra partı́cula) experimenta oscilaciones
armónicas de amplitud A y perı́odo T . Entonces la ecuación (4) describe una onda armónica, de
amplitud A y con longitud de onda λ y perı́odo T que viaja hacia la derecha.
Una forma simplificada de escribir (4) es:
donde k = 2π/λ es el número de onda y ω = 2π/T es la frecuencia angular [rad/s]. Algunas veces se
emplea la frecuencia f = 1/T (Hz=s−1 ). Las ecuaciones (4) y (5) suponen que y(x = 0, t = 0) = 0
lo cual no siempre es asi. Una versión mas general de (4) se escribe como
donde φ se denomina constante de fase. Tanto la amplitud A como la fase φ dependen de las
condiciones iniciales, es decir de las deformaciones impuestas a la cuerda en t = 0.
Las condiciones iniciales también pueden dictar el valor de λ (y k) con lo cual el valor de T (y ω)
queda completamente definido pues:
p
λ/T = ω/k = c = τ /ρ
Alternativamente, la condición inicial puede dictar el valor de T lo cual fija el valor de λ. Las
relaciones anteriores indican que la longitud de onda (o número de onda) NO es independiente del
perı́odo (o frecuencia) cumpliéndose que:
Ondas largas (número de onda pequeño) son ondas de perı́odo largo (baja frecuencia)
Ondas cortas (número de onda grande) son ondas de perı́odo corto (alta frecuencia)
Hasta ahora hemos considerado que la cuerda tiene un largo infinito: −∞ ≤ x ≤ +∞. Veamos
ahora que pasa cuando la cuerda es finita tal que −∞ ≤ x ≤ 0. La condición de borde en x = 0
puede ser de dos tipos:
Supongamos que un pulso u onda se acerca desde la izquierda hacia el punto de empotramiento y
definamos t = 0 cuando la perturbación alcanza x = 0. Entonces:
¿Qué pasa cuando t > 0? Para responder esta pregunta emplearemos el principio de superposición
y supongamos que la cuerda se extiende hacia el infinito también a la derecha de x = 0. Imaginemos
que en el lado imaginario de la cuerda viaja una perturbación idéntica a la real pero invertida y
moviéndose hacia la izquierda:
La fórmula anterior satisface la condición de empotramiento y predice además que pasa para t > 0:
en este caso la perturbación se refleja en x = 0 (es decir comienza a avanzar hacia la izquierda)
invirtiendo su forma pero manteniendo todos los demás parmetros (Figura 3a).
En este caso, la partı́cula en x = 0 puede cambiar su posición pero la tangente a la cuerda siempre
se mantiene horizontal. Si esta condición no se satisface actuarı́a una fuerza transversal finita lo que
generarı́a aceleraciones infinitamente grandes. La condición de borde en este caso es: ∂y(0, t)/∂x =
0. Por analogı́a con el caso anterior se puede demostrar que esta condición de borde es satisfecha
por:
La ecuación anterior indica que nuevamente una perturbación viajera se refleja en x = 0 pero esta
vez comienza a retroceder manteniendo su forma completamente (Fig. 3b).
Supongamos que se generan ondas armónicas en una cuerda finita empotrada en uno de sus ex-
tremos. De acuerdo a la solución general vista en 3.1, en este caso la solución está dada por
y(x, t) = A sin(kx − ωt) − A sin(kx + ωt). Aplicando los teoremas de trigonometrı́a es fácil demos-
trar que:
y(x, t) = 2A sin(ωt) sin(kx) (7)
La ecuación (6) fue deducida de la superposición (suma) de dos ondas viajeras, pero es una onda
estacionaria! (no aparece el término del tipo x − ct). Esta onda estacionaria tiene número de onda
(y λ) igual a las ondas originales (sin(kx)) y oscila en el tiempo al igual que la original (ω o T no
cambian). Los valores k y ω aun satisfacen ω/k = c. La forma de esta onda se muestra en la figura
(4).
Notemos que en los nodos (separados cada λ/2) la amplitud es 0 mientras que en los antinodos
(también separados por λ/2) la amplitud es el doble de la original (2A).
Supongamos ahora que la cuerda esta empotrada en ambos extremos y tiene largo L. Ya no es
posible tener valores arbitrarios de k (u ω) pues los extremos de la cuerda deben ser nodos: y(x =
0, t) = y(x = L, t) = 0 para todo t. La ecuación (6) cumple la primera condición en forma trivial,
pero la segunda requiere
con n un número natural. Además, como c = λ/T = λf → fn = nc/(2L) Los pares (λn , fn ) definen
los modos normales de la cuerda. A medida que n aumenta, disminuye el largo de la onda y aumenta
su frecuencia (Fig. 5).
λn = 4L/(2n − 1) y fn = nc/(4L)
λn = 2L/n y fn = nc/(2L)
Lectura Suplementaria
Los capı́tulos 15 y 16 del libro de Tipler o los capı́tulos 16 y 18 del Serway contienen la materia
necesaria para las unidades 6A y 6B.
Figura 5: Modos normales en una cuerda con sus dos extremos fijos.
7A.1. Introducción
En esta unidad estudiaremos la fuerza que ejerce un conjunto de partículas cuando rebo-
tan en una superficie. Partiremos con nociones simples acerca de rebotes de partículas y
construiremos un modelo más elaborado para cuantificar el rebote de un chorro de ellas.
Lo interesante de este modelo es que nos permite hacerlo extensivo a diversos sistemas
de nuestro entorno. Por ejemplo, entender la fuerza de un gas sobre las paredes de su
contenedor; la fuerza de arrastre del aire sobre un globo cuando se mueve rápidamente; la
fuerza del viento; la fuerza que ejerce una cadena al caer sobre el suelo; el rol de las aspas
de un helicóptero para que vuele; la fuerza que ejerce el chorro de una mangera; etc.
Conviene introducir algunas nociones que nos serán útiles en la discusión que sigue. La
primera de estas es la densidad de partículas, denotada por n, que corresponde al número
de partículas por unidad de volumen. Así, si n(~r) es el número de partículas por unidad de
volumen en la posición ~r, entonces el número de partículas δN en un pequeño volumen de
tamaño δV que abarca ~r está dado por
δN = n(~r) δV .
δV
n ρ
δ N = n δV
δM = ρ δV
Al igual que el caso anterior, podemos definir la emphdensidad de masa, que denotaremos
por ρ. De esta forma, la masa δM contenida en un pequeño volumen de tamaño δV estará
178
Presión Colisional Sistemas Newtonianos 179
dado por
δN = n(~r) δV .
Ciertamente ρ puede depender de la posición, lo que se representa por ρ = ρ(~r).
Si todas las partículas son de igual masa, mo , entonces
ρ = mo n ,
δ~
p = p~f − p~i ,
donde p~i y p~f son los momenta inicial y final, respectivamente. Si el rebote es elástico, con
la partícula incidiendo con rapidez v, entonces la variación de momentum está dada por
δ~
p = 2mo v cos φ n̂ ,
con φ el ángulo de incidencia con respecto a la normal n̂ del plano del rebote y m la masa
de la partícula.
δ p=2mv cos φ n
φ φ
p p
i f
111111111111111111111111111111111111
000000000000000000000000000000000000
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
000000000000000000000000000000000000
111111111111111111111111111111111111
A
φ φ
v δt cos φ
φ
v
δt
δ V =Av δt cos φ
En este caso cuantificaremos el momentum neto δ P~ que imprime la superficie a las partículas
en un lapso δt. En la figura se ilustra un chorro de partículas que incide sobre la superficie.
En ella se muestran las partículas que colisionarán la superficie entre el instante t y t + δt.
El cambio de momentum será igual al cambio experimentado por una de las partículas,
multiplicado por el número de ellas que colisionarán la pared. Tal número, que denotaremos
δN , está dado por
δN = n × Vol = n × A v δt cos φ .
Entonces, el cambio de momentum experimentado por todas las partículas que golpean (en
el lapso δt) es
Identificamos entonces la fuerza media responzable del cambio de momentum de las partí-
culas rebotando:
F~col = 2ρ A v 2 cos2 φ n̂ .
Denotando por A⊥ = A cos φ, la sección transversal del chorro de partículas, entonces
F~col = 2ρ A⊥ v 2 cos φ n̂ .
el cuadrado de la rapidez;
que expresa la fuerza por unidad de área que ejerce el chorro al incidir las partículas. A
esta cantidad se le denomina presión, y dimensionalmente se expresa como fuerza/área.
El el sistema internacional, una unidad de presión se expresa en N/m2 , correspondiente a 1
pascal:
1 N/m2 = 1 pascal = 1 Pa
La presión es una cantidad física muy útil para caracterizar la manifestación mecánica de
un medio fluído. En nuestro modelo es el resultado del cambio de momentum de partículas
en movimiento. En este modelo, las partículas son monoenergéticas (todas chocan con la
misma energía cinética) y colimadas (todas ellas se mueven en la misma direccón).
En el caso del aire, las partículas que golpean un pequeño elemento de superficie provienen
de todas las direcciones dentro de medio hemisferio. Además, las partículas no son monoener-
géticas sino se distribuyen probabilisticamente. El cálculo de la presión en este caso requiere
de herramientas que van más allá de las que se cuentan a este nivel. El origen microscópico
de la presión es el mismo que hemos introducido en esta unidad.
φ
v
T
F
φ
2L
/3
L/
2
Mg R
φ
P
Si el sistema es estático, entonces al imponer torque neto nulo con respecto a P tenemos
2L L 1 2
LT cos φ − F − M g cos φ = 0 → (T − M g) cos φ = F .
3 2 2 3
Si denotamos por To la tensión que ocurre cuanfo F = 0,
1
To = M g ,
2
entonces
3
F = (T − To ) cos φ .
2
Si aceptamos que la fuerza F → Fcol , entonces
3 8
4ρ A⊥ h g cos φ = (T − To ) cos φ → (T − To ) = (ρ A⊥ h) g .
2 3
Este resultado de forma bastante sencilla es consistente con el hecho de que la rapidez del
impacto en la placa es proporcional a v 2 . Además, la ausencia de la dependencia en el ángulo
φ es por la orientacion vertical del cordel que soporta la placa. Las relaciones algebraicas
permiten la cancelación de cos φ provenientes de Fcol . El objetivo de esta unidad es investigar
la validez de esta relación.
Pregunta 4: ¿Cuál sería el ángulo de rebote de una partícula incidiendo con un ángulo φ
si la energía cinética del rebote es λ veces la inicial?
Pregunta 7: Estime la fuerza que ejerce viento de 100 km/h sobre la pared de una casa de
2,5 m de altura y 5 m de ancho.
Pregunta 10: Estime el tiempo que dura la colisión de un puñado de porotos cuando estos
son soltados desde 50 cm de altura.
Ejercicio 2: Un balde con 5 kg de arena se deja caer sobre una balanza durante 1 minuto.
Determine y grafique el peso que registra la balanza hasta los 90 segundos considerando
que los granos de arena no rebotan.
Ejercicio 3: Para efectos de una estimación sencilla, considere las aspas de un helicóptero
formadas por dos placas circunferenciales de aproximadamente 1 m2 cada una, ubicada
a 5 m del eje. Las placas mantienen un ángulo de inclinación entre 25 y 30 gracos con
la horizontal. Estime el período de rotación de la hélice a fin de que con su rotación
se pueda levantar al helicóptero de 500 kg de masa.
1111
0000
1
0
0
1 0000
1111
Ejercicio 5: Cuando una partícula da un rebote elástico sobre una pared, el cambio de
momentum en la colisión es δp⊥ = 2mo v cos φ, con φ el ángulo de incidencia con
respecto a la normal.
Determine δp⊥ en el caso en que el choque ocurra sobre una superficie perfectamente
resbaladiza pero con pérdida de energía, caracterizada por
Ejercicio 6: En la figura se ilustra una cañería que termina en forma de ‘L’, de donde sale
un chorro de agua en forma vertical. Por la cañería pasa agua (ρ=1000 kg/m3 ), a razón
de 2 litros cada minuto.
7B.1. Introducción
Comenzaremos por reconocer que exiten distintos estados de la materia. Todos sabemos que
el agua puede encontrarse en forma de hielo. Es también usual encontrarse que un charco
o pozo de agua líquida puede evaporarse, aunque sea en forma parcial. Estas observaciones
denotan que el agua puede estar en forma líquida, gaseosa o sólida. Estos tres estados de la
materia son muy generales, cualquier material conocido puede bajo condiciones apropiadas
ser transformado de uno de estos estados a otro.
De manera bastante simple podemos decir que un sólido tiene un volumen y una forma
bien definidos. Un líquido en cambio tiene un volumen definido, pero su forma no lo es.
Finalmente, un gas puede ser considerado como un estado que no tiene ni forma ni volumen
bien definidos, pues ambas características dependerán del recipiente que lo contenga.
Tanto los líquidos como gases corresponden a lo que se conoce como un fluido: se trata
de una colección de átomos o moléculas dispuestas aleatoriamente y unidas por fuerzas
cohesivas débiles y por fuerzas ejercidas por la pared del contenedor. Estos fluidos pueden
ser descritos desde el punto de vista de la Mecánica de Newton. Como hemos visto, esta
Mecánica Newtoniana puede ser extendida a un medio compuesto por un gran número de
constituyentes elementales. En el caso de los fluidos, su descripción mecánica se divide en la
Hidroestática y la Dinámica de fluidos. La primera descripción corresponde, como su nombre
lo dice, a la estática de fluidos, y es el tema central de esta unidad. La segunda descripción
corresponde al estudio de cómo las fuerzas que actúan sobre un fluido afectan su movimiento.
También se puede incluir el cómo el movimiento de un fluido afecta el movimiento de un
objeto, como el caso de un avión. Algo de esto ya ha sido abordado al estudiar la fuerza de
arrastre de Rayleigh en la unidad 5B.
En la unidad anterior se estudió el concepto de presión colisional. Se determinó cómo un
chorro de partículas, sean átomos o granos macroscópicos, ejercen una fuerza sobre una
pared que depende, entre otras cosas, de la velocidad del chorro. Este ejemplo es bastante
ilustrativo, y el cálculo es relativamente simple. En particular se supone que todos los átomos
o granos llegan con la misma velocidad a la pared y por lo tanto transfieren la misma cantidad
187
Arquímedes Sistemas Newtonianos 188
de momentum a ella.
Consideremos ahora que pasa con un fluido en reposo, es decir un fluido en una situación es-
tática. ¿Qué significa exactamente que el fluido esté en reposo? ¿Que cada átomo o molécula
que lo constituye está en reposo, es decir con velocidad nula? Esto no es así. Un fluido está en
reposo cuando al considerar cualquier pequeño volumen, pequeño comparado con el volumen
del contenedor pero grande para aún así contener un gran número de átomos o moléculas,
el promedio de las velocidades de todos los constituyentes elementales es cero. Por lo tanto,
es este pequeño volumen el que está en reposo, pero no necesariamento los constituyentes
elementales. Si cada posible pequeño volumen del fluido cumple esta condición, entonces el
fluido completo está en una condición estática.
7B.2. Presión
Es bien conocido que un fluido en reposo ejerce una fuerza sobre las paredes de su contendor.
Basta pensar en un globo. Al inflarlo, éste se llena de gas. Despues de cerrar el globo, el
gas rápidamente llega a un estado en reposo. Podemos deducir esto simplemente por que no
vemos el globo moverse sólo en una dirección. Ahora, recordemos que el globo se mantiene
inflado estando el gas a su interior en reposo. El material que conforma el globo es muy
elástico, y por lo tanto éste trata de volver a su forma inicial (globo desinflado). Pero no
puede, simplemente porque el gas que está a su interior ejerce una fuerza sobre sus paredes
internas. El globo debe finalmente llegar a un compromiso, entre la fuerza que ejerce el
gas y la fuerza elástica del material que lo compone. La fuerza que ejerce el fluido sobre
una pared, dividido por el área de ésta, es una cantidad física que se llama presión, y que
denotaremos P . Es una cantidad escalar y, para un fluido homogéneo, es una cantidad
constante independiente del espacio.
El hecho que un fluido en reposo ejerce una fuerza sobre las paredes de su contendor implica
que efectivamente los constituyentes elementales del gas no están en reposo individualmente,
pues debe existir una transferencia de momentum entre el fluido y las paredes. Como esta
fuerza es para todos efectos prácticos constante, podemos concluir que los átomos o molé-
culas que componen el fluido tienen una rapidez característica. Notemos que es esta rapidez
característica de los átomos o moléculas de un fluido lo que fija su temperatura, T ∼ hv 2 i.
Un fluido más frío tiene átomos o moléculas que se mueven con una rapidez típica menor, y
viceversa. Además, experimentos muestran que en la mayoría de las situaciones, la presión
es proporcional a la temperatura, P ∼ T . Esto es consistente con la presión colisional estu-
diada en la unidad 7A. La diferencia es que la cantidad v (en la dependencia en v 2 ) ahora
no es la velocidad de un chorro o de un flujo, sino que una rapidez promedio.
En resumen, la presión es una cantidad escalar, definida como
|F~ |
P = ,
A
donde F~ es la fuerza que se ejerce sobre una pared de área A. En el sistema internacional
Figura 26: Representación esquemática de un volumen del fluido. El promedio de las veloci-
dades puede se cero, el fluido está en una situación estática. Pero el promedio de la rapidez,
no es cero. Esto define la temperatura de un fluido, T ∼ hv 2 i.
P = Po + ρf gh.
Esta ecuación nos dice que la diferencia de presión entre dos profundidades es proporcional
a la diferencia de profundidades.
Figura 28: Esquema del funcionamiento de una gata hidráulica. Un pistón de área pequeña
se mueve una distancia ∆x1 gracias a una fuerza F1 . Gracias a este movimiento, un pistón
de área más grande A2 se desplaza una distancia x2 ejerciendo el fluido una fuerza F2 sobre
el pistón.
La magnitud de la fuerza de empuje es igual al peso del fluido desplazado por el objeto.
Luego, si un objeto sumergido desplaza una masa M del líquido, el empuje que siente el
objeto es E = M g. Esta fuerza de empuje siempre es en sentido contrario al del peso del
objeto, y por lo tanto reduce el peso aparente de éste. Atención que esta masa es la del
líquido, no la del objeto. Escribiendo lo mismo de otra manera obtenemos
E = ρf Vo g,
donde ρf es la densidad de fluido, y Vo el volumen del objeto. Si no hay más fuerzas presentes,
la ecuación de Newton resulta
Mo~a = E ŷ − Mo g ŷ
donde Mo es la masa del objeto. La aceleración es ~a = ay ŷ, y su signo depende de la
diferencia de densidades del objeto con el fluido: ay = (ρf − ρo )g/ρo . Así, un objeto mas
denso que el fluido tendrá una aceleración negativa, en el mismo sentido que la gravedad, y
se hundirá. En cambio, un objeto de densidad menor al del fluido, tendrá una aceleración
positiva y flotará.
En el caso que un objeto flote, determinemos el volumen del objeto que esta sumergido con
respecto al volumen total del objeto. La fuerza de empuje será en este caso E = ρf gVs ,
donde Vs es el volumen del objeto que está sumergido. El peso del objeto es Mo g = ρo Vo g.
En equilibrio estas fuerzas son iguales, por lo tanto
Vs ρo
= .
Vo ρf
Este resultado dice que si la densidad de objeto es similar a la del fluido, entonces, gran
parte de su volumen se encontrará sumergido. Un caso a considerar es el de un cubo de hielo
de agua sumergido en agua. Siendo las densidades ρo = 917 kg/m3 y ρf = 1000 kg/m3 ¿cuál
es el porcentaje del cubo que se encuentra sumergido?
Pregunta 3: Determinar la fracción del volumen bajo la superficie del agua de un hielo
que flota en equilibrio. La densidad del hielo es 917 kg/m3 . Justificar.
Pregunta 4: Dos vasos idénticos están llenos de agua hasta el mismo nivel. En uno de los
vasos flotan cubos de hielo. Cuál vaso pesa más? Justificar.
Guía de Prácticas
194
Métodos Numéricos Sistemas Newtonianos 195
B. Materiales: Matlab
C. Experiencias
Experiencia 3.- Usando el programa anterior programe la detección del primer cruce en
φ = 0. Con este programa calcule el periodo del péndulo usando g = 9,8m/s2 , L =
π π π
0,5m, ∆t = 0,01s y los ángulos iniciales φ0 = 10 , 4, 2
Complete la siguiente tabla:
Ángulo inicial φ0 Periodo
D. Conclusiones
Presente de manera concisa las conclusiones objetivas de la sesión en general, no debe resumir
otra vez todos los resultados, sino aquellos más importantes. Señale cómo afecta la precisión
del resultado el uso de diferentes valores de ∆t e indique la tendencia que tiene el periodo
al aumentar el ángulo inicial.
E. Lecturas recomendadas
Aprender a usar el programa Measurement & Automation para verficar que el sensor
de fuerza está funcionando correctamente
B. Materiales
C. Experiencias
Ahora verifique que el sensor de fuerza que está conectado con el canal de entrada
análogo número 0 (AI0) está trabajando adecuadamente. Para ello lance la aplicación
Measurement & Automation (Software de Medida y Automatización), cuyo ícono se
encuentra en el escritorio de su PC.
Al abrirse la Barra del Menú Principal, seleccione Configuration → Devices and Inter-
faces → NI-DAQmx → NI USB-6008: ÒDev nÓ, siendo n un número, normalmente 1.
Seleccione Self-Test, la respuesta debe ser The device has passed the self-test, de otra
forma existe un problema de conectividad. Pida ayuda al profesor o un profesor auxiliar
en este caso.
Si el mensaje aparecido es el correcto, presione OK.
La tarjeta requiere +5 V de alimentación para ello, abra Test Panel y haga clic en
Analog Output. Asegúrese que el sensor está siendo alimentado con 5 V. Para ello
ingrese 5 en la casilla Output Value y presione Update. Atención de no cambiar los
parámetros por defecto que se encuentran en la parte superior de esta ventana.
Ahora pruebe que el sensor de fuerza mide correctamente. Seleccione Analog Input.
Use los siguientes parámetros de adquisición:
Mode: Continuous
Max Input Limit: +10 V
Min Input Limit: -10 V
Input Configuration: RSE
Channel Name: AI0
Rate (Hz): 1000
Samples to read: 1000
Lance la medida presionando sobre el botón Start. Como el modo selccionado de ad-
quisición es contínuo, debería ver una medida constante de aproximadamente 2,5 V,
lo que implica que la fuerza sobre el sensor es nula. Puede presionar con su dedo sobre
el gancho del sensor y verá como la señal varía en tiempo real en su pantalla.
Complete la siguiente tabla con los valores aproximados de voltaje que entrega el sensor
de fuerza en dos posiciones, con el gancho apuntando hacia arriba y, hacia abajo:
Posición Voltaje
Arriba
Abajo
D. Conclusiones
Presente de manera concisa las conclusiones objetivas de la sesión en general, no debe resumir
otra vez todos los resultados, sino aquellos más importantes. Enumere las dos fuentes de error
más importantes en su proceso de medición.
Tabla experiencia 3.
B. Materiales
C. Experiencias
Experiencia 2: Considere una barra uniforme de masa M y longitud L que gira con ve-
locidad angular ω en torno a un eje perpendicular. El eje se ubica en la barra, a una
distancia λL de uno de sus extremos. Por lo visto en los apuntes y por simple analisis
dimensional, se tiene que la energía cinética total de la barra se puede expresar como
1
K = γM L2 ω 2 ,
2
donde γ es un coeficiente adimensional que depende de los valores de λ, vale decir,
γ = γ(λ).
D. Análisis e informe
E. Lecturas recomendadas
A. Objetivos
B. Materiales
El montaje experimental es similar al ilustrado en la figura 29, pero con sólo una masa
colgada. El prototipo permite colocar masas de diferentes calibres, a diferentes distancias,
para ejercer fuerzas y torques de distintas magnitudes. El hilo ejercerá la tensión necesaria
para mantener la barra en equilibrio estático. El sensor de fuerza permite medir dicha tensión
para cada configuración de equilibrio.
C. Experiencias
1. Preliminares
Al inicio de cada experimento es necesario verificar el funcionamiento correcto de los
elementos a usar. Hoy utilizarán el mismo sensor de fuerzas que en la Unidad-2. Este es
probable que esté ya instalado, sin embargo antes de utilizarlo deben revisar que esté
instalado correctamente. Recuerden que los cables del sensor deben estar conectados
a la tarjeta de medición en las siguientes puertas:
T = A ∗ V + B[N ]
errorbar(pos,tension_exp,error_tension,’+’).
Noten que hemos utilizado un simbolo + para indicarle a matlab que ponga dicho
símbolo en cada punto graficado y no los una con una línea. Podríamos haber utilizado
algún otro signo. El comando errorbar grafica una línea vertical en cada punto que
indica el error asociado a dicha medida. Estamos suponiendo que no hay error en las
medidas del eje x; en este caso posición.
Ahora buscamos comparar las medidas de la tensión con el valor esperado a partir de
las ecuaciones de la estática. Dado que conocen la geometría del problema (distancias,
ángulos, masas, y g), definan en Matlab un arreglo de posiciones que cubra el rango de
las medidas previamente. El número de puntos no tiene por qué ser el mismo que para
la medida experimental. Para cada posición calculen la tensión del hilo esperada para
mantener el equilibrio. Esto es muy sencillo y pueden traer el programa de Matlab
listo si lo desean. Ojo que es muy común equivocarse con las unidades al escribir este
tipo de programitas. Graficar la tensión calculada como una línea contínua junto al
gráfico anterior.
Ejemplo de gráfico
errorbar(pos,tension_exp,error_tension,’+’)
hold on
plot(pos_calc,tension_calc)
Adjuntar al informe una copia impresa del gráfico bien rotulado con unidades en ambos
ejes, fecha y grupo.
4. Experiencia 3: Variación de la tensión en función del ángulo [2 puntos].
En esta experiencia estudiaremos la variación de la tensión del hilo en función del
ángulo entre el hilo y la barra. Para ello, dado que la polea es fija, se les pide variar
la posición de enganche del hilo sobre la barra. Con una masa fija a una distancia
también fija del eje, deben medir la tensión para al menos tres posiciones distintas del
anclaje.
Para cada posición, mida también el ángulo entre el hilo y la barra, ¿con que precisión
creen que pueden medir dicho ángulo? Graficar la tensión en función del ángulo, de
nuevo se les pide graficar exclusivamente los puntos medidos con el comando errorbar
de Matlab.
En el mismo gráfico dibujar el resultado esperado teóricamente suponiendo que cono-
cen todos los parámetros experimentales de manera exacta.
5. Conclusiones [1 punto]
4B.1. Introducción
En la Unidad 4B comenzamos a estudiar la dinámica de sólidos rígidos. Esta semana estu-
diaremos la dependencia de la energía cinética de rotación con el momento de inercia y la
diferencia de energía potencial de un cuerpo.
Para ello utilizaremos en efecto un péndulo físico (es decir, no puntual) compuesto por una
barra homogénea de longitud 2L = 0,6m, y una barra del mismo material pero largo L
cruzada a la primera (esta última la denominaremos cruceta). Es decir, tenemos una especie
de T sujeta cerca de su extremo por un eje fijo en torno al cual es sistema podrá girar
libremente. Ojo que el mismo sistema lo volverán a utilizar en tres semanas más.
El objetivo de esta experiencia es estudiar el teorema de conservación de energía. La idea
es soltar la barra desde el reposo a una posición inicial fija, y medir la velocidad angular
y el cambio de altura del centro de masa un intervalo de tiempo después. Luego podrán
contrastar el valor de la energía cinética de rotación experimental obtenida a partir de la
medida de velocidad angular, con la energía cinética de rotación teórica calculada a partir
del teorema de conservación de la energía mecánica.
El montaje experimental les permitirá definir cuatro configuraciones distintas (posiciones de
la cruceta), y luego podrán graficar las cantidades medidas y calculadas en función de la
posición de la cruceta. Pongan atención a la forma funcional de dicha dependencia (teórica)
y a los posibles errores que permitan o no obtener dicha función a partir de las mediciones.
En este documento los guiaremos a través de las mediciones que esperamos que realicen y
su modo de presentación. El informe se entregará en un documento aparte más conciso.
A. Objetivos
B. Materiales
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario verificar el funcionamiento correcto de los
elementos a usar. Hoy utilizarán una cámara web para grabar la caída de la T desde
una posición inicial no-vertical y desde el reposo. Deje caer desde el reposo la T desde
diversas posiciones iniciales: ¿Cuánto tiempo demora la T en pasar por la vertical?
¿Depende de la posición inicial? ¿Cuán repetibles son las condiciones iniciales?
La cámara graba 30 cuadros por segundo, luego el tiempo de integración de cada cuadro
es de aproximadamente 1/30 de segundo. El uso de la cámara web está documentado en
el documento anexo Uso-CamaraWeb.pdf. ¿Cuán nítidas son las imágenes obtenidas?
¿Con qué precisión se puede medir el eje o un borde de la T en las imágenes grabadas?
La velocidad angular de la T la medirán a partir de la comparación de dos cuadros
consecutivos en el video que graben; el objetivo es aproximarse a la velocidad angular
instantánea en el momento de pasar por la vertical, luego la precisión de dicha medida
será mejor si la T cae más lento. Por otro lado, es necesario que el centro de masa de
la T baje una distancia apreciable para poder utilizar el principio de conservación de
la energía mecánica.
Elija la posición inicial óptima para realizar este experimento.
a) (1.5 puntos) Para cada posición x buscamos medir la velocidad angular máxima
de la barra (al pasar por la vertical) al soltarla desde la posición inicial definida
por su grupo. Una vez grabado el video, la velocidad angular medida (ωmed ) la
obtendrán a partir de los cuadros más cercanos a la vertical utilizando el software
ImageJ.
Dado que la medida de ω no es demasiado precisa, buscamos mejorar los esti-
madores de la velocidad angular real, repitiendo el experimento 5 o más veces
y calculando el promedio. La desviación estándard de dichas medidas será un
estimador del error aleatóreo de sus mediciones. Se sugiere partir con el valor
máximo de x.
El resultado de este experimento debe ser cuatro valores (uno por cada configu-
ración x) de la velocidad angular ω al pasar por la vertical con su error aleatorio
(ver informe).
b) (1.5 puntos) Ahora buscamos calcular teóricamente la velocidad angular esperada
(ωcalc ). Se sugiere modelar la T como dos barras delgadas y despreciar la distancia
entre el rodamiento y el extremo de la barra larga. Está de más recordarles que
esta parte la pueden preparar antes de llegar a la sala Galileo.
Determine la posición del Centro de masa en función de x.
Determine el momento de inercia de la T en función de x.
Encuentre una ecuación para la velocidad angular esperada ωcalc .
Calcule en Matlab ωcalc (x). Se sugiere definir en Matlab un vector xcalc con
más puntos que las cuatro posiciones medidas, de modo de obtener una curva
más suave en el gráfico, ejemplo x-calc:0.05:0.025:0.55;
c) (1.0 puntos) Finalmente comparar el resultado experimental con el resultado teó-
rico.
Para ello haga un gráfico de ωmed vs x en Matlab en el que cada punto
medido aparece con su barra de error, utilize el comando errorbar. Adjuntar
al mismo gráfico la curva teórica ωcalc vs x con una línea contínua.
Ejemplo de gráfico
errorbar(x_exp,omega_exp,err_omega,’+’)
hold on
plot(x_calc,omega_calc)
Imprimir y Adjuntar al informe una copia impresa del gráfico bien rotulado
con unidades en ambos ejes, fecha y grupo.
3. Conclusiones [1 punto]
4C.1. Introducción
En la Unidad 4C se encontró la ecuación de movimiento para un cuerpo rígido que puede
girar en torno a un eje fijo en un movimiento plano. La ecuación de movimiento relaciona
la aceleración angular con el torque externo y es análoga a la ecuación de Newton, donde el
lugar de la masa es ocupado por el momento de inercia, la aceleración es reemplazada por
la aceleración angular y la fuerza neta por el torque externo.
Para esta experiencia cada grupo
contará con una polea masiva (por
111111111111
000000000000 111111111111
000000000000
lo tanto su momento de inercia no es ω ω
R1 M, I R1 M, I
nulo) que puede girar con poco ro-
ce en torno a un eje ubicado en su R2 R2
O O
centro. La polea tiene canales a dos
radios diferentes donde es posible en-
rollar cuerdas de las que cuelgan los
pesos, tal como se indica en la figu-
ra. Las cuerdas pueden estar enrolla- m1 m1
das en el mismo sentido (figura de la
izquierda) o en sentidos opuestos (fi- m2 m2
gura de la derecha).
El objetivo de esta experiencia es primero comprender las implicancias cualitativas de la
ecuación de movimiento de la polea para lo cual se buscará determinar bajo que condiciones
de distribuciones de masas, la polea gira en uno u otro sentido. Como segunda parte se
buscará que Uds. midan el momento de inercia de la polea usando la ecuación de torques; es
decir, midiendo simultáneamente los torques y la aceleración angular se calculará el momento
de inercia.
En este documento los guiaremos a través de las mediciones que esperamos que realicen y
su modo de presentacióon. El informe se entregará en un documento aparte más conciso.
A. Objetivos
Identificar los efectos del momento de inercia en la dinámica rotacional de los sólidos
rígidos
B. Materiales
Una polea masiva. Los radios externo e interno son R1 = 5 cm y R2 = 2,5 cm. La masa
no es conocida.
Una webcam
Matlab
ImageJ
El montaje experimental es muy sencillo, sólo deberán enrollar pesos de distinta masa en la
polea y medir el giro de la polea mediante una cámara web.
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario verifcar el funcionamiento correcto de los
elementos a usar. En esta experiencia usarán la cámara web y se deben seguir los
mismos procedimientos que la semana pasada para este efecto.
Idealmente la polea debe tener poco roce con el eje. Verifique que es así, para lo cual
hagala girar y deje que se frene libremente; repita la experiencia en ambos sentidos.
Normalmente debería tardar más de 30 segundos. Si no es así consulte con el profesor
o los auxiliares.
Para poder analizar correctamente el movimiento de la polea se debe escribir la ecua-
ción de movimiento de ésta. Para esto considere separadamente las dos geometrías
mostradas en la figura de más arriba y analícelas como se indica. Haga un DCL de
cada una de las masas que cuelga y de la polea. Aplique las ecuaciones cinemáticas
que permiten relacionar la aceleración lineal de cada masa con la aceleración angular
de la polea. A partir de estas ecuaciones se puede despejar la aceleración angular de
la polea en función del momento de inercia, los radios y las masas que cuelgan. ¿Cuál
es esta ecuación? Interprétela.
donde se han absorvido los diferentes términos que tienen t0 en las constantes A y B
y se ha definido C = α/2. En la experiencia ω0 = 0, pero φ0 y t0 serán típicamente
diferentes de cero, de manera que en ecuación (4C.3) los coeficientes A y B serán
diferentes de cero.
Para obtener el valor de α de los valores ti y φi medidos se usa Matlab:
En el informe indique:
Recomendaciones
No use una combinación de pesos que produzca un movimiento muy lento pues
en ese caso se torna dificil medir la aceleracion angular. ¿Por qué?
No use una combinación de pesos que produzca un movimiento muy rápido pues
en ese caso las fotos de la cámara salen borrosas.
Recuerde que los ángulos medidos con ImageJ tienen el cero en la horizontal
hacia la derecha (el eje x habitual), siendo positivos en el sentido antihorario y
negativos en sentido contrario. En el eje x negativo se produce la transición de
−180o a +180o . Luego, en algunos casos Uds deberán sumar o restar 360o para
obtener una curva continua.
No olvide convertir a radianes por segundo cuadrado su aceleración angular.
4. Conclusiones [1 punto]
B1. Materiales:
Tubos de PVC para formar tres anillos de igual diámetro y espesor de paredes.
Plano inclinado
Cuide que la distribución de masa garantice que se cumplan de mejor forma los supuestos
de la descripción teórica de una rueda rodando (material teórico). En particular, que el
centro de masas coincida con el eje del cilindro. Además, que la distribución de masa sea
simétrica con respecto al eje de simetría del cilindro. De no ser así, lo más probable es que el
cilindro tienda a desestabilizarce a medida que incrementa su velocidad. De los tres ejemplos
ilustrados en la figura de abajo, la configuración de la izquierda es aceptable. La del medio
es dudosa y la de la derecha es inaceptable.
C2. Registros
Posando uno de los anillos sobre el plano inclinado, en su extremos superior, comience a
grabar una âĂIJpelículaâĂİ y suelte el anillo para que ruede cuesta abajo.
Detenga la grabación una vez que el anillo llegue al extremo inferior del plano. Cuide que
el anillo baje âĂIJderechoâĂİ sobre el plano. Si se desvía, deseche esa película. Repita esta
grabación 3-4 veces.
Repita lo anterior para las otros dos anillos (produzca 3-4 películas cada vez).
C3. Análisis
Cargue ImageJ e importe cada una de sus películas.
ImageJ tiene una herramienta para medir distancias. La medición está en píxeles, así que
establezca la relación entre pixeles y metros midiendo con ImageJ la banda negra en el borde
del plano inclinado (0.1 m).
Luego mida y tabule en su cuaderno la posición del centro del anillo cada 3 cuadros con
respecto a un punto arbitrario pero fijo durante sus mediciones.
Los cuadros en ImageJ están cada 1/30 segundos. Entonces, para cada película Ud. tiene
una serie de pares (t, x) donde x es la distancia a lo largo del plano desde el extremo superior.
Ingrese esos datos a Matlab y haga un gráfico (t, x). Verifique que obtiene aproximadamente
una parábola y ajuste un polinomio de 2do orden (polyfit(x,t,2)), tal que x(t) = b0 +
b1 t + b2 t2 .
Relacione el coeficiente b2 con la aceleración del anillo (ax ) suponiendo un movimiento unifor-
memente acelerado. De esta forma, para cada anillo, Ud. obtendrá 3-4 valores experimentales
de ax .
A partir de los valores obtenidos experimentalmente para ax , infiera el valor experimental de
Io /M R2 en los tres casos de distribución de masa extudiados. En las situaciones idealizadas
ilustradas en la figura los cuocientes Io /M R2 son 1, 1/2 y 0, respectivamente.
11111
00000
00000
11111
00000
11111
00000
11111
00000
11111
Periferia Uniforme Centrada
C4. Verificación
Finalmente, realice una carrera entre sus tres anillos. Sueltelos desde el reposo en la parte
superior del plano inclinado y grabe su descenso. En este caso, ponga su cámara más arriba
del plano para tener perspectiva de este. Examinando los últimos cuadros de la película
determine en que orden llegan los anillos a la base del plano.
Realice la carrera y anote el orden de llegada unas 6 veces. Comente el resultado de las
carreras y explíquelos en términos físicos.
Identifique las fuentes de errores que afectan el resultado de la experiencia y comente sobre
como minimizarlos.
5A.1. Introducción
En la Unidad 5A estudiaremos las oscilaciones de un péndulo físico. El objetivo fundamental
es iniciar el estudio de movimientos periódicos, en particular el Movimiento Armónico Simple
(MAS). Buscamos que sean capaces de identificar la amplitud, fase, frecuencia angular y
constante de fase de un MAS, en particular para el caso de un péndulo físico.
En el caso particular de la Unidad 5A estudiaremos la posible variación del período de un
péndulo e físico con la amplitud y largo. Para ello utilizaremos en efecto un péndulo físico
compuesto por una barra homogénea de longitud 2L = 0,6m, y una barra del mismo material
pero largo L cruzada a la primera (esta última la denominaremos cruceta). Es decir, tenemos
una especie de T sujeta cerca de su extremo por un eje fijo en torno al cual es sistema podrá
girar libremente. El sistema es idéntico al usado en la unidad 4B.
Para llevar a cabo este experimento, utilizaremos de nuevo una cámara web para filmar
el movimiento del péndulo. Utilizando ImageJ buscamos que ustedes midan el período del
péndulo bajo una variedad de condiciones.
A estas alturas del semestre nos parece importante recalcar algunos conceptos básicos sobre
la medición de cantidades experimentales:
1. Cuando uno mide una cantidad física una gran cantidad de veces, utilizamos el valor
medio (promedio) de todas las medidas como un buen estimador de la cantidad física
que realmente queremos medir, de este modo evitamos que nuestras mediciones están
dominadas por errores experimentales. Por otro lado, una buena medida del error
aleatorio de nuestra medición está dado por la desviación estandard de las mediciones
(σ).
A. Objetivos
B. Materiales
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario verificar el funcionamiento correcto de los
elementos a usar. Hoy deberán medir repetidas veces el período de un péndulo. El
objetivo de esta parte es simplemente probar el equipo y estimar el error con el cual
pueden medir dicho período.
de haberlo olvidado, use el teorema de los ejes paralelos para simplificar su cálculo. Al
llegar a al laboratorio los alumnos ya debieran tener esta información preparada.
Adjuntar al informe una copia impresa del gráfico bien rotulado con unidades en
ambos ejes, fecha y grupo. ¿Es consistente el período medido con la aproximación
de pequeñas oscilaciones? ¿Hasta qué ángulo parece ser válida la aproximación de
pequeñas oscilaciones?
4. Conclusiones [1 punto]
5B.2. Materiales
Un globo e hilo (1 m).
Soporte universal.
Una vez obtenidas las constantes (A, τ , T , to ), verifique la consistencia entre ellas utilizando
la relación
1
Ω2 = ωo2 − 2 ,
4τ
discutida en el Material Teórico de esta unidad. En particular, se espera que con los valores
obtenidos de ω y τ , infiera la longitud del cordel que sostiene al globo. Comente críticamente
sus resultados.
p⇒q ⇐⇒ ∼ q ⇒∼ p .
2. Evite las corrientes de aire mientras hace las filmaciones. Para ello no transite rápida-
mente cerca del globo. Cierre puertas y ventanas si es necesario. Manténgase distante
del globo.
3. Enfoque la cámara hacia el soporte de donde cuelga el globo, cuidando que la extensión
filmada del cordel sea lo suficientemente larga como para hacer óptimas mediciones de
los ángulos. Si es necesario, use un fondo blanco para un buen contraste. El punto de
apoyo debe ser identificable en la imagen. Cuide de no saturar de luz la cámara.
1. El desplazamiento del globo en función del tiempo, según los valores medidos.
2. El desplazamiento del globo en función del tiempo, a si este cayese libremente, sin roce
(∼ gt2 /2).
3. El desplazamiento del globo en función del tiempo, a si este cayese con roce viscoso.
A fin de incluir en este análisis el caso de una fuerza de roce proporcional a v 2 , considere la
siguiente expresión para la velocidad vy (t) mostrada en el Material Teórico
1 − e−gt/Vt
vy (t) = Vt .
1 + e−gt/Vt
Aquí Vt corresponde a la velocidad terminal, el único valor ajustable en esta expresión. Con
esta función y un Vt estimativo (que Ud. debe ajustar para reproducir los tados), calcule
y(t) por recurrencia. Utilice la siguiente construcción en Matlab
%
% CAIDA CON ROCE -beta*v2 (INTEGRACION)
% t es un arreglo uniforme para el tiempo
tf=5 % Tiempo final
g=9.8
dt=0.08
ttt=0:dt:tf
vterm=2.4
v=(1-exp(-g*ttt/vterm))./(1+exp(-g*ttt/vterm));
yn(1)=0;
for i=2:length(ttt)
yn(i)=yn(i-1)+vterm*v(i)*dt
end
%
%
5B.7. El informe:
El informe reportado debe incluir
1. Dos gráficos, uno por cada práctica, en formato PDF. Cualquier otro formato será
rechazado (-2 puntos).
5B.7.1. PREGUNTA
¿Cuál sería el período de oscilación de un péndulo formado por el mismo globo de la primera
parte, pero inflado con un diámetro de unos 15 cm? Verifique experimentalmente.
5C.2. Materiales
Un riel, un carro, dos resortes, un motor y una barra en forma de "L", montado como
se ilustra en la figura.
Lentes de protección.
Matlab.
ω
m
k/2 k/2
M
0 x(t)
La barra en forma de "L" rotará en torno al eje del motor, a una frecuencia máxima de
rotación de aproximadamente 2 vueltas por segundo (2 Hz). Esta barra es la que actúa
como elemento de forzaje para el carro, de manera análoga al modelo simple expuesto en la
guía teórica. Esta barra se encuentra firmemente acoplada al eje de rotacón del motor. Se
exige el uso de lentes de protección durante el uso del equipo experimental.
A diferencia de otras sesiones prácticas, existe un sólo modelo de este experimento por mesa
de trabajo, por lo que un sistema de turnos deberá ser adoptado. La experiencia 1 se reali-
zará una sola vez por mesa, los tres grupos al mismo tiempo. Es una experiencia corta pero
necesaria para el buen desarrollo del resto de la sesión. Las experiencias 2 y 3 se pueden rea-
lizar en cualquier órden entre ellas, pero después de la experiencia 1. Es aquí donde deberán
implementar un sistema de turnos por grupo. La experiencia 4 debe realizarse al final, cada
grupo por su cuenta.
5C.4. Datos
Los datos conocidos del dispositivo experimental son los siguientes:
5C.5. Experiencias
Experiencia 1.- Medición de la frecuencia natural de vibración ωo .
Duración estimada = 15 min
Esta experiencia deberá ser realizada al comienzo de la sesión práctica y una sóla vez
por mesa de trabajo, es decir los tres grupos al mismo tiempo.
Para ello imponga una condición inicial al carro y mida la frecuencia de oscilación (en
unidades de rad/s). Se recomienda una velocidad inicial nula pero una posición inicial
entre 5 y 10 cm con respecto a su posición de equilibrio. En realidad, tratándose de
p
un oscilador amortiguado, lo que se mide es la frecuencia Ω = ωo2 − (1/2τ )2 , siendo
ωo la frecuencia natural de oscilación. Por ahora supondremos que la disipación es
pequeña por lo que Ω ≈ ωo . Esta suposición será validada después.
Para medir la frecuencia de oscilación utilice el cronómetro para medir el tiempo que
toma en realizar 10 oscilaciones. Cada grupo deberá realizar 4 medidas, lo que dará un
conjunto de 12 medidas por mesa, las cuales compartirán. Reporte el valor medio y la
desviación estandar de ωo . Además, determine el valor medio de la constante elástica
total k, suponiendo que los resortes son iguales de constante k/2 (para resortes en
paralelo se suman las constantes elásticas).
E4a ¿Cómo se comparan cualitativamente los valores de la amplitud del estado esta-
cionario obtenidos con el algoritmo de Verlet con lo esperado?
E4b ¿Cómo es la disipación en el experimento? Son realistas las suposiciónes Ω ≈ ωo
y ωo ≈ ωr realizadas en las experiencias 1 y 2 respectivamente?
B. Materiales
Matlab:
Iteraciones
Evaluación de funciones
Gráfico de funciones
Un sistema de varillas para producir ondas de torsión (demostrativo)
Uso de ImageJ para medir distancias
C. Experiencias
1. Preliminares [0 puntos]
Verifique que Matlab funciona y que sabe como correr M-Files.
2. Experiencia 1 [2 puntos]
Se vio en clases que
u(x, t) = f (x − ct) + g(x + ct)
es solución de la ecuación de ondas para cualquier función f y g.
En esta primera experiencia se pide graficar u(x, t) en función de x a medida que el
tiempo avanza para algunas funciones f y g dadas. A partir de los gráficos se pide que
deduzcan el significado de la solución de D’ Alembert.
Caso 1
1cm
f (x) =
(x/5cm)2 + 1
g(x) = 0
f (x) = 0
2cm
g(x) =
exp(x/10cm) + exp(−x/5cm)
c=1;
x=-100:1:100;
for t=0:2:40
y=x-c*t;
u=1./((y./5).^2+1);
plot(x,u)
pause
end
3. Experiencia 2 [1 puntos] Si
Vista en perspectiva
A partir de filmaciones de las oscilaciones del sistema completo, se busca que Uds
determinen el valor de I1 /I2 y lo comparen con la predicción teórica.
En el informe indique el procedimiento que usó para determinar I1 /I2 y el valor obte-
nido. Indique el efecto de tener un mayor o menor I en la propagación de las ondas
5. Conclusiones [1 punto]
1. Guı́a Práctica
A. Objetivos
Visualizar ondas viajeras y estacionarias en una cuerda tensa.
Usar las ecuaciones cinemáticas para ondas viajeras y estacionarias.
Determinar frecuencia y longitud de onda de modos normales de oscilación de una cuerda
finita.
B. Materiales
Video con onda viajera en un arreglo de varillas.
Generador de señales (frecuencia).
Cuerda fija en dos extremos.
Masas para tensar cuerda.
Balanza digital.
1. Medir la longitud de onda (en pixeles) en el sector de varillas largas tantas veces como
sea necesario de modo de obtener un error de medición relativo menor a un 5 %. Para
ello puede medir cuadros distintos del video o bien medir un sectores distintos en cada
cuadro, siempre y cuando sean todas las varillas del mismo largo. Reporte el número de
mediciones Nl , la longitud de onda λl [pixeles] y su desviación estándard σl [pixeles].
1
Semestre 2008-1 Sistemas Newtonianos: Unidad 6B 2
1. Variar lentamente la frecuencia del generador de pulsos, fijar una amplitud pequeña
pero lo suficientemente grande para visualizar la amplitud de las ondas en la cuerda
(viajeras o estacionarias). Determinar la frecuencia de oscilación para los tres primeros
modos normales de la cuerda montada estime su error para anotar en el informe. Mida
la longitud de onda en cada modo normal.
Conclusiones [1 punto]
Se reconozca el valor de la fuerza media (Fm ) como cantidad física para caracterizar
estos sistemas.
7A.2. Materiales
• Placa plana rectangular;
• Soporte universal;
• Hilo (1 m);
• Caja plástica de contención;
• Sensor de fuerzas;
• Transportador;
• Huincha de medir;
• Plomada;
8
Entendemos por pulso al chorro de extensión longitudinal acotada. La duración de ellos es corta en
relación a las escalas de tiempos características de las cantidades físicas a analizar.
(T − To ) = Bh ,
donde
8 ρ A⊥ g
B= .
3
De acuerdo a esto se espera que el incremento T − To crezca linealmente con la altura h
y sea independiente del ángulo φ. Se examina en esta sesión el grado de validez de este
comportamiento y además se busca una estimación de la constante B.
A fin de reducir la carga horaria de esta sesión, quienes estén en los mesones A, B, C, D y J
se encargarán de la Experiencia 1, mientras que quenes se ubiquen en los mesones E, F, G,
H e I, se encargarán de la Experiencia 2.
SENSOR
TUBO
GARBANZOS
HILO h
SOPORTE PLACA INCLINADA
CAJA CONTENCION
En la Experiencia 1 se medirán las fuerzas medias para caídas desde {20 cm, 30 cm, 40 cm,
50 cm}; todas ellas a un ángulo fijo entre 35o y 50o de inclinación. Cada grupo escogerá (e
informará) el valor que le parezca. Para cada una de estas alturas se harán 2 mediciones y
se calculará el promedio de ellas.
En la Experiencia 2 se medirán las fuerzas medias para inclinaciones de {35o , 45o , 55o };
todas ellas a una altura cercana a 40 cm (que debe ser reportada). Para cada uno de estos
ángulos se harán 2 mediciones y se calculará el promedio de ellas.
7A.6. Procedimiento
1. Ajuste y mida la inclinación de la placa;
5. Active el sensor de fuerzas y haga caer los garbanzos. Cerciórese que han impactado la
placa en el lugar correcto, ubicado en la tercera parte (L/3) desde su extremo superior;
7A.8. El informe
Para que el informe se considere completo este debe incluir
4. Indicar al menos dos (2) fuentes de errores en este experimento, sugiriendo como
corregirlas.
5. Proponer al menos una (1) vía para mejorar el modelo colisional presentado.
7A.8.1. No olvide:
F(t)
11111111111
00000000000
00000000000
11111111111 Fm
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
ta t b = t a+ ∆t ta tb
Para simplificar ideas supongamos que en este intervalo se identifican N valores de la fuerza
en instantes uniformemente espaciados (δti = δt), entonces definimos la fuerza media como:
N
X N
X N
1 X
Fm δti = Fi δti ⇒ Fm = Fi .
N
i=1 i=1 i=1
El significado geométrico detrás de esta definición es que el área bajo la curva F (t) vs t,
entre ta y tb , es igual a la de Fm en el mismo intervalo.
Recordando que la suma de variaciones consecutivas es igual a la variación total:
X
δpi = ∆p = pb − pa ,
i
obtenemos
∆p = Fm ∆t .
P
Aquí hemos identificado ∆t = i δti = (t1 − t0 ) + (t2 − t1 ) + . . . + (tN − tN −1 ) = tb − ta .
P
A la cantidad J ≡ i Fi δti se le denomina impulso, que en el contexto de la Segunda Ley
de Newton resulta igual al cambio de momentum. Entonces,
J = Fm ∆t = ∆p ,
>> Media(’colisionH45b.txt’,250,105)
En ella se ve claramente que hay que aumentar el segundo argumento, por ejemplo 250 →
280. Además, hay que achicar el tercer argumento 105 → 85. Lo anterior a fin de tomar el
promedio en la región en que la señal es prominente, que debe coincidir con la curva roja.
Note que en esta figura, al lado derecho se observa una oscilación amortiguada. Lo más
probable es que se trate de oscilaciones de la placa luego de la caída de los grabanzos. Esa
parte no cuenta para el promedio.
>> Media(’colisionH45b.txt’,250,105)
Durac: 0.210 s
Fondo: 2.421 N
Media: 2.355 N
FzaMedia: 0.066 N
10
La palabra output aparece en el Diccionario de la Real Academia de la Lengua Española:
http://www.rae.es/
3. Cuando tengan la configuración experimental lista, realice la medición (corra Run Once).
4. Para exportar los datos, en la ventana Data View copie los datos al Clipboard: con el boton
derecho en la ventana negra seleccione
Export To->Clipboard (Text)
5. Abra el Block de Notas y pegue los datos del Clipboard. Elimine las lı́neas de encabe-
zado y grabe el archivo.
Reconozca que los cuerpos sumergidos experimentan una fuerza de empuje vertical,
en sentido opuesto al peso del mismo.
7B.2. Materiales
Un soporte universal
Matlab y SignalExpress
7B.3. Datos
Aceleración de gravedad en Santiago: g = 9,796 m/s2 .
7B.4. Procedimientos
En esta práctica se realizarán medidas de fuerza (peso), longitud y masa. Las medidas de
fuerza se realizarán usando el sensor, la tarjeta NI-USB-DAQ y el programa SignalExpress.
En U-Cursos encontrará un programa llamado Empuje-SignalExpress.seproj que le permitirá
tomar 1 s de medidas a una velocidad de 200 medidas por segundo. Para cada medida de
fuerza deberá reportar el valor promedio y el error obtenido con esta serie de medidas.
7B.5. Experiencias
Experiencia 1.- Calibración del sensor de fuerza en el rango ±10 N.
Duración estimada = 30 min
El buen desarrollo de esta sesión requiere de una calibración de cada sensor de fuerza.
Esto significa que cada grupo tendrá que obtener las constantes de conversión A y
B de la relación lineal F = A · U + B de su sensor de fuerza. Es importante notar
que la masa del gancho del sensor, que llamaremos mG , no es conocida. En principio
se debe tratar de determinar esta masa. Para ello siga el procedimiento indicado a
continuación.
1. Mida el voltaje correspondiente al peso del gancho del sensor, primero con el sen-
sor apuntándo hacia "arriba", y después, hacia "abajo". Denotaremos al conjunto
voltaje-fuerza (U− , F− ) y (U+ , F+ ) respectivamente.
2. Para al menos 4 masas diferentes, obtenga una medida del valor medio y error
absoluto de cada voltaje asociado. Reporte los datos en una tabla. Agregue a esta
tabla los valores obtenidos en la parte anterior.
3. A través de una regresión lineal obtenga los valores de las constantes A y B. Para
ello siga estos pasos: (1) Dado los valores promedios de cada voltaje asociado a
cada masa se obtiene el siguiente conjunto de ecuaciones
F− = −mG g = A · U− + B,
F+ = +mG g = A · U+ + B,
F1 = (m1 + mG )g = A · U1 + B,
F2 = (m2 + mG )g = A · U2 + B, etc...
donde m1 , m2 , ... son las masas utilizadas. (2) A pesar de que no conocemos
mG se puede realizar una regresión lineal con el conjunto de datos (U1 , m1 g),
(U2 , m2 g), etc.. pero con las constantes A y B 0 ≡ B − mG g. (3) Determinados A
y B 0 , obtenga B y mG con las dos ecuaciones para F− y F+ .
Nota: La regresión lineal se debe hacer con la función polyfit de Matlab, expli-
cada en la guía práctica de la unidad 4C. Si tiene dudas de su utilización escriba
help polyfit en la línea de comando de Matlab.
Figura 30: Esquema del montaje experimental. Para la experiencia 1, donde se calibrará el
sensor, las masas no deben ser sumergidas en la cubeta de agua.
1. Usando la balanza digital, mida por separado la masa de la esfera y del cilindro.
Realice 5 medidas y, reporte un valor promedio y su error absoluto. Si la desviación
estándar es nula, considere la precisión del instrumento.
2. Usando el pie de metro, mida por separado el volumen de la esfera y del cilindro.
Realice 5 medidas y, reporte un valor promedio y su error absoluto. Si la desviación
estándar es nula, considere la precisión del instrumento.
3. Usando los datos obtenidos, determine la densidad de cada uno de los objetos
(ρo ).
4. Con el sensor de fuerza mida el peso de cada objeto, que denotaremos P , sin que
estos estén sumergidos bajo el agua. Compruebe que P = Mo g, donde Mo es la
masa del objeto.
5. Para cada objeto sumergido, mida el empuje E que experimenta. Compare el
comportamiento de sus mediciones E versus V con el Principio de Arquímedes,
es decir E = ρf V g.