El péndulo electrostático

Condensador cargado aislado

Ente las placas de un condensador plano-paralelo, se coloca un péndulo, atando al extremo de un hilo una pequeña esfera de metal normalmente, una bolita que se ha hecho plegando con los dedos papel de aluminio.

Una de las placas del condensador se conecta a un electroscopio, la otra placa se conecta a tierra.

Ponemos en marcha un generador de Van de Graaff y transportamos carga desde la esfera del generador a la placa del condensador conectada al espectroscopio mediante una bola metálica situada en el extremo de una varilla de plástico, o bien, conectamos con un cable esta placa con la esfera del generador.

En el primer caso, la placa se va cargando, la aguja del electroscopio se va desviando de su posición de equilibrio. La bolita metálica es atraída hacia la placa cargada. En un momento dado entra en contacto con ella y el péndulo empieza a moverse rápidamente entre las placas del condensador, descargando la placa unida al electroscopio.

Como vemos en la figura, la placa positiva del condensador está más cerca de la carga negativa inducida en la bolita, que de la carga positiva inducida en la misma. Luego, la fuerza neta será atractiva.

Cuando el campo es suficientemente intenso, la bolita toca la placa positiva del condensador. La bolita se carga con una carga positiva Δq y placa pierde la misma carga que ha adquirido la bolita.

El electroscopio indicará la disminución de carga de la placa, la aguja indicadora formará ahora un ángulo ligeramente más pequeño con la vertical.

Una vez que la bolita se ha cargado positivamente, experimenta una fuerza en el sentido del campo, que la mueve rápidamente hacia la placa opuesta conectada a tierra.

La bolita choca con dicha placa y pierde entonces su carga. Rebota y se mueve hacia la posición de equilibrio, rebasándola debido a su energía cinética y comienza un nuevo ciclo.

El resultado es que la placa positiva va perdiendo una carga Δq, cada vez que la bolita entra en contacto con dicha placa. El electroscopio indicará la disminución paulatina de carga en la placa.

Cuando la carga existente en la placa está por debajo de cierto valor, el campo eléctrico entre las placas es pequeño y la bolita es atraída por la placa positiva con una fuerza que no es suficiente para que entre en contacto con ella. La bolita oscila cada vez con menos amplitud debido al rozamiento con el aire hasta que se detiene.

Actividades

Disponemos de un generador de Van de Graaff para suministrar la carga que transportamos a una de las placas del condensador.

Vemos como se va cargando la placa y se desvía la aguja indicadora del electroscopio. La bolita es atraída hacia la placa positiva, hasta que entra en contacto con ella. En ese momento, la bolita empieza a oscilar y la placa positiva del condensador empieza a perder poco a poco carga en cada oscilación de la bolita.

Para repetir la experiencia se pulsa el botón titulado Nuevo.

Condensador conectado a una batería

En primer lugar, no tendremos en cuenta el peso de la bolita conductora que se moverá horiontalmente hacia adelante y hacia atrás chocando con las placas.

En una segunda aproximación, tendremos en cuenta el efecto de la componente tangencial del peso y calcularemos la frecuencia final en el estado estacionario.

En este análisis no tendremos en cuenta el efecto del rozamiento del la bolita al moverse en el aire

Un condensador plano paralelo se conecta a una batería V0 que suministra una carga Q a las placas. Un pequeño conductor de masa m se suspende de un hilo muy largo y se coloca entre las placas del condensador de área S separadas una distancia d.

El pequeño conductor describe un arco de circunferencia bajo la acción de dos fuerzas: el peso mg y la fuerza que ejerce el campo eléctrico E sobre la carga q.

Inicialmente, el conductor se pone en contacto con la placa positiva y adquiere una carga q del mismo signo. Se mueve hacia la placa negativa, choca con dicha placa cambiando el módulo y el sentido de su velocidad. Mientras está en contacto con la placa negativa adquiere una carga q del mismo signo que la placa. El conductor inicia el viaje de vuelta hacia la placa positiva, choca con dicha placa adquiriendo una carga q del mismo signo, cambiando el módulo y el sentido de su velocidad. Completándose el primer periodo de la oscilación de la bolita cargada.

Movimiento de la bolita cargada en un campo eléctrico uniforme

La fuerza que ejerce el campo eléctrico E existente entre las placas del condensador es qE=qV0/d y la aceleración a=qE/m

Sea q una fracción r de la carga Q de las placas

a= rQ V 0 md = rC V 0 2 md = ε 0 rS m ( V 0 d ) 2

donde C=ε0S/d es la capacidad del condensador plano-paralelo

Primer periodo

Supongamos que la bolita conductora se suelta cuando está en contacto con la placa positiva, es acelerada por el campo eléctrico E una distancia d hasta que llega a la placa negativa

La velocidad inicial, cuando la bolita se pone en contacto con la placa positiva es u0=0. La velocidad final v1 cuando alcanza la placa opuesta es

v 1 = u 0 +at d= u 0 t+ 1 2 a t 2 } v 1 2 = u 0 2 +2ad

Cuando la bolita choca con la placa cambia el sentido de su velocidad, su módulo y el signo de la carga. Si el coeficiente de restitución es ε<1, la velocidad inmediatamente después del choque será menor ε·v1

La bolita cargada con una carga q negativa inicia el viaje de vuelta impulsada por el campo eléctrico E. La velocidad v'1 que tiene al establecer contacto con la placa positiva es

v 1 '2 = ( ε v 1 ) 2 +2ad v 1 '2 = ε 2 u 0 2 +2ad( 1+ ε 2 )

Cuando la bolita choca con la placa cambia el sentido de su velocidad, su módulo y el signo de la carga. La velocidad inmediatamente después del choque será u1=ε·v'1

u 1 2 = ε 4 u 0 2 +2ad( ε 2 + ε 4 )

con lo que se completa el primer periodo

Segundo periodo

La bolita cargada con una carga q positiva inicia el viaje de ida impulsada por el campo eléctrico E. La velocidad v2 que tiene al establecer contacto con la placa negativa es

v 2 2 = u 1 2 +2ad= ε 4 u 0 2 +2ad( 1+ ε 2 + ε 4 )

La velocidad inmediatamente después del choque es ε·v2. La bolita cargada con una carga q negativa inicia el viaje de vuelta impulsada por el campo eléctrico E. La velocidad v'2 que tiene al establecer contacto con la placa positiva es

v 2 '2 = ε 6 u 0 2 +2ad( 1+ ε 2 + ε 4 + ε 6 )

La velocidad inmediatamente después del choque será u2=ε·v'2

u 2 2 = ε 8 u 0 2 +2ad( ε 2 + ε 4 + ε 6 + ε 8 )

con lo que se completa el segundo periodo

Después de n periodos

u n 2 = ε 4n u 0 2 +2ad( ε 2 + ε 4 +...+ ε 4n )

N=10;
e=0.9;
y=zeros(1,N);
for n=1:N
    y(n)=sum((e.^(2:2:4*n)));
end
plot(0:N,[0,y],'o','markersize',4,'markeredgecolor',
'r','markerfacecolor','r')
grid on
xlabel('n')
ylabel('u')
title('Hacia el estado estacionario')

Estado estacionario

Cuando n tiende a infinito, el primer término tiende a cero ya que ε<1. El término entre paréntesis es la suma de infinitos términos de una progresión geométrica cuyo primer término es ε2 y la razón es ε2

u 2 =2ad ε 2 1 ε 2

Una vez que se ha alcanzado el estado estacionario, la velocidad después de de chocar con la placa es u y la velocidad antes de chocar con la placa opuesta es u+a·Δt, donde Δt es el tiempo que tarda la bolita en viajar de una placa a la otra.

( u +aΔt ) 2 = u 2 +2ad Δt= 2d( 1ε ) a( 1+ε )

Comprobamos que en el estado estacionario, el incremento de la velocidad de la bolita durante su viaje entre las dos placas u+a·Δt se pierde en el choque

u =ε( u +aΔt )

El periodo P y la frecuencia f son, respectivamente

P=2Δt f= 1 2 a( 1+ε ) 2d( 1ε ) =k( V 0 d )

La frecuencia es proporcional al campo E=V0/d existente entre las placas del condensador

Actividades

Se introduce

Se representa la velocidad un en función del tiempo t. Vemos que la velocidad final tiene hacia un valor límite constante u. Esta velocidad es proporcional a la diferencia de potencial V0 entre las placas del condensador


La bolita conductora se mueve bajo la acción de su peso y de la fuerza que produce el campo eléctrico

Las fuerzas sobre la bolita cargada son:

Descomponemos las fuerzas en la dirección tangencial y en la dirección normal

La ecuación del movimiento de la bolita es

ml d 2 θ d t 2 =mgsinθ+qEcosθ

Si la separación 2d entre las placas es pequeña comparada con la longitud l del péndulo, la amplitud de la oscilación del péndulo es pequeña, hacemos las siguientes aproximaciones

sinθθ, cosθ≈1

La ecuación del movimiento se escribe

d 2 (lθ) d t 2 + g l x= qE m d 2 x d t 2 + ω 2 x= qE m

La solución de esta ecuación diferencial es

x=Asin( ωt )+Bcos( ωt )+C dx dt =ω( Acos( ωt )Bsin( ωt ) )

La constante C es la solución particular ω2C=qE/m

Primer semiperiodo

Los coeficientes A y B se determinan a partir de las condiciones iniciales: en el instante t=0, la posición es x=-d,dx/dt=u0

x= u 0 ω sin(ωt)( d+δ )cos(ωt)+δ dx dt = u 0 cos(ωt)+ω( d+δ )sin(ωt)

con δ=qE/mω2

Para calcular el tiempo t1 que tarda la bolita en alcanzar la placa negativa x=d utilizamos la equivalencia, Asin(ωt)+Bcos(ωt)=Msin(ωt+φ)

dδ=Msin( ωt+φ ){ Mcosφ= u 0 /ω Msinφ=(d+δ) t 1 = 1 ω ( arcsin( dδ ( u 0 /ω ) 2 + ( d+δ ) 2 )+arctan( d+δ u 0 /ω ) )

La bolita alcanza la velocidad v1 al llegar a la placa negativa

v 1 = u 0 cos(ω t 1 )+ω( d+δ )sin(ω t 1 )

Después del choque la bolita cambia el módulo, el sentido de su velocidad y el signo de la carga, ε·v1

Segundo semiperiodo

Ponemos el reloj a cero t=0, y calculamos el tiempo t2 que tarda la bolita en moverse desde la placa negativa x=d hasta la placa positiva x=-d

Las condiciones iniciales son ahora: el instante t=0, la posición de la bolita es x=d y su velocidad es -ε·v1. Dado que la carga q ha cambiado de signo, la ecuación del movimiento es

x= ε v 1 ω sin(ωt)+( d+δ )cos(ωt)δ dx dt =ε v 1 cos(ωt)ω( d+δ )sin(ωt)

con δ=|q|E/mω2

Calculamos el tiempo t2 que tarda la bolita en llegar a la posición x=-d

t 2 = 1 ω ( arcsin( dδ ( ε v 1 /ω ) 2 + ( d+δ ) 2 )+arctan( d+δ ε v 1 /ω ) )

Inmediatamente después del choque con la placa positiva, la velcidad u1 de la bolita es

u 1 =ε( ε v 1 cos(ω t 2 )ω( d+δ )sin(ω t 2 ) )= ε 2 v 1 cos(ω t 2 )+εω( d+δ )sin(ω t 2 )

Esta es la velocidad inicial para el siguiente periodo

Representamos la posición x de la bolita en función del tiempo t

d=0.05; %mitad de la separación entre placas
e=0.9; %coeficente de restituación
ddp=250000; %fem de la batería
w=sqrt(9.8); % el péndulo tiene una longitud de  1 m.

delta=(ddp/(4*d))^2/(4*pi*9e9*2700*0.02*w^2);
N=10; %periodos
v0=0;
P=0;
hold on
for i=1:N
    t1=(asin((d-delta)/sqrt((v0/w)^2+(d+delta)^2))+atan((d+delta)/(v0/w)))/w; 
    x1=@(t) (v0/w)*sin(w*t)-(d+delta)*cos(w*t)+delta;
    v1=@(t) v0*cos(w*t)+w*(d+delta)*sin(w*t);
    [tt,y]=fplot(x1,[0,t1]);
    plot(tt+P,y, 'r');
    P=P+t1;

    v0=-e*v1(t1);
    t2=(asin((d-delta)/sqrt((v0/w)^2+(d+delta)^2))+atan((d+delta)/(-v0/w)))/w;
    x2= @(t) (v0/w)*sin(w*t)+(d+delta)*cos(w*t)-delta;
    v2= @(t) v0*cos(w*t)-w*(d+delta)*sin(w*t);
    [tt,y]=fplot(x2,[0,t2]);
    plot(tt+P,y, 'b');
    P=P+t2;
    v0=-e*v2(t2);
end
hold off
grid on
xlabel('t')
ylabel('x')
title('Péndulo elestrostático')

Estado estacionario

En el estado estacionario, la bolita sale de la placa positiva x=-d con velocidad u y llega a la placa negativa x=d con velocidad v, algo mayor, debido al fuerza que ejerce el campo eléctrico. Después de chocar con la placa negativa, pierde la velocidad que había ganado, u= ε·v. Las ecuaciones son

d= u ω sin(ωt)( d+δ )cos(ωt)+δ v = u cos(ωt)+ω( d+δ )sin(ωt) u =ε v

Despejamos el tiempo t que tarda la bolita en desplazarse desde la placa positiva a la negativa, que será medio periodo P

cos(ωt)= δ(1+ε)d(1ε) δ(1+ε)+d(1ε)

Calculamos también la velocidad u

u = εω(d+δ)sin(ωt) 1εcos(ωt)

Representamos la velocidad inicial un de la bolita cuando sale de la placa positiva en función del tiempo. Vemos que la velocidad tiende hacia un límite constante u. La distancia horizontal entre dos puntos rojos consecutivos es el periodo, que toma un valor constante P

d=0.05; %mitad de la separación entre placas
e=0.9; %coeficente de restituación
ddp=250000; %fem de la batería
w=sqrt(9.8); % el péndulo tiene una longitud de  1 m.

delta=(ddp/(2*d))^2/(4*pi*9e9*2700*0.02*w^2);
N=10; %periodos
v0=0;
P=0;
for i=1:N
    t1=(asin((d-delta)/sqrt((v0/w)^2+(d+delta)^2))
+atan((d+delta)/(v0/w)))/w; 
    v1=v0*cos(w*t1)+w*(d+delta)*sin(w*t1);
    P=P+t1;

    v0=-e*v1;
    t2=(asin((d-delta)/sqrt((v0/w)^2+(d+delta)^2))+
atan((d+delta)/(-v0/w)))/w;
    v2=v0*cos(w*t2)-w*(d+delta)*sin(w*t2);
    P=P+t2;
    v0=-e*v2;
    v(i)=v0;
    tt(i)=P;
end
plot([0,tt],[0,v],'o','markersize',4,'markeredgecolor',
'r','markerfacecolor','r')
grid on
xlabel('t')
ylabel('v')
title('Péndulo elestrostático')

P=2*acos((delta*(1+e)-d*(1-e))/(delta*(1+e)+d*(1-e)))/w
u_inf=e*w*(d+delta)*sin(w*P/2)/(1-e*cos(w*P/2))

P =    0.2012
u_inf =    0.9341

Actividades

Se introduce

Se representa la velocidad un en función del tiempo t. Vemos que la velocidad final tiene hacia un valor límite constante u. Esta velocidad es proporcional a la diferencia de potencial V0 entre las placas del condensador


Referencias

Amin Rezaeizadeh, Pooya Mameghani.Thermodynamic model for bouncing charged particles inside a capacitor. Am. J. Phys. 81 (8) August 2013, pp. 632-635

Matej Forjan, Marko Marhl, Vladimir Grubelnik. Mathematical modelling of the electrostatic pendulum in school and undergarduate education. Eur. J. Phys. 35 (2014) 015922

Kazutoshi Asano. On the theory of an electrostatic pendulum oscillator. Am. J. Phys. 43 (5) May 1975, pp. 423-427