Dos condensadores con pérdidas en serie

Un condensador con pérdida de carga es equivalente a un condensador ideal de capacidad C que se descarga a través de una resistencia R. Como ya hemos estudiado en la página anterior, la carga del condensador disminuye exponencialmente con el tiempo 

q=Qexp( t RC )

La resistencia R no suele ser constante sino que depende de la diferencia de potencial V=q/C entre las placas del condensador. Sin embargo, en nuestro estudio supondremos que la resistencia R es constante. Hay condensadores que tiene constantes de tiempo RC del orden de minutos. Sin embargo, hay otros como aceites o plásticos especiales cuyas constantes de tiempo se miden en horas o en días

Un condensador con pérdidas

Consideremos un condensador con pérdidas unido a una batería de potencial V0 en serie con una resistencia 'externa' R0.

  1. La intensidad de la corriente i es la suma de dos contribuciones: la intensidad i1 que pasa por la resietncia R y la variación de carga q del condensador ideal, con el tiempo, dq/dt

  2. i= i 1 + dq dt

  3. Diferencia de potencial entre a y b

  4. V ab = i 1 R= q C

  5. Aplicando la regla de Kirchoff al circuito abcd

  6. V 0 =i R 0 + i 1 R

Combinando las tres ecuaciones obtenemos una ecuación diferencial en i1

i= i 1 +RC d i 1 dt V 0 =( i 1 +RC d i 1 dt ) R 0 + i 1 R V 0 =RC d i 1 dt R 0 + i 1 ( R+ R 0 ) d i 1 dt + R+ R 0 RC R 0 i 1 = V 0 RC R 0

Resolvemos esta ecuación diferencial, con la condición inicial t=0, i1=0

La solución particular es i1=B, introducimos en la ecuación diferencial

B= V 0 R+ R 0

La solución completa de la ecuación diferencial es la suma de la particular más la de la homogénea

i 1 =Aexp( R+ R 0 RC R 0 t )+ V 0 R+ R 0

El coeficiente A se determina a partir de las condiciones iniciales, t=0, i1=0

i 1 = V 0 R+ R 0 ( 1exp( R+ R 0 RC R 0 t ) )

Esta es la corriente que pasa por la resistencia R, la corriente total i es

i= i 1 +RC d i 1 dt i= V 0 R+ R 0 ( 1exp( R+ R 0 RC R 0 t ) )+RC V 0 R+ R 0 R+ R 0 RC R 0 exp( R+ R 0 RC R 0 t ) i= V 0 R+ R 0 ( 1+ R R 0 exp( R+ R 0 RC R 0 t ) )

Después de un tiempo muy grande t→∞ se establece el estado estacionario caracterizado por una corriente

i = V 0 R+ R 0

Cuando se conecta la batería, t=0, la intensidad de la corriente es

i 0 = V 0 R 0

Ejemplo

Una línea a trazos señala la inensidad en el estado estacionario

C=5;
R=1;
R0=10;
V0=100.0; %voltios
hold on
in1=@(t) V0*(1-exp(-(R+R0)*t/(R*C*R0)))/(R+R0);
fplot(in1,[0,30])
in=@(t) V0*(1+R*exp(-(R+R0)*t/(R*C*R0))/R0)/(R+R0);
fplot(in,[0,30])
V_inf=V0/(R0+R);
line([0,30], [V_inf,V_inf],'lineStyle','--')
hold off
xlabel('t')
legend('i_1','i','Location','best')
ylabel('i')
grid on

Agrupación de dos condensadores en serie

  1. Diferencia de potencial entre a y b y entre b y c

  2. V ab = i 1 R 1 = q 1 C 1 V bc = i 2 R 2 = q 2 C 2

  3. La intensidad de la corriente i

  4. i= i 1 + d q 1 dt = i 2 + d q 2 dt i= i 1 + C 1 R 1 d i 1 dt = i 2 + C 2 R 2 d i 2 dt

  5. Aplicando la regla de Kirchoff al circuito abcd e

  6. V 0 =i R 0 + i 1 R 1 + i 2 R 2

Intensidad i1

Obtenemos un sistema de dos ecuaciones diferenciales

{ V 0 =( i 1 + C 1 R 1 d i 1 dt ) R 0 + i 1 R 1 + i 2 R 2 V 0 =( i 2 + C 2 R 2 d i 2 dt ) R 0 + i 1 R 1 + i 2 R 2 { d i 1 dt + R 1 + R 0 C 1 R 1 R 0 i 1 + R 2 C 1 R 1 R 0 i 2 = V 0 C 1 R 1 R 0 d i 2 dt + R 2 + R 0 C 2 R 2 R 0 i 2 + R 1 C 2 R 2 R 0 i 1 = V 0 C 2 R 2 R 0

Para resolverlo, derivamos la primera respecto del tiempo

d 2 i 1 d t 2 + R 1 + R 0 C 1 R 1 R 0 d i 1 dt + R 2 C 1 R 1 R 0 d i 2 dt =0

En la segunda ecuación diferencial, despejamos di2/dt e introducimos en ésta

d 2 i 1 d t 2 + R 1 + R 0 C 1 R 1 R 0 d i 1 dt + R 2 C 1 R 1 R 0 ( V 0 C 2 R 2 R 0 R 2 + R 0 C 2 R 2 R 0 i 2 R 1 C 2 R 2 R 0 i 1 )=0 d 2 i 1 d t 2 + R 1 + R 0 C 1 R 1 R 0 d i 1 dt R 1 C 1 C 2 R 1 R 0 2 i 1 R 2 + R 0 C 1 C 2 R 1 R 0 2 i 2 = V 0 C 1 C 2 R 1 R 0 2

En la primera ecuación diferencial, despejamos i2 e introducimos en ésta

d 2 i 1 d t 2 + R 1 + R 0 C 1 R 1 R 0 d i 1 dt R 1 C 1 C 2 R 1 R 0 2 i 1 R 2 + R 0 C 1 C 2 R 1 R 0 2 ( V 0 R 2 C 1 R 1 R 0 R 2 d i 1 dt R 1 + R 0 R 2 i 1 )= V 0 C 1 C 2 R 1 R 0 2 d 2 i 1 d t 2 + 1 R 0 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 ) d i 1 dt + R 0 + R 1 + R 2 C 1 C 2 R 1 R 2 R 0 i 1 = V 0 C 1 C 2 R 1 R 2 R 0

Resolvemos la ecuación de segundo grado, con las siguientes condiciones iniciales, en el instante t=0, i1=0, i2=0, y de la primera ecuación difeencial,

d i 1 dt | t=0 = V 0 C 1 R 1 R 0

La solución particular es la constante i1=D, introducimos en la ecuación difeencial

R 1 + R 2 + R 0 C 1 C 2 R 1 R 2 R 0 D= V 0 C 1 C 2 R 1 R 2 R 0 D= V 0 R 0 + R 1 + R 2

Solución de la homogénea. La ecuación característica es

s 2 + 1 R 0 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 )s+ R 0 + R 1 + R 2 C 1 C 2 R 1 R 2 R 0 =0 s 1 = 1 R 0 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 ) 1 R 0 2 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 ) 2 4 R 0 + R 1 + R 2 C 1 C 2 R 1 R 2 R 0 2 s 2 = 1 R 0 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 )+ 1 R 0 2 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 ) 2 4 R 0 + R 1 + R 2 C 1 C 2 R 1 R 2 R 0 2

Dos raíces reales negativas

La solución completa, es la suma de la particular más la homogénea

i 1 =A e | s 1 |t +B e | s 2 |t + V 0 R 0 + R 1 + R 2

Los coeficientes A y B se determinan a partir de las condiciones iniciales

{ 0=A+B+ V 0 R 0 + R 1 + R 2 V 0 C 1 R 1 R 0 =| s 1 |A| s 2 |B B= V 0 | s 2 || s 1 | ( 1 C 1 R 1 R 0 | s 1 | R 0 + R 1 + R 2 ) A= V 0 | s 2 || s 1 | ( 1 C 1 R 1 R 0 | s 2 | R 0 + R 1 + R 2 )

Intensidad i2

Calculamos la corriente i2 conocido i1. Partimos de la primera ecuación diferencial

d i 1 dt + R 1 + R 0 C 1 R 1 R 0 i 1 + R 2 C 1 R 1 R 0 i 2 = V 0 C 1 R 1 R 0 i 2 = V 0 R 2 R 1 + R 0 R 2 i 1 C 1 R 1 R 0 R 2 d i 1 dt i 2 = V 0 R 2 R 1 + R 0 R 2 ( A e | s 1 |t +B e | s 2 |t + V 0 R 0 + R 1 + R 2 ) C 1 R 1 R 0 R 2 ( A| s 1 | e | s 1 |t B| s 2 | e | s 2 |t ) i 2 = A e | s 1 |t R 2 ( C 1 R 1 R 0 | s 1 | R 1 R 0 )+ B e | s 2 |t R 2 ( C 1 R 1 R 0 | s 2 | R 1 R 0 )+ V 0 R 0 + R 1 + R 2

Comprobamos que para t=0, i2=0. Como puede comprobarse en la segunda ecuación diferencial

d i 2 dt | t=0 = V 0 C 2 R 2 R 0

Para ello, utilizamos las relaciones entre las dos raíces

| s 1 |+| s 2 |= 1 R 0 ( R 1 + R 0 C 1 R 1 + R 2 + R 0 C 2 R 2 ) | s 1 |·| s 2 |= R 0 + R 1 + R 2 C 1 C 2 R 1 R 2 R 0

Después de muchas operaciones, obtenemos los resultados esperados

t=0{ i 2 = V 0 R 0 + R 1 + R 2 + A R 2 ( C 1 R 1 R 0 | s 1 | R 1 R 0 )+ B R 2 ( C 1 R 1 R 0 | s 2 | R 1 R 0 )=0 d i 2 dt =| s 1 | A R 2 ( C 1 R 1 R 0 | s 1 | R 1 R 0 )| s 2 | B R 2 ( C 1 R 1 R 0 | s 2 | R 1 R 0 )= V 0 C 2 R 2 R 0

Intensidad i

Conocido i1, partimos de la ecuación

i= i 1 + C 1 R 1 d i 1 dt i=A e | s 1 |t +B e | s 2 |t + V 0 R 0 + R 1 + R 2 C 1 R 1 ( A| s 1 | e | s 1 |t +B| s 2 | e | s 2 |t ) i=A e | s 1 |t ( 1 C 1 R 1 | s 1 | )+B e | s 2 |t ( 1 C 1 R 1 | s 2 | )+ V 0 R 0 + R 1 + R 2

En el estado estacionario

i = V 0 R 0 + R 1 + R 2

En el instante t=0

i(0)=A( 1 C 1 R 1 | s 1 | )+B( 1 C 1 R 1 | s 2 | )+ V 0 R 0 + R 1 + R 2 = V 0 R 0

Ejemplo

Una línea a trazos señala la inensidad en el estado estacionario

C1=1; %Faradios
R1=5; %Ohm
C2=5;
R2=1;
R0=10;
V0=100.0; %voltios
a=((R1+R0)/(C1*R1)+(R2+R0)/(C2*R2))/R0;
b=sqrt(((R1+R0)/(C1*R1)+(R2+R0)/(C2*R2))^2/R0^2-4*(R0*R1+R2)/(C1*C2*R1*R2*R0));
s1=(a-b)/2;
s2=(a+b)/2;
A=V0*(1/(C1*R1*R0)-s2/(R0+R1+R2))/(s2-s1);
B=-V0*(1/(C1*R1*R0)-s1/(R0+R1+R2))/(s2-s1);

hold on
in1=@(t) A*exp(-s1*t)+B*exp(-s2*t)+V0/(R0+R1+R2);
fplot(in1,[0,30])
in2=@(t) A*exp(-s1*t)*(C1*R1*R0*s1-R1-R0)/R2+B*exp(-s2*t)*
(C1*R1*R0*s2-R1-R0)/R2+V0/(R0+R1+R2);
fplot(in2,[0,30])
in=@(t) A*exp(-s1*t)*(1-C1*R1*s1)+B*exp(-s2*t)*(1-C1*R1*s2)+V0/(R0+R1+R2);
fplot(in,[0,30])
V_inf=V0/(R0+R1+R2);
line([0,30], [V_inf,V_inf],'lineStyle','--')
hold off
xlabel('t')
legend('i_1','i_2','i','Location','best')
ylabel('i')
grid on

Una situación con solución analítica sencilla

Consideremos ahora la situación que se muestra en la figura. Se establece una diferencia de potencial V entre a y c, se ha suprimido la resistencia R0

 

Estado inicial

Inicialmente las cargas de los condensadores son iguales y las diferencias de potencial entre sus placas son respectivamente V1=q/C1 y V2=q/C2. Por lo que cumple que

V 1 V 2 = C 2 C 1

Como V1+V2=V tendremos que

V 2 = V C 1 C 2 + C 1

La corriente comienza a fluir a través de las resistencias R1 y R2 y en general, las cargas en los condensadores serán distintas.

Estado final

Si la diferencia de potencial entre los extremos a y c se mantiene constante, se alcanza un estado estacionario en el que la misma corriente i pasa por las resistencias R1 y R2. Se cumplirá entonces que V1=iR1 y V2=iR2. Por lo que tendremos la relación

V 1 V 2 = R 1 R 2

Como V1+V2=V tendremos que

V 2 = V R 2 R 1 + R 2

La carga final de cada condensador será q1=V1·C1 y q2=C2·V2, de modo que se cumple la relación

q 1 q 2 = C 1 R 1 C 2 R 2

Evolución del estado inicial al final

Supongamos que a los extremos a y c se aplica una diferencia de potencial constante V.

La corriente total i que pasa por a o que sale en b en el instante t, es la suma de dos términos:

tal como hemos visto al estudiar el circuito formado por dos condensadores y una resistencia

Del mismo modo, la corriente i que entra en b o que sale en c  en el instante t, es la suma de dos términos:

d q 1 dt + i 1 = d q 2 dt + i 2 , C 1 d V 1 dt + V 1 R 1 = C 2 d V 2 dt + V 2 R 2

Sustituyendo V2=V-V1 tenemos una ecuación en V2.

( C 1 + C 2 ) d V 2 dt +( 1 R 1 + 1 R 2 ) V 2 = V R 1

o bien

d V 2 dt + 1 T V 2 = V R 1 ( C 1 + C 2 ) , T= R 2 R 1 ( C 1 + C 2 ) R 1 + R 2

Integramos esta ecuación diferencial con las siguientes condiciones iniciales: en el instante t=0, la diferencia de potencial entre las placas del condensador C2 es V02 tal como vimos al principio de este apartado.

V 02 = V C 1 C 2 + C 1

Integrando

V 02 V 2 d V 2 V R 1 ( C 1 + C 2 ) 1 T V 2 = 0 t dt

El resultado final es

V 2 (t)= V R 1 + R 2 ( R 2 + C 1 R 1 C 2 R 2 C 1 + C 2 exp( t T ) )

Comprobación

En el estado inicial t=0  tenemos que

V 2 = V C 1 C 1 + C 2

En el estado final t →∞ tenemos

V 2 = V R 2 R 1 + R 2

Intensidad de la corriente

La intensidad de la corriente i es igual a

i= C 2 d V 2 dt + V 2 R 2 i= V 0 R 1 + R 2 ( 1+ ( C 1 R 1 C 2 R 2 ) 2 R 1 R 2 ( C 1 + C 2 ) 2 exp( t T ) )

La intensidad es máxima en el instante t=0, y tiende hacia un valor constante V 0 R 1 + R 2 para t→∞

Caso particular

En el caso particular de que C1R1=C2R2 las diferencias de potencial V1 y V2 son independientes del tiempo, aunque la corriente sigue fluyendo a través de cada uno de los condensadores

V 2 = V R 2 R 1 + R 2 , V 1 = V R 1 R 1 + R 2 q 2 = C 2 R 2 R 1 + R 2 V, q 1 = C 1 R 1 R 1 + R 2 V

Como C1R1=C2R2 las cargas en los condensadores q1 y q2 son iguales.

Esta condición se puede cumplir eligiendo adecuadamente la constante dieléctrica k y la resistividad ρ del dieléctrico que se coloca entre las placas de un condensador plano paralelo, cuyas placas tiene un área A y están separadas una distancia d.

El producto de ambas magnitudes solamente depende de la constante dieléctrica k y de la resistividad ρ del dieléctrico, RC=ρkε0.

Siempre que de los dieléctricos que separan las placas de los dos condensadores sean tales que cumplan la relación ρ1k12k2 se cumplirá también que C1R1=C2R2

Ejemplo 1:

Representamos: el cociente V1(t)/V2(t) en función del tiempo t, la intensidad i en función del tiempo t

C1=0.2; %micro Faradios
R1=5000; %mega Ohm
C2=0.5;
R2=1000;
V0=100.0; %voltios
T=R2*R1*(C1+C2)/(R1+R2);
V2=@(t) (V0/(R1+R2))*(R2+(C1*R1-C2*R2)*exp(-t/T)/(C1+C2));
subplot(2,1,1)
fplot(@(t) (V0-V2(t))./V2(t), [0,3500])
xlabel('t')
ylabel('V_1/V_2')
grid on
subplot(2,1,2)
inten=@(t) (V0/(R1+R2))*(1+(C1*R1-C2*R2)*(C1*R1-C2*R2)*exp(-t/T)
/((C1+C2)*(C1+C2)*R1*R2));
fplot(inten,[0,3500])
xlabel('t')
ylabel('i')
grid on

En la parte superior, se representa, V1/V2 en función de t. En el instante inicial t=0 la relación V1/V2=2.50. Después de un tiempo suficientemente grande t →∞, V1/V2=5.0

En la gráfica inferior, se representa la intensidad i, que tiende hacia un valor constante V R 1 + R 2 e igual a 0.0167 cuando t →∞

Ejemplo 2:

Se cumple el caso especial de que C1R1=C2R2. La relación V1/V2=5 en la situación inicial (t=0) y en la final (t →∞). Aunque la carga de cada condensador no cambia, la intensidad es 0.0167 no nula.

Referencias

Kirk T. McDonald. Leaky Capacitors (2001)

French A, P. Are the Textbook Writers Wrong about Capacitors?. The Physics Teacher, Vol 31, March 1993, pp156-159