Power System Stability and Control
Table of Contents

Synchronous machines are the backbone of power systems, and understanding their dynamics is crucial. These equations describe how the machine's electrical and mechanical components interact, allowing us to predict and control its behavior during normal operation and disturbances.

The dq0 reference frame simplifies the analysis by transforming three-phase AC quantities into two DC components. This transformation helps us model the machine's response to changes in load, voltage, and frequency, making it easier to design control systems and ensure grid stability.

Synchronous Machine Dynamics in dq0 Frame

Dynamic Equations Formulation

  • Transform the three-phase stator quantities into two DC quantities (d and q axes) and a zero-sequence component using the dq0 reference frame, a rotating reference frame that simplifies the analysis of synchronous machines
  • Align the d-axis with the rotor field winding, while the q-axis is 90 degrees ahead of the d-axis in the direction of rotation
  • Express the stator voltage equations in the dq0 reference frame as:
    • $v_d = -R_si_d - \omega\lambda_q + d\lambda_d/dt$
    • $v_q = -R_si_q + \omega\lambda_d + d\lambda_q/dt$
    • $v_0 = -R_si_0 + d\lambda_0/dt$
  • Represent the stator flux linkage equations in the dq0 reference frame as:
    • $\lambda_d = -L_di_d + L_{md}i_{fd}$
    • $\lambda_q = -L_qi_q$
    • $\lambda_0 = -L_0i_0$

Rotor Equations in dq0 Frame

  • Write the rotor voltage equations in the dq0 reference frame as:
    • $v_{fd} = R_fi_{fd} + d\lambda_{fd}/dt$
    • $v_{kd} = R_{kd}i_{kd} + d\lambda_{kd}/dt$
    • $v_{kq} = R_{kq}i_{kq} + d\lambda_{kq}/dt$
  • Express the rotor flux linkage equations in the dq0 reference frame as:
    • $\lambda_{fd} = L_{fd}i_{fd} + L_{md}i_d$
    • $\lambda_{kd} = L_{kd}i_{kd} + L_{md}i_d$
    • $\lambda_{kq} = L_{kq}i_{kq} + L_{mq}i_q$
  • Use these equations to analyze the dynamic behavior of the synchronous machine under various operating conditions (steady-state, transient, and subtransient)

Stator-Rotor Equation Relationship

Coupling through Mutual Inductances

  • Couple the stator and rotor equations through the mutual inductances $L_{md}$ and $L_{mq}$, which represent the magnetic coupling between the stator and rotor windings
  • Determine the stator d-axis flux linkage ($\lambda_d$) by both the stator d-axis current ($i_d$) and the rotor field current ($i_{fd}$), while the stator q-axis flux linkage ($\lambda_q$) depends only on the stator q-axis current ($i_q$)
  • Express the rotor field flux linkage ($\lambda_{fd}$) and the d-axis damper winding flux linkage ($\lambda_{kd}$) as functions of both the rotor field current ($i_{fd}$) and the stator d-axis current ($i_d$)
  • Relate the rotor q-axis damper winding flux linkage ($\lambda_{kq}$) to both the rotor q-axis damper winding current ($i_{kq}$) and the stator q-axis current ($i_q$)

Dynamic Behavior Analysis

  • Analyze the dynamic behavior of the synchronous machine under various operating conditions using the coupling between the stator and rotor equations
  • Investigate the machine's response to sudden changes in the system (faults, load variations) by examining the relationship between stator and rotor quantities
  • Design appropriate control strategies (excitation systems, power system stabilizers) to enhance the machine's performance and stability based on the coupled stator-rotor equations
  • Optimize the machine's operation (efficiency, power factor) by considering the interaction between stator and rotor quantities

Saliency Impact on Synchronous Machines

Saliency and Inductance Difference

  • Define saliency as the difference between the d-axis and q-axis inductances ($L_d$ and $L_q$) in a synchronous machine
  • Observe that in a salient-pole machine, the d-axis inductance is smaller than the q-axis inductance ($L_d < L_q$) due to the non-uniform air gap caused by the salient poles
  • Note that in non-salient-pole machines (cylindrical rotor machines), the d-axis and q-axis inductances are equal ($L_d = L_q$)

Reluctance Torque

  • Recognize that the presence of saliency introduces an additional torque component called the reluctance torque
  • Calculate the reluctance torque as proportional to the difference between the d-axis and q-axis inductances ($L_d - L_q$) and the product of the d-axis and q-axis currents ($i_di_q$)
  • Observe that the reluctance torque can enhance the overall torque production in salient-pole machines, especially at high load angles
  • Note that in non-salient-pole machines, the reluctance torque is zero due to equal d-axis and q-axis inductances

Impact on Machine Characteristics

  • Understand that saliency affects the machine's transient and subtransient reactances, which determine the machine's response to sudden changes in the system (faults, load variations)
  • Analyze the effect of saliency on the machine's power-angle characteristics, as the reluctance torque contributes to the overall torque production
  • Consider the impact of saliency on the machine's excitation requirements, as the difference in d-axis and q-axis inductances influences the reactive power generation
  • Investigate the role of saliency in the machine's stability, as the reluctance torque can provide additional damping during power system oscillations

Dynamic Response Analysis of Synchronous Machines

Numerical Solution of Dynamic Equations

  • Solve the dynamic equations of a synchronous machine using numerical integration methods (Runge-Kutta, Euler) to determine the machine's response to various disturbances
  • Specify the initial conditions for the machine's state variables (currents, flux linkages) based on the steady-state operating point before the disturbance
  • Simulate disturbances such as changes in the mechanical input power, electrical load, or system voltage and frequency
  • Obtain the time-domain waveforms of the machine's state variables (stator and rotor currents, flux linkages, rotor speed) from the solution of the dynamic equations

Response Analysis and Control Design

  • Analyze the machine's response to determine its stability, damping, and oscillatory behavior under different operating conditions and disturbances
  • Identify the dominant modes of oscillation and their associated damping ratios and frequencies from the time-domain waveforms or eigenvalue analysis
  • Design appropriate control strategies (excitation systems, power system stabilizers) to enhance the machine's performance and stability based on the dynamic response analysis
  • Tune the controller parameters (gains, time constants) to achieve the desired damping, response time, and robustness against system uncertainties and disturbances
  • Validate the control design through simulation studies and sensitivity analysis to ensure satisfactory performance under a wide range of operating conditions