TP Cycle Diesel

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 8

Nom et prénom:

………………………………………………………………………………………………………………………………………………………………
…………………………

/20 Etude du cycle diesel par le logiciel Matlab

I. Objectif du TP

Le but de ce TP est l’étude thermodynamique du cycle diesel à travers l’étude de l’effet de


la variation du rapport de compression volumétrique « r » sur les performances du moteur
essence. Le logiciel EES (Engineering Equation Solver) est utilisé pour cette étude.

Figure 1.

Données du problème

T1=27°C, P1=1 bar, Qin=1800kJ/kg, Taux de compression: r= V1/V2, Taux d’injection,


a=V3/V2, γ=1.4, Cv=0.718 kJ/kg.K, Cp=1.005 kJ/kg.K et R=0.287 kJ/kg.K.

1
II. Travail demandé

Q.1 Ecrire le programme du cycle OTTO (Code MATLAB


clear all
close all
clc

%inputs
D = 0.20; %cylinder diameter or bore of cylinder in [m]
L = 0.30; %stroke length in [m]
conn_rod = 0.42; %Connecting rod length in [m]
Cr = 15; %Compression ratio
cutoff = 8/100;
%Calculate cut off ratio
%COr = cut off ratio
%cut_off = (COr-1)/(Cr-1)
%COr-1 = (Cr-1)*Cutoff
COr = (Cr-1)*cutoff+1;
gamma = 1.4; %cp/cv constant

%State variables
p1 = 1*10^5; %pressure at point 1 in [pascal]
t1 = 27+273; %temp at point 1 in [kelvin]

%Calculate Swept Volume & Clearance Volume


Vs = pi/4*D^2*L; %swept volume
Vc = Vs/(Cr-1); %clearance volume

%Isentropic compression process


%State variables at point 1-2
v1 = Vs+Vc;

constant = p1*v1^gamma;
V_comp = prac(D,L,Cr,conn_rod,180,0);
P_comp = constant./V_comp.^gamma;

%Constant pressure heat addition process


%state variables at point 2-3
v2 = Vc;
%p2*v2^gamma = p1*v1^gamma
%Cr = v1/v2
p2 = p1*(Cr)^gamma;
%t2/t1 = (v2/v2)^(gamma-1)
t2 = t1*(v2/v1)^(gamma-1);

%state variables at point 3-4


%COr = v3/v2
v3 = COr*v2;
p3 = p2; %pressure at point 2 & 3 are remains same
%t3/t2 = v3/v2
t3 = t2*(COr);

constant = p3*v3^gamma;
V_exp = prac(D,L,Cr,conn_rod,0,180);
for i = 1:length(V_exp)
if V_exp(i)<v3
V_exp(i)=v3;
end

2
end
P_exp = constant./V_exp.^gamma;

%state variables at point 4-1


v4 = v1; %volume at point 1 & 4 are remains same
%p3*v3^gamma = p4*v4^gamma
p4 = p3*(v3/v4)^gamma;
%t4/t3 = (v4/v3)^(gamma-1)
t4 = t3*(v4/v3)^(gamma-1);

%efficiency of diesel cycle [Nth]


%Nth = 1-(1/gamma*(Cr)^(gamma-1))*[COr^gamma-1/COr-1]
term1 = 1/(gamma*(Cr)^(gamma-1));
term2 = (COr^gamma-1)/(COr-1);
N = 1-term1*term2;
Nth = N*100;

%plotting
hold on
plot(v1,p1,'*','color','k')
plot(V_comp,P_comp,'color','k')
text(v1,p1,'State 1','color','g')
plot(v2,p2,'*','color','k')
text(v2,p2,'State 2','color','g')
plot(v3,p3,'*','color','k')
plot(V_exp,P_exp,'color','k')
text(v3,p3,'State 3','color','g')
plot(v4,p4,'*','color','k')
text(v4,p4,'State 4','color','g')
line([v2 v3],[p2 p3],'color','k')
line([v4 v1],[p4 p1],'color','k')
title('PV Diagram "Diesel Cycle"')
xlabel('Volume')
ylabel('pressure')

%clear all
%close all
%clc
function V = prac(D,L,Cr,conn_rod,start_crank,end_crank)
%D = 0.20; %bore
%L = 0.30; %stroke
%Cr = 15; %compression ratio
%conn_rod = 0.42; %connecting rod

a = L/2;
R = conn_rod/a;

theta = linspace(start_crank,end_crank,100);

%calculate swept volume and clearance volume


Vs = (pi/4)*D^2*L;
Vc = Vs/(Cr-1);

%calclute formula
%V/Vc = 1+0.5*(Cr-1)[(R+1-cosd(theta)) - (R^2-sin(theta)^2)^0.5]
term1 = 0.5*(Cr-1);
term2 = R+1 - cosd(theta);
term3 = (R^2 - sind(theta).^2).^0.5;

V = (1+term1*(term2-term3)).*Vc;

3
end

Q.2 Remplir les tableaux

r P1 P2 P3 P4 V1 V2 V3 V4
(bar) (bar) (bar) (bar) (m3/kg) (m3/kg) (m3/kg) (m3/kg)
4
6
8
10
12
Tableau 1. Points du cycle P-V

r T1 T2 T3 T4
(K) (K) (K) (K)
4
6
8
10
12

Tableau 2. Points du cycle T-S

r T3 V2 Qin Qout Wnet ηth


(K) ( m3/kg ) (kJ/kg) (kJ/kg) (kJ/kg)

4
6
8
10
12

Tableau 3. Performances du cycle diesel

Avec: r : rapport de compression volumétrique et T1, T2, T3, T4, P1, P2, P3, P4, V1, V2, V3, V4
sont les points du cycle.

4
Q.3 Tracer la variation du rendement du cycle « ηth » en fonction de « r »

Q.4 Tracer la variation du travail net du cycle « Wnet » en fonction de « r »

Q.5 Tracer la variation l’énergie thermique perdue « Qout » en fonction de « r »

5
Q.6 Interpréter les résultats

……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..

6
………………………………………………………………………………………………
……………………………………………………………………………………………….
Q.7 Conclusion générale

……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………

……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………

……………………………………………………………………………………………..
………………………………………………………………………………………………
……………………………………………………………………………………………….
………………………………………………………………………………………………
………………………………………………………………………………………………
……………………………………………………………………………………………..

7
………………………………………………………………………………………………
……………………………………………………………………………………………….

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy