CC BY
20
УДК 531.3
Вестник СПбГУ. Сер. 1. Т. 3(61). 2016. Вып. 1
ЭЛЛИПТИЧЕСКИЙ ПОГРАНСЛОЙ В ОБОЛОЧКАХ ВРАЩЕНИЯ ПРИ УДАРНЫХ ПОВЕРХНОСТНЫХ ВОЗДЕЙСТВИЯХ НОРМАЛЬНОГО ТИПА
И. В. Кириллова, Л. Ю. Коссович
Саратовский государственный университет имени Н. Г. Чернышевского, Российская Федерация, 410012, Саратов, ул. Астраханская, 83
В работе получены асимптотические уравнения, описывающие эллиптический погранс-лой в оболочках вращения при действии ударных нормальных поверхностных нагрузок. Рассматриваемый вид погранслоя имеет место только при поверхностных и торцевых нагрузках нормального типа и сосредоточен в малой окрестности условного фронта поверхностных волн Рэлея. Используются асимптотические методы в комбинациях с символическим методом Лурье. Для погранслоя получены эллиптические уравнения, описывающие поведение погранслоя по толщине оболочки, и гиперболические уравнения, задающие граничные условия на лицевых поверхностях. Библиогр. 18 назв. Ил. 1.
Ключевые слова: оболочки вращения, асимптотические методы, ударные нагрузки, символический метод Лурье, волна Рэлея, погранслой.
Настоящая работа посвящена построению погранслоя в окрестности условного фронта поверхностной волны Рэлея при ударных нормальных поверхностных нагрузках. Рассматриваемый погранслой используется в общей асимптотической схеме расчленения нестационарного напряженно-деформированного состояния (НДС) на составляющие с различными показателями изменяемости [1-7], когда на различных участках фазовой плоскости применяются безмоментная и изгибная составляющие теории Кирхгофа—Лява, квазистатистический погранслой типа Сен-Венана, уравнения квазиплоской и квазиантиплоской задач теории упругости, а также погранслои разных типов.
Представленный вывод асимптотически оптимальных уравнений рассматриваемого погранслоя обобщает исследования, проведенные для случаев полуплоскости, слоя и цилиндрической оболочки [8-17]. Особенностью данного исследования является характер выведенных асимптотических уравнений: эллиптические уравнения описывают изменение погранслоя по толщине, а гиперболические уравнения задают граничные условия и характеризуют волновой процесс на поверхности оболочки.
Рассмотрим оболочку вращения относительно криволинейных координат (а, в, г), изображенную на рисунке. Здесь а — длина дуги вдоль образующей срединной поверхности, в — угол в окружном направлении, 2 — внешняя нормаль к срединной поверхности.
Если мы обозначим напряжения и перемещения как а^ и VI (г = 1, 2, 3), граничные условия на лицевых поверхностях запишутся в виде
азз = -Р (а,г), а 13 =0 при 2 = ±Н, (1)
где Н — полутолщина оболочки. Также рассмотрим однородные начальные условия
vi = ^-=0 (г = 1,2,3) при * = 0. (2)
(¡5 Санкт-Петербургский государственный университет, 2016
Оболочка вращения
Разрешающие динамические уравнения теории упругости для осесимметричного случая запишутся в виде
_2 д2г'1 <92гч 1 д21>3 2/э(1+ г/) а2гч
~ ~ ~~^ О- ]~ч ^ ^ 7-Т . О Ь
да2 дг2 1 — дадг
Е
дг2
+
г ( _2 д2у\ а2гч 2/э(1+ г/) а2г>Л
да2 дг2
К
Е
дг2 )
2 В' ду1
+ I —+ —
. К\ К2
1 1
ду-
3
дгч /3-4 V 1
~д7 + И -2г/Д^ + \-2vRq I да
1 д2у1 д2у1
1 — 2^ дадг ' да2
С ^ д2Уз
, __2Э2г>3 2р(1 + гу)д2Щ 1
I .л I 96 0 0 |
дг2
Е
дг2
+
К1
да2
_2 <92г>з 2/э(1 + г/) 92г'з \ 3 - 4г/ 1 дгч
дг2
Е
дг2 ) 1 — 2^ К1 да
+
1 В' дгч В' д'о-з + 1 - 2г/"В + ~~В~д^
т
2
1 1 \ дп Я~ +Ж )~[Гг
'3
(3)
где р — модуль Юнга, коэффициент Пуассона и плотность материала оболочки, В — расстояние до оси вращения, т2 = (1 — 2v)/(2 — 2v), К1, К2 — радиусы кривизн, а уравнения связи между напряжениями и перемещениями имеют форму
Е
&11
&33
дг> 1 V дг>з г дг>1
V В'
2(1 + V)т2 V дх 1 — V дС К1 дх 1 — V В
У1 +
,1 V 1
+ +
К1 1 — vК2
Е
V ду1 дуз
2{1 + 1у)&2\1~1ь+~д7 1 — и К\ дх ' 1 — V В
V г ду-\ V В'
+ ---+
__
+ 1 -Дй! + д2
Е
&13 =
дг>1 дг>з г>1 г дг>з
2(1 + V) V дг дх К1 К1 дг ) ' Введем безразмерные переменные
У3
У3
(4)
е а г
Ц = Тг' Q=h'
гс2
Е
Т С2 у 2(1 + г/)р'
(5)
Тогда разрешающие уравнения (3)-(4) запишутся в виде
_2d2v\ d2vi 1 d2v3 d2v\ ( ( -2^1 d2v\ d2v\\
Ж ~W + ~W + 1 - 2v д£д( ~ 1h2 + 71 R[ V '86 d^ + ~дС ~ ) + _2 B' dv1 f 1 1 \ dv1 f 3 - 4v 1 1 1 \ dv3
1 д2v1 )2v3 _2 d2v3 д2v3 Z ( ))2v3 _2 )2v3 д2v3\
1 d2v\ d2v3 _2d2v3 d2v3 ( ( d2v3 _2d2v3 d2v3\
1 - 2v д£д(+ ~W + ж + v '+ т ~dC~W2)
3 - 4v 1 dv1 1 B' dv1 B' dv3 _2 f 1 1 \ dv3
E f dv1 v dv3 f Z dv1 v B' f 1 v 1 \ E ( v dvi )v3 ( v Z dvi v B'
a33= 2(1 + v)se2h +~dZ+V\ Т^ЖЖ + T^~B Vl+
E idvi dv3 (vi Z dv3
ai3 = 2(1 + v)h V + Ж ~ v \R~i + ЖЖ
(7)
где малый параметр п = h/R, R — характерное значение радиуса кривизны. Здесь обозначения Ri используются для приведенных значений радиусов кривизн (R* = Ri/R, звездочка опускается).
Анализ разрешающих уравнений показывает, что рассматриваемый тип НДС, который строится при единичных значениях показателей изменяемости и динамичности, может быть классифицирован как коротковолновый и представлен с погрешностью О(п) в виде комбинации симметричной и антисимметричной составляющих, каждая из которых имеет вид
аы = E(a°kl + n^h), vm = Rn(v°m + n^J, k,l,m = 1, 2, 3 (8)
где индекс «0» обозначает асимптотически главную компоненту (симметричную или антисимметричную относительно срединной поверхности), асимптотически второстепенная компонента обозначается индексом «1» (соответственно антисимметричная или симметричная). Тогда разрешающие уравнения для основной компоненты могут быть переписаны в форме
-2 д2v1 d2v1 1 д2v3 d2v1 -2 B' dvh
m ~W + + 1 - 21/ dZdZ ~ 1h2 + ^ ~в Ж =
1 д2v1 )2V3 _2 )2V3 д2v3 f 1 B' дv1 B' )V3\
+ "57Г --515-+ »i i-TT.^^r + ^^r =0;
(9)
1 - 2vдZ,дZ дZ2 дZ2 дт2 ' \,1 - 2v B дZ B ))Z
E f )v1 v )v3 v B'
an = 2(1 + !/)£e2h [Ж + Ж + Г]ТЖ~в""1 E f v )v1 )v3 v B' ,
2(1 + v)&2h \ 1 - v )Z )Z 1 - v B 1 E ( )V1 )V3 °13 = 2(1 + v)h + Ж
Вестник СПбГУ. Сер. 1. Математика. Механика. Астрономия. Т. 3 (61). 2016. Вып. 1
141
(11)
Построение уравнений погранслоя может быть выполнено с помощью символического метода Лурье. Выпишем систему уравнений (9)—(10) в операторной форме:
^ + («г2з? - а? + а£) ..1 + («Г2 =
Е (я , V сЬз , V В'
ап = оп—I—Г^Г + 1--^ + --в1'1
2(1 + v)ж2h\ 1 — V ¿С 1 — V В
Е ( V ¿юз V В' ,
а33 = о?1—I—Г^Г 1-+ 17 + --77^1 , 12
2(1 + v)т2h у 1 — V ¿С 1 — vB 1
Е (¿Ю1 , я
а13 = 2(1+7)^
где операторы , дт выражаются следующим образом: д% = д/ддт = д/дт.
Решим систему (11), рассматривая ее как систему обыкновенных дифференциальных уравнений:
Ю1 =де е-*(1-<)а1 В1 + де е-*(1+<)а1 Б2 + г^е-^1-^?1 Вз — г^е-^1^?1 ВА, ю3 =га1е-1(1-^)а1 В1 — га1е-1(1+^)а1 В2— (13)
где Иг —произвольные операторы, а2 = <9| — аг2д2 + (З2 = <9| — д2 — зг2д2+г]^.
Проанализируем символическое решение, основываясь на прямой аналогии символического метода с методом интегральных преобразований. Решение (13) соответствует изображению двойного интегрального преобразования Лапласа (по времени) и Фурье (по продольной координате) при дт ~ в, ~ гх, где в и х — параметры преобразования Лапласа и Фурье соответственно. В рассматриваемом решении слагаемые с множителями (1 — С) в показателях степени экспоненты соответствуют волнам, инициируемым лицевой поверхностью С =1 и распространяющимся в отрицательном направлении оси С, а слагаемые с множителями (1 + С) —волнам, инициируемым лицевой поверхностью С = —1 и распространяющимся в положительном направлении оси С.
Подставляя решение (13) в граничные условия (1), получим следующую систему алгебраических уравнений для
71^1 + 11е-'2олВ2 + г ^ + гл^ ^ - г (д^ + г]^ В4 = -Р,
а1 В1 — гд% але-12а1 В2 — 121В3 — ^е-2^ В4 = 0
72е-»2а1с£,1 + 72д2 + г ^ + ^^ ^е-^Бз ~ г ^ + = -Р,
(14)
гд% а1 е-аа1 В1 — гд% а1В2 — ^е-^1 В3 — = 0.
Для упрощения вычислений пренебрежем членами с множителями е-г2а1 и е-г2в1 в граничных условиях (14). Это соответствует, в конечном итоге, пренебрежению
членами с множителями ехр(—2^/1 — ее2ее2Ех) и ехр(—2-^/1 — ее2ях) (сн ~ скорость поверхностных волн Рэлея) в асимптотике для первого корня уравнения Рэлея—Лэмба при х ^^ [18]:
= зедХ - ВшхехР ~ + ВшхехР ~ яР&яХ^ >
Бш = 2
же жжд 4жд
+
1 — 1 — ж е 2 — Ж Е
-1
СЕ ^ 1
, аед = — <1. С
(15)
Коэффициент же является корнем известного трансцендентного уравнения
(2 — жЕ)2 — 4^1 — ж^ 1 — ж2жЕ = 0. (16)
Тогда уравнения (14) могут быть преобразованы к виду
Б'4
(17)
¡д^а1В1 — з = 0,
71^2 - « (де + п^ А А = -Р,
—гд5 а1П2 — ~{2хВа = 0.
Рассмотрим большие значения времени, когда т ^ 1 (т.е. когда фронты волн проходят расстояние, много большее толщины оболочки: т = 0(Т),Т ^ 1). Введем малый параметр е = 1/Т и масштабированные переменные у, то:
€ — жет /Ло\ У=---, т0 = £т, (18)
позволяющие построить прифронтовую асимптотику в дальнем поле волны Рэлея. Введем новые операторы ду, дТ0:
ду = д/ду, дт0 = д/дто. (19)
Раскладывая уравнения (14) и выражения для перемещений и напряжений в ряд по малому параметру е, оставим в них асимптотически главные части. Затем вернемся к исходным операторам д%, дт и, соответственно, к дифференциальной форме искомых уравнений. Введем объемные и сдвиговые потенциалы Ламе: ф1,ф2 и '1 ,ф2, где индексом «1» будем обозначать отрицательное, индексом «2» —положительное по С направление распространения возмущений. Тогда потенциальные функции фг,фг в исходных размерных координатах будут определяться эллиптическими уравнениями
—+ —=0
с2)^а2 дг2 ' (зо)
Ся\ д2'^ 1 д 2'г
1 --т) + = 0, г = 1,2,
с2 / да2 дх2
при следующих граничных условиях:
д2Х1 2 д2Х1
_ _ г» _
dt2
д2Х2 dt2
1 да2 д2Х2 1 до2
- K
- еея Кш
&дх1 ' В да дх2
asrВЫЪ дР а да' же Вш b дР
B да
а
д2фг дадг д2ф2 дадг
а2 В дг 1 В'дф2
+
д2'фг |
да2 д2ф2
~9■
да2
"R,
1- mR,
В' дф 1 В' <9V>2
K„ =
да
= h,
= h,
h,
h,
B„
3
1
2
где х® = д'фг/дг. Перемещения и напряжения выражаются через потенциальные функции по формулам
^ _ дер 1 ^ дср2 | |
да да дг дг '
дф2 дф\ дф2
=--1------
дг дг да да
К
~В
(Ф1 + Ф2),
&зз = -3
д2
д2ср2
да2 да2
<Р1
+
д2ф1 д2ф2
дадг дадг
~в
&13
д 2fi д2^2 f д2'ф1 д 2ф2 + т^г- ~9 I -^г +
дадг дадг
да2 да2
д^р1 ¿¥1
да да дг дг
В V да да J
(22)
В качестве примера проанализируем распространение волны Рэлея по лицевой поверхности цилиндрической оболочки, для которой в граничных условиях (1) Р(а,г) = р6(а)И(г), где 6 —дельта-функция, И — функция Хевисайда, а р — заданная постоянная. Использование аппарата интегральных преобразований Фурье по пространственной переменной а и Лапласа по времени дает следующее ожидаемое выражение для перемещения юз на поверхности г = h, задающее граничные условия для малых окрестностей условных фронтов поверхностной волны Рэлея:
V3 =
жеВш p
- 3ln-
2(h + г)2 + (а - cRtf a2(h + г)2 + а2
-grln-
2(h + z)2 + (a + cRtf a2(h + г)2 + а2
+
+ln
b2(h + г)2 + (а - CRt)2 b2(h + г)2 + (а + cR)
2
b2(h + г)2 + а2
ln-
i2(h + г )2 + а2
(23)
Формула (23) иллюстрирует механизм распространения участка возмущения НДС в районе условного фронта поверхностной волны Рэлея (точнее, квазифронта). Поверхностная волна распространяется со скоростью сд вдоль границы и формирует участок возмущения в теле оболочки в малой окрестности рассматриваемого квазифронта.
1
2
а
R
Литература
1. Коссович Л. Ю. Нестационарные задачи теории упругих тонких оболочек. Саратов: Изд-во СГУ, 1986. 176 с.
2. Kaplunov J. D. On the quasi-front in 2D shell theories // C. R. Acad. Sci. Paris, France. Vol.313, N2. P. 731-736 (1991).
3. Kaplunov J. D., Kossovich L. Yu., Nolde E. V. Dynamics of thin walled elastic bodies. San Diego: Academic Press, 1998.
4. Шевцова Ю. В. Погранслой в окрестности квазифронта в трансверсально-изотропной цилиндрической оболочке // Проблемы прочности материалов и конструкций, взаимодействующих с агрессивными средами. Межвуз. науч. сб. Саратов: СГТУ, 2000. С. 114—117.
5. Коссович Л. Ю. Асимптотический анализ нестационарных, упругих волн в тонких оболочках вращения при ударных торцевых воздействиях / Л. Ю. Коссович, Ю. Д. Каплунов // Известия Саратовского университета. Саратов, 2001. T. I, №2. С. 111—131.
6. Кириллова И. В. Области согласования динамического погранслоя и коротковолнового высокочастотного приближения // Математическое моделирование и управление в технических системах: сб. науч. тр. Сарат. гос. агроинженер. ун-т. Саратов: Изд-во СГУ, 1998. С. 3—11.
7. Nigul U. Regions of effective application of the methods of three-dimensional and two-dimensional analysis of transient stress waves in shells and plates // Int. J. Solids Struct. 1969. Vol.5, issue 6. P. 607— 627.
8. Каплунов Ю.Д., Коссович Л.Ю. Асимптотическая модель для вычисления дальнего поля волны Рэлея в упругой плоскости // Докл. АН. 2004. Т. 395. №4. С. 482-484.
9. Kaplunov J., Kossovich L., Zakharov A. An explicit asymptotic model for the Bleustein—Gulyaev wave // C.R. Mecanique. 2004. Vol.332. P. 487-492.
10. Коссович Л. Ю. Поле Релея в бесконечном упругом слое / Л. Ю. Коссович, А. Н. Кушекка-лиев // Математика, механика. Саратов: Изд-во СГУ, 2003. Вып. 5. С. 159-161.
11. Ковалев В. А., Коссович Л. Ю., Таранов О. В. Поле Рэлея в задаче Лэмба для цилиндрической оболочки // Известия вузов. Северо-Кавказский регион. Естественные науки. Спец. вып. 2004. С. 52-54.
12. Ковалев В. А., Коссович Л. Ю., Таранов О. В. Дальнее поле волны Рэлея для упругой полуполосы при действии торцевой нагрузки // Изв. РАН. МТТ. 2005. №5. С. 89-96.
13. Кушеккалиев А. Н. Решение задач о распространении волн в трансверсально-изотропной цилиндрической оболочке при нормальных воздействиях // Механика деформируемых сред. Саратов: Изд-во СГУ, 2002. Вып. 14. С. 106-115.
14. Кушеккалиев А. Н. Волны типа Релея в полубесконечной пластине при нормальном воздействии поперечного типа // Механика деформируемых сред. Саратов: Изд-во СГУ, 2004. C. 66-73.
15. Коссович Л. Ю. Анализ приближений в задаче Лэмба для бесконечного упругого слоя / Л. Ю. Коссович, А. Н. Кушеккалиев // Известия вузов. Северо-Кавказский регион. Естественнные науки. 2003. №9. С. 10-22.
16. Коссович Л. Ю. Расчленение нестационарного НДС в задаче Лэмба для бесконечного слоя на составляющие с разными показателями изменяемости / Л. Ю. Коссович, А. Н. Кушенкалиев // Тр. III Всероссийской конференции по теории упругости с международным участием. Ростов-на-Дону, 2003. C. 232-234.
17. Ковалев В. А., Таранов О. В. Расчленение нестационарного НДС цилиндрических оболочек при ударных торцевых воздействиях нормального типа // Смешанные задачи механики деформированного тела: Материалы V Рос. конф. с международным участием. Саратов: Изд-во СГУ, 2005. С. 191-193.
18. Айнола Л. Я., Нигул У. К. Волновые процессы деформации упругих плит и оболочек // Изв. АН ЭССР, сер. физ.-мат. и техн. наук. 1965. Т. 14, №1. С. 3-63.
Статья поступила в редакцию 22 октября 2015 г. Сведения об авторах
Кириллова Ирина Васильевна — кандидат физико-математических наук, доцент; nano-bio@sgu.ru Коссович Леонид Юрьевич —доктор физико-математических наук, профессор; nano-bio@sgu.ru
ELLIPTIC BOUNDARY LAYER IN SHELLS OF REVOLUTION UNDER NORMAL SHOCK SURFACE LOADING
Irina V. Kirillova, Leonid Yu. Kossovich
Saratov State University, ul. Astrakhanskaya, 83, Saratov, 410012, Russian Federation; nano-bio@sgu.ru
Asymptotic equations for the elliptic boundary layer in shells of revolution at shock normal surface loading are derived in the present paper. Considered boundary layer type has place only at surface and edge shock loading and is localized in the small vicinity of the front of the Rayleigh surface waves. There are
used asymptotic methods together with A. I. Lourier's symbolic integration. There are obtained elliptic equations discribing the boundary layer behaviour with respect to the normal coordinate and hyperbolic equations defined the boundary layer on the shell faces. Refs 18. Figs 1.
Keywords: shell of revolution, asymptotic methods, shock loading, symbolic method Lourier, Rayleigh wave, boundary layer.
References
1. Kossovich L. Yu., Non-stationary Problems in Theory of Thin Elastic Shells (Izd-vo Saratov Univ., Saratov, Russia, 1986) [in Russian].
2. Kaplunov J. D., "On the quasi-front in 2D shell theories", C. R. Acad. Sci. 313(2), 731-736 (Paris, France, 1991).
3. Kaplunov J.D., Kossovich L.Yu., Nolde E. V., Dynamics of thin walled elastic bodies (Academic Press, San Diego, 1998).
4. Shevtsova Yu. V., "Boundary layer in the vicinity of quasi-front in transversal isotropic cylindrical shell", On the strength of the elements of constructions under the action of loads and working media (Izd-vo Saratov Technical Univ. Saratov, Russia, 2000, 114-117) [in Russian].
5. Kossovich L. Yu. "Asymptotic Analysis of the Non-stationary Elastic Waves in Thin-walled Shells of Revolution Initiated by the Edge Shock Load", J. Izvestiya Saratov University 1(2), 111-131 (Izd-vo Saratov Univ., Saratov, Russia, 2001) [in Russian].
6. Kirillova I. V., "Overlap regions of the boundary layer and short-wave high-frequency approximation", Mathematical simulation and control in the technical systems (Izd-vo Saratov Agrarian Univ., Saratov, Russia, 1998, 3-11) [in Russian].
7. Nigul U., "Regions of effective of application of the methods of 3D and 2D of analysis of transient of stress waves in shells and plates", Int. J. Solids Struct. 54, issue 6, 607-627 (1969).
8. Kaplunov Yu.D., Kossovich L.Yu., "An asymptotic model for the Rayleigh wave far-field in the case of an elastic half-plane", J. Russian Phys. Dokl. 395(4), 482-484 (2004) [in Russian].
9. Kaplunov J., Kossovich L., Zakharov A. "An explicit asymptotic model for the Bleustein—Gulyaev wave", C.R. Mecanique 332, 487-492 (2004).
10. Kossovich L. Yu., Kushekkaliev A. N., "Rayleigh's field in the infinite elastic layer", Mathematics. Mechanics issue 5 (Izd-vo Saratov Univ., Saratov, Russia, 2003, 159-161) [in Russian].
11. Kovalev V.A., Kossovich L.Yu. Taranov O.V., "Rayleigh's field in Lamb's problem for the cylindrical shell", J. Izvestiya VUZov. Severo-Kavkazskii region. Natural sciences. Special issue (Rostov-na-Donu, Russia, 2004, 52-54) [in Russian].
12. Kovalev V.A., Kossovich L.Yu., Taranov O. V., "Rayleigh wave far-field in the case of an elastic semi-strip under the edge loading", J. Solid Mechanics 5 (Moscow, Russia, 2005, 89-96) [in Russian].
13. Kushekkaliev A. N., "Solution of the problems about the propagation of waves in the transversal-isotropic cylindrical shell under the normal loading", Continuum mechanics issue 14 (Izd-vo Saratov Univ., Saratov, Russia, 2002, 106-115) [in Russian].
14. Kushekkaliev A. N., "Rayleigh's wave in the semi-infinite plate under the transverse normal loading", Continuum mechanics (Izd-vo Saratov Univ., Saratov, Russia, 2004, 66-73) [in Russian].
15. Kossovich L. Yu., Kushekkaliev A. N., "The analysis of approximations in Lamb's problem for the infinite elastic layer", J. Izvestiya VUZov. Severo-Kavkazskii region. Natural sciences 5 (Rostov-na-Donu, Russia, 2003, 10-22) [in Russian].
16. Kossovich L.Yu., Kushekkaliev A.N., "The separation of non-stationary SSS in Lamb's problem for the infinite layer into the components with the different values of variability", The Theory of Elasticity. Proceedings of the III Russian conference (Izd-vo Rostov Univ., Rostov-na-Donu, Russia, 2003, 232-234) [in Russian].
17. Kovalev V. A., Taranov O. V., "The separation of non-stationary SSS of cylindrical shells under normal type edge loading", Mixed Problems of Mechanics of Solid,. Proceedings of the V Russian conf. with Intern. Participation (ed. Acad. N.F. Morozov; Izd-vo Saratov Univ., Saratov, Russia, 2005, 191-193) [in Russian].
18. Ainola L. Ya., Nigul U.K., "Wave Strain Processes in Elastic Plates and Shells", Izv. AN ESSR 14(1), 3-63 (1965) [in Russian].
