Tổng quan về luận án

Hiện tượng khí động đàn hồi (Aeroelasticity) là một trong những bài toán phức tạp và then chốt nhất trong cơ học công trình và kỹ thuật hàng không vũ trụ hiện đại. Bản chất của hiện tượng này nằm ở sự tương tác phi tuyến, tương hỗ liên tục giữa ba trường lực: lực khí động học (Aerodynamic forces), lực đàn hồi cấu trúc (Elastic forces) và lực quán tính (Inertial forces). Khi kết cấu cánh khí cụ bay hoạt động trong dòng lưu chất, tải khí động gây biến dạng hình học của cánh (uốn, xoắn), sự biến dạng này lập tức làm biến đổi trường phân bố áp suất và tải khí động bao quanh, từ đó kích hoạt một chu trình phản hồi liên tục. Hậu quả cực đoan của hiện tượng khí động đàn hồi động như rung lắc tự kích (Flutter) hay khí động đàn hồi tĩnh như xoắn phá hủy (Torsional Divergence) và đảo chiều tác dụng cánh lái (Control Surface Reversal) có thể dẫn tới sự phá hủy tức thời toàn bộ cấu trúc máy bay.

Mặc dù nền tảng lý thuyết khí động đàn hồi đã được định hình qua nhiều thập kỷ, phần lớn các công trình kinh điển vẫn tập trung vào các dạng hình học cơ bản hoặc giải pháp tuyến tính hóa truyền thống. Thực tế thiết kế khí cụ bay thế hệ mới đòi hỏi việc ứng dụng các biên dạng đặc thù nhằm nâng cao hiệu suất bay và tính tàng hình, tiêu biểu như biên dạng siêu tới hạn (Supercritical airfoil) với lưng cánh phẳng và mép thoát cong cụp xuống, hoặc biên dạng cánh Delta mũi tên góc quét lớn hoạt động ở góc tấn cao. Khoảng trống nghiên cứu (Research gap) cấp thiết hiện nay là sự thiếu hụt các mô hình tính toán tương tác lỏng - rắn hai chiều (Two-way Fluid-Structure Interaction - FSI) phản ánh trung thực bản chất vật lý phi tuyến, song hành với các công cụ số chuyên dụng (như phương pháp biên nhúng - Immersed Boundary Method) và dữ liệu thực nghiệm đối sánh đa chiều cho các hình học cánh đặc thù.

Luận án tiến sĩ chuyên ngành Cơ học (mã số: 9440109) của nghiên cứu sinh Trần Ngọc Khánh, thực hiện tại Đại học Bách khoa Hà Nội dưới sự hướng dẫn khoa học của PGS. Nguyễn Phú Khánh và PGS. Hoàng Thị Kim Dung với tiêu đề: "Nghiên cứu đặc tính khí động đàn hồi cánh khí cụ bay có biên dạng đặc thù", đã giải quyết trực diện bài toán khoa học này. Luận án đặt ra ba câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. RQ1 & H1: Cơ chế tương tác khí động đàn hồi trên cánh có biên dạng siêu tới hạn (Supercritical airfoil) khác biệt như thế nào so với cánh biên dạng đối xứng chuẩn (NACA65A004) khi xét đến sự dịch chuyển của sóng va và phân bố áp suất phi tuyến? Giả thuyết đặt ra là biên dạng siêu tới hạn dịch chuyển tâm khí động (AC) về phía sau, làm thay đổi đáng kể vận tốc tới hạn xoắn phá hủy ($V_{div}$) và tần số dao động Flutter.
  2. RQ2 & H2: Phương pháp tính toán tương tác chất lưu - kết cấu hai chiều (2-way FSI) và phương pháp biên nhúng (IBM) có khả năng mô phỏng chính xác đáp ứng động học và hiện tượng rung lắc cánh Delta (Wing-Rock) ở góc tấn lớn so với dữ liệu thực nghiệm ống khí động hay không?
  3. RQ3 & H3: Hình dạng bình diện cánh (chữ nhật, hình thang, Delta) kết hợp với cấu trúc vật liệu dị hướng (Composite Carbon-Epoxy) tạo ra quy luật phân bố ứng suất và biến dạng đàn hồi như thế nào dưới tải trọng động lực học dòng khí?

Khung lý thuyết của luận án tích hợp tam giác tương tác Collar (1946), phương trình Navier-Stokes cho dòng nhớt kết hợp mô hình rối SST (Shear Stress Transport) và phương trình động lực học kết cấu đàn hồi phi tuyến. Phạm vi nghiên cứu bao quát các mẫu cánh chuẩn AGARD 445.6, cánh phẳng hình chữ nhật và hình thang mang biên dạng NACA65A004, cánh mang biên dạng siêu tới hạn Whitcomb-IL (độ dày cực đại 11% tại 35% dây cung, độ vồng 2.0%), cùng cánh Delta ở vận tốc thấp ($V = 10\text{ m/s}$) và dải góc tấn mở rộng ($\alpha = 0^\circ \div 30^\circ$).

Literature Review và Positioning

Lịch sử nghiên cứu khí động đàn hồi ghi nhận bước ngoặt lý thuyết nền tảng từ sơ đồ tam giác tương tác của Collar (1946), phân định rõ ranh giới giữa Cơ học bay (khí động - quán tính), Dao động cơ học (đàn hồi - quán tính), Khí động đàn hồi tĩnh (khí động - đàn hồi) và Khí động đàn hồi động (khí động - đàn hồi - quán tính). Trước đó, các nghiên cứu giải tích của Prandtl (1918), Ackermann (1918) và Birnbaum (1922) đã phát triển lý thuyết cánh tĩnh sang cánh dao động động lực học dựa trên định lý hoàn lưu Kutta - Joukowsky (1906). Küssner (1929) tiếp tục đặt nền móng giải tích cho bài toán Flutter hai bậc tự do (uốn - xoắn). Những thảm họa lịch sử như sự sụp đổ của cầu treo Tacoma Narrows (Mỹ, 1940) hay tai nạn gãy cánh máy bay ném bom Handley-Page 0/400 (Anh, 1918) và máy bay một tầng cánh của Langley (1903) được Garrick và Reed III (1981) tổng kết đã chứng minh tính chất phá hủy thảm khốc của hiện tượng mất ổn định khí động đàn hồi khi độ cứng chống xoắn ($GJ$) không đủ đáp ứng áp suất động của dòng chảy.

Trong y văn quốc tế, hai trường phái tiếp cận chủ đạo tồn tại nhiều tranh luận:

  • Trường phái giải tích - bán thực nghiệm cổ điển: Sử dụng các phương pháp K, P, P-K kết hợp khí động tuyến tính (Theodorsen, Doublet Lattice Method). Phương pháp này tiết kiệm chi phí tính toán nhưng bộc lộ sai số nghiêm trọng khi áp dụng cho cánh có hình học phức tạp, cánh siêu tới hạn có sóng va (Shock wave) hoặc cánh Delta xuất hiện các cấu trúc xoáy tách dòng phi tuyến (Vortex-induced separation).
  • Trường phái mô phỏng số đa vật lý (Multiphysics FSI): Tích hợp trực tiếp Computational Fluid Dynamics (CFD) và Computational Structural Dynamics (CSD). Các công trình chuẩn của Yates (1987) tại NASA Langley trên mẫu cánh chuẩn AGARD 445.6 đã cung cấp bộ dữ liệu thực nghiệm kinh điển về vận tốc và tần số Flutter từ dải dưới âm ($M = 0.499$) đến siêu âm ($M = 1.141$). Gần đây, các nghiên cứu quốc tế của Badcock et al. và Farhat et al. tập trung vào ghép nối lưới chuyển động (Dynamic Meshing ALE), tuy nhiên chi phí tính toán tái tạo lưới rất lớn và dễ gặp lỗi phân kỳ khi cấu trúc biến dạng lớn.

Luận án của NCS. Trần Ngọc Khánh định vị chính xác vào khoảng trống khoa học này: kết hợp phương pháp FSI hai chiều lặp không tường minh (2-way Iteratively Implicit FSI) trên nền tảng phần mềm công nghiệp ANSYS (Fluent & Mechanical) với việc phát triển một chương trình chuyên dụng tự xây dựng dựa trên Phương pháp biên nhúng (Immersed Boundary Method - IBM). Mô hình này cho phép biểu diễn vật thể biến dạng trên hệ tọa độ phi cấu trúc Lagrange di động lồng trong lưới Euler cố định của chất lưu, loại bỏ hoàn toàn yêu cầu tái tạo lưới đắt đỏ. Nghiên cứu thực hiện đối sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu: thực nghiệm cánh chuẩn AGARD 445.6 của Yates (1987) tại NASA và các nghiên cứu mô phỏng cánh siêu tới hạn của Harris (NASA TM-81907), chứng minh độ tin cậy vượt trội trong việc nắm bắt các hiện tượng phi tuyến cục bộ.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại những bước tiến lý thuyết quan trọng bổ sung vào hệ thống cơ học vật rắn biến dạng và khí động lực học:

  1. Mở rộng lý thuyết xoắn phá hủy tĩnh cho biên dạng 3D đặc thù: Kế thừa lý thuyết giải tích cánh mỏng, tác giả thiết lập công thức xác định áp suất động tới hạn cho cánh hình chữ nhật có phân bố góc xoắn biến thiên phi tuyến dọc sải cánh: $$q_{div} = \frac{\pi^2 GJ}{4 s^2 c e_c a_w} \quad \text{hoặc theo mô hình xấp xỉ bậc một:} \quad q_{div} = \frac{3GJ}{c e_c s^2 a_w}$$ Vận tốc xoắn phá hủy tương ứng được xác định chính xác theo đặc tính hình học và cơ lý tính vật liệu: $$V_{div} = \sqrt{\frac{2 q_{div}}{\rho}} = \sqrt{\frac{\pi^2 GJ}{2 \rho c e_c s^2 a_w}}$$ Trong đó, việc chứng minh tâm khí động ($AC$) dịch chuyển phụ thuộc vào độ vồng phía sau của biên dạng siêu tới hạn đã làm sáng tỏ cơ chế giảm biên độ ổn định xoắn tĩnh so với biên dạng đối xứng.
  2. Cơ chế đảo chiều tác dụng cánh lái: Xác lập biểu thức tương quan giải tích giữa áp suất động đảo chiều ($q_{rev}$) và áp suất động xoắn phá hủy ($q_{div}$): $$\frac{q_w}{q_{rev}} = \frac{e a_C - b_C}{e a_C} \implies q_{rev} = \frac{q_w e a_C}{e a_C - b_C}$$ Bằng chứng toán học chỉ ra hệ số $b_C$ luôn mang giá trị âm, dẫn đến hệ quả vận tốc đảo chiều cánh lái luôn nhỏ hơn vận tốc xoắn phá hủy kết cấu ($V_{rev} < V_{div}$), tạo cơ sở lý thuyết cho việc ấn định dải vận tốc an toàn của các bề mặt điều khiển bay.
  3. Mô hình hóa hiện tượng Wing-Rock phi tuyến trên cánh Delta: Chứng minh cơ chế dao động tự kích chu kỳ giới hạn (Limit Cycle Oscillation - LCO) quanh trục liệng sinh ra từ sự dịch chuyển và bất đối xứng của cặp cuộn xoáy mép vào (Leading-Edge Vortices - LEV) khi góc tấn $\alpha \ge 20^\circ$.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột phương pháp luận:

  • Tích hợp đa lý thuyết: Kết hợp lý thuyết dòng khí nhớt nén được Navier-Stokes, lý thuyết đàn hồi kết cấu Euler-Bernoulli/Timoshenko cho tấm vỏ mỏng, và lý thuyết dao động cơ học đa bậc tự do.
  • Tiếp cận phân tích số kép (Dual Numerical Approach): Sử dụng song song giải pháp FSI thương mại (ANSYS System Coupling) và thuật toán IBM chuyên dụng tự phát triển. Khung tính toán IBM giải phương trình Navier-Stokes trên lưới Euler Cartesian cố định, sử dụng hàm Dirac delta để truyền tương tác lực và vận tốc giữa các điểm nút Lagrange cấu trúc và mạng lưới chất lưu.
  • Điều kiện biên xác định tường minh: Cố định ngàm tại gốc cánh ($y=0 \implies \theta = 0$), biên tự do tại mút cánh ($y=s \implies d\theta/dy = 0$), dòng vào đồng nhất ($V_\infty, P_\infty$), và bề mặt tiếp xúc lỏng - rắn thỏa mãn điều kiện không trượt và liên tục ứng suất ($\mathbf{u}{fluid} = \mathbf{u}{solid}, ; \boldsymbol{\sigma}{fluid} \cdot \mathbf{n} = \boldsymbol{\sigma}{solid} \cdot \mathbf{n}$).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng (Positivism) và chủ nghĩa hiện thực phê phán (Critical Realism), kết hợp chặt chẽ giữa mô phỏng số tiên tiến và thực nghiệm vật lý trong phòng thí nghiệm. Thiết kế đa tầng (Multi-level design) bao gồm:

  • Tầng 1: Kiểm nghiệm độ chính xác (Validation Benchmark): Kiểm chứng mã nguồn CFD và CSD trên các bài toán chuẩn quốc tế gồm dòng chảy Coanda ống đuôi trực thăng NOTAR, dòng chảy qua trụ tròn ở $Re = 100$, phân tích tần số dao động riêng và hiện tượng Flutter trên cánh chuẩn AGARD 445.6.
  • Tầng 2: Khảo sát khí động tĩnh học (Static Aerodynamics): Đánh giá trường vận tốc, phân bố hệ số áp suất ($C_p$), hệ số lực nâng ($C_L$), lực cản ($C_D$) và mô-men ($C_m$) trên các mô hình cánh hình chữ nhật NACA65A004, hình thang NACA65A004, cánh siêu tới hạn Whitcomb-IL và cánh Delta.
  • Tầng 3: Khảo sát tương tác khí động đàn hồi động (Dynamic Aeroelasticity / FSI): Đánh giá đáp ứng chuyển vị thời gian thực, phổ tần số, ứng suất tương đương Von-Mises và phát hiện điểm bifurcation mất ổn định Flutter.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ tiêu chuẩn đo lường và tính toán nghiêm ngặt:

  1. Thiết lập mô phỏng FSI 2 chiều (ANSYS Coupling):
    • Module chất lưu (CFD): Giải phương trình Reynolds-Averaged Navier-Stokes (RANS) kết hợp mô hình rối SST $k-\omega$. Lưới khí động biên dạng lăng trụ tại lớp biên đảm bảo $y^+ \approx 1$.
    • Module kết cấu (CSD): Sử dụng phần tử vỏ (Shell) và khối (Solid) 3D bậc cao. Mô hình hóa đặc tính vật liệu hợp kim nhôm ($E = 70\text{ GPa}, \nu = 0.33, \rho = 2700\text{ kg/m}^3$) và vật liệu composite nhiều lớp Carbon-Epoxy-UD-230GPa-Prepreg ($E_x = 3.31\text{ GPa}, E_y = 3.15\text{ GPa}, G = 0.41\text{ GPa}, \nu = 0.31$).
    • Quy trình lặp không tường minh (2-way Iteratively Implicit): Trao đổi dữ liệu tải áp suất và chuyển vị nút qua giao diện tiếp xúc (Fluid-Solid Interface) sau mỗi bước thời gian $\Delta t = 10^{-3} \div 10^{-4}\text{ s}$ với tiêu chí hội tụ sai số tương đối đạt dưới $10^{-5}$.
  2. Quy trình thực nghiệm ống khí động:
    • Thử nghiệm tại ống khí động dưới âm (Mach $0 \div 0.3$) của Đại học Bách khoa Hà Nội.
    • Hệ thống đo đạc gồm: Bộ gá chuyên dụng đo lực khí động 6 thành phần, cảm biến áp suất vi sai đa kênh nối với các lỗ trích áp phân bố dọc dây cung mút cánh và gốc cánh ($x/c = 0.05 \div 0.95$), hệ thống đo dao động quang học/laser và gia tốc kế áp điện ghi nhận tần số vẫy và biên độ dao động mép vào/mép ra.
+-----------------------------------------------------------------------------------+
|                            QUY TRÌNH NGHIÊN CỨU FSI                               |
+-----------------------------------------------------------------------------------+
|  [CFD: Fluent (RANS/SST k-w)]  <--- Tải khí động (Pressure/Force) ---+            |
|               |                                                      |            |
|     Bước thời gian lặp (dt)                                          |            |
|               v                                                      |            |
|  [Giao diện tiếp xúc FSI]     ---> Dữ liệu chuyển vị (Displacement) -+            |
|               |                                                                   |
|               v                                                                   |
|  [CSD: Mechanical (FEA/Modal)]                                                    |
+-----------------------------------------------------------------------------------+

Data và phân tích

Dữ liệu thu nhận từ mô phỏng và thực nghiệm được xử lý bằng các thuật toán phân tích phổ biến đổi Fourier nhanh (FFT), khớp hàm phi tuyến và phân tích phương thức riêng (Modal Extraction):

  • Độ ổn định và hội tụ lưới: Thực hiện đánh giá độc lập lưới khí động học từ 1.2 triệu đến 4.8 triệu phần tử, chỉ số sai số GCI (Grid Convergence Index) đạt dưới 1.8%.
  • Dữ liệu kiểm nghiệm cánh AGARD 445.6: Tại số Mach cận âm $M = 0.499$, kết quả tính toán tần số dao động uốn mode 1 đạt $9.6\text{ Hz}$ (sai lệch $0.42%$ so với thực nghiệm NASA Yates $9.64\text{ Hz}$), tần số xoắn mode 2 đạt $38.2\text{ Hz}$ (sai lệch $1.03%$ so với thực nghiệm $38.6\text{ Hz}$).
  • Kiểm nghiệm mô hình trụ tròn ($Re = 100$): Hệ số lực cản trung bình $\bar{C}_D = 1.34$ và số Strouhal $St = 0.165$, hoàn toàn trùng khớp với các công bố chuẩn quốc tế của Williamson và Roshko.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Sự khác biệt đột phá về phân bố áp suất trên cánh siêu tới hạn: Dòng chảy qua biên dạng siêu tới hạn Whitcomb-IL tạo ra vùng áp suất âm trải rộng và đồng đều trên toàn bộ lưng cánh phẳng, làm triệt tiêu đỉnh hút áp suất nhọn ở mép vào vốn thấy ở cánh NACA65A004. Phần mép thoát uốn cong xuống tạo thêm lực nâng đuôi (Aft-camber lift), làm tăng hệ số lực nâng cực đại $C_{Lmax}$ thêm $18.4%$ so với cánh đối xứng ở cùng góc tấn $\alpha = 10^\circ$, đồng thời duy trì hệ số lực cản $C_D$ thấp ở dải vận tốc cận âm.
  2. Đáp ứng khí động đàn hồi và dịch chuyển điểm mất ổn định Flutter: Dưới tác dụng của lực khí động ở vận tốc dòng $V = 10\text{ m/s}$ và góc tấn $\alpha = 10^\circ$, cánh hình chữ nhật NACA65A004 có độ võng mút cánh cực đại đạt $14.2\text{ mm}$, trong khi cánh siêu tới hạn cùng kích thước có độ võng lên tới $18.6\text{ mm}$ do tổng lực nâng sinh ra lớn hơn. Ứng suất pháp cực đại tập trung tại vùng ngàm gốc cánh của cánh siêu tới hạn cao hơn $24.7%$ so với cánh NACA65A004. Điều này dẫn tới vận tốc tới hạn Flutter ($V_f$) của cánh siêu tới hạn thấp hơn cánh đối xứng khoảng $8.5%$, đòi hỏi cấu trúc phải được gia cường độ cứng chống xoắn.
  3. Hiện tượng tạo xoáy và ổn định lực nâng ở góc tấn siêu cao trên cánh Delta: Trên cánh Delta, hai cuộn xoáy mép vào (Leading Edge Vortices - LEV) xoay ngược chiều hút dòng khí năng lượng cao xuống lưng cánh, duy trì lực nâng tiếp tục tăng tuyến tính và không bị thất tốc (Stall) đột ngột ngay cả khi góc tấn vượt quá $\alpha = 20^\circ \div 30^\circ$. Tại $\alpha = 20^\circ, \phi = 5^\circ$, hệ số áp suất tại mặt cắt $x/c_r = 0.6$ ghi nhận đáy áp suất $C_p = -2.45$, khẳng định cơ chế tạo lực nâng xoáy vượt trội.
  4. Cơ chế rung lắc cánh tự kích (Wing-Rock) trên cánh Delta: Mô phỏng FSI và thực nghiệm chỉ rõ: ở góc tấn lớn ($\alpha \ge 20^\circ$), sự bất đối xứng về cường độ và vị trí của cặp xoáy LEV khi cánh bị nghiêng một góc xoay $\phi$ tạo ra mô-men lật không tắt dần. Hệ thống rơi vào trạng thái dao động tự kích chu kỳ giới hạn (Limit Cycle Oscillation) với tần số vẫy xác định từ thực nghiệm đạt $f = 2.35\text{ Hz}$ ở $V = 10\text{ m/s}$.
  5. Hiệu năng của phương pháp biên nhúng (IBM) trong bài toán kết cấu chuyển động: Chương trình tính toán IBM tự phát triển tái hiện chính xác trường dòng chảy xoáy phân tầng quanh vật thể chuyển động mà không gặp phải hiện tượng thoái hóa phần tử lưới (Mesh skewness), giảm thời gian xử lý tiền kỳ và tính toán tới $42%$ so với kỹ thuật FSI lưới động thông thường.
+-----------------------------------------------------------------------------------+
|                  SO SÁNH ĐẶC TÍNH KHÍ ĐỘNG VÀ ĐÀN HỒI (V = 10 m/s)                 |
+------------------------------------+--------------------+-------------------------+
| Thông số khảo sát                  | Cánh NACA65A004    | Cánh siêu tới hạn       |
+------------------------------------+--------------------+-------------------------+
| Hệ số lực nâng CL (alpha = 10 deg) | 0.62               | 0.74 (+19.3%)           |
| Chuyển vị mút cánh cực đại         | 14.2 mm            | 18.6 mm (+31.0%)        |
| Ứng suất tương đương tại gốc cánh  | 42.5 MPa           | 53.0 MPa (+24.7%)       |
| Vận tốc tới hạn Flutter dự báo     | Cao hơn chuẩn      | Thấp hơn ~8.5%          |
+------------------------------------+--------------------+-------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Đặt nền tảng cho việc hiệu chỉnh các mô hình khí động đàn hồi phi tuyến, bổ sung số hạng phi tuyến do xoáy tách dòng và hiệu ứng lưng phẳng vào phương trình chuyển động aeroelastic tổng quát.
  • Về mặt phương pháp luận: Khẳng định tính khả thi và độ tin cậy của thuật toán IBM tự xây dựng trong giải quyết các bài toán biên tương tác phức tạp, cung cấp công cụ mã nguồn mở kiểm soát hoàn toàn thuật toán cho cộng đồng nghiên cứu cơ học tính toán.
  • Về mặt thực tiễn thiết kế hàng không: Cung cấp bộ dữ liệu biên dạng và quy tắc thiết kế độ cứng kết cấu ($GJ, EI$) cho các viện thiết kế máy bay không người lái (UAV) trinh sát tầm cao, máy bay phản lực tốc độ cao sử dụng cánh siêu tới hạn hoặc cánh Delta, ngăn ngừa nguy cơ tai nạn do rung lắc cấu trúc.
  • Về mặt chính sách và công nghệ quốc phòng: Đóng góp trực tiếp vào năng lực tự chủ công nghệ tính toán mô phỏng khí cụ bay tiên tiến, giảm phụ thuộc vào các gói phần mềm chuyên dụng quân sự bị bảo mật từ nước ngoài.

Limitations và Future Research

Mặc dù đạt được những đóng góp khoa học đột phá, luận án thẳng thắn thừa nhận các giới hạn nội tại:

  1. Giới hạn số Mach trong thực nghiệm: Các đo đạc thực nghiệm khí động đàn hồi chủ yếu thực hiện trong ống khí động tốc độ thấp ($M < 0.3, V \le 30\text{ m/s}$) do điều kiện trang thiết bị phòng thí nghiệm trong nước, chưa thể tạo lập dòng khí cận âm cao ($M = 0.7 \div 0.85$) để kích hoạt toàn diện sóng va thực tế trên cánh siêu tới hạn.
  2. Mô hình kết cấu tuyến tính hóa vật liệu: Phân tích kết cấu CSD trong FSI chủ yếu giả định vật liệu làm việc trong miền đàn hồi tuyến tính (Hookean), chưa xét đến hiện tượng từ biến, mỏi phi tuyến hoặc phá hủy vi mô cục bộ trong vật liệu composite sau hàng triệu chu kỳ rung lắc.
  3. Bậc tự do thực nghiệm cánh Delta: Mô hình thực nghiệm cánh Delta tập trung khảo sát bậc tự do xoay quanh trục dọc (Roll degree of freedom), chưa tích hợp đầy đủ 6 bậc tự do tự do hoàn toàn trong không gian.

Chương trình nghiên cứu tiếp nối trong giai đoạn 5-10 năm tới bao gồm:

  • Mở rộng thuật toán FSI-IBM cho dòng khí nén được vận tốc cao (Transonic/Supersonic IBM) có xét đến tương tác giữa sóng va cục bộ và lớp biên (Shock-Wave/Boundary-Layer Interaction - SWBLI).
  • Nghiên cứu thực nghiệm hiện tượng Aeroelastic Tailored Wing sử dụng vật liệu composite dị hướng có hướng sợi tối ưu hóa để chủ động triệt tiêu rung lắc Flutter.
  • Phát triển hệ thống điều khiển phản hồi chủ động (Active Flutter Suppression) sử dụng bề mặt lái thông minh bằng vật liệu áp điện hoặc vi cơ điện tử (MEMS).

Tác động và ảnh hưởng

Công trình nghiên cứu tạo ra sức lan tỏa sâu rộng trong cả giới học thuật và công nghiệp:

  • Tác động học thuật: Các kết quả nghiên cứu đã được công bố trên các tạp chí khoa học chuyên ngành uy tín và kỷ yếu hội nghị cơ học toàn quốc. Phương pháp FSI và IBM được trích dẫn và ứng dụng làm tài liệu giảng dạy chuyên sâu cho học viên cao học, nghiên cứu sinh ngành Kỹ thuật Hàng không Vũ trụ và Cơ học Thủy khí.
  • Chuyển giao công nghiệp và ứng dụng thực tế: Các công thức giải tích và quy trình mô phỏng đã được chia sẻ ứng dụng tại các đơn vị nghiên cứu phát triển phương tiện bay không người lái (UAV) chiến thuật, hỗ trợ tối ưu hóa trọng lượng kết cấu cánh và nâng cao trần bay an toàn.
  • Lợi ích kinh tế - kỹ thuật: Giúp rút ngắn chu kỳ thử nghiệm ống khí động tốn kém từ hàng chục lần xuống còn vài lần kiểm chứng trọng điểm, tiết kiệm hàng trăm triệu đồng chi phí chế tạo mẫu thử thực nghiệm trong giai đoạn thiết kế sơ bộ.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới nghiên cứu Cơ học: Tiếp cận nguồn tài liệu học thuật hoàn chỉnh, hệ thống phương trình vi phân chuyển động giải tích tường minh, cùng bộ dữ liệu thực nghiệm chuẩn đối sánh cho cánh NACA65A004, Supercritical và Delta.
  • Kỹ sư R&D Hàng không và Cơ khí Động lực: Nhận được quy trình chuẩn hóa thiết lập bài toán FSI hai chiều lặp không tường minh trên ANSYS và bộ mã nguồn tính toán IBM để ứng dụng trực tiếp vào thiết kế cánh máy bay, tuabin gió, cánh quạt trực thăng.
  • Các cơ quan quản lý và Viện thiết kế quốc phòng: Sở hữu căn cứ khoa học định lượng về biên độ an toàn khí động đàn hồi để xây dựng tiêu chuẩn nghiệm thu kỹ thuật bay cho các khí cụ bay thế hệ mới.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Đó là việc mở rộng thành công lý thuyết xoắn phá hủy và đảo chiều cánh lái cổ điển của Collar và Bisplinghoff sang các hệ biên dạng đặc thù phi đối xứng (Supercritical Whitcomb-IL). Luận án đã giải tích hóa và định lượng được sự dịch chuyển của tâm khí động $AC$ kết hợp với lực nâng đuôi phi tuyến (Aft-camber lift), chứng minh mối quan hệ tất yếu dẫn tới sự suy giảm vận tốc tới hạn Flutter ($V_f$) so với biên dạng đối xứng chuẩn.
  2. Đột phá phương pháp luận so với các nghiên cứu trước đây? So với nghiên cứu FSI một chiều truyền thống của các tác giả trong nước hoặc phương pháp lưới động ALE đắt đỏ của Farhat et al., luận án đã xây dựng thành công chương trình tính toán chuyên dụng phương pháp biên nhúng (IBM) kết hợp kiểm chứng đa chiều với FSI 2 chiều lặp không tường minh (ANSYS Fluent-Mechanical Coupling) và thực nghiệm ống khí động, giải quyết triệt để bài toán biến dạng lớn mà không gây thoái hóa lưới.
  3. Phát hiện nào gây bất ngờ nhất về mặt số liệu? Phát hiện về sự gia tăng tới $31.0%$ chuyển vị biến dạng uốn mút cánh và $24.7%$ ứng suất pháp cục bộ tại ngàm gốc cánh trên cánh siêu tới hạn so với cánh NACA65A004 ở cùng điều kiện dòng vào ($V = 10\text{ m/s}, \alpha = 10^\circ$). Điều này chỉ ra rằng hiệu quả tạo lực nâng vượt trội của cánh siêu tới hạn đi kèm với sự suy giảm đáng kể độ cứng dự trữ khí động đàn hồi.
  4. Quy trình tái lập kết quả (Replication Protocol) có được cung cấp đầy đủ? Hoàn toàn đầy đủ và chi tiết. Luận án công bố tường minh toàn bộ thông số hình học cánh, thông số cơ lý tính vật liệu (hợp kim nhôm, composite nhiều lớp Carbon-Epoxy-UD-230GPa), các điều kiện biên chuyển vị/ngoại lực/gia tốc, mô hình rối SST $k-\omega$, bước thời gian tích phân ($\Delta t = 10^{-3}\text{ s}$), và sơ đồ bố trí các lỗ trích áp trong ống khí động.
  5. Chương trình nghiên cứu 10 năm được định hình ra sao? Chương trình 10 năm tập trung vào 3 mũi nhọn: (1) Nâng cấp mã nguồn IBM cho dòng nén được trên âm có sóng va mạnh; (2) Thiết kế tối ưu hóa cấu trúc Aeroelastic Tailoring đa mục tiêu bằng trí tuệ nhân tạo (AI-driven Optimization); (3) Tích hợp hệ thống kiểm soát chủ động hiện tượng Flutter (Active Flutter Control) trên mô hình bay thực tế.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán nghiên cứu đặc tính khí động đàn hồi trên các mô hình cánh khí cụ bay có biên dạng đặc thù (NACA65A004, siêu tới hạn Whitcomb-IL, cánh Delta) với độ tin cậy khoa học cao.
  2. Xây dựng và chuẩn hóa thành công quy trình mô phỏng tương tác lỏng - rắn hai chiều (2-way FSI) lặp không tường minh, đồng thời phát triển chương trình tính toán chuyên dụng sử dụng phương pháp biên nhúng (IBM) loại bỏ nhu cầu tái tạo lưới phức tạp.
  3. Chế tạo thành công hệ thống mẫu cánh thực nghiệm và thực hiện các đo đạc khí động, dao động vẫy Flutter trong ống khí động, cung cấp bộ dữ liệu kiểm chứng thực nghiệm tin cậy cho các mô hình số.
  4. Phát hiện và làm sáng tỏ cơ chế suy giảm vận tốc Flutter trên cánh siêu tới hạn do hiệu ứng phân bố áp suất phẳng lưng cánh, cùng cơ chế dao động tự kích Wing-Rock sinh ra từ cặp xoáy LEV bất đối xứng trên cánh Delta ở góc tấn cao.
  5. Đóng góp trực tiếp 05 kết quả lý thuyết và thực tiễn đột phá vào kho tàng cơ học hàng không Việt Nam, mở ra 03 hướng nghiên cứu mũi nhọn về khí động đàn hồi phi tuyến, vật liệu thông minh và điều khiển chủ động rung lắc kết cấu bay trong tương lai.