版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
燃燒放熱規(guī)律計(jì)算程序說明燃燒放熱規(guī)律計(jì)算程序說明1任務(wù)說明《內(nèi)燃機(jī)工作過程數(shù)值計(jì)算》教學(xué)程序包括“輸入?yún)?shù).for”給出了發(fā)動(dòng)機(jī)運(yùn)轉(zhuǎn)及結(jié)構(gòu)參數(shù)和參數(shù)賦初值,“零維程序子程序.for”給出計(jì)算所需子程序。
任務(wù):編寫主程序,計(jì)算給定直噴式柴油機(jī)在壓縮、燃燒和膨脹過程中缸內(nèi)狀態(tài)參數(shù)隨曲軸轉(zhuǎn)角的變化(計(jì)算始點(diǎn):進(jìn)氣門關(guān)閉216CA,計(jì)算終點(diǎn):排氣門打開:504CA)。報(bào)告格式:(1)數(shù)學(xué)模型的建立(2)程序流程圖(3)程序清單(4)計(jì)算結(jié)果分析任務(wù)說明《內(nèi)燃機(jī)工作過程數(shù)值計(jì)算》教學(xué)程序包括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編程時(shí)必須加上程序說明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程序說明計(jì)算過量空氣系數(shù)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比內(nèi)能解析式(見第20頁2.3.16)等容比熱解析式(見第20頁2.3.17)C----------------------17C-------------------------RUGKUT()------------------------------SUBROUTINERUGKUT(XVALUE,YVALUE,XSTEP,FUN)EXTERNALFUNDIMENSIONXVALUE(2),YVALUE(2),
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2024外教聘用合同中的工作環(huán)境與安全保障措施3篇
- 《場(chǎng)效應(yīng)管講解》課件
- 2024年度冷鏈物流用地土地使用權(quán)永久轉(zhuǎn)讓與冷鏈服務(wù)合同3篇
- 2024年招投標(biāo)管理創(chuàng)新與實(shí)踐3篇
- 2024年塔吊司機(jī)聘用合同及應(yīng)急響應(yīng)預(yù)案3篇
- 2024年二零二四年度綠色食品代加工保密及質(zhì)量監(jiān)管合作協(xié)議3篇
- 2024年標(biāo)準(zhǔn)舊車交易協(xié)議版B版
- 2025車輛抵押借款合同范本
- 2024年標(biāo)準(zhǔn)游泳館租賃協(xié)議樣稿
- 2024年新一代人工智能軟件銷售及深度定制化服務(wù)合同3篇
- 傳染病護(hù)理學(xué)案例分析報(bào)告總結(jié)
- MOOC 醫(yī)學(xué)基礎(chǔ)-武漢理工大學(xué) 中國大學(xué)慕課答案
- 2024年中考作文十二大高頻熱點(diǎn)主題10-責(zé)任擔(dān)當(dāng)(素材)
- 燃煤采暖鍋爐房煙氣除塵系統(tǒng)設(shè)計(jì)說明書樣本
- 主動(dòng)脈瓣置換手術(shù)
- 2024年四平職業(yè)大學(xué)單招職業(yè)技能測(cè)試題庫及答案解析
- 女生穿搭技巧智慧樹知到期末考試答案2024年
- 作業(yè)成本法在F公司的應(yīng)用研究
- MOOC 大數(shù)據(jù)與法律檢索-湖南師范大學(xué) 中國大學(xué)慕課答案
- 活動(dòng)現(xiàn)場(chǎng)醫(yī)療應(yīng)急方案
- 宏觀經(jīng)濟(jì)學(xué)(山東聯(lián)盟-山東財(cái)經(jīng)大學(xué))智慧樹知到期末考試答案2024年
評(píng)論
0/150
提交評(píng)論