Расчет электрического поля, создаваемого высоковольтными линиями электропередачи ОАО "Костромаэнерго"

program diplom;

uses crt,printer,graph;

const

{Eps0=8.85-12;}

con=1/(2*pi);

e=500;

c1='-------------------------------------------------------------------------------';

label

tabl;

var

st:integer;

g:array[1..2,1..e] of real;

x:array[1..3] of real; {массив координат проводников фаз}

y:array[1..3] of real; {массив координат проводников фаз}

Umod:array[1..3] of real;{массив модулей напряжений фаз}

Uarg:array[1..3] of real; {массив углов напряжений фаз}

Xm,Ym,r0:real; {координаты точки расчета и радиус проводников}

k,i,j:integer; {счетчики циклов}

Ur:array[1..3] of real; {массив действительной части комплексного напряжения}

Ui:array[1..3] of real; { массив мнимой части комплексного напряжения }

delta:array[1..3] of real; {промежуточная переменная}

delta1:array[1..3] of real;{ промежуточная переменная }

Cx:array[1..3] of real; { промежуточная переменная }

Cy:array[1..3] of real; { промежуточная переменная }

Exr,Eyr,Exi,Eyi,Ex,Ey:real;{действительная и мнимая части вертикальной и горизонтальной составляющих напряженности}

LK,LNA,DL:real; {LK-конечная точка расчета, LN- начальная точка расчета, DL-шаг}

Exmod,Eymod,Fix,Fiy:real; {модули вертикальной и горизонтальной составляющих, фазы напряженности}

A,B,Tmax,Tmin,Emax,Emin,Emaxkor2,h:real;

nf,Kei:integer; {nf число проводников в фазе}

DZ,D,rf,rekv:real;

procedure graphika;

var

grDriver: Integer;

grMode: Integer;

ErrCode: Integer;

begin

grDriver := Detect;

InitGraph(grDriver, grMode,'F:\BP\PROGI\D1\');

end;

begin

x[1]:=-11.5;x[2]:=0;x[3]:=11.5;y[1]:=8;y[2]:=8;y[3]:=8;

r0:=0.0124; rf:=0.088; nf:=3;

Umod[1]:=500; Umod[2]:=500; Umod[2]:=500;

Uarg[1]:=0; Uarg[2]:=120; Uarg[3]:=240;

h:=15; Xm:=0; Ym:=2;

LNA:=-20; LK:=20;

clrscr;

for i:=1 to 3 do

begin

Write('X',i,'=');

readln(x[i]);

Write('Y',i,'=');

readln(y[i]);

Write('Umod',i,'=');

readln(Umod[i]);

Write('Uarg',i,'=');

readln(Uarg[i]);

end;

Write('R0=');

readln(R0);

Write('Ln=');

readln(LNA);

Write('Lk=');

readln(LK);

Write('Y=');

readln(Ym);}

for i:=1 to 3 do

begin

Ur[i]:=Umod[i]*cos(Uarg[i]*Pi/180);

Ui[i]:=Umod[i]*sin(Uarg[i]*Pi/180);

writeln('Ur',i,'=',Ur[i]:10:10);

writeln('Ui',i,'=',Ui[i]:10:10);

end;

clrscr;

st:=round(e/15+0.5);

Xm:=1;

DL:=(LK-LNA)/(e-1);

for k:=1 to e do

begin

Xm:=LNA+(k-1)*DL;

DZ:=64/sqrt(50*0.0001);

for i:=1 to 3 do

begin

delta1[i]:=sqr(x[i]-Xm)+sqr(y[i]-Ym);

delta[i]:=sqr(x[i]-Xm)+sqr(y[i]+Ym+DZ);

Cx[i]:=(Xm-x[i])/delta1[i]-(Xm-x[i])/delta[i];

Cy[i]:=(Ym-y[i])/delta1[i]-(Ym+y[i]+DZ)/delta[i];

end;

Exr:=con*(Ur[1]*Cx[1]+Ur[2]*Cx[2]+Ur[3]*Cx[3]);

Exi:=con*(Ui[1]*Cx[1]+Ui[2]*Cx[2]+Ui[3]*Cx[3]);

Eyr:=con*(Ur[1]*Cy[1]+Ur[2]*Cy[2]+Ur[3]*Cy[3]);

Eyi:=con*(Ui[1]*Cy[1]+Ui[2]*Cy[2]+Ui[3]*Cy[3]);

Exmod:=sqrt(sqr(Exr)+sqr(Exi));

Eymod:=sqrt(sqr(Eyr)+sqr(Eyi));

Fix:=arctan(Exi/Exr);

Fiy:=arctan(Eyi/Eyr);

A:=sqr(Exmod)*cos(2*Fix)+sqr(Eymod*cos(2*Fiy));

B:=sqr(Exmod)*sin(2*Fix)+sqr(Eymod*sin(2*Fiy));

Tmax:=arctan((A+sqrt(sqr(A)+sqr(B)))/B);

Tmin:=arctan((A-sqrt(sqr(A)+sqr(B)))/B);

Emax:=sqrt(2)*(sqrt(sqr(Exmod)*sqr(sin(Tmax+Fix))+sqr(Eymod)*sqr(sin(Tmax+Fiy))));

Emin:=sqrt(2)*(sqrt(sqr(Exmod)*sqr(sin(Tmin+Fix))+sqr(Eymod)*sqr(sin(Tmin+Fiy))));

Emaxkor2:=Emax/sqrt(2);

g[1,k]:=Xm;

g[2,k]:=Emaxkor2;

---------------------------------

writeln('Exr=',Exr:20:20);

writeln('Exi=',Exi:20:20);

writeln('Eyr=',Eyr:20:20);

writeln('Eyi=',Eyi:20:20);

writeln('Exmod=',Exmod:20:20);

writeln('Eymod=',Eymod:20:20);

writeln('Fix=',Fix:20:20);

writeln('Fiy=',Fiy:20:20);

writeln('A=',A:2);

writeln('B=',B:20:20);

writeln('Tmax=',Tmax:20:20);

writeln('Tmin=',Tmin:20:20);

writeln('Emax/koren iz 2=',Emaxkor2:20:20);

readln;

---------------------------------

if k>1 then goto tabl;

clrscr;

writeln('Rashet elektricheskogo polia na visote',Ym:2:2,'m ot zemli.');

writeln('');

writeln(c1);

writeln('| X | Emax | Emin | Emax/koren iz 2 | Ex | Ey |');

writeln('| m | kV/m | kV/m | kV/m | kv/m | kV/m |');

writeln(c1);

tabl:

if k mod st =0 then

writeln('| ',Xm:9:3,' | ',Emax:9:3,' | ',Emin:9:3,' | ',Emaxkor2:15:3,' | ',Exmod:9:3,' | ',Eymod:9:3,' |');

end;

writeln(c1);

writeln('Konec rascheta.');

readln;

----------------------------------------------------------------

graphika;

setcolor(13);

line(round(g[1,1]*15+getmaxx/2),round(getmaxy/2)+180,round(20*15+getmaxx/2)+15,round(getmaxy/2)+180);

line(round(g[1,1]*15+getmaxx/2),round(-10*12+getmaxy/2)+180,round(20*15+getmaxx/2),round(-10*12+getmaxy/2)+180);

line(round(g[1,1]*15+getmaxx/2),round(-20*12+getmaxy/2)+180,round(20*15+getmaxx/2),round(-20*12+getmaxy/2)+180);

line(round(g[1,1]*15+getmaxx/2),round(-30*12+getmaxy/2)+180,round(20*15+getmaxx/2),round(-30*12+getmaxy/2)+180);

setcolor(14);

outtextxy(round(g[1,1]*15+getmaxx/2),round(getmaxy/2)+167,'0 kV/m');

outtextxy(round(g[1,1]*15+getmaxx/2),round(-10*12+getmaxy/2)+167,'10 kV/m');

outtextxy(round(g[1,1]*15+getmaxx/2),round(-20*12+getmaxy/2)+167,'20 kV/m');

outtextxy(round(g[1,1]*15+getmaxx/2),round(-30*12+getmaxy/2)+167,'30 kV/m');

setcolor(13);

line(round(g[1,1]*15+getmaxx/2),round(getmaxy/2)+180,round(g[1,1]*15+getmaxx/2),round(-30*12+getmaxy/2)+165);

line(round((g[1,1]+10)*15+getmaxx/2),round(getmaxy/2)+180,round((g[1,1]+10)*15+getmaxx/2),round(-30*12+getmaxy/2)+180);

line(round((g[1,1]+20)*15+getmaxx/2),round(getmaxy/2)+180,round((g[1,1]+20)*15+getmaxx/2),round(-30*12+getmaxy/2)+180);

line(round((g[1,1]+30)*15+getmaxx/2),round(getmaxy/2)+180,round((g[1,1]+30)*15+getmaxx/2),round(-30*12+getmaxy/2)+180);

line(round((g[1,1]+40)*15+getmaxx/2),round(getmaxy/2)+180,round((g[1,1]+40)*15+getmaxx/2),round(-30*12+getmaxy/2)+180);

setcolor(14);

outtextxy(round((g[1,1]+10)*15+getmaxx/2),round(getmaxy/2)+183,'-10');

outtextxy(round((g[1,1]+20)*15+getmaxx/2),round(getmaxy/2)+183,'0');

outtextxy(round((g[1,1]+30)*15+getmaxx/2),round(getmaxy/2)+183,'10');

outtextxy(round((g[1,1]+40)*15+getmaxx/2),round(getmaxy/2)+183,'20');

outtextxy(round((g[1,1])*15+getmaxx/2),round(getmaxy/2)+183,'-20');

setcolor(15);

for i:=2 to e do

line(round(g[1,i-1]*15+getmaxx/2),round(-g[2,i-1]*12+getmaxy/2)+180,round(g[1,i]*15+getmaxx/2),

round(-g[2,i]*12+getmaxy/2)+180);

readkey;

closegraph;

end.