《内燃机工作过程数值计算程序说明》由会员分享，可在线阅读，更多相关《内燃机工作过程数值计算程序说明（18页珍藏版）》请在金锄头文库上搜索。
1、内燃机工作过程数值计算程序说明,任务说明,内燃机工作过程数值计算教学程序包括 “输入参数.for”给出了发动机运转及结构参数和参数赋初值， “零维程序子程序.for”给出计算所需子程序。任务： 编写主程序，计算给定直喷式柴油机在压缩、燃烧和膨胀过程中缸内状态参数随曲轴转角的变化（计算始点：进气门关闭216 CA,计算终点：排气门打开：504 CA）。进气上止点为0 CA 报告格式： （1）数学模型的建立 （2）程序流程图 （3）程序清单 （4）计算结果分析,变量说明,C -BLOCK DATA- C AR-THE ROTATED RADIUM OF CRANKSHAFT,REAL,(M),0
2、C AL-THE LENGTH OF CONNECTING ROD,REAL,(M),0 C D -ENGINE CYLINDER BORE,REAL,(M),0 C S-ENGINE STROKE DISTANCE,REAL,(M),0 C EPS-COMPRESSION RATIO,REAL,0 C AN-ENGINE SPEED,REAL,(R/MIN),0 C R-UNIVERSAL CONSTANT OF GAS,REAL,(MJ/KG),0 C PA-CYLINDER PRSSURE AT I.V.C.,REAL,(Mpa),0 C TA-CYLINDER TEMPERATURE
3、AT I.V.C.,REAL,(K),0 C TW1-WALL TEMPERATURE OF COMBUSTION CHAMBER ON C PISTON,REAL,(K),0 C TW2-WALL TEMPERATURE OF COMBUSTION CHAMBER ON C CYLINDER HEAD,REAL,(K),0 C TW3-WALL TEMPERATURE OF COMBUSTION CHAMBER ON C CYLINDER WALL,REAL,(K),0,变量说明,C AN1-COEFFICIENT OF C HU-HEAT VALUE OF FUEL,REAL,(KJ/KG
4、),0 C AAM-COEFFICIENCY IN WIBES HEAT RELEASE FUNCTION C ETU-COMBUSTION EFFECIENCY IN WIBES HEAT RELEASE FUNCTION C KKK-CONTROL VARITY (-1 OR 0 OR 1) C -1-COMPRESSION PROCESS C 0-COMBUSION PROCESS C +1-EXPANSION PROCESS C +2-EXHAUST AND INDUCTION PROCESSES C ALM-EXCESS AIR COEFFICIENT C AM-QUANTITY O
5、F GAS IN CYLINDER C AML-QUANTITY OF AIR IN CYLINDER C AMB0-QUANTITY OF INJECTION FUEL PER CYCLE C AL0-CONSTANT QUANTITY OF AIR WHICH 1(KG) FULE COMBUSION C IN THEORATY C AMB-QUANTITY OF INJECTION FULE IN CYLINDER C RAD-CONSTANT PI/180. C PI-CONSTANT 3.1415926,变量说明,C H1,H2,H3,H4-THE CALCULATION STEPS
6、 IN PROCESSES C SPRFI-CRANK ANGLE OF INJECTION FUEL,REAL,(CA),0 AND C ENBEFI-CRANK ANGLE OF EARLIY INDUCTION,REAL,(CA),0 AND0 AND 0 AND 0 AND 180 AND 180 AND 0 AND 360 AND 0 C VH-CYLINDER SWEPT VOLUME,REAL,0 C D-ENGINE CYLINDER BORE,REAL,(M),0 C S-ENGINE STROKE DISTANCE,REAL,(M),0 CFUNCTION DVOLU(FI
7、)COMMON /A/AR,AL,D,EPS,AN,SCOMMON /H/RAD,PIFFI=FI*RADALMS=AR/ALVH=PI*D*D*S/4.0SQRD=SQRT(1.0-(ALMS*SIN(FFI)*2)DVOLU=VH/2.*SIN(FFI)*(1.0+ALMS*COS(FFI)/SQRD)*RADRETURNEND,编程时必须加上,C -VOLU()- C VOLU: THIS PROGRAM IS USED TO CALCULATE THE C VOLUME OF CYLINDER C ALMS-HALF OF THE RATIO OF STROKE DISTANCE TO
8、 THE C LENGTH OF CONNECTING ROD,REAL,0 C VH-CYLINDER SWEPT VOLUME CFUNCTION VOLU(FI)COMMON /A/AR,AL,D,EPS,AN,SCOMMON /H/RAD,PIFFI=FI*RADALMS=AR/ALVH=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)RETURNEND,程序说明,C -BURNT()-P45 公式2.6.9 C C BURNT: THIS PR
9、OGRAM IS USED FOR CALCULATION OF COMBUSTION C HEAT RELEASE VALUES AND VALUE OF INJECTED FUEL C DQB-HEAT RELEASE VALUSES by WEIBE function IN COMBUSION c PRECESS C AMB-QUANTITY OF INJECTION FULE IN CYLINDERSUBROUTINE BURNT(FI)COMMON /D/AN1,HU,AAM,ETUCOMMON /F/BEGFI,FIVB,DFI,ENDFICOMMON /G/AM,AML,AMB0
10、,AL0,AMBCOMMON /H/RAD,PICOMMON /J/DQBC=(FI-FIVB)/DFIAA=EXP(-6.908*(FI-FIVB)/DFI)*(AAM+1.)AMB=(1.0-AA)*AMB0DQB=6.908*ETU*AMB0*HU*(AAM+1.)/DFI*(FI-FIVB)/DFI)*AAM*AARETURNEND,程序说明,C -ALANMD()- C ALANMD: THIS PROGRAM IS USED FOR CALCULATION OF EXCESS C AIR COEFFICIENT AND QUANTITY OF GAS C CONTENT: FUNC
11、TION VOLU(FI) C SUBROUTINE BURNT(FI) CSUBROUTINE ALANMD(FI)COMMON /E/KKK,JJJ,ALMCOMMON /G/AM,AML,AMB0,AL0,AMBCOMMON /F/BEGFI,FIVB,DFI,ENDFICOMMON /B/R,PA,TACOMMON /J/DQBAML=PA*VOLU(BEGFI)/R/TASELECT CASE(KKK)CASE(-1) COMPRESSION PROCESSALM=1.0E4AM=AMLCASE(0) COMBUSION PROCESSCALL BURNT(FI)ALM=AML/(AL0*AMB)AM=AML+AMBCASE(1) EXPANSION PROCESS ALM=AML/(AL0*AMB0)AM=AML+AMB0CASE(2) EXHAUST AND INDUCTION PROCESSESAML=AM-AMBALM=AML/(AL0*AMB)CASE DEFAULTWRITE(*,*) “The working process INDEX doesnt match!“END SELECTRETURNEND,
