Неустойчивость упругих цилиндров при одноосном растяжении

Проблема устойчивости механических систем имеет большое научное и прикладное значение, поскольку расчет любых строительных и инженерных конструкций основывается на анализе способности этих конструкций выдерживать различные внешние нагрузки и возмущения. Соответственно, вопрос устойчивости системы под воздействием внешних факторов часто выходит на первый план. Одним из самых первых исследователей устойчивости был Архимед, его работы носили статический характер и относились к твердым телам, погруженным в несжимаемую упругую жидкость. В дальнейшем было предложено большое количество различных определений понятия «устойчивость» и, соответственно, различных теорем об устойчивости, но большинство классических исследований устойчивости упругих систем всё же основывается на том или ином статическом критерии.

Во многих случаях состояние равновесной системы можно определить двумя параметрами — характерным перемещением v и параметром нагрузки Р. Тогда всей совокупности состояний равновесия соответствует некоторая кривая в системе осей v, Р. С помощью этой кривой можно предсказать поведение системы при монотонном возрастании параметра нагрузки (мягкое нагружение) или параметра перемещения (жесткое нагружение), а также отметить критические значения параметров. Возникновение различных методов отыскания критических состояний, соответствующих потере устойчивости, объясняется разнообразием свойств механических систем и кривых v, Р.

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

Во втором варианте статического метода в решение с самого начала вводятся те или иные неидеальности (несовершенства): начальные прогибы, начальные эксцентриситеты или дополнительные внешние силы. При решении часто пользуются линеаризованными уравнениями. Некоторая условность этого подхода состоит в линеаризации задачи, хотя рассматриваемые перемещения не являются малыми. Во многих случаях представляется более правильным учитывать нелинейности, неизбежно проявляющиеся при больших перемещениях. 5.

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

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

В историческом контексте первыми получили развитие различные теории устойчивости конструкций на основе одномерных моделей стержней и двумерных моделей пластин и оболочек [10, 26, 35, 69, 72−74].

Потеря устойчивости тонких и тонкостенных тел при сжатии происходит в области малых деформаций, что позволяет использовать физически линейные определяющие соотношения (закон Гука). Неустойчивость при растяжении наступает обычно при значительных деформациях, что требует полного учета физической нелинейности и трехмерности изучаемых тел. Кроме того, исследование неустойчивости при растяжении, как правило, невозможно в рамках одномерных моделей стержней или двумерных моделей пластин и оболочек. Следует также отметить, что потеря устойчивости при растяжении может рассматриваться как один из механизмов разрушения конструкций [42].

Исследованием вопросов устойчивости на основе трехмерных уравнений нелинейной теории упругости занимались такие отечественные ученые, как B. JL Бидерман [2], В. В. Болотин [4], А. Н. Гузь [12−14],.

B.А. Еремеев [18], Л. М. Зубов [23, 31−33, 35], АЛО. Ишлинский [39], А. И. Лурье [52−54], В. В. Новожилов [60], К. Ф. Черных [79] и другие.

Среди зарубежных авторов следует упомянуть таких ученых, как D. Bigoni [84], М.А. Biot [85], А.Е. Green [92, 93], D. Joseph [40], J.C. Patterson [103], C.E. Pearson [104], S. Reese [106], R.S. Rivlin,.

C.B. Sensenig [111], R.T. Shield [92], C. Truesdell [75, 82], Z. Wesolowski [114,115].

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

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

Существенный вклад в развитие нелинейной теории упругости внесли такие отечественные ученые, как Н. В. Зволинский [19], J1.M. Зубов [24−29, 116], В. А. Левин [49], А. И. Лурье [53, 54], Н. Ф. Морозов [55], В. В. Новожилов [61, 62], Л. А. Толоконников, К. Ф. Черных [79], а также зарубежные ученые: G. Adkins [11], S.S. Antman [81], J.L. Ericksen [90], A.E. Green [11], R. Hill, W. Noll [112], R.W. Ogden [101], R.S. Rivlin [107−110], C. Truesdell [75, 82,112].

Одним из распространенных конструктивных элементов является стержень, работающий на растяжение-сжатие. Проблема устойчивости этого элемента при сжатии хорошо изучена [31, 76, 87, 91], чего нельзя сказать о случае растяжения. Существенный интерес представляет предсказание возникновения неустойчивости при растяжении, а также теоретическое описание закритического поведения растягиваемого стержня, то есть поведения после потери устойчивости, которая выражается в исчерпании стержнем несущей способности с последующим разрушением. Хрупкие материалы (чугун, бетон и др.) разрушаются без видимых предварительных деформаций, тогда как пластичные материалы сталь, медь, некоторые сплавы и др.) перед разрушением испытывают существенные формоизменения, в месте будущего разрушения стержня возникает и развивается шейка.

Процесс образования шейки достаточно хорошо изучен экспериментально [3, 16, 56, 77, 78, 80, 89, 98, 99, Л02, 113], получены формы деформированных образцовзависимость вида шейки от температуры, скорости растяженияхарактер возникновения трещин на начальных стадиях разрушения образца и другие экспериментальные данные. На основе этих данных построен ряд теорий математического описания профиля шейки. Например, в работе [100] дается описание «формы шейки с помощью двух поверхностей вращения. Предложенная модель имеет один параметр, который зависит от материала. Применение указанной модели дает хорошее соответствие фактическим профилям шеек, которые наблюдаются в образцах из мягкой и полутвердой стали, меди, свинца. Разработана также математическая теория образования шейки, согласно которой материал считается наделенным некоторыми идеализированными свойствами [37, 38], а именно считается, что он обладает определенным пределом текучести, после достижения которого деформируется, как вязкое вещество.

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

Из опытов на простое растяжение хорошо известно, что после достижения максимума на диаграмме растяжения стержня процесс однородного деформирования становится неустойчивым, образуется шейка. Здесь, однако, следует четко различать две различных диаграммы растяжения, использующиеся в зависимости от целей конкретного исследования [42, 56, 86, 88, 94]. Это диаграмма растяжения, на которой представлена зависимость условного растягивающего напряжения (то есть растягивающей силы, отнесенной к начальной площади поперечного сечения стержня) от удлинения образцаи диаграмма, на которой представлена зависимость истинного растягивающего напряжения (то есть растягивающей силы, отнесенной к актуальной площади поперечного сечения стержня) от удлинения образца. Первую диаграмму обычно называют условной диаграммой растяжения, вторую называют истинной.

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

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

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

Следует отметить, что существует более строгий метод изучения закритического поведения трехмерных упругих тел, основанный на операторном методе Ляпунова-Шмидта [5] исследования ветвления решений нелинейных задач [20−23, 43, 117]. Указанный метод позволяет точно определить количество решений статических краевых задач нелинейной теории упругости и их характер в окрестности точки бифуркации, однако его применение дает возможность исследовать лишь начальное закритическое поведение деформируемых тел, когда параметр нагружения мало отличается от критического значения. Метод, разработанный в данной работе, позволяет исследовать закритическое поведение трехмерного нелинейно-упругого тела при любых деформациях, что делает его более универсальным.

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

Экспериментальному исследованию некоторых аспектов закритического деформирования посвящены работы П. Бриджмена [3], Р. А. Васина и др. [6], С. Д. Волкова [9], А. А. Лебедева и Н. Г. Чаусова [48], Ф. С. Савицкого и Б. А. Вандышева [67], В. В. Стружанова и В. И. Миронова [71], Я. Б. Фридмана и Б. А. Дроздовского [17].

Различные подходы к теоретическому объяснению данного механического явления представлены в работах С. Д. Волкова [7, 8], В. А. Ибрагимова и В. Д. Клюшникова [36], A.M. Линькова [50, 51], Л. В. Никитина и Е. И. Рыжака [57−59, 64−66], В. В. Стружанова [70] и ряда других ученых.

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

Во многих исследованиях закритического деформирования рассматривается поведение разупрочняющихся материалов, а именно материалов, для которых на определенном этапе имеет место убывающая зависимость истинных напряжений от удлинений. В частности, вопросы, близкие к теме настоящей работы, обсуждались с другой точки зрения в статьях [59, 64−66], где показана физическая осуществимость падающего участка истинной диаграммы нагружения при сохранении свойства сильной эллиптичности материала (строгого неравенства Адамара).

В настоящей работе рассматривается материал, для которого зависимость истинного напряжения от удлинения является монотонно возрастающей, то есть участок разупрочнения отсутствует. Падающий участок наблюдается на диаграмме зависимости условного напряжения от удлинения. Исследованное в [59, 64−66] разупрочнение материала может иметь место на стадии достаточно развитой шейки. Однако в данном исследовании начало образования шейки не связывается с разупрочнением материала, а рассматривается с точки зрения явлений упругой неустойчивости, аналогичных выпучиванию конструкций.

Перейдем к изложению содержания работы. Первая глава посвящена постановке задачи об одноосном растяжении кругового цилиндра. В п. 1.1 рассмотрен общий случай осесимметричной деформации кругового цилиндра из однородного изотропного материала. Получены представления градиента деформации, меры деформации Коши-Грина, логарифмической меры деформации (тензора Генки). Определены собственные значения и собственные векторы указанных мер деформации.

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

В п. 1.2 рассмотрено однородное состояние одноосного растяжения. Представлен вид решения краевой задачи, соответствующий однородному напряженно-деформированному состоянию цилиндра. Показано, что для рассматриваемого материала условная диаграмма растяжения имеет точку максимума, тогда как истинная диаграмма монотонно возрастает.

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

В п. 2.2 вычислено возмущение производной удельной потенциальной энергии деформации по логарифмической мере деформации.

П. 2.3 и п. 2.4 содержат промежуточные преобразования, необходимые для выражения уравнений равновесия тела через три неизвестных функции: добавочные перемещения по осям г и z, а также возмущение давления в несжимаемом теле.

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

В п. 2.7 осуществлен вывод линеаризованных граничных условий на боковой поверхности и на торцах цилиндра. Тем самым получена линейная однородная краевая задача, которая всегда имеет тривиальное решение.

Согласно бифуркационному критерию устойчивости равновесия, исследование устойчивости в малом сводится к нахождению спектра критических значений параметра деформации, при которых указанная краевая задача имеет нетривиальные решения, а также к определению собственных функций — мод выпучивания. В силу «концепции продолжающегося нагружения» бифуркационный критерий пригоден с прикладной точки зрения также и для исследования устойчивости упруго-пластических тел.

Поиск нетривиальных решений краевой задачи производился с использованием специального вида решения, описанного в п. 2.8. Данный вид решения позволяет удовлетворить краевым условиям на торцах цилиндра.

В пп. 2.9, 2.10, 2.11 и 2.12 произведено преобразование соответственно первого уравнения равновесия, второго уравнения равновесия, условия несжимаемости и граничных условий на боковой поверхности цилиндра с учетом специального вида решения. Для обеспечения разрешимости краевой задачи в п. 2.12 в рассмотрение введены два краевых условия на оси цилиндра, которые вытекают из требований неразрывности материала и гладкости решения на оси цилиндра.

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

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

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

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

В третьей главе на основе энергетического метода исследования устойчивости изучено закритическое деформирование растягиваемого цилиндра. В п. 3.1. дана вариационная постановка задачи о растяжении цилиндра. Решение линейной однородной задачи определено во второй главе с точностью до произвольного постоянного множителя. Амплитуды выпучивания и число закритических форм равновесия не могут быть найдены из линеаризованной краевой задачи устойчивости, для анализа закритического поведения цилиндра необходимо рассмотреть нелинейные уравнения, описывающие деформацию упругого тела. Исследование закритического поведения цилиндра основано на применении вариационного принципа Лагранжа, согласно которому выполнение уравнений равновесия, записанных в перемещениях, и граничных условий эквивалентно стационарности потенциальной энергии тела. Основываясь на этом принципе, равновесное состояние считаем устойчивым, пока потенциальная энергия имеет в нём минимальное значение. Если же в процессе нагружения возникает возможность существования равновесного состояния с потенциальной энергией меньшей, чем в текущем состоянии, то текущая форма равновесия теряет устойчивость, сменяясь более выгодной с энергетической точки зрения новой формой равновесия. Такой подход к исследованию устойчивости позволяет, в отличие от линеаризованной задачи, не ограничиваться лишь малыми отклонениями от однородного деформированного состояния.

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

В п. 3.2 описан численный расчет потенциальной энергии и её производных. При вычислении указанных величин удлинение цилиндра и значения амплитуд мод выпучивания считаются фиксированными. Удельная потенциальная энергия и её производные заменяются сеточными функциями по г и 2, и для вычисления потенциальной энергии тела и её производных применяется формула Симпсона.

В п. 3.3 описан итерационный метод поиска решений нелинейной задачи, основанный на замене Ммерных поверхностей, соответствующих первым производным потенциальной энергии, плоскостями, касательными к данным поверхностям в некоторой точке.

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

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

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

По теме диссертации опубликованы работы [30,45−47, 95−97]. Из них работы [30, 95] выполнены и опубликованы в соавторстве с научным руководителем Зубовым JI.M. В указанных работах Зубову JI.M. принадлежат постановка задачи и выбор методов исследования, соискателю принадлежат реализация методов исследования, аналитические выкладки, выполнение и анализ численных расчетов.

Заключение

.

В диссертации получены следующие основные результаты:

1. Рассмотрен общий случай осесимметричной деформации растягиваемого кругового цилиндра из однородного изотропного материала. Получены представления градиента деформации, меры деформации Коши-Грина, логарифмической меры деформации (тензора Генки). Определены собственные значения и собственные векторы указанных мер деформации.

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

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