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

方進小侯爺

#include "math.h"
#include "stdio.h"
9 J3 n) T7 Q( D" }) W% s. q2 P#define PI 3.1415926

* Z2 e5 l; @2 S" b4 D- w  ~double 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];
/*氣缸溫度、壓力、瞬時過量空氣系數 */
/ J( t5 _( t& o! }, [# o        int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
3 n  Q# f' H7 ?. u3 s. K        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
+ n$ ?1 i: t5 y2 y# H        double Cv();
! H1 K4 V' u+ D% E9 L7 z        double dV_dCA();
        double dQw_dCA();
8 K5 d0 X; q) p( I) E) \        double dQB_dCA();
        double U();
% Z6 ]8 n( L1 p/ G( ?* I        double V_CA();
- f" t5 S6 e$ C7 c/ l6 H: f8 _' F        double dU_dnmd();
        double dT_Ys_dCA();
1 s) }& q6 P( a        double dT_Pz_dCA();
& W: L" F% n( |, m& M0 D        double dT_Rs_dCA();
        double dT_Hq_dCA();
" V9 M2 t. Q4 ]( F$ ~3 Z' z9 M( D        double dT_dCA(double an);
: ]' Z8 O; m+ [; K) ~, b! l' L# Z        void Calculate();
' Y- [6 E0 ^3 R. D) g3 |
/* A編寫各子函數程序*/
" a2 X7 z7 t* ^% S4 U        double Cv() /*求定容比熱*/
{
; Y4 f3 C& c$ S* H; `6 e5 |
2 n/ J2 S; i8 S8 p}
5 ^7 G2 m4 ?! y5 `
& J4 Z* B7 Q% W/ ^( P0 {
. z7 ?' A  C+ O8 ~  C, Y% tdouble  dV_dCA() /*求dv/dφ*/
{

! J, I9 B% L, p$ L0 _: T) Z+ i}


$ r! \( y# A2 y# g, Y3 B2 [double V_CA() /*求氣缸容積*/
" {9 h8 C, K% C5 |$ `& ?{
0 K1 X% Z7 ?, s+ I
( m* [, v# u2 e1 |}

2 \& [9 E! ^2 Y# \0 Q  e
double dQw_dCA() /*求dQw/dφ*/
{

}
2 t* D' R7 V* D$ t9 n. S
( \+ Y2 v) d0 X- N
double dQB_dCA() /*求dQB/dφ*/
{
& q+ r: v5 G( S6 i9 H
& \8 `+ E  Z) p: Y" W}

6 \0 ^5 ~0 Z  p! q! p
; \) o/ w: a6 c$ o# ldouble U() /*求u*/
3 b- h; E6 u) ?; e* {- R{

  [6 L0 g. V# G3 @' E2 \, e}


' ]3 R  Z4 q" a$ Y+ |; @- hdouble dU_dnmd()  /*求du/dλ*/
7 ~6 @7 H- u! x1 c  s: P{

}

double dT_Ys_dCA() /*壓縮期*/
* _$ k+ N2 Y# t' Q+ I2 _3 q% O) ^{

}
. @3 n7 c) Q2 t/ r6 z4 O9 B

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

% z' l! ]+ @, K8 q$ c8 r}
) _) L8 G" a4 I9 D- |+ x+ I3 Z

double dT_Rs_dCA() /*燃燒期*/
{
, ]/ y* f# T. c( f" c
, h, C/ u4 @0 x( b% S}
) X4 b/ |1 h$ r( L" O/ t

6 z& e- D9 V$ z' O6 Hdouble dT_dCA(double an) /*求dT/dφ*/
0 d( L% T, z3 q8 e7 n{
        CA=an*PI/180;
        V=V_CA();
' L) A# a- h! i1 F1 g        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
% K& B9 I1 Y  f1 i        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
6 h( m  q- @% a3 V% P4 w% |        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)
# n, q* _2 G+ y3 b2 s- [% s5 Y! a                P=110000-20000*(an>A_SO?an-A_SO:an-A_SO+720)/(A_EC-A_SO+720);
3 [1 S6 e* X3 v+ s  g3 Q3 y( `        if(A_EC<=an && an<=A_SC) P=90000;
  r6 P* c) a, r, W        return 0;
}
& B6 z- @8 l# N3 \& u
. E8 v4 N% D+ r8 E) Y5 n) x3 Gvoid Calculate(void)
{
    double Tpre,K1,K2,an;
7 c8 t* [) S6 Y" E    int index;
3 B! ]. B7 F3 ~0 y+ A0 O        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效壓力bar*/
  J. o, a2 `+ E/ G4 v) N0 S        mBo=Ne*ge/(30000000*n*gnum);   /*單缸循環供油量kg */
1 O( B. ^5 @: I        CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
        V=V_CA();/*進氣關閉時角度下的氣缸容積*/
( \! y, T5 \) e5 ^8 t7 w4 H0 m# q        P=P_SC;
        T=T_SC;
+ a& `& Q- L% {/ G' ?7 ~# e        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;
  \$ _+ ~# V( @6 _5 ^. ]        Tpre=T;
        an=A_SC+1;
) `  P! N# t8 X. T3 U; ]# D
' v! o- Q$ f1 [* Z) b- _        do{
6 K2 I& K& b& P* s( R8 R9 O9 k
  /*B 編寫歐拉法求解常微分方程程序 */
! M2 l/ ]; X+ U, o
        }while(an!=A_SC+1);
}

! W1 X$ }& H3 cvoid main(void)
{
        int i;
& ^" a/ T& c0 T7 b, C        FILE *fp;
" M0 }& J* K! ?* E& G2 G4 [4 U        D=0.11;   /*發動機數據按各自的任務書輸入 */
        S=0.13;
/ `/ P  S! R( I8 k) J3 |9 f        L=0.21;
        ge=236.0;
        n=2400;
1 x0 D. V% O% l! C. v+ |3 B' }        Ne=101568;
7 B: n$ h( D, ^& V7 x5 y        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
- p. T5 u$ K. D. R: X* D3 ?: K9 n& P- T        R=287.08;
( W) Q# M$ }$ s        Pa=101000; /*大氣壓力*/
. [( [6 i1 A5 b% {7 o0 p+ X        Ta=300;    /*大氣溫度*/
        hs=(735.0+R)*Ta;  /*大氣焓值*/
        P_SC=90000;
3 H& `- Z$ C6 a% |        T_SC=313;

' e/ `+ V6 h# T$ }5 p        A_SC=216; /*進氣門關閉*/
        A_CB=350; /*燃燒開始*/
        A_CE=410; /*燃燒結束*/
/ p: `; {1 u; w; x& A2 @% J        A_EO=484; /*排氣門開啟*/
        A_SO=705; /*進氣門開啟*/
        A_EC=15;  /*排氣門關閉*/

        Lo=14.4;

    for(i=0;i<720;i++)
( V; i0 u* z$ ?+ z( C% B* f        {
2 y# s, T- `! F! y: G                m_P[i]=0;
7 O: H6 p$ b/ _0 l8 Y5 W$ y                m_CA[i]=i;
                m_T[i]=0;
  s6 w) f; N' P2 p' `        }
8 Q5 n0 J6 D+ _' w6 C: v! M7 n
# M1 V1 B5 g+ X' Q; f' G2 c        Calculate();

    /*C編寫數據輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數據到outdata.txt文件 */

3 z5 i3 [! T' Q3 }}
4 j! r* X1 V% `3 ]9 `

祥云飛魚

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

一點一點來

這是作業的干活吧……