Commit 045a6ffc authored by Viviane Sapucaia's avatar Viviane Sapucaia
Browse files

corrected mistake in equation of ifd

parent 2d78a293
......@@ -271,7 +271,7 @@ void SynchronGeneratorEMT::stepInPerUnit(Real om, Real dt, Real fieldVoltage, Re
mId_past = mId;
mIq_past = mIq;
mId = ((mLlfd*mLlkd + mLmd*(mLlfd + mLlkd))*mPsid - mLmd*mLlkd*mPsifd - mLlfd*mLmd*mPsikd) / detLd;
mIfd = (mLlkd*mLmd*mPsid - (mLl*mLlkd + mLmd*(mLl + mLlkd))*mPsifd + mLmd*mLl) / detLd;
mIfd = (mLlkd*mLmd*mPsid - (mLl*mLlkd + mLmd*(mLl + mLlkd))*mPsifd + mLmd*mLl*mPsikd) / detLd;
mIkd = (mLmd*mLlfd*mPsid + mLmd*mLl*mPsifd - (mLmd*(mLlfd + mLl) + mLl*mLlfd)*mPsikd) / detLd;
mIq = ((mLlkq1*mLlkq2 + mLmq*(mLlkq1 + mLlkq2))*mPsiq - mLmq*mLlkq2*mPsikq1 - mLmq*mLlkq1*mPsikq2) / detLq;
mIkq1 = (mLmq*mLlkq2*mPsiq - (mLmq*(mLlkq2 + mLl) + mLl*mLlkq2)*mPsikq1 + mLmq*mLl*mPsikq2) / detLq;
......
......@@ -42,8 +42,8 @@ int main(int argc, char* argv[]) {
//SynGenUnitTestBalancedResLoad();
//SynGenUnitTestPhaseToPhaseFault();
//SynGenUnitTestThreePhaseFault();
SynGenDPUnitTestBalancedResLoad();
SynGenUnitTestThreePhaseFault();
//SynGenDPUnitTestBalancedResLoad();
//RXLineResLoad();
//VarFreqRXLineResLoad();
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment