(Luận văn thạc sĩ hcmute) nghiên cứu và xây dựng giải pháp nâng cao ổn định hệ thống treo ôtô

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN NGỌC CHÂU NGHIÊN CỨU VÀ XÂY DỰNG GIẢI PHÁP NÂNG CAO ỔN ĐỊNH HỆ THỐNG TREO Ô TƠ NGÀNH: KỸ THUẬT CƠ KHÍ ĐỘNG LỰC – 60520116 S K C0 3 Tp Hồ Chí Minh, tháng 10/2019 Luan van BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN NGỌC CHÂU NGHIÊN CỨU VÀ XÂY DỰNG GIẢI PHÁP NÂNG CAO ỔN ĐỊNH HỆ THỐNG TREO Ô TÔ NGÀNH: KỸ THUẬT CƠ KHÍ ĐỘNG LỰC – 60520116 Hướng dẫn khoa học: TS NGUYỄN VĂN TRẠNG Tp Hồ Chí Minh, tháng 10 năm 2019 Luan van Luan van Luan van Luan van Luan van Luan van Luan van Luan van Luan van TÀI LIỆU THAM KHẢO [1] Lương Hữu Thoại, nghiên cứu biện pháp tối ưu hóa hệ thống treo buýt sử dụng Tp Hồ Chí Minh, năm 2015 [2] Dương Nguyễn Hắc Lân, nghiên cứu mơ hình hóa mơ hệ thống treo chủ động Matlab/simulink, Đại học sư phạm Tp.HCM năm 2019 [3] Nguyễn Đức Ngọc, Deng Zhaoxiang, phương pháp điều khiển tối ưu hệ thống treo chủ động ô tô, The College of Mechanical Engineering of Chongqing University, China 400030 [4] Hồ Văn Hóa, nghiên cứu mức độ êm dịu xe bus DAEWOO BC 212 MA sử dụng Tp Hồ Chí Minh, năm 2015 [5] Lê Thiện Kim Ngân, nghiên cứu ảnh hưởng tải trọng động đến an toàn chuyển động loại xe khách chuyển động đường, năm 2015 [6] XueMei Sun, Yaxu Chu, Jiuchen Fan, Qiuxiao Yang, Research of Simulation on the Effect of Suspension Damping on Vehicle Ride, International Conference on Future Electrical Power and Energy Systems,2012 [7] John E D Ekoru, Olurotiri A Dahunsi, Jiroh O Pedro PID Control of a Nonlinear Half-Car Active Suspension System via Force Feedback.School of Mechanical, Aeronautical and Industial Engineering, University of the Witwatersrand, Private Bag 03, WITS2050, Johannesburg, South Africa, September 2011 [8] Dirman Hanafi, PID Controller Design for Semi-Active Car Suspension Based on Model from Intelligent System Identification.Department of Mechatronic and Robotic Engineering Universiti Tun Hussein On Malaysia (UTHM),2010 [9] Alexey Kuznetsov, Musa Mammadov, Ibrahim Sultan, Eldar Hajilarov Optimization of a quarter-car suspension model coupled with the driver biomechanical effects,Journal of Sound and Vibration 330,2011 [10] Dr Hayder Sabah Abd AL-Amir, Assistant Lecturer Ali Talib Abd Al Zahra Design a Robust PID Controller of an Active Suspension System, International 94 Luan van Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:14 No:02 april, 2014 [11] Anand Tandel ,A R Deshpande, S P Deshmukh, K R Jagtap, Modeling, Analysis and PID Controller Implementation on Double Wishbone Suspension Using SimMechanics and Simulink, 12th global congress on manufacturing and management, gcmm ,2014 [12] Musa Mohammed Bello, Amir Akramin Shafie ,Raisuddin Md Khan, ElectroHydraulic PID Force Control for Nonlinear Vehicle Suspension System, International Journal of Engineering Research & Technology (IJERT),ISSN : 2278-0181,Vol Issue 01,January-2015 [13] Abdullahi b kunya, atef a ata, half car suspension system integrated with pid controller, Mevlana University, Konya, , Konya –Turkey, april, 2003 [14] Mohammadjavad Zeinali, Intan Zaurah Mat Darus,Fuzzy PID Controller Simulationfor a Quarter-car Semi-active Suspension System Using Magnetorheological Damper, 2012 IEEE Conference on Control, Systems and Industrial Informatics (ICCSII) Bandung, Indonesia, September 23-26, 2012 [15] PGS TS Nơng Văn Vìn, Động Lực Học Thẳng Đứng Và Hệ Thống Treo ô tô, Trường Sư Phạm Kỹ Thuật Hưng Yên, năm 2014 [16] TS Lâm Mai Long, Dao Động Tiếng Ồn Ơ Tơ, NXB.ĐHQG TP.HCM – 2017 [17] Tiêu chuẩn Việt Nam – TCVN 6964 – 1:2001 ; ISO 2631 – 1: 1997 [18] Ikbal Eski, S _ ahin Yıldırım, Vibration control of vehicle active suspension system using a new robust neural network control system,Erciyes University, Faculty of Engineering, Mechanical Engineering Department, Talas Cad, Kayseri 38039, Turkey,2009 [19] Lalitkumar Maikulal Jugulkar, Shankar Singh, Suresh Maruti Sawant Analysis of suspension with variable stiffness and variable damping force for automotive applications, Advances in Mechanical Engineering, 2016 95 Luan van [20] Bùi Quốc Vĩnh, nghiên cứu mức độ êm dịu xe ô tô khách 29 chỗ ngồi sản xuất Việt Nam, Tp Đà Nẵng, năm 2011 [21] Trương Hoàng Tuấn, Xác định tải trọng động tác dụng lên chassi xe sym t880 phương pháp mô phỏng, Trường đại học bách khoa Tp.hcm, năm 2013 [22] VDHS-10 Vehicle Suspension [23] Nguyễn Đức Luyện, Lý thuyết ô tô quân Nhà xuất Quân Đội Nhân Dân Hà Nội, năm 2002 [24] Nguyễn Hữu Thảo, Mô khảo sát dao động ô tô vận tải hàng khách Matlab – Simulink, Đà Nẵng, năm 2012 [25] Nguyễn Phùng Quang, Matlab – Simulink dành cho kỹ sư điều khiển tự động, NXB Khoa Học Kỹ Thuật, Hà Nội 2004 [26] Nguyễn Dỗn Phước, Lý thuyết điều khiển tuyến tính, NXB Khoa Học Kỹ Thuật, Hà Nội 2009 [27] Team 21 Toyota, tài liệu đào tạo kỹ thuật viên http://thacothanhtrung.vn/san-pham/bus/thaco-county/0/86/thaco-hb73s-h140.aspx 96 Luan van PHỤ LỤC Code matlab: %% thong so xe bus p1=50000; i1=400; d1=50000; p=40000; i=200; d=40000; ms=5867;%kg msf=2607;%kg msr=3260;%kg mu=783;% kg muf=348;%kg mur=435;%kg Lf=2.27;% m Lr=1.815;%m B=1.705/2;%m Ix=4264; %kgm2 Iy=19337; % kgm2 klf=596467;% N/m klr=1192935;%N/m ksf=71890;% N/m ksr=130667;% N/m csr= 8256; % Ns/m csf=5476;% Ns/m %% Ham mat duong 97 Luan van h=0.03; % chieu cao map mo (m) v=20; % toc (m/s) lamda=1; % (m) l=4.085;%(m) t=0; phi=0; T=t+10;%(m) omega=2*pi*v/lamda;%(rad/s) qf=h/2*omega*t;% ham tac dung len banh truoc qr=h/2*(omega*t+2*pi*l/lamda); % ham tac dung len banh sau dem=0; k=1000; ghp=40000; ghi=200; ghd=40000; for p=0:k:ghp; for i=0:k/10:ghi; for d=0:k:ghd; dem=dem+1; sim('sosanhketqua_bd_cd'); atb(dem)=sqrt(mean(gtzcga.*gtzcga)); atb_p(dem)=sqrt(mean(gtzcgp.*gtzcgp)); sumax(dem)=max(abs(bdnhipf_a)); bdnhiptruoc_a(dem)=sqrt(mean(bdnhipf_a.*bdnhipf_a)); bdnhiptruoc_p(dem)=sqrt(mean(bdnhipf_p.*bdnhipf_p)); gtgocphitb_a(dem)=sqrt(mean(gtgocphi_a.*gtgocphi_a)); gtgocphitb_p(dem)=sqrt(mean(gtgocphi_p.*gtgocphi_p)); 98 Luan van bdloptruoc_a(dem)=sqrt(mean(bdlopf_a.*bdlopf_a)); bdloptruoc_p(dem)=sqrt(mean(bdlopf_p.*bdlopf_p)); luctruoc(dem=sqrt(mean(luc_ff.*luc_ff)); mta(dem,1)=p; mta(dem,2)=i; mta(dem,3)=d; mta(dem,4)=atb(dem); mta(dem,5)=sumax(dem); mta(dem,6)=gtgocphitb_a(dem); mta(dem,7)=bdloptruoc_a(dem); end end end a=dem; mina=1; mina1=1; for ii=1:1:dem; if (atb(ii)