Аэродинамические характеристики ракеты-носителя "Посейдон С-3"

 














Аэродинамические характеристики ракеты-носителя «Посейдон С-3»


Введение


Ракета-носитель (РН) «Посейдон С-3» [1] - американская баллистическая ракета, размещаемая на подводных лодках. Представляет собой двухступенчатую баллистическую ракету с тандемным расположением ступеней. Обе маршевых ступени оснащены твердотопливными ракетными двигателями. Предназначена для уничтожения площадных незащищённых целей типа городов, а также военных целей, в том числе стартующих баллистических ракет.

Первые испытания РН, проводимые на суше состоялись 16 августа 1968 года. После, 17 июля 1970 года состоялся первый пуск с борта ракетоносца. В течение первых пусков было принято решение сократить наземный цикл испытаний до 20 пусков. В ходе цикла последний пуск был проведён 29 июня 1970 года. Из 20 пусков 13 завершились успешно.

Принятие на вооружение ракет «Посейдон C-3» значительно повысило боевые возможности флота США. При неизменном количестве ракетоносцев размещенное на них количество боеголовок возросло в 2,6 раза. В процессе эксплуатации стартовая надёжность ракет «Посейдон C-3» составила 84%.

Целью исследования является изучение аэродинамики летательного аппарата, а именно величины коэффициента лобового сопротивления и производной коэффициента подъемной силы по углу атаки, которые в дальнейшем позволят рассчитать координату фокуса и сделать заключение об устойчивости данного ракетоносителя.

В данной работе используется метод поэлементного расчета, заключающийся в определении необходимых аэродинамических характеристик отдельно для каждого элемента с последующим определением этих характеристик для «целого» летательного аппарата.


1. Постановка задачи

аэродинамический летательный мах

Для ракеты-носителя «Посейдон С-3» необходимо определить следующие аэродинамические характеристики:

·коэффициент сопротивления трения при нулевом угле атаки Cx тр(M?, h);

·коэффициент сопротивления давления при нулевом угле атаки Cx давл(M?);

·коэффициент аэродинамической продольной силы Cx0(M?) для высоты 10 км и нулевого угла атаки;

·производную коэффициента аэродинамической нормальной силы по углу атаки C?y(M?) для высоты 10 км;

·производную коэффициента аэродинамической подъемной силы по углу атаки C?(M?) для высоты 10 км;

·коэффициент индуктивного сопротивления для различных углов атаки и высоты 10 км;

·коэффициент лобового сопротивления для различных углов атаки и высоты 10 км;

·координату фокуса летательного аппарата .

Значение коэффициентов определить для дискретных значений чисел Маха набегающего потока

M? = (0,1; 0,3; 0,5; 0,7; 0,9; 1; 1,1; 1,3; 1,5; 2; 2,5; 3; 3,5; 4; 4,5; 5),

Высот, км, h = (0,10,20,30,40,60)

и углов атаки, град, ? = (0,2,4,6,8).

Зависимости Cx тр(M?, h), Cx давл.(M?), Cxa(M?, ?), XF (M?) представить в табличном виде и на рисунках.


Рисунок 1 - Чертеж РН «Посейдон С-3»


2. Формирование расчетной схемы ракеты-носителя


В общем случае ракету-носитель можно рассматривать как совокупность корпуса летательного аппарата, одного или нескольких типов ускорителей в виде боковых блоков и одного или нескольких типов крыльев. В данном случае ракета-носитель состоит только из корпуса, представляющего собой сочетание конических и цилиндрических поверхностей.

Формирование расчетной схемы ракетоносителя заключается в уточнении и корректировке масштаба его изображения.

По рисунку 1 находятся две величины, выраженные в миллиметрах, - длина ракеты Lрис и наибольший диаметр корпуса Dрис. В таблице 1 приведены в метрах длина L и наибольший диаметр первой ступени Dдн. По этим данным определяются два масштаба изображения аппарата - масштаб по длине



и масштаб по диаметру



Приведение корпуса аппарата к телу вращения производится за счет исключения из состава конструкции элементов, слабо влияющих на аэродинамические силы, и коррекции формы корпуса введением эквивалентного диаметра для отдельных частей корпуса. В результате формируется расчетная схема летательного аппарата, приведенная на рисунке 2.


Таблица 1 - Основные геометрические параметры РН

НаименованиеОбозначениеРазмерностьВеличинаДлина ракетыLм10,385Носовая частьДлинаLнм2,3ДиаметрDнм1,78Цилиндрический участок 1ДлинаLц1м0,682ДиаметрDц1м1,78Расширяющаяся частьДлинаLюм0,597Диаметр меньшийDю1м1,78Диаметр большийDю2м1,88Цилиндрический участок 2ДлинаLц2м6,223ДиаметрDц2м1,88Сужающаяся частьДлинаLужм0,259Диаметр меньшийDуж1м0,908Диаметр большийDуж2м1,49СоплоДлинаLсм0,324Диаметр меньшийDс1м0,908Диаметр большийDс2м1,102Количествоncшт1

Рисунок 2 - Расчетная схема РН «Посейдон С-3»


3. Геометрические параметры летательного аппарата


Летательный аппарат, схема которого приведена на рисунке 2, имеет следующие геометрические параметры:


- удлинение

площадь миделя


площадь боковой поверхности корпуса



Приведем расчет характеристик РН по каждому элементу:

. Носовая часть


- удлинение носовой части


. Расширяющаяся часть

угол полураствора конуса

длина полного конуса

длина фиктивного конуса

- удлинение

площадь меньшего основания

площадь большего основания


4. Коэффициент силы сопротивления трения ракеты-носителя при нулевом угле атаки


При безотрывном обтекании поверхности корпуса летательного аппарата пограничным слоем, пренебрегая влиянием кривизны поверхности на величину силы трения, а также наклоном отдельных элементов поверхности к продольной оси корпуса, коэффициент трения определяют следующим образом:



где Сf M=0 - коэффициент трения одной стороны плоской пластины в несжимаемом потоке;

?м - коэффициент влияния сжимаемости воздуха на сопротивление трения;

Fкорп. - площадь боковой поверхности корпуса, безотрывно обтекаемая пограничным слоем, м2;

Sм корп. - площадь миделя корпуса, м2.

Удвоенный коэффициент трения 2Сf M=0 определяется по формуле зависимости от числа Рейнольдса Re и относительной координаты перехода ламинарного пограничного слоя в турбулентный

для ламинарного пограничного слоя при значениях чисел Рейнольдса в диапазоне Re < 485000


для турбулентного пограничного слоя при значениях чисел Рейнольдса в диапазоне Re > 107



- для смешанного пограничного слоя при значениях чисел Рейнольдса в диапазоне 485000 < Re < 107



Число Рейнольдса определяется по длине корпуса L



где a? - скорость звука на заданной высоте, м/c;

? - кинематический коэффициент вязкости, м2.

Скорость звука и коэффициент кинематической вязкости определяются по таблице стандартной атмосферы в зависимости от высоты полета летательного аппарата [1].

Относительная координата точки перехода определяется исходя из двух условий, указанных в фигурных скобках



где Lн - длина первой носовой части корпуса летательного аппарата, м;

hw - средняя высота бугорков шероховатости поверхности, м.

Высота бугорков шероховатости выбирается из таблицы 1 [1] для различных поверхностей. Выбор точки перехода осуществляется следующим образом. Если координата попадает в интервал (0,2088; 0,3588), то нужно проводить корректировку и считать, что будет в месте перехода из конической в цилиндрическую часть.

Коэффициент ?м определяется по рисунку 4 [1] для смешанного пограничного слоя. Для полностью ламинарного или турбулентного пограничного слоя коэффициент ?м определяется по формулам

для ламинарного режима течения:



- для турбулентного режима течения:



Расчет коэффициента сопротивления трения РН при полете на высоте

Н=10 км представлен в таблице 2, для других высот в приложении А.

Графики зависимости коэффициента сопротивления трения от числа Маха для различных высот приведены на рисунке 3.


Таблица 2 - Коэффициент сопротивления трения РН при H=10 км

М0,10,08830,99930,00520,05410,30,26480,99400,00520,05330,50,44130,98370,00480,04890,70,61790,96860,00460,04580,80,70610,95950,00450,04450,90,79440,94940,00440,043310,88270,93840,00430,04221,10,97090,92670,00430,04111,31,14750,90110,00420,03901,51,32400,87350,00410,037121,76530,79910,00390,03262,52,20670,72350,00380,028732,64800,65190,00370,02523,53,08930,58670,00370,022243,53060,52890,00360,01974,53,97200,47810,00350,017554,41330,43380,00350,0157

Немонотонность зависимости, изображенной на рисунке 3, при H=30 км, объясняется переходом из ламинарного в смешанный поток, а из него в турбулентный.


Рисунок 3 - Зависимость коэффициента трения РН от числа Маха


5. Коэффициент сопротивления давления ракеты-носителя при нулевом угле атаки


Состав коэффициента сопротивления давления летательного аппарата

Коэффициент продольной силы давления летательного аппарата определяется через его составляющие для отдельных частей аппарата. Однако в данном случае летательный аппарат представляет собой только корпус и коэффициент продольной силы давления



Коэффициент сопротивления давления носовых частей

Сопротивление носовых частей зависит от величины их удлинения и формы контура в продольной плоскости симметрии, проходящей через продольную ось. Значение коэффициента сопротивления давления носовой части Cxa давл.н возрастает по мере увеличения чисел Маха от нулевого значения M? = 0 до значений M? = 1,5, соответствующих присоединению скачка уплотнения. При дальнейшем увеличении числа Маха значения коэффициента монотонно убывают в связи с увеличением угла наклона присоединенного скачка и уменьшением потерь энергии в нем.

Сопротивление носовой части оживальной формы, с удлинением , не приведено в справочных данных учебного пособия. Плоским сечением носовой части оживальной формы является парабола, описываемая уравнением:



варьируется в зависимости от различных чисел Маха и различных удлинений

Исходя из этого, можно написать систему уравнений для 3 удлинений:

оживала с

оживала с

плоского торца

Система будет иметь вид:



Далее находятся коэффициенты и, зная удлинение для исследуемой ракеты, находится коэффициент сопротивления давления носовой части при определенном числе Маха.

Данные, полученные с помощью решения данной выше системы уравнений приведены в таблице 3, а график зависимости изображен на рисунке 4.

Таблица 3 - Коэффициент сопротивления давления носовой части

M0,10,1650,30,1650,50,1490,70,1750,90,25510,3041,10,361,30,4841,50,520,512,50,50230,53,50,48540,4634,50,45850,447

Рисунок 4 - Зависимость коэффициента давления носовой части от числа Маха


Коэффициент носового сопротивления давления переходников в виде

усеченного конуса определяется по формуле

,


где - коэффициент сопротивления давления достроенного конуса [2, рис. 6,7]

- площади миделя цилиндрических частей соответственно с меньшим и большим диаметром, .

Расчет коэффициента сопротивления давления расширяющихся переходных частей РН представлен в таблице 4. Зависимость коэффициента сопротивления давления переходника представлена на рисунке 5.


Таблица 4 - Коэффициент сопротивления давления переходника в виде усеченного конуса

M0,10,00010,30,00050,50,00310,70,00420,90,005210,00781,10,00831,30,00941,50,009820,00722,50,005730,00523,50,004740,00444,50,004250,0041

Рисунок 5 - Зависимость коэффициента давления переходника от числа Маха


Коэффициент сопротивления давления донной части

Донное сопротивление обусловлено возникновением разряжения за тупым основанием тела. Величина разряжения, устанавливающаяся за донным срезом корпуса, зависит от многих факторов: формы кормовой части, наличия или отсутствия хвостового оперения, реактивной струи, длины корпуса, состояния пограничного слоя, температуры поверхности и других факторов.

При числах Маха >0,8 коэффициент донного сопротивления можно определить по формуле


,


где - коэффициент донного давления для тел вращения без сужающейся кормовой части [2, рисунок13];

- коэффициент, учитывающий форму кормовой части [2, рисунок14];

- площадь донного среза.

В случае РН «Посейдон С-3» площадью донного среза будет площадь кольца, образованного миделевым диаметром РН и фермой крепления ЖРД.



Расчет сопротивления донной части при >0,8 представлен в таблице 5.


Таблица 5 - Сопротивление давления донного среза при

M0,90,1480,00610,1950,0081,10,2000,0091,30,1950,0081,50,1850,00820,1570,0072,50,1200,00530,0960,0043,50,0800,00340,0680,0034,50,0590,00350,0500,002

При числах Маха коэффициент донного сопротивления можно определить по формуле


где - коэффициент донного давления для тел вращения без сужающейся кормовой части, определяется по формуле



где - коэффициент сопротивления трения плоской пластины, длина которой равна длине корпуса с учетом влияния сжимаемости, определяется по формуле



Расчет коэффициента сопротивления донной части при полете с

на высоте Н=10 км представлен в таблице 6.


Таблица 6 - Сопротивление давления донного среза при

M0,10,00260,12920,00560,00260,99930,30,00260,13020,00560,00260,99400,50,00240,13600,00590,00240,98370,70,00230,14040,00600,00220,9686

Коэффициент сопротивления давления корпуса

Коэффициент сопротивления давления РН определяется по формуле



Значения коэффициента сопротивления давления РН в диапазоне высот от 0 до 60 км представлены в таблице 7.

Таблица 7 - Сопротивление давления РН

MH, км010203040600,10,15370,15360,15280,15410,15150,15020,30,15450,15400,15380,15300,15300,15130,50,14300,14250,14230,14150,14190,14000,70,16760,16700,16640,16600,16670,16460,80,19800,19800,19800,19800,19800,19800,90,24010,24010,24010,24010,24010,240110,28860,28860,28860,28860,28860,28861,10,33960,33960,33960,33960,33960,33961,30,45150,45150,45150,45150,45150,45151,50,46590,46590,46590,46590,46590,465920,47110,47110,47110,47110,47110,47112,50,46080,46080,46080,46080,46080,460830,45740,45740,45740,45740,45740,45743,50,44280,44280,44280,44280,44280,442840,42230,42230,42230,42230,42230,42234,50,41720,41720,41720,41720,41720,417250,40690,40690,40690,40690,40690,4069

Рисунок 6 - Зависимость коэффициента сопротивления давления РН от числа Маха и высоты

Так как на рисунке 6 невозможно различить зависимости при , на рисунке 7 приведен график для этих зависимостей в увеличенном масштабе.


Рисунок 7 - Зависимость коэффициента сопротивления давления РН от числа Маха и высоты при


Расчет коэффициента продольной силы при нулевом угле атаки

Коэффициент продольной силы летательного аппарата при нулевом угле атаки определяется как сумма коэффициентов сопротивления трения летательного аппарата и коэффициента сопротивления давления летательного аппарата при нулевом угле атаки



Коэффициент продольной силы является функцией числа Маха набегающего потока воздуха. При нулевом угле атаки значения коэффициентов продольной силы и лобового сопротивления совпадают [1].

Значения коэффициента продольной силы РН при нулевом угле атаки в диапазоне высот от 0 до 60 км представлены в таблице 8.

Таблица 8 - Коэффициент продольной силы РН

MH, км010203040600,10,20600,20770,22490,19890,28930,48450,30,19890,20730,21190,23300,23250,34410,50,18400,19140,19480,21270,20330,28900,70,20610,21290,22400,23160,21850,29020,80,23680,24340,25410,26110,24690,31380,90,27660,28340,29440,30130,28560,350510,32420,33080,34150,34790,33170,39311,10,37420,38070,39100,39710,38050,43891,30,48450,49060,50020,50570,53790,54241,50,49730,50310,51220,51710,54690,550120,49880,50370,51160,51560,54050,54282,50,48510,48950,49630,50560,52060,523730,47890,48270,48860,49730,51600,51383,50,46170,46500,47020,47780,49400,493940,43910,44200,44660,45320,46730,46924,50,43210,43470,43870,44460,45700,460650,42030,42260,42610,43130,44230,4473

Рисунок 8 - Зависимость коэффициента осевой силы от числа Маха и высоты

6. Производная по углу атаки коэффициента аэродинамической нормальной силы ракеты-носителя


Состав производной по углу атаки коэффициента нормальной силы летательного аппарата

Производная коэффициента нормальной силы по углу атаки для всего летательного аппарата определяется через его составляющие для отдельных частей аппарата. Однако в данном случае летательный аппарат представляет собой только корпус и коэффициент продольной силы давления



Производная по углу атаки коэффициента нормальной силы изолированного корпуса

Производная коэффициента нормальной силы изолированного корпуса по углу атаки C?y корп определяется для типичной схемы корпуса ракеты-носителя, состоящей из носовой части с продолжающейся цилиндрической частью и кормовой части по рисунку 19 [2]. Согласно теории тонких осесимметричных тел, нормальная сила появляется только на участках корпуса с переменной площадью поперечного сечения S = S(x1) и направление действия этой силы определяется знаком производной площади поперечного сечения dS(x1)/dx1 по координате x1 вдоль продольной оси системы координат Ox1y1z1, ось Ox1 которой противоположна оси Ox связанной системы координат Oxyz и направлена в сторону хвостовой части корпуса. Расширяющиеся части корпуса (в данном случае юбка РН), которым соответствует положительный знак производной, создают положительную нормальную силу, направленная в сторону оси Oy связанной системы координат, которая, в свою очередь, направлена в сторону увеличения угла атаки. Сужающиеся части корпуса, которым соответствует отрицательный знак производной, имеют отрицательную нормальную силу. При этом нормальная сила цилиндрических частей с постоянной площадью поперечного сечения равна нулю.

Однако эксперименты показывают, что при больших дозвуковых и при сверхзвуковых скоростях цилиндрические части корпуса, примыкающие к его расширяющимся частям, также создают некоторую нормальную силу. Этот факт отображен в экспериментальных зависимостях, представленных на рисунке 20 [2].

Величина производной C?y корп зависит от формы корпуса и определяется как сумма производных коэффициентов нормальных сил по углу атаки отдельных частей корпуса, умноженных на отношения характерных площадей этих частей к площади миделя корпуса



где Cy?н1, Cy?н2 - производные коэффициентов нормальных сил, действующих на носовую и переходную часть корпуса, 1/град;

S1, S2 - характерные площади оснований конических и кормовой частей, м2.

Расчет производной Cy?н2 и Cy?н3 для переходной части имеет вид



Производные Cy?н2 и Cy?н2 определяются по рисунку 20 [2]. Производная Cy?н2 находится для псевдоконуса при условии отсутствия цилиндрической части за этим конусом.

Результаты расчетов по определению производной по углу атаки коэффициента аэродинамической нормальной ракетоносителя приведены в таблице 9.

Таблица 9 - Производная коэффициента аэродинамической нормальной и подъемной силы летательного аппарата на высоте h = 10 км

MC?у н1, 1/град, 1/град, 1/град, 1/град, 1/град0,10,03400,03500,03480,00380,03430,30,03420,03530,03480,00410,03450,50,03440,03540,03480,00420,03490,70,03450,03550,03480,00430,03530,80,03470,03650,03480,00500,03750,90,03500,03700,03480,00580,041110,03700,04200,03480,01080,04401,10,04200,04650,03480,01530,05301,30,04550,04610,03480,01510,05581,50,04500,04500,03480,01380,054220,04100,04150,03480,01180,04602,50,03680,04100,03480,00940,042430,03600,04000,03480,00880,04113,50,03560,03900,03480,00780,039840,03540,03800,03480,00660,03804,50,03520,03750,03480,00620,037450,03500,03700,03480,00570,0371

Расчет производной коэффициента аэродинамической подъемной силы летательного аппарата по углу атаки

Коэффициент лобового сопротивления при нулевом угле атаки Cxa0(M?) и производная коэффициента нормальной силы летательного аппарата по углу атаки C?y позволяют определить продольную Х и нормальную Y силы, действующие на летательный аппарат



При нулевом угле скольжения ? = 0 подъемная сила выражается через нормальную и продольную силы. Для малых углов атаки sin? ? ?, cos? ? 1 выражение примет вид


где ? - угол атаки, выраженный в рад.

Представляем подъемную силу Ya через производную коэффициента подъемной силы по углу атаки



В данной работе расчеты производной коэффициента нормальной силы по углу атаки имеют размерность 1/град, в этом случае производная C?ya будет определяться по формуле



График зависимости производной коэффициента подъемной силы от числа Маха в нормальной и скоростной системе координат изображен на рисунке 9.


Рисунок 9 - Зависимость производной коэффициента нормальной силы по углу атаки в с коростной и нормальной системе координат от числа Маха


7. Коэффициент силы лобового сопротивления ракеты-носителя


Состав коэффициента силы лобового сопротивления летательного аппарата

Коэффициент лобового сопротивления летательного аппарата представляется в виде суммы двух коэффициентов - коэффициента лобового сопротивления при нулевом угле атаки Cxa0, зависящего от числа Маха, и коэффициента индуктивного сопротивления Cxai, зависящий от числа Маха и угла атаки. Поэтому коэффициент лобового сопротивления является функцией числа маха и угла атаки



Коэффициент индуктивного сопротивления летательного аппарата, состоящего только из корпуса, определяется следующим образом



где Cxai, Cxai корп - коэффициенты индуктивного сопротивления соответственно летательного аппарата и изолированного корпуса;

Коэффициент индуктивного сопротивления корпуса

В диапазоне малых углов атаки коэффициент индуктивного сопротивления корпуса определяется по формуле



в которой производная C?y корп имеет размерность 1/град, а угол ? град.

Коэффициент, учитывающий влияние перераспределения давления по расширяющимся частям корпуса летательного аппарата, определяется следующим образом



где ?н1, ?н2 - коэффициенты перераспределения давления на носовой, переходной и кормовой частях ЛА;

S1, S2 - характерные площади оснований конических частей, м2.

Коэффициент ?н2 для переходной части корпуса, представляющей собой усеченный конус, определяется по формуле



где ?н2 - коэффициент, учитывающий перераспределение давления по конической носовой части продленного корпуса длиной Lн2.

Коэффициенты ?н1, ?н2 определяются по рисунку 30 [2].

Результаты расчетов по определению коэффициента индуктивного сопротивления ракетоносителя приведены в таблице 10, и на рисунке 10.


Таблица 10 - Коэффициент индуктивного сопротивления корпуса

M?н1? н2? кСxai? = 2o? = 4o? = 6o? = 8o0,1-0,3000-0,0208-0,28900,00170,00680,01520,02700,3-0,3000-0,0208-0,28900,00170,00680,01530,02730,5-0,3000-0,0208-0,28900,00170,00690,01560,02770,7-0,3000-0,0208-0,28900,00180,00700,01580,02820,8-0,3000-0,0208-0,28000,00190,00770,01740,03100,9-0,2900-0,0197-0,25200,00230,00900,02030,036110,0200-0,01760,00030,00310,01230,02760,04911,10,3700-0,00620,32500,00450,01800,04040,07181,30,5000-0,00100,43900,00500,01990,04470,07941,50,57000,00620,51300,00500,02010,04530,080520,70000,02080,64200,00480,01910,04300,07642,50,72000,02490,69000,00460,01860,04180,074330,77000,03010,72000,00460,01850,04160,07403,50,78000,03110,72900,00460,01820,04100,072940,78100,03110,73200,00440,01770,03990,07104,50,78200,03110,73500,00440,01760,03960,070450,79000,03110,73900,00440,01760,03950,0702

Рисунок 10 - Зависимость индуктивного сопротивления от угла атаки и от числа Маха


Результаты расчетов по определению коэффициента силы лобового сопротивления ракетоносителя приведены в таблице 11, и на рисунке 11.

Таблица 11 - Коэффициент лобового сопротивления корпуса

xa? = 2o? = 4o? = 6o? = 8o0,10,20930,21440,22290,23470,30,20900,21410,22260,23450,50,19310,19830,20700,21910,70,21460,21990,22870,24100,80,24530,25110,26080,27430,90,28570,29240,30370,319510,33390,34310,35850,38001,10,38520,39860,42110,45251,30,49550,51040,53520,57001,50,50810,52320,54840,583620,50850,52280,54670,58012,50,49410,50800,53120,563730,48730,50110,52430,55663,50,46960,48320,50600,537940,44640,45970,48190,51304,50,43910,45230,47430,505150,42690,44010,46210,4928

Рисунок 11-Зависимость коэффициента лобового сопротивления ЛА от числа Маха


8. Координата фокуса ракеты-носителя


Определение фокуса летательного аппарата

Фокусом летательного аппарата по углу атаки называют точку на продольной оси, характерную тем, что момент аэродинамических сил относительно поперечной оси Oz, проходящей через эту точку, не зависит от угла атаки. Для ракет-носителей, обладающих продольной плоскостью симметрии и нулевой подъемной силой при нулевом угле атаки, в диапазоне малых углов атаки, для которых нормальная сила линейно зависит от величины угла атаки, фокус является точкой приложений этой нормальной силы. Наличие координаты фокуса и координаты центра масс позволяет определить величину статической устойчивости летательного аппарата.

Состав координаты фокуса летательного аппарата

Координата фокуса летательного аппарата xF определяется по оси Ox1 связанной с корпусом системы координат Ox1y1z1, ось Ox1 которой противоположна оси Ox стандартной связанной системы координат Оxyz и направлена в сторону хвостовой части корпуса, ось Oy1 направлена параллельно оси Oy стандартной связанной системы координат, а ось Oz1 дополняет систему до правой. При ненулевом угле атаки ? ? 0 на корпус действует нормальная сила Yкорп., которая имеет свой центр давления, совпадающий с фокусом и определяемый координатой XF корп. Для ЛА, состоящего только из корпуса



Координата фокуса корпуса летательного аппарата

Координата фокуса корпуса ЛА определяется по формуле



Координата фокуса комбинации носовой оживальной части с цилиндром определяется по теории удлиненных тел с учетом эмпирических поправок

Координата фокуса комбинации носовой оживальной части с цилиндром определяется по теории удлиненных тел с учетом эмпирических поправок


где Lн1 - длина оживальной носовой части, м;

Wн1 - объем оживальной носовой части, м3;

S1 - площадь основания оживальной носовой части, м2;

?XF н - относительное смещение фокуса носовой оживальной части из-за наличия цилиндрической части корпуса, примыкающей к носовой.

Относительная величина смещения фокуса ?XF н зависит от числа Маха М?, удлинений носовой ?н1 и цилиндрической ?ц1 частей и определяется по эмпирической зависимости [2, рисунок 34].

Координата фокуса комбинации усеченного конуса и примыкающего к нему цилиндра определяется следующим образом. Усеченный конус достраивается до полного. Обозначая через X`F Н2 и LН2 координату фокуса и длину полной конической носовой части и соответственно через X`»F Н2 и L»Н2 координату фокуса и длину фиктивного конуса, определяются моменты относительно вершины фиктивного конуса от действия нормальных сил YН2, Y»Н2 и YН2, действующие на полный конус, фиктивный конус и усеченный конус соответственно



Выражая нормальные силы через коэффициенты



координата фокуса усеченного конуса переходной части определяется следующим образом


В этом выражении координата фокуса полного конуса X`F Н2, за которым следует цилиндрическая часть, определяется по формуле



где WH2 - объем полного конуса, м3;

?XFH - относительное смещение фокуса за счет влияния цилиндрической части аппарата;

Lн2 - высота продленного конуса.

Координата фокуса фиктивного конуса X`»F Н2 за которым отсутствует цилиндрическая часть, влияющая на смещение фокуса, определяется по формуле



где W»H2 - объем фиктивного конуса, м3;

Lн2 - высота фиктивного конуса.

Координата фокуса переходной части относительно носка летательного аппарата находится с учетом расстояния вершины фиктивного конуса от носка летательного аппарата А2



Результаты расчетов по определению координаты фокуса ракетоносителя приведены в таблице 12, и на рисунке 12.


Таблица 12 - Определение координаты фокуса ЛА

M?xFн1?x'Fн2X' Fн2, мX'' Fн2, мX` Fн2, мX Fн2, мX Fн1, мX F, м0,10,05000,05008,04327,084015,94008,29501,34092,11140,30,05100,05008,04327,084015,55197,90691,34322,12150,50,05300,05008,04327,084015,18807,54301,34782,09840,70,05300,05008,04327,084014,98697,34191,34782,07700,80,05700,05608,11057,084015,01897,37391,35702,10860,90,06000,07008,26767,084014,61946,97441,36392,028410,06200,09008,49217,084012,55394,90891,36852,23901,10,07000,11008,71657,084012,04184,39681,38692,25741,30,07700,12408,87377,084012,45944,81441,40302,32811,50,08600,14009,05327,084013,51145,86641,42372,553020,08700,12808,91857,084012,64254,99751,42602,42112,50,06300,11808,80637,084014,89067,24561,37082,674330,04900,10008,60437,084013,98536,34031,33862,41013,50,04750,09788,57967,084014,55096,90591,33522,425740,04700,09738,57407,084015,86598,22091,33402,54384,50,04700,09718,57187,084016,21018,56511,33402,545150,04700,09708,57067,084016,84339,19831,33402,5436

Для оценки центра масс ракету необходимо представить как летающую цистерну. Переходные участки не используются для расположения топлива.

Координата центра определяется по формуле



где - вес топлива первого и второго цилиндрических участков, Н;

- вес РН, Н;

- координаты ЦМ участков, м.

Вес топлива определяется по формуле


где - удельный вес топлива, Н/м3;

- объем топлива в баках n ступени, м3.

В качестве горючего в ракете используется полиуретан с , в качестве окислителя - перохлорат аммония с . Для приблизительной оценки положения ЦМ примем осредненное значение удельного веса .

Определим вес топлива

Определим координату ЦМ ЛА

Координата центра масс изображена на рисунке 12.


Рисунок 12-Зависимость изменения координаты фокуса от числа Маха


Заключение


В результате проведенных расчетов для РН «Посейдон C-3» получены аэродинамические характеристики аппарата. Исходя из полученных результатов можно судить, что коэффициент лобового сопротивления достигает своего максимума (Сxa=0,5836, ? = 8o) при значении числа Маха M? = 1,5 после чего начинает убывать. В свою очередь производная по углу атаки коэффициента подъемной силы будет иметь наибольшее значение (=0,0558) при значении числа Маха M? = 1,3. Координата фокуса ракеты-носителя лежит в диапазоне [2,1; 2,7] м в зависимости от чисел Маха. После проведенных расчетов можно сделать вывод о том, что данный летательный аппарат неустойчив (координата фокуса находится перед координатой центра масс летательного аппарата). Для обеспечения устойчивости летательного аппарата необходимо установить крылья, которые смещали бы фокус назад.

Полученные аэродинамические характеристики позволяют определить динамическое воздействие внешней среды на летательный аппарат на активном участке его полета в атмосфере.


Список использованных источников

аэродинамический летательный мах

1.#"justify">2.Расчет аэродинамических характеристик летательных аппаратов [Электронный ресурс]: электрон. учеб пособие /В.В. Васильев, А.Н. Никитин, В.А. Фролов, В.Г. Шахов; Минобрнауки России, Самар. гос. аэрокосм. ун-т им. С.П. Королева (нац. исслед. ун-т). - Электрон. текстовые и граф. дан. (2,315 Мбайт). - Самара, 2012. - 1 эл. опт. диск (CD-ROM).


Аэродинамические характеристики ракеты-носителя «Посейдон С-3» Введение Ракета-носит

Больше работ по теме:

КОНТАКТНЫЙ EMAIL: [email protected]

Скачать реферат © 2017 | Пользовательское соглашение

Скачать      Реферат

ПРОФЕССИОНАЛЬНАЯ ПОМОЩЬ СТУДЕНТАМ