Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

1528

.pdf
Скачиваний:
8
Добавлен:
07.01.2021
Размер:
1.27 Mб
Скачать

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

3

18

Pa =Po - ( 2+ )Wa2 o /(2 106)

PA =PO-BS**2*RO/2E6

19

PrN=0,098+0,539 10-5 nN

 

 

 

 

 

PR1=0,098+0,539*N1*E5

20

P =1,035 P +(PrN -1,035P )(n/nN)2

PR=1,035*PO+(PR1-1,035*РО*

 

r

 

 

 

 

 

o

 

 

 

 

 

 

o

 

 

 

 

(N/N1)**2

21

T

=370(

 

n

)0,633-10(

 

n

 

)+(T

 

 

-360)

TR=370*(N/N1)**0,333-10*(N/N1)+

r

 

 

 

 

 

 

nN

 

 

 

 

 

 

 

nN

 

 

 

rN

 

+(TR1-360)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

T= TN

 

110 0,125n

 

 

 

 

 

 

 

DТ= DТ1*(110+0,0125*N) / (110+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

110 0,0125nN

 

 

 

 

 

 

 

+ 0,0125*N1)

23

д=( N - 0,18)+0,18(n/nN)

 

 

 

 

 

FD=(FN-0,18)+0,18*(N/N1)

 

 

 

 

 

 

24

r=

To T

 

 

Pr

 

 

 

 

 

 

 

Y=(ТО+DТ)*РR/(TR*(PA*S*FD-PR)

 

 

Pa N Pr

 

 

 

 

 

 

 

 

 

 

Tr

 

 

 

 

 

 

 

 

 

25

Ta=

To T cTr

 

 

 

 

 

 

 

 

 

TA = (TO+DT+V*Y*TR)/(1+Y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

КV =ТО*(PA*S*FD-PR)/((DT+TO)*

 

 

 

 

 

 

 

 

 

1 r

 

 

 

 

 

 

 

 

 

26

V =

 

 

To

 

 

 

 

 

Pa N Pr

 

 

 

 

 

 

 

 

 

 

 

 

*PO*(S-1))

To T

 

 

 

 

 

 

 

 

 

 

Po( 1)

 

 

 

 

 

27

(mc )

=21,475+0,00306(t

c

- 400)

MCC=21.475+0,00306*(ТС1-400)

 

 

 

v

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(mcv) c= M(20,234+0,0031tc)+

МРС=(МСО*(20,234+0,0031*ТС1)+

 

+МС*(29.643+0.01323*ТС1)+МН2*

28

+MCO 2(29,643+0,01323tc)+MH 2

(20.619+0.00063ТС1)+МН*

 

(20,619+0,00063tc)+MH

 

O

 

(24.615+0.0054ТС1)+МN*

 

2

 

(20.129+0.00264ТС1))/М2

 

(24,615+0,0054tc) +MN 2

 

 

 

 

 

(20,129+0,00264tc) /M2

 

 

 

 

 

 

29

(mc )

 

CM= (mc ) + (mc )

/(1+ )

МСМ=(МСС+Y*МРС) / (1+Y)

 

 

v c

 

 

 

 

 

 

 

 

v c

 

r

 

v

c

 

 

r

 

 

(mcv) а= MCO(20,808+0,00055ta+

МРА=(МСО*(20.808+0.00055*

 

+MCO 2(27,545+0,02552ta)+

 

 

 

*ТА1)+МС*(27.545+0.02552*ТА1)+

30

+MH

 

(20,302 +0,00188t )+

 

 

 

+МН2*(20.302+0.00318*ТА1)+

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

+МН*(25.184+0.00242*ТА1)+

 

+MH 2O (25,184+0,00242ta)+

 

+МN*(20.704+0.00029*ТА1))/М2

 

+MN

2

(20,704+0,00029ta) / M2

 

31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(mc )

=20,759+0,0008t

a

 

 

 

 

 

МСА = 20.759+0.0008*ТА

 

 

v

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продолжение табл. 9

1

 

 

2

3

32

(mc )

СМ= (mc )

+ ((mc ) /

МАС=(МСА+Y*МРА)/(1+Y)

 

v a

v a

r v a

 

 

/(1+ r)

 

 

33

K =1+8,315

 

КС=1+8.315*(ТС1-ТА1)/(МСМ*

 

*ТС1-МАС*ТА1)

 

1

 

 

tc ta

 

(mc )

CM t

c

(mc )

 

 

CM t

a

 

 

 

 

 

 

v

c

 

 

 

 

 

v

a

 

 

 

 

34

n1 = K1 -0,02

 

 

 

 

 

 

 

 

 

 

 

 

 

NC=KC-0.02

35

Тс = Та n-1

 

 

 

 

 

 

 

 

 

 

 

 

 

ТС=ТА*S**(NC-1)

36

tc=

(Tс 273) tc

 

100

 

 

 

DТС=((ТС-273)-ТС1)*100/(ТС-273))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tc 273

 

 

 

 

 

 

РС=РА*S**NC

37

Pc = Pa

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

0,4

 

 

 

 

 

 

 

 

 

 

 

 

КСI=(N/N1)**0.4-0.49*(N/N1)+

38

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

+(KCI1-0.51)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

0,49

 

 

 

 

 

 

 

nN

 

 

 

 

 

 

nN

 

 

 

 

 

 

 

 

 

0,51

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ТС1 = ТС-273

39

tc = To - 273

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

U c=8,591 + 0,02299(tc -

400)

UCC=8.591+0.02299*(ТС1-400)

 

 

41

U o= MCO(0,02299tc - 0,205)+

UРC =(MCO*(0.2199*ТС1-0.205)+

+МС*(0.04154*ТС1-2.64)+МН2*

 

+ MCO 2(0,0415tc - 2,64)+МН 2

 

(0.02119ТС1-1.28)+МН*(0.02947*

 

(0,02119tc -1,28)+МН

2

О

 

ТС1-1.078)+MN*(0.02252*ТС1-

 

 

-0.534))/М2

 

(0,02947 tc 1,078)N 2

 

 

 

(0,02252tc 0,534) /М2

 

 

 

 

42

= ( + ) / (1+ )

 

 

 

 

 

 

=(В1+Y)/(1+Y)

 

 

 

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

(Hи

Hи)

 

 

 

 

 

 

 

 

 

U z=[

 

 

 

 

 

 

 

 

 

 

 

Uc rUc

 

 

43

 

 

 

 

L

 

 

 

 

UPZ=(KCI*(HU-DHU)/LT+UCC+

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+Y*UPC)/B*(1+Y)

 

(1 r)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

44

А=(0,02889MСО+0,05275MCO 2+

А=(0.02889*МСО+0.05275*МС+

+0,02805МН

2

+0,04848МН

2

О+

+0.02805*МН2+0.04648*МН+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+0.02885*МN)/M2

 

+0,02805 МN

2

)/M2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продолжение табл. 9

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

3

45

B=(7,162МСО+13,598МСО

2

+

 

 

 

 

B2=(7.162*МСО+13.598*МС+9.961*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

МН2+21.965*МН+6.232*МN)/M2

 

+9,961МН 2+21,965МН 2О+

 

 

 

 

 

 

 

 

 

+6,232МN

2

)/M2

 

 

 

 

 

 

 

 

 

 

 

TZ1=(UPZ+B2)/A

46

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t =(U +B)/A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

47

Tz = tz+273

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TZ=TZ1+273

48

Pz= PcTz / Tc

 

 

 

 

 

 

 

 

 

 

 

PZ=PC*B*TZ/TC

49

=Pz / Pc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LM=PZ/PC

50

Hu = 114(1- )Lo

 

 

 

 

 

 

 

 

 

DHU=114*(1-AF)*LO

 

(mcv) z= MCO(23,196+0,00113tz)+

 

MPZ=(MCO*(23,196+0,00113*TZ1)+

51

+MCO 2(41,652+0,00226tz)+MH 2

 

+MC*(41.652+0.00226*TZ1)+MH2*

(19,995+0,00163tz)+MH 2O

 

 

 

(19.995+0.00163*TZ1)+MH*(28.402+

 

 

 

 

+0,00372*TZ1)+MN*(22.735+0.00113*TZ

 

(28,402+0,00372tz)+MN 2(22 +

 

1))/M2

 

+0,00113tz) / M2

 

 

 

 

 

 

 

 

 

 

 

(mcv) b= МСО(21,468+0,0205 tb)+

МРВ=(МСО*(21.488+0.0205*ТВ1)+

 

+MCO

2

(36,682+0,0049tb)+

 

 

 

 

 

+МС*(36.682+0.0049*ТВ1)+МН2*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.632+0.0018*ТВ1)+МН*(29.975+

52

+MH 2(19,632+0,0018tb) +

 

 

 

 

 

 

 

 

 

 

+0.00549*ТВ1)+МN*(20.922+

 

+MH

2

 

О(24,975+0,00549tb) +

 

 

 

+0.00209*ТВ1))/М2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+MN 2(20,922+0,00209tb) / M2

 

 

 

 

53

K2

=

 

 

 

1+8,315

 

 

 

 

 

 

 

 

 

КР =1+8.315*(TZ1-TB1)/(MPZ*

 

 

 

 

 

 

 

 

 

 

 

tz tb

 

 

 

 

 

 

 

 

 

 

 

*TZ1-MPB*TB1)

 

(mcv)ztz (mcv)btb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

N2 = K2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

NP=KP

55

tb=

(Tb 273) tb

100

 

 

 

 

 

 

DТВ=((ТВ-273)-ТВ1)*100)/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tb 273

 

 

 

 

 

 

 

 

 

 

 

/(ТВ-273)

56

P

=P / n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

РВ = PZ/S**NP

 

b

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

= T

 

 

/3

 

 

 

/ Р

 

 

 

 

 

 

 

 

 

 

 

 

57

b

 

 

Р

 

 

 

 

 

 

 

 

 

 

 

 

TR1=TB/(PB/PR)**0.333

r

 

 

 

 

 

 

 

 

 

 

 

b

 

r

 

 

 

 

 

 

 

 

 

58

Tr=(Tr- T r`) 100 / Tr

 

 

 

 

 

 

DТR=(TR-TR1)*100/TR

59

P =

Рс

[

 

 

 

 

(1

 

1

 

)

 

1

 

 

PI1=PC*(LM*(1-1/S**(NP-1))/

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

1 n

2

1

 

 

 

n2 1

 

 

n 1

 

/(NP-1)-(1-1/S**(NC-1))/(NC-1))/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

(1-

 

 

1

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/(S-1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продолжение табл. 9

1

 

 

 

 

 

 

2

 

 

 

 

3

60

Р = P - (P

- P )

 

 

 

PI=MI*PI1-FI*(PR-PA)

 

i

 

 

i

r

a

 

 

 

 

61

gi = 3600 V /(24,4L Pi)

 

 

 

GI=3600*KV/(24.4*LT*PI)

62

i = 3600/(gi Hu)

 

 

 

KI=3600/(GI*HU)

63

Pм =0,039+0,0132 CN

 

 

 

PM=0.039+0.0132*CN

64

P = P –P

 

 

 

 

 

PE=PI-PM

 

e

i

м

 

 

 

 

 

 

65

м = Pe /Pi

 

 

 

 

 

KM=PE/PI

 

 

 

 

 

 

66

ge = gi / м

 

 

 

 

 

GE=GI/KM

67

 

 

=

 

 

 

 

 

 

KE=KI*KM

 

e

i

м

 

 

 

 

 

 

68

N

i

=N /

 

 

 

 

 

NI1=NE/KM

 

 

e

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

69

D =

 

12 Ni

 

 

 

 

DС=SQRT(12*NI1*T/(3.14*PI*

 

ц

 

 

3,14 Pi n i Sp

 

 

 

N*I*SP))

 

 

 

 

 

 

 

 

 

70

iVh=3,14D2ц Sp i 10 / 4

 

 

VHI = 3.14*DС**2*SP*I/4

71

N

 

= P

iV n / (30 )

 

 

 

NE = PE*VHI*N/(30* )

 

e

e

h

 

 

 

 

 

 

72

Me = 3 104Ne / (3,14 n)

 

 

 

ME=3*10**4*NE/(3.14*N)

73

GТ = ge Ne / 1000

 

 

 

GT =GE*NE/1000

74

Q = GТ Hu

 

 

 

 

 

Q=GT*HU

75

Qe = 3,6 Nе

 

 

 

 

 

QE=3.6*NE

76

qe = Qe 100/Q

 

 

 

 

 

QE1=QE 100/Q

77

Qм = 3,6(Ni-Ne)

 

 

 

 

 

QM=3.6*(NI-NE)

 

 

 

 

 

 

78

qм = Qм 100/Q

 

 

 

 

 

QM1=QM*100/Q

 

 

 

 

 

 

79

tr = Tr 273

 

 

 

 

 

TR1 = TR-273

80

to = To 273

 

 

 

 

 

TO1 = TO-273

 

(mcv) r= MCO(20,211+0,00318tr)+

MPR=(MCO*(20.211+0.00318*TR1)+

 

+MCO 2(31,848+0,0095tr)+MH 2

 

+МС*(31.848+0.0095*TR1)+MH2*

81

(20,33+0,00109t )+MH

 

O

*(20.33+0.00109TRI)+MH*

 

 

 

 

 

 

 

r

2

 

 

*(24.253+0.00603*TR1)+MN*

 

(24,253+0,00603tr) +MN 2

 

*(19.976+0.00293TR1))/M2

 

(19,976+0,00293tr) / M2

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 9

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

3

 

82

(mcp) r=(mcv) r+8,315

 

 

 

 

MRP=MPR+8.315

 

83

(mcv) о=20,759+0,0008to

 

 

 

МОС=20.759+0.0008*ТО1

 

84

(mcp) o=(mcv) a+ 8,315

 

 

 

МОР = МОС+8.315

 

85

Qог= GТ L (mcp) r tr

 

 

 

QG=GT*LT*(B*MRP*TR1-MOP*

 

 

 

 

*TO1)/1000

 

(mcp) o to / 1000

 

 

 

 

 

 

 

 

 

 

 

 

86

q=Q100 / Q

 

 

 

 

 

 

 

QG1 = QG*100/Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m H

 

H

 

 

QW = CW*IC*DC**(1+2*MW)*

 

87

Q = C i D

1+2m

n

и

и

 

*N**MW*(HU-DHU)/(AF*HU)

 

ц

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

Hи

 

 

 

 

 

 

 

 

 

 

 

 

QW1=QW*100/Q

 

88

qw = Qw 100 / Q

 

 

 

 

 

 

 

 

 

 

 

 

 

QH=DHU*GT

 

89

Qн = Hu GТ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

QH1=QH*100/Q

 

90

qн = Qн 100 /Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

91

Qs = Q –(Qe+Qм+Q+Qw+Qн)

 

QS =Q-(QE+QM+QC+QW+QH)

 

 

 

 

92

qs = Qs 100 / Q

 

 

 

 

 

 

 

QS1 = QS 100/Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Начало

Ввод

исходных

данных

Вычисление

Sp, Lo, L o, L o

Задание n(ц)

 

 

 

Вычисление А

 

 

 

 

 

 

 

0

0

N

 

 

 

 

f(n)

 

 

 

 

 

 

 

 

 

 

 

Вычисление L, М1,

mТ, МСО, МСО2,

МН2О, МN2, М2, МН2

Вычисление ,

o, Wвп, Ра,

PrN, Pr

Рис. 4. Схема алгоритма теплового расчета двигателя, работающего на сжиженном газе (начало)

Задание (ц)

Вычисление Tr,

TN, д, r, To, V

Вычисление

(mcv) c, (mcv) c, (mcv) a, (mcv)c,

Вычислениеисление

К1, n1,

(mcvv) cа,,(mcс,v) cc,,

Т t

(mcv)(mca, (vmc)аv),c,

 

0

 

2

tc=tc-5

 

tc

 

tc= tc +5

 

 

 

 

 

Вычисление Рс,

Ни, tc, U с, U c, , А, U z, , , B, tz

Вычисление Tz,

Pz, , (mcv) z, К2, (mcv) b, n2, Тb,

tb

 

0

 

2

tb=tb-5

 

tb

 

tb= tb +5

 

 

 

 

 

Продолжение рис. 4

Вычисление

Pb, T r, ∆Tr

< 0

 

> 2

 

TrN= TrN + 5

 

Tr

 

TrN= TrN -5

 

 

 

 

 

Вычисление

P i, Pi, gi, ηi, Pe, Pм, ηм, ge, ηe

n < nN

 

 

n = nN

 

 

n, IR

 

IR 1

 

IR=1

 

 

 

 

 

Вычисление

NiN, Dц, iVh

Вычисление

Ne, Me, GT, G T, Q, Ni, Qe, qe, Qм, qм, tr, tо, (mcv) r

Продолжение рис. 4

Вычисление

Qн, qн, Qs, qs, (mcv) о, (mcp) r,

Qог, qог, (mcp) о,

Qw, qw

 

< 1

> 3

 

m=m+0,005

qs

 

m=m+0,005

 

 

 

 

Печатать

Pz, ge, ηe, Ne

Печатать

Me,GT,qн,qw

Печатать qе, qм, qог, qs

 

Tmix

Tmax

 

 

 

 

 

 

 

Tо = T- T

 

 

 

Tо

 

 

Tо+ T

 

 

 

nmin

 

 

 

n

 

 

n=n- n

 

 

 

 

 

 

Выход

Окончание рис. 4

4.9. Программа теплового расчета двигателя, работающего на сжиженном газе

REAL MT1,N1T,M1T,MC3,MT2,N2T,M2T,MC4,KT1,KT2,N1,KCI1,MI,NE1,IC, *MW,K,LO1,LO2,LO,LT,MT,M1,MCO,MC,MH2,MH,MN,V2,KV,MCC,MPC,MCM, *MCA,MPA,MAC,KC,NC,KCI,LM,MPZ,MPB,KP,NP,KI,KM,KE,NI1,NE,ME, *MPR,MOC,MOP,N3,N,M2

OPEN(5,FILE='ISX_DAN.DAT')

OPEN(6,FILE='REZYL.TXT') C OPEN(7,FILE='CON')

READ (5, 1)MT1,N1T,M1T,MC3,MT2,N2T,M2T,MC4,AF1,KT1,KT2,PO,TO,S READ (5,3) L,M,I

READ (5,5)RB,CN,N1,FK,BS,TRN,FN,DT1,V,TC1,KCI1,HU,TB1,MI,FI *,NE1,T,IC,CW,MW,N3,K

WRITE (6,2)MT1,N1T,M1T,MC3,MT2,N2T,M2T,MC4,AF1,KT1,KT2,PO,TO,S WRITE (6,4)L,M,I

WRITE (6,2)RB,CN,N1,FK,BS,TRN,FN,DT1,V,TC1,KCI1,HU,TB1,MI,FI *,NE1,T,IC,CW,MW,N3,K

1

FORMAT (7(/),(63X,F8.0))

2

FORMAT (1X, F12.5)

3

FORMAT (63X,I8)

4

FORMAT (1X, I 12)

5

FORMAT (63X,F8.0)

C

ПАРАМЕТРЫ РАБОЧЕГО ТЕЛА

SP=30*CN/N1

LO1=(N1T+M1T/4)*MC3/(0.209*MT1)

LO2=(N2T+M2T/4)*MC4/(0.209*MT2)

LO=LO1+LO2

N=N1

9 A=0.7-N/N1

IF (A.GE.0.0) GO TO 10 AF=AF1

GO TO 11

10AF=AF1-0.78*A**2.36

11CONTINUE LT=AF*LO MT=MC3*MT1+MC4*MT2 M1=LT+1/MT

MCO=0.42*(1-AF)*(LO1/(1+KT1)+LO2/(1+KT2)) MC=N1T*MC3/MT1+N2T*MC4/MT2-MCO

MH2=MCO*(KT1+KT2)/2 MH=M1T*MC3/(2*MT1)+M2T*MC4/(2*MT2)-MH2 MN=0.79*AF*LO

M2=MCO+MC+MH2+MH+MN

B1=M2/M1 WRITE (6,60) M2

60FORMAT (1X,'M2=',E11.4)

C

ПРОЦЕСС ВПУСКА

RO=PO*10**6/(RB*TO)

W=0.05433*SP*N1*FK

12PA=PO-BS*W**2*RO/2E6 PR1=0.098+0.539*N1/10**5

13PR=1.035*PO+(PR1-1.035*PO)*(N/N1)**2 IRAB=0

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]