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

книги / Периодические кусочно-однородные упругие структуры

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

§6. Некоторые приложепвя теории гармонического регулярного поля

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

описанную в § 1 структуру.

Если материал каждого компонента структуры имеет плос­ кость упругой симметрии, перпендикулярную оси волокна, то паиряженпос состояние среды при действии средних напряже­ ний <Пц>, <0i2>, <022>, <Oi3>, <02з> распадается на два незави­ симых двумерных состояния — плоской деформации и продоль­

ного сдвига.

В атом последнем случае смещения Щ = U2= 0, из = щ (хи х2);

«отличпы от пуля лишь

сдвиговые

напряжения

On = 013(^1, х2)

Ж 023 — 023 ( X i , Х2) .

 

 

 

 

 

 

 

 

 

Имеют место [3]:

 

 

 

 

 

 

 

 

 

уравнения равновесия

 

 

 

 

 

 

 

( 6.1)

соотпошения Коши

9iOi3 Ч- 62O23= О,

 

 

 

 

 

 

 

 

 

 

 

(6.2)

С\з =

d\U3,

623=

92из,

 

 

 

закон Гука о — е

 

 

 

 

 

 

 

 

 

023 = С44С23 +

045613,

С44 >

О,

С55 >

О,

(6.3)

013 =

054623 +

С55613,

С45 =

С54

 

 

 

либо закоп Гука е — о

 

 

 

 

 

 

 

 

 

023 = 044023 +

045013,

®13 =

054023 +

O5SO13,

(6.4)

причем

 

 

J

 

 

 

 

Cjg

 

Cg

 

 

 

 

 

 

 

Д (о) =

о.ио35 -

<4 = - Щ - >

0.

 

 

Подставляя напряжения из

(6.3)

в

уравнение

равповесня

(6.1), приходим с учетом (6.2) к однородному дифференциаль­ ному уравнению второго порядка относительно из

c6s^iu3 + 2c4fidid2u3 + cu dlu3 = О, Д (с) =

с44с55 — с\ъ > 0.

(6.5)

Решение этого уравнения имеет вид

 

 

Us(xu х2) = Recp(z),

 

(6.6)

тде

»У л (с) —с4

 

z = x 1 + рх2,

 

ц

- « а + гр.

 

Здесь ф (z) — произвольная аналитическая функция переменной z.

181

Вводя

выражение для смещения

из из (6.6) в равенства

(6.2) и

(6.4), получаем

 

 

 

 

<?t3 = R e® (z),

е2з =

Res[ji<D(z)],

 

 

0i3 = УД (с)1ш [цФ (г)],

023 = ~ VA (с)1ш Ф (г),

(6.7)

 

Ф (*)-

*Р(«)

 

 

dz

 

 

Сдвиговые напряжения на площадке с нормалью п (cosi|), sinvf). очевидно, имеет вид

0„з = Oi3 cos г|) + 023 sin t|) = VД (с) 1щ [а (ф) Ф (z) ],

( 6.8)

а (ф) = ц cos ф — sin ф.

 

Переходя теперь к описапнон в § 1 регулярной структуре» потребуем, чтобы в ней имели место средние сдвигающие напря­ жения <01з») и <02з>. Это приводит к условиям

j Опзd s ------ ш^аад),.

[

0»,з ds = <0i3> Н — <023> Л. (6-9)

где h = Re ©г, # = Im сог-

вектора

0 = 013+1023

через произ­

С другой стороны,

поток

вольную дугу АВ определяется, в силу

(6.8), формулой

 

Q =

! orn3ds =

pcj., Im ср (z0) U .

(6.10)

 

А В

 

 

1 0

 

Сравнивая (6.9) и (6.10), находим

 

 

Im {<р (z0 +

£Ojo) — ф (z0)} = —

 

<023> = Im Qlf

(6.11)

 

 

 

 

 

 

Im {ф(z0 +

<в20) — ф (z0)} =

{<оаз> — <а23> ctg а} =

Im Й2.

Наконец, условия сопряжения компонентов среды по пере­ мещениям из и усилиям 0„з, очевидно, совпадают с краевыми условиями (1.5), в которых необходимо положить gh(t) = 0 ,

gkm(t) = 0 1).

Таким образом, продольный сдвиг многокомпонентной сре­ ды, поперечное сечение которой имеет описанную в § 1 регу­ лярную структуру (при условии, что в среде имеют место сред­ ние напряжения (013) и <02з^), полностью описывается теорией регулярного гармонического поля.

') Предполагается, что между компонентами среды соблюдаются усло­ вия идеального контакта.

182

При этом необходимо иметь и виду следующие формулы пе­ рехода:

(Jl = 013, ?2 = 023,

<9l) = <°13),(qs) — <И2з),

К ц — ^55, ^12 = Сц5К. 22-—<?44, ТУ = Ц3,

Xl 1 = Д55, Xi2 = 045,

(6.12)

Х22~ «44,д\W = С|3,

а2ТУ = £23,

q = о,

(/.„ = ОпЗ-

Аналитические функции, описывающие перемещение ц3 в раз­ личных компонентах среды, задаются интегральными представ­ лениями (1.12). Фигурирующие в пих плотности определяются

из системы

(2.1). Уравнения состояния макромоделп имеют вид

 

 

(егз) =

(044X 025) +

(045) <0i3),

(6.13)

 

 

*<Cl3> =

<aS4X<J23i) +

<055><0|3>,

 

 

 

где согласно

(5.6)

 

 

 

 

<»„> =

« ,,(

i - Ц - R e i )

<««> = «45 — X ®44 R e -i,

<0S.|> = «-If, -J,F «55 Re fl2-

<Q55> = «56 (1 —

Re ®lj-

Стандартные

функционалы ah

= 1, 2) определяются соот­

ношениями

(5.4)

и (5.5)

с учетом формул перехода

(6.12), ха­

рактеристическое число Цо вычисляется в (6.5), где с*

(г, & = 4, 5)

отпосятся к матрице.

Результаты расчетов контактных сдвиговых напряжении и макропарамстров <aih> для армированных сред с изотропными компонентами по несколько более простым алгоритмам приведе­ ны в гл. 6.

Поперечная теплопроводность армированных сред (КМ). Пусть анизотропная среда, армированная копгруэптнымп группа­ ми ппородиых анизотропных волокон (которые в свою очередь могут быть армированы инородными волокпамп, так чтобы поперечпое сечеппе такой среды пмело описанную в § 1 структуру), пронизывается нормальным к оси волокна стационарным тепло­ вым потоком со средними значениями в каждой ячейке <qi> и

<q2> (рис. 5.1.1).

из компонентов

структуры

Уравнения состояния в каждом

имеют ппд [2]

Яи>0, A-Z2 > 0 ,

. . . . .

qi = —X u d iT — Х\2д2Т,

 

о

(6.14:)

q2 — — ^'21^1^' — ^22^2^i A M ~ ^п^аг ^та ^

где h j — тензор теплопроводности, причем в силу принципа сим­ метрии кинетических коэффициентов Х\2 = Х2\\ Т = Т (хи х2) — температура. Очевидно, в каждом из компонентов структуры принимают свои постоянные значения, а температура описывает­ ся различными функциями {Т-, Тк, Т ^ ).

183

Обращение закопа

(6.14) дает

 

 

 

 

 

(6.15>

 

01Т = -tu ffi -

41292,

д2Т = -4 2 I(ZI -

422(72»

где

А,,„

 

 

-I

 

 

_

л.

 

 

 

 

^12

 

 

л»

 

Y u

= Т ( X T ’

Y l2

_ Тз1 ”

Д (W ’

722

 

д (W

 

Д (4’)

= Y iiY a e -

Т м = д Щ

>

° -

 

 

 

Условие

div q = О

приводит к

уравнению

(1.4)

относительно-

температур и, следовательно, к формулам (1.7).

Потребовав существование в среде средних потоков <qi>, <g2>, приходим к дополнительным условиям (1.10). Наконец, усло­ вия сопряжения компонентов среды по. температурам и потокам приводят к краевой#задаче (1.6).

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

К\\ = С$5 =

—А,ц,

К\2 =

С45 =

А,12,

К-71 — с44 =

—^-22;

W = из = Т,

(71 =

013,

(72 = 023,

(6.16).

<01^> =

<eia> =

<diT>,

<d2W> = <е23> = <д2Т>,

V-ll = а 55 =

—4u,

Я12 =

045 =

—412,

Я22 = «44 =

—422,

 

 

(7,и =

о „ 3.

 

 

Результаты расчетов осреднеппых теплофизических характе­ ристик армированной среды представлены па рис. 5.6.1—5.6.4. Материал матрицы и включений изотропен в смысле теплофизи­

ческих свойств, причем А»11А п = 20 (А,пД п — коэффициенты теплопроводности матрицы и волокон соответственно). В преде­ лах ячейки имеется одно волокно эллиптического поперечного сечения с полуосями Ri и R 2, ориентированными вдоль осей

Х\ И %%•

На рис. 5.6.1, 5.6.2 покаэапо изменение величин 4п/<4п> и

4гг/<422> для квадратной

решетки

(сог =

icoi)

в

зависимости

от

параметров A, = 2JRI/COI и

А,* = А,//?,. Изменение тех же

вели­

чин для прямоугольной

решетки

(002 = 0,5tcoi)

представлепо

на

рис. 5.6.3, 5.6.4.

 

 

 

 

 

 

 

Поперечное электростатическое поле

в

армпровапных

ди­

электриках. Пусть армированный

диэлектрик,

поперечное

сече­

ние которого имеет описанную в § 1 структуру, пронизывается перпендикулярным оси волокна стационарным электрическим

полем с одинаковыми ■в каждой

ячейке

средними

значениями

электрической индукции

<Di> и

W 2>. Уравнения

состояния в

этом случае имеют вид [5]

 

 

 

 

D\ = ец/?1 + e i2E 2,

D2 =

&2IE I ~\-е22Е 2, 612=

621,

(6.17)

где Е ч и Dv (v = 1, 2) — соответственно

компоненты

вектора

напряженности электрического поля Е и

электрической

индук-

184

7///W

Угг/^Ъг >

Рис. 5.0.1. Зависимость величины 1fn/<lfii> Для тстрагоиальиой решет­ ки (ю2 = (ьм) с подокнами эллипти­

ческого сечения от параметра X = = 2«,/о), для различных значении

А* = А ,/А ,. Кружочками показаны

результаты работы [8], крестикалш— [10], квадратиками в кру­ жочках— [9]

Рис. 5.6.2. Изменение величины Тгг/<Т22> для тетрагональной решет­ ки в зависимости от \ = 2Лг/о>1 н

K = R 2IR,

Рис. 5.G.4. Иэмепеппе ЧнКчпУ для прямоугольной решетки в зависимо­ сти от %= 2Д|/(1)1 нЯ,* = Л /Л

185

ции D, a

eik (i, k = 1,

2) — диэлектрические пропицаемости сре­

ды. Обращение закона

(6.17) дает

 

 

 

 

 

 

£ i = бц1?1 + 6i2-D2)

Ez = 82iDi + 622D2,

 

(6.18)

причем б«, связаны с величинами e ih формулами тина (1

.3).

Вводя

потенциал поля / по

формуле

Е = grad /,

приводим

уравнение

Максвелла

div D = О

к

дифференциальному

уравне­

нию (1.4)

относительно потенциала

/(ад,

яг). Поэтому

для по­

тенциала поля в каждом из компонентов структуры справедли­ вы формулы (1.7).

Характеристическое число |х определяется формулой iУ Ш —е

ц =

----------------—,

А(е) — епем — е*а > 0,

(6.19)

 

®22

 

 

Нормальная

компонента вектора D =D\ + iD2 на площадке

с нормалью тг{созф, sin. имеет в силу (6.19) и (1.7) вид

 

Dn = Di cos ф + D2sin ф = A(e) Im [а(ф) Ф (z) ],

(6.20)

где

 

 

 

a (ф) = p, cos-ф—

Ф (z) = d(^ .

 

Поэтому поток вектора индукции через дугу А В имеет вид

(1.9),

а условия квазипериодичности ср (яо) — (1.10).

Электрические краевые условия на границах раздела диэлект­

рических материалов таковы:

 

 

 

 

 

 

 

 

E t = Е 7,

D$ -

Dn = 4яд,

 

(6.21)

где

Е , — касательная

составляющая

вектора

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

электрического

ноля,

q — плотность

поверхностных

зарядов па

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

 

 

 

 

 

 

В силу (1.3) и (1.9), эти краевые условия сводятся к равенст­

вам

(1.11). Средпие значения

вектора

D имеют

в

соответствии

с (1.13)

вид

 

 

 

 

 

г

 

 

 

*+<о2

 

2

 

 

 

 

< °>>

=

Т Г J

+

[ А.*,

<А> = 4- I

(°-22>

 

 

*

 

1 Z+ttj

 

12+<01

Значения потенциала поля

(7, U; fhm) определяются представ­

лениями (1.7),

(1.12), причем входящие туда плотности вычисля­

ются из системы (2.1).

 

 

 

 

 

 

Средние значения вектора напряженности электрического по­ ля задаем формулами

<а > - <«,/>=

(6.23)

86

Таким образом, поперечное электростатическое поле в арми­ рованных диэлектриках описывается изложенной выше теорией

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

<#1> = <бП><^1> + <6 i2><I>2>,

 

<Е2У= <б2.XZ),) + <622><D2>,

(6-24)

где

 

 

< 6 , 1 - -7- Ве

^ вц Не

 

<62i) — б12— -762,

/’ = ш ,1т ю 2,

 

< * „ > - 8 « [ l - £ R » ( i ) ] .

Формулы перехода, например, к продольному сдппгу армировапной среды таковы:

С55 = еп,

C45 = ei2,

С44 = £32»

 

fl'55 = 611, а45 = 612, Й44 = 622, Dn= 0 „з,

(6.25)

/ = ПЗ, Е 1= б 13,

Е2 ~ С23,

JD| = O13, Z>2 =

023.

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

указанной в § 1 структуры.

 

 

щ, v2 с

Закон Дарси,

связывающий

компопепты скорости

 

традиептами давления д\р, д2р в каждой из областей,

предста­

вим в виде [1, 6]

 

 

 

 

pvi = — andip — a i2d2p,

pv2 = - a 2\dip— а 22д2р,

(6.26)

где a ih— коэффициенты фильтрации, а « > 0 и a i2 = a 2i, р — плот­ ность жидкости. .

Обратный закон Дарси имеет вид

 

д хр =

— Pn^i — Pi2y2«

^ д2Р =

— РаЛ

Ра2у2*

(6.27)

причем

 

 

 

 

 

 

 

 

 

 

 

 

Р1 1 = дЦ р

 

 

Р“ = т Й Г ’

 

 

Д (а) = апа22 — а.{2> О, Д (Р) = рлр22 - pL -

д ^ Г

 

Уравпепие

неразрывности

течения

div(pv) = 0

с

учетом

(6.26)

приводит

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

второго

порядка

(1.4)

относительно функции

давлений р =

~ р (х 1, х2)

в

каждой

из областей структуры. Стало быть, для

функции р

справедливы представления (1.7).

 

 

 

187

Краевые условия па границах раздела сред таковы:

(РУ»)+ =(РУ„)-, Р+ = Г -

(6.28>

Здесь ил = V\cos ф +

sin ф — нормальная

составляющая ско­

рости на площадке с нормалью n = {cos ф, sin ф).

Очевидно, условия

(6.28) сводятся к

краевым условиям

(1.11). Потребовав, чтобы расход жидкости через любую кри­ вую, соединяющую конгруэнтные точки, был постоянным, при­ ходим к условиям кваэипериодичности (1.10).

Осредненпын закон Дарси имеет вид

<0,Р> = - < Pll><pPl> ~ <^12><РУ2>,

<diP> = —<P2^<Pyl^ — <P22><pl’2>, .

где <flifc> берутся из расчета, например, продольного сдвига ар­ мированной среды с учетом следующих формул перехода:

а л

= —Css,

« 1 2 =

-С 5 4 ,

<Х22 =

— С44,

Вп

- 0 5 5 ,

В)2 =

"054,

022 =

"

= Мз,

p v i = 0 iz,

p v z - а г з , pvn - а»з,

 

ei3“

diР,

 

 

( 6.80)

 

 

 

 

<Рп> __<Д55>, -<Pl2> ------- < а« > ,

< 0 2 2 > - - < в 4 4 > ,

< е,з > = < д х р \

< « 2 3 > - < 0 t f > .

Г л а в а 6

ТЕОРИЯ РЕГУЛЯРНО АРМИРОВАННОГО ВОЛОКНИСТОГО КОМПОЗИЦИОННОГО МАТЕРИАЛА

С ИЗОТРОПНЫМИ КОМПОНЕНТАМИ СТРУКТУРЫ

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

Пусть упругие и геометрические свойства трехмерной изо­ тропной кусочно-однородной среды, отнесенной к некоторой пря­ моугольной системе координат О х&^х3 с единичными направляю­ щими векторами эц Э , э3, неизмеппы в направлении 0х3. Распре­ деление упругих и геометрических свойств в плоскости Ох2 предполагаем двоякопериодическим.

Как следует из результатов [1, 2], важную роль при построе­ нии осредненных уравпепий, описывающих упругое поведение такой структуры, играет двоякопериодическое решение, когда компонента деформации е33 не зависит от всех координат, а ос­ тальные компоненты деформации— от координаты х3.

Такое напряженное состояние распадается на два линейно не­ зависимых [18]: обобщенную плоскую деформацию

щ =

щ (х 1, х2),

ejh = ejh(x u хц),

л = ол(х1, х2)

(/, к = 1 , 2 ),

е33 = const,

o33 = 2 |i(l+v)e33 + v(oii + O22),

и продольный сдвиг

е131= б23 — 0i3= O23= О

из (хц х2),

е}з = в]з(xi, х2),

из =

Oj3 = Oj3 (xi, Х2)

(/ = 1 , 2 ),

 

On = е\2 = е 22 =

е33 = О,

Он 012 — 022 — 033 = 0.

Здесь ejk, ой, и, (у, к = 1, 2, 3) — компопеиты тензоров деформа­

ции, папряжения и вектора перемещения, р. и v — модуль сдви­

га и коэффициент Пуассона.

 

 

 

 

Для решения этих

задач

представляется

естественным

ис­

пользовать развитую в

предыдущих главах

схему построения

комплексных

потенциалов,

соответствующих

двоякопериодиче­

скому полю напряжении.

 

в основном содержа­

При изложении § 1—2 мы следовали

нию работы

[101, § 3—7 представляют

собой расширенное

из­

ложение статей [7—9].

 

 

 

 

 

§ 1. Средние компоненты тензоров деформации

 

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

 

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

Будем через й| и й2 обозначать основные периоды структу­ ры, через 1\ и h обозначать |й|| и |й2|. Считаем, что система координат Ох\х2х3 выбрана так, что направление вектора й|

совпадает с положительным направлением оси

Ох\ и что угол

■а. между векторами Й| и й 2, отсчитываемый от

первого против

часовой стрелки, меньше я. Соответствующие векторам й| и й2

основные периоды в комплексной форме имеют вид on = l\,

ю2 =

= /2exp(ia).

 

Параллелограмм, построенный па векторах й| п й 2,

будем

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

Двоякопериодичность

компонент

деформации означает, что

вектор перемещения и в какой-либо

точке с радиусом-вектором

R отличается от вектора перемещения в точке с радиусом-век­

тором R + й, не более чем на вектор

малого жесткого переме­

щения

u(R + й^) =

u (R) -I- Aj + R X Ь^.

 

(1.1)

 

 

Напоминаем, что векторы R,

й ь й 2

 

лежат в плоскости

х\Ох2,

векторы bi, Ь2 перпендикулярны этой плоскости.

 

 

Вычислим при помощи (1.1)

значение

и(Н + й , + й 2)

двумя

способами: идя

через точку

с

радиусами-векторами

R + й ( п

R + й2. Условие

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

дает

 

 

 

Й, X Ь2 = й2 X Ьь

 

 

 

откуда следует, что bi = 0, Ь2 =

0.

 

 

 

 

 

 

Таким образом, из двоякопериодичпости тепзора папряжения

следует квазипериодичпость

вектора

перемещепия.

Запишем

последний в виде

 

 

 

 

 

 

 

 

 

u (R) - U* + *1

-^ Ctg-

+

 

^

па- Д2 + ^ 3Э33,

(1.2)