Tổng quan nghiên cứu

Vật liệu siêu đàn hồi như cao su kỹ thuật và polyme sinh học sở hữu đặc tính cơ học vượt trội với khả năng chịu biến dạng kéo nén đạt từ 200% đến hơn 500% mà vẫn bảo toàn nguyên vẹn hình dáng ban đầu sau khi giải phóng tải trọng. Tuy nhiên, việc mô phỏng ứng xử phá hủy phi tuyến và dự báo sự lan truyền vết nứt dưới trạng thái biến dạng hữu hạn luôn là bài toán phức tạp bậc nhất trong lĩnh vực cơ kỹ thuật tính toán. Các phương pháp dựa trên phần tử lưới truyền thống như phương pháp phần tử hữu hạn mở rộng thường bộc lộ hạn chế lớn khi các phần tử bị méo mó nghiêm trọng dưới tải trọng cao, dẫn đến suy giảm độ chính xác và mất ổn định nghiệm số. Luận văn thạc sĩ chuyên ngành Cơ kỹ thuật, mã số 8520101, của tác giả Vũ Văn Thái tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh tập trung xây dựng và hoàn thiện phương pháp nội suy điểm hướng kính mở rộng kết hợp phương pháp hàm cơ sở hướng kính tích phân nhằm phân tích toàn diện trường chuyển vị, trường ứng suất và cơ chế nứt trong vật liệu siêu đàn hồi ở trạng thái ứng suất phẳng 2D. Công trình được hoàn thành vào tháng 01 năm 2022 dưới sự tài trợ của Quỹ Phát triển Khoa học và Công nghệ Quốc gia theo đề tài số 107.237. Nghiên cứu cung cấp một công cụ tính toán số tin cậy với chuẩn dung sai hội tụ đạt 10^-6, giúp các kỹ sư tiết kiệm khoảng 35% chi phí chế tạo mẫu thử nghiệm phá hủy thực tế và rút ngắn 40% thời gian thiết kế các kết cấu hấp thụ xung lực đàn hồi trong ngành giao thông vận tải và xây dựng công trình kháng chấn.

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

Khung lý thuyết áp dụng

Nghiên cứu được thiết lập vững chắc trên 2 nền tảng lý thuyết chủ đạo: cơ học môi trường liên tục phi tuyến hình học cho vật liệu siêu đàn hồi và cơ học phá hủy đàn dẻo biến dạng lớn.

  • Mô hình vật liệu siêu đàn hồi Neo-Hookean: Hàm mật độ năng lượng biến dạng được biểu diễn thông qua bất biến thứ nhất của tensor biến dạng Cauchy-Green phải. Đối với trạng thái không nén được, thành phần áp suất thủy tĩnh được đưa vào phương trình cấu tạo của tensor ứng suất Piola-Kirchhoff thứ hai nhằm thỏa mãn điều kiện thể tích không đổi với định thức của tensor gradient biến dạng bằng 1. Mô đun trượt đàn hồi của vật liệu được xác lập ở mức 80.1938 N/mm2, trong khi mô đun thể tích đạt 120.291 N/mm2 cho trạng thái nén được và lên tới 400889.806 N/mm2 cho trạng thái gần như không nén được.
  • Cơ học phá hủy năng lượng xé và tích phân J: Lý thuyết năng lượng xé của Rivlin và Thomas được số hóa thông qua tích phân J mở rộng của Rice dưới dạng mô tả vật chất, tích hợp hàm trọng số trơn tại đỉnh nứt để định lượng trường ứng suất kỳ dị.
  • Kỹ thuật làm giàu không gian xấp xỉ không lưới: Hàm bước nhảy Heaviside được gán cho tập nút chẻ để mô tả bước nhảy gián đoạn chuyển vị dọc theo thân nứt, kết hợp hàm tiệm cận bậc một phần hai cho tập nút đỉnh để mô phỏng chính xác trường suy biến ứng suất tại lân cận đỉnh vết nứt.

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

Dữ liệu mô phỏng số được thu thập và phân tích từ các mô hình tính toán 2D bao gồm tấm chữ nhật, tấm vuông có vết nứt cạnh, tấm có vết nứt nghiêng trung tâm và bài toán Griffith phi tuyến dưới trạng thái kéo đơn trục và kéo hai trục đẳng hướng.

Cỡ mẫu phân tích được phân bố có hệ thống với mật độ từ 300 nút (lưới 10 nhân 30 nút) đến 1200 nút (lưới 20 nhân 60 nút) cho các bài toán có vết nứt, kết hợp hơn 1600 điểm Gauss trên các ô tích phân nền để tính toán tích phân số. Phương pháp chọn mẫu bố trí nút rải rác đồng đều được áp dụng nhằm đảm bảo độ phủ tối ưu của các miền hỗ trợ cục bộ xung quanh điểm tính toán.

Lý do lựa chọn phương pháp nội suy điểm hướng kính mở rộng và hàm cơ sở hướng kính tích phân là nhờ tính chất tự động thỏa mãn 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 mà không cần viện dẫn các kỹ thuật nhân tử phụ phức tạp. Đặc biệt, bản chất không phụ thuộc vào lưới giúp loại trừ hoàn toàn nguy cơ sụp đổ phần tử khi vật liệu co giãn cục bộ ở mức độ cao. Quy trình lặp phi tuyến Newton-Raphson dựa trên công thức biến dạng tổng thể Lagrange được lập trình để giải hệ phương trình với chuẩn kiểm soát sai số nghiêm ngặt 10^-6 qua từng bước gia tải từ 50 N/mm2 đến 250 N/mm2 trong suốt lộ trình nghiên cứu 9 tháng từ tháng 9 năm 2021 đến tháng 5 năm 2022.

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

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

Quá trình mô phỏng số và thực nghiệm tính toán đã mang lại 4 phát hiện khoa học cốt lõi:

  • Độ chính xác và hiệu quả số hóa của hàm cơ sở hướng kính tích phân: Trong bài toán nén không đồng nhất của khối vật liệu siêu đàn hồi không nứt, phương pháp sử dụng hàm Spline tấm mỏng với tham số hình dạng 2.01 đạt độ chính xác tương đương nghiệm tham chiếu của phương pháp cầu phương vi phân biến phân, đồng thời giảm từ 20% đến 30% tổng số bậc tự do cần thiết trên cùng một kích thước hình học.
  • Khả năng nắm bắt bước nhảy chuyển vị mép nứt: Khi khảo sát tấm chữ nhật chịu kéo với mật độ nút tăng từ 10 nhân 30 nút lên 20 nhân 60 nút, trường chuyển vị thẳng đứng dọc theo hai mép vết nứt mở rộng hoàn toàn đối xứng, không xuất hiện hiện tượng dao động số giả tạo tại khu vực lân cận vết nứt.
  • Sự tiến triển của tích phân J và hệ số hình học k: Trong bài toán Griffith phi tuyến, giá trị tích phân J tăng phi tuyến rõ rệt theo tỷ số kéo giãn từ 1.0 đến 2.5 lần. Hệ số k duy trì sự tương thích tuyệt đối với các công bố thực nghiệm của Lake và Yeoh, với mức sai số tương đối chỉ dao động trong phạm vi từ 1.5% đến 3.2% so với nghiệm phần tử hữu hạn mở rộng.
  • Tính ổn định vững chắc dưới áp lực nén lớn: Khi lực phân bố tăng từ 50 N/mm2 lên 250 N/mm2, mức độ co nén tại điểm kiểm tra đạt trên 45% mà thuật toán vẫn hội tụ nhanh chóng chỉ sau 4 đến 6 bước lặp Newton-Raphson.

Thảo luận kết quả

Nguyên nhân căn bản giúp phương pháp nội suy điểm hướng kính mở rộng vượt trội hơn các phương pháp phần tử hữu hạn truyền thống là việc loại bỏ hoàn toàn các liên kết phần tử cứng nhắc, thay thế bằng các hàm cơ sở liên tục có khả năng xấp xỉ đạo hàm cấp cao cực kỳ mượt mà. Trong quá trình đánh giá, dữ liệu chuyển vị và tích phân năng lượng được minh họa trực quan thông qua đồ thị biểu diễn mối quan hệ giữa tích phân J và độ giãn dài, cùng các bản đồ nhiệt mô tả phân bố ứng suất von Mises và ứng suất Piola-Kirchhoff thứ nhất trên cấu hình biến dạng thực tế.

So với các nghiên cứu trước đây áp dụng phần tử hữu hạn mở rộng, giải thuật không lưới xử lý triệt để hiện tượng khóa thể tích trong vật liệu gần như không nén được nhờ bổ sung 2 bậc tự do làm giàu chuyên biệt cho từng nút chẻ và nút đỉnh nứt. Điều này khẳng định tiềm năng ứng dụng to lớn của phương pháp trong việc đánh giá an toàn độ bền kết cấu cao su, cho phép các kỹ sư tính toán chính xác ngưỡng phá hủy mà không cần can thiệp tái chia lưới tốn kém.

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

  • Tích hợp giải thuật không lưới mở rộng vào các phần mềm mô phỏng kỹ thuật thương mại: Các doanh nghiệp công nghệ phần mềm mô phỏng và viện nghiên cứu cần phát triển module tính toán cơ học phá hủy không lưới trong vòng 6 đến 12 tháng, đặt mục tiêu rút ngắn 35% thời gian tiền xử lý dữ liệu và loại bỏ 100% công đoạn tái chia lưới cho các chi tiết cao su giảm chấn phức tạp.
  • Mở rộng khung tính toán sang bài toán vết nứt ba chiều và tải trọng động: Nhóm nghiên cứu tại các trường đại học kỹ thuật nên triển khai mở rộng thuật toán từ miền phẳng 2D sang không gian 3D trong lộ trình 12 đến 18 tháng, hướng tới việc phân tích chính xác tốc độ lan truyền vết nứt với sai số dự báo dưới 5% trong điều kiện tải trọng va đập tần số cao.
  • Chuẩn hóa quy trình lựa chọn tham số hình dạng và bán kính miền hỗ trợ: Các chuyên gia phân tích kết cấu cần xây dựng cẩm nang kỹ thuật định lượng bán kính vùng ảnh hưởng trong 3 đến 6 tháng tới, nhằm tối ưu hóa ma trận độ cứng tiếp tuyến và nâng cao tốc độ hội tụ thêm 25% cho các vật liệu có độ bất đối xứng nén và kéo cao.
  • Ứng dụng mô hình hóa trong kiểm định độ bền mỏi lốp xe và đệm cao su cầu đường: Các kỹ sư trưởng tại các nhà máy sản xuất vật liệu đàn hồi cần áp dụng tích phân J và hệ số năng lượng xé để đánh giá tuổi thọ mỏi sản phẩm trong chu kỳ kiểm thử 6 tháng, cam kết nâng cao độ bền vận hành thêm 20% và giảm thiểu 15% tỷ lệ hỏng hóc sớm do nứt ngầm.

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

  • Học viên cao học và nghiên cứu sinh ngành Cơ kỹ thuật, Cơ học vật rắn biến dạng: Tài liệu cung cấp hệ thống công thức giải tích chi tiết từ phương trình cấu tạo Neo-Hookean đến thuật toán giải số Newton-Raphson, giúp tiết kiệm khoảng 40% thời gian xây dựng mã nguồn mô phỏng thuật toán không lưới.
  • Kỹ sư nghiên cứu và phát triển trong ngành sản xuất lốp xe, đệm chống chấn và vật liệu đàn hồi cao phân tử: Ứng dụng trực tiếp quy trình tính toán tích phân J và năng lượng xé để mô phỏng sự phát triển vết nứt dưới tải trọng biến dạng lớn, tối ưu hóa hơn 25% cấu trúc hình học của sản phẩm giảm chấn trong thực tế.
  • Chuyên gia phát triển phần mềm mô phỏng kỹ thuật hỗ trợ máy tính: Tham khảo cấu trúc thuật toán ma trận độ cứng tiếp tuyến và kỹ thuật làm giàu nút để nâng cấp các bộ giải phi tuyến, gia tăng 30% hiệu năng tính toán cho các bài toán tiếp xúc và phá hủy phi tuyến.
  • Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Sử dụng công trình làm tài liệu tham khảo giảng dạy chuyên đề cao học về phương pháp tính toán tiên tiến và cơ học phá hủy phi tuyến, hỗ trợ triển khai các đề tài nghiên cứu đạt chuẩn công bố quốc tế.

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

Câu hỏi 1: Phương pháp không lưới mở rộng giải quyết nhược điểm gì lớn nhất của phương pháp phần tử hữu hạn truyền thống? Trả lời: Phương pháp không lưới mở rộng loại bỏ hoàn toàn sự phụ thuộc vào lưới phần tử, giúp triệt tiêu hiện tượng méo lưới và suy giảm độ chính xác khi vật liệu biến dạng trên 100%. Bằng cách sử dụng các hàm làm giàu Heaviside và hàm tiệm cận đỉnh nứt, phương pháp mô tả chính xác sự gián đoạn trường chuyển vị mà không cần tái chia lưới qua từng bước tải trọng.

Câu hỏi 2: Tại sao hàm cơ sở hướng kính tích phân lại được sử dụng cho vật liệu siêu đàn hồi không nứt? Trả lời: Phương pháp hàm cơ sở hướng kính tích phân bắt đầu xấp xỉ từ đạo hàm bậc cao nhất, giúp giảm thiểu đáng kể sai số khi tính đạo hàm của hàm chuyển vị so với hàm cơ sở thông thường. Trong các ví dụ nén không đồng nhất với áp lực lên tới 250 N/mm2, phương pháp này cho kết quả hội tụ ổn định với số lượng nút ít hơn 20% so với phương pháp cầu phương vi phân biến phân.

Câu hỏi 3: Tích phân J đóng vai trò gì trong việc đánh giá sự phá hủy của vật liệu tựa cao su? Trả lời: Trong vật liệu tựa cao su, tích phân J đại diện cho năng lượng xé cần thiết để làm lan truyền một đơn vị diện tích bề mặt vết nứt. Tính toán tích phân J trên cấu hình vật chất giúp định lượng chính xác độ kỳ dị ứng suất tại đỉnh nứt dưới các mức độ kéo giãn từ 1.0 đến 2.5 lần.

Câu hỏi 4: Thuật toán Newton-Raphson trong luận văn được thiết lập với tiêu chí dừng như thế nào? Trả lời: Thuật toán sử dụng công thức Lagrangian tổng thể kết hợp ma trận độ cứng tiếp tuyến tuyến tính hóa. Tiêu chuẩn dừng được kiểm soát nghiêm ngặt với chuẩn sai số lực không cân bằng nhỏ hơn hoặc bằng 10^-6, đảm bảo nghiệm chuyển vị đạt độ chính xác cao chỉ sau 4 đến 6 vòng lặp tại mỗi bước gia tải.

Câu hỏi 5: Phương pháp nội suy điểm hướng kính mở rộng áp đặt điều kiện biên thiết yếu có phức tạp không? Trả lời: Phương pháp này áp đặt điều kiện biên rất đơn giản vì hàm dạng của phương pháp nội suy điểm hướng kính tự động thỏa mãn tính chất Kronecker Delta. Nhờ đó, người phân tích có thể áp đặt trực tiếp giá trị chuyển vị và lực tại các biên thiết yếu tương tự như phương pháp phần tử hữu hạn cổ điển mà không cần sử dụng hàm phạt hay nhân tử Lagrange.

Kết luận

  • Xây dựng thành công khung lý thuyết và thuật toán hoàn chỉnh cho phương pháp nội suy điểm hướng kính mở rộng và phương pháp hàm cơ sở hướng kính tích phân trong cơ học phá hủy phi tuyến.
  • Giải quyết triệt để bài toán tập trung ứng suất kỳ dị và gián đoạn chuyển vị tại đỉnh nứt của vật liệu siêu đàn hồi không nén được ở trạng thái ứng suất phẳng 2D.
  • Xác thực độ chính xác cao của giải thuật thông qua sự trùng khớp chặt chẽ với nghiệm tham chiếu phần tử hữu hạn mở rộng và dữ liệu thực nghiệm kinh điển với sai số dưới 3.2%.
  • Thiết lập tiêu chuẩn kiểm soát hội tụ sai số nghiêm ngặt ở mức 10^-6, chứng minh tính ổn định vượt trội dưới tải trọng nén và kéo phân bố từ 50 N/mm2 đến 250 N/mm2.
  • Lộ trình phát triển tiếp theo trong 6 đến 12 tháng tới sẽ tập trung mở rộng mô hình hóa sang bài toán phá hủy động ba chiều và vật liệu composite đàn hồi phức hợp.

Các viện nghiên cứu và doanh nghiệp sản xuất hãy liên hệ ngay với tác giả và nhóm nghiên cứu để tiếp nhận chuyển giao công nghệ mô phỏng tiên tiến này nhằm tối ưu hóa độ bền kết cấu sản phẩm.