Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
19
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§2.11. Динамические задачи линейной теории упругости

231

N

оо

 

U = U + w =

+ J2{ank + ank)u (n)) cosu kt -

 

k = 1

n = 1

 

oo

oo

 

- (wfc+ ^ 2 iKk + K k u (n)) sinLVkt) + u + w + ^ 2 (A'n cosLVnt -

T" sino;nt)u(n).

71=1

71— \

(2.11.123)

 

 

Это решение называется вынужденными колебаниями упругой среды под действием полигармонических массовых сил и сводит исходную задачу (2.11.1) к квазистатическим задачам (2.11.107), (2.11.109), (2.11.115) отно­ сительно w^//,_, и и и/^//, а также к задаче (2.11.77) на собственные значе­ ния оо\ и собственные функции U(n), после решения которых по формулам (2.11.119а), (2.11.122) вычисляют а^к и а ^ . Далее по (2.11.118) восстанавли­ вают поля и^" и и^" и по (2.11.116), (2.11.111), (2.11.1146) вычисляют а затем по (2.11.123) находят поле и(х,£).

Формула (2.11.122) имеет важную особенность: если заданная частота вынужденных колебаний одр окажется совпадающей с какой-либо собственной частотой ооп, то соответствующий коэффициент устремляется в бесконеч­ ность, что приводит к неограниченным значениям амплитудной составляющей и (2.11.118) вектора перемещений. Этот эффект называют резонансом, он играет большую роль при расчете конструкций и сооружений на динамиче­ ские нагрузки, поскольку действие сравнительно небольших по амплитуде внешних массовых сил может привести к возникновению неограниченных амплитуд колебаний конструкций и, как следствие, к их разрушению. Что­ бы предотвратить это явление на практике сначала вычисляют собственные частоты ооп тела (конструкции, сооружения), а затем сопоставляют их с ча­ стотами одр вынужденных колебаний; в случае близости значений некоторых из частот ооп и одр принимают меры по изменению, например, геометрии тела или допустимых внешних сил.

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

2.11.9. Установившиеся колебания упругих сред

Важным частным случаем вынужденных колебаний упругой среды явля­ ются установившиеся колебания, для которых:

1)частоты вынужденных колебаний одр существенно меньше первой (низ­ шей) частоты собственных колебаний:

оор < оо\ Vfc;

(2.11.124)

232 Глава 2. Упругие среды с малыми деформациями

2) функции начальных данных UQ, VQ нулевые,

тогда в задаче (2.11.1) можно пренебречь начальными данными, а в выраже­ нии (2.11.123) — соответствующим им полем и.

Поскольку при условии (2.11.124) из (2.11.122) имеем 2 ^ < то амплитудами, соответствующими этим коэффициентам также можно прене­ бречь, и решение (2.11.123) задачи (2.11.1) принимает вид полигармониче-

ской функции:

N

u = 5 > '

cos ujpt — ик sin cDyt) + u,

(2.11.125)

k= 1

 

 

где

 

(2.11.126)

ur

= ur + w r

 

Эти поля можно вычислить непосредственно, решая задачу, являющуюся суперпозицией задач (2.11.107) и (2.11.108), (2.11.109):

f °ршк1# + V . (4С ..

V ® и£") + pf£" = 0,

\ п • • V

(X) и^"|Ест =

Ц",

(2.11.127)

и'-"|Еи = <•";

ГV (4С

V ® й) +

°pf =

0,

} п • •

V 0 uL

=

tne,

(2.11.128)

uL = йе.

2.11.10. Колебания оболочек

Одними из основных объектов приложения рассмотренной выше теории свободных и вынужденных колебаний являются оболочки и пластины. Ди­ намические уравнения теории оболочек могут быть получены из уравнений квазистатического равновесия (2.10.41), выведенных в и. 2.10.6. Для этого плотность массовых сил f a , входящую в выражения (2.10.29) для массовых усилий Fea и массовых моментов Меа, следует заменить на обобщенные

массовые силы f a , включающие силы инерции:

 

-

д2и

(2.11.129)

f a = f a

g2a , а = 1,2,3,

где иа — компоненты вектора перемещений и в ортонормированном физиче­ ском базисе га срединной поверхности оболочки (см. и. 2.10.2).

Подставляя вместо иа выражения (2.10.21), соответствующие модели обо­ лочки Тимошенко, находим

7

f

д2и«

1 о.

7

,

d2w

(2.11.130)

/“ = /“

at’-

х at--' а

U2-

h

h

a ? '

 

§2.11. Динамические задачи линейной теории упругости

233

Если подставить (2.11.130) в (2.10.29), то получим следующие выражения для массовых усилий и моментов:

 

h/2

 

 

 

h/2

 

 

 

 

 

°ph3 d2la

F —

 

°pfa d X

3

= F ea- p h ^ X ,

м еа =

Pf° a d X

3

d X

3

= Mea

-1- е а

J

 

 

 

12 dt2 ’

 

 

 

ot

-h/2

 

 

 

 

 

 

-h/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a = 1,2;

 

 

 

 

 

 

 

 

 

 

Fe3 =

d 2W

 

 

 

 

 

(2.11.131)

 

 

 

 

Fe3- h~ W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя эти выражения в уравнения (2.10.41) вместо Fea и Меа, получаем

динамические уравнения теории оболочек:

La(T) + A xA2 (kaQa + Fea -

p h ^ f ) = 0,

 

La{ M ) - A xA2 {Qa - M eaF ^ ^ - ) =

0 , a =1,2;

(2.11.132)

(A2 QX),X + (AXQ2\ 2 - A xA 2 (k xTxx + k2T22

+ Ap - Fe3 + °p h ^ ~ ) = 0,

где дифференциальные операторы La(T) выражаются по формуле (2.10.41а). Определяющие соотношения (2.10.43) и кинематические соотношения (2.10.26), (2.10.27) остались без изменений, в результате получаем замкнутую систему динамических уравнений теории оболочек Тимошенко относительно пяти неизвестных функций U\, 7ь 72 и W, зависящих от X 1, X 2 и t. Эта система имеет второй порядок производных по координатам и времени и относится к типу гиперболических систем дифференциальных уравнений с

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

После добавления к этой системе граничных условий (2.10.50), а также начальных условий

t = 0: Ua = Uaо, 7с = 7с0, W = W0,

dUa/d t = Ua1, dla /d t = 7сь dW /dt = Wu а = 1,2, (2.11.133)

получаем постановку динамической задачи теории оболочек (2.11.132), (2.11.133), (2.10.26), (2.10.27), (2.10.49) модели Тимошенко, где Uaо, 7^0, Wo, Ua1, 7d и W\ — заданные функции координат Х а.

В задаче о свободных колебаниях оболочек граничные условия (2.10.50) следует принять нулевыми (условия свободного края (2.10.51), жесткой задел­ ки (2.10.52) и шарнирно-закрепленного края (2.10.53) являются примерами таких условий), а также считать отсутствующими все массовые усилия и моменты:

Fea = 0, Меа = 0, РеЪ= 0 , Ар = 0.

234 Глава 2. Упругие среды с малыми деформациями

Решение данной задачи в соответствии с общей теорией ищем в гармони­ ческом виде (2.10.72):

оо

n = Y , ( AnCosujnt - A i; Smujnt)n {n}(X a), Q = {Ua, la ,W } . (2.11.134) П=1

Подставляя (2.11.134) в кинематические соотношения (2.10.26), (2.10.27) и определяющие соотношения (2.10.43), находим, что усилия, моменты, тан­ генциальные деформации и искривления оболочки также имеют гармониче­ ский вид (2.11.134), т. е.

^ = \UQLI ТскМар, ^aph Qa}• (2.11.135)

После подстановки выражений (2.11.134), (2.11.135) в уравнения динами­ ческой теории оболочек (2.11.132) и граничные условия (2.10.50), получаем

следующую

задачу

на собственные значения

ип

и

собственные

функции

Ua(n), 7a(n),

 

 

в теории оболочек Тимошенко:

 

 

 

 

La

) Т ^-1Л-2 (kaQa(n) Т Р^^п^арп))

 

 

 

 

 

<

+

A

\ A 2

( Q a ( n

^ — ( p h /12)o;nyQ/(n))

0,

O L

1,2,

 

, (^2Ql(n)),l +

(^lQ2(n)),2 —^ 1^ 2(^1^ 11(71) + &2^22(rc) —

= 0,

{ Laa(n)

 

Caa^aa(n) + Capeap(n),

T\2 (n) — 2C66ei2(n).

 

Maa(n)

 

 

T ^ск/3^ск/3(п) ’ ^12(n)

 

2Z)00Xi2(n) >

 

Qa(n)

 

^Сб—CK,6—a&a3(n) >

^

 

 

 

 

 

&а,(3(п)

&а/з(иа(n),

H7^)),

Xap(n)

^apiS)a(n)> "^(n) >^/3(n))>

(E7ea, Щ

= 0,

(We,

Q£) = 0,

(7ea,

 

= 0,

a =1,2.

(2.11.136)

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

Подставляя выражения (2.11.134) в начальные условия (2.11.133), полу­ чаем следующие соотношения:

' оо

 

 

Е

К П (п)(Х а) = Ы Х а), По =

{C/Q0, 7a0, Wo},

< n=1oo

wnA"fi(n)(Xa) = П!(Х“),

(2.11.137)

-

E

fij = {C7al) 7al, ТЩ.

k

П=1

 

Система

собственных функций

(при каждом п имеется пять

функций: J7a(n), 7a(n), W(n)), согласно теореме Гильберта — Шмидта, является полной, и поэтому коэффициенты А" разложения функций По(Ха), Hi(Xa) по системе собственных функций всегда существуют — это коэффи­ циенты Фурье.

Аналогичным образом можно рассмотреть задачи о вынужденных и уста­ новившихся колебаниях оболочек.

236

Глава 2. Упругие среды с малыми деформациями

Выражая из (2.11.1406) функцию 7(п)ц и подставляя ее в уравнение (2.11.140а), продифференцированное по X 1, получаем следующую задачу для функции W(ny.

' D n W{n)M n+ L o l ( 2 ^ + Dp)w(n)M-

 

- “2nBp (

' - ^

2n)W(n) = 0,

(2.11.141a)

Щ (п)(0) = Win](l) = о,

w(n)Al(0) = win]M(i) =0.

(2.11.1416)

Решение задачи (2.11.141а, б) ищем в виде

 

 

Щп )

= w n sin

im X 1

n =

1,2,...

(2.11.142)

l

 

 

 

 

 

После подстановки (2.11.142) в (2.11.141a) получаем следующее алгебра­

ическое квадратное уравнение относительно

 

 

 

 

 

 

 

, ,2 ( JD

D pB p 2

|

\

Y

( D D + 2

^

I I

= D

/

п

\ 1

*

(2.11.143)

шп \ В Р ----- лгч— шп" П

+'

j

/

I

\

U '

гч

“ 11

V

7

У

 

 

O55

 

\

I

 

 

O55

J J

 

\

I

/

 

 

Решая это уравнение, находим две последовательности (ветви) собственных частот колебаний балки:

 

1

 

п =

1,2,...,

(2.11.144)

u i = — {Bn ± y B l - A A C n ),

где

 

 

 

 

 

Л— ВрРр

°p2h3

Вп = Вр + (Л!

D0 +

D\\BP

Сп = D n (Trn/l)4.

бъъ

12С55’

Съъ

 

(2.11.145)

 

 

 

 

 

Для большинства упругих сред параметры АСп/В 2 оказываются обычно малыми: АСп/В 2 <С 1, тогда для низшей ветви собственных значений (знак «—» в (2.11.144)) после линеаризации получаем следующее выражение:

<4 = 1вл п

D u ^ n /lf

(2.11.146)

Bp + (im/lf(Dp + DuBp/C55) ’

 

или после подстановки выражений (2.11.145) в (2.11.146) и отбрасывания членов порядка (п//)2, малых по сравнению с 1:

C \ \ h 2 ^7гп^4

а7 _

/7Г/ш \2

С \ \

(2.11.147)

4 = 12p(l+cen) V I

=

\ ~ Г )

12(755*

 

Для балок с конечной сдвиговой жесткостью отношение Сц/(12С55) ~ ~ 1 конечно, и параметры ап оказываются много меньше 1: a n 1, поскольку hjl < 1. В этом случае сдвиговыми свойствами балки можно пренебречь, то­ гда из формулы (2.11.147) получаем следующее выражение для собственных

частот:

,

 

шп = авНтгп/1)2,

ав = у С ц/(12р).

(2.11.148)

 

§2.11. Динамические задачи линейной теории упругости

 

237

Для балок с малой сдвиговой жесткостью

отношение

С ^/С ц

<С 1 и

параметр

ап может иметь значение порядка 1,

а

далее —

много

больше

1. Балку,

для которой выполняется условие ап

1, называют балкой с

очень малой сдвиговой жесткостью, для нее в формуле (2.11.147) можно считать, что 1 + ап рз ап, тогда для собственных частот получаем следующее выражение:

сип = ав ттп/1, ав = \ / С ^ /р .

(2.11.149)

Параметр ав представляет собой скорость поперечных волн в ортотропной среде (см. п. 2.11.2).

Колебания балки с очень малой сдвиговой жесткостью не зависят от ее изгибных упругих свойств (т. е. от Си) и определяются только модулем сдви­ га С55. Аналогичный эффект имеет место для струн (длинных тонких тел), у которых отсутствуют изгибные жесткости, а их колебания определяются продольным нагружением струны сг^. Формула для частоты колебаний струн

совпадает с (2.11.149), где ав ^ а в = \Ja\\/°p.

На практике часто вместо ип используют круговые частоты колебаний:

 

=

(2.11.150)

которые измеряют в герцах (1 Гц =

1 оборот/с).

 

Из (2.11.149) и (2.11.150) следует, что

 

Vn = v\n, п = 1, 2,...;

щ = Л‘ = °Л.

(2.11.151)

Первую круговую частоту щ называют основной частотой (основным тоном) колебаний, а остальные частоты vn (п > 1) — высшими частотами

(обертонами).

Звуковыми колебаниями называют колебания (свободные или вынуж­ денные) с частотами ип в диапазоне 17... 20000 Гц, которые способно вос­ принимать человеческое ухо. Колебания с частотами Vn < 17 Гц называют

инфразвуком, а свыше 20 кГц — ультразвуком.

В силу условий (2.11.73) и (2.11.83) ортонормированности собственных функций (называемых в теории оболочек собственными формами коле­ бании), множитель Wn в формуле (2.11.147) следует выбрать следующим образом: Wn = л/Z jl для всех п.

Из формулы (2.11.142) следует, что собственные формы колебаний балки имеют синусоидальный вид, причем число полусинусоид укладывается на

длине балки целое число раз, т. е.

 

Хп = 21/п, п — 1,2,...,

(2.11.152)

где Хп — длина волны собственной функции W ^ ( X l) (см. рис. 2.11.5).

238

Глава 2. Упругие среды с малыми деформациями

Если же

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

ки (2.11.138)

(при Ме1 = 0, Feз = 0), обусловленных только начальными

условиями (2.11.138г), то, согласно формуле (2.11.72), решение будет иметь следующий вид:

00

( An cos (jjnt — Af/ sin (jjnt) sin

X х

(2.11.153)

W (X l,t) = y 2 / l J 2

— .

n = 1

 

 

Введя обозначения для угла (рп:

 

 

coscpn = А'п/\А п\, sin ipn = А'п/\Ап\, \Ап \2 = А ' 2

+ А"2,

(2.11.153а)

формулу (2.11.153) можно представить в виде

 

 

W { x \t) = У Щ

cos(cont + срп) sin 2тгХ‘

(2.11.154)

^ \ А™

 

 

п = 1

Модули п| называют амплитудами колебаний, а (рп — углами сдвига фазы.

Амплитуды и фазы колебаний зависят от А '£ , которые определяются по

(2.11.85):

 

 

 

 

 

А ' =

ИЦАГЩ т ^ - d X \

А "

= -

V 2

W\ (X 1) sin П Х [ d X [

 

 

 

UJnX i

 

 

 

 

 

(2.11.155)

Например,

если W\ = 0, WQ ф 0,

то

А" =

0

и (рп = 0, т. е. фазы (рп в

упругих средах определяются только начальной скоростью тела (в неупругих средах ситуация иная, см. § 5.6).

Если WQ(X 1) — симметричная функция относительно точки X х = 1/2: WQ(X 1) = WQ(1— X х), т о и з (2.11.155) следует, что А\ Ф 0, а А'п = 0 (п > 1), т. е. колебания происходят только с одной основной частотой щ.

2.11.12. Продольные колебания балки

А. Свободные колебания балки

В качестве еще одного часто встречающегося в приложениях примера

рассмотрим

задачу

о продольных свободных колебаниях балки, оба торца

(X 1 = 0,1)

которой

жестко защемлены (рис. 2.11.5), а в начальный момент

времени t = 0 действует некоторое продольное перемещение

U\ = U\Q(X x)

(скорость dU \/dt = U \\(X x)), удовлетворяющее условиям совместности:

 

 

и ао(0) = и а0(1) = 0, се = 0,1.

(2.11.156)

В этом случае задача (2.11.132), (2.11,133), (2.10.43), (2.10.26), (2.10.52) о свободных колебаниях, подобно квазистатической задаче о растяжении балки