SHARE
TWEET

Untitled

a guest May 16th, 2019 88 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. clc  
  2. clear
  3.  
  4. a=80*10^-3 //alésage m
  5. c=84*10^-3 //course m
  6.  
  7. cylunit=((%pi*a^2)/4)*c //cylindrée unitaire
  8. cyl=cylunit*4
  9.  
  10. //cherchons V1 et V2 connaissant les deux formules de Pic et de la cyl
  11. V1=1782713*10^-9 //m^3
  12. V2=93827*10^-9 //m^3
  13.  
  14. Tc=V1/V2 //taux compression
  15.  
  16. //car Diesel
  17. V4=V1
  18.  
  19. //on a P1=Patm, T1 = temp ambiante (en K)
  20.  
  21. Pa=101325
  22. Y=1.4
  23. Cp=1000
  24. Cv=714
  25. P1=101325
  26. T1=15+273.15 //en K
  27.  
  28. //après calcul
  29.  
  30. P2=(P1*V1^Y)/(V2^Y) //Pascals, 62.5 bars
  31.  
  32. T2=(T1*V1^(Y-1))/V2^(Y-1)
  33.  
  34. P3=P2
  35.  
  36.  
  37. //trouvons T3
  38. r=286.8
  39. PCi=43*10^6
  40. rapport_stoechio=14.7
  41. m_air=(P1*V1)/(r*T1)
  42. m_fuel=m_air/rapport_stoechio
  43. Qin=m_fuel*PCi
  44. T3=Qin/((m_air+m_fuel)*Cp)+T2
  45.  
  46. V3=V2*T3/T2
  47.  
  48. P4=(P3*V3^Y)/(V4^Y)
  49.  
  50. T4=P4*T1/P1
  51.  
  52. Td=V4/V3 //taux détente
  53.  
  54.  
  55. rendement=1-(1/Y)*( 1/Td^Y-1/Tc^Y)/(1/Td-1/Tc)
  56.  
  57.  
  58.  
  59. disp ("Qin="+string(Qin)+" J")
  60.  
  61. //Verifions...
  62.  
  63. Q2_3=m_air*Cp*(T3-T2)
  64. Q4_1=m_air*Cv*(T1-T4)
  65. Qtot=Q2_3+Q4_1
  66.  
  67.  
  68. W2_3=-P2*(V3-V2)
  69. W1_2=(P2*V2-P1*V1)/(Y-1)
  70. W3_4=(P4*V4-P3*V3)/(Y-1)
  71.  
  72. Wtot=W2_3+W1_2+W3_4
  73.  
  74. Qtot
  75.  
  76. rpm=2300
  77.  
  78. //puissancce moyenne lors de la détendre =MEP
  79. //MEP=Mild effective pressure
  80.  
  81.  
  82. MEP=-Wtot/(V4-V3)
  83.  
  84.  
  85. P=MEP*cyl*rpm/120
  86.  
  87. xset('window',1)
  88.  
  89.  
  90. // de 1 à 2
  91.  
  92. cst1=P1*(V1^Y)
  93. volume12=[V2:0.00000001:V1]
  94. pression12=cst1./(volume12^Y)
  95. plot(volume12,pression12,'r')
  96.  
  97.  
  98. // de 2 à 3
  99.  
  100. volume23=[V2,V3]
  101. pression23=[P2,P3]
  102. plot(volume23,pression23,'r')
  103.  
  104. // de 3 à 4
  105. cst2=P3*(V3^Y)
  106. volume34=[V3:0.00000001:V4]
  107. pression34=cst2./(volume34^Y)
  108. plot(volume34,pression34,'r')
  109.  
  110. // de 4 à 1
  111. pression41=[P1,P4]
  112. volume41=[V1,V4]
  113. plot(volume41,pression41,'r')
  114.  
  115.  
  116.  
  117. disp ("Qtot="+string(Qtot)+" J")
  118. disp ("Wtot="+string(Wtot)+" J")
  119.  
  120.  
  121.  
  122. //graphique de la puissance en fonction de la température//////////////////////////////////////////////////
  123.  
  124.  
  125. xset('window',2)
  126. i=210
  127. while (i<320)
  128.    
  129.     i=i+0.25
  130.     T1=i
  131.     T2=(T1*V1^(Y-1))/V2^(Y-1)
  132.      m_air=(P1*V1)/(r*T1)
  133.     m_fuel=m_air/rapport_stoechio
  134.  
  135.     Qin=m_fuel*PCi
  136.     T3=(Qin/(m_air+m_fuel)*Cp)+T2
  137.     V3=V2*T3/T2
  138.     P3=P2
  139.     P4=(P3*V3^Y)/(V4^Y)
  140.     T4=P4*T1/P1
  141.  
  142.    
  143.     W3_4=(P4*V4-P3*V3)/(Y-1)
  144.     W1_2=(P2*V2-P1*V1)/(Y-1)
  145.     W2_3=-P2*(V3-V2)
  146.     Wtot=W2_3+W1_2+W3_4
  147.     MEP=-Wtot/(V4-V3)
  148.     P=MEP*cyl*rpm/120
  149.     plot (i,P,".r")
  150. end
  151.  
  152.  
  153. //remise a niveau////////////////////////////////////////////////////////////////////////////////////////
  154. clc  
  155. clear
  156.  
  157. a=80*10^-3 //alésage m
  158. c=84*10^-3 //course m
  159.  
  160. cylunit=((%pi*a^2)/4)*c //cylindrée unitaire
  161. cyl=cylunit*4
  162.  
  163. //cherchons V1 et V2 connaissant les deux formules de Pic et de la cyl
  164. V1=1782713*10^-9 //m^3
  165. V2=93827*10^-9 //m^3
  166.  
  167. Tc=V1/V2 //taux compression
  168.  
  169. //car Diesel
  170. V4=V1
  171.  
  172. //on a P1=Patm, T1 = temp ambiante (en K)
  173.  
  174. Pa=101325
  175. Y=1.4
  176. Cp=1000
  177. Cv=714
  178. P1=101325
  179. T1=15+273.15 //en K
  180.  
  181. //après calcul
  182.  
  183. P2=(P1*V1^Y)/(V2^Y) //Pascals, 62.5 bars
  184.  
  185. T2=(T1*V1^(Y-1))/V2^(Y-1)
  186.  
  187. P3=P2
  188.  
  189.  
  190. //trouvons T3
  191. r=286.8
  192. PCi=43*10^6
  193. rapport_stoechio=14.7
  194. m_air=(P1*V1)/(r*T1)
  195. m_fuel=m_air/rapport_stoechio
  196. Qin=m_fuel*PCi
  197. T3=Qin/((m_air+m_fuel)*Cp)+T2
  198.  
  199. V3=V2*T3/T2
  200.  
  201. P4=(P3*V3^Y)/(V4^Y)
  202.  
  203. T4=P4*T1/P1
  204.  
  205. Td=V4/V3 //taux détente
  206.  
  207.  
  208. rendement=1-(1/Y)*( 1/Td^Y-1/Tc^Y)/(1/Td-1/Tc)
  209.  
  210.  
  211.  
  212. disp ("Qin="+string(Qin)+" J")
  213.  
  214. //Verifions...
  215.  
  216. Q2_3=m_air*Cp*(T3-T2)
  217. Q4_1=m_air*Cv*(T1-T4)
  218. Qtot=Q2_3+Q4_1
  219.  
  220.  
  221. W2_3=-P2*(V3-V2)
  222. W1_2=(P2*V2-P1*V1)/(Y-1)
  223. W3_4=(P4*V4-P3*V3)/(Y-1)
  224.  
  225. Wtot=W2_3+W1_2+W3_4
  226.  
  227. Qtot
  228.  
  229. rpm=2300
  230.  
  231. //puissancce moyenne lors de la détendre =MEP
  232. //MEP=Mild effective pressure
  233.  
  234.  
  235. MEP=-Wtot/(V4-V3)
  236.  
  237.  
  238. P=MEP*cyl*rpm/120
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247. //graphique de la puissance en fonction de la pression ///////////////////////////////////////////////////////////////////////////////////////
  248.  
  249.  
  250. xset('window',3)
  251.  
  252. i=50000
  253.  
  254. while i<(152000)
  255.    
  256.    
  257.     i=i+400
  258.     P1=i
  259.     P2=(P1*V1^Y)/(V2^Y)
  260.    
  261.     T2=(T1*V1^(Y-1))/V2^(Y-1)
  262.      m_air=(P1*V1)/(r*T1)
  263.     m_fuel=m_air/rapport_stoechio
  264.    
  265.      Qin=m_fuel*PCi
  266.     T3=(Qin/(m_air+m_fuel)*Cp)+T2
  267.     V3=V2*T3/T2
  268.     P3=P2
  269.     P4=(P3*V3^Y)/(V4^Y)
  270.     T4=P4*T1/P1
  271.    
  272.     W3_4=(P4*V4-P3*V3)/(Y-1)
  273.     W1_2=(P2*V2-P1*V1)/(Y-1)
  274.     W2_3=-P2*(V3-V2)
  275.     Wtot=W2_3+W1_2+W3_4
  276.     MEP=-Wtot/(V4-V3)
  277.     P=MEP*cyl*rpm/120
  278.     plot (i,P,".r")
  279. end
  280.  
  281. //remise a niveau ///////////////////////////////////////////////////////////////////////////////////////
  282. clc  
  283. clear
  284.  
  285. a=80*10^-3 //alésage m
  286. c=84*10^-3 //course m
  287.  
  288. cylunit=((%pi*a^2)/4)*c //cylindrée unitaire
  289. cyl=cylunit*4
  290.  
  291. //cherchons V1 et V2 connaissant les deux formules de Pic et de la cyl
  292. V1=1782713*10^-9 //m^3
  293. V2=93827*10^-9 //m^3
  294.  
  295. Tc=V1/V2 //taux compression
  296.  
  297. //car Diesel
  298. V4=V1
  299.  
  300. //on a P1=Patm, T1 = temp ambiante (en K)
  301.  
  302. Pa=101325
  303. Y=1.4
  304. Cp=1000
  305. Cv=714
  306. P1=101325
  307. T1=15+273.15 //en K
  308.  
  309. //après calcul
  310.  
  311. P2=(P1*V1^Y)/(V2^Y) //Pascals, 62.5 bars
  312.  
  313. T2=(T1*V1^(Y-1))/V2^(Y-1)
  314.  
  315. P3=P2
  316.  
  317.  
  318. //trouvons T3
  319. r=286.8
  320. PCi=43*10^6
  321. rapport_stoechio=14.7
  322. m_air=(P1*V1)/(r*T1)
  323. m_fuel=m_air/rapport_stoechio
  324. Qin=m_fuel*PCi
  325. T3=Qin/((m_air+m_fuel)*Cp)+T2
  326.  
  327. V3=V2*T3/T2
  328.  
  329. P4=(P3*V3^Y)/(V4^Y)
  330.  
  331. T4=P4*T1/P1
  332.  
  333. Td=V4/V3 //taux détente
  334.  
  335.  
  336. rendement=1-(1/Y)*( 1/Td^Y-1/Tc^Y)/(1/Td-1/Tc)
  337.  
  338.  
  339.  
  340. disp ("Qin="+string(Qin)+" J")
  341.  
  342. //Verifions...
  343.  
  344. Q2_3=m_air*Cp*(T3-T2)
  345. Q4_1=m_air*Cv*(T1-T4)
  346. Qtot=Q2_3+Q4_1
  347.  
  348.  
  349. W2_3=-P2*(V3-V2)
  350. W1_2=(P2*V2-P1*V1)/(Y-1)
  351. W3_4=(P4*V4-P3*V3)/(Y-1)
  352.  
  353. Wtot=W2_3+W1_2+W3_4
  354.  
  355. Qtot
  356.  
  357. rpm=2300
  358.  
  359. //puissancce moyenne lors de la détendre =MEP
  360. //MEP=Mild effective pressure
  361.  
  362.  
  363. MEP=-Wtot/(V4-V3)
  364.  
  365.  
  366. P=MEP*cyl*rpm/120
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top