CC BY
39
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ
УДК 531.36:534.1
МАТЕМАТИЧЕСКИЕ МЕТОДЫ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ ОБЪЕКТА, ДВИЖУЩЕГОСЯ ПО УПРУГОЙ НАПРАВЛЯЮЩЕЙ
© 2008 г. С.Н. Веричев
Нижегородский филиал Института машиноведения им. А.А. Благонравова РАН
8. verichev@infonet.nnov. ги
Поступнла в редакцню 30.04.2008
Работа посвящена обзору математических методов исследования устойчивости сосредоточенных объектов, движущихся вдоль распределенных упругих направляющих. Перечислятся основные этапы математического анализа устойчивости колебаний (для случая линейных моделей). Затем, все эти этапы демонстрируются на конкретном примере осциллятора, движущегося вдоль балки Бернулли-Эйлера.
Ключевые слова: упругая направляющая, осциллятор, линейная модель.
Математическая постановка задачи о взаимодействии движущегося объекта с упругой направляющей содержит уравнение в частных производных, описывающее колебания направляющей и уравнение(я) колебаний объекта [1—3]. Эти уравнения динамически связаны в движущейся точке контакта. Связь отражает баланс поперечных сил, действующих в контакте. Для описания динамических условий в контакте, в правую часть уравнения колебаний упругой системы вводится обобщенная дельта-функция Дирака, описывающая положение движущегося объекта и вследствие этого зависящая как от времени, так и от пространственной переменной.
При исследовании подобных уравнений очень удобным является введение движущейся системы отсчета, которая сопровождает движущийся объект. В этой системе отсчета дельта-функция перестает зависеть от времени, что позволяет сделать заключение об устойчивости системы в пространстве изображений по Лапласу, не переходя в пространство оригиналов. Далее, применение интегральных преобразований Фурье по пространственной перемен-ной(ым) и Лапласа по времени позволяет перейти от системы дифференциальных уравнений в частных производных к системе алгебраических уравнений и в итоге получить характеристическое уравнение для колебаний объекта движущегося по упругой направляющей. По-
скольку устойчивость системы определяется корнями данного уравнения, то дальнейший анализ связан с изучением этих корней. В силу того, что характеристическое уравнение является интегральным (в него входит обратное преобразование Фурье), решение его «в лоб» довольно затруднительно. Поэтому, для исследования корней характеристического уравнения используются метод Б-разбиения и принцип аргумента [4-6]. Метод ^-разбиения, разработанный Ю.И. Неймарком, используется для разделения пространства параметров системы на области с различным числом «неустойчивых» корней (корней с положительной действительной частью). Последующее применение принципа аргумента используется для нахождения числа данных «неустойчивых» корней в одной из полученных областей (для одного набора параметров). Комбинация этих двух методов является очень удачной и может быть использована для исследования различных, в том числе и трехмерных моделей [7, 8].
Рассмотрим последовательно применение данных методов на простом примере одномассового осциллятора массы т и жесткости к, равномерно и безотрывно движущегося вдоль балки Бернулли-Эйлера [9-13] (см. рис. 1). Данная задача позволит наглядно продемонстрировать вышеупомянутые этапы математического анализа устойчивости рассматриваемых систем.
Рис. 1. Равномерное движение одномассового осциллятора вдоль балки Бернулли-Эйлера
Математическая постановка задачи имеет вид
PFutt с EIuxxxx =-S(x - Vt)
d 2u0
\
m-
с kuО
(1)
= -5Є)
Ґ d 2u (О, t)
m
dt2
с
ku (О, t)
Это позволяет перейти от дифференциального уравнения в частных производных (2) к алгебраическому уравнению:
(pF (s2 - 2ikBsV - V2k2) с EIkB ^ =
2
= - (ms2 с k)Hs (О,s),
(4)
которое может быть переписано в виде
"2 с k)s
D(k,s)uk s = - (ms2 сk us (О,s), D(k, s ) =
u0 (t) = u (Vt, t), lim u (x, t)<ад.
V ' V ' \x-Vt\1 V 71
Здесь u° (t) и u (x, t) - вертикальное отклонение массы m и поперечное смещение балки соответственно; E - модуль Юнга; p и I -погонная плотность материала балки и момент инерции сечения балки на поворот; F - шю-щадь поперечного сечения балки; k - жесткость осциллятора; 5 (...) - дельта-функция.
Заметим, что частные производные в первом уравнении системы (1) являются обобщенными производными.
Введем движущуюся систему отсчета
{ = x-Vt, t = t}. В этой системе уравнения
(1) примут вид
PF (utt - 2VutS с V2u%)с EIum? =
= р(я-2 - 2гк^Г - V2кв2V Е1к4 . (5)
Для нахождения замкнутого уравнения относительно фурье-образа смещения массы й(0, я) применим к (4) обратное преобразова-
нне Фурье по волновому числу кв. Полагая в полученном уравнении £ = 0, придем к сле-
дующему уравнению:
( (
ms2 с k с
й(О, s)
dk
-Л
2^-адD(k, s)
(
11 >/77■ J
dk
2* -адD(k, s)
= О . (6)
(2)
Как видно из (2), переход в движущуюся систему координат позволяет уменьшить число полевых переменных путем исключения из
уравнений смещения массы и0, а также избавиться от временной зависимости в аргументе дельта-функции.
Следующий шаг — применение интегральных преобразований Фурье и Лапласа по переменным £ и і:
ад
{и (§, ^ )}=Л (§, і )}єхр (- ^
0
ад
,, (кв,5)}= I & (^5)}ехр(—ікв£У£, . (3)
Последнее уравнение в (6) представляет собой характеристическое уравнение колебаний осциллятора, равномерно и безотрывно движущегося вдоль балки. Корни этого уравнения определяют устойчивость системы. Достаточным условием возникновения неустойчивости вертикальных колебаний осциллятора является наличие хотя бы у одного корня характеристического уравнения (6) положительной действительной части. Легко видеть, что данное уравнение является интегральным по отношению к переменной Лапласа я, и найти его комплексные корни крайне сложно. Поэтому для дальнейшего исследования удобно использовать метод D-разбнення по одному комплексному параметру. Идея метода состоит в отображении мнимой оси комплексной плоскости (я) на плоскость системных параметров, которые временно будут рассматриваться как комплексные. Закон отображения следует из характеристического уравнения, которое должно быть переписано таким образом, чтобы выразить выбранный параметр явно. Полученная посредством данного отображения линия разделяет пространство параметров на области с различным числом корней характеристического уравнения, имеющих положительную действительную часть. Чтобы отобразить мнимую ось я -плоскости на плоскость некоторого комплексного параметра,
2
например к (забудем на время о его физическом смысле), необходимо подставить 5 = /О в уравнение (6) и затем, используя это уравнение как закон отображения, проварьировать О от —да до +да. Согласно (6) закон отображения на комплексную плоскость (к ) имеет вид
{і да л, V1
к = тО —
2л I
ёк
(7)
2л —'да Б(к, О)
Взаимное расположение кривых .О-разбиения на комплексной плоскости (к) показано на рис. 2.
Как видно из рисунка, возможно два типичных расположения данных кривых. Левый рисунок соответствует случаю, когда скорость осциллятора не превышает минимальную фазовую ско-
( /, ™\1/4 А
рость волн в балке
V < у± = 21 —— т рЛ
пра-
выи же соответствует так называемому «закри-тическому» движению со скоростью V > У1фп . Таким образом, осуществив отображение (7), мы получили разбиение комплексной плоскости (к) на области с различным числом корней характеристического уравнения (6), имеющих положительную действительную часть. Штриховка на линиях .О-разбиения показывает направление перехода из левой в правую полуплоскость комплексной плоскости (5). При пересечении кривых О-разбиения в направлении штриховки число N корней характеристического уравнения с положительной действительной частью (число неустойчивых корней) увеличивается на единицу.
Рис. 2а показывает, что при «докритиче-ском» движении осциллятора (V < Уф) устойчивость системы не зависит от жесткости осциллятора (жесткость осциллятора отражена положительной полуосью оси Яе к). Напротив, при закритическом движении, как видно из рис. 2б, существует критическая жесткость осциллятора к = к *. При переходе к через это значение (в направлении возрастания жесткости), количество корней с положительной действительной частью уменьшается на 2.
Итак, с помощью метода О-разбиения, мы ответили на вопрос: зависит ли количество корней характеристического уравнения с Яе 5 > 0 от жесткости осциллятора, и если да, то каким образом. Для завершения анализа задачи нужно определить число N таких корней для произвольного набора параметров задачи. Это удобно сделать с помощью принципа аргумента. Для его применения уравнение (6) должно быть рассмотрено как некая функция комплексного аргумента Р (5 ) :
( 1 “ л, V1
1 ад
Р (5 )= тя2 + к + —| 2л
ёк
(8)
О(к, 5) j
Далее, на плоскости (5) необходимо выбрать замкнутый контур, и отобразить его на комплексную плоскость Р , используя 5 как параметр. Таким образом, можно получить на
плоскости (Р ) замкнутый контур, число вращений которого в положительном направлении вокруг начала координат равно разнице между
числом корней уравнения Р (5) = 0 и числом
Рис. 2. D-разбиение по комплексному параметру к
полюсов данного уравнения, лежащих внутри выбранного контура. Для определения количества корней с положительной действительной частью, контур на плоскости 5 необходимо выбрать в виде полуокружности, представленной на рис. 3.
Рис. 3. Контур на комплексной плоскости 5, используемый для принципа аргумента
1т Р( 5)
кости Р равно числу корней характеристического уравнения с Яе 5 > 0.
Применение принципа аргумента к отображению (8) показывает, что в «докритическом»
режиме, при к > 0 контур Р (5) , определяемый выражением (8), не включает в себя начало координат комплексной плоскости Р. Это означает, что на рис. 2а N = 0, т.е. колебания осциллятора в «докритическом» режиме всегда устойчивы. Напротив, в закритическом режиме, при
0 < к < к *, контур Р (5) дважды обходит начало координат (см. рис. 4). Из этого следует, что на рис. 2б, N + 2 = 2 ^ N = 0, т.е. колебания осциллятора в «закритическом» режиме устойчивы при к > к и неустойчивы при 0 < к < к *.
Таким образом, применение принципа аргумента завершает математическую процедуру анализа устойчивости колебаний системы движущийся осциллятор — балка.
Работа выполнена при частичной поддержке гранта 1т Р( 5)
ЯеР(5) Яе Р(5)
Рис. 4. Кривая в плоскости (1т Р(5), Яе Р(5) ), полученная с помощью принципа аргумента для «закритиче-
ского» случая движения осциллятора для области 0 < к < к *. Данная кривая дважды обходит начало координат. Это означает, что количество «неустойчивых» корней в данной области N = 2
Данный контур должен быть выбран таким образом, чтобы он лежал в правой полуплоскости комплексной плоскости 5 и имел достаточно большой радиус, чтобы включать в себя все корни и полюса уравнения Р (5) = 0 . Легко показать, что уравнение Р (5) = 0 не имеет полюсов в полуплоскости Яе 5 > 0 и, следовательно, число вращений контура на комплексной плос-
РФФИ, № 08-08-97057-р_поволжье_а.
Список литературы
1. Метрикин А.В. Постановка задачи о взаимодействии движущегося точечного объекта и упругой направляющей с использованием обобщенных функций // Волновые задачи механики: Сб. научн. трудов. Нижний Новгород, 1994. С. 32-37.
2. Весницкий А.И., Каплан Л.Э., Уткин Г.А. Вывод естественных граничных условий для одномерных задач динамики с движущимися нагрузками и закреплениями // Дифференциальные и интегральные уравнения. 1982. С. 75-80.
3. Весницкий А.И., Каплан Л.Э., Уткин Г.А. Законы изменения энергии и импульса для одномерных систем с движущимися закреплениями и нагрузками // Прикладная математика и механика. 1983. Т. 47. С. 863-866.
4. Неймарк Ю.И. Динамические системы и управляемые процессы. М.: Наука, 1978.
5. Араманович И.Г., Лунц Г.Л., Эльсгольц Л.Э. Функции комплексного переменного. Операционное исчисление. Теория устойчивости. М.: Наука, 1968.
6. Verichev S.N., Metrikine A.V. Instability of a moving bogie on a flexibly supported Timoshenko beam // J. of sound and vibrations. 2001. V. 253. P. 653-668.
7. Veritchev S.N. Instability of a vehicle moving on an elastic structure // Delft University press. 2002. Р. 192.
8. Metrikine A.V., Verichev S.N., Blaauwendraad J. Stability of a two-mass oscillator moving on a beam sup-
ported by a visco-elastic half-space // International Journal of Solids and Structures. 2005. V. 42. P. 1187-1207.
9. Денисов Г.Г., Новиков В.В., Кугушева Е.К. К задаче об устойчивости одномерных безграничных систем // Прикл. математика и механика. 1985. Т. 49. С. 691-696.
10. Крысов С.В. Вынужденные колебания и резонанс в упругих системах с движущимися нагрузками: Учебное пособие. Горький: Изд-во ГГУ, 1985.
11. Лисенкова Е.Е. Некоторые общие соотношения для упругих волн, возбуждаемых в направляющей движущимся объектом // Волновые задачи механики: Сб. статей. Нижний Новгород: Нф ИМАШ РАН, 1991. С. 44-52.
12. Метрикин А.В. Неустойчивость поперечных колебаний объекта, равномерно движущегося вдоль упругой направляющей как следствие аномального эффекта Доплера // Акустический журнал. 1994. Т. 40. № 1. С. 99-103.
13. Метрикин А.В. Резонанс в системе упругая направляющая - движущийся осциллятор // Акустический журнал. 1994. Т. 40. № 6. С. 974-978.
MATHEMATICAL METHODS TO STUDY THE STABILITY OF AN OBJECT MOVING ALONG AN ELASTIC GUIDE
S.N. Verichev
Mathematical methods used to study the stability of lumped objects moving along distributed elastic guides have been reviewed. The main stages of mathematical analysis of the stability of an object moving along an elastic guide are listed (for the case of linear models). All these stages are then demonstrated by the example of an oscillator moving along an Euler-Bernoulli beam.
