ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ, 2023, том 59, № 5, с.658-674
УРАВНЕНИЯ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ
УДК 517.958
О СУЩЕСТВОВАНИИ РЕШЕНИЙ
НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
ДЛЯ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
РАВНОВЕСИЯ ОБОЛОЧЕК ТИПА ТИМОШЕНКО
В ИЗОМЕТРИЧЕСКИХ КООРДИНАТАХ
© 2023 г. С. Н. Тимергалиев
Доказывается существование решений краевой задачи для системы пяти нелинейных диф-
ференциальных уравнений с частными производными второго порядка при заданных нели-
нейных граничных условиях, описывающей состояние равновесия упругих пологих неод-
нородных изотропных оболочек с незакреплёнными краями в рамках сдвиговой модели
Тимошенко, отнесённых к изометрическим координатам. Краевая задача сводится к нели-
нейному операторному уравнению относительно обобщённых перемещений в соболевском
пространстве, разрешимость которого устанавливается с использованием принципа сжа-
тых отображений.
DOI: 10.31857/S0374064123050102, EDN: CYVKAG
1. Постановка задачи. В плоской односвязной ограниченной области Ω рассматривается
система нелинейных дифференциальных уравнений вида
(DT)αλ + DGjλμTλμ + DRj = 0, j = 1, 2,
(DTλμw3αμ )αλ + (DTλ3)αλ + DBλμTλμ + DR3 = 0,
(DM)αλ - DTj3 + DGjλμMλμ + DLj = 0, j = 1, 2,
(1)
при выполнении на границе Γ области Ω условий
D(Tj12/ds - Tj21/ds) = Pj(s), j = 1, 2,
D(T132/ds - T231/ds + T1λw3αλ2/ds - T2λw3αλ1/ds) = P3(s),
D(Mj12/ds - Mj21/ds) = Nj (s), j = 1, 2.
(2)
В (1), (2) и ниже используются следующие обозначения:
Tij ≡ Tij(γ) = Dijknλ-1γλ-1kn, Mij ≡ Mij(γ) = Dijknλγλ-1kn,
γ = (γ01), γk = (γk11k12k13k22k23k33), k = 0,1;
Dijknm = Dijknm(α12) =
Bijkn(α123)(α3)m3, m = 0,2, i,j,k,n = 1,2,3;
-h0/2
2
E
νE
E
B1111 = B2222 =
,
B1122 =
,
B1212 =
,
B1313 = B2323 =
;
12
12
2(1 + ν)
2(1 + ν)
γ0jj = wj - Gλjjwλ - Bjjw3 + w23αj /2, j = 1,2,
γ012 = w1α2 + w2α1 - 2Gλ12wλ - 2B12w3 + w3α1 w3α2 , γ1jj = ψj - Gλjjψλ, j = 1,2,
γ112 = ψ1α2 + ψ2α1 - 2Gλ12ψλ, γ0j3 = w3αj + ψj, j = 1,2, γ033 = γ1k3 0, k = 1,2,3;
(3)
658
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
659
остальные Bijkn равны нулю, αj = αj (s), j = 1, 2, - уравнения кривой Γ, s - длина дуги Γ;
нижний индекс αλ в (1)-(3) и далее означает дифференцирование по αλ, λ = 1, 2.
Если выражения Tij, Mij из (3) подставить в систему (1), то относительно функций w1,
w2, w3, ψ1, ψ2 получим систему пяти дифференциальных уравнений с частными производ-
ными второго порядка, линейных относительно w1, w2, ψ1, ψ2 и нелинейных относительно
w3. Совместно с граничными условиями (2) система (1) описывает состояние равновесия упру-
гой пологой изотропной неоднородной оболочки с незакреплёнными краями в рамках сдвиго-
вой модели Тимошенко [1, с. 168-170, 269], отнесённой к криволинейной системе координат.
При этом Tij - усилия, Mij - моменты; γkij, i, j = 1, 2, 3, k = 0, 1, - компоненты дефор-
маций срединной поверхности S0 оболочки гомеоморфной области Ω; wj, j = 1, 2, и w3 -
соответственно тангенциальные и нормальное перемещения точек S0; ψi, i = 1, 2, - углы
поворота нормальных сечений S0; Bij, i, j = 1, 2, - составляющие тензора кривизны поверх-
ности S0, Gλij - символы Кристоффеля, которые в изометрической системе координат даются
формулами [2, c. 18]
G1jj = (-1)j-1Λα1 /(2Λ), G2jj = (-1)jΛα2 /(2Λ), Gj12 = Gj21 = Λα3-j /(2Λ), j = 1,2;
(4)
Rj, Pj, j = 1,2,3, Lk, Nk, k = 1,2, - компоненты внешних сил, действующих на оболочку;
ν - коэффициент Пуассона, E - модуль Юнга, κ2 - коэффициент сдвига, D = D(α12) =
=
A11A22 - A212 - якобиан (Aij - коэффициенты первой квадратичной формы поверхности
S0; в изометрических координатах A11 = A22 = Λ(α12), A12 = 0 и D = Λ), h0 = const -
толщина оболочки; α1, α2 - декартовы координаты точек области Ω.
В (1)-(3) и в дальнейшем по повторяющимся латинским индексам ведётся суммирование
от 1 до 3, по повторяющимся греческим индексам - от 1 до 2.
Задача (1), (2). Требуется найти решение системы (1), удовлетворяющее граничным усло-
виям (2).
В настоящее время достаточно полно изучена разрешимость краевых задач для диффе-
ренциальных уравнений, описывающих состояние равновесия оболочек в рамках простейшей
модели Кирхгофа-Лява [2-5]. В то же время актуальной задачей является исследование подоб-
ных краевых задач в рамках более сложных моделей теории оболочек, не опирающихся на ги-
потезы Кирхгофа-Лява [2, c. 349]. Имеется ряд работ [6-13], в которых в рамках сдвиговой мо-
дели Тимошенко исследована разрешимость и доказаны теоремы существования обобщённых
решений в соболевских пространствах нелинейных задач для пологих оболочек, отнесённых
к евклидовой системе координат. В основе исследований в [6-13] лежат интегральные пред-
ставления для обобщённых перемещений, содержащие произвольные голоморфные функции,
которые находятся таким образом, чтобы обобщённые перемещения удовлетворяли заданным
граничным условиям.
В данной статье метод работ [6-13] развивается на случай пологих неоднородных изотроп-
ных оболочек типа Тимошенко, отнесённых к изометрической системе координат. Переход к
изометрическим координатам, с одной стороны, расширяет класс рассматриваемых оболочек,
а с другой - существенно усложняет систему дифференциальных уравнений.
Краевую задачу (1), (2) будем изучать в обобщённой постановке. Пусть выполнены следу-
ющие условия:
(a) имеют место включения
Bijkn(α123) (W(1)p(Ω)
Cβ(Ω)) × L1[-h0/2,h0/2], i,j,k,n = 1,2,3,
Bλμ(α12) ∈ W(2)p(Ω), λ,μ,k = 1,2;
(b) якобиан D = Λ(α1, α2) > 0 в Ω и имеет место включение Λ(α1, α2)
p (Ω);
(c) компоненты внешних сил Rj , j = 1, 2, 3, и Lk, k = 1, 2, принадлежат пространству
Lp(Ω), а компоненты Pj, j = 1,2,3, Nk, k = 1,2, - пространству Cβ(Γ); внешние силы
самоуравновешены;
(d) Ω - произвольная односвязная область с границей Γ ∈ C1β.
Здесь и далее 2 < p < 4/(2 - β), 0 < β < 1.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
660
ТИМЕРГАЛИЕВ
Определение. Назовём вектор обобщённых перемещений a = (w1,w2,w312) обоб-
щённым решением задачи (1), (2), если a принадлежит пространству
p (Ω), почти всюду
удовлетворяет системе (1) и поточечно граничным условиям (2).
Обозначим
p (Ω) - пространства Соболева. В силу теорем вложения для соболевских
пространств
p (Ω) c p > 2 обобщённое решение a принадлежит пространству Cα(Ω).
Здесь и везде далее α = (p-2)/p. Заметим, что при 2 < p < 4/(2) справедливо неравенство
α < β/2.
Соотношения для компонент деформаций в (3) для удобства в дальнейших исследованиях
запишем в виде
γkij = eksij + ekcij + χkij, i,j = 1,2,3, k = 0,1,
(5)
где приняты обозначения
e0sjj = wj , e0sj3 = γ0j3, e1sjj = ψj , j = 1,2, e0s12 = w1α2 + w2α1 , e1s12 = ψ1α2 + ψ2α1 ,
e0cjj = -Gλjjwλ - Bjjw3, e1cjj = -Gλjjψλ, j = 1,2, e0c12 = -2Gλ12wλ - 2B12w3,
e1c12 = -2Gλ12ψλ, χ0jj = w23αj /2, j = 1,2, χ012 = w3α1 w3α2 ,
χ1ij = χ0j3 = e0s33 = e1sj3 = ekcj3 0, i,j = 1,2,3, k = 0,1.
(6)
2. Построение интегральных представлений для обобщённых перемещений. Вве-
дём в рассмотрение две комплексные функции:
ωj = ωj(z) = D{D1111j-1(w1α1 + w2α2 ) + D1111j(ψ1α1 + ψ2α2 ) +
+ i[D1212j-1 (w2α1 - w1α2 ) + D1212j(ψ2α1 - ψ1α2 )]}, j = 1, 2, z = α1 +2.
(7)
В системе (1) усилия Tjk, моменты Mjk и компоненты деформаций γnjk заменим их вы-
ражениями из (3), (5). Прибавляя после этого к первому уравнению в (1) второе, умноженное
на мнимую единицу i, а к четвёртому уравнению - пятое, умноженное также на i, систему
(1) при помощи функций ωj(z) из (7) представим в удобной для дальнейших исследований
форме
ωjz + hj(a) = fjc(a) + f(a) - Fj(z), j = 1,2,
DD13130(w3α1α1 + w3α2α2 ) + h3(a) = f3c(a) + f3χ(a) - F3(z), z ∈ Ω,
(8)
где приняты следующие обозначения:
ωjz = (ω1 +2)/2, j = 1,2,
hj(a) = (-1)μ-1[(DD1212j+λ-2)α3 νλ2αμ + i(DD1212j+λ-2)αμ νλ1α3 ] - (j - 1)DD13130(γ013 +023)/2,
ν1j = wj, ν2j = ψj, j = 1,2; h3(a) = (DD13130)αλ w3αλ + (DD13130ψλ)αλ ;
fjc(a) = (fc3j-2 + ifc3j-1)/2, f(a) = (fχ3j-2 + ifχ3j-1)/2, j = 1,2,
f3c(a) = fc3(a), f3χ(a) = fχ3(a), fcj(a) = -(DT(ec))αλ - DGjλμTλμ(e),
fc3+j(a) = -(DM(ec))αλ - DGjλμMλμ(e), fχj(a) = -(DT(χ))αλ - DGjλμTλμ(χ),
fχ3+j(a) = -(DM(χ))αλ - DGjλμMλμ(χ), j = 1,2,
fc3(a) = -DBλμTλμ(e), fχ3(a) = -(DTλμw3αλ )αμ - DBλμTλμ(χ),
F1 = D(R1 + iR2)/2, F2 = D(L1 + iL2)/2, F3 = DR3;
e = es + ec, es = (e0s,e1s), ec = (e0c,e1c), eks = (eks11,eks12,eks13,eks22,eks23,eks33),
ekc = (ekc11,ekc12,ekc13,ekc22,ekc23,ekc33), k = 0,1, χ = (χ011012022);
(9)
eksij, ekcij, χkij определены в (6).
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
661
Отметим, что через e и χ обозначены соответственно линейные и нелинейные части ком-
понент деформации γ, поэтому справедливо представление γ = e + χ.
Аналогично граничные условия (2) запишем в виде
Re [(-i)j tωk(t)] + 2(-1)j DD1212k+δ-2νδ3-jαλλ/ds = ϕc3(k-1)+j (a)(t) +
+ ϕχ3(k-1)+j(a)(t) - F3k+j(s), k,j = 1,2,
DD13130[(w3α2 + ψ2)1/ds - (w3α1 + ψ1)2/ds] = ϕχ3(a)(t) - F6(s),
(10)
где
ϕcj (a)(t) = D[Tj2(ec)1/ds - Tj1(ec)2/ds], ϕc3+j(a)(t) = D[Mj2(ec)1/ds -
- Mj1(ec)2/ds], ϕχj(a)(t) = D[Tj2(χ)1/ds - Tj1(χ)2/ds],
ϕχ3+j (a)(t) = D[Mj2(χ)1/ds - Mj1(χ)2/ds], j = 1, 2, ϕc3(a)(t) 0,
ϕχ3(a)(t) = D[(T11(γ)w3α1 + T12(γ)w3α2 )2/ds - (T22(γ)w3α2 + T12(γ)w3α1 )1/ds],
F3+j = -Pj, j = 1,2, F6(s) = P3(s), F6+k = -Nk, k = 1,2,
(11)
усилия Tjk, моменты Mjk определены в (3).
В основе исследования системы уравнений (8) при граничных условиях (10) лежат инте-
гральные представления для обобщённых перемещений wj, j = 1, 2, 3, ψk, k = 1, 2.
Для их вывода рассмотрим уравнения
ωjz = ρj, j = 1,2, DD13130(w3α1α1 + w3α2α2 ) = ρ3,
(12)
где ρ1 = ρ1 +2, ρ2 = ρ4 +5, ρ3 = ρ3 - произвольно фиксированные функции, принадле-
жащие пространству Lp(Ω).
Первые два уравнения в (12) представляют собой неоднородные уравнения Коши-Римана.
Их общие решения даются формулами [14, c. 29]
ωj(z) = Φj(z) +j(z) ≡ ωjj;ρj)(z),
∫∫
1
ρj(ζ)
j(z) = -
dξ dη, j = 1, 2, ζ = ξ + iη,
(13)
π
ζ-z
Ω
где Φj(z) - произвольные голоморфные функции, принадлежащие пространству Cα(Ω).
Известно [14, c. 39-41, 46], что T - вполне непрерывный оператор в пространствах Lp(Ω)
и C(Ω), отображающий их в пространства Cα(Ω) и Ck+1α(Ω) соответственно. Кроме того,
существуют обобщённые производные
∫∫
∂Tf
∂Tf
1
f (ζ)
=f,
≡ Sf = -
dξ dη,
(14)
∂z
∂z
π
(ζ - z)2
Ω
где S - линейный ограниченный оператор в пространствах Lp(Ω), p > 1, и C(Ω).
Представления (13) в свою очередь при помощи функций ω01 = w2 + iw1, ω02 = ψ2 +1
запишем в виде неоднородных уравнений Коши-Римана
ω0j¯z = i(d2j-1[ω1] + d2j[ω2]) ≡ iTjω, j = 1,2, ω = (ω12),
(15)
общие решения которых имеют вид
ω0j(z) = Ψj(z) + iTTjω(z) ≡ ω0jj;ω)(z), j = 1,2.
(16)
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
662
ТИМЕРГАЛИЕВ
В (15), (16) приняты обозначения
d2j+λ-2[ωλ] = d12j+λ-2ωλ + (-1)j+λd22j+λ-2ωλ, j,λ = 1,2,
)
)
1
(D1111
D12124-2k
1
(D12121
D11111
4-2k
dj3k-2 =
+ (-1)j
,
dj2 = dj3 =
+ (-1)j
,
k,j = 1,2,
4D δ0
δ1
4D δ1
δ0
δ0 = D11110D11112 - (D11111)2, δ1 = D12120D12122 - (D12121)2;
(17)
Ψj(z) ∈ C1α(Ω) - произвольные голоморфные функции.
Третье уравнение в (12) запишем в виде
w3zz = ρ3/4,
ρ3 = ρ3/(DD13130), w3z = (w3α1 - iw3α2 )/2,
откуда получим
∫∫
1
z
w3(z) = Re Ψ3(z)
T ρ3 ≡ w33;ρ3)(z),
3 =-
ρ3(ζ)ln
ξ dη,
(18)
1
d
2π
ζ
Ω
где Ψ3(z) ∈ C1α(Ω) - произвольная голоморфная функция.
Соотношения (16), (18) являются искомыми интегральными представлениями для обоб-
щённых перемещений. Для их частных производных первого и второго порядков при помощи
формул (13)-(18) и (8.20) из [14, c. 58] получаем
νjkαk = Im[ω0jz - (-1)kω0jz], νjkαn = Re[ω0jz + (-1)kω0jz], k = n, j,k,n = 1,2,
ω0jz = Ψ′j(z) + iSTjω(z), ω0jz = iTjω, w3αj = 2Re (ij-1w3z), j = 1,2,
w3z = Ψ3(z)/2 + T ρ3(z)/4, νknαjαj = -Re{in[ω0kzz + (-1)j(ω0kzz + ω0kzz)]},
νknα1α2 = Re{in-1(ω0kzz - ω0kzz)}, w3αjαj = 2[w3zz + (-1)j-1Rew3zz], k,n,j = 1,2,
w3α1α2 = -2Im w3zz, ωkzz = Tk1ω + Sk10;ρ0), ωkzz = Tk2ω + Sk20;ρ0),
1
Tkω(τ)
ω0kzz = Ψ′′k(z) +0
(z) -
dτ , k = 1, 2, Φ0 = (Φ1, Φ2), ρ0 = (ρ1, ρ2),
kζζ
2π
(τ - z)2
Γ
w3zz = Ψ′′3(z)/2 +3/4, w3zz = ρ3/4, Tjkω = i[d12j+μ-2,kωμ + (-1)j+μd22j+μ-2,kωμ],
Sjk0;ρ0) = i[d12j+μ-2ωμ,k + (-1)j+μd22j+μ-2ωμ,3-k], ωj,1 ≡ ωjz = Φ′j(z) +j(z),
ωj,2 ≡ ωjz = ρj, djm,1 ≡ djmz, djm,2 ≡ djmz, j,k = 1,2, m = 1,4.
(19)
3. Решение задачи (1), (2). Интегральные представления (16), (18) для обобщённых
перемещений a = (w1, w2, w3, ψ1, ψ2) содержат произвольные голоморфные функции Φj(z),
j = 1,2, Ψk(z), k = 1,2,3, и произвольные функции ρj(z), j = 1,2,3. Их найдём так, чтобы
обобщённые перемещения удовлетворяли системе (8) и граничным условиям (10), при этом
правые части уравнений (8) и граничных условий (10) временно считаем известными. С этой
целью соотношения (16), (18), (19) подставим в левые части системы (8) и граничных условий
(10). В результате система уравнений (8) запишется в виде
ρj(z) + hj1(ρ)(z) + hj2(Φ)(z) = fjc(a)(z) + f(a)(z) - Fj(z), j = 1,2,3, z ∈ Ω,
(20)
где через hj1(ρ)(z) и hj2(Φ)(z) обозначены те части выражения оператора hj (a) в (9), которые
содержат функции ρ = (ρ1, ρ2, ρ3) и Φ = (Φ1, Φ2, Ψ1, Ψ2, Ψ3) соответственно.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
663
Граничные условия (10) с учётом представлений
S(Tj Φ0)+(t) = -(t)2[d12j-1(t1(t) + d12j (t2(t)] + K0j0)(t), Φ0 = (Φ1, Φ2),
d12j+μ-2(t)
ψ(τ, t) - ψ(t, t)
d22j+μ-2(t)
K0j0)(t) = -
Φμ(τ)dτ - (-1)j+μ
×
2πi
τ-t
2πi
Γ
∫∫
ψ(τ, t) - ψ(t, t)
1
d12j+μ-2(ζ) - d12j+μ-2(t)
×
Φμ(τ) dτ -
Φμ(ζ)dξ dη -
τ-t
π
(ζ - t)2
Γ
Ω
∫∫
j+μ
d22j+μ-2(ζ) - d22j+μ-2(t)
(-1)
-
Φμ(ζ)dξ dη, j = 1,2,
π
(ζ - t)2
Ω
ψ(τ, t) = (τ - t)/(τ - t), ψ(t, t) = (t)2,
(21)
получаемых при помощи соотношений (13)-(15), формул (4.7), (4.9) из [14, c. 28] и формул
Сохоцкого [15, c. 66], преобразуются к виду
(-1)j d(t)Re [ij tΦλ(t)] - 2DD1212λ+k-2(t)Re [ij-1tΨ′λ(t)] - 2DD1212λ+k-2(t)Re [ij tK0λ0)(t)] +
+ H3(k-1)+jρ(t) = ϕc3(k-1)+j(a)(t) + ϕχ3(k-1)+j(a)(t) - F3k+j(s), k,j = 1,2,
DD13130(t)Re [itΨ3(t)] + K03(Φ)(t) + H3ρ(t) = ϕχ3(a)(t) - F6(s),
(22)
где приняты следующие обозначения:
H3(k-1)+jρ(t) = Re [(-i)jtk(t)] - 2DD1212k+λ-2(t)Re {ijt(I + S)(Tλ0)+(t)}, k,j = 1,2,
H3ρ(t) = DD13130(t)Re [it(T ρ3(t)/2 + TT20(t))],
K03(Φ)(t) = DD13130(t)Re {t2(t) + iTT2Φ0(t)]},
dkj(t) = (-1)j-1[2(-1)λDD1212λ+k-2(t)d22λ+j-2(t) + 3 - k - j], k,j = 1,2,
(23)
I - тождественный оператор, операторы T, S, Tλ и функции dkj(t) определены в (13), (14),
(15) и (17) соответственно; Φλ(t) Φ+λ(t), t ∈ Γ; символ Φ+λ(t) здесь и далее означает предел
функции Φλ(z) при z → t ∈ Γ изнутри области Ω.
Таким образом, для определения функций ρj ∈ Lp(Ω), j = 1, 2, 3, Φk(z) ∈ Cα(Ω), k =
= 1, 2, Ψj(z) ∈ C1α(Ω), j = 1, 2, 3, получили систему уравнений (20), (22). Голоморфные
функции будем искать в виде интегралов типа Коши с действительными плотностями
Φk(z) = Θ(μ2k)(z) Φk(μ2k)(z), k = 1,2,
1
f (τ)
Ψ′j(z)=i(j-1)(j-2)/2Θ(μ2j-1)(z)Ψ′j(μ2j-1)(z), j=1,2,3, Θ(f)(z)=
,
(24)
2πi
τ(τ -z)
Γ
где μj(t) ∈ Cα(Γ), j = 1, 5, - произвольные действительные функции, τ = dτ/dσ, dσ -
элемент длины дуги кривой Γ.
Для функций Ψj(z), j = 1, 2, 3, имеем представления
Ψj(z) = i(j-1)(j-2)/2Θ0(μ2j-1)(z) + c2j-1 + ic2j Ψj(μ2j-1)(z) + c2j-1 + ic2j, j = 1,2,3,
(
)
1
f (τ)
z
Θ0(f)(z) = -
ln
1-
dτ,
(25)
2πi
τ
τ
Γ
7
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
664
ТИМЕРГАЛИЕВ
где cj , j = 1, 6, - произвольные действительные постоянные, под ln(1 - z/τ) понимается
однозначная ветвь, обращающаяся в нуль при z = 0.
Используя формулы Сохоцкого [15, c. 66], находим Φk(t), k = 1, 2, Ψ′j (t), j = 1, 2, 3, t ∈ Γ.
Подставляя их выражения, а также представления (25) в систему (20), (22), после несложных
преобразований приходим к следующей системе уравнений относительно функций ρ ∈ Lp(Ω)
и μ = (μ12345) ∈ Cα(Γ):
ρj(z) + hj1(ρ)(z) + hj2(μ)(z) = fjc(a)(z) + f(a)(z) + gjc(z) - Fj(z), z ∈ Ω, j = 1,2,3,
]
∑[
μn(τ)
ajn(t)μn(t) + bjn(t)
+ Kjμ(t) + Hjρ(t) =
τ-t
n=1
Γ
= ϕcj(a)(t) + ϕχj(a)(t) + g3+jc(t) - F3+j(t), j = 1,5, t ∈ Γ,
(26)
в которой приняты обозначения
K3(n-1)+jμ(t) = (-1)jd(t){Re [ijtΘ(μ2λ)(t)] - iRe(ij-1)Θ(τμ2λ)(t)} +
+ 2DD1212λ+n-2(t){Re [ij+1tΘ(μ2λ-1)(t)] - iRe (ij )Θ(τμ2λ-1)(t) - Re [ij tK0λ(μ0)(t)]}, n, j = 1, 2,
K3μ(t) = K03(μ)(t) - DD13130(t)Re [tΘ(μ5)(t)];
g2c(z) = DD13130(c4 + ic3)/2, g3c(z) = -c4(DD13130)α1 - c3(DD13130)α2 ,
g6c(t) = DD13130(t)(c42/ds - c31/ds), g1c(z) = g3+jc(t) 0, j = 1,2,4,5;
a3(k-1)+j,2λ(t) = (-1)jd(t)Re (ij)/2, b3(k-1)+j,2λ(t) = (-1)jd(t)Re (ij-1)/(2π),
a3(k-1)+j,2λ-1(t) = -DD1212λ+k-2(t)Re (ij-1), b3(k-1)+j,2λ-1(t) = DD1212λ+k-2(t)Re (ij)/π,
k,j,λ = 1,2, a35(t) = -DD13130(t)/2,
(27)
остальные ajk, bjk равны нулю; здесь hj2(μ)(z) ≡ hj2(Φ(μ))(z), K0j (μ0)(t) ≡ K0j0(μ0))(t),
j = 1, 2; K03(μ)(t) ≡ K03(Φ(μ))(t), Φ(μ) = (Φ1(μ2), Φ2(μ4), Ψ1(μ1), Ψ2(μ3), Ψ3(μ5)), μ0 =
= (μ2, μ4).
Лемма 1. Пусть выполнены условия (a), (b), (c), (d). Тогда:
1) hj1(ρ), j = 1, 2, 3, - линейные вполне непрерывные операторы в пространстве Lp(Ω);
2) hj2(μ), j = 1, 2, 3, - линейные вполне непрерывные операторы из Cν (Γ) для любого
ν ∈ (0,1) в Lp(Ω);
3) Kjμ, j = 1, 5, - линейные вполне непрерывные операторы из Cν (Γ) для любого ν ∈
(0, 1) в Cγ(Γ) при всех γ < β/2;
4) Hjρ, j = 1, 5, - линейные вполне непрерывные операторы из Lp(Ω) в Cα (Γ) при всех
α < α и ограниченные операторы из Lp(Ω) в Cα(Γ);
5) имеют место включения
fjc (a)(z), fχ(a)(z), Fj (z), gc (z) ∈ Lp(Ω), j = 1, 2, 3, ϕcj (a)(t),
ϕχj (a)(t) ∈ Cα(Γ), F3+j (t), gc(t), ajk(t), bjk(t) ∈ Cβ(Γ), j, k = 1, 5.
Доказательство. Известно [14, c. 26-27], что интеграл типа Коши θ(f) в (24) представля-
ет собой ограниченный оператор из Cα(Γ) в Cα(Ω), а его производная θ(f) - ограниченный
оператор из Cα(Γ) в Lq(Ω),
1 < q < 2/(1 - α). Кроме того, нетрудно показать, что θ(f) -
вполне непрерывный оператор из Cα(Γ) в Lp(Ω) при всех p > 1 и в Cα (Ω) для любого
α < α. Учитывая это, а также свойства операторов T, S, определённых в (13), (14), ис-
пользуя представления для производных первого порядка обобщённых перемещений в (19)
и выражения для операторов hj (a) в (9), получаем, что первые утверждения 1), 2) леммы
справедливы.
Так как ψ(τ, t) ∈ Cβ(Γ) × Cβ(Γ) [15, с. 28-32], d(t) ∈ Cβ(Γ), dkj(z),D13130(z),D(z)
∈ Cβ(Ω), то, принимая во внимание следствие 4.3 из [16, с. 124], легко убеждаемся в том,
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
665
что первые два слагаемых в правой части представления для оператора K0j (μ0) в (21) суть
вполне непрерывные операторы из Cν (Γ), ν ∈ (0, 1) в Cγ(Γ) для любого γ < β. Также
нетрудно показать, что третье и четвертое слагаемые этого представления в (21) суть вполне
непрерывные операторы из Cν (Γ), ν ∈ (0, 1) в Cγ(Γ), γ < β. Тогда получаем, что Koj (μ0),
j = 1,2, - линейные вполне непрерывные операторы из Cν(Γ), ν ∈ (0,1) в Cγ(Γ), γ < β.
Аналогично из представления оператора K03(μ) в (23) следует, что K03(μ) - линейный вполне
непрерывный оператор из Cν (Γ) в Cβ(Γ) для всех ν ∈ (0, 1).
Далее два первых слагаемых в правой части формулы для оператора K3(n-1)+j μ в (27)
преобразуем к виду
(
)
}
(-i)j
{ (-1)j - 1
μ2λ(τ) τ - t
(-1)j
μ2λ(τ)
t
d(t)
+
Im
2
2πi
τ
τ-t
π
τ
τ-t
Γ
Γ
Следовательно, с учётом включений τ, d ∈ Cβ(Γ) и равенства
Im[τ/(τ - t)] = k(τ,t)/|τ - t|1-β/2,
где k(τ, t) ∈ Cβ/2(Γ) × Cβ/2(Γ) [15, c. 31-32, 55-56], а также следствий 4.4,4.5 из [16, с. 125]
получаем, что эти слагаемые в выражении для оператора K3(n-1)+j μ в (27) определяют линей-
ный вполне непрерывный оператор из Cν (Γ) в Cγ(Γ) для любых ν ∈ (0, 1) и γ < β/2. Анало-
гично показываем, что третье и четвертое слагаемые в выражении для оператора K3(n-1)+j μ
в (27) обладают таким же свойством. Тогда из представлений операторов Kjμ, j = 1, 5, в (27)
вытекает справедливость утверждения 3) леммы. Справедливость утверждения 4) следует из
представлений операторов Hjρ, j = 1, 5, в (23) с учётом свойств операторов T, S, интеграла
типа Коши и соотношений
)
(
1
1
Tλ0(τ)
S(TλT ρ0)+(t) = T
Tλ0
(t) -
(t)2TλT ρ0(t) -
dτ , λ = 1, 2,
∂ζ
2
2πi
τ-t
Γ
которые получаются с использованием формул (8.20) из [14, c. 58] и формул Сохоцкого. Спра-
ведливость утверждения 5) леммы непосредственно вытекает из формул (9), (11), (27). Лемма
доказана.
Исследуем разрешимость системы уравнений (26) в пространстве Lp(Ω) × Cα (Γ), α < α.
Заметим, что любое решение (ρ, μ) ∈ Lp(Ω)×Cα (Γ) системы (26) в силу леммы 1 принадлежит
пространству Lp(Ω) × Cα(Γ). Используя выражения для ajk(t), bjk(t) из (27), вычисляем
определитель
det [A(t) - πiB(t)] = D3D13130δ1/(32δ0)(a21 - a0a2), an = D1111n + D1122n, n = 0, 1, 2,
где δ0, δ1 определены в (17), а A = (ajk), B = (bjk) - квадратные матрицы пятого порядка.
Итак, det [A(t) - πiB(t)] = 0 на границе Γ и для индекса системы (26) получаем
[
]
1
det (A - πiB)
χ=
arg
=0
2π
det (A + πiB)Γ
(здесь символ [arg ϕ]Γ означает приращение аргумента функции ϕ при обходе кривой Γ один
раз в положительном направлении). Следовательно, к системе (26) применима альтернатива
Фредгольма. Пусть (ρ, μ) ∈ Lp(Ω) × Cα (Γ) - решение системы (26) при нулевой правой части.
Этому решению по формулам (24), (25) с постоянными cj = 0, j = 1, 6, соответствуют голо-
морфные функции Φk(z), Ψj(z), которые в свою очередь по формулам (16), (18) определяют
функции wj , j = 1, 2, 3, ψk, k = 1, 2. Эти функции, как нетрудно видеть, удовлетворяют
однородной системе линейных уравнений (8) (
c
+ fχ - Fj 0, j = 1,2,3) и однородным
линейным граничным условиям (10) (ϕcj + ϕχj - F3+j 0, j = 1, 5). Умножим действи-
тельную и мнимую части первого уравнения однородной системы (8) соответственно на w1
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
7
666
ТИМЕРГАЛИЕВ
и w2, второго уравнения - соответственно на ψ1 и ψ2, а третье уравнение - на w3. После
этого проинтегрируем по области Ω и сложим получившиеся равенства. С учётом однородных
граничных условий (10) получаем, что wj, j = 1, 2, 3, ψk, k = 1, 2, удовлетворяют системе
νj1α1 = 0, νj2α2 = 0, νj1α2 + νj2α1 = 0, w3αj + ψj = 0, j = 1,2, решение которой имеет вид
w1 = -c0α2 + c1, w2 = c0α1 + c2, w3 = -c4α1 - c5α2 + c6, ψ1 = c4, ψ2 = c5,
(28)
где cj - произвольные действительные постоянные.
Так как Ψj (0) = 0, j = 1, 2, 3, w3(0) = 0, то из (28) будем иметь w1 = -c0α2 + c1, w2 =
= c0α1 + c2, w3 = ψ1 = ψ2 0. Тогда ωj(z) = 2ic0DD1212j-1, j = 1,2, и из уравнений (12)
следуют равенства
ρj(z) = 2ic0(DD1212j-1)z, j = 1,2, ρ3(z) 0, z ∈ Ω.
(29)
Используя формулы (13), (16), (18) и представление для ω0jz в (19), находим функции
Φk(z), k = 1,2, Ψ′j(z), j = 1,2,3, подставив которые в (24), получим
μ1(t)/t - c0(t)2 = F-1(t), μ2j(t)/t - 2ic0DD1212j-1(t) = F-2j(t), j = 1,2,
μ2j-1(t)/t = F-2j-1(t), j = 2,3,
где F-j(t) - граничные значения функции F-j(z), голоморфной во внешности Ω и исчезающей
на бесконечности. Следовательно, для функции F-j(z) во внешности области Ω приходим к
задаче Римана-Гильберта с краевым условием Re [itF-j(t)] = f-j(t), j = 1, 5, где
f-1(t) = c0Re(it), f-2j(t) = 2c0DD1212j-1(t)Re t, j = 1,2, f-2j-1(t) = 0, j = 2,3.
Используя решение этой задачи [17, c. 253], для функций μj(t) получаем представления
μj(t) = c0μ0j(t) + β0jμ1j(t), j = 1,2,4, μj(t) = β0jμ1j(t), j = 3,5,
(30)
где μkj(t) - известные действительные функции, принадлежащие пространству Cα(Γ); c0,
β0j - произвольные действительные постоянные.
Решения (29), (30) показывают, что однородная система уравнений (26) имеет шесть ли-
нейно независимых решений. Тогда союзная с ней система уравнений также будет иметь
шесть линейно независимых решений. Для вывода союзной системы действительные и мни-
мые части левых частей уравнений в (20) умножим соответственно на действительные функ-
ции v1, v2, v3, v4, v5 ∈ Lq(Ω),
1/p + 1/q = 1, и проинтегрируем по области Ω, а левые части
уравнений в (22) умножим на действительные функции ν1, ν2, ν3, ν4, ν5 ∈ Cα(Γ) и проинте-
грируем по кривой Γ. После этого их сложим и приравняем к нулю. Заменяя голоморфные
функции Φj(z), Ψk(z), Ψ′k(z) их выражениями из (24), (25) с постоянными, равными нулю,
переставляя порядок интегрирования в полученных повторных интегралах, после несложных,
но достаточно громоздких преобразований приходим к искомой союзной системе уравнений
vj(z) - T3+jv(z) + 2Θ(τνj)(z) = 0, j = 1,2, Re T3v(z) = 0, z ∈ Ω,
Re{i[T3+jv(t) --(τνj)(t)]} = 0, j = 1,2, Re [Tg(v)(t) + Θ-(τDD13130ν3)(t)] = 0,
Re {T[(DD1212λ+j-2)ζvλ](t) --(τDD1212λ+j-2νλ)(t) + (j - 1)[iT0g(v)(t)-
-T(DD13130τν3)(t)]} = 0, j = 1,2, t ∈ Γ;
vj = v3j-2 + iv3j-1, νj = ν3j-2 +3j-1, j = 1,2, v3 = v3, ν3 = ν3.
(31)
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
667
В уравнениях (31) приняты обозначения
T3v(z) = -2Tg(v)(z) + 2DD13130(z)v3(z) - 2Θ(τDD13130ν3)(z), T3+jv(z) =
= 2T dj+2λ-2[Sλv](z) + T d2+j [T3v](z), Sj v(z) = S[(DD1212j+λ-2)ζ vλ](z) - (DD1212j+λ-2)zvλ(z) -
-(τDD1212j+λ-2νλ)(z), j = 1, 2, g(v)(z) = (DD13130)z(z)v3(z) - DD13130(z)v2(z)/4,
∫∫
(
)
(
)
1
ζ
1
τ
T0f(z) = -
f (ζ) ln
1-
dξ dη, Tf(z) = -
f (τ) ln
1-
dσ,
πi
z
2πi
z
Ω
Γ
1
f (τ)
Θ(f)(z) =
,
v = (v1,v2,v3,v4,v5),
(32)
2πi
τ(τ - z)2
Γ
Θ-(f)(t) - граничные значения функции Θ(f)(z) при z → t ∈ Γ извне области Ω; операторы
Tf, Sf, dj[f], Θ(f) определены в (13), (14), (17), (24) соответственно.
Система (31), как отмечено выше, имеет шесть линейно независимых решений. Получим их
явные выражения. Далее в (31) под v ∈ Lq(Ω), 1/p+1/q = 1, ν ∈ Cα(Γ) будем подразумевать
некоторое её решение.
Заметим, что операторы T, T0, T, введённые в (13), (32), определяют функции T f(z),
T0f(z), Tf(z), которые голоморфны во внешности области Ω и обращаются в нуль на бес-
конечности. Этими же свойствами обладает и функция θ(f)(z). Поэтому последние пять ра-
венств на кривой Γ в (31) представляют собой краевые условия задачи Римана-Гильберта с
нулевым индексом для функций, голоморфных вне Ω и исчезающих на бесконечности. Такая
задача, как известно, имеет только нулевое решение. Следовательно, эти пять равенств на
кривой Γ преобразуются к виду
T3+jv(z) - 2Θ(τνj)(z) = 0, j = 1,2, Tg(v)(z) + Θ(τDD13130ν3)(z) = 0,
T [(DD1212λ+j-2)ζ vλ](z) - 2Θ(τDD1212λ+j-2νλ)(z) + (j - 1)[iT0g(v)(z) - T(DD13130τν3)(z)] = 0,
j = 1,2, z ∈ Ω1C \ Ω,
(33)
C - комплексная плоскость.
Из первых трёх равенств в (31) следует, что функции vj , j = 1, 5, принадлежат простран-
ству
q1
(Ω)
Cα(Ω), 1 < q1 < 2/(1 - α). В них перейдём к пределу при z → t ∈ Γ изнутри
области Ω, а в первых трёх равенствах в (33) - извне области Ω, затем последние прибавим
к первым трём соответственно. Принимая во внимание непрерывность функций вида T f(z)
при f ∈ Lp(Ω) на C и используя формулы Сохоцкого, получаем
vj(t) = -2νj(t), j = 1,2, v3(t) = ν3(t), t ∈ Γ.
(34)
Продифференцировав первые два равенства в (31) по z, с учётом (14) получим равенства
vjz = 2dj+2λ-2[Sλv](z) + d2+j[T3v](z), j = 1,2, z ∈ Ω,
откуда, рассматривая их как систему относительно X1 = 2S1v, X2 = 2S2v + T3v и решая её,
будем иметь
Xj = D[(D1111j+λ-2 - D1212j+λ-2)vλz + (D1111j+λ-2 + D1212j+λ-2)vλz], j = 1,2, z ∈ Ω.
(35)
Пусть дополнительно выполнены условия
D1212j, j = 0,1,2, D13130 ∈ W(2)p(Ω).
(36)
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
668
ТИМЕРГАЛИЕВ
Используя соотношения для функций T3v(z), Sjv(z), j = 1, 2, в (32), находим Xjz, j =
= 1, 2, которые, как нетрудно видеть, принадлежат пространству Lq1 (Ω), 1 < q1 < 2/(1 - α).
Теперь эти выражения Xjz, j = 1, 2, подставим в левые части соотношений, полученных
дифференцированием по z равенств (35). Третье равенство в (31) дифференцируем по z и
z. При помощи несложных преобразований полученных соотношений убеждаемся в том, что
вектор-функция v = (v1, v2, 2v3, v4, v5) является решением системы линейных уравнений (8)
при нулевой правой части.
Далее от решения (v, ν) союзной системы уравнений (31) потребуем, чтобы ν(t) ∈ C1α(Γ).
Тогда, как нетрудно видеть, v(z) ∈ C1α(Ω). Теперь в равенствах (35) переходим к пределу
при z → t ∈ Γ изнутри области Ω, при этом левую часть X+j(t) заменим выражением,
полученным с использованием представлений (Sjv)(z), T3v(z) в (32). Затем из них вычтем
соответственно равенства, которые получаются дифференцированием по z последних двух
соотношений в (33) с последующим переходом в них к пределу при z → t ∈ Γ извне Ω.
Далее третьи равенства в (31) и (33) продифференцируем по z, в получившихся равенствах
перейдём к пределу при z → t ∈ Γ соответственно изнутри и извне области Ω и затем вычтем
их друг из друга. При помощи полученных таким образом равенств на кривой Γ, используя
соотношения (34) и формулы
(Sf)+(t) - (Sf)-(t) = -f(t) · (t)2, θ+(τf)(t) - θ′-(τf)(t) = ft + ft · (t)2, t ∈ Γ,
в которых операторы Sf, Θ(f) определены в (14), (32), после несложных преобразований
приходим к тому, что функции v1, v2,
2v3, v4, v5 удовлетворяют также и однородным
линейным граничным условиям в (10). Таким образом, вектор v = (v1, v2, 2v3, v4, v5) являет-
ся решением однородной системы линейных уравнений в (8), удовлетворяющим однородным
линейным граничным условиям в (10). Следовательно, в соответствии с (28) для компонент
вектора v получим следующие представления:
v1 = -c0α2 + c1, v2 = c0α1 + c2, v3 = (-c4α1 - c5α2 + c6)/2, v4 = c4, v5 = c5,
где cj - произвольные действительные постоянные.
Функции νj(t) и vk связаны друг с другом формулами (34). Следовательно, решение
(v, ν)T , v = (v1, v2, v3, v4, v5), ν = (ν1, ν2, ν3, ν4, ν5) союзной системы (31) можно представить
как (v, ν)T = c0γ1 + c1γ2 + c2γ3 + c4γ4 + c5γ5 + c6γ6, где γk = (γk1, γk2, . . . , γk10), k = 1, 6, -
линейно независимые решения системы (31). Тогда для разрешимости системы (26) необходимо
и достаточно, чтобы выполнялись условия
∫∫
{Re [(f1c + f1χ + g1c - F1)(z)(γk1 - iγk2)(z) + (f2c + f2χ + g2c - F2)(z)(γk4 - iγk5)(z)] +
Ω
+(f3c+f3χ+g3c-F3)(z)γk3(z)} dα12+
(ϕcj +ϕχj +g3+jc -F3+j )(t)γk,5+j (t) ds = 0, k = 1, 6,
j=1Γ
которые после несложных преобразований принимают вид
∫∫
∫∫
DRj12 + Pjds +
DGjλμTλμ(γ)12 = 0, j = 1,2,
Ω
Γ
Ω
∫∫
∫∫
D(R1α2 - R2α1)12 + (P1α2 - P2α1)ds +
D(α2G1λμ - α1G2λμ)Tλμ(γ)12 = 0,
Ω
Γ
Ω
∫∫
∫∫
D(αj R3 - Lj )12 + (αj P3 - Nj)ds +
jBλμTλμ(γ)12 -
Ω
Γ
Ω
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
669
∫∫
∫∫
− DTλj(γ)w3αλ12 - DGjλμMλμ(γ)12 = 0, j = 1,2,
Ω
Ω
∫∫
∫∫
DR312 + P3ds +
DBλμTλμ(γ)12 = 0,
(37)
Ω
Γ
Ω
где Rj, Pj, j = 1, 2, 3, Lk, Nk, k = 1, 2, - компоненты внешних сил, γ - произвольно
фиксированный вектор деформации, w3 - произвольно фиксированная функция.
При выполнении условий (37) общее решение системы (26) можно представить в виде
(ρ, μ) = (ρc, μc)(a) + (ρχ, μχ)(a) + (ρ, μ) + (ρF , μF ), (ρc, μc)(a) = Rfc(a),
(ρχ, μχ)(a) = Rfχ(a), (ρ, μ) = Rgc + (ρ, μ), (ρF , μF ) = -RF,
(38)
где fc(a) = (f1c, f2c, f3c, ϕc1, . . . , ϕc5), fχ(a) = (f1χ, f2χ, f3χ, ϕχ1, . . . , ϕχ5), gc = (g1c, . . . , g8c), F =
= (F1, . . . , F8); R = (R1, . . . , R8); Rj , j = 1, 2, 3, и Rk, k = 4, 8, - линейные ограниченные
операторы из Lp(Ω) × Cα(Γ) в Lp(Ω) и в Cα(Γ) соответственно; функции ρ = (ρ1, ρ2, ρ3),
μ=(μ1,...,μ5) определены формулами(29),(30),а
fjc, fχ, ϕck, ϕχk, gc, Fn - формулами
в (9), (11), (27).
Если выражение для вектор-функции μ(t) из (38) подставить в соотношения (24), (25), то
для голоморфной вектор-функции Φ(z) = (Φ0, Ψ), Φ0 = (Φ1, Φ2), Ψ = (Ψ1, Ψ2, Ψ3) получим
представление
Φ(z) = Φc(a)(z) + Φχ(a)(z) + Φ(z) + ΦF (z), z ∈ Ω,
(39)
где
Φc(a)(z) = Φ(μc(a))(z), Φχ(a)(z) = Φ(μχ(a))(z), ΦF (z) = Φ(μF )(z),
Φ(z) = Φ(Rgc)(z) +Φ(z),
Φ(z) = (c0β0(z), c0β1(z), c0γ0(z) + c1 + ic2, 0, 0),
βj(z) = 2iΘ(tDD1212j)(z), j = 0,1, γ0(z) = Θ(tt)(z),
функция Θ(f)(z) определена в (24), cj - произвольные действительные постоянные.
Теперь выражения ρ(z) из (38) и голоморфных функций из (39) подставим в (16), (18).
Тогда задача (1), (2) сведётся к системе нелинейных уравнений относительно вектор-функции
a = (w1,w2,w312), которую представим в виде
ω0j(z) = ω0jc(a) + ω0(a) + ω0j∗(z) + ω0jF (z), j = 1,2,
w3(z) = w3c(a) + w3χ(a) + w3(z) + w3F (z), z ∈ Ω,
(40)
где ω0jc(a) = ω0jjc(a); ωjc(a)), ωc(a) = (ω1c, ω2c), ωjc(a) = ωjjc(a); ρc(a)), j = 1, 2, w3c(a) =
= w33c(a);ρ3c(a)); остальные слагаемые в системе (40) определяются аналогично; операторы
ωjj;ρj), ω0jj;ωj) и w33;ρ3) определены в (13), (16), (18) соответственно.
Отметим, что функции ω0j∗(z), w3(z) и ω0jF (z), w3F (z) зависят соответственно от про-
извольных постоянных и внешних сил, действующих на оболочку. При этом, как нетрудно
заметить, для функций ω01(z) = w2 + iw1, ω02(z) = ψ2 +1, w3(z) имеют место пред-
ставления (28).
Исследуем разрешимость системы (40) в пространстве
p (Ω).
Лемма 2. Пусть выполнены условия (a), (b), (c), (d). Тогда
1) ω0jc(a), j = 1, 2, w3c(a) - линейные вполне непрерывные операторы в
p (Ω);
2) ω0(a), j = 1, 2, w3χ(a) - нелинейные ограниченные операторы в
p (Ω), причём для
любых aj = (wj1, wj2, wj3, ψj1, ψj2)
p (Ω), j = 1, 2, справедливы оценки
∥ω0(a1) - ω0(a2)
,
∥w3χ(a1) - w3χ(a2)
c(∥a1
+ ∥a2
+
p (Ω)
p (Ω)
p (Ω)
p (Ω)
+ ∥w132
+ ∥w232
)∥a1 - a2
,
(41)
p (Ω)
p (Ω)
p (Ω)
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
670
ТИМЕРГАЛИЕВ
где c - известная положительная постоянная, зависящая от физико-геометрических харак-
теристик оболочки;
3) ω0j∗(z), ω0jF (z)
p (Ω), j = 1, 2.
Доказательство. Из представлений для
fjc (a), fχ(a) в (9), ϕcj(a), ϕχj (a) в (11) сле-
дует, что
fjc (a) и ϕcj (a) - линейные вполне непрерывные, а fχ(a) и ϕχj (a) - нелинейные
ограниченные операторы из
(a), ϕχ3(a)
p (Ω) в Lp(Ω) и в Cα(Γ) соответственно; для
χ
справедливы оценки вида (41), а для fχ(a), ϕχj(a), j = 1, 2, - оценки вида
∥f(a1) - f(a2)L
∥ϕχj (a1) - ϕχj (a2)C
+
p(Ω),
α(Γ)c(∥w3W (2)
p (Ω)
+ ∥w23
)∥a1 - a2
,
j = 1,2.
(42)
p (Ω)
p (Ω)
Тогда из (38) с учётом ограниченности операторов Rj получим, что ρc(a) и μkc(a) - линейные
вполне непрерывные, а ρχ(a) и μ(a) - нелинейные ограниченные операторы из
p (Ω) в
Lp(Ω) и в Cα(Γ) соответственно, и для ρχ(a), μ(a) справедливы оценки (41). Следователь-
но, с учётом свойств интеграла типа Коши из (39) будем иметь, что Φkc(a), Ψ′jc(a) - линейные
вполне непрерывные, Φ(a), Ψ′jχ(a) - нелинейные ограниченные операторы из
p (Ω) в
Cα(Ω), причём для нелинейных операторов Φ(a), Ψ′jχ(a) справедливы оценки (41).
Исследуем свойства операторов
Φ′kc(a) = Θ(μ2kc(a)), k = 1,2, Ψ′′jc(a) = i(j-1)(j-2)/2Θ(μ2j-1c(a)), j = 1,2,3,
(43)
где оператор Θ(f) определён в (32).
Заметим, что функции ρc(a)(z), μkc(a)(t), определённые в (38), являются решениями сис-
темы (26) с правой частью
fjc (a)(z), j = 1, 2, 3, ϕck(a)(t), k = 1, 5. Поэтому μc = (μ1c, . . . , μ5c)
можно представить в виде
[
]
μc(a)(τ)
μc(a)(t) = A-1(t) ϕc(a)(t) - B(t)
dτ - Kμc(a)(t) - Hρc(a)(t) ,
(44)
τ-t
Γ
где A-1(t)∈Cβ(Γ) - матрица, обратная к матрице A(t), ϕc =(ϕc1, . . . , ϕc5), K =(K1, . . . , K5),
H = (H1,...,H5), ρc = (ρ1c2c3c).
Если выражение (44) для μc(a)(t) подставить в (43), после этого переставить порядок
интегрирования в повторных интегралах и использовать указанные выше свойства интеграла
типа Коши, операторов T, S, соотношений (4.7), (4.9) из [14, c. 28-29] и лемму 1, то после
несложных, но достаточно громоздких преобразований получим, что операторы Φ′kc(a), k =
= 1, 2, Ψ′′jc(a), j = 1, 2, 3, суть линейные вполне непрерывные операторы из
p (Ω) в Lp(Ω).
При помощи аналогичных рассуждений также будем иметь, что Φ′kχ(a), k = 1, 2, Ψ′′jχ(a) j =
= 1, 2, 3, - нелинейные ограниченные операторы из
p (Ω) в Lp(Ω) и для них справедливы
оценки (41). Теперь, если использовать соотношения (15), (16), (19) и оценки (42), утверждение
леммы становится очевидным. Лемма доказана.
Систему (40) запишем в виде
a - L(a) - G(a) = a + aF,
(45)
где
L = (L1,...,L5), G = (G1,...,G5), a = (w1,w2,w312),
aF = (w1F ,
w2F ,w3F
ψ1F
ψ2F ), ω01 = w2 + iw1, ω02 = ψ2 +1, L3(a) = w3c(a),
G3(a) = w3χ(a), L3(n-1)+j(a) = -Re [ijω0nc(a)], G3(n-1)+j(a) = -Re[ijω0(a)], n,j = 1,2,
wjF = -Re[ijω01F ],
ψjF = -Re [ijω02F ], j = 1,2,
w3F = w3F .
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
671
Отметим, что L(a) - линейный вполне непрерывный, G(a) - нелинейный ограниченный
операторы в
p (Ω), причём для G(a) имеет место оценка (41); aF
p (Ω) - известная
функция, зависящая от внешних сил; компоненты вектора a даются формулами (28).
Уравнение a - L(a) = 0 имеет лишь нулевое решение в пространстве
p (Ω). Дейст-
вительно, если a ∈
p (Ω) - ненулевое его решение, то, как нетрудно заметить, a является
решением системы линейных уравнений равновесия, удовлетворяющим линейным однородным
граничным условиям. Тогда, рассуждая как и в случае системы (26), приходим к тому, что
вектор a удовлетворяет системе
wj - Gλjjwλ - Bjjw3 = 0, j = 1,2,
w1α2 + w2α1 - 2Gλ12wλ - 2B12w3 = 0, ψj - Gλjjψλ = 0, j = 1,2,
ψ1α2 + ψ2α1 - 2Gλ12ψλ = 0, w3αj + ψj = 0, j = 1,2.
(46)
Выведем условия, при выполнении которых система (46) имеет только нулевое решение.
Для этого введём в рассмотрение вспомогательную систему вида
(ω1/Λ)α1 - (ω2/Λ)α2 = c1f1(α1, α2), (ω1/Λ)α2 + (ω2/Λ)α1 = c1f2(α1, α2),
ω1α1 + ω2α2 = c1f3(α12)
(47)
относительно функций ω1(α1, α2), ω2(α1, α2) ∈ C1α(Ω), где fj
p (Ω), j = 1, 2, 3, - произ-
вольно фиксированные функции, c1 - произвольная действительная постоянная.
Решаем систему (47). При помощи комплексной функции ω = ω1/Λ +2/Λ первые два
уравнения в (47) представим в виде неоднородного уравнения Коши-Римана ωz = c1(f1 +
+ if2)/2, общее решение которого даётся формулой ω = ϕ(z) + c1T (f1 + if2)(z)/2, откуда для
ωj получаем представления
ω1 = Λ(α12)Re ϕ(z) + c1g1(α12), ω2 = Λ(α12)Im ϕ(z) + c1g2(α12),
g1 = ΛRe T(f1 + if2)/2, g2 = ΛIm T(f1 + if2)/2, z = α1 +2,
(48)
где ϕ(z) - произвольная голоморфная в области Ω функция, принадлежащая пространству
C1α(Ω), оператор T определён в (13). Из условия (b) и указанных выше свойств оператора T
следует, что функции gj , j = 1, 2, принадлежат пространству
p (Ω).
Выражения решений ω1, ω2 из (48) подставим в третье уравнение системы (47). Тогда
относительно голоморфной функции ϕ(z) получим уравнение вида
Re[ϕ(z) + (ln Λ)zϕ(z)] = c1g3(α12), g3(α12) = (f3 - g1α1 - g2α2 )/(2Λ),
(49)
где функция g3 принадлежит пространству
p (Ω).
Предположим, что в области Ω выполняется условие
Λ0(α12) = (ln Λ)zz = 0,
(50)
т.е. функция ln Λ(α1, α2) в области Ω не является гармонической.
Заметим, что условие (50) в силу уравнения Гаусса [18, c. 193] означает, что гауссова кри-
визна срединной поверхности оболочки не равна нулю.
Найдём решение уравнения (49). Заметим, что функция Re ϕ(z), как действительная
часть голоморфной функции ϕ(z), в области Ω является гармонической, т.е. [Re ϕ(z)]zz =
= 0. Следовательно, в области Ω при любой постоянной c1 должно выполняться условие
[c1g3 - Re ((ln Λ)zϕ)]zz = 0, которое можно записать в виде Re (Λ0ϕ)z - c1g3zz = 0, c1 = const,
z ∈ Ω, откуда с учётом условий (50) и g3zz = 0, z ∈ Ω, получим Re(Λ0ϕ)z = 0, c1 = 0. Итак,
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
672
ТИМЕРГАЛИЕВ
голоморфная функция ϕ(z) удовлетворяет системе уравнений Re (Λϕ)z = 0, Re (Λ0ϕ)z = 0,
которая эквивалентна системе
uα1 + Λα1 u + Λα2 v = 0,
0uα1 + Λ0α1 u + Λ0α2 v = 0,
uα1 = vα2 , uα2 = -vα1, Λ0/Λ = const,
(51)
где ϕ(z) = u(α1, α2) + iv(α1, α2), функция Λ0(α1, α2) определена в (50).
Исключив из первых двух уравнений в (51) производную uα1 , получим
v = Λ1u, Λ1Λ1(α12) = -0/Λ)α1/0/Λ)α2 ∈ W(1)p(Ω).
(52)
Теперь, если выражение функции v из (52) подставить в первое уравнение (51), относи-
тельно функции u получим уравнение
uα1 + Λ2u = 0, Λ2 Λ2(α12) = [(ln Λ)α1 + (ln Λ)α2 Λ1]/2 ∈ W(1)p(Ω),
общее решение которого даётся формулой
[ α1
]
u = c(α23, Λ3Λ3(α12) = exp
-
Λ2(β1, α2)1
∈ W(1)p(Ω),
(53)
α1
0
где c(α2) - произвольная действительная функция переменной α2, принадлежащая прост-
ранству C1α, (α10, α20) Ω - произвольно фиксированная точка.
Тогда для функции v с учётом (52) будем иметь представление
v = Λ1Λ3c(α2).
(54)
Теперь функции u(α1, α2), v(α1, α2) из (53), (54) подставим в третье и четвертое уравнения
системы (51). В результате относительно функции c(α2) почти всюду в области Ω получим
уравнения
Λ3c(α2) + [Λ3α2 + (Λ1Λ3)α1 ]c(α2) = 0, Λ1Λ3c(α2) -3α1 -1Λ3)α2 ]c(α2) = 0,
из которых видно, что c(α2) 0, следовательно, из формул (53), (54) будем иметь ϕ(z) 0
в Ω. Тогда из соотношений (48) получим ω1 = ω2 0 в Ω.
Перейдём к нахождению функций ψ1, ψ2. С этой целью четвертое и пятое равенства
системы (46) сложим и вычтем друг из друга. Принимая во внимание соотношения (4) для
символов Кристоффеля, будем иметь систему
(ψ1/Λ)α1 - (ψ2/Λ)α2 = 0, (ψ1/Λ)α2 + (ψ2/Λ)α1 = 0, ψ1α1 + ψ2α2 = 0,
которая получается из системы (47) при fj 0, j = 1, 2, 3, ω1 = ψ1, ω2 = ψ2. Поэтому в силу
доказанного выше ψ1(z) = ψ2(z) 0 в Ω и из последних двух равенств в (46) следует, что
w3 = c1 = const в Ω. Подставляя это значение функции w3 в первые три равенства системы
(46), при помощи аналогичных рассуждений снова приходим к системе вида (47), в которой
ω1 = w1, ω2 = w2, f1 = (B11 - B22)/Λ, f2 = 2B12/Λ, f3 = B11 + B22. В силу условия (а)
имеем fj
p (Ω), j = 1, 2, 3, при этом условие g3zz = 0 примет вид
[(B11 + B22)/Λ - Re (ln Λ)zRe T f + Im (ln Λ)zIm T f - Re Sf]zz = 0,
f = (B11 - B22 + 2iB12)/Λ,
(55)
где операторы T, S определены в (13), (14). Следовательно, wj 0, j = 1, 2, 3, в Ω. Итак,
уравнение a - L(a) = 0 имеет только нулевое решение в
p (Ω). Таким образом, существует
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
О СУЩЕСТВОВАНИИ РЕШЕНИЙ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ
673
обратный оператор (I - L)-1, ограниченный в
p (Ω), с помощью которого уравнение (45)
сведётся к эквивалентному уравнению
a - G(a) = aF,
(56)
где G(a) = (I - L)-1G(a), aF = (I - L)-1aF .
Отметим, что вектор ac = (I - L)-1a является решением однородной системы линей-
ных уравнений равновесия, удовлетворяющим однородным линейным граничным условиям.
Поэтому в силу доказанного выше ac 0, что и учтено нами при переходе к уравнению (56).
Также отметим, что вектор aF в (56) зависит только от внешних сил и aF = 0, если
внешние силы отсутствуют.
Лемма 3. Пусть выполнены условия (a), (b), (c), (d). Тогда:
1) G(a) - нелинейный ограниченный оператор в
p (Ω), причём для любых векторов
aj = (wj1,wj2,wj3j1j2), j = 1,2, справедлива оценка
∥G(a1)-G(a2)
c(∥a1
+∥a2
+∥w132
+∥w232
)∥a1-a2
,
p (Ω)
p (Ω)
p (Ω)
p (Ω)
p (Ω)
p (Ω)
где c - известная положительная постоянная, зависящая от физико-геометрических ха-
рактеристик оболочки;
2) aF
p (Ω).
Справедливость этой леммы вытекает из леммы 2 с учётом указанных выше свойств опе-
раторов (I - L)-1 и G.
Исследуем разрешимость уравнения (56) в пространстве
p (Ω). Используя лемму 3, для
любых aj
p (Ω), j = 1, 2, принадлежащих шару ∥a∥W (2)
< r, получаем
p (Ω)
∥G(a1) - G(a2)
q∥a1 - a2
,
q = 2cr(1 + r).
p (Ω)
p (Ω)
Предположим, что радиус r шара и внешние силы таковы, что выполняются неравенства
q < 1,
∥aF
< (1 - q)r.
(57)
p (Ω)
Тогда к уравнению (56) можно применить принцип сжатых отображений [19, c. 146], согласно
которому уравнение (56) в шаре ∥a∥
< r имеет единственное решение вида a = R(aF)
p (Ω)
p (Ω), где R - резольвента оператора G.
Заметим, что если внешняя нагрузка отсутствует, то задача (1), (2) имеет только нулевое
решение.
Вернёмся к условиям разрешимости (37), в которых под a = (w1, w2, w3, ψ1, ψ2)
p (Ω)
будем подразумевать решение задачи (1), (2). Используя равенства (1) и (2), убеждаемся в
том, что условия разрешимости (37) выполняются.
Таким образом, доказана следующая
Теорема. Пусть выполнены условия (a), (b), (c), (d), (36), (50), (55) и неравенства (57).
Тогда задача (1), (2) имеет единственное обобщённое решение a = (w1, w2, w3, ψ1, ψ2)
p (Ω), 2 < p < 4/(2 - β).
СПИСОК ЛИТЕРАТУРЫ
1. Галимов К.З. Основы нелинейной теории тонких оболочек. Казань, 1975.
2. Ворович И.И. Математические проблемы нелинейной теории пологих оболочек. М., 1989.
3. Морозов Н.Ф. Избранные двумерные задачи теории упругости. Л., 1978.
4. Карчевский М.М. Исследование разрешимости нелинейной задачи о равновесии пологой незакреп-
лённой оболочки // Уч. зап. Казанского. ун-та. Сер. Физ.-мат. науки. 2013. Т. 155. № 3. С. 105-110.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023
674
ТИМЕРГАЛИЕВ
5. Тимергалиев С.Н. Теоремы существования в нелинейной теории тонких упругих оболочек. Казань,
2011.
6. Тимергалиев С.Н. О существовании решений геометрически нелинейных задач для пологих оболо-
чек типа Тимошенко со свободными краями // Изв. вузов. Математика. 2014. № 3. С. 40-56.
7. Тимергалиев С.Н. К вопросу о существовании решений нелинейной краевой задачи для системы
дифференциальных уравнений с частными производными теории пологих оболочек типа Тимошен-
ко со свободными краями // Дифференц. уравнения. 2015. Т. 51. № 3. С. 373-386.
8. Тимергалиев С.Н., Харасова Л.С. Исследование разрешимости одной краевой задачи для системы
нелинейных дифференциальных уравнений теории пологих оболочек типа Тимошенко // Диффе-
ренц. уравнения. 2016. Т. 52. № 5. С. 651-664.
9. Тимергалиев С.Н. Метод интегральных уравнений в нелинейных краевых задачах для пологих
оболочек типа Тимошенко со свободными краями // Изв. вузов. Математика. 2017. № 4. С. 59-75.
10. Тимергалиев С.Н. К проблеме разрешимости нелинейных задач равновесия пологих оболочек типа
Тимошенко // Прикл. математика и механика. 2018. Т. 82. № 1. С. 98-113.
11. Тимергалиев С.Н. Метод интегральных уравнений исследования разрешимости краевых задач для
системы нелинейных дифференциальных уравнений теории пологих неоднородных оболочек типа
Тимошенко // Дифференц. уравнения. 2019. Т. 55. № 2. С. 239-255.
12. Тимергалиев С.Н. К проблеме разрешимости нелинейных краевых задач для произвольных изо-
тропных пологих оболочек типа Тимошенко со свободными краями // Изв. вузов. Математика.
2021. № 4. С. 90-107.
13. Тимергалиев С.Н. О разрешимости нелинейных краевых задач для системы дифференциальных
уравнений равновесия пологих анизотропных оболочек типа Тимошенко с незакреплёнными краями
// Дифференц. уравнения. 2021. Т. 57. № 4. С. 507-525.
14. Векуа И.Н. Обобщённые аналитические функции. М., 1988.
15. Мусхелишвили Н.И. Сингулярные интегральные уравнения. М., 1962.
16. Пресдорф З. Некоторые классы сингулярных уравнений. М., 1979.
17. Гахов Ф.Д. Краевые задачи. М., 1963.
18. Векуа И.Н. Основы тензорного анализа и теории ковариантов. М., 1978.
19. Красносельский М.А. Топологические методы в теории нелинейных интегральных уравнений. М.,
1956.
Казанский государственный
Поступила в редакцию 17.12.2022 г.
архитектурно-строительный университет
После доработки 13.02.2023 г.
Принята к публикации 16.03.2023 г.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 59
№5
2023