久久久国产一区二区_国产精品av电影_日韩精品中文字幕一区二区三区_精品一区二区三区免费毛片爱

機械社區

標題: 內燃機熱力計算怎么用C語言編程？模板如下，求助。 [打印本頁]

---

作者: 方進小侯爺    時間: 2015-6-28 22:58
標題: 內燃機熱力計算怎么用C語言編程？模板如下，求助。
#include "math.h"
#include "stdio.h"
#define PI 3.1415926

. N( O. G. j# X. n& Wdouble D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
& x- }3 s4 O0 H$ _- _  r8 B/*缸徑、行程、連桿長、比油耗、轉速、功率、氣缸數、壓縮比、低熱值、理論空燃比、Weibe函數常數m、氣體常數、進氣門關閉時壓力和溫度 */
* z* `& i) Y: f1 e: N) c+ [        double m_T[720],m_P[720],m_a[720];
/*氣缸溫度、壓力、瞬時過量空氣系數 */
  U( Y' A: u0 w& @; R3 s        int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
% p; [- q- F) O! x1 s        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
        double Cv();
        double dV_dCA();
        double dQw_dCA();
! {1 V+ v" i) X1 ]* F; L% v        double dQB_dCA();
9 y+ M  M4 m9 k8 r1 c5 h& g0 V4 a# ^* q        double U();
8 w4 L! z+ @+ V) P; ~9 }        double V_CA();
        double dU_dnmd();
6 Y. w9 R- W& {$ z5 S# g! `        double dT_Ys_dCA();
        double dT_Pz_dCA();
        double dT_Rs_dCA();
/ g1 x( T) V% y1 j  j9 o        double dT_Hq_dCA();
        double dT_dCA(double an);
$ X6 j) q8 z0 f9 y( ]        void Calculate();
( j# A% I' W% M% q5 Y" L9 l, @( w6 v
1 V+ e8 p, a* M( y( g; V/* A編寫各子函數程序*/
        double Cv() /*求定容比熱*/
' m9 K. _$ k. i+ o. d& c{

}

! T( D7 k, E$ D9 I7 x
double  dV_dCA() /*求dv/dφ*/
4 U3 ?" j6 C' R0 ?* R8 q{
; X1 P# Z/ J5 A$ \+ m# S3 h
}

0 w8 t7 p9 X( D# U* Q& W8 Z
2 o2 C3 r3 j$ b% I& W, f6 N) Adouble V_CA() /*求氣缸容積*/
# s- Q7 \# d$ }' K' E7 J{
5 U3 o1 Y( z# P4 ~
}


3 u8 N3 Z  X+ ~7 ^double dQw_dCA() /*求dQw/dφ*/
{
3 A" }7 j/ z& v# k: J
}
1 Z7 X& a7 A4 \# c$ p) c& t: E
3 V' X& [% r! g) n, ^4 E. b6 P! e
double dQB_dCA() /*求dQB/dφ*/
3 n; c( ^( [* a1 G{
+ }& [* E& R5 l5 x
2 `7 Q6 P" F/ o! c7 K9 s# T}

2 S! \/ `- s5 r* z. F0 U
6 W1 B& e5 t( T1 v. ]  odouble U() /*求u*/
{
  f8 g7 S- l1 ~. O8 j9 b% Q
}
( y. Q& O! D' e0 C% w2 D
/ Y5 `' j+ v% T5 a6 I* _2 C4 |4 i* K
double dU_dnmd()  /*求du/dλ*/
{

}

' \# E, w# Y) m& pdouble dT_Ys_dCA() /*壓縮期*/
5 o6 f; P, X+ }: C# F0 r1 w* n' \{
' ]" J/ `9 v/ X6 I
8 v! q! Z! q1 X1 Q}


double dT_Pz_dCA() /*膨脹期*/
: t3 p7 U! i2 @) W+ V% V{

. ^' E  V0 Y5 R8 r" m; L6 o}
. G1 T* i: a6 R1 N; S% W" C

double dT_Rs_dCA() /*燃燒期*/
{

# ]7 j2 v+ M  m8 U}
  ~1 m% j2 r+ w8 |+ p
; y) N# Z) u1 k& l2 ?2 T( }8 c
double dT_dCA(double an) /*求dT/dφ*/
9 t$ f1 h9 f" j! w  H6 Y! A* |{
2 e1 T6 Z7 r% ~4 M/ |% f        CA=an*PI/180;
( R& t0 C2 H/ F1 s: Q7 ~6 b! G        V=V_CA();
- r+ A+ B" w0 C2 Q        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
- p9 f6 m3 E, T7 y% ?8 H" U' h        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
' @; A) J! u/ g  e6 o& ]. |' x& U        if(A_EO<an && an<=A_SO)    /*換氣過程不計算，線性賦值，避免零值 */
                P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
# d5 T: Z  x7 k1 I+ B4 R. E    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;
- b* w! E  W$ L, L2 u        return 0;
: {5 `! N( s' C6 K$ P9 ~% n}
1 ], N$ S& R6 Y: C' h# k8 z
void Calculate(void)
1 Y' M7 w9 o. I/ v. O" x2 \  X9 [{
    double Tpre,K1,K2,an;
( c1 j* V6 y2 {0 h% w3 P    int index;
# T" n2 z- s# ?        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效壓力bar*/
        mBo=Ne*ge/(30000000*n*gnum);   /*單缸循環供油量kg */
        CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
        V=V_CA();/*進氣關閉時角度下的氣缸容積*/
7 C2 b6 t! f# z$ S        P=P_SC;
9 D, f* x0 s0 @' J        T=T_SC;
" U9 |- ^/ a% i2 W; ~/ R        index=(int)A_SC;
        m_T[index]=T;
) f; t/ J6 e6 n" i% G) o- K        m_P[index]=P;
        m=P*V/(R*T);/*進氣關閉時角度下的氣缸內工質質量*/
7 h' T" \) c+ s1 H, \        mL=m-0.01*mBo;/*進氣關閉時角度下的氣缸內空氣質量*/
    mf=m-mL;
        Tpre=T;
        an=A_SC+1;
4 E: p0 @8 i" v3 v, i1 H
8 L  N8 D) y* c% {3 I1 z        do{
" F& ^$ z$ ?# R2 ]9 [% t5 m
  /*B 編寫歐拉法求解常微分方程程序 */

        }while(an!=A_SC+1);
}

void main(void)
{
        int i;
        FILE *fp;
        D=0.11;   /*發動機數據按各自的任務書輸入 */
        S=0.13;
. k4 u2 ~3 M! p: r  N  w        L=0.21;
: Z& d* Z" w( M  a4 d/ `        ge=236.0;
        n=2400;
& E& y; D+ {( w+ H& {6 u3 l        Ne=101568;
        gnum=6;  ysb=17;   M=0.7;
        Hu=41868000;
# a: P* h9 x4 u* |* Z* T) \6 e7 }        R=287.08;
        Pa=101000; /*大氣壓力*/
' ?' N" }' B. i        Ta=300;    /*大氣溫度*/
        hs=(735.0+R)*Ta;  /*大氣焓值*/
+ q% ^) H% k  i6 A5 K        P_SC=90000;
        T_SC=313;

5 x- b' `. W3 t, Z        A_SC=216; /*進氣門關閉*/
& n& B1 b/ L+ P3 J        A_CB=350; /*燃燒開始*/
* T2 z3 b# v! V7 I7 A        A_CE=410; /*燃燒結束*/
        A_EO=484; /*排氣門開啟*/
        A_SO=705; /*進氣門開啟*/
+ j5 n/ ~; e( s* k6 z        A_EC=15;  /*排氣門關閉*/
8 y; j, N* D0 L6 q* f( N
4 ]+ p/ d" b2 A* E, X# E7 o' [        Lo=14.4;
# m( A& G' q; V+ @0 \) c
4 m5 r: D+ `/ ^/ f/ R3 q4 N    for(i=0;i<720;i++)
        {
9 _5 R' _: V0 W4 a9 n& f                m_P[i]=0;
3 W; f9 B+ F  j, V7 F- l& m                m_CA[i]=i;
! q4 ]6 h& B6 F/ a/ p                m_T[i]=0;
: @! [0 j8 l8 s0 d        }

' B7 {* m- f2 ^$ ]5 f" V3 c        Calculate();
. {. s/ S; b( |' Y7 Q, z7 X, O
    /*C編寫數據輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數據到outdata.txt文件 */

}
6 l& M; [. Z  Y/ _3 d) J- |0 N

' T7 T0 O( J1 K$ V

---

作者: 祥云飛魚    時間: 2015-6-28 23:45
不難啊  完全是 把 公司 翻譯成 C 又沒有用指針跳轉的。。。

---

作者: 一點一點來    時間: 2015-6-29 01:45
這是作業的干活吧……

---

歡迎光臨 機械社區 (http://www.ytsybjq.com/)

Powered by Discuz! X3.5

感谢您访问我们的网站，您可能还对以下资源感兴趣：

亚洲影视一区二区三区

<li id="yagss"><tbody id="yagss"></tbody></li>

<ul id="yagss"></ul>

<input id="yagss"></input>

<sup id="yagss"><xmp id="yagss"></xmp></sup>

<rt id="yagss"></rt>

<em id="yagss"></em>

<dl id="yagss"><center id="yagss"></center></dl>