(Luận văn thạc sĩ) nghiên cứu đánh giá độ êm dịu của xe thiết giáp BTR 60PB bằng matlab simulink

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Ĩ NGHIÊM XUÂN THÀNH NGHIÊN CỨU ÐÁNH GIÁ ÐỘ ÊM DỊU CỦA XE THIẾT GIÁP BTR – 60PB BẰNG MATLAB – SIMULINK NGÀNH: KỸ THUẬT CƠ KHÍ ÐỘNG LỰC - 60520116 SKC004972 Tp Hồ Chí Minh, tháng 10/2016 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Ỹ NGHIÊM XUÂN THÀNH NGHIÊN CỨU ĐÁNH GIÁ ĐỘ ÊM DỊU CỦA XE THIẾT GIÁP BTR – 60PB BẰNG MATLAB – SIMULINK NGÀNH: KỸ THUẬT CƠ KHÍ ĐỘNG LỰC - 60520116 Tp Hồ Chí Minh, tháng 10/ năm 2016 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Ỹ NGHIÊM XUÂN THÀNH NGHIÊN CỨU ĐÁNH GIÁ ĐỘ ÊM DỊU CỦA XE THIẾT GIÁP BTR – 60PB BẰNG MATLAB – SIMULINK NGÀNH: KỸ THUẬT CƠ KHÍ ĐỘNG LỰC - 60520116 Hướng dẫn khoa học: GVC TS NGUYỄN CHÍ THANH Tp Hồ Chí Minh, tháng 10/ năm 2016 LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƯỢC: Họ tên: Nghiêm Xuân Thành Giới tính: Nam Ngày, tháng, năm sinh: 21/01/1971 Nơi sinh: Vĩnh Phúc Quê quán: Quang Sơn, Lập Thạch, Vĩnh Phúc Dân tộc: Kinh Chỗ riêng địa liên lạc: 116/29 đường Nguyễn Văn Lượng, Phường 17, quận Gị Vấp, thành phố Hồ Chí Minh Điện thoại quan: Điện thoại nhà riêng: 0983414279 Fax: E-mail: nghiemxuanthanhsqktqs@gmail.com II QUÁ TRÌNH ĐÀO TẠO: Trung học chuyên nghiệp: Hệ đào tạo: Chính quy Thời gian đào tạo từ 3/1990 đến 12/1990 Nơi học (trường, thành phố): Ngành học: Trường Hạ sĩ quan Xe tăng II (Vĩnh Phúc) Kỹ thuật Tăng thiết giáp Đại học: * Hệ đào tạo: Cao đẳng quy Thời gian đào tạo từ 8/1992 đến 8/1996 Nơi học (trường, thành phố): Trường Sĩ quan Kỹ thuật Quân sự, Tp HCM Ngành học: Cơ khí động lực Tên đồ án: Khảo sát hệ thống trợ lực lái xe tăng T55 Ngày & nơi bảo vệ đồ án tốt nghiệp: 6/1996 Trường Sĩ quan Kỹ thuật Quân Người hướng dẫn: TS Phạm Văn Trợ * Hệ đào tạo: Đại học chuyên tu Thời gian đào tạo từ 8/1997 đến 5/2000 Nơi học (trường, thành phố): Học viện Kỹ thuật Quân sự, Hà Nội Ngành học: Cơ khí động lực Tên mơn thi tốt nghiệp: Triết học, Kỹ thuật CN tổng hợp Ngày & nơi thi tốt nghiệp: 3/2000, Học viện Kỹ thuật Quân i III Q TRÌNH CƠNG TÁC CHUN MƠN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Công việc đảm nhiệm Cán quản lý giáo dục (có tham gia thỉnh giảng) Cán giảng viên ii LỜI CAM ĐOAN Tôi xin cam đoan luận văn cơng trình nghiên cứu riêng tôi, thực với hướng dẫn khoa học Tiến sĩ Nguyễn Chí Thanh, số liệu kết nghiên cứu nêu luận văn tốt nghiệp hoàn toàn trung thực chưa công bố cơng trình khác Nội dung luận văn hồn toàn phù hợp với tên đề tài đăng ký phê duyệt Hiệu trưởng Trường Đại học Sư phạm Kỹ thuật Tp Hồ Chí Minh Tp Hồ Chí Minh, ngày 25 tháng 10 năm 2016 Tác giả luận văn Nghiêm Xuân Thành iii LỜI CẢM TẠ Trong suốt q trình học tập thực hồn thành luận văn, nhận giúp đỡ quý báu thầy cô bạn bè đồng nghiệp Với kính trọng biết ơn sâu sắc, tơi xin bày tỏ lòng cảm tạ chân thành tới: Ban Giám hiệu, Phòng Đào tạo, Bộ phận sau đại học, Khoa Cơ Khí Động Lực Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, giúp đỡ tạo điều kiện cho suốt trình học tập hồn thành luận văn Giảng viên chính, Tiến sĩ Nguyễn Chí Thanh, Trường Sĩ quan Kỹ thuật Quân Sự tận tình hướng dẫn, động viên giúp đỡ tạo điều kiện cho trình cơng tác, học tập hồn thành luận văn tốt nghiệp Xin chân thành cảm ơn quý thầy Hội đồng chấm luận văn tốt nghiệp (khóa học 2015 - 2017) cho tơi đóng góp q báu để hoàn chỉnh luận văn tiến độ kế hoạch, bảo đảm chất lượng hiệu trình nghiên cứu thực Xin cảm ơn bạn bè đồng nghiệp, gia đình nhiệt tình cổ vũ, giúp đỡ tơi q trình hồn thành luận văn Xin kính chúc Ban Giám hiệu, quý thầy cô mạnh khỏe, hạnh phúc thành công nghiệp trồng người, ánh sáng đường khoa học đầy gian chuân thử thách, dẫn đường cho hệ hôm mai sau vững bước vươn xa, “như sóng dịu êm vỗ bờ, hàng ngày hàng giờ, xô đến bến bờ nhân cách mới” Xin chân thành cảm ơn! iv TÓM TẮT Đề tài: "Nghiên cứu đánh giá độ êm dịu xe thiết giáp BTR - 60PB Matlab - Simulink", nội dung gồm: Thứ nhất, nghiên cứu tổng quan mô tả phương pháp khảo sát dao động hệ thống treo xe bốn cầu phần mềm Matlab - Simulink Thứ hai, xây dựng mơ hình dao động thiết lập hệ phương trình vi phân dao động Thứ ba, xác định giá trị thông số hệ thống dao động Khảo sát đánh giá chất lượng hệ thống treo với ưu điểm nhược điểm Thứ bốn, tối ưu hóa thiết kế hệ thống treo ghế lái theo phương pháp thuật toán di truyền, nhằm nâng cao độ êm dịu cho vị trí ghế lái Phần cuối, kết luận khuyến nghị hướng nghiên cứu Ưu điểm luận văn vừa có tính lý thuyết tính thực tiễn Đề tài quan tâm đến vấn đề dao động ghế lái xe hệ thống tổng thể thân xe, đưa kết hoàn toàn độ tin cậy Chúng ta phát triển đề tài để cải thiện nâng cấp ghế lái hệ thống treo xe bọc thép BTR - 60PB SUMMARY Thesis: "Study assessing the mellow of armored vehicles BTR - 60PB by Matlab - Simulink", the contents are: Firstly, the research described the overview of survey methodology on oscillations of the suspension in four different boxes vehicle by Matlab - Simulink software Secondly, building vehicle model and establishing oscillation equations Next, determining the value of the parameters of the oscillating system Surveying and assessing the quality of the suspension with the advantages and disadvantages After that, optimal seat suspension design of vehicle using genetic algorithms, to enhance the mellow to the driver's seat position Finally, conclusions and recommendations for further studies The advantage of this topic was both theoretical and practical This thesis was concerned about the problems in the driving seat oscillations in overall body system, has launched a completely new results and reliability We can develop this theme to improve the driver's seat and upgrad suspensional system of the armored vehicles BTR - 60PB v MỤC LỤC Tựa trang TRANG Quyết định giao đề tài Lý lịch cá nhân i Lời cam đoan iii Lời cảm tạ iv Tóm tắt v Mục lục vi Danh sách ký hiệu x Danh sách hình vẽ xii Danh sách bảng xiv Chương 1: TỔNG QUAN NGHIÊN CỨU 1.1 Tổng quan hướng nghiên cứu 1.1.1 Tình hình nghiên cứu đề tài 1.1.1.1 Ngoài nước 1.1.1.2 Trong nước 1.1.2 Các hướng nghiên cứu dao động ô tô 1.1.3 Dao động ô tô ảnh hưởng dao động 1.1.3.1 Dao động ô tô 1.1.3.2 Ảnh hưởng dao động ô tô 1.1.3.3 Các yếu tố ảnh hưởng tới dao động ô tô 11 1.1.4 Các tiêu đánh giá độ êm dịu ô tô 12 1.1.4.1 Tần số dao động 13 1.1.4.2 Gia tốc dao động 13 1.1.4.3 Chỉ tiêu dựa số liệu cảm giác theo gia tốc vận tốc dao động 13 1.1.4.4 Đánh giá cảm giác theo công suất dao động 15 1.1.4.5 Đánh giá cảm giác theo gia tốc dao động thời gian tác động 16 1.1.4.6 Đánh giá dựa theo trị số hiệu gia tốc dao động (bình phương trung bình gia tốc dao động) 18 vi 1.1.4.7 Chỉ tiêu an toàn chuyển động tải trọng tác dụng xuống đường 19 1.1.4.8 Tiêu chuẩn độ êm dịu Việt Nam 21 1.1.5 Khái quát mô 22 1.1.5.1 Khái niệm mô 22 1.1.5.2 Khái quát mô dao động ô tô 23 Tính cấp thiết ý nghĩa đề tài 24 2.1 Tính cấp thiết đề tài 24 1.2.2 Ý nghĩa khoa học thực tiễn 25 1.3 Mục đích nghiên cứu đối tượng nghiên cứu 26 1.3.1 Mục đích nghiên cứu 26 1.3.2 Đối tượng nghiên cứu 26 1.4 Nhiệm vụ giới hạn đề tài 26 1.4.1 Nhiệm vụ đề tài 26 1.4.2 Giới hạn đề tài 26 1.5 Phương pháp nghiên cứu 27 Chương 2: CƠ SỞ LÝ THUYẾT THIẾT LẬP HỆ PHƯƠNG TRÌNH VI PHÂN DAO ĐỘNG CỦA XE 28 2.1 Các thành phần hệ dao động 28 2.1.1 Khối lượng treo M 28 2.1.2 Khối lượng không treo m 29 2.2 Mơ hình dao động tương đương xe cầu 29 2.2.1 Mơ hình phẳng dao động liên kết tổng qt “đường – xe – người” 30 2.2.2 Mơ hình dao động phần tư xe 31 2.2.3 Mơ hình phẳng dao động ngang thân xe 31 2.2.4 Mơ hình dao động khơng gian xe 32 2.2.5 Mơ hình phẳng dao động dọc thân xe 33 2.3 Hàm kích thích dao động từ mặt đường 34 2.4 Các giả thiết xây dựng mơ hình khảo sát 36 2.5 Thiết lập hệ phương trình vi phân mơ tả mơ hình dao động 38 vii RMSddZCK{iCb,iKb} =RMSddZ; iKb iCb if max(RMSddzse)>2.5 || max(Maxzsetd)>0.15 OBJ(iCb,iKb) = 50; else OBJ(iCb,iKb) = max(RMSddzse)+max(Maxzsetd); end end end %%%%%% CODE MATLAB GUIDE function buttonMophong_Callback(hObject, eventdata, handles) % hObject handle to buttonMophong (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) %str = get(handles.inputFuntion,'String'); %val = get(handles.inputFuntion,'Value'); %switch val % case % assignin('base','inputFuntion',2); % case % assignin('base','inputFuntion',3); % case % assignin('base','inputFuntion',1); % case % assignin('base','inputFuntion',4); %end msgbox('Bat dau mo phong','Thong bao','help'); tStart = str2num(get(handles.tStart,'String')); tStop = str2num(get(handles.tStop,'String')); temp = who; for iv = 1:length(temp) assignin('base',temp{iv},eval(temp{iv})); end who; l1 = evalin('base','l1'); l2 = evalin('base','l2'); l3 = evalin('base','l3'); l4 = evalin('base','l4'); lse = evalin('base','lse'); S = evalin('base','S'); w0 = evalin('base','w0'); Cb0 = evalin('base','Cb0'); Kb0 = evalin('base','Kb0'); for iCb = 1:length(Cb0) for iKb = 1:length(Kb0) cse = Cb0(iCb); kse = Kb0(iKb); assignin('base','cse',cse); assignin('base','kse',kse); for iw = 1:length(w0) w = w0(iw); assignin('base','w',w); if w == t1 = 0; 85 t2 = 0; t3 = 0; assignin('base','t1',t1); assignin('base','t2',t2); assignin('base','t3',t3); else end t1=2*pi()*(l1-l2)/(w*S); t2=2*pi()*(l1+l3)/(w*S); t3=2*pi()*(l1+l4)/(w*S); assignin('base','t1',t1); assignin('base','t2',t2); assignin('base','t3',t3); opt = simset('InitialStep',1e3,'OutputVariables','txy','MaxStep',0.01); [t,x,y] = sim('mohinh_510', [tStart, tStop], opt); RMSZ(iw) = max(SZ); RMSdZ(iw) = rms(SdZ); RMSddZ(iw) = rms(SddZ); RMSphi(iw) = max(Sphi); RMSdphi(iw) = rms(Sdphi); RMSddphi(iw) = rms(Sddphi); RMSzsetd(iw) = max(Szsetd); RMSzse(iw) = rms(Szse); RMSdzse(iw) = rms(Sdzse); RMSddzse(iw) = rms(Sddzse); z1 = SZ - l1*Sphi; dz1 = SdZ - evalin('base','lse')*Sdphi; ddz1 = SddZ - evalin('base','lse')*Sddphi; RMSz1(iw) = max(z1) RMSdz1(iw) = rms(dz1); RMSddz1(iw) = rms(ddz1); RMSFw1(iw) RMSFw2(iw) RMSFw3(iw) RMSFw4(iw) = = = = rms(Fw1); rms(Fw2); rms(Fw3); rms(Fw4); end %ket thuc iw RMSddZ1CK{iCb,iKb} =RMSddz1 ; RMSddzseCK{iCb,iKb} = RMSddzse; MaxzsetdCK{iCb,iKb} =RMSzsetd ; RMSddZCK{iCb,iKb} =RMSddZ; temp = who; for iv = 1:length(temp) assignin('base',temp{iv},eval(temp{iv})); end %ket thuc iKb end %ket thuc iCb end h = msgbox('Mo phong xong','Trang thai','help'); 86 % - Executes on button press in buttonUpdate function buttonUpdate_Callback(hObject, eventdata, handles) % hObject handle to buttonUpdate (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) assignin('base','S',str2num(get(handles.editS,'String'))); assignin('base','fks',str2num(get(handles.editVks,'String'))); assignin('base','qo',str2num(get(handles.editq0,'String'))); assignin('base','cse',str2num(get(handles.editcse,'String'))); assignin('base','kse',str2num(get(handles.editkse,'String'))); h = msgbox('Da cap nhat','Thong bao','help'); % - Executes on button press in buttonKq function buttonKq_Callback(hObject, eventdata, handles) % hObject % eventdata % handles str_1 = get(handles.typeData,'String'); val_1 = get(handles.typeData,'Value'); GRMSZ = evalin('base','RMSZ'); GRMSdZ = evalin('base','RMSdZ'); GRMSddZ = evalin('base','RMSddZ'); GRMSzse = evalin('base','RMSzsetd'); GRMSdzse = evalin('base','RMSdzse'); GRMSddzse = evalin('base','RMSddzse'); GRMSphi = evalin('base','RMSphi'); GRMSdphi = evalin('base','RMSdphi'); GRMSddphi = evalin('base','RMSddphi'); GRMSz1 = evalin('base','RMSz1'); GRMSdz1 = evalin('base','RMSdz1'); GRMSddz1 = evalin('base','RMSddz1'); GRMSFw1 GRMSFw2 GRMSFw3 GRMSFw4 = = = = evalin('base','RMSFw1'); evalin('base','RMSFw2'); evalin('base','RMSFw3'); evalin('base','RMSFw4'); switch val_1 case set(handles.editGiatoc,'String',num2str(max(GRMSZ))); case set(handles.editGiatoc,'String',num2str(max(GRMSdZ))); case set(handles.editGiatoc,'String',num2str(max(GRMSddZ))); case set(handles.editGiatoc,'String',num2str(max(GRMSphi))); case set(handles.editGiatoc,'String',num2str(max(GRMSdphi))); case set(handles.editGiatoc,'String',num2str(max(GRMSddphi))); case set(handles.editGiatoc,'String',num2str(max(GRMSz1))); 87 case case case 10 case 11 case 12 case 13 case 14 case 15 case 16 set(handles.editGiatoc,'String',num2str(max(GRMSdz1))); set(handles.editGiatoc,'String',num2str(max(GRMSddz1))); set(handles.editGiatoc,'String',num2str(max(GRMSzse))); set(handles.editGiatoc,'String',num2str(max(GRMSdzse))); set(handles.editGiatoc,'String',num2str(max(GRMSddzse))); set(handles.editGiatoc,'String',num2str(round(max(GRMSFw1)))); set(handles.editGiatoc,'String',num2str((round(max(GRMSFw2))))); set(handles.editGiatoc,'String',num2str((round(max(GRMSFw3))))); set(handles.editGiatoc,'String',num2str((round(max(GRMSFw1))))); otherwise end % - Executes on button press in buttonPlot function buttonPlot_Callback(hObject, eventdata, handles) % hObject handle to buttonPlot (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) %Locals funtion grid on; evalin('base','who'); str = get(handles.typeData,'String'); val = get(handles.typeData,'Value'); grid on; axes(handles.axes1); if (get(handles.checkboxHold,'Value')) hold on; else hold off; end str_1 = get(handles.popupmenuChonmau,'String'); val_1 = get(handles.popupmenuChonmau,'Value'); switch val_1 case color = 'r'; case color = 'm'; case color = 'c'; case color ='g'; case color = 'b'; case color = 'w'; case color = 'k'; end str_2 = get(handles.Doday,'String'); val_2 = get(handles.Doday,'Value'); 88 switch val_2 case Linewidth case Linewidth case Linewidth case Linewidth case Linewidth case Linewidth case Linewidth end = 1.5; = 2; = 2.5; = 3; = 3.5; = 4; = 4.5; f0 = evalin('base','f0'); GRMSZ = evalin('base','RMSZ'); GRMSdZ = evalin('base','RMSdZ'); GRMSddZ = evalin('base','RMSddZ'); GRMSzse = evalin('base','RMSzsetd'); GRMSdzse = evalin('base','RMSdzse'); GRMSddzse = evalin('base','RMSddzse'); GRMSphi = evalin('base','RMSphi'); GRMSdphi = evalin('base','RMSdphi'); GRMSddphi = evalin('base','RMSddphi'); GRMSz1 = evalin('base','RMSz1'); GRMSdz1 = evalin('base','RMSdz1'); GRMSddz1 = evalin('base','RMSddz1'); GRMSFw1 GRMSFw2 GRMSFw3 GRMSFw4 = = = = evalin('base','RMSFw1'); evalin('base','RMSFw2'); evalin('base','RMSFw3'); evalin('base','RMSFw4'); RMSddzseCK = evalin('base','RMSddzseCK'); RMSddZ1CK = evalin('base','RMSddZ1CK'); switch val %Z case hZ = plot(f0,GRMSZ,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('Z [m]'); title('DO DICH CHUYEN THAN XE'); axis tight ; axis([0 12 inf inf]); %dZ case 89 plot(f0,GRMSdZ,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('dZ [m/s]'); title('TOC DO DICH CHUYEN THAN XE'); axis tight; axis([0 12 inf inf]); %ddZ case plot(f0,GRMSddZ,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('ddZ [m/s^2]'); title('GIA TOC DICH CHUYEN THAN XE'); axis tight ; axis([0 12 inf inf]); %phi case plot(f0,GRMSphi,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('phi [rad]'); title('GOC DICH CHUYEN THAN XE'); axis tight; axis([0 12 inf inf]); %dphi case plot(f0,GRMSdphi,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('dphi [rad/s]'); title('TOC DO GOC DICH CHUYEN THAN XE'); axis tight; axis([0 12 inf inf]); %ddphi case plot(f0,GRMSddphi,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('ddphi [rad/s^2]'); title('GIA TOC GOC DICH CHUYEN THAN XE'); axis tight; axis([0 12 inf inf]); %z1 case plot(f0,GRMSz1,color,'linewidth',Linewidth); grid on; xlabel('Tan so ylabel('zse=z1 title('DO DICH axis tight; axis([0 12 inf (f) [Hz]'); [m]'); CHUYEN GHE LAI KHONG TREO'); inf]); 90 %dz1 case plot(f0,GRMSdz1,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('dz1 [m/s]'); title('TOC DO DICH CHUYEN GHE LAI KHONG TREO'); axis tight; axis([0 12 inf inf]); %ddz1 case plot(f0,GRMSddz1,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('ddz1 [m/s^2]'); title('GIA TOC DICH CHUYEN GHE LAI KHONG TREO'); axis tight; axis([0 12 inf inf]); %RMSddz1 case 10 plot(f0,GRMSzse,color,'linewidth',Linewidth); grid on; case 11 axis([0 12 inf inf]); xlabel('Tan so (f) [Hz]'); ylabel('zse [m/s^2]'); title('DO DICH CHUYEN TUONG DOI GHE LAI CO TREO'); axis tight; plot(f0,GRMSdzse,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('dzse [m/s^2]'); title('TOC DO DICH CHUYEN GHE LAI CO TREO'); axis tight ; axis([0 12 inf inf]); case 12 plot(f0,GRMSddzse,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('ddzse [m/s^2]'); title('GIA TOC DICH CHUYEN GHE LAI CO TREO'); axis tight; axis([0 12 inf inf]); case 13 plot(f0,GRMSFw1,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('Fw1 [N]'); title('LUC DONG BANH XE CAU SO 1'); axis tight; axis([0 12 inf inf]); 91 %f2 case 14 plot(f0,GRMSFw2,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('Fw2 [N]'); title('LUC DONG BANH XE CAU SO 2'); axis tight; axis([0 12 inf inf]); %f3 case 15 plot(f0,GRMSFw3,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('Fw3 [N]'); title('LUC DONG BANH XE CAU SO 3'); axis tight; axis([0 12 inf inf]); %f4 case 16 plot(f0,GRMSFw4,color,'linewidth',Linewidth); grid on; xlabel('Tan so (f) [Hz]'); ylabel('Fw4 [N]'); title('LUC DONG BANH XE CAU SO 4'); axis tight; axis([0 12 inf inf]); case 17 plot(f0,RMSddzseCK{1,1},' r',f0,RMSddzseCK{1,2},'b',f0,RMSddzseCK{1,3},'-.r','linewidth',1.5) legend('kse = 300','kse = 359','kse = 400'); title('GIA TOC BINH PHUONG TRUNG BINH GHE LAI') xlabel('f [Hz]') ylabel('RMSddzse [m/s^2]') hold on; grid on; axis tight; axis([0 12 inf inf]); case 18 plot(f0,RMSddzseCK{1,2},'-b',f0,RMSddZ1CK{1,2},'.r','linewidth',1.5) legend('Ghe lai co treo','Ghe lai khong treo'); title('SO SANH GIA TOC BINH PHUONG TRUNG BINH GHE LAI') ylabel('RMSddzse [m/s^2]') xlabel('f [Hz]') hold on; grid on; axis tight; axis([0 12 inf inf]); otherwise end function buttonNhapthongso_Callback(hObject, eventdata, handles) % hObject handle to buttonNhapthongso (see GCBO) 92 % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) inputFuntion = 2; %Parameter; Parameter; temp = who; for iv = 1:length(temp) assignin('base',temp{iv},eval(temp{iv})); end set(handles.editc1,'String',num2str(evalin('base','c1'))); set(handles.editc2,'String',num2str(evalin('base','c2'))); set(handles.editc3,'String',num2str(evalin('base','c3'))); set(handles.editc4,'String',num2str(evalin('base','c4'))); %set(handles.editcse,'String',num2str(evalin('base','cse'))); %display l set(handles.editl1,'String',num2str(evalin('base','l1'))); set(handles.editl2,'String',num2str(evalin('base','l2'))); set(handles.editl3,'String',num2str(evalin('base','l3'))); set(handles.editl4,'String',num2str(evalin('base','l4'))); set(handles.editlse,'String',num2str(evalin('base','lse'))); %display k set(handles.editk1,'String',num2str(evalin('base','k1'))); set(handles.editk2,'String',num2str(evalin('base','k2'))); set(handles.editk3,'String',num2str(evalin('base','k3'))); set(handles.editk4,'String',num2str(evalin('base','k4'))); %set(handles.editkse,'String',num2str(evalin('base','kse'))); %display k set(handles.editcw,'String',num2str(evalin('base','cw'))); %display M set(handles.editM,'String',num2str(evalin('base','M'))); %display m %display m set(handles.editm1,'String',num2str(evalin('base','m1'))); set(handles.editm2,'String',num2str(evalin('base','m2'))); set(handles.editm3,'String',num2str(evalin('base','m3'))); set(handles.editm4,'String',num2str(evalin('base','m4'))); set(handles.editmse,'String',num2str(evalin('base','mse'))); %display Jy set(handles.editJy,'String',num2str(evalin('base','Jy'))); %display g set(handles.editg,'String',num2str(evalin('base','g'))); %display che chuyen dong set(handles.editS,'String',num2str(evalin('base','S'))); set(handles.editq0,'String',num2str(evalin('base','qo'))); %set(handles.editVks,'String',num2str(evalin('base','Vks'))); set(handles.editVks,'String',num2str(evalin('base','fks'))); %evalin('base',Parameter.m); temp_1 = evalin('base','k1')+evalin('base','k2')+evalin('base','k3')+evalin('base' ,'k4');%tong he so can giam chan temp_2 = evalin('base','c1')+evalin('base','c2')+evalin('base','c3')+evalin('base' ,'c4'); % tong cung muy_eq = temp_1/(2*sqrt(evalin('base','g')*evalin('base','M')*temp_2)); %muy_eq wr = sqrt(temp_2/(evalin('base','g')*evalin('base','M'))); % tan so dao dong rieng khong giam chan 93 wd = wr*sqrt(1-muy_eq*muy_eq); % tan so dao dong rieng co giam chan assignin('base','wr',wr); assignin('base','wd',wd); set(handles.editTanso,'String',num2str(evalin('base','wd'))); function editFreq_CreateFcn(hObject, eventdata, handles) % hObject handle to editFreq (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) end set(hObject,'BackgroundColor','white'); Ph l c 3: K t qu tính tốn kho sát q0 = [0,04 0,06 0,08] m S = 4,24 m B ng 3.33_1: K t qu tính tốn gia tc BPTB (q0 = 0,04 m, S = 4,24 m) fks [Hz] 0,5 5,5 10 94 B ng 3.33_2: K t qu tính tốn l c ng bánh xe F (q0 = 0,04 m, S = 4,24 m) wi fks [Hz] 0,5 10 B ng 3.33_3: K t qu tính tốn gia tc BPTB (q0 = 0,06 m, S = 4,24 m) fks [Hz] 0,5 5,5 10 B ng 3.33_4: K t qu tính tốn l c ng bánh xe F (q0 = 0,06 m&S = 4, 24 m) wi fks [Hz] 0,5 10 B ng 3.33_5: K t qu tính tốn gia tc BPTB(q0 = 0,08 m, S = fks [Hz] 0,5 5,5 10 B ng 3.33_6: K t qu tính tốn kho sát F (q0 = 0,08 m, S = 4,24 m) wi fks [Hz] 0,5 10 97 ...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Ỹ NGHIÊM XUÂN THÀNH NGHIÊN CỨU ĐÁNH GIÁ ĐỘ ÊM DỊU CỦA XE THIẾT GIÁP BTR – 60PB BẰNG MATLAB – SIMULINK. .. “như sóng dịu êm vỗ bờ, hàng ngày hàng giờ, xô đến bến bờ nhân cách mới” Xin chân thành cảm ơn! iv TÓM TẮT Đề tài: "Nghiên cứu đánh giá độ êm dịu xe thiết giáp BTR - 60PB Matlab - Simulink" ,... lái xe 44 2.5.2.8 Các lực động bánh xe xuống đường 45 Chương 3: KHẢO SÁT ĐỘ ÊM DỊU XE BTR - 60PB BẰNG MATLAB SIMULINK 46 3.1 Giới thiệu hệ thống treo xe BTR - 60PB