eu.itesla_project.iidm.network.util.SV.java Source code

Introduction

Here is the source code for eu.itesla_project.iidm.network.util.SV.java

Source

/**
 * Copyright (c) 2016, All partners of the iTesla project (http://www.itesla-project.eu/consortium)
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at http://mozilla.org/MPL/2.0/.
 */
package eu.itesla_project.iidm.network.util;

import eu.itesla_project.iidm.network.DanglingLine;
import eu.itesla_project.iidm.network.Line;
import eu.itesla_project.iidm.network.TwoWindingsTransformer;
import org.apache.commons.math3.complex.Complex;
import org.apache.commons.math3.complex.ComplexUtils;

/**
 * Utility class to compute the state variables on one side of a branch, knowing
 * the state variables on the other side.
 *
 * @author Geoffroy Jamgotchian <geoffroy.jamgotchian at rte-france.com>
 */
public class SV {

    public static float getRatio(TwoWindingsTransformer twt) {
        float rho = twt.getRatedU2() / twt.getRatedU1();
        if (twt.getRatioTapChanger() != null) {
            rho *= twt.getRatioTapChanger().getCurrentStep().getRho();
        }
        return rho;
    }

    public static float getR(TwoWindingsTransformer twt) {
        float r = twt.getR();
        if (twt.getRatioTapChanger() != null) {
            r *= (1 + twt.getRatioTapChanger().getCurrentStep().getR() / 100);
        }
        return r;
    }

    public static float getX(TwoWindingsTransformer twt) {
        float x = twt.getX();
        if (twt.getRatioTapChanger() != null) {
            x *= (1 + twt.getRatioTapChanger().getCurrentStep().getX() / 100);
        }
        return x;
    }

    public static float getG(TwoWindingsTransformer twt) {
        float g = twt.getG();
        if (twt.getRatioTapChanger() != null) {
            g *= (1 + twt.getRatioTapChanger().getCurrentStep().getG() / 100);
        }
        return g;
    }

    public static float getB(TwoWindingsTransformer twt) {
        float b = twt.getB();
        if (twt.getRatioTapChanger() != null) {
            b *= (1 + twt.getRatioTapChanger().getCurrentStep().getB() / 100);
        }
        return b;
    }

    float p;

    float q;

    float u;

    float a;

    public SV(float p, float q, float u, float a) {
        this.p = p;
        this.q = q;
        this.u = u;
        this.a = a;
    }

    public float getP() {
        return p;
    }

    public float getQ() {
        return q;
    }

    public float getU() {
        return u;
    }

    public float getA() {
        return a;
    }

    public SV otherSide(float r, float x, float g, float b, float ratio) {
        Complex z = new Complex(r, x); // z=r+jx
        Complex y = new Complex(g, b); // y=g+jb
        Complex s1 = new Complex(p, q); // s1=p1+jq1
        Complex u1 = ComplexUtils.polar2Complex(u, Math.toRadians(a));
        Complex v1 = u1.divide(Math.sqrt(3f)); // v1=u1/sqrt(3)

        Complex v1p = v1.multiply(ratio); // v1p=v1*rho
        Complex i1 = s1.divide(v1.multiply(3)).conjugate(); // i1=conj(s1/(3*v1))
        Complex i1p = i1.divide(ratio); // i1p=i1/rho
        Complex i2 = i1p.subtract(y.multiply(v1p)); // i2=i1p-y*v1p
        Complex v2 = v1p.subtract(z.multiply(i2)); // v2=v1p-z*i2
        Complex s2 = v2.multiply(3).multiply(i2.conjugate()); // s2=3*v2*conj(i2)

        Complex u2 = v2.multiply(Math.sqrt(3f));
        return new SV((float) -s2.getReal(), (float) -s2.getImaginary(), (float) u2.abs(),
                (float) Math.toDegrees(u2.getArgument()));
    }

    public SV otherSide(TwoWindingsTransformer twt) {
        return otherSide(getR(twt), getX(twt), getG(twt), getB(twt), getRatio(twt));
    }

    public SV otherSide(Line l) {
        return otherSide(l.getR(), l.getX(), l.getG1() + l.getG2(), l.getB1() + l.getB2(), 1);
    }

    public SV otherSide(DanglingLine dl) {
        return otherSide(dl.getR(), dl.getX(), dl.getG(), dl.getB(), 1);
    }

    @Override
    public String toString() {
        return "p=" + p + ", q=" + q + ", u=" + u + ", a=" + a;
    }

}