(TIỂU LUẬN) bài tập lớn môn học TÍNH TOÁN ĐỘNG cơ đốt TRONG đề tài KIA SORENTO (d)

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP HỒ CHÍ MINH KHOA ĐÀO TẠO CHẤT LƯỢNG CAO BÀI TẬP LỚN MƠN HỌC: TÍNH TỐN ĐỘNG CƠ ĐỐT TRONG Đề Tài: KIA SORENTO (D) GVHD: PGS.TS Lý Vĩnh Đạt Thành viên: Nguyễn Xuân Trường MSSV 19145054 Phạm Duy Minh MSSV 19145268 Nguyễn Quang Sự MSSV 19145301 Tp.HCM, Tháng 12 Năm 2021 CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc PHIẾU NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… …………………………………………………………………………………………………… Tp Hồ Chí Minh, ngày tháng năm 2021 Chữ ký giáo viên hướng dẫn Điểm: MỤC LỤC PHẦN 1: SỐ LIỆU BAN ĐẦU I Các Thông Số Cho Trước Của Động Cơ PHẦN 2: TÍNH TỐN NHIỆT ĐỘNG CƠ ĐỐT TRONG I CHỌN CÁC THƠNG SỐ TÍNH TỐN NHIỆT Áp suất khơng khí nạp: Nhiệt độ khơng khí nạp: Áp suất khí nạp trước xuppap nạp: Nhiệt độ khí nạp trước xuppap nạp: Áp suất cuối trình nạp: Áp suất khí sót: Nhiệt độ khí sót: Độ tăng nhiệt độ khí nạp mới: Chọn hệ số nạp thêm: 10 Chọn hệ số quét buồng cháy: 11 Chọn hệ số hiệu đỉnh tỷ nhiệt: 12 Hệ số nhiệt lợi dụng điểm z: 13 Hệ số nhiệt lợi dụng điểm b: 14 Chọn hệ số dư lượng khơng khí: 15 Chọn hệ số điền đầy đồ thị công: 16 Tỷ số tăng áp: II TÍNH TỐN NHIỆT Quá trình nạp: Quá trình nén: Quá trình cháy: Quá trình giãn nở: 4.1 Hệ số giãn nở sớm ρ: 4.2 Hệ số giãn nở sau δ: 4.3 Chỉ số giãn nở đa biến trung bình n2: 4.4 4.5 Nhiệt độ cuối trình giãn nở Tb: Áp suất cuối trình giãn nở Pb: 4.6 5.Tính Tốn Các Thơng Số Đặc Trưng Của Chu Trình Kiểm nghiệm nhiệt độ khí sót Tr: 5.1Áp suất chỉ thị trung bình tính tốn: 5.2Áp suất chỉ thị trung bình thực tế 5.3Áp suất tổn thất khí Pm 5.4Áp suất có ích trung bình Pe 5.5Hiệu suất giới 5.6Hiệu suất chỉ thị 5.7Hiệu suất có ích 5.8Tính suất tiêu hao nhiên liệu chỉ thị gi 5.9Tính suất tiêu hao nhiên liệu ge 5.10 Tính tốn thơng số kết cấu động 6.Vẽ đồ thị công chỉ thị PHẦN 3: TÍNH TỐN ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC CƠ CẤU PISTON - TRỤC KHUỶU - THANH TRUYỀN I ĐỘNG HỌC CỦA PISTON 1.Chuyển vị Piston 2.Tốc độ Piston 3.Gia tốc Piston II ĐỘNG LỰC HỌC CỦA CƠ CẤU PISTON - KHUỶU TRỤC - THANH TRUYỀN Lực khí thể Lực quán tính chi tiết chuyển động Hệ lực tác dụng cấu trục khuỷu - truyền PHẦN 4: BẢN VẼ CÁC ĐỒ THỊ - CODE ĐỒ THỊ I ĐỒ THỊ 1.Đồ thị công P-V Đồ thị công P-α Đồ thị biểu diễn lực N- α 4.Đồ thị biểu diễn lực T-α Đồ thị biểu diễn lực Z-α 46 Đồ thị T-Z 47 48 49 50 Đồ thị chuyển vị piston x-α Đồ thị vận tốc piston v-α Đồ thị gia tốc piston j-α II CODE ĐỒ THỊ 51 Phân chia công việc Tên nhiệm vụ Tra thông số ban đầu, thông số hiệu chỉnh Chọn thông số tính tốn nhiệt Tính tốn nhiệt Vẽ đồ thị P-V, P5 Tính tốn động lực học Vẽ đồ thị Động cơ: KIA SORENTO 2.2L PHẦN 1: SỐ LIỆU BAN ĐẦU I Các Thông Số Cho Trước Của Động Cơ STT Tên Thông số Loại động Số xylanh Số kỳ Dung tích xylanh Cơng suất động Số vịng quay động Momen xoắn Tỷ số nén Phun nhiên liệu sớm 10 Góc mở sớm xúpap nạp 11 Góc đóng muộn xúpap nạp 12 Góc mở sớm xúpap xả 13 Góc đóng muộn xúpap xả 14 Hành trình piston 15 Đường kính xy lanh PHẦN 2: TÍNH TỐN NHIỆT ĐỘNG CƠ ĐỐT TRONG I CHỌN CÁC THƠNG SỐ TÍNH TỐN NHIỆT Áp suất khơng khí nạp: Áp suất khơng khí nạp chọn áp suất khí quyển: p0 = 0,1 (MN/m2) Nhiệt độ khơng khí nạp: Nước ta thuộc khu vực nhiệt đới, nhiệt độ trung bình ngày chọn Tkk = 29oC, T0 = (tkk + 273)oK = (29 + 273)oK = 302oK Áp suất khí nạp trước xuppap nạp: Động Diesel bốn kì tăng áp ta có pk > p0 pk = 0,2 (MN/m2) Nhiệt độ khí nạp trước xuppap nạp: Vì động kỳ tăng áp nên ta có: T =T( ) k −1 o - T = 355,5 K m Chọn: m = 1,5  T m = 25 oK Áp suất cuối trình nạp: Động Disel bốn kì khơng tăng áp: pa = (0,88 ÷ 0,98).pk (MN/m2), ta chọn: Pa= 0,9.Pk= 0,9.0,2=0,18 (MN/m2) Áp suất khí sót: Đối với động Disel: pr = (0,106 ÷0,115 ), có tốc độ cao nên chọn: Pr = 0,11 (MN/m2) Nhiệt độ khí sót: 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 Disel Tr = (700o ÷ 900)oK, ta chọn: Tr =750o K Độ tăng nhiệt độ khí nạp mới: Khí nạp chuyển động đường ống nạp vào xilanh động Disel tiếp xúc với vách nóng lên sấy nóng lên trị số nhiệt độ ΔT Đối với động Disel ∆T = (20 ÷ 40)oK, ta chọn ∆T = 25oK Chọn hệ số nạp thêm: Hệ số nạp thêm λ1 phụ thuộc chủ yếu vào pha phối khí thông thường chọn giới hạn λ1 = 1,02 ÷ 1,07, ta chọn: λ1 = 1,06 10 Chọn hệ số quét buồng cháy: Đối với động Disel tăng áp ta chọn: λ2 = 0.2 11 Chọn hệ số hiệu đỉnh tỷ nhiệt: Hệ số hiệu đính tỷ nhiệt λt phụ thuộc vào thành phần hỗn hợp ⍺ nhiệt độ khí sót Tr Đối với động Disel có ⍺ = 1,7 thì: λt = 1,11 12 Hệ số nhiệt lợi dụng điểm z: Đối với động Diesel thì, ξ z = (0,65 ÷ 0,85), ta chọn: ξz = 0,8 13 Hệ số nhiệt lợi dụng điểm b: Hệ số lợi dụng nhiệt điểm b phụ thuộc vào tốc độ động tỷ số nén Khi động hoạt động chế độ tồn tải tốc độ động thấp, q trình cháy rớt giảm nên dẫn đến ξblớn, từ ta chọn: ξb = 0,9 14 Chọn hệ số dư lượng khơng khí: Đối với động Disel có buồng cháy dự bị, α = (1,7 ÷ 2,2) , ta chọn 15 α = 1,7 Chọn hệ số điền đầy đồ thị công: Hệ số điền đầy đồ thị công đánh giá phần hao hụt diện tích đồ thị cơng thực tế so với đồ thị cơng tính tốn, có buồng cháy ngăn cách φ d = (0,92 ÷0,96) ta chọn : φd = 0,94 16 Tỷ số tăng áp: Tỷ số tăng áp tỷ số áp suất hỗn hợp khí xi lanh cuối trình cháy & trình nén: Đối với động Diesel = = (1,35 ÷ 2,4) Ta chọn: × =2 Đồ thị biểu diễn lực N- 44 4.Đồ thị biểu diễn lực T- 45 Đồ thị biểu diễn lực Z- 46 Đồ thị T-Z 47 Đồ thị chuyển vị piston x- 48 Đồ thị vận tốc piston v- 49 Đồ thị gia tốc piston j- 50 II CODE ĐỒ THỊ phis = 15; phi1 = 12; phi2 = 28; phi3 = 45; phi4 = 14; S = 0.096; B = 0.854; R = 0.048; Sx = (pi)*(B^2)/4; Vd = Sx*0.096; Vc = Vd/(16-1); Va = Vd + Vc; lamda = 0.25; n = 3800; pa = 0.18; n1 = 1.38; n2 = 1.2; ro = 1.07; Vz = ro*Vc; pz = 16.5; pr = 0.11; %hieu chinh as=[0 10 26]; ps=[0.104 0.114 0.18]; a1=linspace(0,26,1000); X1=R*(1-cosd(a1)+lamda*(1-cosd(2.*a1))./4); 51 V1=X1.*Sx+Vc; p1=interp1(as,ps,a1,'pchip'); %phan nap a2=linspace(26,180,1000); X2=R*(1-cosd(a2)+lamda*(1-cosd(2.*a2))./4); V2=X2.*Sx+Vc; p2=linspace(pa,pa,1000); %phan nen a3=linspace(180,350,1000); X3=R*(1-cosd(a3)+lamda*(1-cosd(2.*a3))./4); V3=X3.*Sx+Vc; p3=pa.*(Va./V3).^n1; %nen - chay as = [350 355 360]; % goc danh lua som toi 360 ps=[6.87 7.7 10];% ap suat tuong ung a4=linspace(350,360,1000); X4=R*(1-cosd(a4)+lamda*(1-cosd(2.*a4))./4); V4=X4.*Sx+Vc; p4=interp1(as,ps,a4,'pchip'); % phan bo cong chay - gian no as=[360 365 370]; ps=[10 14.5 pz]; a5=linspace(360,370,500); X5=R*(1-cosd(a5)+lamda*(1-cosd(2.*a5))./4); V5=X5.*Sx+Vc; p5=interp1(as,ps,a5,'pchip'); %%%%%%%%%% 52 as=[370 375 380]; ps=[pz 15 11.2]; a6=linspace(370,380,500); X6=R*(1-cosd(a6)+lamda*(1-cosd(2.*a6))./4); V6=X6.*Sx+Vc; p6=interp1(as,ps,a6,'pchip'); %phan gian no a7=linspace(380,490,500); X7=R*(1-cosd(a7)+lamda*(1-cosd(2.*a7))./4); V7=X7.*Sx+Vc; p7=pz.*(Vz./V7).^n2; %gian no - thai as=[490 515 540]; ps=[0.7854 0.505 0.3676]; a8=linspace(490,540,500); X8=R*(1-cosd(a8)+lamda*(1-cosd(2.*a8))./4); V8=X8.*Sx+Vc; p8=interp1(as,ps,a8,'pchip'); %%%%%%%%% as=[540 555 575]; ps=[0.3676 0.11 0.11]; a9=linspace(540,572,1000); X9=R*(1-cosd(a9)+lamda*(1-cosd(2.*a9))./4); V9=X9.*Sx+Vc; p9=interp1(as,ps,a9,'pchip'); %ph7an thai a10=linspace(572,720,1000); 53 X10=R*(1-cosd(a10)+lamda*(1-cosd(2.*a10))./4); V10=X10.*Sx+Vc; p10=linspace(pr,pr,1000); %%%%%%%%%% a=[a1 a2 a3 a4 a5 a6 a7 a8 a9 a10];%do X=[X1 X2 X3 X4 X5 X6 X7 X8 X9 X10];%m V=[V1 V2 V3 V4 V5 V6 V7 V8 V9 V10];%m3 p=[p1 p2 p3 p4 p5 p6 p7 p8 p9 p10];%MN/m2 plot(V,p,'k','linewidth',2); title('DO THI P-V'); xlabel('The tich V (m3)'); ylabel('Ap suat P (MN/m2)'); grid on; figure(2) plot(a,p,'k','linewidth',2); title (' DO THI P - PHI '); xlabel('Goc quay truc khuyu (do)'); ylabel('Ap suat P (MN/m2)'); grid on; figure(3) plot(a,X,'k','linewidth',2); grid on title('DO THI CHUYEN VI'); xlabel('Goc quay truc khuyu (do)'); ylabel('Chuyen vi x (m)') ; figure(4) Vp=R*((2*pi*n)/60).*(sind(a)+lamda*(sind(2.*a))./2); 54 plot(a,Vp,'k','linewidth',2); grid on title('DO THI VAN TOC PISTON'); xlabel('Goc quay truc khuyu (do)'); ylabel('Van toc v (m/s) '); figure(5) J=(R*((2*pi*n)/60)^2).*(cosd(a)+lamda.*cosd(2.*a));%m/s2 plot(a,J,'k','linewidth',2); grid on title('DO THI GIA TOC PISTON'); xlabel('Goc quay truc khuyu (do)'); ylabel('Gia toc j (m/s2) '); figure(6) Pj=(-0.0003605*R*((2*pi*n)/60)^2).*(cosd(a)+lamda.*cosd(2.*a)); Pkt=p-0.1; P1=Pkt+Pj;%hop luc plot(a,Pkt,'r','linewidth',2); hold on plot(a,Pj,'k','linewidth',2); hold on plot(a,P1,'g','linewidth',2); grid on 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'); figure(7) 55 T=P1.*sind(a+asind(lamda.*sind(a)))./cosd(asind(lamda.*sind(a))); plot(a,T,'k','linewidth',2); grid on title('DO THI T'); xlabel('Goc quay truc khuyu (do)'); ylabel('T (MN/m2)'); figure(8) Z=P1.*cosd(a+asind(lamda.*sind(a)))./cosd(asind(lamda.*sind(a))); plot(a,Z,'k','linewidth',2); grid on title('DO THI Z'); xlabel('Goc quay truc khuyu (do)'); ylabel('Z (MN/m2)'); figure(9) N=P1.*tand(asind(lamda.*sind(a))); plot(a,N,'k','linewidth',2); grid on title('DO THI LUC NGANG'); xlabel('Goc quay truc khuyu (do)'); ylabel('N (MN/m2)'); figure(10) plot(T,Z,'k','linewidth',2); axis ij; ax = gca; ax.XAxisLocation='top'; ax.YAxisLocation='right'; grid on 56 title('DO THI VECTO PHU TAI'); xlabel('T (MN/m2)'); ylabel('Z (MN/m2)'); figure(11) A1 = [a1 a2]; P1 = [p1 p2]; Pj1=(-0.0003605*R*((2*pi*n)/60)^2).*(cosd(A1)+lamda.*cosd(2.*A1)); Pkt1=P1-0.1; P11=Pkt1+Pj1; A2 = [a3 a4]; P2 = [p3 p4]; Pj2=(-0.0003605*R*((2*pi*n)/60)^2).*(cosd(A2)+1/4.*cosd(2.*A2)); Pkt2=P2-0.1; P12=Pkt2+Pj2; A3 = [a5 a6 a7 a8]; P3 = [p5 p6 p7 p8]; Pj3=(-0.0003605*R*((2*pi*n)/60)^2).*(cosd(A3)+1/4.*cosd(2.*A3)); Pkt3=P3-0.1; P13=Pkt3+Pj3; A4 = [a9 a10]; P4 = [p9 p10]; Pj4=(-0.0003605*R*((2*pi*n)/60)^2).*(cosd(A4)+1/4.*cosd(2.*A4)); Pkt4=P4-0.1; P14=Pkt4+Pj4; T1=P11.*sind(A1+asind(lamda.*sind(A1)))./cosd(asind(lamda.*sind(A1))); T2=P12.*sind(A2+asind(lamda.*sind(A2)))./cosd(asind(lamda.*sind(A2))); T3=P14.*sind(A4+asind(lamda.*sind(A4)))./cosd(asind(lamda.*sind(A4))); 57 T4=P13.*sind(A3+asind(lamda.*sind(A3)))./cosd(asind(lamda.*sind(A3))); plot(A1,T1.*R,'k','linewidth',2); hold on plot(A2-180,T2.*R,'g','linewidth',2); hold on plot(A4-540,T3.*R,'r','linewidth',2); hold on plot(A3-360,T4.*R,'b','linewidth',2); hold on plot(A1,(T1+T2+T3+T4).*R,'y','linewidth',2); grid on title('DO THI MOMEN QUAY TRUC KHUYU'); xlabel('Goc quay truc khuyu (do)'); ylabel('Mq (MN/m2)'); legend('May1','May2','May3','May4','Tong'); 58 ... Chọn thơng số tính tốn nhiệt Tính tốn nhiệt Vẽ đồ thị P-V, P5 Tính tốn động lực học Vẽ đồ thị Động cơ: KIA SORENTO 2.2L PHẦN 1: SỐ LIỆU BAN ĐẦU I Các Thông Số Cho Trước Của Động Cơ STT Tên Thông... Ta có bảng số liệu sau: 24 PHẦN 3: TÍNH TỐN ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC CƠ CẤU PISTON - TRỤC KHUỶU - THANH TRUYỀN I ĐỘNG HỌC CỦA PISTON Hình 3.1 Sơ đồ Động học cấu Piston – Khuỷu trục – Thanh truyền... LIỆU BAN ĐẦU I Các Thông Số Cho Trước Của Động Cơ PHẦN 2: TÍNH TOÁN NHIỆT ĐỘNG CƠ ĐỐT TRONG I CHỌN CÁC THƠNG SỐ TÍNH TỐN NHIỆT Áp suất khơng khí nạp: