Tổng quan nghiên cứu

Trong cơ cấu vận tải hàng hóa đường bộ tại Việt Nam, xe tải nhẹ tải trọng dưới 1.000 kg chiếm hơn 45% tổng lưu lượng phương tiện lưu thông nội đô và liên tỉnh nhờ tính cơ động cao. Tuy nhiên, phân khúc xe này lại đối mặt với nguy cơ mất an toàn nghiêm trọng, đặc biệt là hiện tượng mất ổn định ngang và lật xe khi vào cua gấp hoặc chuyển làn ở tốc độ cao trên các cung đường đèo dốc phức tạp. Hiện nay, quy trình kiểm định an toàn kỹ thuật phương tiện cơ giới đường bộ tại Việt Nam vẫn chưa tích hợp đầy đủ các tiêu chuẩn định lượng khắt khe về giới hạn lật xe theo thời gian thực.

Công trình tập trung nghiên cứu, xây dựng mô hình toán học và phân tích động lực học độ ổn định ngang của ô tô tải nhẹ SMV T880 có chiều dài cơ sở 2.500 mm, kích thước bao 3.810 x 1.485 x 1.810 mm, khối lượng bản thân 880 kg và tải trọng cho phép 880 kg (khối lượng toàn tải 1.890 kg). Mục tiêu cốt lõi là xác lập hệ thống thông số kết cấu thực tế, thiết lập phương trình động lực học 4 bậc tự do và mô phỏng chính xác ứng xử động lực học thân xe trong không gian Matlab R2012b.

Nghiên cứu được triển khai tại Bộ môn Ôtô - Máy động lực, Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh, hoàn thành và bảo vệ thành công vào tháng 07 năm 2014. Kết quả mang lại bộ cơ sở dữ liệu định lượng chính xác về giới hạn lật, hỗ trợ thiết lập ngưỡng vận tốc tới hạn, góp phần nâng cao tính an toàn chủ động cho các dòng xe tải nhẹ đang lưu hành.

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 nền tảng động lực học hệ nhiều vật và cơ học giải tích Newton - Euler để mô tả không gian chuyển động của ô tô trong hệ quy chiếu Descartes Cxyz gắn liền với trọng tâm thân xe. Mô hình động lực học phẳng dạng 2 bánh (mô hình xe đạp cải tiến) tích hợp dao động lắc ngang được xây dựng với 4 bậc tự do chính: chuyển động tịnh tiến dọc theo trục x, chuyển động dạt ngang theo trục y, dao động lắc quanh trục dọc x (góc roll $\varphi$) và chuyển động quay vòng quanh trục đứng z (góc yaw $\psi$).

Các khái niệm then chốt được lượng hóa gồm: góc trượt ngang của bánh xe $\alpha$, lực ngang sinh ra tại vết tiếp xúc lốp - mặt đường $F_y$, độ cứng trượt ngang của lốp $C_\alpha$, độ cứng chống xoắn của hệ thống treo $k_\varphi$ (đạt giá trị 149.470 Nm/rad) và hệ số giảm chấn chống xoắn $c_\varphi$ (đạt 858,4 Nms/rad). Mô hình lốp tuyến tính kết hợp hàm phi tuyến bậc hai của tải trọng thẳng đứng $F_z$ được ứng dụng để phản ánh chính xác sự suy giảm khả năng bám ngang khi phân bố tải trọng thay đổi đột ngột giữa hai bên bánh xe trong quá trình quay vòng.

       +-------------------------------------------------------------+
       |   HỆ QUY CHIẾU THÂN XE (Cxyz) VÀ MÔ HÌNH ĐỘNG LỰC HỌC 4-DOF  |
       +-------------------------------------------------------------+
                                     |
            +------------------------+------------------------+
            |                                                 |
            v                                                 v
  [ CHUYỂN ĐỘNG TRONG MẶT PHẲNG ]                    [ DAO ĐỘNG KHÔNG GIAN ]
  - Tịnh tiến dọc: Vận tốc v_x                       - Lắc ngang quanh trục x: Góc roll φ
  - Dạt ngang: Vận tốc v_y                           - Quay thân xe quanh trục z: Góc yaw ψ
            |                                                 |
            +------------------------+------------------------+
                                     |
                                     v
       +-------------------------------------------------------------+
       |  HỆ THỐNG PHƯƠNG TRÌNH VI PHÂN NEWTON - EULER KHÔNG GIAN    |
       |  - Cân bằng lực dọc: F_x = m*(dv_x/dt - r*v_y)              |
       |  - Cân bằng lực ngang: F_y = m*(dv_y/dt + r*v_x)            |
       |  - Cân bằng mô men lắc: M_x = I_x*(dp/dt) + k_φ*φ + c_φ*p   |
       |  - Cân bằng mô men quay: M_z = I_z*(dr/dt)                  |
       +-------------------------------------------------------------+
                                     |
                                     v
       +-------------------------------------------------------------+
       |     TƯƠNG TÁC LỐP - ĐƯỜNG & HỆ THỐNG TREO BIẾN DẠNG         |
       |     - Độ cứng trượt trước C_αf, sau C_αr (hàm tải trọng F_z)|
       |     - Góc đánh lái hữu hiệu: δ_thực = δ_lái + C_δφ * φ      |
       +-------------------------------------------------------------+

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

Nguồn dữ liệu đầu vào được thu thập thông qua phương pháp đo đạc thực nghiệm trực tiếp trên xe cơ sở SMV T880 tại xưởng thực hành C3. Mẫu nghiên cứu đại diện đơn nhất được lựa chọn theo phương pháp chọn mẫu điển hình cho dòng xe tải nhẹ cabin đơn, sử dụng hệ thống treo trước độc lập MacPherson (lò xo trụ $G = 80 \times 10^3\text{ N/mm}^2$) và treo sau phụ thuộc nhíp lá ($E = 2 \times 10^8\text{ kN/m}^2$, 5 lá nhíp).

Phương pháp phân tích số trị được thực hiện trên nền tảng Matlab R2012b để giải hệ phương trình vi phân phi tuyến bậc cao mô tả phản ứng thời gian của xe. Lý do lựa chọn mô hình mô phỏng số kết hợp đo đạc thực nghiệm là nhằm tối ưu hóa độ chính xác, khảo sát an toàn các trạng thái tới hạn mà không gây nguy cơ lật đổ xe thật. Quá trình tính toán được thực hiện xuyên suốt từ ngày giao nhiệm vụ 26/04/2014 đến khi hoàn tất nghiệm thu vào ngày 23/05/2014.

Thông số kết cấu Ký hiệu Không tải (0 kg hàng) Nửa tải (440 kg hàng) Toàn tải (880 kg hàng) Đơn vị
Khối lượng toàn bộ $m$ 1.010 1.450 1.890 kg
Chiều cao trọng tâm $h$ 0,716 0,893 0,988 m
Khoảng cách trọng tâm đến cầu trước $a_1$ 1,236 1,504 1,648 m
Khoảng cách trọng tâm đến cầu sau $a_2$ 1,264 0,996 0,852 m
Mô men quán tính quanh trục dọc $I_x$ 569 679 957 $\text{kg}\cdot\text{m}^2$
Mô men quán tính quanh trục đứng $I_z$ 1.561 2.102 2.812 $\text{kg}\cdot\text{m}^2$
Độ cứng trượt bánh trước $C_{\alpha f}$ 24.972 28.072 31.132 N/rad
Độ cứng trượt bánh sau $C_{\alpha r}$ 24.438 41.265 56.836 N/rad

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

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

Quá trình mô phỏng trên 3 chế độ tải trọng (không tải 1.010 kg, nửa tải 1.450 kg và toàn tải 1.890 kg) mang lại những phát hiện định lượng quan trọng:

  • Sự dịch chuyển bất lợi của trọng tâm thân xe: Khi chất tải từ không tải lên toàn tải, chiều cao trọng tâm $h$ tăng từ 0,716 m lên 0,988 m (tăng 37,99%), đồng thời tọa độ dọc dịch chuyển về phía cầu sau thêm 0,412 m (khoảng cách $a_1$ tăng 33,33% từ 1,236 m lên 1,648 m). Điều này làm tăng mô men lật quanh trục dọc một cách đáng kể.
  • Gia tăng đột biến mô men quán tính khối lượng: Mô men quán tính lắc ngang $I_x$ tăng từ 569 lên 957 $\text{kg}\cdot\text{m}^2$ (tăng 68,19%), trong khi mô men quán tính xoay đứng $I_z$ tăng từ 1.561 lên 2.812 $\text{kg}\cdot\text{m}^2$ (tăng 80,14%), làm chậm đáng kể khả năng phản hồi và triệt tiêu dao động góc lắc của thân xe.
  • Mất cân bằng độ cứng bám ngang giữa hai cầu: Độ cứng trượt ngang bánh sau $C_{\alpha r}$ tăng vọt 132,57% (từ 24.438 N/rad lên 56.836 N/rad), cao hơn nhiều so với mức tăng 24,67% của bánh trước $C_{\alpha f}$ (từ 24.972 N/rad lên 31.132 N/rad), làm thay đổi căn bản đặc tính quay vòng từ thiếu lái sang thừa lái nguy hiểm ở dải tốc độ cao.
  • Thu hẹp biên độ vận tốc tới hạn khi chuyển làn: Trong bài thử chuyển làn tuần hoàn trên đường bằng, góc lắc cực đại $\varphi_{\max}$ ở trạng thái toàn tải cao gấp 2,14 lần so với không tải cùng tốc độ, khiến vận tốc giới hạn an toàn giảm từ mức 78 km/h xuống chỉ còn dưới 52 km/h.
       +-------------------------------------------------------------+
       |   BIẾN THIÊN CÁC CHỈ SỐ ĐỘNG LỰC HỌC THEO CHẾ ĐỘ TẢI TRỌNG   |
       +-------------------------------------------------------------+
  Chỉ số               Không tải (1.010 kg)          Toàn tải (1.890 kg)      Mức tăng
  ------------------------------------------------------------------------------------
  Chiều cao h (m)      [=== 0.716 m ===]             [===== 0.988 m =====]    +37.99%
  Mô men I_x (kg.m²)   [=== 569 kg.m² =]             [===== 957 kg.m² ===]    +68.19%
  Mô men I_z (kg.m²)   [=== 1561 kg.m² ]             [===== 2812 kg.m² ==]    +80.14%
  Độ cứng C_αr (N/rad) [=== 24438 N/rad]             [===== 56836 N/rad =]    +132.57%
  ------------------------------------------------------------------------------------
  Vận tốc tới hạn      [========== 78 km/h ==========]                        (Không tải)
  chuyển làn an toàn   [===== 52 km/h =====]                                  (Toàn tải: -33.3%)

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự suy giảm nghiêm trọng độ ổn định ngang khi đầy tải là do hiệu ứng cộng hưởng giữa chiều cao trọng tâm dâng cao thêm 0,272 m và sự gia tăng tải trọng tập trung cục bộ lên cầu sau. Khi thực hiện thao tác đánh lái đột ngột với góc $\delta = 0,2\text{ rad}$, lực ly tâm tác dụng lên khối lượng 1.890 kg sinh ra mô men lật $M_x$ vượt quá ngưỡng bù trừ của hệ thống treo nguyên bản ($k_\varphi = 149.470\text{ Nm/rad}$).

So sánh với các tiêu chuẩn an toàn chuyển động quốc tế, việc thân xe đạt góc lắc tĩnh $\varphi_{\text{ss}}$ vượt quá 6 độ ở dải vận tốc 60 km/h là dấu hiệu cảnh báo nguy cơ bánh xe phía trong vệt cua bị tách khỏi mặt đường dẫn đến lật xe hoàn toàn. Dữ liệu phản ứng thời gian được trực quan hóa thông qua các đồ thị dao động tắt dần của vận tốc góc lắc $p(t)$, gia tốc góc xoay $r(t)$ và phân bố lực ngang bánh xe trước - sau $F_{yf}, F_{yr}$, chứng minh rõ nét sự phụ thuộc phi tuyến tính của độ ổn định vào hệ số đàn hồi $k_\varphi$ và hệ số cản $c_\varphi$.

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

Dựa trên các phân tích định lượng, 4 nhóm giải pháp kỹ thuật và quản lý được đề xuất nhằm nâng cao tính năng an toàn ổn định ngang:

  • Nâng cấp độ cứng hệ thống chống nghiêng ngang: Bổ sung thanh cân bằng ngang (anti-roll bar) cho cầu sau và tăng độ cứng chống xoắn $k_\varphi$ thêm tối thiểu 18% (đạt ngưỡng trên 176.000 Nm/rad) nhằm giảm 25% góc nghiêng lắc cực đại khi vào cua. Thời gian triển khai thiết kế và thử nghiệm mẫu: 6 tháng; chủ thể thực hiện: Bộ phận R&D của nhà máy sản xuất ô tô.
  • Quy chuẩn hóa giới hạn xếp tải trọng thùng xe: Khống chế chiều cao xếp dỡ hàng hóa trên thùng xe nhằm đảm bảo tọa độ trọng tâm $h$ không vượt quá 0,85 m khi chở đủ 880 kg hàng, hạn chế triệt để tình trạng chở hàng cồng kềnh vượt chiều cao quy chuẩn. Áp dụng ngay lập tức; chủ thể thực hiện: Doanh nghiệp logistics và lái xe vận tải.
  • Bổ sung chỉ tiêu ổn định ngang vào quy trình kiểm định quốc gia: Xây dựng bộ tiêu chuẩn kiểm định động lực học ngang với ngưỡng gia tốc ngang cực đại cho phép không vượt quá $0,55\text{g}$ trong bài thử chuyển làn mô phỏng. Lộ trình hoàn thiện văn bản quy phạm: 12 đến 24 tháng; chủ thể thực hiện: Cục Đăng kiểm Việt Nam và Bộ Giao thông Vận tải.
  • Ứng dụng hệ thống kiểm soát ổn định điện tử (ESC): Nghiên cứu tích hợp các thuật toán can thiệp phanh chủ động từng bánh xe khi phát hiện góc trượt $\beta$ vượt ngưỡng 0,08 rad hoặc vận tốc góc lắc $p$ tăng đột ngột. Lộ trình nghiên cứu phát triển: 18 tháng; chủ thể thực hiện: Các đơn vị nghiên cứu cơ điện tử ô tô phối hợp cùng nhà sản xuất.

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

Nội dung và bộ dữ liệu của đề tài đem lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế và phát triển ô tô (R&D): Khai thác phương pháp thiết lập phương trình vi phân và dữ liệu độ cứng $k_\varphi$, $c_\varphi$ để tối ưu hóa kết cấu hệ thống treo, bố trí trọng tâm khung gầm cho các dòng xe tải nhẹ thương mại dưới 1.000 kg.
  • Cán bộ quản lý kỹ thuật và đăng kiểm viên: Sử dụng các ngưỡng vận tốc tới hạn và phương pháp mô phỏng số làm căn cứ khoa học để xây dựng bộ tiêu chí đánh giá độ ổn định ngang khi nghiệm thu, cải tạo phương tiện cơ giới.
  • Chủ doanh nghiệp vận tải và tài xế xe tải nhẹ: Nắm bắt các giới hạn vật lý về tải trọng (ngưỡng 880 kg) và vận tốc an toàn khi chuyển làn (dưới 52 km/h ở trạng thái toàn tải) để tổ chức xếp dỡ hàng hóa chuẩn xác, giảm thiểu tai nạn lật xe trong thực tế.
  • Giảng viên, học viên cao học và sinh viên kỹ thuật ô tô: Sử dụng luận văn như tài liệu học thuật mẫu mực về ứng dụng Matlab trong mô phỏng động lực học ô tô 4 bậc tự do và xử lý bài toán tiếp xúc lốp - mặt đường phi tuyến.

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

Vì sao tải trọng lại làm suy giảm nghiêm trọng độ ổn định ngang của xe tải nhẹ?
Khi chở đủ 880 kg hàng hóa, khối lượng toàn bộ tăng 87,13% (lên 1.890 kg) và kéo chiều cao trọng tâm dâng lên 0,988 m (tăng 37,99%). Trọng tâm nâng cao làm gia tăng cánh tay đòn của lực quán tính ly tâm khi quay vòng, tạo ra mô men lật $M_x$ lớn vượt quá khả năng hồi vị của hệ thống treo cơ sở.

Mô hình 2 bánh 4 bậc tự do có đủ tin cậy để đánh giá tính năng lật xe không?
Mô hình toán học 4 bậc tự do mô tả đồng thời chuyển động tịnh tiến, dạt ngang, lắc ngang quanh trục dọc và quay vòng quanh trục đứng. Nhờ tích hợp chính xác đặc tính biến dạng đàn hồi của lốp theo tải trọng thẳng đứng $F_z$, mô hình phản ánh sát thực tế trên 90% động thái chuyển động của xe thật.

Thông số độ cứng chống xoắn $k_\varphi = 149.470\text{ Nm/rad}$ có ý nghĩa gì đối với an toàn?
Giá trị $k_\varphi$ thể hiện khả năng chống lại góc nghiêng lắc thân xe khi chịu mô men ngoại lực. Độ cứng $k_\varphi$ càng lớn thì góc lắc $\varphi$ càng nhỏ, giúp giữ ổn định góc đặt bánh xe và diện tích vết tiếp xúc của lốp với mặt đường, từ đó hạn chế hiện tượng nhấc bánh và lật xe.

Vận tốc giới hạn an toàn của xe SMV T880 khi chuyển làn ở chế độ đầy tải là bao nhiêu?
Kết quả tính toán mô phỏng chuyển làn trên đường bằng xác định vận tốc tới hạn an toàn của xe khi chở toàn tải 1.890 kg là dưới 52 km/h. Vượt qua ngưỡng vận tốc này khi đánh lái chuyển làn khẩn cấp, góc lắc $\varphi$ và vận tốc góc $p$ sẽ rơi vào vùng mất ổn định mất kiểm soát lái.

Có thể áp dụng mô hình toán học trong luận văn cho các loại xe tải khác không?
Hoàn toàn có thể áp dụng. Mô hình được xây dựng dưới dạng giải thuật tổng quát trong môi trường Matlab R2012b. Người nghiên cứu chỉ cần thay đổi tập thông số hình học, phân bố tải trọng, độ cứng lò xo và giảm chấn của loại xe mới để thu được kết quả phân tích động lực học tương ứng.

Kết luận

  • Xây dựng thành công mô hình toán học 4 bậc tự do: Tích hợp đầy đủ các liên kết động học và động lực học giữa thân xe, hệ thống treo và lực tương tác lốp - đường.
  • Lượng hóa chính xác hệ thống thông số xe cơ sở SMV T880: Xác định chi tiết khối lượng, tọa độ trọng tâm và mô men quán tính trên cả 3 chế độ tải trọng (1.010 kg, 1.450 kg và 1.890 kg).
  • Chỉ rõ nguy cơ mất ổn định khi đầy tải: Tốc độ giới hạn an toàn khi chuyển làn bị suy giảm từ 78 km/h xuống dưới 52 km/h do chiều cao trọng tâm dâng lên mức 0,988 m.
  • Xác lập ảnh hưởng cốt lõi của hệ thống treo: Khẳng định vai trò quyết định của độ cứng chống xoắn $k_\varphi$ và hệ số cản $c_\varphi$ trong việc khống chế góc lắc ngang thân xe.
  • Lộ trình phát triển thực tế: Cần sớm tiến hành thử nghiệm kiểm chứng trên đường thực nghiệm (proving ground) và tích hợp các thuật toán điều khiển phanh chủ động ESC trong giai đoạn tiếp theo.

Các cơ quan đăng kiểm, doanh nghiệp sản xuất lắp ráp ô tô và các đơn vị vận tải cần nhanh chóng ứng dụng các kết quả nghiên cứu định lượng này vào thực tiễn thiết kế, kiểm định và vận hành nhằm đảm bảo an toàn giao thông tối đa trên mọi cung đường.