پخش بار نیوتون رافسون در متلب

پخش بار نیوتون رافسون در متلب

با توجه به پیشرفت نرم افزارها، پخش بار یک سیستم چند باسه را می توان با نرم افزارهای متعددی نظیر: دیگسایلنت، ای تپ، متلب و … انجام داد و در تمام آنها نیز شاهد دقت بسیار بالا هستیم.

با توجه به اهمیت پخش بار، امروز برآن شدیم تا برای شما عزیزان یک پروژه پخش بار را جهت استفاده قرار دهیم.

در ابتدا اجازه دهید انواع روش های گرفتن پخش بار را به صورت تئوری برایتان مرور نماییم؛

روش های محاسبه پخش بار:

روش گوس سایدل

روش نیوتن رافسون

روش Decoupled

روش Fast Decoupled

روش پخش بار DC

 

روش گوس سایدل:

گاوس سایدل (Gauss–Seidel method) در جبر خطی عددی، روش تکراری است که برای حل دستگاه معادلات خطی استفاده می‌شود. نام این روش از روی ریاضی‌دانان آلمانی کارل فریدریش گاوس و فیلیپ لودویگ ون سایدل نهاده شده ‌است. اگرچه از این روش می‌توان در هر ماتریسی که دارای درایه قطری صفر نباشد استفاده کرد، اما فقط در صورتی همگرایی تضمین می‌شود که ماتریس مثبت معین یا قطری غالب باشد.

برای یک سیستم مربعی باn  معادلهٔ خطی و دارای مجهولx  داریم:

که در آن:

اگر A را به ماتریس پایین مثلثی و بالا مثلثی L* و U تجزیه کنیم:

معادلات خطی سیستم به شکل زیر بازنویسی خواهند شد:

روش گاوس سایدل از روش تکراری برای حل قسمت چپ عبارت، جهت به دست آوردن x بهره می‌برد، و بدین منظور از مقدار قبلی x در سمت راست عبارت استفاده می‌کند. می‌توان آن را به صورت زیر نوشت:

و با استفاده از خواص ماتریس مثلثی L* می‌توان x(k+1) را به صورت زیر به دست آورد:

در واقع برای محاسبه xi(k+1) به همه عناصر x(k) به جز خود آن به xi(k) دیگر نیاز خواهد شد.

محاسبات تا زمانی ادامه داده می‌شود تا در تکراری خاص به خطایی کمتر از مقدار مورد نظر برسیم.

 

روش نیوتن رافسون (روش مورد استفاده):

این روش نیز یکی از راه های تعین ریشهٔ یک تابع است.

فرض کنید تابعی (نمودار آبی) دارید که می‌خواهید ریشه (محل برخورد تابع با محور xها) را بیابید یا به اصطلاح آن را حل کنید. در روش نیوتن رافسون ابتدا x0 را به عنوان حدس اولیه وارد فرمول زیر می‌کنیم تا x1 بدست آید. به همین ترتیب ادامه می‌دهیم و این بار x1 را در فرمول قرار می‌دهیم.

به همین ترتیب خواهیم داشت:

هر چه تعداد دفعات تکرار بیشتر باشد x بدست آمده به ریشه نزدیک‌تر است.

 

روش Decoupled:

تجربه نشان می دهد که تغییرات P و δ بسیار زیاد بهم وابسته اند. همچنین Q و |V| نیز بهم وابسته اند پس می توان عناصر N و J ماتریس ژاکوبین را برابر صفر قرار داد.به این روش جداسازی یا Decoupled گفته می شود. معادلات توان اکتیو و راکتیو بصورت زیر در می آید:

27

در این روش سرعت محاسبات بطور قابل ملاحظه ای افزایش می یابد اما دقت محاسبات کاهش می یابد.

 

روش Fast Decoupled:

می توان روش Decoupled را برای تسریع در همگرایی باز تقریب زد و با قراردادن اختلاف زاویه بین شین ها برابر صفر معادلات را ساده تر کرد که به این روش Fast Decoupled می گویند. معادلات در این روش بصورت زیر در می آید:

28

 

روش پخش بار DC:

در این روش مقدار تقریبی توان خطوط محاسبه می شود. در این روش از مقاومت خطوط صرف نظر کرده و ولتاژ شین ها را برابر 1puو اختلاف زاویه بین شین ها را کوچک در نظر می گیریم. معادله توان در این روش بصورت زیر در می آید:

29

 

در ابتدا اجازه بدهید چند استاندارد رسمی برای حل سیستم های چند باسه را توضیح دهیم؛

شین های استفاده شده در محاسبات و سیستم معمولاً به سه نوع دسته بندی می شوند:

1) شین مرجع (Slack): این شین که به شین نیروگاهی معروف است به عنوان شین مرجع یا (Slack) در نظر گرفته می شود. زاویه در این شین، صفر درجه و ولتاژ معمولا یک پریونیت در نظر گرفته می شود. پس از محاسبات پخش بار، کمبود تولید و تلفات سیستم توسط این شین باید جبران شود. این شین را معمولا شین شماره 1 در نظر می گیریم و در این شین مقادیر |V| و δ معلوم و مقادیر P و Q مجهول می باشد.

در ادامه بخوانید  آموزش سیمولینک متلب

2) شین های کنترل شده: شین های دارای ژنراتور هستند و به شین های کنترل شده یا PV معروف اند، زیرا در این نوع از شین ها ولتاژ و توان اکتیو مقداری معلوم دارد. در انجام محاسبات، ما باید به دنبال پیدا کردن مقدار فاز و توان راکتیو این شین ها می باشیم.

3) شین های بار: در این نوع شین ها توان اکتیو و راکتیو مشخص است و ولتاژ و فاز قسمت مجهول است که باید محاسبه شود. این شین ها به شین های PQ نیز موسوم هستند.

معمولاً باس شماره یک باس مرجع و باس های 1 تا m باس های کنترل ولتاژ و باس های m+1 تا n باس های مصرفی شماره گذاری می شوند.

در یک سیستم قدرت برای تحلیل، پارامترهای مهمی از هر باس باید در نظر گرفته بشوند که عبارتند از؛

  • اندازه ولتاژ |V|

  • زاویه فاز δ

  • توان اکتیو P

  • توان راکتیو Q

 

اجازه بدهید به سراغ فایل شبیه سازی برویم:

در این پروژه ما به روش نیوتن – رافسون به تحلیل یک شبکه 33 باسه مشخص می پردازیم و در نهایت توان های اکتیو، راکتیو، ولتاژ، توان، زاویه ولتاژ و … را برای هر باس بدست می آوریم.

توضیحات روش استفاده شده در متلب:

در این روش به دلیل اینکه از 3 عدد ام فایل (m-file) استفاده کرده ایم پروژه را کمی پیچیده کرده است در حالی که روش ساده تر هم وجود دارد و آن استفاده از کتابخانه پیش فرض MatPower می باشد که بدون نیاز به نوشتن یک کد خط می توانید داده های باس های خود را (3 باسه، 5 باسه و … ) تشکیل داده و با وارد کردن داده های ضروری در نهایت به راحتی فشردن یک دکمه پخش بار را تحویل بگیرید.

در ادامه به زودی روش ذکر شده را توضیح خواهیم داد.

در ادامه بخوانید  تولباکس MATPOWER v6.0

 

در این پروژه شاهد 5 فایل می باشیم که در زیر برای شما کارکرد آنها را توضیح داده ایم؛

LFYBUS (مخفف لود فلو وای باس): کد تشکیل دهنده ماتریس ادمیتانس برای پخش بار

Busdata33: ماتریس اطلاعات باس برای کد LFYBUS

Linedata33: ماتریس اطلاعات خطوط برای کد LFYBUS

LFNEWTON: کد پخش بار اصلی نیوتن رافسون

BUSOUT: کد نمایش نتایج نهایی

 

 

در زیر برای شما کدهای قسمت های مختلف را قرار داده ایم که می توانید در نرم افزار متلب عینا وارد نمایید؛

LFYBUS:

LFNEWTON:

 

BUSOUT:

 

 

در ویدیو زیر شما نحوه Run کردن این پروژه را خواهید آموخت:

 

 

 

توجه:

به زودی برای شما روش پخش بار استاندارد با MatPower را آموزش خواهیم داد که در آن به راحتی می توانید داده های باس های تعریف شده را تغییر نمایید.

 

 

1415866183_Internet_Download_Manager     دانلود مستقیم فایل پخش بار نیوتون رافسون  | با حجم 1 مگابایت

1415866183_Internet_Download_Manager     دانلود مستقیم ویدیو آموزش پخش بار نیوتون رافسون در متلب – Full HD  | با حجم 12 مگابایت 

1415866190_698841-icon-114-lock-128     پسورد : www.poweren.ir

 

راستي! براي دريافت مطالب جديد در کانال تلگرام PowerEn عضو شويد.

تلگرام

 

پخش بار نیوتون رافسون در متلب
5 (100%) 1 vote
مهندس سیاه تیری
گرایش مورد علاقه ام ماشین های الکتریکی، بخصوص شار محورها هست - عاشق کار با نرم افزار مکسول هستم - به زودی ایده خودم رو که نیروگاه خانگی هست راه اندازی می کنم و تموم موفقیت های داشته و نداشتم رو مدیون کسی هستم که بدون هیچ چشم داشتی کنارم موند. و واقعا خوشحال می شم بتونم کمکتون کنم. دانشجوی کارشناسی ارشد برق-قدرت (ماشین های الکتریکی و الکترونیک قدرت) - دانشگاه صنعتی خواجه نصیر الدین طوسی

2
دیدگاه بگذارید

avatar
1 نظرها
1 پاسخ ها
1 فالورها
 
پر امتیازترین
پر بحث ترین
2 نویسندگان نظرات
مهندس سیاه تیریفلاحی جدیدترین نویسنده ها
  اطلاع رسانی با ایمیل  
جدیدترین قدیمی ترین محبوب ترین
اطلاع از
فلاحی
کاربر
فلاحی

سلام اگر 14 باسه داشته باشیم چطور ماتریس ژاکوپین رو حساب کنم ؟

گرایش رشته تحصیلی
قدرت

نظرت در مورد یه همراه همیشگی چیه!؟

 

 برای شروع این همراهی کافیه ایمیلت رو ثبت کنی - همین

عالی بود حالا از سایت لذت ببر