پروژه گمز2

سلام ....پروژه پخش بار ac برای ieee 33 bus با معادلات خطی رو داخل نرم افزار گمز میتونید پیاده سازی کنید؟ دیتا و معادلات رو خودم واستون میفرستم.این معادلات خطی پخش بار ac خطی برای سیستم ieee 33 bus هست که استفاده خواهد شد.مقدار Sbase برابر 100000 و Vbase برابر 12.665 تا generator وجود داره که در باس های 6 و 18 و 21 و 24 و 30 قرار دارن که همون Pg و Qg در معادلات 22a و 22b هستن.دیتای لود به صورت ساعتیه که به صورت زیر هست: Table E(t, b) "KW"          B1            B2            B3            B4            B*******;که Pg به صورت زیر به دست میاد: e\_Pg(t, i)..   Pg(t, i) =e= sum(B$Bus\_B(B, i), E(t, b)/Sbase);و Bus\_B هم به صورت زیر هست: Sets Bus\_B(B, i) "Mapping between B and Bus" / B1.6, B2.18, B3.21, B4.24, B5.30 /;مقدار Qg هم به صورت زیر هست: e\_Qg(t, i)..   Qg(t, i) =e= sum(B$Bus\_B(B, i), (E(t, b)\*tan(arccos(0.8)))/Sbase);دیتای Pd و Qd به این صورت هست: Table PB\_D(i, \*) "KW and KVAR"         Pd      Qd*******; و به صورت زیر در گمز پیاده سازی میشن: e\_Pd(t, i)..   P\_demand(t, i) =e= PB\_D(i, 'Pd');e\_Qd(t, i)..   Q\_demand(t, i) =e= PB\_D(i, 'Qd'); و دیتای مربوط به لاین ها به صورت زیر هست (که ستون limit اگه تغییر داده بشه مشکلی نیست!): Table LN(i, j, \*)             R             X             limit        G1           B*******; که R و X همون مقاومت و اندوکتانس لاین ها هست. limit حداکثر میزان توان عبوری از لاین ها و G1 و B1 به ترتیب مقدار کندوکتانس و سوسپتانس رو نشون میده که این دو در متلب محاسبه شده هست.باس Slack باس شماره 1 هست و محدوده ولتاژ از 0.95 تا 1.05 هست. زاویه ولتاژ هم از -p تا p هست.