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

方進小侯爺

#include "math.h"
0 J( d' o% N2 f/ t* m9 v, S#include "stdio.h"
#define PI 3.1415926
& k  O* {0 Q8 z9 R& d8 U5 N
" k2 T0 a/ B' O4 D; Ydouble D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
! P! d' |5 I9 X; I  t/*缸徑、行程、連桿長、比油耗、轉速、功率、氣缸數、壓縮比、低熱值、理論空燃比、Weibe函數常數m、氣體常數、進氣門關閉時壓力和溫度 */
* {$ A7 M7 N# g2 U5 T0 U        double m_T[720],m_P[720],m_a[720];
6 B$ f3 N! n5 N% [4 D/*氣缸溫度、壓力、瞬時過量空氣系數 */
        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;
5 m' s& J, g% ^* F$ }  I! s" C" p        double Cv();
        double dV_dCA();
        double dQw_dCA();
        double dQB_dCA();
        double U();
        double V_CA();
) A% n! B0 z; T$ a. K        double dU_dnmd();
8 ^* Q" I4 X* O% F; T- e        double dT_Ys_dCA();
        double dT_Pz_dCA();
# X) b/ E- E; Z7 G" `        double dT_Rs_dCA();
' Q  w  z  ^  C$ H' [6 l        double dT_Hq_dCA();
        double dT_dCA(double an);
        void Calculate();
! g$ A+ _  _3 h8 c, T; `8 L, N; J
, }) ^) ]$ m9 \1 T/* A編寫各子函數程序*/
        double Cv() /*求定容比熱*/
" q" J8 ^* ~0 F/ d, I{

) U3 K% i# O2 b& z}

4 ~, \5 @+ O2 Y
/ W: b# Q' G: P4 c% i9 `  U/ odouble  dV_dCA() /*求dv/dφ*/
{

}

% `/ O  ^/ H1 L! m. y
double V_CA() /*求氣缸容積*/
; u3 T8 Z( w- K{

}
) o" u+ x# V2 F1 Y4 M9 c
0 s: K8 j" M9 T7 K3 ?( a
double dQw_dCA() /*求dQw/dφ*/
{
% W* j1 r* _& I' r
# }9 C2 g) A& m% _6 _: E}

4 w7 s3 V3 ?) q, `
' A  u( X, v1 e3 I) z7 ]double dQB_dCA() /*求dQB/dφ*/
{

}


double U() /*求u*/
: m2 j! F3 i4 k% M0 S( d$ _{
5 U# b7 h$ N+ w3 L$ v5 A
" W+ n* R) a1 C: Q}
' F5 y" C" @/ A, }5 w; d# s
( _# Q3 J6 a5 ]' t$ G. R
. @( i4 h! k3 H- c. O5 N# @. Pdouble dU_dnmd()  /*求du/dλ*/
; ?7 {! S7 p" X" G; Q, F# U{
- O$ z5 A( h& m5 T, u0 B/ A( S
}
5 p) U6 m) ]5 |; z. b
double dT_Ys_dCA() /*壓縮期*/
  H" s& h; c2 E% y4 }( r{
  c: v6 ^# s- N9 V6 A' k+ G
}
( ^8 I, g6 y' `- x( V7 `) a
0 [5 x, u+ A+ P% `4 J) N
double dT_Pz_dCA() /*膨脹期*/
+ o  ?; A' a8 ?  O{

}
  ?8 m; g0 F4 I2 |3 J+ l& n

double dT_Rs_dCA() /*燃燒期*/
{
/ ?* {0 i# e1 L, `, j8 M
9 J6 [1 B1 x7 Y# o6 m& }8 c: g( N}

* N3 x3 N# q2 G9 V
0 y5 S5 _! J) O+ K6 Zdouble dT_dCA(double an) /*求dT/dφ*/
: w* C3 R3 e  \& l, D! u{
# B8 F; g# p) }) K" u0 d0 |/ O8 {        CA=an*PI/180;
        V=V_CA();
( G+ u4 B: K4 c5 c        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
0 Z* i* C5 [. @' k1 I        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)    /*換氣過程不計算，線性賦值，避免零值 */
7 n$ s8 h. j4 S3 j- R% \6 f& y7 P" Y                P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
    if(A_SO<an || an<A_EC)
0 H8 A2 j% D4 ~9 e" u                P=110000-20000*(an>A_SO?an-A_SO:an-A_SO+720)/(A_EC-A_SO+720);
' T9 D: O. m  ?/ @" r9 y        if(A_EC<=an && an<=A_SC) P=90000;
        return 0;
  S" d3 p: c( c! W$ y, Q7 H! U}

void Calculate(void)
{
* _" c+ x& w( p    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 */
# i- E' I0 h2 y5 C8 ^5 H3 ?        CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
        V=V_CA();/*進氣關閉時角度下的氣缸容積*/
        P=P_SC;
- j; N4 D' s  G) b* _        T=T_SC;
& {+ j. }3 J+ I- z# Q; r        index=(int)A_SC;
$ Z' G& W: O+ l! N/ c1 \+ A! O7 Q        m_T[index]=T;
, F4 q2 u! V7 }        m_P[index]=P;
# ], I# ~( J7 j4 L. K6 Z5 ^        m=P*V/(R*T);/*進氣關閉時角度下的氣缸內工質質量*/
        mL=m-0.01*mBo;/*進氣關閉時角度下的氣缸內空氣質量*/
    mf=m-mL;
        Tpre=T;
        an=A_SC+1;

* Z& g3 [' [  n' m        do{
  h  V( E/ ^. u6 v0 g
  /*B 編寫歐拉法求解常微分方程程序 */

        }while(an!=A_SC+1);
}
4 `( `& V+ a  P4 T  R: X
- @; S1 p& [' {2 q4 T7 M. Ovoid main(void)
{
  z6 C% v; {) n) }+ {        int i;
8 H3 ^1 ~( ^" h& c        FILE *fp;
# G$ d) @2 |8 R' j& S        D=0.11;   /*發動機數據按各自的任務書輸入 */
- E5 k& c2 k" ?; J. C        S=0.13;
        L=0.21;
        ge=236.0;
        n=2400;
        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
; m+ U% U+ ~7 Q% W/ P" t        Hu=41868000;
        R=287.08;
        Pa=101000; /*大氣壓力*/
        Ta=300;    /*大氣溫度*/
        hs=(735.0+R)*Ta;  /*大氣焓值*/
3 c/ |3 B% c! k' P/ T$ |        P_SC=90000;
        T_SC=313;
2 C5 N! a' r2 z3 w5 w6 Y; Q
        A_SC=216; /*進氣門關閉*/
        A_CB=350; /*燃燒開始*/
* f# x. p7 {- E: N+ A2 X: s% O        A_CE=410; /*燃燒結束*/
' R+ }: \7 I& l' [        A_EO=484; /*排氣門開啟*/
7 _! K2 Q8 L) i5 m4 g; ?        A_SO=705; /*進氣門開啟*/
        A_EC=15;  /*排氣門關閉*/

$ T# l: {1 @. I9 Q3 d        Lo=14.4;

    for(i=0;i<720;i++)
        {
                m_P[i]=0;
) i5 m" L; R1 M4 D7 E                m_CA[i]=i;
3 G7 k& X3 B! k% j7 F* k% |                m_T[i]=0;
9 E" p4 i5 n& b) B4 T        }
; v8 R5 V. d% W1 k: W
% t9 M: u3 l) G2 g$ Z; ^$ F& l        Calculate();

    /*C編寫數據輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數據到outdata.txt文件 */
+ N8 V, y5 Q6 E1 E. M$ G
; i3 @  Q( b0 w+ V}


% k: ?& d! g) R* `% I

祥云飛魚

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

一點一點來

這是作業的干活吧……