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

книги / Механика композитных материалов. 1982, т. 18, 2

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

Принимая 7 = г/(а + 1р), r)=i|jp и сравнивая (12) и (13), находим

 

/(*/<, у) —

1

Яр=— (Р> 0).

 

— ivy + ixit

 

 

 

 

Значение

постоянной В0 определяем

из

(12) с учетом f{t,y) из (14):

В0 = п/а.

 

 

 

 

Таким

образом, выражение |'i (11)

приобрело полную определен­

ность. Его распространение на отрицательные у легко выполнить с уче­

том тождества f ( - x h, - у )

= - f ( x h,y). В итоге имеем

 

2ni

neWmt] у < 0;

 

L

 

 

/г = —оо

i

(15)

 

У > 0 ,

 

 

где е 0= -g-H е т = 1 при т ^ 1 . Ряд по k для h удовлетворяет уравнению

(6) и является абсолютно сходящимся во всей плоскости Оху, за исклю­ чением точек x = ka (k = 0 , ± 1, ± 2,...), при условии объединения слагае­ мых с одинаковыми положительными и отрицательными индексами k

при его суммировании. Точки x = ka являются полюсами функции gb т. е.

есть внешнее периодическое по х и ограниченное при |#|->-оо решение уравнения (6). Исходя из (15), получаем следующую искомую систему линейно независимых решений:

 

1

 

^птп^

0‘,

 

771= 0

 

(16)

|п =

£

оо

 

 

к = —оо

 

 

 

 

 

 

 

 

771=0

 

 

^тфтп71 1

( - 0 П

(/г= 1, 2, 3,

.).

)71771 — а

( п - 1)!

 

 

С помощью полученной системы (16) составляем еще одно решение урав­

нения (6):

оо

Ф1(2)= Е ^ .

(1?)

71=1

 

Решение (17) в сумме с (10) будет использовано ниже для представ­ ления ф[ (х,у) в (9). Но прежде чем переходить к решению поставлен­

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

и £п к поляр­

ным координатам р, ср.

е±'’*’п? осуществляется представлением

ее в виде

3.

Преобразование

произведения двух экспонент, разложением последних в степенные ряды

и их перемножением. В результате имеем

 

 

 

оо

_____ i

/

 

 

e±j^nl= У1, £-PIP ^ ± nj)71~j/ l+ v 2p^ ег’рф;

е= "|/ 1 —tv

(18)

р = —оо

1

+ iv

 

где Ip — модифицированная цилиндрическая функция первого рода. Для представления в виде ряда Фурье по переменной <р на окружности р запишем выражение для | п в форме

 

 

 

1

1

] -= а0(п) (р) +

ь

л-1

(£-Ао)» +-(£+Ла)«

 

 

 

 

 

+

У , [as<n>(р) e isv + a_s(") (р) e-«'s<p].

(19)

 

S = 1

 

 

 

 

Разложение

в ряд

Фурье

выполняется

в такой

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

ности:

 

 

2п

 

 

 

 

£-" = (д:—vr/)-" = ----[(1 + tv)e“p + (1 -

tv)e-<<p]-n=

 

 

 

Р7‘

 

 

 

 

 

On

p —irup

 

 

 

 

= ---------------- (1 + e2e~2iv) -» =

 

 

 

 

(l+ iv )B

 

 

- (

- /

_L V (_

 

(20)

\

n-> ^ 1

* s i r n - m

 

 

1 + VV

s = 0

 

Здесь принято во внимание неравенство |е| < 1 .

Обозначим через Wn)(p) коэффициенты Фурье суммы по k в выраже­ нии для In (19). Тогда

 

оо

2 л

 

 

Ь-

 

 

 

 

 

оо

2 л

 

 

(21)

 

 

 

2 nbs(n>< р ) - Е

[ J

 

-б/ф+ J

rfcp]

s = 0, 1, 2,

Ь=1

О (Гt,-ka)n

0

(l + ka)n

 

Но

 

 

 

 

 

(%+ka)n=

 

[p(l +Jv)e2,'44-p(l —tv) +2&ae,'4’] n;

(llPka)n= (-^-e-,<p)

[p(H-fv)e2i,*’ + p(l —tv) =F2^ae,<i>] n.

Поэтому между &s<n)(p,v) и 6_5<n'(p,v) имеет место соотношение &s(n)(p,v) = 6_3<n)(p,v). Интегрирование же по <р в правой части первого равенства (21) заменим контурным интегралом вдоль окружности ра­ диусом р в комплексной плоскости z=x + iy:

 

2 л

 

 

 

 

 

.(п)

f

e is(p

,

2 ”

X

z n+*-irfz

6_STA= J —---—— dm= ----- Ф

 

 

Jg

tt=Fka)n

v

»n.

J

 

 

 

 

 

'Ps

 

[ (1+ fa)z2zF2 kaz + p2 (1 —tv) ]"

( 22)

Таким образом, задача свелась к вычислению контурного интеграла от рациональной дроби. Корни знаменателя равны

± k a ± ~J k2a2p2(l+ v 2)

Z\,2m) =

1 +tv

Так как 2i^/l)e2(“,t)= p2e2; 2i(,i)22(ft) = p2e2 и |е |< 1 , то для любого k подын­ тегральное выражение (22) имеет внутри окружности р по крайней мере один полюс. При 2 р < а и |v |^ l * таких полюсов всего один, а именно

ka— У k2a2р2(1 +v2)

, (-*) = ------------------------------ ;

1 +/v

( * = 1 , 2 ,

—ka+ У k2a2 р2(1 +v2)

' = ---------------------------------

1 +iv

) .

Поэтому для вычисления контурного интеграла (22) достаточно восполь­ зоваться теоремой вычетов с учетом того, что подынтегральная функция имеет полюс п-го порядка [9]. В итоге имеем

 

 

 

 

 

1

 

1

 

dn~

 

г

2 n + s - l

 

 

-I

 

 

 

 

 

 

 

 

 

 

 

п- l

 

 

 

 

 

 

b -( l - h =

2 n ^

 

 

 

 

 

 

-----------------

 

2

(-ft);

 

 

 

ps

(l+ tv )n

(п- l ) !

-

 

L

z

zJ

 

 

 

 

 

 

 

dzn~l

(Z— Z2

(~h))\n

 

J Z=Z|

 

 

 

 

2™

1

 

1

 

dn~l

Г

g

 

 

-I

 

 

 

 

 

 

 

zn+s~l

 

 

 

 

 

 

=

2 *

 

 

 

 

 

dzn~l

*-

I n _n.lli)\n

J

 

(ft),

 

 

 

ps (l + tv)n (я—1)!

L

( z - z i ^ ) n

J z=z>

 

Раскрывая

здесь

производные,

нетрудно

установить,

 

что

 

й_5,/4(п)=

= ( —l) n+s6,-s,-ft(n), а коэффициенты Фурье (21) равны

 

 

 

 

 

 

 

 

[ ( - ! ) " + (-!)« ]

(п + з

- 1)!

71-1

 

. 2 i(2 n — j — 2 ) !

 

 

 

 

 

 

 

 

 

2n-'p s(l + t'v)3

[ (/г—1)!]2

j=0

С ^ -

 

 

 

 

 

X

 

 

 

 

 

(n+s —j 1)!

 

 

 

x

oo

/

 

-|/

 

 

 

_

\ 7l+ S_ j _ l

 

 

 

 

 

 

 

 

V

U a -

Г fe2a2 —p2(l+ v 2)/

 

 

 

 

 

 

 

 

 

 

 

 

 

( V ^ - p ^ l + v 2))2" ^ 1

 

 

 

 

 

 

 

 

 

b / n ) =

e -2sb _ s(n)

 

(S= 0 ,1 ,2 ,

).

 

 

 

 

(23)

Принимая теперь во внимание равенства

(23)

и (20), запишем разложе­

ние (19)

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

|„ =

^

as<»>e-»'<n+2s>4>+

6*<">eis<Р;

2р<а;

 

 

 

 

 

 

 

s=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(24)

 

 

 

, . .

(s + n —1)!

/

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

a«(")= ( —l)s—---------- -— е2М ----- — )

----

 

 

 

 

 

 

 

 

 

 

s! (п—1)!

 

'

1+iv-4- '

р71пп

 

 

 

 

4.

 

Возвращаясь

к

решению

краевой

задачи

(6),

(8), представим

<Di (х, у)

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф ) =

А п \ п +

 

В п е -< * п « -и л .)+

 

С пе**п (t-vfts)

 

 

 

 

 

п = 1

 

п=1

 

 

 

 

п=1

 

 

 

 

 

 

(25)

 

 

 

 

 

 

/ii,/i2^>0*

h\-\-fi2 = b.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

выполнения

двух

последних

в

(9) условий

достаточно

воспользо-

*

Неравенство |v | ^ l

можно

считать

всегда

выполненным,

так

 

как

 

И55

 

 

И

 

|v |2= Р44

дело сводится к выбору оси координат х.

ваться соотношениями (16). Опуская промежуточные выкладки, запи­ шем окончательный результат:

 

а0йо=Ьо—5о+ Г;

 

e-WnVlu.bn

 

е-п>„v/n.an

Вп

 

Сп--

 

(26)

e-ii|)nb ctg 0 _ e-i*„v6

 

е~Н>пь ctg 0 _

g-i<>nvb

 

оо

 

оо

 

Выражение (25) с учетом (26) принимает вид

 

оо

оо

 

оо

 

ф 1= Х л п£п+ X dne-Wnt+

Cne^nt;

\y\< b ;

п = 1

77 = 1

 

71 = 1

(27)

 

 

 

gi^nb Ctg 0

 

 

 

 

d n = k n b n \

Сп = A n ^ n J

А п —

~

 

 

 

^Zt|)nv 6 __ei^nbctg 0

Таким образом, получено решение

уравнения (6)

в форме (9), где

Ф\{х,у) определяется равенством (27). Данное решение удовлетворяет всем условиям периодичности рассматриваемой задачи и имеет смысл в любом слое толщиной b, содержащем один ряд волокон. Бесконечная последовательность неопределенных постоянных Ап достаточна для вы­ полнения определенного условия на границе р= R.

Часть Вх + Су решения (9) рассматривается нами в качестве возму­ щающего фактора задачи. Так как Вх + Су= —[В( —х) + С( —у)], то аналогичным свойством должны обладать как Фь так и функция напря­ жений для волокна Ф°. В связи с этим положим в (27) А2п= 0, dn= —сп (п= 1,2,...), а Ф° представим в виде

оо

Ф°=Ф1°+Ф?>; Ф1°= ^ О

2п+1Р2пН(^0);

71 = 0

__________(28)

_____________

Р2П+1 (to) = [to+ %02 - R 2(1 + Vo2) ]2n+4 - [ to - По2- R 2{1 + Vo2) ]2n+1;

to= x - v 0y-

Здесь vo определяется через модули упругости материала волокна теми же формулами, что и v в (10) через модули упругости матрицы. Реше­ ния Pm(to) заимствованы нами из [10].

При выполнении контактных условий (8) необходимо переписать функцию напряжений Ф в полярных координатах р, ф, для чего следует использовать разложения (18) и (24). После некоторых преобразований приходим к следующим алгебраическим уравнениям:

4И'-^)+Ф-^ЧЬ-+

j

ОО

 

 

+

A2h^\Cif-hl) (R) +2e2j+I

Cnhi+i

+ v2) +

/i=0

71 = 1

 

 

CO

 

 

 

+A2n+lb-i2j+l)(R) =-2'^ 2i+l[(l + lA)J^2j+ l(l—t'Vo)2j+1 +

+ 0 —ц )^2Н1(1—iv0)2j+l] ;

(29)

Д= а —2/(YI+ YI) +2р/(у2+У2); 6 = 1+ 217(72—72).

(32)

Таким образом, приведенные модули сдвига выражаются через одну неизвестную бесконечной системой (29). Первое приближение Л,, най­ денное из (29) в результате удержания всего двух уравнений, дает сле­ дующие выражения для 71 и 72:

 

-

-

Г iB + B2 + a 2

 

 

2Ry\k= (1 + р) (х2 + е02Х4)

Р-----------------------

 

-ео 2 ( pH-----+Ц

) ]

+

—I1) («4+ е02хг) X

 

—/р+ Р2+ а 2

 

 

 

 

X [ р -

 

 

 

 

 

а

 

 

 

 

 

2Ry2k = (1 + р) (х2+ ео2х4) [

ip-----— Ье02( i p—

) ]

(1 —р) (х4+ео2хг)

ip-

P - i

-е0‘

 

(33)

а

 

 

 

 

 

 

k = (1 + р)22+ 8о2Х4|2—(1 —р)24+ео2Хг|2'

 

 

 

 

оо

 

 

 

х2= а 0(1>( R ) + b - ,<*» (R) +2е

Д»Ь»Л (iiM l+ v 2/?);

 

 

 

П=1

 

 

 

 

ОО

 

 

 

х4= -&!<*) (R) +

i

X J

nh (iyny l + v 2R ) .

 

 

6 «=.

Данные соотношения позволяют получить в первом приближении эф­ фективные модули сдвига (31) для общего типа анизотропии упругих свойств матрицы и волокон. Для частных случаев анизотропии формулы (31) —(33) принимают более простой вид. Так, например, для изотроп­ ной матрицы Pi = |3 = 0, а, = а= 1 и равенства (33) записываются в форме

2Ryik= —(1 + р) (x2+ eo2x4) [р —1 —ео2(р + 1)] +

 

+ (1 —р) (х4 + 802х2) [ р —1—е02( р + 1 ) ] ;

(34)

 

 

2Ry2ik = (1 + р) (х2+ ео2х4) [р— 1 + ео2(р + 1)] + (1 —р) X

X (х4 + е0'2х2) [р—1 + е02(р + 1)]; Х2 = R

n2R

 

Щ- За2 + 2iR

Дп^тфя

 

71 = 1

 

Задача определения эффективных модулей упругости композитной среды' с изотропной матрицей и анизотропными волокнами рассматрива­ лась в работах [2] и [4] (трансверсально-изотропные волокна) и в работе [И] (ортотропиые волокна). Сравнение первого приближения предлагае­ мой теории с результатами указанных работ представлено на рис. 2—4. При этом эффективные модули вычислялись по формулам

1

1

— р*44= (1-4/7,)->; р:45 = р 54 = 0;

Р 55= 1 —4i/72,

(35)

которые следуют из (31), (32). Структура решетки волокон предполага­ ется гексагональной. Значения упругих модулей сдвига матрицы и воло­

кон выбирались следующими:

G = 1,27 • 104

кгс/см2,

Ц44(0) = !Л55(0)= 4,7Х!

ХЮ5

кгс/см2

([2],

см. рис. 2);

G= 1,61 - 104

кгс/см2,

р,44(0) = р,55(0)= 2,8х

ХЮ5

кгс/см2

([4],

см. рис.

3);

G= 1,15• 104

кгс/см2, ц44(°)= 4,63Х

ХЮ5 кгс/см2, р55(0) = 0,7• 105 кгс/см2

([11], см. рис. 4). Кружки на рис. 3

соответствуют

экспериментальным

данным

[4].

Сплошные линии на

рис. 2—4 отмечают вычисленные значения, а штриховые — результаты работ [2, 4, 11]. Нижняя кривая рис. 4 соответствует отношению 11*55/G. Как следует из рис. 2—4, вычисленные по формулам (35) и приведенные в указанных работах значения эффективных модулей продольного сдвига находятся в хорошем согласии. Следует, однако, заметить, что затрудни­ тельно указать количественные отличия сравниваемых величин, так как в [2, 4, И] результаты приведены лишь в виде графиков. Вместе с тем обращают на себя внимание практически полное совпадение нашей кри­ вой с кривой работы [4] и факт попадания двух экспериментальных то­ чек на эти кривые. Это последнее является довольно неожиданным, если учесть достаточно высокую степень наполнения и точность используемых нами формул (первое приближение предложенного метода).

В заключение заметим следующее. Рассмотренная выше задача о вы­ числении приведенных модулей сдвига в математическом отношении эк­ вивалентна задаче о вычислении эффективных коэффициентов теплопро­ водности для анизотропной в тепловом отношении волокнистой среды. В самом деле, поставив в соответствие перемещению w температуру Г, деформациям ууг и уХ2 компоненты градиента температуры, напряже­ ниям тyz и Txz компоненты вектора плотности потока q, модулям сдвига коэффициенты теплопроводности, получим уравнения стационарной теп­ лопроводности типа (1) — (3). В этой системе уравнения (2) будут выра­ жать закон теплопроводности Фурье, а уравнение (3) — условие ба­ ланса тепла. Кроме этого, выше лишь для упрощения выкладок мы предполагали волокна круглыми. Если поперечные сечения волокон явля­ ются эллипсами, то изложенный метод требует незначительной модифи­ кации. В этом случае необходимо переписать разложения на окружности (18) и (24) в соответствующие разложения на эллипсе и выбрать функ­ ции P2n+i в несколько ином виде [10].

СП И С О К Л И Т Е Р А Т У Р Ы

1.Малмейстер А. К., Тамуж В. П., Тетере Г. А. Сопротивление жестких полимер­ ных материалов. 2-е изд. Рига, 1972. 498 с.

2.Рикарде Р. БЧате А. К. Упругие свойства композита с анизотропными волок­

нами. — Механика полимеров, 1980, № 1, с. 22—29.

3. Хорошун Л. П. Методы теории случайных функций в задачах о макроскопиче­ ских свойствах микронеоднородных сред. — Прикл. механика, 1978, т. 14, № 2, с. 3—17.

4.

Чжень, Чжен. Механические свойства анизотропных волокнистых материа­

лов. —

Прикл. механика (США), 1970, № 1, с. 197— 199.

5.

Шермергор Т. Д. Теория упругости микронеоднородных сред. М., 1977. 399 с.

6.

Головчан В. Т., Гиря М. Г. Распространение упругих волн сдвига в композитной

волокнистой среде. — Механика полимеров, 1979, № 1, с. 146—149.

7.Головчан В. Т., Никитюк Н. И. К решению задачи о сдвиге волокнистой компо­ зиционной среды. — Прикл. механика, 1981, т. 17, № 2, с. 29—35.

8.Бейтмен Г., Эрдейн А. Таблицы интегральных преобразований. М., 1969.

Т.1. 343 с.

9.Смирнов В. И. Курс высшей математики. М., 1969. Т. 3, ч. 2. 672 с.

10.Лехницкий С. Г. Теория упругости анизотропного тела. М., 1977. 415 с.

11.Хант Е. Б. Напряжения в волокне, находящемся в анизотропной матрице, и упругие постоянные материалы с анизотропными волокнами. — Изв. АН СССР. Меха­ ника твердого тела, 1977, № 1, с. 120—130.

Институт сверхтвердых материалов

Поступило в редакцию 28.04.81

АН Украинской ССР, Киев

 

УДК 539.3.001:678.067

И. В. Грушецкий, М. #. Микельсон,, В. П. Тамуж

ИЗМЕНЕНИЕ ЖЕСТКОСТИ ОДНОНАПРАВЛЕННОГО ВОЛОКНИСТОГО КОМПОЗИТА ВСЛЕДСТВИЕ ДРОБЛЕНИЯ ВОЛОКОН

Изменение упругих характеристик изотропного материала вследствие возникновения в нем дискообразных трещин изучено в ряде работ [1—3]. При этом получено достаточно хорошее соответствие между фактически наблюдаемым изменением модуля упругости в результате возникновения повреждений под нагрузкой и теоретически рассчитанным значением из­ менения жесткости [4]. Для анизотропных материалов аналогичный рас­ чет изменения жесткости более затруднителен из-за сложности решения задачи о раскрытии трещин в анизотропной среде. Однако для дискооб­ разной трещины в трансверсально изотропной среде в случае, если нор­ маль к поверхности трещины совпадает с осью симметрии среды, реше­ ние получается в аналитическом виде [5]. Поэтому расчет изменения жесткости такого материала может быть проведен достаточно просто.

Типичным конструктивным трансверсально-изотропным материалом являются однонаправленные композиты. При растяжении такого компо­ зита перед окончательным разрушением в нем накапливается значитель­ ное число разрывов волокон, что приводит к изменению жесткости мате­ риала. Анализ изменения модуля упругости однонаправленного компо­ зита вследствие разрыва волокон дан в статье [6], где результат получен сравнительно громоздкими вычислениями.

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

Рассмотрим описание изменений жесткости однонаправленного ком­ позита, возникших вследствие деформации разрывов волокон. Композит заменяется однонаправленным трансверсально-изотропным материалом, усреднение характеристики жесткости которого определяется согласно формулам теории армирования [7]. Оценим в таком материале радиус дискообразной трещины, эквивалентной разрыву волокна в композите. Пусть модуль упругости матрицы, волокна и композита в направлении армирования равны соответственно £ м, £ в, Ек. Радиус волокна обозна­ чим г и коэффициент армирования р. Если считать £ м<с£в, то £к~р£в- Тогда «размазывая» свойства композита по площади повреждений зоны радиуса а и приравнивая усилия, получим

Екла2= Евтсг2; а = г/]/р.

(1)

Поврежденность композита характеризуем числом разрывов волокон в единице объема. Учитывая, что дробление волокон на части меньше так называемой «неэффективной длины» мало вероятно, удобно считать отрезки волокон длиной б за структурный элемент материала, который может быть в разрушенном или неразрушенном состоянии. Число этих

структурных элементов и, следовательно, максимально возможное коли­

чество разрывов в объеме композита равно N, и если разрушено к kN

структурных элементов, то поврежденность ■P=](jovr' ^ Дальнеишем вво"

дим аналогично изотропному случаю безразмерный параметр

 

kN

 

со= Ра3 = ЮОК

(2)

Поскольку V и N связаны очевидной зависимостью

 

 

N,

(3)

то из (1)—(3) следует

г

 

k

 

100 лбУц

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

о;

(4)

(т , —т2)

2

х = -

(5)

(1+mj) (1+ т 2)

£44(Уи1~ 1/пг)

Постоянные mi,2 и «1,2 связаны с упругими характеристиками материала

следующим образом:

 

^

..

п1,2 = -

 

1-V ,J4,

Г , _Г ,

£ « а -У |га)(1 -У п £ '/£ »)

Г П

1

—vi32Ei/E3

1 1+ 1 ‘ ----------- /?,П -V.. ..12

J I

 

 

Е i ( l — V12V1 3 ) 5

 

 

 

 

C \ \t l i С44

(6)

 

 

 

m t= -

/= 1, 2;

 

 

 

С13+С44

 

 

 

1 - у 132£ , / £ 3

Е^ 3

С44= G13.

Сц = 2G12 I-V 12-V 132Ei/E3

с 1 3 = - 1—V12 — v\32E\IE3

 

Здесь Vjj — коэффициенты Пуассона; йц — модули сдвига; Ei — мо­ дули упругости. Ось 3 совпадает по направлению с осью монотропии. Просуммировав вклад от всех трещин, получаем, что отношение модуля Е3 поврежденного материала к модулю Е3 неповрежденного равно:

Е3 Е3

__________( m i - m 2)__________

(7)

( l + m () ( l + m 2) (У«1 — Ун2)

В случае изотропного материала выражение для %существенно проще и согласно [2]

х=

4(1- v 2)

(8)

 

Ё