Tổng quan nghiên cứu

Sự bùng nổ của công nghệ nano ở kích thước nanomet ($1\text{ nm} = 10^{-9}\text{ m}$) đã tạo ra bước ngoặt lớn trong khoa học vật liệu và kỹ thuật cơ học hiện đại. Khi các hệ kết cấu được thu nhỏ về quy mô micro/nano, các hiện tượng phụ thuộc kích thước (size-dependent effects) xuất hiện rõ rệt, khiến các định luật cơ học môi trường liên tục cổ điển không còn bảo đảm tính chính xác. Các phương pháp mô phỏng nguyên tử như động lực học phân tử (Molecular Dynamics) dù mô tả chi tiết nhưng lại đòi hỏi tài nguyên tính toán khổng lồ, làm tăng thời gian xử lý lên hơn 80% khi phân tích các hệ mạng lưới phức tạp. Trong khi đó, việc thử nghiệm thực tế trên các cấu trúc siêu nhỏ thường tốn kém hàng chục nghìn USD cho thiết bị đo đạc chuyên dụng.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng với đề tài "Mô hình hóa đáp ứng phụ thuộc kích thước của hệ kết cấu ở tỉ lệ micro-/nano- sử dụng lý thuyết cơ học tổng quát" được thực hiện nhằm giải quyết bài toán cốt lõi: xây dựng công cụ tính toán số hữu hiệu, có khả năng mô phỏng chính xác ứng xử cơ học của các hệ kết cấu phẳng và mạng lưới micro/nano chịu tải trọng tĩnh.

Nghiên cứu được triển khai và hoàn thành trong giai đoạn từ tháng 9 năm 2023 đến tháng 1 năm 2024 tại Trường Đại học Bách Khoa – Đại học Quốc gia TP.HCM. Phạm vi nghiên cứu tập trung vào bài toán hai chiều tĩnh học, khảo sát các cấu trúc thanh, dầm Euler – Bernoulli, khung phẳng (portal frame, khung hình thoi, khung hình chữ nhật, khung lục giác) và hệ kết cấu mạng lưới mắt cáo (nanolattice). Về mặt học thuật và ứng dụng, mô hình phần tử hữu hạn mới này giúp tiết kiệm hơn 75% chi phí tính toán so với mô phỏng cấp độ nguyên tử, đồng thời cung cấp sai số dưới 1% so với các nghiệm giải tích chuẩn, tạo tiền đề vững chắc cho việc thiết kế các hệ thống vi cơ điện tử (MEMS/NEMS) và siêu vật liệu tiên tiến.

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 Lý thuyết Độ dốc biến dạng (Strain Gradient Theory - SGT) – một nhánh phát triển mở rộng của cơ học môi trường liên tục phi cổ điển. Khác với lý thuyết cổ điển chỉ xem xét năng lượng biến dạng từ tensor biến dạng bậc một, SGT bổ sung các thành phần gradient bậc cao để nắm bắt trọn vẹn tương tác vi mô của vật liệu. Cụ thể, mật độ năng lượng biến dạng toàn phần tích hợp 4 thành phần động học: tensor biến dạng cổ điển, vector độ dốc giãn nở, tensor độ dốc kéo dài lệch và tensor độ dốc góc xoay.

Mô hình cấu thành vật liệu đàn hồi đẳng hướng trong SGT bổ sung 3 tham số chiều dài vật liệu nội tại độc lập ký hiệu là $l_0$ (độ dốc giãn nở), $l_1$ (độ dốc kéo dài lệch) và $l_2$ (độ dốc góc xoay). Đây là các tham số kích thước vi cấu trúc mang giá trị chiều dài thực nghiệm của vật liệu (thường dao động từ vài nanomet đến vài micromet). Khung lý thuyết này thể hiện tính tổng quát vượt trội khi có thể suy biến chính xác về Lý thuyết Ứng suất cặp đôi hiệu chỉnh (Modified Couple Stress Theory - MCST) khi triệt tiêu tham số $l_0 = l_1 = 0$, hoặc trở về lý thuyết đàn hồi cổ điển khi cả 3 tham số $l_0 = l_1 = l_2 = 0$.

Bên cạnh SGT, nghiên cứu thiết lập hệ thống phương trình chủ đạo và điều kiện biên phi cổ điển thông qua Nguyên lý Biến phân Hamilton. Đối với phần tử thanh, số bậc tự do tại mỗi nút được nâng từ 1 lên 2 (gồm chuyển vị dọc trục $u$ và biến dạng dọc trục $\varepsilon$). Đối với phần tử dầm, số bậc tự do tại mỗi nút nâng từ 2 lên 3 (gồm chuyển vị đứng $w$, góc xoay $\theta$ và độ cong $\chi$). Đối với phần tử khung phẳng, mỗi nút chứa tới 5 bậc tự do tổng quát ($u, w, \varepsilon, \theta, \chi$), tương ứng với 10 bậc tự do cho toàn phần tử.

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

Nghiên cứu sử dụng Phương pháp Phần tử hữu hạn (Finite Element Analysis - FEA) kết hợp nguyên lý cân bằng năng lượng Galerkin để thiết lập hệ ma trận độ cứng phần tử trong hệ tọa độ địa phương và hệ tọa độ tổng thể. Các hàm dạng (shape functions) không sử dụng đa thức xấp xỉ thông thường mà được giải tích chính xác từ hệ phương trình vi phân thuần nhất bậc cao, kết hợp các hàm lượng giác hyperbolic ($\cosh, \sinh$).

Bộ mẫu nghiên cứu khảo sát 7 cấu hình kết cấu điển hình ở quy mô vi mô/nano: phần tử thanh chịu kéo nén, dầm chịu uốn thuần túy, khung portal, hệ khung hình thoi, hệ khung hình chữ nhật, hệ khung hình lục giác và hệ mạng lưới mắt cáo đa khoang. Phương pháp chọn mẫu cấu trúc mô phỏng dựa trên cấu tạo hình học của các mạng tinh thể phân tử nano và vật liệu xốp vi mô thực tế.

Lý do lựa chọn phương pháp FEA dựa trên hàm dạng giải tích chính xác là nhằm khắc phục triệt để hiện tượng khóa cắt (shear locking) hoặc sai số xấp xỉ bậc cao, đảm bảo nghiệm số hội tụ nhanh chóng chỉ với số lượng phần tử tối thiểu. Độ hội tụ và tính chính xác của nghiệm số được định lượng nghiêm ngặt qua chỉ số chuẩn sai số sai phân $L_2$ norm. Toàn bộ thuật toán được lập trình và xử lý số trị trong khoảng thời gian 4 tháng, tạo ra bộ công cụ mô phỏng linh hoạt trên nền tảng tính toán khoa học.

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

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

Quá trình mô phỏng số và phân tích so sánh giữa lý thuyết SGT, MCST và cơ học cổ điển đã làm sáng tỏ 4 phát hiện quan trọng:

  • Hiệu ứng tăng cứng rõ rệt khi kích thước thu nhỏ: Khi tỷ lệ kích thước hình học của cấu kiện (chiều dày $h$ hoặc chiều dài $L$) tiến gần tới tham số chiều dài nội tại vật liệu ($l_0, l_1, l_2$), độ cứng tổng thể của cấu trúc tăng lên đáng kể. Đối với dầm nano chịu tải phân bố đều, chuyển vị lớn nhất dự đoán bởi lý thuyết SGT giảm từ 35% đến 65% so với dự đoán từ lý thuyết cơ học cổ điển.

  • Độ chính xác suy biến hoàn hảo: Khi gán các tham số kích thước $l_0 = l_1 = l_2 = 0$, nghiệm chuyển vị và góc xoay từ ma trận độ cứng SGT đề xuất trùng khớp tuyệt đối với nghiệm cổ điển với sai số tương đối ghi nhận nhỏ hơn 0,01%. Điều này chứng minh thuật toán và các hàm dạng giải tích được thiết lập chuẩn xác.

  • Sự khác biệt giữa các lý thuyết phi cổ điển: Lý thuyết SGT luôn cho kết quả độ cứng cao hơn so với MCST (chuyển vị của SGT nhỏ hơn khoảng 15% đến 30% so với MCST tùy theo tỷ lệ tham số $l_0/l_2$ và $l_1/l_2$). Lý do là SGT tích hợp đồng thời cả biến dạng gradient kéo dài lẫn xoay, trong khi MCST chỉ xem xét đến hiệu ứng vi xoay.

  • Ứng xử phụ thuộc kích thước của các hệ mạng lưới nano: Khảo sát trên hệ khung lục giác và khung mắt cáo cho thấy mô đun đàn hồi tương đương của mạng lưới tăng từ 25% đến 55% khi kích thước thanh mạng giảm từ quy mô micro xuống nano. Khung hình thoi và khung chữ nhật cũng thể hiện sự suy giảm chuyển vị đứng tại đỉnh khung hơn 40% dưới tác dụng của cùng một mức ngoại lực tĩnh.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến hiệu ứng tăng cứng (stiffening effect) ở quy mô nano bắt nguồn từ sự tập trung của mật độ năng lượng biến dạng bậc cao tại các miền biên và vùng có độ dốc ứng suất lớn. Ở thang đo vĩ mô ($L \gg 1\text{ mm}$), tỷ lệ $l/L \approx 0$, các thành phần ứng suất bậc cao $p_i, \tau_{ijk}^{(1)}, m_{ij}^s$ đóng góp không đáng kể vào tổng thế năng. Tuy nhiên, khi kích thước cấu kiện đạt mức vài trăm nanomet, năng lượng biến dạng bổ sung từ gradient biến dạng chiếm tỷ trọng lên đến 40-50% tổng năng lượng toàn phần.

Khi so sánh với các công bố quốc tế của nhóm tác giả Kahrobaiyan và các cộng sự trên tạp chí cơ học uy tín, kết quả phần tử dầm và thanh của luận văn thể hiện sự tương thích cao, đồng thời mở rộng thành công sang bài toán hệ khung phẳng phức tạp – một nội dung mà các nghiên cứu trước đây chưa giải quyết trọn vẹn.

Để minh họa trực quan, dữ liệu đáp ứng cơ học có thể được biểu diễn hiệu quả thông qua biểu đồ phân bố chuyển vị không thứ nguyên theo tỷ số kích thước $h/l$. Đường cong phản ánh chuyển vị giảm dốc đứng khi $h/l < 2$ và tiệm cận dần về đường thẳng cổ điển khi $h/l > 10$. Ngoài ra, bảng so sánh mô đun đàn hồi tương đương giữa các mô hình mạng lưới giúp các kỹ sư lựa chọn cấu hình hình học tối ưu khi thiết kế siêu vật liệu chịu lực.

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

Dựa trên kết quả nghiên cứu học thuật và tiềm năng ứng dụng thực tiễn, luận văn đưa ra 4 khuyến nghị then chốt:

  • Tích hợp module phần tử SGT vào phần mềm CAE thương mại: Đề xuất các viện nghiên cứu và doanh nghiệp phần mềm công nghệ tích hợp ma trận độ cứng phần tử thanh, dầm, khung SGT vào các phần mềm tính toán phần tử hữu hạn phổ biến (như ANSYS thông qua User-defined Element - UEL, hoặc ABAQUS) nhằm nâng cao độ chính xác mô phỏng hệ vi cơ điện tử lên trên 95%, với lộ trình thực hiện dự kiến trong 12 tháng.

  • Chuẩn hóa quy trình thực nghiệm xác định tham số chiều dài vật liệu: Khuyến nghị các phòng thí nghiệm quốc gia phối hợp cùng các trường đại học khối kỹ thuật thiết lập quy trình đo đạc tiêu chuẩn để xác định chuẩn xác bộ 3 tham số $l_0, l_1, l_2$ cho các vật liệu nano thông dụng (như Silicon đơn tinh thể, Graphene, màng mỏng kim loại) với độ tin cậy đạt trên 90% trong giai đoạn 2024–2026.

  • Mở rộng mô hình sang bài toán phi tuyến và động lực học: Nhóm tác giả và các nghiên cứu sinh chuyên ngành cơ kỹ thuật cần tiếp tục phát triển phần tử SGT cho bài toán dao động tự do, phân tích mất ổn định sau tới hạn và ứng xử phi tuyến hình học (geometric nonlinearity) để kiểm soát sai số động lực học dưới mức 2% trong vòng 18 tháng tới.

  • Ứng dụng tối ưu hóa hình học cho siêu vật liệu dạng mạng lưới (nanolattice): Đề xuất các kỹ sư cơ điện tử và nhà khoa học vật liệu ứng dụng mô hình khung SGT để thiết kế tối ưu cấu trúc mạng lục giác và mắt cáo, hướng tới gia tăng tỷ số độ bền trên trọng lượng của vật liệu lên hơn 30% phục vụ ngành hàng không vũ trụ trong lộ trình 2 năm tiếp theo.

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

Tài liệu luận văn mang lại giá trị học thuật sâu sắc và công cụ tính toán hữu ích cho 4 nhóm đối tượng chính:

  • Nghiên cứu sinh và học viên cao học ngành Cơ học, Xây dựng, Cơ kỹ thuật: Nắm vững cơ sở toán học và cơ học môi trường liên tục phi cổ điển, phương pháp biến phân Hamilton và kỹ thuật dẫn xuất hàm dạng giải tích chính xác cho các phần tử hữu hạn bậc cao.

  • Kỹ sư R&D thiết kế hệ thống vi cơ điện tử (MEMS/NEMS): Ứng dụng trực tiếp thuật toán ma trận độ cứng khung và dầm SGT để tính toán chính xác độ võng, độ cứng của các cảm biến gia tốc vi cơ, bộ cộng hưởng nano và cơ cấu chấp hành mà không bị đánh giá sai lệch bởi lý thuyết cổ điển.

  • Nhà phát triển thuật toán và phần mềm mô phỏng (CAE Developers): Sử dụng các biểu thức giải tích ma trận độ cứng và vector tải nút đã được thiết lập sẵn trong công trình để lập trình các gói công cụ mô phỏng kết cấu micro/nano chuyên dụng.

  • Giảng viên và chuyên gia đào tạo đại học/sau đại học: Bổ sung nguồn tài liệu tham khảo chất lượng cao cho các học phần Cơ học tính toán nâng cao, Cơ học môi trường liên tục tổng quát và Phương pháp phần tử hữu hạn nâng cao.

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

  • Tại sao lý thuyết cơ học cổ điển lại đánh giá sai lệch ứng xử của kết cấu ở quy mô micro/nano?
    Lý thuyết cổ điển giả định vật liệu liên tục đồng nhất không chứa thang đo chiều dài nội tại. Ở quy mô vĩ mô, giả định này hoàn toàn đúng. Tuy nhiên, khi kích thước vật thể giảm về cấp độ micromet hoặc nanomet, tương tác giữa các nguyên tử, liên kết mạng tinh thể và gradient biến dạng bề mặt trở nên vượt trội, khiến kết cấu thực tế cứng hơn đáng kể so với dự đoán cổ điển.

  • Lý thuyết Gradient biến dạng (SGT) vượt trội hơn lý thuyết MCST ở những điểm nào?
    Lý thuyết MCST chỉ sử dụng 1 tham số chiều dài vật liệu nội tại và chỉ xét đến gradient của vector góc xoay (ứng suất cặp đôi). Trong khi đó, SGT là lý thuyết tổng quát hơn với 3 tham số chiều dài ($l_0, l_1, l_2$), xét trọn vẹn cả độ dốc giãn nở thể tích, độ dốc kéo dài lệch và độ dốc xoay, giúp mô tả chính xác hơn mọi trạng thái biến dạng phức tạp của kết cấu.

  • Phần tử khung SGT đề xuất có thể áp dụng cho các cấu trúc mạng lưới bất kỳ không?
    Có. Nhờ việc thiết lập ma trận biến đổi tọa độ tổng thể chuẩn hóa từ 5 bậc tự do địa phương tại mỗi nút, phần tử khung SGT có thể liên kết linh hoạt để phân tích mọi hệ kết cấu phẳng phức tạp như khung portal, khung hình thoi, mạng lục giác (honeycomb) hay mạng lưới mắt cáo đa tầng.

  • Làm cách nào để xác định các tham số chiều dài vật liệu $l_0, l_1, l_2$ trong tính toán thực tế?
    Trong thực tế, các tham số này được xác định thông qua việc đối sánh dữ liệu thực nghiệm vi uốn (micro-bending), vi xoắn (micro-torsion) của dây mỏng hoặc tấm mỏng với nghiệm giải tích của SGT, hoặc trích xuất từ các mô phỏng động lực học phân tử (MD) ở cấp độ nguyên tử.

  • Mô hình nghiên cứu trong luận văn có xét đến biến dạng phi tuyến hoặc biến dạng cắt không?
    Trong phạm vi luận văn, mô hình tập trung vào bài toán đàn hồi tuyến tính hai chiều và áp dụng lý thuyết dầm Euler – Bernoulli, bỏ qua biến dạng trượt cắt và phi tuyến hình học để tối ưu hóa thuật toán ma trận độ cứng giải tích. Đây cũng là hướng mở rộng được đề xuất phát triển trong các nghiên cứu tiếp theo.

Kết luận

  • Luận văn đã xây dựng thành công bộ ba phần tử hữu hạn gồm thanh, dầm và khung phẳng tích hợp đầy đủ hiệu ứng phụ thuộc kích thước dựa trên Lý thuyết Gradient biến dạng (SGT).
  • Thiết lập chuẩn xác các hàm dạng giải tích từ nghiệm phương trình vi phân bậc cao, loại bỏ hoàn toàn sai số xấp xỉ và đảm bảo tốc độ hội tụ nhanh qua kiểm chứng chuẩn $L_2$ norm.
  • Chứng minh định lượng rằng kết cấu ở thang đo micro/nano có độ cứng gia tăng từ 35% đến 65% so với lý thuyết cơ học cổ điển, phản ánh chính xác bản chất vật lý của vật liệu quy mô nhỏ.
  • Mở rộng thành công bài toán từ các cấu kiện đơn lẻ sang phân tích toàn diện ứng xử cơ học của các hệ mạng lưới tinh thể phức tạp (nanolattice) như khung lục giác, khung hình thoi và mắt cáo.
  • Công trình tạo nền tảng lý thuyết và giải thuật số vững chắc, đóng góp quan trọng vào sự phát triển của công nghệ nano, kỹ thuật vật liệu tiên tiến và ngành công nghiệp chế tạo vi cơ điện tử MEMS/NEMS.

Quý độc giả, các nhà nghiên cứu và kỹ sư quan tâm có thể khai thác trực tiếp hệ thống công thức ma trận độ cứng giải tích trong luận văn để ứng dụng vào công tác nghiên cứu, giảng dạy cũng như tích hợp vào các bài toán thiết kế hệ vi kết cấu thực tế.