Luận án tiến sĩ 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

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG VIỆN KHOA HỌC VÀ CÔNG NGHỆ QUÂN SỰ TRẦN MẠNH TUÂN 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 LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI - 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHỊNG VIỆN KHOA HỌC VÀ CƠNG NGHỆ QUÂN SỰ TRẦN MẠNH TUÂN 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 Ngành: Cơ kỹ thuật Mã số: 52 01 01 LUẬN ÁN TIẾN SĨ KỸ THUẬT NGƯỜI HƯỚNG DẪN KHOA HỌC PGS TS Nguyễn Văn Chúc TS Lê Đức Hạnh HÀ NỘI - 2022 i LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu Các số liệu, kết luận án trung thực chưa cơng bố cơng trình khác Các liệu tham khảo trích dẫn đầy đủ Hà nội, ngày 26 tháng năm 2022 Người cam đoan NCS Trần Mạnh Tn ii LỜI CẢM ƠN Cơng trình nghiên cứu thực Viện Tên lửa, Viện Khoa học Công nghệ quân - Bộ Quốc Phịng Lời cảm ơn chân thành tơi xin gửi đến hai thầy hướng dẫn PGS TS Nguyễn Văn Chúc TS Lê Đức Hạnh trực tiếp hướng dẫn, tận tình bảo, tạo điều kiện tốt để tơi hồn thành luận án Xin chân thành cảm ơn Ban giám đốc Viện Khoa học Công nghệ Quân sự, Thủ trưởng Phòng đào tạo, Thủ trưởng Viện Tên lửa tạo điều kiện mặt cho tơi q trình học tập nghiên cứu thực luận án Tôi xin chân thành cảm ơn nhà khoa học, cán nghiên cứu Viện Tên lửa, Khoa Hàng không vũ trụ Khoa Vũ khí - Học viện Kỹ thuật quân sự, Viện Cơ khí động lực - Đại học Bách khoa Hà Nội có ý kiến góp ý q giá q trình nghiên cứu hồn thiện luận án Tơi xin bày tỏ biết sâu sắc đến gia đình, người thân bạn bè, đồng nghiệp quan tâm, cổ vũ, động viên tạo điều kiện tốt cho thực luận án NCS Trần Mạnh Tuân iii MỤC LỤC Trang DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT vii DANH MỤC CÁC BẢNG xii DANH MỤC CÁC HÌNH VẼ xiii MỞ ĐẦU Chương TỔNG QUAN CÁC VẤN ĐỀ NGHIÊN CỨU 1.1 Tổng quan tối ưu hóa biên dạng khí động khí cụ bay 1.1.1 Bài toán tối ưu biên dạng khí động khí cụ bay 1.1.2 Phương pháp giải toán tối ưu biên dạng khí động 11 1.1.3 Tổng quan nghiên cứu tối ưu biên dạng khí động khí cụ bay 16 1.2 Tổng quan xu hướng phát triển khí cụ bay điều khiển kênh 18 1.2.1 Tổng quan khí cụ bay điều khiển kênh 18 1.2.2 Các nghiên cứu liên quan đến khí cụ bay điều khiển kênh 21 1.3 Phương pháp nghiên cứu không gian tham số 22 1.3.1 Vùng không gian tham số thiết kế 23 1.3.2 Hàm ràng buộc hàm mục tiêu 23 1.3.3 Giới hạn hàm mục tiêu tập điểm chấp nhận 24 1.3.4 Các điểm hiệu 25 1.3.5 Không gian mục tiêu 25 1.3.6 Sơ đồ thuật toán tối ưu 26 1.4 Những vấn đề tồn hướng nghiên cứu luận án 28 1.4.1 Những vấn đề tồn 28 1.4.2 Hướng nghiên cứu luận án 29 1.5 Kết luận chương 30 iv Chương XÂY DỰNG CHƯƠNG TRÌNH MƠ PHỎNG ĐỘNG LỰC HỌC BAY KHÍ CỤ BAY ĐIỀU KHIỂN MỘT KÊNH 31 2.1 Các hệ tọa độ sử dụng 31 2.2 Xây dựng hệ phương trình mơ tả chuyển động bay khí cụ bay 35 2.3 Lực mô-men tác dụng lên khí cụ bay 39 2.3.1 Trọng lực 40 2.3.2 Lực mơ-men khí động 42 2.3.3 Lực đẩy động 45 2.3.4 Lực căng dây vi cáp điều khiển 48 2.4 Phương pháp dẫn thuật toán điều khiển khí cụ bay 53 2.5 Xây dựng quỹ đạo dẫn 55 2.6 Xây dựng chương trình mơ động lực học bay 58 2.7 Kết luận chương 61 Chương XÂY DỰNG BÀI TOÁN TỐI ƯU BIÊN DẠNG KHÍ ĐỘNG CỦA KHÍ CỤ BAY ĐIỀU KHIỂN MỘT KÊNH 62 3.1 Thiết lập tốn tối ưu biên dạng khí động khí cụ bay điều khiển kênh sử dụng dây vi cáp 62 3.1.1 Xác định tham số thiết kế cần tối ưu 62 3.1.2 Các hàm mục tiêu ràng buộc toán tối ưu 63 3.1.3 Lựa chọn phương pháp giải 68 3.2 Phương pháp giải toán tối ưu 70 3.2.1 Thiết lập bảng kết tính tốn 71 3.2.2 Xác định tập điểm chấp nhận 77 3.2.3 Lựa chọn phương án thiết kế tối ưu 82 3.3 Kết luận chương 85 Chương TỐI ƯU THAM SỐ THIẾT KẾ BỘ CÁNH TRƯỚC CHO v KHÍ CỤ BAY ĐIỀU KHIỂN MỘT KÊNH CẢI TIẾN 87 4.1 Kiểm chứng phương pháp xác định tham số khí động 87 4.1.1 Phương pháp mô số 87 4.1.2 Thực nghiệm ống thổi khí động 90 4.2 Kiểm chứng chương trình mơ động lực học bay 93 4.2.1 Kiểm chứng khí cụ bay CT14M 93 4.2.2 Kiểm chứng phương án khí cụ bay CT14M cải tiến 98 4.3 Tối ưu tham số thiết kế cánh trước cho khí cụ bay điều khiển kênh cải tiến kiểu CT14M 100 4.3.1 Xác định ràng buộc tham số thiết kế 101 4.3.2 Xác định tập phương án chấp nhận 110 4.3.3 Lựa chọn phương án thiết kế tối ưu 114 4.4 Mô quỹ đạo bay phương án thiết kế tối ưu 120 4.4.1 Mục tiêu cự ly 500 m 120 4.4.2 Mục tiêu cự ly 2.500 m 121 4.5 Kết luận chương 122 KẾT LUẬN 123 DANH MỤC CÔNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ 125 TÀI LIỆU THAM KHẢO 126 PHỤ LỤC PL1 Phụ lục 1: Biên dạng khí động khí cụ bay CT14M .PL1 Phụ lục 2: Biên dạng khí động khí cụ bay CT14M cải tiến kiểm chứng phương pháp xác định hệ số khí động PL2 Phụ lục 3: Biên dạng khí động khí cụ bay CT14M cải tiến kiểm chứng chương trình mơ động lực học bay PL3 vi Phụ lục 4: Biên dạng khí động khí cụ bay CT14M cải tiến tối ưu .PL4 Phụ lục 5: Bảng hệ số rj  PL5 l Phụ lục 6: Chương trình tính tốn hệ số khí động khí cụ bay CT14M PL7 Phụ lục 7: Chương trình tính tốn hệ số khí động khí cụ bay CT14M cải tiến PL8 Phụ lục 8: Lực đẩy động khí cụ bay CT14M theo nhiệt độ PL9 vii DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT CD Gia tốc khí cụ bay [m/s2] Dây cung mút cánh [mm] Dây cung gốc cánh [mm] Giới hạn hàm ràng buộc Giới hạn hàm ràng buộc Ma trận cô-sin phương chuyển HTĐ liên kết sang HTĐ vận tốc - Ma trận cô-sin phương chuyển HTĐ cầu mặt đất sang HTĐ mặt đất - Hệ số lực cản Cdk - Đạo hàm hệ số lực điều khiển theo góc lật cánh lái Ceb Cx - Ma trận cô-sin phương chuyển HTĐ liên kết sang HTĐ mặt đất - Hệ số lực nâng - Hệ số lực cản khí động dọc trục Cy - Hệ số lực khí động pháp tuyến Cz - Hệ số lực khí động ngang - Ma trận cô-sin phương chuyển HTĐ vận tốc sang HTĐ mặt đất - Tập điểm chấp nhận - Tập điểm có khả - Độ xác dẫn [m] - Tập điểm hiệu - Tập điểm Pareto - Hàm ràng buộc - Tần số quay quanh trục dọc [Hz] - Tần số dao động riêng kênh chúc ngóc [Hz] - Lực khí động HTĐ vận tốc [N] a bk b0 c* c** Cbw CCe CL Cwe D  D Dmin E E f (A) fx fz Fxa , Fya , Fza  Fb Fd - - Ngoại lực HTĐ liên kết [N] - Lực căng dây vi cáp [N] viii F dk - Lực điều khiển [N] Fcbdk - Lực điều khiển cân [N] dk Fmax - Lực điều khiển cân cực đại [N] Fxedk , Fyedk , Fzedk Fyedk_ tb Fms  FM FP FxbP , FybP , FzbP F PHT FxbPHT , FybPHT , FzbPHT F PP FxbPP , FybPP , FzbPP - Lực điều khiển hệ tọa độ mặt đất [N] - Lực điều khiển trung bình kênh tầm [N] - Lực ma sát tời dây [N] - Lực Mechersky [N] - Lực đẩy động [N] - Lực đẩy động HTĐ liên kết [N] - Lực đẩy ĐCHT [N] - Lực đẩy ĐCHT HTĐ liên kết [N] - Lực đẩy ĐCP [N] - Lực đẩy ĐCP HTĐ liên kết [N] Fxb , Fyb , Fzb - Ngoại lực so với trục HTĐ liên kết [N] Fxba , Fyba , Fzba - Các thành phần lực khí động HTĐ liên kết [N] Fxbd , Fybd , Fzbd - Lực căng dây HTĐ liên kết [N] Fxbgra , Fybgra , Fzbgra - Trọng lực HTĐ liên kết [N] Fxbp , Fybp , Fzbp - Lực đẩy động HTĐ liên kết [N] Fy - Đạo hàm lực pháp tuyến theo góc α [N/độ] Fy - Đạo hàm lực pháp tuyến theo góc δ [N/độ] g - Gia tốc trọng trường [ m/s2] G - Tập điểm thỏa mãn ràng buộc Jx, J y, Jz K Kc Kd Kdk K od - Mơ-men qn tính HTĐ liên kết [ kg.m2] - Hệ số chất lượng khí động Động lượng hệ Hệ số hình dạng ống dây Hệ số tính điều khiển Độ dự trữ ổn định tĩnh 126 TÀI LIỆU THAM KHẢO Tiếng Việt [1] Nguyễn Văn Chúc, Phan Văn Chương, Trần Mạnh Tuân, Lê Quang Thương, Phan Thế Sơn, Đặng Văn Thức, Về giải pháp nâng cao độ xác xác định mơ-men qn tính tên lửa điều khiển, Tạp chí Nghiên cứu KH&CN quân sự, Số 62, 8/2019 [2] Nguyễn Văn Chúc, Phan Văn Chương, Trần Mạnh Tuân, Lê Quang Thương, Trần Phú Hoành, Lê Đức Hạnh, Nghiên cứu phối trí khí động phục vụ cải tiến, thiết kế tên lửa điều khiển kênh, Tạp chí Nghiên cứu KH&CN quân sự, Số 61, 6/2019 [3] Nguyễn Văn Chúc, Nguyễn Văn Sơn, Trần Phú Hoành, Phương pháp xác định hệ số lệnh tên lửa kênh quay quanh trục dọc, Tạp chí Nghiên cứu KH&CN quân sự, Số 45, 10/2016 [4] Nguyễn Văn Chúc, Nguyễn Phú Thắng, Phạm Khắc Lâm, Mô bán tự nhiên thời gian thực tên lửa điều khiển kiểu B72, Tạp chí Nghiên cứu KH&CN quân sự, số Đặc san Tên lửa, 09/2016 [5] Phan Văn Chương, "Nghiên cứu ảnh hưởng số tham số kết cấu đến đặc trưng động lực học hệ thống tên lửa kéo chuỗi nổ mềm liên tục", Luận án Tiến sĩ kỹ thuật, Viện KH-CNQS, 2011 [6] Tô Văn Dực, Nguyễn Văn Sơn, Động lực học bay nguyên lý dẫn khí cụ bay điều khiển kênh, NXB Khoa học kỹ thuật, 2006 [7] Vũ Tùng Lâm, Nghiên cứu phương pháp thiết kế tối ưu kết cấu vỏ khoang KCB đối hải âm, Luận án Tiến sĩ kỹ thuật, Viện KH-CNQS, 2016 [8] Tăng Xuân Long, Trịnh Hồng Anh, Hồ Xuân Vĩnh, Tối ưu thơng số thiết kế máy lái điện – khí tên lửa đối hải phương pháp Monte-Carlo, Tạp chí Nghiên cứu KH&CN quân sự, Số Đặc san KCB, 09/2016 127 [9] Nguyễn Trang Minh, So sánh số phương pháp tìm nghiệm tối ưu xây dựng sở mơ q trình tự nhiên, Tạp chí Nghiên cứu KH&CN quân sự, Số 33, 10/2014 [10] Bùi Phúc Kiển, Quy hoạch đa mục tiêu, Luận văn thạc sĩ Toán học, ĐHSP Thành phố Hồ Chí Minh, năm 2012 [11] Lê Văn Hiệp, Một lớp phương pháp giải toán tối ưu nhiều mục tiêu, Luận văn thạc sĩ tốn học, Đại học Quốc gia Thành phố Hồ Chí Minh, năm 2009 [12] Bùi Thế Tâm, Trần Vũ Thiệu, Các phương pháp tối ưu hóa, Nxb Giao thơng vận tải, Hà Nội, năm 1998 [13] Nguyễn Hải Thanh, Tối ưu hóa, Nxb Đại học Bách khoa, Hà Nội, năm 2006 [14] Nguyễn Đình Thúc, Trí tuệ nhân tạo – Lập trình tiến hóa, Nxb Giáo dục, Hà Nội, 2002 [15] Tổ hợp KCB chống tăng 9K11, Cục kỹ thuật/ Bộ tư lệnh Pháo binh, Hà Nội, 2006 [16] Trần Mạnh Tuân, Nguyễn Phú Thắng, Nguyễn Văn Chúc, Đỗ Tiến Cần, Phạm Khắc Lâm, Xác định hệ số khí động tên lửa B72 quay quanh trục dọc, Tạp chí Nghiên cứu KH&CN quân sự, Số 45, 10/2016 [17] Hồng Quang Tuyến, Bài giảng Tốn tối ưu, Trường ĐH Sự phạm Đà Nẵng, TP Đà Nẵng, 2012 Tiếng Anh [18] Ahmed Z Al-Garni, Ayman H Kassem, Ayman M Abdallah, Aerodynamic-Shape Optimization of Supersonic-Missiles Using Monte-Carlo, International Review of Aerospace Engineering (I.RE.AS.E), 02/2008 [19] ANSYS CFX Release 12.0: “ANSYS CFX-Solver Theory Guide” [20] Anderson M B., Missile Aerodynamic Shape Optimization Using 128 Genetic Algorithms, Journal Of Spacecraft And Rockets, Vol 37, No 5, September–October2000 [21] Colonno M., Palacios F., Economon T D., Lonkar A K., Alonso J J., An Adjoint-Based Aerodynamic Shape Optimization Methodology for Fairing Systems, In 31st AIAA Applied Aerodynamics Conference (p 2649), 2013 [22] Deb Kalyanmoy, Multi-objective optimization using evolutionary algorithms, Chichester, UK: Wiley; 2001 [23] Feyzioğlu E., Roll Characterıstıcs And Shape Optımızatıon Of The Free-To-Rotate Taıl-Fıns On A Canard-Controlled Mıssıle, Doctoral Dissertation, Mıddle East Technıcal Unıversıty, 2014 [24] Holland J H., Adaptation in Natural and Artificial Systems, Cambridge: MIT Press, 1992 [25] Kennedy J., Eberhart R.C , Particle swarm optimization, In Proceedings of the IEEE International Joint [26] Ocokoljic G.J , Aerodynamic shape optimization of guided missile based on wind tunel testing and computational fluid dynamic simulation, “Thermal Science”, vol 21, No3, 2017 [27] Omer Tanrikulu, Veysi Ercan, Optimal external configuration design of unguided missiles, AIAA-97-3725, 1997 [28] Riddle D B., Hartfield R J., Burkhalter J E., & Jenkins R M., Geneticalgorithm optimization of liquid-propellant missile systems, Journal of Spacecraft and Rockets, 46, 151-159, 2009 [29] Runduo C., Xiaobing Z., Multi-objective optimization of the aerodynamic shape of a long-range 2018 [30] Simulink Destop Realtime User’s guide, MathWorks MATLAB, R2015 [31] Sivanandam S N., Algorithms,Springer, 2008 Deepa S N., Introduction to Genetic 129 [32] Srinivas N., Deb K., Multi-Objective function optimization using the non-dominated sorting genetic algorithm, Evolutionary Computaion, vol 2, N°3, pp.221-248 [33] Tanil C., Platin B E., External configuration optimization of missiles in conceptual design, 2009 [34] Usta E., Arslan K., Tuncer I H., Aerodynamic design analysis of missile with strake configuration at supersonic mach numbers, 8th Ankara International Aerospace Conference 10-12 September., 2015 [35] Vidanovic, N., Rasuo, B., Kastratovic, G., Maksimovic, S., Curcic, D., Samardzic M., Aerodynamic–structural missile fin optimization, Aerospace Science and Technology, 65, 26-45, 2017 [36] Wiliam B Blake, Missile Datcom User manual, 1998 [37] Xia C C., Jiang T T., Chen W F., Particle Swarm Optimization of Aerodynamic Shapes With Nonuniform Shape Parameter–Based Radial Basis Function, Journal of Aerospace Engineering, 30(3), 04016089, 2016 [38] Xiaobing Z Runduo C., Multi-objective optimization of the aerodynamic shape of a long-range, 2018 [39] Yang Y R., Jung S K., Cho T H., & Myong R S., Aerodynamic Shape Optimization System of a Canard-Controlled Missile Using TrajectoryDependent Aerodynamic Coefficients, Journal of Spacecraft and Rockets, 49(2), 243-249, 2012 [40] Yugoimport-Sdpr, Upgraded malyutka anti-tank missile family, Belgrade, 2006 Tiếng Nga [41] Артоболевский И И., Генкин М Д., Гринкевич В К., Соболь И М., Статников Р Б., Оптимизацияв теории машин ЛП-поиском, Докл АН СССР,1971, том 200, номер 6, 1287г –1290г [42] Ветров В.В., Грязев М.В., Основы устройства и функционирования 130 противотанковых управляемых рaкет, изд ТулГУ, Тула 2006г [43] Дьяконов В., Круглов В., МАТЛАВ Анализ, идентификация и моделирование систем Изд Пртер С Петербург 2002г [44] Каликин В.В., Численные методы, - М.: Наука, 1978г [45] Кашин В.М., Лифиц А.Л., Ефремов М.И Основы проектирования переносных зенитных ракетных комплексов Москва, Изд-во МГТУ им Н Э Баумана – 2014г [46] Колесников К.С., Динамика ракетa, Машиностроение, 2003г [47] Коростелев О.П., Теоретические основы проектирования ствольных управляемых ракет Defense Express Library Киев 2007г [48] Лебедев А.А Чернобровкин Л.С Динамика полета беспилотных летательных аппаратов Москва Издательство Машиностроение 1973г [49] Лотов А.В., Поспелова И.И., Многокритериальные задачи принятия решений: Учебное пособие – М.: МАКС Пресс, 2008г – 197 с [50] Мануйленко В.Г., Удин Е.Г., Теоретическик основы построения крылатых управляемых, Санкт-Петербург, 2020г [51] Мещерский И.В., Работы по механике тел переменной массы, изд Технико-тоеретической литературы, Москва, 1952г [52] Наземная аппаратура управлерия снарядом 9М14М (9М14) –Пуль управления 4С415 Москва: Военное издательство Министерства обороны СССР, 1972г [53] Назаренко А.М , Эффективный алгоритм многокритериальной суррогатной оптимизации, Выпускная квалификационная работа на степень магистра студента 711гр., Работа выполнена в ИППИ РАН, Москва, 2013г [54] Соболь И М., Статников Р Б ЛП-поиск и задачи оптимального конструирования, Проблемы случайного поиска — Рига: Зинатне, 131 1972г -№ c 117-135 [55] Соболь И.М., Стаников Р.В., Выбор оптимальных параметров в задачах со многими критериями Москва: Дрофа, 2006г [56] Соляник П.Н., М.Л.Сургайло, В.В.Чмовж, Экспериментальная Аэродинамика, Харьков "ХАИ", 2007г [57] Техническое описание наземной аппаратуры управления 9М14 Москва Издательство Машиностроение, 1967г [58] Управляемый снаряд 9М14М Техническое описание Изд Министр Обороны СССР 1966г [59] Федоренко Р П Приближенные решения задач оптимального управления — М.: Наука, 1978г PL1 PHỤ LỤC Phụ lục 1: Biên dạng khí động khí cụ bay CT14M PL2 Phụ lục 2: Biên dạng khí động khí cụ bay CT14M cải tiến kiểm chứng phương pháp xác định hệ số khí động PL3 Phụ lục 3: Biên dạng khí động khí cụ bay CT14M cải tiến kiểm chứng chương trình mơ động lực học bay PL4 Phụ lục 4: Biên dạng khí động khí cụ bay CT14M cải tiến tối ưu PL5 l  Phụ lục 5: Bảng hệ số rj j=1 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 l=1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 1 3 3 7 3 5 7 3 1 3 15 11 13 11 13 13 11 15 15 17 13 15 31 27 15 29 21 23 19 11 25 13 17 25 29 51 61 43 51 59 47 57 35 53 19 51 61 37 33 11 39 63 27 17 15 23 85 67 49 125 25 109 43 89 69 113 47 55 97 37 83 103 27 13 33 115 41 79 255 79 147 141 89 173 43 25 115 97 19 97 197 101 255 29 203 65 195 177 105 17 257 465 439 177 321 181 225 235 103 411 233 59 353 329 463 385 111 475 451 263 19 249 275 10 771 721 1013 759 835 949 113 929 615 157 39 761 169 983 657 647 581 505 833 139 147 203 81 11 1285 823 727 267 833 471 1601 1341 913 1725 2021 1905 375 893 1599 415 605 819 975 915 1715 1223 1367 12 3855 4091 987 1839 4033 2515 579 3863 977 3463 2909 3379 1349 3739 347 387 2381 2821 1873 1959 1929 2389 3251 13 4369 4125 5889 6929 3913 6211 1731 1347 6197 2817 5459 8119 5121 7669 2481 7101 2677 1405 7423 725 2465 471 2887 14 13107 4141 6915 16241 11643 2147 11977 4417 14651 9997 2615 13207 13313 2671 5201 11469 14855 12165 5837 5387 12483 12945 1279 15 21845 28723 16647 16565 18777 3169 7241 5087 2507 7451 13329 8965 19457 18391 3123 11699 721 709 20481 19285 13057 32321 4865 16 65535 45311 49925 17139 35225 35873 63609 12631 27109 12055 35887 9997 1033 31161 32253 15865 26903 41543 12291 5165 28931 29377 64771 17 65537 53505 116487 82207 102401 33841 81003 103445 5205 44877 97323 75591 62487 12111 78043 49173 100419 57545 86017 27985 54019 127427 24321 18 196611 250113 83243 50979 45059 99889 15595 152645 91369 24895 83101 226659 250917 259781 63447 147489 206167 77163 12303 69809 21251 103759 42247 19 327685 276231 116529 252717 36865 247315 144417 130127 302231 508255 320901 187499 234593 36159 508757 81991 241771 357231 299025 128325 62233 472541 338691 20 983055 326411 715667 851901 299009 1032727 685617 775365 172023 574033 810643 628265 308321 232401 974837 802875 987201 378135 774207 164575 248081 1008719 599831 PL6 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 1 3 3 3 3 1 3 1 7 3 5 7 3 13 11 15 15 11 13 11 13 13 11 31 11 23 27 19 21 17 13 15 31 27 15 29 21 23 19 11 25 13 29 21 13 31 25 49 33 19 29 11 19 27 15 25 21 51 19 29 39 37 43 17 29 119 75 73 105 59 65 21 113 61 89 45 107 25 49 111 119 33 67 101 19 39 91 47 169 109 149 15 213 135 253 191 155 175 63 47 23 39 55 71 87 103 119 135 153 169 185 201 217 369 393 167 333 473 469 101 21 451 229 247 297 403 497 61 361 215 393 125 285 501 277 199 301 19 83 351 337 829 989 375 365 131 215 733 451 447 177 57 471 979 197 251 517 137 567 1021 167 877 929 269 327 997 91 663 629 525 469 981 1667 1587 1251 451 481 721 483 1209 1457 415 1435 725 861 41 1619 1579 1701 869 1151 1897 1679 1355 1149 243 3609 1131 1701 143 1339 3497 2499 1571 983 4021 1625 3217 1163 2977 3391 675 3093 1495 3443 557 675 1489 2303 3925 3705 1715 5595 5689 441 3169 4485 6311 3557 483 3781 3195 5213 5085 185 7323 1713 4021 5875 5363 4977 5441 1779 6777 287 6919 1517 1875 187 8133 11819 14471 7615 2981 4081 7223 11843 10799 9277 2031 15371 6603 7563 11617 4129 12061 3471 15919 1097 10369 14343 8475 16139 305 7621 12285 4929 15889 12625 8405 12593 28637 13425 28285 15893 15405 4677 19493 1129 25321 14979 4099 25469 17589 6731 13483 15325 18465 6929 16677 21765 4381 53631 10817 48083 8881 41135 60913 60935 58577 12029 959 19637 26607 56445 36087 52563 5455 12345 47423 50131 43771 58779 33331 63615 46013 34579 45827 9079 110851 8261 67537 34707 106823 15703 94129 69521 86021 19793 87283 20931 26369 66817 37745 68289 102733 29505 33137 23313 36561 118321 43349 52785 120981 91157 94533 4357 189901 63993 85105 107847 26967 109273 217151 217093 213491 186143 54345 27399 98051 81777 209987 21287 124097 98739 151281 116819 59665 30799 75249 239693 113679 37261 153 255947 336469 479495 339031 507907 475921 424277 348165 377941 343297 259163 521499 451841 235347 346179 128115 444613 361365 270519 420599 498897 322567 14035 73299 204881 431301 0.00052 734787 749285 911133 977907 344073 281389 789985 176165 414943 1041185 741087 132383 175361 539895 521289 20689 430923 426737 11187 998391 494137 939017 507165 863545 761911 176455 PL7 Phụ lục 6: Chương trình tính tốn hệ số khí động khí cụ bay CT14M $FLTCON NALPHA=5.,NMACH=1.,MACH=0.35,REN=7600000., ALPHA=-10.,-5.,0.01,2.,10.,$ $REFQ XCG=0.498, SREF = 0.01225, LREF = 0.125, LATREF = 0.125,$ $AXIBOD X0 = 0., NX =6., X(1)=0.000,0.0115,0.175,0.690,0.764,0.778, R(1)=0.0001,0.0085,0.0625,0.0625,0.057,0.057, BNOSE = 0.012, DEXIT = 0.013, TRUNC = FALSE.,$ $FINSET1 SECTYP = HEX, NPANEL=4., PHIF=45.,135.,225.,315., XLE=0.530, SWEEP=45.,16., STA=0., CHORD=0.150,0.150,0.150, SSPAN=0.0625,0.110,0.1965, ZUPPER=0.0063,0.0063,0.0063, LMAXU=0.0313,0.0313,0.0313, LFLATU=0.9375,0.9375,0.9375,$ $DEFLCT DELTA1 = -3.25,-3.25,-3.25,-3.25,$ $PROTUB NPROT = 5., PTYPE = HCYL,BLOCK,HCYL,BLOCK,BLOCK,VCYL, XPROT = 0.5725,0.326,0.350,0.540,0.670,.335, NLOC = 1.,1.,4.,2.,2.,2., LPROT = 0.168,0.056,0.032,0.008,0.012,.015, WPROT = 0.022,0.018,0.019,0.012,0.008,.015, HPROT = 0.025,0.009,0.019,0.012,0.012,.0105, OPROT = 0.,0.,0.,0.,0.,0.,$ DIM M DERIV DEG SAVE DAMP NEXT CASE PL8 Phụ lục 7: Chương trình tính tốn hệ số khí động khí cụ bay CT14M cải tiến $FLTCON NALPHA=7.,NMACH=1.,MACH=0.3,REN=7600000., ALPHA=-10.,-5.,-2.,0.01,2.,5.,10.,$ $REFQ XCG=0.645, SREF = 0.01225,LREF = 0.125, LATREF = 0.125,$ $AXIBOD X0 = 0., NX = 8., X(1)=0.0001,.004,0.014,.230,.310,.880,0.954,.968, R(1)=0.019,.027,0.032,.032,.0625,.0625,.057,0.057, BNOSE = 0.019, DEXIT = 0.013, TRUNC = TRUE.,$ $FINSET1 SECTYP = HEX, XLE=0.239,NPANEL=4., PHIF=0.,90.,180.,270., SWEEP=11., STA=0., CHORD=0.0418,0.0245, SSPAN=0.0443,0.0775, ZUPPER=0.0239,0.0408, LMAXU=0.1195,0.2041, LFLATU=0.761,0.5918,$ $FINSET2 SECTYP = HEX, NPANEL=4., PHIF=45.,135.,225.,315., XLE=0.720, SWEEP=45.,14., STA=0., CHORD=0.15,0.15,0.15, SSPAN=0.0625,0.110,0.1965, ZUPPER=0.0063,0.0063,0.0063, LMAXU=0.0313,0.0313,0.0313, LFLATU=0.9375,0.9375,0.9375,$ $DEFLCT DELTA2 = -3.25,-3.25,-3.25,-3.25,$ $PROTUB NPROT = 5., PTYPE = HCYL,BLOCK,HCYL,BLOCK,BLOCK,VCYL,VCYL, XPROT = 0.7455,0.499,0.523,0.713,0.843,.508,.2215, NLOC = 1.,1.,4.,2.,2.,2.,2., LPROT = 0.168,0.056,0.032,0.008,0.012,.015,0.002, WPROT = 0.022,0.018,0.019,0.012,0.008,.015,0.002, HPROT = 0.025,0.009,0.019,0.012,0.012,.0105,0.0087, OPROT = 0.,0.,0.,0.,0.,0.,0.,$ DIM M DERIV DEG DAMP NEXT CASE PL9 Phụ lục 8: Lực đẩy động khí cụ bay CT14M theo nhiệt độ STT Tên gọi Nhiệt độ Giá trị Đơn vị C -40 +15 +50 Lực đẩy động phóng kg 168 202 227 Thời gian cháy động phóng giây 0,81 0,68 0,61 Lực đẩy động hành trình kg 7,1 8,1 9,0 Thời gian cháy động hành trình giây 30,3 27,1 25,1