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

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

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

Потребуем, чтобы на участке г — ае1Ь% функция имела

аналитическое продолжение, что приводит к функциональному уравне­ нию

2 (Хтла"еый - Х тпа ^ ш ) = 0.

л > 0

Сравнение этого уравнения и (5.37) приводит к определению связи между коэффициентами разложения функций

 

^

= V [ 1 +

(-1. ,2Л*] а"»-»-,кпРта.

(5.39)

Краевое условие для

вспомогательной функции

 

 

откуда

х% (т)-

(т) =

0,

 

 

х = да'®,

(5.40)

 

 

 

 

 

 

 

 

 

 

 

* и ( г ) = * ( 4 )

 

 

<5 -«)

Удовлетворяя второму краевому условию (5.34), приходим к функ­

циональному уравнению

 

 

 

 

 

2 -

а" (<?тпе‘п*-С1тпе-ш ) + сГ (Стеш

- Сте~ш ) = 0,

 

Л>0

 

 

 

 

 

 

 

 

 

 

 

 

Ъа< # < К

 

 

(5-42)

Введем вторую вспомогательную функцию

 

 

 

 

 

 

л>0

 

( 4 ) “ ,

|г |< о ,

 

У т Ы -

 

 

 

 

 

(5.43)

У т ( г ) = 2 ? тп (-Т -)" + Ст Л

I г | > а,

 

 

 

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

 

 

 

Ут (т) — Ут (т) = 0,

 

„ « > « > » .

(5.44)

Отсюда следует

 

 

 

 

 

 

 

 

Ут(*)= У т(-Т-)->

Утп = 4тп-

(5.45)

Вторая

формула

(5.35) — равенство

напряжений на межфазной

границе — приводит к функциональному соотношению

 

[Хга(т )+ аУт (т)]+ - [ Х в (т) + аК- (х)Г -

0,

* , < * < * „ .

(5.46)

Первая формула (5.35) — равенство смещений на участке совершен­ ного контакта — эквивалентна функциональному уравнению для вспомогательных функций

91

[Хга (Т) -

т(Т)]+ + [Хш(т) - „У. (Т)Г = 2 ^ 1 ^ 2

а~"Ы +Ч '">'Х

 

уу /п

П—Ып п

,

п _Л«т

с2л \

(5.47)

 

х

т

+

Ртпа

— ] .

Для определения вспомогательных функций имеем систему функ­ циональных уравнений (5.40), (5.44), (5.46) и (5.47). Соотношение (5.46) определяет аналитическое продолжение комбинации вспомога­ тельных функций в интервале разреза; удовлетворяя последнему пре­ дельному условию (5.35), а также учитывая равенства (5.41) и (5.45), находим

Хт (г) + аУ т (г) = аРт(г),

(5.48)

где

 

Рт(2) = Стгт + Ст^ т

(5.49)

Равенство (5.48) позволяет найти

 

Ут(г) = Рт( г ) - ± Х т(г).

(5.50)

Подставляем (5.50) в неоднородное уравнение (5.47) и сводим определе­ ние Х т (г) к решению неоднородной задачи сопряжения [531, которую решаем методом Карлемана [541:

где

 

^ ( т )

+ Х-(т) = Фт(т),

 

 

(5.51)

2Vа

г.

,_у

,

2ча

 

Р — V

 

 

 

 

 

 

 

 

х

 

 

ч’шМ -

 

 

 

 

 

 

2

 

 

 

X

(Ртпа

- ^ п +

Ртпа‘*Л- ^ )

 

 

(5.52)

Решение задачи (5.51) известно [53]:

 

 

 

 

 

 

 

 

(г) =

И <*)Р т (г) +

 

} ^

С

- г

) '

 

(5'53)

Здесь I — контур межфазной

границы

с

совершенным

контактом:

 

 

 

 

 

 

 

 

 

 

 

| г , <

| г »1-

(а(г) =

V (г — ае‘*а) (г — ае‘0ь)

= |

 

Л“ °

г

,

 

 

 

 

 

 

 

 

 

 

 

 

 

Iг I >

I 21

 

 

 

 

 

[

л=0

 

 

 

 

•о |»

 

 

 

 

 

 

 

 

 

(5.54)

 

 

 

г0= ае1\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*0 = 1,

К = — со5

,

К =

Л» (соз 4-) “

2с08 4 "Рл-» (С05 4 ) +

 

 

Рл-р2 ^соз ”2“^

,

ц ^

2;

 

 

(5.55)

92

Рп(соз -у ) — полином Лежандра первого рода л-го порядка [31]; 60=

=а ^ — угловая координата середины участка с совершенным

контактом; Ф=

— величина угла,

охватываемого

этим участ-

ком (см. рис.

44);

 

 

 

Кт й

= А> +

+ Л*2 + ••• Н—

Н---- ^г- +

•••

Для вычисления интеграла рассмотрим интеграл по замкнутому контуру, охватывающему дугу /. Согласно теории вычетов [29], имеем

<Рт(Т>‘*Т

=

Ке$ (0) +

Кез (оо) + Кез (г),

1 т <г> = -2 ^Г $ - Р(Т)(Х —2)

 

 

 

где Кез (г) — вычет подынтегральной функции в точке г. Для опре­ деления вычета Нез(0) разлагаем функцию в ряд по степеням г, со­ храняя в разложении только главные члены:

 

%»<т>

 

 

2лкх-

1 /1 Г

Т

I

х3 .

 

 

 

V

 

р(т) (т — г)

 

 

 

+

т

+

_? +

 

 

 

 

 

п> 0

 

 

 

 

 

 

а2п

(

Р^соз 4 " )

 

 

 

 

 

 

+

 

 

X

т"

' 1 + -

 

 

 

 

 

 

 

Здесь

использовано разложение

 

 

 

 

 

р ( 2 )

 

л=0

 

 

 

(5.56)

 

 

 

 

 

 

 

 

2 Рп (с05т) (т-)"-1

N>1*1

 

 

Совершая предельный переход т->~0, получаем

 

 

-

^

[

г .

У

 

 

+

 

 

 

 

 

У=1

 

 

 

+

- ^ - 2 2

ггп',а -п,-ч- " +шРяпРп- 1(сое 4 -) - З Г ( т - / 1

 

я>0 /=1

 

 

 

 

 

 

 

Вычет Кез (оо) определяется главным членом с обратным знаком в разложении подынтегральной функции в ряд по большим т->оо:

93

[с -> - + - № ■ I ,

+

+ ^тпЛ‘Х" ^ г ) (I + ~ + 7Г + •••) ! ^ - +

Р4 (с»»-у) +

откуда

 

 

 

 

+

^ р г (с°з4 ) + • • • ) ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ке* (00) =

-

^

 

^

-

(с т 2 * Р;_, (соз ± )

( 4 ) М гпЖ +

 

 

о

 

 

Г

 

 

 

т—2

 

 

 

+ С ^ 6 1 т + ^ -

 

„ - ^ - О + Ь . , ^ + 2 2 а - ^ + ' + ^ Р ^ х

 

 

 

 

 

1

 

 

г»»>0 /=1

 

 

 

 

 

 

 

 

 

*(«4 )(4 ) '^ Ц -

 

 

Стягивая контур интегрирования к берегам

границы

(см. рис. 44),

получаем

искомое значение

интеграла

 

 

 

 

 

1

Г

<Рг о

( т

М

т

 

 

у т ( г )

т

I

 

о * .

 

и+ (т)<т - 2)

 

 

2ц (г)

у + а

Г 1

г2

+

 

т + 2

 

 

 

 

 

 

 

 

 

 

 

 

+ С ,

 

 

( с - I ) ( А ) м ^

+ С » ^ Р „ _ , ( с , * ) X

 

 

 

 

 

 

 

 

 

 

/=1

 

 

 

 

х (-%■)'-$- +

4

> — V |„ - У + 1 + < Х е

 

 

 

 

^ - [ а -

 

 

Ли! +

 

 

 

I

г

/

г0

Т

р + у [

 

 

 

+

2 2

 

 

 

 

 

(соз 4 ) ( А

/ 4

рг +

 

л>0 /*=1

 

 

 

 

 

4

7

4

7

О

+

2 2

а-л<1,+1+й0е2л',/>/_ , /С03 4

) ( 4 ) ' ^

Р

тл]| . (5.57)

 

л>0/=1

 

 

 

 

 

'

/ \

/

 

1)

Здесь 61т — символ Кронекера; р+ (х) +

р“ (т) =

0.

 

 

Заменяем интеграл (5.53) функцией (5.57) и устанавливаем явную зависимость вспомогательной функции Хт (г) от параметров. Расклады­ ваем Х т (г) в ряд по степеням г и приравниваем коэффициенты при оди­ наковых степенях г согласно формуле (5.31). Это позволяет определить коэффициенты Ртп в явном виде и суммировать ряды при значениях аргумента г *= ае™ Если контур области образован кривой, описыва­ емой ортогональными криволинейными координатами, то решение строится путем конформного отображения рассматриваемой области на круг [53].

А

Г Л А В А 6

СРЕДЫ СО СФЕРИЧЕСКИМИ ВКЛЮЧЕНИЯМИ

Однородная матрица со сферическими включениями — простейшая модель дисперсного или обычного упрочненного композиционного ма­ териала. Если размеры включений не превышают доли микрона, то имеем дисперсно-упрочненный материал, а если превышают указанную величину, имеем материал с включениями, например твердосплавные композиции на основе металлов, керметы, строительные материалы, сферопластики и др.

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

§ 1. СТРУКТУРА СРЕДЫ

Реальные материалы имеют нерегулярное размещение частиц в пространстве, различные формы и размеры. Модельные структуры бу­ дем строить исходя из трансляционной и поворотной симметрии регу­ лярной среды с одинаковыми частицами. Распределение центров час­ тиц в пространстве представим в виде последовательности векторов пространственной решетки [67]. Простейшую структуру подобного типа образует бесконечный набор параллелепипедов, угловые точки которых определяются вектором

е = а1тсо1 4- егпю2+

е3ра3,

где ек (к = 1, 2, 3) — единичные векторы,

направленные вдоль сто­

рон параллелепипеда; т , п, р — набор натуральных чисел; сой = = 1, 2, 3) — длины сторон параллелепипеда периодов. В более слож­ ных структурах около каждой угловой точки пространственной решет­ ки находится группа включений, что приводит к различным струк­ турам и группам симметрии [63].

Для построения теории эффективных модулей в первом приближе­ нии первостепенное значение имеют определение относительного объем­ ного наполнения среды частицами С, вид упаковки материала сфера­ ми и циклическая симметрия в его структуре. Первый и второй пара­ метры определяют в регулярной структуре взаимное расстояние меж-

95

Рис. 45

Рис. 46

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

Наиболее простая кубическая решетка = о>3 = ш3 = ш пред­ ставляет собой куб (рис. 45), в углах которого размещены центры одинаковых шаров. Восьмая часть объема каждого шара в узле сетки входит в объем одной ячейки, поэтому на одну ячейку приходится всего один шар. При плотной упаковке шаров радиуса а предельное объем­ ное заполнение достигается, когда сферы касаются друг друга м = = 2а:

(6. 1)

Напряженное состояние в такой структуре будет обладать циклической симметрией в трех взаимно перпендикулярных плоскостях с углом периода а = я/4.

Если в кубическую структуру поместить еще одну сферу в центре куба (рис. 46), так что в одной ячейке будет уже два шара, то плотная упаковка достигается при касании шаров по диагоналям куба. Пре­

дельные размеры определяются из уравнения

 

(4а)2 = со2 + (со УЪ)\ а =

;

предельное объемное содержание

 

 

(6-2)

96

Рис. 47 Рис. 48

Такая структура называется объемно-центрированной кубической. Циклическая симметрия напряженного состояния в этой решетке бу­ дет определяться и наклоном плоскости симметрии, соответственно изменится угол периода.

Если на каждой грани куба в кубической структуре поместить до­ полнительно по одному шару, половина объема которых войдет в дан­ ную ячейку (рис. 47), то в выделенном объеме будет всего четыре вклю­ чения. Предельная упаковка достигается при условии

(4а)2 = 2(о2, а = <о .

Допустимое наибольшее объемное содержаниз включений

(6.3)

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

В гексагональной плотной структуре основание ячейки образова­ но гексагональной сеткой; центры сфер в следующей плоскости сдви­ нуты так, что шары расположены между лежащими в нижней плос­ кости. Каждый шар имеет шесть ближайших соседей в основной плос­ кости ячейки (как в плоской структуре), другие шесть соседей рас­ положены по три — выше и ниже основной плоскости (рис. 48). Пре­ дельная степень упаковки материала достигается при касании шаров друг друга:

(6.4)

97

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

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

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

Введем систему сферических координат г, 0, ср естественно свя­ занной с межфазной границей среды (рис. 49)

XI = Г СОЗ 0,

хг =

Г 51П 0 СОЗ ф, Х3 =

Г ЗШ 6ф,

0 ^ г < о о ,

0 ^

0 ^ тс, 0 < ф <

(6.5)

2л.

Разложение ортогональных единичных орт на сфере еГл еэ и ёф, направленных в сторону возрастания сферических координат, имеет вид

ег = е1соз 0 + е2з т 0 соз ф + е3з т 0 з т ф,

ее =

е1з т 0 4- е2сов 0 соз ф +

е3соз 0 з т ф,

(6.6)

 

вф = — е2з т ф +

еэ соз ф.

 

Первая квадратичная форма и площадь

элемента на поверхности сфе­

ры

йгг +г*4Ф + г2$1П2 б^ф2»

= г2з т 040<йр.

 

^52=

 

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

98

ные свойства не зависят от выбранного направления, поэтому в при­ нятом приближении направим ось хх так, чтобы напряженное^состояние не зависело от угла ср. Деформированное осесимметричное состоя­ ние в сферической системе координат определяется компонентами [47]

 

 

да

 

1

диа

ш

I

 

 

г ,~ И Г '

в0 “ Т 'э Г

+

т и,■

(6.7)

 

 

 

 

 

 

да

дил

1

, Л

1

 

Уге =

1

1

8<р = — и&с% в +

иг,

у

 

+ у г

 

 

 

Д =

ег + е э + еф,

 

 

где иг, ив, //ф — компоненты

вектора

перемещения

 

 

и = ёгиг +

евив + е9и<р.

 

(6.8)

Компоненты тензора напряжений в сферической системе координат

Т = егёгот+ е0ёеа0 + ефеф(Хф+ (егев + евег)аг0.

(6.9)

§ 2. ОСЕСИММЕТРИЧНОЕ СОСТОЯНИЕ

 

Задача теории упругости об осесимметричном напряженно-дефор­

мированном состоянии сферы рассматривалась рядом авторов [50, 72], поэтому для удобств в дальнейшем без доказательств выписыва­

ются основные результаты [50].

 

возрастающие с удалением'.от

 

Компоненты вектора перемещения,

начала координат, имеют вид

 

 

 

 

 

 

= [Апгп+1(гс +

1) (п — 2 +

4у) + Вппгп~х\Рп (соз 0),

(6. 10)

 

 

 

 

 

 

 

(соз е)

 

 

«8° = 1Алг"+1(п +

5 — 4V) +

 

 

 

 

Впгп~1] —у в -----»

 

объемное

расширение

 

 

 

 

 

 

 

 

Д(л) = — 2 (1 — 2V) (2п + 3) (п + 1) Апг*Рп (соз 0),

 

где

Ап, Вп— постоянные; Рп (соз В) — полином Лежандра первого ро­

да

/г-го порядка, я =

0, 1, 2, ...

Общее решение

с этими компонен­

тами

ОО

 

00

00

 

 

 

 

 

 

 

 

 

и г = 2

«(?>.

 

= 2

Д = 2

д".

(б. 11)

 

 

л=0

 

п=1

л=0

 

 

Соответствующие компоненты напряженного состояния

 

о<*> = 20

(п + 1) (п2— п — 2 — 2V) гп + Влп (п— 1) г " -1]Рп{сов 0),

 

о'"е) =

2 0 [Лп(п2 +

2 л -

1 +

2V) т" +

Вп(п — 1 ) г - ’]

^ ,

99

■о*"' = — 2О [Л„ (п2+ 4л + 2

+ 2V) (п + 1) гп + Впп2гл~2] Рп (соз 0) —

— 20 [Ап (п + 5 — 4\) г" + Впгп~2] — л^ -°3 6) с!§ 0,

(6.12)

о<”>= 20 [Ап (п + 1) (п — 2

— 2V — 4пу)гп + Вппгп~2] Рп (соз 0) +

 

йР

(соз 0)

0.

 

+ 20 [Ап (п + 5 — 4у) г* + Впг*-2} -

п^ —

 

Общее напряженное состояние получается суммированием функций (6.12) по всем п. Компоненты вектора перемещения, убывающие .’ с удалением от начала координат, будут

 

 

=

^ / 1 ( п +

3 -

4 ? ) - - ° " Д +

‘-| л ,

(соз 0),

(6.13)

 

 

 

ГС

 

,

 

 

Д

1 арп

(соз 0)

 

 

 

и!*'

 

 

 

(С05

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

объемное

расширение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д,п’ => — 2 (1 — 2V) (2п -

1) пСпг-«+»Рл (соз 0).

(6.14)

Напряженное состояние

 

 

 

 

 

 

 

 

 

 

 

 

сг = 20

 

С п

 

п

п

%

Д , ( я + 1 ) ( я + 2)

 

Рп (С05 0),

- ^ т (я» +

З я - ^ ) +

"

^

 

 

------ -

 

аГ0 =

20

с

 

л

* +

Я * )

Д, (я + 2) 1

й Р п (соз 0)

 

п + 1

\ П

 

+3

-

I

 

п

0

 

 

 

 

 

 

п—

/и2_0-1- 9 ^ ____2—__!

 

 

 

 

 

= 20 - ^

1(п22п — 1 + 2\) — ° п("п+3 1)

| Рп (соз0) —

 

- 2^ 1 ^ т ( - « +

4 - ^

) +

^ з

й Р п (соз 0)

 

 

 

 

 

я

с{§е'

 

О®=

2С[ - ^ (/г + 3 — 4лл> — 2л>) -

А , (я +

1)

Рп (СОЗ 0) +

 

+

20

 

 

( -

п +

4 -

4У) +

- % ,] ' Ч

 

Г

 

°-

 

е -

<6-15)

Условия идеального контакта компонентов при

г =

 

а

 

 

 

 

 

“ в =

“ в*

“? =

“ г.

°г =

а г.

Ог“е =

 

0 ,в.

 

 

(6.16)

Эти условия дополняются заданием напряжений или деформаций, характеризующих однородное взаимодействие включений.

100