Tối ưu hóa hình dạng là gì? Cơ sở toán học, Level-Set và ứng dụng CAD

Tối ưu hóa hình dạng là kỹ thuật toán học ứng dụng và cơ học tính toán tự động nhằm xác định hình dạng biên hình học tối ưu của một cấu trúc cơ khí, cánh khí động hoặc thiết bị nhiệt lỏng để đạt được hiệu năng kỹ thuật cao nhất. Mục tiêu tối ưu thường bao gồm cực tiểu hóa trọng lượng, giảm thiểu nồng độ ứng suất tập trung, triệt tiêu lực cản khí động học hoặc tối đa hóa tốc độ tản nhiệt dưới các ràng buộc cơ lý nghiêm ngặt. Khác với thiết kế truyền thống vốn dựa vào kinh nghiệm thử - sai của kỹ sư, tối ưu hóa hình dạng kết hợp các thuật toán đạo hàm giải tích và mô phỏng số phần tử hữu hạn để tự động biến đổi đường biên. Bài viết này trình bày nền tảng toán học, các phương pháp tính toán, quy trình phân tích độ nhạy và ứng dụng công nghiệp hiện đại.

Phân biệt các cấp độ tối ưu hóa cấu trúc

Theo công trình kinh điển của Bendsøe và Sigmund (2004), tối ưu hóa cấu trúc trong cơ học kỹ thuật được phân cấp thành ba mức độ với không gian biến thiết kế mở rộng dần:

Cấp độ tối ưu	Biến thiết kế chính	Đặc điểm hình học và bậc tự do
Tối ưu hóa kích thước (Sizing Optimization)	Độ dày bản mỏng, diện tích mặt cắt ngang dầm thanh	Hình dạng và tô-pô cấu trúc giữ nguyên cố định
Tối ưu hóa hình dạng (Shape Optimization)	Tọa độ các điểm nút trên đường biên hoặc bề mặt biên	Tô-pô (số lượng lỗ rỗng) giữ nguyên, biên tự do uốn lượn
Tối ưu hóa tô-pô (Topology Optimization)	Phân bố mật độ vật liệu trong toàn bộ miền không gian thiết kế	Tự do tạo ra các lỗ rỗng mới và cấu trúc khung liên kết

Cơ sở toán học và bài toán biến phân biên

Xét một miền vật lý $\Omega(\Gamma)$ được giới hạn bởi đường biên khả vi $\Gamma$. Bài toán tối ưu hóa hình dạng tổng quát được phát biểu dưới dạng cực tiểu hóa phiếm hàm mục tiêu $J$ phụ thuộc vào trường chuyển vị $u(\Gamma)$ và hình dạng biên $\Gamma$:

Tìm đường biên tối ưu $\Gamma^*$ sao cho:

Cực tiểu hóa phiếm hàm mục tiêu $J$ phụ thuộc trường trạng thái $u$ thỏa mãn phương trình trạng thái cơ lý $A(\Gamma)$ với tải trọng $p$ cùng các ràng buộc hình học về biên $\Gamma$.

Theo lý thuyết toán học của Sokolowski và Zolesio (1992), để tìm hướng biến đổi biên làm giảm nhanh nhất hàm mục tiêu $J$, người ta sử dụng phương pháp đạo hàm hình dạng (shape derivative). Khi biên $\Gamma$ chịu một nhiễu loạn vận tốc pháp tuyến $V_n(x)$, đạo hàm Euler của phiếm hàm $J$ được tính theo tích phân trên biên:

Biểu thức đạo hàm hình dạng có dạng tích phân trên biên của tích giữa hàm gradient nhạy $g(x)$ và vận tốc dịch chuyển pháp tuyến $V_n(x)$. Để giảm nhanh nhất giá trị của hàm mục tiêu $J$, hướng dịch chuyển biên tối ưu được chọn ngược hướng với gradient độ nhạy: $V_n = -g(x)$ kết hợp với hệ số bước dịch chuyển lpha.

Các phương pháp tính toán tối ưu hóa hình dạng

1. Phương pháp tham số hóa hình học (CAD-based Parametric Optimization)

Đường biên $\Gamma$ được mô tả thông qua các đường cong tham số mượt như B-spline hoặc NURBS trong môi trường thiết kế CAD. Các biến tối ưu là tọa độ của các điểm điều khiển (control points). Phương pháp này đảm bảo mô hình đầu ra có thể chuyển thẳng vào dây chuyền gia công CNC mà không cần xử lý lại hình học, tuy nhiên số lượng bậc tự do bị giới hạn bởi số điểm điều khiển.

2. Phương pháp phân tích độ nhạy biến dạng nút (Node-based Mesh Morphing)

Trực tiếp coi tọa độ các điểm nút biên trên lưới phần tử hữu hạn là biến thiết kế. Phương pháp này có độ linh hoạt rất cao nhưng dễ dẫn đến hiện tượng méo mó lưới phần tử hữu hạn sau nhiều bước lặp, đòi hỏi thuật toán làm mịn biên và tái chia lưới (remeshing) liên tục.

3. Phương pháp tập mức (Level-Set Method)

Công trình đột phá của Allaire, Jouve và Toader (2004) đã áp dụng phương pháp hàm tập mức (Level-Set) vào tối ưu hóa hình dạng cấu trúc. Biên $\Gamma$ được biểu diễn ẩn dưới dạng đường đẳng mức không của hàm vô hướng nhiều chiều $\phi(x)$. Sự dịch chuyển của biên theo vận tốc pháp tuyến $V_n(x)$ được giải thông qua phương trình đạo hàm riêng Hamilton-Jacobi trên lưới Euler cố định. Phương pháp này giải quyết triệt để vấn đề méo lưới và cho phép kết hợp tối ưu hình dạng với sự thay đổi tô-pô một cách tự nhiên.

Quy trình kỹ thuật mô phỏng và tối ưu hóa

Một chu trình tối ưu hóa hình dạng tự động trong công nghiệp bao gồm các bước tuần tự lặp kín:

• Khởi tạo mô hình CAD/CAE: Thiết lập hình dạng ban đầu của chi tiết và xác định miền biên $\Gamma$ được phép thay đổi.
• Phân tích trạng thái cơ lý (FEA/CFD): Giải hệ phương trình trạng thái (độ bền kết cấu, truyền nhiệt hoặc dòng chảy khí động) bằng phần tử hữu hạn để thu được trường nghiệm $u$.
• Phân tích độ nhạy liên hợp (Adjoint Sensitivity Analysis): Giải bài toán trạng thái liên hợp để tính toán gradient độ nhạy $g(x)$ trên toàn bộ biên mà không tốn kém chi phí tính toán sai phân hữu hạn.
• Cập nhật hình học biên: Dịch chuyển các nút biên theo hướng giảm hàm mục tiêu $J$ với bước lpha và kiểm tra điều kiện hội tụ.

Ứng dụng thực tiễn trong công nghiệp hiện đại

Tối ưu hóa hình dạng là công cụ then chốt trong nhiều lĩnh vực kỹ thuật công nghệ cao:

• Công nghiệp hàng không và ô tô: Tối ưu hóa biên dạng cánh máy bay (airfoil), thân xe hơi và hốc hút gió nhằm giảm thiểu lực cản và tối đa hóa lực nâng khí động học.
• Kỹ thuật cơ khí chính xác: Bo tròn các góc chuyển tiếp (fillet optimization) để triệt tiêu ứng suất tập trung, kéo dài tuổi thọ mỏi của trục khuỷu, tuabin khí và cầu nối chịu lực.
• Thiết bị y sinh: Tối ưu hóa hình dáng stent tim mạch và khớp nhân tạo để phân bố áp lực đồng đều lên mô sinh học xung quanh.

Nghiên cứu và đào tạo kỹ thuật tại Việt Nam

Tại Việt Nam, các phương pháp tối ưu hóa hình dạng và cơ học tính toán đang được giảng dạy và nghiên cứu chuyên sâu tại các trường đại học kỹ thuật hàng đầu như Trường Đại học Bách khoa Hà Nội, Đại học Bách khoa TP.HCM và các viện nghiên cứu chuyên ngành cơ học. Sự kết hợp giữa tối ưu hóa hình dạng số và công nghệ in 3D kim loại (Additive Manufacturing) đang mở ra hướng đột phá mới cho ngành chế tạo máy và công nghiệp quốc phòng trong nước.