Universidad Técnica de Cotopaxi

Nombre: Alex Chicaiza

Ciclo: 8vo “A”
Fecha: 26/12/2022
Materia: Planificación y Confiabilidad de SEP
MODELO DC, TRANSPORTE, HIBRIDO Y DISYUNTIVO
En el presento documento se realizará una comparación de los modelos DC, transporte,
hibrido y disyuntivo, tanto en sus modelos matemáticos como en sus resultados, esto con
el fin de ver si los modelos son capaces de obtener un resultado aproximadamente igual
que el modelo DC.
Para esto se utilizarán los siguientes datos de un sistema de 6:
Tabla 1. Datos de potencia activa de los generadores y las cargas.
Barra Tipo Demanda 𝑷𝒅 Generación Generación
de Máxima Mínima
Barra 𝑴𝑾 𝑷𝒎𝒂𝒙 𝑷𝒎𝒊𝒏
MW
MW
1 3 80 150.00 0.00
2 0 240 000.0 0.00
3 2 40 360.00 0.00
4 0 160 000.0 0.00
5 0 240 000.0 0.00
6 2 0 600.00 0.00

Tabla 2. Datos de las líneas de transmisión.

𝑺𝒎𝒂𝒙𝒊𝒋 𝒓𝒍 𝒙𝒍 𝒃𝒔𝒉𝒍 𝑪𝒍
# línea i j
MVA p.u. p.u. p.u. MUSD
𝒏𝟎
𝒏𝒎𝒂𝒙
1 1 2 100 0.1 0.40 0.018 40 1 2
2 1 3 100 0.09 0.38 0.018 38 0 3
3 1 4 80 0.15 0.60 0.062 60 1 2
4 1 5 100 0.05 0.20 0.05 20 1 2
5 1 6 70 0.17 0.68 0.021 68 0 3
6 2 3 100 0.05 0.20 0.044 20 1 2
7 2 4 100 0.10 0.40 0.046 40 1 2
8 2 5 100 0.08 0.31 0.012 31 0 3
9 2 6 100 0.08 0.30 0.011 30 0 3
10 3 4 82 0.15 0.59 0.029 59 0 3
11 3 5 100 0.05 0.20 0.01 20 1 2
12 3 6 100 0.12 0.48 0.021 48 0 3
13 4 5 75 0.16 0.63 0.024 63 0 3
14 4 6 100 0.08 0.30 0.09 30 0 3
 15 5 6 78 0.15 0.61 0.049 61 0 3

Modelo Matemático
➢ Modelo DC
- Función Objetivo:

𝑀𝑖𝑛 {𝐶𝑇 = ∑ (𝐶𝑙𝑖𝑗 ∗ 𝑛𝑖𝑗 )} (1)

𝑖𝑗∈Ω𝐿

- Restricciones:
Balance de potencia:

𝑃𝑔𝑖 − 𝑃𝑑𝑖 + ∑ 𝑃𝑓𝑗𝑖 − ∑ 𝑃𝑓𝑖𝑗 = 0 Ɐ 𝑖 ∈ Ω𝐵 (2)

𝑗𝑖∈Ω𝐿 𝑖𝑗∈Ω𝐿

Flujos de potencia activa:

(𝜃𝑖 − 𝜃𝑗 )
𝑃𝑓𝑖𝑗 = ∗ (𝑛0𝑖𝑗 + 𝑛𝑖𝑗 ) Ɐ 𝑖𝑗 ∈ Ω𝐿 (3)
𝑥𝑙𝑖𝑗

Limites operativos:

−(𝑛0𝑖𝑗 + 𝑛𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 ≤ 𝑃𝑓𝑖𝑗 ≤ (𝑛0𝑖𝑗 + 𝑛𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (4)

Limites operativos de la maquina:

𝑃𝑚𝑖𝑛𝑖 ≤ 𝑃𝑔𝑖 ≤ 𝑃𝑚𝑎𝑥𝑖 Ɐ 𝑖 ∈ Ω𝐵 (5)

Limites de la cantidad de líneas:

0 ≤ 𝑛𝑖𝑗 ≤ 𝑛𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (6)

➢ Modelo de transporte
- Función Objetivo:

𝑀𝑖𝑛 {𝐶𝑇 = ∑ (𝐶𝑙𝑖𝑗 ∗ 𝑛𝑖𝑗 )} (7)

𝑖𝑗∈Ω𝐿

- Restricciones:
Balance de potencia:
 𝑃𝑔𝑖 − 𝑃𝑑𝑖 + ∑ 𝑃𝑓𝑗𝑖 − ∑ 𝑃𝑓𝑖𝑗 = 0 Ɐ 𝑖 ∈ Ω𝐵 (8)
𝑗𝑖∈Ω𝐿 𝑖𝑗∈Ω𝐿

Límites operativos:

−(𝑛0𝑖𝑗 + 𝑛𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 ≤ 𝑃𝑓𝑖𝑗 ≤ (𝑛0𝑖𝑗 + 𝑛𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (9)

Límites operativos de la maquina:

𝑃𝑚𝑖𝑛𝑖 ≤ 𝑃𝑔𝑖 ≤ 𝑃𝑚𝑎𝑥𝑖 Ɐ 𝑖 ∈ Ω𝐵 (10)

Límites de la cantidad de líneas:

0 ≤ 𝑛𝑖𝑗 ≤ 𝑛𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (11)

➢ Modelo hibrido
- Función Objetivo:

𝑀𝑖𝑛 {𝐶𝑇 = ∑ (𝐶𝑙𝑖𝑗 ∗ 𝑛𝑖𝑗 )} (12)

𝑖𝑗∈Ω𝐿

- Restricciones:
Balance de potencia:

𝑃𝑔𝑖 − 𝑃𝑑𝑖 + ∑ 𝑃𝑓𝑗𝑖 + ∑ 𝑃𝑓0𝑗𝑖 − ∑ 𝑃𝑓𝑖𝑗 − ∑ 𝑃𝑓0𝑖𝑗 = 0 Ɐ 𝑖 ∈ Ω𝐵 (13)

𝑗𝑖∈Ω𝐿 𝑗𝑖∈Ω𝐿 𝑖𝑗∈Ω𝐿 𝑖𝑗∈Ω𝐿

Flujos de potencia activa inicial:

(𝜃𝑖 − 𝜃𝑗 )
𝑃𝑓0𝑖𝑗 = ∗ (𝑛0𝑖𝑗 ) Ɐ 𝑖𝑗 ∈ Ω𝐿 (14)
𝑥𝑙𝑖𝑗

Limites operativos inicial:

−(𝑛0𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 ≤ 𝑃𝑓0𝑖𝑗 ≤ (𝑛0𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (15)

Limites operativos inicial:

−(𝑛𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 ≤ 𝑃𝑓𝑖𝑗 ≤ (𝑛𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (16)

Limites operativos de la maquina:

𝑃𝑚𝑖𝑛𝑖 ≤ 𝑃𝑔𝑖 ≤ 𝑃𝑚𝑎𝑥𝑖 Ɐ 𝑖 ∈ Ω𝐵 (17)
Límites de la cantidad de líneas:
0 ≤ 𝑛𝑖𝑗 ≤ 𝑛𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (18)

➢ Modelo disyuntivo
- Función Objetivo:

𝑀𝑖𝑛 {𝐶𝑇 = ∑ (𝐶𝑙𝑖𝑗 ∗ ( ∑ 𝑊𝑖𝑗 ))} (19)

𝑖𝑗∈Ω𝐿 𝑦 ∈Ω𝑌

- Restricciones:
Balance de potencia:

𝑃𝑔𝑖 − 𝑃𝑑𝑖 + ∑ ∑ 𝑃𝑓𝑗𝑖,𝑦 + ∑ 𝑃𝑓0𝑗𝑖 − ∑ ∑ 𝑃𝑓𝑖𝑗,𝑦

𝑗𝑖∈Ω𝐿 𝑦 ∈Ω𝑌 𝑗𝑖∈Ω𝐿 𝑖𝑗∈Ω𝐿 𝑦 ∈Ω𝑌
Ɐ 𝑖 ∈ Ω𝐵 (20)
− ∑ 𝑃𝑓0𝑖𝑗 = 0
𝑖𝑗∈Ω𝐿

Flujos de potencia activa inicial:

(𝜃𝑖 − 𝜃𝑗 )
𝑃𝑓0𝑖𝑗 = ∗ (𝑛0𝑖𝑗 ) Ɐ 𝑖𝑗 ∈ Ω𝐿 (21)
𝑥𝑙𝑖𝑗

Límite de flujos BIGM_minimo:

Ɐ 𝑖𝑗 ∈ Ω𝐿 (22)
−2 ∗ 𝜃𝑚𝑎𝑥 ∗ (1 − 𝑊𝑖𝑗,𝑦 ) ≤ 𝑃𝑓𝑖𝑗 ∗ 𝑥𝑙𝑖𝑗 − (𝜃𝑖 − 𝜃𝑗 )
Ɐ 𝑦 ∈ Ω𝑌

Límite de flujos BIGM_maximo:

Ɐ 𝑖𝑗 ∈ Ω𝐿 (23)
2 ∗ 𝜃𝑚𝑎𝑥 ∗ (1 − 𝑊𝑖𝑗,𝑦 ) ≥ 𝑃𝑓𝑖𝑗 ∗ 𝑥𝑙𝑖𝑗 − (𝜃𝑖 − 𝜃𝑗 )
Ɐ 𝑦 ∈ Ω𝑌

Limites operativos inicial:

−(𝑛0𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 ≤ 𝑃𝑓0𝑖𝑗 ≤ (𝑛0𝑖𝑗 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (24)

Limites operativos planificados:

Ɐ 𝑖𝑗 ∈ Ω𝐿 (25)
−(𝑊𝑖𝑗,𝑦 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗 ≤ 𝑃𝑓𝑖𝑗,𝑦 ≤ (𝑊𝑖𝑗,𝑦 ) ∗ 𝑆𝑚𝑎𝑥𝑖𝑗
Ɐ 𝑦 ∈ Ω𝑌
Limites secuencial:
Ɐ 𝑖𝑗 ∈ Ω𝐿 (26)
(𝑊𝑖𝑗,𝑦 ) ≤ (𝑊𝑖𝑗,𝑦−1 )
Ɐ 𝑦 ∈ Ω𝑌 (𝑦 > 1)

Limites angulares:

−𝜃𝑚𝑎𝑥 ≤ 𝜃𝑖 ≤ 𝜃𝑚𝑎𝑥 Ɐ 𝑖𝑗 ∈ Ω𝐵 (27)

Limites operativos de la maquina:

𝑃𝑚𝑖𝑛𝑖 ≤ 𝑃𝑔𝑖 ≤ 𝑃𝑚𝑎𝑥𝑖 Ɐ 𝑖 ∈ Ω𝐵 (28)

Límites de la cantidad de líneas:

0 ≤ 𝑊𝑖𝑗,𝑦 ≤ 𝑛𝑚𝑎𝑥𝑖𝑗 Ɐ 𝑖𝑗 ∈ Ω𝐿 (29)

Estructura en AMPL
Tabla 3. Modelo DC
CODIGO AMPL MODELO DC
reset;

# INICIAR EL MODELO

# DECALRACION DE CONJUNTOS
set OB;
set OL within 1..9999 cross OB cross OB;
#k i j

#DECLARACION DE PARAMETROS
#parametros del sistema
param Sb;

#parametro de las barras

param Nombre{OB}symbolic; #nombre de la barra
param TB{OB}; #tipo de barra i #Voltaje inicial
en la barra i
param th0{OB}; #Angulo inicial en la barra i
param Pd{OB}; #Potencia activa demandada en la barra i
param Pmax{OB}; #Potencia activa generada max
param Pmin{OB}; #Potencia activa generada min

#parametros de los RAMALES

param rl{OL}; #resistencia del ramal ij
param xl{OL}; #reactancia del ramal ij
param bshl{OL}; #suceptancia shunt del ramal ij
param Smax{OL}; #cargabilidad del ramal del elemento del ramal ij
param g{OL}; #conductancia del ramal ij inverso de la resistencias
param b{OL}; #suceptancia del
ramal ij inverso de la reactancia
param Cl{OL}; #costo de cada linea
param n0{OL}; #cantidad
de lineas que tiene el sistema actualmente
param nmax{OL}; #cantidad
de lineas maximas que puede tener ese ramal

#decalracion de variables
var th{OB}; #angulo en la barra i con respecto a la barra slack
var Pg{OB};
var Pde{OL}; #Flujo de potencia activa desde la barra i hacia la
barra j en el ramal kij
var n{OL}integer;
#cantidad de lineas que va a tener ese ramal

# funcion objetivo
minimize Costos:
sum{(k,i,j) in OL}Cl[k,i,j]*n[k,i,j];

# RESTRICCIONES
#restriccion de balance de potencia en cada barra
subject to balance{i in OB}:
Pg[i]-Pd[i]+sum{(k,j,i) in OL}(Pde[k,j,i])-sum{(k,i,j) in OL}(Pde[k,i,j])=0;

#flujos de potencia activa

subject to flujo_envio{(k,i,j) in OL}:
Pde[k,i,j]=((th[i]-th[j])/xl[k,i,j])*(n0[k,i,j]+n[k,i,j]);

#limites operativos
subject to limites_operativos_minimos{(k,i,j) in OL}:
-(n0[k,i,j]+n[k,i,j])*Smax[k,i,j]<=Pde[k,i,j];

subject to limites_operativos_maximos{(k,i,j) in OL}:

Pde[k,i,j]<=(n0[k,i,j]+n[k,i,j])*Smax[k,i,j];

#limites operativos de la maquina

subject to limites{i in OB}:
Pmin[i]<=Pg[i]<=Pmax[i];

#limites de la cantidad de lineas

subject to limites_l{(k,i,j) in OL}:
0<=n[k,i,j]<=nmax[k,i,j];
 Tabla 4. Modelo de transporte.
CODIGO AMPL MODELO DE TRANSPORTE
reset;

# INICIAR EL MODELO

# DECALRACION DE CONJUNTOS
set OB;
set OL within 1..9999 cross OB cross OB;
#k i j

#DECLARACION DE PARAMETROS
#parametros del sistema
param Sb;

#parametro de las barras

param Nombre{OB}symbolic; #nombre de la barra
param TB{OB}; #tipo de barra i
param Pd{OB}; #Potencia activa demandada en la barra i
param Pmax{OB}; #Potencia activa generada max
param Pmin{OB}; #Potencia activa generada min

#parametros de los RAMALES

param rl{OL}; #resistencia del ramal ij
param xl{OL}; #reactancia del ramal ij
param bshl{OL}; #suceptancia shunt del ramal ij
param Smax{OL}; #cargabilidad del ramal del elemento del ramal ij
param g{OL}; #conductancia del ramal ij inverso de la resistencias
param b{OL}; #suceptancia del
ramal ij inverso de la reactancia
param Cl{OL}; #costo de cada linea
param n0{OL}; #cantidad
de lineas que tiene el sistema actualmente
param nmax{OL}; #cantidad
de lineas maximas que puede tener ese ramal

#decalracion de variables
var Pg{OB};
var Pde{OL};
var Ppa{OL}; #Flujo de potencia activa desde la barra i hacia la
barra j en el ramal kij
var n{OL}integer;
#cantidad de lineas que va a tener ese ramal

# funcion objetivo
minimize Costos:
sum{(k,i,j) in OL}Cl[k,i,j]*n[k,i,j];

# RESTRICCIONES
#restriccion de balance de potencia en cada barra
subject to balance{i in OB}:
Pg[i]-Pd[i]+sum{(k,j,i) in OL}(Ppa[k,j,i])-sum{(k,i,j) in OL}(Pde[k,i,j])=0;

#limites operativos
subject to limites_operativos_minimos_e{(k,i,j) in OL}:
-(n0[k,i,j]+n[k,i,j])*Smax[k,i,j]<=Pde[k,i,j];

subject to limites_operativos_maximos_e{(k,i,j) in OL}:

Pde[k,i,j]<=(n0[k,i,j]+n[k,i,j])*Smax[k,i,j];

subject to limites_operativos_minimos_l{(k,i,j) in OL}:

-(n0[k,i,j]+n[k,i,j])*Smax[k,i,j]<=Ppa[k,i,j];

subject to limites_operativos_maximos_l{(k,i,j) in OL}:

Ppa[k,i,j]<=(n0[k,i,j]+n[k,i,j])*Smax[k,i,j];

#limites operativos de la maquina

subject to limites{i in OB}:
Pmin[i]<=Pg[i]<=Pmax[i];

#limites de la cantidad de lineas

subject to limites_l{(k,i,j) in OL}:
0<=n[k,i,j]<=nmax[k,i,j];

Tabla 5. Modelo Hibrido.

CODIGO AMPL MODELO HIBRIDO
reset;

# INICIAR EL MODELO

# DECALRACION DE CONJUNTOS
set OB;
set OL within 1..9999 cross OB cross OB;
#k i j

#DECLARACION DE PARAMETROS
#parametros del sistema
param Sb;

#parametro de las barras

param Nombre{OB}symbolic; #nombre de la barra
param TB{OB}; #tipo de barra i #Voltaje inicial
en la barra i
param th0{OB}; #Angulo inicial en la barra i
param Pd{OB}; #Potencia activa demandada en la barra i
param Pmax{OB}; #Potencia activa generada max
param Pmin{OB}; #Potencia activa generada min
#parametros de los RAMALES
param rl{OL}; #resistencia del ramal ij
param xl{OL}; #reactancia del ramal ij
param bshl{OL}; #suceptancia shunt del ramal ij
param Smax{OL}; #cargabilidad del ramal del elemento del ramal ij
param g{OL}; #conductancia del ramal ij inverso de la resistencias
param b{OL}; #suceptancia del
ramal ij inverso de la reactancia
param Cl{OL}; #costo de cada linea
param n0{OL}; #cantidad
de lineas que tiene el sistema actualmente
param nmax{OL}; #cantidad
de lineas maximas que puede tener ese ramal

#decalracion de variables
var th{OB}; #angulo en la barra i con respecto a la barra slack
var Pg{OB};
var Pf0{OL}; #Flujo de potencia activa desde la barra i hacia la
barra j en el ramal kij
var Pf{OL};
var n{OL}integer;
#cantidad de lineas que va a tener ese ramal

# funcion objetivo
minimize Costos:
sum{(k,i,j) in OL}Cl[k,i,j]*n[k,i,j];

# RESTRICCIONES
#restriccion de balance de potencia en cada barra
subject to balance{i in OB}:
Pg[i]-Pd[i]+sum{(k,j,i) in OL}(Pf[k,j,i])+sum{(k,j,i) in OL}(Pf0[k,j,i])-sum{(k,i,j) in
OL}(Pf[k,i,j])-sum{(k,i,j) in OL}(Pf0[k,i,j])=0;

#flujos de potencia activa

subject to flujo_envio{(k,i,j) in OL}:
Pf0[k,i,j]=((th[i]-th[j])/xl[k,i,j])*(n0[k,i,j]);

#limites operativos
subject to limites_operativos_iniciales{(k,i,j) in OL}:
-n0[k,i,j]*Smax[k,i,j]<=Pf0[k,i,j]<=n0[k,i,j]*Smax[k,i,j];

subject to limites_operativos_minimos{(k,i,j) in OL}:

-n[k,i,j]*Smax[k,i,j]<=Pf[k,i,j];

subject to limites_operativos_maximos{(k,i,j) in OL}:

Pf[k,i,j]<=n[k,i,j]*Smax[k,i,j];

#limites operativos de la maquina

subject to limites{i in OB}:
Pmin[i]<=Pg[i]<=Pmax[i];

#limites de la cantidad de lineas

subject to limites_l{(k,i,j) in OL}:
0<=n[k,i,j]<=nmax[k,i,j];

Tabla 6. Modelo Disyuntivo.

CODIGO AMPL MODELO DISYUNTIVO
reset;

# INICIAR EL MODELO

# DECALRACION DE CONJUNTOS
set OB;
set OL within 1..9999 cross OB cross OB;
#k i j
set Y=1..3;

#DECLARACION DE PARAMETROS
#parametros del sistema
param Sb;

#parametro de las barras

param Nombre{OB}symbolic; #nombre de la barra
param TB{OB}; #tipo de barra i
param th0{OB}; #Angulo inicial en la barra i
param Pd{OB}; #Potencia activa demandada en la barra i
param Pmax{OB}; #Potencia activa generada max
param Pmin{OB}; #Potencia activa generada min
param thmax;

#parametros de los RAMALES

param rl{OL}; #resistencia del ramal ij
param xl{OL}; #reactancia del ramal ij
param bshl{OL}; #suceptancia shunt del ramal ij
param Smax{OL}; #cargabilidad del ramal del elemento del ramal ij
param g{OL}; #conductancia del ramal ij inverso de la resistencias
param b{OL}; #suceptancia del
ramal ij inverso de la reactancia
param Cl{OL}; #costo de cada linea
param n0{OL}; #cantidad
de lineas que tiene el sistema actualmente
param nmax{OL}; #cantidad
de lineas maximas que puede tener ese ramal

#decalracion de variables
var W{OL,Y}binary;
var th{OB}; #angulo en la barra i con respecto a la barra slack
var Pg{OB};
var Pf0{OL}; #Flujo de potencia activa desde la barra i hacia la
barra j en el ramal kij
var Pf{OL,Y};
var n{OL}integer;
#cantidad de lineas que va a tener ese ramal

# funcion objetivo
minimize Costos:
sum{(k,i,j) in OL}Cl[k,i,j]*sum{y in Y}W[k,i,j,y];

# RESTRICCIONES
#restriccion de balance de potencia en cada barra
subject to balance{i in OB}:
Pg[i]-Pd[i]+sum{(k,j,i) in OL}(sum{y in Y}Pf[k,j,i,y])+sum{(k,j,i) in OL}(Pf0[k,j,i])-
sum{(k,i,j) in OL}(sum{y in Y}(Pf[k,i,j,y]))-sum{(k,i,j) in OL}(Pf0[k,i,j])=0;

#flujos de potencia activa

subject to flujo_envio{(k,i,j) in OL}:
Pf0[k,i,j]=((th[i]-th[j])/xl[k,i,j])*(n0[k,i,j]);

#flujos de potencia de las lineas proyectadas

#subject to flujo_proyectado{(k,i,j) in OL, y in Y}:
#Pf[k,i,j,y]=((th[i]-th[j])/xl[k,i,j])*(W[k,i,j,y]);

#restriccion de BIGM
subject to BIGM_maximo{(k,i,j) in OL, y in Y}:
(Pf[k,i,j,y]*xl[k,i,j])-(th[i]-th[j])<=2*thmax*(1-W[k,i,j,y]);

subject to BIGM_minimo{(k,i,j) in OL, y in Y}:

(Pf[k,i,j,y]*xl[k,i,j])-(th[i]-th[j])>=-(2*thmax)*(1-W[k,i,j,y]);

#limites operativos
subject to limites_operativos_iniciales{(k,i,j) in OL}:
-n0[k,i,j]*Smax[k,i,j]<=Pf0[k,i,j]<=n0[k,i,j]*Smax[k,i,j];

subject to limites_operativos_minimos{(k,i,j) in OL, y in Y}:

-W[k,i,j,y]*Smax[k,i,j]<=Pf[k,i,j,y];

subject to limites_operativos_maximos{(k,i,j) in OL,y in Y}:

Pf[k,i,j,y]<=W[k,i,j,y]*Smax[k,i,j];

#Limite secuencial
subject to secuencial{(k,i,j) in OL, y in Y: y>1}:
W[k,i,j,y]<=W[k,i,j,y-1];

#otras restricciones
subject to limites_angular{i in OB}:
-thmax<=th[i]<=thmax;
#limites operativos de la maquina
subject to limites_g{i in OB}:
Pmin[i]<=Pg[i]<=Pmax[i];

#limites de la cantidad de lineas

subject to limites_l{(k,i,j) in OL}:
sum {y in Y}W[k,i,j,y]<=nmax[k,i,j];

Resultados
Modelo DC
Knitro 13.0.1: Locally optimal or satisfactory solution.
objective 130; optimality gap -1.13e-08
43 nodes; 60 subproblem solves
suffix feaserror OUT;
suffix opterror OUT;
suffix numfcevals OUT;
suffix numiters OUT;
suffix incumbent OUT;
suffix relaxbnd OUT;
Costo = 130
n :=
1120
2130
3140
4150
5160
6231
7240
8250
9261
10 3 4 0
11 3 5 1
12 3 6 0
13 4 5 0
14 4 6 2
15 5 6 0
;
Pij[k,i,j]*Sbase :=
1 1 2 -11.6948
2130
3 1 4 -10.275
4 1 5 69.445
5160
6 2 3 -148.968
7 2 4 -10.7345
8250
9 2 6 -91.9918
10 3 4 0
11 3 5 170.555
12 3 6 0
13 4 5 0
14 4 6 -181.01
15 5 6 0
;
Pg[i]Sbase [] :=
1 127.475
20
3 359.523
40
50
6 273.001
;
th[i]180/3.14159 [] :=
10
2 2.68026
3 9.07133
4 3.53228
5 -7.95782
6 16.6163

Modelo de transporte
CPLEX 20.1.0.0: optimal integer solution; objective 110
20 MIP simplex iterations
0 branch-and-bound nodes
Costo = 110
n :=
1120
2130
3140
4150
5160
6230
7240
8250
9263
10 3 4 0
11 3 5 1
12 3 6 0
13 4 5 0
14 4 6 0
15 5 6 0
;
Pij[k,i,j]*Sbase :=
1 1 2 -30
2130
3 1 4 60
4 1 5 40
5160
6 2 3 -70
7 2 4 100
8250
9 2 6 -300
10 3 4 0
11 3 5 200
12 3 6 0
13 4 5 0
14 4 6 0
15 5 6 0
;
Pg[i]Sbase [] :=
1 150
20
3 310
40
50
6 300

Modelo hibrido
CPLEX 20.1.0.0: optimal integer solution; objective 110
37 MIP simplex iterations
0 branch-and-bound nodes
: Pf Pf0 n n0 :=
1 1 2 0 0.2 0 1
2 13 0 0 0 0
3 14 0 0 0 1
4 1 5 0 0.4 0 1
5 16 0 0 0 0
6 2 3 0 -1 0 1
7 2 4 0 -0.2 0 1
8 25 0 0 0 0
9 2 6 -1 0 1 0
10 3 4 0 0 0 0
11 3 5 1 1 1 1
12 3 6 0 0 0 0
13 4 5 0 0 0 0
14 4 6 -1.8 0 2 0
15 5 6 0 0 0 0
;

Costos = 110

Modelo disyuntivo
CPLEX 20.1.0.0: optimal integer solution; objective 110
151 MIP simplex iterations
0 branch-and-bound nodes
Pf [*,1,*,1]
: 2 3 4 5 6 :=
1 0 . . . .
2 . 0 . . .
3 . . 0 . .
4 . . . 0 .
5 . . . . 0
 [*,1,*,2]
: 2 3 4 5 6 :=
1 0 . . . .
2 . 0 . . .
3 . . 0 . .
4 . . . 0 .
5 . . . . 0

[*,1,*,3]
: 2 3 4 5 6 :=
1 0 . . . .
2 . 0 . . .
3 . . 0 . .
4 . . . 0 .
5 . . . . 0

[*,2,*,1]
: 3 4 5 6 :=
6 0 . . .
7 . 0 . .
8 . . 0 .
9 . . . 0

[*,2,*,2]
: 3 4 5 6 :=
6 0 . . .
7 . 0 . .
8 . . 0 .
9 . . . 0

[*,2,*,3]
: 3 4 5 6 :=
6 0 . . .
7 . 0 . .
8 . . 0 .
9 . . . 0

[*,3,*,1]
: 4 5 6 :=
10 0 . .
11 . 0.857576 .
12 . . 0

[*,3,*,2]
: 4 5 6 :=
10 0 . .
11 . 0 .
12 . . 0

[*,3,*,3]
: 4 5 6 :=
10 0 . .
11 . 0 .
12 . . 0
 [*,4,*,1]
: 5 6 :=
13 0 .
14 . -0.99798

[*,4,*,2]
: 5 6 :=
13 0 .
14 . -0.99798

[*,4,*,3]
: 5 6 :=
13 0 .
14 . -0.99798

[15,5,6,1] 0

[15,5,6,2] 0

[15,5,6,3] 0
;

Pf0 :=
1 12 0.409091
2 13 0
3 14 -0.393939
4 15 0.684848
5 16 0
6 23 -0.990909
7 24 -1
8 25 0
9 26 0
10 3 4 0
11 3 5 0.857576
12 3 6 0
13 4 5 0
14 4 6 0
15 5 6 0
;

W [*,1,*,1]
: 2 3 4 5 6 :=
1 0 . . . .
2 . 0 . . .
3 . . 0 . .
4 . . . 0 .
5 . . . . 0

[*,1,*,2]
: 2 3 4 5 6 :=
1 0 . . . .
2 . 0 . . .
3 . . 0 . .
4 . . . 0 .
5 . . . . 0

[*,1,*,3]
: 2 3 4 5 6 :=
1 0 . . . .
2 . 0 . . .
3 . . 0 . .
4 . . . 0 .
5 . . . . 0

[*,2,*,1]
: 3 4 5 6 :=
6 0 . . .
7 . 0 . .
8 . . 0 .
9 . . . 0

[*,2,*,2]
: 3 4 5 6 :=
6 0 . . .
7 . 0 . .
8 . . 0 .
9 . . . 0

[*,2,*,3]
: 3 4 5 6 :=
6 0 . . .
7 . 0 . .
8 . . 0 .
9 . . . 0

[*,3,*,1]
: 4 5 6 :=
10 0 . .
11 . 1 .
12 . . 0

[*,3,*,2]
: 4 5 6 :=
10 0 . .
11 . 0 .
12 . . 0

[*,3,*,3]
: 4 5 6 :=
10 0 . .
11 . 0 .
12 . . 0

[*,4,*,1]
: 5 6 :=
13 0 .
14 . 1

[*,4,*,2]
 : 5 6 :=
13 0 .
14 . 1

[*,4,*,3]
: 5 6 :=
13 0 .
14 . 1

[15,5,6,1] 0

[15,5,6,2] 0

[15,5,6,3] 0
;

n0 :=
1 12 1
2 13 0
3 14 1
4 15 1
5 16 0
6 23 1
7 24 1
8 25 0
9 26 0
10 3 4 0
11 3 5 1
12 3 6 0
13 4 5 0
14 4 6 0
15 5 6 0
;

Costos = 110

Análisis de resultados
Tabla 7. Resultados obtenidos por la simulación.
Costo Total
∆%
Modelo DC 130
Modelo de transporte 110 15.38%
Modelo hibrido 110 0%
Modelo Disyuntivo 110 0%

Como se puedo observar el modelo DC es el único que tiene un costo de 130, esto
utilizando el solver knitro en NEMO, el resto de las modelos tienen el mismo resultado
puesto que son modelos lineales con puntos de vista diferentes en la resolución de flujos,
de esta manera se podría mencionar que no existe diferencia en los costos de estos tres
modelos, su costos es de 110.
CONCLUSIONES:
1. Como se pudo observar en la simulación el modelo hibrido, disyuntivo y de
transporte tuvieron el mismo resultado, puesto que los dos primeros modelos son
una adaptación del modelo de transporte, pero más completos, pero aun así los
resultado son los mismo ya que el sistema es pequeño y las decisiones son
aproximadamente las mismas.
2. Estos modelos nos ayudan a ir complementando el sistema DC para el estudio
exclusivo de la implementación de nuevas líneas en sectores ya previsto, como
seria el flujo en esa línea y cuantas se deben implementar en ese periodo de
tiempo.