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

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

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

 

 

о>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) = Ле‘(^ + £°о>