Материал: 1524
<0 >2
<0 >2
Продолжение рис. 5
Вычисление
P i, Pi, ηi, Pм, Pе, ηм, ηе, gе
n ≠ nN |
n = nN |
|
n, IR |
IR≠1 |
IR=1 |
|
Вычисление |
|
NiN, Dц, iVh |
Вычисление
Nе, Mе, Ni, GT, Qе, qе, Qм, qм
Вычисление tr, to, (mcv) r
Вычисление
(mcp) 1r, (mcp) r, (mcp) o
Вычисление
Qог, qог, Qw, qw, Qs, qs
|
|
Продолжение рис. 5 |
|
||
|
|
< 2 |
>5 |
|
|
|
|
|
|
||
|
|
qs |
|
|
|
|
m = m-0,005 |
|
m = m+0,005 |
|
|
|
|
|
|||
|
|
|
|
|
|
Печатать
Pz, gе, ηе, Nе
Печатать
Mе, GT, qе
Печатать
qм, qог, qw,
|
|
>ξmin |
< ξmax |
|
|
|
|
||
|
ξ =ξ-∆ ξ |
ξ |
|
ξ =ξ+∆ ξ |
|
|
|||
|
|
|
|
|
nmin
А |
n = n+∆n |
|
Выход
Окончание рис. 5
5.9. Программа теплового расчета дизельного двигателя с наддувом
C
REAL N1,MT,NK,KCI1,LM,MI,NE1,IC,MW,LO,LT,M1,MC,MH,MN,MO,M2,
*KV,MCC,MCA,KC,NC,KCI,MPZ1,MCZ,MPZ,MPB1,MCB,MPB,KP,NP,KI,KM,
*KE,NI1,NE,ME,NI,MCR,MPR1,MPR,MPO,K,N3,N OPEN(5,FILE='ISX_DAN.DAT')
OPEN(6,FILE='REZYL.TXT')
C OPEN(7,FILE='CON')
READ (5,1)
C,H,O,AF1,N1,MT,TO,PK1,PO,NK,RB,CN,FK,BS,TRN,DT1,S,V,TC1 READ (5,3) L,M,I
READ (5,5) KCI1,HU,LM,TB1,MI,FI,NE1,T,IC,CW,MW,K,N3
WRITE (6,2) C,H,O,AF1,N1,MT,TO,PK1,PO,NK,RB,CN,FK,BS, *TRN,DT1,S,V,TC1
WRITE (6,4) L,M,I
WRITE (6,2) KCI1,HU,LM,TB1,MI,FI,NE1,T,IC,CW,MW,K,N3
1FORMAT (7(/),(63X,F8.0))
2FORMAT (1X,F12.5)
3FORMAT (63X,I8)
4FORMAT (1X,I12)
5FORMAT (63X,F8.0)
C |
ПАРАМЕТРЫ РАБОЧЕГО ТЕЛА |
LO=(C+3*(H-O/8))/(12*0.209) N=N1
10 AF= (AF1-0.6)+(N/N1-0.3)**0.75-0.2357*(N/N1-0.3) LT=AF*LO
M1=LT+1/MT
MC=C/12.0
MH=H/2.0
MN=0.79*AF*LO MO=0.209*(AF-1)*LO M2=MC+MH+MN+MO
B1=M2/M1
WRITE(6,50)M2
50 FORMAT('M2=',E11.4)
C |
ПРОЦЕСС ВПУСКА |
TK1=TO*(PK1/PO)**((NK-1)/NK)
T1K=TK1-273
TK=(T1K-56)+112*(N/N1-0.5)+273
PK=(PK1-0.056)+0.112*(N/N1-0.5)
RK=PK*10**6/(RB*TK)
SP=30*CN/N1
W=0.05433*SP*N*FK
11PA=PK-BS*RK*W**2/2E6 PR1=0.9*PK
PR=(PR1-0.056)+0.07*(N/N1-0.2)
12TR=(TRN+45)-90*(N/N1-0.5)+273
DT=DT1*(110+0.0125*N)/(110+0.0125*N1)
13Y=(TK+DT)*PR/(TR*(PA*S-PR))
14TA=(TK+DT+V*Y*TR)/(1+Y)
15KV=TK*(PA*S-PR)/((TK+DT)*PK*(S-1))
WRITE(6,51)KV
51 FORMAT('KV=',E11.4)
C |
ПРОЦЕСС СЖАТИЯ |
16MCC=22.408+0.00299*(TC1-700)
TA1=TA-273
MCA=20.838+0.00146*(TA1-100) KC=1+8.315*(TC1-TA1)/(MCC*TC1-MCA*TA1)
NC=KC
17TC=TA*S**(NC-1)
DTC=((TC-273)-TC1)*100/(TC-273)
IF(DTC.GT.2.0) GO TO 18
IF(DTC.LT.0.0) GO TO 19
GO TO 20
18TC1=TC1+5
GO TO 16
19TC1=TC1-5
GO TO 16
20CONTINUE PC=PA*S**NC
WRITE(6,52)PC
52 FORMAT ('PC=',E11.4)
C |
ПРОЦЕСС СГОРАНИЯ |
B=(B1+Y)/(1+Y)
IRAB=0
DO 46 JJ=1,M
IRAB= IRAB+1
J=L+(JJ-1)*I
KCI1=J/1000.0
IF(N/N1.GE.1.0) GO TO 21
GO TO 22
21KCI=KCI1 GO TO 60
22KCI=(KCI1-0.2)+(N/N1-0.3)**0.75-0.807*(N/N1-0.3)
60CONTINUE