Giới thiệu dự án

Tai nạn giao thông đường bộ tiếp tục là một trong những thách thức nghiêm trọng nhất đối với xã hội hiện đại. Theo thống kê của Tổ chức Y tế Thế giới (WHO) và Ủy ban An toàn Giao thông Quốc gia, va chạm trực diện (frontal collision) chiếm hơn 50% tổng số các vụ tai nạn đặc biệt nghiêm trọng, dẫn đến tỷ lệ tử vong và đa chấn thương vùng đầu, cổ và lồng ngực cao nhất ở người ngồi trong xe. Để giảm thiểu rủi ro chấn thương, hệ thống an toàn bị động (Passive Safety Systems), đặc biệt là hệ thống túi khí bổ trợ (Supplemental Restraint System - SRS), đóng vai trò sống còn trong việc hấp thụ xung lực và phân tán động năng va chạm.

Tuy nhiên, việc nghiên cứu, phát triển và đánh giá hiệu năng làm việc của túi khí thông qua thử nghiệm phá hủy thực tế (Physical Crash Testing) đòi hỏi chi phí khổng lồ (hàng trăm nghìn USD cho mỗi chu kỳ thử nghiệm), tốn kém thời gian và khó khăn trong việc thu thập trường dữ liệu động lực học vi mô. Đồ án tốt nghiệp "Nghiên cứu mô phỏng hoạt động túi khí trên ô tô" (Chuyên ngành: Công nghệ Kỹ thuật Ô tô, Trường Đại học Văn Lang, Tác giả: Nguyễn Hoàng Vũ, Cán bộ hướng dẫn: PGS. Nguyễn Phụ Thượng Lưu) được triển khai nhằm giải quyết bài toán tối ưu hóa an toàn bị động bằng công nghệ mô phỏng số (Computer-Aided Engineering - CAE).

                      HỆ THỐNG AN TOÀN BỊ ĐỘNG TRÊN Ô TÔ

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

  1. Xây dựng lý thuyết động học và nhiệt động học: Khái quát hóa quá trình kích nổ pyrotechnic, cơ chế sinh khí Nitơ ($N_2$) từ hợp chất Sodium Azide ($NaN_3$) và động lực học giãn nở túi khí.
  2. Thiết lập mô hình 3D CAD/CAE hoàn chỉnh: Xây dựng khoang lái, vô lăng, túi khí người lái (driver airbag module) và tích hợp mô hình người nộm tiêu chuẩn Hybrid III 50th Percentile Male Dummy.
  3. Mô hình hóa phi tuyến phần tử hữu hạn: Thiết lập tương tác tiếp xúc phức tạp (Fluid-Structure Interaction - FSI) và ứng xử cơ học phi tuyến của vải túi khí trong bộ giải LS-DYNA Explicit Solver.
  4. Mô phỏng va chạm trực diện: Thực hiện kiểm thử số kịch bản xe đâm trực diện vào tường cứng ở vận tốc $50\text{ km/h}$ ($13.89\text{ m/s}$) theo tiêu chuẩn an toàn liên bang Hoa Kỳ FMVSS 208.
  5. Định lượng hiệu quả bảo vệ sinh cơ học: Đo lường, trích xuất biểu đồ lực tác động lên ngực, gia tốc đỉnh vùng đầu, phân tích mức giảm thương tổn giữa hai trạng thái: có và không có trang bị túi khí.

Phạm vi và giới hạn của đề tài

  • Đối tượng nghiên cứu: Cụm mô-đun túi khí bảo vệ người lái phía trước (Frontal Driver Airbag System) lắp đặt trong vô lăng xe du lịch con.
  • Kịch bản va chạm: Va chạm trực diện đối xứng $100%$ diện tích đầu xe vào tường cứng không biến dạng (Rigid Wall Barrier) ở dải vận tốc $50\text{ km/h}$.
  • Thời gian phân tích: Khung thời gian mô phỏng giới hạn trong $100\text{ ms}$ – giai đoạn vàng quyết định tính mạng hành khách khi xảy ra va chạm.
  • Giới hạn kỹ thuật: Không can thiệp sâu vào thuật toán lọc tín hiệu số của vi điều khiển ECU/ACU mà giả định cảm biến kích hoạt tức thời tại ngưỡng gia tốc âm $a \ge 4g$ (độ trễ danh định $10\text{--}12\text{ ms}$).

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

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

Trong quy trình phát triển phương tiện, các phương pháp đánh giá an toàn thụ động hiện nay được chia thành ba nhánh tiếp cận chính:

Tiêu chí phân tích Thử nghiệm thực tế (Crash Test) Mô hình giải tích 1D (MATLAB/Simulink) Mô phỏng số 3D CAE (LS-DYNA Explicit)
Độ chính xác vật lý Rất cao (Thực tế $100%$) Trung bình (Rút gọn khối lượng - lò xo) Rất cao (Động lực học phi tuyến 3D)
Chi phí thực hiện Cực kỳ tốn kém ($> 100,000\text{ USD}$/lần) Thấp ($< 500\text{ USD}$) Tiết kiệm ($> 90%$ so với thử nghiệm thật)
Thời gian thiết lập $2\text{--}4$ tuần chuẩn bị $1\text{--}2$ ngày $3\text{--}5$ ngày lập mô hình
Khả năng lặp lại Phá hủy mẫu, không lặp lại nguyên bản Cao, nhưng thiếu hình học không gian Tuyệt đối, dễ dàng thay đổi biến số
Tương tác vải - khí Quan sát ngoài bằng camera tốc độ cao Không mô hình hóa được trường áp suất Chi tiết từng phần tử (FEM/CPM Method)

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must have (Bắt buộc): Mô hình hóa vật liệu màng mỏng phi tuyến *MAT_FABRIC (Nylon 66 phủ silicon); thuật toán tích phân thời gian tường minh (Explicit Integration); mô hình tiếp xúc không xuyên thấu giữa túi khí và dummy; kịch bản va chạm $50\text{ km/h}$.
  • Should have (Nên có): Phương pháp mô hình hóa hạt khí Corpuscular Particle Method (CPM) hoặc mô hình áp suất đồng nhất Uniform Pressure (UP) chính xác; đánh giá tiêu chí chấn thương sọ não (Head Injury Criterion - HIC).
  • Could have (Có thể mở rộng): Tích hợp thêm hệ thống dây đai an toàn căng đai tự động (Seatbelt Pretensioner); mô hình hóa van xả khí (vent holes) điều áp.
  • Won't have (Chưa thực hiện đợt này): Mô phỏng va chạm góc lệch tâm (Small Overlap Barrier $25%$), lật xe (Rollover), hoặc mô phỏng mô mềm người thật (Human Body Model - THUMS).

Thiết kế hệ thống mô phỏng

Kiến trúc quy trình kỹ thuật số (CAE Pipeline) được thiết kế qua 4 giai đoạn chuẩn hóa:

 CAD Modeling (SolidWorks) 
 HyperMesh Meshing (2D Shell / 3D Solid)
 LS-DYNA Explicit Solver (Non-linear Dynamic Engine)
 LS-PrePost (Biomechanical Injury Post-processing)

Technology Stack và Công cụ sử dụng

  • CAD Geometry Engine: SolidWorks 2023 (Dựng hình học khoang lái, cột lái, vô lăng, vỏ cụm túi khí).
  • CAE Pre-processor: Altair HyperMesh 2022 (Chia lưới phần tử hữu hạn 2D Shell, kiểm soát độ méo Jacobian $\ge 0.7$, tỷ lệ co Aspect Ratio $\le 5$).
  • FEA Solver: Ansys LS-DYNA SMP/MPP Version R14.0 (Bộ giải động lực học phi tuyến tường minh).
  • Post-processor: LS-PrePost Version 4.10 (Trích xuất vector gia tốc, ứng suất Von Mises, đường cong phản lực).
  • Biomechanical Target: Mô hình chuẩn hóa Hybrid III 50th Percentile Male Dummy (Khối lượng danh định $75\text{ kg}$, chiều cao $1.75\text{ m}$).

LS-DYNA Keyword Architecture & Data Flow

Cấu trúc file đầu vào (.k / .dyn) được cấu hình phân lớp nghiêm ngặt nhằm đảm bảo độ hội tụ:

*KEYWORD
*CONTROL_TERMINATION
$  ENDTIM    ENDCYC     DTMIN    ENDENG    ENDMAS
    100.0         0       0.0       0.0       0.0
*CONTROL_TIMESTEP
$  DTINIT    TSSFAC      ISDO    TSLIMT     DTMS
     0.00      0.90         0      0.00 -1.0E-06
*MAT_FABRIC
$   MID        RO         E        PR        EA        EB      GAB
      1  1.40E-06     2.000      0.35     1.200     1.200    0.450
*AIRBAG_SIMPLE_AIRBAG_MODEL
$  ABID      CV        CP         T        PE      PAIRC
      1   717.0    1004.0     600.0     1.0E5     1.01E5
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$  SSID      MSID     SSTYP     MSTYP     SBOXID    MBOXID       SPR       MPR
      2         1         3         3          0         0       1.0       1.0
*END

Implementation và kết quả

Development Process & Cơ sở thuật toán

Thuật toán giải chuyển vị trong LS-DYNA dựa trên phương pháp tích phân sai phân trung tâm (Explicit Central Difference Time Integration Scheme). Phương trình cân bằng động lực học phi tuyến tại bước thời gian $t_n$:

$$\mathbf{M} \ddot{\mathbf{u}}_n + \mathbf{C} \dot{\mathbf{u}}_n + \mathbf{F}^{int}(\mathbf{u}_n) = \mathbf{F}^{ext}_n$$

Trong đó:

  • $\mathbf{M}$: Ma trận khối lượng đường chéo (Lumped Mass Matrix).
  • $\ddot{\mathbf{u}}_n, \dot{\mathbf{u}}_n, \mathbf{u}_n$: Vector gia tốc, vận tốc và chuyển vị của các nút phần tử.
  • $\mathbf{F}^{int}, \mathbf{F}^{ext}$: Vector nội lực (tính từ tensor ứng suất Cauchy) và ngoại lực.

Gia tốc nút tại bước $n$ được giải trực tiếp mà không cần đảo ma trận độ cứng:

$$\ddot{\mathbf{u}}_n = \mathbf{M}^{-1} \left( \mathbf{F}^{ext}_n - \mathbf{F}^{int}_n - \mathbf{C} \dot{\mathbf{u}}_n \right)$$

Vận tốc và chuyển vị tại các bước kế tiếp được tích phân qua công thức:

$$\dot{\mathbf{u}}{n+1/2} = \dot{\mathbf{u}}{n-1/2} + \ddot{\mathbf{u}}_n \Delta t_n$$

$$\mathbf{u}{n+1} = \mathbf{u}n + \dot{\mathbf{u}}{n+1/2} \Delta t{n+1/2}$$

Để bài toán không bị mất ổn định số, bước thời gian tính toá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} = \frac{L_e}{\sqrt{\frac{E}{\rho (1-\nu^2)}}}$$

(Trong đó $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, $E$ là mô-đun Young, $\rho$ là mật độ khối lượng, $\nu$ là hệ số Poisson).

Testing và Validation

Mô phỏng được chạy song song trên trạm làm việc hiệu năng cao (HPC Workstation: 16 Cores Intel Xeon, 64GB DDR4 ECC RAM, GPU Acceleration).

                 DIỄN BIẾN ĐỘNG HỌC VA CHẠM (0 - 100 ms)

Kết quả định lượng đạt được

Kết quả mô phỏng số chiết xuất từ LS-PrePost minh chứng rõ nét khả năng bảo vệ người lái vượt bậc của hệ thống túi khí:

Chỉ số đo lường sinh cơ học Trường hợp 1: KHÔNG có túi khí Trường hợp 2: CÓ túi khí (LS-DYNA) Mức độ cải thiện / Giảm thiểu Ngưỡng an toàn FMVSS 208
Lực tác động lên ngực (Chest Load) $4.50\text{ kN}$ $1.80\text{ kN}$ Giảm 60.0% xung lực $\le 3.0\text{ kN}$ (Đạt)
Gia tốc đỉnh vùng đầu (Head Peak Acc.) $62.5\text{ g}$ $30.2\text{ g}$ Giảm 51.6% gia tốc âm $\le 60.0\text{ g}$ (Đạt)
Thời gian bắt đầu tiếp xúc vật lý $25\text{ ms}$ (Đập vào vô lăng) $45\text{ ms}$ (Chạm đệm khí mềm) Kéo dài thời gian va chạm Triệt tiêu xung kích
Khoảng cách dịch chuyển đầu tối đa Va chạm trực diện vành lái Cách vô lăng an toàn $120\text{ mm}$ Triệt tiêu nguy cơ chấn thương nứt sọ An toàn tuyệt đối
Tỷ lệ bảo tồn thể tích khoang lái Biến dạng cột lái $35\text{ mm}$ Hấp thụ biến dạng mềm $8\text{ mm}$ Giảm phá hủy thứ cấp Đạt chuẩn Euro NCAP

Đổi mới và đóng góp

  1. Ứng dụng phương pháp hạt Corpuscular Particle Method (CPM): Thay vì sử dụng mô hình áp suất đồng nhất đơn giản (Uniform Pressure) vốn xem nhẹ tính động học của dòng khí, đề tài đã tiếp cận và triển khai phương pháp hạt khí CPM trong LS-DYNA. Điều này giúp mô tả chính xác áp suất cục bộ và sóng áp suất va đập lên màng vải túi khí trong $15\text{ ms}$ đầu tiên sau kích nổ.
  2. Chuẩn hóa quy trình liên kết CAD-CAE cho an toàn ô tô: Xây dựng hoàn chỉnh chuỗi dữ liệu kỹ thuật từ mô hình hình học SolidWorks, chia lưới tối ưu qua HyperMesh đến thiết lập thẻ từ khóa chuyên sâu trong LS-DYNA, tạo cơ sở thực hành mẫu cho sinh viên và kỹ sư thiết kế ô tô.
  3. Định lượng toàn diện tổn thương sinh cơ học: Khảo sát chi tiết đáp ứng lực và gia tốc của dummy Hybrid III theo thời gian thực ($0\text{--}100\text{ ms}$), cung cấp dữ liệu định lượng chính xác (giảm $60%$ lực ngực, giảm $51.6%$ gia tốc đầu) thay vì chỉ dừng lại ở các phân tích định tính mô tả chuyển động.
  4. Đóng góp học thuật và đào tạo kỹ thuật: Bổ sung nguồn tài liệu chuyên sâu bằng tiếng Việt có cấu trúc bài bản, phục vụ công tác đào tạo kỹ sư công nghệ kỹ thuật ô tô, giảm phụ thuộc vào các tài liệu ngoại ngữ rời rạc.

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

Kịch bản ứng dụng công nghiệp

  • Nghiên cứu & Phát triển (R&D) tại các OEM ô tô: Áp dụng mô hình để tối ưu hóa hình dạng túi khí, vị trí đặt lỗ thoát khí (vent holes), và lựa chọn loại vật liệu dệt túi khí trước khi tiến hành dập khuôn chế tạo hàng loạt.
  • Tối ưu hóa hệ thống an toàn thụ động theo chuẩn NCAP: Hỗ trợ tinh chỉnh thời điểm nổ của ngòi nổ pyrotechnic ($10\text{--}15\text{ ms}$) kết hợp với độ giãn dài của dây đai an toàn nhằm tối đa hóa điểm số thử nghiệm va chạm Euro NCAP và Latin NCAP.
  • Nội địa hóa linh kiện ngành công nghiệp ô tô: Phục vụ các doanh nghiệp phụ trợ tại Việt Nam trong việc thử nghiệm ảo độ bền và khả năng bung của các chi tiết vô lăng, táp-lô chứa mô-đun túi khí.
       LỘ TRÌNH TRIỂN KHAI VÀ CHUYỂN GIAO CÔNG NGHỆ CAE

Yêu cầu phần cứng và triển khai hệ thống

  • CPU: Tối thiểu 8 Cores / 16 Threads (Khuyến nghị AMD EPYC hoặc Intel Xeon chuyên dụng).
  • Bộ nhớ RAM: Tối thiểu 32GB ECC RAM (Khuyến nghị 64GB đối với các mô hình $> 500,000$ phần tử).
  • Lưu trữ: SSD NVMe dung lượng trống $\ge 500\text{ GB}$ để ghi log tệp kết quả động lực học d3plot.
  • Hệ điều hành: Linux RedHat Enterprise / Ubuntu LTS hoặc Windows Server 64-bit.

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

Hạn chế kỹ thuật hiện tại

  • Giả định điều kiện biên cứng: Mô hình giả định tường chắn hoàn toàn cứng (Rigid Wall) và không tính đến sự hấp thụ biến dạng của toàn bộ khung xe phía trước (Crumple Zone), dẫn đến xung lực truyền vào sàn xe có dạng dốc đứng hơn thực tế.
  • Mô hình nhiệt động học khí nén lý tưởng: Bỏ qua tổn thất nhiệt truyền qua bề mặt vải túi khí ra môi trường khoang lái trong quá trình cháy của hạt sinh khí $NaN_3$.
  • Cảm biến va chạm đơn điểm: Kích hoạt túi khí được giả định theo thời gian trễ cố định ($12\text{ ms}$) thay vì xử lý dòng dữ liệu liên tục từ nhiều cảm biến gia tốc và áp suất sườn theo thời gian thực.

Hướng phát triển và mở rộng

  • Tích hợp mô hình cơ thể người sống (Human Body Model - HBM): Chuyển đổi từ Dummy Hybrid III sang mô hình THUMS (Total Human Model for Safety) để đánh giá chi tiết gãy xương sườn, tổn thương dập phổi và xuất huyết nội tạng.
  • Mở rộng các kịch bản va chạm phức tạp: Mô phỏng va chạm lệch tâm góc nhỏ ($25%$ Small Overlap Frontal Crash), va chạm hông xe (Side Pole Impact) và lật xe liên hoàn.
  • Ứng dụng trí tuệ nhân tạo (AI/ML): Sử dụng mạng nơ-ron học sâu để dự đoán biến dạng túi khí tức thời (Surrogate Modeling), rút ngắn thời gian tính toán tối ưu hóa từ nhiều giờ xuống vài giây.

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

  • Sinh viên ngành Kỹ thuật Ô tô & Cơ điện tử: Cung cấp tài liệu mẫu mực về phân tích phần tử hữu hạn phi tuyến, nắm bắt phương pháp thiết lập bài toán va chạm thực tế trong LS-DYNA.
  • Kỹ sư CAE & Thiết kế an toàn xe: Cung cấp bộ thông số tham chiếu chuẩn xác (*MAT_FABRIC, *AIRBAG_PARTICLE) và kỹ thuật kiểm soát năng lượng mô phỏng, hạn chế tối đa lỗi trượt tiếp xúc hoặc mất ổn định bước thời gian.
  • Các doanh nghiệp sản xuất và lắp ráp ô tô: Tiết kiệm hàng triệu USD chi phí thử nghiệm vật lý trong giai đoạn Concept, đẩy nhanh tiến độ đưa sản phẩm đạt chuẩn an toàn ra thị trường.
  • Các nhà nghiên cứu cơ học va chạm: Nền tảng dữ liệu mở để tiếp tục phát triển các mô hình vật liệu dệt thông minh có khả năng thay đổi độ thấm khí theo vận tốc va đập.

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

1. Yêu cầu kỹ thuật tối thiểu để chạy mô phỏng bài toán túi khí trong LS-DYNA là gì?

Hệ thống cần tối thiểu CPU 4 nhân (tần số $\ge 3.0\text{ GHz}$), 16GB RAM và bộ cài LS-DYNA R11 trở lên. Đối với mô hình độ phân giải lưới cao ($> 100,000$ phần tử), khuyến nghị trạm làm việc 16 nhân, 64GB RAM và ổ đĩa SSD NVMe tốc độ cao để đảm bảo thời gian chạy hoàn tất trong vòng $1\text{--}3$ giờ.

2. Làm thế nào để xử lý lỗi mất ổn định bước thời gian (Time Step Drop / Instability) khi túi khí biến dạng lớn?

Nguyên nhân chủ yếu do phần tử vải bị bóp méo quá mức khiến chiều dài đặc trưng $L_e$ tiến về 0. Giải pháp gồm: áp dụng kỹ thuật kiểm soát biến dạng Hourglass Control (*CONTROL_HOURGLASS Type 6 hoặc 8 cho vật liệu vải), sử dụng tính năng Mass Scaling có kiểm soát (*CONTROL_TIMESTEP, giá trị DTMS âm thích hợp) và kích hoạt chế độ tự động định nghĩa lại lưới tiếp xúc (*CONTACT_AUTOMATIC_SINGLE_SURFACE).

3. Phương pháp Corpuscular Particle Method (CPM) ưu việt hơn Uniform Pressure (UP) như thế nào?

Mô hình UP giả định áp suất bên trong túi khí đồng nhất tại mọi điểm ngay khi bung, dẫn đến đánh giá sai lực va đập ban đầu. Phương pháp CPM mô phỏng hàng triệu hạt khí động lực học chuyển động và phản xạ thực tế, giúp tính toán chính xác hiện tượng "túi khí tát vào mặt người lái" (Out-of-Position impact) trong $10\text{--}20\text{ ms}$ đầu tiên.

4. Dữ liệu từ mô phỏng số có đủ độ tin cậy để thay thế hoàn toàn thử nghiệm vật lý không?

Mô phỏng CAE đóng vai trò cắt giảm $80\text{--}90%$ số lượng mẫu thử nghiệm vật lý trung gian trong quá trình phát triển. Tuy nhiên, theo quy định của các cơ quan quản lý giao thông quốc tế (NHTSA, UNECE), thử nghiệm vật lý cuối cùng (Homologation Crash Test) vẫn là bắt buộc để cấp phép lưu hành thương mại cho phương tiện.

5. Chi phí đầu tư cho hệ thống mô phỏng CAE và thời gian hoàn vốn (ROI) ước tính như thế nào?

Đầu tư ban đầu bao gồm bản quyền phần mềm CAE và máy trạm tính toán ($30,000\text{--}70,000\text{ USD}$). Tuy nhiên, chỉ cần giảm thiểu được $2\text{--}3$ mẫu xe thử nghiệm phá hủy thật (chi phí trung bình $50,000\text{ USD}$/xe thử nghiệm), doanh nghiệp đã có thể đạt điểm hòa vốn ngay trong năm đầu tiên phát triển sản phẩm mới.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu mô phỏng hoạt động túi khí trên ô tô" đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra. Thông qua việc ứng dụng công cụ mô phỏng động lực học phi tuyến LS-DYNA và chuẩn hóa mô hình hình nhân Hybrid III, nghiên cứu đã chứng minh một cách thuyết phục:

  • Hệ thống túi khí người lái giúp giảm tới 60.0% xung lực tác động lên lồng ngực (từ $4.5\text{ kN}$ xuống $1.8\text{ kN}$).
  • Giảm 51.6% gia tốc đỉnh vùng đầu (từ $62.5\text{ g}$ xuống $30.2\text{ g}$), đưa các chỉ số thương tổn về sâu dưới ngưỡng nguy hiểm của tiêu chuẩn an toàn liên bang FMVSS 208.

Công trình không chỉ mang ý nghĩa lý thuyết sâu sắc trong việc làm chủ công nghệ mô phỏng số tiên tiến mà còn mở ra tiềm năng ứng dụng thực tiễn to lớn cho ngành công nghiệp ô tô nội địa tại Việt Nam, hướng tới tối ưu hóa an toàn giao thông và giảm thiểu thương vong cho người sử dụng phương tiện. Các kỹ sư, sinh viên và nhà nghiên cứu quan tâm có thể tiếp tục ứng dụng bộ thông số và quy trình trong đồ án để mở rộng cho các bài toán an toàn giao thông đa chiều trong tương lai.