Design of Thin-Walled Single-Curvature Parts for Use in Lightweight Structures
- Authors: Morozov Y.A.1, Belelyubskiy B.F.2
-
Affiliations:
- Bauman Moscow State Technical University
- Moscow Polytechnic University
- Issue: Vol 20, No 4 (2024)
- Pages: 342-354
- Section: Analysis and design of building structures
- URL: https://journals.rudn.ru/structural-mechanics/article/view/41544
- DOI: https://doi.org/10.22363/1815-5235-2024-20-4-342-354
- EDN: https://elibrary.ru/TCWWPU
- ID: 41544
Cite item
Full Text
Abstract
The aim of the study - the purpose of the study was to find the minimum (critical) curvature of sheet material, to which it can be bent without fracture (formation of longitudinal cracks) and which is determined by the combined «play» of two deformational parameters: thinning, responsible for cross-section weakening, and strain hardening of the material, characterized by the intensity of deformations. The existing sheet bending pattern is analyzed with regard to the kinematics of deformational changes in the initial radii of the part due to the continuity of compressive (radial) and tensile (tangential) deformations. Assuming the Bernoulli’s hypothesis in sheet bending conditions, a mathematical model has been developed for estimating the deformational and geometric (thinning) parameters during the formation of a torus surface of various curvatures. The level of radial stresses has been identified taking into account strain hardening and thinning of the bent material, which lead to the exhaustion of its load-bearing capacity (fracture), where the plasticity criterion is the mechanical properties of a particular material obtained in tensile tests (yield and strength limits, relative elongation), approximated by a power law. The obtained results can be applied in the design of lightweight power structures; in modeling the stressstrain state of metal when developing technological processes of sheet stamping (bending) for calculating the magnitude of thinning, assessing the level of radial stresses in metal bending along the end edge of a pressing punch, as well as when designing bending equipment.
Full Text
1. Введение Современное строительство зданий и сооружений предусматривает использование облегченных силовых элементов, среди них светопрозрачные фасадные или навесные системы, образованные стеклянными панелями, упакованными каркасной обвязкой, которую чаще всего выполняют из алюминиевых сплавов или стали (в зависимости от интенсивности действующих нагрузок). Изготовление алюминиевых профилей прессованием позволяет спроектировать любую их конфигурацию, но высокая стоимость инструментальной оснастки также отражается и на стоимости подобных изделий. Производство стальных профилей не требует какого-то инструмента индивидуального исполнения и может осуществляться различными схемами формоизменения (гибка, профилирование и пр.). В этом отношении изготовление гнутых профилей является аналогичным производству тонкостенных деталей одинарной кривизны методами листовой штамповки, обеспечивающей требуемую конфигурацию деталей, которые широко применяются в авиа-, судо-, автопромышленности [1; 2]. При этом актуальным является вопрос определения минимально возможного радиуса изгиба согнутой детали, который в пределе стремится к нулю (рис. 1). Иногда уменьшение предельно допустимого радиуса изгиба диктуется необходимостью снижения материалоемкости изделия, условиями монтажа и др. При этом использование чрезвычайно малых радиусов гибки приближает конфигурацию листовой профильной детали (уголки, швеллеры и др.) к сечению прессованных профилей, т.е. к показателям большей жесткости [3; 4]. Рис. 1. Схема V-образной гибки в штампе И с т о ч н и к: https://ugselmash.ru/uslugi/gibka-metalla/gibka-stali; https://ipmet.ru/razd/metalloobrabotka/gibka-metalla Figure 1. V-die bending S o u r c e: https://ugselmash.ru/uslugi/gibka-metalla/gibka-stali; https://ipmet.ru/razd/metalloobrabotka/gibka-metalla Основной задачей в данном случае является предотвращение разрушения заготовки при ее изгибе, обычно происходящее на выпуклой поверхности радиуса R (по биссектрисе угла), вследствие больших деформаций растяжения и утонения материала (рис. 2)[7]. Рис. 2. Сечение изогнутой заготовки: 1, 2 - зоны тангенциальных напряжений растяжения и сжатия И с т о ч н и к: выполнено Ю.А. Морозовым Figure 2. Section of a curved workpiece: 1, 2 - zones of tangential tensile and compressive stresses S o u r c e: made by Yu.А. Morozov Утонение приводит к снижению прочности детали и является дополнительной причиной концентрации напряжений и усталостного разрушения (обычно допустимое утонение ограничивается 20 %). Технологически процесс гибки (изменение кривизны листового материала) является весьма простой операцией, достаточно освященной в работах отечественных[8] и зарубежных авторов, предлагающих различные методики определения утонения материала и его деформационного упрочнения [5-9]. Однако решение задачи минимального радиуса изгиба требует комплексного рассмотрения кинематики процесса и анализа напряженно-деформированного состояния металла, обусловленного совокупной «игрой» двух деформационных параметров - утонение, приводящее к ослаблению сечения детали, и деформационное упрочнение материала, характеризуемое интенсивностью деформаций [10-14]. Следует отметить, что современное программное обеспечение позволяет решать подобные задачи в диалоговом режиме компьютерного моделирования с использованием численных методов или метода конечных элементов [15; 16]. При этом следует учесть, что основой каждого подобного ПО является математический алгоритм, основанный на определенных моделях и допущениях. В связи с этим рассмотрим поведение металла в указанных условиях формоизменения. 2. Метод Формоизменение гибкой приводит к появлению в материале двух противоположных деформационных процессов: растяжения и сжатия, обуславливающих деформационное изменение первоначальных радиусов R1 и R2 (рис. 3)[9] [17; 18]. Рис. 3. Распределение тангенциальных деформаций по сечению полосы И с т о ч н и к: выполнено Ю.А. Морозовым Figure 3. Distribution of tangential deformations over the strip section S o u r c e: made by Yu.А. Morozov Наружный радиус R1 под действием радиальной деформации εr нар < 0 уменьшается, принимая новое условное значение Rнар. Внутренний же радиус R2 увеличивается до Rвн под действием соответствующей радиальной деформации εr нар > 0 Rнар = +R1 (R1 -ρн)εr нар Rвн =R2 + ρ -( н R2 )εr вн . (1) Величина утонения будет определяться величиной тангенциальных деформаций, действующих в растянутых и сжатых слоях заготовки, вследствие чего первоначальные радиусы R1 и R2 принимают новые продеформированные значения - соответственно Rнар и Rвн с границей раздела в нейтральном сечении ρн , в котором тангенциальные деформации εθ равны нулю: εθнар = Rнар -ρн = Rнар -1 ρн ρн . = ρ -н Rвн = -1 Rвн ε θвн ρн ρн (2) Совместное решение (1) и (2) устанавливает выражения наружного и внутреннего радиусов (при допущении фиксированного внутреннего радиуса R2 = Rвн): +1, (3а) R , (3б) ρн Z1 где Z1 = ρнR1 - коэффициент, определяющий положение нейтральной поверхности. Тогда из условия равенства по модулю тангенциальных деформаций в растянутых и сжатых слоях гнутого элемента получаем деформационную модель изгиба, позволяющую последовательно определить положение нейтральной поверхности Z1 в сечении и установить окончательную толщину полосы для разных значений первоначальных радиусов R1R2 : R2 - - Z εθнар =-εθвн или (2-Z1)- 4Z1 3Z12 = R1 1 . (4) 4Z1-3Z12 -Z1 Z1 С использованием величины утонения при оценке радиальных напряжений изгибаемого материала по торцевой кромке давящего пуансона, можно рассчитать предельный и критический радиус кривизны, допускающий гибку без разрушения гнутого элемента (образование продольных трещин). Рассмотрим гибку детали на кривизну R1R2 =1,1, определяющую следующий характер равновесия тангенциальных деформаций (4): 1 (2 -Z1)- 4Z1 -3Z12 = 1,1 -Z1 . 4Z1 -3Z12 -Z1 Z1 Итерационным перебором определяется коэффициент положения нейтральной поверхности Z1= 0,9535 (рис. 4). Рис. 4. Коэффициент нейтральной поверхности, Z1 И с т о ч н и к: выполнено Ю.А. Морозовым, Б.Ф. Белелюбским Figure 4. Neutral surface factor, Z1 S o u r c e: made by Yu.А. Morozov, B.F. Belelyubskiy Подставляя данное значение в выражения (3а) и (3б), устанавливается относительная величина деформированных радиусов: ρн Z1 0,9535 При исходной толщине листового металла S =0,8 мм и фиксированном внутреннем радиусе гибки Rвнρн =const устанавливается радиус торового скругления давящего пуансона и радиус нейтральной поверхности S 0,8 rп = R2 = Rвн = = = 8 мм; R R1 2 -1 1,1 1- ρ =н Z R1 1 = Z1 R1 R2 = 0,9535 1,⋅ 1 8⋅ = 8,39 мм. R2 Наружный радиус торового скругления R Rнар ρн 1,0466 8,39 8,78 мм. Коэффициент утонения (рис. 5) [19] Рис. 5. Коэффициент утонения листового металла, S'/S И с т о ч н и к: выполнено Ю.А. Морозовым, Б.Ф. Белелюбским Figure 5. Thinning factor of sheet metal, S'/S S o u r c e: made by Yu.А. Morozov, B.F. Belelyubskiy Упрочнение материала при гибке определяется интенсивностью деформаций, т.е. суммарным значением тангенциальных деформаций на наружной и внутренней поверхностях изгибаемой полосы εi = εθнар +εθвн = Rнар - Rвн =1,0466-0,9534 = 0,093 (9,3 %). ρн ρн Принимаем модельным материалом сталь 20[10] [20], пластичность которой будет устанавливаться механическими свойствами, полученными в испытаниях на растяжение и аппроксимированными степенной зависимостью σТ = σ0,2 + Aεin = 245+ 22,4 9,⋅ 30,62 = 334,3 МПа, где σ0,2 = 245 МПа[11] - условный предел текучести стали 20; А, n - коэффициенты упрочнения материала [21]. n; A , где σв = 410 МПа - предел прочности стали 20; δ=25 % - относительное удлинение при разрыве. Радиальное напряжение при V-образной гибке листового материала на угол 90 градусов (α = π2 ) с использованием коэффициента контактного трения f = 0,2 (табл. 1) [22]: σρmax = 2σρТ S S′S e fα = 2 8,39 0,976334,3⋅ 0,8 e0,2π2 = 22,4 МПа, н где σТ - сопротивление пластической деформации. Таблица 1 Напряженно-деформированное состояние металла при гибке Кривизна изогнутого элемента R1R2 Коэффициент нейтральной поверхности Z1 = ρнR1 Коэффициент утонения S′ S Интенсивность деформаций εi , % Радиальные напряжения σρ max , МПа 1,1 0,9535 0,976 9,30 22,40 1,3 0,8779 0,938 24,8 78,50 1,5 0,8192 0,919 37,2 137,9 2,0 0,7171 0,871 60,6 290,0 2,5 0,6514 0,840 77,2 432,1 3,0 0,6055 0,819 89,9 558,2 И с т о ч н и к: выполнено Ю.А. Морозовым Table 1 The stress-strain state of the metal during bending Curvature of a curved element R1R2 Neutral surface factor Z1 = ρнR1 Thinning factor S′ S Deformation intensity εi , % Radial stresses σρ max , МПа 1.1 0.9535 0.976 9.30 22.40 1.3 0.8779 0.938 24.8 78.50 1.5 0.8192 0.919 37.2 137.9 2.0 0.7171 0.871 60.6 290.0 2.5 0.6514 0.840 77.2 432.1 3.0 0.6055 0.819 89.9 558.2 S o u r c e: made by Yu.А. Morozov Предельный радиус кривизны, допускающий гибку без разрушения гнутого элемента, будет при достижении радиальных напряжений начального предела текучести материала σρmax ≤ σ0,2 = 245 МПа. Критический радиус кривизны, характеризуемый началом разрушения гнутого элемента (образование продольных трещин), будет при радиальных напряжениях, сравнимых с пределом прочности данного модельного материала (рис. 6): σρmax < σв = 410 МПа. Рис. 6. Радиальные напряжения при изгибе листового материала (сталь 20) И с т о ч н и к: выполнено Ю.А. Морозовым Figure 6. Radial stresses during bending of sheet material (steel 20) S o u r c e: made by Yu.А. Morozov Аппроксимируя распределение радиальных напряжений (среднеквадратичная ошибка определения R2 =0,999), устанавливается критический радиус кривизны R1R2 ≈ 2,42: 2 R1 σρmax =-10 R2 +324 RR12 -320,6 = = -10 2,42⋅ 2 +324 2,42⋅ -320,6 ≈ 405 МПа. Для оценки влияния пластических свойств различных металлов ниже рассматриваются некоторые стали, используемые для производства листоштампованных изделий и имеющие различные характеристики прочности и пластичности (табл. 2). Таблица 2 Механические свойства материалов для листовой штамповки Параметр Материал Полоса, ГОСТ 1577-931 Прокат, ГОСТ 1050-882 Сталь 08 Сталь 40 Сталь 15 Сталь 20 Условный предел текучестиσ0,2 , МПа 196 335 225 245 Предел прочности σв , МПа 320 570 370 410 Относительное удлинение δ, % 33 19 27 25 Коэффициенты упрочнения A / n 6,35/0,85 60,7/0,46 14,9/0,69 22,4/0,62 И с т о ч н и к: выполнено Ю.А. Морозовым Table 2 Mechanical properties of sheet stamping materials Parameter Material Strip, GOST 1577-93[12] Rolled steel, GOST 1050-88[13] Steel 08 Steel 40 Steel 15 Steel 20 Proof strength σ0,2 , МПа / MPa 196 335 225 245 Tensile strength σв , MPa 320 570 370 410 Relative elongation δ, % 33 19 27 25 Hardening factors, A / n 6.35/0.85 60.7/0.46 14.9/0.69 22.4/0.62 S o u r c e: made by Yu.А. Morozov 3. Результаты и обсуждение Анализ интенсивности распределения радиальных напряжений с учетом пластических свойств рассмотренных материалов устанавливает практически одинаковые критические радиусы Рис. 7. Радиальные напряжения при изгибе листового материала И с т о ч н и к: выполнено Ю.А. Морозовым, Б.Ф. Белелюбским Figure 7. Radial stresses during bending of sheet material S o u r c e: made by Yu.А. Morozov, B.F. Belelyubskiy Найденное значение кривизны, также можно представить фактической величиной критического внутреннего радиуса S S rкр ≥ =R = =(0,71K0,69) S . 4. Заключение В статье проведено исследование минимального (критического) внутреннего радиуса кривизны гнутого элемента, допускающего гибку листового материала без разрушения (образование продольных трещин) на основе математической модели напряженно-деформированного состояния формоизменяемой заготовки. На основании проведенного исследования можно сформулировать следующие выводы: 1. На основе гипотезы плоских сечений представлен метод высчитывания критического радиуса, определяемый интенсивностью деформаций на наружной и внутренней поверхностях изгибаемой заготовки (полосы). 2. Численное моделирование, учитывающее деформационные процессы и механические характеристики прочности и пластичности материала, показало устойчивый рост радиальных напряжений, вызывающих разрушение материала при достижении определенной величины критического радиуса кривизны. Таким образом, предложены математически обоснованные рекомендации для разработки технологических процессов листовой штамповки или проектирования гибочной оснастки с учетом пластических свойств конкретного материала, что позволит снизить материалоемкость изделия с повышением его жесткости. Для исключения значительного объема итерационных расчетов, приводятся соответствующие аппроксимационные зависимости, позволяющие с высокой точностью 2…3 % определить положение нейтральной поверхности гнутого элемента и величину утонения материала.About the authors
Yury A. Morozov
Bauman Moscow State Technical University
Author for correspondence.
Email: akafest@mail.ru
ORCID iD: 0000-0001-9229-7398
SPIN-code: 3189-5426
PhD in Technical Sciences, Associate Professor of the Department of Materials Processing Technologies (MT-13)
Moscow, RussiaBoris F. Belelyubskiy
Moscow Polytechnic University
Email: alib@bk.ru
ORCID iD: 0000-0002-1702-707X
SPIN-code: 2007-1003
PhD in Technical Sciences, Associate Professor of the Department of Metallurgy
Moscow, RussiaReferences
- Vlasov S.V., Yelatontsev N.A. Balans napryazheniy i deformatsiy pri kholodnoy gibke listovoy sudostroitel’noy stali. FEFU: School of Engineering Bulletin. 2021;(1):36–48. (In Russ.) http://www.doi.org/10.24866/2227-6858/2021-1-4
- Dang X., He K., Zhang F., Du R. A new flexible sheet metal forming method of incremental bending. Procedia Manufacturing. 2018;15:1298–1305. https://doi.org/10.1016/j.promfg.2018.07.355
- Morozov Yu.A. Development of the configuration of bent profiles in the design of translucent structures. Informatics and technologies. Information technologies in industry and informatics. Proceedings of the conference. Moscow, April 12–13, 2018. Moscow: RTU MIREA Publ.; 2018;2:733–737. (In Russ.) EDN: YWQWPB
- Morozov Yu.A. Investigation of the deformed state of the material in the production of bent profiles. Informatics and technologies. Information technologies in industry and informatics. Proceedings of the conference. Moscow, April 11–12, 2019. Moscow: RTU MIREA Publ.; 2019;2:288–295. (In Russ.) EDN: IFOSTI
- Ahn K. Plastic bending of sheet metal with tension/compression asymmetry. International Journal of Solids and Structures. 2020;204–205:65–80. https://doi.org/10.1016/j.ijsolstr.2020.05.022
- Barnwal V.K., Lee S.-L., Jisik Choi, Kim J.-H., Barlat F. Fracture assessment in dual phase and transformationinduced plasticity steels during 3-point bending. Theoretical and Applied Fracture Mechanics. 2020;110:102834. https:// doi.org/10.1016/j.tafmec.2020.102834
- Zadpoor A.A., Campoli G., Sinke J., Benedictus R. Fracture in bending — The straining limits of monolithic sheets and machined tailor-made blanks. Materials & Design. 2011;32(3):1229–1241. https://doi.org/10.1016/j.matdes.2010.10.005
- Yoshida M., Yoshida F., Konishi H., Fukumoto K. Fracture limits of sheet metals under stretch bending. International Journal of Mechanical Sciences. 2005;47(12):1885–1896. https://doi.org/10.1016/j.ijmecsci.2005.07.006.
- Romanovskiy V.P. Handbook of Cold Forming. Moscow. Leningrad: Mashinostroyeniye. Publ.; 1979. (In Russ.)
- Li F.F., Zhu J., Zhang W., Fang G. Investigation on the inhomogeneous deformation of magnesium alloy during bending using an advanced plasticity model. Journal of Materials Research and Technology. 2023;25:5064–5075. https:// doi.org/10.1016/j.jmrt.2023.06.264
- Li S., He J., Gu B., Zeng D., Xia Z.C., Zhao Y., Lin Z. Anisotropic fracture of advanced high strength steel sheets: Experiment and theory. International Journal of Plasticity. 2018;103:95–118. https://doi.org/10.1016/j.ijplas.2018.01.003
- Soyarslan C., Malekipour Gharbi M., Tekkaya A.E. A combined experimental-numerical investigation of ductile fracture in bending of a class of ferritic-martensitic steel. International Journal of Solids and Structures. 2012;49(13): 1608–1626. https://doi.org/10.1016/j.ijsolstr.2012.03.009
- Stoughton T.B., Yoon J.W. A new approach for failure criterion for sheet metals. International Journal of Plasticity. 2011;27(3):440–459. https://doi.org/10.1016/j.ijplas.2010.07.004
- Levy B.S., Van Tyne C.J. Predicting breakage on a die radius with a straight bend axis during sheet forming. Journal of Materials Processing Technology. 2009;209(4):2038–2046. https://doi.org/10.1016/j.jmatprotec.2008.04.053
- Bate K., Vilson E. Numerical methods in finite element analysis. Prentice-Hall Publ.; 1976. Available from: https://sciarium.com/file/268214/ (accessed: 02.03.2024).
- Zenkevich O.K. The finite element method in engineering. Moscow: Mir Publ.; 1975. (In Russ.) Available from: https://djvu.online/file/DtUw9BqXrtZCc (accessed: 02.03.2024).
- Lukashkin N.D., Kokhan L.S., Punin V.I., Morozov Yu.A. Bending of profiles on presses and mills. Moscow: MGVMI Publ.; 2005. (In Russ.)
- Kokhan L.S., Roberov I.G., Morozov Yu.A. Investigation into kinematic parameters during bending the sheet materials. Tekhnologiya metallov. 2008;(10):11–13. (In Russ.) EDN: IVMCXK
- Morozov Yu.А. The study of marginal deformations of the leaf extracts with regard to plastic thinning and destruction of the material. Structural Mechanics of Engineering Constructions and Buildings. 2019;15(5):353–359. (In Russ.) https://doi.org/10.22363/1815-5235-2019-15-5-353-359
- Arzamasov B.N., Solovyova T.V., Gerasimov S.A. Handbook of Structural Materials. Moscow: MSTU named after N.E. Bauman Publ.; 2005. (In Russ.)
- Tret’yakov A.V., Zyuzin V.I. Mechanical properties of metals and alloys during pressure treatment. Directory. Moscow: Metallurgiya Publ.; 1973. (In Russ.)
- Isachenkov E.I. Contact friction and lubrication in metal forming. Moscow: Mashinostroyeniye Publ.; 1978. (In Russ.)
Supplementary files










