Updates to the IEEE Microgrid Example

fachif · fachif · commit 42e306bbad78 · 2024-01-15T23:50:20.000-05:00
diff --git a/OpenIPSL/Data/PowerPlant/GenerationGroup/Generator1.mo b/OpenIPSL/Data/PowerPlant/GenerationGroup/Generator1.mo
@@ -130,8 +130,7 @@ equation
           fillPattern=FillPattern.Solid,
           textString="%name"),
                          Line(points={{-60,-20},{-20,20},{20,-20},{60,20}},
-          color={28,108,200}),Ellipse(extent={{-100,-100},{100,100}}, lineColor
-            =
+          color={28,108,200}),Ellipse(extent={{-100,-100},{100,100}}, lineColor=
            {28,108,200})}),
     Documentation(info="<html>
 <p>Generation group for the example that uses the IEEE421.5-based data set.</p>
diff --git a/OpenIPSL/Data/PowerPlant/GenerationGroup/Generator2.mo b/OpenIPSL/Data/PowerPlant/GenerationGroup/Generator2.mo
@@ -127,8 +127,7 @@ equation
           fillPattern=FillPattern.Solid,
           textString="%name"),
                          Line(points={{-60,-20},{-20,20},{20,-20},{60,20}},
-          color={28,108,200}),Ellipse(extent={{-100,-100},{100,100}}, lineColor
-            =
+          color={28,108,200}),Ellipse(extent={{-100,-100},{100,100}}, lineColor=
            {28,108,200})}),
     Documentation(info="<html>
 <p>Generation Group for the example that uses the Anderson-based data set.</p>
diff --git a/OpenIPSL/Electrical/Controls/VariableSpeedDrive/ControllerLogic/VoltsHertzController.mo b/OpenIPSL/Electrical/Controls/VariableSpeedDrive/ControllerLogic/VoltsHertzController.mo
@@ -24,7 +24,7 @@ model VoltsHertzController "Volts/Hertz controller model"
         rotation=180,
         origin={120,40})));
   OpenIPSL.NonElectrical.Continuous.SimpleLag Speed_Sensor(K=1, T=Tr,
-    y_start=Modelica.Constants.eps)
+    y_start=1.9*Modelica.Constants.pi*SysData.fn)
     annotation (Placement(transformation(extent={{-38,-22},{-18,-2}})));
   parameter OpenIPSL.Types.Time Tr=0.01 "Time constant for speed sensor filter"
     annotation (Dialog(group="Control Parameters"));
@@ -36,7 +36,7 @@ model VoltsHertzController "Volts/Hertz controller model"
     annotation (Placement(transformation(extent={{-48,-60},{-28,-40}})));
   Modelica.Blocks.Continuous.Integrator integrator(k=Ki,
     initType=Modelica.Blocks.Types.Init.InitialState,
-    y_start=m0*(2*Modelica.Constants.pi*SysData.fn))
+    y_start=0)
     annotation (Placement(transformation(extent={{-48,-92},{-28,-72}})));
   Modelica.Blocks.Math.Add add1(k1=+1)
     annotation (Placement(transformation(extent={{-14,-76},{6,-56}})));
diff --git a/OpenIPSL/Electrical/Controls/VariableSpeedDrive/PowerElectronics/AC2DCandDC2AC.mo b/OpenIPSL/Electrical/Controls/VariableSpeedDrive/PowerElectronics/AC2DCandDC2AC.mo
@@ -27,7 +27,7 @@ model AC2DCandDC2AC "Phasor based voltage source converter model."
             {-64,-4}})));
   Modelica.Electrical.Analog.Basic.Resistor Resistor(R=Rdc) "DC link resistor"
     annotation (Placement(transformation(extent={{-42,-10},{-22,10}})));
-  Modelica.Electrical.Analog.Basic.Inductor Inductor(i(start=Il0, fixed=true),
+  Modelica.Electrical.Analog.Basic.Inductor Inductor(i(start=Il0, fixed=false),
                                                      L=Ldc) "DC link inductor"
     annotation (Placement(transformation(extent={{-16,-10},{4,10}})));
   Modelica.Electrical.Analog.Ideal.IdealOpeningSwitch switch(Ron=1e-5, Goff=1e-5) "Switch enforcing no inverse current flow in DC link"
@@ -58,6 +58,7 @@ model AC2DCandDC2AC "Phasor based voltage source converter model."
   Modelica.Blocks.Sources.RealExpression Qmotor(y=n.vr*n.ii - n.vi*n.ir) "Reactive power for drawn motor current"
     annotation (Placement(transformation(extent={{24,54},{44,74}})));
   OpenIPSL.Types.PerUnit P "Active Power";
+  Modelica.Units.SI.ActivePower Pdc;
   OpenIPSL.Types.PerUnit Q "Reactive Power";
   OpenIPSL.Types.PerUnit S "Motor Slip";
 
@@ -131,7 +132,8 @@ equation
     Q = (-p.vr*p.ii) + p.vi*p.ir;
     Q = 0;
     S = sqrt(P^2 + Q^2);
-    Resistor.i = smooth(0,(P*S_b)/Vd0.y);
+    Pdc = Vd0.y*Resistor.i;
+    P*S_b = Pdc;
 
     n.vr = vr_m.y;
     n.vi = vi_m.y;
diff --git a/OpenIPSL/Electrical/Machines/PSSE/BaseClasses/baseMotor.mo b/OpenIPSL/Electrical/Machines/PSSE/BaseClasses/baseMotor.mo
@@ -94,8 +94,8 @@ equation
   v = sqrt(p.vr^2 + p.vi^2);
   anglev = atan2(p.vi, p.vr);
   delta = anglev;
-  P = p.vr*p.ir + p.vi*p.ii;
-  Q = (-p.vr*p.ii) + p.vi*p.ir;
+  P = if Ctrl == true then (p.vr*p.ir + p.vi*p.ii)*(we_fix.y/w_b) else (p.vr*p.ir + p.vi*p.ii);
+  Q = if Ctrl == true then ((-p.vr*p.ii) + p.vi*p.ir)*(we_fix.y/w_b) else ((-p.vr*p.ii) + p.vi*p.ir);
   P_motor = P/CoB;
   Q_motor = Q/CoB;
   ns = w_sync/N;
diff --git a/OpenIPSL/Electrical/Renewables/PSSE/BESS.mo b/OpenIPSL/Electrical/Renewables/PSSE/BESS.mo
@@ -18,6 +18,7 @@ extends OpenIPSL.Electrical.Essentials.pfComponent(
   replaceable
     OpenIPSL.Electrical.Renewables.PSSE.InverterInterface.BaseClasses.BaseREGC
     RenewableGenerator(
+    V_b=V_b,
     M_b=M_b,
     P_0=P_0,
     Q_0=Q_0,
@@ -37,6 +38,7 @@ extends OpenIPSL.Electrical.Essentials.pfComponent(
   replaceable
     OpenIPSL.Electrical.Renewables.PSSE.PlantController.BaseClasses.BaseREPC
     PlantController(
+    V_b=V_b,
     M_b=M_b,
     P_0=P_0,
     Q_0=Q_0,
diff --git a/OpenIPSL/Electrical/Renewables/PSSE/PV.mo b/OpenIPSL/Electrical/Renewables/PSSE/PV.mo
@@ -18,6 +18,7 @@ extends OpenIPSL.Electrical.Essentials.pfComponent(
   replaceable
     OpenIPSL.Electrical.Renewables.PSSE.InverterInterface.BaseClasses.BaseREGC
     RenewableGenerator(
+    V_b=V_b,
     M_b=M_b,
     P_0=P_0,
     Q_0=Q_0,
@@ -37,6 +38,7 @@ extends OpenIPSL.Electrical.Essentials.pfComponent(
   replaceable
     OpenIPSL.Electrical.Renewables.PSSE.PlantController.BaseClasses.BaseREPC
     PlantController(
+    V_b=V_b,
     M_b=M_b,
     P_0=P_0,
     Q_0=Q_0,
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/Data/LoadsMicrogrid.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/Data/LoadsMicrogrid.mo
@@ -4,7 +4,7 @@ record LoadsMicrogrid "Data record for load power consumption"
 
   // Load 1
   parameter Real P1=0.02e6;
-  parameter Real Q1=0.005e6;
+  parameter Real Q1=0.0123e6;
 
 
   annotation (Documentation(info="<html>
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/Data/MachinesMicrogrid.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/Data/MachinesMicrogrid.mo
@@ -3,12 +3,22 @@ record MachinesMicrogrid "Data record for power generation"
   extends Modelica.Icons.Record;
 
   // Machine Diesel Unit
-  parameter Real PDT=0.05e6;
-  parameter Real QDT=-0.04e6;
+  parameter Real PDT=0.02e6;
+  parameter Real QDT=-0.0313e6;
+
+  // PV
+  parameter Real PPV=0.035e6;
+  parameter Real QPV=0.0115e6;
+
+  // BESS
+  parameter Real PBESS=0.0055e6;
+  parameter Real QBESS=0;
 
   // Machine Inf
-  parameter Real PInf=-0.02999824e6;
-  parameter Real QInf=-0.1550323e6;
+  parameter Real PInf=-0.0405e6;
+  parameter Real QInf= 0.0121e6;
+
+
 
   annotation (Documentation(info="<html>
 This data record has the active and reactive powers in VA and Var, respectively, that are delivered by the generation units to the system. The data is used to initialized the network in a certain operating point.
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/Data/VoltagesMicrogrid.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/Data/VoltagesMicrogrid.mo
@@ -9,24 +9,29 @@ record VoltagesMicrogrid
   parameter Real A1=0*pi/180;
 
   // Bus LVBus
-  parameter Real V2=1.0000;
-  parameter Real A2=0*pi/180;
+  parameter Real V2=1.00;
+  parameter Real A2=0.001322568*pi/180;
 
   // Bus SubBus
-  parameter Real V3=1.000034;
-  parameter Real A3=-0.0006039513*pi/180;
+  parameter Real V3=1.00;
+  parameter Real A3=0.001845621*pi/180;
 
   // Bus CentralBus
-  parameter Real V4=1.000072;
-  parameter Real A4=-0.001997927*pi/180;
+  parameter Real V4=1.00;
+  parameter Real A4=0.002433209*pi/180;
 
   // Bus LoadBus
-  parameter Real V5=1.000013;
-  parameter Real A5=-0.002251726*pi/180;
+  parameter Real V5=1.00;
+  parameter Real A5=0.002433209*pi/180;
 
   // Bus CapBus
-  parameter Real V6=1.000072;
-  parameter Real A6=-0.001997927*pi/180;
+  parameter Real V6=1.00;
+  parameter Real A6=0.003302665*pi/180;
+
+  // Bus Motor
+  parameter Real V7=1.00;
+  parameter Real A7=0.002433209*pi/180;
+
 
   annotation (Documentation(info="<html>
 This data record has the voltage magnitude, in per unit, and angle, in radians, that are used to initialized the network in a certain operating point.
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/GeneratorGroups/DieselGeneratorUnit.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/GeneratorGroups/DieselGeneratorUnit.mo
@@ -7,53 +7,30 @@ model DieselGeneratorUnit
     Q_0=Q_0,
     v_0=v_0,
     angle_0=angle_0,
-    M_b(displayUnit="kW") = 100000,
+    M_b(displayUnit="V.A") = 50000,
     Tpd0=3.7,
     Tppd0=0.05,
     Tppq0=0.05,
     H=1.07,
     D=0,
-    Xd=1.56,
-    Xq=1.06,
+    Xd=1,
+    Xq=1,
     Xpd=0.296,
     Xppd=0.177,
     Xppq=0.177,
     Xl=0.052,
     S10=0.12,
     S12=0.4,
-    R_a=0.0036,
+    R_a=0,
     Xpq=0.4610,
     Tpq0=0.391)
     annotation (Placement(transformation(extent={{20,-20},{60,20}})));
-  Electrical.Controls.PSSE.ES.ESAC1A eSAC1A(
-    T_R=0.02,
-    T_B=0,
-    T_C=0,
-    K_A=400,
-    T_A=0.02,
-    V_AMAX=14.5,
-    V_AMIN=-14.5,
-    T_E=0.8,
-    K_F=0.03,
-    T_F=1,
-    K_C=0.2,
-    K_D=0.38,
-    K_E=1.0,
-    E_1=4.18,
-    E_2=5.6,
-    S_EE_1=0.1827,
-    S_EE_2=0.255,
-    V_RMAX=6.03,
-    V_RMIN=-5.43)
+  Electrical.Controls.PSSE.ES.SEXS   sEXS
     annotation (Placement(transformation(extent={{-20,-30},{0,-10}})));
   Modelica.Blocks.Sources.Constant zero(k=0) annotation (Placement(transformation(
         extent={{6,6},{-6,-6}},
         rotation=180,
         origin={-66,-16})));
-  Modelica.Blocks.Sources.Constant minusInf(k=-Modelica.Constants.inf) annotation (Placement(transformation(extent={{-60,-74},
-            {-52,-66}})));
-  Modelica.Blocks.Sources.Constant plusInf(k=Modelica.Constants.inf) annotation (Placement(transformation(extent={{-60,-88},
-            {-52,-80}})));
   Electrical.Controls.PSSE.TG.DEGOV dEGOV(
     T1=0.01,
     T2=0.02,
@@ -68,26 +45,26 @@ model DieselGeneratorUnit
 equation
   connect(gENROE.p, pwPin)
     annotation (Line(points={{60,0},{110,0}}, color={0,0,255}));
-  connect(eSAC1A.EFD, gENROE.EFD) annotation (Line(points={{1,-20},{6,-20},{6,
-          -12},{16,-12}}, color={0,0,127}));
-  connect(zero.y, eSAC1A.VOTHSG)
+  connect(sEXS.EFD, gENROE.EFD) annotation (Line(points={{1,-20},{6,-20},{6,-12},
+          {16,-12}}, color={0,0,127}));
+  connect(zero.y, sEXS.VOTHSG)
     annotation (Line(points={{-59.4,-16},{-21,-16}}, color={0,0,127}));
-  connect(minusInf.y, eSAC1A.VUEL) annotation (Line(points={{-51.6,-70},{-14,
-          -70},{-14,-31}}, color={0,0,127}));
-  connect(plusInf.y, eSAC1A.VOEL) annotation (Line(points={{-51.6,-84},{-10,-84},
-          {-10,-31}}, color={0,0,127}));
-  connect(eSAC1A.XADIFD, gENROE.XADIFD) annotation (Line(points={{-2,-31},{-2,
-          -40},{70,-40},{70,-18},{62,-18}}, color={0,0,127}));
-  connect(eSAC1A.EFD0, gENROE.EFD0) annotation (Line(points={{-21,-24},{-26,-24},
+  connect(sEXS.XADIFD, gENROE.XADIFD) annotation (Line(points={{-2,-31},{-2,-40},
+          {70,-40},{70,-18},{62,-18}}, color={0,0,127}));
+  connect(sEXS.EFD0, gENROE.EFD0) annotation (Line(points={{-21,-24},{-26,-24},
           {-26,-44},{74,-44},{74,-10},{62,-10}}, color={0,0,127}));
-  connect(eSAC1A.ECOMP, gENROE.ETERM) annotation (Line(points={{-21,-20},{-30,
-          -20},{-30,-48},{78,-48},{78,-6},{62,-6}}, color={0,0,127}));
-  connect(dEGOV.PMECH, gENROE.PMECH)
+  connect(sEXS.ECOMP, gENROE.ETERM) annotation (Line(points={{-21,-20},{-30,-20},
+          {-30,-48},{78,-48},{78,-6},{62,-6}}, color={0,0,127}));
+  connect(dEGOV.PMECH,gENROE. PMECH)
     annotation (Line(points={{1,12},{16,12}}, color={0,0,127}));
   connect(gENROE.SPEED, dEGOV.SPEED) annotation (Line(points={{62,14},{70,14},{
           70,28},{-24,28},{-24,18},{-18,18}}, color={0,0,127}));
   connect(gENROE.PMECH0, dEGOV.PMECH0) annotation (Line(points={{62,10},{74,10},
           {74,32},{-28,32},{-28,6},{-18,6}}, color={0,0,127}));
+  connect(sEXS.VUEL, zero.y) annotation (Line(points={{-14,-31},{-14,-38},{-36,
+          -38},{-36,-16},{-59.4,-16}}, color={0,0,127}));
+  connect(sEXS.VOEL, zero.y) annotation (Line(points={{-10,-31},{-10,-38},{-36,
+          -38},{-36,-16},{-59.4,-16}}, color={0,0,127}));
   annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(
         coordinateSystem(preserveAspectRatio=false)));
 end DieselGeneratorUnit;
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/IEEEMicrogrid.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/IEEEMicrogrid.mo
@@ -157,16 +157,16 @@ model IEEEMicrogrid
         origin={230,100})));
   Electrical.Buses.Bus CapBus(
     V_b=400,
-    v_0=PF_results.voltages.V6,
-    angle_0=PF_results.voltages.A6,
+    v_0=PF_results.voltages.V5,
+    angle_0=PF_results.voltages.A5,
     displayPF=false) annotation (Placement(transformation(
         extent={{-16,-16},{16,16}},
         rotation=90,
         origin={-130,-70})));
   Electrical.Buses.Bus LoadBus(
     V_b=400,
-    v_0=PF_results.voltages.V5,
-    angle_0=PF_results.voltages.A5,
+    v_0=PF_results.voltages.V6,
+    angle_0=PF_results.voltages.A6,
     displayPF=false) annotation (Placement(transformation(
         extent={{-16,-16},{16,16}},
         rotation=90,
@@ -202,32 +202,57 @@ model IEEEMicrogrid
         origin={140,-2})));
   Electrical.Buses.Bus MotorBus(
     V_b=400,
-    v_0=PF_results.voltages.V5,
-    angle_0=PF_results.voltages.A5,
+    v_0=PF_results.voltages.V7,
+    angle_0=PF_results.voltages.A7,
     displayPF=false) annotation (Placement(transformation(
         extent={{-16,-16},{16,16}},
         rotation=90,
         origin={140,-70})));
   Electrical.Controls.VariableSpeedDrive.PowerElectronics.AC2DCandDC2AC
-    aC2DCandDC2AC
+    aC2DCandDC2AC(
+    V_b=400,
+    Cdc=0.000001,
+    m0=0.9)
     annotation (Placement(transformation(extent={{100,-126},{80,-106}})));
   Electrical.Controls.VariableSpeedDrive.ControllerLogic.VoltsHertzController
-    voltsHertzController
+    voltsHertzController(
+    V_b=400,
+    f_max=60,
+    f_min=0,
+    m0=0.9,
+    Kp=0.5,
+    Ki=0.2)
     annotation (Placement(transformation(extent={{100,-166},{80,-146}})));
-  Electrical.Machines.PSSE.CIM5 Motor1(M_b=50000)
+  Electrical.Machines.PSSE.CIM5 Motor1(
+    V_b=400,
+    M_b=50000,
+    Sup=false)
     annotation (Placement(transformation(extent={{40,-126},{60,-106}})));
   Electrical.Machines.PSSE.CIM5 Motor2(
+    V_b=400,
     M_b=50000,                         Sup=false, Ctrl=false)
     annotation (Placement(transformation(extent={{220,-126},{200,-106}})));
   Modelica.Blocks.Sources.Ramp Sync_Speed(
-    height=0.9*(2*Modelica.Constants.pi*SysData.fn),
+    height=0*(2*Modelica.Constants.pi*SysData.fn),
     duration=5,
-    offset=0.1*(2*Modelica.Constants.pi*SysData.fn),
+    offset=0.9*(1.9*Modelica.Constants.pi*SysData.fn),
     startTime=0)
     annotation (Placement(transformation(extent={{140,-166},{120,-146}})));
-  Electrical.Renewables.PSSE.PV PV
+  Electrical.Renewables.PSSE.PV PV(
+    redeclare OpenIPSL.Electrical.Renewables.PSSE.InverterInterface.REGCA1
+      RenewableGenerator,
+    redeclare OpenIPSL.Electrical.Renewables.PSSE.ElectricalController.REECB1
+      RenewableController,
+    redeclare OpenIPSL.Electrical.Renewables.PSSE.PlantController.REPCA1
+      PlantController)
     annotation (Placement(transformation(extent={{240,120},{220,140}})));
-  Electrical.Renewables.PSSE.BESS BESS
+  Electrical.Renewables.PSSE.BESS BESS(
+    redeclare OpenIPSL.Electrical.Renewables.PSSE.InverterInterface.REGCA1
+      RenewableGenerator,
+    redeclare OpenIPSL.Electrical.Renewables.PSSE.ElectricalController.REECCU1
+      RenewableController,
+    redeclare OpenIPSL.Electrical.Renewables.PSSE.PlantController.REPCA1
+      PlantController)
     annotation (Placement(transformation(extent={{240,60},{220,80}})));
 equation
   connect(Transformer.p, LVBus.p)
@@ -343,6 +368,6 @@ equation
           textString="(7)")}),
     experiment(
       StopTime=20,
-      __Dymola_NumberOfIntervals=5000,
+      __Dymola_NumberOfIntervals=1000,
       __Dymola_Algorithm="Dassl"));
 end IEEEMicrogrid;
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/IEEEMicrogrid_2.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/IEEEMicrogrid_2.mo
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/IEEEMicrogrid_3.mo b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/IEEEMicrogrid_3.mo
diff --git a/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/package.order b/OpenIPSL/Examples/Microgrids/IEEEMicrogrid/package.order
diff --git a/OpenIPSL/Tests/Controls/PSSE/ES/SEXS.mo b/OpenIPSL/Tests/Controls/PSSE/ES/SEXS.mo
