ДИНАМИКА, ПРОЧНОСТЬ, НАДЕЖНОСТЬ |
УДК 621.73
A. Л. Воронцов, В. К. С т р а т ь е в,
B. П. Ступников
ОПРЕДЕЛЕНИЕ НАПРЯЖЕННОГО И ДЕФОРМИРОВАННОГО СОСТОЯНИЙ И УЧЕТ УПРОЧНЕНИЯ ПРИ ВНЕДРЕНИИ ПУАНСОНА В ТЕЛО БОЛЬШИХ ПОПЕРЕЧНЫХ РАЗМЕРОВ. Ч. 2
Определены напряженное и деформированное состояния материала тела больших поперечных размеров при вдавливании в него осесимметричного пуансона. Показана методика учета упрочнения деформируемого материала. Путем сопоставления с независимыми экспериментальными данными подтверждены достоверность и высокая точность полученных теоретических результатов.
E-mail: [email protected]
Ключевые слова: внедрение пуансона в тело больших поперечных размеров, напряжение, напряжение текучести, накопленная деформация, упрочнение.
Продолжим математическое исследование задачи, сформулированной авторами в работе [1], где изложено начало решения.
Используем повышающее оценку силовых параметров (т.е. согласующееся с используемым методом верхней оценки; доказательство см. в работе [2, с. 172-173]) упрощенное условие пластичности
Оу — Op = Os. (1)
Из условия (34) работы [1] следует, что
доу = до?
ду ду
Подставив равенства (1) и (2), а также формулу (33) из работы [1] в уравнения равновесия (2) и (3), приведенные в работе [1], получим
доп F (р) 2os
-ТТ1 + 2-cos у--= 0; (3)
др р р
sin у = 0. (4)
дор +
ду
Продифференцировав уравнение (3) по у, а уравнение (4) по р, придем к уравнениям
— 2 ^sin у = 0; (5)
дрду р
д2
Jp
дрдр
+
рд 2F (р) + 4 dF (р)
sin р = 0.
др2 ' ' др
Вычитая из уравнения (6) уравнение (5), получаем уравнение
д 2F (р) dF (р) F (р)'
р
+ 4-
+ 2-
sin р = 0,
др2 др р которое после умножения на р и сокращения sin ф получает вид
2 д2F (р) дF (р) Р + 4р-^ + 2F(р) = 0.
др
др
(6)
(7)
(8)
Исследуем это однородное дифференциальное уравнение второго порядка, представив его структуру как
р1(р) + Р2(р) ^ + P3fc)F (р) = 0,
др2 др гДе ф1(р) = P2, ф2(р) = 4р, фз(р) = 2. Сначала находим
(9)
-Р1(Р) = 2р- д 2Р1(Р) = 2-
др
Поскольку
др2
др2(р) д 2р1(р)
дР2(р) др
= 4.
= Рз(р^
(10)
(11)
др др2
то, следовательно, уравнение (9) является уравнением в полных дифференциалах и может быть сведено к уравнению первого порядка, имеющему вид [3, с. 145]
дР1(р)"
, ^F(р)
или
др
р
Р2(р) -
др
F (р) = C
-F (р) др
+ 2рF (р) = C,
(12)
(13)
где С — произвольная постоянная.
Разделив уравнение (13) на р2, получим линейное уравнение
^+2 ад = С.
р р р2
Это уравнение имеет структуру
дF (р)
р
+ Фl(р)F (р) = Ф2 (р),
(14)
(15)
где Ф1(р) = 2/р, Ф2(р) = С/р2.
Рис. 1. Расчетная схема вдавливания пуансона с плоским (слева) или сферическим (справа) торцом в полупространство
Решение такого уравнения, как известно, имеет вид [3, с. 36]
Ф2(Р)
F (р) = Ф(р)
dp + C1
где
Ф(р)
ф(р) = e- f *l(p)dp.
С учетом этого можно записать
_ г 2 dp — р J 2 p — p-21n p _
(16)
Ф(р) = e J Р = e
1
F (Р) = Р5
Cdp + Ci
р2
Cp + Ci p2 ■
(17)
(18)
Подставив функцию (18) в равенство (33) из работы [1], найдем касательное напряжение
Cp + C1 . -2-sin у.
ТР¥
Р2
(19)
Произвольные постоянные определяем из граничных условий (рис. 1)
Тру = Р\(73 При у = 90° и р =1; тру = —а8/\/3 при у = 90° и р = Я,
с учетом которых
C = -
R2 + 2pi _
V3(R - 1)'
(20)
С = „ + 2^1). (21)
1 ' ^3(Л - 1)
При необходимости с помощью полученной формулы (19) и уравнения (33) (с учетом равенств (29)-(31)) из работы [1] можно конкретизировать функцию f (р), получив явные поля скоростей, строго согласованные с касательными напряжениями и коэффициентом трения р1.
Подставив равенство (18) в уравнение (3), получим
дар Cp + CI 2а s
-— + 2---cos ф--
dp p3 p
+ 2 r 3 1 cos ф--s = 0, (22)
dp
откуда
ар = 2aslnp + 2Cp + Cl cosp + fi((p) + C2. (23)
P
Подставив теперь равенство (18) в уравнение (4), получим дар 2Cp + C1 .
я =--2-sin P' (24)
p p2
откуда
2Cp + Ci , ч
ар =-2-cos p + f 5 (p) + C2. (25)
p2
Сравнивая выражения (23) и (25), видим, что
f4(p)=0,
(p) = 2as ln p.
С учетом этого
(26)
а„ = 2aslnp + 2Сp + Cl cosp + C2. (27)
p2 2
Из условия пластичности (1)
av = as + 2as ln p + 2Cp + Cl cos p + C2. (28)
p2 2
Произвольную постоянную С2 находим из граничного условия
ау = qc при p = 90° и p = 1:
C2 = qc - as. (29)
Дадим небольшое пояснение. В силу того что выражение (28) зависит от координаты p, удовлетворить условие = q^; во всех точках горизонтальной границы между областями 1 и 2 невозможно. Поэтому с учетом предстоящей минимизации силы процесса в целях нахождения радиуса R, для обеспечения равенства мы выбрали точку p = 1, исключив тем самым вхождение в произвольную постоянную С2 громоздких выражений, содержащих искомый параметр R. В принципе,
использование иной точки соблюдения граничного условия, например, срединной точки p = (R + 1)/2, никаких математических затруднений не вызывает, но, не давая существенного изменения окончательных расчетных результатов, делает завершающие формулы более громоздкими.
С учетом равенства (29) из выражения (27) запишем
ap = 2аs ln p + 2Cp + Cl cos y - as + qc. (30)
P2
На контактной поверхности пуансона нормальное напряжение
ap|p=1 = (2C + Ci) cos у - as + qc, (31)
а касательное напряжение с учетом выражений (19)-(21) имеет вид
т9'-р lp=1 = №ias sin у. (32)
Сила вдавливания определяется выражением
П 2
P = nr2 I ap|p=1 | sin2y + Tp^|p=1 siny^J dy. (33)
o
Интегрируя выражение (33) с учетом формул (31), (32), а также равенств (20), (21) и формулы (17) из работы [1], находим
P = osnr2 < 1 + 0,385
R(R + 2р\) (3R + 5,2^)s R — 1 + R2 — 1 Удельная деформирующая сила определяется как
(34)
q = = о Л 1 + 0,385
nr2
R(R + 2^i) (3R + 5, 2p2)s
(35)
R - 1 R2 - 1
Для анализа деформированного состояния переносную систему координат, связанную с торцом пуансона, считаем условно неподвижной, полагая, что металл, расположенный под нижней границей очага пластической деформации, движется навстречу пуансону со скоростью v0 (показана штриховой линией на рис. 1). Влиянием угловых сдвигов на накопленную деформацию пренебрегаем. В области 2 скорости течения частиц металла задаем наиболее простыми выражениями, удовлетворяющими граничным условиям и условию несжимаемости:
Р - 1
cos у;
(36)
Vp = -vo——- cos у; R1
3р - 2 ' 2(R - 1)
Vy = —— sin у.
Подставляя выражения (36) в равенства (10)-(12) из работы [1], находим скорости деформаций:
с уо
= - -^псой
(37)
с с Уо
Ъ = е. = 2(--1)с08 р
после чего по соотношению (8) [1] находим интенсивность скоростей деформации:
6 = с0й Р (38)
Полагая, что накопленная деформация е^ в явном виде зависит только от координаты р, из известного выражения для определения накопленных деформаций по А.А.Ильюшину [2, с. 175] получаем уравнение
6 = Ур |' (39)
откуда после подстановки выражений (36) и (38) с учетом начального условия е^ = 0 при р = р0 находим
е = 1п ^. (40)
р - 1
Проинтегрировав выражение йр = УрйЬ, учитывая первую формулу из (36) и начальное условие р = р0 при 8 = 0, можно найти зависимость между начальными и текущими координатами материальной точки:
р = 1 + (ро - 1)е-^ с°8* (41)
где е — основание натурального логарифма.
В области 2, где сосредоточены пластические деформации, можно выделить две зоны с различным деформированным состоянием (см. рис. 1, справа): зону 2а, в которую постоянно вовлекаются свежие частицы металла, каждая из которых проходит через нижнюю границу радиуса Я, и зону 2б, в которой сосредоточены частицы, находившиеся в очаге пластической деформации в момент начала вдавливания. Уравнение границы между этими двумя зонами можно найти с помощью выражения (41), подставив в него р0 = Я:
в
р1 = 1 + (Я - 1)е-я-1с°8(42)
На оси симметрии, т.е. при р = 0, положение границы определяется радиусом
ге = 1 + (Я - 1)е- л-1. (43)
С учетом того что для любой частицы металла, расположенной в зоне 2а, начальной координатой являлась координата р0 = R, из формулы (40) находим выражение для определения накопленной деформации в этой зоне:
R1
ei = ln--. (44)
Р - 1
Чтобы определить накопленные деформации в зоне 2б, надо подставить в формулу (40) выражение (41):
s
ei = -5—т cos у. (45)
R - 1
На оси симметрии, т.е. при у = 0, из формулы (45) получаем
s
ei1 ei max Г. (46)
R - 1
Эпюра накопленных деформаций на оси симметрии показана на рис. 1. С достаточной точностью можно принять, что в зоне 2а накопленная деформация изменяется линейно. Тогда средняя накопленная деформация на оси симметрии с учетом формул (43) и (46) будет равна
eii(re — 1) + 0,5ei1(R — re) ( -\ s ,ж ч
ei = -) R_,1 il(-^ = 0,5^1 + e R-ij —-. (47)
Для упрощения примем аналогично работе [4], что средняя накопленная деформация на оси симметрии (47) определяет среднее значение накопленной деформации во всем очаге пластической деформации. С учетом этого можно найти среднее по очагу деформации напряжение текучести и по формуле (35) вычислить удельную силу, необходимую для вдавливания в упрочняющийся материал.
Радиус R, определяющий очаг пластической деформации, находим из условия минимума удельной деформирующей силы. При s = 0 из условия dq/dR = 0 с учетом выражения (35) получаем, что
Ro = 1 + y/T+2^i. (48)
при = 0,5 начальный радиус R0 = 2,414.
Численная минимизация при s > 0 показала, что по мере увеличения хода вдавливания s радиус R, определяющий внешний размер очага пластической деформации, увеличивается по зависимости
R =1 + /1 + 2р1 + 0,48s + ку(2 - e-s - e-5s). (49)
Если пуансон имеет рабочий торец в форме полусферы (см. рис. 1, справа), то в формулы следует подставлять значения р1, руководствуясь рекомендациями раздела 3.2 монографии [2] или работы [5]. При любой другой форме торца пуансона, вписывающейся в полусферу, в упомянутые формулы следует подставлять ¡л1 = 0,5, поскольку, как было указано ранее, в этом случае под торцом пуансона образуется
застойная зона в виде полусферы, являющаяся как бы продолжением пуансона, нивелирующим конкретный вид его торца.
Пример 1. Определить изменение удельной силы по ходу холодного выдавливания полой толстостенной детали из армко-железа с пределом текучести 230 МПа пуансоном с плоским торцом и сравнить расчетные значения с экспериментальными, приведенными в табл. 1 (соответствуют рис.46, с. 250 работы [4]; дополнительно укажем, что данные экспериментальные результаты приведены и на рис. 51, с. 441 справочника [6]).
Таблица 1
Сравнение результатов расчета удельной силы холодного выдавливания полой толстостенной детали из армко-железа с экспериментальными данными
s R ei as, МПа q, МПа q^ МПа 5, %
0,0 2,41 0,000 230 746 750 0,5
0,2 3,28 0,084 327 1130 1050 7,1
1,0 4,44 0,254 472 1882 1850 1,7
2,0 5,15 0,390 550 2429 2450 0,9
3,0 5,71 0,487 590 2817 2850 1,2
4,0 6,22 0,561 614 3128 3000 4,1
Решение. Поскольку торец пуансона плоский, то принимаем = 0,5 и = 0,1. В соответствии с разделом 3.1 монографии [2] используем коэффициент упрочнения ку = 0,95. Напряжение текучести, МПа, армко-железа определяем по аппроксимации
а, = 610 + 520е-е" - 900е-2е*,
полученной для кривой упрочнения 1, показанной на рис. 37, с. 52
[7]. Последовательно используя формулы (49), (47) и (35), находим результаты, представленные в табл. 1.
Пример 2. Определить изменение удельной силы по ходу холодного выдавливания полой толстостенной детали из стали 15 пуансоном с плоским торцом и сравнить расчетные значения с экспериментальными, приведенными в табл. 2 (соответствуют рис. 22, с. 299 справочника
[8]; дополнительно укажем, что данные экспериментальные результаты совпадают и с приведенными на рис. 54, показанном на с. 449 справочника [6]).
Решение. Поскольку торец пуансона плоский, то принимаем ¡л1 = 0,5 и р2 = 0,1. Опять используем коэффициент упрочнения ку = 0,95. Определение напряжения текучести, МПа, стали 15 проводим по аппроксимации
а, = 790 + 260е-е" - 750е-2е*,
Таблица 2
Сравнение результатов расчета удельной силы холодного выдавливания полой толстостенной детали из стали 15 с экспериментальными данными
s R ei as, МПа q, МПа q^ МПа 5, %
0,00 2,41 0,000 300 973 1000 2,7
0,35 3,65 0,124 434 1560 1500 3,9
0,75 4,20 0,210 508 1958 2000 2,1
1,50 4,82 0,329 589 2482 2500 0,7
2,25 5,29 0,417 636 2870 2800 2,4
2,60 5,49 0,452 652 3023 3000 0,8
полученной для кривой упрочнения 4, показанной на рис. 37, с. 52 из работы [7]. Последовательно используя формулы (49), (47) и (35), находим результаты, представленные в табл. 2.
Подчеркнем, что использованные нами в примерах 1 и 2 аппроксимации применялись в целях облегчения автоматизированного выполнения расчетов с помощью ЭВМ. Те же самые результаты можно получить и без данных аппроксимаций, переводя найденные накопленные деформации в относительные и выбирая соответствующие им значения напряжений текучести непосредственно с указанных кривых упрочнения. Продемонстрируем такой путь вычислений в двух следующих примерах.
Пример 3. Определить изменение удельной силы по ходу холодного выдавливания полой толстостенной детали из отожженной стали 10 пуансоном с плоским торцом и сравнить расчетные значения с экспериментальными, приведенными в табл. 3 (соответствуют рис. 51, с. 441 справочника [6]).
Таблица 3
Сравнение результатов расчета удельной силы холодного выдавливания полой толстостенной детали из отожженной стали 10 с экспериментальными
данными
s R ei e as, МПа q, МПа q3, МПа 5, %
0,4 3,74 0,136 0,127 460 1670 1600 4,2
0,8 4,25 0,219 0,197 530 2059 1950 5,3
1,2 4,60 0,286 0,249 570 2328 2250 3,4
2,0 5,15 0,390 0,323 600 2652 2600 2,0
2,8 5,60 0,470 0,375 620 2920 2900 0,7
Решение. Поскольку торец пуансона плоский, то принимаем = 0,5 и = 0,1. Опять используем коэффициент упрочнения ку = 0,95. Определение напряжения текучести отожженной стали 10
Рис. 2. Кривая упрочнения отожженных сталей 10 (а) и 20 (б)
проводим по кривой упрочнения, на рис. 2, а (соответствует кривой упрочнения 5, показанной на рис. 31, с. 62 справочника [6]). Последовательно, используя формулы (49), (47) и (35), а также формулу для условной относительной деформации (для которой была построена кривая упрочнения в справочнике)
e = 1 — e-
(50)
получаем результаты, представленные в табл. 3.
Пример 4. Определить изменение удельной силы по ходу холодного выдавливания полой толстостенной детали из отожженной стали 20 пуансоном с плоским торцом и сравнить расчетные значения с экспериментальными, приведенными в табл. 4 (соответствуют рис. 51, с. 441 справочника [6]).
Решение. Поскольку торец пуансона плоский, то принимаем = 0,5 и = 0,1. Снова используем коэффициент упрочнения ку = 0,95. Определение напряжения текучести отожженной стали 20 проводим по кривой упрочнения, приведенной на рис. 2, б (соответствуют рис. 31, с. 62 справочника [6]). Последовательно используя формулы (49), (47) и (35), а также формулу (50), находим результаты, представленные в табл. 4.
Вдавливание пуансона со сферическим торцом на начальном этапе требует специального рассмотрения (рис. 3). В этом случае форма внедренной части торца в начальный момент близка к плоской и до момента полного внедрения будет дополняться образующейся застойной зоной, в связи с чем можно использовать уже полученные формулы, сделав в них соответствующие изменения.
Таблица 4
Сравнение результатов расчета удельной силы холодного выдавливания полой толстостенной детали из отожженной стали 20 с экспериментальными
данными
s R ei e as, МПа q, МПа q:, МПа 5, %
0,4 3,74 0,136 0,127 500 1816 1750 3,6
0,8 4,25 0,219 0,197 600 2331 2150 7,8
1,2 4,60 0,286 0,249 620 2532 2400 5,2
2,0 5,15 0,390 0,323 650 2873 2800 2,5
2,2 5,26 0,412 0,338 660 2967 2900 2,3
Во-первых, нужно принять во внимание, что на данном этапе область 1 (см. рис. 1) будет отсутствовать, т.е. надо убрать соответствующие члены в формулах (35) и (49), положив в них 8 = 0.
Во-вторых, нужно принять во внимание, что на начальном этапе удельная сила д будет приложена на радиусе гн, в связи с чем удельная сила на этом этапе, приходящаяся на наибольший радиус пуансона г = 1, определяется выражением
дн = — д = гнд. (51)
г
По теореме Пифагора текущее значение ширины
Гн = л/1 - (1 - 8н)2, (52)
где 8н — глубина внедрения на начальном этапе.
В-третьих, с учетом наличия застойной зоны, в принципе, следует подставить в формулы (35) и (49) = 0,5. Однако в случае, если на поверхности самого пуансона коэффициент трения меньше 0,5, такая подстановка приведет к разрыву теоретической кривой изменения удельной силы в момент завершения внедрения торца, т.е. при 8н = 1, из-за нестыковки формул с разными коэффициентами трения: 0,5 и р1. Поэтому для обеспечения гладкой кривой теоретического измене-
Рис. 3. Расчетная схема начального этапа вдавливания пуансона со сферическим торцом в полупространство
ния силы лучше с самого начала и до конца процесса вдавливания использовать одно и то же значение
С учетом этих замечаний формулы (35) и (49) для расчета удельной силы на начальной стадии внедрения примут вид
qн
н = аау/1 - (1 - 8н)2
1 +
0,385R(R + 2^i)
R- 1
R = 1 + ^1 + 2^.
(53)
(54)
После полного внедрения сферического торца, т.е. при зн = 1, следует переходить на расчет по формулам (35) и (49).
Ввиду подтвержденной высокой точности полученных расчетных формул их можно рекомендовать для применения в расчетах технологических операций вдавливания осесимметричного пуансона в заготовки больших поперечных размеров.
СПИСОК ЛИТЕРАТУРЫ
¡.Воронцов А. Л., Стратьев В. К., Ступников В. П. О пробе Бри-нелля и внедрении пуансона в тело больших поперечных размеров // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2012. - № 2. - С. 12-21.
2. В о р о н ц о в А. Л. Технологические задачи теории пластичности. Т. 1. - М.: Машиностроение, 2006. - 474 с.
3. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. -М.: Наука, 1976. - 576 с.
4. С м и р н о в - А л я е в Г. А. Сопротивление материалов пластическому деформированию. - Л.: Машиностроение, 1978. - 368 с.
5. Воронцов А. Л. Выбор коэффициентов трения в расчетах процессов выдавливания // Производство проката. - 2003. - № 4. - С. 2-5.
6. Холодная объемная штамповка: Справочник / Под ред. Г.А. Навроцкого. -М.: Машиностроение, 1973. -496 с.
7. К р о х а В. А. Кривые упрочнения металлов при холодной деформации. - М.: Машиностроение, 1968. - 131 с.
8. К о в к а и штамповка: Справочник / Под ред. Е.И. Семенова. Т. 3. Холодная объемная штамповка / Под ред. Г.А. Навроцкого. - М.: Машиностроение, 1987. -384 с.
Статья поступила в редакцию 18.10.2011
Андрей Львович Воронцов родился в 1955 г., окончил МВТУ им. Н.Э. Баумана в 1978 г. Заслуженный деятель науки РФ, д-р техн. наук, профессор кафедры "Технологии обработки материалов" МГТУ им. Н.Э. Баумана. Автор 420 научных работ в области прикладной теории пластичности, обработки металлов давлением и обработки резанием.
A.L. Vorontsov (b. 1955) graduated from the Bauman Moscow Higher Technical School in 1978. D. Sc. (Eng.), professor of "Technologies of Treatment of Materials" department of the Bauman Moscow State Technical University, Honored Scientist of the Russian Federation. Author of 420 publications in the field of applied theory of plasticity, noncutting shaping of metals and treatment by cutting.
Владимир Петрович Ступников родился в 1954 г., окончил МВТУ им. Н.Э. Баумана в 1980 г. Д-р техн. наук, профессор, заведующий кафедрой "Технологии обработки материалов" МГТУ им. Н.Э. Баумана. Автор 40 научных работ в области обработки металлов давлением.
V.P. Stupnikov (b. 1954) graduated from the Bauman Moscow Higher Technical School in 1980. D. Sc. (Eng.), professor, head of "Technologies of Treatment of Materials" department of the Bauman Moscow State Technical University. Author of 40 publications in the field of noncutting shaping of metals.
Валерий Константинович Стратьев родился в 1951г., окончил МВТУ им. Н.Э. Баумана в 1974 г. Канд. техн. наук, доцент кафедры "Технологии обработки материалов" МГТУ им. Н.Э. Баумана. Автор 50 научных работ в области обработки металлов давлением.
V.K. Strat'ev (b. 1951) graduated from the Bauman Moscow Higher Technical School in 1974. Ph. D. (Eng.), assoc. professor of "Technologies of Treatment of Materials" department of the Bauman Moscow State Technical University. Author of 50 publications in the field of noncutting shaping of metals.