CC BY
11
УДК 621.981.12 В01: 10.18287/2541-7533-2020-19-3-73-84
АНАЛИТИЧЕСКАЯ МОДЕЛЬ ОПРЕДЕЛЕНИЯ ПРИПУСКОВ НА СТЕСНЁННЫЙ ИЗГИБ ЭЛАСТОМЕРОМ ДЕТАЛЕЙ ЛЕТАТЕЛЬНЫХ АППАРАТОВ
© 2020
доктор технических наук, профессор, профессор кафедры производства летательных аппаратов и управления качеством в машиностроении; Самарский национальный исследовательский университет имени академика С.П. Королёва; moiseevvk@mail.ru
кандидат технических наук, доцент кафедры производства летательных аппаратов и управления качеством в машиностроении; Самарский национальный исследовательский университет имени академика С.П. Королёва; pla.gromova@ya.ru
кандидат технических наук, доцент кафедры производства летательных аппаратов и управления качеством в машиностроении; Самарский национальный исследовательский университет имени академика С.П. Королёва; oleg.lomovskoi@ya.ru
ведущий инженер;
Ульяновский филиал КБ «ПАО Туполев», г. Ульяновск; mishaufkb@ya.ru
кандидат технических наук, доцент кафедры производства летательных аппаратов и управления качеством в машиностроении; Самарский национальный исследовательский университет имени академика С.П. Королёва; anplotnikov@ya.ru
кандидат технических наук, доцент кафедры производства летательных аппаратов и управления качеством в машиностроении; Самарский национальный исследовательский университет имени академика С.П. Королёва; aa.sharov@ya.ru
Представлен метод стеснённого изгиба листовых деталей летательных аппаратов с криволинейными бортами эластичной средой. Описан процесс проведения стеснённого изгиба. Изготовление детали осуществляется за два перехода. В результате первого перехода получается деталь с утонением в радиусной части. Второй переход проводят на формблоке сниженной высоты, в результате полученная волна избыточного материала деформируется по радиусу гибочной оправки. В готовой детали наблюдается увеличение толщины заготовки в зоне радиуса гиба. Волна избыточного материала образуется за счёт специального припуска, величина которого является определяющей для достижения увеличения толщины. Аналитическая модель определения припуска базируется на общепринятых допущениях и принципах расчёта технологических параметров. Для определения припуска предлагается аппроксимация формы избыточной волны тангенциальной функцией. Это позволяет получить аналитическую зависимость, связывающую геометрию детали и утолщение стенки с величиной припуска. Величина припуска ограничивается некоторым диапазоном, нижняя граница которого определена из условия начала пластических деформаций, верхняя - из условия возможности потери устойчивости, приводящей к неисправимому браку. Результаты проиллюстрированы графиками зависимости величины минимального и максимального припуска высоты борта от соотношений толщины стенки заготовки, радиуса гиба борта и радиуса борта в плане. Графические и аналитические зависимости представлены в безразмерном виде, где геометрические параметры отнесены к радиусу гибочного инструмента, то есть к радиусу гиба борта детали по внутренней поверхности.
Технология производства летательных аппаратов; штамповка эластомером; стеснённый изгиб; деформации; припуск
Цитирование: Моисеев В.К., Громова Е.Г., Ломовской О.В., Мантусов М.Н., Плотников А.Н., Шаров А.А. Аналитическая модель определения припусков на стеснённый изгиб эластомером деталей летательных аппаратов // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2020. Т. 19, № 3. С. 73-84. DOI: 10.18287/2541-7533-2020-19-3-73-84
В. К. Моисеев
Е. Г. Громова
О. В. Ломовской
М. Н. Мантусов А. Н. Плотников
А. А. Шаров
Введение
В конструкции летательных аппаратов имеется большое множество тонкостенных деталей, изготавливаемых из листовых полуфабрикатов гибкой. При этом борт детали может быть прямолинейным либо криволинейным в плане. Гибка бортов таких деталей осуществляется преимущественно эластомерами и сопровождается нежелательным утонением стенки на радиусе гиба. Утонение регламентируется соображениями равно-прочности и оговаривается в технических условиях на конструкцию детали. Для выполнения задаваемых конструктором условий, а именно для уменьшения величины утонения, предложен технологический процесс штамповки деталей эластичным инструментом со стеснённым изгибом, защищённый патентом, а также разработаны различные схемы его реализации. Результаты проведённых исследований позволяют рассчитывать параметры и проектировать техпроцесс и технологическое оснащение для стеснённого изгиба прямолинейных бортов, тогда как по штамповке криволинейных бортов проведены лишь экспериментальные исследования и проведена единичная промышленная апробация [1-4].
Общая постановка задачи
Схема технологического процесса стеснённого изгиба при штамповке эластомером представлена на рис. 1. На нём условно изображена последовательность реализации стеснённого изгиба эластомером деталей, имеющих угол гиба 90° . Стеснённый изгиб по этой схеме производится за два перехода. Вначале (рис. 1, а) заготовка 1 изгибается по гибочной оправке 2 под действием давления эластичной среды 3, заключённой в контейнер (условно не показан). Затем под торец изогнутого борта заготовки 4 (рис. 1, б) устанавливается опорная пластина 5 для создания определённого припуска борта АН по его высоте между верхней плоскостью гибочной оправки и стенкой заготовки (рис. 1, б, слева). Для предотвращения потери устойчивости изогнутого борта заготовки 4 на опорную пластину 5 ставятся эластичные подпоры 6. Под действием давления эластичной среды стенка заготовки вначале прогибается в центральной части, а затем образуется волна избыточного материала в зоне скругления гибочной оправки с радиусом г0 (рис. 1, б, справа). При дальнейшем увеличении давления эластичной среды волна избыточного материала деформируется по радиусу гибочной оправки (рис. 1, в). В результате толщина заготовки в зоне радиуса гиба увеличивается, а пру-жинение борта уменьшается за счёт изменения схемы деформирования заготовки.
Рис. 1. Схема технологического процесса стеснённого изгиба при штамповке эластомером: а - предварительная гибка; б - начало операции окончательной гибки; в - завершение операции окончательной гибки
Кроме того, появляется возможность уменьшения радиуса изгиба, что также является актуальной задачей, так как при этом повышается жёсткость деталей и снижается их масса в результате уменьшения необходимой для крепления детали высоты борта, например, в клёпаных конструкциях сборочных единиц самолётов.
Технологическим параметром, обеспечивающим реализацию стеснённого изгиба, является припуск высоты борта полуфабриката, величина которого определяет возможность использования такой технологии в конкретном случае и толщину стенки в зоне радиуса гиба. Для назначения размера припуска требуется создание аналитической модели, позволяющей провести расчёт оптимального диапазона возможных припусков для последующего окончательного его выбора при проектировании конкретного технологического процесса.
Основная часть
В случае криволинейного в плане борта детали полуфабрикат перед последней стадией формовки имеет вид обечайки с донцем с торовым скруглением на стенке в зоне кромки оправки-формблока (рис. 2). Для расчёта необходимого припуска на втором переходе вышеописанного процесса штамповки можно выделить две стадии:
1) посадка заготовки на поверхность формблока с образованием волны избыточного материала над кромкой последнего;
2) осадка волны на кромку формблока.
В холодной листовой штамповке принято рассчитывать размеры заготовок, исходя из равенства площадей поверхности заготовки и детали. Для расчёта исходной площади полуфабриката его поверхность следует рассмотреть как состоящую из трёх сегментов:
- донце
S = я( - r )2, (1)
- обечайка
S2 = 2яR (H + А-), (2)
- торовый сегмент (скругление)
я
2
S3 = 2яг1 J[R - r (1 - cos (p)\dy = ж2г1 (R - r1) + 2я;2. (3)
0
Здесь R - радиус кривизны борта в плане; А - припуск высоты борта; r - радиус скругления борта заготовки после первого перехода (рис. 2, а).
Рис. 2. Параметры геометрии полуфабриката на промежуточных стадиях штамповки: а - после первого перехода; б - перед окончательной штамповкой
Суммируя выражения (1) - (3), получим полную площадь поверхности заготовки
Я0 = ж(Я - г )2 + 2жЯ (Н + А- г ) + ж2гг (Я - г) + 2жг2. (4)
Площадь поверхности заготовки в начальной стадии посадки на формблок аналогично имеет три составляющих:
£ = ж(Я - г )2 + 2жЯ (Н - г0) + Я *;
(5)
площадь волны избыточного материала над кромкой формблока, которая имеет вид то-роподобного сегмента (рис. 2, б), определяется соотношением
2
£ * = 2жг01 [ Я - г0 + р(р) cos р]
Р2 (р) +
йр(р) йр
йр.
(6)
Профиль поперечного сечения этой волны из условия гладкого примыкания её ветвей к донцу и борту на основании результатов экспериментов и с достаточной для практики степенью точности может быть аппроксимирован зависимостью
Р(<) =
1 + Иcos2
21 р-
ж
(7)
где г0 - радиус скругления кромки формблока; И - высота волны в долях г0 (рис. 2, б).
Подставляя аппроксимирующую функцию (7) в подкоренное выражение зависимости (6) и разлагая полученную дробно-степенную функцию в ряд по степени И (И < 1) , получим
£ * (И) = 2жг01 [ Я - г0 + р(р) cos р( х
х|1 + Иcos2
(
21 р-
ж
Л
+ 8И2^2
21 р-
ж
Л (
sin
21 р-
ж
(8)
-8И3^4
,21'"
Л (
sin
ж
Л!
2 1р-- I \йр.
После группирования подынтегрального выражения (8) по степеням И и сведения ряда к виду а0 + а1И + а2 И2 + а3И3 последовательно получаются следующие зависимости:
аг
= Я - г0 (1 - cos р).
а
= (Я - г0) cos2 21 р - Ж I + 2г0 cos рcos2 21 р - —
ж
ж
а2 = 8[Я-г0 (1 -соБр)]cos2 2|р- — J ;
Л ( sin2
21 р-
п
л
+ г0 cos рсоБ
а3 = -8
[ Я-Г0 (1 -СОБр)] cos4 21 р-П
V
(
21 р-
п
Sin
Л
,2
(
-г0 СОБ (СОБ
,2р
Л (
Sin
Л"
,2р
= -8(я-г0)соб4I 21 р- —
п
(
Sin
Л
,21Р"
Таким образом, формулу (8) можно представить в виде
5* (И) = 2пг0 (А + Л1И + Л2И2 + Л3И3), где коэффициенты Л0, Л1, Л2 и А3 имеют вид
(9)
Л =(Я - Г0 )0 + Г0/l, Л1 = (Я - Г )2 + 2Г0 ¡3, Л2 = 8 (Я - Г0 )14 + Г0 (815 +16), Л3 =-8 (Я -0 )7.
Здесь коэффициенты ¡i ( = 1,7) являются тригонометрическими интегралами, которые имеют следующие значения:
¡0 = JdР = — , 02
п
2
¡1 = | cosрdр = 1
п
2 (
¡2 = | соб2 2 I р-П I р = —
¡3 = |соб(р)соб2 2\р-п} ^р=15,
п
I (
п 4
¡4 = [соб2 21 р-
п
¡5 =[ СОБ (р)сОБ2 2 \р- —
\\
2 (р-£ ^ с!р=Ж
( ( —Л 2 (
Sin
п
16
Л
21 р-— I dр = —,
и
63
'«= 1
= J cos (pcos
0 п
1 } 4 ( п
17 = Jcos 21 р--
(( п^ 128
21р-^| ар=—
Sin
п) \ , п
21 р--I ар = —.
Г 4 I, 32
Таким образом, коэффициенты А0, А1,А2, А3 примут вид
А =П (Я - г)+^
А1 =п (Я - г0 ) + ^ r0, А2 = | ( - Г0 )+ 31? Г0, Аз =ПП ( - Го )•
После подстановки А0, А1, А2, А3 выражение (9) примет вид
о*/7\ ~ \п п п-2 (п_ 64- 15п
Я (И)= 2пг0Я--г0 +1 — Я +-
0 1 2 2 0 ^ 4 60
Го ) И+(? Я+Г0 ] И2-£ (Я - Г0 )И31.(10)
Далее, подставляя зависимость (10) в выражение (5) и приравнивая результат к начальной площади (4), после преобразований получаем окончательное выражение для припуска высоты борта полуфабриката:
Д = г0 (а- ^
п п -1„/ ,4
2 - 2Р(а +1)
, I п 64 - 15п „
г0 И ^ —+-р-
1 4 60
п+65 е) и-п (1 -Р) },(11)
где
а = ± > 1,
Р = -° < 1. Я
г
0
Величину припуска высоты борта полуфабриката необходимо ограничить некоторым диапазоном, имеющим верхнюю и нижнюю границы. Верхнюю границу припуска (11) следует определить из условия сохранения устойчивости волны избыточного материала при её осадке (обжиме) на кромку формблока, чтобы не допустить её переход в складку. В качестве критерия потери устойчивости предлагается принять равенство поперечного давления и продольного давления на волну, а конкретнее - равенство их равнодействующих Q и Р, возникающих от равномерно распределённого по внешней поверхности волны нормального давления р (рис. 3).
Рис. 3. Схема силового воздействия эластомера на волну избыточного материала
Величины Q и Р пропорциональны суммарным проекциям равномерного нормального давления эластомера р на соответствующие взаимно перпендикулярные направления. Математически это выражается в виде интеграла скалярного произведения вектора единичной нормали к поверхности волны и соответствующей единичной орты. Последние естественным образом определяются как е = (1,0), i = (0,1) .
Вектор единичной нормали при данной ориентации осей имеет координаты
(
n =
Л
d_ dp
(p(p) sin (p) )
d_ dp
(р(р) cos (р))
р2(р)+
d р(р) dp
Р2(Р)+
d р(ф) dp
(12)
J J
Подставляя в выражение (12) аппроксимирующую функцию (7) со сдвигом по уг-
— й Р Л Q
ловой координате р = р--, для приведенных величин Р =-, Q =- получим их
4 Г0Р Г0Р
значения.
Продольный компонент давления определяется следующим образом:
1 4 4
P = — J(n, е) dLp = Jjcos (p([l + hcos2 (2p) J- 4h sin (p) cos (2p) sin (2p) J dp
r0 0 0
Аналогично поперечный компонент давления определяется как
1 4 4
Q (h) = — J(n,Í) dL (p) = J |4h cos (p( cos (2p( sin (2p) + sin (p) [l + hcos2 (2p)]Jdp .
r0 0 0
В результате расчётов определено:
р= 0(И) = 1 *И
Приравнивая полученные значения, получим верхнюю границу относительной высоты волны:
= 42 " 1.
С учётом этого результата из выражения (11) следует, что верхняя граница припуска определяется как
Л тах- = ((Х" 1)
Г0
Ж Ж" 1 ,ч
2 "Ж ?(а+1)
+((- •))*^КЖ*65 " 1)" 4 (1 "Л(3 " 2Г2
Нижнюю границу припуска Лтт можно определить из решения деформационной
краевой задачи теории пластичности. В качестве схемы очага деформации для упрощения задачи допустимо принять вершину волны избыточного материала за цилиндрический сегмент, т.е. представить его в виде толстостенной трубы, осаживаемой опоясывающим давлением на стержневую оправку радиуса г0 при наличии однородного осевого сжатия. Деформированное состояние в очаге деформации будет определяться тремя главными компонентами «р,«р,« . Радиальный и окружный компоненты зависят
только от радиальной координаты и связаны уравнением неразрывности:
«^М^- (13)
Осевой компонент с точностью до малой порядка ^¡К, где Я - толщина заготовки, равен
« =-Нр. (14)
Вторым уравнением служит условие несжимаемости:
= (15)
Граничным условием логично считать контакт заготовки с кромкой формблока в конце процесса посадки волны избыточного материала, т.е. положить радиальное перемещение на внутренней поверхности волны и (г0 (1 + И)) = -г0к . Отсюда следует, что
«.(о**)) и¥#=-т\ (16)
4 ' г0 (1 + И) 1 + И
Подставляя уравнение неразрывности (13) и осевой компонент деформации (14) в условие несжимаемости (15), с учётом выражения (16), получим краевую задачу Коши в виде следующей системы уравнений:
ра-^1 + 2Ер(р)-РК
d р
( (1 + И))
И
1 + И
где г0 (1 + И) < р < г0 (1 + И) + ^.
В результате интегрирования получим
% = р - ^ (1+^ Р
(17)
где % = Р - безразмерная координата.
Условием полного перехода волны избыточного материала в пластическое состояние является выход границы зоны пластичности на внешнюю границу очага деформации р = г0 (1 + И) + % . Согласно критерию Губера-Мизеса в деформациях интенсивность деформаций, определяемая как
л/3
2
(18)
в пластичном состоянии при любой механической схеме напряжений-деформаций
- 77
должна достигать критической величины = = —, где ог и Е -предел текучести и
Е
модуль упругости материала, соответственно.
Подставляя выражения (13), (14), (17) в зависимость (18) и введя обозначение
% = 1 + И + 7], где 7] = — - относительная толщина заготовки, получаем, что для внешней Г0
поверхности волны избыточного материала должно выполняться равенство
И 2 (1+И) 4 (1+1+И Р
Ъ =
(1+И+])
3 РИ + И (1+И) 2 (1++ 2 (1 + И + ] V 2
Е
3
-РИ -
2 (1 + И + ])
И (1 + И) Г +1 + И
[ + И + ])2 V 2
Р
Возводя обе части последнего соотношения в квадрат, после упрощающих преобразований получаем выражение для нижней границы И в виде неявной зависимости:
и2 (1+и)2 г + 1+и
(1 + И + ] V 2
+И р| + 3 „2р2 =3
4 Е2
(19)
Пренебрегая членами, порядок малости которых выше Р, сведем зависимость (19) к квадратному уравнению, решением которого является явное выражение для минимального значения Итт, при котором волна избыточного материала в полном объеме перей-
дет в пластическое состояние
г
0
К
^(1 +г) я-4 р) +(1 +г)2 я(1 +я)-(1 +г) я +1
1 + я
(20)
у/э С
где Я =--- « 1. После упрощающих преобразований выражение (20) для нижней
2 Е
границы примет вид
л/з ^ (1+г)2
К ■ = г- •
ЛИ 2 _ /Т
2 Т^р-^^ (1+г)2
Подставляя найденное значение Ктт в зависимость (11), получим величину минимального припуска высоты борта Лтт, необходимую для перевода в пластическое
состояние волны избыточного материала.
На рис. 4 в безразмерном виде представлены величины минимального и максимального припуска высоты борта в зависимости от соотношений толщины стенки заготовки Я, радиуса гиба борта г0 и радиуса борта в плане Я .
Рис. 4. Зависимости припусков высоты борта от геометрических параметров детали: а - максимально возможный припуск Лтах /г0 ;
б - минимально возможный припуск 100 (Лт;п /г0) при г = Я/г0 = 0,5 ; в - минимально возможный припуск 100 (Лт;п/г0) при г = Б/г0 = 0,05
Из графиков на рис. 4 следует, что припуски находятся в обратной зависимости от кривизны детали в плане. Практические значения припусков высоты борта, например для деталей из алюминиевых сплавов, когда обычно принимается г0 = 3Я, составляют величины порядка 0,5.. .0,9 толщины стенки заготовки.
Заключение
В результате исследований разработанной перспективной технологии стеснённого изгиба листовых деталей летательных аппаратов с криволинейными бортами эластичной средой получена аналитическая модель определения припуска, связывающая геометрию детали и утолщение стенки заготовки с величиной припуска. Величина припуска ограничивается некоторым диапазоном, нижняя граница которого определена из условия начала пластических деформаций, верхняя - из условия возможности потери устойчивости, приводящей к неисправимому браку. Выяснено, что величина припуска находится в обратной зависимости от кривизны детали в плане. Для наиболее распространённых в авиационной и ракетно-космической технике конструкций деталей припуски составляют величины порядка 0,5...0,9 толщины стенки заготовки.
Библиографический список
1. Моисеев В.К., Мантусов М.Н., Плотников А.Н., Ломовской О.В., Шаров А. А. Совершенствование технологии гибки криволинейных бортов деталей самолётов // Известия Самарского научного центра РАН. 2018. Т. 20, № 4 (3). С. 441-444.
2. Мантусов М.Н., Моисеев В.К., Шаров А.А., Громова Е.Г., Рыжаков С.Г. Штамповка листовых деталей с криволинейными бортами эластичным материалом // Известия Самарского научного центра РАН. 2018. Т. 20, № 4 (3). С. 332-336.
3. Гречников Ф.В., Моисеев В.К., Ломовской О.В., Шаров А.А., Мантусов М.Н. Штамповка с тангенциальным сжатием высокоточных листовых деталей // Кузнечно-штамповочное производство. Обработка материалов давлением. 2019. № 7. С. 3-9.
4. Moiseyev V.K., Sharov A.A., Gromova E.G., Mantusov M.N. Results of the sheet parts curved edges constrained bending with elastomer // Key Engineering Materials. 2017. V. 746. P. 285-289. DOI: 10.4028/www.scientific.net/kem.746.285
ANALYTICAL MODEL FOR DETERMINING ALLOWANCES FOR ELASTOMER CONSTRAINED BENDING OF AIRCRAFT PARTS
© 2020
V. K. Moiseev
E. G. Gromova
O. V. Lomovskoi
М. N. Mantusov
A. N. Plotnikov
A. A. Sharov
Doctor of Science (Engineering), Professor; Department of Aircraft Production
and Quality Control in Mechanical Engineering;
Samara National Research University, Samara, Russian Federation;
moiseevvk@mail.ru
Candidate of Science (Engineering), Associate Professor; Department of Aircraft Production and Quality Control in Mechanical Engineering; Samara National Research University, Samara, Russian Federation; pla.gromova@ya.ru
Candidate of Science (Engineering), Associate Professor; Department of Aircraft Production and Quality Control in Mechanical Engineering; Samara National Research University, Samara, Russian Federation; oleg.lomovskoi@ya.ru
Leading Engineer;
Ulyanovsk branch of Tupolev Design Bureau, Ulyanovsk, Russian Federation; mishaufkb@ya.ru
Candidate of Science (Engineering), Associate Professor; Department of Aircraft Production and Quality Control in Mechanical Engineering; Samara National Research University, Samara, Russian Federation; anplotnikov@ya.ru
Candidate of Science (Engineering), Associate Professor; Department of Aircraft Production and Quality Control in Mechanical Engineering; Samara National Research University, Samara, Russian Federation; aa.sharov@ya.ru
In the article, the authors present a method of constrained bending of aircraft sheet parts with curved sides by an elastic medium. The process of conducting constrained bending is described. A schematic description of the process is presented. The manufacture of a part is carried out in two steps. As a result of the first step, a piece is obtained with thinning in the radius part, the second step is accomplished with a unit die of reduced height. The resulting wave of excess material is deformed along the radius of the bending mandrel. In the finished part, an increase in the thickness of the workpiece is observed in the zone of the bending radius. The wave of excess material is formed due to a special allowance, the size of which is the critical factor for achieving this increase in thickness. The analytical model for determining the allowance is based on generally accepted assumptions and principles for calculating process parameters. To determine the allowance, an approximation of the shape of the excess wave by a tangential function is proposed. This allows you to get an analytical dependence that relates the geometry of the part and the thickening of the wall to the size of the allowance. The size of the allowance is limited to a certain range, the lower boundary of which is determined by the onset of plastic deformations, the upper one - by the possibility of loss of stability, leading to irreparable spoilage. The results are illustrated by graphs of the dependence of the minimum and maximum allowance for the side height on the ratios of the workpiece wall thickness, the radius of the side bending and the radius of the side in plan. Characteristic curves and analytical dependencies are presented in a dimensionless form, where the geometric parameters are assigned to the radius of the bending tool, that is, to the radius of bending of the part side on the inner surface.
Aircraft production technology; elastomer forming; constrained bending; deformations; allowance
Citation: Moiseev V.K., Gromova E.G., Lomovskoi O.V., Mantusov M.N., Plotnikov A.N., Sharov A.A. Analytical model for determining allowances for elastomer constrained bending of aircraft parts. Vestnik of Samara University. Aerospace and Mechanical Engineering. 2020. V. 19, no. 3. P. 73-84. DOI: 10.18287/2541-7533-2020-19-3-73-84
References
1. Mosieev V.K., Mantusov M.N., Plotnikov A.N., Lomovskoi O.V., Sharov A.A. Improving the technology of bending curved sides of aircraft parts. Izvestiya Samarskogo Nauchnogo Tsentra RAN. 2018. V. 20, no. 4 (3). P. 441-444. (In Russ.)
2. Mantusov M.N., Moiseev V.K., Sharov A.A., Gromova E.G., Ryzhakov S.G. Stamping of sheet details with curvilinear sides by elastic material. Izvestiya Samarskogo Nauchnogo Tsentra RAN. 2018. V. 20, no. 4 (3). P. 332-336. (In Russ.)
3. Grechnikov F.V., Moiseev V.K., Lomovskoy O.V., Sharov A.A., Mantusov M.N. Stamping with tangential compression of high-precision sheet metal parts. Forging and Stamping Production. Material Working by Pressure. 2019. No. 7. P. 3-9. (In Russ.)
4. Moiseyev V.K., Sharov A.A., Gromova E.G., Mantusov M.N. Results of the sheet parts curved edges constrained bending with elastomer. Key Engineering Materials. 2017. V. 746. P. 285-289. DOI: 10.4028/www.scientific.net/kem.746.285
