УДК 621.81:539.4
О долговечности многослойных матриц для штамповки осесимметричных деталей
Ю.Р. Фаизов, Р.Ш. Валиев, Х.Ш. Газизов
Уфимский государственный авиационный технический университет, 450000, Республика Башкортостан, Уфа, К. Маркса ул., д. 12.
The lifetime of multilayer dies for industrial forging of axisymmetrical parts
Yu.R. Faizov, R.Sh. Valiev, Kh.Sh. Gazizov
Ufa State Aviation Technical University, Karla Marksa str., 12, 450000, Ufa, Bashkortostan, Russian Federation. e-mail: [email protected], [email protected], [email protected]
П Повышение долговечности инструмента для холодной объемной штамповки деталей сложного профиля всегда является актуальной задачей, определяющей его надежность, уровень затрат на изготовление и время эксплуатации. Рассмотрены вопросы напряженно-деформированного состояния многослойных конструкций матриц для штамповки осесимметричных деталей сложной формы. Для уменьшения переменной части максимальных напряжений предложено натяги между слоями сделать переменными по высоте инструмента. Значения натягов определяются известными методами оптимального проектирования. Показано, что предложенный подход позволяет уменьшить долю переменной части напряжений, что существенно повышает долговечность инструмента.
Ключевые слова: долговечность инструмента, надежность, напряженно-деформированное состояние, многослойные конструкции.
D Increasing the lifetime of dies for cold forging of complex parts is always an urgent problem as it determines their reliability, cost, and operating time. The problems of stress-strain state of multilayer dies intended for industrial forging of complex axisymmetrical parts are considered. To reduce variations of maximum stresses, the tension between layers are suggested to make variable along the die height. The tension values are determined by the well-known optimal design methods. It is shown that the proposed approach can reduce the variation of stresses, which significantly improves the die lifetime. It is shown that the height adjustable tension significantly increases not only its carrying capacity, but also durability. A similar problem for cylinders with a uniform thickness under constant internal pressure is referred to in the scientific literature as the academician Gadolin problem.
Keywords: lifetime, reliability, stress-strain state, optimal design, multilayer structure.
В конструкциях матриц, например, для штамповки шаровых пальцев (рис. 1) ввиду сложной формы высаживаемых деталей внутренняя (рабочая) поверхность инструмента подвержена воздействию переменного (по высоте ручья) внутреннего давления. Фрагмент эпюры единичного внутреннего давления на инструмент, полученный в результате численно-
го моделирования процесса штамповки шарового пальца в пакете DEFORM-2D приведен на рис. 2 [1].
Цель работы — создание методики проектирования многослойных матриц для штамповки осесимметричных деталей сложного профиля и их оптимизация за счет переменного по высоте поверхностей контакта слоев натяга.
Рис. 1. Разрез матрицы
К11 К
К 21 К
К (0
К 31
К (0
К 32
К (') К13 ц ю' Х(г)
К (') К 23 < ц 2 . = < 0
К (') К 33 ц ? Х 2г)
II ., п.
(1)
Здесь К;у — блоки (прямоугольные, кроме диагональных) матрицы жесткости слоя г; ц(г) — вектор (матрица-столбец) перемещений узлов на внутренней поверхности слоя г +1 (см. рис. 1); ц 2° — вектор перемещений оставшихся (внутренних) узлов; ц 3 — вектор перемещений узлов на наружной поверхности слоя г; Х(г) — вектор узловых сил на внутренней поверхности слоя; X 2г) — вектор узловых сил на
наружной поверхности слоя; 0 — нулевой вектор; п — число слоев.
Если из уравнений (1) исключить неизвестные ц2г), то получим уравнения суперэлементной формулировки [3].
Учитывая, что
Х(г) =-Х2-1), г = 2, 3, ..., п-1;
Х 2П) = 0 Х(0) = ршахр
-(0)_
(2)
(действие равно противодействию; наружная поверхность матрицы свободна от нагрузок; р — вектор узловых сил от «единичного» (см. рис. 2) внутреннего давления),
,(>)-„(> -1)
ч
Ц3
' + А
г -1>
г = 2, 3, ..., п,
(3)
где А г — вектор радиальных натягов в узлах на поверхности контакта слоев г и г - 1, уравнения (1) можно объединить в одну систему стандартного для МКЭ вида:
К • ц = И. (4)
Здесь ц — вектор неизвестных узловых перемещений,
т и (1)т и (1)т и (1)т и (2)т и (2)т
ц = и
2
и 3 и 2 и 3 ... и 2
(п)т и (п)т
3
;(5)
Рис. 2. Фрагмент эпюры внутреннего давления
В работах [2, 3] предложена методика исследования напряженного состояния многослойных матриц для штамповки осесимметричных деталей методом конечных элементов (МКЭ), причем натяги между слоями полагаются постоянными по высоте инструмента. Суть этой методики заключается в следующем: матричное уравнение равновесия г-го слоя многослойной матрицы в блочном виде записывается в виде
К — глобальная матрица жесткости многослойной конструкции, для трехслойного случая имеет вид
К =
К(1) К11 К (1) К12 К(1) К13 0 0 0 0
К(1) К21 К (1) К 22 К23) 0 0 0 0
К(1) К31 К (1) К 32 к33+К11 К(2) К12 К(2) К13 0 0
0 0 К(2) К21 К(2) К22 К(2) К 23 0 0
0 0 К(2) К31 К(2) К32 К(2) + К(3) К33 + К11 К(3) К12 К(3) К13
0 0 0 0 К(3) К21 К(3) К22 К(3) К23
0 0 0 0 К(3) К31 К(3) К32 К(3) К33
(6)
(вид матрицы для двух и произвольного числа слоев — очевиден); И — модифицированный вектор узловых сил,
Ит = рт||I 0 0 0 0...0 0 0 ||ршах -
- А? || 0 0 К(2)т К (2)т К (1)т... 0 0 0 11--|| 0 0 0 0 К13)т ...0 0 0 11-
-ап-2||0 0 0 0 0...К(1-1)т0 011-
-Ап^0 0 0 0 0...К{?)т К21)т К31)т ||, (7) где I — единичная матрица.
Если натяги по высоте инструмента постоянны, то
Ак =Ак Jk, к =1, 2, ..., п-1. (8)
Здесь Д^ — числовое значение натяга.
Если перемещения узлов в векторах u перечислены в следующем порядке: радиальное и осевое, тогда они имеют размерность 2М(к), где М(к) — число двойных узлов на стыке к,
1к =| 1 0 1 0 1 0 ... 1 0||, к = 1,2,...,п-1; (9)
Сказанное выше остается в силе и в случае, когда вектор радиальных натягов в узлах \к меняется от точки к точке на контактирующих поверхностях. Наиболее технологичным вариантом изменения натягов по высоте поверхностей контакта является кусочно-линейный закон (рис. 3), особенно если внутренний слой составлен из нескольких вставок. В этом случае закон изменения радиальных натягов по высоте инструмента (на стыке к) может быть представлен в следующем виде: пк
Дк(г) = £ Д ыХт(г). (10)
т=1
Здесь Ыт (г) — интерполирующие функции, такие же, что используются в МКЭ с линейными элементами (см. рис. 3), причем
Пк
X Nт (г) = 1 (11а)
т=1
для любого г [4]; Пк — число, на единицу большее, чем число вставок в охватываемом слое к (на рис. 3 три вставки).
Анализ структуры глобальной матрицы жесткости (6) показывает, что ее можно построить обычным для МКЭ способом, как для конструкции из неоднородных материалов. При этом не нужно специально перенумеровывать узлы, что подразумевалось при выводе формул (4) и (6). Модификации подлежит лишь
вектор узловых сил (7). В случае переменных по высоте инструмента натягов вид векторов Т несколько изменяется. Ненулевые элементы, т. е. соответствующие радиальным направлениям элементов на границах соседних вставок, по-прежнему равны 1, уменьшаются до нуля к другой границе (в соответствии с рис. 3) и могут быть представлены в следующем виде: пк
X тт=т. (иб)
т=1
После решения системы линейных алгебраических уравнений (4) неизвестные перемещения узлов на охватывающих поверхностях слоев могут быть вычислены по формуле (3), а затем, обычным способом МКЭ — деформации и напряжения в точках интегрирования конечных элементов и в узлах [3]. В частности, компоненты тензора напряжений в любой точке могут быть записаны в виде
п-1 пк кт
,г) = X X Чт(г,г)Дкт +<з?)(г,г)ртах, (12)
к =1 т=1
где г,) = 1, 2, 3 — соответственно радиальное, осевое и окружное направления; при осесим-метричной деформации ст13 = а23 = 0.
В частном случае постоянных натягов между слоями в (12) следует положить Дкт =Дк, т = 1,2,..., пк.
Однако изложенный выше алгоритм не предусматривает возможность проскальзывания контактирующих поверхностей. В реальных условиях, применительно к многослойным матрицам, где поверхности соприкосновения слоев цилиндрические или близки к таковым, как показывают вычисления, касательные напряжения на поверхностях контакта не превышают сил кулоновского трения, что не является существенным недостатком.
Основная цель создания натягов между слоями толстостенного цилиндра — уменьшение
растягивающих окружных напряжений о33 от внутреннего давления и натягов. При численном решении задачи такая цель математически может быть сформулирована как определение минимума функционала
N I к \2
/ = X ( а33 шах ) к=Г '
>Ш1И,
(13)
где N — число составляющих матрицу элементов; о3з шах — максимальное окружное напряжение на внутренней поверхности вставки к.
Уравнение (13) представляет собой основное соотношение метода наименьших квадратов [5].
Как следует из второго слагаемого (12), функционал (13) зависит от максимального значения внутреннего давления ршах (при заданном законе изменения его по длине цилиндра) и дискретного множества натягов на стыках вставок А гк. Интерес представляет только такое множество переменных проектирования ршах и А гк, которое не только обеспечивает минимум функционалу (13), но одновременно удовлетворяет условиям прочности слоев:
^экв к < ОЦш к , к - 1 2, ... , Л.
(14)
Здесь оэкв к — эквивалентное напряжение для материала к-го слоя конструкции; оцш к — предельное напряжение для материала к-го слоя конструкции. Если решается задача о несущей способности составной конструкции р11ш ^ шах(ршах), в качестве стцш к следует принять пределы текучести материала; если речь идет о выносливости, то — пределы выносливости (сложное нагружение, [6]) и в выражениях (14) знаки неравенства следует заменить знаками равенства.
При наличии ограничений (14) рассчитывается минимум функционала
м
N
У = X (°к3шах)2 + X X к (^экв к - СТ11ш к ) ^ш1п, (15)
к=1 к=1
где X к — множители Лагранжа [5].
Условия минимума У в нашем случае имеют вид
Э У
= 0;
Э У
Э рцш ЭАгк ЭУ
= 0, к = 1,2, ..., п;
(16)
Э X
= 0, г = 1,2, ..., М.
Замкнутая система из
м
N + X пг +1
г=1
уравнений (16) относительно того же числа не-
известных р11ш, Агк, X г является нелинейной, поскольку эквивалентные напряжения оэкв к, вычисленные через компоненты тензора напряжений в какой-либо системе координат (в рассматриваемом случае цилиндрическая), содержат их квадраты [6]. Только в частном случае, когда на поверхностях контакта слоев можно пренебречь касательными напряжениями и за неравенство (14) принять условие текучести Треска — Сен-Венана [7], система уравнений (16) относительно переменных проектирования будет линейной. Если многослойный цилиндр нагружен только внутренним давлением, такое предположение для первого (внутреннего) слоя будет точным. Как показали многочисленные численные эксперименты, погрешность такого упрощения дает несколько заниженные значения эквивалентных напряжений, не превышающие разность их величин, полученных за счет использования вместо условия Треска — Сен-Венана более точного, как принято считать, условия текучести Мизеса [7, 8].
Поскольку главный интерес представляет выносливость матрицы, задачу можно упростить, заменяя неравенства (1) равенствами, так как максимум внутреннего давления имеет место в случае, когда эквивалентные напряжения достигают предельных значений, т. е. предела выносливости. Остается открытым вопрос об опасных точках внутренних контуров слоев. Пока не определены натяги между слоями, их невозможно указать. Для решения этой задачи предлагается следующий итерационный алгоритм:
1) принять в первом приближении (начальный проект) натяги постоянными по высоте контактирующих поверхностей;
2) изложенным выше методом решить задачу о напряженно-деформированном состоянии составной многослойной конструкции;
3) установить положение опасных точек на внутренних контурах слоев;
4) определить оптимальные значения натягов и предельное давление на инструмент;
5) при найденных значениях этих величин вычислить эквивалентные напряжения во внутренних контурах инструмента и, если они не превышают допускаемых перейти к п. 7;
6) вернуться к п. 4 с новыми данными о давлении и положении опасных точек;
7) конец.
Следует отметить, что найденное значение предельного давления зависит не только от натягов, но и от радиусов посадочных поверхностей. Для определения несущей способности
необходим еще один итерационный процесс — поиск этих размеров. Используем для этого градиентный метод [9].
Приведенная конечно-элементная методика определения максимального внутреннего давления piim и узловых значений (на границах вставок) натягов Ark была реализована в виде пакета программ на языке Fortran 90.
Результаты решения тестовой задачи о предельном внутреннем давлении и постоянном натяге двухслойной трубы (r1 — внутренний радиус; r3 — наружный радиус; r2 = (r1r3)1/2 — оптимальное значение радиуса посадочной поверхности, обеспечивающий при радиальном натяге А наибольшее значение внутреннего давления pum) приведены ниже:
Метод plim/tflim ЕА/(«lim »2)
Аналитическое решение [10]......... 0,667 0,303
МКЭ, постоянный натяг................. 0,682 0,312
Как следует из приведенных результатов, наибольшая погрешность при вычислении предельного внутреннего давления составляет 2,2 %, а по натягу — 3 % (уменьшается при измельчении сетки конечных элементов).
По предложенному алгоритму была рассчитана несущая способность реальной конструкции двухслойной матрицы для штамповки шаровых пальцев (см. рис. 1) с собранными с постоянным натягом слоями. Закон изменения внутреннего давления был принят таким же, что и в работе [2]. Эпюра относительных напряжений Мизеса (в числителе в долях от наибольшего из предельных напряжений верхней вставки оцт 1 = 2 300 МПа) верхней части внутреннего контура инструмента представлена
0,90/0,90
0,76/0,79
на рис. 4. Максимальное относительное значение напряжений Мизеса, равное 0,90, имеет место во внутреннем контуре переходной области (поверхность тора, с центром О3, см. рис. 2), где эквивалентные напряжения, вычисленные по гипотезе Треска — Сен-Венана, достигают максимума (относительное значение равно 1,0), что подтверждает правомерность ранее принятого упрощения нелинейной системы (16). На рис. 4 в знаменателе приведены численные результаты расчетов при переменных по высоте инструмента натягах. Как следует из этих данных, использование переменного натяга приводит к некоторому выравниванию напряжений Мизе-са по высоте инструмента в наиболее нагруженном участке, формирующей сферическую часть заготовки.
Полученные результаты подтверждаются и характером разрушения верхней вставки (на рис. 5): усталостная трещина образуется во внутреннем контуре переходной области, именно там, где эквивалентные напряжения достигают максимума.
Напряжения в инструменте складываются из двух частей: постоянная часть, обусловленная натягами, и переменная, вызванная внутренним давлением. Значение последних для инструмента многоразового пользования должны быть таковы, чтобы в сумме с первыми в опасной точке оно не превышало временного предела выносливости материала инструмента.
Числовые значения напряжений Мизеса, МПа, в опасной точке для двух значений внутреннего давления (рабочее значение максимального давления принято равным 0,4рцт) приведены ниже:
Максимальное Постоянный Переменный
давление натяг натяг
1 000 .................................................. 278 904
0 ......................................................... 967 341
0,75/0,54
Рис. 4. Эпюра относительных напряжений Мизеса во внутреннем контуре матрицы
Рис. 5. Диаметральный разрез верхней вставки после разрушения
При практически равных значениях средних напряжений аш = 622 МПа амплитуды переменных напряжений составили: при постоянном натяге ая = 344 МПа, при переменных — а а = 282 МПа, т. е. уменьшились примерно на 18 %. Даже такое, на первый взгляд, незначительное уменьшение амплитуды переменных напряжений должно привести к заметному повышению долговечности матрицы. Как известно, при линейном напряженном состоянии показатель степени кривой усталости имеет достаточно большое значение, например, для сталей он составляет от 12.15. Поэтому уменьшение амплитуды напряжений в 1,18 раза должно привести к увеличению долговечности в 1,1812 раза. Однако случай малоцикловой усталости (высокий уровень напряжений) при объемном напряженном состоянии к настоящему времени все еще мало изучен. Очевидно и здесь следует ожидать повышения долговечности, может быть и в не такой степени.
В ОАО «БелЗАН» (г. Белебей) были изготовлены и проведены опытно-промышленные испытания модернизированной конструкции (с переменными натягами). Они показали, что применение переменных по высоте инструмента натягов приводит к заметному (до 50 %) увеличению долговечности матрицы.
Выводы
1. Показано, что для многослойного инструмента для штамповки осесимметричных деталей сложного профиля за счет переменного по высоте штампа натяга удается заметно уменьшить амплитуды переменных составляющих напряжений.
2. Определены числовые значения оптимальных натягов.
3. Опытно-промышленными испытаниями установлено, что применение предложенного алгоритма позволяет существенно повысить долговечность инструмента.
Литература
[1] Жернаков В.С., Газизов Х.Ш., Валиев Р.Ш., Ермоленко А.Н. Контактное взаимодейст-
вие инструмента и заготовки при штамповке круглых деталей. Вестник УГАТУ, 2006, № 1(14), с. 192-195.
[2] Овчинников С.В., Валиев Р.Ш., Газизов Х.Ш., Фаизов Ю.Р. О напряженном состоянии
составных цилиндров переменной жесткости. Вестник УГАТУ, 2011, т. 15, № 4(44), с. 95-97.
[3] Овчинников С.В., Валиев Р.Ш., Газизов Х.Ш., Фаизов Ю.Р. Напряженное состояние со-
ставных матриц для штамповки сложнопрофильных осесимметричных деталей. Современные методы и технологии создания и обработки материалов. Тр. VIIМеждунар. науч.-техн. конф. Кн. 3. Минск, Изд-во ФТИ НАН Беларуси, 2012, с. 81-86.
[4] Присекин В.Л., Расторгуев Г.И. Основы метода конечных элементов в механике дефор-
мируемых тел. Новосибирск, Изд-во НГТУ, 2010. 238 с.
[5] Губанов В.С. Обобщенный метод наименьших квадратов. Теория и применение в аст-
рометрии. Санкт-Петербург, Наука, 1997. 319 с.
[6] Пантелеев А.В. Вариационное исчисление в примерах и задачах. Москва, Изд-во МАИ,
2000. 230 с.
[7] Горшков А.Г., Старовойтов Э.И., Тарлаковский Д.В. Теория упругости и пластично-
сти. Москва, Изд-во ФИЗМАТЛИТ, 2003. 416 с.
[8] Коллинз Дж. Повреждение материалов в конструкциях. Москва, Мир, 1984. 624 с.
[9] Хог Э., Арора Я. Прикладное оптимальное проектирование. Москва, Мир, 1983. 479 с.
[10] Феодосьев В.И. Сопротивление материалов. Москва, Изд-во МГТУ, 2000. 592 с.
References
[1] Zhernakov V.S., Gazizov Kh.Sh., Valiev R.Sh., Ermolenko A.N. Kontaktnoe vzaimodeistvie instrumenta i zagotovki pri shtampovke kruglykh detalei [Numerical simulation of contact relation of the tool and of the blank while heading of rotational part]. Vestnik UGATU [Herald USATU]. 2006, no. 1(14), pp. 192-195.
[2] Ovchinnikov S.V., Valiev R.Sh., Gazizov Kh.Sh., Faizov Iu.R. O napriazhennom sostoianii
sostavnykh tsilindrov peremennoi zhestkosti [The stressed state of a composite cylinder with variable rigidity]. Vestnik UGATU [Herald USATU]. 2011, vol. 15, no. 4(44), pp. 95-97.
[3] Ovchinnikov S.V., Valiev R.Sh., Gazizov Kh.Sh., Faizov Iu.R. Napriazhennoe sostoianie
sostavnykh matrits dlia shtampovki slozhnoprofil'nykh osesimmetrichnykh detalei [Stress state of composite matrices for punching slozhnoprofilnyh axisymmetric parts]. Sovremen-nye metody i tekhnologii sozdaniia i obrabotki materialov. Trudy 7 mezhdunarodnoi nauch-no-tekhnicheskoi konferentsii [Proceedings of the 7th International Scientific and Technical Conference]. Book 3, Minsk, FTI NAN Belarusi publ., 2012, pp. 81-86.
[4] Prisekin V.L., Rastorguev G.I. Osnovy metoda konechnykh elementov v mekhanike deform-
iruemykh tel [Fundamentals of the finite element method in the mechanics of deformable bodies]. Novosibirsk, NSTU publ., 2010. 238 p.
[5] Gubanov V.S. Obobshchennyi metod naimen'shikh kvadratov. Teoriia i primenenie v astro-
metrii [Generalized least-squares method. Theory and application in astrometry]. St. Petersburg, Nauka publ., 1997. 319 p.
[6] Panteleev A.V. Variatsionnoe ischislenie v primerakh i zadachakh [The calculus of variations
in the examples and problems]. Moscow, MAI publ., 2000. 230 p.
[7] Gorshkov A.G., Starovoitov E.I., Tarlakovskii D.V. Teoriia uprugosti i plastichnosti [Theory
of elasticity and plasticity]. Moscow, FIZMATLIT publ., 2003. 416 p.
[8] Kollinz Dzh. Povrezhdenie materialov v konstruktsiiakh [Damage to materials in the con-
struction]. Moscow, Mir publ., 1984. 624 p.
[9] Khog E., Arora Ia. Prikladnoe optimal'noe proektirovanie [Applied optimal design]. Moscow,
Mir publ., 1983. 479 p.
[10] Feodos'ev V.I. Soprotivlenie materialov [Strength of materials]. Moscow, Bauman Press, 2000. 592 p.
Статья поступила в редакцию 02.06.2014
Информация об авторах
ФАИЗОВ Юрий Рустемович (Уфа) — аспирант кафедры «Сопротивление материалов». Уфимский государственный авиационный технический университет (450000, Республика Башкортостан, Уфа, ул. К. Маркса, д. 12, e-mail: [email protected]).
ВАЛИЕВ Рафаил Шамилевич (Уфа) — кандидат технических наук, доцент кафедры «Сопротивление материалов». Уфимский государственный авиационный технический университет (450000, Республика Башкортостан, Уфа, ул. К. Маркса, д. 12, e-mail: [email protected]).
ГАЗИЗОВ Хатиб Шарфзянович (Уфа) — доктор технических наук, профессор кафедры «Сопротивление материалов». Уфимский государственный авиационный технический университет (450000, Республика Башкортостан, Уфа, ул. К. Маркса, д. 12, e-mail: [email protected]).
Information about the authors
FAIZOV Yuriy Rustemovich (Ufa) — Post-Graduate of «Strength of Materials» Department. Ufa State Aviation Technical University (USATU, Karla Marksa str., 12, 450000, Ufa, Bashkortostan, Russian Federation, e-mail: [email protected]).
VALIEV Rafail Shamilevich (Ufa) — Cand. Sc. (Eng.), Associate Professor of «Strength of Materials» Department. Ufa State Aviation Technical University (USATU, Karla Marksa str., 12, 450000, Ufa, Bashkortostan, Russian Federation, e-mail: [email protected]).
GAZIZOV Khatib Sharfzyanovich (Ufa) — Dr. Sc. (Eng.), Professor of «Strength of Materials» Department. Ufa State Aviation Technical University (USATU, Karla Marksa str., 12, 450000, Ufa, Bashkortostan, Russian Federation, e-mail: [email protected])