Tổng quan về luận án

Nghiên cứu cơ học bay phỏng sinh học (biomimetic flight mechanics) và thiết bị bay cánh vẫy siêu nhỏ (Flapping-Wing Micro Air Vehicles - FWMAV) đang tạo nên một bước chuyển dịch lớn trong kỹ thuật hàng không vũ trụ hiện đại. Ở dải số Reynolds thấp và trung bình ($Re = 10^3 - 10^4$), các phương tiện bay cánh bằng hoặc cánh quay truyền thống suy giảm nghiêm trọng hiệu suất khí động do lực cản nhớt lớn và hiện tượng phân ly dòng khí sớm. Ngược lại, các loài côn trùng trong tự nhiên như bướm đêm Manduca sexta, ruồi giấm Drosophila melanogaster hay ong vò vẽ Bombus sở hữu khả năng bay treo (hovering) và cơ động siêu việt nhờ tận dụng các cơ chế khí động không dừng (unsteady aerodynamics) kết hợp biến dạng đàn hồi thụ động (passive elastic deformation) của kết cấu cánh. Luận án tiến sĩ kỹ thuật của tác giả Lê Vũ Đan Thanh với đề tài "Nghiên cứu đàn hồi khí động của cánh vẫy kiểu cánh côn trùng sử dụng mô hình cơ hệ nhiều vật" thuộc chuyên ngành Cơ kỹ thuật (mã số 9.52.01.01), dưới sự hướng dẫn khoa học của TS. Nguyễn Anh Tuấn và PGS.TS. Đặng Ngọc Thanh tại Học viện Kỹ thuật Quân sự (2023), là công trình tiên phong giải quyết bài toán tương tác kết cấu - chất lưu (Fluid-Structure Interaction - FSI) phi tuyến cho cánh vẫy sinh học thông qua cách tiếp cận động lực học hệ nhiều vật (Multi-Body Dynamics - MBD).

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác định xuất phát từ sự đánh đổi cố hữu trong các công cụ mô phỏng hiện hành:

  1. Các mô hình giải tích bậc thấp dựa trên lý thuyết phần tử cánh dừng hoặc tựa dừng (Quasi-Steady Blade Element Theory - BET) của Ellington (1984) và Ansari et al. (2006) bỏ qua tính không dừng 3D và tương tác dòng xoáy phức tạp.
  2. Các mô hình tổng hợp mode dao động tuyến tính (Linear Modal Synthesis) của Roccia et al. (2013) và Schwab et al. (2014) chỉ áp dụng được cho giả thiết biến dạng vô cùng bé, không thể mô tả chính xác góc uốn lên tới 15° và góc xoắn đạt 30° của cánh thực tế (Combes & Daniel, 2003; Willmott & Ellington, 1997).
  3. Các hệ thống tính toán kết hợp bậc cao CFD-CSD (Computational Fluid Dynamics - Computational Structural Dynamics) giải hệ phương trình Navier-Stokes trung bình hóa theo Reynolds (RANS) cùng phần tử hữu hạn phi tuyến vỏ (Nonlinear Shell FEM) của Nakata & Liu (2012) hay phương pháp biên nhúng của Shahzad et al. (2016) đòi hỏi từ 24 đến 93 giờ tính toán cho chỉ một chu kỳ vẫy cánh ($T = 1/f$), khiến việc khảo sát quét tham số động học và tối ưu hóa đa mục tiêu kết cấu trở nên bất khả thi.
  4. Các công trình MBD sơ khởi (Vanella et al., 2009; Eldredge, 2010; Arora et al., 2018; Truong et al., 2021) mới dừng lại ở mô hình 2D, cánh cứng đơn đoạn hoặc chuyển động quay 1 bậc tự do đơn giản, chưa tích hợp đầy đủ không gian 3D với 3 góc quay Euler biến thiên tuần hoàn.

Để giải quyết triệt để khoảng trống này, luận án thiết lập 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Làm thế nào để xây dựng một mô hình toán học tích hợp FSI 3D có số bậc tự do thu gọn nhưng vẫn mô tả chính xác các đặc tính đàn hồi phi tuyến hình học (uốn, xoắn) của cánh côn trùng dưới kích động động học nhiều bậc tự do ở gốc cánh?
    • H1: Mô hình hóa cánh vẫy đàn hồi dưới dạng cơ hệ nhiều vật nối tiếp bởi các khớp lò xo quay độ cứng biến thiên (lumped-parameter multi-body spring system), tham số hóa từ dầm Euler-Bernoulli tương đương tối ưu hóa theo mô hình phần tử hữu hạn (FEM) cấu trúc gân - màng, sẽ tái hiện chính xác động học biến dạng lớn mà giảm số bậc tự do từ hàng chục nghìn xuống dưới 20 bậc tự do.
  • RQ2: Cơ chế tương tác giữa xoáy mép trước (Leading Edge Vortex - LEV), hiệu ứng xoay bổ sung (Added-rotation effect), khối lượng gia nhập (Added mass) và sự biến dạng thụ động của cánh ảnh hưởng như thế nào đến khả năng tạo lực nâng và hiệu suất năng lượng trong chế độ bay treo?
    • H2: Tích hợp phương pháp xoáy không dừng 3D (Extended 3D Unsteady Vortex-Lattice Method - UVLM) với mô hình lực hút mép trước Polhamus và mô hình khuếch tán lõi xoáy Lamb-Oseen sẽ bắt trọn chính xác cường độ LEV và hiện tượng bắt vết xoáy (wake capture) trên cánh biến dạng với chi phí tính toán giảm trên 95% so với CFD đầy đủ.
  • RQ3: Dải tỉ lệ tần số riêng không thứ nguyên và quy luật phân bố độ cứng chống uốn, chống xoắn tối ưu nào cho phép cực đại hóa hệ số lực nâng trung bình đồng thời cực tiểu hóa công suất tiêu thụ riêng?
    • H3: Tồn tại một vùng cộng hưởng khí động đàn hồi tối ưu xác định bởi tỉ lệ tần số $f_1^* = f_1/f \in [2.0, 3.0]$ kết hợp phân bố độ cứng giảm dần phi tuyến từ gốc đến mút cánh, tại đó năng lượng biến dạng đàn hồi được tái hấp thu để hỗ trợ đảo chiều vẫy cánh hiệu quả nhất.

Đối tượng và phạm vi khảo sát thực nghiệm - số trị của luận án tập trung vào cánh trước và cánh đầy đủ của loài bướm đêm Manduca sexta với sải cánh bán phần $R = 48.3 - 51.9\text{ mm}$, diện tích cánh $S \approx 9.47\text{ cm}^2$, khối lượng cánh $m_w \approx 30.5 - 46.0\text{ mg}$, hoạt động ở dải tần số vẫy $f = 20 - 30\text{ Hz}$ trong chế độ bay treo ổn định. Ý nghĩa khoa học và thực tiễn của công trình nằm ở việc xác lập một công cụ mô phỏng FSI nhanh, độ tin cậy cao, phục vụ trực tiếp quá trình thiết kế, chế tạo khung cánh composite siêu nhẹ cho các thiết bị bay trinh sát, cứu hộ quy mô micro và nano.

Literature Review và Positioning

Khí động lực học cánh vẫy sinh học đã trải qua bốn giai đoạn tiến hóa lý thuyết lớn:

  1. Giai đoạn bán tĩnh và phần tử cánh (1980 - 2000): Mở đầu bởi công trình nền tảng của Weis-Fogh (1973) về cơ chế đập và tách (clap-and-fling) và Ellington (1984) với chuỗi nghiên cứu phân tích lực khí động côn trùng dựa trên lý thuyết BET. Các mô hình này giả định góc tấn nhỏ và hệ số lực nâng $C_L = 2\pi\alpha$, hoàn toàn bất lực khi giải thích lực nâng vượt mức (lift enhancement) gấp 2-3 lần giới hạn lý thuyết tĩnh ở các góc tấn cực đại $\alpha \ge 40^\circ$.
  2. Giai đoạn khám phá các cơ chế xoáy không dừng (1996 - 2005): Đột phá thực nghiệm của van den Berg & Ellington (1997) bằng kỹ thuật khói dòng và Dickinson et al. (1999), Sane & Dickinson (2001, 2002) trên mô hình cánh cơ khí rô-bốt "Robofly" đã định danh 3 hiện tượng vật lý cốt lõi: sự gắn kết bền vững của xoáy mép trước (LEV) nhờ dòng chảy hướng trục ổn định, lực tuần hoàn xoay sớm (hiệu ứng Kramer), và cơ chế tái hấp thụ năng lượng vết xoáy (wake capture). Sun & Tang (2004) bằng các mô phỏng CFD 3D Navier-Stokes đã chứng minh định lượng rằng sự đảo pha xoay cánh sớm khoảng $5%T$ tạo ra xung lực nâng đỉnh cực đại ở đầu mỗi nửa chu kỳ vẫy.
  3. Giai đoạn mô phỏng đàn hồi khí động tuyến tính (2005 - 2015): Khi nhận thức rõ cánh côn trùng không phải vật thể cứng tuyệt đối, hàng loạt tác giả đã ứng dụng lý thuyết dầm dẻo Euler-Bernoulli và phương pháp tổng hợp mode dao động (Nguyen & Han, 2010; Roccia et al., 2013; Schwab et al., 2014). Tuy nhiên, các mô hình này mắc phải điểm nghẽn nghiêm trọng: giả thiết dao động điều hòa biên độ nhỏ tuyến tính khiến ma trận độ cứng không cập nhật được sự thay đổi hình học phi tuyến khi cánh uốn cong lớn, đồng thời các mô hình xoáy UVLM cổ điển bị phân kỳ số trị sau 1-2 chu kỳ vẫy do hiện tượng điểm kỳ dị khi các nút mạng xoáy tiến sát bề mặt cánh biến dạng.
  4. Giai đoạn FSI phi tuyến đầy đủ và mô hình hóa thu gọn (2015 - nay): Nhóm nghiên cứu của Nakata & Liu (2012), Dai et al. (2012) và Masoud & Alexeev (2010) thiết lập chuẩn mực độ chính xác cao khi ghép nối bộ giải CFD lưới trượt 3D với CSD phi tuyến. Dẫu vậy, chi phí tính toán khổng lồ (hàng trăm CPU-hours cho 1 ca tính) ngăn cản việc tối ưu hóa thiết kế. Điều này thúc đẩy sự xuất hiện của cách tiếp cận MBD hạt mềm hoặc hệ vật - lò xo (Vanella et al., 2009; Qi et al., 2019; Truong et al., 2021).
+-----------------------------------------------------------------------------------------+
|                                    LITERATURE MAP                                       |
+-----------------------------------------------------------------------------------------+
|  BET & Quasi-Steady (Ellington, 1984; Ansari, 2006)                                     |
|  --> Đơn giản, tính nhanh, nhưng BỎ QUA tính không dừng 3D & LEV                         |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  High-Fidelity CFD-CSD (Nakata & Liu, 2012; Shahzad, 2016)                              |
|  --> Độ chính xác cao, FSI phi tuyến, nhưng CHI PHÍ TÍNH TOÁN CỰC LỚN (24-93h/chu kỳ)    |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  Linear Modal Synthesis + UVLM (Roccia et al., 2013; Schwab et al., 2014)               |
|  --> Tính toán nhanh, nhưng BỊ GIỚI HẠN bởi giả thiết biến dạng nhỏ, mất ổn định vết xoáy|
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  POSITIONING CỦA LUẬN ÁN LÊ VŨ ĐAN THANH (2023):                                       |
|  --> Tích hợp MBD phi tuyến (Hệ vật - lò xo tương đương) + 3D UVLM cải tiến             |
|      (Polhamus LEV Suction + Lamb-Oseen Core Diffusion + Quasi-Newton FSI Solver)       |
|  ==> Giữ vững bản chất biến dạng lớn phi tuyến 3D, kiểm soát xoáy ổn định đa chu kỳ,   |
|      rút ngắn thời gian tính toán xuống hàng chục phút, mở đường cho tối ưu hóa tham số. |
+-----------------------------------------------------------------------------------------+

Trong bức tranh tổng quan đó, tồn tại hai luồng tranh luận học thuật sâu sắc:

  • Tranh luận 1: Nguồn gốc chủ đạo điều khiển biến dạng cánh. Daniel & Combes (2002) dựa trên các thử nghiệm rung cánh trong môi trường không khí và buồng chân không/khí Helium khẳng định rằng lực quán tính kết cấu (inertial-elastic forces) chiếm đến $90%$ nguyên nhân gây biến dạng uốn, tải khí động chỉ đóng vai trò thứ yếu. Trái lại, Norris (2000) và Mountcastle & Daniel (2016) chứng minh rằng ở dải vận tốc quét cao và biên độ góc xoay lớn, áp lực khí động biến thiên cục bộ dọc sải cánh làm thay đổi cục diện phân bố mô-men uốn và xoắn xoay chiều, đặc biệt là tại thời điểm cánh đạt vận tốc quét cực đại giữa chu kỳ ($t/T = 0.25$ và $0.75$).
  • Tranh luận 2: Bản chất của góc lên - xuống (deviation angle $\theta$). Li et al. (2012) cho rằng chuyển động lên - xuống chỉ là biến dạng thụ động không mong muốn do khớp cánh lỏng lẻo. Ngược lại, Luo et al. (2015) và Phillips & Knowles (2011) chứng minh sự tồn tại của "hiệu ứng xoay bổ sung" (Added-rotation effect): khi $\theta \neq 0$, vận tốc góc quét $\vec{\omega}_z$ sẽ chiếu một thành phần vận tốc góc $\vec{\omega}_y$ dọc theo sải cánh, làm thay đổi góc tấn hiệu dụng và gia tăng lực nâng tức thời lên tới $7.3%$.

Vị thế khoa học (positioning) của luận án: Luận án định vị chính xác tại giao điểm giữa tính chính xác vật lý phi tuyến của CFD-CSD và tốc độ xử lý vượt trội của các mô hình MBD thu gọn. So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  1. So với mô hình CFD-CSD của Nakata & Liu (2012) trên loài bướm Manduca sexta, mô hình MBD-UVLM của luận án tái tạo dạng biến dạng uốn - xoắn và hệ số lực nâng với sai số tương đối dưới $8.5%$, trong khi giảm thời gian giải tích xuống hơn 120 lần.
  2. So với mô hình khí động dầm đàn hồi của Roccia et al. (2013), luận án khắc phục triệt để hiện tượng mất ổn định số trị qua nhiều chu kỳ vẫy nhờ giải thuật khuếch tán bán kính lõi xoáy $r_c(t) = \sqrt{4\nu\delta t + r_0^2}$, cho phép mô phỏng hội tụ ổn định liên tục trên 10 chu kỳ vẫy liên tiếp ở chế độ bay treo.

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

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

Luận án tạo ra những bước tiến lý thuyết đáng kể trong lĩnh vực Cơ học vật rắn và Đàn hồi khí động học phỏng sinh học:

  1. Mở rộng lý thuyết động lực học hệ nhiều vật (Kane & Lagrange Formalism) cho cơ hệ sinh học siêu nhẹ biến dạng lớn: Luận án đã hình thức hóa phương pháp chuyển đổi liên tục từ các trường thông số liên tục của cánh sinh học (phân bố khối lượng $m(r)$, mô-men quán tính $I(r)$, độ cứng chống uốn $EI(r)$, độ cứng chống xoắn $GJ(r)$) sang các thông số rời rạc của hệ $N$ vật rắn ($k = 1 \dots N$) liên kết thông qua các khớp cầu và lò xo xoắn đàn hồi 2 bậc tự do uốn - xoắn độc lập.
  2. Hoàn thiện khung lý thuyết tương tác đàn hồi khí động phi tuyến bậc thấp: Chứng minh giải tích và số trị mối quan hệ ghép nối giữa biến dạng góc uốn $\theta_k$ và góc xoắn $\alpha_k$ dọc sải cánh với ma trận biến đổi tọa độ quay Euler theo chuỗi thứ tự trục quay $(3-1-2)$: $\mathbf{B}_{k0} = \mathbf{B}_k(\alpha_k) \mathbf{B}_k(\theta_k) \mathbf{B}_k(\phi_k)$.
  3. Định lập các mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề P1 (Đồng điều kiện cộng hưởng khí động đàn hồi): Hiệu suất tạo lực nâng cực đại của cánh vẫy mềm đạt được khi và chỉ khi tỉ lệ tần số riêng $f_1^* = f_1/f$ được điều chỉnh trùng khớp với điểm dịch pha biến dạng tối ưu, tại đó biến dạng xoắn thụ động đạt cực đại vào đúng thời điểm vận tốc quét tịnh tiến đạt cực trị.
    • Mệnh đề P2 (Cơ chế suy giảm lực nâng do góc lên - xuống lớn): Sự gia tăng lực nâng nhờ hiệu ứng xoay bổ sung sinh ra bởi góc lệch mặt phẳng vẫy $\theta$ tỷ lệ thuận với $\sin\theta$, nhưng bán kính quét hiệu dụng của cánh lại suy giảm tỷ lệ với $R\cos\theta$. Tồn tại một góc tới hạn $\theta_{\text{crit}} \approx 8.6^\circ - 10.0^\circ$, vượt qua ngưỡng này, suy giảm bán kính hiệu dụng sẽ triệt tiêu hoàn toàn lợi ích khí động học, biến chuyển động lên - xuống thành thành phần tiêu tán năng lượng thuần túy.
    • Mệnh đề P3 (Bảo toàn xoáy mép trước nhờ độ cứng biến thiên): Dạng phân bố độ cứng uốn giảm phi tuyến bậc hai $EI(r) = EI_0(1 - a_{1b}\bar{r} - a_{2b}\bar{r}^2)$ có khả năng triệt tiêu xung áp lực cục bộ tại $75%R$, ngăn ngừa hiện tượng vỡ xoáy sớm (vortex bursting) và duy trì sự ổn định của LEV dọc suốt nửa chu kỳ vẫy.
+-----------------------------------------------------------------------------------------+
|                                THEORETICAL FRAMEWORK                                    |
+-----------------------------------------------------------------------------------------+
|  [ĐỘNG HỌC GỐC CÁNH (3-DOF)]                                                           |
|  - Góc quét phi(t), Góc lên-xuống theta(t), Góc xoay alpha(t)                           |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  [BỘ GIẢI ĐỘNG LỰC HỌC KẾT CẤU MBD (N-BODY SPRING MODEL)]                               |
|  - Tọa độ suy rộng eta_j (uốn theta_k, xoắn alpha_k)                                    |
|  - Ma trận xoay Euler B_kR, Phương trình Lagrange loại II:                              |
|    M(eta) * d2(eta)/dt2 + C(eta, d(eta)/dt) + K * eta = Q_aero + Q_spring               |
+-----------------------------------------------------------------------------------------+
                      |                                       ^
       Vị trí, Vận tốc, |                                       | Tải khí động suy rộng
       Biến dạng tức thời|                                       | Q_j^F, Q_j^M
                      v                                       |
+-----------------------------------------------------------------------------------------+
|  [BỘ GIẢI KHÍ ĐỘNG LỰC HỌC 3D EXTENDED UVLM]                                            |
|  - Mạng lưới khung xoáy Gamma_i trên mặt cánh & vết xoáy tự do (Free wake)             |
|  - Mô hình lực hút mép trước Polhamus (LEV model)                                       |
|  - Mô hình khuếch tán lõi xoáy Lamb-Oseen: r_c(t) = sqrt(4*nu*delta_t + r_0^2)          |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  [THUẬT TOÁN GHÉP NỐI 2 CHIỀU QUASI-NEWTON FSI]                                         |
|  - Hội tụ sai số lực & vị trí: ||R(eta)|| < 10^-5                                       |
|  ==> ĐẦU RA: Lực nâng L_m, Công suất P_m, Suất tiêu hao năng lượng eta_p, eta_z        |
+-----------------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liền mạch ba trụ cột lý thuyết:

  1. Lý thuyết Động lực học hệ nhiều vật (MBD): Sử dụng phương trình vi phân chuyển động Lagrange dạng ma trận cho hệ nhiều vật có liên kết đàn hồi: $$\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{K}\mathbf{q} = \mathbf{Q}^{\text{aero}} + \mathbf{Q}^{\text{spring}}$$ trong đó $\mathbf{q} = [\theta_1, \alpha_1, \theta_2, \alpha_2, \dots, \theta_N, \alpha_N]^T$ là véc-tơ tọa độ suy rộng mô tả chuyển vị góc uốn và góc xoắn của từng đoạn vật thể thứ $k$.
  2. Lý thuyết Độ cứng dầm đàn hồi phi tuyến tương đương: Hàm mục tiêu tối ưu hóa chi phí $J$ được thiết lập dựa trên sai số dạng dao động riêng và giá trị tần số dao động riêng giữa mô hình phần tử hữu hạn chi tiết (ANSYS APDL với phần tử dầm BEAM188 cho gân và phần tử màng SHELL181 cho màng cánh) và mô hình dầm tương đương: $$J = \sum_{i=1}^{M} w_i \left( \frac{f_{\text{FEM}}^{(i)} - f_{\text{beam}}^{(i)}}{f_{\text{FEM}}^{(i)}} \right)^2 + \sum_{i=1}^{M} \gamma_i (1 - \text{MAC}i)$$ với $\text{MAC}i = \frac{(\mathbf{\phi}{\text{FEM}}^{(i)T} \mathbf{\phi}{\text{beam}}^{(i)})^2}{(\mathbf{\phi}{\text{FEM}}^{(i)T} \mathbf{\phi}{\text{FEM}}^{(i)})(\mathbf{\phi}{\text{beam}}^{(i)T} \mathbf{\phi}{\text{beam}}^{(i)})}$ là tiêu chuẩn tương đồng mode dao động (Modal Assurance Criterion).
  3. Lý thuyết Khí động không dừng UVLM mở rộng: Trường vận tốc cảm ứng được tính toán thông qua định luật Biot-Savart giải tích trên từng đoạn xoáy biên và xoáy tự do: $$\vec{V}(\vec{r}) = \frac{\Gamma}{4\pi} \frac{\vec{r}_1 \times \vec{r}_2}{|\vec{r}_1 \times \vec{r}2|^2} \left[ \vec{r}0 \cdot \left( \frac{\vec{r}1}{r_1} - \frac{\vec{r}2}{r_2} \right) \right] \frac{r^2}{r^2 + r_c^2}$$ Lực nâng bổ sung từ xoáy mép trước được lượng hóa thông qua phương pháp tương tự lực hút mép trước của Polhamus: $$\vec{F}{\text{LEV}} = \rho \Gamma{\text{LEV}} (\vec{V}{\text{rel}} \times \vec{t}{\text{LE}})$$ Tải khí động suy rộng $Q_j^F, Q_j^M$ tác dụng lên từng đốt vật cứng được tích phân trực tiếp từ trường áp suất chênh lệch qua bề mặt màng cánh tương ứng.

Điều kiện biên xác định rõ ràng: Mô hình áp dụng cho chế độ bay treo trong không khí tĩnh ($\vec{V}_\infty = 0$), khối lượng riêng không khí $\rho = 1.225\text{ kg/m}^3$, độ nhớt động học $\nu = 1.5 \times 10^{-5}\text{ m}^2/\text{s}$. Cánh được ngàm cứng tại khớp gốc cánh ($r = 0$) và tự do tại mút cánh ($r = R$).

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng thực nghiệm kết hợp cơ học tính toán tất định (deterministic computational mechanics paradigm). Toàn bộ quy trình giải thuật mô phỏng được phân tầng qua ba cấp độ thiết kế:

  • Cấp độ 1 (Micro-Structural Modeling): Xây dựng hình học cánh bướm Manduca sexta chi tiết với 18 nhánh gân chính và phụ dựa trên số liệu giải phẫu sinh học của Combes & Daniel (2003). Mô hình FEM trong ANSYS APDL gồm 1,428 phần tử dầm không gian và 3,856 phần tử màng mỏng đàn hồi tuyến tính trực giao, mô-đun đàn hồi dọc gân $E_{\text{vein}} = 2.5 - 4.5\text{ GPa}$, mô-đun đàn hồi màng $E_{\text{mem}} = 1.0 - 2.0\text{ MPa}$, khối lượng riêng $\rho_{\text{vein}} = 1,150\text{ kg/m}^3$.
  • Cấp độ 2 (Equivalent Beam & Parameter Identification): Trích xuất các phân bố đặc tính quán tính và độ cứng dọc trục đàn hồi không thứ nguyên $\bar{r} = r/R \in [0, 1]$. Bằng phương pháp bình phương tối thiểu có ràng buộc, hàm phân bố khối lượng đơn vị $\bar{m}(\bar{r})$ và mô-men quán tính đơn vị $\bar{I}(\bar{r})$ được xấp xỉ bằng đa thức bậc 4. Hàm độ cứng uốn $EI(r)$ và xoắn $GJ(r)$ được tối ưu hóa đạt độ trùng khớp tần số dao động riêng hai mode đầu tiên ($f_1$ uốn thuần túy, $f_2$ xoắn thuần túy) với chỉ số $\text{MAC} > 0.96$.
  • Cấp độ 3 (Multi-Body Discretization & FSI Solver): Rời rạc hóa dầm thành $N = 5 - 8$ vật cứng liên tiếp. Mỗi vật thứ $k$ có khối lượng $m_k$, ten-xơ quán tính $\mathbf{J}k$, nối với vật $k-1$ bằng cụm lò xo xoắn có độ cứng uốn $K{bk}$ và độ cứng xoắn $K_{tk}$. Chương trình tính toán FSI ghép nối bộ giải MBD với bộ giải UVLM được lập trình trên nền tảng MATLAB và liên kết kiểm chứng chéo với phần mềm thương mại đa vật thể chuyên dụng MSC/ADAMS.
+-----------------------------------------------------------------------------------------+
|                              RESEARCH DESIGN ARCHITECTURE                               |
+-----------------------------------------------------------------------------------------+
|  [BƯỚC 1: DỮ LIỆU SINH HỌC & MÔ HÌNH PHẦN TỬ HỮU HẠN CHI TIẾT (ANSYS APDL)]            |
|  - Cánh bướm Manduca sexta: 18 nhánh gân, màng mỏng Mylar/Chitin                        |
|  - Trích xuất: Dạng dao động riêng f_1 (uốn), f_2 (xoắn), phân bố m(r), I(r)           |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  [BƯỚC 2: TỐI ƯU HÓA THÔNG SỐ DẦM TƯƠNG ĐƯƠNG & RỜI RẠC HÓA HỆ NHIỀU VẬT (MBD)]       |
|  - Tối ưu hóa hàm chi phí J: Đồng pha tần số riêng và Mode shape (MAC > 0.96)          |
|  - Rời rạc hóa thành N = 5 đến 8 vật cứng nối bởi các lò xo uốn K_bk và xoắn K_tk       |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  [BƯỚC 3: MÔ HÌNH KHÍ ĐỘNG UVLM MỞ RỘNG VỚI XOÁY MÉP TRƯỚC (LEV)]                       |
|  - Mạng lưới khung xoáy 3D biến dạng động theo vị trí các vật cứng                     |
|  - Khử kỳ dị bằng bán kính lõi xoáy Lamb-Oseen r_c(t)                                   |
|  - Tính lực hút mép trước Polhamus và cuộn xoáy mép sau                                |
+-----------------------------------------------------------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------------+
|  [BƯỚC 4: BỘ GIẢI GHÉP NỐI FSI HAI CHIỀU BẰNG THUẬT TOÁN QUASI-NEWTON]                 |
|  - Cập nhật biến dạng kết cấu <---> Cập nhật trường áp suất khí động học                |
|  - Kiểm chứng chéo với MSC/ADAMS, dữ liệu thực nghiệm PIV của Lua et al. & Combes-Daniel|
+-----------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu và tam giác đạc (Triangulation) phương pháp được thực thi nghiêm ngặt qua 4 bước kiểm chứng độc lập:

  1. Kiểm chứng bộ giải động lực học MBD: Thiết lập bài toán con hệ 3 con lắc ngược phi tuyến chịu kích động điều hòa. So sánh quỹ đạo chuyển vị không gian của điểm kiểm tra $A$ giữa thuật toán tự phát triển và mô hình ADAMS, độ lệch chuẩn $RMSE < 1.2 \times 10^{-4}\text{ m}$.
  2. Kiểm chứng phân tích dao động đàn hồi cánh: So sánh tần số dao động riêng thứ nhất ($f_1$) và thứ hai ($f_2$) của cánh trước bướm Manduca sexta:
    • Mô hình FEM chi tiết: $f_{1\text{FEM}} = 27.2\text{ Hz}$, $f_{2\text{FEM}} = 76.5\text{ Hz}$.
    • Mô hình dầm tối ưu: $f_{1\text{beam}} = 27.0\text{ Hz}$ (sai số $0.74%$), $f_{2\text{beam}} = 75.8\text{ Hz}$ (sai số $0.92%$).
    • Số chỉ $\text{MAC}_1 = 0.988$, $\text{MAC}_2 = 0.964$, khẳng định tính tương đương động lực học tuyệt đối của mô hình dầm.
  3. Kiểm chứng mô hình khí động UVLM mở rộng: So sánh hệ số lực nâng tức thời $C_L(t)$ với kết quả đo thực nghiệm trên cánh bướm của Lua et al. (2014) trong kênh dẫn lưu chất. Đồ thị $C_L$ ghi nhận chính xác 2 đỉnh tạo lực nâng tại thời điểm đầu và giữa nửa chu kỳ vẫy với hệ số tương quan tuyến tính $R^2 = 0.942$.
  4. Kiểm chứng tích hợp FSI hai chiều: Đối chiếu biên độ góc uốn $\theta(t)$ tại các điểm đo dọc sải cánh ($\bar{r} = 0.35, 0.55, 0.75, 0.95$) với số liệu thực nghiệm quay video tốc độ cao của Combes & Daniel (2003). Góc uốn cực đại tại mút cánh $\theta_{\text{max}} = 14.8^\circ$ từ mô phỏng hoàn toàn trùng khớp với giá trị đo đạc thực nghiệm $15.1^\circ \pm 0.6^\circ$.
+-----------------------------------------------------------------------------------------+
|                                    VALIDATION MATRIX                                    |
+-----------------------------------------------------------------------------------------+
| Cấp độ kiểm chứng         | Phương pháp so sánh        | Tiêu chí định lượng | Kết quả  |
+---------------------------+----------------------------+---------------------+----------+
| 1. Động lực học MBD       | ADAMS Benchmark            | Sai số RMSE         | < 0.12 mm|
| 2. Dao động riêng kết cấu | FEM ANSYS (Beam + Shell)   | Sai số tần số / MAC | < 1% /   |
|                           |                            |                     | > 0.96   |
| 3. Khí động học UVLM      | Thực nghiệm PIV Lua (2014) | Tương quan R^2      | 0.942    |
| 4. Toàn diện FSI 2 chiều  | Thực nghiệm Combes (2003)  | Biên độ uốn mút     | 14.8° so |
|                           |                            |                     | vs 15.1° |
+-----------------------------------------------------------------------------------------+

Dữ liệu và phân tích tham số

Không gian tham số khảo sát toàn diện của luận án bao gồm:

  • Tần số vẫy kích động: $f \in [15, 35]\text{ Hz}$ (bước quét $\Delta f = 1.0\text{ Hz}$).
  • Biên độ góc quét: $\phi_a \in [45^\circ, 75^\circ]$ (bước quét $5^\circ$).
  • Biên độ góc xoay: $\alpha_a \in [30^\circ, 60^\circ]$ (bước quét $5^\circ$).
  • Góc lên - xuống trung bình: $\theta_0 \in [-10^\circ, 10^\circ]$ và biên độ $\theta_a \in [0^\circ, 20^\circ]$.
  • Tỉ lệ tần số dao động riêng không thứ nguyên: $f_1^* = f_1/f \in [1.0, 5.0]$.
  • Các hệ số dạng phân bố độ cứng uốn và xoắn: $a_{1b}, a_{2b}, a_{1t}, a_{2t} \in [-0.8, 0.8]$.

Các đại lượng đánh giá hiệu quả đầu ra được tính toán trung bình theo chu kỳ vẫy ổn định ($t \ge 3T$):

  • Lực nâng trung bình: $L_m = \frac{1}{T}\int_0^T L(t) dt \text{ [mN]}$.
  • Công suất khí động tiêu thụ trung bình: $P_m = \frac{1}{T}\int_0^T \left[ Q_\phi(t)\dot{\phi}(t) + Q_\theta(t)\dot{\theta}(t) + Q_\alpha(t)\dot{\alpha}(t) \right] dt \text{ [mW]}$.
  • Suất tiêu thụ công suất trên một đơn vị lực nâng (Power Economy): $\eta_p = P_p / L_m$ (trường hợp có dự trữ năng lượng đàn hồi) và $\eta_z = P_z / L_m$ (trường hợp tiêu tán hoàn toàn năng lượng âm, không dự trữ) $[\text{W/N}]$.

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

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

  1. Xác lập vùng cộng hưởng khí động đàn hồi tối ưu ($f_1^ = 2.0 - 3.0$):* Khi tỉ lệ $f_1^* < 1.8$ (cánh quá mềm), cánh bị uốn cong quá mức ($\theta_{\text{tip}} > 28^\circ$), làm suy giảm diện tích chiếu hiệu dụng và gây hiện tượng bóc tách xoáy mép trước sớm, dẫn đến lực nâng trung bình $L_m$ suy giảm hơn $35%$. Ngược lại, khi $f_1^* > 4.0$ (cánh quá cứng tiệm cận cánh cứng tuyệt đối), cánh không tận dụng được biến dạng xoắn thụ động để tạo góc tấn tối ưu ở giữa hành trình. Lực nâng cực đại và hiệu suất năng lượng $\eta_p$ đạt giá trị tối ưu toàn cục khi $f_1^* \approx 2.4 - 2.6$, tại đó công suất tiêu thụ giảm $14.2%$ so với cánh cứng trong khi lực nâng trung bình duy trì tương đương.
Lực nâng & Công suất
  ^
  |               /---\  <-- Vùng tối ưu khí động đàn hồi (f1* = 2.0 - 3.0)
  |              /     \     (Lực nâng cao nhất, Tiết kiệm năng lượng 14.2%)
  |             /       \
  |   Cánh mềm /         \------------------  Cánh cứng
  |   (f1* < 1.8)                              (f1* > 4.0)
  |   Lực giảm 35%
  +---------------------------------------------> Tỉ lệ tần số riêng f1* = f1/f
     1.0      1.8     2.4   3.0              5.0

  1. Quy luật gia tăng lực nâng phi tuyến của góc xoay sớm (Phase Advance): Hiện tượng xoay sớm pha $\Delta\psi \approx 5%T$ kích hoạt hiệu ứng bắt vết xoáy (wake capture) cực mạnh: dòng xoáy cuộn từ nửa chu kỳ trước va đập vào mặt dưới của cánh khi cánh vừa đảo chiều ở góc tấn đón gió thuận lợi, tạo ra một xung lực nâng đỉnh đạt giá trị gấp 3.2 lần trọng lượng bướm ($L_{\text{peak}} \approx 3.2 mg$). Nếu xoay chậm pha ($\Delta\psi < 0$), lực nâng trung bình sụt giảm tức thì $42%$ do cánh chuyển động ngược chiều trường vận tốc cảm ứng của vết xoáy cũ.

  2. Lượng hóa ranh giới hiệu dụng của hiệu ứng xoay bổ sung (Added-rotation effect): Luận án làm sáng tỏ một kết quả phản trực giác: việc tăng góc lên - xuống trung bình $\theta_0$ từ $0^\circ$ lên $+8.6^\circ$ giúp tăng lực nâng tức thời thêm $7.1%$ nhờ vận tốc góc cảm ứng $\vec{\omega}_y$ dọc sải cánh hỗ trợ duy trì góc tấn dương. Tuy nhiên, cái giá phải trả về mặt năng lượng là công suất tiêu thụ $P_m$ tăng phi mã tới $28.4%$, khiến chỉ số suất tiêu thụ $\eta_z$ tăng từ $12.1\text{ W/N}$ lên $15.6\text{ W/N}$. Điều này chứng minh côn trùng trong tự nhiên chỉ kích hoạt chuyển động lên - xuống lớn trong các pha cơ động khẩn cấp né tránh vật cản hoặc tấn công, thay vì duy trì thường trực trong chế độ bay treo tiết kiệm năng lượng.

  3. Ưu thế khí động học của dạng phân bố độ cứng phi tuyến: Khảo sát chi tiết ảnh hưởng của hệ số dạng $a_{1b}, a_{2b}$ cho thấy dạng phân bố độ cứng uốn giảm nhanh ở $50%$ sải cánh ngoài mút ($a_{1b} > 0, a_{2b} > 0$) giúp mút cánh uốn cong mượt mà theo dạng parabol, phân phối lại áp suất khí động dọc sải cánh đồng đều hơn, ngăn chặn hiện tượng tập trung ứng suất và giảm tiêu hao năng lượng công suất quán tính $\sim 11.5%$.

+-----------------------------------------------------------------------------------------+
|                                    KEY FINDINGS SUMMARY                                 |
+-----------------------------------------------------------------------------------------+
| Hiện tượng vật lý         | Bằng chứng định lượng cụ thể     | Ý nghĩa cơ học           |
+---------------------------+----------------------------------+--------------------------+
| 1. Tỉ lệ tần số f1*       | Tối ưu tại f1* = 2.4 - 2.6       | Tái hấp thu năng lượng   |
|                           | Giảm 14.2% công suất tiêu thụ    | biến dạng, ổn định LEV   |
+---------------------------+----------------------------------+--------------------------+
| 2. Sớm pha xoay           | Sớm pha 5% chu kỳ vẫy            | Kích hoạt wake capture,  |
|                           | Xung lực nâng gấp 3.2 lần mg     | tăng lực nâng đỉnh       |
+---------------------------+----------------------------------+--------------------------+
| 3. Góc lên - xuống theta  | Lực nâng tăng cực đại ở theta=8.6| Hiệu ứng xoay bổ sung    |
|                           | Công suất tiêu hao tăng vọt 28.4%| chỉ phù hợp khi cơ động  |
+---------------------------+----------------------------------+--------------------------+
| 4. Phân bố độ cứng giảm   | Hệ số phân bố a_1b, a_2b > 0     | Phân bố đều trường áp,   |
|                           | Giảm 11.5% công suất quán tính   | triệt tiêu đỉnh ứng suất |
+---------------------------+----------------------------------+--------------------------+

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp một phương pháp luận chuẩn xác để nghiên cứu các bài toán tương tác chất lưu - kết cấu phi tuyến lớn có số bậc tự do thu gọn, có thể chuyển giao mở rộng cho các hệ thống kỹ thuật khác như cánh quạt trực thăng mềm, cánh tuabin gió nổi ngoài khơi, và vây bơi rô-bốt cá phỏng sinh học dưới nước.
  • Về mặt kỹ thuật chế tạo FWMAV: Đưa ra hướng dẫn thiết kế định lượng trực tiếp cho các kỹ sư: không nên cố gắng chế tạo cánh cứng tuyệt đối bằng cách gia cố quá nhiều thanh carbon làm tăng khối lượng quán tính. Thay vào đó, cần điều chỉnh mật độ gân và chiều dày màng mỏng polymer (Mylar/Kapton) sao cho tần số dao động riêng thứ nhất của cánh gấp từ $2.0$ đến $2.5$ lần tần số phát động của động cơ ($f_1 \approx 2.0 - 2.5 f_{\text{motor}}$).
  • Về mặt an ninh - quốc phòng: Các phát hiện về kiểm soát pha xoay và hiệu ứng xoay bổ sung đặt nền tảng cho việc lập trình thuật toán điều khiển tự động cho các vi thiết bị bay trinh sát ngụy trang kiểu côn trùng (cỡ $< 15\text{ cm}$, khối lượng $< 20\text{ g}$), nâng cao tính cơ động trong không gian hẹp (hang động, đô thị, công trình kín) và kéo dài thời gian hoạt động của pin nguồn.

Limitations và Future Research

Nhằm đảm bảo tính khách quan khoa học, luận án chỉ rõ 4 giới hạn tồn tại:

  1. Giới hạn giả thiết hình học mặt phẳng dây cung: Mô hình dầm MBD hiện thời coi các dải cánh là phân đoạn cứng phẳng theo phương dây cung (chordwise), bỏ qua biến dạng cong dạng lòng máng (chordwise camber deformation, xấp xỉ $2 - 5%$ chiều dài dây cung) và hiện tượng nhăn màng (membrane wrinkling) cục bộ khi chịu lực cắt lớn.
  2. Giới hạn mô hình dòng thế UVLM: Mặc dù đã tích hợp mô hình khuếch tán xoáy và lực hút Polhamus, phương pháp UVLM vẫn dựa trên nền tảng dòng không nén phi nhớt, chưa mô tả được đầy đủ lớp biên nhớt 3D sát bề mặt cánh và hiện tượng tái bám dòng phức tạp ở các số Reynolds cực thấp ($Re < 500$).
  3. Giới hạn điều kiện biên thân cố định: Luận án mới chỉ giải quyết bài toán cân bằng lực nâng trong chế độ bay treo độc lập của cánh, chưa ghép nối 6 bậc tự do chuyển động thân thể (6-DOF body flight dynamics) để khảo sát bài toán ổn định động học dài hạn và dao động ghép thân - cánh (pitch-stroke body oscillation).
  4. Giới hạn chế độ bay: Chưa mở rộng sang chế độ bay tiến (forward flight) với vận tốc dòng tự do $V_\infty > 0$ và góc nghiêng mặt phẳng vẫy biến thiên từ $0^\circ$ đến $90^\circ$.

Chương trình nghiên cứu phát triển giai đoạn tiếp theo (5 - 10 năm tới) định hình 5 hướng đi:

  • Hướng 1: Nâng cấp mô hình mạng nhiều vật 2 chiều (2D multi-body lattice grid) để mô phỏng đồng thời biến dạng uốn - xoắn dọc sải cánh lẫn biến dạng cong lòng máng dọc dây cung.
  • Hướng 2: Tích hợp bộ giải FSI MBD-UVLM vào vòng lặp điều khiển phản hồi kín 6-DOF mô phỏng đầy đủ động lực học bay tự do của cả thân và cánh côn trùng dưới tác động của gió giật (gust perturbation).
  • Hướng 3: Ứng dụng thuật toán tối ưu hóa tiến hóa đa mục tiêu (NSGA-II) kết hợp mạng nơ-ron học sâu (Deep Surrogate Models) để tự động hóa quá trình tối ưu topo cấu trúc gân composite cho cánh FWMAV.
  • Hướng 4: Tiến hành thử nghiệm trong ống khí động đo trường vận tốc PIV đồng bộ cảm biến lực đa thành phần trên mô hình cánh cơ điện tử (MEMS-flapper) để hiệu chỉnh các hệ số thực nghiệm khuếch tán xoáy.
  • Hướng 5: Nghiên cứu vật liệu thông minh (Piezoelectric / Shape Memory Alloys) tích hợp trên khung gân cánh để chủ động điều khiển biến dạng đàn hồi thích ứng theo thời gian thực (Active morphing aeroelastic control).

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

  • Tác động học thuật (Academic Impact): Công trình đóng góp một phương pháp tính toán đột phá cho cộng đồng Cơ học tính toán và Robot phỏng sinh học quốc tế. Với mã nguồn giải thuật MBD-UVLM hiệu năng cao, nghiên cứu mở ra tiền đề trích dẫn mạnh mẽ trong các tạp chí chuyên ngành hàng đầu như Journal of Fluid Mechanics, Bioinspiration & Biomimetics, AIAA Journal, và IEEE Transactions on Robotics.
  • Chuyển đổi công nghệ công nghiệp & Quốc phòng (Industry & Defense Transformation): Kết quả nghiên cứu cung cấp cơ sở dữ liệu kỹ thuật then chốt cho các viện nghiên cứu thuộc Bộ Quốc phòng và các tập đoàn công nghiệp công nghệ cao (Viettel, Z113) trong việc làm chủ quy trình thiết kế, mô phỏng và chế tạo các dòng micro-drone phỏng sinh học phục vụ trinh sát chiến thuật ngầm, thám sát địa hình độc hại và cứu hộ thảm họa sập đổ công trình.
  • Lợi ích xã hội và đào tạo (Societal & Educational Value): Thúc đẩy ứng dụng công nghệ vi cơ điện tử (MEMS) và vật liệu composite tiên tiến tại Việt Nam; cung cấp tài liệu tham khảo chuyên sâu, giáo trình chuẩn mực cho chương trình đào tạo sau đại học ngành Cơ học kỹ thuật, Kỹ thuật Hàng không Vũ trụ và Cơ điện tử.

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

+-----------------------------------------------------------------------------------------+
|                                    STAKEHOLDER BENEFITS                                 |
+-----------------------------------------------------------------------------------------+
| Nhóm đối tượng            | Giá trị & Lợi ích cụ thể nhận được                          |
+---------------------------+-------------------------------------------------------------+
| 1. Nghiên cứu sinh &      | Khung phương pháp luận FSI MBD thu gọn hoàn chỉnh;          |
|    Học giả trẻ            | Công cụ mô phỏng giải quyết bài toán biến dạng lớn phi tuyến|
+---------------------------+-------------------------------------------------------------+
| 2. Nhà khoa học Cơ học    | Làm sáng tỏ các cơ chế vật lý (LEV, Wake capture,           |
|    & Khí động lực học     | Added-rotation) dưới tương tác đàn hồi khí động học 3D      |
+---------------------------+-------------------------------------------------------------+
| 3. Kỹ sư R&D thiết bị     | Công thức & hướng dẫn tối ưu độ cứng cánh (f1* = 2.0 - 3.0) |
|    bay không người lái    | Rút ngắn 90% chu kỳ mô phỏng thử nghiệm khung cánh composite|
+---------------------------+-------------------------------------------------------------+
| 4. Cơ quan quản lý &      | Cơ sở khoa học định hình chiến lược phát triển khí tài trinh|
|    Bộ Quốc phòng          | sát vi mô tự chủ công nghệ cao, nâng cao năng lực tác chiến |
+---------------------------+-------------------------------------------------------------+

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

1. Đóng góp lý thuyết độc đáo và đột phá nhất của luận án là gì?

Đó là việc thiết lập và giải quyết thành công mô hình tương tác đàn hồi khí động học FSI phi tuyến 3D bằng cách tiếp cận cơ hệ nhiều vật (MBD) kết hợp phương pháp xoáy không dừng mở rộng (Extended UVLM). Luận án đã mở rộng lý thuyết dầm dẻo phi tuyến tương đương, bắc cầu nối giải tích giữa cấu trúc gân - màng vi mô phức tạp với hệ các vật cứng lò xo xoắn 2 bậc tự do thu gọn. Đóng góp này phá vỡ thế bế tắc kéo dài nhiều thập kỷ giữa phương pháp giải tích tuyến tính biên độ nhỏ (quá đơn giản, sai lệch vật lý) và phương pháp CFD-CSD đầy đủ (quá nặng nề, không thể quét tham số tối ưu).

2. Sự cách tân trong phương pháp nghiên cứu so sánh với các công trình quốc tế tiêu biểu?

So với nghiên cứu của Roccia et al. (2013) vốn bị phân kỳ tính toán sau 1-2 chu kỳ do kỳ dị mạng xoáy, luận án đã tích hợp thuật toán khuếch tán bán kính lõi xoáy Lamb-Oseen và mô hình lực hút mép trước Polhamus, duy trì tính ổn định số trị tuyệt đối qua hàng chục chu kỳ vẫy. So với mô hình CFD biên nhúng của Shahzad et al. (2016) tiêu tốn 93 giờ cho 1 chu kỳ, phương pháp của luận án giải quyết 1 chu kỳ vẫy chỉ trong 12-15 phút trên máy tính cá nhân tiêu chuẩn (giảm hơn 300 lần thời gian chạy máy) mà vẫn bảo toàn độ chính xác biến dạng góc uốn và xoắn với sai số dưới $8.5%$.

3. Phát hiện nào gây bất ngờ nhất và có bằng chứng định lượng ra sao?

Phát hiện bất ngờ nhất là vai trò thực sự của chuyển động lên - xuống (góc $\theta$). Trước đây nhiều nghiên cứu cho rằng chuyển động này luôn có lợi cho việc tạo lực nâng. Luận án chỉ ra rằng dù góc $\theta = 8.6^\circ$ tăng lực nâng thêm $7.1%$ nhờ hiệu ứng xoay bổ sung, nhưng nó làm tiêu hao công suất khí động vọt lên $28.4%$, làm suy giảm nghiêm trọng hiệu suất kinh tế năng lượng $\eta_z$. Điều này bác bỏ giả thiết cho rằng côn trùng duy trì chuyển động lên - xuống liên tục để tiết kiệm lực trong bay treo.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ. Luận án công bố chi tiết: hệ phương trình vi phân chuyển động Lagrange dạng ma trận, bảng thông số hình học và cơ lý tính của 18 nhánh gân cánh bướm Manduca sexta, các hệ số đa thức xấp xỉ khối lượng $\bar{m}(\bar{r})$ và quán tính $\bar{I}(\bar{r})$, độ cứng lò xo rời rạc $K_{bk}, K_{tk}$, giải thuật Quasi-Newton bước lặp thời gian $\Delta t = T/100$, cùng điều kiện hội tụ sai số dư $\epsilon \le 10^{-5}$, cho phép bất kỳ nhóm nghiên cứu độc lập nào cũng có thể lập trình tái lập chính xác kết quả.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm tập trung vào 3 cột mốc: (1) 2024-2026: Mở rộng mô hình MBD lưới 2D xét biến dạng cong lòng máng dây cung và hiệu ứng nhăn màng; (2) 2027-2029: Ghép nối đầy đủ 6-DOF động lực học bay tự do thân - cánh trong môi trường gió nhiễu loạn thực; (3) 2030-2033: Tích hợp cơ chế vật liệu thông minh biến đổi hình dạng chủ động (Active Morphing) và trí tuệ nhân tạo biên (Edge AI) điều khiển đàn hồi khí động học thời gian thực trên vi thiết bị bay FWMAV thực tế.

Kết luận

  1. Xây dựng thành công phương pháp luận và chương trình mô phỏng tương tác kết cấu - chất lưu (FSI) 3D tiên tiến cho cánh vẫy kiểu côn trùng bằng cách tiếp cận động lực học hệ nhiều vật (MBD) tích hợp phương pháp xoáy không dừng mở rộng (UVLM).
  2. Thiết lập quy trình tối ưu hóa chuyển đổi tương đương chính xác từ mô hình phần tử hữu hạn cấu trúc gân - màng chi tiết sang mô hình dầm dẻo và hệ vật - lò xo thu gọn, đạt tiêu chuẩn tương đồng mode dao động $\text{MAC} > 0.96$.
  3. Khám phá và lượng hóa chính xác vùng cộng hưởng khí động đàn hồi tối ưu tại tỉ lệ tần số riêng $f_1^* = 2.0 - 3.0$, giúp giảm $14.2%$ công suất tiêu thụ năng lượng so với cánh cứng tuyệt đối.
  4. Làm sáng tỏ cơ chế vật lý của hiệu ứng xoay bổ sung sinh ra từ chuyển động lên - xuống, xác định góc giới hạn tối ưu $\theta_{\text{crit}} \approx 8.6^\circ - 10.0^\circ$ và chứng minh sự đánh đổi lớn về mặt năng lượng của góc lệch mặt phẳng vẫy.
  5. Chứng minh định lượng lợi thế khí động học của dạng phân bố độ cứng uốn và xoắn giảm dần phi tuyến về mút cánh, cung cấp cơ sở khoa học trực tiếp cho việc thiết kế, tối ưu hóa kết cấu cánh composite siêu nhẹ cho các phương tiện bay cánh vẫy siêu nhỏ (FWMAV).
  6. Mở ra một hướng tiếp cận tính toán mới có hiệu năng cao, khả thi cho các bài toán tối ưu hóa đa tham số và thiết kế hệ thống điều khiển tự động cho các thiết bị bay tự hành phỏng sinh học thế hệ mới.