۱٫ معرفی
وسایل نقلیه الکتریکی (EVs) به عنوان یک راه حل اساسی در کاهش آلودگی محیط زیست و کاهش CO پدید آمده اند.
۲ انتشار گازهای گلخانه ای در سطح جهان، جایگزین وسایل نقلیه معمولی می شود [
۱,
۲,
۳]. با این حال، محدوده محدود خودروهای برقی عمدتاً از محدودیتهای فعلی در فناوری باتری ناشی میشود. پرداختن به این محدودیت مستلزم توسعه وسایل نقلیه کارآمدتر، به حداکثر رساندن ظرفیت باتری موجود برای افزایش برد آنها است. نقش محوری ماشینهای الکتریکی در EV، که انواع مختلفی از جمله ماشینهای DC، ماشینهای AC و ماشینهای رلوکتانس را در بر میگیرد، بهطور قابل توجهی قابل توجه است. در این میان، موتور سنکرون مغناطیس دائم (PMSM) به دلیل راندمان و عملکرد بالا همراه با ویژگی هایی مانند سادگی در ساخت، در دسترس بودن بالا، کارایی و نگهداری راحت، متمایز می شود و آن را در حوزه های مختلف از جمله کاربردهای فضایی، سیستم های اتوماسیون، همه کاره می کند. دفاع و به ویژه صنعت خودروسازی [
۴,
۵,
۶].
در خودروهای برقی معاصر، برجستگی ماشین های الکتریکی آهنربای دائم به طور قابل توجهی افزایش یافته است. این ماشین های الکتریکی نه تنها به عنوان موتورهای کششی بلکه به عنوان منابع انرژی در طول ترمز احیا کننده نیز خدمت می کنند، برای سرعت و گشتاور از طریق اینورتر کنترل می شوند که ولتاژ DC باتری را به ولتاژ AC تبدیل می کند تا ماشین الکتریکی را به حرکت درآورد. علاوه بر این، باتری علاوه بر سیستم درایو EV، برق تجهیزات کمکی مانند روشنایی، تهویه مطبوع و سیستمهای امنیتی را تامین میکند. در طول انتقال انرژی ترمز احیا کننده به باتری، فعال کردن یک مقاومت ترمز برای اطمینان از تداوم ترمز احیا کننده ضروری است. شارژ باتری از شبکه سه فاز شامل استفاده از یکسو کننده و مدار مبدل تقویت کننده است. [
۷,
۸].
انواع مختلفی از ماشین های الکتریکی در خودروهای الکتریکی کاربرد دارند [
۹,
۱۰,
۱۱,
۱۲]. PMSM ها به دلیل راندمان و قابلیت اطمینان بالا مورد علاقه بسیاری از سازندگان هستند که نیاز به تحریک خارجی را از بین می برند و تلفات مس روتور را به حداقل می رسانند و در نتیجه طراحی فشرده ای دارند. [
۱۳,
۱۴,
۱۵,
۱۶,
۱۷,
۱۸]. با این حال، اتکا به آهنرباهای خاکی کمیاب چالش هایی را در رابطه با ناامنی عرضه، نوسانات قیمت و خطر مغناطیس زدایی حرارتی ایجاد می کند. از سوی دیگر، ماشین القایی (IM) گشتاور و عملکرد قابل مقایسه ای را هنگام استفاده در PMSM ها، به ویژه در سرعت های بالا ارائه می کند. [
۱۹,
۲۰]. از مزایای IM می توان به عدم وجود آهنربا، ساختار مستحکم و مقرون به صرفه بودن اشاره کرد [
۲۱,
۲۲,
۲۳,
۲۴]. با هدف دستیابی به چگالی توان بالا، PMSM را به عنوان تمرکز این مطالعه انتخاب کردیم [
۲۵].
تأثیر متقابل اجزای الکتریکی و مکانیکی در خودروهای الکتریکی به شدت بر مصرف انرژی و عملکرد خودرو تأثیر میگذارد. نیروهای مختلف، از جمله مقاومت چرخ، مقاومت در برابر تپه، مقاومت در برابر شتاب، مقاومت آیرودینامیکی و شناوری، به طور قابل توجهی بر مصرف انرژی و عملکرد خودرو تأثیر میگذارند. [
۲۶,
۲۷].
شکل ۱ به وضوح تغییرات ضرایب مقاومت آیرودینامیکی خودرو را با توجه به طراحی خودرو نشان می دهد.
از لحاظ تاریخی، آیرودینامیک خودرو در طراحی خودرو مورد تاکید قرار نگرفته است. با این حال، افزایش مصرف انرژی و پیگیری عملکرد بالاتر، توجه سازندگان را به این جنبه معطوف کرده است. روش های ریاضی، به ویژه دینامیک سیالات محاسباتی (CFD)، نقش اساسی در طراحی بدنه خودرو، به ویژه در بهینه سازی آیرودینامیک دارند. [
۲۹].
کشش آیرودینامیکی به عنوان یک عامل تعیین کننده عملکرد خودرو است [
۲۶]. حرکت خودرو که عمدتاً از اختلاف فشار در جلو و پشت خودرو ناشی میشود، منجر به افزایش فشار استاتیک در جلو و تشکیل یک گرداب به دلیل جریان نامنظم هوا در عقب میشود. این جریان هوای نامتقارن باعث ایجاد نرخ جریان نسبی در هنگام حرکت خودرو می شود و در نتیجه نیروی مقاومت را به دلیل تغییر در ناحیه جلویی خودرو افزایش می دهد. [
۳۰].
عملکرد خودرو را می توان با در نظر گرفتن نیروی رو به پایین، افزایش قدرت مانور با کاهش نیروی مقاومت و افزایش نیروی رو به پایین افزایش داد. [
۳۱]. زاویه شیب در عقب خودرو به میزان قابل توجهی به مقاومت فشار کمک می کند و بخش قابل توجهی از مقاومت کلی خودرو را نشان می دهد. [
۳۲]. مطالعاتی بر روی ساختار گرداب و تغییر در زاویه عقب خودرو انجام شده است. شناخته شده ترین این مطالعات مدل بدن احمد است [
۳۳]. تلاشها برای بهینهسازی کارایی و ایمنی خودرو شامل به حداقل رساندن گردابهایی است که در نواحی محلی خودرو اتفاق میافتد. سازندگان خودرو اغلب شاسی را صاف می کنند تا این گرداب ها را به حداقل برسانند [
۲۶].
عوامل موثر بر مصرف انرژی خودرو از پارامترهای طراحی خودرو، ویژگیهای اجزا، استفاده از دستگاه کمکی و شرایط خارجی مانند عوامل محیطی، ترافیک، نوع جاده و رفتار رانندگی و غیره ناشی میشوند. [
۳۴]. محاسبه دقیق این عوامل برای تخمین دقیق مصرف انرژی خودرو بسیار مهم است [
۳۵].
سه روش غالب – مدل های تحلیلی، آماری و محاسباتی – برای محاسبه مصرف انرژی EV استفاده می شود. [
۳۶]. مدل های تحلیلی، که بر اساس دینامیک طولی وسیله نقلیه و نقشه های بازده موتور الکتریکی، برآورد تلفات و توان مورد نیاز برای غلبه بر نیروهای وسیله نقلیه است. [
۳۷,
۳۸]. مدلهای آماری بر تجزیه و تحلیل دادههای محرک و توسعه سیستمهای رگرسیون برای تعیین روابط بین مصرف انرژی و عوامل مؤثر متکی هستند. [
۳۹,
۴۰,
۴۱]. مدلهای محاسباتی با استفاده از شبکههای عصبی مصنوعی (ANN) مصرف انرژی را به عنوان تابعی از عوامل ورودی با استفاده از الگوریتمهای آموزشی تخمین میزنند. [
۴۲,
۴۳].
در این مطالعه، یک مدل دینامیکی EV جامع معرفی شده است که پارامترهای مختلفی را در مصرف انرژی در هر دو حوزه مکانیکی و الکتریکی تأثیر میگذارد. این مدل شامل پارامترهای آیرودینامیکی، ضریب اصطکاک لاستیک، مساحت جلوی خودرو، وزن، سرعت خودرو، سرعت باد، شیب جاده، و ویژگیهای موتور الکتریکی، مانند نقشههای بازده در محدوده سرعت و گشتاور، مبدلهای قدرت و باتریها است.
بخش ۲ طراحی و تجزیه و تحلیل مکانیکی خودرو را مورد بحث قرار می دهد و به دنبال آن مدل دینامیکی خودرو را با استفاده از یک رویکرد چند رشته ای اجرا می کند. این بخش همچنین به جزئیات مدل سازی موتور کششی و استراتژی کنترل آن می پردازد. علاوه بر این، الگوریتمی برای تخمین خطوط راندمان موتور نیز در جزئیات ارائه شده است
بخش ۲. این الگوریتم با توجه به حساسیت راندمان موتور به شرایط کاری، توسعه یافته و در مدل دینامیکی EV گنجانده شده است.
بخش ۳ نتایج به دست آمده از مدل را برای دو طرح خودروی مجزا ارائه می دهد. مدل شبیهسازی توسعهیافته امکان دستیابی سریع تمام دادههای مصرف انرژی یک EV را در یک محیط کامپیوتری با استفاده از دادههای جغرافیایی فراهم میکند. در نهایت، این مطالعه با خلاصه کردن بینشها و مفاهیم کلیدی به پایان میرسد.
۲٫ طراحی و مدلسازی خودروی الکتریکی
در این تحقیق، یک مدل EV با دقت ساخته شد تا محاسبه دقیق مصرف انرژی خودرو را تسهیل کند. این مدل به طور پیچیده از بلوک های متمایز از جمله بلوک محاسبه نیرو، بلوک مدل PMSM، بلوک نقشه کارایی PMSM و بلوک محاسبه مصرف انرژی خودرو تشکیل شده است. پارامترهای پارامونت حاکم بر ضریب درگ آیرودینامیکی خودرو (
، ناحیه جلویی (
و جرم وسیله نقلیه (
) از طریق تجزیه و تحلیل های دقیق انجام شده از طریق پلت فرم های نرم افزار CATIA و ANSYS Fluent به دست آمد. پارامترهای مشخص شده در
شکل ۲ به طور قابل توجهی به بلوک محاسبه نیرو کمک می کند، که به طور استراتژیک برای محاسبه نیروی کشش وسیله نقلیه استفاده می شود. این نیروی کشش، که یک جزء ضروری در مدل موتور الکتریکی و بلوک محاسبه مصرف انرژی بعدی است، به عنوان یک عامل محوری در محاسبات کل نگر مدل عمل می کند. علاوه بر این، PMSM به طور جامع با استفاده از پارامترهای به دست آمده از ANSYS Electromagnetics مدلسازی میشود و نمایش دقیقی را در مدل تضمین میکند. [
۲۵]. علاوه بر این، یک الگوریتم اختصاص داده شده به پیشبینی نقشه بازده استفاده میشود که تصویری جامع از ویژگیهای بازده موتور را تسهیل میکند، همانطور که در نشان داده شده است.
شکل ۲. این ادغام بلوک ها و ماژول های الگوریتمی یک چارچوب قوی را تشکیل می دهد که امکان برآورد دقیق مصرف انرژی خودرو را در مدل پیشنهادی EV فراهم می کند.
مصرف انرژی خودرو با توجه به داده های بازده PMSM، نیروی کشش خودرو، راندمان اینورتر و سرعت خودرو محاسبه می شود. انرژی مصرفی آنی خودرو از ظرفیت باتری کم می شود.
۲٫۱٫ آیرودینامیک و تحلیل خودرو
برای تجزیه و تحلیل جریان یک جسم، لازم است معادلات جرم، تکانه و انرژی جریان در تماس با آن جسم در شرایط مرزی فعلی حل شود. معادله (۱) معادله نیروی پسا آیرودینامیکی را نشان می دهد:
جایی که نیروی کشش آیرودینامیکی است، فشار محیط است، فشار اتمسفر است و مساحت جریان عمود بر جریان روی وسیله نقلیه است [۴۴].
بر اساس رابطه (۲) ضریب درگ آیرودینامیکی به صورت زیر بدست می آید:
جایی که ضریب درگ آیرودینامیکی است، چگالی هوا است، سرعت هوا است و قسمت جلوی وسیله نقلیه است [۴۴].
این
ضریب با توجه به شکل وسیله نقلیه و زاویه جریان باد تغییر می کند. بالا
ضریب نامطلوب است و برای رانندگی کارآمد باید تا حد امکان کاهش یابد. کاهش نیروی مقاومت منجر به کاهش مصرف انرژی می شود [
۴۴]. برای بهبود عملکرد خودرو باید آیرودینامیک آن مورد مطالعه قرار گیرد. جریان را می توان با استفاده از معادلات ناویر-استوکس ارائه شده در معادلات (۳) و (۴) تعریف کرد. [
۴۵]:
که در آن سرعت نشان می دهد اجزاء و تنش های آشفته را نشان می دهد. تجزیه و تحلیل با فرض جریان به عنوان حالت پایدار انجام شد. این معادلات را می توان با استفاده از روش معادله ناویر-استوکس با میانگین رینولد (RANS) حل کرد.
تمرکز این روش بر روی اثر آشفتگی بر خواص جریان متوسط است. همچنین معادله ناویر-استوکس را با میانگین زمان تقریبی حل می کند و نمی تواند تمام موقعیت های آشفته را حل کند.
مدلهای تنش رینولدز و رینولدز شناختهشدهترین مدلها در این روش حل هستند. این روش به محاسبات دقیق جریان و منابع محاسباتی قابل توجهی نیاز دارد. بنابراین، این رویکرد پایه اصلی مهندسی در محاسبات جریان در سی سال گذشته بوده است. [
۴۵].
در این مطالعه، قابل تحقق
مدل جریان آشفته از روش حل RANS انتخاب شد. معادلات حمل و نقل مدل شده برای
و
در مدل تحقق پذیر در معادلات (۵) و (۶) نشان داده شده است:
که در آن معادله (۷) است
همانطور که با دیگر
مدل ها، ویسکوزیته گردابی (
) از رابطه (۸) محاسبه شد.
تفاوت بین قابل تحقق
مدل و استاندارد و RNG
مدل ها این است که
دیگر ثابت نیست در عوض، می توان آن را از رابطه (۹) محاسبه کرد:
جایی که در رابطه (۱۰) به صورت زیر ارائه شده است:
جایی که میانگین نرخ چرخش تانسور است که در یک قاب مرجع متحرک با سرعت زاویه ای (). مدل ثابت می شود و در معادله (۱۱) ارائه شده است:
جایی که
می توان مشاهده کرد که
تابعی از میانگین نرخ کرنش و چرخش، سرعت زاویهای چرخش سیستم و میدانهای آشفتگی است.
و
).
در رابطه (۹) می توان نشان داد که مقدار استاندارد ۰٫۰۹ برای یک زیرلایه اینرسی در یک لایه مرزی تعادل را بازیابی می کند. ثابت های مدل به صورت زیر ارائه می شوند:
۲٫۲٫ مدلسازی دینامیک خودرو
مدل دینامیکی EV طراحی شده با استفاده از نرم افزار MATLAB–Simulink ایجاد شد. دادههای مسیر استفادهشده در مدل پویا با مسیری مطابقت دارد که بهطور خاص برای مسابقات نمونه اولیه خودروهای سبک وزن در Korfez، Izmit، از طریق نرمافزار Google Earth 9.0 بهدستآمده است.
شکل ۳یک و
شکل ۳b به ترتیب چیدمان مسیر مسابقه و مشخصات ارتفاعی آن را نشان می دهد.
۲٫۲٫۱٫ مدل PMSM
معادلات ولتاژ مدل PMSM در رابطه (۱۳) ارائه شده است. [
۴۶]:
جایی که ، ، ولتاژهای فاز خنثی هستند، ، ، کل شارهای سیم پیچ فاز هستند و ، ، مقاومت های سیم پیچ استاتور هستند [۴۶].
اگر کل شارهای متعلق به سیم پیچ ها به صورت باز نوشته شوند، معادله (۱۴) به صورت زیر است. [
۴۶]:
جایی که ، ، و خود القایی سیم پیچ های سه فاز هستند و ، ، ، ، ، و اندوکتانس متقابل بین هر سیم پیچ سه فاز هستند. علاوه بر این، ، ، و جریان سیم پیچ های سه فاز هستند و کل شار ایجاد شده توسط آهنرباها در سیم پیچ فاز است [۴۶].
به طور کلی، در PMSM، از آنجایی که آهنرباها روی سطح روتور قرار می گیرند و نفوذپذیری مغناطیسی آهنرباها نزدیک به نفوذپذیری مغناطیسی هوا است، می توان فرض کرد که موقعیت آهنرباها بر اندوکتانس فاز تأثیر نمی گذارد. علاوه بر این، اگر ساختار استاتور متقارن باشد، مقادیر خودالقایی و اندوکتانس متقابل را می توان ثابت در نظر گرفت. [
۴۶]. در این صورت اگر
=
=
=
،
=
=
=
و معادلات شار کل با معادلات ولتاژ جایگزین می شوند، سپس معادله ولتاژ فاز A به معادله (۱۵) تبدیل می شود. [
۴۶].
شکل ماتریسی معادلات مدل سازی PMSM در رابطه (۱۶) نشان داده شده است. [
۴۶].
توان الکترومغناطیسی منتقل شده به روتور را می توان به صورت معادله (۱۷) نشان داد. [
۴۶].
اگر فرض شود که تمام توان الکترومغناطیسی با نادیده گرفتن تلفات مکانیکی به انرژی جنبشی تبدیل می شود، توان مانند رابطه (۱۸) است. [
۴۶].
اگر (۱۷) و (۱۸) با هم ترکیب شوند، سرعت روتور را می توان از رابطه (۱۹) به دست آورد. [
۴۶].
از آنجایی که شکل موج EMF پشت موتور طراحی شده سینوسی است، تابع شکل موج ولتاژ حرکت به عنوان کسینوس تعریف می شود. ولتاژهای حرکت در رابطه (۲۰) نشان داده شده است. [
۴۶].
جزئیات مدل PMSM طراحی شده در نشان داده شده است
شکل ۴.
۲٫۲٫۲٫ نقشه کارایی PMSM
پارامترهای PMSM مدل شده متعلق به موتور طراحی شده و مورد استفاده برای وسیله نقلیه است که در رقابت چالش بهره وری TÜBİTAK شرکت می کند. با انتقال داده های نقشه عملکرد این موتور به محیط MATLAB–Simulink یک جدول جستجو ایجاد شد. از اینجا، راندمان موتور با توجه به سرعت و گشتاور موتور محاسبه شد. سپس از مقدار بازده در بلوک محاسبه مصرف انرژی استفاده شد. مدل نقشه کارایی PMSM در نشان داده شده است
شکل ۵.
۲٫۲٫۳٫ درایور PMSM
تکنیک کنترل جریان باند هیسترزیس برای هدایت PMSM استفاده می شود. در تکنیک کنترل جریان باند هیسترزیس، مقدار جریان مجبور می شود بین حدود خاصی باقی بماند، همانطور که در نشان داده شده است.
شکل ۶.
در مدل درایور PMSM، جریان مرجع مورد نیاز برای شرایط موجود با استفاده از کنترلر PI پیدا شد. با استفاده از بلوک باند هیسترزیس، سیگنال های محرک گیت مورد نیاز برای راه اندازی سوئیچ های بالایی فازهای A، B و C تولید می شوند. مدل درایور گیت در نشان داده شده است
شکل ۷.
۲٫۲٫۴٫ محاسبه نیرو
در این مدل پنج نیروی مختلف موثر بر خودرو محاسبه شده است. با توجه به رابطه (۲۱) نیروی اینرسی وسیله نقلیه محاسبه می شود [
۴۷]:
جایی که نیروی اینرسی خودرو را نشان می دهد، جرم خودرو را نشان می دهد و شتاب خودرو را نشان می دهد. در رابطه (۲۲)، نیروی گرانشی مؤثر بر وسیله نقلیه محاسبه شده است [۴۷]: جایی که نشان دهنده نیروی گرانشی و g نشان دهنده شتاب گرانشی است.
در رابطه (۲۳)، نیروی مقاومت غلتشی چرخ محاسبه شده است [
۴۷]:
جایی که نیروی مقاومت غلتشی چرخ ها را نشان می دهد [N]، ضریب مقاومت غلتشی چرخ را نشان می دهد و به زاویه شیب جاده اشاره دارد.
نیروی مقاومت باد در رابطه (۲۴) محاسبه می شود. [
۴۷].
با توجه به معادله (۲۴)
نیروی مقاومت باد را نشان می دهد [N]،
تراکم هوای خشک را در ۲۰ درجه سانتیگراد نشان می دهد [kg/m
3]،
ضریب درگ آیرودینامیکی را نشان می دهد،
منطقه جلوی خودرو را نشان می دهد [m
2]،
سرعت خودرو را نشان می دهد [m/s]، و
سرعت باد را نشان می دهد [m/s] [
47].
در معادله (۲۵)، بیان نیروی کشش موثر بر خودرو ارائه شده است [
۴۷].
۴٫ نتیجه گیری
هدف اصلی این مطالعه ارزیابی مصرف انرژی طرحهای خودروی برنامهریزیشده با استفاده از رویکرد چند رشتهای و استراتژی کنترل قبل از مرحله ساخت بود. این تحقیق شامل تجزیه و تحلیل دو طرح خودروی متمایز تحت شرایط خارجی یکنواخت مانند شیب جاده و سرعت باد بود. برای تخمین مصرف انرژی، مدلهای خودروی دینامیکی با استفاده از MATLAB-Simulink ساخته شدند که پارامترهای عملکرد حیاتی بهدستآمده از ANSYS Fluent، ANSYS Electronics و CATIA را شامل میشود. ، ، و .
توسعه مدلهای دقیق در MATLAB–Simulink شامل PMSM، الگوریتم کنترل آن، نقشه کارایی آن و نیروی اصطکاک معادل بود. با ارزیابی نیروهای مؤثر بر وسایل نقلیه، محاسباتی برای تعیین نیروی کشش مورد نیاز برای انتقال چرخ در مدلهای دینامیکی انجام شد. این نیروی کششی ارزیابی گشتاور بار موثر بر موتور و محاسبه متعاقب آن مصرف انرژی خودرو را تسهیل کرد. برای تنظیم سرعت خودرو، یک کنترلر PI بر اساس یک چرخه درایو مرجع از پیش تعریف شده در نظر گرفته شده برای سناریوی مسابقه مورد استفاده قرار گرفت. تجزیه و تحلیل دادههای گشتاور بار و سرعت استخراج دادههای راندمان لحظهای از نقشه بازده موتور را امکانپذیر کرد. این پارامترهای کارایی در مدل مصرف انرژی برای محاسبه کل انرژی مصرف شده توسط هر وسیله نقلیه ادغام شدند. پس از آن، مدل باتری مطابق با مصرف انرژی محاسبه شده تخلیه شد.
نتایج مطالعه مصرف انرژی ۱۱۲۴ وات ساعت را برای وسیله نقلیه اول و ۶۳۲ وات ساعت را برای وسیله نقلیه دوم نشان داد، که نشان دهنده کارایی برتر طراحی دوم است. در نتیجه خودروی دوم برای تولید انتخاب و در مسابقه مورد استفاده قرار گرفت. یافتههای تحقیق بر ضرورت یک تحلیل جامع و بینرشتهای برای ارزیابی عملکرد خودروهای الکتریکی سبک وزن، شامل عواملی مانند دینامیک سیالات، اصطکاک تایرها، ضرایب درگ، باتری، موتور، اینورتر و سایر اجزای پیشرانه تاکید میکنند. این مطالعه زمینه را برای تخمین مصرف انرژی در مرحله طراحی EV ها فراهم می کند. نکته قابل توجه، از آنجایی که این وسیله نقلیه برای مقاصد رقابتی در نظر گرفته شده بود، اجزای اضافی مصرف کننده انرژی مانند روشنایی و تهویه مطبوع عمدا حذف شدند. تحقیقات آینده با هدف ارزیابی تاثیر اجزای کمکی بر مصرف انرژی، کمک به بهبود مستمر در بهره وری خودرو خواهد بود. علاوه بر این، در مطالعات آتی، یکپارچهسازی ناوبری در حین رانندگی، پیشبینی محدوده زمان واقعی را با در نظر گرفتن ویژگیهای جغرافیایی مسیر وسیله نقلیه تسهیل میکند.