Refactor VBR models

[dinkelbachjan]

• Reuse the order-specific auxiliar constants in the domain specific VBR models

  • Example 4th order:
    • EMT_Ph3_SynchronGenerator4OrderVBR
    • The same constants are also used by DP_Ph1_SynchronGenerator4OrderVBR and SP_Ph1_SynchronGenerator4OrderVBR
  • Instantiate and calculate the order-specific auxiliar constants already in Base_ReducedOrderSynchronGenerator
  • Same goes for 3rd and 6th order
  • Stick for 3rd/4th order model to nomenclature for auxiliar constants introduced in 6th order model, i.e. use _t as subscript e.g. mAd_t

• Specify and reuse the order-specific resistance matrix constants in the domain specific VBR models

  • Example 4th order
    • DP_Ph1_SynchronGenerator4OrderVBR
    • The same constants can also used by EMT_Ph3_SynchronGenerator4OrderVBR and SP_Ph1_SynchronGenerator4OrderVBR
  • Instantiate and calculate the order-specific resistance matrix constants already in Base_ReducedOrderSynchronGenerator
  • Same goes for 3rd and 6th order

• Remove Ef from auxiliar constants

  • Example 4th order:
    • EMT_Ph3_SynchronGenerator4OrderVBR
  • Remove Ef

• Rename Cq to Dq

  • Example 4th order:
    • EMT_Ph3_SynchronGenerator4OrderVBR

• Reconsider zero component in EMT_Ph3 models

• Use resistance matrix constants A and B in EMT_Ph3_SynchronGenerator4OrderVBR

• Simplify the DP-DQ interface to obtain armature voltage (and armature current), which currently looks as follows
  

  dpsim/models/Source/DP/DP_Ph1_SynchronGeneratorVBR.cpp

  Lines 97 to 102 in 2499325

  	// convert armature voltage into dq reference frame
  	Matrix parkTransform = get_parkTransformMatrix();
  	MatrixComp Vabc_ = mIntfVoltage(0, 0) * mShiftVector * Complex(cos(mNomOmega * mSimTime), sin(mNomOmega * mSimTime));
  	Matrix Vabc = Matrix(3,1);
  	Vabc << Vabc_(0,0).real(), Vabc_(1,0).real(), Vabc_(2,0).real();
  	mVdq = parkTransform * Vabc / mBase_V_RMS;

  
  Instead of applying DP-EMT and EMT-DQ transformation, directly apply DP-DQ transformation, as e.g. done for the voltage behind the reactance
  

  dpsim/models/Source/DP/DP_Ph1_SynchronGenerator3OrderVBR.cpp

  Line 93 in 2499325

  mEvbr = (mKvbr * mEh_vbr * mBase_V_RMS)(0,0);

  
  using
  

  dpsim/models/Source/DP/DP_Ph1_SynchronGeneratorVBR.cpp

  Lines 65 to 67 in 2499325

  	mKvbr = Matrix::Zero(1,2);
  	mKvbr(0,0) = Complex(cos(mThetaMech - mBase_OmMech * mSimTime), sin(mThetaMech - mBase_OmMech * mSimTime));
  	mKvbr(0,1) = -Complex(cos(mThetaMech - mBase_OmMech * mSimTime - PI/2.), sin(mThetaMech - mBase_OmMech * mSimTime - PI/2.));

• Harmonize implementation of DQ-DP and DQ-SP interface implementation
  

  dpsim/models/Source/SP/SP_Ph1_SynchronGeneratorVBR.cpp

  Lines 80 to 82 in 15657fe

  	dqToComplexA <<
  	cos(mThetaMech - mBase_OmMech * mSimTime), -sin(mThetaMech - mBase_OmMech * mSimTime),
  	sin(mThetaMech - mBase_OmMech * mSimTime), cos(mThetaMech - mBase_OmMech * mSimTime);

• Simplify definition of VBR resistance matrix in DP_Ph1 models according to equivalent derivation
  

  dpsim/models/Source/DP/DP_Ph1_SynchronGeneratorVBR.cpp

  Lines 47 to 50 in e6f3776

  	resistanceMatrix(0,0) = mR_const_1ph.real() + mKa_1ph.real() + mKb_1ph.real();
  	resistanceMatrix(0,1) = -mR_const_1ph.imag() - mKa_1ph.imag() + mKb_1ph.imag();
  	resistanceMatrix(1,0) = mR_const_1ph.imag() + mKa_1ph.imag() + mKb_1ph.imag();
  	resistanceMatrix(1,1) = mR_const_1ph.real() + mKa_1ph.real() - mKb_1ph.real();

• Rename DP_Ph1_SynchronGeneratorVBR and SP_Ph1_SynchronGeneratorVBR to be consistent with EMT_Ph3_ReducedOrderSynchronGeneratorVBR

• Make SGModel of type String. Then you can use '6a' and '6b' instead of 6 and 7, which can become confusing in particular if we have models 6c, 6d etc. in the future

• We can add another case distinction during calculation of constants in Base::ReducedOrderSynchronGenerator, since mAd_t and mBd_t are only used within the fourth-order model.
  

  dpsim/dpsim-models/src/Base/Base_ReducedOrderSynchronGenerator.cpp

  Lines 216 to 217 in 2cca71a

  	mAd_t = mTimeStep * Zd_t / (2 * mTq0_t + mTimeStep);
  	mBd_t = (2 * mTq0_t - mTimeStep) / (2 * mTq0_t + mTimeStep);