книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§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, |
|||
\ п • 4С • • V |
(X) и^"|Ест = |
Ц", |
(2.11.127) |
|
и'-"|Еи = <•"; |
||||
ГV • (4С • |
• V ® й) + |
°pf = |
0, |
|
} п • 4С • • |
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 |
/ |
7Гп |
\ 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 7ГП Х [ 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) о свободных колебаниях, подобно квазистатической задаче о растяжении балки