Rather than setting the reactance to Inf, you could simply set the status of the line to zero ...
mpc.branch(k, BR_STATUS) = 0; If you are running an OPF, there should be no overloads (assuming the problem is feasible). If you are using a simple power flow, you can compute the overloads as follows ... r = runpf(mpc); MVAflow_f = sqrt(r.branch(:, PF).^2 + r.branch(:, QF).^2); MVAflow_t = sqrt(r.branch(:, PT).^2 + r.branch(:, QT).^2); MVAflow_max = max([MVAflow_f MVAflow_t], [], 2); overloads = max(MVAflow_max - r.branch(:, RATE_A), 0); -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Feb 9, 2012, at 11:38 AM, iman wrote: > Dear All, > > I have been trying to study line outages in MATpower.By equaivalenting X of > desired branches to INF (infinity) I made them disconntected from bus > bars.The problem is that no matter which branch is disconected, you never can > see the branch overloads.Does anyone know how should I implement it? > > -- > Best regards > Iman >
