SWITCHING SURGE PROGRAM¶

This program was originally written in Fortran in 1989 by David R. Brown and the American Public Power Association.

The following sources provided data on the electrical properties of different conductors:

AL/ACSR data: Aluminum Electrical Conductor Handbook, 1982.
Alumoweld data: Copperweld Steel Company, Southern Division, Fayetteville, TN.
Steel data: Mechanical Data from ASTM Standards, ASTM Standard A475.
CU data: Westinghouse T&D Reference, 1964, ABB Power Systems Inc.

IF SHIELDS ARE SEGMENTED (INSULATED) TO PREVENT 60 HZ SHIELD WIRE CURRENTS, SHIELD WIRE CURRENTS SHOULD BE SET TO ZERO.

{nins : shield currents ? TO INCLUDE SHIELD CURRENTS, TYPE A 1. TO SET SHIELD CURRENTS TO ZERO, TYPE A 2.}

{xlngth: LINE LENGTH IN MILES}

THE SWITCHING SURGE PROGRAM SIMULATES LINE RECLOSING ASSUMING TRAPPED CHARGE ON THE LINE. THE SYSTEM TO WHICH THE LINE WILL BE CONNECTED IS MODELLED ON THE BASIS OF THE SYSTEM THREE-PHASE FAULT CURRENT.

{fltcurr: THREE-PHASE FAULT CURRENT FOR THE SYSTEM IN AMPERES }

{nres : SIMULATION OF INSERTION RESISTORS ? TO USE INSERTION RESISTORS TYPE A 1. TO OMIT RESISTORS TYPE A 2.}

Omit insertion resistance and resistor insertion time if you chose to omit resistors

{rins: insertion resistance(OHMS)}

{tins: resistor insertion time (ms) ? INSERTION TIME IS THE TIME THAT THE INSERTION RESISTORS ARE IN THE CIRCUIT.}

In [1]:

# Click the Blue Plane to preview this as a CrossCompute tool
conductor_table_path = 'conductors.csv'

phase_table_path = 'phases.csv'
shield_table_path = 'shields.csv'
span_table_path = 'spans.csv'

nins = 1
nres = 1
rins = 400
tins = 10
xlngth = 70
fltcurr = 10

target_folder = '/tmp'

In [2]:

nins = 1 if int(nins) == 1 else 2
nres = 1 if int(nres) == 1 else 2

In [3]:

from os.path import join
output_text_path = join(target_folder, 'ssrg.log')
target_text_path = join(target_folder, 'ssrg.txt')

In [4]:

from macros import transmogrify_conductors
conductor_text_path = transmogrify_conductors(join(
    target_folder, 'conductors.dat'), conductor_table_path)

In [5]:

from macros import transmogrify_pylons
pylon_text_path = transmogrify_pylons(join(
    target_folder, 'pylons.dat',
), phase_table_path, shield_table_path, span_table_path)

In [6]:

from subprocess import run

executable = 'ssrg.out'
script = 'ssrg.f'

output_text = run([
    'gfortran',
    '-o',
    executable,
    script,])

In [7]:

from os.path import expanduser
from subprocess import check_output, STDOUT
output_text = check_output([
    join('.', executable),
    conductor_text_path,
    pylon_text_path,
    expanduser(target_text_path),
    str(nins),
    str(nres),
    str(rins),
    str(tins),
    str(xlngth),
    str(fltcurr),
], stderr=STDOUT)
with open(output_text_path, 'wb') as output_file:
    output_file.write(output_text)

In [8]:

print('output_text_path = ' + output_text_path)
print('target_text_path = ' + target_text_path)

output_text_path = /tmp/ssrg.log
target_text_path = /tmp/ssrg.txt

Pay Notebook Creator: Salah Ahmed	0
Set Container: Numerical CPU with TINY Memory for 10 Minutes	0
Total	0

American Public Power Association tools

SWITCHING SURGE PROGRAM¶