CC BY
43
Егорова О.В., Жаворонок С.И., Курбатов А.С. О вариационных уравнениях расширенной теории N-го порядка упругих оболочек и их приложении к некоторым задачам динамики // Вестник ПНИПУ. Механика. - 2015. - № 2. -С. 36-59. DOI: 10.15593/perm.mech/2015.2.03
Egorova O.V., Zhavoronok S.I., Kurbatov A.S. An application of various n-th order shell theories to normal waves propagation problems. PNRPU Mechanics Bulletin. 2015. No. 2. Рр. 36-59. DOI: 10.15593/perm.mech/2015.2.03
ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2015
PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
DOI: 10.15593/perm.mech/2015.2.03 УДК 539.3
0 ВАРИАЦИОННЫХ УРАВНЕНИЯХ РАСШИРЕННОЙ ТЕОРИИ N-ГО ПОРЯДКА УПРУГИХ ОБОЛОЧЕК И ИХ ПРИЛОЖЕНИИ
К НЕКОТОРЫМ ЗАДАЧАМ ДИНАМИКИ
О.В. Егорова1, С.И. Жаворонок2, А.С. Курбатов3
1 Московский авиационный институт (Национальный исследовательский университет)
2
Институт прикладной механики Российской академии наук (ИПРИМ РАН)
3 Исследовательский центр им. М.В. Келдыша
О СТАТЬЕ
АННОТАЦИЯ
Получена: 27 марта 2015 г. Принята: 6 мая 2015 г. Опубликована: 30 июня 2015 г.
Ключевые слова:
оболочки толстостенные, аналитическая механика континуальных систем, формализм Лагранжа, слои упругие, волны нормальные, скорости групповые, явление «обратной волны», теория пластин N-го порядка
Одной из основных задач теории оболочек средней и большей толщины является моделирование высокочастотных колебаний и распространения волн. Во многих практически важных случаях, возникающих в инженерной практике аэрокосмической отрасли, такие процессы не могут быть адекватно описаны традиционными теориями, основанными на методе гипотез. Требуется разработка моделей, учитывающих, помимо векторов перемещения и поворота, высшие степени свободы и, следовательно, взаимодействие движений в тангенциальном и трансверсальном направлениях.
В данной статье рассматриваются некоторые варианты теории оболочек N-го порядка, принадлежащей к классу теорий И.Н. Векуа и предложенной А.А. Амосовым. Применяемые варианты теории основаны на вариационном формализме аналитической механики континуальных систем со связями. Модель оболочки представляет собой поверхность с заданными на ней множеством обобщенных координат (переменных поля) и поверхностной плотностью лагранжиана. Уравнения движения оболочки имеют вид обобщенных уравнений Лагранжа второго рода. Удовлетворение краевых условий на лицевых поверхностях обеспечивается в расширенном варианте теории. Краевые условия переносятся на базисную поверхность и записываются относительно переменных поля, образуя уравнения связей в континуальной механической системе.
Поставлена однородная задача о распространении нормальных волн в плоском упругом слое и рассмотрено ее решение на базе теории пластин N-го порядка. Продолжено исследование описания теориями низших порядков второй распространяющейся продольной моды, отличающейся эффектом «обратной волны» - различными знаками фазовой и групповой скоростей в диапазоне малых волновых чисел. На базе решения спектральной задачи получены формы второй распространяющейся продольной моды при различных значениях волнового числа, изучена эволюция формы, приводятся оценки аппроксимации формы волны на базе теорий пластин различного порядка. Сделаны выводы о точности описания «обратной волны» расширенными и элементарными теориями пластин, основанными на ортогональных разложениях вектора перемещения.
© ПНИПУ
© Егорова Ольга Владимировна - кандидат физико-математических наук, доцент, e-mail: janus_olga@mail.ru Жаворонок Сергей Игоревич - кандидат физико-математических наук, доцент, e-mail: zhavor71@mail.ru Курбатов Алексей Сергеевич - кандидат технических наук, старший научный сотрудник, e-mail: defunt@inbox.ru
Olga V. Egorova - PhD in Physical and Mathematical Sciences, Associate Professor, e-mail: janus_olga@mail.ru Sergey I. Zhavoronok - PhD in Physical and Mathematical Sciences, Associate Professor, e-mail: zhavor71@mail.ru Alexey S. Kurbatov - PhD in Technical Sciences, Senior Researcher, e-mail: defunt@inbox.ru
36
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
THE VARIATIONAL EQUATIONS OF THE EXTENDED N’TH ORDER SHELL THEORY AND ITS APPLICATION TO SOME PROBLEMS OF DYNAMICS
O.V. Egorova1, S.I. Zhavoronok2, A.S. Kurbatov3
Moscow Aviation Institute (National Research University)
Institute of Applied Mechanics of Russian Academy of Sciences
3
Federal State Unitary Enterprise Research Centre named after M.V. Keldysh
ARTICLE INFO
ABSTRACT
Received: 27 March 2015 Accepted: 6 May 2015 Published: 30 June 2015
Keywords:
thick shells, shell theories of Vekua type, analytical mechanics of continua,
Lagrange formalism, elastic layer, normal waves, group velocities, inverse wave effect, plate theory of N'th order
One of topical problems of thick shells theory consists in an accurate modelling of high-frequency shell vibrations and wave propagation, as well as in the application of developed mathematical models to engineering problems of structural dynamics. Many observed processes of both periodic and non-steady shell dynamics cannot be adequately simulated on the groundwork of most traditional hypothesis-based shell theories. Therefore the engineering practice requires new shell models considering higher degrees of freedom besides the displacement and rotation vectors and taking into account the interference of tangential and transversal motions.
Here a variant of the N'th order shell model of I.N. Vekua type is analyzed. The proposed model improves the A. Amosov's shell theory; it is based on the Lagrange's formalism of analytical mechanics extended to continuum systems, on the dimensional reduction approach, and on the biorthogonal expansion technique. The shell is finally represented by a material surface framed by a set of generalized coordinates (or field variables) and a scalar generator function - the surficial density of the Lagrangian. The dynamic equations for the shell are formulated as generalized Lagrange equations of the second kind.
A propagation of normal waves in a plane elastic layer is considered, the corresponding homogeneous problem of N'th order plate theory is formulated, and its solution analysis is extended. In particular, the investigation of the second longitudinal propagating mode and its modelling by approximated plate theories are expanded. The main feature of the second mode is the “inverse wave” effect, i.e. the opposite signs of the phase and group velocities at small wavenumbers. It was shown earlier that the accurate dispersion curves for both phase frequency and group velocity can be obtained on the basis of the fourth- and fifth-order theories. Using the known solution of the spectral problem of N'th order plate theory the eigenfunctions are obtained. The second propagating mode evaluation at the smallest wavenumbers is modelled, the forms of the second propagating wave are investigated, and the approximation of the waveforms at some wavenumbers by plate theories of various orders is estimated. The accuracy of the approximation of the “inverse wave” by some applied theories of various orders based on the orthogonal expansions of displacement vector is discussed.
© PNRPU
Введение
Одним из основных элементов конструкций авиационно-космических систем различного назначения является тонкостенная оболочка, подвергающаяся, как правило, интенсивному силовому и тепловому динамическому воздействию со стороны окружающей среды, собственных энергетических установок и оборудования. Разработка перспективной аэрокосмической техники требует качественного анализа особенностей динамических процессов, протекающих в тонкостенных элементах конструкций при высокочастотных вибрациях и распространении волн. В ряде случаев традиционно используемые модели элементов такого типа, адекватно описывающие низкочастотные процессы и притом обладающие компактной формулировкой, становятся непригодными. Такие ситуации возникают, например, при описании высокочастотных вибраций, а также при моделировании распространения и интерференции волн. По мере ускорения процессов и повышения диапазона частот сначала классические теории, основанные на гипотезах Кирхгоффа-Лява и их уточнения, принимающие во внимание инерцию вращения нормального эле-
37
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
мента, а затем и теории класса Тимошенко-Райсснера-Миндлина, учитывающие сдвиг в трансверсальной плоскости, становятся непригодными. Дальнейшее уточнение модели, как правило, диктует отказ от метода гипотез и перехода к редукции трехмерной задачи механики деформируемого твердого тела.
Исчерпывающее описание как высокочастотных колебаний, так и процессов распространения волн в тонкостенных системах на базе «двумерных» моделей на текущий момент не достигнуто [1, 2]. Обзоры основных работ, посвященных так называемым «неклассическим» моделям оболочек и их приложению к задачам о колебаниях и опубликованных до 1970-х гг., содержатся в работах [3] и [4]. В течение последних 15-20 лет интерес к подобным теориям оболочек возобновился. В работах [5, 6] продемонстрировано существенное влияние трансверсальных нормальных напряжений на динамику слоистых оболочек; статья [7] содержит достаточно подробный аналитический обзор моделей высшего порядка тонкостенных слоистых структур, разработанных до 2002 г, и их практических реализаций на базе метода конечных элементов. Результаты работ [5-7] подтвердили выводы Койтера [8] о необходимости совместного учета трансверсальной сдвиговой и нормальной податливости при уточнении классических моделей оболочек, что представляется возможным в рамках теорий высших порядков. В работе [9] исследовано явление так называемого «толщинного» запирания в теориях оболочек низших порядков, возникающего при предположении об отсутствии деформации изменения толщины оболочки, и показана возможность его исключения при учете высших степеней свободы. В статье [10] обоснована необходимость применения теорий с порядком аппроксимации до четвертого в случае высокочастотных колебаний пластины вне зависимости от толщины. В работе [11], в свою очередь, на основе теорий высших порядков показано, что при высокочастотной вибрации пластины существуют решения типа пограничного слоя, локализующегося вблизи лицевых поверхностей. Эффективность моделей, основанных на разложении по системам специальных функций, продемонстрирована в статье [12] на примере круговой пластины. «Трехмерные» модели на базе степенных разложений были использованы для решения задачи о колебаниях конических оболочек [13]. Другим примером необходимого применения теорий высшего порядка является описание статического деформирования [14], колебаний и потери устойчивости функциональноградиентных пластин [15, 16] и оболочек [17, 18], применяемых в системах тепловой защиты летательных аппаратов при комбинированном воздействии. В [19-22] теории высших порядков использованы при решении обратных задач нестационарной динамики вязкоупругих композиционных тонкостенных элементов.
Помимо исследования высокочастотных колебаний, теории оболочек, учитывающие высшие степени свободы, требуются при решении таких задач, как изучение краевых и кромочных волн; классическая теория пластин способна описать только первую (фундаментальную) краевую волну, тогда как для учета высших волн необходима трехмерная постановка задачи [23]. Следовательно, усовершенствование моделей деформирования оболочек, учитывающих высшие степени свободы и взаимодействие тангенциальных и трансверсальных движений, является актуальным.
Проблема приведения трехмерной задачи к двумерной оболочки не имеет единственного решения [24] и базируется на различных математических методах [25]. «... Вопросы построения приближенных теорий колебаний стержней, пластин и оболочек. представляют собой часть общей проблемы построения моделей механики сплошных сред, получившей развитие в работах Л.И. Седова» [4, с. 6]. Обзору ранних решений проблемы при-
38
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
ведения посвящены работы [4, 24]. Обычно выделяют методы асимптотического интегрирования уравнений теории упругости [26, 27], прямого построения теории материальных поверхностей [28, 29], метод пространственной редукции трехмерной задачи [14-22, 30-43] на базе разложений по степенным [14-18, 30] либо специальным функциям [5-7, 9, 10, 19-22, 33-43], непрерывным либо кусочно-гладким (в случае слоистых оболочек) [57], в комбинации с методом Галеркина ([31, 32, 36-38, 42] и др.) или вариационным подходом [например, 24, 33-35, 39-41 и др.].
В данной работе описано совместное приложение метода И.Н. Векуа [31] пространственной редукции трехмерной задачи теории упругости, редукционного подхода А.А. Амосова, использованного ранее в работах [37-43], и формализма аналитической механики континуальных систем [44]. Вариационный формализм построения трехмерной, т.е. способной описывать напряженно-деформированное состояние в областях неприводимости, теории А. А. Амосова оболочек N-го порядка был предложен в работах [45-49]. Модель оболочки формулируется в терминах двумерного многообразия как континуальная система, заданная конфигурационным пространством с N векторными переменными поля первого рода и определенной на нем поверхностной плотностью функционала Лагранжа. Уравнения движения оболочки строятся на базе формализма [45-48] как обобщенные уравнения Лагранжа второго рода относительно обобщенных усилий [45-48] или переменных поля первого рода [49]. Определяющие уравнения модели N-го порядка анизотропной упругой оболочки приведены в работе [46].
В цитируемых выше статьях [45-49] проблема удовлетворения краевых условий на лицевых поверхностях оболочки рассматривалась на базе «упрощенной», согласно терминологии И.Н. Векуа [31], или «элементарной» теории оболочек N-го порядка [37]. Вариационный метод редукции уравнений трехмерной задачи [44] обеспечивает удовлетворение краевых условий на лицевых поверхностях «в среднем» [37, 45]. В рамках проекционного метода [31] приведения трехмерной задачи теории упругости к двумерной задаче теории оболочек предложены различные приемы, обеспечивающие удовлетворение краевых условий на лицевых поверхностях: метод корректирующего решения [37] или «расширенная» формулировка теории оболочек [37], «метод нормированных моментов» И.Н. Векуа [31], использующего в качестве дополнительных неизвестных остаточные члены рядов компонентов вектора перемещения. Вариационный формализм аналитической механики континуальных систем [44], в свою очередь, позволяет рассматривать краевые условия, перенесенные на базисную поверхность оболочки с ее лицевых поверхностей в процессе редукции модели, в качестве уравнений дополнительных связей, накладываемых на переменные поля первого рода [44], голономных в случае кинематических краевых условий и неголономных в случае силовых краевых условий на лицевых поверхностях. Полученная таким образом вариационная задача с ограничениями в форме равенств решается методом множителей Лагранжа, причем уравнения движения, записанные относительно новых обобщенных сил, включающих слагаемые со множителями, идентичны по форме уравнениям элементарной теории N-го порядка. Разрешение уравнений связей относительно множителей Лагранжа позволяет получить новые определяющие уравнения, учитывающие влияние краевых условий на лицевых поверхностях. Можно показать, что изгибная жесткость пластины в рамках теорий низшего порядка совпадает с цилиндрической жесткостью классической теории без принятия гипотезы о трансверсальной изотропии материала [50, 51] либо применения методов штрафа [9] или преобразования определяющих уравнений [52]. Таким образом, устраняется погреш-
39
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
ность вариационного подхода, согласно [29], где указано, что «...трехмерная энергия учитывает энергию тех движений частиц оболочки, трехмерного тела, которые никак не сказываются на движении собственно двумерной оболочки» (29, с. 15), или эффект «тол-щинного» запирания [9].
Свойства модели, описываемой уравнениями теории N-го порядка [45, 46], исследованы в работах [41, 43, 47, 48, 54-57]. В работах [41] и [43] получены оценки сходимости приближенных решений различных задач статики на базе трехмерной модели N-го порядка к известным точным решениям. Работы [47, 48, 54-57], в свою очередь, посвящены решению различных модельных задач динамики. Так как «.наиболее простым методом исследования уравнений динамической теории стержней, пластин и оболочек является метод дисперсионных уравнений» ([4], с. 25), в работах [47, 54, 55] на базе теории пластин N-го порядка были построены дисперсионные уравнения задачи о распространении нормальных волн в плоском упругом слое. Исследована сходимость спектра при нулевом значении волнового числа (т. е. частот запирания распространяющихся нормальных мод). В статье [56] на основе решения спектральной задачи теории N-го порядка вычислены приближенные формы распространяющихся мод нормальных волн, соответствующие частотам запирания, и их погрешности относительно точных решений [58] по норме евклидова пространства. Кроме того, получены значения фазовых частот в точках равенства соответствующих фазовых скоростей и скорости дилатационной волны с1 [56], а также погрешности приближенных форм нормальных мод, соответствующих данным точкам, относительно точного решения (мод Гудьера-Бишопа) в среднем квадратичном. В работе [56] исследованы погрешности мод Ламе.
При исследовании свойств теорий пластин низших порядков особый интерес представляет их способность правильно описывать вторую нормальную продольную моду, распространяющуюся в плоском упругом слое и отличающуюся существованием диапазона волновых чисел [0, к* ] с положительной фазовой и отрицательной групповой скоростью [58]. Физическая возможность возбуждения «обратных» волн, доказанная экспериментально в работе [59], делает данную задачу весьма интересной как с точки зрения практических приложений, в частности, в акустике машин [60], так и в качестве тестовой для анализа различных приближенных моделей. В работах [47, 49] на базе элементарных теорий 1-5-го порядков, не учитывающих краевые условия на лицевых поверхностях в форме связей и использующих в качестве базисной системы полиномы Лежандра, построен начальный участок второй дисперсионной ветви, соответствующий «обратной» волне. Показано, что адекватное значение групповой скорости достигается на основе теорий не ниже 4-го порядка. В работе [57] исследовано описание формы второй распространяющейся моды продольных волн в области существования эффекта «обратной волны». Эволюция формы волны описана элементарной теорией пластин 5-го порядка.
Ниже приводится анализ второй распространяющейся моды, построенной на основе решения спектральной задачи в рамках элементарной и расширенной теории 2-5-го порядков. Получены зависимости компонентов вектора перемещения от безразмерной нормальной координаты в диапазоне малых волновых чисел, определяющие эволюцию формы «обратной» нормальной волны. Исследована сходимость решения по мере повышения порядка применяемой теории пластин. Полученные результаты являются основой моделирования распространения волн в тонких телах на базе моделей, являющихся частными случаями теории оболочек N-го порядка.
40
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
1. Формулировка элементарной теории оболочек N-го порядка
Параметризация геометрии оболочки и ее базисной поверхности
Пусть оболочка занимает область V е R3, д¥ = S± ® SB ; S± - лицевые, SB - боковые поверхности оболочки. Параметризация оболочки в рамках теорий, описываемых в терминах двумерного многообразия, заключается:
1) в параметризации базисной поверхности S0, являющейся геометрической моделью оболочки и определяющей формулировку двумерной начально-краевой задачи;
2) в параметризации области V е R3, являющейся ключевым пунктом решения «проблемы приведения» [24] трехмерной задачи к двумерному эквиваленту.
Параметризация $0. Введем криволинейные координаты ЕP, p = 1,2 :
VMo е So r(Mо) = r (^Р) = x' (^Р)e,, (1)
где X - декартовы координаты точки М0. Здесь и далее по повторяющемуся индексу подразумевается суммирование; латинские индексы пробегают значения 1, 2, 3; греческие - значения 1, 2. Базис ra и метрика а касательного расслоения TMS0, единичная нормаль S0 и тензор кривизны b определяются обычным образом [31]:
ra (Мо) = даr(ЕР) = J(e,, J( = дх'/дЕ“ ,
aap = ra' rp; rP= a“Pra , ra = Oap^, a “4p = 5P ; a = det (ap) ,
n = (r1 X r2 )Ha ; Ч=-[дП/дЕ“]' rp; Cap= b( 4.
Определим также среднюю H (M0) = 1 b( и гауссову K (M0) = det(b(p) кривизны S0.
Параметризация V. Введем нормальную координату Е3 е [h-, h+ ], так, что
VM е V \ S0 R (M ) = r (M0 ) + Е3п (M0), (2)
где h = h+ - h_ - толщина оболочки. При h = 2(h+ - h_) = 0 S0 - срединная поверхность. Базис Ra, n пространственной системы координат Е1, Е 2, Е3 определен так [31]:
R a
5R(M0,Е3= A(prp(M0); A(p = б£ _Е3Ь(; |A(p| = 1-2ГЯ + (')3,
ra= Apa Rp, =
1 - 2Е3H + ( )2 K 1 1 [8( - Е3 (2H8( - bpa)
Условия однозначности параметризации (2) приведены в [47]. Там же введена система координат г|\ Ч п3, голономным базисом которой являются векторы ra, n :
r =дR _a&paR , ,p. = 5^ p = 8^, , 3
^ A.a ^ a , Aa. ^ a ; ,
a
дЧ дП дЕ'3
V = R a
д
- + n-
д
дпа дП
дп°
=ra a pa
д
- + n-
дЕ
д
3
(3)
41
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
С учетом (3) дисторсия и плотность потенциала записываются в виде [46, 47]
d = (i- ba|A )Г “ + (( + bpua)RP" + d3u„nr “ + д3«3ПП Э3 = , (4)
2W = о: d = s“p (Vp«a - b^) + s3p (Vpu + bp4) + + s33d3u, (5)
о = s “pra Rp + s 3pnRp + s a3ra n + s 33nn, (6)
где Va - ковариантная производная на двумерном многообразии S0; о - симметричный тензор напряжения; ui - компоненты вектора перемещения в базисе r“ , п , не зависящем
г 3
от координаты q .
Формулировка модели N-го порядка анизотропной оболочки
В соответствии с подходом [45-49] модель оболочки N-го порядка - континуальная система на S0 е R3, заданная конфигурационным пространством
Q
{u(k], i = 1...3, k = 0...N, u(kR= u(k 1 (Mo,t), Mо eSo,
(7)
где t e [t0, t1 ] eR+u{0} - время, u(k) - переменные поля I рода [44],
u (M, t) = u(k ) (M0, t )P(k) (C) = [ui*) (M0, t ) (M0)+u3k ) (M0, t R (M0 )]], k) k u(k l(M0, t ) = ( u, (M0, c, t ),p( k > (c))1,
(8)
(9)
где C = 2(q3 -h)/h e [-1,1] - безразмерная нормальная координата [45, 46]; p{k)(C), p(m)(c) - биортогональный базис в гильбертовом пространстве N[-1,1], k, m е[0, N] о Z:
(p k),p m ))1 = J-, P( k) (c)p( m) ( Qd C = 0,
k Ф m.
(10)
Здесь и далее подразумевается суммирование по повторяющимся латинским индексам, заключенным в скобки. С учетом (4), (5), (8)—(10) записываются поверхностная и контурная плотности функционала Лагранжа [45, 46]:
д — д ( u( k) u( k) u(k) u(k) V u( k) V u( kЛ =) Pm)(u“ )k) + u3 u(k)) + P“ u( k)-
ДS = ДS yua , u3 , ua , u3 , Vpu« , Vpu3 J= 2 P(k) \u( m)u« + u( m)u3 J + P(k)ui
(k)
(k))+h-is3i)Dm)
и (11)
'(f)(Vpua)- bpa u3*)) + h s(3kP)^(3 -r"y“a " J(k )^(m.)u3
Дг = qB(k)ui + qB(k)ui ^, (12)
P(m)) (M0 ) = (k),Р('И) )j , s(k) = (hs‘J ,P(k) ) , P(k) = F(k)+% + (-k q± .
Обобщенные силы s(jk) являются коэффициентами биортогонального разложения по системе Р(m)(C) компонентов тензора напряжения, заданного в форме (6). Коэффициенты qaB(k) биортогонального разложения компонентов главного вектора внешних сил на кон-
42
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
туре Г = dS0, так же как и вектора приведенных внешних сил на поверхности S0, введены в работах [45, 46]; q± - компоненты главного вектора сил на лицевых поверхностях S± . Линейные операторы D^?) определяются следующими соотношениями [45]:
D,
.’m)=( d. p p(m))
(k •) { d( p(k ),p )l-
(13)
Постановка начально-краевой задачи теории элементарной теории N-го порядка
В соответствии с вариационным формализмом [45-48] из (11) следуют уравнения
движения и естественные краевые условия в следующем виде:
p(mKm,=^$ - ъа$)- ь-'Dtmm*5+%• p(mKm)=^ *(3f)+vs - h-4m;43,+pi, •
(s(f)vp- qу ))Sa!k 1
M 0еГ
= 0; i = 1,2,3; p = 1,2.
Уравнениям движения (14) соответствуют начальные условия
= V/.
\з)
= U(k); u(k)
(k )
(14)
(15)
(16)
где U(k) (M0), V(lk) (M0) - коэффициенты разложения компонентов векторов начального
смещения и скорости оболочки как трехмерного тела по системе P(k)(Q при t = t0.
Определяющие уравнения теории N-го порядка [45, 46] могут быть записаны относи-
(k) (k)
тельно переменных поля uya , u3 и их ковариантных производных:
s(p) = Cjyjm)+Cju(m); s(?)=qjvj+CkmAm).
(k) (km) 5 j (km) J ’ (k) (km) о j (km) у
(17)
Коэффициенты жесткости, входящие в (17), определены в [45, 46] так:
Clру = ъуСР35 + h-1D("n)ClРу3- Clр3 =-ъ ClРу° + h-1D('n)CiP33
(km) 5 (km) (k«) ( nm) (km) у5 (km) (k •) (nm)
Ci3y = ъуCl335 + h-1D('n')Ci3У3; Ci33 =—h Cl3У5 + h-1D('n')Ci333
(km) °5 (km) (k •) (nm) (km) y5 (km) (k •) (nm)'
cijJkl = (km)
(k),p(m))1; c1 Pj5 = д4р;; cipj3 = рА.р'с1у/3; cl3j5 = c13je = vA;criJ\
p = 1 + H + h2 K-2^(h + hK)2K; H = hH, K =h2K, hi = h/h.
Здесь Clm - физические константы теории N-го порядка, определенные в [45]. Постоянные Cijkl, Cljkl образуются при переходе от основного базиса пространственной системы координат ^а, С, к сопутствующему базису ra , n; H, K - безразмерные средняя и гауссова кривизны поверхности S0; h(M0) - отношение отклонения h базисной поверхности
от срединной к толщине оболочки h. Уравнения (14), (17) с краевыми (15) и начальными условиями (16) формируют начально-краевую задачу трехмерной теории оболочек N-го
43
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
порядка [45-48]. Уравнения динамики (14) есть обобщенные уравнения Лагранжа II рода континуальной системы, определенной на двумерном многообразии S0 [45, 46].
2. Формулировка расширенной теории оболочек N-го порядка
Силовые краевые условия на S± е Sc имеют общий вид о|S • n+ = q+, n+ - вектор единичной внешней нормали лицевой поверхности S+, компоненты которого определяются следующим образом:
n±=N±/N + ; N±=R±xR±, R± = R|±x; ^ N ±=p±Van, N ±=p±Va,
^ n± = v+1n; p± = 1 - h±b\+ h+ det (bp1 ).
Главный вектор внешних сил на S+ записывается так: q = (p±v±) (q“ra + q±n), при этом силовые краевые условия на лицевых поверхностях приводятся к виду
s ±
= q±, i = 1,3, p = 1,2.
(18)
s
При учете определяющих уравнений для s(k) (17), полученных в [45, 46] для случая анизотропной упругой среды, краевые условия (18) переходят в уравнения связей для переменных поля I рода u(k), определенные на S0:
СЙ± (<”>) + )D^u'”» -c&=(C(-h+c(5)p(m)(±1); Ck)
ck)±brSu(k)+CibO1 - q±:
+=(cm - h+c('km3)p(” )(±1>
0,
(19)
Учет краевых условий (18) на лицевых поверхностях S+ в форме связей (19) преобразует экстремальную задачу для функционала (11), (12) в экстремальную задачу с ограничениями, решаемую методом множителей Лагранжа [44]. Начально-краевая задачи теории оболочек N-го порядка, следующая из (11), (12) и (19), имеет вид
P(k)u(“m)=VPs(k) bp“s(k) D(k.) s(”)+ P(k),
p(”Km)=v $+baP $ - Dk^w+pk),
(s(k) vp- qB(k) )§u (k)
M 0еГ
= 0; i = 1,2,3; p = 1,2,
si3 -x±CJi3
s(k) KJC(k)±
p(k )(+1)-q± = 0,
s(p) = C;flvsu(,m) + C(p” )u) + rc(ip + XjC/lP ,
(k) (km) 0 j (km) J J (k)+ J (k)-’
sH) = C;lj0vAm) + ck”и{m] eVc3* + ac? ,
(k) (km) 0 j (km) J J (k)+ J (k)-’
\k)
t=t0 = U(k); ^)
= V(k).
(20)
(21)
(22)
(23)
(24)
44
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
Здесь введены соответствующие лицевым поверхностям оболочки S+ и S_ векторы
множителей Лагранжа X1 = X*r“ + X±n, а также следующие обозначения обобщенных физических постоянных:
(25)
3. Уравнения движения плоского изотропного упругого слоя
В случае плоского однородного изотропного упругого слоя с учетом
ba= bap = 0 ^ Д?. = ^ ^ Sj )=5 fk) sj }=&(k},
p(g =p5(m)); C,jpq = C,,pq = Cjpq = Xgijgpq +v(g,pgjq + giqgJp); C;= CIJpqG{km), (26)
paap, i = a, j = p,
Ci)± = C3,3P(k)g±l) = Pgk) (±0 < X + 2(X 1 = 3, j = 3,
0, i = a, j = 3, i ^ j
уравнения связей преобразуются к виду
aa>(k)(±1)_X±paaPpw _q±a = 0; (5g;3>p(k)(±1)_X±(X + 2v)pN _N = 0 (27)
PN = P(k) (±1)P(k) (±1) = G(km)P(") (±1)P(m) (±1) = G(4P(k) (±1)P(m) (±1)- (28)
Для базиса, образованного полиномами Лежандра, pN запишется так:
Pn = Pgk )(±1)P(k )(±1) = S
N 2
(k) 2
Pv ’
k=0
Pgk)(±1) = (±1) , G G diG
(km) =,
g)
(29)
В силу |aap| ^ 0 уравнения связей разрешимы относительно множителей X± :
X3 =(X+2р)_1 pnN[agk3)P(k)(±1)_q±]; X± = р_рЧр[a(ak3)P(k)(±1)_q±
Определяющие уравнения будут иметь следующий вид:
й (aiP=aa,P+XaaPgX J Pgk )(1) + X_ Pg )(_1)),
5 (ai3)=aai3+^aa,(X; Pg, (1) + X_ Pg- )(-1)).
a ;,3л =
*33 = a(k3) +(X+2p) ((W+X3P(k) (_0 ) •
(30)
(31)
Подстановка Xj в определяющие уравнения и разрешение их относительно <5, приводит к уравнениям, связывающим обобщенные усилия a j) и a j):
*j
(k)
aap = 5 (k) =
aap_____^ aapa33
a(k) X+2^ a a(k)
+^r- a*p R(
X ^apD(m) (k)
X+2^
a(3m)_ PN (q+3Pjm)+ q_3P_m
(32)
45
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
P(m \k)
i a3
_a3
a(~k)=°(k)- Rk)P
a3 —2 / —a + , —a -
°(m)— pN P(m)+ q— P(
d33 = )— ^pN (qlp(t
pN(p(k)Plm)+p—k)P—m)); Rkn)p
+q—3p(
))■
— p\m) —S\m)
-°(k) ; P(k) -°(k)
(33)
(34)
(35)
PN - p(k)p± ) > 0 ■
Начально-краевая задача расширенной теории пластин N-го порядка примет вид
р%)-v(t)—Dkm'5(;()+д )■ (36)
)5u,(k) - 0; ( - 1,2,3; р-1,2,
(б (
V)
Mo еГ
t-to - U(k); Uk )
- V k )■
(37)
(38)
4. Спектральная задача для плоского упругого слоя
Рассмотрим задачу о распространении второй моды продольной волны в плоском упругом слое плотностью р, с параметрами Ламе X, ц, толщиной 2h и со срединной плоскостью х2 - 0. Нормальная волна распространяется в направлении оси Ox1 прямоугольной декартовой системы координат Ox1x2х3. На поверхностях слоя х2 - ±h поставлены однородные краевые условия aa2 (х1, ±h) - 0.
Аналогично [53-57] введем безразмерные переменные
£-Xih—1; т-tcxh~1; - u^K1, р-С2/ c - (1 — 2v)/( 2 — 2v), (39)
где c1 -yj(X + 2д)р 1 , c2 -д/др 1 - скорости распространения волн дилатации и сдвига; v - коэффициент Пуассона.
Однородные уравнения движения в безразмерных переменных, соответствующие расширенной теории N-го порядка [45-48], следуют из (14)-(17) и приводятся к матричной форме записи
^u - (C2 ^ + ДД + C0) u, (40)
u
u
(0)
u
(N) u20)
u
( N)
Ck - Ak + Bk, K = 0,1,2.
Операторы D(('km^ (25) и соответствуют базису, образованному полиномами Лежан-
дра p(k )(Q, С- х2 / h е [—1,1] в соответствии с [37-43, 45-49]:
<)= ( 2m + 1)8(m;21+1),1 <[ К N— m — Матрицы линейных операторов D((mk.))
1)]; СИ2т+3)8(mr2M).1 <[ 1—(• (4 1 )
, D((kmj имеют размерность (N +1) х (N +1).
46
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
Вектор перемещения, соответствующий нормальной волне, может быть представлен в виде [53-57]
u(k= U(kk exp [/ [kxl -Qt)] = U(k^ exp [/ (к£-шт)],
- и(k),
(42)
где U(k) - амплитуда; o = Qh/c2 - безразмерная фазовая частота; к-kh - безразмерное волновое число. С учетом (42) уравнения движения расширенной теории N-го порядка (...) приводятся к 2N + 2 линейному уравнению:
(С2к2 + /CjK + C0 - ш21) U - 0,
(43)
U =
u;o)... u(n) и20)... и
(n )■
, C j = A j + B j , J = 1,2;
A 2 =
A1 =
4|32(1 - 2Р2)5("> 0
Р28(Г)' p2 Dk.k;9
f
0 0
-p2 d;
A0 =
Cp2d(ь4(n\ 0^
~ . n .;
.n) (.;
0 0
B 2
4p2 (1 - 2p2) R(
k)
p2 r;
B1
P2 DM< 44 2 (1 - 2P2) 22
2 (1 - 2P2 Id'^-p2 4><l4n:>
CR2 D(k.) rW D(n.)
0
B0
p2D.k№№) 0
0 Dkk; } R((n)) Dk.:)) ,
(44)
(45)
(46)
(47)
, (48)
(49)
Здесь I - единичная; 0 - нулевая матрицы размерности (2N + 2) х(2N + 2) . Операторы AK соответствуют плоскому напряженному состоянию, операторы BK учитывают влияние связей. Из уравнений движения (43) следует спектральная задача [53-57]
(с0 + 'С + к2С2) -о21
= 0.
(50)
Собственные значения оператора (43) on (к)е R+ определяют на плоскости (к,о) дисперсионные кривые; кеЯсоответствуют распространяющимся модам [58]. Собственные векторы U n (к), n e[1,2N + 2 ]n Z задают формы мод:
47
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
< (C) = U<k )n Pk (C), « = 1,2, k = 0...N, n = 1...N +1. (51)
Система уравнений (43) распадается на независимые подсистемы для симметричных и антисимметричных относительно плоскости C = 0 волн, при этом продольным волнам соответствуют индексы
k,m = {2n,N + 2n + 2}; n e [0,[-2(N +1)]] uZ.
5. Анализ описания формы второй распространяющейся продольной моды нормальных волн теориями различного порядка
Рассмотрим далее вторую распространяющуюся моду продольных нормальных волн: ш2(к), U 2( к), ке[0, к*] [48, 49, 57]. На рис. 1 приведены дисперсионные кривые для фазовой частоты ш2(к) и групповой скорости cg2(к), соответствующие второй распространяющейся моде, построенные на основе результатов [48, 49] и точного решения задачи Рэлея-Лэмба [58].
Рис. 1. Фазовая частота ю (а) и групповая скорость cg (б) второй распространяющейся моды продольной волны в упругом изотропном слое. Теории 1, 2... 5-го порядков [48, 49]
Достаточной для практических приложений [59, 60] точности вычисления групповой скорости cg удается достичь при N = 4 . При N < 3 зависимость групповой скорости в
области «обратной волны» качественно отличается от точного решения [58].
Рассмотрим далее эволюцию формы второй распространяющейся моды при к* < 1,4. Согласно (51) соотношения
u2 (к, C) = Ujk)2 (к)рт(С), k = 0,2... N', N'< N (52)
48
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
описывают зависимость продольного перемещения от нормальной координаты С и волнового числа к ; аналогичное соотношение описывает поперечное перемещение и2:
«22 (к,С) = UK (к)р(,,(С), к = 1,3...N', N'< N, (53)
где выбор индексов (к) соответствует базису, образованному полиномами Лежандра при симметрии перемещения в продольной волне: и1 (С) = и1 (-С) ; иП (С) = ~и2 (-С) •
6. Кинематика второй распространяющейся моды
На рис. 2 приведена зависимость и12 (С, к), к е [0,1,5], С е [-1,1], N = 5 .
Рис. 2. Эволюция формы второй распространяющейся моды продольной нормальной волны в слое в области эффекта «обратной волны». Продольное перемещение и12 (к,С) (52). Теория 5-го порядка
На рис. 3 показаны зависимости безразмерного перемещения и2(С, к), определяемого следующим соотношением [57]:
и« ( к) = иу max
^ к)|
(54)
от нормальной координаты С при характерных значениях волнового числа к, соответствующих cg(к) = min и cg(к) = 0, N = 2...6. Кривые на рис. 3 соответствуют точному
решению, полученному в [58]. Как при минимальной отрицательной скорости с"11" , так
и при cg = 0 в сечении слоя существуют подобласти с противоположными знаками амплитуд продольного перемещения:
с = с)
: С,.,0 еИ-°Л0.5]). Cg = 0: С[.,0 е(*[-0,48,MS]).
Распределения по толщине слоя безразмерного продольного перемещения их2(С, к) в точке дисперсионной кривой, соответствующей cgmin , описываемые теориями порядка
49
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
N > 4, практически совпадают [57]. Зависимости u2^,к), соответствующие N = 2,3, качественно отличаются, представляя собой параболическую аппроксимацию перемещения. Аналогичный вывод следует из анализа зависимости u2((^, к) при к : cg = 0. Следует отметить, что в точке cg)m границы подобласти ux2 <0, вычисленные на базе теорий 2...6 порядков, практически совпадают, несмотря на качественное различие кривых. В точке к: cg = 0 границы подобласти u2 < 0, определяемые теориями 2, 3-го и высших порядков, различаются.
Рис. 3. Распределение по толщине слоя безразмерного продольного перемещения:
(a) cg = min; (б) cg = 0. Теории 2, 3, 4, 5-го порядков
На рис. 4 показана эволюция формы поперечного перемещения u^(£, к) во второй распространяющейся моде, при ке [0,1,5], С,е [-1,1], N = 5 .
Рис. 4. Эволюция формы второй распространяющейся моды продольной нормальной волны в слое в области эффекта «обратной волны». Поперечное перемещение u2 (к,(53). Теория 5-го порядка
50
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
На рис. 5 приведены зависимости безразмерного перемещения и2(С,к), определяемого (54), от координаты С при значениях волнового числа к : cg (к) = min и к : cg (к) = 0 . Аналогично предыдущему случаю N = 2... 6.
Рис. 5. Распределение по толщине слоя безразмерного поперечного перемещения:
(а) cg = min; (б) cg = 0. Теории 2, 3, 4, 5-го порядков
Распределение поперечного безразмерного перемещения и^(С, к) по толщине слоя, вычисленное на базе теории второго порядка (см. рис. 5), соответствует линейной аппроксимации и качественно отличается от зависимости и2, полученной на основе теорий 3, 4, 5-го и высших порядков. При этом в точке к : cg^in кривые для всех порядков теории N> 3 совпадают, в точке к: cg = 0 различие между кривыми, соответствующими 3, 4, 5-му порядкам, становится заметным.
По мере повышения порядка кривизна кривых в областях | С |е [0,5,1] увеличивается. При N > 5 кривые совпадают.
На рис. 6, 7 приведены распределения по толщине слоя продольного (а) и поперечного (б) перемещений при дальнейшем увеличении значения волнового числа. Рис. 6 соответствует 2к / ж = 1, рис. 7 - значению 2к/ ж = 1,5 . Как следует из рис. 1, при таких величинах волнового числа знаки фазовой и групповой скоростей совпадают.
Анализируя приведенные на рис. 6, 7 результаты, можно утверждать, что на ближайшем к области проявления эффекта «обратной волны» отрезке к < 1 точность теории
2
4-5-го порядков достаточна для описания поперечного перемещения и2 .
Для достаточной точности аппроксимации продольного перемещения предпочтительна теория 5-го порядка. При 2к / ж * 1,5 погрешность теории 5-го порядка становится заметной.
51
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
б
а
Рис. 6. Распределение по толщине слоя компонентов безразмерного перемещения при значении волнового числа 2к / п = 1. Теории 2, 3, 4, 5, 6, 7-го порядков
;-№=2 ^ 5^| . . 1. . .
II ГЛ и-"- 1 5,6, 7...
\/ N — 4
I -0,5 0 0,5 1
k22(S»
б
Рис. 7. Распределение по толщине слоя компонентов безразмерного перемещения при значении волнового числа 2к/п = 1,5 . Теории 2, 3, 4, 5, 6, 7-го порядков
Заключение
Получена формулировка расширенной теории N-го порядка нетонких анизотропных оболочек, учитывающая силовые краевые условия на лицевых поверхностях, перенесенных на базовую поверхность оболочки, в форме неголономных уравнений связей. Для плоского упругого слоя получена формулировка уравнений движения в переменных поля первого рода и сформулирована соответствующая спектральная задача. Проведенный анализ аппроксимации формы второй распространяющейся моды нормальных волн в упругом слое на базе теории изотропных пластин N-го порядка И.Н. Векуа - А. А. Амосова [31-33, 37-43] в варианте [45-49], при использовании в качестве базисных функций полиномов Лежандра, соответствующих традиционной формулировке «трехмерных» теорий оболочек [19-22, 31-36], позволяет сделать следующие основные выводы.
52
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
1. Для обнаружения эффекта «обратной волны» в плоском упругом слое, т.е. для получения в результате решения соответствующей спектральной задачи на некотором отрезке
к е [0, к*] отрицательных значений групповых скоростей cg < 0 второй распространяющейся продольной моды, минимально необходимой является теория 2-го порядка [48]. В то же время теория 2-го порядка приводит к существенному занижению к *, а также получению
завышенного значения cmin. Продольное перемещение при этом описывается параболической, поперечное - линейной аппроксимацией. Теория 2-го порядка не может быть использована при решении задач, связанных с учетом явления «обратной волны».
2. Теория 3-го порядка за счет дополнительной степени свободы в трансверсальном направлении, обеспечивающей кубическую аппроксимацию поперечного перемещения и, следовательно, лучшее приближение в области «толщинного резонанса» [53] при малых значениях волнового числа, позволяет получить практически точное значение 2к* / л* 0,51 границы области проявления «эффекта обратной волны», однако приводит
к завышению cmin и смещению точки минимума групповой скорости в сторону больших
волновых чисел. Продольное перемещение в рамках теории 3-го порядка описывается параболической аппроксимацией и близко к аппроксимации на базе теории 2-го порядка. Поперечное перемещение в области малых волновых чисел к < к * описывается теорией 3-го порядка с достаточной точностью.
3. Теории 4-го и 5-го порядка обеспечивают не только близкую к точной форму второй дисперсионной ветви [48], но и при малых значениях волнового числа хорошую аппроксимацию формы продольного перемещения, практически неотличимую от аппроксимаций высших порядков при cmin и cg = 0 . Следовательно, теория 4-го порядка может
считаться достаточно точной приближенной моделью для описания поля перемещений в области проявления эффекта «обратной волны». Теория 5-го порядка при волновых числах, не превышающих к*, не предоставляет заметных преимуществ по сравнению с теорией 4-го порядка.
4. При увеличении волнового числа до значений 2к/л * 1, т. е. при описании «прямой» волны в области, непосредственно прилежащей к диапазону волновых чисел, где проявляется эффект «обратной» волны, точность теорий 4-5-го порядков представляется также достаточной для решения практических задач. При дальнейшем увеличении волнового числа теория 5-го порядка не может описать распределение перемещений, соответствующих второй распространяющейся моде, с достаточной точностью: при 2к/л* 1,5 требуется применение теорий 6-7-го порядков для достижения хорошей аппроксимации продольного перемещения .
На основе перечисленных выводов представляется возможным заключить, что для описания явления «обратной волны» в упругом слое теориями пластин типа И.Н. Векуа -А.А. Амосова, основанных на ортогональных разложениях неизвестных по базису, образованному полиномами Лежандра, минимально необходимыми порядками теории являются 4-5-й порядки, обеспечивающие точность построения как дисперсионной ветви, так и распределения перемещений по толщине слоя при распространении второй продольной моды прогрессивных волн.
Работа выполнена при финансовой поддержке РФФИ (проекты № 13-01-00446_а, № 14-01-00488_а, № 14-01-00890_а).
53
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
Библиографический список
1. Kienzler R., Altenbach H., Otts I. Theories of Plates and Shells: Critical Review and New Applications. - Berlin: Springer-Verlag, 2004.
2. Eremeyev V., Pietraszkiewicz W. Refined theories of plates and shells // Zeitschrift fur Angewandte Mathematik und Mechanik. - 2014. - Vol. 94(1-2). - P. 5-6.
3. Галиньш А. К. Расчет пластин и оболочек по уточненным теориям // Исследования по теории пластин и оболочек. - Вып. 6-7. - Казань: Изд-во Казан. ун-та, 1970. - С. 23-64.
4. Григолюк Э.И., Селезов И.Т. Итоги науки и техники. МДТТ. Т.5. Неклассические теории колебаний стержней, пластин и оболочек. - М.: ВИНИТИ, 1973. - 273 с.
5. Carrera E. A study of transverse normal stress effect on vibration of multilayered plates and shells // J. Sound & Vibration. - 1999. - Vol. 225(5). - P. 803-829.
6. Carrera E. A Reissner’s mixed variational theorem applied to vibration analysis of multilayered shell // Trans. ASME. J. Appl. Mech. - 1999. - Vol. 66. - P. 69-78.
7. Carrera E. Theories and Finite Elements for Multilayered, Anisotropic, Composite Plates and Shells // Arch. Comput. Methods Eng. - 2002. - Vol. 9(2). - P. 87-140.
8. Koiter W.T. A consistent first approximations in the general theory of thin elastic shells // Proc. of Symp. on the Theory of Thin Elastic Shells. August 1959, North Holland, Amsterdam. - P. 12-23.
9. Carrera E., Brischetto S. Analysis of thickness locking in classical, refined and mixed multilayered plate theories // Composite Structures. - 2008. - Vol. 82. - P. 549-562.
10. Brichetto S, Carrera S. Importance of Higher Order Modes and Refined Theories in in Free Vibration Analysis of Composite Plates // J. Appl. Math. - 2010. - Vol. 77(1), 011013 (14).
11. Batra R.C., Aimmanee S. Vibrations of thick isotropic plates with higher order shear and normal deformable plate theories // Comput. and Struct. - 2005. - Vol. 83. - P. 934-955.
12. Three-dimensional vibration analysis of circular and annular plates via the Chebyshev-Ritz method /
D. Zhou, F.T.K. Au, Y.K. Cheung, S.H. Lo // Int. J. Sol. Struct. - 2003. - Vol. 40. - P. 3089-3105.
13. Kang Jae-Hoon, Leissa A. W. Three-Dimensional Vibration Analysis of Thick, Complete Conical Shells // Trans. of ASME. - 2004. - Vol. 71(6). - P. 502-507.
14. Matsunaga H. Stress analysis of functionally graded plates subjected to thermal and mechanical loadings // Compos. Struct. - 2009. - Vol. 87. - P. 344-357.
15. Matsunaga H. Free vibration and stability of functionally graded plates according to a 2D higher-order deformation theory // Compos. Struct. - 2008. - Vol. 82. - P. 499-512.
16. Matsunaga H. Thermal buckling of functionally graded plates according to a 2D higher-order deformation theory // Compos. Struct. - 2009. - Vol. 90. - P. 76-86.
17. Matsunaga H. Free vibration and stability of functionally graded shallow shells according to a 2D higher-order deformation theory // Compos. Struct. - 2008. - Vol. 84. - P. 132-146.
18. Matsunaga H. Free vibration and stability of functionally graded circular cylindrical shells according to a 2D higher-order deformation theory // Compos.Struct. - 2009. - Vol. 88. - P. 519-531.
19. Абросимов Н.А., Баженов В.Г., Куликова Н.А. Идентификация вязкоупругих характеристик композитных материалов по результатам экспериментально- теоретического анализа динамического поведения полусферических оболочек // ПМТФ. - 2006. - Т. 47, № 3. - С. 126-133.
20. Абросимов Н.А., Куликова Н.А. Расчетно-экспериментальный метод идентификации вязкоупругих характеристик композитных материалов в динамически нагруженных оболочках вращения // Механика композит. материалов. - 2007. - Т.43, № 4. - С. 449-464.
21. Абросимов Н.А., Куликова Н.А. Определение параметров моделей нелинейного деформирования изотропных и композитных материалов по результатам расчетно- экспериментального анализа импульсного нагружения круглых пластин // ПМТФ. - 2011. - Т. 52, № 1. - С. 163-172.
22. Абросимов Н.А., Куликова Н.А. Идентификация параметров моделей вязкоупругого деформирования композитных материалов по результатам расчетно-экспериментального анализа импульсного нагружения оболочек вращения // Изв. РАН. МТТ. - 2011. - № 3. - С. 42-57.
54
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
23. Ардзашвили Р.В., Вильде М.В., Коссович Л.Ю. Антисимметричные кромочные волны высшего порядка в пластинах // Изв. Саратов. ун-та. Нов. сер. Математика, механика, информатика. - 2013. - Вып. 1. Ч. 1. - Т. 13. - С. 51-56.
24. Кильчевский Н.А. Основы аналитической механики оболочек. - Киев: Изд-во АН УССР, 1963. - 634 с.
25. Ворович И.И. Некоторые математические задачи теории пластин и оболочек // Тр. II Все-союзн. съезда по теор. и прикл. мех., 1964. - М.: Наука, 1966. - С. 116-136.
26. Гольденвейзер А.Л. Построение приближенной теории оболочек при помощи асимптотического интегрирования уравнений теории упругости // ПММ. - 1963. - Т. 27, № 4. - С. 593-608.
27. Коссович Л.Ю. Асимптотические методы в динамике оболочек при ударных воздействиях // Изв. Сарат. ун-та. Сер. Математика, механика, информатика. - 2008. - Т. 8. - Вып. 2. - С. 12-33.
28. Zhilin P.A. Mechanics of Deformable Cosserat Surfaces // Int. J. Sol. Struct. - 1976. - Vol. 12. -P.635-648.
29. Жилин П.А. Прикладная механика: основы теории оболочек: учеб. пособие. - СПб.: Изд-во С.-Петерб. политехн. ун-та, 2006. - 167 с.
30. Yokoo Y., Matsunaga H. A general nonlinear theory of elastic shells // Int. J. Sol. Struct. -1974. - Vol. 10. - P. 261-274.
31. Векуа И.Н. Некоторые общие методы построения различных вариантов теории оболочек. -М: Наука, 1982. - 282 с.
32. Гуляев В.И., Баженов В.А., Лизунов П.П. Неклассическая теория оболочек и ее приложение к решению инженерных задач. - Львов: Вища Школа, 1978. - 192 с.
33. Хома И.Ю. Общая теория анизотропных оболочек. - Киев: Наук. думка,1986. - 170 с.
34. Абросимов Н.А., Баженов В.Г. Исследование динамического деформирования упругопластических сферических оболочек при тепловом ударе // Изв. АН СССР. МТТ. - 1978. - № 1. -С. 139-143.
35. Абросимов Н.А., Баженов В.Г. Нелинейные задачи динамики композитных конструкций. -Н. Новгород: Изд-во Нижегород. гос. ун-та, 2002. - 400 с.
36. Никабадзе М.У. Развитие метода ортогональных полиномов в механике микрополярных и классических упругих тонких тел. - М.: Изд-во МГУ, 2014. - 515 с.
37. Амосов А.А. Приближенная трехмерная теория толстостенных пластин и оболочек // Строительная механика и расчет сооружений. - 1987. - № 5. - С. 37-42.
38. Амосов А.А., Жаворонок С.И. К проблеме редукции плоской задачи теории упругости к последовательности одномерных краевых задач // Мех. композ. матер. и констр. - 1997. - Т. 3, № 1. - С. 69-80.
39. Амосов А.А., Князев А.А., Жаворонок С.И. О решении некоторых краевых задач о плоском напряженном состоянии криволинейной трапеции // Механика композ. матер. и констр. -1999. - Т. 5, № 1. - С. 60-72.
40. Amosov A.A., Zhavoronok S.I. An approximate high-order theory of thick anisotropic shells // Int. J. Comput. Civil Struct. Eng. - 2003. - Vol. 1. - P. 28-38.
41. Амосов А. А., Жаворонок С.И., Леонтьев К. А. О решении некоторых задач о напряженнодеформированном состоянии анизотропных толстостенных оболочек вращения в трехмерной постановке // Мех. композ. матер. и констр. - 2004. - T. 10, № 3. - С. 301-310.
42. Жаворонок С. И. Модели высшего порядка анизотропных оболочек // Механика композ. матер. и констр. - 2008. - T. 14, № 4. - С. 561-571.
43. Жаворонок С.И., Леонтьев А.Н., Леонтьев К. А. Анализ сходимости решения при расчете толстостенных оболочек вращения произвольной формы // Int. J. Comput. Civil Struct. Eng. -2010. - Vol. 6(1-2). - P. 105-111.
44. Кильчевский Н.А., Кильчинская Г.А., Ткаченко Н.Е. Аналитическая механика континуальных систем. - Киев: Наукова Думка, 1979. - 189 с.
55
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
45. Жаворонок С.И. Вариационные уравнения трехмерной теории анизотропных оболочек // Вестн. Нижегород. ун-та им. Н.И. Лобачевского. - 2011. - № 4(5). - С. 2153-2155.
46. Жаворонок С.И. Обобщенные уравнения Лагранжа второго рода трехмерной теории анизотропных оболочек // Мех. композ. матер. и констр. - 2011. - Т. 17, № 1. - С. 116-132.
47. Zhavoronok S.I. A Vekua-type linear theory of thick elastic shells // Zeitschrift fur Angewandte Mathematik und Mechanik. - 2014. - Vol. 94 (1-2). - P. 164-184.
48. Zhavoronok S.I. Variational formulations of Vekua-type shell theories and some their applications // Shell Structures: Theory and Applications. Vol. 3. CRC Press. - Balkema: Taylor & Francis Gr., 2014. - P. 341-344.
49. Жаворонок С.И. Формулировка начально-краевой задачи приближенной трехмерной теории оболочек N-го порядка в обобщенных перемещениях и ее приложение к задачам стационарной динамики // Мех. композ. матер. и констр. -2012. - Т. 18, № 3. - С. 333-344.
50. Пикуль В.В. Физически корректные модели материала упругих оболочек // Изв. РАН. МТТ. - 1995. - № 2. - С. 103-108.
51. Schwab C., Wright S. Boundary layers in hierarchical beam and plate models // J. Elast. -1995. - Vol. 38. - С. 1-40.
52. Kulikov G.M, Plotnikova S.V. Equivalent single-layer and layerwise shell theories and rigid-body motions. Part I: Foundations // Mech. Adv. Mater. Struct. - 2005. - Vol. 12(4). - P. 275-283.
53. Жаворонок С.И. Исследование гармонических волн в упругом слое на основе трехмерной теории оболочек N-го порядка // Механика композиционных материалов и конструкций. - 2010. -T. 16, № 4/2. - С. 693-701.
54. Жаворонок С.И. Трехмерные модели динамического деформирования толстостенных оболочек // Морские интел. технол. -2011. - № 3 (спец.). - C. 6-10.
55. Жаворонок С.И. Исследование распространяющихся мод гармонических волн в упругом слое на базе трехмерной теории оболочек N-го порядка // Мех. композ. матер. и констр. - 2011. -T. 17, № 2. - С. 278-287.
56. Жаворонок С.И. Исследование кинематики нормальных волн в упругом слое на основе трехмерной теории оболочек N-го порядка для различных значений волновых чисел // Механика композ. матер. и констр. - 2012. - T. 18, № 1. - С. 45-56.
57. Егорова О.В., Жаворонок С.И., Курбатов А.С. О приложении различных вариантов теории оболочек N-го порядка к некоторым задачам о прогрессивных волнах // Изв. Тул. гос. ун-та. Сер. Технические науки. - 2014. - № 11. - Ч. 1. - С. 255-266.
58. Гринченко В.Т., Мелешко В.В. Гармонические колебания и волны в упругих телах. -Киев: Наукова думка, 1981. - 284 с.
59. Meitzler A. H. Backward-wave transmission of stress pulses in elastic cylinders and plates // J. Acoust. Soc. Amer. - 1965. - Vol. 38(5). - P. 835-842.
60. Артоболевский И.И., Бобровницкий Ю.И., Генкин М.Д. Введение в акустическую динамику машин. - М: Наука, 1979. - 306 с.
References
1. Kienzler R., Altenbach H., Otts I. Theories of Plates and Shells: Critical Review and New Applications. Berlin: Springer-Verlag, 2004.
2. Eremeyev V., Pietraszkiewicz W. Refined theories of plates and shells. Zeitschrift fur Angewandte Mathematik und Mechanik, 2014, vol. 94(1-2), pp. 5-6.
3. Galinsh A.K. Raschet plastin i obolochek po utochnennym theoriyam (Computation of plates and shells using advanced theories). Issledovaniia po teorii plastin i obolochek. Kazanskii universitet, 1970, iss. 6-7, pp. 23-64.
4. Grigoliuk E.I., Selezov I.T. Itogi nauki i tekhniki. Mekhanika deformiruemogo tverdogo tela. Vol. 5. Neklassicheskie teorii kolebanii sterzhnei, plastin i obolochek [Resumes of Science and Technics.
56
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
Mechanics of Solids. Vol. 5. Non-classical theories of beams, plates and shells vibrations]. Moscow: VINITI, 1973. 273 p.
5. Carrera E. A study of transverse normal stress effect on vibration of multilayered plates and shells. J. Sound & Vibration, 1999, vol. 225(5), pp. 803-829.
6. Carrera E. A Reissner’s mixed variational theorem applied to vibration analysis of multilayered shell. Trans. ASME. J. Appl. Mech, 1999, vol. 66, pp. 69-78.
7. Carrera E. Theories and Finite Elements for Multilayered, Anisotropic, Composite Plates and Shells. Arch. Comput. Methods Eng., 2002, vol. 9(2), pp. 87-140.
8. Koiter W.T. A consistent first approximations in the general theory of thin elastic shells. Proc. of Symp. on the Theory of Thin Elastic Shells. August 1959. North Holland, Amsterdam, pp. 12-23.
9. Carrera E., Brischetto S. Analysis of thickness locking in classical, refined and mixed multilayered plate theories. Composite Structures, 2008, vol. 82, pp. 549-562.
10. Brichetto S., Carrera S. Importance of Higher Order Modes and Refined Theories in in Free Vibration Analysis of Composite Plates. J. Appl. Math, 2010, vol. 77(1), 011013 (14 pp).
11. Batra R.C., Aimmanee S. Vibrations of thick isotropic plates with higher order shear and normal deformable plate theories. Comput. and Struct, 2005, vol. 83, pp. 934-955.
12. Zhou D., Au F.T.K., Cheung Y.K., Lo S.H. Three-dimensional vibration analysis of circular and annular plates via the Chebyshev-Ritz method. Int. J. Sol. Struct, 2003, vol. 40, pp. 3089-3105.
13. Kang Jae-Hoon, Leissa A.W. Three-Dimensional Vibration Analysis of Thick, Complete Conical Shells. Trans. of ASME, 2004, vol. 71(6), pp. 502-507.
14. Matsunaga H. Stress analysis of functionally graded plates subjected to thermal and mechanical loadings. Compos. Struct, 2009, vol. 87, pp. 344-357.
15. Matsunaga H. Free vibration and stability of functionally graded plates according to a 2D higher-order deformation theory. Compos. Struct, 2008, vol. 82, pp. 499-512.
16. Matsunaga H. Thermal buckling of functionally graded plates according to a 2D higher-order deformation theory. Compos. Struct, 2009, vol. 90, pp. 76-86.
17. Matsunaga H. Free vibration and stability of functionally graded shallow shells according to a 2D higher-order deformation theory. Compos. Struct., 2008, vol. 84, pp. 132-146.
18. Matsunaga H. Free vibration and stability of functionally graded circular cylindrical shells according to a 2D higher-order deformation theory. Compos. Struct., 2009, vol. 88, pp. 519-531.
19. Abrosimov N.A., Bazhenov V.G., Kulikova N.A. Identifikatsiia viazkouprugikh kharakteristik kompozitnykh materialov po rezul'tatam eksperimental'no-teoreticheskogo analiza dinamicheskogo povedeniia polusfericheskikh obolochek [Identification of viscoelastic characteristics of composite materials on the basis of results of an experimental and theoretical analysis of the dynamic behavior of hemispherical shells]. J. Appl. Mech. Techn. Phys, 2006, vol. 47(3), pp. 412-418.
20. Abrosimov N.A., Kulikova N.A. Raschetno-eksperimental'nyi metod identifikatsii viazkouprugikh kharakteristik kompozitnykh materialov v dinamicheski nagruzhennykh obolochkakh vrashcheniia [A computational-experimental method for identifying the viscoelastic properties of composites in dynamically loaded shells of revolution]. Mechanics of Compos. Mater., 2007, vol. 43(4), pp. 449-464.
21. Abrosimov N.A., Kulikova N.A. Opredelenie parametrov modelei nelineinogo deformirovaniia izotropnykh i kompozitnykh materialov po rezul'tatam raschetno-eksperimental'nogo analiza impul'snogo nagruzheniia kruglykh plastin [Determining the model parameters of nonlinear deformation of isotropic and composite materials from the results of an experimental-computational analysis of impulsively loaded circular plates]. J. Appl. Mech. Techn. Phys, 2011, vol. 52(1), pp. 132-140
22. Abrosimov N.A., Kulikova N.A. Parameter identification in viscoelastic strain models for composite materials by analyzing impulsive loading of shells of revolution. Mechanics of Solids: A Journal of Russian Academy of Sciences, 2011, vol. 46(3), pp. 368-379.
23. Ardzashvili R.V., Wilde M.V., Kossovich L.Yu. Antisimmetrichnye kromochnye volny vyschego poriadka v plastinakh [High-order boundary waves in plates]. Izvestiia Saratovskogo universiteta. Mathematika, Mechanika, Informatika, 2013, vol. 13, iss. 1, part 1, pp. 51-56.
57
Egorova O.V., Kurbatov A.S., Zhavoronok S.I. / PNRPU Mechanics Bulletin 2 (2015) 36-59
24. Kilchevskiy N.A. Osnovy analiticheskoi mekhaniki obolochek [Basics of analytical shell mechanics]. Kiev: Ukrainian Acad. Sci., 1963. 634 p.
25. Vorovich I.I. Nekotorye matematicheskie zadachi teorii plastin i obolochek [Some mathematical problems of plate and shell theory]. Proc. of II All-Union Conf. on Theoretical and Applied Mechanics, 1964. Review. Moscow: Nauka, 1966, pp. 116-136.
26. Goldenweiser A.L. Postroenie priblizhennoi teorii obolochek pri pomoshchi asimptoticheskogo integrirovaniia uravnenii teorii uprugosti [Derivation of an approximate theory of shells by means of asymptotic integration of the equations of the theory of elasticity]. Appl. Math. Mech, 1963, vol. 27(4), pp. 903-924.
27. Kosovich L.Yu. Asimptoticheskie metody v dinamike obolochek pri udarnykh vozdeistviakh [Asymptotic Methods in Dynamics of Shells under Shock Loading]. Izvestiia Saratovskogo universiteta. Mathematika, Mechanika, Informatika, 2008, vol. 8, iss. 2, pp. 12-33.
28. Zhilin P.A. Mechanics of Deformable Cosserat Surfaces. Int. J. Sol. Struct., 1976, vol. 12, pp.635-648.
29. Zhilin P.A. Prikladnaia mekhanika: osnovy teorii obolochek [Applied Mechanics: Basics of Shells’ Theory]. St. Petersburgskii polytechnicheskii universitet, 2006. 167 p.
30. Yokoo Y., Matsunaga H. A general nonlinear theory of elastic shells. Int. J. Sol. Struct, 1974, vol. 10, pp. 261-274.
31. Vekua I.N. Nekotorye obshchie metody postroenia razlichnykh variantov teorii obolochek [Some general methods of construction of various variants of theory of shells]. Moscow: Nauka, 1982. 282 p.
32. Guliaev V.I., Bazhenov V.A., Lizunov P.P. Neklassicheskaia teoria obolochek i ee prilozhenie k resheniiu inzhenernykh zadakh [Non-classical shell theory and its application to solving of enrineering problems]. Lwov: Vischa Shkola, 1978. 192 p.
33. Khoma I.Yu. Obschaia teoriia anizotropnykh obolochek [General Theory of Anizotropic Shells]. Kiev: Naukova Dumka, 1986. 170 p.
34. Abrosimov N.A., Bazhenov V.G. Issledovanie dinamicheskogo deformirovaniia uprugoplasticheskikh sfericheskikh obolochek pri teplovom udare [Investigation of the dynamic strain state of elastic-plastic spherical shells under heat impact]. Mechanics of Solids: A Journal of Russian Academy of Sciences, 1978, no. 1, pp. 139-143.
35. Abrosimov N.A., Bazhenov V.G. Nelineinye zadachi dinamiki kompositnykh konstrukcii [Nonlinear Problems of Dynamics of Composite Structures]. Nizhegorodskii universitet, 2002. 400 p.
36. Nkabadze M.U. Razvitie metoda ortogonalnykh polinomov v mekhanike mikropolarnykh i klassicheskikh uprugikh tonkikh tel [Development of the Orthogonal Polynomials Method in the Mechanics of Micropolar and Traditional Elastic Bodies]. Moscovskii gosudarstvennyi universitet, 2014. 515 p.
37. Amosov A.A. Priblizhennaia trekhmernaia teoriia tolstostennykh plastin i obolochek [An approximated three-dimensional theory of thick-walled shells and plates]. Stroitelnaia mekhanika i raschet sooruzhenii, 1987, no. 5, pp. 37-42.
38. Amosov A.A., Zhavoronok S.I. K probleme reduktsii ploskoi zadachi teorii uprugosti k posledovatel'nosti odnomernykh kraevykh zadach [Reduction of the plane problem of elasticity theory to a sequence of one-dimensional boundary-value problems]. Composite Mechanics and Design, 1997, vol. 3(1), pp. 69-80.
39. Amosov A.A., Kniazev A.A., Zhavoronok S.I. On solution of 2D-problem of stressed curvilinear trapezoid. Composite Mechanics and Design, 1999, vol. 5(1), pp. 60-72.
40. Amosov A.A., Zhavoronok S.I. An approximate high-order theory of thick anisotropic shells. Int. J. Comput. Civil Struct. Eng, 2003, vol. 1, pp. 28-38.
41. Amosov A.A., Zhavoronok S.I., Leontiev K.A. O reshenii nekotorykh zadach o napriazhenno-deformirovannom sostoianii anizotropnykh tolstostennykh obolochek vrashcheniia v trekhmernoi postanovke [On the solution of some problems of stress-strain state of thick anisotropic shells of revolution on the basis of the three-dimensional problem’s statement]. Compos. Mech. and Design, 2004, vol. 10(3), pp. 301-310.
58
Егорова О.В., Жаворонок С.И., Курбатов А. С. /Вестник ПНИПУ. Механика 2 (2015) 36-59
42. Zhavoronok S.I. Modeli vysshego poriadka anizotropnykh obolochek [High-order anisotropic shells models]. Composite Mechanics and Design, 2008, vol. 14(4), pp. 561-571.
43. Zhavoronok S.I., Leontiev A.N., Leontiev K.A. Analysis of thick-walled rotation shells based on Legendre polynomials. Int. J. Comput. Civil Struct. Eng., 2010, vol. 6(1-2), pp. 105-111.
44. Kilchevskiy N.A., Kilchinskaia G.A., Tkachenko N.E. Analiticheskaia mekhanika konti-nualnykh system [Analytical mechanics of continual systems]. Kiev: Naukova Dumka, 1979. 189 p.
45. Zhavoronok S.I. Variacionnye uravnenia trekhmernoi teorii anizotropnykh obolochek [Variational equations of the three-dimensional theory of anisotropic shells]. Vestnik Nizhegorodskogo universiteta imeni N.I. Lobachevskogo, 2011, vol. 4(5), pp. 2153-2155.
46. Zhavoronok S.I. Obobshchennye uravneniia Lagranzha vtorogo roda trekhmernoi teorii anizotropnykh obolochek [Generalized Lagrange equations of the second kind of three-dimensional anisotropic shell theory]. Composite Mechanics and Design, 2011, vol. 17(1), pp. 116-132.
47. Zhavoronok S.I. A Vekua-type linear theory of thick elastic shells. Zeitschrift fur Angewandte Mathematik undMechanik, 2014, vol. 94 (1-2), pp. 164-184.
48. Zhavoronok S.I. Variational formulations of Vekua-type shell theories and some their applications. Shell Structures: Theory and Applications, vol. 3. CRC Press. Balkema: Taylor & Francis Gr., 2014. 341-344 p.
49. Zhavoronok S.I. Formulirovka nachal'no-kraevoi zadachi priblizhennoi trekhmernoi teorii obolochek Ngo poriadka v obobshchennykh peremeshcheniiakh i ee prilozhenie k zadacham statsionarnoi dinamiki [A formulation of the three-dimensional approximated shells theory of n-th order using generalized displacements and its application to steady dynamics]. Composite Mechanics and Design, 2012, vol. 18(3), pp. 333-344.
50. Pikul V.V. Fizicheski korrektnye modeli materiala uprugikh obolochek [Physically correct material models for elastic shells]. Mechanics of Solids: A Journal of Russian Academy of Sciences, 1995, no. 2, pp.103-108.
51. Schwab C., Wright S. Boundary layers in hierarchical beam and plate models. J. Elast, 1995, vol. 38, pp. 1-40.
52. Kulikov G.M, Plotnikova S.V. Equivalent single-layer and layerwise shell theories and rigid-body motions. Part I: Foundations. Mech. Adv. Mater. Struct, 2005, vol. 12(4), pp. 275-283.
53. Zhavoronok S.I. Issledovanie garmonicheskikh voln v uprugom sloe na osnove trekhmernoi teorii obolochek N-go poriadka [Investigation of harmonic waves in elastic layer using N-th order threedimensional shells theory]. Composite Mechanics and Design, 2010, vol. 16(4/2), pp. 693-701.
54. Zhavoronok S.I. Trekhmernye modeli dinamicheskogo deformirovania tolstostennykh obolochek [Three-dimensional models of dynamical deformed state of thick-walled shells]. Morskie intellektualnye tekhnologii, 2011, vol. 3, pp. 6-10.
55. Zhavoronok S.I. Issledovanie rasprostraniaiushchikhsia mod garmonicheskikh voln v uprugom sloe na baze trekhmernoi teorii obolochek N-go poriadka [Investigation of propagating modes of harmonic waves in elastic layer using N-th order three-dimensional shells theory]. Comp.Mech.and Design, 2011, vol. 17(2), pp. 278-287.
56. Zhavoronok S.I. Kinematics of normal modes in elastic layer for some wavenumbers investigation based on N-th order three-dimensional shells’ theory. Composite Mechanics and Design, 2012, vol. 18(1), pp. 45-56.
57. Egorova O.V., Zhavoronok S.I., Kurbatov A.S. O prilozhenii razlitchnykh variantov teorii obolochek N-go poriadka k nekotorym zadacham o progrssivnykh volnakh [On the application of different variants of N’th order shell theory to some problems of normal waves propagation]. Part 1. Izvestiia Tulskogo gosudarstvennogo universiteta. Technical Sci., 2014, no. 11, pp. 255-266.
58. Grinchenko V.T., Meleshko V.V. Garmonicheskie kolebaniia i volny v uprugikh telakh [Harmonic waves in elastic bodies]. Kiev: Naukova dumka, 1981. 284 pp.
59. Meitzler A.H. Backward-wave transmission of stress pulses in elastic cylinders and plates. J. Acoust. Soc. Amer., 1965, vol. 38(5), pp. 835-842.
60. Artobolevskiy I.I., Bobrovnitskiy Yu.I., Genkin M.D. Vvedenie v akusticheskuiu dinamiku mashin [Introduction into the acoustical machine dynamics]. Moscow: Nauka, 1979. 306 p.
59
