EMT_Ph3_Switch.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 <array>
 
#include <dpsim-models/Base/Base_Ph3_Switch.h>
#include <dpsim-models/Definitions.h>
#include <dpsim-models/Logger.h>
#include <dpsim-models/MNASimPowerComp.h>
#include <dpsim-models/Solver/MNASwitchInterface.h>
#include <dpsim-models/Solver/MNAVariableCompInterface.h>
 
namespace CPS {
namespace EMT {
namespace Ph3 {

class Switch : public MNASimPowerComp<Real>,
               public Base::Ph3::Switch,
               public SharedFactory<Switch>,
               public MNAVariableCompInterface,
               public MNASwitchInterface {
public:
  enum class SwitchingMode {
    Ideal = 0,
    CurrentZero = 1,
    ExponentialZCSEmulation = 2,
  };

  const Attribute<Bool>::Ptr mOpeningRequested;

  const Attribute<Bool>::Ptr mPoleClosedA;
  const Attribute<Bool>::Ptr mPoleClosedB;
  const Attribute<Bool>::Ptr mPoleClosedC;

  const Attribute<Real>::Ptr mZeroCrossingTimeA;
  const Attribute<Real>::Ptr mZeroCrossingTimeB;
  const Attribute<Real>::Ptr mZeroCrossingTimeC;

  const Attribute<Bool>::Ptr mExponentialTransitionActive;
  const Attribute<Bool>::Ptr mExponentialTransitionClosing;
  const Attribute<Real>::Ptr mExponentialProgress;
  const Attribute<Real>::Ptr mExponentialTransitionStartTime;
  const Attribute<Real>::Ptr mExponentialTransitionEndTime;

  const Attribute<Real>::Ptr mEffectiveResistanceA;
  const Attribute<Real>::Ptr mEffectiveResistanceB;
  const Attribute<Real>::Ptr mEffectiveResistanceC;
 
  Switch(String uid, String name, Logger::Level logLevel = Logger::Level::off);
 
  Switch(String name, Logger::Level logLevel = Logger::Level::off)
      : Switch(name, name, logLevel) {}
 
  SimPowerComp<Real>::Ptr clone(String name) override;

  void setSwitchingMode(SwitchingMode mode);
  SwitchingMode switchingMode() const { return mSwitchingMode; }

  void setZeroCrossingTolerance(Real tolerance);
  Real zeroCrossingTolerance() const { return mZeroCrossingTolerance; }

  void setExponentialSwitchingTime(Real switchingTime);
  Real exponentialSwitchingTime() const { return mExponentialSwitchingTime; }

  void openSwitch() override;

  void closeSwitch() override;
 
  // #### General ####
  void initializeFromNodesAndTerminals(Real frequency) override;
 
  // #### General MNA section ####
  void mnaCompInitialize(Real omega, Real timeStep,
                         Attribute<Matrix>::Ptr leftVector) override;
 
  void mnaCompApplySystemMatrixStamp(SparseMatrixRow &systemMatrix) override;
  void mnaCompApplyRightSideVectorStamp(Matrix &rightVector) override;
 
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaCompUpdateCurrent(const Matrix &leftVector) override;
 
  // #### MNA switch interface ####
  Bool mnaIsClosed() override;

  Bool supportsPrecomputedSystemMatrices() const override;
 
  void mnaCompApplySwitchSystemMatrixStamp(Bool closed,
                                           SparseMatrixRow &systemMatrix,
                                           Int freqIdx) override;
 
  void mnaCompPostStep(Real time, Int timeStepCount,
                       Attribute<Matrix>::Ptr &leftVector) override;
 
  void
  mnaCompAddPostStepDependencies(AttributeBase::List &prevStepDependencies,
                                 AttributeBase::List &attributeDependencies,
                                 AttributeBase::List &modifiedAttributes,
                                 Attribute<Matrix>::Ptr &leftVector) override;
 
  // #### Variable component interface ####
  Bool hasParameterChanged() override;
 
private:
  SwitchingMode mSwitchingMode = SwitchingMode::Ideal;
 
  Real mZeroCrossingTolerance = 1e-6;

  Real mExponentialSwitchingTime = 0.01;

  Real mTimeStep = 0.0;
 
  std::array<Bool, 3> mPoleClosedPrev{{false, false, false}};
  std::array<Real, 3> mEffectiveResistancePrev{{0.0, 0.0, 0.0}};
 
  Matrix mPreviousCurrent = Matrix::Zero(3, 1);
  Real mPreviousCurrentTime = 0.0;
  Bool mPreviousCurrentValid = false;
 
  Bool mResetZeroCrossingHistory = false;
 
  Bool mExponentialTransitionStarted = false;
 
  Bool poleClosed(UInt phase) const;
  void setPoleClosed(UInt phase, Bool closed);
  void setAllPoles(Bool closed);
 
  Bool allPolesClosed() const;
  Bool allPolesOpen() const;
 
  void resetZeroCrossingTimes();
 
  Real effectiveResistance(UInt phase) const;
  void setEffectiveResistance(UInt phase, Real resistance);
  void synchronizeEffectiveResistance(Bool closed);
 
  void resetExponentialTransition();
  void validateExponentialResistanceParameters() const;
  Real exponentialResistance(UInt phase, Real alpha) const;
  void updateExponentialTransition(Real time);
 
  MatrixFixedSize<3, 3> currentConductanceMatrix() const;
 
  Bool currentCrossedZero(Real previousCurrent, Real current) const;
 
  Real interpolateZeroCrossingTime(Real previousCurrent, Real current,
                                   Real previousTime, Real currentTime) const;
 
  void setZeroCrossingTime(UInt phase, Real time);
  void updateZeroCrossingState(Real time);
};
 
} // namespace Ph3
} // namespace EMT
} // namespace CPS