A existência de correntes circulantes em sistemas com inversores operando em paralelo representa um desafio crítico de engenharia. Estas correntes, que circulam entre as unidades sem contribuir para a carga, surgem de desbalanços sutis nos componentes, na temporização das chaves e nos caminhos parasitas. Sua manifestação não apenas reduz a eficiência do sistema, mas também pode provocar sobreaquecimento, estresse térmico nos semicondutores e, em casos extremos, falhas catastróficas. O diagnóstico e a mitigação eficazes dependem de uma compreensão profunda de sua natureza e origem, frequentemente abordada por meio de simulação detalhada.

Para investigar este fenômeno, é comum construir um modelo de simulação dinâmica representando dois inversores conectados em paralelo a uma mesma barra CA e a uma fonte de tensão contínua compartilhada. A precisão do modelo reside na parametrização adequada de elementos que governam o comportamento das correntes parasitas. Parâmetros fundamentais incluem a tensão do barramento CC, a frequência de modulação, as resistências parasitas nos caminhos de corente e, de maneira crucial, os tempos mortos e as indutâncias nos laços de circulação.

% Definição dos parâmetros do sistema
V_barramento = 800; % Tensão no barramento CC (V)
freq_modulacao = 20e3; % Frequência de modulação das portadoras (Hz)
R_parasita = 0.015; % Resistência parasita dos caminhos de corrente (Ohm)
L_circulacao = 45e-6; % Indutância no caminho das correntes circulantes (H)
tempo_morto = 1.8e-6; % Tempo morto configurado nos drivers (s)

Um aspecto crucial da modelagem é a implementação precisa do tempo morto na geração dos sinais de acionamento. O uso de blocos de atraso de transporte puro é essencial, pois blocos de atraso genéricos ou inserções inadequadas podem introduzir erros de fase sistemáticos, distorcendo significativamente os resultados da simulação e invalidando a correlação com medições experimentais.

Após a execução da simulação, a análise foca na extração e decomposição das correntes. A corrente circulante de sequência zero é um indicador chave e pode ser isolada das correntes de linha medidas. O cálculo envolve a média fasorial das correntes trifásicas, seguido da remoção de qualquer componente de valor médio CC que possa ter sido induzida por erros numéricos ou configurações do modelo.

function [I_seq_zero] = extrair_corrente_circulante(Ifase_A, Ifase_B, Ifase_C)
    % Calcula a componente de média das correntes de fase, representando a
    % corrente de sequência zero circulante.
    I_soma_total = Ifase_A + Ifase_B + Ifase_C;
    I_seq_zero = I_soma_total / 3;
    
    % Remove desvios de DC residuais usando filtragem de média móvel
    tamanho_filtro = 100;
    filtro_media = ones(1, tamanho_filtro) / tamanho_filtro;
    offset_cc = filter(filtro_media, 1, I_seq_zero);
    I_seq_zero = I_seq_zero - offset_cc;
end

A análise no domínio da frequência revela informações que a enspeção temporal pode ocultar. A assinatura espectral da corrente circulante, particularmente nas proximidades da frequência de comutação e seus múltiplos, ajuda a identificar a origem do desbalanço. Uma rotina de análise espectral eficiente utiliza métodos de estimativa de densidade espectral de potência com janelamento adequado para capturar as bandas laterais características.

% Análise espectral da corrente circulante
[psd_estimada, vet_freq] = pwelch(I_seq_zero, hanning(4096), 2048, 4096, freq_amostragem);

% Definir bandas de interesse ao redor da frequência de comutação
limite_inferior = 0.85 * freq_modulacao;
limite_superior = 1.15 * freq_modulacao;
indices_harmonicas = find(vet_freq > limite_inferior & vet_freq < limite_superior);

% Verificar a presença de harmônicas relevantes
if max(psd_estimada(indices_harmonicas)) > 1e-5
    disp('Alerta: Harmônicas significativas detectadas na banda de comutação.');
end

Entre as estratégias de mitigação, a otimização da técnica de modulação é uma das mais efetivas. Em modulação por vetor espacial (SVPWM), a inserção deliberada de uma componente de tensão de sequência zero no sinal de referência pode contrabalançar parcialmente as assimetrias que causam a circulação. Esta componente de compensação é calculada a partir das tensões de referência instantâneas.

% Cálculo da tensão de offset para compensação de sequência zero em SVPWM
V_ref_max = max([V_alfa_ref, V_beta_ref]); % Supondo referências em coordenadas alfa-beta
V_ref_min = min([V_alfa_ref, V_beta_ref]);
V_compensacao_zero = -0.5 * (V_ref_max + V_ref_min);
% Este offset é adicionado às referências antes do cálculo dos tempos de chaveamento.

A eficácia de qualquer medida corretiva deve ser validada experimentalmente. Resultados comparativos típicos mostram uma redução progressiva nos picos de corrente circulante: de uma configuração inicial não otimizada, passando por ajustes no tempo morto, até a implementação de algoritmos de compartilhamento de corrente ativo, onde reduções superiores a 70% são possíveis.

Uma lição fundamental na análise de falhas em campo é que a assinatura da corrente circulante pode estar mais evidente no domínio da frequência. A presença de componentes em frequências específicas, como o dobro da frequência da rede, pode apontar diretamente para problemas de desbalanço na impedância de paralelismo, como por exemplo, uma conexão de resistência elevada devido a um contato frouxo no dissipador de calor, um defeito que pode ser quase imperceptível na forma de onda temporal.