Example: ALTERNATOR05
Simple Alternator Stability Model

Description

Simple Model of Alternator and main Network Stability.
Occurence of Fault on Coupling overhead Line and after that Clearing of it.
Transient Action solved from 0 to 2 sec.
Chosen clearing Time is 0.7 sec. In this case Stability was restored after slip.

System Parameters

Pt0166Pt0 = 166[MW]Turbine output Power [MW]
Pmax1125.8Pmax1 = 125.8[MW]Maximal electric Power given by Alternator at usual State
Pmax24.8Pmax2 = 4.8[MW]Maximal electric Power given by Alternator when Fault occurs
Pmax3121.1Pmax3 = 121.1[MW]Maximal electric Power after Fault Clearing
Tporuch0.2Tporuch = 0.2[s]Time of Fault Occurence
Tvyppor0.7Tvyppor = 0.7[s]Time of Fault Clearing
Tf0.35Tf = 0.35[s]Time Constant of Rotor exciting Winding
Omega02PI*50ω0 = 2.π.50[rad/s]Angular synchronous Velocity
Tm1.6Tm = 1.6[s]Alternator mechanical Time Constant
Sng220Sng = 220[MVA]Nominal Alternator Power

Model

[picture]
[end of picture]

Data

Open in DYNAST Shell
  Analyzed system
  Analysis results
*: Simple Alternator Stability Model

*SYSTEM;
Pt0=166;   :: [MW] Turbine output Power [MW]


Pmax1=125.8;   :: [MW] Maximal electric Power given by Alternator at usual State
Pmax2=4.8;     :: [MW] Maximal electric Power given by Alternator when Fault
occurs
Pmax3=121.1;   :: [MW] Maximal electric Power after Fault Clearing
Pa1 /SIN/ B=Pmax1;
Pa2 /SIN/ B=Pmax2;
Pa3 /SIN/ B=Pmax3;

Tporuch=0.2;    :: [s] Time of Fault Occurence
Tvyppor=0.7;   :: [s] Time of Fault Clearing
Tf=0.35;        :: [s] Time Constant of Rotor exciting Winding

: Exciting current with Respect to Impact of overexciting System
Ib=1+1.5*(1-EXP(-(TIME-Tporuch)/Tf))*(TIME>Tporuch); :: Exciting current

: Final mechanical Power without Losses
DeltaP=-Pt0/3+Ib*(
       Pa1(Theta)*(TIME<Tporuch)
      +Pa2(Theta)*(TIME>=Tporuch)*(TIME<Tvyppor)
      +Pa3(Theta)*(TIME>=Tvyppor));

Omega0=2PI*50;   :: [rad/s] Angular synchronous Velocity
Tm=1.6;          :: [s] Alternator mechanical Time Constant
Sng=220;         :: [MVA] Nominal Alternator Power


DeltaM > J Omega = DeltaP/Omega0; :: [N*m]Moment Difference

Jm=Tm*Sng/(Omega0*Omega0); :: [kg*m^2] Inertia Moment
InertiaMoment > C Omega = Jm;

PowerAngle > @Int Omega,Theta; :: [rad] Power Angle

*TR;
init Theta=0.457, Omega=0;
TR 0 2;
PRINT(1001) DeltaP, Omega, Theta, Ib;
RUN;
*END;
:: DeltaP [MW] Final mechanical Power without Losses
:: Omega [rad/s] Angular Velocity
:: Theta [rad] Power Angle
:: Ib [rad/s] Exciting current with Respect to Impact of overexciting System

The following data were generated by DYNAST for the above given system model. You can modify the model parameters and re-submit the data to DYNAST to receive new results.

  

Submodels

[picture]
[end of picture]
[picture]
[end of picture]
[picture]
[end of picture]
[picture]
[end of picture]
[picture]
[end of picture]

Origin

Karel Noháč KEE, FEL, ZČU v Plzni

Last Update

November 13, 2014