Nghiên cứu cải thiện sự hoạt động của tổ hợp tuabin tăng áp khi động cơ diesel làm việc ở chế độ chuyển tiếp

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ GIAO THÔNG VẬN TẢI TRƯỜNG ĐẠI HỌC GIAO THÔNG VẬN TẢI TP.HCM - oOo NGUYỄN VĂN BÌNH NGHIÊN CỨU CẢI THIỆN SỰ HOẠT ĐỘNG CỦA TỔ HỢP TUABIN TĂNG ÁP KHI ĐỘNG CƠ DIESEL LÀM VIỆC Ở CHẾ ĐỘ CHUYỂN TIẾP LUẬN VĂN THẠC SĨ KỸ THUẬT NGÀ NH : KỸ THUẬT CƠ KHÍ ĐỘNG LỰC TP HỒ CHÍ MINH, 07 – 2015 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ GIAO THÔNG VẬN TẢI TRƯỜNG ĐẠI HỌC GIAO THÔNG VẬN TẢI TP.HCM - oOo NGUYỄN VĂN BÌNH NGHIÊN CỨU CẢI THIỆN SỰ HOẠT ĐỘNG CỦA TỔ HỢP TUABIN TĂNG ÁP KHI ĐỘNG CƠ DIESEL LÀM VIỆC Ở CHẾ ĐỘ CHUYỂN TIẾP CHUYÊN NGÀ NH : KỸ THUẬT CƠ KHÍ ĐỘNG LỰC MÃ SỚ : LUẬN VĂN THẠC SĨ KỸ THUẬT NGƯỜI HƯỚNG DẪN KHOA HỌC : TS LÊ VĂN VANG II   LỜI CẢM ƠN Trong suốt thời gian thực luận văn thạc sĩ này, với nỗ lực cá nhân tác giả nhận nhiều giúp đỡ nhiệt tình quý báu Tác giả trân trọng tri ân giúp đỡ Trước tiên, tác giả xin gửi lời cảm ơn sâu sắc tới thầy TS Lê Văn Vang - giảng viên trường Đại học Giao thông vận tải TPHCM Thầy hướng dẫn cung cấp cho nhiều kiến thức chuyên môn, tài liệu chuyên ngành quý báu suốt thời gian thực đề tài Thầy người hướng dẫn tận tình tâm huyết q trình thực hiên cơng trình nghiên cứu đưa ý kiến góp ý giúp tơi hồn tất luận văn Đồng thời, tác giả xin gửi lời cám ơn chân thành tới thầy, cô, đồng nghiệp bạn bè khoa Máy Tàu Thủy nói riêng thầy trường Đại học Giao thơng vận tải TPHCM nói chung truyền đạt cho kiến thức chuyên môn bổ ích q trình tơi học trường Cảm ơn Ban Chủ nhiệm khoa Máy Tàu Thủy giúp đỡ tạo điều kiện thời gian công việc để tơi hồn thành luận văn Cám ơn quan tâm, hỗ trợ động viên gia đình bạn bè giúp tơi có động lực hồn thành cơng việc nghiên cứu luận văn Tp Hồ Chí Minh, ngày 15 tháng 06 năm 2015 Tác giả Nguyễn Văn Bình III   MỤC LỤC LỜI CAM ĐOAN I  LỜI CẢM ƠN II  MỤC LỤC III  DANH MỤC CÁC KÝ HIỆU .VII  DANH MỤC HÌNH VẼ IX  DANH MỤC BẢNG BIỂU XI  LỜI NÓI ĐẦU 1  CHƯƠNG 3  GIỚI THIỆU TỔNG QUAN VỀ ĐỀ TÀI 3  1.1 Giới thiệu chung đề tài 3  1.1.1  Ảnh hưởng trình khai thác động đến môi trường 4  1.1.2 Ảnh hưởng chế độ thay đổi tải đến độ mài mòn chi tiết làm việc động 6  1.1.3 Tình hình khai thác sử dụng động Diesel Việt Nam 7  1.2 Tình hình nghiên cứu liên quan đến đề tài 8  1.3 Phương pháp mục tiêu nghiên cứu 8  CHƯƠNG 9  CƠ SỞ LÝ THUYẾT VÀ QUÁ TRÌNH HOẠT ĐỘNG CHUYỂN TIẾP CỦA ĐỘNG CƠ 9  2.1 Cơ sở lý thuyết tăng áp cho động Diesel 9  2.2 Phân tích chế độ cơng tác động Diesel tăng áp lai phụ tải 13  2.3 Các trình chuyển tiếp thay đổi thơng số động Diesel q trình chuyển tiếp 14  2.3.1 Chế độ chuyển tiếp động Diesel 14  2.3.2 Q trình đóng ngắt tải đột ngột 15  2.3.3 Quá trình tăng tốc 18  IV   2.3.4 Quá trình chuyển tiếp liên hợp 20  2.3.5 Q trình thay đổi tải có tính chu kỳ 21  2.3.6 Quá trình đảo chiều quay trục khuỷu động 22  2.4 Q trình cơng tác động chế độ chuyển tiếp 25  2.4.1 Q trình trao đổi khí động diesel tăng áp tuabin khí xả chế độ chuyển tiếp 26  2.4.2 Q trình hồ trộn hỗn hợp cháy chế độ chuyển tiếp 29  2.4.3 Sự phối hợp công tác tổ hợp tuabin máy nén động diesel trình hoạt động 31  CHƯƠNG 42  CÁC PHƯƠNG PHÁP CẢI THIỆN VÀ TÍNH TỐN SỰ THAY ĐỔI CÁC THÔNG SỐ CÔNG TÁC CỦA TỔ HỢP TUABIN TĂNG ÁP VÀ ĐỘNG CƠ TRONG CHẾ ĐỘ CHUYỂN TIẾP 42  3.1 Giới thiệu chung 42  3.2 Đối với động 43  3.2.1 Hệ thống nhiên liệu 43  3.2.2 Bộ điều tốc 43  3.2.3 Sử dụng chốt tì di động dạng thuỷ lực 44  3.2.4 Lựa chọn thơng số trao đổi khí 45  3.2.5 Cơ cấu phối khí 46  3.3 Đối với tuabin tăng áp 47  3.3.1 Sử dụng tăng áp biến áp 47  3.3.2 Giảm mômen qn tính rơ to tuabin máy nén 48  3.3.3 Điều chỉnh tuabin máy nén 48  3.3.4 Kết cấu tổ hợp tuabin máy nén tăng áp 50  3.3.5 Sử dụng tăng áp kết hợp 51  3.3.6 Sử dụng tăng áp hai cấp 52  V   3.3.7 Sử dụng tăng áp 53  3.3.8 Cấp khí bổ sung 54  3.4 Tính tốn thay đổi thơng số công tác tổ hợp tuabin tăng áp động chế độ chuyển tiếp 59  3.4.1 Lựa chọn thơng số tính tốn 59  3.4.2 Lựa chọn phương pháp tính tốn động 60  3.4.3 Lập thuật tốn chương trình tính 61  3.4.4 Phương trình thay đổi vịng quay hệ trục động - phụ tải 62  3.4.5 Suất tiêu hao khơng khí động 66  3.4.6 Nhiệt độ khí xả 68  3.4.7 Các thông số công tác máy nén 68  3.4.8 Các thơng số cơng tác tuabin khí xả 70  3.4.9 Tính tốn lượng khí cấp bổ sung cho tuabin tăng áp 74  CHƯƠNG 77  LẬP CHƯƠNG TRÌNH TÍNH TỐN CẢI THIỆN SỰ HOẠT ĐỘNG CỦA TỔ HỢP TUABIN TĂNG ÁP KHI ĐỘNG CƠ DIESEL LÀM VIỆC Ở CHẾ ĐỘ CHUYỂN TIẾP BẰNG PHƯƠNG PHÁP CẤP KHÍ BỔ SUNG 77  4.1 Lựa chọn động để tính tốn 77  4.2 Viết chương trình tính 82  4.3 Xây dựng đặc tính máy nén 83  4.4 Tính tốn q trình thay đổi tải 83  4.5 Lập chương trình tính ngôn ngữ Matlab 84  4.5.1 Giới thiệu chung phần mềm Matlab 84  4.5.2 Xây dựng đặc tính hàm lưu động tương tỉ số áp suất px/pk 85  4.5.3 Xây dựng mối quan hệ tỉ số áp suất px/pk khối lượng khí nạp Gkq 86  VI   4.5.4 Xây dựng mối quan hệ suất tiêu hao khơng khí qua động Gk với tỉ số áp suất px/pk 88  4.5.5 Xây dựng mối quan hệ lượng cấp khí bổ sung cơng suất máy nén 89  4.5.6 Xây dựng mối quan hệ lượng cấp khí bổ sung hệ số dư lượng khơng khí α1 90  4.5.7 Xây dựng mối quan hệ lượng cấp khí bổ sung hiệu suất thị động 92  4.5.8 Các tiêu đánh giá chất lượng trình thay đổi tải 93  4.5.9 Kết luận 93  KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 94  Kết luận 94  Kiến nghị hướng phát triển đề tài 94  TÀI LIỆU THAM KHẢO 95  PHỤ LỤC 96          VII   DANH MỤC CÁC KÝ HIỆU Ký hiệu Ý nghĩa Ne Cơng suất có ích động (kw) Nm Công suất giới động (kw) Ni Công suất thị động (kw) Pe Áp suất có ích bình qn (kg/cm2) Pi Áp suất thị (kg/cm2) Pm Áp suất giới (kg/cm2) D Đường kính xylanh (m) S Hành trình pittong (m) n Vòng quay động (v/ph) i Số xy lanh m Hệ số kỳ ε Tỉ số nén λta Tỉ số tăng áp δ Độ khơng đồng vịng quay Cm Tốc độ (m/s) Ta Thời gian chuyển tiếp hệ động (s) J Mơmen qn tính khối lượng quay động phụ tải (kg.ms2) φ Độ thay đổi vòng quay tương đối (rad/s) gct Lượng nhiên liệu cung cấp cho chu trình (kg/s) ge Suất tiêu hao nhiên liệu có ích (kg/kw.h) gi Suất tiêu hao nhiên liệu thị (kg/kw.h) τi Thời gian trì hỗn cháy (s) α Hệ số dư lượng khơng khí VIII   Lo Lượng khơng khí lý thuyết cần đốt cháy hoàn toàn 1kg nhiên liệu (kmol) ηi Hiệu suất thị ηm Hiệu suất giới ηe Hiệu suất có ích Mq Mơmen quay động (N.m) Mc Mômen cản động (N.m) ψ Hàm lưu động tương đối ρk Khối lượng riêng không khí nạp (kg/m3) Lk Cơng đoạn nhiệt máy nén (kj) Nk Công suất máy nén (kw) Pk Áp suất tăng áp (kg/cm2) Px Áp suất khí xả (kg/cm2) kx Chỉ số nén đoạn nhiệt khí xả kkk Chỉ số nén đoạn nhiệt khơng khí Rkk Hằng số trạng thái khơng khí (kj/kg.độ) Rkx Hằng số trạng thái khí xả (kj/kg.độ) Gk Lưu lượng khơng khí qua động (kg/s) Gkd Lượng khơng khí cịn lại xylanh đóng xu páp nạp, thải (kg/s) Gkq Lưu lượng khí nạp (kg/s) Gbs Lượng cấp khí bổ sung (kg/s) Tk Nhiệt độ sau nén có làm mát (0K) To Nhiệt độ môi trường (0K) Tx Nhiệt độ khí xả (0K) Fvan Diện tích cửa van điện từ (m2)   IX   DANH MỤC HÌNH VẼ Hình 2.1: Sơ đồ động tăng áp tua bin khí xả 12  Hình 2.2: Quá trình chuyển tiếp động 16  Hình 2.3: Sự thay đổi cơng suất tiêu thụ phụ thuộc vào vòng quay tăng tốc tàu 19  Hình 2.4: Các đường cong biểu diễn vùng làm việc động 22  Hình 2.5: Sự phụ thuộc mơmen trục khuỷu động vào vịng quay dừng tàu 23  Hình 2.6: Sự phụ thuộc vịng quay động có hộp giảm tốc vào thời gian dừng tàu 24  Hình 2.7: Đặc tính máy nén 32  Hình 2.8: Sự thay đổi thông số công tác tổ hợp tuabin máy nén làm việc trình chuyển tiếp 34  Hình 2.9: Đặc tính cơng tác đồng thời hệ máy nén tăng áp-động 36  Hình 2.10: Sự thay đổi điểm công tác đồng thời hệ máy nén tăng áp-động ổn định 38  Hình 3.1: Chốt tỳ di động thuỷ lực 45 Hình 3.2: Sơ đồ bố trí tổ hợp tuabin tăng áp có đường khí xả 50  Hình 3.3: Sơ đồ bố trí kết cấu động sử dụng tăng áp kết hợp 52  Hình 3.4: Sơ đồ bố trí kết cấu sử dụng tăng áp hai cấp 53  Hình 3.5: Hình dạng thiết kế tăng áp có cấp khí bổ sung 55  Hình 3.6: Sơ đồ bố trí họng cấp khí vào hệ thống điều khiển máy nén 55  Hình 3.7: Ảnh hưởng cấp khí bổ sung q trình chuyển tiếp 56  Hình 3.8 So sánh trình chuyển tiếp sau tăng phụ tải 57  Hình 4.1: Sơ đồ hệ thống tăng áp 79 Hình 4.2: Sơ đồ cấp khí bổ sung 79  Hình 4.3: Sơ đồ điều khiển áp suất khí nạp 80  86   Hình 4.6: Đồ thị biểu thị mối quan hệ hàm lưu lượng tương đối Ψ với tỉ số áp suất px/pk 4.5.3 Xây dựng mối quan hệ tỉ số áp suất px/pk khối lượng khí nạp Gkq Nhiệt độ khơng khí sau máy nén tăng áp: Tk = To = 25.2,25 , Pk = 2,25 bar = 2,25.104 (kg/m2); Tk = 360C = 3090K; = , = 0,254 (kg/m3 ); Gkqdm = giá trị vào công thức (3.17) ta có: = 1,82.0,756 = 1,82 (kg/s) thay tất 87   Chọn kx = 1,4 Px/Pk = 0÷1 ta có Gkq = 1,37 , , Từ chương trình tính Matlab ta lập chương trình tính M.file thể mối quan hệ Gkq Px/Pk viết theo phụ lục Từ hình 4.7 ta thấy Gkq tỉ số px/pk cắt hai điểm (0;1) khối lượng khí nạp đạt giá trị lớn khoảng tỉ số px/pk Khi thêm vào lượng Gkqdm lượng khí nạp tăng lên đáng kể đạt giá trị lớn điểm có tọa độ (0,53;0,9421) Hình 4.7: Đồ thị quan hệ tí số áp suất với khối lượng khí nạp 88   4.5.4 Xây dựng mối quan hệ suất tiêu hao khơng khí qua động Gk với tỉ số áp suất px/pk Từ phương trình (3.15) (3.16) (3.17) ta có: Gk = Gkd + Gkq ; ; , 11,92 dm ), τ = 2, nd = 720 (v/ph) = 24 (rad/s); i = 6; = 1.37 Tính tốn suất tiêu hao khơng khí qua động theo cơng thức: Gk = Gkd+ Gkq = , , , = 0,027 (kg/s) Từ chương trình tính Matlab ta lập chương trình tính M.file thể mối quan hệ Gk px/pk theo phụ lục Hình 4.8: Đồ thị quan hệ suất tiêu hao khơng khí với tỉ số áp suất 89   Từ hình 4.8 ta thấy tăng thêm lượng Gkd lượng tiêu hao khơng khí qua động tăng lên đáng kể điều có nghĩa lượng khơng khí nạp vào động tăng Gk đạt giá trị lớn điểm có tọa độ (0,53;0,9691) 4.5.5 Xây dựng mối quan hệ lượng cấp khí bổ sung cơng suất máy nén Theo cơng thức (3.23) ta có Nk = Lk.Gk trường hợp cấp khí bổ sung ta có Nk = Lk.(Gk+Gbs) = Lk.(0,9691+Gbs) Từ cơng thức (3.22) ta có: 1005 2.25 tỉ số tăng áp (Pra/Pvao) = 25.2,250,44 = 360C; lấy m = 1,8 ∆ ; (m chạy từ 1,5 tới 2); Tk nhiệt độ khơng khí sau máy nén tăng áp Tk=360C suy 31 304 thay vào cơng thức (3.22) ta tính Lk=385,18 (KJ) Công suất máy nén: Nk = Lk.Gk = 385,18.(0,9691+Gbs) Biêu diễn đồ thị Matlab ta viết M.file biểu diễn mối quan hệ Nk Gbs chương trình viết theo phụ lục Qua hình vẽ 4.9 ta thấy việc cấp khí bổ sung q trình thay đổi tải có hiệu quả, cơng suất máy nén tăng qua tăng lưu lượng khí qua máy nén qua cải thiện q trình phối hợp cơng tác tổ hợp tuabin tăng áp động tăng lượng khí Gk nạp vào động trình thay đổi tải 90   Hình 4.9: Đồ thị mối quan hệ công suất máy nén lượng cấp khí bổ sung 4.5.6 Xây dựng mối quan hệ lượng cấp khí bổ sung hệ số dư lượng khơng khí α1 Từ cơng thức (3.15) Gk = Gkd+Gkq với Gkd = 0,027 ta có Gk = 0,027 + Gkq qua hình vẽ 4.7 ta thấy Gkq đạt giá trị lớn khí Px/Pk = 0,5÷0,55 Gkq = 1,37.0,687 = 0,9421 (kg/s) suy Gk = 0,9691(kg/s) theo công thức (3.10): , Ở chế độ cấp khí bổ sung ta thay Gk Gk+Gbs = 0,9691+Gbs Gnl = ge.Ne với ge=194 (g/kW.h) thay vào biểu thức tính Gnl ta có: 91   Gnl = , 780 0,042 (kg/s) từ suy ra: , , Ta biểu thị mối quan hệ Gbs α1 Matlab với Gbs = 0,1÷1 viết chương trình tính theo phụ lục 5: Từ đồ thị ta thấy lượng cấp khí bổ sung nhiều hệ số α1 lớn việc cấp hkis bổ sung có sở làm tăng hệ số dư lượng khơng khí α1 Ta có bảng sau: Bảng 4.1: Giá trị tương ứng lượng cấp khí bổ sung Gbs hệ số α1 Gbs 0,2 0,4 0,5 0,6 0,8 0,9 α1 1,165 1,941 2,282 2,449 2,615 2,949 3,115 3,283 Hình 4.10: Đồ thị mối quan hệ lượng cấp khí bổ sung Gbs hệ số dư lượng khơng khí α1 92   4.5.7 Xây dựng mối quan hệ lượng cấp khí bổ sung hiệu suất thị động Từ công thức (3.9) tính hiệu suất thị ηi ta có: ηi = ηimax 1,065- 0,426 α1 =0,5 1,065- 0,426 α1 =0,5325- 0,213 0,9691+Gbs 0,6 Từ chương trình tính Matlab ta viết M.file biểu thị mối quan hệ Gbs ηi theo phụ lục 6: Hình 4.11: Biểu diễn mối quan hệ lượng cấp khí bổ sung Gbs hiệu suất thị ηi động 93   Qua hình 4.11 ta thấy khơng cấp khí bổ sung hiệu suất thị ηi = 0,4287 Gbs=1 ηi=0,4967 lượng cấp khí bổ sung tăng hiệu suất thị tăng việc cấp khí bổ sung q trình thay đổi tải có sở thực tế giúp cải thiện thông số động 4.5.8 Các tiêu đánh giá chất lượng trình thay đổi tải Nghiên cứu cải thiện hoạt động tổ hợp tuabin tăng áp thực sở tính tốn lý thuyết đánh giá qua mơ Bài tốn đặt điều kiện thực nghiệm tiêu đánh giá chất lượng cải thiện phương án áp dụng tiêu Các tiêu sau để đánh giá trình chuyển tiếp mà hợp lý chúng khẳng định nghiên cứu thực nghiệm tiến hành Độ dài trình chuyển tiếp động Ta thời gian (được tính giây) từ lúc bắt đầu trình chuyển tiếp đến thời điểm mà số vòng quay động thiết lập chế độ ổn định Đây tiêu quan trọng đánh giá chất lượng trình chuyển tiếp động Nó khả thích ứng tổ hợp tuabin-máy nén làm việc ổn định điều tốc với thay đổi tải động Độ thay đổi vòng quay tương đối tính theo vịng quay chế độ tức thời hệ so với vòng quay định mức: Đây tiêu quan trọng đánh giá mức độ làm việc ổn định động 4.5.9 Kết luận - Kết tính tốn mơ phù hợp với qui luật thay đổi thông số mà đề tài nêu - Cấp khí bổ sung cho tổ hợp tuabin tăng áp thời điểm thay đổi tải có sở khoa học cải thiện thông số công tác động hoạt động chế độ 94   KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN Kết luận Những kết nhận từ q trình nghiên cứu tính tốn, mơ việc cấp khí bổ sung cho tổ hợp tuabin tăng áp rút số kết luận sau: - Trong thực tế khai thác, trình thay đổi tải trình làm việc thường xuyên động Do đề tài phân tích đưa giải pháp cải thiện hoạt động động diesel, tổ hợp tuabin tăng áp chế độ chuyển tiếp nghiên cứu phương pháp cụ thể phương pháp cấp khí bổ sung - Phương pháp cấp khí bổ sung cho hệ thống tăng áp động Diesel tăng áp tuabin khí xả giải pháp nhằm tăng tính tin cậy, tính kinh tế giảm thiểu độc tố khí xả thải mơi trường khai thác động Diesel tăng áp hệ cũ, hệ động sử dụng nhiều Việt Nam - Xây dựng mơ hình tính tốn kỹ thuật cho phép tính tốn mơ biến thiên thơng số công tác tổ hợp tuabin tăng áp, động Diesel tải thay đổi Đây sở để nhà chế tạo, nhà quản lý khai thác động đề giải pháp hợp lý giúp cải thiện hoạt động tổ hợp tuabin tăng áp động diesel làm việc chế độ thay đổi tải Kiến nghị hướng phát triển đề tài Hướng nghiên cứu đề tài chế tạo hệ thống tự động cấp khí bổ sung cho hệ thống tăng áp động Diesel làm việc chế độ chuyển tiếp thực nhờ vào tiến khoa học công nghệ thông tin tự động hố Từ áp dụng trực tiếp động Diesel có trang bị hệ thống tăng áp tuabin khí xả 95   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] [9] Constantine D Rakopoulos, Evangelos G Giakoumis,  “Diesel Engine Transient Operation” Springer Publisher (2009) Constantine D Rakopoulos, Evangelos G Giakoumis, “Comparative first- and second-law parametric study of transient diesel engine operation” Elsevire Publisher (2005) 96   PHỤ LỤC LẬP CHƯƠNG TRÌNH TÍNH BẰNG MATLAB PHỤ LỤC % 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 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') PHỤ LỤC % 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 97   % CHON DIEU KIEN BIEN BAN DAU Gkqdm=1.82; % khoi luong quet dinh muc kx=1.4; % he so mu doan nhiet cua xa pk=2.25; % ap suat dong ham truoc tuabin rk=0.254; % khoi luong rieng khong 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 nap Gkq') title ('DO THI QUAN HE TI SO AP SUAT VOI KHOI LUONG KHI NAP') axis ([0 1]) grid on PHỤ LỤC % 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 lai xilanh Gkqdm=1.82; % khoi luong quet dinh muc kx=1.4; % he so mu doan nhiet cua xa pk=2.25; % ap suat dong ham truoc tuabin rk=0.254; % khoi luong rieng khong tai cua nap y=sqrt(2*kx/(kx-1))*(sqrt(x.^(2/kx)-x.^((kx+1)/kx))); Gkq=Gkqdm*y*sqrt(rk*pk) Gk=Gkd+Gkq 98   plot(x,Gk,'linewidth',2) Xlabel('ti so ap suat px/pk') Ylabel('suat tieu hao khong Gk') title ('DO THI QUAN HE TI SO AP SUAT VOI SUAT TIEU HAO KHONG KHI') axis ([0 1]) grid on PHỤ LỤC % 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 bo sung Gbs=0:0.01:1; Nmn=385.18*(0.9691+Gbs) plot(Gbs,Nk,'linewidth',2) Xlabel ('luong cap 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 300 800]) grid on PHỤ LỤC % 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 bo sung x=0:0.01:1; y=(0.9691+x)/0.6 plot(x,y,'linewidth',2) 99   Xlabel ('luong cap 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 4]) grid on PHỤ LỤC % 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 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 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 0.4 0.6]) grid on TP HỒ CHÍ MINH, 06 - 2017 LUẬN VĂN ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC GIAO THƠNG VẬN TẢI THÀNH PHỐ HỒ CHÍ MINH Cán hướng dẫn khoa học : TS PHẠM THỊ NGA Cán chấ m nhận xét 1: PGS.TS HỒ THỊ THU HÒA Cán chấ m nhận xét 2: TS TRƯƠNG QUANG DŨNG Luận văn bảo vệ Hội đồng trường Đại Học Giao Thơng Vận Tải TP Hồ Chí Minh ngày 29 tháng 06 năm 2017 Thành phần Hội đồng đánh giá luận văn Thạc sĩ gồm: TS Nguyễn Văn Khoảng Chủ tịch hội đồng; PGS.TS Hồ Thị Thu Hòa Ủy viên phản biện; TS Trương Quang Dũng Ủy viên phản biện; TS Lý Bách Chấn Ủy viên; TS Lê Phúc Hòa Ủy viên thư ký; Xác nhận Chủ tịch Hội đồng đánh giá luận văn Trưởng Khoa quản lý chuyên ngành sau luận văn sửa chữa CHỦ TỊCH HỘI ĐỒNG TS NGUYỄN VĂN KHOẢNG TRƯỞNG KHOA KINH TẾ TS NGUYỄN VĂN KHOẢNG