Hướng nghiên cứu tiếp theo của đề tài là có thể chế tạo được hệ thống tự động cấp khí bổ sung cho hệ thống tăng áp khi động cơ Diesel làm việc ở chế độ chuyển tiếp được thực hiện nhờ vào tiến bộ của khoa học công nghệ thông tin và tự động hoá. Từ đó có thể áp dụng trực tiếp trên những động cơ Diesel có trang bị hệ thống tăng áp bằng tuabin khí xả.
TÀI LIỆU THAM KHẢO I. Tiếng Việt
[1]. Lê Viết Lượng (2005), “Các Chế Độ Chuyển Tiếp Của Động Cơ Diesel”, NXB Hải Phòng.
[2]. TS. Lê Văn Vang, TS. Trương Thanh Dũng, Ths. Hoàng Văn Sĩ
(2012), “Bài Giảng Động Cơ Diesel Tàu Thủy”, Đại Học Giao Thông Vận Tải Tp.HCM.
[3]. TS. Lê Văn Vang (2013), “Khai Thác Hệ Động Lực Tàu Thủy”, Đại Học Giao Thông Vận Tải Tp.HCM.
[4]. PGS. TS Lê Hữu Sơn,TS (2012), “Động lực học hệđộng lực tàu thủy”. NXB giao thông vận tải.
[5]. Lê Viết Lượng (2000), “Lý Thuyết Động Cơ Diesel”, NXB giáo dục. [6]. Nguyễn Văn Sơn (2002), “Tua Bin Tăng Áp”, Đại Học Giao Thông
Vận Tải Tp.HCM.
[7]. Nguyễn Tất Tiến (2003), “Nguyên Lý Động Cơ Đốt Trong”, NXB giáo dục.
II. Tiếng Anh
[8]. Constantine D. Rakopoulos, Evangelos G. Giakoumis, “Diesel Engine Transient Operation” Springer Publisher (2009).
[9]. Constantine D. Rakopoulos, Evangelos G. Giakoumis, “Comparative first- and second-law parametric study of transient diesel engine operation” Elsevire Publisher (2005)
PHỤ LỤC
LẬP CHƯƠNG TRÌNH TÍNH BẰNG MATLAB 1. PHỤ LỤC 1
% CHUONG TRINH TINH MOI QUAN HE LUU LUONG VOI TI SO AP SUAT % Goi x = px/pk; phi = y x=0:0.01:1; % TRUONG HOP kx1 = 1.4 kx1=1.4; y1=sqrt(2*kx1/(kx1-1))*(sqrt(x.^(2/kx1)-x.^((kx1+1)/kx1))); % TRUONG HOP kx2 = 1.3 kx2=1.3; y2=sqrt(2*kx2/(kx2-1))*(sqrt(x.^(2/kx2)-x.^((kx2+1)/kx2))); % TRUONG HOP kx3 = 1.135 kx3=1.135; y3=sqrt(2*kx3/(kx3-1))*(sqrt(x.^(2/kx3)-x.^((kx3+1)/kx3))); plot(x,y1,x,y2,x,y3,'linewidth',2) grid on
Xlabel ('ti so ap suat (px/pk)[-]')
Ylabel ('ham luu luong khi xa (phi[-])')
title ('DO THI QUAN HE TI SO AP SUAT VOI HAM LUU LUONG') legend ('kx1 = 1.4','kx2 = 1.3','kx3 = 1.1.35','location','southeast')
2. PHỤ LỤC 2
% CHUONG TRINH TINH BIEU DIEN MOI QUAN HE TI SO AP SUAT VA KHOI LUONG KHI NAP
% Gia thiet x=px/pk; x=0:0.01:1
% CHON DIEU KIEN BIEN BAN DAU Gkqdm=1.82; % khoi luong khi quet dinh muc kx=1.4; % he so mu doan nhiet cua khi xa pk=2.25; % ap suat dong ham truoc tuabin
rk=0.254; % khoi luong rieng khong khi tai cua nap y=sqrt(2*kx/(kx-1))*(sqrt(x.^(2/kx)-x.^((kx+1)/kx))); Gkq=Gkqdm*y*sqrt(rk*pk)
plot(x,Gkq,'linewidth',2) Xlabel('ti so ap suat px/pk') Ylabel('khoi luong khi nap Gkq')
title ('DO THI QUAN HE TI SO AP SUAT VOI KHOI LUONG KHI NAP') axis ([0 1 0 1])
grid on
3. PHỤ LỤC 3
% CHUONG TRINH TINH BIEU DIEN MOI QUAN HE TI SO AP SUAT VA SUAT TIEU HAO KHONG KHI
% Gia thiet x=px/pk; x=0:0.01:1
% CHON DIEU KIEN BIEN BAN DAU
Gkd=0.027 % khoi luong khong khi con lai trong xilanh Gkqdm=1.82; % khoi luong khi quet dinh muc
kx=1.4; % he so mu doan nhiet cua khi xa pk=2.25; % ap suat dong ham truoc tuabin
rk=0.254; % khoi luong rieng khong khi tai cua nap y=sqrt(2*kx/(kx-1))*(sqrt(x.^(2/kx)-x.^((kx+1)/kx))); Gkq=Gkqdm*y*sqrt(rk*pk)
plot(x,Gk,'linewidth',2) Xlabel('ti so ap suat px/pk')
Ylabel('suat tieu hao khong khi Gk')
title ('DO THI QUAN HE TI SO AP SUAT VOI SUAT TIEU HAO KHONG KHI')
axis ([0 1 0 1]) grid on
4. PHỤ LỤC 4
% CHUONG TRINH TINH MOI QUAN HE GIUA CONG SUAT MAY NEN VA LUONG CAP KHI BO SUNG
% Goi Nk la cong suat may nen; Gbs la luong cap khi bo sung Gbs=0:0.01:1;
Nmn=385.18*(0.9691+Gbs) plot(Gbs,Nk,'linewidth',2)
Xlabel ('luong cap khi bo sung Gbs') Ylabel ('cong suat may nen Nk')
title ('DO THI QUAN HE HE GIUA CONG SUAT MAY NEN VA LUONG CAP KHI BO SUNG')
axis ([0 1 300 800]) grid on
5. PHỤ LỤC 5
% CHUONG TRINH TINH MOI QUAN HE GIUA HIEU SUAT CHI THI VA LUONG CAP KHI BO SUNG
% Goi y=alpha la hieu suat chi thi; x=Gbs la luong cap khi bo sung x=0:0.01:1;
y=(0.9691+x)/0.6 plot(x,y,'linewidth',2)
Xlabel ('luong cap khi bo sung x') Ylabel ('he so alpha y')
title ('DO THI QUAN HE HE SO ALPHA VA LUONG CAP KHI BO SUNG')
axis ([0 1 1 4]) grid on
6. PHỤ LỤC 6
% CHUONG TRINH TINH MOI QUAN HE GIUA HIEU SUAT CHI THI VA LUONG CAP KHI BO SUNG
% Goi ni la hieu suat chi thi; x=Gbs la luong cap khi bo sung x=0:0.01:1;
ni=0.5325-0.213./(sqrt((0.9691+x)/0.6)).^3; plot(x,ni,'linewidth',2)
Xlabel ('luong cap khi bo sung Gbs') Ylabel ('hieu suat chi thi ni')
title ('DO THI QUAN HE HIEU SUAT CHI THI VA LUONG CAP KHI BO SUNG')
axis ([0 1 0.4 0.6]) grid on