Apuntes
Apuntes
Apuntes
Departamento de Física
Facultad de Ciencias Físicas y Matemáticas
Universidad de Chile
2 Agosto 2010
Indice
1
Información General Sistemas Newtonianos 2
4A.4.Ejemplo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4A.5.Preguntas Conceptuales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
4A.2.Guía Práctica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
4C.5.Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
4C.6.Ejemplos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
5B.8.Ejemplo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180
5B.10.
Ejercicios Semestres Pasados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
5C.10.
Preguntas Conceptuales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 202
5C.11.
Ejercicios Semestres Pasados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
5C.2.Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
5C.4.Datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 205
5C.5.Experiencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 206
7A.2.Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 253
7A.6.Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 255
7B.2.Presión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 262
7B.2.Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 268
7B.3.Datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 268
7B.4.Procedimientos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 269
7B.5.Experiencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 269
5C.2.Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 308
5C.4.Datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 309
5C.5.Experiencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 309
7A.2.Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 319
7A.6.Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 321
7B.2.Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 328
7B.3.Datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 328
7B.4.Procedimientos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 329
7B.5.Experiencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 329
10
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. En algunos casos en que esto sí 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
gracar diversos modelos físicos.
En esta unidad se ilustrará cómo se utiliza el computador, mediante el programa Matlab con
estos nes.
Ejemplo:
Consideremos el ejemplo de una persona parada sobre una
tarima. La persona está sujeta a dos resortes de largo nat-
ural L, cada uno de los cuales está amarrado a la misma
altura de la persona en direcciones opuestas y que inicial-
h
mente no están elongados y en posición 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.
11
Métodos Numéricos Sistemas Newtonianos 12
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
p
−mgh + k( h2 + L2 − L)2 = 0
√
Si esta expresión se expandiera daría lugar a una ecuación de tercer grado en la variable h 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 = 80 kg, L = 10 m , k = 10N /m y g = 9,8 2
m/s , 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
jos. Para eso, lo mejor es hacer un gráco 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 ≈ to ⇒ x(t) ≈ x(to )
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 el 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 ser lmado con una cámara proporciona un muestreo cada
cierta fracción de segundo solamente.
x(t + h) − x(t)
ẋ(t) = lı́m ,
h→0 h
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á jo, 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 − 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 − 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
k
xi+1 = 2xi − xi−1 − xi ∆t2 (1.10)
m
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 cuáles 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 xo y vo . Recordando que la velocidad es la primera derivada de la posición y usando la
expresión (1.6) se tiene
vo = ẋ(0)
x1 − xo
=
∆t
de la que se despeja
x1 = xo + vo ∆t (1.11)
Se tienen así todos los pasos del llamado algoritmo de Verlet (L. Verlet, Physical Review
1967)
• Dados: xo , vo
• Se calcula: x1 = xo + vo ∆t
• Se itera desde i = 1 hasta el tiempo nal:
k
xi+1 = 2xi − xi−1 − xi ∆t2
m
En Matlab, considerando que m=1 kg, k = 0,5N/m, x0 = 5,0 m y v0 = 2,0 m /s, eso se
hace de la siguiente forma:
>> 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) el 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
xi − xi−1
xi+1 = 2xi − xi−1 + F xi , ∆t2 /m (1.13)
∆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.
Pregunta 3: Si x(t) = sen(t) exp(−0,1t), escriba las líneas de Matlab que permitirían gracar
esa función en el intervalo t ∈ [0 : 10] usando un espaciado ∆t = 0,1.
Escriba la ecuación de energía que permite determinar el punto donde se detiene la partícu-
la.
Pregunta 5: Si x(t) = cos(t) exp(−2t), escriba las líneas de Matlab que permitirían gracar 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.
Para eso, resuelva numéricamente la ecuación de Newton que resulta con los parámetros
m = 1 kg y H = 10m con γ = 0; 0,1 kg/s; 0,2 kg/s; . . .; 0,5 kg/s.
Resuelva para m =0,1 kg, Vo =1 m/s y γ= 0,1 kg/s. Compare con la predicción sin roce
H= V02 /2g .
Ejercicio 3: Se desea resolver el movimiento de la Tierra en torno al Sol. Se sabe que en ese
caso la fuerza es la de gravitación universal:
GM m
F~ = − r̂
r2
Con el n de poder tratarla numéricamente, la fuerza se reescribe de la siguiente manera
(considerando el movimiento en el plano x − y)
GM m
F~ = − r̂
r2
GM m
= − ~r
r3
GM m
= − (xî + y ĵ)
r3
GM m
= − 2 (xî + y ĵ)
(x + y 2 )3/2
GM
ẍ = − x (1.18)
(x2 + y 2 )3/2
GM
ÿ = − 2 y (1.19)
(x + y 2 )3/2
Use el método de Verlet visto en clases para resolver estas ecuaciones acopladas. Considere
los siguientes valores de las constantes: G=1 y M = 1. Además considere como condición
inicial para la posición x0 = 1 e y0 = 0 y para la velocidad vx0 = 0, vy0 =0.5, 1.0 y 2.0.
Graque la trayectoria que resulta (plot(x, y)). Compare con los cálculos analíticos que
predicen, para
p los datos del problema, que la velocidad para una órbita circular es Vcirc =
GM/R = 1.
mv̇ = −γv
vi+1 − vi−1 γ
= − vi
2∆t m
Muestre que si resuelve para m = 1, γ = 0,1, ∆t = 0,1, Tfinal = 100 y V (0) = 1 el resultado
no tiene sentido.
Sin embargo, se puede escribir otra discretización centrada, en que el lado derecho se
promedia en dos instantes
vi+1 − vi γ vi+1 + vi
=−
∆t m 2
de la cual se puede despejar vi+1 como
1 − γ∆t/2m
vi+1 = vi
1 + γ∆t/2m
muestre que esta discretización es estable y entrega resultados sensatos.
mẍ = −kx
Resuelva la ecuación de movimiento, es decir calcule x(t), usando el método de Verlet para
el siguiente conjunto de valores: m = 1kg , k = 1N/m, x0 = 3m y v0 = 0.
Una vez que tenga la solución de la ecuación, estudie si numéricamente la energía mecánica
se conserva
1 1
E = mv 2 + kx2
2 2
Para eso escriba de alguna manera discreta la velocidad y evalue la energía en función del
tiempo. ¾Es constante?
Ejercicio 6: Se desea saber cómo disminuye en el tiempo la energía mecánica de un cuerpo que
cae en el aire en presencia de roce turbulento Froce = −γ|v|v , donde |x| es el valor absoluto
de x.
Para eso, considere que se suelta un cuerpo de masa m desde una altura H del piso y se
sigue su evolución hasta que golpea al suelo.
Ejercicio 7: Se desea modelar un billar donde las bolas se mueven en un plano con roce viscoso
sobre la supercie 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 simplicar 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 entrega
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, Santi-
ago es plano).
2. Complete el siguiente programa Matlab que hace el control llamando al método Alarma
que hace sonar la alarma:
for i=1:fin
if(XXXXXXXX)
Alarma;
end
end
g
φ̈ = − sin φ − γ φ̇
L
donde L es el largo del péndulo y γ el coeciente de roce.
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 gura) serán cada vez
menores.
Se busca resolver numéricamente la dinámica del sistema para obtener cómo van dismin-
uyendo 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.
A. Objetivos
1. Conocer las capacidades de los métodos numéricos en la solución y análisis de los sistemas
newtonianos.
B. Materiales: Matlab
C. Experiencias
L
Experiencia 1.- Una partícula de masa m puede deslizar sin roce
por un riel horizontal. La partícula está unida a dos resortes V
idénticos de constante elástica k y largo natural L. Los ex-
tremos de los resortes están jos a una distancia L y 2L del
riel, tal como se indica en la gura. Inicialmente se le da a la
partícula una velocidad V hacia la derecha. Se pide determinar
el alcance de la partícula si L = 0, 1m, m = 0, 1kg, k = 2N/m, 2L
para los siguientes valores de V : 0,3 m/s, 0,5 m/s, 1 m/s.
Velocidad V Alcance
26
Métodos Numéricos Sistemas Newtonianos 27
ϕ̈ = −(g/L) sen(ϕ)
π
Considere que el péndulo se suelta desde el reposo en un ángulo inicial ϕo = 4
Use g = 9, 8m/s2 , L = 0, 5m, ∆t = 0, 01s, 0, 05s, 0, 5s Integre la ecuación hasta T = 10s.
Se pide gracar la solución para los diferentes valores de ∆t y estudiar la conabilidad de
la solución.
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 ϕo = 10 , 4 , 2
Complete la siguiente tabla:
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 sigue el período al aumentar
el ángulo inicial.
¾Cómo gracar 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 .
¾Cómo hacer un gráco de datos con barras de error?
Utilice errorbar(x,y,err,'*'), con err el arreglo de errores asociado a y, en las coorde-
nadas x.
¾Cómo incluir leyenda de los símbolos usados en un gráco?
Utilice legend('leyenda 1, leyenda 2, ... ') después de la última instrucción plot.
El orden de las leyendas sigue el orden de las curvas desplegadas.
F. Lecturas recomendadas
Manual de Matlab
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íco, 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 vericar las predicciones. En el caso de que
éstas no se veriquen, los nuevos resultados experimentales sirven de observación (se vuelve al
paso 1) para poder plantear nuevas hipótesis.
Mediremos fuerzas estáticas y dinámicas, es decir fuerzas constantes o que varíen en el tiempo
respectivamente. Una manera simple de medir fuerzas es usando masas conocidas. Por ejemplo,
30
Métodos Experimentales Sistemas Newtonianos 31
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.
si se quiere aplicar tensión a una cuerda podemos colocar una cuerpo de 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 , supuesto que se desprecie la masa de la cuerda. 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 gauge en inglés, consta de una peque na
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, nalmente, el voltaje medido es
proporcional a la fuerza aplicada.
F = A · U + B,
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 conver-
sió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 sucientemente 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 denir una velocidad. En otras aplicaciones se medirá el ángulo recorrido por un
objeto entre dos cuadros consecutivos, obteniéndose una medida de velocidad angular.
Supongamos que realizamos un conjunto de medidas en forma controlada, bajo las mismas
condiciones experimentales. Estas condiciones se reeren al tipo de instrumento, el tipo de escala
usado con ese instrumento, que determina en general su precisión, como también a las condiciones
de ambiente o del sistema como temperatura, humedad, luminosidad, por nombrar algunas. En
general las condiciones del tipo de instrumento y su uso sí son determinantes en el grado de
repetibilidad de las medidas, pero algunas de las condiciones de ambiente mencionadas (u otras)
puede que no tengan un efecto signicativo.
Tenemos entonces N medidas de una cantidad física C , obtenidas bajo las mismas condiciones
experimentales. Supondremos que estas medidas, que no son iguales, se distribuyen en torno a un
valor. Los valores medidos siguen una distribución de probabilidad, o de ocurrencia. Usualmente,
cuando N es un número grande, digamos sobre 30, los valores siguen en buena aproximación
una distribución particular llamada gaussiana, conocida también como distribución normal. Es-
ta distribución tiene una forma analítica precisa, que no detallaremos por ahora, pero que es
caracterizada por dos parámetros: el valor promedio y el ancho de la distribución. Las N medi-
das experimentales serán entonces representadas por su valor medio hCi y su desviación estándar
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
σ, que es una medida del ancho de la distribución. Es importante notar que el promedio hCi
es una estimación estadística del valor real de la cantidad C. Mientras más medidas tomemos,
mientras más datos tengamos, mejor será esta estimación.
Dos ejemplos de estas distribuciones de ocurrencia son presentadas en la gura 33, con N=
100 y N = 10000 medidas de una cantidad C. El eje horizontal C muestra los valores obtenidos
y el eje vertical muestra la ocurrencia de estos valores, es decir cuántas veces, dentro del número
total de medidas N, se encuentran valores en un rango particular de C.
Si llamamos Ci los valores medidos de C, con i = 1...N , el valor promedio y la desviación
estándar que caracterizan la distribución de estos valores tienen las siguientes deniciones:
N
1 X
hCi = Ci , (2.1)
N
i=1
v
u
u1 X N
σ = t (Ci − hCi)2 (2.2)
N
i=1
De los datos de la gura 33 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 signicativas, 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
Gaussiana de medidas, aproximadamente 68 % de los datos se encuentran en el rango [hCi −
σ, hCi + σ], es decir dentro del ancho dado por 2σ . También se tiene que aproximadamente 95 %
de los datos se encuentran en el intervalo [hCi − 2σ, hCi + 2σ].
Pregunta 2: Usando los datos presentados en la gura 33 con N = 100 ¾Cuál es el porcentaje
de valores que se encuentran fuera del rango [hCi − 3σ, hCi + 3σ]?
Errores aleatorios: como dice su nombre estos errores se producen fortuitamente. Para
minimizar estos errores conviene realizar un número considerable de mediciones y prome-
diar los resultados. Como los errores aleatorios van en ambas direcciones con la misma
probabilidad, al promediar los datos, éstas contribuciones aleatorias se promedian a cero.
La palabra precisión se usa para referirse a una medición con bajo error aleatorio.
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
denirá la dirección vertical de la gravedad.
Los errores aleatorios se pueden presentar por un problema de medición, como por la precisió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 hacemos 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 estatura 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 0, 12 m.
Pero si tomamos un niño para medir su estatura, este valor no será una buena estimación.
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 .
Error relativo: corresponde al cociente entre el error absoluto y el valor medio, de la
forma
∆C
c = .
hCi
Por ejemplo, un error relativo de 0,05 signica que hCi, que es una representación estadística
de la cantidad C, es una buena estimación dentro del 5 % de error.
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
Pregunta 3: Para las cantidades U = 0,55 ± 0,04 V, A = −4,9 ± 0,1 N/V y B = 12,25 ± 0,05
N, obtenga una fuerza con su error asociado. ¾Cuáles son las fórmulas de propagación de errores
que debe utilizar? A y B corresponden a los factores de calibración del sensor pero con errores
asociados.
Deniremos las cifras signicativas 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 signicativas. Si por algún motivo, la medida hubiese dado 14,4±1 cm, entonces lo correcto
es decir que tiene dos cifras siginicativas 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 nal 15 ± 1
cm.
Un error clásico, que lamentablemente se acarrea hasta muy adelante en las carreras cientí-
cas 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 estudiante suele entregar C = 12,23453 ± 0,01236. ¾Se tiene la
precisión para entregar todas esas cifras? ¾Cuáles son las cifras signicativas 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 la forma correcta de expresar el
resultado C = 12,24 ± 0,01, y el resultado tiene cuatro cifras signicativas.
Cuando se realizan operaciones artiméticas con varias cantidades físicas, entonces las cifras
signicativas estarán dadas por el error nal, 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.
Nota: el concepto de cifra signicativa no está realmente bien denido, puesto que este número
depende de la elección de la base del sistema de numeración, que es arbitrario. No obstante es
de suma utilidad prà ctica, ya que es fácil determinar, como en el ejemplo anterior, si una cierta
cantidad de decimales es excesiva en relación a la precisión de un experimento.
>> d = load('datos.txt');
Esta instrucción crea una variable llamada d, que tiene una estructura de vector de dimensiones
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á congurado
Esto se corrige
para grabar los números con comas, lo cual puede provocar errores en Matlab.
cambiando la opción en Conguración regional del sistema operativo.
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 gura 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. Normalmente una
velocidad de 100 a 1000 datos por segundo debería ser suciente 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).
Pregunta 6: Suponga que las constantes de conversión del sensor de fuerza son A = −4,9
N/V y B = 12,25 N. Si se usa la tarjeta en el rango ±5 V, ¾cuál es la resolución nal en unidad
de fuerza? ¾Es mejor o peor que la sensibilidad del sensor de fuerza?
Pregunta 2: Describa las diferencias entre los errores de medición sistemáticos y los aleatorios.
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 conocidos
se puede medir una fuerza.
Pregunta 7: Si se cambia la tarjeta de adquisición de datos por una más barata de 8 bits, ¾Es
posible aun tener una sensibilidad de 0,01 N? La relación fuerza voltaje es F = −4,9U +
12,25 (para rango ± 10), donde U es el voltaje medido y F es la fuerza en N.
2. Dena en lenguaje Matlab las constantes del problema y arreglos para el tiempo,
posición, velocidad, y energía.
P1. [1 pt] Un péndulo simple se usa para medir la aceleración de gravedad, usando
p T =
2π L/g . El período medido es T = 1,51 ± 0,03 s y la longitud, L = 56,7 ± 0,2 cm.
¾Cuál es el valor resultante de g, su error absoluto y relativo?
P2. [1 pt] Se da un valor como 35,562 ± 0,2. Reescríbalo con el número adecuado de cifras
signicativas. 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 un termopar 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 el
termopar está dada por una ecuación a determinar.
El cuadro 1 tabula estas medidas. A partir de estos datos graque la curva temperatura
versus voltaje en el rectángulo cuadriculado que se encuentra en el reverso de esta hoja.
Suguiera y justique 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).
Medida N
◦ T (N) ∆T (N) Observaciones
1 11.7 0.3
2 16.9 0.2
3 12.0 0.5
4 3.5 1.5 Medida mal realizada: hilo se salió de la polea
5 18.2 0.5
6 14.0 0.2
7 15.1 0.1
8 19.4 0.3
9 12.7 0.3
10 10.4 0.4
11 11.2 0.5
12 12.7 0.1
13 14.9 0.2
Ejercicio 4: Un sensor de fuerza tiene una relaación de calibración dada por F = A·U2 +B ·U,
donde U es el voltaje medido y F la fuerza correspondiente. El sensor puede ser utilizado
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 de
las constantes A y B en cada rango de medida. Se sabe que estas constantes son ambas
positivas. ¾En qué se diferencia este sensor del modelo que se usa en este curso?
Nota: Por ahora preocúpese solamente de los valores promedio" 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 la 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
Con el resultado anterior, determine la velocidad media con su error asociado. Se sabe
que la distancia entre estos dos puntos es de L = 800 ± 5 m.
' ?Qué errores aleatorios o sistemáticos pueden ser importantes en este proceso de
medición?
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 signicativas de su resultado
nal.
a/b − c
9,81 · exp(−b · c/10)
a + 3,2 · b − c
Ejercicio 7: Suponga que se realizó el proceso de calibración de un sensor de fuerza (P3), con
lo cual se determinaron las siguientes constantes de calibración en el rango de ±10 N:
A = 4,9 ± 0,1N/V y B = 12,25 ± 0,05 N. Llene la tabla que se muestra a continuación y
determine la fuerza promedio y su error absoluto.
Medida N
o U (V) δU (V) F (N) ∆F (N)
1 3,8 0,1
2 4,5 0,2
3 3,3 0,3
4 3,1 0,1
5 4,1 0,2
a/b + c
5,5a exp(−bc)
a ln(2b) + c
a + 7b − c
abc
a − 3 sen(bc)
Ejercicio 9: Considere los siguientes valores de velocidad entregados por una persona que no
sabe de cifras signicativas ni errores. Corríjalos para que sean consistentes
A. Objetivos
Aprender a usar el programa Measurement & Automation para vercar que el sensor de
fuerza está funcionando correctamente
B. Materiales
C. Experiencias
46
Métodos Experimentaless Sistemas Newtonianos 47
Ahora verique que el sensor de fuerza que está conectado con el canal de entrada análogo
número0 (AI0) está trabajando adecuadamente. Para ello lance la aplicación Measurement
& Automation (programa de medida y automatización), cuyo ícono se encuentra en el
escritorio de su PC.
Al abrirse la Barra del Menú Principal, seleccione Conguration → Devices and Interfaces
→ 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.
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 Conguration: 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 adquisició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
Analog Input y Save to ASCII. Estas son tareas ya preasignadas. Usted puede cambiar los
parámetros pero antes de hacerlo anote los valores preasignados pues en principio están
denidos para un buen uso para la experiencia 3.
Con Analog Output y Analog Input se conguran la salida AO0 y entrada AI0 de una manera
muy similar a lo que se hizo en la experiencia 2.
Con la botonera Save to ASCII puede elegir un nombre de archivo como también un
directorio donde guardarlo. Tenga cuidado de mantener el formato de archivo con ASCII
(formato texto), pues así podrá leer facilmemte los datos con Matlab. Verique que su
conguración de Windows graba los números con punto decimal y no coma.
Ahora pruebe la medida de tensión de corte de un hilo, lanzando la medición con el botón
Run Once. De hecho, esto lanza las tres tareas mencionadas en forma consecutiva.
El hilo debe atarse un extremo al sensor de fuerza y el otro extremo pasarlo por la polea y
anudarlo al dedo de la persona que tirará rme, pero suavemente incrementando la mag-
nitud de la tensión aplicada. Se puede visualizar la adquisición seleccionando la lengueta
Data View y simplemente arrastrando el cursor desde Analog Input hacia esta ventana.
Puede también leer estos datos desde Matlab como se indica en el texto de material teórico
sobre Métodos Experimentales.
Con esta prueba determine el rango en el cual usará el sensor, ±10 N ó ±50 N. Explique
la elección del rango:
Con los datos obtenidos por su grupo realice un histograma que muestre la distribución de
ocurrencia de la tensión de corte del hilo.
Con los datos obtenidos por su mesa realice un histograma que muestre la distribución de
ocurrencia de la tensión de corte del hilo.
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.
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 automóvil,
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 aproximaciones,
y por último, plantear representaciones simplicadas del sistema. Todo esto es en principio de
gran dicultad, 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 está limitada a ciertos
aspectos y ciertos rangos de valores.
En un intento de clasicació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 distancia de
separación puede variar moderadamente ante deformaciones.
iv.- sólidos indeformables, donde las distancias entre moléculas vecinas son invariable en el
49
Sistemas Extendidos Sistemas Newtonianos 50
tiempo.
Figura 4: Choque de dos galaxias interpretado a partir de datos obtenidos por el telecopio espacial
Hubble.
Para estudiar la evolución de un sistema formado por muchos componentes nos valemos de
las leyes de Newton, las que se adaptan de forma muy sencilla a sistemas complejos. Debemos
considerar, provisoriammente, al sistema complejo como un conjunto de sistemas más pequeños,
para los cuales las leyes de movimiento se puede aplicar en forma sencilla. Para esto es necesario
denir el sistema a estudiar, vale decir, el conjunto de partículas distinguibles y enumerables
que uno escoge (cticiamente) para el estudio de un fenómeno dado. El criterio que comunmente
nos guía para denir los componentes del sistema es la simplicidad que nos brinde para respon-
dernos preguntas especícas. Suponemos que las componentes del sistema pueden interactuar
entre sí, y a su vez con el exterior. Es decir, las fuerzas actuando sobre cada constituyente corre-
sponden a la suma de las fuerzas que se ejercen directamente más las fuerzas que las otras partes
ejerzan sobre él. La mecánica newtoniana de sistemas compuestos consiste en una aplicación
de las tres Leyes de Newton, además del Principio de Superposición, postulados para cuerpos
puntuales. Estas leyes se resumen como sigue:
I.- Si un cuerpo en reposo no interactúa con el entorno, entonces este permanecerá en reposo
indenidamente.1
II.- El cambio de momentum δ~p de una partícula es proporcional a la fuerza aplicada y a la
duración δt de su aplicación.2
1
A la inversa, si un cuerpo cambia de su estado de reposo, entonces es porque interactuó con el entorno.
2
Esta relación se expresa δ~p = F~ δt → F~ = d~p/dt = m~a.
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 reere a que las fuerzas son aditivas en el sentido vectorial.
(a)
(b)
+ + + + + + + + + + (c)
(d)
N
~ =
X mi
R ~ri .
M
i=1
Tan solo al mirar la denición podemos interpretar geométricamente el punto ~ como el promedio
R
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
~ = (Rx , Ry , Rz ) ,
R (3.1)
entonces,
N
1 X
Rx = mi xi , (3.3)
M
i=1
N
1 X
Ry = mi yi , (3.4)
M
i=1
N
1 X
Rz = mi zi . (3.5)
M
i=1
n n
1 X 1 X M iL L
Rx = mi xi = × =
M M n+1 n 2
i=0 i=0
111111111111111111
000000000000000000
M, L
0 xi= i d
N PN
i=1 mi yi
X
Ug = mi gyi = M g = M gYCM ,
M
i=1
donde claramente YCM representa la coordenada Y del centro de masas del sistema.
Si R
~ A localiza el centro de masas de un sistema A de masa MA
yR~ B el de un sistema B de masa MB ,
entonces el centro de masas del conjunto está dado por
~ ~
~ = MA RA + MB RB
R
MA + M B
X X X
~ =
MR mi~ri = mi~ri + mi~ri .
todos i∈A i∈B
~ A = (P
MA R ~ B.
Pero, i∈A mi~
ri ), análogamente para R Con ello,
~ = 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 gura 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
3
Homogéneo signica que su composición es idéntica en todas las partes del mismo. En un sentido estricto
no existen materiales homogéneos, son una idealización. En un sentido práctico podemos decir que un material
es homogéneo cuando las diferencias de composición entre sus partes son tan pequeñas que son incapaces de ser
detectadas en un experimento especíco. Es decir: ½Un mismo objeto puede ser homogéneo o no dependiendo del
experimento!
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
N N N
!
~
~ = dR = d
V
1 X
mi~ri =
1 X d~ri
mi =
1 X
mi~vi =
1 ~
P .
dt dt M M dt M M
i=1 i=1 i=1
| {z }
~
P
Entonces
P~ = M V
~ (3.11)
~
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 éste pueda experimentar.
j
i
fj/i
tario, denotemos f~j/i la fuerza que ejerce la componente j sobre la i, naturalmente f~i/i = 0,
i-ésima es N ~
P
entonces la fuerza neta de todas las componentes sobre la j=1 fj/i . La ecuación de
movimiento de Newton para la partícula i es:
N
d~
pi X
= F~i + f~j/i . (3.12)
dt
j=1
N N N X
N
X d~
pi X X
= F~i + f~j/i . (3.13)
dt
|i=1{z } |i=1
{z } i=1 j=1
| {z }
~
dP ~ ~0
F
dt
PN d~
pi ~
dP
la suma de derivadas es igual a la derivada de la suma, por lo cual i=1 dt = dt ;
la suma de todos los pares de fuerzas internas es nulo, en virtud al principio de acción y
~i/j
reacción (f = −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, indepen-
diente 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.
La energía cinética es una cantidad aditiva, vale decir, si representamos el sistema como un
conjunto de N celdas, cada una de masa mi , con i : 1 → N, entonces la energía cinética total
Considerando que la celda i-ésima experimenta describe un trayecto de radio ρi con velocidad
angular ω, entonces su rapidez es ωρi . Sustituyendo en la expresión para la energía cinética
obtenemos "N #
1 X
K= mi ρi ω 2
2
2
| i=1 {z }
I
Aquí hemos denido el momento de inercia del sólido en torno al eje de rotación, que denotamos
por I:
N
X
I= mi ρ2i ,
i=1
con ρi la distancia al eje de la celda i-ésima, de masa mi . Esta denición se extiende a cualquier
sólido distribuido volumétricamente.
Observaciones:
2. No hay restricción a la orientación ni dirección de los ejes con respecto al cual se evalue el
momento de inercia.
ω
ρi ωρi
1
0
0
1
Figura 8: Un cuerpo plano rotando en torno a un eje jo; a la derecha se ilustra una vista desde arriba.
Apéndices
3.5. ¾Altera un espectador pasivo la dinámica de un sistema?
El centro de masas de un sistema, sea este cohesionado, disgregado, amorfo o deformable, es
un punto. Lo interesante es que cual sea la naturaleza del sistema y de sus fuerzas internas, el
movimiento (aceleración, velocidad y trayectoria) de su centro de masas está determinado sólo
por la acción de las fuerzas externas.
Un ejemplo curioso que permite ilustrar este punto lo encontramos en un sistema formado
por dos masas idénticas, A y B, sin que interactúen entre sí y cada una de masa m. Ambas posan
sobre una supercie horizontal lisa y el centro de masas se ubica en el punto medio entre ambas.
Para jar ideas, alineamos ambas partículas en un eje horizontal x, con xA y xB las coordenadas
de A y B, respectivamente.
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
Para calcular el área del círculo de radio R, lo subdividimos en N anillos circulares, cada uno
de ancho d = R/N . Si N es lo sucientemente grande, el área de un anillo de semiradio r es,
en buena aproximación, igual a 2πrd. Esto se ilustra en la gura de más abajo. La construcción
r=d/2:d:1 dene un arreglo de semiradios d/2; 3d/2; ....
El uso de esta misma idea en el cálculo del volumen de una esfera se traduce en subdividirla
en cascarones de igual espesor. El volumen de una cascara de 'semiradio'r y grosor d será 4πr2 d.
El volumen total será la suma de estas contribuciones.
2 πrd
1111111111111
0000000000000
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
r 0000000000000
1111111111111
r
0000000000000
1111111111111
0000000000000
1111111111111
r
0000000000000
1111111111111
0000000000000
1111111111111 2π
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
r=d/2:d:1
d franja extendida
Pregunta 2: Dena el momento de inercia de un sólido e indique al menos una de sus propiedades.
Pregunta 3: En no más de 5 líneas, establezca dos diferencias físicas (½no algebraicas!) entre el
momento de inercia de un sólido y su centro de masas.
Pregunta 4: Excluyendo el caso de un objeto puntual, ilustre con un ejemplo un caso de mo-
mento de inercia nulo.
Pregunta 5: En no más de 5 líneas, establezca dos diferencias físicas (½no algebraicas!) entre el
momento de inercia de un sólido y su momentum.
Pregunta 6: Considere un lápiz de grato 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.
Identique las fuerzas actuando internamente en el sistema y las 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 9: En no mas de 5 líneas, describa dos diferencias físicas entre la masa de un sistema
y su centro de masas.
Pregunta 11: Un canasto con manzanas reposa sobre una mesa. En el contexto del sistema
formado por el canasto y las manzanas, indique dos fuerzas externas y dos fuerzas internas.
Pregunta 12: Describa dos ejemplos de cuerpos cuyos centros de masa se ubican en el exterior
de ellos. Justique su respuesta.
Ejercicio 2: En la gura se ilustra un sólido conformado por dos barras perpendiculares. Una de
ellas es de longitud L y masaM , y la transversal es de longitud L/2 y masa M/2. Ambas
barras se unen a una distancia x del extremo O de la barra mayor. Calcule y graque la
posición del centro de masas del sistema como función de x. Rotule adecuadamente su
gráco, indicando los puntos mas característicos a resaltar.
Ejercicio 3: En la gura se muestra un círculo uniforme de radio R con una perforación cuadra-
da como se indica. La longitud de cada lado del cuadrado es b y su centro dista R/2 del
centro del círculo. Determine la ubicación del centro de masas del círculo perforado.
Ejercicio 4: Se tiene una barra de longitud L homogénea y doblada en V como indica la
gura. El ángulo de doblado de la barra es α.
1. (2pt.) Determine a que distancia del vértice (donde se dobla la barra) se encuentra el
centro de masa del objeto.
A. Objetivos
B. Materiales
C. Experiencias
A Calcular numéricamente las coordenadas (Xcm , Ycm ) del centro de masa de un arco
cerrado (α = 2π ) usando M =1.AA kg y R=1.BB m, con AA y BB los dos últimos
dígitos del RUT de integrantes distintos del grupo. Estudie el caso N = 100.
B Ahora se desea calcular el centro de masa de medio aro, es decir cuando α = π . Calcule
y tabule las coordenads del centro de masas para los casos
p θi = 0, π/4, π/2, 3π/4, π .
Reporte 2 + Y 2 /R
Xcm e interprete. En todos estos casos use N = 100.
cm
65
Sistemas Extendidos Sistemas Newtonianos 66
Experiencia 2: Considere una barra uniforme de masa M y longitud L que gira con velocidad
angular ω en torno a un eje perpendicular. El eje se ubica en la barra, a una distancia λL
(λ ∈ [0 − 1]) de uno de sus extremos. Por lo visto en los apuntes y por análisis 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 coeciente adimensional que depende de los valores de λ, vale decir, γ = γ(λ).
Utilizando la misma estrategia de discretización de la experiencia 1, calcule y graque γ
para los valores λ=0, 1/4, 1/2, 3/4 y 1. Use N = 100. Incluya en su gráco los resultados
para λ = 0, ZZ , con ZZ los dos últimos dígitos de todos los integrantes del grupo.
D. Análisis e informe
E. Lecturas recomendadas
4A.1. Introducción
En este capítulo continuaremos el desarrollo 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, éste puede rotar sin modicar la posición de su CM. Estudiaremos la leyes de la estática,
que nos permitirán evaluar cuándo un cuerpo permanecerá en reposo. Las leyes de la estática
son uno de los conocimientos empíricos más antiguos de nuestra civilización. Este ha permi-
tido la construcción de monumentos formidables que aún, con todos los recursos tecnológicos
disponibles, nos sorprenden. El descubrimiento básico fue que era posible amplicar, 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 .
La validez de estas leyes empíricas encuentra su fundamentación en las leyes de Newton, que
como hemos visto se aplican a objetos puntuales. Veremos que las leyes de la estática aplicadas
a un sólido se resumen en dos restricciones:
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.
67
Estática Sistemas Newtonianos 68
La fuerza necesaria para cerrar una puerta es menor si se ejerce lejos de las bisagras. Por
el contrario mientras la fuerza se ejerza más cerca del borde jo de la puerta, más costará
cerrarla.
Estas observaciones constituyen la motivación fundamental para el concepto que brinda el título
a esta sección: el torque asociado a una fuerza. Si bien aún no hemos denido torque, algo que
corregiremos enseguida, podemos describirlo como fuerza de palanca asociada a una fuerza. La
idea es simple: consideremos un balancín equilibrado como el de la gura. Para simplicar la
discusión supongamos que la barra horizontal no tiene masa (masa despreciable). Los bloques
en los extremos forman parte de lo que denominaremos `sistema', de modo que el contacto de
ellos con la tabla constituye una fuerza interna. Las fuerzas externas al sistema son: 1) el peso
del bloque izquierdo; 2) el peso del bloque derecho y; 3) la fuerza (normal) donde se apoya la
tabla. La primera exigencia es que la suma vectorial de todas ellas sea nula.
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
~ =A
C ~×B
~,
es un vector cuya
B
A
1. ~×A
A ~ = 0, para todo ~;
A
2. ~×B
A ~ = −B
~ ×A
~;
3. ~ × (λB)
A ~ = λ(A
~ × B)
~ ; con λ un escalar;
4. ı̂ × ̂ = k̂ ; ̂ × k̂ = ı̂; y k̂ × ı̂ = ̂;
5. ~×B
A ~ = AB sin θ = A⊥ B = AB⊥ , donde A⊥ es la magnitud de la componente de ~,
A
~.
perpendicular a B
6. ~ × (B
A ~ + C)
~ =A
~×B
~ +A
~×C
~;
4
En ella, ambos vectores se disponen con sus colas coincidentes. Luego, el dedo anular de la mano derecha
toma la dirección y sentido de A~ , con la palma orientada hacia B
~ . El sentido de A ~ coincide con el del dedo
~×B
pulgar.
~ = (Ay Bz − Az By , Az Bx − Ax Bz , Ax By − Ay Bx )
C (4A.1)
~τ = ~r × F~ ,
donde ~r es el vector que une el origen 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
N N
!
X X
~τ = ~ri × (mi~g ) = mi~ri ~ CM × (M~g ) ,
×~g = R
i=1 i=1
| {z }
~ CM
MR
Consideremos un sólido sobre el cual actúan N fuerzas externas: F~1 , F~2 , · · · F~N .
Para que el sistema esté en equilibrio estático es necesario que se cumplan las siguientes condi-
ciones:
III.- ~v = ~0;
IV.- ~ P = ~0.
ω
Donde hemos indicado explícitamente el punto P con respecto al cual se determina los torques
y la velocidad angular.
OBSERVACIONES:
1. La suma de fuerzas es vectorial, de modo que muchas veces resulta conveniente expresar-
la 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 simplique 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 especicar si el torque
entra al plano.
Uno puede convenir el signo (+) del torque para cualquiera de estos casos, quedando los
torques en sentido opuesto con signo ().
4A.4. Ejemplo.
Una carga de masa m reposa sobre un tablón de longitud L y masa M distribuida uniforme-
mente. El tablón está sostenido en sus extremos por dos cuerdas ideales las cuales forman ángulos
θA y θB con respecto la vertical. La carga se ubica en una posición tal que permite que el tablón
se mantenga en forma horizontal. a) Determine la ubicación de la carga con respecto al extremo
A del tablón. b) Verique e interprete concisamente su respuesta para el caso θA = θB .
θ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:
L (m + M ) sin θA cos θB M
x= −
m sin(θA + θB ) 2
L (m + M ) sin θA cos θA M L
x→ − ⇒x→
m 2 sin θA cos θB 2 2
vale decir, si los ∠'s son iguales, la carga debe ubicarse simétricamente en el tablón.
Pregunta 2: Para el diagrama de la gura, identique las fuerzas actuando sobre la barra y
P~
plantee (sin despejar) la ecuación Fi = 0.
Pregunta 3: Una viga de largo L y masa M está rmemente empotrada a una pared vertical
sobre la supercie de la Tierra. Dibujar la viga, un sistema de referencia, y listar los
componentes de la reacción de la pared sobre la viga en dicho sistema.
Pregunta 4: Explique si es posible que una escalera de largo L y masa M descanse en equilibrio
estático sobre una supercie horizontal sin roce y una supercie vertical con roce.
Pregunta 5: Una viga de densidad uniforme tiene masa M y largo L. Si la viga se apoya en un
punto jo al piso a una distancia L/4 de un extremo, ¾Qué fuerza es necesario aplicar y
dónde para que la viga esté 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 gura:
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 gura. 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
conguraciones desde la más eciente a la menos eciente 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 π/4 rad a un
metro de uno de sus extremos y se cuelga de una cuerda como muestra la gura. 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 gura. ¾En cuál caso el centro de masa del sistema adquiere una mayor
rapidez?
en caso (a);
en caso (b);
es igual;
P2. [5 puntos] Un vaso cilíndrico homogéneo (abierto por arriba) de radio basal a y altura
b está posado sobre un plano inclinado en un ángulo β c/r a la horizontal y no resbala
debido a un pequeño tope (de altura despreciable) jo sobre el plano. Suponga que el vaso
se puede modelar como compuesto de una lámina sin espesor y de densidad supercial
(masa por unidad de supercie) 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, 25 m como indica la gura?
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 supercie de corte vertical. A n de
evitar que las semi-esferas se separen, se dispone una cuerda ideal con dos cuerpos dénticos
de masa m cada uno como indica la gura. 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.
4A.1. Introducción
En la Unidad 4A buscamos estudiar las condiciones de equilibrio estático de una barra hor-
izontal sujeta por una cuerda tensa. Las condiciones que se deben cumplir para mantener el
equilibrio se pueden resumir imponiendo que la fuerza neta y el torque neto sobre la barra
sean nulos. En un problema en dos dimensiones, estas condiciones en general se convierten en 3
ecuaciones que se deben cumplir entre las fuerzas y sus puntos y direcciones de aplicación.
Para esta experiencia cada grupo contará con una barra sujeta en un extremo por un tornillo
jo a un soporte vertical. La barra se sostendrá en posición horizontal por un hilo conectado al
sensor de fuerzas que utilizaron hace dos semanas. Por último, se colgarán masas bajo distintas
posiciones de la barra. Cada grupo medirá entonces la tensión de la cuerda bajo una variedad de
conguraciones.
En general, mediciones como las que harán esta semana se utilizan con una variedad de
objetivos en física o ingeniería. A veces se busca encontrar algun tipo de ecuación o ley de la
naturaleza que reproduzca las mediciones realizadas. Por ejemplo podemos buscar si la fuerza
ejercida por un medio elástico es proporcional al estiramiento desde su posición de equilibrio, ¾o
quizás al cuadrado del mismo? También nos puede interesar medir en qué rango de estiramientos
es válida la aproximación matemática encontrada, quizás debido a deformaciones permanentes o
incluso rompimiento del medio elástico en estudio. Otras veces se busca medir alguna propiedad
de los cuerpos tal como podría ser la aceleración gravitacional en el laboratorio, o la tensión de
corte de un hilo como en la Unidad 2 de este curso. En estos últimos casos se busca encontrar
el mejor parámetro que reproduce las observaciones, esto se denomina ajustar una función a los
datos y lo llevarán a cabo más adelante.
El objetivo de esta semana es más sencillo, y simplemente les pedimos comparar las medidas
realizadas por ustedes con los valores esperados de acuerdo a las ecuaciones de la estática. Es
decir, ustedes deberán comparar la medida de la tensión del hilo realizada con el sensor de
fuerzas con la tensión calculada a partir de las leyes de la física (las leyes de estática en este
caso). La conclusión más fundamental a la que esperamos ustedes lleguen es si ¾las leyes de la
estática reproducen bien la tensión medida bajo una variedad de condiciones o no? Sin embargo,
no esperamos simplemente una respuesta positiva o negativa sino una comparación cuantitativa
(diferencia absoluta o relativa) entre la tensión medida y la tensión calculada. En caso de que
79
Estática Sistemas Newtonianos 80
la diferencia encontrada sea superior a lo que ustedes sospechan debiera darles, busquen alguna
explicación plausible para tal efecto.
Es sumamente importante que comprendan que en caso de haber diferencias entre los valores
medidos y los valores calculados, éstas se pueden deber a una gran cantidad de efectos sistemáticos
o aleatorios que pueden ser distintos para cada grupo. Esto hace que las conclusiones a las que
llegue cada grupo no deben ser necesariamente las mismas.
B. Materiales
El montaje experimental es similar al ilustrado en la gura 37, 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 conguración de equilibrio.
C. Experiencias
1. Preliminares
Al inicio de cada experimento es necesario vericar el funcionamiento correcto de los el-
ementos 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é in-
stalado correctamente. Recuerden que los cables del sensor deben estar conectados a la
tarjeta de medición en las siguientes puertas:
Luego Abrir Test Panel e ingrese 5V bajo Analog Output y presione Update. Final-
mente seleccionar Analog Input y vericar que el sensor da una medida cercana a 2,5V
sin fuerza aplicada.
b) Cerrar la aplicación Measurement and Automation y abrir Signal Express para medir
señales de fuerza. Modicaremos el programa Tensión de Corte que se utilizó hace
dos semanas. Hoy mediremos tensión bajo condiciones estáticas, luego no buscaremos
encontrar una tensión mínima o máxima en un intervalo de tiempo, pero sí obtener
la tensión con la mayor precisión posible. Para ello realizaremos una gran cantidad de
mediciones con Signal Express.
Para simplicar el procedimiento no grabaremos los datos obtenidos con Signal Ex-
press sino simplemente obtendremos el valor medio y la desviación estándar de una
gran cantidad de medidas de la tensión.
En esta etapa cada grupo deberá decidir si utiliza el sensor de fuerzas en el rango de
validez de 10 N o de 50 N.
La tensión del hilo y su error se obtienen de la estadística que entrega el programa Signal
Express. Recuerden que para encontrar la tensión T a partir del voltaje promedio medido
V se utilizan las constantes A y B de cada sensor de fuerzas:
T = A ∗ V + B[N ]
Para encontrar la masa de la barra deberá primero hacer un diagrama de cuerpo libre y
resolver las ecuaciones de la estática. ¾Con cuántas cifras signicativas conocen la masa de
la barra? ¾Con respecto a qué punto conviene calcular el torque y porque?
¾Tiene sentido la masa encontrada, es decir, está en un rango más o menos esperado? Si
obtienen valores claramente erróneos tales como 1 mg o 103 kg es casi seguro que tienen
un problema de unidades en sus ecuaciones.
Tomar una masa y medir la tensión de la cuerda al variar la distancia de dicha masa con
respecto al eje del soporte. Medir la tensión en cinco o más posiciones distintas de la carga.
Gracar la tensión en función de la posición junto con su barra de error.
Al hacer un gráco lo más común es unir los puntos con una línea continua, dando la
apariencia de haber tomado un gran número de mediciones. En esta experiencia buscaremos
gracar exclusivamente la tensión medida con su error. En un macro de Matlab (archivo
.m) dena arreglos con las cantidades a gracar, por ejemplo:
errorbar(pos,tension_exp,error_tension,'+').
Noten que hemos utilizado un simbolo + para indicarle a matlab que ponga dicho símbolo
en cada punto gracado y no los una con una línea. Podríamos haber utilizado algún otro
signo. El comando errorbar graca 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 ), denan 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 esperada del hilo 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 programas.
Gracar la tensión calculada como una línea continua junto al gráco anterior.
Ejemplo de gráco
errorbar(pos,tension_exp,error_tension,'+')
hold on
plot(pos_calc,tension_calc)
Adjuntar al informe una copia impresa del gráco bien rotulado con unidades en ambos
ejes, fecha y grupo.
de enganche del hilo sobre la barra. Con una masa ja a una distancia también ja del eje,
deben medir la tensión para al menos tres posiciones distintas del anclaje.
Para cada posición, midan también el ángulo entre el hilo y la barra, ¾con que precisión
creen que pueden medir dicho ángulo? Gracar la tensión en función del ángulo, de nuevo
se les pide gracar exclusivamente los puntos medidos con el comando errorbar de Matlab.
En el mismo gráco dibujar el resultado esperado teóricamente suponiendo que conocen
todos los parámetros experimentales de manera exacta.
5. Conclusiones [1 punto]
Wneto = Kf − Ki (4B.1)
La evaluación del trabajo neto requiere integrar (es decir, sumar con mucho detalle) los dW
entre una posición inicial y otra nal a lo largo de la trayectoria de la partícula. Sin embargo,
existen fuerzas en que el trabajo neto NO depende del camino recorrido sino que exclusivamente
de la posición inicial y nal de la partícula. Estas fuerzas se denominan conservativas e incluyen
al peso y la fuerza elástica. En contraste, la fuerza de roce entre una partícula y una supercie
(usualmente escrita como −µc N ) no es conservativa.
Como para una fuerza conservativa, F~c , el trabajo efectuado es solo función de las coorde-
nadas, se puede denir una función de energía potencial, U, tal que el trabajo realizado es igual
a la disminución en la energía potencial:
Z
Wc = ~ = −∆U = Ui − Uf
F~ dx (4B.2)
Consideremos que sobre una partícula solo actúa una fuerza conservativa F~c . Entonces, com-
binando (4B.1) y (4B.2) obtenemos que ∆K = Wc = −∆U , de donde:
∆(K + U ) = 0 (4B.3)
85
Energía de Rotación Sistemas Newtonianos 86
que es una forma de escribir la ley de la conservación de la energía mecánica. Podemos además
denir la energía mecánica total del sistema como E ≡ K + U.
Como E permanece constante a medida que la partícula se desplaza, cualquier cambio de la
energía cinética ocurre a expensas de igual cambio (pero de signo opuesto) en la energía potencial.
~nc ,
Si sobre el sistema actúa más de una fuerza conservativa y otras fuerzas no-conservativas (F
como el roce), obtenemos una ecuación más general:
X X
W (F~nc ) = ∆(K + U) (4B.4)
Como mencionábamos anteriormente, tanto el peso como la fuerza elástica son fuerzas con-
1 2
servativas, y sus energías potenciales son Ug = mgy y Ue = 2 kδ , respectivamente. Aquí y
corresponde a la altura sobre un nivel de referencia (arbitrario pero jo). Si la partícula está
sobre (bajo) ese nivel y ≥ 0 (y ≤ 0). Además, δ representa la deformación (compresión o defor-
mación) de un resorte adherido a la partícula, g es la magnitud de la aceleración de gravedad y
k la constante elástica del resorte.
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 nal (objeto impacta
el suelo):
1 1
mv 2 + mgH = mvf2 + mg0
2 0 2
de donde obtenemos directamente
1
vf = (2gH) 2
Cada partícula esta sometida a fuerzas externas (por ejemplo, el peso) y fuerzas internas (por
ejemplo, las fuerzas que ejercen las partículas vecinas). Se puede demostrar (se verá con detalle
en el curso Mecánica) que las fuerzas internas en un sólido rígido no realizan trabajo neto.
Es decir, las fuerzas internas pueden transferir energía de una partícula a otra, pero lo que una
gana lo pierde exactamente la otra.
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
que corresponde a la conservaciíon de energía de un sólido. En las siguientes secciones se verá
que esta expresión se simplica aún más.
∆(K + Ug ) = 0 (4B.9)
con K= Iω 2 /2 y Ug = M YCM g .
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.
Figura 12: Momento de Inercia de una lámina, una esfera llena y un cascarón esférico
~ = 1
X
R mi r~i (4B.10)
M
También es conveniente recordar que al considerar un sistema compuesto por varios cuerpos,
podemos calcular primero la posición del CM de cada cuerpo, y luego sumarlos para encontrar
el CM del sistema. En otras palabras, cada cuerpo es tratado como una partícula con la masa
concentrada en su propio CM y luego empleamos la ecuación (4B.10).
Ejemplo: un sistema formado por una barra de largo L y masa M distribuída en forma uni-
forme, 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íca del sistema en este instante. Más importante, sabemos que el CM se ubica a
5
6 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æ̂
Por denición IO =
P
mi x2i , y ICM = Σmi ri2 . Expresamos el vector ~ + r~i ,
x~i = R luego
~ + r~i )2 = R
x~i 2 = (R ~ 2 + 2R
~ · r~i + r~i 2
X
IO = ~ 2 + 2R
mi [R ~ · r~i + r~i 2 ]
X
mi r~i = 0
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
1 2
momento de inercia de una barra de largo L y masa M con respecto a su CM es
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
I B = 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 = 13 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 nal 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 nal y la velocidad nal del centro de masa
ωf = (3g/L)1/2
1 1
vCM f = Lωf = (3gL)1/2
2 2
¾Qué pasa si la barra rota en torno al punto marcado por una x?
Consideremos un cuerpo rígido sobre el cual solo actúa el peso y que puede rotar en torno a un
eje. Si el cuerpo se suelta desde el reposo, podemos evaluar la velocidad angular en cualquier otro
instante (por ejemplo, cuando su CM pasa por su punto más bajo) empleando la conservación
de la energía mecánica
1/2
2M g∆y
ωf =
I
donde ∆y es el cambio en la posición del centro de masa. Como I es proporcional a M, el
resultado anterior es independiente de la masa del cuerpo. De esta forma, la velocidad angular
nal crece con el desplazamiento del centro de masa (se esta transformando energía potencial
gravitacional en energía cinética) pero disminuye con I. Si realizamos el mismo experimento pero
Figura 15: Dos cuerpos de igual masa pero distinta forma, ¾cuál cae más rápido?
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
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
Otro ejemplo interesante es el momento de una placa uniforme de masa M , ancho a y longitud
b. Sabemos que el momento de inercia de una barra de longitud L con respecto a un eje que
pasa por uno de sus extremos, perpendicular a la barra, vale M L2 /3 (revisar Unidad 3, Sistemas
Extendidos). No es difícil convencerse de que este momento de inercia es el mismo al de una
puerta de longitud L, con respecto a un eje que pasa por un eje en su base (ver gura de más
abajo).
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 gura 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’
111111111111
000000000000
z
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’
0000000000000
1111111111111 x’
y
x (xi,y i) 0000000000000
1111111111111
0000000000000
1111111111111
i
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
xi
0000000000000
1111111111111
0000000000000
1111111111111
x y
X X
I⊥ = mi x2i + mi yi2 = Iyy0 + Ixx0 ,
i i
2
P
demostrándose la propiedad. Nótese que hemos denotado Ixx0 = i mi yi . Es un buen ejercicio
convencerse de esta relación.
Abordemos un ejemplo simple. Una placa cuadrada de lados de longitud a está pivoteada en
una de sus esquinas mediante un eje perpendicular al su plano. Inicialmente la placa se suelta
del reposo como se indica en la gura de más abajo. Queremos saber su velocidad angular una
vez que ha girado en 90 grados.
En el estado A, K = 0 y Ug = M ga/2.
ω 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 momento
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.
Determine la velocidad angular del disco cuando el resorte alcanza su largo natural.
4B.1. Introducción
En la Unidad 4B comenzamos a estudiar la dinámica de sólidos rígidos. Esta semana estudi-
aremos 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 jo en torno al cual es sistema podrá girar libremente. Ojo que
el mismo sistema lo volverán a utilizar en tres semanas más.
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.
100
Energía de Rotación Sistemas Newtonianos 101
B. Materiales
El montaje experimental es muy sencillo, sólo deberán cambiar la posición de la cruceta para
realizar las distintas mediciones. Noten que al variar la posición de la cruceta, varía tanto la
posición del centro de masa de la T, como su momento de inercia (en este caso en torno al eje
de rotación del extremo superior).
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario vericar el funcionamiento correcto de los el-
ementos 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?
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 denida 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 estimadores
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 conguració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.
Para ello haga un gráco 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áco la curva teórica ωcalc vs x con una línea contínua.
Ejemplo de gráco:
errorbar(x_exp,omega_exp,err_omega,'+')
hold on
plot(x_calc,omega_calc)
Imprimir y Adjuntar al informe una copia impresa del gráco bien rotulado con
unidades en ambos ejes, fecha y grupo.
3. Conclusiones [1 punto]
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
~˙ = ~r × F~
L (4C.4)
~˙ = ~τ
L (4C.5)
105
Torque y Momento Angular Sistemas Newtonianos 106
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
~˙ = mR2 θ̈k̂
L (4C.7)
2
= mR αk̂ (4C.8)
Como T~ es paralelo al vector posición su torque se anula y sólo queda el peso. Se debe notar que
ésta es una de las comodidades del método de torque, pues se pueden eliminar algunas fuerzas
del análisis (la tensión en este caso). El torque del peso se obtiene de la denición: la magnitud
de ~r es R y el ángulo que forman el peso con el vector posición es θ. Al aplicar la regla de la
mano derecha se obtiene que el torque apunta perpendicularmente al plano, pero hacia adentro.
Es decir
~τ = Rmg sen θ(−k̂) (4C.10)
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.
X
~ =
L ~i
L (4C.13)
i
X
= mi~ri × ~vi (4C.14)
i
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 consideraremos sólo el
movimiento de rotación en torno a un punto jo. En la próxima unidad se verá el caso general.
En el caso de tener un sólido rígido que gira en torno a un punto jo 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 jo es constante ρi , describiendo un movimiento circular de ese radio en torno al punto
jo.
V
ω
111111111111
000000000000
O
ρi vi
ω
De esta forma, el momento angular total del sólido que rota en torno a un punto jo 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 jo. Al pasar de la tercera a la cuarta línea se identicó el momento de inercia
respecto al punto jo O. Además, se denió 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.
~˙ =
X˙
L ~i
L (4C.20)
i
X
= ~τi (4C.21)
i
= ~τtotal (4C.22)
donde ~τi es el torque sobre cada partícula y ~τtotal es la suma de los torques.
~˙ = 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.
~˙ = ~τtotal
IO ω (4C.26)
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 jo, un cuerpo de mayor
momento de inercia tendrá una menor aceleración angular. Es decir, el momento de inercia indica
la dicultad 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 acelerará 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?
Sea un sólido que está descrito como sistema de N partículas de masas mi y posiciones
~ri . En general sobre cada una de las partículas se ejercerán fuerzas provenientes de las otras
partículas del sólido (por ejemplo, las fuerzas moleculares que lo mantienen rígido) y fuerzas que
son ejercidas por otros cuerpos. A las primeras se les llamará fuerzas internas y se denotará por
f~ik la fuerza que la partícula k le ejerce a la partícula i. Al segundo tipo de fuerzas se les llama
fuerzas externas y se les denotará por F~iext . Así la fuerza total sobre la partícula i es
X
F~i = f~ik + F~iext (4C.28)
k
El torque total tiene una componente interna y otra externa. Calculemos primero la compo-
nente interna: XX
int
~τtotal = ~ri × f~ik (4C.31)
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
2 ).
Fik = −Gm1 m2 r̂ik /rik Debido a esa propiedad, los productos cruz en (4C.35) son
todos nulos.
En consecuencia, el torque total de las fuerzas internas es nulo. Esta propiedad es muy
importante porque implica (de acuerdo a la ecuación (4C.26)) que las fuerzas internas no provocan
aceleración angular. Dicho de otra forma, un cuerpo no se pone a girar de manera espontánea.
donde se deja explícita la indicación que el brazo de los torques se mide respecto al punto jo O.
4C.5. Resumen
En resumen, la ecuación de torque para un sólido rígido respecto a un punto jo O se escribe
de las siguientes maneras
~˙ = ~τOext
IO ω (4C.38)
4C.6. Ejemplos
4C.6.1. Movimiento del péndulo físico
Se llama péndulo físico al caso de un sólido rígido de forma arbitraria que puede girar libre-
mente respecto a un punto jo bajo la acción de la gravedad.
θ
G
mg
En la gura se representa el péndulo físico que está sujeto del punto jo 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 jo ~.
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 jo y el signo se obtuvo de la regla de la
mano derecha.
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é?
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 jo 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?
T = mg − mRα (4C.49)
La aceleración de la masa m es
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:
m1
. m2
Tipler y Mosca, Física para la Ciencia y la Ingeniería, Secciones 9.3, 9.4, 9.5, 10.1 (op-
cionales 10.2 y 10.3)
Cualquier libro de Física en los capitulos de Rotación o Dinámica Plana de Sólidos Rígidos.
Pregunta 3: Para el movimiento de un sólido rígido con punto jo, 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.
45 45
Ejercicio 3: Una rueda de masa M, radio R y momento de inercia I está en contacto con una
supercie horizontal rugosa, caracterizada por un coeenete 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 ja a la misma altura del centro de la rueda.
Ejercicio 4: Considere el sistema de la gura, 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
ella en O. En el extremo P de la barra se encuentra un cuerpo de masa m y en el extremo
Q otro de masa 2m. Todo el sistema se mueve bajo la acción de la gravedad.
Ejercicio 6: Considere la T como la de la gura, que puede rotar en torno a un eje que pasa
a distancia d del vértice.
2. (3pt.) Suponga ahora que el eje tiene un roce viscoso, de manera que al girar la barra
ejerce sobre ésta un torque τroce = −γω . Escriba la ecuación de movimiento de la
barra, es decir una ecuación del tipo
θ̈ = f (θ, θ̇)§
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 jo 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.
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.
120
Torque y Momento Angular Sistemas Newtonianos 121
Identicar 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.
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 vericar el funcionamiento correcto de los el-
ementos 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. Verique que es así, para lo cual hágala
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.
En el informe indique las combinaciones probadas, el sentido de giro que se observó (vector
hacia afuera o adentro o ausencia de movimiento) y el signo predicho por la ecuación de
movimiento de la polea. Compare.
Nota: Para que las cuerdas no deslicen por la polea (y así garantizar que se cumplen las
relaciones cinemáticas usadas en el análisis del problema), éstas se deben enrollar con varias
vueltas alrededor de la polea. No haga un lazo pues normalmente quedan ojos y la cuerda
desliza.
Se abre ImageJ y en el menu Import se lee la película como Avi. Hecho esto, es
posible analizar la película cuadro a cuadro.
sale el número de cuadro, que se convierte a tiempo recordando que cada cuadro se
toma cada 1/30 s.
Con esto se obtiene una serie de ángulos para diferentes tiempos {ti , φi }N
i=1 . Para que
la medición sea representativa se deben obtener unos diez valores.
Se sabe que el movimiento de la polea en estas condiciones (ver Material Teórico) ocurre
con aceleración angular constante α. Luego, el ángulo evoluciona en el tiempo como un
movimiento uniformemente acelerado
1
φ = φo + ωo (t − to ) + α(t − to )2 (4C.1)
2
1 2
= A + Bt + αt (4C.2)
2
= A + Bt + Ct2 (4C.3)
donde se han absorbido los diferentes términos que tienen to en las constantes A y B y se
ha denido C = α/2. En la experiencia ωo = 0, pero φo y to serán típicamente diferentes
de cero, de manera que en ecuación (4C.3) los coecientes A y B serán diferentes de cero.
Se pide a Matlab que encuentre la mejor parábola (el mejor polinomio cuadrático)
que pase lo más cercano por los puntos. Este procedimiento se llama ajuste o t y se
hace en Matlab con la función polyfit que entrega los coecientes del polinomio
coef=polyfit(t,phi,2)
coef =
XX YY ZZ
donde XX es el coeciente de t2 (es decir α/2), YY el coeciente de t y ZZ el cociente
constante.
Luego, un ejemplo de código de Matlab para analizar los datos es
plot(t2,phiaju)
alfa=2*coef(1) % El valor buscado de alfa
En el informe indique:
El valor de α medido
El valor de IO calculado.
Recomendaciones
No use una combinación de pesos que produzca un movimiento muy lento pues en ese
caso se torna dicil 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 .
o
Luego, en algunos casos Uds deberán sumar o restar 360 para obtener una curva
continua.
4. Conclusiones [1 punto]
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 identicar dos casos de interés particular:
Ciertamente estos no son los únicos casos concebibles. Podemos agregar movimientos compuestos
de rotación y traslación (como un boomerang en el aire), y de rotación con resbalamiento (como
el neumatico de un vehículo al partir resbalando).
Si pensamos especícamente en una rueda, en la gura de más abajo ilustramos una rueda
dentada que gira, sin deslizar, sobre una supercie igualmente dentada. Cuando 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.
iento o
alam ient
Q resb lam
Q esba
sin sin r
P P
Entre P y Q la velocidad relativa es casi nula Entre P y Q la velocidad relativa es nula
125
Rodadura Sistemas Newtonianos 126
~
dL
~τo = = Io α
~. (4D.1)
dt
Recordar que esta es una relación vectorial. Sin embargo, para cuerpos planos (y si O es escogido
en el plano del papel), entonces los torques resultantes son perpendiculares a este plano. Si k̂ es
un vector unitario saliendo del (o entrando al) papel, entonces podemos expresar
~τo = τo k̂; α
~ = αk̂ ,
τo = Io α .
R
R
T T’ Mg
En la gura de más arriba se ilustran dos casos de rotación en torno a un eje jo: 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 jo
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
el piso está instantáneamente en reposo. Denotamos ese punto por 'C'. Denominando objeto a 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íco,
M~g + f~ + N
~ = M~a ,
M g sin β − f = M ax (4D.3)
Ahora, si aplicamos la ecuación de torques con respecto al punto de contacto, es fácil vericar
τ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.
δx = Rδθ
R δθ
δ s=R δθ
δx δθ
=R ,
δt δt
de modo que al tomar el límite δt → 0,
v = Rω .
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 armar 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 .
Nótese 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. Verique que en tal caso
ax = (5/7)g sin β ≈ 0.71g sin β .
1111
0000
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
H
ESTADO A EVOLUCION A−>B
ESTADO B
El roce: no trabaja puesto que no hay deslizamiento mutuo entre la rueda y el piso, con lo
cual WA→B (roce) = 0.
La energía cinética en cada caso se puede obtener considerando el movimiento de la rueda como
de rotación pura en torno al punto de contacto, es decir
1
K = IC ω 2 .
2
Podemos obtener IC utilizando el teorema de Steiner, IC = M R2 + ICM , con lo cual
2 2 2M gH
ωB = ωA + .
M R2 + ICM
Nuevamente imponemos la condición de rodadura sin resbalamiento, vB = ωB R, vA = ωA R, con
lo cual
2 2 2M gH 2 g sin β
vB = vA + = vA +2 ∆x .
M R2 + ICM 1 + ICM /M R 2
∆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 resul-
tado 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
ω,α
1. Denotamos por ~
R la ubicacion del CM con respecto a P.
V2 ~ .
~a = − ~ ×R
R̂ + α (4D.10)
R
4. Si rotulamos con {1, · · · , i, · · · , N } los puntos de aplicación de las N fuerzas externas sobre
el sólido, y por ~
R al vector que une P con CM, entonces
~ +ρ
~ri = R ~i .
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
~ × (M~a) = IP α
~τCM + R ~.
~ × (M α
~τCM + R ~ = IP α
~ × R) ~.
~τ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 ésta rueda sin resbalar.
x=X+rho*cos(theta)
y=R−rho*sin(theta)
Si tomamos el eje x en la dirección del movimiento del centro de la rueda, y el ángulo θ creciente
en el sentido horario, entonces podemos decir que las coordenadas (X, Y ) del centro del disco
quedan descritas mediante
X = Rθ , Y =R.
Esta construcción garantiza que para θ = 0, el centro de la rueda se ubica en (0, R). Las coorde-
nadas (x, y) de un punto P en la rueda están dadas simplemente por
x = X + ρ cos θ (4D.12)
y = R − ρ sin θ . (4D.13)
con 0 ≤ h ≤ 1. Juegue con la siguiente rutina MatLab, modicándola a gusto, para visualizar lo
que resulta. Al hacer h = 0 se está identicando el centro de la rueda, cuya trayectoria debiera
ser recta. Al hacer h = 1 se identica un punto en la periferia de la rueda. Cualquier valor
intermedio identicará 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 momento
de inercia I) que se muestra en la gura. 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 gura.
Pregunta 5: Dibuje el DCL y escriba Ðsin resolverÐ la ecuaciones dinámicas para el cilindro
(masa M y momento de inercia I) que se muestra en la gura. El cable tira al cilindro desde
su centro.
Pregunta 6: Un cilindro de radio R rueda sin resbalar sobre un plano inclinado en un ángulo
β con respecto a la horizontal. El cilindro, de masa M, tiene un momento de inercia I con
respecto a su eje. Haga un diagrama de cuerpo libre para el cilindro y plantee las ecuaciones
que permitan obtener su aceleración sobre el plano. ( ½no las resuelva!)
Pregunta 7: Explique con detalle el signicado de la condición: RUEDA SIN RESBALAR
Ejercicio 2: Un disco de radio R y masa M rueda sin resbalar sobre una supercie horizontal
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 gura). Haga un gráco 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 supercie
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 vo tal que el resorte se comienza a comprimir.
Calcule la comprensión máxima del resorte y compare con el valor que obten-
dría si en vez del disco se tratara de una partícula de masa M sin roce con la
supercie. 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.
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.
A. Objetivos
1. Visualizar la rodadura de un cuerpo como composición de rotación y traslación.
B1. Materiales:
Tubos de PVC para formar tres anillos de igual diámetro y espesor de paredes.
Plano inclinado.
137
Rodadura Sistemas Newtonianos 138
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; 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.
Apoye la cámara web en la supercie de la mesa y disponga el plano inclinado (en adelante,
la pista) frente a la cámara, a unos 0.5-1.0 m de ella. El eje de la cámara debe ser perpendicular
con el plano sobre el cual ocurre el movimiento sobre la pista. Busque una conguración tal que
se vea la mayor parte de la pista en el video. Si no es posible ver la pista completa, apunte la
cámara tal que se grabe la parte central de la pista
Mida la inclinación de la pista con respecto a la horizontal. Ello se inere de su ángulo con
la vertical, denida por la dirección de una plomada que Ud. puede construir o implementar;
o bien se puede medir con transportador o utilizando ImageJ. La inclinación ideal del plano
probablemente estará en el rango 4-10 .
o
Para elegir la inclinación a usar, es importante probar que los cilindros aceleren suciente-
mente rápido para que sea medible; y que aceleren sucientemente lento de modo de alcanzar a
medir al menos 5 pares (tiempo,posición) en el video grabado (y para que los cilindros aparezcan
nítidos en el video).
½Asegúrese de mantener esta conguración (anillos, plano, cámara) inalterada durante toda
la experiencia!
C2. Registros
Posando uno de los anillos sobre el plano inclinado, en su extremos superior, comience a grabar
una pelí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 derecho 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 (grabe 4 videos 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 aproximadamente cada
3 cuadros con respecto a un punto arbitrario pero jo durante sus mediciones. En cada video
mida un mínimo de 5 pares (cuadro, posición).
Para cada película Ud. tiene una serie de pares (t, x) donde t se mide en [número de cuadros]
x es la distancia a lo largo del plano desde el extremo superior medida en [pixeles]. Ingrese
esos datos a un archivo-m de Matlab,convierta el tiempo de cuadros a segundos, y la posición de
pixeles a metros. La cámara web graba cuadros cada 1/30 segundos. Graque (t, x) y verique que
obtiene aproximadamente una parábola. Ajuste un polinomio de 2do orden (polyfit(t,x,2)),
tal que x(t) = b2 + t2 + b1 t + b0 . Si tiene dudas sobre el comando polyt, escriba help polyfit
en Matlab.
Para uno de los anillos, graque los cuatro ajustes obtenidos con polyt junto con los datos
medidos. Para ello se recomienda dividir la ventana en cuatro para mejor claridad, por ejemplo:
Rotule los grácos con ejes, unidades, fecha y grupo. Adjunte el gráco y el archivo-m
utilizado a su informe.
Relacione el coeciente b2 con la aceleración del anillo (ax ) suponiendo un movimiento uni-
formemente acelerado. De esta forma, para cada anillo, Ud. obtendrá 34 valores experimentales
de ax .
A partir de los valores obtenidos experimentalmente para ax , inera el valor experimental
de Io /M R2 en los tres casos de distribución de masa estudiados. En las situaciones idealizadas
ilustradas en la gura los cuocientes Io /M R2 son 1, 1/2, & << 1 respectivamente.
11111
00000
00000
11111
00000
11111
00000
11111
00000
11111
Periferia Uniforme Centrada
C4. Vericación
Finalmente, realice una carrera entre sus tres anillos. Suéltelos 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 6 veces y anote el orden de llegada. Comente el resultado de las carreras y
explíquelos en términos físicos.
Identique las fuentes de errores que afectan el resultado de la experiencia y comente sobre
como minimizarlos.
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
141
Oscilaciones Sistemas Newtonianos 142
donde hemos utilizado la regla de la mano derecha para expresar el ángulo de forma vectorial;
ds es la longitud del arco de circunferencia entre ambas posiciones, considerado positivo si el
ángulo aumenta. Naturalmente dϕ es adimensional (cuociente entre cantidades con las mismas
unidades), para mayor claridad se dene el radián como una dimensión angular.
Si el movimiento circular es uniforme, entonces por denición la tasa de cambio del ángulo θ
en el tiempo es constante y se denomina velocidad angular:
~
dϕ
ω
~ = = constante, luego
dt
ϕ(t) = ϕo + ω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(ϕo + ωt) ; y(t) = R sen(ϕo + ωt).
R representa entonces la amplitud del movimiento en x o en y . Al ángulo ϕ(t) lo llamamos fase,
y al ángulo inicial ϕo lo llamamos constante de fase. El movimiento es periódico, repitiéndose
cuando la fase cambia en 2π . El período del movimiento lo obtenemos de imponer ωt = 2π :
2π
T = (5A.2)
ω
donde ~xo 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
p
cuyas soluciones son imaginarias. Para simplicar la notación, denimos la constante ωo = k/m
de modo que podemos expresar
b± = ±iωo
Vemos que efectivamente las 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 vericar las siguientes
expresiones:
= B sen(ωo t + ϕs ) (5A.7)
= A cos(ωo t + ϕo ) (5A.8)
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 ϕo . La cantidad ωo = 2π/T la llamamos
frecuencia angular; T es el período del movimiento. Las unidades de la frecuencia angular son
[rad/s]. Notamos que ωo 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 denición del problema.
variable x(t) = A cos(ωo t + ϕo ) que representa posición en el caso de un resorte ideal, pero
puede representar cualquier cantidad física que oscile de forma periódica: ejemplos ángulo
para el movimiento de un péndulo, corriente eléctrica en un circuito de corriente alterna.
Constante de fase ϕo .
Esta es una ecuación muy sencilla, pero al no ser polinomial puede tener un comportamiento
no intuitivo para los alumnos. Les recomendamos gracar 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).
Dada la solución para la posición en función del tiempo x(t) = A cos(ωo t + ϕo ) podemos
encontrar la velocidad derivando c/r al tiempo:
−vo
tg ϕo =
xo ωo
A2 = x2o + (vo /ωo )2
1 1
E = Ue + K = kx2 + mv 2 (5A.9)
2 2
1 2 1
= kA sen2 (ωo t + ϕo ) + mωo2 A2 cos2 (ωo t + ϕo ) (5A.10)
2 2
1 2
= kA (5A.11)
2
La energía mecánica es constante (en la última expresión reemplazamos mωo2 = k ).
Para encontrar una solución para la posición denida 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 torque c/r
dl
= mR2 θ̈
dt
dl
El torque del peso c/r al punto de apoyo del hilo es τg = −Rmg sen θ = dt , luego la ecuación
que describe el movimiento del péndulo se puede expresar como
g
θ̈ + sen θ = 0 (5A.12)
R
Esta ecuación es aparentemente muy similar a la ec. (5A.3) sin embargo el que aparezca
la función sen θ 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 sen θ ≈ θ si θ << 1 rad. Notamos también que las constantes g y R son positivas,
luego podemos denir la constante
g
ωo2 = (5A.13)
R
obteniendo la siguiente ecuación para describir el movimiento del péndulo
θ̈ + ωo2 θ = 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:
donde la amplitud A = θmax es un ángulo. Es muy importante distinguir en este caso la constante
g
ωo = R que representa la frecuencia angular del movimiento de la velocidad angular ω = θ̇ que
para empezar no es constante y tiene un signicado físico distinto. El período del movimiento es
s
2π R
T = = 2π . (5A.16)
ωo g
τg = −dM g sen θ = Iα = I θ̈
dL
donde α es la aceleración angular y la derivada del momento angular c/r al tiempo es
dt = Iα
dado que I es constante. Para pequeñas oscilaciones podemos reemplazar de nuevo sen θ ∼ θ y
obtenemos la ecuación
M gd
θ̈ + θ=0 (5A.17)
I
que describe un movimiento armónico simple con frecuencia angular
M gd
ωo2 = >0
I
s s
I 2L
Tsico = 2π = 2π
mgd 3g
s
L
Tsimple = 2π
2g
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.
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 sec-
ciones 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).
grac
t[s]
Pregunta 2: Determinar el largo que debe tener un péndulo simple sobre la supercie terrestre,
de modo que su período para pequeñas oscilaciones en torno a la posición vertical sea de
60s.
Pregunta 3: 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 4: Determinar el largo que debe tener un péndulo simple sobre la supercie terrestre
de modo que su período sea de 60 s para pequeñas oscilaciones en torno a la vertical local.
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
g
torque obtenemos a = − 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 δ
Solución
La gura ilustra el DCL de la esferita de radio r, las únicas dos fuerzas sobre la esferita son
el peso y la normal. Si calculamos el torque neto c/r al punto de contacto p entre la esferita r y
el cilindro R obtenemos
rmg sin(π − θ) = Ip α
donde el torque apunta hacia afuera del papel si la esferita está al lado derecho del punto más
bajo del cilindro (como en la gura) y hacia adentro del papel si la esferita está al lado izquierdo.
Ahora usamos el hecho que el centro de masa de la esferita se mueve en movimiento circular por
R − r del centro). Luego acm (= rα) representa también la
el interior del cilindro (a una distancia
aceleración tangencial: acm = −(R − r)θ̈ donde los puntos representan derivadas con respecto al
tiempo. Esta es la relación entre la rotación en torno al centro de la esfera y la rotación en torno
al eje del cilindro que es necesario imponer para resolver el problema. Notar el signo negativo: si
la esferita baja desde la posición de la gura entonces α>0 y θ̈ < 0; si la esferita está al lado
izquierdo de la vertical entonces α<0 y θ̈ > 0.
Para pequeñas oscilaciones podemos usar sin θ ∼ θ, llegando a la ecuación
5g
θ̈ + θ=0 ,
7(R − r)
q
5g 2π
que representa un movimiento armónico simple de frecuencia natural ω0 = 7(R−r) = T ; de
donde obtenemos nalmente el período del movimiento como
s
7(R − r)
T = 2π .
5g
1. (1 punto) Un objeto descanza en equilibrio estático desde un resorte colgado del techo
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.
c) (1 punto) Calcule la amplitud de las oscilaciones del resorte una vez que la bolita ha
sido eyectada.
Ejercicio 2:
1. (1 punto) Una persona sentada en un columpio oscila en su frecuencia angular natural ω.
Si la persona se para sobre el mismo columpio, la nueva frecuencia natural del sistema
columpio-persona será:
a) Igual,
b) Mayor,
c) Menor.
c) Expresar la posición x(t) de la masa sujeta al resorte en cada caso suponiendo que el
choque ocurre en t = 0.
1. Determine la ecuación de movimiento para las desviaciones del sistema desde su posición
de equilibrio vertical.
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 identicar la amplitud, fase, frecuencia angular y constante
de fase de un MAS, en particular para el caso de un péndulo físico.
Para llevar a cabo este experimento, utilizaremos de nuevo una cámara web para lmar 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 real-
mente 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 (σ ).
2. Antes de realizar cualquier tipo de medida, es importante conocer el/los instrumentos que
vamos a utilizar y estimar el error probable de las medidas realizadas. Por ejemplo, si
utilizamos una regla graduada, muchas veces la incerteza de una medición de longitud es
similar al intervalo entre las marcas más pequeñas (a menudo 1 mm). Al utilizar un reloj
analógico con segundero el error de cualquier medición es probablemente del orden de 1s.
Al utilizar aparatos más complejos no es tan sencillo estimar a priori el error de cualquier
medida, sin embargo esto siempre es posible de hacer durante el proceso de medición.
¾Con qué precisión pueden medir un ángulo con imageJ? ¾Con qué precisión pueden medir
154
Oscilaciones Sistemas Newtonianos 155
una coordenada? Las respuestas a estas preguntas dependen de la forma de realizar las
mediciones, incluída luminosidad de la sala, la meticulosidad de los alumnos, y el montaje
experimental. Es muy importante que ustedes entiendan y cuantiquen la precisión con la
cual pueden realizar una medición.
3. Para comparar resultados experimetales con resultados teóricos esperados, podemos simple-
mente dibujar los datos medidos y calculados en el mismo gráco. Por otro lado, podemos
determinar uno o más parámetros de una curva teórica de modo que se aproxime lo mejor
posible a los datos medidos (ajuste de parámetros). En cualquier caso la posible conclusión
de si la curva teórica es o no una buena aproximación de las medidas experimentales, de-
penderá de la comparación entre la diferencia entre teoría y observaciones comparada con
la estimación del error de las observaciones. En este caso, el error puede estar dominado
por errores aleatorios (σ ) o por errores sistemáticos que incluyen elementos tales como
aparatos defectuosos o el haber despreciado un fenómeno físico que era relevante para las
mediciones (por ejemplo ¾es importante el roce?).
B. Materiales
El montaje experimental es muy sencillo, sólo deberán cambiar la posición de la cruceta para
realizar las distintas mediciones. Noten que al variar la posición de la cruceta, varía tanto la
posición del centro de masa de la T como su momento de inercia (en este caso en torno al eje de
rotación del extremo superior).
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario vericar el funcionamiento correcto de los ele-
mentos 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.
Para medir el período del péndulo, ustedes grabarán videos con al menos una oscilación
completa del péndulo. Para cada conguración experimental la idea es grabar uno o más
videos y medir el período cuantas veces sea necesario, llegando idealmente a un error menor
a un 10 % en la determinación del período.
Las mediciones serán muy rápidas, ya que simplemente deberán examinar visualmente los
videos cuadro a cuadro para encontrar el número de cuadros que corresponde al período.
Recuerden que el período corresponde al intervalo de tiempo entre dos fases idénticas del
movimiento. La cámara graba 30 cuadros por segundo. En un video ustedes pueden medir
el a período más de una vez, pero jándose de que cada medida sea independiente de
las demás. Por ejemplo, si lmamos un video durante una fase completa del movimiento
(∆t = P ó ∆ϕ = 2π ) pueden medir los intervalos de tiempo entre fase 0 − π/2; π/2 − π ;
π − 3π/2, y 3π/2 − 2π . De ese modo obtienen cuatro medidas independientes de P/4 de
donde pueden obtener un buen valor de P .
Utilizando la T con su longitud máxima, se les pide repetir las medidas del período para
un rango de amplitudes. Elija una serie de 4−6 amplitudes tales como 0.2, 0.4, 0.6,
0.8, 1.0, 1.2 rad. Procure poder establecer con precisión el ángulo inicial de modo que el
experimento sea reproducible. Para cada amplitud mida el período tal como lo hizo en la
experiencia anterior. No es necesario repetir las medidas ya realizadas. Complete la tabla del
informe incluyendo la constante de fase que describe el movimiento. Graque los valores
del período medidos con su error aleatorio en función de la amplitud o con el comando
errorbar en Matlab. Recuerde dibujar sólo los puntos medidos. En el mismo gráco dibuje
el período esperado para un péndulo físico bajo la aproximación de pequeñas oscilaciones.
Para ello debe conocer la aceleración gravitacional local, la posición del centro de masa
y el momento de inercia del péndulo (realmente no es necesario conocer el momento de
inercia, sino el cuociente I/M ). Ustedes ya conocen tanto la posición del centro de masa
como el momento de inercia de la T, en caso de haberlo olvidado, use el teorema de los
ejes paralelos para simplicar 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áco 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]
Consideremos el movimiento del globo en el aire. Los efectos del aire se manifestarán en la
forma de una resistencia del aire sobre el globo, oponiendose al desplazamiento. Estas fuerzas
corresponden a las llamada del tipo roce viscoso.
Esta fuerza que actúa sobre el globo corresponde a una manifestación macroscópica de su
interacción con el gran número de átomos que conforman el volumen de aire que rodea al globo.
El movimiento del globo implica su colisión con dichos átomos y el consiguiente intercambio de
energía y momentum. Dicha complejidad se nos presenta como un problema insoluble, que deja
atrás a las más sosticadas herramientas de cálculo. A pesar de esto, y en realidad gracias a esto,
el efecto neto se puede expresar en forma bastante sencilla. La idea es que la dinámica de todos
los átomos es tan compleja que la única información coherente que se deriva de ella es su efecto
promedio, dejando la información asociada a cada átomo individual como algo indetectable, y en
denitiva, irrelevante. Esta idea, a pesar de su sencillez, puede parecer extremadamente increíble,
casi una trampa. En realidad constituye uno de los principios fundamentales de la física, y en un
modo muy preciso constituye una de las ideas más importantes de la física contemporanea.
Ahora, el efecto neto del movimiento del globo será poner en movimiento un gran número de
átomos, ellos se llevarán energía y el sistema a los ojos del obervador macroscópico parecerá 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
158
Osc. Amortiguadas Sistemas Newtonianos 159
F~ = −f (|~v |)v̂ ,
con −v̂ representando un vector unitario en sentido opuesto al desplazamiento. Con esta expresión
se supone muy poco sobre los mecanismos que determinan f . La forma general de f (v) depende
del regimen de velocidades, densidad del aire, propiedades termodiná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 vericar con gran facilidad, para
f >0 (pero, de otro modo, completamente arbitario) la energía mecánica del sistema es siempre
decreciente. Es decir, esta forma de roce es disipativa.
f (v) → fs (v) ∝ v ,
f (v) → fr (v) ∝ v 2 ,
Fx = −bvx ,
denotando por b al coeciente de roce viscoso (notación Serway). En esta notación es evidente
que el coeciente de roce viscoso tiene unidades de F/v ∼ kg/s. Además, supondremos que el
globo parte (t = 0) con rapidez vo .
11
00 11
00 dx/dt
00
11 00
11
00
11
00
11 00
11
00
11
(t=0) F=− b dx/dt
dvx dvx 1
m = −bvx , ⇒ = − vx ,
dt dt τ
con τ = m/b. Claramente τ tiene dimensiones de tiempo. En esta ecuación uno se plantea la
siguiente pregunta:
¾Cuál es aquella función vx (t), que al derivarla es proporcional a ella misma? La respuesta
la encontramos en la función exponencial, particularmente
vx (t) = A e−t/τ .
Aquí, A representa una constante que está restringida por la condición inicial, vale decir, infor-
mació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áco de esta función se ilustra en la gura 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
t
y(t) = gτ 2 + e−t/τ −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
nalmente
1 − e−2gt/Vt
vy (t) = Vt .
1 + e−2gt/Vt
ẍ + ωo2 x = 0 ,
donde ωo representa la frecuencia natural del sistema. La solución x(t) a esta ecuación es
x(t) = A cos(ωt + ϕo ) .
En ella, las constantes A (amplitud) y ϕo (constante de fase) dependen de las condiciones iniciales
del sistema. En contraste, ωo depende exclusivamente de sus propiedades físicas.
Consideremos esta vez el sistema formado por un bloque de masa m unido a un resorte de
constante elástica k. Como es sabido, en este caso
k
ωo2 = .
m
Además, supongamos la acción de una fuerza de roce viscoso −bẋ actuando sobre el bloque.
k m
O x(t)
b k
mẍ = −kx − bẋ ⇒ ẍ + ẋ + x = 0 .
m m
En forma estándar,
1
ẍ + ẋ + ω 2 x = 0 ,
τ
donde τ representa un tiempo de atenuación. En la asignatura de Mecánica, se estudiará en algún
detalle sistemas que resultan descritos por este tipo de ecuaciones. La solución x(t) en este caso
es
x(t) = A e−t/2τ cos(Ωt + ϕo ) ,
con 2
2 1
Ω = ωo2 − .
2τ
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 ]
Apéndices
5B.4. La fuerza de Stokes
Si consideramos una esfera inmersa en un uido muy viscoso y hacemos que el líquido uya
suavemente, como se ilustra en la gura, entonces la fuerza que el uido ejerce sobre un cuerpo
toma una forma bastante simple. La expresión para tal fuerza fué deducida por Sir George Gabriel
Stokes, en el año 1840. Si la esfera tiene un radio R y la viscosidad del uido es η (de dimensiones
M L−1 T −1 ), entonces la fuerza de arrastre está dada por
FStokes = 6πRη v
donde v es la rapidez del uido. Nótese que la fuerza es proporcional al radio de la esfera, y no
a su área transversal.
1
Fa = ρv 2 Cd A
2
con ρ la densidad del uido, A el área transversal del objeto, y Cd un coeciente de arrastre. El
valor de Cd depende en gran medida de la forma del cuerpo, variando típicamente entre 0.25 y
0.45 para un vehículo.
A modo de estimación, considerando la densidad del aire igual a 1,3 kg/m , un objeto de
3
2
sección transversal de 1 m , a una velocidad de 10 m/s (36 km/h), Cd ∼ 1, entonces la fuerza de
arrastre resulta del orden de
L
(!bL! )
! (T)
m !
(mg) r
Proyectando la ecuación del movimiento exclusivamente según la componente angular φ̂, obten-
emos
1
mLφ̈ = −mg sen φ − bLφ̇ ⇒ φ̈ + ωo2 sen φ + φ̇ = 0 .
τ2
El movimiento resultante es amortiguado pero no armónico: ¾porqué? En el límite de pequeñas
oscilaciones podemos aproximar sen φ ≈ φ, con lo cual
1
φ̈ + ωo2 φ + φ̇ = 0 .
τ2
La ecuación es idéntica a la del resorte con roce viscoso discutido recientemente. En este caso
particular identicamos
g m
ωo2 = , τ= .
L b
5B.8. Ejemplo
5B.8.1. Problema: tiempo de subida y bajada
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 3: Describa en que consiste el M-File en Matlab, como desplegarlo y su uso en esta
Unidad.
Pregunta 4: Indique las precauciones que hay que tener durante el manejo del globo durante
la primera práctica.
1
Fa = ρCd v 2 A ,
2
y explique el contexto en el cual ella aparece en el Material Teórico.
Pregunta 7: Escriba la Fuerza de Stokes y describa el signicado de los términos que la com-
ponen.
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
F = Cd ρ A v 2 .
2
Además suponga la densidad del aire igual a 1,2 kg/m
3 y la del agua igual a 1000 kg/m3 .
2. ¾Podría este paracaidista alcanzar a una lata de soda lanzada antes que él? Considere
en este caso Cd =1,0. Estime la masa de la lata llena y su sección transversal.
1. Estime la fuerza de arrastre del agua sobre el nadador cuando se mueve horizontal-
mente a 2 m/s. Suponga Cd = 0, 7 y el área transversal del orden de 0,1 m .
2
2. Estime la rapidez con que debe mover las manos, con respecto al agua detenida, a
n 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.
Ejercicio 3: Suponga que en su viaje en vehículo consume una cantidad L de combustible para
recorrer una distancia D con rapidez V (100 km/h). Estime el combustible que gastaría
para recorrer la misma distancia a una rapidez λv , con λ=0,8 y 1,2.
Ejercicio 6: ¾En qué unidades queda representado el cociente η/ρ? A esta cantidad se le de-
nomina viscosidad cinemática y es simbolizada por ν . Calcule ν del agua y el aire.
5B.2. Materiales
Un globo e hilo (1 m).
Soporte universal.
170
Oscilaciones Amortiguadas Sistemas Newtonianos 171
Disponga además, en el extremo de la barra, una plomada la cual se usará para calibrar la
vertical.
cordel a n de obtener medidas del ángulo {φi : i = 1, N }, en función del tiempo {ti : i = 1, N }.
La grabación debe contemplar cuatro pasadas por la vertical.
Tabule las mediciones en su cuaderno de apuntes. Ingrese los datos como arreglos en un
programa MatLab ( File > New > M-File) y construya un gráco (x, y) → (t, φ).
Inserte en su programa la función φ(t) dada por
t
φ = A exp(− ) cos(Ω(t − to )) .
2τ
Evalúela para {ti : i = 1, N } y grafíquela en el mismo gráco de los datos. Aquí
2π
Ω= .
T
Por tanteo, busque los valores (A, τ , T , to ) que mejor reproduzcan el comportamiento de los
datos en el intervalo abarcado por los tres primeros nodos. Para encontrar estas constantes
se sugiere:
parta considerando poca atenuación con un τ grande (τ ∼ 50 s), y una amplitud ∼30o .
ejecute su programa Matlab probando la mejor combinación (T ,to ) que reproduzca la posi-
ción de los tres primeros nodos.
Una vez logrado el ajuste de los nodos, proceda a ajustar τ y A de modo que la curva
teórica pase entre los datos.
Una vez obtenidas las constantes (A, τ , T , to ), verique la consistencia entre ellas utilizando
la relación
1
Ω2 = ωo2 − ,
4τ 2
discutida en el Material Teórico de esta unidad. En particular, se espera que con los valores
obtenidos de ω y τ , inera 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 lmaciones. Para ello no transite rápidamente
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
lmada del cordel sea lo sucientemente 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 identicable en la imagen. Cuide de no saturar de luz la cámara.
Deje caer el globo y registre su movimiento. Tabule las mediciones (y del globo en función
del tiempo) en su cuaderno de apuntes. Ingrese los datos como arreglos en un programa MatLab.
Tome el eje y creciente hacia abajo, con su origen en la posición desde la cual se
suelta el globo.
Mediante el uso de Matlab construya un gráco en el cual se representen
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. Para
ello utilice la fórmula dada en el Material Teórico
t
y(t) = gτ 2
+ e−t/τ − 1 .
τ
En este caso utilice el valor de τ obtenido en la primera parte.
A n de incluir en este análisis el caso de una fuerza de roce proporcional a v2, 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 datos), 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
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 inado con un diámetro de unos 15 cm? Verique experimentalmente.
5C.1. Introducción
En las unidades anteriores se estudiaron oscilaciones armónicas mecánicas. Se consideró tam-
bién el efecto del roce viscoso (de algún tipo) sobre estas oscilaciones, lo cual se llama un oscilador
amortiguado. Este oscilador puede representarse grácamente como se muestra en la gura 1,
constando de tres elementos: una masa (inercia), una fuerza restitutiva (resorte) y el roce (amor-
tiguador). Es importante notar que no se ha especicado aún el origen del roce en este sistema.
Puede ser un roce supercial 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 simple-
mente 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 dene 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 signica que
a la imagen anterior agregaremos un forzamiento con una dependencia explícita en el tiempo.
La idea es que si uno fuerza el sistema a una frecuencia particular, el movimiento oscilatorio
será amplicado. Esta frecuencia se denomina 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
forzamiento 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 forzamiento más
que en considerar bien cuál es el mecanismo dominante de roce en estos sistemas.
175
Oscilaciones Forzadas Sistemas Newtonianos 176
k, lo
M
0 x(t)
Figura 19: Esquema de un oscilador amortiguado compuesto por un resorte, una masa y algún
tipo de roce. El forzamiento externo puede ser incluido como una fuerza F~ (t) que actúa sobre la
masa M.
Figura 20: Derrumbe del Hotel Continental ocurrido durante el terremoto de 1985 en cuidad de
México. Este edicio entró en resonancia con la frecuencia de forzamiento del temblor, y sus pisos
superiores colapsaron. Otros edicios contiguos, más altos o más bajos, no sufrieron los mismos
da nos que este hotel.
la gura 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 con-
tinuació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 gura 1 es un comportamiento
universal. Por esto queremos decir que muchos sistemas desarrollan oscilaciones análogas a las
déscritas en la gura 1, independientemente de su naturaleza (sistemas mecánicos, electróni-
cos, electromagnéticos, geológicos, atmosféricos, etc). El simple modelo de la gura 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 continuos y por lo tanto presentan
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 edicios tienen frecuencias naturales de vibración, por muy compleja que sea su estruc-
tura. Si un terremoto excita una frecuencia de resonancia del edicio, este puede llegar a
romperse. Por supuesto que las técnicas actuales reducen esta posibilidad. Un ejemplo de
un colapso real se muestra en la gura 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 forzamiento resonante que produjieron los soldados al marchar.
Se dice que en la actualidad de advierte a los soldados de no marchar sobre un puente.
1
ẍ + ẋ + ωo2 x = 0, (5C.2)
τ
donde τ = M/b representa un tiempo de atenuación, con
p M la masa y b la constante de roce
viscoso (fuerza de roce viscoso = −bẋ), y ωo = k/M es la frecuencia angular de resonancia. La
solución x(t) de esta ecuación es
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 sen(ωt), la ecuación
de Newton resulta
M ẍ = −kx − bẋ + Fo sen(ω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 21: 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)
p = 0. En (a), ω = 7 rad/s y en (b) ω = 9,8 rad/s. La frecuencia
natural de resonancia es ωo = k/M = 10 rad/s.
1 Fo
ẍ + ẋ + ωo2 x = sen(ωt). (5C.6)
τ M
La solución analítica de esta ecuación es mas compleja, pero puede demostrarse que la solución
general es la suma de dos soluciones, la dada por la ecuación (5C.3) más una particular. Por
ahora simplemente daremos la solución:
Fo /M
x(t) = Ae−t/2τ cos(Ωt + ϕo ) + q sen(ωt − δ), (5C.7)
ω 2
(ωo2 − ω 2 )2 + τ
con
ω
tan δ = . (5C.8)
τ (ωo2 − ω2)
Al nal 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
forzamiento 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 forzamiento.
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)
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.
Figura 23: Distintas posibilidades para el incremento en energía mecánica debido al forzamiento.
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.
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 energí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 velocidad 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 pérdida
adicional de energía debido al roce viscoso. Como esta pérdida se incrementa con la rapidez,
Proce = −mv 2 /τ , siempre negativo, la potencia entregada por el forzamiento externo, será even-
tualmente compensada por la disipación. Esto nos permite entender por qué se llega a un régimen
estacionario, donde la partícula oscila con amplitud nita, incluso bajo condiciones de resonan-
cia. Este argumento permite predecir el valor de la amplitud nal. Si la partícula oscila con
amplitud A, la potencia entregada por la fuerza externa es Pext = ωF A sen2 (ωt), mientras que la
2 2 2
potencia disipada es Proce = −mA ω sen (ωt)/τ . De este modo, la condición de compensación
es simplemente: A = (F/M )τ /ω , consistente con el resultado analítico.
El punto de partida entonces es escribir la ecuación de Newton para el centro de masa del
sistema:
d2
M · x + m · (x + y)
(M + m) 2 = −kx − bẋ. (5C.13)
dt M +m
Es importante notar que en este modelo el roce viscoso solo actúa sobre la masa M, y no sobre
la masa pequeña m. Con un poco de álgebra se obtiene
1 m
ẍ + ẋ + ωo2 x = − ÿ, (5C.14)
τ M +m
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 sen(ωt) para todo t, entonces
ÿ = −yo ω 2 sen(ωt). Con esto, la ecuación de Newton resulta
1 myo ω 2
ẍ + ẋ + ωo2 x = sen(ωt), (5C.15)
τ M +m
que básicamente es igual a la ecuación (5C.6) salvo que la constante que va delante desen(ωt),
llamada Fo , ya no es independiente de ω . La solución analítica (5C.7) no cambia, sólo debe
2
cambiarse Fo /M por myo ω /(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 . (5C.16)
2 2 2 ω 2
(ωo − ω ) + τ
y(t)
0' m
k, lo
M
0 x(t)
Figura 24: Esquema de un oscilador amortiguado compuesto por un resorte, una masa y algún
tipo de roce. El forzamiento externo sobre la masa M se obtiene por un movimiento oscilatorio
impuesto a la masa m, es decir con y(t) = yo sen(ω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 25: Curvas de resonancia de la amplitud B como función de ω, en escala lineal (a) y
en escala semilogarítmica (b). La curva segmentada es la misma presentada en la gura 4a. La
curva continua se obtuvo con la ecuación (5C.16) con
p 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.
!
Figura 26: Distintas posibilidades para el incremento en energía mecánica debido al forzamiento.
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.
a. Encuentre la posición del centro de masa del sistema XCM , considerando que la coordenada
del móvil M es x. Luego escriba la ecuación de Newton, o sea:
y ordene la ecuación de modo que en el lado izquierdo hayan sólo términos proporcionales
a x o sus derivadas y en el lado derecho los términos relacionados con y.
b. Compruebe que la función:
1,2ω 2
B(ω) = √ cm
ω 4 − 7ω 2 + 16
5C.9.1. Solución
a. La posición horizontal del centro de M es x(t). La posición horizontal de la bolita m es
x(t) + y(t), luego la posición del centro de masa del sistema es
1 mrω 2
ẍ + ẋ + ω02 x = sen(ωt)
τ M +m
La única dependencia de la solución propuesta x(t) ≡ B(ω) sen(ωt − δ) en el tiempo está
en el numerador luego obtenemos
1 ω
ẋ = B(ω) cos(ωt − δ) y ẍ = −ω 2 x(t) = −ω 2 B(ω) sen(ωt − δ)
τ τ
luego
1 ω
ẍ + ẋ + ω02 x = [−ω 2 + ω02 ]B(ω) sen(ωt − δ) + B(ω) cos(ωt − δ) .
τ τ
Usando sen(ωt−δ) = sen(ωt) cos δ−sen δ cos(ωt), cos(ωt−δ) = cos(ωt) cos δ+sen(ωt) sen δ ;
reemplazando en la ecuación anterior y reagrupando:
1 h ω i hω i
ẍ+ ẋ+ω02 x = (ω02 − ω 2 ) cos δ + sen δ B(ω) sen(ωt)+ cos δ − (ω02 − ω 2 ) sen δ B(ω) cos(ωt)
τ τ τ
ω/τ
factorizando por cos δ y reemplazando tan δ = ω02 −ω 2
se anula el término con cos(ωt) y se
obtiene
ω2
1 2 2 2 1
ẍ + ẋ + ω0 x = (ω0 − ω ) + 2 2 cos δB(ω) sen(ωt)
τ τ (ω0 − ω 2 )
1 1 (ω 2 − ω 2 )
cos δ = √ =q =p 2 0
1 + tan2 δ ω 2 /τ 2 (ω0 − ω 2 )2 + ω 2 /τ 2
1+ (ω02 −ω 2 )2
ω2 mrω 2
1 2 2 2 2 B(ω)
ẍ + ẋ + ω0 x = (ω0 − ω ) + 2 p 2 sen(ωt) = sen(ωt)
τ τ (ω0 − ω 2 )2 + ω 2 /τ 2 M +m
q q
(ω02 − ω 2 )2 + ω 2 /τ 2 = ω 4 + (1/τ 2 − 2ω02 )ω 2 + ω04
2 1 2 √
La amplitud B(ω) es máxima cuando ωmax = ω02 − 2( 2τ ) de donde ωmax = 3,5 ∼ 1,9 ∼
dB(u) 2
ω0 . Alternativamente imponiendo
du = 0 con u = ω se obtiene umax = 32/7 ∼ 4,6 de
donde ωmax ∼ 2,1 ∼ ω0 .
dx d2 x
Nota: ẋ ≡ dt ; ẍ ≡ dt2
Fo /M
x(t) = Ae−t/2τ cos(Ωt + φo ) + q sen(ωt − δ),
ω 2
(ωo2 − ω 2 )2 + τ
P2. Una masa de m, después de caer una distancia h, se adosa a un resorte de constante k. A
partir de ese momento (t = 0) el sistema resultante está descrito por la siguiente ecuación
de movimiento
z̈(t) + 2ω0 ż(t) + ω02 z(t) = C
donde z(t) es la posiciónp
de la masa m (medida hacia abajo desde el punto más alto del
resorte en t = 0), y ω0 = k/m es la frecuencia natural del resorte. Notar que la constante
de amortiguamiento toma un valor crítico = 2ω0 . Para constantes de amortiguamiento
menores que este valor crítico el movimiento corresponde a oscilaciones amortiguadas. En
este caso la solución toma la siguiente forma:
z(t) = (A + Bt)e−ω0 t + D,
donde AyB son constantes que se determinan a partir de las condiciones iniciales (forman
parte de la solución homogénea), y D es una constante que corresponde a la llamada
solución particular (independiente de las condiciones iniciales)
6. ¾Con que frecuencia habría que forzar este sistema para obtener una amplitud de
resonancia máxima?
Reconozca que este sistema se describe por la ecuación de un oscilador amortiguado forzado
Use el método de Verlet para encontrar una solución numérica de la ecuación de Newton
correspondiente.
5C.2. Materiales
Un riel, un carro, dos resortes, un motor y una barra en forma de "L", montado como se
ilustra en la gura.
Lentes de protección.
Matlab.
189
Oscilaciones Forzadas Sistemas Newtonianos 190
ω
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 forzamiento para el carro, de manera análoga al modelo simple expuesto en la guía
teórica. Esta barra se encuentra rmemente acoplada al eje de rotació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 solo modelo de este experimento por mesa
de trabajo, por lo que se debe adoptar un sistema de turnos. La experiencia 1 se realizará 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 realizar en cualquier
órden entre ellas, pero después de la experiencia 1. Es aquí donde deberán implementar un sis-
tema de turnos por grupo. La experiencia 4 debe realizarse al nal, 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 sola 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 un
p
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).
Realice entonces una serie de medidas para al menos 10 valores de ω en torno a ωo ; dicho
de otra manera deberá medir la respuesta del sistema para al menos 10 valores de voltaje
en torno a Vo (±20 %). Recuerde que la respuesta máxima del sistema está dada, según la
p
ecuación (10) de la guía teórica, por ωr = ωo2 − 2(1/2τ )2 . En conclusión, Ω, ωo y ωr son
diferentes, pero para disipación pequeña, ΩE ≈ ωo ≈ ωr .
Se debe medir ω para cada voltaje con el cronómetro y la amplitud de oscilación con la regla
que se encuentra pegada al riel. Recuerde esperar que se amortigue el estado transiente del
sistema pues se quiere obtener la amplitud de oscilación en el estado estacionario. Estime
el error absoluto de su medida de amplitud usando la regla.
Recuerde que para que Matlab ejecute estos programas usted debe elegir como carpeta de
trabajo el directorio donde se encuentran guardados. Antes de usarlos, determine cual es
la masa M y la masa m en este esquema. Se recomienda usar yo = 3 cm. La constante k
(que es la constante total) ya fue determinada en la experiencia 1.
E4a ¾Cómo se comparan cualitativamente los valores de la amplitud del estado estacionario
obtenidos con el algoritmo de Verlet con lo esperado?
Según veremos todas estas ondas presentan propiedades muy similares y tienen una de-
scripción matemática similar. En este curso estudiaremos las ondas que tienen caracterísicas
newtonianas, es decir, cuya dinámica está dada por la ley de Newton y que se propagan en una
dimensión. Más adelante en la carrera se tratará otros tipos de ondas:
193
Ondas Propagativas Sistemas Newtonianos 194
Así una cuerda tensa entre sus extremos tiende a estar recta. Si es deformada, tratará de
volver a su forma original, generando un movimiento que llamaremos ondulatorio.
Al considerar una cuerda tensa, una membrana de tambor o la supercie del agua se observa
que dependiendo de cómo se deforme el sistema se presentan los siguientes fenómenos:
Oscilación colectiva Todo el sistema oscila colectivamente con una amplitud que depende de la
posición, como se observa en una cuerda de guitarra. Este tipo de movimiento se denomina
modo(s) normale(s) y será tratado en la Unidad 6B.
Propagación de pulsos Una perturbación se propaga casi sin deformación, como ocurre en la
supercie del agua, sobre la cual, si se genera una perturbación localizada, ésta se propaga
hacia afuera. De igual manera, si se genera un pulso en una cuerda larga (o una manguera
para regar, por ejemplo), éste se propaga hasta el otro extremo casi sin deformación. Las
ondas de radio también corresponden a este fenómeno.
Ambos tipos de fenómenos se describen de manera unicada como fenómenos ondulatorios, siendo
la propagación de pulsos el que estudiaremos en esta Unidad.
Vista en perspectiva
Tal como se ve en la gura, las varillas pueden girar en torno a su centro, pero, al ser solidarias
al hilo, éste se tuerce ejerciendo torque sobre las varillas vecinas.
Rotulemos con i = 1, . . . , N las varillas y sea θi el ángulo que forma la varilla i-ésima con
la horizontal. Si la varilla que le sigue (la i + 1) se encuentra al mismo ángulo que ésta, el hilo
no está torcido y el torque sobre él es nulo; si la varilla forma un ángulo mayor con la horizonal
(θi+1 > θi ) entonces el hilo está torcido y ejerce torque sobre la varilla i-ésima, el que tiende a
hacer que θi crezca para que los dos ángulos se igualen, es decir un torque positivo (noten que por
acción y reacción el torque que i ejerce sobre i + 1 es negativo); nalmente si θi+1 < θi el torque
será negativo. Considerando la situación en que las diferencias de ángulo no son muy grandes,
se puede considerar entonces como una buena aproximación que el torque sobre i debido a i+1
sea proporcional a la diferencia de ángulos.
τi debido a i+1
= T (θi+1 − θi )
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á dado por la ley de Newton para
sólidos rígidos con un punto jo (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
θi+1 − 2θi + θi−1
2
T (θi+1 − 2θi + θi−1 ) = T ∆ (6A.3)
∆2
2
d θ(x)
= T ∆2 (6A.4)
dx2
d2 θ 2
2d θ
I = T ∆ (6A.5)
dt2 dx2
que se puede escribir como
d2 θ d2 θ
2
= c2 2 (6A.6)
dt dx
donde se ha denido por comodidad
p
c= T ∆2 /I (6A.7)
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 analizado 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á tensada).
Para describir la dinámica de la cuerda, estudiemos lo que pasa en una vecindad de un punto
x de la misma. 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 el trozo de cuerda entre x−∆ y x + ∆. 5 . Al hacer el DCL de ese trozo, las
fuerzas que aparecen son las tensiones por ambos extremos. 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
5
En esta derivación suponemos que ∆ es muy pequeño. El estudiante se preguntará '?pequeño respecto a que?.
Esta es una pregunta dícil de contestar: esencialemente diremos que el cambio en y(x, t) asociado al cambio en
∆ debe ser pequeño, i.e. ∆∂u/∂x y(x, t). Esta discusión es un poco estéril, cuando aprendamos más sobre el
comportamiento ondulatorio sabremos que el criterio realmente útil es que ∆ sea pequeño respecto de la longitud
de onda.
T
θ2
θ1
x−∆ x x+ ∆
La gura muestra una cuerda muy deformada, pero vamos a considerar el caso en que la
deformación es pequeña. Esto signica que 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.
F~ = T (− tg θ1 + tg θ2 )ŷ (6A.13)
dy
tg θ1 = (x − ∆) (6A.14)
dx
dy
tg θ2 = (x + ∆) (6A.15)
dx
Luego, la fuerza es
dy dy
F~ = T (x + ∆) − (x − ∆) ŷ (6A.16)
dx dx
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 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
Esta ecuación tiene la misma forma que la encontrada para el caso de las varillas en torsión,
con la diferencia que la densidad de masa juega el rol de la inercia y la tensión de la cuerda la
fuerza restitutiva.
Nuevamente, si denimos
p
c= T /ρ (6A.21)
d2 y d2 y
= c2 (6A.22)
dt2 dx2
d2 u d2 u
= c2 (6A.23)
dt2 dx2
donde u(x, t) es una variable que describe la deformación relevante en el medio (es decir, u=θ
para el caso de las varillas y u=y para el caso de la cuerda).
Vamos a buscar un tipo de soluciones que se llama de D' Alembert o propagativo. Imaginemos
que tenemos una función f (x) cualquiera (con todas las buenas propiedades de cálculo para poder
calcular las derivadas que sean necesarias). A partir de esta función denimos la función
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 f0 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 cancelación
a ambos lados para cualquier f.
Luego, hemos encontrando una solución de la ecuación de ondas.
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.
es solución (½vericar!).
donde f1 y g1 son funciones arbitrarias. En algunas situaciones resulta más cómodo usar la última
forma.
1. Determine 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 de un pulso depende de la tensión y de la densidad de masa, a través de la
p
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: vf
2 − vi2 = 2a∆. Con
vi = 0, ∆ = y y a = g/2. El tiempo está 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 identicar la dependencia de
la tensión con respecto a la posición. Una vez hecho eso, el resultado es directo si logramos
identicar las ecuaciones resultantes con las de un movimiento acelerado uniformemente. Esta
identicació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 que mientras 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 cor-
responde 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
6A.4.3. Referencias:
Capítulos 14 y 15 deTipler & Mosca.
Capítulo 16 de Serway & Beichman.
Para cada uno de estos ejemplos, indique cuál es la variable que está descrita por la ecuación
de ondas. Es decir, para cada ejemplo indique qué es u en:
d2 u 2
dt2
= c2 ddxu2
d2 u 2
dt2
= c2 ddxu2
P2. [1 pt] Dos cuerdas semi-innitas de densidades distintas están unidas en x = 0. Desde
x>0 viaja hacia x<0 una onda armónica de frecuencia angular ω y longitud de onda λ.
¾Cuáles de las siguientes cantidades cambian al pasar la onda de un medio a otro?
(a) Frecuencia
(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
2
θ(x, t) = A e−(ax+bt) ,
6A.1. Introducción
En el material teórico se estudió la dinámica de las ondas y, en particular, se encontró que la
ecuación de ondas admite una solución particular de la forma
Determinar la relación entre los parámetros físicos de un medio en el cual se propaga una
onda a partir de las oscilaciones del medio.
B. Materiales
Matlab:
Iteraciones
Evaluación de funciones
Gráco de funciones
205
Ondas Propagativas Sistemas Newtonianos 206
C. Experiencias
1. Preliminares [0 puntos]
Verique 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 gracar 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ácos se pide que deduzcan
el signicado de la solución de D' Alembert.
Caso 1
1 cm
f (x) =
(x/5 cm)2 + 1
g(x) = 0
A modo de ejemplo, un programa Matlab que realiza el primer caso con c=1 es
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
Indique en el informe qué se observa y cómo se interpretan las soluciones de D' Alembert.
3. Experiencia 2 [1 puntos] Si
u(x, t) = f (x − ct) + g(x + ct)
1 cm
f (x) = g(x) =
(x/5 cm)2 + 1
Vista en perspectiva
En el laboratorio disponemos de un sistema que tiene varillas de dos largos diferentes, con
momentos de inercia diferentes I1 e I2 .
A partir de lmaciones de las oscilaciones del sistema completo, se busca que Uds deter-
minen 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 obtenido.
Indique el efecto de tener un mayor o menor I en la propagación de las ondas
5. Conclusiones [1 punto]
6B.1. Introducción
Todos los cuerpos exhiben algún grado de exibilidad en cuanto pueden experimentar pe-
queñas deformaciones, sean estas de tipo longitudinal (a lo largo del cuerpo) o transversal (nor-
males 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.
Supondremos el eje x alineado con la cuerda. En este caso la deformación 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
reeren a un elemento innitesimal de la cuerda entre x y x + dx. Al aplicar la 2da ley de Newton
a un elemento innitesimal de la cuerda, obtenemos la ecuación de onda:
∂2y 1 ∂2y
= (6B.1)
∂x2 c2 ∂t2
donde c = ( τρ )1/2 , con τ= tensión de la cuerda y ρ = densidad lineal. Notar que c depende
exclusivamente de las propiedades del medio y no de las condiciones iniciales o amplitud de las
deformaciones.
satisface la ecuación de onda. 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 (es por
eso que el parámetro c se conoce como velocidad de fase).
209
Ondas Estacionarias Sistemas Newtonianos 210
2πx
y(x, 0) = Asen( ) (6B.3)
λ
A representa la amplitud máxima de las deformaciones. Los nodos de la cuerda (y = 0)
2πx
en la condición inicial ocurren en todas las posiciones x que satisfacen
λ = nπ , es decir
x = nλ
2 = {... − 2λ, −3λ −λ λ 3λ
2 , −λ, 2 , 0, 2 , λ, 2 , 2λ, ...}. El gráco siguiente ilustra la forma de la
cuerda en t = 0.
Notar que en x = λ, 2λ, ... 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. Combinando (2) y (3) obtenemos entonces:
π x ct
y(x, t) = Asen(2 (x − ct)) = Asen(2π( − ))
λ λ λ
x t
y(x, t) = Asen(2π( − )) (6B.4)
λ T
λ
donde hemos denido T = c . Interpretemos ahora esta nuevo parámetro. Por simplicidad
consideremos x= 0, de manera que (4) se reduce a y(0, t) = Asen(2π Tt ) cuyo graco se muestra
en la gura a acontinuación.
2π 2π
donde k =
λ = numero de onda y ω = T = frecuencia angular (rad/s). Algunas veces
se emplea la frecuencia f =
1 −1 = Hertz = Hz). Las ecuaciones (4) y (5) suponen que
T (s
y(x = 0, t = 0) = 0, lo cual no siempre es así. Una versión mas general de (4) se escribe como:
donde ϕ se denomina fase. Tanto la amplitud A como la faseϕ dependen de las condiciones
iniciales, es decir las deformaciones impuestas a la cuerda en t = 0.
Las condiciones iniciales también pueden dictar el valor de λ (y entonces k) con lo cual el
valor de T (y entonces ω) queda completamente denido pues:
λ ω τ 1
= = c = ( )2
T k ρ
(Alternativamente, la condición inicial puede dictar el valor de T lo cual ja el valor de λ ).
Las relaciones anteriores indican que la longitud de onda (o número de onda) NO es independi-
ente del período (o frecuencia) cumpliéndose que:
Ondas largas (# de onda pequeño) son ondas de período largo (baja frecuencia)
Ondas cortas (# de onda grande) son ondas de período corto (alta frecuencia)
¾Qué pasa cuando t≥0 ? Para responder esta pregunta emplearemos el principio de super-
posición, suponiendo además que la cuerda se extiende hacia el innito 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 ecuación anterior indica que nuevamente una perturbación viajera se reeja en x = 0 pero
esta vez comienza a retroceder manteniendo su forma completamente (Fig. 30b).
La ecuación (6B.7) 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 que la onda original (sen(kx)) y oscila en el tiempo al igual que la original
ω
(ω o T no cambian). Los valores k y ω aun satisfacen = c. La forma de esta onda se muestra
k
en la gura 31.
λ
Notemos que en los nodos (separados cada
2 ) la amplitud es 0 mientras que en los antinodos
λ
(también separados por ) la amplitud es el doble de la original 2A.
2
nπ 2L
sen(kL) = 0 → kL = nπ → k = → λn =
L n
con n = 1, 2, 3, ...
λ nc
Además, como c= T = λf → fn = 2L
Los pares [λn , fn ] denen los modos normales de la cuerda. A medida que n aumenta, dis-
minuye el largo de la onda y aumenta su frecuencia 32.
4L nc
λn = 2(n−1) y fn = 4L
2L nc
λn = n y fn = 2L
Figura 32: Modos normales en una cuerda con sus dos extremos jos
6B.5.4. Referencias
Capítulo 16 de Tipler & Mosca.
Capítulo 18 de Serway & Beichman.
Pregunta 2: Una onda estacionaria se genera a partir de la superposición de dos ondas viajeras
idénticas en una cuerda tensa. Las ondas viajeras tienen una amplitud A = 0, 005 m y una
frecuencia f = 2 Hz . Determine la amplitud que experimenta la cuerda en los antinodos
de la onda estacionaria y el período con que estos puntos experimentan un ciclo completo
de oscilación.
Pregunta 3: Una onda estacionaria se genera a partir de la superposición de dos ondas viajeras
idénticas en una cuerda tensa. Las ondas viajeras tienen una amplitud A = 0, 005 m y una
longitud de onda λ = 0, 1 m. La velocidad de fase en esta cuerda es c = 20 m/s. Determine
la amplitud que experimenta la cuerda en los antinodos de la onda estacionaria y el período
con que estos puntos experimentan un ciclo completo de oscilación.
Después de un cierto tiempo se vuelve a sacar una fotografía instantánea como se muestra
en la gura 2. El punto A se encuentra a una altura H2 = 0,02 m.
Determine el valor de L (distancia entre nodos), la altura máxima a la cual llegará el punto
A y el tiempo que transcurrirá entre el instante de la fotografía y el momento en que A
alcance su altura máxima.
Ejercicio 2: Se tiene un arreglo como indicado en la gura: un bloque de masa m pasa por una
polea ideal sujeto por una cuerda ideal de densidad lineal de masa ρ = 2 × 10−3 kg/m.
La cuerda está ja a una distancia L=2 m horizontal del eje de la polea a un aparato
que fuerza vibraciones verticales con una frecuencia f ja. Cuando la masa m es 16,0 kg o
25,0 kg se observan ondas estacionarias, y no se producen ondas estacionarias con masas
intermedias.
B. Materiales
Balanza digital.
C. Experimento
219
Ondas Estacionarias Sistemas Newtonianos 220
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 en 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 (en pixeles) y su desviación estándar σl (en
pixeles).
3. Calcular las frecuencias de los tres primeros modos normales para el montaje experi-
mental de la experiencia 2.
7A.1. Introducción
En esta unidad estudiaremos la fuerza que ejerce un conjunto de partículas cuando rebotan en
una supercie. Partiremos con nociones simples acerca de rebotes de partículas y construiremos
un modelo más elaborado para cuanticar 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 que un gas ejerce sobre las paredes de su
contenedor; la fuerza de arrastre del aire sobre un globo cuando se mueve rápidamente; la fuerza
que el viento ejerce sobre un objeto; 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 manguera;
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 denir la densidad de masa, que denotaremos por ρ.
222
Presión Colisional Sistemas Newtonianos 223
De esta forma, la masa δM contenida en un pequeño volumen de tamaño δV estará dada 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 nal, 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
000000000000000000000000000000000000
111111111111111111111111111111111111
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
Una situación interesante ocurre cuando una lluvia de partículas rebota sobre la misma
supercie. Para simplicar la discusión, supondremos que existen n partículas por unidad de
volumen, las cuales se mueven con igual velocidad de incidencia sobre una supercie plana.
A
φ φ
v δt cos φ
φ
v
δt
δ V =Av δt cos φ
En este caso cuanticaremos el momentum neto δ P~ que imprime la supercie a las partículas en
un lapso δt. En la gura se ilustra un chorro de partículas que incide sobre la supercie. En ella
se muestran las partículas que colisionarán la supercie 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
Identicamos entonces la fuerza media responsable 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̂ .
Observamos que la fuerza que la supercie ejerce sobre el chorro es proporcional a:
el cuadrado de la rapidez;
v= 60 km/h≈ 17 m/s, ρ≈ 3
1 kg/m , estimamos F⊥ ≈ 580 N, equivalente a una masa de 58 kg.
Fcol
p≡ = 2ρ v 2 cos2 φ
A
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. En el sistema
2
internacional, una unidad de presión se expresa en N/m , correspondiente a 1 pascal:
1 N/m
2 = 1 pascal = 1 Pa
La presión es una cantidad física muy útil para caracterizar la manifestación mecánica de
un medio uido. 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 dirección).
En el caso del aire, las partículas que golpean un pequeño elemento de supercie provienen de
todas las direcciones dentro de medio hemisferio. Además, las partículas no son monoenergéticas
sino que se distribuyen probabilísticamente. El cálculo de la presión en este caso requiere de
herramientas que van más allá de las con que se cuenta en 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 cuando 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 v2. 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.
7A.4.1. Referencias
Secciones 12.1 a 12.4 de los apuntes de Massman.
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 graque 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 for-
madas por dos placas circunferenciales de aproximadamente 1 m
2 cada una, ubicada a 5
m del eje. Las placas mantienen un ángulo de inclinación entre 25 y 30 grados con la ho-
rizontal. Estime el período de rotación de la hélice a n 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 rebota elásticamente sobre una pared, el cambio de momen-
tum 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 supercie perfectamente
resbaladiza pero con pérdida de energía, caracterizada por
Ejercicio 6: En la gura 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/m ), a razón de 2
3
Determine el incremento de la fuerza normal en la cuña triangular cuando sale agua por la
cañería. Considere el diámetro de la cañería igual a 10 mm.
Para ello se hará caer una carga de garbanzos en forma vertical, se medirá la tensión del cordel que
soporta el plano inclinado y se cuanticará la interacción (fuerza) entre el chorro y la supercie
donde ocurren los impactos. Con esto se espera que:
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;
6
Entendemos por pulso al chorro de extensión longitudinal acotada. La duración de ellos es corta en relación
a las escalas de tiempo características de las cantidades físicas a analizar.
231
Presión Colisional Sistemas Newtonianos 232
• Sensor de fuerzas;
• Transportador;
• Huincha de medir;
• Plomada;
• Tubo (15 cm);
• Un puñado de garbanzos; y
• Programa MATLAB Media.m (disponible en página web del curso)
(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 n 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.
n de evitar la salida de garbanzos. Note que los garbanzos en el piso son un factor
de riesgo.
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 jo entre 35
o 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.
o o o
En la Experiencia 2 se medirán las fuerzas medias para inclinaciones de {35 , 45 , 55 }; 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;
8
B∆=
(ρ A⊥ ∆) g .
3
La cantidad entre paréntesis tiene dimensiones de densidad×volumen=masa. En particular, si
∆ fuese la extensión longitudinal (aproximada) del pulso de garbanzos y ρ fuese más o menos
uniforme, entonces (ρ A⊥ ∆) corresponderá a la masa M de todos los garbanzos. Así entonces,
8M g
∆= .
B
La determinación experimental de B permite inferir la extensión aproximada del pulso de gar-
banzos.
7A.8. El informe
Para que el informe se considere completo este debe incluir
1. Tabulación de los datos de los valores medios de la fuerza en función de la altura o ángulo,
según corresponda. Representar estos datos grácamente. Tanto la tabla como el gráco
puede ser anexados en hojas separadas (corcheteadas al informe). Puede usar papel
milimetrado.
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.
tuvieron la información.
7
Es siempre pertinente citar los autores de trabajos no realizados por uno.
7A.8.1. No olvide:
Rotular debidamente sus grácos;
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 simplicar ideas supongamos que en este intervalo se identican N valores de la fuerza en
instantes uniformemente espaciados (δti = δt), entonces denimos la fuerza media como:
N N N
X X 1 X
Fm δti = Fi δti ⇒ Fm = Fi .
N
i=1 i=1 i=1
El signicado geométrico detrás de esta denició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.
X
δpi = ∆p = pb − pa ,
i
X X
Fi δti = δpi ⇒ Fi δti = δpi = ∆p ,
i i
obtenemos
∆p = Fm ∆t .
P
Aquí hemos identicado ∆t = i δti = (t1 − t0 ) + (t2 − t1 ) + . . . + (tN − tN −1 ) = tb − ta .
J ≡ i Fi δti se le denomina impulso, que en el contexto de la segunda ley
P
A la cantidad 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 n 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 gura 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.
8
Se muestra a continuación la ejecución y salida orrespondiente de Media.m
8
c
>> Media('colisionH45b.txt',250,105)
Durac: 0.210 s
Fondo: 2.421 N
Media: 2.355 N
FzaMedia: 0.066 N
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
denidos. Un líquido en cambio tiene un volumen denido, pero su forma no lo es. Finalmente,
un gas puede ser considerado como un estado que no tiene ni forma ni volumen bien denidos,
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 uido: 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 uidos 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 uidos, su descripción mecánica se divide en la Hidroestática y la Dinámica de
uidos. La primera descripción corresponde, como su nombre lo dice, a la estática de uidos, 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 uido afectan su movimiento. También se puede incluir el cómo el
movimiento de un uido 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.
238
Arquímedes Sistemas Newtonianos 239
Consideremos ahora que pasa con un uido en reposo, es decir un uido en una situación
estática. ¾Qué signica exactamente que el uido 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 uido 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 uido cumple esta condición, entonces el uido completo
está en una condición estática.
7B.2. Presión
Es bien conocido que un uido en reposo ejerce una fuerza sobre las paredes de su contenedor.
Basta pensar en un globo. Al inarlo, é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 inado estando
el gas en 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 desinado). Pero no puede, simplemente porque el
gas que está en su interior ejerce una fuerza sobre sus paredes internas. El globo debe nalmente
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 uido 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
uido homogéneo, es una cantidad constante independiente del espacio.
El hecho que un uido en reposo ejerza 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 uido y las paredes. Como esta fuerza
es para todos efectos prácticos constante, podemos concluir que los átomos o moléculas que
componen el uido tienen una rapidez característica. Notemos que es esta rapidez característica
de los átomos o moléculas de un uido lo que ja su temperatura, T ∼ hv 2 i. Un uido más frío
tiene átomos o moléculas que se mueven con una rapidez típica menor, y viceversa. Además,
los 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 estudiada en la unidad 7A. La
diferencia es que la cantidad v (en la dependencia en v2) ahora no es la velocidad de un chorro
o de un ujo, sino que una rapidez promedio.
|F~ |
P = ,
A
donde F~ es la fuerza que se ejerce sobre una pared de área A. En el sistema internacional de
unidades, la presión tiene como unidad el Pascal, denido como
N
1 Pa = 1 .
m
Figura 33: Representación esquemática de un volumen del uido. El promedio de las velocidades
puede se cero, el uido está en una situación estática. Pero el promedio de la rapidez, no es cero.
Esto dene la temperatura de un uido, 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 34: Representación de un volumen de uido, de altura h, área A, ubicado a una profun-
didad d.
Figura 35: 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 uido una fuerza F2 sobre el pistón.
más grande hay un aumento de presión ∆P , lo que se traduce en una nueva fuerza F2 , la cual
cumple ∆P = F2 /A2 . Finalmente obtenemos
A2
F2 = F1 .
A1
Luego, como A1 < A2 , entonces F2 > F1 . Si la diferencia de áreas es muy grande, entonces la
diferencia de fuerzas también lo es. Esto podria parecer violar alguna ley básica de la Física,
pero no es así. Hay que tener en cuenta que el volumen de uido desplazado es el mismo,
A1 ∆x1 = A2 ∆x2 . Usando esto en la ecuación anterior, obtenemos F1 ∆x1 = F2 ∆x2 , lo cual no
es nada más que la conservación de Energía, pues establece que el trabajo realizado por ambos
pistones es igual.
La magnitud de la fuerza de empuje es igual al peso del uido desplazado por el objeto.
Luego, si un objeto sumergido desplaza una masa M del líquido, el empuje que experimenta 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 uido, 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 uido: ay = (ρf − ρo )g/ρo . Así, un objeto mas denso que el uido
tendrá una aceleración negativa, en el mismo sentido que la gravedad, y se hundirá. En cambio,
un objeto de densidad menor al del uido, tendrá una aceleración positiva y otará.
En el caso que un objeto ote, 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 uido, 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/m
3 y ρf = 1000 kg/m
3 ¾cuál es el
Pregunta 3: Determinar la fracción del volumen bajo la supercie del agua de un pedazo de
hielo que ota en equilibrio. La densidad del hielo es 917 3
kg/m . Justicar.
Pregunta 4: Dos vasos idénticos están llenos de agua hasta el mismo nivel. En uno de los vasos
otan cubos de hielo. ¾Cuál vaso pesa más? Justicar.
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 2
m/s .
245
Principio de Arquímedes Sistemas Newtonianos 246
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 signica 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 sensor
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 yB 0 ≡ B − mG g . (3)
0
Determinados A y B , obtenga B y
mG con las dos ecuaciones para F− y F+ .
Nota: La regresión lineal se debe hacer con la función polyt de Matlab, explicada en
la guía práctica de la unidad 4C. Si tiene dudas de su utilización escriba help polyt
en la línea de comando de Matlab.
Figura 36: 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 calibre, 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 Po , sin que
estos estén sumergidos bajo el agua. Compruebe que Po = Mo g , donde Mo es la masa
del objeto.
5. Con los cuerpos completamente sumergidos en agua, repita las mediciones del peso
de cada objeto (P1 ).
6. Determine la fuerza de empuje E que experimenta cada objeto a partir de las medidas
de Po ± ∆Po P1 ± ∆P1 . Compare sus mediciones de E con la predicción del Principio
y
de Arquímedes, es decir E = ρf V g .
Guía de Prácticas
248
Métodos Numéricos Sistemas Newtonianos 249
1. Conocer las capacidades de los métodos numéricos en la solución y análisis de los sistemas
newtonianos.
B. Materiales: Matlab
C. Experiencias
L
Experiencia 1.- Una partícula de masa m puede deslizar sin roce
por un riel horizontal. La partícula está unida a dos resortes V
idénticos de constante elástica k y largo natural L. Los ex-
tremos de los resortes están jos a una distancia L y 2L del
riel, tal como se indica en la gura. Inicialmente se le da a la
partícula una velocidad V hacia la derecha. Se pide determinar
el alcance de la partícula si L = 0, 1m, m = 0, 1kg, k = 2N/m, 2L
para los siguientes valores de V : 0,3 m/s, 0,5 m/s, 1 m/s.
Velocidad V Alcance
ϕ̈ = −(g/L) sen(ϕ)
π
Considere que el péndulo se suelta desde el reposo en un ángulo inicial ϕo = 4
Use g= 9, 8m/s2 , L = 0, 5m, ∆t = 0, 01s, 0, 05s, 0, 5s Integre la ecuación hasta T = 10s.
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 ϕo = 10 , 4 , 2
Complete la siguiente tabla:
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 sigue el período al aumentar
el ángulo inicial.
¾Cómo gracar 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 .
¾Cómo hacer un gráco de datos con barras de error?
Utilice errorbar(x,y,err,'*'), con err el arreglo de errores asociado a y, en las coorde-
nadas x.
¾Cómo incluir leyenda de los símbolos usados en un gráco?
Utilice legend('leyenda 1, leyenda 2, ... ') después de la última instrucción plot.
El orden de las leyendas sigue el orden de las curvas desplegadas.
F. Lecturas recomendadas
Manual de Matlab
Aprender a usar el programa Measurement & Automation para vercar que el sensor de
fuerza está funcionando correctamente
B. Materiales
C. Experiencias
Ahora verique que el sensor de fuerza que está conectado con el canal de entrada análogo
número0 (AI0) está trabajando adecuadamente. Para ello lance la aplicación Measurement
& Automation (programa de medida y automatización), cuyo ícono se encuentra en el
escritorio de su PC.
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 Conguration: 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 adquisició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
Con Analog Output y Analog Input se conguran la salida AO0 y entrada AI0 de una manera
muy similar a lo que se hizo en la experiencia 2.
Con la botonera Save to ASCII puede elegir un nombre de archivo como también un
directorio donde guardarlo. Tenga cuidado de mantener el formato de archivo con ASCII
(formato texto), pues así podrá leer facilmemte los datos con Matlab. Verique que su
conguración de Windows graba los números con punto decimal y no coma.
Ahora pruebe la medida de tensión de corte de un hilo, lanzando la medición con el botón
Run Once. De hecho, esto lanza las tres tareas mencionadas en forma consecutiva.
El hilo debe atarse un extremo al sensor de fuerza y el otro extremo pasarlo por la polea y
anudarlo al dedo de la persona que tirará rme, pero suavemente incrementando la mag-
nitud de la tensión aplicada. Se puede visualizar la adquisición seleccionando la lengueta
Data View y simplemente arrastrando el cursor desde Analog Input hacia esta ventana.
Puede también leer estos datos desde Matlab como se indica en el texto de material teórico
sobre Métodos Experimentales.
Con esta prueba determine el rango en el cual usará el sensor, ±10 N ó ±50 N. Explique
la elección del rango:
Con los datos obtenidos por su grupo realice un histograma que muestre la distribución de
ocurrencia de la tensión de corte del hilo.
Con los datos obtenidos por su mesa realice un histograma que muestre la distribución de
ocurrencia de la tensión de corte del hilo.
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.
B. Materiales
C. Experiencias
A Calcular numéricamente las coordenadas (Xcm , Ycm ) del centro de masa de un arco
cerrado (α = 2π ) usando M =1.AA kg y R=1.BB m, con AA y BB los dos últimos
dígitos del RUT de integrantes distintos del grupo. Estudie el caso N = 100.
B Ahora se desea calcular el centro de masa de medio aro, es decir cuando α = π . Calcule
y tabule las coordenads del centro de masas para los casos θi = 0, π/4, π/2, 3π/4, π .
p
Reporte 2 + Y 2 /R
Xcm e interprete. En todos estos casos use N = 100.
cm
Experiencia 2: Considere una barra uniforme de masa M y longitud L que gira con velocidad
angular ω en torno a un eje perpendicular. El eje se ubica en la barra, a una distancia λL
(λ ∈ [0 − 1]) de uno de sus extremos. Por lo visto en los apuntes y por análisis 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 coeciente adimensional que depende de los valores de λ, vale decir, γ = γ(λ).
Utilizando la misma estrategia de discretización de la experiencia 1, calcule y graque γ
para los valores λ=0, 1/4, 1/2, 3/4 y 1. Use N = 100. Incluya en su gráco los resultados
para λ = 0, ZZ , con ZZ los dos últimos dígitos de todos los integrantes del grupo.
D. Análisis e informe
E. Lecturas recomendadas
Para esta experiencia cada grupo contará con una barra sujeta en un extremo por un tornillo
jo a un soporte vertical. La barra se sostendrá en posición horizontal por un hilo conectado al
sensor de fuerzas que utilizaron hace dos semanas. Por último, se colgarán masas bajo distintas
posiciones de la barra. Cada grupo medirá entonces la tensión de la cuerda bajo una variedad de
conguraciones.
En general, mediciones como las que harán esta semana se utilizan con una variedad de
objetivos en física o ingeniería. A veces se busca encontrar algun tipo de ecuación o ley de la
naturaleza que reproduzca las mediciones realizadas. Por ejemplo podemos buscar si la fuerza
ejercida por un medio elástico es proporcional al estiramiento desde su posición de equilibrio, ¾o
quizás al cuadrado del mismo? También nos puede interesar medir en qué rango de estiramientos
es válida la aproximación matemática encontrada, quizás debido a deformaciones permanentes o
incluso rompimiento del medio elástico en estudio. Otras veces se busca medir alguna propiedad
de los cuerpos tal como podría ser la aceleración gravitacional en el laboratorio, o la tensión de
corte de un hilo como en la Unidad 2 de este curso. En estos últimos casos se busca encontrar
el mejor parámetro que reproduce las observaciones, esto se denomina ajustar una función a los
datos y lo llevarán a cabo más adelante.
El objetivo de esta semana es más sencillo, y simplemente les pedimos comparar las medidas
realizadas por ustedes con los valores esperados de acuerdo a las ecuaciones de la estática. Es
decir, ustedes deberán comparar la medida de la tensión del hilo realizada con el sensor de
fuerzas con la tensión calculada a partir de las leyes de la física (las leyes de estática en este
caso). La conclusión más fundamental a la que esperamos ustedes lleguen es si ¾las leyes de la
estática reproducen bien la tensión medida bajo una variedad de condiciones o no? Sin embargo,
no esperamos simplemente una respuesta positiva o negativa sino una comparación cuantitativa
(diferencia absoluta o relativa) entre la tensión medida y la tensión calculada. En caso de que
la diferencia encontrada sea superior a lo que ustedes sospechan debiera darles, busquen alguna
explicación plausible para tal efecto.
Es sumamente importante que comprendan que en caso de haber diferencias entre los valores
medidos y los valores calculados, éstas se pueden deber a una gran cantidad de efectos sistemáticos
o aleatorios que pueden ser distintos para cada grupo. Esto hace que las conclusiones a las que
llegue cada grupo no deben ser necesariamente las mismas.
B. Materiales
El montaje experimental es similar al ilustrado en la gura 37, 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 conguración de equilibrio.
C. Experiencias
1. Preliminares
Al inicio de cada experimento es necesario vericar el funcionamiento correcto de los el-
ementos 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é in-
stalado correctamente. Recuerden que los cables del sensor deben estar conectados a la
tarjeta de medición en las siguientes puertas:
Luego Abrir Test Panel e ingrese 5V bajo Analog Output y presione Update. Final-
mente seleccionar Analog Input y vericar que el sensor da una medida cercana a 2,5V
sin fuerza aplicada.
b) Cerrar la aplicación Measurement and Automation y abrir Signal Express para medir
señales de fuerza. Modicaremos el programa Tensión de Corte que se utilizó hace
dos semanas. Hoy mediremos tensión bajo condiciones estáticas, luego no buscaremos
encontrar una tensión mínima o máxima en un intervalo de tiempo, pero sí obtener
la tensión con la mayor precisión posible. Para ello realizaremos una gran cantidad de
mediciones con Signal Express.
Para simplicar el procedimiento no grabaremos los datos obtenidos con Signal Ex-
press sino simplemente obtendremos el valor medio y la desviación estándar de una
gran cantidad de medidas de la tensión.
En esta etapa cada grupo deberá decidir si utiliza el sensor de fuerzas en el rango de
validez de 10 N o de 50 N.
La tensión del hilo y su error se obtienen de la estadística que entrega el programa Signal
Express. Recuerden que para encontrar la tensión T a partir del voltaje promedio medido
V se utilizan las constantes A y B de cada sensor de fuerzas:
T = A ∗ V + B[N ]
Para encontrar la masa de la barra deberá primero hacer un diagrama de cuerpo libre y
resolver las ecuaciones de la estática. ¾Con cuántas cifras signicativas conocen la masa de
la barra? ¾Con respecto a qué punto conviene calcular el torque y porque?
¾Tiene sentido la masa encontrada, es decir, está en un rango más o menos esperado? Si
obtienen valores claramente erróneos tales como 1 mg o 103 kg es casi seguro que tienen
un problema de unidades en sus ecuaciones.
Tomar una masa y medir la tensión de la cuerda al variar la distancia de dicha masa con
respecto al eje del soporte. Medir la tensión en cinco o más posiciones distintas de la carga.
Gracar la tensión en función de la posición junto con su barra de error.
Al hacer un gráco lo más común es unir los puntos con una línea continua, dando la
apariencia de haber tomado un gran número de mediciones. En esta experiencia buscaremos
gracar exclusivamente la tensión medida con su error. En un macro de Matlab (archivo
.m) dena arreglos con las cantidades a gracar, por ejemplo:
errorbar(pos,tension_exp,error_tension,'+').
Noten que hemos utilizado un simbolo + para indicarle a matlab que ponga dicho símbolo
en cada punto gracado y no los una con una línea. Podríamos haber utilizado algún otro
signo. El comando errorbar graca 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 ), denan 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 esperada del hilo 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 programas.
Gracar la tensión calculada como una línea continua junto al gráco anterior.
Ejemplo de gráco
errorbar(pos,tension_exp,error_tension,'+')
hold on
plot(pos_calc,tension_calc)
Adjuntar al informe una copia impresa del gráco bien rotulado con unidades en ambos
ejes, fecha y grupo.
de enganche del hilo sobre la barra. Con una masa ja a una distancia también ja del eje,
deben medir la tensión para al menos tres posiciones distintas del anclaje.
Para cada posición, midan también el ángulo entre el hilo y la barra, ¾con que precisión
creen que pueden medir dicho ángulo? Gracar la tensión en función del ángulo, de nuevo
se les pide gracar exclusivamente los puntos medidos con el comando errorbar de Matlab.
En el mismo gráco dibujar el resultado esperado teóricamente suponiendo que conocen
todos los parámetros experimentales de manera exacta.
5. Conclusiones [1 punto]
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 jo en torno al cual es sistema podrá girar libremente. Ojo que
el mismo sistema lo volverán a utilizar en tres semanas más.
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.
B. Materiales
El montaje experimental es muy sencillo, sólo deberán cambiar la posición de la cruceta para
realizar las distintas mediciones. Noten que al variar la posición de la cruceta, varía tanto la
posición del centro de masa de la T, como su momento de inercia (en este caso en torno al eje
de rotación del extremo superior).
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario vericar el funcionamiento correcto de los el-
ementos 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?
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 denida 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 estimadores
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 conguració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.
Para ello haga un gráco 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áco la curva teórica ωcalc vs x con una línea contínua.
Ejemplo de gráco:
errorbar(x_exp,omega_exp,err_omega,'+')
hold on
plot(x_calc,omega_calc)
Imprimir y Adjuntar al informe una copia impresa del gráco bien rotulado con
unidades en ambos ejes, fecha y grupo.
3. Conclusiones [1 punto]
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.
Identicar 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.
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 vericar el funcionamiento correcto de los el-
ementos 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. Verique que es así, para lo cual hágala
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.
En el informe indique las combinaciones probadas, el sentido de giro que se observó (vector
hacia afuera o adentro o ausencia de movimiento) y el signo predicho por la ecuación de
movimiento de la polea. Compare.
Nota: Para que las cuerdas no deslicen por la polea (y así garantizar que se cumplen las
relaciones cinemáticas usadas en el análisis del problema), éstas se deben enrollar con varias
vueltas alrededor de la polea. No haga un lazo pues normalmente quedan ojos y la cuerda
desliza.
Se abre ImageJ y en el menu Import se lee la película como Avi. Hecho esto, es
posible analizar la película cuadro a cuadro.
sale el número de cuadro, que se convierte a tiempo recordando que cada cuadro se
toma cada 1/30 s.
Con esto se obtiene una serie de ángulos para diferentes tiempos {ti , φi }N
i=1 . Para que
la medición sea representativa se deben obtener unos diez valores.
Se sabe que el movimiento de la polea en estas condiciones (ver Material Teórico) ocurre
con aceleración angular constante α. Luego, el ángulo evoluciona en el tiempo como un
movimiento uniformemente acelerado
1
φ = φo + ωo (t − to ) + α(t − to )2 (4C.1)
2
1 2
= A + Bt + αt (4C.2)
2
= A + Bt + Ct2 (4C.3)
donde se han absorbido los diferentes términos que tienen to en las constantes A y B y se
ha denido C = α/2. En la experiencia ωo = 0, pero φo y to serán típicamente diferentes
de cero, de manera que en ecuación (4C.3) los coecientes A y B serán diferentes de cero.
Se pide a Matlab que encuentre la mejor parábola (el mejor polinomio cuadrático)
que pase lo más cercano por los puntos. Este procedimiento se llama ajuste o t y se
hace en Matlab con la función polyfit que entrega los coecientes del polinomio
coef=polyfit(t,phi,2)
coef =
XX YY ZZ
donde XX es el coeciente de t2 (es decir α/2), YY el coeciente de t y ZZ el cociente
constante.
Luego, un ejemplo de código de Matlab para analizar los datos es
plot(t2,phiaju)
alfa=2*coef(1) % El valor buscado de alfa
En el informe indique:
El valor de α medido
El valor de IO calculado.
Recomendaciones
No use una combinación de pesos que produzca un movimiento muy lento pues en ese
caso se torna dicil 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 .
o
Luego, en algunos casos Uds deberán sumar o restar 360 para obtener una curva
continua.
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; 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.
Apoye la cámara web en la supercie de la mesa y disponga el plano inclinado (en adelante,
la pista) frente a la cámara, a unos 0.5-1.0 m de ella. El eje de la cámara debe ser perpendicular
con el plano sobre el cual ocurre el movimiento sobre la pista. Busque una conguración tal que
se vea la mayor parte de la pista en el video. Si no es posible ver la pista completa, apunte la
cámara tal que se grabe la parte central de la pista
Mida la inclinación de la pista con respecto a la horizontal. Ello se inere de su ángulo con
la vertical, denida por la dirección de una plomada que Ud. puede construir o implementar;
o bien se puede medir con transportador o utilizando ImageJ. La inclinación ideal del plano
probablemente estará en el rango 4-10 .
o
Para elegir la inclinación a usar, es importante probar que los cilindros aceleren suciente-
mente rápido para que sea medible; y que aceleren sucientemente lento de modo de alcanzar a
medir al menos 5 pares (tiempo,posición) en el video grabado (y para que los cilindros aparezcan
nítidos en el video).
½Asegúrese de mantener esta conguración (anillos, plano, cámara) inalterada durante toda
la experiencia!
C2. Registros
Posando uno de los anillos sobre el plano inclinado, en su extremos superior, comience a grabar
una pelí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 derecho 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 (grabe 4 videos 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 aproximadamente cada
3 cuadros con respecto a un punto arbitrario pero jo durante sus mediciones. En cada video
mida un mínimo de 5 pares (cuadro, posición).
Para cada película Ud. tiene una serie de pares (t, x) donde t se mide en [número de cuadros]
x es la distancia a lo largo del plano desde el extremo superior medida en [pixeles]. Ingrese
esos datos a un archivo-m de Matlab,convierta el tiempo de cuadros a segundos, y la posición de
pixeles a metros. La cámara web graba cuadros cada 1/30 segundos. Graque (t, x) y verique que
obtiene aproximadamente una parábola. Ajuste un polinomio de 2do orden (polyfit(t,x,2)),
tal que x(t) = b2 + t2 + b1 t + b0 . Si tiene dudas sobre el comando polyt, escriba help polyfit
en Matlab.
Para uno de los anillos, graque los cuatro ajustes obtenidos con polyt junto con los datos
medidos. Para ello se recomienda dividir la ventana en cuatro para mejor claridad, por ejemplo:
Rotule los grácos con ejes, unidades, fecha y grupo. Adjunte el gráco y el archivo-m
utilizado a su informe.
Relacione el coeciente b2 con la aceleración del anillo (ax ) suponiendo un movimiento uni-
formemente acelerado. De esta forma, para cada anillo, Ud. obtendrá 34 valores experimentales
de ax .
A partir de los valores obtenidos experimentalmente para ax , inera el valor experimental
de Io /M R2 en los tres casos de distribución de masa estudiados. En las situaciones idealizadas
ilustradas en la gura los cuocientes Io /M R2 son 1, 1/2, & << 1 respectivamente.
11111
00000
00000
11111
00000
11111
00000
11111
00000
11111
Periferia Uniforme Centrada
C4. Vericación
Finalmente, realice una carrera entre sus tres anillos. Suéltelos 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 6 veces y anote el orden de llegada. Comente el resultado de las carreras y
explíquelos en términos físicos.
Identique las fuentes de errores que afectan el resultado de la experiencia y comente sobre
como minimizarlos.
Para llevar a cabo este experimento, utilizaremos de nuevo una cámara web para lmar 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 real-
mente 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 (σ ).
2. Antes de realizar cualquier tipo de medida, es importante conocer el/los instrumentos que
vamos a utilizar y estimar el error probable de las medidas realizadas. Por ejemplo, si
utilizamos una regla graduada, muchas veces la incerteza de una medición de longitud es
similar al intervalo entre las marcas más pequeñas (a menudo 1 mm). Al utilizar un reloj
analógico con segundero el error de cualquier medición es probablemente del orden de 1s.
Al utilizar aparatos más complejos no es tan sencillo estimar a priori el error de cualquier
medida, sin embargo esto siempre es posible de hacer durante el proceso de medición.
¾Con qué precisión pueden medir un ángulo con imageJ? ¾Con qué precisión pueden medir
una coordenada? Las respuestas a estas preguntas dependen de la forma de realizar las
mediciones, incluída luminosidad de la sala, la meticulosidad de los alumnos, y el montaje
experimental. Es muy importante que ustedes entiendan y cuantiquen la precisión con la
cual pueden realizar una medición.
3. Para comparar resultados experimetales con resultados teóricos esperados, podemos simple-
mente dibujar los datos medidos y calculados en el mismo gráco. Por otro lado, podemos
determinar uno o más parámetros de una curva teórica de modo que se aproxime lo mejor
posible a los datos medidos (ajuste de parámetros). En cualquier caso la posible conclusión
de si la curva teórica es o no una buena aproximación de las medidas experimentales, de-
penderá de la comparación entre la diferencia entre teoría y observaciones comparada con
la estimación del error de las observaciones. En este caso, el error puede estar dominado
por errores aleatorios (σ ) o por errores sistemáticos que incluyen elementos tales como
aparatos defectuosos o el haber despreciado un fenómeno físico que era relevante para las
mediciones (por ejemplo ¾es importante el roce?).
B. Materiales
El montaje experimental es muy sencillo, sólo deberán cambiar la posición de la cruceta para
realizar las distintas mediciones. Noten que al variar la posición de la cruceta, varía tanto la
posición del centro de masa de la T como su momento de inercia (en este caso en torno al eje de
rotación del extremo superior).
C. Experiencias
1. Preliminares [1 punto]
Al inicio de cada experimento es necesario vericar el funcionamiento correcto de los ele-
mentos 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.
Para medir el período del péndulo, ustedes grabarán videos con al menos una oscilación
completa del péndulo. Para cada conguración experimental la idea es grabar uno o más
videos y medir el período cuantas veces sea necesario, llegando idealmente a un error menor
a un 10 % en la determinación del período.
Las mediciones serán muy rápidas, ya que simplemente deberán examinar visualmente los
videos cuadro a cuadro para encontrar el número de cuadros que corresponde al período.
Recuerden que el período corresponde al intervalo de tiempo entre dos fases idénticas del
movimiento. La cámara graba 30 cuadros por segundo. En un video ustedes pueden medir
el a período más de una vez, pero jándose de que cada medida sea independiente de
las demás. Por ejemplo, si lmamos un video durante una fase completa del movimiento
(∆t = P ó ∆ϕ = 2π ) pueden medir los intervalos de tiempo entre fase 0 − π/2; π/2 − π ;
π − 3π/2, y 3π/2 − 2π . De ese modo obtienen cuatro medidas independientes de P/4 de
donde pueden obtener un buen valor de P .
Utilizando la T con su longitud máxima, se les pide repetir las medidas del período para
un rango de amplitudes. Elija una serie de 4−6 amplitudes tales como 0.2, 0.4, 0.6,
0.8, 1.0, 1.2 rad. Procure poder establecer con precisión el ángulo inicial de modo que el
experimento sea reproducible. Para cada amplitud mida el período tal como lo hizo en la
experiencia anterior. No es necesario repetir las medidas ya realizadas. Complete la tabla del
informe incluyendo la constante de fase que describe el movimiento. Graque los valores
del período medidos con su error aleatorio en función de la amplitud o con el comando
errorbar en Matlab. Recuerde dibujar sólo los puntos medidos. En el mismo gráco dibuje
el período esperado para un péndulo físico bajo la aproximación de pequeñas oscilaciones.
Para ello debe conocer la aceleración gravitacional local, la posición del centro de masa
y el momento de inercia del péndulo (realmente no es necesario conocer el momento de
inercia, sino el cuociente I/M ). Ustedes ya conocen tanto la posición del centro de masa
como el momento de inercia de la T, en caso de haberlo olvidado, use el teorema de los
ejes paralelos para simplicar 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áco 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.
Disponga además, en el extremo de la barra, una plomada la cual se usará para calibrar la
vertical.
cordel a n de obtener medidas del ángulo {φi : i = 1, N }, en función del tiempo {ti : i = 1, N }.
La grabación debe contemplar cuatro pasadas por la vertical.
Tabule las mediciones en su cuaderno de apuntes. Ingrese los datos como arreglos en un
programa MatLab ( File > New > M-File) y construya un gráco (x, y) → (t, φ).
Inserte en su programa la función φ(t) dada por
t
φ = A exp(− ) cos(Ω(t − to )) .
2τ
Evalúela para {ti : i = 1, N } y grafíquela en el mismo gráco de los datos. Aquí
2π
Ω= .
T
Por tanteo, busque los valores (A, τ , T , to ) que mejor reproduzcan el comportamiento de los
datos en el intervalo abarcado por los tres primeros nodos. Para encontrar estas constantes
se sugiere:
parta considerando poca atenuación con un τ grande (τ ∼ 50 s), y una amplitud ∼30o .
ejecute su programa Matlab probando la mejor combinación (T ,to ) que reproduzca la posi-
ción de los tres primeros nodos.
Una vez logrado el ajuste de los nodos, proceda a ajustar τ y A de modo que la curva
teórica pase entre los datos.
Una vez obtenidas las constantes (A, τ , T , to ), verique la consistencia entre ellas utilizando
la relación
1
Ω2 = ωo2 − ,
4τ 2
discutida en el Material Teórico de esta unidad. En particular, se espera que con los valores
obtenidos de ω y τ , inera 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 lmaciones. Para ello no transite rápidamente
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
lmada del cordel sea lo sucientemente 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 identicable en la imagen. Cuide de no saturar de luz la cámara.
Deje caer el globo y registre su movimiento. Tabule las mediciones (y del globo en función
del tiempo) en su cuaderno de apuntes. Ingrese los datos como arreglos en un programa MatLab.
Tome el eje y creciente hacia abajo, con su origen en la posición desde la cual se
suelta el globo.
Mediante el uso de Matlab construya un gráco en el cual se representen
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. Para
ello utilice la fórmula dada en el Material Teórico
t
y(t) = gτ 2
+ e−t/τ − 1 .
τ
A n de incluir en este análisis el caso de una fuerza de roce proporcional a v2, 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 datos), 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
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 inado con un diámetro de unos 15 cm? Verique experimentalmente.
Reconozca que este sistema se describe por la ecuación de un oscilador amortiguado forzado
Use el método de Verlet para encontrar una solución numérica de la ecuación de Newton
correspondiente.
5C.2. Materiales
Un riel, un carro, dos resortes, un motor y una barra en forma de "L", montado como se
ilustra en la gura.
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 forzamiento para el carro, de manera análoga al modelo simple expuesto en la guía
teórica. Esta barra se encuentra rmemente acoplada al eje de rotació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 solo modelo de este experimento por mesa
de trabajo, por lo que se debe adoptar un sistema de turnos. La experiencia 1 se realizará 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 realizar en cualquier
órden entre ellas, pero después de la experiencia 1. Es aquí donde deberán implementar un sis-
tema de turnos por grupo. La experiencia 4 debe realizarse al nal, 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 sola 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 un
p
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).
Realice entonces una serie de medidas para al menos 10 valores de ω en torno a ωo ; dicho
de otra manera deberá medir la respuesta del sistema para al menos 10 valores de voltaje
en torno a Vo (±20 %). Recuerde que la respuesta máxima del sistema está dada, según la
p
ecuación (10) de la guía teórica, por ωr = ωo2 − 2(1/2τ )2 . En conclusión, Ω, ωo y ωr son
diferentes, pero para disipación pequeña, ΩE ≈ ωo ≈ ωr .
Se debe medir ω para cada voltaje con el cronómetro y la amplitud de oscilación con la regla
que se encuentra pegada al riel. Recuerde esperar que se amortigue el estado transiente del
sistema pues se quiere obtener la amplitud de oscilación en el estado estacionario. Estime
el error absoluto de su medida de amplitud usando la regla.
Recuerde que para que Matlab ejecute estos programas usted debe elegir como carpeta de
trabajo el directorio donde se encuentran guardados. Antes de usarlos, determine cual es
la masa M y la masa m en este esquema. Se recomienda usar yo = 3 cm. La constante k
(que es la constante total) ya fue determinada en la experiencia 1.
E4a ¾Cómo se comparan cualitativamente los valores de la amplitud del estado estacionario
obtenidos con el algoritmo de Verlet con lo esperado?
En esta práctica se buscará, mediante análisis numérico, determinar el signicado físico de esta
solución. Para eso se probarán distintas funciones f y g , con distintos valores de c, lo que permitirá
comprender el signicado de esta solución.
Determinar la relación entre los parámetros físicos de un medio en el cual se propaga una
onda a partir de las oscilaciones del medio.
B. Materiales
Matlab:
Iteraciones
Evaluación de funciones
Gráco de funciones
C. Experiencias
1. Preliminares [0 puntos]
Verique 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 gracar 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ácos se pide que deduzcan
el signicado de la solución de D' Alembert.
Caso 1
1 cm
f (x) =
(x/5 cm)2 + 1
g(x) = 0
A modo de ejemplo, un programa Matlab que realiza el primer caso con c=1 es
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
Indique en el informe qué se observa y cómo se interpretan las soluciones de D' Alembert.
3. Experiencia 2 [1 puntos] Si
u(x, t) = f (x − ct) + g(x + ct)
1 cm
f (x) = g(x) =
(x/5 cm)2 + 1
Vista en perspectiva
En el laboratorio disponemos de un sistema que tiene varillas de dos largos diferentes, con
momentos de inercia diferentes I1 e I2 .
A partir de lmaciones de las oscilaciones del sistema completo, se busca que Uds deter-
minen 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 obtenido.
Indique el efecto de tener un mayor o menor I en la propagación de las ondas
5. Conclusiones [1 punto]
B. Materiales
Balanza digital.
C. Experimento
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 en 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 (en pixeles) y su desviación estándar σl (en
pixeles).
3. Calcular las frecuencias de los tres primeros modos normales para el montaje experi-
mental de la experiencia 2.
Para ello se hará caer una carga de garbanzos en forma vertical, se medirá la tensión del cordel que
soporta el plano inclinado y se cuanticará la interacción (fuerza) entre el chorro y la supercie
donde ocurren los impactos. Con esto se espera que:
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;
• Tubo (15 cm);
9
Entendemos por pulso al chorro de extensión longitudinal acotada. La duración de ellos es corta en relación
a las escalas de tiempo características de las cantidades físicas a analizar.
• Un puñado de garbanzos; y
• Programa MATLAB Media.m (disponible en página web del curso)
(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 n 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 jo entre 35
o 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.
o o o
En la Experiencia 2 se medirán las fuerzas medias para inclinaciones de {35 , 45 , 55 }; 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;
8
B∆=
(ρ A⊥ ∆) g .
3
La cantidad entre paréntesis tiene dimensiones de densidad×volumen=masa. En particular, si
∆ fuese la extensión longitudinal (aproximada) del pulso de garbanzos y ρ fuese más o menos
uniforme, entonces (ρ A⊥ ∆) corresponderá a la masa M de todos los garbanzos. Así entonces,
8M g
∆= .
B
La determinación experimental de B permite inferir la extensión aproximada del pulso de gar-
banzos.
7A.8. El informe
Para que el informe se considere completo este debe incluir
1. Tabulación de los datos de los valores medios de la fuerza en función de la altura o ángulo,
según corresponda. Representar estos datos grácamente. Tanto la tabla como el gráco
puede ser anexados en hojas separadas (corcheteadas al informe). Puede usar papel
milimetrado.
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.
obtuvieron la información.
10
Es siempre pertinente citar los autores de trabajos no realizados por uno.
7A.8.1. No olvide:
Rotular debidamente sus grácos;
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 simplicar ideas supongamos que en este intervalo se identican N valores de la fuerza en
instantes uniformemente espaciados (δti = δt), entonces denimos la fuerza media como:
N N N
X X 1 X
Fm δti = Fi δti ⇒ Fm = Fi .
N
i=1 i=1 i=1
El signicado geométrico detrás de esta denició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.
X
δpi = ∆p = pb − pa ,
i
X X
Fi δti = δpi ⇒ Fi δti = δpi = ∆p ,
i i
obtenemos
∆p = Fm ∆t .
P
Aquí hemos identicado ∆t = i δti = (t1 − t0 ) + (t2 − t1 ) + . . . + (tN − tN −1 ) = tb − ta .
J ≡ i Fi δti se le denomina impulso, que en el contexto de la segunda ley
P
A la cantidad 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 n 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 gura 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.
11
c
>> Media('colisionH45b.txt',250,105)
Durac: 0.210 s
Fondo: 2.421 N
Media: 2.355 N
FzaMedia: 0.066 N
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 2
m/s .
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 signica 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 sensor
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 yB 0 ≡ B − mG g . (3)
0
Determinados A y B , obtenga B y
mG con las dos ecuaciones para F− y F+ .
Nota: La regresión lineal se debe hacer con la función polyt de Matlab, explicada en
la guía práctica de la unidad 4C. Si tiene dudas de su utilización escriba help polyt
en la línea de comando de Matlab.
Figura 38: 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 calibre, 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 Po , sin que
estos estén sumergidos bajo el agua. Compruebe que Po = Mo g , donde Mo es la masa
del objeto.
5. Con los cuerpos completamente sumergidos en agua, repita las mediciones del peso
de cada objeto (P1 ).
6. Determine la fuerza de empuje E que experimenta cada objeto a partir de las medidas
de Po ± ∆Po P1 ± ∆P1 . Compare sus mediciones de E con la predicción del Principio
y
de Arquímedes, es decir E = ρf V g .