3.3 Metodos Interactivos Simples

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 4

3.

3 MÉTODOS INTERACTIVOS SIMPLES

Históricamente, debido fundamentalmente a la poca capacidad de cálculo y de memoria de los


computadores primitivos, empezaron usándose métodos que iteraban nudo a nudo sin tener que
manipular grandes matrices, ya fuesen de admitancias [21, 6] o de impedancias [9, 8,3]. Estas
técnicas son aplicaci9ones directas al caso no lineal de los métodos iterativos de solución de
ecuaciones lineales existentes en aquella época, y su interés industrial hoy día es prácticamente
nulo. A pesar de ello, por su simplicidad e interés académico, resumiremos a continuación algunas
delas técnicas más populares que todavía incorporan algunos paquetes comerciales.

3.3.1. MÉTODO DE GAUSS-SEIDEL

Consiste en barrer secuencialmente cada nudo y actualizar su tensión en función delos valores
disponibles en ese momento de todas las tensiones. En general, encontrar el vector x que
satisface el sistema no lineal.

f ( x )=0 (3.19)

Puede formularse como el problema de punto fijo

x=F (x)
Cuya solución, partiendo de un valor inicial x 0, se obtiene iterativamente mediante:

x k+1
i =F i ¿ (3.20)

Obsérvese que cuando se actualiza x i se utilizan los nuevos valores de las variables actualizadas
con anterioridad (i=1,2 , … , i−1).

Para el caso concreto del flujo de cargas, entre las divertidas formas en que puede reescribirse
(3.7), la siguiente parece haberse consolidado como la más eficiente:

1
uki +1= ¿ (3.21)
y ii

i=1,2 , … , n−1
El proceso iterativo se detiene cuando se satisface

max u k+1−u k ≤ ε
i
|i i| (3.22)

Donde ε es un umbral suficientemente pequeño (por ejemplo 0.001).

Aunque el esfuerzo de cálculo por iteración es moderado, la convergencia del método de Gauss-
Seidel es lineal, lo que significa que el número de iteraciones es del orden de n . Esto supone una
limitación importante para sistemas grandes, porque el tiempo de cálculo total crece con n2 . Es
posible disminuir el número de iteraciones, a veces hasta la mitad, mediante un factor de
aceleración α .

uki +1=uki + α (uik+1−uik )


acel
Cuyo valor debe ser menor que 2 para evitar divergencia, los valores óptimos de α están
comprendidos entre 1.4 y 1.6.

Sin embargo, la ecuación (3.21) no puede aplicarse directamente a los nudos PV por los motivos:
esp
1) no se conoce Q i ; 2) el módulo de la tensión resultante no coincide con el especificado. El
método habitualmente empleado para solventar el primer inconveniente consiste en sustituir Q esp
i
por el valor calculado con las tensiones disponibles en cada momento. La segunda limitación se
solventa corrigiendo la tensión obtenida para que tenga el modulo deseado:

uki +1=V esp k +1


corr
i ui /V ik+1

Este mecanismo de corrección no debe llevarse a cabo prematuramente para no perturbar la


convergencia global del proceso, ya de por sí bastante pobre.

Quizás la única utilidad práctica hoy día del método de Gauss-Seidel radique en su utilización como
forma de generar valores iniciales para el método Newton-Raphson, en aquellos raros casos en los
que la convergencia de éste sea problemática partiendo desde el perfil plano.

EJEMPLO 3.2:

Resolveremos la red del ejemplo 3.1 mediante el método de Gauss-Seidel.

Tomando inicialmente Q 3=0 y perfil plano, las tensiones obtenidas después de aplicar dos veces
la ecuación (3.21) a cada uno delos nudos 2 y 3 son:

u2=0.9803−3.369 ∡ ; u3=1.05−2.086 ∡ ;

Los resultados finales obtenidos y el número de iteraciones necesarias se muestran en la siguiente


tabla para dos criterios de convergencia diferentes:

MAX | Δui| 10−4 10−5


Iteraciones 25 33
Nudo V θ V θ
2 0.9931 -9.256 0.9931 -9.262
3 1.05 -5.479 1.05 -5.490
Puede apreciarse la pobre convergencia del método, siendo incluso peor si el nudo 3 fuese de tipo
PQ (33 y 42 iteraciones, respectivamente). En este caso concreto, el uso de factores de aceleración
no aporta ninguna mejora.
Bibliografía:

https://manautomata.files.wordpress.com/2012/10/capitulo5.pdf

https://prezi.com/znotel_t4pyg/soluciones-de-flujo-de-potencia-metodo-gauss-seidel/

https://books.google.com.pe/books?
id=mj13Xk9fR50C&pg=PA126&lpg=PA126&dq=ejercicios+resueltos+de+flujo+de+carga+metodo+g
auss&source=bl&ots=QEFHKkh3i6&sig=WnBuHdOgAIbfh8E4sIGDXexEVqo&hl=es419&sa=X&ved=0
ahUKEwiF5Yaf9MnWAhXBbiYKHV2bA_o4ChDoAQhXMAk#v=onepage&q&f=false

También podría gustarte