Моделирование напряженно-деформированного состояния композиционных оболочек с дефектами

УДК 539.3

БОТ 10.18698/2308-6033-2016-11-1552

Моделирование напряженно-деформированного состояния композиционных оболочек с дефектами

© Ю.В. Захарова, Л.Г. Лохматова МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

Рассмотрена актуальная задача моделирования напряженно-деформированного состояния композиционных оболочек с дефектами. Данная тема является перспективной, так как в настоящее время оболочки из композиционных материалов часто выступают в качестве элементов различных силовых конструкций, а появления дефектов практически невозможно избежать даже на стадии изготовления. Предложена математическая модель деформирования цилиндрических оболочек из композиционных материалов с дефектами, основанная на гипотезе Тимошенко. Рассмотрены дефекты типа непроклея. С помощью данной модели описан алгоритм решения задачи о напряженно-деформированном состоянии цилиндрической оболочки с дефектами. Для численного решения поставленной задачи использовался метод конечных элементов. Представленные результаты численного моделирования позволили установить закономерности влияния локальных несовершенств структуры на напряженно-деформированное состояние цилиндрической оболочки из композиционных материалов с дефектами типа непроклея.

Ключевые слова: напряженно-деформированное состояние, цилиндрическая оболочка, композиционные материалы, модель оболочек Тимошенко, непроклей.

Введение. В современных силовых конструкциях находят широкое применение волокнистые композиционные материалы, представляющие собой твердое полимерное связующее, армированное высокопрочными и высокожесткими волокнами [1, 2]. Элементы конструкций (стержни, пластины, оболочки) из подобных материалов имеют требуемые механические свойства при высокой удельной прочности и жесткости. Такие элементы конструкций применяются в ракетно-космической технике, авиационной промышленности, автомобилестроении и в других отраслях [3-6]. Изучение различных свойств конструкций, изготовленных из подобных материалов, является перспективным направлением.

Кроме того, изготовленные из композиционных материалов изделия существенно легче металлических аналогов, вследствие чего в последнее время сформировалась тенденция увеличения габаритов таких конструкций. Увеличение размеров изделия и его составных элементов непосредственно связано с ростом вероятности появления в них технологических дефектов [7-10], существенно влияющих на прочностные свойства всего изделия.

В настоящее время для решения задач о напряженно-деформируемом состоянии (НДС) изделия используется метод конечных элементов (МКЭ) [11-16]. В данной работе с помощью МКЭ решена задача моделирования НДС цилиндрической оболочки из композиционного материала с дефектами, рассматриваемой в рамках теории Тимошенко [17, 18], под действием равномерного внутреннего давления.

Математическая постановка задачи моделирования НДС композиционной оболочки с дефектами. Рассмотрим цилиндрическую оболочку из композиционного материала, находящегося под действием равномерного внутреннего давления р (рис. 1). В качестве кинематической гипотезы примем гипотезу Тимошенко для многослойного пакета [18, 19].

Рис. 1. Геометрия цилиндрической оболочки

В качестве координатной поверхности (2 = 0) примем среднюю поверхность оболочки. Тогда распределение по толщине оболочки будем вычислять по формулам

Щ,=Ца + , а = Х, у, Ы = Ж,

где иХ, иу — тангенциальные перемещения; Ж — прогибы оболочки; ФХ, Фу — углы поворота сечений.

С учетом принятой гипотезы и геометрии рассматриваемой оболочки деформационные соотношения примут вид

еХ ех + Х , £ у еу + у, £Ху еху +

Ху

где

диХ

-— е

дХ ' еу

дЦ Ж „ -- +—; 2е

ду г

диХ ди

у .

Ху

+; х.

ду дХ

дХ

Ху

дФ

У .

ду

2х

Ху

дФ Х дФ у = Цу

ду + дХ ; Ф дХ ; Фу ~1у2 дХ + Я •

В соотношениях (3) ех, еу, еху — деформации растяжения-сжатия и сдвига срединной поверхности оболочки; хх, ху, хху —

деформации изгиба и кручения.

Для краткости записи деформационные соотношения (3) запишем в матричном виде:

е = Lu .

Здесь е = \ех еу 2е^ х х X у 2Х ху Ъ * Ъ у ]т, и = [Ж Ъ х Ъ у их и у ]т,

где L — матрица дифференциального оператора, полученная из соотношений (3).

Определяющие соотношения для модели оболочек Тимошенко имеют вид

^аа = Саа еа + Сарер + X3

°ху = 2С66еху + 2Х С66Хху, Ту = 2С44Ъу, Тх = 2С55Ъх.

где сар — компоненты тензора напряжений.

Запишем соотношения упругости для внутренних погонных силовых факторов:

Тх = С11ехх + С12еуу; Ту = С22еуу + С21ехх; Тху = 2С66еху; ° = 2С55 У х; 0у = 2С44 У у; Мхх = А1Ххх + А2Xуу; Муу = В22Xуу + А1Ххх; Мху = 2 ВббХ ху,

где Тар — погонные усилия; Мар — моменты; Qа — перерезывающие силы.

Запишем определяющие соотношения в матричном виде:

Т11 ] С11 С12 0 0 0 0 0 0 ехх

Т22 С21 С22 0 0 0 0 0 0 еуу

Т12 0 0 С66 0 0 0 0 0 2еху

М1 0 0 0 А1 А2 0 0 0 х х

М2 0 0 0 А1 О22 0 0 0 х у

М12 0 0 0 0 0 Об6 0 0 2х ху

01 0 0 0 0 0 0 2С55 0 Ъ х

Q2 ] 0 0 0 0 0 0 0 2С44 Ъ у

или

Т = De,

- h5 где Caß = hCan — мембранные жесткости; Daß =— Caß — изгибные

12

жесткости. Вычислим Caß по следующим формулам:

E1h Е2h

C11 =

1 -Hl2^21

C22 =

Сбб = G12h; C55 = Gi

712

55

'13'

1 — ^12^21

C44 = G23^

C12 = ^21C11;

44

„3

D66 =

D11 = C111 11 12 11

D11 = C11'.

22 12 22^

G12 h 12

12 = ^21D11;

Du =

где G12, 023, G13 — модули поперечного сдвига.

Вариационная постановка задачи статики. Используем вариационную постановку задачи в перемещениях на основе принципа Ла-гранжа [20-25]. Это позволит не выписывать явно уравнения совместности деформаций, которые выполняются тождественно.

Разность энергии деформаций и работы внешних сил с учетом выражения деформаций через перемещения (4) дает минимизируемый функционал потенциальной энергии, в котором учтены статические граничные условия (в данной работе в качестве граничных условий принято жесткое закрепление торцев оболочки, что приводит к соотношениям иа = 0; Ж = 0; Фа = 0, а = 1,2):

П =1Л етБейУ -Д итрй$р = - Д

Рис. 2. Разбиение цилиндрической оболочки на треугольные конечные элементы

Минимум функционала необходимо отыскивать в пространстве кинематически допустимых полей. Окончательно, с учетом принятых гипотез при кинематических условиях, определяемых схемой закрепления, задача сводится к минимизации функционала вида

Птп ^ min

Ее приближенное решение может быть получено путем дискретизации модели.

Квадратичный двухмерный конечный элемент. Для численного решения задачи разобьем оболочку на конечные элементы (рис. 2).

В качестве конечного элемента принят шестиузловой треугольник. Его конечно-элементная аппроксимация представлена на рис. 3.

Для конечного элемента с помощью естественных координат были вычислены функции формы:

N1 =(2/4 -1)4; N4 = 4L2Lз;

N2 = 4/1/2; N5 =(2/з -1) /;

N3 =( 2/2 -1) /2; N6 = 4/1/3.

Рис. 3. Расположение узлов в треугольном элементе с квадратичной аппроксимацией

Связь деформаций координатной поверхности с узловыми перемещениями. Выразим перемещения в произвольной точке через перемещения узлов конечного элемента (три линейных и два угла поворота нормали):

и

их

иу Ж

Ъ х

Ъ у

[ N ы2 ы3 ыА ы5 ы6 ]

'51 §2

53

54

55

56

где 51 =[и1; ,и21 ,Wi, Ъ х, Ъ у1 ] — перемещения I -го узла; N — блочная матрица функций форм. При этом

N =[ N1 N2 N3 N.4 N5 N6 ];

" N 0 0 0 0 "

0 N 0 0 0

N = 0 0 N 0 0

0 0 0 N 0

0 0 0 0 N ]

С учетом перемещений (12) связь деформаций координатной поверхности с узловыми перемещениями имеет вид

е = В5,

где В — матрица производных от функций форм, В = LN.

Подставляя полученные аппроксимации в функционал (9), получаем

П

2

где К(е) = | BтDBdS — матрица жесткости элемента; P(e) = | NтpdS —

^

вектор узловых сил.

Суммируя по элементам матрицу и вектор, получаем глобальную матрицу жесткости и вектор. Минимизируя потенциальную энергию, переходим к системе линейных алгебраических уравнений:

К 5 = 0,

решая которую, находим перемещения узлов.

Определение деформаций и напряжений. Для определения деформаций и напряжений вычислим перемещения в точке М, соответствующей середине конечного элемента:

и

Затем вычисляем обобщенные деформации в точке М, используя формулы (3):

е

При численной реализации матрицу В(М) необходимо сформировать при обработке конечного элемента.

Для вычисления усилий в центре конечного элемента необходимо умножить деформации на матрицу коэффициентов приведенных жесткостных характеристик (7):

Т

Далее можно вычислить напряжения, используя формулы (5). Математическое моделирование дефекта оболочки. При

изготовлении оболочек невозможно избежать появления начальных расслоений и непроклеев. Данные виды дефектов моделируются изменением упругих характеристик материала, определяемых формулами (8).

В представленной работе исследовали дефекты типа непроклея прямоугольной формы, характеризуемые уменьшением жесткостных характеристик материала в зоне дефекта в 4 и 9 раз [5]. Рассматривались дефекты, которые или примыкают к краю оболочки, или находятся на ее середине (рис. 4).

Рис. 4. Расположение дефектов в оболочке

Результаты численного моделирования НДС оболочки. Для

численной реализации рассмотрели цилиндрическую оболочку длиной 2 м, радиусом 1 м и толщиной 0,0015 м. Композиционный материал, из которого изготовлена оболочка, имеет следующие характеристики: Ех = Е2 = 200 ГПа; Д12 = Ц21 = 0,3; = ^23 = вп = 76,9 ГПа. Оболочка жестко закреплена с обоих концов и находится под равномерным внутренним давлении р = 0,1 МПа.

Для оценки достоверности результатов, полученных с помощью модели, описанной выше, сравним некоторые расчеты с аналитическими. Согласно работе [17], максимальный прогиб срединной поверхности оболочки можно рассчитать по формуле

1 (рМо +0о),

W =

'' тах

2р3 Н

где М0 — изгибающий момент; Q0 =-— — перерезываю-

щая сила.

Коэффициенты р и Н определяли по формулам:

Н

ЕЬ3

в4

4Я2 Н

12

После подстановки в формулу (20) исходных данных и формул (21) получили:

Жтах = 3,2-10 -4м. Результаты, полученные при численном решении этой задачи, представлены на рис. 5.

0,0003 0,0002 0,0001

( \

О 0,5 1,0 1,5 2,0 /, м Рис. 5. Распределение прогибов по длине оболочки

Также был проведен расчет напряжений по безмоментной теории для оболочки с исходными параметрами, но не имеющей закреплений на торцах, в соответствии с работами [17, 18, 22]:

рЯ

'УУ

h

6,7 -107Па.

При численном решении данной задачи получены результаты, представленные на рис. 6.

вуу, МПа

0,0003 0,0002 0,0001

о

0,5 1,0 1,5 2,0 /, м

Рис. 6. Распределение окружных напряжений в оболочке со свободными концами

В соответствии с представленными на рис. 5 и 6 данными в обоих случаях результаты численных решений довольно близки по значениям к результатам, полученным при аналитических решениях. Относительная погрешность не превышает 7 %. Это доказывает, что разработанная модель и алгоритм работают верно.

Также проведено параметрическое исследование зависимости НДС оболочки от размера, физико-механических свойств дефекта и его расположения в оболочке.

На рис. 7-10 представлены результаты численного моделирования НДС цилиндрической оболочки с дефектами.

Исходя из данных рис. 7-10 можно сделать вывод, что при одинаковых параметрах дефекта материал деформируется сильнее, если дефект примыкает к торцу оболочки.

ж, м

0,0025 0,0020 0,0015 0,0010 0,0005

0

, \\.......

1 \ \ 1 \ 1

' /2

\

Оуу, МПа

800 600 400 200

/4

/х / Ч.)/

! » ^2

А

.... /

0,5 1,0 1,5 2,0 I, м а

0 0,5 1,0 1,5 2,0 /,м б

Рис. 7. Зависимость прогибов (а) и окружных напряжений (б) от размера дефекта, примыкающего к торцу оболочки:

1 — без дефекта; 2-4 — Воб/Бдеф равно 784, 392 и 31 соответственно

0,0020 0,0015 0,0010 0,0005

0

¡1-3

1 2

1

/

г .....I >

<*уу> МПа

800 600 400 200

Л

1 уз

\ ,2

1

/

0,5

1,0 1,5 а

2,0 I, м

О

0,5

1,0 1,5 б

2,0 I, м

Рис. 8. Зависимость прогибов (а) и окружных напряжений (б) от характера изменений жесткостных характеристик дефекта, примыкающего к торцу оболочки:

1 — без дефекта; 2, 3 — уменьшение жесткостных характеристик дефекта в 4 и 9 раз соответственно

Ж, м

0,0025 0,0020 0,0015 0,0010 0,0005

Оуу, МПа

А"-- -4

! 1 \ 11 \

11 \ 1 \ г

1 \

\ * \\ % \

/

800 600 400 200

/ Л -

ь \ ¥

ч I I -4

! I

I I

0 0,5 1,0 1,5 2,0 /,м о 0,5 1,0 1,5 2,0 1,м

а б

Рис. 9. Зависимость прогибов (а) и окружных напряжений (б) от размера дефекта, располагающегося по середине оболочки: 1-4 — см. рис. 7

Щм

0,002 0,001

'V

✓ —

Оуу, МПа

800 600 400 200

3 ,,

1 ¡П

/ 4

V/

0

0,5 1,0 1,5 2,0 /, м

0

0,5 1,0 1,5 2,0 I, м

а б

Рис. 10. Зависимость прогибов (а) и окружных напряжений (б) от характера изменений жесткостных параметров дефекта, располагающегося по середине оболочки:

1-3 — см. рис. 8:

Уменьшение жесткостных свойств в зоне дефекта ведет к более быстрому разрушению материала, чем при увеличении размера дефекта. Численное моделирование показало, что уменьшение жесткостных характеристик в 9 раз приводит к существенному увеличению локальных прогибов и окружных напряжений более чем в 8 раз.

Заключение. Разработана математическая модель расчета НДС цилиндрических оболочек типа Тимошенко. Представлен алгоритм поставленной задачи на основе метода конечных элементов и проведено численное моделирование дефектов типа непроклея. Относительная погрешность результатов сравнения численного моделирова-

ния НДС композиционных оболочек, имеющих дефекты, с аналитическим решением не превысила 7 %. Проведено параметрическое исследование влияния размера, физико-механических характеристик и расположения дефекта на НДС оболочки.

Разработанная математическая модель и ее численная реализация могут быть использованы для дальнейших исследований НДС композиционной оболочки с дефектами.

ЛИТЕРАТУРА

[1] Choudhry R.S., Hassan Syed F., Li S., Day R. Damage in Single Lap Joints of Woven Fabric Reinforced Polymeric Composites Subjected to Transverse Impact Loading. International Journal of Impact Engineering, 2015, vol. 80, pp. 76-93.

[2] Димитриенко Ю.И., Соколов А.П., Шпакова Ю.В., Юрин Ю.В. Моделирование поверхностей прочности композитов на основе микроструктурного конечно-элементного анализа. Наука и образование, 2012, № 11. DOI: 10.7463/1112.0496336 (дата обращения 05.07.2016).

[3] Гузь А.Н., Григоренко Я.М., Бабич И.Ю. и др. Механика композитных материалов и элементов конструкций: В 3 т. Т. 2. Механика элементов конструкций. Киев, Наукова думка, 1983, 464 с.

[4] Сыздыков Е.К., Димитриенко Ю.И., Минин В.В. Моделирование термомеханических процессов в композитных оболочках при локальном нагреве излучением. Механика композиционных материалов и конструкций, 2011, т. 17, № 1, с. 71-90.

[5] Каледин В.О., Аульченко С.М., Миткевич А.Б, Решетникова Е.В., Седова Е.А., Шпакова Ю.В. Моделирование статики и динамики оболочечных конструкций из композитных материалов. Москва, Физматлит, 2014, 196 с.

[6] Aul'chenko S.M., Kaledin V. O., Shpakova Yu.V. Characteristic Features of the Oscillation of Shells of Bodies of Revolution Immersed in a Viscous Fluid Flow. Journal of Engineering Physics and Thermophysics, 2008, vol. 81, iss. 5, pp. 941-944.

[7] Laptev N.I., Mordasov V.I., Poilov V.V., Shulepova O.V. Laser Vibrodiagnostics of Starved Spots of Large-sized Multilayered Honeycomb Shells. Polymer Science. Series D, 2010, vol. 3, no. 4, pp. 263-266.

[8] Cawley P. Low frequency NDT techniques for the detection of disbonds and delaminations. NDT & E International, 1992, vol. 25, iss. 2, pp. 100. D0I:10.1016/0963-8695(92)90556-V (дата обращения 05.07.2016).

[9] Ouinas D., Bouiadjra B. Bachir, Achour T., Benderdouche N. Influence of Disbond on Notch Crack Behaviour in Single Bonded Lap Joints. Materials & Design, 2010, vol. 31, iss. 9, pp. 4356-4362.

[10] Каледин В.О., Аникина Ю.В. Численное моделирование статики и устойчивости подкрепленных оболочек с расслоениями. Вестник Томского государственного университета, 2006, № 19, с. 225-232.

[11] Агапов В.П. Метод конечных элементов в статике, динамике и устойчивости пространственных тонкостенных подкрепленных конструкций: учеб. пособие. Москва, Изд-во АСВ, 2000, 152 с.

[12] Голованов А.И., Тюленева О.И., Шигабутдинов А.Ф. Метод конечных элементов в статике и динамике тонкостенных конструкций. Москва, Физматлит, 2006, 392 с.

[13] Зенкевич О. Метод конечных элементов в технике. Москва, Мир, 1975, 541 с.

[14] Рикардс Р.Б. Метод конечных элементов в теории оболочек и пластин. Рига, 1988, 284 с.

[15] Сегерлинд Л. Применение метода конечных элементов. Москва, Мир, 1979, 392 с.

[16] Каледин В.О., Аникина Ю.В., Бурнышева Т.В., Решетникова Е.В. Математическое моделирование статики сетчатой оболочки с учетом концентрации напряжений. Вестник Томского государственного университета, 2006, № 19, с. 233.

[17] Тимошенко С.П., Войновский-Кригер С., Снитко И.К., ред. Пластинки и оболочки. Москва, Наука, 1966, 636 с.

[18] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 4. Основы механики твердого тела. Москва, Изд-во МГТУ им. Н. Э. Баумана, 2013, 624 с.

[19] Кузьмин М.А., Лебедев Д.Л., Попов Б.Г., Данилова В.Л., ред. Расчеты на прочность элементов многослойных композитных конструкций. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2012, 341 с.

[20] Алфутов Н.А., Зиновьев П.А., Попов Б.Г. Расчет многослойных конструкций и оболочек из композиционных материалов. Москва, Машиностроение, 1984, 264 с.

[21] Попов Б.Г. Расчет многослойных конструкций вариационно-матричными методами. Москва, Изд-во МГТУ им. Н.Э. Баумана, 1993, 294 с.

[22] Биргер И.А., Пановко Я.Г. Прочность, устойчивость, колебания: справочник. В 3 т. Т. 2. Москва, Машиностроение, 1968, 464 с.

[23] Димитриенко Ю.И., Шпакова Ю.В., Богданов И.О., Сборщиков С.В. Моделирование процесса многоуровневой фильтрации жидкого связующего в тканевом композите при RTM методе изготовления. Инженерный журнал: наука и инновации. 2015, вып. 12 (48). URL: http://engjournal.ru/catalog/msm/pmcm/1454.html

[24] Dimitrienko Yu.I., Dimitrienko I.D., Sborschikov S.V. Multiscale Hierarchical Modeling of Fiber Reinforced Composites by Asymptotic Homogenization Method. Applied Mathematical Sciences, 2015, vol. 9, no. 145, рр. 7211-7220. URL: http://dx.doi.org/10.12988/ams.2015.510641

[25] Димитриенко Ю.И., Коряков М.Н., Захаров А.А., Строганов А.С. Численное моделирование сопряженных аэрогазодинамических и термомеханических процессов в композитных конструкциях высокоскоростных летательных аппаратов. Математическое моделирование и численные методы, 2014, № 3, с. 3-24.

Статья поступила в редакцию 26.09.2016

Ссылку на эту статью просим оформлять следующим образом: Захарова Ю.В., Лохматова Л.Г. Моделирование напряженно-деформированного состояния композиционных оболочек с дефектами. Инженерный журнал: наука и инновации, 2016, вып. 11. http://dx.doi.org/10.18698/2308-6033-2016-11-1552

Захарова Юлия Владимировна родилась в 1982 г., окончила Кемеровский государственный университет в 2004 г. Канд. техн. наук, доцент кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 50 научных работ в области механики деформируемого твердого тела, механики композиционных материалов и конструкций. e-mail: shpakovayuliya@bmstu.ru

Лохматова Любовь Григорьевна родилась в 1992 г., окончила МГТУ им. Н.Э. Баумана в 2016 г., кафедра «Вычислительная математика и математическая физика». e-mail: lohmatova0408@gmail.com

Simulation of stress-strain state of defected composite shells

© Yu.V. Zakharova, L.G. Lokhmatova Bauman Moscow State Technical University, Moscow, 105005, Russia

The paper examines an important problem of modeling the stress-strain state of the defected composite shells. The issue is a promising one because presently shells made of composite materials often serve as elements of different bearing constructions and it is often impossible to avoid defects even at the manufacturing stage. We propose a mathematical model of deformation of defected cylindrical shells made of composite materials. The model is based on Timoshenko hypothesis. We examined defects of disbond type and with the help of this model we described an algorithm for solving the problem of stress-strain state of defected cylindrical shells. For the numerical solution of the problem we used the finite element method. The presented results of numerical simulations have established patterns of influence of local imperfections in the structure on the stress-strain state of cylindrical shell made of composite materials with defects such as disbonds.

Keywords: stress-strain state, cylindrical shell, composite materials, Timoshenko shell model, disbonds.

REFERENCES

[1] Choudhry R.S., Hassan Syed F., Li S., Day R. Damage in single lap joints of woven fabric reinforced polymeric composites subjected to transverse impact loading. International Journal of Impact Engineering, 2015, vol. 80, pp. 76-93.

[2] DimitrienkoYu.I., Sokolov A.P., Shpakova Yu.V., Yurin Yu.V. Nauka i obra-zovanie — Science and Education, 2012, no. 11. DOI: 10.7463/1112.0496336 (accessed July 5, 2016).

[3] Guz A.N., Grigorenko Ya.M., Babich I.Yu. et al. Mekhanika kompozitnykh ma-terialov i elementov konstruktsiy. V 3 tomakh. Tom 2. Mekhanika elementov konstruktsii [Mechanics of Composite Materials and structural elements. In 3 vols. Vol. 2. Mechanics of structural elements]. Kiev, Naukova dumka Publ., 1983, 464 p.

[4] Dimitrienko Yu.I., Minin V.V., Syzdykov E.K. Mekhanika kompozitsionnykh materialov i konstruktsiy — Mechanics of composite materials and structures, 2011, vol. 17, no. 1, pp.71-90.

[5] Kaledin V.O., Aulchenko S.M., Mitkevich A.B., Reshetnikova E.V., Sedova E.A., Shpakova Yu.V. Modelirovanie statiki i dinamiki obolochechnykh konstruktsii iz kompozitnykh materialov [Simulation of statics and dynamics of shell structures made of composite materials]. Moscow, Fizmatlit Publ., 2014, 196 p.

[6] Aulchenko S.M., Kaledin V.O., Shpakova Yu.V. Characteristic features of the oscillation of shells of bodies of revolution immersed in a viscous fluid flow. Journal of Engineering Physics and Thermophysics, 2008, vol. 81, iss. 5, pp. 941-944.

[7] Laptev N.I., Mordasov V.I., Poilov V.V., Shulepova O.V. Laser vibrodiagnos-tics of starved spots of large-sized multilayered honeycomb shells. Polymer Science. Series D, 2010, vol. 3, no. 4, pp. 263-266.

[8] Cawley P. Low frequency NDT techniques for the detection of disbonds and de-laminations. NDT & E International, 1992, vol. 25, no. 2, pp. 100. DOI: 10.1016/0963-8695(92)90556-V (accessed July 5, 2016).

[9] Ouinas D., Bouiadjra B. Bachir, Achour T., Benderdouche N. Influence of disbond on notch crack behaviour in single bonded lap joints. Materials & Design, 2010, vol. 31, no. 9, pp. 4356-4362.

[10] Kaledin V.O., Anikina Yu.V. Vestnik Tomskogo gosudarstvennogo universiteta — Tomsk State University Journal, 2006, no. 19, pp. 225-232.

[11] Agapov V.P. Metod konechnykh elementov v statike, dinamike i ustoychivosti prostranstvennykh tonkostennykh podkreplennykh konstruktsiy [The finite element method in statics, dynamics and spatial stability of thin-walled reinforced constructions]. Moscow, ASV Publ., 2000, 152 p.

[12] Golovanov A.I., Tyuleneva O.N., Shigabutdinov A.F. Metod konechnyih ele-mentov v statike i dinamike tonkostennyih konstruktsiy [The finite element method in statics and dynamics of thin-walled structures]. Moscow, Fizmatlit Publ., 2006, 392 p.

[13] Zenkevich A. Metod konechnykh elementov v tekhnike [Finite element method in technology]. Moscow, Mir Publ., 1975, 541 p.

[14] Rickards R.B. Metod konechnykh elementov v teorii obolochek i plastin [Finite element method in the theory of shells and plates]. Riga, 1988, 284 p.

[15] Segerlind L.J. Applied Finite Element Analysis. Wiley, New York, 1976. [In Russ.: Segerlind L. Primenenie metoda konechnykh elementov. Moscow, Mir Publ., 1979, 392 p.].

[16] Kaledin V.O., AnikinaYu.V., Burnysheva T.V., Reshetnikova E.V. Vestnik Tomskogo gosudarstvennogo universiteta — Tomsk State University Journal, 2006, no. 19, pp. 233.

[17] Timoshenko S.P., Woinowskiy-Krieger S., Snitko I.K., eds. Plastinki i oboloch-ki [Theory of plates and shells]. Moscow, Nauka Publ., 1966, 636 p.

[18] Dimitrienko Yu.I. Mehanika sploshnoy sredy. V 4 tomakh. Tom 4. Osnovy me-hanik i tverdyh sred [Continuum Mechanics. In 4 vol. Vol. 4. Fundamentals of mechanics of solid media]. Moscow, BMSTU Publ., 2013, 624 p.

[19] Kuzmin M.A., Lebedev D.L., Popov B.G. Raschety na prochnost elementov mnogosloinykh kompozitnykh konstruktsiy [Strength calculations of elements of multilayer composite structures]. Мoscow, BMSTU Publ., 2012, 341 p.

[20] Alfutov N.A., Zinoviev P.A., Popov B.G. Raschet mnogosloinykh konstruktsiy i obolochek iz kompozitsionnykh materialov [Calculation of multilayer structures and shells made of composite materials]. Moscow, Mashinostroenie Publ., 1984, 264 p.

[21] Popov B.G. Raschet mnogosloinykh konstruktsiy variatsionno-matrichnymi metodami [Calculation of multilayer structures by variation-matrix methods]. Moscow, BMSTU Publ., 1993, 294 p.

[22] Birger I.A., Panovko Ya.G. Prochnost, ustoichivost, kolebaniya [Strength, Stability, Vibrations]. Handbook in 3 vol. Vol. 2. Moscow, Mashinostroenie Publ., 1968, 464 p.

[23] Dimitrienko Yu.I., Shpakova Yu.V., Bogdanov I.O., Sborshikov S.V. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2015, iss. 12 (48). Available at: http://engjournal.ru/catalog/msm/pmcm/1454.html

[24] Dimitrienko Yu.I., Dimitrienko I.D., Sborschikov S.V. Multiscale Hierarchical Modeling of Fiber Reinforced Composites by Asymptotic Homogenization Method. Applied Mathematical Sciences, 2015, vol. 9, no. 145, 7211-7220. Available at: http://dx.doi.org/10.12988/ams.2015.510641 http://www.m-hikari.com/ams/ams-2015/ams-145-148-2015/p/dimitrienkoAMS 145-148-2015.pdf

[25] Dimitrienko Yu.I., Koryakov M.N., Zakharov A.A., Stroganov A.S.

Matematicheskoe modelirovanie i chislennye metody — Маthematical Modeling and Computational Methods, 2014, no. 3 (3), pp. 3-24.

Zakharova Yu.V. (b. 1982) graduated from Kemerovo State University in 2004. Cand. Sci. (Tech.), Assoc. Professor of the Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University. Author of over 50 publications in the field of solid mechanics, mechanics of composite materials and structures. e-mail: shpakovayuliya@bmstu.ru, shpakovayulia@gmail.com

Lokhmatova L.G. (b. 1992) graduated from the Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University in 2016. e-mail: lohmatova0408@gmail.com