МЕХАНИКА ТВЕРДОГО ТЕЛА <3 2008

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ РЕЗАНИЯ УПРУГОВЯЗКОПЛАСТИЧЕСКИХ МАТЕРИАЛОВ В ТРЕХМЕРНОЙ ПОСТАНОВКЕ

В работе методом конечных элементов проведено трехмерное моделирование неустановившегося процесса резания упруговязкопластической пластины (заготовки) абсолютно жестким резцом, движущимся с постоянной скоростью V0 при различных наклонах грани резца а (фиг. 1). Моделирование проводилось на основе связанной термомеханической модели упруговяз-копластического материала. Приведено сравнение адиабатического процесса резания и режима с учетом теплопроводности материала заготовки. Проведено параметрическое исследование процесса резания при изменении геометрии заготовки и режущего инструмента, скорости и глубины резания, а также свойств обрабатываемого материала. Варьировался размер толщины заготовки в направлении оси г. Напряженное состояние изменялось от плосконапряженного H = H/L < 1 (тонкая пластина) до плоскодеформируе-мого H > 1 (широкая пластина), где H - толщина, L - длина заготовки. Задача решалась на подвижной адаптивной лагранжево-эйлеровой сетке методом конечных элементов с расщеплением и использованием явно-неявных схем интегрирования уравнений . Показано, что численное моделирование задачи в трехмерной постановке позволяет исследовать процессы резания с образованием непрерывной стружки, а также с разрушением стружки на отдельные куски. Механизм этого явления в случае ортогонального резания (а = 0) может быть объяснен термическим разупрочнением с образованием адиабатических полос сдвига без привлечения моделей поврежденно-сти. При резании более острым резцом (угол а велик) необходимо привлечение связанной модели термического и структурного разупрочнения. Получены зависимости силы, действующей на резец при разных геометрических и физических параметрах задачи. Показано, что возможны квазимонотонный и осциллирующий режимы и дано их физическое объяснение.

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

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

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

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

Кроме того, как показывают экспериментальные исследования последних лет , при высоких скоростях деформирования е > 105-106 с-1 многие материалы обнаруживают аномальную температурную зависимость, связанную с перестройкой механизма движения дислокаций. Термофлуктуационный механизм сменяется на механизм фононного сопротивления, в результате которого зависимость сопротивления материала от температуры становится прямо противоположной: с увеличением температуры возрастает упрочнение материала. Такие эффекты могут привести к большим неприятностям при высокоскоростном резании. Эти проблемы в литературе до настоящего времени совершенно не изучены. Моделирование высокоскоростного процесса требует развитие моделей, учитывающих сложные зависимости вязкопластического поведения материалов и в первую очередь учета поврежденности и разрушения с образованием трещин и фрагментацией частиц и кусков деформируемого материала. Чтобы учесть все перечис-

8 Механика твердого тела, № 3

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

2. Постановка задачи. 2.1. Геометрия. Принимается трехмерная постановка задачи. На фиг. 1 изображены область и граничные условия в плоскости резания. В направлении перпендикулярном плоскости заготовка имеет конечную толщину И = Н/Ь (Ь - длина заготовки), которая варьировалась в широком диапазоне. Пространственная постановка допускает свободу движения обрабатываемого материала из плоскости резания и более плавный выход стружки, что обеспечивает более благоприятные условия резания.

2.2 Основные уравнения. Полная связанная система уравнений термоупруговязкопла-стичности состоит из уравнения сохранения импульса

рйи/йг = ; (2.1)

закона Гука с температурными напряжениями

йО;/йг = к1 - еы - «М) (2.2) уравнения притока тепла йй

рСе й- = К 0,.. - (3 X + 2ц)а0° е „■ + ко; р (2.3)

где Се - теплоемкость, К - коэффициент теплопроводности, к - коэффициент Куини-Тейлора, учитывающий разогрев материала за счет пластической диссипации.

Имеем также ассоциированный закон пластического течения

ер = Хй^/йо; (2.4)

и условия пластичности

Л, ЕЫ, X;, 9) = Оу (]ЕЫ, X;, 0) < 0 (2.5)

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

йX /йг = яЛк, Хк, 9) (2.6)

2.3 Модель материала. В работе принимается термоупруговязкопластическая модель типа Мизеса - модель пластичности с пределом текучести в виде мультипликативной зависимости (2.7), включающая деформационное и вязкопластическое упрочнение и термическое разупрочнение :

оу (ер, ¿* ,9) = [ а + в (ер)"]

где оу - предел текучести, ер1 - интенсивность пластических деформаций, 0 - относительная температура, отнесенная к температуре плавления 0т: " 0<0*

(0 - 0*) / (0т - 0*), 0* <0<0т

Материал детали принят однородным. В расчетах использовался относительно мягкий материал А12024-Т3 (упругие постоянные: Е = 73 ГПа, V = 0.33; пластические: А = 369 МПа, В = 684 МПа, п = 0.73, е0 = 5.77 ■ 10-4, С = 0.0083, т = 1.7, 9* = 300 К, 9т = 775 К, в = 0.9) и более жесткий 42СгМо4 (Е = 202 ГПа, V = 0.3, А = 612 МПа, В = 436 МПа, п = 0.15, е0 = 5.77 ■ 10-4, С = 0.008, т = 1.46, 9* = 300 К, 9т = 600 К, в = 0.9). Производится сравнение адиабатического процесса резки с решением полной термомеханической задачи.

2.4. Разрушение. Модель разрушения материала основывается на континуальном подходе Майнчена-Сака , основанном на моделировании зон разрушения дискретными частицами. В качестве критерия разрушения принимается критическая величина

интенсивности пластических деформаций ер:

ер = [йх + й2ехр (й311/12)][ 1 + й41п (ёр/ё0)](1 + й59) (2.8)

где й. - константы материала, определяемые из эксперимента.

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

е р = 1.0, скорость движения резца принималась равной 2 м/с и 20 м/с.

2.5. Метод интегрирования уравнений. Для интегрирования приведенной связанной системы уравнений термопластичности (2.1)-(2.8) целесообразно применить метод расщепления, развитый в работе . Схема расщеплений упругопластических уравнений заключается в расщеплении полного процесса на предиктор - термоупругий процесс, в

котором ер = 0 и все связанные с пластической деформацией операторы обращаются в нуль, и корректор - при котором полная скорость деформаций е = 0. На стадии предиктора система (2.1)-(2.6) относительно переменных, обозначенных тильдой примет вид

рйЪ/йг = а]

й аЛ = « - а§«9) рСей9/йг = К.9ц - (3X + 2ц)а90еи

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

АСТАШЕВ В.К., РАЗИНКИН А.В. - 2008 г.

Введение

Глава 1. Общая постановка задачи упруго-пластического деформирования 25

1.1. Кинематика процессов 25

1.2. Определяющие соотношения процессов упругопластического конечного деформирования 32

1.3. Постановка задачи конечного упругопластического деформирования 38

1.4. Постановка процесса разделения 42

Глава 2. Численное моделирование процессов конечного формоизменения 44

2.1. Численная формулировка проблемы 44

2.2. Метод интегрирования разрешающих соотношений 50

2.3. Алгоритмы решения краевых задач упруго-пластичности 51

2.4. Проверка правильности реализации математической модели 54

2.5. Анализ поведения модели при небольших деформациях 57

2.6. Моделирование процесса конечно-элементного разделения материала 58

2.7. Построение модели внедрения жесткого клина в полубесконеч-ное упруго-пластическое тело 60

2.8. Механизм учета трения в модели резания 62

Глава 3. Математическое моделирование процесса резания . 65

3.1. Процесс свободного резания 65

3.2. Факторы, влияющие на процесс стружкообразования 68

3.3. Граничные условия при моделировании 70

3.4. Конечно-элементная реализация процесса резания 74

3.5. Моделирование установившегося режима резания 75

3.6. Итерационный процесе на шаге 77

3.7. Обоснование выбора шага расчета и числа конечных элементов 80

3.8. Сравнение экспериментально найденных и расчетных значений сил резания 83

Список литературы

Введение к работе

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

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

кромки I .

Рис. 1. (а) Схема прямоугольного резания (б) Схема косоугольного резания.

Природа стружкообразования для рассмотренных случаев приблизительно одинакова. Различные авторы разделяют процесс образования стружки как на 4 , так и на 3 вида . Согласно различают три основных вида стружкообразования, показанные на рис. 2: а) прерывистое, включающее периодическое отделение элементов стружки в виде небольших сегментов; б) непрерывное стружкообразование; в) непрерывное с образованием нароста на инструменте.

Введение

По другой концепции еще в 1870 г. И. А. Тиме была предложена классификация типов стружек, образующихся при резании различных материалов. Согласно классификации И. А. Тиме, при резании конструкционных материалов в любых условиях образуются четыре вида стружек: элементная, суставчатая, сливная и надлома. Элементную, суставчатую и сливную стружку называют стружками сдвига, так как их образование связано с напряжениями сдвига. Стружку надлома иногда называют стружкой отрыва, так как ее образование связано с растягивающими напряжениями. Внешний вид всех перечисленных типов стружки изображен на рис. 3.

Рис. 3. Виды стружки по классификации Тиме.

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

Введение 8

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

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

На рисунке Зв - образование сливной стружки. Основным признаком является ее сплошность (непрерывность). Если на пути движения сливной стружки нет никаких препятствий, то она сходит непрерывной лентой, завиваясь в плоскую или винтовую спираль, пока часть стружки не отламывается под действием собственного веса. Поверхность стружки 1 - прилегающую к передней поверхности инструмента, называют контактной поверхностью. Она сравнительно гладкая, а при высоких скоростях резания отполирована в результате трения о переднюю поверхность инструмента. Ее противоположную поверхность 2 называют свободной поверхностью (стороной) стружки. Она покрыта мелкими зазубринками-насечкой и при высоких скоростях резания имеет бархатистый вид. Стружка соприкасается с передней поверхностью инструмента в пределах площадки контакта, ширина которой обозначена через С, а длина равна рабочей длине главного лезвия. В зависимости от рода и свойств обрабатываемого материала и скорости резания ширина площадки контакта в 1,5 - 6 раз больше толщины срезаемого слоя.

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

Введение 9

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

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

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

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

Введение 10

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

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

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

Введение

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

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

Стандартизированы как многие конструкции резцов в целом , так и их отдельные конструктивные элементы. Для унификации конструкций и присоединительных размеров резцедержателей принят следующий ряд сечений стержня, мм: квадратные со стороной а = 4, 6, 8, 10, 12, 16, 20, 25, 32, 40 мм; прямоугольные 16x10; 20x12; 20x16; 25x16; 25x20; 32x20; 21x25; 40x25;40x32;50x32; 50x40; 63x50 (отношение сторон Н:В=1,6 используется для получистовой и чистовой обработки, а Н:В=1,25 - для черновой обработки).

Общероссийский классификатор продукции предусматривает 8 подгрупп резцов с 39 видами в них. На конструкции резцов издано около 60 стандартов и технических условий. Кроме того, стандартизировано 150 типоразмеров пластин из быстрорежущей стали для всех типов резцов, около 500 типоразмеров твердосплавных напаиваемых пластин, 32 вида многогранных неперетачиваемых пластин (свыше 130 типоразмеров). В простейших случаях резец моделируется как абсолютно жесткий клин, без учета многих геометрических параметров .

Основные геометрические параметры резца с учетом сказанного выше.

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

Введение 12

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

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

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

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

Введение 13

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

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

Рис. 6. Форма передней поверхности резцов: а - плоская с фаской; б - криволинейная с фаской

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

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

углом Уф-^~5 - для резцов из быстрорежущей стали и У ф =-5..-25 . для резцов из твердых сплавов, всех видов керамики и синтетических сверхтвердых материалов.

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

Введение

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

Резцы из быстрорежущей стали и твердых сплавов с плоской передней поверхностью без фаски с ^ = 8..15 применяются для обработки хрупких материалов, дающих стружку надлома (чугун, бронза). При малой толщине среза, сравнимой с радиусом округления режущей кромки, величина переднего угла практически не влияет на процесс резания, так как деформирование срезаемого слоя и превращение его в стружку производится закругленной кромкой радиуса. В этом случае передние углы для всех типов инструментальных материалов принимаются в пределах 0...5 0 . Величина переднего угла существенно влияет на стойкость резцов.

Назначение главного угла в плане - изменять соотношение между шириной Ь и толщиной а среза при постоянных глубине резания t и подаче S.

Уменьшение угла повышает прочность вершины резца, улучшает тепло-отвод, повышает стойкость, но увеличивает силы резания P z и, Р у увеличивает

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

Конструкции резцов, в особенности с механическим креплением твердосплавных пластин, предусматривают ряд значений угла#>: 90, 75, 63, 60, 50, 45, 35, 30, 20, 10, что позволяет подобрать угол , наиболее соответствующий конкретным условиям.

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

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

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

T^-\ : " г

Введение

Рис. 7. Контур зоны стружкообразования, установленный путем исследования свободного резания с помощью киносъемки.

Рис. 8. Микрофотография зоны стружкообразования при резании стали с малой скоростью. На микрофотографии очерчены начальная и конечная границы зоны стружкообразования. (100 кратное увеличение)

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

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

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

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

Введение 17

довательно превращается в стружку. Метод пригоден при резании с очень малыми скоростями, не превышающими 0,2 - 0,3 м/мин, и дает только качественное представление о процессе стружкообразования.

Метод скоростной киносъемки. Хорошие результаты дает при съемке с частотой порядка 10 000 кадров в секунду и позволяет выяснить особенности процесса стружкообразования при практически используемых скоростях резания.

Метод делительной сетки. Основан на нанесении точной квадратной делительной сетки с размерами ячейки 0,05 - 0,15 мм. Делительная сетка наносится различными способами: накатыванием типографской краской, травлением, напылением в вакууме, трафаретной печатью, царапанием и т. п. Наиболее точным и простым способом является царапание алмазным индентором на приборе ПМТЗ для измерения микротвердости или на универсальном микроскопе. Для получения неискаженной зоны деформации, соответствующей определенной стадии стружкообразования, применяют специальные приспособления для «мгновенного» прекращения процесса резания, в которых вывод резца из-под стружки осуществляется сильной пружиной или энергией взрыва порохового заряда. На получившемся корне стружки с помощью инструментального микроскопа измеряют размеры ячеек искаженной в результате деформирования делительной сетки. Используя аппарат математической теории пластичности, по размерам искаженной делительной сетки можно определить вид деформированного состояния, размеры и форму зоны деформации, интенсивность деформации в различных точках зоны деформации и другие параметры, количественно характеризующие процесс стружкообразования.

Металлографический метод. Полученный с помощью приспособления для «мгновенного» прекращения резания корень стружки вырезают, тщательно полируют его боковую сторону, а затем протравливают соответствующим реактивом. Полученный микрошлиф корня стружки рассматривают под микроскопом при увеличении в 25-200 раз или делают микрофотографию. Изменение структуры

Введение

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

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

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

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

Введение 19

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

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

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

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

р, = c P f p sy K P

где Ср г - коэффициент, учитывающий влияние на силу некоторых постоянно действующих условий; *р- глубина резания; $^,- продольная подача; К р - обобщенный коэффициент резания; xyz - показатели степени.

Введение 20

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

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

напряжение q F = 0,285^ , где &к - действительное конечное сопротивление разрыву. На этом основании А.А.Розенбергом была получена другая формула для расчета главной составляющей силы резания:

(90-у) " cos/

-- їїдГ + Sin/

P z =0,28S K ab(2,05K a -0,55)

2250QK Qm5{9Q - Y) "

где Ъ - ширина срезаемого слоя.

Недостатком данной формулы является то, что для каждого конкретного

случая расчета сил требуется определение параметров К а и $к экспериментальным путем, который является весьма трудоемким. По данным многочисленных экспериментов было выявлено, что при замене криволинейной линии сдвига прямой, угол У близок к 45 , и следовательно формула примет вид:

d cos У

P z = - "- r + sin^

tg arccos

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

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

r max = гір"х ~ ъ) достигает некоторого определенного значения (для каждого материала своего).

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

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

Введение 23

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

К критериям прочности материалов неодинаково сопротивляющихся растяжению и сжатию, относятся критерий Шлейхера, Баландина, Миролюбова, Ягна . К недостаткам можно отнести сложность применения и плохое подтверждение экспериментальной проверкой.

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

Анализ перечисленных выше подходов к описанию теории разделения позволяет выделить следующие характерные особенности:

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

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

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

Экспериментальные методы способны однозначно предоставить инфор-

Введение

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

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

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

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

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

Основные положения и результаты работы доложены на Всероссийской научной конференции «Современные проблемы математики, механики и информатики» (г. Тула, 2002 г.), а также на зимней школе по механике сплошной среды (г. Пермь, 2003 г.), на международной научной конференции «Современные проблемы математики, механики и информатики» (г. Тула, 2003 г.), на научно-практической конференции «Молодые ученые центра России» (г. Тула, 2003 г.).

Определяющие соотношения процессов упругопластического конечного деформирования

Для индивидуализации точек среды выводится для начальной t - О фиксированной, так называемой, расчетной, конфигурации {KQ }, произвольная система координат 0, с помощью которой каждой частице ставится в соответствие тройка чисел (J,2,3) "приписанных" данной частице и неизменных в течение всего времени движения. Введенную в отсчетной конфигурации систему 0, вместе с базисом, =-г (/ = 1,2,3) называют неподвижной лагранжевой системой координат. Отметим, что в качестве материальных координат могут быть выбраны координаты частиц в начальный момент времени в системе отсчета. Необходимо отметить, что при рассмотрении процессов деформирования среды с зависимыми от истории деформирования свойствами независимо от используемых материальных или пространственных переменных используются две системы координат - одна из лагранжевых и эйлерова.

Как известно , возникновение в теле напряжений порождается деформацией материальных волокон, т.е. изменение их длин и взаимного расположения, поэтому основная задача, решаемая в геометрически нелинейной теории деформаций, заключается в разделении движения среды на трансляционное и "чисто деформационное" и указании мер для их описания. Следует отметить, что такое представление не является однозначным и можно указать несколько подходов к описанию среды, в которых разделение движения на переносное "квазитвердое" и относительное "деформационное" проводится различными способами . В частности, в ряде работ под деформационным движением понимают движение окрестности материальной частицы по отношению к под л вижному лагранжеву базису ёк; в работах в качестве деформационного рассматривается движение по отношению к жесткому базису, трансляционное движение которого определяется тензором ротации, связывающим главные оси левой и правой мер искажения. В данной работе разделение движения окрестности материальной частицы М(рис. 1.1) на трансляционное и деформированное строится на основе естественного представления градиента скорости в виде симметричной и антисимметричной части . В этом случае деформационная скорость определяется как относительная скорость движения частицы относительно жесткого ортогонального триэдра вихревого базиса, вращение которого задается тензором вихря Q. Следует отметить, что в общем случае движения среды главные оси тензора W проходят через разные материальные волокна. Однако, как показано в , для процессов простого и квазипростого нагружения в реальном диапазоне деформаций, исследование деформационного движения в вихревом базисе представляется весьма удовлетворительным. Вместе с тем при построении соотношений, описывающих процесс конечного деформирования среды, выбор мер должен удовлетворять ряду естественных критериев : 1) мера деформации должна быть сопряжена с мерой напряжения через выражение элементарной работы. 2) вращение материального элемента как абсолютно твердого тела не должно приводить к изменению мер деформации и их производных по времени - свойство материальной объективности. 3) при дифференцировании мер должно сохраняться свойство симметрии и условие разделения процессов формоизменения и изменения объема. Последнее требование является весьма желательным.

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

Для определения кривизны и круток траектории используются инварианты

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

Построение модели внедрения жесткого клина в полубесконеч-ное упруго-пластическое тело

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

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

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

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

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

Далее, исходя из предполагаемого характера процесса деформирования, выбирается тип частной теории пластичности (п. 1.2). Обрабатываемые данные экспериментов по одноосному растяжению образцов исследуемого материала формируют конкретный вид определяющих соотношений, используя в соответствии с требованиями п. 1.2 какой-либо один из наиболее распространенных методов аппроксимации экспериментальной кривой. При решении задачи определенный тип теории пластичности принимается неизменным для всего исследуемого объема в течение всего процесса. Справедливость выбора оценивается в последствии по кривизне траектории деформации, вычисляемой в наиболее характерных точках тела. Такой подход использовался при исследовании моделей технологических процессов конечного деформирования трубчатых образцов в режимах простого или близкого к нему внешнего нагружения. В соответствии с выбранной процедурой пошагового интегрирования весь интервал нагружения по параметру t разбивается на ряд достаточно малых этапов (шагов). В дальнейшем решение задачи для типичного шага строится по следующему алгоритму. 1. Для вновь определенной по результатам предыдущего шага конфигурации области вычисляются метрические характеристики деформированного про Глава 2. Численное моделирование процесса конечного формоизменения 53 странства. При первом шаге конфигурация области совпадает с конфигурацией, определенной при t = О. 2. Определяются упругопластические характеристики материала для каждого элемента в соответствии с напряженно-деформированным состоянием, отвечающим концу предыдущего шага. 3. Формируется локальная матрица жесткости и вектора сил элемента. 4. Задаются кинематические граничные условия на контактных поверхностях. При произвольной форме контактной поверхности используется известная процедура перехода к локальной системе координат . 5. Формируется глобальная матрица жесткости системы и соответствующий вектор сил. 6. Решается система алгебраических уравнений, определяется вектор столбец скоростей узловых перемещений. 7. Определяются характеристики мгновенного напряженно деформированного состояния, вычисляются тензоры скорости деформации W, вихря С1, скорости изменения объема 0, вычисляется кривизна траектории деформирования X 8. Интегрируются поля скоростей тензоров напряжения и деформаций, определяется новая конфигурация области. Определяется вид напряженно-деформированного состояния, зоны упругого и пластического деформирования. 9. Определяется достигнутый уровень внешних сил. 10. Проводится контроль выполнения условий равновесия, вычисляются векторы невязок. При реализации схемы без уточняющих итераций осуществля ется сразу переход к п. 1.

Факторы, влияющие на процесс стружкообразования

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

Согласно , непосредственное влияние на процесс стружкообразования при прямоугольном резании оказывают только четыре фактора: угол действия, передний угол инструмента, скорость резания и свойства материала. Все остальные факторы влияют косвенно. Для выявления данных зависимостей выбран процесс свободного прямоугольного резания материала по плоской поверхности, Заготовка разделена на две части линией предполагаемого разделения GA, верхний слой - будущая стружка, толщина снимаемого слоя о, оставшаяся заготовка толщиной h. Точка М - максимальная точка достижения острия резца при врезании, путь пройденный резцом - S. Ширина образца конечна и равна b . Рассмотрим модель процесса резания (Рис.3.1.) Считая, что в начальный момент времени образец является недеформированным, целым, без разрезов. Заготовка из двух поверхностей соединенных очень тонким слоем AG, толщиной 8 .а, где а толщина снимаемой стружки. AG - предполагаемая линия раздела (рис.3.1.). При движении резца возникает контакт по двум поверхностям режущего инструмента. В начальный момент времени разрушения не происходит - внедрение резца без разрушения. В качестве основного материала применяется упруго-пластический изотропный материал. В расчетах рассматривались как пластичные (способность материала получать большие остаточные деформации, не разрушаясь), так и хрупкие (способность материала разрушаться без заметной пластической деформации) материалы. В основу был положен низкоскоростной режим резания, при котором согласно исключается возникновение застойных явлений на передней поверхности. Еще одной особенностью является низкое тепловыделение в процессе резания, что не сказывается на изменении физических характеристик материала и, следовательно, на процессе резания и значении сил резания. Таким образом, возникает возможность как численно, так и экспериментально исследовать процесс резания не осложненного дополнительными явлениями прирезцового слоя.

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

Для решения этих задач была выбрана модель модульного построения комплекса, включающего общую оболочку, как объединяющий элемент, способный управлять подключением различных модулей . Единственным, глубоко интегрированным модулем, стал блок визуализации результатов. Остальные модули разделены на две категории: задачи и математические модели. Допускается не единственность математической модели. В оригинальном проекте их три для двух различных типов элементов. Каждая задача также представляет собой модуль связанный с математической моделью тремя процедурами и с оболочкой одной процедурой вызова модуля, таким образом, интеграция нового модуля сводится к вписыванию четырех строк в проект и перекомпилированию. В качестве инструмента реализации был выбран язык высокого уровня Borland Delphi 6.0 , обладающий всем необходимым для решения поставленной задачи в ограниченные сроки. В каждой задаче возможно применение либо автоматически построенных сеток конечных элементов, либо использования специально подготовленных с помощью пакета AnSYS 5.5.3 , и сохраненных в текстовом формате. Все границы можно разделить на два типа: динамические (где от шага к шагу меняются узлы) и статические (постоянные на всем протяжении расчета). Наиболее сложные при моделировании - динамические границы, если проследить процесс разделения по узлам, то при достижении критерия разрушения в узле принадлежащем границе Ол связь между элементами, которым принадлежит данный узел разрывается путем дублирования узла - добавления нового номера для элементов лежащих ниже линии раздела. Один узел приписывается J- и, а другой 1 із (рис.3.10). Далее из 1 и узел переходит в Ц и делее в Ц. Узел приписанный А п сразу же или через несколько шагов попадает на поверхность резца и переходит в Ц, где может быть откреплен по двум причинам: достижения критерия открепления, либо по достижении точки В, если стружколом определен при решению данной задачи. Далее узел переходит в Г9 , если стоящий перед ним узел уже откреплен.

Сравнение экспериментально найденных и расчетных значений сил резания

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

В основе процесса лежит проверка всех критериев и определение наименьшего коэффициента уменьшения шага, после чего шаг пересчитывается и так до тех пор, пока не станет К 0,99. Часть критериев в ряде задач может быть не задействована, далее описаны все критерии (рис. ЗЛО): 1. Запрет проникновения материала в тело резца - достигается проверкой всех узлов из і \Л 9»! 12 на пересечение границы передней режущей поверхности. Считая перемещение линейным на шаге, находится точка контакта поверхности и узла и определяется коэффициент уменьшения величины шага. Производится пересчет шага. 2. Выявляются элементы перешедшие предел текучести на данном шаге, определяется понижающий коэффициент для шага, так чтобы только несколько элементов "перешли" предел. Производится пересчет шага. 3. Выявляются узлы из некоторой области принадлежащие линии раздела GA, превысившие значение критерия разрушения на данном шаге. Определяется понижающий коэффициент для шага, так чтобы только один узел превысил значение критерия разрушения. Производится пересчет шага. Глава 3. Математическое моделирование процесса резания 4. Запрет проникновения материала в тело резца через заднюю режущую поверхность для узлов из А 6, если данная граница не закреплена. 5. Для узлов 1 8 может быть задано условие открепления и переход вЦ в точке В, если выбрано условие использовать при расчете стружколом. 6. При превышении деформации хотя бы в одном элементе более чем на 25%, производится понижение величины шага до предела 25% деформации. Производится пересчет шага. 7. Определяется минимальный коэффициент уменьшения величины шага, и если он меньше 0.99, то происходит пересчет шага, иначе переход к следующим условиям. 8. Первый шаг считается без трения. После просчета находятся направления движения узлов, принадлежащих А 8 и Ц, добавляется трение и происходит пересчет шага, направление действия силы трения сохраняется в отдельной записи. Если шаг просчитан с трением, то проверяется, не поменялось ли направление движения узлов, на которые действует сила трения. Если изменилось, то данные узлы жестко закрепляются на передней режущей поверхности. Производится пересчет шага. 9. Если осуществляется переход к следующему шагу, а не пересчет, то производится закрепление узлов подошедших к передней режущей по верхности - ПереХОД уЗЛОВ ИЗ і 12 К А 8 10. Если осуществляется переход к следующему шагу, а не пересчет, то для узлов принадлежащих 1 8 осуществляется расчет сил резания и если они отрицательны, то узел проверяется на возможность открепления, т.е. открепление осуществляется лишь в том случае, если он самый верхний. 11. Если осуществляется переход к следующему шагу, а не пересчет, то происходит выявление узла принадлежащего AG, превысившего значение критерия разрушения на данном шаге на допустимое (малое) значение. Включение механизма разделения: вместо одного узла создается два, один принадлежащий - и, другой 1 із; перенумерация узлов тела по специальному алгоритму. Переход к следующему шагу.

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

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

Для исследования чистого процесса резания, без влияния резкого повышения температуры на поверхности взаимодействия и в образце, при котором образуется сливная стружка, без образования нароста на режущей поверхности, согласно необходима скорость резания порядка 0,33 мм/сек. Приняв данную скорость как максимальную, получаем, что для продвижения резца на 1 мм необходимо просчитать 30 шагов (при условии интервала времени 0,1 - который обеспечивает наулучшую устойчивость процесса). При расчете, с использованием тестовой модели, при внедрении резца на 1 мм с учетом использования ранее описанных критериев и без учета трения было получено вместо 30 шагов 190. Это обусловлено уменьшением величины шага продвижения. Однако, в связи с тем, что процесс итерационный, реально было посчитано 419 шагов. Вызвано такое расхождением слишком большим размером шага, что ведет к многократному уменьшению величины шага в силу итерационного характера критериев. Так. при первоначальном увеличении числа шагов до 100 вместо 30, получено расчетное число шагов - 344. Дальнейшее увеличение числа до 150 ведет к увеличению числа расчетных шагов до 390, а следовательно и увеличению времени расчета. Исходя из этого, можно предположить, что оптимальным числом шагов, при моделировании процесса снятия стружки является 100 шагов на 1 мм врезания, при неравномерном разбиении сетки с числом элементов 600-1200. При этом, реальное число шагов, без учета трения будет не менее 340 на 1 мм, а с учетом трения не менее 600 шагов.

«МЕХАНИКА УДК: 539.3 А.Н. Шипачев, С.А. Зелепугин ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ВЫСОКОСКОРОСТНОГО ОРТОГОНАЛЬНОГО...»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2009 Математика и механика № 2(6)

МЕХАНИКА

А.Н. Шипачев, С.А. Зелепугин

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ

ВЫСОКОСКОРОСТНОГО ОРТОГОНАЛЬНОГО РЕЗАНИЯ МЕТАЛЛОВ1

Численно исследованы процессы высокоскоростного ортогонального резания металлов методом конечных элементов в рамках упругопластической модели среды в диапазоне скоростей резания 1 – 200 м/с. В качестве критерия отделения стружки применялось предельное значение удельной энергии сдвиговых деформаций. Выявлена необходимость использования дополнительного критерия стружкообразования, в качестве которого предложена предельная величина удельного объема микроповреждений.

Ключевые слова: высокоскоростное резание, численное моделирование, метод конечных элементов.

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

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

В научной литературе теоретические и экспериментальные исследования процессов высокоскоростного резания материалов представлены крайне недостаточно. Имеются отдельные примеры экспериментально-теоретических исследований влияния температуры на прочностные характеристики материала в процессе высокоскоростного резания . В теоретическом плане проблема резания материалов получила наибольшее развитие в создании ряда аналитических моделей ортогонального резания . Однако сложность проблемы и необходимость более полного учета свойств материалов, тепловых и инерционных эффектов привели к Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проекты 07-08-00037, 08-08-12055), РФФИ и Администрации Томской области (проект 09-08-99059), Минобрнауки РФ в рамках АВЦП «Развитие научного потенциала высшей школы» (проект 2.1.1/5993).

110 А.Н. Шипачев, С.А. Зелепугин использованию численных методов, из которых применительно к рассматриваемой задаче наибольшее распространение получил метод конечных элементов.

– – –

ляется с помощью уравнения состояния типа Ми – Грюнайзена, в котором коэффициенты подбираются на основе констант a и b ударной адиабаты Гюгонио.

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

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

Текущее значение этой энергии вычисляется с помощью формулы:

D Esh = Sij ij (5) dt Критическая величина удельной энергии сдвиговых деформаций зависит от условий взаимодействия и задается функцией начальной скорости удара:

c Esh = ash + bsh 0, (6) c где ash, bsh – константы материала. Когда Esh Esh в расчетной ячейке, эта ячейка считается разрушенной и удаляется из дальнейшего расчета, а параметры соседних ячеек корректируются с учетом законов сохранения. Корректировка заключается в удалении массы разрушенного элемента из масс узлов, принадлежавших этому элементу. Если при этом масса какого-либо расчетного узла становится нулевой, то данный узел считается разрушенным и также удаляется из дальнейшего расчета.

Результаты расчетов Расчеты проводились для скоростей резания от 1 до 200 м/с. Размеры рабочей части инструмента: длина верхней грани 1,25 мм, боковой 3,5 мм, передний угол 6°, задний угол 6°. Обрабатываемая стальная пластина имела толщину 5 мм, длину 50 мм, глубина резания – 1 мм. Материал обрабатываемой заготовки – сталь Ст3, материал рабочей части инструмента – плотная модификация нитрида бора.

Были использованы следующие значения констант материала обрабатываемой заготовки : 0 = 7850 кг/м3, a = 4400 м/с, b = 1,55, G0 = 79 ГПа, 0 = 1,01 ГПа, V1 = 9,2·10–6 м3/кг, V2 = 5,7·10–7 м3/кг, Kf = 0,54 м·с/кг, Pk = –1,5 ГПа, ash = 7·104 Дж/кг, bsh = 1,6·103 м/с. Материал рабочей части инструмента характеризуется константами 0 = 3400 кг/м3, К1 = 410 ГПа, К2 = К3 = 0, 0 = 0, G0 = 330 ГПа, где К1, К2, К3 – константы уравнения состояния в форме Ми – Грюнайзена.

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

– – –

Рис. 1. Форма стружки, заготовки и рабочей части режущего инструмента в моменты времени 1,9 мс (a) и 3,8 мс (б) при движении резца со скоростью 10 м/с Численное моделирование процессов высокоскоростного ортогонального резания 113 вплоть до выполнения критерия отделения стружки. При скоростях резания 10 м/с и ниже в образце возникают области, где несвоевременно срабатывает критерий отделения стружки (рис. 1, a), что свидетельствует о необходимости применения или дополнительного критерия, или замены используемого критерия на новый.

Дополнительно на необходимость корректировки критерия стружкообразования указывает форма поверхности стружки.

На рис. 2 показаны поля температуры (в K) и удельной энергии сдвиговых деформаций (в кДж/кг) при скорости резания 25 м/с в момент времени 1,4 мс после начала резания. Расчеты показывают, что поле температур практически идентично полю удельной энергии сдвиговых деформаций, что свидетельствует о том, что а 1520

– – –

Рис. 3. Поля удельного объема микроповреждений (в см3/г) в момент времени 1,4 мс при движении резца со скоростью 25 м/с Численное моделирование процессов высокоскоростного ортогонального резания 115 Заключение Численно исследованы процессы высокоскоростного ортогонального резания металлов методом конечных элементов в рамках упругопластической модели среды в диапазоне скоростей резания 1 – 200 м/с.

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

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

ЛИТЕРАТУРА

1. Петрушин С.И. Оптимальное проектирование рабочей части режущих инструментов // Томск: Изд-во Том. политехнического университета, 2008. 195 с.

2. Sutter G., Ranc N. Temperature fields in a chip during high-speed orthogonal cutting – An experimental investigation // Int. J. Machine Tools & Manufacture. 2007. No. 47. P. 1507 – 1517.

3. Miguelez H., Zaera R., Rusinek A., Moufki A. and Molinari A. Numerical modelling of orthogonal cutting: Influence of cutting conditions and separation criterion // J. Phys. 2006. V. IV. No. 134.

4. Hortig C., Svendsen B. Simulation of chip formation during high-speed cutting // J. Materials Processing Technology. 2007. No. 186. P. 66 – 76.

5. Campbell C.E., Bendersky L.A., Boettinger W.J., Ivester R. Microstructural characterization of AlT651 chips and work pieces produced by high-speed machining // Materials Science and Engineering A. 2006. No. 430. P. 15 – 26.

6. Зелепугин С.А., Коняев А.А., Сидоров В.Н. и др. Экспериментально-теоретическое исследование соударения группы частиц с элементами защиты космических аппаратов // Космические исследования. 2008. Т. 46. № 6. С. 559 – 570.

7. Зелепугин С.А., Зелепугин А.С. Моделирование разрушения преград при высокоскоростном ударе группы тел // Химическая физика. 2008. Т. 27. № 3. С. 71 – 76.

8. Иванова О.В., Зелепугин С.А. Условие совместного деформирования компонентов смеси при ударно-волновом компактировании // Вестник ТГУ. Математика и механика. 2009. № 1(5).

9. Канель Г.И., Разоренов С.В., Уткин А.В., Фортов В.Е. Исследования механических свойств материалов при ударно-волновом нагружении // Известия РАН. МТТ. 1999. № 5. С. 173 – 188.

10. Зелепугин С.А., Шпаков С.С. Разрушение двуслойной преграды карбид бора – титановый сплав при высокоскоростном ударе // Изв. вузов. Физика. 2008. № 8/2. С. 166 – 173.

11. Горельский В.А., Зелепугин С.А. Применение метода конечных элементов для исследования ортогонального резания металлов инструментом из СТМ с учетом разрушения и температурных эффектов // Сверхтвердые материалы. 1995. № 5. С. 33 – 38.

СВЕДЕНИЯ ОБ АВТОРАХ:

ШИПАЧЕВ Александр Николаевич – аспирант физико-технического факультета Томского государственного университета. E-mail: [email protected] ЗЕЛЕПУГИН Сергей Алексеевич – доктор физико-математических наук, профессор кафедры механики деформируемого твердого тела физико-технического факультета Томского государственного университета, старший научный сотрудник Отдела структурной макрокинетики Томского научного центра СО РАН. E-mail: [email protected], [email protected] Статья принята в печать 19.05.2009 г.

Похожие работы:

«Серия брифингов АПП по правовым вопросам Национальные правозащитные учреждения в роли национальных превентивных механизмов: возможности и проблемы Декабрь 2013 г. Введение Факультативный протокол Конвенции ООН против пыток (ФПКПП) создает систему предупреждения пыток, в основе которой лежит посещение мест содержания под стражей международным органом Подкомитетом, и национальными организациями национальными превентивными механизмами. Государства праве наделять одно или несколько существующих или...»

«Ученый совет: итоги заседания 30 января На заседании Ученого совета СПбГУ 30 января состоялось вручение медали Санкт-Петербургский университет, свидетельств победителей конкурса 2011 года по государственной поддержке молодых российских ученых-кандидатов наук, присвоение званий Почетного профессора СПбГУ, присуждение премий СПбГУ За научные труды, присвоение ученых званий, выборы заведующих кафедрами и конкурс научно-педагогических работников. Проректор по научной работе Николай Скворцов сделал...»

«1. Общие положения В целях выявления и поддержки талантливых молодых исследователей, содействия профессиональному росту научной молодежи, поощрения творческой активности молодых ученых РАН, других учреждений, организаций России и студентов высших учебных заведений России в проведении научных исследований Российская академия наук ежегодно присуждает за лучшие научные работы 19 медалей с премиями в размере 50 000 рублей каждая молодым ученым РАН, других учреждений, организаций России и 19 медалей...»

«ПРАВА ЧЕЛОВЕКА КОМИТЕТ ПО ЛИКВИДАЦИИ РАСОВОЙ ДИСКРИМИНАЦИИ Изложение фактов № 12 Всемирная кампания за права человека Серия Права человека: изложение фактов издается Центром по правам человека Отделения Организации Объединенных Наций в Женеве. В ней отражаются некоторые проблемы прав человека, которые находятся под пристальным вниманием или представляют особый интерес. Публикация Права человека: изложение фактов предназначена для самых широких кругов общественности; ее цель содействовать...»

«Лекция 3 РЫНОК И ГОСУДАРСТВЕННОЕ РЕГУЛИРОВАНИЕ Государство является единственной в своем роде организацией, занимающейся упорядоченным насилием в крупном масштабе. Мюррей Ротбард7 Я всегда отстаивал сбалансированный взгляд на роль государства, признавая ограничения и провалы как рыночного механизма, так и государства, но всегда предполагая, что они действуют совместно, на условиях партнерства. Джозеф Стиглиц8 Основные вопросы: 3.1. Фиаско, или провалы, рынка и необходимость государственного...»

2016 www.сайт - «Бесплатная электронная библиотека - Научные публикации»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам , мы в течении 1-2 рабочих дней удалим его.

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Математика и механика

МЕХАНИКА

А.Н. Шипачев, С.А. Зелепугин

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ВЫСОКОСКОРОСТНОГО ОРТОГОНАЛЬНОГО РЕЗАНИЯ МЕТАЛЛОВ1

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

Ключевые слова: высокоскоростное резание, численное моделирование, метод конечных элементов.

1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проекты 07-08-00037, 08-08-12055), РФФИ и Администрации Томской области (проект 09-08-99059), Минобрнауки РФ в рамках АВЦП «Развитие научного потенциала высшей школы» (проект 2.1.1/5993).

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

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

В численных расчетах используется модель повреждаемой среды, характеризующаяся возможностью зарождения и развития в ней трещин. Общий объем среды Ж составляют ее неповрежденная часть, занимающая объем Жс и характеризующаяся плотностью рс, а также занимающие объем Ж/ трещины, в которых плотность полагается равной нулю. Средняя плотность среды связана с введенными параметрами соотношением р = рс (Жс /Ж). Степень поврежденности среды характеризуется удельным объемом трещин V/ = Ж//(Ж р).

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

где р - плотность, г - время, и - вектор скорости с компонентами и, стгу = - (P+Q)5jj + Бу - компоненты тензора напряжений, Е - удельная внутренняя энергия, - компоненты тензора скоростей деформаций, Р = Рс (р/рс) - среднее давление, Рс - давление в сплошной компоненте (неповрежденной части) вещества, 2 - искусственная вязкость, Бу - компоненты девиатора напряжений.

Моделирование «отрывных» разрушений проводится с помощью кинетической модели разрушения активного типа :

При создании модели предполагалось, что в материале имеются потенциальные очаги разрушения с эффективным удельным объемом V:, на которых образуются и растут трещины (или поры) при превышении растягивающим давлением Рс некоторой критической величины Р = Р)У\/(У\ + V/), которая уменьшается по мере роста образовавшихся микроповреждений. Константы VI, V2, Рк, К/ подбирались путем сопоставления результатов расчетов и экспериментов по регистрации скорости тыльной поверхности при нагружении образца плоскими импульсами сжатия. Один и тот же набор констант материала используется при расчете как роста, так и схлопывания трещин или пор в зависимости от знака Рс.

Давление в неповрежденном веществе считается функцией удельного объема и удельной внутренней энергии и во всем диапазоне условий нагружения опреде-

Постановка задачи

Шу(ри) = 0 ;

0 , если |Рс | < Р* или (Рс > Р* и У^ = 0),

^ = | - я§п (Рс) к7 (Рс | - Р*)(У2 + У7),

если Рс < -Р* или (Рс > Р* и У^ > 0).

ляется с помощью уравнения состояния типа Ми - Грюнайзена, в котором коэффициенты подбираются на основе констант a и b ударной адиабаты Гюгонио.

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

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

Критическая величина удельной энергии сдвиговых деформаций зависит от условий взаимодействия и задается функцией начальной скорости удара:

Esh = ash + bsh U0 , (6)

где ash, bsh - константы материала. Когда Esh > Esch в расчетной ячейке, эта ячейка считается разрушенной и удаляется из дальнейшего расчета, а параметры соседних ячеек корректируются с учетом законов сохранения. Корректировка заключается в удалении массы разрушенного элемента из масс узлов, принадлежавших этому элементу. Если при этом масса какого-либо расчетного узла стано-

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

Результаты расчетов

Расчеты проводились для скоростей резания от 1 до 200 м/с. Размеры рабочей части инструмента: длина верхней грани 1,25 мм, боковой 3,5 мм, передний угол 6°, задний угол 6°. Обрабатываемая стальная пластина имела толщину 5 мм, длину 50 мм, глубина резания - 1 мм. Материал обрабатываемой заготовки - сталь Ст3, материал рабочей части инструмента - плотная модификация нитрида бора. Были использованы следующие значения констант материала обрабатываемой заготовки : р0 = 7850 кг/м3, a = 4400 м/с, b = 1,55, G0 = 79 ГПа, о0 = 1,01 ГПа, V = 9,2-Ю"6 м3/кг, V2 = 5,7-10-7 м3/кг, К= 0,54 м-с/кг, Pk = -1,5 ГПа, ash = 7-104 Дж/кг, bsh = 1,6-10 м/с. Материал рабочей части инструмента характеризуется константами р0 = 3400 кг/м3, К1 = 410 ГПа, К2 = К3 = 0, у0 = 0, G0 = 330 ГПа, где К1, К2, К3 - константы уравнения состояния в форме Ми - Грюнайзена.

Рис. 1. Форма стружки, заготовки и рабочей части режущего инструмента в моменты времени 1,9 мс (а) и 3,8 мс (б) при движении резца со скоростью 10 м/с

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

На рис. 2 показаны поля температуры (в К) и удельной энергии сдвиговых деформаций (в кДж/кг) при скорости резания 25 м/с в момент времени 1,4 мс после начала резания. Расчеты показывают, что поле температур практически идентично полю удельной энергии сдвиговых деформаций, что свидетельствует о том, что

Рис. 2. Поля и изолинии температуры (а) и удельной энергии сдвиговых деформаций (б) в момент времени 1,4 мс при движении резца со скоростью 25 м/с

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

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

0,1201 0,1101 0,1001 0,0901 0,0801 0,0701 0,0601 0,0501 0,0401 0,0301 0,0201 0,0101

Рис. 3. Поля удельного объема микроповреждений (в см /г) в момент времени 1,4 мс при движении резца со скоростью 25 м/с

Заключение

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

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

ЛИТЕРАТУРА

2. Sutter G., Ranc N. Temperature fields in a chip during high-speed orthogonal cutting - An experimental investigation // Int. J. Machine Tools & Manufacture. 2007. No. 47. P. 1507 - 1517.

3. Miguelez H., Zaera R., Rusinek A., Moufki A. and Molinari A. Numerical modelling of orthogonal cutting: Influence of cutting conditions and separation criterion // J. Phys. 2006. V. IV. No. 134. P. 417 - 422.

4. Hortig C., Svendsen B. Simulation of chip formation during high-speed cutting // J. Materials Processing Technology. 2007. No. 186. P. 66 - 76.

5. Campbell C.E., Bendersky L.A., Boettinger W.J., Ivester R. Microstructural characterization of Al-7075-T651 chips and work pieces produced by high-speed machining // Materials Science and Engineering A. 2006. No. 430. P. 15 - 26.

7. Зелепугин С.А., Зелепугин А.С. Моделирование разрушения преград при высокоскоростном ударе группы тел // Химическая физика. 2008. Т. 27. № 3. С. 71 - 76.

9. Канель Г.И., Разоренов С.В., Уткин А.В., Фортов В.Е. Исследования механических свойств материалов при ударно-волновом нагружении // Известия РАН. МТТ. 1999. № 5. С. 173 - 188.

10. Зелепугин С.А., Шпаков С. С. Разрушение двуслойной преграды карбид бора - титановый сплав при высокоскоростном ударе // Изв. вузов. Физика. 2008. № 8/2. С. 166 - 173.

ШИПАЧЕВ Александр Николаевич - аспирант физико-технического факультета Томского государственного университета. E-mail: [email protected]

ЗЕЛЕПУГИН Сергей Алексеевич - доктор физико-математических наук, профессор кафедры механики деформируемого твердого тела физико-технического факультета Томского государственного университета, старший научный сотрудник Отдела структурной макрокинетики Томского научного центра СО РАН. E-mail: [email protected], [email protected]

V 0 z. H/L 1 (широкая пластина), где Н - толщина, L - длина заготовки. Задача решалась на подвижной адаптивной лагранжево-эйлеровой сетке методом конечных элементов с расщеплением и использованием явно-неявных схем интегрирования уравнений ...

В работе методом конечных элементов проведено трехмерное моделирование неустановившегося процесса резания упруговязкопластической пластины (заготовки) абсолютно жестким резцом, движущимся с постоянной скоростью V 0 при различных наклонах грани резца а (фиг. 1). Моделирование проводилось на основе связанной термомеханической модели упруговяз-копластического материала. Приведено сравнение адиабатического процесса резания и режима с учетом теплопроводности материала заготовки. Проведено параметрическое исследование процесса резания при изменении геометрии заготовки и режущего инструмента, скорости и глубины резания, а также свойств обрабатываемого материала. Варьировался размер толщины заготовки в направлении оси z. Напряженное состояние изменялось от плосконапряженного Я = H/L 1 (широкая пластина), где Н - толщина, L - длина заготовки. Задача решалась на подвижной адаптивной лагранжево-эйлеровой сетке методом конечных элементов с расщеплением и использованием явно-неявных схем интегрирования уравнений . Показано, что численное моделирование задачи в трехмерной постановке позволяет исследовать процессы резания с образованием непрерывной стружки, а также с разрушением стружки на отдельные куски. Механизм этого явления в случае ортогонального резания (а = 0) может быть объяснен термическим разупрочнением с образованием адиабатических полос сдвига без привлечения моделей поврежденно-сти. При резании более острым резцом (угол а велик) необходимо привлечение связанной модели термического и структурного разупрочнения. Получены зависимости силы, действующей на резец при разных геометрических и физических параметрах задачи. Показано, что возможны квазимонотонный и осциллирующий режимы и дано их физическое объяснение.

Какие птицы улетают на юг?

Жаворонок хохлатый (galerida cristata)

Муж на час – бизнес для хозяйственных мужчин

Бакалея это какие продукты