Revision history [back]
 | 1 | initial version |
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 metersSimilarly 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 metersSimilarly 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/microFSimilarly 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.
| | 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 metersSimilarly 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 metersSimilarly 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/microFSimilarly 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.
| | 3 | relocated to single answer |
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 metersSimilarly 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 metersSimilarly 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/microFSimilarly 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.
