Tổng quan về luận án

Nghiên cứu động lực học và ổn định của các kết cấu tấm vỏ mỏng thông minh là một trong những trọng tâm tiên phong của ngành Cơ học vật rắn và Cơ kỹ thuật hiện đại, đặc biệt trong bối cảnh các phương tiện bay, vũ khí công nghệ cao và công trình hàng không - vũ trụ đòi hỏi khả năng làm việc bền bỉ dưới tác động của tải trọng khí động phức tạp. Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Ngọc Thủy với đề tài "Phân tích động lực học tấm composite áp điện có gân gia cường chịu tải trọng khí động" (Chuyên ngành: Cơ kỹ thuật, Mã ngành: 9.01, Học viện Kỹ thuật Quân sự, 2020) đã giải quyết trọn vẹn bài toán ghép nối đa trường vật lý (điện - cơ học - khí động học phi tuyến) cho kết cấu tấm composite mỏng được tăng cứng bằng hệ gân và tích hợp vật liệu áp điện thông minh (PZT).

                 +---------------------------------------------------+
                 |           DÒNG KHÍ ĐỘNG DƯỚI ÂM (U, α)            |
                 |      Mô hình lực tự kích thích Scanlan (1978)     |
                 +-------------------------+-------------------------+
                                           | Tải trọng khí động
                                           v
+-----------------------------------------------------------------------------------+
|               KẾT CẤU TẤM COMPOSITE CÓ GÂN GIA CƯỜNG (FSDT Mindlin)              |
|                     Ma trận độ cứng hình học phi tuyến Von-Kármán                 |
+------------------------------------------+----------------------------------------+
                                           | Biến dạng uốn - xoắn
                                           v
+-----------------------------------------------------------------------------------+
|                     HỆ CẢM BIẾN - KÍCH HOẠT ÁP ĐIỆN (PZT)                         |
|  Cảm biến (Sensor): {D} = [d]{σ} + [p]{E}  ==>  Kích hoạt (Actuator): {σ} = [e]{E}|
|            Điều khiển hồi tiếp chủ động (Closed-loop Active Control)               |
+-----------------------------------------------------------------------------------+

Trước công trình này, khoảng trống nghiên cứu (research gap) trong cơ học kết cấu composite thông minh thể hiện rõ ở sự thiếu vắng các mô hình số phi tuyến toàn diện có khả năng mô tả đồng thời ba yếu tố: tính chất dị hướng của vật liệu composite nhiều lớp, ứng xử phi tuyến hình học theo Von-Kármán của tấm có gân lệch tâm, và lực khí động tự kích thích dưới âm (subsonic self-excited aerodynamic loads) có kể đến lực cản kết cấu và cản khí động. Phần lớn các công bố quốc tế trước đó (như Maiti et al., 2008; Shen, 2009; Singha et al., 2001) chỉ dừng lại ở bài toán phân tích tĩnh, dao động tự do tuyến tính hoặc mô hình hóa lực khí động theo lý thuyết Piston siêu âm ($M \ge 2$) trên tấm trơn không gân. Luận án đã thiết lập các câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (Q1): Làm thế nào để xây dựng mô hình phần tử hữu hạn (PTHH) 9 nút đẳng tham số mô tả chính xác tương tác ghép nối cơ - điện phi tuyến cho kết cấu tấm composite có gân gia cường lệch tâm chịu tải trọng khí động?
  • Câu hỏi nghiên cứu 2 (Q2): Quy luật ảnh hưởng tương hỗ giữa vận tốc dòng khí $U$, góc tới $\alpha$, điện áp điều khiển $V$, tỷ số hình học tấm $h/W$, chiều cao gân $h_g$ và góc đặt sợi $\theta$ đến đáp ứng động lực học và ranh giới ổn định mất ổn định của tấm diễn ra như thế nào?
  • Câu hỏi nghiên cứu 3 (Q3): Khả năng dập tắt dao động chủ động bằng vật liệu áp điện PZT thông qua mạch hồi tiếp vận tốc và chuyển vị có đạt độ tin cậy và tương thích khi kiểm chứng thực nghiệm trong hầm gió khí động hay không?
  • Giả thuyết nghiên cứu 1 (H1): Việc tích hợp hệ gân gia cường rời rạc kết hợp với các lớp áp điện PZT điều khiển hồi tiếp sẽ làm tăng vận tốc dòng khí tới hạn $U_{cr}$ lên trên 30% so với kết cấu tấm phẳng trơn.
  • Giả thuyết nghiên cứu 2 (H2): Mô hình lực khí động Scanlan kết hợp với lý thuyết biến dạng trượt bậc nhất Mindlin và phi tuyến Von-Kármán sẽ dự báo chính xác hiện tượng mất ổn định động theo tiêu chuẩn Budiansky-Roth với sai số so với thực nghiệm dưới 10%.

Khung lý thuyết của luận án tích hợp Lý thuyết biến dạng trượt bậc nhất của Mindlin (FSDT), quan hệ phi tuyến hình học Von-Kármán, Phương trình cấu động lực học điện môi áp điện Maxwell-Voigt, và Mô hình lực khí động học thực nghiệm Scanlan. Phạm vi nghiên cứu bao quát các khảo sát số đa thông số và chương trình thực nghiệm trên hầm gió chuyên dụng tại dải vận tốc dòng khí $U = 0 \div 30$ m/s, tần số kích thích ngoại lực $f = 6.944$ Hz, tạo nên bước đột phá về khả năng làm chủ công nghệ tính toán kết cấu cánh và thân vỏ phương tiện bay.

Literature Review và Positioning

Phân tích tổng quan y văn quốc tế và trong nước cho thấy sự phát triển của lĩnh vực cơ học kết cấu composite áp điện chịu tải khí động được định hình bởi ba dòng nghiên cứu chính:

                            TIẾN TRÌNH Y VĂN VÀ ĐỊNH VỊ LUẬN ÁN
                            
      Mô hình lực khí động                    Cơ học tấm composite                 Kết cấu thông minh áp điện
   +-------------------------+             +-------------------------+             +-------------------------+
   | Theodorsen (1935)       |             | Mindlin (1951) - FSDT   |             | Tiersten (1969)         |
   | Scanlan (1978, 1992)    |             | Reddy (2004) - HSDT     |             | Crawley & de Luis (1987)|
   +------------+------------+             +------------+------------+             +------------+------------+
                |                                       |                                       |
                +---------------------------------------+---------------------------------------+
                                                        |
                                                        v
                                   +-----------------------------------------+
                                   | CÁC NGHIÊN CỨU TIỀN NHIỆM:              |
                                   | - Balamurugan et al. (2002): Tấm gân PZT|
                                   |   (chỉ tải điều hòa, không khí động)    |
                                   | - Lee et al. (2002): Khí động Piston    |
                                   |   (chỉ tấm trơn, siêu âm M >= 2)        |
                                   | - Đào Huy Bích, Nguyễn Đình Đức:       |
                                   |   Phương pháp giải tích san đều gân     |
                                   +--------------------+--------------------+
                                                        |
                                                        | Research Gap
                                                        v
                                   +-----------------------------------------+
                                   | LUẬN ÁN NGUYỄN NGỌC THỦY (2020):        |
                                   | Mô hình hóa PTHH 9 nút phi tuyến FSDT   |
                                   | Tấm composite + gân rời rạc + PZT       |
                                   | Khí động Scanlan + Thực nghiệm hầm gió  |
                                   +-----------------------------------------+
  1. Dòng nghiên cứu về kết cấu tấm composite áp điện: Dipak Kumar Maiti (2008), Shen (2009), I. Oh et al. (2000), và Shunqi Zhang & Schmidt (2014) đã đặt nền móng giải thuật PTHH cho tấm composite tích hợp các lớp áp điện chịu tải cơ - nhiệt. Tuy nhiên, các công trình này chủ yếu khảo sát tấm phẳng đẳng hướng hoặc trực giao trực tiếp, bỏ qua sự tham gia của các kết cấu tăng cứng dạng gân rời rạc và hoàn toàn vắng bóng mô hình tải trọng khí động dòng chảy.
  2. Dòng nghiên cứu về khí đàn hồi và ổn định Flutter: Dựa trên mô hình kinh điển Theodorsen (1935) và Scanlan (1978, 1992), các tác giả Lee et al. (2002), Singha et al. (2001), và Makihara et al. (2007) đã phân tích hiện tượng mất ổn định Flutter của tấm cánh máy bay. Hạn chế lớn là các tác giả thường áp dụng lý thuyết Piston tuyến tính vốn chỉ phù hợp với dòng siêu âm vận tốc cao ($M \ge 2$), không phản ánh đúng tương tác xoáy và áp lực phân bố trong dải vận tốc cận âm và dưới âm ($M < 1$).
  3. Dòng nghiên cứu về tấm có gân gia cường: Tại Việt Nam, nhóm nghiên cứu của Đào Huy Bích (2010), Đào Văn Dũng (2013), và Nguyễn Đình Đức (2015) đã ứng dụng xuất sắc phương pháp giải tích kết hợp lý thuyết san đều tác dụng gân của Lekhnitsky. Mặc dù vậy, phương pháp san đều gân chỉ chính xác khi mật độ gân rất dày và kết cấu có điều kiện biên đơn giản (tựa bốn cạnh SSSS), không giải quyết được các bài toán thực tế khi gân bố trí cục bộ tại vị trí xung yếu, biên ngàm phức tạp (CCCC, CFFF) hoặc khi tấm có gắn các phiến áp điện phân tán.

Về mặt học thuật, tồn tại cuộc tranh luận lớn giữa quan điểm mô hình hóa gân rời rạc (discrete stiffener formulation) và quan điểm san đều độ cứng gân (smeared stiffener approach). Luận án của Nguyễn Ngọc Thủy đã định vị chính xác ở việc bảo vệ và hiện thực hóa mô hình gân rời rạc trong không gian phần tử hữu hạn 9 nút đẳng tham số, cho phép xét chính xác độ lệch tâm giữa trục trung hòa của gân và mặt trung bình của tấm thông qua ma trận chuyển đổi $[T_{gx}]$ và $[T_{gy}]$, vượt qua các hạn chế giản lược hóa của trường phái Lekhnitsky.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  • So với công trình của Balamurugan et al. (2002) trên tạp chí Smart Materials and Structures, nghiên cứu của Balamurugan chỉ phân tích dao động tấm có gân dưới tải trọng điều hòa cơ học đơn thuần; luận án của Nguyễn Ngọc Thủy đã tiến xa hơn khi tích hợp trọn vẹn mô hình khí động học Scanlan với các hệ số dẫn xuất thực nghiệm $H_i^, A_i^$.
  • So với công trình của Lee et al. (2002) trên AIAA Journal, vốn sử dụng lý thuyết Piston phi tuyến cho tấm composite trơn chịu dòng khí vượt âm, luận án này đã giải quyết bài toán dòng khí dưới âm có góc tới $\alpha$ biến thiên bất kỳ trên kết cấu tăng cứng thực tế, có kiểm chứng trực tiếp bằng mô hình hầm gió khí động.

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

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

Luận án đã mở rộng Lý thuyết biến dạng trượt bậc nhất Mindlin (FSDT) và nguyên lý biến phân Hamilton tổng quát cho hệ ghép liên tục điện - cơ - khí động học phi tuyến. Phương trình cấu động lực học của hệ được thiết lập từ biểu thức năng lượng:

$$\delta \int_{t_1}^{t_2} \left( T_e + T_{eg} - U_e^M - U_e^E - U_{ge} + W_e^M + W_e^E \right) dt = 0$$

Từ nguyên lý này, luận án đã dẫn xuất tường minh hệ phương trình vi phân phi tuyến toàn cục mô tả dao động của tấm composite áp điện có gân gia cường:

$$[M]{\ddot{q}} + [C]{\dot{q}} + [K(q)]{q} = {R(t)}$$

Trong đó, các ma trận thành phần được cấu thành từ sự kết hợp đa trường:

  • $[M] = [M_e^{uu}]$: Ma trận khối lượng kết cấu tổng thể bao gồm khối lượng tấm và khối lượng gân gia cường.
  • $[C] = [C_R] + [C_A] + [C_{air}]$: Ma trận cản toàn phần, tích hợp đồng thời cản kết cấu Rayleigh ($[C_R] = \alpha_r [M] + \beta_r [K]$), cản áp điện chủ động $[C_A]$ sinh ra từ mạch hồi tiếp vận tốc, và ma trận cản khí động học $[C_{air}]$ từ mô hình Scanlan.
  • $[K] = [K_L] + [K_N(q)] + [K_G] - [K_{air}] + [K_{u\phi}][K_{\phi\phi}]^{-1}[K_{\phi u}]$: Ma trận độ cứng tổng thể, biểu diễn sự tương tác giữa độ cứng đàn hồi tuyến tính $[K_L]$, độ cứng phi tuyến chuyển vị lớn $[K_N]$, độ cứng hình học $[K_G]$, độ cứng khí động $[K_{air}]$, và độ cứng tương tác cơ - điện áp điện.

Luận án đã đóng góp một bước chuyển dịch mô hình (paradigm shift): chứng minh rằng ứng xử phi tuyến hình học Von-Kármán đóng vai trò triệt tiêu hiện tượng phân kỳ biên độ vô hạn dự báo bởi các mô hình tuyến tính, đưa kết cấu về trạng thái dao động chu trình giới hạn (Limit Cycle Oscillations - LCO) ổn định hơn trong vùng vận tốc siêu tới hạn.

       MA TRẬN ĐỘ CỨNG VÀ CẢN TOÀN CỤC CỦA HỆ GHÉP NỐI ĐA TRƯỜNG
       
  +-------------------------------------------------------------------+
  |               MA TRẬN ĐỘ CỨNG HIỆU DỤNG TOÀN PHẦN [K*]           |
  |                                                                   |
  |   [K*] = [K_L] + [K_N(q)] + [K_G] - [K_air] + [K_uphi][K_phiphi]^-1[K_phi u]
  |            |        |        |        |                |          |
  |         Tuyến   Phi tuyến  Hình   Khí động       Ghép nối cơ-điện |
  |          tính   Von-Kármán  học   (Scanlan)      (Cảm biến/Actuator)
  +-------------------------------------------------------------------+
  |               MA TRẬN CẢN TOÀN CỤC [C]                            |
  |                                                                   |
  |   [C]  =  [C_R] (Rayleigh)  +  [C_A] (Áp điện)  +  [C_air] (Khí động)
  +-------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp chặt chẽ 3 trụ cột lý thuyết: (1) Cơ học vật liệu composite nhiều lớp dị hướng, (2) Lý thuyết áp điện Maxwell - Voigt với các quan hệ cấu động cơ - điện thuận nghịch:

$${D} = [d]{\sigma} + [p]{E} \quad \text{và} \quad {\sigma} = [C]{\varepsilon} - [e]^T {E}$$

và (3) Mô hình tải trọng khí động học Scanlan với lực nâng khí động $l_h$ và mô men khí động $m_\alpha$:

$$l_h = \frac{1}{2}\rho_a U^2 B \left[ k H_1^(k)\frac{\dot{h}}{U} + k H_2^(k)\frac{B\dot{\alpha}}{U} + k^2 H_3^(k)\alpha + k^2 H_4^(k)\frac{h}{B} \right]$$

$$m_\alpha = \frac{1}{2}\rho_a U^2 B^2 \left[ k A_1^(k)\frac{\dot{h}}{U} + k A_2^(k)\frac{B\dot{\alpha}}{U} + k^2 A_3^(k)\alpha + k^2 A_4^(k)\frac{h}{B} \right]$$

Điểm sáng tạo cốt lõi trong khung phân tích là việc thiết lập quan hệ chuyển đổi hình học gân lệch tâm. Bằng cách định nghĩa ma trận chuyển đổi $[T_{gx}]$ và $[T_{gy}]$ liên kết các bậc tự do của gân tại tâm hình học với các bậc tự do trên mặt trung bình của tấm composite, luận án bảo toàn được tính liên tục tuyệt đối về chuyển vị và biến dạng trượt tại mặt tiếp xúc mà không cần tăng số lượng nút phần tử một cách nhân tạo.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi quan điểm nhận thức luận thực chứng (positivist paradigm) kết hợp giữa mô phỏng số học chuyên sâu (computational mechanics) và kiểm chứng thực nghiệm vật lý (empirical validation). Mô hình phần tử hữu hạn sử dụng phần tử tấm 9 nút bậc hai đẳng tham số (Lagrange/Serendipity). Mỗi nút phần tử tấm sở hữu 5 bậc tự do cơ học:

$${q_i^M} = [u_{0i}, v_{0i}, w_{0i}, \theta_{xi}, \theta_{yi}]^T \quad (i = 1 \div 9)$$

tương ứng với 45 bậc tự do cơ học cho mỗi phần tử, cùng với 1 bậc tự do điện thế $\phi_i$ cho mỗi nút thuộc lớp áp điện, tạo ra một phần tử đẳng tham số 9 nút đa trường 54 bậc tự do hoàn chỉnh.

       PHẦN TỬ TẤM 9 NÚT ĐẲNG THAM SỐ (54 BẬC TỰ DO ĐA TRƯỜNG)
       
                     s ^ 
                       |
                 7-----6-----5
                 |     |     |   Mỗi nút (i = 1..9):
                 |     |     |   - 5 Bậc tự do cơ học: [u0, v0, w0, θx, θy]
                 8-----9-----4   - 1 Bậc tự do điện: [ϕ]
                 |     |     |
                 |     |     |
                 1-----2-----3 ---> r
                 
      Tích hợp gân gia cường lệch tâm dọc theo các đường nút tham chiếu:
      Ma trận độ lệch tâm: [T_gx], [T_gy] gắn kết đường trung bình gân
      với mặt trung bình tấm, triệt tiêu giả thiết san đều Lekhnitsky.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số được thực hiện qua các bước tính toán lặp phi tuyến nghiêm ngặt:

  1. Rời rạc hóa miền tính toán: Sử dụng tích phân số Gauss-Legendre bậc $3 \times 3$ để tính toán các ma trận độ cứng màng $[A]$, ma trận uốn - xoắn $[B]$, ma trận uốn $[D]$ và độ cứng cắt $[A_s]$, khắc phục hoàn toàn hiện tượng khóa cắt (shear locking).
  2. Giải bài toán phi tuyến động lực học: Sử dụng thuật toán Newmark-$\beta$ gia tốc trung bình ($\beta = 0.25, \gamma = 0.5$) kết hợp phương pháp lặp Newton-Raphson tại mỗi bước thời gian $\Delta t = 10^{-4} \div 10^{-3}$ s để đảm bảo sự hội tụ tuyệt đối của véc tơ dư tải trọng:

$$|{R_{residue}}| \le 10^{-4} \cdot |{R_{ext}}|$$

  1. Xác định ranh giới ổn định: Áp dụng Tiêu chuẩn mất ổn định động Budiansky-Roth. Theo tiêu chuẩn này, đồ thị đáp ứng độ võng cực đại $w_{max}$ được khảo sát theo tham số vận tốc dòng khí $U$. Điểm uốn (inflection point) với độ dốc $dw_{max}/dU \to \infty$ hoặc bước nhảy vọt biên độ dao động xác định chính xác vận tốc mất ổn định tới hạn $U_{cr}$.
     TIÊU CHUẨN MẤT ỔN ĐỊNH ĐỘNG BUDIANSKY - ROTH TRONG LUẬN ÁN
     
      Độ võng cực đại w_max
             ^
             |                                    / (Mất ổn định Flutter)
             |                                   /
             |                                  /
             |                                 /  <-- Điểm uốn (Inflection Point)
             |                                *       dw_max / dU -> vô cùng
             |                              / |       ==> VẬN TỐC TỚI HẠN U_cr
             |                             /  |
             |                            /   |
             |---------------------------/    |
             |  (Vùng dao động ổn định)       |
             +--------------------------------+------------------------->
             0                               U_cr            Vận tốc dòng khí U

Data và phân tích

Toàn bộ thuật toán được tác giả tự lập trình phát triển thành bộ phần mềm độc lập SSP_2018 (Smart Stiffened Plate 2018) với hàng chục nghìn dòng lệnh mã nguồn, được đối chuẩn (benchmarked) với phần mềm thương mại ANSYS và các công bố quốc tế của Gupta (2010) và Pradyumna (2010), cho sai số tần số dao động riêng dưới 1.2%.

Hệ thống thực nghiệm được thiết kế đồng bộ và vận hành tại Phòng thí nghiệm Trọng điểm Quốc gia:

  • Hầm gió khí động: Vận tốc dòng khí điều khiển chính xác từ $U = 0 \div 25$ m/s.
  • Bộ kích thích và khuếch đại: Bộ khuếch đại điện áp Piezo tuyến tính chuyên dụng EPA-104-230 ($V_{out} = 0 \div 200$ V).
  • Thiết bị tạo tín hiệu: Bộ phát tín hiệu chuẩn Onsoku ($f = 0.1 \div 1000$ Hz).
  • Hệ thống đo lường: Cảm biến gia tốc áp điện siêu nhẹ, Bộ thu nhận dữ liệu 2 kênh áp điện HnB75B kết nối máy hiện sóng số GW Instek GDS-2104 (băng thông 100 MHz, tốc độ lấy mẫu 1 GSa/s).
                        SƠ ĐỒ HỆ THỐNG ĐO KIỂM THỰC NGHIỆM HẦM GIÓ
                        
  +--------------------+      Tín hiệu điều hòa      +-----------------------+
  | BỘ PHÁT TÍN HIỆU   | --------------------------> | BỘ KHUẾCH ĐẠI PIEZO   |
  | CHUẨN ONSOKU       |     (Vin = 9.30 V,          | TUYẾN TÍNH EPA-104-230|
  | (f = 6.944 Hz)     |      f = 6.944 Hz)          +-----------+-----------+
  +--------------------+                                         |
                                                                 | Điện áp kích thích (Vout)
                                                                 v
  +----------------------------------------------------------------------------------+
  | HẦM GIÓ KHÍ ĐỘNG HỌC (Vận tốc U = 5 m/s, Góc tới α = 0°..45°)                    |
  |                                                                                  |
  |   Dòng khí U ===> [ TẤM COMPOSITE CÓ GÂN GIA CƯỜNG + TẤM PZT ACTUATOR/SENSOR ]   |
  |                                     |                                            |
  +-------------------------------------+--------------------------------------------+
                                        | Đáp ứng gia tốc / điện thế cảm ứng
                                        v
  +--------------------+      Tín hiệu đo đạc        +-----------------------+
  | MÁY HIỆN SÓNG      | <-------------------------- | BỘ THU NHẬN DỮ LIỆU   |
  | GW INSTEK GDS-2104 |                             | 2 KÊNH ÁP ĐIỆN HnB75B |
  +--------------------+                             +-----------------------+

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

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

Kết quả mô phỏng số và khảo sát thực nghiệm trong luận án đã mang lại 5 phát hiện mang tính đột phá cho ngành Cơ học kết cấu:

+--------------------------------------------------------------------------------------------------+
|                   TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ TỪ MÔ PHỎNG VÀ THỰC NGHIỆM                     |
+---+-----------------------------------+----------------------------------------------------------+
|STT| Tham số khảo sát                  | Kết quả định lượng & Quy luật vật lý                     |
+---+-----------------------------------+----------------------------------------------------------+
| 1 | Hiệu ứng phi tuyến Von-Kármán     | Triệt tiêu phân kỳ vô hạn; biên độ LCO giảm 18.5% so với |
|   |                                   | lý thuyết tuyến tính cổ điển.                            |
| 2 | Kích thước gân gia cường (h_g)    | Tăng h_g từ 0.002m lên 0.006m giúp nâng U_cr thêm 38.2%, |
|   |                                   | độ võng w_max giảm tới 54.6%.                            |
| 3 | Điện áp điều khiển PZT (V)        | Áp đặt V = -100V đến +100V kiểm soát chủ động biên độ    |
|   |                                   | dao động; dập tắt 62.4% năng lượng rung động.            |
| 4 | Góc đặt sợi composite (θ)         | Cấu hình góc đặt sợi [45°/-45°]s cho độ cứng uốn-xoắn    |
|   |                                   | cực đại, tối ưu hóa khả năng kháng Flutter.              |
| 5 | Kiểm chứng thực nghiệm hầm gió    | 5 lần đo lặp tại U = 5 m/s, f = 6.944 Hz cho sai số đáp  |
|   |                                   | ứng gia tốc so với tính toán lý thuyết chỉ 6.8%.         |
+---+-----------------------------------+----------------------------------------------------------+
  1. Hiệu ứng kìm hãm phi tuyến hình học: Khác với phân tích tuyến tính vốn dự báo biên độ dao động tiến tới vô cùng khi vận tốc đạt ngưỡng Flutter, mô hình phi tuyến Von-Kármán chỉ ra rằng lực màng kéo căng (membrane tension) sinh ra khi chuyển vị lớn làm tăng độ cứng hiệu dụng của tấm, kìm hãm biên độ dao động dừng lại ở chu trình giới hạn LCO ổn định, giảm 18.5% giá trị chuyển vị đỉnh so với tính toán tuyến tính.
  2. Tác động vượt trội của gân gia cường lệch tâm: Khi tăng chiều cao gân $h_g$ từ 0.002 m lên 0.006 m, vận tốc mất ổn định tới hạn $U_{cr}$ tăng từ 14.5 m/s lên 20.04 m/s (tăng 38.2%), trong khi độ võng cực đại $w_{max}$ tại tâm tấm giảm 54.6% dưới cùng điều kiện dòng khí. Gân gia cường không chỉ làm tăng độ cứng chống uốn mà còn phân tách phổ tần số riêng, ngăn chặn hiện tượng cộng hưởng cặp mode (mode coalescence).
  3. Cơ chế điều khiển chủ động của điện áp PZT: Việc áp đặt điện thế trái dấu lên các lớp áp điện tạo ra mô men uốn ngược chiều biến dạng cơ học, cho phép hiệu chỉnh độ võng tĩnh và giảm biên độ dao động động học tới 62.4%. Khi vận tốc dòng khí tiệm cận vùng tới hạn $U \approx U_{cr}$, việc kích hoạt mạch hồi tiếp vận tốc $G_v$ đã mở rộng miền ổn định khí đàn hồi thêm 22.8% dải vận tốc an toàn.
  4. Hiệu ứng ghép nối giữa góc đặt cốt $\theta$ và góc tới $\alpha$: Khảo sát trên dải góc sợi $\theta \in [0^\circ, 90^\circ]$ và góc tới của dòng khí $\alpha \in [0^\circ, 90^\circ]$ chứng minh rằng cấu hình sợi $[45^\circ/-45^\circ]_s$ cung cấp độ cứng chống xoắn và độ bền khí đàn hồi tối ưu nhất trước các luồng khí lệch trục. Ngược lại, góc tới $\alpha = 45^\circ$ gây ra trạng thái ứng suất cắt kết hợp uốn phức tạp nhất, làm suy giảm vận tốc tới hạn của tấm.
  5. Độ tin cậy của kiểm chứng thực nghiệm: Kết quả đo dao động trong hầm gió với mẫu tấm composite kích thước chuẩn chịu dòng khí $U = 5$ m/s, tần số kích thích $f = 6.944$ Hz và điện áp đầu vào $V_{in} = 9.30$ V cho thấy dạng sóng và phổ gia tốc thực nghiệm trùng khớp chặt chẽ với kết quả mô phỏng số từ chương trình SSP_2018, sai số tương đối cực đại không vượt quá 6.8% qua 5 lần đo lặp liên tiếp.

Implications đa chiều

  • Về mặt học thuật: Cung cấp cơ sở dữ liệu số và công thức giải tích giải quyết bài toán tương tác đa trường (Multiphysics FSI - Electro-mechanics) trên kết cấu tấm tăng cứng, mở ra hướng nghiên cứu mới cho các cấu trúc thông minh.
  • Về mặt kỹ thuật hàng không - quân sự: Cho phép các kỹ sư thiết kế trực tiếp ứng dụng để tối ưu hóa cánh đuôi tên lửa, cánh máy bay không người lái (UAV), tấm vỏ tàu cao tốc composite nhằm triệt tiêu hiện tượng rung cánh (flutter suppression) mà không làm tăng khối lượng tổng thể của kết cấu.
  • Về mặt chính sách và công nghệ quốc phòng: Đóng góp vào việc tự chủ hóa công nghệ phần mềm tính toán kết cấu chuyên sâu (chương trình SSP_2018), giảm phụ thuộc vào các phần mềm thương mại nước ngoài trong các dự án phát triển khí tài quân sự chiến lược.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu mang tính biên (boundary conditions):

  1. Mô hình khí động học Scanlan áp dụng chủ yếu phản ánh chính xác dòng khí dưới âm ($M < 1$), chưa bao hàm các hiệu ứng sóng xung kích (shock waves) và tương tác nhiệt khí động học trong dải vận tốc siêu thanh ($M > 1$) và cực siêu thanh ($M > 5$).
  2. Giả thiết liên kết giữa lớp composite, gân gia cường và phiến PZT là bám dính tuyệt đối, chưa xét đến hiện tượng phá hủy bóc tách lớp (delamination) và suy giảm đặc tính áp điện do mỏi nhiệt - cơ dài hạn.
  3. Thử nghiệm thực nghiệm mới chỉ thực hiện được tại vận tốc dòng khí dưới 10 m/s do giới hạn công suất của hầm gió phòng thí nghiệm, chưa đo trực tiếp được điểm mất ổn định Flutter thực tế trong ống thổi.

Chương trình nghiên cứu tương lai (Future Agenda):

  • Phát triển mô hình tích hợp Lý thuyết biến dạng cắt bậc cao (HSDT - Reddy) kết hợp mô hình chất lưu số 3D (CFD-FEM FSI coupling) để khảo sát kết cấu trong dải vận tốc chuyển âm và siêu âm.
  • Nghiên cứu cơ chế suy thoái và phá hủy cấu trúc áp điện phi tuyến xét đến hiệu ứng trễ (hysteresis) của vật liệu PZT dưới điện trường cao.
  • Ứng dụng các thuật toán điều khiển thông minh tự thích nghi (Adaptive Neural Network Control, Deep Reinforcement Learning) để điều khiển thời gian thực (real-time active control) hệ thống PZT trên các cánh bay biến hình (morphing wings).

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

Nghiên cứu của tác giả Nguyễn Ngọc Thủy tạo ra những giá trị tác động lan tỏa sâu rộng:

  • Tác động học thuật: Dự báo tiềm năng trích dẫn cao trong các tạp chí chuyên ngành hàng đầu như Composite Structures, Journal of Fluids and Structures, và Smart Materials and Structures, đóng vai trò tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh chuyên ngành Cơ học kỹ thuật và Kỹ thuật Hàng không.
  • Chuyển đổi công nghiệp và quốc phòng: Trực tiếp hỗ trợ Tổng cục Công nghiệp Quốc phòng và các Viện nghiên cứu Hàng không trong việc thiết kế kết cấu vỏ nhẹ chịu áp lực động cho các phương tiện không người lái và khí tài quân sự thế hệ mới.
  • Lợi ích kinh tế - kỹ thuật: Việc chủ động tối ưu hóa chiều cao gân $h_g$ kết hợp bố trí lớp PZT giúp tiết kiệm 20-30% khối lượng vật liệu composite sử dụng, đồng thời kéo dài tuổi thọ mỏi của kết cấu lên hơn 40% nhờ triệt tiêu dao động cộng hưởng biên độ lớn.

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

                              CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI
                              
  +-------------------------------------+     +-------------------------------------+
  |     NCS & HỌC GIẢ CƠ HỌC SỐ         |     |      KỸ SƯ THIẾT KẾ HÀNG KHÔNG      |
  | - Bộ mã nguồn mở thuật toán SSP_2018|     | - Quy trình tối ưu hóa gân & PZT    |
  | - Mô hình phần tử 9 nút FSDT đa pha |     | - Cắt giảm 20-30% khối lượng vỏ cánh|
  +-------------------------------------+     +-------------------------------------+
                                           \ /
                                            X
                                           / \
  +-------------------------------------+     +-------------------------------------+
  |     R&D CÔNG NGHIỆP QUỐC PHÒNG      |     |     CHUYÊN GIA KHÍ ĐÀN HỒI (AERO)   |
  | - Tự chủ công nghệ thiết kế UAV     |     | - Dữ liệu thực nghiệm hầm gió chuẩn |
  | - Triệt tiêu rung động cánh tên lửa |     | - Tiêu chuẩn Budiansky-Roth phi tuyến|
  +-------------------------------------+     +-------------------------------------+
  1. Nghiên cứu sinh và Giảng viên ngành Cơ kỹ thuật: Thụ hưởng một hệ thống lý thuyết chuẩn mực, công thức biến phân đa trường chi tiết và bộ thuật toán phần tử hữu hạn 9 nút hoàn chỉnh có thể mở rộng cho các bài toán vỏ trụ, vỏ nón hoặc vật liệu FGM.
  2. Kỹ sư R&D Hàng không - Vũ trụ và Công nghiệp Quốc phòng: Nắm bắt được phương pháp tính toán và thiết kế tối ưu hệ thống gân tăng cứng kết hợp cảm biến/kích hoạt áp điện, giúp nâng cao độ ổn định khí đàn hồi cho thân cánh phương tiện bay.
  3. Các cơ sở thử nghiệm và kiểm định kết cấu: Tiếp cận quy trình đo đạc thực nghiệm hầm gió kết nối hệ thống DAQ và máy hiện sóng số, làm chuẩn đối sánh cho các thí nghiệm khí đàn hồi vật lý.

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ì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập thành công mô hình phần tử hữu hạn 9 nút đẳng tham số ghép nối toàn phần Điện - Cơ - Khí động học phi tuyến cho tấm composite có gân gia cường lệch tâm. Luận án đã mở rộng Lý thuyết biến dạng trượt bậc nhất Mindlin (FSDT) kết hợp quan hệ hình học phi tuyến Von-Kármán và mô hình khí động học tự kích thích Scanlan, đưa thêm các ma trận chuyển đổi lệch tâm $[T_g]$ để mô tả chính xác tương tác giữa hệ gân rời rạc và tấm vỏ đa lớp, vượt qua giới hạn của lý thuyết san đều độ cứng Lekhnitsky kinh điển.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu tiền nhiệm?

So với các nghiên cứu của Balamurugan et al. (2002) (chỉ xét tải cơ học điều hòa) và Lee et al. (2002) (chỉ xét tấm phẳng trơn với lý thuyết Piston siêu âm), luận án đã:

  • Xây dựng phần tử tấm 9 nút với 54 bậc tự do đa trường có khả năng tích hợp linh hoạt các dải gân và phiến PZT tại các tọa độ bất kỳ.
  • Tích hợp mô hình khí động Scanlan dưới âm với đầy đủ các hệ số dẫn xuất thực nghiệm $H_i^, A_i^$.
  • Ứng dụng tiêu chuẩn Budiansky-Roth trên chuỗi thời gian giải bằng thuật toán Newmark-$\beta$ kết hợp Newton-Raphson để xác định chính xác điểm mất ổn định động phi tuyến.

3. Phát hiện thực nghiệm/mô phỏng gây bất ngờ nhất là gì?

Phát hiện bất ngờ nhất là hiệu ứng tương tác phi tuyến mạnh mẽ giữa góc tới của dòng khí $\alpha$ và góc đặt sợi composite $\theta$ lên ranh giới mất ổn định. Trong khi trực giác thông thường cho rằng góc tới dòng khí $\alpha = 0^\circ$ (dòng khí thổi dọc trục) luôn cho trạng thái ổn định nhất, phân tích số chỉ ra rằng với các tấm có cấu hình góc đặt sợi bất đối xứng, dòng khí có góc tới $\alpha = 45^\circ$ kết hợp với hiệu ứng xoắn của gân lệch tâm có thể làm suy giảm vận tốc tới hạn $U_{cr}$ lên tới 32.4%, tạo ra hiện tượng mất ổn định uốn - xoắn kết hợp sớm hơn nhiều so với dự báo của các mô hình tuyến tính.

4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (replication protocol) không?

Luận án cung cấp giao thức tái lập nghiên cứu hoàn chỉnh 100%:

  • Công bố chi tiết toàn bộ các ma trận phần tử ($[M_e]$, $[K_e]$, $[C_e]$, $[K_{u\phi}]$) và ma trận chuyển đổi gân $[T_{gx}], [T_{gy}]$ trong hệ tọa độ Gauss $(r, s)$.
  • Bảng thông số cơ lý tính vật liệu chi tiết: Tấm composite Graphite/Epoxy ($E_1 = 138$ GPa, $E_2 = 8.96$ GPa, $G_{12} = 7.1$ GPa, $\nu_{12} = 0.3$), Gân gia cường ($E_g = 138$ GPa), Lớp áp điện PZT-G1195N ($E_p = 63$ GPa, $d_{31} = d_{32} = -254 \times 10^{-12}$ m/V).
  • Mã nguồn chương trình SSP_2018 được đính kèm tường minh trong Phụ lục luận án.
  • Sơ đồ mạch đo, thông số thiết bị hầm gió, bộ khuếch đại EPA-104-230 và cảm biến được định lượng chính xác.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?

Luận án định hình lộ trình nghiên cứu 10 năm xoay quanh 3 hướng đột phá:

  1. Giai đoạn 2021-2025: Mở rộng mô hình giải thuật cho tấm và vỏ composite thông minh làm từ vật liệu biến đổi chức năng (FGM) và vật liệu gia cường ống nano carbon (CNTRC) chịu tải đa trường nhiệt - cơ - khí động học.
  2. Giai đoạn 2026-2028: Phát triển hệ thống điều khiển phản hồi thích nghi thời gian thực (Real-time FPGA-based Active Control) tích hợp vật liệu áp điện biến dạng trượt và sợi áp điện Macro Fiber Composite (MFC) cho cánh bay không người lái.
  3. Giai đoạn 2029-2030: Tích hợp mô hình học máy và mạng nơ-ron sâu (Physics-Informed Neural Networks - PINNs) để dự báo suy giảm tuổi thọ mỏi và cảnh báo sớm hiện tượng mất ổn định khí đàn hồi trên các phương tiện bay siêu âm thế hệ mới.

Kết luận

Luận án tiến sĩ kỹ thuật của NCS Nguyễn Ngọc Thủy là một công trình khoa học công phu, mẫu mực, giải quyết trọn vẹn bài toán cơ học đa trường phức tạp từ lý thuyết giải tích, mô phỏng số đến thực nghiệm vật lý. Các đóng góp cốt lõi của công trình bao gồm:

  1. Xây dựng thành công hệ phương trình vi phân phi tuyến toàn cục mô tả dao động và ổn định của tấm composite áp điện có gân gia cường chịu tải trọng khí động dưới âm theo mô hình Scanlan, tích hợp lý thuyết biến dạng trượt bậc nhất Mindlin và phi tuyến Von-Kármán.
  2. Đề xuất và hiện thực hóa mô hình phần tử hữu hạn 9 nút đẳng tham số 54 bậc tự do với ma trận chuyển đổi lệch tâm gân rời rạc $[T_g]$, khắc phục triệt để các nhược điểm của phương pháp giải tích san đều gân Lekhnitsky.
  3. Phát triển hoàn chỉnh bộ chương trình phần mềm tính toán chuyên dụng SSP_2018, cho phép khảo sát sâu sắc đáp ứng động lực học phi tuyến, phân tích chu trình giới hạn LCO và xác định ranh giới ổn định động theo tiêu chuẩn Budiansky-Roth.
  4. Khám phá và định lượng hóa các quy luật tương tác phức hợp giữa các thông số kết cấu ($h_g, h/W, \theta$), điện trường điều khiển ($V$) và dòng khí ($U, \alpha$), chứng minh khả năng tăng vận tốc tới hạn $U_{cr}$ lên trên 38% nhờ hệ gân tối ưu và giảm hơn 60% biên độ rung động nhờ kích hoạt áp điện chủ động.
  5. Thực hiện thành công hệ thống thí nghiệm hầm gió khí động học, kiểm chứng và xác nhận độ chính xác cao của mô hình số với sai số thực nghiệm chỉ 6.8%.
  6. Đặt nền móng khoa học vững chắc cho việc ứng dụng vật liệu và kết cấu composite thông minh vào thực tiễn thiết kế, chế tạo vũ khí công nghệ cao và phương tiện bay trong ngành Công nghiệp Quốc phòng Việt Nam.