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

方进小侯爷

#include "math.h"
  U) ~* y% P4 Y' s% L#include "stdio.h"
#define PI 3.1415926
% L) y# W7 W* ^
9 L3 u. q% X: b! K4 bdouble D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
& l; k5 A  a2 k3 M, B/ ?/*缸径、行程、连杆长、比油耗、转速、功率、气缸数、压缩比、低热值、理论空燃比、Weibe函数常数m、气体常数、进气门关闭时压力和温度 */
        double m_T[720],m_P[720],m_a[720];
/*气缸温度、压力、瞬时过量空气系数 */
        int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
  H# y' l; Q& ]9 v/ d/ B$ I& i- A        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
        double Cv();
% R' U2 I1 `' o; A- w7 E" x        double dV_dCA();
        double dQw_dCA();
$ j/ S# i+ ~1 S5 r( O        double dQB_dCA();
$ r4 T6 T- D, k( ]        double U();
- \8 K, p1 Z/ N% m        double V_CA();
1 J1 d* ^1 X; F+ D" x        double dU_dnmd();
        double dT_Ys_dCA();
0 g% _0 p/ t- }6 L: t; f- \' j        double dT_Pz_dCA();
        double dT_Rs_dCA();
        double dT_Hq_dCA();
        double dT_dCA(double an);
. ^* W2 U3 h! m( G( n/ w  u        void Calculate();

/* A编写各子函数程序*/
        double Cv() /*求定容比热*/
2 f) Y; G* i! n0 I{

# |5 c1 {9 R- C- O5 e}
% u! c3 c- X. u5 Z
3 J, E$ ~9 A% n* F# }; Q$ A
double  dV_dCA() /*求dv/dφ*/
{

$ X: j9 ?; [- @8 @}


( ^0 q' r0 G) T6 F0 V) A9 i2 b$ Ydouble V_CA() /*求气缸容积*/
3 `) Z4 H9 T0 r* x{

# H5 R6 Z2 s2 F& F. {& i% l3 U% @}


3 y* U( N) H# Z9 {3 ?8 adouble dQw_dCA() /*求dQw/dφ*/
) M1 }: o, P& A{
/ K; i  B3 l9 u. u
}
1 p+ H3 Q- L; G+ P3 Z0 {4 u! e

+ U) V- Y* |. h+ m" s6 [double dQB_dCA() /*求dQB/dφ*/
; \4 A+ c* l' D: o{

}

" c* K; E- c% v. F7 Y% v  H, [
double U() /*求u*/
{

! \' P2 h) }* ~. C9 R" e: V, Z. }}
' W0 d; T. K2 x/ ?  {

double dU_dnmd()  /*求du/dλ*/
/ r+ Y: f8 k  C6 }( u- _9 a{

0 t! p7 w% ?4 t/ }, B& L* v  Z+ X}
: x; B; [( G3 D3 X$ b; |" ~0 X
double dT_Ys_dCA() /*压缩期*/
+ H4 N- y. }& ~{

}

, ^5 s3 I2 i5 z! h# d4 L* V
double dT_Pz_dCA() /*膨胀期*/
{

}
% ?" s' ]5 O/ l

double dT_Rs_dCA() /*燃烧期*/
{
. @$ N, C/ W; k/ `& D) L
$ j1 E; D# ?9 r8 C+ |% Z}


double dT_dCA(double an) /*求dT/dφ*/
{
        CA=an*PI/180;
        V=V_CA();
, ]3 G7 }, z: D8 m        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
% m% k# A$ i6 P/ h  _! ^        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
! b9 X  _* `, Y' }2 S) h        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
: ~# ]) U! s) B7 r        if(A_EO<an && an<=A_SO)    /*换气过程不计算，线性赋值，避免零值 */
                P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
! b/ G1 n9 Q( w    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;
1 k, s5 n! P9 q# B        return 0;
}

void Calculate(void)
( @6 z- M& f) t  ]8 A{
    double Tpre,K1,K2,an;
' n6 c3 ?) p8 z# S1 Y$ L" O    int index;
% n" j  y# ?$ N( N$ C6 M8 j* ]( D        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效压力bar*/
! Q) [% t, N" W0 {7 L        mBo=Ne*ge/(30000000*n*gnum);   /*单缸循环供油量kg */
        CA=A_SC*PI/180; /*进气关闭时的角度转弧度*/
        V=V_CA();/*进气关闭时角度下的气缸容积*/
        P=P_SC;
        T=T_SC;
        index=(int)A_SC;
        m_T[index]=T;
- f% m/ p  D- M        m_P[index]=P;
# \: a- ], ]; U2 d5 r# ?        m=P*V/(R*T);/*进气关闭时角度下的气缸内工质质量*/
        mL=m-0.01*mBo;/*进气关闭时角度下的气缸内空气质量*/
& v$ d4 Z: p- e# T( r7 E    mf=m-mL;
        Tpre=T;
        an=A_SC+1;

        do{
! R5 p/ S, h- s. _+ n
. Q' U& ?: m" R9 F4 {$ z2 @. s% s  /*B 编写欧拉法求解常微分方程程序 */

. M# |. `/ D3 g6 Y/ M* l( R        }while(an!=A_SC+1);
: ^2 G$ s$ [6 t8 n' r, Q0 k}

void main(void)
{
# |+ e  h  h5 R$ p) m        int i;
        FILE *fp;
1 E+ c) q% a6 a; c( ]        D=0.11;   /*发动机数据按各自的任务书输入 */
7 D/ A4 @$ q- V9 W2 k) k        S=0.13;
' J( |; I2 x+ x        L=0.21;
        ge=236.0;
        n=2400;
        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
' {. `! |; m$ i, Z2 e: T        Hu=41868000;
) {, {: g2 G" v5 c        R=287.08;
% ^) u5 \7 O$ ]+ X+ p        Pa=101000; /*大气压力*/
        Ta=300;    /*大气温度*/
7 d. }. o2 F2 A2 z* [3 }        hs=(735.0+R)*Ta;  /*大气焓值*/
* K0 |4 T' v/ {* {3 M8 n        P_SC=90000;
        T_SC=313;
" A* m, h* ^, x& O
        A_SC=216; /*进气门关闭*/
        A_CB=350; /*燃烧开始*/
        A_CE=410; /*燃烧结束*/
. M4 Y8 \6 ^, M6 E: r        A_EO=484; /*排气门开启*/
$ Z' N/ |9 q* @+ O) t6 Z: L        A_SO=705; /*进气门开启*/
        A_EC=15;  /*排气门关闭*/

9 n* e: E9 W3 S; O1 U# y' e) o        Lo=14.4;
$ Q) S2 T" i" q$ U* F
. ]. B4 K' ~* f% G# p    for(i=0;i<720;i++)
% A2 N* y+ X+ k0 n( ^        {
                m_P[i]=0;
                m_CA[i]=i;
                m_T[i]=0;
% @+ y$ _' K* T4 p/ F$ D! Z; X        }

        Calculate();
1 h4 b1 f1 Z  N# \* S! h
    /*C编写数据输出程序，计算结束后按顺序输出0-720转角的“缸压，温度”数据到outdata.txt文件 */

7 t( d; u, O% b2 H( `}

祥云飞鱼

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

一点一点来

这是作业的干活吧……