内燃机热力计算怎么用C语言编程？模板如下，求助。

方进小侯爷

#include "math.h"
+ m" s+ h& p! t$ k- V- Q#include "stdio.h"
+ G9 b' C* n( n4 \3 j* w4 j9 l$ r#define PI 3.1415926

/ F" p) b2 }0 }* c" Ddouble D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
, }# M& D$ {" V# t- O- @  x% J/*缸径、行程、连杆长、比油耗、转速、功率、气缸数、压缩比、低热值、理论空燃比、Weibe函数常数m、气体常数、进气门关闭时压力和温度 */
        double m_T[720],m_P[720],m_a[720];
/ M- F; y# `9 U: O' x/ ^/*气缸温度、压力、瞬时过量空气系数 */
% ^7 K: z# W& k# a% q3 H        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;
4 ?/ d8 r3 \: \9 P. f        double Cv();
        double dV_dCA();
" V( n8 z: f: |5 r  n% }        double dQw_dCA();
" u; n" y& x3 a9 Y9 m' r/ }; Y        double dQB_dCA();
" E4 i) y& b) ~2 A5 m' H        double U();
        double V_CA();
* _6 V5 {+ \. T) o4 M9 c$ W- H3 S$ ]        double dU_dnmd();
        double dT_Ys_dCA();
8 F2 W* S$ o% W+ }* _: g8 i        double dT_Pz_dCA();
8 F5 k4 E/ n- f* [        double dT_Rs_dCA();
; I- R, y4 o. u) f# r  s6 M        double dT_Hq_dCA();
        double dT_dCA(double an);
2 P( d# N2 _+ C' p        void Calculate();

) i2 D6 B" ]0 U: E8 q, z2 ]/* A编写各子函数程序*/
        double Cv() /*求定容比热*/
5 x. ~1 B! y9 o& h- l{

}
& |' S6 o' D9 E' V$ v

double  dV_dCA() /*求dv/dφ*/
{
& R. ^8 Y1 {* b8 z2 |) m
/ ^; L: u: f% V2 O! i2 B/ G}
: y( q. i0 B7 ?; D% v

double V_CA() /*求气缸容积*/
{

0 C2 u3 l3 T" A2 u+ I}
2 v) }# N5 S( I, ?5 O7 m5 E
; ?" Z9 d1 ?# z
1 R4 n4 P8 N- V# Y- K5 C/ ldouble dQw_dCA() /*求dQw/dφ*/
{
& f7 L# \' v$ P- F
}
( G5 N6 M% m3 y, G$ l' R, ^
) T6 o% t; c) w; N( G
double dQB_dCA() /*求dQB/dφ*/
{

}

: ?) g- h$ |. E* p5 x& G
3 |" l, K! K( w- N0 Adouble U() /*求u*/
) A& a9 l: |$ S( k2 z{
4 H0 z% k# e9 X: x
, T' K* }: l9 Z2 ?6 \  Z}
% D0 c- a4 d0 W

double dU_dnmd()  /*求du/dλ*/
8 P: ^3 \: y& a; f{
0 Y/ u/ K  o. N/ F7 f! L( C  s3 M$ [
}

/ i( o0 r+ u4 I& m9 V) }double dT_Ys_dCA() /*压缩期*/
8 U% K: B4 F1 A' ^! _{
7 X# Y2 i2 g) L! N. ]8 E) z5 L
$ k. p2 x  h4 M( N}


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

}
1 A) m: J/ p* e% U) S
8 F  U: X2 r) i. n; [* j
; i1 J% x" P$ O/ ~6 A! _# e, hdouble dT_Rs_dCA() /*燃烧期*/
1 o  y7 o, T" w& ?( d' U{
0 X% q6 w0 }9 Q; E, ~# f# A
}

8 u' K1 B7 L% P
double dT_dCA(double an) /*求dT/dφ*/
" K3 w0 J) i4 M" p1 n- ?{
        CA=an*PI/180;
, A' t8 K! B( ?" z7 y        V=V_CA();
        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
; L4 f  \& ]# \- ~$ K$ }        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
0 S  }7 s& s9 ^: |        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
        if(A_EO<an && an<=A_SO)    /*换气过程不计算，线性赋值，避免零值 */
/ t, v5 i) a$ ^" I; i                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;
) f  O6 G7 J5 G6 |4 L}
. h* M4 @' K# T% N6 ]
void Calculate(void)
/ v: O. q8 k! P& {6 v{
" g, A# }$ S' V5 Y3 s' L) ?9 B    double Tpre,K1,K2,an;
; I( E0 {8 a- Y) D    int index;
6 y5 L2 f% ]* ^& E        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效压力bar*/
        mBo=Ne*ge/(30000000*n*gnum);   /*单缸循环供油量kg */
0 K2 _+ n7 L$ m1 T5 L        CA=A_SC*PI/180; /*进气关闭时的角度转弧度*/
        V=V_CA();/*进气关闭时角度下的气缸容积*/
/ K$ G. k( E1 ?        P=P_SC;
        T=T_SC;
        index=(int)A_SC;
: `# H: R0 `. m( q$ H0 }        m_T[index]=T;
        m_P[index]=P;
3 c0 _  i0 Q, V  Q  u        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
& L) E: o. [# q        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
8 s* [  }3 F6 l5 s9 H    mf=m-mL;
: T/ }8 u* M! `        Tpre=T;
        an=A_SC+1;
' `- K, P2 G" ?- Y8 o3 z( q! s6 }
2 T  S7 d, h7 t        do{
2 l7 n+ J+ X9 r) Z4 `: h; n4 R2 I
" r/ v. ], V& \! ^" Q6 U  /*B 编写欧拉法求解常微分方程程序 */

! ]( o, b9 i7 i- J  S; E7 H* t        }while(an!=A_SC+1);
4 X4 l; U: K* A- ~9 x' Q9 y) R}

3 k* @* N1 d8 f+ S2 \( S# pvoid main(void)
{
        int i;
# k' [, k" I  ?7 S        FILE *fp;
        D=0.11;   /*发动机数据按各自的任务书输入 */
5 |7 G1 L- i; A( @" `        S=0.13;
        L=0.21;
# B' \' r2 m2 g+ u4 X+ c) z+ i        ge=236.0;
        n=2400;
        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
        R=287.08;
& \& _* n0 t4 _8 F        Pa=101000; /*大气压力*/
) s! ^) z/ g+ x7 _9 y; ]0 T0 z        Ta=300;    /*大气温度*/
0 f3 W. m  x: Q        hs=(735.0+R)*Ta;  /*大气焓值*/
        P_SC=90000;
        T_SC=313;

( Q6 t% \7 |) p: t. W        A_SC=216; /*进气门关闭*/
7 }( K% j( u/ G. O        A_CB=350; /*燃烧开始*/
- T. q0 M* O9 X6 ?! m4 h# d        A_CE=410; /*燃烧结束*/
, @$ ?) i; B; J# T& A- B        A_EO=484; /*排气门开启*/
# b, e/ E$ m; S+ A        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/

- I+ C' s& S/ ^/ i8 I        Lo=14.4;
% P0 p0 R' }  B9 j+ j, s
+ N* [- {$ n. R2 H    for(i=0;i<720;i++)
' i2 D" d8 Z6 x1 o  E3 H$ o$ M" k        {
                m_P[i]=0;
                m_CA[i]=i;
                m_T[i]=0;
        }

        Calculate();
" d- O6 s9 i3 B2 L; y. s1 Q8 |8 d4 r
    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */

/ g5 I. P: J7 c}


. O2 f: a1 H; A1 \% N( @: f( i2 |1 H

祥云飞鱼

不难啊  完全是 把 公司 翻译成 C 又没有用指针跳转的。。。

一点一点来

这是作业的干活吧……