Размещено на http://www.allbest.ru/

Расчет аэродинамических характеристик летательного аппарата AGM-84B “Harpoon”

Содержание

аэродинамический летательный лобовое сопротивление тангаж

Введение

Подъёмная сила

Определение коэффициента

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

Лобовое сопротивление

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

Коэффициент лобового сопротивления несущих поверхностей при

Момент тангажа

Момент тангажа при

Введение

Аэродинамический расчет является важнейшим элементом аэродинамического исследования ЛА или его отдельных частей (корпуса, крыльев, оперения, управляющих устройств). Результаты такого расчета используются при траекторных вычислениях, при решении задач, связанных с прочностью движущихся объектов, при определении летно-технических характеристик ЛА.

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

Аэродинамические силы и моменты можно определить с использованием аэродинамических коэффициентов.

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

Приведенная методика определения аэродинамических характеристик является приближенной. На рисунке приведена схема ракеты, а все размеры в табл. 3.1, здесь L - длина ЛА, D - диаметр корпуса ЛА, , - длины носовой и цилиндрической частей. , - удлинение носовой и цилиндрической частей корпуса, l - полный размах несущей поверхности,- размах несущей поверхности без учета диаметра корпуса, b , b1 - длины бортовой и концевой хорд консоли, - площади двух консолей передней (задней) несущей поверхности - удлинение передних (задних) консолей, - сужение передних (задних) консолей, - угол стреловидности консолей, измеряемый по линии середины хорд,- угол стреловидности консолей, измеряемый по передней (задней) кромке, - относительная толщина профиля,- длина средней аэродинамической хорды (САХ) крыльев, - координата z средней аэродинамической хорды крыльев,- координата x средней аэродинамической хорды крыльев, - расстояние от передней точки корпуса до передней (задней) консоли.

Для исследования динамики ЛА необходимо учесть действующие силы и моменты, в том числе и аэродинамические. Полную аэродинамическую силу , зависящую от ряда факторов можно представить составляющими по скоростным осям координат (x, y, z) или по связанным (), а полный аэродинамический момент М - разложенным по осям (). В случае симметричного ЛА подъемная сила Y и боковая сила Z имеют одинаковую зависимости соответственно от углов атаки и скольжения , от углов отклонения рулей и .

1. Подъемная сила

Подъемная сила определяется по формуле

,

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

Коэффициент принято определять в скоростной системе координат 0xyz. Наряду с коэффициентом далее рассматривается и коэффициент нормальной силы , определяется в связанной системе координат .

Эти коэффициенты связаны между собой соотношением

. (1)

Представляем ЛА в виде совокупности следующих основных частей: корпуса (фюзеляжа), передних (I) и задних (II) несущих поверхностей. При небольших углах атаки и углах отклонения несущих поверхностей зависимости и близки к линейным, т. е. могут быть представлены в виде

;

;

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

Значения и у беспилотных ЛА в большинстве случаев близки к нулю, поэтому в дальнейшем они не рассматриваются. В качестве органов управления принимаются задние несущие поверхности.

2. Определение коэффициента

Из варианта (I) найдем производную :

Таблица 1

Величина	Значение	Величина	Значение
	3.79		1.932
	0.315		0.798
	0.38		0.273
	3.41		0.147
	1.21		0.1014
	10.83		2.3007
	0.915		1.8571
	0.6		0.2162
	0.735		0.1087
	0.336		3.4449
	0.3213		3.36
	1.12		27
	2.1875		0.5205
	0.561		0.1
	0.1314
	2.0886
	28o
	0.5317
	0.1

. (2)

При малых углах атаки и при можно положить , тогда равенство (2) принимает вид . Представим нормальную силу ЛА в виде суммы трех слагаемых

, (3)

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

;

;

.

Поделив равенство (3) почленно на и изъяв производную по , получим в точке 0

, (4)

где ; - коэффициенты торможения потока;; ; - относительные площади частей ЛА. Рассмотрим подробнее величины, входящие в правую часть равенства (4).

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

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

Величина этой силы выражается через нормальную силу изолированных крыльев (т. е. крыльев, составленных из двух консолей) с помощью коэффициента интерференции k : . Величины и kI подсчитываются при числе Маха .

Третье слагаемое в выражении (4) аналогично второму. Единственное отличие состоит в том, что при определении угла атаки задней несущей поверхности надо учитывать средний угол скоса потока , вызываемого передней несущей поверхностью: . При малых углах атаки зависимость близка к линейной. В том случае и производную можно выразить в виде

. (5)

Все величины, входящие в (5), подсчитываются при числе Маха .

Расчет и содержится в табл. 3.2 и табл. 3.4, где L1 - расстояние от носика корпуса до середины бортовой хорды консолей, Lx - длина хвостовой части корпуса (от конца бортовой хорды до кормового среза корпуса).

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

продифференцируем выражение (1) по углу II :

.

Таблица 2

Величина	M
	0.3	0.95
	0.788	0.548
	0.035	0.048
	1	1
	1.796	1.796
k	0.99	0.965
	0.316	0.316
	0.9924	0.9734
	1	1
	0.9924	0.9734
	0.3	1.2
	0.5955	0.5955
	0.5199	0.5199
	1.068	0.7429
	0.0245	0.028
	0.0274	0.0314
	0.344	0.344
	2.1875	2.1875
	1.842	1.842
	1.37	1.37
	2.2995	2.2995
	102.084	408.336
	0.000014607
	0.0279	0.0313
ж	0.987	0.9856
	7.3	7.3
	0.9896	0.9896
ж	1	1
	1	1
	1.7991	1.7966
	1.3381	1.3362
	0.8	0.785
	0.3975	0.393
	0.3947	0.3947
	0.5385	0.5385
	0.9962	0.985
	Лежат в одной плоскости
	0	0
i	3	3
	1	1
	0.2989	1.1839
	0.3947	0.3947
	1.8571	1.8571
	2.06	2.06
	1.42	1.42
	3.4965	3.4965
	101.7097	402.86
	0.039	0.0439
	0.978	0.9754
	1	1
	11.1	11.1
	0.9984	0.9984
	0	0
	1.7986	1.7938
	1.3377	1.3342
	0.8677	0.9035
	1	1
	4.1245	4.1245
	0.0779	0.0779
	1.3017	1.3017
	0.035	0.048
	0.0493	0.0564
	1.1975	1.1975
	2.1947	2.2542
	1.8571	1.8571
	1.0679	1.0679
	0.02	0.0207
	0.046	0.0476
	0.0109	0.0082
(I)	0.035	0.048
(II)	0.2033	0.2326
(III)	0.0637	0.0715
	0.302	0.3521

При малых углах и это выражение принимает следующий вид :

.

Поделив равенство (3) почленно на qS и взяв производную по , получим

.

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

. Тогда .

Расчет представлен в табл. 3.3, где - угол стреловидности оперения;- коэффициент снижения подъёмной силы из-за щели между рулём и корпусом при отклонении рулей.

Таблица 3

Величина	М=0.3	M=0.95
	1.42	1.42
	0.9788	0.9788
	0.039	0.0443
	0.9794	0.9751
	1	1
	1	1
	1.3907	1.3846
	27.5	27.5
	0.887	0.887
	0.85	0.85
	0.7225	0.7225
	0.0597	0.0603

Таблица 4

Величина	М=0.3	M=0.95
	0.0779	0.0779
	4.2013	4.2013
	3.79	3.79
	0	0
	0.0052	0.0042
	0.98	0.9
	0.0051	0.0038
	0.1402	0.1133
	1.21	1.21
	-0.005	0.24
,м2	0.039	0.039
	0.5	0.5
	0.125	0.275
	0.7071	1
	0.0442	0.1375
	0.1794	0.4908
	0.007	0.0057
	0.98	0.925
	0.0069	0.0053
	0.1	0.1
	1.28	1.28
	0.0088	0.0068
	1.068	0.7429
	1.12	1.12
	0.5199	0.5199
	0.0088	0.0068
	101.71	402.86
	0.2162	0.2162
	0.0087	0.0065
	0.98	0.925
	0.0085	0.006
	1.28	1.28
	0.0109	0.0077
	0.0109	0.0077
	8.249	8.249
	2.6033	2.6033
	0.0726	0.0561
	0.0282	0.0195
	0.2942	0.5947
	0.2969	0.3417

4. Лобовое сопротивление

Сила лобового сопротивления подсчитывается по формуле

.

Коэффициент лобового сопротивления ЛА представим в виде суммы двух слагаемых , где -коэффициент сопротивления при ; - коэффициент индуктивного сопротивления , под которым понимается сопротивление , зависящее от углов , и . Коэффициент ЛА может быть выражен в виде

,

где 1.05 - поправка на неучтенные детали; - отношение суммарной площади всех консолей передней несущей поверхности к характерной площади; - то же для задней несущей поверхности; ,,- коэффициенты изолированных частей ЛА.

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

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

,

где последние три слагаемых представляют собой сопротивление давления.

6. Коэффициент лобового сопротивления несущих поверхностей при

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

Лобовое сопротивление несущей поверхности с заостренными задними кромками при складывается из профильного и волнового сопротивления . В соответствии с этим можно написать

.

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

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

У ЛА с крестообразным расположением крыльев (++) сила лобового сопротивления создается двумя парами передних и задних несущих поверхностей , поэтому коэффициенты и должны перемножаться на соответствующие удвоенные безразмерные площади . Расчет смотри в табл.3.4.

7. Момент тангажа

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

.

Величина аэродинамического момента при данной скорости и высоте полета зависит от ряда факторов и прежде всего от угла атаки и углов отклонения органов управления. Кроме того, на величину момента влияет угловая скорость вращения ЛА , а также скорости изменения угла атаки и отклонения рулей, характеризуемые производными и . Таким образом,

. (6)

При малых значениях аргументов выражение (6) можно представить в виде линейной функции

, (7)

где и т.д. - частные производные момента тангажа по соответствующим параметрам.

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

(7)

Для упрощения записи величин, входящих в выражения (6) и (7),индекс “I” будем в дальнейшем опускать. Кроме того, будем опускать черточки в обозначениях частных производных

8. Момент тангажа при

Рассмотрим величину аэродинамического продольного момента, действующего на ЛА, при условии, что угловая скорость , а угол атаки и углы отклонения органов управления остаются неизменными по времени.

Введем понятие центра давления ЛА. Центр давления - это точка на продольной оси 0x1 , через которую проходит равнодействующая - аэродинамических сил.

Момент аэродинамических сил относительно центра давления можно выразить в виде , а коэффициент момента

,

здесь - координата центра тяжести ЛА, - координата центра давления (отчет производится от носа корпуса).

По аналогии с понятием центра давления всего ЛА введем также понятие центров давления его частей как точек приложения нормальных сил, создаваемых этими частями.

Из условия равновесия имеем

.

Отсюда находим выражение для :

При малых углах атаки и углах отклонения рулей удобно пользоваться понятием аэродинамических фокусов ЛА. Фокусом ЛА по углу атаки называется точка приложения той части нормальной силы, которая пропорциональна углу атаки (т.е. ). Тогда при закрепленных органах управления момент аэродинамических сил относительно оси 0z1 , проходящей через точку фокуса, не зависит от угла атаки. Аналогично можно показать, что момент относительно фокуса по не зависит от , а момент относительно фокуса по не зависит от .

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

,

где

,

,

.

В этих выражениях , , - координаты фокусов по , и .

Расчет содержится в табл.3.5, а расчет - в табл.3.6.

Таблица 5

Величина	М=0.3	М=0.95
	0.38	0.38
	0.0099	0.0099
	0.788	0.548
	8.9504	8.9504
	0.04	0.2
	0.0152	0.076
	0.2681	0.3289
	0.0779	0.0779
	0.2681	0.3289
	2.0886	2.0886
	0.561	0.561
	0.22	0.38
	0.344	0.344
	0.032	0.032
	0.0096	0.0096
	0.5317	0.5317
	2.2069	2.2967
	2.212	2.3018
	1.932	1.932
	0.735	0.735
	0.5955	0.5955
	0.2319	0.3919
	1	1
	1	1
	2.1024	2.22
	2.183	2.2799
	3.4449	3.4449
	0.2162	0.2162
	0.2	0.37
	3.4881	3.5249
	0.5205	0.5205
	0.3947	0.3947
	0.483	0.483
	0.0285	0.0285
	0.0069	0.0069
	3.4845	3.5213
	1.1975	1.1975
	0.224	0.394
	3.36	3.36
	0.273	0.273
	3.4212	3.4676
	3.4703	3.5095
(I)	0.0094	0.0158
(II)	0.4438	0.5303
(III)	0.2211	0.2509
	2.2328	2.2636
	1.6643	1.6951
	-0.0453	-0.0528

Таблица 6

Величина	М=0.3	М=0.95
	0.9586	0.9544
	1.3907	1.3846
	3.4673	3.5071
	-0.0284	-0.0288

9. Момент тангажа, вызванный вращением ЛА вокруг оси Z

Рассмотрим ЛА, летящий со скоростью v и одновременно вращающийся вокруг своей оси (поперечной) с угловой скоростью .

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

Величина очень мала и в расчётах подъёмной силы ею обычно пренебрегают.

Момент существенно влияет на динамические свойства ЛА. Он называется демпфирующим моментом тангажа или продольным демпфирующим моментом.

Величина демпфирующего момента пропорциональна угловой скорости . Поэтому .

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

Представим продольный демпфирующий момент как сумму моментов, создаваемых частями ЛА: . Это выражение можно переписать в соответствии с равенством (9):

.

Сокращая на , получаем:

Где , , .

Расчёт представлен в таблице 7.

Таблица 7

Величина	М=0.3	М=0.95
	3.79	3.79
	12.0317	12.0317
	1.21	1.21
	0.5337	0.5337
	0.4391	0.4473
	-0.44	-0.4382
	2.0886	2.0886
	-0.5025	-0.5177
	0.218	0.218
	0.0594	0.0594
	0.7623	0.7623
	0.013	0.013
	0.0503	0.0477
	0.8679	0.7903
	0.2392	0.2392
	0.7923	0.7923
	0.0519	0.0519
	0.0137	0.0137
	-8.2359	-8.0934
	47.79	46.2412
	0.3382	0.3325
	-14.621	-16.2592
	-0.904	-0.9301
	-26.2986	-29.1658
	-13.1493	-14.5829
	-39.4479	-43.7487
	0.148	0.148
	0.06	0.06
	-0.0817	-0.084
	-0.1841	-0.2023
	-0.7058	-0.7245

Список используемой литературы:

“Динамика полёта” А.А. Лебедев , Л.С. Чернобровкин

Учебное пособие “Траектории. Аэродинамика. Характеристики ЛА”

Internet:

http://www.fas.org/man/dod-101/sys/smart/agm-84.htm

www.danshistory.com

http://faq.guns.ru/agm84d.html

www.time.com/time/daily/iraq/weapons/messiles

http://sterett.org/missile

Размещено на Allbest.ru