DP_Ph3_SynchronGeneratorVBR_Deprecated.h

Go to the documentation of this file.

/* 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_Exciter.h>
#include <dpsim-models/Base/Base_SynchronGenerator.h>
#include <dpsim-models/Signal/TurbineGovernor.h>
#include <dpsim-models/SimPowerComp.h>
#include <dpsim-models/Solver/MNAInterface.h>
 
namespace CPS {
namespace DP {
namespace Ph3 {
class SynchronGeneratorVBR : public SimPowerComp<Complex>,
                             public Base::SynchronGenerator,
                             public MNAInterface,
                             public SharedFactory<SynchronGeneratorVBR> {
 
protected:
  Real mDLmd;
  Real mDLmq;

  Real mLa;
  Real mLb;
 
  //Real mThetaMech;
  Real mThetaMech2;
  Real mTheta0;

  Real mDPsid;
  Real mDPsiq;

  Real mVaRe;
  Real mVaIm;
  Real mVbRe;
  Real mVbIm;
  Real mVcRe;
  Real mVcIm;

  Real mIaRe;
  Real mIaIm;
  Real mIbRe;
  Real mIbIm;
  Real mIcRe;
  Real mIcIm;

  Real mPsimq;
  Real mPsimd;

  Matrix mDqStatorCurrents = Matrix::Zero(2, 1);
 
  // ### Useful Matrices ###
  Matrix mDInductanceMat = Matrix::Zero(3, 3);
  Matrix R_load = Matrix::Zero(6, 6);
  Matrix LD0 = Matrix::Zero(3, 3);
  Matrix L_EQ = Matrix::Zero(6, 6);
  Matrix R_EQ = Matrix::Zero(6, 6);
  Matrix R_eq_DP = Matrix::Zero(6, 6);
  Matrix E_eq_DP = Matrix::Zero(6, 1);

  Matrix mIabc = Matrix::Zero(6, 1);
  Matrix mDVabc = Matrix::Zero(6, 1);
  Matrix mPsikq1kq2 = Matrix::Zero(2, 1);
  Matrix mPsifdkd = Matrix::Zero(2, 1);
  Real mIq_hist;
  Real mId_hist;
  Matrix mKrs_teta = Matrix::Zero(3, 3);
  Matrix mKrs_teta_inv = Matrix::Zero(3, 3);

  Matrix mConductanceMat = Matrix::Zero(6, 6);
  Matrix mISourceEq = Matrix::Zero(6, 1);

  Real c21_omega;
  Real c22_omega;
  Real c13_omega;
  Real c14_omega;
  Matrix K1a = Matrix::Zero(2, 2);
  Matrix K1b = Matrix::Zero(2, 1);
  Matrix K1 = Matrix::Zero(2, 1);
  Matrix K2a = Matrix::Zero(2, 2);
  Matrix K2b = Matrix::Zero(2, 1);
  Matrix K2 = Matrix::Zero(2, 1);
  Matrix h_qdr;
  Matrix K = Matrix::Zero(3, 3);
  Matrix E_r_vbr = Matrix::Zero(3, 1);
  Matrix K_DP = Matrix::Zero(6, 6);

  Real c11;
  Real c12;
  Real c23;
  Real c24;
  Real c15;
  Real c25;
  Real c26;
  Real b11;
  Real b12;
  Real b13;
  Real b21;
  Real b22;
  Real b23;
  Real b31;
  Real b32;
  Real b33;
  Real b41;
  Real b42;
  Real b43;
 
  Real E1_1d;
  Real E2_1d;
 
  Matrix Ea = Matrix::Zero(2, 2);
  Matrix E1b = Matrix::Zero(2, 1);
  Matrix E1 = Matrix::Zero(2, 2);
  Matrix Fa = Matrix::Zero(2, 2);
  Matrix F1b = Matrix::Zero(2, 1);
  Matrix F1 = Matrix::Zero(2, 2);
  Matrix E2b = Matrix::Zero(2, 2);
  Matrix E2 = Matrix::Zero(2, 2);
  Matrix F2b = Matrix::Zero(2, 2);
  Matrix F2 = Matrix::Zero(2, 2);
  Matrix F3b = Matrix::Zero(2, 1);
  Matrix F3 = Matrix::Zero(2, 2);
  Matrix C26 = Matrix::Zero(2, 1);
  Matrix A = Matrix::Zero(6, 6);
  Matrix B = Matrix::Zero(6, 6);
  Matrix Var1 = Matrix::Zero(6, 6);
  Matrix Var2 = Matrix::Zero(6, 6);
  Matrix mVabc = Matrix::Zero(6, 1);
  Matrix E_r_vbr_DP = Matrix::Zero(6, 1);
  Matrix E_r_vbr_DP2 = Matrix::Zero(6, 1);
 
public:
  SynchronGeneratorVBR(String name, Real nomPower, Real nomVolt, Real nomFreq,
                       Int poleNumber, Real nomFieldCur, Real Rs, Real Ll,
                       Real Lmd, Real Lmd0, Real Lmq, Real Lmq0, Real Rfd,
                       Real Llfd, Real Rkd, Real Llkd, Real Rkq1, Real Llkq1,
                       Real Rkq2, Real Llkq2, Real inertia,
                       Logger::Level logLevel = Logger::Level::off);

  void addExciter(Real Ta, Real Ka, Real Te, Real Ke, Real Tf, Real Kf, Real Tr,
                  Real Lad, Real Rfd);
  void addGovernor(Real Ta, Real Tb, Real Tc, Real Fa, Real Fb, Real Fc, Real K,
                   Real Tsr, Real Tsm, Real Tm_init, Real PmRef);
 
  void initialize(Matrix frequencies) override {
    SimPowerComp<Complex>::initialize(frequencies);
  }

  void initialize(Real om, Real dt, Real initActivePower,
                  Real initReactivePower, Real initTerminalVolt,
                  Real initVoltAngle, Real initMechPower);

  void mnaStep(Matrix &systemMatrix, Matrix &rightVector, Matrix &leftVector,
               Real time);

  void stepInPerUnit(Real om, Real dt, Real time, NumericalMethod numMethod);

  void mnaCompPostStep(Matrix &rightVector, Matrix &leftVector, Real time);

  Matrix abcToDq0Transform(Real theta, Real aRe, Real bRe, Real cRe, Real aIm,
                           Real bIm, Real cIm);

  Matrix dq0ToAbcTransform(Real theta, Real d, Real q, Real zero);
 
  void CalculateLandR(Real time, Real dt);
  void CalculateAuxiliarConstants(Real dt);
  void CalculateAuxiliarVariables(Real time);
 
  Matrix &dqStatorCurrents() { return mDqStatorCurrents; }
  Real electricalTorque() { return **mElecTorque * mBase_T; }
  Real rotationalSpeed() { return **mOmMech * mBase_OmMech; }
  Real rotorPosition() { return mThetaMech; }
  Matrix &statorCurrents() { return mIabc; }
};
} // namespace Ph3
} // namespace DP
} // namespace CPS