BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - NGUYỄN KIM KỲ NGHIÊN CỨU ẢNH HƯỞNG CỦA THÔNG SỐ KẾT CẤU ĐẾN DAO ĐỘNG Ô TÔ Chuyên ngành: Kỹ thuật ôtô xe chuyên dụng LUẬN VĂN THẠC SĨ KHOA HỌC Kỹ thuật ôtô xe chuyên dụng NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS LƯU VĂN TUẤN Hà Nội – 2011 LỜI NÓI ĐẦU Ô tô phương tiện vận tải có vai trò quan trọng kinh tế quốc dân, có ưu điểm hẳn loại phương tiện vận tải khác như: thông dụng, đơn giản, dễ sử dụng có tính động cao ứng dụng rộng rãi tất lĩnh vực dân quốc phòng Cùng với phát triển mạnh mẽ khoa học công nghệ, công nghiệp ô tô giới phát triển ngày cao, cho đời nhiều loại ô tô đại phục vụ cho nhu cầu mục đích sử dụng người Việc thiết kế hệ thống, cụm chi tiết ôtô, đánh giá chất lượng làm việc chúng ngày quan tâm nghiên cứu nhằm nâng cao chất lượng ôtô Dao động ô tô vấn đề cần quan tâm nghiên cứu để nâng cao chất lượng khai thác ô tô Hiện phương pháp mô số sử dụng rộng rãi để nghiên cứu hệ, đặc biệt trình mô thể tính ưu việt cao miêu tả hệ phi tuyến áp dụng vào mô hệ dao động ô tô Trong thời gian vừa qua có nhiều tác giả đề cập tới vấn đề nghiên cứu dao động ô tô Tuy nhiên tính chất phức tạp vấn đề nên phần lớn tác giả phải sử dụng giả thiết để đơn giản hoá mô hình như: Chỉ nghiên cứu mô hình phẳng, hệ dao động thực tế chuyển thành hệ tuyến tính phi tuyến yếu… dẫn tới kết nghiên cứu nhiều hạn chế Bên cạnh đó, xe đóng Việt nam đặc biệt xe chở khách chưa nghiên cứu sâu dao động Xuất phát từ yêu cầu trên, lựa chọn đề tài: “Nghiên cứu ảnh hưởng thông số kết cấu đến dao động ô tô” Đề tài tập trung vào nghiên cứu mô hệ dao động ô tô có xét đến yếu tố phi tuyến từ làm sở cho việc chọn thông số kết cấu hợp lý đảm bảo độ êm dịu chuyển động cho xe Đây công việc hữu ích vô quan trọng cho việc nghiên cứu dao động ô tô trình thử nghiệm ban đầu nhằm tiết kiệm thời gian, khối lượng tính toán giảm chi phí cho thí nghiệm đo đạc sản phẩm thử nghiệm thực tế Trong trình thực luận văn, PGS_TS.Lưu Văn Tuấn tập thể cán giảng dạy môn Ôtô Trường Đại Học Bách Khoa Hà Nội bạn đồng nghiệp tận tình giúp đỡ bảo Tôi xin bày tỏ lòng biết ơn sâu sắc PGS_TS.Lưu Văn Tuấn thầy giáo môn Ôtô Xe chuyên dụng Trường Đại Học Bách Khoa Hà Nội tạo điều kiện thuận lợi để hoàn thành đề tài Mặc dù đề tài đạt số kết định, nhiên khả kinh nghiệm hạn chế nên luận văn không tránh khỏi sai sót Rất mong đóng góp ý kiến chân thành Thầy giáo, nhà khoa học, đồng nghiệp bạn bè để đề tài hoàn thiện Xin chân thành cảm ơn! Hà nội ngày 20 tháng 09 năm 2011 Học viên Nguyễn Kim Kỳ MỤC LỤC TRANG BÌA PHỤ LỜI CAM ĐOAN LỜI NÓI ĐẦU MỤC LỤC CÁC KÝ HIỆU VÀ CÁC THÔNG SỐ KỸ THUẬT CỦA XE DANH MỤC BẢNG .13 DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ 14 CHƯƠNG 1: TỔNG QUAN 16 1.1 Sự phát triển ngành công nghiệp ô tô Việt nam 16 1.1.1 Tình hình phát triển, quan điểm phát triển tương lai 16 1.1.2 Chính sánh nội địa hoá phát triển công nghiệp ô tô nước.17 a Quan điểm phát triển ngành công nghiệp ôtô: .17 b Mục tiêu phát triển: 19 c Đặc điểm thị trường ôtô Việt nam: 20 1.2 Tổng quan nghiên cứu dao động ô tô 21 1.2.1 Các nghiên cứu dao động ô tô nước 21 1.2.2 Nghiên cứu dao động ô tô giới 22 a Nghiên cứu biên dạng mặt đường: .23 b Nghiên cứu dao động ô tô: 23 c Nghiên cứu cảm giác người: 24 1.3 Nhiệm vụ phương hướng nghiên cứu luận văn 25 1.3.1 Nghiên cứu tiêu đánh độ êm dịu chuyển động 25 1.3.2 Xây dựng mô hình thiết lập hệ phương trình vi phân mô tả dao động 25 1.3.3 Mô dao động 25 1.4 Các tiêu đánh giá độ êm dịu chuyển động 26 1.4.1 Ảnh hưởng dao động tới thể người xe cộ 26 a Những ảnh hưởng dao động lên thể người 26 b Ảnh hưởng dao động đến độ bền độ an toàn chuyển động .28 1.4.2 Các tiêu đánh giá độ êm dịu chuyển động ô tô 29 a Chỉ tiêu tần số dao động 31 b Chỉ tiêu dịch chuyển tương đối sàn xe 31 c Chỉ tiêu gia tốc dao động 32 d Chỉ tiêu an toàn chuyển động tải trọng tác dụng xuống đường 33 CHƯƠNG 2: XÂY DỰNG MÔ HÌNH VÀ THIẾT LẬP HỆ PHƯƠNG TRÌNH VI PHÂN MÔ TẢ DAO ĐỘNG 35 2.1 Xây dựng mô hình dao động ô tô 35 2.1.1 Mô hình dao động 35 a Mô hình 1/4 35 b Mô hình dao động liên kết .36 c Mô hình 1/2 37 d Mô hình không gian xe 37 2.1.2 Xây dựng mô hình dao động 38 2.1.3 Hàm kích động 39 2.2 Thiết lập hệ phương trình vi phân mô tả dao động 40 2.2.1 Hệ thống treo trước 44 a Khối lượng không treo trước bên trái: .44 b Khối lượng không treo trước bên phải: 45 2.2.2 Hệ thống treo sau 47 a Lực quán tính: 47 b Trọng lực: 47 d Lực truyền từ thân xe xuống cầu qua hệ thống treo trái: 48 e Lực truyền từ thân xe xuống cầu qua hệ thống treo phải: .48 f Phương trình cân lực tác dụng lên cầu sau theo phương thẳng đứng Z: 49 2.2.3 Thân xe 49 a Phương trình cân lực theo phương Z: 49 b Phương trình cân mômen theo toạ độ suy rộng ϕ : .50 c Phương trình cân mômen theo toạ độ suy rộng θ : 50 2.2.4 Hệ phương trình vi phân 51 CHƯƠNG 3: MÔ PHỎNG DAO ĐỘNG 52 3.1 Thư viện tiện ích Simulink 52 3.2 Sơ đồ tổng thể 55 3.3 Sơ đồ khối chức 56 3.3.1 Khối “Thân xe” .56 a Khối “Z”: 57 b Khối “Phi”: 57 c Khối “Teta”: 58 3.3.2 Khối “Hệ thống treo” 59 a Khối “F_Suspension”: 59 b Khối “R_Suspension”: 63 c Khối “gia tri Zi”: 64 3.3.3 Khối “Cầu xe” 65 a Khối “Cau_truoc”: 65 b Khối “Cau_sau”: 67 3.3.4 Khối lốp xe 68 3.3.5 Khối “Duong” 70 CHƯƠNG 4: MỘT SỐ KẾT QUẢ MÔ PHỎNG VÀ KẾT LUẬN 73 4.1 Các thông số xe đưa vào tính toán 73 4.2 Một số kết mô 74 4.2.1 Một số kết mô dao động 74 4.2.2 Ảnh hưởng thông số đến dao động xe 78 a Ảnh hưởng vận tốc xe: 78 b Ảnh hưởng thông số hệ thống treo: 83 c Ảnh hưởng độ cứng lốp: 87 d Ảnh hưởng khối lượng thân xe: 88 e Ảnh hưởng độ cứng ổn định ngang: 90 KẾT LUẬN……………………………………………………………………… 92 TÀI LIỆU THAM KHẢO……………………………………………………… 94 CÁC KÝ HIỆU VÀ CÁC THÔNG SỐ KỸ THUẬT CỦA XE Ký hiệu Giá trị Đơn vị A 1,91 M Khoảng cách từ trọng tâm xe tới cầu trước B 1,64 M Khoảng cách từ trọng tâm xe tới cầu sau B1 1,32 M B2 1,33 M Chiều rộng sở cầu sau B2g 1,16 M Khoảng cách giảm chấn treo sau C1T 29781 N/m Độ cứng hệ thống treo trước trái C1P 29781 N/m Độ cứng hệ thống treo trước phải C2T 46446 N/m Độ cứng hệ thống treo sau trái C2P 46446 N/m Độ cứng hệ thống treo sau phải Cod 89962 Nm/rad Độ cứng chống xoắn ổn định CL1T 152000 N/m Độ cứng lốp trước trái CL1P 152000 N/m Độ cứng lốp trước phải CL2T 310000 N/m Độ cứng lốp sau trái CL2P 310000 N/m Độ cứng lốp sau phải ft 0,16 M Độ võng tĩnh hệ thống treo đầy tải Fqt1T, Fqt1P N Ý nghĩa ký hiệu Chiều rộng sở cầu trước Lực quán tính khối lượng không treo trước trái, phải Fqt2 N Lực quán tính khối lượng không treo sau Ftl1T, Ftl1P N Ftl2 N Trọng lực khối lượng không treo sau Fqt3 N Lực quán tính khối lượng treo Ftl3 N Trọng lực khối lượng treo FC1T, FC1P N Lực đàn hồi hệ thống treo trước trái, phải FK1T, FK1P N Trọng lực khối lượng không treo trước trái, phải Lực cản giảm chấn hệ thống treo trước trái, phải FCL1T, Lực đàn hồi bánh xe trước trái, phải N FCL1P FKL1T, Lực cản giảm chấn bánh xe trước trái, phải N FKL1P FC2T, FC2P N Lực đàn hồi hệ thống treo sau trái, phải FK2T, FK2P N Lực cản giảm chấn hệ thống treo sau trái, phải FCL2T, Lực đàn hồi bánh xe sau trái, phải N FCL2P FKL2T, Lực cản giảm chấn bánh xe sau trái, phải N FKL2P F1T, F1P N Lực truyền từ thân xe xuống cầu trước trái, phải FL1T, FL1P N Lực truyền từ lốp lên cầu trước trái, phải F2T, F2P N FL2T, FL2P N Lực truyền từ thân xe xuống hệ thống treo sau trái, phải Lực truyền từ lốp lên hệ thống treo sau trái, phải g 9,81 m/s Jx 678 Nms2 Jy 23155 Nms2 J2x 41,4 Nms2 Gia tốc trọng trường Mô men quán tính khối lượng treo theo trục quay x Mô men quán tính khối lượng treo theo trục quay y Mô men quán tính khối lượng không treo sau theo trục quay x K1T, K1P Ns/m Hệ số cản giảm chấn trước trái, phải K2T, K2P Ns/m Hệ số cản giảm chấn sau trái, phải K1nm 193 Ns/m Hệ số cản giảm chấn trước nén mạnh K1nn 492,5 Ns/m Hệ số cản giảm chấn trước nén nhẹ K1tm 509 Ns/m Hệ số cản giảm chấn trước trả mạnh K1tn 1517,5 Ns/m Hệ số cản giảm chấn trước trả nhẹ K2nm 245 Ns/m Hệ số cản giảm chấn sau nén mạnh 10 K2nn 625 Ns/m Hệ số cản giảm chấn sau nén nhẹ K2tm 740 Ns/m Hệ số cản giảm chấn sau trả mạnh K2tn 1835 Ns/m Hệ số cản giảm chấn sau trả nhẹ Kođ 14318 N/m Hệ số tác dụng ổn định KL1T 1260 Ns/m Hệ số cản lốp trước trái KL1P 1260 Ns/m Hệ số cản lốp trước phải KL2T 1567 Ns/m Hệ số cản lốp sau trái KL2P 1567 Ns/m Hệ số cản lốp sau phải L 3.55 m Chiều dài sở xe Mxe 2120 kg Khối lượng toàn xe đầy tải M 1950 kg Khối lượng treo đầy tải m1T 40 kg Khối lượng không treo cầu trước trái m1P 40 kg Khối lượng không treo cầu trước phải m2 90 kg Khối lượng không treo cầu sau Mô men quán tính khối lượng treo theo M qtφ Nm M qtθ Nm M1, M2 Nm M1T, M1P Nm Mô men tác truyền từ hệ thống treo trước trái, phải M2T, M2P Nm Mô men tác truyền từ hệ thống treo sau trái, phải MT, MP Nm q1T, q1P m q2T, q2P m t s trục quay y Mô men quán tính khối lượng treo theo trục quay x Mô men tác dụng lên thân xe truyền từ hệ thống treo trước, sau Mô men tác dụng lên thân xe truyền từ hệ thống treo phía trái, phải Kích động mặt đường tác động lên lốp trước trái, phải Kích động mặt đường tác động lên lốp sau trái, phải Thời gian dao động 11 KL1P=KL1T; KL2T=1567; KL2P=KL2T; Kod=Cod/(2*pi*0.1^2); L=3.55; M=2120; m1T=40; m1P=m1T; m2=90; v=60/3.6; - Ảnh hưởng vận tốc xe đến độ lệch quân phương toạ độ trọng tâm thân xe %DolechZ_V: clear all clc nhap; km1; V=10:5:100; nV=length(V); DlZ=ones(1,nV); for j= 1:nV v=V(1,j)/3.6; sim(‘Mo_phong') suZ=sum(Z); nZ=length(Z); phuongsaiZ=0; for i= 2000:nZ Ztb=suZ/nZ; phuongsaiZ=phuongsaiZ+(Z(i,1)-Ztb)^2/nZ; end viii DlZ(1,j)=sqrt(phuongsaiZ); end plot(V,DlZ); xlabel('Van toc xe (Km/h)') ylabel('Sai lech quan phuong toa tam than xe (m)') title('Anh huong cua van toc xe den sai lech quan phuong toa tam than xe Z') grid on - Ảnh hưởng vận tốc xe đến sai lệch quân phương gia tốc trọng tâm thân xe %DolechZ2dot_V: clear all clc nhap; km1; V=10:5:100; nV=length(V); DlZ2dot=ones(1,nV); for j= 1:nV v=V(1,j)/3.6; sim(‘Mo_phong') suZ2dot=sum(Z2dot); nZ2dot=length(Z2dot); phuongsaiZ2dot=0; for i= 2000:nZ2dot Z2dottb=suZ2dot/nZ2dot; phuongsaiZ2dot=phuongsaiZ2dot+(Z2dot(i,1)-Z2dottb)^2/nZ2dot; end DlZ2dot(1,j)=sqrt(phuongsaiZ2dot); end ix plot(V,DlZ2dot); xlabel('Van toc xe (Km/h)') ylabel('Sai lech quan phuong gia toc tam than xe (m/s2)') title('Anh huong van toc xe den sai lech quan phuong gia toc toa than xe') grid on - Ảnh hưởng vận tốc xe đến sai lệch quân phương góc lắc thân xe DPhi, DTeta (Rad): %DolechPhi_V clear all clc nhap; km1; V=10:5:90; nV=length(V); Dlphi=ones(1,nV); for j= 1:nV v=V(1,j)/3.6; sim('Mo_phong') suphi=sum(phi); nphi=length(phi); phuongsaiphi=0; for i= 2000:nphi phitb=suphi/nphi; phuongsaiphi=phuongsaiphi+(phi(i,1)-phitb)^2/nphi; end Dlphi(1,j)=sqrt(phuongsaiphi); end plot(V,Dlphi); xlabel('Van toc xe (Km/h)') ylabel('Sai lech quan phuong goc lac doc than xe (Rad)') x title('Anh huong cua van toc xe den sai lech quan phuong goc lac doc than xe (Phi)') grid on %DolechTeta_V clear all clc nhap; km1; V=10:5:90; nV=length(V);Dlteta=ones(1,nV); for j= 1:nV v=V(1,j)/3.6; sim('Mo_phong') suTeta=sum(Teta); nTeta=length(Teta); phuongsaiTeta=0; for i= 2000:nTeta Tetatb=suTeta/nTeta; phuongsaiTeta=phuongsaiTeta+(Teta(i,1)-Tetatb)^2/nTeta; end DlTeta(1,j)=sqrt(phuongsaiTeta); end plot(V,DlTeta); xlabel('Van toc xe (Km/h)') ylabel('Sai lech quan phuong goc lac ngang than xe (Rad)') title('Anh huong cua van toc xe den sai lech quan phuong goc lac ngang than xe (Teta)') grid on xi - Ảnh hưởng vận tốc xe đến sai lệch toạ độ trọng tâm cầu xe DCsi1T, DCsi1P, DCsi2: %DolechCsi1T_V clear all clc nhap; km1; V=10:5:100; nV=length(V); DlCsi1T=ones(1,nV) for j= 1:nV v=V(1,j)/3.6; sim('Mo_phong') suCsi1T=sum(Csi1T); nCsi1T=length(Csi1T); phuongsaiCsi1T=0; for i= 2000:nCsi1T Csi1Ttb=suCsi1T/nCsi1T; phuongsaiCsi1T=phuongsaiCsi1T+(Csi1T(i,1)-Csi1Ttb)^2/nCsi1T; end DlCsi1T(1,j)=sqrt(phuongsaiCsi1T); end plot(V,DlCsi1T); xlabel('Van toc xe (Km/h)') ylabel('Sai lech quan phuong toa tam cau xe (m)') title('Anh huong cua van toc xe den sai lech quan phuong toa tam cau xe') grid on; hold on; xii - Ảnh hưởng độ cứng lò xo đến sai lệch quân phương toạ độ trọng tâm thân xe DZ (m) gia tốc thân xe DZ&& : %DolechZ_C1_C2 clear all clc nhap; km1; mau=['r' 'b' 'g']; C1=20000:15000:50000; C2=20000:2000:90000; nC2=length(C2); DlZ=ones(1,nC2); for k= 1:3 C1T=C1(1,k); C1P=C1(1,k); for j= 1:nC2 C2T=C2(1,j); C2P=C2(1,j); sim('Mo_phong') suZ=sum(Z); nZ=length(Z); phuongsaiZ=0; for i= 1:nZ Ztb=suZ/nZ; phuongsaiZ=phuongsaiZ+(Z(i,1)-Ztb)^2/nZ; end DlZ(1,j)=sqrt(phuongsaiZ); end plot(C2,DlZ,mau(k)); Hold on xiii end xlabel('Do cung cua lo xo sau C2 (N/m)') ylabel('Sai lech quan phuong toa tam than xe (m)') title('Anh huong cung lo xo C2 den sai lech quan phuong toa tam than xe C1 thay doi') grid on %DolechZdot_C1_C2 clear all clc nhap; km1; mau=['r' 'b' 'g']; C1=20000:15000:50000; C2=20000:2000:90000; nC2=length(C2); DlZ2dot=ones(1,nC2); for k= 1:3 C1T=C1(1,k); C1P=C1(1,k); for j= 1:nC2 C2T=C2(1,j); C2P=C2(1,j); sim('Mo_phong') suZ2dot=sum(Z2dot); nZ2dot=length(Z2dot); phuongsaiZ2dot=0; for i= 1:nZ2dot Z2dottb=suZ2dot/nZ2dot; phuongsaiZ2dot=phuongsaiZ2dot+(Z2dot(i,1)-Z2dottb)^2/nZ2dot; end xiv DlZ2dot(1,j)=sqrt(phuongsaiZ2dot); end plot(C2,DlZ2dot,mau(k)); Hold on end xlabel('Do cung cua lo xo sau C2 (N/m)') ylabel('Sai lech quan phuong gia toc tam than xe (m)') title('Anh huong cung lo xo C2 den sai lech quan phuong gia tốc tam than xe C1 thay doi') grid on - Ảnh hưởng hệ số cản giảm chấn đến sai lệch quân phương toạ độ trọng tâm thân xe DZ (m) gia tốc thân xe DZ&& : %DolechZ_K: clear all clc nhap; km1; hs=0.6:0.1:1.6; nhs=length(hs); DlZ=ones(1,nhs); for j= 1:nhs K1nn=492.5*hs(1,j); K1nm=193*hs(1,j); K1tn=1517.5*hs(1,j); K1tm=609*hs(1,j); K2nn=625*hs(1,j); K2nm=245*hs(1,j); K2tn=1835*hs(1,j); K2tm=740*hs(1,j); sim('Mo_phong') xv suZ=sum(Z); nZ=length(Z); phuongsaiZ=0; for i= 2000:nZ Ztb=suZ/nZ; phuongsaiZ=phuongsaiZ+(Z(i,1)-Ztb)^2/nZ; end DlZ (1,j)=sqrt(phuongsaiZ); end plot(hs,DlZ); xlabel('He so can cua giam chan X*K') ylabel('Sai lech quan phuong toa tam than xe (m/s2)') title('Anh huong can giam chan den sai lech quan phuong toa tam than xe') grid on %DolechZ2dot_K: clear all clc nhap; km1; hs=0.6:0.1:1.6; nhs=length(hs); DlZ2dot=ones(1,nhs); for j= 1:nhs K1nn=492.5*hs(1,j); K1nm=193*hs(1,j); K1tn=1517.5*hs(1,j); K1tm=609*hs(1,j); K2nn=625*hs(1,j); K2nm=245*hs(1,j); xvi K2tn=1835*hs(1,j); K2tm=740*hs(1,j); sim('Mo_phong') suZ2dot=sum(Z2dot); nZ2dot=length(Z2dot); phuongsaiZ2dot=0; for i= 2000:nZ2dot Z2dottb=suZ2dot/nZ2dot; phuongsaiZ2dot=phuongsaiZ2dot+(Z2dot(i,1)-Z2dottb)^2/nZ2dot; end DlZ2dot(1,j)=sqrt(phuongsaiZ2dot); end plot(hs,DlZ2dot); xlabel('He so can cua giam chan X*K') ylabel('Sai lech quan phuong gia toc tam than xe (m/s2)') title('Anh huong can giam chan den sai lech quan phuong gia toc tam than xe') grid on - Ảnh hưởng độ cứng lốp %DolechZ_CLop1_CLop2: clear all clc nhap; km1; Cl2=200000:100000:500000; Cl1=100000:10000:250000; nCl1=length(Cl1); DlZ=ones(1,nCl1); mau=['r' 'b' 'g']; for k= 1:4 xvii CL2T=Cl2(1,k); CL2P=Cl2(1,k); for j= 1:nCl1 CL1T=Cl1(1,j); CL1P=Cl1(1,j); sim('Mo_phong') suZ=sum(Z); nZ=length(Z); phuongsaiZ=0; for i= 2000:nZ Ztb=suZ/nZ; phuongsaiZ=phuongsaiZ+(Z(i,1)-Ztb)^2/nZ; end DlZ(1,j)=sqrt(phuongsaiZ); end plot(Cl1,DlZ,mau(k)); hold on end xlabel('Do cung cua lop truoc (N/m)') ylabel('Sai lech quan phuong toa tam than xe (m)') title('Anh huong cung lop den sai lech quan phuong toa tam than xe') grid on %DolechZdot_CLop1_CLop2: clear all clc nhap; km1; Cl2=200000:100000:500000; Cl1=100000:10000:250000; xviii nCl1=length(Cl1); DlZ2dot=ones(1,nCl1); mau=['r' 'b' 'g']; for k= 1:4 CL2T=Cl2(1,k); CL2P=Cl2(1,k); for j= 1:nCl1 CL1T=Cl1(1,j); CL1P=Cl1(1,j); sim('Mo_phong') suZ2dot=sum(Z2dot); nZ2dot=length(Z2dot); phuongsaiZ2dot=0; for i= 2000:nZ2dot Z2dottb=suZ2dot/nZ2dot; phuongsaiZ2dot=phuongsaiZ2dot+(Z2dot (i,1)-Z2dottb)^2/nZ2dot; end DlZ2dot (1,j)=sqrt(phuongsaiZ2dot); end plot(Cl1,DlZ2dot,mau(k)); hold on end xlabel('Do cung cua lop truoc (N/m)') ylabel('Sai lech quan phuong gia toc tam than xe (m)') title('Anh huong cung lop den sai lech quan phuong gia toc tam than xe') grid on - Ảnh hưởng khối lượng thân xe %DolechZ_M: clear all clc xix nhap; km1; MT=1000:140:2120; nMT=length(MT); DlZ2=ones(1,nMT); for j= 1:nMT M=MT(1,j); b=b-0.0163; a=a+0.0163; sim('Mo_phong') suZ=sum(Z); nZ=length(Z); phuongsaiZ=0; for i= 2000:nZ Ztb=suZ2/nZ; phuongsaiZ=phuongsaiZ+(Z(i,1)-Ztb)^2/nZ; end DlZ(1,j)=sqrt(phuongsaiZ); end plot(MT,DlZ); xlabel('So luong hanh khach (nguoi)') ylabel('Sai lech quan phuong tọa tam than xe ') title('Anh huong khoi luong than xe den sai lech quan phuong toa than xe') grid on %DolechZ2dot_M: clear all clc nhap; km1; MT=1000:140:2120; xx nMT=length(MT); DlZ2dot=ones(1,nMT); for j= 1:nMT M=MT(1,j); b=b-0.0163; a=a+0.0163; sim('Mo_phong') suZ2dot=sum(Z2dot); nZ2dot=length(Z2dot); phuongsaiZ2dot=0; for i= 2000:nZ2dot Z2dottb=suZ2dot/nZ2dot; phuongsaiZ2dot=phuongsaiZ2dot+(Z2dot(i,1)-Z2dottb)^2/nZ2dot; end DlZ2dot(1,j)=sqrt(phuongsaiZ2dot); end plot(MT,DlZ2dot); xlabel('So luong hanh khach (nguoi)') ylabel('Sai lech quan phuong gia toc tam than xe (m/s2)') title('Anh huong khoi luong than xe den sai lech quan phuong gia toc than xe') grid on - Ảnh hưởng ổn định ngang %DolechTeta_Cod: clear all clc nhap; km1; Cod=40000:2000:60000; nCod=length(Cod); DlTeta=ones(1,nCod); for j= 1:nCod xxi kod=Cod(1,j)/(2*pi*0.1^2); sim('Mo_phong') suTeta=sum(Teta); nTeta=length(Teta); phuongsaiTeta=0; for i= 2000:nTeta Tetatb=suTeta/nTeta; phuongsaiTeta=phuongsaiTeta+(Teta(i,1)-Tetatb)^2/nTeta; end DlTeta(1,j)=sqrt(phuongsaiTeta); end plot(Cod,DlTeta); xlabel('Do cung chong xoan cua on dinh ngang (N/m)') ylabel('Sai lech quan phuong goc lac ngang than xe Teta (Rad)') title('Anh huong cua cung on dinh den sai lech quan phuong goc lac ngang than xe') grid on xxii ... chưa nghiên cứu sâu dao động Xuất phát từ yêu cầu trên, lựa chọn đề tài: Nghiên cứu ảnh hưởng thông số kết cấu đến dao động ô tô Đề tài tập trung vào nghiên cứu mô hệ dao động ô tô có xét đến. .. ô tô 21 1.2.1 Các nghiên cứu dao động ô tô nước 21 1.2.2 Nghiên cứu dao động ô tô giới 22 a Nghiên cứu biên dạng mặt đường: .23 b Nghiên cứu dao động ô tô: 23 c Nghiên. .. sỹ tác giả Đặng Việt Hà [10] Nghiên cứu ảnh hưởng thông số kết cấu đến điều kiện làm việc ô tô độ êm dịu chuyển động ý đến ảnh hưởng thông số kết cấu đến dao động xe khách đóng xe tải Nhìn