Saltar al contenido principal

# 5.3: Técnica Iterativa Gauss-Seidel

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$

( \newcommand{\kernel}{\mathrm{null}\,}\) $$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\id}{\mathrm{id}}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\kernel}{\mathrm{null}\,}$$

$$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$

$$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$

$$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\AA}{\unicode[.8,0]{x212B}}$$

$$\newcommand{\vectorA}[1]{\vec{#1}} % arrow$$

$$\newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow$$

$$\newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vectorC}[1]{\textbf{#1}}$$

$$\newcommand{\vectorD}[1]{\overrightarrow{#1}}$$

$$\newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}}$$

$$\newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}}$$

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\avec}{\mathbf a}$$ $$\newcommand{\bvec}{\mathbf b}$$ $$\newcommand{\cvec}{\mathbf c}$$ $$\newcommand{\dvec}{\mathbf d}$$ $$\newcommand{\dtil}{\widetilde{\mathbf d}}$$ $$\newcommand{\evec}{\mathbf e}$$ $$\newcommand{\fvec}{\mathbf f}$$ $$\newcommand{\nvec}{\mathbf n}$$ $$\newcommand{\pvec}{\mathbf p}$$ $$\newcommand{\qvec}{\mathbf q}$$ $$\newcommand{\svec}{\mathbf s}$$ $$\newcommand{\tvec}{\mathbf t}$$ $$\newcommand{\uvec}{\mathbf u}$$ $$\newcommand{\vvec}{\mathbf v}$$ $$\newcommand{\wvec}{\mathbf w}$$ $$\newcommand{\xvec}{\mathbf x}$$ $$\newcommand{\yvec}{\mathbf y}$$ $$\newcommand{\zvec}{\mathbf z}$$ $$\newcommand{\rvec}{\mathbf r}$$ $$\newcommand{\mvec}{\mathbf m}$$ $$\newcommand{\zerovec}{\mathbf 0}$$ $$\newcommand{\onevec}{\mathbf 1}$$ $$\newcommand{\real}{\mathbb R}$$ $$\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}$$ $$\newcommand{\laspan}[1]{\text{Span}\{#1\}}$$ $$\newcommand{\bcal}{\cal B}$$ $$\newcommand{\ccal}{\cal C}$$ $$\newcommand{\scal}{\cal S}$$ $$\newcommand{\wcal}{\cal W}$$ $$\newcommand{\ecal}{\cal E}$$ $$\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}$$ $$\newcommand{\gray}[1]{\color{gray}{#1}}$$ $$\newcommand{\lgray}[1]{\color{lightgray}{#1}}$$ $$\newcommand{\rank}{\operatorname{rank}}$$ $$\newcommand{\row}{\text{Row}}$$ $$\newcommand{\col}{\text{Col}}$$ $$\renewcommand{\row}{\text{Row}}$$ $$\newcommand{\nul}{\text{Nul}}$$ $$\newcommand{\var}{\text{Var}}$$ $$\newcommand{\corr}{\text{corr}}$$ $$\newcommand{\len}[1]{\left|#1\right|}$$ $$\newcommand{\bbar}{\overline{\bvec}}$$ $$\newcommand{\bhat}{\widehat{\bvec}}$$ $$\newcommand{\bperp}{\bvec^\perp}$$ $$\newcommand{\xhat}{\widehat{\xvec}}$$ $$\newcommand{\vhat}{\widehat{\vvec}}$$ $$\newcommand{\uhat}{\widehat{\uvec}}$$ $$\newcommand{\what}{\widehat{\wvec}}$$ $$\newcommand{\Sighat}{\widehat{\Sigma}}$$ $$\newcommand{\lt}{<}$$ $$\newcommand{\gt}{>}$$ $$\newcommand{\amp}{&}$$ $$\definecolor{fillinmathshade}{gray}{0.9}$$

Esta es una de las muchas técnicas para resolver el problema del flujo de carga no lineal. Cabe señalar que esta técnica de solución, aunque sencilla de usar y fácil de entender, tiene una tendencia a utilizar mucho cálculo, particularmente en el trabajo de grandes problemas. También es bastante capaz de converger en soluciones incorrectas (eso es un problema con los sistemas no lineales). Al igual que con otras técnicas iterativas, a menudo es difícil saber cuándo se ha alcanzado la solución correcta. A pesar de estas deficiencias, Gauss-Seidel se puede utilizar para tener una buena sensación de problemas de flujo de carga sin exceso de equipaje de análisis numérico.

Supongamos que tenemos una estimación inicial (ok: guess) para voltajes de red. Podemos particionar la ecuación\ ref {8} como:

$\ \mathbf{S}_{k}=\mathbf{V}_{k} \sum_{j \neq k} \mathbf{Y}_{j k}^{*} \mathbf{V}_{j}^{*}+\mathbf{V}_{k} \mathbf{Y}_{k k}^{*} \mathbf{V}_{k}^{*}\label{9}$

Señalando que$$\ \mathbf{S}_{k}=P_{k}+j Q_{k}$$, podemos resolver por$$\ \mathbf{V}_{k}^{*}$$ y, tomando el complejo conjugado de eso, tenemos una expresión para$$\ \mathbf{V}_{k}$$ en términos de todos los voltajes de la red,$$\ P_{k}$$ y$$\ Q_{k}$$:

$\ \mathbf{V}_{k}=\frac{1}{\mathbf{Y}_{k k}}\left(\frac{P_{k}-j Q_{k}}{\mathbf{V}_{k}^{*}}-\sum_{j \neq k} \mathbf{Y}_{j k} \mathbf{V}_{j}\right)\label{10}$

La ecuación de expresión\ ref {10} es una mejor estimación de la$$\ \mathbf{V}_{k}$$ que comenzamos con. La solución al conjunto de ecuaciones no lineales consiste en llevar a cabo esta expresión, repetidamente, para todos los buses de la red.

Un procedimiento iterativo en el que se calcula una corrección para cada una de las tensiones de la red en un solo paso, y las correcciones aplicadas todas a la vez se denomina Iteración Gaussiana. Si, por otro lado, las variables mejoradas se utilizan inmediatamente, el procedimiento se denomina Iteración Gauss—Seidel.

Tenga en cuenta que la Ecuación\ ref {10} utiliza como sus restricciones$$P$$ y$$Q$$ para el bus en cuestión. Por lo tanto, se puede usar directamente para buses tipo carga. Para otros tipos de restricciones de bus, se requieren modificaciones. Consideramos solo dos de los muchos conjuntos posibles de restricciones.

Para los buses de generadores, generalmente se especifican la potencia real y la magnitud del voltaje del terminal. En cada paso de tiempo es necesario salir con un voltaje terminal de magnitud especificada: el ángulo de fase del voltaje y la potencia reactiva$$Q$$ son las incógnitas. Una forma de manejar esta situación es:

1. Generar una estimación para la potencia reactiva$$Q$$, luego
2. Utilice la ecuación\ ref {10} para generar una estimación para el voltaje del terminal, y finalmente,
3. Al mantener constante el ángulo de fase del voltaje, ajustar la magnitud a la restricción.

En cualquier punto de la iteración, la potencia reactiva es:

$\ Q_{k}=\operatorname{Im}\left\{\mathbf{V}_{k} \sum_{j=1}^{N} \mathbf{Y}_{j k}^{*} \mathbf{V}_{j}^{*}\right\}\label{11}$

Cabe señalar que existen otras formas de hacer este cálculo. Generalmente son más trabajo de configurar pero a menudo convergen más rápidamente. El método y las variaciones de Newton son buenos ejemplos.

Para los buses cargados por impedancia constante, es suficiente agregar la impedancia de carga a la red. Es decir, la admisión de carga va directamente en paralelo con la admitancia del punto de conducción en el nodo en cuestión.

Estos tres tipos de restricción de bus, generador, carga e impedancia constante son suficientes para manejar la mayoría de los problemas de importancia práctica.

This page titled 5.3: Técnica Iterativa Gauss-Seidel is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by James Kirtley (MIT OpenCourseWare) via source content that was edited to the style and standards of the LibreTexts platform.