Home > matpower6.0 > case33bw.m

case33bw

PURPOSE ^

CASE33BW Power flow data for 33 bus distribution system from Baran & Wu

SYNOPSIS ^

function mpc = case33bw

DESCRIPTION ^

CASE33BW  Power flow data for 33 bus distribution system from Baran & Wu
    Please see CASEFORMAT for details on the case file format.

    Data from ...
       M. E. Baran and F. F. Wu, "Network reconfiguration in distribution
       systems for loss reduction and load balancing," in IEEE Transactions
       on Power Delivery, vol. 4, no. 2, pp. 1401-1407, Apr 1989.
       doi: 10.1109/61.25627
       URL: http://doi.org/10.1109/61.25627

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function mpc = case33bw
0002 %CASE33BW  Power flow data for 33 bus distribution system from Baran & Wu
0003 %    Please see CASEFORMAT for details on the case file format.
0004 %
0005 %    Data from ...
0006 %       M. E. Baran and F. F. Wu, "Network reconfiguration in distribution
0007 %       systems for loss reduction and load balancing," in IEEE Transactions
0008 %       on Power Delivery, vol. 4, no. 2, pp. 1401-1407, Apr 1989.
0009 %       doi: 10.1109/61.25627
0010 %       URL: http://doi.org/10.1109/61.25627
0011 
0012 %% MATPOWER Case Format : Version 2
0013 mpc.version = '2';
0014 
0015 %%-----  Power Flow Data  -----%%
0016 %% system MVA base
0017 mpc.baseMVA = 100;
0018 
0019 %% bus data
0020 %    bus_i    type    Pd    Qd    Gs    Bs    area    Vm    Va    baseKV    zone    Vmax    Vmin
0021 mpc.bus = [  %% (Pd and Qd are specified in kW & kVAr here, converted to MW & MVAr below)
0022     1    3    0    0    0    0    1    1    0    12.66    1    1    1;
0023     2    1    100    60    0    0    1    1    0    12.66    1    1.1    0.9;
0024     3    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0025     4    1    120    80    0    0    1    1    0    12.66    1    1.1    0.9;
0026     5    1    60    30    0    0    1    1    0    12.66    1    1.1    0.9;
0027     6    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0028     7    1    200    100    0    0    1    1    0    12.66    1    1.1    0.9;
0029     8    1    200    100    0    0    1    1    0    12.66    1    1.1    0.9;
0030     9    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0031     10    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0032     11    1    45    30    0    0    1    1    0    12.66    1    1.1    0.9;
0033     12    1    60    35    0    0    1    1    0    12.66    1    1.1    0.9;
0034     13    1    60    35    0    0    1    1    0    12.66    1    1.1    0.9;
0035     14    1    120    80    0    0    1    1    0    12.66    1    1.1    0.9;
0036     15    1    60    10    0    0    1    1    0    12.66    1    1.1    0.9;
0037     16    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0038     17    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0039     18    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0040     19    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0041     20    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0042     21    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0043     22    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0044     23    1    90    50    0    0    1    1    0    12.66    1    1.1    0.9;
0045     24    1    420    200    0    0    1    1    0    12.66    1    1.1    0.9;
0046     25    1    420    200    0    0    1    1    0    12.66    1    1.1    0.9;
0047     26    1    60    25    0    0    1    1    0    12.66    1    1.1    0.9;
0048     27    1    60    25    0    0    1    1    0    12.66    1    1.1    0.9;
0049     28    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0050     29    1    120    70    0    0    1    1    0    12.66    1    1.1    0.9;
0051     30    1    200    600    0    0    1    1    0    12.66    1    1.1    0.9;
0052     31    1    150    70    0    0    1    1    0    12.66    1    1.1    0.9;
0053     32    1    210    100    0    0    1    1    0    12.66    1    1.1    0.9;
0054     33    1    60    40    0    0    1    1    0    12.66    1    1.1    0.9;
0055 ];
0056 
0057 %% generator data
0058 %    bus    Pg    Qg    Qmax    Qmin    Vg    mBase    status    Pmax    Pmin    Pc1    Pc2    Qc1min    Qc1max    Qc2min    Qc2max    ramp_agc    ramp_10    ramp_30    ramp_q    apf
0059 mpc.gen = [
0060     1    0    0    10    -10    1    100    1    10    0    0    0    0    0    0    0    0    0    0    0    0;
0061 ];
0062 
0063 %% branch data
0064 %    fbus    tbus    r    x    b    rateA    rateB    rateC    ratio    angle    status    angmin    angmax
0065 mpc.branch = [  %% (r and x specified in ohms here, converted to p.u. below)
0066     1    2    0.0922    0.0470    0    0    0    0    0    0    1    -360    360;
0067     2    3    0.4930    0.2511    0    0    0    0    0    0    1    -360    360;
0068     3    4    0.3660    0.1864    0    0    0    0    0    0    1    -360    360;
0069     4    5    0.3811    0.1941    0    0    0    0    0    0    1    -360    360;
0070     5    6    0.8190    0.7070    0    0    0    0    0    0    1    -360    360;
0071     6    7    0.1872    0.6188    0    0    0    0    0    0    1    -360    360;
0072     7    8    0.7114    0.2351    0    0    0    0    0    0    1    -360    360;
0073     8    9    1.0300    0.7400    0    0    0    0    0    0    1    -360    360;
0074     9    10    1.0440    0.7400    0    0    0    0    0    0    1    -360    360;
0075     10    11    0.1966    0.0650    0    0    0    0    0    0    1    -360    360;
0076     11    12    0.3744    0.1238    0    0    0    0    0    0    1    -360    360;
0077     12    13    1.4680    1.1550    0    0    0    0    0    0    1    -360    360;
0078     13    14    0.5416    0.7129    0    0    0    0    0    0    1    -360    360;
0079     14    15    0.5910    0.5260    0    0    0    0    0    0    1    -360    360;
0080     15    16    0.7463    0.5450    0    0    0    0    0    0    1    -360    360;
0081     16    17    1.2890    1.7210    0    0    0    0    0    0    1    -360    360;
0082     17    18    0.7320    0.5740    0    0    0    0    0    0    1    -360    360;
0083     2    19    0.1640    0.1565    0    0    0    0    0    0    1    -360    360;
0084     19    20    1.5042    1.3554    0    0    0    0    0    0    1    -360    360;
0085     20    21    0.4095    0.4784    0    0    0    0    0    0    1    -360    360;
0086     21    22    0.7089    0.9373    0    0    0    0    0    0    1    -360    360;
0087     3    23    0.4512    0.3083    0    0    0    0    0    0    1    -360    360;
0088     23    24    0.8980    0.7091    0    0    0    0    0    0    1    -360    360;
0089     24    25    0.8960    0.7011    0    0    0    0    0    0    1    -360    360;
0090     6    26    0.2030    0.1034    0    0    0    0    0    0    1    -360    360;
0091     26    27    0.2842    0.1447    0    0    0    0    0    0    1    -360    360;
0092     27    28    1.0590    0.9337    0    0    0    0    0    0    1    -360    360;
0093     28    29    0.8042    0.7006    0    0    0    0    0    0    1    -360    360;
0094     29    30    0.5075    0.2585    0    0    0    0    0    0    1    -360    360;
0095     30    31    0.9744    0.9630    0    0    0    0    0    0    1    -360    360;
0096     31    32    0.3105    0.3619    0    0    0    0    0    0    1    -360    360;
0097     32    33    0.3410    0.5302    0    0    0    0    0    0    1    -360    360;
0098     21    8    2.0000    2.0000    0    0    0    0    0    0    0    -360    360;
0099     9    15    2.0000    2.0000    0    0    0    0    0    0    0    -360    360;
0100     12    22    2.0000    2.0000    0    0    0    0    0    0    0    -360    360;
0101     18    33    0.5000    0.5000    0    0    0    0    0    0    0    -360    360;
0102     25    29    0.5000    0.5000    0    0    0    0    0    0    0    -360    360;
0103 ];
0104 
0105 %%-----  OPF Data  -----%%
0106 %% generator cost data
0107 %    1    startup    shutdown    n    x1    y1    ...    xn    yn
0108 %    2    startup    shutdown    n    c(n-1)    ...    c0
0109 mpc.gencost = [
0110     2    0    0    3    0    20    0;
0111 ];
0112 
0113 
0114 %% convert branch impedances from Ohms to p.u.
0115 [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
0116     VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
0117 [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
0118     TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
0119     ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
0120 Vbase = mpc.bus(1, BASE_KV) * 1e3;      %% in Volts
0121 Sbase = mpc.baseMVA * 1e6;              %% in VA
0122 mpc.branch(:, [BR_R BR_X]) = mpc.branch(:, [BR_R BR_X]) / (Vbase^2 / Sbase);
0123 
0124 %% convert loads from kW to MW
0125 mpc.bus(:, [PD, QD]) = mpc.bus(:, [PD, QD]) / 1e3;

Generated on Fri 16-Dec-2016 12:45:37 by m2html © 2005