SP_Ph1_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_Ph1_Transformer.h>
#include <dpsim-models/CompositePowerComp.h>
#include <dpsim-models/SP/SP_Ph1_Capacitor.h>
#include <dpsim-models/SP/SP_Ph1_Inductor.h>
#include <dpsim-models/SP/SP_Ph1_Resistor.h>
#include <dpsim-models/Solver/MNAInterface.h>
#include <dpsim-models/Solver/PFSolverInterfaceBranch.h>
 
namespace CPS {
namespace SP {
namespace Ph1 {
class Transformer : public CompositePowerComp<Complex>,
                    public Base::Ph1::Transformer,
                    public SharedFactory<Transformer>,
                    public PFSolverInterfaceBranch {
 
private:
  std::shared_ptr<SP::Ph1::Resistor> mSubResistor;
  std::shared_ptr<SP::Ph1::Inductor> mSubInductor;

  std::shared_ptr<SP::Ph1::Resistor> mSubResistor2;
  std::shared_ptr<SP::Ph1::Inductor> mSubInductor2;
  std::shared_ptr<SP::Ph1::Resistor> mSubMagnetizingResistor;
  std::shared_ptr<SP::Ph1::Inductor> mSubMagnetizingInductor;

  Real mRatioAbs = 1;
  Real mRatioPhase = 0;
  Real mNominalOmega;
  Real mReactance;

  Complex mLeakage;

  Real mBaseApparentPower;
  Real mBaseImpedance;
  Real mBaseInductance;
  Real mBaseAdmittance;
  Real mBaseOmega;

  Real mBaseCurrent;

  Real mResistancePerUnit;
  Real mReactancePerUnit;
  Real mInductancePerUnit;
  Complex mLeakagePerUnit;
  Complex mMagnetizingPerUnit;
  Real mRatioAbsPerUnit;
  Complex mRatioPerUnit;
 
  // #### Admittance matrix stamp ####
  MatrixComp mY_element;

  Real mMagnetizingResistance;
  Real 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;
  }

  const Attribute<Real>::Ptr mBaseVoltage;
 
  // #### Power flow results ####
  const Attribute<MatrixComp>::Ptr mCurrent;

  const Attribute<Matrix>::Ptr mActivePowerBranch;

  const Attribute<Matrix>::Ptr mReactivePowerBranch;
  const Attribute<Real>::Ptr mActivePowerInjection;
  const Attribute<Real>::Ptr mReactivePowerInjection;

  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<Complex>::Ptr clone(String name) override;
 
  // #### General ####
  void setParameters(Real nomVoltagePrimary, Real nomVoltageSecondary,
                     Real ratioAbs, Real ratioPhase, Real resistance,
                     Real inductance);
  void setParameters(Real nomVoltagePrimary, Real nomVoltageSecondary,
                     Real ratedPower, Real ratioAbs, Real ratioPhase,
                     Real resistance, Real inductance);
  void createSubComponents() override;
  void initializeParentFromNodesAndTerminals(Real frequency) override;
 
  // #### Powerflow section ####
  Real getNominalVoltagePrimary() const;
  Real getNominalVoltageSecondary() const;
  void setBaseVoltage(Real baseVoltage);
  void calculatePerUnitParameters(Real baseApparentPower, Real baseOmega);
  void pfApplyAdmittanceMatrixStamp(SparseMatrixCompRow &Y) override;
  void updateBranchFlow(VectorComp &current, VectorComp &powerflow);
  void storeNodalInjection(Complex powerInjection);
 
  // #### Getter ####
  MatrixComp Y_element();
 
  // #### 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 Ph1
} // namespace SP
} // namespace CPS