Tổng quan về luận án

Nghiên cứu động lực học bay của các phương tiện phóng và đạn phản lực đóng vai trò nền tảng trong khoa học hàng không vũ trụ và hiện đại hóa kỹ thuật quốc phòng. Trong các thiết kế truyền thống, tên lửa thường được đơn giản hóa như một vật rắn tuyệt đối (rigid body). Tuy nhiên, đối với dòng tên lửa không điều khiển thân mảnh (slender unguided rocket) có tầm bắn dưới 50 km, tải trọng khí động cực lớn khi bay qua tầng khí quyển đậm đặc kết hợp với lực đẩy phản lực sinh ra biến dạng uốn đàn hồi đáng kể dọc theo thân dầm. Biến dạng này làm thay đổi góc tấn cục bộ $\alpha(x, t)$, làm sai lệch hướng véc-tơ lực đẩy tại miệng loa phụt, và làm biến đổi phân bố lực nâng $F_L$ cùng lực cản $F_D$.

Khoảng trống nghiên cứu (research gap) then chốt trong y văn quốc tế và trong nước nằm ở việc phần lớn các mô hình trước đây (chẳng hạn như Xiao Huang & Zeiler, 2006; Xu Yun, Xie Changchuan & Yang Chao, 2008) tiếp cận bài toán cơ hệ biến đổi khối lượng một cách gần đúng: xem lực đẩy phản lực $F_{pl} = \beta U$ như một ngoại lực tập trung độc lập và chỉ cập nhật khối lượng tức thời của tên lửa theo thời gian mà bỏ qua tính phi bảo toàn (non-conservative) nội tại của hệ hở. Cách tiếp cận này triệt tiêu các thành phần lực liên kết động học vuông góc với trục dọc sinh ra do tương tác giữa lưu lượng phụt $\beta$ và vận tốc, gia tốc dao động đàn hồi tại miệng phụt, đồng thời chưa tích hợp đồng thời sự dịch chuyển khối tâm thực tế $G(t)$.

Luận án tiến sĩ của nghiên cứu sinh Đàm Việt Phương (người hướng dẫn khoa học: PGS.TS Vũ Quốc Trụ, Học viện Kỹ thuật Quân sự) tập trung giải quyết 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 (RQ1): Cơ chế toán - cơ học nào mô tả chính xác tương tác giữa chuyển động không gian 3 chiều và biến dạng đàn hồi uốn của cơ hệ dầm biến đổi khối lượng liên tục mang luồng phụt?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ suy giảm lực đẩy hiệu dụng $F_{plhd}$ do góc xoay đàn hồi tại loa phụt và sự trượt khối tâm dọc trục tác động định lượng như thế nào đến quỹ đạo bay, vận tốc, góc chúc ngóc $\vartheta$ và góc tấn $\alpha$?
  • Giả thuyết nghiên cứu 1 (H1): Việc mở rộng nguyên lý Hamilton và phương trình Lagrange cho hệ thống kín gồm thân tên lửa và vi phân khối lượng luồng phụt $dm = \beta dt$ sẽ phát hiện các thành phần lực tương tác mới vuông góc với trục thân mà mô hình vật rắn hoặc hệ Lagrange quy ước bỏ qua.
  • Giả thuyết nghiên cứu 2 (H2): Gia tăng độ cứng chống uốn $EJ$ từ giá trị kết cấu thực $K_{gd}$ lên $100K_{gd}$ sẽ tiệm cận nghiệm vật rắn tuyệt đối, chứng minh dao động uốn là nguyên nhân trực tiếp gây biến thiên góc chúc ngóc cục bộ và tản mát quỹ đạo đạn.

Khung lý thuyết (theoretical framework) của công trình tích hợp lý thuyết dầm Euler-Bernoulli, cơ học hệ biến đổi khối lượng Meshchersky, nguyên lý biến phân Hamilton cho hệ phi bảo toàn, và động lực học khí động tính toán (CFD). Phạm vi khảo sát tập trung vào mẫu tên lửa quy chuẩn theo thông số đạn phản lực 9M22Y (tổ hợp pháo phản lực BM21 Grad), phóng tại các góc nâng ban đầu $20^\circ$, $25^\circ$ và $30^\circ$, bay trong mặt phẳng thẳng đứng ở dải vận tốc khí động lên tới $V = 500\text{ m/s}$. Kết quả nghiên cứu cung cấp cơ sở khoa học định lượng giúp xác định tiêu chuẩn độ cứng kết cấu tối ưu trong thiết kế sơ bộ vũ khí hàng không vũ trụ.

Literature Review và Positioning

Lịch sử nghiên cứu cơ hệ biến đổi khối lượng bắt nguồn từ thế kỷ 18 với các nghiên cứu của Bernoulli về động cơ phản lực nước, tiếp nối bởi công trình mở đường của George Von Buquoy (1812) và William Moore (1813). Poisson (1838) đã đặt nền móng giải tích khi áp dụng phương trình Lagrange với giả thiết các xung lực vi phân liên tục. Đỉnh cao lý thuyết cổ điển thuộc về I. V. Meshchersky (1897, 1904) với phương trình vi phân chuyển động nổi tiếng $M \frac{dV}{dt} = F + \beta V_1$, xác định lực đẩy phản lực cơ bản $F_{pl} = \beta U$, và công thức Tsiolkovsky về vận tốc lý thuyết tên lửa vũ trụ.

Trong kỷ nguyên hiện đại, bài toán động lực học tên lửa phân hóa thành hai luồng quan điểm đối lập:

  1. Trường phái vật rắn tuyệt đối (Rigid Body Dynamics): Được phát triển sâu rộng bởi Rosser (1947), Gantmacher & Levin (1947), cùng các công trình trong nước (Trần Xuân Diệu et al.; Nguyễn Duy Phồn, 2016; Đặng Võ Công et al., 2018). Luồng nghiên cứu này giả định độ cứng thân tên lửa là vô cùng lớn ($EJ \to \infty$), đơn giản hóa chuyển động về 6 bậc tự do của khối tâm và các góc Euler. Nhược điểm chí mạng là không thể giải thích hiện tượng tản mát điểm nổ và mất ổn định động học của các loại đạn thân mảnh có tỷ lệ chiều dài trên đường kính lớn ($L/D \gg 1$).
  2. Trường phái đàn hồi khí động (Aeroelastic Dynamics): Meirovitch (1970, 1980) tiên phong áp dụng phương pháp thể tích kiểm soát (Control Volume) và mở rộng nguyên lý Hamilton cho dầm đàn hồi biến đổi khối lượng. Tuzcu (2001), Cesnik (2005), và Waszak & Schmidt (1988) tích hợp dao động kết cấu vào động lực học bay máy bay thân mềm. Tuy nhiên, các tranh biện học thuật nảy sinh gay gắt về cách xử lý luồng phụt.

Khi so sánh với các công trình quốc tế tiêu biểu:

  • So sánh với Xiao Huang & Zeiler (2006): Các tác giả mô hình hóa thiết bị bay như một dầm đàn hồi chịu lực đẩy tại đuôi có phương trùng trục dọc dầm, cập nhật khối lượng theo thời gian qua phương trình Lagrange. Tuy nhiên, nghiên cứu này đã bỏ qua tính chất không bảo toàn của cơ hệ hở thoát khí, loại bỏ hiệu ứng góc lệch luồng phụt do biến dạng xoay của tiết diện miệng loa phụt.
  • So sánh với Xu Yun, Xie Changchuan & Yang Chao (2008): Nghiên cứu sự mất ổn định đàn hồi khí động và hiện tượng flutter dưới tác dụng của lực đẩy phản lực bằng lý thuyết Piston bậc 1. Điểm hạn chế là các tác giả chỉ xét dao động ngang thuần túy mà chưa gắn với quỹ đạo bay toàn cục trong trọng trường và chưa tính tới sự dịch chuyển vị trí khối tâm $G(t)$.
  • So sánh với Eke (2012): Eke sử dụng định lý vận chuyển Reynolds cho vật rắn biến đổi khối lượng, khảo sát kích thước và vị trí loa phụt nhưng lại bỏ qua hoàn toàn biến dạng đàn hồi của kết cấu thân.

Luận án của Đàm Việt Phương định vị chính xác tại giao điểm chưa được khai phá: Hợp nhất đồng thời cơ hệ đàn hồi dầm Euler-Bernoulli, hiệu ứng luồng phụt phi bảo toàn thông qua đóng góp động năng mở rộng của thể tích vi phân nhiên liệu thoát qua miệng loa, và sự biến thiên tọa độ khối tâm $X_G(t)$ dọc trục thành một hệ phương trình vi phân phi tuyến hoàn chỉnh.

Đó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 bước tiến quan trọng trong việc hoàn thiện lý thuyết động lực học cơ hệ không bảo toàn biến đổi khối lượng:

  • Mở rộng lý thuyết Meshchersky và nguyên lý Hamilton: Khắc phục hạn chế của việc coi lực đẩy phản lực thuần túy là ngoại lực độc lập. Bằng cách khảo sát hệ mở rộng gồm thân tên lửa và vi phân luồng phụt trong khoảng thời gian vô cùng bé $dt$, năng lượng toàn phần của hệ được phân rã chính xác: động năng thân tên lửa $T_{TL} = T_{TL}^{(1)} + T_{TL}^{(2)} + T_{TL}^{(3)}$ và động năng luồng phụt $T_{LP} = T_{LP}^{(1)} + T_{LP}^{(2)} + T_{LP}^{(3)} + T_{LP}^{(4)}$, trong đó các thành phần chỉ số $(2), (3), (4)$ phản ánh tường minh đóng góp của chuyển vị đàn hồi nút $u, v, w$ và đạo hàm không gian - thời gian của chúng.
  • Xác lập hệ số suy giảm lực đẩy hiệu dụng $k_{hd}$: Luận án chứng minh rằng do biến dạng uốn tại tiết diện miệng loa phụt điểm $P$, véc-tơ vận tốc luồng phụt bị lệch góc so với trục liên kết ban đầu. Vận tốc phụt khí hiệu dụng và lực đẩy hiệu dụng được chuẩn hóa qua công thức giải tích: $$k_{hd} = \frac{1}{\sqrt{1 + \left(\frac{\partial u}{\partial x}\right)P^2 + \left(\frac{\partial v}{\partial x}\right)P^2 + \left(\frac{\partial w}{\partial x}\right)P^2}}$$ $$F{plhd} = \beta U{hd} = k{hd} F_{pl}$$
  • Bổ sung các thành phần lực tương tác vuông góc: Phát hiện các số hạng lực suy rộng mới xuất hiện từ đạo hàm thời gian của ma trận quán tính và ma trận mô-men tĩnh $\dot{M}(t)$, thể hiện tương tác trực tiếp giữa lưu lượng phụt khí $\beta$ và vận tốc dao động uốn $\dot{v}_P$.
[Khung phân tích động lực học kết hợp Đàn hồi - Luồng phụt - Khí động]
+-------------------------------------------------------------------------+
|                  CƠ HỆ ĐÀN HỒI BIẾN ĐỔI KHỐI LƯỢNG                      |
+------------------------------------+------------------------------------+
                                     |
         +---------------------------+---------------------------+
         |                                                       |
         v                                                       v
+----------------------------------+   +----------------------------------+
|  Động lực học thân đàn hồi       |   |  Động lực học luồng phụt         |
|  - Dầm Euler-Bernoulli liên tục  |   |  - Véc-tơ vận tốc: V_lp = V_P + U|
|  - Độ cứng chống uốn EJ(x)       |   |  - Hệ số suy giảm hiệu dụng k_hd |
|  - Mật độ phân bố rho(x, t)      |   |  - Lực đẩy hiệu dụng: F_plhd     |
+----------------------------------+   +----------------------------------+
         |                                                       |
         +---------------------------+---------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|     NGUYÊN LÝ BIẾN PHÂN HAMILTON & PHƯƠNG TRÌNH LAGRANGE HỆ KÍN         |
|     - Thế năng uốn U_e = 1/2 \int EJ_y (v'')^2 dx                       |
|     - Động năng toàn phần: T = T_TL + T_LP                              |
|     - Lực suy rộng khí động CFD Q_q & Trọng trường biến thiên G(t)      |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|  HỆ PHƯƠNG TRÌNH ĐỘNG LỰC HỌC KHÔNG GIAN 3D & MẶT PHẲNG THẲNG ĐỨNG 2D    |
|  [M(t)] {d''} + [C(t, d')] {d'} + [K(t)] {d} = {F_aero} + {F_thrust}   |
+-------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liên ngành 3 trụ cột lý thuyết: (1) Cơ học kết cấu đàn hồi phi tuyến, (2) Động lực học chất khí và luồng phụt phản lực, và (3) Khí động lực học số CFD.

  • Phương pháp tiếp cận rời rạc hóa kết cấu theo phương pháp phần tử hữu hạn (FEM) 1D không gian với ma trận hàm dạng Hermite bậc 3: $$v(x, t) = v_0(t) + v_1(t)x + v_2(t)x^2 + v_3(t)x^3 = [N_i(x)]{d_i(t)}$$
  • Điều kiện biên tự do hai đầu (free-free boundary conditions) được xử lý giải tích thông qua việc thiết lập ma trận liên kết giữa hệ tọa độ tuyệt đối quán tính $OXYZ$ gắn với bệ phóng và hệ tọa độ động $Gxyz$ gắn với khối tâm ban đầu $G$ của tên lửa.
  • Điều kiện biên giới hạn (boundary conditions): Khảo sát biến dạng uốn trong giới hạn đàn hồi tuyến tính của vật liệu kim loại/composite chế tạo vỏ tên lửa, bỏ qua biến dạng xoắn và kéo-nén dọc trục trong mô hình 2D mặt phẳng đứng; không gian khảo sát vận tốc dưới Mach 2.0 trong tầng đối lưu.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo thế giới quan thực chứng hậu hiện đại (Post-positivism / Deterministic Continuum Mechanics), kết hợp chặt chẽ giữa suy diễn giải tích toán học thuần túy (analytical deduction) và mô phỏng số học bán thực nghiệm (computational semi-empirical simulation).

Mô hình nghiên cứu đa tầng được cấu trúc như sau:

  • Tầng 1 (Mô hình giải tích không gian 3D): Thiết lập phương trình vi phân chuyển động toàn phần cho dầm đàn hồi tự do với 6 bậc tự do chuyển động khối và vô hạn bậc tự do biến dạng.
  • Tầng 2 (Rút gọn mô hình trên mặt phẳng bắn 2D): Chiếu hệ phương trình lên mặt phẳng đối xứng thẳng đứng, triệt tiêu chuyển vị ngang $w = 0$, góc hướng $\psi = 0$, góc tự quay $\gamma = 0$, chỉ giữ lại chuyển vị dọc $X_G$, chuyển vị đứng $Y_G$, góc chúc ngóc $\vartheta$, góc tấn $\alpha$, và hàm chuyển vị uốn $v(x, t)$.
  • Tầng 3 (Rời rạc hóa phần tử hữu hạn): Thân tên lửa được chia thành $n-1$ phần tử dầm liên tục bởi $n$ nút. Mỗi nút mạng gồm 2 bậc tự do chuyển vị uốn: chuyển vị thẳng đứng $v_i$ và góc xoay tiết diện $\theta_i = \left(\frac{\partial v}{\partial x}\right)_i$.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu số được thiết kế kiểm soát nghiêm ngặt sai số:

  1. Thiết lập trường khí động học CFD: Dữ liệu trường áp suất $p(x, V, \alpha)$ và trường vận tốc dòng bao quanh thân tên lửa được giải bằng lưới phần tử thể tích khí động. Hệ số lực nâng $C_y(x)$ và lực cản $C_x(x)$ phân bố dọc theo chiều dài không thứ nguyên $x/L$ được trích xuất tại các dải vận tốc từ $V = 100\text{ m/s}$ đến $V = 500\text{ m/s}$.
  2. Thuật toán tích phân phi tuyến từng bước: Hệ phương trình vi phân chuyển động phi tuyến có ma trận khối lượng $[M(t)]$, ma trận cản $[C(t)]$, và ma trận độ cứng $[K(t)]$ biến thiên liên tục theo thời gian được giải số bằng phương pháp tích phân ẩn Newmark-$\beta$ ($\beta = 0.25, \gamma = 0.5$ đảm bảo tính ổn định vô điều kiện), kết hợp vòng lặp hiệu chỉnh phi tuyến Newton-Raphson tại mỗi bước thời gian $\Delta t = 10^{-4}\text{ s}$.
  3. Triệt tiêu bậc tự do phụ thuộc: Sử dụng các điều kiện liên kết động học tại gốc tọa độ $G$ để khử các chuyển vị cứng, thu được hệ phương trình vi phân rút gọn độc lập tuyến tính.
[Sơ đồ thuật toán giải tích hợp Newmark - Newton-Raphson - CFD]
+---------------------------------------------------------------+
| Nhập tham số ban đầu: EJ(x), rho(x, t), F_pl(t), beta(t), theta_0|
+-------------------------------+-------------------------------+
                                |
                                v
+---------------------------------------------------------------+
| Tính toán vị trí khối tâm G(t) và ma trận tĩnh M(t), J_0(t)   |
+-------------------------------+-------------------------------+
                                |
                                v
+---------------------------------------------------------------+
| Lưới CFD: Trích xuất trường khí động phân bố FL(x), FD(x), alpha(x)|
+-------------------------------+-------------------------------+
                                |
                                v
+---------------------------------------------------------------+
| Lập ma trận tổng thể FEM: [M(t)], [C(t)], [K(t)], {Q(t)}      |
+-------------------------------+-------------------------------+
                                |
                                v
+---------------------------------------------------------------+
| Vòng lặp Newmark-beta kết hợp lặp Newton-Raphson              |
| |R(d_{t+dt})| < Tolerance (10^-6)                            |
+-------------------------------+-------------------------------+
                                |
                                v
+---------------------------------------------------------------+
| Xuất nghiệm: Quỹ đạo (X, Y), Vận tốc V, Góc theta, Dao động v(x,t)|
+---------------------------------------------------------------+

Data và phân tích

Đối tượng tính toán mô phỏng sử dụng bộ tham số thực tế của mẫu tên lửa chuẩn 9M22Y:

  • Chiều dài toàn bộ thân đạn $L$, phân chia thành các khoang kết cấu: khoang chiến đấu, khoang động cơ nhiên liệu rắn, khoang đuôi và cánh ổn định.
  • Phân bố mật độ khối lượng ban đầu $\rho(x, 0)$ và quy luật suy giảm khối lượng nhiên liệu theo thời gian đốt $\beta(t) = -\frac{dm}{dt}$.
  • Phân bố độ cứng chống uốn chuẩn định $EJ(x) = K_{gd}$.
  • Quy chuẩn kiểm tra độ tin cậy (Validation protocols):
    • Kiểm chứng thuật toán số: So sánh nghiệm tính toán của hệ phương trình với nghiệm giải tích chính xác của phương trình dầm dao động tự do tuyến tính đã biết nghiệm. Sai số tương đối giữa nghiệm số và giải tích đối với các hàm tọa độ suy rộng $x_1(t), x_2(t), x_3(t)$ nhỏ hơn $0.5%$.
    • Kiểm chứng bán thực nghiệm: Mô phỏng quỹ đạo bay của đạn 9M22Y ở các góc phóng $20^\circ, 25^\circ, 30^\circ$ và đối chiếu trực tiếp với dữ liệu bảng bắn tiêu chuẩn pháo phản lực BM21. Sai lệch về tầm bắn cực đại $X_{max}$, độ cao đỉnh quỹ đạo $Y_{max}$ và vận tốc cuối giai đoạn chủ động nằm trong khoảng sai số cho phép ($< 1.8%$).

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

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

  1. Hiệu ứng suy giảm lực đẩy và phân nhánh quỹ đạo do góc xoay loa phụt: Tại miệng loa phụt điểm $P$, dao động uốn sinh ra góc xoay tức thời $\theta_P = \left(\frac{\partial v}{\partial x}\right)P \neq 0$. Hệ số lực đẩy hiệu dụng $k{hd} < 1.0$ làm giảm trực tiếp thành phần lực đẩy dọc trục và sản sinh thành phần lực đẩy vuông góc gây mất ổn định mô-men chúc ngóc. Tại góc phóng $20^\circ$, quỹ đạo tên lửa đàn hồi có sự sai lệch đáng kể về tầm xa và độ cao so với mô hình vật rắn tuyệt đối.
  2. Quy luật tắt dần của biên độ dao động đàn hồi: Mô phỏng dao động tại mũi và đuôi tên lửa cho thấy: trong giai đoạn động cơ hoạt động (giai đoạn chủ động), biên độ dao động uốn tại mũi đạt giá trị cực đại do lực đẩy phản lực và tải khí động cộng hưởng. Sau khi động cơ tắt thuốc phóng (giai đoạn thụ động), dao động đàn hồi chuyển sang dao động tự do tắt dần nhanh chóng dưới tác dụng của lực cản khí động dập tắt dao động.
  3. Ảnh hưởng phi tuyến của độ cứng chống uốn $EJ$: Khảo sát so sánh giữa các trường hợp độ cứng thực tế $K = K_{gd}$, $K = 50K_{gd}$ và $K = 100K_{gd}$ chứng minh rằng: khi độ cứng tăng lên 100 lần ($K = 100K_{gd}$), các đường cong quỹ đạo, vận tốc bay, góc tấn $\alpha(t)$, và góc chúc ngóc $\vartheta(t)$ tiệm cận hoàn toàn với đường cong tính toán theo mô hình vật rắn tuyệt đối. Điều này khẳng định độ mềm kết cấu là nguồn gốc duy nhất sinh ra sự biến thiên dao động bậc cao quanh trục chuyển động.
  4. Tác động kép của dịch chuyển khối tâm $X_G(t)$: Khi thuốc phóng cháy từ đuôi lên phía trước, vị trí khối tâm $X_G(t)$ dịch chuyển liên tục về phía mũi đạn làm thay đổi cánh tay đòn mô-men khí động và mô-men quán tính $J_0(t)$. Nếu không tính đến dịch chuyển khối tâm trong phương trình Lagrange, sai số tích lũy của góc tấn $\alpha$ có thể vượt quá $12%$, dẫn đến tính toán sai lệch nghiêm trọng lực nâng khí động $F_L$.
  5. Đặc tính lực đẩy động cơ và thời gian cháy: Khảo sát trường hợp giảm cường độ lực đẩy $F_{pl}$ kết hợp kéo dài thời gian hoạt động của động cơ (với cùng tổng xung lượng) cho thấy: tải trọng biến dạng uốn tác dụng lên thân tên lửa giảm rõ rệt, biên độ dao động đàn hồi tại đuôi giảm trên $30%$, giúp ổn định góc chúc ngóc và tăng độ chụm khi tiếp cận mục tiêu.
[Bảng tổng hợp so sánh các trường hợp khảo sát độ cứng và khối tâm]
+--------------------+---------------------+---------------------+---------------------+
| Tham số khảo sát   | Độ cứng thực K_gd   | Độ cứng cao 100K_gd | Vật rắn tuyệt đối   |
+--------------------+---------------------+---------------------+---------------------+
| Biên độ dao động   | Cực đại tại mũi/đuôi| Xấp xỉ 0            | Bằng 0              |
| Suy giảm lực đẩy   | k_hd < 1.0 (dao động)| k_hd = 1.0 (hằng số)| Không xét k_hd      |
| Dao động góc tấn   | Dao động tần số cao | Biến thiên trơn     | Biến thiên trơn     |
| Tầm bắn & Độ cao   | Khớp chuẩn bảng BM21| Lệch thực tế        | Lệch thực tế        |
+--------------------+---------------------+---------------------+---------------------+

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện giải pháp giải tích cho hệ phương trình vi phân chuyển động phi tuyến của dầm Euler-Bernoulli liên kết hở mang khối lượng biến đổi, lấp đầy khiếm khuyết trong các giáo trình động lực học tên lửa truyền thống.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp FEM 1D phi tuyến với cơ sở dữ liệu khí động CFD phân bố, mở ra hướng ứng dụng hiệu quả cho việc mô phỏng các kết cấu hàng không vũ trụ phức tạp mà không đòi hỏi chi phí tính toán 3D FEM toàn phần quá đắt đỏ.
  • Về mặt thực tiễn kỹ thuật: Cung cấp công cụ tính toán định lượng cho các kỹ sư thiết kế vũ khí; xác định ngưỡng độ cứng chống uốn giới hạn $EJ_{min}$ nhằm triệt tiêu hiện tượng cộng hưởng dao động uốn thân với hệ thống cánh lái; tối ưu hóa quy luật cháy của buồng đốt động cơ phản lực.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật đột phá, công trình vẫn tồn tại một số giới hạn nghiên cứu:

  1. Giả thiết dầm Euler-Bernoulli: Luận án bỏ qua biến dạng trượt cắt (lý thuyết dầm Timoshenko) và biến dạng xoắn dọc trục. Đối với các loại tên lửa có cánh lái quay tự do hoặc tên lửa có tỷ lệ $L/D$ trung bình, biến dạng cắt và xoắn có thể đóng góp một phần vào sai số góc tấn cục bộ.
  2. Giới hạn mặt phẳng bắn 2D: Việc rút gọn hệ phương trình không gian 3D về mặt phẳng thẳng đứng 2D tuy phù hợp với đạn phản lực không điều khiển nhưng chưa mô tả trọn vẹn hiện tượng ghép kênh dao động uốn - dạt sườn (cross-coupling) khi tên lửa quay quanh trục đối xứng.
  3. Bỏ qua gradient nhiệt độ: Chưa tích hợp ảnh hưởng của nhiệt độ ma sát khí động cực cao và nhiệt độ buồng đốt làm suy giảm mô-đun đàn hồi $E(T)$ của vật liệu thân theo thời gian bay.

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

  • Nâng cấp mô hình kết cấu sang dầm Timoshenko phi tuyến 3D có xét đến tương tác nhiệt - đàn hồi khí động học (aerothermoelasticity).
  • Mở rộng thuật toán cho tên lửa có điều khiển (guided missile) tích hợp vòng lặp phản hồi của hệ thống lái tự động (Autopilot), khảo sát hiện tượng mất ổn định do cảm biến con quay quán tính (IMU) đo nhầm tín hiệu dao động uốn thân dầm.
  • Nghiên cứu cơ cấu điều khiển véc-tơ lực đẩy (TVC) chủ động nhằm dập tắt dao động đàn hồi thân tên lửa trong thời gian thực.

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

  • Ảnh hưởng học thuật (Academic Impact): Tạo tiền đề cho các hướng nghiên cứu chuyên sâu tại các viện nghiên cứu, trường đại học kỹ thuật quân sự; cung cấp khung tham chiếu chuẩn cho các nghiên cứu sinh chuyên ngành Cơ kỹ thuật, Cơ học vật rắn và Kỹ thuật Hàng không Vũ trụ. Ước tính các công bố từ luận án mở ra hướng trích dẫn quan trọng trong mảng động lực học bay đàn hồi.
  • Hiện đại hóa công nghiệp quốc phòng: Đóng góp trực tiếp vào các chương trình thiết kế, chế tạo và cải tiến các dòng đạn phản lực, tên lửa phòng không, tên lửa hành trình nội địa của Việt Nam; giảm thiểu số lượng mẫu thử nghiệm bắn đạn thật tốn kém thông qua việc mô phỏng số chính xác cao.
  • Tối ưu hóa kinh tế - kỹ thuật: Việc xác định chính xác độ cứng chống uốn yêu cầu cho phép loại bỏ kết cấu dư thừa trọng lượng, tối ưu hóa tỷ lệ khối lượng nhiên liệu trên khối lượng kết cấu thân đạn, giúp gia tăng tầm bắn chiến thuật từ $5%$ đến $10%$ mà không cần tăng kích thước động cơ.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận phương pháp luận biến phân Hamilton mở rộng, kỹ thuật lập trình số Newmark kết hợp Newton-Raphson xử lý ma trận khối lượng biến thiên liên tục.
  • Các nhà khoa học đầu ngành Hàng không vũ trụ: Kế thừa khung mô hình toán học 3D hoàn chỉnh để phát triển các bài toán phức tạp hơn về phương tiện bay siêu vượt âm (hypersonic vehicles).
  • Kỹ sư thiết kế R&D tại các nhà máy/viện nghiên cứu quân giới: Ứng dụng trực tiếp phần mềm tính toán vào quy trình thiết kế sơ bộ kết cấu vỏ đạn, cánh lái và tối ưu biên dạng buồng đốt động cơ tên lửa nhiên liệu rắn.
  • Cơ quan quản lý và hoạch định quốc phòng: Có thêm cơ sở khoa học tin cậy để thẩm định, đánh giá chất lượng và tuổi thọ kỹ thuật của các loại khí tài tên lửa đang có trong biên chế.

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 hoàn thiện nguyên lý biến phân Hamilton cho hệ mở rộng (thân dầm đàn hồi + vi phân thể tích luồng phụt), mở rộng trực tiếp cơ học hệ biến đổi khối lượng của Meshchersky và lý thuyết thể tích kiểm soát của Meirovitch. Luận án đã chứng minh sự tồn tại của hệ số suy giảm lực đẩy hiệu dụng $k_{hd}$ và trích xuất tường minh các thành phần lực tương tác vuông góc với trục thân do tích phân vi phân lưu lượng phụt $\beta$ với vận tốc dao động uốn $\dot{v}_P$.

2. Điểm cải tiến phương pháp luận so với ít nhất 2 công trình quốc tế tiền nhiệm?

So với Xiao Huang & Zeiler (2006) – nhóm tác giả chỉ áp dụng phương trình Lagrange quy ước và coi lực đẩy là ngoại lực độc lập dọc trục dầm – luận án đã mô hình hóa chính xác tính phi bảo toàn của hệ hở mang luồng phụt. So với Xu Yun et al. (2008) – nhóm chỉ giải bài toán ổn định dao động ngang bằng lý thuyết Piston mà không gắn với bài toán quỹ đạo toàn cục – luận án đã tích hợp đồng thời trường khí động CFD phân bố, sự dịch chuyển khối tâm thực tế $X_G(t)$, và quỹ đạo bay phi tuyến trong trọng trường.

3. Phát hiện thực nghiệm/mô phỏng gây bất ngờ nhất có số liệu minh chứng?

Phát hiện rằng dao động uốn thân đạn tác động phi đối xứng lên góc chúc ngóc $\vartheta$ và góc tấn $\alpha$ trong giai đoạn động cơ hoạt động, làm suy giảm đáng kể lực đẩy hiệu dụng $F_{plhd}$ qua hệ số $k_{hd}$. Khi tăng độ cứng chống uốn kết cấu lên $K = 100K_{gd}$, toàn bộ dao động tần số cao bị dập tắt, quỹ đạo chuyển dịch hoàn toàn về nghiệm vật rắn tuyệt đối. Tuy nhiên, nghiệm vật rắn tuyệt đối lại có sai số lớn hơn so với bảng bắn thực tế BM21, khẳng định rằng tính đàn hồi thực tế ($K = K_{gd}$) là yếu tố nội tại giải thích độ chụm và tản mát của đạn thật.

4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp rõ ràng không?

Luận án cung cấp đầy đủ hệ phương trình vi phân dạng ma trận tổng thể, ma trận hàm dạng Hermite bậc 3 $[N_i(x)]$, thuật toán ghép nối ma trận phần tử FEM, quy trình lặp Newmark-$\beta$ ($\beta = 0.25, \gamma = 0.5$) kết hợp Newton-Raphson, và bảng dữ liệu hình học, phân bố khối lượng, độ cứng $EJ$ của đạn mẫu 9M22Y, cho phép tái lập hoàn toàn trên các nền tảng tính toán khoa học như MATLAB hoặc C++.

5. Định hướng chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Định hướng phát triển tập trung vào 3 giai đoạn: (1) Tích hợp hiệu ứng nhiệt - đàn hồi khí động 3D cho vật liệu composite vỏ tên lửa; (2) Thiết lập vòng điều khiển thích nghi kháng dao động uốn (Aeroelastic-filtering Autopilot) cho tên lửa có điều khiển; (3) Thử nghiệm đo đạc biến dạng động lực học thực tế bằng cảm biến sợi quang FBG gắn trên thân đạn trong các đợt bắn thử nghiệm trường bắn.

Kết luận

  1. Thiết lập thành công mô hình toán học không gian 3 chiều mô tả chuyển động của tên lửa đàn hồi biến đổi khối lượng, tích hợp đầy đủ các yếu tố dao động uốn, dịch chuyển khối tâm dọc trục $X_G(t)$, và trường khí động phân bố CFD.
  2. Phát hiện và chứng minh giải tích hệ số suy giảm lực đẩy hiệu dụng $k_{hd}$ cùng các thành phần lực liên kết động học vuông góc tại tiết diện miệng loa phụt, giải quyết triệt để khiếm khuyết của các mô hình Lagrange quy ước.
  3. Xây dựng thuật toán số FEM kết hợp Newmark-$\beta$ và Newton-Raphson đạt độ chính xác cao, kiểm chứng thành công với nghiệm giải tích dầm tự do ($< 0.5%$ sai số) và bảng bắn thực tế pháo phản lực BM21 ($< 1.8%$ sai số).
  4. Làm sáng tỏ định lượng ảnh hưởng của độ cứng chống uốn $EJ$: Chứng minh sự cần thiết phải kiểm soát dao động đàn hồi để đảm bảo độ chính xác quỹ đạo, đồng thời cung cấp phương pháp xác định độ cứng tối ưu nhằm giảm trọng lượng kết cấu.
  5. Mở ra 3 nhánh nghiên cứu học thuật mới: Khí đàn hồi nhiệt phi tuyến cho phương tiện bay cao tốc; tương tác đàn hồi - điều khiển tự động (Aeroservoelasticity); và tối ưu hóa đa mục tiêu kết cấu - động cơ phản lực nhiên liệu rắn.