Лабораторная работа: Расчёт зоны плавления

end

else begin

append(fo);

append(f1);

if l=1 then writeln(f1,Zav);

writeln(f1,P:3:3,' ',Zpl:1:5);

end;

for i:=1 to 5 do begin

X:=W*(1-Ph*Z/(2*Ros*Vsz*H*sqrt(W)));

omega:=Ph*sqrt(X);

sigma:=sqrt(X*2*(LamdaM*(Tb-Tm)+m*sqr(Vj)/2)/(Vbx*Rom*(Cs*(Tm-Ts)+Lamda)));

Vsy:=Vbx*Rom*sigma/(2*Ros*X);

if j=true then

begin

for s:=1 to 6 do begin

Ttv:=(Tm-Ts)*exp(Ros*Vsy*Ros/(LamdaS*Cs)*Y)+Ts;

Tg:=m*sqr(Vj)*YY/(2*LamdaM*sigma)*(1-YY/sigma)+(Tb-Tm)*YY/sigma+Tm;

writeln('Y=',' ',Y:1:2,' ','YY=',YY:1:5,' ','Tg=',' ',Tg:3:3,' ','Ttv=',Ttv:3:1);

append(f);

write(f,YY:1:5,' ',Tg:3:1,' ',Y:1:2,' ',Ttv:3:1);

writeln(f);

YY:=YY+sigma/5;

Y:=Y-0.5;

end;

writeln(f);

close(f);

end;

Y:=0;YY:=0;

k:=k+1;