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

книги / Механика композиционных материалов

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
12.36 Mб
Скачать

где

(о — действительное число, называемое круговой

частотой

(см.

1 . 4 . 1 3 ) ) ;

k — волновой

вектор,

характеризующий

направле­

ние распространения волны (k — волновое число):

 

 

 

 

k = kv, А = |Л =

=

1,

(1.6),

длину волны

X и фазовую

скорость

с (скорость

перемещения.

плоскости равных фаз ф=£*дс A-(at в направлении единичного век-

—*

тора v, называемого фазовой нормалью) определим:

А — постоянный вектор, длина которого называется амплитудой:

А = |А |= у AtAt Разумеется, (1.5) будет решением уравне­

ний (1.4) не при всех A, k и со. Для физически линейной среды, оператор определяющих соотношений — линейный, и поэтому уравнения (1.4) будут удовлетворяться при любой амплитуде. В этом случае, чтобы (1.5) было нетривиальным решением урав­ нений (1.4), между волновым числом k и частотой ш должна, во­ обще говоря, существовать связь вида

F((0, 6 )= 0 .

(1.8);

Соотношение (1.8) называется дисперсионным (или частотным) уравнением. Если для всех частот со фазовая скорость распрост­ ранения волны с одинакова, то волна называется недиспергиру­ ющей. В противном случае (фазовая скорость волны зависит от частоты) говорят, что имеет место дисперсия (волна — дисперги­ рующая).

В некоторых случаях волновое число может быть комплекс­ ным:

k = k'—ik'\

(1.9)

(например, при введении в определяющие соотношения комплексных модулей для изучения динамики линейно вязкоупругих материалов; § 4 гл. 1). Тогда дисперсионное соотношение (1.8) записывается в виде двух уравнений с действительными пе­ ременными:

Z7' (со, k't k") =0, F "(со, k', k" ) =0.

(1.10)

В этом случае (1.5) можно записать в виде

 

и = Ае-к'хе«к'х~т<\

(1.11)

причем мы выбрали специальную систему координат, в которой

волновой вектор k направлен по координатной

оси

х\ = х.

Из

(1.11) видно, что величина kf описывает длину

волны

X и

фазо­

вую скорость с'\

 

 

 

а ветчина k" связана с явлением затухания ноли.

Облает^ частот, для которой уравнения (1.10) не имеют дей­ ствительного решения к', называется областью непропускания. Наличие таких областей называется явлением волнового фильт­

ра-

Рассмотрим теперь в качестве примера динамическую задачу -об упругом неоднородном стержне. Для него уравнение движения (1.1) будет иметь вид

(L13)

а граничные условия (1.2) и начальные данные (1.3) — вид

и |х=о — н°, В (X)

1 = S °,

(1.14)

 

ох

 

при /== 0 :гг= U, - ~ =

V.

(1.15)

at

 

 

Так же как и в статической задаче об

упругом стержне

(§ 1

гл. 4), введем быструю координату | и обозначим точкой произ­ водную по этой координате, а штрихом — производную по мед­ ленной координате х.

Ищем решение задачи (1Л3) — (1.15) в виде

 

 

и = v (х, t) + aNl (g) v’ (x, t) + a* [w a (g) v

(x, t)

+ Nt (g)

j +

+ “ 5 [ N, (g) V " + N, (g) ~

] +

«* [ A

(g) nIV +

 

U (0 ^ + ^(0 -g-]+ -

 

 

 

 

 

 

 

 

( и в )

<7=0

P

 

 

 

 

 

 

 

P

 

 

 

 

 

 

где локальные функции Nq{\) таковы, что

 

 

 

 

А ) = 1,

при

q < 20, q < 0,

0<О.

(1.17)

Поэтому суммирование

в (М б )

производится

от 0 = 0, так,

чтобы

выполнялось неравенство q

 

методики

осреднения,

полу­

Производя обычные

процедуры

чим рекурентную последовательность дифференциальных уравне­ ний для определения локальных функций (при q = — 1, 0, 1,...)

(ЕЛГв+2у + (EN,+Iy + £ & + . + EN„ - pN\ = hq

(1.18)

е

где hq — некоторые постоянные, удовлетворяющие условиям

< ENq+1+ ENq - pNq) = hq,

(1.19)

причем, как обычно,

( h = 0; q > 0; В Д = 0.

(1.2 0 )

Тогда уравнение движения при отсутствии «объемных сил» при­ обретает вид

hQvn+ hQ~dF~д2и + а;[v"' + ^ ] +

= S “ ‘' S ^ ' 5 ' t’(,+2_2S)(x’ <) = 0-

(L21)

о

Упражнение 1.1. Показать, что величины hq в (1.18), (1.19) совпадают с величина hq § 1 гл. 4, и поэтому

А0 = й= ll(-y)'' h„= о, <7 > 0 .

(1.22 )

о

Упражнение 1.2. Показать, что локальные функции Nq(l) в (1.16) (1.20) совпадают с локальными функциями Nq(%) § 1 гл. 4, и поэтому для них справедливы формулы (4.1.58) и (4.1.59).

Упражнение 1.3. Показать, что из (1.18) и (1.19) вытекает справедливость следующих соотношений:

к = — (р) =

Ро.

*1 = 0,

 

(1.23)

к (е)=к (?) - (к)>

к

ю-J-^L-fp (?) -

 

 

 

 

 

dr\

 

 

О

 

 

0

0

 

 

 

 

 

 

dr\

(1.24)

 

 

 

 

Я(Л)

 

 

 

 

 

 

Уравнение движения (1.21)

при а = 0 имеет вид

 

 

,

 

дЬ)

 

 

(1.25)

hv =

Р*“ 5 Г '

 

 

 

 

 

где использовано (1.22), (1.23). Сохранение членов при

а в пер­

вой степени благодаря

(1.22), (1.23) ничего не прибавляет к

уравнению (1.25).

 

 

Подставим, согласно

(1.5), выражение

 

 

v (х, t) =

(1.26)

в уравнение движения (1.21), получим

 

Ш2 + р0ш2 + а2(о2 (h2k2 -f Ла(й2) Ч- aHka2 (h3k2 + h3w9) —

— a4©2 (ht&v2+ h4k2(о2+ Л Х ) + . . . = 0 .

(1.27)

Уравнение (1.27) является дисперсионным уравнением вида (1.8). Из него видно, что при а = 0 уравнение движения (1.21), которое в этом случае имеет вид (1.25), дисперсии не описывает, ибо из (1.27) для этого случая имеем

т. е. фазовая скорость для всех частот одинакова и совпадает со скоростью стержневых волн в эквивалентной однородной среде. Из (1.27) также видно: чтобы описать дисперсию, необходимо в

уравнении удержать члены с а2.

задаче

(1.13) — (1.15). Для ее

Возвратимся

к динамической

решения воспользуемся представлением

 

 

 

i / =

£

a >

p.

 

(1.29)

 

 

 

о

 

 

 

Тогда получим

рекуррентную

последовательность

динамических

задач для упругого однородного стержня

(р = 0, 1,....)

(ьз°)

 

+ Х Р ~

Ро

 

 

 

%\z, = up,

hw'o|j, =

Sp,

(1.31)

 

при t = 0 : w. =

Up,

 

= Vp,

(1-32)

 

и

 

y

at

 

где входные данные каждой задачи

(1.30) — (1.32)

определяются

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

' >0:

/•=1 0

= р > о; s0=s»,

 

 

 

 

 

“ р

=

А к .

р >

 

«о = и°;

 

 

 

 

 

 

 

А

=

вр!<=<>•

р > 0;

(/0 =

с/,

 

(1.33)

 

 

 

 

 

 

 

 

 

р > о;

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A e

- I

3 E A

<E> ^ S

- < - i 2- 26)

 

(1.34)

 

 

 

 

 

 

Г=1

Р

 

 

 

 

 

 

§

2.

Динамическая

задача

теории

упругости

и вязкоупругости

(

Рассмотрим

динамическую

задачу

теории

упругости

(1.3.16),

<1.3.19),

(1.2.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[Cijki(x)Uk.t].i + Xt = pWt"

 

 

(2.1)

 

 

 

 

Щ Iz, =

u? ,

CiM .(x)Uk.iij|г, =

Si ,

 

(2.2)

 

 

 

 

 

при t =

0 :u t =

Ui,

tii = Vi.

 

(2.3)

Ищем решение задачи

(2.1) — (2.3)

в виде

 

 

 

 

 

 

 

4=0

2

 

 

 

© - 5 - ° / . * . - * ^ .

(2-4)

 

 

 

 

 

PI

 

 

 

 

 

 

 

 

где суммирование по р происходит от

р = 0

так,

чтобы

выполня­

лось

условие q— 2р> 0

(все величины,

имеющие

отрицательные

индексы, равны нулю).

 

 

 

 

 

 

 

 

 

Для определения локальных функций ,N(9)(P) требуется решить

рекуррентную

последовательность

дифференциальных уравнений

при

q= — 1, 0,

1,...;

р = 0, 1,...

 

 

 

 

 

 

 

 

^iimn^mlk^kqJr2_ 2^ \ i + (CIM TMV^SSi . V . 'V . V 1'

+

 

+ ^ * <7+2-2Pmri^m/t1l.).(.^_f.1_ 2p|n + ^ t'fe<7+2-2Pm^4-r-2p^m/ftlP!.fe(7_2p _

 

 

 

— пЛ^^^НР-1)

 

=

 

 

 

• ДМ9М0) = g .

(2.5)

 

 

 

V . ilk t ...k g+2- 2 f l

 

lk q + 2- 2$l.kq + l - 2 fi'

'4

l ]

'

Как и прежде, локальные функции должны удовлетворять усло­ виям

(NW<P>) =

0,

р >

0; [[N^XP>]] = 0.J

(2.6)

Постоянные величины

Лпл>1(.?>.^+2-2р

определяются

осреднением

по ячейке периодичности следующих выражений:

 

/|(<7)(Р)

=

Л/(<7+П(Р)

4-

ikq+ 2 — 2$lkq+l—2Pft*

'

‘Л</+2—2Pmri т/Л‘ -••*<7+1—2Р*п

+ C’ift(7-|-2_2p^+l-2p ^ тЙ Р). г?-2Р

 

(2.7)

:<7+2-2Р/

В частности, из последнего выражения следует:

*!?>"> = - <Р>

= {СШк + ClimnN M $ ).

(2.8)

В результате введеных соотношений однородные уравнения дви­ жения (при отсутствии объемных сил) примут вид

0,4

•V H -2 P d t2fi V l>k i- ••V f 1-2Э'/ = ® *

p

Для отыскания вектора v нужно воспользоваться нием

V, = J ] apw}p},

/?=О

так что выражение (2.4) с учетом (2.10) примет вид

 

я

а2р

{я- р)^

и{ =

 

р=0 р

■Р-2Р Л2р

..ftp_2p.

я

 

 

^ . 9 )

разложе­

(2. 10)

(2.11)

Тогда получим

рекуррентную последовательность

динамических

задач теории упругости для

однородного тела

(р = 0, 1, 2,...)

 

h (o m J p )

, у { р ) _ / _ \ д 2 гЛ р )

,

/ 9 1 9 v

 

hijkl Wk.ij

+

Ai

— \Р/ “ ГГ ад

(2.12)

 

 

 

 

 

at2

 

 

w\P} I,, = и}'\

tiW u & jn j U. = S<(P),

(2.13)

при / = 0 : wl9) =

u\f ) ,

=

V',(P),

(2.14)

где входные данные каждой задачи

(2.12) — (2.14)

определяются

из рекуррентных соотношений

 

 

 

 

 

р = о,

 

 

 

 

 

 

XIjW .

 

 

 

 

 

 

 

'=1 0

 

 

 

 

 

 

S?,

р = о,

 

 

 

 

 

 

s/p>=

hljlSl..пг+1_ 2р а/2р

Ш«.п1.^лг+1_ 2р

к> Р > 0 ;

Г=1

 

 

 

 

 

 

 

{р}

 

и°г

р =

0,

 

 

 

 

 

 

 

 

 

 

и} ' =

 

 

 

 

 

 

 

 

Qi(P)k ,

Р > 0;

 

 

(U „ p = 0,

1

<?/P^ l/=o,

p > 0;

 

V,.

p = 0

 

 

v / p) =

 

 

(2.15)

- ^ Q i P)\t=o, P >

0;

 

 

a2P

{p-o

^ > s - 1 3 S

■2p

а/2Э Ш/.*|..Л-2Э •

г=1 Э

 

 

 

Нетрудно видеть, что величины h(<7)(0) в (2.5), (2.7) совпадают с величинами h^>, определенными в статической задаче теории уп­ ругости в § 2 гл. 4. Величины h(0)(0) = h являются компонентами эффективного тензора модулей упругости.

По теории эффективного модуля решается динамическая зада­ ча теории упругости для однородной среды:

ktnflk.il +

Х/ =

Po-Jr^f.

Р о= (р> ,

(2.16)

V t

к =

«?,

h t i k f l k . i r i j к

=

S °,

(2.17)

при

< =

0 :» , = £/„ - ^ -

= Vi.

(2.18)

 

 

 

dt

 

 

 

После ее решения можно найти микронапряжения по теории ну­ левого приближения:

о Г = с Ш | ) «*.«(*.<).

(2.19)

гДе Cify — компоненты тензора модулей упругости нулевого приб­ лижения

+

(2.20)

а локальные функции ЛД1)(0) находятся из решения дифференци­ альных уравнений

 

( C

/ / « e >

+

CiW)l(. =

0

 

(2.21)

при выполнении

условий

(2.6). Из

(2.8) и

(2.20)

следует,

что

 

 

 

 

 

 

 

( 2 .22)

В предыдущем

параграфе было

установлено, что

анализ

харак­

тера распространения плоских гармонических волн в данной среде позволяет ответить на вопросы:

а)

затухают ли эти волны в исследуемой среде;

б)

имеет ли место дисперсия;

в)

наблюдается ли явление фильтра.

Для исследования однородных уравнений движения теории эф­ фективного модуля (2.16) подставим в них (1.5). Т^гда получим однородную систему алгебраических уравнений

=

<2-23)

Для того чтобы система (2.23) имела

нетривиальное решение,

ее определитель нужно приравнять нулю. Используя обозначения

(1.6), получим так называемое уравнение Кристоффеля:

 

det |- i -

сг6,„ |= 0.

(2.24)

Решая уравнение (2.24), найдем направление вектороВ перемеще­

ний и и соответствующие фазовые скорости упругих эолн в экви­ валентной однородной упругой среде, характеризующейся эффек­ тивным тензором модулей упругости h. Нетрудно видеть, что волновые числа являются действительными, групповая скорость не зависит от частоты, т. е. волны вида (1.5) являются для рас­ сматриваемой среды недиспергирующими и незатухающими, яв­ ление волнового фильтра отсутствует.

Для того чтобы исследовать дисперсию, необходим0 в уравне­ ниях движения (2.9) сохранить члены при достаточно высоких степенях а.

Для решения динамической задачи линейной теории вязкоуп­ ругости достаточно во всех предыдущих соотношениях этого па­

раграфа заменить величины

С, N(^ ) ,

на соответствующие

вязкоупругие операторы R,

N<?Kp), h ^(|1).

 

Для теории эффективного модуля получаем уравнения движе­

ния

 

 

 

 

 

hijkpk.ij +

Х £ =

р0-^ —- vi* (р) ^ Р о,

 

(2.25)

граничные условия

 

 

 

 

 

vi Iz, =

w? ;

hljklVk,ifij |г, = 5?

 

(2.26)

и начальные данные

(2.18). После решения этой

задачи

можно

найти микронапряжения по теории нулевого приближения

 

 

 

«Г.? = « .% (5)

 

(2.27)

где Ryu — компоненты

тензора-оператора с ядрами релаксации

нулевого приближения

RifkiilJ):

 

 

R\jkl (£, t) =

Rtjki (£, t) + R;jmn

%

(2.28)

причем локальные ядра релаксации N(1)(0)(£, t) находятся из ре­ шения уравнений

[* //« . (Г) Л С Д + R ,M ( I /)],/ - 0

(2.29)

при выполнении условий (2.6). Из (2.28) можно найти эффектив ный тензор ядер релаксации

А|/« (О - A $f0) (0 = <Д|/°Ъ /)>

(2.30)

Особенности динамического поведения вязкоупругих материа­ лов можно рассмотреть на примере однородного вязкоупругого стержня (теория эффективного модуля). Уравнение движения такого стержня имеет вид

d*v

д*и

(2.31)

дх8 =

Ро' ~дР'

К исследованию этого уравнения можно подойти следующим способом. Подставляя (1.26) в уравнение движения упругого стержня, получим

hk2 = р0ы2.

(2.32)

Заменяя в (2.32) модуль упругости h на комплексный модуль Е* (1.4.16), получим

* = j l l L = < » / р7 ~

(2.зз)

VE*(iш)

 

где J* — комплексная податливость (1.4.24). Представляя эту податливость в тригонометрическом виде (1.4.25) и используя выражение (1.9), получим

k’ = u l ^ j 7 T a > s - p А' = ш К рЛ Л "8 1 п -5-.

(2.34)

Поэтому из (1.12) следует выражение для длины волны Я, и фа­ зовой скорости

\ ------------- -----------

— . с = ------------

!-------—

(2.35)

0) / р о | / • | cos

/ р о

I J*I cos —

 

Таким образом, в вязкоупругом стержне, даже однородном, воз­ никает дисперсия (вязкая дисперсия), причем гармонические плоские волны затухают. Явление фильтра отсутствует.

Возникающая в упругих композитах дисперсия носит название геометрической дисперсии.

Упражнение 2.1. Показать, что для малой вязкости, т. е. ма­ лых углов <р (1.4.23), из (2.34) и (1.4.24) следует

ti =г=со Кро («Л»+ J5),

Л' 0) / р , (Л .+ У)

(2.36)

2

с

V?o{J-+Js) -J4JJ'12

§3. Волновой фильтр

Вслоистых композитах наблюдается явление непропускания некоторых коротких волн (или высоких частот) в направлении, ортогональном слоистости. Это явление названо исследователями эффектом волнового фильтра.

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

Пусть длина ячейки периодичности равна 1, «объемная» кон­

центрация одного из компонентов у» а компоненты являются од­ нородными и характеризуются модулями упругости и плотностя­ ми Еа, ра или их комбинациями

Тогда

для каждого компонента

может быть записано

волновое

уравнение

&иа

 

 

 

 

 

с2

 

а =

1,2,

(3.2)

 

дх2

dt*

 

а

 

 

 

а на

границе раздела

компонентов Г (£ = 0 и £ = у) должны быть

выполнены обычные условия идеального контакта

 

 

_ „ р

 

Р

диг

(3 .3 )

 

ULu2t

дх

- Е2——

 

 

 

 

дх

 

После подстановки в уравнение (3.2) выражения (1.26) получим дисперсионное уравнение

COS k = COS ( —

+

——

CD

< '- x)a sin

sin- С - ^ -

(3.4)

\ Cl

С2

}

2x

 

су

c2

 

Спектр недопустимых длин волн Л=2л//г для композита эпок­

сид — бор

 

 

 

 

 

 

 

 

£ . =

0,515.10» — ;

р, =

0,137

КГ-с2

 

 

 

 

 

см2

 

 

см4

 

(3.5)

 

 

 

 

 

 

КГ с2

 

£ 2 =

0,39 - Ю7

— ;

Ра =

0,29

 

 

см4

 

 

 

 

 

СМ2

 

 

 

 

при у = 1/2 приведен в табл. 3.1.

 

 

 

Т а б л и ц а

3.1

 

2,440— 1,048

0,109—0,103

 

 

 

 

0,815—0,524

0,099—0,093

 

 

 

0,458—0,349

0,091—0,087

 

 

 

0,319—0,262

0,081— 0,063

 

 

 

0,244—0,210

0,059—0,055

 

 

 

0,156—0,144

0,053— 0,011

 

 

 

0,138— 0,126

0,009—0,007

 

 

0,122— 0,113