First, for future reference, please use a subject line related to your question 
rather than simply replying to an unrelated thread.

You can simply convert your loads from fixed to dispatchable using 
load2disp()<https://matpower.org/docs/ref/matpower7.1/lib/load2disp.html> and 
run an OPF.

For example with case9, with changes all of the load is served with no problems 
 …

mpc = load2disp(loadcase('case9'));
runopf(mpc)

================================================================================
|     Bus Data                                                                 |
================================================================================
 Bus      Voltage          Generation             Load          Lambda($/MVA-hr)
  #   Mag(pu) Ang(deg)   P (MW)   Q (MVAr)   P (MW)   Q (MVAr)     P        Q
----- ------- --------  --------  --------  --------  --------  -------  -------
    1  1.100    0.000*    89.80     12.86       -         -      24.756     -
    2  1.097    4.892    134.32      0.09       -         -      24.035     -
    3  1.087    3.247     94.19    -22.58       -         -      24.076     -
    4  1.094   -2.464       -         -         -         -      24.756   0.005
    5  1.084   -3.983       -         -       90.00*    30.00*   24.999   0.027
    6  1.100    0.601       -         -         -         -      24.076     -
    7  1.089   -1.198       -         -      100.00*    35.00*   24.254   0.036
    8  1.100    0.904       -         -         -         -      24.035     -
    9  1.072   -4.617       -         -      125.00*    50.00*   24.999   0.112
                        --------  --------  --------  --------
               Total:    318.31     -9.62    315.00    115.00



… but with the line ratings on lines feeding bus 7 reduced, it is forced to 
curtail the bus 7 load …

define_constants;
mpc.branch(5:6, RATE_A) = 40;
runopf(mpc)

================================================================================
|     Bus Data                                                                 |
================================================================================
 Bus      Voltage          Generation             Load          Lambda($/MVA-hr)
  #   Mag(pu) Ang(deg)   P (MW)   Q (MVAr)   P (MW)   Q (MVAr)     P        Q
----- ------- --------  --------  --------  --------  --------  -------  -------
    1  1.100    0.000*    83.52     13.79       -         -      23.375     -
    2  1.094    4.353    110.02     -7.40       -         -      19.903     -
    3  1.083    4.440    100.28    -28.15       -         -      25.570     -
    4  1.094   -2.292       -         -         -         -      23.375   0.008
    5  1.083   -3.523       -         -       90.00*    30.00*   24.624   0.014
    6  1.099    1.609       -         -         -         -      25.570     -
    7  1.094   -0.200       -         -       75.51*    26.43* 4456.8491551.860
    8  1.100    1.077       -         -         -         -      19.903     -
    9  1.071   -4.445       -         -      125.00*    50.00*   22.624   0.091
                        --------  --------  --------  --------
               Total:    293.83    -21.76    290.51    106.43


Hope this helps,

   Ray


On Nov 23, 2020, at 3:33 PM, [email protected]<mailto:[email protected]> 
wrote:
From: [email protected]<mailto:[email protected]>
Subject: Re: Bus, Loads and Generators question
Date: November 23, 2020 at 3:33:40 PM CST


Hi Ray How do you compute optimal amount of load shedding in MATPOWER



Reply via email to