Министерство образования и науки Российской Федерации
Федеральное государственное бюджетное образовательное учреждение
высшего образования
Пермский национальный исследовательский политехнический
университет
Факультет прикладной математики и механики
Кафедра «Динамика и прочность машин»
Направление подготовки: 15.04.03 Прикладная механика
Направление (профиль) образовательной программы:
Динамика и прочность машин, конструкций и механизмов
ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ РАБОТА
(магистерская диссертация)
На тему Особенности распределения остаточных напряжений
и перемещений в упругих телах, вызванных их поверхностной обработкой
Студент Виндокуров Илья Владимирович
Состав ВКР:
1. Пояснительная записка на 38 страницах
2. Графическая часть состоит из 19 рисунков
Допускается к защите
Заведующий кафедрой
д.т.н., профессор
В.П. Матвеенко
« 21 »
июня
2019 г.
Регистрационный номер
Пермь — 2019
Руководитель ВКР
(И. Э. Келлер)
Министерство образования и науки Российской Федерации
Федеральное государственное бюджетное образовательное учреждение
высшего образования
Пермский национальный исследовательский политехнический
университет
Кафедра «Динамика и прочность машин»
УТВЕРЖДАЮ
Зав. кафедрой
В.П. Матвеенко
«5»
июня
2019 г.
ЗАДАНИЕ
на выполнение выпускной квалификационной работы
Фамилия, И. О. Виндокуров Илья Владимирович
Факультет прикладной математики и механики Группа ДПМ-17-1м
Начало выполнения работы
Контрольные сроки просмотра работы кафедрой
Сроки представления на рецензию
Защита работы на заседании ГЭК 27.06.2019
1. Наименование темы Особенности распределения остаточных
напряжений и перемещений в упругих телах, вызванных их поверхностной
обработкой
2. Исходные данные к работе
3. Содержание пояснительной записки
а) Математическая постановка задачи и ее численная формулировка;
б) Исследование стрелы прогиба полосы в зависимости от ее ширины;
в) Исследование особенностей распределения остаточных напряжений
в пластине с локальной областью поверхности, подвергнутой
лазерной проковке;
г) Исследование особенностей распределения остаточных напряжений
вблизи кромки лопатки;
д) Заключение;
4. Дополнительные указания
5. Основная литература
а) Биргер И.А. Остаточные напряжения. М.: Машгиз, 1963. 232с.
Руководитель выпускной квалификационной работы магистра
д.ф.-м.н., профессор кафедры ДПМ
И. Э. Келлер
Задание получил
И. В. Виндокуров
КАЛЕНДАРНЫЙ ГРАФИК ВЫПОЛНЕНИЯ ВЫПУСКНОЙ
КВАЛИФИКАЦИОННОЙ РАБОТЫ
№
1
2
3
4
5
6
Наименование этапа
Разработка основных разделов
диссертации
Оформление диссертации
Разработка и оформление
иллюстративного материала
к защите диссертации
Представление диссертации на
проверку и отзыв научного
руководителя
Представление работы
заведующему кафедрой
Защита на заседании ГЭК
Научный руководитель работы
«5»
июня
2019 г.
Объем, Сроки выполнения Приме
начало
конец
%
чание
75
05.06.19
14.06.19
15
15.06.19
17.06.19
3
18.06.19
19.06.19
5
16.06.19
18.06.19
1
20.06.19
1
27.06.19
27.06.19
(И. Э. Келлер)
Реферат
Современные регламенты проектирования авиационного двигателя со
держат требования прямого учета в расчетах на прочность поверхностного
слоя тангенциальных сжимающих остаточных напряжений. Данный слой
обычно наводится процессом холодной обработки детали путем бомбарди
ровки участка ее поверхности мелкими стальными шариками. Упругость
материала стремится восстановить начальную форму детали вблизи ее
поверхности, что приводит к образованию в поверхностном слое поля
тангенциальных сжимающих остаточных напряжений и необходимого для
равновесия поля остаточных напряжений в детали вне этого слоя. В рабо
те численно исследуются поля остаточных напряжений и перемещений в
некоторых задачах, важных с точки зрения приложений к проектированию
перспективных технологий поверхностной обработки ответственных деталей
авиационного двигателя и методов разрушающего контроля остаточных на
пряжений.
Целью данной работы является:
1. Исследование прогиба пластины, поверхностно упрочненной в зави
симости от ее ширины.
2. Исследование качественного и количественного сравнения распреде
ления остаточных напряжений по глубине в однородно и неоднородно
обработанной пластине.
3. Исследование концентрации растягивающих напряжений в поверх
ностно обработанном затупленном клине.
Существует инженерная проблема, с которой сталкиваются при произ
водстве. Во время работы с различными пластинами возникают вопросы, как
правильно подобрать ширину для различных прочностных испытаний и как
ширина влияет на прогиб. Оказалась что от ширины прогиб не меняется ес
ли его мерить по длине в центральном сечении. Численное моделирование
проводилось в программном комплексе ANSYS APDL. Численный расчет
показал, что пластины выгибаются в направлении границы, вблизи кото
рой заданы пластические деформации. По данным численного эксперимента
была построена зависимость стрелы прогиба от ширины пластины. Увеличе
ние долговечности и надёжности лопаток в двигателе возможна различными
способами, одним из которых является поверхностное упрочнение. В связи
со сложной геометрией и высокими требованиями к качеству поверхности
имеется очень ограниченная номенклатура способов поверхностного пласти
ческого деформирования. Одним из таких методов является бомбардировка
участка поверхности мелкими стальными шариками. Такая обработка при
водит к пластическому растяжению материала в касательной плоскости к
обрабатываемой поверхности в тонком приповерхностном слое. Заданный ис
точник пластических деформаций создает поле сжимающего напряжения по
контуру, а также поле растягивающего напряжения во внутренней области
затупленного клина. Сходимость по максимальным напряжениям показала,
что количество элементов были подобраны верно для оптимального времени
расчета и точности результатов. Основные результаты работы заключают
ся в следующем:
1. Подтверждена независимость стрелы прогиба полосы, подвергну
той односторонней однородной поверхностной обработке пласти
ческим деформированием, от ее ширины. Данный факт важен
при лабораторном исследовании остаточных напряжений методом
Н.Н.Давиденкова, в особенности нестандартных образцов, разработ
ке нового оборудования для этих задач.
2. Показана возможность оценки глубины и максимальной величины
остаточных напряжений с использованием обработки части поверх
ности контрольной пластины. Ошибка по максимальной величине
остаточных напряжений не превышает 8%, при этом для частично
обработанной пластины удовлетворяется с хорошим приближением
интегральное условие равновесия. Данный результат обосновывает
упрощенную методику исследования параметров профиля остаточ
ных напряжений, наведенных различными технологиями обработки.
3. Показано образование опасного концентратора растягивающих на
пряжений вблизи кромки затупленного клина при его поверхностной
обработке пластическими деформациями, способного приводить к
пластическому разрыхлению и образованию усталостных трещин в
данной области.
Оглавление
Стр.
Введение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9
Глава 1. Математическая постановка задачи и ее численная
формулировка . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Глава 2. Исследование стрелы прогиба полосы в
зависимости от ее ширины . . . . . . . . . . . . . . . .
2.1 Постановка задачи . . . . . . . . . . . . . . . . . . . . . . . .
2.2 Численное моделирование в программном комплексе ANSYS
APDL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.3 Результаты численного решения . . . . . . . . . . . . . . . .
2.4 Вывод по главе 2 . . . . . . . . . . . . . . . . . . . . . . . . .
Глава 3. Исследование особенностей распределения
остаточных напряжений в пластине с локальной
областью поверхности, подвергнутой лазерной
проковке . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.1 Постановка задачи . . . . . . . . . . . . . . . . . . . . . . . .
3.2 Численное моделирование в программном комплексе ANSYS
APDL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.3 Результаты численного решения . . . . . . . . . . . . . . . .
3.4 Вывод по главе 3 . . . . . . . . . . . . . . . . . . . . . . . . .
Глава 4. Исследование особенностей распределения
остаточных напряжений вблизи кромки лопатки .
4.1 Постановка задачи . . . . . . . . . . . . . . . . . . . . . . . .
4.2 Численное моделирование в программном комплексе ANSYS
APDL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.3 Результаты численного решения . . . . . . . . . . . . . . . .
4.4 Вывод по главе 4 . . . . . . . . . . . . . . . . . . . . . . . . .
. . 14
. . 14
. . 17
. . 18
. . 20
. . 21
. . 21
. . 23
. . 25
. . 27
. . 28
. . 28
. . 29
. . 31
. . 33
Стр.
Заключение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Список литературы
. . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Список рисунков . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Введение
Современные регламенты проектирования авиационного двигателя со
держат требования прямого учета в расчетах на прочность поверхностного
слоя тангенциальных сжимающих остаточных напряжений [1]. Данный слой
обычно наводится процессом холодной обработки детали путем бомбарди
ровки участка ее поверхности мелкими стальными шариками. Удар каждой
дробинки оставляет на поверхности микроскопическую ямку и в целом по
крытие такими ямками приводит к пластическому растяжению материала в
касательной плоскости к обрабатываемой поверхности в тонком приповерх
ностном слое (при обработке пластины в умеренном количестве это ведет к
ее изгибу с выпуклостью со стороны обработанной поверхности [2]).
Упругость материала стремится восстановить начальную форму детали
вблизи ее поверхности, что приводит к образованию в поверхностном слое по
ля тангенциальных сжимающих остаточных напряжений и необходимого для
равновесия поля остаточных напряжений в детали вне этого слоя. Поскольку
усталостное и коррозионное разрушение, водородное охрупчивание, разруше
ние вследствие износа, выкрашивания, образования задиров и кавитационной
эрозии деталей машин начинается с поверхности, такая ее обработка способ
ствует увеличению их ресурса.
Рассматриваемые технологические процессы поверхностного упрочне
ния разнообразны, чаще всего используются пневмо и гидродробеструйная
обработка [3], дробеструйная обработка и обкатка в сочетании ультразвуком
[4; 5], технологии с применением кавитации [6; 7], обкатка роликом, алмазное
выглаживание. Указанные технологии эффективно применяются для литых
или спеченных из порошка высокопрочных сталей, алюминиевых и титано
вых сплавов [8; 9].
В работе численно исследуются поля остаточных напряжений и пе
ремещений в некоторых задачах, важных с точки зрения приложений к
проектированию перспективных технологий поверхностной обработки от
ветственных деталей авиационного двигателя и методов разрушающего
контроля остаточных напряжений. Использование прямого моделирования
технологического процесса, посредством которого создаются эти остаточные
напряжения, необходимо для исследования закономерностей их формиро
9
вания и выбора рациональных режимов, но для внесения остаточных
напряжений в численную модель прогнозирования ресурса изделия выгля
дит избыточным средством. Для постановки подобных задач вместо расчета
предшествующего процесса поверхностной обработки разумно задавать рас
пределение по глубине приповерхностного слоя пластических деформаций,
являющихся источником остаточных напряжений в упругом теле. В отличие
от поля остаточных напряжений пограничный слой пластических деформа
ций вблизи свободной границы не изменяется при вырезке полосы из изделия,
не должен удовлетворять каким-либо полевым уравнениям и потому предо
ставляет удобный способ введения остаточных напряжений в численную
модель.
Применительно к пластине, подвергнутой односторонней однородной в
ее плоскости обработке, проблема реконструкции распределения пластиче
ской деформаций известна. Все известные варианты ее решения сводятся к
некорректным задачам. При этом используется решение задачи об остаточ
ных напряжениях при известном распределении пластических деформаций
[10—13] либо решению задачи о распределении пластических деформаций в
упруго пластическом полупространстве, вызванном ударом одного шарика
цитата [14—19]. Таким образом, приведенные здесь решения задач пред
ставляются важным элементом методики учета остаточных напряжений в
прогнозных расчетах чувствительных к ним характеристик прочности изде
лий, подвергнутых поверхностной обработке.
Целью данной работы является
1. Исследование прогиба пластины поверхностно упрочненной в зави
симости от ее ширины.
2. Исследование качественного и количественного сравнения распреде
ления остаточных напряжений по глубине в однородно и неоднородно
обработанной пластине.
3. Исследование концентрации растягивающих напряжений в поверх
ностно обработанном затупленном клине.
Объем и структура работы. Работа состоит из введения, двух глав
и заключения. Полный объём работы составляет 38 страниц, включая 19 ри
сунков и 0 таблиц. Список литературы содержит 20 наименований.
10
Глава 1. Математическая постановка задачи и ее численная
формулировка
Рассматривается равновесие свободного от действия внешних сил изо
тропного линейно-упругого тела с полем пластических деформаций, подчи
няющееся уравнениям:
𝜕𝑖 σ𝑖𝑗 = 0,
(1.1)
σ𝑖𝑗 = λε𝑒𝑘𝑘 δ𝑖𝑗 + 2µε𝑒𝑖𝑗 , ε𝑒𝑖𝑗 = ε𝑖𝑗 − ε𝑝𝑖𝑗 , ε𝑝𝑘𝑘 = 0,
(1.2)
1
(1.3)
ε𝑖𝑗 = (𝜕𝑖 𝑢𝑗 + 𝜕𝑗 𝑢𝑖 ),
2
где фигурируют компоненты σ𝑖𝑗 симметричного тензора напряжений,
ε𝑖𝑗 ,ε𝑒𝑖𝑗 ,ε𝑝𝑖𝑗 , - симметричных тензоров полных, упругих и пластических дефор
маций, 𝑢𝑖 - вектора перемещений, 𝜕𝑖 - частные производные по координатам
x, y, z декартовой ортогональной системы, λ, µ - константы Ламе, δ𝑖𝑗
- символ Кронекера, индексы 𝑖, 𝑗, 𝑘 независимо пробегают значения 1, 2, 3 и
принимается соглашение о суммировании по повторяющимся индексам.
Уравнения (1.1) - (1.3) применяются к пластине с пластическими де
формациями, однородно распределенными по ее плоскости и неоднородно
распределёнными по ее толщине:
⎛
⎞
ε𝑝 (𝑧)
0
0
⎜
⎟
ε𝑝𝑖𝑗 = ⎝ 0
ε𝑝 (𝑧)
0
⎠
𝑝
0
0
−2ε (𝑧)
(1.4)
Здесь декартова координата 𝑧 направлена вдоль нормали к плоскости
пластины, а 𝑥 и 𝑦 располагаются в этой плоскости вдоль главных осей тензора
пластических деформаций. Распределение пластических деформаций (1.4) от
ражает характер многих технологических процессов упрочняющей обработки
плоских поверхностей (дробеструйной, ультразвуковой, лазерной, кавитаци
онной, микропластической и т.д.). Полагается, что такое распределение в
свободной пластине вызывает однородное в плоскости пластины распреде
ление напряжений σ𝑖𝑗 (𝑧) и полных деформаций ε𝑖𝑗 (𝑧) Данная гипотеза будет
11
выполняться точно для неограниченной в плане пластины, а для тонких пла
стин с конечным характерным размеров в плане 𝑙, много большим толщины ℎ,
- всюду за исключением областей, примыкающих к краям. Граничные усло
вия задавались в зависимости от типа решаемой задачи и согласно расчетным
схемам изображённым на рисунках (2.5),(3.7),(4.1). Профиль пластических
деформаций имеет вид (Рисунок 1.1)
Рисунок 1.1 — Профиль пластических деформаций
Cтруктура тензора пластических деформации (1.4) имеет трансфер
сальную изотропию в координатах касательной плоскости и нормали обра
ботанной свободной поверхности. Следует заметить, что тензор ε𝑝 (1.4) ,
источник остаточных напряжений, можно трактовать как тензор деформа
ций пластины с неоднородным анизотропным полем тензора коэффициентов
линейного температурного расширения под действием однородного поля тем
пературы:
𝑓 (𝑧) = 𝑘(𝑧)∆𝑇 ,
(1.5)
где 𝑓 (𝑧) - источник остаточных напряжений, 𝑘(𝑧) - коэффициенты линей
ного термического расширения, ∆𝑇 - однородное поле температур
В расчетах источник остаточных напряжений задавался в термоупругой
постановке с однородным полем температуры ∆𝑇 и неоднородным анизотроп
ным полем тензора коэффициентов линейного температурного расширения
(1.5). Реконструкция пластических деформаций показана в работе [8]. Зада
вались следующие упругие свойства стали 65Г:𝐸 =2.15·105 МПа, ν= 0.34.
12
Для решения задач использовался программный комплекс 𝐴𝑁 𝑆𝑌 𝑆
𝐴𝑃 𝐷𝐿, который позволяет задавать анизотропную матрицу, коэффици
ентов линейного термического расширения. Использовалась формула (1.5).
Подставляя в неё известные остаточные напряжения и температуру полу
чали коэффициенты линейного термического расширения. Через команду
𝑚𝑝,𝑎𝑙𝑝𝑖 где 𝑖 − 𝑥,𝑦,𝑧, задавались полученные коэффициенты линейного
термического расширения по нормали к обрабатываемой поверхности.
13
Глава 2. Исследование стрелы прогиба полосы в зависимости от
ее ширины
2.1
Постановка задачи
Основными методами определения остаточных напряжений являются
механические и рентгеновские. Весьма перспективны для промышленного
применения электрофизические методы, при которых остаточные напряже
ния определяются по изменению электромагнитных свойств поверхностного
слоя. Для оптически активных прозрачных материалов или покрытий можно
применять поляризационно-оптические методы фотоупругости и фотопла
стичности. В работе использовался металлический образец, для которого
удобнее всего использовать механический метод. Механические методы опре
деления остаточных напряжений получили наибольшее распространение не
только из-за своей простоты, но и вследствие того, что в них используются
такие же представления о напряжениях и деформациях механики твердого
тела, как и при расчетах деталей на прочность, жесткость и устойчивость.
Это облегчает практическое применение результатов исследований для де
талей из различных материалов (металлов, керамики, пластмасс и т. д.).
Механические методы основаны на предположении, что разрез или удале
ние части детали с остаточными напряжениями эквивалентны приложению к
оставшейся детали, на вновь появившихся поверхностях, напряжений, обрат
ных остаточным. Эти обратные напряжения вызывают деформацию детали
или усилия в устройствах, препятствующих деформации. Измеряя возникшие
деформации (деформационными методами) или силы реакций (силовыми ме
тодами), можно вычислить остаточные напряжения.
Существует инженерная проблема с которой сталкиваются при произ
водстве. Во время работы с различными пластинами возникает вопрос, как
правильно подобрать ширину для различных прочностных испытаний и как
ширина влияет на прогиб, оказалась что от ширины прогиб не меняется если
его мерить по длине в центральном сечении.
Рассмотрим задачу исследование прогиба пластины, поверхностно
упрочненной, в зависимости от ее ширины. В данной работе в качестве об
14
разца рассматривалась калибровочная пластина из стали 65Г, обработанная
с одной стороны. Полевые уравнения (1.1) - (1.4),профиль пластических
деформаций (Рисунок 1.1), размеры пластины меняются по ширине, начало
координат помещено в центр пластины (Рисунок 2.1) −𝑙 ⩽ 𝑥 ⩽ 𝑙, −𝑏 ⩽ 𝑦 ⩽ 𝑏,
0 ⩽ 𝑧 ⩽ ℎ, расчетная схема (Рисунок 2.2). Граничные условия запишутся
как (2.1). Далее фиолетовым цветом будет обозначался объем с пластической
деформацией, а бирюзовым объем без каких либо деформаций.
Рисунок 2.1 — Пластина и связанная с ней система координат
Рисунок 2.2 — Расчетная схема
15
𝑢𝑥 = 𝑢𝑦 = 0
(𝑥 = 𝑦 = 𝑧 = 0)
𝑢𝑧 = 0 (𝑥 = −𝑙; 𝑦 = −𝑏; 𝑧 = 0)
𝑢𝑧 = 0
(𝑥 = −𝑙; 𝑦 = 𝑏; 𝑧 = 0)
𝑢𝑧 = 0
(𝑥 = 𝑙; 𝑦 = 𝑏; 𝑧 = 0)
𝑢𝑧 = 0
(𝑥 = 𝑙; 𝑦 = 𝑏; 𝑧 = 0)
𝑢𝑥 = 0
(𝑥 = 𝑙; 𝑧 = 𝑦 = 0)
𝑢𝑧 = 0
(𝑥 = −𝑙; 𝑧 = 𝑦 = 0)
(2.1)
16
2.2
Численное моделирование в программном комплексе ANSYS
APDL
Были выполнены расчеты в программе ANSYS APDL для упругой пла
стинки размерами 100 х 5h х 1.3мм, 100 х 10h х 1.3мм, 100 х 15h х 1.3мм,
100 х 20h х 1.3мм, 100 х 25h х 1.3мм, 100 х 50h х 1.3мм (размеры l, b и
h соответственно ) со слоем пластических деформаций(1.4). Пластические
деформации вводились в расчет с помощью введения различных коэффи
циентов линейного термического расширения. Пластины были разбиты на
гексагональные структурированные конечные элементы (Рисунок 2.3).
Рисунок 2.3 — Конечно - элементная модель пластинки
17
2.3
Результаты численного решения
Численный расчет показал, что пластины выгибаются в направлении
границы, вблизи которой заданы пластические деформации. Показана одна
из пластин (Рисунок 2.4).
Рисунок 2.4 — Перемещение 𝑢𝑦 пластины, вызванное пограничным слоем
пластических деформаций
18
По данным численного эксперимента была построена зависимость пе
ремещения от ширины пластины 𝑢𝑦 |𝑧=𝑦=0,−𝑙⩽𝑥⩽𝑙 по центральному сечению
вдоль длины (Рисунок 2.5).
�������
|�� -��|/��
�� × ��-�
�� × ��-�
�� × ��-�
�� × ��-�
�������
�
�
�
�
��
�� /��
Рисунок 2.5 — Зависимость стрелы прогиба изогнутой полосы от ее ширины
19
2.4
Вывод по главе 2
Подтверждена независимость стрелы прогиба полосы, подвергнутой од
носторонней однородной поверхностной обработке пластическим деформиро
ванием, от ее ширины. Данный факт важен при лабораторном исследовании
остаточных напряжений методом Н.Н.Давиденкова, в особенности нестан
дартных образцов, разработке нового оборудования для этих задач.
20
Глава 3. Исследование особенностей распределения остаточных
напряжений в пластине с локальной областью поверхности,
подвергнутой лазерной проковке
3.1
Постановка задачи
Рассмотрим задачу о локальном наведении остаточных напряжений
в образце определенной геометрии, приведенной в работе [20]. Некоррект
но говорить о глубине профиля остаточных напряжений, если пластина
обработана неоднородно, так как поля остаточных напряжений могут урав
новешиваться, если пластина достаточно глубокая или, наоборот, если
пластинка обработана однородно глубина может быть достаточно невелика.
Исследуем сравнение распределения остаточных напряжений по глу
бине пластики с размерами 50х20х4 мм неоднородно обработанной пласти
ческой деформацией в центре размерами 15х15 мм и пластинки 50х20х4
мм однородно обработанной пластической деформацией. Полевые уравнения
(1.1) - (1.4), профиль пластических деформаций представлен на рисунке (1.1).
Расчетная схема приведена на рисунках 3.1,3.2. Начало координат задавалось
в центре пластины (Рисунок 2.1), −𝑙 ⩽ 𝑥 ⩽ 𝑙, −𝑏 ⩽ 𝑦 ⩽ 𝑏, 0 ⩽ 𝑧 ⩽ ℎ. Гранич
ные условия одинаковы для задач и записываются как (3.1)
21
Рисунок 3.1 — Расчетные схемы для однородно обработанной пластины
Рисунок 3.2 — Расчетные схемы для неоднородно обработанной пластины
𝑢𝑥 = 𝑢𝑦 = 0
(𝑥 = 𝑧 = 0; 𝑦 = −𝑏)
𝑢𝑥 = 𝑢𝑦 = 0
(𝑥 = 𝑧 = 0; 𝑦 = 𝑏)
𝑢𝑧 = 0 (𝑥 = −𝑙; 𝑦 = −𝑏; 𝑧 = 0)
𝑢𝑧 = 0
(𝑥 = −𝑙; 𝑦 = 𝑏; 𝑧 = 0)
𝑢𝑧 = 0
(𝑥 = 𝑙; 𝑦 = −𝑏; 𝑧 = 0)
𝑢𝑧 = 0
(𝑥 = 𝑙; 𝑦 = 𝑏; 𝑧 = 0)
(3.1)
22
3.2
Численное моделирование в программном комплексе ANSYS
APDL
Пластина однородно обработанная пластической деформацией была
разбита на гексогональные структурированные конечные элементы (Рисунок
3.3). Пластина неоднородно обработанная пластической деформацией была
разбита на гексогональные конечные элементы, структурированные в месте
с неоднородной обработкой и сгущением вблизи этой области (Рисунок 3.4).
Рисунок 3.3 — Конечно - элементная модель однородно обработанной
пластины
23
Рисунок 3.4 — Конечно - элементная модель неоднородно обработанной
пластины
24
3.3
Результаты численного решения
Заданный источник пластических деформаций создает поле сжимаю
щего напряжения вблизи обеих наружных поверхностей пластинки, а также
поле растягивающего напряжения во внутренней области пластинки. Было
получено поле напряжений σ𝑥𝑥 . Поле распределения напряжений показано
на (Рисунок 3.5),(Рисунок 3.6).
Рисунок 3.5 — Распределение напряжений σ𝑥𝑥 неоднородно обработанной
пластины
Рисунок 3.6 — Распределение напряжений σ𝑥𝑥 однородно обработанной
пластины
25
По данным численного эксперимента была построена эпюра распреде
∫︀
ления напряжений σ𝑥𝑥 (Рисунок 3.7). При этом оказалось σ𝑥𝑥 (𝑧)𝑑𝑧 ≈ 0, что
не следует непосредственно из уравнений (1.1) - (1.3).
Рисунок 3.7 — Распределение напряжений σ𝑥𝑥
26
3.4
Вывод по главе 3
Показана возможность оценки глубины и максимальной величины
остаточных напряжений с использованием обработки части поверхности
контрольной пластины. Ошибка по максимальной величине остаточных на
пряжений не превышает 8%, при этом для частично обработанной пластины
удовлетворяется с хорошим приближением интегральное условие равнове
сия. Данный результат обосновывает упрощенную методику исследования
параметров профиля остаточных напряжений, наведенных различными тех
нологиями обработки.
27
Глава 4. Исследование особенностей распределения остаточных
напряжений вблизи кромки лопатки
4.1
Постановка задачи
Имеются свидетельства технологов, что при обработке лопатки поверх
ностными пластическими деформациями появляется опасный концентратор
растягивающих остаточных напряжений, который ведет к появлению пор,
что ведет к образованию трещин из нутри.
Рассмотрим задачу концентрации напряжений в поверхностно обрабо
танном затупленном клине. Начало координат задавалось в центре затуплен
ного клина на нижней грани −𝑏 ⩽ 𝑥 ⩽ 𝑏, 0 ⩽ 𝑦 ⩽ 𝑙. Профиль задавался, в
виде, изображенном на рисунке 1.1. Расчетная схема с закреплением затуп
ленного клина приведены на рисунке 4.1. Вводится параметр затупленности
клина равный отношению h/l, где h - высота затупленного клина, l - мак
симальная высота клина. Значения параметра были выбраны следующими :
0.4; 0.5; 0.6 ; 0.7; 0.8; 0.9; 0.91; 0.92; 0.93; 0.94; 0.95; 0.96; 0.97; 0.98; 0.99
Граничные условия запишутся как (4.1)
𝑢𝑥 = 𝑢𝑦 = 0 (𝑥 = 𝑦 = 0)
𝑢𝑥 = 0
(𝑦 = 0)
(4.1)
Рисунок 4.1 — Расчетная схема
28
4.2
Численное моделирование в программном комплексе ANSYS
APDL
Затупленный клин был разбит на треугольные конечные элементы
сгущающиеся кверху, так как там, предположительно, располагается зона
концентратора напряжений (Рисунок 4.2). Распределение по нормали ε𝑝 для
боковых границ ε𝑝𝑦𝑦 и сверху ε𝑝𝑥𝑥 (Рисунке4.3).
Рисунок 4.2 — Конечно - элементная модель
29
Рисунок 4.3 — Схема затупленного клина с полем пластических деформаций
30
4.3
Результаты численного решения
Заданный источник пластических деформаций создает поле сжимаю
щего напряжения по контуру, а также поле растягивающего напряжения во
внутренней области затупленного клина. Распределение компоненты напря
жения σ𝑦𝑦 показано на (Рисунке 4.4).
Рисунок 4.4 — Распределение компоненты напряжения σ𝑦𝑦
31
По данным численного эксперимента было построено распределение
от параметра затупленного клина (Рису
максимальных напряжений σ𝑚𝑎𝑥
𝑦𝑦
нок 4.5).
Рисунок 4.5 — Распределение максимальных напряжений от параметра
затупленного клина
Так же была проведена сходимость по максимальным напряжениям
σ𝑚𝑎𝑥
𝑦𝑦 . Путём изменения количества элементов записывались максимальные
напряжения до тех пор, пока напряжения не вышли на одно значение (Ри
сунок 4.6).
Рисунок 4.6 — Сходимость максимальных напряжений
32
4.4
Вывод по главе 4
Показано образование опасного концентратора растягивающих напря
жений вблизи кромки затупленного клина при его поверхностной обработке
пластическими деформациями, способного приводить к пластическому раз
рыхлению и образованию усталостных трещин в данной области.
33
Заключение
Основные результаты работы заключаются в следующем:
1. Подтверждена независимость стрелы прогиба полосы, подвергну
той односторонней однородной поверхностной обработке пласти
ческим деформированием, от ее ширины. Данный факт важен
при лабораторном исследовании остаточных напряжений методом
Н.Н.Давиденкова, в особенности нестандартных образцов, разработ
ке нового оборудования для этих задач.
2. Показана возможность оценки глубины и максимальной величины
остаточных напряжений с использованием обработки части поверх
ности контрольной пластины. Ошибка по максимальной величине
остаточных напряжений не превышает 8%, при этом для частично
обработанной пластины удовлетворяется с хорошим приближением
интегральное условие равновесия. Данный результат обосновывает
упрощенную методику исследования параметров профиля остаточ
ных напряжений, наведенных различными технологиями обработки.
3. Показано образование опасного концентратора растягивающих на
пряжений вблизи кромки затупленного клина при его поверхностной
обработке пластическими деформациями, способного приводить к
пластическому разрыхлению и образованию усталостных трещин в
данной области.
34
Список литературы
1. Биргер, И. А. Остаточные напряжения / И. А. Биргер. — М. : Машгиз,
1963. — 232 с.
2. Kopp, R. Flexible sheet forming technology by double-sided simultaneous
shot peen forming / R. Kopp, J. Schulz // CIRP Annals. –– 2002. ––
Vol. 51, no. 1. –– P. 195––198.
3. Dounde, A. A. Study of shot peening process and their effect on surface
properties: A Review / A. A. Dounde, C. Y. Seemikeri, P. R. Tanpure //
International Journal of Engineering, Business and Enterprise Applications
(IJEBEA). –– 2015. –– Vol. 2, no. 12. –– P. 104––107.
4. Simulation of shot dynamics for ultrasonic shot peening: Effects of process
parameters / J. Badreddine [et al.] // International Journal of Mechanical
Sciences. –– 2014. –– Vol. 82. –– P. 179––190.
5. Yu, L. Finite element modeling of ultrasonic surface rolling process / L. Yu,
L. Wang, D. Wang // Journal of Materials Processing Technology. ––
2011. –– Vol. 211, no. 12. –– P. 2106––2113.
6. Кузнецов, Н. Д. Технологические методы повышения надежности дета
лей машин: справочник / Н. Д. Кузнецов, В. И. Цейтлин, В. И. Волков. —
М. : Машиностроение, 1993. — 304 с.
7. Han, B. Compressive residual stress induced by water cavitation peening:
A finite element analysis / B. Han, D. Y. Ju // Materials and Design. ––
2009. –– Vol. 30, no. 8. –– P. 3325––3332.
8. Келлер И. Э., Трофимов В. Н., Владыкин А. В., Плюснин В. В., Пе
тухов Д. С., Виндокуров И. В., К вопросу о реконструкции остаточных
напряжений и деформаций пластины после дробеструйной обработки //
Вестник самарского государственного технического университета. се
рия: физико-математические науки. — 2018. — Т. 22, № 1. — С. 40—64.
9. Champaigne, J. Shot Peening Overview [Электронный ресурс] /
J. Champaigne. — 2001. — URL: https://www.shotpeener.com/library/
pdf/2001012.pdf (дата обр. 19.04.2019).
35
10. Korsunsky, A. M. On the modelling of residual stresses due to surface peen
ing using eigenstrain distributions / A. M. Korsunsky // The Journal of
Strain Analysis for Engineering Design. –– 2005. –– Vol. 40, no. 8. ––
P. 817––824.
11. Korsunsky, A. M. Residual elastic strain due to laser shock peening: Mod
elling by eigenstrain distribution / A. M. Korsunsky // The Journal of Strain
Analysis for Engineering Design. –– 2006. –– Vol. 41, no. 3. –– P. 195––204.
12. Numerical and experimental investigation on shot-peening induced defor
mation. Application to sheet metal forming / F. Cochennec [et al.] //
JCPDS-International Centre for Diffraction. –– 2009. –– P. 511––518.
13. Musinski, W. D. On the eigenstrain application of shot-peened residual
stresses within a crystal plasticity framework: Application to Ni-base super
alloy specimens / W. D. Musinski, D. L. McDowell // International Journal
of Mechanical Sciences. –– 2015. –– Vol. 100. –– P. 195––208.
14. Mechanical approach to the residual stress field induced by shot peen
ing / J. K. Li [et al.] // Materials Science and Engineering. –– 1991. ––
Vol. A147. –– P. 167––173.
15. Shen, S. An analytical model for shot-peening induced residual stresses /
S. Shen, S. N. Atluri // Computers, Materials and Continua. –– 2006. ––
Vol. 4, no. 2. –– P. 75––85.
16. Analytical modelling for residual stresses produced by shot peening /
A. S. Franchim [et al.] // Materials and Design. –– 2009. –– Vol. 30. ––
P. 1556––1560.
17. Davis, J. A study of the residual stress induced by shot peening for an
isotropic material based on Prager’s yield criterion for combined stresses /
J. Davis, M. A. Ramulu // Meccanica. –– 2015. –– Vol. 50. –– P. 1593––1604.
18. Experimental measurement and analytical determination of shot peen
ing residual stresses considering friction and real unloading behavior /
K. Sherafatnia [et al.] // Materials Science and Engineering. –– 2016. ––
Vol. A657. –– P. 309––321.
19. Simulation of shot peening: From process parameters to residual stress Fields
in a structure / D. Gallitelli [et al.] // International Journal of Mechanical
Sciences. –– 2016. –– Vol. 344, no. 4/5. –– P. 355––374.
36
20. Experimental and numerical investigation of residual stresses in laser shock
peened AA2198 / S. Keller [et al.] // Journal of Materials Processing Tech
nology. –– 2018. –– Vol. 255. –– P. 294––307.
37
Список рисунков
1.1
Профиль пластических деформаций . . . . . . . . . . . . . . . . . . 12
2.1
2.2
2.3
2.4
Пластина и связанная с ней система координат . . . . . . . .
Расчетная схема . . . . . . . . . . . . . . . . . . . . . . . . . .
Конечно - элементная модель пластинки . . . . . . . . . . . .
Перемещение 𝑢𝑦 пластины, вызванное пограничным слоем
пластических деформаций . . . . . . . . . . . . . . . . . . . .
Зависимость стрелы прогиба изогнутой полосы от ее ширины
2.5
3.1
3.2
3.3
3.4
3.5
3.6
3.7
4.1
4.2
4.3
4.4
4.5
4.6
. . . 15
. . . 15
. . . 17
. . . 18
. . . 19
Расчетные схемы для однородно обработанной пластины . . . . . .
Расчетные схемы для неоднородно обработанной пластины . . . . .
Конечно - элементная модель однородно обработанной пластины .
Конечно - элементная модель неоднородно обработанной пластины
Распределение напряжений σ𝑥𝑥 неоднородно обработанной
пластины . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Распределение напряжений σ𝑥𝑥 однородно обработанной пластины
Распределение напряжений σ𝑥𝑥 . . . . . . . . . . . . . . . . . . . . .
Расчетная схема . . . . . . . . . . . . . . . . . . . . . . . . . .
Конечно - элементная модель . . . . . . . . . . . . . . . . . . .
Схема затупленного клина с полем пластических деформаций
Распределение компоненты напряжения σ𝑦𝑦 . . . . . . . . . .
Распределение максимальных напряжений от параметра
затупленного клина . . . . . . . . . . . . . . . . . . . . . . . .
Сходимость максимальных напряжений . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
22
22
23
24
25
25
26
28
29
30
31
. . . 32
. . . 32
38
Отзывы:
Авторизуйтесь, чтобы оставить отзыв