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

方进小侯爷

#include "math.h"
#include "stdio.h"
2 ~' _! m- L( w, }#define PI 3.1415926

& d' s" V- Q! Q& }double D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
( V. Z$ |4 f5 \) Q; q/*缸径、行程、连杆长、比油耗、转速、功率、气缸数、压缩比、低热值、理论空燃比、Weibe函数常数m、气体常数、进气门关闭时压力和温度 */
: ~+ H4 g' b! X" I* E        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;
        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
        double Cv();
. S& b8 C! n- y7 W7 ?% w        double dV_dCA();
        double dQw_dCA();
        double dQB_dCA();
        double U();
0 p& u2 G0 F6 h1 y% R6 E        double V_CA();
) L0 N2 _( N5 R1 Y2 ?        double dU_dnmd();
1 a* L# u' z! \8 s; w7 f$ Z        double dT_Ys_dCA();
        double dT_Pz_dCA();
        double dT_Rs_dCA();
4 J! ]. k& |! r' \2 S        double dT_Hq_dCA();
        double dT_dCA(double an);
        void Calculate();
' Q* d# j( G+ J0 h
$ y% A# e& h; q/ |6 H8 L/* A编写各子函数程序*/
3 [# P. j* `# g5 x* Q        double Cv() /*求定容比热*/
5 B1 m7 l$ P$ b5 L1 p{

. K6 j* I$ k! F}
/ a) L: g  {# S: N3 @

double  dV_dCA() /*求dv/dφ*/
* G: E, _' U# t/ `! {! ~# c0 X{

& j1 W3 l6 ~: n7 g  c0 _}


) H4 U9 @# m+ U: b. Cdouble V_CA() /*求气缸容积*/
{
7 y( F& b6 _) O' h( P) b
}

7 ^: i5 D0 [; L9 X4 R' J1 ^
double dQw_dCA() /*求dQw/dφ*/
$ g$ W! _2 I% @9 x; _/ S- N{

* e; B( L. E9 j( c}


$ O& c+ j1 b# |, B6 idouble dQB_dCA() /*求dQB/dφ*/
{
7 U+ Y9 i4 m  u5 E' v% g7 T6 P# y( g
}


- H- C) L* g+ D7 ^6 ~% @% Kdouble U() /*求u*/
{

$ l( `: B; \/ f8 u" t/ }8 T}
" `9 l! D' p1 c

double dU_dnmd()  /*求du/dλ*/
{
7 n$ Z% g; O- y8 i; F; U
+ ^) x2 l( ~9 |, E! g0 t, S}
; L2 `# ~  j: O9 ?8 }( [+ G% G
, [3 H  ~6 P+ Tdouble dT_Ys_dCA() /*压缩期*/
{
  G# b) C' U$ J3 L/ s
# a# X# x) y* v  v  g}


/ V5 m- A" q: l. @, Y/ P8 pdouble dT_Pz_dCA() /*膨胀期*/
( v$ o7 P; k$ T7 `) Y9 _6 V$ C; y5 @{

}
9 k0 ~  m' o2 r! w" e+ A/ z

9 O: D$ O2 C2 @& ]- @1 udouble dT_Rs_dCA() /*燃烧期*/
* Q, `9 M0 X/ ^) n% L4 C{
8 ^5 d7 G( a: @
) R& S8 `) E3 [& Q}


5 w& l* l1 h: d7 @double dT_dCA(double an) /*求dT/dφ*/
+ C' w5 F3 A' G  S& k{
# g5 d: R1 d& t7 p6 U- y        CA=an*PI/180;
        V=V_CA();
3 G5 h( \. m) i" t  }        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
; q$ |* r) A- \8 q  {        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);
3 _! J/ f: \" O8 S9 u2 v    if(A_SO<an || an<A_EC)
- R1 \9 N- S! g" k+ ~" \% f( }                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;
5 [* Y1 n2 F3 O9 x1 _6 F% g        return 0;
% j: m& `5 V( F! P) q; S}

; P( @; T/ r! E( }void Calculate(void)
& K) C4 U6 j8 s* \7 l{
    double Tpre,K1,K2,an;
( c# V$ o& Y, u2 E) }1 E; n/ [    int index;
        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; /*进气关闭时的角度转弧度*/
) P$ u; ~# s' h4 F9 R2 k        V=V_CA();/*进气关闭时角度下的气缸容积*/
4 c. b7 F0 [9 ~/ z" w: j( H        P=P_SC;
" @. \. a9 n( k        T=T_SC;
        index=(int)A_SC;
' ]; Q( R, V; Z* ]        m_T[index]=T;
        m_P[index]=P;
        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
6 I, v! g/ d9 [& o" Z3 l! K        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
) q/ }. w2 _( q3 l3 w    mf=m-mL;
        Tpre=T;
        an=A_SC+1;
! I. ?0 I# P+ Z4 l3 m4 c6 U, r! a
  Z1 I- w+ P7 M  f( [$ r        do{
% K/ n: N0 f3 W; B0 s! W! h
  V' e2 o" l2 B' L  /*B 编写欧拉法求解常微分方程程序 */

        }while(an!=A_SC+1);
) M6 ^( a8 a8 [8 {7 z# h# |4 }}

void main(void)
{
$ _, L4 b0 g; A# y3 \# j4 J8 r        int i;
        FILE *fp;
: s" J) k& X6 ?' Y' X1 @4 }( G        D=0.11;   /*发动机数据按各自的任务书输入 */
        S=0.13;
6 m5 ^! I4 f$ `, C. Y        L=0.21;
        ge=236.0;
        n=2400;
        Ne=101568;
' Z( }3 i" k/ _0 ^+ j! [        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
% r# q2 Z. ]$ z( {6 _5 J        R=287.08;
4 A% y9 g6 j! P2 X4 o7 I1 h6 q* I        Pa=101000; /*大气压力*/
6 ~" v" a& U( O* W4 n        Ta=300;    /*大气温度*/
        hs=(735.0+R)*Ta;  /*大气焓值*/
        P_SC=90000;
        T_SC=313;

# F$ `( b3 w+ k0 g7 T+ e        A_SC=216; /*进气门关闭*/
) G1 G2 t: H9 v5 a+ x        A_CB=350; /*燃烧开始*/
        A_CE=410; /*燃烧结束*/
! p0 Y8 ]" h$ e) g* \8 m% S1 S5 {, J        A_EO=484; /*排气门开启*/
; l' f" R1 J5 Q& F: }( k        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/

        Lo=14.4;

    for(i=0;i<720;i++)
/ b. x" P' j0 b$ s5 }        {
                m_P[i]=0;
                m_CA[i]=i;
; s0 u) Z6 ?4 p6 b                m_T[i]=0;
        }
  @: E( p5 y  P* `$ \2 g4 S( Q5 f2 t
        Calculate();

    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */

& O  l5 c) `+ \}

0 W8 `/ Y+ v$ g& h' T+ I
& x& ]' J9 t) S

祥云飞鱼

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

一点一点来

这是作业的干活吧……