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

方進小侯爺

#include "math.h"
( Z0 m1 S6 G4 {# h& }) {#include "stdio.h"
! O! ], h4 \' _- Q+ k2 w, J#define PI 3.1415926
4 g/ J6 C$ ~7 m$ h
- x# T/ X* j+ ^double D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
/*缸徑、行程、連桿長、比油耗、轉速、功率、氣缸數、壓縮比、低熱值、理論空燃比、Weibe函數常數m、氣體常數、進氣門關閉時壓力和溫度 */
; c  d8 y7 h& K/ b+ O4 j0 ~+ V0 l        double m_T[720],m_P[720],m_a[720];
/*氣缸溫度、壓力、瞬時過量空氣系數 */
9 b' U2 C' L1 T* F: U8 f. g        int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
. k* s+ C1 F) `8 |$ b: T/ C        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
2 h3 O% M4 L0 P7 \2 u2 |3 ~        double Cv();
% l% X' b0 b3 w7 R1 Y        double dV_dCA();
2 h; E- X5 }2 V1 h5 k4 }3 ?' N        double dQw_dCA();
        double dQB_dCA();
, M2 |" ^, W6 ]- q2 D7 M; V& _% N        double U();
        double V_CA();
/ b8 G5 I5 T% V/ n/ B        double dU_dnmd();
4 Y  y& T- x( w/ d        double dT_Ys_dCA();
% W4 P2 t3 Y9 X. F7 b        double dT_Pz_dCA();
        double dT_Rs_dCA();
3 O# i- `! W, \' @+ |+ w        double dT_Hq_dCA();
9 v8 R& J  }; K8 r: X1 E/ \$ B        double dT_dCA(double an);
        void Calculate();

/* A編寫各子函數程序*/
        double Cv() /*求定容比熱*/
( v5 m( I! N3 G  ?5 T2 ]' m1 G9 j{

}
3 F0 a) v4 B  {. o' E
- F$ t, I& D; U: F; d
double  dV_dCA() /*求dv/dφ*/
- m5 _3 n' v# O2 |7 I{

! u; i( j1 x; \( X+ T}

$ x2 U) F- c) b; z# x' @$ g1 L
  ]3 V" O9 B/ P4 Adouble V_CA() /*求氣缸容積*/
+ @4 q: n- E* c; n; D2 i) K4 }- t6 y{
5 p2 M3 I, t, j+ r4 P7 ], a
' o0 n6 i4 U4 k* N  \4 b; e}


double dQw_dCA() /*求dQw/dφ*/
0 D1 J! s+ \  i9 @{
, W/ J3 Z. ~7 t$ W0 M3 S  }
}
& \" @. B( r4 ?6 V9 _
/ ~% |" x$ S8 b; ?: F
double dQB_dCA() /*求dQB/dφ*/
% e" d" D& o- ~* n{

}
% p$ p4 Z$ o" v. ^1 F- O' }
4 N( r; K! G1 r5 W: c( U/ }
+ o( [( c" Z1 j/ e# Udouble U() /*求u*/
+ n1 ?; h3 A% z  b4 v8 z+ J{
- b) M6 F( j6 [9 r  q
}


double dU_dnmd()  /*求du/dλ*/
2 [2 p* T1 M- n8 i' n{
8 w/ Y$ b( i: z
7 Q5 k( G: `# q* ~}

1 S  C- |8 D7 B$ H* s" n) ddouble dT_Ys_dCA() /*壓縮期*/
: P; Y; I8 r9 |( h" H4 V{
  `) t5 V" O  |1 O9 Z3 G4 ]
& f$ c+ l% o0 g% B1 @& {' Z2 l$ p2 G}
) M5 P1 P) x6 y7 F: U6 {
8 U, W' d3 i' R) U  ~0 z
double dT_Pz_dCA() /*膨脹期*/
{

* k0 N' M. K/ y5 i9 c, W- F}
0 g7 j% c8 W! w4 v2 s

7 i$ D7 B' F3 `double dT_Rs_dCA() /*燃燒期*/
6 Q: m7 o/ E% x) i- v3 J  P3 E: U{
% i( U0 ^( j6 r* @( ^4 g
}
" A% E$ v  e) u

double dT_dCA(double an) /*求dT/dφ*/
{
        CA=an*PI/180;
! u& n" A% W& B2 u+ H; t$ n2 ^4 A        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();
% k3 |/ Z. H* T        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
        if(A_EO<an && an<=A_SO)    /*換氣過程不計算，線性賦值，避免零值 */
# R& ~' e1 b3 @8 G+ E  ]5 H                P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
0 m1 v; J0 X3 b9 ]& o: T    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);
& z8 C% f4 M: U2 l! J        if(A_EC<=an && an<=A_SC) P=90000;
        return 0;
' Z* Y3 ?  X' l}

void Calculate(void)
1 s, _* V2 y4 k9 X4 T0 E( V$ V{
! [$ X# P) q( `& H    double Tpre,K1,K2,an;
    int index;
        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效壓力bar*/
        mBo=Ne*ge/(30000000*n*gnum);   /*單缸循環供油量kg */
" D0 }0 h/ Z8 Y  n        CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
" S. j) ]3 u( G- R7 H5 z        V=V_CA();/*進氣關閉時角度下的氣缸容積*/
" M0 \: O: p3 q' A        P=P_SC;
" o0 h9 O2 {) q/ f0 x        T=T_SC;
        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;
1 j$ c1 x/ _0 W5 k        an=A_SC+1;

        do{
" A3 {1 p8 y& f% x: f
" D% i. I/ `: t* B& ?! A( F) S5 i3 }  /*B 編寫歐拉法求解常微分方程程序 */

        }while(an!=A_SC+1);
}
% z& ]/ a8 B) S1 I- q( e- |
void main(void)
{
4 f6 A- c8 {" V9 \. c& G        int i;
+ _/ S4 e) ?. Q) N  m        FILE *fp;
        D=0.11;   /*發動機數據按各自的任務書輸入 */
% E  s* c: s  q( x+ R  U        S=0.13;
6 g: r, O$ X3 ^% f# D        L=0.21;
        ge=236.0;
4 M+ {" D" {9 l; T/ ]2 j        n=2400;
4 a6 K" x. ]$ ]: v5 g5 M        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
" a( r2 L& s/ ~! V0 s- V# ^6 n: E# }        Hu=41868000;
% S+ r7 t6 C1 ]% `2 ~' a        R=287.08;
        Pa=101000; /*大氣壓力*/
        Ta=300;    /*大氣溫度*/
        hs=(735.0+R)*Ta;  /*大氣焓值*/
        P_SC=90000;
        T_SC=313;
3 N/ \# R; m, Y; c6 x4 Q5 n
  v, j1 l  p( _! S$ g6 s3 Z        A_SC=216; /*進氣門關閉*/
        A_CB=350; /*燃燒開始*/
        A_CE=410; /*燃燒結束*/
" ^2 p% H  v# y. F8 T+ C        A_EO=484; /*排氣門開啟*/
        A_SO=705; /*進氣門開啟*/
        A_EC=15;  /*排氣門關閉*/

5 O  o0 r$ y6 P  L" E4 H5 u2 U        Lo=14.4;
: I' Z' H5 P  [* ~0 g3 H
! v; F" V7 s& z; e    for(i=0;i<720;i++)
        {
                m_P[i]=0;
                m_CA[i]=i;
6 I& f& ~! ?5 t' |2 h                m_T[i]=0;
4 y* Y$ c1 f3 s$ c7 @+ z  e( q4 q        }
) p( g& J+ F0 e6 H) L2 a9 N
        Calculate();

    /*C編寫數據輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數據到outdata.txt文件 */
# {6 E5 e' {8 D% e+ ~1 v9 F
}


4 u0 L$ @0 f# r" l3 C$ r

祥云飛魚

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

一點一點來

這是作業的干活吧……