О критических нагрузках сжатой упругой прямоугольной пластины с дислокациями и дисклинациями Текст научной статьи по специальности «Математика»

МЕХАНИКА MECHANICS

УДК 539.3 DOI 10.12737/18157

0 критических нагрузках сжатой упругой прямоугольной пластины с дислокациями

*

и дисклинациями * И. М. Пешхоев1 "

1 Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация

On critical loads of compressed elastic rectangular plate with dislocations and disclinations *** I. M. Peshkhoev1**

1 Don State Technical University, Rostov-on-Don, Russian Federation

Рассматривается задача о критических нагрузках сжатой прямоугольной пластины, содержащей непрерывно распределенные источники собственных напряжений. Анализ задачи проводится на основе модификации системы нелинейных уравнений Кармана для больших прогибов упругих пластин с дислокациями и дисклинациями с различными вариантами краевых условий. Введением замены для функции напряжений задача сводится к исследованию двух задач: линейной краевой задачи относительно функции напряжений, вызванных внутренними источниками и системы нелинейных уравнений относительно прогиба и функции напряжений, вызванных внешними сжимающими нагрузками, которая имеет тривиальное решение. Классическая критическая нагрузка определяется как наименьшее собственное значение линейной краевой задачи, полученной линеаризацией системы нелинейных уравнений относительно тривиального решения. Рассматриваются четыре типа краевых условий: все края подвижно защемлены; все края шарнирно оперты; два противоположных края свободны от напряжений, а два других подвижно защемлены или шар-нирно оперты. Равномерно распределенные сжимающие нагрузки одинаковы на противоположных краях. Установлено, что если мера несовместности является нечетной по одной переменной и четной или нечетной по другой переменной, то напряжения, вызванные только внутренними источниками, не приводят к потере устойчивости плоского равновесного состояния и не влияют на критические значения сжимающих нагрузок.

A problem on critical loads of the compressed rectangular plate containing continuously distributed sources of inherent stress is considered. The task analysis is based on the modification of the Karman nonlinear equations system for large deflections of elastic plates with dislocations and disclinations under different boundary conditions. By the introduction of a replacement for the stress function, the problem reduces to the treatment of two tasks: a linear boundary value problem concerning the stress function caused by internal sources and a system of nonlinear equations concerning the deflection and the stress function caused by external compressive loads, which possesses a trivial solution. The classical critical load is defined as the smallest eigenvalue of the linear boundary value problem obtained by the linearization of the nonlinear equations system relative to the trivial solution. Four types of boundary conditions are treated: all edges are variably restrained; all edges are simply supported; two opposite edges are stress-free, and the other two are either variably restrained or simply supported. Uniformly distributed compressive loads are equal on the opposite edges. It is established that if the measure of inconsistency is odd on one variable and odd or even on another variable, then the stresses caused only by internal sources, do not lead to the loss of the flat equilibrium state and do not affect the critical values of compressive loads.

Keywords: elastic plate, dislocations and disclinations, critical load.

Введение. В работе Л. М. Зубова [1] на основе уравнений Кармана построена система нелинейных уравнений равновесия упругих пластин, содержащих внутренние источники напряжений, вызванных дислокациями и дисклина- ^ циями. Рассмотрена задача об изгибе мембраны вследствие релаксации вызванных дефектами внутренних напряже- щ

ей

ний, которая сведена к уравнению Монжа-Ампера. Найдено несколько точных решений о форме поверхности мем- g браны, содержащей распределенные дисклинации. В [2] рассмотрена задача о влиянии внутренних напряжений на ^

Ключевые слова: упругая пластина, дислокации и дисклина-ции, критическая нагрузка.

* Работа выполнена при финансовой поддержке РФФИ (проект 14-08-00142_А).

** Е-mail: peshkhoev@rambler.ru

*** The research is done with the financial support from RFFI (project 14-08-00142_А).

й о тз

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

В работе В. А. Треногина [3] был развит метод Ляпунова-Шмидта в операторной форме для нелинейных уравнений в банаховых пространствах. В работе [4] Л. С. Срубщик и В. А. Треногин исследовали задачу о влиянии малой поперечной нагрузки на выпучивание и послекритическое поведение пластины произвольной формы под действием параллельных осям координат сжимающих краевых усилий.

В работе Рейсснера [5] рассмотрена задача о влиянии нелинейно-упругого основания на выпучивание и начальное послекритическое поведение безмоментного плоско-напряженного состояния в случае бесконечной пластины с малыми геометрическими несовершенствами, а в случае тонкой сжатой пластины строго выпуклой формы со свободным краем при дополнительном действии малой поперечной нагрузки эта задача решена Л. С. Срубщиком [6]. В работах [7, 8] исследована задача о влиянии малой поперечной нагрузки на послекритическое поведение прямоугольной гибкой пластины, лежащей на нелинейно-упругом основании и равномерно сжатой в продольном направлении, а в [9] эта же задача рассмотрена для пластины с дислокациями и дисклинациями. С помощью операторного метода Ляпунова-Шмидта определено количество решений, соответствующих новым формам равновесий, и для каждого из них строятся асимптотические представления.

1. Постановка задач

1. 1. Система уравнений равновесия и краевые условия. Пусть тонкая упругая прямоугольная пластина, содержащая в плоском состоянии поля непрерывно распределенных краевых дислокаций и клиновых дисклинаций, находится под действием малой нормальной нагрузки интенсивности 40(х,у) и внешнего краевого усилия, составляющие которого вдоль осей X и у равны соответственно Р и Q. Тогда система уравнений равновесия записывается в виде [1, 10]

|бД2Ж = [Т, Ж ]+4 о

|д2Т + ЕН[Ж, Ж^2 = ЕНу. (1.1)

Б = ЕЙ3/12Ц - V2 ), [Т, Ж] = Жхх%т - 2ЖхуТу + Жуу^х

Здесь

ц(х ,У) — скалярная мера несовместности, которая выражается через плотности дислокаций и дисклинаций, Ж (х ,У ) — прогиб пластины, — функция напряжений Эри, Н — толщина пластины, V — коэффициент Пуас-

сона, Е — модуль Юнга, х, У — прямоугольные координаты с началом в центре пластины и осями, параллельными ее краям, 4 — малый числовой параметр.

Систему (1.1) будем исследовать с краевыми условиями одного из следующих типов [9, 10]

|Ж = Жх = Тху = 0, Туу = -Р, |х| = а/2, |У| < Ъ/2

|ж = Жу = Тху = 0,Тхх = |у| = Ъ/2,|х\ <а/2 .

|Ж = Жхх = Тху = 0,Туу = -Р,|х| = а/2,|у| <Ъ/2 (1 з)

[ж = Жуу = Тху = 0, Тхх = -Q, \у\ = Ъ2, |х\ < а/2 .

|Ж = Жх = Тху = 0,Туу = -Р,|х = а/2,\у\ < Ъ/2 (1

|жуу + vЖхх = Жууу + (2-v)Жхху =Тху =Тхх = 0,у = Ъ/2,|х <а/2 '

^ |Ж = Жхх =Тху = 0,Туу =-Р,|х = а2,И <Ъ/2 (15)

' й [Жуу + vЖхх = Жууу + (2 -^Жхху = Тху = Тхх = 0, |у| = Ъ/2,|^ < а/2 .

сл

Краевые условия (1.2) отвечают свободному защемлению, а (1.3) — подвижной шарнирной опоре краев. ¡¿, Сжимающие усилия равномерно распределены по краям х = + а/ 2 и у = + Ъ/2 соответственно. Условия для |у| = Ъ/ 2 в (1.4) и (1.5) соответствуют свободному краю. В этих случаях сжимающие усилия приложены только к краям

х = + а/ 2 .

Пусть Т = Ф- РУ2/2 - QX2 /2. Тогда функция Ф удовлетворяет однородным краевым условиям ФХУ = ФУУ = 0 при |Х| = а/2, ФХУ = Фхх = 0 при |У| = Ь/2 . Заменяя эти условия на эквивалентные Ф = Фх = 0 при Х| = а/2, Ф = ФУ = 0 при |У| = Ь/2 , перепишем систему уравнений (1.1) и краевые условия (1.2) — (1.5) в виде

[вДЖ + РЖхх + QWУУ = [Ф, Ж ]+4 О

(1.6)

(1.7)

(1.8) (1.9)

(1.10)

|д2Ф + ЕИ[Ж, Ж ] 2 = Ehд |Ж = Жх =Ф = Фх = 0,|х| = а/2,У <Ь/2 [ж = ЖУ =Ф = ФУ = 0,|У| = Ь/2,|х| <а/2

Ж = = Ф = Фх = 0, |х| = а2, У < Ь2

= ЖУУ = Ф = ФУ = 0, |У| = Ь/2, |х| < а/2 |Ж = Жх =Ф = Фх = 0,|х| = а/2,У <Ь/2

[ЖУУ + V Жхх = ЖУУУ + (2 - = Ф = ФУ = 0, У = Ь/ 2, | х| < а/ 2

Ж = Жхх = Ф = Фх = 0,|х| = а/2,|У| <Ь/2 [ЖУУ + vWXX = ЖУУУ + (2-v)WXXУ =Ф = ФУ = 0,У = Ь/2,|х| <а/2

В случае краевых условий (1.9) или (1.10) считаем, что Q = 0 в (1.6).

1.2. Переход к безразмерным переменным. В (1.6) — (1.10) перейдем к безразмерным переменным по формулам

х = ах, У = Ьу,Ж (х ,У) = w(x, ур!у, ф(х ,У )= ВЁ (х, у)

Р = Вр/Ь2, Q = Вд/(аа 2 ),О(х ,У ) = hdg (х, у)/(уЬ 4 ) (1.11)

а = Ь2 / а2, ц = цу2/ h2 С учетом (1.11) система (1.6) преобразуется к виду

+ Ртхх + ^уу = ак Ё ]+ 4 g

[д2аЁ + а[w, м>]2 = Ц, IX < 12, |у| < 12, (1.12)

Да О^Оуу + а( )хх , Д2а = ДаДа

Введем числовой параметр г такой, что ц(х, у) = гц1(х, у), где ц (х, у) фиксированная функция. Обозначим через ё решение краевой задачи

Дё = Ц1(х, у), (1.13)

Ё = Ёх = 0,|х = Щу < 12

Ё = Ёу = 0,| у| = 1/2, | х < 12

Тогда решением задачи ДаЁ = ц(х, у) с краевыми условиями (1.14) будет функция ёц = гё1 . Назовем переменную г параметром интенсивности напряжений, вызванных внутренними источниками пластины. Положим Ё = гё1 + Ё . Учитывая (1.11) — (114), выводим краевые задачи

|Д2аw + Р^хх + ^уу = а[w, Ё] + га[w, Ё1 ] + 4 g

[д2аЁ + а^, м]2 = 0, |х| < 1/2, |у| < 1/2, |w = Wx = 0, | х = 1/2, |у| < 12

^ = Wy = 0,|у| = 1/2, | х < 12 ^ = Wxx = 0, | х = 1/2,| у| < 12 ^ = Wyy = 0,| у| = 12, | х < 12

^ = wx = 0, | х = 1/2, |у| < 12

(1.14)

(1.15)

(1.16)

(1.17)

й И К X а X к

w

уу

+ vаwxx = Wyy + (2 - у^хху = а |у| = !2, |х < I2

(1.18)

и=ихх = 0,| x=12,| у <12 (119)

^ + vawxx = + (2 - v)awxxy = 0, |у = V2, |х| < V2 .

Функция Е в (1.15) должна удовлетворять краевым условиям (1.14).

При 4 = 0 краевая задача (1.15) с каждым из краевых условий (1.16) — (119) при всех значениях параметров р, д, г имеет тривиальное решение (и*,Е*)=(0,0).

2. Критические нагрузки выпучивания основного равновесия

2.1. Исследование линеаризованных задач. В работе [9] с использованием результатов [11-13] показано, что критические нагрузки выпучивания основного равновесия определяются из линейной краевой задачи на собственные значения

Аа® + д®уу - гак Е1] = -Р«® XX (21)

с параметрами q, г и одним из следующих краевых условий

г®=®я = 0,| х=12,| у < 12

[® = ®у = 0,|у = 1/2, | X < 12 .

г®=®хх=0,х=У2,н < 12 (23)

[® = Шуу = 0,|у = 1/2, |X < 12 .

г®=®х=0,х=12,| у < V2 (24)

[®уу + vaш хх = ®т + (2 - v)aш„у = 0, |у| = 1/2, | X < 12

|®=®хх=0,1Х=У2,Н <12 (25)

|®уу + va® хх = ®уу +(2 - v)a®хху = 0, |у| = V2, |х < V2 .

Здесь в случае краевых условий (2.4) или (2.5) в уравнении (2.1) полагается д = 0 .

Введем линейные операторы А и В

А® = + д®^ - га[®,Е1], В® = -а®хх (2.6)

Задачи на собственные значения (2.1) с краевыми условиями (2.2-2.5) в случае г = 0 исследованы, где показано, что оператор А формально самосопряжен и ограничен снизу [14]. Формальная самосопряженность оператора А в случае г ф 0 следует из равенства

Ц [<в, Е lpdxdy = Ц [<в, o:>]Fdxdy

□ □

где функция ® удовлетворяют одному из краевых условий (2.2)-(2.5), а функция Е—условиям (1.14). Ограниченность снизу оператора А в случае г ф 0 следует из неравенства

Ц[ю, Е }йdxdy

□

Л [<а, ю]Еdxdy

□

< с

Л гаА^ dxdy

□

которое следует из неравенств Фридрихса [14].

Так как в задаче на собственные значения (2.1) оператор А формально самосопряжен и ограничен снизу, а оператор В положительно определен [13], то можно применять результаты §§39, 40, 44 из [15] и свести задачу к минимизации функционала

2 J(т)= ® , (и,у) = Ц uvdxdy, (2.7)

3 (В®, га) □

и

§ где функция ® удовлетворяет одному из краевых условий (2.2)-(2.5). Приходим к следующему алгоритму [14, 16]. Пусть известны первые собственные значения р1 < р2 < ...рп-1 и соответствующие собственные функции

щ1,т^...,!®^ задачи (2.4) при заданных значениях параметров д и г. Задаем начальное приближение ®^ к соб-

<и

ственной функции ®п . Для I = 0,1... выполняем следующие действия

.&1 1) вычисляем (г +1) - е приближение к собственному значению р

£

п

рП+1) = J (®Пг)) (2.8)

2) находим гаП+1) из уравнения

+1)= р<£+1Ып) (2.9)

при соответствующих краевых условиях из (2.5)-(2.8);

3) ортогонализируем решение 3 к первым п -1 собственным функциям по энергии оператора В

юП+1) = 3П+1) - "1з; (в3П+1), 3; )/(В3;, 3;) (2.10)

] =1

4) вводим функцию 3(П+1)

(х, у ) (2.11)

3п ' = 3п 7 ЦгаП , 3 = тах 3х:

II (х,У)еО

Этот процесс повторяется до тех пор, пока не будет выполнено условие сходимости |рП+1) - рПг)| ^ £ , где е - заданная погрешность.

Критическое значение параметра г0 определяется как наименьшее собственное значение задачи

А^га + ра3 хх + д3 уу = га[га, / ] (2.12)

при заданных значениях параметров р и q с каждым из краевых условий (2.2-2.5) (в случаях (2.4), (2.5) полагается д = 0). Если в (2.12) р = д = 0, то при г = г0 пластина теряет устойчивость плоского равновесного состояния под действием только внутренних источников напряжений.

Для решения задачи (2.12) применяется минимизация функционала (2.7), где операторы А и В определяются по формулам

А3 = А^ш + р агахх + д 3уу, В3 = а[ш, / ] (2.13)

Здесь оператор А ограничен снизу [15], а знакоопределенность оператора В зависит от свойств функции /1, т.е. от свойств меры несовместности ^ (х, у).

В силу симметрии краевых условий (2.2)-(2.5) относительно начала координат при условии четности или нечетности функции /1 по каждой переменной собственные функции 3 задач (2.1) и (2.12) являются соответственно четными или нечетными по каждой переменной. Отсюда следует, что если функция является нечетной по одной переменной и четной или нечетной по другой переменной, то для оператора В в (2.13) справедливо

(В3,3) = а||[3, // ]<odxdy = а||[3, зFldxdy = 0 (2.14)

а а

и для оператора А в (2.6) имеем

= + дзуу -га[<3,/],3)=^^ + д3уу,3) (2.15)

Из (2.14), (2.15) следует, что задача (2.12) имеет только тривиальное решение 3 = 0, а в задаче (2.1) собственные значения не зависят от параметра г. Это означает, что если мера несовместности ^ (х, у) является нечетной по одной переменной и четной или нечетной по другой переменной, то напряжения, вызванные внутренними источниками, не приводят к потере устойчивости плоского равновесного состояния и не влияют на критические значения сжимающих нагрузок р и д.

Если функция // такова, что оператор В в (2.13) положительно определен, то задача (2.12) разрешима методом минимизации функционала (2.7), т.е. критическое значение г = г0 существует.

В данной работе для каждой конкретной функции /, являющейся решением задачи (1.13), (1.14) при заданной мере несовместности ^(х,у) существование решения задачи (2.12) определялось численными экспериментами, т.е. проверялась сходимость алгоритма (2.8)-(2.11), где вместо параметра рп используется гп, а операторы А и В определены формулами (2.13). Полученные при заданных значениях параметров р = р0, д = д0 собственное значение

Г0 и собственная функция го 0 задачи (2.12) сверялись с собственным значением р* и собственной функцией 3 * задаст

и

чи (2.1) при заданных значениях г = г0, д = д0. Во всех случаях сходимости алгоритма (2.8)-(2.11) наблюдалось к

К

р* « р0 и 3* « 30 с погрешностью менее 0,001. ^

и

2.2. Результаты численных расчетов. Для краевых условий (2.3) решения уравнения (2.1) при г = 0 выписы- ^ ваются явно [10]

рт,п

п4 (ат2 + п212- дп 2п2 ^(ап 2т2 ) (2.16)

47

&mn = sin (nm(x + 0,5))sin (nn(y + 0,5)) (2.17)

где m, n — натуральные числа. Для остальных краевых условий при r = 0 и для всех краевых условий при r Ф 0 задача (2.1) решалась по алгоритму (2.8)—(2.11) с применением разностного метода и квадратурной формулы Симпсона.

Формулы (2.16), (2.17) использовались в тестовых расчетах при отладке программ. Из (2.16) следует, что при

r = 0 и q = п 21 - 4а2 )/ а ,0 < а < 1/2 минимальному собственному значению pm\n = Р11 = Р21 = п2(5а + 2)/а отвечают две собственные функции г — четная по обеим переменным, и ш2 1 — нечетная по переменной x, четная по переменной y.

Рассмотрим случай постоянной меры несовместности, т.е. в задаче (1.13) цДх, y) = 1. Тогда решение задачи (1.13), (1.14) — функция F является четной по обеим переменным, что подтверждается численным решением этой задачи разностным методом. На рис.1 для случая квадратной пластины приведены графики функции Fl и её вторых производных F1,xx , F1,xy , F1,yy , которые позволяют судить о распределении напряжений, обусловленных внутренними источниками.

Рис.1. Графики функции Е1 и её вторых производных Е1хх, Е1уу, Е1ху

С помощью численных экспериментов найдены критические значения параметров нагрузок р, д, г которым отвечают одна или две собственные функции при различных значениях геометрического параметра а = Ъ2/а2 и краевых условиях (2.2)-(2.5). Отметим основные результаты для прямоугольных пластин при а = 1/4, 1, 4. Для случая г = г0, р = 0, д = 0 (пластина теряет устойчивость плоского равновесного состояния под действием только внутренних источников напряжений) установлено, что для закрепленной по краям платины (условия (2.5), (2.6)) г0 > 0 , а для пластины с двумя свободными краями (условия (2.4), (2.5)) г0 < 0. На рис.2 приведены графики первой и второй собственных функций задачи (2.12) с условиями (2.4) для квадратной пластины при р = 0, д = 0 .

й о тз

"¡3

и <U

Ü £ -Й

Рис.2. Графики первой и второй собственных функций задачи (2.12) с условиями (2.4)

для квадратной пластины

В табл.1 приведены значения первого и второго собственных значений г1, г2 задачи (2.12) с условиями (2.4) для квадратной пластины при различных значениях параметра р (д = 0).

Таблица 1

Собственные значения г1, г2 задачи (2.12) с условиями (2.4)

p Г\ Г2

0 -3794 3719

7,73 -3282 3334

15,45 -2741 2910

23,17 -2123 2424

30,89 -1371 1812

38,62 0,11 590

В табл. 2 приведены критические значения г0 для р = д = 0 при различных значения параметра а = Ь2/а2 , а = 1 и краевых условиях (2.2) - (2.5).

Таблица 2

Критические значения г0 для р = д = 0

Краевые условия а = 14 а = 1 а = 4

(2.2) 40807 4354 2550

(2.3) 27856 2595 1740

(2.4) -16997 -3794 -3088

(2.5) -10166 -2120 -1708

Заключение.

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

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

Автор благодарит профессоров Л. М. Зубова и Л. С. Срубщика за помощь в работе.

Библиографический список

1. Зубов, Л. М. Уравнения Кармана для упругой пластинки с дислокациями и дисклинациями // Доклады РАН. — 2007. — Т.412, № 3. — С. 343-346.

2. Зубов, Л. М. Сильный изгиб круглой пластинки с непрерывно распределенными дисклинациями / Л. М. Зубов, Т. Х. Фам // Известия вузов. Северо-Кавказский регион. Естественные науки. — 2010. № 4. — С. 28-33.

3. Треногин, В. А. Разветвление решений нелинейных уравнений в банаховом пространстве / В. А. Треногин // Успехи матем. наук. — 1958. — Т. 13. Вып. 4.

4. Срубщик, Л. С. О выпучивании гибких пластин / Л. С. Срубщик, В. А. Треногин // ПММ. — 1968. — Т. 32. Вып.4. — С. 721-727.

5. Reissner E. On Postbuckling Behavior and imperfection sensitivity of Thin Elastic Plates on a Non-linear Elastic

Foundation / E. Reissner // Studies in Appl. Math. — 1970. — Vol. XLIX, N. 1. — P. 45-57. Щ

X

6. Срубщик, Л. С. Краевой эффект и выпучивание тонких пластин на нелинейно-упругом основании /

Л. С. Срубщик // Дифференциальные уравнения. — 1985. — Т. XXI, № 10. — С. 1790-1794. Ц

7. Пешхоев, И. М. Выпучивание и послекритическое поведение сжатой прямоугольной пластины на нелинейно-упругом основании / И. М. Пешхоев, Л. С. Срубщик. — Ростов-на-Дону, 1983. — 17 с. — Деп. в ВИНИТИ 07.83, № 4037-83. 49

8. Баул А. В. Влияние начальных несовершенств на выпучивание продольно сжатых прямоугольных цилиндрических панелей и пластин / А. В. Баул, И. М. Пешхоев, Л. С. Срубщик // Известия вузов. Северо-Кавказский регион. Естественные науки. — 1986. — № 1. — С.34-37.

9. Пешхоев, И. М. Ветвление равновесий сжатой упругой прямоугольной пластины с дислокациями и дис-клинациями / И. М. Пешхоев // XI всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики, сб. докл., Казань, 20 - 24 августа 2015 г., — С. 2989-2991.

10. Тимошенко, С. П. Пластинки и оболочки / С. П. Тимошенко, С. Войновский-Кригер. — Москва : Физматгиз, — 1966. — 636 с.

11. Ворович, И. И. Математические проблемы нелинейной теории пологих оболочек/ И. И. Ворович. — Москва : Наука, 1989. — 376 с.

12. Морозов, Н. Ф. К нелинейной теории тонких пластин / Н. Ф. Морозов // Доклады АН СССР. — 1957. — Т.114, № 5. — С. 968-671.

13. Вайнберг, М. М. Теория ветвления решений нелинейных уравнений / М. М. Вайнберг, В. А. Треногин. — Москва : Наука, 1969. — 528с.

14. Пешхоев, И. М. Асимптотика и ветвление равновесий сжатых упругих прямоугольных пластин и стержней на нелинейно упругом основании : диссерт. ... к-та физ.-мат. наук / И. М. Пешхоев. — Ростов-на-Дону, 1991. — 146с.

15. Михлин, С. Г. Вариационные методы в математической физике / С. Г. Михлин. — Москва : Наука, 1970. — 512с.

16. Bauer, L. Block five diagonal matrices and the fast numerical solution of the biharmonic equation / L. Bauer, E. Reiss // Math. Comput. — 1972. — V.26, 118. — P. 311-326.

References

1. Zubov, LM. Uravneniya Karmana dlya uprugoy plastinki s dislokatsiyami i disklinatsiyami. [Karman equations for an elastic plate with dislocations and disclinations.] Doklady Akademii Nauk, 2007, vol. 412, no. 3, pp. 343-346 (in Russian).

2. Zubov, LM., Pham Thu Huong. Sil'nyy izgib krugloy plastinki s nepreryvno raspredelennymi disklinatsiyami. [Strong bending of circular plate with continuously distributed disclinations.] Izvestiya vuzov. Severo-Kavkazskiy region. Natural Sciences. 2010, no. 4, pp. 28-33 (in Russian).

3. Trenogin, VA. Razvetvlenie resheniy nelineynykh uravneniy v banakhovom prostranstve. [Branching of solutions of non-linear equations in Banach space.] Russian Mathematical Surveys, 1958, vol. 13, iss. 4 (in Russian).

4. Srubshchik, L.S., Trenogin, VA. O vypuchivanii gibkikh plastin. [On buckling of flexible plates.] Journal of Applied Mathematics and Mechanics, 1968, vol. 32, iss. 4, pp. 721-727 (in Russian).

5. Reissner E. On Postbuckling Behavior and imperfection sensitivity of Thin Elastic Plates on a Non-linear Elastic Foundation. Studies in Appl. Math., 1970, vol. XLIX, no. 1, pp. 45-57.

6. Srubshchik, L.S. Kraevoy effekt i vypuchivanie tonkikh plastin na nelineyno-uprugom osnovanii. [The edge effect and buckling of thin plates on nonlinear elastic foundation.] Differential Equations, 1985, vol. XXI, no. 10, pp.1790-1794 (in Russian).

7. Peshkhoev, I.M. Srubshchik, L.S Vypuchivanie i poslekriticheskoe povedenie szhatoy pryamougol'noy plastiny na nelineyno-uprugom osnovanii. [Buckling and post critical behavior of compressed rectangular plate on nonlinear elastic foun-

| dation.] Rostov-on-Don, 1983, 17 p. Dep. in VINITI 07.83, no. 4037-83 (in Russian).

и 8. Baul, А.^, Peshkhoev, I.M., Srubshchik, L.S Vliyanie nachal'nykh nesovershenstv na vypuchivanie prodol'no

С

szhatykh pryamougol'nykh tsilindricheskikh paneley i plastin. [Influence of initial imperfections on the buckling of longitudinally compressed rectangular cylindrical panels and plates.] Izvestiya vuzov. Severo-Kavkazskiy region. Natural Sciences.1986, no. 1, pp. 34-37 (in Russian).

9. Peshkhoev, I.M. Vetvlenie ravnovesiy szhatoy uprugoy pryamougol'noy plastiny s dislokatsiyami i disklinatsiyami. [Branching of equilibria of compressed elastic rectangular plate with dislocations and disclinations.] XI vse-rossiyskiy s''ezd po fundamental'nym problemam teoreticheskoy i prikladnoy mekhaniki, sb. dokl., Kazan', 20 - 24 avgusta 2015 g. [Proc. XI All-Russian Congress on fundamental issues of Theoretical and Applied Mechanics, Kazan, 20 -. 24 August, 2015] pp. 2989-2991 (in Russian).

"¡3

и (U

Ü £ л

10. Timoshenko, S.P., Voynovsky-Krieger, S. Plastinki i obolochki. [Plates and shells.] Moscow: Fizmatgiz, 1966, 636 p. (in Russian).

11. Vorovich, I.I. Matematicheskie problemy nelineynoy teorii pologikh obolochek. [Mathematical problems of nonlinear shallow-shell theory.] Moscow: Nauka, 1989, 376 p. (in Russian).

12. Morozov, N.F. K nelineynoy teorii tonkikh plastin. [On nonlinear theory of thin plates.] Doklady AN SSSR, 1957, vol. 114, no. 5, pp. 968-671 (in Russian).

13. Weinberg, M.M., Trenogin, V.A. Teoriya vetvleniya resheniy nelineynykh uravneniy. [Theory of branching of nonlinear equation solutions] Moscow: Nauka, 1969, 528 p. (in Russian).

14. Peshkhoev, I.M. Asimptotika i vetvlenie ravnovesiy szhatykh uprugikh pryamougol'nykh plastin i sterzhney na nelineyno uprugom osnovanii: dissert. ... k-ta fiz.-mat. nauk. [Asymptotics and branching of equilibria of compressed elastic rectangular plates and rods on nonlinear elastic foundation: Cand.Sci. (Phys.-Math.) diss.]. Rostov-on-Don, 1991, 146 p. (in Russian).

15. Mikhlin, S.G. Variatsionnye metody v matematicheskoy fizike. [Variational methods in mathematical physics] Moscow: Nauka, 1970, 512 p. (in Russian).

16. Bauer, L., Reiss, E. Block five diagonal matrices and the fast numerical solution of the biharmonic equation. Math. Comput., 1972, vol. 26, 118, pp. 311-326.

Поступила в редакцию 28.08.2015 Сдана в редакцию 28.08.2015 Запланирована в номер 22.01.2016

ев И К X св X

К