🌐 Aula 06 — Aplicações do Pensamento Computacional em Engenharia Elétrica

Navegação Rápida

📍 Disciplina: 02 Pensamento Computacional | 🗺️ Área: 06 Computação e Telecomunicações | 🏠 Hub Principal ⬅️ Aula Anterior: Aula 05 - Representação Digital da Informação | 🔝 MOC da Disciplina: Hub 02 Pensamento Computacional

---

📌 Sumário da Aula

1. 🧠 1. Integração Holística dos Quatro Pilares na Engenharia Elétrica
2. ⚡ 2. Redes Elétricas Inteligentes (Smart Grids) e Automação Self-Healing
3. 📐 3. Formulação Matemática do Problema FLISR e Restrições Físicas
4. 🏭 4. Despacho Econômico de Geração e Multiplicadores de Lagrange
5. 🔮 5. Gêmeos Digitais (Digital Twins) e Computação de Borda em Subestações
6. 📊 6. Diagrama de Arquitetura Self-Healing em Smart Grids
7. 📝 7. Questões de Concursos e Problemas Avançados Resolvidos
8. 🔗 8. Conexões e Referências Cruzadas

---

🧠 1. Integração Holística dos Quatro Pilares na Engenharia Elétrica

Ao longo das cinco primeiras aulas, foram dissecados os princípios basilares da ciência computacional: a Decomposição (Aula 01), o Reconhecimento de Padrões (Aula 02), a Abstração (Aula 03), o Projeto de Algoritmos com Análise Assintótica (Aula 04) e a Representação Digital Matricial da Informação (Aula 05).

Na Engenharia Elétrica contemporânea, esses conceitos deixam de ser tópicos teóricos isolados para formar o alicerce operacional dos Sistemas Ciber-Físicos (Cyber-Physical Systems - CPS). A transição da rede elétrica tradicional unidirecional (geração centralizada fóssil/hidrelétrica $\to$ transmissão $\to$ distribuição passiva $\to$ consumidor final) para as Redes Elétricas Inteligentes (Smart Grids) é orientada pela convergência estrita da eletrotécnica pesada com a ciência da computação.

flowchart TD
    Desafio["Desafio de Engenharia:<br/>Operação em Tempo Real de Redes Ativas Complexas"] --> P1["1. Decomposição:<br/>Particionamento em Microrredes e Zonas de Proteção"]
    Desafio --> P2["2. Reconhecimento de Padrões:<br/>Telemetria AMI, IA Preditiva de Carga e Geração Solar/Eólica"]
    Desafio --> P3["3. Abstração:<br/>Modelagem em Grafos Topológicos e Matrizes Nodal Ybus / Jacobiano"]
    Desafio --> P4["4. Projeto de Algoritmos:<br/>Otimização Restrita, FLISR, Despacho Econômico e Self-Healing"]

    P1 --> Operacao["Operação Ciber-Física Segura, Resiliente e de Mínimo Custo"]
    P2 --> Operacao
    P3 --> Operacao
    P4 --> Operacao

---

⚡ 2. Redes Elétricas Inteligentes (Smart Grids) e Automação Self-Healing

Uma das aplicações mais sofisticadas do pensamento computacional reside nos sistemas de Autorrecuperação de Redes de Distribuição (Self-Healing). Em alimentadores aéreos urbanos e rurais de média tensão ($13,8\text{ kV}$ a $34,5\text{ kV}$), descargas atmosféricas, quedas de árvores ou toques de animais provocam curtos-circuitos que desarmam os disjuntores da subestação, interrompendo o fornecimento de energia para milhares de consumidores.

Historicamente, a localização e manobra da falta exigiam o envio de viaturas com eletricistas para inspeção visual ao longo de quilômetros de rede, resultando em índices regulatórios de qualidade de serviço degradados (DEC - Duração Equivalente de Interrupção por Unidade Consumidora e FEC - Frequência Equivalente de Interrupção por Unidade Consumidora).

O sistema FLISR (Fault Location, Isolation, and Service Restoration) automatiza completamente esse processo em menos de $60\text{ segundos}$ mediante a execução coordenada dos quatro pilares:

sequenceDiagram
    autonumber
    participant Falta as Curto-Circuito Físico
    participant IED as Relés Religadores (IED)
    participant SCADA as Algoritmo FLISR Central / Borda
    participant Chaves as Chaves Telecontroladas Motorizadas

    Falta->>IED: Corrente de Falta Ik > 5000 A
    IED->>SCADA: Telemetria Digital de Trip (DNP3 / IEC 61850)
    SCADA->>SCADA: 1. Decomposição: Isolar Trecho Faltoso
    SCADA->>Chaves: Comando Remoto: Abrir Chaves Imediatamente à Montante e Jusante
    Chaves-->>SCADA: Confirmação de Estado Aberto
    SCADA->>SCADA: 2. Otimização: Calcular Caminho de Restauração via Alimentador B
    SCADA->>Chaves: Comando Remoto: Fechar Chave de Fronteira (Tie-Switch)
    Chaves-->>SCADA: Confirmação de Fechamento (Rede Reenergizada)

---

📐 3. Formulação Matemática do Problema FLISR e Restrições Físicas

O problema de Restauração de Serviço pós-falta é formalizado matematicamente como um problema de Programação Linear Inteira Mista (MILP) ou Otimização Combinatória Não Linear com restrições elétricas estritas.

3.1 Função Objetivo

Deseja-se maximizar a quantidade total de potência ativa restabelecida para os blocos de carga desenergizados sadios (${\mathcal{N}}_{\text{sadio}}$), minimizando simultaneamente as perdas elétricas por Efeito Joule nos alimentadores e o número de manobras mecânicas de abertura/fechamento das chaves seccionadoras:

${\max }_{\mathbf{s}}\left[{\sum }_{i\in {\mathcal{N}}_{\text{sadio}}}{w}_i\cdot P_{L,i}\cdot y_i-{\sum }_{k\in \mathcal{E}}{C}_{\text{manobra}}\cdot |\Delta s_k|-\beta {\sum }_{j\in {\mathcal{E}}_{\text{ativo}}}{R}_j\cdot I_j^2\right]$

Onde:

• $y_i\in \{0,1\}$ é a variável binária indicando se a barra de carga $i$ foi restaurada ($1$) ou permanece desenergizada ($0$);
• $w_i$ é o fator de ponderação de prioridade da carga (pesos elevados para hospitais, estações de tratamento de água e telecomunicações);
• $s_k\in \{0,1\}$ é o estado de operação da chave $k$ (aberta ou fechada);
• $|\Delta s_k|=|s_k^{\text{novo}}-s_k^{\text{antigo}}|$ penaliza o desgaste eletromecânico dos contatos da chave;
• $R_j\cdot I_j^2$ computa as perdas dissipativas nos cabos condutores.

3.2 Restrições Físicas e Regulatórias de Operação

1. Restrição de Topologia Radial: A maioria esmagadora das redes de distribuição urbanas deve operar estritamente em topologia radial (estrutura em árvore conexa sem ciclos fechados), para evitar correntes de circulação descontroladas e manter a simplicidade da coordenação de relés de sobrecorrente direcionais: $|{\mathcal{E}}_{\text{ativo}}|=|{\mathcal{V}}_{\text{ativo}}|-1\text{e}\text{Grafo eˊ Conexo}$

2. Capacidade Térmica Máxima dos Condutores (Limite Ampacity): A corrente elétrica em regime permanente fluindo em cada ramo $j$ não pode ultrapassar o limite térmico de aquecimento do cabo condutor estabelecido pela norma ABNT NBR 5422: $I_j\le I_{j,\text{maˊx}}^{\text{teˊrmico}},\forall j\in {\mathcal{E}}_{\text{ativo}}$

3. Limites Regulatórios de Queda de Tensão Nodal (PRODIST Módulo 8 da ANEEL): O perfil de tensão em regime permanente em qualquer barra consumidora $k$ deve permanecer rigorosamente dentro da faixa de tensão adequada por unidade ($\text{pu}$): $\boxed{0,93\text{ pu}\le |V_k|\le 1,05\text{ pu},\forall k\in \mathcal{V}}$

4. Capacidade Nominal do Transformador da Subestação Adjacente: A potência aparente total injetada no alimentador de socorro após o fechamento da chave de interligação (Tie-Switch) não pode provocar sobrecarga no transformador de força da subestação que assume a nova carga: $S_{\text{SE, socorro}}\le S_{\text{SE, nominal}}$

---

🏭 4. Despacho Econômico de Geração e Multiplicadores de Lagrange

Outra aplicação clássica do projeto de algoritmos e modelagem matemática em sistemas elétricos é o Despacho Econômico de Carga (Economic Load Dispatch). Em um sistema interligado contendo $N_g$ usinas geradoras termoelétricas, o operador nacional do sistema (ONS) deve determinar quanto cada usina $i$ deve gerar ($P_{g,i}$) para atender à demanda total requerida pelos consumidores ($P_D$), com o mínimo custo de combustível por hora.

4.1 Formulação Matemática do Custo de Geração

A curva de custo horário de combustível de uma turbina termoelétrica a gás ou carvão é aproximada com excelente precisão por uma parábola quadrática:

$C_i(P_{g,i})={\alpha }_i\cdot P_{g,i}^2+{\beta }_i\cdot P_{g,i}+{\gamma }_i[\text{R\$/h}]$

Onde ${\alpha }_i>0,{\beta }_i>0,{\gamma }_i>0$ são constantes empíricas da planta. O problema de otimização irrestrita (desprezando perdas na linha para formulação analítica) é:

$\min {\sum }_{i=1}^{N_g}{C}_i(P_{g,i})\text{sujeito a}{\sum }_{i=1}^{N_g}{P}_{g,i}-P_D=0$

4.2 Otimização com Multiplicadores de Lagrange

Monta-se a função lagrangiana $\mathcal{L}$:

$\mathcal{L}(P_{g,1},\dots ,P_{g,N_g},\lambda )={\sum }_{i=1}^{N_g}{C}_i(P_{g,i})-\lambda \left({\sum }_{i=1}^{N_g}{P}_{g,i}-P_D\right)$

As condições necessárias de primeira ordem para o ponto ótimo são:

$\frac{\partial \mathcal{L}}{\partial P_{g,i}}=\frac{d{C}_i}{d{P}_{g,i}}-\lambda =0⟹\boxed{\frac{d{C}_i}{d{P}_{g,i}}=\lambda ,\forall i\in \{1,\dots ,N_g\}}$ $\frac{\partial \mathcal{L}}{\partial \lambda }={\sum }_{i=1}^{N_g}{P}_{g,i}-P_D=0$

Onde $\frac{d{C}_i}{d{P}_{g,i}}=2{\alpha }_i{P}_{g,i}+{\beta }_i$ representa o Custo Marginal de Produção da usina $i$.

Teorema Fundamental do Despacho Econômico: A condição de custo global mínimo é alcançada quando todas as usinas em operação operam com exatamente o mesmo custo marginal de geração $\lambda$.

Isolando $P_{g,i}$ em função do multiplicador de Lagrange $\lambda$: $P_{g,i}=\frac{\lambda -{\beta }_i}{2{\alpha }_i}$

Substituindo na equação de balanço de potência: ${\sum }_{i=1}^{N_g}\left(\frac{\lambda -{\beta }_i}{2{\alpha }_i}\right)=P_D⟹\boxed{\lambda =\frac{P_D+{\sum }_{i=1}^{N_g}\frac{{\beta }_i}{2{\alpha }_i}}{{\sum }_{i=1}^{N_g}\frac{1}{2{\alpha }_i}}}$

---

🔮 5. Gêmeos Digitais (Digital Twins) e Computação de Borda em Subestações

A convergência dos 4 pilares atinge o estado da arte na criação de Gêmeos Digitais (Digital Twins) de transformadores de potência de subestações de transmissão ($500\text{ kV}/230\text{ kV}$):

flowchart LR
    TransfReal["Transformador Físico de Potência (300 MVA)<br/>Sensores Ópticos de Fibras no Enrolamento,<br/>Cromatografia de Gases Dissolvidos (DGA)"] -->|Protocolo IEC 61850 SV / GOOSE| Gateway["Edge Computer em Tempo Real (Subestação)<br/>Amostragem Fasorial Sincronizada PMU"]
    Gateway -->|Modelagem Multifísica Discreta| DigitalTwin["Gêmeo Digital (Digital Twin)<br/>Modelo Térmico IEEE C57.91 + Degradação de Papel"]
    DigitalTwin --> Alerta["Predição de Falhas Dielétricas<br/>e Estimativa de Vida Útil Restante (RUL)"]

5.1 Modelo Térmico da Temperatura do Ponto Mais Quente (Hotspot)

O envelhecimento dielétrico do isolamento em papel celulósico kraft imerso em óleo mineral isolante é governado pela temperatura do ponto mais quente do enrolamento (${\theta }_{hs}$). Conforme a norma IEEE C57.91, a elevação de temperatura em regime permanente é calculada em tempo real por equações diferenciais discretizadas executando em microprocessadores na borda:

${\theta }_{hs}={\theta }_{\text{ambiente}}+\Delta {\theta }_{to}+\Delta {\theta }_{hs,to}$

Onde:

• $\Delta {\theta }_{to}$ é a elevação de temperatura do topo do óleo;
• $\Delta {\theta }_{hs,to}$ é a elevação do enrolamento em relação ao óleo.

A taxa de perda acelerada de vida útil do isolamento ($V$) segue a Equação de Arrhenius:

$V=\exp \left(\frac{15000}{110+273,15}-\frac{15000}{{\theta }_{hs}+273,15}\right)$

Para uma temperatura de referência de ${\theta }_{hs}=110^{\circ }\text{C}$, $V=1,0$ (envelhecimento normal). Se a carga elevar a temperatura para ${\theta }_{hs}=122^{\circ }\text{C}$, a taxa salta para $V\approx 2,5$, significando que o transformador envelhece $2,5\text{ dias}$ para cada $24\text{ horas}$ de operação sob essa condição de estresse térmico.

---

📊 6. Diagrama de Arquitetura Self-Healing em Smart Grids

O diagrama abaixo, desenvolvido na estética Catppuccin Mocha, exibe o fluxo algorítmico do Self-Healing para recomposição automática de redes de distribuição:

---

📝 7. Questões de Concursos e Problemas Avançados Resolvidos

Questão 1 (Banca CESGRANRIO - Engenheiro de Eletricidade / Eletrobras)

Duas usinas termoelétricas a gás natural operam em regime de despacho econômico conjunto para alimentar uma carga total demandada de $P_D=800\text{ MW}$. As funções de custo horário de combustível de cada usina são dadas por: $C_1(P_1)=0,004P_1^2+8P_1+1000[\text{R\$/h}]$ $C_2(P_2)=0,006P_2^2+6P_2+1200[\text{R\$/h}]$ Desprezando as perdas elétricas na transmissão e assumindo que nenhuma usina atinge os seus limites mínimos ou máximos de capacidade de geração, determine:

1. O custo marginal de geração ótimo do sistema $\lambda$ em \text{R\/MWh}$.
2. A potência ativa despachada em cada usina ($P_1$ e $P_2$).
3. O custo financeiro total de geração do sistema por hora.

Resolução Passo a Passo:

1. Cálculo dos Custos Marginais Individuais ($\frac{d{C}_i}{d{P}_i}$): $\frac{d{C}_1}{d{P}_1}=2\times 0,004P_1+8=0,008P_1+8$ $\frac{d{C}_2}{d{P}_2}=2\times 0,006P_2+6=0,012P_2+6$

2. Aplicação do Princípio do Despacho Ótimo: $\lambda =0,008P_1+8⟹P_1=\frac{\lambda -8}{0,008}$ $\lambda =0,012P_2+6⟹P_2=\frac{\lambda -6}{0,012}$

3. Balanço Global de Potência Ativa: $P_1+P_2=P_D⟹\frac{\lambda -8}{0,008}+\frac{\lambda -6}{0,012}=800$

   Multiplicando todos os termos por $0,024$ (mínimo múltiplo comum): $3(\lambda -8)+2(\lambda -6)=800\times 0,024$ $3\lambda -24+2\lambda -12=19,2$ $5\lambda -36=19,2⟹5\lambda =55,2$ $\boxed{\lambda =\frac{55,2}{5}=11,04\text{ R\$/MWh}}$

4. Cálculo das Potências Despachadas ($P_1$ e $P_2$): $P_1=\frac{11,04-8}{0,008}=\frac{3,04}{0,008}=380\text{ MW}$ $P_2=\frac{11,04-6}{0,012}=\frac{5,04}{0,012}=420\text{ MW}$

   Verificação: $P_1+P_2=380+420=800\text{ MW}$. Perfeito!

5. Cálculo do Custo Total Horário ($C_{\text{total}}=C_1+C_2$):

   • Para a Usina 1: $C_1(380)=0,004(380{)}^2+8(380)+1000=0,004(144400)+3040+1000$ $C_1(380)=577,60+3040+1000=4617,60\text{ R\$/h}$
   • Para a Usina 2: $C_2(420)=0,006(420{)}^2+6(420)+1200=0,006(176400)+2520+1200$ $C_2(420)=1058,40+2520+1200=4778,40\text{ R\$/h}$

   $C_{\text{total}}=4617,60+4778,40=9396,00\text{ R\$/h}$

$\boxed{\lambda =11,04\text{ R\$/MWh};P_1=380\text{ MW};P_2=420\text{ MW};C_{\text{total}}=9396,00\text{ R\$/h}}$

---

Questão 2 (ENADE Engenharia Elétrica / Concurso Perito Criminal)

Em um alimentador de distribuição urbana com tensão nominal de $13,8\text{ kV}$, um sistema automatizado de Self-Healing isolou com sucesso um trecho faltoso e precisa religar $1200\text{ kVA}$ de carga desenergizada ao alimentador vizinho. Sabe-se que:

• O alimentador vizinho já atende uma carga pré-existente de $3500\text{ kVA}$;
• A capacidade nominal contínua do cabo de saída na subestação é de $300\text{ A}$;
• A tensão de fornecimento de barra é regulada em $1,0\text{ pu}$ ($13,8\text{ kV}$).

O algoritmo de controle do Self-Healing deve autorizar a manobra de fechamento da chave de interligação (Tie-Switch)? Justifique rigorosamente através do cálculo da corrente elétrica resultante.

Resolução Passo a Passo:

1. Cálculo da Carga Total Agregada após a Recomposição: $S_{\text{total}}=S_{\text{preˊ-existente}}+S_{\text{transferida}}=3500\text{ kVA}+1200\text{ kVA}=4700\text{ kVA}$

2. Cálculo da Corrente Trifásica Resultante ($I_{\text{resultante}}$): $I_{\text{resultante}}=\frac{S_{\text{total}}}{\sqrt{3}\times V_{\text{linha}}}=\frac{4700\times 10^3\text{ VA}}{\sqrt{3}\times 13,8\times 10^3\text{ V}}$ $I_{\text{resultante}}=\frac{4700}{1,73205\times 13,8}=\frac{4700}{23,9023}\approx 196,63\text{ A}$

3. Avaliação da Restrição de Ampacidade: $I_{\text{resultante}}=196,63\text{ A}\le I_{\text{maˊx}}^{\text{cabo}}=300\text{ A}$ A corrente solicitada representa apenas $\frac{196,63}{300}\approx 65,54\%$ da capacidade nominal do cabo.

$\boxed{\text{Sim, o algoritmo DEVE autorizar a manobra (Corrente de }196,63\text{ A }<300\text{ A)}}$

Conclusão Operacional: Como a restrição de corrente é plenamente atendida com mais de $34\%$ de margem de segurança e não há risco de sobrecarga térmica no condutor principal, o algoritmo comanda com segurança o pulso de fechamento motorizado da chave de socorro.

---

Questão 3 (Problema Analítico de Engenharia - Envelhecimento Térmico de Transformador)

Um transformador de potência de $50\text{ MVA}$ em uma subestação industrial é submetido a uma sobrecarga de emergência durante 4 horas diárias, elevando a temperatura de seu ponto mais quente para ${\theta }_{hs}=125^{\circ }\text{C}$. Durante as 20 horas restantes do dia, o transformador opera na temperatura nominal de referência ${\theta }_{hs}=110^{\circ }\text{C}$ ($V=1,0$). Calcule a perda relativa diária de vida útil equivalente do transformador utilizando o modelo de Arrhenius da norma IEEE C57.91: $V=\exp \left(\frac{15000}{383,15}-\frac{15000}{{\theta }_{hs}+273,15}\right)$

Resolução Passo a Passo:

1. Cálculo do Fator de Aceleração $V$ para ${\theta }_{hs}=125^{\circ }\text{C}$ ($T=398,15\text{ K}$): $\frac{15000}{383,15}\approx 39,14916$ $\frac{15000}{398,15}\approx 37,67424$ Diferença nos expoentes: $\Delta =39,14916-37,67424=1,47492$

   $V_{125}=\exp (1,47492)\approx 4,3707$

2. Cálculo do Envelhecimento Equivalente Ponderado nas 24 Horas: $V_{\text{diaˊrio, meˊdio}}=\frac{(V_{125}\times 4\text{ h})+(V_{110}\times 20\text{ h})}{24\text{ h}}$ $V_{\text{diaˊrio, meˊdio}}=\frac{(4,3707\times 4)+(1,0\times 20)}{24}=\frac{17,4828+20}{24}=\frac{37,4828}{24}\approx 1,5618$

$\boxed{V_{\text{diaˊrio, meˊdio}}\approx 1,56\text{ (Acelerac¸a˜o de }56,2\%\text{ no consumo de vida uˊtil)}}$

Interpretação: Apenas 4 horas de sobrecarga diária a $125^{\circ }\text{C}$ aceleram o consumo da vida útil do isolamento em mais de $56\%$, fazendo com que um transformador projetado para 30 anos de durabilidade atinja o final de vida em menos de 19 anos se operado continuamente sob esse regime sem controle algorítmico do gêmeo digital.

---

🔗 8. Conexões e Referências Cruzadas

Disciplinas Vinculadas:

• Aula 01: Apresentação do pilar de decomposição.
• Aula 04: Algoritmos em grafos elétricos.
• 05 Distribuição de Energia: Planejamento e operação de redes elétricas urbanas.
• 09 Sistemas de Proteção e Medição: Relés de proteção e esquemas de teleproteção.

Materiais Didáticos Oficiais da Disciplina:

• Arquivo de Curso: 06 - Pensamento Computacional e Aplicações Nas Áreas de Conhecimento.pdf
• Arquivo de Curso: 07 - Praticando.pdf