Mô phỏng phần tử hữu hạn (finite element method, FEM) là phương pháp số học tính toán xấp xỉ dùng để giải các phương trình vi phân đạo hàm riêng chi phối các bài toán cơ học vật rắn, truyền nhiệt, động lực học chất lưu và trường điện từ bằng cách rời rạc hóa miền liên tục phức tạp thành tập hợp các phần tử con hình học đơn giản.
Cơ sở toán học và dạng biến phân yếu
FEM chuyển đổi phương trình vi phân cấp cao (dạng mạnh) thành phương trình tích phân biến phân (dạng yếu) thông qua phương pháp phần dư có trọng số Galerkin hoặc nguyên lý cực tiểu thế năng toàn phần:
Trong đó là hàm nghiệm cần tìm, là hàm kiểm tra (test function), là ngoại lực nguồn và là miền không gian khảo sát.
Nghiệm trên từng phần tử được xấp xỉ qua tổ hợp tuyến tính của các hàm hình dạng (shape functions) :
Tích phân trên toàn bộ các phần tử rời rạc được ghép nối lại để thiết lập hệ phương trình đại số tuyến tính toàn cục:
Trong đó là ma trận độ cứng toàn cục, là véc-tơ chuyển vị tại các nút lưới, và là véc-tơ tải trọng nút tương đương.
| Bước thực hiện | Nội dung tính toán | Mục đích kỹ thuật |
|---|---|---|
| Tiền xử lý | Tạo hình học và chia lưới phần tử | Rời rạc hóa miền liên tục |
| Thiết lập hệ | Tính ma trận độ cứng và véc-tơ tải | Xây dựng hệ phương trình |
| Giải ma trận | Áp đặt điều kiện biên và giải số | Xác định chuyển vị nút |
| Hậu xử lý | Tính trường ứng suất và biến dạng | Đánh giá độ bền kết cấu |
Quy trình thực hiện mô phỏng FEM tiêu chuẩn
Một quy trình phân tích FEM hoàn chỉnh bao gồm ba giai đoạn chính:
- Tiền xử lý (Pre-processing): Xây dựng mô hình hình học CAD, lựa chọn loại phần tử thích hợp (phần tử thanh 1D, phần tử tấm/vỏ 2D, hoặc phần tử khối 3D tứ diện/lục diện), thiết lập đặc tính cơ lý của vật liệu (môđun đàn hồi Young, hệ số Poisson) và áp đặt các điều kiện biên ràng buộc cùng tải trọng tác dụng.
- Xử lý số (Solver Execution): Phần mềm tính toán ma trận phần tử, ghép nối ma trận toàn cục, khử các bậc tự do bị ràng buộc và sử dụng các thuật toán giải hệ ma trận lớn (như phân tích Cholesky, phương pháp gradient liên hợp) để tìm nghiệm chuyển vị tại từng nút lưới.
- Hậu xử lý (Post-processing): Trích xuất kết quả phân bố chuyển vị, tính toán tenxơ biến dạng và ứng suất tương đương von Mises, hiển thị biểu đồ đường mức ứng suất màu và kiểm tra hệ số an toàn chống phá hủy mỏi hoặc dẻo.
Các lớp bài toán kỹ thuật giải quyết bằng FEM
Phương pháp phần tử hữu hạn ngày nay đã trở thành công cụ không thể thay thế trong thiết kế kỹ thuật số hỗ trợ máy tính (CAE):
- Phân tích kết cấu tuyến tính và phi tuyến: Tính toán độ bền tĩnh, va chạm biến dạng lớn của khung xe ô tô, uốn dọc mất ổn định của tháp gió và mở rộng vết nứt mỏi trong cánh máy bay.
- Mô phỏng truyền nhiệt và ứng suất nhiệt: Khảo sát trường nhiệt độ dừng và chuyển tiếp trong động cơ tua-bin khí, lò luyện kim và tản nhiệt chip điện tử.
- Bài toán tương tác đa vật lý: Ghép đôi tương tác chất lưu - kết cấu (FSI), phân tích cơ - nhiệt - điện từ trong các hệ vi cơ điện tử MEMS.