Tổng quan về luận án

Sự bùng nổ của hạ tầng giao thông đường bộ và các tuyến cao tốc huyết mạch tại Việt Nam đã thúc đẩy mạnh mẽ nhu cầu vận tải hành khách liên tỉnh bằng ô tô khách cỡ lớn, với tốc độ vận hành cho phép đạt từ 80 km/h đến 120 km/h. Theo Quyết định số 1211/QĐ-TTg của Thủ tướng Chính phủ về Quy hoạch phát triển giao thông vận tải đường bộ, mục tiêu đến năm 2025 các nhà sản xuất nội địa phải đáp ứng 92% nhu cầu xe khách trong nước với sản lượng 29.102 xe. Tuy nhiên, đặc thù kết cấu của xe khách cỡ lớn là diện tích cản gió mặt bên rất lớn ($A_{side} \gg A_{front}$), trọng tâm xe ($h_g$) đặt cao và hệ số cản hình học phức tạp. Khi vận hành ở dải tốc độ cao qua các vùng địa hình trống trải, cầu vượt biển hoặc trong điều kiện giông bão bất ngờ, tác động cộng hưởng của lực và mô men khí động học có thể làm suy giảm nghiêm trọng độ ổn định chuyển động, gây ra hiện tượng lệch làn đường, trượt ngang hoặc lật xe thảm khốc.

Research gap cốt lõi trong các công trình nghiên cứu động lực học ô tô tại Việt Nam là việc thiếu vắng một công trình tích hợp hoàn chỉnh giữa việc tính toán chính xác trường khí động 3D bằng phương pháp số với mô hình động lực học không gian phi tuyến nhiều bậc tự do. Hầu hết các nghiên cứu nội địa trước đây (như Tô Hoàng Tùng, 2015; Tạ Tuấn Hưng, 2018; Nguyễn Tiến Dũng, 2020) chỉ dừng lại ở việc khảo sát riêng lẻ hệ số cản chính diện $C_x$ nhằm tối ưu hóa tiêu hao nhiên liệu, hoặc sử dụng các hệ số khí động giả định cố định khi phân tích động lực học đoàn xe sơ mi rơ-moóc mà chưa khảo sát tương tác đa chiều của 6 thành phần lực và mô men khí động biến thiên theo góc gió xiên. Trên bình diện quốc tế, các công trình của Baker (1986), Oskar Hellsten (2020) hay Dragan Sekulic (2021) tuy đã phân tích sâu về khí động học nhưng chủ yếu tập trung vào các dòng xe buýt tiêu chuẩn châu Âu hoặc xe buýt hai tầng, chưa phản ánh chính xác đặc tính hình học và thông số phân bố khối lượng của các dòng xe khách 47 chỗ/giường nằm phổ biến tại Đông Nam Á.

Luận án "Nghiên cứu ảnh hưởng của lực khí động đến sự ổn định hướng chuyển động của ô tô khách cỡ lớn" do nghiên cứu sinh Đỗ Tiến Quyết thực hiện dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Trọng Hoan và TS. Trịnh Minh Hoàng tại Đại học Bách khoa Hà Nội (Mã số chuyên ngành: 9520116 - Kỹ thuật cơ khí động lực, 2023) đã giải quyết trọn vẹn khoảng trống học thuật này thông qua hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  • Research Question 1 (RQ1): Quy luật phân bố trường áp suất, trường vận tốc và các hệ số lực/mô men khí động (6 bậc tự do: $C_x, C_y, C_z, C_{Mx}, C_{My}, C_{Mz}$) của xe khách cỡ lớn biến thiên như thế nào theo góc gió xiên tương đương ($\beta_w \in [0^\circ, 90^\circ]$) và vận tốc gió?
  • Research Question 2 (RQ2): Mức độ ảnh hưởng định lượng của các thành phần lực ngang khí động ($F_{wy}$) và mô men xoay thân xe khí động ($M_{wz}$) đến các chỉ tiêu ổn định hướng chuyển động (độ lệch bên $y$, gia tốc bên $a_y$, vận tốc góc xoay $\dot{\psi}$, và phản lực pháp tuyến lốp $F_{zij}$) khi xe chạy ở tốc độ cao (80 - 120 km/h) là bao nhiêu?
  • Hypothesis 1 (H1): Việc áp dụng mô hình dòng rối Navier-Stokes trung bình hóa Reynolds (RANS) kết hợp biến thể Realizable $k-\varepsilon$ và xử lý lớp biên sát thành ($y^+$) cho phép dự báo chính xác hiện tượng tách dòng và hình thành xoáy áp suất âm tại vùng đuôi xe, khớp với dữ liệu thực nghiệm ống khí động học với sai số dưới 8%.
  • Hypothesis 2 (H2): Mô hình động lực học không gian 14 bậc tự do tích hợp mô hình lốp phi tuyến Pacejka sẽ chỉ ra rằng khi xe chạy với vận tốc 120 km/h chịu tác động của gió ngang cấp độ lớn, độ lệch quỹ đạo bên sẽ vượt ngưỡng an toàn làn đường tiêu chuẩn ($y > 0.5\text{ m}$) trong khoảng thời gian phản xạ ngắn ($t < 0.5\text{ s}$), đòi hỏi phải thiết lập hệ thống cảnh báo hoặc can thiệp chủ động.

Đối tượng nghiên cứu thực nghiệm và mô phỏng điển hình là dòng xe khách Thaco HB120SL-H380R-14. Khung lý thuyết của luận án được xây dựng dựa trên sự liên kết giữa Cơ học chất lưu tính toán (CFD - Computational Fluid Dynamics), Động lực học hệ nhiều vật (Multibody System Dynamics), Lý thuyết tiếp xúc lốp - mặt đường phi tuyến (Pacejka's Magic Formula) và Tiêu chuẩn an toàn phương tiện quốc tế ISO 12021 cùng Tiêu chuẩn ổn định khí động Baker.

Literature Review và Positioning

Khảo cứu lịch sử phát triển khí động học và ổn định hướng của phương tiện giao thông đường bộ cho thấy hai trường phái tiếp cận song song nhưng thường bị phân tách: trường phái khí động học thuần túy (Aerodynamics stream) và trường phái động lực học phương tiện (Vehicle dynamics stream).

                      +-------------------------------------------------------------+
                      |         KHOA HỌC ĐỘNG LỰC HỌC KHÍ ĐỘNG Ô TÔ TOÀN CẦU       |
                      +-------------------------------------------------------------+
                                                     |
                    +--------------------------------+--------------------------------+
                    |                                                                 |
                    v                                                                 v
+---------------------------------------+         +-------------------------------------------------------+
|  KHÍ ĐỘNG HỌC & TRƯỜNG DÒNG RỐI (CFD)  |         |        ĐỘNG LỰC HỌC PHƯƠNG TIỆN & ỔN ĐỊNH HƯỚNG       |
+---------------------------------------+         +-------------------------------------------------------+
| - Steady/Unsteady Flow:               |         | - Tiêu chí ổn định & giới hạn rủi ro:                 |
|   Pascal Theissen (2012),             |         |   Baker (1986), Hucho (1998),                         |
|   Makoto Tsubokura (2010),            |         |   Emmelman (1981), Carr (1982).                       |
|   Takuji Nakashima (2013).            |         | - Mô hình hóa hệ thống & mô phỏng số:                 |
| - Tương quan lực/mô men theo góc gió: |         |   Kee (2008), S. Walczak (2012),                      |
|   Youhanna E. Wiliam (2013),          |         |   Milan Batista (2014), Zhang (2015).                 |
|   Huoyue Xiang (2017),                |         | - Nghiên cứu nội địa (thiếu tương tác CFD đầy đủ):    |
|   Oskar Hellsten (2020),              |         |   Đào Đình Tại (2012), Tạ Tuấn Hưng (2018),           |
|   Dragan Sekulic (2021).              |         |   Nguyễn Anh Tuấn (2019), Nguyễn Tiến Dũng (2020).    |
+---------------------------------------+         +-------------------------------------------------------+
                    |                                                                 |
                    +--------------------------------+--------------------------------+
                                                     |
                                                     v
                      +-------------------------------------------------------------+
                      |                 POSITIONING CỦA LUẬN ÁN                     |
                      |   (Đỗ Tiến Quyết, PGS.TS Nguyễn Trọng Hoan, 2023)           |
                      |   - Tích hợp 2 chiều: CFD 3D (RANS Realizable k-epsilon)   |
                      |     và Mô hình động lực học không gian 14 DOFs.             |
                      |   - Đối tượng thực tiễn: Khung vỏ xe Thaco HB120SL.         |
                      |   - Kiểm chứng thực nghiệm ống khí động cân 6 thành phần.   |
                      +-------------------------------------------------------------+

Trong dòng nghiên cứu khí động học, cuộc tranh luận lớn xoay quanh việc lựa chọn mô hình gió ngang ổn định (steady), bán ổn định (quasi-steady) hay bất ổn định (unsteady). Pascal Theissen (2012) khi nghiên cứu tác động của gió giật đã chỉ ra rằng các mô hình quasi-steady có xu hướng đánh giá thấp giá trị tức thời của mô men xoay thân xe ($M_{wz}$) so với mô hình gió giật thực tế, dù lực ngang trung bình có thể xấp xỉ nhau. Makoto Tsubokura (2010) đã sử dụng siêu máy tính SR11000 (Hitachi) tại Đại học Tokyo để mô phỏng dòng rối LES cho xe tải với hàm gió biến thiên hình sin $v = A \sin^2(x/\lambda - ft)$, chứng minh cấu trúc xoáy không gian tạo ra các xung lực ngang có tần số cộng hưởng với hệ thống treo. Ngược lại, Takuji Nakashima (2013) và Youhanna E. Wiliam (2013) lập luận rằng trong bài toán đánh giá đáp ứng ổn định cơ bản, việc quy đổi tác động của vận tốc xe $V_{xe}$ và gió ngang $V_{giongang}$ thành một véctơ vận tốc gió xiên tương đương $V_{gioxien}$ với góc nghiêng $\beta_w = \arctan(V_{giongang}/V_{xe})$ là giải pháp tối ưu về mặt chi phí tính toán nhưng vẫn bảo toàn đầy đủ các đặc trưng áp suất bề mặt tĩnh.

So sánh đối sánh quốc tế cho thấy tính dị biệt sâu sắc về hình thái học phương tiện:

  1. Nghiên cứu của Oskar Hellsten (2020) trên dòng xe khách Volvo Bus tiêu chuẩn châu Âu ghi nhận hệ số cản $C_x$ tăng khi góc nghiêng $\beta_w$ từ $0^\circ$ đến $20^\circ$, sau đó giảm mạnh và chuyển sang giá trị âm khi $\beta_w > 60^\circ$ do áp suất vùng đuôi xe cao hơn áp suất vùng đầu xe bị tách dòng.
  2. Công trình kinh điển của Baker (1986) trên dòng xe buýt hai tầng Leyland Atlantean cho thấy $C_x$ đạt cực đại tại $\beta_w \approx 30^\circ$ và lực ngang $C_y$ đạt giá trị lớn nhất tại $\beta_w = 90^\circ$ do thân xe có dạng khối hộp chữ nhật phẳng.
  3. Nghiên cứu của Dragan Sekulic (2021) trên xe khách liên tỉnh một tầng lại chứng minh $C_y$ tăng tuyến tính từ $0^\circ$ đến $45^\circ$ và đạt bão hòa ở dải góc lớn hơn.

Luận án của Đỗ Tiến Quyết đã định vị chính xác điểm giao thoa này: thiết lập mô hình khí động học 3D cho biên dạng xe khách Thaco đặc thù tại Việt Nam (có kính chắn gió nghiêng, góc bo tròn cạnh trên và bề mặt phẳng kéo dài), từ đó cung cấp bộ hệ số khí động 6 thành phần chính xác làm đầu vào cho bài toán động lực học phi tuyến 14 bậc tự do.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng khung lý thuyết khí động học ứng dụng và động lực học phương tiện thông qua các đóng góp khoa học cụ thể:

  • Phát triển mô hình tương tác khí động - động lực học đa chiều: Thay vì đơn giản hóa lực khí động thành một lực tập trung đặt tại trọng tâm xe, nghiên cứu đã tổng quát hóa hệ 6 thành phần lực và mô men khí động ($F_{wx}, F_{wy}, F_{wz}, M_{wx}, M_{wy}, M_{wz}$) tác động liên tục vào khối lượng được treo, biểu diễn dưới dạng hàm phi tuyến của góc gió xiên tương đương $\beta_w$ và bình phương vận tốc dòng khí tổng hợp $U_\infty^2$.
  • Mở rộng lý thuyết ổn định chuyển động của Baker và Hucho: Nghiên cứu đã xác lập mối quan hệ hàm số giữa vận tốc xe chạy, cấp độ gió ngang và các biến trạng thái mất ổn định (vận tốc góc xoay thân xe quanh trục đứng $\dot{\psi}$, góc lắc ngang của thân xe $\varphi$, góc lắc dọc $\theta$ và gia tốc bên $a_y$).
  • Lý giải cơ chế chuyển dịch tâm áp suất khí động học: Kết quả mô phỏng RANS đã làm sáng tỏ hiện tượng tâm áp suất khí động di chuyển dọc theo chiều dài cơ sở ($l$) khi góc gió xiên thay đổi, giải thích bản chất vật lý của sự xuất hiện mô men xoay thân xe làm mất cân bằng hướng lái ngay cả khi góc đánh lái $\delta_{ij} = 0$.
+----------------------------------------------------------------------------------------------------+
|                                    KHUNG PHÂN TÍCH TỔNG QUÁT                                       |
+----------------------------------------------------------------------------------------------------+
                                                  |
           +--------------------------------------+--------------------------------------+
           |                                                                             |
           v                                                                             v
+------------------------------------+                         +-------------------------------------+
|      PHÂN HỆ KHÍ ĐỘNG HỌC (CFD)     |                         |     PHÂN HỆ ĐỘNG LỰC HỌC PHƯƠNG TIỆN|
+------------------------------------+                         +-------------------------------------+
| - Phương trình RANS:               |                         | - Khối lượng được treo (Sprung Mass)|
|   $\nabla \cdot \vec{u} = 0$       |                         |   6 DOFs: $x, y, z, \varphi, \theta,|
|   Phương trình Navier-Stokes       |                         |   \psi$                             |
|   trung bình hóa.                  |                         | - Khối lượng không được treo (Axles)|
| - Mô hình Realizable $k-\varepsilon$:      |                         |   4 DOFs: Chuyển vị đứng và góc lắc |
|   $k$ (Động năng dòng rối),        |                         |   ngang cầu trước/sau.              |
|   $\varepsilon$ (Hệ số tán xạ).    |                         | - Bánh xe: 4 DOFs chuyển động quay. |
| - Ten-sơ ứng suất Reynolds:        |                         | - Mô hình lốp phi tuyến Pacejka:    |
|   $\tau_{ij}^t = \mu_t \partial_j u_i$                         |   $F_{yij} = f(\alpha_{ij}, F_{zij})|
+------------------------------------+                         +-------------------------------------+
           |                                                                             |
           +--------------------------------------+--------------------------------------+
                                                  |
                                                  v
+----------------------------------------------------------------------------------------------------+
|             BỘ THÔNG SỐ TRUYỀN DẪN VÀ CHỈ BÁO ĐÁNH GIÁ ỔN ĐỊNH HƯỚNG                               |
|   Lực & Mô men khí động: $F_{wx}, F_{wy}, F_{wz}, M_{wx}, M_{wy}, M_{wz}$                           |
|   Chỉ tiêu đánh giá (ISO 12021 & Baker): $y(t), a_y(t), \dot{\psi}(t), F_{zij}(t) \to 0$           |
+----------------------------------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp chặt chẽ 3 nền tảng lý thuyết cốt lõi:

  1. Lý thuyết thủy khí động học chất lưu không nén được (Incompressible Fluid Dynamics): Khai thác hệ phương trình liên tục và phương trình động lượng Navier-Stokes trong điều kiện số Mach $M < 0.3$: $$\frac{\partial u_i}{\partial x_i} = 0$$ $$\frac{\partial u_i}{\partial t} + u_j \frac{\partial u_i}{\partial x_j} = -\frac{1}{\rho}\frac{\partial p}{\partial x_i} + \frac{\partial}{\partial x_j} \left[ \nu \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right) - \overline{u_i' u_j'} \right]$$
  2. Lý thuyết động lực học hệ nhiều vật (Multibody Vehicle Dynamics): Thiết lập hệ phương trình vi phân chuyển động phi tuyến theo nguyên lý D'Alembert và định lý biến thiên động lượng/mô men động lượng cho hệ gồm 1 khối lượng được treo (thân xe) liên kết đàn hồi với 2 khối lượng không được treo (cầu trước và sau).
  3. Lý thuyết tiếp xúc ma sát và biến dạng lốp phi tuyến Pacejka: Mô tả tương tác lực ngang bánh xe $F_{yij}$ phụ thuộc vào góc trượt lệch bên $\alpha_{ij}$, tải trọng động thẳng đứng $F_{zij}$ và hệ số bám $\mu$: $$F_{yij} = D \sin \left( C \arctan \left( B \alpha_{ij} - E (B \alpha_{ij} - \arctan(B \alpha_{ij})) \right) \right)$$

Điều kiện biên xác định: Mặt đường hoàn toàn phẳng, không xét biến dạng gồ ghề; lốp xe hoạt động trong giới hạn bám không bị khóa cứng bánh; không xét phản xạ bù lái chủ động từ người điều khiển để đánh giá độ nhạy thuần túy của phương tiện.

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

Thiết kế nghiên cứu

Luận án tuân thủ triệt để quan điểm bản thể luận thực chứng (Positivism) và nhận thức luận thực nghiệm suy diễn (Deductive Empirical Epistemology). Thiết kế nghiên cứu là sự kết hợp đồng bộ (Mixed-Method Quantitative Approach) giữa mô phỏng số học chuyên sâu độ phân giải cao (High-Fidelity Numerical Simulation) và kiểm chứng thực nghiệm vật lý trong phòng thí nghiệm đạt chuẩn (Wind Tunnel Experimental Verification).

Cấu trúc mô hình động lực học được mô hình hóa với 14 bậc tự do độc lập (14 DOFs):

  • Khối lượng được treo (Thân xe - Sprung mass $M$): 6 bậc tự do gồm 3 chuyển vị tịnh tiến (dọc $x$, ngang $y$, đứng $z$) và 3 chuyển động quay quanh hệ trục tọa độ trọng tâm (lắc ngang $\varphi$, lắc dọc $\theta$, xoay thân xe quanh trục đứng $\psi$).
  • Khối lượng không được treo cầu trước ($m_{A1}$) và sau ($m_{A2}$): 4 bậc tự do gồm chuyển vị thẳng đứng ($z_{A1}, z_{A2}$) và góc lắc ngang của từng cầu quanh trục dọc ($A_1, A_2$).
  • Hệ thống 4 cụm bánh xe: 4 bậc tự do biểu diễn chuyển động quay tự do quanh trục ngang của các bánh xe ($J_{Ay11}, J_{Ay12}, J_{Ay21}, J_{Ay22}$).

Đối tượng mô phỏng và mẫu thí nghiệm tương ứng là dòng xe khách Thaco HB120SL-H380R-14 với các thông số hình học và cơ học chính xác:

  • Chiều dài toàn bộ ($L$): $12.000\text{ mm}$; Chiều rộng toàn bộ ($B$): $2.500\text{ mm}$; Chiều cao toàn bộ ($H$): $3.500\text{ mm}$.
  • Chiều dài cơ sở ($l$): $6.000\text{ mm}$; Chiều cao trọng tâm toàn tải ($h_g$): $1.450\text{ mm}$.
  • Vết bánh xe trước ($2b_1$) và sau ($2b_2$): $2.050\text{ mm}$ và $1.860\text{ mm}$.
  • Khối lượng toàn bộ đầy tải ($M$): $16.000\text{ kg}$ (phân bố cầu trước $M_{c1} = 6.000\text{ kg}$, cầu sau $M_{c2} = 10.000\text{ kg}$).
                                  +-------------------------------------------------------------+
                                  |            QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU RIGOROUS        |
                                  +-------------------------------------------------------------+
                                                                 |
                                +--------------------------------+--------------------------------+
                                |                                                                 |
                                v                                                                 v
+---------------------------------------------------------------+ +---------------------------------------------------------------+
|         MÔ PHỎNG SỐ CFD (ANSYS FLUENT)                        | |        THỰC NGHIỆM ỐNG KHÍ ĐỘNG (WIND TUNNEL)                 |
+---------------------------------------------------------------+ +---------------------------------------------------------------+
| 1. Xây dựng mô hình 3D xe Thaco HB120SL.                      | | 1. Chế tạo mô hình thu nhỏ tỷ lệ hình học chính xác.          |
| 2. Thiết lập miền tính toán khí động (Computational Domain).  | | 2. Lắp đặt trên cân khí động 6 thành phần lực/mô men.         |
| 3. Chia lưới Tetrahedron/Prism, kiểm soát $y^+$ bề mặt.       | | 3. Thu thập dữ liệu qua card cảm biến & phần mềm DasyLab.     |
| 4. Thiết lập mô hình dòng rối Realizable $k-\varepsilon$.             | | 4. Quét dải vận tốc dòng khí và góc xoay mô hình $\beta_w$.   |
| 5. Giải hệ RANS xác định $C_x, C_y, C_z, C_{Mx}, C_{My}, C_{Mz}$.| | 5. Trích xuất lực cản $F_{xkd}$ và lực ngang $F_{ykd}$.       |
+---------------------------------------------------------------+ +---------------------------------------------------------------+
                                |                                                                 |
                                +--------------------------------+--------------------------------+
                                                                 |
                                                                 v
                                +---------------------------------------------------------------+
                                | KIỂM ĐỊNH TƯƠNG QUAN & HIỆU CHUẨN SAI SỐ (< 8%)               |
                                +---------------------------------------------------------------+
                                                                 |
                                                                 v
                                +---------------------------------------------------------------+
                                | MÔ PHỎNG ĐỘNG LỰC HỌC 14 BẬC TỰ DO (ODE NUMERICAL SOLVER)     |
                                | - Nhập ngoại lực khí động biến thiên 6 thành phần.            |
                                | - Tích hợp mô hình lốp phi tuyến Pacejka Magic Formula.       |
                                | - Đánh giá đáp ứng ổn định: $y(t), a_y(t), \dot{\psi}(t), F_z$|
                                +---------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình giải quyết bài toán khí động số trên phần mềm Ansys Fluent được chuẩn hóa nghiêm ngặt:

  • Xây dựng miền tính toán (Domain): Miền không gian bao quanh xe được mở rộng với chiều dài $8L$, chiều rộng $5B$ và chiều cao $4H$ nhằm loại trừ hoàn toàn ảnh hưởng phản xạ sóng biên (Blocking ratio < 3%).
  • Tối ưu hóa chất lượng lưới (Mesh Generation): Sử dụng lưới phi cấu trúc phức hợp (Unstructured Tetrahedral Mesh) kết hợp các lớp lăng trụ biên (Prism inflation layers) sát bề mặt thân vỏ xe. Số lượng phần tử lưới đạt xấp xỉ $3.2 \times 10^6$ nodes. Chỉ số chất lượng lưới trực giao (Orthogonal Quality) đạt trung bình $> 0.85$ và hệ số méo phần tử (Skewness) $< 0.28$. Giá trị thông số lớp biên không thứ nguyên $y^+$ trên toàn bộ diện tích vỏ xe được duy trì trong dải tối ưu ($30 \le y^+ \le 120$) phù hợp với hàm thế vách tiêu chuẩn (Standard Wall Functions).
  • Thiết lập điều kiện biên và thuật giải: Biên vào (Inlet) là Velocity Inlet ($U_\infty = 20 - 40\text{ m/s}$); Biên ra (Outlet) là Pressure Outlet ($p = 0\text{ Pa gauge}$); Mặt đường và bề mặt vỏ xe là No-slip Wall; Thuật toán liên kết áp suất - vận tốc SIMPLEC với lược đồ sai phân bậc hai (Second Order Upwind) cho các phương trình bảo toàn động lượng, động năng rối $k$ và hệ số tán xạ $\varepsilon$.

Data và phân tích

Trong nghiên cứu thực nghiệm, hệ thống thí nghiệm ống khí động học tĩnh hở (Open-jet / Closed-test-section Wind Tunnel) được vận hành với động cơ quạt gió công suất lớn, vận tốc dòng khí hiệu chuẩn từ $10\text{ m/s}$ đến $35\text{ m/s}$. Mô hình thu nhỏ của xe Thaco HB120SL được gá đặt trên cụm cân khí động học tự động 6 thành phần lực/mô men. Tín hiệu điện áp analog từ các cảm biến biến dạng (Strain-gauge load cells) được lọc nhiễu, khuếch đại và số hóa qua bộ chuyển đổi A/D chuyên dụng, sau đó ghi nhận và xử lý thời gian thực bằng phần mềm DasyLab.

Các chỉ số thống kê độ tin cậy của mô hình đo đạc:

  • Hệ số tương quan xác định giữa dữ liệu thực nghiệm và mô phỏng CFD đạt $R^2 = 0.968$ đối với lực cản $F_{xkd}$ và $R^2 = 0.982$ đối với lực ngang $F_{ykd}$.
  • Mức độ sai lệch tương đối giữa mô phỏng và thực nghiệm qua các góc nghiêng $0^\circ, 30^\circ, 60^\circ$ dao động trong khoảng từ $4.2%$ đến $7.8%$, khẳng định tính hội tụ và tính giá trị nội tại (Internal Validity) của mô hình CFD.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Quy luật phi tuyến đảo chiều của hệ số cản $C_x$ và sự bùng nổ của hệ số lực ngang $C_y$: Khi góc gió xiên tương đương $\beta_w$ tăng từ $0^\circ$ đến $30^\circ$, hệ số cản khí động $C_x$ tăng nhẹ từ $0.58$ lên cực đại $0.72$. Tuy nhiên, khi $\beta_w$ tiếp tục tăng từ $30^\circ$ lên $60^\circ$ và $90^\circ$, $C_x$ sụt giảm đột ngột và chuyển sang giá trị âm (tại $\beta_w = 60^\circ$, $C_x \approx -0.15$). Hiện tượng này được giải thích do sự phân bố lại trường áp suất: vùng xoáy áp suất âm bao trùm mặt khuất gió kéo thân xe về phía sau bị lấn át bởi áp suất dương rất lớn đẩy vào sườn xe phía đón gió. Đồng thời, hệ số lực ngang $C_y$ tăng vọt từ $0$ (tại $\beta_w = 0^\circ$) lên cực đại $1.85$ (tại $\beta_w = 60^\circ$), lớn gấp hơn 3 lần hệ số cản chính diện ban đầu.

  2. Sự hình thành mô men xoay thân xe $M_{wz}$ và mô men lật $M_{wx}$ quy mô lớn: Tại vận tốc gió xiên $V_{gioxien} = 30\text{ m/s}$ và góc $\beta_w = 30^\circ$, lực ngang khí động đạt giá trị $F_{wy} = 14.850\text{ N}$ và mô men xoay thân xe đạt $M_{wz} = 18.200\text{ Nm}$. Giá trị mô men này có xu hướng đẩy đầu xe xoay ngược chiều gió trong khi đuôi xe bị trôi, gây mất cân bằng mô men quanh trục thẳng đứng $Oz$.

  3. Mất ổn định quỹ đạo bên ở dải tốc độ cao (80 km/h và 120 km/h): Kết quả giải hệ phương trình vi phân động lực học 14 bậc tự do cho thấy: Khi xe chạy với vận tốc $80\text{ km/h}$ gặp luồng gió ngang $20\text{ m/s}$ ($\approx 72\text{ km/h}$, tương đương gió bão cấp 8-9), sau thời gian $t = 0.8\text{ s}$, độ lệch bên của thân xe đạt $y = 0.62\text{ m}$ và gia tốc bên đạt $a_y = 1.85\text{ m/s}^2$. Khi vận tốc xe tăng lên $120\text{ km/h}$, sau $0.5\text{ s}$ đi vào vùng gió ngang, độ lệch quỹ đạo bên vọt lên $y = 1.15\text{ m}$, vượt xa bề rộng an toàn cho phép của tim làn đường ($y_{allow} \approx 0.5\text{ m}$), trong khi vận tốc góc xoay thân xe đạt $\dot{\psi} = 0.185\text{ rad/s}$.

  4. Hiện tượng suy giảm tải trọng thẳng đứng gây nguy cơ mất bám: Dưới tác động đồng thời của lực nâng khí động $F_{wz}$ ($C_z \approx 0.35$) và mô men lắc ngang $M_{wx}$, tải trọng thẳng đứng tại các bánh xe phía đón gió bên trái giảm tới $42%$ ở vận tốc $80\text{ km/h}$ và giảm tới $68%$ ở vận tốc $120\text{ km/h}$. Mặc dù chưa chạm ngưỡng lật bánh hoàn toàn ($F_{zij} \to 0$), sự suy giảm tải trọng thẳng đứng này làm suy giảm nghiêm trọng khả năng truyền lực bên của lốp theo mô hình Pacejka, khiến xe rơi vào trạng thái trượt ngang không kiểm soát.

+----------------------------------------------------------------------------------------------------+
|                BẢNG TỔNG HỢP SO SÁNH CÁC HỆ SỐ KHÍ ĐỘNG VÀ ĐÁP ỨNG ỔN ĐỊNH HƯỚNG                   |
+----------------------+--------------------+--------------------+-----------------------------------+
| Thông số khảo sát    | Góc $\beta_w = 0^\circ$| Góc $\beta_w = 30^\circ$| Góc $\beta_w = 60^\circ$          |
+----------------------+--------------------+--------------------+-----------------------------------+
| Hệ số cản $C_x$      | $0.58$ (Dương lớn) | $0.72$ (Cực đại)   | $-0.15$ (Đảo chiều âm)            |
| Hệ số lực ngang $C_y$| $0.00$             | $1.24$             | $1.85$ (Cực đại)                  |
| Hệ số mô men $C_{Mz}$| $0.00$             | $0.38$             | $0.52$                            |
| Độ lệch bên ($120$km/h)| $0.00\text{ m}$  | $1.15\text{ m}$ (sau $0.5$s) | $1.82\text{ m}$ (Nguy cơ lật/trượt cực cao) |
| Tải trọng lốp đón gió| $100\%$ định mức   | Giảm $48\%$ tải    | Giảm $68\%$ tải (Trượt phẳng)     |
+----------------------+--------------------+--------------------+-----------------------------------+

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ cơ sở dữ liệu số thực nghiệm chuẩn hóa về 6 hệ số khí động học cho xe khách cỡ lớn dạng thân hộp chữ nhật đặc thù, đóng vai trò tài liệu đối sánh cho các nghiên cứu cơ học động lực học ô tô tại các nước đang phát triển.
  • Về mặt kỹ thuật công nghiệp: Các nhà sản xuất lắp ráp ô tô (như THACO, SAMCO, TRACOMECO) có thể sử dụng kết quả này để tối ưu hóa hình dáng khí động học thân vỏ: bo tròn các cột biên chữ A, hạ thấp trọng tâm cụm mui xe, thiết kế các cánh hướng dòng khuếch tán không khí dưới gầm (Rear Diffuser) nhằm giảm hệ số lực nâng $C_z$ và lực ngang $C_y$.
  • Về mặt an toàn giao thông và chính sách: Cung cấp bằng chứng khoa học thực chứng để Cục Đường bộ và các cơ quan quản lý vận tải xây dựng quy chuẩn vận hành an toàn: Thiết lập biển báo cảnh báo gió ngang và hạ giới hạn tốc độ tối đa từ $120\text{ km/h}$ xuống $70 - 80\text{ km/h}$ đối với ô tô khách khi lưu thông qua cầu vượt biển hoặc đèo dốc trong điều kiện thời tiết gió giật cấp 7 trở lên.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Giả thiết trường gió ổn định (Steady Wind Assumption): Luận án giả định véctơ gió ngang có hướng và độ lớn không đổi theo thời gian tại mỗi bước giải. Trong tự nhiên, gió giật (Turbulent Wind Gusts) có phổ tần số biến thiên ngẫu nhiên có thể kích thích các mode dao động riêng của hệ thống treo.
  2. Loại trừ phản xạ điều khiển của người lái (Closed-loop Driver Absence): Mô hình động lực học được khảo sát ở trạng thái vòng hở (Open-loop system), góc đánh lái cố định $\delta = 0$. Trong thực tế, phản xạ bù lái của tài xế khi nhận biết xe bị trôi làn có thể gây ra hiện tượng đánh lái quá mức (Oversteering), dẫn đến mất ổn định thứ cấp.
  3. Mặt đường phẳng lý tưởng: Chưa tích hợp mấp mô mặt đường ngẫu nhiên theo tiêu chuẩn ISO 8608, vốn có thể làm thay đổi động học góc đặt bánh xe và phản lực lốp liên tục.

Các định hướng nghiên cứu tiếp nối trong giai đoạn 5 - 10 năm tới bao gồm:

  • Mở rộng mô hình CFD sang phương pháp mô phỏng xoáy lớn LES (Large Eddy Simulation) tích hợp dòng gió giật phi dừng đa chiều.
  • Phát triển mô hình điều khiển vòng kín Driver-in-the-loop tích hợp hệ thống cân bằng điện tử chủ động ESP/ESC và hệ thống lái chủ động AFS can thiệp phân bổ mô men phanh độc lập tại các bánh xe để triệt tiêu mô men xoay $M_{wz}$.
  • Nghiên cứu khí động học cho đoàn xe khách chạy bám đuôi tự hành (Platooning) trên cao tốc thông minh.

Tác động và ảnh hưởng

  • Tác động học thuật: Luận án tạo ra một bước chuyển dịch phương pháp luận từ mô phỏng khí động tĩnh đơn thuần sang phân tích tích hợp động lực học nhiều vật phi tuyến. Dự kiến công trình sẽ là nguồn tài liệu trích dẫn quan trọng cho các nghiên cứu sinh, học viên cao học chuyên ngành Kỹ thuật Cơ khí Động lực, Kỹ thuật Giao thông và Cơ học Ứng dụng.
  • Chuyển giao công nghệ công nghiệp: Cung cấp giải pháp phần mềm và quy trình tính toán cho Trung tâm R&D của các tập đoàn sản xuất ô tô trong nước, giúp rút ngắn thời gian thử nghiệm thực tế và giảm chi phí thiết kế khung vỏ xe khách thế hệ mới.
  • Lợi ích xã hội và an toàn cộng đồng: Việc ứng dụng các ngưỡng cảnh báo tốc độ trong điều kiện gió bão sẽ góp phần giảm thiểu trực tiếp các vụ tai nạn giao thông thảm khốc liên quan đến lật xe và trượt lệch làn đường của phương tiện vận tải hành khách khối lượng lớn.

Đối tượng hưởng lợi

+----------------------------------------------------------------------------------------------------+
|                                    MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI                                     |
+-------------------------+--------------------------------------------------------------------------+
| Nhóm đối tượng          | Giá trị thụ hưởng và Ứng dụng thực tiễn                                  |
+-------------------------+--------------------------------------------------------------------------+
| Nghiên cứu sinh &       | - Khung phương pháp luận kết hợp CFD RANS và Động lực học 14 DOFs.       |
| Giảng viên Học thuật    | - Cơ sở dữ liệu 6 thành phần lực/mô men khí động xe khách dạng hộp.      |
+-------------------------+--------------------------------------------------------------------------+
| Kỹ sư R&D Ô tô          | - Bộ thông số khí động học định lượng phục vụ thiết kế kiểu dáng khí động|
| (THACO, SAMCO, VinFast) | - Quy trình mô phỏng kiểm nghiệm ảo trước khi chế tạo dưỡng mẫu thực.    |
+-------------------------+--------------------------------------------------------------------------+
| Cơ quan Quản lý Đường bộ| - Căn cứ khoa học thiết lập ngưỡng tốc độ an toàn theo cấp gió bão.      |
| & CSGT                  | - Quy hoạch lắp đặt hệ thống cảnh báo gió ngang trên cao tốc và cầu lớn. |
+-------------------------+--------------------------------------------------------------------------+
| Doanh nghiệp Vận tải &  | - Cẩm nang kỹ năng điều khiển xe khách cỡ lớn an toàn trong bão gió.     |
| Lái xe Đường dài        | - Nâng cao ý thức kiểm soát tốc độ khi chạy qua vùng địa hình trống trải.|
+-------------------------+--------------------------------------------------------------------------+

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập mô hình toán học giải tích hóa tương tác 6 thành phần lực và mô men khí động phi tuyến biến thiên theo góc gió xiên tương đương $\beta_w$, tích hợp trực tiếp vào hệ phương trình vi phân chuyển động không gian 14 bậc tự do của xe khách cỡ lớn. Công trình đã mở rộng lý thuyết động lực học phương tiện của Wong (2001) và lý thuyết ổn định khí động học của Baker (1986) thông qua việc chứng minh rằng sự kết hợp giữa mô men xoay thân xe $M_{wz}$ và sự suy giảm tải trọng lốp pháp tuyến đón gió (lên tới $68%$) là nguyên nhân trực tiếp kích hoạt hiện tượng lệch làn đường trước khi xảy ra hiện tượng lật xe hình học.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Takuji Nakashima (2013) (chỉ mô phỏng CFD tĩnh 2 nguồn gió tốn kém tài nguyên siêu máy tính) và Zhang (2015) (dùng phần mềm TruckSim với hệ số khí động học mặc định), luận án đã tối ưu hóa mô hình 1 nguồn gió xiên tương đương trong Ansys Fluent kết hợp thuật giải Realizable $k-\varepsilon$, đồng thời tự xây dựng hệ phương trình vi phân động lực học 14 bậc tự do tường minh bản chất vật lý, được hiệu chuẩn chính xác bằng thực nghiệm ống khí động học với cân 6 thành phần chuyên dụng.

3. Phát hiện nào bất ngờ nhất về mặt số liệu thực nghiệm?

Phát hiện bất ngờ nhất là hiện tượng đảo chiều của hệ số cản chính diện $C_x$: khi góc gió xiên vượt qua ngưỡng $57^\circ - 60^\circ$, hệ số cản mang giá trị âm ($C_x \approx -0.15$), nghĩa là lực khí động theo trục dọc không còn cản trở chuyển động tiến của xe mà trở thành lực đẩy tiếp năng lượng dọc, trong khi lực ngang $C_y$ lại tăng vọt lên gấp hơn 3 lần hệ số cản ban đầu, đẩy toàn bộ rủi ro mất ổn định sang phương ngang và mô men quay thân xe quanh trục thẳng đứng.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Toàn bộ thông số kích thước hình học xe Thaco HB120SL, thông số cơ học của hệ thống treo khí nén, thông số độ cứng lốp Pacejka, cấu hình lưới không gian CFD, giá trị điều kiện biên $y^+$, cùng quy chuẩn đo đạc ống khí động theo tiêu chuẩn ISO 12021 được trình bày chi tiết trong các bảng thông số kỹ thuật của luận án, cho phép tái lập hoàn toàn quy trình tính toán trên bất kỳ nền tảng số nào.

5. Khung chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình nghiên cứu phát triển bao gồm: (1) 2024-2026: Tích hợp thuật toán điều khiển mờ thích nghi (Adaptive Fuzzy / MPC) cho hệ thống phanh vi sai chủ động chống lệch làn do gió ngang; (2) 2027-2029: Thử nghiệm trên mô hình xe khách tỷ lệ thật tích hợp cảm biến đo áp suất gió thời gian thực trên đường cao tốc ven biển; (3) 2030-2034: Phát triển chuẩn giao tiếp V2X cảnh báo sớm rủi ro khí động học cho mạng lưới phương tiện tự hành cấp độ 4 (SAE Level 4 Autonomous Buses).

Kết luận

Công trình nghiên cứu của Tiến sĩ Đỗ Tiến Quyết đã khẳng định những giá trị khoa học và thực tiễn vững chắc thông qua 5 đóng góp cốt lõi:

  1. Xây dựng hoàn chỉnh mô hình số học CFD 3D mô phỏng chính xác trường dòng chảy rối bao quanh xe khách Thaco HB120SL-H380R-14, xác định đầy đủ quy luật biến thiên của 6 hệ số lực và mô men khí động học ($C_x, C_y, C_z, C_{Mx}, C_{My}, C_{Mz}$) trong toàn dải góc gió xiên $0^\circ \le \beta_w \le 90^\circ$.
  2. Thiết lập thành công hệ phương trình vi phân động lực học không gian 14 bậc tự do phản ánh trung thực tương tác phi tuyến giữa khối lượng được treo, các cầu xe, hệ thống treo khí nén và mô hình lốp biến dạng Pacejka dưới tác động của ngoại lực khí động học.
  3. Định lượng hóa chính xác mức độ suy giảm ổn định hướng chuyển động: chứng minh rằng ở dải tốc độ cao 120 km/h dưới tác động của gió ngang mạnh ($20\text{ m/s}$), xe sẽ bị lệch làn tới $1.15\text{ m}$ chỉ sau $0.5\text{ s}$, tiềm ẩn nguy cơ va chạm và trượt ngang cực kỳ nguy hiểm nếu không có cảnh báo hoặc can thiệp điều khiển.
  4. Thực hiện kiểm chứng thực nghiệm nghiêm ngặt trong phòng thí nghiệm ống khí động học hiện đại với cân khí động 6 thành phần và phần mềm DasyLab, xác nhận độ tin cậy của mô hình mô phỏng số với sai số nhỏ dưới $8%$.
  5. Mở ra hướng nghiên cứu ứng dụng đột phá trong thiết kế công nghiệp ô tô nội địa và cung cấp luận cứ khoa học thực tiễn giúp các cơ quan quản lý nhà nước hoàn thiện quy chuẩn an toàn giao thông đường bộ trên các tuyến cao tốc hiện đại tại Việt Nam.