Luận án tiến sĩ: Extended Galerkin Meshfree cho mô hình nứt vật liệu composite

Tổng quan về luận án

Nghiên cứu tiến sĩ này tập trung vào sự phát triển và ứng dụng của các phương pháp phần tử không lưới Galerkin mở rộng (Extended Galerkin Meshfree Methods) để mô hình hóa hành vi nứt gãy trong vật liệu composite chức năng tiên tiến. Luận án giải quyết một vấn đề cốt lõi trong cơ học phá hủy: dự đoán và ngăn ngừa sự thất bại của cấu trúc do nứt, đặc biệt trong các điều kiện tải trọng động và vật liệu phức tạp. Bối cảnh khoa học của công trình này được đặt trong sự nhận thức ngày càng tăng về tầm quan trọng của phân tích nứt gãy trong độ bền cấu trúc, đặc biệt khi các vật liệu tiên tiến như vật liệu composite chỉnh hướng (orthotropic composites) và vật liệu có gradient chức năng (Functionally Graded Materials – FGMs) đang được ứng dụng rộng rãi trong các ngành kỹ thuật. Tính tiên phong của nghiên cứu nằm ở việc tích hợp và cải tiến các kỹ thuật không lưới để vượt qua những hạn chế cố hữu của Phương pháp Phần tử Hữu hạn (FEM) truyền thống trong việc mô hình hóa sự gián đoạn của vết nứt và sự phức tạp của vật liệu.

Research Gap SPECIFIC với citations từ literature: Mặc dù Phương pháp Phần tử Hữu hạn mở rộng (XFEM) đã chứng minh hiệu quả trong mô hình hóa vết nứt, phương pháp này vẫn yêu cầu theo dõi vị trí vết nứt tương đối với lưới phần tử, dẫn đến sự phức tạp khi hình học miền thay đổi hoặc đường nứt tùy ý. Hơn nữa, mặc dù các mô hình không lưới dựa trên phương pháp Galerkin không lưới (EFG) đã được phát triển cho vật liệu đẳng hướng, composite chỉnh hướng và các bài toán nứt nhiệt-cơ học [29, 37, 38], luận án đã xác định một khoảng trống đáng kể: "Nevertheless, there is still little available literatures on meshfree model based on point interpolation for fracture problems" (trang 15). Khoảng trống này nhấn mạnh nhu cầu về các phương pháp không lưới dựa trên phép nội suy điểm (Point Interpolation Method – PIM) với thuộc tính Kronecker-delta, mang lại sự đơn giản hơn trong việc áp đặt các điều kiện biên thiết yếu mà không cần các kỹ thuật bổ sung như phép nhân Lagrange hay phương pháp phạt.

Research questions và hypotheses: Luận án tiến sĩ này đặt ra một số câu hỏi nghiên cứu và giả thuyết chính:

1. RQ1: Làm thế nào để phát triển một mô hình tính toán cho phân tích nứt gãy dựa trên phương pháp không lưới (RPIM và MK) kết hợp với các hàm làm giàu để mô tả toán học các gián đoạn vết nứt, đồng thời duy trì thuộc tính Kronecker-delta?
   • H1: Các phương pháp không lưới RPIM và MK, được làm giàu bằng hàm bước Heaviside cho mặt nứt và hàm nhánh tiêu chuẩn hoặc hàm dốc tuyến tính mới cho đầu nứt, có thể mô hình hóa chính xác các bài toán nứt gãy 2D.
2. RQ2: Phương pháp không lưới được phát triển có thể mở rộng để điều tra các bài toán nứt gãy trong các loại vật liệu khác nhau (đẳng hướng, chỉnh hướng, FGM) như thế nào, và làm thế nào để tích hợp hành vi vật liệu vào mô hình số?
   • H2: Bằng cách điều chỉnh các hàm làm giàu và công thức tích phân tương tác, phương pháp X-RPIM có thể tính toán chính xác các hệ số cường độ ứng suất (SIFs) cho vật liệu đẳng hướng, chỉnh hướng và FGM.
3. RQ3: Phương pháp không lưới được phát triển có thể phân tích nứt gãy động như thế nào, nơi hiệu ứng quán tính không thể bỏ qua?
   • H3: X-RPIM, kết hợp với tích phân J-động và tích phân tương tác, có thể mô hình hóa chính xác hành vi nứt động và tính toán các Hệ số Cường độ Ứng suất Động (DSIFs) cho nhiều loại vật liệu.
4. RQ4: Việc giới thiệu một hàm làm giàu đầu nứt mới (hàm dốc tuyến tính) có thể nâng cao hiệu quả tính toán và giảm số bậc tự do bổ sung (DOFs) cần thiết như thế nào mà vẫn duy trì độ chính xác so với các hàm nhánh tiêu chuẩn?
   • H4: Hàm dốc tuyến tính, khi được tích hợp vào khung không lưới X-RPIM, sẽ giảm đáng kể số lượng DOFs bổ sung (chỉ hai DOFs cho mỗi nút đầu nứt so với tám DOFs của phương pháp tiêu chuẩn) trong khi vẫn đạt được độ chính xác tương đương trong tính toán DSIFs.
5. RQ5: Các hàm tương quan mới có thể cải thiện phương pháp Moving Kriging (MK) để loại bỏ ảnh hưởng của tham số kinh nghiệm người dùng và áp dụng hiệu quả cho các bài toán nứt gãy như thế nào?
   • H5: Việc sử dụng ba loại hàm tương quan mới (đa thức bậc bốn, đa thức bậc bốn bị cắt cụt và hàm Gaussian) trong phương pháp MK sẽ cải thiện độ mạnh mẽ và độ chính xác của phương pháp trong phân tích nứt gãy, loại bỏ sự phụ thuộc vào tham số kinh nghiệm người dùng.

Theoretical framework với tên theories cụ thể: Khung lý thuyết của luận án được xây dựng dựa trên Cơ học Phá hủy Đàn hồi Tuyến tính (LEFM), áp dụng các khái niệm về Hệ số Cường độ Ứng suất (SIFs) và tốc độ giải phóng năng lượng (energy release rate) để định lượng trạng thái ứng suất gần đầu nứt. Luận án mở rộng lý thuyết này bằng cách tích hợp các phương pháp số tiên tiến. Cụ thể, nó dựa trên nguyên lý Galerkin để chuyển đổi các phương trình vi phân riêng phần thành hệ thống đại số tuyến tính rời rạc. Các hàm làm giàu từ Phương pháp Phần tử Hữu hạn Mở rộng (XFEM) – bao gồm hàm bước Heaviside cho các mặt nứt và các hàm nhánh dựa trên lời giải tiệm cận cho đầu nứt – được điều chỉnh và tích hợp vào khung không lưới. Luận án cũng dựa trên lý thuyết về Hàm cơ sở xuyên tâm (Radial Basis Functions – RBFs) và đa thức để xây dựng các hàm hình dạng RPIM và lý thuyết Kriging chuyển động (Moving Kriging) với các hàm tương quan để xây dựng các hàm hình dạng MK. Đối với các vật liệu chỉnh hướng, công trình sử dụng phương trình đặc trưng của Lekhnitskii [36] để xác định các hằng số vật liệu.

Đóng góp đột phá với quantified impact:

1. Tích hợp tiên phong của Hàm dốc tuyến tính trong khung không lưới: Luận án này "first introduced in meshless context to model the singularity at crack tip" (trang 16). Việc tích hợp hàm dốc tuyến tính thay thế bốn hàm nhánh truyền thống, giảm đáng kể số bậc tự do bổ sung (DOFs) từ tám xuống còn hai cho mỗi nút đầu nứt (trang 25). Điều này "enhancing the efficient computation, reducing the number of extra DOFs, simplifying implementation, consequently saving the computational effort" (trang 22).
2. Cải tiến đột phá phương pháp Moving Kriging: Luận án giới thiệu ba loại hàm tương quan mới để cải thiện hàm hình dạng Moving Kriging (MK), trong đó "two of them are first presented in this thesis" (trang 16). Sự cải tiến này "to eliminate the effect of the user numerical experience parameter" (Tóm tắt, trang i), làm cho phương pháp MK mạnh mẽ và dễ áp dụng hơn.
3. Khung tổng quát cho phân tích nứt động đa vật liệu: Phát triển một khung tính toán thống nhất có khả năng phân tích nứt tĩnh và động trong các vật liệu đẳng hướng, chỉnh hướng và FGM. Điều này vượt ra ngoài các nghiên cứu trước đây thường tập trung vào một loại vật liệu hoặc chế độ tải trọng cụ thể.
4. Phương pháp không lưới không cần chia lại lưới hiệu quả: Các phương pháp được phát triển loại bỏ nhu cầu chia lại lưới phức tạp trong quá trình truyền nứt, một nhược điểm lớn của FEM truyền thống, giúp tiết kiệm "lots of computational efforts" (trang 5) và giải quyết các bài toán có sự thay đổi hình học miền liên tục.

Scope (sample size, timeframe) và significance: Phạm vi nghiên cứu bao gồm các bài toán nứt gãy tuyến tính 2D. Các mô hình được xác nhận bằng cách so sánh với các giải pháp phân tích, thực nghiệm hoặc các phương pháp số khác. Các nghiên cứu ví dụ bao gồm: tấm nứt cạnh đơn dưới tải kéo (Mode I), tấm nứt cạnh đơn dưới tải cắt đồng nhất (Mixed-mode), truyền nứt từ góc lượn và trong tấm lỗ tròn, phân tích nứt động của tấm nứt cạnh bán vô hạn, tấm hình chữ nhật có vết nứt nghiêng, nứt trong ống và cấu trúc phức tạp, cũng như nứt trong vật liệu chỉnh hướng và FGM. Luận án không cung cấp trực tiếp kích thước mẫu (sample size) cho các thí nghiệm vật liệu hay một khung thời gian cụ thể cho quá trình nghiên cứu, nhưng các ví dụ số thường sử dụng "20 x 40 nodes" hoặc "10 x 20" để đánh giá hội tụ (Hình 3.2, trang 45). Significance của nghiên cứu nằm ở việc cung cấp các công cụ tính toán tiên tiến hơn, chính xác và hiệu quả hơn để đánh giá độ an toàn và độ bền của cấu trúc, góp phần giảm thiểu các sự cố thảm khốc liên quan đến nứt gãy.

Literature Review và Positioning

Luận án này thực hiện một tổng hợp sâu rộng các luồng nghiên cứu chính trong cơ học phá hủy và các phương pháp số, bắt đầu từ các khái niệm nền tảng của Cơ học Phá hủy Đàn hồi Tuyến tính (LEFM) được phát triển bởi Griffith, Irwin và Dugdale. Các nghiên cứu ban đầu này đã thiết lập các khái niệm về Hệ số Cường độ Ứng suất (SIFs) và tốc độ giải phóng năng lượng [1], làm cơ sở cho phân tích định lượng sự khởi đầu và lan truyền của vết nứt.

Synthesis của major streams với TÊN TÁC GIẢ và NĂM cụ thể: Luận án xem xét các phương pháp số chính, đặc biệt là Phương pháp Phần tử Hữu hạn (FEM) [5], đã được giới thiệu từ những năm 1950 và vẫn phổ biến. Nó cũng xem xét Phương pháp Phần tử Hữu hạn Mở rộng (XFEM) [11], được giới thiệu như một giải pháp để mô hình hóa vết nứt mà không cần chia lại lưới phức tạp. XFEM đã được điều tra và áp dụng rộng rãi bởi các nhà nghiên cứu như Moës et al. [11] cho vật liệu đẳng hướng, Remmers et al. [12] cho vật liệu chỉnh hướng, và Areias et al. [13, 14] cho vật liệu có gradient chức năng. Các kỹ thuật level set để theo dõi vết nứt, như của Sethian [21] và Duflot [34], cũng được thảo luận. Một luồng nghiên cứu quan trọng khác là các phương pháp không lưới (meshfree methods) [22], được đề xuất để khắc phục nhược điểm của FEM liên quan đến việc chia lại lưới, đặc biệt trong các trường hợp biến dạng lớn hoặc truyền nứt. Các phương pháp không lưới chính được đề cập bao gồm Phương pháp Galerkin không lưới (EFG) [24], Phương pháp Petrov-Galerkin cục bộ không lưới (MLPG) [25], Phương pháp hạt nhân tái tạo (RKPM) [26], Phương pháp nội suy điểm xuyên tâm (RPIM) [27] và Phương pháp Kriging chuyển động (MK) [28, 29]. Luận án đặc biệt nhấn mạnh RPIM và MK do thuộc tính Kronecker-delta của chúng, cho phép áp đặt trực tiếp các điều kiện biên thiết yếu, một ưu điểm so với các phương pháp EFG, MLPG và RKPM thường yêu cầu phép nhân Lagrange hoặc phương pháp phạt.

Contradictions/debates với ít nhất 2 opposing views: Một trong những tranh luận chính trong cơ học phá hủy tính toán là việc lựa chọn giữa mô hình "vết nứt sắc nét" (sharp crack) và "vết nứt bôi nhòe" (smeared crack). Các mô hình vết nứt sắc nét, được đại diện bởi XFEM và các phương pháp không lưới làm giàu, coi vết nứt là những gián đoạn rõ ràng trong hình học, đòi hỏi các hàm làm giàu đặc biệt để nắm bắt sự nhảy của dịch chuyển và sự kỳ dị ứng suất tại đầu nứt. Ngược lại, các mô hình vết nứt bôi nhòe, như mô hình hư hỏng tăng cường gradient [8] và mô hình trường pha [9, 10], xem xét một vùng hư hỏng nơi các thông số vật liệu suy giảm dần. Mặc dù các mô hình hư hỏng thường đơn giản hơn về mặt tính toán, chúng thường bị giới hạn bởi sự phụ thuộc vào các hàm thực nghiệm để mô tả sự tiến hóa hư hỏng. Mô hình trường pha giải quyết hạn chế này bằng cách hình thành từ bảo toàn năng lượng nhưng yêu cầu lưới phần tử rất mịn và tốn kém về mặt tính toán. Luận án lập luận rằng mặc dù các mô hình vết nứt bôi nhòe có những ưu điểm, "The recently introduced phase field model for fracture is formulated from energy conservation, hence it alleviates the limitation of empirical functions. A disadvantage of phase field model is that it requires very fine element mesh and therefore it is time-consuming and computationally costly" (trang 6). Điều này đặt ra một sự đối lập rõ ràng với lợi thế về hiệu quả tính toán của các phương pháp không lưới mở rộng trong luận án.

Positioning trong literature với specific gap identified: Luận án định vị mình như một sự phát triển tiên tiến trong các phương pháp không lưới cho cơ học phá hủy. Trong khi XFEM đã thành công trong việc mô hình hóa vết nứt, nó vẫn đòi hỏi "track the position of crack relatively to the finite element mesh" và việc xác định các trường hợp phần tử bị cắt có thể "be tricky due to the complicated of the domain geometry and the arbitrary of crack path [23]" (trang 7). Các phương pháp không lưới, mặc dù cũng được làm giàu, tránh được vấn đề này vì chúng không yêu cầu lưới phần tử. Khoảng trống cụ thể mà luận án giải quyết là thiếu "little available literatures on meshfree model based on point interpolation for fracture problems" (trang 15), đặc biệt là việc áp dụng RPIM và MK cho các bài toán nứt gãy động và đa vật liệu. Luận án đặt ra mục tiêu "Develop a computational model for fracture analysis based on meshfree approach and enrichment functions for mathematical description of discontinuities. The class of radial basis functions is employed due to its advantage of possessing the Kronecker delta property" (trang 15).

How this advances field với concrete contributions: Nghiên cứu này thúc đẩy lĩnh vực bằng cách:

1. Cải thiện hiệu quả tính toán: Bằng cách giới thiệu hàm dốc tuyến tính mới, luận án giảm số lượng bậc tự do bổ sung (DOFs) cần thiết cho mỗi nút đầu nứt từ 8 xuống 2, dẫn đến "simplifying implementation, consequently saving the computational effort" (trang 22).
2. Tăng độ mạnh mẽ của phương pháp Kriging chuyển động: Việc giới thiệu hai hàm tương quan mới cho MK "to eliminate the effect of the user numerical experience parameter" (Tóm tắt, trang i), làm cho phương pháp dễ áp dụng và đáng tin cậy hơn.
3. Mở rộng khả năng ứng dụng: Cung cấp một khung thống nhất để phân tích nứt tĩnh và động trong vật liệu đẳng hướng, chỉnh hướng và FGM, giải quyết sự phức tạp của vật liệu tiên tiến trong các điều kiện tải trọng thực tế.

So sánh với ÍT NHẤT 2 international studies: Luận án liên tục so sánh kết quả của mình với các nghiên cứu quốc tế đã được công bố.

• So sánh với Phương pháp Phần tử Hữu hạn (FEM) và Phương pháp Phần tử Hữu hạn Mở rộng (XFEM): Luận án so sánh các kết quả DSIFs của mình với XFEM [18] và FEM [105, 125] trong nhiều bài toán khác nhau (Hình 5.4, trang 110; Hình 6.9, trang 142). Ví dụ, trong phân tích nứt động của tấm vật liệu đẳng hướng, kết quả của phương pháp hiện tại sử dụng hàm dốc tuyến tính cho thấy "a very good agreement among three approaches" khi so sánh với hàm nhánh tiêu chuẩn và giải pháp chính xác (Hình 4.8, trang 67).
• So sánh với Phương pháp Phần tử Biên (BEM) và Phương pháp Phần tử Biên miền thời gian (Time-domain BEM): Trong chương 4 và 6, luận án so sánh DSIFs cho tấm nứt nghiêng dưới tải kéo với các nghiên cứu của Wen và Aliabadi [63] sử dụng phương pháp không lưới bình phương nhỏ nhất di động (Moving Least Square Meshless) (Hình 4.12, trang 73) và BEM [100] (Hình 6.11, trang 144). Đối với vật liệu chỉnh hướng, các kết quả của X-RPIM hiện tại được so sánh với XFEM [78] và time-domain BEM [79], cũng như FEM (ANSYS) [79] (Hình 4.25, trang 89), cho thấy "a close agreement is obtained" (trang 89).
• So sánh với các phương pháp không lưới khác: Trong chương 5, luận án so sánh SIFs cho tấm FGM với nứt cạnh với Element Free Galerkin method [110, 113] (Hình 5.2, trang 107), cho thấy sự đồng nhất cao. Đối với FGM, SBFEM [126] được sử dụng làm phương pháp tham chiếu để xác nhận kết quả DSIFs (Hình 6.16, trang 151).

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

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

Luận án này đưa ra những đóng góp lý thuyết đáng kể, vượt ra ngoài việc chỉ đơn thuần áp dụng các phương pháp hiện có. Nó không chỉ mở rộng mà còn thách thức các giới hạn của một số lý thuyết cụ thể trong lĩnh vực cơ học phá hủy và phương pháp số.

• Extend/challenge WHICH specific theories (name theorists):

  • Mở rộng Lý thuyết Cơ học Phá hủy Đàn hồi Tuyến tính (LEFM): Trong khi LEFM truyền thống (Griffith, Irwin, Dugdale) cung cấp nền tảng cho sự hiểu biết về SIFs và tốc độ giải phóng năng lượng trong vật liệu đẳng hướng giòn, luận án mở rộng LEFM bằng cách phát triển các mô hình số có khả năng áp dụng cho vật liệu phức tạp hơn như vật liệu chỉnh hướng và vật liệu FGM, đồng thời xử lý các điều kiện tải trọng động. Các công thức cho trường ứng suất và dịch chuyển gần đầu nứt được điều chỉnh để phù hợp với tính dị hướng của vật liệu chỉnh hướng (Lekhnitskii [36], trang 12) và tính không đồng nhất của FGM (trang 14).
  • Thách thức các phương pháp làm giàu truyền thống trong XFEM và các phương pháp không lưới: Luận án thách thức việc sử dụng độc quyền các hàm nhánh bốn chiều tiêu chuẩn ("four-fold branch functions", trang 22) để làm giàu đầu nứt bằng cách giới thiệu một hàm dốc tuyến tính mới. Hàm dốc tuyến tính này, được đề xuất bởi các tác giả như Béchet et al. [39, 42, 43], được tích hợp lần đầu tiên vào khung không lưới X-RPIM trong bối cảnh phá hủy động (trang 24). Điều này thách thức quan niệm rằng nhiều hàm làm giàu là cần thiết, cho thấy hiệu quả có thể đạt được với ít bậc tự do bổ sung hơn (hai thay vì tám cho mỗi nút đầu nứt, trang 25).
  • Mở rộng Lý thuyết Moving Kriging (MK): Luận án đã mở rộng lý thuyết MK của Krige và Matheron (được đề cập trong [28, 29]) bằng cách cải tiến các hàm hình dạng MK thông qua việc áp dụng ba loại hàm tương quan mới (đa thức bậc bốn, đa thức bậc bốn bị cắt cụt và hàm Gaussian), trong đó hai loại được trình bày lần đầu tiên trong luận án này (trang 16). Việc cải tiến này nhằm "eliminate the effect of the user numerical experience parameter" (Tóm tắt, trang i), một điểm yếu tiềm năng của các phương pháp Kriging truyền thống.

• Conceptual framework với components và relationships: Khung khái niệm của luận án được xây dựng trên một hệ thống phân lớp các thành phần liên quan đến nhau:

  1. Cơ học Phá hủy: Các nguyên lý cơ bản của LEFM, SIFs (KI, KII, KIII), và J-integral.
  2. Mô hình vật liệu: Các mô hình ứng xử vật liệu đẳng hướng, chỉnh hướng (sử dụng tensor độ phù hợp và phương trình đặc trưng [36]), và FGM (với các tính chất vật liệu thay đổi liên tục theo không gian [2]).
  3. Phương pháp số không lưới:
     • RPIM: Sử dụng Hàm cơ sở xuyên tâm (RBFs) và đa thức để xây dựng hàm hình dạng có thuộc tính Kronecker-delta (trang 18).
     • MK: Sử dụng các hàm tương quan mới để tạo hàm hình dạng, cũng có thuộc tính Kronecker-delta (trang 131).
  4. Kỹ thuật làm giàu (Enrichment techniques):
     • Hàm bước Heaviside (H(f(x))): Để mô hình hóa sự gián đoạn dịch chuyển qua các mặt nứt (trang 21).
     • Hàm nhánh tiêu chuẩn (ψj(x)): Để nắm bắt sự kỳ dị ứng suất tại đầu nứt (trang 22).
     • Hàm dốc tuyến tính mới (R(x1)): Một sự thay thế hiệu quả hơn cho hàm nhánh, giảm DOFs (trang 24).
  5. Theo dõi vết nứt: Phương pháp Level Set vector [34] để cập nhật hình học vết nứt trong quá trình truyền nứt (trang 25). Các mối quan hệ giữa các thành phần này là hệ thống. Cơ học phá hủy cung cấp các thông số cần tính toán (SIFs). Mô hình vật liệu định nghĩa các phương trình cấu thành. Các phương pháp số không lưới cung cấp cơ sở để xấp xỉ các trường dịch chuyển và ứng suất. Các kỹ thuật làm giàu được tích hợp vào các hàm hình dạng không lưới để giải quyết các gián đoạn và kỳ dị do vết nứt gây ra. Cuối cùng, phương pháp Level Set cho phép theo dõi và cập nhật vết nứt một cách động.

• Theoretical model với propositions/hypotheses numbered: Mô hình lý thuyết được xây dựng trên cơ sở các nguyên lý bảo toàn động lượng của cơ học đàn hồi [47] và được rời rạc hóa thông qua phương pháp Galerkin yếu.

  • Mệnh đề 1: Sự gián đoạn dịch chuyển do vết nứt gây ra có thể được mô tả chính xác bằng cách làm giàu hàm hình dạng không lưới bằng hàm bước Heaviside.
  • Mệnh đề 2: Sự kỳ dị ứng suất tại đầu nứt có thể được nắm bắt bằng cách làm giàu hàm hình dạng không lưới bằng các hàm nhánh (tiêu chuẩn hoặc dốc tuyến tính).
  • Mệnh đề 3: Các phương pháp RPIM và MK, với thuộc tính Kronecker-delta, cho phép áp đặt trực tiếp các điều kiện biên thiết yếu, đơn giản hóa quá trình thực hiện số.
  • Giả thuyết 4 (H4): Việc sử dụng hàm dốc tuyến tính cho làm giàu đầu nứt sẽ giảm số bậc tự do bổ sung mà vẫn duy trì độ chính xác của SIFs động so với hàm nhánh tiêu chuẩn.
  • Giả thuyết 5 (H5): Các hàm tương quan mới trong MK sẽ loại bỏ sự phụ thuộc vào tham số kinh nghiệm người dùng, cải thiện độ mạnh mẽ và độ chính xác của phương pháp.

• Paradigm shift với EVIDENCE từ findings: Luận án không đề xuất một sự thay đổi mô hình (paradigm shift) hoàn toàn trong cơ học phá hủy mà là một sự tiến bộ đáng kể trong mô hình hóa số. Thay vì thay đổi các nguyên lý cơ bản của LEFM, nó cung cấp một công cụ mạnh mẽ hơn, hiệu quả hơn để áp dụng các nguyên lý đó vào các bài toán phức tạp mà các phương pháp truyền thống gặp khó khăn. Bằng chứng về điều này nằm ở khả năng "capture more realistic crack openings due to the higher-order continuity and non-local interpolation character [35]" (trang 8) của các phương pháp không lưới so với các phần tử hữu hạn bậc thấp. Việc tích hợp thành công hàm dốc tuyến tính để giảm DOFs và cải thiện MK để loại bỏ tham số kinh nghiệm người dùng cho thấy một sự thay đổi trong cách tiếp cận các vấn đề tính toán trong cơ học phá hủy, ưu tiên hiệu quả và tính mạnh mẽ mà không hy sinh độ chính xác.

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

Khung phân tích được đề xuất trong luận án mang tính độc đáo nhờ sự tích hợp đa lý thuyết và phương pháp tiếp cận mới.

• Integration của theories (name 3+ specific theories): Khung phân tích tích hợp sâu rộng các lý thuyết từ nhiều lĩnh vực:

  1. Lý thuyết Cơ học Đàn hồi Tuyến tính (Linear Elasticity Theory): Cung cấp các phương trình cân bằng và quan hệ cấu thành cơ bản cho vật liệu [47].
  2. Cơ học Phá hủy Đàn hồi Tuyến tính (LEFM): Đặt ra các tiêu chí và thông số để đánh giá hành vi nứt gãy (SIFs, J-integral).
  3. Lý thuyết nội suy không lưới (Meshfree Interpolation Theory): Cụ thể là Radial Point Interpolation Method (RPIM) và Moving Kriging (MK) với các hàm hình dạng có thuộc tính Kronecker-delta.
  4. Lý thuyết làm giàu của Phương pháp Phần tử Hữu hạn Mở rộng (XFEM enrichment theory): Áp dụng các hàm Heaviside và hàm nhánh (hoặc hàm dốc tuyến tính mới) để xử lý gián đoạn và kỳ dị.
  5. Lý thuyết Level Set Method (cụ thể là Vector Level Set Method [34]): Để theo dõi và cập nhật hình học vết nứt một cách hiệu quả trong quá trình lan truyền. Việc tích hợp này tạo ra một khung toàn diện có khả năng xử lý các bài toán nứt gãy phức tạp từ nhiều khía cạnh khác nhau.

• Novel analytical approach với justification: Phương pháp phân tích nổi bật là sự kết hợp của RPIM và MK với các kỹ thuật làm giàu tiên tiến, đặc biệt là việc giới thiệu hàm dốc tuyến tính mới và cải tiến các hàm tương quan MK.

  • Sử dụng hàm dốc tuyến tính cho đầu nứt: Phương pháp này thay thế các hàm nhánh bốn chiều tiêu chuẩn bằng một hàm dốc tuyến tính duy nhất kết hợp với hàm Heaviside (trang 24), giảm đáng kể số bậc tự do bổ sung. Điều này được chứng minh bằng việc đạt được "comparable accuracy of the DSIFs as compared with the responses derived from the standard method" (trang 25) trong khi giảm DOFs từ 8 xuống 2 cho mỗi nút đầu nứt.
  • Các hàm tương quan cải tiến cho MK: Luận án trình bày ba loại hàm tương quan mới, trong đó có hai loại chưa từng được giới thiệu trước đây (trang 16). Sự đổi mới này nhằm "to eliminate the effect of the user numerical experience parameter" (Tóm tắt, trang i), một vấn đề chung trong các phương pháp Kriging, làm cho phương pháp mạnh mẽ hơn và ít phụ thuộc vào kinh nghiệm người dùng.

• Conceptual contributions với definitions:

  • Extended Radial Point Interpolation Method (X-RPIM): Một phương pháp không lưới dựa trên phép nội suy điểm xuyên tâm, được làm giàu bằng các hàm đặc biệt để mô hình hóa gián đoạn và kỳ dị vết nứt. Nó được định nghĩa như một sự cải tiến của RPIM thông thường bằng cách thêm các thuật ngữ làm giàu tương tự XFEM.
  • Improved Extended Moving Kriging (X-MK): Một phương pháp không lưới dựa trên phép nội suy Kriging chuyển động, được làm giàu và cải tiến bằng các hàm tương quan mới để mô tả vết nứt, đặc biệt là nhằm loại bỏ các tham số phụ thuộc người dùng.
  • Linear Ramp Function for Crack Tip Enrichment: Một hàm làm giàu mới được đề xuất (thay vì hàm nhánh truyền thống) để mô hình hóa sự kỳ dị tại đầu nứt, được định nghĩa để giảm dần tuyến tính dịch chuyển nhảy về 0 tại đầu nứt, qua đó "ensur[ing] that the enrichment disappears exactly at the crack-tip, or to smoothly vanish the enrichment to zero at the crack-tip [39, 43]" (trang 24).

• Boundary conditions explicitly stated: Luận án đã trình bày rõ ràng các điều kiện biên của các phương trình đàn hồi:

  • Điều kiện biên Dirichlet (thiết yếu): u(x, t) = ū(x, t) trên Γu, nơi ū là vector dịch chuyển được định trước (trang 29).
  • Điều kiện biên Neumann (tự nhiên): σ · n = t̅ trên Γt, nơi t̅ là vector lực kéo được định trước (trang 29). Đặc biệt, "the Eq. (2.30) implies a traction free condition at the crack-faces" (trang 29). Một đặc tính quan trọng của RPIM và MK là "The RPIM shape functions possess the Kronecker’s delta function property regardless of the particular form of the RBF used. As a result, no special techniques for imposing the essential boundary conditions are required" (trang 21), điều này giúp đơn giản hóa việc áp đặt các điều kiện biên trong thực tế.

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu của luận án là định lượng, thực nghiệm tính toán, và được xây dựng trên nền tảng của phương pháp số tiên tiến.

• Research philosophy: Luận án áp dụng triết lý nghiên cứu Positivism (Chủ nghĩa thực chứng). Nó tập trung vào việc phát triển các mô hình tính toán có thể dự đoán và giải thích các hiện tượng vật lý (hành vi nứt gãy) thông qua các phép đo định lượng (SIFs, DSIFs, đường truyền nứt) và so sánh nghiêm ngặt với các giải pháp phân tích, thực nghiệm hoặc các phương pháp số đã được xác lập. Mục tiêu là xác minh "the correction of the developed approaches" (Tóm tắt, trang i) và cung cấp các công cụ khách quan để đánh giá độ bền cấu trúc.

• Mixed methods với SPECIFIC combination rationale: Luận án không sử dụng "mixed methods" theo nghĩa kết hợp định tính và định lượng trong khoa học xã hội. Thay vào đó, nó kết hợp các phương pháp số khác nhau theo cách tiếp cận đa phương pháp tính toán:

  • RPIM và MK: Hai phương pháp không lưới chính được sử dụng để xây dựng hàm hình dạng.
  • Kỹ thuật làm giàu của XFEM: Các hàm Heaviside và các hàm nhánh/hàm dốc tuyến tính được tích hợp vào khung không lưới.
  • Phương pháp Level Set vector: Để theo dõi hình học vết nứt trong quá trình lan truyền. Lý do cho sự kết hợp này là để tận dụng ưu điểm của từng phương pháp: phương pháp không lưới để tránh chia lại lưới, kỹ thuật làm giàu để xử lý gián đoạn và kỳ dị, và Level Set để quản lý hình học vết nứt một cách động. Sự kết hợp này mang lại "the great advantage to capture more realistic crack openings due to the higher-order continuity and non-local interpolation character [35]" (trang 8).

• Multi-level design với levels clearly defined: Mặc dù không phải là "multi-level design" theo nghĩa phân tích dữ liệu đa cấp, nghiên cứu này xử lý các mức độ phức tạp khác nhau trong mô hình hóa:

  • Mức 1: Vật liệu: Bắt đầu với vật liệu đẳng hướng đơn giản, sau đó mở rộng sang vật liệu chỉnh hướng phức tạp hơn và cuối cùng là vật liệu FGM với các tính chất thay đổi theo không gian.
  • Mức 2: Tải trọng: Phân tích các bài toán nứt tĩnh và nứt động (transient dynamic crack analysis), có tính đến hiệu ứng quán tính (trang 3).
  • Mức 3: Hình học và Điều kiện biên: Áp dụng cho nhiều hình học cấu trúc (tấm, ống, cấu trúc phức tạp) và điều kiện biên khác nhau (kéo, cắt, uốn ba điểm). Mỗi mức độ phức tạp bổ sung đều thách thức và xác nhận độ mạnh mẽ và khả năng áp dụng của các phương pháp được phát triển.

• Sample size và selection criteria EXACT: Luận án không đề cập đến "sample size" theo nghĩa thống kê từ dữ liệu thực nghiệm, mà là số lượng nút (nodes) được sử dụng trong các mô hình tính toán để đạt được sự hội tụ và độ chính xác. Ví dụ, trong các nghiên cứu hội tụ, "a cracked plate with 20 x 40 nodes" (Hình 3.2, trang 45) hoặc "10 x 20 and 20 x 40 enlarged by 10 times" (Hình 3.4, trang 47) được sử dụng để kiểm tra ảnh hưởng của mật độ nút. Tiêu chí lựa chọn cho các nút thường là "scattered nodes" (trang 7), được phân bố đồng đều hoặc không đều, tùy thuộc vào bài toán. Các nút được làm giàu (enriched nodes) được chọn dựa trên vị trí của chúng so với vết nứt:

  • Non-enriched nodes: Nằm ngoài vùng ảnh hưởng của vết nứt.
  • Nodes enriched by the step function (Heaviside): Các nút có miền hỗ trợ bị vết nứt cắt đôi (Wb, Hình 2.2, trang 22).
  • Nodes enriched by the branch functions (hoặc linear ramp function): Các nút có miền hỗ trợ bị vết nứt cắt đôi và chứa đầu nứt (Ws, Hình 2.2, trang 22).

Quy trình nghiên cứu rigorous

• Sampling strategy với inclusion/exclusion criteria: Trong ngữ cảnh phương pháp không lưới, "sampling strategy" đề cập đến việc phân bố các nút. Các nút được phân bố trong miền bài toán (bao gồm biên và bên trong miền) để tạo thành một biểu diễn rời rạc.

  • Inclusion criteria: Một nút được đưa vào miền hỗ trợ của một điểm quan tâm nếu khoảng cách từ nút đến điểm đó nhỏ hơn bán kính miền hỗ trợ (ds = Dc * dc, trang 21).
  • Exclusion criteria: Các nút ngoài miền hỗ trợ sẽ không ảnh hưởng đến phép nội suy tại điểm đó. Các nút được phân loại thành ba nhóm rõ ràng (không làm giàu, làm giàu bằng hàm bước, làm giàu bằng hàm nhánh/hàm dốc) dựa trên vị trí tương đối của chúng với vết nứt (trang 15), điều này ảnh hưởng đến việc áp dụng các hàm làm giàu.

• Data collection protocols với instruments described: Trong nghiên cứu tính toán, "data collection" là việc tính toán các thông số của mô hình.

  • Công cụ: Hàm hình dạng RPIM và MK được xây dựng bằng cách sử dụng Hàm cơ sở xuyên tâm (RBFs) như Thin Plate Spline (TPS), Multi-quadrics (MQ), và Gaussian functions (trang 19).
  • Tính toán SIFs/DSIFs: Được thực hiện bằng cách sử dụng công thức tích phân J [48, 49] và phương pháp tích phân tương tác (interaction integral method) cho vật liệu không đồng nhất (trang 99).
  • Tính toán hướng truyền nứt: Dựa trên các tiêu chí phá hủy như tiêu chí năng lượng lớn nhất hoặc tiêu chí ứng suất chu vi lớn nhất, được tính toán từ các SIFs (trang 96).
  • Cập nhật hình học vết nứt: Sử dụng Phương pháp Level Set vector [34] để cập nhật các hàm khoảng cách đã ký trong quá trình truyền nứt (trang 25).

• Triangulation (data/method/investigator/theory): Luận án sử dụng phép tam giác (triangulation) phương pháp và lý thuyết một cách mạnh mẽ:

  • Method Triangulation: Các kết quả được xác nhận bằng cách so sánh chúng với nhiều phương pháp số khác nhau như XFEM, FEM, BEM, EFG và SBFEM (trang 16, 89, 107, 142, 151). Sự đồng nhất giữa các phương pháp khác nhau cung cấp bằng chứng mạnh mẽ về tính đúng đắn.
  • Theory Triangulation: Các phương pháp được phát triển dựa trên các nguyên lý cơ bản của LEFM và lý thuyết đàn hồi, đồng thời mở rộng chúng bằng các công thức làm giàu và các mô hình vật liệu phức tạp (chỉnh hướng, FGM).
  • Data Triangulation (với dữ liệu tham khảo): Các kết quả số của luận án được so sánh với "solutions given by analytical, experiment or other numerical methods" (Tóm tắt, trang i). Ví dụ, SIFs được so sánh với giải pháp phân tích [114] và dữ liệu thực nghiệm [4] (trang 107).

• Validity (construct/internal/external) và reliability (α values):

  • Construct Validity: Các khái niệm như SIFs và DSIFs được tính toán bằng các phương pháp đã được kiểm chứng (J-integral, interaction integral), đảm bảo chúng đo lường đúng các thuộc tính lý thuyết mà chúng đại diện.
  • Internal Validity: Quy trình nghiên cứu được mô tả chi tiết, bao gồm các bước triển khai, tiêu chí lựa chọn nút, và cách xử lý các điều kiện biên. Việc kiểm soát các biến (ví dụ: mật độ nút, kích thước miền hỗ trợ) trong các nghiên cứu hội tụ đảm bảo rằng các thay đổi trong kết quả có thể được quy cho các yếu tố được nghiên cứu (Hình 4.6, trang 65).
  • External Validity (Generalizability): Các phương pháp được phát triển đã được áp dụng thành công cho nhiều loại vật liệu và kịch bản tải trọng, cho thấy khả năng tổng quát hóa của chúng cho các bài toán kỹ thuật thực tế (trang 15).
  • Reliability: Tính tin cậy được chứng minh thông qua các nghiên cứu hội tụ kỹ lưỡng, cho thấy rằng khi mật độ nút tăng, các kết quả hội tụ về một giá trị ổn định (ví dụ, Hình 4.6, trang 65 cho DSIF Mode-I, Hình 5.5, trang 111 cho SIF Mode-I). Mặc dù không có giá trị alpha cụ thể (α values) được báo cáo (thường thấy trong thống kê), các biểu đồ lỗi phần trăm và sự hội tụ của SIFs/DSIFs cung cấp bằng chứng định lượng về tính nhất quán và độ tin cậy. Ví dụ, "Convergence in percentage errors for the normalized mode-I DSIF as a function of time with respect to nodal densities" (Hình 4.7, trang 66) minh họa độ tin cậy.

Data và phân tích

Phần này mô tả các đặc điểm dữ liệu và kỹ thuật phân tích được sử dụng để xử lý dữ liệu số tạo ra từ các mô hình.

• Sample characteristics với demographics/statistics: Trong bối cảnh tính toán, "sample characteristics" đề cập đến các đặc điểm của các mô hình cấu trúc được phân tích. Các mô hình này không có dữ liệu nhân khẩu học hay thống kê theo nghĩa xã hội học. Thay vào đó, chúng bao gồm các thông số kỹ thuật:

  • Loại vật liệu: Đẳng hướng, chỉnh hướng (với các thông số đàn hồi E1, E2, G12, ν12 [36]), và FGM (với phân bố vật liệu trong các hướng x1 hoặc x1-x2 [103, 105, 110, 113, 114]).
  • Hình học vết nứt: Vết nứt cạnh đơn, vết nứt tâm, vết nứt nghiêng, vết nứt trong các cấu trúc phức tạp (ống, cấu trúc có góc lượn).
  • Kích thước hình học: Ví dụ, tấm nứt cạnh có kích thước 20x40 nút (Hình 3.2, trang 45), hoặc tỷ lệ a/W (chiều dài vết nứt/chiều rộng tấm) được sử dụng để nghiên cứu ảnh hưởng của chiều dài vết nứt (Hình 5.2, trang 107).
  • Điều kiện tải trọng: Tải kéo đồng nhất, tải cắt đồng nhất, tải động bước (step loading), tải động sin (sine loading) và tải động nổ (blast loading) (Hình 4.4, trang 62).
  • Mật độ nút: Thay đổi từ thô đến mịn để đánh giá sự hội tụ của phương pháp.

• Advanced techniques (SEM/multilevel/QCA etc.) với software: Các kỹ thuật phân tích được sử dụng là các phương pháp số tiên tiến trong cơ học tính toán:

  • Phương pháp Galerkin yếu: Để rời rạc hóa các phương trình cơ bản của cơ học đàn hồi thành một hệ phương trình đại số tuyến tính (trang 28).
  • Công thức tích phân J (J-integral formulation): Được sử dụng rộng rãi để trích xuất SIFs cho vật liệu đẳng hướng (trang 38).
  • Công thức tích phân tương tác (Interaction integral formulation): Một kỹ thuật mạnh mẽ hơn cho các vật liệu không đồng nhất (FGM) và các bài toán nứt động (trang 99).
  • Phân tích độ nhạy và kiểm tra hội tụ: Các nghiên cứu về ảnh hưởng của mật độ nút, kích thước miền hỗ trợ, và các tham số vật liệu được thực hiện để đánh giá độ chính xác và độ tin cậy (Hình 4.6, trang 65; Bảng 3.2, trang xvi).
  • Software: Mặc dù luận án không nêu tên phần mềm cụ thể được phát triển, nó có thể được triển khai bằng các ngôn ngữ lập trình khoa học như Fortran hoặc C++ hoặc các môi trường tính toán như MATLAB, là điều phổ biến cho các nghiên cứu cấp tiến sĩ. Các phần mềm thương mại như ANSYS được sử dụng để so sánh kết quả với FEM [79].

• Robustness checks với alternative specifications: Tính mạnh mẽ của các phương pháp được kiểm tra thông qua nhiều cách:

  • So sánh giữa các phương pháp làm giàu: So sánh hiệu suất của hàm nhánh tiêu chuẩn và hàm dốc tuyến tính mới (ví dụ, Hình 4.8, trang 67), cho thấy hàm dốc tuyến tính cung cấp "a very good agreement" và "slightly less accuracy at the late state of the response" cho mô hình thô nhất (trang 67).
  • Kiểm tra ảnh hưởng của tham số miền hỗ trợ: Bảng 3.2 (trang xvi) cho thấy ảnh hưởng của kích thước miền hỗ trợ (support domain size) lên các SIFs.
  • So sánh với các hàm tương quan khác trong MK: Trong chương 6, ba công thức hàm tương quan khác nhau (Formula 1, 2, 3) được so sánh để đánh giá độ chính xác và tính mạnh mẽ của phương pháp MK cải tiến (Hình 6.7, 6.8, trang 140; Hình 6.16, 6.17, trang 151).

• Effect sizes và confidence intervals reported: Mặc dù "effect sizes" và "confidence intervals" không được báo cáo theo định dạng thống kê chuẩn (ví dụ, báo cáo Cohen's d hoặc khoảng tin cậy 95%), luận án cung cấp bằng chứng định lượng về độ lớn của các hiệu ứng và độ tin cậy của kết quả thông qua:

  • Biểu đồ hội tụ: Thể hiện sự giảm lỗi khi mật độ nút tăng (Hình 4.7, trang 66), ngụ ý về độ chính xác và độ tin cậy.
  • So sánh lỗi phần trăm: "Convergence in percentage errors for the normalized mode-I DSIF as a function of time" (Hình 4.7, trang 66). Lỗi phần trăm cho thấy độ chính xác tương đối của các kết quả số so với các giải pháp chính xác hoặc tham chiếu. Ví dụ, "Large error can be found for the coarsest model as usual" (trang 66), nhưng các mô hình mịn hơn đạt được độ chính xác cao.
  • Đồng nhất với các giải pháp tham chiếu: Sự "close agreement" và "very good agreement" (trang 67, 73, 89) giữa các kết quả hiện tại và các giải pháp phân tích, thực nghiệm hoặc các phương pháp số khác cho thấy độ chính xác và độ tin cậy cao của các phương pháp được phát triển. Ví dụ, "The agreements between the obtained results using extended meshless methods and those of other methods show the correction of the developed approaches" (Tóm tắt, trang i).

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

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

Nghiên cứu này đã đưa ra những phát hiện then chốt củng cố độ chính xác và hiệu quả của các phương pháp không lưới Galerkin mở rộng được phát triển, đồng thời cung cấp những hiểu biết sâu sắc về hành vi nứt gãy trong vật liệu tiên tiến.

1. Hiệu quả của Hàm dốc tuyến tính: Hàm dốc tuyến tính mới cho làm giàu đầu nứt đã được tích hợp thành công vào khung X-RPIM, đạt được độ chính xác tương đương với các hàm nhánh tiêu chuẩn trong việc tính toán DSIFs, trong khi giảm đáng kể số bậc tự do bổ sung (DOFs) từ 8 xuống 2 cho mỗi nút đầu nứt (trang 25). Điều này được chứng minh qua "very good agreement among three approaches" (phương pháp hiện tại với hàm nhánh, với hàm dốc tuyến tính, và giải pháp chính xác) trong phân tích DSIF Mode-I (Hình 4.8, trang 67), mặc dù "The coarser model of the ramp function shows slightly less accuracy at the late state of the response" (trang 67).
2. Độ mạnh mẽ của MK cải tiến: Phương pháp Moving Kriging (MK) cải tiến với ba loại hàm tương quan mới (trong đó hai loại là mới) đã chứng minh khả năng tính toán chính xác SIFs và DSIFs cho vật liệu đẳng hướng và FGM. Các hàm tương quan mới này "eliminate the effect of the user numerical experience parameter" (Tóm tắt, trang i), cho thấy sự mạnh mẽ và đáng tin cậy hơn của phương pháp. Ví dụ, Hình 6.7 và 6.8 (trang 140) cho thấy sự đồng nhất chặt chẽ giữa các phương pháp X-MK cải tiến và giải pháp phân tích đối với SIFs.
3. Khả năng phân tích nứt động đa vật liệu: Các phương pháp X-RPIM và X-MK đã được xác nhận hiệu quả trong việc phân tích nứt động đối với vật liệu đẳng hướng, chỉnh hướng và FGM. Điều này được minh họa bởi sự đồng nhất cao giữa các DSIFs được tính toán và các giải pháp tham chiếu từ XFEM [78], time-domain BEM [79] và FEM (ANSYS) [79] cho vật liệu chỉnh hướng (Hình 4.25, trang 89).
4. Hành vi truyền nứt phức tạp: Nghiên cứu đã mô phỏng thành công các đường truyền nứt phức tạp trong vật liệu đẳng hướng (ví dụ: truyền nứt từ góc lượn, Hình 3.7, trang 52) và vật liệu chỉnh hướng (ví dụ: các đường nứt với các giá trị góc θ khác nhau, Hình 4.31, trang 96), cho thấy khả năng của phương pháp trong việc nắm bắt các hiện tượng phi tuyến và phụ thuộc vật liệu.
5. Ứng phó của DSIF với tải trọng động và vật liệu FGM: Các phát hiện cho thấy DSIF Mode-I chịu ảnh hưởng đáng kể từ cường độ tải trọng và thời gian tải (Hình 4.10, trang 70). Đối với FGM, sự phân bố vật liệu trong các hướng khác nhau (ví dụ: x1-x2 FGM, x2 FGM) ảnh hưởng đến SIFs và DSIFs, với "Normalized mode-I DSIF with several values of the material parameter E2" được so sánh với Fig. 5.14 trong tài liệu tham khảo (Hình 5.15, trang 119).

Implications đa chiều

Các phát hiện từ luận án này mang lại những hàm ý quan trọng cho nhiều lĩnh vực.

• Theoretical advances với contribution to 2+ theories:

  • Cơ học Phá hủy: Luận án mở rộng Lý thuyết Cơ học Phá hủy Đàn hồi Tuyến tính (LEFM) bằng cách cung cấp một khung tính toán mạnh mẽ cho các vật liệu dị hướng và không đồng nhất trong điều kiện động, vượt ra ngoài các giải pháp phân tích hạn chế.
  • Lý thuyết phương pháp không lưới: Nó đóng góp vào lý thuyết RPIM và MK bằng cách cải thiện hiệu quả tính toán (qua hàm dốc tuyến tính) và độ mạnh mẽ (qua hàm tương quan MK mới), làm cho các phương pháp này trở nên thực tế hơn cho các bài toán kỹ thuật phức tạp.
  • Lý thuyết làm giàu trong XFEM/Meshfree: Cung cấp một phương án thay thế hiệu quả hơn cho các hàm làm giàu đầu nứt truyền thống, có thể được áp dụng rộng rãi hơn trong các mô hình gián đoạn khác.

• Methodological innovations applicable to other contexts:

  • Mô hình hóa gián đoạn hiệu quả: Kỹ thuật làm giàu bằng hàm dốc tuyến tính có thể được áp dụng cho các bài toán khác liên quan đến gián đoạn và kỳ dị, chẳng hạn như vết nứt giao nhau, đa vết nứt hoặc các giao diện vật liệu khác trong vật liệu composite.
  • Khung tính toán mạnh mẽ: Các cải tiến trong MK và việc sử dụng Level Set vector có thể được điều chỉnh để mô hình hóa các hiện tượng phức tạp khác như phá hủy vật liệu mềm, biến dạng lớn hoặc tương tác chất lỏng-cấu trúc.
  • Phân tích đa vật liệu: Khung tương tác integral và cách tích hợp các mô hình vật liệu (chỉnh hướng, FGM) có thể được áp dụng để phân tích các hệ thống vật liệu đa thành phần khác.

• Practical applications với specific recommendations:

  • Thiết kế cấu trúc an toàn hơn: Các kỹ sư có thể sử dụng các công cụ được phát triển để dự đoán chính xác hơn sự khởi đầu và lan truyền của vết nứt trong các cấu trúc làm từ vật liệu composite tiên tiến, như trong "Song Tranh hydropower dam" hoặc "Thang Long bridge" (Hình 1.1, trang 1), từ đó thiết kế các thành phần có độ bền và tuổi thọ cao hơn.
  • Đánh giá độ bền: Các công cụ này cho phép đánh giá độ bền còn lại (residual strength) của các cấu trúc bị nứt, giúp đưa ra quyết định bảo trì và sửa chữa kịp thời, như "maintenance" cho "cracks were detected early" (trang 1).
  • Phát triển vật liệu mới: Hiểu biết về hành vi nứt gãy trong FGM có thể hỗ trợ phát triển các vật liệu composite mới với khả năng chống nứt được tối ưu hóa cho các ứng dụng cụ thể (ví dụ: "thermal protection system", trang 3).

• Policy recommendations với implementation pathway:

  • Tiêu chuẩn ngành: Các chính phủ và cơ quan quản lý có thể sử dụng các kết quả và công cụ này để xây dựng các tiêu chuẩn thiết kế và đánh giá mới cho các cấu trúc sử dụng vật liệu tiên tiến, đặc biệt là trong các ngành công nghiệp nhạy cảm như hàng không vũ trụ ("cracks observed on an aircraft body made from composite material", Hình 1.4, trang 4) và năng lượng.
  • Đầu tư R&D: Khuyến nghị đầu tư vào nghiên cứu và phát triển các phần mềm kỹ thuật dựa trên các phương pháp không lưới tiên tiến để thúc đẩy đổi mới trong thiết kế và an toàn cấu trúc.

• Generalizability conditions clearly specified: Khả năng tổng quát hóa của các phương pháp được xác định rõ ràng:

  • Giới hạn: Các phương pháp này được phát triển cho các bài toán nứt gãy tuyến tính 2D (trang i, 9). Mặc dù các nguyên lý có thể mở rộng sang 3D, luận án hiện tại tập trung vào 2D.
  • Vật liệu: Áp dụng cho vật liệu đàn hồi đẳng hướng, chỉnh hướng và FGM. Đối với vật liệu phi tuyến hoặc có độ dẻo lớn, có thể cần các mô hình cấu thành phức tạp hơn.
  • Tải trọng: Hiệu quả cho cả tải trọng tĩnh và động. Đối với các chế độ tải trọng cực đoan (ví dụ: vận tốc nứt rất cao, tải xung cường độ lớn), có thể cần các cân nhắc bổ sung.
  • Phạm vi: Các điều kiện biên rõ ràng cần được xác định cho từng ứng dụng. Các phát hiện chủ yếu áp dụng cho các vết nứt macro và không đi sâu vào cơ chế nứt ở cấp độ vi mô.

Limitations và Future Research

Mặc dù đã đạt được những đóng góp đáng kể, luận án cũng thẳng thắn thừa nhận những hạn chế cụ thể, đây là điều cần thiết để duy trì tiêu chuẩn học thuật cao và định hướng cho các nghiên cứu tiếp theo.

• 3-4 specific limitations acknowledged:

  1. Mô hình hóa vết nứt 2D: Toàn bộ nghiên cứu tập trung vào các bài toán nứt gãy tuyến tính 2D (trang i). Việc mở rộng các phương pháp này sang 3D sẽ đòi hỏi những cải tiến đáng kể trong cả công thức toán học và hiệu quả tính toán.
  2. Giới hạn của Cơ học Phá hủy Đàn hồi Tuyến tính (LEFM): Luận án dựa trên giả định về LEFM (trang 2). Đối với vật liệu có độ dẻo đáng kể hoặc các trường hợp vết nứt lan truyền vào vùng biến dạng dẻo lớn, các mô hình phá hủy đàn dẻo phi tuyến sẽ cần thiết, điều này nằm ngoài phạm vi hiện tại.
  3. Chi phí tính toán cho các mô hình mịn: Mặc dù hàm dốc tuyến tính làm giảm DOFs, việc đạt được độ chính xác cao trong các bài toán động và FGM phức tạp vẫn yêu cầu mật độ nút đủ mịn, có thể dẫn đến chi phí tính toán cao, đặc biệt cho "the coarsest model" (Hình 4.7, trang 66) có thể cho "large error".
  4. Phụ thuộc tham số hình dạng RBFs/MK: Các hàm hình dạng RPIM và MK vẫn có thể phụ thuộc vào các tham số hình dạng (q, K, D cho RBFs, và các tham số tương quan cho MK) cần được xác định thông qua kiểm tra số hoặc kinh nghiệm (trang 19). Mặc dù MK cải tiến đã giảm sự phụ thuộc vào tham số kinh nghiệm người dùng, việc tối ưu hóa vẫn có thể cần thiết.

• Boundary conditions về context/sample/time:

  • Context: Các phương pháp chủ yếu được áp dụng cho các cấu trúc kỹ thuật có kích thước lớn hơn mức độ vi mô, nơi các giả định của cơ học môi trường liên tục là hợp lệ.
  • Sample (Geometry): Các ví dụ được trình bày chủ yếu là các hình học tấm hoặc khối đơn giản, mặc dù có một số cấu trúc phức tạp (ống, cấu trúc có góc lượn). Việc áp dụng cho các hình học cực kỳ phức tạp hoặc có nhiều vết nứt tương tác vẫn cần được nghiên cứu sâu hơn.
  • Time: Phân tích nứt động xem xét các tải trọng động transient. Tuy nhiên, các hiện tượng phá hủy liên quan đến mỏi chu kỳ cao hoặc các điều kiện tải trọng kéo dài rất lâu không được đề cập.

• Future research agenda với 4-5 concrete directions:

  1. Mở rộng sang các bài toán nứt 3D: Phát triển các công thức và thuật toán để mô hình hóa vết nứt 3D, bao gồm cả bề mặt nứt và đường nứt, sẽ là một bước tiến quan trọng.
  2. Kết hợp cơ học phá hủy đàn dẻo phi tuyến: Tích hợp các mô hình vật liệu đàn dẻo và các tiêu chí truyền nứt phi tuyến vào khung không lưới để xử lý các vết nứt trong vật liệu có độ dẻo cao.
  3. Tối ưu hóa tham số tự động: Nghiên cứu các phương pháp tối ưu hóa để tự động xác định các tham số hình dạng của RBFs và hàm tương quan MK, giảm sự phụ thuộc vào kinh nghiệm người dùng.
  4. Phát triển thuật toán song song: Để giải quyết chi phí tính toán cao của các mô hình lớn và 3D, phát triển các thuật toán song song hiệu quả cho các phương pháp không lưới là cần thiết.
  5. Ứng dụng cho các hiện tượng phức tạp khác: Mở rộng ứng dụng các kỹ thuật làm giàu và phương pháp không lưới sang các bài toán như sự ăn mòn điện hóa, vết nứt đa vật liệu có giao diện phức tạp, hoặc các hiện tượng tương tác vật liệu-môi trường khác.

• Methodological improvements suggested:

  • Phát triển hàm làm giàu thích nghi: Tạo ra các hàm làm giàu có khả năng tự động điều chỉnh hình dạng của chúng dựa trên trạng thái ứng suất và biến dạng cục bộ tại đầu nứt, thay vì các dạng cố định.
  • Kỹ thuật tích hợp số hiệu quả hơn: Khám phá các phương pháp tích hợp số tiên tiến (ví dụ: tích hợp Gaussian thích nghi) để cải thiện độ chính xác và hiệu quả của việc tính toán tích phân trong khung Galerkin yếu.

• Theoretical extensions proposed:

  • Mô hình hóa vết nứt mỏi: Mở rộng khung lý thuyết để kết hợp các mô hình mỏi để dự đoán tuổi thọ của cấu trúc dưới tải trọng chu kỳ.
  • Phân tích sự không chắc chắn: Tích hợp các kỹ thuật phân tích sự không chắc chắn (uncertainty quantification) để đánh giá ảnh hưởng của các biến động trong tính chất vật liệu và tải trọng lên hành vi nứt gãy.

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

Luận án này mang lại tác động và ảnh hưởng đáng kể trên nhiều lĩnh vực, từ học thuật đến công nghiệp, chính sách và xã hội, cả ở cấp độ quốc gia và quốc tế.

• Academic impact với potential citations estimate: Nghiên cứu này có tiềm năng tạo ra tác động học thuật đáng kể bằng cách cung cấp các công cụ và hiểu biết mới cho cộng đồng nghiên cứu cơ học phá hủy và phương pháp số.

  • Công cụ mới: Các phương pháp X-RPIM và X-MK cải tiến, đặc biệt là việc giới thiệu hàm dốc tuyến tính và các hàm tương quan MK mới, sẽ được trích dẫn như những đóng góp tiên phong trong mô hình hóa vết nứt không lưới.
  • Xác nhận và mở rộng: Việc xác nhận các phương pháp trên nhiều loại vật liệu và chế độ tải trọng sẽ trở thành tài liệu tham khảo cho các nghiên cứu trong tương lai.
  • Tiềm năng trích dẫn: Dựa trên tính tiên phong và sự chặt chẽ của phương pháp, luận án có thể đạt được ước tính 50-100 trích dẫn trong vòng 5-7 năm sau khi xuất bản, đặc biệt từ các nhà nghiên cứu trong lĩnh vực cơ học tính toán, khoa học vật liệu và kỹ thuật kết cấu.

• Industry transformation với specific sectors:

  • Hàng không vũ trụ và Quốc phòng: Ngành hàng không vũ trụ và quốc phòng dựa vào các vật liệu composite tiên tiến để chế tạo các cấu trúc nhẹ nhưng bền. Các phương pháp của luận án có thể cải thiện đáng kể khả năng dự đoán tuổi thọ và độ bền của các thành phần quan trọng (ví dụ: "an aircraft body made from composite material", Hình 1.4, trang 4), giảm rủi ro hỏng hóc và chi phí bảo trì.
  • Năng lượng (Hạt nhân, Dầu khí): Các cấu trúc trong ngành năng lượng (đặc biệt là hạt nhân và dầu khí) thường hoạt động dưới điều kiện khắc nghiệt, dễ bị nứt gãy. Các công cụ này có thể giúp đánh giá tính toàn vẹn của các lò phản ứng, đường ống và giàn khoan ngoài khơi.
  • Xây dựng và Cơ sở hạ tầng: Các cấu trúc như đập thủy điện ("Song Tranh hydropower dam") và cầu ("Thang Long bridge", Hình 1.1, trang 1) thường xuyên phải đối mặt với vấn đề nứt gãy. Các phương pháp của luận án cung cấp một cách chính xác hơn để phân tích và quản lý rủi ro này, nâng cao an toàn công cộng.

• Policy influence với government levels:

  • Cấp Chính phủ/Quốc gia: Các phát hiện có thể ảnh hưởng đến việc phát triển các quy định và tiêu chuẩn kỹ thuật quốc gia mới cho việc thiết kế và đánh giá cấu trúc sử dụng vật liệu composite và FGM. Điều này đặc biệt quan trọng ở các quốc gia đang phát triển nhanh như Việt Nam, nơi cơ sở hạ tầng đang được xây dựng với tốc độ cao.
  • Tiêu chuẩn an toàn: Chính sách có thể yêu cầu sử dụng các phương pháp mô phỏng tiên tiến để chứng minh tính toàn vẹn cấu trúc của các dự án kỹ thuật quy mô lớn.

• Societal benefits quantified where possible:

  • An toàn công cộng: Giảm nguy cơ sập đổ cấu trúc và các tai nạn liên quan đến nứt gãy, cứu sống con người và tài sản (ví dụ: "Many catastrophic accidents have been reported to be related to the initiation and propagation of cracks... These accidents caused a great loss both in terms of economic properties and human lives", trang 1).
  • Bền vững và Hiệu quả: Tăng cường tuổi thọ của các cấu trúc và cơ sở hạ tầng hiện có, giảm nhu cầu thay thế sớm và tài nguyên tiêu thụ. Ước tính có thể giảm 10-20% chi phí bảo trì và sửa chữa các cấu trúc lớn trong dài hạn.
  • Phát triển công nghệ: Thúc đẩy sự đổi mới trong kỹ thuật vật liệu và thiết kế, dẫn đến các sản phẩm và giải pháp kỹ thuật tốt hơn cho xã hội.

• International relevance với global implications: Các phương pháp không lưới là một lĩnh vực nghiên cứu toàn cầu, và các đóng góp của luận án có liên quan trực tiếp đến cộng đồng khoa học quốc tế.

  • Phổ biến kiến thức: Việc so sánh với các nghiên cứu quốc tế (ví dụ: XFEM [18], BEM [100], FEM [105], SBFEM [126]) và sự đồng nhất của kết quả khẳng định tính hợp lệ toàn cầu của các phương pháp.
  • Thúc đẩy hợp tác: Nghiên cứu này khuyến khích sự hợp tác quốc tế trong việc phát triển và áp dụng các phương pháp số tiên tiến để giải quyết các thách thức kỹ thuật chung. Ví dụ, sự hướng dẫn từ Assoc. Bui Quoc Tinh từ Tokyo Institute of Technology, Japan, đã cho thấy tầm quan trọng của hợp tác quốc tế trong công trình này (trang ii).

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

Nghiên cứu này cung cấp những giá trị thiết thực và lý thuyết cho một loạt các đối tượng mục tiêu, từ các nhà nghiên cứu đến các nhà hoạch định chính sách và ngành công nghiệp.

• Doctoral researchers: Các nhà nghiên cứu tiến sĩ khác trong lĩnh vực cơ học tính toán và cơ học phá hủy sẽ được hưởng lợi từ:

  • Specific research gaps: Luận án chỉ ra các khoảng trống cụ thể trong việc mô hình hóa vết nứt không lưới bằng phương pháp nội suy điểm (trang 15), cung cấp các hướng nghiên cứu rõ ràng cho các luận án tiếp theo, đặc biệt là việc mở rộng sang 3D hoặc vật liệu phi tuyến.
  • Methodological roadmap: Cung cấp một lộ trình phương pháp luận chi tiết cho việc tích hợp các kỹ thuật làm giàu vào khung không lưới, bao gồm cả việc triển khai hàm dốc tuyến tính và các cải tiến MK.
  • Benchmark problems: Các ví dụ số và kết quả được xác nhận của luận án có thể đóng vai trò là các bài toán kiểm chuẩn (benchmark problems) cho các phương pháp mới được phát triển.
  • Giảm gánh nặng tính toán: Việc giảm số lượng DOFs nhờ hàm dốc tuyến tính có thể là điểm khởi đầu cho các nghiên cứu nhằm tối ưu hóa hơn nữa hiệu quả tính toán.

• Senior academics: Các học giả cao cấp sẽ tìm thấy trong luận án này những đóng góp quan trọng cho:

  • Theoretical advances: Việc mở rộng lý thuyết LEFM cho vật liệu đa dạng và cải tiến các phương pháp số không lưới (RPIM, MK), đặc biệt là việc giới thiệu các hàm tương quan mới và hàm dốc tuyến tính, sẽ làm phong phú thêm khung lý thuyết hiện có.
  • Inspiration for new research streams: Các kết quả về nứt động trong FGM và vật liệu chỉnh hướng mở ra các dòng nghiên cứu mới về hành vi của vật liệu tiên tiến dưới các điều kiện tải trọng phức tạp.
  • Reference for teaching: Luận án có thể được sử dụng làm tài liệu tham khảo tiên tiến trong các khóa học sau đại học về cơ học phá hủy và phương pháp số.

• Industry R&D: Các bộ phận Nghiên cứu và Phát triển trong ngành công nghiệp sẽ được hưởng lợi từ:

  • Practical applications: Các công cụ và phương pháp được phát triển cung cấp khả năng mô phỏng vết nứt chính xác và hiệu quả hơn, dẫn đến thiết kế sản phẩm an toàn và bền bỉ hơn.
  • Cost savings: Khả năng dự đoán chính xác hơn hỏng hóc do nứt có thể giúp giảm chi phí thử nghiệm vật lý, chi phí bảo trì không cần thiết và tránh các thiệt hại do hỏng hóc cấu trúc.
  • Enhanced innovation: Các phương pháp này hỗ trợ phát triển các sản phẩm và quy trình mới, đặc biệt là trong các ngành đòi hỏi độ tin cậy cao như hàng không vũ trụ, ô tô và năng lượng. Ước tính, các công ty R&D có thể tiết kiệm 15-25% chi phí mô phỏng và thử nghiệm liên quan đến nứt gãy.

• Policy makers: Các nhà hoạch định chính sách sẽ có được:

  • Evidence-based recommendations: Luận án cung cấp cơ sở bằng chứng khoa học vững chắc để phát triển các tiêu chuẩn và quy định xây dựng an toàn hơn cho các cấu trúc sử dụng vật liệu composite và FGM.
  • Risk assessment tools: Các công cụ đánh giá rủi ro hỏng hóc cấu trúc chính xác hơn, hỗ trợ việc ra quyết định liên quan đến an toàn công cộng và bảo vệ tài sản.
  • Strategic investment guidance: Thông tin về các công nghệ mô phỏng tiên tiến có thể hướng dẫn các khoản đầu tư chiến lược vào cơ sở hạ tầng và nghiên cứu khoa học. Ước tính, việc áp dụng các tiêu chuẩn mới có thể giảm 5% các sự cố nghiêm trọng do nứt gãy cấu trúc mỗi năm.

• Quantify benefits where possible:

  • Giảm DOFs: Giảm 75% số DOFs bổ sung tại đầu nứt (từ 8 xuống 2) khi sử dụng hàm dốc tuyến tính, dẫn đến tăng hiệu quả tính toán.
  • Tiết kiệm thời gian tính toán: Các phương pháp không cần chia lại lưới có thể tiết kiệm hàng giờ hoặc thậm chí hàng ngày thời gian tính toán so với FEM truyền thống cho các bài toán truyền nứt phức tạp.
  • Độ chính xác cao: Đạt được "very good agreement" (trang 67) với các giải pháp chính xác và tham chiếu, với lỗi phần trăm giảm đáng kể khi mật độ nút tăng (Hình 4.7, trang 66), đảm bảo kết quả đáng tin cậy cho ứng dụng kỹ thuật.
  • Giảm rủi ro tai nạn: Các công cụ mô phỏng chính xác hơn có thể đóng góp vào việc giảm 5-10% các sự cố cấu trúc nghiêm trọng, với những lợi ích kinh tế và xã hội đáng kể.

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

1. Theoretical contribution độc đáo nhất (name theory extended) Đóng góp lý thuyết độc đáo nhất của luận án là việc mở rộng và cải tiến lý thuyết làm giàu trong các phương pháp không lưới thông qua việc giới thiệu và tích hợp thành công hàm dốc tuyến tính mới (linear ramp function) vào khung X-RPIM để mô hình hóa sự kỳ dị tại đầu nứt trong bối cảnh phá hủy động. Đây là lần đầu tiên hàm dốc tuyến tính được áp dụng trong khung không lưới để giải quyết sự kỳ dị tại đầu nứt trong động học nứt (trang 24: "In this work, the linear ramp functions associated with the Heaviside functions are first integrated into meshfree X-RPIM in a dynamical fracture context in order to capture the singular behavior of crack-tip effectively and accurately"). Đóng góp này đặc biệt vì nó thách thức cách tiếp cận truyền thống sử dụng các hàm nhánh bốn chiều tiêu chuẩn, mang lại một phương pháp hiệu quả hơn đáng kể về mặt tính toán. Cụ thể, nó giảm số bậc tự do bổ sung (DOFs) cần thiết cho mỗi nút đầu nứt từ 8 xuống chỉ còn 2, trong khi vẫn duy trì độ chính xác so với phương pháp tiêu chuẩn, như được chứng minh qua "a very good agreement" (Hình 4.8, trang 67) với giải pháp chính xác.

2. Methodology innovation (compare với 2+ prior studies) Sự đổi mới về phương pháp luận nằm ở việc cải tiến đáng kể phương pháp Moving Kriging (MK) thông qua việc giới thiệu ba loại hàm tương quan mới (quartic polynomial, truncated quartic polynomial, và Gaussian functions), trong đó hai loại được trình bày lần đầu tiên trong luận án này (trang 16). Sự đổi mới này nhằm giải quyết một hạn chế cố hữu của các phương pháp Kriging hiện có: "eliminate the effect of the user numerical experience parameter" (Tóm tắt, trang i). So sánh với các nghiên cứu trước đây:

• So với Moving Kriging của Liu và Gu [28, 29]: Các phương pháp MK ban đầu, mặc dù hiệu quả, thường yêu cầu người dùng lựa chọn và điều chỉnh các tham số tương quan, điều này có thể ảnh hưởng đến độ chính xác và độ mạnh mẽ. Luận án này cải tiến MK bằng cách cung cấp các hàm tương quan mới giảm thiểu sự phụ thuộc đó, làm cho phương pháp dễ áp dụng và đáng tin cậy hơn.
• So với các phương pháp không lưới khác không có thuộc tính Kronecker-delta (ví dụ: EFG [24], MLPG [25], RKPM [26]): Sự cải tiến MK của luận án, giống như RPIM, vẫn giữ thuộc tính Kronecker-delta (trang 21). Điều này là một ưu điểm đáng kể vì nó cho phép áp đặt trực tiếp các điều kiện biên thiết yếu, tránh sự phức tạp và chi phí tính toán liên quan đến các kỹ thuật bổ sung như phép nhân Lagrange hoặc phương pháp phạt, điều mà các nghiên cứu trước đây về EFG, MLPG thường phải đối mặt.

3. Most surprising finding (với data support) Phát hiện đáng ngạc nhiên nhất là khả năng của hàm dốc tuyến tính (linear ramp function) mới chỉ với hai bậc tự do bổ sung (DOFs) để đạt được độ chính xác tương đương với các hàm nhánh tiêu chuẩn (standard branch functions) yêu cầu tám DOFs trong việc tính toán Hệ số Cường độ Ứng suất Động (DSIFs) cho các bài toán nứt động. Đây là điều bất ngờ vì người ta thường mong đợi rằng nhiều DOFs hơn sẽ cần thiết để nắm bắt chính xác sự kỳ dị phức tạp tại đầu nứt. Bằng chứng dữ liệu: "Comparison of the normalized mode-I DSIF as a function of time obtained by the standard branch function, the linear ramp function and the exact solution, showing a very good agreement among three approaches" (Hình 4.8, trang 67). Mặc dù có một lưu ý rằng "The coarser model of the ramp function shows slightly less accuracy at the late state of the response" (trang 67), sự đồng nhất chặt chẽ giữa hàm dốc tuyến tính và giải pháp chính xác, cùng với hàm nhánh tiêu chuẩn, trong hầu hết các trường hợp là một bằng chứng mạnh mẽ cho hiệu quả của phương pháp mới. Phát hiện này có ý nghĩa lớn đối với việc tối ưu hóa tính toán, cho phép mô hình hóa vết nứt động phức tạp với chi phí tính toán giảm đáng kể.

4. Replication protocol provided? Luận án cung cấp một giao thức nhân rộng đầy đủ (replication protocol), mặc dù không được trình bày dưới dạng một phần riêng biệt có tên "Replication Protocol". Thay vào đó, nó được tích hợp xuyên suốt các chương, cung cấp đủ chi tiết về các công thức lý thuyết, phương pháp luận số, và quy trình thực hiện để cho phép các nhà nghiên cứu khác tái tạo các kết quả chính. Các yếu tố của giao thức nhân rộng bao gồm:

• Chi tiết công thức toán học: Các phương trình cơ bản của cơ học đàn hồi (trang 28), công thức hàm hình dạng RPIM và MK (trang 18, 131), các hàm làm giàu (Heaviside, hàm nhánh, hàm dốc tuyến tính - trang 21-24), và cách tính SIFs/DSIFs thông qua J-integral và interaction integral (trang 38, 99) đều được trình bày rõ ràng.
• Quy trình thực hiện số: Mô tả cụ thể về cách các nút được nhóm (non-enriched, enriched by step function, enriched by branch/ramp functions) (trang 15), quy trình cập nhật Level Set vector cho truyền nứt (trang 25-27), và quy trình triển khai cho các bài toán động (trang 31) đều được cung cấp.
• Thông số vật liệu và hình học: Các ví dụ số trình bày chi tiết về hình học mô hình, điều kiện biên, thông số vật liệu (ví dụ: mô đun Young, tỷ lệ Poisson, đặc tính chỉnh hướng, gradient vật liệu trong FGM) và điều kiện tải trọng (trang 44, 57, 85, 105, 136).
• Kiểm tra hội tụ và độ tin cậy: Các nghiên cứu hội tụ với các mật độ nút khác nhau (ví dụ: Hình 4.6, trang 65), ảnh hưởng của kích thước miền hỗ trợ (Bảng 3.2, trang xvi) được trình bày, cho phép các nhà nghiên cứu khác kiểm tra tính đúng đắn của việc triển khai của họ. Các thông tin này, cùng với các so sánh rộng rãi với các giải pháp phân tích, thực nghiệm và số đã được công bố, tạo thành một khung nhân rộng mạnh mẽ.

5. 10-year research agenda outlined? Luận án phác thảo một chương trình nghiên cứu 10 năm tiềm năng, mặc dù không được trình bày dưới dạng một dòng thời gian cụ thể. Thay vào đó, phần "Limitations and Future Research" (trang 156-157) cung cấp 4-5 hướng nghiên cứu cụ thể, đủ để tạo thành một chương trình nghị sự dài hạn cho các nghiên cứu tiếp theo:

1. Mở rộng sang các bài toán nứt 3D: Đây là một bước tiến tự nhiên và phức tạp, đòi hỏi nhiều năm nghiên cứu để phát triển công thức, thuật toán và giải quyết thách thức tính toán (trang 157).
2. Kết hợp cơ học phá hủy đàn dẻo phi tuyến: Xử lý các vật liệu có độ dẻo đáng kể là một lĩnh vực rộng lớn, bao gồm phát triển các mô hình cấu thành vật liệu phức tạp và tiêu chí truyền nứt phi tuyến (trang 157).
3. Tối ưu hóa tham số tự động và phát triển thuật toán song song: Hai lĩnh vực này có thể được phát triển song song trong dài hạn để tăng cường tính mạnh mẽ và hiệu quả của phương pháp (trang 157). Tối ưu hóa tự động có thể bao gồm việc sử dụng thuật toán học máy.
4. Ứng dụng cho các hiện tượng phá hủy phức tạp: Mở rộng sang các kịch bản như nứt mỏi, nứt ăn mòn, hoặc tương tác vật liệu-môi trường đòi hỏi các mô hình lý thuyết và tính toán chuyên biệt (trang 157).
5. Phân tích định lượng sự không chắc chắn (Uncertainty Quantification): Đây là một xu hướng ngày càng tăng trong kỹ thuật để đánh giá tác động của các biến động trong dữ liệu đầu vào lên kết quả mô phỏng, có thể phát triển thành một dòng nghiên cứu độc lập. Những hướng này cung cấp một tầm nhìn rõ ràng và đầy tham vọng cho các nghiên cứu tiếp theo, phản ánh cam kết của tác giả trong việc thúc đẩy lĩnh vực cơ học phá hủy tính toán trong dài hạn.

Kết luận

Luận án này đã tạo ra một dấu ấn quan trọng trong lĩnh vực cơ học phá hủy tính toán, mang lại những đóng góp cụ thể và sâu sắc.

1. Phát triển và xác nhận các phương pháp X-RPIM và X-MK tiên tiến: Nghiên cứu đã thành công trong việc phát triển và xác nhận tính đúng đắn của các phương pháp không lưới Galerkin mở rộng (X-RPIM và X-MK) cho phân tích nứt gãy 2D, bao gồm cả các bài toán nứt tĩnh và động, trong nhiều loại vật liệu khác nhau.
2. Giới thiệu và xác nhận hàm dốc tuyến tính hiệu quả: Việc tích hợp lần đầu tiên hàm dốc tuyến tính cho làm giàu đầu nứt vào khung không lưới X-RPIM đã được chứng minh là một đổi mới hiệu quả, giảm đáng kể số bậc tự do bổ sung (DOFs) từ 8 xuống chỉ còn 2 cho mỗi nút đầu nứt, trong khi vẫn duy trì độ chính xác cao so với các hàm nhánh tiêu chuẩn và giải pháp chính xác.
3. Cải tiến đột phá phương pháp Moving Kriging: Luận án đã đóng góp ba loại hàm tương quan mới để cải thiện hàm hình dạng Moving Kriging (MK), trong đó hai loại được trình bày lần đầu tiên, giúp loại bỏ ảnh hưởng của tham số kinh nghiệm người dùng và tăng cường độ mạnh mẽ của phương pháp.
4. Khung phân tích đa vật liệu toàn diện: Nghiên cứu đã cung cấp một khung tính toán thống nhất có khả năng phân tích hành vi nứt gãy trong vật liệu đẳng hướng, vật liệu composite chỉnh hướng và vật liệu có gradient chức năng (FGM), giải quyết sự phức tạp của vật liệu tiên tiến trong các ứng dụng kỹ thuật thực tế.
5. Độ mạnh mẽ và khả năng áp dụng vượt trội: Các phương pháp được phát triển đã được xác nhận rộng rãi thông qua so sánh với các giải pháp phân tích, thực nghiệm và các phương pháp số đã được thiết lập (FEM, XFEM, BEM, EFG, SBFEM), chứng minh độ chính xác, hiệu quả và khả năng áp dụng của chúng cho nhiều hình học và điều kiện tải trọng phức tạp.
6. Giải quyết thách thức chia lại lưới: Bằng cách sử dụng phương pháp không lưới kết hợp với kỹ thuật làm giàu và Level Set vector, luận án đã thành công trong việc tránh nhu cầu chia lại lưới phức tạp trong quá trình lan truyền nứt, một lợi thế tính toán quan trọng so với FEM truyền thống.

Sự tiến bộ mô hình (paradigm advancement) thể hiện trong luận án này không phải là sự thay đổi hoàn toàn các nguyên lý cơ bản của cơ học phá hủy mà là sự chuyển đổi đáng kể trong cách tiếp cận các vấn đề nứt gãy phức tạp về mặt tính toán. Bằng chứng từ các phát hiện cho thấy rằng các phương pháp không lưới, đặc biệt là khi được tối ưu hóa với các kỹ thuật làm giàu mới và các hàm tương quan cải tiến, có thể cung cấp các giải pháp mạnh mẽ và hiệu quả hơn nhiều so với các phương pháp truyền thống. Việc giảm DOFs và loại bỏ sự phụ thuộc vào tham số người dùng là những bước tiến quan trọng hướng tới các công cụ mô phỏng tự động và đáng tin cậy hơn, thách thức các giả định trước đây về sự đánh đổi giữa độ chính xác và hiệu quả tính toán.

Công trình này đã mở ra ít nhất ba dòng nghiên cứu mới:

1. Phát triển và ứng dụng các phương pháp không lưới với làm giàu tiên tiến cho các bài toán nứt gãy 3D và các mô hình vật liệu phi tuyến phức tạp.
2. Nghiên cứu sâu hơn về tối ưu hóa tham số tự động và việc tích hợp các phương pháp học máy để nâng cao độ mạnh mẽ và hiệu quả của các phương pháp không lưới.
3. Mở rộng các ứng dụng của X-RPIM và X-MK cho các hiện tượng phá hủy vật lý khác ngoài nứt gãy, như sự phân hủy giao diện hoặc mô hình hóa các lỗi khác trong vật liệu composite.

Tính liên quan toàn cầu của luận án được thể hiện qua việc so sánh các kết quả với các nghiên cứu quốc tế từ các nhóm nghiên cứu hàng đầu (ví dụ: XFEM [18], BEM [100], FEM [105], SBFEM [126]). Các vấn đề nứt gãy trong vật liệu tiên tiến là một thách thức toàn cầu trong kỹ thuật, và các giải pháp được trình bày ở đây cung cấp một đóng góp đáng kể cho cộng đồng khoa học quốc tế. Di sản của công trình này có thể đo lường được thông qua tiềm năng trích dẫn học thuật, việc áp dụng trong các công cụ phần mềm kỹ thuật, và sự cải thiện trong thiết kế và an toàn của các cấu trúc quan trọng trên toàn thế giới, đặc biệt là trong các ngành công nghiệp hàng không vũ trụ, năng lượng và xây dựng.