1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
| unit Unit1;
interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, Math;
type
TForm1 = class(TForm)
Button1: TButton;
Edit1: TEdit;
Edit2: TEdit;
Edit3: TEdit;
Edit4: TEdit;
Edit5: TEdit;
Label1: TLabel;
Label2: TLabel;
Label3: TLabel;
Label4: TLabel;
Label5: TLabel;
Button2: TButton;
Label6: TLabel;
procedure calculate(Sender: TObject);
procedure Button2Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;
var
Form1: TForm1;
implementation
{$R *.dfm}
procedure TForm1.calculate(Sender: TObject);
var
A,B,C,D,E: real;
X1,X2: real;
S1,S2,S3,S4: real;
Q1,Q2,Q3,Q4: real;
W1,W2,W3,W4,W5: real;
E1,E2,E3,E4,E5: real;
R1,R2,R3,R4,T1,T2,T3,T4: real;
Y1,Y2,Y3,Y4: real;
U1,U2,U3,U4: real;
I1,I2,I3,I4: real;
O1,O2,O3,O4: real;
P1,P2,P3,P4: real;
G1,G2,G3,G4,F1,F2,F3,F4: real;
begin
DecimalSeparator:= '.';
// Ввод исходных данных.
A:= StrToFloat(Edit1.Text);
B:= StrToFloat(Edit2.Text);
C:= StrToFloat(Edit3.Text);
D:= StrToFloat(Edit4.Text);
E:= StrToFloat(Edit5.Text);
// Вычисление числа ступеней и степени сжатия.
X1:=B/A;
X2:=(ln(X1))/(ln(4));
S1:=exp(3*ln(X1/1));
S2:=exp(3*ln(X1/0.95));
S3:=exp(3*ln(X1/0.9025));
S4:=exp(3*ln(X1/0.857375));
// Определение давления газа по ступеням.
Q1:=A*S1;
Q2:=Q1*S2;
Q3:=Q2*S3;
Q4:=Q3*S4;
// Определение потерь давления между ступенями. дальше идет большой ошибкен
W1:=0.3*(2.6/power(0.26,A*1000000));
W2:=0.3*(2.6/power(0.26,Q1*1000000));
W3:=0.3*(2.6/power(0.26,Q2*1000000));
W4:=0.3*(2.6/power(0.26,Q3*1000000));
W5:=0.3*(2.6/power(0.26,Q4*1000000));
E1:=0.7*(2.6/power(0.26,Q1*1000000));
E2:=0.7*(2.6/power(0.26,Q2*1000000));
E3:=0.7*(2.6/power(0.26,Q3*1000000));
E4:=0.7*(2.6/power(0.26,Q4*1000000));
// Определение давлений по ступеням с учетом потерь.
R1:=(1-W1)*A;
R2:=(1-W2)*Q1;
R3:=(1-W3)*Q2;
R4:=(1-W4)*Q3;
T1:=(1-E1)*Q1;
T2:=(1-E2)*Q2;
T3:=(1-E3)*Q3;
T4:=(1-E4)*Q4;
// Определение коэффициента подачи.
// Определение объемного коэффициента.
Y1:=1-0.08*(power((1/1.2),(T1/R1))-1);
Y2:=1-0.08*(power((1/1.3),(T2/R2))-1);
Y3:=1-0.08*(power((1/1.35),(T3/R3))-1);
Y4:=1-0.08*(power((1/1.35),(T4/R4))-1);
// Определение коэффициента дросселирования.
U1:=1-(((1-0.08)/(23*Y1))*((A-R1)/A));
U2:=1-(((1-0.08)/(23*Y2))*((Q1-R2)/Q1));
U3:=1-(((1-0.08)/(23*Y3))*((Q2-R3)/Q2));
U4:=1-(((1-0.08)/(23*Y4))*((Q3-R4)/Q3));
// Определение коэффициента подогрева.
I1:=1-0.01*((T1/R1)-1);
I2:=1-0.01*((T2/R2)-1);
I3:=1-0.01*((T3/R3)-1);
I4:=1-0.01*((T4/R4)-1);
// Определение коэффициентов плотности для компрессоров в хорошем состоянии.
O1:=0.98;
O2:=0.975;
O3:=0.965;
O4:=0.965;
// Определение коэффициента подачи.
P1:=Y1*U1*I1*O1*1;
P2:=Y2*U2*I2*O2*1;
P3:=Y3*U3*I3*O3*1;
P4:=Y4*U4*I4*O4*1;
// Определение температура сжимаемого газа;
G1:=C;
G2:=C+5;
G3:=G2+5;
G4:=G3+5;
F1:=G1*power((0.32/1.32),(T1/R1));
F2:=G2*power((0.34/1.34),(T2/R2));
F3:=G3*power((0.36/1.36),(T3/R3));
F4:=G4*power((0.38/1.38),(T4/R4));
showmessage('F4');
end;
procedure TForm1.Button2Click(Sender: TObject);
begin
Close;
end;
end. |
(
