Tổng quan nghiên cứu

Vật liệu composite có cơ tính biến đổi (Functionally Graded Materials - FGM) là bước tiến đột phá của khoa học vật liệu hiện đại kể từ khi được khởi xướng tại Viện Sendai, Nhật Bản vào năm 1984. Khác với composite phân lớp truyền thống dễ bị tách lớp hoặc tập trung ứng suất tại mặt phân giới, vật liệu FGM hai thành phần gồm kim loại và gốm sở hữu cấu trúc biến đổi liên tục theo chiều dày. Cấu trúc này cho phép kết cấu vừa duy trì độ dẻo dai cơ học của kim loại, vừa tận dụng khả năng kháng nhiệt vượt trội của gốm ở môi trường làm việc khắc nghiệt từ 1500 độ C đến trên 2500 độ C.

Trong các ứng dụng thực tế của ngành hàng không vũ trụ, năng lượng hạt nhân và công trình biển, tấm composite FGM thường được đặt trên các lớp nền đàn hồi và chịu tác động đồng thời của tải nén dọc trục và tải nhiệt tăng đều. Tuy nhiên, các công bố khoa học trong nước về bài toán ổn định tĩnh phi tuyến và ứng xử sau vồng của tấm FGM đặt trên nền đàn hồi hai thông số vẫn còn rất hạn chế trước giai đoạn năm 2011.

Luận văn thạc sĩ chuyên ngành Cơ học Kỹ thuật của tác giả Đỗ Nam dưới sự hướng dẫn của Phó Giáo sư, Tiến sĩ Nguyễn Đình Đức đã tập trung giải quyết bài toán cốt lõi này. Mục tiêu nghiên cứu là thiết lập hệ phương trình giải tích, xác định chính xác tải trọng tới hạn và khảo sát định lượng toàn diện mối quan hệ phi tuyến giữa tải trọng và độ võng trong giai đoạn sau vồng. Nghiên cứu thực hiện trên mô hình tấm chữ nhật FGM kim loại - gốm (Al/Al2O3 và Ti-6Al-4V/ZrO2) tựa bản lề bốn cạnh trên nền đàn hồi Pasternak và Winkler. Kết quả của luận văn đóng góp hệ thống công thức giải tích chuẩn xác, giúp giảm thiểu rủi ro phá hủy kết cấu và tối ưu hóa chi phí thiết kế trong kỹ thuật công trình hiện đại.

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 xây dựng trên nền tảng lý thuyết tấm mỏng cổ điển kết hợp với quan hệ hình học phi tuyến bậc hai von Kármán nhằm mô tả chính xác các biến dạng lớn trong giai đoạn sau khi kết cấu bị vồng. Ba khái niệm trung tâm tạo nên khung lý thuyết bao gồm: hàm phân bố tỷ lệ thể tích theo luật lũy thừa, mô hình nền đàn hồi hai thông số Pasternak - Winkler và tiêu chuẩn ổn định tĩnh.

Cơ tính hiệu dụng của vật liệu FGM như mô đun đàn hồi E và hệ số giãn nở nhiệt được xác định thông qua quy tắc hỗn hợp theo hàm lũy thừa của tọa độ chiều dày z với chỉ số thể tích N dao động từ 0 đến 1. Khi N bằng 0, tấm trở thành gốm thuần túy với mô đun đàn hồi đạt xấp xỉ 380 GPa; khi N bằng 1, cơ tính phân bố tuyến tính giữa nhôm (70 GPa) và gốm Alumina. Mô hình nền tương tác với kết cấu thông qua tải trọng phản lực kết hợp độ cứng lò xo tuyến tính Winkler K1, độ cứng lớp cắt trượt Pasternak K2 và thành phần phi tuyến bậc ba K3. Tiêu chuẩn cân bằng tĩnh của Euler được áp dụng để xác định trạng thái rẽ nhánh từ cấu hình ban đầu sang dạng mất ổn định.

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

Nguồn dữ liệu tính toán được thiết lập từ các thông số vật lý thực nghiệm chuẩn quốc tế của hệ vật liệu Al/Al2O3 và Ti-6Al-4V/ZrO2 với chiều dày tấm h và tỷ lệ kích thước hình học tiêu chuẩn a/b bằng 1. Quá trình phân tích áp dụng phương pháp giải tích nửa giải tích thông qua việc đưa vào hàm ứng suất Airy để thỏa mãn đồng nhất các phương trình cân bằng nội lực trong mặt phẳng.

Nghiên cứu sử dụng chuỗi lượng giác Navier hai chiều để biểu diễn độ võng và độ bất toàn hình học ban đầu với tham số độ lệch mu dao động từ 0 đến 0.1. Sau đó, thuật toán tích phân xấp xỉ Galerkin được triển khai để chuyển hệ phương trình đạo hàm riêng phi tuyến thành phương trình đại số tường minh xác định tải nén và tải nhiệt tới hạn. Phương pháp này được lựa chọn vì tính chính xác cao, thời gian xử lý nhanh gấp nhiều lần phương pháp phần tử hữu hạn, đồng thời khảo sát chi tiết hơn 90 kịch bản biến thiên của các hệ số nền K1 từ 0 đến 300 và K2 từ 0 đến 250 trong suốt tiến trình hoàn thành luận văn năm 2011.

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

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

Thứ nhất, sự hiện diện của nền đàn hồi Pasternak đóng vai trò quyết định trong việc gia tăng khả năng mang tải của kết cấu tấm FGM. Khi khảo sát tấm hoàn hảo với chỉ số thể tích N bằng 1, việc tăng hệ số Winkler K1 từ 0 lên 300 làm tải nén tới hạn tăng từ 0.4280 GPa lên 0.7575 GPa, tương đương mức tăng trưởng 77.0%. Khi kết hợp thêm hệ số cắt Pasternak K2 bằng 200, tải nén tới hạn đạt mức 1.0075 GPa, tăng hơn 135.4% so với trường hợp không có nền đàn hồi.

Thứ hai, chỉ số tỷ lệ thể tích N có ảnh hưởng nghịch biến rõ rệt đối với độ cứng uốn của tấm FGM. Khi chuyển đổi cấu trúc vật liệu từ dạng phân bố tuyến tính N bằng 1 sang tấm gốm thuần túy N bằng 0, tải nén tới hạn gốc tăng vọt từ 0.4280 GPa lên 0.8586 GPa, tức tăng hơn 100.6%. Tương tự, ở trạng thái nền cứng K1 bằng 300 và K2 bằng 200, giá trị tải vồng của tấm gốm thuần túy đạt đỉnh 1.7697 GPa so với mức 1.0075 GPa của tấm có N bằng 1.

Thứ ba, nghiên cứu làm sáng tỏ ứng xử phi tuyến sau vồng của tấm có độ bất toàn hình học ban đầu mu bằng 0.1. Khi tỷ số độ võng cực đại trên chiều dày W/h phát triển từ 0.5 lên 2.0 trong điều kiện K1 bằng 300 và K2 bằng 100, năng lực chịu tải nén của tấm tăng từ 0.7961 GPa lên 1.6038 GPa, tức tăng hơn 101.4%. Dưới tác dụng của tải nhiệt tăng đều, biến thiên nhiệt độ tới hạn cũng ghi nhận mức tăng tỷ lệ thuận với độ cứng của nền đàn hồi, khẳng định vai trò chống mất ổn định nhiệt của lớp nền đệm.

Thảo luận kết quả

Cơ chế vật lý giải thích cho sự gia tăng vượt bậc của tải vồng tới hạn nằm ở hiệu ứng chống uốn kép do lớp cắt trượt Pasternak K2 mang lại, giúp phân tán đều ứng suất cục bộ thay vì chỉ chịu nén độc lập như mô hình lò xo đơn Winkler K1. Kết quả này hoàn toàn tương thích với các nghiên cứu lý thuyết tấm biến dạng trượt bậc cao của các nhóm tác giả quốc tế nhưng có ưu thế nổi bật về mặt công thức tường minh.

Dữ liệu tính toán trong luận văn được hệ thống hóa trực quan thông qua 95 bảng số liệu chi tiết và hơn 20 đồ thị biểu diễn đường cong quan hệ tải trọng - độ võng (P-W). Đồ thị cho thấy rõ hiện tượng rẽ nhánh đối với tấm phẳng lý tưởng và đường cong phát triển liên tục không có điểm gãy đối với tấm có độ bất toàn ban đầu. Sự so sánh định lượng giữa mô hình phân tích giải tích này với các kết quả phần tử hữu hạn cho thấy độ lệch chỉ dao động trong khoảng 2.0% đến 3.5%, chứng minh độ tin cậy tuyệt đối của phương pháp.

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

Thứ nhất, các đơn vị thiết kế kết cấu hàng không và năng lượng hạt nhân cần tối ưu hóa chỉ số thành phần thể tích N trong khoảng từ 0.2 đến 0.5. Giải pháp này giúp tăng tỷ lệ gốm ở bề mặt chịu nhiệt, qua đó nâng cao giới hạn chịu nhiệt lên trên 1200 độ C và giảm thiểu tới 25% ứng suất nhiệt cục bộ tại các điểm liên kết khung vỏ trong giai đoạn vận hành từ 2 đến 5 năm tới.

Thứ hai, các kỹ sư cơ khí và xây dựng công trình ngầm cần chủ động tích hợp lớp đệm đàn hồi tương đương mô hình Pasternak với hệ số cắt K2 tối thiểu từ 50 đến 100 vào bệ đỡ tấm vỏ. Việc gia cố nền đệm này sẽ giúp cải thiện ngay lập tức từ 30% đến 50% khả năng kháng uốn và ổn định sau vồng của kết cấu dưới tác động của tải nén cơ học bất lợi.

Thứ ba, các cơ sở gia công chế tạo vật liệu FGM phải kiểm soát nghiêm ngặt độ bất toàn hình học ban đầu nhằm duy trì hệ số mu dưới ngưỡng 0.05. Việc kiểm soát dung sai độ phẳng bề mặt bằng công nghệ lắng đọng hơi hóa học hoặc thiêu kết plasma sẽ hạn chế tối đa nguy cơ suy giảm tới hạn tải trọng trong 1000 giờ làm việc đầu tiên.

Thứ tư, các viện nghiên cứu chuyên ngành Cơ học tính toán cần xây dựng các gói công cụ phần mềm tự động hóa tính toán tấm vỏ FGM dựa trên phương pháp giải tích Galerkin đã thiết lập. Kế hoạch này nên hoàn thành trong vòng 12 tháng nhằm phục vụ công tác thẩm định an toàn chịu lực cho các dự án công nghệ cao cấp quốc gia.

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

Nhóm kỹ sư thiết kế kết cấu hàng không vũ trụ và quốc phòng sẽ tìm thấy trong luận văn phương pháp tính toán chuẩn xác khả năng chịu tải nén và tải nhiệt của tấm vỏ FGM trên thân máy bay siêu âm hoặc cánh lái tên lửa, giúp tối ưu tỷ lệ trọng lượng và độ bền.

Nhóm nghiên cứu và phát triển vật liệu composite tiên tiến có thể khai thác các công thức xác định cơ tính hiệu dụng và quy tắc hỗn hợp lũy thừa để định lượng tỷ lệ pha kim loại - gốm, ứng dụng trong việc chế tạo các tấm chắn nhiệt cho lò phản ứng hạt nhân.

Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Cơ học Kỹ thuật, Cơ học Vật rắn biến dạng có thể sử dụng toàn bộ hệ thống phương trình vi phân và phương pháp tích phân Galerkin trong tài liệu như một giáo trình tham khảo mẫu mực cho các đề tài nghiên cứu chuyên sâu.

Các chuyên gia tư vấn giám định công trình công nghiệp nặng có thể ứng dụng các bảng tra cứu hệ số nền đàn hồi Pasternak và Winkler để đánh giá mức độ an toàn của kết cấu sàn đáy và bể chứa chịu nhiệt đặt trực tiếp trên nền đất hoặc nền móng gia cố.

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

Vật liệu FGM khác biệt như thế nào so với composite phân lớp thông thường? Vật liệu FGM sở hữu thành phần biến đổi liên tục từ kim loại sang gốm theo chiều dày, giúp triệt tiêu hoàn toàn mặt phân giới rời rạc. Điều này loại bỏ triệt để hiện tượng bong tách lớp và giảm tập trung ứng suất, nâng cao độ bền nhiệt từ 2000 độ C trở lên trong các điều kiện làm việc khắc nghiệt.

Tại sao mô hình nền đàn hồi Pasternak lại vượt trội hơn mô hình Winkler? Mô hình Winkler chỉ xem nền móng như các lò xo độc lập không có tương tác ngang. Trong khi đó, mô hình Pasternak bổ sung thêm lớp cắt trượt K2, phản ánh chính xác sự liên kết giữa các phần tử đất nền lân cận, giúp kết quả tính toán tải vồng tăng từ 15% đến 40% sát với ứng xử thực tế.

Chỉ số thể tích N ảnh hưởng như thế nào đến độ ổn định của tấm composite FGM? Chỉ số N điều chỉnh tỷ lệ giữa gốm và kim loại trong kết cấu. Khi N tiến về 0, hàm lượng gốm tăng cao làm tăng mô đun đàn hồi tổng thể lên 380 GPa, giúp tải nén tới hạn tăng hơn 100% so với khi N bằng 1, giúp tấm chịu uốn và nén tốt hơn đáng kể.

Độ bất toàn hình học ban đầu mu tác động ra sao đến giai đoạn sau vồng? Độ bất toàn ban đầu làm mất đi điểm rẽ nhánh lý thuyết của hiện tượng vồng, khiến tấm bắt đầu võng ngay khi xuất hiện tải trọng nhỏ. Tuy nhiên, trong giai đoạn sau vồng với tỷ số W/h đạt mức 2.0, kết cấu vẫn giữ được khả năng mang tải tăng thêm nhờ hiệu ứng kéo căng màng.

Phương pháp Galerkin kết hợp nghiệm chuỗi Navier có đảm bảo độ tin cậy không? Phương pháp này cho kết quả giải tích tường minh với độ chính xác rất cao, sai số chỉ dưới 3.0% so với phương pháp phần tử hữu hạn 3D đắt đỏ. Đồng thời, cách tiếp cận này cho phép khảo sát nhanh hàng trăm tổ hợp thông số biên và độ cứng nền trong thời gian ngắn.

Kết luận

  • Luận văn đã thiết lập thành công hệ phương trình phi tuyến giải tích hoàn chỉnh nghiên cứu độ ổn định và ứng xử sau vồng của tấm composite FGM trên nền đàn hồi hai thông số.
  • Chứng minh định lượng rằng độ cứng nền Winkler K1 và lớp cắt Pasternak K2 giúp tăng khả năng chịu tải nén tới hạn của tấm từ 77.0% đến trên 135.4%.
  • Xác định rõ vai trò của chỉ số phân bố thể tích N và độ bất toàn ban đầu mu trong việc kiểm soát độ võng và trạng thái ứng suất sau vồng dưới tải cơ - nhiệt.
  • Đóng góp hệ thống 95 bảng số liệu chuẩn mực cùng phương pháp giải tích tối ưu, tạo tiền đề vững chắc cho việc tính toán thiết kế kết cấu FGM tại Việt Nam.
  • Trong thời gian tới, các đơn vị nghiên cứu cần tiếp tục mở rộng mô hình cho kết cấu vỏ cong FGM đa lớp và tích hợp điều kiện biên phức tạp nhằm ứng dụng trực tiếp vào công nghiệp chế tạo công nghệ cao.