Tổng quan nghiên cứu

Hệ thống giao thông công cộng tại các đô thị lớn của Việt Nam đang chứng kiến sự phát triển mạnh mẽ về quy mô nhằm giải quyết bài toán ùn tắc. Tại Hà Nội, mạng lưới xe buýt đạt trên 1.100 phương tiện hoạt động trên 140 tuyến, đáp ứng khoảng 9% nhu cầu đi lại của người dân. Trong khi đó, Thành phố Hồ Chí Minh vận hành hơn 2.100 xe buýt trên 126 tuyến, phục vụ 9,2% nhu cầu giao thông đô thị và hướng tới mục tiêu 15% đến 20% trong giai đoạn tới. Sự gia tăng này kéo theo lực lượng lao động hơn 6.400 tài xế làm việc luân phiên 2 ca mỗi ngày, với thời gian ngồi lái liên tục từ 8 đến 10 giờ trong điều kiện chịu rung động toàn thân kéo dài.

Rung động truyền từ mặt đường qua sàn xe lên ghế lái trong dải tần số 0,5 đến 100 Hz tác động trực tiếp đến cột sống và hệ thần kinh người điều khiển. Nguy hiểm nhất là vùng tần số cộng hưởng 2 đến 5 Hz theo phương đứng, nơi cơ thể người bị khuếch đại biên độ rung động mạnh nhất, dẫn đến các bệnh nghề nghiệp nghiêm trọng như thoái hóa đốt sống lưng, thoát vị đĩa đệm và suy giảm phản xạ vận hành. Nghiên cứu này tập trung phân tích dao động ghế ngồi người lái xe buýt thông qua việc xây dựng mô hình toán học động lực học 3 bậc tự do (3DOF), kết hợp giữa mô hình 1/4 xe và hệ thống treo ghế gồm bầu hơi đàn hồi phi tuyến, giảm chấn thủy lực phi tuyến cùng cơ cấu dẫn hướng chữ X.

Mục tiêu cốt lõi của luận văn là đánh giá toàn diện độ êm dịu, mức độ phơi nhiễm rung động đối với sức khỏe tài xế và khả năng cách ly rung động (chỉ số SEAT) dưới kích thích điều hòa và ngẫu nhiên mặt đường cấp B, C theo các tiêu chuẩn quốc tế ISO 2631-1:1997 và Chỉ thị Châu Âu 2002/44/EC.

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 động lực học hệ nhiều vật kết hợp với nhiệt động học khí nén và cơ học lưu chất để mô tả chính xác các phản lực phi tuyến trong hệ thống treo ghế. Mô hình tổng thể tích hợp 3 bậc tự do bao gồm: khối lượng phần không được treo của xe, khối lượng thân xe và khối lượng người lái (quy ước 70 kg) đặt trên mặt ghế.

Các khái niệm then chốt được áp dụng trong đánh giá bao gồm:

  • Khả năng cách ly rung động của ghế (SEAT): Tỷ số phần trăm giữa gia tốc hiệu dụng trên mặt ghế so với gia tốc trên mặt sàn xe.
  • Gia tốc trung bình bình phương có trọng số tần số ($a_w$): Đại lượng đặc trưng cho độ êm dịu theo tiêu chuẩn ISO 2631-1:1997, sử dụng hàm trọng số $w_k$ cho dao động theo phương đứng.
  • Giá trị phơi nhiễm rung động hàng ngày tương đương 8 giờ [$A(8)$]: Chỉ số đánh giá tổn hại sức khỏe nghề nghiệp theo tiêu chuẩn 2002/44/EC.
  • Mật độ phổ công suất mặt đường (PSD): Đặc trưng biên dạng mấp mô ngẫu nhiên theo tiêu chuẩn ISO 8608.

Bầu hơi được mô hình hóa dựa trên phương trình trạng thái khí lý tưởng và phương trình trao đổi năng lượng, tính đến sự thay đổi áp suất thiết lập ban đầu trong dải 2 đến 6 bar. Giảm chấn thủy lực ống đơn được thiết lập dựa trên phương trình lưu lượng dòng chảy nhớt qua lỗ tiết lưu có đường kính thay đổi từ 1,5 đến 3,0 mm, phản ánh lực cản nhớt phi tuyến theo vận tốc biến dạng.

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

Nguồn dữ liệu đầu vào sử dụng biên dạng kích thích mặt đường dạng điều hòa trong dải tần số 0,5 đến 20 Hz và mặt đường mấp mô ngẫu nhiên cấp B (mật độ phổ chuẩn từ 2 đến 8 $\times 10^{-6}\text{ m}^3/\text{rad}$) và cấp C (từ 8 đến 32 $\times 10^{-6}\text{ m}^3/\text{rad}$) theo tiêu chuẩn ISO 8608 ở các dải tốc độ xe chạy từ 20 đến 60 km/h.

Nghiên cứu khảo sát 12 tổ hợp thông số vận hành thực nghiệm bao gồm 3 mức áp suất khí ban đầu của bầu hơi (2 bar, 4 bar, 6 bar) kết hợp với 4 mức đường kính lỗ tiết lưu giảm chấn (1,5 mm; 2,0 mm; 2,5 mm; 3,0 mm). Cỡ mẫu 12 tổ hợp này được lựa chọn có chủ đích nhằm bao phủ toàn bộ phạm vi tải trọng thực tế của người lái xe (từ 50 kg đến 90 kg) và các mức cài đặt độ cứng giảm chấn trên các dòng xe buýt thương mại phổ biến.

Phương pháp phân tích số (Numerical Simulation) trên phần mềm MATLAB được lựa chọn thay vì phương pháp tuyến tính hóa truyền thống. Lý do là vì giải thuật tích phân số phi tuyến cho phép giải quyết chính xác sự biến thiên nhiệt độ, áp suất tức thời trong bầu hơi và hiệu ứng phi tuyến tính khi đóng mở van giảm xóc, giúp loại bỏ sai số lên đến 20% so với các mô hình tuyến tính đơn giản. Timeline nghiên cứu được thực hiện trong 10 tháng, bao gồm đo đạc kích thước hình học cơ cấu chữ X thực tế, lập trình thuật toán giải hệ phương trình vi phân và xử lý dữ liệu đáp ứng tần số.

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

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

Phân tích số liệu từ mô phỏng mang lại 4 phát hiện quan trọng về động lực học hệ thống treo ghế:

  1. Hiệu quả cách ly rung động phụ thuộc mạnh vào áp suất bầu hơi: Ở dải tần số kích thích 2 đến 5 Hz, cấu hình áp suất bầu hơi ban đầu 2 bar kết hợp đường kính lỗ tiết lưu 2,0 mm cho chỉ số SEAT thấp nhất, chỉ đạt khoảng 82% đến 85%, nghĩa là ghế đã giảm được 15% đến 18% năng lượng rung động từ sàn xe. Ngược lại, khi tăng áp suất lên 6 bar, chỉ số SEAT cực đại tăng vọt lên mức 125% đến 132%, làm khuếch đại biên độ dao động truyền vào cơ thể tài xế.
  2. Ảnh hưởng của cấp mặt đường và tốc độ xe: Khi xe di chuyển trên đường cấp B ở tốc độ 40 km/h, gia tốc hiệu dụng $a_w$ duy trì trong khoảng 0,38 đến 0,52 m/s², nằm trong ngưỡng cho phép về độ thoải mái. Tuy nhiên, khi chuyển sang mặt đường cấp C ở tốc độ 50 km/h, giá trị $a_w$ tăng hơn 80%, chạm mức 0,85 đến 0,98 m/s², gây cảm giác rất không thoải mái cho tài xế.
  3. Mức độ phơi nhiễm rung động theo thời gian: Giá trị phơi nhiễm $A(8)$ sau 8 giờ làm việc với cấu hình tối ưu đạt mức 0,46 m/s², nằm dưới ngưỡng cảnh báo hành động 0,50 m/s² của Chỉ thị 2002/44/EC. Khi lỗ tiết lưu giảm chấn mở rộng đến 3,0 mm, giá trị $A(8)$ vượt ngưỡng 0,68 m/s², làm tăng nguy cơ tổn thương cột sống lên gấp 1,5 lần.
  4. Đặc tính phi tuyến của cơ cấu dẫn hướng chữ X: Góc nghiêng của hai thanh dẫn hướng chữ X làm thay đổi tỷ số truyền chuyển vị và truyền lực lên bầu hơi từ 10% đến 15% tùy thuộc vào hành trình nhún, điều mà các mô hình 1 bậc tự do dạng thanh trượt đứng thẳng không thể phản ánh được.

Thảo luận kết quả

Nguyên nhân của hiện tượng trên bắt nguồn từ sự tương quan giữa độ cứng động học của khối khí nén và khả năng tiêu tán năng lượng của dòng dầu thủy lực. Khi áp suất thiết lập ban đầu quá cao (6 bar), thể tích khí bị nén chặt làm độ cứng lò xo khí tăng cao, biến hệ thống treo thành một liên kết cứng truyền trực tiếp xung lực tần số thấp lên người ngồi. Về phía giảm chấn, đường kính lỗ tiết lưu quá nhỏ (1,5 mm) gây hiện tượng nghẽn dòng ở vận tốc dao động cao, trong khi đường kính quá lớn (3,0 mm) lại không đủ lực cản để dập tắt dao động cộng hưởng.

So với các nghiên cứu trước đây về ghế đệm mút bị động thông thường vốn chỉ giảm được khoảng 8% đến 10% biên độ gia tốc đỉnh, mô hình hệ thống treo khí nén tích hợp cơ cấu chữ X được hiệu chỉnh thông số tối ưu giúp tăng hiệu quả dập tắt rung động lên hơn 35%.

Để trực quan hóa các mối quan hệ phức tạp này, dữ liệu nghiên cứu có thể được trình bày thông qua biểu đồ bề mặt 3D thể hiện sự phụ thuộc của chỉ số SEAT theo hai trục thông số áp suất và độ mở lỗ tiết lưu, kết hợp bảng phân loại ma trận rủi ro sức khỏe tương ứng với các cấp mặt đường và tốc độ vận hành của phương tiện.

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

Dựa trên kết quả nghiên cứu thực nghiệm mô phỏng, 4 giải pháp cụ thể được đề xuất nhằm bảo vệ sức khỏe người lái và nâng cao hiệu quả khai thác xe buýt đô thị:

  • Hiệu chuẩn áp suất làm việc của bầu hơi theo thể trạng tài xế: Doanh nghiệp vận tải xe buýt cần trang bị đồng hồ đo áp suất và hướng dẫn tài xế tự điều chỉnh áp suất bầu hơi trong dải tối ưu từ 2,0 đến 4,0 bar tương ứng với khối lượng cơ thể từ 55 đến 80 kg, nhằm duy trì chỉ số cách ly rung động SEAT luôn dưới ngưỡng 88%. Thời gian thực hiện kiểm tra định kỳ 1 tháng/lần.
  • Tùy chỉnh độ mở lỗ tiết lưu giảm chấn thủy lực: Bộ phận kỹ thuật bảo dưỡng tại các xí nghiệp xe buýt tiến hành cố định núm xoay điều chỉnh van tiết lưu ở mức đường kính 1,8 đến 2,0 mm. Giải pháp này giúp cắt giảm 25% đến 30% gia tốc đỉnh khi xe vận hành qua các cung đường gồ ghề, hoàn thành triển khai trên 100% đầu xe trong vòng 6 tháng.
  • Xây dựng quy chuẩn giám sát phơi nhiễm rung động nghề nghiệp: Sở Giao thông Vận tải phối hợp cùng Sở Y tế áp dụng tiêu chuẩn TCVN 6964-1:2001 và Chỉ thị 2002/44/EC để đo lường định kỳ mức độ phơi nhiễm rung động $A(8) \le 0,5\text{ m/s}^2$ cho đội ngũ lái xe, đồng thời điều chỉnh thời gian lái xe liên tục không quá 4 giờ mỗi ca, áp dụng trong lộ trình 12 tháng tới.
  • Nâng cấp hệ thống treo ghế bán tích cực (Semi-active): Các đơn vị nghiên cứu và sản xuất phụ tùng ô tô trong nước ứng dụng mô hình toán học của luận văn để tích hợp van điện từ điều khiển liên tục độ giảm chấn theo thời gian thực, hướng tới mục tiêu giảm thêm 20% gia tốc rung động trong giai đoạn 2026 - 2030.

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

Công trình nghiên cứu mang lại giá trị lý thuyết và ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư thiết kế và phát triển sản phẩm ô tô (R&D): Tiếp cận mô hình toán học phi tuyến hoàn chỉnh mô tả bầu hơi, giảm chấn thủy lực và cơ cấu dẫn hướng chữ X, giúp rút ngắn 40% thời gian thử nghiệm và tối ưu hóa hệ thống treo ghế xe thương mại.
  • Cán bộ quản lý kỹ thuật và điều hành doanh nghiệp vận tải: Ứng dụng các bảng tra cứu thông số áp suất bầu hơi và độ mở van giảm xóc để chuẩn hóa quy trình bảo dưỡng, kéo dài tuổi thọ cụm ghế lái thêm 2 đến 3 năm và giảm thiểu tỷ lệ tài xế xin nghỉ ốm do bệnh cột sống.
  • Chuyên gia an toàn lao động và y tế giao thông: Sử dụng phương pháp tính toán gia tốc có trọng số tần số theo ISO 2631-1:1997 và chỉ số $A(8)$ theo 2002/44/EC làm bằng chứng khoa học để đề xuất chính sách bảo hiểm và khám sức khỏe định kỳ cho người lao động.
  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Cơ khí Động lực: Sử dụng mã nguồn mô phỏng trên nền tảng MATLAB và phương pháp giải tích số làm tài liệu tham khảo chất lượng cao cho các nghiên cứu chuyên sâu về động lực học dao động xe.

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

Tại sao cần phải mô hình hóa phi tuyến cho bầu hơi và giảm chấn thay vì dùng mô hình tuyến tính?

Mô hình tuyến tính bỏ qua sự thay đổi nhiệt dung, áp suất khí nén khi thể tích co giãn lớn và đặc tính tiết lưu dòng dầu qua lỗ nhỏ. Mô hình phi tuyến trong luận văn tái hiện chính xác phản lực đàn hồi và cản nhớt ở dải tần số cao lên đến 20 Hz, giảm thiểu sai số tính toán từ 15% xuống dưới 5%.

Mức áp suất ban đầu nào của bầu hơi là tối ưu nhất cho tài xế có cân nặng trung bình?

Đối với người lái có khối lượng từ 65 đến 75 kg, mức áp suất ban đầu tối ưu của bầu hơi là 2,0 đến 3,0 bar. Mức áp suất này giúp đưa tần số dao động riêng của hệ thống ghế về vùng dưới 1,5 Hz, tránh hiện tượng cộng hưởng nguy hiểm tại dải tần 2 đến 5 Hz.

Cơ cấu dẫn hướng chữ X có ảnh hưởng như thế nào đến khả năng giảm chấn của ghế?

Cơ cấu chữ X giữ cho mặt ghế chuyển vị thuần túy theo phương đứng mà không bị lắc ngang, đồng thời hình học của các thanh đòn làm biến đổi phi tuyến hành trình và vận tốc tương đối truyền vào bầu hơi và giảm chấn, tác động trực tiếp đến lực cản dập tắt dao động.

Chỉ số SEAT bao nhiêu thì đạt tiêu chuẩn an toàn cho ghế lái xe buýt?

Chỉ số SEAT phải nhỏ hơn 1 (hoặc dưới 100%) thì hệ thống treo ghế mới đạt yêu cầu cách ly rung động. Trong nghiên cứu này, cấu hình tối ưu đạt chỉ số SEAT từ 82% đến 85%, nghĩa là ghế triệt tiêu được 15% đến 18% năng lượng rung động từ sàn xe truyền lên.

Làm thế nào để kiểm soát mức phơi nhiễm rung động A(8) dưới ngưỡng 0,5 m/s² trong ca làm việc 8 giờ?

Cần kết hợp đồng bộ việc duy trì vận tốc xe dưới 45 km/h trên đường xấu cấp C, thiết lập áp suất bầu hơi ở mức 2 đến 3 bar, đường kính lỗ tiết lưu khoảng 2,0 mm và bố trí thời gian nghỉ ngơi giữa ca ít nhất 15 phút sau mỗi 2 giờ lái xe liên tục.

Kết luận

  • Xây dựng thành công mô hình toán học tích hợp 3 bậc tự do (3DOF) mô tả chi tiết đặc tính phi tuyến của bầu hơi, giảm chấn thủy lực và động học cơ cấu dẫn hướng chữ X trên hệ thống treo ghế lái xe buýt.
  • Xác lập dải thông số vận hành tối ưu với áp suất bầu hơi $P_0 = 2 \div 4\text{ bar}$ và đường kính lỗ tiết lưu giảm chấn $d_h = 2,0\text{ mm}$, giúp hạ chỉ số truyền dao động SEAT xuống dưới mức 85%.
  • Đánh giá định lượng tác động của rung động toàn thân theo tiêu chuẩn ISO 2631-1:1997 và 2002/44/EC, chứng minh giá trị phơi nhiễm $A(8)$ được kiểm soát an toàn dưới ngưỡng 0,5 m/s² khi vận hành đúng thông số.
  • Cung cấp bộ công cụ thuật toán mô phỏng hoàn chỉnh trên môi trường MATLAB, phục vụ hiệu quả cho công tác nghiên cứu, đào tạo và thiết kế kỹ thuật ô tô.
  • Định hướng phát triển tiếp theo trong 24 tháng tới là mở rộng mô hình sang hệ thống treo điều khiển bán tích cực (semi-active suspension) và tiến hành thử nghiệm đo đạc thực địa trên xe buýt chạy tuyến thực tế.

Các doanh nghiệp vận tải hành khách công cộng và đơn vị sản xuất ô tô hãy liên hệ áp dụng ngay giải pháp tối ưu hóa thông số hệ thống treo ghế lái để nâng cao an toàn giao thông và bảo vệ sức khỏe người lao động.