Материал: 1524
m = m + 0,005
Tо = Tо - T
Вычисление Ar,
(mcv)r, (mcp)r, (mcv)о,
(mcp)о, Qог, qог
Продолжение рис. 2
Вычисление
Qн, qн, Qw, qw, Qs
1 |
4 |
m = m – 0,005 |
g3 |
|
|
|
|
|
Печатать Pz, gi,
e, Ne, Me, GT, ge
Печатать qм, qог, qe, qw, qs
Tо max |
Tо max |
|
Tо |
|
Tо = Tо+ T |
|
|
|
n |
n =n - n |
n
Выход
Окончание рис. 2
2.9. Программа теплового расчета карбюраторного двигателя
REAL N1,MT,KT,NK,MPZ1,MPB1,MI,NE1,IC,MPR1,MW,LO,LT,
*M1,MC,MH2,MH,MN,M2,KV,MCC,MCA,K1,NC,KCI,LM,MP1Z,
*MPZ,MP1B,MPB,KP,KI,KM,KE,NE,ME,MP1R,MPR,MRP,
*MCO,MOP,N3,K,N,NI1,KCI1,NP
OPEN(5,FILE='ISX_DAN.DAT')
OPEN(6,FILE='REZYL.TXT')
C OPEN(7,FILE='CON')
READ (5, 1) C,H,AF1,N1,MT,KT,TO
READ (5, 3) L,M,I
READ (5, 5) PO,RB,CM,FK,BS,TRN,FD1,S,V,TC1,DT1,KCI1, *HU,UPZ1,DUPZ,T1Z,MPZ1,MPB1,TB1,MI,
*FI,NE1,T,IC,MPR1,CW,MW,N3,K
WRITE(6, 2) C,H,AF1,N1,MT,KT,TO
WRITE(6, 4) L,M,I
WRITE(6, 2) PO,RB,CM,FK,BS,TRN,FD1,S,V,TC1,DT1,KCI1, *HU,UPZ1,DUPZ,T1Z,MPZ1,MPB1,TB1,MI, *FI,NE1,T,IC,MPR1,CW,MW,N3,K
1FORMAT (7 (/),(63X,F8.0))
2FORMAT (1X, F12.5)
3FORMAT (63X,I8)
4FORMAT (1X, I12)
5FORMAT (63X,F8.0)
C |
ПАРАМЕТРЫ РАБОЧЕГО ТЕЛА |
SP=30*CM/N1
LO=(C+3*H)/(12*0.209)
N=N1
18A=0.7-N/N1
19IF (A.GE.0) GO TO 20
AF=AF1
GO TO 21
20AF=AF1-0.78*A**2.36
21CONTINUE LT=AF*LO
M1=LT+1/MT
MCO=0.418*LO*(1-AF)/(1+KT)
MC=C/12-MCO
|
MH2=KT*MCO |
|
MH=H/2-MH2 |
|
MN=0.791*AF*LO |
|
M2=MCO+MC+MH2+MH+MN |
|
B1=M2/M1 |
C |
ПРОЦЕСС ВПУСКА |
|
IRAB=O |
|
DO 50 JJ=1, M |
|
IRAB= IRAB+1 |
C-------------------------- |
|
|
WRITE(6,100)JJ |
100 FORMAT(1X,'JJ=',I5)
C--------------------------
|
J=L+(JJ-1)*I |
|
TO=J |
|
RO=PO*10**6/(RB*TO) |
|
W=0.05433*SP*N*FK |
|
PA=PO-BS*W**2*RO/2E6 |
22 PR1=0.098+0.539*N/10**5 |
|
|
PR=1.035*PO+(PR1-1.035*PO)*(N/N1)**2 |
C |
IRAB1=O |
24 TR=370*(N/N1)**0.333-10*(N/N1)+(TRN-360) |
|
C-------------------------- |
|
C |
IRAB1=IRAB1+1 |
C |
WRITE(6,103)TRN |
C 103 FORMAT(1X,'TRN=',E11.4) |
|
C |
IF(IRAB1.GT.500) STOP |
C-------------------------- |
|
|
DT=DT1*(110+0.0125*N)/(110+0.0125*N1) |
|
FD=(FD1-0.18)+0.18*(N/N1) |
|
Y=(TO+DT)*PR/(TR*(PA*S*FD-PR)) |
|
TA=(TO+DT+V*Y*TR)/(1+Y) |
|
KV=TO*(PA*S*FD-PR)/((TO+DT)*PO*(S-1)) |
C |
ПРОЦЕСС СЖАТИЯ |
C |
IRAB2=O |
25 MCC=21.475+0.00306*(TC1-400) |
|
C-------------------------- |
|
C |
IRAB2=IRAB2+1 |
C |
WRITE(6,101)TC1 |
C 101 FORMAT(1X,'TC1=',E11.4) |
|
C |
IF(IRAB2.GT.500) STOP |
C--------------------------
TA1=TA-273
MCA=20.758+0.0008*TA1
K1=1+8.315*(TC1-TA1)/(MCC*TC1-MCA*TA1)
NC=K1-0.02
TC=TA*S**(NC-1)
DTC=((TC-273)-TC1)*100/(TC-273)
IF (DTC.GT.2.0) GO TO 26
IF (DTC.LT.0.0) GO TO 27
GO TO 28
26TC1= TC1+5 GO TO 25
27TC1= TC1-5
GO TO 25
28 CONTINUE PC=PA*S**NC
C |
ПРОЦЕСС СГОРАНИЯ |
|
DHU=114*(1-AF)*LO |
|
UCC=8.591+0.02299*(TC1-400) |
|
UP1C=9.1123+0.02459*(TC1-400) |
|
AC=UP1C/9.1123 |
|
UPC=UP1C+AC**1.25*(1.69*(AF-0.7)-0.9565*(AF-0.7)**1.65) |
|
B=(B1+Y)/(1+Y) |
|
KCI=(N/N1)**0.4-0.49*(N/N1)+(KCI1-0.51) |
|
UPZ=(KCI*(HU-DHU)/LT+UCC+Y*UPC)/(B*(1+Y)) |
|
BZ=16*(AF-0.7)-9.075*(AF-0.7)**1.65 |
|
TZ1=(UPZ/(1+BZ/UPZ1)-UPZ1)/DUPZ+T1Z |
|
TZ= TZ1+273 |
|
PZ=B*PC*TZ/TC |
|
LM=PZ/PC |
C |
ПРОЦЕСС РАСШИРЕНИЯ |
|
MP1Z=MPZ1+0.003172*(TZ1-T1Z) |
|
AZ=MP1Z/MPZ1 |
|
MPZ=MP1Z+AZ**1.1*(6.75*(AF-0.7)-3.95*(AF-0.7)**1.75) |