Next Article in Journal
Clustering Analysis of Voltage Sag Events Based on Waveform Matching
Previous Article in Journal
Oxidation during Fresh Plant Processing: A Race against Time
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determining the Drawing Force in a Wire Drawing Process Considering an Arbitrary Hardening Law

by
Sergei Alexandrov
1,
Yeong-Maw Hwang
2,* and
Hiu Shan Rachel Tsui
2
1
Federal State Autonomous Educational Institution of Higher Education, South Ural State University (National Research University), 454080 Chelyabinsk, Russia
2
Department of Mechanical and Electro-Mechanical Engineering, National Sun Yat-sen University, Lien-Hai Rd., Kaohsiung 804, Taiwan
*
Author to whom correspondence should be addressed.
Processes 2022, 10(7), 1336; https://doi.org/10.3390/pr10071336
Submission received: 7 June 2022 / Revised: 29 June 2022 / Accepted: 5 July 2022 / Published: 8 July 2022
(This article belongs to the Section Materials Processes)

Abstract

:
The present paper considers the wire drawing process through a conical die and provides an approximate solution for the drawing force, assuming an arbitrary hardening law. The solution is based on the upper bound theorem. However, this theorem does not apply to the stationary flow of strain-hardening materials. Therefore, an adopted engineering approach is to replace the original material model with a non-homogeneous, perfectly plastic model. The kinematically admissible velocity field is derived from an exact semi-analytical solution for the flow through an infinite channel, which increases the accuracy of the solution. The solution for the homogeneous perfectly plastic material compares with an available solution. The general solution is valid for any die angle. The numerical example focuses on the range of angles used in wire drawing. Since the material model is pressure-independent, it is straightforward to adopt the solution for calculating the force in extrusion.

1. Introduction

Wire drawing is a commonly used metal forming process characterized by its high deformation rate and excellent surface quality [1,2,3]. During the process, metal or alloy billets with high tensile strength, such as copper, aluminum, steel, etc., are reshaped through the conical drawing die. The products are widely used in different industry sectors, such as construction, electrical, chemical, and automotive [4,5,6,7]. Papers [8,9,10,11] have pointed out the importance and influence of some wire drawing process parameters, including the reduction ratio, die angle, drawing speed, and friction.
Several approximate analytical methods are available for calculating wire drawing process parameters. One of these methods employs analytical solutions for material flow in infinite converging conical channels. The first solution of this class has been derived in [12]. This rigid perfectly plastic solution is valid for any isotropic pressure-independent yield criterion. Paper [13] has studied some mathematical features of the solution [12] for the von Mises and Tresca yield criteria. Paper [14] has extended the solution [12] to linear hardening materials, assuming the von Mises yield criterion. An approach to using radial flows to analyze drawing processes has been developed in [15]. The corresponding approximate solutions are valid for any hardening law, but small die angle and friction.
Another widely used method for calculating wire drawing process parameters is based on the upper bound theorem. The first solution of this class has been presented in [16]. Other solutions have been derived in [17,18,19]. A distinguishing feature of the solution [17] is that the kinematically admissible velocity field is continuous. Paper [18] has focused on the formation of a rigid region near the friction surface. The singular velocity field in the vicinity of maximum friction surfaces derived in [20] has been employed in [19]. Upper bound solutions are also available for non-conical dies (see, for example, [21,22,23]).
A disadvantage of the upper bound method is that it does not apply to stationary processes for strain-hardening materials. Several papers have calculated the yield stress at the exit from the die by replacing the equivalent strain in the hardening law with the extrusion ratio [24,25,26]. An extension of this method allows average yield stress in each cross-section of the wire to be calculated using geometric parameters [27,28]. This method allows one to replace a strain-hardening material model with the corresponding model of non-homogeneous perfectly plastic material. The upper bound theorem is valid for the latter model. The present paper adopts this approach to develop a new approximate solution for the drawing force. In addition, the velocity field found in [12] is used as a kinematically admissible velocity field in the plastic region. Therefore, the stress equations are satisfied in this region in the case of homogeneous perfectly plastic material, which increases the accuracy of the solution. A comparison with the solution [16] is made.

2. Statement of the Problem

A schematic diagram of the wire drawing process through a conical die is shown in Figure 1. The initial and final radii of the wire are denoted as R 1 and R 2 , respectively. The die angle is 2 φ 0 , and the drawing speed is V 2 . It is required to calculate the drawing force F.
The material is regarded as rigid/plastic, and incompressible. In the case of perfectly plastic material, the von Mises yield criterion reads:
σ e q = σ 0 .
Here σ e q is the equivalent stress and σ 0 is the yield stress in tension. The latter is a constitutive parameter. The friction stress between the conical die and material is taken as:
τ f = m σ 0 3 .
Here m is the constant friction factor. Its value may vary in the range 0 m 1 . A consequence of the incompressibility equation is:
V 1 R 1 2 = V 2 R 2 2 .
Here V 1 is the velocity of the wire before entering the die.
The solution is based on the upper bound theorem [29]. Therefore, strictly speaking, strain-hardening models cannot be considered because the process is stationary. The von Mises yield criterion of strain-hardening materials is represented as:
σ e q = σ 0 Φ ( ε e q )
where ε e q is the equivalent strain and Φ ( ε e q ) is an arbitrary function of its argument satisfying the conditions Φ ( 0 ) = 1 and d Φ ( ε e q ) d ε e q 0 for all ε e q . An engineering approach to account for strain hardening has been proposed in [28] for linear strain-hardening materials. According to this approach, the equivalent strain in the hardening law is replaced with an average strain defined as:
ε a = l n ( A 0 A )
where A 0 is the cross-sectional area of the wire before entering the die and A is the cross-sectional area of the wire at any distance from the entrance. It is evident from (5) that ε a = 0 at A = A 0 . Therefore, it follows from (4) that σ e q = σ 0 at the entrance to the die. Paper [28] assumes that Φ ( ε e q ) is a linear function of its argument. However, it is straightforward to extend this approach to any function. Then, Equation (4) becomes:
σ e q = σ 0 Φ ( ε a ) .
Since ε a is a known function of position, the corresponding model can be regarded as rigid perfectly plastic with a known nonuniform yield stress distribution. The upper bound theorem applies to such a model. Equation (2) becomes:
τ f = m Φ ( ε a ) 3 .  

3. Kinematically Admissible Velocity Field

It is natural to use a spherical coordinate system ( r , φ , ϕ ) . The solution is independent of ϕ . The shear stress in the spherical coordinate system is denoted as σ r φ . The equation of the die surface is φ = φ 0 (Figure 2). The boundary conditions (2) and (7) become:
σ r φ = m σ 0 3   and   σ r φ = m σ 0 Φ ( ε a ) 3
for φ = φ 0 , respectively. Axial symmetry demands
σ r φ = 0
for φ = 0 . The spherical coordinate system and the general structure of the kinematically admissible velocity field are shown in Figure 2. The kinematically admissible velocity field consists of a plastic region and two rigid regions. The rigid/plastic boundaries are velocity discontinuity lines in meridian planes.

3.1. Kinematically Admissible Velocity Field in the Plastic Region

Choosing a kinematically admissible velocity field that satisfies the stress equations when combined with the associated flow rule is advantageous. Semi-analytical solutions to boundary value problems supply such velocity fields. In particular, the solution presented in [12] can be used in the case under consideration. This solution is valid for any pressure-independent yield criterion. The present paper employs the solution for the von Mises yield criterion. The dimensionless shear stress τ = 3 σ r φ σ 0 satisfies the following Equation [12]:
d τ d φ + τ c o t φ + 2 3 1 τ 2 = c .
Here c is constant. Its value should be found from the solution. It is convenient to put:
τ = s i n ψ   and   1 τ 2 = c o s ψ .
Equations (10) and (11) combine to give:
d ψ d φ = c s i n ψ c o t φ 2 3 c o s ψ c o s ψ .
The first equation in (8) and (9) become τ = m for φ = φ 0 and τ = 0 for φ = 0 , respectively. Using (11), one can transform these boundary conditions to:
ψ = 0
for φ = 0 and:
ψ = ψ f = a r c s i n   m
for φ = φ 0 . It is seen from (13) that the term s i n ψ c o t φ involved in (12) reduces to the expression 0 0 at φ = 0 . Assume that ψ is a linear function of φ in the vicinity of φ = 0 . This function must satisfy (13). Substituting it into (12) and expanding the right-hand side of the resulting equation in a power series, one obtains
ψ = ( c 2 3 ) 2 φ + o ( φ )
as φ 0 . Equation (12) should be solved numerically using (15). The solution depends on c. The boundary condition (14) supplies the equation to determine c. Figure 3 illustrates the dependence of c on m and φ 0 . In what follows, it is assumed that ψ is a known function of φ .
The only non-vanishing velocity component is the radial velocity. The latter is given by [12]
u = β V R 2 r 2 e x p [ 2 3 φ 0 φ τ 1 τ 2 d γ ] .
Here β > 0 is constant whose value should be found from the solution. In (16) and what follows, the quantity V R 2 may be represented as V 1 R 1 2 or V 2 R 2 2 due to (3). Using (11), one transforms (16) to
u = β V R 2 r 2 q ( φ )
where
q ( φ ) = e x p [ 2 3 φ 0 φ t a n ψ d γ ] .
A consequence of the incompressibility equation is:
π V R 2 = 2 π 0 φ 0 u r 2 s i n φ d φ .
Eliminating here u using (17), one arrives at the following equation for β :
1 = 2 β 0 φ 0 q ( φ ) s i n φ d φ .
In what follows, it is assumed that β has been calculated. Therefore, the kinematically admissible velocity field in the plastic region has been determined.

3.2. Velocity Discontinuity Lines

There are the two velocity discontinuity lines shown in Figure 4 (BB1 and AA1). Two orthogonal systems of unit vectors associated with line BB1 are depicted in this figure. A similar system of vectors may be associated with line AA1. The derivation below is valid for both lines.
The base vectors of the spherical coordinate system are denoted as e r and e φ . The vectors directed along the normal and tangent to the velocity discontinuity line are denoted as n and τ , respectively. Moreover, i is the unit vector directed along the axis of symmetry. It follows from the geometry of Figure 4 that the relationships between the unit vectors are:
n = e r c o s γ + e φ s i n γ ,   τ = e r s i n γ + e φ c o s γ ,   i = e r c o s φ e φ s i n φ .
Here γ is the inclination of the vector n to the r-axis, measured clockwise. The velocity vectors in the plastic and rigid regions are represented as:
U p = u e r   and   U r = V i ,  
respectively. Here, and in what follows, V should be replaced with V1 for line BB1 and with V2 for line AA1. The normal velocity must be continuous across the velocity discontinuity line. Therefore, U p n = U r n . Upon substitution from (21) and (22), this equation becomes:
u V = c o s ( γ + φ ) c o s γ .
The amount of velocity jump across the velocity discontinuity line is [ u ] = U r τ U p τ . Upon substitution from (21) and (22), this equation becomes:
[ u ] V = u V s i n γ + s i n ( γ + φ ) .
Eliminating here u using (23), one obtains:
[ u ] V = s i n φ c o s γ .
An infinitesimal length element of the velocity discontinuity line is:
d l = ( d r ) 2 + ( r d φ ) 2 = r ( d r r d φ ) 2 + 1 d φ .
It follows from the geometry of Figure 4 that:
t a n γ = d r r d φ .
Equations (26) and (27) combine to give:
d l = r c o s γ d φ .
Using (17), one can transform Equation (23) to:
β R 2 r 2 q ( φ ) = c o s ( γ + φ ) c o s γ .
Employing trigonometric identities and (27), one rewrites this equation as:
d s d φ = 2 s i n φ [ β q ( φ ) s c o s φ ]
where s = r 2 R 2 . Here R = R 1 for line BB1 and R = R 2 for line AA1. Equation (30) determines the velocity discontinuity line. This line must pass through point B or A (Figure 4). Therefore, the boundary condition to Equation (30) is:
s = 1 s i n 2 φ 0
for φ = φ 0 . The general solution of Equation (30) is:
s = 1 s i n 2 φ [ 2 β 0 φ s i n χ q ( χ ) d χ + s 0 ] .
Here s 0 is constant. The value of s must be finite at φ = 0 . Otherwise, the velocity discontinuity line does not reach the symmetry axis. Therefore, it is necessary to put s 0 = 0 . Then, Equation (32) becomes:
s = 2 β s i n 2 φ S ( φ ) ,   S ( φ ) = 0 φ s i n χ q ( χ ) d χ .
This solution reduces to the expression 0 0 at φ = 0 . Applying L’Hospital’s rule leads to:
s = β q 0
at φ = φ 0 . Here q 0 is the value of q at φ = 0 . It follows from (18) that:
q 0 = e x p [ 2 3 0 φ 0 t a n ψ d γ ] .
It is seen from (20) and (33) that (31) is satisfied.
An infinitesimal area element of the velocity discontinuity surface is determined from (28) as:
d Σ = r 2 s i n φ c o s γ d φ d ϕ .
Using the definition for s and (33), one transforms this equation to:
d Σ = 2 β R 2 S ( φ ) c o s γ s i n φ d φ d ϕ .

4. Plastic Work Rate

It is necessary to calculate the plastic work rate (i) in the plastic region, (ii) at two velocity discontinuity surfaces, and (iii) at the friction surface.

4.1. Plastic Work Rate in the Plastic Region

The infinitesimal volume element is d Ω = r 2 s i n φ d r d φ d ϕ . Using the definition for s, one transforms this equation to:
d Ω = R 3 s s i n φ 2 d s d φ d ϕ .
The plastic work rate is defined as:
W Ω = Ω σ e q ξ e q d Ω .
Here ξ e q is the equivalent plastic strain rate. In the case under consideration, it is given by:
ξ e q = 2 3 ξ r r 2 + ξ φ φ 2 + ξ ϕ ϕ 2 + 2 ξ r φ 2 .
The strain rate components are expressed through the radial velocity as:
ξ r r = u r ,   ξ φ φ = ξ ϕ ϕ = u r , ξ r φ = 1 2 r u φ .
Equations (17), (18) and (41) combine to give:
ξ r r = 2 β V R 2 r 3 q ( φ ) ,   ξ φ φ = ξ ϕ ϕ = β V R 2 r 3 q ( φ ) ,   ξ r φ = 3 β V R 2 r 3 q ( φ ) t a n ψ .
Substituting (41) into (40) and using the definition for s yields:
ξ e q = 2 β V q ( φ ) R s s c o s ψ .
Equations (38), (39) and (43) combine to give:
w Ω = W Ω π V R 2 σ 0 = 2 β 0 φ 0 s 1 ( φ ) s 2 ( φ ) σ e q σ 0 q ( φ ) s i n φ s c o s ψ d s d φ .
This equation involves both velocity discontinuity lines. Therefore, it is necessary to adopt the same definition for s. Assuming that R = R 2 and using (33) results in:
s 1 ( φ ) = 2 β s i n 2 φ 0 φ s i n χ q ( χ ) d χ   and   s 2 ( φ ) = 2 β s i n 2 φ ( R 1 R 2 ) 2 0 φ s i n χ q ( χ ) d χ .
The first equation determines the velocity discontinuity line through A and the second through B (Figure 4). The equations in (45) allow the integral in (44) to be evaluated numerically. In the case of homogeneous properties, Equation (1) is valid, and Equations (44) and (45) yield:
w Ω = 4 β l n ( R 1 R 2 ) 0 φ 0 q ( φ ) s i n φ c o s ψ d φ .
In the case of non-homogeneous properties, Equation (6) is valid, and Equation (44) yields:
w Ω = 2 β 0 φ 0 s 1 ( φ ) s 2 ( φ ) Φ ( ε a ) q ( φ ) s i n φ s c o s ψ d s d φ .
In the spherical coordinate system, Equation (5) becomes:
ε a = 2 l n ( R 1 r s i n φ 0 ) .
Using the definition for s, one can transform this equation to:
ε a = 2 l n ( R 1 R 2 s i n φ 0 s ) .
Eliminating ε a in (47) using (49), one represents the integrand as a function of s and φ .

4.2. Plastic Work Rate at the Velocity Discontinuity Surfaces

The plastic work rate at each velocity discontinuity surface is given by:
W d = 1 3 Σ σ e q [ u ] d Σ .
Using (25) and (37), one transforms this equation to:
w d = W d π V R 2 σ 0 = 4 β 3 0 φ 0 σ e q σ 0 c o s 2 γ S ( φ ) d φ .
This equation is valid for both velocity discontinuity surfaces. In the case of homogeneous properties, Equation (51) becomes:
w d = 4 β 3 0 φ 0 S ( φ ) c o s 2 γ d φ
for each velocity discontinuity surface. In the case of non-homogeneous properties, Equation (51) becomes:
w d = 4 β 3 0 φ 0 Φ ( ε a ) c o s 2 γ S ( φ ) d φ .
Here ε a should be regarded as a function of φ . One derives this function by replacing s in (49) with s 1 ( φ ) or s 2 ( φ ) . Therefore, the plastic work rate is not the same at the two velocity discontinuity surfaces.
Equations (52) and (53) involve c o s 2 γ . It follows from (27) and (30) that:
t a n γ = s c o s φ β q ( φ ) s s i n φ .
The right-hand side of this equation is a known function of φ due to (33). Therefore, c o s 2 γ is determined as a function of φ using (54) and the identity c o s 2 γ = ( 1 + t a n 2 γ ) 1 .

4.3. Plastic Work Rate at the Friction Surface

The radial velocity at the friction surface is determined from (17) and (18) as:
u f = β V R 2 r 2 .
The plastic work rate is:
W f = 2 π m 3 R 2 s i n φ 0 R 1 s i n φ 0 σ e q | u f | r s i n φ 0 d r .
Equations (55) and (56) combine to give:
w f = W f π V R 2 σ 0 = 2 m β s i n φ 0 3 R 2 s i n φ 0 R 1 s i n φ 0 σ e q σ 0 d r r .
In the case of homogeneous properties, this equation becomes:
w f = 2 m β s i n φ 0 3 l n ( R 1 R 2 ) .
In the case of non-homogeneous properties, Equation (57) becomes:
w f = 2 m β s i n φ 0 3 R 2 s i n φ 0 R 1 s i n φ 0 Φ ( ε a ) r d r .
Here one should eliminate ε a using (48).

5. Drawing Force

Let F be the drawing force. It follows from the upper bound theorem that:
F V 2 W Ω + W d 1 + W d 2 + W f .
Here W d 1 is W d calculated at the velocity discontinuity line through A, and W d 2 is W d calculated at the velocity discontinuity line through B (Figure 4). Using (44), (51) and (57), one transforms Equation (60) to:
f u = F u π R 2 2 σ 0 = w Ω + w d 1 + w d 2 + w f .
Here F u is the upper bound on F and f u is its dimensionless representation. The right-hand side of Equation (61) has been calculated in the previous section for both homogeneous and non-homogeneous materials. Quantitative results can be obtained after prescribing the function Φ ( ε e q ) and evaluating the integrals numerically.

6. Numerical Examples

The right-hand side of Equation (61) has been calculated using the corresponding build-in commands in Wolfram Mathematica (version 11.3). The present section illustrates these results and emphasizes the effect of the reduction ratio, the die angle, and the friction factor on the drawing force. Even though the general solution is valid for any values of these parameters, the quantitative results focus on their typical ranges in the drawing of thin wires. In the case of strain-hardening material, the stress–strain curve of AISI-316 stainless steel is approximated by the equation:
Φ ( ε e q ) = ( 1 + 84.75 ε e q ) 0.52 .
The constitutive parameter σ 0 involved in (4) and required to calculate the drawing force F after determining f is σ 0 = 119 MPa.
Figure 5 compares the present solution for rigid perfectly plastic material and the solution provided in [16]. The former is shown by the solid lines, and the latter by the broken lines. The difference is invisible in the figure. It is worthy of note that the solution [16] has been specifically derived for small die angles. The present solution is more general and is valid for any die angle. Figure 5 demonstrates that this advantageous feature of the solution does not affect its accuracy.
Figure 6, Figure 7, Figure 8 and Figure 9 show the effect of process parameters on the dimensionless drawing force. Each figure shows the effect of the die angle for several friction factor values. The effect of the reduction ratio is revealed by comparing the results depicted in all these figures. The effect of strain hardening is revealed by comparing the results depicted in Figure 5 and Figure 6. All results coincide with physical expectations. In particular, the processing of strain-hardening materials requires a larger force than materials with no hardening. In the case under consideration, this difference is considerable. It is because of the hardening law (62). It is seen from each figure that an optimal die angle (i.e., the drawing force attains a minimum at this angle) exists if the other process parameters are kept constant. Its value depends on these other process parameters and may be outside the typical ranges considered. The drawing force increases as the friction factor increases.

7. Discussion

The present paper has developed an engineering approach that allows one to use the upper bound technique for analyzing the stationary metal forming processes for strain-hardening materials. This development is important because the upper bound solutions are in good agreement with finite element solutions when the upper bound theorem is valid. It has been demonstrated in [30] in the case of wire drawing processes. It is unsurprising because the upper bound theorem provides an accurate result, even from comparatively crude kinematically admissible velocity fields [31]. The present solution is based on the kinematically admissible velocity field that, in conjunction with the associated flow rule, satisfies all the fundamental and constitutive equations in the plastic region in the case of rigid perfectly plastic material. This kinematically admissible velocity field feature increases the solution’s accuracy.
The solution found shows the effect of process parameters on the drawing force. It is straightforward to use this solution for other applications. Examples of such applications are the prediction of central bursting defects using the approaches proposed in [32,33], and the prediction of fracture initiation using the extended Bernoulli’s theorem [34] in conjunction with uncoupled ductile fracture criteria. A review of such criteria is available, for example, in [35].

8. Conclusions

The process of wire drawing through a conical die has been investigated. A new theoretical solution for the drawing force has been derived, assuming an arbitrary strain- hardening law. The solution is based on the upper bound theorem, although it is not an upper bound solution for the standard model of rigid plastic strain-hardening material. The reason for it is that the theorem does not apply to the stationary flow of rigid plastic strain-hardening materials. An engineering approach has been used to overcome this difficulty. In particular, the equivalent strain in the strain-hardening law is replaced with an average strain in each cross-section of the wire. The latter is calculated using geometric parameters of the die. As a result, the original material model is replaced with the rigid perfectly plastic model with a non-homogeneous yield stress distribution. The upper bound theorem applies to such a model.
The solution has been reduced to several ordinary integrals that should be evaluated numerically. Section 6 illustrates the effect of process parameters on the dimensionless drawing force, assuming the hardening law (62). The behavior of the solution agrees with physical expectations. The solution practically coincides with solution [16] (Figure 5). The latter has been developed for small die angles. The new solution is valid for any die angle. An advantage of this solution is that the kinematically admissible velocity field is taken from the exact solution for material flow through an infinite channel [12].
Since the material is pressure-independent, it is straightforward to adopt the solution provided for calculating the extrusion force. To this end, it is necessary to add the same hydrostatic stress to all normal stresses such that the drawing force vanishes. Front and back forces can be included similarly.

Author Contributions

Conceptualization, S.A. and Y.-M.H.; formal analysis, H.S.R.T.; writing—review and editing, S.A. and Y.-M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of the Republic of China, grant number MOST 108-2923-E-110-002-MY3.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, Y.C.; Kim, B.M. A study on the corner filling in the drawing of a rectangular rod from a round bar. Int. J. Mach. Tools Manuf. 2000, 40, 2099–2117. [Google Scholar] [CrossRef]
  2. Felder, E.; Levrau, C.; Mantel, M.; Truong Dinh, N.G. Identification of the work of plastic deformation and the friction shear stress in wire drawing. Wear 2012, 286–287, 27–34. [Google Scholar] [CrossRef]
  3. Martinez, G.A.S.; Qian, W.L.; Kabayama, L.K.; Prisco, U. Effect of Process Parameters in Copper-Wire Drawing. Metals 2020, 10, 105. [Google Scholar] [CrossRef] [Green Version]
  4. Vega, G.; Haddi, A.; Imad, A. Investigation of process parameters effect on the copper-wire drawing. Mater. Des. 2009, 30, 3308–3312. [Google Scholar] [CrossRef]
  5. Bitkov, V. Nonuniformity of deformation in drawing corrosionresistant steel and alloy wire. Steel Transl. 2010, 40, 169–173. [Google Scholar] [CrossRef]
  6. Huang, X.; Wang, J.; Zhao, S.; Yao, Z.; Liu, C. High-resolution multiscale modeling of mechanical behavior of cold-drawn pearlitic steels. J. Mater. 2021, 15, 5920–5935. [Google Scholar] [CrossRef]
  7. Liu, X.Y.; Zhang, S.H. The design of a drawing die based on the logistic function for the energy analysis of drawing force. Appl. Math. Model. 2022, 109, 833–847. [Google Scholar] [CrossRef]
  8. Dixit, U.S.; Dixit, P.M. An analysis of the steady-state wire drawing of strain hardening materials. J. Mater. Process Technol. 1995, 47, 201–229. [Google Scholar] [CrossRef]
  9. Cetlin, P.R.; Correa, E.C.S.; Aguilar, M.T.P. The effect of strain path on the work hardening of austenitic and ferritic stainless steels in axi-symmetric drawing. Metall. Mater. Trans. A 2003, 34, 589–601. [Google Scholar] [CrossRef]
  10. Haddi, A.; Imad, A.; Vega, G. Analysis of temperature and speed effects on the drawing stress for improving the wire drawing process. Mater. Des. 2011, 32, 4310–4315. [Google Scholar] [CrossRef]
  11. Tintelecan, M.; Sas-Boca, I.M.; Ilutiu-Varvara, D.A. The influence of the dies geometry on the drawing force for steel wires. Procedia Eng. 2017, 181, 193–199. [Google Scholar] [CrossRef]
  12. Shield, R.T. Plastic flow in a converging conical channel. J. Mech. Phys. Solids 1955, 3, 246–258. [Google Scholar] [CrossRef]
  13. Bonitskaya, O.V.; Krasavin, R.V.; Markin, A.A. The steady axisymmetric flow of ideally plastic materials in a conical channel. J. Appl. Math. Mech. 2010, 74, 494–500. [Google Scholar] [CrossRef]
  14. Durban, D. Axially symmetric radial flow of rigid/linear-hardening materials. J. Appl. Mech. 1979, 46, 322–328. [Google Scholar] [CrossRef]
  15. Durban, D. Radial flow simulation of drawing and extrusion of rigid/hardening materials. Int. J. Mech. Sci. 1983, 25, 27–39. [Google Scholar] [CrossRef]
  16. Avitzur, B. Analysis of wire drawing and extrusion through conical dies of small cone angle. J. Eng. Ind. 1963, 85, 89–95. [Google Scholar] [CrossRef]
  17. Lambert, E.R.; Kobayashi, S. A theory on the mechanics of axisymmetric extrusion through conical dies. J. Mech. Eng. Sci. 1968, 10, 367–380. [Google Scholar] [CrossRef]
  18. Tirosh, J. On the dead-zone formation in plastic axially-symmetric converging flow. J. Mech. Phys. Solids 1971, 19, 39–47. [Google Scholar] [CrossRef]
  19. Alexandrov, S.; Mishuris, G.; Mishuris, W.; Sliwa, R.E. On the dead zone formation and limit analysis in axially symmetric extrusion. Int. J. Mech. Sci. 2001, 43, 367–379. [Google Scholar] [CrossRef]
  20. Alexandrov, S.; Richmond, O. Singular plastic flow fields near surfaces of maximum friction stress. Int. J. Non-Linear Mech. 2001, 36, 1–11. [Google Scholar] [CrossRef]
  21. Gordon, W.A.; Van Tyne, C.J.; Sriram, S. Extrusion through spherical dies—An upper bound analysis. J. Manuf. Sci. Eng. 2002, 124, 92–97. [Google Scholar] [CrossRef]
  22. Chang, D.-F.; Wang, J. Optimized upper bound analysis of axisymmetric extrusion using spherical velocity field. J. Manuf. Sci. Eng. 2006, 128, 4–10. [Google Scholar] [CrossRef]
  23. Zhang, S.H.; Chen, X.D.; Zhou, J.; Zhao, D.W. Upper bound analysis of wire drawing through a twin parabolic die. Meccanica 2016, 51, 2099–2110. [Google Scholar] [CrossRef]
  24. Gur’yanov, G.N. Effect of the deformation parameters on the forces in the tool during wire-drawing. Metallurgist 2008, 52, 197–201. [Google Scholar] [CrossRef]
  25. Gur’yanov, G.N. Evaluating the maximum draft and safety factor in wire-drawing. Metallurgist 2008, 52, 582–587. [Google Scholar] [CrossRef]
  26. Gur’yanov, G.N.; Zhelezkov, O.S.; Platov, S.I.; Terent’ev, D.V. Influence of the drawing angle on wire quality. Steel Transl. 2015, 45, 18–24. [Google Scholar] [CrossRef]
  27. Kazeminezhad, M. Combination of the upper bound and potts models for simulation of microstructure in wire drawing and annealing processes. J. Mater. Eng. Perform. 2009, 18, 26–31. [Google Scholar] [CrossRef]
  28. Panteghini, A.; Genna, F. An engineering analytical approach to the design of cold wire drawing processes for strain-hardening materials. Int. J. Mater. Form. 2010, 3, 279–289. [Google Scholar] [CrossRef]
  29. Hill, R. The Mathematical Theory of Plasticity; Oxford University Press: New York, NY, USA, 1950. [Google Scholar]
  30. Luis, C.J.; Leon, J.; Luri, R. Comparison between finite element method and analytical methods for studying wire drawing processes. J. Mater. Process. Technol. 2005, 164–165, 1218–1225. [Google Scholar] [CrossRef]
  31. Hill, R. New horizons in the mechanics of solids. J. Mech. Phys. Solids 1956, 5, 66–74. [Google Scholar] [CrossRef]
  32. Avitzur, B. Analysis of central bursting defects in extrusion and wire Drawing. J. Eng. Ind. 1968, 90, 79–90. [Google Scholar] [CrossRef]
  33. Parghazeh, A.; Haghighat, H. Prediction of central bursting defects in rod extrusion process with upper bound analysis method. Trans. Nonferr. Metal. Soc. 2016, 26, 2892–2899. [Google Scholar] [CrossRef]
  34. Richmond, O.; Alexandrov, S. Extension of Bernoulli’s theorem on steady flows of inviscid fluids to steady flows of plastic solids. C. R. L’académie Sci.-Ser. IIB–Mech. 2000, 328, 835–840. [Google Scholar] [CrossRef]
  35. Atkins, A.G. Fracture in forming. J. Mater. Process. Technol. 1996, 56, 609–618. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the process.
Figure 1. Schematic diagram of the process.
Processes 10 01336 g001
Figure 2. General structure of the kinematically admissible velocity field.
Figure 2. General structure of the kinematically admissible velocity field.
Processes 10 01336 g002
Figure 3. Variation of c with the die semi-angle at several friction factors.
Figure 3. Variation of c with the die semi-angle at several friction factors.
Processes 10 01336 g003
Figure 4. Geometry of velocity discontinuity lines.
Figure 4. Geometry of velocity discontinuity lines.
Processes 10 01336 g004
Figure 5. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.07 and several friction factors for perfectly plastic material. The solid curves correspond to the present solution and the broken curves to the solution [16].
Figure 5. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.07 and several friction factors for perfectly plastic material. The solid curves correspond to the present solution and the broken curves to the solution [16].
Processes 10 01336 g005
Figure 6. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.07 and several friction factors for strain-hardening material.
Figure 6. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.07 and several friction factors for strain-hardening material.
Processes 10 01336 g006
Figure 7. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.1 and several friction factors for strain-hardening material.
Figure 7. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.1 and several friction factors for strain-hardening material.
Processes 10 01336 g007
Figure 8. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.13 and several friction factors for strain-hardening material.
Figure 8. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.13 and several friction factors for strain-hardening material.
Processes 10 01336 g008
Figure 9. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.16 and several friction factors for strain-hardening material.
Figure 9. Variation of the dimensionless drawing force with the die semi-angle at R 1 R 2 = 1.16 and several friction factors for strain-hardening material.
Processes 10 01336 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alexandrov, S.; Hwang, Y.-M.; Tsui, H.S.R. Determining the Drawing Force in a Wire Drawing Process Considering an Arbitrary Hardening Law. Processes 2022, 10, 1336. https://doi.org/10.3390/pr10071336

AMA Style

Alexandrov S, Hwang Y-M, Tsui HSR. Determining the Drawing Force in a Wire Drawing Process Considering an Arbitrary Hardening Law. Processes. 2022; 10(7):1336. https://doi.org/10.3390/pr10071336

Chicago/Turabian Style

Alexandrov, Sergei, Yeong-Maw Hwang, and Hiu Shan Rachel Tsui. 2022. "Determining the Drawing Force in a Wire Drawing Process Considering an Arbitrary Hardening Law" Processes 10, no. 7: 1336. https://doi.org/10.3390/pr10071336

APA Style

Alexandrov, S., Hwang, Y. -M., & Tsui, H. S. R. (2022). Determining the Drawing Force in a Wire Drawing Process Considering an Arbitrary Hardening Law. Processes, 10(7), 1336. https://doi.org/10.3390/pr10071336

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop