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

[復制鏈接]

方進小侯爺

電梯直達

1^#

發表于 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;
/*缸徑、行程、連桿長、比油耗、轉速、功率、氣缸數、壓縮比、低熱值、理論空燃比、Weibe函數常數m、氣體常數、進氣門關閉時壓力和溫度 */
double m_T[720],m_P[720],m_a[720];
/*氣缸溫度、壓力、瞬時過量空氣系數 */
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編寫各子函數程序*/
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); /*單缸循環供油量kg */
CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
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);/*進氣關閉時角度下的氣缸內工質質量*/
mL=m-0.01*mBo;/*進氣關閉時角度下的氣缸內空氣質量*/
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; /*發動機數據按各自的任務書輸入 */
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; /*燃燒結束*/
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編寫數據輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數據到outdata.txt文件 */

}

內燃機, C語言

分享到: QQ好友和群 QQ空間

收藏0 轉播分享 淘帖0 問題專業，描述清楚0 伸手黨/灌水/看不懂0

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

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

相關帖子