Tổng quan về luận án

Công trình nghiên cứu khoa học "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 của GS.TS Nguyễn Thái Chung tại Học viện Kỹ thuật Quân sự (2022), chuyên ngành Cơ kỹ thuật (Mã ngành: 9.01), thuộc nhóm các công trình tiên phong trong lĩnh vực cơ học công trình phòng thủ và kết cấu chịu tải trọng đặc biệt. Trong bối cảnh kỹ thuật quân sự và xây dựng hiện đại, các dạng kết cấu vỏ mỏng chịu lực phân bố dạng vòm che chắn (như hầm trú ẩn, công sự kiên cố, vòm bảo vệ khí tài tên lửa, mái kho đạn dược) thường xuyên phải đối mặt với tải trọng nổ cường độ cao. Phần lớn các công trình nghiên cứu trước đây (như công trình của Trần Anh Dũng, Nguyễn Đức Thắng, Türkmen) chủ yếu tiếp cận bài toán sóng xung kích đơn (single shock wave) hoặc khảo sát trên vỏ trơn không gân, liên kết biên lý tưởng hóa (ngàm cố định hoặc khớp tựa cứng tuyệt đối).

Khoảng trống nghiên cứu (research gap) trọng yếu mà luận án tập trung giải quyết là: Sự thiếu hụt một mô hình toán học - cơ học hoàn chỉnh và công cụ số chính xác cao để dự báo đáp ứng động phi tuyến của vỏ trụ cong có hệ gân gia cường trực giao đặt trên gối tựa đàn hồi biến dạng dưới tác động của chuỗi nhiều xung kích nổ liên tiếp (shock wave series). Hiện tượng chồng lấn tải trọng, cộng hưởng pha và tích lũy động năng giữa các vụ nổ liên tiếp đặt ra bài toán tương tác phức tạp mà các phương pháp giải tích truyền thống không thể giải quyết tường minh.

Các câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi bao gồm:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật tương tác động lực học phi tuyến của kết cấu vỏ trụ gân - gối đàn hồi biến thiên như thế nào khi chịu chuỗi sóng xung kích có khoảng trễ thời gian $\Delta t_0$?
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ ảnh hưởng tương hỗ giữa cấu hình gân gia cường ($h_g/b_g$), độ vồng $f_0$, bán kính cong $R$, lỗ khoét kỹ thuật và độ cứng liên kết lò xo $k$ đến sự phân bố lại trường ứng suất và biến dạng là gì?
  3. Giả thuyết khoa học 1 (H1): Việc kết hợp hệ gân trực giao tối ưu với liên kết đàn hồi biên sẽ kích hoạt cơ chế tiêu tán năng lượng sóng nổ, làm giảm chuyển vị đỉnh từ 25% đến 40% so với kết cấu liên kết cứng.
  4. Giả thuyết khoa học 2 (H2): Tồn tại các dải thời gian trễ tới hạn $\Delta t_{0,\text{crit}}$ giữa các xung nổ dẫn tới hiện tượng tích lũy biên độ gia tốc và ứng suất cực đại vượt ngưỡng phá hủy cục bộ dù cường độ đỉnh của từng sóng riêng lẻ nằm trong giới hạn đàn hồi.

Khung lý thuyết của luận án được xây dựng dựa trên: Lý thuyết vỏ thoải Reissner - Mindlin (tính đến biến dạng trượt bậc nhất FSDT), quan hệ biến dạng - chuyển vị phi tuyến hình học von Kármán, nguyên lý biến phân Hamilton, phương pháp phần tử hữu hạn (PTHH) sử dụng phần tử vỏ cong kết hợp phần tử thanh 3D cong, và thuật toán tích phân từng bước Newmark kết hợp lặp Newton - Raphson. Phạm vi nghiên cứu bao quát từ mô hình số hóa giải thuật đến thực nghiệm hiện trường tại bãi thử nổ với lượng nổ TNT thật, cung cấp bộ dữ liệu kiểm chứng có giá trị khoa học và độ tin cậy vượt trội.

Literature Review và Positioning

Tổng quan tài liệu chuyên ngành cho thấy sự phát triển của bài toán kết cấu chịu tải nổ trải qua ba dòng nghiên cứu chính:

Dòng thứ nhất tập trung vào lý thuyết tải trọng nổ và sóng xung kích trong không khí, khởi nguồn từ mô hình thực nghiệm kinh điển của Viện sĩ M.A. Sadovsky (1952) với công thức tính siêu áp mặt sóng tới: $$\Delta P_\Phi = 0.84\frac{\sqrt[3]{C_{tp}}}{R} + 2.7\left(\frac{\sqrt[3]{C_{tp}}}{R}\right)^2 + 7.0\left(\frac{\sqrt[3]{C_{tp}}}{R}\right)^3$$ kết hợp lý thuyết sóng phản xạ trên bề mặt vật cản của Henrych (1979) và Baker (1983): $$\Delta P_{fx} = 2\Delta P_\Phi + \frac{6\Delta P_\Phi^2}{\Delta P_\Phi + 7.2}$$ Dòng nghiên cứu này đã thiết lập chuẩn mực xác định tải trọng mặt sóng nhưng phần lớn chỉ áp dụng cho một xung nổ đơn lập.

Dòng thứ hai là các nghiên cứu giải tích và bán giải tích về tấm, vỏ chịu tải trọng động của các nhóm tác giả Đào Huy Bích, Đào Văn Dũng, Vũ Hoài Nam, Nguyễn Đình Đức (2010–2018). Các tác giả đã giải quyết thành công bài toán dao động phi tuyến của vỏ trụ FGM, vỏ composite có gân gia cường trên nền đàn hồi Pasternak chịu tải điều hòa hoặc xung lực va chạm. Tuy nhiên, hạn chế cố hữu của phương pháp giải tích (như phương pháp Galerkin, Runge-Kutta) là chỉ giải quyết được các hình học đối xứng đơn giản, không thể bao quát kết cấu có lỗ khoét kỹ thuật, cấu hình gân bất đối xứng hoặc điều kiện biên thực tế phức tạp.

Dòng thứ ba là các nghiên cứu số và thực nghiệm hiện đại trên thế giới. Điển hình, Gabriele Imbalzano, Phuong Tran, Tuan D. Ngo (2016, 2017) đã khảo sát tấm composite sandwich có lõi tăng trưởng (auxetic unit cell) chịu sóng xung kích nổ, chứng minh khả năng tiêu tán năng lượng nhờ hệ số Poisson âm nhưng chưa mở rộng cho kết cấu vỏ trụ cong có gân gia cường trên gối tựa biến dạng. Nghiên cứu của Lin Jing và cộng sự (2013, 2014) khảo sát thực nghiệm và mô phỏng số phá hủy của vỏ trụ sandwich lõi nhôm xốp chịu tải nổ, khẳng định vai trò của lớp vỏ ngoài trong việc hấp thu xung lượng nổ nhưng chỉ áp dụng tải trọng nổ đơn. Tại Việt Nam, tác giả Nguyễn Đức Thắng (2016) nghiên cứu vỏ trụ thoải trên gối đàn hồi chịu sóng xung kích đơn, mở ra hướng tiếp cận mới nhưng chưa giải quyết bài toán gia cường tăng cứng và chuỗi đa sóng nổ.

Vị trí của luận án được định vị tại giao điểm đột phá: Tích hợp mô hình phần tử vỏ cong phi tuyến và phần tử dầm 3D cong đồng thời trong môi trường thuật toán tự phát triển CPAS_2018 (Compute Plate And Shell 2018) liên kết ANSYS Mechanical và MATLAB, giải quyết hoàn chỉnh bài toán chuỗi hệ sóng xung kích nổ liên tiếp tương tác với kết cấu có gối tựa đàn hồi và kiểm chứng độc lập bằng thực nghiệm nổ TNT quy chuẩn.

Đó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 khung lý thuyết động lực học công trình thông qua việc xây dựng hệ phương trình vi phân mô tả dao động phi tuyến hoàn chỉnh của hệ vỏ trụ - gân gia cường - gối đàn hồi chịu chuỗi tải trọng động ngẫu nhiên:

$$[M]{\ddot{q}} + [C]{\dot{q}} + ([K_L] + [K_{NL}(q)]){q} = {R(t)}$$

Trong đó, $[M]$, $[C]$, $[K_L]$, $[K_{NL}]$ lần lượt là các ma trận khối lượng, ma trận cản nhớt Rayleigh ($[C] = \alpha_{ra}[M] + \beta_{ra}[K]$), ma trận độ cứng tuyến tính và ma trận độ cứng phi tuyến được suy biến từ thế năng biến dạng toàn phần $U_e$ và động năng $T_e$.

Mô hình lý thuyết đã mở rộng lý thuyết tấm vỏ Reissner - Mindlin khi giải quyết sự tương thích chuyển vị tuyệt đối tại bề mặt tiếp xúc giữa phần tử vỏ cong 2D và phần tử dầm gân cong 3D (kể đến độ lệch tâm của gân gia cường so với mặt trung hòa của vỏ). Véc tơ ngoại tải ${R(t)}$ được mô hình hóa dưới dạng chuỗi tổng quát của $N$ xung sóng nổ có bước trễ thời gian $\Delta t_i$: $$p(t) = \sum_{i=1}^{N} \Delta P_{m,i} \left(1 - \frac{t - t_i}{\tau_i}\right) \exp\left(-\frac{\alpha (t - t_i)}{\tau_i}\right) \cdot H(t - t_i)$$ với $H(t - t_i)$ là hàm bước nhảy Heaviside.

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

Khung phân tích của luận án kết hợp đồng thời ba cấu phần lý thuyết:

  1. Lý thuyết phần tử vỏ cong bậc cao: Khắc phục triệt để hiện tượng khóa cắt (shear locking) và mô tả chính xác độ cong liên tục $R$ mà không gây ra sai số hình học như khi dùng phần tử vỏ phẳng ghép (flat shell elements).
  2. Mô hình tương tác gân - vỏ lệch tâm: Kể đến toàn bộ mô men xoắn, lực cắt và lực màng truyền qua lại giữa hệ gân trực giao và mặt vỏ.
  3. Mô hình liên kết đàn hồi biên phân bố: Sử dụng ma trận độ cứng lò xo $[K_{pill}]$ dọc theo chu vi, cho phép mô phỏng sát thực tế độ mềm của công sự đặt trên nền đất đá hoặc khung đỡ biến dạng.

Biên điều kiện áp dụng (boundary conditions) được thiết lập chặt chẽ cho cả 4 cạnh: tựa đơn đàn hồi (Elastically Supported Simple Edges), ngàm đàn hồi (Elastically Clamped Edges) và tựa cứng truyền thống.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái nhận thức thực chứng kết hợp hiện thực phản biện (Critical Realism). Thiết kế nghiên cứu đa phương pháp (multi-method design) kết hợp giữa mô hình hóa giải tích - số bậc cao (Computational Mechanics) và thử nghiệm phá hủy/đo đạc sóng nổ thực địa (Field Blast Testing).

Mô hình tính toán được cụ thể hóa thông qua việc chia lưới phần tử hữu hạn tối ưu: Vỏ trụ được rời rạc hóa bằng $N_s$ phần tử vỏ cong 8 nút hoặc 4 nút có 6 bậc tự do tại mỗi nút ($u, v, w, \theta_x, \theta_y, \theta_z$), hệ gân gồm $N_{sb}$ phần tử thanh 3D cong tương thích nút hoàn toàn với vỏ, và $N_{pill}$ phần tử liên kết đàn hồi lò xo tại chu vi.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số và thực nghiệm tuân thủ các bước:

  1. Rời rạc hóa không gian & thời gian: Thiết lập ma trận phần tử trong hệ tọa độ tự nhiên, chuyển đổi sang hệ tọa độ tổng thể thông qua ma trận biến đổi tọa độ $[T_{tr}]$ chứa các cosin chỉ phương $c_{ij}$.
  2. Thuật toán giải phi tuyến: Sử dụng phương pháp Newmark-$\beta$ ($\gamma = 0.5, \beta = 0.25$) tích phân trực tiếp theo thời gian với bước thời gian siêu nhỏ $\Delta t = 10^{-5} \div 10^{-6}\text{ s}$ để đảm bảo điều kiện hội tụ Courant-Friedrichs-Lewy (CFL). Tại mỗi bước thời gian, hệ phương trình phi tuyến được giải lặp bằng phương pháp Newton - Raphson hiệu chỉnh thông qua ma trận độ cứng tiếp tuyến hiệu quả $[K_t^{(i-1)}]$.
  3. Thiết kế thực nghiệm hiện trường: Thí nghiệm được tổ chức tại bãi thử nổ chuyên dụng của Học viện Kỹ thuật Quân sự. Mẫu vỏ trụ thực tế chế tạo bằng thép hợp kim có gân gia cường hàn trực giao, đặt trên hệ khung thép cứng liên kết qua các lò xo đàn hồi đã được chuẩn hóa độ cứng $k$ trên máy nén thủy lực MTS-810 Landmark.

Data và phân tích

Thiết bị đo đạc và thu thập dữ liệu hiện trường bao gồm:

  • Hệ thống đo tải trọng nổ: Đầu đo áp suất áp điện Kistler Type 6233AA1000 chịu áp lực cao gắn tại vị trí đón sóng để ghi nhận chính xác profile xung nổ $\Delta P(t)$.
  • Hệ thống đo đáp ứng kết cấu: Cảm biến gia tốc áp điện độ nhạy cao (PCB Piezotronics) và các tấm điện trở dán đo biến dạng (Strain Gauges) gắn tại các điểm trọng yếu (điểm trung tâm A, điểm gân giao nhau B, và vị trí mép biên).
  • Hệ thống thu thập xử lý tín hiệu: Máy phân tích động học đa kênh LMS Scadas Mobile kết hợp phần mềm LMS Test.Lab với tần số lấy mẫu lên đến $204.8\text{ kHz/kênh}$, bảo toàn dạng sóng tín hiệu biến dạng và gia tốc tức thời.
  • Nguồn nổ: Khối thuốc nổ TNT chuẩn hóa (khối lượng 50g đến các mức tăng dần) kích nổ bằng kíp điện vi sai, treo tại các cao độ $H$ xác định để tạo sóng tới đơn và chuỗi sóng tới kép.

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

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

  1. Hiệu ứng cộng hưởng xung nổ tích lũy (Blast Wave Pulse Superposition): Khi khoảng cách thời gian giữa hai đợt nổ $\Delta t_0$ xấp xỉ bằng bội số nửa chu kỳ dao động riêng thứ nhất của kết cấu ($\Delta t_0 \approx (2k+1)T_1/2$), đáp ứng gia tốc $\ddot{W}A$ và ứng suất $\sigma{Ax}$ tại tâm vỏ không bị triệt tiêu mà cộng hưởng cực đại, làm biên độ dao động tăng từ 135% đến 168% so với trường hợp chỉ chịu một xung nổ đơn cùng tổng năng lượng.
  2. Quy luật hình học gân tối ưu ($h_g/b_g$ Ratio): Khảo sát biến thiên tỷ số chiều cao trên bề rộng gân $h_g/b_g$ từ 1.0 đến 5.0 chỉ ra rằng hiệu quả gia cường đạt điểm tối ưu kỹ thuật - kinh tế trong khoảng $h_g/b_g = 2.5 \div 3.0$. Khi vượt quá ngưỡng 3.0, độ cứng chống uốn cục bộ tăng nhưng xuất hiện hiện tượng tập trung ứng suất tiếp trượt tại chân gân, dễ dẫn đến phá hủy trượt liên kết giữa gân và vỏ.
  3. Cơ chế tiêu tán năng lượng của gối tựa đàn hồi: So với liên kết ngàm cứng tuyệt đối, sự hiện diện của liên kết đàn hồi có độ cứng $k$ thích hợp giúp hấp thụ động năng ban đầu của sóng xung kích, biến đổi một phần cơ năng thành thế năng biến dạng của lò xo, làm giảm giá trị đỉnh của ứng suất chính màng $\sigma_{max}$ trên vỏ tới 38.5% và giảm gia tốc đỉnh $\ddot{W}_{max}$ tới 42.1%.
  4. Ảnh hưởng phi tuyến của lỗ khoét kỹ thuật ($a_1 \times a_1$): Vỏ có lỗ khoét tại tâm làm suy giảm độ cứng tổng thể, gây dịch chuyển tần số dao động riêng cơ bản về miền tần số thấp, đồng thời tạo ra vùng tập trung ứng suất viền lỗ với hệ số tập trung ứng suất $K_\sigma > 2.8$, đòi hỏi phải bố trí gân bo viền lỗ chuyên biệt.
  5. Kiểm chứng thực nghiệm xuất sắc: Dữ liệu thực nghiệm đo đạc từ hệ thống LMS với lượng nổ TNT 50g cho thấy độ sai lệch giữa mô phỏng số bằng chương trình CPAS_2018 và dữ liệu thực tế về tần số dao động riêng, biên độ gia tốc cực đại và biến dạng tương đối dao động trong biên độ rất nhỏ (từ 6.8% đến 11.4%), khẳng định độ tin cậy tuyệt đối của mô hình phần tử vỏ cong kết hợp thanh 3D cong.

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 phi tuyến kết cấu vỏ thoải gia cường chịu tác dụng của hệ tải trọng động ngẫu nhiên, lấp đầy khoảng trống tính toán đa xung kích nổ trong cơ học công trình phòng thủ.
  • Về mặt phương pháp luận: Khẳng định tính ưu việt của việc kết hợp mô hình số phần tử vỏ cong 3D với thuật toán Newmark lặp Newton-Raphson tự lập trình, cho phép can thiệp sâu vào các tham số vi mô mà các phần mềm thương mại đóng gói không hỗ trợ.
  • Về ứng dụng thực tiễn: Cung cấp bộ công cụ tính toán và biểu đồ thiết kế tra cứu nhanh cho các kỹ sư công trình quân sự khi thiết kế hầm trú ẩn, vòm che pháo binh, hầm chỉ huy dã chiến chịu được các đợt tập kích hỏa lực liên tiếp.
  • Về mặt chính sách và tiêu chuẩn: Đóng góp luận cứ khoa học thực nghiệm để bổ sung và hiệu chỉnh Tiêu chuẩn Thiết kế Công trình Quốc phòng (TCVN/BQP) trong điều kiện vũ khí công nghệ cao tấn công chính xác nhiều lần vào cùng một mục tiêu.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Mô hình vật liệu trong thuật toán số mới dừng lại ở vật liệu đàn hồi đẳng hướng tuyến tính vi phân, chưa xét đến tính đàn - dẻo dẻo quánh (elasto-viscoplasticity) hoặc mô hình phá hủy rách vỏ nứt vỡ dưới tốc độ biến dạng cực cao (strain rate effects Johnson-Cook).
  2. Chưa xét đến trường nhiệt độ nổ đi kèm trong không gian vụ nổ gần (tương tác cơ - nhiệt đồng thời).
  3. Thí nghiệm hiện trường bị khống chế về quy mô lượng nổ TNT (50g) do điều kiện an toàn bãi thử, chưa thể kích nổ các khối nổ quy mô hàng trăm kilogam để đánh giá biến dạng dư lớn.

Chương trình nghiên cứu 10 năm tiếp theo đề xuất:

  • Mở rộng thuật toán CPAS cho vật liệu FGM (Functionally Graded Materials) và composite đa lớp gia cường sợi carbon dệt nano.
  • Phát triển mô hình tương tác chất lưu - kết cấu (Fluid-Structure Interaction - FSI) 3D mô phỏng sự biến dạng của môi trường khí nổ tương tác với vỏ mềm.
  • Nghiên cứu ứng dụng các cơ cấu giảm chấn thông minh (như vật liệu nhớt từ Magnetorheological hoặc siêu hợp kim nhớ hình SMA) vào hệ gối tựa biên để điều khiển thích nghi đáp ứng động lực học của vỏ.

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

Nghiên cứu tạo ra bước nhảy vọt trong tính toán an toàn công trình chịu nổ:

  • Tác động học thuật: Dự kiến đóng góp các công bố quốc tế chất lượng cao trên các tạp chí thuộc danh mục ISI/Scopus (Q1/Q2) như Thin-Walled Structures, International Journal of Mechanical Sciences, Journal of Fluids and Structures, đóng vai trò tài liệu tham khảo chuẩn mực cho các luận án tiến sĩ chuyên ngành Cơ kỹ thuật và Xây dựng công trình ngầm.
  • Chuyển đổi công nghiệp & quốc phòng: Trực tiếp chuyển giao quy trình tính toán cho Viện Thiết kế Bộ Quốc phòng, Viện Khoa học và Công nghệ Quân sự để tối ưu hóa kết cấu mái vòm bảo vệ khí tài hiện đại, tiết kiệm 15–20% khối lượng vật liệu chế tạo mà vẫn nâng cao khả năng sống còn của công trình trước các cuộc tập kích đường không.
  • Lợi ích xã hội: Ứng dụng trong việc thiết kế các công trình công cộng lưỡng dụng (nhà ga tàu điện ngầm, hầm giao thông xuyên núi, bể chứa hóa chất công nghiệp nguy hiểm) chống chịu các nguy cơ cháy nổ công nghiệp hoặc khủng bố.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận thuật toán số mở, khung phương trình vi phân tương tác vỏ cong - dầm 3D cong và phương pháp giải Newmark phi tuyến hoàn chỉnh.
  • Kỹ sư kết cấu công trình ngầm & quốc phòng: Nhận được các khuyến cáo định lượng về tỷ số gân $h_g/b_g$, khoảng cách gân, độ cứng gối lò xo $k$ để đưa vào bản vẽ thi công thực tế.
  • Các cơ quan quản lý & phê duyệt tiêu chuẩn xây dựng: Có căn cứ khoa học kiểm chứng bằng thực nghiệm để ban hành các quy chuẩn kỹ thuật phòng chống nổ cho hạ tầng trọng yếu quốc gia.

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ì?

Đóng góp độc đáo nhất là việc mở rộng thành công lý thuyết Reissner-Mindlin phi tuyến hình học von Kármán cho hệ kết cấu vỏ trụ có gân trực giao lệch tâm đặt trên gối tựa đàn hồi chịu tải chuỗi đa sóng nổ $p(t) = \sum p_i(t - \Delta t_i)$, thiết lập ma trận độ cứng tiếp tuyến hiệu quả $[K_t]$ giải quyết trọn vẹn hiện tượng phi tuyến hình học kết hợp tải trọng động gián đoạn.

2. Đổi mới phương pháp luận so với các nghiên cứu trước đây?

So với nghiên cứu của Türkmen (2003) dùng phương pháp giải tích Runge-Kutta chỉ áp dụng cho vỏ trơn đối xứng và nghiên cứu của Nguyễn Đức Thắng (2016) chỉ tính cho 1 xung nổ đơn trên vỏ không gân, luận án đã tự phát triển trọn vẹn phần mềm CPAS_2018 sử dụng phần tử vỏ cong kết hợp thanh 3D cong, cho phép tính toán kết cấu có hình học bất kỳ, có lỗ khoét, có gân gia cường đa hướng và chịu tải nổ liên tiếp ngẫu nhiên.

3. Phát hiện gây bất ngờ nhất có số liệu minh chứng?

Phát hiện về sự trượt pha của chuỗi sóng nổ: Một chuỗi gồm 2 sóng nổ cường độ nhỏ ($p_m = 0.5\text{ kG/cm}^2$) có khoảng trễ $\Delta t_0 = 0.005\text{ s}$ có thể gây ra ứng suất uốn đỉnh lớn hơn 28% so với một sóng nổ đơn có cường độ gấp đôi ($p_m = 1.0\text{ kG/cm}^2$). Điều này chứng minh rằng việc thiết kế công trình chỉ dựa trên áp lực đỉnh của một vụ nổ là không an toàn trước chiến thuật tấn công hỏa lực dồn dập.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) không?

Có. Luận án cung cấp toàn bộ sơ đồ thuật toán, mã nguồn chương trình tính toán số trong phần phụ lục, bảng thông số vật liệu thép ($E = 2.1 \times 10^{11}\text{ N/m}^2, \nu = 0.3, \rho = 7850\text{ kg/m}^3$), kích thước hình học chi tiết của mẫu thử nghiệm, tọa độ bố trí đầu đo áp suất Kistler 6233AA1000, cảm biến gia tốc và strain gauges, cho phép tái lập hoàn toàn quy trình số và thực nghiệm.

5. Khung chương trình nghiên cứu 10 năm tới?

Chương trình 10 năm tập trung vào 3 trục: (1) Tích hợp mô hình phá hủy cơ học phi tuyến dẻo - trượt phân lớp cho vật liệu FGM/Sandwich; (2) Ứng dụng mạng nơ-ron nhân tạo (AI/Machine Learning) dự báo nhanh đáp ứng hư hại kết cấu theo thời gian thực (Real-time Structural Health Monitoring); (3) Tối ưu hóa đa mục tiêu hình học gân và gối đỡ hấp thụ năng lượng bằng thuật toán di truyền di động.

Kết luận

  1. Xây dựng thành công hệ phương trình vi phân và thuật toán phần tử hữu hạn phi tuyến giải quyết trọn vẹn bài toán động lực học vỏ trụ thoải có gân gia cường trên liên kết đàn hồi chịu chuỗi sóng xung kích nổ trong không khí.
  2. Phát triển và làm chủ hoàn toàn chương trình máy tính CPAS_2018 kết hợp giữa ANSYS và MATLAB, tối ưu hóa độ chính xác tính toán vượt trội so với các phần tử vỏ phẳng thông thường.
  3. Làm sáng tỏ bản chất cơ học của hiện tượng cộng hưởng tải trọng đa xung nổ, xác lập dải thông số tối ưu cho tỷ số kích thước gân ($h_g/b_g = 2.5 \div 3.0$) và vai trò triệt tiêu xung lực của gối đàn hồi.
  4. Thu thập và công bố bộ cơ sở dữ liệu thực nghiệm hiện trường quý giá về gia tốc, ứng suất và áp lực sóng nổ TNT bằng hệ thống đo hiện đại LMS và cảm biến Kistler.
  5. Mở ra ba nhánh nghiên cứu đột phá: Kết cấu vỏ thông minh tiêu tán năng lượng, phương pháp đánh giá độ tin cậy công trình phòng thủ dưới đòn tấn công chính xác, và tiêu chuẩn hóa thiết kế kết cấu chống va nổ nhiều tầng.