Tổng quan nghiên cứu

Sự bùng nổ của công nghệ nano và công nghệ vi cơ điện tử đang định hình lại toàn bộ nền tảng kỹ thuật hiện đại, với tốc độ tăng trưởng thị trường toàn cầu dự kiến vượt 12% mỗi năm trong giai đoạn 2020 – 2030. Khi các cấu trúc kỹ thuật thu nhỏ về kích thước từ $10^{-6}$ mét (micrometer) đến $10^{-9}$ mét (nanometer), tỷ lệ diện tích bề mặt trên thể tích tăng vọt, kích hoạt các hiệu ứng cơ học phụ thuộc kích thước mà lý thuyết cơ học môi trường liên tục cổ điển không thể giải thích được. Các mô hình cổ điển của Euler – Bernoulli hay Navier thường đánh giá thấp độ cứng thực tế của kết cấu siêu nhỏ, dẫn đến sai số dự báo chuyển vị và ứng suất lên tới 30% đến trên 50%, tiềm ẩn nguy cơ phá hủy nghiêm trọng cho các thiết bị vi cơ điện tử (MEMS/NEMS).

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học và thuật toán phần tử hữu hạn tiên tiến dựa trên Lý thuyết độ dốc biến dạng (Strain Gradient Theory – SGT) nhằm mô phỏng chính xác đáp ứng tĩnh học của các phần tử thanh, dầm và hệ kết cấu khung/mạng lưới phẳng ở quy mô micro/nano. Nghiên cứu được triển khai thực hiện tại Trường Đại học Bách Khoa – ĐHQG TP.HCM trong khoảng thời gian từ tháng 09/2023 đến tháng 01/2024, tập trung vào phạm vi bài toán hai chiều với vật liệu đàn hồi tuyến tính chịu tải trọng tĩnh không đổi.

Về mặt ý nghĩa khoa học và thực tiễn, công trình cung cấp một công cụ tính toán số tin cậy, giúp tiết kiệm hơn 80% thời gian và tài nguyên tính toán so với phương pháp động lực học phân tử (Molecular Dynamics). Mô hình này đóng vai trò quan trọng trong việc tối ưu hóa thiết kế các cảm biến siêu nhạy, bộ dao động nano, và các dòng vật liệu nhân tạo metamaterials có cấu trúc tinh thể siêu nhẹ với độ bền cơ học vượt trộ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 vận dụng đồng thời ba khung lý thuyết cơ học môi trường liên tục: Lý thuyết độ dốc biến dạng (SGT), Lý thuyết ứng suất cặp đôi hiệu chỉnh (MCST), và Cơ học môi trường liên tục cổ điển. Trong đó, SGT là lý thuyết cơ học tổng quát nhất, kế thừa phát triển từ các công trình của Lam và Kahrobaiyan, tích hợp đồng thời ba tham số chiều dài vật liệu nội tại độc lập: tham số độ dốc giãn nở ($l_0$), tham số độ dốc kéo dài lệch ($l_1$), và tham số độ dốc góc xoay ($l_2$).

Năng lượng biến dạng tích lũy trong vật liệu không chỉ phụ thuộc vào tensor biến dạng cổ điển mà còn chứa đựng các thành phần năng lượng của vector độ dốc giãn nở, tensor độ dốc kéo dài lệch và tensor độ dốc góc xoay. Thông qua nguyên lý Hamilton và phương pháp biến phân, phương trình vi phân chủ đạo bậc 4 cho phần tử thanh và bậc 6 cho phần tử dầm được thiết lập chặt chẽ cùng hệ thống điều kiện biên phi cổ điển.

Năm khái niệm then chốt xuyên suốt mô hình bao gồm:

  • Tham số chiều dài đặc trưng vật liệu nội tại ($l_0, l_1, l_2$).
  • Tensor ứng suất bậc cao đại diện cho năng lượng biến dạng vi mô.
  • Hàm dạng giải tích chính xác suy xuất trực tiếp từ nghiệm phương trình vi phân thuần nhất.
  • Phần tử khung phẳng siêu vi tích hợp 5 bậc tự do tại mỗi nút.
  • Mô đun đàn hồi tương đương của cấu trúc mạng lưới tinh thể phẳng (Lattice Metamaterials).

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

Nguồn dữ liệu của nghiên cứu là dữ liệu số mô phỏng được phát triển qua các thuật toán phần tử hữu hạn tự lập trình, đối chiếu trực tiếp với các nghiệm giải tích và dữ liệu thực nghiệm đã công bố trên các tạp chí quốc tế uy tín của Lam (2003), Park (2008), và Kahrobaiyan (2013).

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu đã khảo sát toàn diện 7 mô hình kết cấu tiêu biểu bao gồm: thanh chịu kéo nén đúng tâm, dầm Euler – Bernoulli chịu uốn tĩnh, hệ khung cổng (portal frame), khung hình thoi, khung hình chữ nhật, khung mạng lục giác (honeycomb) và hệ khung mắt cáo (lattice frame). Quá trình khảo sát thực hiện trên hơn 50 kịch bản biến thiên thông số hình học, với tỷ số kích thước chiều cao trên chiều dài nội tại ($h/l$) thay đổi liên tục từ 0,5 đến 20 nhằm quét trọn vẹn dải chuyển tiếp từ thang nano sang thang macro.

Lý do lựa chọn phương pháp phân tích: Nghiên cứu lựa chọn phương pháp phần tử hữu hạn trực tiếp kết hợp hàm dạng giải tích chính xác (exact shape functions) thay cho các hàm dạng đa thức bậc thấp truyền thống. Giải pháp này loại bỏ hoàn toàn hiện tượng khóa số (locking phenomena), đảm bảo ma trận độ cứng đạt độ chính xác tuyệt đối mà không cần chia lưới phần tử quá mịn. Sự hội tụ và độ chính xác của trường chuyển vị được kiểm soát nghiêm ngặt bằng chuẩn sai số năng lượng Norm L2 với ngưỡng sai số kiểm định nhỏ hơn $10^{-6}$. Toàn bộ quy trình tính toán được thực hiện trong khung thời gian từ ngày 04/09/2023 đến 18/12/2023.

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

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

Phát hiện thứ nhất: Hiệu ứng tăng cứng cấu trúc (stiffening effect) thể hiện vô cùng mạnh mẽ ở quy mô micro/nano. Khi tỷ lệ kích thước tiết diện tiệm cận chiều dài nội tại của vật liệu ($h/l \approx 1$), độ cứng uốn và độ cứng dọc trục của dầm và khung tăng vọt. Kết quả tính toán cho thấy chuyển vị cực đại của dầm theo lý thuyết SGT giảm từ 35% đến 65% so với dự báo của lý thuyết cổ điển, khẳng định vật liệu thực tế cứng hơn nhiều so với các giả định truyền thống.

Phát hiện thứ hai: Tính tổng quát và khả năng suy biến hoàn hảo của phần tử SGT đề xuất. Khi cho các tham số chiều dài vật liệu $l_0, l_1, l_2$ dần tiến về 0, ma trận độ cứng của phần tử thanh (2 bậc tự do/nút), phần tử dầm (3 bậc tự do/nút) và phần tử khung (5 bậc tự do/nút) tự động suy giảm chính xác 100% về ma trận độ cứng theo lý thuyết MCST và lý thuyết cơ học cổ điển, với sai số số học dưới 0,01%.

Phát hiện thứ ba: Đáp ứng phụ thuộc kích thước của các hệ kết cấu khung phức tạp. Tại nút đỉnh của hệ khung hình thoi và khung cổng portal, chuyển vị ngang và góc xoay giảm từ 40% đến 58% khi kích thước thanh cấu thành giảm từ 100 micromet xuống 10 nanomet dưới cùng một tỷ lệ tải trọng chuẩn hóa.

Phát hiện thứ tư: Sự biến đổi của mô đun đàn hồi tương đương trong các cấu trúc mạng lưới. Đối với cấu trúc mạng lục giác và mạng mắt cáo, mô đun đàn hồi hiệu dụng tăng từ 1,5 đến 2,8 lần khi xét hiệu ứng gradient biến dạng so với khi tính toán bằng cơ học môi trường liên tục tiêu chuẩn, chứng minh tiềm năng vượt trội của các cấu trúc tinh thể nhân tạo ở cấp độ phân tử.

Thảo luận kết quả

Nguyên nhân vật lý dẫn đến hiện tượng tăng cứng là do ở quy mô nanometer, gradient biến dạng (đạo hàm bậc cao của trường biến dạng) xuất hiện với mật độ cực lớn tại các vùng biên và vị trí tập trung ứng suất. Năng lượng dự trữ trong các liên kết vi mô tạo ra hệ ứng suất bậc cao và các mômen cặp đôi chống lại sự uốn và giãn nở của kết cấu.

Khi so sánh với các nghiên cứu của Park (2008) và Kahrobaiyan (2013), kết quả của luận văn thể hiện sự tương đồng hoàn toàn về mặt định tính và định lượng, đồng thời mở rộng phạm vi áp dụng từ các thanh đơn lẻ sang các hệ kết cấu khung và mạng lưới hoàn chỉnh.

Về mặt hiển thị trực quan, các kết quả tính toán có thể được biểu diễn trực quan qua đồ thị đường cong tương quan giữa tỷ số chuyển vị chuẩn hóa $w_{SGT} / w_{Classical}$ theo tỷ số kích thước $h/l$. Đồ thị này minh họa một đường cong phi tuyến dốc đứng khi $h/l < 2$ và tiệm cận về đường nằm ngang giá trị 1 khi $h/l > 10$. Đồng thời, các bảng so sánh chuyển vị nút và góc xoay tại nút số 2 của khung portal giữa ba lý thuyết Cổ điển, MCST và SGT phản ánh rõ nét bước nhảy độ cứng khi kích thước hình học thu nhỏ.

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

Thứ nhất, chuẩn hóa quy trình phân tích và thiết kế hệ thống vi cơ điện tử MEMS/NEMS. Các doanh nghiệp công nghệ cao và viện nghiên cứu thiết bị bán dẫn cần tích hợp mô hình phần tử khung SGT 5 bậc tự do vào phần mềm mô phỏng kỹ thuật, nhằm giảm sai số dự báo tần số dao động và độ võng từ mức 40% hiện nay xuống dưới 5% trong giai đoạn 2025 – 2026.

Thứ hai, thiết lập cơ sở dữ liệu quốc gia về tham số chiều dài vật liệu nội tại. Bộ Khoa học và Công nghệ cùng các phòng thí nghiệm trọng điểm cần thực hiện các thí nghiệm uốn và xoắn vi mô trên thiết bị nano-indentation chuyên dụng, nhằm định lượng chính xác bộ ba tham số $l_0, l_1, l_2$ cho ít nhất 10 loại vật liệu nano phổ biến (như Silicon đơn tinh thể, Graphene, màng mỏng kim loại) trước năm 2027.

Thứ ba, mở rộng mô hình tính toán sang bài toán phi tuyến hình học và động lực học dao động. Các nhóm nghiên cứu chuyên sâu tại các trường đại học kỹ thuật cần tiếp tục nâng cấp thuật toán SGT để phân tích đáp ứng quá độ, hiện tượng mất ổn định đàn hồi (buckling) và dao động phi tuyến cho các hệ khung không gian 3 chiều trong lộ trình 12 đến 24 tháng tới.

Thứ tư, ứng dụng mô hình mạng lưới micro/nano vào chế tạo vật liệu siêu nhẹ (Metamaterials). Các kỹ sư cơ khí hàng không và công nghệ y sinh cần tận dụng công thức mô đun đàn hồi tương đương của khung lục giác và mắt cáo để thiết kế các khung xương nhân tạo và vỏ bọc chịu lực, giúp tăng khả năng hấp thụ năng lượng lên 30% trong khi giảm 25% khối lượng kết cấu giai đoạn 2026 – 2028.

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

Nhóm 1: Các nhà nghiên cứu và học viên cao học chuyên ngành Cơ học tính toán, Cơ kỹ thuật và Kỹ thuật xây dựng. Luận văn cung cấp toàn bộ quy trình thiết lập biến phân, nghiệm giải tích hàm dạng và ma trận độ cứng tổng thể của phần tử thanh, dầm, khung theo SGT, đóng vai trò là tài liệu tham khảo nền tảng để phát triển các phần tử tấm, vỏ vi mô bậc cao.

Nhóm 2: Kỹ sư thiết kế hệ thống vi cơ điện tử (MEMS/NEMS) và cảm biến nano. Giúp các kỹ sư nắm bắt chính xác quy luật phụ thuộc kích thước, từ đó tính toán chính xác độ võng, lực tĩnh điện kéo sập (pull-in voltage) và ứng suất phá hủy của các vi cần gạt (micro-cantilevers) và màng cảm biến áp suất.

Nhóm 3: Các chuyên gia phát triển phần mềm mô phỏng kỹ thuật (CAE). Cung cấp các công thức ma trận độ cứng giải tích tường minh và thuật toán chuyển đổi tọa độ địa phương - tổng thể, giúp lập trình viên dễ dàng tích hợp module phân tích cơ học phi cổ điển vào các phần mềm thương mại hoặc mã nguồn mở.

Nhóm 4: Các doanh nghiệp và kỹ sư nghiên cứu phát triển vật liệu tiên tiến (Advanced Materials). Hỗ trợ tính toán và tối ưu hóa các thông số hình học của mạng tinh thể nhân tạo, cấu trúc tổ ong và khung mắt cáo vi mô nhằm chế tạo vật liệu metamaterials có hệ số Poisson âm hoặc độ cứng riêng biệt cho ngành hàng không vũ trụ và cấy ghép y học.

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

Lý thuyết độ dốc biến dạng SGT khác biệt như thế nào so với cơ học môi trường liên tục cổ điển?

Lý thuyết SGT bổ sung thêm gradient của tensor biến dạng và ba tham số chiều dài vật liệu nội tại ($l_0, l_1, l_2$) vào hàm thế năng biến dạng. Trong khi cơ học cổ điển giả định vật liệu hoàn toàn đồng nhất và không xét đến kích thước tuyệt đối, SGT nắm bắt đầy đủ các hiệu ứng cấu trúc vi mô, giúp giải thích hiện tượng tăng cứng rõ rệt khi kích thước kết cấu nhỏ hơn 1 micromet.

Tại sao kết cấu ở tỉ lệ micro/nano lại có độ cứng cao hơn so với tính toán cổ điển?

Hiện tượng tăng cứng xuất phát từ mật độ gradient biến dạng cục bộ rất lớn khi kích thước hình học tiệm cận thang đo vi mô của mạng tinh thể. Các trường ứng suất bậc cao sinh ra các mômen cặp đôi kháng lại chuyển vị, làm cho độ cứng uốn thực tế của dầm nano có thể cao hơn gấp 2 đến 3 lần so với dự đoán của công thức Euler – Bernoulli truyền thống.

Phần tử khung SGT đề xuất trong luận văn có bao nhiêu bậc tự do tại mỗi nút?

Mỗi nút của phần tử khung SGT sở hữu 5 bậc tự do: chuyển vị dọc trục $u$, chuyển vị đứng $w$, biến dạng dọc trục $\varepsilon$, góc xoay $\theta$ và độ cong $\chi$. Cấu trúc này mở rộng vượt bậc so với 3 bậc tự do thông thường, tương ứng trực tiếp với 5 thành phần lực nút gồm lực dọc, lực cắt, mômen uốn và hai thành phần lực bậc cao.

Phương pháp phần tử hữu hạn SGT có ưu thế gì so với mô phỏng động lực học phân tử (MD)?

Mô phỏng động lực học phân tử đòi hỏi tài nguyên tính toán khổng lồ và bị giới hạn ở kích thước mẫu dưới vài triệu nguyên tử. Ngược lại, phương pháp phần tử hữu hạn SGT giải quyết các hệ khung mạng lưới phức tạp hàng nghìn phần tử chỉ trong vài giây trên máy tính cá nhân thông thường, tiết kiệm trên 90% thời gian mà vẫn đảm bảo độ chính xác cơ học liên tục.

Mô hình trong luận văn có thể áp dụng cho những dạng hình học kết cấu nào?

Luận văn đã kiểm chứng thành công trên 7 dạng hình học từ thanh, dầm đơn lẻ đến khung cổng portal, khung hình thoi, khung chữ nhật, khung lục giác tổ ong và hệ mạng mắt cáo. Thuật toán ma trận độ cứng tổng thể có tính khái quát cao, dễ dàng mở rộng cho các hệ dàn và mạng tinh thể không gian đa chiều.

Kết luận

  • Xây dựng thành công hệ phương trình vi phân chủ đạo và ma trận độ cứng giải tích chính xác cho phần tử thanh, dầm và khung phẳng dựa trên Lý thuyết độ dốc biến dạng (SGT).
  • Thiết lập phần tử khung phẳng vi mô với 5 bậc tự do tại mỗi nút, sử dụng hàm dạng giải tích chính xác giúp triệt tiêu hoàn toàn hiện tượng khóa số mà không cần tăng mật độ chia lưới.
  • Kiểm chứng độ tin cậy tuyệt đối của phương pháp với sai số năng lượng Norm L2 nhỏ hơn $10^{-6}$ và suy biến hoàn hảo về lý thuyết MCST và cổ điển khi triệt tiêu tham số kích thước.
  • Định lượng chi tiết hiệu ứng phụ thuộc kích thước trên 7 hệ cấu trúc phẳng, chứng minh mức độ tăng cứng làm giảm chuyển vị từ 35% đến 65% ở quy mô nanometer.
  • Khảo sát và thiết lập quy luật biến thiên của mô đun đàn hồi tương đương cho các cấu trúc mạng lục giác và mắt cáo, mở ra hướng ứng dụng thiết kế vật liệu metamaterials siêu bền nhẹ.

Về kế hoạch phát triển tiếp theo, nhóm nghiên cứu dự kiến hoàn thiện bộ phần tử tấm vỏ 3D và tích hợp phân tích phi tuyến hình học trong giai đoạn 2024 – 2025. Để khai thác và áp dụng toàn bộ các thuật toán ma trận độ cứng và mã nguồn phân tích kết cấu vi mô phục vụ cho các dự án nghiên cứu và phát triển sản phẩm MEMS/NEMS, hãy kết nối ngay với nhóm tác giả tại Trường Đại học Bách Khoa – ĐHQG TP.HCM để nhận tài liệu chuyên sâu và hỗ trợ kỹ thuật toàn diện.