Лабораторные работы / Программы (Pascal) - 1994 / GIDRA / Dina1
.docDina1
program dinam;
const m=2 (*количество решаемых ОДУ*); np=8; nk=5; nv=11;
type areal=array [1..m] of real;
var del (* шаг*),k,x0 (*начальное значение аргумента *),x1 :real;
y0 (*начальные значения функций *),y1 :areal;
n (* число шагов *), kv (* кратность вывода *),i,kl,kl1,kl2 :integer;
vm,v,ak:array[1..nk] of real; p:array[1..np] of real;
hg,h,s,pr:areal; ro,pn,g,x:real;
label 10,20;
FUNCTION SGN(ab:real):integer;
begin if ab>0 then sgn:=1; if ab<0 then sgn:=-1; if ab=0 then sgn:=0; end;
procedure dydx(x:real;y:areal; var pr:areal);
(*расчет производных(pr) для текущих x,y*)
label 300,400; var vm:array[1..5] of real; i:integer;
begin
h:=y; p[7]:=pn * hg[1] / (hg[1] - h[1]); p[8]:=pn*hg[2]/(hg[2]-h[2]);
p[5]:=p[7] + ro * g * h[1] *0.000001; p[6]:=p[8]+ro*g*h[2]*0.000001;
v[1]:=ak[1] * SGN(p[1] - p[5]) * SQRt(ABS(p[1] - p[5]));
v[3]:=ak[3] * SGN(p[5] - p[3]) * SQRt(ABS(p[5] - p[3]));
v[5]:=ak[5]*sgn(p[5]-p[6])*sqrt(abs(p[5]-p[6]));
v[2]:=ak[2] * SGN(p[2] - p[6]) * SQRt(ABS(p[2] - p[6]));
v[4]:=ak[4] * SGN(p[6] - p[4]) * SQRt(ABS(p[6] - p[4]));
pr[1]:=(v[1]-v[3]-v[5])/s[1]; pr[2]:=(v[2] + v[5] - v[4])/s[2];
FOR i:=1 TO 5 do vm[i]:=v[i] * ro;
IF kl = 0 THEN goto 400; IF kl = 1 THEN goto 300;
writeln; writeln(' p(5-7) vm');
writeln(' ', p[5]:nv:5, vm[1]:nv:5);
writeln(' ', p[6]:nv:5, vm[2]:nv:5);
writeln(' ', p[7]:nv:5, vm[3]:nv:5);
writeln(' ', vm[4]:nv:5);
writeln(' ', vm[5]:nv:5);
300: writeln('x = ', x:nv:2, ' y[1]=',y[1]:nv:5,' y[2]=',y[2]:nv:5);
writeln('pr[1]=',pr[1],' pr[2]=',pr[2]);
400:
end;
procedure step(x:real;y:areal; del:real;var x1:real; var y1:areal);
(*шаг интегрирования:по X и Y расчет X1 и y1*)
var pr,y12:areal; i:integer ;
begin
dydx(x,y,pr); for i:=1 to m do y12[i]:=y[i]+del*pr[i]/2;
dydx(x+del/2,y12,pr); x1:=x+del; for i:=1 to m do y1[i]:=y[i]+del*pr[i];
end;
begin
writeln('высоты емкостей (1,2) /м/ = ');read( hg[1], hg[2]);
writeln('площади оснований (1,2) /м^2/ =');read(s[1],s[2]);
writeln('плотность жидкости /кг/м^3/ = ');read(ro);
writeln('начальное давление /МПА/ = ');read(pn);
writeln('давления (1-4) /МПА/ = '); FOR i := 1 TO 4 do read(p[i]);
writeln('коэффициенты сужающих устройств(1-5) = ');
FOR i := 1 TO 5 do read(ak[i]);
10: writeln('н.у: x0,y0[1],y0[2]'); read(x0,y0[1],y0[2]);
20: writeln('число шагов,кратрость вывода ,del'); read(n,kv,del);
writeln('вывод промежуточ.: нет - 0 частич. - 1 полн.- 2 ');read(kl1);
writeln('вывод результатов: частич. - 1 полн.- 2 ');read(kl2);
for i:=1 to n do
begin
kl:=kl1; step(x0,y0,del,x1,y1);
x0:=x1; y0:=y1;
if (i div kv) =(i/kv) then
begin
write('ПРИ i=',i:3);
if kl2=1 then writeln(' x=',x0:10:5,' y[1]=',y0[1]:nv:5,' y[2]=',y0[2]:nv:5);
if kl2=2 then begin kl:=kl2; dydx(x0,y0,pr);end
end
end;
writeln('прекращение счета - 0; продолжение - 1; счет с новыми н.у.- 2 ');
read(kl);
if kl=1 then goto 20; if kl=2 then goto 10;
end.