New Die-Compaction Equations for Powders as a Result of Known Equations Correction: Part 1–Review and Analysis of Various Die-Compaction Equations

: The well-known equations for the powder compaction process (PCP) in a rigid die published from the beginning of the last century until today were considered in this review. Most of the considered equations are converted into the dependences of densification pressure on the powder’s relative density. The equations were analyzed and their ability to describe PCP was assessed by defining the coefficient of determination when approximating experimental data on the compaction of various powders. It was shown that most of the equations contain two constants the values of which are determined by fitting the mathematical dependence to the experimental curve. Such equations are able to describe PCP with high accuracy for the compaction of powders up to a relative density of 0.9–0.95. It was also shown that different equations can describe PCP in the density range from the initial density to 0.9 with the same high accuracy, but when the process of compaction is extrapolated to higher values of density, the curves diverge. This indicates the importance of equations that can unambiguously describe PCP to a relative density equal to or close to 1.0. For an adequate description of PCP for relative density greater than 0.95, equations containing three or four constants have proven useful.


Introduction
The powder compaction process (PCP) in a rigid die is one of the main processes in final product fabrication from various powders.Thus, the comprehensive study and description of PCP has attracted the attention of many specialists in the field of powder metallurgy.There are many equations that describe PCP, and various approaches have been used to derive them.Initially, to describe PCP, simple mathematical functions were used, e.g., exponential [1][2][3], power [4][5][6], and logarithmic [7][8][9][10][11].Later the proposed equations took into account various physical phenomena occurring in PCP such as contact interaction of particles, local and general deformation, hardening, shear between particles and within particles, friction, etc. [12][13][14][15][16][17].To derive the equations with physical constants, some researchers considered a powder body as a combination of a large number of individual particles contacting each other [14,18,19].An approach based on the evolution of the contact interaction between individual particles under the effect of pressure can be designated as a discrete PCP theory.There is also an approach that considers a powder sample as a quasi-continuous two-phase body that can decrease in volume under the effect of pressure.In this case, the compaction of the quasi-continuous body in a rigid die makes it possible to obtain the corresponding equations for describing PCP [20][21][22][23][24][25].This approach can be considered a continuous one.
However, the equations for PCP in a rigid die are not able to describe the entire compaction process with high accuracy, since the process itself is complex and multi-stage.The known equations are able to adequately describe only one or two of the three available Powders 2024, 3 112 stages [18].In addition, there are fundamental discrepancies between the equations, which are associated with the description of PCP upon approaching a pore-free state, where some equations give a finite value of the compaction pressure while other ones give infinite pressure.It is therefore difficult to choose an equation that allows one to determine the compaction pressure accurately enough to achieve a relative density close to 1.0.
Certain difficulties are associated with the determination of the true plastic deformation of particles during PCP in a rigid die because of particle rearrangement that can take place and that depends on the shape of particles and on their yield strength.Many equations for PCP do not explicitly reflect the effect of particle rearrangement on the increase in sample density.Some authors [17,26] have made an attempt to take into account the degree of density change due to particle rearrangement.These authors believe that rearrangement takes place throughout almost the entire compaction process.However, according to other researchers [27], rearrangement ends at the initial stage of the compaction process.This implies the existence of a porosity threshold beyond which no rearrangement occurs.The presence of a threshold or critical phenomena during PCP in a rigid die can be described based on percolation theory [28][29][30], but the modification of the powder compaction equations using this theory has not yet led to adequate equations.
Another issue occurs when describing the compaction of "high" samples, i.e., samples with a height-to-diameter ratio exceeding 1.0.It refers to the non-uniform distribution of density along both the height and the diameter of the sample.To solve the problem of determining the density of a sample at its different points the finite element method known for metal forming is used.This method allows one to determine the density at any point of a sample with known parameters of powder particle strength, using the hardening law and the friction forces between the particles and the die wall [31][32][33].However, the task of determining the true strength of particles and the law of hardening during their deformation is not easy.
To determine the properties of a specific powder, experiments are required on triaxial powder compression and on the determination of the yield surface depending on the sample density.Such experiments are relatively complex and must be carried out in special laboratories [34][35][36], which makes it difficult to obtain quick results for different powders.In addition, information about the properties of a powder can also be obtained from experiments on compaction in a rigid die of "low" samples in which, after compaction, the height becomes half the diameter.In these samples the effect of friction forces can be neglected.It should be noted that most of the experiments are performed with low samples, and that compaction equations are tested mainly on such samples.Therefore, before solving the problems of determining the density distribution in high samples (height/diameter ≥ 1) and samples of complex shape (stepped), it is necessary to have an equation that allows one to describe the compaction process of low samples with the highest accuracy and with a change in density from initial density to close to 1.0.
The adequacy of an equation for PCP in rigid die is determined by assessing the accuracy of fitting the equation to the experimental data, that is, by finding the determination coefficient R 2 .However, in many cases the adequacy of equations is unknown.To confirm this thesis, two review papers on the equations of PCP in a rigid die, published in 2007 and 2017, can be noted.The author of one of them is T. Çomo glu (2007) [37], who considered the equations proposed by Walker, Bal'shin, Heckel, Kawakita and Lüdde, Cooper and Eaton, Leuenberger, Shapiro, Sonnergaard.In that review, it was noted that the PCP in a rigid die typically consists of three stages and it is difficult or impossible to describe it by a single equation.The author analyzes known equations with various constant parameters and tries to relate these parameters to the particular stage of compaction and to the mechanical properties of materials.He also contrasts or compares some equations and shows their similarities and differences.An attempt was made to evaluate the equations for their ability to describe the compaction process of plastic and ceramic powders, and a possibility was also considered of describing the compaction of powder mixtures consisting of plastic and brittle powders.However, in that review there is no graphical representation of the experimental and theoretical curves of the compaction process, and there is also no quantitative assessment of the accuracy of the powder compaction process description by different equations.
In the second review, by Popescu and Vidu (2017) [38], the equations of Shapiro-Kolthoff and Konopicky, Bal'shin, Heckel, Cooper and Eaton, Kawakita and Lüdde, Ge Rong-de, Panelli and Filho, Parilak and Dudrova, Castagnet and Leal Neto, Gerdemann and Jablonski were considered.That review compares different equations describing the compaction process of powder mixtures consisting of plastic and non-plastic powders, and also establishes, in some cases, the influence of the shape and size of brittle particles on the constants of equations.However, in that review, as in the previous one, there are no graphical examples (with the exception of one figure ) showing the correspondence of theoretical curves to experimental data, and there is also no quantitative assessment of the agreement degree.This does not make it possible to understand how a wider range of density changes affects the accuracy of the PCP description, and also does not allow us to evaluate the compaction process at the stage of powder compression in the region of relative density close to unity.
Therefore, the purpose of this communication is to review the known equations for PCP in a die.An attention will be focused on assessing the accuracy of PCP description by various equations, as well as on a graphical representation (visualization) of correspondence theoretical curves to experimental data.In this case, it seems to be appropriate to convert the known equations that show the dependence of relative or absolute density on the compacting pressure into the dependence of pressure on relative density, since such dependence can have a simpler form and thus makes it easy to determine the specific compaction energy of the powder using the formula: where w is the specific energy in J/m 3 or in Mpa, p(ρ) is the experimental or theoretical external pressure (on the punch) dependence on the relative density, and ρ 0 and ρ are the initial and current relative densities of the sample, respectively.The assessment of accuracy of PCP description by various equations was performed by computer approximation with these equations of the experimental data on compaction different powders.For approximation, the Russified program "Wolfram Mathematica 10.4" was used.

Accuracy Assessment of PCP Description Using Equations Obtained by Selecting Mathematical Functions
Equations in the form of simple mathematical functions were among the first formulated ones, and they appeared at the beginning of the 20th century.Most researchers point out that the first powder compaction equation was proposed by E. Walker (1923) [1].Then, almost the same equation was proposed by M. Balshin (1948) [2] and H. Lipson (1950) [3].If we take the compaction pressure as a function and the relative density as an argument, then the compaction equation for the powders, by these authors, will be an exponential function: where p is the compaction pressure, β = 1/ρ is the relative volume of the powder sample, ρ is the relative density, k 1 and k 2 are the constants, and a 1 = Exp(k 1 ).
Other researchers, e.g., W. Rutkowski and H. Rutkowska (1949) [4], C. Agte and M. Petrdlik (1951) [5], and G. Meyerson (1962) [6] used a power function to describe PCP in the form of the following relationship between compaction pressure and relative density: Powders 2024, 3 114 where a 2 and b are constant parameters determined by approximating experimental data on powder compaction.The presented mathematical functions (exponential (2) and power (3)) are quite simple, but there is no clear understanding what range of density change can be accurately described by them.It is only clear that the smaller this range, the more accurately it is described by these equations.The accuracy of approximation of real experimental data on various PCPs by these exponential and power dependences is shown by approximation of experimental data on the compaction of iron, copper, and nickel powders from publications [39][40][41], respectively.The experimental data on the compaction of these powders are given in Table 1.This makes it possible to use them to checkup the approximation carried out and approximate with equations by other researchers.The results of the approximation of experimental data by exponential (2) and power (3) equations are shown in Figure 1, and the values of the constants of these equations and the coefficient of determination R 2 are presented in Table 2.As seen, the exponential (2) and power (3) dependences describe PCP quite accurately from the initial relative density to a density of 0.95, and power dependence (3) describes PCP with a higher accuracy that was also noted in another research [6].
where  2 and b are constant parameters determined by approximating experimental data on powder compaction.The presented mathematical functions (exponential (2) and power (3)) are quite simple, but there is no clear understanding what range of density change can be accurately described by them.It is only clear that the smaller this range, the more accurately it will be is described by these equations.The accuracy of approximation of real experimental data on various PCPs by these exponential and power dependences is shown by approximation of experimental data on the compaction of iron, copper, and nickel powders from publications [39][40][41], respectively.The experimental data on the compaction of these powders are given in Table 1.This makes it possible to use them to checkup the approximation carried out and approximate with equations by other researchers.The results of the approximation of experimental data by exponential (2) and power (3) equations are shown in Figure 1, and the values of the constants of these equations and the coefficient of determination R 2 are presented in Table 2.As seen, the exponential (2) and power (3) dependences describe PCP quite accurately from the initial relative density to a density of 0.95, and power dependence (3) describes PCP with a higher accuracy that was also noted in another research [6].In the above equations, PCP is conceded to be a monotonous process without any indication that, in reality, it proceeds in several stages with their own mechanisms of powder densification.Most experts agree on the existence of three stages of PCP, although according to some authors, there are 4 stages [42,43].The three stages are characterized as follows.The first stage is the rearrangement of particles when elastic and slight plastic deformation can occur.The second stage is the compaction due to the local plastic deformation of the particles during their contact.The third stage is interpreted ambiguously: it is either the localization of plastic deformation near the pore or the general plastic deformation of particles.These stages are clearly manifested in case of constructing a graphical dependence of the experimental relative density on the densification pressure in logarithmic coordinates [44].Moreover, each stage can be described, as the authors [44] suggested, by a power equation with its own constants: where m is a constant different for each stage, ρ* and p* are the minimum relative density and pressure, respectively, for a particular compaction stage.We note that the PCP staging was also established when considering one of the most common equations for PCP in a rigid die.This equation was proposed by different authors, e.g., L.F. Athy (1930) [7], I. Shapiro and I. M. Kolthoff (1947) [8], K. Konopicky (1948) [9], T. N. Znatokova and V. I. Likhtman (1954) [10], and R.W. Heckel (1961) [11], and in the original record it has the form: where k and a are constants.This equation is often referred to as the Konopicky's equation or as the Heckel's equation.If the experimental data on the compaction of powders are presented in coordinates "ln(1/1ρ)" and "p", then a broken line is obtained, indicating different stages of compaction with their own values k and a.The Equation (5) often met in literature [45,46], and after transformation, where the pressure is a function and the relative density is an argument, It takes the following form: where It is important to note that the Equation ( 6) allows one to describe only the intermediate or second stage of powder compaction.Therefore, in most cases, when experimental data on powder compaction are limited by low density (ρ ≤ 0.9), the Equation (6) describes the compaction process without an initial stage.It should also be noted that Equation ( 6) can be written otherwise if the value of the constant a is determined by substituting pressure equal to zero (p = 0) into this equation.In this case, Equation (6) will look as: Here, the parameter ρ ′ 0 means the conditional initial density that is greater which the actual initial density of the powder ρ 0 .This indicates that Equations ( 6) and ( 7) cannot take into account the initial process of powder compaction.
Furthermore, Equation ( 5) has been modernized several times.In one case, modernization was proposed by M. Kuntz and H. Leuenberger [30].They transformed the Heckel's equation, i.e., Equation (6), to the form: where C is a constant, ρ c is the conditional initial density of the powder or the second constant.
The use of Equation ( 8) for approximation of the experimental process of various powders compaction showed that it was more adequate than the Equation ( 6).In particular, the experimental data on the compaction of iron, copper, and nickel powders presented in Table 1 were used for approximation by Equations ( 6) and (8).The results of approximation are shown in Figure 2, and the values of the equations constants and the coefficient of determination R 2 are presented in Table 3.In addition, Table 3 shows the value of the conditional initial density ρ ′ 0 , which is used in Equation ( 7) and is also determined by fitting the curve to the experimental data.
Powders 2024, 3, FOR PEER REVIEW 6 can be written otherwise if the value of the constant  is determined by substituting pressure equal to zero (p = 0) into this equation.In this case, Equation (6) will look as: Here, the parameter  0 ′ means the conditional initial density that is greater which the actual initial density of the powder ρ0.This indicates that Equations ( 6) and ( 7) cannot take into account the initial process of powder compaction.
Furthermore, Equation ( 5) has been modernized several times.In one case, modernization was proposed by M. Kuntz and H. Leuenberger [30].They transformed the Heckel's equation, i.e., Equation (6), to the form: where C is a constant, ρc is the conditional initial density of the powder or the second constant.
The use of Equation ( 8) for approximation of the experimental process of various powders compaction showed that it was more adequate than the Equation ( 6).In particular, the experimental data on the compaction of iron, copper, and nickel powders presented in Table 1 were used for approximation by Equations ( 6) and ( 8).The results of approximation are shown in Figure 2, and the values of the equations constants and the coefficient of determination R 2 are presented in Table 3.In addition, Table 3 shows the value of the conditional initial density  0 ′ , which is used in Equation ( 7) and is also determined by fitting the curve to the experimental data.Figure 2 and Table 3 show that the new Equation ( 8) describes PCP with higher accuracy than the original Heckel's equation (6).In this case, the value of the conditional initial density ρc in Equation ( 8), obtained as a result of the approximation, turns out to be lower than the actual initial density ρ0 that makes it possible to accurately describe the initial stage of the compaction process.Table 3. Constants and coefficient of determination R 2 obtained via approximation of experimental data on PCP for iron, copper, and nickel using Equations ( 6) and ( 8).   3 show that the new Equation ( 8) describes PCP with higher accuracy than the original Heckel's equation (6).In this case, the value of the conditional initial density ρ c in Equation ( 8), obtained as a result of the approximation, turns out to be lower than the actual initial density ρ 0 that makes it possible to accurately describe the initial stage of the compaction process.

Parameter
A higher accuracy of approximation of experimental data on PCP compared to the equation of K. Konopicky or R. W. Heckel can be reached by another modernization of the equation, which leads to the appearance of the so-called double logarithmic form of the equation for powder compaction in a rigid die.Three different teams of authors proposed such an equation independently.First, the work of Ge Rong-de (1991) [47] should be noted, in which he proposed a new differential equation for PCP in a rigid die: where K, n, and m are constants.At the same time, Ge Rong-de claimed that at n = 0, the integration of Equation ( 9) leads to an equation that describes PCP at low and high pressures with high accuracy.This equation in the original record is as follows: where A and B are constants.
In another case, the same equation was proposed by other authors-A.B. Yu and Z. H. Gu (1993) [48], who used their own differential equation in the form: where k, a, b are constants.Upon integrating Equation (11), these authors obtained an equation for PCP in the form: where K = k/(1 + b) and n = 1 + b are new constants.In the third case, L. Parilak and E. Dudrova (1994) [49] proposed an equation that in its original form was written as where θ and θ 0 are the current porosity and the initial porosity, respectively, of the powder in the die, p is the compaction pressure, K and n are constants.Equation ( 13), after making double logarithm and replacement of porosity by relative density, takes the form of Equation ( 12) into which Equation (10) can also be easily converted.Consequently, the analysis of PCP by at least three groups of authors led to the derivation of the same equation, which contains the powder initial density ρ 0 , i.e., the powder density after filling the die, and two constants that can be determined by fitting the curve of Equation (12) to experimental data on PCP.If in the Equation ( 12) we interchange the function and the argument, i.e., we take the pressure as a function and the relative density as an argument, then we obtain a simpler equation: where m = 1/n and B = (1/K) m are new constants.
Powders 2024, 3 It is of interest to check the accuracy of PCP description for different powders by Equation ( 12), or rather by Equation (14).Prior to this, we need to remember one more equation which also contains the powder initial density and two constants, and which was proposed much earlier than the Equation (12).We are talking about the equation of K. Kawakita and K.H. Ludde (1970) [50].That equation is rather often referred in the literature, and in original form it is as: where C is the degree of volume reduction, V 0 is the initial volume of the powder, V is the volume of the powder under pressure P, a and b are constants characterizing the powder.
If we transform it in such a way that the compacting pressure p is a function, and the relative density ρ is an argument, then it takes the form: where ρ 0 is the initial relative density of the powder filled into the die, a and b ′ = 1/b are constants.
To check the approximation accuracy of the experimental data on PCP by the equation of K. Kawakita and K.H. Lüdde ( 16) and the equation of Ge Rong-de, A. B. Yu and Z. H. Gu, L. Parilak and E. Dudrova ( 14), we used the experimental data on the compaction of powders of iron, copper and nickel presented in Table 1.The results of approximation are presented in the form of graphs in Figure 3, and the values of the constants of Equations ( 16) and ( 14) and the coefficient of determination R 2 are shown in Table 4.
the curve of Equation (12) to experimental data on PCP.If in the Equation ( 12) we interchange the function and the argument, i.e., we take the pressure as a function and the relative density as an argument, then we obtain a simpler equation: where m = 1/n and B = (1/K) m are new constants.It is of interest to check the accuracy of PCP description for different powders by Equation ( 12), or rather by Equation (14).Prior to this, we need to remember one more equation which also contains the powder initial density and two constants, and which was proposed much earlier than the Equation (12).We are talking about the equation of K. Kawakita and K.H. Ludde (1970) [50].That equation is rather often referred in the literature, and in original form it is as: where C is the degree of volume reduction, V0 is the initial volume of the powder, V is the volume of the powder under pressure P, a and b are constants characterizing the powder.
If we transform it in such a way that the compacting pressure p is a function, and the relative density ρ is an argument, then it takes the form: where ρ0 is the initial relative density of the powder filled into the die, a and b′ = 1/b are constants.
To check the approximation accuracy of the experimental data on PCP by the equation of K. Kawakita and K.H. Lüdde ( 16) and the equation of Ge Rong-de, A. B. Yu and Z. H. Gu, L. Parilak and E. Dudrova (14), we used the experimental data on the compaction of powders of iron, copper and nickel presented in Table 1.The results of approximation are presented in the form of graphs in Figure 3, and the values of the constants of Equations ( 16) and ( 14) and the coefficient of determination R 2 are shown in Table 4.  Table 4. Constants and coefficient of determination R 2 obtained via approximation of experimental data on the compaction of powders of iron, copper, and nickel using Equations ( 14) and ( 16).Accuracy assessment of the experimental data approximation on the compaction of iron, copper and nickel powders by Equations ( 14) and (16) show that these equations describe PCP very accurately and practically in the same way, since the coefficient of determination differs slightly.In this case, the coefficient of determination by Equation ( 14) for copper and nickel is slightly higher than that by Equation ( 16), and for iron it is slightly lower (Table 3).In addition, these equations take into account the initial density of the powder and accurately describe the initial stage of powder compaction.But the fundamental difference between them consists in the fact that at a relative density ρ = 1, the compaction pressure by Equation ( 16) has a finite value, while by Equation ( 14), it is equal to infinity.It is still difficult to say, which equation is more adequate, but the higher accuracy of Equation ( 14) in two cases out of three may indicate that the equation in which the pressure tends to infinity as the relative density increases up to 1.0 seems to be more valid.

Parameter
The above comparison of two equations shows how important it is to accurately describe the entire PCP, including the process of compaction to a relative density close to 1.0.Alongside with this, it should be noted that there are very few experimental data in the literature on the compaction of metal powders with a final density that exceeds 0.95.Such experimental data, in particular, on the compaction of various iron powders, are available in the book by R. Kieffer and W. Hotop (1948) [51].They also showed that it is difficult to achieve a relative density ρ = 1.0, even at a pressure of 3000 MPa.In this regard, of great interest are the equations which are able to describe PCP with high accuracy when the pressure tends to infinity at a relative density ρ < 1.0.In practice, it occurs when metal powders with low plasticity or very hard powders, such as ceramics, are compacted.To describe such a case, the J. Secondi's equation ( 2002) [52] that contains the relative density parameter at infinite pressure is used.This parameter is denoted as ρ ∞ , and it means that the infinite pressure can occur at a relative density ρ < 1.0.In the original record, the J. Secondi's equation [52] had the form: where ρ ∞ is the relative density at which the compaction pressure tends to infinity, ρ 0 is the initial relative density of the powder, p is the compaction pressure, K and n are constant parameters which control the hardening and plasticity of the powder material.With taking the compacting pressure p as a function and the relative density ρ as an argument, the J. Secondi's Equation ( 17) is converted to the form: where K ′ = 1/K m and m = 1/n are new constants.One should pay attention to the fact that Equation ( 17) turns into Equation ( 14), if the limiting relative density ρ ∞ , will be substituted by the limiting density ρ ∞ = 1.0.Therefore, the Secondi's equation can describe a wider class of powder materials including both hard or brittle and plastic powders.As an example of the approximation by this equation of experimental data on compaction of hard-to-deform plastic and almost nondeformable brittle powders, we used the experimental data on compaction of titanium carbide powder [53] and three iron powders [51] that were compacted to an extremely high density.The values of relative density at different compaction pressures for titanium carbide and three iron powders are given in Table 5.
The results of approximation by the Secondi's Equation ( 18) of these powders are shown in Figure 4, and the values of constant parameters and the coefficient of determination are listed in Table 6.In addition, Figure 4 and Table 6 also show for comparison the results of these powders approximation by Equation (14).As seen, the Equation ( 18) describes PCP for brittle and ductile powders with higher accuracy than Equation ( 14) does.

Material
Table 5. Experimental data on compaction of titanium carbide powders [53] and some iron powders [51].* Iron powders correspond to the following methods of their preparation: FeKH3-vibration grinding, FeKH6-carbonyl, and FeKH9 -steel.

Material
The results of approximation by the Secondi's Equation ( 18) of these powders are shown in Figure 4, and the values of constant parameters and the coefficient of determination are listed in Table 6.In addition, Figure 4 and Table 6 also show for comparison the results of these powders approximation by Equation ( 14).As seen, the Equation ( 18) describes PCP for brittle and ductile powders with higher accuracy than Equation ( 14) does.One of the factors that cause the increased accuracy of the Secondi's Equation ( 18) when describing the compaction process to an extremely high density is probably the presence of three constants.Another important advantage of the Secondi's equation is the fact that it provides for a real powder compaction process, where the compaction pressure tends to infinity at a relative density of significantly below 1.0.The equation of the authors Ge Rong-de, A. B. Yu and Z. H. Gu, L. Parilak and E. Dudrova (14) does not provide for such a possibility.It follows that the equations for PCP with two empirical constants cannot describe this process to the relative density level of 0.96-0.99 with high accuracy.[53] and three iron powders from [51] using Equations ( 18) and (14).It was important to analyze another modified Konopicky-Heckel's equation describing PCP to a high density.We mean the equation of R. Panelli and F. Ambrosio Filho (1998) [54], which in the original record had the form:

Parameter
where A and B are constants, and which can also be written in a shorter form: where ρ ′ 0 = 1 − Exp(−B) is also constant.If pressure is taken as a function, then the Equation ( 19) takes the form: where To evaluate the accuracy of the PCP description by this equation, we have approximated the experimental data on the compaction of two groups of powders.In one group the experimental data show a final relative density less than 0.95, and in the other group-more than 0.95.That is, for the first group we use the data on iron, copper, and nickel powders from Table 1, and for the latter one-iron powders from Table 5.The results of the approximation of PCP for these powders by Equation ( 21) are shown in Figure 5, and the values of the constant parameters and the coefficient of determination R 2 -in Table 7.As seen from the results of this approximation, Equation ( 21) describes PCP quite accurately when the powders are compacted to a density of less than 0.95.If the powders are compacted to a density close to 1.0, then the Equation ( 21) describes the PCP with low accuracy, at which the beginning and middle of PCP are distorted.In this case, it is impossible to correctly assess the features of the powder and its mechanical and deformation properties.It follows that in order to determine the real properties of a powder during its compaction in a die, it is necessary to have, firstly, experimental data on the compaction of a particular powder to a density close to 1.0 and, secondly, an equation that allows to describe such a process with high accuracy.In this regard, the equations for PCP that contain some physical characteristics of the materials used for powder production are of interest.

Equations in Which the Physical Characteristics of Compact Materials Are Used
The above equations for PCP description have basically two constants.In some cases, attempts were made to establish the physical meaning of these constants.e.g., two constants, K and n, in Equations ( 13) and ( 17) characterize, according to the authors, the powder material plasticity and the hardening work during powder deformation.These constants are similar to the coefficients characterizing the plasticity and strain hardening of compact materials.In this regard, some researchers have tried to obtain an equation for As seen from the results of this approximation, Equation ( 21) describes PCP quite accurately when the powders are compacted to a density of less than 0.95.If the powders are compacted to a density close to 1.0, then the Equation ( 21) describes the PCP with low accuracy, at which the beginning and middle of PCP are distorted.In this case, it is impossible to correctly assess the features of the powder and its mechanical and deformation properties.It follows that in order to determine the real properties of a powder during its compaction in a die, it is necessary to have, firstly, experimental data on the compaction of a particular powder to a density close to 1.0 and, secondly, an equation that allows to describe such a process with high accuracy.In this regard, the equations for PCP that contain some physical characteristics of the materials used for powder production are of interest.

Equations in Which the Physical Characteristics of Compact Materials Are Used
The above equations for PCP description have basically two constants.In some cases, attempts were made to establish the physical meaning of these constants.e.g., two constants, K and n, in Equations ( 13) and ( 17) characterize, according to the authors, the powder material plasticity and the hardening work during powder deformation.These constants are similar to the coefficients characterizing the plasticity and strain hardening of compact materials.In this regard, some researchers have tried to obtain an equation for PCP taking the strength and plasticity characteristics of the material used for powder production from reference books.
An attempt to relate one of the equation constants to the yield strength of the powder material was made by S. Torre (1948) [12] and A. N. Nikolaev (1962) [13].They proposed the following equations: where σ S is the yield strength of the powder material and C = 2.5 ÷ 3 is a constant.

123
The Equation ( 23) was modified by G. M. Zhdanovich (1999) [14] as follows: Equation ( 24) does not have the disadvantage inherent in Equation (23), where the pressure becomes negative when the relative density is less than 0.5, but the parameter C is not known in advance and must be determined, as G. M. Zhdanovich points out, from the experiment.The above equations are difficult to use for description of PCP, since the yield strength of particles can differ significantly from the yield strength of an absolutely dense material.
Of interest is an equation that takes into account the features of both ductile and brittle powders, as well as the effect of powder friction against the die wall during PCP.Such an equation was proposed in the work of Li S., Khosrovabadi P.B., Kolster B.H. (1994) [16].Moreover, these authors presented experimental data on the compaction of not only plastic Ni powder but also a mixture of plastic (Ni) and ceramic (Al 2 O 3 ) powders, as well as nickel-coated ceramic (Al 2 O 3 , SiC) powders.They took into account the forces of powder friction against the walls of the die and, as a result, proposed the following equation (in original form): where P is the external pressure, D is the absolute density of the powder blank at pressure P, D 0 is the initial density of the powder, D m is the theoretical density of the powder material, M 0 is the compaction modulus for a dense sample, K is a dimensionless parameter associated with the friction coefficient of the powder (against the die wall) and workpiece geometry, and m and n are empirical constants.In this case, the parameters M 0 and K are also constant.
If the above equation is potentiated and the absolute density is converted into relative density, then the following equation is obtained which associates pressure with the relative density: where the three constants A = M/K, m and n can be determined by fitting the equation to the experimental curve.The authors of this equation showed that the compaction of nickel powder, equivolume mixture of nickel and aluminum oxide powders, as well as aluminum oxide and silicon carbide powder coated with aluminum can be described by Equation (25) with high accuracy [16].However, it should be kept in mind that the experimental data obtained by authors apply to the materials with a final relative density less than 0.8.For approximation by Equation (26) the experimental data for higher final density, e.g., the data on the compaction of iron and copper powders from Table 1, as well as unusual data on the compaction of coarse and fine iron and copper powders from the book by F. V. Lenel [55], (p.96, Figures 3-25) were of interest.The experimental data on compaction of coarse and fine Fe and Cu powders, presented in [55], were obtained by digitizing the experimental points on the corresponding curves, and these data are given in Table 8.Note that in [55] only the bulk or apparent density of coarse and fine iron and copper powders is indicated without information on the sizes of the powders.
The results of approximation by Equation (26) of experimental data on the compaction of various iron and copper powders are shown in Figure 6, and the values of the constant parameters of Equation ( 26) and the coefficient of determination R 2 are listed in Table 9. on the corresponding curves, and these data are given in Table 8.Note that in [55] only the bulk or apparent density of coarse and fine iron and copper powders is indicated without information on the sizes of the powders.The results of approximation by Equation ( 26) of experimental data on the compaction of various iron and copper powders are shown in Figure 6, and the values of the constant parameters of Equation ( 26) and the coefficient of determination R 2 are listed in Table 9.Table 9. Constant parameters of Equation ( 26) obtained via approximation of experimental data on the compaction of iron and copper powders from Table 1 and iron and copper powders from Table 8.The presented results of approximation of experimental data on the compaction of various iron and copper powders by Equation (26) show that the equation proposed by previously mentioned authors [16] allows description of PCP quite accurately in cases where the powder is compressed to a low density.And this happens despite the fact that the proposed equation has three constants that implicitly take into account the resistance of the powder material to deformation, the powder friction against the die walls, and the compacted powders geometry.In addition, a characteristic feature of Equation ( 26) is the fact that it provides for obtaining a final pressure when the relative density of the powder reaches 1.0.Moreover, the value of this pressure (see p max in Table 9) is unexpectedly low that confirms the impossibility of this equation to accurately describe PCP at the final stage.

Parameter
Noteworthy is another equation for PCP presented in one of the works by G. Aryanpour and M. Farzaneh (2015) [17].When deriving the equation, the authors analyzed various mechanisms of PCP, in particular-particles rearrangement and plastic deformation.With that, they made the assumption that these two mechanisms work together up to a sufficiently high relative density of 0.95, at which, as they believe, there are no more open pores, and the particle rearrangement mechanism does not operate.To describe compaction due to the particle rearrangement mechanism, the authors suggested the following equation: where b and a are constants determined from the experiment, ρ 0 is the initial density of the powder in the die.After a simple transformation and taking logarithm, the Equation ( 27) takes a form: where c = 1/a is a constant.
To describe the powder compaction due to plastic deformation, the researchers [17] used the Heckel's Equation [11] written in the form: where k is a constant, but ρ ′ 0 is the conditional initial density of the powder.According to G. Aryanpour and M. Farzaneh, these two mechanisms are summed up during compaction of the powder, and as a result, the following equation was proposed (in this case, the conditional initial density ρ ′ 0 in Equation ( 29) is replaced by the actual initial density ρ 0 ): The Equation (30), as noted by its authors, can be used to describe PCP in the density range from ρ 0 to ρ = 0.95.It takes into account the action of two mechanisms-particles rearrangement and particles deformation.In this case, at the beginning of the compaction process, the rearrangement mechanism prevails, and at the end of the compaction, the deformation mechanism does.Unfortunately, this equation cannot be transformed in such a way that it could be solved with regard to the pressure p and an appropriate approximation can be made.However, to estimate the accuracy of approximation of various experimental data by Equation (30), it is necessary to transform it so that the function is the relative density ρ.In this case, it will acquire the following form: But the Equation (31) hides compaction mechanisms, making it less useful.At the same time, this equation makes it possible to approximate experimental data on powder compaction and determine the values of three constants-b, c, and k.Taking into consideration that this equation has three constants, its adequacy in describing PCP not only up to a relative density of 0.95 but also up to a higher density is of interest.Therefore, to approximate the experimental data on PCP by this equation, we chose the experimental Powders 2024, 3 126 results previously used by us on the compaction of iron, copper, and nickel powders, presented in Table 1, as well as on iron powders compacted to high density presented in Table 5.The results of the approximation of experimental data on the compaction of these powders are shown in Figure 7.The values of the constants of Equation ( 31) and the coefficient of determination R 2 are given in Table 10.
But the Equation ( 31) hides compaction mechanisms, making it less useful.At the same time, this equation makes it possible to approximate experimental data on powder compaction and determine the values of three constants-b, c, and k.Taking into consideration that this equation has three constants, its adequacy in describing PCP not only up to a relative density of 0.95 but also up to a higher density is of interest.Therefore, to approximate the experimental data on PCP by this equation, we chose the experimental results previously used by us on the compaction of iron, copper, and nickel powders, presented in Table 1, as well as on iron powders compacted to high density presented in Table 5.The results of the approximation of experimental data on the compaction of these powders are shown in Figure 7.The values of the constants of Equation ( 31) and the coefficient of determination R 2 are given in Table 10.  1, as well as iron powders (d) FeKH3, (e) FeKH6, and (f) FeKH9 from Table 5.  1, as well as iron powders (d) FeKH3, (e) FeKH6, and (f) FeKH9 from Table 5.
Table 10.Constant parameters in Equation ( 31) resulted from the approximation of experimental data on the compaction of iron, copper, and nickel powders from Table 1 and iron powders FeKH3, FeKH6, and FeKH9 from Table 5.As seen from presented results, the accuracy of the approximation of experimental data on the compaction of various powders by Equation ( 31) is high for cases of powder compaction up to a density of 0.95.This may indicate the adequacy of the hypothesis that allows the simultaneous existence of the mechanisms of rearrangement and deformation of particles during the compaction process up to a relative density of 0.95.Alongside with this, the Equation (31) can, in some cases, describe with high accuracy the PCP up to a relative density exceeding 0.95 (Table 10, FeKH9 powder).It is also important that this equation allows one to estimate the degree of density change during powder compaction due to both particles rearrangement (∆ρ r ) and plastic deformation (∆ρ d ).

Parameter Powder
In order to evaluate the change in density by two mechanisms, it was suggested to determine the derivatives of the relative density with regard to pressure separately for Powders 2024, 3 127 the rearrangement and plastic deformation mechanisms [17].If we fully expand these derivatives, then they will look as: -for particles rearrangement (32) and -for plastic deformation (33) The areas under the curves ( 32) and ( 33) within the pressure range from zero to a value corresponding to a given density will show the degree of change in density for the corresponding compaction mechanism.Numerical integration of dependences (32) and (33) for the examples shown in Figure 7, within the pressure range from zero to the pressure corresponding to a density of 0.95, led to the results presented in Table 11.In addition, this Table 11 presents the true value of particles plastic deformation, which, for the case of reaching the relative density of 0.95, was calculated using the following formula: Table 11.Degree of change in density due to rearrangement and plastic deformation of particles during compaction of the powder from the initial density to a density of 0.95, as well as the true degree of plastic deformation of the particles and the degree of change in the relative density from initial to 0.95.

Density Change Mechanism Powders
Fe [ By the data in Table 11, the degree of change in the powder density due to rearrangement during compaction in a rigid die can exceed the degree of change in density through plastic deformation of the particles.In relatively hard powders, such as FeKH3, FeKH6, and FeKH9, densification occurs practically due to the particles rearrangement.There cannot be excluded the case when description of solid powders particles rearrangement with the Equation ( 27) is very approximate.In this regard, I would also like to draw attention to the fact that the Equation ( 27) exactly corresponds to the equation of M. Yu.Bal'shin (1972) [18] (p.163, Eq.V.47) which he proposed to describe the third stage of the densification process, i.e., the stage where practically there is no rearrangement of the particles.
It should also be noted that the process of particle rearrangement during compaction of powders in a rigid die is explicitly reflected in the equation proposed as far back as in the middle of the last century (1962) by A. R. Cooper and L. E. Eaton [26].In the original record this equation is as follows: where V* is the degree of powder sample volume change under pressure, V 0 , V and V ∞ are initial volume at zero pressure, current volume at pressure P, and volume at infinite pressure, respectively, P is the compaction pressure, a 1 and k 1 are constants characterizing the rearrangement process, and a 2 and k 2 are constants characterizing the process of plastic deformation.
Powders 2024, 3 128 When compacting the powder, it is advisable to operate with the relative density ρ = V com /V (V com is the volume of a pore-free sample), then the Equation ( 35) can be converted into the dependence of relative density (ρ) on pressure (p): where ρ 0 is the initial density of powder, ρ ∞ ≤ 1 is the density at infinite pressure.This equation was proposed in order to describe the compaction process of ceramic powders.Therefore, it can be used to approximate the compaction process of any hard powders.In this regard, the applicability of Equation ( 36) was tested for describing the process of compaction of titanium carbide powders and iron powders presented in Table 5.The results are shown in Figure 8, and the values of the constant parameters of Equation ( 36) and the coefficient of determination R 2 are listed in Table 12.The approximation was performed at a given value of the parameter ρ ∞ for each powder.
where V* is the degree of powder sample volume change under pressure, V0, V and V∞ are initial volume at zero pressure, current volume at pressure P, and volume at infinite pressure, respectively, P is the compaction pressure, a1 and k1 are constants characterizing the rearrangement process, and a2 and k2 are constants characterizing the process of plastic deformation.
When compacting the powder, it is advisable to operate with the relative density ρ = Vcom/V (Vcom is the volume of a pore-free sample), then the Equation ( 35) can be converted into the dependence of relative density (ρ) on pressure (p): where ρ0 is the initial density of powder, ρ∞ ≤ 1 is the density at infinite pressure.This equation was proposed in order to describe the compaction process of ceramic powders.Therefore, it can be used to approximate the compaction process of any hard powders.In this regard, the applicability of Equation ( 36) was tested for describing the process of compaction of titanium carbide powders and iron powders presented in Table 5.The results are shown in Figure 8, and the values of the constant parameters of Equation ( 36) and the coefficient of determination R 2 are listed in Table 12.The approximation was performed at a given value of the parameter ρ∞ for each powder.5.
Table 12.Constant parameters in Equation ( 36) resulted from the approximation of experimental data on the compaction of titanium carbide and iron powders FeKH3, FeKH6, and FeKH9 from Table 5.   5.It follows from the presented approximation results that Equation (36) that contains four constant parameters (the fifth parameter ρ ∞ is predetermined) is capable of describing the compaction process of hard metal and ceramic powders with high accuracy.The most important result of approximating the process of compaction of hard powders by this equation is the fact that the process of particle rearrangement affects the change in porosity to a much greater extent than the process of powder deformation.This is evidenced by the value of the constant a 1 in Table 12, which reaches 80% of the entire compaction process.A similar result follows when describing the process of compaction of hard iron powders by Equation (31), Table 11.2017) [56] is worth of attention.As said above, these authors previously proposed an equation with two constants, which in many cases allows one to describe PCP quite accurately.In new article, they proposed a novel equation based on their former one (13) which was transformed into an equation with one constant due to the relation established between the constants K and n in the form ln(K) = 1.2952 − 7.3349n.With this, Equation ( 13) after taking a double logarithm acquires the form:

Modern Works
where P 0 and P are the initial and current porosities of the powder sample.
After potentiation, the Equation ( 37) is transformed into the following expression: from which one can get an equation where the compaction pressure is a function, and porosity or relative density (ρ = 1 − P) is an argument: where m = 1/n is a new constant.
There is an ambiguous attitude to this new equation for PCP.On one hand, with an equation containing one constant parameter it is easier to describe the compaction of different powders, but, on the other hand, such an equation is not capable of taking into account the great diversity in the morphology and properties of powders.All differences in the behavior of powders during compaction are averaged in one parameter.This may reduce the accuracy of the densification process description of certain powders.Indeed, the approximation by Equation (39) of experimental data on the compaction of iron, copper and nickel powders (Table 1) showed a lower accuracy in two cases out of three (Table 13) as compared with the approximation by two-parameter equation ( 14) (Table 4).With a reduced accuracy in describing PCP, we will get to know with reduced accuracy the degree of plastic deformation of particles, the work of plastic deformation, the magnitude of the resistance of particles to deformation, as well as other characteristics of the powder.Therefore, a one-parameter equation can be used to approximate assessment of the compaction process and powder properties.
Another recent work on compressibility of powders during compaction in a die deserves an attention, namely, the article of J.M. Montes, F.G. Cuevas, J. Cintas et al. (2018) [57].When developing their equation, the authors aimed at minimizing the number of experimental constants in the equation and using physical parameters, which are characteristic of the powder material and known in advance from reference books.Unfortunately, in many cases, such a desire cannot lead to obtaining an adequate compaction equation due to significant differences in the strength and plastic properties of powders and bulk material they are made of.The researches [57] proposed the following equation to describe PCP: where P N is the external pressure during compaction of the powder, Θ and Θ M are the current and initial porosity of the powder sample, ξ is the coefficient that takes into account the friction of the powder against the die wall, k and n are the Hollomon equation parameters which characterize the law of hardening during deformation of a dense sample obtained from compacted powder.
According to the authors of [57], the parameters k and n are very close to the reference values for the powder material.Therefore, in Equation ( 40), the parameters k and n are known in advance, and the parameters Θ M and ξ are determined by fitting the curve according to Equation (40) to the experimental compaction curve of a particular powder.In the authors' opinion, the advantage of Equation ( 40) is the use of previously known parameters k and n with a clear physical meaning.However, the proposed equation does not take into account the phenomenon of particle rearrangement during powder compaction, as well as the possible discrepancy between the strain hardening of particles themselves and the strain hardening of a compact sample made from these particles.
In reality, such factors as the rearrangement of particles and the specific nature of the particle strain hardening, especially in the compaction of fine-grained particles, are of fundamental importance in describing the compaction process and for the adequacy of the corresponding equation.Therefore, the free status of the parameters k and n in Equation ( 40) can take into account both the particle rearrangement factor and the specific nature of the strain hardening of these particles.In this regard, an attempt was made to use Equation (40) to approximate the experimental curves for the compaction of various powders, provided that the parameters k and n are free, and the initial porosity is known.To perform the approximation, we transformed Equation ( 40) into the dependence of the compaction pressure p on the relative density ρ: where K ≈ k (since ( √ 3/2) × ξ ≈ 1) and n are constant parameters obtained from the experiment, ρ 0 is the initial relative density of the powder in the die.
To carry out the approximation, the experimental data for iron, copper and nickel powders from Table 1 as well as hard-to-compact iron powders from Table 5 were used.The results of the approximation are shown in Figure 9, and the values of the constants are listed in Table 14.
where K ≈ k (since (√3/2) × ξ ≈ 1) and n are constant parameters obtained from the experiment, ρ0 is the initial relative density of the powder in the die.
To carry out the approximation, the experimental data for iron, copper and nickel powders from Table 1 as well as hard-to-compact iron powders from Table 5 were used.The results of the approximation are shown in Figure 9, and the values of the constants are listed in Table 14.The results of the approximation of various powders by Equation (41) show that in some cases this equation allows one to describe very accurately the PCP, despite the change in some principles underlying the derivation of Equation ( 40).However, hard   41) obtained via approximation of experimental data on the compaction of iron [39], copper [40], and nickel [41] powders, as well as FeKH3, FeKH6, and FeKH9 iron powders from the book [51].The results of the approximation of various powders by Equation (41) show that in some cases this equation allows one to describe very accurately the PCP, despite the change in some principles underlying the derivation of Equation (40).However, hard powders compacted to a density greater than 0.95 cannot be adequately described by this equation.

Parameter
A. Molinari et al. (2018) [27] also proposed a new equation for powder compaction in a rigid die.The authors studied the compaction of low-alloyed iron powders (alloyed with CuMo, Mo, CrMo, or CuMoNi) in a mixture with graphite and lubricant, and suggested describing PCP by the absolute density dependence on the average compaction pressure in the form of a power law: where ρ 0 and ρ are the absolute initial and current density of the powder in the die (the dash above the symbol distinguishes the absolute density from the relative density, which we denote as ρ), respectively, P m is the average compaction pressure, A and B are constant parameters obtained from the experiment.It should be noted that a feature of the experimental PCP in article [27] is a relatively narrow range of compaction pressure up to 540 MPa, as well as the use of the absolute Powders 2024, 3 132 density of the material instead of the relative density.If we transform Equation ( 42) into the dependence of pressure on relative density, then we obtain power law: where ρ 0 and ρ are relative initial and current densities, respectively, A ′ = 1/A (1/B) and B ′ = 1/B are constant parameters obtained from the experiment.Note that using of the power law to describe PCP is not a new solution (see Equation (3)), but Equation ( 42) takes into account the real fact that at a density equal to the initial density ρ 0 , the compaction pressure is zero.However, in the long-known equation ( 3), it is easy to take into account this fact if instead of the current relative density ρ we substitute the expression (ρ − ρ 0 )/(1 − ρ 0 ), that allows taking into account not only the initial density but also the entire range of density changes from ρ 0 to ρ = 1.Then Equation (3) will acquire the form: where c and d are constant parameters obtained from the experiment.
Comparison of the adequacy of Equations ( 43) and (44) showed that in describing the compaction of iron, copper, and nickel powders that we used the Equation ( 44) was more accurate.This follows from the approximation results presented in Figure 10 and in Table 15.As seen from Figure 10 and Table 15, the power-law dependence of pressure on relative density does not provide a high accuracy of the PCP description, especially in the initial stage of compaction, which is described by Equation ( 43) less accurately.Obviously by this reason the authors [27] divided PCP into two stages, in one of which particles are rearranging, and in the other, plastic deformation occurs.In this case, each stage of compaction was approximated by a power equation with its own constant parameters.However, such a method for PCP description was proposed by researches [44] as early as in 1975.
The ability to describe almost the entire compaction process with a single equation raises the question of how one equation can describe the different stages of PCP.In this   43) and ( 44) after approximation of experimental data on the compaction of iron, copper, and nickel powders from Table 1.As seen from Figure 10 and Table 15, the power-law dependence of pressure on relative density does not provide a high accuracy of the PCP description, especially in the initial stage of compaction, which is described by Equation ( 43) less accurately.Obviously by this reason the authors [27] divided PCP into two stages, in one of which particles are rearranging, and in the other, plastic deformation occurs.In this case, each stage of compaction was approximated by a power equation with its own constant parameters.However, such a method for PCP description was proposed by researches [44] as early as in 1975.

Parameter
The ability to describe almost the entire compaction process with a single equation raises the question of how one equation can describe the different stages of PCP.In this case, two possibilities may exist: either the compaction process occurs according to one yet unknown mechanism, or the equation must contain a number of constants, which take into account, to varying degrees, the manifestation of different mechanisms.In the next two messages of this article (Part 2 and Part 3), new equations for powders compaction in a rigid die are proposed, which contain three and four constants.But in one case (Part 2), to obtain such equations, various equations of M. Yu.Bal'shin were corrected, and in another case (Part 3), the plasticity equations of a porous body proposed by Kuhn and Downey or Skorokhod, Martynova and Shtern were transformed.
According to the results of analysis concerning the equations presented in the form of simple mathematical functions or in the form of semi-empirical equations, which take into account some physical constants of the powder material, the following conclusion can be formulated.

1.
The considered equations for the compaction of powders in a rigid die in most cases contain two constants, which depend on the powder material and the type of equation and are determined from the experiment.At the same time, the accuracy of describing the process of powder compaction by equations with two constants is different.In addition, some equations lead to a finite value of compaction pressure upon reaching the pore-free state of the sample, while the others lead to an infinite pressure.

2.
High accuracy in describing the process of powder compaction using different equations is observed in the case of using a shortened range of density changes, i.e., up to 0.8-0.9, but such a description of the compaction process is not complete, since fundamental differences in the behavior of powders are observed when a relative density close to unity is reached.When compacting powders to a density greater than 0.95, the highest accuracy in describing the compaction process can be provided by equations with three constants, for example, such as the J. Secondi's equation.

3.
Increasing the number of constants in the equation to three or four allows us to identify compaction mechanisms such as particle rearrangement and plastic deformation.This is shown by discussing the equations of A. R. Cooper and L. E. Eaton or G. Aryanpour and M. Farzaneh, which consider that the rearrangement and plastic deformation mechanisms operate in parallel throughout the compaction of powders to a relative density of 0.95.4.
In a number of proposed equations, there was an attempt to take into account both the strength and plasticity characteristics of a powder compact material and the shape and size of particles, as well as the internal friction forces between particles and external friction forces between particles and the die wall.But this came down to setting the appropriate coefficients, which are constant throughout the entire compaction process.Moreover, in most cases, the value of these coefficients must be determined in the experiment and this neutralizes the desire of researchers to work without using empirical constants.

Figure 1 .
Figure 1.Approximation dependences of experimental data on PCP for (a) iron, (b) copper, and (c) nickel using exponential Equation (2) (dashed line) and power Equation (3) (solid line).Here and below, light circles are the experimental data.

Figure 1 .
Figure 1.Approximation dependences of experimental data on PCP for (a) iron, (b) copper, and (c) nickel using exponential Equation (2) (dashed line) and power Equation (3) (solid line).Here and below, light circles are the experimental data.

Figure 2
Figure 2 and Table3show that the new Equation (8) describes PCP with higher accuracy than the original Heckel's equation(6).In this case, the value of the conditional initial density ρ c in Equation (8), obtained as a result of the approximation, turns out to be lower than the actual initial density ρ 0 that makes it possible to accurately describe the initial stage of the compaction process.A higher accuracy of approximation of experimental data on PCP compared to the equation of K. Konopicky or R. W. Heckel can be reached by another modernization of the equation, which leads to the appearance of the so-called double logarithmic form of the equation for powder compaction in a rigid die.Three different teams of authors proposed such an equation independently.First, the work of Ge Rong-de (1991)[47] should be noted, in which he proposed a new differential equation for PCP in a rigid die:
Describing the Process of Powder Compaction in a Rigid Die Recently published papers with new equations for PCP are of increased interest.First of all, a new article by L. Parilak, E. Dudrova et al. (

Figure 9 .
Figure 9. Approximation dependences according to Equation (41) as a result of processing experimental data on the compaction of (a) iron, (b) copper, and (c) nickel, as well as hard-to-compress iron powders (d) FeKH3, (e) FeKH6, and (f) FeKH9.

Figure 9 .
Figure 9. Approximation dependences according to Equation (41) as a result of processing experimental data on the compaction of (a) iron, (b) copper, and (c) nickel, as well as hard-to-compress iron powders (d) FeKH3, (e) FeKH6, and (f) FeKH9.

Figure 10 .
Figure 10.Approximation curves according to Equation (43) (dotted line) and Equation (44) (solid line) as a result of processing experimental data on the compaction of (a) iron, (b) copper, and (c) nickel powders.

Figure 10 .
Figure 10.Approximation curves according to Equation (43) (dotted line) and Equation (44) (solid line) as a result of processing experimental data on the compaction of (a) iron, (b) copper, and (c) nickel powders.

Table 2 .
Constants and coefficient of determination R 2 obtained via approximation of experimental data on PCP for iron, copper, and nickel using exponential (2) and power (3) equations.

Table 6 .
Constants and coefficient of determination R 2 obtained via approximation of experimental data on the compaction of titanium carbide powder

Table 7 .
(21)tants and coefficient of determination R 2 obtained via approximation of experimental data on the compaction of powders of iron, copper, nickel, and FeKH3, FeKH6 and FeKH9 using Equation(21).

Table 7 .
Constants and coefficient of determination R 2 obtained via approximation of experimental data on the compaction of powders of iron, copper, nickel, and FeKH3, FeKH6 and FeKH9 using Equation (21).

Table 8 .
Experimental data on the compaction of coarse and fine iron and copper powders[55].

Table 8 .
Experimental data on the compaction of coarse and fine iron and copper powders[55].
It follows from the presented approximation results that Equation (36) that contains four constant parameters (the fifth parameter ρ∞ is predetermined) is capable of describing the compaction process of hard metal and ceramic powders with high accuracy.The most

Table 12 .
Constant parameters in Equation (36) resulted from the approximation of experimental data on the compaction of titanium carbide and iron powders FeKH3, FeKH6, and FeKH9 from Table5.

Table 13 .
(39)value of the constant m in Equation(39)and the coefficient of determination R 2 obtained after the approximation of experimental data on compaction of iron, copper, and nickel powders from Table1.

Table 14 .
Constant parameters in Equation (

Table 15 .
(44)tants in Equations (43) and(44)after approximation of experimental data on the compaction of iron, copper, and nickel powders from Table1.