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

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

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

причем

* 1 1 - дw *и = лЩ '

Условие сохранения div q = О приводит к однородному диф­ ференциальному уравнению эллиптического типа относительно

функции W в каждой из областей <25, ЯА) dhm (к = 1, 2, . . . г;

т =

— 1 ,2 .........rfc) :

 

K n d\W + 2Kn dld.1W + K j i W = 0.

(1.4)

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

всякого з ^ 2D имеют место равенства

 

| qnds = ± l/A(/QlmQv = cons* (v = 1, 2),

(1.5)

где qn — нормальная составляющая вектора q на дуге, соединяю­ щей точки z и z + cD.v. Верхний знак соответствует Кц > 0, К 22 > > 0, пггжппй К\\ < 0, К 22 < 0 .

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

Вкаждой из областей Об, B h, dhm построить регулярные реше­ ния уравнения (1.4), удовлетворяющие дополнительным услови­ ям (1.5) и следующим краевым условиям на границах раздела компонентов среды:

w (t) = wk (/) + gh (<),

Чп(t)

= qn (/),

t<=Lh

(/c = 1 , 2 , . . . ,

r),

ivh (t) = whn(t) +

ghm(l),

qhn(t) =

<£"(0.

(1.6)

t e

h n {m = 1, 2,

. . . , /л).

 

Здесь gh{t), ghm(t)— заданные соответственно на Lk, ZAm произ­ вольные непрерывные по Гельдеру функции; под ними можно понимать, например, скачки температур на границах раздела сред.

Преобразуем описанную краевую задачу к более удобному для исследования виду. С этой целью выразим общее решение уравнения (1.2) через аналитические функции. Имеем

u? = Re<p(zo),

 

= Re фл (гл),

whm = Re флт(zAlH) ,

20 =

£i + ЦоФ,

zK = x 1+ |XAI 2,

zAm= X\+ pAm#2-

Здесь

 

 

 

- K ^ ± iV M K )

h> =

a +

*P’

Ц(2) =

Pit =

GCk +

*P/o

K .

Pfcm = a hm +

 

(1-7)

г е й ) , z eR ft, z e= dhrnt

171

причем верхний

знак берется при Кгч > 0, а пижний при

Кгг < 0 .

^

Для обозначения образов точек, кривых и ооластеи, получаю­ щихся прп аффинном отображении физической плоскости z на плоскости Zo, Zft, Zftm, будем приписывать к обозначениям одно­ именных прообразов в плоскости z соответственно пуль, штрих

и два штриха.

Соответственно этому запишем основные периоды в плос­

кости ZQ:

 

 

 

сою cot, ©го =

Re о>2 + ро

сог =

h + щ Н ho 4* iHo,

ho = h + aH,

Но =

( 1.8)

$Н.

Вычислим поток вектора q через произвольную кривую, соеди­

няющую точки А и В

в области 3).

 

 

Находим, учитывая (1.7) и (1.2),

 

Q = \ Чпds =

f (<7Хdx2 — q2dxy) = ptf221ш ф(z0) |д®. (1.9)

А В

А В

 

 

В силу (1.9), дополнительные условия (1.5) приобретают вид

Im (<p(z0 + Шлю)— <p(zo)) = Im

(v = 1, 2).

(1.10)

Из (1.10) следует квазипериодичность аналитической в области &>о функции <p(zo) по периодам сою и сего.

Краевые условия (1.6), учитывая (1.7) и (1.9), можно пред­

ставить в следующем равносильном виде:

 

 

ф (*о) =

елфл (h) + Е/*Фи (h) + gh (0.

t <= L h,

t0 <= L l,

th €= L'k.

(1.11)

фЛ {h) = ekm(Phm (hm) + Ейтф/ют{hm) +

ghm (0»

t e

l]lm>

th S

Ihvu

hm S

him.

Здесь

 

 

 

 

 

 

 

1 + h

 

*

1- Ч

 

 

=

2

'

=

~ 2

 

 

Н -Я,17П

*

1 - W

 

Р/|тп^22

2

 

Е/т------- 2

"w

Р А

 

 

 

{к =

1,

2,

... , г; т = 1, 2,

 

г*).

Постоянные интегрирования, которые, вообще говоря, долж­

ны фигурировать в

правых частях

(1.11),

включены в искомые

функции фЛ, фЛт.

Таким образом, приходим к следующей краевой задаче: опре­ делить квазпперподическую в Фо функцию ф(го) и регулярные

в областях Вп и dhm

соответственно фушщии фл(гл) и ф;,т (2Лт)

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

(1.11) и дополнительным условиям (1.10).

172

При этом подразумевается, что квазиперподлиность <p(zo) дости­ гается за счет выбора для пее специального представления.

Положим

Ф Ы

=

- t i I

[е W Р W “

8* (О PW1С («о -

2о) + А .

е 0 О,

 

 

гь

«

______

Л

 

фЛ Ы

=

2^]

j

(0 :— ®ктРкт (01

 

 

 

 

 

+ 2^jJ

~т~-----d£ft,

ZfeSBjI,

Z k m ^ d h m -

(1*12)

 

 

 

 

 

Lk h ~ Zb

 

 

 

 

 

 

 

 

 

<phm (z/tm) =

Лт f

t PkT^ }

dtkm

(A — 1,2,

, r; m —1, 2,

. . . , rA),

 

 

J

‘fern

 

 

hm

 

 

 

 

 

 

 

 

 

 

причем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<0=

0,5/ (1 — i(Jt0) +

0,5/ (1 +

ip0),

i e L ,

Z s i 0,

 

 

 

/л =

0,5/ (1 — ijXft) +

0,5/ (1 +

/рь),

Zs

Z/д

(J

ta,

tk e

L k U Z/tmi

Zfcin =

0,5/ (1 — tp/un) +

0,5/ (1 + ipfcm).

ZS

Ihmr

tkm ^

^ftm>

 

Здесь

^(zo)— дзета-функция

Вейерштрасса,

построенная

на

периодах

©ю и ©го,

p (t)=

{ph[p)t

t ^ L k\t

e ( t ) = i e ht t ^ L h)

и

e*(Z) = {e£, Ze

J;

Л — подлежащая

определению

константа.

Обход при интегрировании

вдоль

hm> h

ведется

по

часовой

стрелке, а

при

интегрировании

по

Lk,

L h — против

часовой

стрелки. Очевидно, функция q>(zo) в

(1.12)

квазипериодична в 2£)Q.

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

(1.5)

дают возможность ввести средние по ячейке потоки <qi> и

<?2>

(рис. 5.1.1)

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<</i>=^

J

 

 

 

 

|

 

 

 

 

 

 

 

 

 

 

 

\

 

 

 

1

г+"1

 

 

 

 

(1.13)

 

 

<q2> =

 

j*

 

qnds,

Н =

Im ©.,,

а =

arg ©,.

 

 

 

 

 

1 z+ffl

 

 

 

 

 

 

 

 

 

 

 

Соотношения (1.13) с учетом (1.9) и (1.10) связывают средние потоки с величинами Im Qi, Im £?г:

<?i>------ № г 1ш F ) ' (?!> — Pr = IlIl| <114)

Всоответствии с этим в дальнейшем будем считать, что поле

вструктуре вызывается действующими в пей средними потоками

173

Подставляя в дополнительные условия (1.10) приращения: <p(z0) из (1.12), приходим к уравнениям

1т(Л<й«о — я&ю)=

 

(v = 1, 2),

(1.15)

W

,

о

'

 

 

 

 

a - E T j

[e (0 / » W -e *(t)p (0 ]*o .

 

 

 

Ч

 

 

 

 

 

6vo =

С (zo +

“ vo) —

?> (zo) =

2С ( “^ ) -

 

Решение этих уравнений с учетом

(1.14)

дает

 

где

 

А = А Р + АЧ,

 

 

(1.16)

 

 

 

 

 

 

 

= ^

а “

7 ? Кеа’

 

=

<■>,„ Im (о2П,

 

А

=

д(8) {<9i> +

<<72>1-

 

Таким образом, представления (1.12) определяют кваанперподическую в £>о функцию <p(zo), удовлетворяющую дополнитель­ ным условиям (1.10).

§ 2. Интегральные уравнения регулярного поля

Подставляя предельные значения функций (1.12) в краевые условия (1.11), приходим к системе интегральных уравнений от­ носительно функций рк, phm:

Рн(т ) - Мкiph(t), phm(t), т) = F h(т), т е Lh (к = 1, 2, . .. , г),

 

 

 

 

 

 

(2.1)

Рхш(т)— Мкт iphm(t), ph{t), х\ — F hm(т), x ^ L hm (m = 1, 2,

... ,

?*h),

+ 2iir h

J

Iе (f) P (*) -

e* (0 P (01 C (<0 - To) dt0 +

 

 

+ Т0Лр +

2

л

J

l

Phm (t) — 4m Phm (0) -

 

 

 

 

l/{m

1

 

 

 

 

_

J _

V

J" l^hmPhm (t) fytmPhm(f) 1j

^ '»

 

 

 

2nt ^

 

 

 

 

m-i Jfiw

 

 

174

jVkm -

2 J

(<) — ehsPh* (<)1

 

lhs

 

 

Штрпх пад суммой в выражении для Мкт означает, что слагаемое с номером s = т необходимо опустить.

Легко видеть, что ядра в (2.1) могут обладать не более чем слабыми особенностями, т. е. система (2.1) фредгольмова. Ниже будет доказана ее разрешимость.

 

 

§ 3.

Теоремы единственности

 

 

Будем предполагать,

что

решение

краевой

задачи

(1.11),

ДНО)

существует.

 

 

 

 

 

 

 

 

Т е о р е м а

3.1.

Между любыми двумя решениями

краевой

задачи

(1.11)

ф(1) (z„),

 

(zft),

(zhm) и

ср(2) (г0), <p(ft2

Ы , Фа« fom),

каждое

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

(1.10), имеют место соотношения

 

 

 

 

 

Ф(0) =

ф(2) (z0) -

Ф(1)(20)= С , Chm=

Re С+ j- i — b i С,

 

 

 

 

 

 

 

 

ЛЬЛЛт71

 

•фк0) =

ф!2) (Z0) -

(z0) = Сл

(ft — 1,2.........?•;

HI =

1,2.......... г*),

 

 

 

 

 

 

 

 

 

 

(3.1)

 

Цкт=

ЦШ (Zhm) ~

Cf>im(Zftm) = Cfcmi

Ch. =

Re C + j- Im C.

Для доказательства заметим прежде всего, что для любого регулярного в области В решения уравнения (1.4) имеет место

«энергетическое» равенство

 

 

j j |d2W - pdJY I2 dx1 dx.2 =

f Wqnds.

(3.2)

в

22 LB

 

Здесь L B — граница области В , qn — нормальная составляющая потока q, введенного в (1.2); обход при интегрировании — против стрелки часов.

Формула (3.2) выводится при помощи обычных преобразова­ ний, вполне аналогичных тем, которые приводят к интегральной

175

формуле Дирихле в теории потенциала

[4]. При К п = Я 22, К Х2=

= 0 она и совпадает с последней.

 

 

 

 

 

Применяя

(3.2)

к нашей многокомпонентной структуре, полу­

чим с учетом

(1.11)

и (1.7)

 

 

 

 

 

 

/ =

К 22 f J

|d2W -

n dJV |2 dx, dx„ +

 

 

 

 

 

b v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ jS

* и И

|

-

tydjWj |2 dx, dx., +

 

 

 

7*

rj

 

 

 

 

 

 

 

 

 

 

 

+

2^ ^2^ K {2 ( j

I d2Wjm

 

|- dx, dx., =

 

 

= JWgn ds +

2

f

 

(0 ds +

2 2

f

cjU’jm (t) ds.

(3.3)

 

Г

 

 

J - * L j

 

 

 

 

4 *

 

 

В (3.3) контуры lj„,

Lj

и Г

(граница основной ячейки По)

обхо­

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

— г -Ь 1-связная область с грани­

цей Г U L.

 

 

 

 

 

 

 

 

 

 

 

 

В

силу

квазипериодпчности <p(zo)

находим

с учетом

(1.7),

(1.5)

и (1.9)

 

f wqn ds = Im(йД,) / д,

 

(3.4)

 

 

 

 

 

где

 

 

 

г

 

 

 

 

 

 

 

 

 

 

Qv =

<p(z0 +

C0v0) - <p(z0)

(v = 1,

2).

 

 

 

 

 

Подставляя (3.4)

в

равенство

(3.3)

и применяя затем

(3.3)

к разности двух решений краевой задачи

(1.11), удовлетворяю­

щих условию

(1.10), приходим к требуемому утверждению.

 

Функции Ф , фл , фля»

можно интерпретировать как решения

однородной краевой задачи

(1.11), соответствующей

 

 

 

 

&(*)-0; g hm(t) = 0 ,

<?,>=0

<q2>=0.

 

Т е о р е м а

3.2.

 

Пусть

%i(t)

и сД£) — граничные значения

функций %f(zo) и Oj(zj), регулярных соответственно в конечной области на плоскости zo, ограниченной кривой L j, и в дополне­

нии

до полной плоскости zj.

Если

в окрестности

бесконечно

удаленной точки Oj(z}) =

0(lz*l-1), то краевая задача

 

 

Ъ (О = W i (0 +

e*o7(ij,

t<=Lj

(/ =

1,2,

. . . ,

/•),

имеет лини» тривиальное решение.

 

 

 

 

Для доказательства

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

(3.2)

к

двухкомпо­

нентной области, ограниченной окружностью CR достаточно боль­ шого радиуса.

176

Имеем

Ро^22 |* j I X (Zo) I” d x i dX2 + Pj^-22 \ j

I Gj (Zj) |2 C?Zj C?Zo =

 

 

 

=

f <7n Re a,- (/j) <fc,

(3.5)

 

 

 

cR

 

S)n — двухсвязная область с границей L3UCn.

 

Неограниченно

увеличивая R,

приходим к требуемому

ре­

зультату..

 

 

 

 

§ 4.

Разрешимость системы

(2.1)

 

Докажем, что при сделанных выше предположениях относи­ тельно граничных лнппй н функций gh(t), gAm(f) система инте­ гральных уравнепий (2.1) всегда разрешима. Для этого рассмот­ рим соответствующую ей одпородпую систему.

Очевидно, для того чтобы /',,(т)=0, Fhm(x) = 0, необходимы и достаточны условия

Ы = 0, ghm(т) = 0 =

1 , 2, ...,

г; т = 1 , 2, ..., гА)

<2,> = 0.

<<7г> = 0.

(4,1)

Таким образом, однородная система соответствует однородной краевой задаче (1.11) при однородных дополнительных усло­

виях

(1.10).

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

решения

одпороднон системы

(2.1) через

Ph(t)T.

P U

I ).

 

 

 

 

 

 

 

 

 

 

 

 

На основании теоремы единственности 3.1 запишем

 

 

ф(0) Ы = ^ - j Н О Ро (0 -

е* (t)

 

£ (t0 -

ZQ) dt0 + A0z„ =

C,

 

 

L °

 

 

 

 

 

 

 

 

 

 

 

 

2,

£

cpffi Ы

-

{

 

Я ш =

C„„

(4.2)

 

 

 

 

 

 

n

l hm

"Ant

 

 

 

 

 

 

 

 

 

 

lkm

 

 

 

 

 

 

фк0)(Zk= 2лiJ

 

 

 

[ leAm7*Am(0—

 

 

 

 

Llt

 

 

 

lhm

 

 

 

 

 

 

 

 

 

&umPkm (t) ) - - -

-

= C;t,

Zfc e

2?;t,

z,tm e

« W

 

 

 

 

 

SA

 

 

 

 

 

 

 

Вычисляя приращения

функции

q)°(zo)

 

в

первой

формуле

(4.2)

при переходе от точки z

к

конгруэнтной точке

г + a>v, по­

лучаем

 

 

 

 

 

 

 

 

 

 

 

 

 

2^ f

U (9 Ро (9 — е* (<) /й*)} С(^о — zo) dlo =

 

 

(4-3)

 

ь°

4 о

= 0,

а0 =

0,

Zo е 0 ° .

 

 

 

 

 

 

 

 

 

 

 

 

12 э. И. Грнголгок, Л. А. ФилыцтнискшЧ

 

 

 

 

 

 

 

 

177

 

 

 

 

 

 

 

 

 

Из (4.3) следует, что функция t(t)p 0(t)~ e*(t)p 0(t) является граничным значением некоторых функций, регулярных в коноч­ ных областях плоскости zo, ограниченных контурами

= 1,2, . . г). Поэтому интеграл в (4.3) исчезает, и полу­ чаем на основанпп (3.1)

£ = 0, Ск = 0, Скт= О (& = 1, 2, ... , г; m = 1, 2, ..., гА) . (4.4)

Введем в рассмотрение функции

i%n(I)= ehPh(t) —£*kPk(0.

ton V) -

Ph (t),

t<= Lh

= 1,2, . . . , r),

(4.5)

*XAm (0 =

£hmPhm (0 — ehmPhrn(£)»

 

^Am (0 =

Pkm(t)i

t S Z/ITJI

(//1= 1, 2, . . ., ?'/i).

 

Из (4.2) с учетом (4.4) н (4.5) заключаем, что %k(t) является граничным значением функций %K(ZQ), регулярных в коночных

областях, ограниченных кривыми L k; 0 /,(t) — граничное значение

функций Ofc(zft), регулярных вне Ж>к и исчезающих па бесконеч­ ности; Xhm(t) — граничное значение функций xkm(zhm), регулярных

в конечных областях, ограниченных кривыми 1кт,' ы, пакопец, <yhm(t) является граничным значением функций afcm(z„m), регу­ лярных вне областей dkm п исчезающих на бесконечности.

Исключая из (4.5) функции Ph(t), pkm(i). ‘приходим к спсте-

ме независимых краевых задач

 

 

Ха(0 = ehah(t) + e,!oft(/), /.<= L h (к =

1 , 2 , . . . , ?•),

 

_____ :

 

(4.6)

Хат(0 = ehmakm (t) + etm(jhm (t), t <= lkm

(m = 1, 2,

. . . , rh).

Функции oft(zA) n ahm(zhm) в окрестности бесконечно удален­

ной точки затухают не медленнее, чем lzfth l и lzAm|-1 соответ­ ственно.

В силу теоремы 3.2 имеем

М * ) =

Х * ( 0 = 0 ,

Oum(t) =

x hm( t ) = 0

 

(к =

1, 2,

... , г; т = 1,

2.........

г„).

(4*7)

Отсюда, на основании (4.5), заключаем

 

 

Ph (0 = 0, p lm(l) = о

(Лс =

1, 2, . . . ,

г;

/л = 1 ,2 .........

rh). (4‘.8)

Таким образом, система уравнений (2.1) всегда имеет реше­ ние и притом единственное.

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

478

Устремляя периоды к бесконечности, получим решение крае­ вой задачи (4.5) для мыогосвязной области с конечным числом комнопсптон.

 

§ 5. Модель регулярного поля (регулярной структуры)

ды

Положим ?а(0==0, <7«п(0в 0- В каждом из компонентов сре­

имеет место скалярное поле, определяемое функциями ш,

wky

whm. Поток в каждой точке поля задается формулой (1.2),

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

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

с этим мы определили средипе потоки <д\>

и <</ > формулами

(1.13).

Нод макромоделыо регулярной структуры будем понимать здесь однородную анизотропную (в определенном смысле) сре­ ду. уравнения состояния которой совпадают с законом связи между средними потоками <д(>, <д2> и средними градиентами <д\u>>, <dow> в структуре.

Средине градиенты определяются соотношениями

 

co sily )

= R cQ ,,

hidiiv) + #<<?2ц;> =

Re Q2.

(5.1)

Выражая Re Qv

(v = 1,

2) через средние

потоки,

находим,

учитывая (1.12) и

(1.16),

 

 

 

<?!> — тг-Нев,

B . Q . - - &

* < * > + “

 

 

 

 

 

-

М

 

" - “ ■■<^>— -^-Re(a(o20).

 

 

А (К)

 

А (1с)

 

 

 

 

 

 

Вставляя эти значения Re fiv

(v =

l,

2) в

 

формулы (5.1), по­

лучаем

 

К» . .

 

 

 

 

 

 

 

 

<5дш> =

^

<

92> _ ^ R e a ,

 

 

<<7i> _

 

 

 

 

 

 

 

 

 

 

 

 

(5.3)

< » . » > -

- т а г < «i>

+ т ш < & > -

т

1 R e №

) •

 

Пусть

р(1) (t) — решение

 

системы

(2.1)

при

<gi> =

p_/f22,

<д2> = 0,

а p l2\(t) — решение

этой

системы

при <gi> = 0, <g2> ==

= ВКм- Тогда

<?i>

 

< 0

 

 

 

 

 

 

 

 

 

 

 

 

 

a= =K 7 a i + f e r a-

 

 

 

(5-4>

где функционал ah (к = 1, 2) определен на стандартном решении системы p w (i) формулой

0 1 . - 2 г !

(5.5)

L\

12*

17&

Подставляя выражение (5.4) в соотношения (5.3) и учитыиая (1.3), приходим к уравнениям состояния макромодели структуры

<dl W> =

.< Я ц Х д | > +

<Х|2><52>,

 

 

(5.6)

<д2гиУ =

<X2i><gi> +

<«22

 

 

 

 

 

 

 

<xu> = xu (l -

Ц - Re ttj),

<Х 12> =

хи -

^

Х и Re о2,

<««> =

«г1 [ l

-

 

Ио ( ^ - ) ].

 

 

 

Коэффициенты

<ха»>

назовем

макроскопическими

(эффек­

тивными, средними) параметрами структуры.

 

 

 

 

Т е о р е м а 5.1.

Макроскопические параметры

<х<щ>

образу­

ют симметричную положительно определенную матрицу х.

Для доказательства подставим

в

правую часть

(3.4)

вместо

Refiy, lm&v их выражения из (1.14) и

(5.2). Тогда (3.3) при­

мет вид (суммирование по повторяющимся индексам)

 

7 = < g e><5eU>>

F

 

( а = 1 , 2).

 

 

(5.7)

Вводя стандартные решения wy\ и>у), w}tm(v =

1, 2)

по формулам

w = <<7V) w(v\

wh = <<7У> w{u\

whm =

<r/v>

 

(5.8)

и подставляя их в левую часть

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

 

 

 

 

I = I vll<qM qll> = F < qa><daw>

 

(л>, ц =

1,

2),

(5.9)

причем I \2 = h\ и смысл

ясен из кЬнтекста.

 

 

 

 

В силу произвольности <gv> приходим к

закону (5.6). Мат­

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

ложительной определенности функционала

I

в . (3.3)

d e tx > 0 .

Функционалы ai и а2, фигурирующие

в

<хар>,

учитываю^

мшфоструктуру ячейки. Если структура отсутствует, т. е. сре­ да однородна, то а х = а2= 0 и (5.6) совпадает с уравнениями состояния (1.3).

Отметим, что построенная теория справедлива и в том слу­ чае, когда материал какого-либо компонента структуры изотро­

пен. В этом

случае соответствующее характеристическое

число

ц = £. Результаты этого параграфа суммирует

следующая

теоре­

ма, определяющая решение проблемы осреднения.

 

 

Т е о р е м а

5.2. Если в регулярной

структуре

имеет

место

гармопическое поле W с уравнениями состояния (1.2) или (1.3),

вызванное

средними в

смысле (1.5),

(1.13)

потоками,

то су­

ществует

макромодель

структуры, управляемая

законом

(5.6).

180