CC BY
54
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2017. Т. 21, № 3. С. 496-506 188№ 2310-7081 (опИпе), 1991-8615 (ргт!) с! http://doi.org/10
УДК 539.3
Интегро-дифференциальные уравнения второй краевой задачи линейной теории упругости. Сообщение 1. Однородное изотропное тело
В. В. Стружанов
Институт машиноведения УрО РАН,
Россия, 620049, Екатеринбург, ул. Комсомольская, 34.
Аннотация
Система уравнений второй краевой задачи линейной теории упругости для однородного изотропного тела сведена к двум разным интегро-дифференциальным уравнениям фредгольмовского типа, что позволило для их исследования применить теоремы Фредгольма. Определены спектральные радиусы соответствующих операторов и доказано существование и единственность решения второй краевой задачи. Также установлено, что решение второго интегро-дифференциального уравнения можно найти методом последовательных приближений и представить его сходящимся со скоростью геометрической прогрессии рядом Неймана. Применение методики проиллюстрировано на примере расчета остаточных напряжений в закаленном цилиндре.
Ключевые слова: вторая краевая задача, однородное изотропное тело, интегро-дифференциальное уравнение, спектральный радиус, последовательные приближения.
Получение: 12 июля 2017 г. / Исправление: 23 августа 2017 г. / Принятие: 18 сентября 2017 г. / Публикация онлайн: 22 сентября 2017 г.
Введение. Классическая теория упругости сохраняет свое почетное место в науке о поведении деформируемого твердого тела. Ее исходные определения являются общими для всех разделов этой науки, а методы постановки и решения задач служат для нее образцами.
Несмотря на то, что линейная теория упругости является полностью замкнутой математической теорией, положение в этой области постоянно меняется. Это связано в первую очередь с использованием все более мощного формального аппарата. В этом случае система уравнений изучается сама по себе на принятом в математике уровне строгости, т.е. исследуются чисто
Краткое сообщение
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Стружанов В.В. Интегро-дифференциальные уравнения второй краевой задачи линейной теории упругости. Сообщение 1. Однородное изотропное тело // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2017. Т. 21, № 3. С. 496-506. doi: 10.14498/vsgtu1555. Сведения об авторе
Валерий Владимирович Стружанов А http://orcid.org/0000-0002-3669-2032 доктор физико-математических наук, профессор; главный научный сотрудник; лаб. микромеханики материалов; e-mail: stru@imach.uran.ru
.14498/vsgtu1555
математические свойства, что позволяет раскрыть все еще неизвестные возможности теории и получить новые эффективные методы решения краевых задач, в том числе и приближенных, имеющих большое значение для решения важных технических задач.
В данной работе показано, что систему уравнений второй краевой задачи линейной теории упругости однородного изотропного тела можно свести к одному из двух интегро-дифференциальных уравнений, которые относятся к классу уравнений Фредгольма второго рода. Это позволило применить теоремы Фредгольма, т. е. свести задачу к проблеме собственных чисел соответствующих интегро-дифференциальных операторов. Для каждого оператора определены их спектральные радиусы, использование которых позволило привести новое доказательство существования и единственности решения второй краевой задачи. При исследовании второго интегро-дифференциального уравнения было установлено, что его решение (следовательно, и решение второй краевой задачи) можно найти методом последовательных приближений и представить его сходящимся рядом Неймана. Определены условия, при выполнении которых ряд Неймана возможно свернуть. В качестве примера решена задача об определении остаточных напряжений в закаленном цилиндре.
1. Первое интегро-дифференциальное уравнение второй краевой задачи теории упругости для однородного изотропного тела.
Система уравнений краевой задачи линейной теории упругости, записанная в инвариантной форме, имеет вид [1,2]:
V- а + д = 0, е = ёе! V, а = С ■-(е — е*), «|г = (1)
Требуется найти тензоры напряжений а и деформаций е, вектор перемещений V внутри упругого тела V в состоянии равновесия, если известны перемещения vг точек его границы (кусочно-гладкой поверхности Г).
В системе (1) первая группа уравнений — это уравнения равновесия (д — вектор объемных сил, V — набла-оператор Гамильтона [1], точкой обозначено скалярное произведение вектора Гамильтона на симметричный тензор второго ранга), вторая группа — соотношения Коши, третья — закон Гука (С — однородный изотропный тензор четвертого ранга модулей упругости, двумя точками обозначено двойное скалярное произведение тензоров [3], е* — тензор первоначальных деформаций свободных от связей элементов тела V, возникающих при нагреве, фазовых превращениях и т. п. [4]).
В уравнениях (1) произведем замену V = и + и*, где и — неизвестная вектор-функция, на границе равная нулю, а и* — известная вектор-функция, на границе равная и*|г = vг, например, гармоническая вектор-функция
* ( г№ и = — V —— аГ, 7г ап
где С — функция Грина оператора Лапласа А для области V, п — внешняя нормаль к поверхности Г [5].
Подставляя теперь закон Гука в уравнения равновесия и заменяя деформации соотношениями Коши, получаем уравнения Навье—Ляме [6]:
— (^Аи + (Л + и) = р, и|г = 0, (2)
где Л = ((1 + — )) \ ^ = Е(2(1 + V)) 1 — постоянные Ляме; Е, V — соответственно модуль Юнга и коэффициент Пуассона;
р = д + V ■ а* + ц,Аи* + (Л + grad ё1у и*;
а* = С ■ •£* —формальный тензор псевдонапряжений.
Уравнения (2) будем рассматривать как некоторое отображение пространства ) в пространство Ь2^) (и е ), р е Ь2^)). Здесь ) — вещественное полное сепарабельное гильбертово пространство вектор-функций, компоненты которых интегрируемы с квадратом, а ) — вещественное полное сепарабельное гильбертово пространство вектор-функций, компоненты которых обращаются в нуль на Г и принадлежат ¿2 (V) вместе со своими обобщенными производными до второго порядка включительно [7,8]. Пусть область V и граница Г таковы, что оператор Лапласа А устанавливает биекцию между пространствами Ж^^) и Ь2^) и существует функция Грина оператора Лапласа. Таким образом, оператор (—А) имеет обратный [7,9]:
(—А-1 X) = / СХ^, X е ^2^), А-1Х е
■IV ,
Оператор (—А) — положительно определенный оператор в ) с дискрет-
ным спектром, то есть его собственные числа вещественны, положительны и дискретны [10]. Следовательно, оператор (—А-1) является вполне непрерывным.
Применяя к обеим частям уравнения (2) оператор (—А-1), получаем инте-гро-дифференциальное уравнение (первое интегро-дифференциальное уравнение второй краевой задачи)
и = Пи + f. (3)
Здесь
Пи = —тА 1 graddiv и, f = — А 1р, т =
^ ^ 1 — 2v
Оператор П действует из ) в ^220(V), / е ^220(V). Дифференциаль-
ный оператор grad diу — линейный непрерывный (ограниченный) оператор из пространства в пространство Ь2(V) [11,12]. Тогда оператор П вполне
непрерывен как произведение вполне непрерывного и ограниченного операторов [13]. Таким образом, уравнение (3)—уравнение Фредгольма второго рода [14].
2. Спектральный радиус оператора И и существование решения второй краевой задачи. Для уравнения (3) справедлива альтернатива Фредгольма [14]. Поэтому вопрос существования и единственности решения сводится к проблеме собственных чисел оператора П. Для определения его спектрального радиуса р(П) рассмотрим уравнение
—(Аи + тgraddivи) — к(—Аи) = 0, и е к е (—те, +те). (4)
Выражение в первых скобках — оператор теории упругости. В случае закрепленной границы он положительно определенный с дискретным спектром [10]. Таков же и оператор (—Д) [10]. Тогда собственные числа k уравнения (4) вещественны, положительны и лежат в пределах inf A(u) = ki ^ k ^ k2 = = sup A(u) [10], где
, . (-[Ди + m graddiv u], u) (— graddiv u, u)
A( u) =--—---- = I + m--:—г---. (5)
v 7 (-Дщ, u) (-Дu, u) v '
Здесь круглыми скобками обозначены скалярные произведения в ). Далее, применяя формулу Остроградского—Гаусса и учитывая равенство [15]
Дu = grad div u — rot rot u, (6)
находим
(— graddivu,u) = (divu,divu) ^ 0, (—Д^u) = (gradu,gradu) ^ 0,
(—Д^u) — (— graddivu,u) = (rotrotu,u) = (rotu,rotu) ^ 0. Отсюда
(—Au,u) ^ (— graddivu,u) ^ 0.
Тогда, используя выражение (5), находим, что собственные числа k уравнения (4) лежат в пределах I ^ k ^ I + m.
Применяя к равенству (4) оператор (—Д-1), получаем эквивалентное уравнение su = Hu (s = I — k). Собственные числа s оператора П лежат в отрезке [—m, 0]. Для вполне непрерывного линейного оператора спектральный радиус равен наибольшему по модулю собственному значению оператора [16]. Тогда р(П) = m. Если 0 ^ v ^ 0.5, то m е (I, те) и р(П) > I. Отсюда оператор П не является оператором сжатия [16] и решение уравнения (3) невозможно определить методом последовательных приближений.
Отметим, что k = 0 не является собственным числом уравнения (4). Отсюда s = I также не является собственном числом оператора П. Согласно теореме Фредгольма [14] решение уравнения (3) существует и единственно. Таким образом, решение второй краевой задачи теории упругости также существует и единственно при 0 < v < 0.5. Если v = 0.5, то m = те и вопрос о существовании и единственности решения требует специального исследования.
3. Второе интегро-дифференциальное уравнение второй краевой задачи. Используя равенство (6), из уравнения (3) получаем второе интегро-дифференциальное уравнение
u = Qu + h, (7)
где
1 I m 1
— 1 ^ Я 7
Qu = —1Д rotrot u, h =-f, l =
1 + т 1 + т 2(1 — V)'
Здесь оператор Q также действует из ) в Ж^о^) и является вполне
непрерывным, т. е. уравнение (7) — уравнение Фредгольма второго рода.
Для отыскания спектрального радиуса оператора Q применим прием, изложенный выше. Рассмотрим уравнение
-^u + l rot rot u) — y(—Д^ = 0, u е W220(V), 7 е (—те, +те). (8)
Собственные числа 7 уравнения (8) вещественны, положительны и лежат в промежутке inf B (u) = 71 ^ 7 ^ 72 = sup B (u), где
(—^u + l rot rot u], u) (— rot rot u, u)
B(u) =--—----= I + l—-—---—. (9)
v ' (-Дu, u) (—Д^ u) v 7
Далее имеем оценки
(rot rot u, u) = (rot u, rot u) ^ 0, (—Д^ u) + (— rot rot u, u) = (— grad div u, u) ^ 0, 0 ^ (— rot rot u, u) ^ — (—Д^ u).
Теперь, используя выражение (9), находим, что собственные числа 7 уравнения (8) лежат в промежутке (I — l) ^ 7 ^ I.
Применяя оператор (—Д-1) к (8), получаем эквивалентное уравнение Au = = Qu (A = I — 7). Собственные числа оператора Q лежат в отрезке [0, l]. Спектральный радиус p(Q) = l. При 0 ^ v < 0.5 имеем p(Q) < I, то есть оператор Q является оператором сжатия. Тогда решение уравнения (7) можно найти методом последовательных приближений и представить его сходящимся по норме пространства W|0(V) рядом Неймана u = Qnh. Причем ряд сходится со скоростью геометрической прогрессии со знаменателем, сколь угодно близким к l [16].
Отметим, что 7 = 0 (v < 0.5) не является собственным числом уравнения (8) и A = I не является собственным числом оператора Q. Отсюда из теоремы Фредгольма вытекает существование и единственность решения уравнения (7).
Если v = 0.5 (материал несжимаемый), то l = I и A = I —собственное число оператора Q. В данном случае вопрос о существовании и единственности решения уравнения (7), а также второй краевой задачи теории упругости, требует специального исследования [14].
4. Потенциальные и соленоидальные поля. Согласно теореме Гельм-гольца, всякое непрерывное векторное поле можно представить в виде суммы потенциального и соленоидального полей [9]:
w = grad ^ + rot A.
Здесь функция ^ — скалярный потенциал, а вектор A — векторный потенциал. В нашем случае w = grad ^>+rot A е W22 0(V). Рассмотрим соленоидальное поле w = rot A. Тогда из равенства (6) вытекает
Дгс^ A = — rot rot rot A. (10)
Таким образом, в соленоидальном поле Д = — rot rot. Возьмем теперь вектор Д-1 rot A. Используя равенство (6), получаем
grad div Д-1 rot A — rot rot Д-1 rot A = rot A,
т. е.
grad div A-1 rot A = 0
и вектор A-1 rot A принадлежит соленоидальному полю. Если w = grad то A grad = grad div grad Отсюда A = grad div. Возьмем вектор A-1 grad Из равенства (6) следует, что
grad div A-1 grad ^ — rot rot A-1 grad ^ = grad
т. е. A-1 grad ^ — вектор потенциального поля.
Представим теперь вектор h в уравнении (7) суммой h = grad ф + rot B. Тогда, используя равенство (10), можно свернуть ряд Неймана. В результате получаем решение
u = grad ф + --- rot B.
1 — l
5. Остаточные напряжения в закаленном цилиндре. В качестве примера определим напряжения в закаленном длинном стальном круговом цилиндре. В результате закалки часть зерен аустенита в приповерхностных слоях перешла в мартенситное состояние. Пусть P — объемное содержание мартенсита:
Г 0, 0 ^ r < а; P = < ^ r — а
P | Pot-, а < r < b,
I b — а
где b — радиус основания цилиндра, (b — а) — глубина закалки, Po ^ 1 — объемное содержание мартенсита в поверхностном слое. Отсюда после закалки цилиндр состоит из двух однородных изотропных компонентов, каковыми являются аустенит и мартенсит. В первом приближении полагаем их свойства одинаковыми. Так как зерна мартенсита имеют несколько больший объем [17], элементарные объемы материала в приповерхностных слоях находятся в стесненном состоянии, что вызывает появление остаточных (закалочных) напряжений.
Рассмотрим произвольный свободный от связей кубический элемент материала. После закалки в нем образуются равномерно распределенные по объему зерна мартенсита с объемным содержанием P. В результате объем элемента увеличится с сохранением кубической формы. В этом случае компоненты деформации
s j = aP5ij, i, j = 1,2,3.
Здесь 5ij — символ Кронекера, а — параметр свободной структурной деформации мартенсита [17].
Деформации sjjj не удовлетворяют условиям совместности и не могут быть реализованы в сплошном теле. Для сохранения сплошности к деформациям sj необходимо добавить такие деформации sj, чтобы суммарные деформации si j = sj + sj уже удовлетворяли условиям совместности. Это означает, что к каждому элементу должны быть приложены усилия, определяемые напряжениями а- (собственные, остаточные напряжения). Напряжения aj и деформации sj связаны законом Гука
aij = Cij«esaв = Cij«e(s«e — s«e), i, j, а, в = 1, 2, 3.
Здесь Cijae — компоненты тензора модулей упругости (суммирование по повторяющимся индексам).
Стесненный компонент расположен симметрично относительно оси цилиндра. Поэтому точки поверхности получают постоянные по величине радиальные перемещения:
vr I , = v° = const.
' lr=o r
Таким образом, задача по определению закалочных напряжений осесиммет-ричная и цилиндр находится в плоском деформированном состоянии. Так как
^11 = 4(r) ^22 = £l(r) £3з = е;12 = £1з = 4з = 0,
£ И = ^22 = £3з = aP (r) £ 1 2 = £ 1 3 = £2з = 0
уравнения закона Гука имеют вид
оГ' = + A(4 + 4) - a(3A + 2p)P,
of = 2^4 + A(£r + 4) - a(3A + 2p)P, Oz = A(£r + 4) - a(3A + 2p)P. Подставляя сюда соотношения Коши
j dvr f Vr
£r = ~T~, £l =
dr r
и затем полученное выражение в уравнение равновесия
dor + or ol = 0
dr r
после известных преобразований получаем уравнение
. , 3A + 2u . .
Дv + l rot rot v = a——-— grad P, v = vr (r). A + 2p
Затем, делая замену
v0r r
v = u + u , u = , u = ur (r),
имеем
3A + 2u .
Дu + l rot rot u = a—-grad P, ur I , = 0.
A + 2p ' r|r=b
Данное уравнение эквивалентно интегро-дифференциальному уравнению
. . _1 3A + 2u 1
u = —l Д rot rot u + —-аД grad P.
A + 2p 6
Здесь оператор (—Д-1) определен с использованием функции Грина для круга [18].
Как следует из рассуждений, приведенных выше, решение интегро-диф-ференциального уравнения определяет ряд Неймана. Так как вектор
А 1 grad Р принадлежит потенциальному полю, ряд Неймана имеет только один первый член:
ЗА + 2и _1
ur = a-
Л + 2^
Тогда
Здесь
■A-i grad P.
3Л + 2u
vr = u* + ur = vOt- + a——:—— A grad P. b Л + 2^
P
о
A-1 grad P =
,3„
br ar a3r\
b — a V 3
P0 (3(b — r) —
3 + 6b2^'
a3
— r2)),
0 ^ r ^ a, a ^ r ^ b.
Ь — а V Зч" 7 6Ь2г
Подставляя уг в соотношения Коши и полученные деформации в закон Гука, вспоминая, что цилиндр не нагружен, то есть = 0, вычисляем значение
перемещения точек границы V0. Переходя снова к напряжениям, получаем выражения для закалочных напряжений:
aPoE
it
or =
^ =
(1 — v )(b — - a)b2 V
aPoE (1
(1 — v )(b — -a) V3
aPoE (
(1 — v)(b — a)b2
aPoE ( ^
(b3 _ ^ 2
a
+ ar
a3 ( 1
6 b2
b2a a3 4 + "6"
'"> + £ ( b2
1
0 ^ r ^ a, a ^ r ^ b;
0 ^ r ^ a, , a ^ r ^ b.
(1 — V )(Ь — а) \3
Напряжение вдоль оси цилиндра определим по формуле, которая является следствием равенства нулю деформации :
О" =
v« + of), 0 < r < a, v(о? + о#) - Ее*, a < r < b.
(11)
Напряжение (11) возникает в том случае, когда торцевые поверхности закреплены от осевого перемещения. Если они свободны, то на напряжение (11) следует наложить равномерное напряжение Ee" (е" = const). При этом деформация е" подбирается так, чтобы равнодействующая напряжений, распределенных по торцевой поверхности, обращалась в нуль:
, 2aP0 (b a a3 ) е" = b—a v 3 — 2 + 6b2 / .
В результате получаем выражение для о"', при котором граничные условия на торцах удовлетворяются в смысле принципа Сен—Венана:
о"' = —Ее" + О" =
2aP0E (b a a3 )
(b — a)(1 — v) ( 3 — 2 + 6b2 ): 2aP0E (b r a3 )
(b — a)(1 — v) U — 2+6b2J,
0 ^ r ^ a, a ^ r ^ b.
2
1
2
r
3
Заключение. Приведенные результаты показывают, что применение нетрадиционного для теории упругости формального аппарата функционального анализа позволяет не только исследовать чисто математические свойства второй краевой задачи, но и разработать оригинальный метод последовательных приближений для нахождения ее решения.
Конкурирующие интересы. У меня нет конкурирующих интересов. Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Исследование не имело финансирования.
Библиографический список
1. Лурье А. И. Теория упругости. М.: Наука, 1970. 939 с.
2. Елисеев В. В. Механика упругих тел. СПб.: СПбГПУ, 2002. 341 с.
3. Димитриенко Ю. И. Тензорное исчисление. М.: Высш. шк., 2001. 575 с.
4. Timoshenko S. P., Goodier J. N. Theory of elasticity/ Engineering Societies Monographs. International Student Edition. New York: McGraw-Hill Book Comp., 1970. xxiv+567 pp.
5. Владимиров В. С. Уравнения математической физики. М.: Наука, 1971. 512 с.
6. Hahn H. G. Elastizitätstheorie. Grundlagen der linearen Theorie und Anwendungen auf eindimensionale, ebene und räumliche Probleme / Leitfäden der Angewandten Mathematik und Mechanik. vol. 62. Stuttgart: B. G. Teubner, 1985. 332 pp.
7. Ладыженская О. А., Уральцева Н. Н. Линейные и квазилинейные уравнения эллиптического типа. М.: Наука, 1973. 576 с.
8. Соболев С. Л. Некоторые применения функционального анализа в математической физике. М.: Наука, 1988. 334 с.
9. Кошляков Н. С., Глинер Э. Б., Смирнов М. М. Уравнения в частных производных математической физики. М.: Высш. шк., 1970. 712 с.
10. Михлин С. Г. Вариационные методы в математической физике. М.: Наука, 1970. 512 с.
11. Треногин В. А. Функциональный анализ. М.: Наука, 1980. 496 с.
12. Функциональный анализ / Справочная математическая библиотека / ред. С. Г. Крейн. М.: Наука, 1972. 544 с.
13. Люстерник Л. А., Соболев В. И. Элементы функционального анализа. М.: Наука, 1965. 520 с.
14. Кантарович Л. В., Акилов Г. П. Функциональный анализ. М.: Наука, 1977. 741 с.
15. Коренев Г. В. Тензорное исчисление. М.: МФТИ, 2000. 240 с.
16. Красносельский М. А., Вайникко Г. М., Забрейко П. П. Приближенное решение операторных уравнений. М.: Высш. шк., 1969. 455 с.
17. Юрьев С. Ф. Удельные объемы фаз в мартенситном превращении аустенита. М.: Металлургиздат, 1950. 48 с.
18. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М.: Наука, 1966. 724 с.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 3, pp. 496-506
ISSN: 2310-7081 (online), 1991-8615 (print) d http://doi.org/10.14498/vsgtu1555
MSC: 74C10
Integro-differential equations the second boundary value problem of linear elasticity theory. Message 1. Homogeneous isotropic body
V. V. Struzhanov
Institute of Engineering Science, Ural Branch of RAS,
34, Komsomolskaya st., Ekaterinburg, 620049, Russian Federation.
Abstract
The system of equations of the second boundary value problem of the linear theory of elasticity for homogeneous isotropic bodies is reduced to two separate integro-differential equations of Fredholm type, which allowed to apply for their research the theorem of Fredholm. The spectral radii of the corresponding operators are determined and the existence and uniqueness of the solution of the second boundary value problem are proved. It is also established that the decision of the second integro-differential equation can be found by successive approximations and presented convergent with a geometric rate close to Neumann. The method application is illustrated on the example of calculation of residual stresses in a quenched cylinder.
Keywords: second boundary-value problem, homogeneous isotropic body, integro-differential equation, spectral radius, successive approximation.
Received: 12th July, 2017 / Revised: 23rd August, 2017 / Accepted: 18th September, 2017 / First online: 22nd September, 2017
Competing interests. I have no competing interests.
Author's Responsibilities. I take full responsibility for submitting the final manuscript in print. I approved the final version of the manuscript. Funding. The research has not had any sponsorship.
Short Communication
3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Struzhanov V. V. Integro-differential equations the second boundary value problem of linear elasticity theory. Message 1. Homogeneous isotropic body, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 3, pp. 496-506. doi: 10.14498/vsgtu1555 (In Russian). Author's Details:
Valery V. Struzhanov A http://orcid.org/0000-0002-3669-2032
Dr. Phys. & Math. Sci., Professor; Chief Researcher; Lab. of Material Micromechanics;
e-mail: stru@imach.uran.ru
References
1. Lurie A. I. Theory of Elasticity, Foundations of Engineering Mechanics. Berlin, Heidelberg, Springer, 2005, 1050 pp. doi: 10.1007/978-3-540-26455-2.
2. Eliseev V. V. Mekhanika uprugikh tel [Mechanics of Elastic Bodies]. Saint-Petersburg, SPbGPU, 2002, 341 pp. (In Russian)
3. Dimitrienko Yu. I. Tenzornoe ischislenie [Tensor Calculus]. Moscow, Vyssh. shk., 2001, 575 pp. (In Russian)
4. Timoshenko S. P., Goodier J. N. Theory of elasticity, Engineering Societies Monographs. International Student Edition. New York, McGraw-Hill Book Comp., 1970, xxiv+567 pp.
5. Vladimirov V. S. Uravneniia matematicheskoi fiziki [Equations of mathematical physics]. Moscow, Nauka, 1971, 512 pp. (In Russian)
6. Hahn H. G. Elastizitätstheorie. Grundlagen der linearen Theorie und Anwendungen auf eindimensionale, ebene und räumliche Probleme, Leitfäden der Angewandten Mathematik und Mechanik, vol. 62. Stuttgart, B. G. Teubner, 1985, 332 pp.
7. Ladyzhenskaia O. A., Ural'tseva N. N. Lineinye i kvazilineinye uravneniia ellipticheskogo tipa [Linear and quasilinear equations of elliptic type]. Moscow, Nauka, 1973, 576 pp. (In Russian)
8. Sobolev S. L. Nekotorye primeneniia funktsional'nogo analiza v matematicheskoi fizike [Some applications of functional analysis in mathematical physics]. Moscow, Nauka, 1988, 334 pp. (In Russian)
9. Koshlyakov N. S., Gliner E. B., Smirnov M. M. Uravneniia v chastnykh proizvodnykh matem-aticheskoi fiziki [Equations in partial derivatives of mathematical physics]. Moscow, Vyssh. shk., 1970, 712 pp. (In Russian)
10. Mikhlin S. G. Variatsionnye metody v matematicheskoi fizike [Variational Methods in Mathematical Physics]. Moscow, Nauka, 1970, 512 pp. (In Russian)
11. Trenogin V. A. Funktsional'nyi analiz [Functional analysis]. Moscow, Nauka, 1980, 496 pp. (In Russian)
12. Functional analysis, Wolters-Noordhoff Series of Monographs and Textbooks on Pure and Applied Mathematics, ed. S. G. Krejn. Groningen, Netherlands, Wolters-Noordhoff Publ., 1972, xv+379 pp.
13. Lyusternik L. A., Sobolev V. I. Elementy funktsional'nogo analiza [Elements of Functional Analysis]. Moscow, Nauka, 1965, 520 pp. (In Russian)
14. Kantorovich L. V., Akilov G. P. Funktsional'nyi analiz [Functional Analysis]. Moscow, Nauka, 1977, 741 pp. (In Russian)
15. Korenev G. V. Tenzornoe ischislenie [Tensor Calculus]. Moscow, Mosk. Fiz.-Tekhn. Inst., 2000, 240 pp. (In Russian)
16. Krasnosel'sky M. A., Vainikko G. M., Zabreiko P. P., Rutitsky Ya. B., Stetsenko V. Ya. Pri-blizhennoe reshenie operatornykh uravnenii [Approximate Solution of Operator Equations]. Moscow, Nauka, 1969, 455 pp. (In Russian)
17. Yuriev S. F. Udel'nye ob"emy faz v martensitnom prevrashchenii austenita [Specific Volumes of Phases in Martensite Transformation of Austenite]. Moscow, Metallurgizdat, 1950, 48 pp. (In Russian)
18. Tikhonov A. N., Samarskii A. A. Equations of mathematical physics. New York, A Pergamon Press Book The Macmillan Co., 1963, xvi+765 pp.
