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

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

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

[ Qmn = т + *я, a £

 

 

 

 

 

суммирование распростра-

j

 

 

т.п

 

 

 

 

 

 

 

 

 

 

 

 

няется

по

всем

целым

т и п,

кроме

т = п= 0.

Кроме того,

J]

означает суммирование

по всем

целым

т и п ,

 

 

 

т,п

2* — суммирова­

ние по

нечетным

значениям

индекса. Коэффициенты а0,

b0,

ah,

bk, ch, dh (k=\, 3, 5,...)

действительные числа.

 

 

 

 

Функции

V

- —

— — (k >

3) являются

двояко-периодическими,

 

 

т,п

(* — Ртл)Л

 

_

 

 

 

 

 

 

а функции с (2),

( 2)

и

fk(z) = V -— д”1

 

(k ^

3)

обладают сле-

дующими свойствами:

 

 

т,п

(2утп)к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S ( 2 + 1) — S ( Z ) = Я ,

 

 

 

 

 

 

 

 

 

 

£(г+ 1 )-£ (г) = - « ,

 

 

 

 

 

 

 

 

Q ( Z +

l)-Q(z) = - £ '( z ) -

5S4

 

 

 

 

 

 

 

Q(z +

i)-< 3 (z ) =

( r ( z ) - i - ^ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.4)

 

 

/*(* + 0 — /*(z)= — tJ ]

 

 

 

 

 

 

где S.

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя условие периодичности функций Uj и U2, получим,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

До =

 

~х^

1 ■ ь„ =

а, ^х, +

 

V.

 

(2.5)

Условия (1.23) и (1.24) можно записать для

/еГоо

в виде

 

 

<Pi (0 +

*фГ(0 + ЫО + Yi (Р) *"+ Y2 (Р) t = Ф2 (0 + *фГ(0 + ф^),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 2.6)

х* [х1ф1 (0 -

Щ Г) -

^М01 = * 2<Р2 (0 -

'фГ(0 -

ФЛО-

(2.7)

Здесь

 

 

 

 

 

 

K* =

\bi\x,ь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Yi (Р) = (Сад —Сад —Сад + СпРр)/2,

(2.8)

Y2 (Р) = (Спрр—Спрр + Сад — Сад)/2.

 

 

 

X* fa +

1)

 

 

4n*+2<I„ (A > 1 ),

 

(2. 11)

 

Ck+ 2 =

 

X2

 

 

 

 

 

 

x* +

П=1

 

 

 

 

 

 

4

X»(xt +

1)

ak +

- ^ y A W R b i k -

 

 

 

 

 

X*— 1

 

 

 

 

 

 

 

 

 

 

 

x*fa + l)

(£ +

2)

^

T]nft+2an.

 

 

 

 

 

 

X* +

x2

 

 

 

n=l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Упражнение 2.1. Доказать, что свободные члены системы

(2.9)

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

 

 

 

 

 

 

 

 

 

I.

 

 

 

 

1*14*41*1

 

 

(2Л2)

 

 

k,n=1

Л,п=1

 

 

А,п=1

 

 

 

сходятся, т. е. система (2.9)

является системой нормального

типа.

Упражнение 2.2. Доказать, используя условие периодичности

функций U 1, U2,что

 

 

 

 

 

 

 

 

 

 

 

h\\р + Лггрр = (Сцрр + С22рр) —2nR

clf

 

 

 

 

 

 

 

 

 

 

 

H-i

x

 

 

^22pp — h\\pp =

 

( C 22pp —

Cupp)

+

 

[(1 + Щ) a i

+ Yi (P) # ] • ( 2 - 1 3 )

Упражнение 2.3. Показать, что для задачи Жр? вектор объем­

ных псевдосил

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xi(pq)= Cijpq.j

 

 

 

(2Л4)

имеет представленный в табл. 2.1 вид.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

2.1

 

Hpq)

 

 

х х

 

 

 

 

у

 

у

 

(pq)

'

 

.

 

 

 

 

 

 

 

 

Аз

 

 

 

 

 

 

 

 

 

 

 

 

(И )

 

(Я. + 2р)|4

 

 

 

А/,о

 

0

 

(22)

 

 

Ki

 

 

 

 

А

+ 21*).,

 

0

 

(33)

 

 

^.1

 

 

 

 

^,i

 

0

 

(12)

 

 

И,2

 

 

 

 

 

0

 

(13)

 

 

0

 

 

 

 

 

0

 

 

 

(23)

 

 

0

 

 

 

 

 

0

 

Ц,2

 

В табл. 2.2 приведены значения относительных эффективных

модулей

Лцц/(Ai—f—2fXj), /11122А 1,

^пззАь

Лзззз/A i“b2m)

в зависи­

мости от

х* при

разных

объемных

концентрациях

волокна у.

В расчете принималось v! = 0,39; ^2 = 0,2.

V

0 ,4

0 ,5 5

0 ,7

0 ,7 5

0 ,7 8

 

h m i

6 » , ,

^п»>

А»МЭ

X*

*1+ 2ц,

к

К

U + 2ц»

 

6

1 ,4 2

1,11

1 ,0 3

1 ,7 5

20

1 ,7 7

1 ,3 6

1 ,2 6

4 ,2 6

120

1 ,9 6

1 ,4 8

1 ,3 9

2 1 ,6 2

400

1 ,9 9

1 ,5 0

1,41

7 0 ,1 0

6

1 ,6 7

1 ,1 2

1 ,0 4

2 ,0 3

20

2 ,4 2

1 ,5 0

1 ,4 6

5 ,5 1

120

2 ,9 3

1 ,6 8

1 ,7 2

2 9 ,4 0

400

3 ,0 3

1 ,7 0

1 ,7 7

9 6 ,0 6

6

1 ,9 9

1 ,1 4

1 ,0 7

2 ,3 1

20

3 ,6 8

1 ,7 5

1 ,8 4

6 ,7 9

120

5,71

1 ,9 3

2 ,6 2

3 7 ,3 2

400

6 ,2 4

1 ,8 7

2 ,8 0

1 2 2 ,2 0

6

2 ,1 2

1 ,1 6

1 ,0 7

2 ,4 1

20

4 ,4 4

2 ,0 2

2 ,0 8

7 ,2 4

120

8 ,7 9

2 ,3 7

3 ,6 5

4 0 ,1 4

400

10,51

2 ,1 0

4 ,1 7

131 ,17

6

2 ,2 0

1 ,1 7

1 ,0 8

2 ,4 6

20

5 ,0 7

2 ,3 2

2 ,3 0

7 ,5 3

120

14 ,8 4

4 ,0 5

5 ,8 3

4 2 ,2 3

400

2 9 ,9 5

3 ,8 7

8 ,7 0

1 3 7 ,5 0

§ 3. Решение антиплоской задачи

Задача Ж 1 3 (1 .26) — (1 .28)

(см. табл. 2 .1)

называется

плоской задачей теории упругости.

будем опускать

В этом параграфе для

сокращения записи

индексы (1 3 ). Функции W ,

W

будем искать в виде

t/S= M

t ) + i ! r Re'Pi(z)’

 

(3.1 )

u h

------- Re ф2 (z),

 

Ра

где

 

ф1 (г) = а М (*)

+ j Е °* ( т л и " )

(3.2 )

* ( * ) = 2 4 ( -* = £ = -)* k=l

Константу А определяем

из

условия периодичности функции С/3:

 

А = ------— nRzat.

 

(3.3)

Тогда

 

 

 

 

 

 

-

Ч> = -

n F t a i +

Я»; (*).

 

°и -

iah =%(*)■

 

 

 

Условия (1.27) и (1.28) на границе Г<ю можно записать

в виде

А - ^ - + —

R e«M 0 = —

 

R e«h(0.

(3.5)

IK

(X!

 

p2

 

 

Im q>! (t) + ( HI

 

+ Hi -

ц„)

 

= Im <p2 (t).

(3.6)

Используя разложение в ряд Лорана функции q>i(z) в окрестно­ сти Г<хь получим следующую систему уравнений для определения коэффициентов а^:

e» - x £ *

^ . = Jr ^ - * 6 u ,

(3.7)

П=1

 

 

р - f c L

l Rk+iTk+ u n = \ ,

(3.8)

 

 

*l4 r t . » > з .

_

(

я,

* = 1,

(3.9)

*+‘

l s * f l ,f t > 3 ,

 

v __

1—X*

(3.10)

 

 

 

 

1+Х* •

Коэффициенты ch выражаются через ak следующим образом:

Ck = *'(ak + £ Pnifln)-

(3.11)

П—1

 

Можно показать, что система (3.7) является системой нормаль­ ного типа.

Упражнение 3.1. Доказать, что

 

 

^1313 _ J

%Уа1

(ЗЛ2)

Pi

R\it '

 

В таблице 3.2 приведены значения относительного эффекти ного модуля /1131з/ц1 при разных объемных концентрациях волок-

V

0 ,4 0

0 ,5 5

0 ,7 0

0 ,7 5

0 ,7 8

0 ,7 8 5

ч* Чжень

6

1 ,805

20

2 ,1 4 7

120

2 ,3 1 4

400

2 ,3 4 0

сю

 

6

2 ,3 2 6

20

3 ,1 8 4

120

3 ,5 5 5

400

3 ,5 7 8

сю

 

6

3 ,1 7 6

20

5 ,2 2 2

120

6 ,9 4 5

400

7,291

 

6

3 ,6 2 0

20

7 ,0 0 6

120

1 1 ,1 7 0

400

1 2,5 2 3

сю

 

6

 

20

-

120

400

 

оо

 

6

-

20

Адамс

Ван Фо Фы

(3.12)

1 ,7 9 6

1 ,8 00

1 ,8 04

2 ,1 3 7

2 ,1 3 4

2 ,1 4 5

2 ,3 0 5

2 ,2 9 7

2 ,3 1 3

 

2 ,3 2 2

2

,3 3 9

2 ,3 4 3

2 ,3 3 3

2 ,3 5 0

2 ,3 0 4

2 ,2 9 4

2 ,3 2 5

3 ,0 4 4

2,98 1

3 ,0 7 7

3 ,4 6 9

3 ,3 5 6

3 ,5 0 6

 

3 ,4 1 7

3 ,5 7 7

3 ,5 7 0

3 ,4 4 4

3 ,6 0 9

3 ,1 6 3

3 ,0 0 0

3 ,1 7 3

5 ,1 8 7

4 ,4 5 4

5 ,2 1 3

6 ,8 7 8

5 ,4 1 9

6 ,9 2 9

 

5 ,5 8 9

7 ,2 7 3

7 ,3 7 2

5 ,6 6 6

7 ,4 3 2

3 ,6 4 6

3 ,3 0 7

3 ,6 1 9

7 ,0 0 5

5 ,2 2 2

7 ,0 0 4

11,035

6 ,6 2 2

1 1 ,1 6 4

 

6 ,8 8 2

1 2 ,2 2 6

12,540

7 ,0 0 0

12,751

 

3 ,5 1 6

3 ,9 7 7

-

5 ,7 9 6

9 ,4 2 7

7 ,5 8 7

2 3 ,6 7 0

 

 

7 ,9 3 2

3 1,022

 

8 ,0 9 0

35 ,9 3 3

-

3 ,5 5 5

4 ,0 4 7

 

5 ,9 1 0

10,227

на v, вычисленные различными авторами, а

в последнем

столб-

це — по

формуле

(3.12).

В

предпоследнем

столбце вычислены

значения по приближенной формуле

 

 

 

 

 

 

 

Лип

1-XV

 

 

 

(3.13)

 

 

 

Pi

1 + XY ’

 

 

 

 

предложенной Ван Фо Фы.

рассмотрим

решение

 

^

В заключение

параграфа

задачи Ж 12,

которое мало чем

отличается

от решения задач Ж рр

(Р —

А

Индексы

(12) для краткости

будем опускать. Функции и к (а,

К = 1, 2)

ищем в виде

 

 

 

_____

 

 

 

2ра ( ^ + ЯА) = ^афа(^) - Z $ a (z) — $a(z),

 

(3-Н)

где

 

 

 

 

 

 

 

(3.15)

 

 

ф а (г) =

1'ф а (z), фа (z) =

*'фа (z),

 

 

 

 

 

а Фа (а) и г|)а(2)

имеют

вид

(2.3). Из

условия

периодичности

функций £/ь U2 находим, что

 

 

 

 

 

 

 

 

 

 

а0 =

--------60=

- a

1 (x 1- - ^

- U -

я2

/

 

(3.16)

 

 

Xi т 1

 

 

 

 

V

 

 

 

 

Условия идеального контакта (1.32), (1.33) примут вид

 

ФИО —<ф!(0 — ФИО—(Hi— Иг)^= ч>«(0 — <ч>2(0 — ♦«(О. (3-17)

« 'K < P i(0 +

<ф|(0 + ЫО]= Ф»(0 + <фг(0 +

’МО-

(3.18)

Для определения

коэффициентов

ah получим

 

следующую

нор­

мальную систему уравнений:

 

 

 

 

 

 

 

 

 

 

a k + ^ o a A f t

+

n=i

nk +

B G nk) =

E

(Hi —

Иг)

 

(3 .1 9 )

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л0 = в ( - ^ -

— хЛ v; B =

 

x>~

!

 

 

 

°

 

 

 

 

 

V

 

 

 

1+K*:;*xt

 

 

 

л =

 

в

X2 —X*Xi

;

E =

-

1

 

 

 

(3.20)

 

 

 

 

X* + x 2

 

 

 

1-ГХ*Х!

 

 

 

Коэффициенты 6/t,

, dh выражаются

через

ак следующим об­

разом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&i =

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

Ь/е+2 = kak +

^

 

Л/г к+2ап (& ^

1 )>

 

 

 

 

 

 

X* (XI +

1)

 

 

 

 

 

 

(3.21)

 

 

 

 

 

Х2 +

1

Ё

-Т1"А "

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П=1

 

 

 

 

 

 

Ск+2 =

X* (X! + 1)

ОО

 

 

 

 

 

 

 

 

Т\пк+2С1п (£ > 1 ),

 

 

 

х* + х2

п=

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

dk= j q * ± l L a k

 

X* (Цх — ц2)

 

 

 

 

 

1 — X*

М и -

 

 

 

 

 

1—X*

 

 

 

 

 

 

 

 

 

 

X* (х, +

1)

(^ +

2) ^

 

т]п. к+2fln.

 

 

 

 

 

X* + х 2

 

 

 

 

 

 

 

 

 

 

Упражнение 3.1. Показать, что

 

 

 

 

 

 

 

 

 

 

Иг

 

 

 

 

 

 

лгр-1

 

 

 

(3.22)

 

 

 

 

 

 

 

 

 

 

 

 

В таблице 3.2 приведены значения относительного эффективного

модуля Л12121 при

разных

объемных

концентрациях

волокна

Y (vi = 0,39;

V 2 = 0,2)

 

 

 

 

 

 

 

 

Т а б л и ц а 3.2

X*

 

 

V

 

 

0.4

0,55

0.7

0,75

0,78

 

6

1,5 82

1,951

2 ,6 4 4

3 ,0 8 2

3 ,4 9 3

20

1 ,7 4 3

2 ,2 9 0

3 ,6 4 6

4 ,9 4 6

7 ,1 0 2

120

1,8 1 0

2 ,4 4 8

4,281

6 ,6 1 6

14T261

400

1,8 2 0

2 ,4 7 2

4 ,3 9 2

6 ,9 6 9

1 7 ,1 0 5

00

1 ,824

2 ,4 8 3

4 ,4 4 2

7 ,1 3 4

18,770

В силу того что в выражения компонент эффективного тензо­ ра модулей упругости (2.13), (3.12), (3.22) входят только кон­ станты аь сь можно дать явные выражения компонент эффек­ тивного тензора модулей упругости (точное аналитическое пред­ ставление) :

^1111 — ^2222 — (Cllll)

+

 

(

У2( 1)

+ a 2v ,(l)j

 

Y 1ai

 

-

 

 

 

 

 

 

l

 

Pi

 

 

 

,

 

.

1

 

Y1(D

 

 

 

^1122— (C1122) +

Y |ai

Pi

‘ a2Vl(l)} ■

 

 

 

 

1

 

Y*(1)Y.(3)

 

 

 

 

 

 

\

 

»

 

"1133 — "2233 — (^1122/

~r Ya i

Pi

 

 

 

 

 

 

 

 

 

 

(3.23)

 

 

 

 

 

 

 

Y? (3)

 

 

^зззз —

( Q i i i )

+

Yai

 

 

 

Pi

 

 

 

 

 

 

 

 

 

 

 

 

 

^1212

 

 

(xi +

1) XiY _

 

 

 

 

1*1

1

x .- x ? n v - V .

 

 

 

_ __ ^2323

11'

 

 

 

2XY

 

 

 

Pi

Pi

 

 

 

1 + % y - X 9NiY % '

 

где

 

 

 

 

 

 

 

 

 

 

x2 — 1

|Y

1 — x,i(«i +

1)

*»-

- N J - ' N A ,

 

2ao

2a0

 

 

 

 

XI + 1

 

 

 

 

 

 

(i+ x * x 1) ( x , - X i W

'_1£/«)

 

 

 

„ _

x2 — X*Xi

 

 

_

 

1 X*

 

 

(3.24)

a3 —

-----“

 

» M —

1+ X*Xi

 

 

x* + x2

 

 

 

 

 

X, = 1 + Xj | ( * ! +

Y + Gn — < V n j .

% a = 1 + ^1 |

-------^7~)

Y

Оц

a 3ru j ■

Компоненты матриц бесконечного порядка

(T)|s), К, (a|s), (Л (tt*s),£/a(«ft) .

^(rij,), V2(vl), Y(yls), Z(zls), V(vls), U(uls) выражаются по формулам

Ills —

Л 14 s - b

y !s —

G l 4s+l

+

« З г 14s+li

 

wjs =

G\4s+i — a3rj 4S+1',

=

Л4/-1 i;

 

 

 

vtt ~

^4t+11+

«3r4/+l i;

 

 

 

 

 

=

G*t+l 1--- a3r4t+l Ь

 

 

 

Ун —& 14s—i — XV4t—14s—i — X2/ ;

(3.25).

2/s = ^4f—1 4s—1

+ X i ( G ^ —] 4s—1

 

 

 

1-I- K*Ki -

 

&3r 4t—l 4s—1 '

an//

 

 

 

 

 

 

 

 

 

 

 

U*s =

^ 4/+ l 4s+l —

X j (G 4/-I-1

4s+l

+

« 3 Г4 Ж

4 s+ l);

 

Uts — ^ 4/+ l 4s+l

+

Xx (G 4/+ I

4S+1 —

<*3r 4t+\

4 s+ l);

 

I

=

T\4t-l

1T]l 4s—1 (t, S =

1,2, 3, . . .).

 

Аналогично можно получить результаты и для неидеального кон­ такта. Приведем здесь только выражение для модуля /i!3i3 в слу­ чае, когда волокно проскальзывает относительно матрицы

 

=

1 --------------

^ ----------

(3.26).

 

(*i

l+ Y -^ J V 'A T ,

 

где компоненты матрицы Yo(y0is) имеют вид

 

 

= б«-1 4s—1— Г4t—l 4S-1 + /.

(3.27)'

В таблице 3.3

приведены

значения

эффективного модуля

(3.26)

в зависимости

от объемной

концентрации волокна у.

 

V

Т а б л и ц а 3.3

0,78

Л1313/ Н-1

0 ,4 2 5

0 ,2 7 7

0 ,1 3 4

0 ,0 7 8

0 ,0 2 7

Нетрудно провести подобный анализ и для произвольного па­ раллелепипеда периодов.

§ 4. Модельные задачи

Рассмотрим применение численного метода, основанного на использовании матрицы А. А. Ильюшина, к решению задачи тео­