Составление математической модели теоретических обводов крыла - конспект -  Астрономия, Конспект из Астрономия
filizia
filizia11 June 2013

Составление математической модели теоретических обводов крыла - конспект - Астрономия, Конспект из Астрономия

PDF (605.7 KB)
16 страница
526количество посещений
Описание
Rybinsk State Academy of Aviational Technology. Лекции и рефераты по Астрономии. Классификация несущих поверхностей Все многообразие проектируемых несущих поверхностей можно классифицировать следующим образом: линейча...
20очки
пункты необходимо загрузить
этот документ
скачать документ
предварительный показ3 страница / 16
это только предварительный показ
консультироваться и скачать документ
это только предварительный показ
консультироваться и скачать документ
предварительный показ закончен
консультироваться и скачать документ
это только предварительный показ
консультироваться и скачать документ
это только предварительный показ
консультироваться и скачать документ
предварительный показ закончен
консультироваться и скачать документ

8 Составление математической модели теоретических обводов крыла

8.1 Классификация несущих поверхностей

Все многообразие проектируемых несущих поверхностей можно

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

существенно нелинейчатые и интегральные. При этом в основу

классификации положен скорее не теоретический подход, а некоторая

практическая характеристика, которую можно назвать

геометротехнологической.

Чаще всего, как и в нашем случае, при математическом описании

несущих поверхностей применяются линейчатые поверхности, образованные

путем перемещения прямолинейной образующей по двум криволинейным

пространственным направляющим.

Для однозначности определения положения в пространстве

прямолинейной образующей необходимо задать закон ее перемещения. Этот

закон может задаваться в виде направления в пространстве, например,

плоскости параллелизма, или третьей направляющей. Если же две

направляющие являются плоскими кривыми, лежащими в параллельных

плоскостях, а в качестве третьей направляющей выбрана прямая,

параллельная указанным плоскостям, то образованная в этом случае

поверхность будет называться поверхностью с пропорциональной разбивкой.

Действительно, если мы рассмотрим проекцию направляющих и образующих

на плоскость, перпендикулярную прямолинейной образующей, нетрудно

видеть, что отрезки проекций криволинейных направляющих, отсекаемые

образующими, будут пропорциональны.

Иногда такие поверхности называют линейчатыми поверхностями с

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 64

docsity.com

процентной разбивкой.

Частным случаем линейчатых поверхностей является развертываемая

поверхность, отличающаяся тем, что прямолинейная образующая,

соединяющая две точки на направляющих, и касательные в них

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

направляющих. Исходя из этого углы наклона касательных, как в начале, так

и в конце используемых отрезков двух кривых должны быть равны между

собой, а изменение углов наклона вдоль кривых должно быть гладким и

непрерывным.

Нелинейчатой будем называть такую поверхность, у которой способ

перехода от сечения к сечению в параллельных плоскостях не обеспечивает

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

линиями. Существенно нелинейчатая поверхность — это поверхность такого

крыла, геометрические параметры которого (форма профиля, толщина и

вогнутость его средней линии и другие) значительно изменяются вдоль

размаха крыла. Кроме того, форма в плане описывается криволинейными

передней и задней кромками. Следует отметить, что такая поверхность

позволяет существенно повысить аэродинамические характеристики крыла.

Примером такого крыла является, например, крыло сверхзвукового

пассажирского самолета Ту-144.

Для получения дополнительного выигрыша в аэродинамических

характеристиках (интерференция) фюзеляжа и крыла, а также бол ее полного

использования компонуемого объема в последние годы получило широкое

распространение объединение поверхности фюзеляжа и крыла в гладкую

единую поверхность с обеспечением их плавного сопряжения. Такая

компоновка крыла и фюзеляжа получила название интегральной. Это

решение реализовано при проектировании американского самолета В1-А.

8.2 Основные геометрические характеристики крыла

Геометрические характеристики крыла в основном можно определить

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 65

docsity.com

по его форме в плане. Вообще говоря, хорда крыла определяется как

условная линия, соединяющая точки передней и задней кромок крыла,

полученные в результате их пересечения плоскостью, параллельной

плоскости симметрии самолета.

Хорда, взятая в произвольном по размаху месте крыла, называется

местной, ее длина )'(zb равна:

   2кпкз2кпкз zyzyzxzxzb )()()()()( ....  , (8.1) где )(. zx кп , )(. zy кп - координаты передней кромки крыла; )(. zx кз ,

)(. zy кз - координаты задней кромки крыла.

Хорда, определяемая при z = 0 в системе координат самолета,

называется центральной (корневой) 0b • Бортовая хорда — это хорда крыла

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

треугольного в плане крыла концевая хорда вырождается в ноль.

При расчете аэродинамических характеристик крыла чаще пользуются

геометрическими параметрами его проекции на базовую или строительную

плоскости.

Базовая плоскость крыла (БПК) — это плоскость, перпендикулярная

плоскости симметрии самолета и проходящая через корневую хорду крыла.

Строительной плоскостью крыла (СПК) называют плоскость,

проходящую через хорду одного из сечений крыла (чаще всего корневого

или бортового) и точку, лежащую на хорде концевого сечения. Тогда при

прямолинейности задней кромки крыла СПК будет определяться двумя

пересекающимися линиями — корневой хордой и задней кромкой крыла.

При нулевом значении угла поперечного V крыла базовая и

строительная плоскости совпадают. Поэтому будем считать, что крыло в

плане ограничено проекциями линий передней и задней кромок на СПК,

корневой и концевой .хордами. Площадь, ограниченную этими линиями,

будем называть проекционной площадью крыла.

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 66

docsity.com

Местный угол стреловидности передней кромки крыла кп. - угол

между касательной к линии передней кромки в заданной точке и плоскостью,

перпендикулярной к корневой хорде крыла. Аналогично определяется угол

стреловидности задней кромки кз. и линии четверти хорды крыла.

Удлинение крыла определяется как отношение квадрата полного

размаха l к его площади S :

S l2  . (8.2)

Другой важной характеристикой формы крыла в плане является

сужение, которое определяется как отношение корневой хорды и концевой:

к

0

b b  (8.3)

В ряде случаев из конструктивных соображений или по

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

обрезают. В этом случае для определения сужения крыла исходную форму в

плане заменяют фиктивным трапециевидным крылом равной площади с

совпадающими передней и задней кромками. Концевая хорда кb' такого

крыла определяется из условия равенства площадей по формуле:

    

  

 

   

    кзкп0кзкп10к tgtg2

1btgtgl 2 1bb ....' , (8.3)

а сужение крыла определяется так:

к

0

b b '

 . (8.4)

При этом следует отметить, что полученное фиктивное крыло нельзя

использовать для расчета таких характеристик, как удлинение и средняя

аэродинамическая хорда.

Средняя аэродинамическая хорда (САХ) определяется как хорда

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 67

docsity.com

прямоугольного крыла, равного по размаху и площади исходному. САХ

является одним из важнейших геометрических параметров несущей

поверхности, используемых при расчетах аэродинамических и динамических

характеристик, и рассчитывается на основании приведенного выше

определения так:

 2

1

z

z

2 a dzzbS 2b )( . (8.5)

Формулой (8.5) пользуются для определения САХ сложного по форме

в плане крыла. Однако в большинстве случаев форму крыла в плане можно

привести к одной или нескольким трапециям. В этом случае САХ

рассчитывается по известным геометрическим формулам:

площадь крыла

2 lbbS к0 )(  ; (8.6)

положение САХ по размаху

к0

к0 a bb

b2b 6 1z

 

 ; (8.7)

длина САХ

)()( aк0кa z2 l

l bb1

bb  

 ; (8.8)

положение носка САХ относительно начала координат крыла

l z2xx aкa  . (8.9)

Если воспользоваться такими характеристиками крыла, как сужение и

удлинение, то формулы (8.7) и (8.8) принимают вид

)( 1 11

6 1za   ; (8.10)

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 68

docsity.com



  

 

 S

1 11

3 4b 2a )(

)( . (8.11)

Для крыла, составленного из двух трапеций, САХ и ее положение

определяются по формулам:

21

aa aa SS

zz zz 12 1 

  ; (8.12)

21

21 2

aa

aa aa zz

bb bb

  ; (8.13)

12

12 1

aa

aa aa zz

xx xy

  . (8.14)

Здесь индексом "1" обозначены параметры внутренней, а индексом "2"

- внешней секции крыла.

Для многосекционного крыла, состоящего из n трапеций, длина САХ

определяется по формуле:

 n

1i i

n

1i ia

a

S

Sb b

i

. (8.15)

Приведенные выше зависимости для определения геометрических

характеристик крыла справедливы и для других несущих поверхностей,

таких, как вертикальное и горизонтальное оперение, с той лишь разницей,

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

определяется как сумма длин консолей , т.е.

кl2l  . (8.16)

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 69

docsity.com

8.3 Геометрические характеристики аэродинамического профиля

Аэродинамический профиль является основой построения поверхности

крыла и определяет основные его характеристики. В общем случае профилем

крыла следует считать плоский замкнутый контур, полученный в результате

пересечения поверхности крыла плоскостью, перпендикулярной

строительной плоскости крыла и пересекающей переднюю и заднюю кромки

крыла.

Задачу об обтекании крыла потоком теоретическая аэродинамика

разделяет на две: задачу об обтекании прямоугольного недеформированного

крыла заданной толщины и задачу об обтекании деформированной пластины

нулевой толщины. При решении задачи обтекания поверхность крыла и

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

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

и сдвигом срединной поверхности, т.е. ординаты точек поверхности

определяются в виде:

).,(),(),(

);,(),(),(

zxyzxyzxy

zxyzxyzxy

дсн

дсв



 (8.17)

где ),( zxyв - ордината верхней поверхности крыла; ),( zxyн - ордината

нижней поверхности крыла; ),( zxyc - положительная ордината

симметричной поверхности крыла; ),( zxy д - ордината деформированной

срединной поверхности.

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

которая определяется как расстояние между крайними его точками,

являющимися точками вертикальных касательных.

В местной системе координат, начало которой лежит в носке профиля,

а ось x направлена вдоль его хорды, ординаты профиля можно представить в

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 70

docsity.com

виде:

).()()( );()()( xyxyxy xyxyxy

fсн

fсв



 (8.18)

где )(xyв - ордината верхнего контура профиля; )(xyн - ордината

нижнего контура профиля; )(xyс - положительная ордината симметричной

части профиля; )(xy f - ордината средней линии профиля.

Преобразуя (8.18), получаем:

  2

xyxyxy нвc )()()(  ; (8.19)

  2

xyxyxy нвf )()()(  . (8.20)

Для удобства сравнения профилей различных форм и размеров были

введены безразмерные, или относительные, координаты:

.

;)(

b xx

b xyy

 (8.21)

где b – хорда аэродинамического профиля.

Основными геометрическими характеристиками аэродинамического

профиля являются: максимальная относительная толщина симметричной

части профиля c и ее положение на единичной хорде cx , максимальная

кривизна f и ее положение fx .

В практике проектирования несущих поверхностей широко

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

относительную толщину и кривизну по формулам:

.

;

и

3 ff

и

3 cc

f f

yy

c c

yy

и3

и3

(8.22)

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 71

docsity.com

где индексом «3» отмечены параметры и координаты искомого профиля, а

индексом «и» - исходного профиля.

Важной характеристикой формы профиля является также

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

представляющий собой значение радиуса кривизны контура в точке 0x  . В полете под действием аэродинамических сил происходит

деформация крыла: изгиб вдоль размаха и закрутка сечений относительно

продольной оси крыла. В результате закрутки сечений происходит

увеличение местного угла атаки профиля, причем это увеличение нарастает к

концам крыла. На больших углах атаки полета самолета в концевых частях

крыла возникает срыв потока и уменьшение подъемной силы.

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

(поворот) каждого текущего сечения крыла на угол  относительно

принятой оси и отгиб носовой части профиля на угол  . Относительную величину носовой части профиля, на которую распространяется деформация

отгиба, обозначим н0x .

Излом по обводу профиля, особенно в носовой части, недопустим.

Поэтому в точке н0x необходимо обеспечить как минимум первый порядок

гладкости стыковки. В этом случае при заданном в сечении z угле отгиба  деформация будет определяться формулой:

)()()( zbxx x ztgy 20

0 z н

н

 

 . (8.23)

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

н0xx  , т.е. 0y 2 '' , то формула для определении деформации при отгибе

принимает вид:

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 72

docsity.com

)()()( zbxx x ztgy 302

0 z н

н

 

 . (8.24)

За ось крутки сечений крыла обычно принимается задняя кромка. В

этом случае деформация крутки рассчитывается по формуле:

)()(sin)( zbzx1yy 22t  , (8.25)

где y - ордината несимметричного профиля без крутки.

В ряде случаев для упрощения расчетов ввиду малости углов )(z

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

сдвига и величину деформации крутки ty рассчитывают по формуле:

)()()( zbztgx1y t  , (8.26)

Из технологических соображений удобнее бывает задавать положение

передней кромки закрученного крыла )( . zy

кпt . В этом случае формула

(8.26) преобразуется к виду:

))(( .

x1zyy кпtt  . (8.27)

Аэродинамический профиль является исходной информацией при

проектировании крыла летательного аппарата, и требование выдерживания

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

первый план по сравнению с требованиями компоновки, технологичности и

т.д. Поэтому вопросам описания обводов типа аэродинамический профиль

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

в самолетостроении.

В зависимости от назначения профиля предъявляются

соответствующие требования к его форме и геометрическим

характеристикам.

1. Дозвуковые профили характеризуются утолщенной носовой

частью, смещением максимальной толщины профиля вперед и

плавными сходами к хвостовой части (рис. 2.1, а).

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 73

docsity.com

2. Околозвуковые профили отличаются несколько более

заостренной носовой частью, смещением максимальной

толщины в более заднее положение и более плавными

формами в районе максимальной толщины (рис. 2.1,6).

3. Появившиеся в последние годы трансзвуковые, или

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

верхней линией профиля и значительным искривлением

хвостовой части (рис. 2.1, в).

4. Сверхзвуковые профили обычно представляют собой обводы с

заостренными носовой и хвостовой частями (рис. 2.1, г).

5. Гиперзвуковые профили отличаются заостренной носовой

частью и резко затупленной хвостовой частью, а также

значительным смещением максимальной толщины назад (рис.

2.1, д).

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

достаточно условна. Выбор формы профиля диктуется конкретными

задачами.

По методам описания обводов аэродинамические профили делятся на

две группы:

Рисунок 8.1. Профили.

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 74

docsity.com

1) профили, обводы которых имеют аналитическое описание;

2) профили, обводы которых заданы дискретным массивом координат.

Обводы профилей первой группы могут быть получены по заданным

аэродинамическим характеристикам в результате решения задачи обтекания

кругового цилиндра с использованием конформного отображения (профили

Жуковского, Кармана-Трефтца, Мизеса, Карафоли). При этом уравнения

контура профиля достаточно сложны. Поэтому были предложены способы

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

простого вида (степенными, строфоидами и т.п.) с последующим

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

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

30-40-х гг. В последующие годы более широкое распространение получили

профили, обводы которых получены путем численного решения

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

доводкой на основании исследований в аэродинамических трубах с целью

получения заданных характеристик.

Информация об обводах профилей второй группы обычно

представляется в виде таблицы значений координат точек, принадлежащих

контуру.

Поэтому одной из задач проектирования поверхности крыла является

задача описания обвода, заданного дискретным точечным рядом.

Прежде чем приступить к выбору функции, аппроксимирующей

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

должна отвечать эта функция.

Эти требования определяются, с одной стороны, условиями работы

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

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

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

удобства проектировщику при работе с ним, представляя собой

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 75

docsity.com

неотъемлемую часть автоматизированной системы проектирования

поверхности. Таким образом, аппроксимирующая функция должна

удовлетворять следующим основным требованиям:

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

порядка;

2) обеспечивать по возможности описание наибольшего количества типов

профилей;

3) обеспечивать гладкую аппроксимацию профиля без предварительного

графического сглаживания исходных данных;

4) обладать минимальным, но достаточным числом параметров,

варьируемых для управления формой профиля.

В настоящее время при проектировании плоских контуров типа

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

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

степенные уравнения, уравнения специальных контуров, сплайн-функции.

8.4 Проектирование поверхности линейчатого крыла

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

применение получили линейчатые поверхности, что обусловлено простотой

алгоритма их построения и высокой степенью технологичности. И если в

общем случае линейчатая поверхность не является разворачиваемой, как,

например, гиперболоид вращения, то для поверхностей крыльев можно

получить развертку с достаточно высокой степенью точности. Это важно при

изготовлении обшивки крыла, особенно на участках кессона, где ее толщина

велика и, следовательно, в качестве технологического процесса изготовления

может быть использована только гибка в случае изготовления из металлов.

Особенно широкое распространение линейчатые поверхности

получили также в силу простоты их увязки и построения графическим

способом. Однако в последние годы из-за повышения требований к

технологии подготовки производства и к точности изготовления оснастки и

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 76

docsity.com

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

автоматизации конструкторских работ все большее распространение

получают математические методы описания линейчатых поверхностей. При

этом применяются алгоритмы проектирования поверхностей на основе как

традиционного точечного задания профилей, так и аналитического описания

профилей и поверхностей.

Рассмотрим алгоритм расчета линейчатого крыла, направляющие

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

переменных:

),( zxfy IIтеор  ; ),( zxfy IIIIтеор  (8.28)

Пусть задана точка А на плановой проекции крыла с координатами Tx ,

Tz . Необходимо определить третью координату этой точки Ty . Рассмотрим

каждый этап этого алгоритма. На первом этапе по задан ной координате Tz

определяем значения координат передней и задней кромок, которые заданы

как функции одной переменной,

)(... Tкпткп zfx  ; (8.29)

)(... Tкзткз zfx  . (8.30)

Используя эти величины, можно легко определить длину текущей

хорды, координаты x точки в местной системе координат с началом на

передней кромке крыла, а также значение относительной координаты x точки

по следующим формулам:

ткпткзz хxb T ....  ; (8.31)

ткпттпр хxх ...  ; (8.32)

Tz

тпр тпр b х

х ..  . (8.33)

На следующем этапе определяем координаты точек первого и второго

теоретических сечений с равнопроцентными координатами x:

docsity.com

теортеор прIIпрIтпр хxх . . (8.34)

Умножая каждую относительную координату на величину хорды,

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

Находим координаты у точек на первом и втором теоретических сечениях:

IтпрпрI bxх . ; (8.35)

IIтпрпрII bxх . ; (8.36)

),( IпрIIпрI zxfy  ; (8.37)

),( IIпрIIIIпрII zxfy  . (8.38)

Для перехода в систему координат агрегата прибавим к прIх и прIIх

значения соответствующих координат передней кромки:

IкппрII хyх . ; (8.39)

IIкппрIIII хxх . . (8.40)

Таким образом, нам известны координаты х двух точек образующей

линейчатого крыла, которая проходит через точку А.

Для определения неизвестной координаты у точки А необходимо

подставить известные ее координаты в уравнение проекции образующей:

III

IT

III

IT

III

IT

xx xx

zz zz

yy yy

 

  

  

, (8.41)

откуда:

I III

ITIII T yzz

zzyy y 

 

))((

, (8.42)

или:

I III

ITIII T yzz

xxyy y 

 

))((

. (8.43)

При расчете сечений поверхности линейчатого крыла формулы (8.42) и

(8.43) можно использовать при взаимной перпендикулярности плоскости

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 77

docsity.com

сечения и плоскости хорд крыла. Если эти плоскости не перпендикулярны, то

неизвестные координаты точки А определяются из решения системы

линейных уравнений:

.coscoscos ;coscoscos

;coscoscos

0Pzyx 0Pzyx 0Pzyx

33T3T3T

22T2T2T

11T1T1T

 

 (8.44)

где 1P , 2P и 3P - некоторые постоянные коэффициенты.

Для решения этой системы можно использовать известный метод

Крамера. Рассмотренный алгоритм определения координат произвольной

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

заданной стреле прогиба, т.к. при аналитическом задании профилей можно

определить произвольную точку на поверхности.

Лист

Изм Лист № докум. Подп. Дата ДП 1301.02.07.00.00.00 ПЗ 78 9

docsity.com

комментарии (0)
не были сделаны комментарии
Напиши ваш первый комментарий
это только предварительный показ
консультироваться и скачать документ
Docsity не оптимизирован для браузера, который вы используете. Войдите с помощью Google Chrome, Firefox, Internet Explorer 9+ или Safari! Скачать Google Chrome