First time here? We are a friendly community of Power Systems Engineers. Check out the FAQ!

Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

answered Mar 9 '12

amaity gravatar image

Now perform the spacing calculations.

D12 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D21 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D13 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters
D31 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters

Similarly compute for D14 ... D81, D23 ... D82, D34 ... D83, D45 ... D84, D56 ... D85, D67 ... D87, D78, D87.

Then compute to the vertical distances.

d12 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d21 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d13 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters
d31 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters

Similarly compute for d14 ... d81, d23 ... d82, d34 ... d83, d45 ... d84, d56 ... d85, d67 ... d87, d78, d87.

Calculate the Potential Coefficients.

P11 = 17.975109*math.log(2*y1av/GMrceq); # in km/microF
P22 = 17.975109*math.log(2*y2av/GMrceq); # in km/microF
P33 = 17.975109*math.log(2*y3av/GMrceq); # in km/microF
P44 = 17.975109*math.log(2*y4av/GMrceq); # in km/microF
P55 = 17.975109*math.log(2*y5av/GMrceq); # in km/microF
P66 = 17.975109*math.log(2*y6av/GMrceq); # in km/microF
P77 = 17.975109*math.log(2*y7av/GMre); # in km/microF
P88 = 17.975109*math.log(2*y8av/GMre); # in km/microF
P12 = 17.975109*math.log(D12/d12); # in km/microF
P21 = 17.975109*math.log(D12/d12); # in km/microF

Similarly compute for P13 ... P81, P23 ... P82, P34 ... P83, P45 ... P84, P56 ... P85, P67 ... P87, P78, P87.

Calculate the phase impedance.

Z11 = (Rc/nb)+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/(4*nb))+math.log(Derc/GMRceq));
Z22 = Z11; Z33 = Z11; Z44 = Z11; Z55 = Z11; Z66 = Z11;
Z77 = Re+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/4)+math.log(Derc/GMRe));
Z88 = Z77;
Z12 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);
Z21 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);

Similarly compute for Z13 ... Z81, Z23 ... Z82, Z34 ... Z83, Z45 ... Z84, Z56 ... Z85, Z67 ... Z87, Z78, Z87.

click to hide/show revision 2
No.2 Revision

Now perform the spacing calculations.calculations. Start off with the horizontal distances.

D12 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D21 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D13 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters
D31 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters

Similarly compute for D14 ... D81, D23 ... D82, D34 ... D83, D45 ... D84, D56 ... D85, D67 ... D87, D78, D87.

Then compute to the vertical distances.

d12 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d21 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d13 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters
d31 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters

Similarly compute for d14 ... d81, d23 ... d82, d34 ... d83, d45 ... d84, d56 ... d85, d67 ... d87, d78, d87.

Calculate the Potential Coefficients.

P11 = 17.975109*math.log(2*y1av/GMrceq); # in km/microF
P22 = 17.975109*math.log(2*y2av/GMrceq); # in km/microF
P33 = 17.975109*math.log(2*y3av/GMrceq); # in km/microF
P44 = 17.975109*math.log(2*y4av/GMrceq); # in km/microF
P55 = 17.975109*math.log(2*y5av/GMrceq); # in km/microF
P66 = 17.975109*math.log(2*y6av/GMrceq); # in km/microF
P77 = 17.975109*math.log(2*y7av/GMre); # in km/microF
P88 = 17.975109*math.log(2*y8av/GMre); # in km/microF
P12 = 17.975109*math.log(D12/d12); # in km/microF
P21 = 17.975109*math.log(D12/d12); # in km/microF

Similarly compute for P13 ... P81, P23 ... P82, P34 ... P83, P45 ... P84, P56 ... P85, P67 ... P87, P78, P87.

Calculate the phase impedance.

Z11 = (Rc/nb)+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/(4*nb))+math.log(Derc/GMRceq));
Z22 = Z11; Z33 = Z11; Z44 = Z11; Z55 = Z11; Z66 = Z11;
Z77 = Re+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/4)+math.log(Derc/GMRe));
Z88 = Z77;
Z12 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);
Z21 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);

Similarly compute for Z13 ... Z81, Z23 ... Z82, Z34 ... Z83, Z45 ... Z84, Z56 ... Z85, Z67 ... Z87, Z78, Z87.

click to hide/show revision 3
relocated to single answer

updated Mar 16 '12

JervisW gravatar image

Now perform the spacing calculations. Start off with the horizontal distances.relocated to single answer

D12 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D21 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D13 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters
D31 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters

Similarly compute for D14 ... D81, D23 ... D82, D34 ... D83, D45 ... D84, D56 ... D85, D67 ... D87, D78, D87.

Then compute the vertical distances.

d12 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d21 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d13 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters
d31 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters

Similarly compute for d14 ... d81, d23 ... d82, d34 ... d83, d45 ... d84, d56 ... d85, d67 ... d87, d78, d87.

Calculate the Potential Coefficients.

P11 = 17.975109*math.log(2*y1av/GMrceq); # in km/microF
P22 = 17.975109*math.log(2*y2av/GMrceq); # in km/microF
P33 = 17.975109*math.log(2*y3av/GMrceq); # in km/microF
P44 = 17.975109*math.log(2*y4av/GMrceq); # in km/microF
P55 = 17.975109*math.log(2*y5av/GMrceq); # in km/microF
P66 = 17.975109*math.log(2*y6av/GMrceq); # in km/microF
P77 = 17.975109*math.log(2*y7av/GMre); # in km/microF
P88 = 17.975109*math.log(2*y8av/GMre); # in km/microF
P12 = 17.975109*math.log(D12/d12); # in km/microF
P21 = 17.975109*math.log(D12/d12); # in km/microF

Similarly compute for P13 ... P81, P23 ... P82, P34 ... P83, P45 ... P84, P56 ... P85, P67 ... P87, P78, P87.

Calculate the phase impedance.

Z11 = (Rc/nb)+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/(4*nb))+math.log(Derc/GMRceq));
Z22 = Z11; Z33 = Z11; Z44 = Z11; Z55 = Z11; Z66 = Z11;
Z77 = Re+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/4)+math.log(Derc/GMRe));
Z88 = Z77;
Z12 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);
Z21 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);

Similarly compute for Z13 ... Z81, Z23 ... Z82, Z34 ... Z83, Z45 ... Z84, Z56 ... Z85, Z67 ... Z87, Z78, Z87.