Method Article

Estudo Comparativo de Simulação de Aumento de Temperatura na Unidade Principal do Anel

DOI:

10.3791/66643

July 5th, 2024

In This Article

Summary

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Este artigo aborda o problema de aumento de temperatura da unidade principal do anel, estabelecendo um modelo simplificado e conduzindo uma análise comparativa em dois módulos de resolução de campo de temperatura.

Abstract

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

A Unidade Principal em Anel (RMU) é um dispositivo crítico em sistemas de distribuição de energia usados para conectar e distribuir eletricidade. No entanto, devido à sua estrutura interna compacta e alta carga de corrente, os problemas de dissipação de calor são particularmente proeminentes. Para resolver esse problema, este estudo propõe de forma inovadora um modelo RMU simplificado, empregando métodos de simulação de elementos finitos para resolver com precisão as perdas ôhmicas de condutores em condições reais de operação e obter dados de perda ôhmica para vários componentes. Esta é a primeira investigação aprofundada do problema de aumento de temperatura da RMU usando uma abordagem tão abrangente. Posteriormente, o campo de temperatura foi resolvido usando dois módulos diferentes de análise de campo de temperatura, com uma comparação detalhada e análise dos resultados da simulação para identificar semelhanças, diferenças e tendências na distribuição de temperatura. Os resultados indicam que o modelo de solução de campo de temperatura, que considera a transferência de calor por convecção, é mais preciso e se alinha com as condições reais de operação. Esta pesquisa fornece uma abordagem inovadora e soluções práticas para o projeto e otimização de RMUs. Pesquisas futuras podem explorar ainda mais os métodos de análise de acoplamento multifísico para abordar questões de projeto estrutural e validação obrigatória para RMUs de alta e ultra-alta tensão e outros equipamentos elétricos, fornecendo assim informações importantes para o projeto de engenharia.

Introduction

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

A unidade principal do anel é um grupo de aparelhagem de alta tensão montada em um gabinete de metal de aço ou feita de unidade de fonte de alimentação de rede de anel espaçado montada de equipamentos elétricos. A estrutura geral do interruptor de carga e do circuito condutor consiste no circuito condutor, que inclui vários componentes que compõem o núcleo principal da unidade de anel. No entanto, devido à sua estrutura interna compacta, a unidade principal em anel enfrenta desafios na dissipação de calor. Isso pode levar à deformação térmica e envelhecimento ao operar por longos períodos em ambientes de alta temperatura. Esses problemas não afetam apenas a vida útil da unidade, mas também afetam suas propriedades isolantes, apresentando riscos à segurança. Em particular, danos ao equipamento e acidentes elétricos tornam-se mais prováveis, representando riscos significativos à segurança.

Em diferentes áreas de pesquisa, os estudiosos realizaram uma série de estudos sobre o aumento da temperatura dos painéis de distribuição de linhas aéreas e analisaram vários fatores que afetam a distribuição de temperatura1. Em Polykrati et al.2, é apresentado um modelo matemático para a estimativa do aumento de temperatura de componentes instalados na rede de distribuição durante uma falha de curto-circuito. O modelo foi aplicado às chaves seccionadoras comuns da rede, e as características dos resultados foram plotadas de acordo com as diferentes formas da parte assimétrica da forma de onda da corrente de curto-circuito e o valor inicial do componente de corrente CC de curto-circuito. Guan et al., por outro lado, levaram em consideração a resistência de contato e a repulsão eletromagnética construindo uma ponte de contato equivalente para simular a interface de contato e analisaram ainda mais o campo de acoplamento eletromagnético-térmico e o experimento de aumento de temperatura3. Além disso, os pesquisadores investigaram o campo de temperatura e a distribuição de tensões térmicas dos contatos dinâmicos e estáticos dentro da unidade principal do anel por simulação de elementos finitos, que forneceu uma base para o estudo da vida útil do disjuntor4. Finalmente, Mueller et al. concentraram-se nas características geométricas dos dissipadores de calor e avaliaram os efeitos da seleção de materiais, área de superfície total, uniformidade de temperatura e temperatura máxima da superfície no desempenho térmico5. Esses estudos fornecem informações e métodos valiosos para melhorar o desempenho e a confiabilidade do painel de distribuição, reduzir o aumento da temperatura e prolongar a vida útil do equipamento. Wang et al. propuseram um MiNET Deep Learning Model (MDLM) no ambiente UPIOT com o objetivo de detectar o diagnóstico de falhas em gabinetes de anéis elétricos, que foi validado para ter uma precisão de identificação de 99,1%, que é significativamente maior do que a de outros métodos6. Lei et al. estudaram o desempenho térmico de um barramento GIS em estado estacionário usando o método de análise de acoplamento magneto-fluido-térmico, otimizando assim o diâmetro do condutor e do tanque com base nos resultados da simulação de aumento de temperatura7. Ouerdani et al. usaram o modelo de simulação de aumento de temperatura da RMU para determinar o aumento de temperatura em locais críticos dentro dela, fixando assim a duração da sobrecarga máxima para os componentes dentro da RMUde acordo 8. Zheng et al. descreveram um barramento retangular convencional em um modelo de painel de alta corrente, construindo um modelo bidimensional e aplicando o método dos elementos finitos (FEM) para cálculos de campo eletromagnético. Isso permitiu que eles obtivessem a distribuição da densidade de corrente do condutor do barramento e da perda de potência. Um barramento irregular foi projetado após considerar os efeitos do efeito de proximidade e efeito de pele. Este design de barramento irregular melhorou o desempenho do barramento retangular convencional9.

Quanto ao aspecto do uso da simulação do icepak, Wang et al. realizaram uma simulação de aumento de temperatura por meio de teorias de campo de vórtice, campo de fluxo de ar e campo de temperatura e descobriram que o aumento de temperatura da unidade principal do anel era mais grave sob convecção natural. Eles reduziram com sucesso o nível de aumento de temperatura adicionando resfriamento de ar forçado e fazendo melhorias na estrutura de contato interno10. Zhu et al.11 utilizaram o icepak para simular um modelo térmico a fim de comparar o efeito da presença de vias térmicas na PCB e a presença de dissipadores de calor na temperatura dos dispositivos de potência. Por fim, a análise teórica é comparada com os resultados da simulação para verificar a exatidão da análise teórica. Mao et al.12 estudaram a temperatura e a distribuição do fluxo de ar interno em condições de operação no verão por simulação térmica baseada no software CAE na simulação icepak. O problema de como melhorar a eficiência de resfriamento e controlar o aumento de temperatura de vários contatos prateados é dado, e os contornos de temperatura e fluxo de ar interno capturados na simulação estabelecerão as bases para o projeto do esquema de resfriamento para os seis contatos prateados montados na unidade de vedação. Por outro lado, no uso de um módulo térmico de estado estacionário, os métodos de modelagem de Zhang13 são discutidos para resolver a rede térmica de uma bucha de alta pressão usando um procedimento transitório alternativo. Os resultados do teste e da simulação estão de acordo com o estado estacionário térmico e os estados transitórios da bucha. Os resultados transitórios são então usados para avaliar a capacidade de sobrecarga da bucha. Vaimann et al.14 desenvolveram e analisaram um modelo térmico analítico de um motor de relutância síncrona para prever a temperatura de seus diferentes componentes e a rede térmica do parâmetro total definido.

Com o avanço contínuo da pesquisa em equipamentos elétricos, como unidades principais de anel, os testes convencionais de aumento de temperatura e os métodos de produção são relativamente ineficientes. Portanto, ao utilizar a tecnologia de elementos finitos combinada com testes off-line, não apenas os problemas de custo do projeto são resolvidos, mas ajustes e otimizações podem ser feitos prontamente para problemas do mundo real com base em simulações. Com base no progresso da pesquisa mencionado acima, o uso do ANSYS Icepak e do acoplamento térmico de estado estacionário para análise comparativa raramente é mencionado. Portanto, o protocolo descreve a pesquisa de mecanismos de elementos finitos, usa combinações numéricas e morfológicas para estabelecer um modelo de simulação de aumento de temperatura de elementos finitos para o recinto e discute o modelo de simulação de aumento de temperatura de elementos finitos com base nos resultados dos dois módulos analíticos, comparando os resultados dos dois módulos de simulação. Através da comparação entre os dois módulos de simulação, obteremos as características da tendência de aumento de temperatura da unidade principal do anel e encontraremos o método mais aplicável para fornecer a base necessária e ideias de pesquisa para uma estratégia de mitigação do aumento de temperatura da unidade principal do anel.

Protocol

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

1. Modelo

NOTA: Devido à estrutura complexa da unidade principal do anel (Figura 1A), um software de design on-line foi escolhido para simplificar a operação da unidade principal do anel.

  1. Simplificação da modelagem
    1. Simplifique parcialmente o modelo, preservando a seção da caixa de ar da RMU enquanto remove ou simplifica outros componentes, como eixos isolantes, parafusos de fixação, porcas, componentes de vedação e suportes de suporte de pressão. A versão simplificada é mostrada na (Figura 1B).
      1. No processo de simplificação da unidade principal de anel tipo 630A, remova o eixo isolado que conecta a sala do disjuntor com a caixa de instrumentos e muitos parafusos e porcas fixos. Retire as peças de vedação e o suporte de retenção de pressão e conecte os contatos estáticos da viga estática isolada com o barramento do ramal inferior sob a premissa de garantir que toda a configuração tenha a mesma corrente condutora e apenas o disjuntor a vácuo, a placa de fixação do disjuntor e os contatos estáticos e o disjuntor a vácuo sejam retidos.
      2. Guarde apenas o disjuntor a vácuo, a placa de fixação do disjuntor, o contato estático e a placa de bloqueio do disjuntor a vácuo. No geral, remova os parafusos e juntas do modelo, preencha os orifícios após remover os parafusos com sólidos, reduza o número de peças de malha e otimize as formas irregulares das peças. Remova os instrumentos para operação do painel, placas de montagem, suportes e outras peças operacionais, como caixas de instrumentos, que não afetam o processo de simulação de aumento de temperatura.
      3. A remoção dos mancais isolados de alguns componentes pode ser desconsiderada na simulação, pois tem pouco efeito nos resultados da simulação. Além disso, as chaves de aterramento que não afetam o uso do equipamento durante a operação normal as removem e mantêm a sala do disjuntor para simulação.
    2. Para excluir qualquer seção, basta selecioná-la e clicar na opção Excluir .

2. Solução de campo parasita

  1. Configurações de pré-processamento
    NOTA: A emulação de campo de corrente parasita é a base para realizar a solução de campo de temperatura, que requer a análise subsequente da fonte de calor resolvida como uma carga no campo de temperatura.
    1. Consulte a documentação do equipamento da unidade principal em anel e os manuais relevantes para obter informações sobre as propriedades físicas e parâmetros de cada componente da unidade principal em anel. Defina os atributos físicos e parâmetros dos componentes da unidade principal do anel em Maxwell com base nas informações obtidas, conforme detalhado na Tabela 1.
    2. Defina a corrente de carga de campo da corrente parasita em 630 A com uma frequência de 50 Hz. No software Maxwell, selecione Um lado dos braços de saída superior e inferior, entre no módulo de excitação e defina a magnitude da corrente para 630 A. Na seção de configurações da solução, escolha uma frequência de 50 Hz.
      NOTA: No circuito condutor de uma unidade principal em anel, o caminho formado por todos os componentes do braço de saída superior ao braço de saída inferior é conhecido como sequência de fases. Portanto, neste artigo, as fases A, B e C são organizadas da esquerda para a direita.
    3. Os parâmetros do material dos componentes da unidade principal do anel são mostrados na Tabela 2.
    4. Direcione a corrente através dos braços de linha de saída, conexões flexíveis, barramentos, disjuntores, barramentos de suporte de contato estático e barramentos de derivação para cada fase. O objetivo é realizar um caminho de corrente que permita que os componentes completem a carga.
    5. Utilize a malha adaptativa de Maxwell para completar o controle de grade para o modelo. Use o método de particionamento de malha adaptável Maxwell para componentes maiores e o refinamento de malha local para componentes internos menores.
      NOTA: Maxwell pode melhorar continuamente a precisão da grade durante o processo de resolução, eliminando a necessidade de clicar em Operações de malha para particionamento de malha adicional.
    6. Defina o tamanho da etapa da solução. Clique em Análise na árvore do modelo, abra as configurações da Etapa de resolução e defina o Número máximo de passagens como 10. Mantenha outras configurações em seus valores padrão sem fazer alterações.
  2. Princípio do cálculo de campo de correntes parasitas 15,16.
    1. Use a primeira equação de Maxwell, que descreve a ação da carga na geração de um campo elétrico17.
      figure-protocol-1(1)
      onde ρ representa a densidade de carga; ε0 representa a constante dielétrica do vácuo.
    2. Use a segunda equação de Maxwell, que descreve a relação entre um campo magnético variável e um campo elétrico e o efeito de um campo magnético no movimento de uma carga.
      figure-protocol-2(2)
      onde figure-protocol-3 representa a intensidade do campo magnético. Esta equação descreve que um campo magnético variável produz um campo elétrico de vórtice, ou seja, o spin do campo elétrico do vórtice é igual ao negativo da taxa de variação do campo magnético com o tempo.
    3. Use a terceira equação de Maxwell, que descreve o efeito da carga magnética na produção de um campo magnético.
      figure-protocol-4(3)
      Esta equação descreve o campo magnético produzido por uma carga magnética como sendo passivo, ou seja, não há monopolos no campo magnético.
    4. Use a quarta equação de Maxwell, que descreve a relação entre um campo elétrico variável e um campo magnético e o efeito de uma corrente elétrica em um campo magnético.
      figure-protocol-5(4)
      onde figure-protocol-6 representa a densidade de corrente e μ0 representa a permeabilidade ao vácuo. Esta equação descreve que um campo elétrico variável produz um campo magnético de vórtice, ou seja, o spin do campo magnético do vórtice é igual à soma da densidade de corrente e da taxa de variação do campo elétrico com o tempo.
    5. Com base nas equações acima, use o Maxwell 3D usando o módulo solucionador de correntes parasitas para resolver a perda ôhmica gerada pelo circuito condutor na RMU, que fornece uma fonte de calor para a análise de simulação térmica subsequente. Sua expressão matemática é dada como18
      figure-protocol-7(5)
      onde σ denota a condutividade do material do loop condutor; J é a densidade de corrente no loop.
  3. Resultados do cálculo
    1. Clique na opção Maxwell 3D na interface e abra a verificação de validação para revisar todas as configurações em busca de erros. Se não houver erros, clique em Analisar tudo para iniciar o processo de resolução.
    2. Utilize a calculadora de pós-processamento de Maxwell para calcular e representar graficamente as perdas ôhmicas no campo de correntes parasitas da unidade principal do anel, conforme mostrado na Tabela 3.

3. Solução de campo de temperatura

NOTA: Para fins comparativos, divida o campo de temperatura em Icepak e térmico de estado estacionário. Configure e resolva cada um separadamente para obter uma análise comparativa.

  1. Configuração do modelo Icepak
    1. Defina as propriedades do material da seguinte forma: designe todos os materiais sólidos do circuito como-Pure, com superfícies usando superfície polida com. Para os componentes do painel, selecione o material Aluminum6061-T6 , com um revestimento de superfície de superfície Paint-AL com uma emissividade de 0,35. Consulte a Tabela 4 para obter detalhes. Clique com o botão direito do mouse no componente selecionado, clique em Editar e vá para Propriedades para definir o material para a superfície e os materiais sólidos.
    2. Selecione o modelo e clique em Definir no menu Editar e, em seguida, escolha Nível de malha multinível para ajustar as configurações de malha. Defina o gabinete externo para um nível de malha de 2 e todos os limites para um nível de malha de 2. Para todos os outros componentes, defina o nível de malha como 3. Por fim, abra o controle de malha e clique em Gerar para criar a malha.
    3. Para garantir a precisão e eficiência da simulação, independentemente do tamanho da grade, é necessária a validação da independência da rede. Importe o modelo geométrico do invólucro do campo de temperatura, estabelecido usando o software de projeto para geração de malhas.
    4. Conforme ilustrado na Figura 2, as curvas de polarização dos quatro conjuntos de grade estão bem alinhadas. Com uma tensão de trabalho de 0,5 V, as densidades de corrente para os quatro conjuntos de grade são 2,357 A/cm2, 2,358 A/cm2, 2,356 A/cm2 e 2,454 A/cm2, respectivamente, com o erro entre as densidades máxima e atual sendo inferior a 1%. Para equilibrar eficiência e precisão, determine o tamanho da grade que vem a ser 987924.
  2. Configuração da solução
    1. Defina as direções do Gabinete do domínio da solução para Abertura.
    2. No software, selecione Etapa do problema. Em Parâmetros básicos, verifique o Modelo de radiação superfície-superfície, escolha Equação zero para Regime de fluxo turbulento, selecione a opção Gravidade para Convecção natural e defina a temperatura ambiente para 20 °C.
    3. Nas configurações de arquivo, escolha Perdas de calor volumétricas para mapeamento EM e selecione Todos os objetos mostrados para concluir as configurações de perda.
  3. Cálculo do campo de temperatura
    1. No icepak, aplique três equações principais de conservação de energia: equação de conservação de massa, equação de conservação de momento e equação de conservação de energia. Especificamente, use a equação de conservação do momento, que é a seguinte19:
      figure-protocol-8(6)
      Equação de Conservação de Energia:
      figure-protocol-9(7)
      Equações de conservação de massa:
      figure-protocol-10(8)
      Equação de transferência de energia para transferência de calor de uma fonte de calor sólida:
      figure-protocol-11(9)
      ρ representa a densidade do fluido; v representa o vetor velocidade do fluxo; T representa a temperatura; p é a pressão; τ é a força viscosa na superfície do micrometabólito; κ é o coeficiente de transferência de calor; Shé a fonte de calor do corpo; h é a entalpia específica do fluido e F é a força corporal do micrometabólito.
      NOTA: Os resultados dos cálculos do campo de temperatura são mostrados na Figura 3A e na Figura 4A.
  4. Configuração do modelo térmico de estado estacionário
    1. Mantenha as propriedades do material de acordo com a Tabela 3 nas configurações do material. Gere as perdas ôhmicas resultantes da análise de simulação de campo de correntes parasitas no módulo Térmico de Estado Estacionário clicando em Geração de Carga Térmica .
    2. Clique no valor da temperatura por convecção e defina-o para 20 °C, com um coeficiente de convecção de 5 (W/m²°C) aplicado às paredes internas do gabinete, componentes e gabinete externo. Aplique as configurações e gere. Defina a saída para resolver a temperatura clicando em Resolver > Resultados da Saída.
      NOTA: O principal campo de temperatura que rege as equações no princípio de cálculo do campo de temperatura térmica em estado estacionário 20,21,22 é geralmente derivado da lei da condução de calor (lei de condução de calor de Fourier). No caso unidimensional, a equação de transferência de calor do campo de temperatura pode ser expressa como20:
      figure-protocol-12(10)
      Nesta equação, T representa a temperatura dentro do objeto, t é o tempo, x são as coordenadas espaciais e α é a difusividade térmica. Esta equação descreve a variação da temperatura em relação ao tempo e ao espaço, onde o lado direito expressa a relação entre a taxa de condução de calor e o gradiente de temperatura. Em um cenário tridimensional mais geral, a equação de condução de calor para o campo de temperatura pode ser expressa da seguinte forma:
      figure-protocol-13(11)
      ρ representa a densidade do objeto, c é a capacidade térmica específica, K é a condutividade térmica e Q é o termo da fonte de calor dentro do volume. Esta equação descreve a variação do campo de temperatura, influenciada pela condução de calor, fontes de calor e capacitância térmica.
    3. Os resultados do cálculo do campo de temperatura são mostrados na Figura 3. Compare os valores de temperatura resumidos na Tabela 5 e na Tabela 6.

Results

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Com base nos dados da Tabela 3, as seguintes conclusões podem ser tiradas: As perdas gerais para as Fases A, B e C são relativamente semelhantes. Especificamente, as perdas totais para a Fase A são de 16.063 W/m³, a Fase B é de 16.12 W/m³ e a Fase C é de 19.57 W/m³. Os locais com maiores perdas podem estar nas conexões de vários componentes. Isso ocorre principalmente porque a resistência de contato e a resistência do condutor normalmente existem nesses pontos de conexão. Quando a corrente passa por essas conexões, um calor significativo é gerado, levando a um aumento da temperatura e maiores perdas nessas áreas.

Os braços de saída superior e inferior da unidade principal do anel suportam algumas perdas, especialmente ao transportar a carga principal. As perdas nos braços de saída das três fases são aproximadamente as mesmas. Isso ocorre porque essa parte da corrente é relativamente concentrada e, devido à forma regular, o valor da resistência é pequeno. Portanto, as perdas nessas partes são aproximadamente iguais.

As perdas no barramento do ramal são relativamente altas, principalmente devido às suas inúmeras curvas e à presença de seções angulares. A maior parte da corrente está concentrada na área da curva e perto dos cantos. Além da seção do disjuntor a vácuo, as perdas na parte do tubo de cobre do barramento também são relativamente altas, totalizando 22,32 W/m³, representando 40% das perdas ôhmicas totais das três fases. As perdas na parte de contato estático são relativamente pequenas em comparação com as perdas gerais.

Outros componentes, como a placa fixa do disjuntor, o defletor do disjuntor e o invólucro externo da unidade principal do anel, têm perdas menores. Por não participarem diretamente da carga, suas perdas são geradas principalmente pelos componentes internos por condução, resultando em menores valores de perda. Nos cálculos de pós-processamento, as perdas desses componentes não são particularmente detalhadas. Em resumo, as principais características da distribuição de perdas foram descritas, fornecendo informações valiosas para uma maior otimização do projeto.

Combinando os resultados da Figura 3 e da Figura 4, juntamente com a Tabela 5 e a Tabela 6, foi realizada uma comparação completa dos resultados da solução de campo de temperatura entre os módulos Icepak e Thermal Steady-state. No processo de análise, foram observadas diferenças e semelhanças entre os dois módulos.

Na faixa de temperatura de 20-31,24 °C, conforme mostrado na Figura 2A, as áreas de temperatura mais alta estão concentradas nas regiões onde o invólucro RMU entra em contato com os componentes internos. A principal razão é que essas áreas de contato normalmente servem como caminhos críticos para a condução de calor. Durante a operação, os componentes internos geram calor, que é conduzido para o invólucro através das superfícies de contato com boa condutividade térmica, causando assim um aumento da temperatura nessas regiões. Essas áreas incluem os pontos de contato do braço de linha de saída superior e as partes do braço de linha de saída inferior em contato com seu invólucro próximo, com uma faixa de temperatura de aproximadamente 24-27 °C. Comparado com a Figura 3B, a faixa de aquecimento das conchas em ambos os casos é aproximadamente a mesma, e as mudanças de temperatura se propagam dos pontos quentes para as áreas mais frias. Examinando a Figura 4A, B, as tendências gerais de temperatura em ambos os casos também são aproximadamente semelhantes. A principal fonte de aumento de temperatura é a condução de calor através do braço da linha de saída inferior, que passa pelos condutores internos e se dissipa do braço da linha de saída superior. A tendência de aquecimento está concentrada em áreas próximas ao braço inferior da linha de saída, incluindo os barramentos de ramificação de cada fase e os barramentos de suporte conectados. Além disso, a partir da observação da distribuição de temperatura nos circuitos internos na Figura 4, seja usando o Icepak ou o módulo térmico de estado estacionário para resolver, a temperatura geral da Fase B é consistentemente mais alta do que as outras duas fases. Isso sugere principalmente que a Fase B, durante o processo de carga, não apenas suporta o calor gerado pela corrente dessa fase, mas também, devido à estrutura interna excessivamente compacta da RMU, o calor produzido por cada fase não pode se dissipar prontamente. Devido à troca de calor, a temperatura da Fase B permanece mais alta do que as outras duas fases, independentemente do solucionador usado. As tendências gerais de temperatura também exibem um alto grau de semelhança, com o principal aumento de temperatura decorrente da condução de calor através do braço da linha de saída inferior e da dissipação de calor do braço da linha de saída superior. Essa tendência parece altamente consistente entre os dois módulos. Além disso, as posições dos pontos quentes nos mapas de distribuição de temperatura de ambos os módulos são altamente consistentes. Particularmente, na seção de barramento de ramificação da Fase C, ambos os solucionadores mostram a mesma temperatura mais alta e as diferenças são quase insignificantes. Isso implica que, independentemente do solucionador usado, a identificação precisa das posições de temperatura mais quentes é crucial para o gerenciamento térmico da RMU.

Em segundo lugar, pode-se ver na tabela que a mesma corrente produz a mesma perda, mas o valor da temperatura de cada componente de fase na Tabela 6 é geralmente menor do que na Tabela 5; por exemplo, o valor da temperatura do barramento de ramificação na área do ponto quente de toda a RMU é a temperatura mais alta na Tabela 6 a 30,91 °C e, por outro lado, a temperatura mais alta do barramento de ramificação na Tabela 5 é de 31,24 °C. O motivo é a lógica da solução do Icepak: a RMU como fonte de calor continuará a trocar calor e dissipar calor, e quando um fluido é definido como meio para transferência de calor, a temperatura no domínio da solução será gradualmente emitida através do meio, o que resulta em uma diminuição na temperatura. O módulo térmico de estado estacionário coloca menos ênfase na convecção de calor com o ar circundante e, em vez disso, se concentra em uma solução baseada em condução de calor. Em comparação com os modelos que levam em conta a convecção, essa abordagem resulta em uma solução de campo de temperatura menos abrangente. Consequentemente, em áreas onde a temperatura é mais alta, os pontos quentes tornam-se mais pronunciados. No processo do teste de aumento de temperatura real, além de considerar as características do próprio modelo, a transferência de calor ao mesmo tempo também considera a influência do meio de dissipação de calor, como ar, etc., na temperatura geral. Combinado com as necessidades práticas e considerações de operação experimental, a RMU como fonte de calor continua a realizar convecção e dissipação de calor; A simulação de temperatura do ICEPAK está mais alinhada com as necessidades reais. No caso em que a unidade principal do anel é usada como fonte de calor com convecção e dissipação contínua de calor, a simulação de temperatura do icepak está mais alinhada com as necessidades práticas. Pelo contrário, o módulo térmico de estado estacionário se concentra principalmente na condução de calor, que pode não ser capaz de atender às necessidades reais em alguns casos.

figure-results-1
Figura 1: Modelo da unidade principal em anel. (A) Modelo geral da unidade principal do anel (B) Modelo simplificado da unidade principal do anel. Clique aqui para ver uma versão maior desta figura.

figure-results-2
Figura 2: Curvas de polarização para quatro conjuntos de malhas. A malha adaptativa Maxwell é usada para peças maiores e o refinamento de malha local é usado para peças menores. No software de projeto, um modelo 3D simplificado do gabinete principal do anel é criado e importado para o Maxwell, e o modelo é gerado usando seu módulo Mesh. A figura representa a verificação das grades do modelo. A densidade de corrente dos quatro grupos de grades é de 2,357 A/cm2, 2,358 A/cm2, 2,356 A/cm2 e 2,354 A/cm2 sob a carga de tensão, e o erro relativo entre a densidade de corrente máxima e a densidade de corrente mínima é inferior a 1%, e para levar em consideração a eficiência e precisão dos cálculos, O número final das grades é determinado como 1169091. Clique aqui para ver uma versão maior desta figura.

figure-results-3
Figura 3: Modelo de análise de campo de temperatura do gabinete externo. (A) Distribuição da temperatura do invólucro com solução Icepak. (B) Distribuição da temperatura do invólucro com solução térmica em estado estacionário. Clique aqui para ver uma versão maior desta figura.

figure-results-4
Figura 4: Modelo de análise de campo de temperatura do circuito interno. (A) Distribuição de temperatura do circuito condutor com solução Icepak (B) Distribuição de temperatura do circuito condutor com solução térmica em estado estacionário. Clique aqui para ver uma versão maior desta figura.

MaterialArChapa de aço galvanizadoCobreAlumínio
Capacidade calorífica específica (J/(kg· K))1007500500897
Densidade (kg/m³)1.1614803089002689
Permeabilidade relativa10.30.9999991/
Permeabilidade magnética relativa1250011
Condutividade térmica (W / (m · K))0.02616386237
Emissividade/0.650.30.1
Condutividade elétrica (S/m)/0.85.80E+07/

Tabela 1: Parâmetros físicos de alguns materiais na unidade principal do anel.

ComponentesPermeabilidade relativaCondutividade/sMaterial
Circuitos condutores de corrente (barramento, barramento de derivação, etc.)03,00,00,000Liga de alumínio
Disjuntor a vácuo0.99992,00,00,000Ligas de cobre
Corpo20011,00,000Aço estrutural

Tabela 2: Lista de materiais para cada componente.

Perda ôhmica W/m³UmBC
Saia deixe o braço3.783.723.73
Barramento de ramificação2.12.092.1
Conexão flexível1.31.31.3
Disjuntor a vácuo1.0230.950.98
Contatos estáticos0.360.360.36
barramento alimentador1.331.351.32
Tubos de cobre para barramentos de ramificação6.196.359.78

Tabela 3: Valores de perda ôhmica para os componentes das fases A, B e C.

Número12
MateriaAlumínio6061-T6-Puro
Condutividade térmica (W / m · K)167387.6
Densidade (kg/m³)89332700
Capacidade calorífica específica/(kg· K)/896
Material da superfícieSuperfície Paint-AL-Lustrado-superfície
Emissividade0.350.052

Tabela 4: Configurações de parâmetros de material.

Ponto de monitoramento de temperatura/temperatura (°C)UmBC
Braço de saída superior25.4825.7925.63
Barramento principal25.9326.2826.13
Barramento de ramificação2929.1830.01
Tubo de cobre de barramento de ramificação31.0431.1831.24
Braço de saída inferior26.526.9826.92

Tabela 5: Valores de temperatura de cada fase da unidade principal do anel sob o módulo de resolução de campo de temperatura térmica em estado estacionário.

Ponto de monitoramento de temperatura/temperatura (°C)UmBC
Braço de saída superior23.7323.8223.81
Barramento principal25.1525.1725.35
Barramento de ramificação27.7628.0429.07
Tubo de cobre de barramento de ramificação28.4229.3130.91
Braço de saída inferior24.9524.8526.33

Tabela 6: Valores de temperatura nos pontos de monitoramento de cada fase do gabinete do anel sob o módulo de solução de campo de temperatura Icepak

Discussion

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Este artigo é uma análise de simulação comparativa do aumento de temperatura do gabinete do anel com base em software de modelagem de engenharia e software de elementos finitos, e a solução mais adequada para a situação real de aumento de temperatura é analisada por dois módulos de solução de campo de temperatura de elementos finitos. O gerenciamento térmico também é descrito no Icoz23 como um componente crítico e essencial para manter a alta eficiência e confiabilidade dos componentes eletrônicos. A importância de realizar uma análise comparativa é resumida com base no trabalho de Steiner24: uma análise comparativa foi realizada usando COMSOL e ANSYS Mechanical. Portanto, no processo de fabricação da RMU real, sua distribuição de temperatura pode ser analisada por software de simulação finita, o que pode economizar muito mão de obra e custos de produção.

O modelo tridimensional da RMU pode ser criado usando o SolidWorks, conforme mostrado na Figura 1A. Na Etapa 1.1, a simplificação do modelo é enfatizada como uma etapa crucial na análise de elementos finitos23. Dados os inúmeros componentes da RMU, que podem afetar a precisão da solução e os resultados computacionais, as peças desnecessárias são removidas após considerar apenas os componentes condutores. As principais peças de trabalho são mantidas, conforme ilustrado na Figura 1B.

No estágio de pré-processamento do campo de correntes parasitas, um aspecto crítico envolve a adição e verificação da excitação da carga. Conforme descrito na Etapa 2.1, é essencial adicionar com precisão a corrente de carga aos braços de saída superior e inferior, garantindo que todo o circuito condutor forme um caminho para a corrente. Isso facilita o cálculo do mapa de nuvens de distribuição do campo de correntes parasitas. A dificuldade reside na necessidade de verificar o circuito antes do cálculo da solução, pois um caminho de corrente incompleto no circuito pode levar à não convergência durante a resolução ou situações em que Maxwell não consegue calcular. A análise subsequente do campo de temperatura depende inteiramente das perdas obtidas com a resolução desse campo de correntes parasitas, tornando a inspeção do caminho da corrente uma etapa necessária.

Na solução de campo de temperatura, a Etapa 3 apresenta dois módulos diferentes para resolver o modelo. No entanto, garantir condições iniciais idênticas durante o processo de solução representa um desafio. Como os módulos têm ênfases diferentes em suas soluções, é necessário aproximar ou definir temperaturas ambientais, coeficientes convectivos e outras condições de solução idênticas para garantir a consistência sob a fonte de calor exclusiva fornecida pela Maxwell. Além disso, controlar os intervalos para solução de campo de temperatura em ambos os modelos para serem consistentes permite uma comparação direta das exibições de temperatura para os mesmos componentes na mesma faixa de temperatura, destacando as diferenças nas soluções e facilitando conclusões intuitivas.

Este trabalho apresenta um método inovador para medir e analisar com precisão a distribuição de temperatura de equipamentos elétricos durante a operação. Enfrentando o desafio de medir com precisão as perdas ôhmicas com métodos tradicionais, este estudo utiliza o solucionador de correntes parasitas Maxwell para cálculos precisos, que servem como base para resolver o campo de temperatura. Posteriormente, o módulo de resolução de campo de temperatura é usado para exibir visualmente a distribuição de temperatura de cada componente, aumentando significativamente a eficiência e a precisão dos testes de engenharia. A novidade deste artigo reside não apenas em fornecer um método eficiente para resolver o campo de temperatura, mas também em demonstrar como usar esses resultados para otimização estrutural e análise de dissipação de calor. Isso oferece novos meios técnicos para o projeto e otimização de equipamentos elétricos. No entanto, houve limitações no processo de pesquisa e, em trabalhos futuros, novas pesquisas serão realizadas nas seguintes áreas. Em primeiro lugar, refinamento do modelo e acoplamento de vários campos. O modelo será refinado para refletir com precisão as características térmicas do gabinete de rede em anel durante a operação real. Em segundo lugar, estudar campos de temperatura transientes. Pesquisas aprofundadas serão realizadas sobre as mudanças de campo de temperatura do gabinete de rede em anel sob condições de carga variadas para fornecer suporte adicional para a operação estável do gabinete de rede em anel em condições operacionais complexas.

Disclosures

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Os autores não têm interesses conflitantes.

Acknowledgements

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Os autores agradecem ao Sr. Wu, MS Sun, Sr. Wang, Sr. Mu e Sr. Li por sua ajuda. Este estudo foi apoiado pela China Postdoctoral Science Foundation (2022M721604) e pelo Wenzhou Key Science and Technology Tackling Programmer (ZG2023015).

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
Ar//Gases convencionais
Alumínio//Materiais de liga
Cobre//Materiais de liga
IcepakANSYS companyANSYS 2021R1A CFD software de simulação térmica
Hospedagem de PC/CPU Intel(R) Core(TM) i5-13500F de 12ª geraçãoEquipamento
de computador hostSolidWorksSubsidiária da Dassault SystèmesSolidWorks2021Uma ferramenta de desenho de software de engenharia
Empresa ANSYS térmica de estado estacionárioANSYS 2021R1Uma ferramenta de solução de simulação térmica

References

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Temperature rise test and analysis of high current switchgear in distribution system. J Engg. , 754-757 (2019).">Xia, H., et al. Temperature rise test and analysis of high current switchgear in distribution system. J Engg. , 754-757 (2019).
  2. Thermal effect on electric power network components under short-circuit currents. Electric Power Syst Res. 72 (3), 261-267 (2004).">Polykrati, A. D., Karagiannopoulos, C. G., Bourkas, P. D. Thermal effect on electric power network components under short-circuit currents. Electric Power Syst Res. 72 (3), 261-267 (2004).
  3. Multiphysics analysis of plug-in connector under steady and short circuit conditions. IEEE Trans Comp Packag Manu Technol. 5 (3), 320-327 (2015).">Guan, X., Shu, N., Kang, B., Zou, M. Multiphysics analysis of plug-in connector under steady and short circuit conditions. IEEE Trans Comp Packag Manu Technol. 5 (3), 320-327 (2015).
  4. Simulation analysis on the impact of different filling gases on the temperature rise of C-GIS. IEEE Trans Comp Packag Manu Technol. 9 (10), 2055-2065 (2019).">Wang, L., Wang, R., Li, X., Jia, S. Simulation analysis on the impact of different filling gases on the temperature rise of C-GIS. IEEE Trans Comp Packag Manu Technol. 9 (10), 2055-2065 (2019).
  5. Numerical design and optimization of a novel heatsink using ANSYS steady-state thermal analysis. 2020 27th International Workshop on Electric Drives: MPEI Department of Electric Drives 90th Anniversary (IWED. , 1-5 (2020).">Mueller, A., et al. Numerical design and optimization of a novel heatsink using ANSYS steady-state thermal analysis. 2020 27th International Workshop on Electric Drives: MPEI Department of Electric Drives 90th Anniversary (IWED. , 1-5 (2020).
  6. Optimizing GIS partial discharge pattern recognition in the ubiquitous power internet of things context: A MiNET deep learning model). Int J Electrical Power Energy Sys. 125, 106484(2021).">Wang, Y., Yan, J., Yang, Z., Zhao, Y., Liu, T. Optimizing GIS partial discharge pattern recognition in the ubiquitous power internet of things context: A MiNET deep learning model). Int J Electrical Power Energy Sys. 125, 106484(2021).
  7. A 3-D steady-state analysis of thermal behavior in. EHV GIS Busbar. J Electr Engg Tech. 11 (3), 781-789 (2016).">Lei, J., et al. A 3-D steady-state analysis of thermal behavior in. EHV GIS Busbar. J Electr Engg Tech. 11 (3), 781-789 (2016).
  8. Temperature rise simulation model of RMU with switchfuse combinations for future load profiles. CIRED 2021 - The 26th International Conference and Exhibition on Electricity Distribution. , 360-364 (2021).">Ouerdani, Y., et al. Temperature rise simulation model of RMU with switchfuse combinations for future load profiles. CIRED 2021 - The 26th International Conference and Exhibition on Electricity Distribution. , 360-364 (2021).
  9. Multi-physical field coupling simulation and thermal design of 10 kV-KYN28A high-current switchgear. Thermal Sci Engg Prog. 43, 101954(2021).">Zheng, W., Jia, X., Zhou, Z., Yang, J., Wang, Q. Multi-physical field coupling simulation and thermal design of 10 kV-KYN28A high-current switchgear. Thermal Sci Engg Prog. 43, 101954(2021).
  10. Electromagnetic-thermal-flow field coupling simulation of 12-kV medium-voltage switchgear. IEEE Trans Comp Packag Manufact Technol. 6 (8), 1208-1220 (2016).">Wang, L., et al. Electromagnetic-thermal-flow field coupling simulation of 12-kV medium-voltage switchgear. IEEE Trans Comp Packag Manufact Technol. 6 (8), 1208-1220 (2016).
  11. Thermal analysis and design of GaN device of energy storage converter based on Icepak. 2022 4th International Conference on Smart Power & Internet Energy Systems (SPIES. , 762-767 (2022).">Zhu, Y., et al. Thermal analysis and design of GaN device of energy storage converter based on Icepak. 2022 4th International Conference on Smart Power & Internet Energy Systems (SPIES. , 762-767 (2022).
  12. Thermal simulation of high-current switch cabinet based on Icepak). Electr Ener Mgmt Technol. , 1-7 (2018).">Ye Mao, Thermal simulation of high-current switch cabinet based on Icepak). Electr Ener Mgmt Technol. , 1-7 (2018).
  13. Evaluation of thermal transient and overload capability of high-voltage bushings with ATP. IEEE Trans Power Delivery. 24 (3), 1295-1301 (2009).">Zhang, S. Evaluation of thermal transient and overload capability of high-voltage bushings with ATP. IEEE Trans Power Delivery. 24 (3), 1295-1301 (2009).
  14. Steady-state thermal model of a synchronous reluctance motor. 2018 IEEE 59th International Scientific Conference on Power and Electrical Engineering of Riga Technical University (RTUCON. , 1-5 (2018).">Ghahfarokhi, P. S., et al. Steady-state thermal model of a synchronous reluctance motor. 2018 IEEE 59th International Scientific Conference on Power and Electrical Engineering of Riga Technical University (RTUCON. , 1-5 (2018).
  15. Temperature field and power loss calculation with coupled simulations for a medium-voltage simplified switchgear). Electrica. 23 (1), 107-120 (2021).">Şeker, E. A., Çelik, B., Yildirim, D., Sakaci, E. A., Deniz, A. Temperature field and power loss calculation with coupled simulations for a medium-voltage simplified switchgear). Electrica. 23 (1), 107-120 (2021).
  16. Research and application of temperature load of switchgear. J Physics: Conf Series. 2378 (2022), (2022).">Ruibo, Y., et al. Research and application of temperature load of switchgear. J Physics: Conf Series. 2378 (2022), (2022).
  17. Simulation for impact of nanofluid spectral splitter on efficiency of concentrated solar photovoltaic thermal system. Sust Cities Soc. 101, 105139(2024).">Sheikholeslami, M., Khalili, Z. Simulation for impact of nanofluid spectral splitter on efficiency of concentrated solar photovoltaic thermal system. Sust Cities Soc. 101, 105139(2024).
  18. Environmental and energy assessment of photovoltaic-thermal system combined with a reflector supported by nanofluid filter and a sustainable thermoelectric generator. J Cleaner Prod. 438, (2024).">Sheikholeslami, M., Khalili, Z., Scardi, P., Ataollahi, N. Environmental and energy assessment of photovoltaic-thermal system combined with a reflector supported by nanofluid filter and a sustainable thermoelectric generator. J Cleaner Prod. 438, (2024).
  19. Solar photovoltaic-thermal system with novel design of tube containing eco-friendly nanofluid. Renewable Ener. , 222(2024).">Sheikholeslami, M., Khalili, Z. Solar photovoltaic-thermal system with novel design of tube containing eco-friendly nanofluid. Renewable Ener. , 222(2024).
  20. Environmental and energy analysis for photovoltaic-thermoelectric solar unit in existence of nanofluid cooling reporting CO2 emission reduction. J Taiwan Inst Chem Eng. 156, 105341(2024).">Sheikholeslami, M., Khalili, Z. Environmental and energy analysis for photovoltaic-thermoelectric solar unit in existence of nanofluid cooling reporting CO2 emission reduction. J Taiwan Inst Chem Eng. 156, 105341(2024).
  21. Research on the temperature rise characteristics of medium-voltage switchgear under different operation conditions. IEEJ Trans Elect Electr Engg. 17 (5), 654-664 (2022).">Zhao, L., et al. Research on the temperature rise characteristics of medium-voltage switchgear under different operation conditions. IEEJ Trans Elect Electr Engg. 17 (5), 654-664 (2022).
  22. Influence of heat source location on air temperatures in sealed MV switchgear. CIRED - 24th International Conference on Electricity Distribution. , 1-5 (2017).">Fjeld, E., Rondeel, W., Vaagsaether, K., Attar, E. Influence of heat source location on air temperatures in sealed MV switchgear. CIRED - 24th International Conference on Electricity Distribution. , 1-5 (2017).
  23. Light weight high performance thermal management with advanced heat sinks and extended surfaces. IEEE Trans Comp Pack Technol. 33 (1), 161-166 (2010).">Icoz, T., Arik, M. Light weight high performance thermal management with advanced heat sinks and extended surfaces. IEEE Trans Comp Pack Technol. 33 (1), 161-166 (2010).
  24. High temperature steady-state experiment for computational radiative heat transfer validation using COMSOL. and ANSYS. Results Engg. 13, 100354(2022).">Steiner, T. R. High temperature steady-state experiment for computational radiative heat transfer validation using COMSOL. and ANSYS. Results Engg. 13, 100354(2022).

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

Ring Main UnitTemperature RiseFinite Element SimulationOhmic LossesTemperature Field AnalysisConvective Heat TransferSteady State ThermalMesh GenerationNatural ConvectionElectrical Equipment

Related Articles