data16m.m

潮流计算

	31  62  0.      0.026   0.      1.04   0.;
	32  63  0.      0.013   0.      1.04   0.;
	36  64  0.      0.0075  0.      1.04   0.;
	37  65  0.      0.0033  0.      1.04   0.;
	41  66  0.      0.0015  0.      1.     0.;
	42  67  0.      0.0015  0.      1.     0.;
	52  68  0.      0.0030  0.      1.     0.;
	 1  27  0.032   0.32    0.41    1.     0.];


% Machine data format
%       1. machine number,
%       2. bus number,
%       3. base mva,
%       4. leakage reactance x_l(pu),
%       5. resistance r_a(pu),
%       6. d-axis sychronous reactance x_d(pu),
%       7. d-axis transient reactance x'_d(pu),
%       8. d-axis subtransient reactance x"_d(pu),
%       9. d-axis open-circuit time constant T'_do(sec),
%      10. d-axis open-circuit subtransient time constant
%                T"_do(sec),
%      11. q-axis sychronous reactance x_q(pu),
%      12. q-axis transient reactance x'_q(pu),
%      13. q-axis subtransient reactance x"_q(pu),
%      14. q-axis open-circuit time constant T'_qo(sec),
%      15. q-axis open circuit subtransient time constant
%                T"_qo(sec),
%      16. inertia constant H(sec),
%      17. damping coefficient d_o(pu),
%      18. dampling coefficient d_1(pu),
%      19. bus number
%      20. saturation factor S(1.0)
%      21. saturation factor S(1.2)
% note: all the following machines use subtransient reactance model

mac_con = [
1 53  100 0.0125 0.0  0.1    0.031   0.025    10.2  0.05 0.069  0.028   0.025   1.5   0.035 42.   4.0  0  53   0 0;
2 54  100 0.035  0.0  0.295  0.0697  0.05     6.56  0.05 0.282  0.060   0.05    1.5   0.035 30.2  9.75 0  54   0 0;
3 55  100 0.0304 0.0  0.2495 0.0531  0.045    5.7   0.05 0.237  0.050   0.045   1.5   0.035 35.8  10   0  55   0 0;
4 56  100 0.0295 0.0  0.262  0.0436  0.035    5.69  0.05 0.258  0.040   0.035   1.5   0.035 28.6  10   0  56   0 0;
5 57  100 0.027  0.0  0.33   0.066   0.05     5.4   0.05 0.31   0.060   0.05    0.44  0.035 26.   3    0  57   0 0;
6 58  100 0.0224 0.0  0.254  0.05    0.04     7.3   0.05 0.241  0.045   0.04    0.4   0.035 34.8  10   0  58   0 0;
7 59  100 0.0322 0.0  0.295  0.049   0.04     5.66  0.05 0.292  0.045   0.04    1.5   0.035 26.4  8    0  59   0 0;
8 60  100 0.028  0.0  0.29   0.057   0.045    6.7   0.05 0.280  0.050   0.045   0.41  0.035 24.3  9    0  60   0 0;
9 61  100 0.0298 0.0  0.2106 0.057   0.045    4.79  0.05 0.205  0.050   0.045   1.96  0.035 34.5  14   0  61   0 0;
10 62 100 0.0199 0.0  0.169  0.0457  0.04     9.37  0.05 0.115  0.045   0.04    1.5   0.035 31.0  5.56 0  62   0 0;
11 63 100 0.0103 0.0  0.128  0.018   0.012    4.1   0.05 0.123  0.015   0.012   1.5   0.035 28.2  13.6 0  63   0 0;
12 64 100 0.022  0.0  0.101  0.031   0.025    7.4   0.05 0.095  0.028   0.025   1.5   0.035 92.3  13.5 0  64   0 0;
13 65 200 0.0030 0.0  0.0296 0.0055  0.004    5.9   0.05 0.0286 0.005   0.004   1.5   0.035 248.0 33   0  65   0 0;
14 66 100 0.0017 0.0  0.018  0.00285 0.0023   4.1   0.05 0.0173 0.0025  0.0023  1.5   0.035 300.0 100  0  66   0 0;
15 67 100 0.0017 0.0  0.018  0.00285 0.0023   4.1   0.05 0.0173 0.0025  0.0023  1.5   0.035 300.0 100  0  67   0 0;
16 68 200 0.0041 0.0  0.0356 0.0071  0.0055  7.8    0.05 0.0334 0.006   0.0055 1.5   0.035 225.0 50   0  68   0 0] ;

% exciter data dc12 model
%     1 - exciter type (1 for DC1, 2 for DC2)
%     2 - machine number
%     3 - input filter time constant T_R
%     4 - voltage regulator gain K_A
%     5 - voltage regulator time constant T_A
%     6 - voltage regulator time constant T_B
%     7 - voltage regulator time constant T_C
%     8 - maximum voltage regulator output V_Rmax
%     9 - minimum voltage regulator output V_Rmin
%    10 - exciter constant K_E
%    11 - exciter time constant T_E
%    12 - E_1
%    13 - saturation function S_E(E_1)
%    14 - E_2
%    15 - saturation function S_E(E_2)
%    16 - stabilizer gain K_F
%    17 - stabilizer time constant T_F
exc_con1 = [...
1 1 0 5.0 0.06 0 0 10.0  -2.0  -0.0485 0.25 3.5461 0.080 4.7281 0.260 0.040 1.0;
1 2 0 6.2 0.05 0 0 10.0  -2.0  -0.0633 .405 0.9183 0.660 1.2244 0.880 0.057 0.5;
1 3 0 5.0 0.06 0 0 10.0  -2.0  -.0198 0.50 2.3423 0.130 3.1230 0.340 0.080 1.0;
1 4 0 5.0 0.06 0 0 10.0  -2.0  -.0525 0.50 2.8681 0.080 3.8241 0.314 0.080 1.0;
1 5 0 40. 0.02 0 0 10.  -10.    1.0   .785 3.9267 0.070 5.2356 0.910 0.030 1.0;
1 6 0 5.0 0.02 0 0 10.0  -2.0  -.0419 .471 3.5868 0.064 4.7824 0.251 .0754 1.246;
1 7 0 40. 0.02 0 0 6.5  -6.5    1.0   .730 2.8017 0.530 3.7356 0.740 0.030 1.0;
1 8 0 5.0 0.02 0 0 10.0  -2.0  -0.047 .528 3.1915 0.072 4.2553 0.282 .0854 1.26;
1 9 0 40. 0.02 0 0 10.5 -10.5   1.0   1.40 3.2568 0.0620 4.6757 0.850 0.030 1.0 ];

% modified exc_con;
exc_con = [...
1 1 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 2 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 3 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 4 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 5 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 6 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 7 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
1 8 0 40. 0.02 0 0 10.  -10.  1.0    .785 3.9267 0.070 5.2356 0.910 0.030 1.0 0 0 0;
0 9 0.01 200.0  0     0     0    5.0  -5.0...
    0    0      0     0     0    0    0      0      0    0   0];
pss_con = [...
1 9 100.0  10.0  0.2  0.1  0.2  0.1  0.2 -0.05];

%Switching file defines the simulation control
% row 1 col1  simulation start time (s) (cols 2 to 6 zeros)
%       col7  initial time step (s)
% row 2 col1  fault application time (s)
%       col2  bus number at which fault is applied
%       col3  bus number defining far end of faulted line
%       col4  zero sequence impedance in pu on system base
%       col5  negative sequence impedance in pu on system base
%       col6  type of fault  - 0 three phase
%                            - 1 line to ground
%                            - 2 line-to-line to ground
%                            - 3 line-to-line
%                            - 4 loss of line with no fault
%                            - 5 loss of load at bus
%       col7  time step for fault period (s)
% row 3 col1  near end fault clearing time (s) (cols 2 to 6 zeros)
%       col7  time step for second part of fault (s)
% row 4 col1  far end fault clearing time (s) (cols 2 to 6 zeros)
%       col7  time step for fault cleared simulation (s)
% row 5 col1  time to change step length (s)
%       col7  time step (s)
%
%
%
% row n col1 finishing time (s)  (n indicates that intermediate rows may be inserted)

sw_con = [...
0    0    0    0    0    0    0.01;%sets intitial time step
0.1  29    28    0    0    0    0.005; %apply three phase fault at bus 29, on line 29-28
0.15 0    0    0    0    0    0.005; %clear fault at bus 
0.20 0    0    0    0    0    0.005; %clear remote end
0.50 0    0    0    0    0    0.01; % increase time step 
1.0  0    0    0    0    0   0.02; % increase time step
5.0  0    0    0    0    0    0]; % end simulation
%
%
% non-conforming load
% col 1           bus number
% col 2           fraction const active power load
% col 3           fraction const reactive power load
% col 4           fraction const active current load
% col 5           fraction const reactive current load

% non-conforming load on all load buses
load_con = [...
1  0 0 .5 0;
3  0 0 .5 0;
4  0 0 .5 0;
7  0 0 .5 0;
8  0 0 .5 0;
9  0 0 .5 0;
15 0 0 .5 0;
16 0 0 .5 0;
18 0 0 .5 0;
20 0 0 .5 0;
21 0 0 .5 0;
23 0 0 .5 0;
24 0 0 .5 0;
25 0 0 .5 0;
26 0 0 .5 0;
27 0 0 .5 0;
28 0 0 .5 0;
29 0 0 .5 0;
33 0 0 .5 0;
37 0 0 .5 0;
39 0 0 .5 0;
40 0 0 .5 0;
41 0 0 .5 0;
42 0 0 .5 0;
44 0 0 .5 0;
45 0 0 .5 0;
46 0 0 .5 0;
47 0 0 .5 0;
48 0 0 .5 0;
49 0 0 .5 0;
50 0 0 .5 0;
51 0 0 .5 0;
52 0 0 .5 0];