Теория пластин, основанная на методе асимптотических разложений

УДК 539.3

Теория пластин, основанная на методе асимптотических разложений

© С.В. Шешенин, К. А. Скопцов

МГУ им. М.В. Ломоносова, Москва, 119991, Россия

Приведено сравнение результатов асимптотического анализа поперечного изгиба многослойной пластины под воздействием поверхностной нагрузки с классическими теориями тонких и толстых пластин. Слои пластины полагаются составленными из однородных упругих ортотропных материалов.

Ключевые слова: асимптотический метод, метод осреднения, слоистая пластина, теории тонких и толстых пластин, пластины Кирхгофа и Миндлина — Рейснера.

Одним из способов упрощения трехмерной краевой задачи является понижение размерности, которое может быть проведено, например при учете малости каких-либо параметров, входящих в задачу. Поскольку толщина пластины мала по сравнению с ее другими характерными размерами, пластину можно моделировать плоскими конечными элементами, функции формы для которых получают исходя из решения вспомогательной задачи деформирования тонкого тела. Понижение размерности при этом может быть осуществлено при помощи методов асимптотического анализа [1, 2], заменяющих традиционные гипотезы теории пластин.

Исследование асимптотического поведения уравнений равновесия при помощи рядов специального вида, применяющихся при осреднении уравнений с периодическими быстроосциллирующими коэффициентами [1-3], позволяет при некоторых допущениях о граничных условиях на боковой поверхности построить решение, асимптотически сходящееся к решению задач трехмерной теории упругости. Для случая пластин применение асимптотического разложения в ряд началось с работы [4]. Обзор литературы можно найти в [2]. Анализ первых трех приближений приведен в [5, 6]. Исследование пологих оболочек дано в [7]. Сравнение асимптотической теории пластин с трехмерной теорией можно найти в [8, 9].

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

Цель статьи состоит в том, чтобы показать, каким образом из асимптотического разложения возникает бигармоническое уравнение теории пластин Кирхгофа — Лява [10] и уравнения теорий толстых

пластин [11-12]. Частично этот вопрос уже рассмотрен в [13]. Также в статье описана общая процедура получения локальных функций и уравнений произвольного приближения для изгиба слоистой анизотропной пластины. Отметим, что при выводе уравнений не используются гипотезы, свойственные классическим теориям пластин. Рассмотрение касается только случая линейной теории. Асимптотическому анализу нелинейной проблемы посвящена работа [14], а также работы [15, 16], в которых изучена задача изгиба тонких пластин из физически нелинейного материала.

Постановка задачи. Рассматривается краевая задача линейной теории упругости: уравнение равновесия

в области О = ш[-к /2, к /2] с границей Е = Еп иЕа с граничными условиями

на части границы Еп и щ = 0 на части границы Ел.

Асимптотическое исследование проводится при к ^ 0. Примем, что модули упругости Сф1 зависят только от вертикальной координаты. Удобно считать, что область определения модулей упругости нормирована толщиной пластины, т. е. модули упругости в слое х3 суть Суы (х3/к). Границу области ш, которую занимает тело в плоскости Ох1 х2, обозначим у. Если не сказано иного, будем полагать, что нагружена лишь верхняя граница (Х3 = к /2) нагрузкой э (х1, х2 ) = = э35г-3. Начнем с замечания относительно гипотезы Кирхгофа — Лява.

Гипотеза Кирхгофа — Лява и линейная часть поля перемещений. Гипотеза Кирхгофа — Лява состоит в том, что поле перемещений в тонкой однородной изотропной пластине, находящейся в состоянии изгиба, представляется в виде

где w (х1, Х2) — прогиб пластины. Здесь и далее прописными латинскими буквами будем обозначать индексы, принимающие значения 1, 2 в отличие от индексов, обозначаемых строчными латинскими буквами, принимающими значения 1, 2, 3.

Покажем, что указанные представления суть главные линейные части поля перемещений трехмерной задачи. Не снижая общности,

= С!]к1щк, ¡п] = Э

щ (х1, ^ х3) = - хэ^ (х1, х2); и3 (хЬ х2, х3) = w (хь х2), I = 1, 2,

считаем, что aj (хь 0) = 0 (состояние изгиба без растяжения-сжатия отсчетной поверхности). Из уравнений равновесия следует

a13,3

откуда a13 (x1, x2, x3) = a13 (x1, x2, 0) + O(x2) при x3 ^ 0. Так как aI3 (x1, x2, h/2) = sI = 0, то aI3 (x1, x2, x3) = O (h2) при h ^ 0 и x3/ h = const. В силу ортотропности пластины, как следствие, имеем sI3 (x1, x2, x3) = O(h2). Значит,

ui,3

Разложим uI в ряд по x3, тогда

Ul

Обозначив w (x1, x2) = U3 (x1, x2, 0), с учетом соотношения (5) и ui (x1, x2, 0) = 0 получаем из (6) следующее выражение

ui = -x3U3, i

что в первом приближении соответствует гипотезе Кирхгофа — Лява.

Для вертикальной компоненты поля перемещений можно провести похожие рассуждения. А именно в силу уравнений равновесия имеем на свободной границе

a33,3

откуда a33 (x1, x2, x3) = O (h2) при h ^ 0. Значит, u3 3 (x1, x2, 0) = = O (h2 ), следовательно, получаем

U3

что также соответствует гипотезе Кирхгофа — Лява. Поэтому в работах [5, 6] первые члены асимптотического ряда для представления поля перемещений в трехмерной задаче теории упругости были выбраны в виде (7), (9).

Асимптотическое понижение размерности задачи. Представленное рассмотрение может быть упрощено и обобщено при помощи асимптотического анализа. Будем искать решение задачи (1), (2) в более общем виде, чем в [5, 6]:

щ

^ 0

в котором х = х1, х2, и суммирование проводится по всевозможным непустым мультииндексам а.

Уравнения равновесия (1) предполагаются выполненными при подстановке 2, = х3/к. Таким образом, Е, играет роль «быстрой» координаты, как в стандартном методе осреднения [1, 3].

Подстановка ряда (10) в выражение о у = С^щ, 1 для компонент

тензора напряжений и приведение подобных по степеням параметра к дает асимптотический ряд

Оу = 2 кlа-1PjkаWk,а, (11) Н> 0

в котором

Рук аЬ = СугЬ—гк а + Суг 3 —Гк аЬ и —2к = ^гк.

Здесь и далее штрихом обозначена производная соответствующей функции от Е, по

Дальнейшая процедура понижения размерности состоит в разделении задачи теории упругости на задачу определения локальных функций N к а, зависящих от «быстрой» вертикальной координаты, и задачу определения функций wk, зависящих от «медленных» пла-нарных координат. При этом п-м приближением мы будем называть асимптотический анализ при учете всех слагаемых рядов (10), (11) при И ^ п +1 (номер приближения соответствует максимальной степени к в ряде (11)).

Нахождение локальных функций. При помощи выбора локальных функций —ка можно асимптотически обратить дифференциальное уравнение равновесия (2) в тождество. Подстановка ряда (11) в уравнение равновесия (2) дает формальное равенство

к-1Р'3кл,Ь +2 кН-1

Н> 0

Если потребовать выполнения равенств

р3кЬ = 0 и р3каЬ + Рьк а = 0,

то левая часть (13) станет асимптотически эквивалентна нулю. Эти равенства являются обыкновенными дифференциальными уравнениями второго порядка на функции N\ка, решения которых приобре-

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

Предположим, что найдены все локальные функции Pijka и Nika.

Покажем, как найти все функции PiJkaL и Nik aL следующего приближения. Уравнения (14) позволяют получить все Pi3kaL с точностью до постоянных слагаемых, которые можно фиксировать дополнительным условием. Из соотношений (12) (при J = 3) можно выразить Nik aL:

%

NikaL

0

где M (%) — матрица, обратная матрице чисел Ct3r3 (%). Остается воспользоваться соотношениями (12) при J = 1,2 для нахождения оставшихся функций PiJk aL:

PJkaL = CiJrLNrk a + CJr3Mrs

Как видно из соотношений (15), кроме условий на Pi 3kaL, для однозначности локальных функций требуются условия, позволяющие фиксировать NikaL (0). О них и пойдет речь далее.

Получение двухмерной задачи. Покажем, как можно удовлетворить нагрузке si, распределенной по верхней лицевой поверхности тела x3 = h /2 (в линейной задаче можно без потери общности считать, что si = 0 на нижней лицевой поверхности). Наложим на функции Pjka условия P3ka(-1/ 2) = 0. Они обеспечивают равенство

нулю вектора напряжений на свободной границе. Тогда граничные условия при x3 = h /2 можно записать следующим образом:

Z h0-1P3ka

ю> 0

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

В классических теориях пластин принято считать прогиб w вертикальным перемещением серединной поверхности (в теории Мин-

длина [11]) или же усредненным вертикальным перемещением (в теории Рейсснера [12]). Аналогичным физическим смыслом обладает поле wk при выборе условий вида

NikLa

или

1/2

J NikLa

-1/2

В первом случае функции wk имеют смысл перемещений серединной поверхности, а во втором — усредненных перемещений вертикальных отрезков тела.

Уравнение (17) (формально бесконечного порядка) может быть стандартным образом преобразовано к последовательности уравнений конечного порядка подстановкой

<х>

wk = ^ hnwkn

n=0

и приравниванием членов при одинаковых степенях h . При этом получаются уравнения вида

2 Dkawkn, a

\a=m

в которых m — наименьший порядок ненулевых Dika, а fi выражаются через si и wkl р при l< m и |р| = 2m -l. Заметим, что в данном пункте описана общая процедура получения двумерного уравнения, в том числе и для случая sj Ф 0. Далее подробно изучается случай поперечного изгиба, когда sj = 0.

Локальные функции нулевого приближения. Покажем, как при помощи вышеописанной методики получить локальные функции нулевого приближения. В силу условия Pi3kL (-1/2) = 0 и уравнений

(14) имеем Рж (4) = Pi3kL (4) = 0. Последнее равенство эквивалентно уравнению

C 3kL + C Ъг 3 Kl = 0.

Соотношения (15) принимают вид

4

NkL

0

Теперь можно выразить Р1ЛЬ через найденные функции согласно (16):

Р1ЛЬ = СикЬ - Сиг3МгрСр3кЬ-

Заметим, что при к = 3 формула (23) приобретает вид

N3Ь

что с учетом любого из условий вида (18) или (19) превращается в N31 (4) = -5^ 4- Это значит, что нулевое приближение асимптотического метода при нулевых и отсутствии растяжения-сжатия от-счетной поверхности дает поле перемещений вида

Щ

соответствующее гипотезе Кирхгофа — Лява. Еще раз отметим, что в [5, 6] это соотношение в нулевом приближении принималось априори, а последующие приближения получались с помощью асимптотического разложения-

Уравнения третьего приближения в частном случае. Продемонстрируем полную процедуру получения уравнений третьего приближения в частном случае симметричных (относительно плоскости х3 = 0) свойств пластины.

Для начала заметим, что коэффициенты Р3ка(1/2) двумерной задачи (17) могут быть выражены через функции Р 3к р меньших порядков:

1/2 1/2

Р3каХ (1/2) = | Р3ка! (4) ¿4 = - { Рька (4) ¿4- (27) -1/2 -1/2

Если 1 = 3, то такое преобразование может быть проведено еще один раз:

1/2

Рр33каЬК (1/2) = - | Р3Кка1 (4)^4 = -1/2

1/2

= -Р3Кк*ь

-1/2

Применив к первому слагаемому правой части (28) преобразование (27), получим в итоге

1/2 ( 1 ^

Рр33ка£Я

-1/2 ^ 2 '

Для простоты дальнейшего изложения примем = 0 и Щ (х1, х2) = Щ (х1, х2) и будем искать лишь Ру3а и 3а. Нетрудно

показать, что для симметричных пластин при отсутствии горизонтально направленной составляющей нагрузки эти предположения справедливы. При к = 3 соотношение (24) принимает вид Ри 3Ь = 0. Поэтому Р133щ = 0, соответственно Р^33Ьд = Р3г-3щ = 0, а вот функции Ри 3Ьд, определяемые из соотношений (16), вообще говоря, отличны от нуля. Это значит, что в силу (14) Р1331дх также отличны от нуля. При этом Р333щк по-прежнему равны нулю. Отсюда можно

сделать вывод, что первое асимптотическое приближение позволяет получить аппроксимацию компонент а и, второе — компонент а 13, третье — а33.

Итак имеем, что все Р333а(1/2) при |а| < 4 равны нулю. Однако функции Р333щр$ (1/2) отличны от нуля и могут быть выражены с помощью (29) в виде

1/2 ( 1 ^

Р333ьдвз

Если материал пластины изотропен, то вычисления показывают,

что

Р333ЩЖ

где Б — стандартная изгибная жесткость пластины, а первое (и единственное) из уравнений (21) при условии £3 ~ И3 принимает вид

В^Ои

что соответствует уравнению для прогиба в теории Кирхгофа — Ля-ва. Подчеркнем еще раз, что это уравнение получается в третьем приближении, т. е. при учете всех слагаемых ряда (11) для а у порядка И3 и ниже. При этом заметим, что для вычисления а ¿у и получения уравнения (32) не обязательно искать функции а при |а| > 2, а достаточно найти все функции только первого приближения, так как функции Рука разных приближений связаны между собой по формулам (27), (30).

Решив уравнение (32) с граничными условиями на боковой поверхности [1, 2, 8], можно по формулам (11), (12) найти напряжения

Оу. При этом в разложении (11) имеет смысл рассматривать лишь

первые члены с ненулевыми локальными функциями (т. е. аш могут

быть найдены с погрешностью порядка И2, а 13 — с погрешностью

порядка И3, а а33 — с погрешностью порядка И4).

Рассмотрение же всех уравнений (21), соответствующих пятому приближению (т. е. при учете всех слагаемых ряда для а^ порядка

И5 и ниже), дает при том же предположении относительно £3 (т. е. £3 порядка И3) соотношение вида

характерное для теорий толстых пластин (Рейснера — Миндлина). Поле w, найденное по уравнению (33), имеет погрешность порядка

И3, поэтому формулы

позволяют найти а ш с погрешностью порядка И4, а 13 — с погрешностью порядка И5, а33 — с погрешностью порядка И6.

Сравнение с теорией толстых пластин. В качестве типичной теории толстых пластин рассмотрим теорию Миндлина [10], [11], иногда называемую теорией Тимошенко. В теории Миндлина поле перемещений (в отсутствие растяжения-сжатия) представлено в виде

и

При этом считается, что w (Х1, Х2) есть прогиб серединной поверхности. Поле деформаций представлено в виде

в1Ш

6/3

Множитель к в формуле (36) обычно принимается отличным от 1 (в то время как дифференцирование соотношений (34) дает к = 1). Это связано с тем, что распределение касательных деформаций по х3 в трехмерной задаче близко к параболическому, а в (36) в 13 не зависят от х3. Потому к обычно выбирают так, чтобы внутренняя энергия деформации соответствовала тем или иным соображениям адекватности теории [10], [17]. Наиболее распространенные значения: к = 5/6 и к - 4/5.

Теория Миндлина при р (хь х2 ) = 53 (хь х2И3 приводит к

уравнениям для прогиба щ и углов ф отклонения нормали от вертикального направления:

Уравнение (33), аналогичное (37), получается при применении асимптотического метода в пятом приближении. При наложении на локальные функции условия (18) функции щ из разложений (10) и (34) приобретают одинаковый смысл щ (хь х2, 0). Поэтому множитель к может быть найден из уравнения пятого приближения (33), в

, 5 29 V котором к = — + — ———. к 6 6 12 - 7v

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

-с поправками порядка И3, пропорциональными кубу вертикальной координаты, а в 13 с точностью до поправок распределены параболически по толщине пластины. В теории Миндлина ви линейны по толщине, но с другим угловым коэффициентом (отличающимся от -Щу на величину порядка И2, которая задается уравнениями (38)), а в 13 постоянны по толщине.

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

Достоинством предложенного метода является отсутствие кинематических гипотез (типа гипотезы Кирхгофа — Лява или Рейссне-ра), обычно используемых для анализа тонких и толстых пластин. Наоборот, кинематические гипотезы можно считать следствием результатов асимптотического анализа.

Работа выполнена частично при поддержке гранта РФФИ

13-01-00688.

ЛИТЕРАТУРА

[1] Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. Математические задачи механики композиционных материалов. Москва, Наука, 1984.

[2] Levinski T., Telega J.J. Plates, Laminates and Shells. Asymptotic Analysis and Homogenization. New York, World Scientific Publishing Co., 2010.

[3] Победря Б.Е. Механика композиционных материалов. Москва, Изд-во Моск. ун-та, 1984.

[4] Kohn R.V., Vogelius M. A new model of thin plates with rapidly varying thickness. Int. J. Solids and Struct, 1984, pp. 333-350.

[5] Шешенин С.В. Асимптотический анализ периодических в плане пластин. Изв. РАН. Механика твердого тела, 2006, № 6, с. 71-79.

[6] Шешенин С.В. Применение метода осреднения к пластинам, периодическим в плане. Вестник Моск. ун-та. Сер. 1. Математика и механика, 2006, № 1, с. 47-51.

[7] Скопцов К., Шешенин С. Асимптотический анализ слоистых пластин и пологих оболочек. Изв. РАН. Механика твердого тела, 2011, № 1, с. 161171.

[8] Димитриенко Ю.И., Яковлев Д.О. Сравнительный анализ решений асимптотической теории многослойных тонких пластин и трехмерной теории упругости. Инженерный журнал: наука и инновации, 2013, № 7(19). URL: http://engjournal.ru/articles/899/899.pdf

[9] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Асимптотическая теория конструктивно-ортотропных пластин с двухпериодической структурой. Математическое моделирование и численные методы, 2014, № 1, с. 36-57.

[10] Reddy J.N. Theory and analysis of elastic plates. Philadelphia, Taylor and Francis, 1999.

[11] Mindlin R.D. Influence of rotary inertia and shear on flexural motions of isotropic, elastic plates. ASME Journal of Applied Mechanics, 1951, vol. 18, pp. 31-38.

[12] Reissner E. The effect of transverse shear deformation on the bending of elastic plates. ASME Journal of Applied Mechanics, 1945, vol. 12, pp. 68-77.

[13] Скопцов К., Шешенин С.В. Асимптотический метод получения уравнений теории пластин Рейсснера — Миндлина. Вестник Моск. ун-та. Сер. 1. Математика и механика, 2013, № 2, с. 65-67.

[14] Ghosh S., Lee K., Moorthy S. Two scale analysis of heterogeneous elastic-plastic materials with asymptotic homogenization and Voronoi cell finite element model. Comput. Methods Appl. Mech. Enrgr, 1996, no. 132, pp. 63-116.

[15] Шешенин С.В., Савенкова М.И. Осреднение нелинейных задач в механике композитов. Вестник Моск. ун-та. Сер. 1. Математика, механика. 2012, № 5, с. 58-61.

[16] Шешенин С.В., Савенкова М.И. Об осреднении композитов при наличии нелинейности. Упругость и неупругость. Доп. материалы Международного научного симпозиума по проблемам механики деформируемых тел, посвященного 100-летию со дня рождения А.А. Ильюшина. Москва, Изд-во Моск. ун-та, 2012, с. 260-269.

[17] Yu Liu, Chee-Kiong Soh. Shear correction for Mindlin type plate and shell elements. Int. J. Numer. Meth. Engng, 2007, no. 69, pp. 2789-2806.

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

Ссылку эту на статью просим оформлять следующим образом: Шешенин С.В., Скопцов К.А. Теория пластин, основанная на методе асимптотических разложений. Математическое моделирование и численные методы, 2014, № 2, с. 49-61.

Шешенин Сергей Владимирович — д-р физ.-мат. наук, профессор, зам. заведующего кафедрой «Механика композитов» механико-математического факультета МГУ им. М.В. Ломоносова. e-mail: sheshenin@mech.math.msu.su

Скопцов Кирилл Александрович — аспирант кафедры «Механика композитов» механико-математического факультета МГУ им. М.В. Ломоносова. e-mail: arbrk1@gmail.com

Theory of plates based on the method of asymptotic decompositions

© S.V. Sheshenin, K.A. Skoptsov

Lomonosov Moscow State University, Moscow, 119991, Russia

The paper presents comparison of asymptotic analysis of a composite plate lateral bending under the impact of surface load against classical theories of thin and thick plates. Layers of the plate are assumed to be homogeneous linear elastic orthotropic materials.

Keywords: asymptotic method, averaging method, laminated plate, theory of thin plates, theory of thick plates, Kirchhoffplate, Mindlin — Reissner plate.

REFERENCES

[1] Bakhvalov N.S., Panasenko G.P. Osrednenie protsessov v periodicheskikh

sredakh. Matematicheskie zadachi mekhaniki kompozitsionnykh materialov [The averaging of processes in periodic media. Mathematical problems of mechanics of composite materials]. Moscow, Nauka Publ., 1984.

[2] Levinski T., Telega J.J. Plates, Laminates and Shells. Asymptotic Analysis and

Homogenization. N.J., World Scientific Publishing Co., 2nd ed., 2010.

[3] Pobedrya B.E. Mekhanika kompozitsionnykh materialov [Mechanics of Com-

posite Materials]. Moscow, Moscow State University Publ., 1984.

[4] Kohn R.V., Vogelius M. A new model of thin plates with rapidly varying

thickness. Int. J. Solids and Struct, 1984, no. 20, pp. 333-350.

[5] Sheshenin S.V. Izvestiya Rossiyskoy ahdemii iauk. Mkhanika Verdoгo

tela — Mechanics of Solids. Journal of the Russian Academy of Sciences, 2006, no. 6, pp. 71-79.

[6] Sheshenin S.V. Vestnik Moskovskogo Universiteta. Matem. Mekhan. — MSU

Bulletin. Mathematics.Mechanics, 2006, no. 1, pp. 47-51.

[7] Skoptsov K., Sheshenin S. Izvestiya Rossiyskoy akademii nauk. Mekhanika

tverdoгo tela — Mechanics of Solids. Journal of the Russian Academy of Sciences, 2011, no. 1, pp. 161-171.

[8] Dimitrienko Yu.I., Yakovlev D.O. Enzhinerny zhurnal: nauka i innovatsii —

Engineering Journal: Science and Innovation, 2013, 7(19). Available at: http://engjournal.ru/articles/899/899.pdf

[9] Dimitrienko Yu.I., Gubareva E.A., Sborshchikov S.V. Matematicheskoye

modelirovaniye i chislennye metody — Mathematical Modelling and Numerical Methods, 2014, no. 1, pp. 36-57.

[10] Reddy J.N. Theory and analysis of elastic plates. Taylor and Francis, Philadelphia, 1999.

[11] Mindlin R.D. Influence of rotary inertia and shear on flexural motions of isotropic, elastic plates. ASME Journal of Applied Mechanics, vol. 18, 1951, pp. 31-38.

[12] Reissner E. The effect of transverse shear deformation on the bending of elastic plates. ASME Journal of Applied Mechanics, 1945, vol. 12, pp. 68-77.

[13] Skoptsov K., Sheshenin S.V. Vestnik Moscovskogo Universiteta. Seriya 1. Matematika i Mechanika — MSU Bulletin. Series 1. Mathematics and Mechanics, 2013, no. 2, pp. 65-67.

[14] Ghosh S., Lee K., Moorthy S. Two scale analysis of heterogeneous elastic-plastic materials with asymptotic homogenization and Voronoi cell finite element model. Comput. Methods Appl. Mech. Eng., 1996, no. 132, pp. 63116.

[15] Sheshenin S.V., Savenkova M.I. VestnikMoskovskogo Universiteta. Seriya 1. Matematika i Mechanika — MSU Bulletin. Series 1. Mathematics and Mechanics, 2012, no. 5, pp. 58-61.

[16] Sheshenin S.V., Savenkova M.I. Ob osrednenii kompozitov pri nalichii nelineynosti.Uprugost' i neuprugost': dopolnitel'nye materialy Mezhdunarod-nogo nauchnogo simpoziuma po problemam mekhaniki deformiruemykh tel, posvyashchennogo 100-letiyu so dnya rozhdeniya A.A. Il'yushina. [On the Averaging of Composites with Nonlinearity. Elasticity and Inelasticity: the Additional Materials of the International Scientific Symposium on Problems of Mechanics of Deformable Bodies, dedicated to the 100th anniversary of the birthday of A. A. Ilyushin]. Moscow, Moscow State University Publ., 2012, pp. 260-269.

[17] Yu Liu, Chee-Kiong Soh. Shear correction for Mindlin type plate and shell elements, Int. J. Numer. Meth. Eng., 2007, 69, pp. 2789-2806.

Sheshenin S.V., Dr. Sci. (Phys.&Math.), professor, deputy head of the Department of Composite Mechanics, Faculty of Mechanics and Mathematics, at Lomonosov Moscow State University. e-mail: sheshenin@mech.math.msu.su

Skoptsov K.A., graduate student of the Department of Composite Mechanics, Faculty of Mechanics and Mathematics, at Lomonosov Moscow State University. e-mail: arbrk1@gmail.com