Péndulo compuesto

El péndulo compuesto es un sólido en rotación, alrededor de un eje fijo perpendicular a la varilla que pasa por O
La ecuación de la dinámica de rotación se escribe
IO·α=-mgxsinθ
Donde x es la distancia entre el centro de masa y el centro de oscilación O.
IO es el momento de inercia del cuerpo respecto del eje de rotación que pasa por O.
Cuando la varilla se separa un ángulo θ de la posición de equilibrio y se suelta, sobre el sólido actúa el momento del peso, que tiene signo contrario al desplazamiento.
Expresamos la ecuación de la dinámica de rotación en forma de ecuación diferencial
Esta no es la ecuación diferencial de un Movimiento Armónico Simple. Si la amplitud es pequeña aproximamos el seno del ángulo al ángulo medido en radianes sinθ≈θ. La ecuación diferencial se escribe entonces
Esta es la ecuación diferencial de un M.A.S. de frecuencia angular ω y periodo P
Por el teorema de Steiner
IO=IC+mx2=mR2+mx2
R se denomina radio de giro, para una varilla R2=l2/12, siendo l la longitud de la varilla. El periodo se escribe
Cuando se representa P en función de x. Aparecen dos curvas simétricas con respecto a la posición de centro de masas. El periodo alcanza un valor infinito para x=0, es decir, cuando coincide el centro de masa con el centro de oscilación O. La curva presenta un mínimo para un cierto valor de x que se puede calcular derivando P respecto de x e igualando a cero.
%péndulo compuesto de l=1 m de longitud fplot(@(x) 2*pi*sqrt((1/12+x.^2)./(9.8*x)),[0.05,0.45]) xlim([0,0.5]) grid on xlabel('x (m)') ylabel('P(s)') title('Periodo de un péndulo compuesto')
Actividades
Se mide el periodo de las oscilaciones para cada una de las posiciones del centro de oscilación.
El péndulo compuesto es una varilla de longitud l=1 m en la que se han hecho agujeros equidistantes 5 cm. El péndulo aparece oscilando en el primer agujero.
Se selecciona el centro de oscilación en el control Posición, se pulsa el botón titulado Nuevo
Momento de inercia y aceleración de la gravedad
Elevando al cuadrado la fórmula del periodo P del péndulo compuesto
Medimos el periodo Pi de péndulo para cada posición xi, completando una tabla con N pares de datos
x (cm) | P (s) |
---|---|
5 | 2.620 |
10 | 1.936 |
15 | 1.668 |
20 | 1.568 |
25 | 1.520 |
30 | 1.512 |
35 | 1.536 |
40 | 1.576 |
45 | 1.600 |
MATLAB nos permite ajustar los datos experimentales del cuadrado del periodo y=P2/(4π2) para cada posición x, a la función no lineal y=a/x+bx mediante la función
x=(5:5:45)/100; y=[2.620 1.936 1.668 1.568 1.520 1.512 1.536 1.576 1.600].^2/(4*pi^2); hold on %representa los datos experimentales plot(x,y,'bo','markersize',6,'markerfacecolor','b') %modelo de función f_ajuste =@(a,x) a(1)./x+a(2)*x; a0=[0.3 4]; %valor inicial de los parámetros af=nlinfit(x,y,f_ajuste,a0) %representa la función x=linspace(0.04,0.45,100); y=f_ajuste(af,x); plot(x,y,'r') title('Péndulo compuesto') xlabel('x (m)') ylabel('P^2/(4·\pi^2)') hold off
El valor de a es el primer elemento del vector
af = 0.0084 0.1020
Los coeficientes a y b
La aceleración de la gravedad vale g=1/b=9.8 m/s2. Conocida la aceleración de la gravedad g, obtenemos el momento de inercia del péndulo respecto de un eje que pasa por el centro de masas, mR2=0.0084·9.8·m=0.0823·m kg·m2
Fuerzas en el centro de oscilación

En el centro de oscilación actúa una fuerza F desconocida que forma un ángulo φ con la dirección radial (la que une el centro de oscilación con el centro de masa)
Ya hemos deducido al principio de esta página, la ecuación de la dinámica de rotación alrededor de un eje perpendicular que pasa por el centro de oscilación
Tomando momentos respecto del centro de masa, la ecuación de la dinámica de rotación se escribe
Ic es el momento de inercia respecto de un eje perpendicular al plano de oscilación del péndulo que pasa por el c.m.. El teorema de Steiner relaciona los dos momentos de inercia.
Aplicamos la ecuación de la dinámica del movimiento circular uniforme, para describir el movimiento en la dirección radial, de forma similar al péndulo simple
Tomando el nivel cero de energía potencial en el centro de oscilación, el principio de conservación de la energía se escribe
Despejamos las componentes, Fcosφ y Fsinφ
En el péndulo simple, solamente tenemos componente en la dirección radial, la tensión de la cuerda
Las componentes rectangulares de la fuerza F en el centro de oscilación son
Horizontal, Fx
Vertical, Fy
La componente Fx es similar al péndulo simple salvo el factor de escala que multiplica al peso mg. La componente Fy es similar al péndulo simple salvo el factor de escala que multiplica al peso y el término adicional Ic/(md2)
Representamos , el término proporcional a Fx entre corchetes, en función de θ para varias amplitudes
hold on for th=[30,60,90,120,150,179]*pi/180 %amplitudes fplot(@(x) -(3*cos(x)-2*cos(th)).*sin(x),[-th,th], 'displayName',num2str(th*180/pi)) end hold off legend('-DynamicLegend','location','best') xlabel('\theta') ylabel('F_x') grid on title('Componente F_x')
Representamos , el término proporcional a Fy sin el término constante, en función de θ para varias amplitudes
hold on for th=[30,60,90,120,150,179]*pi/180 %amplitudes fplot(@(x) (3*cos(x)-2*cos(th)).*cos(x),[-th,th], 'displayName',num2str(th*180/pi)) end hold off legend('-DynamicLegend','location','best') xlabel('\theta') ylabel('F_y') grid on title('Componente F_y')
El lector interesado, puede estudiar en detalle estas funciones, máximos, mínimos y los puntos de corte con el eje horizontal, es decir, los ángulos para los cuales la componente Fx se anula
Referencias
Peter F Hinrichsen. Compound pendulum pivot force. Eur. J. Phys. 42 (2020) 025001