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

方进小侯爷

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

( l& S% k4 A0 ]1 A5 }+ ydouble 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;
) V4 m6 Z, v. L1 p        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
' ?" X. t( g/ y  S        double Cv();
        double dV_dCA();
        double dQw_dCA();
        double dQB_dCA();
        double U();
        double V_CA();
  o$ q8 d& T9 M) ?& [  m7 {        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();

& E' d- F5 I+ o6 G/* A编写各子函数程序*/
        double Cv() /*求定容比热*/
, M3 p6 k4 g4 d9 q* p, y{

8 h- C( |1 _4 M5 c! ^}


double  dV_dCA() /*求dv/dφ*/
{
; y5 g" D) ~  r0 g
}

4 l8 h+ _* X) d. L' M: k$ Y6 N
double V_CA() /*求气缸容积*/
{
$ Z. W4 v  |- P- p9 X, |
}
1 M& i7 |* y/ P0 P8 z: \
4 _  q3 g4 C3 y( }* \
5 z" ^2 J& z  u* U1 d4 b9 |6 Tdouble dQw_dCA() /*求dQw/dφ*/
+ P8 _& j7 }+ n4 G, G8 @! ]{

}
7 |& S2 s, C/ |. F9 N  ]2 h$ S
$ |* _. i$ _$ u5 A" A
$ A# C/ y, Z, U0 [; p* S/ Pdouble dQB_dCA() /*求dQB/dφ*/
" d7 @0 p9 U6 }% y' R+ }& j{

  g8 k  C- h9 Z2 i}


double U() /*求u*/
{

3 d3 a: |4 ~, h3 K. x8 N/ L}


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

}

0 p; |+ F# T) I/ Edouble dT_Ys_dCA() /*压缩期*/
{
, N8 u: ~. O, E" b, `. @2 [3 k# U
" @  j  |+ n" I4 x: j. K& @1 p}
; a5 o! z& E+ r/ K

double dT_Pz_dCA() /*膨胀期*/
{
6 k& W2 A6 T" X/ a) A
+ ]' Y) g6 X5 K3 b9 p. n# a. O4 r# _}

/ D! d( J- ^* G5 T3 t& n
# S: B! i9 D4 M: u; }2 \) Zdouble dT_Rs_dCA() /*燃烧期*/
{

& q( e0 B8 f/ V* G* Y7 [0 u! G4 f}


+ p7 K% t1 Q4 ^8 ]) m" Adouble dT_dCA(double an) /*求dT/dφ*/
8 x6 u# a1 K; [6 X{
        CA=an*PI/180;
        V=V_CA();
$ n& u' M/ u) O3 a$ ~* B, u        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);
) b0 f" g- h8 ?1 D3 d    if(A_SO<an || an<A_EC)
5 A% j( D, j* u, ~                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;
9 a( X, Y$ \( s  f: u' u        return 0;
}

! z% c/ n! j2 ~% n( Fvoid Calculate(void)
{
    double Tpre,K1,K2,an;
    int index;
7 I+ l  D/ \0 x1 o. ~7 l        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效压力bar*/
% p9 \( c; i/ \: D        mBo=Ne*ge/(30000000*n*gnum);   /*单缸循环供油量kg */
        CA=A_SC*PI/180; /*进气关闭时的角度转弧度*/
: j! |$ }; J' S4 r        V=V_CA();/*进气关闭时角度下的气缸容积*/
        P=P_SC;
1 \& k; m9 U+ ]        T=T_SC;
        index=(int)A_SC;
        m_T[index]=T;
) M; r% Z: M: R' N        m_P[index]=P;
        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
0 N( @7 W! _5 S% `, g( Q% C0 f    mf=m-mL;
( U! t& U) Q: D% e7 v. ^# _        Tpre=T;
        an=A_SC+1;
! L3 U6 M8 Y% F6 Z& @% T' U6 f
        do{
# m2 p0 V$ [0 x6 a( @& Z& b9 _% K8 f
, p$ h1 a" A# X: b! b! L( [  H  /*B 编写欧拉法求解常微分方程程序 */
1 @! G  F4 i2 Q1 S" M
        }while(an!=A_SC+1);
" M/ f: J7 Q$ m/ X6 z9 f}

# I" I$ U* w: N/ p) Q/ k) uvoid main(void)
% k- O0 }) D" P; e% o{
        int i;
. w% w' l' X# c  t$ T' R' |        FILE *fp;
, q2 O! `+ g' p( ~- c        D=0.11;   /*发动机数据按各自的任务书输入 */
9 P0 E- h) y& X" m7 @        S=0.13;
7 ]4 y, u3 ~- z$ `1 \) o/ A        L=0.21;
9 |; ^% _% h0 I        ge=236.0;
% d( c; [  x8 O- n; R0 `        n=2400;
: q* g' L/ F  \% B( k        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
        R=287.08;
        Pa=101000; /*大气压力*/
        Ta=300;    /*大气温度*/
1 h4 t( n6 M* ~        hs=(735.0+R)*Ta;  /*大气焓值*/
        P_SC=90000;
        T_SC=313;

        A_SC=216; /*进气门关闭*/
0 d3 ?' i" g; |1 b# {: B# |! Y4 Y        A_CB=350; /*燃烧开始*/
        A_CE=410; /*燃烧结束*/
        A_EO=484; /*排气门开启*/
+ L# W8 X+ B& C8 f" {        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/

, k1 k. Q1 m6 M1 b0 N, I! F, z/ S3 d9 O        Lo=14.4;
- c. {+ f' ^- _
, k4 t! S  L) R7 }9 t7 {) j- T( f    for(i=0;i<720;i++)
        {
                m_P[i]=0;
                m_CA[i]=i;
                m_T[i]=0;
( W9 t0 ]' Y4 _- T/ x& C        }
( G$ F8 D  t# }% N
. }# b% ?# B5 C) G6 J        Calculate();
; I$ ]* x  ?1 F1 c4 G) Q
    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */
' R/ ?0 V% j* B+ X' E% I% C# [
0 ?; _" m5 m" U0 |5 b$ Z}
9 S* X1 I" Z6 a3 V, n% M+ x
) \/ g/ v! g3 w$ b4 u; K
6 i& C5 v2 N- M8 R8 B$ a8 `9 C

祥云飞鱼

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

一点一点来

这是作业的干活吧……