EMT_Ph3_SynchronGeneratorDQ.h

/* Copyright 2017-2021 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/Base/Base_SynchronGenerator.h>
#include <dpsim-models/MNASimPowerComp.h>
#include <dpsim-models/Solver/MNAInterface.h>
 
namespace CPS {
namespace EMT {
namespace Ph3 {
class SynchronGeneratorDQ : public MNASimPowerComp<Real>,
                            public Base::SynchronGenerator {
protected:
  Matrix mCompensationCurrent;

  SynchronGeneratorDQ(String name, String uid,
                      Logger::Level logLevel = Logger::Level::off);
  SynchronGeneratorDQ(String name, Logger::Level logLevel = Logger::Level::off);

  Matrix abcToDq0Transform(Real theta, Matrix &abc);

  Matrix dq0ToAbcTransform(Real theta, Matrix &dq0);
 
public:
  const std::vector<String> attrParamNames = {
      "Rs",   "Ll",   "Ld",    "Lq",    "Ld_t",  "Ld_s",
      "Lq_t", "Lq_s", "Td0_t", "Td0_s", "Tq0_t", "Tq0_s"};
 
  virtual ~SynchronGeneratorDQ();

  void setParametersFundamentalPerUnit(
      Real nomPower, Real nomVolt, Real nomFreq, Int poleNumber,
      Real nomFieldCur, Real Rs, Real Ll, Real Lmd, Real Lmq, Real Rfd,
      Real Llfd, Real Rkd, Real Llkd, Real Rkq1, Real Llkq1, Real Rkq2,
      Real Llkq2, Real inertia, Real initActivePower, Real initReactivePower,
      Real initTerminalVolt, Real initVoltAngle, Real initMechPower);

  void setParametersOperationalPerUnit(Real nomPower, Real nomVolt,
                                       Real nomFreq, Int poleNumber,
                                       Real nomFieldCur, Real Rs, Real Ld,
                                       Real Lq, Real Ld_t, Real Lq_t, Real Ld_s,
                                       Real Lq_s, Real Ll, Real Td0_t,
                                       Real Tq0_t, Real Td0_s, Real Tq0_s,
                                       Real inertia);

  void setInitialValues(Real initActivePower, Real initReactivePower,
                        Real initTerminalVolt, Real initVoltAngle,
                        Real initMechPower);

  void applyParametersOperationalPerUnit();

  void initializeFromNodesAndTerminals(Real frequency) override;
  void initializeMatrixAndStates();

  void initialize(Matrix frequencies) override;
  Real electricalTorque() const;
  Real rotationalSpeed() const;
  Real rotorPosition() const;

  void step(Matrix &voltage, Real time);
 
  // #### MNA Functions ####

  void mnaCompApplyRightSideVectorStamp(Matrix &rightVector) override;
  void mnaCompApplySystemMatrixStamp(SparseMatrixRow &systemMatrix) override;

  virtual void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaCompPostStep(Real time, Int timeStepCount,
                       Attribute<Matrix>::Ptr &leftVector) override;

  void
  mnaCompAddPostStepDependencies(AttributeBase::List &prevStepDependencies,
                                 AttributeBase::List &attributeDependencies,
                                 AttributeBase::List &modifiedAttributes,
                                 Attribute<Matrix>::Ptr &leftVector) override;
};
} // namespace Ph3
} // namespace EMT
} // namespace CPS