Laboratorio Final de Circuitos Eléctricos

Programa de aplicación - Circuitos trifásicos.

Integrantes:
Samuel Martínez
Esteban Gallego
Juan Fernando Bermúdez

Sabado, 19 de Julio del 2025

1. RESUMEN:

El mundo moderno está en constante cambio. Es fundamental ser conscientes de este hecho y adaptarnos a los
cambios, haciendo uso de las herramientas que las nuevas tecnologías brindan. Una de estas herramientas es la
simulación, que permite a los ingenieros tener un acercamiento con que pasaría si se implementa alguna idea o
solución que tienen en la realidad, con la facilidad de solamente dar unos clics y sin tener que arriesgarse a errores
humanos a la hora de realizar cálculos. La idea de este laboratorio es utilizar la herramienta de Matlab con el fin
de simular un escenario que podría darse en la realidad y mediante este informe se busca exponer los resultados
obtenidos y lo que hizo falta para llegar a ellos.

2. OBJETIVOS:

Entender el funcionamiento de Matlab como un programa de simulación útil para modelar la realidad

Simular las conexiones de 4 cargas trifásicas a un sistema trifásico de generación y entender que efecto tendrán
sobre el sistema.

Aprender como se realiza el cálculo de la facturación de energía eléctrica para un hogar y para una industria.

Comprender como se mejora el factor de potencia de una industria con el fin de que no se le cobre más debido
a penalización

3. MARCO TEÓRICO:

La ley de Ohm establece que:

V = IZ (1)

Donde V es la tensión, I es la corriente y Z es la impedancia. Está fórmula resulta ser de mucha utilidad siempre que
se analizan circuitos eléctricos, ya que a partir de ella es posible despejar V, I o Z según el caso lo requiera.

1
2 Universidad Nacional de Colombia - Sede Manizales

Sea Z la impedancia, la definición de impedancia varia dependiendo del elemento que se analice, para los resistores,
la impedancia es sencillamente igual a la resistencia, pero para los inductores la impedancia es:

Z = jwL (2)

Donde j es la unidad imaginaria, L es la inductancia en Henrios y ω es la frecuencia angular la cual se calcula así:

w = 2π f (3)

Siendo f la frecuencia en hertz de la red eléctrica que alimente el circuito, en el caso de Colombia 60Hz
Para un capacitor, la impedancia viene dada por

1
ZC = (4)
jωC

La potencia en DC viene dada por:

P = VI (5)

Por su parte, la potencia en AC viene dada por el triángulo de potencias:

Figura 1: Triángulo de potencias

A partir de este triángulo y la definición clásica de potencia se obtiene que lo que originalmente era potencia se
vuelve la potencia aparente, es decir:

S = VI (6)

En AC, tanto V que es la tensión como I que es la corriente se trabajan en su valor RMS, que es lo que se mide
mediante un multímetro clásico, y es lo que se obtendría al pasar la señal AC a DC. Para este laboratorio se va a
necesitar también la potencia activa la cuál se puede deducir del triangulo de potencias como:

P = V ICos(ϕ) (7)

En esta ecuación ϕ es el ángulo de desfase entre V e I, tomando a V como referencia. Luego la potencia reactiva es:

Q = V IS in(ϕ) (8)
3 Universidad Nacional de Colombia - Sede Manizales

La diferencia entre estas 2 potencias radica en que, mientras P es la potencia útil en el sistema, Q es una potencia
que utilizan los componentes como inductores o capacitores para mantener sus campos magnéticos.
Un circuito trifásico tiene 4 combinaciones posibles: Y-Y, Y-∆, ∆-Y y ∆-∆. Por su parte, la configuración Y-Y,
despreciando sus impedancias de linea y de fuente luce asi:

Figura 2: Circuito trifásico Y-Y

Un método para analizar este circuito es por nodos, si se hacen las ecuaciones del nodo n, tomando el neutro (el
medio entre las 3 fuentes) como referencia, se obtiene que:

Va − Vn Vb − Vn Vc − Vn Vn
+ + = (9)
Za Zb Zc Z0

Si el circuito está equilibrado, las impedancias Za, Zb y Zc son iguales por lo que es posible nombrarlas de igual
manera, con lo que se obtiene que:

Va + Vb + Vc
!
3 1
= Vn + (10)
Z Z Z0

Suponiendo que la fuente trifásica este equilibrada, Va, Vb y Vc van a tener misma magnitud y tendrán desfase entre
ellas de 120 grados, luego:

Va + Vb + Vc = V∠0° + V∠120° + V∠240° = 0 (11)

Por lo tanto:
!
3 1
Vn + =0 (12)
Z Z0

Para que esto se cumpla Vn debe ser igual a 0, y como la diferencia de potencial entre el Vn y la referencia es 0
idealmente, se supone que por Z0 la corriente que pasa es 0. Gracias a este hecho es posible cambiar el valor de
Z0 por cualquier valor sin miedo a afectar los cálculos. Esto es importante, ya que lo siguiente que se hace para
analizar un circuito trifásico equilibrado es pasarlo a su equivalente monofásico, el cual se obtiene de solo analizar
el comportamiento en una fase así:
4 Universidad Nacional de Colombia - Sede Manizales

Figura 3: Equivalente monofásico del circuito trifásico Y-Y

En este caso se cambia la impedancia Z0 por un corto circuito, pero tal y como se vio anteriormente esto no afectara
los cálculos. Gracias a que las cargas y las tensiones son las mismas, este análisis será valido para las 3 fases, y en lo
único que diferirá con hacer el análisis con el diagrama completo es en la corriente que pasa por el neutro.
Como el valor de la impedancia Z0 es irrelevante para el cálculo, es posible poner el neutro para el análisis aún
cuando no halla conexión física entre el neutro de la fuente y el neutro del circuito, ya que en realidad siempre habrá
una impedancia, aunque sea muy alta por medio del ambiente.
La tensión en un circuito trifásico en forma Y-Y tiene 2 maneras de medirse. La llamada tensión de fase que se mide
tomando neutro como referencia y la tensión de linea que se mide tomando otra linea como referencia. Estas se
relacionan matemáticamente así:

VAB = VAN − VBN (13)

Donde VAB es la tensión medida entre la linea A y la linea B, VAN la tensión medida entre la linea A y el neutro
y VBN la tensión medida entre la linea B y el neutro. Las tensiones medidas entre una y otra linea son llamadas
tensiones de linea, mientras que las tensiones medidas entre el neutro y una linea son llamadas tensiones de fase.
Ahora bien, si VAN y VBN pertenecen a la misma fuente trifásica equilibrada, entonces tendrán la misma magnitud y
un desfase de 120 grados entre si de esta manera:

VAN = V∠0° (14)

VBN = V∠240° (15)

Luego la tensión entre las 2 lineas será:

√
VAB = VAN − VBN = V∠0° − V∠240° = 3V∠30° (16)

Una vez establecidas las tensiones de linea y tensiones de fase, es posible hallar la corriente de linea para esta
configuración, la cual viene dada por la fórmula:

VAN
IA = (17)
Za
Donde VAN es la tensión de fase y Za es la impedancia de la rama. Una característica de los circuitos trifásicos en
configuración Y-Y es que la corriente de linea es igual a la de fase.
5 Universidad Nacional de Colombia - Sede Manizales

En caso de ser necesario analizar un circuito trifásico balanceado en delta, es posible pasar este a un circuito trifásico
balanceado en Y mediante la expresión:

Z∆
ZY = (18)
3
Como la potencia que absorbe un capacitor es puramente reactiva, es posible decir que:

S 2 = P2 + Q 2 = Q 2 (19)

Pero S viene dado por:

S =V∗I (20)

Y si se despeja I por ley de Ohm se obtiene que:

S = V 2Z (21)

Es decir que para un capacitor o un inductor siempre se cumple que:

Q = V 2Z (22)

4. PROCEDIMIENTO:

La guía enviada por el monitor brinda los siguientes datos del estado actual de la red:

Tensión linea-linea: 208V

Frecuencia del sistema: 60Hz

Potencia activa total del sistema: 1MW

Carga industrial (Modelada en delta - balanceada):

1. Empresa No. 1:
• Potencia reactiva: 9, 3KV Ar
• Factor de potencia: 0, 85
2. Empresa No. 2:
• Potencia reactiva: 20KV Ar
• Factor de potencia: 0, 87

Carga Residencial (Modelada en estrella - balanceada)

• Número de viviendas: 60
• Factor de potencia: 0,9
• Demanda residencial promedio: 6kW al día

Para poder hacer la simulación es necesario empezar por pasar la tensión que dan a tensión de fase.
6 Universidad Nacional de Colombia - Sede Manizales

V f ase = 208
√ ≈ 120V
3

Luego se calcula la potencia total consumida por las 60 casas:

Ptotal60casas = 6kW ∗ 60 = 360kW

Está es la potencia activa total consumida por las 60 casas. Para poder realizar el análisis, se pasa el circuito a su
equivalente monofásico, para esto se divide la potencia activa total entre 3 para saber la potencia activa por cada fase:

360
3 = 120kW

Con lo que el equivalente monofásico quedaría:

Figura 4: Equivalente monofásico de las 60 casas

Luego es posible calcular la corriente que pasa por la impedancia Za usando la fórmula de potencia activa y
despejando:

I= P
VCos(θ) = 120kW
120(0,9) ≈ 1111A

Está es la magnitud de la corriente, ahora hace falta saber su ángulo de desfase. Sabiendo que el ángulo θ usado en
el calculo del factor de potencia es el desfase entre tensión y corriente, es posible despejar el ángulo asi:

θ = Cos−1 (0, 9) ≈ 25, 84°

Ahora, lo más común en cargas de la vida real es que sean principalmente inductivas, es decir, que la corriente
atrasara a la tensión, y si se toma la tensión como referencia, se obtendrá que los fasores de corriente y tensión para
este circuito son:

V = 120∠0°
I = 1111A∠ − 25, 84°

Con estos datos es posible despejar de la ley de Ohm el valor de la impedancia:

7 Universidad Nacional de Colombia - Sede Manizales

Z= V
I = 120∠0°
1111∠−25,84° ≈ 0, 097 + j0, 047

Ya que se supone una carga mayormente inductiva, de la parte compleja de esta impedancia es posible despejar la
inductancia, así:

j0,047
L= ZL
j2π f = j2π∗60Hz ≈ 0, 00012467H

Por otra parte, para los datos de la empresa 1 se tiene que su potencia reactiva Q es de 9, 3kV Ar, es decir que por
fase tiene una potencia reactiva de:

9,3kV Ar
3 = 3, 1kV Ar

El factor de potencia de esta empresa según la guía es de 0, 85, por lo que el angulo θ es:

θ = Cos−1 (0, 85) ≈ 31, 8°

Es posible, al igual que con la carga de las casas pasar con los datos de tensión de fase antes calculados y con la
potencia reactiva por fase el circuito a su equivalente monofásico.
Luego por la fórmula de potencia reactiva se puede despejar la corriente I, con lo que se obtiene:

Q
I= VS en(θ) = 3,1kV Ar
120V∗S en(31,8°) ≈ 49A

Nuevamente, si se supone carga inductiva la corriente esta en retraso con respecto a la tensión, y tomando a la
tensión como referencia se obtiene que:

V = 120∠0°
I = 49∠ − 31, 8°

Con estos valores se despeja de la ley de Ohm la impedancia y se obtiene que:

Z= 120∠0°
49∠−31,8° ≈ 2, 081 + j1, 29

Pero está impedancia se corresponde con ZY y la guía pide que se modele la empresa 1 como una carga balanceada
en delta, por lo que es necesario convertir esta impedancia a una impedancia Z∆ por medio de la fórmula:

Z∆ = 3ZY = 6, 244 + j3, 871

Nuevamente con esta impedancia se calcula la inductancia:

j3,871
L= ZL
j2π f = j2π∗60Hz ≈ 0, 010268H

Ahora, se realizan los cálculos para la empresa 2. Para esta empresa el factor de potencia es de 0,87 por lo que el
angulo θ es:

θ = Cos−1 (0, 87) ≈ 29, 54°

8 Universidad Nacional de Colombia - Sede Manizales

La potencia reactiva de esta empresa es de 20kV Ar, luego por fase se tendría que:

20kV Ar
3 ≈ 6, 67kV Ar

Luego por la fórmula de potencia reactiva se despeja la corriente y se obtiene que:

Q
I= VS en(θ) = 6,67kV Ar
120V∗S en(29,54°) ≈ 112, 74A

Nuevamente se toma la tensión como referencia, y debido a que es una carga inductiva la corriente atrasa a la tensión
con lo que se obtienen los fasores:

V = 120∠0°
I = 112, 74∠ − 29, 54°

Luego la impedancia por ley de Ohm es:

Z= 120∠0°
112,74∠−29,54° ≈ 0, 926 + j0, 525

Esta es la impedancia ZY , con la cual se despeja la impedancia Z∆ :

Z∆ = 3ZY ≈ 2, 778 + j1, 574

Luego la inductancia es:

j1,574
L= ZL
j2π f = j2π∗60Hz ≈ 0, 004175H

Con estos datos calculados se comienza a montar el circuito en Matlab en la herramienta de Simulink, para elegir
un componente se da doble clic y se busca por el nombre. Para empezar se buscan 3 fuentes AC, y se unen por
el terminal negativo para simular una fuente trifásica balanceada en configuración Y. Luego se busca una carga
trifásica balanceada en Y, se le da doble clic luego de colocarla en el tablero, en la interfaz que despliega se le
colocan conectores expandidos de tres fases, en parametrización se le dice que pida especificar los valores de los
componentes directamente, y en estructura de componente se le dice que es una serie RL. Luego en resistencia y en
inductancia se colocan los valores calculados para la carga de las 60 casas y se cierra esta ventana.
Para la tensión de las fuentes, se calcula la tensión pico, ya que la que se da en la guía es RMS mediante la formula:
√
V = VRMS ∗ 2 ≈ 169, 7V

Este es el valor de tensión que se coloca a las 3 fuentes, luego se desfasan entre si en secuencia positiva.
Nuevamente se da doble clic al tablero y se busca el componente de una carga en delta, a la cual se configura igual
que con la carga en Y, colocando los respectivos valores de resistencia e inductancia calculados para la empresa 1.
Se hace lo mismo para la empresa 2.
Ahora se coloca un sensor de tensión y corriente trifásico mediante el procedimiento de siempre de dar doble clic a
la pantalla y buscar el componente. Este se conecta a la fuente trifásica, es decir, al positivo de las 3 fuentes AC que
se colocaron antes. Luego se coloca un sensor de potencia trifásico para cada una de las cargas, y se conectan las 3
9 Universidad Nacional de Colombia - Sede Manizales

fases de cada carga a su respectivo sensor de potencia. Por último se conectan los sensores de potencia al sensor de
tensión y corriente.
Adicionalmente, para saber la potencia total consumida por la red se conecta un sensor de potencia entre el sensor de
corriente y tensión y la fuente trifásica balanceada.
A cada sensor se le conecta un PS Simulink Converter, y a este se le conecta un display para poder ver los valores
nominales que se necesitan. Entre el Converter y el display de los sensores de tensión y corriente se conecta un RMS
Measurement. Al sensor de corriente y tensión se le conectan un Converter más tanto a corriente como a tensión, al
cual se le pone un Scope osciloscopio. Para que todo funcione se conecta una tierra al neutro de la fuente y a esta
tierra se conecta un solver. A ese neutro se le conecta un sensor de corriente, que está en serie con el neutro de la
conexión en Y de las casas. La simulación queda así:

Figura 5: Simulación del estado actual de la red

Se toman los datos de Fasores de corriente y ondas de corriente de las lineas A, B, C y neutro, estos datos se
consignan en la sección de registro de medidas tomadas en el laboratorio.
Ahora la guía pide corregir el factor de potencia de las empresas en caso de ser necesario. En Colombia se penaliza a
las industrias si tienen un factor de potencia menor a 0, 9, por lo que, con el fin de hacer tan barata como sea posible
esta solución se buscara que ambas empresas tengan este factor de potencia.
En este momento la empresa 1 tiene una potencia reactiva de 9, 3kV Ar, gracias a su factor de potencia de 0, 85, lo
que le da un angulo θ de −31, 8° lo que significa una potencia activa de:

P= 9,3kV Ar
tan(−31,8°) ≈ 14, 999kW

El factor de potencia buscado es de 0, 9 es decir que el ángulo entre la potencia activa y reactiva debe ser de

θ = Cos−1 (0, 9) ≈ 25, 84°

Lo cual, manteniendo la misma potencia activa, significaría una potencia reactiva menor, de tan solo:

Q = P ∗ T an(25, 84°) = 14, 999kW ∗ T an(25, 84°) ≈ 7, 26kV Ar

Es decir, que la potencia reactiva que los capacitores deben cubrir en total es de:
10 Universidad Nacional de Colombia - Sede Manizales

9, 3kV Ar − 7, 26kV Ar = 2, 036kV Ar

Ya que es un sistema trifásico, cada linea o cada fase debe aportar un tercio de esta potencia reactiva:

2,036kV Ar
3 ≈ 0, 6787kV Ar

Pero la potencia reactiva Q de un capacitor viene dada por:

Q = V 2 ZC

Luego, como se quiere que la configuración del banco de capacitores sea en delta al igual que las impedancias, se
realiza el calculo con la tensión de linea, ya que es la que pasara por cada capacitor, así:

Q = Vlinea
2 Z
C

Se despeja la impedancia y se obtiene que:

2
Vlinea (208V)2
ZC = Q = 0,6787kV Ar ≈ 63, 75

Este valor es un número complejo, del cual se puede despejar la capacitancia así:

C= 1
2π f ZC = 1
2π∗60Hz∗63,75 ≈ 4, 16 ∗ 10−5 F

Con esto en mente se implementa el banco de capacitores para la empresa 1 en la simulación, los datos obtenidos
están en la sección de registro de medidas.
Ahora se hacen los cálculos para el banco de capacitores de la empresa 2, la cual tiene una potencia reactiva de
20kV Ar con un factor de potencia de 0, 87, esto da un ángulo θ de 29, 54° lo cual implica una potencia activa de:

P= 20kV Ar
T an(29,54°) ≈ 35, 29kW

Nuevamente se busca un factor de potencia de 0, 9 que implica un ángulo de 25, 84°, esto con la misma potencia
activa quiere decir que la potencia reactiva tendría que ser de tan solo:

Q = P ∗ T an(25, 84°) = 35, 29kW ∗ T an(25, 84°) ≈ 17, 09kV Ar

Es decir, que la potencia reactiva que los capacitores deben cubrir en total es de:

20kV Ar − 17, 09kV Ar ≈ 2, 909kV Ar

Ya que es un sistema trifásico, cada linea o cada fase debe aportar una potencia reactiva de:

2,909kV Ar
3 ≈ 0, 97kV Ar

Por la fórmula de la potencia reactiva Q de un capacitor se despeja la impedancia y se obtiene que:

11 Universidad Nacional de Colombia - Sede Manizales
2
Vlinea (208V)2
ZC = Q = 0,97kV Ar ≈ 44, 602

Luego el valor de la capacitancia debe ser de:

C= 1
2π f ZC = 1
2π∗60Hz∗44,602 ≈ 5, 95 ∗ 10−5 F

Con esto en mente se implementa el banco de capacitores en el simulador, se toman los datos antes y después y se
registran en la siguiente sección.
Ahora la guía pide que se haga facturación de cada vivienda, para este fin se consultan las tarifas actuales de CHEC,
para este mes de Julio de 2025 se publicaron los siguientes valores:

Figura 6: Tarifa CHEC Julio 2025

Para simplificar los cálculos se harán con el precio del kWh sin subsidiar. Así usando el bloque Gain en Simulink
se calcula el precio de la factura para cada vivienda y se toman los resultados. Estos se exponen en el registro de
medidas tomadas en el laboratorio

5. REGISTRO DE MEDIDAS TOMADAS EN EL LABORATORIO:

De la simulación del estado actual de la red se toman la siguiente medida de corriente:

Figura 7: Corriente medida en el estado actual de la red

Las ondas de corriente por su parte dan de la siguiente forma:

12 Universidad Nacional de Colombia - Sede Manizales

Figura 8: Ondas de corriente en el estado actual de la red

Por otra parte la corriente del neutro da:

Figura 9: Corriente del neutro en el estado actual de la red

Debido a que la corriente del neutro es tan pequeña (casi cero) es posible ver que el sistema está balanceado, por
lo tanto las corrientes van a tener todas la misma magnitud, es decir que la medida tomada antes de 1274A es la
magnitud RMS de las 3 corrientes, las cuales están desfasadas entre si 120 grados tal como las fuentes en secuencia
positiva.
El estado actual de la red puede evidenciarse en la potencia que se mide mediante el sensor de potencia:

Figura 10: Estado actual de la red

Aquí se ve que la potencia activa en el sistema actualmente es de 410800W aproximadamente, es decir que la red no
está dando ni la mitad de su capacidad, lo que en primera instancia es una buena señal que podría indicar que quizás
sea plausible cargar los 80 apartamentos a la red.
Antes de implementar el banco de capacitores se toman las siguientes medidas de potencia:
13 Universidad Nacional de Colombia - Sede Manizales

Figura 11: Potencia de la empresa 1 antes del banco de capacitores

Con estos datos el factor de potencia de la empresa 1 es de:

FP = Cos(T an−1 ( 14990

9298
)) ≈ 0, 8498

Luego de implementar el banco de capacitores, se obtienen los siguientes datos:

Figura 12: Potencia de la empresa 1 luego del banco de capacitores

Con lo cual el nuevo factor de potencia de la empresa 1 es de:

FP = Cos(T an−1 ( 14990

7265
)) ≈ 0, 9

Es decir que se cumplió con el objetivo.

Ahora en la empresa 2 se miden los siguientes datos antes de implementar el factor de potencia:

Figura 13: Potencia de la empresa 2 antes del banco de capacitores

Con estos datos se calcula el siguiente factor de potencia:

FP = Cos(T an−1 ( 35300

20010
)) ≈ 0, 87
14 Universidad Nacional de Colombia - Sede Manizales

Se implementa el banco de capacitores y los resultados cambian a ser:

Figura 14: Potencia de la empresa 2 luego del banco de capacitores

Con estos resultados el factor de potencia de la empresa 2 pasa a ser:

FP = Cos(T an−1 ( 17110

35300 )) ≈ 0, 9

Es decir que se cumplió con el objetivo.

Mediante Matlab se hizo también el calculo del precio de la factura para cada vivienda aproximadamente, y se
obtuvo lo siguiente como resultado:

Figura 15: Precio factura de cada residencia

6. CUESTIONARIO:

7. CONCLUSIONES:

Referencias
[1] Dorf, R. C., & Svoboda, J. A. (2010). Circuitos eléctricos (8.ª ed.). México: Pearson Educación.

[2] Nilsson J. W. & Riedel S. A. (2005). Circuitos eléctricos (7.ª ed.). Madrid: Pearson Educación.

[3] Michael Deyvid Estacio Ramirez (2025). Capitulo 7: Programa de aplicación . circuitos trifásicos