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 download by : skknchat@gmail.com 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 Điểm: Chữ ký giáo viên hướng dẫn download by : skknchat@gmail.com 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: download by : skknchat@gmail.com 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-α download by : skknchat@gmail.com Đồ 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 download by : skknchat@gmail.com 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ị download by : skknchat@gmail.com Độ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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com Đồ thị biểu diễn lực N- 44 download by : skknchat@gmail.com 4.Đồ thị biểu diễn lực T- 45 download by : skknchat@gmail.com Đồ thị biểu diễn lực Z- 46 download by : skknchat@gmail.com Đồ thị T-Z 47 download by : skknchat@gmail.com Đồ thị chuyển vị piston x- 48 download by : skknchat@gmail.com Đồ thị vận tốc piston v- 49 download by : skknchat@gmail.com Đồ thị gia tốc piston j- 50 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com 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 download by : skknchat@gmail.com ... download by : skknchat@gmail.com 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... 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ị download by : skknchat@gmail.com Động cơ: KIA SORENTO 2.2L PHẦN 1: SỐ LIỆU BAN... 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: