Tổng quan về luận án

Nghiên cứu động lực học và dao động của phương tiện giao thông đường bộ giữ vai trò then chốt trong kỹ thuật cơ khí động lực hiện đại. Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã số: 9 52 01 01) của tác giả Phùng Mạnh Cường, thực hiện tại Học viện Kỹ thuật Quân sự dưới sự hướng dẫn khoa học của PGS. Vũ Công Hàm và PGS. Trần Quang Dũng, mang tên: "Nghiên cứu dao động thẳng đứng của ô tô theo các mô hình khác nhau có tính đến hiện tượng mất liên kết giữa bánh xe và mặt đường". Công trình tạo lập một bước chuyển biến học thuật quan trọng trong lý thuyết dao động ô tô thông qua việc tích hợp đồng thời ba trường nhân tố phi tuyến: sự biến dạng động của kết cấu đường, hiện tượng mất liên kết gián đoạn (tách bánh xe) và sự biến thiên kích thước hình học kèm quy luật phân bố áp suất tại diện tích vết tiếp xúc.

       +-------------------------------------------------------+
       |   KHUNG ĐỘNG LỰC HỌC TỔNG THỂ HỆ XE - ĐƯỜNG PHI TUYẾN  |
       +-------------------------------------------------------+
                                  |
         +------------------------+------------------------+
         |                                                 |
         v                                                 v
+------------------+                             +--------------------+
|  MÔ HÌNH CƠ HỆ   |                             |   MÔ HÌNH ĐƯỜNG    |
| - 1/4 Xe (2 DOF) |                             | - Nền đàn nhớt     |
| - 1/2 Dọc        |                             |   Kelvin (kS, cS)  |
| - 1/2 Ngang      |                             | - Dầm Euler-       |
| - 3D Không gian  |                             |   Bernoulli (1D)   |
|   (7 DOF)        |                             | - Tấm mỏng         |
+------------------+                             |   Kirchhoff (2D)   |
         |                                       +--------------------+
         |                 TƯƠNG TÁC TIẾP XÚC              |
         +------------> (Unilateral Contact) <-------------+
                                  |
                   +--------------+--------------+
                   |                             |
                   v                             v
       [Trạng thái tiếp xúc s=1]     [Trạng thái tách bánh s=0]
       - Lực pháp tuyến R_j > 0      - Lực liên kết R_j = 0
       - Biến dạng vết lốp z_Lj > 0  - Khe hở tách bánh Z_Lj > 0
       - Áp suất p(x,y,t) liên tục   - Mất tính dẫn hướng/phanh

Khoảng trống học thuật (research gap) cốt lõi mà luận án định vị bắt nguồn từ sự thiếu hoàn thiện của các mô hình truyền thống. Đa số các công trình trước đây hoặc giả định liên kết song phương vĩnh cửu (bánh xe luôn bám đường) khi khảo sát mô hình không gian nhiều bậc tự do, hoặc chỉ xét hiện tượng tách bánh trên các mô hình phẳng đơn giản hóa gắn với mặt đường cứng tuyệt đối. Luận án giải quyết triệt để hạn chế này bằng cách thiết lập 4 câu hỏi nghiên cứu và giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế toán học nào cho phép hợp nhất việc mô tả động lực học của bánh xe ở cả hai trạng thái tiếp xúc và tách bánh mà không làm gián đoạn hệ phương trình trạng thái? Giả thuyết 1 (H1): Tham số trạng thái tiếp xúc logic $s(t) \in {0, 1}$ kết hợp với biến số khe hở tách bánh $Z_L(t)$ cho phép mô tả chính xác lực pháp tuyến đơn phương và chuyển đổi liên tục cấu trúc ma trận của hệ phương trình vi phân.
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ ảnh hưởng tương hỗ giữa độ đàn nhớt của nền đường và dao động thẳng đứng của xe diễn ra như thế nào khi xuất hiện hiện tượng mất liên kết? Giả thuyết 2 (H2): Nền đường có hệ số độ cứng $k_S$ thấp làm gia tăng độ võng cục bộ $w_D(t)$, kéo dài tổng thời gian mất liên kết và kích hoạt hiện tượng va đập thứ cấp khi tiếp xúc trở lại.
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật phân bố áp suất (đều, parabol, tam giác) và kích thước vết tiếp xúc biến thiên ảnh hưởng như thế nào đến phản lực mặt đường so với giả định tải trọng tập trung? Giả thuyết 3 (H3): Phân bố áp suất liên tục trên diện tích tiếp xúc hữu hạn $A_j = d_{cj} \times b_L$ làm mượt phổ ứng suất tiếp xúc và thay đổi đáng kể giá trị cực đại của phản lực động so với mô hình điểm tiếp xúc tập trung.
  • Câu hỏi nghiên cứu 4 (RQ4): Sự khác biệt về đáp ứng động lực học giữa các mô hình phân cấp (1/4, 1/2 dọc, 1/2 ngang và không gian 3D) thể hiện ra sao trong cùng điều kiện kích thích mặt đường tiền định? Giả thuyết 4 (H4): Mô hình không gian phản ánh chính xác hiện tượng lệch tải và hiện tượng tách bánh bất đối xứng giữa các bánh xe mà các mô hình phẳng không thể dự báo.

Phạm vi nghiên cứu tập trung vào ô tô chuyển động thẳng đều với vận tốc danh nghĩa $V \in [10, 120] \text{ km/h}$, chịu kích thích động học tiền định (sóng sin nhiều chu kỳ, xung đơn parabol/cung tròn/hình thang, gờ giảm tốc tiêu chuẩn). Kết cấu mặt đường biến dạng được mô hình hóa bằng dầm Euler-Bernoulli (cho mô hình 1/4 và 1/2) và tấm mỏng Kirchhoff (cho mô hình không gian) đặt trên nền đàn nhớt Kelvin ($k_S, c_S$). Ý nghĩa khoa học của luận án thể hiện ở việc lượng hóa chính xác tổng thời gian mất liên kết, giá trị hiệu dụng gia tốc (RMS) và phản lực tiếp xúc cực đại, tạo cơ sở dữ liệu nền tảng cho việc thiết kế hệ thống treo chủ động và đánh giá tuổi thọ công trình đường bộ.

Literature Review và Positioning

Lịch sử nghiên cứu tương tác động lực học xe - đường được định hình bởi nhiều trường phái khoa học trên thế giới và tại Việt Nam. Phân tích tổng quan tài liệu cho thấy hai luồng quan điểm và phương pháp tiếp cận chính:

+-------------------------------------------------------------------------+
|                  TIẾN TRÌNH PHÁT TRIỂN CỦA Y VĂN QUỐC TẾ                |
+-------------------------------------------------------------------------+
|                                                                         |
|  [Trường phái Liên kết Song phương]    [Trường phái Liên kết Đơn phương]|
|  - Vũ Công Hàm & Nguyễn Đình Dũng      - Mahmoud Rababad & Atanu Bhuyan |
|    (2011, 2012): Mô hình 3D 7 DOF,       (2018): Cảnh báo tách bánh     |
|    bánh xe luôn bám đường.               dưới điều khiển Skyhook.       |
|  - Yang Shaopu et al. (2015): Xe -     - Cheng et al. (2001), Zhu (2014)|
|    đường tấm 2 lớp nền Kelvin,           Mô hình xe - dầm/cầu, tiếp xúc |
|    bỏ qua mất liên kết lốp.              đơn phương, thuật toán LCP.    |
|                                                                         |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|                ĐỊNH VỊ ĐỘT PHÁ CỦA LUẬN ÁN (PHÙNG MẠNH CƯỜNG, 2022)      |
|  - Thiết lập chuỗi mô hình phân cấp: 1/4 -> 1/2 Dọc -> 1/2 Ngang -> 3D. |
|  - Kết hợp đồng thời: Nền Kelvin + Mất liên kết + Vết tiếp xúc 2D.      |
|  - Khử phương trình đạo hàm riêng bằng Galerkin giải số trên MATLAB.     |
+-------------------------------------------------------------------------+

Trường phái thứ nhất tập trung vào mô hình hóa dao động ô tô nhiều bậc tự do nhưng giữ nguyên giả định liên kết song phương lý tưởng. Vũ Đức Lập (2006) trong giáo trình kinh điển về dao động ô tô đã phân tích định tính nguy cơ mất an toàn khi bánh xe tách khỏi mặt đường, song chưa lượng hóa được động học tiếp xúc phi tuyến. Nguyễn Đình Dũng (2014) phát triển mô hình không gian 7 bậc tự do kết hợp dầm/tấm đàn hồi trên nền Kelvin chịu kích thích tiền định; tuy nhiên, mô hình vẫn giả thiết bánh xe luôn tiếp xúc liên tục với mặt đường. Nhóm tác giả Vũ Công Hàm và Nguyễn Đình Dũng (2011, 2012) khảo sát lực tương tác xe - đường và ghi nhận dấu của lực liên kết bị đảo chiều, qua đó kết luận điều kiện tách bánh xảy ra rất phổ biến nhưng chưa thiết lập hệ phương trình vi phân chuyển động cho giai đoạn tách rời.

Trường phái thứ hai tiếp cận bài toán tiếp xúc đơn phương (unilateral contact) và tương tác xe - kết cấu công trình. Cheng cùng các cộng sự (2001) khảo sát dao động dầm cầu liên tục dưới tác dụng của dòng phương tiện chuyển động có xét hiện tượng mất liên kết, giải hệ phương trình tách biến từng chặng nhưng bỏ qua tính đàn hồi và giảm chấn của lốp xe. Zhu và các cộng sự (2014) áp dụng bài toán bù tuyến tính (Linear Complementarity Problem - LCP) kết hợp phương pháp phần tử hữu hạn (FEM) để mô phỏng sự tách bánh và va chạm của xe khi qua mố cầu. Gamaleddine Elnashar, Rama B. Bhat và Ramin Sedaghati (2014) khảo sát mô hình 1/4 xe trên dầm Euler-Bernoulli nền Kelvin chịu kích thích điều hòa bằng phương pháp Galerkin, chỉ ra độ cứng nền đường thấp tác động mạnh đến gia tốc thân xe, đồng thời khuyến nghị phải tích hợp hiện tượng mất liên kết và hình học vết tiếp xúc. Yang Shaopu cùng cộng sự (2015) nghiên cứu hệ xe - đường dạng tấm hai lớp trên nền đàn nhớt phi tuyến nhưng vẫn dựa trên mô hình tiếp xúc điểm.

Sự đối nghịch giữa hai luồng nghiên cứu nằm ở việc cân bằng giữa tính xác thực vật lý của liên kết bánh xe và độ phức tạp tính toán của bài toán biên chuyển động. Các nghiên cứu quốc tế về tiếp xúc đơn phương thường giới hạn ở mô hình phẳng 1/4 hoặc kết cấu dầm cầu nhịp đơn. Luận án của Phùng Mạnh Cường đã định vị chính xác khoảng trống này: Thiết lập một hệ thống hoàn chỉnh từ mô hình 1/4, mô hình 1/2 phẳng (dọc và ngang) đến mô hình không gian 3D của hệ xe - đường kết hợp có tính đến biến dạng nền móng, hiện tượng mất liên kết và phân bố áp suất diện tích tiếp xúc. Đây là bước tiến vượt bậc so với công trình của Gamaleddine Elnashar (2014) và Yang Shaopu (2015).

Đó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 và làm sâu sắc lý thuyết động lực học công trình và cơ học tiếp xúc phi tuyến thông qua việc hoàn thiện hệ phương trình trạng thái mô tả hệ thống xe - đường. Công trình thách thức giả định truyền thống về tính liên tục của liên kết lốp - đường vốn được áp dụng phổ biến trong các tiêu chuẩn thiết kế hệ thống treo và mặt đường giao thông.

       +-------------------------------------------------------------+
       |   QUY TRÌNH BIẾN ĐỔI TOÁN TỬ VÀ RỜI RẠC HÓA HỆ PHƯƠNG TRÌNH  |
       +-------------------------------------------------------------+
                                      |
                                      v
+----------------------------------------------------------------------------+
| [PTVP Đạo hàm riêng của Tấm Kirchhoff]:                                    |
|   D \nabla^4 w(x,y,t) + \rho h_P \ddot{w} + c_S \dot{w} + k_S w = p(x,y,t) |
+----------------------------------------------------------------------------+
                                      |
                                      v [Phép khai triển Fourier kép Galerkin]
+----------------------------------------------------------------------------+
|   w(x,y,t) = \sum_{m=1}^\infty \sum_{n=1}^\infty T_{mn}(t) W_{mn}(x,y)     |
|   với W_{mn}(x,y) = \sin(m\pi x / L_P) \sin(n\pi y / B_P)                  |
+----------------------------------------------------------------------------+
                                      |
                                      v [Chuyển tổng 2 chỉ số -> 1 chỉ số k]
+----------------------------------------------------------------------------+
| [Hệ PTVP Vi phân thường ma trận đồng nhất]:                                |
|   [M] \ddot{q}(t) + [C] \dot{q}(t) + [K] q(t) = F(t, q, \dot{q})           |
+----------------------------------------------------------------------------+

Khung lý thuyết của luận án thiết lập mệnh đề toán - cơ học về liên kết đơn phương gián đoạn:

  • Mệnh đề 1 (Trạng thái liên kết bám đường, $s_j(t) = 1$): Lực tiếp xúc pháp tuyến $R_j(t) > 0$, khoảng cách tách bánh $Z_{Lj}(t) = 0$. Phản lực tiếp xúc được tính bởi phương trình tích phân diện tích: $$R_j(t) = \iint_{A_j} p_j(x, y, t) , dA = k_{Lj} z_{Lj}(t) + c_{Lj} \dot{z}{Lj}(t)$$ trong đó $z{Lj}(t) = u_{cj}(t) - w_{Dj}(t) - r_{Dj}(t)$ là độ biến dạng hướng tâm của lốp xe thứ $j$.
  • Mệnh đề 2 (Trạng thái mất liên kết tách bánh, $s_j(t) = 0$): Khi biến dạng của lốp triệt tiêu và chuyển sang trạng thái kéo danh nghĩa ($F_{Lj} \le 0$), bánh xe tách khỏi nền đường: $$R_j(t) = 0 \quad \text{và} \quad Z_{Lj}(t) = r_{Dj}(t) + w_{Dj}(t) - u_{cj}(t) > 0$$

Mô hình lý thuyết này chứng minh rằng hiện tượng tách bánh không đơn thuần là sự gián đoạn lực cơ học, mà tạo ra sự thay đổi tức thời cấu trúc không gian trạng thái của cơ hệ nhiều vật (multibody system), làm suy giảm đột ngột ma trận cản và độ cứng tiếp xúc của toàn hệ thống.

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

Khung phân tích của luận án là sự hợp nhất đồng bộ của ba nền tảng lý thuyết:

  1. Lý thuyết dầm Euler-Bernoulli và tấm mỏng Kirchhoff: Biểu diễn kết cấu mặt đường dưới dạng phần tử liên tục có xét quán tính khối lượng $\rho$, chiều dày $h_P$, độ cứng uốn $D = E h_P^3 / [12(1-\nu^2)]$.
  2. Lý thuyết nền móng đàn nhớt Kelvin: Mô hình hóa lớp đệm đất nền thông qua hệ lò xo độ cứng $k_S$ phân bố kết hợp phần tử cản nhớt $c_S$ tiêu tán năng lượng.
  3. Cơ học tiếp xúc động lực học: Khảo sát diện tích tiếp xúc chữ nhật $A_j = d_c \times b_L$ với 4 quy luật phân bố áp suất: phân bố đều (Uniform), phân bố hình sin (Sinusoidal), phân bố parabol bậc hai (Parabolic) và phân bố tam giác (Triangular).

Điều kiện biên xác định (boundary conditions) được chuẩn hóa: Dầm và tấm chịu liên kết tựa đơn (simply supported) trên toàn bộ chu vi, phản ánh sát thực điều kiện làm việc của các tấm mặt đường bê tông xi măng có khe co giãn hoặc kết cấu dầm cầu tiêu chuẩ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ủ nghiêm ngặt thế giới quan thực chứng (positivism), kết hợp chặt chẽ giữa giải tích cơ học giải pháp giải tích (analytical formulation) và mô phỏng số học trên hệ thống phân cấp 4 cấp độ mô hình:

[Mô hình 1/4 Xe] --------> [Mô hình 1/2 Dọc] --------> [Mô hình 1/2 Ngang] -------> [Mô hình 3D Không gian]
  (2 DOF + 1D Beam)          (4 DOF + 1D Beam)           (4 DOF + 1D Beam)            (7 DOF + 2D Plate)
  - Khảo sát đáp ứng         - Dao động nhấp nhô         - Dao động lắc ngang         - Khảo sát tương tác
    cơ bản, biến dạng          và lắc dọc thân xe,         thân xe và cầu,              toàn diện 4 bánh,
    vết lốp, loại xung.        lệch pha trước/sau.         nguy cơ lật bên.             lệch pha đa chiều.

Thiết kế đa tầng này cho phép cô lập từng thành phần chuyển động (nhấp nhô, lắc dọc pitch, lắc ngang roll) trước khi tổng hợp vào không gian trạng thái 3D toàn phần 7 bậc tự do (gồm chuyển vị đứng $u_b$, góc lắc dọc $\varphi_b$, góc lắc ngang $\theta_b$ của thân xe; chuyển vị đứng $u_{c1}, u_{c2}$ và góc lắc ngang $\theta_{c1}, \theta_{c2}$ của hai cầu xe).

Quy trình nghiên cứu rigorous

Quy trình giải quyết hệ phương trình vi phân đạo hàm riêng phi tuyến được tiến hành theo các bước chuẩn hóa:

  • Bước 1: Thiết lập hệ phương trình vi phân liên tục: Biểu diễn chuyển động của hệ xe và phương trình dao động của dầm/tấm chịu tải trọng di động $p(x,y,t)$.
  • Bước 2: Rời rạc hóa không gian bằng phương pháp Galerkin: Khai triển chuyển vị thẳng đứng của mặt đường theo chuỗi hàm trực giao: $$w(x, y, t) = \sum_{m=1}^\infty \sum_{n=1}^\infty T_{mn}(t) \sin\left(\frac{m\pi x}{L_P}\right) \sin\left(\frac{n\pi y}{B_P}\right)$$ với $L_P, B_P$ lần lượt là chiều dài và chiều rộng tính toán của tấm đường.
  • Bước 3: Chuyển đổi chỉ số tổng kép: Sử dụng thuật toán ánh xạ song ánh quy đổi tổng kép $(m, n)$ về tổng đơn $k = 1 \div N_{\text{mod}}$ nhằm tối ưu hóa cấu trúc ma trận khối lượng $[M]$, ma trận cản $[C]$ và ma trận độ cứng $[K]$.
  • Bước 4: Tích hợp thuật toán điều kiện tiếp xúc: Tại mỗi bước thời gian $\Delta t$, thuật toán kiểm tra đồng thời điều kiện lực $F_{Lj}(t)$ và điều kiện hình học $Z_{Lj}(t)$. Nếu $F_{Lj} > 0 \Rightarrow s_j = 1$; nếu $F_{Lj} \le 0 \Rightarrow s_j = 0$, ngắt liên kết đàn nhớt của lốp khỏi mặt đường.
  • Bước 5: Giải số tích phân phi tuyến: Lập trình giải hệ PTVP thường trên môi trường MATLAB bằng các thuật toán biến đổi bước lặp thích nghi (ODE suite: ode45, ode15s chuyên dụng cho hệ stiff).
          +----------------------------------------------------+
          |      THUẬT TOÁN ĐIỀU KIỆN TIẾP XÚC ĐƠN PHƯƠNG      |
          +----------------------------------------------------+
                                   |
                                   v
                      +--------------------------+
                      |  Tính toán bước thời gian|
                      |     t = t + \Delta t     |
                      +--------------------------+
                                   |
                                   v
          +--------------------------------------------------+
          | Kiểm tra điều kiện lực lò xo - giảm chấn bánh xe |
          |      F_{Lj}(t) = k_{Lj} z_{Lj} + c_{Lj} \dot{z}   |
          +--------------------------------------------------+
                                   |
                 +-----------------+-----------------+
                 |                                   |
                 v [F_{Lj} > 0]                      v [F_{Lj} <= 0]
      +---------------------+             +---------------------+
      |  TRẠNG THÁI BÁM     |             |  TRẠNG THÁI TÁCH    |
      | - s_j(t) = 1        |             | - s_j(t) = 0        |
      | - R_j(t) = F_{Lj}(t)|             | - R_j(t) = 0        |
      | - Z_{Lj}(t) = 0     |             | - Z_{Lj}(t) > 0     |
      | - Đường chịu tải p_j|             | - Đường tự do p_j=0 |
      +---------------------+             +---------------------+
                 |                                   |
                 +-----------------+-----------------+
                                   |
                                   v
          +--------------------------------------------------+
          | Cập nhật ma trận [M], [C], [K] và giải số bước   |
          | tiếp theo bằng phương pháp Runge-Kutta thích nghi|
          +--------------------------------------------------+

Data và phân tích

Toàn bộ tham số vật lý của cơ hệ được thu thập từ các xe tham chiếu thực tế (xe tải 2 cầu và xe khách tiêu chuẩn):

  • Khối lượng thân xe $m_b = 8500 \text{ kg}$, mômen quán tính $J_{bx} = 2800 \text{ kg}\cdot\text{m}^2$, $J_{by} = 12500 \text{ kg}\cdot\text{m}^2$.
  • Khối lượng cầu trước $m_{c1} = 650 \text{ kg}$, cầu sau $m_{c2} = 1100 \text{ kg}$.
  • Độ cứng hệ treo trước $k_{Tf} = 220 \text{ kN/m}$, hệ treo sau $k_{Tr} = 380 \text{ kN/m}$.
  • Độ cứng lốp xe $k_{Lf} = 850 \text{ kN/m}$, $k_{Lr} = 1700 \text{ kN/m}$.
  • Thông số đường dầm/tấm: Môđun đàn hồi bê tông $E = 3.1 \times 10^{10} \text{ N/m}^2$, khối lượng riêng $\rho = 2500 \text{ kg/m}^3$, hệ số Poisson $\nu = 0.18$, chiều dày tấm $h_P = 0.25 \text{ m}$. Độ cứng nền Kelvin khảo sát trong dải $k_S = 1.0 \times 10^7 \div 1.0 \times 10^8 \text{ N/m}^3$.
  • Kích thích mặt đường: Xung đơn chiều cao $h_E = 0.05 \div 0.08 \text{ m}$, chiều dài $L_E = 1.0 \div 3.0 \text{ m}$; sóng sin liên tục bước sóng $L_E = 5 \text{ m}$.

Chỉ tiêu đánh giá thống kê tập trung vào: Giá trị cực đại ($\text{Max}$) và giá trị hiệu dụng ($\text{RMS}$) của gia tốc thẳng đứng thân xe $\ddot{u}_b$, gia tốc cầu $\ddot{u}c$, lực tiếp xúc $R_j(t)$ và tổng thời gian mất liên kết tích lũy $\sum \Delta t{\text{MLK}}$.

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

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

Trích dẫn trực tiếp từ văn bản luận án xác nhận tính cấp bách của hiện tượng:

"Sự mất tiếp xúc của bánh xe với mặt đường sẽ làm mất tính điều khiển (phanh, lái, tốc độ), có thể gây mất an toàn cho người, hàng hóa và chính phương tiện." (Trang 2) "Chu trình 'dao động - phá hỏng - dao động - v.v.' cứ lặp đi lặp lại với mức độ ngày càng tăng." (Trang 2, 7) "Xét trên phương diện cơ học thì mô hình dao động của ô tô có kể đến hiện tượng tách bánh là mô hình của cơ hệ có liên kết thay đổi theo thời gian. Do liên kết là một bộ phận cấu thành của mô hình nên sự thay đổi của liên kết kéo theo sự thay đổi của mô hình." (Trang 28)

Dữ liệu mô phỏng số đã mang lại 4 phát hiện đột phá:

Tham số / Hiện tượng Kết quả khi KHÔNG tính MLK Kết quả khi CÓ tính MLK (Phát hiện của Luận án) Ý nghĩa cơ học & An toàn
Phản lực tiếp xúc cực tiểu $R_{\min}$ Xuất hiện lực âm giả tạo ($R < 0$, lực kéo dính) $R(t) \equiv 0$ khi tách bánh ($Z_L > 0$) Phản ánh chính xác bản chất cơ học tiếp xúc đơn phương.
Phản lực tiếp xúc cực đại $R_{\max}$ Bị đánh giá thấp hơn thực tế từ $22% - 38%$ Tăng vọt tức thời khi tiếp xúc trở lại (va đập) Làm rõ nguyên nhân gây phá hủy mỏi nhanh mặt đường và ổ trục.
Gia tốc thân xe $\text{RMS}(\ddot{u}_b)$ Tăng tuyến tính đều theo vận tốc Biến thiên phi tuyến, xuất hiện các đỉnh cộng hưởng cục bộ Đánh giá chính xác mức độ suy giảm độ êm dịu theo tiêu chuẩn ISO 2631-1.
Ảnh hưởng của vận tốc $V$ Bánh xe luôn giữ tiếp xúc Vượt ngưỡng $V_{cr} \approx 45 \text{ km/h}$, thời gian MLK tăng vọt Cảnh báo mất kiểm soát lái và phanh ở dải tốc độ cao qua mấp mô. 
    Lực tiếp xúc R(t) [kN]
       ^
       |        /\               /\  (Va đập khi tiếp xúc lại)
       |       /  \             /  \
  Tĩnh + - - -/- - \- - - - - -/- - \- - - - - - - - - - - - -
       |     /      \         /      \
     0 +----+--------+-------+--------+------------------------> Thời gian t [s]
       |              \_____/          <--- VÙNG MẤT LIÊN KẾT (MLK)
       |          (Lực âm giả tạo          [Luận án: R = 0, Khe hở Z_L > 0]
       |         trong mô hình cũ)
  1. Hiệu ứng va đập tái tiếp xúc (Re-contact impact phenomenon): Trong mô hình tuyến tính cũ, khi bánh xe đi qua gờ mấp mô lớn, lực liên kết tính toán bị âm (mặt đường "hút" bánh xe). Khi đưa điều kiện đơn phương vào, lực tiếp xúc triệt tiêu ($R=0$), bánh xe chuyển động tự do trong không gian dưới tác dụng quán tính và lực bung của lò xo treo. Khi rơi trở lại mặt đường, hiện tượng va đập xuất hiện làm phản lực động cực đại $R_{\max}$ tăng vọt từ $1.35$ đến $1.85$ lần so với tải trọng tĩnh, giải thích hiện tượng nứt vỡ kết cấu bê tông mặt đường tại các vị trí sau gờ giảm tốc.
  2. Tác động phi tuyến của độ cứng nền đường $k_S$: Khi $k_S$ giảm (nền đất yếu), độ võng động $w_D(t)$ của đường tăng lên, kéo theo sự gia tăng của khe hở tách bánh $Z_L(t)$. Thời gian mất liên kết tại cầu sau tăng đến $32%$ khi xe chạy qua mấp mô dạng sóng hình sin trên nền đất yếu so với nền đường cứng.
  3. Sự nhạy cảm với quy luật phân bố áp suất: Quy luật phân bố parabol và phân bố đều trên diện tích tiếp xúc $d_c \times b_L$ cho giá trị gia tốc RMS thân xe thấp hơn từ $4.5% - 7.2%$ so với mô hình tập trung điểm, chứng minh tính năng lọc mấp mô hình học của lốp xe thực tế.
  4. Hiện tượng tách bánh bất đối xứng trong không gian 3D: Khác với mô hình 1/2 dọc dự báo hai bánh trên cùng một cầu tách bánh đồng thời, mô hình không gian 3D chỉ ra rằng dưới kích thích mấp mô lệch pha giữa hai vệt bánh trái - phải, chuyển vị góc lắc ngang $\theta_b$ kết hợp với chuyển vị góc cầu $\theta_c$ làm một bánh xe bị tách hoàn toàn trong khi bánh đối diện chịu tải trọng tăng gấp đôi, làm tăng nguy cơ lật xe (rollover risk).

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp khung giải tích toán - cơ học chuẩn mực cho hệ đa vật thể có liên kết biến đổi theo thời gian, đặt nền móng cho việc mô hình hóa các hệ thống động lực phức tạp khác như tàu cao tốc - ray hay máy bay khi hạ cánh trên đường băng mấp mô.
  • Về mặt phương pháp luận: Khẳng định sự bắt buộc phải áp dụng quy trình biến đổi Galerkin kết hợp giải thuật kiểm soát tiếp xúc phi tuyến khi phân tích động lực học kết cấu giao thông.
  • Về mặt ứng dụng công nghiệp ô tô: Cung cấp thông số ngưỡng tới hạn (critical velocity thresholds) cho các nhà sản xuất xe tải và xe quân sự để tối ưu hóa đặc tính phi tuyến của giảm chấn thủy lực, ngăn chặn hiện tượng tách bánh kéo dài.
  • Về mặt quy hoạch và bảo trì cầu đường: Cung cấp công cụ tính toán chính xác hệ số tải trọng động (Impact Factor) phục vụ thiết kế kết cấu áo đường mềm và mặt đường bê tông xi măng, giúp tối ưu hóa chu kỳ duy tu sửa chữa đường bộ.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn chỉ rõ 4 giới hạn nội tại:

  1. Giả thiết kích thích mặt đường thuần túy tiền định (deterministic profiles) chưa bao quát hết tính chất ngẫu nhiên bất quy tắc (stochastic stationary process) của biên dạng mặt đường thực tế theo phân loại ISO 8608 (cấp mặt đường A, B, C, D, E).
  2. Các cụm lò xo và giảm chấn trong hệ thống treo và lốp xe được giả thiết có đặc tính đàn hồi và cản tuyến tính, chưa xét đến tính phi tuyến phi đối xứng (hành trình nén khác hành trình trả) của giảm chấn thủy lực thực tế.
  3. Bỏ qua sự mất mát động năng và vận tốc do va chạm cơ học tại thời điểm bánh xe tái lập liên kết với mặt đường (hệ số phục hồi va chạm $e = 1$).
  4. Khảo sát giới hạn ở chế độ chuyển động thẳng đều với vận tốc không đổi ($V = \text{const}$), chưa tích hợp quá trình phanh khẩn cấp, tăng tốc đột ngột hay quay vòng chuyển làn.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi chiến lược:

  • Hướng 1: Mở rộng mô hình không gian sang kích thích ngẫu nhiên sử dụng phương pháp mật độ phổ công suất (Power Spectral Density - PSD) và hàm truyền phi tuyến.
  • Hướng 2: Tích hợp mô hình lốp phi tuyến nâng cao (như Magic Formula của Pacejka hoặc mô hình lốp chùm dải đàn hồi Fiala) kết hợp tổn hao va chạm phi đàn hồi.
  • Hướng 3: Khảo sát bài toán điều khiển chủ động (Active Suspension Control) dựa trên thuật toán điều khiển thích nghi $\mathcal{H}_\infty$ và MPC nhằm chủ động triệt tiêu hiện tượng mất liên kết.
  • Hướng 4: Tiến hành thực nghiệm kiểm chứng quy mô đầy đủ (Full-scale experimental validation) sử dụng xe thí nghiệm chuyên dụng gắn cảm biến đo áp lực bánh xe không dây và hệ thống quét laser mặt đường 3D.

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

Luận án của Phùng Mạnh Cường tạo lập giá trị học thuật và thực tiễn sâu rộng:

  • Tác động học thuật: Dự báo đóng góp trích dẫn quan trọng cho các nghiên cứu thuộc nhóm ngành Kỹ thuật Giao thông, Cơ học Động lực học và Kỹ thuật Cơ sở hạ tầng trong nước và quốc tế. Khung toán tử của luận án thiết lập chuẩn mực mới cho các đề tài nghiên cứu sinh về tương tác phương tiện - kết cấu công trình.
  • Chuyển đổi công nghiệp ô tô & quốc phòng: Đóng góp trực tiếp cho các viện nghiên cứu thiết kế xe chuyên dụng quân sự (thuộc Bộ Quốc phòng) và các trung tâm R&D sản xuất ô tô dân dụng trong việc tối ưu hóa hệ thống treo chịu tải cao di chuyển trên địa hình gồ ghề phức tạp.
  • Chính sách & Tiêu chuẩn giao thông: Cung cấp cơ sở khoa học định lượng để Bộ Giao thông Vận tải và Cục Đường bộ chuẩn hóa quy chuẩn kích thước hình học của gờ giảm tốc trên các tuyến quốc lộ, hạn chế tối đa các xung lực va đập nguy hiểm cho phương tiện.

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

  • Nghiên cứu sinh & Nhà khoa học chuyên sâu: Thụ hưởng phương pháp luận giải tích chuyển đổi PDE sang ODE bằng phương pháp Galerkin và thuật toán xử lý liên kết biến đổi theo thời gian.
  • Kỹ sư R&D công nghiệp ô tô: Sử dụng bộ chương trình tính toán và mô phỏng trên MATLAB làm công cụ kiểm nghiệm ảo các thông số độ cứng $k_T, k_L$ và hệ số cản $c_T, c_L$ trước khi chế tạo mẫu thử nghiệm vật lý.
  • Kỹ sư thiết kế công trình đường và cầu: Ứng dụng mô hình tương tác tấm trên nền Kelvin có xét lực va đập động để tính toán chính xác chiều dày áo đường và cốt thép gia cường mặt đường sân bay, bến bãi hạng nặng.
  • Cơ quan quản lý an toàn giao thông: Tiếp cận các dữ liệu thực nghiệm mô phỏng về tốc độ giới hạn an toàn khi xe lưu thông qua các đoạn đường có chất lượng bề mặt xấu.

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ì?

Đóng góp độc đáo nhất là việc thiết lập thành công mô hình toán học giải tích hợp nhất cho cơ hệ xe - đường kết hợp có liên kết đơn phương gián đoạn (unilateral intermittent contact). Luận án đã mở rộng lý thuyết tấm Kirchhoff trên nền đàn nhớt Kelvin bằng cách nhúng trực tiếp biến trạng thái logic $s(t)$ và khe hở động $Z_L(t)$ vào hệ phương trình vi phân dao động không gian 3D, khắc phục triệt để nghịch lý "lực tiếp xúc âm" của cơ học cổ điển khi nghiên cứu ô tô qua mấp mô.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu trước đây?

So sánh với nghiên cứu của Vũ Công Hàm & Nguyễn Đình Dũng (2011) vốn chỉ dừng lại ở việc phát hiện sự đổi dấu của lực tiếp xúc trên mô hình liên kết cứng, và nghiên cứu của Gamaleddine Elnashar et al. (2014) chỉ giải bài toán 1/4 xe trên dầm phẳng, luận án của Phùng Mạnh Cường đã thực hiện bước nhảy vọt:

  • Mở rộng từ mô hình dầm 1D sang mô hình tấm 2D không gian toàn phần.
  • Chuyển đổi thành công hệ phương trình đạo hàm riêng liên tục phi tuyến thành hệ phương trình vi phân thường ma trận nhờ tích hợp phép khai triển Galerkin với thuật toán ánh xạ chỉ số ma trận đồng nhất, giải trực tiếp trên MATLAB mà không cần đơn giản hóa liên kết.

3. Phát hiện bất ngờ nhất kèm minh chứng số liệu là gì?

Phát hiện bất ngờ nhất là sự suy giảm đột ngột và phi tuyến của tính ổn định khi xe đạt dải vận tốc $V = 50 \div 70 \text{ km/h}$ qua mấp mô dạng sóng sin nhiều chu kỳ. Tại dải tốc độ này, tần số kích thích từ mặt đường tiệm cận tần số dao động riêng của khối lượng không được treo (cầu xe), khiến thời gian mất liên kết chiếm tới $18% - 25%$ tổng thời gian chuyển động qua đoạn mấp mô. Đồng thời, lực va đập khi liên kết trở lại tạo ra đỉnh phản lực gấp $1.85$ lần tải trọng tĩnh, gây mất kiểm soát lái hoàn toàn ở mô hình không gian 3D trong khi mô hình 1/4 chỉ dự báo mức tăng $1.25$ lần.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập nghiên cứu:

  • Hệ thống ký hiệu toán học nhất quán (danh mục 42 ký hiệu La-tinh và Hy Lạp).
  • Tường minh hóa toàn bộ các ma trận khối lượng $[M]$, ma trận cản $[C]$, ma trận độ cứng $[K]$ và véc-tơ lực kích thích ${F}$ cho cả 4 mô hình (từ 1/4 đến 3D không gian).
  • Quy trình giải thuật toán rời rạc hóa Galerkin và điều kiện chuyển đổi trạng thái tiếp xúc $s(t)$ được trình bày dưới dạng sơ đồ khối logic, sẵn sàng cho việc lập trình độc lập trên bất kỳ ngôn ngữ tính toán số nào (MATLAB, Python, C++).

5. Lộ trình nghiên cứu 10 năm (10-Year Research Agenda) được phác thảo ra sao?

Lộ trình phát triển học thuật giai đoạn 2022–2032 được xác lập qua 3 giai đoạn:

  • Giai đoạn 1 (2022–2025): Hoàn thiện mô hình tiếp xúc đơn phương cho hệ thống treo phi tuyến (phi tuyến hình học và phi tuyến giảm chấn thủy lực) dưới kích thích ngẫu nhiên mặt đường theo tiêu chuẩn ISO 8608.
  • Giai đoạn 2 (2025–2028): Tích hợp động lực học phanh ngặt (ABS/EBD) và hệ thống lái chủ động trong điều kiện mất liên kết cục bộ; phát triển phần mềm mô phỏng đồng thời (co-simulation) giữa động lực học hệ nhiều vật và cơ học phá hủy áo đường.
  • Giai đoạn 3 (2028–2032): Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) và mạng nơ-ron vật lý (PINNs) để dự báo thời gian thực hiện tượng mất liên kết và điều khiển chủ động hệ thống treo trên các phương tiện tự hành (Autonomous Vehicles) di chuyển trên địa hình tự nhiên phức tạp.

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Phùng Mạnh Cường đã khẳng định vị thế một công trình nghiên cứu hàn lâm xuất sắc trong lĩnh vực Cơ kỹ thuật và Động lực học ô tô. Những đóng góp then chốt được đúc kết toàn diện:

  1. Xây dựng hoàn chỉnh hệ thống 4 mô hình dao động phân cấp (1/4, 1/2 dọc, 1/2 ngang và không gian 3D) của hệ xe - đường có tính đến đồng thời biến dạng nền đàn nhớt Kelvin, hiện tượng mất liên kết đơn phương và sự biến thiên diện tích vết tiếp xúc.
  2. Thiết lập thuật toán giải tích kết hợp phương pháp số học Galerkin giải quyết triệt để hệ phương trình vi phân đạo hàm riêng phi tuyến mô tả tương tác tấm mỏng Kirchhoff và thân xe nhiều bậc tự do.
  3. Lượng hóa chính xác cơ chế cơ học của hiện tượng va đập tái tiếp xúc, chứng minh phản lực động cực đại tăng từ $22% - 38%$ so với các mô hình bỏ qua hiện tượng tách bánh.
  4. Xác lập mối quan hệ định lượng giữa vận tốc chuyển động, độ cứng nền đường $k_S$, quy luật phân bố áp suất vết lốp với tổng thời gian mất liên kết và độ êm dịu chuyển động theo chuẩn ISO 2631-1.
  5. Cung cấp bộ công cụ mô phỏng số học tin cậy trên nền tảng MATLAB, phục vụ trực tiếp công tác nghiên cứu, đào tạo sau đại học và ứng dụng thiết kế kỹ thuật trong công nghiệp ô tô và công trình giao thông tại Việt Nam.