Tổng quan về luận án

Nghiên cứu tối ưu hóa khí động học cho các phương tiện bay quân sự có điều khiển đóng vai trò cốt lõi trong việc nâng cao hiệu quả tác chiến và hiện đại hóa vũ khí trang bị. Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã số: 9 52 01 01) của nghiên cứu sinh Trần Mạnh Tuân, với đề tài "Nghiên cứu tối ưu biên dạng khí động của khí cụ bay điều khiển một kênh", được thực hiện tại Viện Tên lửa thuộc Viện Khoa học và Công nghệ quân sự - Bộ Quốc phòng, dưới sự hướng dẫn khoa học của PGS. TS Nguyễn Văn Chúc và TS Lê Đức Hạnh (Hà Nội, 2022). Công trình giải quyết bài toán cấp thiết: tối ưu hóa biên dạng hình học khí động cho tổ hợp khí cụ bay (KCB) điều khiển một kênh sử dụng dây vi cáp (điển hình như mẫu KCB CT14M thuộc tổ hợp CT11) nhằm đáp ứng đồng thời các yêu cầu khắt khe về chất lượng khí động, độ ổn định động học và khả năng cơ động bám quỹ đạo dẫn.

Khoảng trống nghiên cứu (research gap) lớn trong y văn hiện đại xuất phát từ việc các nghiên cứu khí động học trước đây phần lớn tập trung vào phương tiện bay điều khiển đa kênh truyền thống hoặc tên lửa không dây, bỏ qua các ràng buộc đặc thù của KCB điều khiển một kênh tự quay quanh trục dọc và tải trọng phi tuyến từ dây vi cáp. Như văn bản nguồn khẳng định: "Chưa có tài liệu nào công khai trình bày một cách đầy đủ, khoa học, có hệ thống về xây dựng mô hình toán mô phỏng động lực học bay cũng như tối ưu tham số biên dạng khí động cho đối tượng KCB điều khiển một kênh sử dụng dây vi cáp."

Luận án thiết lập 4 câu hỏi nghiên cứu (Research Questions - RQ) và các giả thuyết tương ứng (Hypotheses - H):

  • RQ1: Làm thế nào để mô hình hóa chính xác chuyển động 6 bậc tự do (6-DOF) của KCB điều khiển một kênh chịu tác động đồng thời của biến thiên khối lượng, lực phản lực Meshchersky và lực căng phi tuyến của dây vi cáp tời ra liên tục?
    • H1: Tải trọng cơ học và mô-men cản do lực căng dây vi cáp $F_d$ kết hợp với độ dịch chuyển khối tâm $x_T(t)$ làm suy giảm đáng kể độ dự trữ ổn định tĩnh $K_{od}$ và độ chính xác dẫn nếu không được bù trừ khí động học.
  • RQ2: Phương pháp tối ưu hóa nào giải quyết hiệu quả bài toán đa mục tiêu, đa ràng buộc với không gian tham số hình học phi tuyến mà không làm bùng nổ khối lượng tính toán?
    • H2: Phương pháp nghiên cứu không gian tham số (Parameter Space Investigation - PSI) kết hợp chuỗi phân bố chuẩn đều $LP_\tau$ cho phép xác định tập nghiệm Pareto tối ưu toàn cục với số bước lặp hội tụ thấp hơn đáng kể so với các thuật toán tiến hóa cổ điển.
  • RQ3: Các tham số thiết kế hình học cánh trước (vị trí gốc cánh $x_a$, sải cánh $L$, dây cung gốc $b_0$, dây cung mút $b_k$, góc mũi tên $\chi$) tác động tương hỗ ra sao đến hệ số chất lượng khí động $K = C_L / C_D$ và hệ số điều khiển $K_{dk}$?
  • RQ4: Phương án cấu hình cánh trước tối ưu có triệt tiêu được hiện tượng cộng hưởng tựa chu kỳ giữa tần số quay trục dọc $f_x$ và tần số dao động riêng kênh chúc ngóc $f_z$ trong toàn dải cự ly bắn từ 500 m đến 3.000 m không?

Khung lý thuyết của luận án tích hợp cơ học vật rắn biến khối lượng Meshchersky, khí động học dòng chảy nén được (Computational Fluid Dynamics - CFD), lý thuyết điều khiển tự động dẫn bắn 3 điểm (MCLOS/SACLOS) và lý thuyết tối ưu hóa đa mục tiêu Pareto. Nghiên cứu thực hiện khảo sát toàn diện trên mẫu KCB CT14M cải tiến trong không gian tham số đa chiều, giải bài toán quét không gian với các cỡ mẫu $N = 1024$ phương án trong miền tìm kiếm $P$, sau đó thu hẹp tinh chỉnh ở các miền $P_1$ ($N = 64$) và $P_2$ ($N = 64$ và $N = 128$), đem lại giá trị định lượng vượt trội về độ ổn định và góc tấn cân bằng.

Literature Review và Positioning

Nghiên cứu tối ưu hóa hình học khí động học cho phương tiện bay trên thế giới đã trải qua quá trình phát triển sâu rộng với nhiều trường phái phương pháp luận:

  1. Trường phái thuật toán di truyền (GA) và tiến hóa (EA): Tanil và cộng sự (2017) đã tích hợp phần mềm Missile Datcom (MD) và MATLAB Simulink để tối ưu hóa cấu hình khí động KCB với hàm mục tiêu cực tiểu hóa khối lượng phóng dưới các ràng buộc về tầm bay và vận tốc cực đại. Riddle (2004) áp dụng GA kết hợp phần mềm Aerodsn và MD nhằm cực tiểu khối lượng KCB nhiên liệu lỏng. Vidanovic và cộng sự (2012) ứng dụng thuật toán di truyền đa mục tiêu kết hợp mô phỏng số CFD để tối ưu hóa cánh nâng tên lửa siêu âm nhằm cực đại hóa tầm bắn và cực tiểu khối lượng cánh. Anderson (2001), Yang và cộng sự (2015), Runduo & Xiaobing (2016) mở rộng GA để tối ưu hóa đồng thời hệ số lực pháp tuyến $C_y$, hệ số lực cản $C_x$ và độ dự trữ ổn định tĩnh $K_{od}$.
  2. Trường phái tối ưu hóa bầy đàn (PSO) và Gradient: Usta và cộng sự (2014) sử dụng thuật toán PSO kết hợp phương pháp giải phương trình Navier-Stokes và Missile Datcom để tối ưu cánh lái KCB siêu âm, đạt trạng thái hội tụ sau 67 bước lặp. Xia và cộng sự (2019) áp dụng PSO giúp giảm hệ số lực cản 14% cho cánh cận âm và 15% cho KCB siêu âm. Feyzioglu (2017) sử dụng phương pháp gradient ngược bậc nhất kết hợp thuật toán tìm kiếm đường (Line Search) để cực tiểu vận tốc góc cánh lái.
  3. Trường phái thực nghiệm ống thổi và tham số hóa: Ocokoljic và cộng sự (2014) thực hiện tối ưu hình học phần mũi KCB cải tiến dựa hoàn toàn vào chuỗi thử nghiệm đắt đỏ trong ống thổi khí động học nhằm so sánh hệ số chất lượng khí động và độ ổn định so với mẫu CT14M nguyên bản.

Văn bản chỉ rõ hai luồng quan điểm đối nghịch trong tối ưu hóa khí động KCB: Một luồng quan điểm (đại diện bởi các phương pháp dựa trên gradient hoặc Monte Carlo truyền thống như Ahmed, 2010; Colonno, 2012) ưu tiên tốc độ hội tụ toán học cục bộ nhưng phụ thuộc sâu vào đạo hàm giải tích và điểm khởi tạo ban đầu, khiến việc tìm nghiệm tối ưu toàn cục trong không gian nhiều chiều bị hạn chế. Luồng quan điểm thứ hai (như Holland, 1960; Goldberg, 1989; Kennedy & Eberhart) ủng hộ các thuật toán phỏng sinh học (GA, PSO), tuy tìm được nghiệm toàn cục nhưng đòi hỏi duy trì quần thể lớn qua nhiều thế hệ, tiêu tốn tài nguyên tính toán khổng lồ và không thể tận dụng dữ liệu tính toán khi muốn mở rộng cỡ mẫu.

Luận án của NCS Trần Mạnh Tuân định vị chính xác ở điểm giao thoa giữa cơ học bay thực tế và tối ưu hóa hệ thống hiện đại: Luận án không tách rời bài toán khí động học thuần túy khỏi bài toán động lực học điều khiển, mà liên kết trực tiếp việc tính toán hệ số khí động (CFD, Missile Datcom) với hệ phương trình vi phân chuyển động 6-DOF hoàn chỉnh trong môi trường MATLAB Simulink. Bằng việc ứng dụng phương pháp nghiên cứu không gian tham số PSI kết hợp chuỗi phân bố chuẩn đều $LP_\tau$ của Sobol & Statnikov, nghiên cứu giải quyết triệt để bài toán đa mục tiêu phức tạp mà các công trình quốc tế như của Tanil (2017) hay Ocokoljic (2014) chưa bao quát được đối với KCB điều khiển một kênh tự quay.

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

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

Luận án tạo bước tiến khoa học quan trọng bằng việc kế thừa và mở rộng hệ thống lý thuyết động lực học bay cổ điển của Kolesnikov (1995), Korosteleb (1993) và Lebedev (1987) sang miền bài toán KCB điều khiển một kênh thế hệ I và II (tổ hợp CT11/CT14M).

Đóng góp lý thuyết cốt lõi được thể hiện qua các mệnh đề và mô hình toán sau:

Mệnh đề 1 (Động lực học biến khối lượng chịu lực căng dây): Chuyển động khối tâm và chuyển động quay quanh khối tâm của KCB điều khiển một kênh sử dụng dây vi cáp tuân theo hệ phương trình vector vi phân cơ học biến thiên: $$\tilde{F}_b = m(t)\left(\frac{d\tilde{V}_b}{dt} + \tilde{\omega}_b \times \tilde{V}_b\right) - \tilde{F}_M$$ $$\tilde{M}_b = \frac{d\tilde{L}_c}{dt} + \tilde{\omega}_b \times \tilde{L}_c$$ Trong đó, ngoại lực tổng hợp $\tilde{F}_b$ và ngoại mô-men $\tilde{M}b$ trong hệ tọa độ liên kết $Ox_by_bz_b$ chứa thành phần lực căng dây vi cáp $\tilde{F}{bd}$, mô-men do dây $\tilde{M}_d$, lực đẩy động cơ phóng $\tilde{F}_b^{PP}$, lực đẩy động cơ hành trình $\tilde{F}_b^{PHT}$, lực phản lực Mechersky $\tilde{F}M = \frac{dm{day}}{dt}\tilde{V}_b$, và các thành phần khí động $\tilde{F}_b^a, \tilde{M}^a$.

Mệnh đề 2 (Mô hình hóa lực căng và ma sát dây vi cáp): Lực cản tời dây và lực ma sát dọc ống dây được xác định qua hàm phi tuyến của vận tốc bay $V$, chiều dài dây tời $l_{day}$, khối lượng riêng dây $\rho_d$, đường kính dây $d$ và hệ số hình dạng ống dây $K_d(\lambda_d)$: $$F_d = F_{ms} + \rho_d V^2 + f(l_{day}, \lambda_d)$$ Mô hình chứng minh lực căng dây $F_d$ tác dụng lệch trục khối tâm tạo ra mô-men $M_d$ liên tục gây nhiễu loạn góc chúc ngóc và hướng, đòi hỏi hệ thống khí động phải có độ dự trữ ổn định tĩnh tối ưu $K_{od} > 0$.

Mệnh đề 3 (Độ dự trữ ổn định tĩnh và điều khiển cân bằng): Văn bản xác lập rõ ràng công thức định lượng tiêu chuẩn ổn định tĩnh: $$K_{od} = \frac{x_T - x_F}{L_{ref}}$$ với $x_T$ là tọa độ khối tâm, $x_F$ là tọa độ tâm áp tính từ đỉnh mũi KCB, $L_{ref}$ là chiều dài đặc trưng. Độ nhạy điều khiển được chuẩn hóa qua biểu thức: $$K_{dk} = \frac{\alpha_{cb}}{\delta_{cb}}$$ trong đó $\alpha_{cb}$ là góc tấn cân bằng và $\delta_{cb}$ là góc lệch cánh lái (hoặc góc lệch loa phụt tương đương $\delta_{cbtd}$).

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

Khung phân tích của luận án tích hợp liên hoàn ba lý thuyết nền tảng: (1) Khí động học tính toán dòng chảy nén dưới âm/trên âm; (2) Động lực học hệ thống điều khiển một kênh bám quỹ đạo 3 điểm (MCLOS/SACLOS); (3) Phương pháp quy hoạch toán học không gian tham số PSI ($LP_\tau$).

Cách tiếp cận độc đáo này tạo ra một vòng lặp kín (closed-loop optimization): Từ một điểm thiết kế $A(x_a, L, b_0, b_k, \chi)$ trong không gian cấu hình $P$, mô hình khí động tự động kết xuất bảng hệ số lực và mô-men ($C_x, C_y, C_z, m_x, m_y, m_z$), đưa trực tiếp vào bộ tích phân vi phân 6-DOF trong MATLAB Simulink để mô phỏng quỹ đạo bay thực dưới tác động của dây vi cáp, lực đẩy động cơ thay đổi theo thời gian và hệ thống tự quay. Dữ liệu đầu ra lập tức được kiểm tra qua các điều kiện biên nghiêm ngặt:

  • Ràng buộc động học chống cộng hưởng: Tần số quay quanh trục dọc $f_x \ge 2 f_z$ (với $f_z$ là tần số dao động riêng kênh chúc ngóc).
  • Ràng buộc quá tải động: Hệ số quá tải pháp tuyến $n_y \ge n_{yVmin}$ nhằm đảm bảo nâng trọng lượng KCB và duy trì chuyển động nằm ngang ở pha bay hành trình.
  • Giới hạn góc tấn tới hạn: $\alpha \le \alpha_{th}$ để ngăn ngừa hiện tượng mất lực nâng dòng khí.

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 thực nghiệm (positivism paradigm), kết hợp chặt chẽ giữa mô hình hóa toán học giải tích, tính toán mô phỏng số độ trung thực cao và thực nghiệm vật lý kiểm chứng. Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm:

  • Cấp độ 1 (Khí động vi mô): Tính toán phân bố trường áp suất bề mặt và các hệ số khí động $C_L, C_D, C_y^\alpha, m_z^\alpha, m_z^{\omega_z}$ thông qua phần mềm khí động chuyên dụng Missile Datcom và mô phỏng số CFD 3D (giải hệ phương trình Reynolds-Averaged Navier-Stokes).
  • Cấp độ 2 (Động lực học hệ thống vĩ mô): Mô phỏng thời gian thực chuyển động 6 bậc tự do của KCB với đầy đủ các pha: pha phóng (động cơ phóng ĐCP, thời gian cháy $t_{dcp}$, góc lệch loa phụt phóng $50'$ tạo mô-men quay ban đầu), pha giữ chậm $t_{gc}$, và pha hành trình (động cơ hành trình ĐCHT cháy trong thời gian $t_{dcht}$, cánh ổn định đuôi đặt lệch $3^\circ 15'$ duy trì vận tốc góc $\omega_x$, loa phụt ĐCHT lật đổi hướng tạo lực điều khiển vuông góc với trục dọc).
  • Cấp độ 3 (Tối ưu hóa đa không gian): Giải bài toán tối ưu hóa đa mục tiêu bằng phương pháp PSI.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ nghiêm ngặt chuẩn mực khoa học kỹ thuật quân sự:

  1. Mô hình hóa lực và mô-men: Hệ tọa độ sử dụng bao gồm Hệ tọa độ mặt đất $O_ex_ey_ez_e$, Hệ tọa độ liên kết $O x_b y_b z_b$, Hệ tọa độ vận tốc $O x_w y_w z_w$, và Hệ tọa độ cầu mặt đất $O_e r \varepsilon \sigma$. Ma trận chuyển đổi cô-sin chỉ phương $C_e^b, C_b^w, C_w^e$ được thiết lập chuẩn xác theo các góc Euler ($\psi, \vartheta, \gamma$) và góc khí động ($\alpha, \beta, \gamma_c$).
  2. Kỹ thuật thu thập và đối chiếu dữ liệu (Triangulation): Luận án kết hợp 3 nguồn kiểm chứng độc lập:
    • Dữ liệu tính toán lý thuyết và bán thực nghiệm bằng phần mềm chuyên dụng Missile Datcom.
    • Mô phỏng số khí động 3D CFD trên lưới phi cấu trúc với kiểm soát lớp biên nghiêm ngặt.
    • Thử nghiệm vật lý thực tế trên mô hình thổi khí động trong ống thổi khí động học (Wind Tunnel Testing) kết hợp đo lực căng dây vi cáp bằng cảm biến tải chuyên dụng.
  3. Độ tin cậy của thuật toán PSI: Sử dụng bảng hệ số phát sinh chuỗi chuẩn đều số học $LP_\tau$ (Sobol sequences). Điểm vượt trội của chuỗi $LP_\tau$ so với việc tạo số ngẫu nhiên Monte Carlo thuần túy là độ phân tán cực kỳ đồng đều trong không gian $n$ chiều, loại bỏ hoàn toàn hiện tượng vón cục hoặc khoảng trống dữ liệu.

Data và phân tích

  • Đặc tính mẫu khảo sát: Mô hình KCB CT14M nguyên bản có khối lượng phóng ban đầu $m_0$, đường kính thân, chiều dài đặc trưng $L_{ref}$, tầm bắn hiệu quả từ 500 m đến 3.000 m. Vận tốc góc quay trục dọc được thiết kế quanh ngưỡng đảm bảo ổn định con quay.
  • Không gian biến thiết kế cánh trước: 5 biến hình học độc lập được khảo sát trong các miền giới hạn kỹ thuật:
    • Vị trí điểm đầu dây cung gốc cánh: $x_a \in [x_a^*, x_a^{**}]$ (mm)
    • Sải cánh trước: $L \in [L^*, L^{**}]$ (mm)
    • Dây cung gốc cánh: $b_0 \in [b_0^*, b_0^{**}]$ (mm)
    • Dây cung mút cánh: $b_k \in [b_k^*, b_k^{**}]$ (mm)
    • Góc mũi tên cánh: $\chi \in [\chi^*, \chi^{**}]$ (độ)
  • Công cụ phân tích và phần mềm: Chương trình mô phỏng quỹ đạo và động lực học 6 bậc tự do được lập trình hoàn chỉnh trong môi trường MATLAB & Simulink, tích hợp tự động với bộ giải Missile Datcom và các script xử lý số liệu tối ưu PSI.

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

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

Quá trình phân tích dữ liệu và mô phỏng thực nghiệm mang lại 4 phát hiện khoa học mang tính bước ngoặt:

  1. Quy luật phân bố Pareto trong không gian mục tiêu: Khi quét 1024 điểm phân bố $LP_\tau$ trong không gian tham số ban đầu $P$, phần lớn các phương án bị loại bỏ do vi phạm điều kiện tự quay chống cộng hưởng ($f_x < 2 f_z$) hoặc độ dự trữ ổn định tĩnh âm ($K_{od} < 0$). Tập điểm thỏa mãn ràng buộc $G$ và tập các điểm chấp nhận được $D$ chỉ chiếm một tỷ lệ nhỏ, phân bố trên một dải hẹp dạng biên cong Pareto.
  2. Hiện tượng phi tuyến của vị trí cánh trước đối với độ dự trữ ổn định tĩnh: Khi dịch chuyển vị trí cánh trước $x_a$ về phía mũi KCB, tâm áp $x_F$ dịch chuyển mạnh về phía trước nhanh hơn sự dịch chuyển khối tâm $x_T$, làm giảm $K_{od}$. Tuy nhiên, điều này lại làm tăng đáng kể hệ số tính điều khiển được $K_{dk} = \alpha_{cb} / \delta_{cb}$, giúp giảm bán kính cong quỹ đạo $R = mV^2/F_{dk}$ khi tấn công mục tiêu ở cự ly gần (500 m).
  3. Sự đánh đổi (trade-off) giữa sải cánh $L$ và lực cản hành trình $C_x$: Tăng sải cánh $L$ và diện tích cánh trước giúp tăng hệ số lực nâng $C_L$, tạo lực điều khiển cân bằng $F_{cbdk}$ lớn hơn, nhưng đồng thời làm tăng hệ số lực cản cảm ứng và lực cản ma sát, dẫn đến suy giảm vận tốc cuối hành trình ở cự ly cực đại 3.000 m.
  4. Xác lập phương án thiết kế tối ưu $B_2$: Thông qua 3 bước lặp thu hẹp không gian tìm kiếm từ $P \rightarrow P_1 \rightarrow P_2$ ($N = 128$), phương án thiết kế $B_2$ được lựa chọn sở hữu bộ thông số vượt trội:
    • Hệ số chất lượng khí động $K = C_L / C_D$ tối ưu hóa ở dải góc tấn làm việc.
    • Độ dự trữ ổn định tĩnh $K_{od}$ luôn duy trì dương và ổn định trong khoảng an toàn dọc theo toàn bộ hành trình tiêu hao nhiên liệu của cả ĐCP và ĐCHT.
    • Thỏa mãn tuyệt đối điều kiện tần số dao động: $f_z$ đầu hành trình đạt xấp xỉ 2,17 Hz, cuối hành trình đạt 2,87 Hz, trong khi tần số quay quanh trục dọc $f_x$ luôn lớn hơn $2 f_z$, dập tắt hoàn toàn dao động ghép kênh.

Implications đa chiều

  • Về mặt lý thuyết học thuật: Công trình hoàn thiện cơ sở lý thuyết khí động lực học cho lớp phương tiện bay tự quay điều khiển một kênh có liên kết cơ học mềm (dây vi cáp), bổ sung một nhánh nghiên cứu chuyên sâu vào ngành Cơ kỹ thuật và Kỹ thuật Hàng không vũ trụ quân sự.
  • Về mặt phương pháp luận: Chứng minh tính ưu việt vượt bậc của phương pháp PSI sử dụng chuỗi $LP_\tau$ trong các bài toán kỹ thuật quân sự phức tạp: cho phép mở rộng gấp đôi số lượng mẫu thử ($N \rightarrow 2N$) mà không cần tính toán lại các điểm cũ, đồng thời giữ nguyên vai trò quyết định của chuyên gia thiết kế trong việc lựa chọn giới hạn hàm mục tiêu $\Phi_v^{**}$.
  • Về mặt thực tiễn và quốc phòng: Cung cấp trực tiếp bộ hồ sơ thiết kế hình học cánh trước hoàn chỉnh phục vụ dự án nâng cấp, cải tiến các tổ hợp KCB chống tăng hiện có trong biên chế Quân đội nhân dân Việt Nam, nâng cao độ chính xác bắn trúng mục tiêu di động và cố định ở dải cự ly từ 500 m đến 3.000 m với chi phí sản xuất thấp.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (limitations) mang tính khoa học:

  1. Giả thiết về tính đồng nhất của dòng khí quyển: Mô hình mô phỏng động lực học 6-DOF hiện tại chưa xét đến trường gió giật cục bộ phi tuyến mạnh hoặc sự thay đổi đột ngột của mật độ không khí theo thời tiết khắc nghiệt.
  2. Mô hình lực căng dây vi cáp: Lực căng dây $F_d$ và ma sát tời dây $F_{ms}$ được mô hình hóa theo các hàm bán thực nghiệm dựa trên các thử nghiệm đo đạc ở điều kiện tiêu chuẩn, chưa mô phỏng chi tiết tương tác tiếp xúc động lực học đa vật thể (multibody contact dynamics) của từng vòng dây trong cuộn cáp.
  3. Phạm vi khảo sát hình học: Luận án tập trung tối ưu hóa 5 thông số hình học của cụm cánh trước cải tiến, giữ nguyên hình dạng thân trụ cơ sở và cụm cánh đuôi nguyên bản của mẫu CT14M.

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

  • Hướng 1: Mở rộng bài toán tối ưu hóa đồng thời toàn bộ cấu hình khí động (thân, cánh trước, cánh ổn định đuôi và loa phụt) dưới điều kiện dòng chảy ba chiều siêu âm cao.
  • Hướng 2: Tích hợp thuật toán điều khiển thích nghi thông minh (Adaptive / Robust Control) trên khoang để tự động điều chỉnh góc lật loa phụt theo thời gian thực khi dây cáp bị đứt hoặc gặp gió nhiễu loạn mạnh.
  • Hướng 3: Nâng cấp mô hình dây vi cáp sang liên kết điều khiển quang cáp tốc độ cao, cho phép truyền dữ liệu hình ảnh độ nét cao về đài điều khiển bán tự động SACLOS thế hệ mới.
  • Hướng 4: Nghiên cứu ứng dụng các thuật toán tối ưu hóa đa mục tiêu lai ghép (Hybrid PSI-GA/PSO) nhằm tự động hóa hoàn toàn quá trình tìm kiếm biên Pareto trên máy tính hiệu năng cao (HPC).

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

  • Tác động học thuật: Luận án đóng góp hệ thống phương pháp luận mô phỏng - tối ưu hóa hoàn chỉnh, làm tài liệu tham khảo mẫu mực cho các chương trình đào tạo sau đại học chuyên ngành Cơ kỹ thuật, Cơ khí động lực, Vũ khí và Hàng không vũ trụ tại Viện Khoa học và Công nghệ quân sự, Học viện Kỹ thuật Quân sự và Đại học Bách khoa Hà Nội. 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 trong nước và quốc tế.
  • Tác động công nghiệp quốc phòng: Tạo tiền đề kỹ thuật vững chắc cho Viện Tên lửa và các nhà máy quốc phòng thuộc Tổng cục Công nghiệp Quốc phòng làm chủ công nghệ thiết kế, chế tạo, cải tiến các loại vũ khí có điều khiển chính xác cao, giảm phụ thuộc vào việc nhập khẩu trang bị nước ngoài.
  • Lợi ích kinh tế - xã hội và an ninh: Việc tận dụng và cải tiến thành công các tổ hợp khí cụ bay thế hệ I/II sẵn có giúp tiết kiệm hàng trăm tỷ đồng ngân sách quốc phòng so với mua sắm mới các tổ hợp tên lửa thế hệ thứ III đắt đỏ, đồng thời bảo đảm khả năng sẵn sàng chiến đấu bảo vệ chủ quyền lãnh thổ.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp thiết lập hệ phương trình vi phân chuyển động 6-DOF biến khối lượng và kỹ thuật tối ưu hóa đa mục tiêu PSI bằng chuỗi $LP_\tau$.
  • Các nhà khoa học và chuyên gia khí động học: Thụ hưởng mô hình toán học đã được kiểm chứng bằng thực nghiệm ống thổi và CFD đối với KCB tự quay có tải trọng dây vi cáp.
  • Kỹ sư thiết kế và R&D công nghiệp quốc phòng: Sở hữu ngay bộ thông số hình học cánh trước tối ưu ($B_2$) và chương trình mô phỏng MATLAB Simulink để áp dụng trực tiếp vào quy trình sản xuất vũ khí.
  • Cơ quan quản lý và chỉ huy quân sự: Có căn cứ khoa học xác đáng để phê duyệt các dự án nâng cấp, hiện đại hóa vũ khí trang bị kỹ thuật quân sự trong giai đoạn 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ì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng hệ phương trình động lực học bay 6 bậc tự do cổ điển của Kolesnikov và Lebedev cho lớp phương tiện bay tự quay điều khiển một kênh bằng việc bổ sung đồng thời vector lực biến khối lượng Meshchersky $\tilde{F}M$ và vector lực căng dây vi cáp phi tuyến $\tilde{F}d$, giải quyết chính xác sự tương tác giữa độ dự trữ ổn định tĩnh $K{od} = (x_T - x_F)/L{ref}$ và hệ số điều khiển cân bằng $K_{dk} = \alpha_{cb}/\delta_{cb}$.
  2. Phương pháp nghiên cứu có điểm gì đột phá so với ít nhất 2 công trình quốc tế tiêu biểu? Trả lời: So với nghiên cứu của Tanil và cộng sự (2017) (dùng GA tốn kém thời gian tính toán và dễ mắc kẹt) và Ocokoljic và cộng sự (2014) (chỉ dựa vào thử nghiệm ống thổi thuần khí động), luận án đã kết hợp phương pháp PSI sử dụng chuỗi chuẩn đều $LP_\tau$ với mô phỏng động lực học bay 6-DOF khép kín, cho phép khảo sát toàn diện 1024 phương án với khả năng nhân đôi cỡ mẫu mà không cần chạy lại các bước tính toán trước đó.
  3. Phát hiện nào mang tính bất ngờ nhất và có bằng chứng số liệu ra sao? Trả lời: Phát hiện bất ngờ nhất là việc tăng sải cánh trước không đồng biến với việc cải thiện quỹ đạo dẫn ở cự ly ngắn (500 m) do sự gia tăng đột biến của lực cản ma sát và sự trôi dịch tâm áp $x_F$, làm thay đổi tần số dao động riêng kênh chúc ngóc $f_z$. Chỉ khi phối hợp tối ưu đồng thời cả 5 thông số ($x_a, L, b_0, b_k, \chi$) tại phương án $B_2$, KCB mới vừa dập tắt nguy cơ cộng hưởng ($f_x > 2 f_z$), vừa bám sát đường ngắm 3 điểm.
  4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Có. Toàn bộ hệ phương trình vi phân chuyển động 6-DOF, ma trận cosin chỉ phương, bảng hệ số phát sinh chuỗi $LP_\tau$, cấu trúc chương trình tính toán trên MATLAB Simulink, cùng bản vẽ biên dạng hình học cánh trước và bảng số liệu khí động từ Missile Datcom đều được trình bày chi tiết trong phần chính văn và 8 phụ lục của luận án.
  5. Chương trình nghiên cứu 10 năm tới định hình những hướng đi nào? Trả lời: Phát triển hệ thống dẫn bắn bán tự động SACLOS truyền tín hiệu qua sợi quang học; tích hợp cánh lái điều khiển khí động kết hợp loa phụt xung lực; và ứng dụng trí tuệ nhân tạo (Deep Reinforcement Learning) trong dẫn bắn bám bắt mục tiêu cơ động cao.

Kết luận

Luận án tiến sĩ của NCS Trần Mạnh Tuân đã hoàn thành xuất sắc toàn bộ mục tiêu đề ra với các kết luận khoa học cô đọng:

  1. Xây dựng thành công mô hình toán học và chương trình mô phỏng số 6 bậc tự do hoàn chỉnh cho KCB điều khiển một kênh sử dụng dây vi cáp, tích hợp đầy đủ lực căng dây, lực biến khối lượng Meshchersky và lực điều khiển loa phụt.
  2. Thiết lập bài toán tối ưu hóa đa mục tiêu biên dạng khí động học với các hàm mục tiêu cốt lõi: chất lượng khí động $K$, độ dự trữ ổn định tĩnh $K_{od}$, tính điều khiển được $K_{dk}$ và bán kính cong cơ động $R$.
  3. Ứng dụng thành công phương pháp nghiên cứu không gian tham số PSI kết hợp chuỗi phân bố chuẩn đều $LP_\tau$, sàng lọc hiệu quả qua các không gian $P, P_1, P_2$ với cỡ mẫu lên tới $N = 1024$.
  4. Xác định thành công cấu hình cánh trước tối ưu $B_2$, nâng cao vượt bậc chất lượng bay, bảo đảm ổn định con quay $f_x \ge 2 f_z$, dập tắt dao động và đảm bảo độ chính xác dẫn bắn ở mọi cự ly tác chiến từ 500 m đến 3.000 m.
  5. Kiểm chứng độ tin cậy tuyệt đối của mô hình thông qua sự trùng khớp chặt chẽ giữa tính toán giải tích, mô phỏng số CFD và thực nghiệm thổi khí động học trong phòng thí nghiệm.
  6. Mở ra hướng nghiên cứu và ứng dụng thực tiễn to lớn trong việc nâng cấp, làm chủ công nghệ chế tạo các tổ hợp vũ khí có điều khiển chính xác cao, phục vụ đắc lực cho sự nghiệp xây dựng quân đội và củng cố quốc phòng an ninh quốc gia.