Luận án tiến sĩ: Nghiên cứu đặc tính khí động đàn hồi cánh máy bay

Tam giác Collar biểu diễn mối liên hệ giữa ba lực chính tác động lên cánh. Lực khí động phát sinh từ chuyển động tương đối giữa cánh và không khí. Lực đàn hồi xuất hiện khi kết cấu cánh bị biến dạng. Lực quán tính liên quan đến gia tốc của khối lượng cánh. Ba lực này tương tác liên tục, tạo nên hiện tượng khí động đàn hồi. Tương tác khí động đàn hồi mạnh nhất ở vận tốc cao và kết cấu mềm. Mô phỏng CFD giúp quan sát chi tiết dòng khí quanh cánh. Phương pháp phần tử hữu hạn phân tích ứng xử kết cấu dưới tải trọng.

Đàn hồi kết cấu là khả năng cánh biến dạng và phục hồi dưới tác dụng của lực. Vật liệu composite hiện đại cho phép thiết kế cánh nhẹ nhưng bền. Độ cứng chống xoắn GJ và độ cứng chống uốn là thông số quan trọng. Biến dạng khí động gồm biến dạng uốn và biến dạng xoắn. Biến dạng uốn làm cánh võng xuống hoặc võng lên. Biến dạng xoắn thay đổi góc tấn dọc theo sải cánh. Cơ học kết cấu nghiên cứu phản ứng của cánh với tải trọng phức tạp.

Cánh có biên dạng đặc thù như cánh siêu tới hạn mang lại ưu điểm khí động vượt trội. Thiết kế đặc biệt giảm lực cản, tăng hiệu suất bay. Tuy nhiên, biên dạng phức tạp tạo ra thách thức về khí động đàn hồi. Phân tích chuyên sâu cần thiết để tránh hiện tượng nguy hiểm. Flutter cánh máy bay có thể xảy ra ở tốc độ nhất định. Xoắn phá hủy và đảo chiều tác dụng cánh lái là rủi ro cần kiểm soát. Nghiên cứu này góp phần nâng cao độ tin cậy thiết kế.

Xoắn phá hủy xảy ra khi vận tốc bay vượt quá vận tốc tới hạn Vdiv. Áp suất động q tăng theo bình phương vận tốc. Lực nâng tạo mô-men xoắn quanh trục dọc cánh. Góc xoắn θ tăng làm góc tấn tăng, sinh thêm lực nâng. Vòng phản hồi dương này dẫn đến mất ổn định. Độ cứng chống xoắn GJ càng lớn, Vdiv càng cao. Khoảng cách giữa tâm khí động và tâm cứng ec ảnh hưởng mạnh. Phân tích FSI xác định vận tốc an toàn khai thác.

Đảo chiều xảy ra khi cánh lái mất hiệu quả do biến dạng đàn hồi. Khi mở cánh lái góc β, tạo ra lực khí động bổ sung. Lực này gây xoắn cánh theo hướng giảm góc tấn hiệu dụng. Ở vận tốc thấp, hiệu ứng cánh lái chiếm ưu thế. Ở vận tốc cao, biến dạng xoắn triệt tiêu tác dụng cánh lái. Hệ số giảm tác dụng điều khiển ξ đo lường hiện tượng này. Thiết kế phải đảm bảo vận tốc đảo chiều cao hơn vận tốc khai thác. Tương tác khí động đàn hồi phức tạp hơn ở biên dạng đặc thù.

Flutter là dao động tự kích kết hợp giữa chuyển động uốn và xoắn. Năng lượng từ dòng khí cung cấp cho dao động. Khi tần số dao động uốn wh gần bằng tần số xoắn wα, nguy cơ flutter tăng. Vận tốc flutter phụ thuộc vào khối lượng, độ cứng và hình dạng cánh. Hệ số giảm chấn khí động γ quyết định ổn định hay không ổn định. Phương pháp Modal phân tích tần số và dạng dao động riêng. Động lực học chất lỏng phi tuyến làm phức tạp bài toán. Mô phỏng CFD phi định thường cần thiết cho dự báo chính xác.

Cánh siêu tới hạn giảm cường độ sóng xung kích bằng cách phân bố lại áp suất. Bề mặt trên phẳng hơn làm chậm gia tốc dòng khí. Vận tốc cục bộ tăng chậm hơn, trì hoãn đạt số Mach 1. Sóng xung kích khi xuất hiện yếu hơn và ở vị trí xa hơn về phía sau. Phần dưới cong nhiều hơn bù lại lực nâng. Hệ số lực nâng CL duy trì ở mức cao. Hệ số lực cản CD giảm 10-15% so với biên dạng thông thường. Khí động lực học phức tạp đòi hỏi phân tích chi tiết.

Biên dạng phức tạp tạo phân bố tải trọng không đồng đều. Biến dạng khí động ảnh hưởng mạnh đến hiệu suất cánh. Vùng sóng xung kích di chuyển khi cánh biến dạng. Tương tác khí động đàn hồi mạnh hơn ở vận tốc xuyên âm. Phương pháp tuyến tính không còn chính xác. Phân tích FSI phi tuyến cần thiết để dự báo đúng. Độ cứng kết cấu phải được tối ưu hóa cẩn thận. Quá cứng làm tăng trọng lượng, quá mềm gây mất ổn định.

Cánh siêu tới hạn được sử dụng rộng rãi trên máy bay thương mại hiện đại. Boeing 787 và Airbus A350 đều áp dụng công nghệ này. Tiết kiệm nhiên liệu 5-10% mang lại lợi ích kinh tế lớn. Tầm bay tăng hoặc tải trọng tăng với cùng nhiên liệu. Thiết kế đòi hỏi nghiên cứu khí động đàn hồi kỹ lưỡng. Thử nghiệm bay và mô phỏng CFD xác nhận an toàn. Cơ học kết cấu composite hỗ trợ thiết kế tối ưu. Công nghệ này đại diện cho tiến bộ quan trọng trong hàng không.

Mô phỏng CFD giải các phương trình Navier-Stokes mô tả chuyển động chất lưu. Phương pháp thể tích hữu hạn chia miền tính toán thành các ô lưới. Mô hình rối như k-epsilon hoặc k-omega SST mô phỏng dòng rối. Điều kiện biên xác định vận tốc vào, áp suất ra và tường không trượt. Số Reynolds Re và số Mach M là thông số quan trọng. Hệ số áp suất Cp cho biết phân bố áp suất trên bề mặt cánh. Lực nâng L và lực cản D tính từ tích phân áp suất và ứng suất cắt. Độ chính xác phụ thuộc vào chất lượng lưới và mô hình rối.

Phương pháp phần tử hữu hạn chia kết cấu thành các phần tử nhỏ. Mỗi phần tử có ma trận độ cứng [K] và ma trận khối lượng [M]. Tải trọng khí động áp dụng lên các nút phần tử. Phương trình chuyển động [M]{ẍ} + [K]{x} = {F} được giải. Vec-tơ chuyển vị {q} cho biết biến dạng tại mỗi nút. Vật liệu composite cần mô hình trực hướng với mô-đun đàn hồi khác nhau. Mô-đun đàn hồi Ex, Ey và mô-đun trượt G là thông số vật liệu. Phân tích modal xác định tần số riêng và vector riêng [Φ].

Phân tích FSI kết nối mô phỏng CFD và FEM qua giao diện. Lực khí động từ CFD truyền sang kết cấu làm biến dạng. Biến dạng kết cấu thay đổi hình dạng cánh, cập nhật lưới CFD. Quá trình lặp tiếp tục đến khi hội tụ. Phương pháp ghép mạnh giải đồng thời cả hai miền. Phương pháp ghép yếu giải luân phiên từng miền. Tương tác khí động đàn hồi được mô phỏng chính xác. Phương pháp này cần thiết cho biên dạng đặc thù và vận tốc cao. Thời gian tính toán lớn nhưng kết quả đáng tin cậy.

Bài toán CFD kiểm nghiệm với cánh NACA65A004 hình chữ nhật. Số Reynolds Re = 1.5×10⁶ tương ứng điều kiện bay thực tế. Hệ số lực nâng CL sai lệch dưới 3% so với thực nghiệm. Hệ số lực cản CD có độ chính xác tương tự. Phân bố áp suất Cp trên bề mặt cánh trùng khớp tốt. Phương pháp Modal xác định tần số riêng với sai số dưới 2%. Dạng dao động riêng phù hợp với lý thuyết. Phân tích FSI cho biến dạng chính xác trong phạm vi 4%.

Cánh siêu tới hạn cho hệ số lực nâng cao hơn cánh thông thường 8%. Hệ số lực cản giảm 12% ở số Mach 0.75. Sóng xung kích xuất hiện muộn hơn và yếu hơn. Phân bố áp suất trên bề mặt trên phẳng hơn. Vùng áp suất âm kéo dài nhưng cường độ thấp hơn. Hệ số mô-men chúc ngóc Cm ổn định hơn. Hiệu suất khí động lực học vượt trội rõ rệt. Kết quả mô phỏng CFD phù hợp với đo đạc thực nghiệm.

Biến dạng uốn tăng tuyến tính với vận tốc bay. Biến dạng xoắn phức tạp hơn do phân bố tải không đều. Cánh siêu tới hạn có biến dạng xoắn lớn hơn 15% so với cánh thông thường. Tương tác khí động đàn hồi làm thay đổi phân bố lực nâng. Hiệu ứng phản hồi dương xuất hiện ở vận tốc cao. Vận tốc xoắn phá hủy Vdiv xác định chính xác. Vận tốc flutter dự báo với độ tin cậy cao. Kết quả cho thấy cần thiết tối ưu độ cứng kết cấu.

Phương pháp nghiên cứu áp dụng trực tiếp cho thiết kế cánh máy bay mới. Dự báo chính xác vận tốc tới hạn đảm bảo an toàn bay. Tối ưu hóa độ cứng kết cấu giảm trọng lượng mà vẫn đảm bảo độ bền. Phân tích FSI giúp phát hiện sớm vấn đề tiềm ẩn. Thời gian phát triển rút ngắn nhờ mô phỏng chính xác. Chi phí thử nghiệm giảm khi số lần thử ít hơn. Kết quả nghiên cứu góp phần nâng cao năng lực thiết kế trong nước.

Thiết kế cánh cần cân bằng nhiều mục tiêu đối kháng. Tăng lực nâng thường đi kèm tăng lực cản. Giảm trọng lượng có thể làm giảm độ cứng. Tối ưu hóa đa mục tiêu tìm giải pháp tốt nhất. Thuật toán di truyền và tối ưu hóa bầy đàn hiệu quả. Không gian thiết kế rộng lớn với biên dạng đặc thù. Tích hợp phân tích khí động đàn hồi vào vòng lặp tối ưu. Kết quả là thiết kế vượt trội về nhiều khía cạnh.

Nghiên cứu mở rộng sang cánh thích ứng có hình dạng thay đổi. Vật liệu thông minh cho phép điều khiển biến dạng chủ động. Phân tích khí động đàn hồi phi tuyến hoàn toàn cần phát triển thêm. Mô phỏng đa tỷ lệ kết nối từ vi mô đến vĩ mô. Trí tuệ nhân tạo hỗ trợ dự báo nhanh và tối ưu hóa. Thực nghiệm với mô hình tỷ lệ lớn xác nhận kết quả. Hợp tác quốc tế mở ra cơ hội tiếp cận công nghệ tiên tiến. Nghiên cứu khí động đàn hồi tiếp tục là lĩnh vực quan trọng.