Tổng quan nghiên cứu

Vật liệu composite hiện đại đang chiếm hơn 60% tỉ trọng ứng dụng trong các ngành công nghiệp đòi hỏi yêu cầu kỹ thuật khắt khe như hàng không vũ trụ, chế tạo tàu thủy, công nghiệp dầu khí và sản xuất vi mạch sinh học. Trong điều kiện làm việc thực tế, các kết cấu dạng tấm mỏng thường xuyên phải chịu tải trọng nhiệt lớn và gradient nhiệt độ biến đổi qua chiều dày từ 20 mm đến 50 mm, dẫn đến hiện tượng cong vênh, dạn nứt hoặc suy giảm độ bền kết cấu nghiêm trọng. Luận văn thạc sĩ khoa học chuyên ngành Cơ học vật thể rắn của tác giả Nghiêm Thị Thu Hà, dưới sự hướng dẫn khoa học của Giáo sư Nguyễn Đình Đức tại Đại học Quốc gia Hà Nội vào năm 2011, đã giải quyết bài toán cơ học nhiệt cấp thiết này.

Mục tiêu cốt lõi của nghiên cứu là thiết lập mô hình toán học giải tích nhằm tính toán chính xác độ võng của tấm composite mỏng hình chữ nhật cốt hạt Titan nền polymer PVC tựa bản lề bốn cạnh, đặt trong trường nhiệt độ dừng và không dừng. Nghiên cứu thực hiện mô phỏng số cho tấm kích thước tiêu chuẩn chiều dài 2.0 m, chiều rộng 1.0 m và chiều dày 0.02 m, chịu tác động trao đổi nhiệt đối lưu giữa hai mặt với nhiệt độ môi trường chênh lệch từ 300 K đến 330 K. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc tìm ra nghiệm giải tích tường minh dạng chuỗi lượng giác kép, giúp nâng cao hơn 25% độ chính xác trong dự báo ứng xử uốn nhiệt và xác định tỉ lệ hạt độn tối ưu, hỗ trợ đắc lực cho các kỹ sư tối ưu hóa kết cấu vật liệu composite chịu nhiệt.

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 kết hợp chặt chẽ giữa cơ học môi trường liên tục, lý thuyết đàn hồi nhiệt tuyến tính và lý thuyết tấm mỏng cổ điển Kirchhoff-Love. Hai khung lý thuyết vật lý nền tảng được áp dụng bao gồm:

  • Định luật dẫn nhiệt Fourier và nguyên lý nhiệt động lực học môi trường biến dạng đàn hồi thông qua hàm năng lượng tự do Helmholtz, phục vụ việc thiết lập phương trình vi phân truyền nhiệt tổng quát.
  • Lý thuyết tấm mỏng Kirchhoff dựa trên giả thiết pháp tuyến thẳng và bỏ qua ứng suất pháp theo phương vuông góc với mặt phẳng trung hòa, giúp quy đổi trạng thái ứng suất - biến dạng không gian 3 chiều về hệ phương trình vi phân 2 chiều cho mặt giữa.

Bên cạnh đó, nghiên cứu vận dụng các mô hình vi cơ học tiên tiến để xác định các đặc trưng cơ lý vĩ mô tương đương của vật liệu composite cốt hạt đẳng hướng. Mô đun đàn hồi thể tích hiệu dụng và mô đun trượt tương đương được tính toán theo mô hình tương tác hạt - nền đàn hồi cầu tự phối hợp. Hệ số giãn nở nhiệt thể tích của composite hạt được nội suy thông qua phương trình tương quan cơ nhiệt chính xác. Các khái niệm trung tâm xuyên suốt công trình bao gồm: trường nhiệt độ dừng, trường nhiệt độ không dừng biến thiên theo thời gian, mômen uốn nhiệt tương đương, lực kéo nén mặt phẳng và độ cứng trụ uốn của tấm mỏng.

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

Nghiên cứu sử dụng nguồn dữ liệu thông số vật liệu thực nghiệm tiêu chuẩn của composite hạt Titan phân tán trong nền Polyvinyl Clorua (PVC). Các thông số đầu vào xác định gồm: pha nền PVC có mô đun đàn hồi 3.0 GPa, hệ số Poisson 0.16; pha cốt hạt Titan có mô đun đàn hồi 100.0 GPa, hệ số Poisson 0.10. Hệ số truyền nhiệt đối lưu bề mặt trên và mặt dưới lần lượt là 60 W/(m²·K) và 40 W/(m²·K).

Cỡ mẫu khảo sát mô phỏng số bao gồm 4 mức tỉ lệ thể tích hạt độn Titan đại diện: 10%, 15%, 20% và 30%, kết hợp với 3 mốc thời gian chuyển tiếp không dừng là 120 giây, 600 giây và 1200 giây trước khi tiến tới trạng thái dừng hoàn toàn. Phương pháp chọn mẫu thông số dựa trên nguyên tắc quét biên kỹ thuật nhằm khảo sát đầy đủ từ vùng pha loãng đến vùng nồng độ hạt bão hòa.

Quy trình phân tích sử dụng phương pháp giải tích giải tích kết hợp tính toán số trên phần mềm Matlab:

  • Sử dụng phép biến đổi tích phân Laplace để giải phương trình đạo hàm riêng truyền nhiệt phi dừng qua chiều dày tấm, sau đó giải phương trình siêu việt bằng thuật toán chia đôi để thu được trường nhiệt độ chính xác.
  • Áp dụng chuỗi Fourier kép lượng giác Navier để giải phương trình vi phân bậc bốn uốn tấm. Phương pháp giải tích được lựa chọn thay vì phương pháp phần tử hữu hạn thông thường nhằm cung cấp công thức nghiệm hiển minh, loại bỏ hoàn toàn sai số làm tròn của lưới số và cho phép phân tích độ nhạy của mọi tham số hình học - vật liệu một cách tức thì.

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

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

Phân tích số liệu và biểu thức nghiệm giải tích đã mang lại bốn phát hiện khoa học mang tính đột phá về ứng xử cơ nhiệt của tấm composite:

Thứ nhất, quy luật phân bố nhiệt độ qua chiều dày tấm thể hiện tính phi tuyến rõ rệt trong giai đoạn đầu (thời gian dưới 120 giây). Khi thời gian tiến dần đến 1200 giây, trường nhiệt độ chuyển dịch hoàn toàn sang dạng tuyến tính ổn định. Tỉ lệ thể tích hạt Titan tăng từ 10% lên 30% làm tăng hệ số dẫn nhiệt tương đương, giúp giảm độ chênh lệch nhiệt độ giữa hai bề mặt tấm khoảng 18%.

Thứ hai, nghiên cứu chứng minh một quy luật quan trọng: gradient nhiệt độ truyền thuần túy theo chiều dày không gây ra bất kỳ chuyển vị nào theo phương trục tọa độ trong mặt phẳng giữa (chuyển vị u = 0, v = 0). Điều này khẳng định mặt giữa của tấm hoàn toàn không phát sinh biến dạng nhiệt phẳng, giúp giản lược cấu trúc phương trình uốn tấm.

Thứ ba, tỉ lệ thể tích hạt Titan trong khoảng từ 10% đến 20% giúp cải thiện vượt trội khả năng chống uốn của kết cấu. Khi tăng tỉ lệ hạt từ 10% lên 20%, độ võng cực đại tại tâm tấm giảm hơn 35%, đồng thời phân bố độ võng trên toàn bộ bề mặt trở nên phẳng và đồng đều hơn.

Thứ tư, khi tiếp tục tăng hàm lượng hạt Titan lên mức 30%, ứng xử uốn nhiệt của tấm xuất hiện sự biến đổi đột ngột. Độ võng tại các điểm tăng vọt khoảng 40% đến 45% so với mức nồng độ 20%, kèm theo hiện tượng tập trung độ võng không đồng đều giữa tâm và biên tấm.

Thảo luận kết quả

Nguyên nhân vật lý sâu xa của hiện tượng trên nằm ở sự tương tác cơ nhiệt giữa hai pha vật liệu. Pha hạt Titan có độ cứng cao hơn pha nền PVC xấp xỉ 33 lần (100 GPa so với 3 GPa). Khi phân tán hạt ở nồng độ hợp lý dưới 20%, các hạt cầu đóng vai trò như các chốt cơ học gia cường vi mô, làm tăng mô đun đàn hồi khối tương đương của composite, nâng cao độ cứng trụ uốn và cản trở biến dạng dãn nở nhiệt. Tuy nhiên, khi hàm lượng hạt đạt ngưỡng 30%, tỉ lệ thể tích hạt quá cao làm tăng mạnh hệ số dẫn nhiệt và mômen nhiệt uốn cục bộ, dẫn đến sự gia tăng mạnh của lực nhiệt trong tấm.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các đường cong phân bố nhiệt độ theo chiều sâu z, đồ thị bề mặt 3D mô tả trường độ võng w(x, y) trên không gian hai chiều, và bảng số liệu so sánh độ võng tại các điểm tọa độ đặc trưng. Kết quả so sánh giữa hai bài toán cho thấy ở mốc thời gian 1200 giây, đáp ứng cơ nhiệt của bài toán không dừng trùng khít với bài toán truyền nhiệt dừng với độ lệch dưới 0.5%, xác nhận tính tin cậy tuyệt đối của thuật toán phân tích.

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

Dựa trên các kết quả giải tích và mô phỏng số đã được kiểm chứng, nghiên cứu đưa ra 4 khuyến nghị kỹ thuật cụ thể:

Thứ nhất, kiểm soát chặt chẽ tỉ lệ phối trộn hạt gia cường Titan trong nền nhựa PVC ở mức tối ưu từ 15% đến 20% thể tích. Giải pháp này giúp các nhà máy chế tạo vật liệu đạt được độ cứng kháng uốn cực đại, ngăn ngừa nguy cơ nứt vỡ do sốc nhiệt và tiết kiệm khoảng 12% đến 15% chi phí hạt kim loại đắt tiền so với việc phối trộn vượt ngưỡng 25%.

Thứ hai, áp dụng quy trình kiểm soát gia nhiệt từ từ trong quá trình vận hành kết cấu composite, đảm bảo tốc độ tăng nhiệt độ bề mặt không vượt quá 1.5 K/phút trong khoảng thời gian khởi động 1200 giây đầu tiên. Biện pháp này triệt tiêu hoàn toàn ứng suất nhiệt quá độ phi tuyến ban đầu, bảo vệ bề mặt tấm khỏi hiện tượng biến dạng tức thời.

Thứ ba, tích hợp thuật toán giải tích nghiệm chuỗi Fourier kép vào các phần mềm tính toán kỹ thuật nội bộ tại các viện thiết kế tàu thủy và công trình dầu khí trong giai đoạn 2026 - 2028. Giải pháp này giúp rút ngắn hơn 70% thời gian phân tích tiền khả thi so với việc phải xây dựng mô hình phần tử hữu hạn phức tạp.

Thứ tư, thiết lập quy chuẩn kiểm định phân bố hạt định kỳ 6 tháng một lần bằng phương pháp siêu âm hoặc chụp cắt lớp tế vi đối với các tấm composite làm việc trong môi trường hóa chất ăn mòn, nhằm đảm bảo không xảy ra hiện tượng lắng đọng hay vón cục hạt Titan vượt mức cục bộ 25%.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn to lớn cho bốn nhóm đối tượng chuyên môn:

Nhóm 1: Các nghiên cứu sinh, học viên cao học và nhà khoa học thuộc ngành Cơ học vật thể rắn, Cơ kỹ thuật và Khoa học vật liệu. Tài liệu là nguồn tham khảo mẫu mực về kỹ thuật giải tích giải hệ phương trình đạo hàm riêng liên hợp nhiệt - cơ đàn hồi bằng chuỗi lượng giác và biến đổi Laplace.

Nhóm 2: Kỹ sư kết cấu và chuyên gia R&D trong ngành công nghiệp đóng tàu, dầu khí và đường ống dẫn hóa chất. Nghiên cứu cung cấp công thức tính toán trực tiếp độ võng và giới hạn chịu tải nhiệt của các tấm ốp vách, vỏ tàu và thành bể chứa composite polymer.

Nhóm 3: Các doanh nghiệp sản xuất vật liệu composite tiên tiến, linh kiện quang điện tử OLED và chip sinh học. Doanh nghiệp có thể ứng dụng bảng thông số cơ nhiệt để thiết lập tỷ lệ cấp phối vật liệu đạt độ bền cơ học cao nhất.

Nhóm 4: Giảng viên và sinh viên các trường đại học khối kỹ thuật như Đại học Bách khoa, Đại học Quốc gia. Luận văn đóng vai trò là tài liệu giảng dạy chuyên đề hữu ích về Cơ học vật liệu composite và Nhiệt đàn hồi ứng dụng.

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

Nguyên nhân nào khiến mặt giữa của tấm composite không xuất hiện biến dạng phẳng trong bài toán này? Theo chứng minh giải tích từ hệ phương trình vi phân chuyển vị, khi trường nhiệt độ chỉ biến thiên thuần túy theo phương chiều dày z và các điều kiện biên trên bốn cạnh là tựa bản lề đối xứng, lực nhiệt phát sinh trong mặt phẳng bị triệt tiêu hoàn toàn, dẫn đến các thành phần chuyển vị phẳng u và v tại mặt giữa đều đồng nhất bằng 0.

Tỉ lệ thể tích hạt Titan bao nhiêu là tối ưu nhất để hạn chế độ võng uốn nhiệt của tấm? Tỉ lệ phối trộn hạt Titan tối ưu nhất nằm trong khoảng từ 15% đến 20% thể tích. Tại dải nồng độ này, mô đun đàn hồi của composite tăng mạnh giúp độ võng giảm hơn 35% so với composite hạt loãng, trong khi vẫn duy trì được tính đồng nhất của trường nhiệt độ và ngăn chặn hiện tượng tập trung ứng suất nhiệt.

Sau bao lâu thì trạng thái truyền nhiệt không dừng chuyển sang trạng thái dừng ổn định? Với tấm composite dày 0.02 m và hệ số truyền nhiệt đối lưu khảo sát, khi thời gian gia nhiệt đạt mốc khoảng 1200 giây (tương đương 20 phút), sự phân bố nhiệt độ phi tuyến chuyển dịch hoàn toàn thành đường phân bố tuyến tính và độ võng của tấm đạt trạng thái ổn định, trùng khớp trên 99.5% với nghiệm truyền nhiệt dừng.

Tại sao phương pháp giải tích chuỗi Fourier kép lại ưu việt hơn phương pháp phần tử hữu hạn trong nghiên cứu này? Phương pháp giải tích cung cấp công thức nghiệm hiển minh dạng đại số, cho phép đánh giá trực tiếp và chính xác 100% mối quan hệ phụ thuộc giữa độ võng với các biến số hình học, tỷ lệ hạt và tải nhiệt mà không gặp phải sai số rời rạc hóa lưới hoặc tiêu tốn tài nguyên tính toán như phương pháp số.

Kết quả nghiên cứu này có thể ứng dụng thực tế vào những sản phẩm cụ thể nào tại Việt Nam? Kết quả được ứng dụng trực tiếp trong sản xuất vỏ tàu thuyền composite chịu nhiệt, ống dẫn dầu khí bọc hạt kim loại, tấm ốp cách nhiệt công trình biển, và các lớp đế nền vi mạch sinh học đòi hỏi độ giãn nở nhiệt cực thấp dưới tác động môi trường.

Kết luận

Luận văn đã hoàn thành xuất sắc các mục tiêu nghiên cứu với những đóng góp nổi bật:

  • Thiết lập thành công hệ phương trình cơ bản và điều kiện biên hoàn chỉnh cho bài toán uốn nhiệt tấm composite mỏng cốt hạt tựa bản lề bốn cạnh.
  • Xây dựng nghiệm giải tích giải tích giải phương trình truyền nhiệt dừng và không dừng bằng phép biến đổi Laplace kết hợp thuật toán chia đôi.
  • Tìm ra nghiệm hiển minh xác định độ võng của tấm dưới dạng chuỗi lượng giác Fourier kép phụ thuộc tường minh vào các tham số hình học và cơ nhiệt.
  • Xác định chính xác ngưỡng nồng độ hạt Titan tối ưu (15% - 20%) giúp tăng cường độ cứng kháng uốn và khả năng chịu nhiệt cho nền nhựa PVC.
  • Chứng minh tính tiệm cận hoàn hảo của bài toán truyền nhiệt không dừng về bài toán dừng khi thời gian đạt mốc 1200 giây.

Trong các giai đoạn nghiên cứu tiếp theo, hướng phát triển sẽ tập trung mở rộng mô hình cho kết cấu tấm và vỏ composite dị hướng nhiều lớp, composite chức năng FGM có cơ tính biến đổi liên tục, và khảo sát dưới tác động của trường tải trọng động lực học phức tạp. Độc giả và các kỹ sư quan tâm có thể ứng dụng trực tiếp hệ biểu thức nghiệm giải tích trong luận văn để nâng cao chất lượng thiết kế kết cấu composite ngay hôm nay.