Next Article in Journal
Pilot-Scale Anaerobic Digestion of Pig Manure with Thermal Pretreatment: Stability Monitoring to Improve the Potential for Obtaining Methane
Previous Article in Journal
A Metallic Fracture Estimation Method Using Digital Image Correlation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ideal Flow Design of Plane-Strain Bending Driven by Springback

by
Sergei Alexandrov
1,2,
Marina Rynkovskaya
2 and
Yeong-Maw Hwang
3,*
1
Ishlinsky Institute for Problems in Mechanics RAS, 101-1 Prospect Vernadskogo, 119526 Moscow, Russia
2
Department of Civil Engineering, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St., 117198 Moscow, Russia
3
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(8), 1601; https://doi.org/10.3390/pr10081601
Submission received: 16 July 2022 / Revised: 7 August 2022 / Accepted: 9 August 2022 / Published: 13 August 2022
(This article belongs to the Section Materials Processes)

Abstract

:
The present paper is concerned with the design of the plane-strain bending of wide sheets driven by the sheet’s geometric shape after unloading. The solution is based on the ideal flow theory. It determines the loading path (the dependence between the dimensionless bending moment and the dimensionless tensile force) that produces the desirable shape. The ideal flow theory is valid for the rigid perfectly plastic material model. A method proposed in the literature combines a rigid plastic solution at loading and an elastic solution for springback. The final design solution is practically analytical. A numerical method is only necessary to solve a system of transcendental equations. An illustrative example is provided. This example shows that, in some cases, it is sufficient to use the solution at loading to determine an accurate loading path that produces the desirable sheet’s shape after springback.

1. Introduction

All sheet metal forming processes incorporate some bending. There is a vast amount of literature on springback for such processes. A comprehensive review of several issues related to springback was provided in [1]. Most theoretical studies have focused on the springback predictions using analytical and numerical methods. Several works have proposed procedures for reducing springback [2,3,4,5,6]. Paper [2] recommended reheating before post-stretching to reduce the elastic recovery. Optimal temperature, punch speed, holding time, and sheet orientation relative to the rolling direction were found in [3] for anisotropic sheets of Inconel 625 alloy. The effect of different heating methods on springback in the dieless bending of mild steel sheets has been studied in [4,5]. Paper [6] has employed regression analysis to develop an analytical model for optimizing the V-bending process of JSC-590 sheets. The present paper provides a design procedure that aims to achieve the desired shape after four-point bending followed by unloading. The design is based on the ideal flow theory [7].
The ideal flow theory was invented in [8], as the planar stationary flow of rigid perfectly plastic material was considered. The theory has been extended to three-dimensional stationary and non-stationary flows in [9]. It is used for the preliminary design of metal forming processes [7]. However, all available non-trivial bulk ideal flow theory applications are restricted to stationary flows. A plane strain strip drawing design was proposed in [8]. A more general solution for an arbitrary flow with a symmetry axis was derived in [10]. An axisymmetric drawing design was proposed in [11]. A review of ideal flow design solutions for drawing and extrusion has been provided in [12]. All these solutions determine a die shape that produces an ideal flow. Additionally, solution [11] supplies the die of minimum length.
The present paper combines the ideal flow theory for non-stationary planar processes and the method for analyzing plane-strain bending proposed in [13]. A review of the latter is presented in [14]. The design objective is to determine an ideal flow loading path that results in desirable geometric parameters after unloading. It is worthy of note that the ideal flow theory is based on a rigid plastic model. Therefore, the approach proposed in [15] is adopted for calculating springback. It has been shown in [16] that the effect of neglecting elasticity on loading on springback in bending is negligible. It is worthy of note that the new design disposes of the reversed loading for fibers in a central part of the sheet. Therefore, the Bauschinger effect, which may affect springback [17], is immaterial.
The solution at loading is purely analytical. Assuming that unloading is purely elastic, paper [18] derives a general solution for the final shape. This solution is used in conjunction with the new solution developed in the present paper.

2. Statement of the Problem

A wide sheet of initial thickness H and initial width 2L is subject to bending under tension. Its initial and intermediate/final configurations are shown in Figure 1a,b, respectively. The bending moment is denoted as M and the tensile force as F. Both are in per unit length. The current thickness, radius of AB, and radius of CD are denoted as h, rAB, and rCD, respectively. It is required to find a loading path in the MF space such that
r C D = r f and h = h f
after unloading. Additionally, the process must be an ideal flow process. It is assumed that the material obeys the von Mises yield criterion and its associated flow rule.
The general concept of ideal flow design was described in detail in [7]. The present section applies this concept to the process of plane-strain bending under tension. The general structure of the rigid/plastic solution for this process is illustrated in Figure 2a. All fibers in Region 1 are extended, and all fibers in Region 2 are compressed in the circumferential direction throughout the deformation process. All fibers in Region 3 are first compressed and afterward extended in the circumferential direction. Region 3 does not exist at the initial instant. The ideal flow condition is satisfied in Regions 1 and 2. The boundary between Regions 2 and 3 propagates to surface CD and moves through the material. For this reason, the ideal flow condition is not satisfied in Region 3. Therefore, the ideal flow design should determine such a dependence between the bending moment and the tensile force that the initial boundary between Regions 1 and 2 does not move through the material. Thus, Region 3 does not exist in the ideal flow process of bending under tension throughout the deformation process (Figure 2b).
It is worthy of note that Region 2 does not exist at the initial instant if the tensile force is large enough. Such processes will be considered separately.
The ideal plastic solutions are rigid/plastic. Nevertheless, the determination of residual stresses and springback is possible using the method proposed in [15].

3. Solution at Loading

3.1. Kinematics

With no loss of generality, one can choose the origin of the Cartesian coordinate system (x, y) shown in Figure 2b at the intersection of the axis of symmetry of the process and the outer surface of the sheet. Then, the process is symmetric relative to the x-axis. The kinematics of plane strain bending of incompressible material was described in [13]. This description is independent of the constitutive equations. It is also valid for plane strain bending under tension [18]. It follows from [13] that
x H = ζ a + s a 2 cos ( 2 a η ) s a and y H = ζ a + s a 2 sin ( 2 a η ) .
Here, (ζ, η) is the Lagrangian coordinate system satisfying the condition ζ = x / H and η = y / H at the initial instant, a is a time-like parameter, and s is a function of a. This function should be found from the solution. At the initial instant:
a = 0 and s = 1 4 .
It is also convenient to introduce a polar coordinate system ( r ,   θ ) by the transformation equations:
r cos θ H = x H + s a and r sin θ H = y H .
Then, r = r A B on AB, r = r C D on CD, θ = θ 0 on BC, and θ = θ 0 on AD (Figure 2b). Equations (2) and (4) combine to give
r H = ζ a + s a 2 and θ = 2 a η .
It follows from the definition for the Lagrangian coordinates that ζ = 0 on AB, ζ = 1 on CD, η = L / H on BC, and η = L / H on AD (Figure 2b). Then, Equation (5) supplies
r A B H = s a , r C D H = s a 2 1 a and h H = r A B r C D H = s s a a .
The equations in (2) show that the principal strain rate trajectories coincide with the coordinate curves of the Lagrangian coordinate system. Moreover, the quadratic invariant of the strain rate tensor is determined from these equations as
ξ e q = | ζ + d s / d a | 3 ( ζ a + s ) d a d t .
Then, ξ e q = 0 at
ζ = ζ 0 = d s d a .
This equation is valid at the boundary between Regions 1 and 2 at the initial instant (Figure 2b).
One should distinguish two cases: ζ 0 1 and ζ 0 < 1 at the initial instant. Since the coordinate ζ is Lagrangian, the ideal flow condition formulated in Section 2 requires that ζ 0 is independent of a if ζ 0 1 at the initial instant. The line ζ = ζ 0 lies outside the domain occupied by the material if ζ 0 < 1 . Therefore, it is not necessary that ζ 0 is independent of a, though it is important to verify that the condition ζ 0 1 is satisfied throughout the deformation process. However, the condition that ζ 0 is independent of a allows one to determine an ideal flow process if ζ 0 < 1 . This process is considered below. Then, independently of the value of ζ 0 at the initial instant, solving (8) together with the initial condition (3) yields
s = 1 4 ζ 0 a .
Let a l , h l , r A B ( l ) and r C D ( l ) be the values of a, h, r A B and r C D at the end of loading, respectively. Then, it follows from (6) and (9) that
h l H = 1 4 ζ 0 a l 1 4 a l ( ζ 0 + 1 ) 2 a l , r A B ( l ) H = 1 4 a l ( ζ 0 + 1 ) 2 a l , r C D ( l ) H = 1 4 a l ( ζ 0 + 1 ) 2 a l .

3.2. Stress Solution

Let σ r and σ θ be the radial and circumferential stresses, respectively, referred to the ( r θ ) system introduced in (5). It has been shown above that the principal strain rate trajectories coincide with the coordinate curves of the Lagrangian coordinate system. The material model is coaxial. Therefore, the principal stress trajectories also coincide with the coordinate curves of the Lagrangian coordinate system. Then, Equation (5) implies that the radial and circumferential stresses are the principal stresses.

3.2.1. Case ζ 0 1

The stress solution is discontinuous across the curve ζ = ζ 0 . The yield criterion is
σ r σ θ = 2 σ 0 3
in the region 1 ζ ζ 0 and
σ r σ θ = 2 σ 0 3
in the region ζ 0 ζ 0 . Here, σ 0 is the yield stress in tension, a material constant. The only stress equilibrium equation that is not identically satisfied is
σ r r + σ r σ θ r = 0 .
Using (5), (9), (11), and (12), one transforms this equation into
σ r ζ = 4 a σ 0 3 [ 1 + 4 a ( ζ ζ 0 ) ]
in the region 1 ζ ζ 0 and
σ r ζ = 4 a σ 0 3 [ 1 + 4 a ( ζ ζ 0 ) ]
in the region ζ 0 ζ 0 .
Surface AB is traction free (Figure 2b). Therefore,
σ r = 0
for ζ = 0 . The tensile force is assumed to be balanced by a uniform pressure P applied over surface CD. For equilibrium, it is evident that P r C D = F . Using (6) and (9), one can rewrite this equation as
p = P σ 0 = 2 f a 1 4 a ( ζ 0 + 1 ) , f = F σ 0 H .
Then,
σ r σ 0 = p
for ζ = 1 . Moreover, the radial stress must be continuous across the curve ζ = ζ 0 :
[ σ r ] = 0
for ζ = ζ 0 . Here, the square brackets denote the amount of jump in an enclosed quantity.
The solution of Equations (14) and (15) satisfying conditions (16) and (18) is
σ r σ 0 = 1 3 ln [ 1 + 4 a ( ζ ζ 0 ) 1 4 a ζ 0 ]
in the region ζ 0 ζ 0 and
σ r σ 0 = 1 3 ln [ 1 4 a ( 1 + ζ 0 ) 1 + 4 a ( ζ ζ 0 ) ] p
in the region 1 ζ ζ 0 . Substituting (20) and (21) into (19) and using (17) determines f as a function of a at a given value of ζ 0 :
f = 1 4 a ( ζ 0 + 1 ) 2 3 a ln { [ 1 4 a ( 1 + ζ 0 ) ] ( 1 4 a ζ 0 ) } .
The right-hand side of this equation reduces to the expression 0 / 0 as a 0 . Applying l’Hospital’s rule, one can find that
f = 2 ( 1 + 2 ζ 0 ) 3
at a = 0 .
Equations (11), (12), (20), and (21) combine to give
σ θ σ 0 = 1 3 ln [ 1 + 4 a ( ζ ζ 0 ) 1 4 a ζ 0 ] + 2 3
in the region ζ 0 ζ 0 and
σ θ σ 0 = 1 3 ln [ 1 4 a ( 1 + ζ 0 ) 1 + 4 a ( ζ ζ 0 ) ] p 2 3
in the region 1 ζ ζ 0 . The bending moment referred to the central fiber is determined as
M = r C D r A B ( σ θ F h ) r d r .
Using (5), (6), and (17), one can transform this equation to
m = 2 3 M σ 0 H 2 = 3 a 1 0 σ θ σ 0 d ζ 3 f ( s s a ) .
The last term is a known function of a due to (9) and (22). One can evaluate the integral in (27) using (24) and (25). As a result,
1 0 σ θ σ 0 d ζ = 3 ln { [ 1 4 a ( 1 + ζ 0 ) ] ( 1 4 a ζ 0 ) } 4 3 a [ 1 + 2 ζ 0 + 3 p ( 1 + ζ 0 ) ] 12 a .
Here, p is a function of a due to (17) and (22). Then, substituting (28) into (27) determines m as a function of a:
m = 4 a + 8 a ζ 0 + ( 1 4 a ζ 0 ) [ 1 4 a ( 1 + ζ 0 ) ] ln { ( 1 4 a ζ 0 ) [ 1 4 a ( 1 + ζ 0 ) ] } 4 a 2 .
The right-hand side of this equation reduces to the expression 0 / 0 as a 0 . Applying l’Hospital’s rule, one can find that
m = 4 ζ 0 ( 1 + ζ 0 )
at a = 0 .

3.2.2. Case ζ 0 < 1

In this case, Equation (20) is valid in the range 0 ζ 0 1 . Equation (21) is immaterial. It follows from (18) and (20) that
p = 1 3 ln [ 1 4 a ζ 0 1 4 a ( 1 + ζ 0 ) ] .
This equation and (17) combine to give
f = 1 4 a ( ζ 0 + 1 ) 2 3 a ln [ 1 4 a ζ 0 1 4 a ( 1 + ζ 0 ) ] .
Applying l’Hospital’s rule, one can find that f = 2 / 3 at a = 0 .
Equation (27) is valid, and Equation (24) is valid in the range 0 ζ 0 1 . Then,
1 0 σ θ σ 0 d ζ = 1 3 [ 1 4 a ( 1 + ζ 0 ) ] 4 3 a ln [ 1 4 a ( 1 + ζ 0 ) 1 4 a ζ 0 ] .
Substituting the right-hand side of this equation into (27) yields
m = 1 a + ( 1 4 a ζ 0 ) [ 1 4 a ( 1 + ζ 0 ) ] 4 a 2 ln [ 1 4 a ( 1 + ζ 0 ) 1 4 a ζ 0 ] .
Applying l’Hospital’s rule, one can find that f = 0 at a = 0 .

4. Design Procedure at Loading

At a given value of ζ 0 , the solution in Section 3 determines the path in the (f, m)−space in parametric form, with a being the parameter. This path ensures that the process is an ideal flow process. A restriction imposed on the solution is that p 0 . It is seen from (17) that this inequality is equivalent to
f 0
if a 0 . Equation (32) shows that the inequality in (34) is always satisfied if ζ 0 < 1 .
Consider the case 0 ζ 0 1 . It follows from (22) and (31) that
1 4 a ( ζ 0 + 1 ) 2 3 a ln { [ 1 4 a ( 1 + ζ 0 ) ] ( 1 4 a ζ 0 ) } 0 .
Consider the equation
1 4 a ( ζ 0 + 1 ) 2 3 a ln { [ 1 4 a ( 1 + ζ 0 ) ] ( 1 4 a ζ 0 ) } = 0 .
It is seen from (10) that 1 4 a ( ζ 0 + 1 ) > 0 . Then, solving Equation (36) for a and using (35) leads to
a a c ( ζ 0 ) = 1 + 2 ζ 0 4 ζ 0 ( ζ 0 + 1 ) .
Let ρ C D and β be desirable values of ( r C D / H ) 2 and h / H at the end of loading. The inequality in (37) restricts these values. It follows from Equation (10) that
ρ C D = 1 4 a ( ζ 0 + 1 ) 4 a 2 and β = a 1 ( 1 4 ζ 0 a 1 4 a ( ζ 0 + 1 ) ) .
Solving these equations for ζ 0 and a leads to
a l = 1 β ( β + 2 ρ C D ) and ζ 0 = β ( β + 2 ρ C D ) 4 ρ C D β ( β + 2 ρ C D ) 1 .
Substituting (37) into (38) yields
ρ C D ( ζ 0 ) = 1 4 a c ( ζ 0 ) ( ζ 0 + 1 ) 4 a c 2 ( ζ 0 ) , β ( ζ 0 ) = 1 2 a c ( ζ 0 ) ( 1 4 ζ 0 a c ( ζ 0 ) 1 4 a c ( ζ 0 ) ( ζ 0 + 1 ) ) .
These equations determine a curve in the ( ρ C D ,   β ) space. This curve shown in Figure 3 corresponds to the condition f = 0 . The domains of the validity and invalidity of the solution are also indicated in this figure. Choosing ρ C D and β in the domain of validity, one can find ζ 0 and a l from (39). Equations (22) and (29) supply the loading path for this design if 0 ζ 0 1 , and Equations (32) and (33) if ζ 0 < 1 .

5. Unloading and Final Design

Unloading is assumed to be purely elastic. This assumption is verified a posteriori. Let ρ A B be the value of ( r A B / H ) 2 at the end of loading. The general solution for unloading is provided in [18]. Using this solution, one can find
R 0 = ρ C D { 1 + C 1 [ 1 ( 1 ρ 0 2 ) 2 ln ρ 0 ] 3 k p 2 ( 1 + ρ 0 2 2 ln ρ 0 ) } , R 1 = ρ A B { 1 + C 1 [ ρ 0 2 ( 1 ρ 0 2 ) 2 ln ρ 0 ] 3 k p ρ 0 2 2 ( 1 + 1 2 ln ρ 0 ) } .
Here, R 0 is the value of r C D / H after springback, R 1 is the value of r A B / H after springback, ρ 0 = ρ C D / ρ A B , k = σ 0 / ( 3 G ) and G is the shear modulus of elasticity. Additionally,
C 1 = B 1 A 1 , B 1 = 3 k f ln ρ 0 β 3 k p ρ 0 4 2 + k ρ 0 2 2 [ 3 p + 2 ln ρ 0 ( 3 m R 0 2 3 f β 3 p + 6 p ln ρ 0 ) ] , A 1 = ( 1 ρ 0 2 ) 2 4 ρ 0 2 ln 2 ρ 0 .
The quantities m, p, and f involved in (41) and (42) are understood to be calculated at the end of loading. The dimensionless strip’s thickness after springback is determined as
β = R 1 R 0 .
The solution in Section 3 supplies m, p and f as functions of ζ 0 and a l . Equation (10) provides ρ C D and ρ A B as functions of the same parameters. Since R 0 and β are prescribed as the design objectives, Equations (41) and (43) constitute a system for determining ζ 0 and a l . This system should be solved numerically.

6. Numerical Example

The numerical solutions were found for k = 0.0035 . This value of k is typical for ultra-high-strength steels [19]. The two design solutions shown in Table 1 are compared. Design (a) requires that the radius of CD (Figure 1) and strip thickness take prescribed values after springback. Design (b) requires that these quantities take prescribed values at the end of loading. For illustrative purposes, it is assumed that R 0 = 1 in Design (a) and ρ C D = 1 in Design (b).
Design (b) does not require numerical efforts other than solving the first equation in (40), from which the maximum value of ζ 0 is determined. This maximum value is ζ 0 0.661 . The corresponding value of β is determined from the second equation in (40) as β = β m 0.947 . The loading paths that require ζ 0 > 0.661 violate the inequality in (34). Several paths calculated using the procedure in Section 4 are depicted in Figure 4, Figure 5 and Figure 6. The curves in these figures are for β = β m , β = 0.9 , β = 0.8 , and β = 0.7 . The corresponding values of ζ 0 for the last three cases are ζ 0 = 0.73 ,   ζ 0 = 0.886 , and ζ 0 = 1.057 . Figure 4 and Figure 5 show the dependencies of f and m on a, respectively. Due to (5), these curves can be regarded as dependencies on θ 0 at any given value of η . Figure 6 depicts the loading paths in the (f, m)−space.
Design (a) requires a numerical method for solving Equations (41) and (43). This solution has shown that both designs lead to practically the same dependence of f on a. Therefore, this dependence is not illustrated for Design (a). The variation of m with a slightly depends on the design objectives. Figure 7 compares the two designs ( β = β m and β = β m in Figure 7a, β = 0.9 and β = 0.9 in Figure 7b, β = 0.8 and β = 0.8 in Figure 7c, and β = 0.7 and β = 0.7 in Figure 7d).
The solid curves correspond to Design (a) and the broken curves to Design (b). The difference decreases as the desired final thickness decreases and becomes invisible in Figure 7d.

7. Conclusions

A design solution for plane-strain bending was presented. The solution is based on the ideal flow theory and determines the loading path that ensures desirable geometric parameters after unloading. The solution is practically analytical. A numerical method is only necessary to solve a system of transcendental equations. Section 6 provides several loading paths that illustrate the general solution.
The two design objectives specified in Table 1 were compared. One of these designs requires that the process produces desirable geometric parameters of the sheet after springback. The other requires that these geometric parameters be obtained at the end of loading. Subsequent unloading changes the sheet’s geometry. The illustrative example in Section 6 showed that the difference between these two designs is very small. However, this conclusion is not universal and may depend on the design objectives.
The general structure of the found solution shows that it may be extended to strain-hardening materials, even though the general ideal flow theory is not available for this class of material models. This extension should allow for the ideal flow design of real sheet metal forming processes driven by springback to be provided.

Author Contributions

Conceptualization, S.A. and Y.-M.H.; formal analysis, M.R.; 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.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This paper was supported by the RUDN University Strategic Academic Leadership Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wagoner, R.H.; Lim, H.; Lee, M.-G. Advanced issues in springback. Int. J. Plast. 2013, 45, 3–20. [Google Scholar] [CrossRef]
  2. Cao, F.; Zeng, Y.; Liu, B.; Wu, W. A novel method for reducing springback of hot stretch bending parts. ICAACE 2019 IOP Conf. Ser. J. Phys. Conf. Ser. 2019, 1213, 052035. [Google Scholar] [CrossRef]
  3. Badrish, A.; Morchhale, A.; Kotkunde, N.; Singh, S.K. Parameter optimization in the thermo-mechanical V-bending process to minimize springback of inconel 625 alloy. Arab. J. Sci. Eng. 2020, 45, 5295–5309. [Google Scholar] [CrossRef]
  4. Wei, B.; Weib, Y.; Zhang, F.; He, K.; Dang, X.; Du, R. Influence of different heating methods on springback of mild steel plate during dieless bending process. Proc. Manuf. 2020, 50, 318–323. [Google Scholar] [CrossRef]
  5. Wei, B.; Wei, Y.; Zhang, F.; He, K.; Dang, X.; Du, R. Springback control and plastic deformation of metal plates with large curvature in heat-assisted incremental bending process. Int. J. Adv. Manuf. Technol. 2021, 112, 1483–1500. [Google Scholar] [CrossRef]
  6. Wasif, M.; Fatima, A.; Ahmed, A.; Iqbal, S.A. Investigation and optimization of parameters for the reduced springback in JSC-590 sheet metals occurred during the V-bending process. Trans. Indian. Inst. Met. 2021, 74, 2751–2760. [Google Scholar] [CrossRef]
  7. Chung, K.; Alexandrov, S. Ideal flow in plasticity. Appl. Mech. Rev. 2007, 60, 316–335. [Google Scholar] [CrossRef]
  8. Richmond, O.; Devenpeck, M.L. A die profile for maximum efficiency in strip drawing. Proc. 4th U.S. Natl. Congr. Appl. Mech. ASME 1962, 2, 1053–1057. [Google Scholar]
  9. Richmond, O.; Alexandrov, S. The theory of general and ideal plastic deformations of Tresca solids. Acta Mech. 2002, 158, 33–42. [Google Scholar] [CrossRef]
  10. Hill, R. A remark on diagonal streaming in plane plastic strain. J. Mech. Phys. Solids 1966, 14, 245–248. [Google Scholar] [CrossRef]
  11. Richmond, O.; Morrison, H.L. Streamlined wire drawing dies of minimum length. J. Mech. Phys. Solids 1967, 15, 195–203. [Google Scholar] [CrossRef]
  12. Richmond, O. Theory of streamlined dies for drawing and extrusion. Mech. Solid State 1968, 154–167. [Google Scholar] [CrossRef]
  13. Alexandrov, S.; Kim, J.-H.; Chung, K.; Kang, T.-J. An alternative approach to analysis of plane-strain pure bending at large strains. J. Strain Anal. Eng. Des. 2006, 41, 397–410. [Google Scholar] [CrossRef]
  14. Alexandrov, S.; Lyamina, E.; Hwang, Y.-M. Plastic bending at large strain: A review. Processes 2021, 9, 406. [Google Scholar] [CrossRef]
  15. Lee, Y.; Dawson, P.R. Obtaining residual stresses in metal forming after neglecting elasticity on loading. ASME. J. Appl. Mech. 1989, 56, 318–327. [Google Scholar] [CrossRef]
  16. Alexandrov, S.; Lyamina, E.; Hwang, Y.-M. Finite pure plane strain bending of inhomogeneous anisotropic sheets. Symmetry 2021, 13, 145. [Google Scholar] [CrossRef]
  17. Zhai, R.; Zhao, Z.; Yang, J.; Ma, B.; Yu, G. Analytical prediction of stretch-bending springback based on the proportional kinematic hardening model. Symmetry 2021, 13, 2389. [Google Scholar] [CrossRef]
  18. Alexandrov, S.; Manabe, K.; Furushima, T. A general analytic solution for plane strain bending under tension for strain-hardening material at large strains. Arch. Appl. Mech. 2011, 81, 1935–1952. [Google Scholar] [CrossRef]
  19. Liu, X.; Cao, J.; Huang, S.; Yan, B.; Li, Y.; Zhao, R. Experimental and numerical prediction and comprehensive compensation of springback in cold roll forming of UHSS. Int. J. Adv. Manuf. Technol. 2020, 111, 657–671. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the process: (a) initial configuration; and (b) intermediate and final configurations.
Figure 1. Schematic diagram of the process: (a) initial configuration; and (b) intermediate and final configurations.
Processes 10 01601 g001
Figure 2. Structure of the solutions: (a) general solution; and (b) ideal flow solution.
Figure 2. Structure of the solutions: (a) general solution; and (b) ideal flow solution.
Processes 10 01601 g002
Figure 3. Domains of validity and invalidity of the solution.
Figure 3. Domains of validity and invalidity of the solution.
Processes 10 01601 g003
Figure 4. Variation of the dimensionless force with a for several β values.
Figure 4. Variation of the dimensionless force with a for several β values.
Processes 10 01601 g004
Figure 5. Variation of the dimensionless bending moment with a for several β values.
Figure 5. Variation of the dimensionless bending moment with a for several β values.
Processes 10 01601 g005
Figure 6. Loading paths in the (f, m) space for several β values.
Figure 6. Loading paths in the (f, m) space for several β values.
Processes 10 01601 g006
Figure 7. Effect of the design objectives in Table 1 on the dependence of the dimensionless bending moment on a. The solid curves correspond to Design (a) and the broken curves to Design (b).
Figure 7. Effect of the design objectives in Table 1 on the dependence of the dimensionless bending moment on a. The solid curves correspond to Design (a) and the broken curves to Design (b).
Processes 10 01601 g007aProcesses 10 01601 g007b
Table 1. Design objectives.
Table 1. Design objectives.
Design (a) R 0 = 1 β = β m ,   or   β = 0.9 , or   β = 0.8 , or   β = 0.7
Design (b) ρ C D = 1 β = β m ,   or   β = 0.9 , or   β = 0.8 , or   β = 0.7
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.; Rynkovskaya, M.; Hwang, Y.-M. Ideal Flow Design of Plane-Strain Bending Driven by Springback. Processes 2022, 10, 1601. https://doi.org/10.3390/pr10081601

AMA Style

Alexandrov S, Rynkovskaya M, Hwang Y-M. Ideal Flow Design of Plane-Strain Bending Driven by Springback. Processes. 2022; 10(8):1601. https://doi.org/10.3390/pr10081601

Chicago/Turabian Style

Alexandrov, Sergei, Marina Rynkovskaya, and Yeong-Maw Hwang. 2022. "Ideal Flow Design of Plane-Strain Bending Driven by Springback" Processes 10, no. 8: 1601. https://doi.org/10.3390/pr10081601

APA Style

Alexandrov, S., Rynkovskaya, M., & Hwang, Y. -M. (2022). Ideal Flow Design of Plane-Strain Bending Driven by Springback. Processes, 10(8), 1601. https://doi.org/10.3390/pr10081601

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