Docsity
Docsity

Prepare-se para as provas
Prepare-se para as provas

Estude fácil! Tem muito documento disponível na Docsity


Ganhe pontos para baixar
Ganhe pontos para baixar

Ganhe pontos ajudando outros esrudantes ou compre um plano Premium


Guias e Dicas
Guias e Dicas

Solução de Sistemas Lineares por Métodos Iterativos, Resumos de Engenharia Elétrica

A solução de sistemas lineares por métodos iterativos, como o método de gauss-seidel e o método de newton-raphson, com base em um exemplo de engenharia elétrica. O processo iterativo é detalhado passo a passo, demonstrando a convergência dos métodos e a solução final do sistema linear.

Tipologia: Resumos

2022

Compartilhado em 22/11/2022

eduardo-henrique-cabral
eduardo-henrique-cabral 🇧🇷

2 documentos

1 / 71

Toggle sidebar

Esta página não é visível na pré-visualização

Não perca as partes importantes!

bg1
Centro Universitário Avantis
Curso de Graduação em Engenharia Elétrica
Sistemas de Potência
Métodos iterativos de solução de fluxo de potência
Parte III
Prof. Victor Freitas, Dr.
victor.freitas@uniavan.edu.br
Prof. Victor Freitas (UniAvan) Módulo 8 2022/2 1/ 30
pf3
pf4
pf5
pf8
pf9
pfa
pfd
pfe
pff
pf12
pf13
pf14
pf15
pf16
pf17
pf18
pf19
pf1a
pf1b
pf1c
pf1d
pf1e
pf1f
pf20
pf21
pf22
pf23
pf24
pf25
pf26
pf27
pf28
pf29
pf2a
pf2b
pf2c
pf2d
pf2e
pf2f
pf30
pf31
pf32
pf33
pf34
pf35
pf36
pf37
pf38
pf39
pf3a
pf3b
pf3c
pf3d
pf3e
pf3f
pf40
pf41
pf42
pf43
pf44
pf45
pf46
pf47

Pré-visualização parcial do texto

Baixe Solução de Sistemas Lineares por Métodos Iterativos e outras Resumos em PDF para Engenharia Elétrica, somente na Docsity!

Centro Universitário Avantis Curso de Graduação em Engenharia Elétrica

Sistemas de Potência

Métodos iterativos de solução de fluxo de potência

Parte III

Prof. Victor Freitas, Dr. victor.freitas@uniavan.edu.br

Sumário

(^1) Expressões Gerais dos Fluxos de Potência

(^2) Métodos de Solução do Problema

(^3) Solução pelo Método Gauss/Gauss-Seidel

(^4) Métodos Baseados na Matriz Y

(^5) Solução pelo Método de Newton-Raphson

Expressões Gerais dos Fluxos de Potência

Equações básicas são obtidas impondo-se a conservação das potências ativa e reativa em cada nó; Isto equivale a se impor a Lei de Kirchhoff das Correntes (LKC) aos nós do sistema. Equações para linha de transmissão

a) Quais as expressões de fluxo de potência ativa e reativa entre as barras k e m? b) Como podem ser calculadas as perdas ativas e reativas da linha;

Expressões Gerais dos Fluxos de Potência

Equações básicas são obtidas impondo-se a conservação das potências ativa e reativa em cada nó; Isto equivale a se impor a Lei de Kirchhoff das Correntes (LKC) aos nós do sistema. Equações para linha de transmissão

a) Quais as expressões de fluxo de potência ativa e reativa entre as barras k e m? b) Como podem ser calculadas as perdas ativas e reativas da linha;

Formulação básica do problema

Para que o problema seja determinístico:

Duas variáveis são conhecidas (dadas) e as outras duas devem ser calculadas (incógnitas). Dependendo de quais são dadas e de quais são incógnitas, são definidos três tipos de barras:

Tipo Dados Incógnitas Características PQ Pk , Qk Vk , θk Barras de carga PV Pk , Vk Qk , θk Barras de Geração incluindo compensadores síncronos

Referência Vk , θk Pk , Qk

Barra de geração (geral- mente uma unidade geradora de grande capacidade (V θ, slack, swing)

Formulação básica do problema

Aplicação de lei das corrente de Kirchhoff para cada barra corresponde ao balanço de potência na barra:

São realizados os balanços das potências ativa e reativa

Pk = ∑

m∈Ωk

Pkm(Vk , Vm, θk , θm)

Qk + Qshk (Vk ) = ∑

m∈Ωk

Qkm(Vk , Vm, θk , θm)

Métodos de Solução do Problema

Ao longo dos anos, vários métodos de solução do FC foram propostos;

Para cada aplicação existem os métodos mais apropriados;

Os fatores considerados na escolha são: Tipos de Solução Precisa Aproximada Sem controle de limites Com controle de limites Off-line On-line Caso simples Casos múltiplos

Métodos de Solução do Problema

Ao longo dos anos, vários métodos de solução do FC foram propostos;

Para cada aplicação existem os métodos mais apropriados;

Os fatores considerados na escolha são: Tipos de Solução Precisa Aproximada Sem controle de limites Com controle de limites Off-line On-line Caso simples Casos múltiplos

Métodos de Solução do Problema

1 Métodos iterativos baseados na matriz Ybarra Baseados na resolução de sistemas de equações lineares I = Y · E; Ex: Gauss/Gauss-Seidel, Glimn-Stagg, Ward-Hale, etc.

2 Método iterativo de Newton Jacobiano completo Métodos desacoplados

Métodos de Solução do Problema

1 Métodos iterativos baseados na matriz Ybarra Baseados na resolução de sistemas de equações lineares I = Y · E; Ex: Gauss/Gauss-Seidel, Glimn-Stagg, Ward-Hale, etc.

2 Método iterativo de Newton Jacobiano completo Métodos desacoplados

Solução pelo Método Gauss/Gauss-Seidel

Gauss

Considere o sistema de equações algébricas Ax = b   

10 x 1 + 2 x 2 + x 3 = − 7 x 1 + 5 x 2 + x 3 = − 8 2 x 1 + 3 x 2 + 10 x 3 = − 6

O sistema pode ser reescrito em termos de x 1 , x 2 e x 3 :

x 1 = (− 7 − 2 x 2 − x 3 )/ 10 x 2 = (− 8 − x 1 − x 3 )/ 5 x 3 = (− 6 − 2 x 1 − 3 x 2 )/ 10

Solução pelo Método Gauss/Gauss-Seidel

Gauss

Atribuindo-se um valor inicial x = [0 0 0]T^ , o seguinte processo iterativo pode ser estabelecido: Iteração 1:

x 1 = (− 7 − 2 x 2 − x 3 )/ 10 = (− 7 − 2 x 0 − 0 )/ 10 = − 0 , 7 x 2 = (− 8 − x 1 − x 3 )/ 5 = (− 8 − 0 − 0 )/ 5 = − 1 , 6 x 3 = (− 6 − 2 x 1 − 3 x 2 )/ 10 = (− 6 − 2 x 0 − 3 x 0 )/ 10 = − 0 , 6

Iteração 2 (com x = [− 0 , 7 − 1 , 6 − 0 , 6 ]T^ ):

x 1 = (− 7 − 2 x(− 1 , 6 ) − (− 0 , 6 ))/ 10 = − 0 , 32 x 2 = (− 8 − (− 0 , 7 ) − (− 0 , 6 ))/ 5 = − 1 , 34 x 3 = (− 6 − 2 x(− 0. 7 ) − 3 x(− 1 , 6 ))/ 10 = 0 , 02

Solução pelo Método Gauss/Gauss-Seidel

Gauss

Atribuindo-se um valor inicial x = [0 0 0]T^ , o seguinte processo iterativo pode ser estabelecido: Iteração 1:

x 1 = (− 7 − 2 x 2 − x 3 )/ 10 = (− 7 − 2 x 0 − 0 )/ 10 = − 0 , 7 x 2 = (− 8 − x 1 − x 3 )/ 5 = (− 8 − 0 − 0 )/ 5 = − 1 , 6 x 3 = (− 6 − 2 x 1 − 3 x 2 )/ 10 = (− 6 − 2 x 0 − 3 x 0 )/ 10 = − 0 , 6

Iteração 2 (com x = [− 0 , 7 − 1 , 6 − 0 , 6 ]T^ ):

x 1 = (− 7 − 2 x(− 1 , 6 ) − (− 0 , 6 ))/ 10 = − 0 , 32 x 2 = (− 8 − (− 0 , 7 ) − (− 0 , 6 ))/ 5 = − 1 , 34 x 3 = (− 6 − 2 x(− 0. 7 ) − 3 x(− 1 , 6 ))/ 10 = 0 , 02

Solução pelo Método Gauss/Gauss-Seidel

Gauss

Continuando o processo iterativo, temos: Iteração 6: x = [− 0 , 3894 − 1 , 5048 − 0 , 0663 ]T

Iteração 7: x = [− 0 , 3894 − 1 , 5048 − 0 , 0707 ]T

Critério de Parada : Solução com erro máximo (Norma Infinita ou L 1 )

ε = max |xυ+^1 − xυ| < 0 , 01

|xυ+^1 − xυ|T^ = [ 0. 0031 0. 0041 0. 0044 ]

x = [− 0 , 3894 − 1 , 5048 − 0 , 0707 ]T