内燃机工作过程数值计算程序说明.ppt

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/K

4、G),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

5、OF 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 STEP

6、S 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 540,BLOCK DATA OPRATION COMMON /A/AR,AL,D,EPS,AN,S COMMON /B/R,PA,TA COMMON /C/TW1,TW2,TW3 COMMON /D/AN1,HU,AAM,ETU COMMON /E/K

7、KK,JJJ,ALM COMMON /F/BEGFI,FIVB,DFI,ENDFI COMMON /G/AM,AML,AMB0,AL0,AMB COMMON /H/RAD,PI COMMON /J/DQB COMMON /K/H1,H2,H3,H4 COMMON /L/SPRFI,ENBEFI,ENAFFI,EXBEFI,EXAFFI COMMON /M/TE,TB,PE,PB COMMON /N/W,变量说明,有名公用区，给不同程序单位的变量分配同一存储单元,COMMON语句，各程序单位数据交换的一种方式,数据块子程序，它的作用是初始化有名公用区中的变量,DATA AR,AL,D,EPS,A

8、N,S/0.0575,0.190,0.10,18.0,2300.,0.115/ DATA R,PA,TA/0.28706E-3,0.09,330./ DATA TW1,TW2,TW3/543.,493.,473./ DATA AN1,HU,AAM,ETU/1.35,44100.,0.5,1./ DATA FIVB,DFI/353.,70./ DATA AMB0,AL0/3.3E-5,14.4/ DATA RAD,PI/0.0174532,3.1415926/ DATA H1,H2,H3,H4/4.,0.4,4.0,2.0/ DATA SPRFI,ENBEFI,ENAFFI,EXBEFI,EXA

9、FFI/25.,16.,36.,36.,14./ DATA TE,TB,PE,PB/300.,300.,0.1,0.1/ END,变量说明,程序说明,C -DVOLU()- C C DVOLU: THIS PROGRAM IS USED FOR THE CALCULATION OF C CYLINDER VOLUME VARY WITH THE CRANK ANGLE C ALMS-HALF OF THE RATIO OF STROKE DISTANCE TO THE C LENGTH OF CONNECTING ROD,REAL,0 C VH-CYLINDER SWEPT VOLUME,RE

10、AL,0 C D-ENGINE CYLINDER BORE,REAL,(M),0 C S-ENGINE STROKE DISTANCE,REAL,(M),0 C 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 END,编程时必须加上,C -VOLU()- C VOLU

11、: THIS PROGRAM IS USED TO CALCULATE THE C VOLUME OF CYLINDER C ALMS-HALF OF THE RATIO OF STROKE DISTANCE TO THE C LENGTH OF CONNECTING ROD,REAL,0 C VH-CYLINDER SWEPT VOLUME C 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(FF

12、I)*2) VOLU=VH/2.*(2./(EPS-1.)+1.0-COS(FFI)+(1.0-SQRD)/ALMS) RETURN END,程序说明,C -BURNT()-P45 公式2.6.9 C C BURNT: THIS PROGRAM 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 INJE

13、CTION FULE IN CYLINDER 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,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

14、,程序说明,C -ALANMD()- C ALANMD: THIS PROGRAM IS USED FOR CALCULATION OF EXCESS C AIR COEFFICIENT AND QUANTITY OF GAS C 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

15、/DQB AML=PA*VOLU(BEGFI)/R/TA SELECT CASE(KKK) CASE(-1) COMPRESSION PROCESS ALM=1.0E4 AM=AML CASE(0) COMBUSION PROCESS CALL BURNT(FI) ALM=AML/(AL0*AMB) AM=AML+AMB CASE(1) EXPANSION PROCESS ALM=AML/(AL0*AMB0) AM=AML+AMB0 CASE(2) EXHAUST AND INDUCTION PROCESSES AML=AM-AMB ALM=AML/(AL0*AMB) CASE DEFAULT

16、 WRITE(*,*) The working process INDEX doesnt match! END SELECT RETURN END,计算过量空气系数,C -HECOEF()- C C HECOEF: THIS PROGRAM IS USED FOR CALCULATION OF HEAT TRANSFER C COEFFICIENT CYLINDER C WSCHINIS FROMULA C CONTENT: C FUNCTION VOLU(FI) C FUNCTION PRES(FI,T) C SUBROUTINE ALANMD(FI) C FUNCTION HECOEF(F

17、I,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) VH=PI/4.*D*D*S CM=AN*S/30. 活塞平均速度 m/s C=0.,程序说明,程序说明,SELECT CASE(KKK) CASE(-1) C1=2.28 CASE(0,1) C2=3.24E-3 P

18、0=(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! 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()- C C HETRAN: THIS PROGRAM IS USED FOR CAL

19、CUATION OF VALUES OF HEAT C TRANFER IN CYLINDER VARY WITH CRANK ANGLE C FUNCTION HETRAN(FI,T) COMMON /A/AR,AL,D,EPS,AN,S COMMON /C/TW1,TW2,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,C -PROPER()- C C PROPER: T

20、HIS PROGRAM IS USED FOR CALCULATION OF THERMAL C PROPERTY OF GAS C WRITTEN BY: DIVISION OF INTERNAL COMBUSTION ENGINE XIAN JIAOTONG UNIVERSITY C SUBROUTINE PROPER(FI,T,CV,U,DULM) COMMON /E/KKK,JJJ,ALM A1=0.0975 A2=0.0485A0=0.14455 B1=7.768 B2=3.36 C1=489.6 C2=46.4 D1=1356.8 T1=(T-273)*1.E-2 CALL ALA

21、NMD(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-B2*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) RETURN END,比内能解析式 （见第20页2.3.16） 单位：kJ/kg 等容比热解析式

22、（见第20页2.3.17） 单位：kJ/(kg.K),C -FOMU()- C FOMU: THIS PROGRAM IS USED FOR SIMULATION OF COMPRESSION, C COMBUSTION AND EXPANSION PROCESSES C 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) COMMON

23、/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 单位转化为 kJ/CA WORK=PRES(FI,T)*DVOLU(FI)*1.0E3 IF(KKK.EQ.0) THEN CALL BURNT(FI) BURN=DQB ENDIF FOMU=BURN+HEAT-WORK-U*BURN/HU FOMU=1.0/(AM*CV)*(FOMU+AM*DULM*AML/(AL0*AMB*2*HU)*BURN) RE