Von Karman Wind Turbulence Model (Continuous)

Сгенерируйте непрерывную турбулентность ветра со скоростными спектрами Фон Карман

  • Библиотека:
  • Aerospace Blockset / Среда / Ветер

  • Von Karman Wind Turbulence Model (Continuous) block

Описание

Блок Von Kármán Wind Turbulence Model (Continuous) использует Фон Карман спектральное представление, чтобы добавить турбулентность в космическую модель путем передачи ограниченного полосой белого шума через соответствующие фильтры формирования. Этот блок реализует математическое представление в Военной Спецификации MIL-F-8785C и Военном Руководстве MIL-HDBK-1797. Для получения дополнительной информации см. Алгоритмы.

Ограничения

  • Замороженное полевое предположение турбулентности допустимо для случаев скорости среднего ветра.

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

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

    • Шероховатость ландшафта

    • Уровень ошибки

    • Сдвиги ветра

    • Средняя величина ветра

    • Другие метеорологические фракции (кроме высоты)

Порты

Входной параметр

развернуть все

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

Типы данных: double

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

Типы данных: double

Матрица направляющего косинуса в виде 3х3 матрицы.

Типы данных: double

Вывод

развернуть все

Скорости турбулентности, возвращенные как трехэлементный сигнал, в заданных модулях.

Типы данных: double

Турбулентность угловые уровни в виде трехэлементного вектора, в радианах в секунду.

Типы данных: double

Параметры

развернуть все

Единицы скорости ветра из-за турбулентности в виде:

МодулиСкорость ветраВысотаВоздушная скорость
Metric (MKS) Метры/секундаМетрыМетры/секунда
English (Velocity in ft/s) Ноги/секундаФутыНоги/секунда
English (Velocity in kts) УзлыФутыУзлы

Программируемое использование

Параметры блоков: units
Ввод: символьный вектор
Значения: 'Metric (MKS)' | 'English (Velocity in ft/s)' | 'English (Velocity in kts)'
Значение по умолчанию: 'Metric (MKS)'

Военная ссылка, которая влияет на приложение длин шкалы турбулентности в боковых и вертикальных направлениях в виде MIL-F-8785C, MIL-HDBK-1797, или MIL-HDBK-1797B.

Программируемое использование

Параметры блоков: spec
Ввод: символьный вектор
Значения: 'MIL-F-8785C' | 'MIL-HDBK-1797' | 'MIL-HDBK-1797B'
Значение по умолчанию: 'MIL-F-8785C'

Модель турбулентности ветра в виде:

Continuous Von Karman (+q -r)

Используйте непрерывное представление скоростных спектров Фон Карман с положительными вертикальными и отрицательными боковыми угловыми спектрами уровней.

Continuous Von Karman (+q +r)

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

Continuous Von Karman (-q +r)

Используйте непрерывное представление скоростных спектров Фон Карман с отрицательными вертикальными и положительными боковыми угловыми спектрами уровней.

Continuous Dryden (+q -r)

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

Continuous Dryden (+q +r)

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

Continuous Dryden (-q +r)

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

Discrete Dryden (+q -r)

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

Discrete Dryden (+q +r)

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

Discrete Dryden (-q +r)

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

Выборы Continuous Von Kármán соответствуют описаниям передаточной функции.

Программируемое использование

Параметры блоков: model
Ввод: символьный вектор
Значения: 'Continuous Von Karman (+q +r)' | 'Continuous Von Karman (-q +r)' | 'Continuous Dryden (+q -r)' | 'Continuous Dryden (+q +r)' | 'Continuous Dryden (-q +r)' | 'Discrete Dryden (+q -r)' | 'Discrete Dryden (+q +r)' | 'Discrete Dryden (-q +r)'
Значение по умолчанию: 'Continuous Von Karman (+q +r)'

Измеренная скорость ветра на высоте 20 футов (6 метров) в виде действительного скаляра, который обеспечивает интенсивность для низковысотной модели турбулентности.

Программируемое использование

Параметры блоков: W20
Ввод: символьный вектор
Значения: действительный скаляр
Значение по умолчанию: '15'

Измеренное направление ветра на высоте 20 футов (6 метров) в виде действительного скаляра, который является углом, чтобы помочь в преобразовании низковысотной модели турбулентности в тело, координирует.

Программируемое использование

Параметры блоков: Wdeg
Ввод: символьный вектор
Значения: действительный скаляр
Значение по умолчанию: '0'

Вероятность интенсивности турбулентности, превышаемой в виде 10^-2 - Light, 10^-1, 2x10^-1, 10^-3 - Moderate, 10^-4, 10^-5 - Severe, или 10^-6. Выше 2 000 футов интенсивность турбулентности определяется из интерполяционной таблицы, которая дает интенсивность турбулентности как функцию высоты и вероятность превышаемой интенсивности турбулентности.

Программируемое использование

Параметры блоков: TurbProb
Ввод: символьный вектор
Значения: '2x10^-1'| '10^-1' | '10^-2 - Light' | '10^-3 - Moderate'| '10^-4' | '10^-5 - Severe'| '10^-6'
Значение по умолчанию: '10^-2 - Light'

Длина шкалы турбулентности выше 2 000 футов в виде действительного скаляра. Эта длина принята постоянная.

Из военных технических требований 1 750 футов рекомендуются для продольной длины шкалы турбулентности спектров Драйдена.

Примечание

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

Программируемое использование

Параметры блоков: L_high
Ввод: символьный вектор
Значения: действительный скаляр
Значение по умолчанию: '762'

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

Программируемое использование

Параметры блоков: Wingspan
Ввод: символьный вектор
Значения: действительный скаляр
Значение по умолчанию: '10'

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

Программируемое использование

Параметры блоков: ts
Ввод: символьный вектор
Значения: действительный скаляр
Значение по умолчанию: '0.1'

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

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

Программируемое использование

Параметры блоков: Seed
Ввод: символьный вектор
Значения: четырехэлементный вектор
Значение по умолчанию: '[23341 23342 23343 23344]'

Чтобы сгенерировать сигналы турбулентности, установите этот флажок.

Программируемое использование

Параметры блоков: T_on
Ввод: символьный вектор
Значения: 'on' | 'off'
Значение по умолчанию: 'on'

Алгоритмы

развернуть все

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

 MIL-F-8785CMIL-HDBK-1797
Продольный

Φu(ω)

2σu2LuπV1[1+(1.339LuωV)2]56

2σu2LuπV1[1+(1.339LuωV)2]56

Φp(ω)

σw2VLw0.8(πLw4b)131+(4bωπV)2

σw22VLw0.8(2πLw4b)131+(4bωπV)2

Ответвление

Φv(ω)

σv2LvπV1+83(1.339LvωV)2[1+(1.339LvωV)2]116

2σv2LvπV1+83(2.678LvωV)2[1+(2.678LvωV)2]116

Φr(ω)

(ωV)21+(3bωπV)2Φv(ω)

(ωV)21+(3bωπV)2Φv(ω)

Вертикальный

Φw(ω)

σw2LwπV1+83(1.339LwωV)2[1+(1.339LwωV)2]116

2σw2LwπV1+83(2.678LwωV)2[1+(2.678LwωV)2]116

Φq(ω)

±(ωV)21+(4bωπV)2Φw(ω)

±(ωV)21+(4bωπV)2Φw(ω)

Переменная b представляет размах крыла самолета. Переменные Lu, Lv, Lw представляют длины шкалы турбулентности. Переменные σu, σv, σw представляют интенсивность турбулентности:

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

pg=wgy

qg=wgx

rg=vgx

pg=wgy

qg=wgx

rg=vgx

pg=wgy

qg=wgx

rg=vgx

Изменения влияют только на вертикаль (qg) и ответвление (rg) турбулентность угловые уровни.

Следует иметь в виду, что продольная турбулентность угловой спектр уровня, Фp (ω), является рациональной функцией. Рациональная функция выведена из подбора кривых комплексная алгебраическая функция, не вертикальный скоростной спектр турбулентности, Фw (ω), умноженный на масштабный коэффициент. Поскольку турбулентность, которую угловые спектры уровня вносят меньше в ответ порыва самолета, чем скоростные спектры турбулентности, он может объяснить изменения их определений.

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

ВертикальныйОтветвление

Фq (ω)

Фq (ω)

−Фq (ω)

−Фr (ω)

Фr (ω)

Фr (ω)

Чтобы сгенерировать сигнал с правильными характеристиками, модульным отклонением, ограниченный полосой белый шумовой сигнал передается посредством формирования фильтров. Формирующиеся фильтры являются приближениями скоростных спектров Фон Карман, которые допустимы в области значений нормированных частот меньше чем 50 радианов. Эти фильтры могут быть найдены и в Военном Руководстве MIL-HDBK-1797 и в ссылке Ли и Чаном.

Следующие две таблицы показывают передаточные функции.

 MIL-F-8785C
Продольный

Hu(s)

σu2πLuV(1+0.25LuVs)1+1.357LuVs+0.1987(LuV)2s2

Hp(s)

σw0.8V(π4b)16Lw13(1+(4bπV)s)

Ответвление

Hv(s)

σv1πLvV(1+2.7478LvVs+0.3398(LvV)2s2)1+2.9958LvVs+1.9754(LvV)2s2+0.1539(LvV)3s3

Hr(s)

sV(1+(3bπV)s)Hv(s)

Вертикальный

Hw(s)

σw1πLwV(1+2.7478LwVs+0.3398(LwV)2s2)1+2.9958LwVs+1.9754(LwV)2s2+0.1539(LwV)3s3

Hq(s)

±sV(1+(4bπV)s)Hw(s)

 MIL-HDBK-1797
Продольный

Hu(s)

σu2πLuV(1+0.25LuVs)1+1.357LuVs+0.1987(LuV)2s2

Hp(s)

σw0.8V(π4b)16(2Lw)13(1+(4bπV)s)

Ответвление

Hv(s)

σv1π2LvV(1+2.74782LvVs+0.3398(2LvV)2s2)1+2.99582LvVs+1.9754(2LvV)2s2+0.1539(2LvV)3s3

Hr(s)

sV(1+(3bπV)s)Hv(s)

Вертикальный

Hw(s)

σw1π2LwV(1+2.74782LwVs+0.3398(2LwV)2s2)1+2.99582LwVs+1.9754(2LwV)2s2+0.1539(2LwV)3s3

Hq(s)

±sV(1+(4bπV)s)Hw(s)

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

Примечание

Тот же результат передаточных функций после оценки длин шкалы турбулентности. Различия в длинах шкалы турбулентности и балансе передаточных функций турбулентности возмещены.

Ссылки

[1] Американское военное руководство MIL-HDBK-1797, 19 декабря 1997.

[2] Американская военная спецификация MIL-F-8785C, 5 ноября 1980.

[3] Мел, Чарльз, Т.П. Нил, Т.М. Харрис, Фрэнсис Э. Притчар и Роберт Дж. Вудкок. "Справочная информация и Руководство пользователя для MIL-F-8785B (ASG), 'Военные Управляющие Спецификацией Качества Пилотируемых Самолетов", AD869856. Buffalo, Нью-Йорк: Авиационная лаборатория Корнелла, август 1969.

[4] Hoblit, Фредерик М., нагрузки порыва на самолет: Концепции и приложения. Образовательный ряд AIAA, 1988.

[5] Ly, U. и И. Чан. "Расчет временного интервала ковариационных матриц порыва самолета". Бумага AIAA 80-1615. Представленный на атмосферной конференции бортмехаников, Дэнверзе, MA, 11-13 августа 1980.

[6] Макруер, Дуэн, Ирвинг Ашкенас и Данстан Грэм. Динамика самолета и автоматическое управление. Принстон: Издательство Принстонского университета, июль 1990.

[7] Moorhouse, Дэвид Дж. и Роберт Дж. Вальдшнеп. "Справочная информация и руководство пользователя для MIL-F-8785C, 'Военные управляющие спецификацией качества пилотируемых самолетов". ADA119421. Рейс динамическая лаборатория, июль 1982.

[8] Макфарлэнд, R. "Стандартная кинематическая модель для симуляции рейса в NASA-Ames". НАСА CR-2497. Computer Sciences Corporation, январь 1975.

[9] Tatom, Франк Б., Стивен Р. Смит и Джордж Х. Фичтл. "Симуляция атмосферных турбулентных порывов и градиентов порыва", бумага AIAA 81-0300. Космическая научная встреча, Сент-Луис, MO., 12-15 января 1981.

[10] Yeager, Джесси. "Реализация и тестирование моделей турбулентности для симуляции F18-HARV". НАСА CR-1998-206937. Хэмптон, ВА: Lockheed Martin Engineering & Sciences, март 1998.

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью Simulink® Coder™.

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