El péndulo giratorio de Pohl
El péndulo de Pohl es un sistema oscilante que consta de un anillo de cobre unido a un muelle helicoidal que puede girar alrededor de un eje horizontal.
El disco se frena mediante las corrientes de Foucault que genera el campo magnético producido por una bobina en el anillo de cobre. Como se estudia en la página titulada “Corrientes de Foucault” el momento de las fuerzas que ejerce el campo magnético sobre las corrientes inducidas es proporcional a la velocidad angular de rotación y de sentido contrario a ésta.
La intensidad del campo magnético es proporcional a la corriente i que pasa por la bobina, la fuerza sobre dichas corrientes es también proporcional al campo magnético. El momento de frenado es proporcional, por tanto, al cuadrado de la intensidad de la corriente que pasa por la bobina.
La fuerza oscilante se proporciona mediante un motor de velocidad variable, que dispone de una rueda impulsora y una excéntrica unida a una biela. La biela se atornilla a una varilla que puede girar alrededor del mismo eje y cuyo extremo está unido al muelle helicoidal. La varilla dispone de una ranura que permite ajustar la amplitud de la oscilación forzada. La varilla impulsora y el disco giran independientemente uno del otro, solamente están conectados por el muelle helicoidal.
Oscilaciones amortiguadas
Se desplaza el disco de la posición de equilibrio y se suelta
La ecuación de la dinámica de rotación del anillo de cobre es
Iα=-kθ- λω
-
α es la aceleración angular del disco cuando el indicador del péndulo se encuentra en la posición angular θ.
-
I es el momento de inercia del anillo respecto del eje de rotación
-
k es la constante del muelle helicoidal
-
λω es el momento de rozamiento proporcional a la velocidad angular ω de rotación, debidos a las corrientes de Foucault producidas en el disco de cobre por el campo magnético de la bobina.
Expresamos la ecuación del movimiento en forma de ecuación diferencial
ω0 es la frecuencia natural o propia del oscilador y γ es la constante de amortiguamiento.
El péndulo de Pohl de la marca Leybold-Heraeus que disponemos en el laboratorio de Física de la Escuela de Ingeniería de Eibar tiene una frecuencia de aproximadamente f0=0.5 Hz, un periodo P0=1/f0=2 s, la frecuencia angular propia es ω0=2πf0=π rad/s en ausencia de rozamiento, es decir, cuando no se conecta la bobina a la fuente de alimentación de corriente continua.
Oscilación amortiguada. γ<ω0
La solución de la ecuación diferencial es
ω es la frecuencia angular de la oscilación amortiguada, que difiere poco de la frecuencia propia ω0, si el amortiguamiento γ es pequeño. No se debe de confundir ω la frecuencia de la oscilación amortiguada, con la velocidad angular de rotación del disco.
La velocidad angular de rotación dθ/dt es
Las condiciones iniciales determinan los valores de las constantes A y φ.
En la experiencia simulada, se desplaza el disco de cobre 150º, θ0=5π/6 de la posición de equilibrio y se suelta, momento en el que empieza a contar el tiempo t=0.
Las sucesivas posiciones de máximo desplazamiento angular, para las cuales la velocidad angular del disco es cero, se denominan posiciones de retorno y se calculan poniendo dθ/dt=0
Medida de la constante de amortiguamiento
En el instante tm=mπ/ω, el máximo desplazamiento del indicador del péndulo es
El cociente
Ejemplo. En el instante t=0, la posición de máximo desplazamiento es la posición inicial θ0=150º, en el instante t9=9.0 el desplazamiento máximo medido es -48.3º.
intensidad=0.3; gamma=1.396*intensidad^2; %amortiguamiento w0=pi; %frecuencia propia w=sqrt(w0^2-gamma^2); %frecuencia de las oscilación amortiguadas phi_0=5*pi/6; %ángulo inicial t=linspace(0,10,200); x=phi_0*exp(-gamma*t).*(gamma*sin(w*t)/w+cos(w*t)); plot(t,x); grid on set(gca,'YTick',-pi:pi/6:pi) set(gca,'YTickLabel',{'-\pi','-5\pi/6','-2\pi/3','-\pi/2','-\pi/3' ,'-\pi/6','0','\pi/6','\pi/3','\pi/2','2\pi/3','5\pi/6','\pi' }) xlabel('t') ylabel('\theta') title('Oscilaciones amortiguadas')
La constante de amortiguamiento γ vale
Comprobación
t9=9·π/ω=9.007 s
Oscilación crítica. γ=ω0
La solución de la ecuación diferencial es
La velocidad angular de rotación dθ/dt es
Las condiciones iniciales determinan los valores de las constantes A y B.
En la experiencia simulada, se desplaza el disco de cobre 150º, θ0=5π/6 de la posición de equilibrio y se suelta, momento en el que empieza a contar el tiempo t=0.
θ0=B
0=A-γB
Despejamos del sistema de dos ecuaciones con dos incógnitas A y B.
El oscilador tiende a la posición de equilibrio θ=0, después de un tiempo t→∞, sin oscilar
Oscilación sobreamortiguada. γ>ω0
La solución de la ecuación diferencial es
La velocidad angular de rotación dθ/dt es
Las condiciones iniciales determinan los valores de las constantes A y B.
En la experiencia simulada, se desplaza el disco de cobre 150º, θ0=5π/6 de la posición de equilibrio y se suelta, momento en el que empieza a contar el tiempo t=0.
θ0=A+B
0=β(B-A)-γ(A+B)
Despejamos del sistema de dos ecuaciones con dos incógnitas A y B.
El oscilador tiende a la posición de equilibrio θ=0, después de un tiempo t→∞, sin oscilar.
intensidad=2; gamma=1.396*intensidad^2; %amortiguamiento w0=pi; %frecuencia propia beta=sqrt(gamma^2-w0^2); phi_0=5*pi/6; %ángulo inicial t=linspace(0,10,200); x=(phi_0/2)*((1-gamma/beta)*exp(-beta*t) +(1+gamma/beta)*exp(beta*t)).*exp(-gamma*t); plot(t,x); grid on set(gca,'YTick',0:pi/6:pi) set(gca,'YTickLabel',{'0','\pi/6','\pi/3','\pi/2','2\pi/3','5\pi/6','\pi' }) xlabel('t') ylabel('\theta') title('Oscilaciones sobreamortiguadas')
Actividades
Se introduce,
-
La intensidad i de la corriente que pasa por la bobina, lo que equivale a una constante de amortiguamiento γ=1.396·i2
-
La frecuencia propia o natural del oscilador se mantiene fija e igual a ω0=π rad/s
Se pulsa el botón titulado Nuevo.
Se observa las oscilaciones amortiguadas del péndulo de Pohl
Oscilaciones forzadas
En un oscilador amortiguado, la oscilación desaparece al cabo de un cierto tiempo, teóricamente infinito. En esta sección se describe el mecanismo que permite mantener la oscilación.
Disponemos de un motor de velocidad angular variable ωf. El eje M del motor está situado a una distancia horizontal c y vertical a del eje O del anillo de cobre donde situamos el origen. Una excéntrica de radio r=|MB| está unida a una biela AB. El extremo A de la biela de longitud b está unida a una varilla excitadora AP que puede girar alrededor del eje O. La distancia R=OA se puede modificar mediante un tornillo que sujeta el extremo A de la biela a una ranura existente en la varilla. Modificando R se modifica la amplitud de la fuerza oscilante.
Como el motor gira con velocidad angular constante, en un determinado instante el ángulo que forma la excéntrica con la horizontal θ=ωf.·t. Dado el ángulo θ, calcularemos el ángulo φ, que forma la varilla con la vertical.
Sea α el ángulo que forma la biela con la horizontal,
R·cosφ-b·sinα=a-r·sinθ
R·sinφ+b·cosα=c+r·cosθ
Eliminamos el ángulo α de este sistema de dos ecuaciones
(Rcosφ-a+rsinθ)2+(c+rcosθ-Rsinφ)2=b2
y teniendo en cuenta que φ es un ángulo pequeño, por lo que sinφ≈φ y cosφ≈0
En la figura, se muestra la gráfica del ángulo de la fuerza φ en función de θ, para los siguientes datos r=1.0, b=22.0, a=10.0, c=22.0,
-
La curva roja se ha trazado para R=7.0
-
La curva azul se ha trazado para R=12.0
r=1; b=22; a=10; c=22; theta=(0:360)*pi/180; hold on; for R=[7,12] phi=(R^2+r^2+a^2+c^2-b^2-2*R*a+2*R*r*sin(theta) -2*a*r*sin(theta)+2*c*r*cos(theta)). /(2*c*R+2*r*R*cos(theta)); plot(theta,phi) end hold off legend('7','12','location','southwest') grid on set(gca,'XTick',0:pi/6:2*pi) set(gca,'XTickLabel',{'0','\pi/6','\pi/3','\pi/2','2\pi/3' ,'5\pi/6','\pi','7\pi/6', '4\pi/3','3\pi/2','5\pi/3','11\pi/6','2\pi'}) xlim([0,2*pi]) xlabel('t') ylabel('\theta') title('Fuerza oscilante')
Ambas curvas difieren de un MAS, aunque pueden hacerse próximas a un MAS con una elección apropiada de los parámetros.
Como vemos en la figura, cambiando R, la distancia entre el eje de rotación del anillo de cobre O y el extremo A de la biela cambia la amplitud de la fuerza oscilante.
En la simulación supondremos que la fuerza oscilante produce un momento respecto del eje de rotación del anillo, descrito por la función armónica
Mf=M0f·cos(ωf·t)
siendo M0f· la amplitud y ωf la frecuencia angular.
La ecuación del movimiento del disco se escribe ahora
Iα=-kθ- λω+Mf
-
α es la aceleración angular del disco cuando el indicador del péndulo se encuentra en la posición angular θ.
-
I es el momento de inercia del anillo respecto del eje de rotación
-
k es la constante del muelle helicoidal
-
λω es el momento de rozamiento proporcional a la velocidad angular ω de rotación, debidos a las corrientes de Foucault producidas en el disco de cobre por el campo magnético de la bobina.
-
Mf es el momento de la fuerza oscilante producida por el motor respecto del eje de rotación O.
Expresamos la ecuación del movimiento en forma de ecuación diferencial
- donde ω0 es la frecuencia natural o propia del oscilador
- ωf es la frecuencia angular del momento producido por la fuerza oscilante de amplitud M0f
- γ es la constante de amortiguamiento,
Como mostramos en la página titulada “El estado transitorio y su evolución hacia el estado estacionario”, la solución de la ecuación diferencial con las condiciones iniciales θ=0, y dθ/dt=0 en el instante t=0, el indicador del anillo parte del origen con velocidad angular inicial nula, es
intensidad=0.3; gamma=1.396*intensidad^2; %constante amortiguamiento w0=pi; %frecuencia propia rFuerza=9; %amplitud frec=0.5; %frecuencia de la fuerza oscilante wf=2*pi*frec; A=((17.0-rFuerza)/10)/((w0^2-wf^2)*(w0^2-wf^2)+4*gamma^2*wf^2); w=sqrt(w0^2-gamma^2); t=linspace(0,20,400); x=A*((w0^2-wf^2)*(cos(wf*t)-exp(-gamma*t).*cos(w*t)) +2*gamma*wf*(sin(wf*t)-((w0^2+wf^2)/(2*w*wf))*exp(-gamma*t).*sin(w*t))); plot(t,x); grid on set(gca,'YTick',-pi:pi/6:pi) set(gca,'YTickLabel',{'-\pi','-5\pi/6','-2\pi/3','-\pi/2' ,'-\pi/3','-\pi/6','0','\pi/6','\pi/3','\pi/2','2\pi/3','5\pi/6','\pi' }) xlabel('t') ylabel('\theta') title('Oscilaciones forzadas')
Actividades
Se introduce,
-
La Intensidad i de la corriente que pasa por la bobina, lo que equivale a una constante de amortiguamiento γ=1.396·i2
-
La frecuencia del momento de la fuerza oscilante f, en el control titulado Frecuencia fuerza, ωf=2πf rad/s.
-
La frecuencia propia o natural del oscilador se mantiene fija e igual a, f0=0.5 Hz, ω0=π rad/s
-
La amplitud del momento de la fuerza oscilante M0f/I, en el control titulado Amplitud.
Se pulsa el botón titulado Nuevo.