Trabajo final: curso de fortran

Instituto Politécnico Nacional
SEPI – ESIME
Unidad Zacatenco

Curso de FORTRAN - Proyecto Final

Tema: Análisis de Flujos de Potencia por el Método de Gauss Seidel

Elaboro: Gladis Pacheco Bautista
Índice

Análisis De Flujos De Potencia Por El Método De Gauss Seidel 2

Algoritmo Computacional 2

Diagramas de Rutinas 3

Codificación 4

Simulación 14

Conclusión 19

AnálisisDe Flujos De Potencia Por El Método De Gauss Seidel

Los estudios de flujo de potencia, son sumamente importantes para evaluar el funcionamiento de los sistemas de potencia, su control y planificación para expansiones futuras. Un estudio de flujo de potencia define principalmente las potencias activa y reactiva y el vector de tensión en cada bus en el sistema.

Entonces la necesidadsistemática de cálculos numéricos requiere que se ejecuten por medio de un procedimiento iterativo; dos de los normalmente más usados son el método Gauss-Seidel y el método Newton-Raphson.

El método de Gauss-Seidel es un procedimiento iterativo para resolver simultáneamente ecuaciones no lineales, los valores obtenidos son utilizados inmediatamente después de haber sido calculados aunque no haya terminadola iteración en curso para que el sistema converja con mayor rapidez.

Algoritmo Computacional

1.- Método de Análisis nodal es el más empleado.

2.- Construir la Matriz de Admitancias de barras: YBUS.

3.- Ajustarse a los tipos de buses (carga, tensión controlada, referencia)

4.- Modelización de las Líneas mediante su esquema en π.

5.- Calculo de Voltajes por Gauss Seidel

6.-Calculo de Flujos de Potencia

Diagramas de Rutinas

El proyecto llamado FLUJOSGS consta de un programa principal del mismo nombre y de 6 subrutinas que se describen a continuación.

[pic]
Diagrama De Bloques Del Proyecto FLUJOSGS

FLUJOSGS Programa principal, desde este se llaman las subrutinas.

PRESENTA Es la “Portada” del proyecto, y también designa el archivo de datos y el deresultados.

LEEDATOS Lee los datos del sistema que se va a analizar.

MATRIZYBUS Construye la matriz YBUS por inspección de la red.

PAUX Calcula los parámetros auxiliares KLp & YLpq para el cálculo de voltajes nodales.

VOLT Calcula los voltajes nodales por el método de Gauss-Seidel.

POTENCIA Calcula los flujos de potencia en las línea, perdida total y balance reactivo del sistema.Codificación

Consideraciones Básicas

Los elementos básicos del programa consisten en:

- Especificar las variables a utilizar.
- Apertura de archivos de datos y de resultados.
- Lectura de datos de conectividad de elementos e impedancias primitivas.
- Determinación de los componentes de la matriz YBUS.
- Cálculo de parámetros auxiliares.
- Cálculo devoltajes nodales.
- Cálculo de flujos de potencia y pérdidas.
- Vaciado de resultados en archivo correspondiente.

Programa Principal FLUJOSGS.F90

PROGRAM FLUJOSGS

IMPLICIT NONE

! Llamado de subrutinas
CALL PRESENTA
CALL LEEDATOS
CALL MATRIZYBUS
CALL PAUX
CALL VOLT
CALL POTENCIA

PAUSE

PRINT *, '>> TERMINO FLUJOS'
END PROGRAM FLUJOSGSSub-rutina PRESENTA.F90

SUBROUTINE PRESENTA

IMPLICIT NONE

! ZONA DE DECLARACION DE VARIABLES: REAL, INTEGER, COMPLEX, CHARACTER,...
CHARACTER(len=30) ::NOMBREALUMNO
CHARACTER (len=25) :: archdat,archres

PRINT *,'PROGRAMA2: FLUJOS POR GAUSS SEIDEL'
PRINT *,'TECLEA TU NOMBRE:'

READ(5,5) NOMBREALUMNO
5 FORMAT(A30)

! SALIDA CON FORMATO
WRITE(6,10) NOMBREALUMNO10 FORMAT(3/,5X,'PROGRAMA2: FLUJOS POR GAUSS SEIDEL',2/,5X, 'ALUMNO:', & A30,/,5X,'SEPI ESIME IPN')

PRINT *,'DAME EL NOMBRE DEL ARCHIVO DE DATOS: '
READ(5,15) archdat

15 FORMAT(A25)

OPEN(UNIT=7,FILE=archdat)

PRINT *,'DAME EL NOMBRE DEL ARCHIVO DE RESULTADOS: '
READ(5,15) archres

OPEN(UNIT=8,FILE=archres)...