Hello everyone, I need your help. I am currently modelling a distribution network (13 buses) using Julia software via PowerModelDistribution with a MATHEMATICAL model. The problem I am encountering is that when calculating the power flow, the voltages do not vary and remain at 1.0 PU, except for those of the Slack Bus. I have tried everything and realise that my network is disconnected, even though the lines and other components of the network are well defined. For this reason, I would like to know if anyone has ever had this kind of problem and how I can solve it.
I remain at your disposal for further information on the problem.
Thank you.
Hello and welcome to the forum.
Please start here:
For the average person (I mean not domain experts) to be able to help, we would need to first see some code.
###### Here is my main code where I connect the network and perform the calculations.
using PowerModelsDistribution
using Ipopt
using Printf
const PMD = PowerModelsDistribution
cd(@__DIR__)
bus_dict = include(joinpath(@__DIR__,"Bus.jl"))
load_dict = include(joinpath(@__DIR__,"Lasten.jl"))
line_dict = include(joinpath(@__DIR__, "Linien.jl"))
math_data = Dict{String, Any}(
"bus" => bus_dict["bus"],
"load" => load_dict["load"],
"branch" => line_dict["branch"],
"source" => Dict{String, Any}(
"1" => Dict{String, Any}("connections" => [1, 2, 3],
"bus" => "632",
"index" => 1,
"status" => 1,
"vm" => [1.0210, 1.0420, 1.0174],
"va" => [-2.49*(pi/180), -121.72*(pi/180), 117.83*(pi/180)]
)
),
"gen" => Dict{String, Any}(
"1" => Dict{String, Any}(
"connections" => [1, 2, 3],
"gen_bus" => "632", "index" => 1,
"gen_status" => 1,
"gen_type" => 3,
#"vm" => [ 1.0210, 1.0420, 1.0174],
#"va" => [ -2.49*(pi /180), -121.72*(pi /180), 117.83*(pi /180)], #radiant winckel
"status" =>1,
#"index" =>1,
"pg" => [5.0, 5.0, 5.0],
"qg" => [0.1, 0.1, 0.1],
"pmax" => [10.0, 10.0, 10.0],
"pmin" => [0.0, 0.0, 0.0],
"qmax" => [10.0, 10.0, 10.0],
"qmin" => [-10.0, -10.0, -10.0]
)
),
"switch" => Dict{String, Any}(
"1" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "692",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"state" => 1,
"status" => 1,
"index" => 1,
"rs" => [0.001, 0.001, 0.001],
"xs" => [0.0, 0.0, 0.0],
"angmin" => [-1.0, -1.0, -1.0], "argmax"=>[1.0,1.0,1.0]
)
),
"transformer" => Dict{String, Any}(
"1" => Dict{String, Any}(
"f_bus" => "633",
"t_bus" => "634",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"configuration" => PMD.WYE,
"tm_nom" => 1.0,
"tm_set" => [1.0, 1.0, 1.0],
"polarity" => 1,
"vbase" => [4.16, 0.48],
"xsc" => [0.8, 0.8, 0.8],
"rw" => [0.1, 0.1, 0.1],
"index" => 1,
"status" => 1,
"tm_status" => 1
)
),
"settings" => Dict{String, Any}("sbase" => 1.0, "vbase" => 4.16),
"per_unit" => true, "data_model" => PMD.MATHEMATICAL
)
###############################################################################
println("Nombre de bus chargés : ", length(math_data["bus"]))
println("Nombre de branches chargées : ", length(math_data["branch"]))
println("Nombre de charges chargées : ", length(math_data["load"]))
println("Nombre de générateurs chargés : ", length(get(math_data, "gen", Dict())))
println("Nombre de transformateur chargés : ", length(get(math_data, "transformer", Dict())))
println("Nombre de Breacker chargés : ", length(get(math_data, "switch", Dict())))
id_tr = first(keys(math_data["transformer"]))
println("Transformateur reliant le bus $(math_data["transformer"][id_tr]["f_bus"]) au bus $(math_data["transformer"][id_tr]["t_bus"])")
println("Branches reconnues par PMD : ", length(math_data["branch"]))
println("Charges reconnues par PMD : ", length(math_data["load"]))
###
opt = optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6, "print_level"=>0)
resultat_pf = PMD.solve_mc_pf(math_data, PMD.ACPUPowerModel, opt)
# Vérifie immédiatement si la branche 1 existe dans la solution
if haskey(resultat_pf["solution"], "branch") && haskey(resultat_pf["solution"]["branch"], "1")
flux = resultat_pf["solution"]["branch"]["1"]["pf"]
println(">>> Flux de puissance sur la branche 1 (632-633) : ", flux)
else
println(">>> ALERTE : La branche 1 est INEXISTANTE dans la solution finale !")
end
# 3. AFFICHAGE ET CALCUL DES INJECTIONS
if resultat_pf["termination_status"] in [PMD.LOCALLY_SOLVED, PMD.ALMOST_LOCALLY_SOLVED]
# On récupÚration des données de sol
ma_sol = resultat_pf["solution"]
sbase_kw = get(math_data["settings"], "sbase", 1.0) * 1000.0
println("\n" * "â"^95)
println(" RAPPORT DE FLUX DE PUISSANCE : ĂTAT DES BUS")
println("â"^95)
@printf("%-5s | %-5s | %-10s | %-10s | %-15s | %-15s\n",
"BUS", "Ph", "V (pu)", "Angle (°)", "P_net (kW)", "Q_net (kVAR)")
println("-"^95)
# Tri des bus par ordre numérique
bus_ids = sort(collect(keys(ma_sol["bus"])), by=x->parse(Int, x))
for b_id in bus_ids
# Données de tension du bus
b_sol = ma_sol["bus"][b_id]
vm = b_sol["vm"]
va = rad2deg.(b_sol["va"])
#CALCUL DES INJECTIONS NETTES PAR PHASE...
p_net = zeros(length(vm))
q_net = zeros(length(vm))
# production des Générateurs
for (g_id, g_res) in get(ma_sol, "gen", Dict())
if string(math_data["gen"][g_id]["gen_bus"]) == b_id
pg_vals = g_res["pg"]
qg_vals = g_res["qg"]
for i in 1:min(length(p_net), length(pg_vals))
p_net[i] += pg_vals[i]
q_net[i] += qg_vals[i]
end
end
end
# 2. Soustraire la consommation (Charges)
for (l_id, l_info) in get(math_data, "load", Dict())
if string(l_info["load_bus"]) == b_id
pd_vals = l_info["pd"]
qd_vals = l_info["qd"]
for i in 1:min(length(p_net), length(pd_vals))
p_net[i] -= pd_vals[i]
q_net[i] -= qd_vals[i]
end
end
end
# Conversion finale en unités physiques (kW)
p_final = p_net .* sbase_kw
q_final = q_net .* sbase_kw
# Affichage de chaque phase pour le bus
for i in 1:length(vm)
@printf("%-5s | %d | %12f | %9.2f | %12.2f | %12.2f\n",
b_id, i, vm[i], va[i], p_final[i], q_final[i])
end
println("-"^95)
end
println("â"^95)
else
println(" Le solveur n'a pas convergé. Statut : ", resultat_pf["termination_status"])
end
############# Here I have defined the different branches of the network.
#Erstellung der Linien
# 09 branch
using LinearAlgebra
# Constantes de conversion
mi2Km = 1.60934
ft2km = 0.0003048
ms2F = 1/2/pi/60*1e-6
# Per unit umrechnen: Z_base = (V_base)^2/S_base
V_base = 4.16
Z_base = (V_base)^2 /1.0
# Configuration 601 (Triphasée - Standard)
R_601 = [0.3465 0.1560 0.1580; 0.1560 0.3375 0.1535; 0.1580 0.1535 0.3414]./Z_base
X_601 = [1.0179 0.5017 0.4236; 0.5017 1.0478 0.3849; 0.4236 0.3849 1.0348] ./Z_base
B_601 = [6.2998 -1.9958 -1.2595; -1.9958 5.9597 -0.7417; -1.2595 -0.7417 5.6386].*Z_base
# Convert for SPS
R_601_sps = R_601/mi2Km
L_601_sps = X_601/mi2Km/2/pi/60
C_601_sps = B_601/mi2Km*ms2F
# Configuration 602 (Triphasée)
R_602 = [0.7526 0.1580 0.1560; 0.1580 0.7475 0.1535; 0.1560 0.1535 0.7436]./Z_base
X_602 = [1.1814 0.4236 0.5017; 0.4236 1.1983 0.3849; 0.5017 0.3849 1.2112]./Z_base
B_602 = [5.6990 -1.0817 -1.6905; -1.0817 5.1795 -0.6588; -1.6905 -0.6588 5.4246].*Z_base
# Convert for SPS
R_602_sps = R_602/mi2Km
L_602_sps = X_602/mi2Km/2/pi/60
C_602_sps = B_602/mi2Km*ms2F
# Configuration 603
R_603 = [0 0 0 ; 0 1.3294 0.2066; 0 0.2066 1.3238]./Z_base
X_603 = [0 0 0; 0 1.3471 0.4591; 0 0.4591 1.3569] ./Z_base
B_603 = [0 0 0; 0 4.7097 -0.8999;0 -0.8999 4.6658].*Z_base
# Convert for SPS
R_603_sps = R_603/mi2Km
L_603_sps = X_603/mi2Km/2/pi/60
C_603_sps = B_603/mi2Km*ms2F
# Configuration 604
R_604 = [1.3238 0 0.2066;0 0 0; 0.2066 0 1.3294]./Z_base
X_604 = [1.3569 0 0.4591;0 0 0; 0.4591 0 1.3471]./Z_base
B_604 = [4.6658 0 -0.8999;0 0 0; -0.8999 0 4.7097].*Z_base
# Convert for SPS
R_604_sps = R_604/mi2Km
L_604_sps = X_604/mi2Km/2/pi/60
C_604_sps = B_604/mi2Km*ms2F
#Configuration 605
R_605 = [0 0 0;0 0 0;0 0 1.3292]./Z_base
X_605 = [0 0 0;0 0 0;0 0 1.3475]./Z_base
B_605 = [0 0 0;0 0 0;0 0 4.5193].*Z_base
# Convert for SPS
R_605_sps = R_605/mi2Km
L_605_sps = X_605/mi2Km/2/pi/60
C_605_sps = B_605/mi2Km*ms2F
# Configuration 606
R_606 = [0.7982 0.3192 0.2849; 0.3192 0.7891 0.3192; 0.2849 0.3192 0.7982]./Z_base
X_606 = [0.4463 0.0328 0.0143; 0.0328 0.4041 0.0328; 0.0143 0.0328 0.4463]./Z_base
B_606 = [96.8897 -1e-6 -1e-6; -1e-6 96.8897 -1e-6; -1e-6 -1e-6 96.8897].*Z_base
# Convert for SPS
R_606_sps = R_606/mi2Km
L_606_sps = X_606/mi2Km/2/pi/60
C_606_sps = B_606/mi2Km*ms2F
# Configuration 607
R_607 = [1.3425 0 0;0 0 0;0 0 0]./Z_base
X_607 = [0.5124 0 0;0 0 0;0 0 0]./Z_base
B_607= [88.9912 0 0;0 0 0;0 0 0].*Z_base
# Convert for SPS
R_607_sps = R_607/mi2Km
L_607_sps = X_607/mi2Km/2/pi/60
C_607_sps = B_607/mi2Km*ms2F
# --- DĂFINITION DES LIGNES ---
Dict(
"branch" => Dict{String, Any}(
"1" => Dict{String, Any}( #632_633
"f_bus" => "632",
"t_bus" => "633",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_602_sps.*500*ft2km,
"br_x" => X_602.*500*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_602.*500*ft2km/2, "b_to" => B_602.*500*ft2km/2,
"br_status" => 1, "index" => 1,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"2" => Dict{String, Any}( #632_645
"f_bus" => "632",
"t_bus" => "645",
"f_connections" => [2, 3],
"t_connections" => [2, 3],
"br_r" => R_603_sps[2:3,2:3].*500*ft2km,
"br_x" => X_603[2:3,2:3].*500*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_603[2:3,2:3].*500*ft2km./2, "b_to" => B_603[2:3,2:3].*500*ft2km./2,
"br_status" => 1, "index" => 2,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"3" => Dict{String, Any}( #645_646
"f_bus" => "645",
"t_bus" => "646",
"f_connections" => [2, 3],
"t_connections" => [2, 3],
"br_r" => R_603_sps[2:3,2:3].*300*ft2km,
"br_x" => X_603[2:3,2:3].*300*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_603[2:3,2:3].*300*ft2km./2, "b_to" => B_603[2:3,2:3].*300*ft2km./2,
"br_status" => 1, "index" => 3,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"4" => Dict{String, Any}( #632_671
"f_bus" => "632",
"t_bus" => "671",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_601_sps.*2000*ft2km,
"br_x" => X_601.*2000*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_601.*2000*ft2km/2, "b_to" => B_601.*2000*ft2km/2,
"br_status" => 1, "index" => 4,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"5" => Dict{String, Any}( #692_675
"f_bus" => "692",
"t_bus" => "675",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_606_sps.*500*ft2km,
"br_x" => X_606.*500*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_606.*500*ft2km/2, "b_to" => B_606.*500*ft2km/2,
"br_status" => 1, "index" => 5,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"6" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "680",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_601_sps.*1000*ft2km,
"br_x" => X_601.*1000*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_601.*1000*ft2km/2, "b_to" => B_601.*1000*ft2km/2,
"br_status" => 1, "index" => 6,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"7" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "684",
"f_connections" => [1, 3],
"t_connections" => [1, 3],
"br_r" => R_604_sps[[1, 3], [1, 3]].*300*ft2km,
"br_x" => X_604[[1, 3], [1, 3]].*300*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_604[[1, 3], [1, 3]].*300*ft2km ./ 2, "b_to" => B_604[[1, 3], [1, 3]].*300*ft2km./ 2,
"br_status" => 1, "index" => 7,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"8" => Dict{String, Any}(
"f_bus" => "684",
"t_bus" => "652",
"f_connection" => [1],
"t_connection" => [1],
"br_r" => R_604_sps[1: 1, 1: 1].*800*ft2km,
"br_x" => X_604[1: 1, 1: 1].*800*ft2km,
"g_fr" => zeros(1,1), "g_to" => zeros(1,1),
"b_fr" => B_604[1: 1, 1: 1].*800*ft2km ./ 2, "b_to" => B_604[1: 1, 1: 1].*800*ft2km ./ 2,
"br_status" => 1, "index" => 8,
"angmin" => [-3.14], "angmax" => [3.14]
),
"9" => Dict{String, Any}(
"f_bus" => "684",
"t_bus" => "611",
"f_connection" => [3],
"t_connection" => [3],
"br_r" => R_605_sps[3: 3, 3: 3].*300*ft2km,
"br_x" => X_605[3: 3, 3: 3].*300*ft2km,
"g_fr" => zeros(1,1), "g_to" => zeros(1,1),
"b_fr" => B_605[3: 3, 3: 3].*300*ft2km ./ 2, "b_to" => B_605[3: 3, 3: 3].*300*ft2km ./ 2,
"br_status" => 1, "index" => 9,
"angmin" => [-3.14], "angmax" => [3.14]
),
)
)
please format your code correctly, this is also explained in Please read: make it easier to help you
You donât have to make a new post, just edit the original one.
My code is structured into four main parts: one for buses, lines, loads, and the main code.
using PowerModelsDistribution
using Ipopt
using Printf
const PMD = PowerModelsDistribution
cd(@__DIR__)
bus_dict = include(joinpath(@__DIR__,"Bus.jl"))
load_dict = include(joinpath(@__DIR__,"Lasten.jl"))
line_dict = include(joinpath(@__DIR__, "Linien.jl"))
math_data = Dict{String, Any}(
"bus" => bus_dict["bus"],
"load" => load_dict["load"],
"branch" => line_dict["branch"],
"source" => Dict{String, Any}(
"1" => Dict{String, Any}("connections" => [1, 2, 3],
"bus" => "632",
"index" => 1,
"status" => 1,
"vm" => [1.0210, 1.0420, 1.0174],
"va" => [-2.49*(pi/180), -121.72*(pi/180), 117.83*(pi/180)]
)
),
"gen" => Dict{String, Any}(
"1" => Dict{String, Any}(
"connections" => [1, 2, 3],
"gen_bus" => "632", "index" => 1,
"gen_status" => 1,
"gen_type" => 3,
#"vm" => [ 1.0210, 1.0420, 1.0174],
#"va" => [ -2.49*(pi /180), -121.72*(pi /180), 117.83*(pi /180)], #radiant winckel
"status" =>1,
#"index" =>1,
"pg" => [5.0, 5.0, 5.0],
"qg" => [0.1, 0.1, 0.1],
"pmax" => [10.0, 10.0, 10.0],
"pmin" => [0.0, 0.0, 0.0],
"qmax" => [10.0, 10.0, 10.0],
"qmin" => [-10.0, -10.0, -10.0]
)
),
"switch" => Dict{String, Any}(
"1" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "692",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"state" => 1,
"status" => 1,
"index" => 1,
"rs" => [0.001, 0.001, 0.001],
"xs" => [0.0, 0.0, 0.0],
"angmin" => [-1.0, -1.0, -1.0], "argmax"=>[1.0,1.0,1.0]
)
),
"transformer" => Dict{String, Any}(
"1" => Dict{String, Any}(
"f_bus" => "633",
"t_bus" => "634",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"configuration" => PMD.WYE,
"tm_nom" => 1.0,
"tm_set" => [1.0, 1.0, 1.0],
"polarity" => 1,
"vbase" => [4.16, 0.48],
"xsc" => [0.8, 0.8, 0.8],
"rw" => [0.1, 0.1, 0.1],
"index" => 1,
"status" => 1,
"tm_status" => 1
)
),
"settings" => Dict{String, Any}("sbase" => 1.0, "vbase" => 4.16),
"per_unit" => true, "data_model" => PMD.MATHEMATICAL
)
###############################################################################
println("Nombre de bus chargés : ", length(math_data["bus"]))
println("Nombre de branches chargées : ", length(math_data["branch"]))
println("Nombre de charges chargées : ", length(math_data["load"]))
println("Nombre de générateurs chargés : ", length(get(math_data, "gen", Dict())))
println("Nombre de transformateur chargés : ", length(get(math_data, "transformer", Dict())))
println("Nombre de Breacker chargés : ", length(get(math_data, "switch", Dict())))
id_tr = first(keys(math_data["transformer"]))
println("Transformateur reliant le bus $(math_data["transformer"][id_tr]["f_bus"]) au bus $(math_data["transformer"][id_tr]["t_bus"])")
println("Branches reconnues par PMD : ", length(math_data["branch"]))
println("Charges reconnues par PMD : ", length(math_data["load"]))
###
opt = optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6, "print_level"=>0)
resultat_pf = PMD.solve_mc_pf(math_data, PMD.ACPUPowerModel, opt)
# Verification de la branche 1
if haskey(resultat_pf["solution"], "branch") && haskey(resultat_pf["solution"]["branch"], "1")
flux = resultat_pf["solution"]["branch"]["1"]["pf"]
println(">>> Flux de puissance sur la branche 1 (632-633) : ", flux)
else
println(">>> ALERTE : La branche 1 est INEXISTANTE dans la solution finale !")
end
# 3. AFFICHAGE ET CALCUL DES INJECTIONS
if resultat_pf["termination_status"] in [PMD.LOCALLY_SOLVED, PMD.ALMOST_LOCALLY_SOLVED]
# On récupÚration des données de sol
ma_sol = resultat_pf["solution"]
sbase_kw = get(math_data["settings"], "sbase", 1.0) * 1000.0
println("\n" * "â"^95)
println(" RAPPORT DE FLUX DE PUISSANCE : ĂTAT DES BUS")
println("â"^95)
@printf("%-5s | %-5s | %-10s | %-10s | %-15s | %-15s\n",
"BUS", "Ph", "V (pu)", "Angle (°)", "P_net (kW)", "Q_net (kVAR)")
println("-"^95)
# Tri des bus par ordre numérique
bus_ids = sort(collect(keys(ma_sol["bus"])), by=x->parse(Int, x))
for b_id in bus_ids
# Données de tension du bus
b_sol = ma_sol["bus"][b_id]
vm = b_sol["vm"]
va = rad2deg.(b_sol["va"])
#CALCUL DES INJECTIONS NETTES PAR PHASE...
p_net = zeros(length(vm))
q_net = zeros(length(vm))
# production des Générateurs
for (g_id, g_res) in get(ma_sol, "gen", Dict())
if string(math_data["gen"][g_id]["gen_bus"]) == b_id
pg_vals = g_res["pg"]
qg_vals = g_res["qg"]
for i in 1:min(length(p_net), length(pg_vals))
p_net[i] += pg_vals[i]
q_net[i] += qg_vals[i]
end
end
end
# 2. Soustraire la consommation (Charges)
for (l_id, l_info) in get(math_data, "load", Dict())
if string(l_info["load_bus"]) == b_id
pd_vals = l_info["pd"]
qd_vals = l_info["qd"]
for i in 1:min(length(p_net), length(pd_vals))
p_net[i] -= pd_vals[i]
q_net[i] -= qd_vals[i]
end
end
end
# Conversion finale en unités physiques (kW)
p_final = p_net .* sbase_kw
q_final = q_net .* sbase_kw
# Affichage de chaque phase pour le bus
for i in 1:length(vm)
@printf("%-5s | %d | %12f | %9.2f | %12.2f | %12.2f\n",
b_id, i, vm[i], va[i], p_final[i], q_final[i])
end
println("-"^95)
end
println("â"^95)
else
println("â Le solveur n'a pas convergĂ©. Statut : ", resultat_pf["termination_status"])
end
#Erstellung der Lasten / Load
# 15 Lasten
using LinearAlgebra
# Constantes de conversion
mi2Km = 1.60934
ft2km = 0.0003048
ms2F = 1/2/pi/60*1e-6
# Per unit umrechnen: Z_base = (V_base)^2/S_base
V_base = 4.16
Z_base = (V_base)^2 /1.0
# Configuration 601 (Triphasée - Standard)
R_601 = [0.3465 0.1560 0.1580; 0.1560 0.3375 0.1535; 0.1580 0.1535 0.3414]./Z_base
X_601 = [1.0179 0.5017 0.4236; 0.5017 1.0478 0.3849; 0.4236 0.3849 1.0348] ./Z_base
B_601 = [6.2998 -1.9958 -1.2595; -1.9958 5.9597 -0.7417; -1.2595 -0.7417 5.6386].*Z_base
# Convert for SPS
R_601_sps = R_601/mi2Km
L_601_sps = X_601/mi2Km/2/pi/60
C_601_sps = B_601/mi2Km*ms2F
# Configuration 602 (Triphasée)
R_602 = [0.7526 0.1580 0.1560; 0.1580 0.7475 0.1535; 0.1560 0.1535 0.7436]./Z_base
X_602 = [1.1814 0.4236 0.5017; 0.4236 1.1983 0.3849; 0.5017 0.3849 1.2112]./Z_base
B_602 = [5.6990 -1.0817 -1.6905; -1.0817 5.1795 -0.6588; -1.6905 -0.6588 5.4246].*Z_base
# Convert for SPS
R_602_sps = R_602/mi2Km
L_602_sps = X_602/mi2Km/2/pi/60
C_602_sps = B_602/mi2Km*ms2F
# Configuration 603
R_603 = [0 0 0 ; 0 1.3294 0.2066; 0 0.2066 1.3238]./Z_base
X_603 = [0 0 0; 0 1.3471 0.4591; 0 0.4591 1.3569] ./Z_base
B_603 = [0 0 0; 0 4.7097 -0.8999;0 -0.8999 4.6658].*Z_base
# Convert for SPS
R_603_sps = R_603/mi2Km
L_603_sps = X_603/mi2Km/2/pi/60
C_603_sps = B_603/mi2Km*ms2F
# Configuration 604
R_604 = [1.3238 0 0.2066;0 0 0; 0.2066 0 1.3294]./Z_base
X_604 = [1.3569 0 0.4591;0 0 0; 0.4591 0 1.3471]./Z_base
B_604 = [4.6658 0 -0.8999;0 0 0; -0.8999 0 4.7097].*Z_base
# Convert for SPS
R_604_sps = R_604/mi2Km
L_604_sps = X_604/mi2Km/2/pi/60
C_604_sps = B_604/mi2Km*ms2F
#Configuration 605
R_605 = [0 0 0;0 0 0;0 0 1.3292]./Z_base
X_605 = [0 0 0;0 0 0;0 0 1.3475]./Z_base
B_605 = [0 0 0;0 0 0;0 0 4.5193].*Z_base
# Convert for SPS
R_605_sps = R_605/mi2Km
L_605_sps = X_605/mi2Km/2/pi/60
C_605_sps = B_605/mi2Km*ms2F
# Configuration 606
R_606 = [0.7982 0.3192 0.2849; 0.3192 0.7891 0.3192; 0.2849 0.3192 0.7982]./Z_base
X_606 = [0.4463 0.0328 0.0143; 0.0328 0.4041 0.0328; 0.0143 0.0328 0.4463]./Z_base
B_606 = [96.8897 -1e-6 -1e-6; -1e-6 96.8897 -1e-6; -1e-6 -1e-6 96.8897].*Z_base
# Convert for SPS
R_606_sps = R_606/mi2Km
L_606_sps = X_606/mi2Km/2/pi/60
C_606_sps = B_606/mi2Km*ms2F
# Configuration 607
R_607 = [1.3425 0 0;0 0 0;0 0 0]./Z_base
X_607 = [0.5124 0 0;0 0 0;0 0 0]./Z_base
B_607= [88.9912 0 0;0 0 0;0 0 0].*Z_base
# Convert for SPS
R_607_sps = R_607/mi2Km
L_607_sps = X_607/mi2Km/2/pi/60
C_607_sps = B_607/mi2Km*ms2F
# --- DĂFINITION DES LIGNES ---
Dict(
"branch" => Dict{String, Any}(
"1" => Dict{String, Any}( #632_633
"f_bus" => "632",
"t_bus" => "633",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_602_sps.*500*ft2km,
"br_x" => X_602.*500*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_602.*500*ft2km/2, "b_to" => B_602.*500*ft2km/2,
"br_status" => 1, "index" => 1,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"2" => Dict{String, Any}( #632_645
"f_bus" => "632",
"t_bus" => "645",
"f_connections" => [2, 3],
"t_connections" => [2, 3],
"br_r" => R_603_sps[2:3,2:3].*500*ft2km,
"br_x" => X_603[2:3,2:3].*500*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_603[2:3,2:3].*500*ft2km./2, "b_to" => B_603[2:3,2:3].*500*ft2km./2,
"br_status" => 1, "index" => 2,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"3" => Dict{String, Any}( #645_646
"f_bus" => "645",
"t_bus" => "646",
"f_connections" => [2, 3],
"t_connections" => [2, 3],
"br_r" => R_603_sps[2:3,2:3].*300*ft2km,
"br_x" => X_603[2:3,2:3].*300*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_603[2:3,2:3].*300*ft2km./2, "b_to" => B_603[2:3,2:3].*300*ft2km./2,
"br_status" => 1, "index" => 3,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"4" => Dict{String, Any}( #632_671
"f_bus" => "632",
"t_bus" => "671",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_601_sps.*2000*ft2km,
"br_x" => X_601.*2000*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_601.*2000*ft2km/2, "b_to" => B_601.*2000*ft2km/2,
"br_status" => 1, "index" => 4,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"5" => Dict{String, Any}( #692_675
"f_bus" => "692",
"t_bus" => "675",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_606_sps.*500*ft2km,
"br_x" => X_606.*500*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_606.*500*ft2km/2, "b_to" => B_606.*500*ft2km/2,
"br_status" => 1, "index" => 5,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"6" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "680",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_601_sps.*1000*ft2km,
"br_x" => X_601.*1000*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_601.*1000*ft2km/2, "b_to" => B_601.*1000*ft2km/2,
"br_status" => 1, "index" => 6,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"7" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "684",
"f_connections" => [1, 3],
"t_connections" => [1, 3],
"br_r" => R_604_sps[[1, 3], [1, 3]].*300*ft2km,
"br_x" => X_604[[1, 3], [1, 3]].*300*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_604[[1, 3], [1, 3]].*300*ft2km ./ 2, "b_to" => B_604[[1, 3], [1, 3]].*300*ft2km./ 2,
"br_status" => 1, "index" => 7,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"8" => Dict{String, Any}(
"f_bus" => "684",
"t_bus" => "652",
"f_connection" => [1],
"t_connection" => [1],
"br_r" => R_604_sps[1: 1, 1: 1].*800*ft2km,
"br_x" => X_604[1: 1, 1: 1].*800*ft2km,
"g_fr" => zeros(1,1), "g_to" => zeros(1,1),
"b_fr" => B_604[1: 1, 1: 1].*800*ft2km ./ 2, "b_to" => B_604[1: 1, 1: 1].*800*ft2km ./ 2,
"br_status" => 1, "index" => 8,
"angmin" => [-3.14], "angmax" => [3.14]
),
"9" => Dict{String, Any}(
"f_bus" => "684",
"t_bus" => "611",
"f_connection" => [3],
"t_connection" => [3],
"br_r" => R_605_sps[3: 3, 3: 3].*300*ft2km,
"br_x" => X_605[3: 3, 3: 3].*300*ft2km,
"g_fr" => zeros(1,1), "g_to" => zeros(1,1),
"b_fr" => B_605[3: 3, 3: 3].*300*ft2km ./ 2, "b_to" => B_605[3: 3, 3: 3].*300*ft2km ./ 2,
"br_status" => 1, "index" => 9,
"angmin" => [-3.14], "angmax" => [3.14]
),
)
)
#Erstellung der Lasten / Load
# 15 Lasten
# Dien Leistungen sind in kWatt und durch sbase = 1000 haben wir dies Per Unit
Dict(
"load"=> Dict{String,Any}(
"1" =>Dict{String, Any}( #632_Yq_PQ
"load_bus"=> "632",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [17, 66, 117]./2000.0,
"qd" => [10, 38, 68]./2000.0,
"status" => 1,
"index" => 1
),
"2" =>Dict{String, Any}( #634_Yg_PQ
"load_bus"=>"634",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [160, 120, 120]./1000.0,
"qd" => [110, 90, 90]./1000.0,
"status" => 1,
"index" => 2
),
"3" =>Dict{String, Any}( #645_PQ
"load_bus"=> "645",
"connections" =>[2],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
#"pd" => [0.0, 170, 0.0]./1000.0,
#"qd" => [0.0, 125, 0.0]./1000.0,
"pd" => [170]/1000.0,
"qd" => [125]/1000.0,
"status" => 1,
"index" => 3
),
"4" =>Dict{String, Any}( # a verifier completement #646_Z
"load_bus"=> "646",
"connections" =>[2,3],
"configuration" => PMD.DELTA, #DELTA weil sie zwichen 2 phasen ist
"model" => PMD.IMPEDANCE,
#"pd" => [0.0, 230, 0.0]./1000.0, #1 und 2 dann 2 und 3 danach 3 und 1 wiwr haben nur 2 und 3 andere sind 0
#"qd" => [0.0, 132, 0.0]./1000.0,
"pd" => [230, 0.0]./1000.0, #1 und 2 dann 2 und 3 danach 3 und 1 wiwr haben nur 2 und 3 andere sind 0
"qd" => [ 132, 0.0]./1000.0,
"status" => 1,
"index" => 4
),
"5" =>Dict{String, Any}( #671_Yq_PQ
"load_bus"=> "671",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [17, 66, 117]./2000.0,
"qd" => [10, 38, 68]./2000.0,
"status" => 1,
"index" => 5
),
"6" =>Dict{String, Any}( #671_D_PQ
"load_bus"=> "671",
"connections" =>[1,2,3],
"configuration" => PMD.DELTA,
"model" => PMD.POWER,
"pd" => [385, 385, 385 ]./1000.0,
"qd" => [220, 220, 220]./1000.0,
"status" => 1,
"index" => 6
),
"7" =>Dict{String, Any}( #675_Yg_Z
"load_bus"=> "675",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.IMPEDANCE,
"pd" => [0.0, 0.0, 0.0],
"qd" => [-200, -200, -200]./1000.0,
"status" => 1,
"index" => 7
),
"8" =>Dict{String, Any}( #675_Yg_PQ
"load_bus"=> "675",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [485, 68, 290]./1000.0, #en Watt
"qd" => [190, 60, 212]./1000.0,
"status" => 1,
"index" => 8
),
"9" =>Dict{String, Any}( #692_D_I
"load_bus"=> "692",
"connections" =>[1,2,3],
"configuration" => PMD.DELTA,
"model" => PMD.CURRENT,
"pd" => [1e-6, 1e-3, 170]./1000.0, # en watt
"qd" => [0.0, 0.0, 151]./1000.0, # en Var
"status" => 1,
"index" => 9
),
"10" =>Dict{String, Any}( #652_Z
"load_bus"=> "652",
"connection" =>[1],
"configuration" => PMD.WYE,
"model" => PMD.IMPEDANCE,
"pd" => [128]/1000.0,
"qd" => [86]/1000.0,
"status" => 1,
"index" => 10
),
"11" =>Dict{String, Any}( #611_Z
"load_bus"=> "611",
"connection" =>[3],
"configuration" => PMD.WYE,
"model" => PMD.IMPEDANCE,
"pd" => [0.0],
"qd" => [-100]/1000.0,
"status" => 1,
"index" => 11
),
"12" =>Dict{String, Any}( #611_I
"load_bus"=> "611",
"connection" =>[3],
"configuration" => PMD.WYE,
"model" => PMD.CURRENT,
"pd" => [170]/1000.0,
"qd" => [80]/1000.0,
"status" => 1,
"index" => 12
),
)
)
#Erstellung der Knoten / Bussen
# 13 Bussen
Dict(
"bus" => Dict{String, Any}(
"632" => Dict{String,Any}(
"bus_type" => 3,
"vbase" => 4.160/sqrt(3), #kV
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [true, true,true],
"index" => 1,
"vmin" => [0.1, 0.1, 0.1],
"vmax" => [1.1, 1.1, 1.1],
"vm" => [1.0210, 1.0420, 1.0174],
"va" => [-2.49*(pi/180), -121.72*(pi/180), 117.83*(pi/180)]
),
"633" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 2,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"634" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 0.480/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 3,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"645" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [2,3],
#"terminals" => [1,2,3],
"grounded" => [false,false],
"index" => 4,
"vmin" => [0.0, 0.0],
"vmax" => [1.5, 1.5]
#"vmin" => [0.0, 0.0, 0.0],
#"vmax" => [1.5, 1.5, 1.5]
),
"646" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [2,3],
"grounded" => [false,false],
"index" => 5,
#"vmin" => [0.0, 0.0, 0.0],
#"vmax" => [1.5, 1.5, 1.5]
"vmin" => [0.0, 0.0],
"vmax" => [1.5, 1.5]
),
"671" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 6,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"675" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 7,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"692" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 8,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"680" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 9,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"684" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,3],
"grounded" => [false,false],
"index" => 10,
"vmin" => [0.0, 0.0],
"vmax" => [1.5, 1.5]
),
"652" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1],
"grounded" => [false],
"index" => 11,
"vmin" => [0.0],
"vmax" => [1.5]
),
"611" => Dict{String,Any}(#Erstellung der Knoten / Bussen
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [3],
"grounded" => [false],
"index" => 12,
"vmin" => [0.0],
"vmax" => [1.5]
)
)
)
Hi,
A 13-bus distribution network is a small case. I think itâs a good starting point to use a toy example like that model, if you are relatively inexperienced. Since itâs small, you can directly draw that network as a graph on a paper so you can check whether it is a proper network. And after solution you can analyze the power flowsâŠ
From my personal experience about numeric computation (especially optimization), the data can generally be generated via random number in proper ranges. Therefore all datasets are not mandatory for me. I would borrow only the graph structure (i.e. the branch-node connection pattern) from the standard cases. And if the existing graph structure is not standard (e.g. there are isolated nodes, or there are redundant parallel branches), I would manually modify it so it can serve as a proper numeric test case that I can write algorithms on.
1 Like
Thank you very much.
I did that, but the problem persists.
Hi @Curtis96, welcome to the forum 
Iâve edited your previous posts to format them correctly.
You needed triple back-ticks, like this:
```julia
1 + 1
```
Unfortunately I donât know if I have the domain expertise the understand your problem. Youâve probably chosen the wrong combination of solve_mc_pf and the ACPUPowerModel.
@ccoffrin may be able to point us in the direction of someone who can help.
Hereâs a reproducible example for future reference
using PowerModelsDistribution
using Ipopt
using Printf
const PMD = PowerModelsDistribution
using LinearAlgebra
# Constantes de conversion
mi2Km = 1.60934
ft2km = 0.0003048
ms2F = 1/2/pi/60*1e-6
# Per unit umrechnen: Z_base = (V_base)^2/S_base
V_base = 4.16
Z_base = (V_base)^2 /1.0
# Configuration 601 (Triphasée - Standard)
R_601 = [0.3465 0.1560 0.1580; 0.1560 0.3375 0.1535; 0.1580 0.1535 0.3414]./Z_base
X_601 = [1.0179 0.5017 0.4236; 0.5017 1.0478 0.3849; 0.4236 0.3849 1.0348] ./Z_base
B_601 = [6.2998 -1.9958 -1.2595; -1.9958 5.9597 -0.7417; -1.2595 -0.7417 5.6386].*Z_base
# Convert for SPS
R_601_sps = R_601/mi2Km
L_601_sps = X_601/mi2Km/2/pi/60
C_601_sps = B_601/mi2Km*ms2F
# Configuration 602 (Triphasée)
R_602 = [0.7526 0.1580 0.1560; 0.1580 0.7475 0.1535; 0.1560 0.1535 0.7436]./Z_base
X_602 = [1.1814 0.4236 0.5017; 0.4236 1.1983 0.3849; 0.5017 0.3849 1.2112]./Z_base
B_602 = [5.6990 -1.0817 -1.6905; -1.0817 5.1795 -0.6588; -1.6905 -0.6588 5.4246].*Z_base
# Convert for SPS
R_602_sps = R_602/mi2Km
L_602_sps = X_602/mi2Km/2/pi/60
C_602_sps = B_602/mi2Km*ms2F
# Configuration 603
R_603 = [0 0 0 ; 0 1.3294 0.2066; 0 0.2066 1.3238]./Z_base
X_603 = [0 0 0; 0 1.3471 0.4591; 0 0.4591 1.3569] ./Z_base
B_603 = [0 0 0; 0 4.7097 -0.8999;0 -0.8999 4.6658].*Z_base
# Convert for SPS
R_603_sps = R_603/mi2Km
L_603_sps = X_603/mi2Km/2/pi/60
C_603_sps = B_603/mi2Km*ms2F
# Configuration 604
R_604 = [1.3238 0 0.2066;0 0 0; 0.2066 0 1.3294]./Z_base
X_604 = [1.3569 0 0.4591;0 0 0; 0.4591 0 1.3471]./Z_base
B_604 = [4.6658 0 -0.8999;0 0 0; -0.8999 0 4.7097].*Z_base
# Convert for SPS
R_604_sps = R_604/mi2Km
L_604_sps = X_604/mi2Km/2/pi/60
C_604_sps = B_604/mi2Km*ms2F
#Configuration 605
R_605 = [0 0 0;0 0 0;0 0 1.3292]./Z_base
X_605 = [0 0 0;0 0 0;0 0 1.3475]./Z_base
B_605 = [0 0 0;0 0 0;0 0 4.5193].*Z_base
# Convert for SPS
R_605_sps = R_605/mi2Km
L_605_sps = X_605/mi2Km/2/pi/60
C_605_sps = B_605/mi2Km*ms2F
# Configuration 606
R_606 = [0.7982 0.3192 0.2849; 0.3192 0.7891 0.3192; 0.2849 0.3192 0.7982]./Z_base
X_606 = [0.4463 0.0328 0.0143; 0.0328 0.4041 0.0328; 0.0143 0.0328 0.4463]./Z_base
B_606 = [96.8897 -1e-6 -1e-6; -1e-6 96.8897 -1e-6; -1e-6 -1e-6 96.8897].*Z_base
# Convert for SPS
R_606_sps = R_606/mi2Km
L_606_sps = X_606/mi2Km/2/pi/60
C_606_sps = B_606/mi2Km*ms2F
# Configuration 607
R_607 = [1.3425 0 0;0 0 0;0 0 0]./Z_base
X_607 = [0.5124 0 0;0 0 0;0 0 0]./Z_base
B_607= [88.9912 0 0;0 0 0;0 0 0].*Z_base
# Convert for SPS
R_607_sps = R_607/mi2Km
L_607_sps = X_607/mi2Km/2/pi/60
C_607_sps = B_607/mi2Km*ms2F
line_dict = Dict(
"branch" => Dict{String, Any}(
"1" => Dict{String, Any}( #632_633
"f_bus" => "632",
"t_bus" => "633",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_602_sps.*500*ft2km,
"br_x" => X_602.*500*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_602.*500*ft2km/2, "b_to" => B_602.*500*ft2km/2,
"br_status" => 1, "index" => 1,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"2" => Dict{String, Any}( #632_645
"f_bus" => "632",
"t_bus" => "645",
"f_connections" => [2, 3],
"t_connections" => [2, 3],
"br_r" => R_603_sps[2:3,2:3].*500*ft2km,
"br_x" => X_603[2:3,2:3].*500*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_603[2:3,2:3].*500*ft2km./2, "b_to" => B_603[2:3,2:3].*500*ft2km./2,
"br_status" => 1, "index" => 2,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"3" => Dict{String, Any}( #645_646
"f_bus" => "645",
"t_bus" => "646",
"f_connections" => [2, 3],
"t_connections" => [2, 3],
"br_r" => R_603_sps[2:3,2:3].*300*ft2km,
"br_x" => X_603[2:3,2:3].*300*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_603[2:3,2:3].*300*ft2km./2, "b_to" => B_603[2:3,2:3].*300*ft2km./2,
"br_status" => 1, "index" => 3,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"4" => Dict{String, Any}( #632_671
"f_bus" => "632",
"t_bus" => "671",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_601_sps.*2000*ft2km,
"br_x" => X_601.*2000*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_601.*2000*ft2km/2, "b_to" => B_601.*2000*ft2km/2,
"br_status" => 1, "index" => 4,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"5" => Dict{String, Any}( #692_675
"f_bus" => "692",
"t_bus" => "675",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_606_sps.*500*ft2km,
"br_x" => X_606.*500*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_606.*500*ft2km/2, "b_to" => B_606.*500*ft2km/2,
"br_status" => 1, "index" => 5,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"6" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "680",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"br_r" => R_601_sps.*1000*ft2km,
"br_x" => X_601.*1000*ft2km,
"g_fr" => zeros(3,3), "g_to" => zeros(3,3),
"b_fr" => B_601.*1000*ft2km/2, "b_to" => B_601.*1000*ft2km/2,
"br_status" => 1, "index" => 6,
"angmin" => [-3.14, -3.14, -3.14], "angmax" => [3.14, 3.14, 3.14]
),
"7" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "684",
"f_connections" => [1, 3],
"t_connections" => [1, 3],
"br_r" => R_604_sps[[1, 3], [1, 3]].*300*ft2km,
"br_x" => X_604[[1, 3], [1, 3]].*300*ft2km,
"g_fr" => zeros(2,2), "g_to" => zeros(2,2),
"b_fr" => B_604[[1, 3], [1, 3]].*300*ft2km ./ 2, "b_to" => B_604[[1, 3], [1, 3]].*300*ft2km./ 2,
"br_status" => 1, "index" => 7,
"angmin" => [-3.14, -3.14], "angmax" => [3.14, 3.14]
),
"8" => Dict{String, Any}(
"f_bus" => "684",
"t_bus" => "652",
"f_connection" => [1],
"t_connection" => [1],
"br_r" => R_604_sps[1: 1, 1: 1].*800*ft2km,
"br_x" => X_604[1: 1, 1: 1].*800*ft2km,
"g_fr" => zeros(1,1), "g_to" => zeros(1,1),
"b_fr" => B_604[1: 1, 1: 1].*800*ft2km ./ 2, "b_to" => B_604[1: 1, 1: 1].*800*ft2km ./ 2,
"br_status" => 1, "index" => 8,
"angmin" => [-3.14], "angmax" => [3.14]
),
"9" => Dict{String, Any}(
"f_bus" => "684",
"t_bus" => "611",
"f_connection" => [3],
"t_connection" => [3],
"br_r" => R_605_sps[3: 3, 3: 3].*300*ft2km,
"br_x" => X_605[3: 3, 3: 3].*300*ft2km,
"g_fr" => zeros(1,1), "g_to" => zeros(1,1),
"b_fr" => B_605[3: 3, 3: 3].*300*ft2km ./ 2, "b_to" => B_605[3: 3, 3: 3].*300*ft2km ./ 2,
"br_status" => 1, "index" => 9,
"angmin" => [-3.14], "angmax" => [3.14]
),
)
)
#Erstellung der Lasten / Load
# 15 Lasten
# Dien Leistungen sind in kWatt und durch sbase = 1000 haben wir dies Per Unit
load_dict = Dict(
"load"=> Dict{String,Any}(
"1" =>Dict{String, Any}( #632_Yq_PQ
"load_bus"=> "632",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [17, 66, 117]./2000.0,
"qd" => [10, 38, 68]./2000.0,
"status" => 1,
"index" => 1
),
"2" =>Dict{String, Any}( #634_Yg_PQ
"load_bus"=>"634",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [160, 120, 120]./1000.0,
"qd" => [110, 90, 90]./1000.0,
"status" => 1,
"index" => 2
),
"3" =>Dict{String, Any}( #645_PQ
"load_bus"=> "645",
"connections" =>[2],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
#"pd" => [0.0, 170, 0.0]./1000.0,
#"qd" => [0.0, 125, 0.0]./1000.0,
"pd" => [170]/1000.0,
"qd" => [125]/1000.0,
"status" => 1,
"index" => 3
),
"4" =>Dict{String, Any}( # a verifier completement #646_Z
"load_bus"=> "646",
"connections" =>[2,3],
"configuration" => PMD.DELTA, #DELTA weil sie zwichen 2 phasen ist
"model" => PMD.IMPEDANCE,
#"pd" => [0.0, 230, 0.0]./1000.0, #1 und 2 dann 2 und 3 danach 3 und 1 wiwr haben nur 2 und 3 andere sind 0
#"qd" => [0.0, 132, 0.0]./1000.0,
"pd" => [230, 0.0]./1000.0, #1 und 2 dann 2 und 3 danach 3 und 1 wiwr haben nur 2 und 3 andere sind 0
"qd" => [ 132, 0.0]./1000.0,
"status" => 1,
"index" => 4
),
"5" =>Dict{String, Any}( #671_Yq_PQ
"load_bus"=> "671",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [17, 66, 117]./2000.0,
"qd" => [10, 38, 68]./2000.0,
"status" => 1,
"index" => 5
),
"6" =>Dict{String, Any}( #671_D_PQ
"load_bus"=> "671",
"connections" =>[1,2,3],
"configuration" => PMD.DELTA,
"model" => PMD.POWER,
"pd" => [385, 385, 385 ]./1000.0,
"qd" => [220, 220, 220]./1000.0,
"status" => 1,
"index" => 6
),
"7" =>Dict{String, Any}( #675_Yg_Z
"load_bus"=> "675",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.IMPEDANCE,
"pd" => [0.0, 0.0, 0.0],
"qd" => [-200, -200, -200]./1000.0,
"status" => 1,
"index" => 7
),
"8" =>Dict{String, Any}( #675_Yg_PQ
"load_bus"=> "675",
"connections" =>[1,2,3],
"configuration" => PMD.WYE,
"model" => PMD.POWER,
"pd" => [485, 68, 290]./1000.0, #en Watt
"qd" => [190, 60, 212]./1000.0,
"status" => 1,
"index" => 8
),
"9" =>Dict{String, Any}( #692_D_I
"load_bus"=> "692",
"connections" =>[1,2,3],
"configuration" => PMD.DELTA,
"model" => PMD.CURRENT,
"pd" => [1e-6, 1e-3, 170]./1000.0, # en watt
"qd" => [0.0, 0.0, 151]./1000.0, # en Var
"status" => 1,
"index" => 9
),
"10" =>Dict{String, Any}( #652_Z
"load_bus"=> "652",
"connection" =>[1],
"configuration" => PMD.WYE,
"model" => PMD.IMPEDANCE,
"pd" => [128]/1000.0,
"qd" => [86]/1000.0,
"status" => 1,
"index" => 10
),
"11" =>Dict{String, Any}( #611_Z
"load_bus"=> "611",
"connection" =>[3],
"configuration" => PMD.WYE,
"model" => PMD.IMPEDANCE,
"pd" => [0.0],
"qd" => [-100]/1000.0,
"status" => 1,
"index" => 11
),
"12" =>Dict{String, Any}( #611_I
"load_bus"=> "611",
"connection" =>[3],
"configuration" => PMD.WYE,
"model" => PMD.CURRENT,
"pd" => [170]/1000.0,
"qd" => [80]/1000.0,
"status" => 1,
"index" => 12
),
)
)
#Erstellung der Knoten / Bussen
# 13 Bussen
bus_dict = Dict(
"bus" => Dict{String, Any}(
"632" => Dict{String,Any}(
"bus_type" => 3,
"vbase" => 4.160/sqrt(3), #kV
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [true, true,true],
"index" => 1,
"vmin" => [0.1, 0.1, 0.1],
"vmax" => [1.1, 1.1, 1.1],
"vm" => [1.0210, 1.0420, 1.0174],
"va" => [-2.49*(pi/180), -121.72*(pi/180), 117.83*(pi/180)]
),
"633" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 2,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"634" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 0.480/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 3,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"645" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [2,3],
#"terminals" => [1,2,3],
"grounded" => [false,false],
"index" => 4,
"vmin" => [0.0, 0.0],
"vmax" => [1.5, 1.5]
#"vmin" => [0.0, 0.0, 0.0],
#"vmax" => [1.5, 1.5, 1.5]
),
"646" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [2,3],
"grounded" => [false,false],
"index" => 5,
#"vmin" => [0.0, 0.0, 0.0],
#"vmax" => [1.5, 1.5, 1.5]
"vmin" => [0.0, 0.0],
"vmax" => [1.5, 1.5]
),
"671" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 6,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"675" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 7,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"692" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 8,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"680" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,2,3],
"grounded" => [false, false,false],
"index" => 9,
"vmin" => [0.0, 0.0, 0.0],
"vmax" => [1.5, 1.5, 1.5]
),
"684" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1,3],
"grounded" => [false,false],
"index" => 10,
"vmin" => [0.0, 0.0],
"vmax" => [1.5, 1.5]
),
"652" => Dict{String,Any}(
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [1],
"grounded" => [false],
"index" => 11,
"vmin" => [0.0],
"vmax" => [1.5]
),
"611" => Dict{String,Any}(#Erstellung der Knoten / Bussen
"bus_type" => 1,
"vbase" => 4.160/sqrt(3),
"vnom_kv" => 4.16 / sqrt(3),
"terminals" => [3],
"grounded" => [false],
"index" => 12,
"vmin" => [0.0],
"vmax" => [1.5]
)
)
)
math_data = Dict{String, Any}(
"bus" => bus_dict["bus"],
"load" => load_dict["load"],
"branch" => line_dict["branch"],
"source" => Dict{String, Any}(
"1" => Dict{String, Any}("connections" => [1, 2, 3],
"bus" => "632",
"index" => 1,
"status" => 1,
"vm" => [1.0210, 1.0420, 1.0174],
"va" => [-2.49*(pi/180), -121.72*(pi/180), 117.83*(pi/180)]
)
),
"gen" => Dict{String, Any}(
"1" => Dict{String, Any}(
"connections" => [1, 2, 3],
"gen_bus" => "632", "index" => 1,
"gen_status" => 1,
"gen_type" => 3,
#"vm" => [ 1.0210, 1.0420, 1.0174],
#"va" => [ -2.49*(pi /180), -121.72*(pi /180), 117.83*(pi /180)], #radiant winckel
"status" =>1,
#"index" =>1,
"pg" => [5.0, 5.0, 5.0],
"qg" => [0.1, 0.1, 0.1],
"pmax" => [10.0, 10.0, 10.0],
"pmin" => [0.0, 0.0, 0.0],
"qmax" => [10.0, 10.0, 10.0],
"qmin" => [-10.0, -10.0, -10.0]
)
),
"switch" => Dict{String, Any}(
"1" => Dict{String, Any}(
"f_bus" => "671",
"t_bus" => "692",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"state" => 1,
"status" => 1,
"index" => 1,
"rs" => [0.001, 0.001, 0.001],
"xs" => [0.0, 0.0, 0.0],
"angmin" => [-1.0, -1.0, -1.0], "argmax"=>[1.0,1.0,1.0]
)
),
"transformer" => Dict{String, Any}(
"1" => Dict{String, Any}(
"f_bus" => "633",
"t_bus" => "634",
"f_connections" => [1, 2, 3],
"t_connections" => [1, 2, 3],
"configuration" => PMD.WYE,
"tm_nom" => 1.0,
"tm_set" => [1.0, 1.0, 1.0],
"polarity" => 1,
"vbase" => [4.16, 0.48],
"xsc" => [0.8, 0.8, 0.8],
"rw" => [0.1, 0.1, 0.1],
"index" => 1,
"status" => 1,
"tm_status" => 1
)
),
"settings" => Dict{String, Any}("sbase" => 1.0, "vbase" => 4.16),
"per_unit" => true, "data_model" => PMD.MATHEMATICAL
)
opt = optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6)
resultat_pf = PMD.solve_mc_pf(math_data, PMD.ACPUPowerModel, opt)