Emt syngen trstab

dpsim-models/include/dpsim-models/Base/Base_Ph3_Switch.h

@@ -36,6 +36,9 @@ namespace Ph3 {
 		}
 		void closeSwitch() { **mSwitchClosed = true; }
 		void openSwitch() { **mSwitchClosed = false; }
+
+		/// Check if switch is closed
+		Bool isClosed() { return **mSwitchClosed; }
 	};
 }
 }

dpsim-models/include/dpsim-models/Components.h

@@ -115,6 +115,7 @@
 #include <dpsim-models/EMT/EMT_Ph3_Transformer.h>
 #include <dpsim-models/EMT/EMT_Ph3_NetworkInjection.h>
 #include <dpsim-models/EMT/EMT_Ph3_Switch.h>
+#include <dpsim-models/EMT/EMT_Ph3_varResSwitch.h>
 #include <dpsim-models/EMT/EMT_Ph3_SynchronGeneratorTrStab.h>
 
 #include <dpsim-models/Signal/DecouplingLine.h>

dpsim-models/include/dpsim-models/EMT/EMT_Ph3_SynchronGeneratorTrStab.h

@@ -42,7 +42,11 @@ namespace Ph3 {
 		Matrix mStates;
 		/// Nominal system angle
 		Real mThetaN = 0;
-
+		/// Flag for usage of attribute of w_ref (otherwise mNomOmega is used)
+		Bool mUseOmegaRef = false;
+		/// Flag for usage of actual mechanical speed for torque conversion (otherwise mNomOmega is used)
+		Bool mConvertWithOmegaMech = true;
+
 	public:
 		// #### Model specific variables ####
 		/// emf behind transient reactance
@@ -73,6 +77,8 @@ namespace Ph3 {
 		Matrix getParkTransformMatrixPowerInvariant(Real theta);
 
 		// #### General Functions ####
+		/// Flags to modify model behavior
+		void setModelFlags(Bool convertWithOmegaMech);
 		///
 		void setInitialValues(Complex elecPower, Real mechPower);
 		/// \brief Initializes the machine parameters
@@ -105,6 +111,8 @@ namespace Ph3 {
 		void mnaParentAddPreStepDependencies(AttributeBase::List &prevStepDependencies, AttributeBase::List &attributeDependencies, AttributeBase::List &modifiedAttributes) override;
 		void mnaParentAddPostStepDependencies(AttributeBase::List &prevStepDependencies, AttributeBase::List &attributeDependencies, AttributeBase::List &modifiedAttributes, Attribute<Matrix>::Ptr &leftVector) override;
 
+		void setReferenceOmega(Attribute<Real>::Ptr refOmegaPtr, Attribute<Real>::Ptr refDeltaPtr);
+
 		class AddBStep : public Task {
 		public:
 			AddBStep(SynchronGeneratorTrStab& generator) :

dpsim-models/include/dpsim-models/EMT/EMT_Ph3_varResSwitch.h

@@ -0,0 +1,89 @@
+/* Copyright 2017-2020 Institute for Automation of Complex Power Systems,
+ *                     EONERC, RWTH Aachen University
+ *
+ * This Source Code Form is subject to the terms of the Mozilla Public
+ * License, v. 2.0. If a copy of the MPL was not distributed with this
+ * file, You can obtain one at https://mozilla.org/MPL/2.0/.
+ *********************************************************************************/
+
+#pragma once
+
+#include <dpsim-models/Solver/MNAInterface.h>
+#include <dpsim-models/MNASimPowerComp.h>
+#include <dpsim-models/Solver/MNASwitchInterface.h>
+#include <dpsim-models/Solver/MNAVariableCompInterface.h>
+#include <dpsim-models/Definitions.h>
+#include <dpsim-models/Logger.h>
+#include <dpsim-models/Base/Base_Ph3_Switch.h>
+
+namespace CPS {
+namespace EMT {
+namespace Ph3 {
+	/// \brief
+	///
+	/// Switch with variable resistance to avoid numerical oscillations, when an inductive current is suddenly interrupted.
+	/// It is useful as a fault switch especially on a faulted generator or transformer winding.
+
+	/// The switch resistance changes at a defined fixed rate by multiplying previous resistance value with a factor for the rate of change
+	/// The MNA variable component interface is used to recompute the system Matrix after each change.
+	class varResSwitch :
+		public MNASimPowerComp<Real>,
+		public Base::Ph3::Switch,
+		public MNAVariableCompInterface,
+		public MNASwitchInterface,
+		public SharedFactory<varResSwitch> {
+
+	protected:
+
+		Bool mPrevState=false;
+		Real mDeltaResClosed = 0; 
+		Real mDeltaResOpen = 1.5;
+		Matrix mPrevRes; // previous resistance value to multiply with rate of change
+		// because we change the base value mClosedResistance & mOpenResistance to recompute the system Matrix
+		// we need to save the initialisation values to use them as target values in the transition
+		Matrix mInitClosedRes;
+		Matrix mInitOpenRes;
+
+
+
+	public:
+		/// Defines UID, name, component parameters and logging level
+		varResSwitch(String uid, String name, Logger::Level logLevel = Logger::Level::off);
+		/// Defines name, component parameters and logging level
+		varResSwitch(String name, Logger::Level logLevel = Logger::Level::off)
+			: varResSwitch(name, name, logLevel) { }
+
+		SimPowerComp<Real>::Ptr clone(String name);
+
+		// #### General ####
+		void setInitParameters(Real timestep);
+		/// Initializes component from power flow data
+		void initializeFromNodesAndTerminals(Real frequency);
+
+		// #### General MNA section ####
+		void mnaCompInitialize(Real omega, Real timeStep, Attribute<Matrix>::Ptr leftVector);
+		/// Stamps system matrix
+		void mnaCompApplySystemMatrixStamp(SparseMatrixRow& systemMatrix);
+		/// Stamps right side (source) vector
+		void mnaCompApplyRightSideVectorStamp(Matrix& rightVector);
+		/// Update interface voltage from MNA system result
+		void mnaCompUpdateVoltage(const Matrix& leftVector);
+		/// Update interface current from MNA system result
+		void mnaCompUpdateCurrent(const Matrix& leftVector);
+
+		// #### MNA section for switches ####
+		/// Check if switch is closed
+		Bool mnaIsClosed() { return isClosed(); }
+		/// Stamps system matrix considering the defined switch position
+		void mnaCompApplySwitchSystemMatrixStamp(Bool closed, SparseMatrixRow& systemMatrix, Int freqIdx);
+
+		void mnaCompPostStep(Real time, Int timeStepCount, Attribute<Matrix>::Ptr &leftVector) override;
+
+		/// Add MNA post step dependencies
+		void mnaCompAddPostStepDependencies(AttributeBase::List &prevStepDependencies, AttributeBase::List &attributeDependencies, AttributeBase::List &modifiedAttributes, Attribute<Matrix>::Ptr &leftVector) override;
+
+		Bool hasParameterChanged();
+	};
+}
+}
+}

dpsim-models/src/CMakeLists.txt

@@ -98,6 +98,7 @@ list(APPEND MODELS_SOURCES
 	EMT/EMT_Ph3_SynchronGeneratorIdeal.cpp
 	# EMT/EMT_Ph3_SynchronGeneratorVBRStandalone.cpp
 	EMT/EMT_Ph3_SynchronGeneratorTrStab.cpp
+	EMT/EMT_Ph3_varResSwitch.cpp
 
 	SP/SP_Ph1_VoltageSource.cpp
 	SP/SP_Ph1_Capacitor.cpp

dpsim-models/src/DP/DP_Ph1_SynchronGeneratorTrStab.cpp

@@ -176,15 +176,15 @@ void DP::Ph1::SynchronGeneratorTrStab::initializeFromNodesAndTerminals(Real freq
 	mSubVoltageSource->setVirtualNodeAt(mVirtualNodes[1], 0);
 	mSubVoltageSource->initialize(mFrequencies);
 	mSubVoltageSource->initializeFromNodesAndTerminals(frequency);
-	addMNASubComponent(mSubVoltageSource, MNA_SUBCOMP_TASK_ORDER::TASK_BEFORE_PARENT, MNA_SUBCOMP_TASK_ORDER::TASK_BEFORE_PARENT, true);
+	addMNASubComponent(mSubVoltageSource, MNA_SUBCOMP_TASK_ORDER::TASK_AFTER_PARENT, MNA_SUBCOMP_TASK_ORDER::TASK_BEFORE_PARENT, true);
 
 	// Create sub inductor as Xpd
 	mSubInductor = DP::Ph1::Inductor::make(**mName + "_ind", mLogLevel);
 	mSubInductor->setParameters(mLpd);
 	mSubInductor->connect({mVirtualNodes[0], terminal(0)->node()});
 	mSubInductor->initialize(mFrequencies);
 	mSubInductor->initializeFromNodesAndTerminals(frequency);
-	addMNASubComponent(mSubInductor, MNA_SUBCOMP_TASK_ORDER::TASK_BEFORE_PARENT, MNA_SUBCOMP_TASK_ORDER::TASK_BEFORE_PARENT, true);
+	addMNASubComponent(mSubInductor, MNA_SUBCOMP_TASK_ORDER::TASK_AFTER_PARENT, MNA_SUBCOMP_TASK_ORDER::TASK_BEFORE_PARENT, true);
 
 	SPDLOG_LOGGER_INFO(mSLog, "\n--- Initialize according to powerflow ---"
 				"\nTerminal 0 voltage: {:e}<{:e}"