psse2psat

这是一个很适合研究和学习用的电力系统仿真软件

# -----------------------------------------------------------------------
# write Bus.con
# -----------------------------------------------------------------------
if ($nbus >= 0) {
    print OUT "Bus.con = [ ...\n";
    $format = "%4d %8.2f %8.5f %8.5f %2d %2d;\n";
    for ($i = 0; $i <= $nbus-1; $i++) {
	printf OUT $format,$busidx[$i],$kvb[$i],$ein[$i],$ang[$i],
	$kae[$i],$kzo[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write SW.con
# -----------------------------------------------------------------------
if ($nsw >= 0) {
    print OUT "SW.con = [ ...\n";
    $format = "%4d %7.2f %6.2f %8.5f 0.00 %8.5f %8.5f 1.1 0.9 %8.5f 1];\n\n";
    printf OUT $format,$swbus,$swsb,$kvbase{$swbus},
    $swvs,$swqt,$swqb,$swpg;
}

# -----------------------------------------------------------------------
# write PV.con
# -----------------------------------------------------------------------
if ($npv >= 0) {
    printf OUT "PV.con = [ ...\n";
    $format = "%4d %7.2f %6.2f " . "%8.5f " x 4 . "1.1 0.9 1;\n";
    for ($i = 0; $i <= $npv; $i++) {
	printf OUT $format,$pvbus[$i],$pvsb[$i],$kvbase{$pvbus[$i]},
	$pvpg[$i],$pvvs[$i],$pvqt[$i],$pvqb[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write PQ.con
# -----------------------------------------------------------------------
if ($npq >= 0) {
    printf OUT "PQ.con = [ ...\n";
    $format = "%4d $pbas %8.2f %8.5f %8.5f 1.1 0.9 1;\n";
    for ($i = 0; $i <= $npq; $i++) {
	$h = $pqbus[$i];
	printf OUT $format,$h,$kvbase{$h},$pqpl[$i],$pqql[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Pl.con
# -----------------------------------------------------------------------
if ($npl >= 0) {
    print OUT "Pl.con = [ ...\n";
    $format = "%4d $pbas %8.2f 60 %8.5f %8.5f 0 %8.5f %8.5f 0 0;\n";
    for ($i = 0; $i <= $npl; $i++) {
	printf OUT $format,$plbus[$i],$kvbase{$plbus[$i]},$plyp[$i],
	$plip[$i],$plyq[$i],$pliq[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Shunt.con
# -----------------------------------------------------------------------
if ($nsh >= 0) {
    print OUT "Shunt.con = [ ...\n";
    $format = "%4d $pbas %8.2f 60 %8.5f %8.5f;\n";
    for ($i = 0; $i <= $nsh; $i++) {
	if ($psh[$i] != 0 or $qsh[$i] != 0) {
	    $h = abs($idxsh[$i]);
	    printf OUT $format,$h,$kvbase{$h},$psh[$i],$qsh[$i];
	}
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Line.con
# -----------------------------------------------------------------------
if ($nline >= 0) {
    print OUT "Line.con = [ ...\n";
    $format = "%0.1s %4d %4d %8.2f %8.2f 60 0 %5.2f " . "%8.5f " x 5 . 
	" %8.3f %8.3f %8.3f;\n";
    for ($i = 0; $i <= $nline; $i++) {
	if ($busidw{$bustt[$i]}) {
	    if ($kvbase{$bustt[$i]} != 0) {
		$k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};  
	    } else { 
		$k = 0;
		print "WARNING: Check voltage rating of bus ".$busidw{$bustt[$i]}."\n";
	    }    
	    if ($k == 1) { $k = 0; }
	    if ($stt[$i] == 1) { 
		$stt = ' ';
	    } else {
		$stt = '%';
	    }
	    printf OUT $format,$stt,$busfr[$i],$bustt[$i],$lbas[$i],
	    $kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],$ratio[$i],$phas[$i],
	    $ratea[$i],$rateb[$i],$ratec[$i];
	} else {
	    print "WARNING: Bus ".$bustt[$i]." does not exists -> Line ".$i." ignored.\n";
	}
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Ltc.con
# -----------------------------------------------------------------------
if ($nltc >= 0) {
    print OUT "Ltc.con = [ ...\n";
    $format = "%4d %4d %8.2f %8.2f 60 1 0.001 0.1 " . 
	"%8.5f " x 6 . " %4d %4d;\n";
    for ($i = 0; $i <= $nltc; $i++) {
	printf OUT $format,$ltc_busfr[$i],$ltc_busto[$i],$ltc_bas[$i],
	$kvbase{$ltc_busfr[$i]},$ltc_vma[$i],$ltc_vmi[$i],$ltc_dm[$i],
	$ltc_ref[$i],$ltc_x[$i],$ltc_r[$i],$ltc_busct[$i],$ltc_ctr[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Phs.con
# -----------------------------------------------------------------------
if ($nphs >= 0) {
    print OUT "Phs.con = [ ...\n";
    $format = "%4d %4d %8.2f %8.2f %8.2f 60 0.001 10 50 " . 
	"%8.5f " x 5 . ";\n";
    for ($i = 0; $i <= $nphs; $i++) {
	printf OUT $format,$phs_busfr[$i],$phs_busto[$i],$phs_bas[$i],
	$kvbase{$ltc_busfr[$i]},$kvbase{$ltc_busto[$i]},$phs_ref[$i],
	$phs_r[$i],$phs_x[$i],$phs_vma[$i],$phs_vmi[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Twt.con
# -----------------------------------------------------------------------
if ($cmtwt > 0) {
    print OUT "% Commented lines are disconnected 3 winding transformers\n";
}  
if ($flag_version == V29 && $ntw >= 0) {
    print OUT "Twt.con = [ ...\n";
    $format = "%0.1s %4d %4d %4d $pbas 60 %5.2f %5.2f %5.2f " . "%8.5f " x 6 . 
	" 1 " . "%7.3f " x 9 . ";\n";
    for ($i = 0; $i <= $ntw; $i++) {
	if ($stt[$i] == 1) {
	    $stt = ' ';
	} else {
	    $stt = '%';
	}
	printf OUT $format,$stt,$bustw1[$i],$bustw2[$i],$bustw3[$i],
	$kvbase{$bustw1[$i]},$kvbase{$bustw2[$i]},$kvbase{$bustw3[$i]},
	$r12[$i],$r13[$i],$r23[$i],$x12[$i],$x13[$i],$x23[$i],
	$rt1a[$i],$rt1b[$i],$rt1c[$i],$rt2a[$i],$rt2b[$i],$rt2c[$i],
	$rt3a[$i],$rt3b[$i],$rt3c[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Svc.con
# -----------------------------------------------------------------------
if ($nsvc >= 0) {
    print OUT "Svc.con = [ ...\n";
    $format = "%4d $pbas %8.2f 60 1 1 10 " . "%8.5f " x 3 . ";\n";
    for ($i = 0; $i <= $nsvc; $i++) {
	printf OUT $format,$svcbus[$i],$kvbase{$svcbus[$i]},
	$svcref[$i],$svcbma[$i],$svcbmi[$i];
    }
    print OUT "   ];\n\n";
}

# -----------------------------------------------------------------------
# write Varname.bus
# -----------------------------------------------------------------------
for ($i = 0; $i <= $nbus; $i++) {
    $busname[$i] =~ s/\'//g;
    $busname[$i] = deblank ($busname[$i]);    
}
$nbus >= 0 && print OUT "Varname.bus = { ...\n";
$h = ($nbus+1) % 5;
if ($h == 0) {$h = 5;}
if (($nbus+1) > 5) {
    for ($i = 0; $i <= $nbus-$h; $i+=5) {
	print OUT "  '$busname[$i]'; '$busname[$i+1]'; " . 
	    "'$busname[$i+2]'; '$busname[$i+3]'; '$busname[$i+4]';\n";
    }
}
print OUT "  ";
for ($i = $nbus-$h+1; $i <= $nbus-1; $i++) {
    print OUT "'$busname[$i]'; ";
}
print OUT "'$busname[$nbus]'};\n\n";

# -----------------------------------------------------------------------
# close output data file
# -----------------------------------------------------------------------
close(OUT) || die "cannot close $ARGV[1]: $!\n";
print "Coversion completed assuming PSS/E V". $flag_version ." data format.\n";
   
#########################################################################
#########################################################################
   
# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub assign {
    my $param;
    $param = substr($_[0],$_[1],$_[2]);
    $param =~ s/\s*//g;
    unless ($param) {$param = "0";}
    unless ($param =~ /[1-9]/) {$param = $_[3];}
    return $param;
}

# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub deblank {
    my $name = $_[0];
    $name =~ s/( *)$//;
    return $name;
}

# -----------------------------------------------------------------------
# function for writing a separator
# -----------------------------------------------------------------------
sub carret {
    $verbose && print "-" x 50 . "\n";
}

# -----------------------------------------------------------------------
# function for reading comma separated data
# -----------------------------------------------------------------------
sub mysplit {
    if ($_ =~  /\s*,\s*/) {
	my $string = $_[0];
	
	# this is to parse line data in PSS/E V23,
	# V24 and V25 data files
	if (/,,,/)  { 
	    $string =~ s/\/(.*)//;
	    if ($V26_up == 0) { 
		$string =~ s/,,,/ /;
		my @mydata = split /\s+/, $string; 
		if ($mydata[0] eq "") { shift(@mydata); } 
		return @mydata; 
	    } else {
		# this is to parse line data in PSS/E V26 
		# data files	   
		$string =~ s/,,,/,/;
		$string =~ s/\s*//g;
		my @mydata = split /,/, $string; 
		if ($mydata[0] eq "") { shift(@mydata); } 
		return @mydata;    
	    }
	}
   
	# this is to parse data in PSS/E V23,
	# to V30 data files
	$string =~ s/\/(.*)//; 
	my @mydata = split /\s*,\s*/, $string; 
	if ($mydata[0] eq "") { shift(@mydata); } 
	return @mydata; 
    } elsif ($_ =~  /\s+/)  {
	my $string = $_[0];
	$string =~ s/\/(.*)//; 
	$string =~ s/\'//g; 
	my @mydata = split /\s+/, $string; 
	if ($mydata[0] eq "") { shift(@mydata); } 
	return @mydata; 
    } 
}

sub mysplitdata { 
    my $nn = 0;
    my @localdata;
    $_[0] =~ s/://g;
    $_[0] =~ s/,//g;
    $_[0] =~ s/\///g;
    last if /^\s*0/;
    
    $nn = 0;
    while ($_[0]) {
	$_[0] =~ s/\s*-?\d+[\.,]?\d*\s*|\s*\'.*?\'\s*//;
	$localdata[$nn] = $&;
	$localdata[$nn] =~ s/^\s*\'?//;
	$localdata[$nn] =~ s/,/\./;
	$localdata[$nn] =~ s/\s*\'?\s*$//;
	$nn++;
	if ($nn > 100) {die "Bad input string (infinite regexp loop...)";}
    }
    return @localdata;
}
# -----------------------------------------------------------------------
# function for finding multiple PQ loads at the same bus
# -----------------------------------------------------------------------
sub pqidx {
    my $mypqbus = $_[0];
    my $check = -1;
    
    for ($i = 0; $i <= $npq; $i++) {
	if ($pqbus[$i] == $mypqbus) {
	    $check = $i; 
	    last;
	}	    
    }	
    return $check;
}

# -----------------------------------------------------------------------
# function for finding multiple PV generators at the same bus
# -----------------------------------------------------------------------
sub pvidx {
    my $mypvbus = $_[0];
    my $check = -1;
    
    for ($i = 0; $i <= $npv; $i++) {
	if ($pvbus[$i] == $mypvbus) {
	    $check = $i; 
	    last;
	}	    
    }	
    return $check;
}

# -----------------------------------------------------------------------
# function for determining if the bus is isolated
# -----------------------------------------------------------------------
sub ididx {
    my $myisbus = $_[0];
    my $check = 0;
    
    for ($i = 0; $i <= $nisl; $i++) {
	if ($islbus[$i] == $myisbus) {
	    $check = 1; 
	    last;
	}	    
    }	
    return $check;
}