sub_flow.cpp

电力系统发输电系统基于非序贯的蒙特卡罗抽样的风险评估程序

#include<math.h>
#include<stdlib.h>

#include"datastruct.h"
#include"cut.h"
//计算子系统潮流程序，并进行过载、线路过载判断，负荷削减等操作
void sub_flow(struct Busdata *tempbusdata,struct Branchdata *tempbranchdata,struct Gendata *tempgendata,
			  int Busnum,int Branchnum,int Gennum,int *igen,
			  struct Busdata *sub_busdata,struct Branchdata *sub_branchdata,
			  int sub_busnum,int sub_branchnum)
{
	int *balanceno;//定义平衡节点号
	int judge_balanceno=0;//判断是否有平衡节点，0为没有，1为有平衡节点
	balanceno=new int[1];
	if(sub_busnum==1)///子系统为孤立节点
	{
		if(sub_busdata[0].PL>sub_busdata[0].PG)
			sub_busdata[0].PL=sub_busdata[0].PG;
		//*(sub_curtailment)+=sub_busdata[0].PL-sub_busdata[0].PG;
		//printf("孤立节点：%d\n",sub_busdata[0].busno);
	}
	else if(sub_busnum>1)///非孤立节点系统
	{///1
		double PG=0,PL=0;
		for(int i=0;i<sub_busnum;i++)
		{
			PG+=sub_busdata[i].PG;
			PL+=sub_busdata[i].PL;
		}
		if(PL>0)//子系统有负荷
		{
			if(PG==0)//子系统没有发电机
			{
				for(i=0;i<sub_busnum;i++)
					sub_busdata[i].PL=0;
			}
			//*sub_curtailment+=PL;
			else if(PG>0)//子系统有发电机
			{
				/*先找到平衡节点，没有平衡节点则设置平衡节点*/
				for(i=0;i<sub_busnum;i++)
				{
					if(sub_busdata[i].busflag==-1)
					{
						*balanceno=i;
						judge_balanceno=1;
						break;
					}
				}
				if(judge_balanceno==0)//No balance node;
				{
					double maxPG=0;//如果没有平衡节点，则设置发电机出力最大的节点为平衡节点
					for(i=0;i<sub_busnum;i++)
					{
						if(maxPG<sub_busdata[i].PG)
						{
							maxPG=sub_busdata[i].PG;
							*balanceno=i;
						}
					}
				}
				curtailment(tempbusdata,tempbranchdata,tempgendata,Busnum,Branchnum,Gennum,igen,//原系统数据
				 sub_busdata,sub_branchdata,//子系统数据
				 sub_busnum,sub_branchnum,balanceno);
			}
		}
	}///1
	for(int i=0;i<Busnum;i++)
	{
		for(int j=0;j<sub_busnum;j++)
		{
			if(sub_busdata[j].busno==tempbusdata[i].busno)
			{
				tempbusdata[i].PL=sub_busdata[j].PL;
				break;
			}
		}
	}
	delete []balanceno;
	balanceno=NULL;
}