Полуобратный метод анализа технологических операций обработки металлов давлением с использованием несимметричного тензора напряжений тема диссертации и автореферата по ВАК РФ 05.16.05, доктор технических наук Базайкин, Владимир Ильич

• Базайкин, Владимир Ильич
• доктор технических наук
• 2000

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

В Л Состояние вопроса : методы теоретического анализа.5

В.2 Состояние вопроса : моделирование конкретных процессов.10

В.З Обзор работ последних лет.14

В.4 Общая характеристика работы.25

1 Обоснование использования пространственного тензора напряжений при анализе конечного формоизменения.33

1.1 Некоторые замечания общего характера.33

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

1.3 Объективность тензоров, характеризующих локальную ротационную деформацию.40

1.4 Вопросы локального силового состояния, сопряжённого с кинематикой, допускающей ротационное деформирование.42

1.5 Адекватность моделей деформирования с полузаданной геометрией деформирования.45

1.6 Выводы по разделу 1 .49

2 Моделирование процесса волочения круглого профиля.52

2.1 Кинематика осесимметричного течения в зоне волоки.52

2.2 Напряжения, уравнения движения и определяющие соотношения.65

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

2.4 Решение задачи для пластических напряжений в локальном пространственном базисе.78

2.5 Аналитическое ( приближённое ) решение задачи для напряжений в пространственных координатах с учётом пластичности и вязкости материала.85

2.6 Напряжения, действующие в очаге деформации и на контакте металл-матрица .92

2.7 Приложение результатов к электростимулированию процесса волочения.95

2.8 Выводы по разделу 2.99

3 Моделирование операций холодной высадки стержневых изделий.101

3.1 Методологические аспекты теоретического анализа.101

3.2 Моделирование истечения материала в зоне деформирования.103

3.3 Определение поля тензора напряжений и усилий на пуансонах для идеального жёсткопластического материала заготовки.111

3.4 Напряжения и усилия высадки с учётом вязкости материала.124

3.5 Элементы локального вихревого движения в зонах деформации.132

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

3.7 Сопоставление теоретических распределений макродеформаций с данными экспериментальных исследований на субструктурном уровне.140

3.8 Выводы по разделу 3.144

4 Моделирование процесса холодной продольной прокатки листов.148

4.1 Построение деформирующего отображения.148

4.2 Определение полей скоростей и кинематических тензоров.156

4.3 Задача для напряжений при прокатке идеального жёсткопластического материала.165

4.4 Расчёт силовых параметров прокатки.180

4.5 Выводы по разделу 4.

5 Исследование способов удаления окалины с поверхности термообработанной проволоки.

5.1 Остаточные напряжения после калибровки.

5.2 Миграция водорода в процессе и после травления проволоки в кислотной ванне до калибровки .

5.3 Распределение водорода после калибровки.

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

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

5.6 Выводы по разделу 5.

Заключение.

Основные выводы.

Список литературных источников.

ВВЕДЕНИЕ

В.1. Состояние вопроса : методы теоретического исследования

Новые технологии обработки металлов давлением (ОМД) нуждаются в возможно более точном описании кинематических и силовых характеристик процесса деформирования , так как наблюдаются тенденции и к усложнению формы очага деформации и условий на его границах , и к переходу к материалам с особыми свойствами. Одно из направлений продвижения к такому описанию заключается в сближении методов механики деформируемого твёрдого тела с методами теории ОМД. Имеется в виду переход от приближённых уравнений равновесия , приспособленных для конкретных процессов , к балансовым равенствам в дивергентной форме. Соответственно от усреднённых скоростей , напряжений , деформаций необходимо переходить к определению тензорных полей этих величин. В результате обеспечивается учёт локальных особенностей характеристик процесса в любой части очага деформации.

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

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

II) Постановка задач с частично предугадываемым внутренним состоянием [3] , где априори задаются видом распределения некоторых внутренних параметров ( скоростей , напряжений ). Сами задачи формулируются в классическом виде с применением симметричного тензора напряжений. Известны решения задач для кручения , некоторых осесимметричных состояний , частных случаев изгиба. Как правило , здесь используется тензор малых линейных деформаций , мало пригодный для описания процессов ОМД. Предпринятые в работах последнего времени попытки энергетического обоснования полуобратных методов [4] подразумевают такие уровни пластических деформаций , которые не приводят к заметным изменениям формы. В то же время установлено [5] , что представления о выпуклости поверхности текучести совместимы с конечностью формоизменения только в случае идеального жёсткопластического тела. Этот факт был впоследствии учтён при моделировании конкретных технологических процессов [6,7].

III) С учётом последних замечаний в задачах плоской деформации ( при анализе ковки , выдавливания , прокатки ) давно и успешно применяется метод линий скольжения. При больших пластических деформациях большая часть их энергии выделяется в виде теплоты. Возникающий градиент температур обусловливает наличие термических напряжений , которые здесь не учитываются. При течении материала напряженное состояние локально определяется постоянным напряжением сдвига и переменным гидростатическим давлением. Уравнения равновесия и условие текучести переходят в интегралы Генки. Если граничные условия формулируются только для напряжений , то задача для напряжений является статически определимой. Таким образом , решение задачи по методу линий скольжения может сопрягаться с различными кинематиками , для выбора единственного решения необходимо применение критерия минимума мощности деформирования. гу) Инженерный метод расчёта , направленный на оценку нагрузки на инструмент , необходимой для пластического деформирования материала , и основывающийся на экстремальных принципах теории пластичности [3] -- метод верхней оценки такой нагрузки. Погрешность метода определяется степенью соответствия действительного и принятого в расчёте кинематически возможного поля скоростей. В практической реализации метод часто предполагает разбиение очага деформации на движущиеся друг относительно друга жёсткие блоки. Их конструкция определяется формой поверхности инструмента (матрицы). у) Метод упругих решений в двух формах : в. форме дополнительных нагрузок и в форме переменных параметров упругости [54]. Решение исходной нелинейной задачи сводится к решению последовательности линейных упругих задач с локально изменяющимися модулем упругости и коэффициентом Пуассона , причём начальное приближение является линейной задачей для однородно упругого тела. Найденная на каждом предыдущем этапе приближённая зависимость между интенсивностями напряжений и деформаций служит для определения коэффициентов упругости последующего этапа. Ограниченность этого метода заключается в том же недостатке , который присущ в целом деформационной теории пластичности : он применим в случаях незначительного изменения формы обрабатываемого материала ( "помнящего" начальную форму).

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

Однако интерес к возможным следствиям , вытекающим из факта существования поверхностных микромоментов-пар , привёл к необходимости использовать пространственный тензор напряжений общего вида ( несимметричный ) [11,12]. Было показано , что антисимметричная часть тензора напряжений компенсируется поверхностной плотностью моментов-пар , а сама эта плотность энергетически сопряжена с тензором кривизны-кручения. Эти обстоятельства потребовали более широкого привлечения аппарата дифференциальной геометрии при теоретическом изучения обычного и ротационного деформирования , например , в объёме [13].

Многие исследователи вспомнили работу братьев Коссера начала века [14] , в которой рассматривался упругий материал с необычным свойством : каждая его материальная точка могла не только смещаться, вызывая обычную деформацию растяжения-сжатия и сдвига окрестности , но и вращаться вокруг собственной оси. Физический смысл ротационной свободы не рассматривался , она просто постулировалась. Мерой ротационной деформации являлся тензор кривизны-кручения , сопряжённый тензору плотности моментов-пар ( моментных напряжений ). Так как евклидовость материального континуума заставляет тензор кривизны-кручения быть нулевым , то был сделан вывод о связи поля поверхностных моментов-пар с некоторой континуальной неевклидовой структурой , сосуществующей с материальным континуумом , именно - со структурой микродефектов внутри материального континуума. Теория локальных проявлений этих дефектов в виде дислокаций и дисклинаций была уже достаточно развита [15,16].

Так возник следующий этап анализа , на котором разрешённые несимметричным тензором напряжений моменты-пары были привязаны к континуальному распределению структуры микродефектов , находящихся и развивающихся в условиях макродеформации некоторого объёма. Для анализа параметров этих двух сосуществующих сред привлечён современный теоретический аппарат : лагранжианы и гамильтонианы [17,18] , калибровочные поля [19]. Естественно , здесь понадобилась техника дифференцирования тензорных полей достаточно высокого уровня [20,21]. Следует подчеркнуть , что ротационное деформирование и соответствующее ему силовое описание связывается не с геометрией деформирования , а с дефектной структурой. Далее , в этих теориях для возможности использования энергетических потенциалов необходимым элементом описания являются упругие макродеформации , а уровень пластических (дефектных) микродеформаций по умолчанию имеет порядок упругих деформаций.

В работах [22 — 24] предприняты попытки описания локального ротационного деформирования на основе антисимметричной части градиента пространственного поля скоростей ; кинематика течения металла в очаге деформации при волочении и высадке была частично задана. Вопрос о том , какой кинематический тензор является сопряжённым по мощности диссипации энергии деформирования тензору моментных напряжений , не ставился. При этом были использованы обоснования применения полуобратного метода , изложенные в [25 - 29] , а также результаты обсуждения правомочности применения общих термодинамических принципов к моделированию необратимого конечного формоизменения и роли при этом упругих деформаций [30 - 35].

В.2. Состояние вопроса : моделирование конкретных процессов

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

Опираясь на классические результаты по анализу процесса волочения [36 - 39], связанные с именами И.Л.Перлина, В.Л.Колмогорова, П.И.Минина, С.И.Губкина, Н.З.Днестровского, Н.Г.Решетникова, было предложено деформирующее отображение [40], определяющее конечное формоизменение локально несжимаемого материала ; использован известный математический аппарат [41 - 43]. В полном объёме предложенный анализ волочения представлен в [44].

В [45] предложена модификация определяющего соотношения для вяз-копластического тела, учитывающая влияние внешнего энергетического воздействия (стимуляции волочения импульсами электрического тока) на уровень вязких напряжений. При этом использована концепция определяющего соотношения как зависимости между инвариантами тензоров кинематических и напряжений, подчиняющейся аксиомам механики сплошной среды (локальности, объективности,.) [46]. Попытка применения ослабленного условия текучести Мизеса [47] потребовала применения метода конечных разностей [48] для расчёта напряжений. Относительные веса квазистатических пластических, вязких и чисто динамических напряжений определены в [49]. В [50] оценена роль чисто ротационного микродеформирования. Его мощность не превышает 1 - 2% общей мощности деформирования.

Обсуждения имеющихся данных о методах и результатах исследований процессов холодных осадки и высадки представлены в работах Г.А.Навроцкого, Г.А.Смирнова-Аляева, Ю.А.Алюшина, ГЛ.Гуна, М.В.Сторожева, Е.А.Попова, Я.М.Охрименко, В.А.Тюрина, В.Г.Паршина, О.С.Железкова, Б.В. Кучеряева [51 - 65]. Эти данные позволили заключить об отсутствии моделей связанной кинематики многопозиционной высадки даже в случае осевой симметрии процесса. Далее, не имеются данные о полях напряжений и конечных деформаций, возникающих при высадке, обычно используются усреднённые по объёму их оценки. Исходя из известных и собственных экспериментальных данных, для двухпозиционной высадки головки болта (винта) была частично задана геометрия деформирования, а из условия локальной несжимаемости материала заготовки определена кинематика процесса. Необходимо подчеркнуть, что в отличие от волочения высадка - нестационарный процесс. Тем не менее, используя в качестве временеподобной координаты текущую координату положения границы зон деформации, были определены история пластических и вязких составляющих напряжений и распределения компонент накопленного тензора конечных деформаций [66,67]. При этом была использована та же схема вязкопластического тела Бингама, что и при волочении [68]. Учёт вязких свойств деформируемого материала позволил показать зависимость напряжений и усилий от скорости движения пуансонов [69]. Для определения доли вязких составляющих в общих значениях напряжений использовано представление коэффициента вязкости в виде произведения второго инварианта девиатора тензора напряжений на время запаздывания сдвига в микромеханизмах деформации. Это представление ранее успешно использовалось при моделировании напряжений в процессе волочения [70].

Основные результаты теоретического анализа геометрии и кинематики деформирвания при продольной прокатке полос, а также формулирование и решение задач для действующих напряжений связаны с именами А.И.Цели-кова, И.Я.Тарновского, И.М.Павлова, В.С.Смирнова, В.Л.Колмогорова, Г.Я.Гуна, В.Н.Выдрина, В.П.Севрденко, Е.С.Рокотяна, Д.Бленда, Г.Форда [39,71 - 85]. Однако несмотря на то, что многие исследователи рассмотрели такие тонкие факты как внеконтактное деформирование [82,83], закон уширения по очагу деформации [73,74], площадь контактной поверхности [78], картины неоднородности деформирования и полей скоростей в очаге деформации изучены недостаточно. Среди наиболее важных достигнутых результатов по геометрии и кинематике продольной прокатки - вывод о том, что так называемый геометрический фактор и неодинаковая вытяжка горизонтальных слоев может определённым образом изгибать исходно прямолинейные поперечные сечения [79] ; следует указать также на представления о форме и положении критической поверхности [74,83,85]. Неоходимо обратить вримание на вывод об отставании полосы ( по скоростям) практически по всей контактной поверхности, сделанный в [82,83].

Силовое описание процесса прокатки основывается на решениях стати-тически определимых задач для приближённого уравнения равновесия в напряжениях [71,74] или на балансе мощностей : формоизменения, контактного трения, на бочке валка, динамической [78]. Как в том, так и в другом случае используются усреднённые величины (касательные напряжения выражаются через неопределённый коэффициент трения, при определении мощностей объёмные интегралы вычисляются по теореме о среднем значении). Надо отдельно остановиться на роли трения в описании процесса. Сила трения привязывается к нормальному давлению соотношениями Амонтона-Кулона или Зибеля через коэффициент трения, который является внешним по отношению к задаче о напряжениях параметром. Поэтому любые выводы о напряжениях в очаге деформации, об усилиях и мощности прокатки зависят от распределения коэффициента трения вдоль контактной поверхности. Но в такой постановке любые выводы о закономерностях трения, а значит, и о напряжениях, основываются на экспериментальных данных, которые часто трактуются неоднозначно. Этим обстоятельством обусловлен пристальный интерес к контактному трению при прокатке в работах [86 - 90].

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

Предлагаемый в данной работе подход к моделированию кинематики и напряжений был применён для задачи, призванной решить конкретную технологическую проблему. Оценивались методы удаления окалины с проволоки в технологической последовательности операций над ней : волочение -- термообработка - удаление окалины - калибровка под последующую высадку. Способ удаления окалины - кислотный или механический ( дробеструйный) - оказался тесно связанным с остаточными напряжениями после калибровки ; последняя в работе рассматривается как волочение с малыми обжатиями (и, следовательно, усилиями). При кислотном удалении в проволоке за счёт диффузии атомарного водорода из кислотной ванны возникает повышенная концентрация водорода. Последующая калибровка перераспре-ляет его таким образом, что в испытаниях на осадку заготовки показывают значительное снижение предельной пластичности [91,92]. Воспользовавшись принципом упругой разгрузки калибрующейся проволоки, опираясь на известные результаты по упругому деформированию стержней и методы математической физики, соотнесясь с данными по диффузии водорода в металлах [93 - 105], авторы работы [106] построили модель влияния водорода на технологическую пластичность стержневых заготовок и негативно оценили кислотный метод удаления окалины. Опираясь на общепринятые воззрения на фазовые превращения водорода в металлах [107,108], .авторы работы [109] предложили механизм воздействия остаточных макронапряжений на функционирование "ловушек" атомарного водорода.

Анализ возможного напряжённого состояния проволоки после дробеструйной обработки проведён в форме моделирования этого состояния остаточными эффективными напряжениями, создающимися в конструкции из упругого цилиндрического стержня, содержащегося в пластической оболочке. Используя некоторые решения упругих задач [110,111], в [112] была построена модель эффективных напряжений в зоне дробеструйной обработки, использующая принципы разграничения упругих и пластических напряжений и деформаций и упругой разгрузки [113 - 116]. Показано, что распределение указанных остаточных напряжений является благоприятным для смягчения схемы напряжённого состояния при последующей калибровке. При этом обсуждены существующие подходы к истолкованию самого понятия остаточных напряжений [117-121].

В.З. Обзор работ последних лет а) Модели деформируемых сред и новые подходы к решению задач.

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

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

Для теоретического анализа процессов ОМД работа [125] носит обзорный характер. Рассмотрены типы определяющих соотношений в связи с их использованием в принципах виртуальной мощности и виртуальных скоростей и напряжений. Высказывается утверждение, что пренебрежение упругими деформациями, являясь грубым приближением, и вызывает необходимость применения методов верхней и нижней оценок силы деформирования. Ссылаясь на свои результаты [126], автор предлагает новый тип определяющих соотношений, имеющих значение для ОМД. Вывод этого соотношения основан на разложении в текущем лагранжевом базисе именно ко-вариантных компонент тензора полных деформаций на упругие и пластические составляющие. Предложена модификация ассоциированнного закона течения, в которой скорость пластической деформации зависит не только от девиатора напряжений, но и от скорости его изменения. Полагая, что среднее напряжение надо вычислять не как первый инвариант тензора напряжений, а через относительное изменение упругого объёма, автор получает отличное от принятого .понятие девиатора напряжений.

Претендует на фундаментальность подход В.Л.Колмогорова с сотрудниками, показанный в [127]. Задаваясь полем перемещений частиц, находят связи между отсчётными и актуальными положениями частиц, строят базисы сопутствующей системы координат, в которых рассматривают ко - и конт-равариантные компоненты векторов скоростей и ускорений, тензора скоростей деформации. Из шести компонент тензора напряжений три (недиагональные) выбирают произвольно, но с учётом граничных условий, три нормальные компоненты получают из уравнений движения, рассматриваемых в сопутствующем базисе. Из последних получают решения, содержащие неизвестные коэффициенты-функции. Подставляя решения в функционал, переводят задачу о его экстремуме в задачу для обыкновенного дифференциального уравнения, решаемого одним из численных методов. Для случая осадки параллелепипеда отмечено, что существенное снижение уровня напряжений в конце деформации отражает роль инерционных нагрузок.

В [128] представлен алгоритм расчёта ускорений и инерционных сил при трёхмерной деформации. Он основан на использовании принципа наименьшего принуждения Гаусса в полях скоростей из жёстких блоков для оценки ускорений и динамических параметров процесса деформации [129], обоснование такого подхода изложено в [130]. В этой работе указывается, что для безвихревых полей скоростей в несжимаемом материале интенсивность деформации сдвига определяется только дивергенцией вектора ускорения ( на примере однородной осадки ). В то же время уравнения движения предполагают отсутствие ускорения в условиях однородной деформации. Далее, ускорения могут претерпевать разрывы даже в зонах с непрерывными скоростями. Эти затруднения преодолеваются методом верхней оценки с учётом ускорений. Предложена методика расчёта мощности сил инерции. В [128] эта методика была обоснована для использования при верхней оценке мощности деформирования и усилий на основе кинематически возможных полей скоростей из непрерывно деформируемых областей. При линейном распределении скоростей внутри каждой области компоненты вектора ускорения определяются конвективными и локальными составляющими определённого вида. Мощности инерционных сил находятся с учётом разрывности поля скоростей. Отмечена трудность реализации локальной производной по времени, которая должна быть разрешена для каждого конкретного процесса. Общие соотношения и алгоритмы для областей прямоугольного типа изложены в [131], для осесимметричных областей - в [132]. Показано, что при малых скоростях деформирования (около 1 м/с) инерционными силами можно пренебречь, но мощность, определённая только по полю скоростей, и мощность ускорений сопоставимы при скоростях, больших 50 м/с.

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

Следует также отметить проект автоматизации метода линий скольжения и разработанный на его основе пакет программ МЕЛИСА [134], позволяющий : рассчитывать напряжённо-деформированное состояние ; определять усилия деформации, поля деформации, поля показателя повреждённости. Приведены примеры плоской деформации. б) Исследования процесса волочения. В первую очередь необходимо отметить статью В.И.Тарновского [135]. Анализируются деформированное состояние и усилия волочения полосы из идеального жёсткопластического материала через клиновидную волоку. Касательные напряжения трения на контактной поверхности определяются как предел текучести на сдвиг, умноженный на коэффициент, учитывающий условия трения. Факт радиальности течения в плоскости сечения полосы предопределяет использование полярной системы координат. Конструкция предполагаемой скорости истечения весьма жестка : множителем при постоянной входной скорости служит отношение функции, зависящей только от полярного угла, к полярному радиусу. Напряжения в полярной системе координат в тригонометрической форме" входят в систему уравнений равновесия типа [113]. В соответствующей форме используется условие текучести Мизеса. Получено замкнутое аналитическое, но неявное описание напряжённо-деформированного состояния. Показано, что относительная удельная сила волочения линейно зависит от логарифма вытяжки с коэффициентом, определяемым углом волоки и условиями трения.

Условиям трения при волочении посвящена и работа [136]. В ней дано аналитическое описание изотермического процесса захвата сухой технологической смазки (мыльного порошка) гладкими поверхностями проволоки и одинарной волоки. Учитываются многократность волочения и влияние тоннельного эффекта в слое порошка смазки вокруг колеблющейся заготовки перед первой волокой при значительной скорости волочения. Определена толщина слоя смазки при наличии этих факторов. Использован аппарат гидродинамики вязкой жидкости.

В продолжение этих исследований в [137] исследованы факторы, влияющие на нормированный коэффициент проявления тоннельного эффекта (имеющего негативный характер) : скорость волочения, влажность мыла, его дисперсность. Показано, что тоннельный эффект проявляется максимально при уменьшении на 25 - 50 % расчётного значения толщины слоя сухой смазки. в) Теория осадки и высадки. В работе [138] исследуется осадка параллелепипеда из линейно вязкого материала. Поле скоростей определяется через две неизвестные функции, зависящие от относительной координаты в направлении осадки. Так как материал вязкий, минимизируемая полная мощность деформирования (без трения, материал прилипает к плитам) зависит только от второго инварианта тензора скоростей деформации. Минимум функционала энергии приводит к системе двух обыкновенных дифференциальных уравнений для функций, входящих в выражения компонент скоростей. Естественные начальные условия позволяют полностью определить поля скоростей. Давление осадки находится из баланса мощностей внешних и внутренних сил (с постоянной величиной сопротивления деформации). Таким образом, используется метод верхней оценки. Для высоких цилиндров отмечен факт немонотонности изменения по высоте скорости поперечной деформации, что трактуется как описание двойной бочкообразности.

Серия статей М.Я.Бровмана посвящена осесимметричной деформации. В [139] рассмотрен класс кинематически допустимых скоростей для задачи о вдавливании с постоянной скоростью в идеальное жёсткопластическое полупространство кольцевого штампа. Границы пластических областей аппроксимируются отрезками парабол. Определяются мощности формоизменения в пластически« зонах, среза на границах областей. В качестве варьируемых параметров избраны размеры пластических областей. Через их значения выражаются три функции, определяющие среднее давление на поверхности штампа. Получены формулы для верхней и нижней оценок усилий вдавливания. В работе [140] проанализированы два случая для цилиндрического пуансона. В первом из них материал проскальзывает вдоль пуансона, во втором материал движется совместно с пуансоном без скольжения. По методу верхней оценки определены мощности формоизменения, среза, трения под пуансоном. Подчёркнуто существенное значение ~ил трения (до 20 % общей мощности). Построено упрощенное статически допустимое поле напряжений, получена нижняя оценка коэффициента напряжённого состояния.

Интерес автора вышеуказанных работ к возможно более общим постановкам отражён в [141]. Здесь задачи рассматриваются в криволинейных ортогональных координатах, используются уравнения движения, но процесс осе-симметричного течения - стационарный, одно из семейств координатных линий служит линиями тока. Связь между инвариантами тензоров напряжений и скоростей деформации принята в виде степенной функции с множителем - пределом текучести на сдвиг. Это позволяет варьировать определяющие соотношения от идеального жёсткопластического типа до линейно-вязкого. Фактически применён полуобратный метод ; применён симметричный тензор Пиолы-Кирхгофа в локальном базисе деформированной конфигурации ; к сожалению, компоненты этого тензора в примере рассматриваются как физические напряжения. Отмечено, что динамические составляющие напряжений могут достигать 23 % общих их значений при скоростях течения около 30 м/с, при этом значения касательных напряжений сохраняются на статических уровнях.

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

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

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

Среди работ, анализирующих конкретные операции, необходимо отметить [145,146]. В [145] с использованием вариационного метода в дискретной постановке разработана методика определения энергосиловых параметров процесса высадки головок стержневых изделий вращающимся валком. Определены усилия деформирования и момент прокатки. Материал заготовки представлялся жёсткопластическим с нелинейным упрочнением (по Г.А.Смирнову-Аляеву). В работе [146] на основе расчёта силовых параметров с кинематически возможными полями скоростей из жёстких блоков при различных расположениях заусенцев, образующихся в процессе высадки головки болта, указаны рациональные схемы деформирования. Во внимание принимались возможные усилия на заключительной стадии процесса.

Интересна для возможного теоретического анализа новая конструкция штампа [147] для реализации высадочных операций при формировании головок стержневых изделий. Конструкция обеспечивает плавный наборный переход к конечной форме. Детали поддерживающих элементов образуют вращательную, поступательную и клиновую пары. Конструкция особенно эффективна при больших длинах высаживаемой части заготовки. г) Работы по холодной продольной прокатке. В них ставятся и решаются специальные задачи теории прокатки. Авторы [148] , пользуясь гипотезой плоских сечений, задаются распределением вдоль оси прокатки продольных составляющих скоростей течения в очаге деформации в виде полиномов второй, третьей или пятой степеней (без обоснования). Поставив и решив для идеального жёсткопластического материала соответствующую вариационную задачу (по методу верхней оценки), авторы показали, что степень полинома растёт с увеличением степени обжатия. При обжатиях в диапазоне 1,2. 1,8 оптимально описание полиномом третьей степени, в этом случае можно пренебречь влиянием коэффициента опережения.

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

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

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

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

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

Автор [155] предлагает инженерную методику оценки желобчатости полос при прокатке-волочении - предельном случае прокатки с кинематической несимметрией. Классическая постановка : приближённое уравнение равновесия, закон трения по Зибелю, прямолинейность линии контакта, статически определимая задача - позволяет разрешить напряжённо-деформированное состояние в приконтактных к ведущему и к ведомому валкам слоях полосы. Интегрированием приращений деформаций вдоль осей валков находятся различные конечные ширины приконтактных слоёв. По их значениям и находится радиус желобчатости.

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

Задавая искажения семейства прямых, параллельных оси прокатки, в очаге деформации, исходя из условия локальной несжимаемости, автор [157] построил поля тензора скоростей деформации и тензор-спина. На основе этой кинематики в [158] определено поле пространственного тензора напряжений, распределения касательных и радиальных контактных напряжений ; материал полосы в очаге деформации - идеальный жёсткопластический. Исползуя результаты двух последних работ и перейдя к традиционным угловым переменным, в [159] по распределению контактных касательных напряжений определён крутящий момент на валке. д) Анализ некоторых проблем, сопряжённых с технологией волочения.

Как отмечается в [160] , при производстве осесимметричных изделий в них после деформирования реализуется схема плоского деформированного состояния (отсутствуют осевые деформации). В силу симметрии уравнения равновесия для упругого материала упрощаются. Методом Ритца решается вариационная задача для потенциальной энергии упругих деформаций, обусловленных остаточными напряжениями. Фактически решается упругая задача с заданием нулевых радиальных напряжений на поверхности изделия ; не учитывается схема предыдущего пластического нагружения.

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

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

25

В.4. Общая характеристика работы.

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

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

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

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

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

Настоящая работа проводилась в соответствии с межвузовской программой НИР "Металл", раздел 09.16 ; координационным планом НИР АН СССР на 1986 - 1990 гг., раздел 13.2 ; Программой фундаментальных исследований "Повышение надёжности систем машина - человек - среда" АН СССР на 1989 - 2000 гг., раздел 3 ; Федеральной программой "Интеграция" на 1997 - 2000 гг.; грантами, Министерства образования РФ по фундаментальным проблемам металлургии (1992 - 2000).

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

Для этого ста- вились и решались следующие задачи.

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

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

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

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

Основные положения , выносимые на защиту.

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

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

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

4. Разделение чисто пластических, вязких и динамических составляющих напряжений и оценка роли двух последних типов напряжений в операциях обработки металлов давлением.

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

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

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

Научная новизна.

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

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

3. Выполнено разделение пластических, вязких и динамических составляющих напряжений при волочении и произведена оценка их относительных долей .

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

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

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

Практическая значимость. Полученные в работе результаты позволяют не только определить силовые интегральные параметры процессов ОМД , но и проанализировать локальные особенности распределений напряжений. Сопоставление этих особенностей при симулировании волочения импульсами электрического тока с данными оптической и электронной микроскопии позволило выделить механизмы микродеформирования, чувствительные к токовому воздействию , и предложить технологические рекомендации по внедрению электростимуляции в процесс волочения (ОАО "ЗСМК", г. Новокузнецк, имеется акт об использовании результатов диссертации).

Расчёты усилий на пуансонах при двухпозиционной высадке головок болтов и накопленных в них деформаций позволили разработать и внедрить оптимальные схемы формоизменения при производстве крепёжных изделий (ОАО "Автонормаль", г. Белебей, имеется акт об использовании результатов диссертации).

Результаты моделирования остаточных напряжений, образующихся после травления проволоки в кислотной ванне и последующей калибровки , показали их негативное влияние на предельную пластичность при осадке, подтверждённое данными испытаний в условиях ОАО "ЗСМК", г. Новокузнецк.

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

Научные результаты, полученные на основании полуобратного метода анализа кинематики и напряжений, рекомендованы к использованию в специальных курсах, читаемых в Сибирском государственном индустриальном университете, в Томском государственном архитектурно-строительном университете. УМО по металлургии Министерства образования РФ монографии "Электростимулированное волочение : структура и анализ" (авторы В.Е.Громов, В.Я.Целлермаер, В.И.Базайкин) присвоен гриф "Учебное пособие".

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

I Всесоюзной конференции "Действие электромагнитных полей на пластичность и прочность металлов и сплавов", Юрмала, 1987 г.; Всесоюзном семинаре "Пластическая деформация материалов в условиях внешних энерге-тичесих воздействий", г.Новокузнецк, 1991 г.; III Международной конференции "Прочность и пластичность материалов в условиях внешних энергетических воздействий", г.Новокузнецк, 1993 г. ; Российской научно-технической конференции "Новые материалы и технологии", г.Москва, 1994 г.; III Международной конференции "Действие электромагнитных полей на пластичность и прочность материалов", г.Воронеж, 1994 г.; I Международной конференции "Актуальные проблемы прочности", г.Новгород, 1994 г.; XIV Международной конференции "Физика прочности и пластичности материалов", г.Самара, 1995 г.; III Международной конференции "Структурно-морфологические основы модификации материалов .", г.Обнинск, 1995 г.; IV Международной конференции "Прочность и пластичность материалов в условиях внешних энергетических воздействий", г.Новокузнецк, 1995 г.; I конференции "Материалы Сибири", г.Новосибирск, 1995 г.; IV Международной конференции "Компьютерное конструирование перспективных материалов и технологий",г.Томск,

1995 г.; VII Международной конференции "Структура дислокаций и технологические свойства металлов и сплавов", г.Екатеринбург, 1996 г. ; II Сибирском конгрессе по прикладной и индустриальной математике, г.Новосибирск,

1996 г.; Международной конференции "Математические модели и численные методы механики сплошных сред", г.Новосибирск, 1996 г.; III Международной школе-семинаре "Эволюция дефектных структур в конденсированных средах", г.Барнаул, 1996 г.; IV Международной конференции "Действие электромагнитных полей на пластичность и прочность материалов", г.Воронеж, 1996 г.; Симпозиум "Синэргетика, структура и свойства материалов, самоорганизующиеся технологии" г.Москва, 1996 г.; Научно-технической конференции "Структурная перестройка металлургии : экономика, экология, управление, технология", г.Новокузнецк, 1996 г.; IV Мехсдународной конференции "Структурные основы модификации материалов .", г.Обнинск, 1997 г.; Международной конференции "Всесибирские чтения по математике и механике", г.Томск, 1997 г.; Научно-практической конференции "Современные проблемы и пути развития металлургии", г.Новокузнецк, 1997 г.; I Международном семинаре "Актуальные проблемы прочности", г.Москва, 1997 г. ; V Международной конференции "Водородное материаловедение и химия гидридов металлов", г. Ялта, 1997 г.; III Сибирском конгрессе по прикладной и индустриальной математике, г.Новосибирск, 1998 г. ; II Международной конференции "Водородная

31 обработка материалов", г.Донецк, 1998 г.; Международной научно-практической конференции "Современные проблемы и пути развития металлургии", г.Новокузнецк, 1998 г.; межвузовской научной конференции "Численно-ана-литческие методы решения краевых задач", г.Новокузнецк, 1998 г.; VI Международной научно-технической конференции "Актуальные проблемы материаловедения", г.Новокузнецк, 1999 г.; 11 и 12 зимних школах по механике сплошных сред , г. Пермь , 1997, 1999 гг.

Публикации. Результаты работы опубликованы в двух монографиях и более 80 других работах. Список основных из них (3 8 публикаций в научных журналах и в сборниках научных трудов) приведён в конце автореферата.

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

Заключение диссертации по теме «Обработка металлов давлением», Базайкин, Владимир Ильич

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

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

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

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

225 ограниченных значениях вязких компонент (до 10 % общих значений) имеются объёмы со сжимающими всеми компонентами тензора напряжений (в отсутствии противонатяжения).

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

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

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

9. С определением поля тензора напряжений найдены: распределения его компонент в очаге деформации; векторы полного, касательного и нормального напряжений на дуге контакта полосы с валком , коэффициент трения ; крутящий момент и давление на валки ; положение точки нейтрали на дуге контакта зависит от значений переднего и заднего натяжений.

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

ЗАКЛЮЧЕНИЕ

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

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

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

Изгибное деформирование означает депланацию площадки с некоторой неподвижной осью в ней.

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

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

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

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

Список литературы диссертационного исследования доктор технических наук Базайкин, Владимир Ильич, 2000 год

1. Колтунов М.А., Кравчук A.C., Майборода В.П. Прикладная механика деформируемого твёрдого тела. М.: Высшая школа, 1983. - 349 с.

2. Клюшников В.Д. Математическая теория пластичности. М.: Изд - во МГУ 1979.-206 с.

3. Буи Х.Д. Введение в теорию обратных задач механики материалов (перевод с англ.). Караганда : Изд - во КарГУ, 1997. - 378 с.

4. Циглер Г. Экстремальные принципы термодинамики необратимых процессов и механика сплошной среды. М.: Мир, 1966. - 134 с.

5. Базайкин В.И., Закиров Д,М., Громов В.Е., Целлермаер В.Я. Напряжённо --деформированное состояние при высадке головки болта.// Изв.вуз.Чёрная металлургия. 1996.№12. с.28 30.

6. Годунов С.К. Элементы механики сплошной среды. М.: Наука, 1978. -303 с.

7. Трусдел К. Первоначальный курс рациональной механики сплошных сред. --М.: Мир, 1975.-592 с.

8. Ю.Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1978. - 736 с.11 .Кувшинский Е.В., Аэро Э.Л. Континуальная теория асимметрической упру-гости.Учёт "внутреннего" вращения.// ФТТ. 1963. т.5, в.9. с. 2591 2598.

9. Аэро Э.Л., Булыгин А.Н., Кувшинский Е.В. Асимметрическая гидромеханика.// ПММ. 1965. т.29, в.2. с. 138 145.

10. Схоутен Я.А. Тензорный анализ для физиков. M.: Наука, 1965. - 455 с.

11. Cosserat Е., Cosserat F. Theorie des Corps deformable. Paris : Hermann, 1909. -324 p.

12. Лихачёв В.А., Хайров Р.Ю. Введение в теорию дисклинаций. JL: Изд - во ЛГУ, 1975.- 183 с.

13. Ландау Л.Д., Лифшиц Е.М. Теория упругости. М.: Наука, 1965. - 202 с.

14. Попов В.Л. Моментные напряжения в кристаллической среде с дислокациями.// Изв.вуз.Физика. 1994.№8. с.117 — 118.

15. Попов В.Л. Динамическая теория пластических поворотов в кристаллах. //Изв.вуз.Физика. 1994.№8. с.119- 121.

16. Гриняев Ю.В., Чертова Н.В. Калибровочные теории пластической деформации в механике сплошных сред.// Изв.вуз.Физика. 1990.№2. с.36 50.

17. Лурье А.И. Нелинейная теория упругости. М.: Наука, 1980. - 512 с.

18. Победря Б.Е. Лекции по тензорному анализу. М.: Изд - во МГУ, 1986. --262 с.

19. Базайкин В.И., Громов В.Е., Целлермаер В.Я., Базайкина Т.В. О полярности деформируемого материала при волочении проволоки.// Изв.вуз.Чёрная металлургия. 1994. №12. с.21 -24.

20. Базайкин В.И., Громов В.Е. Перспективы применения полуобратного метода для анализа напряжений в процессах ОМД.// Перспективы горно-металлургической индустрии. Новосибирск : Сибирские огни, 1999. - 371.

21. Эриксен Дж. Исследования по механике сплошных сред. М.: Мир, 1977. - 245 с.

22. Adkins J.E. A recinrocal property of the finite plane strain equations .//J. Mech. and Phys. Solids. 1958.№6. p.267 275.

23. Shield R.T. Inverse deformation results in finite elasticity.// Z. angrew. Math, und Phys. 1967. №18. p.490 500.

24. Гелей Ш. Расчёт усилий и энергии пластической деформации металлов. -М.: Металлургиздат, 1958. 289 с.

25. Carlson D.E., Shield R.T. Inverse deformation results for elastic materials.// Z.angrew. Math, und Phys. 1969. №20. p.261 263.

26. Gurtin M.E., Murdoch I. A continuum theory of elastic material surfaces.// Arch. Ration. Mech. And Analysis. 1975.№57. p.291 323.

27. Sewell M.J. On the calculation of potential functions defined on curved boundaries.// Proc. Roy. Soc. London. 1965. A286. p.402 411.

28. Wang C.C., Truesdell C. Introduction to rational elasticity. Leuden : Noordhof International, 1973.-320 p.

29. Эриксен Дж.Л. Некоторые проблемы устойчивости в нелинейной теории упругости.// Сб. "Успехи механики деформируемых сред". М.: Наука, 1975. -с.559 - 565.

30. Stolz С. Sur les equations generates de la dynamique des milieux continus anela-tiques.//Comptes R. Acad. Sci. 1991. №2. p.307 316.

31. Базайкин В.И., Громов B.E., Кравченко П.Е. Использование модели жёстко-пластического тела при анализе течения и напряжений в теории ОМД.// Изв. вуз.Чёрная металлургия. 1997. №6. с.57 59.

32. Минин П.И. Исследование волочения прутков и проволоки. М.: Машгиз, 1948.-324 с.

33. Павлов И.М. Теория прокатки. М.: Металлургиздат, 1950. - 430 с.

34. Перлин И.Л., Ерманок М.З. Теория волочения. М.: Металлургия, 1971. -- 446 с.

35. Джонсон У., Меллор Л. Теория пластичности для инженеров. М.: Машиностроение, 1979. - 566 с.

36. Базайкин В.И., Громов В.Е., Кузнецов В.А., Перетятько В.Н. Волочение круглого профиля с внешней энергетической стимуляцией.Определяющее соотношение и кинематика.//Изв.вуз.Чёрная металлургия. 1989.№6.с.51 54.

37. Шварц Л. Анализ.Том 1.-М.:Мир, 1972.-824 с.

38. Фрейденталь А., Гейрингер X. Математические теории неупругой сплошной среды. М.: Физматгиз, 1962. - 432 с.

39. Качанов Л.М. Основы теории пластичности М.: Наука, 1969. - 420 с.

40. Физика и механика волочения и объёмной штамповки (под ред. В.Е.Громова и Э.В.Козлова). М.: Недра, 1997. - 289 с.

41. Базайкин В.И., Громов В.Е., Кузнецов В.А., Перетятько В.Н. Волочение круглого профиля с внешней энергетической стимуляцией. Напряжения и анализ.// Изв.вуз.Чёрная металлургия. 1989. №8. с.76 80.

42. Жермен П. Курс механики сплошных сред. М.: Высшая школа, 1983. -- 399 с.

43. Базайкин В.И.,Громов В.Е., Полторацкий Л.М., Перетятько В.Н. Моделирование напряжений при волочении проволоки с токовой стимуляцией.// Изв. вуз.Чёрная металлургия. 1993. №2. с.33 36.

44. Годунов С.К., Рябенький B.C. Разностные схемы. М.: Наука, 1973. - 400 с.

45. Базайкин В.И., Громов В.Е., Целлермаер В.Я. Разделение пластических и вязких напряжений при волочении круглого профиля.// Вестник горно-металлургической секции АЕН РФ.Отделение металлургии. Вып. 1. 1994.с.93 98.

46. Базайкин В.И., Целлермаер В.Я., Громов В.Е., Кравченко П.Е. Анализ напряжённого состояния при волочении проволоки из материала с варьируемым определяющим соотношением.//Изв.вуз.Чёрная металлургия. 1996. №8. с.17 30.

47. Навроцкий Г.А. Кузнечно-штамповочные автоматы. -М.: Машиностроение, 1965. 418 с.

48. Кучерляев Ь.В., Потанков H.A. Механика сплошных сред (раздел : Математические методы решения задач ОМД). М.: МИСиС, 1992. - 164 с.

49. Ильюшин A.A. Пластичность. -М.: Изд. АН СССР, 1963. 270 с.

50. Александров A.B., Потапов В.Д. Основы теории упругости и пластичности. М.: Высшая школа, 1990. 398 с.

51. Смирнов-Аляев Г.А. Механические основы пластической обработки металлов. JI: Машиностроение, 1968. - 270 с.

52. Гун Г.Я. Теоретические основы обработки металлов давлением. М.: Металлургия, 1980. - 456 с.

53. Алюшин Ю.А. Теория обработки металлов давлением. Ростов-на-Дону : Изд-во РИСХМ, 1977. - 85 с.

54. Алюшин Ю.А. Исследование процессов обработки металлов давлением с помощью кинематически возможных полей скоростей. Ростов-на-Дону Изд-во РИСХМ, 1978. - 96 с.

55. Сторожев М.В., Попов Е.А. Теория обработки металлов давлением. М.: Машиностроение, 1977. - 422 с.

56. Охрименко Я.М., Тюрин В А. Теория процессов ковки. М.: Высшая школа, 1977.-294 с.

57. Холодная объёмная штамповка : Справочник (под ред. Г.Н.Навроцкого). -М.: Машиностроение, 1973. -273 с.

58. Паршин В.Г., Поляков М.1^., Железков О.С. Метод определения усилий холодной высадки болтов и винтов.//Бюл. ин-та "Черметинформация". 1975. №12. с.48 49.

59. Паршин В.Г., Железков О.С. Определение усилий холдной объёмной штамповки осесимметрических деталей.// Изв.вуз.Чёрная металлургия. 1980. №3 с.86 89.

60. Железков О.С., Артюхин В.И. Измерение усилий в процессе холодной штамповки головок стержневых изделий.// Межвуз.сб. "Теория и практика производства метизов". Магнитогорск, 1989. - с.97 - 101.

61. Кучеряев Б.В., Попов В.Д. Расчёт технологических параметров высадки головки гвоздя на кривошипно-шатунных автоматах.// Сб. "Пластическая деформация сталей и сплавов". -М.: МИСиС, 1996. с. 179 - 187.

62. Базайкин В.И., Громов В.Е. Об учёте вязкости материала при анализе напряжений в процессах холодной объёмной штамповки.// Тезисы 3-го Сибирского конгресса по прикладной и индустриальной математике. Новосибирск, июнь 1998. с.85.

63. Закиров Д.М., Базайкин В.И., Громов В.Е., Целлермаер В.Я. Напряжения и деформации при формировании головки болта.// Изв.вуз.Чёрная металлургия. 1997. №2. с.22 29.

64. Малинин H.H. Прикладная теория пластичности и ползучести. М.: Маши ностроение, 1975. - 398 с.

65. Громов В.Е., Целлермаер В.Я., Базайкин В.И. Электростимулированное волочение : структура и анализ. М.: Недра, 1996. - 160 с.

66. Смирнов B.C. Теория прокатки. М.: Металлургия, 1967. - 459 с.

67. Чижиков Ю.М. Моделирование процесса прокатки. М. : Металлургиздат, 1963.-362 с.

68. Целиков А.И. Основы теории прокатки. М.: Металлургия, 1965. - 247 с.

69. Целиков А.И., Гришков А.И. Теория прокатки. М.: Металлургия, 1970. -358 с.

70. Томлёно в А. Д. Теория пластического деформирования металлов. М.: Металлургия, 1972. -408 с.

71. Полухин П.И., Гун Г.Я., Галкин A.M. Сопротивление деформации металлов и сплавов. М.: Металлургия, 1976. - 487 с.

72. Колмогоров В.Л. Напряжения, деформации, разрушение. М.: Металлургия, 1970. - 230 с.

73. Выдрин В.Н., Федосиенко A.C., Крайнов В.И. Процесс непрерывной прокатки. М.: Металлургия, 1970. -455 с.

74. Павлов И.М. Теория прокатки, ч.1. -М.: Металлургиздат, 1950. 380 с.

75. Рокотян С.Е. Теория прокатки и качество металла. М.: Металлургия, 1981. 222 с.

76. Лаврик В.И., Савенков В.Н. Справочник по конформным отображениям. Киев : Наукова думка, 1970. 251 с.

77. Смирнов B.C., Кирицэ В.В. //Сб."Обработка металлов давлением".Л.: ЛПИ, 1965. Вып. 243. с.79-84.

78. Тарновский И.Я. и др. Деформация металла при прокатке. М.: Металлург-издат, 1956. -349-с.

79. Головин Ф. Прокатка.Часть 2. М.: ОНТИ, 1934. - 263 с.

80. Целиков А.И. О распределении деформации по сечению прокатываемого металла.// Труды ЦНИИТМАШ. 1955. Вып. 73. М.: Машгиз. 211 с.

81. Полухин П.И. Николаев В.А., Полухин В.П. и др. Контактное взаимодействие металла и инструмента при прокатке. М. - Металлургия, 1974. - 200 с.

82. Грудев А.П. О методах определения коэффициента трения при прокатке // Сб. "Прокатное и трубное производство". М.: Металлургиздат, 1958. с.47 59.

83. Грудев А.П. Современные методы исследования внешнего трения при прокатке.// Сб. "Теория прокатки".М.: Металлургиздат, 1962. с. 118 — 127.

84. Грудев А.П. Коэффициент трения на двух стадиях процесса прокатки.// Сб. "Обработка металлов давлением". Вып. 4. М.: Металлургиздат, 1956.с.132 140.

85. Гончаров Ю.В., Прокофьев В.Н. Исследование экспериментальных методов определения коэффициентов внешнего трения при прокатке.// Труды Днепропетровского металлургического института.Вып. 49. М.: Металлургия, 1965. с.83 89.

86. Полторацкий Л.М., Громов В.Е., Чинокалов В.Я., Целлермаер В.Я. Влияние водорода на пластичность проволоки при холдной осадке.// Изв.вуз. Чёрная металлургия. 1991. №4. с.56 58.

87. Полторацкий Л.М., Громов В.Е., Базайкин В.И. Влияние остаточных напряжений на пластичность проволоки, полученной из травленной катанки. // Изв.вуз.Чёрная металлургия. 1994. №2. с.84 85.

88. Амензаде Ю.А. Теория упругости. М.: Высшая школа, 1971. - 286.

89. Кодес Е.С. Водород и несовершенства структуры металлов. М.: Металлургия, 1980.-221 с.

90. Гельд П.В.,Рябов P.A., Мохрачёва Л.П. Водород и физические свойства металлов и сплавов. М.: Наука, 1985. - 232 с.

91. Калачёв Б.А. Водородная хрупкость металлов. М.: Металлургия, 1985. --216 с.

92. Спивак Л.В., Скрябина Н.Е., Кац М.Я. Водород и механическое последействие в металлах и сплавах. Пермь : Изд-во Перм. унив-та, 1993, - 344 с.

93. Уманский Я.С., Скаков Ю.А. Физика металлов. М.: Атомиздат, 1978. -- 350 с.

94. Араманович И.Г., Левин В.И. Уравнения математической физики. М.: Наука, 1969.-286 с.

95. Джефрис Г., Свирлс Б. Методы математической физики. Вып. 3. М.: Мир, 1970.-343 с.

96. Кошляков М.С., Глинер Э.В., Смирнов М.М. Уравнения в частных производных математической физики. М.: Высшая школа, 1970. - 710 с.

97. Справочник по специальным функциям (под ред. М.Абрамовича и И.Стегана). М.: Наука, 1979. - 830 с.

98. Кузнецов Д.С. Специальные функции. М.: Высшая школа, 1965. - 420 с

99. Морозов А.Н. Водород и азот в стали. М.: Металлургия, 1968. - 280 с.

100. Физическое металловедение. Под ред. Р.Кана. М.: Мир, 1968. - 490 с.

101. Базайкин В. И., Кравченко П. Е., Целлермаер В.Я., Громов В.Е., Заки-ров Д.М. Влияние водорода на технологическую пластичность при волочении.// Изв.вуз.Чёрная металлургия. 1996. №8. с.41 -44.

102. Саменков В.А., Шильштейн С.Ш. Фазовые превращения водорода в металлах. М.: Наука, 1978. - 80 с.

103. Уманский Я.С. Карбиды твёрдых сплавов. М.: Металлургиздат, 1947. -324 с.

104. Закиров Д.М., Громов В.Я., Базайкин В.И., Кравченко П.Е.// Тезисы докладов 4-ой Международной конф. "Структурные основы модификации2.3^1материалов". Обнинск : МОПО РФ, 1997. с.85. 1 Ю.Лурье А.И. Теория упругости. - М.: Наука, 1970. - 939 с.

105. Папкович П.Ф. Теория упругости. М.: Оборонгиз, 1939. - 420 с.

106. Базайкин В.И., Громов В.Е., Закиров Д.М. Эффективные напряжения после удаления окалины с термически обработанной проволоки.// Изв.вуз. Чёрная металлургия. 1998. №10. с.55 58.

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

108. Аркулис Г.Э., Дорогобид В.Г. Теория пластичности. М.: Металлургия, 1987. - 351 с.

109. Безухов Н.И. Примеры и задачи по теории упругости, пластичности и пол зучести. -М.: Высшая школа, 1965. 318 с.

110. Алюшин Ю.А., Базайкин В.И, Милинис О.В. и др. Определение остаточных напряжений в тонких листах малоуглеродистой стали.// Заводская лаборатория. 1977. №4. с.485 -487.

111. Ильюшин А.А., Ломакин В.А., Шмаков А.П. Задачи и упражнения по механике сплошной среды. М.: Изд-во МГУ, 1978. - 162 с.

112. Contensou P.G. Note sur la cinematique du mobile dirige.// Bulletin ATMA. 1949. v.45.p.341 -347.

113. Andriex S. The Thin Shell Approuch for 3D Engineering Inverse Problems. // Inverse Problems in Engineering Mechanics (Eds. M. Tanaka and H.D.Bui). --Berlin : Springer Yerlag, 1993. p.539 - 551.

114. Kupradze V.D. Dynamical Problems in Elasticity.//Progress in Solid Mech. (Eds. I.N.Sneddon and R.Hill). Amsterdam : North Holland, 1963. p.3 - 17.

115. Gao Z., Mura T. On the Inversion of Residual Stresses from Surface Displacement.//J. Appl. Mech. 1969. №56. p.508 521.

116. Аэро Э.Л., Бессонов H.M., Булыгин A.H. Динамика моментной анизотропной жидкости.// ГТММ. 1996. т.60 , вып. 5. с. 778 785.

117. Булатов Э.А. К вопросу о моментной теории упругости. Плоская деформация. Сообщение 1.//Проблемы прочности. 1998. №3. с. 103 113.

118. Булатов Э.А. К вопросу о моментной теории упругости. Плоская деформация. Сообщение 3.// Проблемы прочности. 1998. №5. с.'82 87.

119. Коновалов A.B. Некоторые особенности моделирования напряжённо-деформированного состояния в процессах ОМД.// Кузнечно-штамповочное производство. 1998. №12. с.З 8.

120. Коновалов A.B. Определяющие соотношения упруго-пластической среды при больших пластических деформациях.//МТТ. 1997. №5. с.139 147.

121. Колмогоров B.JL, Федотов В.П., Горшков A.B. Трёхмерный анализ напряжённо-деформированного состояния в процессах обработки металлов давлением.//Кузнечно-штамповочное производство. 1998. №8. с.23-28.

122. Ерастов В.В. Алгоритм расчёта ускорений и инерционных сил при трёхмерной деформации.// Изв. вузов.Чёрная металлургия. 1997. №6. с.59 -61.

123. Ерастов В.В. Использование принципа наименьшего принуждения Гаусса для расчёта ускорений в полях скоростей из жёстких областей.// Изв. вузов.Чёрная металлургия. 1996. №8. с.З 1 33.

124. Алюшин Ю.А., Ерастов В.В. Ускорение при описании процессов пласти ческой деформации кинематически возможными полями из жёстких блоков.// Изв. вузов.Чёрная металлургия. 1990. №6. с.49 52.

125. Ерастов В.В., Пятайкин Е.М., Сергеев Н.М. Общий алгоритм расчёта про. цессов трёхмерной и плоской деформации.// Изв. вузов.Чёрная металлургия. 1997. №2. с.29-31.

126. Ерастов В.В., Пятайкин Е.М., Мезенцев В.В. Алгоритм расчётов процессов осесимметричной деформации с учётом инерционных сил на основе полей скоростей из непрерывно деформируемых областей.// Изв. вузов.Чёрная металлургия. 1997. №4. с.42 -44.

127. Базайкин В.В., Громов В.Е. Роль вязкости в силовой оценке процессов обработки металлов давлением.// Материалы Международной научно-практической конференции "Современные проблемы и пути развития металлургии". Новокузнецк. 1998. с. 125- 126.

128. Герасимов В.В., Логазяк H.A., Вайсбурд Р.Л., Колмогоров В.Л. Пакет программ МЕЛИСА для решения технологических задач ОМД методом линий скольжения.//Изв. вузов.Чёрная металлургия. 1997. №9. с.32 36.

129. Тарновский В.И. Исследование волочения как жёстко-пластического течения в сходящемся канале.// Изв. вузов.Чёрная металлургия. 1998. №5. с.31 -35.

130. Должанский A.M. Теоретическое определение толщины сухой изотермической смазки при волочении.//Изв. вузов.Чёрная металлургия. 1997. №1 с.47 50.

131. Должанский A.M. Исследование "тоннельного эффекта" в сухой смазке при волочении проволоки.//Изв. вузов.Чёрная металлургия. 1997. №3.с.31 -34.

132. Тарновский В.И., Безручко М.А. Объёмное деформированное состояние высоких заготовок при осадке.// Изв. вузов.Чёрная металлургия. 1998. №1. с.37 43.

133. Бровман М.Я. Определение усилий при осесимметричной деформации. // Изв. вузов.Чёрная металлургия. 1992. №2. с.31 32.

134. Бровман М.Я. Определение усилий при вдавливании цилиндрического пуансона в заготовку.//Изв. вузов Чёрная металлургия. 1997.№1. С.44 47.

135. Бровман М.Я. Динамические задачи осесимметричной пластической деформации.// Изв. вузов.Чёрная металлургия. 1996. №5. с.27 30.

136. Логинов Ю.И. Вариационное решение задачи формоизменения пористого цилиндра при осадке с прилипанием.// Изв. вузов.Чёрная металлургия. 1997. №11. с.38 41.

137. Смирнов О.М., Ершов А.Н., Чумаченко С.Е., Кропотов В.А. Анализ напряжённо-деформированного состояния заготовки в процессах осесимметричной штамповки осадкой с кручением.// Кузнечно-штамповочное производство. 1998. №6. с.9 12.

138. Каплунов Б.Г. Математическая модель расчёта пластической деформации в процессах штамповки заготовок тел вращения.// Изв. вузов.Чёрная металлургия. 1994. №2. с.20 25.23 7

139. Железков О.С. Исследование энергосиловых параметров процесса высадки прокаткой головок стержневых изделий.// 1994. №6. с.33 36.

140. Ерастов В.В. Барыльников. Расчёт силовых параметров и выбор рацио нальных способов высадки головки болта на основе метода верхней оценки.// Изв. вузов.Чёрная металлургия. 1996. №10. с.35 38.

141. Хван Д.В., Мозгалин B.JI. Штамп для высадки головок стержневых изделий.// Кузнечно-штамповочное производство. 1998. №3. с.33 34.

142. Илюкович Б.М., Измайлова М.К. Исследование процесса прокатки прямоугольной полосы на гладкой бочке вариационным методом.// Изв. вузов.Чёрная металлургия. 1996. №11. с.33 35.

143. Николаев В.А. Определение силы при холодной прокатке по экспериментальной дуге контакта.//Изв. вузов.Чёрная металлургия. 1997. №5.с.32 34.

144. Миленин А.А. О реализации граничных условий в напряжениях при моделировании процесса прокатки методом граничных элементов.//Изв. вузов.Чёрная металлургия. 1997. №4. с.28 -31.

145. Malinovski Z., Pitrzyk М., Lenard J.G.// J. Of Materials Processing Technology. 1993. V. 33. P.373 387.

146. Челышев H.А. Контактные напряжения и коэффициент плеча момента от усилия прокатки.// Изв.вузов.Чёрная металлургия. 1997. №8. с.38 39.

147. Максименко О.П., Подберезный Н.П. Уточнение модели контактно-гидродинамической смазки при прокатке.// Изв. вузов.Чёрная металлургия. 1997. №10. с.53 56.

148. Фёдоров НН., Перетятько В.Н., Фёдоров Н.А. Интенсивность пластической деформации при периодической асимметричной прокатке в эксцентричных валках.//Изв. вузов.Чёрная металлургия. 1997. №10.с.45 50.

149. Выдрин А.В. Инженерная методика оценки желобчатости полос при прокатке-волочении.// Изв. вузов.Чёрная металлургия. 1997. №5. с.40- 43.

150. Стеблянко В.Л., Литичевская Л.А., Литичевский В.М. Моделирование38формоизменения биметаллического профиля при прокатке в четырёх-валковом стрельчатом калибре.//Изв. вузов.Чёрная металлургия. 1997. №3. с.34 37.

151. Базайкин В.И. Полуобратный метод моделирования процесса холодной продольной прокатки. Кинематика.//Изв. вузов.Чёрная металлургия. 1998. №2. с.24 27.

152. Базайкин В.И., Громов В.Е., Лоскутов Д.Р. Полуобратный метод моделирования процесса холодной продольной прокатки. Напряжения.//Изв. вузов. Чёрная металлургия. 1998. №4. с.20 23.

153. Колмогоров С.Л., Хрущёв Р.И. Остаточные напряжения в изделиях, полученных осесимметричным пластическим деформированием.//Изв. вузов. Чёрная металлургия. 1993. №8. с.18-20.

154. Колмогоров Г.Л., Куратова H.A., Каменев С.А. Остаточные напряжения и предельная деформируемость при волочении осесимметричных изделий. //Изв. вузов.Чёрная металлургия. 1996. №5. с.23 -25.

155. Тепловые процессы при обработке металлов и сплавов давлением /Н.И.Яловой, М.А.Тылкин, П.И.Полухин, Д.И.Васильев /-М. : Высшая школа , 1973. 631 с.

156. Гарбер Э.А., Кузнецов С.А., Виноградов А.И., Семёнов С.Ю. Новая технология подготовки стальной катанки к волочению.//Изв. вузов.Чёрная металлургия. 1999. №3. с.46 48.239