內燃機熱力計算怎么用C語言編程？模板如下，求助。

方進小侯爺

#include "math.h"
+ ^" A4 w) U6 g, l7 |% o#include "stdio.h"
#define PI 3.1415926

double D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
# x! |, f* j; i9 C9 v, t/*缸徑、行程、連桿長、比油耗、轉速、功率、氣缸數、壓縮比、低熱值、理論空燃比、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;
        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
$ z& j- a' p. x9 H6 a8 Z' @        double Cv();
        double dV_dCA();
/ Z' G9 ~8 G  r" s8 E% N        double dQw_dCA();
        double dQB_dCA();
        double U();
        double V_CA();
) o) v0 R3 f5 K4 V' q        double dU_dnmd();
        double dT_Ys_dCA();
8 k6 A: f3 V; @9 U, k! `% F: M        double dT_Pz_dCA();
        double dT_Rs_dCA();
6 C. }) m( K1 P        double dT_Hq_dCA();
        double dT_dCA(double an);
        void Calculate();
7 y1 B+ @# C5 p. T6 P
' J- f7 W8 `! J: M2 \/* A編寫各子函數程序*/
        double Cv() /*求定容比熱*/
  S; I, u- P# D6 Y0 p& _{

2 }2 S& B) V- o/ o/ A2 _}
: ?( O& A6 I3 P! {; u% ~% n
' e1 _. }" S* i1 j0 g8 y8 A" |, r. H
( J% i: x* z( r9 c* M9 edouble  dV_dCA() /*求dv/dφ*/
{

4 U2 h; v6 q$ y" g% [% W- F}

% U8 d- `( z6 e4 ]4 K# H
* E/ w+ W# |2 Y( C7 Q( gdouble V_CA() /*求氣缸容積*/
. x$ W3 N  r+ _3 @+ R5 f{
" l5 W: M& [! I4 n+ S+ P
' c6 Z. u. f: _' K4 l% X! U7 w}


double dQw_dCA() /*求dQw/dφ*/
$ |& O- E3 s0 M9 e! L$ n! v{
; v  W, x2 i+ [+ L7 ?
4 p0 @' N' u! t}


0 k7 Z1 }" q# a" J0 A- B& \double dQB_dCA() /*求dQB/dφ*/
{
: V0 I% e$ l" L9 P  _1 f* k2 `
}


double U() /*求u*/
{

}
$ O2 W2 \5 H; ]# i% F4 i

- F% j8 ]% e: i/ C: |double dU_dnmd()  /*求du/dλ*/
{

* {/ F' J) z) k" O1 E  L}

double dT_Ys_dCA() /*壓縮期*/
{
' d9 Y' x% [; a1 B8 e4 a" x, h
}


double dT_Pz_dCA() /*膨脹期*/
{

. z+ V+ {* }! z+ C}


) M) k& U* l9 x, z, ^- [double dT_Rs_dCA() /*燃燒期*/
2 K. g4 _' B- h7 d: n; K{

}
$ `5 K) R+ K  s7 t

double dT_dCA(double an) /*求dT/dφ*/
6 g- C9 I, m" p5 _( z, I, M5 B( f{
) g& R+ ^0 F  J        CA=an*PI/180;
        V=V_CA();
        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
7 \5 `1 ~5 {. C% x3 a        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
* e( r7 F  r( X" j5 u8 C        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
9 ]7 r! b6 B7 h# L        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)
+ G7 ~0 C  E6 Y5 Z6 P  Q7 @                P=110000-20000*(an>A_SO?an-A_SO:an-A_SO+720)/(A_EC-A_SO+720);
/ o0 O* Z" M! ^. i! c  a* W        if(A_EC<=an && an<=A_SC) P=90000;
        return 0;
}
: w- Q) S9 [, ]* z
void Calculate(void)
{
    double Tpre,K1,K2,an;
1 n9 x/ g) x5 x- k    int index;
* b/ T- K. N! k9 ^        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效壓力bar*/
        mBo=Ne*ge/(30000000*n*gnum);   /*單缸循環供油量kg */
  o' K/ t* I) p$ [! \2 a8 x7 H        CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
& X6 `& t6 S3 _* y        V=V_CA();/*進氣關閉時角度下的氣缸容積*/
# n/ n* e# q( B$ o; o% M6 E        P=P_SC;
        T=T_SC;
        index=(int)A_SC;
+ ^# L2 n" b# r  S) ?8 ]        m_T[index]=T;
) S8 J6 H8 u  T6 t$ e; ?        m_P[index]=P;
: |4 z- Q" j0 f% I$ w3 p        m=P*V/(R*T);/*進氣關閉時角度下的氣缸內工質質量*/
        mL=m-0.01*mBo;/*進氣關閉時角度下的氣缸內空氣質量*/
. _  L, ?' r$ z6 w: R- y    mf=m-mL;
+ ^. X) F, ^% Y2 a- {1 \3 a+ l        Tpre=T;
! Z8 U( L9 F& T4 S        an=A_SC+1;

        do{
& k  R  y3 \: z+ p0 B4 R
  /*B 編寫歐拉法求解常微分方程程序 */

1 k' W8 [8 R/ ^* \        }while(an!=A_SC+1);
4 w1 k( k' x. L5 o* Q9 A}
# S, E0 O4 G3 l
void main(void)
1 c; r& S2 e2 o/ P{
" D+ Z% ~2 D4 p        int i;
        FILE *fp;
        D=0.11;   /*發動機數據按各自的任務書輸入 */
5 [; J  N& `. a/ n: s0 @2 b. @        S=0.13;
        L=0.21;
        ge=236.0;
        n=2400;
        Ne=101568;
) k" a1 e& L, C0 U% s        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
        R=287.08;
        Pa=101000; /*大氣壓力*/
: w) t* v: J+ m1 \        Ta=300;    /*大氣溫度*/
        hs=(735.0+R)*Ta;  /*大氣焓值*/
! N+ |, L5 H4 ^        P_SC=90000;
        T_SC=313;
4 T( q0 o+ \% m2 ~+ L, [! @
2 O+ `  }9 s, V2 M# _" c        A_SC=216; /*進氣門關閉*/
        A_CB=350; /*燃燒開始*/
* C; a; v( K) B) ^' h+ P7 `2 w        A_CE=410; /*燃燒結束*/
0 T& a( X& o# D+ S        A_EO=484; /*排氣門開啟*/
        A_SO=705; /*進氣門開啟*/
        A_EC=15;  /*排氣門關閉*/
2 v! _& g( _0 c0 y' h
/ P5 P, P2 F- f, p        Lo=14.4;
% Y% r1 |; g2 C- d
    for(i=0;i<720;i++)
        {
                m_P[i]=0;
                m_CA[i]=i;
  f3 |& j9 g" J: j6 Y* [! ]                m_T[i]=0;
8 o$ s1 A0 F& X* p" O5 `5 R6 Z- b        }
8 a  v$ m1 I1 j9 }/ J
7 j( b9 C2 J: i        Calculate();
) \+ C6 Y" I2 @1 f' K
: l5 W1 D6 s, k5 L% w/ b    /*C編寫數據輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數據到outdata.txt文件 */

9 r+ I/ `4 b* c1 V& N( L}

- L* b3 f7 ?8 ~8 h

祥云飛魚

不難啊  完全是 把 公司 翻譯成 C 又沒有用指針跳轉的。。。

一點一點來

這是作業的干活吧……