Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 44 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
44
Dung lượng
566,41 KB
Nội dung
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 : PHAN TRẦN KHẢI MSSV: 18145370 SVTH : NGUYỄN THÀNH ĐẠT MSSV: 18145334 SVTH : NGUYỄN HỮU LƯƠNG MSSV: 18145392 GVHD: TS NGUYỄN VĂN TRẠNG Tp Hồ Chí 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 : PHAN TRẦN KHẢI MSSV: 18145370 SVTH : NGUYỄN THÀNH ĐẠT MSSV: 18145 SVTH : NGUYỄN HỮU LƯƠNG MSSV: 18145392 GVHD: TS NGUYỄN VĂN TRẠNG Tp Hồ Chí 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 Độc lập - Tự – Hạnh phúc Bộ mơn Tính Toá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ư MSSV: 18145465 Hội đồng:………… Họ tên sinh viên: Phan Trần Khải MSSV: 18145370 Hội đồng:………… Họ tên sinh viên: Nguyễn Thành Đạt MSSV: 18145334 Hội đồng:………… Họ tên sinh viên: Nguyễn Hữu Lương MSSV: 18145392 Hội đồng:………… Tên đề tài:Tí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, tí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 Điểm Điểm đạt tối đa 30 Đúng format với đầy đủ cả hình thức và nô ̣i dung của các mục 10 Mục tiêu, nhiệm vụ, tổng quan đề tài 10 Tính cấp thiết đề tài 10 Nội dung 50 Khả ứng dụng kiến thức toán học, khoa học kỹ thuật, khoa học xã hơ ̣i… Khả thực hiện/phân tích/tổng hợp/đánh giá 10 Khả thiết kế chế tạo hệ thống, thành phần, quy trình 15 đá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 15 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 10 Sản phẩm cụ thể 10 Tổng điểm 100 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 Độc lập – Tự – Hạnh phú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ư MSSV: 18145465 Hội đồng:………… Họ tên sinh viên: Phan Trần Khải MSSV: 18145370 Hội đồng:………… Họ tên sinh viên: Nguyễn Thành Đạt MSSV: 1814 Hội đồng:………… Họ tên sinh viên: Nguyễn Hữu Lương MSSV: 18145392 Hội đồng:………… Tên đề tài:Tí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: (Mã 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, tí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 Điểm Điểm đạt tối đa 30 Đúng format với đầy đủ cả hình thức và nô ̣i dung của các mục 10 Mục tiêu, nhiệm vụ, tổng quan đề tài 10 Tính cấp thiết đề tài 10 Nội dung 50 Khả ứng dụng kiến thức toán học, khoa học kỹ thuật, khoa học xã hơ ̣i… Khả thực hiện/phân tích/tổng hợp/đánh giá 10 Khả thiết kế, chế tạo hệ thống, thành phần, quy trình 15 đá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 15 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 10 Sản phẩm cụ thể 10 Tổng điểm 100 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 vii 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 plot(0.1126,-0.1913,'b.'); plot(0.1686,-0.3209,'b.'); plot(0.1451,-0.2753,'b.'); plot(1.2173,-0.3718,'b.'); plot(1.1583,-0.4406,'b.'); hold off; %%%%%%%%%%%% TRON % % % % % la la la la la thetaS=15 theta1=26 theta4=22 theta2=48 theta3=56 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... HCM CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Độc lập - Tự – Hạnh phú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... HCM CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Độc lập – Tự – Hạnh phú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