digsilent2psat

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

	    unless ($pvsb[$nopv]) { $pvsb[$nopv] = $pbase{$symtn{$symname[$nopv]}}; }
	    unless ($pvsb[$nopv]) { $pvsb[$nopv] = 100; }
	    $pvpg[$nopv] = assign("OP_ElmSym",'pgini',1)/$pbas;
	    $pvqg[$nopv] = assign("OP_ElmSym",'qgini',1)/$pbas;
	    $pvqt[$nopv] = assign("OP_ElmSym",'q_min',1)/$pbas;
	    $pvqb[$nopv] = assign("OP_ElmSym",'q_max',1)/$pbas;
	    $pvvs[$nopv] = assign("OP_ElmSym",'usetp',1);
	    # assign slack generator
	    $nsw = $busidxn{$pvidx{$symname[swgen($pvpg[0])]}};
	    $swbus = assign("OP_ElmSym",'Name',1);
	    $symstt[$nsw] = assign("OP_ElmSym",'outserv',1);
	    $opmode[$nsw] = assign("OP_ElmSym",'iv_mode',1);
	    $swsb = $pbase{$symtn{$symname[swgen($pvpg[0])]}};
	    unless ($swsb) { $swsb = 100; }
	    $swpg = $pvpg[swgen($pvpg[0])];
	    $swvs = $pvvs[swgen($pvpg[0])];
	    $swqt = $pvqt[swgen($pvpg[0])];
	    $swqb = $pvqb[swgen($pvpg[0])];
	}
    }
}

# <$$OP_ElmTr2>; data card
# two winding transformer operational data 
$j = matching("OP_ElmTr2");
if ($j == $n && $flgchk < 0) { 
    while (<IN>) { 
        next if /^\*+/;
        next if /^\*+\s+\w+/;
        next if /^.$/;
        last if /^\$\$/;
        @data = mysplit($_);
    }
}

# <$$OP_ElmTr3> data card
# three winding operational data
$j = matching("OP_ElmTr3");
if ($j == $n && $flgchk < 0) { 
    while (<IN>) { 
        next if /^\*+/;
        next if /^\*+\s+\w+/;
        next if /^.$/;
        last if /^\$\$/;
        @data = mysplit($_);
    }
}

# <$$OP_ElmXnet> data card
$j = matching("OP_ElmXnet");
if ($j == $n && $flgchk < 0) { 
    while (<IN>) { 
        next if /^\*+/;
        next if /^\*+\s+\w+/;
        next if /^.$/;
        last if /^\$\$/;
    }
}

# <$$OP_StaSwitch> data card
$j = matching("OP_StaSwitch");
if ($j == $n && $flgchk < 0) { 
    while (<IN>) { 
        next if /^\*+/;
        next if /^\*+\s+\w+/;
        next if /^.$/;
        last if /^\$\$/;
    }
}
}

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

# -----------------------------------------------------------------------
# open output data file
# -----------------------------------------------------------------------
if ($flgchk < 0) {
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 DIGSILENT data file.\n";
print OUT "% "."-" x 72 . "\n";
print OUT "% Author:   Juan Carlos Morataya\n";
print OUT "% E-mail:   jc.morataya\@ieee.org\n";
print OUT "% "."-" x 72 . "\n";
for ($j = 0; $j <= $ntitle; $j++) {
print OUT "% $title[$j]\n";
}
print OUT "% "."-" x 72 . "\n\n";


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

for ($i = 0; $i <= $nbus-1; $i++) {
    printf OUT "$format;\n",$busidx[$i],$kvb[$i];
}
$i = $nbus;
$nbus >= 0 && printf OUT "$format];\n\n",$busidx[$i],$kvb[$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", $nsw,$pbas,$kvbase{$nsw},
    $swvs,$swqt,$swqb,$swpg;
    $nopv = $nopv - 1;
}
    
# -----------------------------------------------------------------------
# write PV.con
# -----------------------------------------------------------------------
$nopv >= 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 <= $nopv-1; $i++) {
	$j = $symname[$i];
    $h = $busidxn{$pvidx{$j}};
    if ($h != $nsw) {  
    printf OUT "$format;\n",$h,$pbas,$kvbase{$h},$pvpg[$i],
    $pvvs[$i],$pvqt[$i],$pvqb[$i]; 
    }
}
$i = $nopv;
$j = $symname[$i];
$h = $busidxn{$pvidx{$j}};
$pb = $pbase{$symtn{$j}};
if ($h != $nsw) {  
    printf OUT "$format];\n\n",$h,$pbas,$kvbase{$h},$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 <= $nopq-1; $i++) {
    $h = $pqname[$i];
    $k = $busidxn{$pqidx{$h}};
    printf OUT "$format;\n",$k,$kvbase{$k},$pqpl[$i],$pqql[$i];
}
$h = $pqname[$nopq];
$k = $busidxn{$pqidx{$h}};
printf OUT "$format];\n\n",$k,$kvbase{$k},$pqpl[$i],$pqql[$i];

# -----------------------------------------------------------------------
# write Mn.con
# -----------------------------------------------------------------------
if ($nmn >= 0) {
    printf OUT "Mn.con = [ ...\n";
    $format = "%4d $pbas %8.2f  1  1  %2.2f %2.2f 1";
    for ($i = 0; $i <= $nopq-1; $i++) {
	$h = $pqname[$i];
	$j = $pqtype{$pqidx{$h}};
	$k = $busidxn{$pqidx{$h}};
	printf OUT "$format;\n",$k,$kvbase{$k},$ktep{$j},$kteq{$j};
    }
    $h = $pqname[$nopq];
    $j = $pqtype{$pqidx{$h}};
    $k = $busidxn{$pqidx{$h}};
    printf OUT "$format];\n\n",$k,$kvbase{$k},$ktep{$j},$kteq{$j};
}

# -----------------------------------------------------------------------
# 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++) {
    $h = $shbus[$i];
    $j = $busidxn{$h};
    $k = $pbas/$kvbase{$j}**2;
    printf OUT "$format;\n", $busidxn{$h},$kvbase{$j},
    $k*$pshnt[$i],$k*$qshnt[$i];
}
$i = $nsh;
$h = $shbus[$i];
$j = $busidxn{$h};
$k = $pbas/$kvbase{$j}**2;
printf OUT "$format];\n\n", $busidxn{$h},$kvbase{$j},
$k*$pshnt[$i],$k*$qshnt[$i];

# -----------------------------------------------------------------------
# write Line.con
# -----------------------------------------------------------------------
$nline >= 0 && print OUT "Line.con = [ ...\n";
$format = "%4d %4d %8.2f %8.2f 60 0 %5.2f " . "%8.5f " x 5 . 
    " %8.3f %8.3f %8.3f";
for ($i = 0; $i <= $nline-1; $i++) {
	$h = $linename[$i];
	$k = $kvbase{$busidxn{$busfr[$i]}}/$kvbase{$busidxn{$bustt[$i]}};
	$j = $pbase{$h}/$kvbase{$busidxn{$busfr[$i]}}**2;
    if ($k == 1) {$k = 0;} 
    printf OUT "$format;\n",$busidxn{$busfr[$i]},$busidxn{$bustt[$i]},
    $pbase{$h},$kvbase{$busidxn{$busfr[$i]}},$k,$j*$rest{$h},$j*$reat{$h},
    0.000001*$j*$susc{$h},0,$phsf{$h},$lnerate{$h},0,0;
}
    $i = $nline;
	$h = $linename[$nline];
	$k = $kvbase{$busidxn{$busfr[$i]}}/$kvbase{$busidxn{$bustt[$i]}};
	$j = $pbase{$h}/$kvbase{$busidxn{$busfr[$i]}}**2;
    if ($k == 1) {$k = 0;} 
    printf OUT "$format];\n\n",$busidxn{$busfr[$i]},$busidxn{$bustt[$i]},
    $pbase{$h},$kvbase{$busidxn{$busfr[$i]}},$k,$j*$rest{$h},$j*$reat{$h},
    0.000001*$j*$susc{$h},0,$phsf{$h},$lnerate{$h},0,0;

# -----------------------------------------------------------------------
# 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  0  0 0 0 0 0 " . "%7.3f " x 9;
for ($i = 0; $i <= $ntw-1; $i++) {
	 $h = $tr3wname[$i];
     printf OUT "$format;\n",$busidxn{$bustw1[$i]},$busidxn{$bustw2[$i]},
     $busidxn{$bustw3[$i]},$kvbase{$busidxn{$bustw1[$i]}},
     $kvbase{$busidxn{$bustw2[$i]}},$kvbase{$busidxn{$bustw3[$i]}},
     $r12{$h},$r13{$h},$r23{$h},$x12{$h},$x13{$h},$x23{$h},
     $rt3a{$h},$rt3b{$h},$rt3c{$h};
}
$i = $ntw;
$h = $tr3wname[$i];
     printf OUT "$format];\n\n",$busidxn{$bustw1[$i]},$busidxn{$bustw2[$i]},
     $busidxn{$bustw3[$i]},$kvbase{$busidxn{$bustw1[$i]}},
     $kvbase{$busidxn{$bustw2[$i]}},$kvbase{$busidxn{$bustw3[$i]}},
     $r12{$h},$r13{$h},$r23{$h},$x12{$h},$x13{$h},$x23{$h},
     $rt3a{$h},$rt3b{$h},$rt3c{$h};

# -----------------------------------------------------------------------
# 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";
} else {
print "WARNING: Not enough arguments to parse data in "."\"$nmcrds[$crd]\""." data card...\n";
print "Conversion failed.\n";;
}
	
# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub assign {
	my $param;
	my $idxnum = 0;
    my $str == ();
    my $idxspt = 0;
	my $mycard = $_[0];
	if ($crdord{$_[0]} =~ /$_[1]/) {
    $idxnum = index $crdord{$_[0]},$_[1];
    $str = substr ($crdord{$_[0]},0,$idxnum);
    $idxspt = $str =~ tr/;// -1 ;
    $param = $data[$idxspt];
    #if (!/\w*\W*\s*\w\W*/) { $param =~ s/\s*//g; }
    #$param =~ s/\s*//g; #conflicto con el nombre
    unless ($param) { $param = $_[2]; }
    #if ($param =~ /[\w+\W+]/) { $param = $param; }
    unless($param =~ /[\w+\W+]/) { $param = $_[2]; }
    } else {
    $param = $_[2]; 
    }
    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 {
	my $string = $_[0];
	my @mydata = split /;/,$string; 
	if ($mydata[0] eq "") { shift(@mydata); } 
	return @mydata; 
}

# -----------------------------------------------------------------------
# function for enable data cards according to input data file
# -----------------------------------------------------------------------
sub matching {
    my $mycard = $_[0];
    my $check = -1;

    for ($i = -1; $i <= $ncrd; $i++) {
	if ($nmcrds[$i] eq "$mycard") {
	    $check = $i;
 	    last;
	}	    
    }
    return $check;
    
}

# -----------------------------------------------------------------------
# function for finding multiple PQ loads at the same bus
# -----------------------------------------------------------------------
sub pqidx {
    my $mypqbus = $_[0];
    my $check = -1;
    for ($i = 0; $i <= $nopq; $i++) {
	if ($pqidx{$pqname[$i]} eq "$pqidx{$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 <= $nopv; $i++) {
	if ($pvidx{$symname[$i]} eq "$pvidx{$mypvbus}") {
	    $check = $i; 
	    last;
	}	    
    }	
    return $check;
}


# -----------------------------------------------------------------------
# function for finding multiple PQ loads at the same bus
# -----------------------------------------------------------------------
sub lneidx {
    my $mylne = $_[0];
    my $check = -1;

    for ($i = 0; $i <= $nsec; $i++) {
	if ($linename[$i] eq "$mylne") {
	    $check = $i; 
	    last;
	}	    
    }	
    return $check;
}

# -----------------------------------------------------------------------
# function for finding the slack generator
# -----------------------------------------------------------------------
sub swgen {
    my $swpgen = $_[0];
    my $swpos = 0;
    for ($k = 0; $k <= $nopv; $k++) {
    if ($swpgen >= $pvpg[$k]) { 
    $swpgen = $swpgen;
    } else {
	$swpos = $k;
    $swpgen = $pvpg[$k];
    }
    }
    return $swpos;
}