1、燃烧放热规律计算程序说明燃烧放热规律计算程序说明任务说明任务说明内燃机工作过程数值计算内燃机工作过程数值计算教学程序包括教学程序包括“输入参数输入参数.for”给出了发动机运转及结构参数和参数赋初值，给出了发动机运转及结构参数和参数赋初值，“零维程序子程序零维程序子程序.for”给出计算所需子程序。给出计算所需子程序。 任务：任务：编写主程序，计算给定直喷式柴油机在压缩、燃烧和膨胀过程编写主程序，计算给定直喷式柴油机在压缩、燃烧和膨胀过程中缸内状态参数随曲轴转角的变化（计算始点：进气门关闭中缸内状态参数随曲轴转角的变化（计算始点：进气门关闭216 CA,计算终点：排气门打开：计算终点：排气门打
2、开：504 CA）。）。报告格式：报告格式：（1）数学模型的建立）数学模型的建立（2）程序流程图）程序流程图（3）程序清单）程序清单（4）计算结果分析）计算结果分析变量说明变量说明C -BLOCK DATA-C AR-THE ROTATED RADIUM OF CRANKSHAFT,REAL,(M),0C AL-THE LENGTH OF CONNECTING ROD,REAL,(M),0C D -ENGINE CYLINDER BORE,REAL,(M),0C S-ENGINE STROKE DISTANCE,REAL,(M),0C EPS-COMPRESSION RATIO,REAL,0C
3、 AN-ENGINE SPEED,REAL,(R/MIN),0C R-UNIVERSAL CONSTANT OF GAS,REAL,(MJ/KG),0C PA-CYLINDER PRSSURE AT I.V.C.,RAEL,(Mpa),0C TA-CYLINDER TEMPERATURE AT I.V.C.,REAL,(K),0C TW1-WALL TEMPERATURE OF COMBUSTION CHAMBER ON C PISTON,RAEL,(K),0 C TW2-WALL TEMPERATURE OF COMBUSTION CHAMBER ONC CYLINDER HEAD,REAL
4、,(K),0C TW3-WALL TEMPERATURE OF COMBUSTION CHAMBER ONC CYLINDER WALL,RAEL,(K),0变量说明变量说明C AN1-COEFFICIENT OF C HU-HEAT VALUE OF FUEL,REAL,(KJ/KG),0C AAM-COEFFICIENCY IN WIBES HEAT RELEASE FUNCTIONC ETU-COMBUSTION EFFECIENCY IN WIBES HEAT RELEASE FUNCTIONC KKK-CONTROL VARITY (-1 OR 0 OR 1)C -1-COMPRES
5、SION PROCESSC 0-COMBUSION PROCESSC +1-EXPANSION PROCESSC +2-EXHAUST AND INDUCTION PROCESSESC ALM-EXCESS AIR COEFFICIENTC AM-QUANTITY OF GAS IN CYLINDERC AML-QUANTITY OF AIR IN CYLINDERC AMB0-QUANTITY OF INJECTION FUEL PER CYCLEC AL0-CONSTANT QUANTITY OF AIR WHICH 1(KG) FULE COMBUSIONC IN THEORATY C
6、AMB-QUANTITY OF INJECTION FULE IN CYLINDERC RAD-CONSTANT PI/180.C PI-CONSTANT 3.1415926变量说明变量说明C H1,H2,H3,H4-THE CALCULATION STEPS IN PROCESSESC SPRFI-CRANK ANGLE OF INJECTION FULE,RAEL,(CA),0 ANDC ENBEFI-CRANK ANGLE OF EARLIY INDUCTION,REAL,(CA),0 AND0 AND 0 AND 0 AND 180 AND 180 AND 0 AND 360 AND
7、0C VH-CYLINDER SWEPT VOLUME,REAL,0C D-ENGINE CYLINDER BORE,REAL,(M),0C S-ENGINE STROKE DISTANCE,REAL,(M),0C FUNCTION DVOLU(FI) COMMON /A/AR,AL,D,EPS,AN,S COMMON /H/RAD,PI FFI=FI*RAD ALMS=AR/AL VH=PI*D*D*S/4.0 SQRD=SQRT(1.0-(ALMS*SIN(FFI)*2) DVOLU=VH/2.*SIN(FFI)*(1.0+ALMS*COS(FFI)/SQRD)*RAD RETURN EN
8、D180sin1cossinsin222sshVddV编程时必须加上编程时必须加上C -FOMU()-C FOMU: THIS PROGRAM IS USED FOR SIMULATION OF COMPRESSION,C COMBUSTION AND EXPANSION PROCESSESC CONTENTS: FUNCTION HETRAN(FI,T)C FUNCTION PRES(FI,T)C FUNCTION DVOLU(FI)C SUBROUTINE PROPER(FI,T,CV,U,DULM)C SUBROUTINE BURNT(FI)C FUNCTION FOMU(FI,T) C
9、OMMON /D/AN1,HU,M,ETU COMMON /G/AM,AML,AMB0,AL0,AMB COMMON /E/KKK,JJJ,ALM COMMON /J/DQB AMB=1.0 BURN=0. CALL PROPER(FI,T,CV,U,DULM) HEAT=HETRAN(FI,T)/1.0E3 WORK=PRES(FI,T)*DVOLU(FI)*1.0E3 IF(KKK.EQ.2) THEN CALL BURNT(FI) BURN=DQB ENDIF FOMU=BURN+HEAT-WORK-U*BURN/HU FOMU=1.0/(AM*CV)*(FOMU-AM*DULM*AML
10、/(AL0*AMB*2*HU)*BURN) RETURN ENDddumddmuddmhddmhddVpddQddQmcddTAEWB21C -VOLU()-C VOLU: THIS PROGRAM IS USED TO CALCULATE THEC VOLUME OF CYLINDERC ALMS-HALF OF THE RATIO OF STROKE DISTANCE TO THE C LENGTH OF CONNECTING ROD,REAL,0C VH-CYLINDER SWEPT VOLUMEC FUNCTION VOLU(FI) COMMON /A/AR,AL,D,EPS,AN,S
11、 COMMON /H/RAD,PI FFI=FI*RAD ALMS=AR/AL VH=PI*D*D*S/4. SQRD=SQRT(1.0-(ALMS*SIN(FFI)*2) VOLU=VH/2.*(2./(EPS-1.)+1.0-COS(FFI)+(1.0-SQRD)/ALMS) RETURN END22sin111cos1122)(sshVV程序说明程序说明C -VOLU()-C VOLU: THIS PROGRAM IS USED TO CALCULATE THEC VOLUME OF CYLINDERC ALMS-HALF OF THE RATIO OF STROKE DISTANCE
12、TO THE C LENGTH OF CONNECTING ROD,REAL,0C VH-CYLINDER SWEPT VOLUMEC FUNCTION VOLU(FI) COMMON /A/AR,AL,D,EPS,AN,S COMMON /H/RAD,PI FFI=FI*RAD ALMS=AR/AL VH=PI*D*D*S/4. SQRD=SQRT(1.0-(ALMS*SIN(FFI)*2) VOLU=VH/2.*(2./(EPS-1.)+1.0-COS(FFI)+(1.0-SQRD)/ALMS) RETURN END22sin111cos1122)(sshVV程序说明程序说明C -BURN
13、T()-CC BURNT: THIS PROGRAM IS USED FOR CALCULATION OF COMBUSTION C HEAT RELEASE VALUES AND VALUE OF INJECTED FUELC DQB-HEAT RELEASE VALUSES by WEIBE function IN COMBUSION PRECESSC SUBROUTINE BURNT(FI) COMMON /D/AN1,HU,AAM,ETU COMMON /F/BEGFI,FIVB,DFI,ENDFI COMMON /G/AM,AML,AMB0,AL0,AMB COMMON /H/RAD
14、,PI COMMON /J/DQB C=(FI-FIVB)/DFI AA=EXP(-6.908*(FI-FIVB)/DFI)*(AAM+1.) AMB=(1.0-AA)*AMB0 DQB=6.908*ETU*AMB0*HU*(AAM+1.)/DFI*(FI-FIVB)/DFI)*AAM*AA RETURN END程序说明程序说明10908. 6exp) 1(908. 6mVBmVBuBuBmHmddQC -ALANMD()-CC ALANMD: THIS PROGRAM IS USED FOR CALCULATION OF EXCESSC AIR COEFFICIENT AND QUANTIT
15、Y OF GASC CONTENT: FUNCTION VOLU(FI)C SUBROUTINE BURNT(FI)C SUBROUTINE ALANMD(FI) COMMON /E/KKK,JJJ,ALM COMMON /G/AM,AML,AMB0,AL0,AMB COMMON /F/BEGFI,FIVB,DFI,ENDFI COMMON /B/R,PA,TA COMMON /J/DQB AML=PA*VOLU(BEGFI)/R/TA SELECT CASE(KKK) CASE(-1) ALM=1.0E4 AM=AML CASE(0) CALL BURNT(FI) ALM=AML/(AL0*
16、AMB) AM=AML+AMB CASE(1) ALM=AML/(AL0*AMB0) AM=AML+AMB0 CASE(2) AML=AM-AMB ALM=AML/(AL0*AMB) CASE DEFAULT WRITE(*,*) The working process INDEX doesnt match! END SELECT RETURN END程序说明程序说明计算过量空气系数计算过量空气系数C -HECOEF()-CC HECOEF: THIS PROGRAM IS USED FOR CALCULATION OF HEAT TRANSFERC COEFFICIENT CYLINDERC
17、 WSCHINIS FROMULAC CONTENT: C FUNCTION VOLU(FI)C FUNCTION PRES(FI,T)C SUBROUTINE ALANMD(FI)C FUNCTION HECOEF(FI,T) COMMON /A/AR,AL,D,EPS,AN,S COMMON /B/R,PA,TA COMMON /D/AN1,HU,AAM,ETU COMMON /E/KKK,JJJ,ALM COMMON /F/BEGFI,FIVB,DFI,ENDFI COMMON /G/AM,AML,AMB0,AL0,AMB COMMON /H/RAD,PI CALL ALANMD(FI)
18、 VH=PI/4.*D*D*S CM=AN*S/30. C=0. 程序说明程序说明K)W/(m,)(82028 . 00212 . 053. 08 . 0ppVpVTCCCDTpaaham程序说明程序说明SELECT CASE(KKK) CASE(-1) C1=2.28 CASE(0,1) C2=3.24E-3 P0=(VOLU(BEGFI)/VOLU(FI)*AN1*PA C=C2*VH/(AM*R)*(PRES(FI,T)-P0) CASE(2) C1=6.18 CASE DEFAULT WRITE(*,*) The working process INDEX doesnt match!
19、END SELECT C=(C1*CM+C)*0.8 P=PRES(FI,T) HECOEF=820./T*0.53*PRES(FI,T)*0.8/D*0.2*C RETURN END用绝热过程代替倒拖压力用绝热过程代替倒拖压力C -HETRAN()-CC HETRAN: THIS PROGRAM IS USED FOR CALCUATION OF VALUES OF HEATC TRANFER IN CYLINDER VARY WITH CRANK ANGLEC FUNCTION HETRAN(FI,T) COMMON /A/AR,AL,D,EPS,AN,S COMMON /C/TW1,TW
20、2,TW3 COMMON /H/RAD,PI AA1=PI*D*D/4. AA2=AA1 AA3=4.*VOLU(FI)/D HETRAN=HECOEF(FI,T)*(TW1-T)*AA1+(TW2-T)*AA2+(TW3-T)*AA3)/6./AN RETURN END)(61)(1313131TTAnTTAddQddQWiiiWiiiiWiWC -PROPER()-CC PROPER: THIS PROGRAM IS USED FOR CALCULATION OF THERMALC PROPERTY OF GASC WRITTEN BY: DIVISION OF INTERNAL COMB
21、USTION ENGINE XIAN JIAOTONG UNIVERSITYC SUBROUTINE PROPER(FI,T,CV,U,DULM)COMMON /E/KKK,JJJ,ALMA1=0.0975A2=0.0485A0=0.14455B1=7.768B2=3.36C1=489.6C2=46.4D1=1356.8T1=(T-273)*1.E-2CALL ALANMD(FI)U=A0*(-(A1+A2/ALM*0.75)*T1*3+(B1+B2/ALM*0.8)*T1*2+(C1+C2 $ /ALM*0.93)*T1+D1)DULM=A0*(0.75*A2/ALM*1.75*T1*3-B
22、2*0.8/ALM*1.8*T1*2-0.93* $ C2/ALM*1.93*T1)CV=A0*(-3*(A1+A2/ALM*0.75)*T1*2*1.E-2+2*(B1+B2/ALM*0.8)* $ T1*1.E-2+(C1+C2/ALM*0.93)*1.E-2)RETURNEND比内能解析式比内能解析式 （见第（见第20页页2.3.16）等容比热解析式（见第等容比热解析式（见第20页页2.3.17）C -RUGKUT()-SUBROUTINE RUGKUT(XVALUE,YVALUE,XSTEP,FUN) EXTERNAL FUN DIMENSION XVALUE(2),YVALUE(2)
23、,AUX(5) X =XVALUE(1) XEND=XVALUE(2) Y =YVALUE(1) DO WHILE(X+XSTEP).LE.XEND) K=1 X1=X Y1=Y DO WHILE(K.LE.4) AUX(K)=FUN(X,Y) SELECT CASE(K) CASE(1,2) COER=2.0 CASE(3) COER=1.0 CASE(4) K=5 AUX(K)=(AUX(1)+2.0*AUX(2)+2.0*AUX(3)+AUX(4)/6.0 END SELECT X=X1+XSTEP/COER Y=Y1+AUX(K)*XSTEP/COER K=K+1 END DO END
24、 DO YVALUE(2)=Y RETURN END主程序主程序 PROGRAM MAIN EXTERNAL FOMUCOMMON /E/KKK,JJJ,ALM COMMON /F/BEGFI,FIVB,DFI,ENDFICOMMON /B/R,PA,TA BEGFI=216. ENDFI=504. FI=216. T=330. 11 IF(FI.GT.353.)THEN IF(FI.GE.423.)THEN IF(FI.GT.504.)THEN STOP END ELSE kkk=1 end if ELSE kkk=0 END IF ELSE kkk=-1 END IF CALL RUGKUT(FI,T,1.,FOMU) CALL PRINT(FI,T,H,W) IF(FI.LE.503.) GOTO 11 END
