Questo articolo affronta il problema dell'aumento di temperatura dell'unità principale dell'anello stabilendo un modello semplificato e conducendo un'analisi comparativa in due moduli di risoluzione del campo di temperatura.
Articolo metodologico
Studio comparativo della simulazione dell'aumento di temperatura nell'unità principale dell'anello
In questo articolo
Sommario
Abstract
L'unità principale ad anello (RMU) è un dispositivo critico nei sistemi di distribuzione dell'energia utilizzati per il collegamento e la distribuzione dell'elettricità. Tuttavia, a causa della sua struttura interna compatta e dell'elevato carico di corrente, i problemi di dissipazione del calore sono particolarmente importanti. Per affrontare questo problema, questo studio propone in modo innovativo un modello RMU semplificato, che impiega metodi di simulazione agli elementi finiti per risolvere con precisione le perdite ohmiche dei conduttori in condizioni operative reali e ottenere dati sulle perdite ohmiche per vari componenti. Questa è la prima indagine approfondita sul problema dell'aumento della temperatura dell'RMU utilizzando un approccio così completo. Successivamente, il campo di temperatura è stato risolto utilizzando due diversi moduli di analisi del campo di temperatura, con un confronto dettagliato e un'analisi dei risultati della simulazione per identificare somiglianze, differenze e tendenze nella distribuzione della temperatura. I risultati indicano che il modello di soluzione del campo di temperatura, che considera il trasferimento di calore convettivo, è più accurato e si allinea con le condizioni operative effettive. Questa ricerca fornisce un approccio innovativo e soluzioni pratiche per la progettazione e l'ottimizzazione delle RMU. La ricerca futura può esplorare ulteriormente i metodi di analisi dell'accoppiamento multifisico per affrontare la progettazione strutturale e i problemi di convalida obbligatoria per RMU ad alta e altissima tensione e altre apparecchiature elettriche, fornendo così importanti informazioni per la progettazione ingegneristica.
Introduzione
L'unità principale ad anello è un gruppo di quadri ad alta tensione montati in un armadio metallico in acciaio o costituiti da un'unità di alimentazione di rete ad anello distanziata assemblata di apparecchiature elettriche. La struttura complessiva dell'interruttore di carico e del circuito conduttivo è costituita dal circuito conduttivo, che comprende una serie di componenti che costituiscono il nucleo principale dell'unità ad anello. Tuttavia, a causa della sua struttura interna compatta, l'unità principale ad anello deve affrontare sfide nella dissipazione del calore. Ciò può portare a deformazione termica e invecchiamento quando si opera per lunghi periodi in ambienti ad alta temperatura. Questi problemi non solo influiscono sulla durata dell'unità, ma influiscono anche sulle sue proprietà isolanti, ponendo rischi per la sicurezza. In particolare, i danni alle apparecchiature e gli incidenti elettrici diventano più probabili, ponendo notevoli rischi per la sicurezza.
All'interno di diverse aree di ricerca, gli studiosi hanno condotto una serie di studi sull'aumento della temperatura dei quadri delle linee aeree e hanno analizzato vari fattori che influenzano la distribuzione della temperatura1. In Polykrati et al.2, viene presentato un modello matematico per la stima dell'aumento di temperatura dei componenti installati sulla rete di distribuzione durante un guasto di cortocircuito. Il modello è stato applicato ai sezionatori comuni della rete e le caratteristiche dei risultati sono state tracciate in base alle diverse forme della parte asimmetrica della forma d'onda della corrente di cortocircuito e al valore iniziale della componente di corrente CC in cortocircuito. Guan et al., d'altra parte, hanno preso in considerazione la resistenza di contatto e la repulsione elettromagnetica costruendo un ponte di contatto equivalente per simulare l'interfaccia di contatto e hanno ulteriormente analizzato il campo di accoppiamento elettromagnetico-termico e l'esperimento di aumento della temperatura3. Inoltre, i ricercatori hanno studiato il campo di temperatura e la distribuzione delle sollecitazioni termiche dei contatti dinamici e statici all'interno dell'unità principale ad anello mediante simulazione agli elementi finiti, che ha fornito una base per lo studio della durata dell'interruttoreautomatico 4. Infine, Mueller et al. si sono concentrati sulle caratteristiche geometriche dei dissipatori di calore e hanno valutato gli effetti della selezione del materiale, dell'area superficiale totale, dell'uniformità della temperatura e della temperatura superficiale massima sulle prestazioni termiche5. Questi studi forniscono preziose informazioni e metodi per migliorare le prestazioni e l'affidabilità dei quadri elettrici, ridurre l'aumento della temperatura e prolungare la durata delle apparecchiature. Wang et al. hanno proposto un modello di apprendimento profondo MiNET (MDLM) nell'ambiente UPIOT con lo scopo di rilevare la diagnosi dei guasti degli armadi elettrici ad anello, che è stato convalidato per avere una precisione di identificazione del 99,1%, che è significativamente superiore a quella di altri metodi6. Lei et al. hanno studiato le prestazioni termiche di una sbarra GIS in uno stato stazionario utilizzando il metodo di analisi dell'accoppiamento magneto-fluido-termico, ottimizzando così il diametro del conduttore e del serbatoio in base ai risultati della simulazione dell'aumento di temperatura7. Ouerdani et al. hanno utilizzato il modello di simulazione dell'aumento della temperatura dell'RMU per determinare l'aumento della temperatura in punti critici al suo interno, fissando così la durata del sovraccarico massimo per i componenti all'interno dell'RMU di conseguenza8. Zheng et al. hanno descritto una sbarra rettangolare convenzionale in un modello di quadro ad alta corrente costruendo un modello bidimensionale e applicando il metodo degli elementi finiti (FEM) per i calcoli dei campi elettromagnetici. Ciò ha permesso loro di ottenere la distribuzione della densità di corrente e della perdita di potenza del conduttore del bus. Una sbarra irregolare è stata progettata dopo aver considerato gli effetti dell'effetto di prossimità e dell'effetto pelle. Questo design irregolare delle sbarre ha migliorato le prestazioni delle tradizionali sbarre rettangolari9.
Per quanto riguarda l'aspetto dell'utilizzo della simulazione icepak, Wang et al. hanno effettuato una simulazione dell'aumento di temperatura attraverso le teorie del campo del vortice, del campo del flusso d'aria e del campo di temperatura e hanno scoperto che l'aumento della temperatura dell'unità principale dell'anello era più grave sotto convezione naturale. Hanno ridotto con successo il livello di aumento della temperatura aggiungendo un raffreddamento ad aria forzata e apportando miglioramenti alla struttura di contatto interna10. Zhu et al.11 hanno utilizzato l'icepak per simulare un modello termico al fine di confrontare l'effetto della presenza di vie termiche sul PCB e la presenza di dissipatori di calore sulla temperatura dei dispositivi di potenza. Infine, l'analisi teorica viene confrontata con i risultati della simulazione per verificare la correttezza dell'analisi teorica. Mao et al.12 hanno studiato la temperatura e la distribuzione del flusso d'aria interno in condizioni operative estive mediante simulazione termica basata sul software CAE nella simulazione icepak. Viene presentato il problema di come migliorare l'efficienza di raffreddamento e controllare l'aumento di temperatura di più contatti argentati e i contorni della temperatura e del flusso d'aria interno catturati nella simulazione getteranno le basi per la progettazione dello schema di raffreddamento per i sei contatti argentati montati nell'unità di tenuta. Al contrario, nell'uso di un modulo termico allo stato stazionario, vengono discussi i metodi di modellazione di Zhang13 per risolvere la rete termica di una boccola ad alta pressione utilizzando una procedura transitoria alternativa. I risultati dei test e delle simulazioni sono in buon accordo con lo stato stazionario termico e gli stati transitori della boccola. I risultati transitori vengono quindi utilizzati per valutare la capacità di sovraccarico della boccola. Vaimann et al.14 hanno sviluppato e analizzato un modello termico analitico di un motore sincrono a riluttanza per prevedere la temperatura dei suoi diversi componenti e la rete termica a parametro totale impostato.
Con il continuo progresso della ricerca sulle apparecchiature elettriche come le unità principali ad anello, i test convenzionali di aumento della temperatura e i metodi di produzione sono relativamente inefficienti. Pertanto, utilizzando la tecnologia degli elementi finiti combinata con i test offline, non solo vengono affrontati i problemi dei costi di progettazione, ma è possibile apportare prontamente modifiche e ottimizzazioni ai problemi del mondo reale sulla base di simulazioni. Sulla base dei progressi della ricerca sopra menzionati, l'uso di ANSYS Icepak e dell'accoppiamento termico allo stato stazionario per l'analisi comparativa è raramente menzionato. Pertanto, il protocollo descrive la ricerca dei meccanismi degli elementi finiti, utilizza combinazioni numeriche e morfologiche per stabilire un modello di simulazione dell'aumento della temperatura degli elementi finiti per l'involucro e discute il modello di simulazione dell'aumento della temperatura degli elementi finiti basato sui risultati dei due moduli analitici confrontando i risultati dei due moduli di simulazione. Attraverso il confronto tra i due moduli di simulazione, otterremo le caratteristiche dell'andamento dell'aumento di temperatura dell'unità principale ad anello e troveremo il metodo più applicabile in modo da fornire le basi necessarie e le idee di ricerca per una strategia per mitigare l'aumento di temperatura dell'unità principale ad anello.
Protocollo
1. Modello
NOTA: A causa della complessa struttura dell'unità principale ad anello (Figura 1A), è stato scelto un software di progettazione online per semplificare il funzionamento dell'unità principale ad anello.
- Semplificazione della modellazione
- Semplifica parzialmente il modello, preservando la sezione dell'air box dell'RMU e rimuovendo o semplificando altri componenti come alberi isolanti, bulloni di fissaggio, dadi, componenti di tenuta e staffe di supporto della pressione. La versione semplificata è mostrata in (Figura 1B).
- Nel processo di semplificazione dell'unità principale ad anello di tipo 630A, rimuovere l'albero isolato che collega la stanza dell'interruttore automatico con la scatola degli strumenti e molti bulloni e dadi fissi. Estrarre le parti di tenuta e la staffa di mantenimento della pressione e collegare i contatti statici della trave statica isolata con la sbarra di derivazione inferiore con la premessa di garantire che l'intera configurazione abbia la stessa corrente conduttiva e che vengano mantenuti solo l'interruttore a vuoto, la piastra di fissaggio dell'interruttore automatico e i contatti statici e l'interruttore a vuoto.
- Conservare solo l'interruttore a vuoto, la piastra di fissaggio dell'interruttore automatico, il contatto statico e la piastra di blocco dell'interruttore a vuoto. Nel complesso, rimuovere i bulloni e le guarnizioni dal modello, riempire i fori dopo aver rimosso i bulloni con solidi, ridurre il numero di parti mesh e ottimizzare le forme irregolari delle parti. Rimuovere gli strumenti per il funzionamento del pannello, le piastre di montaggio, le staffe e altre parti operative, come le scatole degli strumenti, che non hanno alcun effetto sul processo di simulazione dell'aumento della temperatura.
- La rimozione degli alloggiamenti isolati di alcuni componenti può essere ignorata durante la simulazione in quanto hanno scarso effetto sui risultati della simulazione. Inoltre, gli interruttori di messa a terra che non hanno alcun effetto sull'uso dell'apparecchiatura durante il normale funzionamento li rimuovono e mantengono la stanza dell'interruttore automatico per la simulazione.
- Per eliminare qualsiasi sezione, è sufficiente selezionarla e fare clic sull'opzione Elimina .
- Semplifica parzialmente il modello, preservando la sezione dell'air box dell'RMU e rimuovendo o semplificando altri componenti come alberi isolanti, bulloni di fissaggio, dadi, componenti di tenuta e staffe di supporto della pressione. La versione semplificata è mostrata in (Figura 1B).
2. Soluzione in campo parassito
- Impostazioni di pre-elaborazione
NOTA: L'emulazione del campo delle correnti parassite è la base per l'esecuzione della soluzione del campo di temperatura, che richiede la successiva analisi della fonte di calore risolta come carico sul campo di temperatura.- Fare riferimento alla documentazione dell'apparecchiatura dell'unità principale ad anello e ai relativi manuali per raccogliere informazioni sulle proprietà fisiche e sui parametri di ciascun componente dell'unità principale ad anello. Impostare gli attributi fisici e i parametri dei componenti dell'unità principale dell'anello in Maxwell in base alle informazioni ottenute, come descritto in dettaglio nella Tabella 1.
- Impostare la corrente di carico del campo a correnti parassite a 630 A con una frequenza di 50 Hz. Nel software Maxwell, selezionare un lato dei bracci in uscita superiore e inferiore, accedere al modulo di eccitazione e impostare l'ampiezza della corrente su 630 A. Nella sezione delle impostazioni della soluzione, scegli una frequenza di 50 Hz.
NOTA: Nel circuito conduttivo di un'unità principale ad anello, il percorso formato da tutti i componenti dal braccio di uscita superiore al braccio di uscita inferiore è noto come sequenza di fase. Pertanto, in questo documento, le fasi A, B e C sono disposte da sinistra a destra. - I parametri del materiale dei componenti dell'unità principale ad anello sono mostrati nella Tabella 2.
- Dirigere la corrente attraverso i bracci di linea in uscita, i collegamenti flessibili, le sbarre, gli interruttori automatici, le sbarre di supporto dei contatti statici e le sbarre di derivazione per ogni fase. Lo scopo è quello di realizzare un percorso di corrente che permetta ai componenti di completare il carico.
- Utilizza la mesh adattiva di Maxwell per completare il controllo della griglia per il modello. Utilizzare il metodo di partizionamento adattivo della mesh Maxwell per i componenti più grandi e il perfezionamento della mesh locale per i componenti interni più piccoli.
NOTA: Maxwell è in grado di migliorare continuamente la precisione della griglia durante il processo di risoluzione, eliminando la necessità di fare clic su Operazioni mesh per ulteriori partizioni mesh. - Impostare le dimensioni del passaggio della soluzione. Fare clic su Analisi nell'albero del modello, aprire le impostazioni del passo di risoluzione e impostare il numero massimo di passate su 10. Mantieni le altre impostazioni ai valori predefiniti senza apportare modifiche.
- Principio del calcolo del campo delle correnti parassite15,16.
- Usa la prima equazione di Maxwell, che descrive l'azione della carica sulla generazione di un campo elettrico17.
(1)
dove ρ rappresenta la densità di carica; ε0 rappresenta la costante dielettrica del vuoto. - Usa la seconda equazione di Maxwell, che descrive la relazione tra un campo magnetico variabile e un campo elettrico e l'effetto di un campo magnetico sul moto di una carica.
(2)
dove
rappresenta l'intensità del campo magnetico. Questa equazione descrive che un campo magnetico variabile produce un campo elettrico a vortice, cioè lo spin del campo elettrico a vortice è uguale al negativo della velocità di variazione del campo magnetico nel tempo. - Usa la terza equazione di Maxwell, che descrive l'effetto della carica magnetica sulla produzione di un campo magnetico.
(3)
Questa equazione descrive il campo magnetico prodotto da una carica magnetica come passivo, cioè non ci sono monopoli nel campo magnetico. - Usa la quarta equazione di Maxwell, che descrive la relazione tra un campo elettrico variabile e un campo magnetico e l'effetto di una corrente elettrica su un campo magnetico.
(4)
dove
rappresenta la densità di corrente e μ0 rappresenta la permeabilità al vuoto. Questa equazione descrive che un campo elettrico variabile produce un campo magnetico a vortice, cioè lo spin del campo magnetico a vortice è uguale alla somma della densità di corrente e della velocità di variazione del campo elettrico nel tempo. - Sulla base delle equazioni di cui sopra, utilizzare il modulo risolutore Maxwell 3D utilizzando correnti parassite per risolvere la perdita ohmica generata dal circuito conduttivo nell'RMU, che fornisce una fonte di calore per la successiva analisi della simulazione termica. La sua espressione matematica è data come18
(5)
dove σ denota la conduttività del materiale dell'anello conduttore; J è la densità di corrente nel circuito.
- Usa la prima equazione di Maxwell, che descrive l'azione della carica sulla generazione di un campo elettrico17.
- Risultati del calcolo
- Fare clic sull'opzione Maxwell 3D nell'interfaccia e aprire il controllo di convalida per rivedere tutte le impostazioni alla ricerca di errori. Se non ci sono errori, procedi facendo clic su Analizza tutto per avviare il processo di risoluzione.
- Utilizzare il calcolatore di post-elaborazione di Maxwell per calcolare e tracciare le perdite ohmiche nel campo delle correnti parassite dell'unità principale dell'anello, come mostrato nella Tabella 3.
3. Soluzione del campo di temperatura
NOTA: A scopo comparativo, dividere il campo di temperatura in Icepak e termica allo stato stazionario. Imposta e risolvi ciascuno separatamente per ottenere un'analisi comparativa.
- Configurazione del modello Icepak
- Impostare le proprietà del materiale come segue: designare tutti i materiali solidi del circuito come Cu-Pure, con superfici che utilizzano Cu-polished-surface. Per i componenti del pannello, selezionare il materiale Aluminium6061-T6 , con un rivestimento superficiale di superficie Paint-AL con un'emissività di 0,35. Vedere la Tabella 4 per i dettagli. Fare clic con il pulsante destro del mouse sul componente selezionato, fare clic su Modifica e quindi andare su Proprietà per impostare il materiale sia per la superficie che per i materiali solidi.
- Seleziona il modello e fai clic su Imposta nel menu Modifica, quindi scegli Livello di mesh multilivello per regolare le impostazioni della mesh. Impostare l'armadio esterno su un livello di mesh pari a 2 e tutti i limiti su un livello di mesh pari a 2. Per tutti gli altri componenti, impostare il livello di mesh su 3. Infine, apri il controllo Mesh e fai clic su Genera per creare la mesh.
- Per garantire l'accuratezza e l'efficienza della simulazione, indipendentemente dalle dimensioni della rete, è necessaria la convalida dell'indipendenza dalla rete. Importa il modello geometrico dell'involucro del campo di temperatura, stabilito utilizzando il software di progettazione per la mesh.
- Come illustrato nella Figura 2, le curve di polarizzazione dei quattro insiemi di griglie sono ben allineate. A una tensione di lavoro di 0,5 V, le densità di corrente per i quattro set di reti sono rispettivamente di 2,357 A/cm2, 2,358 A/cm2, 2,356 A/cm2 e 2,454 A/cm2, con un errore tra la densità massima e quella di corrente inferiore all'1%. Per bilanciare efficienza e precisione, determinare la dimensione della griglia che viene 987924.
- Configurazione della soluzione
- Impostare le indicazioni del dominio della soluzione Cabinet su Opening.
- Nel software, selezionare Passaggio problema. In Parametri di base, controllare il modello di radiazione superficie-superficie, scegliere Equazione zero per Regime di flusso turbolento, selezionare l'opzione Gravità per Convezione naturale e impostare la temperatura ambiente su 20 °C.
- Nelle impostazioni del file, scegli Perdite di calore volumetriche per la mappatura EM e seleziona Tutti gli oggetti mostrati per completare le impostazioni di perdita.
- Calcolo del campo di temperatura
- In icepak, applicare tre principali equazioni di conservazione dell'energia: equazione di conservazione della massa, equazione di conservazione della quantità di moto ed equazione di conservazione dell'energia. In particolare, utilizzare l'equazione di conservazione della quantità di moto, che è la seguente19:
(6)
Equazione di conservazione dell'energia:
(7)
Equazioni di conservazione della massa:
(8)
Equazione del trasferimento di energia per il trasferimento di calore da una fonte di calore solida:
(9)
ρ rappresenta la densità del fluido; v rappresenta il vettore della velocità del flusso; T rappresenta la temperatura; p è la pressione; τ è la forza viscosa sulla superficie del micrometabolita; κ è il coefficiente di scambio termico; Shè la fonte di calore corporeo; h è l'entalpia specifica del fluido e F è la forza corporea del micrometabolita.
NOTA: I risultati dei calcoli del campo di temperatura sono mostrati nella Figura 3A e nella Figura 4A.
- In icepak, applicare tre principali equazioni di conservazione dell'energia: equazione di conservazione della massa, equazione di conservazione della quantità di moto ed equazione di conservazione dell'energia. In particolare, utilizzare l'equazione di conservazione della quantità di moto, che è la seguente19:
- Configurazione del modello termico in stato stazionario
- Mantenere le proprietà del materiale come indicato nella Tabella 3 nelle impostazioni del materiale. Generare le perdite ohmiche risultanti dall'analisi della simulazione del campo delle correnti parassite nel modulo termico a stato stazionario facendo clic su Generazione carico termico .
- Fare clic sul valore della temperatura convettiva e impostarlo su 20 °C, con un coefficiente convettivo di 5 (W/m²°C) applicato alle pareti interne dell'armadio, dei componenti e dell'armadio esterno. Applicare le impostazioni e generare. Imposta l'output per risolvere la temperatura facendo clic su Risolvi > Risultati output.
NOTA: Il principale campo di temperatura che governa le equazioni nel principio di calcolo del campo di temperatura termicaallo stato stazionario 20,21,22 è solitamente derivato dalla legge di conduzione del calore (legge di Fourier della conduzione del calore). Nel caso unidimensionale, l'equazione di trasferimento del calore del campo di temperatura può essere espressa come20:
(10)
In questa equazione, T rappresenta la temperatura all'interno dell'oggetto, t è il tempo, x sono le coordinate spaziali e α è la diffusività termica. Questa equazione descrive la variazione della temperatura rispetto al tempo e allo spazio, dove il lato destro esprime la relazione tra la velocità di conduzione del calore e il gradiente di temperatura. In uno scenario tridimensionale più generale, l'equazione di conduzione del calore per il campo di temperatura può essere espressa nella seguente forma:
(11)
ρ rappresenta la densità dell'oggetto, c è la capacità termica specifica, K è la conducibilità termica e Q è il termine della fonte di calore all'interno del volume. Questa equazione descrive la variazione del campo di temperatura, influenzata dalla conduzione del calore, dalle fonti di calore e dalla capacità termica. - I risultati del calcolo del campo di temperatura sono mostrati nella Figura 3. Confrontare i valori di temperatura riassunti nella Tabella 5 e nella Tabella 6.
Risultati
Sulla base dei dati della Tabella 3, si possono trarre le seguenti conclusioni: Le perdite complessive per le fasi A, B e C sono relativamente simili. In particolare, le perdite totali per la Fase A sono di 16,063 W/m³, la Fase B è di 16,12 W/m³ e la Fase C è di 19,57 W/m³. Le posizioni con perdite più elevate possono essere in corrispondenza delle connessioni di vari componenti. Ciò è dovuto principalmente al fatto che in questi punti di connessione esistono tipicamente la resistenza di contatto e la resistenza del conduttore. Quando la corrente passa attraverso queste connessioni, viene generato un calore significativo, che porta ad un aumento della temperatura e a maggiori perdite in quelle aree.
I bracci uscenti superiore e inferiore dell'unità principale ad anello sopportano alcune perdite, soprattutto quando si trasporta il carico principale. Le perdite nei bracci uscenti delle tre fasi sono approssimativamente le stesse. Questo perché questa parte della corrente è relativamente concentrata e, a causa della forma regolare, il valore di resistenza è piccolo. Pertanto, le perdite in queste parti sono approssimativamente uguali.
Le perdite nella sbarra di derivazione sono relativamente elevate, principalmente a causa delle sue numerose curve e della presenza di sezioni angolate. La maggior parte della corrente è concentrata nella zona di curvatura e in prossimità degli angoli. Oltre alla sezione dell'interruttore a vuoto, anche le perdite nella parte in tubo di rame della sbarra sono relativamente elevate, per un totale di 22,32 W/m³, pari al 40% delle perdite ohmiche complessive delle tre fasi. Le perdite nella parte di contatto statico sono relativamente piccole rispetto alle perdite complessive.
Altri componenti, come la piastra fissa dell'interruttore, il deflettore dell'interruttore e il guscio esterno dell'unità principale ad anello, hanno perdite minori. Poiché non partecipano direttamente al carico, le loro perdite sono generate principalmente dai componenti interni per conduzione, con conseguenti valori di perdita più piccoli. Nei calcoli di post-elaborazione, le perdite di questi componenti non sono particolarmente dettagliate. In sintesi, sono state delineate le principali caratteristiche della distribuzione delle perdite, fornendo preziose informazioni per un'ulteriore ottimizzazione del design.
Combinando i risultati della Figura 3 e della Figura 4, insieme alla Tabella 5 e alla Tabella 6, è stato condotto un confronto approfondito dei risultati della soluzione del campo di temperatura tra i moduli Icepak e Steady-state Thermal. Nel processo di analisi, sono state osservate differenze e somiglianze tra i due moduli.
Nell'intervallo di temperatura compreso tra 20 e 31,24 °C, come mostrato nella Figura 2A, le aree di temperatura più elevate sono concentrate nelle regioni in cui l'involucro RMU entra in contatto con i componenti interni. Il motivo principale è che queste aree di contatto fungono tipicamente da percorsi critici per la conduzione del calore. Durante il funzionamento, i componenti interni generano calore, che viene condotto all'involucro attraverso le superfici di contatto con una buona conduttività termica, provocando così un aumento della temperatura in queste regioni. Queste aree includono i punti di contatto del braccio di linea uscente superiore e le porzioni del braccio di linea uscente inferiore a contatto con il guscio vicino, con un intervallo di temperatura di circa 24-27 °C. Rispetto alla Figura 3B, l'intervallo di riscaldamento delle conchiglie in entrambi i casi è più o meno lo stesso e le variazioni di temperatura si propagano dai punti caldi alle aree più fredde. Esaminando la Figura 4A, B, anche le tendenze generali della temperatura in entrambi i casi sono approssimativamente simili. La fonte primaria di aumento della temperatura è la conduzione del calore attraverso il braccio di linea uscente inferiore, che passa attraverso i conduttori interni e si dissipa dal braccio di linea uscente superiore. La tendenza al riscaldamento è concentrata nelle aree vicine al braccio inferiore della linea in uscita, comprese le sbarre di derivazione di ciascuna fase e le sbarre di supporto collegate. Inoltre, dall'osservazione della distribuzione della temperatura nei circuiti interni nella Figura 4, sia che si utilizzi Icepak o il modulo termico allo stato stazionario per la risoluzione, la temperatura complessiva della Fase B è costantemente più alta rispetto alle altre due fasi. Ciò suggerisce principalmente che la Fase B, durante il processo di carico, non solo sopporta il calore generato dalla corrente di quella fase ma anche, a causa della struttura interna eccessivamente compatta dell'RMU, il calore prodotto da ciascuna fase non può dissipare tempestivamente. A causa dello scambio di calore, la temperatura della Fase B rimane più alta rispetto alle altre due fasi indipendentemente dal solutore utilizzato. Anche l'andamento generale della temperatura mostra un alto grado di somiglianza, con l'aumento principale della temperatura derivante dalla conduzione del calore attraverso il braccio di linea in uscita inferiore e dalla dissipazione del calore dal braccio di linea in uscita superiore. Questa tendenza sembra molto coerente tra i due moduli. Inoltre, le posizioni degli hotspot nelle mappe di distribuzione della temperatura di entrambi i moduli sono altamente coerenti. In particolare, nella sezione delle sbarre di derivazione della Fase C, entrambi i solutori mostrano la stessa temperatura massima e le differenze sono quasi trascurabili. Ciò implica che, indipendentemente dal solutore utilizzato, l'identificazione accurata delle posizioni di temperatura più calde è fondamentale per la gestione termica dell'RMU.
In secondo luogo, dalla tabella si può vedere che la stessa corrente produce la stessa perdita, ma il valore di temperatura di ciascun componente di fase nella Tabella 6 è generalmente inferiore a quello della Tabella 5; ad esempio, il valore della temperatura del bus di derivazione nell'area del punto caldo dell'intera RMU è la temperatura più alta nella Tabella 6 a 30,91 °C e, d'altra parte, la temperatura più alta del bus di derivazione nella Tabella 5 è 31,24 °C. Il motivo è la logica della soluzione di Icepak: l'RMU come fonte di calore continuerà a scambiare calore e dissipare calore, e quando un fluido viene impostato come mezzo per il trasferimento di calore, la temperatura nel dominio della soluzione verrà gradualmente emessa attraverso il mezzo, il che si traduce in una diminuzione della temperatura. Il modulo termico allo stato stazionario pone meno enfasi sulla convezione del calore con l'aria circostante e si concentra invece su una soluzione basata sulla conduzione del calore. Rispetto ai modelli che tengono conto della convezione, questo approccio si traduce in una soluzione di campo di temperatura meno completa. Di conseguenza, nelle zone in cui la temperatura è più elevata, i punti caldi diventano più pronunciati. Nel processo del test di aumento della temperatura effettivo, oltre a considerare le caratteristiche del modello stesso, il trasferimento di calore considera allo stesso tempo anche l'influenza del mezzo di dissipazione del calore, come l'aria, ecc., sulla temperatura complessiva. In combinazione con le esigenze pratiche e le considerazioni sul funzionamento sperimentale, l'RMU come fonte di calore continua a svolgere la convezione e la dissipazione del calore; La simulazione della temperatura di IcePak è più in linea con le reali esigenze. Nel caso in cui l'unità principale ad anello venga utilizzata come fonte di calore con convezione e dissipazione continua del calore, la simulazione della temperatura di icepak è più in linea con le esigenze pratiche. Al contrario, il modulo termico stazionario si concentra principalmente sulla conduzione del calore, che in alcuni casi potrebbe non essere in grado di soddisfare le effettive esigenze.
Figura 1: Modello dell'unità principale ad anello. (A) Modello generale dell'unità principale ad anello (B) Modello semplificato dell'unità principale ad anello. Clicca qui per visualizzare una versione più grande di questa figura.
Figura 2: Curve di polarizzazione per quattro serie di mesh. La mesh adattiva Maxwell viene utilizzata per le parti più grandi e la finitura della mesh locale viene utilizzata per le parti più piccole. Nel software di progettazione, viene creato e importato in Maxwell un modello 3D semplificato dell'armadio principale ad anello e il modello viene meshato utilizzando il modulo Mesh. La figura rappresenta la verifica delle griglie del modello. La densità di corrente dei quattro gruppi di griglie è 2,357 A/cm2, 2,358 A/cm2, 2,356 A/cm2 e 2,354 A/cm2 sotto il carico di tensione e l'errore relativo tra la densità di corrente massima e la densità di corrente minima è inferiore all'1% e per tenere conto dell'efficienza e dell'accuratezza dei calcoli, Il numero finale delle griglie è determinato in 1169091. Clicca qui per visualizzare una versione più grande di questa figura.
Figura 3: Modello di analisi del campo di temperatura dell'armadio esterno. (A) Distribuzione della temperatura del guscio con soluzione Icepak. (B) Distribuzione della temperatura del guscio con soluzione termica allo stato stazionario. Clicca qui per visualizzare una versione più grande di questa figura.
Figura 4: Modello di analisi del campo di temperatura del circuito interno. (A) Distribuzione della temperatura dell'anello conduttore con soluzione Icepak (B) Distribuzione della temperatura dell'anello conduttore con soluzione termica allo stato stazionario. Clicca qui per visualizzare una versione più grande di questa figura.
| Materiale | Aria | Lamiera d'acciaio zincato | Rame | Alluminio | ||
| Capacità termica specifica (J/(kg· K)) | 1007 | 500 | 500 | 897 | ||
| Densità (kg/m³) | 1.1614 | 8030 | 8900 | 2689 | ||
| Permeabilità relativa | 1 | 0.3 | 0.9999991 | / | ||
| Permeabilità magnetica relativa | 1 | 2500 | 1 | 1 | ||
| Conducibilità termica (W/(m·K)) | 0.026 | 16 | 386 | 237 | ||
| Emissività | / | 0.65 | 0.3 | 0.1 | ||
| Conducibilità elettrica (S/m) | / | 0.8 | 5.80E+07 | / | ||
Tabella 1: Parametri fisici di alcuni materiali nell'unità principale dell'anello.
| Componenti | Permeabilità relativa | Conducibilità/s | Materiale |
| Circuiti che trasportano corrente (sbarre, sbarre derivate, ecc.) | 0 | 3,00,00,000 | Lega di alluminio |
| Interruttore automatico in vuoto | 0.9999 | 2,00,00,000 | Leghe di rame |
| Corpo | 200 | 11,00,000 | Acciaio strutturale |
Tabella 2: Elenco dei materiali per ogni componente.
| Perdita ohmica W/m³ | Un | B | C |
| Braccio di uscita | 3.78 | 3.72 | 3.73 |
| Sbarre di derivazione | 2.1 | 2.09 | 2.1 |
| Collegamento flessibile | 1.3 | 1.3 | 1.3 |
| Interruttore automatico in vuoto | 1.023 | 0.95 | 0.98 |
| Contatti statici | 0.36 | 0.36 | 0.36 |
| sbarre di alimentazione | 1.33 | 1.35 | 1.32 |
| Tubi di rame per sbarre di derivazione | 6.19 | 6.35 | 9.78 |
Tabella 3: Valori di perdita ohmica per i componenti di fase A, B e C.
| Numero | 1 | 2 |
| Materia | in alluminio6061-T6 | Cu-Puro |
| Conducibilità termica (W/m·K) | 167 | 387.6 |
| Densità (kg/m³) | 8933 | 2700 |
| Capacità termica specifica/(kg· K) | / | 896 |
| Materiale di superficie | Superficie Paint-AL | Cu-Superficie lucidata |
| Emissività | 0.35 | 0.052 |
Tabella 4: Impostazioni dei parametri del materiale.
| Punto di monitoraggio della temperatura/temperatura (°C) | Un | B | C |
| Braccio di uscita superiore | 25.48 | 25.79 | 25.63 |
| Sbarra principale | 25.93 | 26.28 | 26.13 |
| Sbarre di derivazione | 29 | 29.18 | 30.01 |
| Tubo di rame per sbarre di derivazione | 31.04 | 31.18 | 31.24 |
| Braccio di uscita inferiore | 26.5 | 26.98 | 26.92 |
Tabella 5: Valori di temperatura di ciascuna fase dell'unità principale ad anello nel modulo di risoluzione del campo di temperatura termica stazionaria.
| Punto di monitoraggio della temperatura/temperatura (°C) | Un | B | C |
| Braccio di uscita superiore | 23.73 | 23.82 | 23.81 |
| Sbarra principale | 25.15 | 25.17 | 25.35 |
| Sbarre di derivazione | 27.76 | 28.04 | 29.07 |
| Tubo di rame per sbarre di derivazione | 28.42 | 29.31 | 30.91 |
| Braccio di uscita inferiore | 24.95 | 24.85 | 26.33 |
Tabella 6: Valori di temperatura nei punti di monitoraggio di ciascuna fase della cabina ad anello sotto il modulo di soluzione sul campo di temperatura Icepak
Discussione
Questo documento è un'analisi di simulazione comparativa dell'aumento di temperatura dell'armadio ad anello basata su software di modellazione ingegneristica e software agli elementi finiti e la soluzione più adatta per l'effettiva situazione di aumento della temperatura viene analizzata da due moduli di soluzione del campo di temperatura agli elementi finiti. La gestione termica è anche descritta in Icoz23 come una componente critica ed essenziale per mantenere l'elevata efficienza e affidabilità dei componenti elettronici. L'importanza di condurre un'analisi comparativa è riassunta attingendo al lavoro di Steiner24: un'analisi comparativa è stata effettuata utilizzando COMSOL e ANSYS Mechanical. Pertanto, nel processo di produzione dell'RMU effettiva, la sua distribuzione della temperatura può essere analizzata da un software di simulazione finito, che può far risparmiare notevolmente manodopera e costi di produzione.
Il modello tridimensionale dell'RMU può essere creato utilizzando SolidWorks, come mostrato nella Figura 1A. Nel passaggio 1.1, la semplificazione del modello è enfatizzata come un passaggio cruciale nell'analisi degli elementi finiti23. Dati i numerosi componenti dell'RMU, che possono influire sull'accuratezza della soluzione e sui risultati computazionali, le parti non necessarie vengono rimosse dopo aver considerato solo i componenti conduttivi. Le parti di lavoro principali vengono mantenute, come illustrato nella Figura 1B.
Nella fase di pre-elaborazione del campo a correnti parassite, un aspetto critico riguarda l'aggiunta e la verifica dell'eccitazione del carico. Come descritto nel passaggio 2.1, è essenziale aggiungere con precisione la corrente di carico ai bracci in uscita superiore e inferiore, assicurandosi che l'intero circuito conduttivo formi un percorso per la corrente. Ciò facilita il calcolo della mappa delle nuvole di distribuzione del campo a correnti parassite. La difficoltà risiede nella necessità di controllare il circuito prima del calcolo della soluzione, poiché un percorso di corrente incompleto nel circuito può portare a una non convergenza durante la risoluzione o a situazioni in cui Maxwell non è in grado di calcolare. La successiva analisi del campo di temperatura si basa interamente sulle perdite ottenute dalla risoluzione di questo campo di correnti parassite, rendendo l'ispezione del percorso della corrente un passaggio necessario.
Nella soluzione del campo di temperatura, il passaggio 3 presenta due diversi moduli per la risoluzione del modello. Tuttavia, garantire condizioni iniziali identiche durante il processo di soluzione rappresenta una sfida. Poiché i moduli hanno enfasi diverse nelle loro soluzioni, è necessario approssimare o impostare temperature ambientali, coefficienti convettivi e altre condizioni di risoluzione identici per garantire la coerenza con l'unica fonte di calore fornita da Maxwell. Inoltre, il controllo coerente degli intervalli per la soluzione del campo di temperatura in entrambi i modelli consente un confronto diretto delle visualizzazioni di temperatura per gli stessi componenti nello stesso intervallo di temperatura, evidenziando le differenze nelle soluzioni e facilitando conclusioni intuitive.
Questo documento presenta un metodo innovativo per misurare e analizzare con precisione la distribuzione della temperatura delle apparecchiature elettriche durante il funzionamento. Affrontando la sfida di misurare con precisione le perdite ohmiche con i metodi tradizionali, questo studio utilizza il risolutore a correnti parassite Maxwell per calcoli accurati, che funge quindi da base per risolvere il campo di temperatura. Successivamente, il modulo di risoluzione del campo di temperatura viene utilizzato per visualizzare visivamente la distribuzione della temperatura di ciascun componente, migliorando significativamente l'efficienza e l'accuratezza dei test ingegneristici. La novità di questo articolo risiede non solo nel fornire un metodo efficiente per risolvere il campo di temperatura, ma anche nel dimostrare come utilizzare questi risultati per l'ottimizzazione strutturale e l'analisi della dissipazione del calore. Ciò offre nuovi mezzi tecnici per la progettazione e l'ottimizzazione delle apparecchiature elettriche. Tuttavia, ci sono state limitazioni nel processo di ricerca e nel lavoro futuro saranno condotte ulteriori ricerche nelle seguenti aree. In primo luogo, l'affinamento del modello e l'accoppiamento multi-campo. Il modello sarà ulteriormente perfezionato per riflettere accuratamente le caratteristiche termiche dell'armadio di rete ad anello durante il funzionamento effettivo. In secondo luogo, lo studio dei campi di temperatura transitori. Verranno condotte ricerche approfondite sulle variazioni del campo di temperatura dell'armadio della rete ad anello in condizioni di carico variabili per fornire un supporto aggiuntivo per il funzionamento stabile dell'armadio della rete ad anello in condizioni operative complesse.
Dichiarazioni
Gli autori non hanno interessi contrastanti.
Ringraziamenti
Gli autori ringraziano il signor Wu, MS Sun, il signor Wang, il signor Mu e il signor Li per il loro aiuto. Questo studio è stato supportato dalla China Postdoctoral Science Foundation (2022M721604) e dal Wenzhou Key Science and Technology Tackling Programmer (ZG2023015).
Materiali
| Nome | Azienda | Numero di catalogo | Commenti |
|---|---|---|---|
| Aria | / | / | Gas convenzionali |
| Alluminio | / | / | Lega Materiali |
| Rame | / | / | Lega |
| Materiali Icepak | ANSYS azienda | ANSYS 2021R1 | A CFD software di simulazione termica |
| PC hosting | / | Intel(R) Core(TM) i5-13500F CPU di 12a generazione | Apparecchiature |
| informatiche hostSolidWorks | Filiale di Dassault Systèmes | SolidWorks2021 | Uno strumento di disegno di software di ingegneria |
| ANSYS | termica allo stato stazionario | ANSYS 2021R1 | Uno strumento di soluzione di simulazione termica |
Riferimenti
- Xia, H., et al. Temperature rise test and analysis of high current switchgear in distribution system. J Engg. , 754-757 (2019).
- 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).
- 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).
- 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).
- 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).
- 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).
- 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).
- 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).
- 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).
- 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).
- 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).
- Ye Mao, Thermal simulation of high-current switch cabinet based on Icepak). Electr Ener Mgmt Technol. , 1-7 (2018).
- Zhang, S. Evaluation of thermal transient and overload capability of high-voltage bushings with ATP. IEEE Trans Power Delivery. 24 (3), 1295-1301 (2009).
- 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).
- Ş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).
- Ruibo, Y., et al. Research and application of temperature load of switchgear. J Physics: Conf Series. 2378 (2022), (2022).
- 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).
- 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).
- Sheikholeslami, M., Khalili, Z. Solar photovoltaic-thermal system with novel design of tube containing eco-friendly nanofluid. Renewable Ener. , 222(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).
- 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).
- 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).
- 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).
- Steiner, T. R. High temperature steady-state experiment for computational radiative heat transfer validation using COMSOL. and ANSYS. Results Engg. 13, 100354(2022).
Ristampe e permessi
