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

方进小侯爷

#include "math.h"
#include "stdio.h"
#define PI 3.1415926
) x. j, n' D( I
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;
* |; r) c* L+ @& t+ p. Z        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
$ z" W' L  ~2 Y, G        double Cv();
, C$ X, q1 ^" h, K7 C: N- R        double dV_dCA();
        double dQw_dCA();
# e) k: I" O" H4 ~6 u) G$ j        double dQB_dCA();
        double U();
        double V_CA();
        double dU_dnmd();
        double dT_Ys_dCA();
$ {. y; j- A" z- n3 Y0 `        double dT_Pz_dCA();
2 S5 X& b4 w. Y& _, S1 I        double dT_Rs_dCA();
0 C5 W, z3 p1 P2 s1 n9 \% Z        double dT_Hq_dCA();
1 {( D: y" s, f+ a. e7 ^        double dT_dCA(double an);
6 R. M& D& i  A7 {3 d  M        void Calculate();

% S( n2 \0 T- M$ j0 D" V2 m/* A编写各子函数程序*/
* E) u3 A) j: P        double Cv() /*求定容比热*/
{
/ _  T  V- T# m8 p6 l
6 n/ {5 o1 _' W# Z. l}


, K$ B% ]6 z1 N6 Qdouble  dV_dCA() /*求dv/dφ*/
4 n( l. I$ K" h{
6 O5 n' P) b% D$ }, @; A- M
! _+ F2 p1 n& r3 k- ]! z}
$ N# ^: v3 U1 s, V/ Y
- N- i2 B! |0 P- L; J4 `9 d
. a' d% B: x: [+ j& {. M! ddouble V_CA() /*求气缸容积*/
{

}


/ I3 _2 R7 W1 c( ~. W! K% }' \double dQw_dCA() /*求dQw/dφ*/
' z! g/ Y# m4 u7 Z) u' @' X6 N: U{
9 ~$ I/ T) c3 ^& W  ]
}

7 I. r& f- T- {' g. L
double dQB_dCA() /*求dQB/dφ*/
{

  m2 O* i- a' T  P3 C}


double U() /*求u*/
2 R+ k1 t3 e% I6 _" V. c# q  X1 {+ k{
+ L& j' G6 }) I, g' X5 v
}
2 @! b' Y" l6 x. J

  Z! X9 _9 g# x( U# A# mdouble dU_dnmd()  /*求du/dλ*/
  ~! Q& Z3 R! Q5 e: T{

3 B" F7 J+ J5 c- H}
. g( Q: Z; \; ]5 `9 Z
double dT_Ys_dCA() /*压缩期*/
: Q5 g% Y$ a4 V% S8 \2 H{
" Z8 Y& l8 b  y+ m6 m7 X; g/ i2 G  l
}


# u* l. R% X% Ldouble dT_Pz_dCA() /*膨胀期*/
{

* W( c* G9 N3 A: X}
8 d6 x. c/ i; E

4 [( b/ C# }; d: I2 vdouble dT_Rs_dCA() /*燃烧期*/
{

}


" C  {8 Z( U8 z8 [& ndouble dT_dCA(double an) /*求dT/dφ*/
  f: O" J9 J' P% u) ~/ `' Y{
        CA=an*PI/180;
6 m' d- O+ S! w6 V( q        V=V_CA();
$ ?* c+ w1 M, y! x        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
7 j7 B: D) P* r8 n* S0 p/ f        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
# `1 {6 \5 c6 ~        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
# u( s- y) @# U1 N+ b! s        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;
* r! r% n7 J6 `        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效压力bar*/
7 P) ]; e/ L* _+ c: T        mBo=Ne*ge/(30000000*n*gnum);   /*单缸循环供油量kg */
        CA=A_SC*PI/180; /*进气关闭时的角度转弧度*/
        V=V_CA();/*进气关闭时角度下的气缸容积*/
/ }) [# X: D8 m% O6 L( H        P=P_SC;
        T=T_SC;
; z/ v+ B0 Z2 N        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{

! S; v" o5 ]0 O) T7 ]  /*B 编写欧拉法求解常微分方程程序 */

5 K$ y8 A; u. W! Y( ^- x2 j3 o& G        }while(an!=A_SC+1);
}

- l0 ?  h- `# V0 o# gvoid main(void)
9 \" H+ h1 K6 O{
        int i;
        FILE *fp;
; s6 D& _; I+ }5 G' }4 O. d  v. v        D=0.11;   /*发动机数据按各自的任务书输入 */
        S=0.13;
        L=0.21;
( }+ l9 m& U" Q. i# j3 R  f- C        ge=236.0;
        n=2400;
        Ne=101568;
" _4 }0 \) C9 I6 e2 l  H* A        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
, o2 ~0 e/ v# s, V5 h0 w        R=287.08;
7 p* q1 r. p; }# j3 K% s2 T        Pa=101000; /*大气压力*/
        Ta=300;    /*大气温度*/
7 |0 Y& w4 E% n  n3 e        hs=(735.0+R)*Ta;  /*大气焓值*/
        P_SC=90000;
5 R/ }& f0 {2 T3 k+ m+ K( W1 c        T_SC=313;

        A_SC=216; /*进气门关闭*/
- l- |( m- A6 @6 J        A_CB=350; /*燃烧开始*/
: X6 I  `+ }+ _0 I+ f& r% }+ `        A_CE=410; /*燃烧结束*/
' a( M) Q' g$ L, |$ L3 m' Z  u2 f        A_EO=484; /*排气门开启*/
! w! E4 `. o# s9 C5 u! K2 c' W        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/
, W1 l$ Q( N' h" E' v7 D% T
. P4 R% Q3 M8 ^' c        Lo=14.4;
7 S& E0 G4 l0 }: G
$ }7 d; P: D) f; Y; O7 @    for(i=0;i<720;i++)
9 ]: U5 P2 z# `        {
, v7 A0 A2 _; m+ V! K                m_P[i]=0;
                m_CA[i]=i;
" m& J) ^' x8 E4 Y; M0 G3 i                m_T[i]=0;
        }
: T8 G7 D" h! Q+ b3 X
: c) u2 a, E3 z2 [, z        Calculate();
6 x7 K( D7 y, N/ B1 q& X1 X
    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */
: ^, @& B7 a' b0 U5 k
}

. ~# O4 Q: B, J% F
  I1 \6 l* `- E$ g( A6 i- f: x& N

祥云飞鱼

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

一点一点来

这是作业的干活吧……