Giải Mã Động Lực Học Phá Hủy: Mô Phỏng Vụ Nổ Vật Liệu LX-17 Trong Hình Trụ 3D Bằng Bộ Giải CESE Của LS-DYNA

Trong lĩnh vực mô phỏng động lực học va chạm và nổ hiện nay, một trong những bài toán kinh điển và thử thách nhất đối với các bộ giải số (solver) là mô phỏng quá trình một khối thuốc nổ hình thanh (slab) được bao bọc bởi một lớp vỏ kim loại, phát nổ và gây biến dạng vỏ chứa. Đây là nền tảng để hiểu rõ về hiệu suất của đầu đạn, kiểm soát bán kính sát thương, và thiết kế các cấu trúc chịu nổ.

Một trường hợp kết quả đánh giá điển hình trong lĩnh vực này là sóng xung kích trong ống. Điểm đặc biệt là thuốc nổ được kích nổ tại một đầu bởi một kíp nổ mồi (booster), tạo ra sóng nổ di chuyển (propagate) dọc theo chiều dài khối thuốc. Áp lực từ phía sau của sóng nổ đủ mạnh để khiến cho lớp vỏ kim loại bọc bên ngoài bị biến dạng dẻo, phồng rộng, thậm chí có thể dẫn đến vỡ vụn.

Bài toán này là phép thử hoàn hảo để kiểm tra năng lực của một công nghệ giải số tiên tiến trong LS-DYNA, sử dụng bộ giải CESE (Conservation Element & Solution Element), một phương pháp lưới cấu trúc dành riêng cho các bài toán động lực học chất lưu có sóng xung kích mạnh.

1. Thách thức Kỹ thuật: Sóng Nổ Di Động và Phản Ứng Vật Liệu

Mô phỏng quá trình sóng nổ lan truyền trong khối thuốc và tương tác với vỏ kim loại đặt ra ba thách thức mang tính định luật:

• Tính di động của sóng nổ (Detonation Wave Propagation): Sóng nổ không phải là một vụ nổ tức thời toàn khối. Nó bắt đầu từ một điểm kích thích (booster) và di chuyển với vận tốc siêu thanh trong lòng khối thuốc. Bộ giải cần mô phỏng chính xác khoảnh khắc này để xác định vùng áp suất cực cao tác dụng lên vỏ chứa tại từng thời điểm.
• Động lực học phản ứng hóa học: Không phải cứ có nổ là áp suất tăng ngay lập tức. Cần phải mô hình hóa "phản ứng cháy nổ" — quá trình chuyển đổi từ vật liệu chưa phản ứng sang sản phẩm khí. Bài toán sử dụng định luật tỷ lệ phản ứng Ignition and Growth (I&G) để mô phỏng giai đoạn khởi cháy và giai đoạn phát triển của vụ nổ.
• Phương pháp tiếp cận đa vật lý (Multiphase Flow): Vụ nổ tạo ra một dòng chảy hai pha phức tạp (khí nổ và mảnh vỡ). Sử dụng phương pháp Lagrange truyền thống cho khối thuốc nổ sẽ thất bại do lưới bị xoắn vặn cực độ. Cần một bộ giải có khả năng xử lý dòng chảy Euler có phản ứng hóa học.

2. Giải pháp Công nghệ trong LS-DYNA: Bộ Giải DUALCESE

Để giải quyết bài toán này, LS-DYNA cung cấp công nghệ DUALCESE — một bộ giải dòng chảy đa pha với hệ thống từ khóa riêng biệt, được tối ưu hóa để mô phỏng sóng xung kích trong chất lưu và khí:

Mô hình hóa Thuốc nổ và Vỏ chứa với CESE

Thay vì các phần tử hữu hạn khối thông thường, bài toán sử dụng DUALCESE_EOS_JWL, và quan trọng nhất là sử dụng mô hình tỷ lệ phản ứng Ignition and Growth (I&G) (*DUALCESE_REACTION_RATE_IG). Mô hình này mô phỏng cách vật liệu chưa nổ chuyển hóa thành sản phẩm khí sau khi sóng xung kích đi qua một lưới. Khối thuốc nổ được khai báo là Part đa pha (Multiphase Part) thông qua từ khóa *DUALCESE_PART_MULTIPHASE. Chất nổ (thuốc LX-17) được định nghĩa bởi phương trình trạng thái JWL (Jones-Wilkins-Lee).

Phương pháp Điều khiển Thời gian và Xuất Kết quả

Sóng nổ di chuyển cực nhanh, đòi hỏi bước thời gian cực nhỏ. LS-DYNA sử dụng *DUALCESE_CONTROL_TIMESTEP với hệ số CFL (Courant-Friedrichs-Lewy) để tự động điều chỉnh bước nhảy, đảm bảo sự ổn định của thuật toán.

3. Kết quả Mô phỏng: Tiến trình Sóng Nổ theo Thời gian

Dựa trên cấu hình mẫu, quá trình xảy ra theo các pha như sau:

1. Pha kích nổ (T = 0 µs): Kíp nổ mồi (booster) được kích hoạt tại đầu trái của khối thuốc, giải phóng một lượng lớn năng lượng cục bộ.
2. Pha lan truyền (0 – 7.5 µs): Sóng nổ hình thành và di chuyển dọc theo thanh thuốc về phía bên phải với vận tốc siêu thanh. Phía sau mặt sóng là vùng khí nổ có áp suất lên tới hàng chục Gigapascal.
3. Pha gia tải cho vỏ (> 7.5 µs): Khi sóng nổ đi qua từng đoạn của vỏ hình trụ, áp suất cao tác dụng lên thành trong. Vỏ kim loại bắt đầu phình ra, chịu ứng suất kéo cực lớn trong một khoảng thời gian cực ngắn.
4. Pha biến dạng tối đa (15 µs – Kết thúc mô phỏng): Vỏ hình trụ bị phồng rộng hoàn toàn. Dữ liệu xuất ra dạng *DATABASE_BINARY_D3PLOT cho phép dựng lại hình ảnh 3D của vỏ bị biến dạng sau vụ nổ.

4. Lợi ích Chiến lược cho R&D và An ninh

Thông qua bài toán kích nổ trong ống, LS-DYNA mang lại những giá trị cốt lõi vượt trội:

• Hiểu sâu về cơ chế lan truyền nổ: Các kỹ sư có thể hình dung rõ ràng hình dạng sóng xung kích với mô hình thuốc nổ phức tạp, từ đó thiết kế hệ thống kích nổ (mạng lưới dây nổ) hiệu quả nhất.
• Tối ưu hóa hiệu suất gây nổ: Xác định chính xác áp lực tác động lên vỏ đạn tại từng thời điểm, giúp tối ưu hóa độ dày vỏ chứa, tiết kiệm vật liệu trong khi vẫn đảm bảo phân bố mảnh vỡ như mong muốn.
• Nền tảng cho mô phỏng cấp cao: Kết quả thu được từ bộ giải CESE có thể được sử dụng làm đầu vào cho các bài toán thứ cấp khác trong LS-DYNA như tương tác giữa mảnh vỡ và kết cấu (FSI) hoặc tối ưu hóa đa mục tiêu.

5. Kết luận

Mô phỏng thử nghiệm Slab Stick 3D bằng bộ giải CESE trong LS-DYNA là một minh chứng rõ ràng cho sức mạnh của phương pháp số trong việc giải quyết các bài toán động lực học phá hủy phức tạp nhất. Từ việc mô phỏng tỷ lệ phản ứng hóa học I&G, phương trình trạng thái JWL cho khí nổ, đến việc sử dụng lưới cấu trúc để tối ưu tốc độ tính toán, công nghệ này đã giúp các kỹ sư mô phỏng thành công hiện tượng sóng nổ di động. Nhờ đó, các nhà phát triển khí tài quân sự không chỉ hiểu sâu sắc về cơ chế nổ bên trong lòng kết cấu, mà còn có thể định hình được những giải pháp phòng thủ và công nghệ tiên tiến, chính xác hơn cho tương lai.