CC BY
30
УДК 539.3
И. Г. Емельянов 12, А. В. Кузнецов 1
1 Институт машиноведения Уральского отделения РАН (Екатеринбург, Россия)
2 Уральский государственный университет путей сообщения (Екатеринбург, Россия)
ОПРЕДЕЛЕНИЕ НАПРЯЖЕННОГО СОСТОЯНИЯ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ С ИСПОЛЬЗОВАНИЕМ МЕТОДОВ ТЕОРИИ ОБОЛОЧЕК
Дата поступления 22.09.2017 Решение о публикации 17.10.2017
Аннотация:
Введение. С развитием численных методов и вычислительных комплексов произвести оценку напряженного состояния тонкостенных конструкций, в виде тел вращения достаточно просто. Однако при решении таких задач методом конечных элементов необходимо выбрать такую сетку конечных элементов, чтобы «схватить» все возможные сингулярности напряженного состояния. Чтобы корректно их учитывать, необходимо уменьшать размер конечных элементов. Уменьшение размеров элементов приводит к увеличению требуемых вычислительных мощностей.
Постановка проблемы. При решении прикладных задач даже при достаточно грубой сетке количество элементов может превышать сотни тысяч. При решении задач для реальных конструкций в трехмерной постановке, объем вычислений может быть достаточно велик и далеко не всякий даже суперкомпьютер сможет обработать такое решение.
Цель работы. Цель данной работы использовать известный подход, применяемый в теории оболочек, который позволяет свести трехмерную задачу к решению одномерной задачи, что существенно снижает требования к вычислительным мощностям.
Метод (методология). Рассматривается задача об определении напряженного состояния оболочечных конструкций в виде тел вращения. Подход основан на интегрировании уравнений теории оболочек и разложении функций в ряды Фурье для разделения переменных. В работе использовано разложение в дискретный ряд Фурье по косинусам и синусам, которое описывает произвольные несимметричные механические нагрузки.
Результаты. Рассмотрена тонкостенная цилиндрическая конструкция, шарнирно закрепленная по торцам. Конструкция нагружена в трех местах распределенной силой, действующей по нормали к поверхности оболочки. После интегрирования системы уравнений для оболочки найденное напряженно-деформированное состояние оболочки определяется компонентами напряжений на внешней и внутренней поверхностях оболочки и компонентами перемещений. В работе приводится сравнение результатов расчета с помощью предлагаемой методики и метода конечных элементов.
Заключение. Показано, что использование методов теории оболочек, и предложенное разложение разрешающих функций и нагрузки в ряд Фурье, позволяет решать задачи с использованием небольших вычислительных ресурсов. При этом
обеспечивается необходимая точность вычисления по всем компонентам напряженно-деформированного состояния конструкции.
Ключевые слова: напряженное состояние, оболочка вращения, метод дискретной ортогонализации, ряд Фурье, метод конечных элементов.
UDC 539.3
L. G. Emel'yanov 12, A. V. Kuznetsov 1
1 Federal State Budgetary Scientific Institution "Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences", (Ekaterinburg, Russia)
2 Ural State University of Railway Transport, (Ekaterinburg, Russia) DETERMINATION OF THE STRESSED STATE OF THIN-CONSTRUCTION STRUCTURES USING THE METHODS OF THE THEORY OF SHELLS
Annotation:
Introduction. With the development of numerical methods and computational complexes, it is quite easy to evaluate the stress state of thin-walled structures in the form of rotation bodies. However, when solving such problems by the finite element method, it is necessary to choose such finite element grid to "grasp" all possible singularities of the stressed state. To correctly take them into account, you must reduce the size of the finite elements. Reducing the size of the elements leads to an increase in the required computing power.
Formulation of the problem. When solving applied problems, even with a sufficiently coarse grid, the number of elements can exceed hundreds of thousands. When solving problems for real constructions in a three-dimensional setting, the amount of computation can be quite large and not every supercomputer can even handle such a solution.
Objective. The purpose of this paper is to use the well-known approach used in shell theory, which allows us to reduce the three-dimensional problem to the solution of a one-dimensional problem, which substantially reduces the requirements for computing power.
Method (methodology). The problem of determining the stress state of shell structures in the form of bodies of revolution is considered. The approach is based on the integration of the equations of the theory of shells and the expansion of functions into Fourier series for separation of variables. The expansion into a discrete Fourier series in cosines and sines is used in this paper, which describes arbitrary asymmetric mechanical loads.
Results. A thin-walled cylindrical structure hinged at the ends is considered. The structure is loaded in three places by a distributed force acting normal to the surface of the shell. After integrating the system of equations for the shell, the found stress-strain state of the shell is determined by the stress components on the outer and inner surfaces of the shell and the displacement components. The paper compares the calculation results with the proposed methodology and the finite element method.
The conclusion. It is shown that the use of methods of shell theory, and the proposed expansion of resolving functions and loads in a Fourier series, allows solving problems using small computing resources. At the same time, the necessary accuracy of calculation for all components of the stress-strain state of the structure is ensured.
Key words: stress state, shell of rotation, discrete orthogonalization method, Fourier series, finite element method.
Введение
В настоящее время при определении напряженного состояния различных конструкционных элементов широко используются вычислительные комплексы, основанные на методе конечных элементов. Произвести оценку напряженного состояния тонкостенных конструкций, в виде тел вращения используя такие комплексы как ANSYS, COSMOS [1, 2] и др. не представляет труда даже в геометрически и физически нелинейной постановке. Напряженное состояние может быть найдено в одномерной, двумерной или трехмерной постановке в зависимости от типа и цели задачи и располагаемой вычислительной мощности. Однако существует класс задач, где заложенная в данных вычислительных комплексах универсальность, теряет свою эффективность. Это происходит при решении задач с подвижными границами, т.е. контактных задач для конструкций, взаимодействующих с жесткими или упругими телами. При решении задач методом конечных элементов для тонкостенных конструкций, взаимодействующих с другими телами, необходимо выбрать такую сетку конечных элементов, чтобы «схватить» все возможные сингулярности напряженного состояния, но при этом использовать минимально возможную размерность задачи. Этого не всегда легко добиться, поскольку, например, при контакте оболочечных конструкций возникают локализованные нагрузки на границе области контакта [3, 4, 5]. Чтобы корректно их учитывать, необходимо уменьшать размер конечных элементов. Уменьшение размеров элементов приводит к увеличению требуемых вычислительных мощностей. При решении реальных прикладных задач даже при достаточно грубой сетке количество элементов может превышать сотни тысяч. Например, при определении не осесимметричного напряженного состояния (двумерная задача) оболочки вагона-цистерны потребовалось около 100 тыс. оболочечных конечных элементов типа shell [6]. При решении такой задачи в одномерной постановке (осесимметричное напряженное состояние) необходимо около 1 тыс. элементов. При решении данной задачи в трехмерной постановке (в рамках теории упругости) необходимо не менее 2 х106 объемных конечных элементов типа solid для случая если толщина оболочки будет аппроксимирована одним конечным элементом (при необходимом соблюдении соотношений между размерами конечного элемента). Если для увеличения точности расчета толщина оболочки будет аппроксимирована тремя элементами, то число элементов возрастет до 1х1°7 конечных элементов. Если толщина оболочки будет аппроксимирована шестью элементами необходимо не менее 8 х1°7 элементов. Если для первого расчетного случая для получения решения необходимо около суток, то для третьего случая речь пойдет уже о месяцах, к тому же объем вычислений
будет достаточно велик и далеко не всякий компьютер сможет такое решение обработать. Эти оценки приведены для решения задачи для оболочки вагона-цистерны на суперкомпьютере «Уран» ИММ УрО РАН. Следовательно, решение такой задачи в трехмерной постановке становится затруднительным даже при использовании суперкомпьютеров. Поэтому вопрос об уменьшении размерности задачи для решения прикладных задач для тонкостенных конструкций является весьма актуальным.
Постановка задачи
Одним из эффективных методов понижения размерности задачи для тонкостенных конструкций в виде тел вращения является использования методов теории оболочек и применение рядов Фурье. Известно, что в методе конечных элементов для тел вращения при аппроксимации внешних нагрузок также используют разложение в ряд Фурье по окружной координате [7, 8]. Однако при этом размерность задачи не понижается.
В данной работе рассматривается подход, который позволяет свести трехмерную задачу к решению к одномерной задачи, что существенно снижает требования к вычислительным мощностям. Рассматривается задача об определении напряженного состояния тонкостенных оболочечных конструкций в виде тел вращения. Подход основан на интегрировании уравнений теории оболочек (понижение размерности задачи на единицу) и разложении разрешающих функций в ряды Фурье (понижение размерности еще на единицу). В работе использовано разложение в дискретный ряд Фурье по косинусам и синусам, которое описывает произвольные несимметричные механические нагрузки. Такой путь решения позволяет не наносить на поверхность оболочки сетку конечных элементов, от которых зависит объём задачи (размер матрицы жесткости конструкции). В работе используется сетка виртуальных элементов [9], количество которых не увеличивает объём задачи.
Если для описания тонкостенной конструкции использовать оболочечные модели, основанные на различных гипотезах, то напряженное состояние в ортогональной криволинейной системе координат а OP, описывается системой дифференциальных уравнений в частных производных [10, 11, 12]
-=/ Л(а)-+ f(a,P)• (1)
да mV ' dpm f P V '
Здесь n -порядок производных, зависящий от принятых допущений; Л - матрицы, элементы которых определяются через геометрические и
жесткостные параметры оболочки; f - вектор, компоненты которого зависят от приложенных к оболочке поверхностных нагрузок и интегральных характеристик температурного поля. Для классической теории оболочек
искомый вектор N имеет восьмой порядок, а для уточненных моделей -десятый (модель с учетом гипотезы Тимошенко), двенадцатый (модель с учетом сдвига и обжатия) и др.
Система (1) дополняется граничными условиями на контуре а = const. Поскольку рассматриваются оболочки вращения замкнутые в направлении Op, то граничные условия в этом направлении заменяются условиями периодичности. Краевая задача системы (1) за счет периодичности позволяет для всех искомых функций представить решения в виде рядов Фурье по координате p.
Материалы и методы исследования
В данной работе будем использовать классическую теорию оболочек, основанную на гипотезах Кирхгофа-Лява [10]. Координатную поверхность оболочки вращения отнесем к криволинейной ортогональной системе 5,9. На рис. 1 показаны: 5 - длина дуги меридиана, 9 - центральный угол в круге, перпендикулярном оси вращения 2. Следовательно, задача определения напряженного состояния оболочки сводится к системе [13]
-У= £ А (5 +/ М+т (5,9) ,
^ д9 (2)
У = { N,N ,М, ,иг иг, V }, (3)
Яо = К, qz,, 0, 0, 0, 0, 0}Г , Д(в,9еО) = 1, Л(з,0^О) = 0, где У - вектор разрешающих функций; Nr, ^ - радиальное и осевое усилия; иг, и? - радиальное и осевое перемещения; § - сдвигающее усилие;
м5 - меридиональный изгибающий момент; V - окружное перемещение; -угол поворота нормали; Е - единичная матрица, ^ - область контакта.
В отличие от уравнения (1) в правую часть уравнения (2) добавлено еще одно слагаемое, которое учитывает действие контактной
нагрузки <с (по радиальному, осевому и окружному направлению), если оболочка контактирует с основанием.
Внешнюю нагрузку, действующею на оболочку вращения можно представить в виде компонент распределенной нагрузки действующих по
касательной к образующей < и направляющей <9 и по нормали <г к поверхности оболочки. Поскольку обычно для оболочки вращения легко можно подобрать поперечную ось, относительно которой внешние нагрузки будут симметричны или антисимметричны, то их компоненты можно представить в виде разложения [7, 10].
ад
<5 (5,9) = £ <5 (5)С08 к9
I S '
k=0
qe (s,e) = ]Г qe( s)sin kd
k=1
qr (s,0) = £ qr (s) cos кв
' r'
к=0
где ^ - номер гармоники.
В силу периодичности компонент поверхностной нагрузки все функции вектора У можно разложить в ряд Фурье по окружной координате в. Для симметричных компонент вектора ^г,^, м,иг,и2 имеем
ад
ф(5, в) = (5)008кв
к=0
а, для антисимметричных компонент S,v
(5)
к компонент S,v
Фans (s,e) = Уф Г (s)sin кв
k=i (6)
Индексом sim отмечены симметричные а, индексом ans -антисимметричные компоненты вектора.
Учитывая (4) - (6) система в частных производных (2) сводится к ряду систем обыкновенных дифференциальных уравнений в нормальной форме восьмого порядка [10] dY - —
— = Лк (s)Yk + gk (s), (k = 0, 1, 2 ...), ds
AK(s)= ||a(k)(s| (i, j = 1, ... 8), so < s < sL, gk ={gik, ...g8k},
(7)
с граничными условиями
BiY (so ) = b
BY (sL ) = Ä2 . _ _
Здесь Bl, B2 - заданные матрицы, b - заданные векторы. Вектор
нагрузок gk состоит из компонентов внешней и контактной нагрузки. Компоненты контактной нагрузки должны быть заранее найдены из решения контактной задачи. Методы решения одномерных и двумерных
контактных задач и нахождение вектора g k подробно рассмотрены в работе [13].
Для решения краевой задачи системы (7) будем использовать метод дискретной ортогонализации [14]. Поскольку система (7) содержит амплитудные значения разрешающих функций, то все нагрузки и функции (4) - (6) будем строить на дискретном множестве точек, т.е. применять дискретные ряды Фурье [15].
Таким образом, для интегрирования уравнения (7) учитывая (4) значение нагрузки должно быть представлено в виде разложения в дискретный ряд Фурье [15, 16].
N
2
к=1
Поскольку используем подход, основанный на применении дискретных рядов Фурье для функций, то для этого на поверхность оболочки нанесем криволинейную сетку с равным шагом по меридиану и окружности. Таким образом, получим множество виртуальных оболочечных элементов. На рис. 1 показан виртуальный элемент с
а а ^
размерами * и в по меридиану и окружности, нагруженный
распределенной нагрузкой ^. Таким образом, будем считать, что на любом виртуальном элементе известно значение компонент внешней и контактной нагрузки._
Рис. 1. Цилиндрическая оболочка, нагруженная внешней нагрузкой q распределенной
на виртуальном элементе
При решении задач с использованием теории оболочек [10] обычно рассматриваются нагрузки симметричные относительно поперечной оси. Следовательно, при вертикальной оси симметрии X, показанной на рис.1, можно использовать разложение (4).
При решении некоторых задач для тел вращения могут встречаться нагрузки не симметричные относительно поперечной оси. Примером таких нагрузок могут быть ветровые нагрузки, силы инерции при землетрясениях, когда земля получает ускорение и другие виды нагрузок. Следовательно, необходимо аппроксимировать нагрузки, распределенные произвольным образом, которые можно представлять в виде [8].
/
ад ад
qs (s, в) = У q7 (s)cos кв+У qT (s)sin кв
к=0 к=1
ад
¿5,0) = -У qgm (s) sin к в + У q™ (s) cos кв
, (9)
к=1 к=0
ад
q, (s, в) = У q;m (s)cosкв +У q;ns(s)sin кв
к=0 к=1
Разрешающие функции должны при этом быть дополнены мыми - для симметричн
ад
ф (5, в) -уф (фт кв
слагаемыми - для симметричных компонент Nr, Nz, M, wr, Mz
к=1 ,
S V
(10)
и, для антисимметричных компонент
ад
Фsim (s, в) = Уф fm (s)cos к в
к=0
(11)
Следовательно, напряженно-деформированное состояние оболочки будет определяться суммой компонент вектора разрешающих функций (3)
Y = Yi + Y2, (12)
где Y1 - вектор разрешающих функций при разложении (5), (6), а Y2 при разложении (10), (11).
Таким образом, поскольку функция нагрузки известна на множестве виртуальных элементов, периодическая с периодом N, то ее необходимо разложить в ряд Фурье по косинусам и синусам. Среди разнообразных разложений в ряд Фурье дискретных функции будем использовать коэффициенты, дающие хорошую точность аппроксимации [17] .
1 N1 2Ш
ак =Т7 У q cos^r ,
N i=0 N (13)
1 N-1 • , N
b, =— У qi sin-, k = 0, ... —
к N X N 2
где qi - значение нагрузки на каждом виртуальном элементе.
Для демонстрации рассмотренного выше подхода рассмотрим тонкостенную цилиндрическую конструкцию, шарнирно закрепленную по торцам. Конструкция нагружена в трех местах распределенной силой q, действующей по нормали к поверхности оболочки. На рис. 1. показан один нагруженный виртуальный элемент.
При расчете принималось - длина оболочки L=2.40 м, внешний радиус
R = 1,50 м, толщина h = 0,006 м, модуль упругостиE =2,1x10 МПа, коэффициент Пуассона v= 0,3. Граничные условия приняты
Nr=Nz = Ms=S =0, при s =0,
Nr=Nz = Ms=S =0, при s =L.
ад
ад
Окружность разбивалась углом равным в =2°, который стягивает
виртуальный элемент длиной а& , т.е. длина окружности оболочки аппроксимирована 180 виртуальными элементами. Длина оболочки
аппроксимирована 10 виртуальными элементами с длиной а =0.24 м.
Рассмотрен случай, когда на оболочку действуют три силы р =1000 р р
Н, р=1000 Н и р =3000 Н, равномерно распределенных на трех виртуальных элементах. Координаты центра нагруженных элементов 5 =Ь/2 ив=к/3, 5 =Ь/2 и в=ж/2, 5 =Ь/2 и в=2ж/3 соответственно. При действии такой нагрузки не представляется возможным подобрать поперечную ось оболочки, чтобы использовать стандартную процедуру разложения в ряд (4) при действии симметричных нагрузок. Поэтому для решения этой задачи необходимо использовать суперпозицию двух решений (12), рассмотренную выше. При этом численным экспериментом установлено, что в разложении в ряд Фурье (8) необходимо удерживать к = N/2 гармоник.
После интегрирования системы (7) найденное напряженно-деформированное состояние оболочки определяется компонентами
напряжений (а* - меридиональное, °в -окружное, т*в -касательное напряжение) на внешней и внутренней поверхностях оболочки и компонентами перемещений срединной линии [10].
Точность решения задачи зависит от используемой модели оболочки (1), от точности интегрирования системы уравнений (7) и от точности аппроксимирующей нагрузки (13). Эта точность может быть определена сопоставлением с решением аналогичной задачи другим численным методом или с экспериментом. В настоящее время, как отмечено выше, наиболее универсальным численным методом исследования тонкостенных конструкции является метод конечных элементов. Поэтому результаты различных компонентов напряженно-деформированного состояния для исследуемой оболочки сравнивались с результатами, полученными с помощью вычислительного комплекса АКБУБ [1], основанного на методе конечных элементов.
Результаты
На рис.2 сплошной линией показано изменение по окружности
радиального перемещения оболочки иг в зоне приложения сил. На рис.2
штриховой линией показано изменение радиального перемещения иг, полученное при использовании программы АКБУБ. Максимальное
5=\Аи | хм _1 х100% ,
отклонение результатов расчета, вычисленное как 1 г 1 гср (где
и
гср - среднее значение перемещения) не превышает 6.7%. На рисунке видны
три характерных локальных прогиба оболочки в местах приложения
нагрузки прив=^/3, в=ж/2, в =2 п/3. При нагрузкеРз прогиб, естественно, увеличенный.
Рис. 3. Изменение меридионального напряжения а,, на внутренней поверхности оболочки
На рис.3 сплошной линией показано изменение меридионального
напряжения по окружности на внутренней поверхности оболочки а5 в зоне
приложения сил. На рис.3 штриховой линией показано изменение
меридионального напряжения, полученное при использовании программы
АКБУБ. Максимальное отклонение результатов расчета, вычисленное как
5 = Да; | ха; ;1 х 100% . а; ч
ср (где ср - среднее значение напряжения) не превышает
5.5%.
На рис.4 сплошной линией показано изменение окружного
напряжения на внутренней поверхности оболочки а<в по окружности в зоне
приложения сил. На рис.4 штриховой линией показано изменение
окружного напряжения на внутренней поверхности оболочки, полученного
при использовании программы АКБУБ. Максимальное отклонение
5 =| Да в | ха в ; х 100% . а в, ч
в в ср (где в ср - среднее значение напряжения) не
превышает 0.6%.
Рис. 4. Изменение окружного напряжения а на внутренней поверхности оболочки
Из рисунков 2 - 4 видно, что характер напряженно-деформированного состояния исследуемой оболочки вычисленный, с использованием двух различных численных методов полностью совпадает.
Обсуждение результатов
При расчете методом конечных элементов окружность аппроксимирована 180 четырех узловыми элементами БЬеИ-181. Поскольку точки выдачи напряженного состояния в методе конечных элементов предусмотрена по узлам, а в предложенной методе в середине элемента, то происходит сдвижка в представлении информации, что приводит к расхождению результатов в некоторых точках до 6.7%. Следует так же отметить, что данное расхождение находится в пределах допустимой погрешности (7-8%), характерной для метода конечных элементов [1].
Заключение
Таким образом, использование методов теории оболочек, и предложенное разложение разрешающих функций и нагрузки в ряд Фурье, позволяет понизить размерность задачи и, следовательно, решать задачи с использованием небольших вычислительных ресурсов. При этом обеспечивается необходимая точность вычисления по всем компонентам напряженно-деформированного состояния для тонкостенных конструкций.
Библиографический список
1. Басов К. А. АКБУБ: Справочник пользователя / К. А. Басов.
- М.: ДМК Пресс, 2005. - 640 с.
2. Алямовский А. А. SolidWorks / СОЗМОЗ^Ъгкв. Инженерный анализ методом конечных элементов. - М.: ДМК Пресс, 2004. - 432 с.
3. Григолюк Э. И. Контактные задачи теории пластин и оболочек / Э. И. Григолюк, В. М. Толкачев. - М.: Машиностроение, 1980. - 411 с.
4. Джонсон К. Механика контактного взаимодействия / К. Джонсон.
- М.: Мир, 1989. - 510 с.
5. Моссаковский В. И. Контактное взаимодействие элементов оболочечных конструкций / В. И. Моссаковский, В. С. Гудрамович, Е. М. Макеев. - Киев: Наукова думка, 1988. - 288 с.
6. Емельянов И. Г. Определение напряженного состояния и ресурса оболочечной конструкции / И. Г. Емельянов, В. И. Миронов, А. В. Кузнецов // Проблемы машиностроения и надежности машин. - 2007.
- №5. - С. 57-65.
7. Зенкевич О. К. Метод конечных элементов в теории сооружений и механике сплошных сред / О. К. Зенкевич, И. Чанг. - М.: Недра, 1974.
- 240 с.
8. Галлагер Р. Метод конечных элементов. Основы / Р. Галлагер.
- М.: Мир, 1984. - 428с.
9. Емельянов И. Г. Применение виртуальных элементов при определении напряженного состояния оболочек вращения / И. Г. Емельянов, А. В. Кузнецов // Вычислительная механика сплошных сред. - 2014. - Т. 7. - № 3. - С. 245-252.
10. Григоренко Я. М. Методы расчета оболочек: Теория оболочек переменной жесткости / Я. М. Григоренко, А. Т. Василенко. - Киев: Наукова думка, 1981. - Т. 4. - 544 с.
11. Григоренко Я. М. Численно-аналитическое решение задач механики оболочек на основе различных моделей / Я. М. Григоренко, Г. Г. Влайков, А. Я. Григоренко. - Киев: Академпериодика, 2006. - 472 с.
12. Василенко А. Т. Расчет параметров напряженного состояния конструктивных элементов из композиционных материалов на основе оболочечных моделей / А. Т. Василенко, Г. П. Голуб, Я. М. Григоренко // Расчеты на прочность. Вып. №30. - М.: Машиностроение, 1989. - С. 87-96.
13. Емельянов И. Г. Контактные задачи теории оболочек / И. Г. Емельянов. - Екатеринбург: УрО РАН , 2009. - 185с.
14. Годунов С. К. О численном решении краевых задач для систем линейных обыкновенных дифференциальных уравнений / С. К. Годунов // Успехи математических наук. - 1961. - Вып. 16. - № 3. - С. 171-174.
15. Григоренко Я. М. Решение краевых задач о напряженном состоянии упругих тел сложной геометрии и структуры с применением дискретных рядов Фурье / Я. М. Григоренко // Прикладная механика. - 2009. - Вып. 45. - №5. - С. 3-51.
16. Корн Г. Справочник по математике для научных работников и инженеров / Г. Корн, Т. Корн. - М.: Наука, 1968. - 720 с.
17. Ба-хуссейн А. А. Дискретное преобразование Фурье. - URL: http://ilab.xmedtest.net/?q=node/3740 (01.02.2015).
References
1. Basov K. A. ANSYS: Spravochnik pol'zovatelya. [ANSYS: User's reference]. Moscow, 2005. 640 p.
2. Alyamovskiy A. A. SolidWorks . COSMOSWorks. Inzhenernyy analiz metodom konechnykh elementov. [SolidWorks. COSMOSWorks. Engineering analysis by the finite element method]. Moscow, 2004. 432 p.
3. Grigolyuk E. I. & Tolkachev V. M. Kontaktnyye zadachi teorii plastin i obolochek [Contact problems of the theory of plates and shells]. Moscow, 1980. 411 p.
4. Dzhonson K. Mekhanika kontaktnogo vzaimodeystviya. [The mechanics of contact interaction]. Moscow, 1989. 510 p.
5. Mossakovskiy V. I., Gudramovich V. S. & Makeyev Ye. M. Kontaktnoye vzaimodeystviye elementov obolochechnykh konstruktsiy.[ Contact interaction of elements of shell structures]. Kiyev, 1988. 288 p.
6. Yemel'yanov I. G., Mironov V. I. & Kuznetsov A. V. Problemy mashinostroyeniya i nadezhnosti mashin - Problems of machine building and reliability of machines, 2007, no. 5, pp. 57-65.
7. Zenkevich O. K. & Chang I. Metod konechnykh elementov v teorii sooruzheniy i mekhanike sploshnykh sred [The finite element method in the theory of structures and the mechanics of continuous media]. Moscow, 1974. 240 p. 8.
Gallager R. Metod konechnykh elementov. Osnovy [Finite element method. Basics]. Moscow, 1984. 428p.
9. Yemel'yanov I. G. & Kuznetsov A. V. Vychislitel'naya mekhanika sploshnykh sred - Computational mechanics of continuous media, 2014, vol. 7, no. 3, pp. 245-252.
10. Grigorenko Ya. M. & Vasilenko A. T. Metody rascheta obolochek: Teoriya obolochek peremennoy zhestkosti [Methods for calculating shells: The theory of shells of variable rigidity]. Kiyev, 1981, T. 4. 544 p.
11. Grigorenko Ya. M., Vlaykov G. G. & Grigorenko A. Ya. Chislenno-analiticheskoye resheniye zadach mekhaniki obolochek na osnove razlichnykh modeley [Numerically-analytical solution of the problems of shell mechanics on the basis of different models]. Kiyev, 2006. 472 p.
12. Vasilenko A. T., Golub G. P. & Grigorenko Ya. M. Raschety na prochnost' - Strength calculations, Moscow, 1989, vol. 30, pp. 87-96.
13. Yemel'yanov I. G. Kontaktnyye zadachi teorii obolochek [Contact problems in the theory of shells]. Yekaterinburg, 2009. 185 p.
14. Godunov S. K. Uspekhi matematicheskikh nauk - Success in Mathematical Sciences, 1961. vol. 16, no. 3, pp. 171-174. (In Russian)
15. Grigorenko Ya. M. Prikladnaya mekhanika - Applied Mechanics, 2009, no. 5, pp. 3-51.
16. Korn G. & Korn T. Spravochnik po matematike dlya nauchnykh rabotnikov i inzhenerov [A Handbook of Mathematics for Scientists and Engineers]. Moscow, 1968. 720 pp.
17. Ba-khusseyn A. A. Diskretnoyepreobrazovaniye Fur'ye [Discrete Fourier Transform]. URL: http://ilab.xmedtest.net/?q=node/3740 (01/02/2015).
Сведения об авторах:
ЕМЕЛЬЯНОВ Игорь Георгиевич, главный научный сотрудник Института машиноведения УрО РАН, профессор Уральского государственного университета путей сообщения. Тел.: +7(922)2073447. E-mail: emelyanov@imach.uran.ru
КУЗНЕЦОВ Алексей Викторович, научный сотрудник Института машиноведения УрО
РАН. Тел.: +7(922)2912357.
E-mail: Alekseikuz031082@rambler.ru
Information about authors:
Igor G. EMELIANOV, chief research officer of Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences, professor of the Ural State University of Railway Transport. Tel.: +7 (922) 2073447. E-mail: emelyanov@imach.uran.ru
Alexey V. KUZNETSOV, Research Fellow, Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences. Tel.: +7 (922) 2912357. E-mail: Alekseikuz031082@rambler.ru
© EME^kSHOB HT., KY3HE^B A. B., 2017
