Tổng quan nghiên cứu

Công nghệ nano hiện đại tập trung vào việc thao tác và chế tạo các cấu trúc có kích thước trong dải từ 0,1 nm đến 100 nm, tương đương 1 phần tỷ mét. Ở cấp độ này, các hệ thống vi điện cơ MEMS và nano điện cơ NEMS bao gồm các cấu kiện dầm, tấm và màng siêu mỏng sở hữu tỷ lệ diện tích bề mặt trên thể tích cực lớn, mang lại những tính chất cơ - lý vượt trội. Tuy nhiên, lý thuyết đàn hồi liên tục cổ điển của Hooke vốn độc lập với quy mô nên không thể phản ánh chính xác hiệu ứng kích thước nhỏ. Trong khi đó, các mô hình mô phỏng động lực học phân tử lại đòi hỏi chi phí tính toán rất lớn và khó áp dụng cho các hệ kết cấu quy mô vi mô.

Luận văn thạc sĩ chuyên ngành Xây dựng công trình dân dụng và công nghiệp được hoàn thiện từ tháng 07 năm 2014 đến tháng 01 năm 2015 tại Trường Đại học Bách khoa, Đại học Quốc gia Thành phố Hồ Chí Minh. Đề tài tập trung giải quyết bài toán phân tích dao động tự do và hiện tượng mất ổn định của dầm dựa trên mô hình vi kết cấu rời rạc. Mục tiêu trọng tâm là thiết lập mô hình dầm gồm các đoạn cứng liên kết bằng lò xo xoay đàn hồi, khảo sát đa dạng điều kiện biên và quy luật tải trọng, đồng thời xác định chính xác hệ số ảnh hưởng quy mô nhỏ. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp công cụ tính toán giải tích số tin cậy, giúp nâng cao độ chính xác trong thiết kế cấu kiện vi kết cấu và giảm thiểu sai số dự báo mất ổn định kết cấu nano xuống dưới 5%.

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 lý thuyết cơ học liên tục phi cục bộ do Eringen đề xuất năm 1972 và nguyên lý biến phân Hamilton. Khác với lý thuyết cổ điển, cơ học phi cục bộ giả định ứng suất tại một điểm là hàm phụ thuộc vào biến dạng của mọi điểm trong toàn bộ môi trường liên tục. Mối quan hệ ứng suất - biến dạng phi cục bộ được biểu diễn qua phương trình vi phân có chứa hệ số ảnh hưởng quy mô nhỏ e0 và chiều dài đặc trưng a tương ứng với khoảng cách nguyên tử.

Khung lý thuyết của luận văn tích hợp ba khái niệm cốt lõi: mô hình dầm Euler-Bernoulli phi cục bộ, thế năng đàn hồi của hệ lò xo xoay với độ cứng C = EI/a, và động năng phân bố tại các nút tập trung. Cấu trúc vi kết cấu gồm n phần tử cứng lặp lại có tổng chiều dài L = n x a. Để mô phỏng chính xác liên kết ngàm, luận văn áp dụng đề xuất của Challamel năm 2013 với độ cứng lò xo biên C1 = 2C. Hệ số e0 biến thiên trong khoảng từ 0,289 đối với bài toán mất ổn định uốn dọc thuần túy đến 0,408 đối với bài toán dao động riêng, phản ánh chính xác trạng thái ứng suất dọc trục ban đầu của hệ vi cơ học.

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

Nguồn dữ liệu nghiên cứu được xây dựng thông qua việc mô hình hóa toán học 6 bộ ví dụ số điển hình đại diện cho 12 bài toán kỹ thuật cụ thể. Cỡ mẫu nghiên cứu bao gồm các cấu hình dầm có chiều dài tiêu chuẩn 4 m, diện tích mặt cắt 2 x 10^-3 m^2, mô đun đàn hồi Young 2 x 10^11 N/m^2 và khối lượng riêng 7850 kg/m^3, được chia rời rạc từ 5 đến 50 phần tử hữu hạn vi cấu trúc. Phương pháp chọn mẫu tập trung vào việc bao quát toàn diện các điều kiện biên thực tế trong xây dựng và cơ khí chính xác, bao gồm: dầm hai đầu khớp, dầm hai đầu ngàm, dầm công-xôn một đầu ngàm một đầu tự do, dầm một đầu ngàm một đầu khớp, cùng các dạng dầm có tiết diện thay đổi bậc và thay đổi đều.

Luận văn lựa chọn phương pháp phân tích giải tích số kết hợp nguyên lý Hamilton vì phương pháp này cho phép thiết lập trực tiếp ma trận độ cứng tổng thể K từ các biểu thức thế năng đàn hồi, thế năng ngoại lực và động năng. Thuật toán giải phương trình đặc trưng định thức Det(K) = 0 được lập trình hoàn chỉnh trên môi trường Matlab R2011a trong khoảng thời gian nghiên cứu 6 tháng. Lý do lựa chọn phương pháp này là khả năng kiểm soát số bậc tự do linh hoạt, dễ dàng trích xuất đồng thời tần số dao động riêng và tải trọng tới hạn, đồng thời thuận tiện so sánh kiểm chứng trực tiếp với phần mềm phần tử hữu hạn thương mại SAP2000.

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

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

Thứ nhất, kết quả khảo sát sự hội tụ cho thấy khi tăng số lượng phần tử rời rạc từ 5 lên 50 phần tử, lực tới hạn và tần số dao động riêng tăng dần và đạt trạng thái ổn định tiệm cận hoàn hảo tại mốc n = 50 phần tử. Mức độ thay đổi giá trị khi tăng từ 45 lên 50 phần tử là không đáng kể, khẳng định mật độ 50 phần tử là ngưỡng chia tối ưu về mặt kinh tế tính toán.

Thứ hai, việc mô hình hóa độ cứng lò xo tại liên kết ngàm C1 đóng vai trò quyết định đến độ chính xác của bài toán dầm hai đầu ngàm. Khi thiết lập C1 = C, lực tới hạn đạt 7.579.490 N với sai số âm 3,87% so với lời giải chính xác Euler. Khi tăng độ cứng lên C1 = 2C, sai số tiệm cận về mức tối thiểu, trùng khớp hoàn toàn với lý thuyết giải tích. Khi độ cứng biên tăng lên mức rất lớn C1 = 100C, sai số lực tới hạn chỉ dừng ở mức 4,93%, chứng minh tính ổn định cao của mô hình.

Thứ ba, trong phân tích phi tuyến từng bước với 25 cấp tăng tải, số gia mỗi cấp đạt 315.795 N đối với lực dọc trục và 31.579,5 N/m đối với tải phân bố đều, mô hình vi kết cấu với C1 = 2C cho đường cong chuyển vị tại điểm giữa dầm trùng khớp tuyệt đối với kết quả phân tích từ phần mềm SAP2000 qua 23 cấp tải đầu tiên trước khi đạt điểm mất ổn định tới hạn.

Thứ tư, đối với dầm có tiết diện biến thiên với diện tích từ 6 x 10^-3 m^2 giảm xuống 2 x 10^-3 m^2 và mô men quán tính 5,6 x 10^-5 m^4, mô hình vi kết cấu phản ánh chính xác dạng dao động riêng mode 1 và mode 2, đồng thời xác định chính xác vị trí tập trung ứng suất uốn lớn nhất.

Thảo luận kết quả

Nguyên nhân tạo nên sự tương đồng chặt chẽ giữa mô hình dầm vi kết cấu và mô hình dầm phi cục bộ liên tục là do bản chất rời rạc của hệ lò xo xoay đã mô phỏng trực tiếp tương tác nội tại giữa các hạt nguyên tử liền kề. So sánh với các nghiên cứu trước đây của Wang năm 2013 và Challamel năm 2013, mô hình vi kết cấu trong luận văn đã mở rộng thành công phạm vi áp dụng từ dầm khớp đơn giản sang các điều kiện biên ngàm phức tạp và dầm tiết diện thay đổi.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các đồ thị quan hệ lực - chuyển vị, biểu đồ tần số dao động riêng dạng sóng và các bảng so sánh sai số định lượng. Đường cong chuyển vị uốn giữa dầm trên biểu đồ thể hiện rõ xu hướng suy giảm biên độ khi tăng độ cứng lò xo ngàm từ C lên 4C. Phát hiện này mang ý nghĩa to lớn đối với kỹ thuật tính toán, cung cấp cơ sở để các kỹ sư thiết kế cảm biến vi cơ MEMS có thể dự đoán chính xác hiện tượng dính bám mất ổn định mà không cần dựa vào các phần mềm mô phỏng phân tử đắt đỏ.

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

Thứ nhất, chuẩn hóa giá trị độ cứng lò xo liên kết ngàm C1 = 2C vào các tiêu chuẩn tính toán thiết kế kết cấu vi mô. Các viện nghiên cứu và phòng phát triển công nghệ MEMS/NEMS cần áp dụng ngay công thức này trong giai đoạn 2025 - 2026 để giảm thiểu sai số mô phỏng kết cấu ngàm xuống dưới 1%.

Thứ hai, phát triển mô-đun phần mềm tính toán chuyên dụng tích hợp thuật toán dầm vi kết cấu trên nền tảng ngôn ngữ lập trình mã nguồn mở như Python hoặc Matlab. Các nhóm nghiên cứu cơ học tính toán cần hoàn thành bộ công cụ này trong vòng 12 tháng, hướng tới mục tiêu rút ngắn 40% thời gian xử lý dữ liệu ma trận so với phương pháp phần tử hữu hạn thông thường.

Thứ ba, mở rộng phạm vi nghiên cứu từ dầm một chiều sang mô hình tấm và vỏ vi kết cấu hai chiều chịu tải trọng động học phức tạp. Bộ môn Công trình và Cơ học ứng dụng tại các trường đại học kỹ thuật nên triển khai các đề tài nghiên cứu chuyên sâu trong giai đoạn 2026 - 2027 nhằm đáp ứng yêu cầu tính toán của vật liệu cấu trúc nano hai chiều graphene.

Thứ tư, phối hợp với các phòng thí nghiệm vật liệu nano để tiến hành các thử nghiệm uốn tĩnh và dao động cộng hưởng trên mẫu dầm vi mô thực tế có kích thước dưới 100 nm. Hoạt động kiểm chứng thực nghiệm cần được thực hiện trong thời hạn 18 tháng nhằm hiệu chỉnh chính xác hệ số chiều dài đặc trưng a theo từng loại vật liệu chế tạo cụ thể.

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

Thứ nhất, kỹ sư thiết kế hệ vi cơ điện tử MEMS và NEMS: Tài liệu cung cấp cơ sở giải tích chuẩn xác để tối ưu hóa kích thước dầm cộng hưởng, màng cảm ứng vi mô và vi công tắc với độ dung sai gia công dưới 5 nm.

Thứ hai, các nhà nghiên cứu trong lĩnh vực cơ học tính toán và vật liệu tiên tiến: Luận văn là tài liệu tham khảo giá trị về cách ứng dụng nguyên lý biến phân Hamilton, biến đổi vi phân giả và khai triển Pade để giải quyết các hệ phương trình vi phân phi cục bộ bậc cao.

Thứ ba, học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng, Cơ kỹ thuật: Đây là biểu mẫu chuẩn mực về phương pháp xây dựng thuật toán Matlab, kỹ thuật phân tích hội tụ số và phương pháp kiểm chứng kết quả mô phỏng với phần mềm SAP2000.

Thứ tư, giảng viên đại học phụ trách các môn học Động lực học công trình và Cơ học kết cấu nâng cao: Công trình cung cấp hệ thống bài tập lớn, ví dụ số chi tiết và giáo trình chuyên đề về mối liên hệ giữa cơ học rời rạc và cơ học môi trường liên tục.

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

Tại sao lý thuyết đàn hồi cổ điển không thể áp dụng chính xác cho kết cấu nano từ 0,1 nm đến 100 nm?

Lý thuyết cổ điển xem vật chất là môi trường liên tục đồng nhất và giả định ứng suất tại một điểm chỉ phụ thuộc duy nhất vào biến dạng tại điểm đó. Khi kích thước cấu kiện giảm xuống dưới 100 nm, lực tương tác giữa các nguyên tử kế cận trở nên vượt trội, khiến lý thuyết cổ điển đánh giá sai độ cứng và tần số dao động của kết cấu vi mô.

Vai trò chính của hệ số ảnh hưởng quy mô nhỏ trong tính toán dầm vi kết cấu là gì?

Hệ số e0 đóng vai trò là tham số hiệu chỉnh mối quan hệ ứng suất - biến dạng phi cục bộ. Giá trị e0 dao động từ 0,289 khi phân tích uốn dọc đến 0,408 khi phân tích dao động, giúp mô hình hóa chính xác sự suy giảm độ cứng tổng thể của cấu kiện dầm vi mô dưới tác dụng của tải trọng nén dọc trục ban đầu.

Vì sao số lượng 50 phần tử được lựa chọn làm ngưỡng hội tụ tối ưu trong nghiên cứu?

Khảo sát số trên dầm 4 m cho thấy kết quả lực tới hạn và tần số dao động riêng biến thiên rõ rệt từ 5 đến 30 phần tử, nhưng bắt đầu đi vào ổn định từ 40 phần tử và tiệm cận không đổi tại 50 phần tử. Việc chia nhiều hơn 50 phần tử chỉ làm tăng khối lượng tính toán ma trận mà không cải thiện độ chính xác.

Giá trị độ cứng lò xo ngàm C1 = 2C mang lại ưu điểm vượt trội nào so với các mức độ cứng khác?

Độ cứng C1 = 2C thiết lập sự cân bằng cơ học hoàn hảo tại biên ngàm, triệt tiêu sai số góc xoay cục bộ. Kết quả kiểm chứng qua 23 bước tăng tải trên SAP2000 chứng minh C1 = 2C cho đường cong chuyển vị và tải trọng tới hạn trùng khớp hoàn toàn với phương pháp phần tử hữu hạn, trong khi C1 = C gây sai số 3,87%.

Mô hình dầm vi kết cấu có khả năng ứng dụng cho dầm tiết diện thay đổi không?

Mô hình hoàn toàn thích ứng tốt với dầm có tiết diện thay đổi bậc hoặc thay đổi đều. Thông qua việc gán ma trận độ cứng lò xo và khối lượng riêng biệt cho từng đoạn phần tử, thuật toán đã giải quyết chính xác bài toán dao động của dầm tiết diện biến thiên từ diện tích 6 x 10^-3 m^2 về 2 x 10^-3 m^2.

Kết luận

  • Luận văn đã xây dựng thành công mô hình dầm vi kết cấu hoàn chỉnh dựa trên nguyên lý Hamilton và lý thuyết phi cục bộ Eringen để phân tích toàn diện dao động riêng và mất ổn định đàn hồi.
  • Thiết lập quy tắc xác định độ cứng lò xo liên kết ngàm tối ưu C1 = 2C, giải quyết triệt để bài toán mô phỏng điều kiện biên ngàm trong cơ học vi cấu trúc với sai số dưới 1%.
  • Chứng minh tính hội tụ bền vững của thuật toán với ngưỡng chia 50 phần tử rời rạc, đồng thời kiểm chứng độ tin cậy tuyệt đối qua 12 bài toán so sánh với SAP2000 và nghiệm giải tích.
  • Đóng góp giải pháp tính toán hiệu quả cho dầm tiết diện thay đổi, mở ra phương pháp xác định hệ số quy mô nhỏ e0 mà không cần phụ thuộc vào các thí nghiệm cơ học phức tạp.
  • Lộ trình phát triển tiếp theo trong 12 - 24 tháng tới là mở rộng thuật toán cho kết cấu tấm nano hai chiều và tích hợp thư viện tính toán vi cơ học vào các phần mềm CAE thương mại. Bạn hãy tải toàn văn luận văn để tham khảo chi tiết mã nguồn Matlab và ứng dụng ngay vào các dự án nghiên cứu kết cấu vi mô của mình.