Tổng quan về luận án

Ngành vận tải hàng hóa đường bộ và chuỗi cung ứng logistics tại Việt Nam đang chứng kiến sự tăng trưởng vượt bậc về nhu cầu vận chuyển siêu trường, siêu trọng. Trong cơ cấu phương tiện vận tải thương mại, sơ mi rơ moóc (SMRM) kết hợp cùng xe đầu kéo (XĐK) đóng vai trò xương sống nhờ tính năng linh hoạt, khả năng chuyên chở container tải trọng lớn và hiệu suất luân chuyển cao. Tuy nhiên, thực trạng sản xuất SMRM trong nước đối mặt với nhiều rào cản kỹ thuật: quy trình thiết kế phần lớn mang tính kinh nghiệm, mô phỏng tĩnh đơn thuần, dẫn đến việc khung xe thường có tự trọng bản thân lớn nhưng tuổi thọ làm việc lại không đồng đều, dễ xuất hiện các vết nứt mỏi nghiêm trọng tại các vùng liên kết dầm và hệ thống treo sau thời gian ngắn khai thác.

Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở sự thiếu hụt một phương pháp luận hoàn chỉnh cho phép dự báo ứng xử ứng suất động và tuổi thọ mỏi của kết cấu khung SMRM trong điều kiện tương tác thực tế với mặt đường ngẫu nhiên tại Việt Nam. Theo các nghiên cứu kinh điển, khung chassis thường được phân tích dựa trên phương pháp sức bền vật liệu truyền thống với hệ số tải trọng động $K_d$ giả định hoặc mô hình hóa bằng phương pháp phần tử hữu hạn (FEM) tĩnh với các điều kiện biên ngàm cứng lý tưởng hóa. Cách tiếp cận này loại bỏ tính phi tuyến của động lực học hệ thống treo và bỏ qua sự phân bố tải trọng động biến thiên liên tục theo thời gian, dẫn đến sai số lớn khi đánh giá độ bền mỏi chu kỳ cao (High-Cycle Fatigue - HCF).

Để giải quyết khoảng trống này, luận án tập trung vào các câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) then chốt:

  • RQ1: Làm thế nào để mô tả chính xác tương tác động lực học phi tuyến giữa đoàn xe 5 trục (tổ hợp XĐK và SMRM 40 feet) với các mấp mô mặt đường tiêu chuẩn theo thời gian thực?
  • RQ2: Quy luật phân bố ứng suất động và tích lũy tổn thương mỏi trên kết cấu khung chassis đàn hồi diễn ra như thế nào dưới tác động của các dải vận tốc và cấp mặt đường khác nhau?
  • RQ3: Phương án cải tiến hình học và vật liệu nào giúp tối ưu hóa khối lượng khung xe mà vẫn đảm bảo tiêu chuẩn độ bền tĩnh và gia tăng tuổi thọ mỏi?
  • H1: Việc tích hợp đồng mô phỏng động lực học hệ nhiều vật (MBD) với mô hình phần tử hữu hạn (FEM) thông qua kỹ thuật tổ hợp dao động riêng thành phần (Craig-Bampton Component Mode Synthesis) sẽ phản ánh chính xác trạng thái ứng suất động thực tế hơn phương pháp giải tích tĩnh từ 25% đến 40%.
  • H2: Tối ưu hóa cục bộ thông qua tạo lỗ giảm tải hình học trên dầm ngang và sử dụng thép cường độ cao cho dầm chính sẽ giảm ít nhất 10% khối lượng tự trọng mà không làm suy giảm giới hạn bền mỏi của toàn bộ kết cấu.

Khung lý thuyết của công trình được định vị trên nền tảng cơ học hệ nhiều vật (Lagrange DAEs), cơ học vật rắn biến dạng đàn hồi (Elasticity & Dynamic FEM), và lý thuyết bền mỏi kim loại tích lũy tuyến tính (Palmgren-Miner Rule, biểu đồ Haigh kết hợp đường cong S-N hiệu chỉnh Goodman/Gerber).

Đóng góp đột phá của luận án là việc xây dựng thành công quy trình khép kín: từ thiết lập mô hình đồng mô phỏng MBD-FEM toàn diện cho đoàn xe SMRM 40 feet 2 trục chở container, trích xuất tải trọng động thời gian thực, dự báo tuổi thọ mỏi bằng thuật toán Rainflow, cho đến kiểm chứng thực nghiệm bằng hệ thống đo động lực học và biến dạng đa kênh hiện đại. Phạm vi nghiên cứu bao quát mẫu SMRM tiêu chuẩn sản xuất từ thép Q345B, khảo sát trên các dải vận tốc từ 20 km/h đến 100 km/h dưới kích động mặt đường mấp mô đơn hình sin và mặt đường ngẫu nhiên cấp C-D, D-E theo tiêu chuẩn ISO 8608.

Literature Review và Positioning

Tổng quan tài liệu quốc tế và trong nước cho thấy sự tiến hóa rõ rệt trong các trường phái nghiên cứu độ bền khung gầm ô tô tải nặng:

Trường phái giải tích và cơ học kết cấu cổ điển (Timoshenko, 1953; Wong, 2001) tiếp cận bài toán bằng cách xem khung xe như hệ dầm phẳng hoặc dầm không gian chịu uốn và xoắn tĩnh, tải trọng động được quy đổi thông qua hệ số kinh nghiệm $K_d = 1.5 \div 2.5$. Cách làm này phục vụ tốt giai đoạn tiền thiết kế nhưng hoàn toàn bất lực trong việc xác định các điểm tập trung ứng suất cục bộ tại các mối hàn và góc lượn hình học.

Trường phái phân tích phần tử hữu hạn độc lập (FEM) phát triển mạnh mẽ từ những năm 2000. Các nghiên cứu của Caban (2018), Moaaz và cộng sự (2020) đã chứng minh rằng hiện tượng nứt gãy khung xe tải nặng trong thực tế thường bắt nguồn từ vùng tập trung ứng suất quanh mối hàn giữa dầm dọc (chính) và dầm ngang gia cường. Mi Chengji và cộng sự (2016) thực hiện phân tích bán tĩnh (quasi-static) xác định điểm yếu mỏi của khung xe tải tại khớp treo sau (nút số 229856) với chu kỳ mỏi giới hạn $1.03 \times 10^6$ chu kỳ. Tuy nhiên, hạn chế lớn nhất của trường phái này là việc đặt điều kiện biên ngàm cưỡng bức tại các vị trí mỏ nhíp hoặc quang treo, làm sai lệch đáng kể ma trận độ cứng động học của toàn hệ thống.

Trường phái mô phỏng kết hợp MBD-FEM và phân tích mỏi chu kỳ cao đại diện cho xu hướng tiên tiến nhất hiện nay. Điển hình là nghiên cứu của Horn (2014) thực hiện khảo sát trên 10 mẫu khung SMRM khác nhau, kết hợp đo đạc hơn 40 thông số động lực học (gia tốc thẳng đứng, áp suất bầu khí, ứng suất động) để hiệu chỉnh mô hình số. Abhijeet S. Hange (2018) hợp tác cùng Altair và TATA ứng dụng nền tảng CAE tích hợp nhằm trích xuất lịch sử tải trọng động để phân tích mỏi theo nguyên lý tích lũy tổn thương, giúp giảm 85% thời gian tính toán so với thử nghiệm truyền thống. Nghiên cứu của Abhishek Sharma và cộng sự (2020) tối ưu hóa dầm chassis chữ C bằng thép hợp kim AISI 4130, đạt mức giảm khối lượng 13% trong khi vẫn bảo toàn tần số dao động xoắn riêng.

Tại Việt Nam, các nghiên cứu của Nguyễn Phúc Hiểu (1996), Ngô Thế Quân (2007) và Phạm Lê Tiến (2011) đã đặt nền móng cho việc khảo sát dao động ngẫu nhiên và lan truyền vết nứt mỏi trên phương tiện giao thông. Giai đoạn gần đây, Vũ Tuấn Đạt (2015, 2017) và Trần Phúc Hòa (2016) đã ứng dụng ANSYS và nCode DesignLife để dự báo tuổi thọ mỏi cho dầm cầu và khung xe tải tự đổ. Luận án của Lương Văn Vạn (2020) sử dụng mô hình động lực học trên Matlab kết nối với FEM trên ANSYS nhưng vẫn dừng lại ở việc ngàm cứng các gối treo khi phân tích ứng suất. Các nghiên cứu về SMRM của Nguyễn Thanh Tùng (2016), Nguyễn Tiến Dũng (2017) và Tạ Tuấn Hưng (2019) chủ yếu tập trung vào động lực học phanh, ổn định hướng và chống lật ngang của đoàn xe, chưa đi sâu vào độ bền mỏi của kết cấu khung mang tải.

       TIẾN TRÌNH PHÁT TRIỂN CÁC PHƯƠNG PHÁP NGHIÊN CỨU KHUNG CHASSIS
┌──────────────────────────────────────────────────────────────────────────┐
│ Cơ học cổ điển (Timoshenko, Wong)                                        │
│ - Dầm giải tích tĩnh, hệ số tải Kd kinh nghiệm                           │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│ Phương pháp Phần tử hữu hạn độc lập (FEM) (Caban, Moaaz, Mi Chengji)     │
│ - Khảo sát tập trung ứng suất, nhưng biên ngàm tĩnh làm sai lệch động học│
└────────────────────────────────────┬─────────────────────────────────────┘
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│ Động lực học hệ nhiều vật rời rạc (MBD) (Nguyễn Tiến Dũng, Tạ Tuấn Hưng) │
│ - Đánh giá ổn định, lật ngang, phanh; xem khung là vật rắn tuyệt đối    │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│ VỊ TRÍ LUẬN ÁN: Đồng mô phỏng MBD - Flexible FEM - Thực nghiệm đa kênh  │
│ - Khung đàn hồi (CMS), tải động ngẫu nhiên ISO 8608, Rainflow/Goodman    │
└──────────────────────────────────────────────────────────────────────────┘

Luận án này định vị chính xác tại điểm giao thoa giữa động lực học hệ nhiều vật phi tuyến và cơ học phá hủy mỏi của kết cấu đàn hồi. Điểm đột phá so với các công trình quốc tế và trong nước là việc thiết lập mô hình SMRM đàn hồi hoàn chỉnh không bị ràng buộc cưỡng bức, kết hợp kiểm chứng thực nghiệm đồng bộ giữa biến dạng cơ học và thông số động học trên xe thật tại điều kiện đường sá Việt Nam.

Đó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à hoàn thiện khung lý thuyết động lực học kết cấu phương tiện vận tải hạng nặng thông qua các đóng góp cụ thể:

  1. Mở rộng lý thuyết Động lực học hệ nhiều vật (MBD) cho đoàn xe tổ hợp 5 trục: Thiết lập hệ phương trình vi phân đại số (Differential-Algebraic Equations - DAEs) mô tả chuyển động tương đối không gian 3 chiều giữa xe đầu kéo và sơ mi rơ moóc qua khớp cầu mâm xoay (chốt kéo Kingpin) kết hợp các chuỗi góc xoay Euler ($\phi, \theta, \psi$), ma trận chuyển đổi tọa độ từ hệ quy chiếu gắn liền vật (BCS) sang hệ tọa độ cố định mặt đất (GCS).
  2. Ứng dụng lý thuyết Tổ hợp dao động riêng thành phần (Component Mode Synthesis - CMS): Sử dụng phương pháp Craig-Bampton để giảm bậc tự do (modal reduction) cho khung SMRM đàn hồi. Cấu trúc ma trận khối lượng và độ cứng tổng quát được phân rã thành các bậc tự do biên ($u_B$) và bậc tự do nội ($u_I$), bảo toàn chính xác các dạng dao động riêng (mode shapes) tần số thấp mang năng lượng kích động chính từ mặt đường:

$$\mathbf{M}{gen} = \begin{bmatrix} \mathbf{M}{tt} & \mathbf{M}{tr} & \mathbf{M}{tm} \ \mathbf{M}{tr}^T & \mathbf{M}{rr} & \mathbf{M}{rm} \ \mathbf{M}{tm}^T & \mathbf{M}{rm}^T & \mathbf{M}{mm} \end{bmatrix}, \quad \mathbf{K}{gen} = \begin{bmatrix} \mathbf{0} & \mathbf{0} & \mathbf{0} \ \mathbf{0} & \mathbf{0} & \mathbf{0} \ \mathbf{0} & \mathbf{0} & \mathbf{\hat{K}}{NN} \end{bmatrix}$$

  1. Chuyển đổi hệ hình mẫu (Paradigm Shift) trong cơ học phá hủy khung chassis: Chuyển từ quan điểm "đánh giá độ bền tĩnh tương đương có điều kiện biên giả định" sang quan điểm "đánh giá tổn thương mỏi tích lũy động lực học phi tuyến liên tục". Kết quả nghiên cứu chứng minh rằng biên ngàm cứng truyền thống làm triệt tiêu mô men xoắn tự do của dầm chính khi xe chuyển động qua mấp mô lệch pha, dẫn đến ước tính sai lệch ứng suất von Mises cục bộ lên đến 32.6%.

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

Khung phân tích của luận án là sự tích hợp hữu cơ của 3 lý thuyết nền tảng: Lý thuyết Động lực học nhiều vật phi tuyến, Lý thuyết Phần tử hữu hạn đàn hồi phi tuyến hình học, và Lý thuyết Mỏi đa trục chu kỳ cao.

                  KHUNG PHÂN TÍCH ĐỒNG MÔ PHỎNG TÍCH HỢP
┌──────────────────────────────────────────────────────────────────────────┐
│                   MÔ HÌNH HÓA KÍCH ĐỘNG MẶT ĐƯỜNG                        │
│   - Mấp mô hình sin (Obstacle bump: H = 65mm, 200mm, 300mm)              │
│   - Mấp mô ngẫu nhiên PSD ISO 8608 (Cấp C-D, D-E; V = 20 - 100 km/h)    │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│       MÔ HÌNH ĐỘNG LỰC HỌC ĐOÀN XE (MBD) TRONG MOTIONSOLVE               │
│   - Xe đầu kéo 3 trục + SMRM 2 trục (Treo cân bằng nhíp lá liên động)    │
│   - Khớp nối mâm xoay Kingpin, ma trận biến đổi Euler 3D                 │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     │ (Tương tác lực & chuyển vị thời gian thực)
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│       MÔ HÌNH PHẦN TỬ HỮU HẠN ĐÀN HỒI (FEM) TRONG OPTISTRUCT             │
│   - Khung SMRM 40 feet (Craig-Bampton Component Mode Synthesis)          │
│   - Lưới CQUAD4/CTRIA3 + RBE2 (Kingpin) + RBE3 (Tải container phân bố)   │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│       PHÂN TÍCH MỎI TÍCH LŨY TRONG HYPERLIFE                             │
│   - Thuật toán đếm chu trình Rainflow (Rainflow Cycle Counting)          │
│   - Hiệu chỉnh ứng suất trung bình Gerber / Goodman trên biểu đồ Haigh   │
│   - Đường cong mỏi S-N Basquin của vật liệu thép Q345B                   │
└────────────────────────────────────┬─────────────────────────────────────┘
                                     ▼
┌──────────────────────────────────────────────────────────────────────────┐
│       KIỂM CHỨNG THỰC NGHIỆM ĐA KÊNH & TỐI ƯU HÓA KẾT CẤU                │
│   - Dewesoft DAQ + Cảm biến lá Strain Gauge + Quang học HF-750C + S-Motion│
│   - Đề xuất 4 phương án cải tiến: Khoét lỗ dầm ngang, tối ưu dầm chính   │
└──────────────────────────────────────────────────────────────────────────┘

Điều kiện biên và ràng buộc (Boundary Conditions): Khung xe được giải phóng hoàn toàn khỏi các ngàm cố định nhân tạo. Các điểm tiếp nhận lực bao gồm: mâm xoay (ràng buộc qua chùm phần tử cứng RBE2 truyền lực kéo/phanh và phản lực thẳng đứng $F_{z2}$), các gối đỡ quang treo và cân bằng của hệ thống treo nhíp lá trục 4 và trục 5 (ràng buộc động học với ma trận độ cứng $C_{4}, C_{5}$ và hệ số cản $K_{4}, K_{5}$), cùng tải trọng hàng hóa giả lập container phân bố đều thông qua mạng phần tử nội suy lực RBE3.

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ủ thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp luận phức hợp đa tầng (Triangulation Method):

  • Tầng mô hình hóa số trị (Numerical Simulation): Xây dựng mô hình 3D chính xác từng chi tiết của khung SMRM 40 feet 2 trục trong môi trường Altair HyperMesh, thiết lập động học đoàn xe trong MotionView và giải bài toán động lực học phi tuyến trong MotionSolve.
  • Tầng tính toán bền mỏi (Durability & Fatigue Analysis): Sử dụng HyperLife để xử lý các chuỗi ứng suất biến thiên theo thời gian trích xuất từ mô hình động lực học.
  • Tầng thực nghiệm hiện trường (Field Experimental Testing): Tiến hành đo đạc trực tiếp trên đoàn xe thực tế tại bãi thử tiêu chuẩn nhằm đối soát, hiệu chỉnh và xác thực mô hình số.

Đối tượng nghiên cứu thực nghiệm là tổ hợp đoàn xe gồm: Xe đầu kéo 3 trục (công thức bánh xe $6 \times 4$) kết nối với Sơ mi rơ moóc chở container 40 feet 2 trục do Việt Nam chế tạo. Vật liệu chế tạo dầm khung chính là thép kết cấu hợp kim thấp cường độ cao Q345B (giới hạn chảy $\sigma_{ch} = 345\text{ MPa}$, giới hạn bền kéo $S_u = 510\text{ MPa}$, mô đun đàn hồi $E = 2.1 \times 10^5\text{ MPa}$, hệ số Poisson $\nu = 0.3$).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được chuẩn hóa qua các bước thực nghiệm và tính toán nghiêm ngặt:

  1. Mô hình hóa phần tử hữu hạn khung SMRM:

    • Kết cấu dầm dọc chữ I biến tiết diện, dầm ngang, thanh gia cường và chốt kéo được chia lưới với kích thước phần tử trung bình 10 mm. Sử dụng chủ đạo phần tử vỏ tấm 2D bậc một 4 nút CQUAD4 và 3 nút CTRIA3 nằm trên mặt trung hòa (Mid-surface), kết hợp phần tử thanh dầm CROD mô tả dầm giằng ngang.
    • Kiểm soát chất lượng lưới gắt gao: Tỷ lệ co dãn (Aspect ratio) $< 5$, góc xiên (Skewness) $< 45^\circ$, độ cong vênh (Warpage) $< 15^\circ$, đảm bảo 100% phần tử vượt qua tiêu chuẩn kiểm định phần tử hữu hạn khắt khe của ngành ô tô.
    • Mô phỏng đường hàn liên kết giữa bản cánh (flange) và bản bụng (web) của dầm chữ I bằng dải phần tử CQUAD4 tương đương với đặc tính cơ lý mối hàn đạt chuẩn JIS D34.

  2. Thiết lập hệ thống đo thực nghiệm đa kênh:

    • Đo ứng suất/biến dạng: Sử dụng các cảm biến lá điện trở (Strain gauges) gồm loại đơn phương (1 chiều) và tam giác 3 phương (Triaxial Strain Gauge Rosettes) dán tại 12 vị trí xung yếu trên dầm chính, gần mâm xoay và khu vực quang treo trục sau. Cảm biến được kết nối theo sơ đồ mạch cầu Wheatstone (Quarter-bridge và Half-bridge) có bù nhiệt độ.
    • Đo chuyển vị động: Sử dụng cảm biến quang học phi tiếp xúc HF-750C với dải đo $\pm 750\text{ mm}$, độ phân giải cao để ghi nhận chuyển vị tương đối của khối lượng được treo tại cầu 4 và cầu 5.
    • Đo vận tốc và gia tốc: Sử dụng cảm biến quang học đo vận tốc dài S-Motion kết hợp cảm biến gia tốc áp điện 3 trục bố trí tại trọng tâm đoàn xe và trên các trục cầu.
    • Hệ thống thu thập dữ liệu (DAQ): Bộ xử lý tín hiệu DEWEsoft SIRIUS đồng bộ hóa tín hiệu với tần số lấy mẫu $f_s = 1000\text{ Hz}$, lọc nhiễu thông thấp Butterworth trước khi phân tích phổ.
                    SƠ ĐỒ HỆ THỐNG THỰC NGHIỆM ĐA KÊNH
┌──────────────────────────────────────────────────────────────────────────┐
│                           ĐOÀN XE SMRM THỰC TẾ                           │
│  [Xe đầu kéo 6x4] + [Mâm xoay / Chốt kéo] + [Khung SMRM 40ft Thép Q345B] │
└──────┬──────────────────────┬──────────────────────┬─────────────────────┘
       │                      │                      │
       ▼                      ▼                      ▼
┌──────────────┐      ┌──────────────┐      ┌──────────────────────────────┐
│ CẢM BIẾN LÁ  │      │ CẢM BIẾN     │      │ CẢM BIẾN QUANG & GIA TỐC     │
│ STRAIN GAUGE │      │ CHUYỂN VỊ    │      │ S-MOTION & ACCELEROMETER     │
│ (12 vị trí)  │      │ OPTICAL      │      │ (Đo vận tốc dài, gia tốc z)  │
│ Mạch cầu     │      │ HF-750C      │      │                              │
│ Wheatstone   │      │ (Cầu 4 & 5)  │      │                              │
└──────┬───────┘      └──────┬───────┘      └──────────────┬───────────────┘
       │                      │                             │
       └──────────────────────┼─────────────────────────────┘
                              ▼
┌──────────────────────────────────────────────────────────────────────────┐
│              BỘ THU THẬP & XỬ LÝ TÍN HIỆU ĐỒNG BỘ DEWESOFT               │
│               (Tần số lấy mẫu fs = 1000 Hz, Lọc Butterworth)             │
└─────────────────────────────┬────────────────────────────────────────────┘
                              ▼
┌──────────────────────────────────────────────────────────────────────────┐
│           SO SÁNH & XÁC THỰC MÔ HÌNH SỐ (MBD - FEM HYPERWORKS)           │
│         Hệ số tương quan Pearson rXY > 0.90; Sai số cực trị < 12%        │
└──────────────────────────────────────────────────────────────────────────┘
  1. Giao thức thử nghiệm hiện trường:
    • Thí nghiệm 1 (Tải tĩnh): Đặt tải cọc bê tông khối lượng chuẩn phân bố đều, xe đỗ trên mặt đường phẳng; tiếp đó cho bánh sau bên phải leo lên khối mấp mô chuẩn $H = 65\text{ mm}$ để đo trạng thái xoắn tĩnh.
    • Thí nghiệm 2 (Động lực học qua mấp mô): Đoàn xe chạy qua dải mấp mô hình sin (chạy đồng pha 2 vệt bánh và lệch pha 1 vệt bánh) ở các cấp tốc độ 20 km/h, 30 km/h và 40 km/h.

Data và phân tích

Dữ liệu thực nghiệm và mô phỏng được xử lý qua chuỗi thuật toán số học chuẩn hóa:

  • Đánh giá độ tin cậy mô hình: Hệ số tương quan tuyến tính Pearson ($r_{XY}$) giữa tín hiệu ứng suất mô phỏng và đo thực tế tại các điểm đo then chốt đều đạt giá trị từ $0.902$ đến $0.948$. Sai số cực trị chuyển vị tương đối của khối lượng được treo tại cầu 4 và cầu 5 dao động trong khoảng $6.8% \div 11.4%$, khẳng định mô hình MBD-FEM có độ tin cậy vượt trội.
  • Xây dựng đường cong mỏi S-N: Ước lượng đường cong S-N của thép Q345B theo phương trình Basquin hiệu chỉnh:

$$S_a = S_f' (2N_f)^b$$

Với hệ số độ bền mỏi $S_f' = 878.5\text{ MPa}$, số mũ chịu mỏi $b = -0.095$, giới hạn mỏi $S_e = 212\text{ MPa}$ tại $2 \times 10^6$ chu kỳ.

  • Thuật toán Rainflow & Biểu đồ Haigh: Quá trình trích xuất ứng suất động trên miền thời gian $\sigma(t)$ được phân rã thành ma trận biên độ ứng suất ($S_a$) và ứng suất trung bình ($S_m$). Ảnh hưởng của ứng suất trung bình kéo được hiệu chỉnh theo tiêu chuẩn Goodman và Gerber:

$$\frac{S_a}{S_e} + \left(\frac{S_m}{S_u}\right) = 1 \quad \text{(Goodman)}, \qquad \frac{S_a}{S_e} + \left(\frac{S_m}{S_u}\right)^2 = 1 \quad \text{(Gerber)}$$

Tổn thương mỏi lũy tích ($D$) được tính theo quy tắc Palmgren-Miner:

$$D = \sum_{i=1}^{k} \frac{n_i}{N_i}$$

Khi $D \ge 1.0$, kết cấu được coi là bị phá hủy mỏi hoàn toàn.

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

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

Kết quả mô phỏng và thực nghiệm mang lại 4 phát hiện khoa học mang tính bước ngoặt:

  1. Sự khuếch đại tải trọng động phi tuyến theo cấp đường và vận tốc: Khi đoàn xe di chuyển trên đường ngẫu nhiên tiêu chuẩn ISO 8608 loại C-D với vận tốc tăng từ 60 km/h lên 100 km/h, phản lực thẳng đứng động tại chốt kéo ($F_{z2}$) tăng 138%, trong khi tại trục cân bằng cầu 4-5 ($F_{z45}$) tăng 162%. Đặc biệt trên đường loại D-E ở vận tốc 40 km/h, biên độ tải trọng động đạt cực đại, tạo ra các xung lực va đập thẳng đứng gấp $2.14$ lần tải trọng tĩnh định mức.

  2. Quy luật phân bố ứng suất và nhận diện điểm nút nguy hiểm nhất: Ứng suất von Mises lớn nhất tập trung tại phần tử số #261150, nằm tại vị trí chuyển tiếp tiết diện (cổ ngỗng dầm dọc) và liên kết gối đỡ nhíp cầu trước SMRM. Dưới điều kiện vận hành trên đường D-E ở tốc độ 40 km/h, biên độ ứng suất động đạt $287.4\text{ MPa}$ (tiệm cận giới hạn chảy $345\text{ MPa}$). Ngược lại, tại trạng thái tĩnh thuần túy, ứng suất tại điểm này chỉ đạt $112.3\text{ MPa}$. Điều này chứng minh phương pháp tính tĩnh truyền thống đã bỏ sót hoàn toàn nguy cơ dẻo hóa cục bộ do tải động.

       SO SÁNH ỨNG SUẤT TẠI NÚT NGUY HIỂM #261150 THEO CÁC ĐIỀU KIỆN
┌──────────────────────────────────────────────────────────────────────────┐
│ Tĩnh định mức (TH01): 112.3 MPa                                          │
│ ██████████                                                               │
├──────────────────────────────────────────────────────────────────────────┤
│ Tĩnh xoắn mấp mô 30cm (TH03): 184.6 MPa                                  │
│ ████████████████                                                         │
├──────────────────────────────────────────────────────────────────────────┤
│ Động: Đường C-D @ 100 km/h: 224.8 MPa                                    │
│ ████████████████████                                                     │
├──────────────────────────────────────────────────────────────────────────┤
│ Động: Đường D-E @ 40 km/h: 287.4 MPa (Gần ngưỡng chảy 345 MPa!)          │
│ ██───────────────────────────────────────────────███████████████████████ │
└──────────────────────────────────────────────────────────────────────────┘

  1. Bản đồ tuổi thọ mỏi và sự suy giảm chu kỳ phá hủy: Trên mặt đường loại C-D ở tốc độ 60 km/h, vùng làm việc của khung có số chu kỳ mỏi đạt trên $10^7$ chu kỳ (vùng bền vô hạn). Tuy nhiên, khi vận tốc đạt 100 km/h hoặc khi xe chạy trên đường D-E (40 km/h), phổ chu kỳ mỏi tại vùng liên kết dầm chính và thanh giằng cầu sau tụt dốc nghiêm trọng xuống còn $1.26 \times 10^5 \div 4.85 \times 10^5$ chu kỳ, hoàn toàn khớp với thực tế xuất hiện các vết nứt chân chim tại mối hàn sau 1-2 năm vận hành khắc nghiệt.

  2. Hiệu quả vượt trội của các phương án cải tiến kết cấu (PACT): Luận án đã khảo sát 4 phương án thiết kế cải tiến kết cấu:

    • PACT 1 (Thay đổi mác thép dầm chính sang thép cường độ cao HSS): Giảm chiều dày bản cánh dầm chính từ 14 mm xuống 12 mm, ứng suất cực đại giảm 8.4%, tuổi thọ mỏi tăng 2.3 lần.
    • PACT 2 (Tối ưu hình học tiết diện dầm I): Tăng chiều cao thân dầm tại đoạn giữa thêm 50 mm, giảm chiều dày bản bụng, giúp mô men chống uốn $W_u$ tăng 14.2%.
    • PACT 3 (Khoét 3 lỗ giảm tải hình ô van trên dầm ngang): Giảm được 185 kg khối lượng kết cấu dầm ngang mà không làm tăng ứng suất uốn cục bộ quá 3.5%.
    • PACT 4 (Tích hợp tối ưu toàn diện): Giảm tự trọng bản thân khung xe $12.8%$ (tương đương giảm hơn 420 kg), phân bố lại dòng ứng suất đều trên toàn bộ khung gầm, nâng hệ số an toàn bền mỏi ($SF_{fatigue}$) từ $1.18$ lên $1.65$.
     MA TRẬN SO SÁNH CÁC PHƯƠNG ÁN THIẾT KẾ CẢI TIẾN KHUNG SMRM
┌───────────┬────────────────────┬───────────────────┬─────────────────────┐
│  Phương án│ Tự trọng khung (kg)│ Ứng suất max (MPa)│ Chu kỳ mỏi cực tiểu │
├───────────┼────────────────────┼───────────────────┼─────────────────────┤
│ Nguyên bản│ 3,280 (Chuẩn)      │ 287.4             │ 1.26 x 10^5 chu kỳ  │
│ PACT 1    │ 3,090 (-5.8%)      │ 263.2             │ 2.90 x 10^5 chu kỳ  │
│ PACT 2    │ 3,120 (-4.9%)      │ 251.0             │ 3.45 x 10^5 chu kỳ  │
│ PACT 3    │ 3,095 (-5.6%)      │ 278.1             │ 1.85 x 10^5 chu kỳ  │
│ PACT 4    │ 2,860 (-12.8%)     │ 242.6             │ 5.12 x 10^5 chu kỳ  │
└───────────┴────────────────────┴───────────────────┴─────────────────────┘

Implications đa chiều

  • Về mặt học thuật: Thiết lập khung phương pháp luận chuẩn xác cho bài toán tương tác động lực học giữa kết cấu đàn hồi và cơ cấu truyền động nhiều vật, đóng vai trò tài liệu tham khảo cốt lõi cho các nghiên cứu chuyên sâu về xe chuyên dụng và cơ khí động lực.
  • Về mặt công nghệ và công nghiệp: Chuyển giao trực tiếp quy trình thiết kế số hóa bằng bộ công cụ HyperWorks cho các tập đoàn sản xuất ô tô trong nước (như THACO, SAMCO, Tân Thanh). Doanh nghiệp có thể loại bỏ giai đoạn chế tạo thử nghiệm tốn kém, rút ngắn chu kỳ R&D từ 6-9 tháng xuống còn vài tuần với độ tin cậy được chứng thực.
  • Về mặt kinh tế - xã hội: Việc giảm 420 kg tự trọng cho mỗi khung SMRM cho phép tăng tải trọng hàng hóa hữu ích tương ứng, tiết kiệm hàng triệu lít nhiên liệu tiêu thụ và cắt giảm lượng phát thải carbon đáng kể trên quy mô toàn bộ đội xe vận tải quốc gia.
  • Về mặt chính sách: Cung cấp cơ sở khoa học xác đáng cho Cục Đăng kiểm Việt Nam và Bộ Giao thông Vận tải trong việc bổ sung, cập nhật Quy chuẩn kỹ thuật quốc gia QCVN 11:2011/BGTVT, chuyển từ kiểm tra bền tĩnh đơn thuần sang yêu cầu đánh giá độ bền mỏi và mô phỏng động lực học bắt buộc đối với phương tiện siêu trường siêu trọng.

Limitations và Future Research

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

  1. Giới hạn về trường tải trọng: Nghiên cứu chủ yếu tập trung vào kích động tải trọng thẳng đứng và tải trọng xoắn do mấp mô mặt đường gây ra. Chưa xét đến sự kết hợp đồng thời của tải trọng phanh khẩn cấp cường độ cao (quán tính dọc cực hạn) và tải trọng lực ly tâm khi xe quay vòng gấp trên đường đèo dốc (mất ổn định lật ngang).
  2. Giới hạn về mô hình hóa mối hàn: Mối hàn được mô phỏng thông qua dải phần tử vỏ tương đương CQUAD4 trên mặt trung hòa mà chưa đi sâu vào phân tích ứng suất tập trung tế vi tại chân mối hàn theo phương pháp dải ứng suất cấu trúc (Structural Hot-Spot Stress) hoặc cơ học nứt tế vi đa trục.
  3. Môi trường vận hành: Nghiên cứu chưa tích hợp ảnh hưởng của môi trường ăn mòn nhiệt đới ẩm, muối biển và dao động nhiệt độ làm suy giảm giới hạn mỏi theo thời gian (Corrosion-Fatigue interaction).

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: Thu thập dữ liệu tải trọng đường thực tế (Road Load Data Acquisition - RLDA) thông qua mạng lưới cảm biến gắn trên xe chạy liên tục trên các tuyến quốc lộ huyết mạch Bắc - Nam.
  • Hướng 2: Ứng dụng thuật toán tối ưu hóa hình thái đa mục tiêu (Multi-objective Topology and Topometry Optimization) kết hợp trí tuệ nhân tạo (Machine Learning Surrogate Models) để tự động hóa hoàn toàn việc thiết kế khung chassis thế hệ mới.
  • Hướng 3: Nghiên cứu độ bền mỏi đa trục của các vật liệu tiên tiến siêu nhẹ như thép hợp kim boron xử lý nhiệt, nhôm hợp kim chuyên dụng hoặc vật liệu composite composite nền sợi carbon/kim loại cho xe chuyên dụng.
  • Hướng 4: Xây dựng mô hình bản sao số thời gian thực (Digital Twin) kết hợp cảm biến thông minh IoT tích hợp trên khung xe để theo dõi sức khỏe kết cấu (Structural Health Monitoring - SHM) và cảnh báo nứt mỏi sớm trong quá trình khai thác.

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

  • Tác động học thuật: Công trình tạo ra bước đột phá phương pháp luận, mở ra hướng nghiên cứu liên ngành giữa Động lực học hệ nhiều vật (MBD), Cơ học tính toán (CAE) và Cơ học phá hủy mỏi tại Việt Nam. Các bài báo công bố từ luận án trên các tạp chí khoa học uy tín trong nước và quốc tế tạo tiền đề gia tăng chỉ số trích dẫn trong lĩnh vực kỹ thuật cơ khí động lực.
  • Chuyển đổi công nghiệp: Thúc đẩy các doanh nghiệp cơ khí ô tô nội địa chuyển dịch từ phương thức gia công theo mẫu sang năng lực tự chủ thiết kế và tối ưu hóa sản phẩm công nghệ cao. Tăng cường tỉ lệ nội địa hóa thực chất của sơ mi rơ moóc lên trên 60%, nâng cao năng lực cạnh tranh sòng phẳng với các sản phẩm nhập khẩu từ Trung Quốc và ASEAN.
  • Lợi ích xã hội và môi trường: Nâng cao hệ số an toàn vận tải, giảm thiểu rủi ro tai nạn giao thông thảm khốc do gãy đứt khung chassis khi chở hàng quá tải trên đường đèo dốc. Tối ưu hóa trọng lượng giúp giảm tiêu hao nhiên liệu khoảng $0.8 \div 1.2\text{ lít}/100\text{ km}$ cho mỗi đoàn xe, đóng góp trực tiếp vào mục tiêu Net-Zero phát thải của quốc gia.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ khí Động lực: Tiếp cận một tài liệu tham khảo học thuật chuẩn mực về quy trình đồng mô phỏng phức hợp MBD-FEM và phương pháp thực nghiệm đo đạc động lực học phương tiện.
  • Kỹ sư NC&PT (R&D) tại các nhà máy ô tô: Sở hữu một cẩm nang hướng dẫn thực hành từng bước trên nền tảng HyperWorks (HyperMesh, MotionView, MotionSolve, OptiStruct, HyperLife) để ứng dụng trực tiếp vào việc phát triển các dòng xe mới.
  • Doanh nghiệp vận tải và Logistics: Được sử dụng các sản phẩm SMRM sản xuất trong nước có tự trọng nhẹ hơn, tải trọng hàng hóa cao hơn, độ bền bỉ vượt trội và thời gian bảo hành kéo dài.
  • Cơ quan quản lý Nhà nước (Cục Đăng kiểm, Bộ GTVT): Có luận cứ khoa học vững chắc để hoàn thiện hệ thống văn bản quy chuẩn, tiêu chuẩn kỹ thuật quốc gia về kiểm định an toàn phương tiện cơ giới đường bộ.

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?

Luận án đã mở rộng lý thuyết Động lực học hệ nhiều vật đàn hồi phi tuyến thông qua việc kết hợp thuật toán khử bậc tự do Craig-Bampton (CMS) với phương trình vi phân đại số Lagrange (DAEs) cho tổ hợp đoàn xe 5 trục. Đóng góp độc đáo nhất là việc loại bỏ hoàn toàn các ràng buộc ngàm cứng cưỡng bức phi thực tế trong các mô hình FEM truyền thống, thiết lập mô hình khung chassis ở trạng thái "tự do đàn hồi hoàn toàn" dưới tác động của dòng tải trọng động biến thiên liên tục theo thời gian trích xuất từ mô hình MBD.

2. Điểm đổi mới về mặt phương pháp luận so với ít nhất 2 nghiên cứu quốc tế là gì?

So với nghiên cứu của Mi Chengji và cộng sự (2016) vốn chỉ sử dụng phân tích ứng suất quasi-static tĩnh rời rạc trên ANSYS để ước lượng mỏi gối treo, luận án đã thực hiện đồng mô phỏng động lực học thời gian thực (Full Dynamic Transient Co-simulation) phản ánh chính xác tương tác chuyển động không gian 3 chiều. So với nghiên cứu của Horn (2014) khảo sát thực nghiệm thuần túy trên các xe mẫu, luận án đã xây dựng quy trình số hóa khép kín từ mô phỏng MBD-FEM đến thuật toán đếm chu trình Rainflow/Goodman, đồng thời tối ưu hóa giảm $12.8%$ khối lượng khung mà nghiên cứu của Horn chưa thực hiện.

3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?

Điểm bất ngờ nhất là ứng suất cục bộ tại vị trí nút #261150 (cổ ngỗng dầm dọc) trong điều kiện xe chạy động qua đường gồ ghề D-E ở tốc độ 40 km/h đạt tới $287.4\text{ MPa}$, cao gấp $2.56$ lần so với giá trị tính toán trong điều kiện chịu tải tĩnh định mức ($112.3\text{ MPa}$). Điều này chỉ ra rằng các tính toán an toàn tĩnh truyền thống với hệ số $K_d = 1.5$ vẫn đánh giá thấp tải trọng động thực tế tới $41%$, giải thích triệt để nguyên nhân nhiều khung xe bị nứt gãy sớm dù đạt chuẩn kiểm định tĩnh.

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

Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ tham số đầu vào: bảng thông số hình học dầm I, ma trận độ cứng lò xo nhíp ($C_1, C_{23}, C_4, C_5$), hệ số cản giảm chấn ($K$), đặc tính lốp kép ($C_L, K_L$), tiêu chuẩn chất lượng lưới FEM (kích thước, loại phần tử CQUAD4/CTRIA3), thông số cơ lý và đường cong mỏi Basquin của thép Q345B, cùng sơ đồ kết nối mạch cầu Wheatstone và cấu hình tần số lấy mẫu trên phần mềm DEWEsoft. Bất kỳ phòng thí nghiệm cơ khí động lực nào cũng có thể tái lập chính xác quy trình mô phỏng và thực nghiệm này.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm tập trung vào 3 giai đoạn: (1) Giai đoạn 2024-2027: Chuẩn hóa bộ dữ liệu phổ tải trọng mặt đường Việt Nam (RLDA) và thương mại hóa gói phần mềm thiết kế tối ưu SMRM cho doanh nghiệp nội địa; (2) Giai đoạn 2027-2030: Ứng dụng vật liệu composite/kim loại nhẹ và công nghệ in 3D kim loại cho các khớp nối khung xe chuyên dụng; (3) Giai đoạn 2030-2035: Tích hợp hệ thống cảm biến quang sợi Bragg (FBG) vào khung xe để hình thành kết cấu khung thông minh tự chẩn đoán nứt mỏi theo thời gian thực (Self-monitoring Smart Chassis).

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán cấp thiết về đánh giá độ bền và dự báo tuổi thọ mỏi cho khung sơ mi rơ moóc sản xuất tại Việt Nam bằng phương pháp tiếp cận hiện đại, đồng bộ và khoa học.
  2. Thiết lập thành công mô hình đồng mô phỏng tích hợp giữa Động lực học hệ nhiều vật (MBD) và Phần tử hữu hạn đàn hồi (Flexible FEM), loại bỏ triệt để các sai số do điều kiện biên cưỡng bức của phương pháp tĩnh truyền thống.
  3. Nhận diện chính xác vị trí tập trung ứng suất nguy hiểm nhất (nút #261150) và cơ chế phá hủy mỏi chu kỳ cao dưới kích động của mấp mô mặt đường ngẫu nhiên theo tiêu chuẩn ISO 8608.
  4. Xây dựng hệ thống thực nghiệm đo động lực học và ứng suất đa kênh đồng bộ thời gian thực với độ tin cậy rất cao, hệ số tương quan Pearson $r_{XY} > 0.90$ so với mô hình số.
  5. Đề xuất thành công phương án cải tiến kết cấu tổng thể PACT 4, giúp cắt giảm $12.8%$ khối lượng tự trọng khung xe (giảm 420 kg) trong khi gia tăng tuổi thọ mỏi và hệ số an toàn kết cấu lên $1.65$.
  6. Đóng góp của luận án mở ra một kỷ nguyên mới trong năng lực làm chủ công nghệ thiết kế, tính toán và chế tạo phương tiện cơ giới hạng nặng, đóng góp thiết thực vào sự phát triển bền vững của ngành công nghiệp ô tô Việt Nam.