Научная статья на тему 'Асимптотические формулы для частот осесимметричных колебаний оболочки вращения'

Асимптотические формулы для частот осесимметричных колебаний оболочки вращения Текст научной статьи по специальности «Физика»

CC BY
35
8
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Асимптотические формулы для частот осесимметричных колебаний оболочки вращения»

Статья публикуется по приглашению редакции

Асимптотические формулы для частот осесимметричных колебаний оболочки вращения.

Асланян А.Г.(1), Лидский В.Б. (2)

(1)-Московский институт радиотехники, электроники и автоматики (2)-Московский физико-технический институт

Рассматриваются осесимметричные колебания тонкой упругой оболочки вращения, которые описываются системой уравнений (см.[1,с. 104-136]) - и' + аи + а2 и + а3 и' + а4 и = Яп

± И 2(и(у) + Ь1 ит + Ь2 и" + Ь3 и') + схи + с2 и + с3 и = Яи

(1)

Здесь дифференцирование ведется по длине дуги меридиана оболочки (¿1 < ^ < 52); п(5), и(5) - компоненты вектора перемещения точки срединной поверхности (п(5) - перемещение вдоль меридиана, и(5) - по нормали); Я -спектральный параметр, отличающийся постоянным множителем от квадрата собственной частоты О, И - толщина оболочки (малый параметр).

Функции а(5), Ь(5), с}- (5) определяются геометрией оболочки; их явный

вид можно найти, например, в [1, с. 136].

Дифференциальный оператор, порожденный левой частью системы (1), формально самосопряжен. На каждой граничной параллели 5 = 51 и 5 = 52 ставятся три условия, каждое из которых совпадает с одним из условий жесткого закрепления, либо свободного края. Подробности см . [1].

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

Так возникают 64 самосопряженные краевые задачи, порожденные системой (1), каждая из которых имеет отрицательный дискретный спектр, а соответствующие собственные функции (с.ф.) образуют ортонормированный базис [1]. Оператор, порожденной системы (1) и одним из вариантов указанных граничных условий, обозначим через 1к. Численное определение собственных значений с.з. и с.ф. оператора 1И является не простой задачей, особенно при малых значениях толщины оболочки И. Дело в том, что спектр оператора 1к

образует неоднородное множество со сложной структурой, при этом с уменьшением толщины оболочки И плотность распределения частот возрастает.

Выделяют три области изменения спектрального параметра Я, в каждой из которых структура спектра может быть описана (см. [1]): [0,а1) - область 1;

[а1, ] - область 11; (Д - область 111. Здесь [а1, ] - отрезок значений функций

^(5)=(1 -а2У-2(5),< 5 < ^2, (2)

о - коэффициент Пуассона, У2 (5) - один из главных радиусов кривизны оболочки (он равен расстоянию по нормали от точки срединной поверхности до оси вращения).

В области 1 с.з., как правило, расположены редко и слабо зависят от к. Однако, могут иметь место и крайние случаи: в указанной области с.з. могут отсутствовать вовсе или их (при малых к) может быть достаточно много [2]. Это определяется геометрией оболочки и свойствами безмоментной задачи, в которую вырождается оператор ¡к при к ^+0 (подробности см. [1,3]). Безмоментная задача проще, чем исходная моментная: порядок системы на четыре меньше и отсутствует малый параметр при старшей производной. Отметим, что безмоментный оператор 10 также является самосопряженным (см. [1]).

В области 1 вырождение моментной задачи в безмоментную регулярно в смысле М.И. Вишика - Л.А. Люстерника (см. [4]). Это позволяет написать асимптотические формулы для Х^ (к), в которых главные члены суть с.з.

безмоментной задачи.

Сложнее устроен спектр в областях 11 и 111. Так, 11 - зона непрерывного спектра безмоментного оператора 10. В ней при малых к наблюдается наивысшая концентрация частот моментной задачи, соответствующих преимущественно изгибным (квазипоперечным) колебаниям (см. [1, с. 107]). Зона 111 содержит бесконечно много с. з. Они отвечают как квазипоперечным, так и квазитангенциальным колебаниям. Последние слабо реагируют на уменьшение к, в то время как «квазипоперечные» с.з. с уменьшением к движутся налево пропорционально к1 2, заполняя в пределе зону 11 (см. [1, с. 107]). Разумеется, вырождение оператора ¡к в зонах 11 и 111 не регулярно. Тем не менее при фиксированном Х > Д удается для функции распределения моментной задачи

п (Х)= £1 (3)

0<Х^ <Х

найти асимптотически точную оценку. Введем следующее обозначение

о(5, Х) = | V Х-ф1 (|(4)

Пусть Л - число удовлетворяющее при фиксированном Х неравенствам Л > Х, Л > 7, где

1 -о2

4 = ^РТ^ > 5 < 5 2 , (5)

В (5 )

В(5) - расстояние от меридиана до оси вращения. Рассмотрим по аналогии с [6] функцию а(5,Х) - характеристический определитель оператора ¡0, рассматриваемого на отрезке [, 5], 51 < 5 < 52.

Обозначим через р0 (Л) - число 5- нулей а( Х) на интервале (, 52), а через р1 (Х, Л), число ¿-нулей а(52,1) на интервале (Х, Л).

Для оператора ¡к определим числа т и V согласно табл.1. В первом столбце и в первой строке таблицы указан тип граничных условий в точках ^ = 51 и ^ = s2 соответственно. В каждой клетке таблицы слева дано число т , справа V. Таблица симметрична относительно граничных условий в точках

5 = 5 и 5 = 52

Таблица 1 .

(1,1,1) (1,1,2) (1,2,1) (2,1,1) (2,2,1) (2,1,2) (1,2,2) (2,2,2)

(1,1,1) -0.5, 0

(1,1,2) -0.25,0 0, 0

(1,2,1) 0.25, 0 0.5, 0 0, 1

(2,1,1) -0.5, 0 -0.25,0 0.25, 0 -0.5, 0

(2,2,1) 0.25, 0 0.5, 0 0, 1 0.25, 0 0, 2

(2,1,2) -0.25,0 0, 0 0.5, 0 -0.25,0 0.5, 1 0, 1

(1,2,2) 0.5, 0 -0.25,0 0.25, 0 0.75, 0 0.25, 1 -0.25,0 -0.5, 2

(2,2,2) 0.5, 0 -0.25,0 0.25, 0 0.5, 1 -0.75,2 -0.25,0 -0.5, 2 0.5, 3

Ниже через [х] и {х} обозначены целая и дробная части х. Справедлива следующая Теорема 1. Пусть в1 + £<Я<Я0,£> 0

Пусть

|а( 2, Я}>£ (6)

Пусть Я и [ всюду таковы, что £<{т + 5(2,Я)/[п}< 1 — £. Тогда существует такое [0, что при всех 0 < [ < [0 справедлива формула

п

(Я)

т +—— 2, Я)

+ Р0 (Л)-Р1 (Я, Л)+у (7)

Здесь [4 = к2 /12 - малый параметр, а фазовая поправка т и целое число У,0 <У < 3, определяются путем решения безмоментной задачи на малом отрезке (см. табл. 1 ).

Впервые асимптотическая формула (7) была получена в [5] при условиях жесткого защемления торцов. Впоследствии формула была обобщена на некоторые типы граничных условий [6] .

Как и в цитированных работах, при выводе формулы для пк (Я) существенно используется осцилляционная теорема, основанная на том, что при сужении отрезка [, 5 2 ] с.з. задачи, расположенные праве некоторого фиксированного у (см. (5)) возрастают. Последнее утверждение при произвольных граничных условиях требует специального рассмотрения.

Литература.

1. Гольденвейзер А.Л., Лидский В.Б., Товстик П.Е. Свободные колебания тонких упругих оболочек. М.; Наука, 1979.

2. Асланян А.Г., Лидский В.Б. Асимптотические формулы для частот нижней серии в теории оболочек вращения.//Изв. АН Фрм. ССз. Сер. Матем. 1971, Т. 6, №2-3, стр. 113-130.

3. Асланян А.Г. Связь моментной задачи с безмоментной в теории колебаний тонких упругих оболочек. //Изв. АН СССР. Механика твердого тела. 1 977, №5, стр. 118-124.

4. Вишик М.И., Люстерник Л.А. Регулярное вырождение и пограничный слой для линейных дифференциальных уравнений с малым параметром. //Успехи матем. Наук. 1957, т. 12, вып. 5(77), стр. 3-122.

5. Асланян А.Г., Лидский В.Б. Формула для числа частот осесимметричных колебаний оболочки вращения. //Дифференц. Ур-ия. 1977, т. 13, № 8, стр. 1355-1365.

6. Асланян А. Г. Формулы для числа частот осесимметричных колебаний оболочки вращения при различных граничных условиях. //Функц. анализ и его прлож. 1978, т. 12, вып. 3, стр. 61-63.

7. Асланян А.А., Асланян А.Г., Лидский В.Б. Асимптотические формулы для частот осесимметричных колебаний оболочки вращения. //Журнал вычислительной математики и математической физики. 1998, т. 38, №2, стр. 298-309

i Надоели баннеры? Вы всегда можете отключить рекламу.