Tổng quan về luận án

Công trình nghiên cứu tiến sĩ của nghiên cứu sinh Phạm Minh Phúc với đề tài "Nghiên cứu ảnh hưởng của vết nứt tới ổn định và dao động của kết cấu tấm FGM", thực hiện dưới sự hướng dẫn khoa học của GS.TSKH. Nguyễn Đình Đức và TS. Đoàn Hồng Đức tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội (chuyên ngành Cơ kỹ thuật, mã số: 9520101.01, năm 2023), là một công trình khoa học chuyên sâu giải quyết bài toán cơ học phi tuyến phức tạp trong cơ học vật rắn biến dạng và cơ học phá hủy hiện đại.

Bối cảnh khoa học của đề tài bắt nguồn từ sự phát triển vượt bậc của vật liệu có cơ tính biến thiên (Functionally Graded Materials – FGM) kể từ khi được khởi xướng bởi các nhà khoa học vật liệu tại Viện Sendai (Nhật Bản) vào năm 1984. Cấu trúc lai ghép giữa pha gốm (chịu nhiệt, chống oxy hóa, mô đun đàn hồi cao) và pha kim loại (dẻo dai, cản trở nứt gãy) cho phép vật liệu FGM hoạt động bền bỉ trong môi trường khắc nghiệt như vỏ khí động học tên lửa, kết cấu hàng không vũ trụ và buồng phản ứng nhiệt hạch. Tuy nhiên, trong quá trình chế tạo và vận hành dưới tải trọng cơ - nhiệt chu kỳ, sự hình thành các khuyết tật cấu trúc vi mô và vết nứt tế vi xuyên suốt chiều dày là không thể tránh khỏi. Sự xuất hiện của vết nứt làm suy giảm độ cứng cục bộ nghiêm trọng, kích hoạt quá trình mất ổn định đàn hồi (buckling) sớm và làm thay đổi toàn bộ phổ tần số dao động tự do của kết cấu tấm.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định rõ: phần lớn các nghiên cứu trước đây (như Praveen & Reddy [11], Thai & Kim [10], Natarajan [64]) chủ yếu tập trung vào kết cấu tấm FGM lý tưởng không khuyết tật hoặc sử dụng các mô hình nứt rời rạc cổ điển (như phương pháp đóng vết nứt ảo VCCT, phần tử hữu hạn mở rộng XFEM) vốn đòi hỏi kỹ thuật theo dõi đường nứt nghiêm ngặt và tái chia lưới (remeshing) phức tạp tại vùng kỳ dị ứng suất. Luận án đã tiên phong ứng dụng lý thuyết trường pha (Phase-Field Fracture Theory) kết hợp với các lý thuyết tấm biến dạng cắt bậc cao để giải quyết đồng thời bài toán mất ổn định tới hạn và phổ dao động tự do của tấm FGM có vết nứt với chiều dày biến đổi, đặt trên nền đàn hồi Pasternak hai tham số và chịu trường nhiệt độ cao.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được thiết lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa trạng thái suy giảm độ cứng liên tục do vết nứt gây ra trong tấm FGM có chiều dày thay đổi mà không gặp phải kỳ dị toán học và không cần tái lập lưới phần tử?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật biến thiên hàm mũ thể tích (Power-law index $n$), chiều dài tương đối ($c/H$), góc nghiêng vết nứt ($\alpha$) và gradient chiều dày $h(x)$ tương tác như thế nào đến tải trọng mất ổn định tới hạn ($P_{cr}$)?
  • Câu hỏi nghiên cứu 3 (RQ3): Tác động hiệp đồng của nền đàn hồi hai tham số (Winkler $k_w$, Pasternak $k_s$) và gradient nhiệt độ ($T_c, T_m$) chi phối phổ tần số và các dạng mode dao động tự do của kết cấu tấm FGM nứt ra sao?
  • Giả thuyết khoa học 1 (H1): Việc đưa biến trạng thái trường pha $s \in [0, 1]$ vào thế năng biến dạng toàn phần cho phép mô phỏng chính xác quá trình chuyển đổi "phá hủy – mềm hóa – nguyên vẹn", loại bỏ hoàn toàn hiện tượng nghẽn cắt khi kết hợp với lý thuyết tấm bậc cao của Shi.
  • Giả thuyết khoa học 2 (H2): Nền đàn hồi Pasternak có thể bù trừ đáng kể phần độ cứng bị suy giảm do vết nứt gây ra, giúp nâng cao tải trọng tới hạn và tần số dao động riêng của tấm FGM.

Khung lý thuyết tổng quát tích hợp lý thuyết biến dạng cắt bậc nhất Mindlin (FSDT), lý thuyết biến dạng cắt bậc ba của Shi (TSDT), cơ học phá hủy Griffith thông qua biến phân trường pha biến liên tục và mô hình nền đàn hồi hai tham số Pasternak. Công trình tạo đột phá định lượng khi xác lập mức độ suy giảm tải trọng tới hạn lên tới trên 32% khi tỷ lệ vết nứt tăng từ $c/H = 0.1$ đến $0.5$, đồng thời chứng minh góc nghiêng vết nứt $\alpha = 60^\circ$ khôi phục khả năng chịu lực nén đơn trục so với vết nứt trực giao $\alpha = 0^\circ$. Phạm vi nghiên cứu bao quát các dạng tấm FGM chữ nhật đồng chất và phân lớp $P\text{-FGM}$ ($\text{Al}/\text{Al}_2\text{O}_3$, $\text{SUS304}/\text{Si}_3\text{N}_4$, $\text{Ti}-6\text{Al}-4\text{V}/\text{Aluminum oxide}$), chiều dày không đổi hoặc biến thiên tuyến tính/hàm mũ, làm việc dưới các tổ hợp điều kiện biên đối xứng và bất đối xứng.

Literature Review và Positioning

Tổng quan y văn trong lĩnh vực cơ học tính toán tấm FGM phản ánh quá trình tiến hóa từ các mô hình giải tích cổ điển đến các phương pháp số hóa cao cấp. Nghiên cứu của Yang [8], Praveen & Reddy [11], và Reddy [12] đã đặt nền móng cho việc phân tích ứng xử tĩnh và động lực học phi tuyến của tấm FGM phẳng dựa trên lý thuyết tấm biến dạng cắt bậc nhất (FSDT) và bậc cao (HSDT). Tiếp đó, Thai & Kim [10] đã tích hợp nguyên lý Hamilton với lý thuyết biến dạng cắt bậc ba để khảo sát độ ổn định của tấm FGM trên nền đàn hồi Pasternak. Đối với bài toán tấm có chiều dày thay đổi, Nerantzaki & Katsikadelis [28] đã phát triển phương pháp phương trình tương tự (AEM), trong khi Wittrick & Ellen [27] chỉ ra tính chất phi liên tục của hiện tượng mất ổn định khi độ dày thay đổi theo không gian.

Trong lĩnh vực cơ học phá hủy kết cấu, hai dòng quan điểm đối lập đã tồn tại trong nhiều thập kỷ:

  1. Trường phái tiếp cận gián đoạn cục bộ (Discrete crack models): Đại diện bởi các nghiên cứu sử dụng XFEM của Moës [19], Natarajan [64], và Liu et al. [17]. Nhóm tác giả này mô hình hóa vết nứt bằng cách bổ sung các hàm làm giàu gián đoạn Heaviside và hàm kỳ dị đỉnh vết nứt. Nhược điểm lớn của hướng tiếp cận này là sự phức tạp tột độ trong việc xử lý điều kiện biên vết nứt ba chiều, khó khăn khi xử lý các vết nứt giao nhau hoặc phân nhánh, và phụ thuộc nặng nề vào việc theo dõi hình học đường nứt.
  2. Trường phái biến phân làm mịn liên tục (Continuous diffuse interface / Phase-field models): Bắt nguồn từ lý thuyết biến phân Griffith của Francfort & Marigo [74], được phát triển tính toán bởi Borden [26, 76], Kuhn [24], và Nguyễn Đình Đức cùng cộng sự [20, 21]. Hướng tiếp cận này chuyển đổi mặt gián đoạn hình học sắc nét thành một dải khuếch tán được đặc trưng bởi tham số trật tự trường pha biến thiên trơn $s \in [0, 1]$, cho phép tích phân toàn miền giải tích mà không cần xử lý kỳ dị toán học cục bộ.

Vị trí của luận án trong hệ thống y văn được định vị rõ nét: đây là công trình tiên phong kết nối thành công lý thuyết trường pha phá hủy với lý thuyết tấm biến dạng cắt bậc ba của Shi [86] trên nền tảng phần tử hữu hạn để giải quyết bài toán giá trị riêng (eigenvalue problems) kép: xác định đồng thời tải trọng mất ổn định đàn hồi tới hạn và tần số dao động tự do của tấm FGM có chiều dày thay đổi.

So sánh trực tiếp với hai nghiên cứu thực nghiệm và số hóa quốc tế tiêu biểu khẳng định tính chuẩn xác vượt trội của luận án:

  • So sánh với nghiên cứu thực nghiệm của Seifi et al. [80]: Trên mẫu tấm vuông nứt đồng chất kích thước $240 \times 240 \times 12\text{ mm}$ ($E = 70\text{ GPa}, \nu = 0.33$), với vết nứt tâm $c/H = 0.3, \alpha = 0^\circ$, kết quả tải trọng tới hạn của thực nghiệm Seifi đạt $1531\text{ N}$, trong khi mô hình số của luận án tính toán ra $1537\text{ N}$, độ sai khác chỉ $0.39%$. Ở góc nghiêng $\alpha = 30^\circ$, giá trị thực nghiệm là $1551\text{ N}$ so với $1626\text{ N}$ của luận án (sai khác $4.83%$). Toàn bộ 9 mẫu thử nghiệm đều có sai số dưới $8.97%$, hoàn toàn nằm trong dung sai gia công hình học và khuyết tật tế vi của mẫu chế tạo thực tế.
  • So sánh với nghiên cứu của Hamid Foroughi [82] và Liu et al. [17]: Khi đối chiếu hệ số ổn định $k_{cr}$ của tấm FGM có vết nứt đặt trên nền đàn hồi Pasternak ($k_w, k_s$), sai lệch giữa thuật toán Phase-Field của luận án và phương pháp XFEM/IGA của các tác giả quốc tế duy trì ổn định dưới mức $1.2%$, chứng minh tính nhất quán cao độ và năng lực triệt tiêu hiện tượng nghẽn cắt của mô hình đề xuấ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 và phát triển các hệ lý thuyết cơ học kinh điển thông qua việc xây dựng mô hình toán học giải tích - số tích hợp:

  • Phát triển lý thuyết tấm biến dạng cắt bậc ba của Shi [86]: Khắc phục triệt để nhược điểm phụ thuộc vào hệ số hiệu chỉnh cắt $\kappa = 5/6$ trong lý thuyết Mindlin (FSDT), tự động thỏa mãn điều kiện ứng suất tiếp tuyến bằng không tại hai bề mặt tự do trên và dưới của tấm ($z = \pm h/2$), đồng thời biểu diễn chính xác trường chuyển vị phi tuyến khi chiều dày tấm biến thiên không gian $h = h(x)$.
  • Mở rộng lý thuyết trường pha phá hủy (Phase-Field Fracture Mechanics): Tích hợp tham số trật tự trường pha $s(x,y)$ vào phiếm hàm năng lượng toàn phần của kết cấu tấm composite phân lớp P-FGM, thiết lập mối quan hệ giáng cấp độ cứng đàn hồi thông qua toán tử suy giảm $(1-s)^2 + \eta$ (với hệ số ổn định số $\eta \ll 1$).
  • Mô hình hóa tương tác nền hai tham số Pasternak: Thiết lập hệ phương trình vi phân chủ đạo liên kết phản lực nền trượt chống cắt $k_s \nabla^2 w$ và phản lực đứng Winkler $k_w w$ với trường biến dạng cắt bậc cao của tấm FGM bị suy giảm liên tục do vết nứt.

Mô hình lý thuyết được lượng hóa thông qua hệ phương trình cân bằng năng lượng biến phân tổng quát:

$$\Pi(d, s) = U(d, s) + U_{crack}(s) + V_{foundation}(d) + V_{thermal}(d) - W_{ext}(d)$$

Trong đó, năng lượng bề mặt vết nứt Griffith được xấp xỉ hóa theo hàm mật độ trường pha của Borden & Francfort:

$$U_{crack}(s) = \int_{\Omega} \left[ \frac{G_c}{2\ell}(1-s)^2 + \frac{G_c \ell}{2} |\nabla s|^2 \right] h(x) , d\Omega$$

với $G_c$ là tốc độ giải phóng năng lượng tới hạn, $\ell$ là tham số độ dài đặc trưng điều hòa vết nứt, và $s=0$ biểu thị trạng thái nứt hoàn toàn, $s=1$ biểu thị trạng thái vật liệu nguyên vẹn.

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

Khung phân tích được xây dựng trên sự tích hợp đa tầng giữa: (1) Quy luật phân bố thành phần thể tích hàm mũ P-FGM qua chiều dày $V_c(z) = (z/h + 1/2)^n$; (2) Lý thuyết nhiệt đàn hồi với các thông số vật liệu $E(z,T), \alpha(z,T), K(z,T)$ phụ thuộc phi tuyến vào nhiệt độ theo công thức Touloukian; (3) Lý thuyết trường pha phá hủy giải quyết gián đoạn hình học; và (4) Phương pháp phần tử hữu hạn đẳng tham số.

Khái niệm "vùng chuyển tiếp mềm hóa" (softening zone) được định nghĩa chuẩn xác như một trạng thái hình thành vi nứt liên tục, cho phép toán tử độ cứng cục bộ của phần tử hữu hạn chuyển đổi êm thuận mà không làm suy biến định thức ma trận Jacobian. Điều kiện biên của khung phân tích bao gồm các liên kết tựa đơn (SSSS), ngàm cứng (CCCC), tự do (FFFF) và hỗn hợp (CFCF), áp dụng cho các tấm chữ nhật có tỷ lệ kích thước $L/H \in [0.5, 2.0]$, tỷ lệ chiều dày $h_L/h_0 \in [0.2, 1.0]$, và dải nhiệt độ khảo sát $T \in [300\text{ K}, 600\text{ K}]$.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị theo hệ hình thực chứng (positivist paradigm), sử dụng phương pháp luận tất định dựa trên mô hình hóa toán học chính xác và mô phỏng số học (deterministic numerical simulation), được kiểm chứng chéo nghiêm ngặt bằng dữ liệu thực nghiệm chuẩn quốc tế. Thiết kế nghiên cứu đa cấp độ (multi-level design) được triển khai qua 3 tầng phân tích:

  • Cấp độ 1 (Micro/Meso-scale): Đồng nhất hóa vi cơ học xác định cơ tính hiệu dụng của vật liệu P-FGM dọc theo trục chiều dày $z$ và tính toán hàm suy giảm năng lượng phá hủy trường pha.
  • Cấp độ 2 (Element-scale): Xây dựng phần tử tấm tam giác 6 nút đẳng tham số với 5 bậc tự do tại mỗi nút ($u_0, v_0, w_0, \theta_x, \theta_y$), tổng cộng 30 bậc tự do cho mỗi phần tử, tích hợp ma trận vi phân hàm dạng bậc cao khử triệt để hiện tượng nghẽn cắt.
  • Cấp độ 3 (Structural-scale): Lập trình giải bài toán giá trị riêng tổng quát xác định tải trọng tới hạn mất ổn định $\det([K(s)] - \lambda [K_G]) = 0$ và tần số dao động tự do $\det([K(s)] - \omega^2 [M(s)]) = 0$.

Quy trình nghiên cứu rigorous

Quy trình tính toán số được thiết lập với độ tin cậy và tính lặp lại tuyệt đối thông qua 4 bước đồng bộ:

  1. Khởi tạo và làm mịn lưới cục bộ: Sử dụng ngôn ngữ lập trình phần tử hữu hạn mã nguồn mở chuyên sâu FreeFem++-cs kết hợp các thuật toán tiền xử lý. Lưới tam giác được phân bố thích nghi (adaptive refinement), làm giàu mật độ cực cao xung quanh vùng vết nứt để đảm bảo kích thước phần tử $h_e \le \ell / 2$, đáp ứng tiêu chuẩn hội tụ của trường biến pha.
  2. Giải trường pha phi tuyến: Giải phương trình Euler-Lagrange xác định phân bố biến trường pha $s(x,y)$ trên toàn miền tấm từ trường hình học vết nứt quy định bởi hàm khoảng cách $\psi(d)$.
  3. Thiết lập ma trận hệ thống: Tích phân số Gauss đa điểm trên từng phần tử để lập ma trận độ cứng đàn hồi $[K^e(s)]$, ma trận độ cứng hình học $[K_G^e]$, ma trận phản lực nền Pasternak $[K_{foundation}^e]$, và ma trận khối lượng nhất quán $[M^e(s)]$.
  4. Kiểm tra độ hội tụ lưới: Khảo sát chi tiết sự biến thiên của tải trọng tới hạn $P_{cr}$ trên tấm vuông đồng chất ($L = H = 0.24\text{ m}, h = 0.012\text{ m}, c/H = 0.3, \alpha = 0^\circ$) qua 6 mức mật độ lưới: từ 761 phần tử ($P_{cr} = 1537.99\text{ N}$), 2973 phần tử, 4129 phần tử, 4944 phần tử, 5487 phần tử đến 5942 phần tử. Kết quả hội tụ tiệm cận hoàn hảo tại mức 4944 phần tử với thời gian tính toán tối ưu, được chọn làm cấu hình chia lưới chuẩn cho toàn bộ các phân tích tiếp theo.

Data và phân tích

Dữ liệu tính toán bao quát các tổ hợp vật liệu composite tiên tiến: Nhôm – Gốm Alumina ($\text{Al}/\text{Al}_2\text{O}_3$), Thép không gỉ – Silicon Nitride ($\text{SUS304}/\text{Si}_3\text{N}_4$) và Hợp kim Titan – Zirconia. Toàn bộ các đặc trưng vật liệu phụ thuộc nhiệt độ được tính toán chính xác theo công thức:

$$P(T) = P_0 (P_{-1}T^{-1} + 1 + P_1 T + P_2 T^2 + P_3 T^3)$$

Thuật toán giải bài toán giá trị riêng quy mô lớn sử dụng phương pháp lặp không gian con (subspace iteration) và thuật toán Lanczos tích hợp trong FreeFem++-cs. Kiểm tra độ nhạy (robustness checks) được thực hiện bằng cách thay đổi độ dày dải làm mịn trường pha $\ell \in [0.005L, 0.02L]$, biến thiên tham số phân rã năng lượng và kiểm chứng tính độc lập của nghiệm số đối với hướng chia lưới tam giác.

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

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

Công trình đã mang lại 5 phát hiện khoa học cốt lõi có giá trị thực chứng và lý thuyết sâu sắc:

  1. Quy luật suy giảm độ ổn định theo chiều dài vết nứt: Tải trọng tới hạn mất ổn định $P_{cr}$ giảm phi tuyến tính khi tăng chiều dài vết nứt tương đối $c/H$. Trên tấm P-FGM ($n=1$), khi $c/H$ tăng từ $0.1$ lên $0.5$, hệ số ổn định $k_{cr}$ sụt giảm từ $21.4%$ đến $35.8%$ tùy thuộc vào điều kiện biên. Mức độ suy giảm này diễn ra mạnh nhất ở các tấm có liên kết ngàm bốn cạnh (CCCC) do sự tập trung ứng suất uốn tại vùng liên kết biên tương tác tiêu cực với vùng mềm hóa của vết nứt trung tâm.
  2. Hiệu ứng phục hồi độ cứng theo góc nghiêng vết nứt: Khi góc nghiêng vết nứt $\alpha$ tăng từ $0^\circ$ (trực giao với phương nén đơn trục $y$) lên $30^\circ$ và $60^\circ$, tải trọng tới hạn tăng trở lại tương ứng $4.8% - 12.6%$. Nguyên nhân cơ học là do diện tích hình chiếu của vết nứt lên phương vuông góc với trục tải trọng nén bị thu hẹp, làm giảm tốc độ giải phóng năng lượng biến dạng cục bộ và tái phân bố dòng truyền lực dọc trục hữu hiệu.
  3. Ảnh hưởng phức hợp của gradient chiều dày biến thiên: Tấm FGM có chiều dày giảm dần theo quy luật tuyến tính hoặc hàm mũ ($h(x) = h_0(1 - c_x x/L)$) biểu hiện hiện tượng mất đối xứng trong các mode mất ổn định và mode dao động. Vùng có chiều dày mỏng hơn ($h_L < h_0$) trở thành miền khởi phát biến dạng uốn cục bộ sớm, làm suy giảm tần số dao động tự do cơ bản ($\omega_1$) tới $41.2%$ so với tấm có chiều dày không đổi cùng khối lượng danh định.
  4. Hiệu ứng bù trừ độ cứng của nền Pasternak: Hệ số kháng cắt của nền Pasternak ($k_s$) có vai trò vượt trội so với hệ số đứng Winkler ($k_w$) trong việc kiềm chế sự suy giảm tần số dao động của tấm nứt. Khi bổ sung nền Pasternak với tham số không thứ nguyên $K_s = k_s L^2 / D_0 = 10$, tần số dao động tự do bậc 1 của tấm FGM nứt tăng trung bình từ $35%$ đến $58%$, làm thay đổi trật tự phân kỳ của các mode dao động bậc cao.
  5. Sự suy thoái nhiệt động lực học đối với tần số riêng: Trong môi trường nhiệt độ đồng nhất hoặc gradient nhiệt ($T_c = 400\text{ K} - 600\text{ K}, T_m = 300\text{ K}$), ứng suất nén nhiệt tiền mất ổn định sinh ra bên trong tấm kết hợp với sự suy giảm mô đun đàn hồi $E(T)$ làm tần số riêng giảm đột ngột. Khi nhiệt độ tiến sát ngưỡng nhiệt độ tới hạn $T_{cr}$, tần số riêng thứ nhất triệt tiêu ($\omega_1 \to 0$), đánh dấu sự chuyển tiếp sang trạng thái mất ổn định nhiệt động (thermal buckling bifurcation).

Implications đa chiều

  • Đóng góp lý thuyết: Bổ sung cơ sở toán học chuẩn xác cho cơ học phá hủy tấm composite FGM, chứng minh tính khả thi của việc giải bài toán giá trị riêng trên miền liên tục hóa trường pha mà không cần sử dụng các toán tử gián đoạn phức tạp.
  • Đổi mới phương pháp luận: Thiết lập quy trình số hóa hoàn chỉnh trên nền tảng mã nguồn mở FreeFem++-cs, mở ra phương thức tính toán chi phí thấp nhưng đạt độ chính xác tương đương các phần mềm thương mại đắt tiền (Abaqus, ANSYS) trong việc xử lý bài toán nứt kết cấu phức tạp.
  • Ứng dụng kỹ thuật thực tiễn: Cung cấp biểu đồ thiết kế kỹ thuật (design charts) và công thức số giúp kỹ sư kết cấu hàng không, tàu thủy và năng lượng hạt nhân đánh giá chính xác hệ số an toàn dư (residual safety margin), dự báo tuổi thọ mỏi và xác định ngưỡng tải trọng cho phép khi phát hiện vết nứt trong quá trình kiểm tra không phá hủy (NDT).
  • Khuyến nghị chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học để cập nhật các tiêu chuẩn thiết kế kết cấu FGM chịu nhiệt, khuyến nghị bổ sung hệ số suy giảm độ cứng bắt buộc $\gamma_{crack} \in [0.65, 0.85]$ trong quy chuẩn tính toán an toàn công trình chịu nhiệt độ cao.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn học thuật cần tiếp tục hoàn thiện:

  1. Giả thiết vết nứt tĩnh xuyên suốt (Through-crack assumption): Mô hình hiện tại khảo sát vết nứt xuyên suốt chiều dày với biên dạng thẳng, chưa bao quát được các dạng vết nứt bề mặt không xuyên (part-through cracks), vết nứt cong phức tạp hoặc hiện tượng tách lớp (delamination) cục bộ giữa các pha vật liệu.
  2. Khuôn khổ đàn hồi tuyến tính vi phân nhỏ: Nghiên cứu giả định biến dạng nhỏ và vật liệu đàn hồi tuyến tính trước mất ổn định, chưa xét đến tính phi tuyến vật liệu (chảy dẻo kim loại ở nhiệt độ cao) và phi tuyến hình học Von Kármán hoàn chỉnh trong bài toán hậu mất ổn định (post-buckling).
  3. Mô hình trường pha tĩnh: Chưa giải bài toán liên kết động lực học hoàn chỉnh giữa sự lan truyền vết nứt theo thời gian thực ($da/dt$) với đáp ứng dao động quá độ phi tuyến.

Chương trình nghiên cứu 10 năm tiếp theo được định hình qua 5 hướng đi đột phá:

  • Hướng 1: Mở rộng mô hình Phase-Field 3D đầy đủ để mô phỏng vết nứt bề mặt và tách lớp trong tấm FGM dày và vỏ thoải composite FGM.
  • Hướng 2: Tích hợp phương pháp Isogeometric Analysis (IGA) dựa trên Non-Uniform Rational B-Splines (NURBS) kết hợp trường pha để nâng cao độ trơn của trường chuyển vị và giảm triệt để số bậc tự do tính toán.
  • Hướng 3: Nghiên cứu bài toán tối ưu hóa đa mục tiêu topology và phân bố gradient vật liệu P-FGM/S-FGM nhằm tối đa hóa tải trọng mất ổn định khi có mặt vết nứt ngẫu nhiên.
  • Hướng 4: Khảo sát đáp ứng động lực học phi tuyến và hiện tượng phân nhánh động (dynamic bifurcation) của tấm FGM nứt dưới tác dụng của tải trọng cơ học di động và sóng nổ va đập.
  • Hướng 5: Phát triển mô hình dự báo khuyết tật dựa trên học máy (Physics-Informed Neural Networks – PINNs) kết hợp dữ liệu trường pha để nhận dạng và định vị vết nứt theo thời gian thực từ phổ dao động đo được.

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

  • Tác động học thuật: Các bài báo khoa học xuất bản từ luận án trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Composite Structures, Thin-Walled Structures, Applied Mathematical Modelling) đóng góp nguồn tài liệu trích dẫn phong phú, nâng cao vị thế của cơ học tính toán Việt Nam trên trường quốc tế với ước tính hàng trăm lượt trích dẫn trong vòng 5 năm tới.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp trong ngành công nghiệp chế tạo tấm chắn nhiệt tàu vũ trụ, cánh tuabin khí hàng không, thân vỏ tên lửa hành trình và các tấm lát đáy bể phản ứng nhiệt hạch thế hệ mới, giúp giảm chi phí bảo trì và ngăn ngừa thảm họa phá hủy kết cấu đột ngột.
  • Tác động xã hội và môi trường: Nâng cao độ tin cậy và hệ số an toàn của các phương tiện giao thông vận tải tốc độ cao và công trình năng lượng hạt nhân, gián tiếp giảm thiểu rủi ro tai nạn công nghiệp và tối ưu hóa hiệu suất tiêu hao vật liệu chế tạo.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận mẫu mực về ứng dụng lý thuyết Phase-Field và lý thuyết tấm bậc cao của Shi, khai thác bộ mã nguồn FreeFem++-cs mở để mở rộng cho các bài toán vỏ, dầm và tấm nano.
  • Giáo sư và nhà nghiên cứu cấp cao: Sử dụng hệ thống dữ liệu số chuẩn hóa và các công thức giải tích trong luận án làm bộ dữ liệu đối chuẩn (benchmark data) đáng tin cậy cho các phương pháp số mới (SBFEM, XIGA, Peridynamics).
  • Kỹ sư R&D trong công nghiệp quốc phòng và hàng không: Sử dụng trực tiếp các quy trình tính toán để tối ưu hóa thiết kế kết cấu tấm chịu nhiệt, xác định kích thước vết nứt giới hạn cho phép vận hành.
  • Các cơ quan quản lý và thiết lập tiêu chuẩn kỹ thuật: Tham khảo dữ liệu thực chứng để chuẩn hóa các quy trình kiểm định an toàn kết cấu composite tiên tiến chịu tải liên kết cơ - nhiệt.

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

1. Đóng góp lý thuyết độc đáo và đột phá nhất của luận án là gì? Đóng góp đột phá nhất là việc mở rộng thành công lý thuyết biến dạng cắt bậc ba của Shi [86] kết hợp với lý thuyết biến phân trường pha phá hủy (Phase-Field Fracture Theory) trên nền tảng phần tử hữu hạn để giải bài toán mất ổn định và dao động của tấm P-FGM có chiều dày thay đổi. Mô hình đã loại bỏ hoàn toàn hệ số hiệu chỉnh cắt nhân tạo, triệt tiêu hiện tượng nghẽn cắt và mô phỏng chính xác sự suy giảm độ cứng cục bộ do vết nứt mà không xuất hiện kỳ dị ứng suất toán học.

2. Điểm mới về phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm? So với phương pháp XFEM của Moës [19] và Natarajan [64] (phải theo dõi mặt gián đoạn và thêm hàm làm giàu kỳ dị) hay phương pháp IGA của Loc et al. [68], phương pháp trường pha của luận án mô tả vết nứt qua trường biến liên tục $s \in [0, 1]$. Toàn bộ bài toán trường pha và cơ học kết cấu được giải đồng bộ trên nền tảng mã nguồn mở FreeFem++-cs với kỹ thuật chia lưới thích nghi cục bộ, đảm bảo tính ổn định số tuyệt đối và giảm thiểu đáng kể chi phí thuật toán.

3. Phát hiện bất ngờ và thú vị nhất thu được từ dữ liệu số là gì? Phát hiện bất ngờ nhất là hiện tượng "tự phục hồi độ cứng hình học" khi vết nứt bị nghiêng một góc $\alpha = 60^\circ$ dưới tải nén đơn trục dọc trục $y$. Tải trọng tới hạn mất ổn định không những không giảm sâu mà còn tăng cao hơn đáng kể so với góc nghiêng $\alpha = 0^\circ$ và $30^\circ$, do sự thay đổi hướng của trường ứng suất chính và giảm tốc độ giải phóng năng lượng biến dạng đàn hồi tại đỉnh vết nứt.

4. Luận án có cung cấp quy trình và giao thức tái lập nghiên cứu (Replication Protocol) không? Có. Toàn bộ hệ phương trình biến phân, ma trận độ cứng phần tử $[K^e]$, ma trận khối lượng $[M^e]$, thuật toán giải trường pha, các thông số vật liệu phụ thuộc nhiệt độ, quy chuẩn lưới phần tử tam giác 6 nút (4944 phần tử tối ưu), và sơ đồ khối thuật toán chi tiết đều được công bố tường minh trong Chương 2 và Chương 3, cho phép cộng đồng khoa học lập trình tái tạo chính xác $100%$ các kết quả số.

5. Chương trình nghị sự nghiên cứu 10 năm tới được vạch ra như thế nào? Chương trình tập trung vào: (1) Phát triển mô hình Phase-Field 3D phi tuyến hình học Von Kármán cho tấm FGM có gân gia cường; (2) Tích hợp cơ chế phá hủy dẻo nhiệt và lan truyền vết nứt động lực học theo thời gian; (3) Ứng dụng trí tuệ nhân tạo (Physics-Informed Neural Networks) để nhận dạng và kiểm chuẩn khuyết tật vết nứt từ phổ tần số dao động thực nghiệm.

Kết luận

  1. Xây dựng thành công mô hình toán học giải tích – số tiên tiến tích hợp lý thuyết biến dạng cắt bậc ba của Shi và lý thuyết trường pha phá hủy, giải quyết triệt để bài toán mất ổn định đàn hồi và dao động tự do của tấm FGM có khuyết tật nứt.
  2. Thiết lập quy trình số hóa tin cậy trên phần mềm FreeFem++-cs với phần tử tam giác 6 nút đẳng tham số, kiểm chứng đạt độ chính xác cao với sai số dưới $0.39% - 8.97%$ so với thực nghiệm quốc tế của Seifi et al. [80].
  3. Định lượng hóa toàn diện ảnh hưởng của các tham số then chốt: tỷ lệ chiều dài vết nứt $c/H$, góc nghiêng vết nứt $\alpha$, chỉ số hàm mũ vật liệu $n$, quy luật thay đổi chiều dày $h(x)$, độ cứng nền Pasternak ($k_w, k_s$) và gradient nhiệt độ ($T_c, T_m$) tới tải trọng tới hạn $P_{cr}$ và tần số dao động tự do $\omega$.
  4. Phát hiện quy luật phục hồi độ cứng kết cấu khi xoay nghiêng góc vết nứt và cơ chế bù trừ độ cứng vượt trội của lớp kháng cắt nền Pasternak đối với kết cấu tấm FGM bị suy giảm liên kết do nứt.
  5. Mở ra 3 hướng nghiên cứu học thuật mới đầy triển vọng trong cơ học tính toán: mô phỏng trường pha 3D cho vật liệu chức năng biến thiên đa chiều, tối ưu hóa topology kết cấu FGM kháng nứt, và chẩn đoán khuyết tật thông minh dựa trên dao động phi tuyến.
  6. Đóng góp vào kho tàng khoa học kỹ thuật một công cụ tính toán mạnh mẽ, kinh tế và chuẩn xác, phục vụ đắc lực cho công tác thiết kế, kiểm định và nâng cao tuổi thọ của các kết cấu composite tiên tiến trong các ngành công nghiệp công nghệ cao.