Research Article

Cascading Failures Analysis Considering Extreme Virus Propagation of Cyber-Physical Systems in Smart Grids

Algorithm 1

Dynamic power flow method.
Input: Electrical network information
Output: Power flows of branches, frequency violation
Step 1: Unbalanced power: Employ Equation (13) to
   calculate the unbalanced power of the
   network.
Step 2: System frequency characteristic: IF , employ
   Equation (14) to calculate and go to the next
   step; ELSE, IF , employ Equation (15) to
   calculate and go to the next step; ELSE, go to Step 4.
Step 3: Frequency calculation: Calculate variation of
   frequency of the network by Equation (12) solved
   using improved Euler method (Δt = 0.01 s).
Step 4: Generator frequency detection: IF is beyond
   limits of generator i, cut off generator i and return
   to Step 3; ELSE, go to the next step.
Step 5: Frequency violation: IF is beyond limits
   the network, end the algorithm and send warnings of
   of instability; ELSE, go to the next step.
Step 6: Power flow calculation: Employ the P-Q power
   flow to calculate the power flows of every branch.