Nghiên cứu đặc tính bôi trơn nhiệt thủy động của ổ có dạng đầu to thanh truyền257

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI Phạm Trung Thiên NGHIÊN CỨU ĐẶC TÍNH BƠI TRƠN NHIỆT THỦY ĐỘNG CỦA Ổ CĨ DẠNG ĐẦU TO THANH TRUYỀN LUẬN ÁN TIẾN SĨ KỸ THUẬT CƠ KHÍ Hà Nội – 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI Phạm Trung Thiên NGHIÊN CỨU ĐẶC TÍNH BƠI TRƠN NHIỆT THỦY ĐỘNG CỦA Ổ CÓ DẠNG ĐẦU TO THANH TRUYỀN Ngành : Kỹ thuật khí Mã số : 9520103 LUẬN ÁN TIẾN SĨ KỸ THUẬT CƠ KHÍ NGƯỜI HƯỚNG DẪN KHOA HỌC TS TRẦN THỊ THANH HẢI TS PHẠM MINH HẢI Hà Nội – 2020 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Cơng trình thực Bộ mơn Máy Ma sát học - Viện Cơ khí, Trường Đại học Bách khoa Hà Nội hướng dẫn TS Trần Thị Thanh Hải TS Phạm Minh Hải Các số liệu kết nghiên cứu luận án trung thực chưa cơng bố cơng trình khác Hà Nội, ngày THAY MẶT TẬP THỂ HƯỚNG DẪN TS Trần Thị Thanh Hải tháng năm Người cam đoan Phạm Trung Thiên i LỜI CẢM ƠN Để hoàn thiện Luận án Tiến sĩ này, bên cạnh cố gắng thân, Tôi nhận động viên giúp đỡ nhiều Thầy Cô giáo, nhà Khoa học, đồng nghiệp bạn bè Tôi xin bày tỏ lòng biết ơn sâu sắc đến TS Trần Thị Thanh Hải TS Phạm Minh Hải Th.S Lưu Trọng Thuận người tận tình hướng dẫn, định hướng, đào tạo giúp đỡ suốt q trình nghiên cứu hồn thành luận án Xin chân thành cảm ơn Thầy, Cô giáo Bộ môn Máy Ma sát học Trường ĐHBKHN giảng dạy, bảo tạo điều kiện giúp đỡ thời gian làm nghiên cứu hoàn thành luận án Xin cảm ơn đồng nghiệp Khoa Cơ khí, lãnh đạo trường Đại học Kinh tế - Kỹ thuật Công nghiệp tạo điều kiện giúp đỡ thời gian ủng hộ để tơi hồn thành luận án Cuối cùng, tơi xin bày tỏ kính trọng, biết ơn lòng yêu thương tới đại gia đình, bạn bè thực động viên, giúp đỡ suốt thời gian học tập Trường Đại học Bách khoa Hà Nội Hà Nội, ngày tháng năm Nghiên cứu sinh Phạm Trung Thiên ii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT TT Ký hiệu, viết tắt Ý nghĩa C Khe hở bán kính DAQ Bộ xử lý tín hiệu d1 Đường kính đầu nhỏ d2 Đường kính ngồi đầu nhỏ ĐCT Điểm chết ĐCD Điểm chết e Độ lệch tâm n Số khoảng chia theo chiều dài Vector pháp tuyến c bề mặt ổ 10 m Số khoảng chia theo chu vi 11 lđt Chiều dài đầu to truyền 12 p Áp suất thủy động 13 PLC Bộ lập trình 14 R Bán kính vịng quay trục khuỷu 15 Rb Bán kính bạc 16 Rt Bán kính trục 17 s chiều dài hành trình pít-tơng 18 W 19 20 x, z Tọa độ theo phương chu vi chiểu dài φ Góc chất tải 21 θ Toạ độ trụ hệ Oxyz 22 ω t, ωb Tải trọng Vận tốc góc trục bạc iii MỤC LỤC LỜI CAM ĐOAN I LỜI CẢM ƠN II DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT III MỤC LỤC IV DANH MỤC BẢNG BIỂU VII DANH MỤC HÌNH ẢNH, ĐỒ THỊ VIII PHẦN MỞ ĐẦU 1 TÍNH CẤP THIẾT CỦA ĐỀ TÀI ĐỐI TƯỢNG NGHIÊN CỨU PHẠM VI NGHIÊN CỨU MỤC TIÊU NGHIÊN CỨU NỘI DUNG NGHIÊN CỨU PHƯƠNG PHÁP NGHIÊN CỨU Ý NGHĨA KHOA HỌC VÀ THỰC TIỄN NHỮNG ĐIỂM MỚI CỦA ĐỀ TÀI CẤU TRÚC LUẬN ÁN CHƯƠNG TỔNG QUAN VỀ BÔI TRƠN THỦY ĐỘNG CỦA Ổ CÓ DẠNG ĐẦU TO THANH TRUYỀN 1.1 Ổ ĐẦU TO THANH TRUYỀN 1.1.1 Khái niệm 1.1.2 Các Hiện tượng, nguyên nhân hư hỏng 1.2 ĐỘNG CƠ ĐỐT TRONG 1.3 TÌNH HÌNH NGHIÊN CỨU TRONG NƯỚC VÀ TRÊN THẾ GIỚI 1.3.1.Tình hình nghiên cứu giới 1.3.1.1 Nghiên cứu lý thuyết bôi trơn ổ đầu to truyền 1.3.1.2 Nghiên cứu th ực nghiệm bôi trơn ổ đầu to truyền 14 1.3.2.Tình hình nghiên cứu nước 18 KẾT LUẬN CHƯƠNG 20 CHƯƠNG CƠ SỞ LÝ THUYẾT BÔI TRƠN THỦY ĐỘNG Ổ ĐẦU TO THANH TRUYỀN VÀ MÔ PHỎNG SỐ NHIỆT ĐỘ MÀNG DẦU Ổ ĐẦU TO THANH TRUYỀN TRONG THIẾT BỊ THỰC NGHIỆM 22 iv 2.1 LÝ THUYẾT BÔI TRƠN THỦY ĐỘNG 22 2.1.1 Phương trình Reynolds 22 2.1.2 Phương trình chiều dày màng dầu 23 2.1.3 Phương trình Reynolds cho ổ đỡ thủy động 24 2.1.4 Phương trình cân tải 26 2.1.5 Phương trình lượng 27 2.2 MÔ PHỎNG NHIỆT ĐỘ MÀNG DẦU TRONG Ổ ĐẦU TO THANH TRUYỀN 28 2.2.1 Điều kiện biên 28 2.2.2 Mơ hình phần tử hữu hạn cho nhiệt độ màng dầu 30 2.2.3 Mô trường nhiệt độ màng dầu bôi trơn ổ đầu to truyền 32 2.2.4 Kết mô trường nhiệt độ ổ 35 KẾT LUẬN CHƯƠNG 39 CHƯƠNG PHƯƠNG PHÁP NGHIÊN CỨU VÀ HỆ THỐNG THIẾT BỊ THỰC NGHIỆM 40 3.1 PHƯƠNG PHÁP NGHIÊN CỨU 40 3.2 THIẾT BỊ THỰC NGHIỆM 41 3.2.1 Nguyên lý hoạt động 41 3.2.2 Các cụm chi tiết điển hình 44 3.2.2.1 Thanh truyền 44 3.2.2.2 Cụm pít-tơng truyền dẫn 46 3.2.2.3 Cơ cấu tạo tải 48 3.2.2.4 Hệ thống thủy lực 50 3.3 HỆ THỐNG ĐIỀU KHIỂN VÀ HỆ THỐNG ĐO 51 3.3.1 Hệ thống điều khiển 51 3.3.1.1 Biến Tần 51 3.3.1.2 Động 53 3.3.1.3 Tủ điều khiển 53 3.3.2 Hệ thống đo 55 3.3.2.1 Hệ thống đo lực 55 3.3.2.2 Hệ thống đo áp suất 58 3.3.2.3 Hệ thống đo nhiệt độ màng dầu 64 KẾT LUẬN CHƯƠNG 69 CHƯƠNG KẾT QUẢ NGHIÊN CỨU 70 4.1 KẾT QUẢ THỰC NGHIỆM 70 4.1.1 Xử lý số liệu thực nghiệm 70 v 4.1.2 Kết đo tải tác dụng lên truyền 72 4.1.3 Kết đo áp suất màng dầu ổ đầu to truyền 75 4.1.4 Kết đo nhiệt độ màng d ầu ổ đầu to truyền 78 4.2 KẾT QUẢ MÔ PHỎNG SỐ NHIỆT ĐỘ MÀNG DẦU Ổ ĐẦU TO THANH TRUYỀN 84 4.2.1 Kết mô độ chênh nhiệt độ ổ với dầu Besil F100 84 4.2.2 Kết mô độ chênh nhiệt độ màng dầu ổ với dầu Atox 320 86 4.3 SO SÁNH NHIỆT ĐỘ MÀNG DẦU MÔ PHỎNG VÀ THỰC NGHIỆM 89 4.3.1 So sánh kết thực nghiệm với kết mô nhiệt độ màng dầu với dầu Besil F100 89 4.3.2 So sánh kết thực nghiệm với kết mô dầu Atox 320 91 KẾT LUẬN CHƯƠNG 94 KẾT LUẬN VÀ KIẾN NGHỊ 95 DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 97 TÀI LIỆU THAM KHẢO 98 PHỤ LỤC 105 Phụ lục 1: Quy trình đo tải tác dụng lên truyền 105 Phụ lục 2: Quy trình đo áp suất màng dầu ổ đầu to truyền 107 Phụ lục 3: Quy trình đo nhiệt độ màng dầu ổ đầu to truyền 109 Phụ lục 4: Code chương trình tính tốn mơ số 111 vi DANH MỤC BẢNG BIỂU Bảng 1: Thống kê phương pháp thực giới khoảng 20 năm gần Bảng 1: Thông số ổ đỡ 32 Bảng 2: Thông số dầu 33 Bảng 1: Các thông số kỹ thuật thiết b ị thực nghiệm 43 Bảng 2: Thông số kỹ thu ật cảm biếp áp suất XCQ-062 60 Bảng Giá trị nhiệt độ 20 mẫu cảm biến góc 00 truyền 100 vg/ph góc 3600 trục khuỷu 70 Bảng Nhóm mẫu số liệu phân tích phương sai ANOVA 72 Bảng Tổng hợp số liệu thực nghiệm nhiệt độ 93 vii DANH MỤC HÌNH ẢNH, ĐỒ THỊ Hình 1.1: Các phận truyền Hình 1.2: Một số hư hỏng thường gặp ổ đầu to truyền Hình 1.3: Cấu tạo động đốt trong[1] Hình 1.4: Mơ hình rời r ạc ổ đầu to truyền [34] 11 Hình 1.5: Ứng suất truyền chịu tác dụng lực siết bu-lơng [35] 11 Hình 1.6: Biến dạng ổ đầu to truyền chịu tác dụng lực siết bu-lông [36] 11 Hình 1.7: Hình dạng ban đầu bạc lót lực xiết bu lơng [37] 12 Hình 1.8: Hình dạng bạc lót chu kỳ hoạt động [35] 12 Hình 1.9: Mơ hình ổ đầu to truyền 2D [41] 13 Hình 1.10: Sự dịch chuyển bạc lót ổ đầu to truyền [43] 13 Hình 1.11: (a) Tiếp xúc trịn xoay (b) Mơ hình ổ đầu to truyền [44] 14 Hình 1.12: ổ đầu to truyền Pierre-Eugene[49] 15 Hình 1.13: Thiết bị nghiên cứu bôi trơn[50] 16 Hình 1.14: Thanh truyền hệ thống đo đặc tính bơi trơn đầu to [59] 17 Hình 1.15: Băng thử để khảo sát bôi trơn ổ đầu đo truyền[40] 17 Hình 1.16: Sơ đồ lực tác dụng lên ổ đầu to truyền [62] 18 Hình 1.17: Ứng suất truyền [63] 18 Hình 1.18: Trượt tương đối bạc lót ổ [63] 18 Hình 2.1: Hệ trục tọa độ 22 Hình 2.2: Mặt cắt ổ đỡ 24 Hình 2.3: Miền khai triển ổ 25 Hình 2.4: Sơ đồ cân lực tác dụng lên truyền 26 Hình 2.5: Hệ trục toạ độ 27 Hình 2.6: Điều kiện biên 29 Hình 2.7: Miền tích phân màng dầu 30 Hình 2.8: Phép chuyển hệ toạ độ 30 viii [49] Piere-Eugene J.(1983), “Contribution l’Etude de la Déformation Elastique d’un Coussinet de Tête de Bielle en Fonctionnement Hydrodynamique Permanent”, Thèse de Doctorat de l’Université de Poitiers [50] Optasanu V.(2000), “Modélisation Expérimentale et Numérique de la Lubrification des Paliers Compliants sous Chargement Dynamique”, Thèse de Doctorat de l’Université de Poitiers [51] Hoang L.V.,(2002), “Modélisation Expérimentale de la Lubrification Thermoélastohydrodynamique des Paliers de Tête de Bielle Comparaison entre les Résultats Théoriques et Expérimentaux”, Thèse de Doctorat de l’Université de Poitiers [52] Hoang L V., Souchet D., Bonneau D.(2002), "Connecting-rod big end bearing thermoelastohydrodynamic lubrication (TEHD) - Comparison between theory and experiment", Int Journal of Applied Mech and Eng., vol 7, p 231-236 [53] M’hammed El Gadari ,n, Aurelian Fatu ,, Mohamed Hajjam.(2016), " Shaft roughness effect Experimentation on elasto-hydrodynamic lubrication and numerical simulation", of rotary journal lip seals: homepage: www.elsevier.com/locate/triboint [54] Cooke W.L., (1965-1966), “Dynamic Displacement in a Diesel Engine Main Bearing”, Proceeding Lubrication and Wear Second Convention, Instn Mech Engrs., Vol 23 [55] Rosenberg R.C., (1973), “A Method for Determining the Influence of Multigrade oils on Journal Bearing Performance”, SEA TRANS Paper 730483, Vol 82 [56] Goodwin G., Holmes R., (1975), “Determination of the Oil Film Thickness in a Crankshaft Main Bearing”, The Journal of Automotive Engineering, Instn, Mech, Engrs., 1975 [57] Bates T.W., Evans P.G., (1985), “Effect of Oil Rheology on Journal Bearing Performance: Part Instrumentation of the Big-End Bearing of a Fired Engine”, Proc Of the JSLE International Tribology Conference, 8-10 juillet, Tokyo, Japon, 1985 [58] Bates s T.W., Benwell S., Evans P.G.(1987), “Effect of Oil Rheology on Journal Bearing Performance : Part - Oil Film Thickness in the Big-End Bearing of an Operating Engine”, Proc 4th SAE Int Pacific Conference on Automotive Engineering, Melbourne, Australia, Paper No 871272 [59] Bates T.W., Benwell S.(1988),“Effect of Oil Rheology on Journal Bearing Performance : Part - Newtonian Oils in the Connecting -Rod Bearing of an Operating Engine”, SAE Paper No 880679 [60] Dutfoy L., Ligier J.-L.(2011), "Wear prediction of connecting rod bearing " 102 (Prédiction de l usure d un coussinet de tête de bielle), 10th EDF/Pprime Workshop, Futuroscope, October & [61] Moreau H.(2001), “Mesures des Epaisseurs du Film d’Huile dans les Paliers de Moteur Automobile et Comparaisons avec les Résultats Théoriques”, Thèse de Doctorat de Université de Poitiers [62] Ligier J.-L., Dutfoy L.(2010), "Modeling and prediction of a simplified seizure mechanism occurring in conrod bearing " (Modélisation et prévision d un mécanisme simplifié de grippage de palier de tête de bielle), 9th EDF/Pprime Workshop, Futuroscope, October & [63] Ligier J.-L., Dutfoy L.(2011), "Fatigue prediction for engine bearing" (Prévision de la tenue mécanique d un coussinet), 10th EDF/Pprime Workshop, Futuroscope, October & [64] Nguyễn Phương (1990), "Nghiên cứu thiết kế, tính tốn thực nghiệm ổ chặn", Luận án phó tiến sĩ, trường ĐH Bách khoa Hà Nội [65] Trần Thị Thanh Hải (1999), "Nghiên cứu tính toán ổ đỡ thuỷ động xây dựng thiết bị khảo sát đặc tính bơi trơn ổ", luận văn thạc sĩ, trường ĐH Bách khoa Hà Nội [66] Phan Thạch Hổ, Nguyễn Anh Tuấn, Phạm Văn Hùng (2001), " Một vài phương pháp tính tốn ổ trục chịu tải trọng động", Hội nghị khoa học lần thứ 19, Đại học Bách Khoa Hà Nội [67] Nguyễn Xuân Toàn, Trần Thị Thanh Hải, Dương Minh Tuấn (2001), "Thiết lập chương trình tính tốn bơi trơn thủy động (ổ đỡ, ổ chặn) có tính đến sai số hình học, hình dạng", Hội nghị khoa học lần thứ 19, Đại học Bách Khoa Hà Nội [68] Phan Thạch Hổ (2002), "Nghiên cứu thiết kế, tính tốn ổ trục khuỷu động NT855", Luận án tiến sĩ, trường ĐH Bách khoa Hà Nội [69] PGS.TS Phạm Văn Hùng (2005), "Nghiên cứu thiết kế chế tạo thiết bị đo áp suất ổ thủy động dùng phương pháp biểu thị đánh giá đại " Đề tài cấp giáo dục & đào tạo [70] TS.Trần Thị Thanh Hải (2014), "Mô hình hóa bơi trơn ổ đầu to truyền động nhiệt" Đề tài cấp trường, Trường ĐH Bách Khoa Hà Nội [71] TS Trần Thị Thanh Hải cộng (2016), "Nghiên cứu thiết kế chế tạo thiết bị thực nghiệm bôi trơn ổ đầu to truyền" Đề tài cấp bộ, Bộ Khoa học công nghệ [72] Nguyễn Đình Tân (2018) “Nghiên cứu mơ số bôi trơn thủy động ổ đầu to 103 truyền động đốt trong”, Luận án tiến sĩ, trường ĐH Bách khoa Hà Nội [73] Holmberga K., Andersson A P., Erdemirb A (2012), "Global energy consumption due to friction in passenger cars", Tribology International, vol 47, p 221-234 [74] Kragelsky, Alisin August (1968), "Cativation in lubricating films", Tribology, Volume 1, Issue 3, , Page 153-156 [75] Nguyễn Xuân Toàn (2007), Công nghệ bôi trơn NXB Bách khoa – Hà Nội, trang 180-182 [76] ThS Nguyễn Đình Tân, Lưu Trọng Thuận, TS Trần Thị Thanh Hải (2014), ‘’Xây dựng chương trình tự động tính tốn bơi trơn ổ đỡ có tính đến tượng gián đoạn màng dầu’’ Tạp chí khí Việt Nam, số năm 2014, Trang 54 – 60 [77] Đinh Tan Nguyen, Thi Thanh Hai Tran, Aurelian FATU (2016) ’’Effect of the Bolt Preload on Elastohydrodynamic Lubrication Analysis of a Connecting Rod Big – End Bering’’ Kỷ yếu hội nghị khoa học cơng nghệ tồn quốc khí – động lực 2016 tập 2, CK6 305, Trang 462- 467 [78] Nguyễn Đình Tân, Trần Thị Thanh Hải, Lưu Trọng Thuận (2018), ‘’Ảnh hưởng khe hở bán kính tới quỹ đạo tâm trục ổ đầu to truyền động 5S-FE.’’ Tạp chí Khoa học & Công nghệ trường đại học kỹ thuật [79] Trần Thị Thanh Hải, Nguyễn Đình Tân, Lưu Trọng Thuận (2018), ‘’Ảnh hưởng khe hở bán kính tới phân bố áp suất ổ đầu to truyền động 5S-FE’’ Tạp chí Khoa học & Cơng nghệ trường đại học kỹ thuật [80] Michael M Khonsari, E Richard Booserauth.(2008), M J Neale, M Priest, G Stachowiakeds "Applied Tribology Bearing Design and Lubrication", Second Edition [81] T.zeghloul, T.T.H.Tran, D.Bonneau (2002), "Modelisation numerique des interaction entre les differents coprs d’un palier de tete de bielle", Journal of Tribology [82] G Knoll, J Lang, A Rieacker (1996), "Transient EHD connecting rod analysis: Full dynamic versus quái-static deformation", Journal of Tribology, Vol 118, pp 349355 [83] T.T.H.Tran , T.Zeghloul, D.Bonneau (2007), "Experimental study of the interaction between the diferent bodies of a connecting rod big end bearing ", 12 th IFToMM world congress, Besancon, June 18-21, 2007 [84] T.T.H.Tran (2018), "A solution for creating the simulating load on connecting-rod in the experimental device for lubrication condiction of the connecting-rod big end bearing ", journal of science and technology, september, 2018 [85]Trần Thị Thanh Hải Lưu Trọng Thuận (2016)," Thermal Effect for Hydrodynamic Lubrication", Hội nghị KH&CN toàn quốc khí động lực 104 PHỤ LỤC Phụ lục 1: Quy trình đo tải tác dụng lên truyền Hình ảnh minh họa Bước thao tác Nội dung Bước Cấp nguồn cho module NI-9219 máy tính Kiểm tra kết nối thiết bị - module – máy tính, xác nhận kết nối qua cổng USB 2.0 Bước Khởi động chương trình lập trình đo tải phần mềm LabView Bước Cấu hình thiết bị đo, số mẫu, thời gian lấy mẫu, calib thiết bị đo Bước Điều chỉnh thông số hoạt động thiết bị: + Chế độ chạy (vô cấp/cố định) + Tốc độ quay trục khủy + 105 Bước Đo liệu tải + Trên giao diện LabView, click biểu tướng Start + Quan sát đồ thị hiển thị kết thông số tải Bước Kết thúc trình đo + Nhấn nút Stop giao diện LabView + Điều chỉnh tốc độ động “0” Tắt nguồn điện cho thiết bị Bước Lưu liệu Click chuột phải vào vị trí đồ thị, chọn Export => Data to Excel để lưu kết đo dạng file Excel xử lý số liệu 106 Phụ lục 2: Quy trình đo áp suất màng dầu ổ đầu to truyền Bước thao tác Nội dung Hình ảnh minh họa Bước Khởi động phát tín hiệu + Bật nút nguồn phát + Quan sát đèn tín hiệu (màu đỏ) phát, xác nhận thiết bị cấp đủ nguồn, hoạt động bình thường Bước Kết nối thu tín hiệu với máy tính qua cổng USB + Cắm cáp kết nối + Quan sát đèn tín hiệu (màu đỏ) thu sáng, xác nhận thiết bị hoạt động bình thường Bước Kiểm tra kết nối thu máy tính Bước Khởi động chương trình lập trình đo áp suất LabView Cài đặt thông số đầu vào (địa cổng COM) cho thu 107 Bước Điều chỉnh thông số hoạt động thiết bị: + Chế độ chạy (vô cấp/cố định) + Tốc độ quay trục khủy + Bước Đo liệu áp suất + Trên giao diện LabView, click biểu tượng Start + Nhấn chuột trái vào nút Read để bắt đầu nhận giá trị áp suất mà thu thu Bước Kết thúc trình đo + Nhấn nút Stop giao diện LabView + Điều chỉnh tốc độ động “0” Tắt nguồn điện cho thiết bị 108 Phụ lục 3: Quy trình đo nhiệt độ màng dầu ổ đầu to truyền Bước thao tác Nội dung Hình ảnh minh họa Cấp nguồn cho module NI-9213 máy tính Bước Kiểm tra kết nối thiết bị - module – máy tính, xác nhận kết nối qua cổng USB 2.0 Bước Khởi động chương trình lập trình đo tải LabView Cấu hình thiết bị đo, điều chỉnh Bước thông số thời gian lấy mẫu, số mẫu calib thiết bị đo Điều chỉnh thông số hoạt động thiết bị: Bước + Chế độ chạy (vơ cấp/cố đính) + Tốc độ quay trục khủy + 109 Thực đo nhiệt độ + Trên giao diện LabView, click Bước biểu tượng Start + Quan sát đồ thị hiển thị kết thơng số nhiệt độ Kết thúc q trình đo + Nhấn nút Stop giao diện Bước LabView + Điều chỉnh tốc độ động “0” Tắt nguồn điện cho thiết bị Lưu liệu Click chuột phải vào vị trí đồ thị, Bước chọn Export => Data to Excel để lưu kết đo dạng file Excel xử lý số liệu 110 Phụ lục 4: Code chương trình tính tốn mơ số PROGRAM THERMO HYDRODINAMIC LUBRICATION CONECTING ROD BIG END BEARING DIMENSION X(200),Y(200),P(20000),H(20000),T(20000),EDA(20000),POLD(20000) DATA EDA0,RO,C,AJ,B,R,RATIO,AN,EPSON/0.05,890.0,1870.0,4.184,1.0E-1,6.0E2,0.003,3.0E3,0.7/ OPEN(70,FILE='FILM.DAT',STATUS='UNKNOWN') OPEN(8,FILE='PRESSURE.DAT',STATUS='UNKNOWN') OPEN(9,FILE='TEM.DAT',STATUS='UNKNOWN') OPEN(10,FILE='EDA.DAT',STATUS='UNKNOWN') PI=3.1415926 N=61 M=41 DX=2.0*PI/FLOAT(N-1) DY=1./FLOAT(M-1) ALFA1=R/B C0=RATIO*R U=R*AN*2.0*PI/60.0 ALFA=(R/B*DX/DY)**2 A=U*R*EDA0/(2.0*AJ*RO*C*C0**2) T0=303.0 KG=0 CALL INIT(N,M,DX,DY,A,T0,EPSON,X,Y,H,T,EDA) 10 CALL SUBP(N,M,DX,EPSON,ALFA,H,P,EDA) CALL THERM(N,M,A,ALFA1,DX,DY,EDA0,T0,X,Y,P,H,T,EDA) IF(KG.EQ.0)THEN CALL PSAVE(KG,N,M,P,POLD,ERO) KG=1 GOTO 10 ENDIF CALL PSAVE(KG,N,M,P,POLD,ERO) WRITE(*,*)'EROP,KG=',ERO,KG 111 KG=KG+1 IF(ERO.GT.1.E-6.AND.KG.LT.10)GOTO 10 CALL OUTPUT(N,M,A,T0,X,Y,H,P,T,EDA) STOP END SUBROUTINE INIT(N,M,DX,DY,A,T0,EPSON,X,Y,H,T,EDA) DIMENSION X(N),Y(M),H(N,M),T(N,M),EDA(N,M) DO I=1,N X(I)=(I-1.0)*DX ENDDO DO J=1,M Y(J)=-0.5+(J-1)*DY DO I=1,N H(I,J)=1.0+EPSON*COS(X(I)) T(I,J)=T0/A EDA(I,J)=1.0 ENDDO ENDDO RETURN END SUBROUTINE SUBP(N,M,DX,EPSON,ALFA,H,P,EDA) DIMENSION H(N,M),P(N,M),EDA(N,M) DO I=2,N-1 DO J=2,M-1 P(I,J)=0.5 ENDDO ENDDO DO J=1,M P(1,J)=0.0 P(N,J)=0.0 ENDDO DO I=1,N 112 P(I,1)=0.0 P(I,M)=0.0 ENDDO IK=0 10 C1=0.0 ALOAD=0.0 DO I=2,N-1 I1=I-1 I2=I+1 DO J=2,M-1 PD=P(I,J) J1=J-1 J2=J+1 A1=(0.5*(H(I2,J)+H(I,J)))**3/(0.5*(EDA(I2,J)+EDA(I,J))) A2=(0.5*(H(I,J)+H(I1,J)))**3/(0.5*(EDA(I,J)+EDA(I1,J))) A3=ALFA*H(I,J)**3/(0.5*(EDA(I,J2)+EDA(I,J))) A4=ALFA*H(I,J)**3/(0.5*(EDA(I,J)+EDA(I,J1))) A5=A1*P(I2,J)+A2*P(I1,J)+A3*P(I,J2)+A4*P(I,J1) A6=A1+A2+A3+A4 P(I,J)=(-DX*(H(I2,J)-H(I1,J))+A5)/A6 P(I,J)=0.7*PD+0.3*P(I,J) IF(P(I,J).LT.0.0)P(I,J)=0.0 C1=C1+ABS(P(I,J)-PD) ALOAD=ALOAD+P(I,J) ENDDO ENDDO IK=IK+1 C1=C1/ALOAD IF(C1.GT.1.E-6)GOTO 10 RETURN END SUBROUTINE PSAVE(KG,N,M,P,POLD,ERO) 113 DIMENSION P(N,M),POLD(N,M) IF(KG.EQ.0)GOTO 10 ERO=0.0 EROMAX=-1.0 W=0.0 DO I=1,N DO J=1,M AE=ABS(P(I,J)-POLD(I,J)) ERO=ERO+AE IF(AE.GT.EROMAX)THEN II=I JJ=J EROMAX=AE ENDIF W=W+ABS(P(I,J)) ENDDO ENDDO ERO=ERO/W 10 DO I=1,N DO J=1,M POLD(I,J)=P(I,J) ENDDO ENDDO RETURN END SUBROUTINE THERM(N,M,A,ALFA1,DX,DY,EDA0,T0,X,Y,P,H,T,EDA) DIMENSION X(N),Y(M),H(N,M),P(N,M),T(N,M),EDA(N,M) KG=0 10 ERT=0.0 DO I=2,N DO J=1,M/2+1 TOLD=T(I,J) 114 EDA(I,J)=EXP((ALOG(EDA0)+9.67)*(((A*T(I,J)-138.0)/(T0-138.0))**(-1.1)-1.0)) IF(I.NE.N)DPDX=0.5*(P(I+1,J)-P(I-1,J))/DX IF(I.EQ.N)DPDX=0.5*(P(1,J)-P(I-1,J))/DX QX=0.5*H(I,J)-0.5*H(I,J)**3*DPDX DPDY=ALFA1*(P(I,J+1)-P(I,J))/DY IF(J.EQ.M/2+1)DPDY=0.0 QY=-0.5*H(I,J)**3*DPDY AA=QX/DX*T(I-1,J)-ALFA1*QY*T(I,J+1)/DY AB=2.0*EDA(I,J)/H(I,J) AC=6.0*H(I,J)/EDA(I,J)*(DPDX**2+ALFA1**2*DPDY**2) BB=QX/DX-ALFA1*QY/DY T(I,J)=(AA+AB+AC)/BB IF(A*T(I,J).GE.403.)THEN WRITE(*,*)'T OVER THE LIMIT 100' WRITE(*,*)I,J,T(I,J) STOP 00001 ENDIF T(I,J)=0.7*TOLD+0.3*T(I,J) ERT=ERT+ABS(T(I,J)-TOLD)/303 ENDDO ENDDO ERT=A*ERT/((N-1)*(M-1)) KG=KG+1 DO J=1,M T(1,J)=0.5*(T(1,J)+T(N,J)) ENDDO WRITE(*,*)ERT,'KG=',ERT,KG IF(ERT.GT.1.E-3)GOTO 10 DO I=1,N DO J=1,M/2 T(I,M-J+1)=T(I,J) EDA(I,M-J+1)=EDA(I,J) 115 ENDDO ENDDO RETURN END SUBROUTINE OUTPUT(N,M,A,T0,X,Y,H,P,T,EDA) DIMENSION X(N),Y(M),H(N,M),P(N,M),T(N,M),EDA(N,M) TMAX=0.0 DO I=1,N DO J=1,M T(I,J)=(A*T(I,J)-T0) IF(T(I,J).GT.TMAX)TMAX=T(I,J) END DO ENDDO WRITE(*,*)'TAMX=',TMAX i=1,M j=1,N WRITE(9,30)X(j),Y(i),T(j,i) enddo !$$$$$$ WRITE(70,30)X(1),(Y(J),J=1,M) !$$$$$$ WRITE(8,30)X(1),(Y(J),J=1,M) !$$$$$$ WRITE(9,30)X(1),(Y(J),J=1,M) !$$$$$$ WRITE(10,30)X(1),(Y(J),J=1,M) !$$$$$$ DO I=1,N !$$$$$$ WRITE(70,30)X(I),(H(I,J),J=1,M) !$$$$$$ WRITE(8,30)X(I),(P(I,J),J=1,M) !$$$$$$ WRITE(9,30)X(I),(T(I,J),J=1,M) !$$$$$$ WRITE(10,30)X(I),(EDA(I,J),J=1,M) ENDDO 30 FORMAT(3F10.5) RETURN END 116 ... giới nghiên cứu nhiều vấn đề bôi trơn ổ đầu to truyền Tuy nhiên Vi ệt Nam chưa có nhiều cơng trình nghiên cứu bơi trơn ổ đầu to truyền, nghiên cứu bổ sung thêm cho nghiên cứu bôi trơn ổ đầu to. .. nhiều nghiên c ứu bôi trơn ổ đầu to truyền mô số thực nghiệm Hơn nữa, nghiên cứu chưa tính tới hiệu ứng nhiệt ổ Vì mà em lựa chọn đề tài ‘? ?Nghiên cứu đặc tính bơi trơn nhiệt thủy động ổ có dạng đầu. .. hình nghiên cứu giới Từ đầu năm 1969 đến có nhiều nghiên cứu lý thuyết thực nghiệm nghiên cứu bôi trơn cho ổ đầu to truyền Các nghiên cứu bôi trơn ổ đầu to truyền theo hai hướng nghiên cứu lý