sample.piecewise.IEEE14BusSample.java Source code

Introduction

Here is the source code for sample.piecewise.IEEE14BusSample.java

Source

/*
 * @(#)IEEE14BusSample.java   
 *
 * Copyright (C) 2006 www.interpss.org
 *
 * This program is free software; you can redistribute it and/or
 * modify it under the terms of the GNU LESSER GENERAL PUBLIC LICENSE
 * as published by the Free Software Foundation; either version 2.1
 * of the License, or (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * @Author Mike Zhou
 * @Version 1.0
 * @Date 04/15/2016
 * 
 *   Revision History
 *   ================
 *
 */

package sample.piecewise;

import java.util.Hashtable;
import java.util.function.Function;

import org.apache.commons.math3.complex.Complex;
import org.interpss.CorePluginObjFactory;
import org.interpss.IpssCorePlugin;
import org.interpss.fadapter.IpssFileAdapter;
import org.interpss.numeric.datatype.ComplexFunc;
import org.interpss.piecewise.CuttingBranch;
import org.interpss.piecewise.impl.PiecewiseAlgorithmImpl;

import com.interpss.CoreObjectFactory;
import com.interpss.core.aclf.AclfBus;
import com.interpss.core.aclf.AclfLoadCode;
import com.interpss.core.aclf.AclfNetwork;
import com.interpss.core.algo.LoadflowAlgorithm;

public class IEEE14BusSample {

    public static void main(String args[]) throws Exception {
        // initialize InterPSS plugin
        IpssCorePlugin.init();

        // =========== Aclf Network data preparation =====================//
        /*
         * The network diagram can be found at https://drive.google.com/drive/folders/0BzjeDvtdQBeyZTA5YTE4NGItMDIwMS00OGYzLTkyNGItOTNlZjM1YTEwMTEw
         */

        // Load the sample network
        AclfNetwork net = CorePluginObjFactory.getFileAdapter(IpssFileAdapter.FileFormat.IpssInternal)
                .load("testdata/ieee14.ipssdat").getAclfNet();

        // calculate loadflow
        LoadflowAlgorithm algo = CoreObjectFactory.createLoadflowAlgorithm(net);
        algo.loadflow();
        //System.out.println(net.net2String());

        // Turn all loads to Constant-Z load to create a linear network
        net.getBusList().forEach(bus -> {
            if (bus.isLoad())
                bus.setLoadCode(AclfLoadCode.CONST_Z);
        });

        // ============ Solve network voltage by piecewise algo ====================//

        // define the cutting branch set
        CuttingBranch[] cuttingBranches = new CuttingBranch[] { new CuttingBranch("4->71(1)"),
                new CuttingBranch("4->91(1)"), new CuttingBranch("5->61(1)") };

        // define bus injection current calculation function
        Function<AclfBus, Complex> injCurrentFunc = bus -> { // this function calculates bus injection current
            // The bus injection current is based on gen bus load flow results.
            return bus.isGen() ? bus.getNetGenResults().divide(bus.getVoltage()) : new Complex(0.0, 0.0);
        };

        // define a piecewise algo object and calculate the network bus voltage
        Hashtable<String, Complex> voltages = new PiecewiseAlgorithmImpl(net).calculateNetVoltage(cuttingBranches,
                injCurrentFunc);

        // output network bus voltage
        voltages.forEach((busId, voltage) -> {
            System.out.println("Bus id:" + busId + ", v = " + ComplexFunc.toStr(voltage));
        });
        /* you should see output like the following
         Bus id:10, v = 0.84642 + j0.39398
         Bus id:1, v = 0.93198 + j0.68774
         ...
         */
    }
}