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

方进小侯爷

#include "math.h"
#include "stdio.h"
; M) a- ]4 ^1 y7 ?6 l( |3 ~#define PI 3.1415926
9 [; g) o1 B' x5 U  c4 ?
" b- E7 F1 D" x, E5 C: Kdouble 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];
9 X& K1 C, E) n" Y. S( l" K& R+ V/*气缸温度、压力、瞬时过量空气系数 */
        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();
        double dV_dCA();
        double dQw_dCA();
        double dQB_dCA();
; _0 W* D  t- m8 P9 D        double U();
; c! Z: y8 s* `& v/ G! F        double V_CA();
# m- e' M; W# X9 k) g6 C        double dU_dnmd();
) M- p. o# m) Y6 _8 `5 v5 |9 m3 E        double dT_Ys_dCA();
        double dT_Pz_dCA();
% C& Y  B3 Y& K9 z3 z        double dT_Rs_dCA();
        double dT_Hq_dCA();
        double dT_dCA(double an);
3 Z( |' }7 |9 @* \# g        void Calculate();
5 c: X/ Q$ g  G1 i
/* A编写各子函数程序*/
        double Cv() /*求定容比热*/
: k2 V! Y. D& u) X9 B$ D( o. Q) ~{

6 ]$ }/ O0 f7 A- D& X2 a% U}
' M- H  Y: h) `9 U

# I' C- s. q9 edouble  dV_dCA() /*求dv/dφ*/
{

}


double V_CA() /*求气缸容积*/
4 {0 e: V" K& K0 t* X6 U+ d9 F{

}

+ G( A1 u8 V% i4 e& [$ X
double dQw_dCA() /*求dQw/dφ*/
{

- L  s' |4 t4 N9 q" p}
2 A: k7 {( g. m4 p

double dQB_dCA() /*求dQB/dφ*/
{
, @( `& r1 Z: d
- W& t6 F: h4 i) b; r$ s+ Y}
* r. O; N! \6 y% D+ @0 p' C
: C" P2 w! ?2 \& k5 J% B& W( X
double U() /*求u*/
{
1 ^! {+ x6 d: G2 C% p
0 l' {# u) ]3 P, y- O( ~' A. x}
% B1 `6 x8 `2 l- d' ^

* v* Q3 n- ^7 ^double dU_dnmd()  /*求du/dλ*/
{
+ ^, {! V, A. i% ^- h
# r1 Y- I3 e! j2 |* h3 z: i}
: T9 ^% g. g1 n, y: Z
double dT_Ys_dCA() /*压缩期*/
{

  M4 R! `4 Y5 Y% B7 S% b}


double dT_Pz_dCA() /*膨胀期*/
{
( x) M/ Z. i$ ]/ t' ~
}


( N: U8 i- c  Ldouble dT_Rs_dCA() /*燃烧期*/
{
. N9 O5 w% ^' P; I6 Y# @& U8 X/ [
; F3 G* |( T# T; r0 ?' r}
: J  X6 N8 N0 f; @
* L- T% ]. Y3 I! R1 C  i
double dT_dCA(double an) /*求dT/dφ*/
: w' O% d" p5 `. Y# C8 s" V; w{
, j8 M: O8 I7 T# S. f8 C        CA=an*PI/180;
' u  B, r6 r% U  o2 K( C8 |        V=V_CA();
        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
2 P! ?7 \4 W6 C+ ^0 [, D        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
8 L5 Q2 X' p0 Z' W9 B! c- M3 y5 M        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
        if(A_EO<an && an<=A_SO)    /*换气过程不计算，线性赋值，避免零值 */
/ m* }/ X( d: c                P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
    if(A_SO<an || an<A_EC)
3 }5 W! ^! B' ^/ J- ^8 M                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;
) X; Z% @$ K3 M' _& e5 }/ Y$ @5 ?        return 0;
}

void Calculate(void)
9 _4 o. b: `7 d+ Q# W2 ]3 I{
    double Tpre,K1,K2,an;
    int index;
        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效压力bar*/
* Z: `/ ~& U6 w        mBo=Ne*ge/(30000000*n*gnum);   /*单缸循环供油量kg */
- ~* R' V/ G3 w& `# k$ o( u        CA=A_SC*PI/180; /*进气关闭时的角度转弧度*/
, s1 _5 z9 w6 m/ V8 O- M; g        V=V_CA();/*进气关闭时角度下的气缸容积*/
        P=P_SC;
8 r' m' a5 c; ?* p' s        T=T_SC;
& w+ y) m7 C/ k/ A: C7 K        index=(int)A_SC;
% J$ ], A: z* E$ E        m_T[index]=T;
4 l: }8 ?& X! R. g) o  f        m_P[index]=P;
2 Q6 G: S) x2 K4 n1 R        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
4 o! A, e) t) H/ Z* G+ t8 w        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
    mf=m-mL;
        Tpre=T;
6 _. s' |6 C! F9 x2 I+ q        an=A_SC+1;
% ]' V& S: x+ Y+ \  x
. Q+ |- g. j4 I& d$ s8 q2 }6 a        do{
3 r- H2 g5 q1 S$ U1 s
  /*B 编写欧拉法求解常微分方程程序 */
# P9 ~3 P: ~9 ]2 T
        }while(an!=A_SC+1);
}
: u# C+ e1 x# C5 y
, x: v, d4 r# c" [1 r+ m$ k# Z8 pvoid main(void)
{
  M) x/ M7 E1 e9 n        int i;
        FILE *fp;
* V; C  z' v" U1 I. z" {: @5 @        D=0.11;   /*发动机数据按各自的任务书输入 */
        S=0.13;
( a% J: A1 i% k! G0 s2 F2 m# X4 T        L=0.21;
        ge=236.0;
        n=2400;
" \/ m  Q* t7 j. y$ S        Ne=101568;
- b3 X, q7 z- u/ i' W        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
# f0 g3 S* a+ G, V        R=287.08;
/ A: r2 m) }; k. O7 Y        Pa=101000; /*大气压力*/
        Ta=300;    /*大气温度*/
2 i( a  C4 B: l9 {8 Y; u9 D        hs=(735.0+R)*Ta;  /*大气焓值*/
4 Z6 K7 ?  f2 m8 N        P_SC=90000;
        T_SC=313;
' L4 l. |" D/ ^7 M3 t0 N
8 q6 j3 F5 ^, |1 M" C% P        A_SC=216; /*进气门关闭*/
% G/ p0 W0 H# D6 D! V, W        A_CB=350; /*燃烧开始*/
. T  t2 \3 w& e; t1 j  V5 h: K: r        A_CE=410; /*燃烧结束*/
0 _& o' ?4 R  k, {- p        A_EO=484; /*排气门开启*/
        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/
' q7 S2 r0 ~2 g  R
        Lo=14.4;

: R/ Y; E! w' L) n! ?$ |    for(i=0;i<720;i++)
; d6 g- S5 d3 _9 |        {
% g/ g0 e" w' {& _" N8 o) m, W                m_P[i]=0;
                m_CA[i]=i;
, V1 h- l& V! }0 q+ a& e: \( F                m_T[i]=0;
        }
  R/ i4 e9 }; y7 k" R: `1 {
        Calculate();

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

}


6 R" t4 m2 c3 b5 G6 G- \

祥云飞鱼

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

一点一点来

这是作业的干活吧……