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

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

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

§ 4.7. Колебания линейно-вязкоупругих сред

451

преобразования (4.7.11), то для стационарных колебаний получим, что тензор деформации снова имеет гармонический вид (4.7.9), а амплитуды ст*п^ и связаны соотношением

£ (п) =

4тт* (

• • <т ( п ) >

(4.7.31)

* Ы

где 4П *{шп) — тензор комплексных податливостей, имеющий вид

 

 

ап

а

 

 

 

4п * Ы = £ ПЩ ип)

ааа(3-' +

*—'

Паа(шп) Га ,

(4.7.32)

а,/3=1

 

 

а = т

-\-1

 

ЩЦЦп) = 1Щ - 11Щ,

(4.7.33)

4П *Ц ) = 4П' - i4n".

(4.7.34)

Спектральные комплексные податливости П*Дщп) аналогично форму­ лам (4.7.19) и (4.7.14) представляют собой предельные значения при t —>+оо:

К р Ы = п а/3 + Nal3( x ) e - ^ X dx,

(4.7.35)

о

аопределяют аналогично (4.7.25):

°

оо

 

Г

 

П Щ Д , шп) = и а/3 +

Ма/3(х)еШпХsin Д ж dx,

 

=

Х ар(х)еШпХsinш'пх dx.

(4.7.36)

О

В силу взаимной обратимости исходных интегральных определяющих со­ отношений (4.6.4) и (4.6.6), будут взаимнообратными и соотношения (4.7.17) и (4.7.31), поэтому тензоры 4П* и 4R* также являются взаимнообратными:

4П *

(4.7.37)

Подставляя (4.7.29) и (4.7.34) в (4.7.37) и разделяя действительные и мнимые части, можно записать это соотношение в эквивалентном виде:

4П' • • 4R ' + 4П" • 4R" = А, - 4П" • • 4R' + 4lT • 4R" = 0.

(4.7.38)

Для экспоненциальных ядер ползучести Nap{t), имеющих вид (4.2.107), интегралы (4.7.36) также легко вычисляются (см. упр. 4 к § 4.7), в резуль­ тате получаем следующие выражения для действительной и мнимой частей комплексных спектральных податливостей:

 

N

л(7)/i

Па/з — Пар + Е

ЛавО шпгав)

7

~Р\ (1 “ ^пЩ )2 + (я'пТрр)2

§ 4.7. Колебания линейно-вязкоупругих сред

453

Методом упругой аналогии можно вычислить собственные частоты и собственные функции (формы) колебаний вязкоупругих оболочек, используя для этого формулировку соответствующей задачи (2.11.136) на собственные значения для упругих оболочек.

4.7.4.Свободные колебания вязкоупругой балки

Вкачестве примера рассмотрим задачу о свободных поперечных колеба­ ниях ортотропной вязкоупругой балки, жестко защемленной по обоим торцам (см. рис. 2.11.5 из п. 2.11.12). Используя метод упругой аналогии, с помощью формулы (2.11.149) находим выражение для собственных частот колебаний вязкоупругой балки с очень малой сдвиговой жесткостью:

<4 = f ( т ) 2'

(4.7.44)

О

о

Здесь учтено, что для ортотропных сред Са/з = Сар, поэтому замена Са/з —►

—►

правомерна.

 

 

 

Разделяя в формуле (4.7.44) действительную и мнимую части, находим

 

, /2 _

, // 2 (ттп)2R 55(u)'n!

, и ” )

(т гп)2Щ 5 (и 'п , и 'Р

(4.7.45)

 

^п

^п

 

pi2

 

 

Pi2

 

 

Поскольку R'I'( и R '^ зависят от и ^ (см., например, (4.7.45)), то формулы (4.7.45) представляют собой систему двух нелинейных алгебраических урав­ нений относительно и'п и CJ", решая которую, находим uJn и CJ".

Наличие мнимой составляющей CJ", как отмечалось выше, приводит к затуханию свободных колебаний вязкоупругой балки, в результате для про­ гиба W (X l,t) балки с учетом (4.7.22) и (2.11.153) получаем следующее выражение:

оо

W ( x \ t ) = Re ( J 2 Aneiu}nt)Wn(X l) = П=1

ОО

= ^2(Л'п cosUj'nt - А'п sinш'У)е~ш'^Шп(Х 1), (4.7.46) П=1

где Wn(X l) — собственные формы колебаний вязкоупругой балки (2.11.142), являющиеся вещественными; коэффициенты И* = А'п + iИ" вычисляют с помощью системы (4.7.20), которая в силу ортогональности функций Wn является конечной и имеет решение

А'п = лД[1 W0{ X x) s m ^I- d X \

(4.7.47)

о

 

454

Глава 4. Вязкоупругие среды

 

 

I

 

 

n X l

 

 

W\ (X 1) sin

WQ(X 1) sin

dX

о

T ~

 

0

 

 

Формулу (4.7.46) можно представить также в виде, аналогичном (2.11.154), если ввести модули \Ап\ и фазы рп по формулам (2.11.153а):

^

I А„ I г.пя(о/ i -А- (п„V- 1

XXI

 

 

W (X l, t) = J j

\An\cos(Jnt + <Pn)e ^

sin ЖпХ

 

 

 

n = 1

I

 

 

 

 

iff

2

 

 

 

 

 

(4.7.48)

cosifn = A'n/\A n\, sm<pn = A "/\A n\, \А?п\ = А!% + А.

 

*

 

 

in

 

 

Если для упругой балки углы сдвига фазы рп были обусловлены только ненулевой начальной скоростью W\ балки и обращались в нуль при W\ = = 0, то для вязкоупругой балки даже при W\ = 0 ввиду наличия мнимой составляющей частоты CJ" углы сдвига фазы всегда отличны от нуля: рп ф 0.

Как отмечалось в и. 4.7.1, из-за наличия комплексной составляющей ча­ стот CJ" свободные колебания вязкоупругой балки всегда будут затухающими по амплитуде во времени. Эффект затухания свободных изгибных колебаний консольных балок хорошо знаком каждому из его практического опыта. Этот же опыт подсказывает, что все твердые среды будут обладать вязкоупругими свойствами, поскольку свободные колебания будут затухающими для балок, изготовленных из любых материалов. Вывод этот действительно справед­ лив — разные твердые материалы различаются только быстротой затухания колебаний: наиболее быстро они затухают у резин — материалов с ярко выраженными вязкоупругими свойствами, несколько дольше у пластмасс, также обладающих существенными вязкоупругими свойствами, и сравни­ тельно долго будут затухать свободные колебания металлических балок, поскольку металлы обладают достаточно слабо выраженными вязкоупругими свойствами.

4.7.5. В ы нуж денны е колебания

Рассмотрим задачу (4.6.12) о вынужденных колебаниях линейно-вязко- упругих тел под действием полигармонических нагрузок

N

 

 

п = £

+ n, а = {f, tne, ue},

(4.6.49)

k=1

 

 

где Щ = Qfk + i заданные комплексные амплитуды колебаний; одр > 0 — заданные частоты колебаний (вещественные числа); Cl — статически действу­ ющие нагрузки.

Общий ход решения задачи (4.6.12), (4.7.49) такой же, как и для анало­ гичной задачи (2.11.1), (2.11.96) о вынужденных колебаниях упругих сред

§ 4.7. Колебания линейно-вязкоупругих сред

455

(см. п. 2.11.8): ее решение представляют в виде суммы

и = и + W.

(4.7.50)

Здесь W (x,t) — решение квазистатической задачи вязкоупругости, анало­ гичной задаче (2.11.97), (2.11.98); u (x,t) — решение динамической задачи (4.6.12) с нулевыми граничными условиями. Ищем эти функции в виде

N N

u = ^ u £ e i° fct + u + u(a;,i)I k=1

W = y ^ W *kelQkt + W,

(4.7.51)

k=1

 

причем для

+ iW^ и W получаем следующие квазистатические

задачи:

 

 

V • (4К*Щ ) • • V (g) W*k) = 0,

(4.7.52)

п • (4R*(wfc) • • V eg) W£)ECT= t*nek,

W*|Eu= u eV k = l , . . . , N ;

 

(V • (4R • • V <g>W) = 0,

(4.7.53)

 

jn • (4R • • V 0 W )sCT= t ne,

 

W |Su= u e,

для = ufk + iuk — задачу об установившихся колебаниях, а для и — квазистатическую задачу, обе с нулевыми граничными условиями, но с ненулевыми массовыми силами:

Рш1и1 + v

а к + рЦ -

о. Ч - Ч +

 

< сг*к = *В,*(йк) - - Т к,

4

= (1/2)(V®2^ + V ® s r ) ,

(4.7.54)

kn • <rfe|E^ =

0, ufelsu

=

 

fv .(* R ..

V 0 fi) + pf = O

{ n - (4R ..

V <8>G)ECT = 0, a + = o.

При формулировке задач (4.7.52) и (4.7.53) использован метод ста­ ционарных колебаний, изложенный в и. 4.7.1, согласно которому отклик линейно-вязкоупругой среды на гармонические возмущения (т. е. оператор определяющих соотношений (4.7.11)—(4.7.16)) также ищется в гармоническом виде. Тензор комплексных модулей упругости 4R*(5+) имеет такой же вид (4.7.14), (4.7.15), но йоу — вещественные.

Для функций и получаем задачу о собственных колебаниях вида (4.7.1) со специальными начальными данными (2.11.111):

N

N

 

t = 0: u = u0 = u0- u - W - ^ 2 (й 'к+ Wfc),

u = VQ+ y^tDfc(Wfc+ Ufc).

к=1

k~l

(4.7.56)

 

Решение этой задачи имеет вид (4.7.3).

§ 4.7. Колебания линейно-вязкоупругих сред

457

Разделяя в (4.7.57а) действительные и мнимые части, получаем уравнения

для вычисления оо'

и ш":

 

 

, /2 _

2 _

2

2

 

Кп гУ л / . //\

о. / . // _ Кп тэН( / //л

(4.7.58)

 

 

 

 

Используя формулу (2.11.180), находим перемещение вязкоупругой балки

при продольных вынужденных колебаниях

 

 

 

 

 

 

 

 

 

 

 

Ui = Re{{U* + \¥*)еш + U)

 

 

 

(4.7.59)

 

 

 

 

 

 

 

 

 

 

 

ОО

тгпХ1

 

 

 

 

 

VF* =

 

 

 

17* = U*

 

К sm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П=1

 

 

 

 

 

 

ь.-. = ^

+ <

=

(-1 )"+|

( Ц

2р и и 2 ии2

 

 

 

 

 

R * M ~ co 2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция свободных колебаний U, согласно (2.11.177), имеет вид

 

t/ = V y i Y . A> " nt sin

 

 

 

 

= <

+ Ь4",

 

(4.7.60)

 

 

n=l

 

 

 

 

 

 

 

 

 

 

 

 

 

где коэффициенты

 

вычисляют

по

начальным

данным

по аналогии с

(2.11.178), но с учетом того, что fr*

и сип комплексные:

 

 

 

 

4 i

= \ / f

7 0 sin

 

 

= л ; - (ь; + Сп) ^

+ а д ',

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- «

Д

+ Л > ") = Л

^0 sin

 

 

dX 1= А" - Д бЩ '' + Ь"Д ),

 

 

Uo = Uю - U ' - W ',

щ = Un ~ й(Ц" + W").

 

(4.7.61)

Разделяя в (4.7.59) действительную и мнимую части, получаем

 

U\ =

|U* + W* | cos(<Dt + <£) +

 

УУ Д п| cos(uj!nt + фп) е ~ ^ sin 7ГП'^ ,

 

 

 

 

 

 

 

 

П=1

 

 

 

 

 

 

(4.7.62)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ylrp

 

 

 

 

 

 

 

 

 

1 7 " + Т Р "

 

U‘

 

 

 

 

°°

 

 

. 7 Г п Х

tg <Р= -~7

 

 

 

 

s + £ ( В Д

- №

 

 

 

 

 

' + Ю = X

-

sm ‘

1

 

 

U' + W'

 

 

 

 

 

 

 

 

 

п = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х ‘С/"

 

 

 

 

 

 

жпХ1

 

 

 

 

+ Ж" =

 

 

+Ё(ВД+№e ; sm -

^

 

 

 

 

 

 

 

 

 

П=1

 

 

 

 

 

 

I чЗ

 

 

=

 

 

рУ -

р'У'

 

 

,

 

 

 

 

р

 

 

 

 

 

 

 

/2, // 2

 

к п = 2(—l)n+1 (— )

|Д ||’

 

| Д П 1 ( Д 2 - Д

 

' 2 - й

; 2 ) 2

+ 4

 

 

 

 

V '

7

 

 

Д У

 

 

 

 

 

7ГП

458

 

 

Глава 4. Вязкоупругие среды

 

ь" =

 

p'q"

-P"q'

=

q" = - 2 J nJ ',

I | (w„ -

-

u r f + 4w„ Ц

-

 

 

 

 

 

//

= R \lcj/ncd" + R u i^n

~ u " 2)- (4.7.63)

 

 

 

n

Из этих формул следует, что коэффициенты Ь^", в отличие от упругого реше­ ния (см. п. 2.11.12), ни при каких положительных а/, иоп и од не обращаются в бесконечность, т. е. амплитуды колебаний в (4.7.62) остаются ограниченными и резонанс отсутствует. Однако в стержне, как отмечалось выше, может возникнуть квазирезонанс, если oj"/oufn <С 1 и й рз uJn.

Отметим также, что при U" = 0 в случае упругой балки угол сдвига фазы (угол задержки колебаний) р = 0, так как в этом случае Ь" = О, U" + W " = = 0. Это означает, что вынужденные колебания упругой балки с частотой возмущения О происходят в одной и той же фазе во всех сечениях балки с координатой X х.

Для вязкоупругой балки ф ф 0 даже при U" = 0; более того, ф является функцией координаты X х, т. е. вынужденные колебания вязкоупругой балки происходят в разных ее сечениях с различным сдвигом фазы ф (Хх).

4 .7 .7 . Уст ановивш иеся и квазист ат ические колебания

Как и для упругих сред (см. п. 2.11.9), часто рассматривают случай

установившихся колебаний вязкоупругих сред, для которых: 1) частоты вынужденных колебаний йоу существенно меньше первой (основной) частоты собственных колебаний: йоу < ио\ Vfc; 2) функции начальных данных UQ, VQ нулевые. Тогда начальными данными в задаче (4.6.12) можно пренебречь, а в выражениях (4.7.50), (4.7.51) — соответствующим им полем й (х ,£).

Оставшиеся два поля в (4.7.51) и и W (первое соответствует ненулевым граничным условиям, а второе — ненулевым массовым силам) в данном слу­ чае можно рассматривать совместно, в результате решение (4.7.50), (4.7.51) можно представить в виде

U = У ] 11*ке

^ + й.

(4.7.64)

к=1

 

 

Здесь

 

 

и* = й* + W£,

а = u + W .

(4.7.65)

Объединяя задачи (4.7.52) и (4.7.54), а также (4.7.53) и (4.7.55), получаем для комплексных амплитуд и средних значений й перемещений следующие задачи:

'°рй2кик +

V •

+ °рЦ = °>

 

 

< <т*к = 4Щйк)--е%,

e* =

(l/2)(Vcx)u* + V cx )u n ,

(4.7.66)

n /-Г* I

_ -f*

П* I

_

п**

 

KU а к1ест

Uk\zu ~

Ue’

 

§ 4.7. Колебания линейно-вязкоупругих сред

459

V • (4R • •

V

(X) й) + °р{ = О,

(4.7.67)

п • (4R • •

V

0 u )sCT= t ne,

 

Задача (4.7.66) — это задача об установившихся вынужденных колебаниях вязкоупругих сред (где Ok заданы и вещественны), а (4.7.67) — квазистатическая задача вязкоупругости.

Если же все частоты вынужденных колебаний Ok много меньше основной

частоты собственных колебаний:

 

 

сгк < сщ

Vfc,

(4.7.68)

то инерционным членом в задаче (4.7.66) можно пренебречь, в результате приходим к постановке задачи о квазистатических колебаниях вязкоупру­ гих тел:

Отметим, что эта же задача может быть непосредственно получена из постановки (4.6.13) квазистатической задачи теории вязкоупругости, если подставить в нее решение (4.7.64) и рассмотреть стационарный режим вязкоупругих колебаний.

Пример решения квазистатической задачи (4.7.69) будет рассмотрен в и. 4.7.9.

4 .7 .8 . Связанны е задачи т ермовязкоупругост и

при уст ановивш ихся колебаниях

Как отмечалось в пп. 4.4.1 и 4.4.2, вследствие неидеальности вязкоупру­ гих сред, функция диссипации отлична от тождественного нуля, что может приводить к саморазогреву тела даже без притока тепла к нему извне. Эффект саморазогрева наиболее ярко проявляется при колебаниях вязкоупру­ гих тел. Поскольку процесс деформирования тел в этом случае становится неизотермическим, то необходимо рассматривать постановку (4.6.9) связан­ ной динамической (или квазистатической (4.6.11) при выполнении условия (4.7.68)) задачи теории термовязкоупругости, предполагая, что внешние на­ грузки П = {if, t fne,u e} в них изменяются по гармоническому закону (4.7.49). Дальнейшее преобразование постановок (4.6.9) и (4.6.11) в задачах типа (4.7.66) и (4.7.69) при установившихся колебаниях осуществляется аналогич­ но рассмотренному выше в и. 4.7.7 изотермическому случаю, но при этом принимаются дополнительные допущения:

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

uok = ки , и = и 1, к = 1,2, ... ;

(4.7.70)