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

方进小侯爷

#include "math.h"
#include "stdio.h"
#define PI 3.1415926
5 {% Q6 U, l9 E1 D! j. o' U8 s" D
* N: Q/ Q6 f( J" Kdouble D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
* S! u0 k) G) @/*缸径、行程、连杆长、比油耗、转速、功率、气缸数、压缩比、低热值、理论空燃比、Weibe函数常数m、气体常数、进气门关闭时压力和温度 */
2 h1 P7 r* K$ F$ ~% r        double m_T[720],m_P[720],m_a[720];
/*气缸温度、压力、瞬时过量空气系数 */
, b( V) v6 [0 S7 e" u5 [# T        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();
! S) u" s' u* A        double dV_dCA();
+ b; B/ |$ c$ Q6 X' I& ^; V        double dQw_dCA();
        double dQB_dCA();
3 S, o/ g& `+ _2 Z( X        double U();
        double V_CA();
+ Z6 b6 \+ q; ^. C        double dU_dnmd();
5 B- h2 N, q. e# M  ~/ N: w        double dT_Ys_dCA();
, c9 ^! y4 t. j4 e& x8 r# h6 N        double dT_Pz_dCA();
        double dT_Rs_dCA();
        double dT_Hq_dCA();
, E0 d! u$ F( V0 x( Q, c; [6 \        double dT_dCA(double an);
; l; h5 O  _+ D/ p7 [$ |        void Calculate();
/ [* }' l% x$ a" D8 Q
0 a4 d8 ]  R9 z9 e/* A编写各子函数程序*/
        double Cv() /*求定容比热*/
: s* H  J. n7 f/ s- Q& \( K) f' Q9 e{
! y7 E. `4 r8 G* F
/ A6 p; `. w1 c) g; C) h}

8 N0 h8 M' c7 k6 L+ w" F8 o
7 ]6 w, Q+ M7 M5 ^3 a0 @9 `, b9 |; Kdouble  dV_dCA() /*求dv/dφ*/
8 ?- y! d" R9 G; I; A{
+ G3 N. F" W: s: e& p1 y1 d
}

1 w- I0 J. f- v) ^2 }6 F  w9 w
# p/ I2 v2 Y3 [double V_CA() /*求气缸容积*/
{
* ~, m" A- c6 j% G( O
- a# y! Q+ ]$ ^6 O6 v}

! o9 T7 s' ]: P) H, b  q1 `
double dQw_dCA() /*求dQw/dφ*/
{

) @8 L% Z: X7 g8 F1 o. P}


double dQB_dCA() /*求dQB/dφ*/
{

: U  [) C5 c" m/ V- [) k& O0 `5 t}
/ E3 J1 \  Q$ g3 U/ q1 L) m0 a

double U() /*求u*/
{

}


9 f( `! D3 N' h# e" _) {double dU_dnmd()  /*求du/dλ*/
4 W+ A: S" d, W4 A; W6 R{
% z9 o2 Y0 A. p$ q8 j/ X
. C4 b  x- s& {: m}
  b. Q$ I( {9 _1 t
4 m# k* |+ h; f6 qdouble dT_Ys_dCA() /*压缩期*/
  x" q1 p8 l- C/ v, m' L{

}


double dT_Pz_dCA() /*膨胀期*/
+ r  z1 \* {5 F; J4 \8 B2 J0 g{
7 a. Y2 s; ^- h6 b
! `+ i3 L. P. J  s# `}


$ s' W  |9 ?9 X4 O7 \- ?7 Y% ydouble dT_Rs_dCA() /*燃烧期*/
1 c0 e9 L0 G0 M8 P{
% h, ~8 N' I4 z0 Q% X. k
}
- b2 @% B. b2 E# |7 C8 T$ U

double dT_dCA(double an) /*求dT/dφ*/
{
$ r$ R: p* y* W1 ?7 O  B( A& Q. _        CA=an*PI/180;
        V=V_CA();
$ U' S4 b/ R' l6 Y. z; d        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
! x7 L7 _* s6 p        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
0 M/ i. f5 c! F8 n9 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)
                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;
/ Y! {' \, z  p        return 0;
) P4 D0 \& f. h2 }0 ^}
- i# Q$ R6 t. |  S
1 o" p7 x% O; ^* |/ qvoid Calculate(void)
2 G) W5 I, l, F6 T% {5 k( H$ z{
& _3 K8 v. T& s7 h9 B    double Tpre,K1,K2,an;
    int index;
        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效压力bar*/
7 j0 a! Y: p( u: ]0 @4 P        mBo=Ne*ge/(30000000*n*gnum);   /*单缸循环供油量kg */
        CA=A_SC*PI/180; /*进气关闭时的角度转弧度*/
6 s- O% d, Z3 r0 u% Z) ~        V=V_CA();/*进气关闭时角度下的气缸容积*/
        P=P_SC;
        T=T_SC;
        index=(int)A_SC;
0 N& A* C- R/ ]% P7 }, m3 n        m_T[index]=T;
. f: w/ U: R  S. m- Y        m_P[index]=P;
        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
: U' [7 m) J! ?+ ?$ B        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
    mf=m-mL;
        Tpre=T;
        an=A_SC+1;

        do{

) w% y, \0 D/ o  /*B 编写欧拉法求解常微分方程程序 */
! W4 ~, y. C! f2 I  o
# r( b( }4 t2 |1 Z        }while(an!=A_SC+1);
8 L( }9 Q7 I; e4 M7 X}
- Q% [; f$ s2 i- w& D
+ N2 t. y+ Z! i) E4 D9 X+ dvoid main(void)
6 B9 Y  H+ w- l! c$ V{
3 Y6 F9 P! y1 a; t        int i;
        FILE *fp;
        D=0.11;   /*发动机数据按各自的任务书输入 */
        S=0.13;
        L=0.21;
        ge=236.0;
+ ~6 H; W+ G9 [7 P        n=2400;
        Ne=101568;
% a# }: c% K' a  m8 P4 o, h" B) N        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
2 D" F* s) ]) ?6 p: |4 T        R=287.08;
  N3 `) f3 \4 g4 B        Pa=101000; /*大气压力*/
, P3 h7 N2 k1 Y3 \3 C* U# Q        Ta=300;    /*大气温度*/
# Y3 u, A* o; D1 m+ ^        hs=(735.0+R)*Ta;  /*大气焓值*/
# e! ^0 f! d' b3 \        P_SC=90000;
        T_SC=313;

        A_SC=216; /*进气门关闭*/
7 G& g1 k& t" X! g7 r- K2 U' L4 Y7 \        A_CB=350; /*燃烧开始*/
        A_CE=410; /*燃烧结束*/
        A_EO=484; /*排气门开启*/
% v$ |# u8 q: n+ s1 K; _7 ]: h        A_SO=705; /*进气门开启*/
( \" N7 e# g  r" U6 \- V6 D: {        A_EC=15;  /*排气门关闭*/
# C# F- z& x; v$ e# \
; s7 }6 l; r! m( H" t8 {        Lo=14.4;

    for(i=0;i<720;i++)
, p* o3 y9 I% @( I        {
                m_P[i]=0;
: t4 V3 i& A& k$ H                m_CA[i]=i;
                m_T[i]=0;
3 K. i4 ~6 x1 G        }

        Calculate();
6 n; L6 p2 K. Q) h
    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */

}

# M; c  W9 f! P% N: l

祥云飞鱼

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

一点一点来

这是作业的干活吧……