книги / Микромеханика композиционных материалов
..pdfПотребуем, чтобы на участке г — ае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) |
|
|
|
|
|
|
|
|
||
|
1Г 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(п2— 2п — 1 + 2\) — ° п("п+3 1) |
| Рп (соз0) — |
|
||||||||||||||
- 2^ 1 ^ т ( - « + |
4 - ^ |
) + |
^ з |
й Р п (соз 0) |
|
|
|
||||||||||
|
|
я |
— |
с{§е' |
|
||||||||||||
О®= |
2С[ - ^ (/г + 3 — 4лл> — 2л>) - |
А , (я + |
1) |
Рп (СОЗ 0) + |
|
||||||||||||
+ |
20 |
|
|
( - |
п + |
4 - |
4У) + |
- % ,] ' Ч |
|
Г |
|
°- |
|
е - |
<6-15) |
||
Условия идеального контакта компонентов при |
г = |
|
а |
|
|
|
|||||||||||
|
|
“ в = |
“ в* |
“? = |
“ г. |
°г = |
а г. |
Ог“е = |
|
0 ,в. |
|
|
(6.16) |
Эти условия дополняются заданием напряжений или деформаций, характеризующих однородное взаимодействие включений.
100