%TINH TOAN DONG CO DOT TRONG

%-----------------------------------------------------------------------
S = 1.049; % hanh trinh piston
D = 0.954; % duong kinh piston
B = 0.954; % duong kinh xi lanh
R =S/2; % ban kinh truc khuyu
lamda =1/3.5; % thong so ket cau
Sp =(pi*(D^2))/4; % dien tich thiet dien mat cat vuong goc truc xi lanh
vd = 1.141; % the tich cong tac
vc =vd/(16.3-1); % the tich buong chay
va =(vd + vc); % dung tich xy lanh
n1 =1.37; % chi so nen da bien trung binh
n2 =1.27; % chi so gian no da bien trung binh
% phía trên tính theo đơn vị (dm)
w = (2*pi*3600)/60;
vz = 1.3*vc; %(do dong co diesel)
vb = va; pa = 0.1056; %MN/m2
pc = 4.835; %MN/m2
pz = 9.669; %MN/m2
pb = 0.421; %MN/m2
pr = 0.109; %MN/m2
%QUA TRINH NAP (doan rr'')
ahc1=[0 7 14]; %xuppap thải đóng muộn góc 14 độ
phc1=[pr (pr+pa)./2 pa];
a1=linspace(0,14,141);
x1=R.*(1-cosd(a1)+(lamda/4).*(1-cosd(2*a1)));
v1=(x1.*Sp)+vc;
p1=interp1(ahc1,phc1,a1,'pchip');
12
j1=(R/10)*(w^2).*(cosd(a1)+lamda.*cosd(2.*a1));
%QUA TRINH NEN (doan ac')
a3=linspace(180,348, 1681); % góc phun nhiên liệu sớm 12 độ
x3=R.*(1-cosd(a3)+(lamda/4).*(1-cosd(2*a3)));
v3=(x3.*Sp)+vc;
p3= pa.*(va./v3).^n1;
j3=(R/10)*((w)^2).*(cosd(a3)+lamda.*cosd(2.*a3));
%QUA TRINH NAP (doan r''a)
a2=linspace(14,180,1661);
x2=R.*(1-cosd(a2)+(lamda/4).*(1-cosd(2*a2)));
v2=(x2.*Sp)+vc;
p2=linspace(pa,min(p3),1661); % Nhằm cho đường p2 và p3 giao nhau tại 180 độ
j2=(R/10)*(w^2).*(cosd(a2)+lamda.*cosd(2.*a2));
%QUA TRINH CHAY VA GIAN NO (doan c'z")
%%%VE DOAN (c'c")
%%%%xac dinh diem c'
vc1 = min (v3); pc1 = max (p3);
%%%%xac dinh toan do diem c"
pcz1 = pz - pc;
pc2 = pcz1/3 + pc;
vc2 = vc; ahc4=[348 354 360];
phc4=[pc1 (pc1+pc2)./2 pc2];
a4=linspace(348,360, 121);
x4=R.*(1-cosd(a4)+(lamda/4).*(1-cosd(2*a4)));
v4=(x4.*Sp)+vc;
p4=interp1(ahc4,phc4,a4,'pchip');
j4=(R/10)*((w)^2).*(cosd(a4)+lamda.*cosd(2.*a4));
%%%VE DOAN AP SUAT CUC DAI
13
a22=linspace(380,492, 1121);
xz2=R.*(1-cosd(a22)+(lamda/4).*(1-cosd(2*a22)));
vz2=(xz2.*Sp)+vc;
pz2=pz.*(vz./vz2).^n2;
pz22=max(pz2);
ahc5=[360 370 380];
phc5=[pc2 pz pz22];
a55=linspace(360,380,201);
x55=R.*(1-cosd(a55)+(lamda/4).*(1-cosd(2*a55)));
v55=(x55.*Sp)+vc;
p55=interp1(ahc5,phc5,a55,'pchip');
j55=(R/10)*((w)^2).*(cosd(a55)+lamda.*cosd(2.*a55));
%QUA TRINH CHAY VA GIAN NO
a7=linspace(380,492,1121); % góc mở sớm xuppap thải 48 độ
x7=R.*(1-cosd(a7)+(lamda/4).*(1-cosd(2*a7)));
v7=(x7.*Sp)+vc;
p7=pz.*(vz./v7).^n2;
j7=(R/10)*((w)^2).*(cosd(a7)+lamda.*cosd(2.*a7));
%QUA TRINH CHAY VA GIAN NO (doan b'->b")
%%%xac dinh diem b'
a8=linspace(492,540, 481);
x8=R.*(1-cosd(a8)+(lamda/4).*(1-cosd(2.*a8)));
v8=(x8.*Sp)+vc; pB1=pz.*( vz./v7).^n2;
pb1=min(pB1);
vb1=max(v7); %%%xac dinh diem b" pb2=((pb+pa)/2);
vb2=va;
%%%
xb3=R.*(1-cosd(ab3)+(lamda/4).*(1-cosd(2*ab3)));
14
vb3=(xb3.*Sp)+vc;
pb3=pz.*(vz./vb3).^n2;
ahc2=[492 516 540];
phc2=[pb1 pb3 pb2];
p8=interp1(ahc2,phc2,a8,'pchip');
j8=(R/10)*((w)^2).*(cosd(a8)+lamda.*cosd(2.*a8));
%QUA TRINH THAI 1 (b''->r')
a9=linspace(540,580, 401); % góc mở sớm xuppap nạp 12 độ
x9=R.*(1-cosd(a9)+(lamda/4).*(1-cosd(2.*a9)));
v9=(x9.*Sp)+vc;
ahc3=[540 560 580];
phc3=[pb2 (pb2+pr)./2 pr];
p9=interp1(ahc3,phc3,a9,'spline');
j9=(R/10)*((w)^2).*(cosd(a9)+lamda.*cosd(2.*a9));
%QUA TRINH THAI 2 (r'->r)
a10=linspace(580,720,1401);
x10=R.*(1-cosd(a10)+(lamda/4).*(1-cosd(2.*a10)));
v10=(x10.*Sp)+vc;
p10=linspace(pr,pr,1401);
j10=(R/10)*((w)^2).*(cosd(a10)+lamda.*cosd(2.*a10));
%CAC LENH VE VA CHINH SUA DO THI
a =[a1 a2 a3 a4 a55 a7 a8 a9 a10];
x =[x1 x2 x3 x4 x55 x7 x8 x9 x10];
v =[v1 v2 v3 v4 v55 v7 v8 v9 v10];
p =[p1 p2 p3 p4 p55 p7 p8 p9 p10];
j =[j1 j2 j3 j4 j55 j7 j8 j9 j10];
figure %DO THI P-V
plot(v,p,'r','linewidth',1.5)
15
hold on;
title('DO THI P-V');
xlabel('The tich(lit)');
ylabel('ap suat P(MN/m2)');
grid on
figure %DO THI CHUYEN VI PISTON THEO GOC
plot(a,x,'r','linewidth',1.5)
hold on;
title('DO THI CHUYEN VI');
xlabel('Goc quay truc khuyu (do)');
ylabel('Chuyen vi x (dm)');
grid on
figure
%DO THI VAN TOC
vt=(R)*(w).*(sind(a)+lamda/2.*sind(2.*a));
plot(a,vt,'r','linewidth',1.5)
hold on;
title('DO THI VAN TOC PISTON');
xlabel('Goc quay truc khuyu (do)');
ylabel('Van toc v (dm/s) ');
grid on
figure %DO THI GIA TOC
plot(a,j,'r','linewidth',1.5)
hold on;
title('DO THI GIA TOC PISTON');
xlabel('Goc quay truc khuyu (do)');
ylabel('Gia toc J (dm/s2) ');
grid on
16
figure %DO THI Pkt Pj P1
pkt=(p-0.1013);
plot(a,pkt,'b','linewidth',1.5)
hold on
grid on
pj=-14.75*R*(w.^2)*(10.^(-6))*(cosd(a)+lamda*(cosd(2*a)));
plot(a,pj,'g','linewidth',1.5)
hold on
grid on
p1=pkt+pj;
plot(a,p1,'r','linewidth',1.5)
title('DO THI Pkt Pj P1');
xlabel('Goc quay truc khuyu(do)');
ylabel('Pkt(MN/m2) Pj(MN/m2) P1(MN/m2)');
legend('Pkt','Pj','P1')
%lay so lieu P V Pkt Pj P1 x vt j
xlswrite('BsolieuPVPktPjP1.xlsx',[a',pkt',pj',p1']);
xlswrite('solieuxvj.xlsx',[a',x',vt',j',v',p']);

17