CC BY
78
UDC 621.45
Doi: 10.31772/2587-6066-2019-20-2-243-250
For citation: Ushakova E. S. [Modeling of the stress-strain state of rocket-space technology structural elements manufactured by using additive technologies]. Siberian Journal of Science and Technology. 2019, Vol. 20, No. 2, P. 243-250. Doi: 10.31772/2587-6066-2019-20-2-243-250
Для цитирования: Ушакова Е. С. Моделирование напряженно-деформированного состояния конструкций ракетно-космической техники, изготовленных с использованием аддитивных технологий // Сибирский журнал науки и технологий. 2019. Т. 20, № 2. С. 243-250. Doi: 10.31772/2587-6066-2019-20-2-243-250_
MODELING OF THE STRESS-STRAIN STATE OF ROCKET-SPACE TECHNOLOGY STRUCTURAL ELEMENTS MANUFACTURED BY USING ADDITIVE TECHNOLOGIES
E. S. Ushakova
Bauman Moscow State Technical University 5, 2-ya Baumanskaya St., Moscow, 105005, Russian Federation E-mail: ellizaweta@gmail.com
One of the promising areas for improving the methods of manufacturing structural elements of rocket and space technology is the use of selective laser melting technology which represents a unique opportunity to manufacture metal products by melting powder and producing a one-piece solid phase structure. However, pores and other structural defects can appear in the formed element during laser sintering which causes a decrease in the strength characteristics of the parts produced. An important step in the additive technologies introduction is the development of methodology for the preliminary prediction of the strength characteristics of manufactured structural elements under the influence of mechanical loads with the help of mathematical modeling. The methodology for estimating the material strength reduction of a rocket-space technology element obtained using additive technologies by simulating a porous structure and calculating the characteristics of the stress-strain state is presented.
The proposed mathematical model and the methodology for calculating the specimen loading on the basis of the distortion energy theory allow calculating the stress-strain state in the process of numerical simulation for different values of the pore diameter. The reduction in yield strength due to the material porosity of the part is estimated using a coefficient equal to the ratio of equivalent stresses arising when a load is applied to a specimen manufactured using traditional and additive technologies. The value of the introduced coefficient characterizes the structure of the grown product and is considered as a function of the random arrangement of pores in the specimen under study. The appearance of pores is the result of a combination of factors: the composition and dispersion of the original metal powder, feed rate, removal distance and laser power during sintering, part orientation and sintering direction, the height of the level of powder deposited on a special base before sintering, etc.
The paper evaluates the reduction in strength for the working part of a series of tensile test specimens grown from metal powder of different dispersity. The non-linear nature of the dependence of the yield strength on the particle diameter of the original metal powder is established. The maximum value of the yield strength corresponds to the specimen with the minimum value of the total surface area of the pores.
Keywords: additive technologies, liquid rocket engine, combustion chamber, porosity, stress-strain state, yield strength.
МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ КОНСТРУКЦИЙ РАКЕТНО-КОСМИЧЕСКОЙ ТЕХНИКИ, ИЗГОТОВЛЕННЫХ С ИСПОЛЬЗОВАНИЕМ АДДИТИВНЫХ ТЕХНОЛОГИЙ
Е. С. Ушакова
Московский государственный технический университет имени Н. Э. Баумана Российская Федерация, 105005, г. Москва, 2-я Бауманская ул., 5 Е-mail: ellizaweta@gmail.com
Одним из перспективных направлений совершенствования методов изготовления конструктивных элементов ракетно-космической техники является применение технологии селективного лазерного плавления, которая дает уникальную возможность изготавливать изделия из металла посредством расплавления порошка и получения сплошной твердофазной структуры. Однако при лазерном спекании в формируемом элементе могут образовываться поры и прочие дефекты структуры, что вызывает снижение прочностных характеристик изготавливаемых деталей. Важным этапом при внедрении аддитивных технологий является разработка ме-
тодов предварительного прогнозирования прочностных характеристик изготавливаемых элементов конструкции в условиях воздействия механических нагрузок с помощь математического моделирования. Представлена методика оценки снижения прочности материала элемента конструкции ракетно-космической техники, полученного с использованием аддитивных технологий с помощью моделирования пористой структуры и расчета характеристик напряженно-деформированного состояния.
Предложенные математическая модель и методика расчета нагружения образца, основанная на теории энергии формоизменения, позволяют провести расчёт напряженно-деформированного состояния в процессе численного моделирования для разных значений диаметра пор. Снижение предела текучести в связи с пористостью материала детали оценивается с помощью коэффициента, равного отношению эквивалентных напряжений, возникающих при приложении нагрузки к образцу, изготовленному посредством традиционных и аддитивных технологий. Значение введенного коэффициента характеризует структуру выращенного изделия и рассматривается как функция случайного расположения пор в исследуемом образце, появление которых является результатом влияния совокупности факторов: состав и дисперсность исходного металлического порошка, скорость подачи, расстояние выведения и мощность лазера при спекании, ориентация детали и направление спекания, высота уровня порошка, нанесенного на специальное основание перед спеканием и др.
В работе проведена оценка снижения прочности для рабочей части серии образцов для испытаний на растяжение, выращенных из металлического порошка различной дисперсности. Установлен нелинейный характер зависимости предела текучести от диаметра частиц исходного металлического порошка. Максимальное значение предела текучести соответствует образцу с минимальным значением суммарной площади поверхности пор.
Ключевые слова: аддитивные технологии, жидкостной ракетный двигатель, камера сгорания, пористость, напряженно-деформированное состояние, предел текучести
Introduction. Currently, improving manufacturability and reducing the cost of manufacturing rocket-space technology (RST) structural elements are important scientific and technical challenges. The disadvantage of traditional methods of manufacturing the most heavily loaded units of a liquid-propellant rocket engine (LRE) - combustion chambers (CCs) - is the problem of ensuring the reliability of structures and controlling their quality due to the use of a large number of soldered joints in manufacturing as well as the need to manufacture expensive tooling. For example, only the process of electroerosive deposition of the "artificial" roughness creating turbulence in the coolant flow on the bottom of the cooling tract channels of the CC inner shell lasts 900 n / h [1].
One of the promising directions of improving the methods of manufacturing the rocket-space technology structural elements is the use of additive technologies (AT) [2]. Their main advantages include: manufacturing objects of complex shape with high accuracy, optimizing time spent on manufacturing, using a combination of metal powders (for example, BrH08 + stainless steel) in order to obtain strength characteristics corresponding to the operating conditions of the grown product [3].
The development of additive technologies, in particular, the technology of selective laser melting (SLM), represents a unique opportunity to manufacture metal products by melting powder and obtaining a one-piece solid phase structure [4]. The use of SLM technology allows increasing the material utilization rate to almost 99 % and, thereby, reducing the cost of production [5]. However, during SLM pores and other structural defects (non-melts, cracks, inclusions, etc.) may appear in the formed element, which causes a decrease in the strength characteristics of the parts produced. One of the objects of RST where the introduction of AT seems to be promising is the LRE chamber.
Since the working process in the LRE is characterized by relatively large values of pressure and temperature in
the CC [6], an important step in the introduction of AT is the development of methods for preliminary prediction of the strength characteristics of manufactured structural elements under the influence of mechanical loads using mathematical modeling.
The purpose of this work is to develop a methodology for evaluating the strength reduction of the material obtained using additive technologies by modeling the porous structure and calculating the characteristics of the stressstrain state (SSS).
Mathematical model. The methodology for assessing the strength reduction of parts manufactured using additive technologies on the basis of the SSS determination in the process of numerical simulation using the distortion energy theory (Mises-Huber-Genki theory) [7] has been developed. According to this methodology it is assumed that the specimen begins to break (or to become unac-ceptably deformed) under the condition aeq > oyield, where <5eq is von Mises equivalent stress arising in the part under the action of a given load, Gyield is the yield strength of the material of the part.
A safety factor - the ratio of yield strength to equivalent stress Gyield determined, for example, according to Mises [7] - can be used as a characteristic of the strength reliability of a RST product:
K.
yield
The reduction of the yield strength GyieU.add of the RST structural element material grown from metal powder is taken into account using the coefficient n (<1):
} yield.add n ^yield'
Then the safety factor for parts made using additive technology will be determined by the formula:
K,
yield .add
a
n • a
yield
(i)
a
eq
In this formulation, the value of the coefficient n characterizes the structure of the grown product and is considered as a function of the random arrangement of pores in the specimen under study. The appearance of pores is the result of the influence of a combination of factors: the composition and dispersion of the original metal powder, feed rate, removal distance and laser power during sintering, part orientation and sintering direction, the height of the level of the powder deposited on a special base before sintering, etc. [8-12].
Pores (as well as cracks, non-melts, etc.) in the structure of any material are local sharp changes in the uniformity of shape and rigidity of the product structure and lead to local increase in the value of internal stresses. Consequently, the value of the equivalent stress Geq.pore caused by a given load for a part with pores will be higher than the value for a monolithic part ceq made of material with the same physical and mechanical properties. The safety factors for the cases considered are determined by the formulas:
K,
.solid
K,
K.
u yield
CT„„
yield
(2)
K > K
s.solid pore '
By presenting a part grown from metal powder as a monolithic one with the presence of local inhomoge-neities in the material structure it is possible to consider the values of safety factors from formulas (1) and (2): Ks.add = Kspore to be equal. Then, in accordance with the idea of the physical meaning of the coefficient n, its value will be equal to the ratio of equivalent voltages: n -ct„
yield
yield
a
a
n = -
Consequently, the decrease in the tensile strength of a metal part made by the method of SLM is due to increase in the value of equivalent stresses as compared to a metal part made using traditional methods.
The porosity of the part material is determined by the formula:
S pore = 1 _ klay ,
where: klay= 1 —
Vp
V
is the packing factor, V is part
volume, Vpore is the pore volume in the part material.
Calculation model of the specimen loading. Structural cryogenic steel AISI 316L the mechanical properties [13] of which are listed in Table 1 is considered as a material of the structural RST element.
The porous structure of the specimen created by the method of selective laser sintering from ASIS 316L steel was studied in [14]. The relationship between the structure of the grown material and the sintering and the dispersion modes of the original metal powder was experimentally established in this work. A complex configuration of pores can be approximated by a set of spherical surfaces on the basis of a number of images obtained on the surface of thin sections of material specimens. In the first approximation we will consider the diameter and distance between pores commensurate with the diameter of a metal particle dp (fig. 1):
dp
'' 2 ' dp , hpore
-2 • dp
In this paper it is expedient to evaluate the reduction in strength for a specimen of a material and not for the design of a LRE chamber as a whole. For this purpose a series of specimens for tensile test were modeled according to GOST 1497-84 [15]. The dimensions the specimens correspond to type 3, number 9 (fig. 2).
The results of the stress-strain state calculations.
Assuming that the applied load is evenly distributed over the cross section the strength calculation was carried out directly for a cylindrical fragment of the specimen working part (h = 0.5 mm, V = 3.534 mm3, fig. 2) grown from AISI 316L powder of the following series of diameters: dpi = 60; 90; 120; 150; 180; 210; 240; 270; 300 microns (fig. 3). The calculated values of the packing factors klayi for each element are presented in tab. 2 A fragment of the working part was fixed on the surface of the rear end; a tensile load was applied along the normal to the surface of the front end. An example of a finite-element model of a calculation object with dp = 60 ^m in the form of a computational grid using tetrahedral elements with a total number of ~ 6 ■ 104 is shown in fig. 4. Verification of the calculated diagram of porous specimens tension (fig. 5) was carried out using an experimental diagram for steel [15]. Comparison of dependences shows satisfactory coincidence of the qualitative nature of the calculated and experimental diagrams for elastic and the beginning of plastic deformations.
Table 1
Steel AISI 316L mechanical properties at T = 20 oC
Density Yield strength Tensile strength Young's modulus Poisson's ratio
P kg m ° yield ,Mpa [Vtens ]' MPa E, GPa
7860 300 570 200 0'26
eq. pore
a
eq
a
eq. pore
а б
Fig. 1. Scheme of distribution (a) and geometrical characteristics (b) of pores in the object under study
Рис. 1. Схема распределения (а) и геометрические характеристики пор (b) в объекте исследования
Fig. 2. 3D-model of the object under study Рис. 2. 3D-модель объекта исследования
5
6
7
8
Fig. 3. 3-D models of porous specimens grown from metal powder with dpi: 60 ^m (1), 90 ^m (2), 150 ^m (3), 180 ^m (4), 210 |xm (5), 240 ^m (6), 270 ^m (7), 300 |xm (8)
Рис. 3. 3-D модели пористых образцов выращенных из металлического порошка с dpi: 60 мкм (1), 90 мкм (2), 150 мкм (3), 180 мкм (4), 210 мкм (5), 240 мкм (6), 270 мкм (7), 300 мкм (8)
Fig. 4. A finite-element model in the form of a structured grid of the calculation object at dp = 60 ^m
Рис. 4. Конечно-элементная модель в виде структурированной сетки объекта расчета при ¿р = 60 мкм
1
2
3
4
The value of the equivalent stress oeq.porei (tab. 2) arising from the application of tensile load equal to F = 2000 N is obtained as a result of a series of calculations for each specimen. The values of the coefficient ni and the yield strength Gyield.porei are presented in tab. 2. The dependence of the value of the yield strength on the particle diameter (hence, on the pore diameter) is nonlinear (fig. 6). The existence of an extremum can be explained by the opposite influence of geometrical parameters on the value of the surface area of the pores in the entire volume of the specimen. The total area of spherical surfaces is proportional to the number of pores in the
specimen and at the same time to the diameter of each pore dporei. The smaller the specified diameter of the initial metal powder particle is, the smaller the diameter and surface area of the pore formed in the SLM are, but the smaller the distance between the pores (hpore ~ 2-dpore) is, the greater their number is, and, therefore, the larger the total surface area of pores is.
The maximum value of the yield strength (cyield.p.max = = 291.9 MPa) corresponds to the specimen with the highest value of the packing coefficient (klay max = 0.956), that is, the specimen with the minimum total surface area of pores.
Fig. 5. The calculated diagram of the specimen stress-strain state with: 1 - dp = 0 |im; 2 - dp = 180 |im; 3 - dp = 300 |im
Рис. 5. Расчетная диаграмма напряженно-деформированного состояния образца при: 1 - dp = 0 мкм; 2 - dp = 180 мкм; 3 - dp = 300 мкм
Table 2
The results of the yield strength calculation for porous specimens
№ ^ eq. porei ni ® yield. porei klayi
mkm MPa - MPa
1 60 254.4 0.917 275.0 0.929
2 90 251.3 0.928 278.4 0.929
3 150 240.8 0.968 290.4 0.955
4 180 239.6 0.973 291.9 0.956
5 210 240.7 0.969 290.7 0.944
6 240 244.5 0.954 286.2 0.939
7 270 249.1 0.936 280.9 0.929
8 300 255.1 0.914 274.2 0.914
0 (monolithic) 233.2 1 300.0 1
265
О 50 100 150 200 250 300 350
dp, [Ш1
Fig. 6. The dependence of the specimen yield strength value on the metal powder particle diameter
Рис. 6. Зависимость значения предела текучести образца от диаметра частиц металлического порошка
Conclusion. The paper considers the methodology of evaluating the reduction of the strength of the RST structural element material obtained by the SLM method by modeling the porous structure and calculating the characteristics of the stress-strain state. During the calculation it was found that the value of the total surface area of the pores in the material due to the size of the original metal powder particles is directly proportional to the value of the equivalent stress aeq,porei arising in the object under study when tensile load is applied.
The dispersion value of the metal powder dp is determined using tensile testing of a series of metal specimens. The value determined provides the minimum total surface area of pores in the material structure resulting from laser beam sintering, thereby causing greater value of the yield strength Gyield.pore of the part. The optimum value of the powder dispersion for steel AISI 316L in the preparation of specimens for tensile is dp opt = 150...200 microns.
References
1. Artemov A. L. [Development of design and technology solutions for additive manufacturing of prototype inner lining for combustion chamber of multifunctional liquid-propellant rocket engine]. Kosmicheskaya tekhnika i tekhnologii. 2017, No. 1, P. 50-62 (In Russ.).
2. Abdullin M. I., Basyrov A. A., Nikolaev A. V. [Metal filled compositions for 3D printing]. Universum: khimiia i biologiia. 2015, No. 11(18). Available at: http://7universum.com/ru/nature/ archive/item/2701 (accessed 05 July 2017).
3. Grigoryanc A. G., Kazaryan M. A., Lyabin N. A. Lazernaya precizionnaya mikroobrabotka materialov [Laser Precision Micromachining]. Moscow, Fizmatlit Publ., 2017, 413 p.
4. Elistratova A. A., Korshakevich I. S., Tikhonenko D. V. [3D-printing technologies: advantages and disadvantages]. Reshetnevskie chteniia. Sbornik materialov
XVIII Mezhdunarodnoi nauchno-prakticheskoi konferent-sii [Reshetnev readings. The collection of materials of XVIII International scientific and practical conference]. Krasnoyarsk, 11-14 November 2014, SibSAU Publ., 2014, P. 557-559.
5. Zlenko M. A., Popovich A. A., Mutylina I. N. Additivnye tekhnologii v mashinostroenii [Additive technologies in mechanical engineering]. Saint-Petersburg, Izd-vo Politekhnicheskogo universiteta Publ., 2013, 223 p.
6. Mosolov S. V., Sidlerov D. A. [Comparative analysis of the characteristic features of operational process in liquid rocket engine combustion chamber with jet-centrifugal and cebtrifugal-centrifugal injectors by means of numerical simulation]. Vestnik MGTU im. N. E. Baumana. Ser. Mashinostroenie. 2016, No. 2, P. 60-71 (In Russ.).
7. Meiz Dzh. Teoriia i zadachi mekhaniki sploshnykh sred [Theory and problems of continuum mechanics]. Moscow, LKI Publ., 2007, 320 p.
8. Galimova L. A., Atroshchenko V. V., Smirnov V. V. et al. [Structure and properties of stainless steel specimens received by method of selective sintering]. Vestnik Bash-kirskogo universiteta: matematika i mekhanika. 2016, No. 2, P. 258-263 (In Russ.).
9. Hendrickson J. W. Use of Direct Metal Laser Sintering for Tooling in High Volume Production. USU Library, Logan, Utah, 2015. 35 p.
10. Torrado A. R., Roberson D. A. Failure Analysis and Anisotropy Evaluation of 3D-Printed Tensile Test Specimens of Different Geometries and Print Raster Patterns. Journal of Failure Analysis and Prevention. 2016, No. 1, P. 154-164.
11. Delgado J., Ciurana J., Rodriguez C. A. Influence of process parameters on part quality and mechanical properties for DMLS and SLM with iron-based materials. The International Journal of Advanced Manufacturing Technology. 2012, Vol. 60, No. 5-8, P. 601-610.
12. Smith D. H. et al. Microstructure and mechanical behavior of direct metal laser sintered Inconel alloy 718. Materials Characterization. 2016, Vol. 113, P. 1-9.
13. Features of Local Corrosion of AISI316L Steel Manufactured by Selective Laser Melting. Kuznetsov P. A., Krasikov A. V., Staritsyn M. V., Mushnikova S. Y., Parmenova O. N. Protection of Metals and Physical Chemistry of Surfaces. 2018, Vol. 54, No 3, P. 484-489.
14. Gu D., Shen Y. Processing conditions and microstructural features of porous 316L stainless steel components by DMLS. Applied Surface Science. 2008, Vol. 255, No. 5, P. 1880-1887.
15. GOST1497-84 Metally. Metody ispytanij na rast-yazhenie [State Standard 1479-84. Metals. Tensile test methods]. Moscow, Standartiform, 2005, 24 p.
Библиографические ссылки
1. Отработка конструктивных и технологических решений для изготовления опытных образцов внутренней оболочки камеры сгорания многофункционального жидкостного ракетного двигателя с использованием аддитивных технологий / Артемов А. Л. [и др.] // Космическая техника и технологии. 2017. № 1. С. 50-62.
2. Абдуллин М. И., Басыров А. А., Николаев А. В. Металлополимерные композиции для 3D печати // Universum: химия и биология. 2015. № 11(18). URL: http//7universum.com/ru/nature/archive/item/2701 (дата обращения 20 декабря 2018).
3. Григорьянц А. Г., Казарян М. А., Лябин Н. А. Лазерная прецизионная микрообработка материалов : монография. М. : Физматлит, 2017. 413 с.
4. Елистратова А. А., Коршакевич И. С., Тихоненко Д. В. Технологии 3D-печати: преимущества и недостатки // Решетневские чтения : сб. матер. XVIII Меж-дунар. науч.-практ. конф. (Красноярск, 11-14 ноября 2014 г.) / СибГАУ. Красноярск, 2014. С. 557-559.
5. Зленко М. А., Попович А. А., Мутылина И. Н. Аддитивные технологии в машиностроении. СПб. : Изд-во Политехн. ун-та, 2013. 223 с.
6. Мосолов С. В., Сидлеров Д. А. Анализ особенностей рабочего процесса в камерах сгорания ЖРД со
струйно-центробежными и центробежно-центробеж-ными форсунками // Вестник МГТУ им. Н. Э. Баумана. Сер. «Машиностроение». 2016. № 2. C. 60-71.
7. Мейз Д. Теория и задачи механики сплошных сред. М. : Изд-во ЛКИ, 2007. 320 с.
8. Галимова Л. А., Атрощенко В. В., Смирнов В. В., Чуракова А. А., Гундеров Д. В., Заманова Г. И. Структура и механические свойства образцов из нержавеющей стали, полученных методом селективного спекания // Вестник Башкирского ун-та: математика и механика. 2016. № 2. С. 258-263.
9. Hendrickson J. W. Use of Direct Metal Laser Sintering for Tooling in High Volume Production. USU Library, Logan, Utah, 2015. 35 p.
10. Torrado A. R., Roberson D. A. Failure Analysis and Anisotropy Evaluation of 3D-Printed Tensile Test Specimens of Different Geometries and Print Raster Patterns // Journal of Failure Analysis and Prevention. 2016. No. 1. Р. 154-164.
11. Delgado J., Ciurana J., Rodriguez C. A. Influence of process parameters on part quality and mechanical properties for DMLS and SLM with iron-based materials // The International Journal of Advanced Manufacturing Technology. 2012. Vol. 60, No. 5-8. P. 601-610.
12. Smith D. H. et al. Microstructure and mechanical behavior of direct metal laser sintered Inconel alloy 718 // Materials Characterization. 2016. Vol. 113. P. 1-9.
13. Features of Local Corrosion of AISI316L Steel Manufactured by Selective Laser Melting / P. A. Kuznetsov, A. V. Krasikov, M. V. Staritsyn [et al.] // Protection of Metals and Physical Chemistry of Surfaces. 2018. Vol. 54, No. 3. P. 484-489.
14. Gu D., Shen Y. Processing conditions and microstructural features of porous 316L stainless steel components by DMLS // Applied Surface Science. 2008. Vol. 255, No. 5. P. 1880-1887.
15. ГОСТ 1497-84 Металлы. Методы испытаний на растяжение. М. : Стандартиформ, 2005. 24 с.
© Ushakova E. S., 2019
Ushakova Elizaveta Sergeevna - is a technician of NIIEM and an employee of the Department of Rocket Engines; Bauman Moscow State Technical University. E-mail: ellizaweta@gmail.com.
Ушакова Елизавета Сергеевна - техник НИИЭМ, сотрудник кафедры ракетных двигателей, Московский государственный технический университет имени Н. Э. Баумана. E-mail: ellizaweta@gmail.com.
