Computer modeling of hot isostatic pressing process of porous blank
Автор: Barkov L.A., Cheskidov P.A., Samodurova M.N., Latfulina Yu.S.
Журнал: Вестник Южно-Уральского государственного университета. Серия: Металлургия @vestnik-susu-metallurgy
Рубрика: Порошковая металлургия и композиционные материалы
Статья в выпуске: 3 т.16, 2016 года.
Бесплатный доступ
Existing mathematical models for the process of HIP porous blanks are based on the solution of approximately differential equations of the equilibrium of quasi-continuous medium. The continual approach is used for computer modeling of the compressible body HIP process. Partial differential equations of motion for a quasi-continuous packing medium and physical equations for a viscous-plastic isotropic porous material subjected to work hardening are taken as the basis for the mathematical model by means of these simultaneous equations. Besides the equations of motion and the rheological equation, the equation of continuity deformation and the equation of heat flow are used. Numerical calculations for additional packing of a high-speed steel powder billet preliminary pressed in the hydrostat are performed. Numerical calculation of the problem for hardening the cylindrical high-speed steel billet with a mild steel shell is done using Lagrange’s method by means of the difference scheme of continuous calculation of the Wilkins’ type. Computer modeling allows to control the process of hardening and changing the form of a porous body during the HIP process.
Continuum, mathematical model, powder blank, different equation, rheology
Короткий адрес: https://sciup.org/147157033
IDR: 147157033 | DOI: 10.14529/met160318
Текст научной статьи Computer modeling of hot isostatic pressing process of porous blank
Hot isostatic pressing (HIP) is one of the most suitable techniques for producing the high-speed steel cemented carbide, superalloys, soft ferrites and composites at a simultaneous application of high-pressure and high temperature [1–3]. The governing of this advanced process is performed, as a rule, by computer. Therefore, development of more perfect mathematical models of the process of plastic deformation is the most important task of further improving the HIP technology and increasing the quality of new materials.
Two basic approaches to the construction of the mathematical model adequately describing the HIP process are generally used. The first method, a discrete one, is developed in the works [4–8]. Here, they consider different mechanisms of packing at HIP such, as plastic deformation, creeping, and diffusion described with kinetic equations. In the article [8] a porous material is modeled as a sphere with the spherical pore subjected to creeping deformation under isostatic loading conditions. However, introduction of the central symmetry considerably delays the packing process. It should be noted that there is a somewhat difficulty in application of the first method when modeling the HIP of workpieces having a complex configuration and shells which distort the uniformity of a stress-strain state of a porous body.
The second method, a continual one, is developed in [9–17]. Here, a porous body is considered as a complex packed medium with an irreversible volumetric strain. This allows to solve multidimensional problems of elasto-viscous-plastic flow of a porous body and a shell. In the work [17] the analyses of HIP technological modes based on the application of both the first and the second methods are given.
We have developed the continual approach. Unlike the works of other authors, the present investigation deals with the differential equations of motion introduced into the mathematical model instead of equilibrium equations. This allows to employ Lagrange’s method in numerical realization of the model by computer. The equations, in a special case, are written in the cylindrical coordinate system.
1. Basic equations for computer modeling
The HIP process, in a general case, is a threedimensional non-stationary non-isothermal problem. Following the works [15–17] we have constructed simultaneous equations of motion describing the flow of a viscous-plastic isotropic compressible material subjected to strain hardening in the right-hand rectangular coordinate system x i . Denote in the equations: p is the density, t is the temperature, т is the time, v i is the component of the flow velocity vector, σ i j , s i j , ξ i j are the components of stress tensors, a stress deviator, and strain rates. The basic equations have a form:
-
– the equation of motion ^^
-
» ц:=сил;(1)
-
– the equation of continuity
^+ри = 0;(2)
-
p dx uv 7
-
– the equation of heat conductivity
-
— сP dl = (A grad г)(а + VOij^ij,(3)
where c is the specific heat, λ = λ(t, ρ) is the coeffi- cient of heat conductivity. The differentiation according to time is done along the path of motion of a material particle.
T h e c om pon en t s o f th e v e l oc ity v e ct or a n d t h ose o f t h e defor m a ti on rate a re de t e r m i n e d w i th t h e kn o w n relations:
idp + du ? + dv - _ p dx dz dr
-
dV r ds— । d^ rz । P dx dr dz
v- r ’
Srr-See + 3a r dr ’
dV y ds -z | ^^ zz | ^ rz | da dx dr dz r dz’
dx, ;
^lj = - (v t,t + v i,i ). (5)
The stress tensor, as usual, con sists of a sph e ri c al constituent and a deviator:
a tj = a& tj + s tj ,
drdz
Tx_V r , Tx_V z ,
Srr _ 2^(dVL + ^^) rr dr3p dx
^_2Ц(^ + ^^).
zz dz3p dx
Srz_^(dV : + dV z ), rz dz ^r ’
,
,
where σ i s t he me a n n or m al stre ss, a n d δ ij is Kronecker’s symbol. The equations (1)– ( 6) a re c om pl e t e d w i th the rheological equations:
S tj = 2vGlJ-;W (7)
a = 3k,
v 1 dp. s ee _2^(7- + ^pd; ),
_ к dp p dx ’
T 2 = 2 2 -^ 4- ?2
i s rr s zz s ee s rz ,
F(T,a,as,q t ) _ 0.
where ξ i s th e ra t e of re l a t iv e c h a ng e o f volume, μ, k are coe f f icie n ts of sh ea r an d v ol ume t ri c v i sc os i ty de p e n de n t u pon de nsi ty a n d t em pe ra t u re .
The l o ad in g sur f ace, i n ge ner al cas e, ha s t he fo rm :
F (T,o,os,q i ,^,qn) _ 0, (9)
The behavior of the metallic shell of the porous cylinder corresponds to the elasto-plastic flow of the material [18]. The differential equations have to satisfy the following initial and boundary conditions:
where T = ^st j St j i s th e i n t en s i ty o f t a ng en t i a l stresses, σ s is t h e y i e ld s t ress o f a sol i d p h a se o f t h e ma terial a n d i t is th e te m pera t u re f un c t i on , q i are work hardening parameters.
2. Mathematical model of HIP process
C on sider th e ma th e m ati c a l mode l of t h e HIP is o t h e rma l proce ss for a cyli n d ri c a l porou s body w i th the radius R and the height H w i th in a s h ell . Th e body su r f ac e i s su bje c t e d t o t h e pre ss u re of gas p = pG (τ) at a st a ti on a ry te m pera t u re duri n g th e pre ssin g cy c l e . The simultaneous equations (1) – (9) describing the viscous-pla sti c fl o w o f a cyli n drica l porou s c o mpressible body is written as follows:
at t = 0, p = p o , Sy = 0, о = 0, v= 0;
at t ≥ 0 on the surface of the shell
σ = – pG (τ), srr = szz = sθθ = srz = 0.
Within the boundary of the porous body and the shell the condition of continuity of the velocity vector ῡ is assigned.
Solution of the equations is done by means of different schemes of continuous calculations of Wilkins’ type [18].
3. Results of numerical calculations
In numerical calculations the HIP process of a porous compressible high speed steel billet pressed in the hydrostat was modeled. The powder billet within

Fig. 1. Dynamics of a cylindrical porous body at HIP: p G = 0–200 MPa, at a velocity of the pressure increase w , equal to 10 MPa/s

Fig. 2. Dynamics of a porous body at HIP: p G = 0–200 MPa, w = 500 MPa/s, ρ0 = 0,7 ρ c
Порошковая металлургия и композиционные материалы
the mild steel shell was subjected to heating and then to hot gasostatic pressing.
The initial density of the porous body at HIP is assumed to be equal to 0.93ρ c , where ρ c is the density of non-porous high speed steel, and the process temperature being as high as 1300, 1400 and 1500 °C. The dependence pG (τ) is taken as a piecewise linear one. The stress surface is assigned:
| г2+Д(р)а2 = ^(рШО^)2, (11)
where p = — is the relative density,
P c
/ i (p) = Q i (1 — p) ”1 ,
A®^”2, f (0 = ®® + 1).
In Figs. 1 and 2 there are shown the cylindrical porous body at the initial instant (1) and final instant (2) in the HIP process as well as the difference grid on which the calculation is done. The coefficient values are taken from the works [7, 11].
Conclusion
The HIP process of a porous compressible body is described with the mathematical model based on differential equations of motion for a quasi-conti-nuous medium having an irreversible volumetric and shear deformation. The equations for a viscous- plastic isotropic porous hardening material are used as rheological ones. For the flows with an axial symmetry of the billet material after CIP the constitutive equations are written in the cylindrical coordinates. Numerical calculation of the problem for hardening the cylindrical high-speed steel billet within a mild steel shell is done using Lagrange’s method by means of the difference scheme of continuous calculation of Wilkins’ type.
Computer modeling allows to control the process of hardening and changing the form of a porous body at HIP in the course of time.
Список литературы Computer modeling of hot isostatic pressing process of porous blank
- Advanced Technology for New Materials HIP, Brochure No. M3C2010. -Kobe Steel Ltd., Japan, 1969. -21 р.
- Quintus Hot Isostatic Presses. Pamphlet A08 -400 8E. -ABB Asea Brown Boveri, Sweden, 1971. -16 р.
- Your Partner in High Pressure Technology National Forge. -National Forge, USA, 1972. -32 p.
- Arzt, E. Practical Applications of Hot-Isostatic Pressing Diagrams: Four Case Studies/E. Arzt, M.F. Ashby//MetallurgicalTransactions A. -1983. -Vol. 14A, no. 1. -P. 211-221 DOI: 10.1007/BF02651618
- Helle, A.S. Hot-Isostatic Pressing Diagrams: New Developments/A.S. Helle, K.E. Easterling, M.F. Ashby//Acta Metallurgica. -1985. -Vol. 33, no. 12. -P. 2163-2174 DOI: 10.1016/0001-6160(85)90177-4
- Li, W.-B. On Densification and Shape Change During Hot Isostatic Pressing/W.-B. Li, M.P. Ashby//Acta Metallurgica. -1987. -Vol. 35, no. 12. -Р. 2831-2842 DOI: 10.1016/0001-6160(87)90282-3
- Лаптев, А.М. Параметры горячего изостатического прессования пористых материалов/А.М. Лаптев, В.Н. Самаров, С.В. Подлесный//Известия АН СССР. Металлы. -1988. -№ 5. -С. 92-99.
- Duva, J.M. The Densification of Powders by Power-Law Creep During Hot Isostatic Pressing/J.M. Duva, P.D. Crow//Acta Metallurgica et Materialia. -1992. -Vol. 40, no. 1. -P. 31-35 DOI: 10.1016/0956-7151(92)90196-L
- Numerical Simulation of the Densification Behaviour of Metal Powder During Hot Isostatic Pressing/A. Nohara, T. Nakagawa, T. Soh, T. Shinke//International Journal for Numerical Methods in Engineering. -1988. -No. 25. -P. 213-225 DOI: 10.1002/nme.1620250117
- Nohara, A. Computer Simulation of Densification Behaviour of Metal Powder During HIP/A. Nohara, T. Nakagawa, T. Soh//Bulletin of the Japan Institute of Metals. -1989. -Vol. 28, no. 11. -P. 922-926 DOI: 10.2320/materia1962.28.922
- Shima, S. Plasticity Theory for Porous Metals/S. Shima, M. Oyane//International Journal of Mechanical Sciences. -1976. -Vol. 18, no. 6. -P. 285-291 DOI: 10.1016/0020-7403(76)90030-8
- Abouaf, M. Elasto-Viscous-Plastic Flow of a Porous Body/M. Abouaf//Mémoires et Études Scientifiques de la Revue de Métallurgie. -1986. -Vol. 83, no. 6. -P. 325-328.
- Nomoto, T. Analytical Studies on Sintering Deformation Behavior During Hot Isostatic Pressing/T. Nomoto, H. Yoshimura//Quarterly Journal of the Japan Welding Society. -1990. -Vol. 8, no. 1. -P. 32-38 DOI: 10.2207/qjjws.8.32
- Nakagawa, T. Simulation of Cold and Hot Isostatic Pressing/T. Nakagawa, A. Nohara//Kobelko Technology Review. -1991. -No. 12. -P. 10-14 DOI: 10.1016/0026-0657(92)91574-4
- Скороход, В.В. Теория нелинейно-вязкого и пластического поведения пористых материалов/В.В. Скороход, М.Б. Штерн, И.Ф. Мартынова//Порошковая металлургия. -1987. -№ 8. -С. 23-30.
- Лаптев, А.М. Расчёт процессов уплотнения пористых материалов на основе теории ползучести/А.М. Лаптев//Технология лёгких сплавов. -1982. -№ 7. -С. 16-20.
- Фролов, А.А. Система проектирования технологических режимов горячего изостатического прессования порошковых материалов. 1. Постановка задачи. Описание модели верхнего уровня/А.А. Фролов, О.Б. Садыхов, Г.Я. Гун//Порошковая металлургия. -1991. -№ 6. -С. 12-17.
- Уилкинс, М.Л. Расчёт упруго-пластических течений: пер. с англ./М.Л. Уилкинс//Вычислительные методы в гидродинамике/под ред. Б. Олдера, С. Фернбаха, М. Ротенберга. -М.: Мир, 1967. -С. 212-263.