(TIỂU LUẬN) báo cáo CUỐI kỳ TÍNH TOÁN ĐỘNG cơ đốt TRONG

TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA CƠ KHÍ ĐỘNG LỰC - - BÁO CÁO CUỐI KỲ TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG SVTH : BÙI CÔNG THƯ MSSV: 18145465 SVTH : MSSV: SVTH : MSSV: SVTH : MSSV: GVHD: PHAN TRẦN KHẢI 18145370 NGUYỄN THÀNH ĐẠT 18145334 NGUYỄN HỮU LƯƠNG 18145392 TS NGUYỄN VĂN TRẠNG Tp Hồ Chí Minh, tháng 12 năm 2020 TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA CƠ KHÍ ĐỘNG LỰC - - BÁO CÁO CUỐI KỲ Chuyên ngành: Cơng nghệ Kỹ thuật Ơ tơ Tên đề tài TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG SVTH : BÙI CƠNG THƯ MSSV: 18145465 SVTH : MSSV: SVTH : MSSV: SVTH : MSSV: GVHD: PHAN TRẦN KHẢI 18145370 NGUYỄN THÀNH ĐẠT 18145 NGUYỄN HỮU LƯƠNG 18145392 TS NGUYỄN VĂN TRẠNG Tp Hồ Chí Minh, tháng 12 năm 2020 TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP HCM CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Bộ mơn Tính Tốn Động Cơ Đốt Trong PHIẾU NHẬN XÉT BÁO CÁO CUỐI KỲ (Dành cho giảng viên hướng dẫn) Họ tên sinh viên: Bùi Công Thư Họ tên sinh viên: Phan Trần Khải Họ tên sinh viên: Nguyễn Thành Đạt Họ tên sinh viên: Nguyễn Hữu Lương Tên đề tài:Tíí́nh tốn động đốt Ngành đào tạo: Cơng nghệ Kỹ thuật Ơ tơ Họ tên GV hướng dẫn: TS Nguyễn Văn Trạng Ý KIẾN NHẬN XÉT Nhận xét tinh thần, thái độ làm việc sinh viên (không đánh máy) Nhận xét kết thực (không đánh máy) 2.1.Kết cấu, cách thức trình bày: 2.2 Nội dung: (Cơ sở lý luận, tíí́nh thực tiễn khả ứng dụng đồ án, hướng nghiên cứu tiếp tục phát triển) i 2.3.Kết đạt được: 2.4 Những tồn (nếu có): Đánh giá: TT Mục đánh giá Hình thức kết cấu Đúng format vơi đu ca hinh thưc va nôịdung cua cac muc Mục tiêu, nhiệm vụ, tổng quan đề tài Tíí́nh cấp thiết đề tài Nội dung Khả ứng dụng kiến thức toán học, khoa học kỹ thuật, khoa hoc xã hơị… Khả thực hiện/phân tíí́ch/tổng hợp/đánh giá Khả thiết kế chế tạo hệ thống, thành phần, quy trình đáp ứng yêu cầu đưa với ràng buộc thực tế Khả cải tiến phát triển Khả sử dụng công cụ kỹ thuật, phần mềm chuyên ngành… Đánh giá khả ứng dụng đề tài Sản phẩm cụ thể Tổng điểm ii Kết luận:  Được phép bảo vệ  Không phép bảo vệ TP.HCM, ngày tháng 12 năm 2020 Giảng viên hướng dẫn ((Ký, ghi rõ họ tên) iii TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP HCM CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Bộ mơn Tính Tốn Động Cơ Đốt Trong PHIẾU NHẬN XÉT BÁO CÁO CUỐI KỲ (Dành cho giảng viên phản biện) Họ tên sinh viên: Bùi Công Thư Họ tên sinh viên: Phan Trần Khải Họ tên sinh viên: Nguyễn Thành Đạt Họ tên sinh viên: Nguyễn Hữu Lương Tên đề tài:Tíí́nh tốn động đốt Ngành đào tạo: Cơng nghệ Kỹ thuật Ơ tơ Họ tên GV phản biện: (Mãã̃ GV) Ý KIẾN NHẬN XÉT Kết cấu, cách thức trình bày: Nội dung: (Cơ sở lý luận, tíí́nh thực tiễn khả ứng dụng đồ án, hướng nghiên cứu tiếp tục phát triển) Kết đạt được: iv Những thiếu sót tồn tại: Câu hỏi: Đánh giá: v TT Mục đánh giá Hình thức kết cấu Đúng format vơi đu ca hinh thưc va nôịdung cua cac muc Mục tiêu, nhiệm vụ, tổng quan đề tài Tíí́nh cấp thiết đề tài Nội dung Khả ứng dụng kiến thức toán học, khoa học kỹ thuật, khoa hoc xã hôị… Khả thực hiện/phân tíí́ch/tổng hợp/đánh giá Khả thiết kế, chế tạo hệ thống, thành phần, quy trình đáp ứng yêu cầu đưa với ràng buộc thực tế Khả cải tiến phát triển Khả sử dụng công cụ kỹ thuật, phần mềm chuyên ngành… Đánh giá khả ứng dụng đề tài Sản phẩm cụ thể Tổng điểm Kết luận:  Được phép bảo vệ  Không phép bảo vệ TP.HCM, ngày tháng 12 năm 2020 Giảng viên phản biện ((Ký, ghi rõ họ tên) vi 17 Pm = A + B*Vp + (Pr - Pa); % Xác ??nh áp su?t có ích trung bình Pe Pe = Pi - Pm; % Hi?u su?t c? gi?i Nm = - Pm/Pi; % Xác ??nh hi?u su?t Ni = 8.314*(M1*Pi*Tk)/(QH*Pk*Nv); % Hi?u su?t có ích Ne = 8.314*(M1*Pe*Tk)/(QH*Nv*Pk); % Tiêu hao nhiên li?u ch? th? Gi = 3600/(QH*Ni)*10^3; % Tiêu hao nhiên li?u có ích Ge = 3600/(QH*Ne)*10^3; % TINH THONG SO KET CAU CUA ??NG C? % lamdaphay lamdaphay = 1/3.5; % Th? tích cơng tác Vh Vh = (30*4*120)/(Pe*n*8); Vh = double(Vh); % ???ng kính piston D = nthroot((4*Vh)/(pi*StrenD),3); D = double(D); % Hành trình piston S = (StrenD)*D; %%%%%%%%%%%%%%% Vc = Vh/(ep - 1); Va = Vh + Vc; Vz = Vc*p; Vb=Va; %CAC DIEM DAC BIET %-diem a y0 = Pa; x0 = Va; %-diem c y1= Pc; x1 = Vc; %-diem z y2 = Pz; x2 = p*Vc; %-diem b y3 = Pb; x3 = Va; %-diem r y4 = Pr; x4 = Vc; % Ve tung diem dac biet x = [x0 x1 x2 x3 x4]; y = [y0 y1 y2 y3 y4]; plot(x,y,'b.'); hold on; %%%%%%%%%%% % DIEM XAC DINH TREN DO THI plot(Vc,Pz,'b.')%Z' 18 plot((Vc+x2)/2,Pz,'b.')%Z'' plot(x1,(Pz-y1)/3+y1,'b.')%C'' plot(x0,y0+(y3-y0)/2,'b.')%B'' lc= Pa*(Va.^n1./0.1126.^n1); plot(0.1126,lc,'b.')%C' lb=(Pz)*(Vz.^n2./1.1583.^n2); plot(1.1583,lb,'b.')%B' %%%%%%%%%%%%%%%%% plot(1.2173,y0,'b.')%A' plot(0.1451,y0,'b.')%R' plot(0.1686,y4,'b.')%R'' %%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%% % DUNG DUONG CONG NEN phu1 = (Va - 0.09012)/0.001; Vxn = linspace(Va,0.09012,phu1); Pxn = Pa*(Va.^n1./Vxn.^n1); plot(Vxn,Pxn,'r-'); % DUNG DUONG DAN NO phu2= (1.1583 -0.1401)/0.001; Vxg = linspace(0.1401,1.1583,phu2); Pxg = (Pz)*(Vz.^n2./Vxg.^n2); plot(Vxg,Pxg,'r-'); ylim([-1 7]) xlim([0 1.5]) axis fill; %%%%%%%%%%%%%%%% %%% doan x= 0.08112:0.001:0.09013; y= ( 3.42e-09)*x.^-8.13+2.345; plot (x,y,'r-'); plot ([x1, x1],[ (Pz-y1)/3+y1 , 5.958],'r-'); x= 0.08612:0.001:0.1403; y= 408.2*x.*exp(-9.958*x) -8.364; plot (x,y,'r-'); x= 0.08112:0.001:0.08612; y= (-1.716e-15)*x.^-13.56+7.022 plot (x,y,'r-'); %%% doan plot ([x4,1.15 ],[ y4 , y4]); y= 0.105:0.001:0.4649; x= 8.804*y.*exp(-4.097*y) +0.5488; plot (x,y,'r-'); %%% doan x= 0.08112:0.001:0.1451; y= 0.04215*x.^-0.3634; plot (x,y,'r-'); plot ([0.1451,x0 ],[ y0 , y0],'r-'); % VE DUONG NUA DUONG TRON BEN DUOI Ophay = (lamdaphay*(Vh/2))/2; xCenter_2 = (Va+Vc)/2; yCenter_2 = 0; radius = (Va-Vc)/2 ; plot(xCenter_2,0,'.'); plot(xCenter_2 + Ophay,0,'r.'); fxy =@(x,y) (x - xCenter_2).^2 + (y - yCenter_2).^2 - radius.^2 +(y>=0) ; 19 fimplicit(fxy); %%%%%%%%%%%% % DIEM XAC DINH TREN DUONG TRON plot(0.1126,0.1913,'b.'); % la thetaS=15 plot(0.1686,-0.3209,'b.'); % la theta1=26 plot(0.1451,-0.2753,'b.'); % la theta4=22 plot(1.2173,-0.3718,'b.'); % la theta2=48 plot(1.1583,-0.4406,'b.'); % la theta3=56 hold off; %%%%%%%%%%%% title('Do thi cong P - V') ylabel('P'); xlabel('V'); close all; clc; %TINH TOAN DONG LUC HOC CO CAU TRUC KHUYU THANH TRUYEN %Dong Hoc Cua Piston %Chuyen vi cua piston R= S/2*10^-1; goc_alpha = : pi/30 : 2*pi; goc_beta= asin(lamdaphay*sin(goc_alpha)); Sp = R*((1- cos(goc_alpha))+ (lamdaphay/4)*(1 - cos(2*goc_alpha))); subplot ( 3,1,1); plot(goc_alpha*180/pi,Sp,'r-'); title('Chuyen vi cua piston') legend('Sp','Sp1','Sp2') ; ylabel('Sp'); xlabel('alpha'); hold on omega = pi*n/30; Vp = R*omega*(sin(goc_alpha)+ lamdaphay/2*sin(2*goc_alpha)); Vp1 = R*omega*sin(goc_alpha); Vp2 = R*omega*(lamdaphay/2)*sin(2*goc_alpha); subplot(3,1,2); plot(goc_alpha*180/pi, Vp, 'r-'); hold on plot(goc_alpha*180/pi, Vp1, 'b-'); plot(goc_alpha*180/pi, Vp2, 'g-'); title('Van toc cua piston') legend('Vp','Vp1','Vp2') ; ylabel('Vp'); xlabel('alpha'); Jp = R*omega^2*(cos(goc_alpha) + lamdaphay*cos(2*goc_alpha)); Jp1 = R*omega^2*cos(goc_alpha); Jp2 = Jp - Jp1; subplot(3,1,3); plot(goc_alpha*180/pi, Jp , 'r-'); hold on plot(goc_alpha*180/pi, Jp1 , 'b-'); plot(goc_alpha*180/pi, Jp2 , 'g-'); title('Gia toc cua piston') legend('Jp','Jp1','Jp2') ; ylabel('Jp'); 20 xlabel('alpha'); maphay = 0.35*mttphay; mbphay = 0.65*mttphay; mjphay = mpphay + maphay; mrphay = mkphay + mbphay; mj = mjphay*10.^-3*pi*(D*10)^2/4; mr = mrphay*pi*D^2/4; hold on; ylim([-7 7]) xlim([-4 4]) axis fill; i=0:1:180; Pkt= Pa - 0.1; Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i=180:1:345; j=1:0.0934:16.5; Pkt= (Pa*j.^n1 - 0.1); Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i=345:2:361; Pkt= 6.622*exp(-((i-359.3)/19.28).^2);; Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i=360:1:518; length(i); j=1:0.0618:delta; length(j); Pkt= (Pz./(j.^n2) - 0.1); Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; 21 goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i= 518:1:530; Pkt= 4.436e+112*i.^-41.7; Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i= 530:1:570; Pkt= 6.368e+114*i.^-42.49; Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i=570:1:720; Pkt= Pr-0.1; Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); title('Do thi phu tai tac dung len chot khuyu') ylabel('- Z'); xlabel('alpha"'); hold on; ylim([-3 7]) xlim([0 800]) axis fill; i=0:1:180; Pkt= Pa - 0.1; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i=180:1:345; 22 Pkt= 6.38e-59*i.^ 23.16; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i=345:2:361; Pkt= 6.622*exp(-((i-359.3)/19.28).^2);; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i=360:1:518; length(i); j=1:0.0618:delta; length(j); Pkt= (Pz./(j.^n2) - 0.1); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i= 518:1:660; Pkt= 0.3162*exp(-((i-492.4)/89.75).^2); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i=660:1:720; Pkt= Pr-0.1; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); title('Luc phap tuyen gay uon truc khuyu') ylabel('Z (MN)'); xlabel('GQTK'); hold on; ylim([-3 7]) xlim([0 800]) axis fill; i=0:1:180; Pkt= Pa - 0.1; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; 23 goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i=180:1:345; Pkt= 6.38e-59*i.^ 23.16; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i=345:2:361; Pkt= 6.622*exp(-((i-359.3)/19.28).^2);; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i=360:1:518; length(i); j=1:0.0618:delta; length(j); Pkt= (Pz./(j.^n2) - 0.1); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i= 518:1:660; Pkt= 0.3162*exp(-((i-492.4)/89.75).^2); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i=660:1:720; Pkt= Pr-0.1; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); title('Luc tiep tuyen lam quay truc khuyu') ylabel('T (MN)'); xlabel('GQTK'); 24 %plot(359,Pz-0.1,'r.'); hold on; ylim([-3 7]) xlim([0 800]) axis fill; i=0:15:180; Pkt= Pa - 0.1; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot([0 180],[ (Pa - 0.1) (Pa - 0.1)],'g-'); i=180:1:345; Pkt= 6.38e-59*i.^ 23.16; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); i=345:2:361; Pkt= 6.622*exp(-((i-359.3)/19.28).^2); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); i=360:1:518; length(i); j=1:0.0618:delta; length(j); Pkt= (Pz./(j.^n2) - 0.1); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); 25 i= 518:1:660; Pkt= 0.3162*exp(-((i-492.4)/89.75).^2); Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); i=660:1:720; Pkt= Pr-0.1; Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot([660 720],[ Pr-0.1 Pr-0.1],'g-'); title('Do thi luc the Pkt, luc quan tin Pj, luc tong hop') legend('Pj','Ptt','Pkt') ; ylabel('P'); xlabel('alpha'); 26 27 28 29 30 31 ... PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA CƠ KHÍ ĐỘNG LỰC - - BÁO CÁO CUỐI KỲ Chuyên ngành: Công nghệ Kỹ thuật Ơ tơ Tên đề tài TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG SVTH : BÙI CÔNG THƯ MSSV: 18145465... KỸ THUẬT TP HCM CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Bộ mơn Tính Tốn Động Cơ Đốt Trong PHIẾU NHẬN XÉT BÁO CÁO CUỐI KỲ (Dành cho giảng viên hướng dẫn) Họ tên sinh viên: Bùi... KỸ THUẬT TP HCM CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Bộ mơn Tính Tốn Động Cơ Đốt Trong PHIẾU NHẬN XÉT BÁO CÁO CUỐI KỲ (Dành cho giảng viên phản biện) Họ tên sinh viên: Bùi