内燃机工作过程数值计算程序说明精讲课件

燃烧放热规律计算程序说明燃烧放热规律计算程序说明1任务说明《内燃机工作过程数值计算》教学程序包括“输入参数.for”给出了发动机运转及结构参数和参数赋初值，“零维程序子程序.for”给出计算所需子程序。

任务：编写主程序，计算给定直喷式柴油机在压缩、燃烧和膨胀过程中缸内状态参数随曲轴转角的变化（计算始点：进气门关闭216CA,计算终点：排气门打开：504CA）。报告格式：（1）数学模型的建立（2）程序流程图（3）程序清单（4）计算结果分析任务说明《内燃机工作过程数值计算》教学程序包括2变量说明C------------------------BLOCKDATA-----------------------------CAR--------THEROTATEDRADIUMOFCRANKSHAFT,REAL,(M),>0CAL--------THELENGTHOFCONNECTINGROD,REAL,(M),>0CD--------ENGINECYLINDERBORE,REAL,(M),>0CS---------ENGINESTROKEDISTANCE,REAL,(M),>0CEPS-------COMPRESSIONRATIO,REAL,>0CAN--------ENGINESPEED,REAL,(R/MIN),>0CR---------UNIVERSALCONSTANTOFGAS,REAL,(MJ/KG),>0CPA-------CYLINDERPRSSUREATI.V.C.,RAEL,(Mpa),>0CTA-------CYLINDERTEMPERATUREATI.V.C.,REAL,(K),>0CTW1-----WALLTEMPERATUREOFCOMBUSTIONCHAMBERONCPISTON,RAEL,(K),>0CTW2-------WALLTEMPERATUREOFCOMBUSTIONCHAMBERONCCYLINDERHEAD,REAL,(K),>0CTW3-------WALLTEMPERATUREOFCOMBUSTIONCHAMBERONCCYLINDERWALL,RAEL,(K),>0变量说明C------------------3变量说明CAN1-------COEFFICIENTOFCHU--------HEATVALUEOFFUEL,REAL,(KJ/KG),>0CAAM-------COEFFICIENCYINWIBE'SHEATRELEASEFUNCTIONCETU----COMBUSTIONEFFECIENCYINWIBE'SHEATRELEASEFUNCTIONCKKK-------CONTROLVARITY(-1OR0OR1)C-1--------------COMPRESSIONPROCESSC0--------------COMBUSIONPROCESSC+1--------------EXPANSIONPROCESSC+2--------------EXHAUSTANDINDUCTIONPROCESSESCALM-------EXCESSAIRCOEFFICIENTCAM--------QUANTITYOFGASINCYLINDERCAML-------QUANTITYOFAIRINCYLINDERCAMB0------QUANTITYOFINJECTIONFUELPERCYCLECAL0-------CONSTANTQUANTITYOFAIRWHICH1(KG)FULECOMBUSIONCINTHEORATYCAMB-------QUANTITYOFINJECTIONFULEINCYLINDERCRAD-------CONSTANTPI/180.CPI--------CONSTANT3.1415926变量说明CAN1-------COEFFICI4变量说明CH1,H2,H3,H4---THECALCULATIONSTEPSINPROCESSESCSPRFI-----CRANKANGLEOFINJECTIONFULE,RAEL,(CA),>0ANDCENBEFI----CRANKANGLEOFEARLIYINDUCTION,REAL,(CA),>0AND<100CENAFFI---CRANKANGLEOFLATERINDUCTION,REAL,(CA),>0AND<100CEXBEFI----CRANKANGLEOFEARLIYEXHAUSTION,RAEL,(CA),>0AND<100CEXAFFI----CRANKANGLEOFLATEREXHAUSTION,REAL,(CA),>0AND<100CBEGFI-----CRANKANGLEOFTHEBEGINNINGOFCALCULATION,REAL,C(CA),>180AND<360CFIVB------CRANKANGLEINWIBE'SHEATRELEASE

FUNCTION,RAEL,(CA)C>180AND<360CDFI-------CRANKANGLEINWIBE'SHEATRELEASEFUNCTION,RAEL,(CA)C>0AND<100CENDFI-----CRANKANGLEATE.V.O.,RAEL,(CA),>360AND<540变量说明CH1,H2,H3,H4---THE5 BLOCKDATAOPRATION 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/KKK,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

变量说明 BLOCKDATAOPRATION变量说明6DATAAR,AL,D,EPS,AN,S/0.0575,0.190,0.10,18.0,2300.,0.115/DATAR,PA,TA/0.28706E-3,0.09,330./DATATW1,TW2,TW3/543.,493.,473./DATAAN1,HU,AAM,ETU/1.35,44100.,0.5,1./DATAFIVB,DFI/353.,70./DATAAMB0,AL0/3.3E-5,14.4/DATARAD,PI/0.0174532,3.1415926/DATAH1,H2,H3,H4/4.,0.4,4.0,2.0/DATASPRFI,ENBEFI,ENAFFI,EXBEFI,EXAFFI/25.,16.,36.,36.,14./DATATE,TB,PE,PB/300.,300.,0.1,0.1/END变量说明DATAAR,AL,D,EPS,AN,S/0.0575,07程序说明C----------------DVOLU()------------------------------------CCDVOLU:THISPROGRAMISUSEDFORTHECALCULATIONOFCCYLINDERVOLUMEVARYWITHTHECRANKANGLECALMS------HALFOFTHERATIOOFSTROKEDISTANCETOTHECLENGTHOFCONNECTINGROD,REAL,>0CVH--------CYLINDERSWEPTVOLUME,REAL,>0CD---------ENGINECYLINDERBORE,REAL,(M),>0CS---------ENGINESTROKEDISTANCE,REAL,(M),>0CFUNCTIONDVOLU(FI)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----------------DV8C-------------------------FOMU()--------------------------------CFOMU:THISPROGRAMISUSEDFORSIMULATIONOFCOMPRESSION,CCOMBUSTIONANDEXPANSIONPROCESSESCCONTENTS:FUNCTIONHETRAN(FI,T)CFUNCTIONPRES(FI,T)CFUNCTIONDVOLU(FI)CSUBROUTINEPROPER(FI,T,CV,U,DULM)CSUBROUTINEBURNT(FI)CFUNCTIONFOMU(FI,T)COMMON/D/AN1,HU,M,ETUCOMMON/G/AM,AML,AMB0,AL0,AMBCOMMON/E/KKK,JJJ,ALMCOMMON/J/DQBAMB=1.0BURN=0.CALLPROPER(FI,T,CV,U,DULM)HEAT=HETRAN(FI,T)/1.0E3WORK=PRES(FI,T)*DVOLU(FI)*1.0E3IF(KKK.EQ.2)THENCALLBURNT(FI)BURN=DQBENDIFFOMU=BURN+HEAT-WORK-U*BURN/HUFOMU=1.0/(AM*CV)*(FOMU-AM*DULM*AML/(AL0*AMB**2*HU)*BURN)RETURNENDC----------------------9C---------------VOLU()--------------------------------------CVOLU:THISPROGRAMISUSEDTOCALCULATETHECVOLUMEOFCYLINDERCALMS------HALFOFTHERATIOOFSTROKEDISTANCETOTHECLENGTHOFCONNECTINGROD,REAL,>0CVH--------CYLINDERSWEPTVOLUMECFUNCTIONVOLU(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---------------VOLU()-10C---------------VOLU()--------------------------------------CVOLU:THISPROGRAMISUSEDTOCALCULATETHECVOLUMEOFCYLINDERCALMS------HALFOFTHERATIOOFSTROKEDISTANCETOTHECLENGTHOFCONNECTINGROD,REAL,>0CVH--------CYLINDERSWEPTVOLUMECFUNCTIONVOLU(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---------------VOLU()-11C------------------------BURNT()--------------------------------CCBURNT:THISPROGRAMISUSEDFORCALCULATIONOFCOMBUSTIONCHEATRELEASEVALUESANDVALUEOFINJECTEDFUELCDQB-------HEATRELEASEVALUSESbyWEIBEfunctionINCOMBUSIONPRECESSCSUBROUTINEBURNT(FI)COMMON/D/AN1,HU,AAM,ETUCOMMON/F/BEGFI,FIVB,DFI,ENDFICOMMON/G/AM,AML,AMB0,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----------------------12C------------------------ALANMD()-------------------------------CCALANMD:THISPROGRAMISUSEDFORCALCULATIONOFEXCESSCAIRCOEFFICIENTANDQUANTITYOFGASCCONTENT:FUNCTIONVOLU(FI)CSUBROUTINEBURNT(FI)CSUBROUTINEALANMD(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/TASELECTCASE(KKK)CASE(-1)ALM=1.0E4AM=AMLCASE(0)CALLBURNT(FI)ALM=AML/(AL0*AMB)AM=AML+AMBCASE(1)ALM=AML/(AL0*AMB0)AM=AML+AMB0CASE(2)AML=AM-AMBALM=AML/(AL0*AMB)CASEDEFAULTWRITE(*,*)"TheworkingprocessINDEXdoesn'tmatch!"ENDSELECTRETURNEND程序说明计算过量空气系数C----------------------13C------------------------HECOEF()-------------------------------CCHECOEF:THISPROGRAMISUSEDFORCALCULATIONOFHEATTRANSFERCCOEFFICIENTCYLINDERCWSCHINI'SFROMULACCONTENT:CFUNCTIONVOLU(FI)CFUNCTIONPRES(FI,T)CSUBROUTINEALANMD(FI)CFUNCTIONHECOEF(FI,T)COMMON/A/AR,AL,D,EPS,AN,SCOMMON/B/R,PA,TACOMMON/D/AN1,HU,AAM,ETUCOMMON/E/KKK,JJJ,ALMCOMMON/F/BEGFI,FIVB,DFI,ENDFICOMMON/G/AM,AML,AMB0,AL0,AMBCOMMON/H/RAD,PICALLALANMD(FI)VH=PI/4.*D*D*SCM=AN*S/30.C=0.

程序说明C----------------------14程序说明SELECTCASE(KKK)CASE(-1)C1=2.28CASE(0,1)C2=3.24E-3P0=(VOLU(BEGFI)/VOLU(FI))**AN1*PAC=C2*VH/(AM*R)*(PRES(FI,T)-P0)CASE(2)C1=6.18CASEDEFAULTWRITE(*,*)"TheworkingprocessINDEXdoesn'tmatch!"ENDSELECTC=(C1*CM+C)**0.8P=PRES(FI,T)HECOEF=820./T**0.53*PRES(FI,T)**0.8/D**0.2*CRETURNEND用绝热过程代替倒拖压力程序说明SELECTCASE(KKK)用绝热过程代替倒拖压15C------------------------HETRAN()-------------------------------CCHETRAN:THISPROGRAMISUSEDFORCALCUATIONOFVALUESOFHEATCTRANFERINCYLINDERVARYWITHCRANKANGLECFUNCTIONHETRAN(FI,T)COMMON/A/AR,AL,D,EPS,AN,SCOMMON/C/TW1,TW2,TW3COMMON/H/RAD,PIAA1=PI*D*D/4.AA2=AA1AA3=4.*VOLU(FI)/DHETRAN=HECOEF(FI,T)*((TW1-T)*AA1+(TW2-T)*AA2+(TW3-T)*AA3)/6./ANRETURNENDC----------------------16C------------------------PROPER()-------------------------------CCPROPER:THISPROGRAMISUSEDFORCALCULATIONOFTHERMALCPROPERTYOFGASCWRITTENBY:DIVISIONOFINTERNALCOMBUSTIONENGINEXIANJIAOTONGUNIVERSITYC SUBROUTINEPROPER(FI,T,CV,U,DULM) COMMON/E/KKK,JJJ,ALM A1=0.0975 A2=0.0485 A0=0.14455 B1=7.768 B2=3.36 C1=489.6 C2=46.4 D1=1356.8 T1=(T-273)*1.E-2 CALLALANMD(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）等容比热解析式（见第20页2.3.17）C----------------------17C-------------------------RUGKUT()------------------------------SUBROUTINERUGKUT(XVALUE,YVALUE,XSTEP,FUN)EXTERNALFUNDIMENSIONXVALUE(2),YVALUE(2),