Tổng quan luận án

Phân tích đáp ứng của thiết diện cánh máy bay là bài toán cơ học kỹ thuật quan trọng phục vụ quá trình thiết kế, chế tạo, vận hành và bảo dưỡng các phương tiện bay trong lĩnh vực dân sự và quân sự. Khi cánh máy bay chuyển động trong dòng khí, sự tương tác giữa lực khí động, lực quán tính và lực đàn hồi làm phát sinh các hiện tượng khí đàn hồi động lực học phức tạp, tiêu biểu là dao động tự kích flutter và dao động vòng giới hạn (Limit Cycle Oscillation - LCO). Các hệ phương trình vi phân mô tả chuyển động của thiết diện cánh thường mang tính phi tuyến mạnh. Phương pháp tuyến tính hóa tương đương kinh điển khi áp dụng cho các hệ phi tuyến này thường gặp hạn chế về độ chính xác của nghiệm. Trong khi đó, cách tiếp cận đối ngẫu trong cơ học phi tuyến do GS. Nguyễn Đông Anh khởi xướng ban đầu chỉ mới được nghiên cứu cho dao động ngẫu nhiên và hệ điều khiển giảm chấn khối lượng TMD, chưa từng được phát triển và áp dụng cho bài toán phân tích đáp ứng phi tuyến và mất ổn định flutter của thiết diện cánh máy bay.

Luận án xác định 4 mục tiêu nghiên cứu cụ thể được đánh số như sau:

  1. Phát triển phương pháp luận cho cách tiếp cận đối ngẫu trong phương pháp tuyến tính hóa tương đương áp dụng trong bài toán phân tích đáp ứng của thiết diện cánh chịu lực khí động.
  2. Xây dựng các cải tiến có hiệu quả cho tiêu chuẩn tuyến tính hóa tương đương đối ngẫu, áp dụng cho bài toán flutter phi tuyến của thiết diện cánh.
  3. Tăng độ chính xác cho nghiệm của bài toán ổn định flutter thiết diện cánh bằng cách áp dụng tiêu chuẩn đối ngẫu được cải tiến.
  4. Thu được các kỹ thuật tính toán theo tiêu chuẩn đối ngẫu cho bài toán flutter của thiết diện cánh.

Đối tượng nghiên cứu của luận án là thiết diện cánh máy bay theo mô hình hai chiều chịu tác động của lực khí động. Phạm vi nghiên cứu tập trung vào mô hình mặt cắt điển hình 2 bậc tự do (chuyển vị uốn $h$ và góc xoắn $\alpha$) chuyển động trong dòng khí không nén được; lực khí động được tính toán dựa trên lý thuyết cánh mỏng dừng và tựa dừng quy về lực tập trung tại tâm khí động; khảo sát tính phi tuyến kết cấu tập trung dạng đa thức bậc ba và khe hở tự do (free-play); phân tích hiện tượng đáp ứng dao động vòng giới hạn (LCO) và mở rộng khảo sát hệ thiết diện cánh có cánh nhỏ (aileron) điều khiển bằng thuật toán PID.

Tổng quan tài liệu và vị trí của luận án

Luận án điểm lại hệ thống các công trình nghiên cứu quốc tế và trong nước về lĩnh vực khí đàn hồi và phương pháp tuyến tính hóa tương đương:

Nghiên cứu quốc tế:

  • Nền tảng khí đàn hồi được Collar (1978) hệ thống hóa qua tam giác tương tác ba lực (lực khí động, lực quán tính, lực đàn hồi). Lịch sử và ứng dụng thiết kế máy bay được khảo sát bởi Garrick và Reid (1981), Flomenhoft (1997), Friedmann (1999), Livne (2003), Hodges và Pierce (2002), Dowell và cộng sự (2003, 2015).
  • Mô hình thiết diện cánh 2 chiều được chuẩn hóa bởi Fung (1993), Wright và cộng sự (2007), Strganac và cộng sự (2000), kết hợp mô hình FSI giữa tính toán động lực học chất lưu CFD và mô hình kết cấu tính toán CSM của Henshaw và cộng sự (2007).
  • Về phân tích phi tuyến và dao động vòng giới hạn (LCO): Yang và Zhao (1988) nghiên cứu lý thuyết kết hợp thực nghiệm mô hình cánh 2 chiều; Liu và Zhao (1992) sử dụng phương pháp cân bằng điều hòa và lý thuyết đa tạp trung tâm; Yang (1995) dùng phương pháp tiệm cận Krylov-Bogoliubov-Mitropolskii (KBM); Shahrzad và Mhazoon (2002) dự báo rẽ nhánh Hopf; Ding và Wang (2006) khảo sát flutter ở vận tốc trên âm; Chen Feixin và cộng sự (2012, 2013) phân tích flutter bằng tuyến tính hóa tương đương; Lee và cộng sự (1997, 1999, 2005) khảo sát giải pháp tính toán số cho thiết diện cánh; Zhao và Hu (2004) nghiên cứu kết hợp phi tuyến bậc 3 và khe hở tự do; Liu và cộng sự (2000, 2002), Chung và cộng sự (2007), Tang và cộng sự (2000, 2004a, 2004b) thực hiện mô phỏng số Runge-Kutta và thí nghiệm hầm gió; Djayapertapa và cộng sự (2001a, 2001b), Allen và cộng sự (2005) nghiên cứu phi tuyến khí động bằng CFD; Pines (1958) chứng minh điều kiện không xảy ra flutter khi trọng tâm nằm trước tâm đàn hồi; Wei và cộng sự (2014) thiết lập tiêu chuẩn ổn định LCO đối với nhiễu biên độ.

Nghiên cứu trong nước:

  • Trong lĩnh vực khí động và flutter: Lã Hải Dũng (2005) khảo sát flutter uốn - xoắn cánh; Hoàng Thị Bích Ngọc và cộng sự (2009) nghiên cứu hiện tượng đàn hồi cánh chịu lực khí động; Phạm Duy Hòa và cộng sự (2014) khảo sát ổn định flutter mặt cắt cầu treo có khe slot; Trần Thế Văn (2012) nghiên cứu tấm composite chịu tải khí động bằng phương pháp phần tử hữu hạn; Trần Ngọc An (2014) phát triển phương pháp bước lặp Matsumoto tính flutter dầm chủ cầu treo có thiết bị giảm chấn.
  • Về cách tiếp cận đối ngẫu: GS. Nguyễn Đông Anh (2010), Anh và cộng sự (2012a, 2012b, 2013) đề xuất kỹ thuật cực tiểu bình phương đối ngẫu cho hệ dao động ngẫu nhiên kích động ồn trắng, hệ nhiều bậc tự do và hệ TMD; TS. Nguyễn Ngọc Linh (2015) phát triển tiêu chuẩn đối ngẫu có trọng số cho hệ 1 bậc tự do chịu kích động ngẫu nhiên dừng Gauss.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết: Luận án của Nguyễn Ngọc Linh (2015) giới hạn trong dao động ngẫu nhiên phi tuyến, trong khi luận án này phát triển và mở rộng cách tiếp cận đối ngẫu sang lớp dao động tuần hoàn và hiện tượng mất ổn định flutter phi tuyến của thiết diện cánh máy bay, đồng thời đề xuất 3 kỹ thuật xác định tham số trọng số mới nhằm nâng cao độ chính xác so với phương pháp tuyến tính hóa kinh điển.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên các khung lý thuyết, khái niệm và công cụ phân tích sau:

  • Lý thuyết khí đàn hồi động học: Dựa trên tam giác Collar, mô tả sự tương tác qua lại giữa ba trường lực: lực đàn hồi (hệ số độ cứng lò xo $k_h, k_\alpha$), lực quán tính (khối lượng $m$, mô men tĩnh $S$, mô men quán tính $I$) và lực khí động.
  • Lý thuyết cánh mỏng dừng và tựa dừng (Thin Airfoil Theory): Sử dụng định lý Joukowsky, điều kiện biên vận tốc triệt tiêu $v/U = 0$, điều kiện thoát dòng Kutta $\gamma(c) = 0$, phép đổi biến Glauert $x = \frac{c}{2}(1 - \cos\theta)$ và phân bố xoáy vận tốc $\gamma(\theta) = 2U [A_0 \cot(\theta/2) + \sum A_n \sin(n\theta)]$ để thiết lập công thức tính lực nâng $L$ và mô men uốn $M$. Xác định vị trí tâm khí động nằm tại điểm $1/4$ dây cung tính từ biên dẫn.
  • Lý thuyết cơ học phi tuyến và tuyến tính hóa đối ngẫu: Xây dựng tiêu chuẩn sai số bình phương trung bình đối ngẫu kết hợp hàm trọng số $p(\mu)$ phụ thuộc vào mức độ phụ thuộc tuyến tính $\mu$, thiết lập hệ tương đương tuyến tính thay thế hệ phi tuyến dạng Duffing và khe hở tự do.
  • Lý thuyết ma trận và đại số đa thức: Sử dụng định thức ma trận Sylvester (kết thức $R_{p,q} = 0$) tách từ phần thực $p(\omega)$ và phần ảo $q(\omega)$ của đa thức đặc trưng để tìm vận tốc tới hạn flutter $U_f$ và tần số dao động tới hạn $\omega$.
  • Phương pháp thực tế:
    • Phương pháp giải tích: Xây dựng hệ phương trình vi phân chuyển động 2 bậc tự do, phương trình trạng thái 5 biến cho hệ có điều khiển PID, thiết lập công thức nghiệm giải tích cho vận tốc flutter và tần số dao động.
    • Phương pháp số: Sử dụng thuật toán Runge-Kutta giải hệ phương trình vi phân phi tuyến để đối chứng nghiệm; sử dụng phần mềm ANSYS Fluent giải hệ phương trình Navier-Stokes trung bình Reynolds (RANS) mô phỏng dòng khí qua profile cánh.
    • Phương pháp tối ưu hóa: Áp dụng thuật toán lập trình toàn phương liên tiếp (Sequential Quadratic Programming - SQP) để tối ưu hóa hình dạng khí động học của thiết diện cánh.
  • Nguồn tư liệu và dữ liệu: Thông số cơ lý hợp kim nhôm 7075-T6; hình học profile cánh NACA 2412 và Eppler 66; trường số liệu mô phỏng CFD tại số Reynolds $Re = 5 \times 10^5$ với các góc tới $\alpha = 0^\circ, 5^\circ, 10^\circ, 18^\circ$; dữ liệu đối sánh từ công trình của Li và cộng sự (2011).

Nội dung chính theo từng chương

Chương 1: Tổng quan về bài toán phân tích đáp ứng của thiết diện cánh

Chương 1 hệ thống hóa các khái niệm cơ bản về khí đàn hồi thông qua tam giác Collar. Tác giả phân tích các cấp độ mô hình hóa cánh máy bay: mô hình 2 chiều mặt cắt điển hình gắn trên hai lò xo, mô hình cánh dạng tấm ngàm một đầu, mô hình khí động học tính bằng CFD và mô hình kết hợp FSI (CFD-CSM). Luận án phân tích cơ chế hình thành dao động vòng giới hạn (LCO) và hiện tượng rẽ nhánh Hopf do tính phi tuyến kết cấu tập trung gây ra (gồm phi tuyến độ cứng bậc ba và phi tuyến khe hở tự do). Chương này làm rõ nguồn gốc của phương pháp tuyến tính hóa đối ngẫu và xác định phạm vi nghiên cứu của luận án là tập trung vào mô hình 2 bậc tự do tựa dừng với luật điều khiển PID.

Chương 2: Mô hình cơ học của thiết diện cánh chuyển động trong dòng khí

Chương 2 tập trung thiết lập các phương trình cơ học chi tiết:

  • Dẫn xuất lực khí động dừng theo lý thuyết dải xoáy mỏng: $$L = \pi \rho c U^2 \alpha$$ $$M = \frac{1}{4} \pi \rho c^2 U^2 \alpha$$ Biểu thức mô men tại vị trí cách biên dẫn khoảng cách $x$ là $M(x) = (x - c/4)\pi \rho c U^2 \alpha$, chứng minh tâm khí động nằm ở vị trí $1/4$ dây cung ($x = c/4$).
  • Mở rộng sang mô hình khí động tựa dừng với chuyển vị thẳng đứng $h$ và góc xoay $\alpha$: $$L = \frac{1}{2} \rho c U^2 \left[ C_L + c_l \left( \frac{\dot{h}}{U} + \left(\frac{3c}{4} - x_0\right)\frac{\dot{\alpha}}{U} \right) \right]$$ $$M = \frac{1}{4} \rho c^2 U^2 \left[ C_M + c_m \left( \frac{\dot{h}}{U} + \left(\frac{3c}{4} - x_0\right)\frac{\dot{\alpha}}{U} \right) \right]$$
  • Thiết lập hệ phương trình vi phân chuyển động 2 bậc tự do uốn - xoắn: $$\begin{cases} m\ddot{h} + S\ddot{\alpha} + k_h h = -L \ I\ddot{\alpha} + S\ddot{h} + k_\alpha \alpha = M \end{cases}$$
  • Phân tích điều kiện mất ổn định:
    • Mất ổn định 1 bậc tự do: hiện tượng stall khi tăng góc tới làm giảm lực nâng gây mất ổn định uốn; mất ổn định xoắn tĩnh xảy ra khi $b^2 U^2 c_m \ge 2 k_\alpha$.
    • Mất ổn định 2 bậc tự do: thiết lập phương trình đa thức trùng phương $a_4 p^4 + a_2 p^2 + a_0 = 0$. Ranh giới ổn định xác định bởi điều kiện $a_2^2 - 4a_4 a_0 = 0$, dẫn đến phương trình bậc 4 trùng phương theo vận tốc $U$: $$b_4 \left(\frac{U^2}{2}\right)^2 + b_2 \left(\frac{U^2}{2}\right) + b_0 = 0$$ Vận tốc tới hạn flutter được xác định chính xác theo nghiệm giải tích: $$U_f = \sqrt{2 \cdot \frac{-b_2 - \sqrt{b_2^2 - 4b_0 b_4}}{2b_4}}$$ Khi mô men tĩnh $S \le 0$ (trọng tâm cánh nằm trước trục đàn hồi), nghiệm $U_f^2$ âm, chứng minh không xuất hiện flutter (theo kết quả Pines 1958).
  • Thiết lập mô hình không gian trạng thái 5 chiều cho hệ cánh 2 chiều tích hợp điều khiển PID qua cánh nhỏ: $$\mathbf{M}\dot{\mathbf{z}} = \mathbf{A}(U)\mathbf{z}$$ Đa thức đặc trưng bậc 5 được tách thành 2 phương trình thực và ảo theo $\omega$: $$\begin{cases} (m_2^2 - m_1 m_3)\omega^4 + (k_4 k_{10} - m_2 k_2 - m_2 k_6 - m_1 k_7 - k_5 k_9 - k_1 m_3)\omega^2 + (k_6 k_2 - k_3 k_9 - k_8 k_4 + k_1 k_7) = 0 \ (m_2 k_9 - k_4 m_3 - m_2 k_5 - m_1 k_{10})\omega^4 + (k_7 k_4 - k_1 k_{10} - k_6 k_5 + m_1 k_8 + m_2 k_3 - k_2 k_9)\omega^2 - (k_1 k_8 - k_6 k_3) = 0 \end{cases}$$ Vận tốc tới hạn $U_f$ được xác định bằng cách giải phương trình triệt tiêu kết thức Sylvester giữa hai đa thức trên; tần số dao động LCO $\omega$ được tính bằng cách khử $\omega^4$.
  • Phân tích CFD trên profile NACA 2412 và Eppler 66 bằng ANSYS Fluent tại $Re = 5 \times 10^5$, phân tích trường áp suất, trường vận tốc, ứng suất tương đương trên vật liệu hợp kim nhôm 7075-T6 và tối ưu hóa biên dạng cánh bằng thuật toán SQP.

Chương 3: Phát triển kỹ thuật đối ngẫu cho bài toán dao động phi tuyến

Chương 3 trình bày cơ sở toán học của phương pháp tuyến tính hóa tương đương:

  • So sánh các tiêu chuẩn kinh điển: tiêu chuẩn sai số phương trình (minS), tiêu chuẩn sai số thế năng, tiêu chuẩn tương đương điều chỉnh.
  • Xây dựng tiêu chuẩn đối ngẫu có trọng số: tích hợp hệ số trọng số $p$ kết hợp giữa tiêu chuẩn kinh điển và tiêu chuẩn đối ngẫu.
  • Đề xuất 3 cải tiến lựa chọn hàm trọng số $p(\mu)$ phụ thuộc vào mức độ phụ thuộc tuyến tính $\mu$:
    • Cải tiến 1: Xây dựng hàm trọng số dựa trên đặc tính sai số cực tiểu.
    • Cải tiến 2: Xây dựng hàm trọng số tuyến tính từng đoạn theo mức độ phụ thuộc tuyến tính $\mu$.
    • Cải tiến 3: Xây dựng hàm trọng số giải tích tối ưu hóa năng lượng tiêu tán.
  • Kiểm chứng các cải tiến trên phương trình phi tuyến dạng Duffing bậc cao dao động tự do và dao động ngẫu nhiên, chứng minh sai số tần số và độ lệch chuẩn dịch chuyển của các cải tiến đối ngẫu đều nhỏ hơn tiêu chuẩn kinh điển khi so sánh với nghiệm giải tích chính xác.

Chương 4: Áp dụng kỹ thuật tuyến tính hóa đối ngẫu cho bài toán phân tích đáp ứng phi tuyến của thiết diện cánh

Chương 4 áp dụng kỹ thuật tuyến tính hóa đối ngẫu cải tiến để giải hệ phương trình vi phân chuyển động của thiết diện cánh 2 bậc tự do có tính phi tuyến kết cấu và hệ điều khiển PID. Tác giả khảo sát chi tiết 5 ví dụ tính toán:

  • Ví dụ 1: Xác định biên flutter phi tuyến và đường cong biên độ - vận tốc cho hệ cánh 2 bậc tự do phi tuyến bậc ba. Kết quả giải tích theo kỹ thuật đối ngẫu trùng khớp với kết quả mô phỏng số Runge-Kutta và đường biên flutter trong công trình của Li và cộng sự (2011).
  • Ví dụ 2: Khảo sát đường cong tần số - vận tốc và so sánh vận tốc flutter ở các mức biên độ dao động khác nhau.
  • Ví dụ 3: Phân tích ảnh hưởng của tham số điều khiển tỉ lệ ($K_p = -0.5$ và $K_p = 0.5$). Khi $K_p = -0.5$, vận tốc flutter giảm; khi $K_p = 0.5$, vận tốc flutter được nâng cao đáng kể.
  • Ví dụ 4: Phân tích tham số điều khiển tích phân ($K_i = -2\text{ s}^{-1}$ và $K_i = 2\text{ s}^{-1}$). Luận án chỉ ra khi $K_i = 2\text{ s}^{-1}$ với điều kiện đầu $h(0)=0\text{ m}, \alpha(0)=0.1\text{ rad}$, hệ xuất hiện sự phân kỳ của dao động xoắn.
  • Ví dụ 5: Khảo sát tham số điều khiển vi phân ($K_d = -0.5\text{ s}$ và $K_d = 0.5\text{ s}$). Luận án chứng minh hệ số $K_d > 0$ làm tăng mức cản tương đương, giúp mở rộng miền vận tốc ổn định của cánh.

Tổng hợp các bảng thông số và dữ liệu tính toán chính trong luận án:

Ký hiệu bảng trong luận án Nội dung số liệu tổng hợp
Bảng 1 Các đặc điểm cơ lý của Hợp kim nhôm 7075-T6
Bảng 2 Các đặc tính khí động lực trên thiết diện cánh tại $Re = 5 \times 10^5$ (lực nâng, lực cản theo góc tới)
Bảng 4 So sánh tần số dao động tính toán theo các tiêu chuẩn tuyến tính hóa với tần số chính xác
Bảng 5 So sánh các dịch chuyển bình phương trung bình
Bảng 6, 7, 8, 9, 10 Các bộ thông số đầu vào cho 5 ví dụ tính toán (Ví dụ 1 đến Ví dụ 5)
Bảng 11, 12, 13, 14 So sánh vận tốc flutter tính theo tiêu chuẩn kinh điển, tiêu chuẩn đối ngẫu và phương pháp số cho Ví dụ 1 và Ví dụ 2
Bảng 15, 16 So sánh vận tốc flutter trong Ví dụ 3 với $K_p = -0.5$ và $K_p = 0.5$
Bảng 17, 18 So sánh vận tốc flutter và các giá trị riêng tính theo các phương pháp tuyến tính hóa trong Ví dụ 4 ($K_i = -2\text{ s}^{-1}$ và $K_i = 2\text{ s}^{-1}$)
Bảng 19, 20 So sánh vận tốc flutter trong Ví dụ 5 với $K_d = -0.5\text{ s}$ và $K_d = 0.5\text{ s}$

Kết quả và những đóng góp mới

Luận án tự tổng hợp các kết quả và đóng góp mới về mặt khoa học và thực tiễn như sau:

Về mặt lý luận và khoa học:

  • Mở rộng thành công phương pháp luận của cách tiếp cận đối ngẫu và tiêu chuẩn tuyến tính hóa tương đương đối ngẫu có trọng số sang lớp bài toán dao động tuần hoàn và phân tích mất ổn định flutter phi tuyến của thiết diện cánh máy bay chịu lực khí động.
  • Đề xuất 3 kỹ thuật cải tiến lựa chọn tham số trọng số $p(\mu)$ cho tiêu chuẩn tuyến tính hóa đối ngẫu, khắc phục nhược điểm của tiêu chuẩn tuyến tính hóa kinh điển và làm tăng độ chính xác của nghiệm giải tích đối với các hệ phi tuyến mạnh.
  • Thiết lập quy trình giải tích xác định vận tốc tới hạn flutter và tần số dao động vòng giới hạn (LCO) cho hệ cơ học tự dao động mở rộng 5 biến trạng thái (kết hợp điều khiển PID) thông qua việc giải kết thức ma trận Sylvester.

Về mặt thực tiễn và ứng dụng:

  • Cung cấp các công thức giải tích tường minh cho phép dự báo nhanh ranh giới flutter và đáp ứng LCO của cánh máy bay mà không cần thực hiện các chuỗi mô phỏng số phương trình vi phân phức tạp và tốn kém thời gian.
  • Khảo sát định lượng ảnh hưởng của các hệ số điều khiển PID ($K_p, K_i, K_d$) đối với cánh nhỏ, xác định được các vùng tham số gây mất ổn định phân kỳ (như trường hợp $K_i = 2\text{ s}^{-1}$) và các vùng tham số giúp mở rộng giới hạn vận tốc bay an toàn.
  • Cung cấp dữ liệu mô phỏng CFD chi tiết (phân bố áp suất, phân bố vận tốc, ứng suất) và kết quả tối ưu hóa hình học cánh bằng thuật toán SQP trên các profile NACA 2412 và Eppler 66 phục vụ thiết kế các mô hình máy bay không người lái (UAV).

Kiến nghị giải pháp:

  • Ứng dụng kỹ thuật tuyến tính hóa đối ngẫu cải tiến vào giai đoạn thiết kế sơ bộ kết cấu cánh hàng không nhằm đánh giá nhanh độ ổn định khí đàn hồi.
  • Tích hợp cánh nhỏ điều khiển PID với thông số vi phân $K_d > 0$ và tỉ lệ $K_p > 0$ hợp lý để kiểm soát và dập tắt dao động flutter.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế:

  • Mô hình cơ học của luận án mới dừng lại ở thiết diện cánh 2 chiều phẳng (2 bậc tự do uốn - xoắn), chưa mô hình hóa đầy đủ hình học cánh 3 chiều có tính đến biến dạng uốn dọc sải cánh và hiệu ứng đầu mút cánh.
  • Lực khí động chỉ khảo sát trong phạm vi lý thuyết dừng và tựa dừng cho dòng khí không nén được, chưa xét đến mô hình khí động không dừng tổng quát ở dải tần số cao hoặc dòng khí nén được (cận âm, siêu âm).
  • Tính phi tuyến kết cấu chỉ xem xét ở dạng phi tuyến tập trung (bậc ba và khe hở tự do), chưa mở rộng cho phi tuyến phân bố dọc sải cánh và phi tuyến vật liệu composite.
  • Thuật toán điều khiển mới dừng lại ở mức điều khiển PID cổ điển.

Hướng nghiên cứu tiếp:

  • Phát triển kỹ thuật đối ngẫu cho mô hình cánh 3 chiều kết hợp tương tác lưu chất - kết cấu đầy đủ (FSI 3D).
  • Nghiên cứu mô hình lực khí động không dừng trong dòng khí nén được.
  • Ứng dụng các thuật toán điều khiển phi tuyến hiện đại (như điều khiển thích nghi, điều khiển trượt) để kiểm soát flutter phi tuyến.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo hữu ích cho các đối tượng và lĩnh vực cụ thể:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Cơ kỹ thuật, Cơ học vật rắn, Kỹ thuật Hàng không và Cơ điện tử: Tham khảo phương pháp luận tuyến tính hóa đối ngẫu có trọng số trong Chương 3 để giải quyết các bài toán dao động phi tuyến phức tạp trong cơ học kết cấu.
  • Kỹ sư thiết kế khí động và kết cấu hàng không (đặc biệt trong lĩnh vực UAV): Tham khảo quy trình tính toán vận tốc flutter giải tích, kỹ thuật thiết lập kết thức Sylvester cho hệ điều khiển PID trong Chương 2 và Chương 4, cũng như phương pháp mô phỏng CFD kết hợp tối ưu hóa SQP biên dạng cánh.
  • Các nhà nghiên cứu về cầu treo và kết cấu nhịp lớn: Tham khảo mô hình phân tích tương tác khí đàn hồi 2 bậc tự do để khảo sát hiện tượng mất ổn định khí động của dầm chủ cầu dây văng, cầu treo nhịp lớn.

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

1. Tại sao theo lý thuyết cánh mỏng dừng, tâm khí động của thiết diện cánh luôn nằm ở vị trí $1/4$ chiều dài dây cung tính từ biên dẫn? Theo lý thuyết cánh mỏng dừng của Glauert, khi phân bố xoáy $\gamma(\theta)$ được xác định và tích phân lực nâng $L = \pi \rho c U^2 \alpha$ cùng mô men tại biên dẫn $M_{l.e.} = -\frac{1}{4}\pi \rho c^2 U^2 \alpha$, mô men khí động tại điểm cách biên dẫn khoảng cách $x$ có dạng $M(x) = xL + M_{l.e.} = (x - c/4)\pi \rho c U^2 \alpha$. Tại vị trí $x = c/4$, mô men này triệt tiêu hoàn toàn ($M(c/4) = 0$), do đó vị trí $1/4$ dây cung được định nghĩa là tâm khí động của thiết diện cánh.

2. Điều kiện nào về mặt cơ học đảm bảo hệ thiết diện cánh 2 bậc tự do không bao giờ xuất hiện hiện tượng mất ổn định flutter? Dựa trên công thức giải tích vận tốc flutter $U_f$ rút ra từ phương trình trùng phương bậc 4, khi mô men tĩnh $S = \int r dm \le 0$ (nghĩa là trọng tâm của thiết diện cánh nằm trùng hoặc nằm phía trước trục đàn hồi theo chiều dòng khí), tất cả các nghiệm của phương trình vận tốc flutter đều âm hoặc ảo. Khi đó, sự ghép nối năng lượng giữa dao động uốn và dao động xoắn không tạo ra cản âm, do đó hiện tượng flutter không xảy ra (kết quả phù hợp với chứng minh của Pines 1958 và Dowell 2015).

3. Phương pháp kết thức Sylvester được áp dụng như thế nào để xác định vận tốc flutter và tần số LCO trong hệ có điều khiển PID? Hệ thiết diện cánh 2 chiều tích hợp điều khiển PID được biểu diễn qua phương trình trạng thái 5 chiều $\mathbf{M}\dot{\mathbf{z}} = \mathbf{A}(U)\mathbf{z}$. Đa thức đặc trưng bậc 5 tại trạng thái tới hạn (giá trị riêng thuần ảo $i\omega$) được tách thành hai đa thức thực $p(\omega)$ và ảo $q(\omega)$. Điều kiện để hai đa thức có nghiệm chung $\omega$ là định thức của ma trận Sylvester tương ứng (kết thức $R_{p,q}$) phải bằng 0. Phương trình $R_{p,q}(U) = 0$ là phương trình đơn biến xác định vận tốc tới hạn $U_f$. Tần số $\omega$ sau đó được xác định trực tiếp bằng cách khử $\omega^4$ giữa hai phương trình thực và ảo.

4. Luận án này có điểm gì khác biệt cốt lõi so với luận án của tác giả Nguyễn Ngọc Linh (2015)? Luận án của Nguyễn Ngọc Linh (2015) nghiên cứu phát triển tiêu chuẩn đối ngẫu có trọng số cho bài toán dao động ngẫu nhiên phi tuyến của hệ 1 bậc tự do chịu kích động ồn dừng Gauss trung bình không. Ngược lại, luận án này mở rộng kỹ thuật đối ngẫu sang lớp bài toán dao động tuần hoàn và hiện tượng mất ổn định flutter phi tuyến của thiết diện cánh 2 bậc tự do chịu lực khí động, đồng thời đề xuất 3 cách tiếp cận mới để xác định tham số trọng số $p(\mu)$ cho hệ tự dao động có điều khiển.

5. Tiêu chuẩn ổn định của dao động vòng giới hạn (LCO) đối với các nhiễu biên độ dao động được xác định bằng điều kiện nào? Theo tiêu chuẩn ổn định của Wei và cộng sự (2014) được luận án áp dụng, một dao động LCO có biên độ $A$ là ổn định nếu tích giữa sự biến thiên phần thực của giá trị riêng $\Delta\sigma$ và độ nhiễu biên độ $\Delta A$ thỏa mãn bất đẳng thức: $$\Delta\sigma \cdot \Delta A < 0$$ Điều kiện này chỉ ra rằng khi biên độ tăng ($\Delta A > 0$), phần thực của giá trị riêng trở nên âm ($\Delta\sigma < 0$) làm tiêu tán năng lượng đưa biên độ về vị trí cân bằng, và ngược lại.

Kết luận

Luận án tiến sĩ "Phân tích đáp ứng của profile cánh máy bay theo cách tiếp cận đối ngẫu" của tác giả Nguyễn Minh Triết (chuyên ngành Cơ kỹ thuật, mã số 62 52 01 01) đã phát triển và mở rộng thành công phương pháp tuyến tính hóa tương đương đối ngẫu cho bài toán dao động phi tuyến tuần hoàn và hiện tượng mất ổn định flutter của thiết diện cánh máy bay. Thông qua việc xây dựng 3 cải tiến cho hàm trọng số và ứng dụng kết thức Sylvester, luận án đã thiết lập được các biểu thức giải tích chính xác để dự báo vận tốc flutter tới hạn và tần số dao động vòng giới hạn cho hệ 2 bậc tự do có điều khiển PID. Các kết quả giải tích đã được kiểm chứng tin cậy thông qua so sánh với mô phỏng số Runge-Kutta, các công bố quốc tế và phân tích khí động CFD trên phần mềm ANSYS Fluent, đóng góp phương pháp luận có giá trị cho ngành cơ học tính toán và kỹ thuật hàng không.