BÀI TẬP LỚN MÔN HỌC TÍNH TOÁN ĐỘNG CƠ đốt TRONG đề TÀI FORD ECOSPORT

TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP.HCM KHOA CƠ KHÍ ĐỘNG LỰC BỘ MƠN ĐỘNG CƠ BÀI TẬP LỚN MƠN HỌC TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG ĐỀ TÀI: FORD ECOSPORT Giảng viên hướng dẫn: PGS.TS Lý Vĩnh Đạt Lớp: Sáng thứ 2, 3; Tiết 1,2,3 Nhóm sinh viên thực hiện: Võ Minh Duy - 19145354 Nguyễn Minh Đức - 19145367 Nguyễn Ngọc Lâm Huy - 19145112 Tp Hồ Chí Minh, Tháng 8, năm 2021 TRƯỜNG ĐẠI HỌC SP KT TP.HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Độc Lập - Tự Do - Hạnh Phúc BỘ MÔN ĐỘNG CƠ NHIỆM VỤ BÀI TẬP LỚN MƠN TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG Họ Tên SV: Nguyễn Minh Đức Võ Minh Duy MSSV: 19145367 MSSV: 19145354 Nguyễn Ngọc Lâm Huy MSSV: 19145112 LỚP: Sáng Thứ 2,3 Tiết 123 Nhóm 11  Số liệu ban đầu Loại động cơ: 1.0L Ecoboost DOHC 12 van i3, Turbo Số kỳ, : kỳ Cơng suất có ích, Pmax (kW): 125Hp 93,2 kWSố vòng quay, n (vòng/phút): 6675 v/ph Tỉ số nén, : 10:1 Hệ số dư lượng khơng khí,(): 0.29 Làm mát bằng: dung dịch  Nội dung thuyết minh 2.1 Tính tốn nhiệt xây dựng giản đồ cơng thị động 2.2 Tính toán động lực học cấu piston – trục khuỷu – truyền  Nội dung vẽ 3.1 Bản vẽ đồ thị công thị P – V 3.2 Bản vẽ đồ thị P – , PJ, P1 3.3 Bản vẽ đồ thị quảng đường SP, vận tốc VP, gia tốc JP piston Ngày giao nhiệm vụ : tháng 7/2021 Ngày hoàn thành : 8/8/2021 GIÁO VIÊN HƯỚNG DẪN CHÍNH (Ký ghi rõ họ tên) PGS.TS Lý Vĩnh Đạt THỜI GIAN THỰC HIỆN  Tiến độ thực Tính tốn nhiệt: tuần - Chọn thơng số đầu vào: tuần - Tính tốn nhiệt: tuần - Vẽ đồ thị công thị: tuần Tính tốn động học, động lực học vẽ vẽ động học động lực học: tuần - Tính tốn động học động lực học: tuần - Vẽ vẽ động học động lực học: tuần Viết thuyết minh chuẩn bị báo cáo: tuần  Nội dung qui định: Bản vẽ đồ thị công thị, động lực học, Thuyết minh viết khổ giấy A4 (in file) Nội dung thuyết minh gồm: - Nhiệm vụ tập lớn GV hướng dẫn thông qua - Phần số liệu tính tốn nhiệt, động học động lực học có kèm theo đồ thị - Kết luận tập lớn MỤC LỤC MỞ ĐẦU NỘI DUNG Chương 1: TÍNH TỐN NHIỆT A Các thông số cho trước động cơ: B Chọn thơng số tính tốn nhiệt: Áp suất khơng khí nạp ( Po ) Nhiệt độ khơng khí nạp ( T0 ): 3 Nhiệt độ khí nạp trước xupap nạp ( Pk ): Nhiệt độ khí nạp trước xuppap nạp ( Tk ): Áp suất cuối trình nạp ( Pa ): Chọn áp suất khí sót ( Pr ) Nhiệt độ khí sót ( Tr ): Độ tăng nhiệt độ khí nạp ( T ): .4 Chọn hệ số nạp thêm ( 1 ): 10 Chọn hệ số quét buồng cháy ( 2 ): 11 Chọn hệ số hiệu đính tỷ nhiệt ( t ): .5 12 Hệ số lợi dụng nhiệt điểm Z (  Z ): 13 Hệ số lợi dụng nhiệt điểm b ( b ): 14 Hệ số dư lượng khơng khí (  ): 15 Chọn hệ số điền đầy đồ thị công ( d ): 16 Tỷ số tăng áp: B Tính tốn nhiệt: Quá trình nạp: Quá trình nén: Quá trình cháy: C Tính thơng số đặc trưng chu trình: 13 Áp suất thị trung bình tính tốn: 13 Áp suất thị trung bình thực tế: 13 Áp suất tổn thất khí: .13 Áp suất có ích trung bình: 13 Hiệu suất giới: 14 Hiệu suất thị: 14 Hiệu suất có ích: 14 Tính suất tiêu hao nhiên liệu thị: 14 Tính suất tiêu hao nhiên liệu: 14 D Vẽ đồ thị công: 14 Chương 2: TÍNH TỐN ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC 19 A Động học piston: 19 Chuyển vị piston: 19 Tốc độ piston: 19 Gia tốc piston: .19 B Động lực học cấu trục khuỷu truyền: 20 Lực khí thể: 20 Chương 3: ĐỒ THỊ 22 Chương 4: CODE 25 MỞ ĐẦU Ơ tơ ngày sử dụng rộng rãi nước ta phương tiện lại cá nhân vận chuyển hành khách , hàng hoá phổ biến Sự gia tăng nhanh chóng số lượng ơtơ xã hội , đặc biệt loại ôtô đời kéo theo nhu cầu đào tạo lớn nguồn nhân lực phục vụ ngành công nghiệp ôtô linh vực thiết kế Sau học xong giáo trình ‘Động đốt trong2’ chúng em giao nhiệm vụ làm đồ án môn học Qua đồ án giúp sinh viên chúng em nắm lực tác dụng , công suất động điều kiện đảm bảo bền số nhóm chi tiết ơtơ , máy kéo Vì thiết thực với sinh viên nghành cơng nghệ kỹ thuật ơtơ Động nhóm chúng em lựa chọn Động 1.0L Ecoboost DOHC 12 van i3, Turbo Ford Một hãng xe thông dụng lựa chọn nhiều nước ta với động Ecosport tiên tiến với tiện ích tiết kiệm nhiên liệu Ecosport với công nghệ tiên tiến kỹ sư Ford trang bị lên động như:phun nhiên liệu trực tiếp, trục cam biến thiên, turbo tăng áp…Bộ điều khiển trung tâm giúp động vận hành êm Nhằm để hiểu công suất động điều kiện bền nên nhóm em định chọn Ecosport cho việc tính tốn Vì bước đầu làm quen với cơng việc tính tốn , thiết kế ơtơ nên khơng tránh khỏi bỡ ngỡ vướng mắc.Nhưng với quan tâm, động viên , giúp đỡ, hướng dẫn tận tình thầy giáo hướng dẫn, giáo viên giảng dạy thầy giáo khoa nên chúng em cố gắng để hoàn thành đồ án thời gian giao Nhóm em xin chân thành cảm ơn ạ! NỘI DUNG Chương 1: TÍNH TỐN NHIỆT A Các thông số cho trước động cơ: - Loại động cơ: Động 1.0L Ecoboost DOHC 12 van i3, Turbo - Kiểu, loại động cơ, số kỳ  : Động thẳng hàng, động xăng, tăng áp, kỳ - Số xilanh: xilanh - Cách bố trí xilanh: 1-2-3 - Đường kính xilanh, D(mm): 71,9mm - Hành trình piston, S(mm): 82mm - Dung tích động cơ: 998cc - Công suất thiết kế, N e (kW): 125Hp= 93,2kW - Số vòng quay thiết kế, n (v/ph): 6675v/ph - Tỷ số nén  : 10:1 - Kiểu buồng cháy phương pháp tạo hỗn hợp: Buồng cháy kiểu SprayGuide, tạo hỗn hợp phân lớp - Kiểu làm mát: Làm mát dung dịch - Chiều dài truyền, L(mm): 137mm - Khối lượng nhóm piston, mnp : 109,3 kg / m - Khối lượng nhóm truyền, m tt : 136,7 kg / m - Thể tích công tác xilanh: Vh  - Thông số kết cấu   0,998  0,33(l ) R  0, 29 L B Chọn thơng số tính tốn nhiệt: Áp suất khơng khí nạp ( Po ) Áp suất khơng khí nạp chọn áp suất khí quyển, giá trị Po phụ thuộc vào độ cao so với mực nước biển Càng lên cao Po giảm khơng khí lỗng, độ cao so với mực nước biển: Po  0,1013 MN / m Nhiệt độ khơng khí nạp ( T0 ): Nhiệt độ khơng khí nạp phụ thuộc chủ yếu vào nhiệt độ trung bình mơi trường, nơi xe sử dụng Điều khó khăn xe thiết kế để sử dụng vùng có khoảng biến thiên nhiệt độ ngày lớn Miền Nam nước ta thuộc vực nhiệt đới, nhiệt độ trung bình ngày chọn t kk  29 C cho khu vực miền Nam, đó: T0  (tkk  273)  302 K Nhiệt độ khí nạp trước xupap nạp ( Pk ): Động bốn kỳ tăng áp: p k áp suất khí nạp nén sơ cấp trước máy nén tăng áp bơm quét khí Pk  (0,14  0, 4) ta chọn Pk  0,14 Mpa Nhiệt độ khí nạp trước xuppap nạp ( Tk ): Đối với động bốn kỳ tăng áp khơng có làm mát trung gian Tk xác định công thức: P  Tk  T0  k   P0  m1 m Trong đó: m - số nén đa biến trung bình khí nén, phụ thuộc vào loại máy nén ( m  1,5 1,65 ), thông thường chọn 1,4 Vậy ta có: Tk  302  0,14   0,1013  1,4 1 1,4  331 K Áp suất cuối trình nạp ( Pa ): Đối với động tăng áp: Pa   0, 88  0, 98  Pk Chọn Pa  0,88 Pk  0,88.0,14  0,1232 Mpa Chọn áp suất khí sót ( Pr ) Pr  (1, 05  1,1) Pth Pth  1, 02  1, 04  P0 Chọn Pth  1, 02 P0  1, 02.0,1013  0,1033 Mpa Chọn Pr  1, 05 Pth  1, 05.0,1033  0,1085 Mpa Nhiệt độ khí sót ( Tr ): Nhiệt độ khí sót Tr phụ thuộc vào chủng loại đơng cơ.Nếu q trình giản nở triệt để, nhiệt độ Tr thấp: Đối với động xăng: Tr  900  1000K Chọn Tr  900K Độ tăng nhiệt độ khí nạp ( T ): Khí nạp chuyển động đường ống nạp vào xylanh động tiếp xúc với vách nóng nên sấy nóng lên trị số nhiệt độ T Khi tiến hành tính tốn nhiệt động người ta thường chọn trị số T vào số liệu thực nghiệm 23 24 Chương 4: CODE A Code đồ thị P-V động lực học: TSN = 10; S = 82; R = S/2; lamda = 0.29; Vd = 1.0/3; Vc = Vd/(TSN-1); Va = TSN*Vc; Fp = (pi*0.719^2)/4; % qua trinh nap a1 = linspace(6,180,100); X1 = 0.01*R*((1 - cosd(a1)) + (lamda*(1 - cosd(2*a1)))/4); Va1 = X1*Fp + Vc; Pa1 = linspace(0.1232,0.1232,100); %plot(Va1,Pa1) % qua trinh nen a den c' a2 = linspace(180,346,100); X2 = 0.01*R*((1 - cosd(a2)) + (lamda*(1 - cosd(2*a2)))/4); Vc1 = X2*Fp + Vc; Pc1 = 0.1232.*((Va./Vc1).^1.369); % c' den c'' 25 a3 = linspace(346,360,100); Vcn = 0.01*R*((1 - cosd(360)) + (lamda*(1 - cosd(2*360)))/4)*Fp + Vc; Pcn = 0.1232.*((Va./Vcn).^1.369); Pz1 = Pcn*3.5; % diem c' V1c = 0.01*R*((1 - cosd(346)) + (lamda*(1 - cosd(2*346)))/4)*Fp + Vc; P1c = 0.1232.*((Va./V1c).^1.369); % diem c'' P2c = (Pz1 - Pcn)/3 + Pcn ; V2c = Vc; % tong diem P3c = [P1c,2.88,P2c]; V3c = [V1c,0.0402,V2c]; Vxc = linspace(V1c,V2c,100); Pxc = interp1(V3c,P3c,Vxc,'spline'); % chay hon hop c'' den z'' a31 = linspace(360,376.85,100); Vz2 = 0.01*R*((1 - cosd(376)) + (lamda*(1 - cosd(2*376)))/4)*Fp + Vc; % diem z' P1z = Pz1; V1z = (Vz2 - Vc)/2 + Vc; % diem z'' V2z = 0.01*R*((1 - cosd(376.85)) + (lamda*(1 - cosd(2*376.85)))/4)*Fp + Vc; 26 P2z = Pz1.*((Vz2./V2z).^1.21); % tong diem z V3z = [Vc,0.0385,0.04,V1z,0.043,V2z]; P3z = [P2c,7.1,8.37,P1z,9.9,P2z]; Vxz = linspace(Vc,V2z,100); Pxz = interp1(V3z,P3z,Vxz,'spline'); % chay gian no z'' den b' a4 = linspace(376.85,488,100); X4 = 0.01*R*((1 - cosd(a4)) + (lamda*(1 - cosd(2*a4)))/4); Vb = X4*Fp + Vc; Pb = Pz1.*((Vz2./Vb).^1.21); % qua trinh xa b' den b'' a5 = linspace(488,540,100); Pb1 = Pz1.*((Vz2/Va).^1.21); % diem b' V1b = 0.01*R*((1 - cosd(488)) + (lamda*(1 - cosd(2*488)))/4)*Fp + Vc; P1b = Pz1.*((Vz2./V1b).^1.21); % diem b'' V2b = Va; P2b = (Pb1 - 0.1085)/2 + 0.1085; % tong diem V3b = [V1b,0.3,0.4,0.5,0.6,V2b]; P3b = [P1b,0.8,0.7,0.6,0.5,P2b]; Vxb = linspace(V1b,V2b,100); 27 Pxb = interp1(V3b,P3b,Vxb,'spline'); % b'' den r' a6 = linspace(540,560,100); % diem r' V1r = 0.01*R*((1 - cosd(560)) + (lamda*(1 - cosd(2*560)))/4)*Fp + Vc; P1r = 0.1085; V4b = [V2b,0.37,0.368,V1r]; P4b = [P2b,0.3,0.2,P1r]; Vx1b = linspace(V2b,V1r,100); Px1b = interp1(V4b,P4b,Vx1b,'spline'); % duong Pr Pr1 = linspace(0.1085,0.1085,100); ar1 = linspace(560,720,100); Xr1 = 0.01*R*((1 - cosd(ar1)) + (lamda*(1 - cosd(2*ar1)))/4); Vr1 = Xr1*Fp + Vc;% r den Pa % diem r V2r = Vc; P2r = 0.1085; % diem r'' V3r = 0.01*R*((1 - cosd(5)) + (lamda*(1 - cosd(2*5)))/4)*Fp + Vc; P3r = 0.17; % diem a' V1a = 0.01*R*((1 - cosd(7)) + (lamda*(1 - cosd(2*7)))/4)*Fp + Vc; 28 P1a = 0.2; % tong diem Vra = [V2r,V3r,V1a]; Pra = [P2r,P3r,P1a]; Vxra = linspace(V2r,V1a,100); Pxra = interp1(Vra,Pra,Vxra,'spline'); figure(1) V = [Va1,Vc1,Vxc,Vxz,Vb,Vxb,Vx1b,Vr1,Vxra]; p = [Pa1,Pc1,Pxc,Pxz,Pb,Pxb,Px1b,Pr1,Pxra]; plot(V,p); axis square title('DO THI P - V'); xlabel('THE TICH V (dm3)'); ylabel('AP SUAT P (MN/m2)'); hold on grid on figure % dong luc hoc a1 = linspace(0,180,100); mt = 15*95.4/80 + (1/3)*(477/16); P = [Pa1,Pc1,Pxc,Pxz,Pb,Pxb,Px1b,Pr1]; Pkt = P - 0.1; alpha = [a1,a2,a3,a31,a4,a5,a6,ar1]; 29 plot(alpha,Pkt,'b'); title('DO THI Pkt - Pj - P1'); xlabel('GOC QUAY TRUC KHUYU (DO)'); ylabel('Pkt - Pj - P1 (MN/m2)'); hold on Pj = -mt*R*((pi*4000/30))^2*0.00000001*(cosd(alpha) + lamda*cosd(2*alpha)); plot(alpha,Pj,'g'); hold on P1 = Pkt + Pj; plot(alpha,P1,'r'); legend('Pkt','Pj','P1'); grid on ax = gca; ax.XAxisLocation = 'origin'; ay = gca; ay.YAxisLocation = 'origin'; figure(3) 30 B Code động học: TSN = 17.5; S = 82; B = 71.9; Vd = 999/4; Vc = Vd/(TSN-1); Va = TSN*Vc; n1 = 1.369; n2 = 1.21; pa = 0.1232; pc = 2.88; pb = 0.623; pr = 0.1033; pz = 10.1; lamda = 0.29; R = S/2000; %mm Fp = (pi*(B/1000)^2)/4; %m^2 n = 6675; %v/p w = 2*pi*n/60; mnp = 109.3; %kg/m^2 mtt = 136.7; %kg/m^2 mk = 136.7; %kg/m^2 mA = mtt/3; 31 mB = 2*mtt/3; mt = mnp + mA; mr = mk + mB; %Hieu chinh a tu 0-6 a01 = linspace(0,6,30); p01j = -mt*R*w^2*(cosd(a01)+lamda*cosd(2*a01))*10^-6; %MN/m^2 xhc = R*((1-cosd(3))+lamda*(1-cosd(2*3))/4); Vhc = (xhc*pi*B^2/4)+Vc; Xr1 = R*((1-cosd(6))+lamda*(1-cosd(2*6))/4); Vr1 = (Xr1*pi*B^2/4)+Vc; pr1 = pa; ahc = [0,3,6]; phc = [0.095,pr,pr1]; pkt01 = interp1(ahc,phc,a01,'spline'); p101 = p01j + pkt01; %Qua trinh nap a1 = linspace(6,180,170); p1j = -mt*R*w^2*(cosd(a1)+lamda*cosd(2*a1))*10^-6; pkt1 = linspace(pa,pa,170); p11 = p1j + pkt1; %Qua trinh nen: a2 = linspace(180,346,160); p2j = -mt*R*w^2*(cosd(a2)+lamda*cosd(2*a2))*10^-6; 32 xc1 = R*((1-cosd(a2))+lamda*(1-cosd(2*a2))/4);%m vc1 = (xc1*pi*B^2/4)+Vc; pkt2 = (pa*(Va.^n1*vc1.^-n1)); p12 = p2j + pkt2; %Hieu chinh a tu 340-360 a23 = linspace(346,360,40); p23j = -mt*R*w^2*(cosd(a23)+lamda*cosd(2*a23))*10^-6; Xc1 = R*((1-cosd(346))+lamda*(1-cosd(2*346))/4); Vc1 = (Xc1*pi*B^2/4)+Vc; pc1 = pa*(Va.^n1*Vc1.^-n1); pz1 = pz; pcz1 = pz1 - pc; pcc2 = pcz1/3; pc2 = pcc2 + pc; ahc1 = [346,350,353,357,360]; phc1 = [pc1,7.9,8.75,9.5,pc2]; a11 = linspace(346,360,40); pkt11 = interp1(ahc1,phc1,a11,'spline'); p123 = p23j + pkt11; %Qua trinh chay a3 = linspace(360,364,10); p3j = -mt*R*w^2*(cosd(a3)+lamda*cosd(2*a3))*10^-6; Vz = 1.3*Vc; 33 ro = 1; Vz1 = Vz/ro; pz2 = 16; pkt3 = linspace(pc2,pz2,10); p13 = p3j + pkt3; %Hieu chinh a tu 364-376.85 a34 = linspace(364,376.85,20); p34j = -mt*R*w^2*(cosd(a34)+lamda*cosd(2*a34))*10^-6; ahc2 = [364,370,376.85]; phc2 = [pz2,pz1,pz2]; a22 = linspace(364,376.85,20); pkt22 = interp1(ahc2,phc2,a22,'spline'); p134 = p34j + pkt22; %Qua trinh gian no a4 = linspace(376.85,460,90); p4j = -mt*R*w^2*(cosd(a4)+lamda*cosd(2*a4))*10^-6; xz2 = R*((1-cosd(370))+lamda*(1-cosd(2*370))/4); vz2 = (xz2*pi*B^2/4)+Vc; xb1 = R*((1-cosd(a4))+lamda*(1-cosd(2*a4))/4); vb1 = (xb1*pi*B^2/4)+Vc; pkt4 = (pz2*vz2.^n2*vb1.^-n2); p14 = p4j + pkt4; %Hieu chinh a tu 460-540 34 a45 = linspace(460,540,80); p45j = -mt*R*w^2*(cosd(a45)+lamda*cosd(2*a45))*10^-6; Xb1 = R*((1-cosd(460))+lamda*(1-cosd(2*460))/4); Vb1 = (Xb1*pi*B^2/4)+Vc; pb1 = pz2*vz2.^n2*Vb1.^-n2; ahc3 = [460,490,505,520,540]; phc3 = [pb1,0.37,0.34,0.32,pr]; a33 = linspace(460,540,80); pkt33 = interp1(ahc3,phc3,a33,'spline'); p145 = p45j + pkt33; %Qua trinh thai a5 = linspace(540,720,200); p5j = -mt*R*w^2*(cosd(a5)+lamda*cosd(2*a5))*10^-6; pkt5 = linspace(pr,pr,200); p15 = p5j + pkt5; figure(2); % Do thi chuyen vi cua piston alpha = [a01,a1,a2,a23]; X = R*((1-cosd(alpha))+(lamda*(1-cosd(2*alpha)))/4); plot(alpha,X,'linewidth',1); hold on; X1 = R*(1-cosd(alpha)); plot(alpha,X1,'linewidth',1); 35 hold on; X2 = R*((lamda*(1-cosd(2*alpha)))/4); plot(alpha,X2,'linewidth',1); xlabel('Alpha (Do)'); ylabel('Sp (m)'); title('DO THI CHUYEN VI CUA PISTON'); legend('Sp','Sp1','Sp2'); ax = gca; ax.XAxisLocation = 'origin'; grid on figure(3); % Do thi van toc cua piston V = R*w*(sind(alpha)+lamda*sind(2*alpha)/2); plot(alpha,V,'linewidth',1); hold on; V1 = R*w*sind(alpha); plot(alpha,V1,'linewidth',1); hold on; V2 = R*w*lamda*sind(2*alpha)/2; plot(alpha,V2,'linewidth',1); xlabel('Alpha (Do)'); ylabel('Vp (m/s)'); title('DO THI VAN TOC CUA PISTON'); 36 legend('V','V1','V2'); ax = gca; ax.XAxisLocation = 'origin'; grid on figure(4); % Do thi gia toc cua piston J = R*w^2*(cosd(alpha)+lamda*cosd(2*alpha)); plot(alpha,J,'linewidth',1); hold on; J1 = R*w^2*cosd(alpha); plot(alpha,J1,'linewidth',1); hold on; J2 = R*w^2*lamda*cosd(2*alpha); plot(alpha,J2,'linewidth',1); hold on xlabel('GQTK (Do)'); ylabel('jp (m/s^2)'); title('DO THI GIA TOC CUA PISTON'); legend('jp','jp1','jp2'); ax = gca; ax.XAxisLocation = 'origin'; grid on 37 ... DOHC 12 van i3, Turbo Ford Một hãng xe thông dụng lựa chọn nhiều nước ta với động Ecosport tiên tiến với tiện ích tiết kiệm nhiên liệu Ecosport với công nghệ tiên tiến kỹ sư Ford trang bị lên động... trình thải - Dựng đường cong nén: Trong hành trình nén, khí xilanh bị nén với số nén đa biến trung bình n1 , từ phương trình: Pa Va n1  Pxn Vxn n1  const Trong đó: Pa ,Va áp suất thể tích khí... ĐỘNG LỰC Độc Lập - Tự Do - Hạnh Phúc BỘ MÔN ĐỘNG CƠ NHIỆM VỤ BÀI TẬP LỚN MƠN TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG Họ Tên SV: Nguyễn Minh Đức Võ Minh Duy MSSV: 19145367 MSSV: 19145354 Nguyễn Ngọc Lâm Huy MSSV: