
книги / Оболочки и пластины
..pdfк замкнутой, ее колебания удлинений переходят в аналогичные колеба ния для замкнутой оболочки.
Как и в случае цилиндрической оболочки, можно показать, что колебания не мо гут совершенно не сопровождаться деформациями удлинений и поправка, необходимая для удовлетворения условиям на краях, значительнее той, которая нужна для того,
чтобы удовлетворить уравнениям |
движения. |
|
Отсюда можно заключить, что на свободном краю удлинения сравнимы по вели |
||
чине с деформациями изгиба и |
практически эти |
удлинения ограничиваются только |
узкой полосой у краев. |
|
характера колебаний, возникающие |
Если представить себе постепенные изменения |
с возрастанием кривизны, начиная от -плоокой пластины и кончая замкнутой сфериче ской оболочкой, то класс колебаний, протекающих практически без удлинений, полно стью исчезнет. Основание для этого нужно было бы искать в быстром росте частоты всех колебаний, принадлежащих этому классу, при значительном уменьшении отвер стия в оболочке.
§ 5. ПОЛОГАЯ ОБОЛОЧКА
Дифференциальное уравнение собственных колебаний пологих обо лочек и слегка искривленных пластин в наиболее общем случае можно представить в виде [7]
(4,5,1)
Здесь
a F = F ( x , |
у , / ) — функция |
напряжений; у — объемный вес материала |
оболочки; |
С — коэффициент |
упругого основания; g — ускорение силы |
тяжести.
Дифференциальное уравнение (4,5,1) .с соответствующими гранич ными условиями, присоединенными к нему, полностью описывают коле бания пологих оболочек как обыкновенных (с = 0), так и лежащих на упругом основании, находящихся под действием заданной внешней по
верхностной нагрузки Z = Z ( x y у , t ) и испытывающих |
заданные внут |
ренние осевые (нормальные и сдвигающие) силы № \у N%t |
S°. |
Эти силы могут быть как статические (не изменяющиеся во време |
ни), так и динамические (изменяющиеся во времени по заданному, на пример, периодическому закону).
Рассмотрим задачу о собственных колебаниях пологой оболочки, имеющей в плане форму прямоугольника со сторонами а и b и опираю щейся на контур шарнирно. Предполагая в целях общности решения, что оболочка испытывает по двум взаимно перпендикулярным направ
лениям заданные сжимающие силы N°i = const, N% = const, будем иметь дифференциальное уравнение
Введем дифференциальный оператор L :
L “ 0 ,1 |
+ 4 D “ |
а л Г |
+ 2 ( 0 ,! + 2£>*') |
+ |
||
|
+ 4D " - £ F + d * £ ; ' |
|
<4 Д 6 > |
|||
Пользуясь (4,6,6), перепишем уравнение (4,6,1): |
|
|
||||
|
— |
+ - ^ - L w |
= 0 . |
|
-(4,6,7) |
|
|
dt2 |
|
yh |
|
|
4 |
Представим решение этого уравнения в виде произведения |
||||||
|
w = (A cos p t + |
В sin pi) W (x, |
у). |
$ , 6,8) |
||
Здесь |
|
|
|
|
|
|
есть частота собственного колебания пластины. |
|
|
||||
Внося (4,6,8) |
в (4,6,7), получим для |
W уравнение |
|
|||
|
L W |
- p 2 ^ - W |
= 0. |
|
(4,6,9) |
|
|
|
|
8 |
|
|
|
Найдем решение этого уравнения, |
удовлетворяющее |
граничным усло |
||||
виям. Граничные |
условия для |
функции |
W (х, у ) |
в случае заделанного, |
опертого и свободного края не отличаются от известных условий для прогиба w ( x , y ) при равновесии пластины.
Представим решение уравнения (4, 6, 9) в виде ряда Фурье, содер жащего произвольные постоянные.
Выполнив требования, что функция W (х, у ) удовлетворяет гранич ным условиям и является решением уравнения (4,6,9), получим си стему однородных уравнений для неизвестных произвольных постоян ных; эта система имеет решения, отличные от нуля, если ее определитель А ( р) равен нулю. Отсюда получаем уравнение частот
Д(Р) = 0. |
(4,6,10) |
Уравнение (4,6,10) имеет бесконечное множество решений, |
кото |
рые составляют спектр частот данной пластины. Частоты зависят от
двух параметров т и п ( т = 1 , 2, 3, ..., |
п = |
1 , 2, 3, |
...). Наименьшая ча |
стота называется частотой основного тона, |
остальные — частотами выс |
||
ших порядков, или частотами обертонов. |
функция |
W m n { x , y ) , которая |
|
Каждой частоте р тп соответствует |
определится с точностью до постоянного множителя. Функции W m n ( x , y ) , называемые собственными функциями, определяют формы изогнутой поверхности, соответствующие частотам р тп.
Определение прогиба w ( x , y , t ) в любой точке и в любой момент времени осуществляют так: разлагают заданные начальный прогиб w 0 и начальную скорость v 0 в ряды по собственным функциям
Частота основного |
тона |
равна |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
ла т г ~ т .. т а * |
|
T v* |
|
|
|
|
<4’«'2о> |
||||||||
|
|
p" |
~ |
^ |
V |
|
- ^ |
V |
|
|
|
|
|
|
|
|
||
В случае |
квадратной |
пластины |
со стороной |
а частоты |
будут соответственно |
|||||||||||||
|
|
Рт п = - J |
] |
/ |
у Г |
^ Dim* + 2D3m W + |
D&*, |
|
(4,6,21) |
|||||||||
|
|
|
Pii = |
-^ - |
j |
/ |
" |
V D l + 2D3 + Di . |
|
|
(4,6,22) |
|||||||
Для прямоугольной изотропной пластины частоты будут соответственно |
|
|||||||||||||||||
|
|
|
Р т п = — |
1 / |
- з г Ц — |
j * п * | , |
|
|
(4,6,23) |
|||||||||
|
|
|
Ри - S /f [OfM |
|
|
(4,6,24) |
||||||||||||
|
|
|
|
|
|
|||||||||||||
В случае изотропной квадратной пластины со стороной а частоты будут соответ |
||||||||||||||||||
ственно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рт п = |
j t |
л / ~ |
g£_ |
(т2 + п * ). |
|
|
|
(4,6,25) |
||||||
|
|
|
|
|
|
|
а2 |
У |
уh |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
2я2 |
, |
f~ g D |
|
|
|
|
|
(4,6,26) |
||
|
|
|
|
|
|
P il~ |
|
а2 |
V |
|
yh |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Форму колебаний, соответствующую данным частотам рЖп, определяют собствен |
||||||||||||||||||
ные функции |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V /т п — sin |
тпх |
|
ппу |
|
|
|
|
(4,6,27) |
||||||
|
|
|
|
-------sin |
Ь |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
а |
|
|
|
|
|
|
|
||
Полный прогиб пластины |
представится |
|
как |
результат |
наложения бесконечного |
|||||||||||||
ряда прогибов |
вида |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ - |
|
|
|
, |
|
|
|
|
ч |
тпх |
. |
ппу |
, |
(4,6,28) |
||
|
и>тп — (Атп cos Pmnt -t- Втп sin pmnt) sin |
|
sin |
- |
||||||||||||||
а именно: |
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
OO |
00 |
|
|
|
|
|
|
Pmn . |
|
Д . |
mnx |
. |
nny |
|
|||
W = |
|
|
a m/i COS Pmnt + |
|
(4,6,29) |
|||||||||||||
|
|
------ |
Sin pmnt. Js i n ------- |
Sin— — . |
||||||||||||||
|
m=1££n= (1 |
|
|
|
|
|
pmn |
|
|
/ 0 |
. |
|
b |
|
||||
Здесь Ctm n, Pmn — коэффициенты разложения |
начального |
прогиба и начальной скоро |
||||||||||||||||
сти в ряды по собственным |
функциям; |
коэффициенты |
рядов определяются |
по фор |
||||||||||||||
мулам |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
а |
Ь |
|
|
тпх |
|
ппу |
, |
, |
|
|
|
||
|
|
|
|
г |
С |
|
|
|
|
|
|
|||||||
|
|
a m n ~ ~ a b |
V |
) |
^ o s i n — |
|
sin — ^ — dxdy, |
|
|
|
||||||||
|
|
|
|
|
|
о |
о |
|
|
|
|
|
|
|
|
|
|
(4,6,30) |
|
|
|
|
|
|
а |
Ь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ппу |
J |
|
|
|
|
||
|
|
Л |
|
4 |
Г |
Г |
|
• |
« я * |
. |
J |
|
|
|
||||
|
|
Pmn = |
— Г- |
\ |
\ |
«О S14 |
------- |
Sin— ;— |
dxdy. |
|
|
|
||||||
|
|
|
|
ab |
J |
.J |
|
|
|
a |
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
о |
о |
|
|
|
|
|
|
|
|
|
|
|
Обозначим во |
втором из |
(4,6,35) подынтегральное |
выражение через |
|
V { W ) ; тогда |
|
|
|
|
|
v |
= - £ - ^ V { W ) d x d y . |
(4,6,36) |
|
Уравнение |
Лагранжа |
(4,6,34) |
после внесения в |
него выражений |
Т и V принимает вид |
|
|
|
|
|
|
q + |
p 2q = О, |
(4.6.37) |
где р — частота, определяемая формулой |
|
|||
|
pi = J g _ |
tt'V(W)dxdy |
(4.6.38) |
|
|
|
УЬ |
Я W2dxdy |
|
|
|
|
||
Выражение для прогиба получается в виде |
|
|||
|
w = |
(A cos р/ + В sin pt) W. |
(4.6.39) |
Точность определения частот в значительной мере зависит от удач ного выбора выражения для W В случае простых контуров подбор этих функций, как правилб, не представляет труда. В качестве первого приближения для функции W можно взять выражение, пропорциональ ное статическому прогибу пластины с теми же условиями закрепления края под действием равномерно распределенной нормальной нагрузки.
Выбор формы изогнутой поверхности является подменой рассматри ваемой пластины, близкой к ней, но более жесткой. При приближенном способе определения частот это дает завышение их значения.
Определим функцию W (х, у ) так, чтобы она давала минимум вы ражения (4,6,38). Для приближенного решения этой вариационной задачи воспользуемся методом Ритца.
Зададим выражение для W в виде
^ |
= |
аЕ |
» £.л |
|
тп |
п |
|
где Wmn — непрерывные функции, зависящие от двух параметров
т и п и удовлетворяющие всем |
граничным* условиям. Тогда (4,6,38) |
||
представится в виде дроби |
|
|
|
_2 _ |
М (Атп) |
(4,6,41) |
|
Р ~ |
N (Атп) ’ |
||
|
числитель и знаменатель которой — однородные квадратичные функции коэффициентов А тп. Найдем минимум этой функции. Составим частные производные д р 21дАтп и приравняем их нулю. Очевидно,
ар2 |
1 |
дм |
М |
dN |
1 |
Г дм |
2 dN \ |
dAmn |
N |
дАтп |
N* |
дАтп |
N |
V дАтп |
дАщп ) |
= 1 Г - 1 Г - ( М - Р 2^ - |
<4’6’42> |
N дАтп |
|
Следовательно, задача отыскания минимума дроби (4,6,41) равно сильна задаче отыскания минимума выражения
S = M — p2N , |
(4,6,43) |
или
5=И[^ F{W)~p2W]dx y■ |
(4,6>4) |
Внося сюда значение W по формуле (4,6,40), получаем 5 в виде одно родной квадратичной функции коэффициентов А тп. Отыскание мини мума функции S приводит к системе однородных уравнений первой степени относительно коэффициентов Атп\ из условия равенства нулю определителя этой системы получим уравнение частот:
Л(р) = 0. |
(4,6,45) |
Наименьшее решение этого уравнения, отличное от нуля, и будет при ближенным значением частоты основного тона.
Для иллюстрации метода Релея—Ритца рассмотрим пример. Опре делим частоту основного тона прямоугольной ортотропной пластины, у которой все четыре стороны жестко защемлены. Точное решение этой задачи нам неизвестно. Приближенное решение получим, полагая
W mn = ( > - ~ ~ У ( у 2 - -£-)* *"Ч/П, |
(4,6,46) |
где т , п — целые числа.
Функции W mn будут удовлетворять граничным условиям:
|
I |
о. |
|
|
|
|
|
|
|
|
при х — ± — |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
Wmn = |
= |
0, |
|
(4,6,47) |
при у = |
± |
Y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
w mn = |
= |
0. |
|
(4,6,48) |
|
|
|
|
|
|
°У |
|
|
|
|
В |
качестве |
первого |
приближения |
возьмем первый |
член |
суммы |
||||
(4,6,40): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W = А ('х 2 - -^-)2 ( у 2 - - j -у |
|
(4 ,6 ,4 9 ) |
|||
Внося это в (4 ,6 ,3 8 ) и производя интегрирование, получим |
|
|
||||||||
|
|
f t . = |
- Т |
Г - 1 |
+ |
|
° ’571 |
|
t 4 ’6 ' 50» |
|
Для квадратной |
изотропной пластины с жесткостью |
D получим |
||||||||
|
|
|
|
|
|
|
|
|
|
(4Д 51) |
Можно указать и иное выражение для W mn: |
|
|
||||||||
|
|
W |
= |^1— (— 1)'»cos 2тпх |
£l — (— l)"cos —jp -J, |
(4,6,52) |
|||||
|
|
w тп |
|
|
|
|
|
|
|
|