S2.1-1FLDS120, Section 2.1, June 2002 РАЗДЕЛ 2.1 Теория аэродинамики. - презентация

1 S2.1-1FLDS120, Section 2.1, June 2002 РАЗДЕЛ 2.1 Теория аэродинамики

2 S2.1-2FLDS120, Section 2.1, June 2002

3 S2.1-3FLDS120, Section 2.1, June 2002 Теории аэродинамики Нелинейно вязкая, сжимаемая среда Навье-Стокса Модель турбулентности CFD - вычислительная гидродинамика Нелинейная, невязкая, сжимаемаемая среда Эйлера Безвихревая, невязкая, сжимаемая среда Теория потенциала: панельный метод Теория тонкого тела

4 S2.1-4FLDS120, Section 2.1, June 2002 Теория тонкого тела Теория тонкого тела – это упрощенная теория потенциала. Допускается что тело вызывает только малые, вдоль оси х, возмущения в равномерном потоке. Возмущение скорости, которое было внесено в невозмущенный поток, является причиной возникновения скорости касательной к поверхности тела, называемой normal wash (перпендикулярный поток) или downwash (скос потока вниз).

5 S2.1-5FLDS120, Section 2.1, June 2002 Аэродинамическая теория в MSC.Nastran Все аэродинамические теории в MSC.Nastran основываются на теории тонкого тела. Стандартные методы Метод дипольных решеток (DLM) для дозвуковых задач Метод гармонического градиента (Зона 51) для сверхзвуковых задач

6 S2.1-6FLDS120, Section 2.1, June 2002 Метод дипольных решеток Линия диполей Точка скоса потока Только несущие поверхности тонкое тело

7 S2.1-7FLDS120, Section 2.1, June 2002 Метод Зона 51 Линия Маха Постоянное давление Точка скоса потока Зона 51

8 S2.1-8FLDS120, Section 2.1, June 2002 Дополнительные теории аэродинамики в MSC.Nastran Давление р Теория поршня Метод «Маховых» панелей Теория «ленты»

9 S2.1-9FLDS120, Section 2.1, June 2002 Приведенная частота Приведенная частота это ключевой параметр в нестационарной аэродинамике Определяется как Где циклическая частота колебаний относительная длина скорость полета приведенная скорость расстояние пройденное за один период колебания

10 S2.1-10FLDS120, Section 2.1, June 2002 Основные аэродинамические матрицы Матрица коэффициентов аэродинамического влияния может быть найдена через скос потока и коэффициент давления Основная матрица дифференциалов – через скос потока и перемещения Матрица интегралов – через аэродинамические нагрузки и коэффициент давления

11 S2.1-11FLDS120, Section 2.1, June 2002 Матрица аэродинамической жесткости Матрица аэродинамической жесткости определяется через аэродинамические нагрузки и перемещения Она может быть вычислена через три основные матрицы

12 S2.1-12FLDS120, Section 2.1, June 2002 Метод дипольных решеток: источники Метод дипольных решеток – это расширенный метод стационарных вихревых решеток для колеблющегося потока. Источники: Hedman, S. G., "Vortex-Lattice Method for Calculation of Quasi-Steady State Loadings on Thin Elastic Wings," Aeronautical Research Institute of Sweden, Report 105, October Albano, E. and Rodden, W.P. A Doublet Lattice Method for Calulating Lift Distributions on Oscillating Surfaces in Subsonic Flows. AIAA J., Vol. 7 No. 2, pp , 1969, and Vol. 7, No. 11, p. 2192, 1969 Rodden, W.P., Taylor, P.F. and McIntosh, S.C., Further Refinement of the Subsonic Doublet-Lattice Method, AIAA J. Vol. 35, No. 5, 1998 Blair, Max, A Compilation of the Mathematics Leading to the Doublet-Lattice Method, NASA Report WL-TM , 1994

13 S2.1-13FLDS120, Section 2.1, June 2002 Метод дипольных решеток (МДР) МДР может быть использован для расчета несущих поверхностей в дозвуковом потоке. Предполагается что все несущие поверхности расположены чуть-чуть не параллельно потоку, из-за наличия малых возмущений. Используется линейная аэродинамическая теория. Каждая несущая поверхность делится на маленькие трапецеидальные элементы (панели). Панели располагаются в форме ленты (одна за другой), параллельной направлению скорости свободного потока. Линии перегиба и шарнирные линии необходимо располагать на границе панелей. Опция симметрии необходима для уменьшения размерности задачи.

14 S2.1-14FLDS120, Section 2.1, June 2002 Метод дипольный решеток: теория Метод дипольных решеток получен путем дискретизации интегрального выражения. Дискретное уравнение где подынтегральная функция получена из нестационарной дозвуковой теории аэродинамики. Интеграл вычисляется с помощью высокоточной аппроксимации подынтегральной функции.

15 S2.1-15FLDS120, Section 2.1, June 2002 Метод дипольный решеток : интегрирование Исходный программный код МДР использует порабольческую апроксимацию подынтегральной функции. Она обозначается N5KA. В MSC.Nastran 70.6 была применена высокоточная аппроксимация 4-й степени. Она обозначается N5KQ. Апроксимация 4-й степени выбирается путем установки системной переменой cell 270 (QUARTICDLM) в значение 1. Пораболическая аппроксимация используется по умолчанию.

16 S2.1-16FLDS120, Section 2.1, June 2002 Метод дипольных решеток: опция N5KQ Эта опция обеспечивает высокоточную аппроксимацию разброса значений подынтегральной функции. Сравнение действительной части коэффициентов подъемной силы. Отношение геометрических размеров панели (Box Aspect Ratio ) = 10.0

17 S2.1-17FLDS120, Section 2.1, June 2002 Метод дипольных решеток Принципы моделирования Обозначение аэродинамических панелей начинается с названия макроэлемента обозначенного CAERO1 и возрастает на единицу для последующей. Обозначение следует сначала вдоль потока а потом вдоль размаха крыла Нумерация аэродинамических узлов производится независимо от нумерации структурных узлов, скалярных и особых точек, так что допускается дублирование номеров структурной и аэродинамической моделей. Аэродинамические узлы не могут иметь одинаковые обозначения Панель. Номер панели = CAERO_ID + 0 … n Угловой узел Номер узла = CAERO_ID n Макроэлемент или CAERO

18 S2.1-18FLDS120, Section 2.1, June 2002 Метод дипольных решеток Принципы моделирования Угловому узлу аэродинамической панели присваивается наименьший номер из (NCHORD+1) x (NSPAN+1). При создании второй аэродинамической сетки, порядок номеров возрастает на Рекомендованный порядок нумерации автоматически выполняется во Flightloads я аэродинамическая сетка 2 - я аэродинамическая сетка

19 S2.1-19FLDS120, Section 2.1, June 2002 Метод дипольных решеток Принципы моделирования Панели обозначаются как k-множество степеней свободы. Каждая панель имеет 2 степени свободы. Одна панель представляет собой тонкую пластину трапецевидной формы. Два ребра должны быть обязательно параллельны направлению потока. Одно параллельное ребро может быть вырождено в точку. Панели могут иметь отношение геометрических размеров меньше 3 при использовании аппроксимации N5KA и меньше 6 при N5KQ. Возможно, что в зависимости от конфигурации, большие отношения геометрических параметров могут быть использованы в N5KA и N5KQ. При использовании больших значений рекомендуется исследовать сходимость. N5KA используется по умолчанию. N5KQ может быть активирован при помощи NASTRAN команды: NASTRAN QUARTICDLM = ON

20 S2.1-20FLDS120, Section 2.1, June 2002 Метод дипольных решеток Принципы моделирования Пусть будет число панелей в длине волны минимум 12,5. Число панелей в длине волны: Однако, должно быть использовано не меньше 4-х панелей и рекомендуется всегда исследовать сходимость результатов. Во Flightloads эти расчеты будут выполняться самостоятельно. x < 0.08 v / f Панели должны быть сконцентрированы в области неоднородности среды, таких например как передняя кромка крыла, кромка оперения и линия шарнира. Используйте узкие панели на конце крыла, где нагрузка уменьшается. Доступные группы аэродинамической интерференции снижают затраты времени на генерацию матриц и изучение влияния поверхностей друг на друга. = местная хорда = базовая хорда = приведенная частота = циклическая частота v = скорость

21 S2.1-21FLDS120, Section 2.1, June 2002 Метод дипольных решеток Принципы моделирования В случае сложных почти компланарных поверхностей не допускайте что бы средние линии одних панелей находились на одной линии с ребрами других панелей НЕ ДОПУСКАЕТСЯ!

22 S2.1-22FLDS120, Section 2.1, June 2002 Зона 51 Зона 51 может быт использована для расчета несущих поверхностей в сверхзвуковом потоке Эта опция в MSC.Nastran называется Aeroelasticity II. Вычисляется матрица [A jj ] Теория метода Зона 51 аписа в ниже приведенных источниках: Chen, P. C. and D. D. Liu, "A Harmonic Gradient Method for Unsteady Supersonic Flow Calculations," Journal of Aircraft, Vol.22, No. 5, May 1985, pp Liu, D. D., D. K. James, P. C. Chen, and A. S. Pototsky, "Further Studies of Harmonic Gradient Method for Supersonic Aeroelastic Applications," Paper , DGLR/AAAF/RAeS European Forum on Aeroelasticity and Structural Dynamics, Aachen, FRG, April Предполагается что все несущие поверхности расположены чуть-чуть не параллельно потоку, из-за наличия малых возмущений. Используется линейная аэродинамическая теория. Каждая несущая поверхность делится на маленькие трапецеидальные элементы (панели).

23 S2.1-23FLDS120, Section 2.1, June 2002 Зона 51 Панели располагаются в форме ленты, параллельной направлению скорости свободного потока. Линии перегиба и шарнирные линии необходимо располагать на границе панелей. Опция симметрии необходима для уменьшения размерности задачи. Эффектом толщины пренебрегают, но можно будет приближенно посчитать коэффициенты коррекции, используя теорию сверхзвукового аэродинамического профиля второго порядка.

24 S2.1-24FLDS120, Section 2.1, June 2002 Зона 51. Теория Для сверхзвуковых скоростей Неизвестное давление принимается как постоянное для каждой панели Известные условия скоса потока заданы на расстоянии 95% от длины хорды панели в центре пролета. В MSC.Nastran аэродинамический узел задается на средней хорде в середине пролета Линия Маха Постоянное давление Скос потока, расположеный в точке находящейся на расстоянии 95% от носка панели

25 S2.1-25FLDS120, Section 2.1, June 2002 Зона 51. Теория Для n панелей В матричной форме (как и в МДР) {w} = [A]{C p } Интеграл вычисляется при помощи высокоточной аппроксимации подынтегральной функции w i a j C pj K xK x i y i x j y j,,, d j=1j=1 n = Постоянная Подинтегральная функция, получена из нестационарной сверхзвуковой теории аэродинамики

26 S2.1-26FLDS120, Section 2.1, June 2002 Метод постоянного даления Эта опция в пакете Aero 1 предназначена для вычисления матрицы A ji при сверхзвуковых расчетах. Такая же опция устанавливается для метода Зона 51 Такие же принципы принципы Не часто применялся, поэтому будте осторожны при использовании Установки: PARAM, SUPAERO, CMP or ZONA (default). Источники : Appa, K and Smith, M.J.C., Evaluation of the Constant Pressure Panel AIAA/ASME/ASCE/AHS 29 th Structures, Structural Dynamics and Materials Conference, Williamsburg Virginia, April 1988.

27 S2.1-27FLDS120, Section 2.1, June 2002 Дозвуковое крыло – теория интерференции тел Теория тонкого тела для изолированного тела имеет большую историю. Существует свыше двухсот статей по этой тематике, список литературы можно посмотреть в журнале Revell, J. D., "Second-Order Theory for Steady or Unsteady Subsonic Flow Past Slender Lifting Bodies of Finite Thickness," AIAA Journal, Vol. 6, No. 6, June 1969, pp В применяемой теории тонкого тела отображаются несущие характеристики каждого тела. Нагрузка пропорциональна коэффициенту площади поперечного сечения. Взаимное влияние тел использует метод суперпозиции сингулярностей и их подобия. Тела имеют продольную (y) и вертикальную (z) степень свободы. Первоначально интерференция крыло-тело расчитывалась для систем подобия с вихрями из МДР и диполями на цилиндрическом интерференционном теле, аписанном как узкое тело (Метод подобия Гисенга).

28 S2.1-28FLDS120, Section 2.1, June 2002 Дозвуковое крыло – теория интерференции тел Однородное уравнение для скоса потока где w w = скос потока на панели крыла расположеный на расстоянии 0.75 длины хорды от носка w l = 0=скос потока для инетерферирующих элементов вносит сингулярность одних тонких элементов на другие w s = скос потока для элемента узкого тела C w = силы, расположеные вдоль панели крыла на расстоянии 0.25 хорды r = диполи остаточного течения s = диполи тонкого тела

29 S2.1-29FLDS120, Section 2.1, June 2002 Дозвуковое крыло – теория интерференции тел Уравнение для силовых факторов где S ww = площадь панели G, S sw S sr S ss = матрицы, использующиеся для получения массовых сил (смотри раздел of Giesing, J.P., Kalman, T.P., and Rodden, W.P., "Subsonic Unsteady Aerodynamics for General Configurations; Part II, Vol. 1, - Application of the Doublet-Lattice method and the Method of Images to Lifting-Surface/Body Interference," Air Force Flight Dynamics Laboratory Report No. AFFDL-TR-71-5, Part II, Vol. 1, April 1972)

30 S2.1-30FLDS120, Section 2.1, June 2002 Дозвуковое крыло – теория интерференции тел Рисунок 3-1a. Изображены панели и элементы тонких тел. N5KA Бомбардировщик. 3 поверхности, 10 панелей

31 S2.1-31FLDS120, Section 2.1, June 2002 Дозвуковое крыло – теория интерференции тел Рисунок 3-2b. Изображены интерфереционные элементы. N5KA Бомбардировщик 3 поверхности, 10 панелей, 2 тела, 9 тонкотельных элементов и 7 интерфереционных элементов

32 S2.1-32FLDS120, Section 2.1, June 2002 Альтернатива для тел Сложность моделирования замкнутых тел может быть устранена заменой их крестообразными несущими поверхностями Весовые коэффициенты могут быть использованы для понижения эффективности несущих поверхностей Кольцеобразные крылья могут быть использованы вместо гондол со сквозным отверстием Эти технические приемы мотивируются тем фактом что удельный массовый расход через тело равен 0.7