במחקר זה פותחה מערכת טעינה לרכבים חשמליים (EV) הניזונה ממערכת פוטו-וולטאיות (PV) דו-צדדית עצמאית וממערכת אחסון אנרגיה בסוללה (BESS), כפי שמוצג ב-איור 1, תוך שימוש בשיטת MPPT הנשלטת על ידי רשת עצבית מלאכותית (ANN). מעגל שווה-ערך בעל דיודה אחת משמש למידול המערך הפוטו-וולטאי הדו-צדדי, תוך התחשבות בטמפרטורה, התנגדות טורית, התנגדות מקבילית, קרינה סולארית וקרינה מהצד האחורי. הפלט הפוטו-וולטאי מסופק לתחנת הטעינה של רכב החשמל באמצעות ממיר DC–DC boost, בעוד שה-BESS תומכת בוויסות מתח אפיק ה-DC. אלגוריתם אופטימיזציה של ציפור מזכיר (SBOA) משמש לאופטימיזציה של הגב리를 של בקר ה-PI ושל פרמטרי המסנן. מידול רכיבי המערכת והמתודולוגיה המוצעת מתוארים בסעיפי המשנה הבאים.
מידול של רכיבים
1. מערכת PV דו-צדדית
המודולים הפוטוולטאיים (PV) מסייעים בהשגת פלט המתח והזרם הנדרשים. במחקר הנוכחי, שני מודולי PV חד-פניתיים שימשו לסימולציה של מאפייני ייצור האנרגיה של מערכים דו-פניתיים תחת תנאי אלבדו שונים של פני השטח9. הרווח הדו-פניתי של הקרינה BGg נקבע לפי משוואה (1)
(1)
כאן, הקרינה הקדמית היא Gf, והקרינה האחורית היא Gr. זרם ה-PV הביפסיאלי המקסימלי הוא Imp, והמתח הוא Vmp. לבסוף, התפוקה המיוצרת על ידי המערכת הביפסיאלית עם מספר המודולים הנבחר בטור Ns ובמקביל Np ניתנת על ידי משוואה (2). הבקרת ה-PV מוצגת ב-איור 2
(2)
2. מערכת אגירת סוללות (BSS)
סוללות הן אמצעי קריטי לאגירה של אנרגיה חשמלית18. עבודה זו משתמשת בסוללות Li-ion כמערכת לאגירת אנרגיה. מבין אפשרויות הסוללות הניתנות לטעינה, עבודה זו מתמקדת בסוללות Li-ion, הנחשבות לאופטימליות עבור רכבים חשמליים (EVs) בשל יעילותן וצפיפות ההספק שלהן. משוואות (3) ו-(4) מתארות את סוללת ה-Li-ion עם מגבלות.
(3)
המפרטים של ה-PV, סוללת האגירה וההמרה מסוכמים בטבלה 2, בעוד שטבלה 3 מציגה את הקצאת זרימת ההספק בין מערך ה-PV, יחידת אגירת האנרגיה והעומסים המחוברים.
(4)
3. ממיר DC–AC
ממירי המתח בעלי נקודת ניטרל מקובעת (NPC) מאומצים באופן נרחב במערכות המרה של הספק עבור אנרגיה מתחדשת. בהשוואה למבנה הדו-רמתי המסורתי, תצורת ה-NPC מציעה עיוות הרמוני כולל (THD) מופחת, הפרעה אלקטרומגנטית נמוכה יותר ומאמץ מתח נמוך יותר על התקני המוליכים למחצה. חיבור ה-DC של טופולוגיה זו מחולק לשלוש רמות פוטנציאל +Vdc/2, 0, ו- -Vdc/2, באמצעות קבלים זוגיים ודיודות קיבוע המהווים את צומת הניטרל. על ידי בקרה מתאימה של המיתוג, הממיר מייצר שלושה מצבי מתח בדידים בטרמינלי הפלט, ובכך משפר את איכות גל המתח ואת יעילות המרת ההספק. נקודת הניטרל של פס ה-DC בממיר NPC ממלאת תפקיד מרכזי באיזון מתחי הקבלים ובהבטחת פעולה תקינה של הממיר. סכימת הבקרה של הממיר מוצגת ב-איור 3.
4. ממיר מתח an-DC למתח DC (מיר מעלה)
ממירי DC-DC מעלים את מתח הכניסה (DC) למתח מוצא גבוה יותר על ידי בקרת מחזור העבודה (duty cycle) של רכיב מיתוג בתדר גבוה, תוך שמירה על נצילות. לממיר שני מצבי פעולה. במהלך מצב ה-ON, הדיודה נמצאת בממתח אחורי (reverse-biased), ומקור הכניסה מספק אנרגיה לסליל, האוגר אותה כשדה מגנטי, בעוד קבל המוצא מספק אנרגיה לעומס. במהלך מצב ה-OFF, הדיודה עוברת לממתח קדמי (forward-biased). בשלב זה, האנרגיה האגורה בסליל משתחררת דרך הדיודה אל הקבל המחובר למוצא ולעומס, ובשילוב עם מקור הכניסה היא מייצרת מתח מוצא הגבוה ממתח הכניסה. על ידי שינוי מחזור העבודה של ה-PWM, מתח המוצא של הממיר משתנה בהתאם לקשר
, כאשר D הוא מחזור העבודה (duty cycle). הוא משמש במערכות סולאריות, תאי דלק, רכבים חשמליים (EV) ומערכות אגירה.
סכמת בקרה של רשת נוירונים מלאכותית (ANN) עבור MPPT
בעבודה זו, נבחר MPPT מבוס ANN כדי לאסוף את ההספק המקסימלי המופק ממערכת הסולארית17מערכי נתונים של קרינה וטמפרטורה ששימשו לאימון ה-ANN הופקו ב-MATLAB/Simulink תחת תנאי תפעול סביבתיים משתנים, התואמים לרמות קרינה סולארית של 80–10 W/m²2 וסיכומי טמפרטורה של 20–25 °C, עם מחזור עבודה (D). רשת העצבים המלאכותית (ANN) מאומנת כדי למזער את השגיאה הריבועית הממוצעת (MSE) בין הפלטים שהתקבלו לפלטים הנדרשים (Op,
כדי להפיק את ההספק המקסימלי ממערכת הפוטו-וולטאיקה (PV) הדו-צדדית. כאן, הקרינה הסולארית והטמפרטורה משמשות כקלטים לרשת העצבית מלאכותית (ANN), ומחזור העבודה (duty cycle) הוא הפלט לשליטה בממיר המעלה (boost converter). מבנה מודל ה-ANN שפותח עבור MPPT, כפי שמוצג ב- איור 4, נלקח בחשבון בעבודה. טבלה משלימה 1 מציג את היתרונות והחסרונות של השיטה המוצעת בהשוואה לשיטות סטנדרטיות אחרות.
רשת עצבית מלאכותית (ANN) מורכבת משלוש שכבות עיקריות: שכבת פלט (OPL), שכבת קלט (IPL) ושכבה נסתרת (HIL). מבנה זה מסייע בהעברת הנתונים בין ה-IPL ל-HIL. לאחר מכן, אלו מופקים על ידי המשקולות בקישורים המחברים בין ה-IPL ל-HIL. בהקשר זה, חישובים מבוצעים עם הטיה (bias) מסוימת המוחלת על משתנה ה-HL, והתוצאות המתקבלות נצברות במשתנה ה-OPL. כאן נבחרה ANN מסוג LMBP [17]. משקולות החיבור מכווננות במהלך האימון על ידי מדידת השגיאה כדי להגיע לפלט המטרה. כאן נעשה שימוש ב-LMBP לאימון רשתות ANN כאשר ה-MSE משמש כפונקציית הביצועים. אלגוריתם ה-LMBP משתמש בנגזרות הגזורות כדי לעדכן את המשקולות, המציגות תכונות של למידה יעילה והסגרה (convergence) מואצת
לכל נוירון ברשת פרספטרון רב-שכבית יש פונקציית סכימה ופונקציית הפעלה. עם זאת, קיימים משקולות מספריות מסוימות (wpk) המחברות בין נוירונים אלו לאורך השכבות. כאשר קלטים מוכפלים במשקולות כפי שמפורט במשוואה (5), נלקחת בחשבון פונקציית סיגמיד לא-ליניארית עם מזעור של MSE כפי שמוצג במשוואה (6).
(5)
(6)
SBOA ביצעה אופטימיזציה למסנן ולפרמטרי הבקרה של המערכת המוצעת
בימים אלה, אלגוריתמים מטא-היוריסטיים ממלאים תפקיד מרכזי בפתרון בעיות הנדסיות. איור 5 מציג סיווג של אלגוריתמים. אלגוריתם SBOA23 מחקה את התנהגותו של ה-SB כדי לשרוד בבית גידולו הטבעי. שלבי החקירה (exploration) והניצול (exploitation) של התנהגות הציד של ה-Secretary Bird ממודלים ב-SBOA. שלב החקירה של האלגוריתם מחקה את התנהגות ה-SB בתפיסת נחשים, בעוד ששלב הניצול שלו מחקה את התנהגותם בהתחמקות מטורפים כגון נשרים. תהליך הפתרון מתחיל ביצירת ערכים אקראיים בתוך הגבולות המתאימים עבור כל SB באוכלוסייה, ובחישוב ערך פונקציית המטרה על ידי הרצת מודל Simulink, תוך התייחסות לערכו של כל SB כפרמטר תכנוני. בהתבס על ערך המטרה, כל SB משנה את מיקומו כדי לייצג פתרון טוב יותר. המידול המתמטי של SBOA נדון להלן:
שלב ההכנה הראשוני
ראשית, יש לזהות את הפתרונות הראשוניים ששימשו להתחלת החיפוש אחר מינימום טיפוסי של פונקציית המטרה F(Y). במקרה זה, האוכלוסייה האקראית הראשונית Y = [Y1, Y2, ......YN] של SB‘s נוצרה באמצעות N פתרונות ראשוניים. משוואה (7) מאתחלת את האוכלוסייה, כאשר Yi מייצגת פתרון.
(7)
כאשר הגבולות התחתון והעליון של משתני הבחירה מסומנים על ידי lb ו-ub. בטווח [0, 1], r הוא מספר אקראי. N הוא ממד הבעיה. יתר על כן, ערך הכושר של הפתרון Fi = F(Yi) משמש לכימות האיכות של Yi.
אסטרטגיית הציד של SB’s
איתור, אכילה ותקיפת הטרף הם שלושת השלבים העיקריים של הציד של ה-SB. תהליך הציד כולו חולק לשלושה פרקי זמן שווים, t < 1/3T, 1/3T < t < 2/3T, ו-2/3T < t < T. כאן, t הוא האיטרציה הנוכחית, ו-T הוא מספר האיטרציות המקסימלי. פרקים אלו תואמים לשלושת שלבי הטורפה של ה-SB: חיפוש, צריכה ותקיפת טרף. חלוקות אלו מבוסות על הסטטיסטיקה הביולוגית שלבי הציד של הציפור ועל משכי הזמן של כל שלב. כתוצאה מכך, להלן האופן שבו מומדל כל שלב של ה-SBOA:
חיפוש אחר טרף
על ה-SB לצוד טרף ממרחק בטוח במהלך שלב זה. כדי לאסוף נתונים מספיקים על פני כל אזור החיפוש, השלב הראשון באלגוריתמי אופטימיזציה דורש חקירה (exploration) חזקה יותר. באמצעות שימוש במיקומם של שני ה-SB-ים האחרים כמדריך, ה-SB יכול לבחון אזורים אפשריים נוספים. כדי להגביר עוד יותר את הגיוון של האלגוריתם, מוצגים תהליכי מוטציה דיפרנציאליים. כאשר t קטן מ-1/3 מ-T, משתמשים במשוואה (8) כדי לעדכן את המיקום של כל פרט. הוא יחפש, יצרוך ויתקוף טרף. כתוצאה מכך, כך ממודל כל שלב ב-SBOA:
(8)
כאשר, Yi(t) הוא המיקום הנוכחי של ה-SB ה-ith, Ynewi(t) מייצג את הפתרון המעודכן, ושני פרטים Yi(t), Yr2(t) נבחרו באופן אקראי מהאוכלוסייה הקיימת. וקטור אקראי, עבור אסטרטגיית בחירת הפרמטרים, R1 בעל גודל 1 × N נבחר אקראית מתוך [0, 1] למען איזון יעיל בין שלבי הניצול (exploitation) והסריקה (exploration). הפתרון המעודכן נשקל רק אם הוא מוביל לשיפור בערך פונקציית המטרה בהשוואה לפתרון שהתקבל קודם לכן.
צריכת טרף
ה-SB משתמש בזמינות ובגמישות ותמרן להמראה סביב הנחש לאחר זיהוי טרף פוטנציאלי. סבלנותו של הטרף תנוצל להחליש את הגנאות שלו באמצעות תצפית ופיתוי של היריבים בעודם מקיפים אותו. SB-ים אחרים מתאימים את תנוחתם בקרבת הטרף. הצלחת הציד תעלה באופן משמעותי בגישה זו. משוואה (9) ממחישה כיצד תנועת בראוניאנית (B) משמשת לסימולציה של תנועתו האקראית של ה-SB כאשר 1/3 T < t < 2/3 T.
(9)
כאשר B מומש ב-MATLAB באמצעות התפלגות נורמלית סטנדרטית על ידי הווקטור r(1,M) שנוצר באופן אקראי. לאחר מכן, ה-SB מעדכנים את מיקומם באמצעות משוואה (10).
(10)
כאשר, Ybest(t) הוא הפתרון הטוב ביותר, ומקדם scaling האקספוננציאלי exp((t/T)4) נבחר כדי להגביר את יכולת הניצול (exploitation) בכל איטרציה.
תקיפת טרף
לאחר צריכה רציפה, הקורבן יסיים את מאגרי האנרגיה שלו. עצי השכרונים צריכים להשיק את ההתקפה כעת. כאן, שיטות תקיפה שונות, כגון צעדים בלתי פוסקים / קפיצות ארוכות ספורדיות לאורך פרק זמן קצר, סומללו באמצעות טכניקת טיסת לוי (Lévy flight). המאפיינים מתוארים על ידי משוואות (1) ו-(12). האפשרויות המועמדות קרובות לתשובה האופטימלית ברגע זה מכיוון ש-SB מתקרבים במהירות לטרף. תכנית זו תיושם כאשר t > 2/3T.
(11)
(12)
היכן, LF מייצג את אסטרטגיית טיסות לבי (Levy flight) ו-
, s = 0.01, ו-φ = 1.5 ערכים קבועים. μ ו- נראה שהקלט שסיפקת אינו מכיל טקסט באנגלית לתרגום. אנא ספק את הטקסט המדעי שברצונך לתרגם לעברית. הם ערכים שרירותיים הנמצאים בטווח [0, 1]. כאשר, τ מסמן את גמא ו- η = 0.5, נבחר מקדם השקלול 0.5, כפי שהוצע בניסוח המקורי של SBOA, כדי להסדיר את גודל הצעד ולשמור על איזון בין קפיצות חקירה גדולות לבין זיקוק מקומי. טבלה 4 מציג את הערכים שבהם נעשה שימוש עבור ה-ANN ואלגוריתם ה-SB במערכת שפותחה.
אסטרטגיית בריחה עבור SB
בעת ציד טרף אחר בטבע, ציפורי מזכיר (SBs) מסתכנות בכך שיהיו הן הטרף. עיטים, נצים, שועלים ובניקלבים הם היריבים העיקריים עימם עליהן להתמודד. הן חייבות להשתמש במגוון טכניקות התחמקות כדי להגן על עצמן או על טרפן כאשר הן מזהות סכנה. אלגוריתם זה מדמה את שיטות הבריחה על ידי מידול מצבי ריצה (D2) והסוואה (D1). כאשר הן ניצבות מול יריבים, ציפורי המזכיר נטמעות תחילה בסביבתן כדי להישאר בטוחה. כאן, ציפורי המזכיר שואפות להימנע מאופטימום מקומי באלגוריתמים על ידי עדכון מיקומן סביב הטרף (הפרט הטוב ביותר). הן ישתמשו בטקטיקות של תעופה או ריצה מהירה כדי להישאר בטוחות אם לא יצליחו להתחמק מהאוייב. לצורך התייחסות, פרט אקראי Yrand נבחר כמוביל כדי למנוע הגבלה לאופטימום מקומי. ציפורי המזכיר משתמשות במשוואות (13) ו-(14) כדי לעדכן את מיקומן. תרשים הזרימה של SBOA מוצג ב-איור 6.
(13)
(14)
ייצוג של משתני תכנון
משתני הבעיה במחקר זה כוללים את פרמטרי התכנון, כגון הגב리를 ה-PI של ממירים DC-DC, כמו גם את ההתנגדויות וההשראויות של המסננים והמירים. משוואה (15) מציגה את ייצוג משתני הבקרה. גבולות פרמטרי הבקרה מפורטים ב-טבלה משלימה 2.
(15)
פונקציית כשירות (FF)
במחקר זה נבחרה המינימיזציה של THD כפונקציית המטרה (Obj). הערך המקסימלי של FF הוא הערך המינימלי של Obj, כמפורט במשוואות (16) ו-(17).
(16)
(17)