книги / Механика композиционных материалов
..pdf
|
|
о>1 |
tON |
|
|
|
|
|
|
|
|
|
|
|
|
j |
J |
Q(a>lt |
,M [ /(« 1 * |
, ЮлО — |
|
|
|
|
|||
M |
|
w; |
co^ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
— J ] |
|
|
(“ l> |
. . . , <Dw)] фрK , |
, 0)*) |
|
|
, d(0Ny |
|
|
|||
<7=1 |
|
|
|
|
p = |
1,2, |
, Af. |
|
|
|
|
|
(3.9) |
|
|
|
|
|
|
|
|
|
|
||||
После^решения этой системы можно заменить в (3.5) |
о)? на опера |
||||||||||||
торы (о<7, |
1 |
, |
7V и решение задачи |
(3.3), |
(3.4) |
записать |
в виде |
||||||
|
|
|
_ |
|
м |
^ |
|
|
|
|
|
|
|
|
|
|
f< (*. 0 |
= 2 |
С,% (“ 1 ..........М |
и?, |
|
|
|
(3.10) |
|||
|
|
|
|
|
<7=1 |
|
|
|
|
|
|
|
|
а так как ядро ф<;(£—т) |
каждого канонического |
оператора |
уя мо |
||||||||||
жет быть найдено экспериментально, то |
|
|
|
|
|
|
|||||||
|
|
|
Щ(X, О = 2 |
СчJ % V ~ т) du° W - |
|
|
|
|
(З .п ) |
||||
|
|
|
|
|
<7=1 |
0 |
|
|
|
|
|
|
|
Описанный метод представляет собой аналитический |
метод |
ап |
|||||||||||
проксимации |
в теории |
вязкоупругих |
композитов. |
Если |
задачу |
||||||||
теории упругости, соответствующую задаче (3.3), (3.4), |
аналити |
||||||||||||
чески решить |
не удается, можно воспользоваться |
методом |
чис |
ленной реализации упругого решения, который представляет со бой численный метод аппроксимации в теории вязкоупругих ком
позитов. |
|
|
|
|
|
|
|
Предположим, |
что |
|||
Его существо заключается в следующем. |
||||||||||||
заданные |
на |
границе перемещения |
(3.4) |
можно представить в |
||||||||
виде |
|
|
|
|
|
|
|
|
|
|
|
|
«»(х, |
0 = |
|
|
а = |
1 ,2 ,3 ; |
/ = |
1..........L. |
(3.12) |
||||
|
|
!=i |
|
|
|
|
|
|
|
|
|
|
Положим |
теперь Va (t) = 0, / = |
2 , 3 , . . . , L; |
a Va (0 = 1 |
и |
|
Ре" |
||||||
шим соответствующую |
задачу |
при |
различных |
наборах |
|
чисел |
||||||
4 V), 0)2V)' ••• ’ ©Sr» У = |
1, |
.Г . |
Тогда, |
воспользовавшись |
ап |
|||||||
проксимацией для каждой компоненты а вектора |
у, |
|
|
|
||||||||
|
|
м |
С&<р, (а,. |
, cow), / = |
1,2. |
, L, |
|
|
(3.13) |
|||
|
VW = £ |
|
|
|||||||||
|
|
<7=1 |
|
|
|
|
|
|
|
|
|
|
в каждой точке тела получаем |
систему |
алгебраических |
|
уравне |
||||||||
ний |
м |
|
|
|
|
|
|
|
|
|
|
|
VWW= |
С® ф, (о/Д |
. (0%'); |
I = |
1, |
, L; у = 1. |
. Г |
|
(3.14) |
||||
£ |
|
для определения постоянных Cqa и, |
возможно, |
|
некоторых |
|||||||||||
«плавающих» |
констант |
операторов |
y q |
(поэтому |
Г > М ). |
После |
||||||||
определения |
всех этих |
констант полагаем |
в |
(3.12) |
VaW(t) = \, |
a |
||||||||
|
|
|
|
|
|
|
остальные |
|
Va{l\ / = 1, |
3, |
||||
|
|
|
|
|
|
|
4, |
|
L равными |
нулю |
и |
|||
|
|
|
|
|
|
— |
снова |
решаем |
систему |
|||||
|
|
|
|
|
|
алгебраических |
|
уравне- |
||||||
|
|
|
|
|
|
—1- |
нцй (3.14). Заметим, что |
|||||||
|
|
|
|
|
|
-+-Р |
параметры |
со/М; |
/= 1 ,..., |
|||||
|
|
|
|
|
|
—*- |
N; |
у= 1, |
..., |
Г |
не обяза- |
|||
|
|
|
|
|
|
-► |
тельно в |
соответствии |
с |
|||||
|
|
|
|
|
|
— |
их физическим |
|
смыслом |
|||||
|
|
|
|
|
|
|
выбирать |
между |
нулем |
и |
||||
|
Рис. |
65. |
|
|
|
единицей. Более того, как |
||||||||
|
|
|
|
показывает |
|
численный |
||||||||
мы «раздвигаем область изменения этих |
эксперимент, чем |
больше |
||||||||||||
параметров, тем точнее |
||||||||||||||
получаются результаты |
счета. |
После исчерпания |
всех |
VJlK / = |
||||||||||
= 1,2,..., L получим решение задачи теории |
вязкоупругости (3.3), |
|||||||||||||
(3.4) в квадратурах |
(по времени): |
|
|
|
|
|
|
|
|
|
||||
|
|
|
М |
L |
t |
|
|
|
|
|
|
|
|
|
|
Иа (X , 0 |
= |
£ |
£ |
C fa j |
Ф, (/ - Т) dV® (Т). |
|
|
(3.15) |
|||||
|
|
|
<7=1 |
*=1 |
О |
|
|
|
|
|
|
|
|
|
В качестве примеров рассмотрим две задачи о простом слоис том двухкомпонентном композите. В первой из них рассматрива ется бесконечная ортотропная пластинка с круговым отверстием. На бесконечности в направлении оси х действует растягивающая сила Р (рис. 65).
Из решения соответствующей упругой задачи следует, что мак
симальное напряжение |
а0 |
возникает |
в точках, |
показанных на |
|||
рис. 65 кружочком, причем* оно равно |
|
|
|
||||
|
((Те)max= (l+ * i + X2)P, |
(3.16) |
|||||
где xi и х2 — корни |
бигармонического уравнения |
|
|||||
__j _ |
1 |
2vjfи \ |
л, |
, |
1 |
(3.17) |
|
—— |
(—-----* 0 * - ] |
х* + |
- L - = о, |
||||
Ех |
\ |
Gxy |
Ех |
) |
|
Еу |
|
причем |
|
|
|
|
|
|
|
*4 + Х2 |
|
/ |
2 ( / t |
- v- ) + - |
(3.18) |
||
|
|
где Ех, Еу, GXy, vxy — модули упругости ортотропной пластинки. Компоненты тензора эффективных ядер ползучести известны. Они
выражаются через операторы 1, |
я, |
где |
Кч (2 -(- М]) |
||
Ер -- |
1+ Р© |
1 +9Д4 |
1 - Y |
М = |
|
|
■ i |
|
|
9/CIGL)! |
(3.1 9)
Подставляя в (3.18) вместо модулей упругости соответствующие операторы, получим
F (со) = |
+ |
>с2 = 2 11 — m (1 — v) ^ |
1 + |
т [2 + |
(1 — Ш ) со]2 |
|
+ |
|||||||
|
|
|
|
|
|
|
|
|
|
9/Иш(2 + ю) |
|
|
|
|
|
+ |
т |
2 + (1 —9М) а + 27М2 (1 + у) со2 |
] 1/2 |
|
(3.20) |
||||||||
|
3М |
|
|
|
со (2 + |
ш) |
|
|
|
|
||||
где |
|
т==у( 1—у), со= |
0)2, |
v = vj. |
|
|
(3.21) |
|||||||
|
|
|
|
|||||||||||
Запишем аппроксимацию функции F {Q ) |
в |
виде |
|
|
|
|
||||||||
|
|
|
F(со) = Л 1+ Л2я + Л 3^ Р. |
|
|
(3.22) |
||||||||
Чтобы иметь возможность интегрировать в (3.7) от 0 до |
1, |
вве |
||||||||||||
дем весовую функцию й = со. Тогда |
|
|
|
|
|
|
|
|
||||||
|
|
|
со/7(со) = В\-\-В2( |
й |
- |
\ |
- |
(3.23) |
||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В, = |
Аг + - ^ - , |
Вг = А г, Ва |
|
|
|
|
(3.24) |
|||||
Для нахождения величин Bi |
(i= 1, 2, 3) |
имеем |
систему |
алгебраи |
||||||||||
ческих уравнений |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
LiiBi= Nil |
|
|
|
|
(3.25) |
||||
где |
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Lii = |
J ф/фу-dco, |
Ni = |
j* со/7 (со) ф4- (со) dco, |
|
|
|
||||||
|
|
|
о |
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
ф ! = 1, |
ф2 = С0, |
фЗ = |
£ ц . |
|
|
(3.26) |
||||
После решения системы (3.25) получим, |
учитывая (3.24), |
для |
||||||||||||
коэффициента |
концентрации |
k = |
^ 9^тах- : |
|
|
|
|
|||||||
|
|
|
k = 1+ А 1+ Л2л; (0 |
|
(0 • |
|
|
|
(3.27) |
|||||
В частности, для ядра |
со(0 |
в виде- |
|
|
|
|
|
|
||||||
|
|
|
|
to (t)=a+be~at |
|
|
|
|
(3.28) |
|||||
имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6---- е а+Ь |
|
|
|
||
|
k (t)= 1 |
+ |
А + |
4 |
[ ^ |
- - - (а + |
Ь) |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
ad+pg) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
’ 1+Ра+РЬ |
|
|
(3.29) |
1+ Ра (1 + N (1 + Ра + РЬ)
На рис. 66 показана зависимость k(t) |
при |
следующих |
парамет |
|||||||
рах: |
|
|
|
|
|
|
|
|
|
|
•0,4; |
Е1= |
10е — |
; К, = |
Ю6 — |
; m = |
0, 1; |
(3.30) |
|||
|
а = |
0,01; 6 = 0,99; |
а = 1. |
|
|
|
||||
Рассмотрим теперь задачу о слоистой трубе из вязкоупругого |
||||||||||
простого композита под действием нагрузки |
(6.6.1), показанной |
|||||||||
|
на рис. 27 (с. 177, |
|3= я/16). Для такой зада |
||||||||
|
чи |
может |
быть |
построено |
точное |
решение. |
||||
|
Ядро |
со(0 |
для |
связующего |
было |
выбрано |
||||
|
в виде |
(3.28), где |
|
|
|
|
|
|||
|
|
|
|
а = 0,1; |
6 = 0,9; а = 1 . |
(3.31) |
||||
|
Аппроксимация (3.13) выбиралась в виде |
|||||||||
|
va — AiaCO -Ь Ага + |
' |
■+ |
1+ |
м |
|
-р2й) |
|||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
(3.32) |
Для определения |
неизвестных |
коэффициентов |
|
А и, |
А2а, А3а, |
|||||
А 4а, А5а и «плавающих» |
параметров |
|
(Зг решалась система семи |
алгебраических уравнений (мы отбросили для |
простоты |
записи |
|||||
индекс |
а ) : |
А4 |
|
|
|
|
|
|
*'> = А ^ п + А2 + |
|
|
|
(3.33) |
||
|
1 + р 1(й(/) |
|
1+ М '-) |
||||
|
№ |
|
|
|
|||
Е сли в выражении (3.32) зафиксировать один |
из |
параметров |
(3, |
||||
например рь оставив плавающим только |
то в системе |
(3.33) |
|||||
будет только 6 неизвестных, и она является |
линейной. Если |
же |
|||||
оставить «плавающими» оба параметра (3i и |
(32, то |
система |
(3.33) |
||||
будет |
нелинейной, и ее можно решать методом |
итераций. |
На |
||||
рис. 67 и 68 показаны графики изменения |
со |
временем |
макси |
||||
мальных напряжений ве/р в вязкоупругом слое |
(рис. 67) и в упру |
||||||
гом слое (рис. 68) слоистой трубы, состоящей |
из |
5 пакетов (у = |
|||||
= 1/2) |
при различных значениях |
параметра |
М (3.19). Буквой э |
помечены кривые, соответствующие решению по теории эффек тивного модуля, буквой т— кривые, соответствующие точному ре шению, нулем — кривые, соответствующие решению по теории нулевого приближения. Штриховой линией показаны асимптоты соответствующих кривых при о.
§ 4. Нелинейные задачи
Связанная задача термовязкоупругости даже для линейной теории является нелинейной за счет функции рассеивания (1.4.39). Эта задача для квазистатического случая заключается в решении трех уравнений равновесия
[ & /« Й (ukii - a ki(x)T)},i + X;= 0, |
(4.1) |
где под Т понимается разность между текущей температурой и температурой актуального состояния, и уравнения теплопровод ности (1.2.31)
р (х) ср Сх) |
*L = [ХГ/ (х) Т.,1, + |
Рq + w |
- |
|
at |
|
|
- Г * |
[а,7 (х) Rm (икл- |
аш (х) Г)] |
(4.2) |
at |
|
|
|
при выполнении граничных условий |
|
|
|
и, |s, = lie, |
Rijik (х) (ик,1 — « « (х) Т) П/Is,Si, |
||
(aMXj,(x)Tjnt + b w T)Uq = xM |
(4.3) |
||
и начальных данных (1.2.33) |
|
|
|
|
при t= 0 :Т — Т°(х). |
|
(4.4) |
Решение этой задачи можно искать описанным ранее методом. Методика определения эффективных характеристик: эффективных тензоров ядер релаксации, теплопроводности, теплового расшире ния уже обсуждалась.
Поэтому для теории эффективного модуля получаем связанную задачу термовязкоупругости в виде
hi jkiVk.u — Р//Ф,/ + |
= |
0 , |
(4 .5 ) |
||
(Рср) ^ 7“ = |
А , / |
+ р<7 + W |
------- |
[р,,и,-./ — РО]. |
(4 .6 ) |
v, Ь, = |
u°t, |
(huktVk.i — р./'О) Щ[т, = S°t, |
(4.7) |
||
(а^Л ^Ф ./Я,+ 6,,>ф)|г, = |
т<«>; |
|
|||
|
при t = 0: ft=T°, |
|
(4.8) |
||
причем введены следующие обозначения: |
|
|
|||
f>u = h ilktau, Р = |
Pi;a i7,. |
(4.9) |
где а*н — компоненты эффективного тензора теплового расшире ния, Л// — компоненты тензора теплопроводности.
Функция рассеивания W* в теории эффективного модуля вы
бирается в виде |
|
|
|
|
|
|
¥• = б?,-J-(v,i - |
* ■ $ )- ± |
Нг,ы (0) ± - |
(erf.о/,), |
(4.10) |
||
где |
|
|
|
|
|
|
|
|
erf/ = |
hijkPkA — |
|
(4.11) |
|
После решения |
задачи |
(4.5) — (4.8) |
можно найти |
напряжения по |
||
теории нулевого приближения |
|
|
|
|||
|
|
с $ |
|
|
|
(4.12) |
Рассмотрим |
теперь |
квазистатическую задачу |
главной |
квази |
линейной теории вязкоупругости. Определяющие соотношения для
случая |
|
изотропных компонентов композита |
запишем |
согласно |
(1.4.44), |
(1.4.45): |
|
|
|
|
|
t |
|
|
|
(Уи = j |[г(Я t — Т) — Гф(х, t — т) ф(х, е, 0)] ги (т) + |
|||
|
|
о |
|
|
+ |
[ г , Й / - т ) - Г фЙ * - т ) ф ( х , е, 0 )----- L |
( r C M |
- T ) - |
|
|
|
— Гф(х, t — T)<p(x, е, 0) J 0(т) б,,- j |
da = |
|
|
|
[Г (х) — ГФ(?) ср (х, е, 0)] е,7 + J Г, (х) — |
|
— Гф (х) ф (х, е, 0)-----(Г (х) — Гф (х) ф (х, е, 0))J 06,7. |
(4.13) |
Если шаровые части тензоров напряжений и деформаций связаны по закону теории упругости
|
|
<т=К(*)6, |
(4.14) |
|
то соотношения |
(4.13) |
примут вид (1.4.48): |
|
|
|
<Уц= |
[Г (X) — Гф (х) ф (х, е)] е,, + |
|
|
+ |
| А" (х )-----[Г (*) — Г, (х) ср (х, е)] j 68,,. |
(4.15) |
||
Рассматривая теперь соотношения (4.13) или (4.15) как опе |
||||
раторные |
|
|
1 |
|
|
|
<*ц = |
е), |
(4.16) |
можно повторить все выкладки § 2 предыдущей главы. |
|
|||
Тогда для теории |
эффективного модуля в случае1 слоистых |
|||
вязкоупругих композитов получим задачу |
|
|||
|
|
Gii.i + |
-X/ — 0» |
(4.17) |
|
|
n, k = Ub |
<Ji,n,k = S°, |
(4.18) |
причем тензор напряжений (индекс «э» будем опускать) связан с тензором деформаций соотношениями
|
|
|
°V/ = |
— бз/^з/) + 0Г3363А / + |
Pii + 2Qtj, |
|
(4.19) |
||||
где |
|
|
|
|
|
|
|
_ |
|
|
|
~ |
|
|
|
|
|
|
Lfi (t |
|
|
||
* = j* [7,х (t — т) 6 (т) + L2 {t — г) е33 (т) |
т) fl (0, е33,!р2, q2)] drt |
||||||||||
|
о |
|
|
|
|
|
|
|
|
|
|
|
*33 — J [7.2 (^ |
Т) ® (Т) “Ь 7.3 (^ |
Т) е33 (Т) |
Lf2(t |
т) X |
|
|||||
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
/ 2 (0, е33, р\ q2)] dr, |
|
|
|
||
|
Ри'= j |
[7.4 (* — т) — 7,f3(t — t) /3 (0, e33,p2, q2)] pu (T ) dr, |
(4.20) |
||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
Q/j = |
j |
[LB(^ — T) — Lfi(t — т )/ 4(?, e33, p2, q2)]qu {rydr, |
||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
-а |
аргументы |
функций |
fa, a = l, 2, |
3, 4, |
зависят от |
времени т, |
|||||
тензоры |
р, |
q и величины р2, q2 определены |
формулами |
(7;3.20), |
|||||||
(7.3.21), |
|7.3724), |
(7.3.25). К (4.17) — (4.20) |
следует добавить еще |
||||||||
я |
соотношения Коши (1.2.1). Заметим, что функции fa, a = l,2 , 3,4, |
независимыми не являются.
Для упрощенной главной |
квазилинейной теории |
упругости |
|
трансверсально изотропной среды следует положить |
|
||
W |
s |
= 0, |
(4.21) |
/з = /з(р2, <72). |
U - W , ? ) , |
(4.22) |
так что определяющие соотношения имеют вид
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
<*п = |
J |
{in « |
- |
г) в (г) + |
L2 (/ - |
х) е33 (т)-+ |
- L |
[ Ц 0 - |
т) - |
||||
|
О |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
— |
1 ,3 0 — |
Т) /з ] [е ,1 (т) — е 22 (т)] J d x , |
|
||||||||
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
сгаа = |
| |
( l 2 0 |
— |
т) 6 (т ) |
+ |
I * |
0 — |
т ) е 33 (т ) + |
- J |
[L 6 (t — |
x) — |
||
|
|
— |
L , 3 ( t — |
x ) f t ] |
[e22 (т) — ец(т)] j |
d x , |
|
|
|||||
|
o 33 = j [ L 3 0 — T ) 9 (T ) + L 3 ( t — т ) e 33 (T )] d x , |
( 4 .2 3 ) |
|||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
j |
{ U ( t ~ |
x ) |
— |
L f 3 ( t |
— |
x ) /3] e„(T) d x , |
|
|||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
<*13 = |
j |
[ L 3 0 - |
T ) |
- |
Z.,4 0 |
- |
T ) / , ] e 13 ( T ) |
d r, |
|
||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
<J23 = |
J [^5 |
T) |
|
Ln (t |
|
т) / 4] e23 (T) dr. |
|
Для «простейшей» главной квазилинейной теории вязкоупру гости трансверсально изотропной среды в (4.22) и (4.23) следует положить
/з = /з ( Р 2) , / 4 = Ы <72) . |
(4 .2 4 ) |
Н Е К О Т О Р Ы Е Л И Т Е Р А Т У Р Н Ы Е У К А З А Н И Я
Обзор по теории вязкоупругих композитов имеется в работе [96].
§ 1 . Метод осреднения к решению квазистатических и динамиче ских задач линейной теории вязкоупругости применен в [84].
§2. Эффективные вязкоупругие характеристики простых компо зитов рассмотрены в '[84], непростых композитов — в [73]. Метод канонических операторов предложен в [85].
§ 3. О методе численной |
реализации упругого решения говорится: |
|
в [84]. Решение задачи о концентрации напряжений в анизо |
||
тропной пластинке с |
отверстием в упругом случае дано |
в |
[54], в вязкоупругом |
— в [2]. Эффективные вязкоупругие |
|
характеристики в плоском случае приведены в [79]. Задача о* |
||
слоистой вязкоупругой трубе под действием локальной |
на |
|
грузки решена в работе [28]. |
|
§4. Постановка связанной задачи термовязкоупругости для ани зотропных сред дана в работе [77], различные нелинейные теории вязкоупругости рассмотрены в [38, 78].
Г л а в а 9
КОЛЕБАНИЯ И ВОЛНЫ
Описывается характер распространения плоских гармониче ских волн в неограниченной среде. На примере неоднородного уп ругого стержня демонстрируется техника осреднения в динамиче ских задачах. Далее эта техника применяется к пространственной динамической задаче теории упругости и линейной вязкоупруго сти. Описывается явление волнового фильтра. Обсуждаются неко торые вопросы разрушения композитов.
§ 1. Динамическая задача об упругом неоднородном стержне
Рассмотрим динамическую задачу МДТТ для композита. Она заключается в решении трех уравнений движения (1.2.10)
( ы ) •относительно 3-х компонент вектора перемещений при удовлетво
рении граничным условиям |
(1.2.9) |
|
|
|
|
« ( Is. = “ ( . 3 ц (Х,~и) п , |z, = S; |
|
|
( 12) |
||
и начальным данным |
|
|
|
|
|
при |
t = |
0 : |
|
|
(1-3) |
К задаче (1.1) — (1.3) |
может быть, как и для |
квазистатиче- |
|||
•ского случая, применена техника осреднения. Прежде |
чем |
это |
|||
сделать, заметим, что |
при изучении динамики МДТТ |
часто инте |
|||
ресуются характером |
распространения плоских |
гармонических |
|||
волн в неограниченной среде. Для этого рассматриваются |
одно |
||||
родные уравнения движения (1.1) |
|
|
|
||
|
|
|
|
|
(1.4) |
и ищется их частное решение вида
и(х, t) = Ле‘(^ + £°о>