1FORMAT (7(/),(63X,F8.0))
2FORMAT (1X, F12.5)
3FORMAT (63X,I8)
4FORMAT (1X, I12)
5FORMAT (63X,F8.0)
C |
ПАРАМЕТРЫ РАБОЧЕГО ТЕЛА |
SP=30*CN/NN
HU1=35.8*CH4+63.4*C2H6+91*C3H8
LO1=(N1+M1/4)*CH4/0.21
LO2=C2H6*(N2+M2/4)/0.21
LO3=C3H8*(N4+M3/4)/0.21
LO=LO1+LO2+LO3
N=NN
9 A=0.7-N/NN
IF (A.GE.0.0) GO TO 10 AF=AF1
GO TO 11
10 AF=AF1-0.78*(0.7-N/NN)**2.36
11 CONTINUE
C WRITE (6, 60) M2,AF LT=AF*LO
MI1=1+LT K1=(1.515*HC1-0.1727)**0.3225 K2=(1.515*HC2-0.1727)**0.3225 K3=(1.515*HC3-0.1727)**0.3225
MCO=0.42*(1-AF)*(LO1/(1+K1)+LO2/(1+K2)+LO3/(1+K3)) MH2=MCO*(CH4*K1+C2H6*K2+C3H8*K3) MH=(M1*CH4+M2*C2H6+M3*C3H8)/2-MH2 MC=N1*CH4+N2*C2H6+N4*C3H8-MCO MN=0.79*AF*LO
M2=MCO+MC+MH2+MH+MN
B1=M2/MI1
CWRITE (6, 60) M2,MI1,B2,MCO,MC,MH2,MH,MN WRITE (6, 60) M2
60 FORMAT (' M2= ', (E11.4))
C |
ПРОЦЕСС ВПУСКА |
RO=PO*10**6/(RB*TO)
SP=30*CN/NN
W=0.054*SP*NN*FK
12PA=PO-BS*RO*W**2/(2*10**6) PR1=0.098+0.539*NN/10**5
13PR=1.035*PO+(PR1-1.035*PO)*(N/NN)**2 IRAB=0
C WRITE (6, 160) RO,SP,W,PA,PR1,PR
C160 FORMAT (' RO= ', (E11.4)) DO 55 JJ=1, M
IRAB= IRAB+1 J=L+(JJ-1)*I S=J/10.0 IRAB1=0
14 TR=370*(N/NN)**0.333-10*(N/NN)+(TRN-360)
C-----------------------------------
C IRAB1=IRAB1+1
C WRITE (6,553) TRN
C 553 FORMAT(' TRN=',E11.4)
C IF (IRAB1.GT.500) STOP
C-----------------------------------
|
DT=DT1*(110+0.0125*N)/(110+0.0125*NN) |
|
|
FN=(FN1-0.18)+0.18*(N/NN) |
|
|
Y=(TO+DT)*PR/(TR*(PA*S-PR)) |
|
|
TA=(TO+DT+V*Y*TR)/(1+Y) |
|
C |
KV=(TO*(PA*S*FN-PR))/((TO+DT)*PO*(S-1)) |
|
WRITE (6,61) Y,KV,TO,DT,PR,TR,PA,S,FN |
||
C |
61 FORMAT (' Y=', E11.4,' KV=',(E11.4)) |
|
C |
STOP |
|
C |
ПРОЦЕСС СЖАТИЯ |
|
|
IRAB2=0 |
|
|
MCC=21.475+0.00306*(TC1-400) |
|
|
DC1=TC1-400 |
|
|
18 MPC=(MCO*(21.475+0.0031*DC1)+MC*(29.643+0.0132*DC1)+ |
|
|
*MH2*(20.619+0.00063*DC1)+MH*(24.615+0.0054*DC1)+ |
|
C |
*MN*(21.186+0.00264*DC1))/M2 |
|
IRAB2=IRAB2+1 |
||
C |
||
C |
WRITE (6,554) TC1 |
|
C |
554 FORMAT(/,' TC1=',E11.4) |
|
C |
IF (IRAB2.GT.500) STOP |
|
C----------------------------------- |
MCM=(MCC+Y*MPC)/(1+Y) |
|
|
||
|
TA1=TA-273 |
|
|
MCA=20.759+0.0008*TA1 |
|
|
MPA=(MCO*(20.808+0.00055*TA1)+MC*(27.545+0.0255*TA1)+ |
|
|
*MH2*(20.302+0.00318*TA1)+MH*(25.184+0.00242*TA1)+ |
|
|
*MN*(20.704+0.00029*TA1))/M2 |
|
|
MAC=(MCA+Y*MPA)/(1+Y) |
|
|
KC=1+8.315*(TC1-TA1)/(MCM*TC1-MAC*TA1) |
|
|
NC=KC-0.01 |
|
|
TC=TA*S**(NC-1) |
|
C |
DTC=((TC-273)-TC1)*100/(TC-273) |
|
WRITE (6,651) DTC,MPC,MCM,MCA,MAC,MPA |
||
C |
||
C |
651 FORMAT(' DTC=',E10.4,' MPC=',E10.4,' MCM=',E10.4,' MCA=',E10.4, |
|
C |
*' MAC=',E10.4,/,' MPA=',E10.4) |
|
C |
WRITE (6,751) KC,TA,S,TC |
|
C |
751 FORMAT(' KC=',E10.4,' TA=',E10.4,' S=',E10.4,' TC=',E10.4) |
|
C |
WRITE (6,752) MCO,MC,MH2,MH,MN |
|
C |
752 FORMAT(' MCO=',E10.4,' MC=',E10.4,' MH2=',E10.4,' MH=',E10.4, |
|
C |
*' MN=',E10.4) |
|
C----------------------------------- |
IF (DTC.GT.2.0)GO TO 20 |
|
|
||
|
IF (DTC.LT.0.0)GO TO 21 |
|
|
GO TO 22 |
20TC1=TC1+5 GO TO 18
21TC1=TC1-5 GO TO 18
22 CONTINUE PC=PA*S**NC WRITE (6, 62) PC
62 FORMAT (' PC=',E11.4)
C |
ПРОЦЕСС СГОРАНИЯ |
DHU=112.5*(1-AF)*LO
TC1=TC-273
KCI=(N/NN)**0.4-0.49*(N/NN)+(KCI1-0.51) UCC=8.591+0.02299*(TC1-400) UPC=(MCO*(0.02199*TC1-0.205)+MC*(0.0415*TC1-2.641)+MH2* *(0.02119*TC1-1.28)+MH*(0.029*TC1-1.078)+MN*(0.0225*TC1-0.534))/M2 B=(B1+Y)/(1+Y)
UPZ=(KCI*(HU-DHU)/MI1+UCC+Y*UPC)/(B*(1+Y)) A=(0.02889*MCO+0.05275*MC+0.02805*MH2+0.04648* *MH+0.02805*MN)/M2 B2=(7.162*MCO+13.598*MC+9.961*MH2+21.995*MH+6.232*MN)/M2 TZ1=(UPZ+B2)/A
TZ=TZ1+273
PZ=B*PC*TZ/TC
LM=PZ/PC
CWRITE (6, 63)TZ,UPZ,HU,DHU,MI1,UCC,Y,UPC,B2,A,B,KCI WRITE (6, 63)TZ
C |
63 FORMAT (' TZ=', (E11.4)) |
|
ПРОЦЕСС РАСШИРЕНИЯ |
||
|
IRAB3=0 |
|
|
MPZ=(MCO*(23.196+0.00113*TZ1)+MC*(41.652+0.00226*TZ1)+MN* |
|
|
*(22.735+0.00113*TZ1)+MH2*(19.995+0.00163*TZ1)+MH* |
|
|
*(28.402+0.00372*TZ1))/M2 |
|
|
30 MPB=(MCO*(21.468+0.00205*TB1)+MC*(36.682+0.0049*TB1)+MH2* |
|
|
*(19.632+0.0018*TB1)+MH*(24.972+0.00549*TB1)+MN*(20.992+0.00209* |
|
C |
*TB1))/M2 |
|
IRAB3=IRAB3+1 |
||
C |
||
C |
WRITE (6,555) TB1 |
|
C |
555 FORMAT(/,' TB1=',E11.4) |
|
C |
IF (IRAB3.GT.500) STOP |
|
C----------------------------------- |
KP=1+8.315*(TZ1-TB1)/(MPZ*TZ1-MPB*TB1) |
|
|
||
|
NP=KP-0.01 |
|
|
TB=TZ/S**(NP-1) |
|
C |
DTB=((TB-273)-TB1)*100/(TB-273) |
|
WRITE (6,755) DTB,TZ1,MPZ,MPB |
||
C |
||
C |
755 FORMAT(' DTB=',E11.4,' TZ1=',E11.4,' MPZ=',E11.4,' MPB=',E11.4) |
|
C |
WRITE (6,655) TZ,S,KP,TB |
|
C |
655 FORMAT(' TZ=',E11.4,' S=',E11.4,' KP=',E11.4,' TB=',E11.4) |
|
C----------------------------------- |
IF (DTB.GT.2.0) GO TO 31 |
|
|
||
|
IF (DTB.LT.0.0) GO TO 32 |
|
|
GO TO 33 |
|
|
31 TB1=TB1+5 |
|
|
GO TO 30 |
|
|
32 TB1=TB1-5 |
|
|
GO TO 30 |
|
|
33 CONTINUE |
|
|
PB=PZ/S**NP |
|
|
WRITE (6,64) PB |
|
|
64 FORMAT (' PB=',E11.4) |
|
|
TRI=TB/(PB/PR)**0.333 |
|
|
DTR=(TR-TRI)*100/TR |
|
|
IF(DTR.GT.2.0) GO TO 34 |
|
|
IF(DTR.LT.0.0) GO TO 35 |
|
|
GO TO 36 |
|
|
34 TRN=TRN-5 |
|
GO TO 14 |
|
35 TRN=TRN+5 |
|
GO TO 14 |
C |
36 CONTINUE |
ПОКАЗАТЕЛИ РАБОТЫ ЦИКЛА |
|
|
PI1=PC*(LM*(1-1/S**(NP-1))/(NP-1)-(1-1/S**(NC-1))/(NC-1))/(S-1) |
|
PI=MI*PI1-FI*(PR-PA) |
|
GI=3.6*KV/(MI1*PI) |
|
KI=3.6/(HU1*GI) |
|
PM=0.039+0.0132*CN |
|
PE=PI-PM |
|
KM=PE/PI |
|
37 GE=GI/KM |
|
38 KE=KI*KM |
|
IF(N.EQ.NN.AND.IRAB.EQ.1) GO TO 39 |
|
GO TO 40 |
|
39 NI1=NE1/KM |
C |
DC=SQRT(1.2*NI1*T/(3.14*PI*NN*SP)) |
WRITE (6,66) KE,DC,NI1,T,PI,NN,SP |
|
|
VHI=3.14*SP*10*IC*DC**2/4 |
|
40 CONTINUE |
|
NE=PE*VHI*N/(30*T) |
|
ME=3*10**4*NE/(3.14*N) |
|
GT=GE*NE1 |
|
GT1=GE/22.5 |
|
NI=NE/KM |
CWRITE (6,66) KE,KI,KM,PI,KV,MI1,GI,GT,DC WRITE (6,66) KE
66 FORMAT (' KE=',(E11.4))
C |
ТЕПЛОВОЙ БАЛАНС ДВИГАТЕЛЯ |
Q=GT*HU1
QE=3.6*NE
QE1=QE*100/Q
QM=3.6*(NI-NE)
QM1=QM*100/Q
TR1=TR-273
TO1=TO-273
MPR= (MCO*(20.211+0.00318*TR1)+MC*(31.848+0.0095*TR1)+ *MH2*(20.33+0.00109*TR1)+MH*(24.253+0.00603*TR1)+ *MN*(19.976+0.00293*TR1))/M2
MRP=MPR+8.315
MOC=20.759+0.0008*TO1
MOP=MOC+8.315
QG=(GT/22.5)*(B*MRP*TR1-MOP*TO1)*LT/1000 QG1=QG*100/Q
QH=DHU*GT/22.5
QH1=QH*100/Q
C WRITE (6,67) Q,GT,HU1,QE,QM,QG,QH,MRP,MOP,TR1,TO1,LT
C67 FORMAT (' Q=',(E11.4)) IRAB4=0
41 QW=CW*IC*(DC**(1+2*MW))*(N**MW)*(HU-DHU)/(AF*HU)
C----------------------------------- |
IRAB4=IRAB4+1 |
C |
|
C |
WRITE (6,556) MW,CW,IC,DC,N,HU,DHU,AF,HU |
C |
WRITE (6,556) MW |
C |
556 FORMAT(' MW=',(E11.4)) |
C |
IF (IRAB4.GT.500) STOP |
C-----------------------------------
QW1=QW*100/Q
QS=Q-(QE+QM+QG+QH+QW)
QS1=QS*100/Q
C--------------------------
C WRITE (6,756) QS1,QW,QS
C 756 FORMAT(' QS1=',E11.4,' QW=',E11.4,' QS=',E11.4)
C----------------------------
IF (QS1.GT.3.0) GO TO 42
IF (QS1.LT.1.0) GO TO 43
GO TO 44
42MW=MW+0.005 GO TO 41
43MW=MW-0.005 GO TO 41
44CONTINUE WRITE(6,45)
45FORMAT ( 1X,30('*')) WRITE (6,70) S,N
70FORMAT (' S=', E11.4,' N=', E11.4) WRITE (6, 46) PZ, GE, KE
WRITE (6, 47) NE, ME, GT WRITE (6, 48) QE1, QM1, QG1 WRITE (6, 49) QH1, QW1, QS1
46FORMAT (3X,'PZ=' ,F5.3,2X,'GE=',F5.1,2X,'KE=' ,F6.4)
47FORMAT (3X,'NE= ',F6.1,2X,'ME=',F5.1,2X,'GT=' ,F5.2)
48FORMAT (3X,'QE1=',F6.3,2X,'QM=',F5.3,2X,'QG1=',F6.3)
49FORMAT (3X,'QH1=',F5.2,2X,'QW=',F5.2,2X,'QS1=',F4.2)
55CONTINUE N=N-K
IF (N.GE.N3) GO TO 9 STOP
END
Контрольные вопросы
1.Целесообразность использования природного газа в качестве топлива для двигателей внутреннего сгорания.
2.Порядок определения давления и температуры в конце наполнения.
Вкаких пределах лежат эти параметры для двигателей внутреннего сгорания?
3.Как изменяется токсичность отработавших газов в связи с переводом бензинового двигателя на природный газ?
4.Какое влияние оказывает степень сжатия на мощностные и экономические показатели двигателя?
5.Назовите статьи внешнего теплового баланса.
6.Влияние какого параметра на мощностные и экономические показатели двигателя определяется в составленной программе?
7.Какие статьи теплового баланса выводятся на печать в составленной программе?