燃燒放熱規(guī)律計算程序說明燃燒放熱規(guī)律計算程序說明1任務說明《內燃機工作過程數值計算》教學程序包括“輸入參數.for”給出了發(fā)動機運轉及結構參數和參數賦初值，“零維程序子程序.for”給出計算所需子程序。

任務：編寫主程序，計算給定直噴式柴油機在壓縮、燃燒和膨脹過程中缸內狀態(tài)參數隨曲軸轉角的變化（計算始點：進氣門關閉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),