Re: magnitude and the phase angle of the current in the branches

Sat, 31 Mar 2018 12:41:47 -0700

Assuming that the system data is in structure 'mpc', and that the buses are numbered consecutively starting from 1,  this should give you the injection currents:

% Run power flow
results = runpf(mpc);
% Get column pointers
define_constants;
% Compute vector of complex bus voltages
V = results.bus(:, VM) .* exp(sqrt(-1)*results.bus(:, VA)*pi/180);
% Compute nodal admittance matrix
Ybus = makeYbus(results);
% Compute vector of complex bus injection currents
Ibus = Ybus * V;

Carlos.

Mostafa Mohammadpourfard wrote:
Thank you, Carlos. May I ask you to write the code here? I am not a power system student. Thanks


On Sat, Mar 31, 2018 at 9:41 AM, Carlos E Murillo-Sanchez <[email protected]> wrote:
You have the bus voltages; now you just need to multiply the system's nodal admittance matrix by the complex voltage vector to get the complex currents vector.  See the makeYbus function in MATTPOWER.

Carlos.

Arkan Arkan wrote:
Thank you, Ilias. But, unfortunately, there is not any direct information about current in the Manpower like you mentioned and we should calculate based on the information obtained from runpf.

 On Fri, Mar 30, 2018 at 9:06 AM, Ilias Sarantakos <[email protected] <mailto:[email protected]m>> wrote:

    Hi Arkan,

    If you type ''result = runpf(mpc)'', then you will get the power
    flow results of your network, stored in the ''result'' variable.
    Now that you have the voltages (magnitudes and angles) at all
    buses, you can calculate the current at each branch. I think it
    would be helpful to check Appendix B in MATPOWER manual, which can
    guide you how to access the results of your network. For example,
    ''result.bus(1,8)'' and ''result.bus(1,9)'', contain the magnitude
    and angle of voltage at bus 1.

    Hope this helps.

    Kind regards,

    Ilias Sarantakos

    2018-03-29 2:25 GMT+01:00 Arkan Arkan <[email protected]
     <mailto:[email protected]>>:

        Hi everyone,

        I am wondering how can I get the magnitude and the phase angle
        of the current in the branches of any case file in the Matpower.

        Thanks for your help.







Reply via email to