diff options
Diffstat (limited to '3035/CH4/EX4.18/Ex4_18.sce')
| -rwxr-xr-x | 3035/CH4/EX4.18/Ex4_18.sce | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/3035/CH4/EX4.18/Ex4_18.sce b/3035/CH4/EX4.18/Ex4_18.sce new file mode 100755 index 000000000..9adde50f3 --- /dev/null +++ b/3035/CH4/EX4.18/Ex4_18.sce @@ -0,0 +1,34 @@ + +// Variable Declaration +A = 0.9*exp(%i*1.0*%pi/180) //Line constant +B = 143.0*exp(%i*84.5*%pi/180) //Line constant(ohm) +V_R = 220.0 //Receiving end voltage(kV) +V_S = 240.0 //Sending end voltage(kV) +P = 100.0 //Power(MVA) +pf = 0.8 //Power factor lagging + +a = phasemag(A)*%pi/180 //Phase angle of A(radian) +b = phasemag(B)*%pi/180 //Phase angle of B(radian) + +// Calculation Section +P_R = P * pf //Active power at receiving end(MW) +cos_b_delta = (P_R*abs(B)/(V_R*V_S))+(abs(A)*V_R/V_S)*cos(b-a) //cos(b-delta)[in radians] +delta_1 = (b - acos(cos_b_delta)) +Q_R = (V_R*V_S/abs(B))*sin(b-delta_1)-(abs(A)*V_R**2/abs(B))*sin(b-a) //Reactive power at receiving end(MVAR) +P_Re = P *(1-pf**2)**0.5 //Reactive power(MVAR) +rating = P_Re - Q_R //Rating of phase modifier(MVAR) + +delta_2 = b //Maximum power is received when delta = b +P_Rmax = (V_R*V_S/abs(B))-(abs(A)*V_R**2/abs(B))*cos(b-a) //Maximum power at receiving end(MW) +Q_R = -(abs(A/B)*V_R**2)*sin(b-a) //Reactive power at receive end(MVAR) +P_S = (V_S**2*abs(A/B))*cos(b-a)-(V_S*V_R/abs(B))*cos(b+delta_2) //Sending end power(MW) +n_line = (P_Rmax/P_S)*100 //Line efficiency(%) + +// Result Section +printf('Case(a) :') +printf('Rating of phase modifier = %.3f MVAR' ,rating) +printf('Power angle , delta = %.2f°' ,(delta_1*180/%pi)) +printf('\nCase(b) :') +printf('Maximum power at receive end , P_Rmax = %.2f MW' ,P_Rmax) +printf('Reactive power available , Q_R = %.2f MVAR' ,Q_R) +printf('Line efficiency = %.2f percent' ,n_line) |