$title Garver 6 Bus Network $onText Tochi Nwachukwu - Transmission Expansion Planning $offText Scalar alpha /1000/ Sbase /100 /; Variables f12,f14,f15,f23,f24,f26,f35,f46; Variables delta1,delta2,delta3,delta4,delta5,delta6; Variables g3,r2,r4,r5,g6,r1,Z; Integer Variables n12,n14,n15,n23,n24,n26,n35,n46; Equation const23, const24, const25, const26, const27, const28, const29, const30,const31,const32,const33,const34,const35,const36,const37,const38, const39,const40,const41,const42,const43,const44,const45,const46,const47, const48,const49,const50,const51,const52,const53; *********************************************************************** const23.. (-f12-f14-f15)*1/sbase =e= 0.30; const24.. (f12 -f23 - f26 + r2)*1/sbase =e= 2.40; const25.. (f23 -f35 + g3)*1/sbase =e=0; const26.. (f14 + f24-f46+r4)*1/sbase =e=1.60; const27.. (f15 + f35 + r5)*1/sbase =e=2.40; const28.. (f26 + f46 + g6)*1/sbase =e=0; const29.. (f12)*1/sbase -5/2*(1+n12)*(delta1-delta2) =e= 0; const30.. (f14)*1/sbase - 5/3*(1+n14)*(delta1-delta4) =e= 0; const31.. (f15)*1/sbase - 5*(1+n15)*(delta1-delta5) =e= 0; const32.. (f23)*1/sbase - 5*(1+n23)*(delta2-delta3) =e= 0; const33.. (f24)*1/sbase - 5/2*(1+n24)*(delta2-delta4) =e= 0; const34.. (f26)*1/sbase - 10/3*(n26)*(delta2-delta6) =e= 0; const35.. (f35)*1/sbase - 5*(1+n35)*(delta3-delta5) =e= 0; const36.. (f46)*1/sbase - 10/3*(n46)*(delta4-delta6) =e= 0; const37.. (f12)*1/sbase =l= (1+n12); const38.. (f12)*1/sbase =g= -(1+n12); const39.. (f14)*1/sbase =l= 0.8*(1+n14); const40.. (f14)*1/sbase =g= -0.8*(1+n14); const41.. f15*1/sbase =l= (1+n15); const42.. f15*1/sbase =g= -(1+n15); const43.. f23*1/sbase =l= (1+n23); const44.. f23*1/sbase =g= -(1+n23); const45.. f24*1/sbase =l= (1+n24); const46.. f24*1/sbase =g= -(1+n24); const47.. f26*1/sbase =l= n26; const48.. f26*1/sbase =g= -n26; const49.. f35*1/sbase =l= (1+n35); const50.. f35*1/sbase =g= -(1+n35); const51.. f46*1/sbase =l= n46; const52.. f46*1/sbase =g= -n46; const53.. Z =g= 40000*n12+ 60000*n14 + 20000*n15 + 20000*n23 + 40000*n24 + 30000*n26 + 20000*n35 + 30000*n46 + alpha*( r2*0.01 + r4*0.01 + r5*0.01); ; delta1.fx=0; delta2.up=3.142; delta2.lo=-3.142; delta3.up=3.142; delta3.lo=-3.142; delta4.up=3.142; delta4.lo=-3.142; delta5.up=3.142; delta5.lo=-3.142; delta6.up=3.142; delta6.lo=-3.142; g3.up = 1.25; g3.lo = 0; g6.up = 5.45; g6.lo =0; r1.up = 0.30; r1.lo =0; r2.up = 2.40; r2.lo =0; r4.up = 1.60; r4.lo =0; r5.up = 2.40; r5.lo =0; MODEL TEP_DET /ALL/; SOLVE TEP_DET USING MINLP MINIMIZING Z;