Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, các tấm composite nhiều lớp truyền thống thường đối mặt với nguy cơ phá hủy nghiêm trọng do hiện tượng tập trung ứng suất cắt và tách lớp tại bề mặt tiếp xúc khi làm việc dưới chênh lệch nhiệt độ lớn lên đến 2000K. Để khắc phục nhược điểm này, vật liệu phân lớp chức năng (Functionally Graded Material - FGM) với cơ chế biến thiên thành phần liên tục theo chiều dày đã trở thành bước đột phá vượt bậc. Vấn đề cốt lõi đặt ra là việc mô phỏng chính xác ứng xử cơ học của tấm sandwich FGM đòi hỏi sự cân bằng giữa độ chính xác và chi phí tính toán.

Lý thuyết tấm cổ điển (CPT) bỏ qua hoàn toàn biến dạng cắt ngang nên chỉ phù hợp với tấm rất mỏng có tỷ lệ cạnh trên chiều dày lớn hơn 100. Ngược lại, các lý thuyết biến dạng cắt bậc cao (HSDT) chứa từ 9 đến 13 ẩn số khiến phương trình vi phân chuyển động trở nên vô cùng phức tạp, đòi hỏi thời gian xử lý ma trận lớn. Trong khi đó, lý thuyết biến dạng cắt bậc nhất (FSDT) tiêu chuẩn với 5 biến số lại mặc định gán hệ số hiệu chỉnh cắt bằng giá trị cố định 5/6 (tương đương 0,833), dẫn đến sai số vượt mức 12% khi phân tích vật liệu không đồng nhất FGM.

Mục tiêu trọng tâm của nghiên cứu là phát triển mô hình tấm sandwich chức năng FGM cải tiến dựa trên lý thuyết biến dạng cắt bậc nhất rút gọn chỉ còn 4 ẩn số, đồng thời xây dựng công thức xác định hệ số hiệu chỉnh cắt thích ứng theo phân bố cơ tính thực tế. Nghiên cứu tập trung khảo sát dao động tự do và ổn định đàn hồi của các cấu hình tấm sandwich FGM kim loại - gốm với tỷ lệ chiều dày 1-2-1, 2-2-1 và chỉ số mũ phân bố từ 0 đến 10. Kết quả thu được giúp giảm thiểu 20% số lượng biến chuyển vị so với FSDT truyền thống, duy trì độ chính xác cao với sai số dưới 1,5% so với nghiệm giải tích ba chiều chính xác.

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 ứng dụng nền tảng của Lý thuyết biến dạng cắt bậc nhất cải tiến (Refined First-Order Shear Deformation Theory - RFSDT) kết hợp với Mô hình đồng hóa vi mô Voigt. Trường chuyển vị của tấm được phân tách độc lập thành hai thành phần uốn và cắt, biểu diễn qua 4 biến chuyển vị cơ bản thay vì 5 hoặc nhiều hơn như các mô hình cũ. Phương trình chuyển động và các điều kiện biên được thiết lập chặt chẽ thông qua việc áp dụng Nguyên lý biến phân Hamilton cho hệ cơ học bảo toàn năng lượng.

Khung lý thuyết của đề tài xoay quanh các khái niệm then chốt:

  • Tấm sandwich chức năng FGM: Kết cấu gồm 3 lớp với lõi đồng nhất liên kết cùng hai lớp bề mặt FGM có hàm mật độ thể tích biến thiên liên tục theo quy luật hàm mũ.
  • Hệ số hiệu chỉnh cắt thích ứng: Đại lượng hiệu chỉnh trường ứng suất cắt ngang thỏa mãn điều kiện triệt tiêu ứng suất tại hai bề mặt tự do trên và dưới của tấm.
  • Tải trọng tới hạn mất ổn định đàn hồi: Giá trị lực nén màng lớn nhất mà kết cấu có thể chịu đựng trước khi chuyển sang trạng thái biến dạng mất ổn định.
  • Tần số dao động tự do không thứ nguyên: Thông số động lực học phản ánh đặc tính dao động riêng của kết cấu theo các bậc dao động khác nhau.

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

Nghiên cứu sử dụng nguồn dữ liệu số chuẩn hóa từ 6 cấu hình hình học tấm sandwich FGM điển hình bao gồm: 1-0-1, 1-1-1, 1-2-1, 2-1-1, 2-1-2 và 2-2-1. Bộ dữ liệu vật liệu khảo sát gồm 3 cặp vật liệu gốm - kim loại tiêu chuẩn quốc tế: Nhôm - Alumina ($Al/Al_2O_3$ với mô đun đàn hồi $E_m = 70$ GPa, $E_c = 380$ GPa), Nhôm - Silicon Carbide ($Al/SiC$ với $E_c = 427$ GPa) và Nhôm - Tungsten Carbide ($Al/WC$ với $E_c = 717$ GPa), cùng hệ số Poisson cố định ở mức 0,3.

Phương pháp giải tích Navier dựa trên chuỗi Fourier kép được lựa chọn để giải phương trình vi phân chuyển động cho tấm chữ nhật có 4 biên tựa đơn (SSSS). Đối với các điều kiện biên phức tạp khác như 4 biên ngàm (CCCC), ngàm - tựa (CCSS) và ngàm - tự do (FFCC), phương pháp xấp xỉ dạng sóng lượng giác Sobhy được áp dụng nhằm đảm bảo tính hội tụ nhanh và ổn định nghiệm. Toàn bộ thuật toán được lập trình tự động hóa trên nền tảng phần mềm Matlab với hơn 1000 lượt mô phỏng tham số trong thời gian triển khai từ tháng 6/2013 đến tháng 11/2013.

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

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

Thứ nhất, hệ số hiệu chỉnh cắt của tấm sandwich FGM không phải là một hằng số 5/6 như các nghiên cứu kinh điển giả định. Với tấm vuông $Al/Al_2O_3$ cấu hình 1-2-1, hệ số hiệu chỉnh cắt biến thiên phi tuyến rõ rệt trong khoảng từ 0,582 đến 0,845 tùy thuộc vào chỉ số mũ phân bố thể tích $p$ và tỷ số tương phản độ cứng giữa mặt gốm và kim loại. Khi tỷ số mô đun $E_t/E_b$ tăng cao, hệ số này giảm mạnh, phản ánh sự thay đổi lớn trong phân bố ứng suất cắt ngang qua chiều dày tấm.

Thứ hai, chỉ số mũ hàm mật độ thể tích $p$ đóng vai trò chi phối toàn diện đến độ cứng tương đương của tấm sandwich. Khi chỉ số $p$ tăng từ 0 (hoàn toàn là gốm) lên 10 (thành phần kim loại chiếm ưu thế), tần số dao động tự nhiên của 3 mode đầu tiên giảm trung bình từ 35,4% đến 48,2%. Đồng thời, tải trọng mất ổn định tới hạn không thứ nguyên chịu lực nén đơn trục suy giảm hơn 52,6% do độ cứng uốn tổng thể của hệ kết cấu bị suy giảm tương ứng.

Thứ ba, ảnh hưởng của hiệu ứng biến dạng cắt ngang phụ thuộc chặt chẽ vào độ mảnh tấm ($b/h$). Đối với tấm dày ($b/h = 5$), việc sử dụng hệ số hiệu chỉnh cắt cố định 5/6 gây ra sai số từ 8,5% đến 12,3% đối với tải trọng mất ổn định so với mô hình hiệu chỉnh thích ứng. Tuy nhiên, khi tấm trở nên mỏng hơn với $b/h \ge 50$, sai số này nhanh chóng giảm xuống dưới 1,2%, chứng minh tính nhạy cảm đặc biệt của hiệu ứng cắt trong kết cấu tấm sandwich dày.

Thứ tư, điều kiện liên kết chu vi làm thay đổi căn bản khả năng chịu lực và tần số dao động riêng. Kết quả mô phỏng chỉ ra rằng tấm 4 biên ngàm (CCCC) có tải trọng mất ổn định tới hạn cao gấp 2,8 đến 3,2 lần và tần số dao động tự nhiên cao hơn khoảng 1,8 lần so với tấm 4 biên tựa đơn (SSSS) có cùng kích thước hình học và tỷ lệ vật liệu.

Thảo luận kết quả

Nguyên nhân căn bản của các hiện tượng trên bắt nguồn từ sự phân bố vật chất không đồng nhất qua chiều dày tấm FGM, làm lệch trục trung hòa hình học ra khỏi trục trung hòa cơ học. Điều này tạo ra sự tương tác phức tạp giữa biến dạng màng và biến dạng uốn được thể hiện qua ma trận độ cứng tương tác $B_{ij}$. Mô hình 4 biến cải tiến đã mô tả chính xác trường ứng suất tiếp thông qua việc tích phân trực tiếp các phương trình cân bằng ứng suất phẳng, giúp phản ánh chân thực năng lượng biến dạng cắt của hệ.

So sánh đối chiếu với các công bố quốc tế sử dụng lý thuyết biến dạng cắt bậc cao (HSDT của Reddy và Zenkour) cho thấy độ lệch kết quả của mô hình đề xuất chỉ dao động trong phạm vi từ 0,6% đến 1,4%, nhưng cấu trúc ma trận giải tích được tinh giản đáng kể. Dữ liệu nghiên cứu được biểu diễn trực quan qua biểu đồ tương tác giữa tần số dao động và tải trọng nén màng, thể hiện quy luật suy giảm tuyến tính của bình phương tần số dao động khi lực nén màng tiệm cận về ngưỡng mất ổn định tới hạn.

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

Thứ nhất, tích hợp mô hình biến dạng cắt bậc nhất cải tiến 4 biến vào các phần mềm phần tử hữu hạn thương mại như ANSYS hoặc ABAQUS thông qua việc xây dựng các phần tử người dùng (UEL), nhằm giảm thiểu ít nhất 30% thời gian giải ma trận độ cứng cho các bài toán phân tích tấm sandwich FGM quy mô lớn.

Thứ hai, áp dụng quy trình tính toán hệ số hiệu chỉnh cắt thích ứng theo phân bố hàm mũ $p$ thay thế hoàn toàn hệ số mặc định 5/6 trong các tiêu chuẩn thiết kế kết cấu composite nhiều lớp, đặc biệt đối với các tấm có tỷ lệ cạnh trên chiều dày nhỏ hơn 20 làm việc trong điều kiện tải trọng động.

Thứ ba, tối ưu hóa chỉ số mũ phân bố vật liệu $p$ trong khoảng thiết kế từ 0,5 đến 2,0 và ưu tiên lựa chọn tỷ lệ chiều dày lớp 1-2-1 đối với các vách ngăn chịu nhiệt trong công nghiệp hàng không vũ trụ, nhằm đạt tỷ số tối ưu giữa độ cứng uốn và khối lượng kết cấu với mục tiêu tăng 25% khả năng chịu tải tới hạn.

Thứ tư, triển khai các chương trình thử nghiệm thực nghiệm kiểm chứng trên mẫu tấm sandwich FGM chế tạo bằng công nghệ luyện kim bột hoặc in 3D laser kim loại - gốm trong lộ trình 12 đến 24 tháng tới, do các viện nghiên cứu vật liệu chuyên ngành chủ trì thực hiện nhằm hoàn thiện cơ sở dữ liệu thực nghiệm.

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

  • Kỹ sư thiết kế kết cấu hàng không, hàng hải và xây dựng công trình đặc biệt: Sử dụng trực tiếp các bảng tra cứu tải trọng tới hạn và tần số dao động để tính toán các tấm vỏ tàu thủy, vách ngăn chống cháy tàu vũ trụ và sàn kết cấu chịu tải cơ nhiệt phức tạp.
  • Nhà nghiên cứu và học viên cao học ngành Cơ học tính toán: Tiếp cận phương pháp thiết lập phương trình chuyển động từ nguyên lý Hamilton và kỹ thuật giải tích Navier/Sobhy cho các mô hình vật liệu phân lớp chức năng nâng cao.
  • Chuyên gia phát triển phần mềm mô phỏng kỹ thuật (CAE): Ứng dụng thuật toán 4 biến rút gọn để tối ưu hóa module phần tử tấm vỏ trong các bộ giải kết cấu, giúp tiết kiệm bộ nhớ RAM và tăng tốc độ xử lý bài toán động lực học.
  • Doanh nghiệp sản xuất vật liệu composite và vật liệu tiên tiến: Khai thác các phân tích về tỷ lệ phân bố thể tích và độ tương phản vật liệu để định hình quy trình chế tạo tấm sandwich FGM kim loại - gốm đạt độ bền uốn và độ bền nhiệt tối ưu.

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

Vật liệu phân lớp chức năng FGM có ưu điểm gì vượt trội so với composite nhiều lớp truyền thống?
Vật liệu FGM loại bỏ hoàn toàn sự gián đoạn cơ tính tại mặt tiếp xúc nhờ cơ chế biến thiên thành phần liên tục theo chiều dày. Điều này ngăn ngừa triệt để hiện tượng tập trung ứng suất tiếp và chống tách lớp khi làm việc ở nhiệt độ lên tới 2000K, mang lại độ bền cơ - nhiệt vượt trội cho kết cấu.

Tại sao mô hình 4 biến cải tiến lại hiệu quả hơn lý thuyết biến dạng cắt bậc cao HSDT?
Mô hình 4 biến phân tách chuyển vị ngang thành thành phần uốn và cắt độc lập, giúp loại bỏ một biến chuyển vị xoay so với FSDT và giảm từ 5 đến 9 biến so với HSDT. Mô hình mang lại độ chính xác tương đương HSDT với sai số dưới 1,5% nhưng giải thuật đơn giản và khối lượng tính toán giảm hơn 30%.

Hệ số hiệu chỉnh cắt 5/6 truyền thống gây ra sai số như thế nào đối với tấm FGM?
Hệ số 5/6 chỉ đúng cho tấm đồng nhất đẳng hướng. Đối với tấm sandwich FGM có mô đun đàn hồi thay đổi liên tục, hệ số hiệu chỉnh thực tế dao động từ 0,582 đến 0,845. Việc áp dụng cố định giá trị 5/6 sẽ gây sai số vượt quá 12% khi tính toán tải trọng mất ổn định của tấm dày.

Chỉ số mũ hàm mật độ thể tích $p$ ảnh hưởng ra sao đến tần số dao động và tải trọng mất ổn định?
Khi chỉ số $p$ tăng từ 0 đến 10, tỷ lệ thể tích kim loại trong tấm tăng dần làm giảm độ cứng uốn tổng thể. Kết quả tính toán cho thấy tần số dao động tự nhiên của tấm giảm từ 35,4% đến 48,2% và tải trọng mất ổn định tới hạn suy giảm hơn 50% so với trạng thái gốm nguyên chất.

Phương pháp giải tích nào được sử dụng để xử lý các điều kiện biên khác nhau của tấm?
Nghiên cứu sử dụng lời giải Navier dạng chuỗi Fourier kép cho tấm có 4 cạnh tựa đơn (SSSS) và áp dụng phương pháp xấp xỉ sóng lượng giác Sobhy cho các điều kiện biên ngàm (CCCC), ngàm - tựa (CCSS), tự do - ngàm (FFCC), đảm bảo nghiệm số hội tụ chính xác tuyệt đối.

Kết luận

  • Xây dựng thành công mô hình tấm sandwich chức năng FGM 4 biến cải tiến, tinh gọn số lượng phương trình chuyển động nhưng vẫn bảo toàn đầy đủ ứng xử cơ học của hệ.
  • Thiết lập công thức giải tích tính toán hệ số hiệu chỉnh cắt thích ứng, khắc phục triệt để hạn chế sai số lớn của hệ số 5/6 truyền thống trên vật liệu FGM.
  • Cung cấp bộ dữ liệu giải tích toàn diện về tần số dao động tự nhiên và tải trọng mất ổn định tới hạn cho 6 cấu hình tấm và 4 dạng điều kiện biên khác nhau.
  • Chứng minh độ chính xác vượt trội của mô hình với sai số dưới 1,5% so với lý thuyết HSDT và lời giải giải tích ba chiều chính xác.
  • Mở ra hướng nghiên cứu mở rộng sang phân tích phi tuyến hình học, ứng xử nhiệt - cơ kết hợp và nứt gãy động lực học trong kế hoạch 12 đến 18 tháng tới. Các kỹ sư và nhà nghiên cứu kết cấu được khuyến nghị ứng dụng ngay khung tính toán này vào quy trình thiết kế và tối ưu hóa kết cấu composite tiên tiến.