книги / Механика композиционных материалов
..pdfгде |
(о — действительное число, называемое круговой |
частотой |
||||
(см. |
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)(Р) |
2Р |
= (С |
Л/(<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 ) |
|
UL— u2t |
дх |
- — Е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 |
|
|
|
|
|