国产精品乱码一区-性开放网站-少妇又紧又爽视频-西西大胆午夜人体视频-国产极品一区-欧美成人tv-四虎av在线-国产无遮挡无码视频免费软件-中文字幕亚洲乱码熟女一区二区-日产精品一区二区三区在线观看-亚洲国产亚综合在线区-五月婷婷综合色-亚洲日本视频在线观看-97精品人人妻人人-久久久久久一区二区三区四区别墅-www.免费av-波多野结衣绝顶大高潮-日本在线a一区视频高清视频-强美女免费网站在线视频-亚洲永久免费

<blockquote id="dsql3"></blockquote>

<blockquote id="dsql3"></blockquote>

機械社區(qū)

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

---

作者: 方進小侯爺    時間: 2015-6-28 22:58
標題: 內燃機熱力計算怎么用C語言編程？模板如下，求助。
#include "math.h"
0 x& s/ v* G2 B8 C7 ^1 u#include "stdio.h"
' K$ w. i5 I0 |6 t#define PI 3.1415926
: v6 ?7 r) \, }  x; M! r% X9 x
double D,S,L,ge,n,Ne,gnum,ysb,Hu,Lo,M,R,P_SC,T_SC;
' z1 K. l1 I- R! ?5 ?0 Y/*缸徑、行程、連桿長、比油耗、轉速、功率、氣缸數(shù)、壓縮比、低熱值、理論空燃比、Weibe函數(shù)常數(shù)m、氣體常數(shù)、進氣門關閉時壓力和溫度 */
        double m_T[720],m_P[720],m_a[720];
/*氣缸溫度、壓力、瞬時過量空氣系數(shù) */
7 W! m) @+ c6 t& q        int m_CA[720],A_SC,A_CB,A_CE,A_EO,A_SO,A_EC;
! j: p7 i; e0 v# s5 P        double Pa,Ta,hs,P,V,T,CA,m,mf,nmd,Pe,mBo,mL,dXY;
( _3 M" T2 R0 i        double Cv();
" k# Z$ h1 y2 `        double dV_dCA();
        double dQw_dCA();
7 q- W: X2 [! g        double dQB_dCA();
        double U();
        double V_CA();
0 \+ j' s0 u, h        double dU_dnmd();
) Y  l* y7 ~' ^( e& `7 V* D        double dT_Ys_dCA();
        double dT_Pz_dCA();
6 M& ^/ i4 ~, v! D8 l8 d4 {        double dT_Rs_dCA();
. c6 s  r6 u. `% ^  u) |        double dT_Hq_dCA();
        double dT_dCA(double an);
        void Calculate();

+ _. q8 G6 F* J: |+ K- Z1 l/* A編寫各子函數(shù)程序*/
        double Cv() /*求定容比熱*/
{
4 \) O# [6 ]6 L6 C* A& X
: Z; m2 V! I( X# F}
0 Z8 v7 F( K" J  H0 G4 u, N
& K8 Z# Y/ B2 g5 o3 T
double  dV_dCA() /*求dv/dφ*/
7 _; @) q) R# \+ S6 ]6 R8 X) I+ y7 B4 d{
9 w' X" H* j' K5 s7 E0 e
}

; ]( R$ b/ `$ I0 ?3 C
; o- }) y$ K- b$ j: z1 n! W  Odouble V_CA() /*求氣缸容積*/
{

; f2 s; ]$ z) l4 J}

+ z- t* {3 X1 n# x
double dQw_dCA() /*求dQw/dφ*/
{

}
" y- ?1 E' N6 M( q$ b
& ]8 C' i+ G3 m$ q  m6 D; ~
double dQB_dCA() /*求dQB/dφ*/
{

}

4 z4 l6 R* }: k( p0 b
double U() /*求u*/
{
( ]. [& @9 f( ~* H1 a
# K# J; U# `1 S1 J}


double dU_dnmd()  /*求du/dλ*/
{
* b& t1 u- o5 ~4 p7 B0 u. V9 S
}
6 A- d5 m) ]5 K9 B; L
/ g9 `  h) c. |double dT_Ys_dCA() /*壓縮期*/
* y% m  _1 d  m4 y2 G# j{
/ f. x8 n. x9 u6 m& C6 N% a
}


7 |. b9 ]6 P" Z( I8 J6 x' B4 s5 g3 |double dT_Pz_dCA() /*膨脹期*/
- ]( z9 D  m+ q* T# J{
) \/ d& R) c( X
4 k* u3 z# f2 S/ x8 S& K}
6 G- }# G5 r$ M' C2 p6 N* r+ [4 G# C

double dT_Rs_dCA() /*燃燒期*/
$ c# p- ^: ]  A9 n4 r{
( U/ O  v7 f. Z- L, U- t7 O- U  V
}
6 _9 R; }3 r# P* t: U- [3 M

$ b2 \7 w5 m- b1 ]# ndouble dT_dCA(double an) /*求dT/dφ*/
% n& y% g& k. a. y! f4 `% p{
. O% s* x7 T2 w  g        CA=an*PI/180;
        V=V_CA();
: K% M! H$ M# B" G$ t6 V        if(A_SC<an && an<A_CB)   return dT_Ys_dCA();
        if(A_CB<=an && an<=A_CE) return dT_Rs_dCA();
        if(A_CE<an && an<=A_EO)  return dT_Pz_dCA();
        if(A_EO<an && an<=A_SO)    /*換氣過程不計算，線性賦值，避免零值 */
+ q8 x2 Q" W& g1 J8 q7 o$ b                P=110000+(m_P[(int)A_EO]-110000)*(A_SO-an)/(A_SO-A_EO);
& e6 R6 p2 ~$ v) c1 [    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;
        return 0;
}
9 j7 q  V* _' f9 p! b* n6 n
void Calculate(void)
- E' Z$ D: j3 ~! ^4 ~3 u& T- t. y{
    double Tpre,K1,K2,an;
7 S6 \' D/ R6 [4 ]( u! t& W/ D    int index;
        Pe=30*4.0*Ne/(gnum*100000*n*S*PI*D*D/4);  /*平均有效壓力bar*/
        mBo=Ne*ge/(30000000*n*gnum);   /*單缸循環(huán)供油量kg */
        CA=A_SC*PI/180; /*進氣關閉時的角度轉弧度*/
        V=V_CA();/*進氣關閉時角度下的氣缸容積*/
        P=P_SC;
        T=T_SC;
        index=(int)A_SC;
, f* c. B" @* Q( `. l& `3 ~$ D        m_T[index]=T;
8 m* _3 ^) D3 @& |4 [& L        m_P[index]=P;
( Q4 m2 R( d& ~( x, V        m=P*V/(R*T);/*進氣關閉時角度下的氣缸內工質質量*/
5 S8 O* m# O; Q, ]3 N- }        mL=m-0.01*mBo;/*進氣關閉時角度下的氣缸內空氣質量*/
    mf=m-mL;
        Tpre=T;
) V, s3 P/ g# d* Q9 i9 J2 O        an=A_SC+1;
) {) N- H! P* s' ~: U" N) Z
        do{
- `. C) I/ U! b& h, H( S, C
! ~( Y% F( P; \( x6 Q  /*B 編寫歐拉法求解常微分方程程序 */

$ ]% g0 q: i/ i4 u7 [& i        }while(an!=A_SC+1);
}

void main(void)
{
0 v4 p6 ~* ]+ E8 z        int i;
0 A9 e1 D& f+ V! G8 C2 i3 t        FILE *fp;
/ ^$ n3 m" c0 `" t$ M        D=0.11;   /*發(fā)動機數(shù)據(jù)按各自的任務書輸入 */
        S=0.13;
        L=0.21;
& i2 a9 A" ]' O& [6 Y3 a        ge=236.0;
        n=2400;
2 l  v" n0 a, x! k/ u4 M        Ne=101568;
* C8 P$ P& d* u+ ~) R3 U        gnum=6;  ysb=17;   M=0.7;
9 i: `4 _7 x$ `8 R        Hu=41868000;
        R=287.08;
/ x, _$ a& ]9 y  T0 I        Pa=101000; /*大氣壓力*/
        Ta=300;    /*大氣溫度*/
        hs=(735.0+R)*Ta;  /*大氣焓值*/
        P_SC=90000;
        T_SC=313;
" G+ k' T( O4 j; v, N* Y  u9 k2 d, N
        A_SC=216; /*進氣門關閉*/
* j1 p9 |' ^2 ^. y& u* K        A_CB=350; /*燃燒開始*/
  W  \& ?- v! t3 }2 a        A_CE=410; /*燃燒結束*/
# l* f0 v* P# U2 q& r" O" _- {        A_EO=484; /*排氣門開啟*/
* Y+ W; W8 |: L! l        A_SO=705; /*進氣門開啟*/
( _: ]5 o& [, e2 h3 s: j( S# b  R, e        A_EC=15;  /*排氣門關閉*/
+ E" G" t/ H' Z' u
6 ^  t/ H% u4 r' \, ^) n9 X        Lo=14.4;

9 a0 _! q7 G. q* N( N    for(i=0;i<720;i++)
5 T" @( z* j1 j& R( ^. }        {
                m_P[i]=0;
' c9 S1 }3 t2 s" K7 ?. e7 F+ [' S                m_CA[i]=i;
/ t7 ^+ i9 N# d+ e4 v; c5 e9 l                m_T[i]=0;
        }

        Calculate();

6 d& E) v* W, ^% G. f  R    /*C編寫數(shù)據(jù)輸出程序，計算結束后按順序輸出0-720轉角的“缸壓，溫度”數(shù)據(jù)到outdata.txt文件 */

}
) b0 x4 F; y9 r4 i. w" {% W
+ ~1 a9 U! f- j  B7 T
8 `6 Q4 @' a% P, ?' @

---

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

---

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

---

歡迎光臨 機械社區(qū) (http://www.whclglass.com.cn/)

Powered by Discuz! X3.5

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

国产精品乱码一区-性开放网站-少妇又紧又爽视频-西西大胆午夜人体视频-国产极品一区-欧美成人tv-四虎av在线-国产无遮挡无码视频免费软件-中文字幕亚洲乱码熟女一区二区-日产精品一区二区三区在线观看-亚洲国产亚综合在线区-五月婷婷综合色-亚洲日本视频在线观看-97精品人人妻人人-久久久久久一区二区三区四区别墅-www.免费av-波多野结衣绝顶大高潮-日本在线a一区视频高清视频-强美女免费网站在线视频-亚洲永久免费

<style id="n0pgu"></style>

<style id="n0pgu"><i id="n0pgu"><strike id="n0pgu"></strike></i></style>

<style id="n0pgu"></style>

<style id="n0pgu"><optgroup id="n0pgu"></optgroup></style>