УДК 539.3
Численное моделирование и экпериментальное исследование деформирования упругопластических пластин при смятии
© Ю.И. Димитриенко1, Е.А. Губарева1, С.В. Сборщиков1,
2 2 В С. Ерасов , Н О. Яковлев
1 МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
2
2 ФГУП «ВИАМ» ГНЦ РФ, Москва, 119049, Россия
Предложена методика численного конечно-элементного решения задачи овализа-ции, которую используют при экспериментальной отработке новых материалов для авиационной промышленности, в целях определения сопротивления деформированию элементов конструкций с наличием концентраторов напряжений, главным образом, соединительных элементов. Методика основана на трехмерном конечно-элементном решении задачи упругопластического деформирования пластин с отверстием при смятии и предназначена для сокращения экспериментальных исследований путей замены их на численные эксперименты. Используется модель малых упругопластических деформаций Ильюшина. Представлены результаты численного моделирования трехмерного напряженно-деформированного состояния упругопластических пластин при смятии, а также результаты экспериментальных исследований деформирования пластин из алюминиевого сплава 163. Показано, что результаты численного и экспериментального моделирования деформирования пластин при смятии достаточно хорошо совпадают.
Ключевые слова: задача овализации, смятие пластин, численное моделирование, метод конечных элементов, упругость, пластичность, диаграммы деформирования.
Введение. При экспериментальной отработке новых материалов для авиационной промышленности в целях определения сопротивления деформированию элементов конструкций с наличием концентраторов напряжений, главным образом соединительных элементов, в промышленности применяют стандарт ЛБТМ Б 5961/Б 5961М, в соответствии с которым проводят испытания образца в виде пластины с отверстием. Пластина одним краем жестко закреплена, на другом крае в отверстие вставлен цилиндрический стержень, к которому приложена нагрузка (рис. 1). При определенной нагрузке Р пластина начинает деформироваться — возникает эффект овализации отверстия, или происходит смятие пластины. По стандарту ЛБТМ Б 5961/Б 5961М используется пластина прямоугольного сечения шириной 36 ± 1 мм, длиной не менее 135 мм и толщиной 3...5 мм. Применяется цилиндрический стальной (сталь Р6М5) стержень диаметром 5,99 ± 0,02 мм. Испытание проводят до деформации 5 %
или момента разрушения, если оно наступает ранее. В процессе испытания записывают диаграмму деформирования Р = Р(5) «напряжение смятия — деформация смятия 5 (овализация отверстия)». Диаграмма деформирования является характеристикой материала в составе конструкции и характеризует сопротивляемость материала воздействию нагрузок в зоне концентраторов напряжений, главным образом соединительных элементов конструкций.
Цель исследования — разработка численного алгоритма решения задачи смятия пластин на основе трехмерного конечно-элементного анализа напряженно-деформированного состояния (НДС) пластины вплоть до разрушения.
Рис. 1. Схема проведения испытания и вид образца при испытании на смятие
Алгоритм основан на трехмерном конечно-элементном решении задачи упруго-пластического деформирования пластин с отверстием при смятии и предназначен для возможности проведения численных экспериментов. Используется модель малых упругопластических деформаций Ильюшина [1]. Решение задачи овализации отверстия пластины при растяжении рассматривается в рамках трехмерной постановки задачи теории пластичности. С помощью этой теории вычисляют поля напряжений, деформаций и перемещений в пластине, а также строят диаграмму деформирования.
Математическая задача деформирования упругопластических пластин при смятии. Рассмотрим постановку трехмерной задачи деформационной теории пластичности при активном нагруже-нии [2]:
=0 в у ;
ау = ) в V
Ц = 2(Щ, у + у) в у; (1)
Щ = Щ на 2И1; Щ = 0 та Ъи2;
ОуПу = 0 на ,
где а у, гк1 — декартовы компоненты тензоров напряжений и деформаций; щ — компоненты вектора перемещений; Г у ) — нелинейная тензорная функция, описывающая модель пластичности; щ у —
е
частные производные; Щ — заданные перемещения на части поверхности тела; часть поверхности ЕЩ 2 полагаем жесткозащемленной, на ней заданы нулевые перемещения. Для плоской пластины, представляющей собой тонкий параллелепипед с отверстием, ребра которого ориентированы по осям координат, эти перемещения были заданы на половине поверхности отверстия пластины и направлены по
оси 0%1 , поэтому
ие = и (Г) = су ; и2 = мэ = 0. (2)
Задаваемое перемещение и (^) изменяется по линейному закону во времени (постоянная скорость нагружения Сд в «жестком» режиме). После решения задачи (1) для каждого момента времени I вычисляем силу, действующую на полуповерхности отверстия,
Рх{г) = |ои щёЪ 0, (3)
Е0
где 1 — полуповерхность отверстия, к которой приложено нагружающее устройство (стержень). Далее определяем среднее напряжение ап и деформацию овализации отверстия §1 в пластине:
5„ =1. Р{,); 5, = М , (4)
Е0 и0
где Б0 — начальный диаметр отверстия в пластине; ) — текущее
значение диаметра отверстия в направлении оси растяжения.
Материал пластины полагаем изотропным, подчиняющимся деформационной теории малых упругопластических деформаций Иль-
юшина для активного нагружения [1]. В этом случае функция ¥ . (в^) имеет следующий вид:
а. = Ру (вм) = К - 3 0[1 - ш (Г2 (в))]вкк5. + 20[1 - ш (У2 (в))]в., (5)
где -^2(в) = 11 -в22)2 + (в 22-взз)2 + (взз -вП)2 + 6(в^2 +в1з2 +в2з2) —
второй инвариант (интенсивность) [з] тензора деформации; К — модуль объемного сжатия материала; G — модуль сдвига; ш — функция пластичности Ильюшина, которую представим в степенном виде [2]:
ю(*2(в)) =
0, ^;
1 -(-^У-"; >вт, (6)
Ч(в) ' ^ т'
где п — показатель упрочнения; вт — деформация начала текучести
материала. Такой выбор функции пластичности характерен для материалов со степенным упрочнением диаграммы деформирования, например для некоторых марок алюминиевых сплавов и сталей. Вводя тензор четвертого ранга
С(У2(в))т = (К - з 0(1 -ш(72(в))))5й 5. + 0(1 -ш(72(в)))(5гк 5 . + 5, 5 .к), (7)
определяющее соотношение (5), записываем в следующем псевдолинейном виде:
= С (в)) ук1 вк1. (8)
для изотропного материала: а. = (К — 0)вкк5. + 20в. .
Это соотношение в частном случае ш = 0 переходит в закон Гука
2 з
Численный метод решения задачи деформирования при смятии. Для решения задачи теории пластичности (1)-(8) применим один из вариантов метода упругих решений [1, 4], согласно которому определяющее соотношение (8) линеаризуем с помощью итерационного алгоритма
а .) = С ^(в(т-1))) 1]к1 вк1 >, (9)
где а(т) и вкт) — тензоры напряжений и деформаций на т-м шаге итерации.
В качестве критерия выбора числа необходимых итераций выбираем условие достижения заданного значения А — относительного отклонения напряжений
И°(Г}(**)-°Г-1)( **)
(т-1),
тах
хк
и
г ,3
_(т-1)
( хк )|
< А.
(10)
С помощью итерационного процесса (9) задачу (1) можно представить в виде последовательности линейных задач
т(т) -1 г] ,3
о3. ■ 3 = 0 в V;
о
(т) _
= С (72(е(т ~-1))ум
,(т) '-к!
в V ^Е;
(11)
3 _ к«3 + Ц]?) в V;
и(т) _ Ц на Е^; иг(т) _ 0 на Еи2; о^пу _ 0 на Еа.
Для численного решения каждой задачи (11) на каждой итерации был применен метод конечных элементов (МКЭ). Вариационную формулировку задачи (11) можно записать стандартным образом:
|5в(т)То(т)(IV _ | 5и(т)Т5(т)(Е,
(12)
дV
куда введены координатные столбцы для компонент вектора перемещений, тензоров напряжений и деформаций:
и(т) _ [ц(т),и2т),и3т)]т;
о(т) _ [о(т) о(т) о(т) о(т) о(т) о(т)]т.
а -[а11 ,а22 ,а33 ,°12 ,а23 ,а31 ] ;
о(т) _ [о(т) о(т) о(т) 2р(т) 2о(т) 2о(т)]т
11 22 33 12 23 13
а 5(т)=[^1(т), 52(т), 53(т)]т — координатный столбец вектора поверхностных сил. Запишем определяющие соотношения и соотношения Коши в системе (11) в матричном виде: о(т) _ С(т)в(т); в(т) _ Би(т),
где С(т) — матрица размером [6*6], составленная из компонент тензора С (У2(г( т~1)))ук1; Б — матрица размером [3*6] дифференциальных операторов [5, 6].
Тогда вариационное уравнение (12) примет следующий вид:
| (П8и (т))т С (т)Ви(т)ёУ = |би(т)т 5 (т)ёЪ. (13)
V ж
Разбивая всю область У на конечные элементы, аппроксимируем перемещения и(т) в каждом КЭ в следующем виде: и(т) = Ф( хк )д(т), где Ф(хк) — матрица функций формы размером [3*п] ( п = 3т, где т — число узлов конечного элемента); q(т) — перемещения в узлах КЭ. Для деформаций и напряжений имеем следующие аппроксимации: а(т) = с(т)В^т); в(т) = £>Ф(хк ^(т}; В = ЯФ(х,). Тогда вариационное уравнение (13) приводим к системе линейных алгебраических уравнений (СЛАУ) для каждого КЭ: К(ет)q(т) = /е, где К(ет) =
= УкВТС(т)В; /е = {Ф(х,)Т 5(т)ёЪ. Путем стандартной процедуры
Ук
ансамблирования СЛАУ собираются в глобальную систему К(т)q(m) = у . Для решения такой СЛАУ большой размерности применяют метод сопряженных градиентов, который в совокупности с разреженным методом хранения СЛАУ [7] обеспечивает хорошую скорость решения.
В рамках проведенного исследования применены конечные элементы в виде тетраэдра с четырьмя узлами (т = 4, п = 12) в вершинах, которые обеспечивают линейную аппроксимацию перемещений и(т) и приводят к постоянным напряжениям а( т) и деформациям в(т) в каждом КЭ. Программная реализация разработанного алгоритма осуществлялась на базе программного обеспечения Научно-образовательного центра «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» МГТУ им. Н.Э. Баумана.
Экспериментальные диаграммы деформирования алюминиевых сплавов. Тестирование разработанного алгоритма численного моделирования осуществлялось путем сравнения расчетных и экспериментальных данных в задаче смятия. На рис. 2 представлены экспериментальные диаграммы деформирования алюминиевого сплава 1163 при растяжении и сжатии, использованного в качестве тестового материала.
Эти диаграммы были использованы для определения констант пластичности в рамках степенной модели (6), (7) Е, V, вт, К, п . Значения этих параметров для сплава 1163: Е = 70 ГПа, К = 63 ГПа, О = 26,3 ГПа, V = 0,33, вт = 0,004, п = 0,31.
ои,МПа 450
360
270
180
90
0 4,4 8,8 13,2 17,6 8П,°Ь
Рис. 2. Диаграммы деформирования алюминиевого сплава 1163, полученные на разных образцах при сжатии (С) и растяжении (Р)
Результаты численного моделирования НДС в задаче овали-зации. Для численного моделирования НДС в задаче овализации при смятии с помощью программы SolidWorks была создана геометрия пластины с отверстием, затем с помощью программы NetGEN была сгенерирована конечно-элементная сетка для решения данной задачи (рис. 3).
Рис. 3. Конечно-элементная сетка, построенная для решения задачи овализации
Параметры и значения сетки:
Используемый тип КЭ..............................Четырехузловой тетраэдр
Число КЭ в сетке....................................41 044
Число узлов в сетке................................. 25 088
Точность определения решения................. 0,00001
Максимальное число итераций в методе
решения СЛАУ..................................... 250
Размер шага по методу упругих решений ... 0,0001
Число точек расчета................................30
Среднее число итераций для вычисления одной точки...........................................25
На рис. 4-6 представлены некоторые из результатов расчетов: поля нормальных напряжений оп и о22, касательного напряжения
о12 в пластине в задаче овализации. Максимальные значения всех
напряжений достигаются непосредственно на поверхности отверстия пластины, в той ее части, где эта поверхность контактирует с растягивающим устройством. От отверстия в сторону передней свободной поверхности (на рис. 4 справа) направлены линии — условные границы максимальных напряжений в пластине. Особенно заметны эти линии на рис. 4 для касательных напряжений.
Определим параметр повреждаемости упругопластического материала
z - Ml) h
( \ о
V0т у
+YHo) h
(
о ос
л
+(^0 о + Bi о2 + B2 Y2 (о))
1 - h
V у
W ( 1 - h
V° у у
с уу
где коэффициенты
1
B0-■
B = °c-0т
0с0т
3о2
°с° т
B2 -
3о2
Y2 (0) ((ОП -022)2 + (0п -033)2 + (022 -°33)2
-б(о:
12
23
"013))
интенсивность тензора напряжений; о = ой — первый инвариант тензора напряжений; о т, о ^, ос — пределы прочности материала при одноосном растяжении, сдвиге и сжатии соответственно; И (у) — функция Хевисайда ( И (у) = 1 при у > 0, и И (у) = 0 при у < 0). Разрушение
материала начинается, когда параметр повреждаемости г = 1.
Такой параметр повреждаемости соответствует критерию прочности Писаренко — Лебедева [5, 6]. При о >от и о < -ос параметр 2
1
s
s
представляет собой отношение интенсивности напряжений к соответствующему пределу прочности при растяжении или сжатии, что соответствует критерию прочности Мизеса. А при -ас < а < ат параметр повреждаемости соответствует критерию прочности Ягна.
И 189 881 (-97.9065 ., 189,881)
173.893 157,905
125.928 Н 109.94 Н 93,9518 Н 77.9636
1 61,9754 45.9873 29.9991 14,0109 Н -1,97733 Н -17,9655 Н -33.9537
1 -49.9419 -65.9301 -81.9183 -97,9065
ь
II
7
I 5 4
I 2
I
I
¿Г.
Рис. 4. Поля напряжений а11 и а12, МПа, в пластине в задаче овализации
Рис. 5. Поля поперечного нормального напряжения ст 22, МПа, и продольного перемещения щ, м, в пластине в задаче овализации
Рис. 6. Поле параметра повреждаемости материала в пластине при двух значениях растягивающего перемещения и (/)
На рис. 6 показано поле параметра повреждаемости материала в пластине при двух значениях растягивающего перемещения и(^), максимум этого параметра достигается на поверхности контакта растягивающего устройства и отверстия пластины, в передней ее точке.
В этой точке и начинается разрушение материала, когда параметр повреждаемости г = 1.
На рис. 7 показаны экспериментальная и расчетная диаграммы деформирования ап ^ 51 алюминиевого сплава 1163 при смятии.
Видно, что совпадение экспериментальных и расчетных результатов достаточно хорошее, что позволяет говорить об адекватности разработанной вычислительной модели.
ап,МПа 500
400 300 200 100
0 0Л 0,2 0,3 0,4 0,5 8!,%
Рис. 7. Экспериментальная (Э) и расчетная (Р) диаграммы деформирования алюминиевого сплава 1163 при смятии
э чр
Разработанный алгоритм может быть применен для расчета НДС тонких пластин [8-14], содержащих конструктивные отверстия, соединительные элементы и другие неоднородности конструкции.
Выводы. Разработан численный конечно-элементный алгоритм решения трехмерной задачи овализации отверстия при смятии образца для упругопластических материалов, деформирование которых подчиняется теории малых упругопластических деформаций Ильюшина.
Численным путем получены поля нормальных напряжений а33
и а22 и касательного напряжения а23 в пластине из алюминиевого
сплава 1163, а также поле параметра повреждаемости материала в пластине. Установлено, что максимум этого параметра достигается на поверхности контакта растягивающего устройства и отверстия пластины, в передней ее точке. В этой точке и начинается разрушение материала, когда параметр повреждаемости достигает единицы.
Численным путем получены диаграммы деформирования при смятии образца с отверстием из алюминиевого сплава 1163. Показано, что имеет место достаточно хорошее совпадение эксперимен-
тальных и расчетных диаграмм овализации, что позволяет говорить об адекватности разработанной вычислительной модели.
ЛИТЕРАТУРА
[1] Ильюшин А.А. Механика сплошной среды. Москва, УРСС, 2014, 320 с.
[2] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 4: Основы механики твердого тела. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2013, 624 с.
[3] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 1: Тензорный анализ. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2011, 463 с.
[4] Биргер И. А., Мавлютов Р. Р. Сопротивление материалов. Москва, Наука,
1986, 560 с.
[5] Димитриенко Ю.И., Яковлев Н.О., Ерасов В.С., Федонюк Н.Н., Сборщиков С.В., Губарева Е.А., Крылов В. Д., Григорьев М.М., Прозоровский А.А. Разработка многослойного полимерного композиционного материала с дискретным конструктивно-ортотропным заполнителем. Композиты и наноструктуры, 2014, т. 6, № 1, с. 32-48.
[6] Димитриенко Ю.И., Федонюк Н.Н., Губарева Е.А., Сборщиков С.В., Прозоровский А.А. Многомасштабное конечно-элементное моделирование трехслойных сотовых композитных конструкций. Наука и образование: электронное научно-техническое издание, 2014, № 10.
ао1: 10.7463/1014.0730105
[7] Димитриенко Ю.И., Соколов А.П. Многомасштабное моделирование упругих композиционных материалов. Математическое моделирование, 2012, т. 24, № 5, с. 3-20.
[8] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Асимптотическая теория конструктивно-ортотропных пластин с двухпериодической структурой. Математическое моделирование и численные методы, 2014, № 1, с. 36-57.
[9] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Конечно-элементное
моделирование эффективных вязкоупругих свойств однонаправленных композиционных материалов. Математическое моделирование и численные методы, 2014, № 2, с. 28-48.
[10] Димитриенко Ю.И., Соколов А.П. Разработка системы автоматизированного вычисления эффективных упругих характеристик композитов. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2008, № 2, с. 57-67.
[11] Димитриенко Ю.И., Соколов А.П. Система автоматизированного прогнозирования свойств композиционных материалов. Информационные технологии, 2008, № 8, с. 31-38.
[12] Димитриенко Ю.И., Сборщиков С.В., Соколов А.П., Садовничий Д.Н., Гафаров Б.Р. Численное и экспериментальное моделирование прочностных характеристик сферопластиков. Композиты и наноструктуры, 2013, № 3, с. 35-51.
[13] Димитриенко Ю.И., Губарева Е.А., Юрин Ю.В. Асимптотическая теория термоползучести многослойных тонких пластин. Математическое моделирование и численные методы, 2014, № 4, с. 18-36.
[14] Димитриенко Ю.И., Губарева Е.А., Яковлев Д.О. Асимптотическая теория вязкоупругости многослойных тонких композитных пластин. Наука и образование: электронное научно-техническое издание, 2014, № 10. ао1: 10.7463/1014.0730105
Статья поступила в редакцию 31.01.2015
Ссылку на эту статью просим оформлять следующим образом: Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В., Ерасов В.С., Яковлев Н.О. Численное моделирование и экспериментальное исследование деформирования уп-ругопластических пластин при смятии. Математическое моделирование и численные методы, 2015, № 1, с. 67-82.
Димитриенко Юрий Иванович родился в 1962 г., окончил МГУ им. М.В. Ломоносова в 1984 г. Д-р физ.-мат. наук, профессор, директор Научно-образовательного центра «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» МГТУ им. Н.Э. Баумана; заведующий кафедрой «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 300 научных работ в области механики сплошных сред, вычислительной механики, механики и термомеханики композитов, математического моделирования в науке о материалах, вычислительной газодинамики. e-mail: dimit.bmstu@gmail.com
Губарева Елена Александровна родилась в 1982 г., окончила МГУ им. М.В. Ломоносова в 2004 г. Канд. физ.-мат. наук, доцент кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор 30 научных работ в области механики сплошных сред, механики контактного взаимодействия, математического моделирования, механики композитов. e-mail: gubareva_ea@pochta.ru
Сборщиков Сергей Васильевич родился в 1989 г., окончил МГТУ им. Н.Э. Баумана в 2012 г. Аспирант кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор 12 научных работ в области вычислительной механики композитов. e-mail: servasbor@gmail.com
Ерасов Владимир Сергеевич родился в 1950 г., окончил МИФИ в 1974 г. Канд. техн. наук, начальник лаборатории «Прочность и надежность материалов воздушного судна» ФГУП «ВИАМ» ГНЦ РФ. Автор более 105 научных работ в области механики материалов и физики прочности.
Яковлев Николай Олегович родился в 1984 г., окончил МИСиС в 2007 г. Канд. техн. наук, научный сотрудник лаборатории «Прочность и надежность материалов воздушного судна» ФГУП «ВИАМ» ГНЦ РФ. Автор более 20 научных работ в области прочности и надежности материалов, методов механических испытаний полимерных, полимерных композиционных и металлических материалов.
Computational modeling and experimental investigation of elastic-plastic plates deforming under crushing
© Yu.I. Dimitrienko1, E.A. Gubareva1, S.V. Sborschikov1, V.S. Erasov2, NO. Yakovlev2
1 Bauman Moscow State Technical University, Moscow, 105005, Russia
2 Federal State Unitary Enterprise "All-Russian Scientific Research Institute of Aviation Materials", Moscow, 105005, Russia
The article presents a suggested method of numerical finite-element solving the 'hole ovalization' problem. This method can be applied for experimental development of advanced aviation materials with the aim of determining structure element resistance against deforming with stress concentrators, mainly, connectors. The method is based on three-dimensional finite element solution of the problem of elastoplastic deformation of plates with a hole under crushing. It is appropriate for reduction of experimental studies
and replacing them by the numerical experiments. The Ilyushin model of small elasto-plastic deformations has been used. The results of numerical simulation of a three-dimensional stress-strain state of elastoplastic plates under crushing are presented as well as results of experimental investigations of deforming plates of Al-alloy 163. It is shown that the results of numerical and experimental modeling for deforming plates under crushing agree quite well.
Keywords:' hole ovalization' problem, plate crumpling, numerical modeling, finite element method, elasticity, plasticity, deforming diagram.
REFERENCES
[1] Ilyushin A.A. Mekhanika sploshnoy sredy [Continuum Mechanics]. Moscow, URSS Publ., 2014, 320 p.
[2] Dimitrienko Yu.I. Mekhanika sploshnoi sredy [Continuum Mechanics]. In 4 vols. Vol. 4. Osnovy mekhaniki tverdogo tela [Fundamentals of Solid Mechanics]. Moscow, BMSTU Publ., 2013, 624 p.
[3] Dimitrienko Yu.I. Mekhanika sploshnoi sredy [Continuum Mechanics]. In 4 vols. Vol. 1. Tenzornyy analiz [Tensor analysis]. Moscow, BMSTU Publ., 367 p.
[4] Birger I.A., Mavlyutov R.R. Soprotivlenie materialov [ Strength of Materials]. Moscow, Nauka Publ., 1986, 560 c.
[5] Dimitrienko Yu.I., Yakovlev N.O., Erasov V.S., Fedonuk N.N., Sborschikov S.V., Gubareva E.A., Krylov V.D., Grigoryev M.M., Prozorovskiy A.A. Kompozity i nanostruktury - Composites andNanostructures, 2014, no. 1, vol. 6, pp. 32-48.
[6] Dimitrienko Yu.I, Fedonuk N.N., Gubareva E.A., Sborschikov S.V, Prozorovsky A.A. Nauka i obrazovanie. Elektronnoe nauchno-tekhnicheskoe izdanie — Science and Education. Electronic Scientific and Technical Journal, 2014, no. 10. doi: 10.7464/1014.0730105
[7] Dimitrienko Yu.I., Sokolov A.P. Matematicheskoe modelirovanie -Mathematical Modeling, 2012, no. 5, vol. 24, pp. 3-20.
[8] Dimitrienko Yu.I., Gubareva E.A., Sborschikov S.V. Matematicheskoe modelirovanie i chislennye menody - Mathematical modeling and Numerical Methods, 2014, no. 1, pp. 36-57.
[9] Dimitrienko Y.I., Gubareva E.A., Sborschikov S.V. Matematicheskoe modelirovanie i chislennye menody - Mathematical modeling and Numerical Methods, 2014, no. 2, pp. 28-48.
[10] Dimitrienko Yu.I., Sokolov A.P. Vestnic MGTU im. N.E. Baumana. Seria Estestvennye nauki - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2008, no. 2, pp. 57-67.
[11] Dimitrienko Yu.I., Sokolov A.P. Informatsionnye Tekhnologii - Information Technologies, 2008, no. 8, pp. 31-38.
[12] Dimitrienko Yu.I., Sborschikov S.V., Sokolov A.P., Sadovnichiy D.N., Gafarov B.R. Kompozity i nanostruktury - Composites and Nanostructures, 2013, no. 3, pp. 35-51.
[13] Dimitrienko Yu.I., Gubareva E.A., Yurin Yu.V. Matematicheskoe modelirovanie i chislennye menody - Mathematical modeling and Numerical Methods, 2014, no. 4, pp.28-48.
[14] Dimitrienko Yu.I., Gubareva E.A., Yakovlev D.O. Nauka i obrazovanie. Elektronnoe nauchno-tekhnicheskoe izdanie — Science and Education. Electronic Scientific and Technical Journal, 2014, no. 10. doi: 10.7463/1014.0730105.
Dimitrienko Yu. I. (b.1962) graduated from Lomonosov Moscow State University in 1984. Dr. Sci. (Phys. & Math.), Professor, Head of the Department of Computational Mathematics and Mathematical Physics, Director of Scientific-Educational Center of Supercomputer Engineering Modeling and Program Software Development at Bauman Moscow State Technical University. Member of the Russian Academy of Engineering Science. The author of more than 300 publications in the field of computational mechanics, gasdynamics, thermomechanics of composite materials, mathematical simulations in material science. e-mail: dimit.bmstu@gmail.com.
Gubareva E.A. (b.1982) graduated from Lomonosov Moscow State University in 2004. Candidate of sciences (Phys. & Math.), associated professor at the Department of Computational Mathematics and Mathematical Physics, at Bauman Moscow State Technical University. The author of more than 30 publications in the field of continuum mechanics, trybology, mathematical modeling and composite mechanics. e-mail: gubareva_ea@pochta.ru.
Sborschikov S.V. (b.1989) graduated from Bauman Moscow State Technical University in 2012. A Ph.D. student at the Department of Computational Mathematics and Mathematical Physics, at Bauman Moscow State Technical University. The author of 12 publications in the field of continuum mechanics of composites. e-mail: servasbor@gmail.com.
Erasov V.S. (b.1950) graduated from Moscow Engineering Physics Institute in 1974. Candidate of Engineering Sciences, head of the laboratory of "Durability and reliability of the aircraft materials" at the Federal State Unitary Enterprise "All-Russian Scientific Research Institute of Aviation Materials". The author of more than 105 publications in the field of mechanics of materials and physics of strength.
Yakovlev N.O. (b.1984) graduated from Moscow Institute of Steels and Alloys in 2007. Candidate of Engineering Sciences, researcher at the laboratory of "Durability and reliability of the aircraft materials" at the Federal State Unitary Enterprise "All-Russian Scientific Research Institute of Aviation Materials". The author of more than 20 publications in the field of durability and reliability of the materials, techniques of mechanical testing polymer, polymer composite and metallic materials.