УДК 533.6.011
DOI 10.18698/2308-6033-2017-4-1610
Аэродинамические характеристики тонких конических оболочек при сверхзвуковых скоростях набегающего
потока
© А.Ю. Луценко, Д.К. Назарова, М.А. Фомин МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
Рассмотрены аэродинамические характеристики тонких отделяемых от ракет-носителей оболочек конической формы. Проведено математическое моделирование обтекания при сверхзвуковых скоростях набегающего потока, получены аэродинамические характеристики для различных конических оболочек. Построены графики зависимостей аэродинамических коэффициентов от угла атаки, отображающие влияние геометрических параметров исследуемых моделей на их аэродинамические характеристики. Показано, что результаты расчетов хорошо согласуются с имеющимися экспериментальными данными.
Ключевые слова: аэродинамические характеристики, ракета-носитель, отделяемые элементы, головной обтекатель, конические оболочки, моделирование обтекания, SolidWorks Flow Simulation
Введение. При достижении высоких слоев атмосферы осуществляется сброс головного обтекателя ракеты-носителя путем разделения его на несколько створок, которые достигают поверхности Земли. Траектория падения отделившегося элемента и размеры районов падения зависят от его аэродинамических характеристик (АДХ). Головные обтекатели большинства ракет-носителей, как правило, представляют собой тонкие оболочки и имеют цилиндрическую и коническую части [1].
Исследованию острых и затупленных конических тел как элементов конструкций сверхзвуковых летательных аппаратов уделялось большое внимание. Имеются многочисленные результаты как экспериментальных исследований, проведенных в ЦАГИ, ЦНИИмаше [2, 3] и др., так и численных расчетов [4, 5]. В то же время АДХ оболочек исследованы мало, получены только некоторые АДХ головных обтекателей конкретных ракет-носителей без проведения параметрических исследований [6-9].
В данной работе проанализированы АДХ конических оболочек, полученные путем численного моделирования в пакете SolidWorks Flow Simulation.
Постановка задачи. Рассмотрены тонкие конические оболочки с различными геометрическими характеристиками, представленными на рис. 1 и в таблице.
Аэродинамические коэффициенты оболочек определены в системе координат, представленной на рис. 1. При вычислении коэффициента Сх продольной силы, действующей на оболочку, использовалась
площадь = п ()2 (где Ь — длина модели; рк — угол полураствора конуса), а при расчете коэффициента Су нормальной силы —
площадь проекции модели на плоскость 0x2. Углы атаки а отсчи-тывались против хода часовой стрелки в плоскости 0ху.
Геометрические характеристики исследуемых моделей
Модель Угол полураствора конуса Рк, град Угол раскрытия оболочки у, град Формула расчета площади
1 12,5
2 18,0 360 ^ = ьЬ 1мРк
3 24,0
4 30,0
5 12,5
6 18,0 180 ^ = ьЬ 1мРк
7 24,0
8 30,0
9 18,0 90 42 2 2
В качестве тестового был проведен расчет АДХ конической оболочки с углом раскрытия у = 180° и углом полураствора конуса рк = = 12,5о (модель 5) при числе Маха набегающего потока Мда = 4 и углах атаки а в диапазоне значений 0.. .360° и выполнено сравнение по-
лученных результатов с экспериментальными данными [2]. На рис. 2, а и б представлены зависимости соответственно коэффициентов продольной и нормальной сил от угла атаки для модели 5 конической оболочки.
Сх
О 60 120 180 240 300 а, град
a
0 60 120 180 240 300 а, град
б
Рис. 2. Аэродинамические характеристики модели 5 при Мм = 4, полученные в результате моделирования (— ) и эксперимента ( • )
Результаты тестового расчета отличаются от экспериментальных данных менее чем на 10 %, что является приемлемым для оценки аэродинамических сил, действующих на оболочку в полете. Наибольшее расхождение данных наблюдается при поперечном обтекании оболочки как с вогнутой (а = 90°), так и с выпуклой (а = 270°) стороны, что может быть обусловлено наличием развитой отрывной зоны или влиянием поддерживающих устройств в эксперименте.
Анализ результатов расчета. Проведены параметрические расчеты АДХ различных моделей оболочек, представленных в таблице, исследовано влияние угла рк полураствора конуса и угла у раскрытия оболочки.
На рис. 3, а и б представлены зависимости соответственно коэффициентов продольной и нормальной сил моделей круговых конических оболочек (модели 1, 3 и 4) с различными углами полураствора при Мю = 3.
90
а
а,град
90
б
а,град
Рис. 3. Аэродинамические характеристики тонких конических оболочек с у = 360° и различными углами рк при Мм = 3
Поскольку все исследуемые модели симметричны относительно плоскости 0x2, расчеты проведены для углов атаки а = 0.. .180°.
Коэффициент Сх продольной силы оболочки возрастает по модулю при увеличении угла рк, так как увеличивается площадь наветренной поверхности. В донной части конического тела имеет место застойная область, давление в которой не изменяется при увеличении угла рк полураствора конуса, соответственно, увеличивается коэффициент волнового сопротивления, а вклад трения в общий коэффициент Сх тела при сверхзвуковых скоростях незначителен [8]. Коэффициент Сх конической оболочки принимает нулевые значения при а = 100.110°.
Коэффициент Су нормальной силы возрастает монотонно до значений а = 50...70° в зависимости от конфигурации исследуемого тела, после чего уменьшается. С увеличением угла рк максимальное значение Су возрастает по модулю и достигается при меньшем угле атаки. Это связано с тем, что при сверхзвуковых скоростях нормальная сила определяется в основном силами давления на наветренную сторону оболочки.
Сравнение АДХ оболочки и сплошного конуса со сходными геометрическими характеристиками [9] (см. рис. 3) показало, что значения коэффициентов нормальной и продольной сил оболочки и сплошного конуса при углах атаки до 70° практически совпадают. Это связано с тем, что при а < 90° АДХ определяются в основном формой внешней поверхности конуса (оболочки). Некоторые расхождения наблюдаются для зависимости Су(а) при а > 30°.
На рис. 4, а и б представлены зависимости соответственно коэффициентов продольной и нормальной сил от угла атаки для моделей разрезных оболочек с углом раскрытия у = 180° (модели 5-8) и различными углами полураствора конуса при Мда = 4.
Сх
-2
0 60 120 180 240 300 а, град
а
Су
-6
о 60 120 180 240 300 а, град
б
Рис. 4. Аэродинамические характеристики разрезных оболочек с у = 180° и различными углами рк при Мм = 4
Графики зависимостей Сх(а) и Су(а) для различных углов рк имеют схожий сложный характер и несимметричны относительно а = = 180°. С увеличением угла рк увеличивается и площадь поверхности оболочки, соответственно, повышенное давление с наветренной стороны создает большую силу, при этом с подветренной стороны образуется обширная отрывная область, давление в которой практически не зависит от рк.
Коэффициент Сх оболочек обращается в нуль при углах атаки, чуть больших рк, в интервалах 30.45° и 210.230°. Коэффициент Су также имеет два нулевых значения при углах атаки а1 ~ рк и а2 ~ 180° + + рк. Максимальные по модулю значения коэффициента нормальной силы оболочки получены при а3 ~ 90° + рк и а4 ~ 270° + рк, причем при обтекании со стороны вогнутой поверхности (углы атаки а = = 0.180°) коэффициенты Сх и Су по модулю больше, чем при обтекании со стороны выпуклой поверхности.
На рис. 5, а и б представлены зависимости соответственно Сх(а) и Су(а) для оболочек с углом полураствора конуса рк = 18° и углами раскрытия у = 360, 180 и 90° (модели 2, 6 и 9 соответственно).
180
а
а,град
180
б
а, град
Рис. 5. Аэродинамические характеристики конических оболочек с различными углами раскрытия у при Мм = 3
На рис. 5, a видно, что характер зависимостей коэффициента Сх разрезных конических оболочек от угла атаки а существенно отличается от характера зависимости для круговой оболочки с у = 360°. Значение Cx практически постоянно при а = 150...220°, это связано с особенностями обтекания внутренней полости со стороны дна тонкого конуса. Зависимость Су(а) круговой оболочки имеет схожий характер с зависимостями Су(а) разрезных оболочек, однако в случае круговой оболочки экстремум соответствует меньшему углу атаки а.
При углах атаки а > 210° коэффициенты нормальных сил оболочек практически совпадают, так как Су в этом случае определяется в большей степени формой наветренной поверхности.
При а < 210° значения Су и Сх оболочки с углом у = 180° по модулю больше значений соответствующих коэффициентов оболочки с углом у = 90°, что, как и при увеличении угла ß^ можно объяснить увеличением площади наветренной поверхности.
Наименьшие по модулю значения Cx имеет оболочка с углом раскрытия у = 90°, так как при уменьшении у исследуемая модель приближается по форме к треугольному крылу. Соответственно, зависимость Су(а) модели 9 также приближается к зависимости Су(а) для треугольного крыла [10, 11].
Выводы. В ходе комплексных исследований процесса обтекания тонких конических оболочек выявлены закономерности изменения АДХ в зависимости от углов атаки, полураствора конуса и раскрытия оболочки.
Для того чтобы подтвердить возможность применения пакета SolidWorks Flow Simulation к задачам определения АДХ тонких конических оболочек, был проведен тестовый расчет и выполнено сравнение полученных данных с имеющимися экспериментальными данными ЦАГИ.
Выявлено, что АДХ тонких оболочек зависят не только от числа Маха набегающего потока и угла атаки, но и в довольно большой степени от углов полураствора конуса и раскрытия оболочки. Аэродинамические характеристики круговой оболочки принципиально отличаются от соответствующих характеристик разрезных оболочек.
Полученные зависимости могут быть полезны при оценке АДХ головных обтекателей ракет-носителей для определения их траекторий и расчета зон падения, а также при решении задач стабилизации отделяемых частей.
ЛИТЕРАТУРА
[1] Потапов А.М., Коваленко В.А., Кондратьев А.В. Сравнение головных обтекателей существующих и перспективных отечественных ракет-носителей и их зарубежных аналогов. Авиационно-космическая техника и технология, 2015, № 1 (118), с. 35-43.
[2] Петров К.П. Аэродинамика тел простейших форм. Москва, Факториал, 1998, 432 с.
[3] Артонкин В.Г., Леутин П.Г., Петров К.П., Столяров Е.П. Аэродинамические характеристики острых и притупленных конусов при дозвуковых и сверхзвуковых скоростях. Тр. ЦАГИ. Москва, Изд. отдел ЦАГИ, 1972, № 1413, 93 с.
[4] Коваленко В.В., Кравцов А.Н., Мельничук Т.Ю. Сопротивление конических носовых частей при сверхзвуковом обтекании. Ученые записки ЦАГИ, 2011, № 1, с. 31-36.
[5] Савкина Н.В., Биматов В.И., Христенко Ю.Ф. Расчет обтекания и аэродинамических характеристик острого конуса на основе решения прямой задачи нелинейной аэробаллистики. Вестник Томского государственного университета, 2014, № 1 (27), с. 110-116.
[6] Харитонова А.Н., Шахов В.Г. Особенности решения пространственных задач обтекания ракет-носителей с надкалиберными головными обтекателями с использованием программного комплекса ANSYS FLUENT. Вестник Самарского государственного аэрокосмического университета, 2012, № 4 (35), с. 116-123.
[7] Xuechang Z., Xiaojing Y., Yan H. Aerodynamic Characteristics of Fairing Separation at Initial Opening Angle. Proceedings of the 1st International Conference on Mechanical Engineering and Material Science, 2012, pp. 259-262. DOI: 10.2991/mems.2012.160
[8] Дядькин А.А., Луценко А.Ю., Назарова Д.К. Математическое моделирование обтекания тонкостенных конструкций в до- и трансзвуковом диапазоне скоростей. Научный вестник Московского государственного технического университета гражданской авиации, 2016, № 223 (1), с. 45-50.
[9] Дядькин А.А., Крылов А.Н., Луценко А.Ю., Михайлова М.К., Назарова Д.К. Особенности аэродинамики тонкостенных конструкций. Космическая техника и технологии, 2016, № 3 (14), с. 15-25.
[10] Wood R. Supersonic Aerodynamics of Delta Wings. NASA Technical Paper. Virginia, Landley Research Center Hampton, 1988, no. 2771, 104 p.
[11] Scallings R., Lamb Jr., Lamb M. Wing-Alone Aerodynamic Characteristics for Hight Angles of Attack at Supersonic Speeds. NASA Technical Paper. Virginia, Landley Research Center Hampton, 1981, no. 1889, 188 p.
Статья поступила в редакцию 20.02.2017
Ссылку на эту статью просим оформлять следующим образом: Луценко А.Ю., Назарова Д.К., Фомин М.А. Аэродинамические характеристики тонких конических оболочек при сверхзвуковых скоростях набегающего потока. Инженерный журнал: наука и инновации, 2017, вып. 4. http://dx.doi.org/10.18698/2308-6033-2017-4-1610
Статья подготовлена по материалам доклада, представленного наXLIАкадемических чтениях по космонавтике, посвященных памяти академика С.П. Королёва и других выдающихся отечественных ученых - пионеров освоения космического пространства. Москва, МГТУ им. Н.Э. Баумана, 24-27января 2017 г.
Луценко Александр Юрьевич — канд. техн. наук, доцент кафедры «Динамика и управление полетом ракет и космических аппаратов» МГТУ им. Н.Э. Баумана. Автор более 70 научных публикаций в области аэродинамики струйных и отрывных течений, проектирования органов управления полетом. e-mail: aulutsenko@mail.ru
Назарова Динара Камилевна — аспирант кафедры «Динамика полета ракет и космических аппаратов» МГТУ им. Н.Э. Баумана. Область научных интересов — аэродинамика струйных и отрывных течений. e-mail: dknazarova@mail.ru
Фомин Максим Андреевич — студент 6-го курса кафедры «Динамика полета ракет и космических аппаратов» МГТУ им. Н.Э. Баумана. e-mail: fominm92@gmail.com
Aerodynamic characteristics of thin conical shells at supersonic speeds of the ram airflow
© A.Yu. Lutsenko, D.K. Nazarova, M.A. Fomin
Bauman Moscow State Technical University, Moscow, 105005, Russia
The article considers aerodynamic characteristics of thin conical shells separable from carrier rockets. Mathematical modeling of flow at the supersonic velocities of the ram airflow is performed; aerodynamic characteristics for various conical shells are obtained. The graphs of aerodynamic coefficients vs. angles of attack are plotted. They display the influence of the geometric parameters of the models under study on their aerodynamic characteristics. It is shown that the results of the calculations are in good agreement with the available experimental data.
Keywords: aerodynamic characteristics, launch vehicle, separable elements, nose fairing, conical shells, flow simulation, SolidWorks Flow Simulation
REFERENCES
[1] Potapov A.M., Kovalenko V.A., Kondratyev A.V. Aviatsionno-kosmicheskaya tekhnika i tekhnologiya — Aerospace engineering and technology, 2015, no. 1 (118), pp. 35-43.
[2] Petrov K.P. Aerodinamika tel prosteyshikh form [Aerodynamics of bodies of the simplest forms]. Moscow, Faktorial Publ., 1998, 432 p.
[3] Artonkin V.G., Leutin P.G., Petrov K.P., Stolyarov E.P. Aerodinamicheskie kharakteristiki ostrykh i prituplennykh konusov pri dozvukovykh i sverkhzvukovykh skorostyakh [Aerodynamic characteristics of sharp and blunted cones at subsonic and supersonic speed].Trudu TsAGI [Proceedings of TsAGI]. Moscow, TsAGI Publ., 1972, no. 1413, 93 p.
[4] Kovalenko V.V., Kravtsov A.N., Melnichuk T.Yu. Uchenye zapiski TsAGI — TsAGI Scientific Journal, 2011, no. 1, pp. 31-36.
[5] Savkina N.V., Bimatov V.I., Khristenko Yu.F. Vestnik Tomskogo gosudar-stvennogo universiteta - Tomsk State University Journal, 2014, no. 1 (27), pp. 110-116.
[6] Kharitonova A.N., Shakhov V.G. Vestnik Samarskogo gosudarstvennogo universiteta — Vestnik of the Samara State Aerospace University, 2012, no. 4 (35), pp. 116-123.
[7] Xuechang Z., Xiaojing Y., Yan H. Aerodynamic Characteristics of Fairing Separation at Initial Opening Angle. Proceedings of the 1st International Conference on Mechanical Engineering and Material Science, 2012, pp. 259-262. DOI: 10.2991/mems.2012.160
[8] Dyadkin A.A., Lutsenko A.Yu., Nazarova D.K. Nauchnyy vestnikMoskovskogo gosudarstvennogo tekhnicheskogo universiteta grazhdanskoy aviatsii — Bulletin of Moscow State Technical University of Civil Aviation, 2016, no. 223 (1), pp. 45-50.
[9] Dyadkin A.A., Krylov A.N., Lutsenko A.Yu., Mikhaylova M.K., Nazarova D.K. Aviatsionno-kosmicheskaya tekhnika i tekhnologiya — Aerospace engineering and technology, 2016, no. 3 (14), pp. 15-25.
[10] Wood R. Supersonic Aerodynamics of Delta Wings. NASA Technical Paper. Virginia, Landley Research Center Hampton, 1988, no. 2771, 104 p.
[11] Scallings R., Lamb Jr., Lamb M. Wing-Alone Aerodynamic Characteristics for Hight Angles of Attack at Supersonic Speeds. NASA Technical Paper. Virginia, Landley Research Center Hampton, 1981, no. 1889, 188 p.
Lutsenko A.Yu., Cand. Sc. (Eng.), Associate Professor, Department of Space Flight Dynamics and Control, Bauman Moscow State Technical University, author of over 70 scientific publications in the field of aerodynamics of jet and separation flows, the design of flight controls. e-mail: aulutsenko@mail.ru
Nazarova D.K., Postgraduate (Ph.D.) student, Department of Space Flight Dynamics and Control, Bauman Moscow State Technical University. Research interests: aerodynamics of jet and separation flows. e-mail: dknazarova@mail.ru
Fomin M.A., 6th year student, Department of Space Flight Dynamics and Control, Bauman Moscow State Technical University. e-mail: fominm92@gmail.com