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

方进小侯爷

#include "math.h"
#include "stdio.h"
1 N6 o& M8 w9 M8 q9 R#define PI 3.1415926

& q' o! X% U' K5 Adouble 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];
' ]- W; C+ E( _/*气缸温度、压力、瞬时过量空气系数 */
+ n) x9 _/ M+ w3 A! H7 Z        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();
8 |- k% \2 K6 E8 x; n% V        double dV_dCA();
# S( Q, x+ J6 p8 @& ]7 D2 F        double dQw_dCA();
+ Z" Q. @) C" d: g        double dQB_dCA();
4 h7 s7 Y/ z5 @        double U();
8 A0 C+ ~1 S- A        double V_CA();
+ c! I- n  }! }8 \" j5 ~9 c; d        double dU_dnmd();
        double dT_Ys_dCA();
        double dT_Pz_dCA();
" r- y" J2 L. y2 \) S# l) X8 ?: z        double dT_Rs_dCA();
) o; H4 Z  m0 d; g        double dT_Hq_dCA();
        double dT_dCA(double an);
) U4 s  E6 h6 L1 T. b        void Calculate();

/* A编写各子函数程序*/
        double Cv() /*求定容比热*/
{
- V5 K& f# [8 |
  F! l2 Y3 M6 V}
( w* A7 @& [$ `) s
' N* S3 ?/ g, f3 f- i$ _4 S0 Z
) i. Y& m+ U# l8 J9 h, U9 u3 Qdouble  dV_dCA() /*求dv/dφ*/
{

}


. h+ C: b9 r* m; jdouble V_CA() /*求气缸容积*/
- q% o7 u9 ?/ ], _' P. e* g{
) u$ J' P- y7 g
3 A( l8 e/ H0 c}


double dQw_dCA() /*求dQw/dφ*/
{
0 J+ S6 u# h7 R6 v8 f
9 z7 J' D7 l1 W1 t& w/ S" M1 r}
: N5 L$ C8 T3 P; D/ H
. Y3 P  u- k& Y8 Z" O, O
6 e7 F# ~( N: l: i/ Vdouble dQB_dCA() /*求dQB/dφ*/
{

}

( P6 y  ?9 L' ?% x# I) P) R$ ^, L
- |$ ~+ W$ {3 J1 I0 D; ydouble U() /*求u*/
" I% B7 z. q% f) F{
# R' B2 B- @  K$ p
}


5 v( p% b$ i' _( `: d$ odouble dU_dnmd()  /*求du/dλ*/
{
9 ~  r, S+ Z& m8 |7 z4 m/ z
}
7 e6 b1 I+ ?* W9 A, D% l2 x
4 |0 v8 P, ?0 G+ s* ^- vdouble dT_Ys_dCA() /*压缩期*/
4 l1 E5 M/ _& v2 K2 A+ Q; y{

}


  Q5 i, r% _# l# i2 jdouble dT_Pz_dCA() /*膨胀期*/
4 p6 ]) M, e5 Q8 L{
, D8 P" \2 W/ r2 T3 b4 H
( G$ N3 {0 O1 o& A6 Y* n6 y  G}


double dT_Rs_dCA() /*燃烧期*/
5 W) ~3 a- y2 |3 Z- ^/ o$ R9 k1 H3 j{
) H, S* o0 }4 ]$ H* d
9 E" e/ x  i1 @; Y9 D1 s( V}
7 S$ P. O2 r7 L! c

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();
' B3 T! @5 C& @  B4 @        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);
' F* z6 F. w) N    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);
) ~, Z# |. a# d7 v        if(A_EC<=an && an<=A_SC) P=90000;
: M+ k7 U3 }4 }. p3 y$ |- n) j( {5 F' X        return 0;
" z1 n, ~/ v; H: Y. ~7 y}

( w, ?5 L8 D) r' _7 s  e3 Z5 o& n6 _+ avoid Calculate(void)
{
$ n* t' i6 x7 I0 b4 i3 B    double Tpre,K1,K2,an;
    int index;
+ q; k& e5 O9 }8 [        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; /*进气关闭时的角度转弧度*/
' A; i8 L( o: y9 k  ]9 ^1 p! a0 x% w        V=V_CA();/*进气关闭时角度下的气缸容积*/
        P=P_SC;
2 S# x, H' P- A& a        T=T_SC;
        index=(int)A_SC;
3 D9 G- Z; s# a) N5 `+ J+ p, E        m_T[index]=T;
        m_P[index]=P;
        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
' P; C9 H& T: d1 m. g/ c6 k        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
    mf=m-mL;
        Tpre=T;
        an=A_SC+1;
* }2 l1 G, B' i/ c( {% f
        do{
8 V  a& z: @  N1 `
  /*B 编写欧拉法求解常微分方程程序 */
4 @% n. F. A9 A5 t
1 ^0 W3 a" N& ]! C* T1 r        }while(an!=A_SC+1);
8 n7 i3 P" m$ ]  T2 b7 H}

! {  B" u' J" E' R$ o& a8 svoid main(void)
0 B9 t+ w, W3 R% H. c. e. I  K& I. f{
" m7 R4 d0 }! v' c7 ^' Y+ D" U+ q        int i;
; M8 y' g7 h! ~# i9 }8 J1 A        FILE *fp;
        D=0.11;   /*发动机数据按各自的任务书输入 */
) t$ J# a. J6 z, q2 U0 G! ]' z, }        S=0.13;
3 @# `5 S1 O( V3 P, D4 P        L=0.21;
        ge=236.0;
        n=2400;
        Ne=101568;
' R# x& u) c4 ~4 v+ {        gnum=6;  ysb=17;   M=0.7;
/ P7 {# m6 ]) e4 p$ o        Hu=41868000;
$ D2 }8 S) f1 [& |( B+ t. |6 s        R=287.08;
4 ]9 g  V' U0 k+ F+ o) c% D) Q( k& B        Pa=101000; /*大气压力*/
7 f; Z; }2 F% `$ U; N        Ta=300;    /*大气温度*/
        hs=(735.0+R)*Ta;  /*大气焓值*/
+ X( k; t0 j0 p4 G6 Q        P_SC=90000;
        T_SC=313;
' [1 j; {8 N; ]1 Y, ^3 P
        A_SC=216; /*进气门关闭*/
        A_CB=350; /*燃烧开始*/
' c3 N- E" T* o9 J* ]        A_CE=410; /*燃烧结束*/
! a9 @  d  E1 K" E) ~        A_EO=484; /*排气门开启*/
        A_SO=705; /*进气门开启*/
4 O$ O4 Q, G0 w8 y  U1 F        A_EC=15;  /*排气门关闭*/
6 }' G+ x! V; J
        Lo=14.4;
1 ~% g9 M1 Y8 _) L& }
    for(i=0;i<720;i++)
        {
                m_P[i]=0;
                m_CA[i]=i;
2 Z; T& u$ w; Q6 M9 ~# T                m_T[i]=0;
        }
/ Q3 R) Z; G  {# C% v# S9 B8 y
( o: E% B1 }$ _1 q( f6 `% Q7 \        Calculate();
" c2 |- c# D5 G% |
    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */
( o# B: s0 ]/ ^
}

' z4 M7 p- n$ O2 y8 E9 ?; x

祥云飞鱼

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

一点一点来

这是作业的干活吧……