Scribd
Upload
20 views2 pages

Code Matlab

cc

Uploaded by

hoanglinhvippro1234
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
20 views2 pages

Code Matlab

cc

Uploaded by

hoanglinhvippro1234
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

R = 0.

488; %dm
D = 0.805; %dm
S = 0.976; %dm
l = 1.43; %dm
Tsn = 13;
Vd = 0.496; %lit
Sp = (pi*D^2)/4; % dien tich piston
lamda = R/l;
Pa = 0.09; %MN/m
Pr = 0.12;
Pz = 9.499;
Pc = 3.061;
ro = 1;
Vc = Vd/(Tsn -1); %dm^3
Vz = ro*Vc;
n1 = 1.37;
n2 = 1.295;
n = 6600; %vong/phut
w = (pi*n)/30; %rad/s
Xa = R*(1-cosd(180) + (lamda/4).*(1-cosd(2.*180)));
Va =Xa*Sp + Vc; %(lit)
%Góc xupap nap mo som: 1 do Truoc DCT
%Góc xupap nap dong muon: 65 do Truoc DCD
%Góc xupap xa mo som: do Sau DCD
%Góc xupap xa dong muon: 43 do Sau DCT
%Góc danh lua som: 20 do
%HIEU CHI THAI – NAP
ahc1 = [0 20 43];
phc1 = [Pr 0.1 Pa];
a1 = linspace(0,43,1000);
X1 = R*(1-cosd(a1) + (lamda/4).*(1-cosd(2.*a1)));
V1 = X1*Sp + Vc;
Pl = interp1(ahc1,phc1,a1, 'spline');
%QUA TRINH NAP
a2 = linspace(43,180,1000);
X2 = R*(1-cosd(a2)+(lamda/4)*(1-cosd(2*a2)));
V2 = X2*Sp+Vc;
P2 = linspace(Pa,Pa,1000);%MN/m2
% QUA TRINH NEN
a3 = linspace(180,340,1000);
X3 = R*(1-cosd(a3)+(lamda/4)*(1-cosd(2*a3)));
V3 = X3*Sp+Vc;
P3 = Pa*((Va./V3).^n1);%MN/m2
% HIEU CHINH NEN - CHAY
Pcc =(Pz-Pc)/3 + Pc;
ahc4 = [340 350 360];
phc4 = [max(P3) 2.35 Pcc];
a4 = linspace(340,360,1000);
X4 = R*(1-cosd(a4)+(lamda/4)*(1-cosd(2*a4)));
V4=X4*Sp+Vc;
P4 = interp1(ahc4,phc4,a4, 'spline');
% Diem z”
X41 = R*(1-cosd(370) + (lamda/4).*(1-cosd(2.*370)));
V41 = X41*Sp + Vc;
% QUA TRINH GIAN NO
a6 = linspace(370,497,1000);
X6 = R*(1-cosd(a6) + (lamda/4).*(1-cosd(2.*a6)));
V6 = X6.*Sp + Vc;
P6 = Pz.*(Vz./V6).^n2;
% HIEU CHINH CHAY - GIAN NO
ahc5=[360 365 370];
phc5=[Pcc 5.91 max(P6)];
a5=linspace(360,370,1000);
X5 = R*(1-cosd(a5) + (lamda/4).*(1-cosd(2.*a5)));
V5 = X5.*Sp + Vc;
P5 = interp1(ahc5,phc5,a5,'spline');
% HIEU CHINH GIAN NO - THAI
Pb = 0.5;
% Diem b’
Pb = min(P6);
% Diem b”
Vb = Va;
Pbb =(Pb-Pr)/2+Pr;
ahc7 = [497 540 580];
phc7 = [Pb Pbb Pr];
a7 = linspace(497,580,1000);
X7 = R*(1-cosd(a7) + (lamda/4).*(1-cosd(2.*a7)));
V7 = X7.*Sp + Vc;
P7 = interp1(ahc7,phc7,a7,'spline');
% QUA TRINH THAI
a8 = linspace(580,720,1000);
X8 = R*(1-cosd(a8) + (lamda/4).*(1-cosd(2.*a8)));
V8 = X8*Sp + Vc;
P8 = linspace(Pr,Pr,1000);
V = [V1 V2 V3 V4 V5 V6 V7 V8];
P = [Pl P2 P3 P4 P5 P6 P7 P8];
a0 = [a1 a2 a3 a4];
ap= [a1 a2 a3 a4 a5 a6 a7 a8];
X = [X1 X2 X3 X4 X5 X6 X7 X8];
% DO THI P-V
figure(1)
plot(V,P);
title('DO THI P-V');
xlabel('The tich V (dm3)');
ylabel('Ap suat P (MN/dm2)');
grid on
hold on

You might also like