Added f_grid - grid frequency to Data, corrected reference speed for normalization of speed (to p.u.) in Power2Torque and TorqueEquation. Added test case for TorqueEquation in OpenHPLTest

Author: boerrebj

OpenHPL/Data.mo

@@ -19,6 +19,7 @@ record Data "Provides a data set of most common used settings"
     annotation (Dialog(group = "Waterway properties"));
   parameter SI.Compressibility beta_total = 1 / (rho*1000^2) "Total compressibility"
    annotation (Dialog(group = "Waterway properties"));
+  parameter SI.Frequency f_grid=50 "Grid frequency"  annotation (Dialog(group = "System properties"));
   parameter Boolean SteadyState=false "If checked, simulation starts in steady state"
     annotation (choices(checkBox = true), Dialog(group="Initialization"));
   parameter SI.VolumeFlowRate Vdot_0 = 0 "Initial volume flow rate through the system"

OpenHPL/ElectroMech/BaseClasses/Power2Torque.mo

@@ -67,7 +67,7 @@ partial model Power2Torque "Converts a power signal to a torque in the rotationa
         extent={{-6,-6},{6,6}},
         rotation=180,
         origin={-50,-40})));
-  Modelica.Blocks.Math.Gain w_m2pu(k=(p/2)/(2*C.pi*data.f_0))
+  Modelica.Blocks.Math.Gain w_m2pu(k=(p/2)/(2*C.pi*data.f_grid))
     annotation (Placement(transformation(extent={{66,-46},{78,-34}})));
   Modelica.Blocks.Nonlinear.Limiter torqueLimit(uMax=Pmax/f_0)
     annotation (Placement(transformation(
@@ -101,7 +101,7 @@ partial model Power2Torque "Converts a power signal to a torque in the rotationa
   Modelica.Blocks.Interfaces.RealInput f_in if enable_f_in and not enable_nomSpeed
     "Speed input of the unit [pu]"
     annotation (Placement(transformation(extent={{-20,-20},{20,20}},rotation=90,origin={-80,-120})));
-  Modelica.Blocks.Math.Gain pu2w_s(k=2*C.pi*data.f_0) if enable_f_in or enable_nomSpeed
+  Modelica.Blocks.Math.Gain pu2w_s(k=2*C.pi*data.f_grid) if enable_f_in or enable_nomSpeed
     annotation (Placement(transformation(extent={{40,-90},{60,-70}})));
 
 equation

OpenHPL/ElectroMech/BaseClasses/TorqueEquation.mo

@@ -44,7 +44,7 @@ partial model TorqueEquation
     Placement(transformation(extent = {{0, 60}, {20, 40}})));
   Modelica.Mechanics.Rotational.Components.Fixed fixed annotation(
     Placement(transformation(extent = {{20, 50}, {40, 70}})));
- Modelica.Blocks.Math.Gain w_m2pu(k = (p/2)/(2*Modelica.Constants.pi*data.f_0)) annotation(
+  Modelica.Blocks.Math.Gain w_m2pu(k = (p/2)/(data.f_grid*2*Modelica.Constants.pi)) "Convert from rad/s to pu" annotation(
     Placement(transformation(extent = {{66, -46}, {78, -34}})));
   Modelica.Blocks.Interfaces.RealOutput f if enable_f "Speed output of the unit [pu]" annotation(
     Placement(transformation(extent = {{100, -50}, {120, -30}}), iconTransformation(extent = {{100, -50}, {120, -30}})));
@@ -64,12 +64,12 @@ partial model TorqueEquation
     Placement(transformation(extent = {{-10, -80}, {10, -60}})));
   Modelica.Blocks.Interfaces.RealInput f_in if enable_f_in and not enable_nomSpeed "Speed input of the unit [pu]" annotation(
     Placement(transformation(extent = {{-20, -20}, {20, 20}}, rotation = 90, origin = {-80, -120})));
-  Modelica.Blocks.Math.Gain pu2w_s(k = 2*Modelica.Constants.pi*data.f_0) if enable_f_in or enable_nomSpeed annotation(
+  Modelica.Blocks.Math.Gain pu2w_s(k = 2*Modelica.Constants.pi*data.f_grid) if enable_f_in or enable_nomSpeed annotation(
     Placement(transformation(extent = {{40, -90}, {60, -70}})));
   protected
    Modelica.Mechanics.Rotational.Sources.Torque torque_transfer annotation(
     Placement(transformation(extent = {{-36, -6}, {-24, 6}})));
-  
+   
 equation
   connect(w, speedSensor.w) annotation(
     Line(points = {{110, 40}, {40, 40}, {40, -42}, {10, -42}, {10, -39}}, color = {0, 0, 127}, pattern = LinePattern.Dash));
@@ -105,6 +105,6 @@ equation
     Line(points = {{-80, 30}, {-94, 30}, {-94, 0}, {-38, 0}}, color = {0, 0, 127}));
   annotation(
     Icon(graphics = {Text(visible = enable_w, extent = {{80, 50}, {100, 30}}, textColor = {0, 0, 0}, textString = "w"), Text(visible = enable_f, extent = {{80, -30}, {100, -50}}, textColor = {0, 0, 0}, textString = "f"), Text(visible = enable_f_in, extent = {{-100, -70}, {-60, -90}}, textColor = {0, 0, 0}, textString = "f_in")}),
-Documentation(info = "<html><head></head><body><p>Abstract (partial) base class for including the torque equation:<br> $$ J\frac{d\omega}{dt}=T $$ <br>In the future this base class can replace <font face=\"Courier\">Power2Torque</font> in the turbine models to avoid the issue at zero speed. This class is also better suited for fundamental or mechanistic turbine models.</p></body></html>"));
+Documentation(info = "<html><head></head><body><p>Abstract (partial) base class for including the torque equation:<br> $$ J\\frac{d\\omega}{dt}=T $$ <br>In the future this base class can replace <font face=\"Courier\">Power2Torque</font> in the turbine models to avoid the issue at zero speed. This class is also better suited for fundamental or mechanistic turbine models.</p></body></html>"));
 
 end TorqueEquation;

OpenHPLTest/TorqueEquation.mo

@@ -2,4 +2,36 @@ within OpenHPLTest;
 
 package TorqueEquation
   extends Modelica.Icons.ExamplesPackage;
+  import SI = Modelica.Units.SI;
+  //
+
+  model TorqueElement
+     extends OpenHPL.Icons.ElectroMech;
+     extends OpenHPL.ElectroMech.BaseClasses.TorqueEquation;
+  equation
+
+  end TorqueElement;
+
+  model TorqueTest
+    extends Modelica.Icons.Example;
+    //
+    parameter SI.Torque shaftTorque0 = 1.e+03;
+    SI.Torque shaftTorque;
+    inner OpenHPL.Data data annotation(
+      Placement(transformation(origin = {-52, 74}, extent = {{-10, -10}, {10, 10}})));
+    TorqueElement te1(J = 10, f_0 = 0, torque(y = shaftTorque), enable_f = true, p = 10) annotation(
+      Placement(transformation(origin = {-30, 52}, extent = {{-10, -10}, {10, 10}})));
+    TorqueElement te2(J = 10, f_0 = 0, torque(y = shaftTorque), enable_f = true, p = 22) annotation(
+      Placement(transformation(origin = {-30, 24}, extent = {{-10, -10}, {10, 10}})));
+    equation
+    if (time > 0.1 and time < 0.4) then
+      shaftTorque = shaftTorque0;
+    elseif (time > 0.5 and time < 0.8) then
+      shaftTorque = -shaftTorque0;
+    else
+      shaftTorque = 0.0;
+    end if;
+  annotation(
+      Diagram(coordinateSystem(extent = {{-80, 80}, {-20, 0}})));
+end TorqueTest;
 end TorqueEquation;