example network changed

change to a robust example network

docs/src/networks.md

@@ -6,28 +6,25 @@ The `Net` module provides functionality for creating and manipulating power syst
 # Import the Net module
   using Sparlectra
   using BenchmarkTools
-  # Create a network
-  # Slack bus B1 at 220 kV with 1.0 pu voltage and 0.0 pu angle
-  # PQ bus B2 at 220 kV with 1.0 pu voltage and 0.0 pu angle
-  # PQ bus B3 at 22 kV with 1.0 pu voltage and 0.0 pu angle
-  # Shunt at bus B3 with 0.0 kW and 150.0 kVar
-  # AC line from bus B1 to B2 at 220 kV with 100 km length, 0.0653 ohm/km resistance, 0.398 ohm/km reactance, 9.08 nF/km capacitance, and 0.0 power factor
-  # 2W transformer from bus B2 to B3 with 1000 MVA rating, 13.0% voltage ratio, 0.28% resistance, 20.0 kW power factor, and 0.06% no-load current
-  # Energy consumer at bus B3 with 285.0 kW and 200.0 kVar
-  # Synchronous machine at bus B1 with 1.02 pu voltage and 0.0 pu angle
-  #
-  # G->1* ----L---- 2 --T-- 3->L
-  #                 x
-  #   
-  net = Net(name = "testnet", baseMVA = 100.0)
-  addBus!(net = net, busName = "B1", busType = "Slack", vn_kV = 220.0, vm_pu = 1.0, va_deg = 0.0)
-  addBus!(net = net, busName = "B2", busType = "PQ", vn_kV = 220.0, vm_pu = 1.0, va_deg = 0.0)
-  addBus!(net = net, busName = "B3", busType = "PQ", vn_kV = 22.0, vm_pu = 1.0, va_deg = 0.0)
-  addShunt!(net = net, busName = "B3", pShunt = 0.0, qShunt = 150.0)
-  addACLine!(net = net, fromBus = "B1", toBus = "B2", length = 50.0, r = 0.0653, x = 0.398, c_nf_per_km = 9.08, tanδ = 0.0)
-  add2WTrafo!(net = net, fromBus = "B2", toBus = "B3", sn_mva = 1000.0, vk_percent = 13.0, vkr_percent = 0.28, pfe_kw = 20.0, i0_percent = 0.06)
-  addProsumer!(net = net, busName = "B3", type = "ENERGYCONSUMER", p = 50.0, q = 200.0)
-  addProsumer!(net = net, busName = "B1", type = "SYNCHRONOUSMACHINE", referencePri = "B1", vm_pu = 1.02, va_deg = 0.0)
+  net = Net(name = "case5", baseMVA = 100.0)
+  addBus!(net = net, busName = "B1", busType = "PQ", vn_kV = 110.0, vm_pu = 1.0, va_deg = 0.0)
+  addBus!(net = net, busName = "B2", busType = "PQ", vn_kV = 110.0, vm_pu = 1.0, va_deg = 0.0)
+  addBus!(net = net, busName = "B3", busType = "PQ", vn_kV = 110.0, vm_pu = 1.0, va_deg = 0.0)
+  addBus!(net = net, busName = "B4", busType = "PQ", vn_kV = 110.0, vm_pu = 1.0, va_deg = 0.0)
+  addBus!(net = net, busName = "B5", busType = "Slack", vn_kV = 110.0, vm_pu = 1.0, va_deg = 0.0)
+
+  addACLine!(net = net, fromBus = "B1", toBus = "B2", length = 25.0, r = 0.2, x = 0.39)
+  addACLine!(net = net, fromBus = "B1", toBus = "B3", length = 25.0, r = 0.2, x = 0.39)
+  addACLine!(net = net, fromBus = "B2", toBus = "B4", length = 25.0, r = 0.2, x = 0.39)
+  addACLine!(net = net, fromBus = "B3", toBus = "B4", length = 25.0, r = 0.2, x = 0.39)
+  addACLine!(net = net, fromBus = "B4", toBus = "B5", length = 25.0, r = 0.2, x = 0.39)
+
+  addProsumer!(net = net, busName = "B1", type = "ENERGYCONSUMER", p = 1.0, q = 2.0)
+  addProsumer!(net = net, busName = "B2", type = "ENERGYCONSUMER", p = 1.0, q = 2.0)
+  addProsumer!(net = net, busName = "B3", type = "ENERGYCONSUMER", p = 1.0, q = 2.0)
+
+  addProsumer!(net = net, busName = "B5", type = "SYNCHRONMASCHINE", referencePri = "B5", vm_pu = 1.0, va_deg = 0.0)
+  addProsumer!(net = net, busName = "B1", type = "GENERATOR", p = 1.1, q = 2.0)
 
   # Run power flow
   tol = 1e-6