371型車用汽油機(jī)的熱力學(xué)、動(dòng)力學(xué)設(shè)計(jì)計(jì)算與分析

華中科技大學(xué)文華學(xué)院畢業(yè)設(shè)計(jì)論文題目：371型車用汽油機(jī)的熱力學(xué)、動(dòng)力學(xué)設(shè)計(jì)計(jì)算與分析目錄摘要····················································2關(guān)鍵詞··················································2Abstract················································3Keywords···············································3前言····················································41．本文主要研究?jī)?nèi)容與方法································51.1研究?jī)?nèi)容·········································51.2研究方法·········································52.熱力學(xué)、動(dòng)力學(xué)設(shè)計(jì)計(jì)算································52.1內(nèi)燃機(jī)的熱力學(xué)計(jì)算·······························52.2內(nèi)燃機(jī)的運(yùn)動(dòng)學(xué)計(jì)算·······························62.3內(nèi)燃機(jī)的動(dòng)力學(xué)計(jì)算······························102.4計(jì)算結(jié)果分析····································142.5改進(jìn)方案········································15參考文獻(xiàn)···············································16附錄··················································17致謝··················································28371型車用汽油機(jī)的熱力學(xué)、動(dòng)力學(xué)設(shè)計(jì)計(jì)算與分析摘要本文簡(jiǎn)要介紹了汽油機(jī)的熱力學(xué)、動(dòng)力學(xué)設(shè)計(jì)計(jì)算與分析，根據(jù)內(nèi)燃機(jī)原理和內(nèi)燃機(jī)機(jī)構(gòu)運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)知識(shí)，建立設(shè)計(jì)計(jì)算數(shù)學(xué)模型。并且完成全部設(shè)計(jì)計(jì)算，繪制發(fā)動(dòng)機(jī)的熱力學(xué)、運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)曲線；并對(duì)結(jié)果進(jìn)行實(shí)事求是的分析、評(píng)估和提出改進(jìn)方案。關(guān)鍵詞：內(nèi)燃機(jī)；熱力學(xué)；動(dòng)力學(xué)；運(yùn)動(dòng)學(xué)；設(shè)計(jì)計(jì)算與分析。Type371gasolineengineforvehicledynamicsandthermodynamics,thedesigncalculationandanalysisAbstractThispaperbrieflyintroducesthegasolineenginethermodynamic,kineticdesignandanalysis,accordingtotheprinciplesofinternalcombustionengineandinternalcombustionenginemechanismkinematics,dynamics,establishesthemathematicalmodelfordesignandcalculation.Andcompletethedesigncalculation,renderingenginethermodynamics,dynamicsandkinematicscurves;andtheresultsofseektruthfromfactsanalysis,assessmentandputforwardtheimprovementscheme.KeyWords:Internalcombustionengine;thermodynamics;kinetics;kinematics;designcalculationandanalysis.前言隨著人類環(huán)保意識(shí)的不斷增強(qiáng),對(duì)車輛動(dòng)力的排放要求也越來(lái)越嚴(yán)格。這種發(fā)展趨勢(shì)勢(shì)必要求發(fā)動(dòng)機(jī)研制、生產(chǎn)單位開發(fā)出排放更清潔、效率更高的動(dòng)力裝置,才能迎接未來(lái)更加激烈的市場(chǎng)競(jìng)爭(zhēng)。綜合評(píng)述了發(fā)動(dòng)機(jī)系統(tǒng)的熱力學(xué)模擬發(fā)展和現(xiàn)狀。在熱力學(xué)第一定律和第二定律的基礎(chǔ)上,展示了國(guó)外在增壓和非增壓汽油機(jī)、帶有催化轉(zhuǎn)換器的汽油機(jī)以及渦輪增壓柴油機(jī)在穩(wěn)態(tài)和瞬態(tài)工況下所建立的分析模型,并介紹了研究發(fā)動(dòng)機(jī)氣缸缸內(nèi)調(diào)諧和暖機(jī)特性的模擬模型，但是,實(shí)現(xiàn)上述目標(biāo)需要在發(fā)動(dòng)機(jī)的整個(gè)開發(fā)過程中實(shí)行多學(xué)科的共同努力。所以選擇對(duì)車用汽油機(jī)的熱力學(xué)計(jì)算分析后對(duì)汽油機(jī)的經(jīng)濟(jì)效益以及環(huán)保問題有益。由于近年來(lái)歐美遇到了經(jīng)濟(jì)危機(jī)，汽車工業(yè)發(fā)展逐漸衰退；而我國(guó)在新世紀(jì)進(jìn)行了產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，人民生活水平不斷提高，購(gòu)買力增強(qiáng)，對(duì)汽車的需求明顯提高，我國(guó)汽車工業(yè)迎來(lái)了一個(gè)發(fā)展的春天。工信部裝備工業(yè)司預(yù)計(jì)，到2020年我國(guó)汽車保有量將超過2億輛，這意味著屆時(shí)我國(guó)將真正成為汽車超級(jí)大國(guó)。綜上所述，針對(duì)于我國(guó)汽車工業(yè)的發(fā)展環(huán)境，當(dāng)前和今后很長(zhǎng)一段時(shí)間內(nèi)，提高車用汽油機(jī)的性能，改善燃油經(jīng)濟(jì)性，降低汽車所引起的污染物和溫室氣體排放是解決當(dāng)前所有問題的關(guān)鍵。而要實(shí)現(xiàn)這幾點(diǎn)，就需要找到適合于我國(guó)具體國(guó)情的車用汽油機(jī)產(chǎn)品類型。我國(guó)的車用汽油機(jī)事業(yè)值得期待。主要研究方案：1.根據(jù)內(nèi)燃機(jī)原理和內(nèi)燃機(jī)機(jī)構(gòu)運(yùn)動(dòng)學(xué)，動(dòng)力學(xué)知識(shí)，建立設(shè)計(jì)計(jì)算數(shù)學(xué)模型；2.運(yùn)用C預(yù)言編制設(shè)計(jì)計(jì)算程序，然后，根據(jù)計(jì)算結(jié)果對(duì)數(shù)學(xué)模型和計(jì)算程序進(jìn)行多次調(diào)試，光順，以求得出接近實(shí)際情況的曲線形態(tài)；3.在此基礎(chǔ)上，確定計(jì)算程序是適用的，進(jìn)而完成全部設(shè)計(jì)計(jì)算；4.選用適用，有效的數(shù)據(jù)處理工具，繪制發(fā)動(dòng)機(jī)的熱力學(xué)，運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)曲線；5.最后，輸出設(shè)計(jì)計(jì)算的結(jié)果，并進(jìn)行實(shí)事求是的分析，評(píng)估和提出改進(jìn)方案；6.用“AutoCAD”繪制該機(jī)型的《活塞—連桿—曲柄機(jī)構(gòu)》的部件工程圖。1本文主要研究?jī)?nèi)容及研究方法1.1研究?jī)?nèi)容根據(jù)內(nèi)燃機(jī)原理和內(nèi)燃機(jī)機(jī)構(gòu)運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)知識(shí)、建立設(shè)計(jì)計(jì)算數(shù)學(xué)模型，以便計(jì)算與分析。1.2研究方法運(yùn)用C語(yǔ)言編制設(shè)計(jì)計(jì)算程序，然后，根據(jù)計(jì)算結(jié)果對(duì)數(shù)學(xué)模型和計(jì)算程序進(jìn)行多次調(diào)試，光順，以求得出接近實(shí)際情況的曲線形態(tài)；在此基礎(chǔ)上，確定計(jì)算程序是適用的，進(jìn)而完成全部設(shè)計(jì)計(jì)算；選用適用，有效的數(shù)據(jù)處理工具，繪制發(fā)動(dòng)機(jī)的熱力學(xué)，運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)曲線；最后，輸出設(shè)計(jì)計(jì)算的結(jié)果，并進(jìn)行實(shí)事求是的分析，評(píng)估和提出改進(jìn)方案；用“AutoCAD”繪制該機(jī)型的《活塞—連桿—曲柄機(jī)構(gòu)》的部件工程圖。2熱力學(xué)、動(dòng)力學(xué)設(shè)計(jì)計(jì)算2.1內(nèi)燃機(jī)的熱力學(xué)計(jì)算所介紹的熱力學(xué)計(jì)算，其考慮的影響因素是最基本的，沒有AVL的BOOST那么全面，其功能非常有限。因此，只能適于對(duì)主要技術(shù)指標(biāo)及相關(guān)參數(shù)作粗略的估算，運(yùn)行程序后的結(jié)果如下：TheCalculationResultsof371CMforNormalEngine~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~PetrolEngineType:SQR371CMRatedPower:51.0kW/6000r/minBoreDiameter:71mm,Crankthrow:41.95mm,Displacement:0.9965LED=0.008mm,SS=83.9721mm,TD=2.506deg,BD=184.633degCP=1.280,HP=6.748,nc=1.379,ne=1.096Pa=0.096MPa,Pc=1.880MPa,Pz=5.828MPa,Pb=0.718MPaTa=325.6K,Tc=737.9K,Tz=2749.3K,Tb=2312.2KVs=0.9974L,El=10.500,Ee=8.641,AF=0.8458Mi=109.266N.m,Me=81.290N.m,Nm=17.58kW,Ne=51.08kWhm=0.7440ge=265.9g/kW.h,Pe=1.024MPa,he=0.3080注意：由于氣缸中心線正向偏置8毫米，導(dǎo)致活塞的位移曲線不對(duì)稱，且最大位移為84.0485mm、并出現(xiàn)在下止點(diǎn)前2°CA處；實(shí)際排量為0.9983L。2.2內(nèi)燃機(jī)的運(yùn)動(dòng)學(xué)計(jì)算運(yùn)動(dòng)學(xué)計(jì)算主要是求解活塞的運(yùn)動(dòng)位移、速度和加速度，為熱力學(xué)、動(dòng)力學(xué)和性能指標(biāo)的計(jì)算奠定基礎(chǔ)。從圖中幾何關(guān)系可以得出：影響內(nèi)燃機(jī)結(jié)構(gòu)的一個(gè)重要特性參數(shù)。現(xiàn)代內(nèi)燃機(jī)的連桿比值一般在(即0.26一0.31)之間。采用較大連桿比優(yōu)點(diǎn)：可使機(jī)體高度減小，重量減輕，缺點(diǎn)：使活塞加速度和連桿擺角有所增大，慣性力增加2.2.1活塞的位移將上式按二項(xiàng)式定理展開后，忽略高階項(xiàng)可得近似公式活塞的位移隨連桿比的變化2.2.2活塞的速度整理，得：活塞的速度活塞的速度在曲柄旋轉(zhuǎn)一周中，時(shí)快時(shí)慢地變化著，他的平均速度可以表示為：Cm=sn/30(m/s)式中，s—活塞行程，單位為m；n—曲軸轉(zhuǎn)速，單位為r/min?；钊骄俣菴m雖然只能粗略地估計(jì)活塞運(yùn)動(dòng)的快慢，但是它是表征內(nèi)燃機(jī)性能指標(biāo)的重要參數(shù)之一。2.2.3活塞的加速度λ值越大，活塞速度最大值也越大，相應(yīng)達(dá)到活塞最大值的曲柄轉(zhuǎn)角α越小。其三種圖形均不對(duì)稱。2.3內(nèi)燃機(jī)的動(dòng)力學(xué)計(jì)算在運(yùn)動(dòng)學(xué)、熱力學(xué)計(jì)算的基礎(chǔ)上，進(jìn)行動(dòng)力學(xué)計(jì)算。2.3.1動(dòng)力計(jì)算是分析發(fā)動(dòng)機(jī)各主要部件的受力情況，為強(qiáng)度計(jì)算提供數(shù)據(jù)。其內(nèi)容如下：1、曲柄連桿機(jī)構(gòu)運(yùn)動(dòng)質(zhì)量的確定初步設(shè)計(jì)時(shí)常用經(jīng)驗(yàn)數(shù)據(jù)來(lái)粗略估算。連桿的質(zhì)量還需要通過連桿重心的確定，再用雙質(zhì)量系統(tǒng)替代之。從而求出往復(fù)質(zhì)量和迥轉(zhuǎn)不平衡質(zhì)量。2、P-φ示功圖的求取將熱計(jì)算得的P-V示功圖，用發(fā)動(dòng)機(jī)運(yùn)動(dòng)學(xué)公式將其展開，即得P-φ示功圖。將活塞的位移轉(zhuǎn)換成對(duì)應(yīng)的曲柄轉(zhuǎn)角，以α代表曲柄轉(zhuǎn)角，由此可得各曲柄轉(zhuǎn)角α下的氣體力值Pg（α），單位為MPa。3、往復(fù)慣性力Pj（α）計(jì)算（MPa）或（MPa）式中：mj—往復(fù)運(yùn)動(dòng)質(zhì)量，kg；R—曲柄半徑，mm；D—?dú)飧字睆?，mm；ω—曲軸旋轉(zhuǎn)角速度，rad/s；β—連桿擺角，rad。4、總作用力P（α）式中：Pβ—活塞底部氣體壓力，。5、活塞側(cè)推力PH（α）（MPa）6、連桿力PC（α）（MPa）7、法向力PN（α）（MPa）8、切向力PT（α）（MPa）9、總切向力對(duì)四沖程曲柄均勻排列的情況：（MPa）式中：z—?dú)飧讛?shù)。最后求出平均切力(Σpr)cr，其單位為：MPa。10、曲柄銷負(fù)荷RB（α）計(jì)算水平分量RΒh(α)=PT(α)(MPa)垂直分量RBv(α)=PN(α)-Prβ(MPa)合力(MPa)式中：(MPa)mB—連桿組則算到大端的質(zhì)量，kg。在完成動(dòng)力計(jì)算后，應(yīng)繪出Pg(α)、Pj(α)、PH(α)、PN(α)、PT(α)、ΣPT(α)、RB(α)各曲線圖。最后按總切力曲線作準(zhǔn)確性校核：由總切力得到的指示功率：（kW）式中，—平均切力，MPa；Fp—活塞面積，mm2；R—曲柄半徑，mm；ω—曲軸角速度，rad/s，；n—曲軸轉(zhuǎn)速度，r/min。將Ni與熱計(jì)算得的Ni比較：誤差：△應(yīng)小于5%2.3.2作用于曲柄連桿機(jī)構(gòu)中的力和力矩2.3.3活塞與活塞銷的相互作用力（Mi—109.266N.m指平均指示扭矩，波浪狀為瞬時(shí)指示扭矩）2.3.4活塞與連桿小頭\連桿大頭與曲柄銷之間的相互作用力2.4計(jì)算結(jié)果分析由以上結(jié)果數(shù)據(jù)可知，基于當(dāng)今汽油機(jī)的燃料燃燒技術(shù)，燃燒效率ηc可達(dá)98％左右。以SQR371為例，若其有效效率ηe＝0.3081、機(jī)械效率ηm＝0.7440。根據(jù)以上分析，即使SQR371的循環(huán)熱效率ηt為0.4224，其有效效率（總效率）僅為0.3080。因此，在平常的口頭交流中，不能將熱效率ηt與有效效率（總效率）ηe混為一談，以免引起誤解。2.5改進(jìn)方案在發(fā)動(dòng)機(jī)的機(jī)械損失中，活塞組的摩擦損失占據(jù)40～50%。當(dāng)發(fā)動(dòng)機(jī)的機(jī)械效率為50%時(shí)，若活塞組的摩擦損失減少25%，燃油消耗率則可以降低9%。常用的低轉(zhuǎn)速工況下，其機(jī)械損失扭矩可達(dá)全負(fù)荷扭矩的10％左右；在高速工況下，則可達(dá)最大扭矩的20％以上。發(fā)動(dòng)機(jī)的機(jī)械損失隨著轉(zhuǎn)速的提高而增加；而負(fù)荷的變化對(duì)機(jī)械損失影響不大。因此，發(fā)動(dòng)機(jī)的負(fù)荷率越低，其機(jī)械效率相對(duì)也越低。進(jìn)氣閥關(guān)幾何壓縮比實(shí)際壓縮比機(jī)械效率有效扭矩標(biāo)定功率耗油率指示壓力有效效率℃AN.MkWg/Kw.hMPa23910.58.6340.744081.29051.08265.91.0240.308025010.57.9210.750284.26552.95235.31.0620.3480根據(jù)內(nèi)燃機(jī)原理的基本公式，由熱力學(xué)、動(dòng)力學(xué)的計(jì)算結(jié)果，很容易算出標(biāo)定工況下的功率、扭矩、燃油消耗率、平均有效壓力及發(fā)動(dòng)機(jī)的總效率（有效效率）ηe等主要指標(biāo)。眾所周知：配氣凸輪的型線及配氣相位對(duì)發(fā)動(dòng)機(jī)的性能指標(biāo)影響很大。如前所述，實(shí)際壓縮比直接與進(jìn)氣終點(diǎn)所對(duì)應(yīng)的活塞位置（構(gòu)成的氣缸容積Va）有關(guān)。若SQR371按米勒循環(huán)的做法、使進(jìn)氣終點(diǎn)設(shè)定在下止點(diǎn)后70（原為59）°CA，當(dāng)不改變其它條件時(shí)，其主要性能指標(biāo)的變化如下表所列：可見，若進(jìn)氣閥比原設(shè)計(jì)推遲11℃A關(guān)閉，可是發(fā)動(dòng)機(jī)的有效效率提高四個(gè)百分點(diǎn)。參考文獻(xiàn)1.周建立《發(fā)動(dòng)機(jī)原理實(shí)驗(yàn)教程》中國(guó)電力出版社出版社2.周建立《發(fā)動(dòng)機(jī)原理實(shí)驗(yàn)教程》中國(guó)電力出版社出版社3.楊有安等.《C語(yǔ)言》.人民郵電出版社.2009年2月第1版4.朱建元等《船舶內(nèi)燃機(jī)》.人民交通出版社.2008年9月第1版5.王元慶TCR發(fā)動(dòng)機(jī)運(yùn)動(dòng)學(xué)、熱力學(xué)和動(dòng)力學(xué)的分析計(jì)算華中科技大學(xué)6.長(zhǎng)春汽車研究所《汽車實(shí)驗(yàn)技術(shù)手冊(cè)》吉林科學(xué)出版社7.吳克剛《發(fā)動(dòng)機(jī)測(cè)試技術(shù)》人民交通出版社附錄一：《活塞—連桿—曲柄機(jī)構(gòu)》的部件工程圖曲軸部件圖活塞部件圖附錄二：C語(yǔ)言程序如下：/*TheThermodynamicsandMechanicsCalculationfortheEccentricI.C.Engine*/#include"math.h"#include"conio.h"#include"stdio.h"#include"graphics.h"FILE*c71,*ydx,*dlx,*mbc,*zcs;main(){inti,jz,ig,nv,ec,ic,kr,kz,kp,eo,io,ib;floatpe,ge,ee,ui,ue,tq,bt,tne;floatpy=3.14159265,rd=.01745329,hs[721],hv[721],ha[721],pg[721],pj[721];floatpp[721],pd[721],ph[721],po[721],prn[721],pt[721],rph[721],rpv[721];floatrkh[721],rkv[721],pob[721],ptb[721],tqi[721],nhi[721],tb,pb,ss,net;floatcm=.87,hm=.126,om=.004,hl=43961.,zc=1.11,yn=1.36,fi=.9462,zt=2108.;floatez=.7,pn=1.28,gr=.02,eb=.85,db,rs,cr,ed,rl,el,wp,wc,ws,wk,ej,pz,pk;floattk,dp,dt,tr,af,em,otq,oei,one,pf,snf,fh,tc,pc,gl,gs,bl,bs,xz,bz,br;floatog,pl,td,bd,rt,sa,sb,sc,pa,ta,es,ev,cn,ys,bo,cb,wj,wl,fp,prl,ei,gi;floatsalf,calf,sbet,cbet,be,tqt,tqm,hit,snh,nht,pr,ri,cz,tz,hu,cp,hp,en;floatpit,pet,git,get,het,hmt,prk,pi,tqe,nm,vs;c71=fopen("c:\\my_files\\c71.inp","r");fscanf(c71,"%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%f,%f,%f,%f,%f,%f," "%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f", &jz,&ig,&nv,&ec,&ic,&kr,&kz,&kp,&eo,&io,&db,&rs,&cr,&wp,&wc, &ws,&wk,&ej,&pz,&pk,&tk,&dp,&dt,&tr,&af,&pf,&snf,&fh,&ed);fclose(c71);rl=rs/cr;el=ed/cr;og=nv*py/30;td=atan(el/sqrt((1+rl)*(1+rl)-el*el))/rd;/*"+"reversal*/bd=atan(el/sqrt((1-rl)*(1-rl)-el*el))/rd+180;/*"+"reversal*/bt=bd-td;rt=rd*td;ydx=fopen("c:\\my_files\\ydx.c71","w");/*c-e*/for(i=0;i<721;i++)/*721*/{ri=rd*i;salf=sin(ri-rt);calf=cos(ri-rt);sbet=rl*salf+el;cbet=sqrt(1-sbet*sbet);hs[i]=1e3*(cr*(sqrt((1+rl)*(1+rl)-el*el)-cbet)-rs*calf);if(hs[i]>hs[i-1])ss=hs[i];hv[i]=rs*og*(salf+(rl*salf+el)*calf/cbet);ha[i]=rs*og*og*(calf+(rl*(calf*calf-salf*salf)-el*salf)/cbet+(rl* salf*calf+el*calf)*(rl*rl*salf*calf+el*rl*calf)/pow(cbet,3));fprintf(ydx,"%4d%9.4f%8.4f%9.2f\n",i,hs[i],hv[i],ha[i]);printf("%4d%9.4f%8.4f%9.2f\n",i,hs[i],hv[i],ha[i]);}fclose(ydx);ib=(int)bt;if((bt-ib)>=.5)ib+=1;if(hs[ib]!=ss)ss=ss;sa=hs[ic];sb=hs[eo];sc=ss/(ej-1);es=(sa+sc)/sc;vs=db*db*py*ig*ss/4;for(i=0;i<=ic;i++){pg[i]=pk-dp*(i+720-io)/(ic+720-io);}pa=pk-dp;ta=(tk+dt+zc*gr*tr)/(1+gr);for(cn=1+8.315/(19.259+.0025*ta*(1+pow(es,(yn-1))));fabs(cn-yn)/cn>.001;yn=.5*(cn+yn),cn=1+8.315/(19.259+.0025*ta*(1+pow(es,(yn-1)))));tc=ta*pow(es,(cn-1));pc=pa*pow(es,cn);gl=(cm/12+hm/4-om/32)/.21;gs=af*gl;bl=(gs+hm/4+om/32)/(gs+(1./286.));bs=(bl+gr)/(1+gr);xz=ez/eb;bz=1+(bl-1)*xz/(1+gr);for(i=ic+1;i<=kr;i++){pg[i]=pa*pow((hs[ic]+sc)/(hs[i]+sc),cn);}ys=pz/pc;bo=1.67*tr;key_1:cb=((1.064+gr)*(20.47+.0036*bo)+(af-1)*(1-gr)*(19.259+.0025*bo))/(af*(1+gr)+.064);hu=hl+gs*(1+gr)*(bs*cb-(19.259+.0025*bo))*tk;key_2:cz=((1.064*xz+gr)*(20.47+0.0036*zt)+(af*(1+gr)-(xz+gr))*(19.259+.0025*zt))/(af*(1+gr)+.064*xz);tz=(ez*hu/gs+(19.259+.0025*zt+8.315*ys)*tc+gr*(cb+8.315*ys)*tc)/bz/(1+gr)/(cz+8.315);zt=.5*(tz+zt);if(fabs(tz-zt)>.5){zt=.5*(tz+zt);gotokey_2;}cp=bs*tz/ys/tc;hp=es/cp;key_3:br=tz*bs/bz/pow(hp,(pn-1));if(fabs(br-bo)>1.){bo=.5*(bo+br);gotokey_1;}key_4:en=1+6.528*(bz*tz/(bz+.0002)-br)/(hu*(eb-ez)/(gs*(1+gr)*bs)+bz*cz*tz/(bz+.0002)-cb*br);if((fabs(en-pn)/pn)>.001){pn=.5*(en+pn);gotokey_3;}tb=tz*bs/bz/pow(hp,(en-1));if(fabs(tb-br)>1.){br=.5*(br+tb);gotokey_4;}pb=pz/pow(hp,en);for(i=kr+1;i<=kz;i++){pg[i]=pz-.00176*pow((kz-i),2.598);}for(i=kz;i<=kp+17;i++){pg[i]=pg[kz]-.00118*pow((i-kz),2.356);}for(i=kp+18;i<=kp+25;i++){pg[i]=pg[396]*pow((hs[396]+sc)/(hs[i]+sc+.218),(1.528));}for(i=kp+26;i<=eo;i++){pg[i]=pb*pow((hs[eo]+sc)/(hs[i]+sc),en);}for(i=(eo+1);i<=(eo+22);i++){pg[i]=pg[eo]-pow((i-eo),1.48)*(pg[eo]-pa)/(io-eo);}for(i=(eo+23);i<=(eo+48);i++){pg[i]=pg[eo+22]-.0098*pow((i-eo-22),.668);}for(i=(eo+49);i<=io;i++){pg[i]=pg[eo+48]-(i-eo-48)*(pg[eo+48]-pa)/(io-eo-48);}for(i=(io+1);i<721;i++){pg[i]=pk-dp*(i-io)/(ic+720-io);}pl=pc*(ys*(cp-1)+cp*ys*(1-1/pow(hp,(en-1)))/(en-1)-(1-1/pow(es,(cn-1)))/(cn-1))/(es-1);pi=fi*(pl*(es-1)+(pb-pa)*(ej-es))/(ej-1);ev=es*tk*pa/(ej-1)/ta/pk/(1+gr);gi=433283.4*ev*pk/gs/tk/pi;ei=3.6e6/gi/hl;pe=pi*em;ge=gi/em;ee=ei*em;ui=pi*ig*nv*py*ss*db*db/120/jz;ue=ui*em;tq=3e4*ue/nv/py;wj=wp+ws;wl=wc-ws;fp=py*db*db/4;prl=wl*rs*og*og/fp/1e6;prk=wk*rs*og*og/fp/1e6;pr=prl+prk;dlx=fopen("c:\\my_files\\dlx.c71","w");/*c-e*/tqt=0.;nht=0.;for(i=0;i<721;i++)/*721*/{ri=rd*i;salf=sin(ri-rt);calf=cos(ri-rt);sbet=rl*salf+el;cbet=sqrt(1-sbet*sbet);be=atan(sbet/cbet);pj[i]=-wj*ha[i]/fp/1e6;pp[i]=pg[i]+pj[i]-pf;ph[i]=pp[i]*sbet/cbet;nhi[i]=1e3*fh*fp*fabs(ph[i]*hv[i]);pd[i]=pp[i]/cbet;po[i]=pd[i]*cos(ri-rt+be);pt[i]=pd[i]*sin(ri-rt+be);prn[i]=po[i]-prl;tqi[i]=pt[i]*rs*fp*1e6;rph[i]=prl*sin(ri-rt+be);rpv[i]=pd[i]-prl*cos(ri-rt+be);rkh[i]=(pd[i]*sin(be)+pr*salf)/2;rkv[i]=(pd[i]*cos(be)-pr*calf)/2;pob[i]=(po[i]-pr)/2;ptb[i]=pt[i]/2;tqt+=tqi[i];tqt=tqt-tqi[0];nht+=nhi[i];nht=nht-nhi[0];fprintf(dlx,"%4d%7.4f%8.4f%7.4f%7.4f%7.4f%7.4f%7.4f%7.4f%8.3f\n", i,pg[i],pj[i],pd[i],ph[i],rkv[i],rkh[i],pob[i],ptb[i],tqi[i]);printf("%4d%7.4f%8.4f%7.4f%7.4f%7.4f%7.4f%7.4f%7.4f%8.3f\n", i,pg[i],pj[i],pd[i],ph[i],rkv[i],rkh[i],pob[i],ptb[i],tqi[i]);}fclose(dlx);c71=fopen("e:\\my_files\\c71.tbc","w");/*c-e*/for(i=0;i<721;i++)/*721*/{fprintf(c71,"%3d%8.4f%8.3f%9.2f%7.4f%8.4f%8.4f%8.4f%8.4f\n", i,hs[i],hv[i],ha[i],pg[i],ph[i],pd[i],pob[i],ptb[i]);}fclose(c71);mbc=fopen("c:\\my_files\\mbc.btc","w");/*c-e*/for(i=0;i<721;i++)/*721*/{fprintf(mbc,"%3d%8.4f%8.4f\n",i,pob[i],ptb[i]);}fclose(mbc);tqm=tqt/720;pit=tqm*py*jz/fp/ss/1000;/*ss@mm*/git=433283.*ev*pk/gs/tk/pit;hit=3.6e6/git/hl;snh=ig*nht/720;tne=ig*py*tqm*nv/3e4;net=tne-snh-snf;hmt=net/tne;het=hit*hmt;pet=pit*hmt;get=git/hmt;tqm=ig*tqm;tqe=tqm*hmt;nm=snf+snh;zcs=fopen("c:\\my_files\\zcs.c71","w");/*c-e*/fprintf(zcs,"TheCalculationResultsof371CMforNormalEngine\n""~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\n""PetrolEngineType:SQR371CMRatedPower:51.0kW/6000r/min\n""BoreDiameter:71mm,Crankthrow:41.95mm,Displacement:0.9965L\n""\n""ED=%4.3fmm,SS=%6.4fmm,TD=%5.3fdeg,BD=%7.3fdeg\n""CP=%5.3f,HP=%5.3f,nc=%5.3f,ne=%5.3f\n""Pa=%5.3fMPa,Pc=%5.3fMPa,Pz=%5.3fMPa,Pb=%5.3fMPa\n""Ta=%5.1fK,Tc=%5.1fK,Tz=%6.1fK,Tb=%6.1fK\n""Vs=%6.4fL,El=%6.3f,Ee=%6.3f,AF=%6.4f\n""Mi=%6.3fN.m,Me=%6.3fN.m,Nm=%5.2fkW,Ne=%5.2fkW\n""hm=%6.4fge=%5.1fg/kW.h,Pe=%5.3fMPa,he=%6.4f\n""\n",1e3*ed,ss,td,bd,cp,hp,cn,en,pa,pc,pz,pb,ta,tc,tz,tb,vs,ej,es,af,tqm,tqe,nm,net,hmt,get,pet,het);fclose(zcs);//closegraph();return(0);}/*TheCalculationResultsof371CMforNormalEnginePetrolEngineT