LC en paralelo con la resonancia del circuito R

He estado haciendo un experimento que implica medir el voltaje a través del capacitor a diferentes frecuencias de conducción. Considere el siguiente circuito:ingrese la descripción de la imagen aquí

Un bosquejo de la curva experimental. ( R 1 < R 2 < R 3 ) con $(L,C)$ constante:ingrese la descripción de la imagen aquí

Sin embargo, todavía estoy tratando de entender por qué se produce la resonancia en primer lugar. Mi lógica es la siguiente:

V L y V C están en antifase entre sí (el inductor está en fase con el condensador), por lo que el voltaje de la combinación LC y, por lo tanto, el voltaje del circuito es:

V = V L V C
V está en fase con I R . La corriente del circuito es
I = I L C 2 + I R 2
V Z = ( V Z L C ) 2 + ( V R ) 2
. Ya que V = V L C = V L V C , I L C Z L C = I L C X L I L C X C ,
Z L C = X L X C
1 Z = 1 ( X L X C ) 2 + 1 R 2 Z = | X L X R | R R 2 + ( X L X C ) 2
por resonancia, X L = X C Z = 0 . Pero Z = 0 daría lugar a I así que algo no está bien con mi derivación para Z (impedancia del circuito)?

Además, me han dado la amplitud del pico de resonancia:

A = k ( ω 2 ω 0 2 ) 2 + ( γ ω ) 2
donde k es algún factor de normalización, ω = 2 π F y γ es el coeficiente de amortiguamiento. Para A aumentar con R , ¿estoy en lo cierto al suponer que γ 1 R ? ¿Existe una ecuación diferencial para describir el circuito?

EDITAR: La reactancia interna de la fuente de voltaje ~ 600 Ω , ( R 1 , R 2 , R 3 ) = ( 1 Ω , 2 Ω , 3 Ω ) , /$L=1mH, C=100nF$/ Obtuve la siguiente ecuación diferencial:

I ¨ L C ( t ) + 1 L R o tu t pags tu t R R + R o tu t pags tu t I ˙ L C ( t ) + I L C ( t ) L C = R V 0 ω C o s ( ω t ) ( R + R o tu t pags tu t ) L
, donde ω es la frecuencia angular impulsora y V 0 es el voltaje de conducción pico. Esto significa γ = R o tu t pags tu t R L ( R + R o tu t pags tu t ) incluso si se tiene en cuenta la resistencia de salida, pero una R más alta no conduciría a una mayor γ y por lo tanto menor amplitud?

¿Es su fuente de voltaje ideal o tiene una resistencia interna que no se muestra en su diagrama?
Plantea bien tu pregunta por favor. Como lo tiene actualmente, "R" no tendrá ningún efecto y mucho menos R1, R2 y R3 (donde sea que estén).
Acabo de darme cuenta de que la fuente de voltaje (picoscopio) tiene una reactancia interna de 600 Ω ( R 1 , R 2 , R 3 ) = 1 Ω , 2 Ω , 3 Ω ?
Por favor, vea las ediciones, obtuve una ecuación diferencial. pero no tiene sentido.
@floccinaucinihilipilificator ¿Cuáles son sus valores para L y C ?
@Jan L=1mH, C=100nF

Respuestas (2)

La impedancia equivalente que calculó solo es válida si la mide a través R , sin tener en cuenta la fuente, pero su salida es a través C . Para calcular la función de transferencia, primero debe tener en cuenta la impedancia (considerémosla una resistencia) de la fuente de voltaje, que se puede modelar como una resistencia en serie. De lo contrario, no tiene sentido, ya que una fuente de voltaje se considera cero, o una impedancia muy baja, por lo tanto, una resistencia a través de ella que es más grande que su fuente no contribuirá con nada.

Considere el voltaje a través R ser V 1 . Entonces, la impedancia vista desde el lado de la fuente es Z 1 :

(1) Z L C = s L + 1 s C Z 1 = 1 1 Z L C + 1 R (2) = R s 2 + 1 L C s 2 + R L s + 1 L C (3) V o tu t = V 1 Z 1 Z 1 + R I norte (4) H ( s ) = R R I norte + R s 2 + 1 L C s 2 + R I norte R I norte + R R L s + 1 L C

( 4 ) muestra un paso bajo de segundo orden atenuado y amortiguado, que debería tener una respuesta similar a esta (para algunos valores aleatorios de L = 1 mH , C = 1 m F , R = 100 Ω , R I norte = 50 Ω ):

prueba

Pero sus gráficos muestran solo un pico que se asemeja a un paso de banda. Entonces, en el caso de que sus conexiones estén de acuerdo con su esquema, supongo que solo está midiendo un ancho de banda finito alrededor del pico, algo como esto:

rango

Si este es el caso, entonces necesita aumentar el rango; de lo contrario, es posible que tenga algo diferente, y solo se puede determinar mostrando los valores reales de los elementos (incluidos los parásitos, para L, al menos), la impedancia de salida de la fuente , punto de medición, cómo se mide, rango de frecuencias y amplitudes/magnitudes que estás midiendo.


Dada la edición reciente, aún no ha dicho cuál es el rango de frecuencias que está midiendo, pero aquí está la respuesta modificada de su circuito con sus valores:

nuevo

Esto refuerza aún más mi suposición de que está midiendo la salida en un ancho de banda limitado, creando la ilusión de que es un paso de banda, en lugar de un paso bajo muy subamortiguado.

En cuanto a la amortiguación, como muestra su fórmula recién encontrada, es exactamente como se muestra en mi derivación, lo que significa que lo habitual R / L es ahora R mi q / L , donde R mi q = R I norte | | R . En resumen, ambos R y R I norte contribuir. Y si L tenía una resistencia en serie, eso también habría estado allí. Si te lo pone más fácil, aquí está la función de transferencia reescrita:

H ( s ) = R R I norte + R s 2 + 1 L C s 2 + R I norte | | R L s + 1 L C

Su nueva edición pareció dejar más claro cuál es su pregunta, vea si la actualización ayuda.

Primero, presentaré un método que usa Mathematica para resolver este problema. Cuando estaba estudiando estas cosas, usaba el método todo el tiempo (sin usar Mathematica, por supuesto). Además de que la respuesta de @un ciudadano preocupado es excelente.

Bueno, estamos tratando de analizar el siguiente circuito:

esquemático

simular este circuito : esquema creado con CircuitLab

Cuando usamos y aplicamos KCL , podemos escribir el siguiente conjunto de ecuaciones:

(1) I I = I 1 + I 2

Cuando usamos y aplicamos la ley de Ohm , podemos escribir el siguiente conjunto de ecuaciones:

(2) { I I = V I V 1 R 1 I 1 = V 1 R 1 I 2 = V 1 V 2 R 2 I 2 = V 2 R 3

Sustituir ( 2 ) dentro ( 1 ) , para obtener:

(3) { V I V 1 R 1 = V 1 R 1 + V 1 V 2 R 2 V I V 1 R 1 = V 1 R 1 + V 2 R 3

Ahora, podemos resolver para V 2 :

(4) V 2 = R 1 R 3 V I R I ( R 2 + R 3 ) + R 1 ( R I + R 2 + R 3 )

Donde usé el siguiente código de Mathematica para encontrar ( 4 ) :

In[1]:=FullSimplify[
 Solve[{Ii == I1 + I2, Ii == (Vi - V1)/Ri, I1 == V1/R1, 
   I2 == (V1 - V2)/R2, I2 == V2/R3}, {Ii, I1, I2, V1, V2}]]

Out[1]={{Ii -> ((R1 + R2 + R3) Vi)/((R2 + R3) Ri + R1 (R2 + R3 + Ri)), 
  I1 -> ((R2 + R3) Vi)/((R2 + R3) Ri + R1 (R2 + R3 + Ri)), 
  I2 -> (R1 Vi)/((R2 + R3) Ri + R1 (R2 + R3 + Ri)), 
  V1 -> (R1 (R2 + R3) Vi)/((R2 + R3) Ri + R1 (R2 + R3 + Ri)), 
  V2 -> (R1 R3 Vi)/((R2 + R3) Ri + R1 (R2 + R3 + Ri))}}

Entonces, la función de transferencia viene dada por:

(5) H := V 2 V I = R 1 R 3 R I ( R 2 + R 3 ) + R 1 ( R I + R 2 + R 3 )


Ahora, aplicando esto a su circuito, debemos usar (a partir de ahora, uso las letras minúsculas para la función en el dominio s 'complejo' donde usé la transformada de Laplace ):

  • (6) R 2 = SL
  • (7) R 3 = 1 Carolina del Sur

Entonces, podemos reescribir ( 5 ) :

h ( s ) = R 1 1 Carolina del Sur R I ( SL + 1 Carolina del Sur ) + R 1 ( R I + SL + 1 Carolina del Sur ) =
(8) R 1 CL ( R I + R 1 ) s 2 + RC I R 1 s + R I + R 1

Usando el hecho de que para voltajes sinusoidales sabemos que s := j ω (donde j 2 = 1 y ω = 2 π F ), podemos escribir para la amplitud:

| h _ ( j ω ) | = | R 1 CL ( R I + R 1 ) ( j ω ) 2 + RC I R 1 j ω + R I + R 1 | =
| R 1 | | R I + R 1 CL ( R I + R 1 ) ω 2 + RC I R 1 ω j | =
(9) R 1 ( R I + R 1 CL ( R I + R 1 ) ω 2 ) 2 + ( RC I R 1 ω ) 2

Ahora, usando el hecho de que R I = 600   Ω , C = 100 10 9   F , y L = 1 10 3   H . Entonces obtenemos:

(10) | h _ ( j ω ) | = 10 10 R 1 ( R 1 + 600 ) 2 ω 4 + 2 10 10 ( R 1 ( 17 R 1 1200 ) 360000 ) ω 2 + 10 20 ( R 1 + 600 ) 2

Donde usé el siguiente código de Mathematica para encontrar ( 10 ) :

In[2]:=Clear["Global`*"];
x = R1/(R1 + Ri + c R1 Ri s + c L (R1 + Ri) s^2);
s = I*\[Omega];
Ri = 600;
c = 100*10^(-9);
L = 1*10^(-3);
FullSimplify[
 Sqrt[(ComplexExpand[Re[x]])^2 + (ComplexExpand[Im[x]])^2], 
 Assumptions -> R1 > 0 && \[Omega] >= 0]

Out[2]=(10000000000 R1)/Sqrt[
100000000000000000000 (600 + R1)^2 + 
 20000000000 (-360000 + R1 (-1200 + 17 R1)) \[Omega]^2 + (600 + 
    R1)^2 \[Omega]^4]

El máximo se produce cuando (suponiendo que debe ser positivo):

ω 0 = 2 π F 0 = 100000 R 1 ( 1200 17 R 1 ) + 360000 R 1 + 600    
(11) F 0 = 100000 R 1 ( 1200 17 R 1 ) + 360000 2 π ( R 1 + 600 ) = 50000 R 1 ( 1200 17 R 1 ) + 360000 π ( R 1 + 600 )

Donde usé el siguiente código de Mathematica para encontrar ( 11 ) :

In[3]:=FullSimplify[
 Solve[{D[%2, \[Omega]] == 0, R1 > 0 && \[Omega] >= 0}, \[Omega]]]

Out[3]={{\[Omega] -> 
   ConditionalExpression[0, 
    0 < R1 < 600/17 (1 + 3 Sqrt[2]) || 
     R1 > 600/17 (1 + 3 Sqrt[2])]}, {\[Omega] -> 
   ConditionalExpression[(100000 Sqrt[360000 + (1200 - 17 R1) R1])/(
    600 + R1), 0 < R1 < 600/17 (1 + 3 Sqrt[2])]}}