Численные решения двумерных технологических задач теории пластичности тема диссертации и автореферата по ВАК РФ 01.02.04, кандидат физико-математических наук Песков, Александр Владимирович

• Песков, Александр Владимирович
• кандидат физико-математических наук
• 1983

Оглавление диссертации кандидат физико-математических наук Песков, Александр Владимирович

ВВЕДЕНИЕ

Глава I. Построение математической модели решения плоской задачи теории пластичности

1.1. Постановка задачи и модель изотропно упрочняющегося вязко-пластического мате -риала

1.2. Численный метод решения задачи

1.3. Построение численной схемы

1.3.1. Построение системы уравнений

1.3.2. Алгоритм решения уравнений пластичности. 27 1.3<3. Доказательство единственности решения . 32 1.3.4. Исследование сходимости решения

Глава 2. Задача о кристаллизации стального двутав рового профиля

2.1. Кристаллизация в свободном пространстве.

2.2. Кристаллизация в замкнутом пространстве.

2.3. Кристаллизация в частично замкнутом пространстве

Глава 3. Задача о деформации стальной полосы в неоднородном температурном поле наклон -ными штампами

3.1. Деформация однородной по структуре полосы

- з

3.2. Деформация полосы, имеющей дефектные образования

3.3. Деформация полосы штампом сложной конфигурации

Глава 4. Задача о деформации стальной неоднород -ной полосы под действием системы плоскопараллельных штампов

4.1. Деформация полосы под действием несим -метрично приложенной нагрузки

4.2. Внедрение штампов различной конфигурации в ограниченное полупространство

4.3. Деформация биметаллической полосы под действием системы штампов

Глава 5. Применение полученных решений к анализу некоторых технологических процессов

Введение диссертации (часть автореферата) на тему «Численные решения двумерных технологических задач теории пластичности»

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

Большой вклад в развитие теории упругости и теории пластичности внесли А.А.Ильюшин [28], Д.Д.Ивлев [ 5], Л.С.Лейбен-зон [50], Н.И.Мусхелишвили [64], А.И.Лурье [51], А.Надаи [55], [56], В.В.Новожилов [9,52], В.Прагер [4,57], Ю.Н.Работнов [8], Л.И.Седов [53,54j, В.В.Соколовский [7], Р.Хилл [3J, Ф.Ходж [4] и другие ученые.

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

Одним из первых методов, с помощью которого были решены многие важные практические задачи, был метод характеристик или линий скольжения. Он был разработан для решения задач пластического течения при использовании модели жестко-пластического тела. Основные теоретические положения метода разра -ботаны Г.Генки и Л.Прандтлем [ 1,2]. Большой вклад в разви -тие метода внесли Р.Хилл [3], В.Прагер [4], Ф.Ходж [41, А.Фрейденталь [II], Х.Гейрингер [II] и другие ученые. Общее решение плоской задачи получено С.А.Христиановичем [ 10]. В случае плоской деформации, напряженное состояние в точке полностью определено, если известно направление линий скольжения и величина среднего напряжения. Для частных задач, учиты -вая свойства линий скольжения и граничные условия, можно оп -ределить поля линий скольжения, которые являются характерис -тиками дифференциальных уравнений плоской деформации. Сравнительная легкость, с которой существующие методы анализа могли быть применены в этом случае, послужила широкому внедрению метода линий скольжения к решению двумерных технологических задач. Одной из первых была рассмотрена задача о начальном те чении при вдавливании плоского жесткого штампа в полубеско -нечное тело. Полученные Л.Прандтлем и Р.Хиллом решения [2,3] отличаются полем скоростей, но дают одинаковое значение дав -ления штампа. В работах В.Прагера [4] и Г.И.Быковцева [13] по строены решения, являющиеся комбинацией решений Л.Прандтля и Р.Хилла. Это говорит о том, что при использовании схемы идеального жесткопластического тела возможна неоднозначность решения.

В дальнейшем метод линий скольжения получил развитие при исследовании таких технологических процессов как листовая вытяжка, выдавливание, прокатка, ковка, волочение в работах Г.И. Быковцева [6,14,15,25,26], Б.А.Друянова [16,17,19-23], Д.Д. Ивлева [5,6], Р.И.Непершина [24], В.Прагера [4], В.В.Соколовского [ 7,18], Р.Хилла [3], Ф.Ходжа [4], А.Д.Томленова [ 27] и других ученых.

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

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

В указанных вше работах получены решения для начального, установившегося и неустановившегося пластических течений, для однородного и неоднородного материала, с учетом и без учета трения на границе инструмента с деформируемым материалом. Основными и наиболее простыми методами применяемыми для построения полей линий скольжения, являются численные, основанные на переходе к конечно-разностным соотношениям и использовании свойств линий скольжения. Различные варианты таких построений приведены в работах Р.Хилла [з], В.Прагера [4] В.В.Соколовского [7]. Несмотря на широкое развитие, метод имеет ряд допущений при рассмотрении поведения реальных материалов: не учитываются действительные свойства материалов--температура, степень и скорость деформации; большие затруд нения вызывает построение линий скольжения при исследовании сложной деформации.

Изучение деформации материала с отражением его действительных свойств может быть осуществлено с применением вариационных методов. Основоположниками их применения к задачам теории пластичности и упругости являются А.А.Гвоздев [751, А.А.Ильюшин [28], А.А.Марков [29], С.Г.Михлин [30], В.В.Но -вожилов [9], В.Драгер и Ф.-Ходж [4], Р.Хилл [з] и другие ученые.

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

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

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

Развитием вариационного метода в этом направлении является метод локальных вариаций,предложенный Ф.Л.Черноусько [ 31] и развитый в работах [32-35]. Идея этого метода заключается в том, что область, ограниченная контуром, разбивается на прямоугольные ячейки.Уеловия на контуре переносятся на граничные точки полученной прямоугольной области. Интеграл по области переписывается в виде суммы интегралов по каждой ячейке. Этот интеграл приблизительно определяется через усредненные по ячейке значения функций скорости перемещения и ее производных, представленных в конечно-разностном виде. При варьировании функции скорости в узле изменяется не вся сумма интегралов, а только те интегралы, которые включают эту функцию. Поэтому варьирование носит локальный характер.

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

Дня решения дифференциальных уравнений в частных производных был разработан метод конечных разностей - работы В.Ва зова и Дж.Форсайта Г 38], Ш.Е.Микеладзе [37], Л.Коллатца [Зб1 А.А.Самарского [39] и других ученых. Наибольшее распростра -нение этот метод получил при решении плоских задач [40-42] . Внимание к этому методу в последнее время обусловлено внед -рением в практику расчетов ЭВМ и успешным использованием аппарата матричной алгебры. Это привело к упрощению записи алгоритма задач и возможности решения задач, имеющих боль -шой объем вычислений.

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

К недостаткам метода следует отнести неоднородность ко

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

В последнее время, наверное, основным численным методом решения прикладных задач стал метод конечных элементов. Возникший в начале 30-х годов, как инженерный метод расчета на прочность металлических конструкций, он получил математиче -ское обоснование в трудах Дж.Аргириса [44,65], О.Зенкевича [43,45,46], Р.Клафа [47], В.Г.Корнеева [58,67]. Согласно это му методу, область, занимаемая телом, разбивается на элементы простой формы. Для плоского случая это чаще всего тре -угольники, а для пространственного - тетраэдры. По каждому элементу задаются некоторые функции, позволяющие определить перемещения внутри элементов по перемещениям в узлах. Определяя затем деформации и напряжения, получают выражение для энергии, как функцию узловых смещений, при вариации ко -торых формируется система алгебраических уравнений. Из ее ре шения определяются искомые перемещения узлов. Необходимо отметить, что распределенная нагрузка заменяет^сяГэквивалентной узловой. Метод конечных элементов часто трактуется как метод Ритца. Различие между ними заключается в выборе системы ко -ординатных функций. В методе Ритца эти функции задаются для всей рассматриваемой области, а в методе конечных элементов-дпя каждого элемента и через множество этих функций опреде -ляется состояние всей системы. В первом случае варьируют по параметрам, имеющимся в каждом члене ряда, а во втором - по перемещениям узлов. Необходимо отметить, что метод конечных элементов хорошо разработан для решения упругих и упруго-пла стических задач - работы В.Г.Корнеева [59], Л.А.Розина [60],

Дж.Аргириса [44], О.Зенкевича [46], Р.Клафа [471 Особенности постановки и конечно-элементной аппроксимации для задач с жестко-пластической и вязко-пластической средами рассмотрена в работах В.Н»Сегала и Г.П.Свирида 48,49 , В.Г.Корнеева и С.Е.Пономарева 68 ,

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

Эти недостатки отсутствуют в численном методе решения дифференциальных уравнений пластического течения,предложен -ным В.И.Одиноковым [61]. Суть его состоит в следующем. Рас -сматриваемая область делится на элементы конечных размеров, для каждого элемента записывается полная система дифференциальных уравнений пластического течения в конечно-разностной форме через значения скоростей перемещений и напряжений по граням элемента. При записи полной системы уравнения состояния получаются из условия коаксиальности девиаторов напряжений и скоростей деформаций,на наклонных площадках, примыкающих к поверхности области течения, записываются уравнения Ко-ши. С учетом граничных условий получается замкнутая опреде -ленная система алгебраических уравнений. Данный метод позво -ляет построить такие рекуррентные соотношения, при которых значительно уменьшается количество неизвестных (в 6-7 раз). В полученной эквивалентной системе коэффициенты вычисляются численным методом. Весь процесс состоит из нескольких итераций. В каждой из них коэффициент пропорциональности между девиато-рами напряжений и скоростей деформаций, который является функ цией от температуры тела, степени и скорости деформации, принимается постоянным. Для решения системы линейных уравнений используется метод Гаусса. Результат решения - поля напряже -ний и скоростей перемещений по граням элементов. Решение этим методом ряда модельных и прикладных задач показывает соответствие полученных результатов экспериментальным - работы В.И. Одинокова, Б.Г.Каплунова, Е.И.Макеранца [63, 69-74].

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

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

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

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

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

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

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

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

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

В приложении представлены акты о внедрении результатов теоретических исследований и программа решения задач на языке Фортран.

Заключение диссертации по теме «Механика деформируемого твердого тела», Песков, Александр Владимирович

ЗАКЛЮЧЕНИЕ

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

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

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

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

Разработан общий алгоритм решения плоской задачи при наличии общих граничных условий. Составлена программа на языке "Фортран".

Алгоритм решения применим к исследованию деформаций не одаородных по структурным и температурным характеристикам сред с учетом истории нагружения.

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

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

4. Применение численной схемы рассмотрено на ряде за -- дач, решение некоторых из них стало возможным только при разработке данной численной схемы: а. Задача об определении напряженно-деформированного состояния при кристаллизации двутаврового профиля в свободном, замкнутом и частично замкнутом пространствах. Результаты ее решения позволили исследовать технологический процесс получения двутавровых заготовок на машинах непрерывного литья. Определена наилучшая форма кристаллизатора из критерия наименьших растягивающих напряжений, возникающих в затвердевающей корочке. Рекомендации по данной работе выданы НИИтяжмашу.Ожидаемый экономический эффект 100 тыс.руб. б. Задача о деформации стальной полосы наклонными штам -пами, совершающими плоско-параллельное движение позволила определить энергосиловые параметры и кинематику течения металла в области деформирования для нового разрабатываемого агрегата периодической деформации. Рассмотрена задача о деформации полосы, имеющей дефектные образования в поверхностных слоях. в. Задача о деформации полосы под действием несимметрично приложенной нагрузки дала возможность рассмотреть процесс резания полосы и определить зависимость между технологическими параметрами этого процесса. Результаты решения задачи позволили выработать рекомендации, которые внедрены на ЧМЗ. Го -довой экономический эффект составил 200 тыс.руб. г. Рассмотрена задача о внедрении плоского и фигурного штампов в ограниченную область. Ее решение использовано при исследовании процесса пробивки фигурного паза в рельсовой под кладке. Рекомендации по данному процессу использованы Урал -НИИЧМ при проектировании штампа на КМК. Ожидаемый экономиче -ский эффект 200 тыс.руб.

Список литературы диссертационного исследования кандидат физико-математических наук Песков, Александр Владимирович, 1983 год

1. Генки Г. О некоторых статически определимых случаях равновесия в пластических телах. - В кн.: Теория пластичности. М., Изд.иностр.лит., 1948, с.80-101.

2. Прандтль Л. О твердости пластических материалов и сопротив -лении резанию. В кн.: Теория пластичности* М., Изд.иностр. лит., 1948, с.70-79.

3. Хилл Р. Математическая теория пластичности. М., ГИТТЛ, 1956, 407 с,

4. Прагер В,, Ходж Ф» Теория идеально пластических тел. М.,Изд. иностр.лит., 1956, 398 с*

5. Ивлев Д.Д, Теория идеальной пластичности. М.» Наука, 1966, 232 с.

6. Ивлев Д.Д., Быковцев Г.И. Теория упрочняющего пластического тела. М*, Наука, 1971, 231 с.

7. Соколовский В.В* Теория пластичности. М., Высшая школа, 1969 608 с.

8. Работнов Ю.Н. Механика деформируемого твердого тела. М.,Наука, 1979, 744 с.

9. Новожилов В.В. Теория упругости. Л., Судпромгиз, 1958, 370 с.

10. Христианович С.А. Плоская задача математической теории пла -стичности при внешних силах, заданных на замкнутом контуре. Матем.сб., 1936, тД, № 4, с.511-543.

11. Друянов Б,А, Метод решения статически неопределимых задач плоского течения идеально-пластических тел, ДАН СССР,1964 т.143, № 4, с,808.

12. Ильюшин А.А» Пластичность. Гостехиздат, 1948, 376 с»

13. Марков А.А. О вариационных методах в теории пластичности. ПММ, АН СССР, 1947, т.Х1

14. Микеладзе Ш.Е. Численные методы интегрирования дифференци -альных уравнений в частных производных. М., Изд-во АН СССР, 1936, 108 с.

15. Вазов В., Форсайт Дж, Разностные методы решения дифферен -циапьных уравнений в частных производных. М., ИЛ, 1963, 487 с.

16. Самарский А.А. Введение в теорию разностных схем. Наука, М., 1972, 552 с.

17. Варвак П,М. Развитие и приложение метода сеток к расчету пластинок. Изд.АН УССР, ч.1, 1949, 136 е.; ч.И, 1952, 116 с.

18. Длугач М.И. Метод сеток в смешанной плоской задаче теории упругости. Киев, "Наукова думка", 1964, 260 с.

19. Дятловицкий Л.И. Напряжения в гравитационных плотинах на нескальных основаниях. Киев, Изд.АН УССР, 1959, 339 с.

20. Сегал В,М., Свирид Г.П, Исследование стационарного течения жестко-пластического материала численным методом конечных элементов. Прикладная механика, 1973, т.IX, вып.4, с.76-80

21. Сегал В.М., Свирид Г.П, Численное исследование асимметричных вязко-пластических течений методом конечных элементов. Доклад АН БССР, 1973, т.ХУП, № 9.

22. Лейбензон Л,С. Курс теории упругости. М., Гостехиздат,1947 464 с.

23. Лурье А.И. Теория упругости. Наука, 1970, 939 с,

24. Новожилов В,В. Основы нелинейной теории упругости. Гостех-издат, 1948, 211 с,

25. Седов Л.И, Механика сплошных сред. Наука, 1976, т.1, 536 с

26. Седов Л,И. Механика сплошных сред. Наука, 1976, т.2, 576 с

27. Надаи А. Пластичность и разрушение твердых тел, т.1, ИЛ., 1954, 648 с.

28. Надаи А. Пластичность и разрушение твердых тел. т.2, ИЛ,, 1954, 863 с.

29. Прагер В. Введение в механику сплошных сред. ИЛ., 1963, 311 с.

30. Корнеев В,Г, Сопоставление метода конечных элементов с вариационно-разностным методом решения задач теории упруго -сти. Изв.ВНИИГ им.Б.Е.Веденеева", 1967, т.87, с,286-307.

31. Корнеев В.Г. 0 методе конечных элементов для решения задач упругого равновесия» В сб.: "Строительная механика соору жений". Изд.ЛПИ им.М.И.Калинина, М», 1971, с,28-46.

32. Розин Л.А, Расчет гидротехнических сооружений на ЭЦВМ. Метод конечных элементов. Изд-во "Энергия", 1971, 214 с,

33. Одиноков В,И. Численный метод решения дифференциальных уравнений пластического течения. Прикладная механика, 1973 т.IX, вып.12, с.64-70»- 102

34. Одиноков В.И. Численное решение некоторых задач о деформаг ции несжимаемого материала. Прикладная механика. 1974, т.Х вып.1, с.84-91.

35. Макеранец Е.И., Одиноков В.И. Расчет пластического течения полых овальных цилиндров неограниченной длины. Изв.АН СССР, МТТ, 1976, № 2, с.ЮЗ-ПО.

36. Мусхелишвили Н.И. Некоторые основные задачи математической теории упругости. М», Наука, 1966, 707 с,

37. Агдсъи J.H. Energy theorems and stzvctutalanalysis. -tJiiczaft £лд." volp. /25~-/Я

38. Couzant R. Variational mettiods /ог tfie solution of p го Hems o^ e^uililzium and vi£ zaiions.1. Bale. Amez. . Soc. 23.

39. Корнеев В.Г. Схемы метода конечных элементов высоких поряд ков точности. Л., Изд-во Ленинградского ун-та, 1977, 208 с.

40. Корнеев В.Г., Пономарев С.Е, Решение задач теории пласти -ческого течения методом конечных элементов* Шурн.вычисл. матем.и матем.физ,", 1977, т.17, № 2, с.437-453.

41. Одиноков В.И., Каплунов Б»Г. Расчет пластического течения цилиндрического тела при осесимметричной деформации. При -кладная механика, 1978, 14, № 5, с. 16-24.

42. Одиноков В.IT.,Макеранец Е.И. К расчету пластического течения полых криволинейных стержней конечной длины. Приклад -ная механика, 1981, 17, № 9, с.93-98.

43. Каплунов Б.Г., Одиноков В.И. Силовые параметры горячего обратного выдавливания. Изв.вузов Черная металлургия, 1982, № 4, с.55-58.

44. Гвоздев А.А. Расчет несущей способности конструкций по методу предельного равновесия. Стройиздат, М., 1949,280с.

45. Одиноков В.И», Бубнов Э»А», Песков А.В» Численная схема расчета плоских криволинейных конструкций в условиях пластической деформации, В сб.Механика деформируемых сред» Куйбышевский госуниверситет, 1978, с»148-157.

46. Одиноков В.И,, Песков А.В», Лехов О.С,, Химич Т.Г. Теоретическое исследование процесса резания металла в горячем состоянии. В сб,Механика деформируемых сред. - Куйбышев ский госуниверситет, 1981, с.83-86.

47. Одиноков В.И., Лехов О.С., Химич Т.Г., Песков А.В. К воп -росу об определении усилий резания металла в горячем со -стоянии. Изв.вузов.Черная металлургия, 1981, № 8, с.142-145.

48. Одиноков В.И., Песков А.В. Расчет напряженно-деформированного состояния при кристаллизации балочной заготовки, отливаемой на МНЛЗ. Изв.вузов.Черная металлургия, 1983, № I, с.46-51.

49. Соколов Л.Д. Сопротивление металлов пластической деформа -ции. Металлургиздат, 1963.

50. Поздеев А.А., Тарновский В.И., Еремеев В.И., Баакашвили B.C. Применение теории ползучести при обработке металлов давлением. Изд-во "Металлургия", 1973, 192 с*

51. Песков А.В. Математическое моделирование некоторых техно -логических процессов, связанных с пластическим течением материала, Тезисы докладов конференции молодых ученых "Пути экономии черных и цветных металлов", Свердловск, 1982, с.76