WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


На правах рукописи

Кондратьева Людмила Александровна

ОБРАТНЫЕ ЗАДАЧИ ДИНАМИКИ ПОДВОДНЫХ И ЛЕТАТЕЛЬНЫХ АППАРАТОВ

Специальность 01.02.01 – Теоретическая механика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Москва – 2012

Работа выполнена на кафедре «Математическая кибернетика» в Московском авиационном институте (национальном исследовательском университете).

Научный руководитель доктор физико-математических наук, профессор Галиуллин Ильяс Абдэльхакович

Официальные оппоненты: доктор физико-математических наук, профессор Гребеников Евгений Александрович доктор физико-математических наук, профессор Мухарлямов Роберт Гарабшевич

Ведущая организация: ФГБОУ ВПО Московский государственный технологический университет «СТАНКИН»

Защита состоится « 27 » апреля 2012 г. в 12 ч. 00 мин. на заседании Диссертационного совета Д 212.125.14 Московского авиационного института (национального исследовательского университета) по адресу: 125993, Москва, А-80,ГСП-3, Волоколамское шоссе, 4, Ученый совет МАИ.

С диссертацией можно ознакомиться в библиотеке Московского авиационного института (национального исследовательского университета).

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

Автореферат разослан «__»_________2012 г.

Ученый секретарь Диссертационного совета Д 212.125.кандидат физико-математических наук, доцент ______________ Гидаспов В.Ю.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ



Актуальность работы.

Обратные задачи занимают важное место в исследовании теоретикомеханических моделей. Математически строгая формулировка понятия обратных задач динамики была дана А.С. Галиуллиным. Соответствующая тематика получила интенсивное развитие в работах его последователей, прежде всего, И.А.

Мухаметзянова и Р.Г Мухарлямова, причем, начиная с работы И.А. Галиуллина, стало возможным исследовать подобные задачи не только в евклидовых пространствах, но и на произвольных дифференцируемых многообразиях.

Многие обратные задачи динамики связаны с условиями программного движения аэрогидродинамических или космических аппаратов, т.е. с выбором функций управления или параметров аппарата, обеспечивающих его движение по траектории с заданными свойствами. Типичным примером такого рода является решённая В.Т. Грумондзом задача о движении по винтовой линии центра динамически симметричного подводного аппарата.

Важной проблемой является устойчивость соответствующего программного движения относительно параметров процесса. Разные классы подобного рода обратных задач механики рассматривались в работах А.С. Галиуллина, Р.Г. Мухарлямова, О.М. Алифанова, Е.А. Гребеникова и Ю.А. Митропольского. В частности, задача о движении геостационарного спутника решена Е.А. Гребениковым, Ю.А. Митропольским и Ю.А. Рябовым. Общим вопросам динамики космических аппаратов посвящены работы В.В. Белецкого. Устойчивость движения спутников изучалась А.П. Маркеевым и О.В. Холостовой. В монографии Ю.А. Митропольского, О.Б. Лыковой определяется эволюция свободных (при отсутствии всех возмущений кроме влияния силы тяжести) орбит спутников и исследуется устойчивость этих орбит. Перспективное направление, связанное с малыми космическими аппаратами как эволюционной ступенью перехода к микро и наноспутникам в последнее время развивается О.М. Алифановым.

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

Большое значение имеет задача аналитического приближения программного движения и оценки его погрешности, решаемая с помощью тех или иных методов анализа обыкновенных дифференциальных уравнений. Одной из наиболее эффективных процедур для решения подобных задач является версия метода гармонического баланса, описанная в работах Б. Деламотта и Д. Поланда. Весьма важной является задача выбора управления, обеспечивающего устойчивое периодическое движение механического объекта в том или ином фазовом пространстве, что в математической формулировке означает существование устойчивого предельного цикла для соответствующей системы дифференциальных уравнений. Методы решения данной задачи на плоскости хорошо известны, однако уже в размерности 3 возникают существенные трудности. С этой точки зрения весьма полезен восходящий к Н.М. Крылову, Н.Н. Боголюбову метод интегральных многообразий и, как его обобщение, развитый Р. Смитом, М. Миклавчичем, А.В.

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

Ряд обратных задач аэродинамики сводится к интегральным или интегродифференциальным уравнениям. Эффективным инструментом решения таких уравнений, служит, как отмечено в работах А.И. Задорина, метод сплайнколлокации.

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

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

1) применительно к движению подводных аппаратов по заданной траектории;

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

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

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

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

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

Научная новизна. В диссертации представлены следующие основные результаты, имеющие научное и прикладное значение.

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

2. На основе понятия инерциального многообразия решена обратная задача выбора функций управления спутником (представляемым материальной точкой), обеспечивающих его полёт вокруг небесного объекта в заданной плоскости с полярными координатами (r,) по замкнутой устойчивой (в обобщённых & & координатах (r,r,) ) траектории. Показано, что в реальных пространственных координатах соответствующее движение будет периодическим или условнопериодическим в зависимости от соотношений между параметрами системы.

3. На основе одной из версий метода гармонического баланса найдены аналитические приближения для замкнутых устойчивых орбит космического & & аппарата в фазовых координатах (r,r,), а также оценки периода обращения.

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

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

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

• при проектировании подводных аппаратов;

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

• при создании систем слежения и выбора устойчивых траекторий облета космического объекта, например, управление спутником-инспектором для диагностики и устранения причин выхода из строя других космических аппаратов;

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

• при изучении вопросов допустимой аппроксимации программного движения, а также вопросов его устойчивости.

Апробация работы. Основные результаты работы докладывались и обсуждались на следующих научных конференциях и симпозиумах:

• Sixth International Symposium on Classical and Celestial Mechanics, Moscow – Velikie Luki, 2007;

• XI, XII, XIV-XVI Международный научно-технический семинар «Современные технологии в задачах управления, автоматики и обработки информации», Алушта: 2002, 2003, 2005-2007.

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

Публикации Основные результаты работы опубликованы в статьях [1-4] в журналах, входящих в Перечень ВАК, в других изданиях [5, 6], а также в материалах научных конференций [7-12].

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, заключения, списка литературы (78 источников), 12 рисунков и 1 таблицы. Объем диссертации - 108 м.п.с.

СОДЕРЖАНИЕ РАБОТЫ

Во введении даётся обзор различных направлений развития обратных задач динамики. Данная теория имеет свою историю, восходящую к работам Ньютона, Гельмгольца и Биркгофа, российских ученых Суслова, Горячева и Жуковского, французских ученых Пуанкаре и Картана.

В формулировке А.С. Галиуллина обратные задачи динамики - это задачи по определению приложенных к механической системе активных сил и моментов, а также параметров данной системы и наложенных на неё связей, при которых становится возможным одно из движений с заданными свойствами.

Классическая постановка обратной задачи связана с определением параметров механической системы, обеспечивающих движение с заданными свойствами. На этой основе в диссертации рассматривается винтовое движение подводного аппарата.

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





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

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

В разделе 1.1 описывается модель движения твердого тела в жидкости1 при отсутствии сил вязкости. Рассматриваются следующие векторные уравнения движения аппарата в предположении, что его моделью служит твёрдое тело, а формой и движением Земли можно пренебречь:

d Q + Q = G + A + T, d t (1.1) d LO G A T + LO +VO Q = M + M + M.

d t Q, LO - соответственно количество движения и кинетический момент тела и Здесь:

идеальной жидкости. В правой части первого уравнения учитываются внешние силы, а именно (см. рис.1): сила тяжести с координатами точки приложения G xc, yc, zc ; гидродинамическая архимедова сила A; сила тяги T. Векторные величины в правой части второго уравнения представляют собой моменты сил относительно центра величины O.

Рис.Проецируя векторные уравнения (1.1) на связанные оси координат, имеем шесть соотношений, содержащих следующие величины:

u,v,w VO компоненты вектора скорости точки О;

p, q, r компоненты угловой скорости тела ;

Ki j -безразмерные коэффициенты, которые вычисляются через симметричную 6 ( ) -матрицу присоединенных масс { }, характеризующих инерционность i j и форму аппарата;

cx,cy,cz,mx,my,mz - аэрогидродинамические характеристики.

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

Грумондз В.Т., Яковлев Г.А. Алгоритмы аэрогидробаллистического проектирования. М.: Изд-во МАИ, 1994.

коэффициентов вращательных производных компонент q c,c,cr,cz,c,cz аэрогидродинамической силы ;

y y y z коэффициентов вращательных производных аэрогидродинамических моментов r m,m,mz,my,m,mq ;

z z y y углов отклонения рулей Г, В, Э, соответственно, горизонтальных, вертикальных и элеронов.

Кроме того, в проекции уравнений (1.1) входят углы ,, , определяющие положение связанных осей относительно земных (см. рис.1), так что & & p = + sin, & q = cos cos +& sin, & r = - cos sin +& cos.

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

При этом предполагается, что:

1) угол атаки и угол скольжения, определяющие вектор скорости в связанных осях, равны нулю;

2) шаг винтовой линии постоянный, т.е. = 0 = const ;

=3) угол крена аппарата постоянный, т.е. ;

& = e = const 4) угловая скорость прецессии аппарата постоянна, т.е. ;

5) углы отклонения рулей Г, В, Э постоянны;

y z 6) моменты инерции вдоль осей и полагаются равными и в силу динамической симметрии аппарата имеют место соотношения m = my, m = m, mr = mq, z z y z y q c = -cz, c = -c, cr = -cz.

y y z y В этой постановке обратной задачи задана траектория движения аппарата, требуется найти соотношения, определяющие такое движение.

В разделе 1.4 решена задача определения стационарного движения аппарата со - смещенным центром масс при допущениях 1) 6). В принятых допущениях 2 p0, q0, rкомпоненты угловой скорости и постоянны.

0 = = p0 + q0 + rВ этом случае проецирование уравнений (1.1) на связанные оси даёт систему p0, q0, r0, VOВ, , , алгебраических уравнений относительно неизвестных.

Г Э Пусть выполняются следующие условия:

tg 0 = yc / zc 1) ;

2) cx = mx b, что достигается подбором угла Э, от которого зависят величины сx и mx.

Тогда можно найти величину VO и значение sin0 (0 < 0) через компоненты p0, q0, r0 при некотором фиксированном значении Э. С учетом угловой скорости динамической симметрии аппарата находится, также, постоянная угловая скорость & 0. По условию = r =, следовательно, центр аппарата будет обращаться вокруг вертикальной оси цилиндра, на котором расположена винтовая линия, с угловой & скоростью = e = 0. При этом радиус и шаг пространственной спирали определяются соотношениями:

l l R =VO cos0 h = VO sin ,, n sin 0 n sin где постоянная величина n зависит от параметров механической системы.

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

Во второй главе объектом исследования служит летательный аппарат, совершающий облет космического тела по плоской орбите и моделируемый материальной точкой. Как известно, интеграл кинетического момента определяет особенность траекторий материальной точки в каком-либо центральном (в том числе гравитационном) поле, состоящую в том, что они необходимо являются плоскими кривыми. Предполагается, что сам объект, вокруг которого совершает облёт летательный аппарат, имеет пренебрежимо малую массу. Это может быть другой летательный аппарат или небесное тело небольших размеров, например, недавно открытые спутники Урана. Заметим, что своего рода «центр орбиты» может вообще не представлять собой материальный объект, а быть фиксированной геометрической точкой в пространстве, некоторым навигационным центром.

В данной работе рассматривается класс систем трёх обыкновенных дифференциальных уравнений (ОДУ) первого порядка, для которых удаётся установить существование устойчивого предельного цикла. Показано, что при подходящем выборе управляющих сил уравнения такого типа описывают плоское движение летательного аппарата и, таким образом, позволяют установить наличие & & устойчивой периодической орбиты в фазовых переменных (r,r,), где (r,) - исходная полярная система координат в фиксированной плоскости обращения спутника.

Существование устойчивого предельного цикла для рассмотренного класса систем дифференциальных уравнений доказывается с помощью адаптированной к конечномерному случаю техники, основанной на понятии инерциального многообразия. В общем случае инерциальное многообразие (ИМ)2 гладкого Miklavcic M. A sharp condition for existence of an inertial manifold // J. Dyn. Differ. Eq. 1991. V.3, №3. P.437-456.

диссипативного полупотока в (конечномерном или бесконечномерном) евклидовом фазовом пространстве X – это гладкая или липшицева инвариантная (интегральная) конечномерная поверхность H X, экспоненциально и равномерно притягивающее при большом времени каждое ограниченное множество D X.

Данное многообразие содержит в себе глобальный аттрактор, т.е. совокупность всех ограниченных целых (существующих при любом времени t (-,+) ) траекторий соответствующей динамической системы. В частности, ИМ содержит все точки покоя и замкнутые траектории данной системы.

& & Показано, что в реальных пространственных координатах (r,r,) -устойчивое периодическое движение летательного аппарата оказывается периодическим или условно-периодическим в зависимости от определённых соотношений между параметрами механической системы.

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

В разделе 2.2 применен метод инерциальных многообразий для трёхмерных систем дифференциальных уравнений с выделенной симметричной линейной частью.

Обычно, инерциальные многообразия рассматривались в связи с изучением финальной динамики нелинейных параболических уравнений в частных производных. Между тем, это понятие оказалось полезным и при изучении систем обыкновенных дифференциальных уравнений первого порядка, так как позволяет сводить задачу описания финальной динамики системы с n степенями свободы к m соответствующей задаче для индуцированной системы с степенями свободы ( m m < n, где размерность ИМ). Долгое время единственной работой, в которой рассматривался подобный подход, была статья Р. Смита3. В данной работе теория инерциальных многообразий была использована для сведения изучаемой там трёхмерной динамической системы на двумерное инерциальное многообразие, что дало возможность использовать теорию Пуанкаре Бендиксона для обнаружения устойчивых предельных циклов. Технику инерциальных многообразий в конечномерных задачах динамики можно рассматривать как плодотворное развитие хорошо известных методов теории интегральных многообразий4.

Для класса рассматриваемых в данном разделе систем ОДУ в случае нечётного векторного поля доказана центральная симметричность периодических орбит (лемма 2.1).

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

1. Согласно интегралу кинетического момента материальная точка в центральном (в том числе, гравитационном) поле движется по плоской кривой.

Smith R.A. Poincare index theorem concerning periodic orbits of differential equations // Proc. London Math. Soc., (3).

1984. V.48, №2. P.341-362.

Митропольский Ю.А., Лыкова О.Б. Интегральные многообразия в нелинейной механике. М.: Наука, 1973.

2. Спутник массы m моделируется материальной точкой, на которую действуют F Fr управляющие силы с радиальной и трансверсальной компонентами.

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

4. Относительное движение объекта исследования либо криволинейное с малой скоростью и малой кривизной траектории, либо прямолинейно (с постоянной скоростью), так что можно пренебречь переносной и кориолисовой силами инерции.

В предположениях 1-4 динамика летательного аппарата описывается уравнениями && & && & & m r - r2 = Fr, m r + 2r = F ( ) ( ). (2.1) При этом ставится обратная задача выбора радиальной и трансверсальной управляющих сил, обеспечивающих устойчивое периодическое движение спутника & & в фазовых переменных (r,r,), где (r,) - полярная система координат в плоскости его обращения.

& & & & F (r,r,) определены соотношениями Fr (r,r,) Пусть управляющие силы и 2 - r2 -( ) ( ) & & Fr = m µ2r & µ1 + µ2 (r + µ2r) + g , & & F = m - µ3r) + r (2r & µ1,µ2,µ3 > 0 g() с постоянными и некоторой гладкой функцией. Тогда, произведя линейную замену переменных & & x1 = r + µ2r x2 = r x3 = ,,, (2.2) исходные уравнения движения спутника (2.1) можно записать в виде трёхмерной системы & -µ1x1 + g(x3), x1 = -µ2x2 + x1, &(2.3) x = & x3 = -µ3x3 + x2, которая является основным объектом исследования во второй и третьей главах работы.

µ1, µ2, µ3 и допустимой Каждому выбору управляющих параметров & g() управляющей функции соответствует вариант управления спутником.

В разделе 2.4 сформулирована основная теорема данной главы. В ней g(x3) получены аналитические условия на управляющую функцию и управляющие µ1,µ2,µпараметры, при которых для системы (2.3) существует устойчивый предельный цикл. Тем самым установлена возможность устойчивого & & периодического движения спутника в координатах (r,r,).

g(x3) ТЕОРЕМА 2.2. Пусть гладкая управляющая функция в системе уравнений (2.3) удовлетворяет условиям:

0 < g(x3) mи g x( ) 1для x (-,+);

3 g( ) > 0 : g ( ) = -hи< 0 = µ1µ2µ3.

µ1, µ2, µТогда, если для упорядочения 1 2 < 3 управляющих параметров справедливы оценки h + 123 > 1 + 2 + 3 1 2 + 1 3 + 2 3 3 - 2 > ( ) ( ),, (2.4) то система (2.3) обладает устойчивым предельным циклом. Таким образом, соответствующие радиальная и трансверсальная управляющие силы 2 - r2 - ( ) ( ) & Fr = m, µ2r & µ1 + µ2 (r& + µ2r )+ g & & F = m(2r - µ3r) + r2 обеспечивают для спутника с массой m существование устойчивой замкнутой & & орбиты в фазовом пространстве (r,r,).

При этом последовательно доказывается:

g(v) g(v) системы (2.3);

x* =,, v а) единственность точки покоя µ1 µ1µ2 б) её неустойчивость (критерий Рауса Гурвица);

в) существование двумерного инерциального многообразия в R с последующим применением теории Пуанкаре Бендиксона.

На рис.2,3 представлены изображения двух различных приближающихся к устойчивому предельному циклу траекторий движения спутника в центрированных Рис.2 Рис.фазовых переменных & & x1 = r + µ2r - x2 = r - x3 = (2.5) 2µ1, 2µ1µ2, 2µ1µ2µс управляющей функцией & & g() = arcctg( - ) =,, (2.6) 2µ1µ2µединственной стационарной точкой x = (0,0,0) и значениями управляющих µ1,µ2 = 0.05 µ2 = 2.параметров,.

Перечислены также возможные альтернативные по сравнению с (2.6) варианты g(x3) выбора управляющей функции в системе уравнений (2.3).

Движение спутника рассмотрено в трёх различных системах координат, а & & именно: в физически значимых координатах (r,r,), обобщённых координатах (2.2) и (для случая управляющей функции (2.6)) в центрированных координатах (2.5).

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

& & В разделе 2.5 рассматривается вопрос о том как выглядит (r,r,) -устойчивое T -периодическое движение летательного объекта в реальных полярных x2(t) = r(t) координатах (r,), действующих в его плоскости обращения. Поскольку в системе уравнений (2.3), то соответствующая радиальная компонента данного движения T -периодична.

ТЕОРЕМА 2.4. Пусть движение спутника в координатах (r, ) отвечает его 2 = устойчивому T -периодическому (с частотой ) движению в переменных T & & (r,r,), существующему при условиях теоремы 2.2 на управляющие силы & & & & & Fr = mµ2r2 - r2 -( ) ( ) µ1 + µ2 (r + µ2r) + g , F = m(2r - µ3r) + r2.

T = Если число, где = (t) dt - средняя за период угловая скорость & T движения спутника, иррационально, то его плоское движение будет условно периодическим. В случае рационального числа данное движение будет k k,l = периодическим. Если при этом, где целые положительные взаимно l lT k простые числа, то за время космический аппарат совершит оборотов вокруг lT центра обращения и его плоское движение будет периодическим с периодом.

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

g(x3) > > Условия, в условиях основной теоремы 2.2 нужны лишь для того, чтобы обеспечить (необходимое в данной физической задаче, связанной с & вращением спутника) сохранение положительности переменных r(t), (t) при возрастании времени. Если же абстрагироваться от этой ситуации и ставить вопрос о существовании устойчивого периодического движения у произвольной трёхмерной системы уравнений вида (2.3), то можно сформулировать следующий результат.

g(x3) ТЕОРЕМА 2.5. Предположим, что гладкая функция в системе дифференциальных уравнений (2.3) удовлетворяет условиям:

g(x3) m g (x3) и на R ;

R : g ( ) = -h < 0 µ1µ2µ3 = g() и.

1 2 < 3 µ1, µ2, µТогда, если для упорядочения величин справедливы оценки (2.4), то система (2.3) обладает устойчивым предельным циклом.

Система уравнений (2.3) сводится к скалярному уравнению третьего порядка y + (1 + 2 + 3 ) y +( + + ) y +12 y =( y) g (2.7) 1 2 1 3 2 3 y(t) = x3 (t) 1, 2, 3 > относительно неизвестной с параметрами.

ТЕОРЕМА 2.6. В условиях теоремы 2.5 на положительные параметры 1 2 < 3 g(x3) и гладкую функцию уравнение вида (2.7) обладает периодическим y(t) решением.

В главе 3 для описанной выше математической модели движения спутника с управляющей функцией вида (2.6) решена задача аналитической аппроксимации его & & устойчивой замкнутой орбиты в фазовых координатах (r,r,) и оценки периода обращения. С этой целью использован метод5, 6 аналитической локализации периодических решений дифференциальных уравнений с полиномиальными нелинейностями, а также систем таких уравнений. Суть метода состоит в замене искомого периодического решения начальным отрезком ряда Фурье с неопределёнными коэффициентами, а также значением периода, и последующим нахождением указанных параметров из соответствующих систем нелинейных алгебраических уравнений. Данная процедура соответствует общей идеологии метода гармонического баланса и даёт хорошие результаты в случаях, когда удаётся сравнить приближённые и точные периодические решения.

Соответствующий подход обладает рядом преимуществ по сравнению с другими возможными: например, c методом Крылова–Боголюбова– Митропольского, вариационными методами, а также различными методами теории возмущений и теории бифуркаций. Во-первых, метод гармонического баланса не связан с наличием в системе малого параметра и позволяет находить периодические решения с большой амплитудой, а во-вторых, оказывается весьма простым в применении, хотя при его реализации возникают достаточно сложные системы нелинейных алгебраических уравнений относительно коэффициентов Фурье и периода приближённого периодического решения.

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

Delamotte B. Nonperturbative (but approximate) method for solving differential equations and finding limit cycles // Physical Review Letters. 1993. V.70, №22. P.3361-3364.

Poland D. Loci of limit cycles // Physical Review E. 1994. V.49, №1. P.157-165.

Ранее метод гармонического баланса обычно применялся к нелинейным дифференциальным уравнениям второго порядка. В данной работе исходная система уравнений (2.3) сводится к одному уравнению третьего порядка, которое после дифференцирования превращается в уравнение четвёртого порядка с квадратичной нелинейностью. Такое повышение порядка приводит к существенному усложнению возникающих систем алгебраических уравнений относительно коэффициентов Фурье и периода искомого приближения. С другой стороны показано, что в аппроксимативных отрезках ряда Фурье достаточно учитывать лишь нечётные гармоники и, тем самым, ускорить сходимость соответствующей вычислительной процедуры. На третьем шаге метода возникает система из семи алгебраических уравнений шестой степени и полученная аппроксимация при выбранных значениях управляющих параметров движения спутника даёт удовлетворительную точность 3.6%.

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

В разделе 3.3 после перехода к центрированным координатам (2.5) уравнения движения спутника (2.3) с управляющей функцией (2.6) приобретают вид & x1 = - µ1x1 - arctg x3, x = - µ2x2 + x1, &(3.1) x = - µ3x3 + x2.

& Данная система сводится к скалярному дифференциальному уравнению третьего порядка, которое после дифференцирования превращается в уравнение четвёртого порядка 1+ y2 y(4) + 1y + 2y + 3y + y = (3.2) ( )( ) y(t) = x3(t) с коэффициентами относительно 1 = µ1 + µ2 + µ3 2 = µ1µ2 + µ1µ3 + µ2µ3 3 = µ1µ2µ,, (3.3) и квадратичной нелинейностью.

В разделе 3.4 описывается применение метода гармонического баланса к уравнению (3.2). Показано, что в данной ситуации можно ограничиться начальными N отрезками ряда Фурье с нечётными гармониками. На шаге приближённое периодическое решение уравнения (3.2) ищется в виде N y(t) = bk sin (2k -1)t + ck cos(2k - 1)t ( ).

k =(N ) x1(t), x2(t) Компоненты приближённого периодического решения X (t) системы x3(t) = y(t) (3.1) легко восстанавливаются по третьей компоненте. В качестве меры правдоподобия данного решения предлагается относительная невязка (t) (N ) o = , где (t) - невязка X (t) относительно (3.1) и равномерная y(t) норма периодических функций. Фактически, (t) - это невязка первого из уравнений (3.1), тогда как два других уравнения удовлетворяются точно.

В разделе 3.5 процедура аналитической аппроксимации искомого периодического движения спутника численно реализована при значениях µ1,µ2 = 0.05, µ2 = 2.1, обеспечивающих согласно теореме управляющих параметров 2.2 существование устойчивого предельного цикла для системы уравнений (2.3), а значит и для индуцированной системы (3.1). В качестве пороговой точности o = 0.05 N =1,2,принимаем. В результате трёх шагов метода ( ) получаем (1) (2) (3) o последовательные приближения X (t) X (t) X (t). Заданная точность,, ( ( удовлетворяется лишь на третьем шаге, а именно:,, o1) = 0.150 o2) = 0.0(3) (3) o = 0.036. Приближённое T -периодическое решение X (t) системы уравнений (3.1) имеет вид x1 t = -0.352sin t +3.029cos t -0.006sin 3 t +0.176cos3t -0. 0005sin 5t + 0.027cos5t ( ) ( ) x t = 6.471sin t +1.480cos t +0.129sin 3 t +0.009cos3 t +0. 012sin 5t + 0.0005cos5t x t = 3.089sin t +0.033cos t +0.045sin 3 t -0.025cos3 t +0.003sin 5t - 0.003cos5t ( ) (1) (2) (3) = 0.4с и T = 2 =13.760. Периоды приближений X (t) X (t) X (t),, практически не зависят от номера шага. Именно: T1 =13.630, T2 =13.688, T3 =13.7, так что значение T3 можно принять за удовлетворительную оценку периода X (t) точного периодического решения системы (3.1).

(3) X t На рис.4 приведен график невязки (t) приближённого решения ( ).

Рис.На рис.5 представлены три последовательные аналитические приближения (1) (2) (3),, X (t) X (t) X (t), а также предельный цикл X (t), построенный численными методами в пакете MAPLE 11. Видно, что амплитуда аппроксимирующих орбит увеличивается с ростом номера шага.

(1) (2) (3) Рис.5. Решения X (t) X (t) X (t) и предельный цикл X (t).

,, (3) В реальных полярных координатах (r,) приближённому решению X (t) системы (3.1) соответствует аппроксимативная траектория летательного аппарата:

r(t) = 628.32 + 6.471sint +1.480cost + 0.129sin3t+ 0.00 9cos3+ t + 0.012sin5t + 0.0005cos5t, (t) = 6.793 + 299.2 t - 6.759cos t + 0.072sin t - 0.0328cos3 t - 0.0182sin3t - 0.0013cos5t- 0.0013sin5t = 0.4со значением частоты и (0) = 0.

1 = 2 = 3 = 2 + 2 0 < < = 0.0В разделе 3.6 полагаем,,, для µ1,µ2,µупорядочения 1 2 < 3 управляющих параметров. В этих условиях теорема 2.2 гарантирует существование устойчивого предельного цикла для системы уравнений (3.1). Первый шаг (N =1) метода гармонического баланса для (1) (3.1) приводит к приближённому периодическому решению X (t) вида x1(t) = -bµ1(µ2 + µ3 )sin t + b( µ2 + µ3)cos t x (t) = bµ3 sin t + bcos t (3.4) x (t) = bsin t - = с и b = 2 -1, где =12 -3, коэффициенты 1,2,3 определены соотношениями (3.3) и 0 < <1. Кривая (3.4) представляет собой эллипс в (1) центрированных координатах (2.5). При этом оказывается, что для X (t) o 0 относительная невязка как при так и при. Итак, в данных µ1,µ2,µпредположениях на поведение величины (а значит, и на параметры ) уже на первом шаге метода получаем хорошее (эллиптическое) приближение искомого предельного цикла.

В разделе 3.8 показано, что метод гармонического баланса для приближённого аналитического определения устойчивого периодического движения спутника можно применить в случае двух альтернативных по сравнению с (2.6) вариантов g(x3) выбора управляющей функции в системе уравнений (2.3):

c = g(x3) =c -sin(x3 - ) c >1) с, ;

µ1µ2µ8 g(x3) = = = 2) с,.

3 3(1 + (x3 - + ) ) 3µ1µ2µg(x3) В случае 1) функция эффективно приближается многочленами Маклорена, а в g (x3) случае 2) производная рациональна. Таким образом, в обеих ситуациях систему (2.3) можно свести к скалярному дифференциальному уравнению с полиномиальной нелинейностью. При выбранных значениях , устойчивый предельный цикл системы (2.3) существует для наиболее широкого диапазона µ1,µ2,µ3 > управляющих параметров.

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

В разделе 4.2 рассматривается обратная задача определения значений интенсивности слоя диполей на поверхности крыла по заданной нормальной компоненте результирующей скорости потока на данной поверхности (условие непротекания). Граничное условие на несущей поверхности для потенциала x, y, z возмущённых скоростей ( ) с учетом постоянной скорости набегающего V потока имеет вид:

+V sin = или, для крыла расположенного в плоскости O x y :

n z x - x1 (x1, y1) dx1 dy1 ( ) -V sin = y - y1) + z2 ( x - x1) - ( y - y1) - z, (4.1) 2 2 2 2 n x ( где - часть поверхности крыла (рис.6) лежащая внутри обратного конуса Маха с 2 x, y, z ( ) вершиной в точке, = M -1, M - число Маха набегающего потока.

Таким образом, задача сводится к решению интегро-дифференциального уравнения7 (4.1) относительно функции распределения интенсивности слоя диполей на поверхности крыла.

Градиент потенциала возмущенных скоростей z x - x1 (x1, y1) dx1 dy1 ( ) x, y, z = ( ), y - y1)2 2 x + z2 x - x1 2 - y - y1 2 - z( ( ) ( ) если найдена функция, определяет скорость течения газа в каждой точке области возмущенного потока внутри конуса Маха (рис.6), а также на поверхности крыла.

Рис.Окончательный шаг в решении обратной задачи состоит в определении p(x, y, z) по формуле Бернулли распределения давления 2 p - p = -V - + ( ) ( ) , x 2 y z p где плотность газа, давление в бесконечно удалённой точке.

Величина коэффициента давления на крыле определяется из линеаризованного Cp = уравнения Бернулли.

x В разделе 4.3 интегро-дифференциальное уравнение (4.1) решается методом сплайн-коллокации, т.е. аппроксимируется системой линейных алгебраических Сирс У.Р. (ред.). Общая теория аэродинамики больших скоростей М.: Воениздат, 1962.

уравнений. Эта система получается разделением поверхности на малые Pi j прямоугольные элементы в предположении, что скачок потенциала представляется в виде многочлена 2-го порядка на каждом элементе. При этом значения производной потенциала в каждом узле коллокации можно получить в n Pi j аналитическом виде, суммируя по всем элементам (внутри обратного конуса Маха), на которые разделена несущая поверхность.

Точно удовлетворяя уравнение (4.1) в узлах коллокации, получаем m n систему линейных алгебраических уравнений. Условия построения двумерного C параболического сплайна класса дают дополнительные соотношения и замыкают систему.

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

Распределение давления (формула Бернулли), поле скоростей и коэффициент давления могут быть получены через градиент потенциала возмущенных скоростей (x, y, z) после построения сплайна для функции распределения интенсивности слоя диполей на поверхности крыла.

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ, ВЫНОСИМЫЕ НА ЗАЩИТУ В работе представлены следующие основные результаты, имеющие научное и прикладное значение.

1. С использованием классических методов обратных задач динамики решена задача о движении подводного аппарата по винтовой линии в предположении, что его центр масс смещён относительно оси динамической симметрии. Найдены условия существования желаемой траектории. Показано, что заданное движение достигается при определенных углах отклонения элеронов.

2. С помощью техники инерциальных многообразий установлено существование устойчивого предельного цикла для класса систем нелинейных обыкновенных дифференциальных уравнений и показано, что при соответствующем выборе функций управления такие системы описывают плоское движение летательного аппарата вокруг некоторого центра вращения. Тем самым установлена возможность облета спутником исследуемого объекта по плоской траектории в реальном трёхмерном пространстве. Данное движение является T & & периодическим и устойчивым в переменных (r,r,), где (r,) - полярная система координат в плоскости вращения летательного аппарата. В исходных пространственных координатах движение спутника оказывается периодическим или условно-периодическим в зависимости от явно выписанных аналитических соотношений между величиной периода T и средней за период угловой скоростью.

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

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

Публикации в журналах перечня ВАК 1 Галиуллин И.А., Кондратьева Л.А. Движение подводного аппарата по пространственной спирали // Вестник Московского авиационного института.

2007. Т.14, №2. С.41-46.

2 Кондратьева Л.А. Обратные краевые задачи на многообразиях // Вестник Российского университета дружбы народов, Серия Математика. Информатика.

Физика. 2010. №1. С.34-38.

3 Галиуллин И.А., Кондратьева Л.А. Спутниковые инерциальные многообразия и предельные циклы // Космонавтика и ракетостроение. 2011. №3(64). С.73-76.

4 Кондратьева Л.А. Приближённое аналитическое вычисление устойчивой периодической орбиты спутника // Вестник Московского авиационного института. 2012. Т.19, №1. С.75-80.

Публикации в других изданиях 5 Кондратьева Л.А. Расчёт аэродинамических характеристик при сверхзвуковом обтекании прямоугольного крыла // Сб. научн. трудов: «Экспериментальное и теоретическое исследование аэродинамических характеристик ЛА и его частей».

М.: Изд-во МАИ, 1983. С.26-28.

6 Кондратьева Л.А. Аппроксимационная обратная задача для предельных циклов // Качественное и численное исследование математических моделей динамических систем. Межвуз. сб. научн. трудов. М.: Изд-во РГОТУПС, 2006. С.72-75.

Доклады на научных конференциях 7 Галиуллин И.А., Кондратьева Л.А. Тезисы доклада «К вопросу о программном движении подводного аппарата», Труды XIV международного научнотехнического семинара «Современные технологии в задачах управления, автоматики и обработки информации». Алушта, 2005 г.

8 Кондратьева Л.А. Тезисы доклада «Численные оценки сложности для нелинейных параболических уравнений», Труды XI международного научнотехнического семинара «Современные технологии в задачах управления, автоматики и обработки информации». Алушта, 2002 г.

9 Кондратьева Л.А. Тезисы доклада «Конечномерная динамика эволюционных уравнений и нелинейный метод Галеркина», Труды XII международного научнотехнического семинара «Современные технологии в задачах управления, автоматики и обработки информации». Алушта, 2003г.

10 Кондратьева Л.А. Тезисы доклада «Аппроксимационная обратная задача для предельных циклов», Труды XV международного научно-технического семинара «Современные технологии в задачах управления, автоматики и обработки информации». Алушта, 2006 г.

11 Кондратьева Л.А. Тезисы доклада «Обратные задачи в теории сплайн-функций», Труды XVI международного научно-технического семинара «Современные технологии в задачах управления, автоматики и обработки информации».

Алушта, 2007 г.

12 Кондратьева Л.А. Тезисы доклада «On the Loci of Limit Cycles for the Spacecraft», Sixth International Symposium on Classical and Celestial Mechanics, Moscow – Velikie Luki, 2007 г.






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.