Trả lời trực tiếp: hội tụ lưới (mesh convergence) là quy trình tăng dần mật độ phần tử và so sánh kết quả để xác nhận giá trị đã ổn định — nhưng có một trường hợp không bao giờ hội tụ: tại các điểm kỳ dị toán học (stress singularity, ví dụ góc vuông trong hoặc điểm tải/ràng buộc tập trung), ứng suất tính được sẽ tăng liên tục không giới hạn khi lưới càng mịn. Đây không phải lỗi phần mềm mà là bản chất toán học của mô hình lý tưởng hóa — không nhận ra điều này là sai lầm phổ biến nhất khi đọc kết quả FEA. cokhi.net đã có bài mô phỏng ứng suất bằng Ansys giải thích quy trình thực hành chung — bài này đào sâu riêng chủ đề lưới và hội tụ.
Nội dung thuộc chuyên trang CAD – CAM – CAE của cokhi.net, cụm Mô phỏng CAE.
Loại phần tử và bậc phần tử
| Đặc điểm | Tứ diện (tetrahedral) | Lục diện (hexahedral) |
|---|---|---|
| Tạo lưới tự động | Dễ, phù hợp mọi hình dạng phức tạp | Khó tự động hóa, thường cần hình học đơn giản/chia khối thủ công |
| Độ chính xác cùng số nút | Thấp hơn, cần nhiều phần tử hơn để đạt cùng độ chính xác | Cao hơn — là lựa chọn ưu tiên khi hình học cho phép |
| Mức dùng phổ biến | Mặc định trong đa số phần mềm cho mô hình phức tạp | Dùng khi cần độ chính xác cao trên hình dạng đơn giản (dầm, tấm, trục) |
Bên cạnh loại phần tử, bậc phần tử cũng ảnh hưởng độ chính xác: phần tử bậc 1 (linear) có cạnh thẳng, xấp xỉ chuyển vị tuyến tính trong phần tử; phần tử bậc 2 (quadratic) có thêm nút giữa cạnh, cạnh cong theo hình dạng thật, xấp xỉ chính xác hơn nhiều cho vùng uốn/cong nhưng tốn tài nguyên tính toán hơn đáng kể (nhiều bậc tự do hơn trên cùng một lưới).
Quy trình chạy hội tụ lưới đúng cách
- Chạy mô phỏng ở mức lưới thô ban đầu, ghi lại giá trị quan tâm (thường là ứng suất đỉnh tại vị trí nguy hiểm nhất).
- Làm mịn lưới toàn bộ hoặc cục bộ tại vùng quan tâm, chạy lại, ghi giá trị mới.
- Lặp lại 3-4 mức mật độ lưới tăng dần, vẽ đồ thị giá trị quan tâm theo số phần tử (hoặc số nút).
- Khi giá trị thay đổi nhỏ (thường dưới 5%) giữa 2 mức liên tiếp, coi là đã hội tụ — dùng giá trị ở mức đó, không cần làm mịn thêm (tốn tài nguyên tính toán vô ích).
Trường hợp KHÔNG BAO GIỜ hội tụ: điểm kỳ dị ứng suất
Tại các vị trí sau, lý thuyết đàn hồi cổ điển cho ứng suất tiến tới vô hạn về mặt toán học:
- Góc vuông trong sắc nét (re-entrant corner): góc 90° lõm vào trong hình học, không có bo góc.
- Điểm tải hoặc điểm ràng buộc tập trung tại 1 điểm/1 cạnh (thay vì phân bố trên một diện tích).
Khi mô phỏng các vị trí này, càng làm mịn lưới, giá trị ứng suất càng tăng — không bao giờ dừng lại ở một số ổn định. Đây KHÔNG phải "chưa hội tụ cần làm mịn thêm" mà là bản chất mô hình không cho hội tụ. Sai lầm phổ biến nhất của người mới dùng FEA là cố làm mịn lưới liên tục hòng "tìm ra giá trị đúng" tại các điểm này, trong khi giá trị đúng về mặt vật lý thực tế không tồn tại ở dạng đó (vì chi tiết thật luôn có bo góc nhỏ, tải luôn phân bố trên diện tích chứ không phải 1 điểm toán học).
Cách xử lý đúng:
- Nếu chi tiết thật có bo góc (hầu hết trường hợp thực tế), mô hình hóa đúng bán kính bo đó thay vì để góc vuông sắc trong CAD — khi đó ứng suất sẽ hội tụ về giá trị hữu hạn thật.
- Nếu không thể tránh singularity (một số điều kiện biên lý tưởng hóa), chỉ dùng giá trị tại đó để so sánh xu hướng tương đối giữa các phương án thiết kế, không dùng làm giá trị tuyệt đối để so sánh với giới hạn chảy vật liệu.
FAQ
Vì sao 2 lần chạy mesh khác nhau ra kết quả khác nhau đáng kể? Có 2 khả năng: (1) chưa hội tụ — cần làm mịn thêm và kiểm tra lại theo quy trình ở trên; (2) đang ở vùng singularity — giá trị sẽ tiếp tục tăng mãi dù làm mịn bao nhiêu, cần kiểm tra hình học tại vị trí đó (góc vuông sắc? tải tập trung 1 điểm?).
Phần tử bậc 2 luôn tốt hơn bậc 1? Chính xác hơn cho cùng số phần tử, đặc biệt ở vùng uốn/cong, nhưng tốn tài nguyên tính toán hơn nhiều (nhiều bậc tự do hơn). Với hình dạng đơn giản và lưới đủ mịn, phần tử bậc 1 vẫn có thể cho kết quả đủ tin cậy với chi phí tính toán thấp hơn.
Làm sao phân biệt "chưa hội tụ" với "đang gặp singularity"? Chạy 3-4 mức mịn dần và nhìn xu hướng: nếu giá trị tăng chậm dần rồi chững lại (hội tụ) — dùng được; nếu giá trị tăng liên tục không có dấu hiệu chững, đặc biệt tại góc vuông sắc hoặc điểm tải/ràng buộc tập trung — khả năng cao là singularity, cần xử lý theo hướng ở mục trên thay vì tiếp tục làm mịn.
Thuật ngữ
- h-refinement: làm mịn lưới bằng cách tăng số lượng phần tử (giảm kích thước từng phần tử).
- p-refinement: làm mịn lưới bằng cách tăng bậc đa thức xấp xỉ của phần tử (ví dụ chuyển bậc 1 sang bậc 2) mà không tăng số phần tử.
- Singularity (điểm kỳ dị): vị trí lời giải lý thuyết cho giá trị tiến tới vô hạn, không hội tụ khi làm mịn lưới.
- Re-entrant corner (góc lõm sắc nét): góc trong của hình học không có bo, nguồn gốc phổ biến nhất của singularity trong phân tích ứng suất.
Tài liệu tham khảo
- NAFEMS, A Designer's Guide to Simulation Governance and Best Practice (quy trình hội tụ lưới thực hành): https://www.nafems.org/
- R.D. Cook, D.S. Malkus, M.E. Plesha, R.J. Witt, Concepts and Applications of Finite Element Analysis — giáo trình chuẩn ngành về phần tử hữu hạn, chương hội tụ và kỳ dị ứng suất.
- Ansys Help, Understanding and Avoiding Singularities in Structural Analysis: https://www.ansys.com/