⭐ 欢迎来到虫虫下载站! | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站

📄 psse2psat

📁 电力系统的psat
💻
📖 第 1 页 / 共 3 页
字号:
}
$i = $nbus;
$nbus >= 0 && printf OUT "$format];\n\n", $busidx[$i], $kvb[$i], 
    $ein[$i],$ang[$i],$kae[$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++) {
    #$h = $idxpv[$i];
    printf OUT "$format;\n", $pvbus[$i],$pvsb[$i],$kvbase{$pvbus[$i]},
    $pvpg[$i],$pvvs[$i],$pvqt[$i],$pvqb[$i];
}
#$h = $idxpv[$npv];
$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 = abs($idxsh[$i]);
       printf OUT "$format;\n", $h,$kvbase{$h},$psh[$i],$qsh[$i];
    }
}
$h = $idxsh[$nsh];
$i = $nsh;
$nsh >= 0 && printf OUT "$format];\n\n",$h,$kvbase{$h},$psh[$i],$qsh[$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 3 . 
#    " 0 %8.5f %8.3f %8.3f %8.3f";
$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++) {
    if ($busidw{$bustt[$i]}) {
	if ($kvbase{$bustt[$i]} != 0) {
	    $k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};  
	} else { # print warnings
	    $k = 0;
	    print "WARNING: Check voltage rating of bus ".$busidw{$bustt[$i]}."\n";
	}    
	if ($k == 1) {$k = 0;} 
	printf OUT "$format;\n", $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";	
    }
}
$i = $nline;
if ($nline >= 0) {
    if ($busidw{$bustt[$i]}) {
	if ($kvbase{$bustt[$i]} != 0) {
	    $k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};  
	} else { # print warnings
	    $k = 0;
	    print "WARNING: Check voltage rating of bus ".$busidw{$bustt[$i]}."\n";
	}    
	if ($k == 1) {$k = 0;}
	printf OUT "$format];\n\n",$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";	
    }
}

# -----------------------------------------------------------------------
# write Ltc.con
# -----------------------------------------------------------------------
$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";
for ($i = 0; $i <= $nltc-1; $i++) {
    printf OUT "$format;\n", $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];
}
$i = $nltc;
if ($nltc >= 0) {
    printf OUT "$format];\n\n", $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];
}

# -----------------------------------------------------------------------
# write Phs.con
# -----------------------------------------------------------------------
$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;
for ($i = 0; $i <= $nphs-1; $i++) {
    printf OUT "$format;\n",$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];
}
$i = $nphs;
if ($nphs >= 0) {
    printf OUT "$format];\n\n",$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];
}

# -----------------------------------------------------------------------
# write Twt.con
# -----------------------------------------------------------------------
if ($flag_version == V29) {
    $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", $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];
}

    $i = $ntw;
if ($ntw >= 0) {
    printf OUT "$format];\n\n", $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];
}
}
#my ($nsvc,@svcbus,@svcbma,@svcbmi,@svcref);
# -----------------------------------------------------------------------
# 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",$svcbus[$i],$kvbase{$svcbus[$i]},
    $svcref[$i],$svcbma[$i],$svcbmi[$i];
}
$i = $nsvc;
if ($nsvc >= 0) {
    printf OUT "$format];\n\n",$svcbus[$i],$kvbase{$svcbus[$i]},
    $svcref[$i],$svcbma[$i],$svcbmi[$i];
}


# -----------------------------------------------------------------------
# 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";
if ($pti > 20) {$flag_version = V23;} 
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";}
    #$param =~ s/^0*\.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];
	if (/,,,/)  { 
            # this is to parse branch data only in PSS/E v26 data files
	    $string =~ s/\/(.*)//; 
	    $flag_version = V26;
	    $string =~ s/,,,//;
	    my @mydata = split /\s+/, $string; 
	    if ($mydata[0] eq "") { shift(@mydata); } 
            #print "PSS/E V26 data parsed correctly...\n"; 
	    return @mydata; 
	} else { }
	
	$string =~ s/\/(.*)//; 
	my @mydata = split /\s*,\s*/, $string; 
	if ($mydata[0] eq "") { shift(@mydata); } 
        #print "Data parsed correctly...\n"; 
        return @mydata; 
    } elsif ($_ =~  /\s+/)  {
	my $string = $_[0];
	$string =~ s/\/(.*)//; 
	# if ($_ =~  /\'/) {
	#        my $diff = substr"$_",0,18;
	#   if ($diff =~  /\'/) {
	#	   $flag_version = V26;
	#	   print "$diff\n";
	#   } else {
	#   $flag_version = 3;
	#   print "$diff V23\n";
	#   }
	# }
        $string =~ s/\'//g; 
        my @mydata = split /\s+/, $string; 
        if ($mydata[0] eq "") { shift(@mydata); } 
        #print "Data parsed correctly...\n"; 
        return @mydata; 
    } 
}
# -----------------------------------------------------------------------
# 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;
}

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -