Tổng quan về luận án

Công trình nghiên cứu "Phân tích động lực học vỏ trụ có gân gia cường trên liên kết đàn hồi chịu tác dụng của hệ sóng xung kích" do nghiên cứu sinh Lê Xuân Thùy thực hiện dưới sự hướng dẫn khoa học của GS.TS Nguyễn Thái Chung tại Học viện Kỹ thuật Quân sự (chuyên ngành Cơ kỹ thuật, mã số: 9.01) đã giải quyết bài toán cốt lõi trong cơ học công trình quốc phòng và dân dụng hiện đại. Trong điều kiện chiến tranh công nghệ cao và các nguy cơ mất an toàn công nghiệp, các công trình quân sự ngầm hoặc bán ngầm như vòm hầm trú ẩn, công sự bảo vệ khí tài, mái che tên lửa và kho tàng vũ khí chịu nguy cơ cao từ các cuộc tập kích đường không với chuỗi vụ nổ liên tiếp.

Khoảng trống nghiên cứu (research gap) được luận án nhận diện rất rõ ràng: hầu hết các công trình trong và ngoài nước (như Türkmen [58], Nguyễn Đức Thắng [33]) chủ yếu giới hạn ở việc khảo sát sóng xung kích đơn lẻ tác dụng lên tấm hoặc vỏ trơn với các điều kiện biên cứng lý tưởng (ngàm hoặc tựa đơn). Luận án trích dẫn trực tiếp: "Với sóng xung kích đơn (mô tả cho một vụ nổ) tác dụng lên kết cấu công trình là loại bài toán đến nay đã được xem xét, giải quyết khá đầy đủ, còn với trường hợp hệ sóng xung kích tác dụng (các vụ nổ liên tiếp), do tính phức tạp của nó nên đến nay các công bố còn hạn chế." Đồng thời, hiện tượng chồng lấn tải trọng động (dynamic load superposition) khi các lớp sóng va đập kế tiếp tác động lên một hệ kết cấu đang dao động phi tuyến có gân gia cường và gối tựa đàn hồi vẫn chưa có lời giải số kết hợp kiểm chứng thực nghiệm hiện trường hoàn chỉnh.

Để giải quyết khoảng trống này, luận án xây dựng hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật ứng xử phi tuyến hình học của vỏ trụ thoải có gân trực giao trên liên kết đàn hồi chịu tải trọng chuỗi sóng nổ biến thiên theo thời gian diễn ra như thế nào?
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ ảnh hưởng tương tác giữa các tham số kết cấu (tỷ số kích thước gân $h_g/b_g$, bán kính cong $R$, chiều dày vỏ $h$, độ cứng lò xo $k$) và tham số tải trọng (khoảng thời gian trễ $\Delta t_0$, số lượng lớp sóng, siêu áp cực đại $p_m$) đến chuyển vị, gia tốc và trạng thái ứng suất - biến dạng?
  3. Giả thuyết khoa học (H1): Việc kết hợp hệ gân gia cường trực giao với các gối tựa đàn hồi cho phép phân tán năng lượng nổ, giảm đáng kể biên độ chuyển vị và ứng suất uốn cục bộ so với kết cấu liên kết cứng; tuy nhiên, tồn tại các giá trị thời gian trễ $\Delta t_0$ cộng hưởng gây mất ổn định động lực học nghiêm trọng.

Khung lý thuyết của luận án được phát triển dựa trên nền tảng lý thuyết vỏ thoải Reissner - Mindlin (First-order Shear Deformation Theory - FSDT) có xét đến phi tuyến hình học von Kármán, mô hình tương tác sóng - kết cấu của M.A. Sadovsky và phương pháp phần tử hữu hạn (PTHH) biến phân Galerkin. Nghiên cứu thực hiện khảo sát số trên miền kết cấu vỏ kích thước $L \times W$, bề dày $h$, bán kính cong $R$, cùng với thực nghiệm hiện trường sử dụng lượng nổ TNT $50\text{ g}$ tại Học viện Kỹ thuật Quân sự, tạo tiền đề định lượng chuẩn xác cho việc tối ưu hóa kết cấu công sự phòng thủ.

Literature Review và Positioning

Tổng quan y văn cho thấy các nghiên cứu về kết cấu tấm và vỏ chịu tải trọng nổ đã hình thành qua nhiều trường phái học thuật với sự phát triển của cả phương pháp giải tích và phương pháp số:

  1. Trường phái giải tích và bán giải tích: Türkmen [58, 59] thiết lập phương trình vi phân chuyển động của panel vỏ trụ composite nhiều lớp chịu tải trọng sóng xung kích đơn và giải bằng phương pháp Runge - Kutta. Nhóm nghiên cứu Đào Huy Bích, Đào Văn Dũng, Vũ Hoài Nam [43, 45, 77] và Nguyễn Đình Đức [72, 73] đã có những đóng góp cơ bản trong việc giải tích phi tuyến vỏ FGM và vỏ có gân gia cường trên nền đàn hồi Pasternak/Winkler. Tuy nhiên, phương pháp giải tích bộc lộ điểm yếu cố hữu khi không thể mô hình hóa chính xác các kết cấu vỏ có độ cong bất kỳ, có lỗ khoét hình học cục bộ, hoặc chịu chuỗi tải trọng động không tuần hoàn.
  2. Trường phái phương pháp phần tử hữu hạn: Trên thế giới, Anqi Chen, Luke A. Louca và Ahmed Y. Elghazouli [39] áp dụng mô hình Euler - Lagrange (CEL) trên phần mềm số để phân tích thùng chứa chịu nổ. Lin Jing và cộng sự [67, 68, 69] nghiên cứu thực nghiệm và mô phỏng vỏ sandwich lõi xốp nhôm. Yonghui Wang và cộng sự [97, 100] khảo sát phi tuyến vỏ thép - bê tông (SCS) và bể chứa giảm chấn bằng chất lỏng. Tại Việt Nam, Nguyễn Đức Thắng [33] đã kết hợp giải tích và PTHH cho vỏ thoải trên liên kết đàn hồi nhưng chỉ dừng lại ở sóng xung kích đơn và vỏ trơn không gân.

Cuộc tranh luận học thuật then chốt tồn tại giữa hai quan điểm:

  • Quan điểm 1: Sử dụng liên kết ngàm cứng tuyệt đối để tối đa hóa độ cứng tổng thể nhằm cản trở biến dạng tức thời.
  • Quan điểm 2 (Quan điểm của luận án): Ứng dụng liên kết đàn hồi biên (hệ lò xo phân bố) kết hợp hệ gân tăng cứng trực giao để chuyển đổi năng lượng xung kích thành biến dạng đàn hồi có kiểm soát của gối tựa, qua đó giảm đột biến ứng suất pháp $\sigma_x, \sigma_y$ và ngăn ngừa hiện tượng phá hủy giòn kết cấu.

Vị trí học thuật của luận án xác lập một bước tiến đột phá: tiên phong xây dựng mô hình phần tử vỏ cong kết hợp phần tử thanh 3D cong mô tả gân gia cường rời rạc, tích hợp hệ lò xo đàn hồi biên chịu tác dụng của chuỗi nhiều lớp sóng nổ liên tiếp, đồng thời thực hiện kiểm chứng đối sánh với hai nghiên cứu quốc tế độc lập của Steeve Chung Kim Yuen [91] và Türkmen [58], khẳng định tính ưu việt về độ chính xác và khả năng ứng dụng thực tế.

Đó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 trực tiếp lý thuyết tấm và vỏ Reissner - Mindlin sang bài toán động lực học phi tuyến của vỏ trụ thoải có gân gia cường cục bộ đặt trên hệ lò xo đàn hồi chịu tải trọng va đập đa xung kích:

                             KHUNG MÔ HÌNH LÝ THUYẾT VÀ TỔNG HỢP NỘI LỰC
 
        Tải trọng chuỗi sóng xung kích: P(t) = P1(t) + P2(t - Δt0) + ... (Sadovsky / Phản xạ)

  1. Mô hình động học và quan hệ biến dạng - chuyển vị: Trường chuyển vị tại một điểm bất kỳ trong vỏ được xác định theo mô hình 5 bậc tự do: $$u(x,y,z,t) = u_0(x,y,t) + z \theta_x(x,y,t)$$ $$v(x,y,z,t) = v_0(x,y,t) + z \theta_y(x,y,t)$$ $$w(x,y,z,t) = w_0(x,y,t)$$ Xét đến tính phi tuyến hình học theo von Kármán, véc tơ biến dạng toàn phần ${\varepsilon}$ gồm thành phần tuyến tính ${\varepsilon^L}$ và phi tuyến ${\varepsilon^N}$: $$\varepsilon_x = \frac{\partial u_0}{\partial x} - \frac{w_0}{R} + \frac{1}{2}\left(\frac{\partial w_0}{\partial x}\right)^2 + z \frac{\partial \theta_x}{\partial x}$$ $$\varepsilon_y = \frac{\partial v_0}{\partial y} + \frac{1}{2}\left(\frac{\partial w_0}{\partial y}\right)^2 + z \frac{\partial \theta_y}{\partial y}$$ $$\gamma_{xy} = \frac{\partial u_0}{\partial y} + \frac{\partial v_0}{\partial x} + \frac{\partial w_0}{\partial x}\frac{\partial w_0}{\partial y} + z \left(\frac{\partial \theta_x}{\partial y} + \frac{\partial \theta_y}{\partial x}\right)$$ $$\gamma_{xz} = \theta_x + \frac{\partial w_0}{\partial x}, \quad \gamma_{yz} = \theta_y + \frac{\partial w_0}{\partial y}$$

  2. Thiết lập phương trình cân bằng năng lượng: Áp dụng nguyên lý biến phân Hamilton $\delta \int_{t_1}^{t_2} (T_e - U_e + W_e) dt = 0$, ma trận độ cứng phi tuyến và ma trận khối lượng được tích phân chính xác trên từng phần tử vỏ và gân: $$[M_{sb}^e]{\ddot{q}^e} + [C_{sb}^e]{\dot{q}^e} + \left([K_{sb}^e] + [K_{NL}^e(q)]\right){q^e} = {f_{sb}^e(t)}$$

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

Khung phân tích của luận án tích hợp đồng thời 3 trường lý thuyết:

  • Lý thuyết cơ học phá hủy và sóng nổ khí động: Áp dụng công thức kinh nghiệm Viện sĩ M.A. Sadovsky (Xađôvxki) để xác định siêu áp sóng tới $\Delta P_\Phi$ tại khoảng cách rút gọn $\bar{R} = R / \sqrt[3]{C_{tp}}$: Luận án trích dẫn trực tiếp công thức: $$\Delta P_\Phi = 0{,}84 \frac{\sqrt[3]{C_{tp}}}{R} + 2{,}7 \frac{\sqrt[3]{C_{tp}^2}}{R^2} + 7{,}0 \frac{C_{tp}}{R^3} \quad (\text{kG/cm}^2)$$ và quy luật sóng phản xạ pháp tuyến trên bề mặt vỏ: $$\Delta P_{fx} = 2\Delta P_\Phi + \frac{6\Delta P_\Phi^2}{\Delta P_\Phi + 7{,}2} \quad (\text{kG/cm}^2)$$
  • Lý thuyết tương thích chuyển vị vỏ - gân rời rạc: Thay vì dùng phương pháp làm đều (smear out approach) vốn san phẳng độ cứng và bỏ qua hiệu ứng tập trung ứng suất tại chân gân, luận án sử dụng phương pháp mô hình hóa rời rạc (discrete stiffener model). Gân dọc và gân vành được biểu diễn bằng phần tử thanh 3D cong, bảo toàn hoàn toàn tính bất liên tục hình học và tương thích biến dạng tại mặt phân giới.
  • Lý thuyết động lực học kết cấu phi tuyến: Tích hợp ma trận liên kết đàn hồi biên $[K_{pill}^e]$ mô phỏng hệ lò xo đàn hồi đa hướng với độ cứng $k$, cho phép thiết lập hệ phương trình vi phân phi tuyến toàn cục: $$\mathbf{M} \ddot{\mathbf{q}} + \mathbf{C} \dot{\mathbf{q}} + \left[\mathbf{K}L + \mathbf{K}{NL}(\mathbf{q})\right] \mathbf{q} = \mathbf{R}(t, \Delta t_0, p_m)$$

Điều kiện biên danh định giới hạn trong miền vỏ trụ thoải ($f_0 / \min(L, W) \le 1/5$), vật liệu đẳng hướng đàn hồi tuyến tính trong suốt quá trình chịu tải va đập cực ngắn.

Phương pháp nghiên cứu tiên tiến

                                QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU TOÀN DIỆN

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (positivism paradigm), kết hợp hài hòa giữa mô hình hóa toán học giải tích số và thực nghiệm vật lý hiện trường. Thiết kế nghiên cứu đa cấp độ (multi-level design) bao gồm:

  • Cấp độ 1 (Vật liệu): Kiểm định cơ tính vật liệu chế tạo vỏ thép và độ cứng thực tế của hệ lò xo nén trên máy kéo nén vạn năng thủy lực tự động.
  • Cấp độ 2 (Phần tử & Kết cấu số): Rời rạc hóa kết cấu vỏ cong và hệ gân 3D trên môi trường số học tự lập trình.
  • Cấp độ 3 (Thực địa nổ thực tế): Đo đạc sóng phản xạ, gia tốc đỉnh và đáp ứng biến dạng động lực học tức thời dưới tác động nổ của các khối nổ TNT thực.

Quy trình nghiên cứu rigorous

  1. Giao thức thu thập dữ liệu thực nghiệm: Thí nghiệm hiện trường sử dụng lượng nổ TNT chuẩn $50\text{ g}$ kích nổ bằng kíp điện số 8, treo ở độ cao xác định theo quy chuẩn nổ không khí ($H > 15 \div 18 r_0$).
  2. Thiết bị đo lường độ nhạy cao:
    • Đầu đo áp suất sóng nổ chuyên dụng Kistler Type 6233AA1000 chịu áp lực xung kích cực lớn.
    • Đầu thu biến dạng điện trở (strain gauges) dán tại các điểm xung yếu (điểm B1 chân gân và mặt vỏ).
    • Cảm biến gia tốc áp điện trở (accelerometers G1, G2) đo đáp ứng rung động thẳng đứng.
    • Hệ thống thu nhận và phân tích tín hiệu động lực học tốc độ cao LMS (LMS Data Acquisition System).
    • Máy kiểm tra cơ tính vật liệu tĩnh - động MTS-810 Landmark xác định mô đun đàn hồi $E = 2{,}1 \times 10^{11}\text{ N/m}^2$, hệ số Poisson $\nu = 0{,}3$, khối lượng riêng $\rho = 7850\text{ kg/m}^3$.

Data và phân tích

  • Kỹ thuật số hóa và phần mềm phát triển: Luận án phát triển thành công bộ chương trình tính toán số chuyên biệt CPAS_2018 (Compute_PlateAndShell_2018) viết bằng ngôn ngữ MATLAB kết nối truyền dữ liệu trực tiếp với ANSYS Mechanical. Thuật toán giải hệ vi phân phi tuyến sử dụng tích phân từng bước Newmark - $\beta$ ($\gamma = 0{,}5; \beta = 0{,}25$) kết hợp vòng lặp điều chỉnh tiếp tuyến ma trận Newton - Raphson: $$\mathbf{K}t^{(i-1)} \Delta \mathbf{q}^{(i)} = \mathbf{R}{t+\Delta t} - \mathbf{F}_{int}^{(i-1)}$$
  • Kiểm tra độ tin cậy (Robustness checks): Kết quả tính toán từ chương trình CPAS_2018 được kiểm chứng chéo với các công trình chuẩn quốc tế. Cụ thể, khi so sánh dạng dao động riêng và gia tốc đỉnh với kết quả của Türkmen [58] và mô phỏng trên ANSYS Shell281, sai số tương đối không vượt quá $4{,}8%$. Đối với thực nghiệm hiện trường, sai lệch giữa gia tốc đỉnh tính toán và đo đạc thực tế chỉ dao động trong khoảng $5{,}2% \div 9{,}4%$, chứng minh độ hội tụ và độ tin cậy tuyệt đối của giải thuật.

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

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

                  TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CỦA NGHIÊN CỨU
  1. Hiệu ứng cộng hưởng chồng lấn tải trọng: Phát hiện phản trực giác quan trọng nhất là mức độ tàn phá của chuỗi vụ nổ không phụ thuộc tuyến tính vào tổng năng lượng thuốc nổ mà phụ thuộc mật thiết vào khoảng thời gian trễ $\Delta t_0$. Khi $\Delta t_0$ xấp xỉ chu kỳ dao động riêng cơ bản thứ nhất ($T_1$) của vỏ, chuyển vị đỉnh $W_{\max}$ và ứng suất pháp $\sigma_x$ tăng vọt từ $145%$ đến $182%$ so với tải nổ đơn có cùng tổng đương lượng TNT do hiện tượng cộng hưởng pha nén.
  2. Quy luật tỷ số hình học gân tối ưu ($h_g/b_g$): Khi tăng tỷ số chiều cao trên bề rộng gân $h_g/b_g$ từ $1{,}0$ lên $3{,}5$ (giữ nguyên diện tích tiết diện gân), độ võng cực đại tại điểm giữa vỏ ($W_A$) giảm $42{,}6%$, và ứng suất uốn cực đại giảm $31{,}2%$. Tuy nhiên, khi $h_g/b_g > 4{,}0$, hiệu quả tăng cứng bị bão hòa và xuất hiện nguy cơ mất ổn định uốn - xoắn cục bộ của bản gân.
  3. Vai trò giảm chấn của gối tựa đàn hồi: Khác với ngàm cứng làm tập trung ứng suất cắt cực đại tại mép biên ($Q_x, Q_y$), liên kết lò xo đàn hồi ($k = 10^5 \div 10^7\text{ N/m}$) biến đổi một phần động năng xung kích thành thế năng biến dạng của gối đỡ, giúp hạ thấp đỉnh ứng suất von Mises trong lòng vỏ tới $36{,}8%$.
  4. Ảnh hưởng của độ vồng ($f_0$) và lỗ khoét kỹ thuật ($a_1 \times a_1$): Vỏ có bán kính cong $R$ càng nhỏ (độ vồng $f_0$ càng lớn) thì lực màng $N_x, N_y$ phát triển càng mạnh, triệt tiêu nhanh chuyển vị đứng. Sự xuất hiện của lỗ khoét trung tâm làm gia tăng ứng suất tập trung tại góc lỗ lên $2{,}15$ lần, tạo ra vùng tập trung biến dạng nguy hiểm đòi hỏi phải có gân vành bao quanh.

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện lý thuyết phân tích động lực học kết cấu mỏng thành mỏng có gân rời rạc trong miền biến dạng lớn, cung cấp hệ phương trình tổng quát mô tả sự tương tác giữa kết cấu đàn hồi và chuỗi tải trọng va đập phức tạp.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn tích hợp giữa mô hình tính toán số PTHH cải tiến và công nghệ đo động lực học hiện trường bằng cảm biến áp điện trong điều kiện nổ thực tế, có khả năng chuyển giao cho các bài toán phân tích va chạm tàu thuyền, máy bay rơi hoặc nổ trong hầm lò.
  • Về mặt ứng dụng và chính sách: Đề xuất hướng dẫn kỹ thuật phục vụ Cục Tác chiến, Viện Thiết kế Bộ Quốc phòng trong việc thiết kế kết cấu vòm hầm bán ngầm: bắt buộc bố trí gân tăng cứng trực giao dạng chữ T hoặc chữ nhật với tỷ số $h_g/b_g \in [2{,}5; 3{,}5]$ và tích hợp gối đỡ cao su/lò xo giảm chấn thay cho chân ngàm bê tông trực tiếp.

Limitations và Future Research

  1. Giới hạn mô hình vật liệu: Luận án giả định vật liệu làm việc trong miền đàn hồi tuyến tính đẳng hướng ($E, \nu = \text{const}$). Dưới áp lực nổ siêu cao khu gần ($p_m > 50\text{ kG/cm}^2$), vật liệu thép thực tế sẽ chuyển sang trạng thái đàn - dẻo và phụ thuộc tốc độ biến dạng (viscoplasticity / Johnson - Cook model).
  2. Giới hạn mô hình liên kết đàn hồi: Hệ gối đỡ mới chỉ khảo sát dạng lò xo đàn hồi tuyến tính đơn thuần, chưa xét đến tính phi tuyến phi đối xứng khi kéo/nén và hiện tượng tiêu tán năng lượng của liên kết đàn - nhớt (viscoelastic damper).
  3. Định hướng nghiên cứu tương lai:
    • Mở rộng thuật toán cho vỏ làm bằng vật liệu cơ tính biến thiên (FGM) và vật liệu composite nhiều lớp gia cường hạt nano graphene.
    • Tích hợp mô hình vật liệu đàn dẻo phá hủy (elasto-plastic damage mechanics) có xét đến hiệu ứng tốc độ biến dạng cao ($\dot{\varepsilon}$).
    • Nghiên cứu tối ưu hóa đa mục tiêu topology hệ gân gia cường dưới tác động của trường nhiệt - cơ kết hợp do nổ nhiệt hạch.

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

  • Tác động học thuật: Luận án mở ra hướng nghiên cứu mới về cơ học va đập đa xung kích tại Việt Nam. Các kết quả đã được công bố trên 8 công trình khoa học chuyên ngành uy tín, bao gồm các bài báo trên tạp chí khoa học quốc gia và kỷ yếu hội nghị cơ học toàn quốc.
  • Tác động công nghiệp và quốc phòng: Trực tiếp chuyển giao dữ liệu số học và quy trình tính toán cho các đơn vị tư vấn thiết kế công trình quân sự, nâng cao khả năng sống còn của các công trình bảo vệ mục tiêu chiến lược quốc gia lên trên $40%$ trước nguy cơ tấn công bằng vũ khí dẫn đường chính xác công phá liên tiếp.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật/Xây dựng: Tiếp cận hệ thống mã nguồn CPAS_2018, thuật toán Newmark kết hợp Newton - Raphson và phương pháp thiết lập phần tử thanh 3D cong tương thích vỏ.
  • Kỹ sư kết cấu và chuyên gia R&D công trình phòng thủ: Nhận được các khuyến cáo định lượng cụ thể về kích thước gân, độ cứng gối đàn hồi để ứng dụng trực tiếp vào hồ sơ thiết kế công trình ngầm.
  • Cơ quan quản lý và hoạch định tiêu chuẩn kỹ thuật: Bộ Quốc phòng, Bộ Xây dựng có cơ sở khoa học để ban hành quy chuẩn tính toán công trình chịu tải trọng nổ liên tiếp.

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?

Đóng góp độc đáo nhất là việc mở rộng thành công Lý thuyết vỏ Reissner - Mindlin kết hợp phi tuyến hình học von Kármán cho bài toán chuỗi sóng xung kích liên tiếp tác động lên vỏ trụ có gân gia cường rời rạc trên gối đàn hồi. Luận án trích dẫn trực tiếp: "tính toán kết cấu với loại tải trọng động tác dụng liên tiếp (hệ sóng xung kích, hệ vật rơi liên tiếp, v.v.) là vấn đề phức tạp vì có sự xuất hiện chồng lấn tải trọng trong quá trình làm việc của hệ." Công trình đã loại bỏ hoàn toàn hạn chế của mô hình làm đều gân truyền thống, thiết lập ma trận tương thích chuyển vị giữa phần tử vỏ cong và phần tử thanh 3D cong một cách tường minh.

2. Điểm mới về phương pháp luận so với ít nhất 2 công trình quốc tế tiền nhiệm?

So với nghiên cứu của Türkmen [58] (chỉ giải tích sóng nổ đơn trên vỏ biên tựa đơn) và Anqi Chen et al. [39] (chỉ mô phỏng số vỏ trơn chứa chất lỏng không có gân), luận án đã:

  • Tự phát triển hoàn toàn công cụ phần mềm CPAS_2018 trên nền tảng tích hợp MATLAB - ANSYS Mechanical mà không phụ thuộc vào các hộp đen thương mại thuần túy.
  • Thiết kế thành công bài toán biên đàn hồi tổng quát kết hợp kiểm chứng hiện trường bằng lượng nổ TNT thực tế và hệ thu nhận dữ liệu LMS tốc độ cao.

3. Phát hiện thực nghiệm/tính toán nào gây bất ngờ nhất?

Phát hiện về vùng trễ thời gian nguy hiểm ($\Delta t_0$ cộng hưởng): khi khoảng cách thời gian giữa hai vụ nổ liên tiếp rơi vào dải $\Delta t_0 = (0{,}9 \div 1{,}1) T_1$, đáp ứng gia tốc $\ddot{W}$ và ứng suất tức thời tại chân gân tăng đột biến hơn $180%$, cao hơn đáng kể so với việc gộp chung hai khối thuốc nổ thành một vụ nổ duy nhất có cùng năng lượng ban đầu.

4. Giao thức tái lập nghiên cứu (Replication protocol) có được cung cấp đầy đủ không?

Luận án cung cấp tường minh toàn bộ quy trình tái lập bao gồm:

  • Hệ phương trình vi phân và dạng hiển ma trận phần tử vỏ cong, thanh cong 3D, lò xo biên trong chương 2 và phụ lục mã nguồn.
  • Thông số chi tiết mẫu thử nghiệm: thép kết cấu ($E = 2{,}1 \times 10^{11}\text{ N/m}^2$, $\nu = 0{,}3$, $\rho = 7850\text{ kg/m}^3$), lò xo độ cứng $k$, lượng nổ $50\text{ g}$ TNT, vị trí tọa độ gắn cảm biến Kistler Type 6233AA1000 và các đầu đo biến dạng B1, gia tốc G1, G2.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình nghị sự 10 năm tập trung vào 3 trụ cột:

  1. Phát triển mô hình giải tích - số cho bài toán tương tác nổ - nhiệt - kết cấu (Thermal-Blast FSI) có xét đến cháy nổ thứ cấp.
  2. Ứng dụng vật liệu siêu nhẹ Auxetic (vật liệu có hệ số Poisson âm) và vật liệu composite cấu trúc tổ ong vào lớp đệm bảo vệ vỏ.
  3. Ứng dụng thuật toán Machine Learning / Deep Neural Operators để dự báo tức thời thời gian thực (real-time damage prediction) trạng thái tổn thương của công trình ngầm khi bị tập kích hỏa lực.

Kết luận

Luận án của tác giả Lê Xuân Thùy là một công trình khoa học xuất sắc, mẫu mực trong lĩnh vực Cơ kỹ thuật và Động lực học công trình công sự với 6 đóng góp cốt lõi:

  1. Thiết lập thành công mô hình toán học và phương trình động lực học phi tuyến hoàn chỉnh cho vỏ trụ thoải có hệ gân gia cường trực giao trên liên kết đàn hồi.
  2. Xây dựng thuật toán và đóng gói phần mềm chuyên dụng CPAS_2018, giải quyết triệt để bài toán chuỗi sóng xung kích do các vụ nổ liên tiếp trong không khí gây ra.
  3. Làm sáng tỏ bản chất cơ học của hiện tượng chồng lấn xung lực nổ và chỉ ra dải thời gian trễ $\Delta t_0$ cộng hưởng nguy hiểm cần tránh trong thiết kế.
  4. Xác lập bộ thông số tối ưu hóa hình học gân gia cường ($h_g/b_g \in [2{,}5; 3{,}5]$) và độ cứng gối đàn hồi nhằm triệt tiêu sóng phản xạ phá hoại.
  5. Cung cấp bộ số liệu thực nghiệm hiện trường vô cùng quý giá bằng nổ TNT thật, xác nhận độ tin cậy của thuật toán với sai số dưới $9{,}4%$.
  6. Định hình chuẩn mực tính toán mới cho các công trình quốc phòng trọng yếu, đưa nghiên cứu cơ học nổ công trình tại Việt Nam tiệm cận trình độ khoa học tiên tiến trên thế giới.