Derivation of the Coefﬁcients in the Coulomb Constant Shear Friction Law from Experimental Data on the Extrusion of a Material into V-Shaped Channels with Different Convergence Angles: New Method and Algorithm

: The combined Coulomb constant shear friction law is widely used in commercial and research software for the ﬁnite-element analysis (FEA) of metalworking and is naturally more ﬂexible and hence, more relevant to real-life manufacturing than the individual Coulomb and constant shear friction laws. In this work, a new mathematical model of coefﬁcients in the Coulomb constant shear friction law for extruding a metal through narrow V-shaped channels with small convergence angles has been developed and evaluated and compared with laboratory measurements. The extrusion of the model material (lead) through narrow V-shaped channels with small convergence angles varying from 0 to 3.5 degrees has been studied. The Coulomb friction coefﬁcient µ and the constant friction factor m appear to be independent of the dimension ratio and are inﬂuenced mostly by roughness and range from µ = 0.363 (with lubricant) to µ = 0.488 (without lubricant) and from m = 0.726 (with lubricant) to 0.99 (without lubricant). The relative length dominated by the Coulomb friction law is less than 1%, and the Coulomb’s coefﬁcient of friction can be approximated as 12 the constant shear friction factor for all tested cases. The developed method and algorithm can be used in both FEA of manufacturing processes and efﬁciency tests for lubricants used in metalworking.


Introduction
The critical importance of friction for metal-forming processes and its adverse impacts on wear resistance of dies are well-established. Mathematical modeling of thermomechanical and chemical interactions at the workpiece-tool interface and other relevant phenomena directly affecting such key process parameters as productivity, quality, and scrap level relies heavily on the finite element method (FEM) [1][2][3][4] widely used for the analysis (FEA) and optimization of various processes and technologies. FEM, a vital component of commercial and research software, is generally capable of delivering great results, but we need a precise description of the behavior of visco-plastic materials at elevated temperatures and accurate assignment of boundary conditions, especially those related to friction [5][6][7][8][9][10][11].
However, despite impressive progress achieved in both theory and over past decades, there is still significant uncertainty in friction models and their parameters in the case when friction occurs at high temperatures [12]. In most cases, non-linear FEM simulations employ phenomenological friction models [13,14], allowing us to uniquely set the boundary conditions for stresses at the workpiece-tool interface via the friction coefficients that greatly improve the convergence of numerical solutions for strongly deformed parts of complex geometries at reasonable computational costs. In view of these circumstances, the correct choice of the friction model and accurately determining the corresponding friction coefficients are a must in the case when a real-life manufacturing process is being simulated.
In the modeling of bulk metal forming, there are a few options to choose from, such as the Coulomb friction law [15][16][17], the constant shear (Siebel) friction law [18][19][20][21][22][23][24], both Coulomb and constant shear friction laws [25,26], general ones [27][28][29] or the integrated law, i.e., Coulomb constant shear friction law [9,14,[30][31][32][33]. Being naturally more flexible and, hence more relevant to real-life manufacturing than conventional Coulomb and constant shear friction laws, the Coulomb constant shear friction law is widely used in FEA of metalworking. Here, it is important to note that a fairly common approach to deduce the Coulomb constant shear friction law coefficients by equating calculated and measured total power-energy properties can lead to large uncertainty in the flow stress and strain rates that in turn may result in fatal errors in the design of dies.
In this work, a new mathematical model of the Coulomb constant shear friction coefficients for extruding a metal through narrow V-shaped channels with small convergence angles was derived based on an approximate solution of the equilibrium equations of elasticity, the kinematically admissible velocity field, and the distribution of tangential stresses in the plastic deformation zone meeting realistic physics-based boundary conditions. Based on the above-mentioned model, the algorithm of the derivation of the Coulomb constant shear friction law coefficients from laboratory measurements was developed.
In order to determine the friction coefficients, classical experiments on the upsetting of a ring-shaped workpiece between flat plates have been performed [25,26,34,35]. In [36], a method of the derivation of the friction coefficients from measurements on the upsetting of a cylindrical workpiece at high strain rates was reported. Different tests to determine friction coefficients under different forming conditions have been developed on the basis of extrusion processes [8,9,21,23,24,32,33], especially when other methods, such as the ring compression test, are not suitable for metal-forming processes where the surface expansion is high [18]. Some examples are the double-cup extrusion test [19], the boss and rib test [20], backward extrusion [37], or combined forward-backward extrusion [23], among others.
However, it is important to note that in experiments when both the metal flow and the geometry of the deformation zone are simple, the deformation pressure is low, and the free surface formed due to the deformation is fairly small, and, hence, the friction coefficients/factors are unreliable and cannot be used to simulate real-life manufacturing processes [38]. According to [8,38], measurements of forward and backward extrusion provide much more realistic estimates.
In addition to the Coulomb friction law, the constant shear (Siebel) friction law, general ones or the integrated Coulomb constant shear friction law used in Deform [13] and Abaqus [14] software, other laws have been developed and used in bulking-forming modelling, such as the Levanov friction law used in QForm software [39]. The Levanow friction model represents a generalization of the shear friction model and Coulomb friction model. The similar Wanheim and Bay [2,[27][28][29] friction law assumes friction to be proportional to the normal contact pressure at low normal pressure, but progressively approaching a constant value at high normal pressure.
In [38,40], experiments on extrusion into narrow lateral radial channels were used to study friction in the case when the contact area between the tool and the workpiece is large.
In a previous study [38], a comparative analysis of the lubricant's composition influence on the frictional forces was proposed by analyzing the filling of the V-shaped channels of the die. The analysis was carried out by the similarity theory method. In the present study, the Coulomb constant shear friction parameters were derived from laboratory measurements on the extrusion of a model material (lead) into V-shaped channels with different convergence angles using the lower bound limit analysis approach. The flowing into narrow channels formed by die surfaces is critically important for pressure processing of metals, with flowing into the groove of the die streams during hot stamping and into the matrix during drawing being the most typical examples of such processes.
The large contact area and low local deformations typical of these processes have a large impact of friction on the flow pattern. Since friction in the large contact area may differ greatly depending on the conditions of the die surface and changes in the contact pressure, the Coulomb constant shear friction law is better suited for the description of friction than either the Coulomb friction law or the constant shear friction law, and thus, the friction coefficients derived from analysis of the occupancy of narrow V-shaped channels in the present study are more relevant to the real-life manufacturing than those obtained using other methods.
The core idea of the developed method and the main component of its novelty involve the use of a comparative analysis of work of filling channels with different angles of convergence calculated from the experimental data on the volumes of metal extruded into these channels. Depending on the friction coefficients values, different ratios of the lengths of material extruded into the channels are observed. A particular advantage of this method over competitors is its independence of the yield strength, leading to its invariance with respect to the material being processed and independence of the processing conditions due to the simultaneous extrusion into all radial channels during experiments. The developed mathematical model was used to study the influence of the dimension ratio on the extrusion into the channels and Coulomb constant shear friction coefficients. Experiments on bulkforming process through matrices of special design were carried out with lead chosen for its high visco-elasticity, allowing us to simulate the metal flowing in the stamp during the hot stamping under laboratory conditions. The developed model was evaluated by comparing modeled and experimental friction coefficients by analyzing extrusion into two dies of identical shape but with different, by a dimensional ratio of two, linear dimensions of the forming cavity.

Theoretical Simulation of Material Flow into V-Shape Narrow Channels
In this section, easy-to-use analytical equations describing the Coulomb constant shear coefficients as a function of the channel properties and other relevant variables describing extrusion conditions and properties of a material being extruded are derived. The channel height is assumed to be much larger than wider; therefore, the metal flow in the channel can be considered two-dimensional.
The most common friction models used to describe metalworking processes are [11]: the Coulomb friction law: the constant shear friction model: the generalized friction model [27][28][29]: where τ is the friction stress, µ is the friction coefficient, p is the normal stress, k = σ s √ 3 -is the shear stress, f is the friction coefficient under real contact conditions, and α is the ratio of actual and theoretical contact areas. The tangential stress along the V-shaped-generating lines of the channel is determined by the constant shear friction law, and the Coulomb friction law τ = µσ n at σ n < σ p Here, σ s is the yield strength, σ n is the normal stress on the contact surface, and σ p is the maximum of σ n when the Coulomb friction law changes to the constant shear friction law. First, a cylindrical coordinate system O(r,θ,z) is introduced in which the plane θ = 0 is the mid-plane of the generating lines of the converging channel (see Figure 1), where θ = α/2, α is the convergence angle of the generating lines, r 0 is the radius of the generatrix of the V-shaped channel with entrance thickness s (the radius of the element extruded into the channel with a base length l), r 1 is the apex radius of the element extruded into the channel, and r 1 ≤ r ≤ r 0 is the current radius.
The tangential stress along the V-shaped-generating lines o mined by the constant shear friction law,

= √3
at ≥ and the Coulomb friction law

= at <
Here, σs is the yield strength, σn is the normal stress on the contac maximum of σn when the Coulomb friction law changes to the con First, a cylindrical coordinate system O(r,θ,z) is introduced in is the mid-plane of the generating lines of the converging channel = α/2, α is the convergence angle of the generating lines, r0 is the r of the V-shaped channel with entrance thickness s (the radius of the the channel with a base length l), r1 is the apex radius of the elem channel, and r1 ≤ r ≤ r0 is the current radius.  According to (1), there is a radius r = rp, at which σn = σp(rp) a condition is met: Figure 1. Scheme of the determination of the Coulomb constant shear friction law coefficients. α-convergence angle, r 0 -radius of the generatrix of the V-shaped channel, r 1 -apex radius of the element extruded into the channel, r p -transition area radius from constant to Coulomb friction law. λ = r 0 r 1 is the elongation ratio of the extruded element, and λ l = l s is the relative length of the extruded element to make a transition to dimensionless quantities.
For small angles α (in rad.), the following relationships are used: According to (1), there is a radius r = r p , at which σ n = σ p (r p ) and, thus, the following condition is met: For narrow channels, i.e., for small α values, the linear distribution law for tangential stress τ rθ along the angle θ is used. By making use of (4), (5), (9) and the boundary conditions τ rθ = 0 at θ = 0; τ rθ = τ, at θ = α/2, we obtain: The equilibrium equation in cylindrical coordinates in the projection onto the axis θ has the following form: 1 r ∂σ θθ ∂θ where σ θθ may be expressed as a product of two functions, σ(θ) and (σ(r) + C 0 ).
After inserting the expressions for the distribution law of the tangential stress (4) and σ θθ (6) into the equilibrium Equation (11), Equation (13) is obtained: After the integration of (13) over θ, Equation (14) is obtained: In the case where θ ≤ α/2 is small and θ 2 is the second-order infinitesimal O(10 −2 ), σ(θ) ≈ C θ is obtained if the following condition, proved below, is met. From (12), with accounting for (14) and (15), Equation (16) is obtained: The typical velocity field of a metal flow in a converging channel is described by the following dependency: and, therefore, the effective strain rate is expressed as: which at small angles θ = 0 . . . α/2 reduces to: After making use of the constitutive equations: where σ is the hydrostatic pressure, and the continuity condition for the metal flow is: Finally, Equation (22) is obtained: The stress law is found: where, according to (10): The constant C r in (24) is derived using the following boundary condition on the free surface σ rr (r = r 1 ) = 0 (27) from which the following are obtained: From (28), it follows that in (12): After inserting expressions for σ(r) and C 0 from (28) into (15), the following inequity is shown: if r 0 ≥ r > r 1 holds identically, assumption (15) is proved to be true. The upper limit of C is shown in Figure 2.
As can be seen in Figure 2, the value of C asymptotically tends to zero when the elongation coefficient increases.
Contact pressure σ n decreases with decreasing radius r. After the radius reaches r = r p corresponding to λ p = r p r 1 , a transition from the constant shear friction law to the Coulomb friction law occurs: and therefore:  Contact pressure σn decreases with decreasing radius r. After the radius rea rp corresponding to = 1 , a transition from the constant shear friction law to t lomb friction law occurs: (1 + (1 + ) ) 2 = √3 2 and therefore: The channel numbers associated with the index are denoted as i = 1…N.
Since μi is independent of angle αi, the following condition The channel numbers associated with the index are denoted as i = 1 . . . N.
Since µ i is independent of angle α i , the following condition should be met at any i, j = 1 . . . N. Figure 3 presents the Coulomb friction coefficient vs. the constant shear friction coefficient at λ l1 = 10.
Elongation coefficients λ pi for each channel can be found from: which follows from (37) with accounting for (36). It is clear that for channel 1 α 1 = 0, λ p1 = 1.  Elongation coefficients λpi for each channel can be found from: which follows from (37) with accounting for (36). It is clear that for channel 1 α1 = 0, λp1 = 1. Figure 4 presents the dependence of the elongation coefficients     Elongation coefficients λpi for each channel can be found from: which follows from (37) with accounting for (36). It is clear that for channel 1 α1 = 0, λp1 = 1. Figure 4 presents the dependence of the elongation coefficients = 1 on the convergence angle of channel i and the constant shear friction coefficient at λl1 = 10.  As can be seen in Figure 4, the Coulomb friction zone increases with increasing convergence angle of each channel.
Since the plastic deformations zones at the entrance to each channel i with identical dimensions s and h of the entry hole size are equal, the work of A 0i of forcing the metal through the plastic deformation zone at the entrance to channel i is proportional to volume W i of the material extruded into channel i.
where σ m is equal for all channels. The work regarding the plastic deformation associated with the extrusion into channel i is expressed by Equation (40).
where the first term on the right hand side of (40) is the work of deformation before the entrance to channel i and the second is the work of metal forming in channel i. The case where the length of the extruded section exceeds l = r 0 − r p is considered. Since at r = r p , the constant shear friction law transforms into the Coulomb friction law, the work of the formation of the section extruded into the radial channel has, according to (40), the following equation.
Integrating (42) and taking into account the last expression, Equation (43) is obtained.
Here, it is important to note that due to the minimum energy principle, the work of forming of each radial element/section is equal. Channel i (for example, with minimum α) is filled with micro-volumes of metal ∆W i until the work of forming A i (∆W i ) becomes greater than A j (∆W j ) in the other channel j, which fills with the metal until A j (∆W j ) > A m (∆W m ), and therefore, the filling of channel m begins. In this way, all channels are sequentially filled with metal at the equilibrium and the balance of the work conducted and thus, the work of bulk metal formation.
for any i and j.
The master Equation (43) of the developed mathematical model express the work of filling a narrow V-shaped channel by a material under friction following the constant shear Coulomb law with the effective friction coefficient µ i determined by Equation (42). Based on the developed mathematical model, an algorithm of the derivation of friction coefficients µ and m from experimental data has been developed.

Algorithm
The algorithm for calculating the friction coefficients µ and m for the Coulomb constant shear friction law (Algorithm 1) consists of the following steps: The friction coefficients µ and m for the Coulomb constant shear friction law 1.
The length of the extruded radial elements l i ,i = 1 . . . N is measured.

2.
The dimensions of the plastic deformation zone and elongation ratio are calculated: The initial conditions that define zero approximations for µ and m are set: Elongation ratio for radius r pi is calculated using Equation (38): µ 0i is calculated using (42) at [r 1i , r pi ], where friction follows the Coulomb friction law,

6.
The work of A i is calculated using (43). 7.
The average value and the relative variance for the work of filling the radial channels is computed: The

Experiments: Methods, Instruments, and Materials
In the experiments, the dependencies of metal volumes extruded in radial channels on the channels' geometry, length, and dimension ratio were investigated. The model material (led) was extruded into the die equipped with two removable blocks with inserts of eight radial channels. Inserts, shown in Figure 5a,b, provide simultaneous extrusion in eight channels with converging angles varying from 0 to 3.5 • . The forming cavities of the blocks provide stamped forgings with dimension ratios of 1 and 1 2 . The average surface roughness Ra of the forming cavities amounts to 2.5 µm for both blocks. Extrusion of the model material into the forming cavities of each of the two dies was carried out under two different friction boundary conditions, with and without lubricant represented by the oil suspension of graphite powder. Pb was used as a model material because of the lack of significant deformation hardening, which would have introduced an additional error into the experiment.
The experiments were carried out using a hydraulic press PSU-250 (RSCIM, Moscow, Russia) with the maximum force of 250 tons ( Figure 6). The press is equipped with a torsion force meter. The process pressure is displayed on the analogous dual-band dial. The loading module of the press is equipped with two vertical screwed columns and an electrically driven mobile traverse. Extrusion of the model material into the forming cavities of each of the two dies was carried out under two different friction boundary conditions, with and without lubricant represented by the oil suspension of graphite powder. Pb was used as a model material because of the lack of significant deformation hardening, which would have introduced an additional error into the experiment.
The experiments were carried out using a hydraulic press PSU-250 (RSCIM, Moscow, Russia) with the maximum force of 250 tons ( Figure 6). The press is equipped with a torsion force meter. The process pressure is displayed on the analogous dual-band dial. The loading module of the press is equipped with two vertical screwed columns and an electrically driven mobile traverse.
(a) (b) Figure 5. Removable blocks with inserts for the simultaneous extrusion in eight channels with converging angles varying from 0 to 3.5° and dimension ratios of (а) 1 and (b) 1/2.  The die tooling is a block-style die made up of eight inserts ( Figure 5). The block was connected to the die bottom die via a special support. The upper part of the die with the slideway was fixed with eight bolts. The force was transmitted through the upper press plate to the punch, which presses on the billet along the slideway (Figure 7). The die tooling is a block-style die made up of eight inserts ( Figure 5). The block was connected to the die bottom die via a special support. The upper part of the die with the slideway was fixed with eight bolts. The force was transmitted through the upper press plate to the punch, which presses on the billet along the slideway (Figure 7). Measured dimensions of the forgings are summarized in Table 1.  The deforming force, billet height, and depth of filling were measured. The lengths of the extruded sections of forgings were measured in a minimum of seven sections along the height and averaged before being used in calculations of the friction coefficients.

Results and Discussion
The forgings photo and the details of the forgings measurements are given in Figure  8 and Table 2. Measured dimensions of the forgings are summarized in Table 1. The deforming force, billet height, and depth of filling were measured. The lengths of the extruded sections of forgings were measured in a minimum of seven sections along the height and averaged before being used in calculations of the friction coefficients.

Results and Discussion
The forgings photo and the details of the forgings measurements are given in Figure 8 and Table 2.    The friction coefficients were calculated based on the algorithm developed in Section 2. The code for the computations of friction coefficients was developed based on the algorithm of nonlinear models of linear programming [41,42] and implemented in the MatLab 2014b (MathWorks, Natick, MA, USA) environment using the f minimax solver.
The analysis of the computations was carried out accounting for λ l = λ l1 → 0, i.e., r p → r 0 . In this case, metal being extruded into radial channels follows the Coulomb friction law because the size of the constant shear friction zone is very small.
Here, friction during the extrusion is dominated by the constant shear friction law and, thus, µ dominates friction while m is negligible.
In this case, the constant shear friction law with coefficient µ controls the range of [r 0 , r p ] while the Coulomb friction law with coefficient m is dominant in the interval [r p , r 1 ].
The results of the channel-filling simulations for different dimension ratios are shown in Table 3. The dependencies of the relative pressure on the convergence angle for the dimension ratios DR = 1 and DR = 0.5 with and without lubrication are presented in Figure 9. m/μ ≈2 Figure 9. Dependence of the relative pressure at the entrance to the radial channel on the convergence angle.
The experimental data clearly indicated that the extrusion of the material into radial channels with small convergence angles was dominated by the constant shear law because the relative length of the Coulomb friction law zone was negligible (see Table 3). In this case, friction was proportional to the normal contact pressure on the channel walls. For dimension ratio DR = 0.5, the higher number of high values of the relative pressures at the entrance to the radial channels in the case of extrusion with lubricant vs. extrusion without lubricant is explained only by larger volumes of extruded material in the first case (see Figure 9).
The Coulomb friction coefficients calculated for the model material appeared to be independent of the dimension ratio, influenced mostly by roughness and ranging from 0.363 to 0.488. The lubrication provided more stable extrusion conditions as evidenced by the lower relative error than in the case when no lubricant was applied. However, lubrication did not significantly affect the length of the constant shear friction law zone. Contact pressure played a key role here, in agreement with the Coulomb law friction. High values of the friction coefficient can be explained by the specific properties of the model material such as low yield stress and high viscosity.
Recent studies tried to establish an approximate relationship between μ and m [15,16,[43][44][45][46]. Molaei et al. [44], proposed the following relation based on FEM simulations of double cup extrusion and barrel compression tests, which provided better accuracy than some other relationships: Figure 9. Dependence of the relative pressure at the entrance to the radial channel on the convergence angle.
The experimental data clearly indicated that the extrusion of the material into radial channels with small convergence angles was dominated by the constant shear law because the relative length of the Coulomb friction law zone was negligible (see Table 3). In this case, friction was proportional to the normal contact pressure on the channel walls. For dimension ratio DR = 0.5, the higher number of high values of the relative pressures at the entrance to the radial channels in the case of extrusion with lubricant vs. extrusion without lubricant is explained only by larger volumes of extruded material in the first case (see Figure 9).
The Coulomb friction coefficients calculated for the model material appeared to be independent of the dimension ratio, influenced mostly by roughness and ranging from 0.363 to 0.488. The lubrication provided more stable extrusion conditions as evidenced by the lower relative error than in the case when no lubricant was applied. However, lubrication did not significantly affect the length of the constant shear friction law zone. Contact pressure played a key role here, in agreement with the Coulomb law friction. High values of the friction coefficient can be explained by the specific properties of the model material such as low yield stress and high viscosity.
Recent studies tried to establish an approximate relationship between µ and m [15,16,[43][44][45][46]. Molaei et al. [44], proposed the following relation based on FEM simulations of double cup extrusion and barrel compression tests, which provided better accuracy than some other relationships: Dixit and Narayanan [45] obtained the die-pressure distributions in the plane-strain upsetting of a specimen of height h and width w for Coulomb and constant friction models. Based on this study, Dixit et al. [46] after several simplifications, obtained the following result: for (µw/h) 1, m ≈ 2µ [46]. Zhou et al. [39], in FE simulation of sideways extrusion, adopted the Levanov friction model. When normal contact pressure was low ( σ n σ s < 0.5) Levanov's friction stresses were close to those of the Coulomb friction model with a friction factor m ≈ 2µ. Here, m = 0.9 was assumed for the Levanov friction law in the simulation as no lubrication was applied in extrusion practice. Zhang et al. [16] used a load-curve slope in a T-shape compression test to determine the friction factor/coefficient of the constant friction/Coulomb friction model under lubricated conditions and found that m ≈ 2.3µ.
The comparable results of [2,16,[44][45][46] in experiments for the analytical calculation of friction law parameters for extruding metal into narrow V-shaped channels are presented in Table 3, with dependence (32). This analysis showed that it possible to extend the relationship between m and its equivalent µ to the process of filling narrow channels with a large contact area, which, for example, allows one to calculate the forging flash formation when designing finishing dies. The importance of converting these factors to each other is specifically highlighted to introduce the frictional conditions in some professional and commercial finite element software.

Conclusions
An accurate assignment of friction boundary conditions is critically important for FEM analysis of metal forming. The accurate description of friction laws and their coefficients is key for the accurate assignment of the friction boundary conditions. In this work, we developed a mathematical model of extrusion into narrow V-shaped channels with friction following the Coulomb constant shear friction law and an algorithm, implemented in the MatLab environment, for the derivation of the friction coefficient from experimental data. Since friction in the large contact area may differ greatly depending on the conditions of the die surface and changes in the contact pressure, the Coulomb constant shear friction law is better suited for the description of friction during metal forming than the individual Coulomb friction and constant shear friction laws, and thus, the friction coefficients derived from experiments on extrusion into narrow V-shaped channels are relevant directly to real-life manufacturing.
The main results of the present work are given below: • A new mathematical model of coefficients in the Coulomb constant shear friction law for extruding a metal through narrow V-shaped channels with small convergence angles has been developed and evaluated against laboratory measurements. The model is based on an approximate solution of the two-dimensional elasticity problem for equilibrium equations, the kinematically admissible velocity field, and the distribution of tangential stresses in the plastic deformation zone that meet certain boundary conditions. For the zone where a transition from the Coulomb friction law to the constant shear friction law occurs, analytical formulas relating the corresponding friction coefficients have been derived. • An algorithm of the derivation the Coulomb constant shear friction law coefficients from experimental data on the extrusion of a material through narrow V-shaped channels with small convergence angles has been developed.

•
Laboratory experiments on the extrusion of the model material (led) into narrow V-shaped channels with small convergence angles (0-3.5 • ) have been performed. It was found that in this particular case, friction largely obeys the constant friction law, and contact pressure plays a key role in the metal forming. The Coulomb friction coefficients calculated for the model material appear to be independent of the dimension ratio are are influenced mostly by roughness and range from 0.363 to 0.488. The lubrication provides more stable extrusion conditions as evidenced by the lower relative error than in the case when no lubricant is applied. However, lubrication does not significantly affect the length of the constant shear friction law zone.
The developed mathematical model and algorithm can be used to derive friction coefficients, to assign the friction boundary conditions in FEM simulations of a number of processes such as bulk metal forming, extrusion, and drawing, and for the performance tests of lubricants.
A feature of the present study was the flow formation conditions of structural elements of forging with a large ratio of the metal-die interface area to the extruded volume. This circumstance extends the results of the study to forgings with thin ribbon sections or the flash formation for the complete cavity-filling stage in the finishing dies. The friction plays a significant role in these cases and obeys a constant friction law for most narrow V-shaped channels with a large surface area. Determining the transition region from the constant friction to Coulomb's law, which depends on the convergence angle of the Vshaped geometry of the forging flash, allows us to more accurately design dies. In this case, a combined law of friction is set when FE-simulating the metal flow.
It should be noted that the considered laws of friction are of a common nature and, for example, can be applied to solving problems of contact interaction in relation to problems of surface plastic deformation.