Это издание является перепечаткой страниц 1, 3, 6-13, 54-58 из отчета
Боткин Н.Д., Пацко В.С., Турова В.Л. Разработка численных методов построения экстремальных ветровых возмущений, действующих на самолет на этапе посадки. Разработка алгоритмов построения экстремальных ветровых возмущений. Отчет о научно-исследовательской работе (промежуточный). Свердловск, 1987.Эти страницы посвящены моделированию микровзрыва. Хотя мое имя и не вошло в число авторов в работе по моделированию я принимал непосредственное участие.N ГР01880003467, Инв. N 02880054701. 58 с.
Я буду благодарен всем, кто сообщит о замеченных опечатках (которые могли появиться при распознавании оригинала) и неточностях (адрес электронной почты: iagsoft@imm.uran.ru).
Последние версии этого текста, приложение (текст программы) в электронном виде и другие сопутствующие документы размещаются по адресу: http://iagsoft.imm.uran.ru/ElePub/wind.html.
Составитель выражает благодарность Г.Ф.Корниловой и Н.Н.Моргуновой за помощь в "борьбе" с LATEX'ом. А также Институту математики и механики УрО РАН, на чьей технике производились работы.
Отчет 58 стр., 31 рисунок, 12 источников
СДВИГ ВЕТРА, МИКРОВЗРЫВ, ЭКСТРЕМАЛЬНОЕ ВЕТРОВОЕ ВОЗМУЩЕНИЕ, ЗАДАЧА ПОСАДКИ,
ДИФФЕРЕНЦИАЛЬНАЯ ИГРА, ЧИСЛЕННЫЕ МЕТОДЫ
Исследуется задача посадки самолета в условиях ветрового возмущения, обусловленного микровзрывом.
Цель работы - создание на ЭВМ программы микровзрыва, стыковка ее с комплексом программ моделирования процесса посадки.
Разработана программа микровзрыва, проведено моделирование процесса посадки для различных вариантов расположения центра микровзрыва относительно глиссады снижения и при нескольких вариантах параметров микровзрыва. Использовались два способа управления: традиционный и минимаксный (игровой), основанный на методах теории дифференциальных игр.
Полученные результаты могут быть полезны при совершенствовании комплексных самолетных тренажеров, а также при испытании новых методов автоматической посадки.
Отчет является непосредственным продолжением работы /I/. Как и в /I/,
рассматривается задача посадки самолета. В /I/ основное внимание
было уделено вопросу имитации экстремальных ветровых возмущений при помощи
методов теории дифференциальных игр /2/. В этом отчете изучается действие
сдвига ветра, вызванного локализованным в пространстве явлением
микровзрыва
Микровзрыв образуется за счет нисходящего потока воздуха, который
ударяется о поверхность земли и растекается затем с образованием
вихря. При прохождении самолетом на взлете или посадке зоны микровзрыва
происходит резкое изменение скорости ветра. Традиционные способы управления
не справляются с таким возмущением. В связи с этим в настоящее время
делаются попытки построить новые способы управления
Используемая в отчете модель микровзрыва взята из статьи /3/. Микровзрыв задается тремя параметрами: скорость ветра в центральной части, высота центральной части, радиус вихря. Для конкретизации расположения микровзрыва относительно глиссады снижения следует выбрать координаты центра микровзрыва по осям xg, zg земной нормальной системы координат. Ось xg направлена по оси взлетно-посадочной полосы (ВПП), ось zg ортогональна к ней в плоскости земли, ось yg лежит в вертикальной плоскости.
Расчетные формулы и программа микровзрыва приведены в разделе I. Текст программы помещен в приложении. Программа разработана студентами Уральского госуниверситета им. А.М.Горького А.Г.Ивановым и Е.А.Рязанцевой.
Во втором разделе выписаны соотношения, определяющие действующий (принятый в настоящее время) способ управления самолетом на посадке. Действующий способ, наряду с минимаксным (игровым), использовался при моделировании процесса посадки. Игровой способ описан в /I/.
Раздел 3 посвящен результатам моделирования процесса посадки при наличии микровзрыва. В проведении расчетов участвовала техник отдела динамических систем ИММ Э.Е.Тарасьева.
Работа выполнена в соответствии с договором
I.I РАСЧЕТНЫЕ ФОРМУЛЫ
Математическая модель микровзрыва взята из статьи /З/. Стационарная стадия микровзрыва идеализируется в виде трехмерного осесимметричного вихревого поля, в котором выделяется тороидальная область ("ядро"), где скорость, начиная от нулевой в центре, линейно увеличивается по радиусу до границы ядра. Вне ядра вихревое поле задается функцией потока через приближенно вычисляемые комбинации полных эллиптических интегралов. Численное дифференцирование функции потока дает радиальную и вертикальную составляющие скорости ветра. Первая из них затем раскладывается на две компоненты: параллельно и перпендикулярно оси ВПП.
Параметры микровзрыва (рис.I.I):
V0 - скорость в центральной части микровзрыва, [м/с],
H0 - высота центральной части, [м],
R0 - радиус вихря, [м].
Радиус RC ядра вихря полагается равным 0.8 H0.
Геометрическое положение микровзрыва определяется координатами его центра в горизонтальной плоскости. При расчете скорости ветра в заданной точке e используются величины (рис.1.2):
Y - высота точки e,
R - расстояние от точки e до центральной оси микровзрыва,
r1 - минимальное расстояние от точки e до основного кольца вихря,
r2 - минимальное расстояние от точки e до зеркального кольца вихря,
r3 - максимальное расстояние от точки e до основного кольца вихря,
r4 - максимальное расстояние от точки e до зеркального кольца вихря
Кроме того:
k1 - модуль эллиптических интегралов основного вихря,
k2 - модуль эллиптических интегралов зеркального вихря,
ap - аппроксимация комбинации эллиптических интегралов основного вихря,
am - аппроксимация комбинации эллиптических интегралов зеркального вихря,
STR - функция потока,
CIR - циркуляционная сила вихря,
WR - радиальная скорость ветра в точке e,
WY - вертикальная скорость ветра в точке e,
WX, WZ - разложение радиальной скорости параллельно и перпендикулярно оси ВПП.
Приведем последовательность расчетных формул. Циркуляционная сила
CIR
связана с
V0,
H0,
R0 соотношением
Расстояния ri:
Далее:
Значение функции потока в точке e:
Радиальная и вертикальная составляющие скорости ветра вне ядра находятся
численным дифференцированием (b = O.2):
Разложение радиальной составляющей:
Внутри ядра скорость ветра находится пересчетом с границы пропорционально
расстоянию от центра ядра.
1.2. Описание программы
Модель микровзрыва реализуется подпрограммами STR, WAP, VV на языке Фортран (текст в приложении).
I) Подпрограмма-функция
STR.
Назначение: вычисление функции потока
STR.
Входные параметры:
H0,
R0,
Y,
R,
CIR.
Результат:
STR.
Используемые подпрограммы: нет.
Форма обращения:
STR(
H0, R0, Y, R).
Примечание. Циркуляционная сила CIR подсчитывается в основной программе, реализующей нелинейную модель движения самолета, и передается в подпрограмму STR через COMMON-блок.
2) Подпрограмма
WAP.
Назначение: вычисление составляющих
WR,
WY
на основе функции
STP.
Входные параметры:
V0,
H0,
R0,
Y,
R.
Используемые подпрограммы:
STR.
Форма обращения:
CALL
WAP(
V0,
H0,
R0, Y, R,
WR,
WY).
3) Подпрограмма
VV.
Назначение: вычисление составляющих
WX,
WY,
WZ.
Входные параметры:
V0,
H0,
R0,
X,
Y,
Z.
Выходные параметры:
WX,
WY,
WZ.
Используемые подпрограммы:
WAP
Форма обращения:
CALL
VV(
V0,
H0,
R0, X, Y, Z,
WX,
WY,
WZ).
Параметры микровзрыва V0, H0, R0, координаты его центра по осям xq, zq задаются в основной программе движения самолета. Там же на каждом шаге интегрирования просчитываются величины X, Y, Z.
В отчете приведены результаты моделирования процесса посадки в условиях сдвига ветра, возникающего при прохождении самолета зоны микровзрыва. Применяемая модель микровзрыва весьма проста и может быть использована в летных тренажерах, а также в устройствах имитации ветрового возмущения при испытании новой аппаратуры.
Текст программы, приведенный здесь не соответствует тексту отчета (в нем были ошибки), а взят из более поздних разработок.
В отличие от отчета циркуляционная сила CIR вычисляется не в основной программе, а в подпрограмме VV. Т. е. приводимый фрагмент полностью автономен и не требует дополнительной поддержки в основной программе.
В электронном виде Вы можете скачать этот фрагмент по адресу: http://home.imm.uran.ru/iagsoft/ElePub/windprog.for.
C ..................................................... C : МОДЕЛЬ МИКРОВЗРЫВА : C : Xg,Yg,Zg - КООРДИНАТЫ ИНТЕРЕСУЮЩЕЙ НАС ТОЧКИ : C : X=Xg-DX, Y=Yg-DY, Z=Zg-DZ - КООРДИНАТЫ ОТНОСИ- : C : ТЕЛЬНО ЦЕНТРА : C : МИКРОВЗРЫВА : C : : C : WX,WY,WZ - СКОРОСТИ В ЭТОЙ ТОЧКЕ ; : C : V0,H0,R0 - ХАРАКТЕРИСТИКИ МИКРОВЗРЫВА. : C :...................................................: SUBROUTINE VV (X,Y,Z,WX,WY,WZ,V0,H0,R0) C ................................................... C П/ПРОГРАММА ОПРЕДЕЛЕНИЯ СКОРОСТИ ВЕТРА : C ПРОИЗВОЛЬНОЙ ТОЧКЕ : C ..................................................: COMMON/CIR/CIP CIP=2*V0*R0/(1.-1./(1. +(2*H0/R0)**2)**1.5) C ................................................... C RC - РАДИУС ЯДРА МИКРОЗРЫВА : C RD - РАССТОЯНИЕ ОТ ЦЕНТРА БЛИЖАЙШЕГО ЯДРА : C ДО ИССЛЕДУЕМОЙ ТОЧКИ : C ..................................................: RC=0.8*H0 R=SQRT(X**2+Z**2) RD=SQRT((ABS(R)-R0)**2+(ABS(Y)-H0)**2) C ................................................... C ПРОВЕРКА: НАХОДИТСЯ ЛИ ИССЛЕДУЕМАЯ ТОЧКА В ЯДРЕ : C ..................................................: IF(RD-RC) 21,22,22 C ................................................... C РЕАЛИЗАЦИЯ ВАРИАНТА - ТОЧКА В ЯДРЕ : C ..................................................: 21 Q=RD/RC IF(Q)24,23,24 C ................................................... C (R1, H1) - ВЕКТОР: : C ЦЕНТР ЯДРА - ИССЛЕДУЕМАЯ ТОЧКА : C ..................................................: 24 R1=R-SIGN(R0,R) Y1=Y-SIGN(H0,Y) CALL WAP(V0,H0,R0,SIGN(H0,Y)+Y1/Q,SIGN(R0,R)+R1/Q,WR1,WY1) 23 WR=Q*WR1 WY=Q*WY1 GO TO 25 C ................................................... C РЕАЛИЗАЦИЯ ВАРИАНТА - : C ИССЛЕДУЕМАЯ ТОЧКА ВНЕ ЯДРА : C ..................................................: 22 IF(ABS(R).LE.0.5) R=0.5 CALL WAP(V0,H0,R0,Y,R,WR,WY) 25 WX=WR*X/R WZ=WR*Z/R WY=-WY RETURN END FUNCTION STR(H0,R0,Y,R) C ................................................. C ВЫЧИСЛЕНИЕ ВИХРЕВОЙ ФУНКЦИИ STR : C ................................................: COMMON/CIR/CIP C ................................................. C ВЫЧИСЛЕНИЕ РАССТОЯНИЙ R : C ................................................: R1=SQRT((Y-H0)**2+(R-R0)**2) R2=SQRT((Y+H0)**2+(R-R0)**2) R3=SQRT((Y-H0)**2+(R+R0)**2) R4=SQRT((Y+H0)**2+(R+R0)**2) C ................................................ C ВЫЧИСЛЕНИЕ КОЭФФИЦИЕНТОВ RK : C ...............................................: RK1=((R3-R1)/(R3+R1))**2 RK2=((R4-R2)/(R4+R2))**2 C ................................................ C ВЫЧИСЛЕНИЕ ПРИБЛИЖЕНИЙ ЭЛЛИПТИЧЕСКИХ : C ИНТЕГРАЛОВ : C ...............................................: AR=0.788*RK1/(0.25+0.75*SQRT(1.-RK1)) AM=0.788*RK2/(0.25+0.75*SQRT(1.-RK2)) C ................................................. C ВЫЧИСЛЕНИЕ ВИХРЕВОЙ ФУНКЦИИ : C ................................................: STR=-CIP*((R1+R3)*AR-(R2+R4)*AM)/6.2831852 RETURN END SUBROUTINE WAP(V0,H0,R0,Y,R,WR,WY) C ................................................ C ОПРЕДЕЛЕНИЕ НАПРАВЛЕНИЯ ВЕТРА ПО : C ВИХРЕВОЙ ФУНКЦИИ ВНЕ ЯДРА : C D - ДЕЛЬТА, ПРИМЕНЯЕМАЯ ПРИ ВЫЧИСЛЕНИИ : C ПРОИЗВОДНОЙ : C ...............................................: D=0.2 C ................................................ C ВЫЧИСЛЕНИЕ РАДИАЛЬНОЙ И ВЕРТИКАЛЬНОЙ : C СОСТАВЛЯЮЩИХ СКОРОСТИ ВЕТРА : C ...............................................: STR1=STR(H0,R0,Y,R) WY=(STR1-STR(H0,R0,Y,R+D))/(R*D) WR=(STR(H0,R0,Y-D,R)-STR1)/(R*D) RETURN END
© сканирование, верстка, IAGSoft (iagsoft@imm.uran.ru), 2000
Версия этой страницы для распечатки в формате PostScript
Программа построения графического распределения потоков ветра в сечении микровзрыва
Домашняя страница | Электронные публикации | Печатные работы |