Программа PotyagAS

 

PROGRAM POTYAG; {$M $4000,0,0}

USES DOS,CRT,GRAPH;

CONST NN=25;

LABEL 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22;

TYPE Dl=ARRAY [1..NN] OF INTEGER;

Uk=ARRAY [1..NN] OF REAL;

Rds=ARRAY [1..NN] OF INTEGER;

DlRez=ARRAY [1..NN] OF INTEGER;

UkRez=ARRAY [1..NN] OF REAL;

DRd=ARRAY [1..NN] OF INTEGER;

PUk=ARRAY [1..NN] OF REAL;

VAR FV:text;

FO:text;

FC:text;

FR:text;

FD:text;

FK:text;

FT:text;

FL:text;

FE:text;

FA:text;

NAME:string;

S,st1,st2,st3:String;

i,j,k,n,m,z,d,di,x1,y1,x2,y2,amax,amin,Vrem,kk,Pr,Kdv,

p,Ns,Nr,m1,Wvf,ObKolVag,Vag4,Vag6,Vag8,ftek,kkk,vv,g:integer;

LM:Word;

bmin,bmax,

OsudsoVa,Put,Skor,SkorPr,UdSopLoc,UdSopSos,SopLoc,SopSos,SopPo,

UdSopPo,Tyag,UdTyag,SilOs,Ksi,LocMas,MasSos,MasPoez,DlSos,DlLoc,

DlPoez,Kob,kppr,Kst,RasEner,TyagaPr,SkorZad,UklUch,Usk,UskZad,

delta,Tyaga,DlinOb,PutProm,UdSopTor,TorSila,SkorTor,PrUkl,TorDel,

TorSilaUd,UklOb,PutTek,Prom,mu,Inf,Pnd,fn,kpd,cofi,lsras,Dsr,rs,

rrp,Unl,X0,Xs,Xrp,MasSosPr,Tau,Tau1,Tau2,f1,f2,U,Un,UnStr,sn,

alpha,betta,bettak,Kgen,bs,cs,ds,es,Mk,Mm,chis1,chis2,znam1,znam2,

Potok,Bd,Bz,Hd,kmu,Fdel,SumF,Imu,Irp,Idv,Mn,Kper,n2,Ptep,Kred,Dk,

delf,f2k,f2n,Kdelta,bzub,Md,Pst,PotokN,Krrp,f1max,Idv2,Ia,Ix,Itep,

TyagaTor,EnRec,LTor,Uw,LTMex,Ipr,Itep1,U1,U2,U3,U4,U5,I1,I2,I3,I4,

I5,Mp,delfx,Ux,f1x,Udv,Ptep2,Ptep1,Ptep3,Ptep4,Ptep5,Ptep6,Ptep7,

obkol,obkol2,obkol3,obkol4,obkol5,obkol6,topcikl,RasTop,Pdiz,Igen,

Itep2,Skor1,Skor2,zz,Mx,Md5,Md4,Md3,Md2,Md1,tyaga5,tyaga4,tyaga3,

tyaga2,tyaga1,Pte5,Pte4,Pte3,Pte2,Pte1,I11,I12,I13,I14,I15,Pte,

tyagav,Udvv,Ite2,Itep2v,PteS:real;

L:Dl;

Ukl:Uk;

Rad:Rds;

RezL:DlRez;

RezUkl:UkRez;

DlRds:DRd;

PrUklRez:PUk;

c:char;

{ Информационное окно }

PROCEDURE okno1;

Const

LU=#201;

RU=#187;

LD=#200;

RD=#188;

H=#205;

V=#186;

Var

k:integer;

begin

X1:=2; Y1:=2; X2:=78; Y2:=17;

TextBackground(LightBlue);

TextColor(Yellow);

Window(X1,Y1,X2,Y2);

Write(LU);

For k:=X1+1 to X2-1 do Write(H);

Write(RU);

For k:=Y1+1 to Y2-1 do

begin

GotoXY(1,k-Y1+1);

Write(V);

GotoXY(X2-X1+1,WhereY);

Write(V);

end;

Write(LD);

Window(X1,Y1,X2,Y2+1);

GotoXY(2,Y2-Y1+1);

For k:=X1+1 to X2-1 do Write(H);

Write(RD);

Window(X1+1,Y1+1,X2-1,Y2-1);

ClrScr;

end;

{ Диалоговое окно }

PROCEDURE okno2;

Const

LU=#201;

RU=#187;

LD=#200;

RD=#188;

H=#205;

V=#186;

Var

k:integer;

begin

X1:=2; Y1:=19; X2:=78; Y2:=24;

TextBackground(Cyan);

TextColor(Blue);

Window(X1,Y1,X2,Y2);

Write(LU);

For k:=X1+1 to X2-1 do Write(H);

Write(RU);

For k:=Y1+1 to Y2-1 do

begin

GotoXY(1,k-Y1+1);

Write(V);

GotoXY(X2-X1+1,WhereY);

Write(V);

end;

Write(LD);

Window(X1,Y1,X2,Y2+1);

GotoXY(2,Y2-Y1+1);

For k:=X1+1 to X2-1 do Write(H);

Write(RD);

Window(X1+1,Y1+1,X2-1,Y2-1);

ClrScr;

end;

{ Процедура диалога при вводе данных }

Procedure Okna(var s1,s2,s3:string);

begin

Okno2;

Writeln(s1);

Writeln(s2);

Writeln(s3);

Okno1;

end;

{ Процедура расчета параметров двигателя }

Procedure Raschet(Var U,f1,Skor,n2,betta:real);

begin

{ if (f1>fn) and (Skor<SkorZad-5) then begin f1:=f1+delf; U:=U+Un/5; end;}

if U>Un then U:=Un;

{ if (RezUkl[i]<0) and (Skor>SkorZad-2) then U:=U-Un/20;}

if U<=10 then U:=10;

U5:=U4;

U4:=U3;

U3:=U2;

U2:=U1;

U1:=U;

U:=(U1+U2+U3+U4+U5)/5;

(* write({FX,}f1:4:2,' ',U:5:2,' ',Skor:3:2,' ',Put:7:2,' ');*)

Potok:=0.96*U/(4.44*f1*Wvf*Kob);

Bd:=0.7073*p*Potok/(Dsr*lsras);

Bz:=3.1417*Bd*Dsr/(bzub*Ns*Kst);

if Bz>1.6 then

Hd:=-85650+38000*Bz

else

if (Bz>=1.2) and (Bz<=1.6) then

Hd:=Bz/(0.002982*Bz*Bz-0.00109*Bz*Bz*Bz-0.00037*Bz*Bz*Bz*Bz)

else

Hd:=610*Bz;

kmu:=Hd/(610*Bz);

Fdel:=0.8*Bd*delta*Kdelta*1000;

SumF:=2*Fdel*kmu;

Imu:=m1*SumF/(0.9*p*Wvf*Kob);

if f1>fn then

Md:=Potok*Potok*m1*4.44*4.44*Wvf*Wvf*Kob*Kob*fn*(Rrp/betta)/

(6.2834*(Xrp*Xrp+Rrp*Rrp/(betta*betta))) else Md:=Mn;

if U<=10 then Md:=50;

if RezUkl[i]<0 then Md:=50;

Irp:=1.4142*Md/(p*m1*Wvf*Kob*Potok*0.985);

Ia:=Pst*(Potok/PotokN)*(Potok/PotokN)*(f1/fn)*sqrt(f1/fn)/(m1*Un)+Irp*0.985;

Ix:=Imu+Irp*0.174;

Idv:=sqrt(Ia*Ia+Ix*Ix);

if (RezUkl[i]<0) and (Skor>SkorZad-2) then Idv:=10;

cofi:=Ia/Idv;

Itep:=m1*Idv*Kdv{10}/(2*cofi*Krrp);

{ Itep2:=(Ipr+Itep)/4;}

if Itep2>10000 then Itep2:=10000;

Ipr:=Itep2;

Itep2:=(Ipr+Itep)/4;

if Itep2<0 then Itep2:=0;

PteS:=U*Itep2/(Kgen*1000);

Tyaga:=mu*Kdv*2{*10}*Md/Dk;

UdTyag:=Tyaga/(LocMas+MasSos);

UdSopLoc:=19.64+0.098*Skor+0.0029*Skor*Skor;

OsUdSoVa:=(6.87+(29.43+0.98*Skor+0.0245*Skor*Skor)/21.4)*

Vag4/ObKolVag+(6.87+(78.48+0.9*Skor+0.0206*Skor*Skor)/21.4)

*Vag6/ObKolVag+(6.87+(58.9+0.373*Skor+0.0206*Skor*Skor)/

21.4)*Vag8/ObKolVag;

SopLoc:=UdSopLoc*LocMas;

SopSos:=OsUdSoVa*MasSos;

SopPo:=SopLoc+SopSos;

UdSopPo:=SopPo/(LocMas+MasSos);

ksi:=0.28+3/(50+20*Skor)-0.0007*Skor;

TyagaPr:=ksi*LocMas*9810/(LocMas+MasSos);

SkorPr:=Skor;

if Idv>15*Inf then Skor:=SkorPr

else

Skor:=Skor+0.0034*(UdTyag-UdSopPo-9.81*RezUkl[i]-PrUklRez[i]*9.81+

TyagaTor);

Put:=Put+(SkorPr+Skor)*10/(2*36);

n2:=Skor*mu/(0.188*Dk);

if n2<=0 then n2:=0;

if Skor<=0 then Skor:=0;

f2:=f1-n2*p/60;

bettak:=rrp*(1+Tau1)/(Xs*(1+Tau2)+Xrp);

f2k:=bettak*fn;

betta:=f2/fn;

Ptep:=PteS/2;

case g of

100: begin

if Itep<0 then Ptep:=200; if Ptep>2000 then Ptep:=2000;

Ptep2:=Ptep*Ptep; Ptep3:=Ptep*Ptep2; Ptep4:=Ptep*Ptep3;

Ptep5:=Ptep*Ptep4; Ptep6:=Ptep*Ptep5; Ptep7:=Ptep*Ptep6;

obkol:=3.6606e-19*Ptep7-2.9586e-15*Ptep6+9.7182e-12*Ptep5-

1.6666e-8*Ptep4+1.5933e-5*Ptep3-0.0083656*Ptep2+

2.4169*Ptep+157.01;

if obkol<200 then begin Ptep:=200; obkol:=200; end;

obkol2:=obkol*obkol; obkol3:=obkol*obkol2;

obkol4:=obkol*obkol3; obkol5:=obkol*obkol4;

topcikl:=(-1.4498e-12*obkol5+4.6782e-9*obkol4-5.792e-6

*obkol3+0.0034975*obkol2-1.1101*obkol+322.85)/3600;

if topcikl<0 then topcikl:=0;

end;

200: begin

if Itep2<0 then Ptep:=150; if Ptep>1350 then Ptep:=1350;

Ptep2:=Ptep*Ptep; Ptep3:=Ptep*Ptep2; Ptep4:=Ptep*Ptep3;

Ptep5:=Ptep*Ptep4; Ptep6:=Ptep*Ptep5; Ptep7:=Ptep*Ptep6;

obkol:=2.5808e-18*Ptep7-1.5114e-14*Ptep6+3.655e-11*Ptep5-

4.7202e-8*Ptep4+3.5151e-5*Ptep3-0.01518*Ptep2+

3.8756*Ptep-19.226;

if obkol<300 then begin Ptep:=150; obkol:=300; end;

obkol2:=obkol*obkol; obkol3:=obkol*obkol2;

obkol4:=obkol*obkol3; obkol5:=obkol*obkol4;

obkol6:=obkol*obkol5;

topcikl:=(3.4165e-17*obkol6-1.1263e-13*obkol5+1.7834e-9*obkol4-

4.1611e-6*obkol3+0.004023*obkol2-1.8692*obkol+501.01)/3600;

if topcikl<0 then topcikl:=0;

end;

300: begin

if Itep2<0 then Ptep:=100; if Ptep>909 then Ptep:=909;

Ptep2:=Ptep*Ptep; Ptep3:=Ptep*Ptep2; Ptep4:=Ptep*Ptep3;

Ptep5:=Ptep*Ptep4; Ptep6:=Ptep*Ptep5; Ptep7:=Ptep*Ptep6;

obkol:=-2.2502e-12*Ptep5+6.3453e-9*Ptep4-6.3992e-6*Ptep3+

0.002872*Ptep2-0.053651*Ptep+295.97;

if obkol<200 then begin Ptep:=100; obkol:=200; end;

obkol2:=obkol*obkol; obkol3:=obkol*obkol2;

obkol4:=obkol*obkol3; obkol5:=obkol*obkol4;

topcikl:=(-2.9687e-12*obkol5+8.6719e-9*obkol4-1.0107e-5*obkol3

+0.0060552*obkol2-1.9864*obkol+528.22)/3600;

if topcikl<0 then topcikl:=0;

end;

400: begin

if Itep2<0 then Ptep:=100; if Ptep>588 then Ptep:=588;

Ptep2:=Ptep*Ptep; Ptep3:=Ptep*Ptep2; Ptep4:=Ptep*Ptep3;

Ptep5:=Ptep*Ptep4; Ptep6:=Ptep*Ptep5; Ptep7:=Ptep*Ptep6;

obkol:=3.2615e-15*Ptep7-8.1075e-12*Ptep6+8.2991e-9*Ptep5-

4.5065e-6*Ptep4+0.0013916*Ptep3-0.24249*Ptep2+

23.171*Ptep-516.12;

if obkol<300 then begin Ptep:=100; obkol:=300; end;

obkol2:=obkol*obkol; obkol3:=obkol*obkol2;

obkol4:=obkol*obkol3; obkol5:=obkol*obkol4;

topcikl:=(9.779e-13*obkol5-3.3931e-9*obkol4+4.2183e-6*obkol3

-0.0020784*obkol2+0.16168*obkol+263.43)/3600;

if topcikl<0 then topcikl:=0;

end;

end;

(* writeln({FX,}f2:4:2,' ',cofi:2:4,' ',Potok:2:4,' ',Bd:2:4,' ',Bz:2:4,' ',

kmu:2:4,' ',Imu:5:2,' ',Irp:5:2,' ',Idv:6:2,' ',Iel:7:2,' ',

Md:6:2,' ',UdTyag:5:2,' ',UdSopPo:5:2,' ',f2k:3:2,' ',

n2:4:2,' ',TyagaPr:5:2);

writeln(' U=',U:7:2); readln;*)

 

end;

 

{ Основная программа }

BEGIN

TextBackGround(Black);

LM:=LastMode;

21: ClrScr;

{ Считывание данных поезда }

assign(fc,'C:\PoDizAs\DanPoezd.dat');

reset(fc);

readln(fc,ObKolVag,Vag8,Vag6,Vag4);

readln(fc,MasSos,SkorZad,UskZad);

Close(fc);

{ Считывание данных электровоза }

assign(fe,'C:\PoDizAs\Electra.dat');

reset(fe);

readln(fe,mu,Kdv,DlLoc,LocMas);

close(fe);

{ Считывание данных электродвигателя }

assign(fa,'C:\PoDizAs\ElDvig.dat');

reset(fa);

readln(fa,Unl,Inf,Pnd,fn,p,Ns,Dk,sn);

readln(fa,kpd,cofi,lsras,Dsr,bzub,rs,rrp);

readln(fa,Xrp,m1,delta,Wvf,Xs,X0,f2n,Mn);

{ Определение длины поезда }

DlSos:=20.5*Vag8+17*Vag6+14.6*Vag4;

DlPoez:=DlLoc+DlSos;

{ Считывание профиля пути }

st1:=' Нажмите клавишу, чтобы осуществить движение: ';

st2:=' <Y> - по неспрямленному профилю пути; ';

st3:=' <N> - по спрямленному профилю пути';

Okna(st1,st2,st3);

{ Открытие файла профиля пути }

repeat

c:=readkey

until (c in ['y','Y','n','N']);

case c of

'y','Y': begin

assign(fv,'C:\PoDizAs\profil.dat');

reset(fv);

readln(fv,m);

i:=1; kkk:=1;

while not EOF(fv) and (i<=m) do

begin

readln(fv,l[i],rad[i],DlRds[i],ukl[i]);

RezL[i]:=l[i]; RezUkl[i]:=Ukl[i];

If DlRds[i]<=0

then PrUklRez[i]:=0.0

else

if DlPoez<=DlRds[i]

then PrUklRez[i]:=700/rad[i]

else

PrUklRez[i]:=(700/Rad[i]*DlRds[i])/DlPoez;

inc(i)

end;

close(fv);

end;

'n','N': begin

assign(fo,'C:\PoDizAs\SprProfA.dat');

reset(fo);

readln(fo,m);

i:=1; kkk:=0;

while not EOF(fo) and (i<=m+1) do

begin

readln(fo,RezL[i],RezUkl[i],PrUklRez[i]);

inc(i)

end;

close(fo);

end;

end;

Okno2;

Writeln(' Укажите мощность используемого дизеля ');

Writeln(' Для этого нажмите: 1 - для выбора дизеля 1Д80Б ');

Writeln(' 2 - 4Д80Д, 3 - 4Д80Б, 4 - 10Д80А ');

repeat c:=readkey

until (c in ['1','2','3','4']);

case c of

'1': begin

writeln(' Выбран дизель типа 1Д80Б ');

writeln(' Номинальная мощность дизеля 2080 кВт ');

writeln(' Максимальные обороты коленчатого вала 850 об/мин ');

g:=100;

end;

'2': begin

writeln(' Выбран дизель типа 4Д80Д ');

writeln(' Номинальная мощность дизеля 1390 кВт ');

writeln(' Максимальные обороты коленчатого вала 750 об/мин ');

g:=200;

end;

'3': begin

writeln(' Выбран дизель типа 4Д80Б ');

writeln(' Номинальная мощность дизеля 909 кВт ');

writeln(' Максимальные обороты коленчатого вала 750 об/мин ');

g:=300;

end;

'4': begin

writeln(' Выбран дизель типа 10Д80А ');

writeln(' Номинальная мощность дизеля 588 кВт ');

writeln(' Максимальные обороты коленчатого вала 1000 об/мин ');

g:=400;

end;

end;

readln;

{ Расчет скорости, пути и времени движения }

Writeln;

Writeln(' Результаты расчета');

Writeln(' Для прекращения счета и выхода в программу',

' нажмите клавишу ВВОД');

Readln;

Okno1;

kk:=0;

20: If kk=0 then

begin

Assign(FD,'C:\PoDizAs\RezPoezA.dat');

Assign(FR,'C:\PoDizAs\rezulta.dat');

Rewrite(FD);

Rewrite(FR);

Okno2;

Writeln;

Writeln(' Расчет параметров',

' движения');

Writeln(' Движение поезда от',

' станции А');

Okno1;

{ end

else

begin

Close(FD);

Close(FR);

Assign(FK,'C:\PoDizAs\RezPoezB.dat');

Assign(FT,'C:\PoDizAs\rezultb.dat');

Assign(FL,'C:\PoDizAs\sprprofb.dat');

Rewrite(FK);

Rewrite(FT);

Rewrite(FL);

for i:=1 to m div 2 do

begin

Pr:=RezL[i]; RezL[i]:=RezL[m+1-i];

RezL[m+1-i]:=Pr; Prom:=RezUkl[i];

RezUkl[i]:=-RezUkl[m+1-i];

RezUkl[m+1-i]:=-Prom; Prom:=PrUklRez[i];

PrUklRez[i]:=PrUklRez[m+1-i];

PrUklRez[m+1-i]:=Prom;

end;

if ((m) mod 2)<>0 then

RezUkl[((m) div 2)+1]:=-RezUkl[((m) div 2)+1];

writeln(FL,m);

For i:=1 to m do

writeln(FL,RezL[i]:12,RezUkl[i]:10:2,

PrUklRez[i]:10:2);

Okno2;

Writeln;

Writeln(' Расчет параметров',

' движения');

Writeln(' Движение поезда от',

' подстанции Б');

Okno1; }

end;

Vrem:=0; Put:=0; Skor:=0.01; SilOs:=17.5; Tyaga:=0; RasEner:=0; U:=0;

EnRec:=0.0; DlinOb:=0; UklOb:=0; UklUch:=50; Put:=0.0; Ipr:=0.0; vv:=0;

For i:=1 to m+2 do DlinOb:=DlinOb+RezL[i]; Md:=0; f1:=0; TyagaTor:=0;

i:=1;

PutProm:=RezL[i];

{ Обмоточный коэффициент }

Kob:=0.91;

{ Кпд генератора }

Kgen:=0.93;

{ Кпд преобразователя частоты }

Krrp:=0.9;

{ Фазное напряжение двигателя }

Un:=Unl/1.73205;

{ Номинальный магнитный поток двигателя }

PotokN:=0.109;

{ Коэффициент заполнения стали }

Kst:=0.95;

{ Коэффициент Картера }

Kdelta:=1.38;

{ Потери в стали статора }

Pst:=900;

{ Предельная масса состава }

MasSosPr:=15000;

{ Расчет вспомогательных переменных }

Tau1:=Xs/X0;

Tau2:=Xrp/X0;

Tau:=Tau1+Tau2+Tau1*Tau2;

bs:=rs*(1+Tau2);

cs:=X0*Tau;

ds:=rs/X0;

es:=1+Tau1;

chis1:=cs*cs*sn/rrp+es*es*rrp/sn;

znam1:=(bs*bs+cs*cs)*sn/rrp+(ds*ds+es*es)*rrp/sn+2*rs;

UnStr:=Un*sqrt(chis1/znam1);

Skor:=0.02;

UdSopLoc:=19.64+0.098*Skor+0.0029*Skor*Skor;

OsUdSoVa:=(6.87+(29.43+0.98*Skor+0.0245*Skor*Skor)/21.4)*

Vag4/ObKolVag+(6.87+(78.48+0.9*Skor+0.0206*Skor*Skor)/21.4)

*Vag6/ObKolVag+(6.87+(58.9+0.373*Skor+0.0206*Skor*Skor)/

21.4)*Vag8/ObKolVag;

SopLoc:=UdSopLoc*LocMas;

SopSos:=OsUdSoVa*MasSosPr;

SopPo:=SopLoc+SopSos;

Md:=SopPo*Dk/(16*2*mu);

f1:=0.2;

Idv:=8000;

While Idv>2.5*Inf do

begin

alpha:=f1/fn; Mp:=Mn; Md:=Mn;

betta:=alpha;

chis2:=(bs*bs+cs*cs*alpha*alpha)*betta/rrp+

(ds*ds+es*es*alpha*alpha)*rrp/betta+2*rs*alpha;

znam2:=(bs*bs+cs*cs)*betta/rrp+(ds*ds+es*es)*rrp/betta+2*rs;

U:=UnStr*sqrt((chis2/znam2)*Md/Mp);

f2:=betta*fn;

Raschet(U,f1,Skor,n2,betta);

Mp:=Mp-Mp/8;

f1:=f1+0.1;

Skor:=0;

end;

LTor:=2000{0.006*SkorZad*(LocMas+MasSos)*8/6.5Kdv};

Repeat

if (DlinOb-LTor)>Put then

begin

If Put>PutProm then

begin

i:=i+1;

PutProm:=PutProm+RezL[i];

UklUch:=UklOb;

end;

PutTek:=Put+RezL[i]-PutProm;

UklOb:=UklUch+RezUkl[i]*PutTek/1000;

ClrScr;

if f1<=fn then

begin

delf:=0.957*UskZad*mu/Dk;

if Skor>1 then

begin

Md:=Md+Mn/5;

if Md>Mn then Md:=Mn;

end else Md:=Mn;

Ux:=U; f1x:=f1; delfx:=delf;

if Skor>SkorZad+2 then

begin

delfx:=delf;

delf:=0;

Ux:=U; f1x:=f1; vv:=0;

Skor2:=Skor1;

Skor1:=Skor;

if Skor2>Skor1 then

begin

zz:=zz+10;

Md:=Md/zz;

Mx:=Md;

end

else

begin

zz:=zz-10;

if zz<=1 then zz:=1;

Md:=Mx/zz;

end;

{ Md5:=Md4;

Md4:=Md3;

Md3:=Md2;

Md2:=Md1;

Md1:=Md;

Md:=(Md5+Md4+Md3+Md2+Md1)/5; }

if Skor>SkorZad+3 then

begin

TorDel:=PutProm-Put;

if TorDel>1010 then TorDel:=1000;

if TorDel<900 then TorDel:=1000;

TyagaTor:=-TyagaTor-abs((-40.85*(Skor*Skor-SkorZad*SkorZad)/TorDel

-(UdSopPo+PrUklRez[i]*9.81+RezUkl[i]*9.81)));

end;

end

else

begin

if vv=0 then begin U:=Ux; f1:=f1x; delf:=delfx; end;

vv:=vv+1;

if delf>2*f2n then delf:=2*f2n;

f1:=n2*p/60+delf;

TyagaTor:=0;

end;

alpha:=f1/fn;

betta:=f2/fn;

chis2:=abs((bs*bs+cs*cs*alpha*alpha)*betta/rrp+

(ds*ds+es*es*alpha*alpha)*rrp/betta+2*rs*alpha);

znam2:=abs((bs*bs+cs*cs)*betta/rrp+(ds*ds+es*es)*rrp/betta+2*rs);

U:=UnStr*sqrt(chis2/znam2)*sqrt(abs(Md/Mn));

if U>Un then U:=Un;

Raschet(U,f1,Skor,n2,betta);

end

else

begin

delf:=0.957*UskZad*mu/Dk;

{ Md:=Potok*Potok*m1*4.44*4.44*Wvf*Wvf*Kob*Kob*fn*(Rrp/betta)/

(6.2834*(Xrp*Xrp+Rrp*Rrp/(betta*betta)));}

if Skor>SkorZad+2 then

begin

delfx:=delf;

delf:=0;

Ux:=U; f1x:=f1; vv:=1;

Skor2:=Skor1;

Skor1:=Skor;

if Skor2>Skor1 then

begin

zz:=zz+10;

Md:=Md/zz;

Mx:=Md;

end

else

begin

zz:=zz-10;

if zz<=1 then zz:=1;

Md:=Mx/zz;

end;

{ Md5:=Md4;

Md4:=Md3;

Md3:=Md2;

Md2:=Md1;

Md1:=Md;

Md:=(Md5+Md4+Md3+Md2+Md1)/5;

if (RezUkl[i]<0) and (Iel2<0) then U:=U+Un/10;}

if Skor>SkorZad+5 then

begin

TorDel:=PutProm-Put;

if TorDel>1010 then TorDel:=1000;

if TorDel<900 then TorDel:=1000;

TyagaTor:=-TyagaTor-abs((-40.85*(Skor*Skor-SkorZad*SkorZad)/TorDel

-(UdSopPo+PrUklRez[i]*9.81+RezUkl[i]

*9.81)));

end;

end

else

begin

if vv=0 then begin U:=Ux; f1:=f1x; delf:=0.957*UskZad*mu/Dk; end

else delf:=0.957*UskZad*mu/Dk;

vv:=vv+1;

delf:=0.957*UskZad*mu/Dk;

if delf>{2*}f2k{n} then delf:={2*}f2k{n};

f1:=n2*p/60+{3*}delf;

TyagaTor:=0;

if Skor<SkorZad-5 then begin U:=U+Un/10; delf:=0.957*UskZad*mu/Dk; end;

if U>Un then U:=Un;

end;

if UdTyag>0.95*TyagaPr then delf:=-delf/6;

if Md>Mn then Md:=Mn;

alpha:=f1/fn;

betta:=f2/fn;

bettak:=rrp*sqrt((ds*ds+es*es*alpha*alpha)/(bs*bs+cs*cs*alpha*alpha));

f2k:=bettak*fn;

if Ptep>2100 then f1:=f1-delf else f1:=f1+delf;

Raschet(U,f1,Skor,n2,betta);

end;

Okno1;

Writeln(#10);

write(' Время движения ');

TextColor(LightGray);

Writeln(Vrem+1:6,' c;');

TextColor(Yellow);

Write(' Движение по ');

TextColor(LightGray);

Write(i);

TextColor(Yellow);

If kkk>=1 then

Write(' участку общей длиной ')

else

Write(' спрямленному участку общей длиной ');

TextColor(LightGray);

Writeln(PutProm:6:2,' м');

TextColor(Yellow);

Write(' Путь пройденный поездом ');

TextColor(LightGray);

Writeln(Put:8:2,' м;');

TextColor(Yellow);

Write(' Скорость ');

TextColor(LightGray);

Writeln(Skor:4:1,' км/ч;');

TextColor(Yellow);

U5:=U4;

U4:=U3;

U3:=U2;

U2:=U1;

U1:=U;

Udvv:=(U1+U2+U3+U4+U5)/5;

Write(' Напряжение на электродвигателе ');

TextColor(LightGray);

Writeln(Udvv:4:2,' B;');

TextColor(Yellow);

Write(' Частота тока статора ');

TextColor(LightGray);

Writeln(f1:4:2,' Гц;');

TextColor(Yellow);

if Ptes>4400 then PteS:=4400;

Pte5:=Pte4;

Pte4:=Pte3;

Pte3:=Pte2;

Pte2:=Pte1;

Pte1:=PteS;

Pte:=(Pte5+Pte4+Pte3+Pte2+Pte1)/5;

if Pte<=50 then Pte:=50;

Write(' Мощность дизеля тепловоза ');

TextColor(LightGray);

writeln(Pte:4:2,' кВт');

TextColor(Yellow);

if Itep2>5000 then Itep2:=5000;

I5:=I4;

I4:=I3;

I3:=I2;

I2:=I1;

I1:=Itep2;

Itep2v:=(I1+I2+I3+I4+I5)/10;

if Itep2v<=0 then Itep2v:=0;

{ if Itep2v>2000 then Itep2v:=2000;}

if Itep2v>5000 then Itep2v:=5000 else Itep2v:=Itep2;

Write(' Ток тягового генератора ');

TextColor(LightGray);

Writeln(Itep2v:8:0,' A');

TextColor(Yellow);

Write(' Цикловая подача топлива ');

TextColor(LightGray);

Writeln(topcikl:2:6,' г/КВт с');

TextColor(Yellow);

RasTop:=RasTop+topcikl*Ptep;

Write(' Расход топлива дизелем тепловоза ');

TextColor(LightGray);

Writeln(RasTop/1000:12:2,' кг');

TextColor(Yellow);

Write(' Обороты коленчатого вала дизеля тепловоза ');

TextColor(LightGray);

Writeln(obkol:12:2,' об/мин');

TextColor(Yellow);

Vrem:=Vrem+1;

if tyaga<20 then tyaga:=20;

tyaga5:=tyaga4;

tyaga4:=tyaga3;

tyaga3:=tyaga2;

tyaga2:=tyaga1;

tyaga1:=tyaga;

tyagav:=(tyaga5+tyaga4+tyaga3+tyaga2+tyaga1)/5;

If kk=0 then

begin

Writeln(FD,Vrem:7,Itep2:15:2,Put:16:2);

Delay(50);

Writeln(FR,Vrem:7,Skor*35:16:2,Put*0.1:16:2,UklOb*20:16:2,

Itep2v*2:16:2,Pte*1.8:16:2,tyagav/42:16:2,RasTop/16:16:2);

Delay(50);

{ end

else

begin

Writeln(FD,Vrem:7,IEl2:15:2,Put:16:2);

Delay(50);

Writeln(FR,Vrem:7,Skor*35:16:2,Put*0.1:16:2,UklOb*20:16:2,

Itep2v*5:16:2,Pte*1.8:16:2,tyaga/42:16:2,RasTop/16:16:2);

Delay(50); }

end;

if Skor<=0 then Skor:=0.02;

If Skor<=0.0 then

begin

Okno2;

Writeln;

Writeln(' Состав данной массы не может',

' проследовать ');

Writeln(' по данному перегону ');

Readln;

Okno1;

Goto 17;

end;

end

else

begin

Okno2;

Writeln;

Writeln(' Торможение ');

Okno1;

SkorTor:=Skor;

LTMex:=SkorTor;

Repeat

SkorPr:=SkorTor;

UdSopLoc:=23.54+0.108*SkorTor+0.0034*SkorTor*SkorTor+9.81*

RezUkl[i];

OsUdSoVa:=(6.87+(29.43+0.98*SkorTor+0.0245*SkorTor*SkorTor)

/21.4+9.81*RezUkl[i])*Vag4/ObKolVag+(6.87+(78.48+

0.9*SkorTor+0.0245*SkorTor*SkorTor)/21.4+9.81*

RezUkl[i])*Vag6/ObKolVag+(6.87+(58.9+0.373*

SkorTor+0.0206*SkorTor*SkorTor)/21.4+9.81*

RezUkl[i])*Vag8/ObKolVag;

UdSopTor:=(UdSopLoc*LocMas+OsUdSoVa*MasSos)/(LocMas+MasSos);

TorSilaUd:=-(40.85*SkorTor*SkorTor/(DlinOb-Put)+UdSopTor);

TorSila:=TorSilaUd*(LocMas+MasSos);

SkorTor:=SkorTor-0.0034*abs(UdSopTor+abs(TorSilaUd));

If SkorTor<0 then SkorTor:=0; { writeln(TorSilaUd:4:4,' ',UdSopTor:4:4); readln;}

Put:=Put+(SkorPr+SkorTor)*10/(2*36);

Usk:=(SkorTor-SkorPr)/3.6;

Itep2:=0;

ClrScr;

Writeln(#10);

TextColor(Yellow);

Write(' Время движения ');

TextColor(LightGray);

Writeln(Vrem,' c');

TextColor(Yellow);

Write(' Скорость ');

TextColor(LightGray);

Writeln(SkorTor:4:1,' км/ч;');

TextColor(Yellow);

Write(' Путь ');

TextColor(LightGray);

Writeln(Put:8:2,' м;');

TextColor(Yellow);

Write(' Ток генератора тепловоза ');

TextColor(LightGray);

Writeln(Itep2:12:2,' А ');

TextColor(Yellow);

Write(' Расход топлива на движение ');

TextColor(LightGray);

Writeln(RasTop/1000:12:2,' кг');

TextColor(Yellow);

Vrem:=Vrem+1;

Writeln;

Ptep:=100; tyaga:=0;

If kk=0 then

begin

Writeln(FD,Vrem:7,Itep2:15:2,Put:16:2);

Delay(50);

Writeln(FR,Vrem:7,SkorTor*35:16:2,Put*0.1:16:2,UklOb*20:16:2,

Itep2:16:2,Ptep*3.5:16:2,tyaga/42:16:2,RasTop/16:16:2);

Delay(50);

{ end

else

begin

Writeln(FD,Vrem:7,IEl2:15:2,Put:16:2);

Delay(50);

Writeln(FR,Vrem:7,SkorTor*35:16:2,Put*0.1:16:2,UklOb*20:16:2,

Itep2:16:2,Ptep*3.5:16:2,tyaga/42:16:2,RasTop/16:16:2);

Delay(50); }

end;

until (SkorTor<=0) or KeyPressed;

If SkorTor<=0 then

begin

Okno2;

Writeln;

Writeln(' Поезд прибыл на конечную',

' станцию ');

Writeln(' Для продолжения нажмите',

' клавишу ВВОД');

Readln;

Okno1;

Goto 17;

end;

end;

until (Put>DlinOb) or KeyPressed;

17: { kk:=kk+1;

If kk<=1 then begin U:=0; Md:=Mn; Put:=0.0; Skor:=0.02; goto 20; end;}

Close(FD);

Close(FR);

TextMode(LM);

END.