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

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

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

 

э

4

8

 

16

 

 

64

128

N \

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

7,85

7,88

7,93

8,07

 

8,23

8,99

 

(157)

(78.5)

(39,3)

(19,6)

(9,8)

(4,9)

 

 

10

 

4,27

4,28

4,30

4,35

 

4,41

4,74

 

(314)

(157)

(78,5)

(39,3)

(19,6)

(9,8)

 

 

20

 

2,87

2,88

2,90

2,90

 

2,96

3,14

 

(628)

(314)

(157)

(78,5)

(39,3)

(19.6)

 

 

Погрешность Лсг< точности нулевого приближения можно оце-

нить по формуле

 

maxlq9— q^0) |

 

 

 

 

 

 

А Ое =

 

 

(6. 10)

 

 

 

 

шах |q0 |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г.0

 

 

 

 

причем

оба

максимума достигаются

в сече­

Р 6ff, /д

 

нии 0 = 0, причем для знаменателя — на внут­

 

 

 

ренней поверхности трубы т—1, а для числи­

ч

----------

N=S

теля — «а

границе раздела

компонентов

8 4

 

внешнего пакета. На рис. 34 приведены гра­

7-

 

 

фики погрешности теории нулевого приближе­

6~

 

 

ния в зависимости от числа ячеек периодич­

5 \

 

N=10

ности (пакетов)

N и степени локализации на­

4 -

 

 

грузки

0,.

 

 

 

 

процен­

--------------_

N=20

Значения этой же погрешности (в

2 ‘

 

 

тах) приведены

в табл. 6.1, причем в скобках

 

 

указано, сколько пакетов укладывается на ду­

/ -

 

 

ге— 0 < 0 < 0 .

 

 

 

 

—I__1 1 1_i i i i i

Когда угол 0 уменьшается

настолько, что

0

5101520253035404?

 

 

J3, град

участок действия локальной

нагрузки

стано­

 

 

 

вится

порядка

двух-трех длин

пакетов, по­

 

Рис.

34.

грешность резко возрастает. Поэтому в таких случаях целесообразно использовать метод выделения особен- но'сгпи.

§ 7. Внутренние напряжения в трубе при ее намотке

Рассмотрим слоистую трубу, получаемую в результате намот­ ки ленты толщиной б, рассматриваемой как слоистый композит, на цилиндрическую оправку радиуса R0 (рис. 35).

Мы предполагаем, что при намотке сделано достаточно боль­ шое число оборотов, так что число образованных при этом паке-

лоненты, причем одна

из

них,

h m 3, в

решении

задачи

не

прини­

мает участия). Продифференцируем (7.1) по параметру

R:

 

dR

\

+

_ и Л _ - * 2 _

_*L =

0.

 

(7.4)

 

 

г

)

г2

dR

 

 

 

Уравнение (7.4)

имеет решение вида

 

 

 

 

 

 

- g - = C 1(i?)r-P + C2(R)rP,

 

 

(7.5)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7-6)

а функции

Ci(R)

и

C2(R)

должны

определяться из

граничных

условий (7.2).

 

 

 

 

 

 

 

 

 

 

В частности, из

(7.2)

следует, что

 

I

= 0 .

Поэтому

 

 

 

 

 

 

 

 

dR \r = R 0

 

 

 

 

C1(R)R0^ + C2(R)RO? = 0.

 

 

(7.7)

Из (III.2)

(приложение

III)

следует, что уравнение

равновесия

для рассматриваемого случая имеет вид

 

 

 

 

 

 

 

 

„ ; +

 

 

= о .

 

 

(7.8)

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

1 II о:

ь°|^

 

•Однако на внешнем радиусе

 

 

J |

___ даг

_

а0

 

dR

 

R •

Из (7.10), (7.5) и закона

Гука:

 

 

(7.9)

(7.10)

------Ь М '

(7.11)

г

следует, что

h [C i(R )R ~ * -'+ C,{R)Rf-'} +

+

[ -

С, (R) R-*~' + С2 (R) /?»-'] =

-

(7.12)

 

 

 

 

 

R

 

Из (7.7) и (7.12)

находим функции C\(R) и C2(R):

 

Q / т

__________________ а°

________________

 

1

 

R ^

(к + рад -

R~» (Х6 -

рх5)

'

(7.13)

с2(Л) = _____________ MR)_______________

+ - К 02р« - р(Х ,-Р ?Ч1)

Тогда решение поставленной задачи по теории эффективного мо­ дуля получается в квадратурах

v(R, /•)= | [С, (R) г-» + Ct (R) г » ] dR.

(7.14)

Г

 

Заметим, что если натяжение ленты постоянно, т. е.

 

Oo(R) = а 0= const,

(7.15)

то интеграл (7.14) вычисляется в элементарных функциях. В са­ мом деле, полагая

rvq0

-,

= _

г-*о0

ЯоР(*б- Р * з)

В = —

Я«Б--

 

 

(7.16)

’^Б +

я;? а.-рад

получаем из (7.14)

 

v (R ,r)= | И ( г ) + В ( г ) ] — R^dR

(7.17)

aR

 

Каждое конкретное число р можно аппроксимировать рациональ­ ным:

(7.18)

п

где т, п — целые числа. Тогда, делая в интеграле

! _

г

fPdR

 

1

x'lfdx

~

J

1 _ а « 2(>

Р

)

замену

 

 

 

 

1

 

 

 

 

 

 

 

II

>

* II

3

приходим к интегралу

 

 

 

 

 

 

j _

« 1/rt

 

,

 

 

?

хт+п-\ ^

 

 

 

J

1ax2n

 

 

 

r\/n

 

 

(7.19)

(7.20)

(7.21)

который берется элементарно.

После этого легко находятся напряжения теории нулевого при­ ближения (7.1):

_________ 1______ v• I <Я/(Я - f - 2ц)) _v___

г

<1/(х + ад>

(1/(х + 2р)>

г

г

(0) =

М£)___________ гУ_

+ [ М|) + 2(i (|)

 

Л / f \

I О..

 

/1 Ч 1 1_

 

 

 

 

 

 

 

Я (6)+2*1 (6)

(l/(X + 2|i)>

 

 

 

 

 

х*а)is;______,I______ M#v уSь)/

W(X+2(Q)

 

— ,

(7.22)

Ж6) + 2ц (g)

 

M i) +2(1(0

(l/(l + 2|i)>

 

Г

 

(0) _

' M i)

 

______*

о +

[M i)

 

 

2

<1/*>\ I

Q.. /fc\

/ 1// “1

 

 

 

M i) +

2(i(i)

<1/(1V(l +

2(1))

 

 

 

 

Xa(l)

 

 

я. (6)

a/(X +

2ii))

I

_u_

 

А-(Б) +2|* (Б)

+

M 6)+2p(6)

(l/(X + 2|i)>

J

г

где Л(£) и |x(£) — постоянные Ламе, являющиеся разрывными функциями | (0 < | < 1 ).

Формулами (7.22) и описываются внутренние напряжения, возникающие в каждом компоненте композиционной ленты при ее намотке на жесткую оправку.

Упражнение 7.1. Показать, что в случае, если эффективный тензор модулей упругости композита является изотропным (мак­ роскопически изотропный материал) и выполняется условие (7.15), то решение задачи (7.1), (7.2) по теории эффективного модуля имеет вид

 

 

 

(2 +

g>)K» + 3o)ftg

.

(7.23)

 

 

 

In—

 

" ' "

Г

°

 

 

 

(2 + <о)/* +3aRl

 

 

 

 

 

 

\ ^

(2+(0)(y + 3aRg

 

 

 

г2

I

(2 + о)/-1 +

3mR§

 

 

 

 

 

 

 

 

 

(7.24)

 

(R, г)

= <Jr(R, г) +

«г,,

 

 

 

 

 

где К — модуль сжатия, <о — безразмерный параметр:

 

 

 

 

 

 

 

 

 

 

 

 

(7.25)

Упражнение

7.2. Показать, что

из

(7.24)

следует,

что

при

и (о<1

на оправке

появляются

сжимающие

 

окружные

напряжения авэ< 0 и что при со = 1

на оправке

всегда

Овэ=ао.

§8. Численное решение пространственных задач

Впоследнее время в вычислительной механике широкое рас­ пространение получили разностные методы и их разновидности (например, метод конечных элементов). Однако непосредственное их применение к решению задач МДТТ для композитов наталки­

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

телом. Ведь каждый

компонент композита должен быть разбит

на достаточное число

«вычислительных»

ячеек, что приводит к

чрезмерно большому

числу алгебраических

уравнений,

решить

которые на современных ЭВМ не удается.

 

 

С помощью метода осреднения в ряде случаев все же можно

преодолеть указанную

трудность. Если

ЭВМ

позволяет

решить

задачу с п «вычислительными» ячейками, то можно это разбие­

ние использовать

для решения

задач

серии

Ж,

т. е. на

 

одной

 

ячейке периодичности

композита,

после

г2

чего

эти

п

«вычислительных»

ячеек мо­

 

жно

использовать

для

решения

 

задач

 

серии Д. Задачи этой серии

являются

 

задачами МДТТ для однородной среды,

 

и поэтому число «вычислительных» яче­

 

ек

никак

не

должно

согласовываться

 

с числом ячеек периодичности и «вычис­

 

лительная»

ячейка

может

содержать

 

достаточно

большое

число

структурных

 

элементов

(ячеек

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

 

 

Дело

обстоит

как

раз

так,

что чем

Рис. 36

труднее

непосредственное

 

применение

разностных методов из-за сложности

 

 

проблемы

дискретизации

(т.

е.

чем

больше структурных элементов содержит композит), тем лучше работает метод осреднения и проще его реализация (в преды­ дущих параграфах было установлено, что при большом числе ячеек периодичности теория! нулевого приближения достаточно точна)ц

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

На примере предыдущих параграфов этой главы было видно, что для слоистых композитов метод осреднения иногда позволяет точно удовлетворить уравнениям равновесия и граничным усло­ виям на поверхностях Гг, эквидистантных слоям. На «слоистых»

поверхностях Гь ограничивающих

композит (т.

е. составленных

из различных компонентов слоистого композита)

граничные усло­

вия могут удовлетворяться только

интегрально

(рис. 36).

Поскольку для слоистых композитов локальные функции пер­ вого уровня и эффективные тензоры модулей упругости и упру­ гой податливости определены, решение задачи Д (0) по теории эффективного модуля, т. е. для однородной анизотропной упру­ гой среды, позволяет построить решение по теории нулевого при­ ближения, т. е. получить микроперемещения или микронапря­ жения.

Для численного решения задачи Д (0) воспользуемся вариаци­ онно-разностным методом, который заключается в том, что вы­ писывается лагранжиан (2.1.8) для задачи в перемещениях или функционал (2.3.3) для задачи в напряжениях в разностном

виде, т. е. производные заменяются разностными, а интегралы в (2.1.8) и (2.3 .3)— соответствующими суммами. Тогда и лагран­ жиан (2.1.8) и функционал (2.3.3) становятся функциями конеч­ ного числа переменных — перемещений или напряжений в узловых точках. Система разностных уравнений получается приравнива­ нием нулю производных этих функций по каждой переменной.

Рассмотрим, например, задачу теории упругости в перемеще­ ниях о слоистом параллелепипеде.

Минимизацией лагранжиана

£ Й = - у j Ч

/fti »<,/»*,! dV— ^X'O'dV —

(8. 1)

V

V

г.

получим решение v по теории эффективного модуля, а по нему найдем микронапряжения теории нулевого приближения:

 

 

 

=

(!)% ./(*).

 

(8.2)

где тензор О 0>определяется по формуле (1.11).

 

Поместив начало

координат в одну из вершин параллелепи­

педа и направив

оси

вдоль

 

ребер, введем

равномерную сетку с

числом

узлов Na на

ребре

 

длиной /а, направленном

вдоль оси

ха, a = lt

2, 3 (хз — направление, перпендикулярное

слоям). Се­

точная

функция

uh(i,

k), i=l,...,A V ,

/= 1,...,М 2;

k=\ , ...,N3

представляет собой дискретный аналог и(х). Заменим частные производные соответствующими разностными производными, а интегралы — конечными суммами. Разностный аналог лагранжи­

ана L(v) обозначим

через

Lh(vh). Тогда

условие

стационарности

лагранжиана

(2.1.9)

6L(v, 6и)=0

перейдет в условия

 

dLh

= 0, £ = 1,

. « а

/ = i.

,

k = 1

N3-,

dva(i, /\ k)

 

 

 

а = 1, 2, 3.

 

 

(8.3)

Поскольку L/l— квадратичная функция от va{i, j, k), условия {8.3) представляют собой разностную систему 3NlN2N3 линейных алгебраических уравнений, которую можно записать в виде

Avh+ fh = 0.

(8.4)

Рассмотрим решениесистемы (8.4) с помощью итерационного метода переменных направлений с чебышевским набором итера­ ционных параметров {тп}

 

B S + i - <

= A j. + fkt

(8.5)

 

Tn+1

 

 

п

 

, Я = 0 , 1,

, N — 1

Jt(2rt— 1)

(8.6)

Ъ +

Yi — (Уг — Yi) cos

2N

(YI и У2 определены в

(8.8)),

где факторизованный

оператор

 

з

 

/ — Ла

0

0

\

(8.7)

В = | | (£ +

^а)*

Ra — |

0

— Ла

О

| ’

а=1

~

\

О

0

— Ла/

 

обращается тремя последовательными прогонками по направле­

ниям Х\,

Х2, Хз. (Ла — одномерный

разностный

оператор Лап­

ласа.)

 

итерационного процесса

(8.5)

решения

Скорость сходимости

системы

(8.4) будет тем

выше, чем

больше отношение

у\/у2 по­

стоянных энергетической эквивалентности операторов —А и В:

 

у1В < - Л < у 2В>

(8.8)

 

~ з

~

~

 

 

Введя регуляризатор

R = ^ R a,

энергетически

эквивалентньГ

 

а=1

 

 

 

 

операторам —А и В с постоянными си с2 и уД у20

 

 

ci R < — A < c2Rt

(8.9)

 

 

 

 

 

( 8. 10)

получим у 1= у 10Сь у22°С2. Известно,

что при

СО 1= 0 2 = 0)3 — 0)

шах У1°/72°= З г]2/3( 1 + 2п) достигается при

 

 

со= а,0= ( У б Д

( У Ь + У Г ) ) - \

 

где

б/А; 6 = min ба; А =

max Да,

 

т] =

 

 

 

а

 

 

 

а ба и Да определяются из неравенств ба / < R a < Да / Ц — тождест­ венный оператор). Поскольку оператор А анизотропный, естест­ венно в данном случае использовать оператор В (8.7) с различ­

ными (Da.

В трехмерном случае наилучший набор о)а неизвестен, но мо­ жет быть найден численно пошаговым поиском от точки o i= ( 02= =<оз = (Оо в направлении возрастания yi°/V2°. Если при этом извест­ ны величины cPa, а 1, 2, 3; 0 = 1, 2, такие, что

з3

 

^

ClaRa ^

А ^ £

°2а 5а,

( 8. 11)

то

а = 1

~

а=1

 

 

 

 

 

3

 

2

 

 

 

Yi — min

с1а Х(Х

 

 

2 с2а Уа

а=1_________ -,

у2 =

max

а = 1 ___________

 

3

 

 

&а<Уа<*а

3

 

П (I +®а*а)

 

П (1 + < W

 

а=1

 

 

 

а=1

 

 

 

 

 

( 8. 12)

и можно непосредственно искать максимум у\/у2 по соа.

188

Выбирая различный масштаб длины по осям хи х2 и х3, мож­ но приблизить анизотропный оператор А к изотропному, для ко­

торого скорость сходимости итерационного процесса (8.5) с опе­ ратором В вида (8.7) выше. Это эквивалентно решению задачи

для тела с соответствующим образом измененными отношениями линейных размеров и модулей упругости (последние становятся близкими к модулям упругости некоторого изотропного тела). Описанные приемы позволяют в случае ярко выраженной анизо­ тропии в несколько раз ускорить сходимость итерационного про­ цесса (8.5).

В случае анизотропного тела в главных осях ортотропии по­

лучим оценку (8.9) с постоянными

 

 

 

с, =

min (ftull — - lm Лп

 

hji22^223!

 

 

^2233

 

1Лпзз

 

fa

^1133 ^2233

2К,

2^2323j>

(8.13>

 

 

С2

гпах (hnll + hu22 + Л1133, ^ 2 2 2 2

+ ^ 1 1 2 2 +

 

+ ^2233» ^3333 “Ь ^1133 “Н^2233* ^1212* 2Лl3ig, 2/ta323).

(8*14}

Итерационный процесс (8.5) прекращается, как только удовлетворяется неравенство

IK+.-41I

(8.15)

И - 4 1 1

 

где под II •II понимается длина вектора.

Рассмотрим численный пример. Пусть на гранях х3 = 0, х3=1 заданы нагрузки, а на гранях *i = 0 ; *i = 8 ; *2= 0; *2 = 8 — усло­ вия жесткого защемления (рис. 37)

S?|*,:=о = 0, / = 1,2; = - о о а ~ ^ х

s?U=, = o, 1=1,2, з,

«°1 = 0 при х, = 0,

хг = 8; хг = 0 ,

хг = 8.

(8.16)

Пусть ячейка периодичности

укладывается

на отрезке

0 < х 3<1

целое, достаточно большое, число раз. Предположим, что компо­

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

компонент —

однородный

и изотропный с коэффициентами Пуассона vi и v2 и

отношением

модулей Юнга к= Е2/Еi и объемной

концентрацией

у= 1/2.

Счет

каждого варианта

на ЭВМ БЭСМ-6 для N\=N2= 9,

N3=5 (т.

е. 1215 разностных

уравнений) составлял 5—8 минут.

По результатам счета были получены напряжения теории ну­ левого приближения (8.2).

На рис. 39—41 представлены графики распределения микро­ напряжений, отнесенных к величине q (8.16) для каждого из компонентов (1) и (2) по высоте параллелепипеда х3 (рис. 37).

Ввиду разного масштаба, принятого для различных комподентов, все графики изображены непрерывными линиями. В гра­

фиках, соответствующих компоненту (1), нужно исключить все отрезки, принадлежащие компоненту (2), а в графиках, соответ­ ствующих компоненту (2),— исключить отрезки, принадлежащие компоненту (1).

Было рассмотрено несколько случаев с различными отноше­ ниями модулей Юнга компонентов к и коэффициентами Пуассо­ на vi, V2- Поскольку напряженно-деформированное состояние об­

ладает очевидной

симметрией, достаточно

 

рассмотреть область

0 < * i< 4 , x i< x 2<4, 0 < ^ з< 1 (рис.

38).

 

 

 

 

 

 

Расчеты показывают, что в случае %=1 максимальные растя­

гивающие

напряжения итп возникают

у

обоих

компонентов

в

окрестности точки Х\(0, 4, 0).

 

 

 

Т а б л и ц а 8.1 (х =

1)

 

 

 

 

 

 

 

 

 

 

Компонент (1)

 

 

 

Компонент (2)

 

 

 

напряжения

 

 

 

 

напряжения

 

Vi

нормальные

касательные

v2

нормальные

касательные

 

 

 

 

0,01

aiV Й ) = 14,00 <*13 Й ) =<*1зЙ )=

0,25

а<2> Й ) =

15,80 а(12з) Й )= с т (123)( й =

 

 

 

=

6,00

 

 

 

 

= 6,00

 

0,01

о '1,» Й

= 9 ,4 2

oft* (*Т) = 6,74

0,49

а(2> Й ) =

20,50

а(2) (х2) = 6,74

0,25

a ft Й )

= 16,50 <*13 Й ) =

а 11з Й ) =

0,49

o ff Й

) =

28,90 а(2) Й ) = o f f f c ) -

 

 

 

=

5,45

 

 

 

 

= 5,45