книги / Периодические кусочно-однородные упругие структуры
..pdf§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) |