Tổng quan nghiên cứu

Hệ thống giao thông công cộng tại các đô thị lớn đang phát triển mạnh mẽ nhằm giải quyết tình trạng ùn tắc giao thông. Tại Thành phố Hồ Chí Minh, mạng lưới xe buýt đã mở rộng với 152 tuyến xe và hơn 3.200 phương tiện hoạt động liên tục trên 148 tuyến đường. Theo Quyết định số 2545/QĐ-UBND, thành phố đã phê duyệt đề án đầu tư và thay mới 1.680 xe buýt giai đoạn 2014 - 2017 nhằm hiện đại hóa đội ngũ phương tiện công cộng. Tuy nhiên, sự gia tăng mật độ phương tiện cũng đi kèm áp lực lớn đối với đội ngũ vận hành. Người lái xe buýt đô thị thường xuyên phải làm việc từ 10 giờ mỗi ngày, 6 đến 7 ngày mỗi tuần, với quãng đường vận hành trung bình đạt 4.472,52 km mỗi tháng trong điều kiện giao thông đông đúc và đường sá nhiều gờ mấp mô.

Tình trạng phơi nhiễm rung động kéo dài với tần số bất lợi gây ảnh hưởng nghiêm trọng đến cột sống, làm tăng mức độ mệt mỏi và tiềm ẩn nguy cơ gây mất an toàn giao thông. Mục tiêu của nghiên cứu là xây dựng mô hình toán học và phương pháp tính toán dao động của ghế lái trên xe buýt đô thị SAMCO BGQ1. Nghiên cứu tập trung giải quyết bài toán tối ưu hóa hệ số đàn hồi và hệ số giảm chấn của hệ treo ghế lái dựa trên phương pháp Lagrange và công cụ mô phỏng số. Đề tài được triển khai thực hiện từ tháng 1 năm 2015 đến tháng 6 năm 2015 tại Trường Đại học Bách Khoa Thành phố Hồ Chí Minh. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở lý thuyết chuẩn mực để giảm hơn 40% gia tốc dao động thẳng đứng, đưa chỉ số rung động toàn thân về ngưỡng an toàn dưới 0,315 m/s² theo tiêu chuẩn hiện hành, từ đó cải thiện rõ rệt điều kiện làm việc cho tài xế xe buýt.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng lý thuyết cơ học dao động phi tuyến và phương pháp giải tích động lực học Lagrange để thiết lập hệ phương trình vi phân chuyển động bậc hai cho hệ thống dao động đa bậc tự do. Hệ dao động được mô tả thông qua các ma trận khối lượng, ma trận độ cản tiêu tán Rayleigh và ma trận độ cứng đàn hồi. Đánh giá độ êm dịu và tác động sinh học của dao động lên cơ thể người lái dựa trên tiêu chuẩn quốc tế ISO 2631-1:1997, tiêu chuẩn rung động Việt Nam TCVN 6964-1:2001 và tiêu chuẩn hướng dẫn kỹ thuật VDI 2057:2002 của Đức.

Năm khái niệm cốt lõi được sử dụng xuyên suốt bao gồm: mô hình động lực học một phần tư xe với 3 bậc tự do; gia tốc bình phương trung bình trọng số; dải tần số nguy hiểm đối với cơ thể người từ 4 đến 8 Hz theo phương thẳng đứng; hệ số độ cứng lò xo k; và hệ số cản giảm chấn c của hệ thống treo. Các chỉ tiêu đánh giá độ êm dịu chuyển động dựa trên phản ứng sinh lý học, trong đó mức gia tốc thẳng đứng vượt ngưỡng 0,63 m/s² được xác định là gây nguy hại trực tiếp đến cấu trúc cột sống của người lái khi tiếp xúc liên tục từ 4 đến 8 giờ mỗi ngày.

Phương pháp nghiên cứu

Nguồn dữ liệu đầu vào được thu thập trực tiếp từ các thông số kết cấu thực tế của xe khách đô thị 50 chỗ SAMCO BGQ1 do Tổng Công ty Cơ khí GTVT Sài Gòn sản xuất. Các thông số chính gồm: khối lượng toàn bộ 8.500 kg, khối lượng phân bố lên cầu trước là 2.935 kg, khối lượng một phần tư thân xe quy đổi đạt 1.467,5 kg, khối lượng phần không được treo là 124 kg, khối lượng ghế và người lái là 81 kg. Độ cứng của lốp xe xác định theo công thức bán thực nghiệm đạt 564.271,2 N/m và độ cứng nhíp lá cầu trước đạt 65.544 N/m.

Về phương pháp chọn mẫu, nghiên cứu áp dụng phương pháp chọn mẫu phương tiện đại diện cho dòng xe buýt đô thị kích thước trung bình phổ biến nhất tại Việt Nam. Kích động mặt đường được xây dựng theo mô hình hàm bán bình phương hình sin với chiều dài bước sóng 1 m, chiều cao mấp mô 0,25 m ở vận tốc vận hành 40 km/h. Phương pháp phân tích số trên nền tảng phần mềm MATLAB kết hợp Simulink được lựa chọn vì khả năng giải chính xác các phương trình vi phân phi tuyến bậc cao, giúp rút ngắn hơn 80% thời gian khảo sát tham số và tiết kiệm chi phí thử nghiệm thực địa. Toàn bộ quá trình tính toán, mô phỏng và phân tích tối ưu được thực hiện chặt chẽ trong thời gian 5 tháng, từ ngày 19 tháng 01 năm 2015 đến ngày 14 tháng 06 năm 2015.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Khảo sát trên mô hình nguyên bản chưa tối ưu với độ cứng ban đầu của ghế k = 15.000 N/m và hệ số cản c = 150 N.s/m cho thấy biên độ dao động của ghế lái bị khuếch đại đáng kể khi xe di chuyển qua mấp mô ở dải tốc độ 20 đến 60 km/h. Tại vận tốc 40 km/h, gia tốc dao động lớn nhất truyền lên người lái vượt ngưỡng 0,63 m/s², nằm trong khu vực cảnh báo gây mệt mỏi nghiêm trọng và thoái hóa cột sống theo tiêu chuẩn TCVN 6964-1:2001.

Thông qua quá trình quét đa biến độc lập trên phần mềm mô phỏng, nghiên cứu đã xác định được vùng thông số tối ưu cho hệ treo ghế lái:

  • Độ cứng tối ưu của hệ treo ghế lái được xác định ở dải từ 8.500 N/m đến 10.000 N/m. Khi độ cứng giảm từ 15.000 N/m xuống mức tối ưu này, biên độ chuyển vị đỉnh của ghế giảm 35,4% và lực động truyền lên cơ thể người giảm rõ rệt.
  • Hệ số giảm chấn tối ưu của ghế được xác định trong khoảng từ 350 N.s/m đến 450 N.s/m, tăng gấp hơn 2,3 lần so với giá trị ban đầu là 150 N.s/m. Hệ số này giúp tăng tốc độ tiêu tán năng lượng dao động lên 40,8%, dập tắt hoàn toàn rung lắc sau 0,75 giây thay vì kéo dài hơn 1,8 giây như hệ thống cũ.
  • Chỉ số gia tốc bình phương trung bình trọng số của ghế sau khi tối ưu hóa đồng thời cả k và c giảm hơn 42,5%, đưa mức rung động từ trạng thái khá không thoải mái xuống vùng an toàn và dễ chịu dưới 0,315 m/s².

Thảo luận kết quả

Nguyên nhân của sự cải thiện vượt bậc này xuất phát từ việc điều chỉnh độ cứng k và hệ số cản c giúp dịch chuyển tần số dao động riêng của cụm ghế và người lái ra khỏi dải tần số cộng hưởng nguy hiểm của cơ thể người từ 4 đến 8 Hz. Hệ số cản được tăng cường giúp hấp thụ hiệu quả các xung lực đột ngột truyền từ mặt đường qua hệ thống nhíp lá cầu trước có độ cứng 65.544 N/m. Kết quả này tương đồng và khẳng định tính đúng đắn khi so sánh với các nghiên cứu động lực học 4 bậc tự do của Ashraf Salah và nghiên cứu mô hình cơ thể người 7 bậc tự do của Stefan Tabacu.

Dữ liệu mô phỏng được biểu diễn trực quan và sinh động thông qua hệ thống đồ thị dạng sóng 2D biểu diễn đáp ứng chuyển vị và vận tốc theo thời gian thực từ 0 đến 2 giây. Bên cạnh đó, các biểu đồ không gian 3D thể hiện tương quan bề mặt giữa vận tốc xe từ 10 đến 60 km/h, độ cứng k và gia tốc cực đại đã minh chứng rõ nét điểm cực tiểu của gia tốc tại vùng thông số tối ưu. Bảng tổng hợp so sánh các giá trị đỉnh trước và sau khi tối ưu hóa cung cấp bằng chứng định lượng rõ ràng cho việc thiết kế chế tạo hệ thống treo thực tế.

Đề xuất và khuyến nghị

Nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật và chính sách mang tính khả thi cao nhằm cải thiện điều kiện làm việc của tài xế xe buýt:

  1. Chế tạo cụm treo giảm chấn độc lập cho ghế lái xe buýt đô thị: Các doanh nghiệp sản xuất ô tô như Tổng Công ty Cơ khí GTVT Sài Gòn SAMCO cần ứng dụng bộ thông số độ cứng k khoảng 9.000 N/m và hệ số giảm chấn c khoảng 400 N.s/m vào quy trình dập khuôn và lắp ráp ghế lái. Mục tiêu giảm tối thiểu 40% gia tốc dao động thẳng đứng trong vòng 6 tháng triển khai thử nghiệm.
  2. Ban hành quy chuẩn kỹ thuật bắt buộc về độ êm dịu ghế lái: Cục Đăng kiểm Việt Nam và Sở Giao thông Vận tải cần xây dựng khung tiêu chuẩn kiểm định định kỳ độ êm dịu của ghế tài xế xe công cộng, quy định gia tốc rung động hiệu dụng không được vượt quá 0,315 m/s² theo TCVN 6964-1:2001, hoàn thành lộ trình ban hành trong vòng 12 tháng.
  3. Hoán cải và nâng cấp hệ thống ghế ngồi trên đội xe hiện hữu: Các hợp tác xã và doanh nghiệp vận tải hành khách công cộng tại Thành phố Hồ Chí Minh cần triển khai thay thế cụm giảm chấn thủy lực hoặc lắp đệm khí nén bổ trợ cho hơn 3.200 xe buýt đang lưu hành, phấn đấu đạt tỷ lệ nâng cấp trên 85% phương tiện trong thời gian từ 18 đến 24 tháng.
  4. Phát triển hệ thống treo bán tích cực thông minh: Các viện nghiên cứu và trường đại học chuyên ngành cần tiếp tục nghiên cứu tích hợp giảm chấn từ biến từ tính với thời gian phản hồi dưới 15 mili-giây, cho phép tự động thay đổi độ cản thích ứng theo vận tốc vận hành từ 20 đến 60 km/h, thực hiện chuyển giao công nghệ trong giai đoạn 2 đến 3 năm tới.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo giá trị cho nhiều nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế và chuyên gia R&D ngành ô tô: Tiếp cận phương pháp thiết lập phương trình Lagrange và quy trình mô phỏng số trên MATLAB để tính toán tối ưu hóa hệ thống treo cho các dòng xe thương mại mới.
  • Doanh nghiệp sản xuất, lắp ráp xe khách và xe buýt: Khai thác bộ thông số kỹ thuật tối ưu để trực tiếp thiết kế, chế tạo các mẫu ghế lái giảm chấn chất lượng cao, giúp tiết kiệm từ 30% đến 50% chi phí so với nhập khẩu nguyên cụm từ nước ngoài.
  • Cơ quan quản lý giao thông và tổ chức bảo hộ lao động: Sử dụng cơ sở đánh giá thực nghiệm và tiêu chuẩn ISO 2631-1:1997 để xây dựng chính sách chăm sóc sức khỏe nghề nghiệp, giảm thiểu tai nạn lao động cho đội ngũ tài xế.
  • Giảng viên, học viên cao học và sinh viên kỹ thuật giao thông: Sử dụng đề tài làm tài liệu giảng dạy và học tập điển hình về giải quyết bài toán dao động phi tuyến và cơ học ứng dụng trong kỹ thuật ô tô hiện đại.

Câu hỏi thường gặp

Tại sao dải tần số từ 4 đến 8 Hz lại gây nguy hiểm nhất cho tài xế xe buýt? Tần số từ 4 đến 8 Hz trùng với tần số dao động tự nhiên của các cơ quan nội tạng và cột sống con người theo tiêu chuẩn VDI 2057:2002. Khi xe vận hành qua đường mấp mô tạo ra kích động trong dải này, hiện tượng cộng hưởng sinh học sẽ xảy ra, làm gia tăng biên độ chấn động và gây tổn thương cột sống nếu lái xe liên tục trên 4 giờ mỗi ngày.

Mô hình một phần tư xe có phản ánh chính xác dao động thực tế của xe buýt không? Mô hình một phần tư xe với 3 bậc tự do phản ánh chính xác đến hơn 85% các dao động thẳng đứng tác động trực tiếp lên người lái. Do dao động thẳng đứng chiếm tỷ trọng năng lượng lớn nhất gây mệt mỏi so với dao động lắc ngang hay lắc dọc, việc sử dụng mô hình này giúp đơn giản hóa bài toán vi phân nhưng vẫn đảm bảo độ tin cậy kỹ thuật cao.

Việc tối ưu hóa độ cứng k và hệ số cản c mang lại lợi ích kinh tế nào? Việc xác định chính xác cặp thông số tối ưu giúp các xưởng chế tạo trong nước sản xuất trực tiếp hệ thống treo ghế lái phù hợp địa hình Việt Nam, giúp giảm giá thành sản phẩm từ 30% đến 50% so với nhập khẩu nguyên bộ, đồng thời giảm thiểu chi phí bồi thường bệnh nghề nghiệp và tai nạn giao thông.

Phần mềm MATLAB và Simulink có độ tin cậy ra sao trong bài toán dao động ô tô? MATLAB và Simulink là công cụ tiêu chuẩn quốc tế trong ngành kỹ thuật ô tô, cho phép giải các hệ phương trình vi phân Lagrange với độ chính xác cao và sai số thuật toán dưới 1%. Mô phỏng số cho phép khảo sát liên tục hàng trăm biến thể tốc độ và tải trọng mà không tốn kém chi phí thử nghiệm phá hủy.

Tiêu chuẩn TCVN 6964-1:2001 quy định ngưỡng rung động an toàn là bao nhiêu? TCVN 6964-1:2001 quy định mức gia tốc rung động bình phương trung bình dưới 0,315 m/s² là trạng thái thoải mái không gây khó chịu. Mức từ 0,5 đến 1,0 m/s² gây cảm giác rất khó chịu, và trên 2,0 m/s² là mức cực kỳ nguy hiểm, không được phép duy trì trong môi trường lao động kéo dài.

Kết luận

  • Luận văn đã xây dựng thành công mô hình toán học động lực học một phần tư xe với 3 bậc tự do bằng phương pháp Lagrange, phản ánh chính xác dao động thẳng đứng của tài xế xe buýt.
  • Ứng dụng thành công phần mềm mô phỏng MATLAB để khảo sát toàn diện ảnh hưởng của độ cứng k và hệ số giảm chấn c trên xe khách SAMCO BGQ1.
  • Xác định chính xác thông số tối ưu của ghế lái với độ cứng k đạt khoảng 9.000 N/m và hệ số cản c đạt khoảng 400 N.s/m, giúp giảm hơn 42,5% gia tốc dao động thẳng đứng.
  • Đưa mức rung động toàn thân của người lái về ngưỡng an toàn dưới 0,315 m/s² theo tiêu chuẩn quốc tế ISO 2631-1:1997 và TCVN 6964-1:2001.
  • Đóng góp cơ sở lý thuyết và thực tiễn vững chắc phục vụ chế tạo nội địa hóa ghế lái giảm chấn cho các phương tiện vận tải công cộng tại Việt Nam.

Lộ trình triển khai ứng dụng thực tế dự kiến sẽ được thử nghiệm trên các tuyến buýt nội đô trong thời gian từ 6 đến 12 tháng tới. Các viện nghiên cứu, doanh nghiệp sản xuất ô tô và đơn vị vận tải hãy liên hệ hợp tác chuyển giao công nghệ để chung tay nâng cao sức khỏe người lao động và xây dựng hệ thống giao thông công cộng an toàn, văn minh.