Tổng quan về luận án

Trong kỹ thuật hiện đại, kết cấu vỏ làm từ vật liệu composite nhiều lớp (laminated composite) và vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) đóng vai trò then chốt trong các ngành hàng không vũ trụ, hàng hải, công trình ngầm và công nghiệp quốc phòng nhờ tỷ thế độ bền trên trọng lượng vượt trội, khả năng kháng ăn mòn, chịu nhiệt và tùy biến hướng sợi tối ưu hóa ứng suất. Tuy nhiên, việc mô hình hóa và phân tích số học chính xác đáp ứng tĩnh, dao động riêng, phi tuyến hình học và đặc biệt là động lực học phi tuyến dưới tải trọng xung kích nổ trong môi trường nước đặt ra thách thức lớn đối với cơ học tính toán đương đại. Luận án tiến sĩ kỹ thuật "Nghiên cứu đáp ứng tĩnh và động của các kết cấu vỏ composite bằng phương pháp phần tử hữu hạn trơn (SFEM) kết hợp với phần tử vỏ MITC3" của tác giả Phạm Quốc Hòa (chuyên ngành Cơ kỹ thuật, Học viện Kỹ thuật Quân sự, Hà Nội - 2019) là công trình nghiên cứu tiên phong phát triển và ứng dụng phương pháp phần tử hữu hạn trơn trên cạnh kết hợp kỹ thuật nội suy ten xơ hỗn hợp cho phần tử tam giác (Edge-based Smoothed Finite Element Method with Mixed Interpolation of Tensorial Components - ES-MITC3).

Khoảng trống nghiên cứu (research gap) mà luận án giải quyết xuất phát từ hạn chế cố hữu của phương pháp phần tử hữu hạn (FEM) truyền thống. Các phần tử vỏ bậc thấp (như tam giác 3 nút T3) thường gặp hiện tượng "khóa cắt" (shear locking) nghiêm trọng khi chiều dày kết cấu rất mỏng ($a/h \ge 100$) và hiện tượng "khóa màng" (membrane locking), đồng thời ma trận độ cứng thường quá cứng (overly stiff), dẫn đến tốc độ hội tụ chậm và đánh giá sai lệch chuyển vị. Ngược lại, các phần tử bậc cao (như 8 nút, 9 nút Lagrangian) hoặc phương pháp không lưới (Meshfree methods) đòi hỏi chi phí tính toán cực lớn, khó xử lý điều kiện biên phức tạp. Tại Việt Nam và trên thế giới, hầu hết các nghiên cứu về vỏ composite/FGM chịu tải trọng nổ hoặc sóng xung kích ngầm chủ yếu dựa trên phương pháp giải tích (Nguyễn Đình Đức, 2014; Đào Huy Bích, 2013) vốn bị giới hạn ở các dạng hình học đơn giản và biên chuẩn hóa. Luận án đã lấp đầy khoảng trống này bằng việc thiết lập một công cụ số bậc thấp linh hoạt, hội tụ nhanh và triệt tiêu hoàn toàn hiện tượng khóa cắt.

Mục tiêu nghiên cứu được cụ thể hóa thông qua 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để xây dựng mô hình phần tử tam giác 3 nút bậc thấp cho vỏ composite/FGM vừa đảm bảo khả năng thích ứng lưới phi cấu trúc, vừa loại bỏ hiện tượng khóa cắt và cải thiện độ mềm của hệ thống?
  • H1: Việc tích hợp kỹ thuật làm trơn biến dạng trên cạnh (ES-FEM) với kỹ thuật nội suy biến dạng cắt hỗn hợp MITC3 (tạo thành phần tử ES-MITC3) sẽ làm mềm ma trận độ cứng $\tilde{\mathbf{K}}$, cho phép chuyển vị và tần số dao động tự do tiệm cận nghiệm giải tích chính xác với số bậc tự do tối thiểu.
  • RQ2: Ứng xử phi tuyến hình học và hiện tượng mất ổn định uốn bẻ (snap-through, snap-back) của vỏ composite/FGM dưới tải trọng tập trung được mô tả như thế nào qua thuật toán điều khiển số?
  • H2: Sự kết hợp giữa ES-MITC3 và thuật toán chiều dài cung (Arc-length) lặp Newton-Raphson sẽ vượt qua các điểm giới hạn (limit points) và điểm chuyển tiếp kỳ dị trên đường cong tải trọng - chuyển vị một cách ổn định.
  • RQ3: Mô hình ES-MITC3 có khả năng mô phỏng chính xác đáp ứng động lực học phi tuyến của vỏ composite chịu tác dụng của sóng xung kích nổ trong môi trường nước (underwater shock wave) hay không?
  • H3: Thuật toán tích phân Newmark kết hợp vòng lặp Newton-Raphson trên nền phần tử ES-MITC3 sẽ phản ánh chân thực trường chuyển vị tức thời và phân bố ứng suất theo thời gian khi chịu tải trọng nổ ngầm biến thiên phi tuyến.
  • RQ4: Các thông số hình học (độ cong Gauss, tỷ lệ bán kính trên chiều dày $R/h$) và thông số vật liệu (góc đặt sợi composite, hệ số mũ thể tích $n$ của FGM) ảnh hưởng định lượng ra sao đến đáp ứng tĩnh và động của kết cấu?
  • H4: Tối ưu hóa cấu trúc lớp xếp sợi và kiểm soát hệ số gradient $n$ sẽ làm giảm đáng kể biên độ chuyển vị động và ứng suất đỉnh khi kết cấu làm việc trong điều kiện va đập thủy động lực học.

Khung lý thuyết chủ đạo của luận án được xây dựng dựa trên Lý thuyết biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) của Reissner-Mindlin, định luật Hooke tổng quát cho vật liệu trực giao nhiều lớp và vật liệu FGM biến thiên liên tục qua chiều dày, nguyên lý năng lượng cực tiểu Lagrange, và kỹ thuật làm trơn gradient của G.R. Liu. Phạm vi nghiên cứu bao gồm vỏ thoải hai độ cong ($R_x = R_y = R$, độ cong Gauss dương), vỏ trụ ($R_x = R, R_y = \infty$), vỏ hyperbol ($R_x = -R_y$, độ cong Gauss âm) và tấm phẳng, chịu các điều kiện biên ngàm (Clamped - C), tựa đơn (Simply supported - S), biên tự do (Free - F), dưới tác dụng của tải trọng tĩnh phân bố đều, tải trọng tập trung và tải trọng sóng xung kích nổ thuốc nổ TNT ($M = 10\text{ kg}$, cự ly nổ $S = 3\text{ m}$) trong môi trường nước.

Literature Review và Positioning

Lịch sử tính toán kết cấu tấm vỏ ghi nhận các mốc phát triển lý thuyết then chốt từ lý thuyết vỏ mỏng cổ điển Kirchhoff-Love [125] vốn bỏ qua biến dạng trượt ngang nên chỉ đúng với kết cấu rất mỏng ($a/h > 100$). Khi kết cấu có chiều dày vừa và lớn ($a/h \le 20$), lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin [79] được áp dụng rộng rãi với hệ số hiệu chỉnh cắt $\kappa = 5/6$. Để mô tả chính xác hơn phân bố ứng suất cắt dạng parabol qua chiều dày mà không cần hệ số $\kappa$, Reddy [42] đề xuất lý thuyết biến dạng cắt bậc cao (HSDT), tiếp nối bởi công thức đồng nhất Carrera (CUF) của Tornabene và cộng sự [129], lý thuyết màng Zig-zag của Oh và Cho [91], và mô hình vi phân cầu phương tổng quát (GDQ) của Viola và cộng sự [140, 141].

Tuy nhiên, khi số hóa bằng phương pháp phần tử hữu hạn (FEM), lý thuyết Reissner-Mindlin gặp hiện tượng "khóa cắt" khi $a/h \to \infty$. Để khắc phục, Bathe và Dvorkin [20] đã giới thiệu phần tử tứ giác MITC4, Lee và Bathe [67] cùng Jeon và cộng sự [51] phát triển phần tử tam giác MITC3 và MITC3+, Bletzinger và cộng sự [26] đề xuất phương pháp rời rạc sai phân trượt (Discrete Shear Gap - DSG3), trong khi Tessler và Hughes [120] đưa ra phương pháp biến dạng tự nhiên (MIN3, ANS).

Đầu thập niên 2000, G.R. Liu và cộng sự [72, 73] đã tạo bước đột phá khi đề xuất phương pháp phần tử hữu hạn trơn (Smoothed Finite Element Method - S-FEM), bao gồm làm trơn trên miền phần tử (CS-FEM), trên nút (NS-FEM) và trên cạnh (ES-FEM). Tại Việt Nam, trường phái nghiên cứu S-FEM phát triển mạnh mẽ với các công trình của Nguyễn Thời Trung và cộng sự [122, 123] (phần tử CS-DSG3), Nguyễn Xuân Hùng và cộng sự [148, 150, 151] (ES-DSG3 cho tấm FGM và áp điện), Phùng Văn Phúc và cộng sự [135, 137], Trần Vĩnh Lộc [130], Thái Hoàng Chiến [121]. Mặt khác, các nghiên cứu chuyên sâu về vỏ FGM của Đào Huy Bích và cộng sự [23, 24, 25], Nguyễn Đình Đức và cộng sự [37, 38, 39, 40], Hoàng Văn Tùng [132], Trần Quốc Quân [105], Phạm Hồng Công [32] chủ yếu khai thác các phương pháp giải tích Runge-Kutta bậc 4, Galerkin kết hợp hàm ứng suất Airy.

Trong bức tranh tổng quan đó, luận án xác lập vị trí học thuật độc bản bằng cách dung hòa ưu điểm của hai nhánh nghiên cứu: tính linh hoạt khử khóa cắt của MITC3 và hiệu ứng làm mềm hệ thống của ES-FEM. Tác giả luận án khẳng định trong phần tổng quan: "việc kết hợp phương pháp phần tử hữu hạn trơn trên cạnh (ES-FEM) và phương pháp nội suy hỗn hợp các thành phần ten xơ (MITC3) nhằm nâng cao tính chính xác và khả năng hội tụ của kết quả phân tích tĩnh và động cho vỏ composite hiện nay chưa được công bố ở bất kỳ tài liệu nào trong nước cũng như trên thế giới".

So sánh với các nghiên cứu quốc tế tiêu biểu:

  1. So với nghiên cứu của Lee và Bathe [67], vốn chỉ áp dụng MITC3 thuần túy trên phần tử tam giác cho vỏ đẳng hướng tuyến tính, mô hình ES-MITC3 của luận án tạo ra bước nhảy vọt khi mở rộng ma trận biến dạng trơn trên miền cạnh $\Omega_k$, nâng cao độ hội tụ cho cả bài toán phi tuyến hình học và động lực học xung nổ trên vật liệu bất đẳng hướng phức tạp (Laminated composite và FGM).
  2. So với nghiên cứu của Arciniega và Reddy [18], vốn sử dụng phần tử Lagrange 7 tham số bậc cao để chống khóa màng cho vỏ FGM phi tuyến (chi phí bậc tự do ma trận toàn cục rất lớn), phương pháp ES-MITC3 của luận án chỉ sử dụng phần tử tam giác 3 nút đơn giản (mỗi nút 6 bậc tự do), tiết kiệm tài nguyên bộ nhớ nhưng vẫn mô tả chính xác hiện tượng snap-through và snap-back.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp quan trọng vào cơ học vật rắn biến dạng và lý thuyết kết cấu vỏ thông qua việc mở rộng lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin vào không gian làm trơn gradient. Cụ thể:

  • Mở rộng quan hệ chuyển vị - biến dạng phi tuyến von Kármán cho kết cấu vỏ thoải hai độ cong. Định nghĩa trường chuyển vị tổng quát 6 bậc tự do tại mỗi nút $\mathbf{d}_i = [u_0, v_0, w_0, \beta_x, \beta_y, \beta_z]^T$, giải quyết triệt để vấn đề bậc tự do xoay quanh trục pháp tuyến (drilling DOF $\beta_z$) bằng cách thiết lập hệ tọa độ ảo và ma trận xoay cục bộ, loại bỏ hiện tượng suy biến ma trận độ cứng (singularity).
  • Thiết lập mô hình giải tích cấu tạo cho hai hệ vật liệu tiên tiến: Vật liệu composite nhiều lớp với ma trận độ cứng liên kết $\mathbf{A}, \mathbf{B}, \mathbf{D}, \mathbf{C}$ phụ thuộc vào góc xoay trục sợi $\theta$, và vật liệu composite có cơ tính biến thiên (FGM) cấu thành từ Gốm (Ceramic - $ZrO_2, Al_2O_3$) và Kim loại (Metal - $Ti-6Al-4V, Al$) với hàm phân bố tỷ lệ thể tích biến thiên liên tục theo chiều dày: $$V_c(z) = \left(\frac{z}{h} + \frac{1}{2}\right)^n$$ với $n \ge 0$ là hệ số mũ thể tích, $z \in [-h/2, h/2]$.
  • Chứng minh về mặt toán học và cơ học tính toán rằng kỹ thuật làm trơn trên cạnh $\Omega_k^S$ tạo ra hiệu ứng "softening effect", bù trừ một cách tự nhiên độ cứng cưỡng bức của phương pháp phần tử hữu hạn tiêu chuẩn theo nguyên lý năng lượng cực tiểu Rayleigh-Ritz, đưa cận dưới của nghiệm năng lượng hội tụ nhanh về nghiệm giải tích thực.

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết:

  1. Lý thuyết vỏ thoải FSDT: Mô tả trường biến dạng gồm biến dạng màng $\hat{\boldsymbol{\varepsilon}}_m$, biến dạng uốn $\hat{\boldsymbol{\kappa}}$, và biến dạng trượt ngang $\hat{\boldsymbol{\gamma}}$ với hệ số hiệu chỉnh cắt trích dẫn chính xác từ văn bản luận án: "ở đây $\kappa = 5/6$ là hệ số hiệu chỉnh cắt làm tăng độ chính xác trong quá trình tính toán biến dạng cắt của vỏ".
  2. Kỹ thuật nội suy ten xơ biến dạng cắt MITC3: Khử triệt để hiện tượng khóa cắt bằng cách nội suy trường biến dạng trượt ngang $\gamma_{xz}, \gamma_{yz}$ tại trung điểm các cạnh tam giác thay vì tính trực tiếp từ đạo hàm hàm dạng nút.
  3. Kỹ thuật làm trơn trên cạnh (Edge-based Smoothing Domain): Miền tính toán $\Omega$ được phân chia thành các miền trơn $\Omega_k^S$ gắn liền với từng cạnh $k$ của lưới tam giác ($k = 1, 2, \dots, N_e$). Miền trơn $\Omega_k^S$ được tạo thành bằng cách nối hai nút đầu mút của cạnh $k$ với các trọng tâm của hai phần tử tam giác kề nhau chứa cạnh đó. Diện tích miền trơn được xác định: $$A_k = \frac{1}{3} A_1 + \frac{1}{3} A_2$$ Ma trận biến dạng trơn $\tilde{\mathbf{B}}k$ trên miền $\Omega_k^S$ được tính tích phân đường dọc theo biên $\Gamma_k$: $$\tilde{\mathbf{B}}k = \frac{1}{A_k} \int{\Omega_k^S} \mathbf{B}(\mathbf{x}) d\Omega = \frac{1}{A_k} \sum{j=1}^{N_e^k} \frac{1}{3} A_j \mathbf{B}j$$ Ma trận độ cứng trơn toàn cục được lắp ráp trực tiếp từ các ma trận độ cứng miền cạnh: $$\tilde{\mathbf{K}} = \sum{k=1}^{N_e} \tilde{\mathbf{K}}k = \sum{k=1}^{N_e} \int_{\Omega_k^S} \tilde{\mathbf{B}}_k^T \mathbf{D} \tilde{\mathbf{B}}k d\Omega = \sum{k=1}^{N_e} \tilde{\mathbf{B}}_k^T \mathbf{D} \tilde{\mathbf{B}}_k A_k$$

Khung phân tích này hoạt động trong điều kiện biên hình học: Vỏ thoải có tỷ lệ giữa chiều dài cạnh nhỏ nhất và bán kính cong của vỏ thỏa mãn trích dẫn từ nguồn: "Vỏ thoải là vỏ có độ cong nhỏ, với tỉ lệ giữa chiều dài cạnh nhỏ nhất và bán kính cong của vỏ không lớn hơn 1/5 [104]".

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng lập trường nhận thức luận thực chứng (Positivism) kết hợp phương pháp mô phỏng số học (Computational Numerical Modeling) dựa trên thuật toán tất định. Mô hình thiết kế đa cấp độ (Multi-level numerical design):

  • Cấp độ 1: Rời rạc hóa cấu trúc vi mô vật liệu qua quan hệ vi phân đàn hồi Hooke cho từng lớp composite và lớp FGM.
  • Cấp độ 2: Rời rạc hóa không gian hình học 2D trên mặt trung hòa thành hệ thống lưới tam giác 3 nút phi cấu trúc.
  • Cấp độ 3: Chuyển đổi tọa độ vi phân cục bộ $(\hat{x}, \hat{y}, \hat{z})$ của từng phần tử về tọa độ Decartes toàn cục $(x, y, z)$ thông qua các ma trận chuyển trục $\boldsymbol{\Lambda}$.
  • Cấp độ 4: Tích hợp phi tuyến hình học và tích phân động lực học thời gian thực.

Quy trình nghiên cứu rigorous

Quy trình giải bài toán số được chuẩn hóa thành ba nhánh thuật toán độc lập:

  1. Thuật toán phân tích tĩnh và dao động riêng tuyến tính: Lắp ráp ma trận độ cứng trơn $\tilde{\mathbf{K}}$ và ma trận khối lượng trơn $\tilde{\mathbf{M}}$. Giải bài toán tĩnh $\tilde{\mathbf{K}}\mathbf{d} = \mathbf{F}$ để tìm trường chuyển vị và ứng suất; giải bài toán trị riêng $(\tilde{\mathbf{K}} - \omega^2 \tilde{\mathbf{M}})\boldsymbol{\Phi} = \mathbf{0}$ bằng thuật toán Subspace Iteration để trích xuất 6 mode dao động đầu tiên cùng các tần số không thứ nguyên tương ứng $\bar{\omega} = \omega (a^2/h)\sqrt{\rho/E_2}$.
  2. Thuật toán phân tích phi tuyến tĩnh (Geometrically Nonlinear Static): Sử dụng phương pháp điều khiển tham số chiều dài cung (Arc-length method) kết hợp vòng lặp điều chỉnh hiệu chỉnh Newton-Raphson. Thuật toán này cho phép vượt qua các điểm uốn tới hạn cực đại/cực tiểu (snap-through) và các nhánh hồi quy phức tạp (snap-back) trên đồ thị tải trọng - chuyển vị mà phương pháp tải trọng thuần túy (Load control) hoặc chuyển vị thuần túy (Displacement control) bị phân kỳ.
  3. Thuật toán động lực học phi tuyến chịu sóng xung kích nổ ngầm: Tải trọng sóng xung kích trong môi trường nước sinh ra từ vụ nổ thuốc nổ TNT được mô hình hóa theo công thức thực nghiệm Cole: $$P(t) = P_m e^{-(t - t_a)/\theta}$$ với $P_m = 52.4 \times 10^6 \left(\frac{M^{1/3}}{S}\right)^{1.13} \text{ (Pa)}$ là áp suất đỉnh, $\theta = 0.084 M^{1/3} \left(\frac{M^{1/3}}{S}\right)^{-0.23} \times 10^{-3} \text{ (s)}$ là hằng số thời gian phân rã, $M$ là khối lượng thuốc nổ TNT (kg), $S$ là khoảng cách từ tâm nổ đến tâm vỏ (m). Phương trình động lực học phi tuyến: $$\tilde{\mathbf{M}}\ddot{\mathbf{d}}{t+\Delta t} + \mathbf{F}{int}(\mathbf{d}{t+\Delta t}) = \mathbf{F}{ext}(t+\Delta t)$$ được giải chính xác bằng phương pháp tích phân ẩn từng bước Newmark ($\gamma = 1/2, \beta = 1/4$) kết hợp lặp cân bằng phi tuyến Newton-Raphson tại mỗi bước thời gian $\Delta t = 10^{-5}\text{ s}$.

Data và phân tích

Toàn bộ hệ thống thuật toán và mã nguồn được tác giả tự lập trình hoàn toàn trong môi trường phần mềm MATLAB.

  • Kích thước và mật độ lưới: Khảo sát độ hội tụ lưới với các phân bố từ $4 \times 4, 8 \times 8, 12 \times 12, 16 \times 16$ đến $20 \times 20$ phần tử tam giác (tương đương 800 phần tử, 441 nút).
  • Tham số cơ học vật liệu:
    • Composite nhiều lớp đẳng hướng trực giao: $E_1/E_2 = 25$ hoặc $40$; $G_{12}/E_2 = G_{13}/E_2 = 0.5$; $G_{23}/E_2 = 0.2$; $\nu_{12} = 0.25$.
    • Composite cơ tính biến thiên FGM (Gốm Zirconia: $E_c = 151\text{ GPa}, \rho_c = 3000\text{ kg/m}^3$; Kim loại Nhôm: $E_m = 70\text{ GPa}, \rho_m = 2707\text{ kg/m}^3, \nu = 0.3$).
  • Kiểm tra tính vững (Robustness Checks): So sánh đối chuẩn kết quả tính toán của ES-MITC3 với các lời giải giải tích chính xác (Navier, 3D Elasticity), phương pháp phần tử hữu hạn tiêu chuẩn (T3, MITC3 gốc, MITC4, FEM Q9), và các công bố quốc tế của Reddy [107], Bathe [20], Kapania và Mohan [56], Tornabene [127].

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

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

  1. Hiệu năng hội tụ vượt trội và triệt tiêu hoàn toàn hiện tượng khóa cắt: Trong bài toán phân tích tĩnh vỏ thoải hai độ cong chịu tải trọng phân bố hình sin với tỷ lệ độ dày từ dày ($a/h = 10$) đến cực mỏng ($a/h = 1000$), phần tử ES-MITC3 đạt độ chính xác chuyển vị không thứ nguyên $\bar{w} = 1000 w(a/2, b/2, 0) / (h^3 E_2 / P a^4)$ tiệm cận nghiệm giải tích 3D với sai số dưới $0.5%$ ngay ở lưới thô $12 \times 12$. Trong khi phần tử tam giác cổ điển T3 bị khóa cắt hoàn toàn ở $a/h = 1000$ (chuyển vị triệt tiêu về gần 0), ES-MITC3 phản ánh chính xác độ võng đàn hồi mà không cần kỹ thuật giảm bậc tích phân nhân tạo.
  2. Độ chính xác cao trong trích xuất tần số dao động riêng: Tần số không thứ nguyên $\bar{\omega} = \omega (a^2/h)\sqrt{\rho/E_2}$ của vỏ cầu composite đối xứng $[0^\circ/90^\circ/90^\circ/0^\circ]$ và vỏ hyperbol ngàm các cạnh tính bằng ES-MITC3 cho kết quả trùng khớp với phương pháp vi phân cầu phương GDQ của Tornabene [129] và HSDT của Ferreira [42], với thời gian chạy chương trình giảm hơn $35%$ so với phần tử tứ giác 9 nút (FEM Q9).
  3. Mô phỏng chính xác hiện tượng snap-through và snap-back phi tuyến: Khi phân tích vỏ trụ FGM (Aluminum/Zirconia) chịu tải trọng tập trung tại tâm với các chiều dày $t = 6\text{ mm}$ và $t = 12\text{ mm}$, ES-MITC3 kết hợp Arc-length đã bắt trọn điểm giới hạn tới hạn và nhánh uốn bẻ snap-back phức tạp. Kết quả chỉ ra rằng khi tăng hệ số mũ thể tích $n$ (tăng hàm lượng kim loại), độ cứng uốn của vỏ giảm làm tăng độ võng và hạ thấp tải trọng tới hạn mất ổn định tĩnh.
  4. Quy luật đáp ứng động lực học dưới tải trọng sóng xung kích nổ trong nước: Mô phỏng đáp ứng chuyển vị thời gian thực của vỏ hai độ cong composite dưới tác dụng vụ nổ $10\text{ kg}$ thuốc nổ TNT ở khoảng cách $3\text{ m}$ trong nước chứng minh:
    • Áp suất đỉnh tác động tức thời gây biến dạng uốn cực đại trong khoảng $t \in [0.002, 0.005]\text{ s}$ đầu tiên trước khi chuyển sang pha dao động suy giảm tự do.
    • Góc xếp lớp cốt sợi ảnh hưởng mang tính quyết định: Cấu trúc lớp $[45^\circ/-45^\circ]$ cho khả năng phân tán năng lượng va đập thủy động và giảm chuyển vị đỉnh tại tâm vỏ tốt hơn cấu trúc $[0^\circ/90^\circ]$ khoảng $14.2%$.
    • Bán kính cong của vỏ $R$ tỷ lệ nghịch với khả năng chịu nổ: Vỏ có độ cong lớn ($R$ nhỏ) tăng cường hiệu ứng vòm màng, làm giảm đáng kể chuyển vị cực đại so với tấm phẳng cùng kích thước và độ dày.

Implications đa chiều

  • Về mặt lý thuyết: Khẳng định tính đúng đắn và sự ưu việt của việc tích hợp không gian làm trơn trên cạnh (ES-FEM) với lý thuyết ten xơ biến dạng cắt (MITC3) trong cơ học kết cấu phi tuyến; cung cấp một hệ thống công thức giải tích vi phân trơn hóa hoàn chỉnh cho các nhà nghiên cứu cơ học tính toán.
  • Về mặt phương pháp luận: Mở ra phương pháp tiếp cận số hiệu quả cho các bài toán kết cấu phức tạp: Sử dụng phần tử tam giác bậc thấp (dễ tự động hóa chia lưới cho các hình học phức tạp) nhưng đạt độ chính xác tương đương phần tử bậc cao, khắc phục nhược điểm phụ thuộc vào cấu trúc lưới hình học chữ nhật của MITC4 hoặc GDQ.
  • Về mặt ứng dụng thực tiễn: Cung cấp bộ công cụ tính toán và các biểu đồ đáp ứng định lượng phục vụ trực tiếp cho các kỹ sư thiết kế kết cấu vỏ tàu ngầm, thân tàu thủy vỏ composite, đáy tàu chống mìn, vách ngăn áp lực trong công trình ngầm và vỏ bảo vệ thiết bị thăm dò biển sâu chịu tải va chạm nổ thủy lôi.
  • Về mặt chính sách và công nghệ quốc phòng: Hỗ trợ quy trình kiểm chuẩn an toàn cấu trúc trang thiết bị quân sự, tối ưu hóa vật liệu FGM nhằm nâng cao khả năng sống còn của tàu chiến và phương tiện quân sự thế hệ mới trước các mối đe dọa nổ ngầm.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật và phạm vi áp dụng:

  • Giới hạn về mô hình hình học: Nghiên cứu chỉ giới hạn trong phạm vi lý thuyết vỏ thoải (shallow shells) với tỷ lệ độ cong $a/R \le 1/5$. Mô hình chưa áp dụng trực tiếp cho các kết cấu vỏ cong sâu (deep shells) hoặc vỏ có hình dạng bất đối xứng không gian lớn với độ cong thay đổi đột ngột.
  • Giới hạn về mô hình vật liệu: Vật liệu composite và FGM được giả định làm việc trong miền đàn hồi tuyến tính phi tuyến hình học; chưa xét đến quá trình phá hủy dẻo kim loại, nứt gãy vi mô của gốm, hoặc hiện tượng tách lớp liên pha (delamination) giữa các lớp vật liệu khi chịu tải trọng nổ cường độ lớn.
  • Giới hạn về tương tác chất lỏng - kết cấu (Fluid-Structure Interaction - FSI): Tải trọng sóng xung kích nổ trong nước được đưa vào dưới dạng hàm áp suất tác dụng trực tiếp qua công thức Cole giải tích mà chưa giải bài toán ghép đôi tương tác hai chiều đầy đủ (Fully Coupled FSI) bao gồm hiện tượng xâm thực (cavitation) và xung lực do bọt khí nổ co dãn (gas bubble pulsation).

Các hướng nghiên cứu phát triển tiếp theo bao gồm:

  1. Phát triển phần tử vỏ trơn ES-MITC3 trên nền tảng lý thuyết vỏ tổng quát không gian 3 chiều (Continuum-based Shell Elements) để phân tích vỏ cong sâu và vỏ phức hợp nhiều nhánh.
  2. Tích hợp mô hình hư hại dẻo phá hủy (Damage Mechanics) và mô hình phần tử dính kết Cohesive Zone Model (CZM) để mô phỏng sự tách lớp của composite dưới tác dụng của sóng xung kích nổ.
  3. Mở rộng giải thuật mô phỏng tương tác chất lỏng - kết cấu FSI bằng cách kết hợp ES-MITC3 với phương pháp thủy động lực học hạt mịn (Smoothed Particle Hydrodynamics - SPH) hoặc phương pháp thể tích hữu hạn (FVM).
  4. Khảo sát đáp ứng nhiệt - cơ ghép phi tuyến (Thermo-mechanical nonlinear response) của vỏ FGM trong môi trường nhiệt độ biến thiên khắc nghiệt kết hợp tải trọng va đập tốc độ cao.

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

Luận án tạo ra tác động học thuật và công nghệ sâu rộng:

  • Tác động học thuật: Đặt nền móng cho hướng phát triển các phần tử hữu hạn trơn cải tiến tại Việt Nam; đóng góp các bài báo khoa học chất lượng cao trên các tạp chí khoa học chuyên ngành trong nước và quốc tế, tạo nguồn tài liệu tham khảo giá trị cho các nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học vật rắn và Xây dựng công trình biển.
  • Chuyển giao công nghiệp và quốc phòng: Bộ chương trình tính toán số viết trên nền tảng MATLAB có khả năng tích hợp vào các gói phần mềm phân tích kết cấu nội địa, phục vụ trực tiếp cho các viện thiết kế tàu thủy quân sự (Viện Kỹ thuật Hải quân, Viện Thiết kế Tàu quân sự), các cơ sở đóng tàu và chế tạo vật liệu hàng không vũ trụ.
  • Hiệu quả kinh tế - kỹ thuật: Việc sử dụng vật liệu FGM và composite tối ưu hóa theo kết quả tính toán của luận án giúp giảm trọng lượng thân vỏ từ $20 - 30%$ so với thép hợp kim truyền thống trong khi vẫn đảm bảo độ bền uốn và khả năng hấp thụ năng lượng sóng nổ, nâng cao tuổi thọ công trình biển.

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

  • Nghiên cứu sinh và Giảng viên Cơ kỹ thuật: Tiếp cận phương pháp luận số học tiên tiến (ES-MITC3), nắm bắt quy trình xây dựng ma trận độ cứng trơn, xử lý bậc tự do drilling và thuật toán giải bài toán phi tuyến phức tạp.
  • Các nhà khoa học tính toán (Computational Mechanics): Khai thác các phát hiện về hiệu ứng làm mềm hệ thống của ES-FEM để mở rộng sang các bài toán cơ học đa trường (nhiệt - điện - từ - cơ).
  • Kỹ sư R&D công nghiệp quốc phòng và đóng tàu: Sử dụng dữ liệu định lượng về phân bố ứng suất, chuyển vị và các mode dao động để thiết kế tối ưu các vách ngăn chịu lực, vỏ tàu ngầm mini, và panel bảo vệ chống nổ.
  • Các nhà hoạch định chính sách khoa học - công nghệ quân sự: Có cơ sở khoa học xác đáng để xây dựng tiêu chuẩn thiết kế, thử nghiệm và đánh giá độ bền của các hệ thống vũ khí, khí tài quân sự làm việc trong môi trường nước biển và chịu tải trọng va đập cực hạn.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo và có tính đột phá nhất của luận án là gì?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công mô hình phần tử vỏ phẳng tam giác ES-MITC3 bằng cách tích hợp phương pháp làm trơn biến dạng trên cạnh (ES-FEM) với kỹ thuật nội suy biến dạng trượt hỗn hợp MITC3 dựa trên Lý thuyết biến dạng cắt bậc nhất FSDT. Mô hình này giải quyết đồng thời hai nghịch lý cơ bản của cơ học tính toán: Triệt tiêu hoàn toàn hiện tượng khóa cắt ở kết cấu vỏ siêu mỏng ($a/h \ge 100$) mà không gây ra hiện tượng mất ổn định ma trận hay mode dao động giả (spurious zero-energy modes), đồng thời làm mềm ma trận độ cứng giúp kết quả hội tụ đơn điệu về phía nghiệm giải tích chính xác.

2. Sự cải tiến phương pháp luận của ES-MITC3 so với các công trình quốc tế tiền nhiệm thể hiện ở điểm nào?

So với phần tử MITC3 chuẩn của Lee và Bathe (2004) vốn bị đánh giá là vẫn còn quá cứng trên lưới tam giác thô, ES-MITC3 đã thay thế phép tích phân Gauss cục bộ trên từng phần tử bằng phép tích phân biến dạng trơn trên các miền cạnh $\Omega_k^S$. So với phần tử ES-DSG3 của Nguyễn Xuân Hùng và cộng sự (2010), ES-MITC3 thể hiện độ ổn định hình học vượt trội khi chia lưới tam giác có tỷ lệ cạnh méo mó (mesh distortion).

3. Phát hiện số học nào bất ngờ và có giá trị khoa học cao nhất trong luận án?

Phát hiện bất ngờ nhất là khả năng bắt trọn điểm kỳ dị uốn bẻ kép (snap-through và snap-back) trên đường cong tải trọng - chuyển vị của vỏ trụ FGM dưới tải trọng tập trung bằng giải thuật Arc-length kết hợp ES-MITC3. Phần tử tam giác bậc thấp thông thường thường bị phân kỳ hoặc bỏ qua nhánh snap-back do ma trận tiếp tuyến trở nên suy biến cục bộ, nhưng cấu trúc ma trận làm trơn $\tilde{\mathbf{K}}$ đã duy trì được tính liên tục trơn tru của trường gradient biến dạng, cho phép thuật toán lặp hội tụ ổn định ngay tại các bước tải có độ dốc âm cực đại.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?

Có. Luận án trình bày tường minh toàn bộ quy trình toán học: Từ biểu thức ma trận vật liệu $\mathbf{A}, \mathbf{B}, \mathbf{D}, \mathbf{C}$, ma trận biến dạng $\tilde{\mathbf{B}}_k$, sơ đồ chia miền trơn $\Omega_k^S$, thuật toán chuyển đổi hệ tọa độ ảo để khử suy biến bậc tự do drilling $\beta_z$, đến sơ đồ khối chi tiết của các chương trình tính tĩnh, dao động riêng, Arc-length phi tuyến và tích phân động lực học Newmark trên MATLAB. Mọi thông số về hằng số vật liệu, kích thước hình học và bước thời gian $\Delta t$ đều được ghi nhận chi tiết, cho phép tái lập độc lập kết quả tính toán với độ sai lệch xấp xỉ $0%$.

5. Định hướng chương trình nghiên cứu 10 năm phát triển từ luận án này được phác thảo ra sao?

Chương trình 10 năm tập trung vào 4 giai đoạn chiến lược:

  • Giai đoạn 1 (Năm 1-2): Mở rộng ES-MITC3 cho bài toán vỏ composite/FGM nứt gãy và tách lớp liên hợp sử dụng kỹ thuật phần tử hữu hạn mở rộng (XFEM/S-XFEM).
  • Giai đoạn 2 (Năm 3-5): Xây dựng phần tử vỏ trơn khối 3 chiều (Smoothed Solid-Shell Elements) phân tích các cấu trúc vỏ dày phức tạp chịu tải trọng va chạm đạn đạo tốc độ cao.
  • Giai đoạn 3 (Năm 6-8): Phát triển gói phần mềm tính toán tương tác chất lỏng - kết cấu (FSI) quy mô lớn mô phỏng quá trình tạo bọt khí và sóng đập xâm thực của các vụ nổ ngầm dưới biển.
  • Giai đoạn 4 (Năm 9-10): Ứng dụng trí tuệ nhân tạo (Physics-Informed Neural Networks - PINNs) kết hợp dữ liệu từ ES-MITC3 để tối ưu hóa thời gian thực cấu trúc vi mô vật liệu FGM tự phục hồi cho các phương tiện không người lái dưới nước (UUV).

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Phạm Quốc Hòa là một công trình khoa học mẫu mực, có tính thời sự cao, giải quyết trọn vẹn và chuyên sâu bài toán cơ học kết cấu vỏ composite và FGM. Tóm lược 5 đóng góp cốt lõi của luận án:

  1. Phát triển thành công phần tử vỏ tam giác trơn trên cạnh ES-MITC3 dựa trên lý thuyết FSDT, loại bỏ triệt để hiện tượng khóa cắt và nâng cao vượt bậc tốc độ hội tụ số.
  2. Xây dựng mô hình tính toán chuyển vị và ứng suất chính xác cho các dạng vỏ thoải hai độ cong, vỏ trụ, vỏ hyperbol làm bằng composite nhiều lớp và FGM.
  3. Thiết lập thuật toán giải phi tuyến hình học Arc-length kết hợp Newton-Raphson, mô tả chính xác các hiện tượng mất ổn định uốn bẻ snap-through và snap-back của vỏ mỏng.
  4. Ứng dụng thành công phương pháp tích phân Newmark phi tuyến mô phỏng đáp ứng tức thời của vỏ composite chịu tải trọng sóng xung kích nổ trong môi trường nước theo công thức Cole.
  5. Xây dựng bộ chương trình tính toán số hoàn chỉnh trong môi trường MATLAB, được kiểm chuẩn nghiêm ngặt và đối chứng độ tin cậy cao với các công bố quốc tế uy tín.

Công trình không chỉ làm phong phú thêm kho tàng lý thuyết của trường phái phần tử hữu hạn trơn (S-FEM) mà còn mở ra những triển vọng to lớn trong việc ứng dụng vật liệu và kết cấu tiên tiến phục vụ sự nghiệp công nghiệp hóa, hiện đại hóa và củng cố tiềm lực quốc phòng - an ninh biển đảo của đất nước.