simpow2psat

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

    $mstt[$ndyn] = 2 * assign($_,'H',1,10);
    $damp[$ndyn] = assign($_,'DAMP',4,0);
        
    # read third line
    $_ = <IN>;
    chomp;
    $td0b[$ndyn] = assign($_,'TD0S',4,0);
    $tq0b[$ndyn] = assign($_,'TQ0S',4,0);
    $mstt[$ndyn] = 2 * assign($_,'H',1,10);
    $damp[$ndyn] = assign($_,'DAMP',4,0); 
    
    # read fourth line
    $_ = <IN>;
    chomp;
    $mstt[$ndyn] = 2 * assign($_,'H',1,10);
    $damp[$ndyn] = assign($_,'DAMP',4,0);
    $gbas[$ndyn] = assign($_,'SN',2,0);    
    $genbas{$synidx[$ndyn]} = $gbas[$ndyn];
}

# -----------------------------------------------------------------------
# close DYNPOW/SIMPOW data file
# -----------------------------------------------------------------------
close(IN) || die "cannot close $ARGV[0]: $!\n";

# -----------------------------------------------------------------------
# open output data file
# -----------------------------------------------------------------------
print "Writing PSAT file \"$ARGV[1]\"...\n";
carret();
open(OUT,">$ARGV[1]") || die "cannot open $ARGV[1]: $!\n";

# -----------------------------------------------------------------------
# write output data file
# -----------------------------------------------------------------------
print OUT "% File generated by PSAT from SIMPOW data file.\n";
print OUT "% "."-" x 78 . "\n";
print OUT "% Author:   Juan Carlos Morataya\n";
print OUT "% E-mail:   jmorataya\@eegsa.net\n";
print OUT "% "."-" x 78 . "\n";
print OUT "% $title\n";
for ($i = 0; $i <= $ntitle; $i++) { print OUT "% $title[$i]\n"; }
print OUT "\n";

# -----------------------------------------------------------------------
# write Bus.con
# -----------------------------------------------------------------------
$format = "%4d  %8.5f  1   0   1  %2d";
$nbus >= 0 && print OUT "Bus.con = [ ...\n";

for ($i = 0; $i <= $nbus-1; $i++) {
    printf OUT "$format;\n",$busidx[$i],$kvb[$i],
    $kzo[$i];
}

$nbus >= 0 && printf OUT "$format];\n\n",$busidx[$i],$kvb[$i],
$kzo[$i];

# -----------------------------------------------------------------------
# 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";
    printf OUT "$format];\n\n", $swbus,$swsb,$kvbase{$swbus},
    $swvs,$swqt,$swqb,$swpg;
}

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

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

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

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

# -----------------------------------------------------------------------
# write Line.con
# -----------------------------------------------------------------------
if ($nline >= 0) {
    print OUT "Line.con = [ ...\n";
    $format = "%4d %4d " . " %8.5f " x 2 ."  60   0 " . " %8.5f " x 4 . 
    "" . " %8.5f " x 5 . "";
    for ($i = 0; $i <= $nline-1; $i++) {
    $k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};
    if ($k == 1) { $k = 0; }
    else { $pbas = $tfsb[$i]; }
    printf OUT "$format;\n", $nameidx{$busfr[$i]},$nameidx{$bustt[$i]},
    $pbas,$kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],
    $ratio[$i],$phas[$i],$ratea[$i],$rateb[$i],$ratec[$i];
    }
    $i = $nline; 
    $k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};
    if ($k == 1) { $k = 0; }
    else { $pbas = $tfsb[$i]; }
    printf OUT "$format];\n\n",$nameidx{$busfr[$i]},$nameidx{$bustt[$i]},
    $pbas,$kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],
    $ratio[$i],$phas[$i],$ratea[$i],$rateb[$i],$ratec[$i];
}    

# -----------------------------------------------------------------------
# write Twt.con
# -----------------------------------------------------------------------
$ntw >= 0 && print OUT "Twt.con = [ ...\n";
$format = "%4d %4d %4d $pbas 60 %5.2f %5.2f %5.2f " . "%8.5f " x 6 . 
" 1 " . "%7.3f " x 9;
for ($i = 0; $i <= $ntw-1; $i++) {
    printf OUT "$format;\n", $nameidx{$bustw1[$i]},$nameidx{$bustw2[$i]},
    $nameidx{$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];
}

$i = $ntw;
if ($ntw >= 0) {
    printf OUT "$format];\n\n", $nameidx{$bustw1[$i]},$nameidx{$bustw2[$i]},
    $nameidx{$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];
}

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

$i = $nsvc;
if ($nsvc >= 0) {
    printf OUT "$format];\n\n",$nameidx{$idxsvc[$i]},
    $kvbase{$idxsvc[$i]},$einsch{$idxsvc[$i]},$svcbma[$i],$svcbmi[$i];
}

# -----------------------------------------------------------------------
# write Syn.con
# -----------------------------------------------------------------------
$ndyn >= 0 && print OUT "Syn.con = [ ...\n";
$format = "%4d %6.2f %6.2f $fsys 6 %6.2f " . "%8.5f " x 13 . 
" " . " 0 " x 2 . " " . " 1 " x 3 . " " . " 0 " x 3;
for ($i = 0; $i <= $ndyn-1; $i++) {
    printf OUT "$format;\n", $nameidx{$synidx[$i]},$gbas[$i],
    $kvbase{$synidx[$i]},$xl[$i],$ra[$i],$xd[$i],$xdp[$i],$xdb[$i],
    $td0p[$i],$td0b[$i],$xq[$i],$xqp[$i],$xqb[$i],$tq0p[$i],
    $tq0b[$i],$mstt[$i],$damp[$i];
}

$i = $ndyn;
if ($ndyn >= 0) {
    printf OUT "$format];\n\n", $nameidx{$synidx[$i]},$gbas[$i],
    $kvbase{$synidx[$i]},$xl[$i],$ra[$i],$xd[$i],$xdp[$i],$xdb[$i],
    $td0p[$i],$td0b[$i],$xq[$i],$xqp[$i],$xqb[$i],$tq0p[$i],
    $tq0b[$i],$mstt[$i],$damp[$i];
}

# -----------------------------------------------------------------------
# write Varname.bus
# -----------------------------------------------------------------------
$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 "Conversion completed.\n";



# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub assign {
    my $param;
    my $tmp;
    my $idxnum = 0;
    my $idxcut = 0;
    my $check;
    $_[0] =~ s/UB=/VOLT=/g; # these are to avoid badly catch values
    $_[0] =~ s/U=/V=/g;
    $_[0] =~ s/R=/REST=/g;
    $_[0] =~ s/X=/REAC=/g;
    $_[0] =~ s/XDB=/XDS=/g;
    $_[0] =~ s/XQB=/XQS=/g;
    $_[0] =~ s/TD0B=/TD0S=/g;
    $_[0] =~ s/TQ0B=/TQ0S=/g;
    $_[0] =~ s/B=/SUSC=/g;
    $_[0] =~ s/L=/LENGTH=/g;
    $_[0] =~ s/RTYP/RTY/g;
    $_[0] =~ s/TYPE/TY/g;
    $_[0] =~ s/P=/PPOW=/g;
    $_[0] =~ s/Q=/QPOW=/g;
    $_[0] =~ s/SW/1/g;
    $_[0] =~ s/UP/2/g;
    $_[0] =~ s/PQ/3/g;
    $_[0] =~ s/D=/DAMP=/g;

    $_[0] =~ s/\s*=\s*//g; 
    
    $idxnum = (index $_[0],$_[1]) + $_[2];
    $check = $idxnum - $_[2];
    for ($i = 1; $i < 40; $i++) {
	     $tmp = substr($_,$idxnum+($i-1),1);
	      if ($tmp =~ /\s/ ) { 
		      $idxcut = $i;
		      last;
	      }
     }
    $param = substr($_[0],$idxnum,$idxcut); 
    $param =~ s/\s*//g;
    unless ($param) {$param = "0";} 
    unless ($param =~ /[1-9]/) { $param = $_[3]; }
    if ( $check < 0 ) { $param = $_[3]; }  # this assigns a determined value
    return $param;                         # when a field/value isn't founded
}

# -----------------------------------------------------------------------
# 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 scanning the file up to the data entries
# -----------------------------------------------------------------------
sub splitdata {
	my $string = $_[0];
	#$string =~ s/\w+=/ /g;
	my @mydata = split /\s+/, $string; 
}

# -----------------------------------------------------------------------
# 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;
}