Tổng quan nghiên cứu

Sự phát triển của công nghệ vật liệu nano trong những thập kỷ gần đây đã mở ra bước đột phá cho ngành kỹ thuật vật liệu, đặc biệt là sự xuất hiện của ống nano carbon (Carbon Nanotubes - CNT). Ống nano carbon sở hữu mô đun đàn hồi vượt trội lên đến xấp xỉ 1000 GPa cùng độ bền kéo cao gấp 50 đến 100 lần so với thép kết cấu thông thường, trong khi trọng lượng riêng chỉ bằng khoảng một phần sáu. Mặc dù vậy, việc khai thác hiệu quả cơ tính của CNT trong các hệ vật liệu nền polymer vẫn đối mặt với nhiều thách thức lớn. Vật liệu polymer thông thường có mô đun đàn hồi tương đối thấp, dao động trong khoảng 2 GPa đến 5 GPa, giới hạn khả năng chịu lực trong các kết cấu kỹ thuật cao.

Mục tiêu trọng tâm của nghiên cứu là xây dựng mô hình số hóa, tính toán và mô phỏng chi tiết các đặc trưng cơ học cơ bản của vật liệu nanocomposite nền polymer cốt ống nano carbon, bao gồm mô đun đàn hồi dọc trục, mô đun trượt và hệ số Poisson. Nghiên cứu tập trung giải quyết bài toán ứng xử cơ học tại vùng tiếp xúc pha giữa cốt CNT và nền polymer khi chịu các trạng thái tải trọng kéo và xoắn thuần túy.

Phạm vi nghiên cứu được thực hiện tại Trường Đại học Bách Khoa Hà Nội trong giai đoạn 2009 đến 2011, ứng dụng phương pháp phần tử hữu hạn trên phần tử thể tích đại diện (RVE) với hàm lượng thể tích CNT khảo sát từ 1% đến 5%. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp công cụ mô phỏng dự báo chính xác độ tăng cứng cơ học, giúp tiết kiệm tới 70% chi phí thử nghiệm mẫu thực nghiệm, đồng thời mở rộng tiềm năng ứng dụng vật liệu siêu nhẹ, độ bền cao cho ngành hàng không, vũ trụ và ô tô thế hệ mớ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 của cơ học môi trường liên tục và cơ học vật liệu composite đa thang. Để dự báo và đánh giá cơ tính của vật liệu nanocomposite, đề tài áp dụng hai mô hình giải tích kinh điển kết hợp mô hình tiếp xúc hiện đại:

Thứ nhất, quy luật hỗn hợp (Rule of Mixtures) được sử dụng làm chuẩn đối sánh lý thuyết giới hạn trên đối với mô đun đàn hồi dọc trục, xác định mối quan hệ trực tiếp giữa tỷ lệ thể tích và đặc trưng đàn hồi của các pha thành phần.

Thứ hai, mô hình bán thực nghiệm Halpin - Tsai được đưa vào để tính toán đặc tính đàn hồi của composite cốt sợi ngắn hoặc sợi liên tục, có xét đến tỷ lệ kích thước hình học (chiều dài trên đường kính) của ống nano carbon.

Thứ ba, lý thuyết vùng kết dính (Cohesive Zone Model) được tích hợp nhằm mô tả chính xác tương tác phi tuyến và sự suy giảm liên kết tại bề mặt ranh giới giữa ống CNT và chất nền polymer. Các khái niệm học thuật then chốt xuyên suốt đề tài bao gồm: ống nano carbon đơn vách (SWCNT), ống nano carbon đa vách (MWCNT), phần tử thể tích đại diện (Representative Volume Element - RVE), mô đun đàn hồi (E), mô đun trượt (G) và hệ số Poisson.

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

Nghiên cứu sử dụng phương pháp phân tích phần tử hữu hạn (Finite Element Method - FEM) thông qua phần mềm mô phỏng thương mại MSC.Marc kết hợp bộ tiền/hậu xử lý Mentat. Dữ liệu đầu vào về cơ tính của vật liệu nền polymer được chuẩn hóa với mô đun đàn hồi ban đầu đạt 5 GPa và hệ số Poisson bằng 0,30, trong khi ống nano carbon đơn vách SWCNT được gán mô đun đàn hồi dọc trục đạt xấp xỉ 1000 GPa.

Kích thước mẫu mô phỏng được xác định thông qua việc thiết lập phần tử thể tích đại diện RVE hình trụ có tính chất tuần hoàn, chứa ống nano carbon đặt đồng trục với thể tích nền xung quanh. Lý do lựa chọn mô hình RVE kết hợp FEM là vì phương pháp này cho phép kiểm soát chính xác tỷ lệ thể tích từ 1% đến 5%, thiết lập điều kiện biên đối xứng chặt chẽ và giải quyết bài toán tiếp xúc phi tuyến mà phương pháp giải tích thông thường không thể hiện thực hóa đầy đủ.

Quy trình mô phỏng thực hiện hai bài toán thử nghiệm ảo: thí nghiệm kéo đẳng tâm dọc trục để xác định mô đun đàn hồi và hệ số Poisson; thí nghiệm xoắn thuần túy thanh tiết diện tròn để xác định mô đun trượt. Toàn bộ chu trình tính toán được thực hiện theo tiến độ khoa học 24 tháng, bảo đảm kiểm chứng độc lập giữa mô hình phần tử hữu hạn và các công thức giải tích.

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

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

Nghiên cứu đã định lượng hóa sự gia tăng vượt trội về cơ tính của vật liệu composite nền polymer khi được gia cường bằng ống nano carbon qua các kịch bản liên kết giao diện khác nhau:

Thứ nhất, trong mô hình giả thiết liên kết lý tưởng giữa cốt CNT và nền, mô đun đàn hồi của vật liệu tăng đột biến. Khi hàm lượng thể tích CNT chỉ chiếm 1%, mô đun đàn hồi của nanocomposite tăng gấp 3 lần, từ mức ban đầu 5,0 GPa lên xấp xỉ 15,0 GPa (tương đương mức tăng 200%). Khi nâng tỷ lệ CNT lên 5%, mô đun đàn hồi tăng vọt hơn 10 lần, đạt giá trị 53,14 GPa (tăng trưởng 962,8% so với polymer nguyên bản). Hệ số Poisson trong mô hình lý tưởng duy trì ổn định ở mức tương đương hai vật liệu thành phần là xấp xỉ 0,30.

Thứ hai, trong mô hình liên kết sử dụng phần tử tiếp xúc (Cohesive Zone) phản ánh điều kiện trượt ranh giới thực tế, khả năng gia cường vẫn thể hiện rõ rệt nhưng ở mức độ vừa phải hơn. Cụ thể, với 1% CNT, mô đun đàn hồi của hệ vật liệu tăng gấp 2 lần, từ 5,0 GPa lên 10,0 GPa (tăng 100%). Khi tỷ lệ CNT đạt 5%, mô đun đàn hồi tăng gấp 6 lần, đạt giá trị xấp xỉ 30,0 GPa (tăng 500%).

Thứ ba, sự thay đổi của hệ số Poisson trong mô hình phần tử liên kết dao động trong khoảng từ 0,250 đến 0,285. Khi tỷ lệ phần trăm thể tích của CNT tăng từ 1% lên 5%, hệ số Poisson của vật liệu nanocomposite có xu hướng giảm đều khoảng 12,3%, minh chứng cho việc biến dạng ngang của cấu trúc bị triệt tiêu đáng kể khi vật liệu chịu kéo dọc trục.

Thứ tư, sự hội tụ giữa mô phỏng số và lý thuyết được kiểm chứng chặt chẽ: kết quả mô hình liên kết lý tưởng hoàn toàn trùng khớp với định luật hỗn hợp, trong khi kết quả từ mô hình vùng kết dính bám sát đường cong dự báo của mô hình Halpin - Tsai với độ lệch trung bình dưới 5%.

Thảo luận kết quả

Cơ chế gia cường của CNT bắt nguồn từ tỷ lệ diện tích bề mặt trên thể tích cực lớn và khả năng truyền ứng suất cắt vượt trội dọc theo thân ống nano. Khi chất nền polymer chịu lực kéo, ứng suất được phân tán nhanh chóng vào mạng lưới CNT có độ cứng siêu cao. Sự chênh lệch đáng kể giữa mô hình liên kết lý tưởng (53,14 GPa tại 5% CNT) và mô hình vùng kết dính (30,0 GPa tại 5% CNT) phản ánh rõ nét hiện tượng suy giảm ứng suất cắt do trượt cục bộ tại bề mặt phân cách pha polymer/CNT.

Dữ liệu tính toán trong luận văn có thể được trực quan hóa thông qua bảng so sánh đa thông số và đồ thị phân bố ứng suất von Mises dọc theo chiều dài phần tử RVE. Đồ thị biểu diễn mối quan hệ giữa hàm lượng thể tích CNT và mô đun đàn hồi cho thấy một hàm tăng phi tuyến, trong đó đường cong mô phỏng phần tử liên kết nằm giữa đường biên trên của định luật hỗn hợp và đường thực nghiệm hiệu chỉnh Halpin - Tsai. Điều này khẳng định tính đúng đắn của việc đưa vùng kết dính Cohesive Zone vào phân tích kết cấu nano, cung cấp cơ sở tin cậy cho việc thiết kế vật liệu composite có độ chính xác cao.

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

Dựa trên các kết quả tính toán số trị, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật nhằm tối ưu hóa việc chế tạo và ứng dụng vật liệu nanocomposite polymer/CNT trong thực tế:

Thứ nhất, tăng cường xử lý hóa tính bề mặt của ống nano carbon (Surface Functionalization) bằng các nhóm chức hóa học phù hợp nhằm nâng cao độ bền liên kết tiếp xúc giữa cốt và nền. Mục tiêu đặt ra là gia tăng ứng suất cắt ranh giới tiếp xúc lên trên 50 MPa, hạn chế tối đa hiện tượng trượt bề mặt trong vòng 12 tháng nghiên cứu tiếp theo do các phòng thí nghiệm vật liệu đảm nhiệm.

Thứ hai, mở rộng mô hình tính toán đa thang kết hợp động lực học phân tử (Molecular Dynamics) và phần tử hữu hạn đa cấu trúc nhằm dự báo chính xác các hiệu ứng bất đẳng hướng và sự uốn lượn tự nhiên của CNT trong nền polymer. Mục tiêu kiểm soát sai số mô phỏng dưới mức 3% đối với các phân bố không gian ngẫu nhiên, thực hiện trong khung thời gian 18 tháng bởi nhóm nghiên cứu cơ học tính toán.

Thứ ba, hoàn thiện công nghệ khuấy trộn siêu âm phân tán năng lượng cao kết hợp dung môi bay hơi trong quy trình đúc composite, nhằm ngăn chặn 95% hiện tượng vón cục của CNT khi tăng tỷ lệ thể tích vượt quá 3%, do các kỹ sư quy trình vật liệu triển khai trong thời hạn 6 đến 9 tháng.

Thứ tư, xây dựng bộ tiêu chuẩn đánh giá và cơ sở dữ liệu số hóa về đặc trưng cơ tính nanocomposite cho tối thiểu 10 loại nền polymer thông dụng (như Epoxy, PEEK, Polypropylene) gia cường nano carbon, phục vụ thiết kế cơ khí kết cấu. Đề án tiêu chuẩn hóa này cần được hoàn thành trong vòng 24 tháng dưới sự chủ trì của các trường đại học kỹ thuật và cơ quan quản lý tiêu chuẩn đo lường chất lượng.

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

Luận văn là tài liệu tham khảo học thuật và kỹ thuật giá trị dành cho các nhóm độc giả sau:

Thứ nhất, học viên cao học, nghiên cứu sinh và giảng viên thuộc các chuyên ngành Cơ học kỹ thuật, Khoa học vật liệu và Cơ kỹ thuật tính toán. Luận văn cung cấp phương pháp luận chuẩn mực về xây dựng phần tử thể tích đại diện RVE và kỹ thuật mô phỏng tiếp xúc trên phần mềm MSC.Marc.

Thứ hai, kỹ sư nghiên cứu và phát triển (R&D) trong các ngành công nghiệp hàng không, vũ trụ, thiết bị quân sự và sản xuất ô tô. Nhóm đối tượng này có thể ứng dụng các bộ thông số cơ tính để thiết kế vỏ thân máy bay, cánh quạt tuabin nhẹ hơn từ 20% đến 40% so với kết cấu kim loại truyền thống.

Thứ ba, chuyên gia phân tích kết cấu CAE và tính toán mô phỏng phần tử hữu hạn (FEA Engineers). Tài liệu mang lại hướng tiếp cận rõ ràng trong việc thiết lập điều kiện biên tuần hoàn và mô hình hóa vùng kết dính Cohesive Zone trong các bài toán tiếp xúc phức tạp.

Thứ tư, doanh nghiệp sản xuất vật liệu composite tiên tiến, vật liệu thể thao và linh kiện điện tử tản nhiệt cao, giúp hoạch định tỷ lệ pha trộn phụ gia CNT tối ưu theo bài toán chi phí và hiệu năng cơ học.

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

Ống nano carbon giúp cải thiện độ cứng của polymer như thế nào? Ống nano carbon sở hữu mô đun đàn hồi xấp xỉ 1000 GPa, khi phân tán vào nền polymer sẽ tạo thành mạng lưới truyền lực chịu ứng suất chính. Ngay cả với hàm lượng thấp từ 1% đến 5%, vật liệu nanocomposite đã có thể gia tăng mô đun đàn hồi từ 2 đến 10 lần so với mức 5 GPa của polymer nguyên sinh.

Vì sao kết quả của mô hình phần tử liên kết lại thấp hơn mô hình liên kết lý tưởng? Mô hình liên kết lý tưởng giả định sự dính bám tuyệt đối không trượt tại ranh giới pha. Trong khi đó, mô hình phần tử liên kết xét đến sự biến dạng cắt, hiện tượng trượt vi mô và suy giảm truyền lực tại vùng giao tiếp phân tử, phản ánh chính xác các giới hạn vật lý của tương tác Van der Waals trong thực tế.

Hàm lượng ống nano carbon bao nhiêu phần trăm là tối ưu cho vật liệu composite? Nghiên cứu cho thấy tỷ lệ thể tích từ 1% đến 3% mang lại hiệu quả gia cường kinh tế và kỹ thuật cao nhất. Ở nồng độ 1%, mô đun đàn hồi đã tăng 100% đối với mô hình kết dính, trong khi vượt quá 5% vật liệu dễ gặp hiện tượng vón cục hạt nano trong thực tế chế tạo.

Phương pháp phần tử hữu hạn trên MSC.Marc có đáng tin cậy trong cơ học nano không? Phần mềm MSC.Marc có khả năng xử lý bài toán phi tuyến hình học, phi tuyến vật liệu và phi tuyến tiếp xúc rất mạnh mẽ. Kết quả mô phỏng RVE kéo và xoắn trong đề tài có độ sai lệch dưới 5% khi đối chiếu với mô hình giải tích Halpin - Tsai và quy luật hỗn hợp đã được thừa nhận rộng rãi.

Ý nghĩa của việc hệ số Poisson giảm từ 0,285 xuống 0,250 khi tăng tỷ lệ CNT là gì? Hệ số Poisson giảm phản ánh cấu trúc nanocomposite trở nên cứng vững hơn theo phương ngang khi chịu kéo dọc trục. Sự hiện diện của mạng lưới ống nano hạn chế sự co thắt tiết diện của chất nền polymer, giúp kết cấu duy trì độ ổn định hình học tốt hơn dưới tải trọng kéo lớn.

Kết luận

Nghiên cứu đã hoàn thành xuất sắc các mục tiêu khoa học đề ra với những đóng góp quan trọng:

  • Xây dựng thành công mô hình phần tử thể tích đại diện RVE chuẩn hóa cho vật liệu nanocomposite nền polymer cốt ống nano carbon chịu tải kéo và xoắn.
  • Định lượng hóa mức tăng mô đun đàn hồi từ 2 lần đến 10 lần khi thay đổi tỷ lệ thể tích CNT từ 1% đến 5% trên hai mô hình liên kết lý tưởng và vùng kết dính.
  • Làm sáng tỏ quy luật biến thiên của hệ số Poisson trong khoảng 0,250 đến 0,285 dưới tác động của hiệu ứng cản trở biến dạng ngang từ ống CNT.
  • Chứng minh tính tương thích và hội tụ cao giữa phương pháp mô phỏng số MSC.Marc với mô hình giải tích Halpin - Tsai và định luật hỗn hợp.
  • Thiết lập cơ sở lý thuyết và quy trình tính toán số phục vụ định hướng chế tạo thực nghiệm vật liệu composite nano tính năng cao.

Các nhóm nghiên cứu và đơn vị phát triển vật liệu có thể khai thác ngay khung phương pháp mô phỏng này để rút ngắn chu kỳ thử nghiệm sản phẩm mới trong giai đoạn 12 đến 24 tháng tới, thúc đẩy ứng dụng vật liệu cấu trúc nano tiên tiến vào các ngành công nghiệp mũi nhọn.