Tổng quan luận án
Tính cấp thiết và khoảng trống nghiên cứu
Trong lĩnh vực khai thác và điều khiển phương tiện bay, chất lượng bay phụ thuộc đồng thời vào đặc tính khí động của máy bay, tác động điều khiển của người lái và các yếu tố nhiễu động từ môi trường bên ngoài như nhiệt độ, áp suất, nhiễu động không khí. Trong các yếu tố môi trường, gió cạnh là một dạng nhiễu làm lệch dòng không khí, biến đổi trường dòng chảy bao đối xứng quanh máy bay thành dòng chảy phi đối xứng. Sự xuất hiện của gió cạnh làm phân bố lại áp suất trên bề mặt khí động và làm mất tính đối xứng trong các phép tính toán đặc tính khí động (ĐTKĐ).
Quá trình hạ cánh bao gồm các giai đoạn bay diễn ra ở độ cao thấp và vùng vận tốc nhỏ. Đây là thời điểm các đặc tính khí động của máy bay chịu tác động rõ rệt từ gió cạnh – một hiện tượng ngẫu nhiên, khó dự báo và tiềm ẩn nguy cơ uy hiếp an toàn bay. Theo thống kê an toàn bay được tác giả trích dẫn, có tới khoảng 85% các sự cố và tai nạn bay xảy ra có liên quan đến tác động của gió cạnh trong quá trình hạ cánh. Tuy nhiên, tại Việt Nam, các nghiên cứu trước đây chưa dành nhiều sự quan tâm đến ảnh hưởng của gió cạnh tới ĐTKĐ của máy bay và chưa đánh giá đầy đủ gió cạnh như một nguyên nhân kỹ thuật trực tiếp gây mất an toàn. Do đó, đề tài "Nghiên cứu ảnh hưởng của gió cạnh tới đặc tính khí động của máy bay trong quá trình hạ cánh" được thực hiện nhằm giải quyết khoảng trống nghiên cứu này, phục vụ yêu cầu đảm bảo an toàn bay, thiết kế chế tạo và khai thác sử dụng máy bay.
Mục tiêu nghiên cứu
Mục đích tổng quát của luận án là xem xét mức độ ảnh hưởng của gió cạnh đến các ĐTKĐ của máy bay trong quá trình hạ cánh ở vùng tốc độ nhỏ. Để đạt được mục đích đó, luận án xác định 3 nội dung nhiệm vụ cụ thể:
- Xây dựng mô hình toán máy bay theo sơ đồ không gian với cấu hình cất, hạ cánh trong dòng chảy bao dừng, phi đối xứng.
- Xây dựng phương pháp nghiên cứu xác định mức độ ảnh hưởng của gió cạnh đến ĐTKĐ của máy bay trong quá trình hạ cánh ở vùng tốc độ nhỏ.
- Đề ra một số giải pháp làm giảm ảnh hưởng của gió cạnh đến ĐTKĐ của máy bay cũng như nâng cao chất lượng bay, đảm bảo an toàn bay.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Máy bay có phần tạo lực nâng cố định, được cơ giới hóa (hệ thống cánh tà) và trang bị các thiết bị hạ cánh (hệ thống càng) để hoạt động trên đường băng. Đối tượng áp dụng tính toán và kiểm chứng cụ thể là máy bay huấn luyện hai buồng lái L-39.
- Phạm vi nghiên cứu: Khảo sát ảnh hưởng của gió cạnh đến một số đặc tính khí động học cơ bản (thuộc chuyên ngành hẹp khí động học) trong quá trình hạ cánh ở vùng vận tốc nhỏ, trọng tâm là khả năng tạo lực nâng của máy bay trong điều kiện dòng chảy bao dừng, phi tuyến.
- Giả thiết tính toán: Gió cạnh được giả thiết thổi trong mặt phẳng nằm ngang với góc thổi $\delta = 90^\circ$ và có cường độ không đổi.
Thông tin tổng quan luận án
| Tiêu chí | Thông tin chi tiết |
|---|---|
| Tên luận án | Nghiên cứu ảnh hưởng của gió cạnh tới đặc tính khí động của máy bay trong quá trình hạ cánh |
| Nghiên cứu sinh | Nguyễn Đình Sơn |
| Người hướng dẫn khoa học | 1. PGS. TS Phạm Vũ Uy 2. GS. TSKH Nguyễn Đức Cương |
| Cơ sở đào tạo | Viện Khoa học và Công nghệ Quân sự (Bộ Quốc phòng) |
| Chuyên ngành | Kỹ thuật cơ khí động lực |
| Mã số chuyên ngành | 62 52 01 16 |
| Năm bảo vệ | 2014 |
| Quy mô tư liệu | 100 trang nội dung chính, 9 bảng số liệu, 48 hình vẽ và đồ thị, danh mục tài liệu tham khảo và các phụ lục |
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu quốc tế
Nghiên cứu khí động lực học của dòng chảy bao ở vận tốc nhỏ đã được triển khai từ thế kỷ 20 với các mốc công trình tiêu biểu:
- Nghiên cứu thực nghiệm vận tốc nhỏ ban đầu: Năm 1937, A. Kovalev (Nga) tiến hành các thí nghiệm trong ống khí động (OKĐ) và ghi nhận hiện tượng lực cản tăng nhẹ ở vùng vận tốc nhỏ. Nhà khoa học F. Schmitz (Đức) nghiên cứu sự ảnh hưởng của độ rối dòng chảy trong ống thổi đến các hệ số khí động qua các thí nghiệm mô hình ở dải vận tốc giảm dần từ 30–50 m/s xuống 4–12 m/s, chỉ ra sự phụ thuộc của phép đo lực vào đặc tính dòng, thiết bị đo và nhiệt độ môi trường. Các nghiên cứu lý thuyết và thực nghiệm giai đoạn 1939–1970 tập trung vào việc lựa chọn profil cánh cho quạt gió và mô hình ở số Reynolds ($Re$) nhỏ.
- Nghiên cứu dòng dưới âm và vượt âm: Từ cuối thập niên 1970, nhiều công trình nghiên cứu dòng chảy bao dưới âm cho dòng nén được và không nén được ở vùng số $Re$ lớn (bỏ qua lực nhớt động). Tác giả nhận định các tài liệu công bố hầu như thiếu vắng các kết quả và đồ thị khí động ở dải số Mach $M \in (0 \div 0,2)$, nơi mà tác động của lực nhớt động không thể bỏ qua.
- Nghiên cứu về tác động của gió cạnh: Các công trình nghiên cứu gió cạnh tập trung vào công trình xây dựng, phương tiện giao thông mặt đất và khí cụ bay. Về mặt khí động cánh và thân:
- S. Belotserkovsky (Белоцерковский C.): Ứng dụng lý thuyết lực nâng Zhukovsky để xây dựng phương pháp xoáy rời rạc (XRR) tính toán ĐTKĐ cho cánh nâng mỏng trong dòng dừng và không dừng (giả thiết khí lý tưởng, bỏ qua tính nén và độ nhớt).
- G. Chernov (Чернов Г.) và các tác giả: Khảo sát dòng chảy bao thân tròn, cấu trúc xoáy phụ thuộc hình dạng thân, cường độ xoáy và tương tác khí động thân - cánh có tính đến độ nhớt ở các số $Re$ khác nhau.
- Will E. Graf, Yann Colin: Khảo sát ảnh hưởng của gió cạnh đến dòng chảy trước cửa thu khí động cơ và cánh quạt đuôi trực thăng, chỉ ra sự suy giảm hiệu suất và phát sinh mô men lệch.
- Yoshiharu Tanahashi (MHI - HUST): Công bố kết quả chỉ ra ở vùng cận âm, khi thử nghiệm OKĐ chỉ đảm bảo đồng dạng số $M$ mà không đảm bảo số $Re$, hệ số lực nâng cực đại $C_{y\max}$ của mô hình thu được nhỏ hơn giá trị khi bay thử thực tế, sự sai lệch càng lớn khi $Re$ càng nhỏ.
Phần lớn các nghiên cứu quốc tế đã công bố chủ yếu giải quyết bài toán động lực học bay (thay đổi quỹ đạo, kỹ thuật lái khi có gió cạnh) mà chưa đi sâu vào việc xác định sự suy giảm ĐTKĐ tổng thể của máy bay ở cấu hình hạ cánh.
Tình hình nghiên cứu trong nước
Tại Việt Nam, một số công trình nghiên cứu đã bước đầu tiếp cận phương pháp XRR phi tuyến để tính toán khí động học:
- Các nghiên cứu đã mô phỏng cấu trúc vết xoáy sau máy bay trong điều kiện dòng chảy không đối xứng (khi góc tấn $\alpha = 10^\circ$, góc trượt cạnh $\beta = 10^\circ$).
- Xây dựng chương trình tính toán cho máy bay có cấu hình không gian phức tạp bằng cách tổ hợp các mặt nâng cơ sở mỏng.
Tuy nhiên, các nghiên cứu trong nước trước đó chủ yếu tập trung vào dòng chảy bao đối xứng ở cấu hình bay bằng (bay thẳng), chưa mô hình hóa đầy đủ cấu hình cất/hạ cánh (khi thả cánh tà và thả càng) và chưa giải quyết bài toán ĐTKĐ khi có gió cạnh thổi ngang.
Khoảng trống luận án lựa chọn giải quyết
Luận án tập trung giải quyết bài toán xác định các ĐTKĐ chủ yếu ($C_y, C_x, C_z, m_x, m_y, m_z$) của máy bay L-39 ở cấu hình cất/hạ cánh hoàn chỉnh trong vùng vận tốc nhỏ có tính đến góc trượt cạnh $\beta$ và góc tấn $\alpha$. Đồng thời, luận án xây dựng phương pháp luận khắc phục khó khăn về kiểm chứng thực nghiệm bằng việc kết hợp đồng thời hai phương pháp số độc lập.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết khí động học
- Đặc trưng dòng chảy vận tốc nhỏ: Khi máy bay chuyển động ở vận tốc nhỏ, thành phần biến thiên đối lưu của vận tốc giảm, thành phần lực nhớt động (khuếch tán) trở nên đáng kể và không thể triệt tiêu. Thực nghiệm chỉ ra ở vận tốc dưới 12 m/s, dòng chảy bao quanh mô hình bắt đầu mất ổn định và hiện tượng mất ổn định gia tăng khi dòng chảy phi đối xứng (xuất hiện góc trượt cạnh $\beta$).
- Lý thuyết lực nâng Zhukovsky và định luật Biot-Savart: Sử dụng để xác định trường vận tốc cảm ứng và phân bố tải khí động học dạng không thứ nguyên từ mạng lưới xoáy liên kết và xoáy tự do.
- Tích phân Cauchy-Lagrange: Dùng để tính toán trường áp suất và lực khí động tác dụng lên các mặt phẳng nâng mỏng.
- Hệ phương trình Navier-Stokes: Mô tả chuyển động của chất lưu nhớt thực, giải gần đúng bằng phương pháp thể tích hữu hạn trên phần mềm Ansys, kết hợp các mô hình chảy rối ($k$-$\varepsilon$, $k$-$\Omega$) và thuật toán liên kết áp suất - vận tốc (SIMPLE, SIMPLEF, SIMPLEN).
Phân tích điều kiện đồng dạng và hạn chế của ống khí động
Luận án làm rõ lý do không sử dụng thực nghiệm OKĐ truyền thống trong điều kiện Việt Nam thông qua bài toán tính toán điều kiện đồng dạng đồng thời về số Mach ($M = V/a$) và số Reynolds ($Re = VL/\nu$):
$$\text{Xét tỷ số: } S = \frac{a}{\nu}$$
Giả thiết bài toán tính toán so sánh:
- Máy bay thực khi bay thử (BT): Sải cánh $L_{BT} = 5\text{ m}$, vận tốc $V_{BT} = 50\text{ m/s}$, nhiệt độ $t_{BT} = 35^\circ\text{C}$, áp suất $p = 101.325\text{ Pa}$, vận tốc âm thanh $a_{BT} = 352\text{ m/s}$, $M_{BT} = 0,142$, tỷ số $S_{BT} = 2,137 \times 10^5\text{ m}^{-1}$.
- Mô hình thử nghiệm trong OKĐ (TN): Sải cánh $L_{TN} = 2\text{ m}$ (tỷ lệ hình học $k = L_{BT}/L_{TN} = 2,5$).
- Để thỏa mãn đồng thời $M_{TN} = M_{BT} = 0,142$ và $Re_{TN} = Re_{BT}$, giá trị $S_{TN}$ trong ống thổi phải đạt:
$$S_{TN} = k \cdot S_{BT} = 2,5 \times 2,137 \times 10^5 = 5,343 \times 10^5\text{ m}^{-1}$$
Để đạt được giá trị này ở áp suất khí quyển, môi chất trong OKĐ bắt buộc phải được làm lạnh sâu xuống nhiệt độ:
$$t_{TN} = -116,5^\circ\text{C}$$
Khi đó, vận tốc âm thanh tương ứng là $a_{TN} = 250,95\text{ m/s}$ và vận tốc dòng thổi là $V_{TN} = 35,63\text{ m/s}$. Nếu thử nghiệm ở nhiệt độ môi trường bình thường ($35^\circ\text{C}$), để giữ nguyên $Re$ thì số Mach trong OKĐ phải tăng lên $M_{TN} = 0,355$ ($V_{TN} = 124,96\text{ m/s}$), làm phá vỡ tính đồng dạng về khả năng nén của dòng khí. Do các OKĐ tại Việt Nam chưa có hệ thống làm lạnh sâu chuyên dụng, phương pháp thực nghiệm OKĐ không đảm bảo độ tin cậy khoa học cho bài toán này.
Phương pháp nghiên cứu kết hợp hai công cụ số
Để khắc phục trở ngại về thực nghiệm, luận án áp dụng giải pháp "phòng thí nghiệm số" với sự phối hợp của hai phương pháp số độc lập:
- Phương pháp Xoáy rời rạc (XRR) phi tuyến: Thuật toán tính toán nhanh, cấu trúc toán học tường minh, phù hợp cho việc khảo sát định hướng và xác định miền hội tụ nghiệm ban đầu.
- Phần mềm tính toán động lực học dòng chảy Ansys (CFD): Giải hệ phương trình Navier-Stokes đầy đủ trên lưới không gian 3 chiều có bản quyền, mô phỏng chính xác các hiệu ứng nhớt, hình học vỏ thân, cánh tà và càng hạ cánh.
Độ tin cậy của kết quả nghiên cứu được xác lập khi sai lệch tính toán giữa hai phương pháp số độc lập nằm trong ngưỡng cho phép ($\Delta \le 15%$), đồng thời các kết quả ở chế độ bay cơ sở được đối chiếu kiểm chứng trực tiếp với số liệu trong Tài liệu kỹ thuật (TLKT) của máy bay L-39.
Nội dung chính theo từng chương
Chương 1. Tổng quan về ảnh hưởng của gió cạnh đến đặc tính khí động của máy bay, các phương pháp nghiên cứu
Chương 1 hệ thống hóa cơ sở lý luận về dòng chảy bao quanh máy bay ở vùng vận tốc nhỏ khi có tác động của gió cạnh:
- Làm rõ bản chất vật lý của hiện tượng gió cạnh: gây mất tính đối xứng của trường áp suất, thúc đẩy quá trình chuyển tiếp lớp biên từ chảy tầng sang chảy rối và làm xuất hiện hiện tượng tách dòng sớm trên bề mặt cánh nâng.
- Phân tích ưu và nhược điểm của các phương pháp xác định ĐTKĐ: phương pháp giải tích (bị giới hạn đối với các kết cấu hình học phức tạp), phương pháp thực nghiệm (vướng rào cản kỹ thuật về điều kiện đồng dạng $M$ và $Re$ trong ống khí động lạnh) và phương pháp số (phương pháp panel, XRR, CFD).
- Xác lập cơ sở phương pháp luận: lựa chọn phương pháp XRR làm công cụ định hướng và phần mềm Ansys làm công cụ tính toán chính xác, kết hợp đối soát với tài liệu kỹ thuật của máy bay L-39.
Chương 2. Mô hình bài toán xác định đặc tính khí động của máy bay
Chương 2 tập trung xây dựng mô hình toán học và thiết lập quy trình tính toán cho máy bay ở chế độ bay thẳng (chưa xét trượt cạnh):
- Mô hình XRR cánh phẳng đơn: Bề mặt cánh được chia thành $N$ dải xoáy dọc và $n$ đoạn xoáy ngang, tạo thành lưới gồm $N \times n$ phần tử tứ giác. Toàn bộ hệ thống nâng được thay thế bằng hệ xoáy liên kết ngang ($\Gamma_{kk1}$), xoáy liên kết dọc ($\Gamma_{k1k}$) và các hệ xoáy tự do vùng I (mép sau), vùng II (mép phải) thỏa mãn điều kiện Chaplygin-Zhukovsky. Thiết lập hệ phương trình đại số tuyến tính $N \times n$ ẩn số dựa trên điều kiện không chảy thấu qua bề mặt:
$$\sum_{k=1}^N \sum_{k_1=1}^n \Gamma_{kk_1} a_{kk_1}^{pp_1} = -2\pi \sin\alpha$$
- Mô hình XRR cấu hình không gian phức tạp: Thân máy bay được mô hình hóa bằng hai mặt nâng phẳng vuông góc (mặt phẳng ngang và mặt phẳng đứng); cánh và cánh đuôi ngang nằm trong mặt phẳng ngang; cánh đuôi đứng nằm trong mặt phẳng đứng. Tải khí động trên từng phần tử được xác định theo phương pháp tuyến và quy đổi về hệ tọa độ gốc $Oxyz$.
- Mô hình số trên Ansys: Mô hình hóa hình học thân, profil cánh và đuôi máy bay L-39; thiết lập miền tính toán và chia lưới không gian; áp dụng phương pháp thể tích hữu hạn với thuật toán SIMPLE để giải trường dòng dừng.
- Đánh giá và kiểm chứng chế độ bay thẳng: So sánh kết quả tính toán hệ số lực nâng $C_y = f(\alpha)$ giữa XRR, Ansys và TLKT máy bay L-39.
Bảng so sánh hệ số lực nâng Cy của máy bay L-39 ở chế độ bay thẳng (trích xuất từ Bảng 2.1 và 2.2):
| Góc tấn $\alpha$ (độ) | $C_y$ (XRR) | $C_y$ (Ansys) | $C_y$ (TLKT) | Chênh lệch XRR vs TLKT (%) | Chênh lệch Ansys vs TLKT (%) |
|---|---|---|---|---|---|
| $0^\circ$ | 0,082 | 0,085 | 0,080 | +2,50% | +6,25% |
| $2^\circ$ | 0,235 | 0,241 | 0,230 | +2,17% | +4,78% |
| $4^\circ$ | 0,388 | 0,396 | 0,380 | +2,11% | +4,21% |
| $6^\circ$ | 0,540 | 0,548 | 0,530 | +1,89% | +3,40% |
| $8^\circ$ | 0,690 | 0,698 | 0,680 | +1,47% | +2,65% |
| $10^\circ$ | 0,838 | 0,845 | 0,830 | +0,96% | +1,81% |
Sai số tương đối giữa kết quả tính toán của hai phương pháp số so với TLKT đều nhỏ hơn $7%$, khẳng định tính đúng đắn và sự hội tụ của mô hình.
Chương 3. Xây dựng mô hình toán xác định đặc tính khí động của máy bay với cấu hình cất, hạ cánh
Chương 3 phát triển mô hình toán học để phản ánh trạng thái hạ cánh thực tế:
- Mô hình XRR cánh nâng cơ sở có trượt cạnh: Xây dựng mô hình không gian bài toán chảy bao khi có góc trượt cạnh $\beta \ne 0$. Hệ thống xoáy tự do được bổ sung thêm vùng II' tại mép trái của cánh nâng nhằm phản ánh sự bất đối xứng của dòng bao. Thiết lập giải thuật vector xác định vector pháp tuyến $\vec{n}_m$ của tứ giác khảo sát $PQRT$ và vector chỉ phương pháp tuyến của các đoạn xoáy liên kết trong không gian 3 chiều.
- Mô hình hóa cơ cấu hạ cánh:
- Cánh tà ở vị trí thả: Sử dụng nguyên tắc ghép nối không đồng phẳng để mô hình hóa cánh tà lệch góc $25^\circ$ (chế độ cất cánh/chuẩn bị hạ cánh) và $44^\circ$ (chế độ hạ cánh hoàn toàn).
- Càng hạ cánh ở vị trí thả: Mô hình hóa hình học càng và tích hợp vào hệ thống mặt nâng cản.
- Kết quả mô phỏng cấu hình hạ cánh ($\beta = 0^\circ$): Khảo sát đường đặc tính lực nâng $C_y = f(\alpha)$ khi thả cánh tà $25^\circ$ và $44^\circ$ bằng cả hai phương pháp XRR và Ansys. Kết quả cho thấy khi thả cánh tà $44^\circ$, hệ số lực nâng tăng đáng kể trên toàn bộ dải góc tấn tuyến tính, độ sai lệch giữa XRR và Ansys so với số liệu kỹ thuật duy trì dưới ngưỡng $10%$.
Chương 4. Nghiên cứu ảnh hưởng của gió cạnh đến đặc tính khí động của máy bay trong quá trình hạ cánh
Chương 4 trình bày các kết quả phân tích định lượng về tác động của gió cạnh tới ĐTKĐ và động lực học hạ cánh:
- Sự thay đổi các hệ số khí động theo góc trượt cạnh $\beta$: Khảo sát trên máy bay L-39 ở cấu hình càng thả, cánh tà thả $44^\circ$, góc tấn $\alpha = 10^\circ$, góc trượt cạnh $\beta$ biến thiên từ $0^\circ$ đến $20^\circ$:
- Hệ số lực nâng $C_y$ suy giảm khi $\beta$ tăng lên.
- Hệ số lực cản $C_x$ và hệ số lực cạnh $C_z$ tăng mạnh.
- Xuất hiện mô men nghiêng cánh $m_x$ và mô men hướng $m_y$ gây mất ổn định không gian.
- Vị trí tiêu điểm khí động theo trục dọc $x_F$ và trục ngang $z_F$ bị dịch chuyển đáng kể theo sự thay đổi của $\beta$.
Bảng biến thiên hệ số khí động của máy bay L-39 theo góc trượt cạnh (α = 10 độ, cánh tà 44 độ, trích xuất Chương 4):
| Góc trượt cạnh $\beta$ | Hệ số lực nâng $C_y$ | Hệ số lực cản $C_x$ | Hệ số lực cạnh $C_z$ | Hệ số mô men nghiêng $m_x$ | Hệ số mô men hướng $m_y$ |
|---|---|---|---|---|---|
| $0^\circ$ | 1,420 | 0,215 | 0,000 | 0,000 | 0,000 |
| $5^\circ$ | 1,405 | 0,224 | -0,065 | 0,018 | -0,012 |
| $10^\circ$ | 1,368 | 0,248 | -0,132 | 0,038 | -0,026 |
| $15^\circ$ | 1,302 | 0,286 | -0,205 | 0,062 | -0,043 |
| $20^\circ$ | 1,210 | 0,339 | -0,285 | 0,091 | -0,063 |
- Mối liên hệ cân bằng giữa các góc lái điều khiển với góc trượt cạnh: Thiết lập tương quan giữa góc lệch cánh lái hướng $\delta_H$ và cánh lái liệng $\delta_L$ cần thiết để triệt tiêu mô men xoay $M_y(\beta)$ và mô men nghiêng $M_x(\beta)$ do gió cạnh sinh ra.
- Xác định góc trượt cạnh tới hạn $\beta^*$ và miền vận tốc hạ cánh an toàn: Luận án chỉ ra rằng khi vận tốc bay $V$ giảm, góc trượt cạnh tới hạn $\beta^$ cho phép cũng giảm theo. Vượt quá ngưỡng $\beta^$, lực nâng bị suy giảm đột ngột và hệ thống điều khiển lái không còn đủ hiệu lực để giữ cân bằng máy bay.
Bảng góc trượt cạnh tới hạn β* theo vận tốc hạ cánh của máy bay L-39 (trích xuất Bảng 4.6):
| Vận tốc hạ cánh $V$ (km/h) | Vận tốc hạ cánh $V$ (m/s) | Góc trượt cạnh tới hạn $\beta^*$ (độ) | Trạng thái an toàn |
|---|---|---|---|
| 160 | 44,4 | $10,5^\circ$ | Giới hạn nguy hiểm cao |
| 170 | 47,2 | $12,0^\circ$ | Miền cảnh báo |
| 180 | 50,0 | $13,8^\circ$ | Miền an toàn mở rộng |
| 190 | 52,8 | $15,5^\circ$ | Miền an toàn định mức |
| 200 | 55,6 | $17,2^\circ$ | Miền an toàn cao |
- Khảo sát quỹ đạo và biện pháp xử lý: Phân tích quỹ đạo hạ cánh có trượt cạnh qua 2 giai đoạn: bay xuống ($h = 20,6\text{ m}$) và cải bằng ($h = 3\text{ m}$). Đề xuất các kỹ thuật phối hợp lái nghiêng cánh và bẻ lái hướng để duy trì tiếp đất dọc theo trục tâm đường băng.
Kết quả và những đóng góp mới
Những đóng góp mới của luận án
Theo công bố của tác giả, luận án mang lại 4 đóng góp mới về mặt khoa học và kỹ thuật:
- Phát triển mô hình tính toán khí động: Mở rộng mô hình không gian cho cánh nâng cơ sở bằng phương pháp xoáy rời rạc phi tuyến, cho phép tính toán dòng chảy bao dừng phi đối xứng có trượt cạnh ($\beta \ne 0$) và ghép nối không đồng phẳng với hệ thống cánh tà hạ và càng hạ cánh.
- Phát hiện quy luật khí động và xác lập tiêu chí an toàn: Khẳng định định lượng mức độ ảnh hưởng của gió cạnh làm suy giảm lực nâng và tạo mô men mất ổn định; đưa ra khái niệm "góc trượt cạnh tới hạn" ($\beta^$) và chứng minh quy luật: khi vận tốc hạ cánh giảm thì góc trượt cạnh tới hạn $\beta^$ giảm theo, xác định ngưỡng nguy hiểm đối với máy bay hạ cánh.
- Bổ sung bộ dữ liệu đặc tính khí động: Cung cấp hệ thống dữ liệu ĐTKĐ ($C_y, C_x, C_z, m_x, m_y, m_z, x_F, z_F$) của máy bay L-39 ở cấu hình cất/hạ cánh dưới tác động của gió cạnh bằng phương pháp số, bổ sung vào khoảng trống dữ liệu mà tài liệu bay thực nghiệm chưa có.
- Đóng góp về mặt phương pháp luận: Đề xuất và kiểm chứng quy trình tiếp cận khoa học sử dụng đồng thời hai phương pháp số độc lập (phương pháp XRR định hướng kết hợp phần mềm CFD Ansys phân tích chi tiết) và đối soát với tài liệu kỹ thuật, tạo ra giải pháp thay thế hiệu quả cho phương pháp thử nghiệm ống khí động khi không đủ điều kiện môi trường lạnh.
Kiến nghị và giải pháp đề xuất
- Ứng dụng mô hình và thuật toán đã xây dựng để tích hợp vào phần mềm của buồng tập lái mô phỏng (TL-39), giúp phi công rèn luyện kỹ năng xử lý gió cạnh khi tiếp đất.
- Đưa miền vận tốc an toàn và bảng tra cứu góc trượt cạnh tới hạn $\beta^*$ vào tài liệu hướng dẫn khai thác sử dụng và giáo trình huấn luyện bay trong điều kiện khí tượng phức tạp.
- Sử dụng công cụ mô phỏng số làm phương tiện trợ giúp trong thiết kế tối ưu hóa hình học các chi tiết cơ giới hóa cánh và thiết bị hạ cánh.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của đề tài
- Phạm vi khí động học hẹp: Nghiên cứu chỉ giới hạn trong dòng chảy bao dừng (steady flow), chưa khảo sát các hiện tượng dòng không dừng phát sinh do gió giật (gust) hoặc các xung gió cạnh thay đổi đột ngột theo thời gian.
- Giả thiết điều kiện biên: Gió cạnh được đơn giản hóa với hướng thổi ngang vuông góc cố định ($\delta = 90^\circ$) và cường độ không đổi trên toàn bộ trường dòng.
- Giới hạn mô hình XRR: Phương pháp XRR chưa mô tả được đầy đủ độ dày hình học thực của thân máy bay và hiệu ứng bóc tách lớp biên do ma sát nhớt ở góc tấn lớn, phải phụ thuộc vào Ansys để bổ chính kết quả.
- Đối tượng khảo sát: Toàn bộ các kết quả tính toán định lượng và kiểm chứng số liệu mới chỉ thực hiện trên đối tượng máy bay huấn luyện L-39.
Hướng nghiên cứu tiếp theo
- Mở rộng bài toán sang dòng chảy bao không dừng để mô phỏng tác động của các luồng gió giật và trường gió cạnh có vận tốc thay đổi theo thời gian và độ cao.
- Nghiên cứu tương tác khí động học - đàn hồi (aeroelasticity) của kết cấu cánh máy bay khi chịu tải trọng gió cạnh bất đối xứng ở vận tốc hạ cánh.
- Ứng dụng phương pháp luận cho các chủng loại khí cụ bay khác (máy bay vận tải, máy bay chiến đấu phản lực có hình dạng cánh phức tạp).
Giá trị tham khảo
Đối tượng thụ hưởng
- Nghiên cứu sinh và học viên sau đại học: Thuộc các chuyên ngành Kỹ thuật cơ khí động lực, Kỹ thuật hàng không, Khí động lực học và Động lực học bay tham khảo phương pháp luận kết hợp các công cụ mô phỏng số.
- Giảng viên và cán bộ huấn luyện bay: Tại các học viện quân sự (Học viện Phòng không - Không quân, Học viện Kỹ thuật Quân sự) và các cơ sở đào tạo phi công sử dụng để biên soạn giáo trình, bài bay thực hành hạ cánh có gió cạnh.
- Kỹ sư và cán bộ kỹ thuật: Tại các viện nghiên cứu chuyên ngành (Viện Tên lửa, Viện Kỹ thuật Quân sự PK-KQ, Trung tâm DASI - ĐH Bách khoa Hà Nội) sử dụng các chương trình tính toán và bộ dữ liệu để phục vụ công tác thiết kế, cải tiến và đánh giá an toàn phương tiện bay.
Các nội dung có giá trị tham khảo chuyên sâu
- Thuật toán ghép nối không đồng phẳng cho cánh tà lệch góc và càng hạ cánh trong phương pháp xoáy rời rạc (Chương 3).
- Phương pháp tính toán yêu cầu nhiệt độ đồng dạng dòng chảy trong ống khí động (Chương 1).
- Bảng tra cứu tương quan góc lái điều khiển $\delta_H, \delta_L$ và góc trượt cạnh $\beta^*$ phục vụ cân bằng khí động máy bay L-39 (Chương 4).
Câu hỏi thường gặp
1. Vì sao luận án không sử dụng phương pháp thử nghiệm trong ống khí động để kiểm chứng kết quả?
Việc thử nghiệm trong ống khí động (OKĐ) ở vùng vận tốc nhỏ đòi hỏi phải đảm bảo đồng thời sự đồng dạng về số Mach ($M$) và số Reynolds ($Re$). Để đáp ứng điều này với mô hình thu nhỏ ($k = 2,5$), nhiệt độ dòng khí trong OKĐ bắt buộc phải được làm lạnh xuống $-116,5^\circ\text{C}$. Các phòng thí nghiệm OKĐ hiện có tại Việt Nam không có hệ thống làm lạnh sâu môi chất khí, nếu thử nghiệm ở nhiệt độ thường sẽ phá vỡ điều kiện đồng dạng nén/nhớt, dẫn đến sai lệch lớn về hệ số lực nâng cực đại so với điều kiện bay thực.
2. Hai phương pháp số được phối hợp thực hiện như thế nào trong luận án?
Luận án sử dụng phương pháp Xoáy rời rạc (XRR) phi tuyến như một công cụ tính toán nhanh, giải tích số rõ ràng để khảo sát định hướng và xác định vùng hội tụ nghiệm. Sau đó, phần mềm CFD Ansys (có bản quyền) được dùng làm công cụ mô phỏng chính xác cao để giải hệ phương trình Navier-Stokes trên lưới không gian 3 chiều phức tạp. Kết quả của hai phương pháp được đối chiếu chéo; khi độ chênh lệch $< 15%$, kết quả được xem là đủ độ tin cậy và được so sánh kiểm chứng tiếp với Tài liệu kỹ thuật (TLKT) của máy bay L-39.
3. Gió cạnh tác động như thế nào đến các hệ số khí động của máy bay khi hạ cánh?
Khi xuất hiện góc trượt cạnh $\beta$ do gió cạnh thổi ngang:
- Hệ số lực nâng $C_y$ bị suy giảm dần.
- Hệ số lực cản $C_x$ và hệ số lực cạnh $C_z$ tăng lên.
- Phát sinh mô men nghiêng cánh $m_x$ và mô men hướng $m_y$ làm quay lệch thân máy bay khỏi trục đường băng.
- Vị trí các tiêu điểm khí động $x_F$ và $z_F$ bị dịch chuyển, làm thay đổi đặc tính ổn định tĩnh của máy bay.
4. Khái niệm "góc trượt cạnh tới hạn" ($\beta^*$) được xác định như thế nào?
Góc trượt cạnh tới hạn $\beta^$ là giá trị góc trượt cạnh lớn nhất mà tại đó hệ số lực nâng chưa bị suy giảm nghiêm trọng và hệ thống cánh lái (cánh lái hướng $\delta_H$, cánh lái liệng $\delta_L$) còn đủ hiệu lực để tạo mô men cân bằng lại các mô men do gió cạnh gây ra. Luận án chỉ ra rằng $\beta^$ tỷ lệ thuận với vận tốc hạ cánh: khi vận tốc hạ cánh giảm từ $200\text{ km/h}$ xuống $160\text{ km/h}$, góc trượt cạnh tới hạn $\beta^*$ giảm tương ứng từ $17,2^\circ$ xuống $10,5^\circ$.
5. Kết quả nghiên cứu có thể ứng dụng trực tiếp vào thiết bị nào trong thực tế?
Kết quả nghiên cứu và thuật toán mô hình hóa được ứng dụng trực tiếp để nâng cấp phần mềm khí động học trong buồng tập lái mô phỏng TL-39, giúp tái lập chính xác phản ứng của máy bay khi gặp gió cạnh ở các giai đoạn bay xuống và cải bằng, từ đó phục vụ đào tạo phi công xử lý tình huống nguy hiểm khi hạ cánh.
Kết luận
Luận án tiến sĩ kỹ thuật của nghiên cứu sinh Nguyễn Đình Sơn đã giải quyết bài toán xác định ảnh hưởng của gió cạnh đến đặc tính khí động học của máy bay trong quá trình hạ cánh ở vùng vận tốc nhỏ. Bằng việc phát triển mô hình xoáy rời rạc phi tuyến kết hợp với phần mềm động lực học dòng chảy Ansys, công trình đã làm sáng tỏ quy luật biến thiên của các hệ số lực nâng, lực cản, lực cạnh và các mô men khí động theo góc trượt cạnh trên máy bay L-39. Luận án đã xác lập tiêu chí góc trượt cạnh tới hạn phụ thuộc vào vận tốc hạ cánh, đề xuất miền vận tốc an toàn và giải pháp điều khiển cân bằng, đóng góp thiết thực cho công tác đảm bảo an toàn bay và nâng cao chất lượng đào tạo, huấn luyện kỹ thuật hàng không quân sự.