Напряженно-деформированное состояние континуума "жесткий индентор - упругопластическая среда" при динамических нагружениях тема автореферата и диссертации по механике, 01.02.04 ВАК РФ
|
Овчинникова, Наталья Владимировна
АВТОР
|
||||
|
кандидата технических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
|
Саратов
МЕСТО ЗАЩИТЫ
|
||||
|
2013
ГОД ЗАЩИТЫ
|
|
01.02.04
КОД ВАК РФ
|
||
|
|
||||
На правах рукописи
ОВЧИННИКОВА Наталья Владимировна
НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ КОНТИНУУМА «ЖЕСТКИЙ ИНДЕНТОР - УПРУГОПЛАСТИЧЕСКАЯ СРЕДА» ПРИ ДИНАМИЧЕСКИХ НАГРУЖЕНИЯХ
Специальность 01.02.04 - Механика деформируемого твердого тела
Автореферат диссертации на соискание ученой степени кандидата технических наук
005544986
Саратов 2013
005544986
Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Саратовский государственный технический университет имени Гагарина Ю.А.»
Научный руководитель - доктор технических наук, профессор,
Лауреат государственной премии СССР в области науки и техники Чеботаревский Юрий Викторович
Официальные оппоненты: Белосточный Григорий Николаевич
доктор технических наук, профессор, ФГБОУ ВПО «Саратовский государственный университет им. Н. Г. Чернышевского», профессор кафедры «Математическая теория упругости и биомеханика»
Тарлаковский Дмитрий Валентинович доктор физико-математических наук, профессор, ФГБОУ ВПО «Московский государственный университет им. М. В. Ломоносова», Научно-исследовательский институт механики, заведующий лабораторией динамических испытаний
Ведущая организация - Федеральное государственное учреждение
науки Институт проблем точной механики и управления Российской академии наук (г. Саратов)
Защита состоится «3» декабря 2013 г. в 15:00 часов на заседании диссертационного совета Д 212.242.06 при ФГБОУ ВПО «Саратовский государственный технический университет имени Гагарина Ю.А.» по адресу: 410054, Саратов, ул. Политехническая, 77, ауд. 319/1.
С диссертацией можно ознакомиться в научно-технической библиотеке ФГБОУ ВПО «Саратовский государственный технический университет имени Гагарина Ю.А.».
Автореферат разослан « 1 » ноября 2013 г.
/'
Ученый секретарь с /
диссертационного совета — Попов B.C.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. В последнее время для улучшения эксплуатационных характеристик различных механизмов получили применение комбинированные методы упрочнения рабочих поверхностей их деталей из металла, объединяющих одновременно несколько способов силового воздействия на обрабатываемый материал. В частности, один из них основан на использовании в процессе поверхностного упрочнения металлов импульсного нагружения рабочего инструмента в сочетании с ультразвуковым воздействием на него. Применение такого режима обработки позволяет существенно повысить качество обрабатываемой поверхности и обеспечить более высокую износоустойчивость готовых деталей. В её процессе проявляются новые динамические эффекты воздействия ультразвука на материал, прогнозирование которых без проведения предварительных аналитических и численных исследований весьма затруднительно. Для осознанного управления технологическим процессом поверхностного упрочнения необходимо определить комплекс геометрических, механических и физических параметров, его характеризующих и исследовать их влияние на результаты обработки. Механические процессы, происходящие при взаимодействии рабочего инструмента и поверхности детали, настолько сложны и разнообразны, что учесть все их специфические особенности в рамках решения единой задачи не представляется возможным. Косвенными характеристиками, позволяющими судить об уровне упрочнения, могут служить возникающие в его процессе остаточные пластические напряжения и деформации в обрабатываемом материале, а также его плотность, являющаяся функцией интенсивности остаточных деформаций. В связи с этим весьма актуальным представляется исследование напряженно-деформированного состояния упрочняемого материала при комбинированном динамическом силовом и ультразвуковом воздействиях.
С точки зрения механики сплошных сред задача об исследовании напряженно-деформированного состояния упрочняемого материала при таком воздействии является контактной и должна рассматриваться в динамической постановке с учетом упругопластического характера поведения обрабатываемого материала. Решению контактных задач посвящены пионерские работы Г. Герца, И. Я. Штаермана, Г. Лоренца, А. Н. Динника, Н. М. Беляева. Дальнейшее развитие это направление получило в работах Н. И. Мусхелишвили, Л. А. Галина, В. М. Александрова и многих других российских и зарубежных ученых. Подавляющее большинство опубликованных ими работ посвящено исследованию НДС рассматриваемых объектов при статическом нагружении одного из контактирующих тел. Однако класс задач механики деформируемого твердого тела, рассматриваемых в данной работе, до сих пор является малоизученным. Поэтому разработка новых подходов к решению контактных задач для упругопластических тел, находящихся под воздействием описанных выше сложных динамических нагрузок, представляется весьма актуальной.
С учетом этого целью настоящей диссертационной работы является разработка новых подходов к решению задачи о напряженно-деформированном
состоянии континуума «жесткий индентор - упругопластическая среда» при различных видах статического и динамического нагружения индентора и неупругом поведении материала среды.
Задачи исследования.
1. Постановка динамической задачи о напряженно-деформированном состоянии континуума «жесткий индентор - упругопластическая среда».
2. Получение вариационного уравнения движения континуума «индентор -деформируемая среда», а также применение метода штрафных функций и метода множителей Лагранжа для учета условий контакта непосредственно в вариационных соотношениях, описывающих движение континуума.
3. Исследование напряженно-деформированного состояния «жесткий индентор - упругопластическая среда» с помощью метода конечных элементов.
4. Разработка упрощенной физико-математической модели для решения задачи о напряженно-деформированном состоянии континуума «абсолютно жесткий индентор - упругопластическая среда».
5. Изучение динамики поведения абсолютно твердого индентора, контактирующего с упругопластической средой.
6. Исследование напряженно-деформированного состояния упругопластической среды под действием эквивалентной нагрузки.
Методы исследований. При выполнении работы использованы аппарат дифференциального и интегрального исчисления, вариационного исчисления, численные методы математического анализа, включая метод конечных элементов.
Для описания упругопластических свойств материала среды в рамках теории течения принята модель изотропно кинематического упрочнения.
Достоверность полученных результатов обеспечивается корректностью математической постановки рассматриваемых задач, использованием фундаментальных принципов механики, а также их совпадением с известными аналитическими решениями в некоторых частных случаях. Полученные численные результаты качественно подтверждаются имеющимися в периодической литературе экспериментальными данными и соответствуют физике исследуемых процессов.
На защиту выносятся следующие основные положения:
1. Вариационные уравнения движения континуума «жесткий индентор -деформируемая среда», полученные на основе вариационного принципа Даламбера - Лагранжа с применением методов штрафных функций и множителей Лагранжа и позволяющие учесть искривление деформируемых вследствие контакта поверхностей без наложения ограничений на их форму и физические уравнения материалов контактирующих тел (упругое, вязко-упругое, упруго-пластическое);
2. Упрощенная физико-математическая модель процесса контактного взаимодействия тел, образующих континуум «абсолютно твердый индентор -деформируемая среда», позволяющая исследовать напряженно-деформированное состояние среды с меньшими затратами вычислительных ресурсов:
- при динамическом нагружении индентора и идеально упругом поведении материала среды;
— при статическом нагружешш индентора и упруго - пластическом поведении материала среды.
3. Результаты численного исследования напряженно-деформированного состояния тел, образующих континуум, полученные для физической модели изотропно кинематического упрочнения материала среды при различных видах приложения нагрузки к индентору (статическая, динамическая, импульсная, с приложением ультразвуковых воздействий).
Научная новизна работы заключается в следующем:
— выполнена классическая постановка динамической упруго пластической задачи о напряженно-деформированном состоянии континуума «жесткий индентор - деформируемая среда» в виде, позволяющем учесть нарушение первоначальной геометрической формы деформируемых в процессе контакта поверхностей без наложения ограничений на их форму;
— на основе принципа Даламбера — Лагранжа с применением методов штрафных функций и множителей Лагранжа выведены вариационные уравнения движения континуума без наложения, каких-либо ограничений на форму деформированных поверхностей контактирующих тел, физические уравнения их материалов и геометрические соотношения;
— с целью применения метода конечных элементов для получения численных результатов проведена дискретизация рассматриваемой задачи в вариационной постановке;
— получены уравнения для определения координат точек, лежащих на границе контактной области в случае применения метода штрафных функций и при использовании метода множителей Лагранжа;
— на базе физической модели изотропно кинематического упрочнения материала среды проведено численное исследование напряженно-деформированного состояния тел, образующих континуум при различных видах приложения нагрузки к индентору (статическая, динамическая, импульсная, с приложением ультразвуковых воздействий);
— сформулирована упрощенная инженерная расчетная модель изучаемого процесса контактного взаимодействия тел, образующих континуум «абсолютно твердый индентор - деформируемая среда», позволяющая в некоторых частных случаях нагружения индентора исследовать напряженно-деформированное состояние среды при упругом и упруго-пластическом поведении её материала с меньшими затратами вычислительных ресурсов;
— определены границы применимости предложенной упрощенной расчетной модели континуума «абсолютно твердый индентор — деформируемая среда» для исследования напряженно-деформированного состояния среды при различных видах нагружения индентора.
Практическая ценность работы состоит в использовании её результатов при проектировании технологического оборудования в ООО «Научно-производственное предприятие нестандартных изделий машиностроения» и в учебном процессе при подготовке инженеров по специальности 151701.65 «Проектирование технологических машин и комплексов» в Саратовском государственном техническом университете имени Гагарина Ю.А.
Апробация работы. Основные положения и результаты работы докладывались на 5 научно-технических конференциях и 3 научных семинарах,
в том числе: XV и XVI Международных симпозиумах «Динамические и технологические проблемы механики конструкций и сплошных сред» им. А. Г. Горшкова (Ярополец, 16-20 февраля 2009 г, 15-19 февраля 2010г.); Всероссийской научно-технической конференции «Совершенствование техники, технологий и управления в машиностроении» (Саратов, 20-24 октября 2009 г.); Всероссийской научно-практической конференции молодых ученых «Инновации и актуальные проблемы техники и технологий», (Саратов, 26-29 октября 2010 г.); XXIV Международной научной конференции «Математические методы в технике и технологиях» ММТТ-24 (Саратов, 2011 г.), научном семинаре «Механика деформируемого твердого тела» СГТУ (Саратов, июнь 2013 г.).
Публикации. Основные результаты работы опубликованы в 12 печатных работах, в том числе 7 работ в изданиях, рекомендованных ВАК РФ для публикации материалов кандидатских и докторских диссертаций.
Структура и объем работы. Диссертационная работа состоит из введения, 4 глав, заключения, списка использованной литературы и приложения.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснованы актуальность темы диссертационной работы, ее научная новизна и практическая ценность. Сформулирована цель и поставлены задачи исследования.
В первой главе приводится обзор литературных источников, тем или иным образом связанных с изучаемой проблемой. На основе их анализа отмечается, что ни один из предлагаемых в изученных литературных источниках подходов без существенных изменений и дополнений не может быть положен в основу решения рассматриваемой в данной работе проблемы. Обосновывается выбор направления исследований. Формулируется математическая модель изучаемого технологического процесса в виде континуума, включающего цилиндрический упругий жесткий индентор со сферической рабочей поверхностью и упругопластическую среду, размеры и форма которой в силу локальности изучаемого процесса особого значения не имеют. Силовое воздействие рабочего инструмента на обрабатываемый материал моделируется путем приложения к нему направленного вдоль общей оси симметрии индентора и среды комбинированного усилия с постоянно действующей и переменной составляющими:
F(f) = F0+/;0(í) + ^v ятаг, (1)
где Р0, Рв О) и ^ - соответственно постоянная и переменная составляющие усилия, и амплитуда переменной составляющей, колеблющейся с ультразвуковой частотой а, г - время. Наличие постоянной составляющей усилия в (1) диктуется необходимостью обеспечения в реальном технологическом процессе гарантированного натяга между рабочим инструментом и обрабатываемой поверхностью и предотвращения тем самым его отрыва от неё в процессе обработки. Закон изменения во времени переменных составляющих обусловлен видом технологического воздействия и может иметь как импульсный, так и характер ультразвуковьк колебаний или их комбинацию.
При постановке задачи о напряженно-деформированном состоянии континуума для удобства формулирования граничных условий в среде выделяется конечная область в виде правильного кругового цилиндра, к центру одной из торцевых поверхностей которого прикладывается технологический инструмент (рис. 1). Её размеры - радиус К2 и толщина /?2 произвольны, ничем не ограничены и могут принимать любые значения. Нижняя торцевая поверхность выделенного объема
среды
жестко закреплена. Боковые
Г77Шпш^1П1ПП1
Рис. 1
поверхности индентора и среды Г5 и Г2, а также верхняя торцевая поверхность среды Г3 и рабочая поверхность индентора Г4, за исключением зоны контакта Гг (на рисунке не показана), свободны от внешней нагрузки. Зона контакта определяется как пересечение контактирующих поверхностей = Г, п Г4. С учетом того, что индентор по условиям своего закрепления может совершать только возвратно-поступательное движение вдоль общей оси симметрии с объемом Ог, исследование напряженного и деформированного состояния континуума проводится в осесимметричной постановке. В предположении, что упругие и пластические деформации материалов среды и индентора являются малыми, и процесс деформирования является изотермическим, определение НДС континуума в рамках теории течения сведено к решению следующей контактной задачи:
Ъо[г | <71гг-а'т | да'ГТ .¿У Э^ . 1 Э , .¿У
Э г
°гг = <*л =0 при (г,в,г)е Г2 и Г5. (3)
=-р(0> агг =0 при (г,е,г)е Г6, (4)
иг= 0 при (г,е,г)е Г,, (5)
ояя = аях = 0 при (г,0, г)е Г3 \ Гг и (г, 0, е Г4 \ Г,., (6)
о'* =<*™ и = о"г = 0 при (г, 0, г)е Гс, (7)
"я + и" - 8 < 0, (8)
, , ^ (1и[ с1и[ иг=и,= 0, —- =—- = 0 при Г = 0 (9)
Ж Ж
В соотношениях (1) - (9): а'„И,] = г,в,г) и н' = н,'(г,г,г) (/'= г,г) (/ = /,//)-компоненты тензора напряжений и вектора перемещений соответственно; />' (/ = /,//) - плотность материала среды или индентора; Г - время; />(г) -равномерно распределенная нагрузка по верхней торцевой поверхности индентора, эквивалентная приложенной к нему сосредоточенной силе:
^) = 2я- \pitydr-,
о
а'т и (1=1,11) — нормальные и касательные компоненты вектора
напряжений на площадках перпендикулярных внешним нормалям деформированных поверхностей индентора и среды в области контакта и её окрестностях, связанные с компонентами тензора напряжений соотношениями:
< (Ю)
<т1г=»1{<гУ.+<»,г)-»'г{<п'.+<п,г) (1 = 1,Л), (И)
где«' и п[ (/ = /, 7/) —направляющие косинусы главных внешних нормалей
п' и п" деформированных поверхностей индентора и среды, включая и область контакта. Здесь и в дальнейшем верхний индекс 1 = 1 или I - II означает, что данная величина независимо от её физической природы относится соответственно либо к индентору, либо к среде.
Компоненты полной деформации е9> состоящие из упругих е^' и пластических е? частей (составляющих)
, 02)
где -б^аЛ (Л7=^,г), (13)
и
связаны с компонентами вектора перемещений соотношениями Коши: ди, _ 8иг _иг _1 (ди диА
£rr > ezz --Г—. Еее - —. епг -- — + > (14)
or dz г 2\or dz J
фЦ, ) - пластический потенциал.
В силу непреодолимых математических трудностей получить аналитическое решение контактной задачи в приведенной выше постановке не представляется возможным.
Поэтому во второй главе диссертации для упрощения её решения предложена замена исходной задачи в классической постановке её вариационным аналогом, получаемым на основе дифференциального вариационного принципа Даламбера - Лагранжа для распределенной системы материальных точек, подчиненной неудерживающим связям. С этой целью наряду с радиальной иг и осевой иг составляющими вектора перемещения^ рассмотрение введены соответствующие им компоненты вектора скорости:
dur . du,
и их вариации Svr и Svz. После проведения ряда сложных и громоздких преобразований вариационное уравнение движения континуума «индентор-среда» получено в виде
*2 h «I X, (г.') Я|
<5Ф{r,z,t)= \ ]fnrdzdr+\ j/'rdzdr- \dv'lprdr=0 , (16)
где
дг г дг Л (1 = 1П) (1?)
дг дг т
Х: ('%') (1-1,П) - функции, входящие в уравнения деформированных вследствие контакта поверхностей индентора и среды
Х,{г,2,1) = г-Х,(г,1)=0 (1 = 1,11). (18)
Уравнение (16) описывает движение континуума «индентор -деформируемая среда» с учетом граничных условий на поверхностях тел его образующих и обобщенного контактного условия:
(19)
но без учета требования <т'т < 0 (1 = 1,11) в области контакта и условия отсутствия взаимопроникновения («условия непроникновения»):
Г„=<+У'п'< 0. (20)
В соотношениях (19) и (20)
(1 = 1,11)- (21)
составляющие векторов скорости принадлежащих индентору и среде точек возможного контакта в направлении соответствующей нормали.
Для введения условия «непроникновения» (20) и требования а'т < 0 (1 = 1,11) в области контакта непосредственно в вариационное соотношение (16), использованы два метода: штрафных функций и множителей Лагранжа. С применением первого из них вариационное уравнение движения континуума (16) преобразовано к виду:
Я2 Ы Л, х;(г,0
8П(г, 2,/)= | \frdzdr + | \frdzdr -
°Х„М °Х,(Л|.'Н (22)
- \bvlprdr + Г— + =
О О
где = {^(г)},
[1
(23)
г' и г" — координаты точек возможного контакта, лежащих на поверхности среды и рабочей поверхности индентора соответственно (Рис.2). Задание штрафной функции в виде (22) с учетом (23) позволило сделать её положительно определенной, что обеспечило выполнение требования ет" < 0 в области контакта, и отличной от нуля только в точках виртуального проникновения тела индентора в материал среды.
Система
уравнении для определения начальных координат ги и гя граничных точек контактной области получена исходя из условия равенства координат г1а,г" и г,', z" точек пересечения поверхностей индентора и среды при наличии виртуального проникновения (рис. 2):
(16)
>о2 "(гЯ-Яо + = и»{г"Л/) .
С применением метода множителей Лагранжа уравнение преобразовано к нестрогому вариационному равенству:
) " + } '}/"гсЫг- '}&>; ргс1г + З[я{у'„ + )}гс1г> 0, (24)
°-Г,(*..')-*. О л-„(г.,) О 0Иг
где Л - множитель Лагранжа, удовлетворяющий в каждой точке контактной поверхности соотношению Я, и равный нулю за её пределами.
Система уравнений для определения начальных координат г' и г"
граничных точек контактной области. Ввиду громоздкости полученных соотношений и ограниченности объема реферата эти уравнения здесь не приводятся.
Вариационные соотношения (22) и (24) универсальны в том смысле, что они получены без наложения каких либо ограничений на форму деформированных поверхностей индентора и среды вследствие контакта, геометрические соотношения и на физические уравнения, определяющие состояние материала тел, образующих приведено доказательство того, что определение точек континуума на основе вариационных соотношений (22) или (24), обеспечивающих в каждый момент времени безусловный минимум соответствующего функционала, соответствует решению исходной задачи (2) - (9).
В третьей главе проведена дискретизация задачи в вариационной постановке путем разбиения континуума на конечные элементы. Связь между компонентами вектора скорости любой точки континуума и их значениями в узловых точках разбиения задана соотношениями:
континуум, скоростей
Рис. 2
В работе движения
= (/ = г,г, / = /,//), (25)
где п, (1 = 1,11) - количество узлов в объемах индентора и среды соответственно, Ы'к(г,г) - функции формы. С учетом (25) на базе вариационного уравнения (22) получена следующая система из 2(п, + п„) обыкновенных дифференциальных уравнений для определения составляющих скоростей в узловых точках разбиения:
= "\flrdzdr-и)гс1г, и = иО
1, л у1 Я] Х;М "1
М.-Т- = -/ \firdzdr-\—М[п,\у{г,и)г<1г + \Nlprdr,
т=1 ^ О X/0«, . о
Ё = -] )/"гйгАг - /-1 Л^"п>(/-, г, О^г, (26)
1>ХнМ 0 "г
При использовании нестрогого вариационного равенства (24) дополнительно к (25) введена связь между значениями множителя Лагранжа в любой точке контактной поверхности и его значениями в узловых точках контакта следующим образом:
Я(лг,/) = |Х(г,г)Л(г). (27)
где ЛДг,г) - функции формы множителей Лагранжа; л, - общее число узлов, принадлежащих индентору и среде, находящихся в контакте. Смешанная система из 2(п, + п„) обыкновенных дифференциальных уравнений и пс нестрогих алгебраических неравенств для определения составляющих скоростей в узловых точках разбиения континуума и значений множителей Лагранжа в узловых точках индентора и среды, находящихся в контакте, получена после подстановки (25) и (27) в соотношение (24) и выполнения ряда преобразований в виде
«, Л у' Я] Х| (<■•') щ а , _
а1 ПХЛЧЫУ-Щ т=1„Пг
л, Л Я1 Х/(г,г) Л, "г " 1
'"=' О м=1 (I П,
"и Л ,>" Й2 "г " 1
«=1 ас о Х;/(г,г) т=1 о П2
(28)
и=,0"г т='°"г * = (29)
+ё > о
т=1 о '"=' О
При решении конкретных задач системы уравнений (26) и (28) и нестрогие неравенства (29) дополнены геометрическими соотношениями (14), устанавливающими зависимости между компонентами вектора перемещений и тензора деформаций, и физическими уравнениями (13), описывающими поведение материалов контактирующих объектов в процессах нагружения и разгрузки. Материал индентора задан идеально упругим в силу его относительно высокой жесткости. Для материала среды приняты условие текучести Губера — Мизеса и модель изотропно-кинематического упрочнения, в соответствии с которой пластический потенциал ф(с$-) определен соотношением
Ф = /(ст,а)-ст° =0, (30)
где
/(ст, а) = ^(агг - ат )2 + (ат - оа )2 + (а2г - агг )2 + 6ап 2 (31)
( ¡,] = г,в,г) (32)
В выражениях (31) и (32)
5,=<т,-а,, а,=С±;{<г,-а,)в*-1га,ё* Ц,] = Г,в,г), (33)
ад — компоненты тензора микронапряжений; <т0 и <т° - начальное и текущее значения предела текучести; Си/— начальное значение модуля кинематического упрочнения, и параметр, характеризующий его убывание с ростом пластических деформаций; 0„ и Ъ — соответственно величина максимально возможного изменения предела текучести материала среды и параметр, характеризующий быстроту его изменения в зависимости от скорости интенсивности пластических деформаций.
В качестве критериев, позволяющих определять переход материала в каждой точке среды из упругого состояния в пластическое, а также различать процессы нагружения и разгрузки при её упругопластическом поведении, приняты следующие:
а) в процессе нагружения; б) в процессе разгрузки.
—Лт„ >0; Ф = 0; £■'" >0 — Лт„ <0; Ф = 0; =0 (34)
да, дау
Как частный рассмотрен случай абсолютно жесткого индентора. На основе аналитического решения, полученного А. Динником для задачи о
12
статическом сжатии двух соприкасающихся тел при упругом поведении материала, дана оценка погрешности такого допущения. Представление индентора в виде твердого тела позволило сократить число уравнений в системах (26) и (28) до 2(1 + и соответственно понизить требования к объему вычислительных ресурсов и уменьшить время счета при получении численных результатов.
Конечноэлементная схема численного решения обоих вариантов вариационной задачи реализована с помощью программного комплекса АВАрШ. В силу локальности рассматриваемого процесса деформирования для оптимизации процесса вычислений при разбиении объектов континуума использована нерегулярная сетка: более густая вблизи области контакта и более редкая в областях, удаленных от контактирующих поверхностей. Расчеты проводились для двух типов конечных элементов - треугольных и четырехугольных, и двух типов функций формы - линейных и квадратичных двумя способами на основе методов Эйлера и Ньютона. В качестве начального принято невозмущенное состояние континуума. Точность расчетов при учете динамики нагружения континуума оказалась в значительной степени, зависящей от размера шага вычислений по времени. При его выборе как критерий использовано условие: отличие расчетных данных для двух соседних и всех последующих значений шага по времени не должно превышать 1% по отношению к предыдущему. С учетом динамики нагружения индентора наибольшее значение шага по времени, позволившее получить удовлетворительные результаты, составило 10"6 с. Для определения скоростей движения и перемещений точек континуума по найденным значениям их ускорений ввиду малости шага по времени использован закон равнопеременного движения. В качестве упрочняемого материала рассмотрена углеродистая сталь, а за материал рабочей поверхности индентора принят синтетический (технический) алмаз. Использование квадратичных функций формы обеспечило большую точность, чем линейных, но время счета при этом увеличилось. Процедура вычислений на базе метода Ньютона оказалась предпочтительней как по точности получаемых результатов, так и по затратам машинного времени. Использование сетки, состоящей из треугольных элементов и исходных уравнений на базе метода множителей Лагранжа, дало результат наиболее близкий к упомянутому выше аналитическому решению. В точках наибольших значений интенсивности напряжений отличие расчетных данных от аналитических находилось в пределах 1%.
В качестве примера на рис. 3 и 4 приведены графики поведения интенсивности напряжений с течением времени в центральной точке контактной поверхности и в точке её наибольших значений ( г = 0, 2 = ЮЛи) при практически мгновенном возрастании нагрузки, прилагаемой к абсолютно жесткому индентору, от нуля до заданного значения Г0 =19,62Н за промежуток времени Д/ = 10"6с.
Из них следует, что процесс деформирования носит колебательный характер и наибольшего значения интенсивность напряжений в центральной точке достигает в момент времени превышающий время возрастания прилагаемой к индентору нагрузки до заданного значения почти в 25 раз (рис. 36).
Рис. 3. Поведение интенсивности напряжений в центральной точке контакта (г = 0, г = 0) стечением времени
Рис. 4. Поведение интенсивности напряжений в точке наибольших значений (г = 0, г = 1 ■ 10"4 м) с течением времени
Затем амплитуда интенсивности напряжений убывает, колебательный процесс стремится к квазистационарному (рис. Зв) и с течением времени становится похожим на гармонический (рис. Зг). При заданном уровне нагружения интенсивность напряжений достигает максимального значения не в области контакта, а на некотором расстоянии от неё под центром давления (рис. 4), в то время как осевое напряжение о, имеет наибольшее значение на её поверхности в центре контура давления. Поэтому пластические деформации, впервые возникают не на поверхности контакта, а на некотором расстоянии от неё под центром давления, в точке с координатами г = 0, 7,5 • 10~5./и. Время, по истечении которого интенсивность напряжений впервые достигает предела текучести сг0 = 21 (Г На, составляет 2-Ю-3 с. Зарождаясь в этой внутренней точке среды, область пластических деформаций в течение промежутка времени 2 10"5-5-4-10~5с продолжает неравномерно расширяться по всем направлениям, достигая по его истечении предельного объема. Начиная с момента времени 4-КГ5с дальнейшие изменения полей напряжений и деформаций на уровень интенсивности пластических деформаций и форму занимаемой ими области влияния не оказывают. При этом граница области пластических деформаций до поверхности контакта не доходит, а зона наибольших значений их интенсивности смещается относительно точки их зарождения вглубь среды. Существенное влияние на процесс формирования пластической зоны оказывает скорость нагружения индентора до заданного постоянного значения силы 19,62//. При увеличении промежутка
времени нагружения А? качественная картина формирования зоны пластических деформаций и её форма сохраняются, а занимаемый её объем и уровень интенсивности пластических деформаций уменьшаются. При возрастании времени нагружения до А? =10 ~4с и Д? = 10~2с уровень наибольших значений интенсивности пластических деформаций падает до
2,443-10"4 и 1 • 1соответственно, что более чем в два и более чем в пять раз меньше, чем при нагружении за промежуток времени А? = 1(Г''с. При этом расстояние от поверхности контакта до упрочненного слоя во всех приведенных выше случаях практически остается тем же, равным 5 ■ 10~5л*.
С увеличением прикладываемого за тот же промежуток времени Лг = 10"6с к индентору усилия в два раза толщина поверхностного упругого слоя остается практически той же, область существования пластических деформаций растет в основном вглубь среды, точка наибольших установившихся значений интенсивности пластических деформаций смещается вниз и сами значения возрастают в 1,78 раза. Качественный характер поведения интенсивностей напряжений и деформаций с течением времени остается тем же, но вместе с тем появляется едва заметная тенденция к изменению формы области существования пластических деформаций. При дальнейшем увеличении задаваемых значений усилия происходят изменения как количественных, так и качественных характеристик рассматриваемого механического процесса. При значениях усилия в пределах 19,62Н - 117,72Н во всех точках области существования пластических деформаций наблюдается монотонный рост их уровня с течением времени без выхода её границы на поверхность среды. С превышением этого интервала объем занимаемый пластической зоной существенно изменяет свою форму, а сами пластические деформации выходят на поверхность среды в окрестностях границы пятна контакта.
Объем занимаемый областью пластических деформаций достигает своего предельного положения в довольно таки короткие промежутки времени порядка (4-5-14)- 10~5с, в связи с чем целесообразно использование импульсного характера силового воздействия в сочетании с последующим линейным перемещением индентора по обрабатываемой поверхности, с выбором при этом длительности импульса и интенсивности силы исходя из промежутка времени необходимого для формирования предельного положения пластической зоны и требуемой толщины упрочненного слоя. При изучении одиночных импульсных силовых воздействий длительностью КГ4с-И 0~3с показано, что после прекращения их действия значения интенсивности пластических деформаций стабилизируются, амплитуда колебаний интенсивности напряжений становится пренебрежимо малой.
Ввиду ограниченности объема реферата результаты численного анализа поведения среды при воздействии серии импульсов усилия, а также их влияния в сочетании с ультразвуковым воздействием здесь не приводятся.
В четвертой главе предложен упрощенный инженерный способ исследования напряженного деформированного состояния среды, в основе которого лежит физико-математическая модель изучаемого процесса, основанная на замене реального контактного воздействия рабочего инструмента на материал обрабатываемой среды, действием эквивалентной ему, виртуальной распределенной по определенному закону нагрузки (рис. 5).
q(r,t)
Рис. 5. Схема нагружения индентора и среды в рамках предлагаемой модели
Движение абсолютно жесткого индентора приведенной массы т под действием приложенной к нему активной силы F(t) предложено рассматривать отдельно от движения упруго-пластической среды, заменяя её действие на индентор равнодействующими реактивных N, (t) и диссипативныхМ, (t) сил контактного взаимодействия. При определении равнодействующей реактивных сил учтено, что рабочий орган может совершать только поступательное перемещение, обрабатываемая поверхность плоская и граница зоны контактного взаимодействия имеет форму окружности некоторого радиуса, а также то, что при вдавливании абсолютно жесткой сферы в упругое полупространство радиус границы области контакта а и распределение напряжений по её поверхности q(r) связаны с перемещением w центральной точки контакта в направлении оси Oz соотношениями:
: = -Jr. ■ w q(r) =
IE
wjl —
(35)
тсТЛ^-у7)" V аС учетом того, что при динамической постановке контактной задачи и упругопластическом поведении материала среды входящие в (35) характеристики контактного взаимодействия являются функциями времени, выражение для равнодействующей сил реакций связи получено в виде
W,(/) = \\q{r,t)rdrdq = cnpw2{t)-
(36)
где с р =4л/Л£/з([ -V2) - приведенный коэффициент жесткости материала среды, Е — его модуль упругости при отсутствии пластических деформаций или секущий модуль при их наличии, V - коэффициент Пуассона. Для учета влияния рассеяния энергии при циклическом упругопластическом деформировании среды на движение индентора принята условная вязкоупругая модель Кельвина — Фойхта с коэффициентом вязкости, энергетически эквивалентным материалу среды в отношении его поглощающих свойств, согласно которой выражение для равнодействующей диссипативных сил представлено в виде
АГ2=-Ьнр*. (37)
Здесь bnp - ~jRw0 E\\>/яю (l - v2 ) - приведенный коэффициент вязкости материала среды, у - его коэффициент поглощения энергии, w0 - некоторое статическое смещение, которое получил бы индентор под действием постоянной составляющей приложенной к нему активной силы F0, а о -круговая частота колебаний её переменной составляющей. С учетом принятых допущений и соотношений (37), (38) изучение движения индентора при нулевых начальных условиях сведено к решению следующей краевой задачи Коши:
3
mw + bnpw + cnpw2(t) = F(t)
w = 0 , w - О при t = 0.
Присутствие в левой части уравнения (38) слагаемого, содержащего перемещение центральной точки контакта индентора в дробной степени с четным знаменателем подтверждает необходимость наличия в обязательном порядке в составе равнодействующей активных сил F(t) стационарной (постоянной) составляющей F0, обеспечивающей необходимый рабочий натяг между рабочим инструментом и обрабатываемой средой.
Найденные путем решения краевой задачи (38) значения перемещения центральной точки рабочей поверхности индентора w(f) являются исходными данными для определения в каждый момент времени интенсивности распределенной виртуальной эквивалентной нагрузки q(r,t) и радиуса её приложения a(t) к среде по формулам (35), причем при идеально упругом поведении материала среды задача о движении индентора рассматривается независимо от задачи о напряженно-деформированном состоянии среды. При появлении в среде пластических деформаций обе задачи решаются совместно в связи с необходимостью перерасчета секущего модуля материала среды в задаче о движении индентора на каждом шаге вычислений.
Фактическое существование периодических решений «существенно» нелинейной задачи о движении индентора (38), описывающих состояние колебательного равновесия системы в окрестностях некоторого стационарного (статического) положения подтверждено результатами исследования их устойчивости с применением метода фазовой плоскости и теоремы Ляпунова. При исследовании движения индентора обнаружено явление нелинейного многопикового резонанса в диапазоне частот, кратных частоте собственных колебаний соответствующей (38) линеаризованной системы.
Задача о напряженно-деформированном состоянии среды в классической постановке аналогично тому, как это было сделано в главах 2 и 3 диссертации, сведена к решению системы обыкновенных дифференциальных уравнений первого порядка для определения скоростей в узловых точках ее разбиения на элементы:
Я, h
m=l at 0X!,(r,t)
= + } г, О^г (39)
Соответствие предложенной упрощенной инженерной модели исходной задаче подтверждено совпадением полученных с ее помощью результатов с аналитическим решением А. Н. Динника в случае статического нагружении индентора при упругом поведении материала среды. Определены пределы применимости предложенной модели с точки зрения сокращения времени счета. Установлено, что упрощенную модель целесообразно использовать при динамическом нагружении и идеально упругом состоянии материалы среды или при статическом нагружении и упругопластическом его поведении. При динамическом нагружении индентора и учете появления пластических деформаций применение упрощенной модели нецелесообразно в виду необходимости многократного перерасчета секущего модуля в приведенном коэффициенте жесткости и, как следствие эквивалентной нагрузки, на каждом временном шаге вычислений.
В заключении приведены основные результаты диссертации.
В приложении представлены документы, подтверждающие использование результатов работы ООО «Научно-производственное предприятие нестандартных изделий машиностроения» и в учебном процессе.
Основные результаты работы:
1) Сформулирована математическая модель процесса поверхностного упрочнения металлического изделия «алмазным карандашом» в виде континуума, состоящего из цилиндрического упругого жесткого индентора со сферической рабочей поверхностью оказывающего силовое воздействие на упругопластическую среду. На её базе в рамках теории течения поставлена динамическая контактная упругопластическая задача о напряженно-деформированном состоянии континуума в виде, позволяющем учесть нарушение первоначальной геометрии деформируемых вследствие контакта поверхностей без наложения ограничений на их форму.
2) На основе принципа Даламбера - Лагранжа с применением методов штрафных функций и множителей Лагранжа построен вариационный аналог поставленной задачи. Вариационные уравнения движения континуума выведены без наложения, каких-либо ограничений на форму деформированных поверхностей контактирующих тел, физические уравнения их материалов и
геометрические соотношения.
3) С использованием метода конечных элементов исследование НДС континуума сведено к определению полей скоростей узловых точек путем решения системы из 2(п, +п„) обыкновенных дифференциальных уравнений при использовании метода штрафных функций и смешанной системы, включающей в себя 2(п, +п„) обыкновенных дифференциальных уравнений и пс нестрогих алгебраических неравенств, обеспечивающих безусловный минимум соответствующих вариационных функционалов.
4) Проведен численный анализ напряженно-деформированного состояния упругопласгической среды при различных вариантах приложения нагрузки к индентору (статическая, динамическая, импульсная, с приложением
ультразвуковых воздействий). На его основе сформулированы рекомендации по выбору силовых воздействий, обеспечивающих при заданном уровне упрочнения сохранение по возможности однородного поверхностного упругого слоя, играющего роль своего рода амортизатора, снижающего негативное влияние динамических нагрузок на готовое изделие в процессе эксплуатации.
5) Разработана упрощенная инженерная расчетная модель изучаемого процесса, позволяющая проводить численные исследования напряженно-деформированного состояния среды с меньшими затратами вычислительных ресурсов. С учетом физических свойств материала среды и возможных вариантов нагружения индентора даны рекомендации по целесообразности её применения с точки зрения уменьшения временных затрат при проведении вычислений.
Основные положения диссертации опубликованы в следующих работах:
В изданиях, рекомендованных перечнем ВАК РФ
1. Овчинникова Н. В. Модельная задача для исследования процессов поверхностного упрочнения пластическим деформированием с применением ультразвуковых воздействий / Н. В. Овчинникова, Д. Г. Павлов, Ю. В. Чеботаревский // Вестник Саратовского государственного технического университета. - 2007. -ЛЬ 4 (28). - Вып. 1.-С. 14-18.
2. Овчинникова Н. В. О некоторых особенностях применения метода конечных элементов к решению контактной задачи на базе программного комплекса АВАСШБ / Н. В. Овчинникова, Ю. В. Чеботаревский // Известия Саратовского университета. Новая серия. Серия: Математика. Механика Информатика. - 2009. - Т. 9. - Вып. 2. - С. 82-88.
3. Овчинникова Н. В. К расчету напряженно-деформированного состояния упругопластического полупространства, контактирующего с абсолютно жестким индентором / Н. В. Овчинникова, Ю. В. Чеботаревский // Вестник Саратовского государственного технического университета - 2010 -№4(51). -Вып. 3. - С. 10-17.
4. Овчинникова, Н. В. О движении абсолютно жесткого индентора, взаимодействующего с упругопластической средой / Н. В. Овчинникова, Ю. В. Чеботаревский // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2010. -№ 3 (15). — С. 151-164.
5. Овчинникова Н. В. Вариационное уравнение движения континуума «жесткий индентор - деформируемая среда» / Н. В. Овчинникова, Ю. В. Чеботаревский // Вестник Саратовского государственного технического университета. - 2011. - №4(60). - Вып. 2. - С. 48-57.
6. Овчинникова Н. В. Применение метода множителей Лагранжа к решению контактной задачи о взаимодействии деформируемой среды с относительно жестким индентором / Н. В. Овчинникова, Ю. В. Чеботаревский // Вестник Саратовского государственного технического университета. - 2012. -№4(68). - Вып. 2. - С. 36-43.
7. Овчинникова Н. В. Применение метода штрафных функций к выводу вариационного уравнения движения континуума «индентор - деформируемая среда» / Н. В. Овчинникова, Ю. В. Чеботаревский // Вестник Самарского государственного технического университета. - Сер. Технические науки -2013.-№ 1 (37).-С. 127-134.
В других изданиях
8. Овчинникова Н. В. К решению задачи о поверхностном упрочнении металлов с применением ультразвука / Н. В. Овчинникова, Д. Г. Павлов, Ю. В. Чеботаревский // Динамические и технологические проблемы механики конструкций и сплошных сред: материалы XV Междунар. симпозиума им.
A.Г.Горшкова, Ярополец, 16-20 февраля 2009 г. - М.: Изд-во «Типография «ПАРАДИЗ», 2009. - Т. 1. - С. 121-122.
9. Овчинникова Н. В. Динамика абсолютно жесткого индентора, взаимодействующего с упругопластической средой / Н. В. Овчинникова,
B. Э. Джашитов, Ю. В. Чеботаревский // Динамические и технологические проблемы механики конструкций и сплошных сред: материалы XVI Междунар. симпозиума им. А. Г. Горшкова, Ярополец, 15-19 февраля 2010 г. - Чебоксары: ГУП «ИГПС «Чувашия», 2010. - Т. 1. - С. 63-65.
10. Овчинникова Н. В. Об одном аналитическом решении задачи о нелинейных колебаниях материальной точки // Инновации и актуальные проблемы техники и технологий: материалы Всерос. науч.-практ. конф. молодых ученых, Саратов, 26-29 октября 2010 г. Саратов: Сарат. гос. техн. ун-т, 2010.-Т. 1,-С. 201-203.
11. Овчинникова Н. В. Расчет напряженно-деформированного состояния упругопластической среды, контактирующей с абсолютно жестким индентором / Н. В. Овчинникова, Ю. В. Чеботаревский // Инновации и актуальные проблемы техники и технологий: материалы Всерос. науч.-практ. конф. молодых ученых, Саратов, 26-29 октября 2010 г. Саратов: Сарат. гос. техн. ун-т, 2010. - Т. 2. - С. 122-125.
12. Овчинникова Н.В. Инженерный способ расчета НДС упругопластической среды, контактирующей с абсолютно жестким индентором / Н. В. Овчинникова, Ю. В. Чеботаревский // Участники школы молодых ученых и программы У.М.Н.И.К.: сб. тр. XXIV Междунар. науч. конф. ММТТ-24. -Саратов, 2011. - С. 234-235.
Подписано в печать 17.10.13 Формат 60x84 1/16
Бум. офсет. Усл. печ. л. 1,0 Уч.-изд. л. 1,0
Тираж 100 экз. Заказ 160 Бесплатно
Саратовский государственный технический университет
410054, Саратов, Политехническая ул., 77 Отпечатано в Издательстве СГТУ. 410054, Саратов, Политехническая ул., 77 Тел.: 24-95-70; 99-87-39, e-mail: izdat@sstu.ru
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Саратовский государственный технический университет имени Гагарина Ю.А.»
На правах рукописи
04201451247
Овчинникова Наталья Владимировна
НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ КОНТИНУУМА «ЖЕСТКИЙ ИНДЕНТОР - УПРУГОПЛАСТИЧЕС*САЯ СРЕДА» ПРИ ДИНАМИЧЕСКИХ НАГРУЖЕНИЯХ
Специальность 01.02.04 - Механика деформируемого твердого тела
Диссертация
на соискание ученой степени кандидата технических наук
Научный руководитель: доктор технических наук, профессор Чеботаревский Ю. В.
Саратов -2013
Оглавление
Оглавление................................................................................... 2
Введение..................................................................................... 4
Глава 1. Задача о напряженно-деформированном состоянии континуума «индентор-деформируемая среда»
1.1. Обоснование выбора направления исследования.............................. 11
1.2 Постановка динамической задачи о напряженно-деформированном состоянии континуума «жесткий индентор - упругопластическая среда»... 20
Глава 2. Задача о движении континуума «индентор-деформируемая среда» в вариационной постановке
2.1 Вариационное уравнение движения континуума «жесткий индентор -
деформируемая среда».................................................................... 27
2.2 Применение метода штрафных функций для учета ограничений в вариационном уравнении движения континуума................................... 35
2.3 Использование метода множителей Лагранжа для учета ограничений в вариационном уравнении движения континуума................................... 45
Глава 3. Исследование напряженно-деформированного состояния «жесткий индентор - упругопластическая среда»
3.1 Применение метода конечных элементов для решения задачи в вариационной постановке................................................................ 51
3.2 Выбор рациональной схемы численного решения поставленной задачи
с использованием программного комплекса АВАС>и8........................... 61
3.3 Исследование напряженно - деформированного состояния упругопластической среды, контактирующей с абсолютно жестким индентором... 67
Глава 4. Решение задачи о движении континуума «жесткий индентор -
упругопластическая среда» на основе упрощенной модели
4.1 Упрощенная физико-математическая модель взаимодействия жесткого
индентора и упругопластической среды..........................................................................................84
4.2 Динамика абсолютно твердого индентора, контактирующего с деформируемой средой....................................................................................................................................87
4.3 Задача о напряженно-деформированном состоянии упругопластической среды под действием эквивалентной нагрузки............................................................................97
4.4 Определение пределов применимости упрощенной модели..................................100
Заключение....................................................................................................................................................................106
Список использованной литературы......................................................................................................108
Приложение А............................................................................................................................................................121
Введение
В настоящее время активно осваиваются различные новые подходы к повышению эксплуатационных характеристик деталей машин за счет упрочнения их рабочих поверхностей различными способами. К их числу принадлежат технологии, использующие для этих целей механизмы нестандартных динамических воздействий различной физической природы, как механической, так и не механической. К наиболее распространенным методам поверхностного упрочнения с немеханическими источниками энергии относятся: термохимический, электронно-лучевой, плазменный, лазерный, ионная имплантация, детонационное нанесение покрытий, а также термомагнитная и импульсная обработка в магнитном поле. Среди «статических» механических методов поверхностного упрочнения обрабатываемого объекта на сегодняшний день наибольшее распространение получили: калибровка отверстий, алмазное выглаживание, дорнование, обкатка (раскатка) и многие другие. Наряду с ними, в машиностроении достаточно широко используются методы поверхностного пластического деформирования, основанные на ударном динамическом воздействии рабочего инструмента на обрабатываемую поверхность. К ним относятся: виброударная, ультразвуковая, дробеструйная, центробежно -шариковая обработки, чеканка, виброобкатка и вибровыглаживание.
К сожалению, ни один из перечисленных методов не является универсальным, способным дать положительный эффект при всех возможных условиях эксплуатации упрочненных деталей готового изделия. Поэтому в последнее время для улучшения качества обработки получило применение комбинированных методов упрочнения, объединяющих одновременно несколько способов силового воздействия на обрабатываемый материал. В частности, один из них, основан на использовании в процессе поверхностного упрочнения металлов импульсного нагружения рабочего инструмента в сочетании с ультразвуковым на него воздействием. Применение такого режима обработки позволяет существенно повысить качество обрабатываемой поверхности и
обеспечить более высокую износоустойчивость готовых деталей. В её процессе проявляются новые динамические эффекты воздействия ультразвука на материал, прогнозирование которых без проведения предварительных аналитических и численных исследований весьма затруднительно. Для осознанного управления технологическим процессом поверхностного упрочнения необходимо определить комплекс геометрических, механических и физических параметров, его характеризующих и исследовать их влияние на результаты обработки. Механические процессы, происходящие при взаимодействии рабочего инструмента и поверхности детали, настолько сложны и разнообразны, что учесть их все специфические особенности в рамках решения единой задачи не представляется возможным. Одними из косвенных характеристик, позволяющими судить об уровне упрочнения, могут служить возникающие в его процессе остаточные пластические напряжения и деформации в обрабатываемом материале, а также его плотность, являющаяся функцией интенсивности остаточных деформаций. В связи с этим весьма актуальным представляется исследование напряженно-деформированного состояния упрочняемого материала при комбинированном динамическом силовом и ультразвуковом воздействиях.
В основе механического упрочнения пластическим деформированием лежат процессы силового воздействия рабочего инструмента на поверхность обрабатываемого материала. Поэтому с точки зрения механики сплошных сред задача об исследовании напряженно деформированного состояния упрочняемого материала является контактной. Кроме того, с учетом специфики прикладываемого к рабочему инструменту силового воздействия она должна рассматриваться в динамической постановке с учетом упругопластического характера поведения обрабатываемого материала. Однако данный класс задач механики деформируемого твердого тела до сих пор является мало изученным. Поэтому разработка новых подходов к решению контактных задач для упругопластических тел, находящихся под воздействием описанных выше сложных динамических нагрузок, представляется весьма актуальной.
С учетом этого целью настоящей диссертационной работы является разработка новых подходов к решению задачи о напряженно-деформированном состоянии континуума «жесткий индентор - упругопластическая среда» при различных видах статического и динамического нагружения индентора и неупругом поведении материала среды.
Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка использованной литературы и приложения.
В первой главе проводится обзор и анализ состояния рассматриваемой проблемы, обосновывается выбор предмета и направления исследований. Для принятой математической модели процесса упрочнения формулируется осесимметричная динамическая контактная задача о напряженном и деформированном состоянии континуума «жесткий индентор - деформируемая среда». Постановка задачи проводится в виде, позволяющем учесть нарушение первоначальной геометрической формы деформируемых в процессе контакта поверхностей и появление необратимых деформаций в среде.
Во второй главе на основе принципа Даламбера - Лагранжа выводится исходное вариационное уравнение, описывающее движение континуума. Для введения в него ограничений, накладываемых на движение тел составляющих континуум, и сведения решения поставленной задачи к поиску безусловного минимума соответствующего функционала используются методы штрафных функций и множителей Лагранжа.
В третьей главе с применением метода конечных элементов на базе полученных вариационных соотношений проводится дискретизация исходной задачи. Её численное решение сводится к отысканию либо скоростей движения узловых точек разбиения континуума из системы обыкновенных дифференциальных уравнений, либо их ускорений из системы алгебраических уравнений. Формулируются условия и выводятся уравнения для определения границы области контакта в случаях применения метода штрафных функций и метода множителей Лагранжа. Проводится выбор рациональной схемы
численного решения поставленной задачи с использованием программного комплекса АВА(^и8. Приводятся результаты численных расчетов и их анализа.
В четвертой главе на основе предлагаемой физико - математической модели разрабатывается упрощенный инженерный подход к исследованию напряженно - деформированного состояния континуума, основанный на замене реального взаимодействия контактирующих тел воздействием на среду эквивалентной ему виртуальной нагрузки. Исследуется динамика абсолютно твердого индентора контактирующего с упругопластической средой. Выводится вариационное уравнение движения среды под действием эквивалентной нагрузки. С применением метода конечных элементов проводится дискретизация динамической упругопластической задачи о движении среды. Определяются границы рациональной применимости предложенной упрощенной модели.
В заключении формулируются основные выводы, обобщающие результаты проведенных исследований.
В приложении приводятся документы, подтверждающие использование полученных результатов в научно - производственной и образовательной деятельности.
Научная новизна работы заключается в следующем:
- выполнена постановка динамической упругопластической задачи о напряженно - деформированном состоянии континуума «жесткий индентор -деформируемая среда» в виде, не накладывающем ограничений на форму области контакта и позволяющем учесть нарушение первоначальной геометрической формы деформируемых в процессе контакта поверхностей;
- на основе принципа Даламбера - Лагранжа с применением методов штрафных функций и множителей Лагранжа выведены вариационные уравнения движения континуума без наложения, каких - либо ограничений на форму деформированных поверхностей контактирующих тел, физические уравнения их материалов и геометрические соотношения;
- с целью применения метода конечных элементов для получения численных результатов проведена дискретизация рассматриваемой задачи в вариационной постановке;
- получены уравнения для определения координат точек, лежащих на границе контактной области в случае применения метода штрафных функций и при использовании метода множителей Лагранжа;
- на базе физической модели изотропно кинематического упрочнения материала среды проведено численное исследование напряженно -деформированного состояния тел образующих континуум при различных видах приложения нагрузки к индентору (статическая, динамическая, импульсная, с приложением ультразвуковых воздействий).
- сформулирована упрощенная инженерная расчетная модель изучаемого процесса контактного взаимодействия тел, образующих континуум «абсолютно твердый индентор - деформируемая среда», позволяющая в некоторых частных случаях нагружения индентора исследовать напряженно - деформируемое состояние среды при упругом и упругопластическом поведении её материала с меньшими затратами вычислительных ресурсов;
- определены границы применимости предложенной упрощенной расчетной модели континуума «абсолютно твердый индентор - деформируемая среда» для исследования напряженно - деформированного состояния среды при различных видах нагружения индентора.
Достоверность полученных результатов обеспечивается корректностью математической постановки рассматриваемых задач, использованием известных аналитических и численных методов, в том числе метода конечных элементов, а также их совпадением с известными аналитическими решениями в некоторых частных случаях. Полученные численные результаты качественно
подтверждаются имеющимися в периодической литературе экспериментальными данными и соответствуют физике исследуемых процессов.
Практическая ценность работы состоит в использовании её результатов при проектировании технологического оборудования в ООО «Научно -производственное предприятие нестандартных изделий машиностроения» и в учебном процессе при подготовке инженеров по специальности №№№ «Проектирование технологических комплексов» в Саратовском государственном техническом университете.
Апробация работы. Основные положения и результаты работы докладывались на 5 научно - технических конференциях и 3 научных семинарах, в том числе: на XV Международном симпозиуме «Динамические и технологические проблемы механики конструкций и сплошных сред» им. А. Г. Горшкова (Ярополец, 16-20 февраля 2009г.); на Всероссийской научно-технической конференции «Совершенствование техники, технологий и управления в машиностроении» (Саратов, 20-24 октября 2009г.); на XVI Международном симпозиуме «Динамические и технологические проблемы механики конструкций и сплошных сред» им. А. Г. Горшкова (Ярополец, 15-19 февраля 2010г.); на Всероссийской научно-практической конференции молодых ученых «Инновации и актуальные проблемы техники и технологий», (Саратов 2629 октября 2010г.); на XXIV Международной научной конференции «Математические методы в технике и технологиях» ММТТ-24 (Саратов, 2011г.).
Публикации. Основные результаты работы опубликованы в 12 печатных работах, в том числе 7 работ в изданиях рекомендованных ВАК РФ для публикации материалов кандидатских и докторских диссертаций.
На защиту выносятся следующие основные положения:
1. Вариационные уравнения движения континуума «жесткий индентор -деформируемая среда», полученные на основе вариационного принципа Даламбера - Лагранжа с применением методов штрафных функций и множителей Лагранжа и позволяющие учесть искривление деформируемых вследствие контакта поверхностей без наложения ограничений на их форму и физические
уравнения материалов контактирующих тел (упругое, вязкоупругое, у пру гопл астическое);
2. Упрощенная физико - математическая модель процесса контактного взаимодействия тел, образующих континуум «абсолютно твердый индентор -деформируемая среда», позволяющая исследовать напряженно - деформируемое состояние среды с меньшими затратами вычислительных ресурсов:
- при динамическом нагружении индентора и идеально упругом поведении материала среды;
- при статическом нагружении индентора и упругопластическом поведении материала среды.
3. Результаты численного исследования напряженно - деформированного состояния тел, образующих континуум, полученные для физической модели изотропно кинематического упрочнения материала среды при различных видах приложения нагрузки к индентору (статическая, динамическая, импульсная, с приложением ультразвуковых воздействий).
Глава 1. Задача о напряженно-деформированном состоянии континуума «индентор-деформируемая среда»
1.1 Обоснование выбора направления исследования
Подавляющее большинство технологических задач обработки металлов давлением с точки зрения механики сплошных сред относятся к контактным. Во многих случаях при рассмотрении этих задач о взаимодействии двух тел локальными явлениями в зоне контакта пренебрегают и саму задачу приводят к смешанной задаче теории упругости с граничными условиями в виде заданных усилий и перемещений. Однако такое упрощение не всегда оправдано. В реальных технологических процессах закон распределения истинных контактных давлений оказывает существенное влияние на напряженно-деформированное состояние взаимодействующих тел. В этих случаях необходимо учитывать, что размеры и форма зоны контакта, а также условия взаимодействия на ней нелинейно зависят от приложенной нагрузки и могут быть определены только в процессе решения задачи. При больших значениях внешней нагрузки напр
