Nghiên cứu phương pháp meshfree mở rộng cho vật liệu siêu đàn hồi có vết nứt

Khám phá phương pháp meshfree mở rộng cho vật liệu siêu đàn hồi có vết nứt trong luận văn thạc sĩ cơ kỹ thuật. Tìm hiểu ứng dụng và hiệu quả.

• Tác giả: Vũ Văn Thái
• Người hướng dẫn: TS. Nguyễn Thanh Nhã
• Trường học: Trường Đại Học Bách Khoa - ĐHQG - HCM
• Chuyên ngành: Cơ Kỹ Thuật
• Đề tài: Phương Pháp Không Lưới Mở Rộng Cho Vật Liệu Siêu Đàn Hồi Có Vết Nứt
• Loại tài liệu: Luận Văn Thạc Sĩ
• Năm xuất bản: 2022
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI MỞ ĐẦU
• 1. CHƯƠNG 1: INTRODUCTION
• 1.1. State of the art
• 1.2. Scope of the study
• 1.3. Research objectives
• 1.4. Author’s contributions
• 1.5. Thesis outline
• 2. CHƯƠNG 2: METHODOLOGY
• 2.1. Constitutive equations of hyperelastic material
• 2.2. Fracture analysis of hyperelastic material
• 2.3. Meshfree shape functions construction
• 2.3.1. Radial Point Interpolation Method (RPIM)
• 2.3.2. Integrated Radial Basis Functions Method (iRBF)
• 2.4. The XRPIM for crack problem in hyperelastic bodies
• 2.4.1. Enriched approximation of the displacement field by XRPIM
• 2.4.2. Weak form for nonlinear elastic problem and discrete equations
• 2.4.2.1. Numerical implementation procedure
• 2.4.2.2. Computation procedure of K matrix and fint matrix
• 2.4.2.3. Computation procedure of B matrix and O matrix
• 3. CHƯƠNG 3: NUMERICAL EXAMPLES
• 3.1. Non-cracked hyperelastic solid
• 3.1.1. Inhomogeneous compression problem
• 3.1.2. Curved beam problem
• 3.2. Cracked hyperelastic solid
• 3.2.1. Rectangular plate with an edge crack under prescribed extension
• 3.2.2. Square plate with an edge crack under prescribed extension
• 3.2.3. Nonlinear Griffith problem
• 3.2.4. Square plate with an inclined central crack
• 4. CHƯƠNG 4: IMPLEMENTATION
• 5. CHƯƠNG 5: CONCLUSION AND OUTLOOK
• LIST OF PUBLICATIONS
• REFERENCES

Tóm tắt

I. Giới thiệu về phương pháp meshfree

Phương pháp meshfree đã trở thành một công cụ quan trọng trong phân tích các vấn đề liên quan đến vật liệu siêu đàn hồi. Đặc điểm nổi bật của phương pháp này là khả năng loại bỏ sự phụ thuộc vào lưới, cho phép các nút trong miền tính toán có thể được thêm, bớt hoặc thay đổi vị trí một cách linh hoạt. Điều này đặc biệt hữu ích trong các bài toán mà hình dạng miền thay đổi liên tục. Việc áp dụng phương pháp meshfree giúp tăng độ chính xác trong việc mô hình hóa các hiện tượng bất liên tục và các trường ứng suất xung quanh đỉnh vết nứt trong vật liệu siêu đàn hồi. Theo nghiên cứu của Mai et al., việc sử dụng các hàm cơ sở hướng kính tích hợp (iRBF) có thể cải thiện độ chính xác cho các phép xấp xỉ đạo hàm của một hàm số.

1.1. Tính toán trường ứng suất và chuyển vị

Việc tính toán trường ứng suất và chuyển vị trong vật liệu siêu đàn hồi là một phần quan trọng trong nghiên cứu này. Các phương pháp như phân tích cấu trúc và mô hình hóa hóa vật liệu được áp dụng để xây dựng các phương trình mô tả mối quan hệ giữa ứng suất và biến dạng. Phương pháp meshfree cho phép mô hình hóa các trường ứng suất mà không cần lưới, từ đó giảm thiểu sai số do biến dạng lớn trong quá trình tính toán. Phương pháp này cũng cho phép tính toán hiệu quả hơn trong các bài toán có vết nứt, nơi mà các ứng suất có thể thay đổi đột ngột.

II. Phân tích vết nứt trong vật liệu siêu đàn hồi

Phân tích vết nứt trong vật liệu siêu đàn hồi là một thách thức lớn trong cơ học tính toán. Các vết nứt có thể gây ra sự suy giảm đáng kể về khả năng chịu lực của vật liệu. Việc áp dụng phương pháp meshfree cho phép mô hình hóa chính xác hơn các hiện tượng này. Đặc biệt, phương pháp XRPIM (Extended Radial Point Interpolation Method) đã được chứng minh là hiệu quả trong việc phân tích sự phát triển của vết nứt trong các vật liệu này. Nghiên cứu cho thấy rằng việc sử dụng hàm Heaviside và hàm nhánh trong XRPIM giúp mô tả chính xác hơn các trường ứng suất xung quanh đỉnh vết nứt.

2.1. Mô hình hóa và giải pháp phân tích

Mô hình hóa vết nứt trong vật liệu siêu đàn hồi yêu cầu một phương pháp phân tích đặc biệt. Phương pháp XRPIM kết hợp giữa phương pháp nội suy điểm hướng kính và các hàm làm giàu để mô tả các bất liên tục trong trường chuyển vị. Phương pháp này không chỉ giúp tính toán chính xác trường ứng suất mà còn cho phép theo dõi sự phát triển của vết nứt một cách hiệu quả. Nghiên cứu đã chỉ ra rằng các kết quả thu được từ phương pháp XRPIM có độ chính xác cao khi so sánh với các giải pháp tham khảo, cho thấy giá trị thực tiễn của phương pháp này trong phân tích vật liệu siêu đàn hồi có vết nứt.

III. Ứng dụng thực tiễn của phương pháp meshfree

Phương pháp meshfree mở rộng có nhiều ứng dụng thực tiễn trong ngành công nghiệp, đặc biệt là trong lĩnh vực vật liệu composite và các cấu trúc chịu tải trọng lớn. Việc áp dụng phương pháp này giúp cải thiện độ tin cậy của các mô hình phân tích, từ đó tối ưu hóa thiết kế và giảm thiểu rủi ro trong các ứng dụng thực tế. Các nghiên cứu cho thấy rằng phương pháp này có thể được áp dụng hiệu quả trong việc phân tích các cấu trúc phức tạp, nơi mà các phương pháp truyền thống gặp khó khăn trong việc xử lý. Ngoài ra, phương pháp cũng có thể được mở rộng để áp dụng cho các loại vật liệu khác nhau, bao gồm cả vật liệu thông minh.

3.1. Tương lai của nghiên cứu và phát triển

Tương lai của nghiên cứu phương pháp meshfree hứa hẹn sẽ mang lại nhiều tiến bộ trong lĩnh vực cơ học tính toán. Các nghiên cứu tiếp theo có thể tập trung vào việc cải thiện độ chính xác của các mô hình, cũng như mở rộng khả năng ứng dụng của phương pháp này trong các lĩnh vực khác nhau. Việc tích hợp các công nghệ mới như trí tuệ nhân tạo và học máy vào phương pháp meshfree có thể mở ra những hướng đi mới trong nghiên cứu và ứng dụng thực tiễn, giúp giải quyết các bài toán phức tạp một cách hiệu quả hơn.

Nội dung chính

Tổng quan nghiên cứu

Vật liệu siêu đàn hồi, điển hình là các vật liệu giống cao su, được ứng dụng rộng rãi trong nhiều lĩnh vực kỹ thuật nhờ khả năng biến dạng lớn dưới tải trọng nhỏ và khả năng phục hồi hình dạng ban đầu sau khi gỡ tải. Theo ước tính, các vật liệu này có thể chịu biến dạng lớn gấp nhiều lần so với vật liệu đàn hồi thông thường, đồng thời giữ được tính chất không nén được hoặc gần như không nén được. Tuy nhiên, việc mô phỏng sự phá hủy biến dạng lớn trong vật liệu siêu đàn hồi vẫn là một thách thức lớn trong cơ học tính toán do tính phi tuyến và sự phức tạp của trường ứng suất, đặc biệt khi xuất hiện vết nứt.

Luận văn tập trung nghiên cứu phương pháp không lưới mở rộng (Extended Radial Point Interpolation Method - XRPIM) để giải quyết bài toán nứt trong vật liệu siêu đàn hồi ở trạng thái không nén được, trong điều kiện ứng suất mặt phẳng. Mục tiêu chính là xây dựng và phát triển phương pháp số nhằm mô phỏng chính xác trường chuyển vị, ứng suất, tích phân J và hệ số k đặc trưng cho sự phá hủy của vật liệu. Ngoài ra, phương pháp tích hợp hàm cơ sở hướng kính (integrated Radial Basis Functions - iRBF) cũng được áp dụng để phân tích vật liệu siêu đàn hồi không có vết nứt ở cả trạng thái nén được và gần như không nén được.

Phạm vi nghiên cứu bao gồm các mô hình 2D, với các ví dụ số minh họa như tấm chữ nhật có vết nứt cạnh, tấm vuông có vết nứt trung tâm nghiêng, và các bài toán nén không đồng nhất. Kết quả thu được được so sánh với các lời giải tham khảo từ các phương pháp tiên tiến như XFEM, cho thấy độ chính xác và hiệu quả của phương pháp đề xuất. Nghiên cứu có ý nghĩa quan trọng trong việc phát triển các công cụ tính toán không lưới cho bài toán phá hủy phi tuyến, góp phần nâng cao khả năng thiết kế và đánh giá độ bền của các cấu kiện sử dụng vật liệu siêu đàn hồi.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn dựa trên các lý thuyết và mô hình sau:

• Mô hình vật liệu siêu đàn hồi Neo-Hookean không nén được: Mô hình này sử dụng hàm năng lượng biến dạng $\psi = \frac{\mu}{2}(\mathrm{tr}(C) - 3)$, trong đó $C$ là tensor biến dạng Cauchy-Green phải, $\mu$ là mô đun cắt. Căng thẳng được xác định qua đạo hàm bậc nhất của hàm năng lượng theo biến dạng, đồng thời áp dụng điều kiện không nén được thông qua áp suất thủy tĩnh $p$.

• Phân tích phá hủy phi tuyến dựa trên tích phân J: Tích phân J được sử dụng để đánh giá năng lượng xé và sự tập trung ứng suất quanh đầu vết nứt, mở rộng từ cơ học phá hủy tuyến tính sang phi tuyến. Công thức J được biểu diễn trong hệ tọa độ vật liệu, tích phân trên đường bao quanh đầu nứt.

• Phương pháp không lưới Radial Point Interpolation Method (RPIM): RPIM sử dụng các hàm cơ sở hướng kính để xây dựng hàm dạng (shape functions) không phụ thuộc vào lưới phần tử, giúp mô phỏng biến dạng lớn và sự thay đổi hình học phức tạp.

• Phương pháp tích hợp hàm cơ sở hướng kính (iRBF): iRBF bắt đầu từ đạo hàm bậc cao của hàm cần xấp xỉ, sau đó tích phân để thu được hàm gốc, giúp cải thiện độ chính xác của đạo hàm và giảm sai số trong mô phỏng biến dạng lớn.

• Phương pháp mở rộng XRPIM: Kết hợp RPIM với các hàm làm giàu (enrichment functions) như hàm Heaviside để mô tả sự gián đoạn chuyển vị tại vết nứt và hàm branch để mô tả trường ứng suất có điểm kỳ dị tại đầu vết nứt, tương tự như trong Extended Finite Element Method (XFEM).

Các khái niệm chính bao gồm: tensor biến dạng Cauchy-Green phải $C$, tích phân J, hệ số k đặc trưng cho trạng thái biến dạng quanh đầu nứt, hàm Heaviside $H(x)$, hàm branch $L(x)$, và các ma trận biến dạng, ứng suất trong công thức Lagrange tổng.

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

Nguồn dữ liệu nghiên cứu chủ yếu là các mô hình số được xây dựng và giải bằng phương pháp không lưới XRPIM và iRBF. Cỡ mẫu gồm các nút phân bố rải rác trong miền tính toán, số lượng nút thay đổi tùy theo bài toán (ví dụ từ 10×30 đến 20×60 nút cho tấm có vết nứt). Phương pháp chọn mẫu là phân bố nút đều hoặc theo yêu cầu độ chính xác tại vùng đầu nứt.

Phân tích được thực hiện theo trình tự:

• Xây dựng hàm dạng RPIM và iRBF, bao gồm tính toán ma trận moment và giải hệ phương trình để xác định hệ số.

• Áp dụng hàm làm giàu cho các nút nằm trên hoặc gần vết nứt để mô tả chính xác sự gián đoạn và điểm kỳ dị.

• Sử dụng công thức Lagrange tổng để mô tả biến dạng lớn, tính toán ma trận độ cứng tiếp tuyến và lực nội tại.

• Tích phân số được thực hiện bằng phương pháp Gaussian Quadrature trên các ô nền (background cells).

• Giải hệ phương trình phi tuyến bằng phương pháp Newton-Raphson với sai số ngưỡng $10^{-6}$.

• Tính toán các đại lượng đặc trưng như trường chuyển vị, ứng suất, tích phân J, hệ số k, và so sánh với các lời giải tham khảo.

Thời gian nghiên cứu kéo dài từ tháng 9/2021 đến tháng 5/2022, thực hiện tại Trường Đại học Bách Khoa – ĐHQG TP.HCM.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

1. Độ chính xác và hội tụ của phương pháp iRBF trong bài toán nén không đồng nhất: Với các lực phân bố $f = 50, 100, 150, 200$ N/mm² cho trạng thái nén được và $f = 100, 150, 200, 250$ N/mm² cho trạng thái gần như không nén được, kết quả biến dạng tại điểm M hội tụ tốt khi tăng số lượng nút. Sai số so với phương pháp tham chiếu Variational Differential Quadrature (VDQ) dưới 2%, trong khi số nút cần thiết ít hơn khoảng 20-30%.

2. Khả năng mô phỏng biến dạng lớn và ứng suất trong vật liệu siêu đàn hồi không nứt: Phương pháp iRBF cho thấy khả năng mô phỏng chính xác trường ứng suất Piola-Kirchhoff đầu tiên và biến dạng lớn, với sai số dưới 3% so với các lời giải tham khảo. Tốc độ hội tụ của trạng thái nén được nhanh hơn trạng thái gần như không nén được.

3. Mô phỏng trường chuyển vị và ứng suất quanh vết nứt bằng XRPIM: Các ví dụ số như tấm chữ nhật có vết nứt cạnh và tấm vuông có vết nứt trung tâm nghiêng cho thấy XRPIM mô phỏng chính xác sự gián đoạn chuyển vị và tập trung ứng suất tại đầu nứt. Giá trị tích phân J và hệ số k biến đổi phù hợp với các kết quả từ Extended Finite Element Method (XFEM), với sai số dưới 5%.

4. Ảnh hưởng của góc nghiêng vết nứt và tải trọng lên tích phân J: Kết quả cho thấy tích phân J tăng theo tải trọng và biến đổi không tuyến tính theo góc nghiêng vết nứt, phản ánh đúng đặc tính phi tuyến của vật liệu siêu đàn hồi. Biểu đồ so sánh kết quả XRPIM với các nghiên cứu trước minh họa sự tương đồng cao.

Thảo luận kết quả

Nguyên nhân của độ chính xác cao đến từ việc XRPIM kết hợp hiệu quả các hàm làm giàu để mô tả chính xác sự gián đoạn và điểm kỳ dị ứng suất, đồng thời loại bỏ nhược điểm méo lưới trong các phương pháp phần tử hữu hạn truyền thống. Việc sử dụng công thức Lagrange tổng và tích phân Gaussian giúp mô phỏng biến dạng lớn một cách ổn định.

So với các nghiên cứu trước, kết quả của luận văn khẳng định tính khả thi và ưu việt của phương pháp không lưới trong phân tích phá hủy phi tuyến của vật liệu siêu đàn hồi. Đặc biệt, phương pháp iRBF cải thiện độ chính xác trong mô phỏng vật liệu không nén được, vốn là thách thức lớn trong cơ học vật liệu.

Dữ liệu có thể được trình bày qua các biểu đồ hội tụ biến dạng theo số nút, biểu đồ biến thiên tích phân J theo tải trọng và góc nghiêng vết nứt, cũng như bảng so sánh giá trị k với các phương pháp khác. Các biểu đồ này minh họa rõ ràng sự ổn định và chính xác của phương pháp đề xuất.

Đề xuất và khuyến nghị

1. Mở rộng ứng dụng XRPIM cho mô phỏng phá hủy vật liệu siêu đàn hồi trong 3D: Đề xuất phát triển thuật toán và chương trình để áp dụng phương pháp cho các bài toán 3 chiều, nhằm nâng cao tính thực tiễn và khả năng ứng dụng trong công nghiệp. Thời gian thực hiện dự kiến 1-2 năm, do các nhóm nghiên cứu cơ học tính toán và phát triển phần mềm đảm nhận.

2. Tích hợp mô hình vật liệu phức tạp hơn như Ogden hoặc Mooney-Rivlin: Để mô phỏng chính xác hơn các vật liệu cao su thực tế, cần xây dựng và kiểm định các mô hình vật liệu phi tuyến đa tham số trong khuôn khổ XRPIM. Mục tiêu cải thiện độ chính xác của dự báo ứng suất và biến dạng, thực hiện trong vòng 12 tháng.

3. Phát triển công cụ tự động hóa phân bố nút và làm giàu hàm số: Đề xuất xây dựng thuật toán tự động chọn vị trí nút và hàm làm giàu phù hợp với hình học vết nứt và biến dạng lớn, nhằm tối ưu hóa hiệu suất tính toán và độ chính xác. Chủ thể thực hiện là nhóm nghiên cứu phần mềm, thời gian 6-9 tháng.

4. Áp dụng phương pháp vào thiết kế và đánh giá độ bền các cấu kiện cao su trong công nghiệp ô tô và xây dựng: Khuyến nghị phối hợp với các doanh nghiệp để thử nghiệm và hiệu chỉnh mô hình trên các cấu kiện thực tế, từ đó hỗ trợ thiết kế tối ưu và dự báo tuổi thọ sản phẩm. Thời gian triển khai 1-3 năm, với sự tham gia của các kỹ sư thiết kế và chuyên gia vật liệu.

Đối tượng nên tham khảo luận văn

1. Nhà nghiên cứu và giảng viên cơ học kỹ thuật: Luận văn cung cấp phương pháp số tiên tiến trong phân tích phá hủy phi tuyến, giúp mở rộng kiến thức và ứng dụng trong giảng dạy, nghiên cứu chuyên sâu về vật liệu siêu đàn hồi và phương pháp không lưới.

2. Kỹ sư thiết kế và phát triển sản phẩm cao su, vật liệu đàn hồi: Các kết quả và phương pháp mô phỏng giúp tối ưu hóa thiết kế sản phẩm chịu biến dạng lớn, giảm chi phí thử nghiệm thực nghiệm và nâng cao độ tin cậy sản phẩm.

3. Chuyên gia phát triển phần mềm mô phỏng cơ học tính toán: Luận văn trình bày chi tiết thuật toán XRPIM và iRBF, cung cấp cơ sở để phát triển hoặc cải tiến các phần mềm mô phỏng không lưới, đặc biệt trong lĩnh vực phá hủy vật liệu phi tuyến.

4. Doanh nghiệp và tổ chức nghiên cứu ứng dụng vật liệu mới: Các kết quả nghiên cứu có thể hỗ trợ đánh giá độ bền, dự báo tuổi thọ và phát triển vật liệu siêu đàn hồi mới, góp phần nâng cao năng lực cạnh tranh và đổi mới công nghệ.

Câu hỏi thường gặp

1. Phương pháp XRPIM khác gì so với XFEM trong mô phỏng vết nứt?
   XRPIM là phương pháp không lưới, không phụ thuộc vào lưới phần tử như XFEM, giúp tránh méo lưới khi biến dạng lớn. XRPIM sử dụng các nút phân bố rải rác và hàm làm giàu tương tự XFEM để mô tả gián đoạn và điểm kỳ dị, nhưng linh hoạt hơn trong xử lý hình học phức tạp.

2. Tại sao chọn mô hình Neo-Hookean cho vật liệu siêu đàn hồi?
   Mô hình Neo-Hookean đơn giản, phù hợp với vật liệu cao su ở biến dạng lớn và trạng thái không nén được. Nó cung cấp công thức năng lượng biến dạng dễ xử lý và được sử dụng rộng rãi trong nghiên cứu cơ bản, làm nền tảng cho các mô hình phức tạp hơn.

3. Làm thế nào để tính tích phân J trong bài toán phi tuyến?
   Tích phân J được tính theo công thức mở rộng trong hệ tọa độ vật liệu, sử dụng tích phân đường bao quanh đầu vết nứt hoặc tích phân thể tích với hàm trọng số. Phương pháp Gaussian Quadrature được áp dụng để tính tích phân số chính xác.

4. Phương pháp iRBF có ưu điểm gì so với RPIM?
   iRBF bắt đầu từ đạo hàm bậc cao của hàm cần xấp xỉ, sau đó tích phân để giảm sai số đạo hàm, cải thiện độ chính xác trong mô phỏng biến dạng lớn. iRBF phù hợp với các bài toán không có vết nứt, trong khi RPIM được mở rộng thành XRPIM để xử lý vết nứt.

5. Phương pháp này có thể áp dụng cho vật liệu nén được không?
   Có, luận văn cũng nghiên cứu vật liệu siêu đàn hồi ở trạng thái nén được và gần như không nén được. Phương pháp iRBF được sử dụng để mô phỏng vật liệu nén được, cho kết quả chính xác và ổn định trong các bài toán biến dạng lớn.

Kết luận

• Phương pháp không lưới mở rộng XRPIM đã được xây dựng và áp dụng thành công cho bài toán nứt trong vật liệu siêu đàn hồi không nén được, mô phỏng chính xác trường chuyển vị, ứng suất và tích phân J.
• Phương pháp iRBF được phát triển để phân tích vật liệu siêu đàn hồi không có vết nứt ở cả trạng thái nén được và gần như không nén được, cho kết quả hội tụ và chính xác cao.
• Các kết quả số được so sánh với các lời giải tham khảo như XFEM và VDQ, cho thấy độ chính xác và hiệu quả vượt trội của phương pháp đề xuất.
• Luận văn góp phần mở rộng ứng dụng của các phương pháp không lưới trong cơ học tính toán, đặc biệt trong phân tích phá hủy phi tuyến của vật liệu siêu đàn hồi.
• Các bước tiếp theo bao gồm mở rộng mô hình sang 3D, tích hợp các mô hình vật liệu phức tạp hơn và phát triển công cụ phần mềm hỗ trợ thiết kế kỹ thuật.

Để tiếp tục nghiên cứu và ứng dụng, độc giả và các nhà nghiên cứu được khuyến khích áp dụng phương pháp XRPIM và iRBF trong các bài toán thực tế, đồng thời phát triển thêm các thuật toán tối ưu hóa và mô hình vật liệu mới.