Media TensiónIEC 60287AmpacidadElectromagnetismoEfecto Pelicular

Efeito Pelicular e Efeito de Proximidade em Cabos de Média Tensão: Impacto Térmico e Normas

Análise técnica do efeito pelicular e de proximidade em cabos de média tensão. Impacto na ampacidade segundo a norma IEC 60287.

Ing. Francisco Ramírez

Fundamentos Eletrodinâmicos e Teorema de Poynting em Condutores de Média Tensão

A análise rigorosa do comportamento eletromagnético em cabos de média tensão exige o abandono das aproximações circuitais de parâmetros concentrados e o aprofundamento na resolução estrita das Equações de Maxwell em meios materiais lineares, isotrópicos e homogêneos. Quando uma corrente alternada de frequência angular \omega circula por um condutor cilíndrico maciço, a lei da indução eletromagnética de Faraday e a lei circuital de Ampère-Maxwell impõem uma distribuição não uniforme da densidade de corrente \vec{J}. Este fenômeno, conhecido como efeito pelicular ou skin effect, é regido pela equação de difusão magnética derivada diretamente do rotacional do campo magnético \vec{H} e do campo elétrico \vec{E}:

×H=J+ϵEt\nabla \times \vec{H} = \vec{J} + \epsilon \frac{\partial \vec{E}}{\partial t}

Desprezando a corrente de deslocamento em relação à corrente de condução em metais de alta condutividade, como o cobre recozido ou o alumínio de grau elétrico, e aplicando a relação constitutiva de Ohm microscópica \vec{J} = \sigma \vec{E}, obtém-se a equação de Helmholtz para o campo elétrico harmônico no interior do condutor:

2E=jωμσE\nabla^2 \vec{E} = j \omega \mu \sigma \vec{E}

Onde \mu = \mu_0 \mu_r representa a permeabilidade magnética absoluta do condutor e \sigma a condutividade elétrica. O parâmetro fundamental que governa a profundidade de penetração das ondas eletromagnéticas no metal é a profundidade de skin (\delta), definida matematicamente como:

δ=2ωμσ=1πfμσ\delta = \sqrt{\frac{2}{\omega \mu \sigma}} = \sqrt{\frac{1}{\pi f \mu \sigma}}

Para um cabo de média tensão típico operando em frequência industrial de f=50Hzf = 50 Hz ou 60Hz60 Hz, com condutor de cobre onde \sigma \approx 5.8 \times 10^7 S/m a 20°C20 °C e \mu_r \approx 1, a profundidade de skin resulta na ordem de 9.35mm9.35 mm a 50Hz50 Hz. À medida que a seção transversal do condutor é incrementada para atender a elevadas capacidades de condução de corrente em redes de distribuição subterrânea, o raio do condutor r0r_0 supera consideravelmente \delta. Isso faz com que a densidade de corrente decaia exponencialmente da periferia para o núcleo do condutor, segundo a função de Bessel modificada de ordem zero:

J(r)=J0ber(2rδ)+jbei(2rδ)ber(2r0δ)+jbei(2r0δ)J(r) = J_0 \frac{ber(\sqrt{2}\frac{r}{\delta}) + j bei(\sqrt{2}\frac{r}{\delta})}{ber(\sqrt{2}\frac{r_0}{\delta}) + j bei(\sqrt{2}\frac{r_0}{\delta})}

Esse confinamento espacial da energia de corrente alternada gera um acréscimo direto na resistência efetiva do condutor em relação à resistência ôhmica em corrente contínua RdcRdc. O fator de incremento de resistência por efeito pelicular é habitualmente denotado por ysy_s, de modo que a resistência em corrente alternada é expressa por Rac = Rdc(1 + y_s). Da mesma forma, o fluxo de energia eletromagnética, quantificado pelo vetor de Poynting \vec{S} = \vec{E} \times \vec{H}, experimenta uma distorção vetorial significativa, concentrando o fluxo de potência ativa e reativa através do meio isolante circundante (polietileno reticulado XLPE ou etileno-propileno EPR) de maneira heterogênea.

Interação Eletromagnética Multicondutora e Efeito Proximidade

Enquanto o efeito pelicular é um fenômeno intrínseco de um condutor isolado, em sistemas trifásicos de média tensão — seja em configurações de ternas compactas de três cabos unipolares ou em um cabo tripolar blindado — opera simultaneamente o efeito proximidade (proximity effect). Esse fenômeno surge porque o campo magnético alternado gerado pelas correntes que circulam nos condutores adjacentes induz correntes parasitas (correntes de Foucault) e redistribui as linhas de corrente dentro da seção transversal do condutor analisado.

A formulação analítica do efeito proximidade exige o emprego de coordenadas bipolares e a expansão em séries de funções harmônicas cilíndricas para resolver as equações de campo na presença de múltiplas fontes. O fator de incremento de resistência devido ao efeito proximidade é denotado como ypy_p. De acordo com as formulações normatizadas internacionais da norma IEC 60287, o fator total de perdas suplementares por corrente alternada no condutor é calculado somando-se ambos os efeitos sob condições de temperatura de operação especificada:

Rac=R0[1+ys+yp]Rac = R0 \left[ 1 + y_s + y_p \right]

Onde R0R_0 é a resistência em corrente contínua do condutor na temperatura máxima de operação. O termo ypy_p depende fortemente da geometria de disposição dos cabos (disposição em trevo ou em formação plana horizontal), da distância entre centros dos condutores ss, do raio do condutor rr e da profundidade de penetração \delta. Para condutores tubulares ou maciços circulares, o fator de proximidade é aproximado por expressões da forma:

yp=xp2(rs)2Fp(x)y_p = x_p^2 \left( \frac{r}{s} \right)^2 F_p(x)

Onde xpx_p é o parâmetro de espessura de pelicularidade corrigido e F_p(x) é uma função analítica tabulada em normativas técnicas. Em arranjos em formação plana com espaçamento reduzido, os campos magnéticos se somam vetorialmente de maneira assimétrica, fazendo com que a densidade de corrente se desloque para as zonas externas ou internas do condutor conforme a fase relativa das correntes circundantes. Isso gera gradientes térmicos severos ao longo da circunferência do condutor, acelerando a degradação termomecânica das camadas semicondutoras extrudadas e do isolamento principal de XLPE.

Matriz Paramétrica e Limites Normativos IEEE / IEC

O dimensionamento e a avaliação operacional dos condutores de média tensão exigem rigorosa conformidade com as normas internacionais IEC 60287 e IEEE Std 835 / IEEE Std 399. A matriz comparativa a seguir detalha os parâmetros críticos, limites operacionais e as consequências dielétricas associadas ao comportamento das perdas por efeito pelicular e proximidade.

Parâmetro Eletrotécnico Limite Padrão IEC 60287 Limite Padrão IEEE (Std 835 / 399) Condição Crítica de Falha Consequências Operacionais e Dielétricas
Fator de Incremento por Skin (ysy_s) Máx. 0,15 para seções \le 630 mm ^2 a 50Hz50 Hz Máx. 0,18 sob fator de potência nominal Sobrecarga harmônica acumulativa (THD_i > 8\%) Perdas Joule não contabilizadas, pontos quentes (hot-spots) locais e contração térmica do XLPE.
Fator de Incremento por Proximidade (ypy_p) Dependente do espaçamento s/d \ge 2,5 Restrição rigorosa em dutos fechados e circuitos em formação plana Circuitos em formação plana sem espaçamento adequado Deformação assimétrica do campo elétrico, estresse tangencial na blindagem semicondutora.
Temperatura Máxima do Condutor 90°C90 °C (XLPE) / 105°C105 °C (EPR) em regime permanente 90°C90 °C contínua, 130°C130 °C emergência Excursão térmica por curto-circuito ou distorção harmônica elevada Termofusão das blindagens semicondutoras, perda irreversível da rigidez dielétrica do polímero.
Seção Transversal Crítica >500mm2> 500 mm ^2 em cobre / >800mm2> 800 mm ^2 em alumínio >1000kcmil> 1000 kcmil Uso de cabo maciço ou compactado padrão sem segmentação Impedância de sequência positiva incorreta, instabilidade na atuação de proteções de distância.

Análise Forense de Engenharia de Falhas em Sistemas de Média Tensão

As falhas catastróficas em redes subterrâneas de média tensão (operando tipicamente em níveis de 13.8kV13.8 kV, 23kV23 kV, 34.5kV34.5 kV ou até 69kV69 kV) frequentemente têm sua causa raiz na análise inadequada dos fenômenos de alta frequência associados ao efeito pelicular e proximidade. Quando uma subestação industrial alimenta cargas altamente não lineares (conversores de frequência de média tensão, fornos a arco, inversores solares fotovoltaicos ou eólicos de grande porte), o espectro harmônico da corrente introduz componentes de frequências múltiplas da fundamental (h \cdot f).

Considerando a equação da profundidade de skin, torna-se evidente que, para a harmônica de ordem hh, a profundidade de penetração se reduz drasticamente:

δh=δ1h\delta_h = \frac{\delta_1}{\sqrt{h}}

Para uma harmônica de ordem h=25h = 25 (1250Hz1250 Hz em um sistema de 50Hz50 Hz ou 1500Hz1500 Hz em 60Hz60 Hz), a profundidade de penetração é reduzida por um fator de 5. Isso confina a corrente harmônica a uma camada extremamente delgada na periferia do condutor. A resistência efetiva do condutor para essa componente harmônica se eleva exponencialmente, provocando aquecimento localizado severo (Efeito Joule concentrado). Esse fenômeno desencadeia os seguintes mecanismos forenses de falha:

  • Degradação Térmica da Isolação XLPE: O calor gerado no interior do condutor não é dissipado com rapidez suficiente devido à baixa condutividade térmica do próprio isolamento polimérico (\lambda \approx 0.4 W/m \cdot K ). Isso gera um gradiente térmico radial abrupto, elevando a temperatura da interface entre o condutor e a camada semicondutora interna acima dos limites de projeto de reticulação do polímero.
  • Formação de Microcavidades e Arborescências Elétricas (Electrical Treeing): Os ciclos térmicos contínuos induzidos pelas perdas suplementares geram esforços termomecânicos de expansão e contração diferencial entre o metal do condutor e o isolamento extrudado. Isso causa a delaminação micromecânica da blindagem semicondutora, criando vazios de ar onde se iniciam descargas parciais (DP). As descargas parciais atuam como mecanismo precursor para o crescimento de arborescências elétricas (electrical treeing), culminando na perfuração dielétrica e no curto-circuito franco para a terra.
  • Falhas em Emendas e Terminações: Os conectores mecânicos ou de compressão utilizados em emendas e terminações de média tensão são projetados sob a premissa de uma distribuição de corrente relativamente homogênea. O efeito pelicular severo e a assimetria do efeito proximidade nas extremidades das conexões concentram as linhas de corrente nas paredes de contato dos conectores, provocando a fusão do composto inibidor de corrosão, a oxidação das superfícies de contato, o aumento progressivo da resistência de contato (efeito em cascata) e a destruição térmica do acessório.

Estratégias de Projeto Avançado, Mitigação e Condutores Segmentados

Para mitigar os efeitos adversos do efeito pelicular e de proximidade em condutores de alta capacidade de condução de corrente em média tensão, a engenharia moderna implementa soluções construtivas avançadas no projeto do condutor e nas configurações geométricas da instalação.

Condutores Segmentados tipo Milliken

Quando as seções transversais superam 800mm2800 mm ^2 em cobre ou 1200mm21200 mm ^2 em alumínio, o uso de condutores cilíndricos compactados convencionais torna-se inviável devido ao elevado fator ys+ypy_s + y_p. A solução padrão da indústria é o emprego de condutores segmentados tipo Milliken. Esse projeto consiste em dividir a área transversal total do condutor em quatro ou mais setores isolados eletricamente entre si por meio de uma fina camada de fita semicondutora ou isolante. Além disso, cada setor é encordoado com um passo helicoidal que transpõe periodicamente a posição física dos setores ao longo do comprimento do cabo.

Eletrodinamicamente, essa segmentação força cada setor a conduzir uma fração equitativa da corrente total, neutralizando os caminhos preferenciais das correntes de Foucault induzidas pelo efeito pelicular e pelo campo magnético de proximidade. Matematicamente, a resistência equivalente em corrente alternada do condutor Milliken é reduzida expressivamente, aproximando-se de forma assintótica da resistência em corrente contínua RdcRdc:

Rac,MillikenRdc(1+ys,sectorn2)R_{ac,Milliken} \approx Rdc \left( 1 + \frac{y_{s,sector}}{n^2} \right)

Onde nn é o número de setores independentes (tipicamente 4, 5 ou 6). Essa tecnologia permite reduzir as perdas por efeito Joule em mais de 40% em cabos de grande seção nominal.

Otimização Geométrica de Ternas Subterrâneas

Em instalações de cabos unipolares de média tensão em dutos, canaletas ou diretamente enterrados, a disposição geométrica das fases define a magnitude do efeito proximidade. A disposição em trevo (triangular compacta) minimiza o enlace de fluxo magnético assimétrico em comparação com disposições em formação plana horizontal ou vertical. Contudo, a disposição em trevo intensifica o acoplamento térmico mútuo entre os cabos, reduzindo a capacidade de dissipação térmica para o solo circundante. O cálculo rigoroso da ampacidade deve resolver simultaneamente o sistema de equações térmicas acopladas de Neher-McGrath:

Δθ=I2RacT1+I2Rac(1+ys)T2+[I2Rac(1+ys)+Wd]T3\Delta \theta = I^2 Rac T_1 + I^2 Rac (1 + y_s) T_2 + \left[ I^2 Rac (1 + y_s) + W_d \right] T_3

Onde T1,T2,T3T_1, T_2, T_3 representam as resistências térmicas do isolamento, da cobertura/duto e do solo circundante, respectivamente, e WdW_d expressa as perdas dielétricas do isolamento.

Aplicação Prática e Análise Computacional via Vexten Suite

Para ilustrar a aplicação dos princípios teóricos apresentados, demonstra-se um estudo de engenharia desenvolvido por meio dos módulos de cálculo da plataforma de engenharia Vexten Suite, aplicando as normas internacionais IEC 60909 para dimensionamento a esforços de curto-circuito e IEC 60287 / NEC 310 para a capacidade de condução de cabos sob fatores de desclassificação térmica por harmônicos.

Definição do Cenário Industrial

Analisa-se uma rede subterrânea de distribuição industrial operando em tensão nominal de 34.5kV34.5 kV trifásica, na frequência fundamental de 60Hz60 Hz. A subestação principal alimenta uma carga retificadora de grande porte que injeta expressivo conteúdo harmônico de corrente no sistema. Os parâmetros de projeto do circuito crítico são apresentados a seguir:

  • Tensão Nominal de Operação: Vn=34.5kVV_n = 34.5 kV (Fase-Fase)
  • Corrente da Frequência Fundamental: I1=650AI_1 = 650 A
  • Espectro Harmônico de Corrente: h_5 = 18\%,\, h_7 = 12\%,\, h11 = 7\%,\, h13 = 5\%
  • Distorção Harmônica Total de Corrente (THDiTHD_i): 23.47\%
  • Cabo Selecionado: Cabo unipolar com condutor de Cobre Segmentado tipo Milliken, isolamento em XLPE, blindagem metálica em fios de cobre e cobertura externa em PEAD (HDPE).
  • Seção Transversal Inicialmente Avaliada: 630mm2630 mm ^2
  • Temperatura Ambiente do Solo: 30°C30 °C
  • Resistividade Térmica do Solo: \rho_t = 1.0 K \cdot m/W

Execução do Algoritmo de Cálculo no Vexten Suite

O núcleo de cálculo eletromagnético e térmico do Vexten Suite processa individualmente cada ordem harmônica, calculando a resistência em corrente alternada específica para cada componente de frequência por meio da formulação de Bessel:

Irms=h=1nIh2=6502+(0.18×650)2+(0.12×650)2+(0.07×650)2+(0.05×650)2=670.85AIrms = \sqrt{\sum_{h=1}^{n} I_h^2} = \sqrt{650^2 + (0.18 \times 650)^2 + (0.12 \times 650)^2 + (0.07 \times 650)^2 + (0.05 \times 650)^2} = 670.85 A

Em seguida, o software avalia o fator total de perdas suplementares ponderado pelo espectro harmônico. Para o condutor convencional compactado de 630mm2630 mm ^2, o fator de skin a 60Hz60 Hz é ys,1=0.12y_{s,1} = 0.12; todavia, para a componente de 13a13^{ a } ordem (h=13h = 13, correspondente a 780Hz780 Hz), a profundidade de penetração cai para \delta13 = 2.59 mm , elevando o fator de skin harmônico para ys,13=0.89y_{s,13} = 0.89.

O algoritmo do Vexten Suite aplica o Fator de Redução Harmônica (Harmonic Derating Factor - HDF) formulado com base nas diretrizes das normas IEEE Std 539 e IEC 60287-2-2:

HDF=I12+h=2nIh2I12+h=2nIh2(Rac(h)Rac(1))HDF = \sqrt{\frac{I_1^2 + \sum_{h=2}^{n} I_h^2}{I_1^2 + \sum_{h=2}^{n} I_h^2 \left(\frac{Rac(h)}{Rac(1)}\right)}}

Ao se parametrizar um condutor convencional de cobre compactado de 630mm2630 mm ^2, o HDF calculado pela plataforma resulta no valor crítico de 0.780.78. Isso indica que a ampacidade do condutor sofre uma redução de 22\% decorrente unicamente da intensificação do efeito pelicular e de proximidade sob regime harmônico. A temperatura de regime do condutor atinge 96.4°C96.4 °C, excedendo o limite normativo contínuo de 90°C90 °C estabelecido para o isolamento em XLPE.

Ação Corretiva Implementada no Vexten Suite

Diante da não conformidade térmica, o motor de otimização paramétrica do Vexten Suite reespecifica a solução para um condutor de Cobre Segmentado tipo Milliken de 800mm2800 mm ^2. Sob essa nova topologia construtiva:

  • O fator de incremento por efeito pelicular na frequência fundamental reduz-se para ys,1=0.03y_{s,1} = 0.03.
  • A influência das harmônicas de ordem superior é atenuada de forma expressiva em função da segmentação dos setores, reduzindo o fator de perdas suplementares harmônicas em 65\%.
  • O Fator de Redução Harmônica (HDF) eleva-se para 0.950.95.
  • A temperatura de operação em regime permanente estacionário estabiliza-se em 82.1°C82.1 °C, atendendo plenamente aos critérios da norma IEC 60287 e assegurando uma vida útil projetada superior a 40 anos, isenta de riscos de formação de arborescências elétricas (electrical treeing) ou envelhecimento térmico prematuro da isolação polimérica.

Essa análise evidencia a necessidade indispensável da integração de modelos eletromagnéticos avançados baseados em equações de campo e difusão no dimensionamento de sistemas elétricos industriais e de distribuição, superando as limitações das tabelas simplificadas e assegurando máxima confiabilidade operacional em redes de média tensão sob condições reais de operação.