Материал: 1524
n=n+7.9.Программа теплового расчета газодизельного двигателя
REAL M1,MC,MH,MN,MO,M2,C,H,O,MD,AF1,N1,NB,MB,M1B,MT,MB1,KCI1, *NE1,MI,IC,CW,MW,LO1,LO2,LO,AF,LT,RO,PA,PR,TR,KV,N,ME,KI, *MCA,MCC,KC,NC,KCI,UCC,UPC1,UPC,MPZ,MPD1,MPB,KP,NP,TB,PB,
*PI,KM,KE,GE,NI,NE,GT,QE1,QM1,QG1,QW1,QS1,LM,K,N3,MCB,MPB1 OPEN(5,FILE='ISX_DAN.DAT')
OPEN(6,FILE='REZYL.TXT') C OPEN(7,FILE='CON')
READ (5,1)C,H,O,MD,AF1,N1,NB,MB,M1B,MT,MB1,TO,PO,RB,CN, *FK,BS,TRN,DT1,S,V,TC1,KCI1,HU,LM,TB1,MI,FI,NE1,T,IC,CW,
*MW,K,N3 READ (5,3) L,M,I
WRITE(6,2)C,H,O,MD,AF1,N1,NB,MB,M1B,MT,MB1,TO,PO,RB,CN,
*FK,BS,TRN,DT1,S,V,TC1,KCI1,HU,LM,TB1,MI,FI,NE1,T,IC,CW,
*MW,K,N3
WRITE(6,4)L,M,I
1FORMAT (7(/),(63X,F8.0))
2FORMAT (1X,F12.5)
3FORMAT (63X,I8)
4FORMAT (1X,I12)
5FORMAT (63X,F8.0)
C |
ПАРАМЕТРЫ РАБОЧЕГО ТЕЛА |
|
N=N1 |
|
LO1=(C+3*(H-O/8))/(0.21*12) |
|
LO2=(NB+MB/4)/(0.21*MB1) |
|
LO=MD*LO1+M1B*LO2 |
|
N=N1 |
C------- |
10 AF=AF1*(1.672*(1+N/N1)+(N/N1)**2) |
10 AF=AF1
|
LT=AF*LO |
|
M1=LT+MD/MT+M1B/MB1 |
|
MC=MD*C/12+M1B*NB/MB1 |
|
MH=MD*H/2+M1B*MB/(2*MB1) |
|
MN=0.79*AF*LO |
|
MO=0.21*(AF-1)*LO |
|
M2=MC+MH+MN+MO |
|
B1=M2/M1 |
C |
ПРОЦЕСС ВПУСКА |
|
RO=PO*10**6/(RB*TO) |
|
SP=30*CN/N1 |
|
W=0.05433*SP*N1*FK |
|
IRAB=0 |
|
DO 46 JJ=1,M |
IRAB=IRAB+1
J=L+(JJ-1)*I
BS=J/10
11 PA=PO-BS*W**2*RO/(2*10**6)
PR1=0.098+0.539*N1/(10**5) PR=1.035*PO+(PR1-1.035*PO)*(N/N1)**2
C IRAB1=0
12 TR=(TRN+45)-90*(N/N1-0.5)+273
DT=DT1*(110+0.0125*N)/(110+0.0125*N1)
13 Y=(TO+DT)*PR/(TRN*(PA*S-PR))
14 TA=(TO+DT+Y*V*TR)/(1+Y)
15 KV=TO*(PA*S-PR)/((TO+DT)*PO*(S-1)) C------------------------
|
WRITE(6,51)KV |
51 FORMAT('KV=',E10.4) |
|
C------------------------ |
|
C |
ПРОЦЕСС СЖАТИЯ |
|
TA1=TA-273 |
|
MCA=20.758+0.0008*TA1 |
C |
IRAB2=0 |
C |
DTC=0.0 |
C |
16 MCC=22.09+0.00318*(TC1-600) |
|
|
C |
IRAB2=IRAB2+1 |
C |
WRITE (6,551) TC1 |
C 551 FORMAT(/,' TC1=',E10.4) |
|
C |
IF (IRAB2.GT.500) STOP |
C----------------------------------- |
KC=1+8.315*(TC1-TA1)/(MCC*TC1-MCA*TA1) |
|
|
|
NC=KC |
17 TC=TA*S**(NC-1) |
|
C |
DTC=((TC-273)-TC1)*100/(TC-273) |
|
|
C |
WRITE (6,651) DTC,TA1,MCC,MCA |
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.5.0) GO TO 18
IF(DTC.LT.0.0) GO TO 19 GO TO 20
18 TC1=TC1+5
GO TO 16
19 TC1=TC1-5
GO TO 16
20 CONTINUE
PC=PA*S**NC
C-------------------------
|
WRITE(6,52)PC |
52 FORMAT ('PC=',E10.4) |
|
C------------------------- |
|
C |
ПРОЦЕСС СГОРАНИЯ |
B=(B1+Y)/(1+Y)
KCI=(KCI1-0.15)+(N/N1-0.3)**0.75-0.878*(N/N1-0.3) TC1=TC-273
UCC=13.255+0.0243*(TC1-600)
UPC1=14.779+0.02806*(TC1-600)
UPC=UPC1/AF+UCC*(AF-1)/AF H1=((KCI*HU/M1+UCC+Y*UPC)/(1+Y)+0.008315*LM*TC)/B TZ1=(H1+(3.59+2.942/AF))/(0.036575+0.00524/AF)
21TZ=TZ1+273
22PZ=LM*PC
R=B*TZ/(LM*TC)
D=S/R
C------------------------
|
WRITE(6,53)TZ |
53 FORMAT('TZ=',E10.4) |
|
C------------------------ |
|
C |
ПРОЦЕСС РАСШИРЕНИЯ |
|
MPZ1=28.638+0.0023*(TZ1-1800) |
|
MCZ=25.003+0.001615*(TZ1-1800) |
|
MPZ=MPZ1/AF+MCZ*(AF-1)/AF |
C |
IRAB3=0 |
C |
DTB=0.0 |
23 MPB1=25.498+0.00376*(TB1-800)
C-----------------------------------
C IRAB3=IRAB3+1
C WRITE (6,552) TB1
C 552 FORMAT(/,' TB1=',E10.4)
C IF (IRAB3.GT.500) STOP
C-----------------------------------
MCB=22.713+0.00293*(TB1-800)
MPB=MPB1/AF+MCB*(AF-1)/AF
KP=1+8.315*(TZ1-TB1)/(MPZ*TZ1-MPB*TB1)
NP=KP
TB=TZ/D**(NP-1)
DTB=((TB-273)-TB1)*100/(TB-273)
C-----------------------------------
C WRITE (6,652) MCB,MPB,MPB1,AF
C 652 FORMAT(' MCB=',E11.4,' MPB=',E11.4,' MPB1=',E11.4,' AF=',E11.4) C WRITE (6,752) KP,TB,DTB
C 752 FORMAT(' KP=',E11.4,' TB=',E11.4,' DTB=',E11.4) C-----------------------------------
IF(DTB.GT.7.0)GO TO 24
IF(DTB.LT.0.0)GO TO 25 GO TO 26
24TB1=TB1+5 GO TO 23
25TB1=TB1-5 GO TO 23
26CONTINUE PB=PZ/D**NP TRI=TB/(PB/PR)**0.333
DTR=(TR-TRI)*100/TR IF(DTR.GT.5.0)GO TO 27 IF(DTR.LT.0.0)GO TO 28 GO TO 29
27TRN=TRN-5
GO TO 12
28TRN=TRN+5 GO TO 12
29CONTINUE
C--------------------------
|
WRITE (6, 54) PB |
54 FORMAT ('PB=', E10.4) |
|
C-------------------------- |
|
C |
ПОКАЗАТЕЛИ РАБОТЫ ЦИКЛА |
PI1=PC*(LM*(R-1)+LM*R*(1-1/D**(NP-1))/(NP-1)-(1-1/S**(NC-1))/ *(NC-1))/(S-1)
PI=MI*PI1-FI*(PR-PA)
GI=3600*KV/(24.4*LT*PI)
KI=3600/(GI*HU)
PM=0.089+0.0118*CN PE=PI-PM KM=PE/PI KE=KI*KM
31GE=GI/KM IF(N.EQ.N1.AND.IRAB.EQ.1.0)GO TO 32
GO TO 33
32NI1=NE1/KM