Проведение аналитических исследований на ПЭВМ с использованием прикладной программы

(3.17)

Формулы прогонки включают в себя коэффициенты прогонки

(3.18)

  (3.19)

Для i- ых узлов

(3.20)

(3.21)

При i=N условия 3- го рода

(3.22)

Недостатком рассмотренной схемы является рост погрешности расчета при увеличении шага по времени.

4. Разработка алгоритма, блок-схемы решения задачи и прикладной программы

5. Разработка и отладка прикладной программы для ПЭВМ

unit Unit1;

interface

uses

Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,

Dialogs, StdCtrls,Math;

const

dtau=5;

dx=0.01;

Tpmin= 1000;

Tpmax=1280;

tau1=36000;

zm=4;

type

TForm1 = class(TForm)

Button3: TButton;

Memo3: TMemo;

Memo1: TMemo;

Button1: TButton;

Memo2: TMemo;

Button2: TButton;

procedure Button3Click(Sender: TObject);

procedure Button1Click(Sender: TObject);

procedure Button2Click(Sender: TObject);

private

{ Private declarations }

public

{ Public declarations }

end;

var

Form1: TForm1;

Tn,ff,b,tau,Tsr,razn,a,z,Tp,al,Tpe:real;

T: array [1..19] of real= (20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20);

T1: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

alfa: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

betta: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

c: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

ro: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

l: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

f: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

lamplys: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

fplys: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

lammin: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

fmin: array [1..19] of real= (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0);

i: integer;

implementation

{$R *.dfm}

procedure TForm1.Button3Click(Sender: TObject);

begin

tau:=0;

Tsr:=20;

z:=0;

while z<zm do

begin

repeat

if tau<3600 then Tp:= Tpmin else

Tp:=Tpmin+(Tpmax-Tpmin) *tau/tau1;

l[1]:=44.2-0.016*T[1];

l[2]:=44.2-0.016*T[2];

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,

l[19]:=44.2-0.016*T[19];

lamplys[1]:=-0.016*0.5*(T1[1]+T1[2])+44.2;

fplys[1]:= lamplys[1]*dtau/(c[1]*ro[1]*dx*dx);

lamplys[2]:=-0.016*0.5*(T1[2]+T1[3])+44.2;

lammin[2]:= -0.016*0.5*(T1[2]-T1[1])+44.2;

fmin[2]:= lammin[2]*dtau/(c[2]*ro[2]*dx*dx);

fplys[2]:= lamplys[2]*dtau/(c[2]*ro[2]*dx*dx);

c[1]:=0.198*t1[1]+463.6;

ro[1]:= -0.308*T1[1]+7836;

lamplys[1]:=-0.016*0.5*(T1[1]+T1[2])+44.2;

fplys[1]:= lamplys[1]*dtau/(c[1]*ro[1]*dx*dx);

alfa[1]:= 2*fplys[1]/(1+2*fplys[1]);

beta[1]:= Tsr/(1+2*fplys[1]);

c[2]:=0.198*T1[2]+463.6;

ro[2]:= -0.308*T1[2]+7836;

lamplys[2]:=-0.016*0.5*(T1[2]+T1[3])+44.2;

lammin[2]:= -0.016*0.5*(T1[2]-T1[1])+44.2;

fmin[2]:= lammin[2]*dtau/(c[2]*ro[2]*dx*dx);

fplys[2]:= lamplys[2]*dtau/(c[2]*ro[2]*dx*dx);

alfa[2]:= fplys[2]/(1+fplys[2]+fmin[2]*(1-alfa[1]));

beta[2]:= (Tsr+fmin[2]*beta[1])/(1+fplys[2]+fmin[2]*(1-alfa[1]));

c[3]:=0.198*T1[3]+463.6;

ro[3]:= -0.308*T1[3]+7836;

lamplys[3]:=-0.016*0.5*(T1[3]+T1[4])+44.2;

lammin[3]:= -0.016*0.5*(T1[3]-T1[2])+44.2;

fmin[3]:= lammin[2]*dtau/(c[3]*ro[3]*dx*dx);

fplys[3]:= lamplys[3]*dtau/(c[3]*ro[3]*dx*dx);

alfa[3]:= fplys[3]/(1+fplys[3]+fmin[3]*(1-alfa[2]));

beta[3]:= (Tsr+fmin[3]*beta[3])/(1+fplys[3]+fmin[3]*(1-alfa[2]));

c[4]:=0.198*T1[4]+463.6;

ro[4]:= -0.308*T1[4]+7836;

lamplys[4]:=-0.016*0.5*(T1[4]+T1[5])+44.2;

lammin[4]:= -0.016*0.5*(T1[4]-T1[3])+44.2;

fmin[4]:= lammin[4]*dtau/(c[4]*ro[4]*dx*dx);

fplys[4]:= lamplys[4]*dtau/(c[4]*ro[4]*dx*dx);

alfa[4]:= fplys[4]/(1+fplys[4]+fmin[4]*(1-alfa[3]));

beta[4]:= (Tsr+fmin[4]*beta[4])/(1+fplys[4]+fmin[4]*(1-alfa[3]));

c[5]:=0.198*T1[5]+463.6;

ro[5]:= -0.308*T1[5]+7836;

lamplys[5]:=-0.016*0.5*(T1[5]+T1[6])+44.2;

lammin[5]:= -0.016*0.5*(T1[5]-T1[4])+