DecouplingLineEMT_Ph3.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/Definitions.h"
#include "dpsim-models/EMT/EMT_Ph3_ControlledCurrentSource.h"
#include <vector>
 
#include <dpsim-models/EMT/EMT_Ph3_Resistor.h>
#include <dpsim-models/SimSignalComp.h>
#include <dpsim-models/Task.h>
 
namespace CPS {
namespace Signal {
class DecouplingLineEMT_Ph3 : public SimSignalComp,
                              public SharedFactory<DecouplingLineEMT_Ph3> {
protected:
  Real mDelay;
  Matrix mResistance = Matrix::Zero(3, 3);
  Matrix mInductance = Matrix::Zero(3, 3);
  Matrix mCapacitance = Matrix::Zero(3, 3);
  Matrix mSurgeImpedance;
 
  std::shared_ptr<EMT::SimNode> mNode1, mNode2;
  std::shared_ptr<EMT::Ph3::Resistor> mRes1, mRes2;
  std::shared_ptr<EMT::Ph3::ControlledCurrentSource> mSrc1, mSrc2;
  Attribute<Matrix>::Ptr mSrcCur1, mSrcCur2;
 
  // Ringbuffers for the values of previous timesteps
  // TODO make these matrix attributes
  Matrix mVolt1, mVolt2, mCur1, mCur2;
  UInt mBufIdx = 0;
  UInt mBufSize;
  Real mAlpha;
 
  Matrix interpolate(Matrix &data);
 
public:
  typedef std::shared_ptr<DecouplingLineEMT_Ph3> Ptr;
 
  const Attribute<Matrix>::Ptr mSrcCur1Ref;
  const Attribute<Matrix>::Ptr mSrcCur2Ref;

  const Attribute<Matrix>::Ptr mStates;
 
  DecouplingLineEMT_Ph3(String name,
                        Logger::Level logLevel = Logger::Level::info);
 
  void setParameters(SimNode<Real>::Ptr node1, SimNode<Real>::Ptr node2,
                     Matrix resistance, Matrix inductance, Matrix capacitance);
  void initialize(Real omega, Real timeStep);
  void step(Real time, Int timeStepCount);
  void postStep();
  Task::List getTasks();
  IdentifiedObject::List getLineComponents();
 
  class PreStep : public Task {
  public:
    PreStep(DecouplingLineEMT_Ph3 &line)
        : Task(**line.mName + ".MnaPreStep"), mLine(line) {
      mPrevStepDependencies.push_back(mLine.mStates);
      mModifiedAttributes.push_back(mLine.mSrc1->mCurrentRef);
      mModifiedAttributes.push_back(mLine.mSrc2->mCurrentRef);
    }
 
    void execute(Real time, Int timeStepCount);
 
  private:
    DecouplingLineEMT_Ph3 &mLine;
  };
 
  class PostStep : public Task {
  public:
    PostStep(DecouplingLineEMT_Ph3 &line)
        : Task(**line.mName + ".PostStep"), mLine(line) {
      mAttributeDependencies.push_back(mLine.mRes1->mIntfVoltage);
      mAttributeDependencies.push_back(mLine.mRes1->mIntfCurrent);
      mAttributeDependencies.push_back(mLine.mRes2->mIntfVoltage);
      mAttributeDependencies.push_back(mLine.mRes2->mIntfCurrent);
      mModifiedAttributes.push_back(mLine.mStates);
    }
 
    void execute(Real time, Int timeStepCount);
 
  private:
    DecouplingLineEMT_Ph3 &mLine;
  };
};
} // namespace Signal
} // namespace CPS