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

方进小侯爷

#include "math.h"
# ^" n; ?/ p* F" p/ O#include "stdio.h"
# a6 u% k3 K& t#define PI 3.1415926

double D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
3 {$ x' O# ^3 P, [( y' S/*缸径、行程、连杆长、比油耗、转速、功率、气缸数、压缩比、低热值、理论空燃比、Weibe函数常数m、气体常数、进气门关闭时压力和温度 */
- q# h' y8 `; O! _        double m_T[720],m_P[720],m_a[720];
" m2 B7 _" `# _/*气缸温度、压力、瞬时过量空气系数 */
        int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
4 s. \8 y3 w* A1 r; d" `        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
# d# x! \: ?1 K+ F! Y/ h        double Cv();
        double dV_dCA();
: D; v8 b0 C% D9 z, H* d" @' ]        double dQw_dCA();
        double dQB_dCA();
        double U();
1 h' a' z9 E& h        double V_CA();
        double dU_dnmd();
( g4 e% b8 W# z2 w* L! P6 k        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() /*求定容比热*/
{

}
' O* O  i6 h2 f/ a: C

double  dV_dCA() /*求dv/dφ*/
0 {, N& c' }, |8 _* ?9 l{
" a2 w; T. V  J( m
}


+ v8 G' G1 L# H( R" `4 D1 Gdouble V_CA() /*求气缸容积*/
{

/ @9 @6 v$ c8 W( A; a}
0 B. b. L7 ]4 V8 P
& h3 x& ]4 m! b3 p+ Q  i: r, ^, B5 s
3 ]* Q6 L, |4 Zdouble dQw_dCA() /*求dQw/dφ*/
{
% t/ U: h! d/ }
7 o+ P, F$ r; f2 s  M" y8 B}


! E) ^9 E# @3 vdouble dQB_dCA() /*求dQB/dφ*/
$ d6 `: Q4 i1 ?. b' c0 ]( s- {$ v( b{

7 b/ ~1 a+ U- {! Q# ]}

# k4 N7 G) B' I6 t; X0 t4 f
4 J2 W) ]5 {2 idouble U() /*求u*/
{
* G& [. E- c, K0 K
}

" S2 x: n6 T3 w) c& v8 m  q
% Z9 |: k& c3 l: z$ Z: b5 m4 Kdouble dU_dnmd()  /*求du/dλ*/
{

}
3 C( f+ t0 p0 i0 i5 _" f% }* P7 s
0 P! t8 v6 A  i3 f! u* }& [double dT_Ys_dCA() /*压缩期*/
9 F" A( t/ w5 Q' y5 ~- n0 b{

2 \' q3 T* s" p( D" Y}

/ t- R7 n# G0 L9 c& p0 @
" T$ ^9 o* a$ d5 j. X: mdouble dT_Pz_dCA() /*膨胀期*/
{

; a0 S% H% {# M8 j' [- W}


. ^1 z+ L2 M8 A) ~double dT_Rs_dCA() /*燃烧期*/
; X. J) }& r" b2 a5 s{

2 P, p1 Y4 ^# a( {6 u}

& o* }  {+ o; @, M/ j/ V1 D4 X% |( C
( j% z; ?9 _* d1 R& Ddouble dT_dCA(double an) /*求dT/dφ*/
{
        CA=an*PI/180;
        V=V_CA();
2 t% G" g! G. [2 M9 J7 X. n        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();
" H2 D4 Q. x# T5 \        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);
; s/ E0 \( g; \; s7 }0 O0 z5 s        if(A_EC<=an && an<=A_SC) P=90000;
        return 0;
}
% c( R5 w. h) c9 B
void Calculate(void)
+ U1 ?2 k0 ~8 ^% s* X/ h{
    double Tpre,K1,K2,an;
- m# q% w0 |! t/ X; p' H5 }    int index;
7 F4 y8 b- S) e3 e        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; /*进气关闭时的角度转弧度*/
+ z5 H8 S5 M9 S9 g6 n  N+ P7 \        V=V_CA();/*进气关闭时角度下的气缸容积*/
        P=P_SC;
        T=T_SC;
4 b7 I2 h  K+ Z$ b2 y        index=(int)A_SC;
+ V( \" U. ]+ I9 \        m_T[index]=T;
        m_P[index]=P;
9 @+ [3 d7 b" Q  _; r# k        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
    mf=m-mL;
        Tpre=T;
        an=A_SC+1;
( N$ G6 Q% O  ^1 K' f$ [% o' \: S
7 I; p1 g3 B4 X$ B2 R        do{
$ D4 Q3 y$ v9 h. g2 T- B: [
1 ?* c' i1 y/ z2 N  n7 B  /*B 编写欧拉法求解常微分方程程序 */

: C& P8 `  A- @1 G, d, W+ ^" p        }while(an!=A_SC+1);
* J& z. Y  ~" B$ m6 e2 h# c}
7 \! o8 d3 Y. X
void main(void)
{
* D) G* v/ f" _8 _$ y# n! Z        int i;
        FILE *fp;
        D=0.11;   /*发动机数据按各自的任务书输入 */
        S=0.13;
% G/ H- B, h5 d5 @% w* [' O$ ^        L=0.21;
        ge=236.0;
+ c5 x4 A2 h& H. a  a        n=2400;
% s( g, ?6 s! W) t* p8 a9 Y' C% ^        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
) `* b8 @8 d9 b- m; F* _, v  M        Hu=41868000;
        R=287.08;
        Pa=101000; /*大气压力*/
        Ta=300;    /*大气温度*/
9 ^" w6 o9 ~& M6 b- H6 k5 R        hs=(735.0+R)*Ta;  /*大气焓值*/
        P_SC=90000;
4 e( c/ u; X) _' S2 {        T_SC=313;

' v. c9 ?+ u, `4 e. ^( N$ t4 }        A_SC=216; /*进气门关闭*/
        A_CB=350; /*燃烧开始*/
  J, R. C6 M9 s" |% p- ^, j        A_CE=410; /*燃烧结束*/
8 o7 J$ T& Q4 {; ]- k) ?7 q, c        A_EO=484; /*排气门开启*/
; s1 W0 k  z' N' J; B; I        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/
3 N8 M' @; K. p. }% F' u2 ~0 e
* D" i1 [' i4 y: B# U        Lo=14.4;
2 p9 J$ ~  U6 f2 j# O$ q+ n+ r% H
    for(i=0;i<720;i++)
7 K" R/ b. }) C* F        {
8 h) F% \$ n+ c; n/ N' m                m_P[i]=0;
, @! K6 M6 K! L* r% M' W                m_CA[i]=i;
                m_T[i]=0;
        }

% J5 N) C. M9 J6 W        Calculate();

; T  x- i/ A+ M9 L2 W6 m    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */

}

6 ~  V: @$ ~& q$ y
* W4 S9 l) ]; V" o+ Z' q

祥云飞鱼

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

一点一点来

这是作业的干活吧……