In questa ricerca, viene sviluppato un sistema autonomo di ricarica per veicoli elettrici (EV) alimentato da un impianto fotovoltaico bifacciale e da un sistema di accumulo energetico a batteria (BESS), come mostrato nella Figura 1, utilizzando un metodo MPPT controllato da una rete neurale artificiale (ANN). Per modellare l'array fotovoltaico bifacciale si utilizza un circuito equivalente a singolo diodo, tenendo conto della temperatura, della resistenza in serie, della resistenza in parallelo, dell'irradianza solare e dell'irradianza sul lato posteriore. L'uscita del fotovoltaico viene fornita alla stazione di ricarica del veicolo elettrico attraverso un convertitore CC–CC rialzatore, mentre il BESS supporta la regolazione della tensione del bus in corrente continua. L'algoritmo di ottimizzazione ispirato al comportamento del serpentario (Secretary Bird Optimization Algorithm, SBOA) viene utilizzato per ottimizzare i guadagni del controllore PI e i parametri del filtro. La modellizzazione dei componenti del sistema e la metodologia proposta sono descritte nei paragrafi seguenti.
Modellizzazione dei Componenti
1. Sistema fotovoltaico bifacciale
I moduli fotovoltaici contribuiscono a raggiungere l'output di corrente e tensione richiesto. Nello studio attuale, sono stati utilizzati due moduli fotovoltaici monofacciali per emulare le caratteristiche di generazione energetica di array bifacciali in diverse condizioni di albedo superficiale9. Il guadagno bifacciale dell'irradianza BGg è determinato dall'Equazione (1)
(1)
In questo caso, l'irradiazione frontale è Gf, e l'irradianza posteriore è Gr. La corrente massima del fotovoltaico bifacciale è Imp, e la tensione è Vmp. Infine, l'output prodotto dal sistema bifacciale con il numero scelto di moduli in serie Ns e in parallelo Np è dato dall'Equazione (2). Il controllo del fotovoltaico è mostrato nella Figura 2
(2)
2. Sistema di accumulo dell'energia (BSS)
Le batterie sono un mezzo fondamentale per immagazzinare energia elettrica18. Questo studio utilizza batterie agli ioni di litio come sistema di accumulo energetico. Tra le opzioni disponibili di batterie ricaricabili, questo lavoro si concentra sulle batterie agli ioni di litio, considerate ottimali per i veicoli elettrici grazie alla loro efficienza e densità di potenza. Le equazioni (3) e (4) descrivono la batteria agli ioni di litio con vincoli.
(3)
Le specifiche del pannello fotovoltaico, della batteria di accumulo e del convertitore sono riassunte nella Tabella 2, mentre la Tabella 3 presenta la ripartizione del flusso di potenza tra il campo fotovoltaico, l'unità di accumulo dell'energia e i carichi collegati.
(4)
3. Convertitore CC – CA
L'inverter a punto neutro clamped (NPC) è ampiamente utilizzato nei moderni sistemi di conversione dell'energia per le fonti rinnovabili. Rispetto alla struttura convenzionale a due livelli, la configurazione NPC offre una ridotta distorsione armonica totale (THD), un minore disturbo elettromagnetico e una minore sollecitazione di tensione sui dispositivi semiconduttori. Il collegamento in corrente continua (DC link) di questa topologia è suddiviso in tre livelli di potenziale: +Vdc/2, 0 e -Vdc/2, mediante coppie di condensatori e diodi di clamp che formano il nodo neutro. Controllando opportunamente la commutazione, l'inverter genera tre stati di tensione discreti ai terminali di uscita, migliorando così la qualità dell'onda e l'efficienza di conversione dell'energia. Il punto neutro del bus in corrente continua in un inverter NPC svolge un ruolo fondamentale nel bilanciare le tensioni dei condensatori e garantire un corretto funzionamento dell'inverter. Lo schema di controllo dell'inverter è mostrato in Figura 3.
4. Convertitore CC-CC rialzatore
I convertitori CC-CC innalzano la tensione di ingresso in una tensione di uscita più elevata controllando il ciclo di lavoro di un dispositivo di commutazione ad alta frequenza, mantenendo nel contempo l'efficienza. Dispone di due modalità operative. Durante lo stato ON, il diodo è polarizzato inversamente e la sorgente di ingresso fornisce energia all'induttore, che la immagazzina sotto forma di campo magnetico, mentre il condensatore di uscita fornisce energia al carico. Durante lo stato OFF, il diodo diventa polarizzato direttamente. In questa fase, l'energia immagazzinata nell'induttore viene rilasciata attraverso il diodo al condensatore collegato all'uscita e al carico, sommandosi alla sorgente di ingresso per produrre una tensione di uscita maggiore rispetto a quella di ingresso. Variando il ciclo di lavoro del PWM, la tensione di uscita del convertitore varia secondo la relazione
, dove D è il ciclo di lavoro. È utilizzato nei sistemi solari, a celle a combustibile, nei veicoli elettrici e nei sistemi di accumulo.
Schema di controllo ANN per MPPT
In questo studio, viene selezionato un MPPT basato su reti neurali artificiali (ANN) per raccogliere la potenza massima generata dal sistema solare17. I dataset di irraggiamento e temperatura utilizzati per l'addestramento della rete ANN sono stati generati in MATLAB/Simulink in condizioni operative ambientali variabili, corrispondenti a livelli di irraggiamento solare compresi tra 800 e 1000 W/m2 e intervalli di temperatura tra 20 e 25 °C, con un duty cycle (D). La rete ANN viene addestrata per minimizzare l'errore quadratico medio (MSE) tra le uscite ottenute e quelle richieste (Op,
) al fine di estrarre la potenza massima dal sistema fotovoltaico bifacciale. In questo caso, l'irraggiamento solare e la temperatura vengono utilizzati come ingressi per la rete ANN, mentre il duty cycle rappresenta l'uscita per controllare il convertitore boost. La struttura del modello ANN sviluppato per il MPPT, mostrata in Figura 4, è quella considerata nel presente lavoro. Tabella supplementare 1 presenta i vantaggi e gli svantaggi del metodo proposto rispetto ad altri metodi standard.
Una rete neurale artificiale (ANN) è composta da tre strati principali: uno strato di uscita (OPL), uno strato di ingresso (IPL) e uno strato nascosto (HIL). Questo consente il trasferimento dei dati tra IPL e HIL. Successivamente, i dati vengono prodotti dai pesi nei collegamenti che uniscono IPL e HIL. In questo contesto, i calcoli vengono eseguiti applicando un certo bias alla variabile HL, e i risultati ottenuti vengono accumulati nella variabile OPL. In questo caso, viene selezionata un'ANN di tipo LMBP [17]. I pesi dei collegamenti vengono regolati durante l'addestramento misurando l'errore al fine di raggiungere l'uscita desiderata. Qui, LMBP viene utilizzato per addestrare le ANN con l'MSE come funzione di prestazione. L'algoritmo LMBP utilizza le derivate ottenute per aggiornare i pesi, che presentano le caratteristiche di un apprendimento efficace e di una convergenza accelerata.
Ogni neurone in una rete di perceptron multistrato ha una funzione di somma e una funzione di attivazione. Tuttavia, esistono determinati pesi numerici (wpk) che collegano questi neuroni tra i diversi livelli. Quando gli ingressi vengono moltiplicati per i pesi come specificato nell'Equazione (5). Viene considerata la funzione sigmoidea non lineare con la minimizzazione dell'errore quadratico medio data nell'Equazione (6).
(5)
(6)
SBOA ha ottimizzato i parametri del filtro e di controllo del sistema proposto
Oggi, gli algoritmi metaeuristici svolgono un ruolo fondamentale nella risoluzione di problemi ingegneristici. Figura 5 mostra una classificazione degli algoritmi. L'algoritmo SBOA23 imita il comportamento del SB per sopravvivere nel suo habitat naturale. Le fasi di esplorazione e sfruttamento del comportamento di caccia del Gruccione del Serpente sono modellate nell'SBOA. La fase di esplorazione dell'algoritmo imita il comportamento del SB nell'acchiappare serpenti, mentre la fase di sfruttamento imita il suo comportamento nell'evitare predatori come aquile. Il processo di soluzione inizia generando valori casuali entro i rispettivi limiti per ogni SB nella popolazione e calcolando il valore della funzione obiettivo eseguendo il modello Simulink, trattando ciascun valore del SB come un parametro di progetto. Sulla base del valore obiettivo, ogni SB modifica la propria posizione per rappresentare una soluzione migliore. La modellizzazione matematica dell'SBOA è discussa di seguito:
Fase iniziale di preparazione
Innanzitutto, è necessario identificare le soluzioni iniziali utilizzate per avviare la ricerca di una tipica minimizzazione della funzione obiettivo F(Y). In questo caso, la popolazione iniziale casuale Y = [Y1, Y2, ......YN] di SB viene formata tramite N soluzioni iniziali. L'equazione (7) inizializza la popolazione e Yi rappresenta una soluzione.
(7)
Dove i limiti inferiore e superiore delle variabili di scelta sono indicati rispettivamente con lb e ub. In [0, 1], r è un numero casuale. N è la dimensione del problema. Inoltre, il valore di adattamento della soluzione Fi = F(Yi) viene utilizzato per quantificare la qualità di Yi.
Strategia di caccia degli SB
Individuare, mangiare e attaccare la preda sono le tre fasi principali della caccia degli SB. L'intero processo di caccia è stato suddiviso in tre intervalli di tempo uguali, t < 1/3T, 1/3T < t < 2/3T e 2/3T < t < T. Qui, t rappresenta l'iterazione corrente, mentre T indica il numero massimo di iterazioni. Questi intervalli corrispondono alle tre fasi della predazione da parte degli SB: ricerca, consumo e attacco alla preda. Tale suddivisione si basa sulle statistiche biologiche relative alle fasi di caccia dell'uccello e sulla durata temporale di ciascuna fase. Di conseguenza, ogni fase dell'SBOA viene modellata nel modo seguente:
Ricerca della preda
L'SB deve cacciare la preda da una distanza di sicurezza durante questa fase. Per raccogliere dati sufficienti su tutta la regione di ricerca, il primo passo negli algoritmi di ottimizzazione richiede un'esplorazione più intensa. Utilizzando le posizioni degli altri due SB come guida, l'SB può esplorare ulteriori regioni potenziali. Per aumentare ulteriormente la diversità dell'algoritmo, vengono quindi introdotti processi di mutazione differenziale. Quando t è inferiore a 1/3 di T, si utilizza l'Equazione (8) per aggiornare la posizione di ciascun individuo. Esso cercherà, consumerà e attaccherà la preda. Di conseguenza, ogni fase dell'SBOA è modellata nel modo seguente:
(8)
Dove, Yi(t) è la posizione attuale dell'iesimo SB, Ynewi(t) fornisce la soluzione aggiornata, due individui Yi(t), Yr2(t) sono stati scelti casualmente dalla popolazione esistente. Un vettore casuale, per la strategia di selezione dei parametri, R1 di dimensione 1 × N è scelto casualmente nell'intervallo [0, 1] per garantire un efficace equilibrio tra le fasi di sfruttamento ed esplorazione. La soluzione aggiornata viene considerata solo se porta a un valore obiettivo migliore rispetto alla soluzione ottenuta in precedenza.
Consumo della preda
SB utilizza un'agilità rapida e si muove per indugiare intorno al serpente dopo aver individuato una possibile preda. La pazienza della preda verrà sfruttata per abbassarne le difese, osservando e attirando gli avversari mentre si muovono in cerchio. Altri SB aggiustano la propria postura vicino alla preda. Il successo della caccia aumenterà in modo significativo con questo approccio. L'equazione (9) illustra come il moto browniano (B) venga utilizzato per simulare il movimento casuale di SB quando 1/3 T < t < 2/3 T.
(9)
Dove B è implementato in MATLAB utilizzando la distribuzione normale standard tramite il vettore generato casualmente r(1,M). I SB utilizzano quindi l'Equazione (10) per aggiornare le loro posizioni.
(10)
Dove, Ybest(t) è la soluzione migliore e il fattore di scala esponenziale exp((t/T)4) è scelto per aumentare la capacità di sfruttamento ad ogni iterazione.
Attacco alla preda
Dopo aver consumato energia in modo continuo, la vittima esaurirà le proprie risorse. I serpenti del Capo dovrebbero lanciare l'attacco ora. Qui, diversi metodi di attacco, come passi incessanti o salti lunghi sporadici in un breve periodo di tempo, vengono simulati utilizzando la tecnica del volo di Lévy. Le caratteristiche sono descritte dalle Equazioni (11) e (12). In questo momento, le opzioni candidate si trovano vicino alla soluzione ottimale poiché i SB si avvicinano rapidamente alla preda. Questo piano verrà attuato quando t > 2/3T.
(11)
(12)
Dove, LF rappresenta la strategia di volo di Levy e
, con s = 0,01 e φ = 1,5 valori fissi. μ e υ sono valori arbitrari compresi nell'intervallo [0, 1]. Inoltre, τ indica la funzione gamma e η = 0,5; il fattore di peso 0,5 è stato adottato, come proposto nella formulazione originale dell'SBOA, per regolare l'ampiezza del passo e mantenere un equilibrio tra grandi salti esplorativi e affinamenti locali. Tabella 4 presenta i valori utilizzati per la rete neurale (ANN) e l'algoritmo SB nel sistema sviluppato.
Strategia di fuga per SB
Nella ricerca di altre prede in natura, i SB rischiano a loro volta di essere cacciati. Aquile, falchi, volpi e sciacalli sono i principali avversari con cui devono confrontarsi. Devono utilizzare una varietà di tecniche di evasione per proteggere se stessi o le loro prede quando percepiscono un pericolo. Questo algoritmo simula i metodi di fuga modellando le modalità di corsa (D2) e il mimetismo (D1). Quando affrontano avversari, i sergentari inizialmente si confondono con l'ambiente circostante per restare al sicuro. In questo contesto, i sergentari cercano di evitare ottimi locali nell'algoritmo aggiornando le proprie posizioni intorno alla preda (l'individuo migliore). Se non riescono a eludere il nemico, ricorreranno a tattiche di volo o di corsa rapida per mettersi in salvo. A tal fine, un individuo casuale Yrand viene scelto come leader per evitare di rimanere intrappolato in un ottimo locale. Le equazioni (13) e (14) vengono utilizzate dai sergentari per aggiornare le proprie posizioni. Il diagramma di flusso dello SBOA è mostrato in Figura 6.
(13)
(14)
Rappresentazione delle variabili di progetto
Le variabili del problema in questo studio comprendono i parametri di progettazione, come i guadagni del controllore PI dei convertitori CC-CC, nonché le resistenze e le induttanze dei filtri e dei convertitori. L'equazione (15) fornisce la rappresentazione delle variabili di controllo. I limiti dei parametri di controllo sono elencati nella Tabella Supplementare 2.
(15)
Funzione di adattamento (FF)
La minimizzazione della THD è stata scelta come funzione obiettivo (Obj) in questo studio. Il valore massimo di FF corrisponde al valore minimo di Obj, indicato nelle equazioni (16) e (17).
(16)
(17)