EMT_Ph3_Transformer.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_Ph3_Transformer.h>
#include <dpsim-models/CompositePowerComp.h>
#include <dpsim-models/EMT/EMT_Ph3_Capacitor.h>
#include <dpsim-models/EMT/EMT_Ph3_Inductor.h>
#include <dpsim-models/EMT/EMT_Ph3_Resistor.h>
#include <dpsim-models/Solver/MNAInterface.h>
 
namespace CPS {
namespace EMT {
namespace Ph3 {
class Transformer : public CompositePowerComp<Real>,
                    public SharedFactory<Transformer>,
                    public Base::Ph3::Transformer {
private:
  std::shared_ptr<EMT::Ph3::Resistor> mSubResistor;
  std::shared_ptr<EMT::Ph3::Inductor> mSubInductor;

  std::shared_ptr<EMT::Ph3::Resistor> mSubResistor2;
  std::shared_ptr<EMT::Ph3::Inductor> mSubInductor2;
  std::shared_ptr<EMT::Ph3::Resistor> mSubMagnetizingResistor;
  std::shared_ptr<EMT::Ph3::Inductor> mSubMagnetizingInductor;

  Matrix mMagnetizingResistance;
  Matrix mMagnetizingInductance;
  Real mNoLoadCurrent = 0.01;
  Real mNoLoadLoss = 1e-3;

  Bool mWithResistiveLosses;

  UInt mReferenceTerminal = 0;
  UInt nonReferenceTerminal() const { return 1 - mReferenceTerminal; }
  Complex mRatioFromReference;
  Real mOrientationSign = 1.;
  WindingReference mReferenceWinding = WindingReference::Auto;
  void resolveWindingRoles();
 
public:
  void setMagnetizingBranch(Real noLoadCurrent, Real noLoadLoss) {
    mNoLoadCurrent = noLoadCurrent;
    mNoLoadLoss = noLoadLoss;
  }

  void setReferenceWinding(WindingReference referenceWinding) {
    mReferenceWinding = referenceWinding;
  }

  Transformer(String uid, String name,
              Logger::Level logLevel = Logger::Level::off,
              Bool withResistiveLosses = false);
  Transformer(String name, Logger::Level logLevel = Logger::Level::off)
      : Transformer(name, name, logLevel) {}
 
  SimPowerComp<Real>::Ptr clone(String name) override;
 
  // #### General ####
  void setParameters(Real nomVoltagePrimary, Real nomVoltageSecondary,
                     Real ratedPower, Real ratioAbs, Real ratioPhase,
                     Matrix resistance, Matrix inductance);
  void createSubComponents() override;
  void initializeParentFromNodesAndTerminals(Real frequency) override;
 
  // #### MNA section ####
  void mnaParentInitialize(Real omega, Real timeStep,
                           Attribute<Matrix>::Ptr leftVector) override;
  void mnaCompApplySystemMatrixStamp(SparseMatrixRow &systemMatrix) override;
  void mnaCompUpdateCurrent(const Matrix &leftVector) override;
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaParentPreStep(Real time, Int timeStepCount) override;
  void mnaParentPostStep(Real time, Int timeStepCount,
                         Attribute<Matrix>::Ptr &leftVector) override;
  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;
};
} // namespace Ph3
} // namespace EMT
} // namespace CPS