О ПРИМЕНЕНИИ МЕТОДА ГИПЕРСИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ К РЕШЕНИЮ ЗАДАЧ ДЛЯ УПРУГОЙ ПЛОСКОСТИ С КОЛЛИНЕАРНОЙ СИСТЕМОЙ ТРЕЩИН

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

УДК 539.3 DOI 10.18522/1026-2237-2020-2-72-83

О ПРИМЕНЕНИИ МЕТОДА ГИПЕРСИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ К РЕШЕНИЮ ЗАДАЧ ДЛЯ УПРУГОЙ ПЛОСКОСТИ С КОЛЛИНЕАРНОЙ СИСТЕМОЙ ТРЕЩИН

© 2020 г. С.М. Мхитарян1&2

ON THE APPLICATION OF THE METHOD OF HYPERSINGULAR INTEGRAL EQUATIONS TO SOLVING PROBLEMS FOR AN ELASTIC PLANE WITH A COLLINEAR SYSTEM OF CRACKS

S.M. Mkhitaryan1&2

institute of Mechanics, National Academy of Sciences of the Republic of Armenia, Yerevan, Armenia, 2National University of Architecture and Construction of Armenia, Yerevan, Armenia

Мхитарян Сурен Манукович - доктор физико-математических наук, член-корреспондент НАН Армении, профессор, заведующий отделом механики упругих и вязкоупругих тел, Институт механики Национальной академии наук Республики Армения, пр. Маршала Баг-рамяна, 24б, г. Ереван, 0019, Республика Армения; заведующий кафедрой высшей математики, Национальный университет архитектуры и строительства Армении, ул. Теряна, 105, г. Ереван, 0009, Республика Армения, e-mail: smkhitaryan39@rambler. ru

Suren M. Mkhitaryan - Doctor of Physics and Mathematics, Corresponding Member, National Academy of Sciences of Armenia, Professor, Head of the Department of Mechanics of Elastic and Viscoelastic Bodies, Institute of Mechanics, National Academy of Sciences of the Republic of Armenia, Marshal Bagramyan Ave., 24b, Yerevan, 0019, Republic of Armenia; Head of the Department of Higher Mathematics, National University of Architecture and Construction of Armenia, Teryan St., 105, Yerevan, 0009, Republic of Armenia, e-mail: smkhitaryan39@rambler. ru

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

In the present paper, using the method of hypersingular integral equations, based on the formulas of the inversion of the corresponding singular integral equations, the exact quadrature solution of the classical problems of the mechanics of an elastic plane with a collinear system of cracks is constructed. The elastic plane is in a state of antiplane deformation or plane deformation; in case of antiplane deformation, crack edges are symmetrically loaded by tangential forces, while in case of plane deformation, they are again loaded symmetrically but by normal forces. Mixed boundary-value problems for an elastic half-plane equivalent to these problems are formulated. Under plane deformation, the mixed boundary-value problem for an elastic half-plane is discussed as well when the plane boundary is reinforced by two similar and symmetrically located semiinfinite stringers between which a system of an arbitrarily final number of stringers is situated. It is considered that the stringers are absolutely rigid for expansion and compression and absolutely flexible for bending. A particular case of two

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

similar symmetrically located cracks is considered more in detail. In this case, the exact solution to the problem is also constructed by the method of Chebyshev orthogonal polynomials.

Введение

Вслед за теорией сингулярных интегральных уравнений (СИУ) и связанных с ней методов краевых задач аналитических функций в последние десятилетия интенсивно развиваются исследования по теории гиперсингулярных интегральных уравнений (ГСИУ) и методам вычисления гиперсингулярных интегралов в смысле Адамара. ГСИУ находят многочисленные приложения в задачах механики разрушения, геомеханики, гидромеханики, акустики, в теории теплопроводности тел, теории потенциала, в разнообразных областях математической физики и прикладной математики [1—5]. ГСИУ в идейном и методологическом аспектах тесно связаны с СИУ и теорией краевых задач аналитических функций, которые широко используются в исследованиях по контактным и смешанным задачам теории упругости и механики разрушения [6-11].

Из часто встречающихся в приложениях ГСИУ наиболее глубоко как в теоретическом, так и в вычислительном плане разработана теория классического ГСИУ с ядром — X) 2 на одном интервале, допускающего точное решение [12, 13]. Это же ГСИУ на двух симметрических интервалах рассмотрено в [14], где решение представлено сложным двухкратным интегралом с внутренним сингулярным интегралом от ядра Коши. Класс точно решаемых ГСИУ расширен в [15], где рассмотрены гиперсингулярные ядра, представляющие аналоги ядра Гильберта и родственные с ним.

В настоящей статье, следуя методике работы [15], на основании формулы обращения классического СИУ с ядром Коши построено в квадратурах

точное решение ГСИУ с ядром (^ — х) 2 на нескольких интервалах. К нему сводятся простейшие задачи о напряженном состоянии линейно -упругой изотропной плоскости с коллинеарной системой трещин. Предполагается, что упругая плоскость находится в состоянии антиплоской или плоской деформации. Берега трещин нагружены симметрически, соответственно, касательными или нормальными силами. Излагается вывод этого определяющего ГСИУ. Одновременно формулируются эквивалентные задачам о трещинах смешанные граничные задачи для упругой полуплоскости. Входящие в решение ГСИУ интегралы вычисляются по квадратурным формулам Гаусса по чебышевским узлам. Более подробно рассматривается частный случай двух симметрично расположенных трещин (интервалов). В этом случае решение ГСИУ разбивается на симметрическую и кососимметрическую части. Они выражаются аналитическими формулами довольно простой структуры. Одновременно методом ортогональных многочленов Чебышева получены точные решения разбираемого ГСИУ для его симметрической и кососимметрической частей. В данном частном случае полученные здесь результаты по сравнению с результатами работы [14] выгодно отличаются своей простотой и компактностью.

Постановка задач и вывод

определяющих уравнений

Рассмотрим следующие простейшие задачи механики трещин. Пусть отнесенное к правой прямоугольной системе координат Охуг линейно-упругое изотропное пространство с модулем сдвига О в плоскости у = 0 содержит систему О из произвольного конечного числа сквозных (туннельных) бесконечно-ленточных трещин (дк:

= {У = 0, ak < х < bk, -œ <z <œj

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

Ty\\y=±0 =Т{Х\\ х eQ , где zyz компонента касауцу=±0 " V/& " &---& &у2

тельных напряжений.

Считается, что под действием таких сил упругое пространство находится в состоянии антиплоской деформации (продольного сдвига) в направлении оси Ог с базовой плоскостью Оху. Тогда единственной отличной от нуля компонентой упругих перемещений будет компонента в направлении оси Ог: = ы2 (х, у), которая в каждой из полуплоскостей у > 0 и у < 0 является гармонической функцией.

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

Ввиду симметрии относительно оси Ох задача о трещинах математически эквивалентна смешанной граничной задаче теории упругости для нижней полуплоскости у < 0:

Au, + = 0, - œ < х < œ; - œ < y < 0,

= r(4 х g L : L = U [ak,bk ]

Приступив к определению функции /(х), обратимся к представлению в виде интеграла Пуассона для нижней полуплоскости у < 0 [16, с. 224, ф-ла (5)]

у р и(5 ^

= 0, х е Ь , Ь = , (1)

^ 0 пРи Х& + У2 ^ Ю

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

Отметим, что мы опять придем к граничной задаче (1), если рассмотрим задачу о системе трещин Ь на горизонтальной оси Ох в кусочно-однородной плоскости Оху, считая, что верхняя полуплоскость у > 0 абсолютно жесткая и касательными силами нагружены только нижние берега трещин.

В граничной задаче (1) требуется определить перемещения на системе отрезков Ь , т.е. функцию

f (х) = и2 (х, у) у=_о, х е Ь.

Предполагается, что f (х) непрерывна на Ь и

обладает производной, удовлетворяющей на Ь условию Гёльдера с некоторым показателем а,

f к ) = Ж ) = 0, к = \\Гп. (2)

Раскрытие трещин и плотности дислокаций на берегах определяются, соответственно, формулами Д( х) = _2 f (х), А & (х) = _2 f& (х), х е Ь.

Коэффициенты интенсивности напряжений и разрушающие напряжения вне системы трещин Ь на оси Ох находятся по формулам из [10].

Отметим, что смешанную граничную задачу (1) можно интерпретировать также следующим образом: каков должен быть режим перемещений f (х)

при антиплоской деформации на системе Ь неза-щемленных отрезков границы нижней упругой полуплоскости у < 0, чтобы на них могли действовать наперед заданные касательные силы интенсивности г(х ).

(х, y )=-y ir^

- œ < х < œ, y < 0,

ж JJ s -х)2 + y2

и(х) =

-¡f (х), хg L,

" 0, х g L &.

Исходя из (3), по закону Гука вычислим касательные напряжения Tyz :

тт = G—2 yz cy

_ G r (s - х)2 ■y

ж L [(s - х)2 + y2

_ ю < х < ю, у < 0. Отсюда предельным переходом у ^ _0 получим

I О г f

Му=_0 =--|Т-у , _Ю< х <Ю ,

Ж Ь _ х)

где интеграл на Ь при 5 = х , х е (ак, Ьк) понимается в смысле Адамара:

Ьк f (х + е)+f (х _е)

г^+0

i fiSf +i

(s - х)2 х+Д* - х)2 Ь

к = 1, n, или в эквивалентной форме

Ы-+0

f(s)ds | ) f(s)ds 2f(х)

(s - х)2 х^С* - х)

k = 1, n.

При помощи (4) реализуем первое граничное условие задачи (1). В результате придем к следующему определяющему ГСИУ поставленной задачи:

■ з I ^

Ж Ь (* _ х)2 О

Его решение должно удовлетворять условиям (2).

После определения функции f (х) найдём касательные напряжения вне системы трещин Ь на оси Ох по уравнению (4) при х е Ь&.

Отметим, что уравнением (5) в соответственных перемещениях и напряжениях описываются с точностью до упругой постоянной аналогичные задачи

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

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

Рассмотрим еще одну смешанную граничную задачу для упругой полуплоскости при плоской деформации, сводящуюся к решению ГСИУ (5). Пусть отнесенная к прямоугольной системе координат 0xy нижняя упругая полуплоскость у < 0

обладает модулем Юнга E, коэффициентом Пуассона V и находится в условиях плоской деформации.

Пусть далее эта полуплоскость на своей границе у = 0 по двум полубесконечным отрезкам

[— го,—а] и [а, да] и по системе отрезков

L=Jk, bk ],

u(x, y )

y =-0 =^y (X

(х, y )

х g [-œ, -a

Tyx(X y) y=-0 =т(х), х g h

y=-0" y=-0

] J [a, œ]J L,

= 0, x g l,

^x ^yx ^ 0 ПРИ x + y •

Здесь l = J[ck,dk],

С = -a,

ck+i = bk

Для вывода определяющего ГСИУ задачи (6) воспользуемся известным решением первой граничной задачи для нижней упругой полуплоскости [6, с. 406-408].

du(x,-0) œ T (s)ds

T (x )=^yx(x, y )| y=-0 =■

u(x,-0) = u(x, y )

-œ<x<œ, 3 =

r(x), x g l,

r0(x), x g R/l, [f(x), x g l,

= и [а

—а < а < К < ••• < а < К < ••• < К < а,

заключенной в отрезке [— а, а], усилена абсолютно жесткими на растяжение-сжатие и абсолютно гибкими на изгиб стрингерами, а на остальной части

отрезка [ а, а] действуют наперед заданные касательные силы интенсивности г(х). В рамках известной модели Мелана [17] описанную задачу можно сформулировать в виде следующей смешанной граничной задачи для нижней упругой полуплоскости:

По формуле обращения Гильберта из (7) получим следующее ключевое уравнение задачи:

T (x)=--к J<

-œ<x<œ.

Это уравнение рассмотрим на системе отрезков I и реализуем второе граничное условие задачи (6). Получим следующее СИУ:

Ж* S - x

g l, f(ck ) = f(dk ) = 0, k = 1, n +1,

= -73 r(x), x g l,

которое после интегрирования по частям дает определяющее ГСИУ обсуждаемой задачи:

ж\\ — х)2

Ж) = Ж) = 0, к = 1Й+1. (9)

ГСИУ (9) с точностью до упругой постоянной совпадает с ГСИУ (5).

Если же ключевое уравнение (8) рассматривать вне системы I, то получим касательные напряжения под стрингерами:

Г. (х)=--^ | , X Е И/,.

^ = а, к = 1, п, = а - система отрезков в [— а,а], свободная от стрингеров; Ы,х,у,Гух горизонтальные перемещения точек полуплоскости и компоненты напряжений.

Требуется определить перемещения граничных точек упругой полуплоскости на системе отрезков

I на оси Ох, т.е. функцию I(х) = и(х, у) у=—0 ,

X Е /, I(Ск) = I(йк) = 0, к = 1пт1, I(х)е с1&2(I).

Ниже займемся решением ГСИУ (5) при условиях (2). После введения безразмерных величин

£ = х/а, п = Vа, (£) = I(а£)/а, g(^) = т(a^)/О придем к следующему ГСИУ:

(р(ар)= 0, (p(ßp)= 0, p = 1,n,

J (ap ß ap = ap/a > ßp = bp/a.

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

Здесь а - характерный параметр длины, например а = |а|, если а ^ 0.

Решение ГСИУ (10)

При помощи интегрирования по частям уравнение (10) можно привести к виду

!j^=-gj, ^gL0.

жL л-ь

Решение СИУ (11) представляется формулой [7, с. 50, ф-лы 244-246].

_ (- 1)n+p = Ш

- S (- Г ÎVn^ Pjg)

ap <J<ßp, p =1,n,

Pn-1(^) = C0 + cj + ... + C^j1 ,

Пп (^)=n^-am| J-ßm| ,

где с ■, 7 = 0, п _ 1, - пока неизвестные коэффициенты многочлена (£).

В (12) £ заменим на Т и по Т произведем интегрирование от а^ до £, а < £ < Рр, Р = 1, п. Получим

((£)= (13)

n+p+1

S (- 1)n - m i^np (ЬлХ/ПЛ) g №n

n-1 J zjdz л/Пn Л)

+ SCj Jj=0 a

, ap , p =1,n

Кпр(£,1)= | , , Т-т.

п п (т)(т_л)

Очевидно, что ((а ) = 0 , р = 1,п . Полагая в (13) £ = Рр и принимая во внимание условие ((р )= 0, для определения постоянных с^ полуS j = Ap), p=1 n C

=1,n C = C-l,

_rj-1dz. (z)

, j =1, n

A;)=S (- 1)n-m+1 J KnP (ßP ,лХ/п;(Л) g (л)л,

Разрешимость СЛАУ (14) может быть показана (12) аналогично тому, как это было сделано в [6, 7].

Так как входящие в (13) и (14) интегралы

Кпр(£,^) и А(п) аналитически не вычисляются,

для их вычисления воспользуемся известными [18] квадратурными формулами типа Гаусса по чебы-шевским узлам, которые после перевода на интервал (а, Р) будут иметь следующий вид:

= жMyf ß-at +ß + a

V(z- a)(ß- z) MS ^ 2 9 2 t = cos((2q -1)ж/ 2M ), q = 1, M,

JV(r-aXß-r) h(z)d&

ß-аЛ ж

• 2 f жт sin I-|X

M +1

Jß-a ß + a x h| —-z + +

z = cos I

m = 1, M .

М +1

Здесь М - любое натуральное число; ^ и Тт -чебышевские узлы, т.е. корни многочленов Чебы-шева первого рода ТМ () и второго рода им (т) соответственно. Полагая

ТЩТ) = ТПрЩт_а/, )(Рр _т) ,

Ппр)(т)= ПТ_«т|т_Рт| ,

m=1,m^ p

-р/ ~ •> ~------г—----------------------- ;

чим следующую систему линейных алгебраических при помощи формул (15) получим уравнений (СЛАУ):

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2 M 1 1

K (ß r(p))=—У

^npXHmi &"m ) дгУ

t( p W p >

(p ) ßp-ap ßp +ap

q 2 q r(p ) = ßp -a

pr ßp + ap 1 m +

A(n ) = Ap =

n-m+1 ,

m+1 tr ^ 2 . ,=.

• 21 —7 sin21 ^

,M +1

nnm)(r)= Пк-a.lk-ß^,

r=1,r Ф m

Knp(^)=—T

Далее рассмотрим важный частный случай, когда L = [— a,—b\\^\\b,a]. Пусть a = |aj . В этом случае будем иметь

а1 = —1, Д = —k; а2 = k, f32 = 1, k = b/a,

Lo =(— 1,—k)U (k,l), П2 (?) = ( — ? № — k2 Л k <?< 1.

В разбираемом частном случае элементы матрицы СЛАУ (14) вычисляются аналитически.

Введём для сокращения дальнейших построений функции = V(1 — rf )р2 — k2 ),

Го( = = V(1 — ?fe2 — k2). Тогда

Lï = frdr

—1Г0 (г) к Гс(0 Далее, воспользовавшись выражением известного интеграла из [19, с. 260, ф-ла 3.153.10], найдём

= К (к&) = К&, к& = ^¡\\-кI, где К (к&) - полный эллиптический интеграл первого рода модуля к&. Обращаясь к , запишем

L(2 ) = L1 О г rdr _ г J1/0(^) f Уа(г)

Приняв во внимание формулы (15) и (16), представим (14) в виде

Г4(2С + т{2С = л{2) 1411 С1 ^ ¿12 С2 = л1 ,

| И2)С + И2)С = Л(2) 1421С1 ^ 422С2 = Л2 ,

где Л(2) = 1

= | [К21(— к п (— п) — К21(— к, п)Е (п)] X

х*Лл\\ 2 ; m +1 У Ы- k (-Г )-K,1(- k ,rf)g (r<1&)]>

sin2[ M7+I № + r?M& + k)

ri" =

Г a22) =

y [a,-r? ))g (-r")-K22(1,r& >)g r ))]> 7

^ M + 1-k

TJV(1 + ri!))W2) + k ).

—r +— 2 7 2

K21 (- *.,) = -f ^

У0 (r)(r + Л) 1 1

V(1+tqi})(k+tf) -q

<щ t «

tq4+v

K22(1,?) = f

M (1+tq2l)(tq2,+k ) (q &-л) &

/0(г)(г-л) 1

t <1) =

t (2) = &q

Решение ГСИУ (10) в случае двух симметрически Этот элементарный интеграл берется подста- расположенных интервалов 40 = (— 1,—к)и (к,1)

новкой 5 =Г . В результате =—Ж/В построим в других аналитических формах более остальных двух случаях простой структуры. С этой целью решение уравне¿(2) = А2) = К& 4-2) = —4,2) =ж/ 2 ния (10) разобьём на симметрическую и кососим21 ? 22 М2 / т"> ^

метрическую части. В случае симметрической чаk

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

Для J£,л, к)

сти g(_£)= g(£) и, следовательно, (_£) = (£), ((_£) = _((£), £е[_1,_к]у[к,1], причем

((± к) = ((±1) = 0. В этом случае формула обращения (12) для соответствующего СИУ (11) преобразуется к виду

г(£)= 1 г + _Со

имеем

J J, л; k )=-^-л -1

п1ж ,^(л), k&l-П

жк0(£) k k<£<1.

Гс(£)

Л > - ,а(л), k&

к <£< 1, ,

где П(р, n, к ) - неполный эллиптический интеграл

Для определения постоянной С0 воспользуемся третьего рода [19, с. 919], откуда

| ( (£^£ = 0, которое после простых 3 (£,±1; к) = 3 (£,1; к ) =

условием J ç

преобразований даёт 1

g0 (z)= g УCt + d ) 1

I0 (k ,z)= J

I J, k )= К & - F k <£< 1.

[k 2П(щ,1, k )+ k2 F Щ, k& )],

щ = arcsin

= arcsin

с =-, d =-.

Далее в (17) £ заменим на и и по и произведем интегрирование от и до £ , к < £ < 1. Получим

((£) = -13(£,); к) g(л)л7o(л)dл+ СоI(£,к), Ж к

к <£< 1, (19)

£ £ I(£,к)=|-^, 3{£,П, к)=( /V 2).

к Уо(и) к Уо(и)\\п _и)

Вычислим интегралы I(£, к) и 3(£,), к) из (19). Имеем

I (£, к ) =

-Г Л _х{Р

к<£<1.

Так как ((1) = 0, то из этих формул вытекает, что постоянная С0, кроме формулы (18), может

быть выражена также формулой 1

со = ,2 , \\/3(1,Пк)g{l)nYo(Л)dl, I (1, к ) = К,

J (1,л, k ) =

k2nl 1,1, k 1+k&2К&

, л = 1.

„V(1 - t2 )(1 - k&2t2 ) J0 V(1 -12 )(1 - k&2t2 )

= F (p, k&),

где F(p, k&) - неполный эллиптический интеграл первого рода [19, с. 918]. В результате

Таким образом, симметрическая часть решения ГСИУ (10) в обсуждаемом частном случае выражается формулами (19)-(22).

Симметрическую часть решения указанного ГСИУ построим также методом многочленов Че-бышева. С этой целью воспользуемся спектральным соотношением [20, 21]

k \\£-л У0(л) n = 0,1,2,..., k <£< 1,

X = cos 3,

= МЛ (X ),

К& J V(u2 -1)(1 - k2u2 )

Y = cos ç,

, k < £, л < 1,

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

K f V(u2 -1)(1 - k2u2 ) &

n = 1,2,...

Коэффициенты у определим при помощи

условий ортогональности многочленов Чебышева второго рода с видоизмененным аргументом:

K & ( 7 nK

K, n = 0.

f^m-1(X )Un-1 (X )

л/1 - X2 dÇ Ytâ)

—, m = n, 2

которые дают

Здесь 1п (.л ) - многочлены Чебышева первого рода аргумента X; 3 при помощи подстановки и = 1/ / в интеграле будет выражаться неполным эллиптическим интегралом первого рода [19, с. 260, ф-ла 10]:

( ( 1,2 ^

Уп = — cth(7nKjK&)gn, n = 1,2,...

gn =f g(ÇU-1(X)dç •

, k& = V 1 - k2 .

V V " У y

Условия ортогональности многочленов Чебышева Tn (X ):

(K&, да = n = 0,

Далее обе части (27) проинтегрируем от к до £, предварительно заменив £ на П . Получим

((Ç)=y0 fvcri +

k /0(л) n=1

У Уп f

Tn (Y)dn

К&/2, т = п * 0, 0, т * п.

Обе части спектрального соотношения (23) продифференцируем по £• После элементарных преобразований получим

к Г0(п)

Положим £ = 1 и воспользуемся формулой (24). В результате находим, что у0 = 0^ Для вычисления второго интеграла из (29) перейдем к переменной ( . Согласно третьей формуле из (23)

—П = — Следовательно,

У0(п) Ж

Tn (Y Ул_ K &

— f cos(n()d( ■ 7 Î

n = 1,2,...

К&) *,(£)

f У0(л)

K& ■ ( п\\ K& sin(n3) .

= — sin(n3) =--1—- sin3:

К& I= К^/1 — X2ип—1(х), п = ¡Д^ тп

В результате

((£) = КлЯ—Х2 У ^иАX),

Ж п=1 п

где коэффициенты уп даются формулами (28). Перейдем к кососимметрической части решения ГСИУ (10) при 4 =[— 1,—к и (к,1). В этом случае g{—£) = —g £), ((—£)=—((£) и ((—£) = ((£), £е т0. Формула обращения (12) с неизвестными коэффициентами уп. Этот ряд для соответствующего СИУ (11) преобразуется к подставим в (26). Воспользовавшись соотношением виду

= -g(ç), k <Ç< 1.

к п —£

Решение СИУ (26) представим в форме бесконечного ряда

((£) =— У упТп(X) , к <£< 1, (27)

/0(£) п=0

(25), получим

„/l _ œ

^-X- Y.VnynU-(X) = g (ç) , k <Ç< 1 •

У0(£) n=1

((Ç) = 2^rf 2 «2

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

В этой формуле £ заменим на и и по и произведем интегрирование от к до £, к <£< 1. В результате получим

ç(£) = - J ln

g {ri)dri +

ç£) = - J j£ л; k Клк Шл + СЛ(£, k ),

—+ arcsin 2

r 2£2 - k2 - Л

V 1 -k2 I

, k<£<1.

,(£, k )=/■

k Ï0(u)

Jl(£,r; k ) = Judu

к 7о(и)()2 _ и2)

Для вычисления интеграла I (£, к) положим в

нем и2 = I. После элементарных выкладок находим

I (£, k ) = 1 arcsin

^2£2 - k2 -1Л ж

к<£<1.

В случае интеграла J1 (£,^; n) воспользуемся,

= 41 — к,2t2 . Тогда dt

Очевидно, что ((к)= 0. С другой стороны, удовлетворив условию ((1) = 0, обнаружим, что С = 0. В итоге симметрическая часть решения ГСИУ (10) для двух симметрических отрезков выражается формулой

((£) = 1 |1п ,У(£п) + у(),£) g())d),

ж к \\г(£,)) _у(п,£)| к <£< 1.

Кососимметрическую часть решения ГСИУ (10) построим также методом ортогональных многочленов Чебышева. С этой целью обратимся к известному соотношению

как и выше, подстановкой u 1

I fln—L- Tn (s d 77" J

Jl(£,r;k) = —j J

k2 (t2 - A2)

ж\\ \\х - s -1

и в нем положим

ln2, n = 0, -1 < х < 1,

Л2 = bf, /(£)=-^л/1—£.

Далее применяем третью подстановку Эйлера

л/Г—t2 = v(1 +1 t = 1—4 ,v = • 11—7

х =-:-, S =---,

k <£, л< 1, £ =4сх + d, л = 4es + d

приводят исходный интеграл к виду

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

следующее:

k&2 (1 - A2 ) J (v2 - B2 )(v2 - B&2 ),

(1 + v2 )dV

ж Jlnb^ T

^2л-k2 -1Л л/л

в2 = Т_А, В& = I, А = 1_). 1 + А В к&2

Вычисляя элементарные интегралы от рациональной функции в правой части, после несложных выкладок окончательно находим

\\£2 -л

f 2£2 - k2 V 1 -k2 I

n = 1,2,...,

n = 0, k <£< 1.

гументом имеют вид

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

y>(Ç)

Ж—1 (5 — х V1 —

Если в его левой части произвести интегрирование по частям, то получим следующее соотношение:

, f ^-М^=-nUn-1 (x ),

^ J (s - x)

n = 1,2,...; -1 <x < 1.

нии получим нужное нам следующее соотношение:

f2\\ - к2 -1 ^

-2- Wo\\)

V 1 к y

^2Ç2 - k2 -1 ^ 1 - k2

, n = 1,2,..., k <Ç< 1.

Условия ортогональности входящих в (36) многочленов Чебышева второго рода с видоизмененным аргументом имеют вид

^ 2Ç2 - k2 -1 ^ v 1 -k 2 ;

X 2Çy0(Ç) dÇ =

m = n, 1 - k 2 & c =.(37)

Далее ГСИУ (10), когда L0 = [- 1,-k]|J (k,l), в

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

—, m = n = 0, 2

— , m = n Ф 0, 4

Обе части спектрального соотношения (34) продифференцируем по £ , а затем по (33) перейдем к переменным х, В результате получим известное соотношение [19, с. 847, ф-ла 7.344.1]

,k<Ç<1.

! \\ -е )2 = 2е

Его решение представим посредством бесконечного ряда

((ç) = y0(ç)y xnUn_,

^2Ç2 - k2 -1 ^ v 1 -k2 ;

к<£<1, с неизвестными коэффициентами х .

Подставляя (39) в (38), меняя порядок интегрирования и суммирования и используя соотношение (36), получим

-y nx„U„ ,

^^ n n-1

^2Ç2 - k2 -1 ^ v 1 -k2 ;

, k<Ç< 1.

Отсюда по условиям ортогональности (37) находим

Интеграл в (35) при s = x понимается в смысле Адамара. Это соотношение, которое может быть интерпретировано как своеобразное спектральное

выражение для ядра 1/|^(s — x)2 ], приведено в [5]. Теперь в (35) по формулам (33) обратно перейдём к переменным rf. После очевидных преобразование

Л=f g «ц ;

n = 1,2,...

Yo(Ç) dÇ,

Последний интеграл можно вычислить по второй квадратурной формуле (15).

Заключение

К решению ГСИУ на нескольких интервалах с классическим гиперсингулярным ядром Коши здесь применены методы СИУ с обычным ядром Коши. К этому ГСИУ, как в работах [18-20], можно применить и методы краевых задач теории аналитических функций. Однако в обоих случаях возникают значительные трудности при аналитическом вычислении нужных сингулярных и обычных интегралов. В плане преодоления этих трудностей эффективнее представляется метод квадратурных формул типа Гаусса, с вычислительной - метод ортогональных многочленов. Оба метода частично использованы в настоящей работе. При помощи квадратурных формул можно провести численный анализ аналитических формул, полученных здесь методами ГСИУ и ортогональных многочленов Чебышева. Изложенные результаты могут быть использованы в смежных исследованиях.

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

Литература

References

ISSN 1026-2237 BULLETIN OF HIGHER EDUCATIONAL INSTITUTIONS. NORTH CAUCASUS REGION. NATURAL SCIENCE. 2020. No. 2

Поступила в редакцию /Received_19 марта 2020 г. /March 19, 2020