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

Министерство науки, высшей школы и технической политики

Казанский ордена Трудового Красного Знамени и ордена Друябы народов Государственный технический университет им. А.Н.Туполева (КГТУ),

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

ХУЗЯТОВ ШЙФИК ИАЕХОВИЧ

УСТОЙЧИВОСТЬ СИСТЕМ С РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ С ЗАПАЗДЫВАЮЩИМ АРГУМЕНТОМ И ПРИЛОЖЕНИЕМ К ИССЛЕДОВАНИЮ УСТОЙЧИВОСТИ УПРУГИХ КОЛЕБАНИЙ КРНЛй

Специальность: 05,13.01 - Управление в технических системах

Автореферат диссертации на соискание ученой степени кандидата технических наук

Казань - 1993

Работа выполнена в Казанском ордена Трудового Красного Знамени и ордена Дружбы Народов Государственном техническом университете им. й,Н.Туполева.

Научный руководитель - Заслуженный деятель науки и техники

РФ и РТ, академик АНТ, доктор технических наук, профессор Т.К.Сираэетдинов

Официальные оппоненты: доктор технических наук,

профессор Ф.Д.Байрамов, кандидат технических наук, доцент Г.Г.Бильченко,

Ведущее предприятие : Опытно-конструкторское бюро "Сокол",

Защита состоится "¿1" г. в {2. часов

на заседании специализированного совета К.063.043.03 Казанского ордена Трудового Красного Знамени и ордена Дружбы Народов Государственного технического университета им. А.Н.Туполева, 420084, Казань, ул. К.Маркса, д. 10.

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

Автореферат разослан _______199^_г.

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

Л.0.Шацилло

ОБЩЙЯ ХАРАКТЕРИСТИКА ДИССЕРТАЦИОННОЙ РАБОТЫ

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

Одним из основных методов исследования устойчивости процессов с распределенными параметрами и с запаздывающим аргументом является метод функций Ляпунова, Разработкой метода функций Ляпунова и исследованием вопросов устойчивости систем с распределенными параметрами с запаздывающим аргументом занимались H.H. Красовский, В.М. Матросов, Б.С. Разумихин, В,В. Румянцев, A.A. Мовчан, Т.К. Сиразетдинов, Г,Л. Дегтярев, П.К. Ванг, ü.M. Зайцев, П.К. Семенов, Ф.Д, Байрамов, В,Б, Кол-мановский и многие другие исследователи. Но тем не менее исследования устойчивости как систем с распределенными параметрами, так и с запаздывающим аргументом остается актуальной научной проблемой при проектировании современных инженерных объектов.

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

параметров.

Таким образом, недостаточная разработанность метода Функций Ляпунова для исследования устойчивости таких систем задерживает решение многих прикладных задач, что и обусловило актуальность темы диссертационной работы,

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

Научная новизна заключается в следующем:

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

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

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

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

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

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

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

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

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

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

Реализация и внедрение. Результаты диссертационной работы использованы при улучшении динамических характеристик двигателей внутреннего сгорания в НТЦ КамАЗа (г.Наб. Челны),' при исследовании колебаний упругих тел и оболочек в институте механики и маииностроения РАН, в учебном процессе в КамПЙ, о чем имеются соответствующие акты.

Апробация работы. Основные результаты диссертации были долоаены и обсуждены на и-й Республиканской научно-технической конференции.КамАЗ-КамПИ (г,Наб. Челны, 1986 г.), Международной математической конференции "Ляпуновские чтения" (г.Харьков, 1992 г.), научно-технической конференции "Научный потенциал вузов - программе "Конверсия" (г.Казань, 1993 г.). Отдельные результаты работы были использованы при составлении двух отчет.ов плановых научно-исследовательских работ, проводимых в АНТ.

Публикации. По результатам выполненных исследований опубликованы 3 научные статьи, 3 тезиса доклада научно-технических конференций, 3 научно-технических отчета.

Структура и обьем работы. Диссертация состоит из введения, трех глав, заклвчения, двух приложений, списка литературы из 98 наименования и 7 рисунков. Полный объем диссертации

составляет страниц.

содержание работ а.

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

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

Рассматривается мноаество F/i ~ Рц (-В) измеримых вектор-функций <f(x)— . • • , (/>„М)Т с равномерным

приращением, т.е. удовлетворяющие неравенству:

sup I<f. 1*0 - у. IX"JI < м (t'bnj

где Л - ограниченная односвязная область в евклидовом пространстве Rs , М - заданное положительное число. Этому классу принадлежат класс равномерно ограниченных функций, класс равностепенно непрерывных функций. Сам класс №) принадлежит пространству /г(2)) суммируемых функций.

Для исследования устойчивости возмущенного процесса вводится мера отклонения процесса;

f Iftxti = fix) + <f*(rj , ^ If] = / <>[f(r)]cl* .

Рассматривается функция такая, что l}(0, x, 0;

xe JJ, и множество вектор-функций <f(x) e F«

таких, что функция

является суммируемой в J7 при всех f-> А» , На множестве функций <f>txj введем

функционал:

Ztyj]^ f mm x,i] dx

Определение 1. Функция J) i] называется определенно

положительной по мере $[<f(xj] почти всюду в D , если для лвбого положительного числа £ 4 h, h>0 и для всех векторов-функций , удовлетворяющих неравенству

4. h, хе 2? , каким бы малым числом yi>0 не задавались, существует число <Г = Т(£,р)>0 такое, что мера множества

X6U} #[(f(X), <

удовлетворяет неравенству (Uts-)</* ПРИ всех

Справедлива следующая теорема.

Теорема 1, Для определенно половительности функционала V3 fyi] по мере f [f] для функций tf(x) , принадлежащих классу Ft> необходимо и достаточно, чтобы функция iJ была

определенно полоаительна по мере f>L<f(*)] почти всшду в Л для всех функций </(х/ из области f^l^J^R, R > О-

Далее рассматривается некоторое множество вектор-функций <sj = [(/> fas), (fatysj) , заданных на множестве G-X^S

где $ = о] . Предполоким, что при каждом фиксированном se £ вектор-функция. f(x,s) принадлежит классу функций с

равномерным приращением и при каждом фиксированном хе Л принадлежит классу /f равностепенно непрерывных функций, т.е. yiz,s)e />х/(.

На множестве х К функций <f(*,s) введем функционал о

Кг i(f> ^^ ^ <t* ds

'г Л

Теорема 2. Для определенно положительности функционала ^[tf,!] по меРе Для Функций tf!^s)e Ft)-у К достаточно,

чтобы функция г^/у, х, s/ /J была определенно положительна по мере q [if] почти всюду в С = Лх £ для всех функций <f(*>s) из области К , R>0.

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

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

Теорема 5. Для определенно положительности квадратичной Формы у по мере почти всюду в Л

необходимо, а в случае ограниченности элементов I) и достаточно, чтобы существовала положительная почти всюду в 3) функция ^(х,1) такая, что при всех выполнялись неравенства:

Ъ; <*, V ' 2 О-"Л,

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

Формулам:

¿„<т,4)=±ыиы)-2: См)"*,

И С/ 1)~ 0 ПРИ С - ¿,п ; ¡<С .

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

рекуррентным критериям.

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

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

(Л.« Й + $* «>

Е (ем + !=, Ы) + Шч + ?. Му = О, (г)

хе Т, 10,

где - п-мерный вектор фазовых функций, ^ =

- м-мерный вектор фазовых функций, кк(х), Вк (х! , £к(х]^ Рк(х)

0; (х) матрицы, элементы которых принадлежат

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

Граничные условия задаются на некоторой части "ЗД, границы И Л области Л в виде:

Л !р(х, {/ у>(х, о = О , С 21) Ш

где р = р Иг) - матрицы, элементы которых

непрерывные ограниченные функции времени.

При исследовании устойчивости нулевого решения системы У) - (з) методом функций Ляпунова имеется два различных способа построения функций Ляпунова.

Первый способ построения функций Ляпунова. Для решения вышепоставленной задачи используем функции Ляпунова в виде интегральной квадратичной формы

И, ЦЬ V] = / I) * (х) у (X, I) с1х Л

Вычислим производную сН/з /сИ в силу (V - (5) некоторых преобразований приведем его к виду

п..г.... 01

с(1

(V

и после

= УТ:)с1г (5)

где

Л

5

0 Щ^АЬ)_ е.)-(¿А. >Р<

к=/ о

(в)

Здесь и Рг матрицы, которые используются при

вычислении производной с(Ц, /с( I для учета системы

уравнений (2) , которые не содержат производных по времени от Фазовых функций.

Согласно теоремы об устойчивости, нулевое решение системы (У - (3) устойчиво, если функционал Уз [у] определенно положителен по мере и производная с1]/3/с11

неположительна на множестве состояний, удовлетворяющих неравенству : 4 .

Здесь для оценки производной /¿I в силу

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

Лемма. Билинейная форма Ц>г при условии

допускает следующую оценку;

С использованием этой леммы условие неполовительности производной dV3/cH записывается в виде:

Vrt!(*)<t'/¿^ЧшУ]2 J

max --iïi--L—¿i , xeJJ (4)

f (rowf r„\.*<krts*

Максимум этой функции зависит только от направления вектора <f . Поэтому для определения максимума достаточно рассмотреть единичные вектора. При заданных коэффициентах системы ¡1), (2J и элементов матрицы определение

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

Для определения области устойчивости в пространстве параметров по оценке (Ч) на языке Turbo Paskal была составлена программа.

При этом размерность пространства параметров считалась произвольной. Локальный максимум (минимум) в оценке (Ч) определялся методом скорейшего спуска, область предполагаемого глобального максимума определялась методом сеток. Время расчета при решении различных задач на машине PC/fiT 5-40 минут.

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

квадратичных форм, то оценка (Ч) представляется в виде:

- (v

где if - максимальные собственные

числа матриц l^fxj и M'Ux) ^ (х) в точке sea Л .

Заметим, что эта оценка является более удобной в практическом использовании, но несколько грубой оценкой. Как показывают примеры, области устойчивости, построенные с использованием оценки (Ч) и (8) шире, чем области, построенные с использованием оценок, предложенных ранее в других работах.

Второй способ построения .функций Ляпунова. Для получения достаточных условий устойчивости нулевого решения системы (1) - (3) используется функционал:

[у] = Ут(хМ(х) у(х, {¡о[х + £ / ^/(х^х^Ш^ (9)

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

7Г ^ Vх м

/

где и)е = (О - 27 и). .

Пусть квадратичная форма <-рт 1?(х) у определенно

полоиительна почти всюду в Л , квадратичные формы У'ьЗ,-!*)^ ~ 1, ¿) неотрицательны, а квадратичная

форма

/ /

№ • • V у -¿я у ^ * £ ч- ^ м

П I { + 1) ~ переменных является неотрицательной при всех хеЛ • Тогда выполняется условие теоремы об устойчивости нулевого решения системы ¡1) ~{3).

Для определения области устойчивости в пространстве параметров разработаны алгоритмы и программы. Условие неотрицательности квадратичной формы , ..^ }

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

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

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

^ - -Ь-г)

оЦг (Ц

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

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

В третьей главе исследуется устойчивость изгибно-крутиль-ных колебаний стреловидного крыла при действии аэродинамических сил, сил упругости и сил воздействия органов управления, некоторые компоненты которых содержат запаздывание. При некоторых предположениях, крыло заменяется эквивалентной балкой; ось которой совпадает с осью жесткости крыла. Изгибно-крутиль-ные колебания крыла рассматриваются как изгиб оси жесткости крыла Ц> = у /2, и закручивания сечений крыла вокруг

оси жесткости на угол I) , где 2 - координата

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

(Ю

Здесь т - погонная масса; (> - расстояние от центра тяжести сечения крыла до оси жесткости; 1Ж - распределенный момент инерции сечения относительно оси жесткости; б

- жесткость крыла на изгиб и на кручение; Уа / Ма распределенные аэродинамические нагрузки: подьемная сила крыла, момент относительно оси жесткости крыла; 1Л - угол отклонения управляющей поверхности крыла; 1У [2) и хв (¿^ задают соответственно распределение управляющих аэродинамических сил и моментов при отклонении органа управления, где г е(г4/¿¿) с (011) • 1^(2) Ы

управляющая сила; И - управлявший момент. Заметим,

что здесь все переменные приведены в безразмерном виде.

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

лявтся в виде линейных комбинаций от

Ъу/д?, управлаящие силы и моменты представляются в виде линейных комбинаций значений б, ^ ^ 2^/р/ в момент

времени 1-т.

Разрешая систему (4 £) относительно высших производных повремени '2гу/9{г, 0*9/21* и вводя новые переменные

* = « * = 4 , ^ = ^ а, _ с г 5 6 ^ к 9*' 7' Ы *№)> ъ

эту систему запишем в виде системы уравнений (1) , (2) . Ери выполнении этих преобразований введена новая система обозначений коэффициентов уравнений, согласно которой коэффициенты уравнения представлены в виде рядов от расстояния . Эта система обозначений значительно облегчает построение функций Ляпунова, удовлетворяющих условиям теоремы об устойчивости и физическую интерпретацию коэффициентов, входящих в уравнения колебаний и функции Ляпунова.

Используя функционалы вида (4) И I^[у] (9)

построены функции Ляпунова,, удовлетворяющие условиям устойчивости, Из условий представления производной в виде (9) непосредственно определяются некоторые элементы матрицы т/(х) или выражаются через другие элементы матрицы Мх): - 1>ц( = ^О, ¡/чч = * /Е1И ,

-*1"' Ф*/!«, Далее, если предполагать

Е1н ~ £ X • 70 из этих же Условий получается

= м , = Заметим, что эти

элементы матрицы и г/^ совпадают с коэффициентами

полной энергии консервативной системы, т.е. системы, описываемой уравнениями (12) в случае Уа — Д/в = Ц Остальные элементы матрицы с некоторой степенью

произвольности определяйся из условий определенно пслокитель-ности матриц гЛ и оЭ .

Используя алгоритмы и программы, разработанные в главе 2. ,

для некоторых численных значений параметров крыла и потока построены области устойчивости в пространстве параметров {¿¿о,^, • где ¿го ~ коэффициент статической устойчивости крыла, // , - параметры, характеризующие структуру управления.

Основные результаты работы :

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

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

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

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

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

В.' Разработаны алгоритмы и программы построения области устойчивости в пространстве параметров уравнений, описывающих колебания упругого крыла.

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

1. Сиразетдинов Т.К., Хузятов 11.В. Об устойчивости систем с распределенными параметрами. ПММ.-1993,-т.5?,вып.6.-с,14-21.

2. Сиразетдинов Т.К., Хузятов 11.1. Необходимые и достаточные условия определенно положительности некоторых функционалов (Часть I) //йвтоматика-1993-т.24, N 4,-с.102-110,

3. Сиразетдинов Т.К., Хузятов в, и. Устойчивость многомерных процессов с распределенными параметрами и запаздывающим аргументом //Автоматика-1993, N 5. -с.65-74,

4. Сиразетдинов Т.К., Хузятов И. И. Устойчивость процессов с распределенными параметрами с запаздывающим аргументом. Международная математическая конференция "Ляпуновские чтения": Тез, докл., Харьков, lS92.-c.146.

5. Хузятов II.И. К устойчивости систем с распределенными параметрами с запаздывающим аргументом. Научно-техническая конференция "Научный потенциал вузов-программе "Конверсия"": Тез. докл., Казань, 1993,-с.82.

6. Хузятов 1Л., Саляхиев И,Г. Аналитическое построение управления для нагрева тонких тел. Научно-техническая конференция "КамЙЗ-КамПЙ": Тез, докл., г.Наб.Челнн, 1980. - с.205.

Формат 60x84 1/16. Бумага д/шож. аппар. Пешть офсетная. Печ.л.1,0. Усл.печ.л.0,93, Усл.кр.-отт.1,16. Уч.-изд.л,1,04. Тираж 100. Заказ 353/Р.

¡казанский государственный технический университет

имени А.Н.Туполева Ротапринт Казанского государственного технического университета имени А.Н.Туполева 420111, Казань, К.Маркса, 10