內(nèi)燃機熱力計算怎么用C語言編程,？模板如下,，求助,。

方進小侯爺 · 發(fā)表于 2015-6-28 22:58:34

#include "math.h"
#include "stdio.h"
#define PI 3.1415926

double D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
/*缸徑,、行程,、連桿長、比油耗,、轉(zhuǎn)速,、功率、氣缸數(shù),、壓縮比,、低熱值、理論空燃比,、Weibe函數(shù)常數(shù)m,、氣體常數(shù)、進氣門關閉時壓力和溫度 */
double m_T[720],m_P[720],m_a[720];
/*氣缸溫度,、壓力,、瞬時過量空氣系數(shù) */
int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
double Cv();
double dV_dCA();
double dQw_dCA();
double dQB_dCA();
double U();
double V_CA();
double dU_dnmd();
double dT_Ys_dCA();
double dT_Pz_dCA();
double dT_Rs_dCA();
double dT_Hq_dCA();
double dT_dCA(double an);
void Calculate();

/* A編寫各子函數(shù)程序*/
double Cv() /*求定容比熱*/
{

}

double  dV_dCA() /*求dv/dφ*/
{

}

double V_CA() /*求氣缸容積*/
{

}

double dQw_dCA() /*求dQw/dφ*/
{

}

double dQB_dCA() /*求dQB/dφ*/
{

}

double U() /*求u*/
{

}

double dU_dnmd()  /*求du/dλ*/
{

}

double dT_Ys_dCA() /*壓縮期*/
{

}

double dT_Pz_dCA() /*膨脹期*/
{

}

double dT_Rs_dCA() /*燃燒期*/
{

}

double dT_dCA(double an) /*求dT/dφ*/
{
CA=an*PI/180;
V=V_CA();
if(A_SC<an && an<A_CB) return dT_Ys_dCA();
if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
if(A_EO<an && an<=A_SO) /*換氣過程不計算,，線性賦值，避免零值 */
P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
if(A_SO<an || an<A_EC)
P=110000-20000*(an>A_SO?an-A_SO:an-A_SO+720)/(A_EC-A_SO+720);
if(A_EC<=an && an<=A_SC) P=90000;
return 0;
}

void Calculate(void)
{
double Tpre,K1,K2,an;
int index;
Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效壓力bar*/
mBo=Ne*ge/(30000000*n*gnum); /*單缸循環(huán)供油量kg */
CA=A_SC*PI/180; /*進氣關閉時的角度轉(zhuǎn)弧度*/
V=V_CA();/*進氣關閉時角度下的氣缸容積*/
P=P_SC;
T=T_SC;
index=(int)A_SC;
m_T[index]=T;
m_P[index]=P;
m=P*V/(R*T);/*進氣關閉時角度下的氣缸內(nèi)工質(zhì)質(zhì)量*/
mL=m-0.01*mBo;/*進氣關閉時角度下的氣缸內(nèi)空氣質(zhì)量*/
mf=m-mL;
Tpre=T;
an=A_SC+1;

do{

  /*B 編寫歐拉法求解常微分方程程序 */

}while(an!=A_SC+1);
}

void main(void)
{
int i;
FILE *fp;
D=0.11; /*發(fā)動機數(shù)據(jù)按各自的任務書輸入 */
S=0.13;
L=0.21;
ge=236.0;
n=2400;
Ne=101568;
gnum=6;  ysb=17; M=0.7;
Hu=41868000;
R=287.08;
Pa=101000; /*大氣壓力*/
Ta=300; /*大氣溫度*/
hs=(735.0+R)*Ta;  /*大氣焓值*/
P_SC=90000;
T_SC=313;

A_SC=216; /*進氣門關閉*/
A_CB=350; /*燃燒開始*/
A_CE=410; /*燃燒結(jié)束*/
A_EO=484; /*排氣門開啟*/
A_SO=705; /*進氣門開啟*/
A_EC=15;  /*排氣門關閉*/

Lo=14.4;

for(i=0;i<720;i++)
{
m_P[i]=0;
m_CA[i]=i;
m_T[i]=0;
}

Calculate();

/*C編寫數(shù)據(jù)輸出程序,，計算結(jié)束后按順序輸出0-720轉(zhuǎn)角的“缸壓,，溫度”數(shù)據(jù)到outdata.txt文件 */

}

祥云飛魚 · 發(fā)表于 2015-6-28 23:45:10

不難啊完全是把公司翻譯成 C 又沒有用指針跳轉(zhuǎn)的,。,。。

一點一點來 · 發(fā)表于 2015-6-29 01:45:20

這是作業(yè)的干活吧……

		自動登錄	找回密碼
密碼			注冊會員

內(nèi)燃機熱力計算怎么用C語言編程,？模板如下,，求助,。

相關帖子