EMT_Ph3_SSN_GFL.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2026 Institute for Automation of Complex Power Systems, EONERC, RWTH Aachen University
// SPDX-License-Identifier: MPL-2.0
 
#pragma once
 
#include <vector>
 
#include <dpsim-models/EMT/EMT_Ph3_TwoTerminalVTypeVariableSSNComp.h>
 
namespace CPS {
namespace EMT {
namespace Ph3 {

class SSN_GFL final : public TwoTerminalVTypeVariableSSNComp,
                      public SharedFactory<SSN_GFL> {
private:
  static constexpr Int mStateSize = 14;
  static constexpr Int mInputSize = 3;
  static constexpr Int mOutputSize = 3;
 
  enum StateIndex : Int {
    ThetaPLL = 0,
    PhiPLL = 1,
    PFiltered = 2,
    QFiltered = 3,
    PhiD = 4,
    PhiQ = 5,
    GammaD = 6,
    GammaQ = 7,
    VcA = 8,
    VcB = 9,
    VcC = 10,
    IfA = 11,
    IfB = 12,
    IfC = 13
  };
 
  // Electrical parameters
  Real mLf;
  Real mCf;
  Real mRf;
  Real mRc;
 
  // PLL parameters
  Real mOmegaN;
  Real mKpPLL;
  Real mKiPLL;
 
  // Power measurement and controller parameters
  Real mOmegaCutoff;
  Real mPRef;
  Real mQRef;
  Real mKpPowerCtrl;
  Real mKiPowerCtrl;
  Real mKpCurrCtrl;
  Real mKiCurrCtrl;
 
  // Logging attributes
  const Attribute<Real>::Ptr mVcD;
  const Attribute<Real>::Ptr mVcQ;
  const Attribute<Real>::Ptr mIrcD;
  const Attribute<Real>::Ptr mIrcQ;
  const Attribute<Real>::Ptr mPInst;
  const Attribute<Real>::Ptr mQInst;
  const Attribute<Real>::Ptr mOmegaPLL;
 
  Matrix getParkTransformMatrix(Real theta) const;
  Matrix getInverseParkTransformMatrix(Real theta) const;

  void evaluateStateDerivative(const Matrix &x, const Matrix &u,
                               Matrix &stateDerivative) const;

  void evaluateOutput(const Matrix &x, const Matrix &u, Matrix &output) const;

  void calculateAnalyticalJacobians(const Matrix &x, const Matrix &u, Matrix &A,
                                    Matrix &B, Matrix &C, Matrix &D) const;

  void buildStateSpaceModel(const Matrix &x, const Matrix &u, Matrix &A,
                            Matrix &B, Matrix &C, Matrix &D, Matrix &E,
                            Matrix &F) const;
 
protected:
  Bool updateComponentParameters() override final;
  void updateLogAttributes(const Matrix &u) const override final;
 
public:
  using SharedFactory<SSN_GFL>::make;
 
  SSN_GFL(String uid, String name, Logger::Level logLevel = Logger::Level::off);
 
  SSN_GFL(String name, Logger::Level logLevel = Logger::Level::off)
      : SSN_GFL(name, name, logLevel) {}
 
  std::vector<String> getLocalStateNames() const override final;
 
  std::vector<SSNComp::LocalAbcStateBlock>
  getLocalAbcStateBlocks() const override final;

  void setParameters(Real lf, Real cf, Real rf, Real rc, Real omegaN,
                     Real kpPLL, Real kiPLL, Real omegaCutoff, Real pRef,
                     Real qRef, Real kpPowerCtrl, Real kiPowerCtrl,
                     Real kpCurrCtrl, Real kiCurrCtrl);
 
  void initializeFromNodesAndTerminals(Real frequency) override final;
 
  Matrix getState() const;
  Matrix getStateDerivative() const;
  Matrix getInterfaceVoltage() const;
  Matrix getInterfaceCurrent() const;
};
 
} // namespace Ph3
} // namespace EMT
} // namespace CPS