(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường

(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường(Luận văn thạc sĩ) Nghiên cứu ảnh hưởng của tải trọng động đến an toàn chuyển động của loại xe khách khi chuyển động trên đường

LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu Các số liệu, kết nêu luận văn trung thực chƣa đƣợc cơng bố cơng trình khác Tp Hồ Chí Minh, ngày … tháng … năm 201… (Ký tên ghi rõ họ tên) Lê Thiện Kim Ngân trang LỜI CẢM ƠN Tôi xin chân thành cảm ơn quý thầy cô trƣờng Đại học Sƣ Phạm Kỹ Thuật Tp.HCM, đặc biệt quý thầy cô khoa Cơ khí Động Lực giảng dạy, tạo điều kiện cho tơi hồn thành khóa học Tơi xin chân thành cảm ơn Thầy PGS.TS Nguyễn Văn Phụng, ngƣời hết lịng tận tình hƣớng dẫn tơi thực luận văn Tôi xin gửi lời cám ơn đến ban lãnh đạo Trung Tâm Ứng dụng Tiến Khoa học Công Nghệ - Sở Khoa học Công Nghệ tỉnh Đồng Nai tạo điều kiện cho suốt trình học tập Mặc dù luận văn hồn thành nhƣng chắn cịn nhiều thiếu sót, mong đƣợc đóng góp ý kiến quý thầy để luận văn đƣợc hồn thiện Học viên Lê Thiện Kim Ngân trang TÓM TẮT Việc xác định tải trọng động tác dụng lên thân xe có ý nghĩa vai trị vơ quan trọng Đề tài đƣợc tiến hành nghiên cứu theo phƣơng pháp đƣợc áp dụng giới nghiên cứu dao động ô tô vận dụng tối đa thành tựu công nghệ thông tin để rút ngắn thời gian tính tốn, tăng độ tin cậy suất tính tốn Luận văn tính tốn xuất phát từ việc xây dựng mơ hình tốn học, thiết lập phƣơng trình vi phân mơ tả chuyển động tơ đƣợc mơ hình hóa, tốn dao động ô tô đƣợc giải chƣơng trình máy tính cho kết nhanh chóng xác Luận văn sử dụng phần mềm Matlab để giải tốn dao động tơ Mục đích luận văn nghiên cứu ảnh hƣởng tải trọng động đến an toàn chuyển động xe, dựa vào kết tính tốn ta xác định đƣợc điều kiện để xe chạy an toàn đƣờng, xác định đƣợc vận tốc tới hạn độ cao mấp mô nguy hiễm mà xe bị lái, phanh lực kéo SUMMARY The determination of dynamic loads acting on the vehicle body significance and crucial role Subject is studied by the method being applied worldwide ranged from a study of car is the maximum use of information technology achievements to shorten calculation time, increased reliability and yield calculations Thesis calculations derived from the construction of mathematical models, set the differential equations describing the motion of the cars have been modeled, automotive oscillation problem will be solved by a computer program the results quickly and accurately This paper uses Matlab software to solve the problem ranged cars The aim of the thesis was to study the effects of dynamic loads to safe movement of the car, based on the results of calculations can determine the conditions for safe movement of vehicles on the road, determine the velocity to trang term and dangerous bumpy altitude at which the vehicle can lose steering, braking or loss took scissors trang MỤC LỤC Trang tựa Quyết định giao đề tài .2 Lý lịch cá nhân Lời cam đoan .4 Lời cảm ơn Tóm tắt Mục lục Danh mục hình……………………………………………………………… 11 DANH MỤC CÁC BẢNG 14 CHƢƠNG TỔNG QUAN 15 1.1 Đặt vấn đề 15 1.2 Tình hình nghiên cứu ngồi nƣớc: 16 1.3 Mục đích nghiên cứu 17 1.4 Đối tƣợng nghiên cứu: 17 1.5 Phạm vi nghiên cứu 18 1.6 Giá trị thực tiễn đề tài 18 1.7 Phƣơng pháp nghiên cứu 18 1.8 Nội dung đề tài 19 CHƢƠNG CƠ SỞ LÝ THUYẾT 20 2.1 Kết cấu hệ dao động ôtô 20 2.1.1 Khối lƣợng đƣợc treo 20 2.1.2 Khối lƣợng không đƣợc treo 20 2.1.3 Hệ thống treo 20 2.1.4 Lốp 21 2.2 Thông số kết cấu, tác dụng phận hệ thống treo 21 2.2.1 Thành phần đàn hồi 21 2.2.2 Bộ phận giảm chấn 22 trang 2.2.3 Bộ phận dẫn hƣớng 23 2.2.4 Sự đàn hồi lốp 23 2.3 Chỉ tiêu an toàn chuyển động tải trọng tác dụng xuống đƣờng 23 2.5 Các khái niệm tƣơng đƣơng 29 2.6 Các giả thiết 29 2.7 Mơ hình khơng gian 29 2.8 Các chế độ tải trọng tác dụng 31 2.9 Xác định moment quán tính khối lƣợng 32 2.9.1 Xác định moment quán tính khối lƣợng phần đƣợc treo Ix, Iy 32 2.9.2 Tính Ixra monen quán tính khối lƣợng cụm bánh xe sau cầu sau 33 2.10 Tính tốn độ cứng nhíp: 34 2.11 Xác định độ cứng lốp xe ktf , k tr 35 2.12 Biên dạng mặt đƣờng dạng bán bình phƣơng hàm sin 36 CHƢƠNG 3THƠNG SỐ TÍNH TỐN MƠ PHỎNG 37 3.1 Kết cấu xe Hyundai County 37 3.2 Thông số kỹ thuật xe Hyundai County 37 3.3 Tính tốn độ cứng nhíp 39 3.4 Xác định độ cứng lốp xe ktf , k tr 40 3.5 Tính momen quán tính khối lƣợng phần đƣợc treo 41 3.6 Tính momen quán tính khối lƣợng I xra cầu sau 42 3.7 Thông số tính tốn xe Huyndai County 43 3.8 Tính tốn ổn định 44 3.8.1 Xác định bán kính quay vịng ô tô 44 3.8.2 Tính ổn định dọc tơ 46 3.8.3 Tính ổn định ngang ô tô 47 3.9 Áp dụng phƣơng trình Lagrange II để tính tốn tải trọng động 49 Động K hệ: 49 Thế V hệ: 50 Hàm tiêu tán D hệ: 50 trang 10 3.10 Phƣơng pháp xác định tải trọng động 53 3.11 Kết tính tốn cho trƣờng hợp biên dạng mặt đƣờng dạng bán bình phƣơng hàm sin 53 Tần số tự nhiên dạng dao động: 53 Động lực học: 60 KẾT LUẬN VÀ HƢỚNG PHÁT TRIỂN ĐỀ TÀI 65 4.1 Kết luận 65 4.2 Hƣớng phát triển đề tài………………………………….………………… 66 trang 11 DANH MỤC HÌNH HÌNH TRANG Hình 1.1 Bố trí chung tơ Huyndai County 16 Hình 2.1 Quan hệ tần số dao động riêng n phần đƣợc treo với độ võng tĩnh ft .20 Hình 2.2 Sự biến thiên tải trọng thẳng đứng theo thời gian 20 Hình 2.3 Sự phân bố xác xuất tải trọng thẳng đứng đặt lên bánh xe .22 Hình 2.4Mơ hình khơng gian hệ giao động tơ khách 28 Hình 2.5 Chế độ tải trọng uốn 29 Hình 2.6 Chế độ tải trọng xoắn 30 Hình 2.7 Khối hộp thân xe chữ nhật đồng chất 31 Hình 2.8 Mơ hình giả định tính tốn momen qn tính khối lƣợng cầu sau 32 Hình 2.9 Thơng số kích thƣớc lò xo 32 Hình 2.10 Mặt cắt lốp xe 33 Hình 2.11 Biên dạng mặt đƣờng dạng bán bình phƣơng hàm sin 34 Hình 3.1 Kết cấu xe Hyundai County 35 Hình 3.2 Thơng số kích thƣớc lị xo 37 Hình 3.3 Mặt cắt lốp xe 38 Hình 3.4 Khối hộp thân xe đồng chất .39 Hình 3.5 Mơ hình giả định tính momen quán tính khối lƣợng cầu sau 40 Hình 3.6 Sơ đồ động học quay vịng ô tô 45 Hình 3.7 Sơ đồ tính tốn ổn định dọc xe lên dốc .46 Hình 3.8 Sơ đồ tính tốn ổn định ngang ô tô 47 Hình 3.9 Các dạng dao động thứ : 3,7,4 55 Hình 3.10 Biên độ gia tốc x theo tần số lực kích động, 56 Hình 3.11 Gia tốc theo tần số lực kích động, 56 Hình 3.12 Biên độ gia tốc theo tần số lực kích động 57 Hình 3.13 Gia tốc x theo thời gian 58 Hình 3.14 Gia tốc theo thời gian 59 trang 12 Hình 3.15 Gia tốc theo thời gian 60 Hình 3.16 Xung biên độ tải trọng động .61 trang 13 DANH MỤC CÁC BẢNG BẢNG TRANG Bảng 2.1 Liệt kê số hàm biểu diễn profil phổ biến 27 Bảng 3.1 Thông số xe Hyundai County 37 Bảng 3.2 Thơng số lị xo .38 Bảng 3.3 Thơng số tính tốn xe Huyndai County 43 Bảng 3.4 Tổng hợp tần số dao động tự nhiên hệ tọa độ tƣơng ứng 55 trang 14 cc(7,5)=cc(5,7); cc(7,6)=cc(6,7); cc(7,7)=2*cr; C=cc; % -dieu kien bien y1=0; y2=y1; y3=0;y4=0; % -F(1,1)=0;F(2,1)=0;F(3,1)=0;F(4,1)=0; F(5,1)=y1*ktf;F(6,1)=y2*ktf;F(7,1)=(y3+y4)*ktr; dy=zeros(14,1); dy(1)=y(8); dy(2)=y(9); dy(3)=y(10); dy(4)=y(11); dy(5)=y(12); dy(6)=y(13); dy(7)=y(14); dy(8)=-1/M(1,1)*(C(1,:)*y(8:14)+K(1,:)*y(1:7)-F(1,1)); dy(9)=-1/M(2,2)*(C(2,:)*y(8:14)+K(2,:)*y(1:7)-F(2,1)); dy(10)=-1/M(3,3)*(C(3,:)*y(8:14)+K(3,:)*y(1:7)-F(3,1)); dy(11)=-1/M(4,4)*(C(4,:)*y(8:14)+K(4,:)*y(1:7)-F(4,1)); dy(12)=-1/M(5,5)*(C(5,:)*y(8:14)+K(5,:)*y(1:7)-F(5,1)); dy(13)=-1/M(6,6)*(C(6,:)*y(8:14)+K(6,:)*y(1:7)-F(6,1)); dy(14)=-1/M(7,7)*(C(7,:)*y(8:14)+K(7,:)*y(1:7)-F(7,1)); end %gia tốc x theo thời gian trang 77 clear clc;close all fn=1.45; %chon gia tri lon nhat d1=1;d2=0.25; vx=2*d1*fn/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1/vx; t_step=0.001; t_smooth_end=4; t_bump_state=t_start:t_step:t_end-t_step; %gan het bien dang t_smooth_state=t_end:t_step:t_smooth_end; x_initial=zeros(14,1); %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]= ode45(@ptvp_xe_chiu_uon,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]= ode45(@ptvp_xe_chiu_xoan,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]= ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_smooth_state, x_initial_bending); [T_smooth_state Y_smooth_state_torsion]= ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_smooth_state, x_initial_torsion); ddx_bend=[Y_bump_state_bending(:,8);Y_smooth_state_bending(:,8)]; ddx_tor=[Y_bump_state_torsion(:,8);Y_smooth_state_torsion(:,8)]; time=[T_bump_state;T_smooth_state]; axes('FontSize',16); trang 78 hold on plot(time,ddx_bend,'-b','LineWidth',4); plot(time,ddx_tor,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('t, s'); ylabel('Accel of x, m/s^2 '); title('GIA TOC DDX THEO THOI GIAN'); %gia tốc theo thời gian clear clc;close all fn=1.45; %chon gia tri lon nhat d1=1;d2=0.25; vx=2*d1*fn/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1/vx; t_step=0.001; t_smooth_end=4; t_bump_state=t_start:t_step:t_end-t_step; %gan het bien dang t_smooth_state=t_end:t_step:t_smooth_end; x_initial=zeros(14,1); %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]= ode45(@ptvp_xe_chiu_uon,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]= ode45(@ptvp_xe_chiu_xoan,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]= trang 79 ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_smooth_state, x_initial_bending); [T_smooth_state Y_smooth_state_torsion]= ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_smooth_state, x_initial_torsion); ddphi_bend=[Y_bump_state_bending(:,9);Y_smooth_state_bending(:,9)]; ddphi_tor=[Y_bump_state_torsion(:,9);Y_smooth_state_torsion(:,9)]; time=[T_bump_state;T_smooth_state]; axes('FontSize',16); hold on plot(time,ddphi_bend,'-b','LineWidth',4); plot(time,ddphi_tor,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('t, s'); ylabel('Accel of phi, rad/s^2 '); title('GIA TOC DDPHI THEO THOI GIAN'); %gia tốc theo thời gian clear clc;close all fn=1.45; %chon gia tri lon nhat d1=1;d2=0.25; vx=2*d1*fn/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1/vx; t_step=0.001; t_smooth_end=5; t_bump_state=t_start:t_step:t_end-t_step; %gan het bien dang t_smooth_state=t_end:t_step:t_smooth_end; x_initial=zeros(14,1); trang 80 %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]= ode45(@ptvp_xe_chiu_uon,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]= ode45(@ptvp_xe_chiu_xoan,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]= ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_smooth_state, x_initial_bending); [T_smooth_state Y_smooth_state_torsion]= ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_smooth_state, x_initial_torsion); ddtheta_bend=[Y_bump_state_bending(:,10); Y_smooth_state_bending(:,10)]; ddtheta_tor=[Y_bump_state_torsion(:,10); Y_smooth_state_torsion(:,10)]; time=[T_bump_state;T_smooth_state]; axes('FontSize',16); hold on plot(time,ddtheta_bend,'-b','LineWidth',4); plot(time,ddtheta_tor,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('t, s'); ylabel('Accel of x, rad/s^2 '); title('GIA TOC DDTHETA THEO THOI GIAN'); %gia tốc x theo tần số lực kích động trang 81 clear clc;close all d1=1;d2=0.25; vx=linspace(0.2,10,30)/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1./vx; t_step=0.001; t_smooth_end=10; for i=1:length(vx) t_bump_state=t_start:t_step:t_end(i)-t_step; t_smooth_state=t_end(i):t_step:t_smooth_end; x_initial=zeros(14,1); x_initial(15,1)=vx(i); %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]=ode45( @ptvp_xe_chiu_uon_1,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]=ode45( @ptvp_xe_chiu_xoan_1,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]= ode45(@ptvp_xe_chay_tren_mat_duong_phang_1,t_smooth_state, x_initial_bending); [T_smooth_state Y_smooth_state_torsion]= ode45(@ptvp_xe_chay_tren_mat_duong_phang_1,t_smooth_state, x_initial_torsion); ddx_bend=[Y_bump_state_bending(:,8);Y_smooth_state_bending(:,8)]; max_ddx_bend(i)=max(ddx_bend); ddx_tor=[Y_bump_state_torsion(:,8);Y_smooth_state_torsion(:,8)]; trang 82 max_ddx_tor(i)=max(ddx_tor); end axes('FontSize',16); hold on plot(3.6*vx/(2*d1),max_ddx_bend,'-b','LineWidth',4); plot(3.6*vx/(2*d1),max_ddx_tor,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('f, Hz'); ylabel('Accel magnitude of x, m/s^2 '); title('GIA TOC DDX THEO TAN SO LUC KICH DONG'); %gia tốc theo tần số lực kích động clear clc;close all d1=1;d2=0.25; vx=linspace(0.2,10,30)/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1./vx; t_step=0.001; t_smooth_end=10; for i=1:length(vx) t_bump_state=t_start:t_step:t_end(i)-t_step; t_smooth_state=t_end(i):t_step:t_smooth_end; x_initial=zeros(14,1); x_initial(15,1)=vx(i); %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]= ode45(@ptvp_xe_chiu_uon_1,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]= ode45(@ptvp_xe_chiu_xoan_1,t_bump_state,x_initial); trang 83 %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]= ode45(@ptvp_xe_chay_tren_mat_duong_phang_1,t_smooth_state, x_initial_bending); [T_smooth_state Y_smooth_state_torsion]= ode45(@ptvp_xe_chay_tren_mat_duong_phang_1,t_smooth_state, x_initial_torsion); ddx_bend=[Y_bump_state_bending(:,9);Y_smooth_state_bending(:,9)]; max_ddx_bend(i)=max(ddx_bend); ddx_tor=[Y_bump_state_torsion(:,9);Y_smooth_state_torsion(:,9)]; max_ddx_tor(i)=max(ddx_tor); end % time=[T_bump_state;T_smooth_state]; axes('FontSize',16); hold on plot(3.6*vx/(2*d1),max_ddx_bend,'-b','LineWidth',4); plot(3.6*vx/(2*d1),max_ddx_tor,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('f, Hz'); ylabel('Accel magnitude of phi, m/s^2 '); title('GIA TOC DDPHI THEO TAN SO LUC KICH DONG'); %gia tốc theo tần số lực kích động clear clc;close all d1=1;d2=0.25; vx=linspace(0.2,10,30)/3.6; %doi tu km/h sang m/s trang 84 t_start=0; t_end=0.5*d1./vx; t_step=0.001; t_smooth_end=10; for i=1:length(vx) t_bump_state=t_start:t_step:t_end(i)-t_step; t_smooth_state=t_end(i):t_step:t_smooth_end; x_initial=zeros(14,1); x_initial(15,1)=vx(i); %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]= ode45(@ptvp_xe_chiu_uon_1,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]= ode45(@ptvp_xe_chiu_xoan_1,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]= ode45(@ptvp_xe_chay_tren_mat_duong_phang_1,t_smooth_state, x_initial_bending); [T_smooth_state Y_smooth_state_torsion]= ode45(@ptvp_xe_chay_tren_mat_duong_phang_1,t_smooth_state, x_initial_torsion); ddx_bend=[Y_bump_state_bending(:,10);Y_smooth_state_bending(:,10)]; max_ddx_bend(i)=max(ddx_bend); ddx_tor=[Y_bump_state_torsion(:,10);Y_smooth_state_torsion(:,10)]; max_ddx_tor(i)=max(ddx_tor); end % time=[T_bump_state;T_smooth_state]; axes('FontSize',16); trang 85 hold on plot(3.6*vx/(2*d1),max_ddx_bend,'-b','LineWidth',4); plot(3.6*vx/(2*d1),max_ddx_tor,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('f, Hz'); ylabel('Accel magnitude of theta, m/s^2 '); title('GIA TOC DDTHETA THEO TAN SO LUC KICH DONG'); %tải trọng động theo thời gian clear all clc;close all m=6670; fn=1.45; %chon gia tri lon nhat d1=1;d2=0.25; vx=2*d1*fn/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1/vx; t_step=0.001; t_smooth_end=2; t_bump_state=t_start:t_step:t_end-t_step; %gan het bien dang t_smooth_state=t_end:t_step:t_smooth_end; x_initial=zeros(14,1); %============BUMP_STATE================== [T_bump_state Y_bump_state_bending]=ode45(@ptvp_xe_chiu_uon,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]=ode45(@ptvp_xe_chiu_xoan,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); trang 86 [T_smooth_state Y_smooth_state_bending]=ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_s mooth_state,x_initial_bending); [T_smooth_state Y_smooth_state_torsion]=ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_sm ooth_state,x_initial_torsion); ddx_bend=[Y_bump_state_bending(:,8);Y_smooth_state_bending(:,8)]; ddx_tor=[Y_bump_state_torsion(:,8);Y_smooth_state_torsion(:,8)]; time=[T_bump_state;T_smooth_state]; f_dynamic_tor=ddx_tor*m; f_dynamic_bend=ddx_bend*m; axes('FontSize',16); hold on plot(time,f_dynamic_bend/1000,'-b','LineWidth',4); plot(time,f_dynamic_tor/1000,' r','LineWidth',4); grid on h=legend('bending','torsion',1); xlabel('t, s'); ylabel('F Dy, KN '); title('DO THI TAI TRONG DONG THEO THOI GIAN'); %%%%%%%%%%%%%%% tinh bien va rong temp1=max(abs(f_dynamic_tor)); temp2=max(abs(f_dynamic_bend)); if temp1>temp2 biendo=temp1; flag=[1;0]; else biendo=temp2; trang 87 flag=[0;1]; end f_dynamic=[f_dynamic_tor,f_dynamic_bend]; f_can_tim=f_dynamic*flag; vitri=find(abs(f_can_tim)==biendo,1);i=vitri;j=vitri; while f_can_tim(i)*f_can_tim(i-1)>0 i=i-1; end while f_can_tim(j)*f_can_tim(j+1)>0 j=j+1; end do_rong_xung=time(j)-time(i) F_dynamic_max=biendo F_static=m*9.81; ti_so_taitrongdong=F_dynamic_max/F_static*100 %Giá trị tuyệt đối tải trọng động clear all clc;close all m=6670; fn=1.45; d1=1;d2=0.25; vx=2*d1*fn/3.6; %doi tu km/h sang m/s t_start=0; t_end=0.5*d1/vx; t_step=0.001; t_smooth_end=2; t_bump_state=t_start:t_step:t_end-t_step; %gan het bien dang t_smooth_state=t_end:t_step:t_smooth_end; x_initial=zeros(14,1); %============BUMP_STATE================== trang 88 [T_bump_state Y_bump_state_bending]=ode45(@ptvp_xe_chiu_uon,t_bump_state,x_initial); [T_bump_state Y_bump_state_torsion]=ode45(@ptvp_xe_chiu_xoan,t_bump_state,x_initial); %============SMOOTH_STATE================== x_initial_bending=Y_bump_state_bending(end,:); x_initial_torsion=Y_bump_state_torsion(end,:); [T_smooth_state Y_smooth_state_bending]=ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_s mooth_state,x_initial_bending); [T_smooth_state Y_smooth_state_torsion]=ode45(@ptvp_xe_chay_tren_mat_duong_phang,t_sm ooth_state,x_initial_torsion); ddx_bend=[Y_bump_state_bending(:,8);Y_smooth_state_bending(:,8)]; ddx_tor=[Y_bump_state_torsion(:,8);Y_smooth_state_torsion(:,8)]; time=[T_bump_state;T_smooth_state]; f_dynamic_tor=ddx_tor*m; f_dynamic_bend=ddx_bend*m; F_static=m*9.81; axes('FontSize',16); hold on plot(time,-(f_dynamic_bend-F_static)/1000,'-b','LineWidth',4); plot(time,-(f_dynamic_tor-F_static)/1000,' r','LineWidth',4); set(gca,'xlim',[0 2],'ylim',[0 75]); grid on h=legend('bending','torsion',1); xlabel('t, s'); ylabel('F Dy, KN '); trang 89 title('DO THI GIA TRI TUYET DOI TAI TRONG DONG THEO THOI GIAN'); trang 90 S K L 0 ... Lagrange Kết nghiên cứu đề tài sở để đánh giá mức độ ảnh hƣởng tải trọng động đến an toàn chuyển động xe di chuyển đƣờng Xác định đƣợc vận tốc tới hạn độ cao mấp mô làm cho xe chuyển động an toàn. .. nhân gây tải trọng động trực tiếp tác dụng lên thân xe Có nhiều nghiên cứu liên quan tải trọng động tác dụng lên thân xe, tính tốn góc xoay chuyển vị thân xe, tính tốn tải trọng tĩnh động tác... nghiên cứu ảnh hƣởng tải trọng động đến an toàn chuyển động xe, dựa vào kết tính tốn ta xác định đƣợc điều kiện để xe chạy an toàn đƣờng, xác định đƣợc vận tốc tới hạn độ cao mấp mô nguy hiễm mà xe