Library logo

Buildings.Electrical.AC.ThreePhasesUnbalanced.Validation.IEEETests.Test4NodesFeeder.BaseClasses

This package contains the base classes used by the IEEE 4 nodes test feeder

Information

This package contains base classes used by the models that are part of the package Buildings.Electrical.AC.ThreePhasesUnbalanced.Validation.IEEETests.Test4NodesFeeder.

Extends from Modelica.Icons.BasesPackage (Icon for packages containing base classes).

Package Content

Name Description
Buildings.Electrical.AC.ThreePhasesUnbalanced.Validation.IEEETests.Test4NodesFeeder.BaseClasses.IEEE4 IEEE4 Base model of the IEEE 4 nodes test feeder

Buildings.Electrical.AC.ThreePhasesUnbalanced.Validation.IEEETests.Test4NodesFeeder.BaseClasses.IEEE4 Buildings.Electrical.AC.ThreePhasesUnbalanced.Validation.IEEETests.Test4NodesFeeder.BaseClasses.IEEE4

Base model of the IEEE 4 nodes test feeder

Buildings.Electrical.AC.ThreePhasesUnbalanced.Validation.IEEETests.Test4NodesFeeder.BaseClasses.IEEE4

Information

This is a partial model that is extended by all the other validation test cases. This model defined replaceable probes and transformer so they can be easily changed when implementing the different tests.

image

More information can be found in each model that extends this one.

Reference results

The reference results for the tests are saved as parameter of the model and compared to the simulated ones. The error between the results (herein called Xmodel) and the references (herein called Xref) are computed in both absolute and relative way. Note that Xmodel and Xref can be either voltage amplitudes or phase angles.

Errabs = Xmodel - Xref

Err% = Errabs / Xref

The variables that store the results of the comparison are listed in the table below

Variable Description Unit
err_V2[3] Error between simulated voltage at node 2 and reference results [V]
err_V3[3] Error between simulated voltage at node 3 and reference results [V]
err_V4[3] Error between simulated voltage at node 4 and reference results [V]
err_Theta2[3] Error between simulated phase angle at node 2 and reference phase angle [rad], displayed as [deg]
err_Theta3[3] Error between simulated phase angle at node 2 and reference phase angle [rad], displayed as [deg]
err_Theta4[3] Error between simulated phase angle at node 2 and reference phase angle [rad], displayed as [deg]
err_V2_percent[3] Relative error between simulated voltage at node 2 and reference results [%]
err_V3_percent[3] Relative error between simulated voltage at node 3 and reference results [%]
err_V4_percent[3] Relative error between simulated voltage at node 4 and reference results [%]
err_Theta2_percent[3] Relative error between simulated phase angle at node 2 and reference phase angle [%]
err_Theta3_percent[3] Relative error between simulated phase angle at node 2 and reference phase angle [%]
err_Theta4_percent[3] Relative error between simulated phase angle at node 2 and reference phase angle [%]

Extends from Modelica.Icons.Example (Icon for runnable examples).

Parameters

TypeNameDefaultDescription
VoltageVLL_side112.47e3Voltage line to line side 1 [V]
VoltageVLL_side24.16e3Voltage line to line side 2 [V]
ApparentPowerVARbase6000e3Base VA power of the transformer [V.A]
Booleanline1_use_Z_ytrueChoose between Zy or Zd impedance matrix for line 1
Booleanline2_use_Z_ytrueChoose between Zy or Zd impedance matrix for line 2
VoltageV2_ref[3]{7107,7140,7121}Reference RMS voltage node 2 - IEEE results [V]
VoltageV3_ref[3]{2247,2269,2256}Reference RMS voltage node 3 - IEEE results [V]
VoltageV4_ref[3]{1918,2061,1981}Reference RMS voltage node 4 - IEEE results [V]
AngleTheta2_ref[3]{-0.3,-120.3,119.6}Reference voltage phase angle node 2 - IEEE results [rad]
AngleTheta3_ref[3]{-3.7,-123.5,116.4}Reference voltage phase angle node 3 - IEEE results [rad]
AngleTheta4_ref[3]{-9.1,-128.3,110.9}Reference voltage phase angle node 4 - IEEE results [rad]
GeneralizedProbenode1node1(perUnit=false, V_nomin...Probe at source
GeneralizedProbenode2node2(perUnit=false, V_nomin...Probe at the primary side of the transformer
GeneralizedProbenode3node3(perUnit=false, V_nomin...Probe at the secondary side of the transformer
GeneralizedProbenode4node4(perUnit=false, V_nomin...Probe at the load

Modelica definition

partial model IEEE4 "Base model of the IEEE 4 nodes test feeder" extends Modelica.Icons.Example; parameter Modelica.Units.SI.Voltage VLL_side1=12.47e3 "Voltage line to line side 1"; parameter Modelica.Units.SI.Voltage VLL_side2=4.16e3 "Voltage line to line side 2"; parameter Modelica.Units.SI.ApparentPower VARbase=6000e3 "Base VA power of the transformer"; parameter Boolean line1_use_Z_y = true "Choose between Zy or Zd impedance matrix for line 1"; parameter Boolean line2_use_Z_y = true "Choose between Zy or Zd impedance matrix for line 2"; parameter Modelica.Units.SI.Voltage V2_ref[3]={7107,7140,7121} "Reference RMS voltage node 2 - IEEE results"; parameter Modelica.Units.SI.Voltage V3_ref[3]={2247,2269,2256} "Reference RMS voltage node 3 - IEEE results"; parameter Modelica.Units.SI.Voltage V4_ref[3]={1918,2061,1981} "Reference RMS voltage node 4 - IEEE results"; parameter Modelica.Units.SI.Angle Theta2_ref[3](each displayUnit="deg") = {-0.3, -120.3,119.6} "Reference voltage phase angle node 2 - IEEE results"; parameter Modelica.Units.SI.Angle Theta3_ref[3](each displayUnit="deg") = {-3.7, -123.5,116.4} "Reference voltage phase angle node 3 - IEEE results"; parameter Modelica.Units.SI.Angle Theta4_ref[3](each displayUnit="deg") = {-9.1, -128.3,110.9} "Reference voltage phase angle node 4 - IEEE results"; Modelica.Units.SI.Voltage err_V2[3]=node2.V - V2_ref "Error on voltage at node 2"; Modelica.Units.SI.Voltage err_V3[3]=node3.V - V3_ref "Error on voltage at node 3"; Modelica.Units.SI.Voltage err_V4[3]=node4.V - V4_ref "Error on voltage at node 4"; Modelica.Units.SI.Angle err_Theta2[3](each displayUnit="deg") = node2.theta - Theta2_ref "Error on voltage at node 2"; Modelica.Units.SI.Angle err_Theta3[3](each displayUnit="deg") = node3.theta - Theta3_ref "Error on voltage at node 3"; Modelica.Units.SI.Angle err_Theta4[3](each displayUnit="deg") = node4.theta - Theta4_ref "Error on voltage at node 4"; Real err_V2_percent[3] = 100*{err_V2[i]/V2_ref[i] for i in 1:3} "Error in RMS voltage at node 2 -- percent"; Real err_V3_percent[3] = 100*{err_V3[i]/V3_ref[i] for i in 1:3} "Error in RMS voltage at node 3 -- percent"; Real err_V4_percent[3] = 100*{err_V4[i]/V4_ref[i] for i in 1:3} "Error in RMS voltage at node 4 -- percent"; Real err_Theta2_percent[3] = 100*{err_Theta2[i]/Theta2_ref[i] for i in 1:3} "Error in voltage phase angle at node 2 -- percent"; Real err_Theta3_percent[3] = 100*{err_Theta3[i]/Theta3_ref[i] for i in 1:3} "Error in voltage phase angle at node 3 -- percent"; Real err_Theta4_percent[3] = 100*{err_Theta4[i]/Theta4_ref[i] for i in 1:3} "Error in voltage phase angle at node 4 -- percent"; Buildings.Electrical.AC.ThreePhasesUnbalanced.Sources.FixedVoltage source( f=60, V=VLL_side1) "Voltage source"; Buildings.Electrical.AC.ThreePhasesUnbalanced.Lines.TwoPortMatrixRL line1( Z11=L1*(if line1_use_Z_y then Z11_y else Z11_d), Z12=L1*(if line1_use_Z_y then Z12_y else Z12_d), Z13=L1*(if line1_use_Z_y then Z13_y else Z13_d), Z22=L1*(if line1_use_Z_y then Z22_y else Z22_d), Z23=L1*(if line1_use_Z_y then Z23_y else Z23_d), Z33=L1*(if line1_use_Z_y then Z33_y else Z33_d), V_nominal=VLL_side1) "Line at primary side"; Buildings.Electrical.AC.ThreePhasesUnbalanced.Lines.TwoPortMatrixRL line2( Z11=L2*(if line2_use_Z_y then Z11_y else Z11_d), Z12=L2*(if line2_use_Z_y then Z12_y else Z12_d), Z13=L2*(if line2_use_Z_y then Z13_y else Z13_d), Z22=L2*(if line2_use_Z_y then Z22_y else Z22_d), Z23=L2*(if line2_use_Z_y then Z23_y else Z23_d), Z33=L2*(if line2_use_Z_y then Z33_y else Z33_d), V_nominal=VLL_side2) "Line at secondary side"; Buildings.Electrical.AC.ThreePhasesUnbalanced.Loads.Inductive loadRL( pf=0.9, V_nominal=VLL_side2, mode=Buildings.Electrical.Types.Load.VariableZ_P_input, use_pf_in=true) "Load"; replaceable Buildings.Electrical.AC.ThreePhasesUnbalanced.Sensors.BaseClasses.GeneralizedProbe node1(perUnit=false, V_nominal=VLL_side1) "Probe at source"; replaceable Buildings.Electrical.AC.ThreePhasesUnbalanced.Sensors.BaseClasses.GeneralizedProbe node2(perUnit=false, V_nominal=VLL_side1) "Probe at the primary side of the transformer"; replaceable Buildings.Electrical.AC.ThreePhasesUnbalanced.Sensors.BaseClasses.GeneralizedProbe node3(perUnit=false, V_nominal=VLL_side2) "Probe at the secondary side of the transformer"; replaceable Buildings.Electrical.AC.ThreePhasesUnbalanced.Sensors.BaseClasses.GeneralizedProbe node4(perUnit=false, V_nominal=VLL_side2) "Probe at the load"; protected parameter Real L1 = 2000*(1.0/5280.0) "Length line 1 in miles"; parameter Real L2 = 2500*(1.0/5280.0) "Length line 2 in miles"; parameter Modelica.Units.SI.Impedance Z11_d[2]={0.4013,1.4133} "Element [1,1] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z12_d[2]={0.0953,0.8515} "Element [1,2] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z13_d[2]={0.0953,0.7266} "Element [1,3] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z22_d[2]={0.4013,1.4133} "Element [2,2] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z23_d[2]={0.0953,0.7802} "Element [2,3] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z33_d[2]={0.4013,1.4133} "Element [3,3] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z11_y[2]={0.4576,1.078} "Element [1,1] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z12_y[2]={0.1559,0.5017} "Element [1,2] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z13_y[2]={0.1535,0.3849} "Element [1,3] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z22_y[2]={0.4666,1.0482} "Element [2,2] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z23_y[2]={0.158,0.4236} "Element [2,3] of impedance matrix"; parameter Modelica.Units.SI.Impedance Z33_y[2]={0.4615,1.0651} "Element [3,3] of impedance matrix"; equation connect(source.terminal, line1.terminal_n); connect(line2.terminal_p, loadRL.terminal); end IEEE4;