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
|
DO 55 JJ = 1, M |
|
|
IRAB=IRAB+1 |
|
|
J=L+(JJ-1)*I |
|
|
S=J/10.0 |
|
C |
IRAB1=0 |
|
C |
14 TR=370*(N/N1)**0.333-10*(N/N1)+(TRN-360) |
|
IRAB1=IRAB1+1 |
||
C |
||
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*N1) |
|
|
||
|
FD=(FN-0.18)+0.18*(N/N1) |
|
|
15 Y=(TO+DT)*PR/(TR*(PA*S*FD-PR)) |
|
|
16 TA=(TO+DT+V*Y*TR)/(1+Y) |
|
|
17 KV=TO*(PA*S*FD-PR)/((TO+DT)*PO*(S-1)) |
|
|
WRITE (6,61) KV |
|
C |
61 FORMAT(1X,'KV=',E11.4) |
|
ПРОЦЕСС СЖАТИЯ |
||
C |
IRAB2=0 |
|
|
18 MCC=21.475+0.00306*(TC1-400) |
|
|
MPC=(MCO*(20.234+0.0031*TC1)+MC*(39.643+0.01323*TC1)+MH2* |
|
|
*(20.619+0.0063*TC1)+MH*(24.615+0.0054*TC1)+MN*(20.129+ |
|
C |
*0.00264*TC1))/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 |
|
|
MPA=(MCO*(20.808+0.00055*TA1)+MC*(27.545+0.02552*TA1)+MH2* |
|
|
*(20.302+0.00318*TA1)+MH*(25.184+0.00242 * |
|
|
*TA1)+MN*(20.704+0.00029*TA1))/M2 |
|
|
MCA=20.759+0.0008*TA1 |
|
|
MAC=(MCA+Y*MPA)/(1+Y) |
|
|
KC=1+8.315*(TC1-TA1)/(MCM*TC1-MAC*TA1) |
|
|
NC=KC-0.02 |
|
|
TC=TA*S**(NC-1) |
|
C |
DTC=((TC-273)-TC1)*100/(TC-273) |
|
WRITE (6,651) DTC,TA1,MCC,MCA |
||
C |
||
C |
651 FORMAT(' DTC=',E11.4,' TA1=',E11.4,' MCC=',E11.4,' MCA=',E11.4) |
|
C |
WRITE (6,751) KC,TA,S,TC |
|
C |
751 FORMAT(' KC=',E11.4,' TA=',E11.4,' S=',E11.4,' TC=',E11.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
22CONTINUE
|
23 PC = PA*S**NC |
|
WRITE (6, 62) PC |
C |
62 FORMAT (' PC = ',E11.4 ) |
ПРОЦЕСС СГОРАНИЯ |
|
|
DHU=114*(1-AF)*LO |
|
KCI=(N/N1)**0.4-0.49*(N/N1)+(KCI1-0.51) |
|
TC1=TC-273 |
|
UCC=8.591+0.02299*(TC1-400) |
|
UPC=(MCO*(0.02199*TC1-0.205)+MC*(0.0415*TC1-2.64)+MH2* |
|
*(0.02119*TC1-1.28)+MH*(0.02947*TC1-1.078)+MN* |
|
*(0.02252*TC1-0.534))/M2 |
|
B=(B1+Y)/(1+Y) |
|
UPZ=(KCI*(HU-DHU)/M1+UCC+Y*UPC)/(B*(1+Y)) |
|
A=(0.02889*MCO+0.05273*MC+0.02805*MH2+0.04648*MH+0.02805*MN)/M2 |
|
B2=(7.162*MCO+13.598*MC+9.961*MH2+21.955*MH+6.232*MN)/M2 |
|
TZ1=(UPZ+B2)/A |
24TZ=TZ1+273
25PZ=PC*B*TZ/TC
LM=PZ/PC WRITE (6,63) TZ
63FORMAT(1X,'TZ=',E11.4)
C |
ПРОЦЕСС РАСШИРЕНИЯ |
|
MPZ=(MCO*(23.196+0.00113*TZ1)+MC*(41.652+0.00226*TZ1)+MH2* |
|
*(19.995+0.00163*TZ1)+MH*(28.402+0.00372*TZ1)+MN* |
|
*(22.735+0.00113*TZ1))/M2 |
C |
IRAB3=0 |
|
30 MPB=(MCO*(21.468+0.00205*TB1)+MC*(36.682+0.0049*TB1)+ |
|
*MH2*(19.532+0.0018*TB1)+MH*(24.975+0.00549*TB1)+ |
|
*MN*(20.922+0.00209*TB1))/M2 |
|
KP=1+8.315*(TZ1-TB1)/(MPZ*TZ1-MPB*TB1) |
|
NP=KP |
|
TB=TZ/S**(NP-1) |
|
DTB=((TB-273)-TB1)*100/(TB-273) |
IF (DTB.GT.2.0) GO TO 31
IF (DTB.LT.0.0) GO TO 32
GO TO 33
31TB1=TB1+5 GO TO 30
32TB1=TB1-5 GO TO 30
33CONTINUE PB=PZ/S**NP WRITE (6, 64) PB
64FORMAT (' PB = ', E11.4) TRI=TB/(PB/PR)**0.333 DTR=(TRN-TRI)*100/TRN IF (DTR.GT.2.0) GO TO 34 IF (DTR.LT.0.0) GO TO 35 GO TO 36
34TRN=TRN-5
GO TO 14
35TRN=TRN+5 GO TO 14
36CONTINUE
C |
ПОКАЗАТЕЛИ РАБОТЫ ЦИКЛА |
|
|
PI1=PC*(LM*(1-1/S**(NP-1))/(NP-1)-(1-1/S**(NC-1))/(NC-1))/(S-1) |
|
|
40 PI=MI*PI1-FI*(PR-PA) |
|
|
GI=3600*KV/(24.4*LT*PI) |
|
|
KI=3600/(GI*HU) |
|
|
PM=0.039+0.0132*SP*N/30 |
|
|
PE=PI-PM |
|
|
KM=PE/PI |
|
|
42 GE=GI/KM |
|
|
43 KE=KI*KM |
|
|
IF (N.EQ.N1.AND.IRAB.EQ.1) GO TO 44 |
|
|
GO TO 45 |
|
|
44 NI1=NE1/KM |
|
|
DC=SQRT(12*NI1*T/(3.14*PI*N1*IC*SP)) |
|
|
VHI=3.14*DC**2*SP*IC*2.5 |
|
|
45 CONTINUE |
|
|
NE=PE*VHI*N/(30*T) |
|
|
ME=3*10**4*NE/(3.14*N) |
|
|
GT=GE*NE/1000 |
|
|
NI=NE/KM |
|
|
WRITE (6, 66) KE |
|
C |
66 FORMAT ( ' KE = ' , E11.4 ) |
|
ТЕПЛОВОЙ БАЛАНС ДВИГАТЕЛЯ |
||
|
Q=GT*HU |
|
|
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*LT*(B*MRP*TR1-MOP*TO1)/1000 |
|
|
QG1=QG*100/Q |
|
|
QH=DHU*GT |
|
|
QH1=QH*100/Q |
|
C |
IRAB4=0 |
|
C |
QS1=0.0 |
|
C |
46 QW=CW*IC*DC**(1+2*MW)*N**MW*(HU-DHU)/(AF*HU) |
|
IRAB4=IRAB4+1 |
||
C |
||
C |
WRITE (6,556) MW,QS1 |
|
C |
556 FORMAT(' MW=',E11.4,' QS1='E11.4) |
|
C |
IF (IRAB4.GT.500) STOP |
|
C----------------------------------- |
QW1=QW*100/Q |
|
|
||
|
QS=Q-(QE+QM+QG+QH+QW) |
|
|
QS1=QS*100/Q |
|
|
IF (QS1.GT.3.0) GO TO 47 |
|
|
IF (QS1.LT.1.0) GO TO 48 |
|
|
GO TO 49 |
|
|
47 MW=MW+0.005 |
GO TO 46
48MW=MW-0.005 GO TO 46
49CONTINUE
WRITE(6,50)
50 FORMAT (1X,30('*')) WRITE (6,70) S,N
70 FORMAT (1X,'S=',E11.4,' N=',E11.4) WRITE (6,51) PZ, GE, KE
WRITE (6,52) NE, ME, GT WRITE (6,53) QE1, QM1, QG1 WRITE (6,54) QH1, QW1, QS1
51 FORMAT (3X,'PZ=' ,F5.3,2X,'GE=' ,F5.1,2X,'KE=' ,F6.4)
52 FORMAT (3X,'NE=' ,F6.1,2X,'ME=' ,F5.1,2X,'GT=' ,F5.2)
53 FORMAT (3X,'QE1=',F6.3,2X,'QМ1=',F5.3,2X,'QG1=',F6.3)
54 FORMAT (3X,'QH1=',F5.2,2X,'QW1=',F5.2,2X,'QS1=',F4.2)
55 CONTINUE N=N-K
IF (N. GE. N3) GO TO 9 STOP
END
Контрольные вопросы
1.Экономическая эффективность использования сжиженного газа в двигателях с воспламенением рабочей смеси от электрической искры.
2.Как изменится токсичность двигателя при переводе его на питание
сбензина на сжиженный газ?
3.Достоинства и недостатки двигателя, работающего на сжиженном
газе.
4.Чем ограничивается степень сжатия при использовании легких топлив (бензина, сжиженного газа, природного газа)?
5.Что называется литровой мощностью двигателя, и какие факторы влияют на её величину?
6.Влияние какого параметра на мощностные и экономические показатели работы двигателя определяется в программе?
7.Какие показатели работы двигателя выведены на печать в составленной программе?