Departamento de Energía y Automatización
TITULO
DESCRIPCIÓN
FECHA
i
Dedicatoria
ii
Agradecimiento
iii
Resumen
Palabras Clave:
iv
Índice General
v
Índice de Tablas y Figuras
Contenido
1. Introducción........................................................................................1
1.1. Presentación...................................................................................................1
2. Principios Básicos...............................................................................1
2.1. Fenómenos.....................................................................................................1
3. Montaje experimental........................................................................3
3.1. Presentación...................................................................................................3
3.1.1 Análisis del circuito RLC bajo régimen de resonancia permanente..........3
4. Resultados Experimentales y su Análisis.......................................10
4.1. Presentación.................................................................................................10
5. Conclusiones......................................................................................11
6. Bibliografía........................................................................................12
7. Apéndices..........................................................................................13
vi
Lista de Símbolos y Abreviaturas
ε: Fem o fuerza electromotriz.
: Constante dieléctrica del medio.
Φ: Flujo del campo magnético.
1
1. Introducción
1.1. Presentación
2. Principios Básicos
2.1. Fenómenos
Un sistema de calentamiento inductivo consta de tres partes principales: un generador de
corriente alterna, una bobina de calentamiento y la pieza a calentar, que se sitúa en el interior
de la bobina tal y como se puede ver en la figura 2.1.
Fig. 2.1: Diagrama básico de un sistema de calentamiento inductivo
2
12
10
Profundidad de Penetración (mm) 8
60
500
0
3000
10000
280 200 )
160 (Hz
30000
120 80 40 ci a
100000
In tensidad de 10 en
Ca mpo (A/m r e cu
) F
Fig. 2.2: Variación de la profundidad de penetración () en función de la intensidad de campo
magnético (H) y de la frecuencia a temperatura de 21 C.
La tabla 2.2
21 Frecuencia (Hz)
C 60 500 3000 10000 30000 100000
10 2,5 0,88 0,36 0,2 0,11 0,06
40 4,7 1,63 0,67 0,36 0,21 0,12
80 6,3 2,2 0,9 0,49 0,28 0,16
H (A/m)
d (mm)
120 7,76 2,69 1,1 0,6 0,35 0,19
160 8,76 3,03 1,24 0,68 0,39 0,21
200 9,63 3,33 1,36 0,75 0,43 0,24
280 11,2 3,89 1,59 0,87 0,5 0,27
Tabla 2.2. Variación de la profundidad de penetración (mm) para acero al carbono 1045 a
21ºC en función de la intensidad de campo magnético y de la frecuencia.
3
A B C
Fig. 2.3: Distribución de corrientes por efecto proximidad en conductores cilíndricos: A) un
solo cable, B) dos cables con corrientes opuestas, C) dos cables con corrientes en el mismo
sentido
3. Montaje experimental
3.1. Presentación
V 1 = jωL p I 1 − jω MI 2 (4.1)
−V 2 =− jω MI 1 + jωLs I 2 (4.2)
donde j es j2 f, siendo f la frecuencia de trabajo y j=√−1 .
Tomando en consideración que I1 señala hacia el extremo con punto e I2 hacia el extremo sin
punto, la convención de punto dicta los signos menos para los términos de inductancia mutua,
por lo que I2 puede colocarse en términos de I1:
jωM
I 2= I
R s + jωL s 1 (4.3)
V 2
R s=
Donde I2. Utilizando 4.3, para eliminar I2 de 4.1, se tiene
[
V 1 = jωL p +
ω2 M 2
]I
R s + jωL s 1
(4.4)
3.1.1 Análisis del circuito RLC bajo régimen de resonancia permanente
Una de las características mas relevantes del presente trabajo es que el sistema debe funcionar
en régimen de resonancia permanente, por lo que se hace necesario calcular que sucede con un
4
sistema serie al cual se le excita con una señal V de frecuencia igual a la de resonancia. Para
efectos del cálculo, se tomarán los valores L, C y R totales definidos en la figura 4.5, dado que
todos esos parámetros influyen en la frecuencia de resonancia del sistema. Usando la Ley de
Kirchoff de tensiones en el dominio de Laplace para el circuito mostrado en la figura 4.5, se
tiene
V−
I (s )
sC [
−s⋅LI ( s )−R⋅I ( s )=0 ⇒V =I ( s )
1
sC
+ sL+ R
] (4.12)
Por otro lado, el valor que se quiere deducir finalmente es la potencia real sobre el sistema y
ésto implica obtener el voltaje (o la corriente) sobre la resistencia óhmica del sistema, por lo
tanto se puede decir
V R=I (s ) R (4.13)
luego, la función de transferencia del sistema, referenciado a la resistencia es
V R( s) R
=
V (s) 1 + sL+ R
sC (4.14)
Arreglando los términos y normalizando la ecuación, se tiene
R
V R( s) s
L
=
V (s) 2 R 1
s +s +
L LC (4.15)
Al buscar las raíces de la ecuación de segundo orden del denominador de 4.15, se puede
factorizar la misma y por expansión en fracciones parciales obtener la solución temporal de la
ecuación, luego
R
s=− ±
2L √( 2L )
R 2 1
−
LC (4.16)
Ahora bien, es necesario que el sistema tenga una respuesta compleja, ya que se quiere que sea
oscilatoria, por lo tanto se debe cumplir que
5
( )
R 2 1
2L
<
LC (4.17)
luego, 4.15 puede escribirse como
R R
V R( s) s s
L L
= =
V ( s ) ( s +( a+ jω ) )( s +( a− jω )) (( s+ a)2 +ω2 ) (4.18)
donde
a=
R
2L
y jω=
√( 2L )
R 2 1
−
LC . Hay que hacer notar que la definición de frecuencia
1
√
de resonancia difiere ligeramente de la definición clásica LC , debido a que el sistema
trabaja precisamente en el punto de resonancia y no es factible realizar aproximaciones.
Ahora bien, el sistema debe alimentarse con una señal de naturaleza tal que alguna de sus
armónicas coincida con la frecuencia de resonancia del sistema, esto es, V =V i sen(ωt ) por lo
tanto
R R
s s
L L ω
V R ( s )= 2 2
V ( s )= 2 2
⋅V i 2 =
(( s+ a) + ω ) ( ( s+ a ) +ω ) ( s + ω2 )
R sω
Vi
L ( ( s +a )2 +ω 2 )( s 2 +ω 2 ) (4.19)
Realizando la expansión en fracciones parciales y antitransformando 4.19, se obtiene,
finalmente
V R (t )=V i
R 1
L (a2 +4 ω2 ) [
(1−ωe−at
)cos(ωt )+
2 ω 2−(a2 +2 ω2 )e−at
a
sen (ωt ) ] (4.20)
R
2 a=
Dado que L , 4.20 puede escribirse como
2
V R (t )=V i 2 2
[ a (1−ωe−at )cos(ωt )+sen (ωt )(2 ω2−(a2+2ω 2 )e−at )]
(a +4 ω ) (4.21)
6
Figura 4.4. Evaluación, mediante simulación AC, del sistema de salida resonante: (I Rp)
freuencia 36,224KHz, Corriente 22,072A, (IRT) frecuencia 36,224KHz, Corriente 22,150A.
0.8
0.6
0.4
Amplitud (V)
0.2
-0.2
-0.4
-0.6
-0.8
-1
0 0.5 1 1.5
t (s) -3
x 10
Figura 4.5. Respuesta analítica del sistema a la excitación sen(ωt )u(t ) evaluada mediante
7
MATLAB.
R2 R5
63.4k 54.9k
C1 C3
150p 68p
R1 R3 R6
4
63.4k 110k 154k
Va 2
Vee
R4 6
Vee
V-
V-
- -
Vi 1
Vo 54.9k
7
1Vac C2 O U T U 1A C4 O U T U 1B
0Vdc 470p 3 680p 5
V+
V+
+ +
Vcc Vcc
0 0 0 0 0
8
Figura 5.6. Circuito del filtro pasa-bajo propuesto
Figura 5.7. Diagrama de Bode del filtro pasa-bajo Simulado Mediante SPICE.
8
Dado que el filtro está compuesto de dos etapas idénticas, salvo los valores de sus
componentes, se procederá a obtener la función de transferencia de la primera etapa y
posteriormente se multiplicará por otra igual con los valores ajustados.
Nuevamente, haciendo uso de la teoría de amplificadores operacionales y el equivalente de los
componentes reactivos en el dominio de Laplace y usando las tensiones de nodo Vi, Va y Vo
mostradas en la figura 5.7, se pueden plantear las ecuaciones de nodo del filtro
V i−V a V a V a V a −V O
− − − =0
R1 1 R3 R2
sC 2 (5.20)
V a V a −V O
−
R3 1
sC 1 (5.21)
De 5.21 se puede despejar Va
V O sR 3 C 1
V a =−
1−sR 3 C 1 (5.22)
y de 5.20 se puede despejar Vi
( )
R 1 R2 + R1 R 3 + R2 R 3 + sR1 R 2 R 3 C 2 R1 V O
V i =V a R1 −
R 1 R2 R3 R2
(5.23)
Sustituyendo 5.22 en 5.23
( )
V O sR 3 C 1 R 1 R 2 +R1 R3 +R2 R3 +sR 1 R2 R 3 C2 R1 V O
V i =− −
1−sR 3 C1 R 2 R3 R2
(5.24)
Simplificando y sacando factor común Vo
Vi
VO
=−
[
( R1 R2 + R1 R3 + R 2 R3 + sR 1 R2 R3 C 2 ) sC 1 R1
( 1−sR 3 C1 ) R2
+
R2 ] (5.25)
Sumando las fracciones e invirtiendo
9
( )
VO
Vi 1
=−
[ ( 1−sR 3 C 1 ) R 2
( R1 R 2 + R1 R3 + R2 R3 + sR 1 R2 R3 C 2 ) sC 1 + ( 1−sR 3 C1 ) R 1 ] (5.26)
dado que la ganancia de la segunda etapa es funcionalmente idéntica a la primera, se puede
inferir directamente que
( )
VO
Vi 2
=−
[ ( 1−sR6 C3 ) R5
( R 4 R5 + R4 R 6 + R 5 R6 + sR 4 R 5 R6 C 4 ) sC 3 + ( 1−sR6 C3 ) R4 ] (5.27)
Los subíndices de 5.26 y 5.27 se refieren a la etapa particular del filtro.
Sustituyendo los valores dados para el filtro y multiplicando 5.26 y 5.27, se obtiene
( ) [ ]
VO 0 ,6014 s 2−9 ,388⋅104 s +3 , 481⋅10 9
=
V i T 1 , 020⋅10−9 s 4 +1 , 292⋅10−4 s3 + 7 ,766 s 2 +2 , 436⋅10 5 s +3 , 481⋅10 9
(5.28)
El diagrama de Bode de esta función de transferencia, evaluado analíticamente usando Matlab,
puede verse en la figura 5.9, la cual se compara favorablemente con la figura 5.8 obtenida
mediante la simulación SPICE.
-20
Magnitud (dB)
-40
-60
-80
-100
360
270
180
Fase (deg)
90
-90
-180
2 3 4 5 6
10 10 10 10 10
Frecuencia (Hz)
10
Figura 5.8. Diagrama de Bode Analítico del filtro pasa-bajo.
11
4. Resultados Experimentales y su Análisis
4.1. Presentación
12
5. Conclusiones
1. El sistema
13
6. Bibliografía
[1]. Bell, Eric; “A valve-maintained high-frequency induction furnace and some notes on the
performance of induction furnaces”, Proc. Phys. Soc. 40, pp.193-205, 1927.
[2]. Clark D.; Varney, W.; “Physical Metallurgy for Engineers, Second Edition”, pp 132 ff,
Van Nostrand Company, Princeton, New Jersey, 1965.
[3]. Espí, J.M.; “Inversor Resonante de Tres Elementos L-LC con Características
Cortocircuitables para Aplicaciones de Calentamiento por Inducción”, Tesis Doctoral,
Universitat deValència, Valencia, 1999, España
14
7. Apéndices
7.1.
