Giới thiệu dự án

Tai nạn giao thông đường bộ là một trong những thách thức toàn cầu nghiêm trọng nhất đối với sức khỏe cộng đồng và kinh tế xã hội. Theo thống kê của Tổ chức Y tế Thế giới (WHO), mỗi năm có khoảng 1,3 triệu người tử vong và hơn 50 triệu người bị thương tật vĩnh viễn do tai nạn giao thông. Tại Việt Nam, theo báo cáo của Cục Cảnh sát Giao thông trong Quý I/2024, toàn quốc xảy ra 6.552 vụ tai nạn giao thông, làm 2.723 người tử vong và 5.246 người bị thương (trung bình mỗi ngày xảy ra 72 vụ, cướp đi sinh mạng của 30 người). Đáng chú ý, tỷ lệ tử vong do tai nạn giao thông tại Việt Nam đứng thứ hai trong khu vực Đông Nam Á (khoảng 30 ca tử vong/100.000 dân vào năm 2020), chỉ xếp sau Thái Lan.

+--------------------------------------------------------------------------------+
|                         THỰC TRẠNG AN TOÀN GIAO THÔNG                          |
|  - Toàn cầu (WHO): 1,3 triệu tử vong/năm | 50 triệu thương tật vĩnh viễn/năm   |
|  - Việt Nam (Q1/2024): 6.552 vụ tai nạn | 2.723 người chết | 5.246 người bị thương|
|  - Thách thức R&D: Chi phí thử nghiệm va chạm thực tế (Destructive Crash Test)  |
|    rất lớn (150.000 - 300.000 USD/mẫu), thiếu dữ liệu động lực học phi tuyến.   |
+--------------------------------------------------------------------------------+

Bài toán đặt ra cho ngành công nghiệp sản xuất ô tô và các viện nghiên cứu là: Làm thế nào để dự đoán chính xác mức độ biến dạng kết cấu khung vỏ, sự suy giảm năng lượng động lực học và bảo vệ khoang hành khách (Cabin Safety Cell) mà không phải tiêu tốn hàng triệu USD cho các bài thử nghiệm va chạm phá hủy thực tế (Destructive Crash Testing)?

Đề tài tốt nghiệp chuyên ngành Công nghệ Kỹ thuật Ô tô với tiêu đề "Ứng dụng phần mềm LS-DYNA mô phỏng và đánh giá mức độ biến dạng của ô tô khi va chạm" do sinh viên Nguyễn Thành Danh và Trần Quốc Hiếu thực hiện dưới sự hướng dẫn của TS. Nguyễn Văn Trạng (Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh) đã giải quyết triệt để bài toán này thông qua phương pháp phần tử hữu hạn phi tuyến (Non-linear Finite Element Method).

Mục tiêu nghiên cứu cụ thể:

  1. Xây dựng và chuẩn hóa mô hình phần tử hữu hạn (FEM) 3D cho các phân khúc xe phổ biến: Bán tải (Pickup Truck - Chevrolet C2500), Sedan (Toyota Yaris 2010) và Crossover/SUV (Toyota RAV4).
  2. Thiết lập chính xác các điều kiện biên động lực học (Dynamic Boundary Conditions), đặc tính vật liệu đàn - dẻo phụ thuộc tốc độ biến dạng (Strain-rate Dependent Elasto-plasticity) và ma sát tiếp xúc (Contact Formulations).
  3. Khảo sát chi tiết quá trình hấp thu năng lượng, biến dạng kết cấu khung gầm ở dải vận tốc va chạm từ 56 km/h, 64 km/h đến 80 km/h theo các tiêu chuẩn quốc tế ASEAN NCAP, Euro NCAP và IIHS.
  4. Đánh giá tính tương thích va chạm (Crash Compatibility) giữa các dòng xe có sự chênh lệch khối lượng lớn và đề xuất giải pháp cải tiến kết cấu thanh dầm dọc phía trước (Front Crash Box/Rail) nhằm tối ưu hóa khả năng hấp thu xung lực.

Phạm vi và giới hạn đề tài:

  • Đối tượng khảo sát: Mô hình kết cấu khung vỏ xe Chevrolet C2500, Toyota Yaris 2010 và Toyota RAV4.
  • Các kịch bản va chạm: Va chạm trực diện 100% diện tích đầu xe vào tường cứng tuyệt đối (Full Frontal Rigid Barrier) ở 64 km/h và 80 km/h; Va chạm đối đầu trực diện giữa 2 xe bán tải ở 56 km/h; Va chạm vuông góc $90^\circ$ giữa 2 xe bán tải ở 56 km/h; Va chạm trực diện giữa xe bán tải (khối lượng lớn) với xe Sedan và SUV.
  • Giới hạn mô hình: Khảo sát tập trung vào biến dạng kết cấu cơ học khung vỏ kim loại; chưa tích hợp mô hình hình nhân sinh cơ học (Biomechanics Anthropomorphic Test Dummy) để đo chỉ số chấn thương chi tiết HIC (Head Injury Criterion).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trong quá trình phát triển phương tiện giao thông, việc đánh giá độ an toàn thụ động (Passive Safety) đòi hỏi sự kết hợp giữa các phương pháp vật lý và mô phỏng số.

Tiêu chí so sánh Thử nghiệm thực tế (Crash Test) Mô phỏng động học vật thể cứng (MBD) Mô phỏng phần tử hữu hạn phi tuyến (LS-DYNA)
Chi phí thực hiện Cực kỳ cao (150.000 - 300.000 USD/xe) Rất thấp (chỉ cần phần mềm) Trung bình (chi phí bản quyền & máy trạm)
Thời gian thiết lập Vài tuần đến vài tháng Vài giờ Vài ngày (Pre-processing & Meshing)
Độ chi tiết kết quả Giới hạn tại các điểm gắn cảm biến Chỉ phân tích chuyển vị tổng thể Chi tiết từng phần tử: Ứng suất, Biến dạng, Năng lượng
Khả năng lặp lại Rất khó lặp lại chính xác 100% Hoàn hảo Hoàn hảo (Deterministic Numerical Solver)
Khảo sát phi tuyến Trực tiếp nhưng phá hủy mẫu Kém (không tính biến dạng dẻo sâu) Xuất sắc (Phi tuyến hình học, vật liệu, tiếp xúc)

Phân loại yêu cầu hệ thống theo mô hình MoSCoW:

  • Must have (Bắt buộc phải có): Bộ giải tường minh động lực học phi tuyến (Explicit Dynamic Solver); Đường cong ứng suất - biến dạng thực (True Stress-Strain Curve); Thuật toán chống xuyên thấu phần tử (Penalty Contact Formulation); Định luật bảo toàn năng lượng với sai số Hourglass $< 5%$.
  • Should have (Nên có): Mô hình hóa chi tiết dầm hấp thu xung lực (Crash Box) dạng ống mỏng; Khảo sát biến thiên gia tốc giảm tốc (Deceleration Pulse).
  • Could have (Có thể mở rộng): Tối ưu hóa đa mục tiêu hình học thanh dầm chịu lực (Trigger beads/Grooves); Khảo sát va chạm đa phương (Oblique/Side impact).
  • Won't have (Chưa thực hiện đợt này): Mô hình hóa biến dạng xé rách kèm nổ túi khí và hệ thống căng đai khẩn cấp chủ động.

Thiết kế hệ thống

Quy trình mô phỏng va chạm ô tô được thực hiện theo cấu trúc khép kín từ tiền xử lý (Pre-processing), giải toán số (Explicit Solving) đến hậu xử lý (Post-processing).

flowchart TD
    A["Mô hình CAD 3D Hình học"] --> B["Tiền xử lý (LS-PrePost 4.8)"]
    B --> C["Chia lưới FEM: Belytschko-Tsay Shell Elements"]
    C --> D["Gán vật liệu: Mat_024 Piecewise Linear Plasticity"]
    D --> E["Thiết lập điều kiện biên & Vận tốc ban đầu (*INITIAL_VELOCITY)"]
    E --> F["Khai báo tiếp xúc (*CONTACT_AUTOMATIC_SINGLE_SURFACE)"]
    F --> G["Bộ giải LS-DYNA Explicit Solver (LS-Run)"]
    G --> H{"Kiểm tra cân bằng năng lượng"}
    H -- "Hourglass > 5% hoặc Lỗi tiếp xúc" --> B
    H -- "Bảo toàn năng lượng hợp lệ" --> I["Hậu xử lý (LS-PrePost 4.8)"]
    I --> J["Xuất biểu đồ: Vận tốc, Gia tốc, Stress Von Mises, Biến dạng dẻo"]

Công nghệ và công cụ sử dụng:

  • Bộ giải số (FEA Solver): LS-DYNA R11.1 / R12.0 (Ansys / LSTC) - Tối ưu hóa tính toán song song SMP/MPP.
  • Môi trường tiền & hậu xử lý: LS-PrePost v4.8 / v4.10 (Hỗ trợ định dạng đầu ra d3plot, binout, glstat, matsum).
  • Quản lý phiên tính toán: LS-Run 2024.
  • Nền tảng thiết kế cấu trúc: CATIA V5 / SolidWorks 2023.

Cơ sở toán học và cơ học va chạm:

  1. Phương trình chuyển động động lực học tổng quát: $$\mathbf{M}\ddot{\mathbf{u}}(t) + \mathbf{C}\dot{\mathbf{u}}(t) + \mathbf{K}\mathbf{u}(t) = \mathbf{F}^{ext}(t)$$ Trong đó $\mathbf{M}$ là ma trận khối lượng đường chéo (Lumped Mass Matrix), $\mathbf{C}$ là ma trận cản, $\mathbf{K}$ là ma trận độ cứng, $\mathbf{u}$ là vector chuyển vị và $\mathbf{F}^{ext}$ là ngoại lực tác dụng.

  2. Tích phân thời gian tường minh (Explicit Central Difference Integration): Bước thời gian tính toán giới hạn $(\Delta t)$ phải thỏa mãn điều kiện ổn định Courant-Friedrichs-Lewy (CFL): $$\Delta t \le \Delta t_{crit} = \frac{L_e}{c} = L_e \sqrt{\frac{\rho}{E}}$$ Với $L_e$ là chiều dài đặc trưng của phần tử nhỏ nhất, $c$ là vận tốc truyền sóng âm trong vật liệu, $\rho$ là khối lượng riêng và $E$ là mô đun đàn hồi Young.

  3. Mô hình vật liệu dẻo đàn hồi Johnson-Cook: $$\sigma_y = \left(A + B \varepsilon_p^n\right) \left(1 + C \ln \frac{\dot{\varepsilon}_p}{\dot{\varepsilon}_0}\right) \left[1 - \left(T^\right)^m\right]$$ Trong đó: $A$ là ứng suất chảy ban đầu, $B$ là hệ số biến cứng, $n$ là số mũ biến cứng, $C$ là hệ số nhạy tốc độ biến dạng, $\dot{\varepsilon}_p/\dot{\varepsilon}_0$ là tỷ số tốc độ biến dạng dẻo chuẩn hóa, và $T^$ là nhiệt độ tương đối.

  4. Định luật biến thiên năng lượng va chạm đối đầu giữa 2 xe: Năng lượng hấp thu biến dạng dẻo cực đại ($\Delta E$) của hệ 2 xe có khối lượng $m_1, m_2$ và vận tốc tương đối $v_r = v_{01} - v_{02}$ được tính theo công thức: $$\Delta E = \frac{m_1 m_2}{2(m_1 + m_2)} v_r^2$$

Phương pháp triển khai (Methodology)

Dự án áp dụng phương pháp nghiên cứu thực nghiệm số (Numerical Experimental Methodology) kết hợp đối chuẩn (Benchmarking) theo các mốc thời gian chặt chẽ:

+-----------------------------------------------------------------------------------+
|                           LỘ TRÌNH THỰC HIỆN DỰ ÁN                                |
|  [Tháng 1-2]: Thu thập dữ liệu hình học, chuẩn hóa lý thuyết biến dạng dẻo        |
|  [Tháng 3]: Xây dựng lưới FEM & gán thuộc tính vật liệu *MAT_024                  |
|  [Tháng 4]: Thiết lập bài toán biên va chạm đơn và đa xe trên LS-PrePost          |
|  [Tháng 5]: Chạy mô phỏng trên LS-Run, kiểm tra cân bằng năng lượng (glstat)      |
|  [Tháng 6]: Đánh giá biến dạng khoang lái, tối ưu hóa thiết kế dầm Crash Box      |
+-----------------------------------------------------------------------------------+

Ma trận đánh giá rủi ro kỹ thuật và giải pháp xử lý:

Rủi ro kỹ thuật (Technical Risk) Mức độ Biện pháp giảm thiểu & Giải pháp kỹ thuật
Lỗi phần tử âm (Negative Volume) Cao Sử dụng công thức tích phân giảm điểm Belytschko-Tsay; Kích hoạt tùy chọn tự hủy phần tử khi biến dạng vượt ngưỡng phá hủy (*MAT_ADD_EROSION).
Hiện tượng biến dạng ảo (Hourglass Mode) Cao Áp dụng kiểm soát năng lượng Hourglass Type 4/6 (Flanagan-Belytschko stiffness/viscous form) với hệ số $QM = 0.10$.
Xuyên thấu tiếp xúc (Contact Penetration) Trung bình Khai báo *CONTACT_AUTOMATIC_SINGLE_SURFACE với thuật toán Penalty Stiffness Scaling Factor $SLSFAC = 0.9$.
Bước thời gian quá nhỏ ($\Delta t < 10^{-7} s$) Trung bình Áp dụng kỹ thuật tăng khối lượng cục bộ có kiểm soát (Mass Scaling) qua thẻ *CONTROL_TIMESTEP, khống chế độ tăng khối lượng toàn hệ thống $< 1%$.

Thực thi và Kết quả nghiên cứu

Quá trình thực thi

Mô hình phần tử hữu hạn của các dòng xe được cấu thành từ hàng trăm nghìn phần tử vỏ (Shell Elements) dạng tứ giác và tam giác, được phân bố mật độ dày tại khu vực đầu xe (Crumple Zone) và thưa dần về phía đuôi xe để tối ưu hóa thời gian tính toán.

*KEYWORD
*TITLE
MOH PHONG VA CHAM TRUC DIEN O TO - LS-DYNA SIMULATION
*CONTROL_TERMINATION
$  ENDTIM    ENDCYC     DTMIN    ENDENG    ENDMAS
    0.120         0       0.0       0.0       0.0
*CONTROL_TIMESTEP
$  DTINIT    TSSFAC      ISDO    TSLIMT     DTMS
      0.0       0.9         0       0.0 -1.20E-06
*CONTROL_HOURGLASS
$     IHQ        QH
        4      0.10
*MAT_PIECEWISE_LINEAR_PLASTICITY
$   MID       RO          E        PR      SIGY      ETAN      FAIL      TDEL
      1 7.85E-09     210000      0.30     310.0    1200.0      0.35       0.0
*CONTACT_AUTOMATIC_SINGLE_SURFACE
$  SSID      MSID     SSTYP     MSTYP    SBOXID    MBOXID       SPR       MPR
      0         0         0         0         0         0         1         1
$    FS        FD        DC        VC       VDC    PENCHK    BT/DT       SFS
   0.30      0.20       0.0       0.0      20.0         0       0.0       1.0
*INITIAL_VELOCITY_GENERATION
$   SID      STYP      OMEGA        VX        VY        VZ     IVATN      ICID
      1         2          0   17777.8       0.0       0.0         0         0
*END

Đoạn mã trên mô tả thẻ dữ liệu Keyword điều khiển bước thời gian tính toán, khai báo vật liệu thép kết cấu dẻo, điều kiện tiếp xúc ma sát bề mặt và vận tốc ban đầu $17.778\text{ mm/s}$ (tương đương $64\text{ km/h}$).

Thử nghiệm và kiểm chứng mô hình

Tính đúng đắn của mô hình tính toán số được kiểm chứng thông qua định luật bảo toàn năng lượng trong tập tin đầu ra glstat.

==================================================================================
KIỂM CHỨNG CÂN BẰNG NĂNG LƯỢNG TOÀN HỆ THỐNG (ENERGY BALANCE VERIFICATION)
==================================================================================
Tổng năng lượng:       E_total = E_kinetic + E_internal + E_contact + E_hourglass
Độ lệch năng lượng:    |(E_total(t) - E_total(0)) / E_total(0)| <= 1.5% (ĐẠT CHUẨN)
Tỷ lệ năng lượng ảo:   E_hourglass / E_total <= 3.8% (TIÊU CHUẨN QUỐC TẾ: < 5.0%)
Khối lượng gia tăng:   Delta_Mass / Mass_total = 0.28% (TIÊU CHUẨN: < 1.0%)
==================================================================================

Phân tích chi tiết các kịch bản va chạm thực nghiệm:

+-----------------------------------------------------------------------------------+
|                        CÁC TRƯỜNG HỢP MÔ PHỎNG VA CHẠM                            |
| 1. Chevrolet C2500 vs Tường cứng tuyệt đối @ 64 km/h (Chuẩn Euro NCAP/IIHS)       |
| 2. Chevrolet C2500 vs Tường cứng tuyệt đối @ 80 km/h (Khảo sát quá tải tốc độ)    |
| 3. Xe Bán tải vs Xe Bán tải: Va chạm đối đầu 100% diện tích @ 56 km/h             |
| 4. Xe Bán tải vs Xe Bán tải: Va chạm góc vuông 90 độ @ 56 km/h                    |
| 5. Xe Bán tải (2.2 tấn) vs Xe Sedan Toyota Yaris (1.1 tấn) @ 56 km/h              |
| 6. Xe Bán tải (2.2 tấn) vs Xe SUV Toyota RAV4 (1.6 tấn) @ 56 km/h                 |
+-----------------------------------------------------------------------------------+

  1. Va chạm xe Bán tải vào tường cứng ở 64 km/h so với 80 km/h:

    • Ở vận tốc $64\text{ km/h}$, toàn bộ động năng ban đầu chuyển hóa thành nội năng biến dạng dẻo trong khoảng thời gian $\Delta t = 0.08\text{ s} - 0.10\text{ s}$. Phần đầu xe co rúm hoàn toàn tại khu vực dầm cản trước và dầm dọc, khoang hành khách (A-pillar và khung cửa) không bị xâm lấn nghiêm trọng.
    • Khi tăng vận tốc lên $80\text{ km/h}$, động năng tăng $56.25%$ ($E_k \propto v^2$). Biến dạng dẻo vượt qua giới hạn hấp thu của dầm trước, lan sâu vào vách ngăn động cơ (Firewall), trục lái bị đẩy lùi về phía sau gây nguy hiểm trực tiếp đến chi dưới và lồng ngực của người lái.

  2. Va chạm tương thích giữa các phân khúc xe (Bán tải vs Sedan & SUV):

    • Trong va chạm đối đầu giữa Chevrolet C2500 (khối lượng $\approx 2.2\text{ tấn}$) và Toyota Yaris (khối lượng $\approx 1.1\text{ tấn}$) ở vận tốc đối đầu $56\text{ km/h}$, xe Sedan chịu mức gia tốc giảm tốc (Deceleration Peak) lên tới hơn $38g$, trong khi xe Bán tải chỉ chịu khoảng $18g$.
    • Do sự chênh lệch độ cao khung gầm (Ride Height Mismatch), thanh dầm của xe Bán tải có xu hướng trèo lên đầu xe Sedan (Overriding Effect), ép bẹp nắp capo và đe dọa trực tiếp đến không gian sinh tồn bên trong khoang lái của chiếc Yaris.

Kết quả đạt được

Kịch bản va chạm Vận tốc Gia tốc đỉnh ($g$) Năng lượng hấp thu biến dạng ($kJ$) Độ lùi trụ A (A-Pillar Intrusion) Đánh giá an toàn theo IIHS
Pickup vs Tường cứng 64 km/h 24.5 420.5 45 mm Tốt (Good)
Pickup vs Tường cứng 80 km/h 36.8 655.8 168 mm Kém (Poor)
Pickup vs Pickup (Đối đầu) 56 km/h 22.1 512.4 52 mm Chấp nhận được (Acceptable)
Pickup vs Pickup (Góc $90^\circ$) 56 km/h 28.4 388.7 142 mm (Cửa xe) Biên giới hạn (Marginal)
Pickup vs Sedan Yaris 56 km/h Pickup: 18.2 / Yaris: 38.6 485.2 Yaris: 195 mm Kém cho xe nhỏ (Poor for Sedan)
Pickup vs SUV RAV4 56 km/h Pickup: 20.4 / SUV: 27.5 502.1 SUV: 78 mm Chấp nhận được (Acceptable)

Đổi mới và Đóng góp kỹ thuật

Cải tiến kết cấu dầm hấp thu xung lực (Front Longitudinal Rail)

Điểm sáng tạo cốt lõi của đề tài là việc nghiên cứu cải tiến hình học của dầm chịu lực đầu xe (Crash Box) - bộ phận đóng vai trò hấp thu phần lớn động năng va chạm thông qua cơ chế sụp đổ tiến triển (Progressive Folding Mode).

MẪU THIẾT KẾ DẦM CHỊU LỰC ĐẦU XE (CRASH BOX PROFILES):
--------------------------------------------------------------------------------
1. Mẫu 01: Dầm ống thép tiết diện vuông nguyên bản (Standard Square Tube)
2. Mẫu 02: Dầm tạo gân tăng cứng cục bộ (Stiffener Ribs)
3. Mẫu 03: Dầm đột lỗ giảm tải và định hướng gập (Trigger Holes)
4. Mẫu 04: Dầm dạng ống lượn sóng định hình biến dạng (Corrugated Profile)
--------------------------------------------------------------------------------

graph LR
    subgraph "Dầm nguyên bản (Mẫu 01)"
        A1["Lực đỉnh F_max rất lớn"] --> B1["Biến dạng uốn không ổn định (Global Buckling)"]
        B1 --> C1["Gia tốc truyền vào cabin cao"]
    end
    subgraph "Dầm cải tiến tạo rãnh định hướng (Mẫu 04)"
        A2["Lực kích hoạt F_initial được kiểm soát"] --> B2["Biến dạng nếp gấp đều đặn (Progressive Folding)"]
        B2 --> C2["Hấp thu năng lượng tăng +28.4%, giảm chấn thương"]
    end

Bảng so sánh hiệu suất hấp thu năng lượng giữa các kết cấu dầm:

Mẫu thử nghiệm kết cấu Khối lượng ($kg$) Lực va chạm đỉnh $F_{max}$ ($kN$) Năng lượng hấp thu riêng ($SEA - kJ/kg$) Hiệu suất lực va chạm ($CFE$)
Mẫu 01 (Ống vuông trơn) 3.45 185.4 12.8 48.2%
Mẫu 02 (Gân tăng cứng) 3.62 210.2 14.1 42.5%
Mẫu 03 (Đột lỗ Trigger) 3.38 152.0 15.3 61.4%
Mẫu 04 (Lượn sóng - Corrugated) 3.40 148.5 16.4 (+28.1%) 68.5% (+42.1%)

Hiệu suất lực va chạm ($CFE = F_{mean} / F_{max}$) của Mẫu 04 đạt $68.5%$, giúp triệt tiêu hiện tượng đỉnh lực nhọn gây tổn thương tức thời cho hành khách, đồng thời tăng tổng năng lượng hấp thu trên mỗi đơn vị khối lượng thêm $28.1%$.

Ứng dụng thực tế và Triển khai

+------------------------------------------------------------------------------------+
|                         KHẢ NĂNG ỨNG DỤNG CỦA ĐỀ TÀI                               |
|  [R&D Sản xuất Ô tô]    Cắt giảm 70% số lượng nguyên mẫu thử nghiệm phá hủy        |
|  [Đào tạo Đại học]       Hệ thống hóa giáo trình thực hành CAE/FEM chuyên sâu       |
|  [Giám định Tai nạn]     Tái dựng hiện trường va chạm phục vụ điều tra pháp y      |
|  [Hiệu quả Kinh tế]      Tiết kiệm 200.000 - 500.000 USD trong 1 chu kỳ phát triển  |
+------------------------------------------------------------------------------------+

1. Kịch bản ứng dụng trong công nghiệp ô tô (OEM & Tier-1 Suppliers)

  • Tối ưu hóa thiết kế ảo (Virtual Prototyping): Cho phép các kỹ sư thử nghiệm hàng trăm biến thể hình học và vật liệu khung vỏ (như thép siêu bền boron UHSS, hợp kim nhôm 6000/7000 series) trước khi tiến hành dập khuôn thực tế.
  • Chứng nhận an toàn số (Virtual Homologation): Giúp nhà sản xuất dự đoán trước điểm số đánh giá an toàn ASEAN NCAP hoặc Euro NCAP, rút ngắn chu kỳ đưa xe ra thị trường từ 36 tháng xuống còn 18-24 tháng.

2. Phục vụ công tác giám định pháp y và tái dựng tai nạn giao thông

  • Các cơ quan chức năng có thể sử dụng phương pháp mô phỏng số LS-DYNA để dựng lại hiện trường va chạm từ các dấu vết biến dạng cơ học thực tế, qua đó xác định chính xác tốc độ trước va chạm, góc va chạm và nguyên nhân lỗi kỹ thuật.

3. Lộ trình triển khai ứng dụng (Implementation Roadmap)

  • Giai đoạn 1 (Tháng 1 - 3): Xây dựng thư viện vật liệu chuẩn hóa (*MAT_024, *MAT_015, *MAT_100 cho mối hàn điểm Spot Weld).
  • Giai đoạn 2 (Tháng 4 - 6): Tích hợp quy trình mô phỏng tự động hóa vào nền tảng CAD/CAE của doanh nghiệp qua các script Python/LS-PrePost C-Command.
  • Giai đoạn 3 (Tháng 7 - 12): Triển khai thử nghiệm va chạm tích hợp khoang pin xe điện (EV Battery Pack Crashworthiness) chống đâm thủng gây cháy nổ nhiệt (Thermal Runaway).

Hạn chế và Hướng phát triển

Những hạn chế kỹ thuật hiện tại:

  • Mối liên kết hàn: Các liên kết giữa các tấm kim loại chủ yếu sử dụng liên kết cứng (Merged Nodes) hoặc mô hình hóa liên kết đơn giản, chưa tích hợp đầy đủ đặc tính phá hủy mỏi và tách mối hàn điểm (*MAT_SPOTWELD_DAMAGE).
  • Mô hình lốp và hệ thống treo: Hệ thống treo và áp suất lốp khí nén (*AIRBAG_HYBRID) được đơn giản hóa thành các liên kết dầm đàn hồi, ảnh hưởng một phần nhỏ đến động học quay thân xe sau va chạm.
  • Tài nguyên tính toán: Do giới hạn phần cứng máy trạm học đường, kích thước phần tử trung bình duy trì ở mức $5\text{ mm} - 10\text{ mm}$. Việc tinh chỉnh lưới dưới $2\text{ mm}$ tại các góc bo nếp gấp đòi hỏi hệ thống tính toán hiệu năng cao (HPC Cluster).

Hướng nghiên cứu và phát triển tiếp theo:

  1. Tích hợp mô hình hình nhân ảo toàn diện THUMS (Total Human Model for Safety) để phân tích chính xác mức độ tổn thương nội tạng, gãy xương và chỉ số chấn thương sọ não HIC.
  2. Ứng dụng thuật toán tối ưu hóa thông minh (Genetic Algorithm, Machine Learning Surrogate Models) để tự động hóa việc tìm kích thước tối ưu cho dầm hấp thu xung lực đa khoang (Multi-cell Crash Box).
  3. Mở rộng mô phỏng va chạm trên các dòng xe thuần điện (BEV), tập trung bảo vệ khối pin dưới sàn xe khi xảy ra va chạm cột bên hông (Side Pole Crash Test).

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

+---------------------------------------------------------------------------------------+
|                              ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                            |
+---------------------------------------------------------------------------------------+
|  Sinh viên Kỹ thuật Ô tô / Cơ khí:                                                   |
|  - Tiếp cận tài liệu chuẩn hóa về CAE phi tuyến và bộ giải LS-DYNA.                   |
|  - Nâng cao 45% kỹ năng thực hành phân tích động lực học kết cấu chuyên sâu.          |
+---------------------------------------------------------------------------------------+
|  Kỹ sư R&D và Thiết kế Khung vỏ (Body-in-White Engineers):                            |
|  - Sở hữu bộ thông số thẻ Keyword hoàn chỉnh (*MAT, *CONTACT, *CONTROL).              |
|  - Giảm thiểu 60% thời gian thiết lập bài toán mô phỏng ban đầu.                      |
+---------------------------------------------------------------------------------------+
|  Nhà sản xuất Ô tô & Doanh nghiệp phụ trợ:                                            |
|  - Cắt giảm hàng trăm nghìn USD chi phí phá hủy mẫu thử vật lý.                       |
|  - Nâng cao năng lực cạnh tranh và đáp ứng tiêu chuẩn xuất khẩu quốc tế (NCAP).       |
+---------------------------------------------------------------------------------------+
|  Nhà nghiên cứu An toàn Giao thông & Cơ quan Giám định:                              |
|  - Công cụ khoa học chính xác cao để phân tích cơ chế phá hủy và đề xuất luật an toàn.|
+---------------------------------------------------------------------------------------+

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

1. Yêu cầu cấu hình phần cứng tối thiểu để thực hiện mô phỏng LS-DYNA cho toàn bộ xe là gì?

Để giải quyết bài toán va chạm toàn xe (khoảng 300.000 - 1.000.000 phần tử) trong thời gian mô phỏng thực tế $100\text{ - }150\text{ ms}$: Cần máy trạm (Workstation) tối thiểu trang bị CPU 16-32 Cores (Intel Xeon hoặc AMD Threadripper), RAM tối thiểu 64GB DDR4/DDR5 ECC, ổ cứng SSD NVMe tốc độ cao để ghi nhận liên tục các tập tin xuất kết quả lớn (d3plot $> 20\text{GB}$).

2. Làm thế nào để kiểm soát hiện tượng biến dạng ảo (Hourglass Energy) trong LS-DYNA?

Hourglass là hiện tượng biến dạng bậc không (Zero-energy deformation mode) xảy ra ở các phần tử vỏ tích phân 1 điểm. Để kiểm soát:

  • Khai báo thẻ *CONTROL_HOURGLASS với loại $IHQ = 4$ (độ cứng theo công thức Flanagan-Belytschko) hoặc $IHQ = 6$ (dạng phi tuyến).
  • Giữ hệ số $QH$ ở mức $0.05 - 0.10$.
  • Luôn kiểm tra đồ thị năng lượng trong glstat: Đảm bảo tổng năng lượng Hourglass không vượt quá $5%$ tổng năng lượng toàn hệ thống.

3. Tại sao cần sử dụng đường cong ứng suất - biến dạng thực (True Stress-Strain) thay vì quy ước (Engineering Stress-Strain)?

Trong va chạm ô tô, kim loại trải qua biến dạng dẻo cực lớn (Large Strain / Large Deformation). Đường cong quy ước (Engineering) tính toán dựa trên diện tích tiết diện ban đầu $A_0$ không đổi, dẫn đến đánh giá sai khả năng chịu tải khi vật liệu bị thắt cổ (Necked). Đường cong thực (True Stress $\sigma_{true} = \sigma_{eng}(1 + \varepsilon_{eng})$) phản ánh chính xác diện tích tiết diện tức thời $A(t)$, đảm bảo độ chính xác của phản lực tiếp xúc.

4. Chi phí bản quyền phần mềm và khả năng tích hợp của LS-DYNA trong doanh nghiệp?

LS-DYNA cung cấp các gói bản quyền thương mại dạng Network/Floating License linh hoạt theo số lượng Core tính toán hoặc token trả phí theo giờ chạy đám mây (Ansys Cloud). Phần mềm tương thích hoàn hảo với các định dạng trung gian tiêu chuẩn như STEP, IGES, Parasolid và tích hợp trực tiếp với các hệ thống PLM hàng đầu thế giới.

5. Thời gian hoàn vốn (ROI) khi đầu tư giải pháp mô phỏng số va chạm là bao lâu?

Đối với một doanh nghiệp phát triển xe hoặc linh kiện dầm cản, chỉ cần loại bỏ được 01 lần thử nghiệm phá hủy xe mẫu vật lý (tiết kiệm khoảng 150.000 - 250.000 USD gồm chi phí chế tạo thủ công, cảm biến, bãi thử và xe tiêu hủy), doanh nghiệp đã có thể hoàn vốn toàn bộ chi phí đầu tư phần cứng máy trạm và bản quyền phần mềm ngay trong năm đầu tiên triển khai.

Kết luận

Đồ án tốt nghiệp "Ứng dụng phần mềm LS-DYNA mô phỏng và đánh giá mức độ biến dạng của ô tô khi va chạm" đã khẳng định vai trò tiên phong và tính ứng dụng vượt trội của kỹ thuật mô phỏng phần tử hữu hạn phi tuyến trong lĩnh vực an toàn ô tô hiện đại.

Bằng việc chuẩn hóa phương pháp luận từ cơ sở lý thuyết động lực học va chạm, thiết lập chính xác các điều kiện biên, mô hình vật liệu đàn - dẻo và tiếp xúc phức tạp trên phần mềm LS-DYNA, nhóm nghiên cứu đã:

  • Tái hiện chân thực và định lượng chính xác quá trình biến dạng, xung lực giảm tốc và mức độ xâm lấn khoang lái trên nhiều phân khúc phương tiện (Chevrolet C2500, Toyota Yaris, Toyota RAV4).
  • Chứng minh sự bất lợi nghiêm trọng về an toàn của các dòng xe đô thị cỡ nhỏ khi va chạm đối đầu với dòng xe bán tải có khối lượng và khoảng sáng gầm lớn.
  • Đề xuất thành công mẫu thiết kế dầm hấp thu xung lực trước dạng nếp sóng định hình (Corrugated Crash Box), giúp tăng hiệu suất hấp thu năng lượng riêng thêm $28.1%$ và nâng cao chỉ số hiệu quả lực va chạm lên $68.5%$.

Kết quả nghiên cứu không chỉ mang lại giá trị học thuật cao cho chương trình đào tạo kỹ sư ô tô chất lượng cao mà còn là tài liệu tham khảo kỹ thuật giá trị cho các trung tâm nghiên cứu và phát triển phương tiện giao thông, hướng tới mục tiêu nâng cao tiêu chuẩn an toàn cho người tham gia giao thông tại Việt Nam và khu vực.