CC BY
123
УДК004.925:539.3:62-41
Кудин А. В.1, Чопоров С. В.2
Преподаватель, Запорожский национальный университет, Украина, E-mail: [email protected] 2Канд. техн. наук, доцент, Запорожский национальный университет, Украина
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ИЗГИБА КРУГЛОЙ ТРЕХСЛОЙНОЙ ПЛАСТИНЫ С ИСПОЛЬЗОВАНИЕМ _АНАЛИТИЧЕСКОГО И ЧИСЛЕННОГО ПОДХОДОВ_
Предложены уравнения изгиба трехслойных круглых пластин симметричного строения с изотропными наружными слоями и нелинейно-упругим заполнителем. На основе полученных уравнений разработана компьютерная модель в системе Maple. Построена трехмерная конечно-элементная модель круглой трехслойной пластинки. В качестве примера, рассмотрена задача осесимметричного поперечного изгиба трехслойной круглой пластинки в линейной постановке, а также задача осесимметричного поперечного изгиба однослойной пластины в линейной и нелинейной постановке.
Ключевые слова: компьютерная модель, круглая трехслойная пластина, аналитическое решение, метод конечных элементов.
ВВЕДЕНИЕ
Компьютерное моделирование позволяет заменить исследование физического образца вычислительным экспериментом и является одним из методов исследования сложных систем. При этом в основе компьютерной модели, как правило, лежат математические уравнения и системы уравнений. Основные подходы при получении определяющих соотношений модели: аналитический подход на базе распространенных систем компьютерной алгебры, численные исследования. Следовательно, актуальным является сопоставление результатов таких исследований. Совместное использование аналитического и численного подхода позволяет рассмотреть объект с разных сторон, используя разные гипотезы, сформировать область допустимых решений.
Трехслойные элементы конструкций широко применяются в авиа- и судостроении, космической промышленности, гражданском строительстве, радиоэлектронике и др. отраслях промышленности. Поэтому актуальной является проблема разработки эффективных методов расчета напряженно-деформированного состояния трехслойных элементов конструкций, а также обобщения классических теорий с применением уточненных моделей, отражающих поведение современных материалов.
В настоящей работе приводится вариант уравнений изгиба трехслойных пластин симметричного строения с изотропными наружными слоями и нелинейно-упругим по [1] материалом заполнителя; построена компьютерная модель на базе пакета символьных вычислений Maple (аналитическая модель).
Предложена трехмерная конечно-элементная модель трехслойной пластинки, на основе которой получены численные оценки параметров состояния трехслойной конструкции. На примере задачи осесимметричного поперечного изгиба трехслойной круглой пластинки проводится сравнение полученных моделей с известны-
ми работами [2, 3, 4]. В качестве частного случая решена задача осесимметричного поперечного изгиба однослойной пластины в линейной и нелинейной постановке; результаты сравниваются с работами [1, 5].
1 АНАЛИТИЧЕСКАЯ МОДЕЛЬ
Рассмотрим круглую трехслойную пластинку, которая подвергается воздействию поперечной нагрузки д(г), распределеной симметрично относительно оси, проходящей через центр пластины. Внешние слои пластинки толщиной 8 изготовлены из изотропного материала, подчиняющегося закону Гука. Средний слой пластины толщиной 2Н выполнен из нелинейно-упругого по [1] изотропного материала.
Деформированное состояние пластинки определяется радиальным перемещением и (г) и прогибом w(г).
Дифференциальные уравнения равновесия круглой трехслойной пластины в общем случае произвольной нагрузки приведены в [5]. Принимая во внимание симметрию деформированного состояния пластины и отсутствие угловых перемещений в данной постановке задачи, нелинейные дифференциальные уравнения равновесия с переменными коэффициентами имеют вид
Л{ ^ ггг + Лу2м> гг + Л13и гг + Л[4w г + Л15и г + + Л16и (г ) + Ф1(и 3, w3,...) = 0,
(1)
Л2\W гггг + Л22w ггг + Л23и ггг + Л24w гг + Л25и гг +
+Л26 wг + Л27и г + Л28и (г)-щ + Ф2(и3, w3,...) = 0.
Нелинейные члены уравнений (1.4) Ф^3, м?,..) и
Ф2(и3,w3,...), а также коэффициенты Лп..Л^, Л21.Л28 приведены в [6], здесь не приводятся ввиду их громоздкости.
© Кудин А. В., Чопоров С. В., 2014 DOI 10.15588/1607-3274-2014-1-11
Граничные условия для свободного опирания имеют следующий вид [4]:
w(r) = 0, w,rr + — Wr = 0, ur + — u(r) = 0 при г = R. r ' r
Граничные условия при защемлении на контуре имеют следующий вид [4]:
w(r) = 0, w,r = 0, u(r) = 0 при r = R.
Одним из подходов к решению системы нелинейных уравнений (1) является метод разложения по малому параметру [1, 7, 8], который позволяет рассматривать нелинейные задачи как последовательность уточняющих друг друга решений систем линейных уравнений. В свою очередь, рассматриваемые системы линейных уравнений в данном случае могут быть сведены к уравнению Бесселя, решение которого известно [2].
Построение компьютерной модели с помощью пакета Maple проводилось в несколько этапов: описание напряженного состояния пластины; описание деформированного состояния пластины; получение функционала энергии; вывод дифференциальных уравнений равновесия на основе принципа Лагранжа аналогично [9]. Приведем команды пакета Maple, описывающие напряженное состояние пластины.
Закон Гука для внешних слоев пластины [10]:
sigma[1,r]:=E/(1-muA2)*(epsilon[1,r]+mu*epsilon[1,phi]);
sigma[1,phi]:=E/(1-muA2)*(epsilon[1,phi]+mu*epsilon[1,r]);
tau[1,rphi]:=G12*epsilon[1,rphi];
sigma[2,r]:=E/(1-muA2)*(epsilon[2,r]+mu*epsilon[2,phi]);
sigma[2,phi]:=E/(1-muA2)*(epsilon[2,phi]+mu*epsilon[2,r]);
tau[2,rphi]:=G12*epsilon[2,rphi];
Закон упругости для среднего слоя [4]:
tau[3,rphi]:=G[3]*gamma1 *epsilon[3,rphi]; tau[3,rz]:=G[3]*gamma1*epsilon[3,rz]; tau[3,phiz]: =G[3] *gamma 1*epsilon[3,phiz]; sigma[3,r]:=3*K[3]*chi1*epsilon[0]+2*G[3]*gamma1*(epsilon[3,r]-epsilon[0]);
sigma[3,phi]:=3*K[3]*chi1*epsilon[0]+2*G[3]*gamma1*(epsilon[3,phi]-epsilon[0]);
Использование системы компьютерной алгебры Maple существенно упрощает оперирование с громоздкими выражениями и позволяет автоматизировать процесс получения дифференциальных уравнений задачи из основных уравнений теории упругости. 2 КОНЕЧНО-ЭЛЕМЕНТНАЯ МОДЕЛЬ Рассмотрим решение задачи исследования напряженно-деформированного состояния равномерно нагруженной круглой трехслойной пластинки в конечно-элементной постановке. Учитывая симметричность поставленной задачи достаточно исследовать лишь сектор пластинки (рис. 1), у которого для узлов, соответствующих грани OA, запрещено перемещение в направлении 2, а узлам на гране OB запрещено перемещение в на-
2 Jk .1
В
О Л
--?
1
Рис. 1. Геометрическая модель пластинки
правлении 1. При этом учет кусочной неоднородности (слоистости) пластинки приводит к необходимости исследования конструкции в трехмерной постановке.
Используя аппарат теории Я-функций В. Л. Рвачева [11, 12] геометрическая модель может быть формализована при помощи выражения
/ (х,у,2) = (Г2 - X2 - у2 )л X Л у Л (2 + Н + §1) л(Н + §2 - 2),
где г - радиус пластинки; х, у, 2 - пространственные ко-
^2 2 Х1 + Х2 -
Я-конъюнкция.
Функция / (х,у,2) принимает значения больше нуля во внутренних точках сектора пластинки, равна нулю на его границе и меньше нуля во внешних точках (боковые грани ограничены плоскостями х = о и у = 0, нижняя грань плоскостью г = -(Н + §2), а верхняя 2 = Н + §1). Следовательно, / (ху,?} может быть использована в качестве предиката (управляющей функции) для автоматической генерации системы конечных элементов.
Для обеспечения большей точности моделирования при построении дискретной модели необходимо учесть структуру слоев: границы слоев должны аппроксимироваться границами элементов (т. е. вдоль границы слоя материала должна проходить граница слоя конечных элементов). Удовлетворение этого требования совместно с учетом ограничения на соотношение размеров сторон при построении сетки шестигранных элементов [13-15] приводит к относительно большому числу конечных элементов. Например, при г = 0,4, §1 = §2 = 0,001 и 2Н = 0,016 сетка состоит из 102277 узлов и 94032 конечных элементов. Опытным путем установлено, что если соотношение сторон элементов превышает 10, то вычислительный процесс дает значительно большую погрешность.
С учетом трехслойности, формула для вычисления локальной матрицы жесткости [Ке] примет вид
[ ке ] = \\\вт овахауа2,
где
В
д_
дх
0 0
д_
ду
0 — 0
0 0 —
А д.
ду дх
0
0
Л
дг
дг ду
0
д_
дх.
#1 # 2 . оо ¡4 0 0 . . 0 0 0 . . 0 "
0 0 . . 0 #1 #2 . . #8 0 0. . 0
0 0. . 0 0 0 . .. 0 #1 #2 . . #8.
Б =
ц ц
1
1 -ц 1 -ц ц 1 Ц
1 -ц 1 -ц
ц ц
Е(-ц) .
(1 + ц)(1 - 2ц)
0 0
1 -ц 1 -ц
0 0
0 0
0 0
1
0
0
1 - 2ц
2 (1 -ц)
0 0
0
0
0
1 - 2ц
2 (1 -ц)
0
0 0 0 0 0
1 - 2ц
2 (1 -ц)
Е =
(Е^, г > к или г < -к,
|Е, - к < г < к,
ц^2 ,г> к или г < -к,
ц, - к < г < к,
3 СОПОСТАВЛЕНИЕ ЗНАЧЕНИИ ПОПЕРЕЧНОГО ИЗГИБА КРУГЛЫХ ОДНОСЛОЙНЫХ ПЛАСТИН
В качестве частного случая рассмотрим поперечный изгиб круглой однослойной пластины (толщина пластинки к = 18 -10 м, радиус Я = 0,4 м; модуль сдвига и коэффициент Пуассона материала - О = 8 -104 МПа и
ц = 0,27 соответственно) и сопоставим полученные значения с решением Тимошенко [5] для равномерно нагруженной круглой пластинки. В табл. 1 приводятся результаты такого сопоставления. Через ^шах обозначен максимальный прогиб в центре пластины; q - величина равномерно распределенной поперечной нагрузки. Модели, используемые при расчетах, обозначены в табл. 1 следующим образом: 1 - аналитическая модель, 2 - конечно-элементная модель, 3 - модель Тимошенко круглой однослойной пластины [5].
Из сопоставления значений прогибов в табл. 1 можно сделать вывод о соответствии построенных моделей решению Тимошенко [5]. Прогибы, полученные по аналитической модели, с точностью до шестого знака после запятой совпадают с прогибами модели Тимошенко. Для свободного опирания отклонения прогибов конечно-элементной модели от прогибов Тимошенко находятся в пределах 2,7 %-2,8 %; для защемления - 3,7 %-3,9 %. Полученные результаты позволяют сделать вывод об адекватности построенных моделей модели Тимошенко и возможности их применения для дальнейшего исследования прогибов слоистых пластин. Разница численных значений прогибов конечно-элементной модели и модели [5] объясняется использованием гипотезы Кирхгофа-Лява в модели Тимошенко, в отличие от конечно-элементной модели.
Рассмотрим задачу изгиба однослойной круглой пластины из нелинейно-упругого материала. Сравним полученные значения с результатами Каудерера [1]. Примем следующие параметры пластины: толщина пластин-
-3
ки к = 18 -10 м, радиус пластинки Я = 0,4 м; модуль
где Е12 и ц12 - модуль упругости и коэффициент Пуассона внешних слоев; Е и ц - модуль упругости и коэффициент Пуассона внутреннего слоя.
Матрица жесткости вычисляется для каждого конечного элемента и затем производится ее ансамблирова-ние (суммирование по специальным правилам) в глобальную матрицу жесткости, которая совместно с вектором узловых нагрузок определяет СЛАУ для вычисления перемещений в узлах.
При этом найденные перемещения в третьем направлении будут соответствовать прогибам пластинки.
Таблица 1. Изгиб круглой однослойной пластинки в линейной постановке
д , МПа 3 ^шах = 10 м
Тип закрепления
Свободное опирание Защемление
Модель Модель
1 2 3 1 2 3
0,05 0,779 0,800 0,779 0,188 0,195 0,188
0,07 1,090 1,121 1,090 0,263 0,273 0,263
0,09 1,402 1,441 1,402 0,338 0,351 0,338
0,11 1,714 1,761 1,714 0,413 0,429 0,413
х
X
сдвига и модуль объемной деформации материала -О = 2,77 -104 МПа, К = 6 -104 МПа, параметр у2 =-3,878-105 [1, 16]. Вычисления производятся для случая защемления контура пластинки. Модели, используемые при расчетах, обозначены в табл. 2 следующим образом: 1 - аналитическая модель, 2 - решение Каудерера [1].
Таблица 2. Изгиб круглой однослойной пластинки в нелинейной постановке
д , МПа ^шах = ^(0) , 103 м
Модель
1 2
0,05 0,529 0,523
0,07 0,756 0,740
0,09 1,000 0,977
0,11 1,275 1,329
При вычислении значений прогибов табл. 2 учитывались три приближения в разложении по методу малого параметра.
Как видно из табл. 2, отклонения значений прогибов аналитической модели от прогибов [1] находятся в пределах 1 %-4 %, что говорит о возможности применения аналитической модели для решения нелинейных задач теории пластин.
4 СОПОСТАВЛЕНИЕ ЗНАЧЕНИЙ ПОПЕРЕЧНОГО ИЗГИБА КРУГЛЫХ ТРЕХСЛОЙНЫХ ПЛАСТИН
Поскольку нулевое приближение дает наиболее существенный вклад в разложение искомых перемещений по малому параметру, исследуем далее задачу о поперечном изгибе круглой трехслойной пластины в линейной постановке и проведем сравнительный анализ полученных решений с известными результатами.
Рассмотрим поперечный изгиб круглой трехслойной пластины со следующими параметрами: толщина среднего слоя 2к = 16 -10 м, толщина внешних слоев 8 = 1.10_3 м, радиус пластинки Я = 0,4 м; модуль сдвига и коэффициент Пуассона внешних слоев - О = 8 • 104 МПа
и ц = 0,27 соответственно, модуль сдвига и модуль объемной деформации заполнителя - О = 2,77 -104 МПа,
К = 6 • 104 МПа. Тангенциальные напряжения изменяются линейно по толщине среднего слоя.
В табл. 3 приводятся значения максимального прогиба в центре пластины в линейной постановке при различных граничных условиях. Значения максимального прогиба в центре пластинки, полученные из работ [2-4], сопоставляются со значениями прогибов полученными из решения уравнений (1) и прогибами конечно-элементной модели. Модели, используемые при расчетах, обозначены в табл. 3 следующим образом: 1 - аналитическая модель, 2 - конечно-элементная модель; 3 - модель И. А. Михайлова [3], 4 - модель А.П. Прусакова [4], 5 - модель А. Г. Горшкова и др. [2].
Из сравнения полученных значений максимальных прогибов приведенных в табл. 3 можно сделать вывод о соответствии результатов полученных на основе построенных в работе аналитической и конечно-элементной моделей результатам известных моделей [2-4]. В силу различных упрощающих гипотез значения прогибов табл. 3 можно рассматривать как разные оценки области точного решения. Наблюдаются следующие значения отклонений рассматриваемых моделей: аналитическая модель и конечно-элементная модель - 8 % -11 %; модель И. А. Михайлова [3] - 1 %-45 %; модель А. П. Прусакова -30 %-80 %; модель А. Г. Горшкова [2] - 16%-55% для случая защемления.
На рис. 2 приведены прогибы срединной плоскости трехслойной пластинки для случая свободного опира-ния (рис. 2, а) и для случая защемления на контуре
(рис. 2, б); значение нагрузки принимается д = 0,11 МПа.
Отклонения табл. 3 объясняются тем, что при построении аналитической модели и моделей [2-4] применялись различные упрощающие гипотезы и как следствие, в тензорах напряжений и деформаций учитывалось различное количество компонент. Так, при построении конечно-элементной модели учитывались все компоненты тензора напряжений и тензора деформаций (трехмерная постановка задачи), а в аналитической модели - в зависимости от слоя от (трех до пяти).
Таблица 3. Изгиб трехслойной круглой пластинки в линейной постановке
д , МПа 3 ^шах = ^С0^ 10 м
Тип закрепления
Свободное опирание Защемление
Модель Модель
1 2 3 4 1 2 3 4 5
0,05 1,606 1,449 1,430 2,621 0,388 0,357 0,348 0,635 0,291
0,07 2,248 2,029 2,001 3,670 0,543 0,499 0,487 0,889 0,407
0,09 2,891 2,609 2,573 4,718 0,698 0,642 0,626 1,142 0,523
0,11 3,533 3,188 3,145 5,767 0,853 0,784 0,765 1,396 0,639
\ч(г)
0.005
0.004
0.003
0.002
0.001
Г 'ч.
\ \
\ Ч \
\ х \ Чч '
% Ч \ \
ч>(г)
0.0012 0.0010 0.0008 0.0006 0.0004 0.000. 0
Ч.
\ \
ч \
\ \\ \
Ч N \ Л 1
1А._ \
0 0.1 0.2 0.3
г
0.4
а)
- аналитическая модель,
конечно-элементная модель,
0 0.1 0.2 0.3 г
б)
- модель Михайлова, • » • ■
0.4
модель
Прусакова, - модель Горшкова
Рис. 2. Прогиб срединной плоскости в случае распределенной нагрузки
ВЫВОДЫ
Разработана компьютерная модель на базе нового варианта дифференциальных уравнений равновесия круглых трехслойных пластин. Пластины считаются симметричными по толщине с изотропными наружными слоями и нелинейно-упругим по [1] материалом заполнителя. Построена трехмерная конечно-элементная модель.
В качестве примера рассмотрена задача осесиммет-ричного поперечного изгиба однослойной круглой пластинки в линейной и нелинейной постановке. Также рассмотрена задача осесимметричного поперечного изгиба трехслойной круглой пластинки в линейно-упругой постановке.
Отметим, что приведенные в табл. 3 модели используют различные упрощающие гипотезы, чем обусловливается разница полученных значений прогибов. Так, при построении конечно-элементной модели использовалась трехмерная постановка задачи. В модели А. Г. Горшкова [2] и И. А. Михайлова [3] для внешних слоев принимается гипотеза Кирхгофа-Лява, для заполнителя же принимаются различные гипотезы о распределении тангенциальных напряжений по толщине слоя. Существенные отклонения результатов расчетов по уравнениям модели А. П. Прусакова [4] объясняются тем, что здесь заполнитель считается легким, в отличие от других вышеупомянутых моделей.
Сравнение полученных значений прогибов на основе аналитической и конечно-элементной моделей с работами [1-5] свидетельствует об адекватности построенных моделей и возможности решения на их основе более широкого класса задач статики и динамики трехслойных пластин.
На основе аналитической модели могут быть получены компактные расчетные формулы для применения в инженерной практике. В дальнейшем, при решении систем дифференциальных уравнений первого и последующего приближений разложения по малому параметру, становится возможным учет влияния нелинейной упругости материала заполнителя на деформированное состояние трехслойной пластинки.
Недостатком трехмерной конечно-элементной модели является требование относительно большого объема вычислений, по сравнению с приведенными аналитическими моделями. Однако такой подход может использоваться для расчета многослойных пластин нестандартной формы и для проверки моделей использующих упрощающие гипотезы. Сравнение результатов аналитических решений с численными позволяет делать заключение о правильности выбранных гипотез.
Перспективы дальнейшего исследования связаны с рассмотрением задач нелинейной устойчивости и динамики трехслойных элементов конструкций, а также созданием гибридных аналитико-вычислительных подходов на базе конечных элементов.
СПИСОК ЛИТЕРАТУРЫ
1. Каудерер, Г. Нелинейная механика / Каудерер Г. ; пер. с нем. - М. : Изд-во иностр. лит., 1961. - 777 с.
2. Горшков, А. Г. Механика слоистых вязкоупругопласти-ческих элементов конструкций / А. Г. Горшков, Э. И. Ста-ровойтов, А. В. Яровая. - М. : ФИЗМАТЛИТ, 2005. -576 с.
3. Михайлов, И. П. Некоторые задачи осесимметричного изгиба круглых трехслойных пластин с жестким заполнителем // Труды Ленинградского кораблестроительного института. - 1969. - Вып. 66. - С. 125-131.
4. Прусаков, А. П. Некоторые задачи изгиба круглых трехслойных пластин с легким заполнителем // Тр. конф. по теор. пластин и оболочек. - 1961. - № 1. - С. 293-297.
5. Тимошенко, С. П. Пластинки и оболочки / С. П. Тимошенко, С. Войновский-Кригер. - М. : Физматгиз, 1964. -636 с.
6. Кудин, А. В. Применение метода малого параметра при моделировании изгиба симметричных трехслойных пластин с нелинейно-упругим заполнителем / А. В. Кудин, Ю. Н. Тамуров // Вюник Схщноукрашського нацюналь-ного ушверситету iменi Володимира Даля. - 2011. - № 11 (165). - С. 32-40.
7. Тамуров, Ю. Н. Вариант обобщенной теории трехслойных пологих оболочек с учетом обжатия физически нелинейного заполнителя // Прикладна механика. - 1990. - Т. 26, № 12. - С. 39-45.
8. Тамуров, Ю. Н. Уравнения изгиба и устойчивости трехслойных оболочек с ортотропными и нелинейно-упругими свойствами материалов // Исслед. по теор. пластин и оболочек. - 1990. - № 20. - С. 102-112.
9. Liu, Renhuai Nonlinear Bending of Circular Sandwich Plates / Liu Renhuai // Applied Mathematics and Mechanics. English Edition. - 1981. - Vo1. 2 N°. 2. - p. 189-208. Стаття надшшла до редакци 28.03.2014.
Кудш О. В.1, Чопоров С. В.2
'Викладач, 3anopi3bKM нацюнальний ушверситет, Укра!на
2Канд. техн. наук, доцент Запорiзький нацюнальний ушверситет, Украша
КОМП'ЮТЕРНЕ МОДЕЛЮВАННЯ ВИГИНУ КРУГЛО1 ТРИШАРОВО1 ПЛАСТИНИ З ВИКОРИСТАННЯМ АНАЛ1ТИЧНОГО ТА ЧИСЕЛЬНОГО П1ДХОД1В
Запропоновано рiвняння вигину тришарових круглих пластин симетрично! будови з iзотропними зовшшшми шарами i нелшшно-пружним заповнювачем. На основi отриманих рiвнянь розроблено комп'ютерну модель у системi Maple. Побудована тривимiрна сюнченно-елементна модель кругло! тришарово! пластинки. В якост прикладу, розглянуто задачу вююиметричного поперечного вигину тришарово! кругло! пластинки в лшшнш постановщ, а також задача вююиметричного поперечного вигину одношарово! пластини в лшшнш i нелшшнш постановщ.
Ключовi слова: комп'ютерне моделювання, кругла тришарова пластина, аналтчний розв'язок, метод скiнченних еле-ментiв.
Kudin A. V.1, Choporov S. V.2
'Assistant Professor, Zaporizhzhya National University, Ukraine
2Ph.D., Associate professor, Zaporizhzhya National University, Ukraine
COMPUTER SIMULATION OF BENDING CIRCULAR SANDWICH PLATES USING ANALYTICAL AND NUMERICAL APPROACHES
The paper describes equations of circular symmetrical sandwich plates bending for plates with isotropic outer layers and nonlinear elastic filler. Two computer models of circular sandwich plates are proposed: 1) analytical model on the basis of Maple CAS, 2) three-dimensional FEM-model. We compare results of analytical model with results of numerical simulation using two problems: 1) the problem of axisymmetric transverse bending of a sandwich plate, 2) the problem of axisymmetric transverse bending of a single-layer plate with linear and nonlinear formulation.
The first section of the article describes basic differential equations and analytical model constructing method. FEM-model is described in the second section. Numerical results of computer simulation are showed in the last section of the article.
Keywords: computer simulation, sandwich symmetrical plate, circular plate, nonlinear-elastic core, axisymmetric bending, FEMmodel.
10. Амбарцумян, С. А. Теория анизотропных пластин: Прочность, устойчивость и колебания. - 2-е изд., перераб. и доп. - М. : Наука, 1987. - 360 с.
11. Рвачев, В. Л. Новые подходы к построению уравнений трехмерных локусов с помощью R-функций / [Рвачев В. Л., Толок А. В., Уваров Р. А., Шейко Т. И.] // В юник Запорiзького державного ушверситету. - 2000. - № 2. -С. 119-130.
12. Рвачев, В. Л. Введение в теорию R-функций / В. Л. Рвачев, Т. И. Шейко // Проблемы машиностроения. - 2001. -Т. 4, № 1-2. - С. 46-58.
13. Гоменюк, С. И. Дискретизация трехмерных областей, заданных R функциями, на шестигранные конечные элементы / С. И. Гоменюк, С. В. Чопоров // Вестник Херсонского национального технического университета. - 2011. -№ 3 (42). - С. 146-153.
14. Schneiders, R. A Grid-based Algorithm for the Generation of Hexahedral Element Meshes / R. Schneiders // Engineering with Computers. - 1996. - No. 12. - P. 168 177.
15. Schneiders, R. Octree-based Generation of Hexahedral Element Element Meshes / R. Schneiders, R. Schindler, F. Weiler // 5th Annual International Meshing Roundtable. -1996. - P. 205-216.
16. Цурпал, И. А. Расчет элементов конструкций из нелинейно-упругих материалов. К. : Техника, 1976. - 176 с.
REFERENCES
1. Kauderer G. Nelinejnaya mexanika, Per. s nem. Moscow, Izd-vo inostr. lit., 1961, 777 p.
2. Gorshkov A. G., Starovojtov E'. I., Yarovaya A. V. Mexanika sloistyx vyazkouprugoplasticheskix e'lementov konstrukcij. Moscow, FIZMATLIT, 2005, 576 p.
3. Mixajlov I. P. Nekotorye zadachi osesimmetrichnogo izgiba kruglyx trexslojnyx plastin s zhestkim zapolnitelem, Trudy Leningradskogo korablestroitel'nogo instituta, Vyp. 66, 1969, pp. 125-131.
4. Prusakov A. P. Nekotorye zadachi izgiba kruglyx trexslojnyx plastin s legkim zapolnitelem, Tr. konf. po teor. plastin i obolochek, 1961, No. 1, pp. 293-297.
5. Timoshenko S. P., Vojnovskij-Kriger S. Plastinki i obolochki. Moscow, Fizmatgiz, 1964, 636 p.
6. Kudin A. V., Tamurov Yu. N. Primenenie metoda malogo parametra pri modelirovanii izgiba simmetrichnyx trexslojnyx plastin s nelinejno-uprugim zapolnitelem, Visnik Sxidnoukraïns'kogo nacional'nogo universitetu imeni Volodimira Dalya. 2011, No. 11(165), pp. 32-40.
7. Tamurov Yu. N. Variant obobshhyonnoj teorii tryoxslojnyx pologix obolochek s uchyotom obzhatiya fizicheski nelinejnogo zapolnitelya, Prikl. mexanika, 1990, Vol. 26, No. 12, pp. 39-45.
8. Tamurov Yu. N. Uravneniya izgiba i ustojchivosti trexslojnyx obolochek s ortotropnymi i nelinejno-uprugimi svojstvami
materialov, Issled. po teor. plastin i obolochek, 1990, No. 20, pp. 102-112.
9. Liu Renhuai Nonlinear Bending of Circular Sandwich Plates, Applied Mathematics and Mechanics, English Edition, 1981, Vol. 2, No. 2, pp. 189-208.
10. Ambarcumyan S. A. Teoriya anizotropnyx plastin: Prochnost', ustojchivost' i kolebaniya. 2-e izd., pererab. i dop. Moscow, Nauka, 1987, 360 p.
11. Rvachev V. L., Tolok A. V., Uvarov R. A., Shejko T. I. Novye podxody k postroeniyu uravnenij trexmernyx lokusov s pomoshh'yu R funkcij, Visnik Zaporiz'kogo derzhavnogo universitetu, 2000, No. 2, pp. 119-130.
12. Rvachev V. L., Shejko T. I. Vvedenie v teoriyu R-funkcij, Problemy mashinostroeniya, 2001, Vol. 4, No. 1-2, pp. 4658.
13. Gomenyuk S. I., Choporov S. V. Diskretizaciya trexmernyx oblastej, zadannyx R funkciyami, na shestigrannye konechnye e'lementy, VestnikXersonskogo nacional 'nogo texnicheskogo universiteta, 2011, No. 3(42), pp. 146-153.
14. Schneiders R. A Grid-based Algorithm for the Generation of Hexahedral Element Meshes, Engineering with Computers, 1996, No. 12, pp. 168 177.
15. Schneiders R., Schindler R., Weiler F. Octree-based Generation of Hexahedral Element Element Meshes, 5th Annual International Meshing Roundtable, 1996, pp. 205-216.
16. Curpal I.A. Raschet e'lementov konstrukcij iz nelinejno-uprugix materialov. Kiev, Texnika, 1976, 176 p.
