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

книги / Метод эффективного поля в механике композитных материалов

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

561

J* ^а^сфЛ.ц {x)(Kl = ^ea{x)AapX tl{x )d r .

(ПЗ.1.9)

n

г

 

П3.2. Производные потенциала двойного слоя статической теории упругости

Пусть вначале Q - замкнутая поверхность Ляпунова, огра­

ничивающая односвязную конечную область V в R3. Рассмот­ рим потенциал

e j x ) = |K ^ , ( x - x ' ) C ^ ( x % ( x ’ )dCr, (ПЗ.2.1)

Q

К а ^ ( х ) = > (ПЗ.2.2)

где G( x) - функция Грина статической теории упругости для однородной среды с тензором модулей С°. Из структуры под­ ынтегрального выражения следует равенство £(;t)=defw(;t), в котором и(х) - потенциал двойного слоя статической теории упругости, определенный соотношениями (1.1.26), (1.1.27). Для любой непрерывной плотности b (х) поля и(Х) и е(х) непрерывны всюду за исключением поверхности Q , причем разрыв и(х) на О определяется соотношением (1.3.28).

Рассмотрим интеграл (ПЗ.2.1) при Z>=const. Применяя фор­ мулу Гаусса и учитывая, что функция G(x) удовлетворяет уравнению (1.1.6)

V a C ^ V x G jx ) = -S fiA x ) ,

(П3.2.3)

получим равенства

 

£afi(X) = J

aflXM(X ~ X')C°X/Jvpdx'bp =

(ПЗ.2.4)

V

 

 

= JУ[ад(х- x'Yix^dx'b^ щ а{х)ьлП(х),

562

где П(лс) - дельта-функция, сосредоточенная на поверхности Q , п (х) - внешняя нормаль к Q .

Пусть теперь Ь(х) в (ПЗ.2.1) есть линейная функция вида

ba(x) = Dafi{ x - x ') fi,

(ПЗ.2.5)

где D - постоянный тензор. Рассмотрим предельные значения

потенциала (ПЗ.2.1) при х—»xoeQ.

Обозначим через Пр часть поверхности Q, расположен­

ную внутри сферы малого радиуса р с центром в точке хо, а

через Qpдополнение Qp до Q. Введя локальную декартову

систему координат ух,у2,уг с началом в точке хо и осью у3,

направленной вдоль нормали п(хо)=по, представим интеграл (ПЗ.2.1) в форме

y')C°n{y')b{y')dCl' =

(ПЗ.2.6)

Q

=| К {y-y')C°n{y')D{y-y')dQ'+ J K{y-y')C°n{y')D(y-y')dQ'.

п,

Интеграл по Qp есть, очевидно, непрерывная функция в точке у= 0 (х=хо).

Найдем предел инте1рала по О р при у —>0. Поскольку в

дальнейшем р-У 0, можно считать, что Qp - плоская круговая

область и п(у')=по. Введя замену переменных у,=\у\4,, будем иметь

lim JК(^ - y')C°noD(y - y ’)dQ! =

(ПЗ.2.7)

^ Ч

563

где интеграл справа вычисляется по всему пространству R3 и учтено, что К(дг) - четная однородная функция степени (-3).

Представим функцию «„£(«„•£') под интегралом в (ПЗ.2.7) в форме

» . А « . - д = ч Л п . - ? ) , я (')= {о 1 ',<о>

(ПЗ'2'8)

и, интецжруя по частям, найдем

 

| к ( £ - ?)C-n.D(Z- ? ) # « . •? =

(ПЗ.2.9)

= -Н (п . 4 )D + | K (f- ? )Н (п . ? ) d ? C D .

 

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

JК (£- £ )Н {п ? ) d ? =(2я) 3J}С {к )Н '{к )е -‘кЧ к =

= 1 [Г (0) + sign(«0•£ ) К » ] .

(ПЗ.2.10)

Здесь учтено, что преобразование Фурье функции Н(по £) имеет вид

H*(k) = H*{kl,k2,ki) = {2 я)2 ^ к х)^ к 2 )[л £ к г)+ ik? ].

(ПЗ.2.11) Обозначим через J+(no) значение предела (ПЗ.2.7) при

стремлении у к нулю со стороны нормали по, а через J_(no) - с противоположной стороны. Очевидно, что эти пределы не

зависят от размера области С1р и в силу ПЗ.2.9, ПЗ.2.10 имеют

вид

y+(/iJ4[K*(o)+K*(/ie)]ceD-z), J (HJ 4 [ K *(O)-K > ,) ] C °D .

(ПЗ.2.12)

При замене по на (-«„) интеграл (ПЗ.2.7) меняет знак, тог­ да как предельные значения его по-прежнему определяются

соотношениями (ПЗ.2.12). Отсюда следует равенство ./+(ио)=

564

=-J_(no), которое приводит к соотношению К*(0)=(С°) 1 (сравни с 1.2.21).

Запишем теперь выражения для предельных значений по­ тенциала (ПЗ.2.1) с плотностью (ПЗ.2.5) в точках поверхности Q . Устремляя р к нулю в (ПЗ.2.6) и учитывая (ПЗ.2.12), полу­

чим, что при х— предельные значения рассматриваемого потенциала со стороны нормали е+и с противоположной сто­ роны Е имеют вид

п

(П3.2.13) где символ / означает интеграл в смысле главного значения по Коши, который существует в силу четности К (х ),

Л сфл)л (И) = ^

— Е арХм ,

(ПЗ.2.14)

Ех - единичный четырехвалентный тензор.

Пусть теперь в (ПЗ.2.1) функция Ь ( Х ) бесконечно диффе­

ренцируемая на Q (Z>(x)eCco(Q)). Представим потенциал f(jc ) в форме

£ -Jx)=[K am{x-x')r^lx')[b^x')-b^x)-dJ>p{x){x'-x)^ £ l’+

Q

+|

- x']C ^ X x')d n 'b f(x) +

(ПЗ.2.15)

Q

 

 

+J "x')C°Xllvtp v{ ^ ) ^ bP{x){x' ~ x)dQ!. a

Здесь первый интеграл является непрерывной функцией во всем пространстве и сходится абсолютно при всех х. При

х^1 его можно понимать в смысле главного значения и представить в виде суммы двух интегралов с плотностями b(x’)-b(x) и cb(x)(x'-x) соответственно, каждый из кото­

565

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

дельные значения потенциала £(лс) на Q имеют вид (х ^ 2 )

= | К afiXM( x - x ’)ClMv/iv{x') х

(ПЗ.2.16)

О

х [^ (х ') -6 р( х )]б /О Т 1 Л ^ (и )^ /1(х),

где тензор Л (И) определен соотношением (ПЗ.2.14). Интеграл здесь существует и при менее жестких шраничениях на b (X) [158].

Рассмотрим теперь потенциал <г(л;) вида

= jS ^ ix - x '^ ix '^ ix ^ d Q ',

(ПЗ.2.17)

О

S„02„(.л

где Q - по-прежнему замкнутая поверхность в R3.

Очевидно, что потенциалы (ПЗ.2.1) и (ПЗ.2.17) связаны

между собой соотношением

 

° J x ) =

(ПЗ.2.18)

Поэтому из (ПЗ.2.16) следует, что его предельные значения на О определяются соотношением

*) = l Safix^x-x,)nx{x ,)[btl{x ') -b tt{x)\iQ: +

Q

* т £ й>)<?а М>

sr(n)=ctr(«)c--c.

 

(ПЗ.2.19)

566

Учитывая выражение для К*(п) (1.1.35), можно убедиться,

что вектор паиар(х) - непрерывен при переходе через поверх­

ность Q (naS*apX)i(n)=0). При х e Q выражение для вектора

<тар(Х) имеет вид

 

па(х)аар{х) =

 

Q

 

Т<ф(х’х') = -n xi^ S x a frix -x ')^ ')-

(ПЗ.2.20)

Если поверхность О не замкнута, томожно считать, что интегрирование в (ПЗ.2.16), (ПЗ.2.19) проводится по любой гладкой замкнутой поверхности, включающей Q, b (Х)=0 вне

П. Пусть Q - дополнение П до замкнутой поверхности. Учи­ тывая представление (1.2.4) для S(x)

SafiXft{x) = rot^ rot^ Z ^ J x ),

(ПЗ.2.21)

где Z (x )

- тензор Грина для внутренних напряжений [80], и

теорему Стокса (П3.3.4), нетрудно получить равенство

 

=

(ПЗ.2.22)

о

г

 

где dTM- векторный элемент длины на контуре Г - границе

Q , ориентация которого согласована по обычному правилу с

ориентацией Q . Таким образом, в случае незамкнутой повер­ хности соотношение (ПЗ.2.19) принимает вид

0*efi(x) = fS a fix J x -X 'M x 'lfo M -h J x fy n ' + Q

+f rotvxZapv,{x - х')сЛГ[Ьм(х) + \S'af}Xfi (n)<9xbM{x),

г

(П3.2.23) где положительный обход контура Г согласован с ориентаци­ ей поверхности О . Отсюда и из ПЗ.2.16, ПЗ.2.17 следует, что

567

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

* * (* ) = lKapX^X- X')Cl»VM X'\ h (X')-bp(X)]d& +

Q

+(с°)^Л/1irot^Jx - x')dT£Xx)+1 («)% (*) •

г

(ПЗ.2.24) Выражение для непрерывного на Q векторного поля

п(х)<т(х), где сг(х) потенциал (ПЗ.2.17), в точках незамкну­ той поверхности Q имеет следующий вид:

nfi{x)<rfia(x) = - fT afi{x,x’)[bfi{x')-b fi(x)'^in' + Q

+ f «/>(*)го\уХ2а^(х-х')с1Г'хЬ^х).

(П3.2.25)

ПЗ.З. Потенциал с плотностью, являющейся тензором поверхности Q

Рассмотрим потенциал типа (ПЗ.2.1)

еар(х) = \KafiXft{ x - x ,)qXfl{x')da, ,

(П3.3.1)

о

 

с плотностью q ( X ) 9 являющейся тензором поверхности Q:

ПаЯа(>(Х) = в, ®afiXM(X)<IxM(X) = <lafi(X)> * s Q . (П3.3.2)

Здесь проектор ®(дг) имеет вид (3.2.4)

 

0(дс) = ®(и) = Е' - 2 Е5(п) + Е6(п),

(ПЗ.З.З)

ОСч

П(х) - нормаль к

(W) - элементы основного базиса П1.1.1,

Используя теорему Гаусса для поверхности (ПЗ.1.9), най­

дем

568

= - JV^V('aG^(x - x'УдXflVp{x')qvp(x')dQ.^ = (П3.3.4)

Q

=-J ^ [v ('aG^(x-x')^(x')]^Q'+J ^[aGpJ<x -x ,)ffxqJJ<x,)d a= Q Q

=J^ (с£*р)^х~X )^хЯхм{х }d£l —^V(aG^)/i(x—^x )<Jx/i(x ')ex(x )dr', а г

где ел(х) - нормаль к Г, определенная в п.3.3.

Здесь учтено, что 0 ^ ,( х ) У Л^ _.„(*) = © ^ ( х ) ^ Дх),

да - оператор дифференцирования по поверхности (ПЗ.1.8).

Первый инте!рал в этом соотношении представляет собой потенциал (1.3.1), предельные значения которого на поверх­ ности определяются соотношением (1.3.25). Таким образом, предельные значения потенциала (П3.3.1) на Q имеют вид

S±J Kx) =

-

(П3.3.5)

 

О

 

Г

Л аДц(«) =И( а ^ ( « ) , « = П(х) , Х€П.

Инте1рал в смысле главного значения в этом соотношении можно переписать в виде

/ V ;A J * - x 't e l j x ') * * ' =

(П3.3.6)

a

=l^ {aG ^ {x - x')dx[qXtl{x')-q lfl{x,x')]dQ.\

а

где - тензорное поле на Q , удовлетворяющее уравнению

569

ЗД Д *,*') = 0, х , х ’ е П

(П3.3.7)

исовпадающее с q(X) при х' = х.

Спомощью теоремы Гаусса для поверхности П3.19. интег­ рал (П3.3.6) приводится к виду

=

(П3.3.8)

Q

= - 1 ^ V 'aG ^ (x - x ^ q j x 1) - ^ ( х ,х ') р П ' +

Q

+ | V ( ' ( * - *')[<7л„( Х ' ) - Я х М( Х ’ * ') \ Ю

Г

Отсюда следует регулярное представление предельных зна­ чений потенциала (П3.3.1) в точках поверхности О в форме

< P{x) = l K afiXfl{x -x')[q X)i{x’) - q Xtl{>г,х')]Л2'- (П3.3.9)

Q

~j У[авлм(х '~ х')я1м{х,х')ех{х')dT' ± \ K ^{n)dxqXtt{x).

г

Заметим, что из определения тензора А° в (П3.3.5) следу­ ет, что "касательная" составляющая тензора е ^ х ) на поверх­

ности Q офХ^{х)еXfl(X)) непрерывна на П и определяется следующим соотношением

® «**(*)**,(*) = l U ^ { x tx,\q^{x,)-q lJ <xyx,)]fiV -

О

- &архм(хЦ ^х^^(х - х')я°РХх,х’)еР(х’)<^, (П3.3.10)

г

У ф г(Х>Х') = Qafivpi^K^ix-X ')®тЗХм(х').

ПРИЛОЖЕНИЕ 4

МАТРИЦЫ И ВЕКТОРЫ ПЕРЕХОДА ЧЕРЕЗ ГРАНИЦЫ СЛОЕВ В ЗАДАЧАХ

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

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

Матрица перехода Г' через границу /-го и (/+1)-го слоев в

точке r=at (/=1,2....N ) имеет размерность (12x12) и пред­

ставляется в форме (® - прямое (декартово) произведение матриц)

Г' = Г / 0 Г '0 Г ' ,

(П4.1.1)

где матрицы Г,' (8x8) и (2x2) зависят от модулей упругос­

ти /-го и (/+1)-го слоев и имеют вид

 

1

0

0

0

0

0

0

0

 

 

0

1

0

0

0

0

0

0

 

 

Г31

Г32

Г33

0

Г35

Гза

0

0

 

Г/ =

Г41

Г42

0

Г*

Г45

0

0

0

, (П4.1.2)

0

0

0

0

1

0

0

 

0

0

0

0

0

1

0

0

 

 

Г71

Г72

0

0

Г75

Г7б

Г77

0

 

 

Г81

 

0

0

Г85

Г86

0

^88

 

_ 1 г г

= г = Mi

г = —г

31>

г 31= -

, г 32 2 1 31>1 33

1 44

Г35 = 136

3

Мм

 

Мм

 

 

 

Соседние файлы в папке книги