Tổng quan nghiên cứu

Vật liệu phân lớp chức năng (FGMs) lần đầu tiên được giới thiệu vào năm 1984 nhằm đáp ứng nhu cầu chế tạo các cấu kiện chịu nhiệt độ cao trên 1000 độ C trong ngành hàng không vũ trụ và năng lượng hạt nhân. Điểm ưu việt của vật liệu FGMs là cơ tính biến thiên liên tục theo chiều dày, chẳng hạn mô đun đàn hồi biến thiên từ 70 GPa ở mặt kim loại nhôm lên đến 380 GPa ở mặt gốm oxide nhôm, giúp loại bỏ hoàn toàn hiện tượng tập trung ứng suất và bong tách lớp thường thấy ở vật liệu composite truyền thống. Mặc dù sở hữu nhiều ưu điểm vượt trội, việc mô phỏng chính xác ứng xử cơ học của tấm FGMs vẫn là một thách thức lớn trong cơ học tính toán. Các phần tử tam giác tiêu chuẩn thường gặp hiện tượng khóa cắt nghiêm trọng khi tỷ số chiều dày trên chiều dài tấm giảm xuống dưới 0.01, làm phát sinh sai số tính toán độ võng vượt quá 30% so với thực tế.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng được thực hiện tại Thành phố Hồ Chí Minh vào tháng 4 năm 2019 tập trung giải quyết triệt để rào cản tính toán này. Mục tiêu cốt lõi của đề tài là xây dựng công thức phần tử hữu hạn tam giác 3 nút trơn trên miền cạnh kết hợp kỹ thuật nội suy ten-xơ hỗn hợp cải tiến (ES-MITC3+) dựa trên lý thuyết biến dạng cắt bậc cao của Reddy. Nghiên cứu thực hiện phân tích ứng xử tĩnh đối với kết cấu tấm phẳng và tấm xiên FGMs chịu tải trọng phân bố đều và tải trọng hình sin trong miền đàn hồi tuyến tính. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc nâng cao độ chính xác tính toán với sai số dưới 1.5% so với lời giải giải tích 3D, đồng thời tối ưu hóa khoảng 35% thời gian tạo lưới phần tử hữu hạn cho các kết cấu có biên dạng phức tạp.

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

Khung lý thuyết áp dụng

Khung lý thuyết của nghiên cứu được xây dựng trên sự kết hợp chặt chẽ giữa mô hình vật liệu phân lớp chức năng và các lý thuyết cơ học môi trường liên tục tiên tiến:

  • Quy luật phân bố vật liệu FGMs: Tấm FGMs được cấu tạo từ hai thành phần kim loại và gốm với sự thay đổi liên tục theo phương chiều dày z từ âm h/2 đến dương h/2. Hàm phân phối thể tích pha gốm tuân theo quy luật hàm mũ với chỉ số vật liệu n biến thiên từ 0 (vật liệu gốm thuần nhất) đến vô cùng (vật liệu kim loại thuần nhất). Mô đun đàn hồi E(z) và hệ số Poisson biến thiên tương ứng theo tỷ lệ thể tích.
  • Lý thuyết biến dạng cắt bậc ba của Reddy (HSDT): Mô hình sử dụng trường chuyển vị 7 thành phần độc lập (gồm 3 chuyển vị màng u0, v0, w0, 2 góc xoay quanh trục và 2 hàm độ vênh). Lý thuyết này phản ánh đúng phân bố phi tuyến bậc 2 của ứng suất cắt ngoài mặt phẳng dọc theo chiều dày và tự động thỏa mãn điều kiện ứng suất cắt bằng 0 tại hai bề mặt tấm mà không cần sử dụng hệ số hiệu chỉnh lực cắt nhân tạo như trong lý thuyết bậc nhất (FSDT).
  • Kỹ thuật khử khóa cắt MITC3+: Trường biến dạng cắt ngoài mặt phẳng được xấp xỉ lại dựa trên giá trị biến dạng tính toán tại 6 điểm buộc hình học xác định trên phần tử tam giác có bổ sung nút nổi tại trọng tâm, loại bỏ triệt để hiện tượng khóa cắt khi tấm chuyển sang trạng thái mỏng.
  • Phương pháp phần tử hữu hạn trơn trên cạnh (ES-FEM): Trường biến dạng uốn trong mặt phẳng được làm trơn thông qua phép lấy tích phân đường trên các miền làm trơn quanh cạnh chung của hai phần tử kề nhau, giúp làm mềm độ cứng giả tạo của mô hình phần tử hữu hạn tiêu chuẩn.

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

Nghiên cứu áp dụng quy trình tính toán số hiện đại với các thông số mẫu mô phỏng được chuẩn hóa:

  • Quy mô mẫu lưới và kỹ thuật rời rạc: Kết cấu tấm được rời rạc hóa thành các cấp độ lưới có cấu trúc gồm 8x8x2, 12x12x2, 16x16x2, 20x20x2 và 24x24x2 phần tử tam giác 3 nút, tương ứng với số lượng phần tử dao động từ 128 đến 1152 phần tử và số bậc tự do tối đa lên tới hơn 3500 phương trình đại số tuyến tính.
  • Phương pháp chọn mẫu phân tích: Lựa chọn mô hình lưới tam giác phân chia đều theo phương đường chéo của các ô chữ nhật. Phương pháp này đảm bảo tính đối xứng của trường ứng suất và phản ánh trung thực gradient biến dạng tại các vùng tập trung nội lực.
  • Lý do lựa chọn phương pháp phân tích: Việc kết hợp kỹ thuật ES-FEM và MITC3+ cho phép phần tử tam giác 3 nút đạt được độ chính xác tương đương hoặc vượt trội so với phần tử tứ giác 4 nút bậc cao, trong khi vẫn duy trì sự linh hoạt tối đa khi tạo lưới cho các miền hình học nghiêng hoặc xiên. Toàn bộ thuật toán được lập trình và giải số trên phần mềm MATLAB trong giai đoạn nghiên cứu từ năm 2018 đến tháng 4 năm 2019.

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 và kiểm chứng trên các bài toán chuẩn đã mang lại 4 phát hiện quan trọng:

  • Khắc phục hoàn toàn hiện tượng khóa cắt ở dải tỷ lệ chiều dày rộng: Khi tỷ số chiều dày trên cạnh tấm h/L giảm từ 0.25 (tấm dày) xuống 0.1 (tấm dày trung bình) và đạt mức 0.001 (tấm siêu mỏng), độ võng không thứ nguyên wc của tấm vuông Al/ZrO2-1 tựa đơn giữ ổn định ở mức 0.168. Độ võng tuyệt đối tăng tỷ lệ nghịch với chiều dày mà không hề xảy ra hiện tượng suy giảm độ võng giả tạo do khóa cắt.
  • Tác động rõ nét của chỉ số vật liệu n đến độ võng kết cấu: Đối với tấm vuông Al/ZrO2-1 chịu tải phân bố đều có h/L = 0.2, khi chỉ số n tăng từ 0 (100% gốm, mô đun đàn hồi 200 GPa) lên n = 2 (tăng hàm lượng kim loại nhôm, mô đun đàn hồi 70 GPa), độ võng không thứ nguyên tại tâm tấm tựa đơn tăng từ 0.3093 lên 0.5892, tương ứng với mức tăng hơn 90.5% do độ cứng chống uốn tổng thể của tấm bị suy giảm.
  • Quy luật làm việc của kết cấu tấm xiên FGMs: Khi góc xiên psi của tấm tăng từ 15 độ lên 75 độ đối với tấm tựa đơn (n = 0, h/L = 0.1), độ võng không thứ nguyên wc tại tâm tấm giảm mạnh từ 0.7125 xuống còn 0.0085, tương đương mức giảm tới 98.8%. Điều này chứng minh hình học xiên làm tăng đáng kể độ cứng tổng thể của tấm.
  • Phân bố ứng suất cắt phi tuyến chính xác: Tại bài toán tấm Al/Al2O3 chịu tải hình sin (h/L = 0.1), trường ứng suất tiếp tau-xy và ứng suất cắt tau-xz phân bố hoàn toàn trơn tru dọc theo chiều dày, tự động triệt tiêu về 0 tại các mặt ngoài z = dương h/2 và z = âm h/2, khớp chính xác với kết quả từ lời giải giải tích Quasi-3D với sai số dưới 1.2%.

Thảo luận kết quả

Về mặt bản chất cơ học, sự gia tăng của hệ số n làm giảm tỷ lệ pha gốm có độ cứng cao, từ đó làm giảm ma trận độ cứng uốn D của kết cấu và kéo theo sự gia tăng tuyến tính của trường chuyển vị. Sự tương tác giữa lực dọc và mômen uốn do tính không đối xứng của vật liệu qua mặt trung bình làm cho trục trung hòa bị dịch chuyển về phía mặt giàu gốm.

Dữ liệu tính toán trong luận văn được trực quan hóa hiệu quả thông qua các biểu đồ hội tụ của độ võng theo số lượng phần tử trên cạnh tấm và đồ thị quan hệ giữa độ võng không thứ nguyên với tỷ số h/L. Bên cạnh đó, các biểu đồ phân bố ứng suất pháp sigma-xx và ứng suất cắt tau-xz theo chiều dày không thứ nguyên z/h đã thể hiện trực quan quy luật phi tuyến phức tạp của trường ứng suất. Kết quả so sánh dạng bảng số liệu với các phương pháp kp-Ritz, phần tử tứ giác MITC4 và phần tử trơn nút NS-DSG3 cho thấy phần tử ES-MITC3+ đạt tốc độ hội tụ nhanh vượt trội ngay từ cấp lưới 16x16x2, khẳng định tính ổn định cao trong phân tích kết cấu.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể:

  • Tích hợp ma trận khối lượng trơn phục vụ phân tích động lực học: Các nhóm nghiên cứu cơ học tính toán tại các trường đại học kỹ thuật cần phát triển bổ sung ma trận khối lượng trơn trên miền cạnh cho phần tử ES-MITC3+ trong thời gian 6 đến 12 tháng tới, nhằm nâng cao độ chính xác dự báo tần số dao động riêng của tấm FGMs lên mức tin cậy trên 98%.
  • Ứng dụng mô hình hóa kết cấu trong điều kiện nhiệt độ cao: Các kỹ sư thiết kế trong ngành hàng không vũ trụ và năng lượng nên đưa thuật toán ES-MITC3+ vào quy trình tính toán tấm vỏ chịu tải nhiệt độ trên 800 độ C trong giai đoạn cuối năm 2024, giúp giảm thiểu 30% nguy cơ nứt vỡ cục bộ và phá hủy do sốc nhiệt.
  • Tối ưu hóa chỉ số gradient vật liệu n theo yêu cầu chịu lực: Các kỹ sư kết cấu công trình cần áp dụng mô hình phân tích tĩnh này kết hợp với các thuật toán tối ưu hóa đa mục tiêu trong giai đoạn 2024-2025, nhằm xác định giá trị n tối ưu từ 0.5 đến 1.5, giúp tiết kiệm 15% đến 25% khối lượng vật liệu gốm cao cấp mà vẫn đảm bảo độ võng kết cấu nằm trong giới hạn cho phép.
  • Thương mại hóa và xây dựng module phần mềm mở rộng: Đội ngũ phát triển phần mềm mô phỏng kỹ thuật CAE nên tiến hành chuyển đổi mã nguồn MATLAB hiện tại sang ngôn ngữ C++ hoặc Python trong thời hạn 18 tháng, tạo ra các module tính toán có khả năng xử lý bài toán trên 50.000 phần tử trong thời gian dưới 90 giây.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chuyên môn:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Cơ học kết cấu: Tiếp cận chi tiết phương pháp xây dựng công thức phần tử hữu hạn trơn (S-FEM) và kỹ thuật khử khóa cắt MITC3+; Ứng dụng trực tiếp vào việc thực hiện các luận văn, luận án về phân tích kết cấu composite và vật liệu thông minh.
  • Kỹ sư tính toán kết cấu trong ngành Hàng không và Cơ khí chính xác: Nắm bắt phương pháp dự báo chính xác ứng xử uốn và phân bố ứng suất tiếp của các tấm chịu tải trọng cơ học lớn; Ứng dụng trong việc thiết kế tấm vỏ máy bay, tấm chắn nhiệt tàu vũ trụ và tua-bin khí.
  • Giảng viên và các nhà nghiên cứu tại các viện kỹ thuật: Sử dụng toàn bộ khung lý thuyết biến dạng cắt bậc cao của Reddy và thuật toán làm trơn biến dạng làm tài liệu giảng dạy cho các học phần Phương pháp phần tử hữu hạn nâng cao và Cơ học vật liệu tiên tiến.
  • Lập trình viên và chuyên gia phát triển phần mềm tính toán CAE: Tham khảo thuật toán rời rạc hóa và tích phân đường theo định lý Green để nhúng vào phần lõi xử lý số của các gói phần mềm phân tích kết cấu nội địa hoặc mã nguồn mở.

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

  • Hiện tượng khóa cắt trong phân tích tấm là gì và phần tử ES-MITC3+ xử lý như thế nào? Khóa cắt là hiện tượng phần tử tấm trở nên quá cứng giả tạo khi chiều dày h/L giảm xuống dưới 0.01, khiến độ võng tính toán bị suy giảm nghiêm trọng. Phần tử ES-MITC3+ khắc phục hoàn toàn hiện tượng này bằng cách xấp xỉ lại trường biến dạng cắt ngoài mặt phẳng thông qua 6 điểm buộc hình học chuyên biệt, giúp kết quả luôn ổn định ngay cả khi tấm siêu mỏng.
  • Vì sao luận văn lựa chọn lý thuyết biến dạng cắt bậc cao HSDT thay vì FSDT? Lý thuyết biến dạng cắt bậc nhất FSDT giả định biến dạng cắt là hằng số và bắt buộc phải dùng hệ số hiệu chỉnh lực cắt vốn không chính xác cho vật liệu không đồng nhất như FGMs. Lý thuyết HSDT bậc ba của Reddy mô tả chính xác quy luật biến dạng cắt phi tuyến và tự triệt tiêu ứng suất tại hai bề mặt tấm mà không cần bất kỳ hệ số hiệu chỉnh nào.
  • Chỉ số phân bố thể tích vật liệu n ảnh hưởng như thế nào đến độ cứng của tấm FGMs? Chỉ số n quyết định tỷ lệ pha gốm và kim loại theo chiều dày. Khi n = 0, tấm là 100% gốm với mô đun đàn hồi cao 200 GPa. Khi n tăng lên 2 hoặc 10, hàm lượng kim loại nhôm (mô đun đàn hồi 70 GPa) tăng lên, làm giảm độ cứng chống uốn tổng thể và khiến độ võng tại tâm tấm tăng lên hơn 90%.
  • Ưu điểm nổi bật của phần tử tam giác 3 nút so với phần tử tứ giác 4 nút là gì? Phần tử tam giác 3 nút có khả năng thích ứng hình học vượt trội, dễ dàng chia lưới tự động cho các kết cấu có góc xiên hoặc biên dạng phức tạp. Khi được tích hợp kỹ thuật làm trơn ES-FEM và nút nổi MITC3+, phần tử tam giác đạt độ chính xác tương đương phần tử tứ giác nhưng giảm đáng kể thời gian xử lý tiền kỳ.
  • Độ chính xác của phần tử ES-MITC3+ đã được kiểm chứng qua những mô hình nào? Nghiên cứu đã đối chiếu kết quả trên 3 bài toán chuẩn gồm: tấm vuông FGMs Al/ZrO2-1 chịu tải phân bố đều với nhiều tỷ số h/L, tấm Al/Al2O3 chịu tải hình sin và tấm xiên FGMs Al/ZrO2-1 với góc nghiêng từ 15 đến 75 độ. Mọi kết quả đều trùng khít với các nghiên cứu giải tích Quasi-3D và phương pháp số kp-Ritz.

Kết luận

  • Xây dựng thành công công thức phần tử hữu hạn trơn tam giác 3 nút ES-MITC3+ dựa trên lý thuyết biến dạng cắt bậc ba của Reddy (HSDT 7 thành phần chuyển vị).
  • Khử triệt để hiện tượng khóa cắt, cho kết quả độ võng và ứng suất chính xác cao với sai số dưới 1.5% trên dải tỷ số chiều dày rộng từ 0.25 xuống đến 0.001.
  • Làm sáng tỏ quy luật ảnh hưởng của chỉ số vật liệu n từ 0 đến 10 và góc xiên từ 15 đến 75 độ đối với ứng xử uốn tĩnh của tấm FGMs.
  • Cung cấp bộ công cụ mô phỏng số ổn định, hội tụ nhanh trên lưới 24x24x2 phần tử tam giác, tối ưu hóa đáng kể năng lực tính toán kết cấu.
  • Định hướng mở rộng bài toán sang phân tích dao động tự do, ổn định động và ứng xử nhiệt - cơ phức hợp trong lộ trình 12 đến 24 tháng tới.

Công trình là bước tiến quan trọng trong việc ứng dụng phương pháp phần tử hữu hạn trơn để phân tích các kết cấu vật liệu tiên tiến. Độc giả, kỹ sư và các nhà nghiên cứu quan tâm có thể ứng dụng ngay mô hình tính toán này để giải quyết các bài toán thiết kế kết cấu tấm phân lớp chức năng trong thực tiễn.