Next Article in Journal
Complex Concentrated Alloys for Substitution of Critical Raw Materials in Applications for Extreme Conditions
Next Article in Special Issue
The Relation between Concrete, Mortar and Paste Scale Early Age Properties
Previous Article in Journal
Characterization of a New Lightened Gypsum-Based Material Reinforced with Fibers
Previous Article in Special Issue
Sizing and Topology Optimization of Trusses Using Genetic Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating of Bending Force and Curvature of the Bending Plate in a Three-Roller Bending System Using Finite Element Simulation and Analytical Modeling

1
Department of Mechanical Engineering, Dunarea de Jos University of Galati, Domneasca 47, 800008 Galati, Romania
2
Center of Excellence Polymer Processing, Dunarea de Jos University of Galati, Domneasca 47, 800008 Galati, Romania
*
Author to whom correspondence should be addressed.
Materials 2021, 14(5), 1204; https://doi.org/10.3390/ma14051204
Submission received: 3 February 2021 / Revised: 22 February 2021 / Accepted: 25 February 2021 / Published: 4 March 2021

Abstract

:
Many industries such as shipbuilding require steel bending plates in a wide range of radii, thus bending machines are often designed and produced on a custom basis in shipyards. From a design perspective, however, the bending force and the radius of the bending plate as a function of vertical displacement of the upper roller must be known. In this paper, a hybrid numerical–analytical approach is proposed to investigate the three-roller bending process for two plates of steel used in the naval industry. Firstly, the bending process is modeled using the finite element (FE) method and regression models for the bending force as a function of plate thickness and vertical displacement of the upper roller were constructed. Then, based on the findings from FE analysis, using the bent bar theory, two analytical expressions for the bending force were derived. Using geometric and deformation compatibilities, analytical expressions for the vertical displacement of the upper roller as a function of the curvature of the bending plate were also developed. The FE results suggest that the cross section of the plate is practically a plastic hinge in the tangent area of the upper roller and that the deformation compatibilities must be considered in order to estimate the curvature radius of the bending plate using analytical formulations. These results are of practical importance in designing rolling machines to estimate the setting parameters.

1. Introduction

Cylindrical and conical steel shells are important components used in various engineering applications such as bilge planking as a part of the naval hull, cylindrical tanks, pressure vessels, etc. [1]. Rolling machines with three or four rollers are generally considered for manufacturing shells (bending plates) with different curvatures [1,2,3,4]. In order to design a rolling machine, it is necessary to calculate the bending force and the curvature radius of the bending plate for given material properties and geometrical parameters. Therefore, efforts have been devoted to the prediction of the bending force or curvature radius by experimental, theoretical, and numerical approaches.
In general, the bending force is estimated using an analytical approach based on the one-dimensional beam theory combined with experimental analysis for validation. For example, Chudasama and Raval established an analytical model for the prediction of the bending force for the three-roller conical multi-pass bending process [4]. Furthermore, Chudasama and Raval proposed an analytical model for force prediction in the three-roller conical bending process, taking into account the effects of various material properties and geometrical parameters [5], which can also be used to estimate the roller plate-interface friction. Padgan et al. have conducted an analytical and experimental analysis of bending forces for five materials, such as aluminum alloy, copper alloy, stainless steel, gray cast iron, and magnesium alloy [6]. However, the analytical models are not able to fully describe the complex deformation behavior of the materials during the bending process due to oversimplified assumptions [2,7]. Therefore, in the attempt to optimize the roll bending process, in the last years, the finite element method (FEM) was applied to simulate the bending process and extract different information, such as bending force, the radius of curvature, stress-strain state in the contact area under various conditions, and the finite element (FE) environment [8,9,10,11,12,13]. For example, Tran investigated the forming process of an asymmetrical roll bending machine using a three-dimensional (3D) dynamic FE model in the Ansys/LS-DYNA environment and the results were validated by experiments [8]. Feng and Champliaud reported on the numerical simulation of cylindrical roll bending to predict the position of the lateral roll [9] and on the conical roll bending [10] using Ansys/LS-DYNA environment with explicit time integration. Ktari et al. carried out two-dimensional (2D) modeling of the rolling process in three roll bending pyramid systems using Abaqus explicit module [11]. Taylor et al. investigated the influence of various parameters of the three-roll bending process on the final radius of curvature of the bent sheet using an elastoplastic explicit dynamic FE method under the LS-DYNA environment [12]. Neto et al. developed an FE model to analyze the stress–strain in the vicinity of the contact area where the plastic deformation increases due to the forming tool [13]. Moreover, several attempts have been made to combine the FE simulation with the analytical modeling of different bending systems. Patel et al. proposed a moment curvature-based model for elastoplastic micro-beam bending to solve a micro-cantilever subjected to a normal load undergoing large deflection [14]. Pandit and Srinivasan presented an explicit numerical approach for three-point bending of a thin elastoplastic beam undergoing large deflection supported on cylindrical rollers with a radius comparable to the deflection [15]. Using both analytical and finite element approaches, Zemin et al. studied the three-roll bending of a large workpiece using a model based on one dimension beam theory to determine the inner radius displacement and the Abaqus FE model for optimizing the process parameters [16]. It was shown that, in general, the FE analysis provides more accurate results for the bending force and radius of the curvature compared with the experimental data than the oversimplified analytical models based on bending beam theory.
In summary, even though a number of researches have been carried out to improve the roll bending technology, estimation of the bending force and the curvature radius of the bending plates as a function of the displacement of the upper roller is still difficult at this stage, especially for large sheets and mainly relays on expensive physical experiments. Thus, the three-roller bending process needs further investigation both numerically and analytically. Moreover, the FE analysis, in combination with existing analytical theories, may offer a cheap and easy alternative to the trial-and-error, experimental approach.
The goal of the paper is to derive simple yet effective analytical expressions for calculating the bending force and the vertical displacement of the upper roller as a function of the radius of curvature, taking into account the findings from the FE analysis of the three-roller bending process. The novelty of this paper consists of coupling the plastic hinge condition observed during the FE analysis with the bent bar theory.

2. Numerical Simulation of the Three-Roller Banding Process

2.1. Research Methodology

Considering the complexity of the deformation during the rolling process, a hybrid numerical–analytical approach is proposed to investigate the three-roller bending process, as shown in Figure 1. First, the three-roller bending process is modeled using the FE method and Ansys Workbench Static Structural program (version 19.0, 2019, Ansys, Inc., Canonsburg, PA, USA), under plane strain conditions. In addition, regression models for the bending force as a function of plate thickness and vertical displacement of the upper roller are derived for two types of structural steel commonly used in the shipbuilding industry. Second, based on the findings from FE analysis, according to which the cross section of the plate is practically a plastic hinge in the tangent area of the upper roller, two analytical approaches for estimating the bending force in different hypotheses (frictional and frictionless) were derived based on the bent bar theory and compared with the FE results. Third, using geometric and deformation compatibilities, analytical expressions for the vertical displacement of the upper roller as a function of the radius of curvature, by taking into account the thickness of the sheet, were derived and verified against the FE results.

2.2. Finite Element Model

Figure 2 shows a schematic representation of the three-roller bending process with identical cylindrical rollers. This system employs one upper roller and two lower rollers as a forming tool. The radius of curvature of the bending plate is controlled by changing the vertical position of the upper roller. To reach a desire (final) curvature radius, the upper roller moves down in a vertical direction pressing down the plate.
Three design parameters, i.e., the vertical displacement of the upper roller, the force applied to the upper roller, and the torsional moments at the axes of the rollers, are required for design a three-roller bending system according to the imposed specifications, such as the plate thickness, final radius of the bending plate, and material properties.
To determine the design parameters for the three-roller bending system, in this paper, a static finite element analysis was carried out using the Ansys Workbench Static Structural program [17] with specific nonlinear settings. It should be pointed out that modeling large plastic deformations, which is the case of the three-rolling process, is based on the use of logarithmic deformations (Hencky), the Jaumann derivative of the Cauchy stress tensor, and the von Mises flow condition [17].
The bending system, particularly used in the naval shipyards, consists of three identical rollers with an outer diameter of 189 mm and a distance between the centers of the lower rollers of 320 mm. A previous study [18] indicated that, for the displacements of the upper roller equal to the roller radius, the applied force increases asymptotically and, therefore, the moments at the axes of rollers increase asymptotically, indicating that the technological solution may be inappropriate for the specified conditions.
To simulate the deformation behaviors, because the analysis of the entire three-roller bending process is very complicated, a two-dimensional (2D) model was considered. Because the plate’s width to thickness ratio is more than 200, and with the assumption of uniform loading along the width of the plate, the plane strain model is suitable. The simulation process was broken down into three stages as a function of time (as a pseudo-variable), which are (i) stage 1: The time varies in the range of 0–1 s in which the upper roller has an imposed vertical downward displacement; (ii) stage 2: The time varies in the range of 1–2 s in which the rollers have an imposed rotation; and (iii) stage 3: The time varies in the range of 2–3 s in which the upper roller has an imposed vertical upward displacement.
Figure 3a shows the 2D plane strain model for the three-roller bending process corresponding to simulation stages 1 and 2. In addition, a symmetrical 2D plane strain model corresponding to stage 1 (i.e., the vertical displacement of the upper roller) was developed to gain information that could be used in the analytical approach (on a symmetrical model) for the bending force calculation (Figure 3b). In this particular case, the center of the bottom roller was fixed, and the force is applied in the center of the upper roller along the axis of symmetry.
The models were drawn in the xy plane using the Design Modeler (version 19.0, 2019, Ansys, Inc., Canonsburg, PA, USA) included in the Static Structural module [17]. The origin of the reference system is in the center of the upper roller, whereas the x and y axes define the plane (in plane strain), where x is the horizontal axis and y is the vertical axis. The plate was discretized with 2D eight-node quadrilateral higher order elements (type PLANE 183 (Ansys, Inc., Canonsburg, PA, USA) [17], with large deflection and contact conditions included. These elements account for the nonlinearities introduced by the nonlinear behavior of the material, large displacements, and the use of contact elements between the roller surfaces and the plate [17]. The finite element size was 1 mm, resulting in 10 elements with that thickness for the 10-mm plate thickness in Figure 4. A total number of 15,000 elements were used to model the plate.
The contact between the plate and the rollers was considered frictional in the sense of Coulomb with a coefficient of friction of 0.3. In the 2D model, for the surface-to-surface contact, the “contact” and “target” elements are CONTA192 (Ansys, Inc., Canonsburg, PA, USA) and TARGE169 (Ansys, Inc., Canonsburg, PA, USA) for contact bodies and target bodies, respectively. The augmented Lagrangian formulation with asymmetric behavior was considered and the small sliding was considered off. The model worked without activating the interface treatment Adjust to Touch [17].
The simulation of the bending process was performed by increasing the vertical displacement of the upper roller from 0 to the maximum physical displacement imposed by the geometry of the system (see Table 1). However, to avoid numerical instability (i.e., the applied force increases asymptotically with increasing the radius of curvature [18]), the vertical displacement was limited to 50 mm. The simulations were carried out for two types of structural naval steels, for which an elastoplastic material model with isotropic hardening was considered, as in Table 1. It should be noted that the spring-back of the bending plate was not taken into account, and the rollers were assumed to be deformable in the linear-elastic domain.

2.3. Numerical Simulation Results

Figure 5 shows the variation of the normal stress σ X for the 10-mm plate thickness of S235JR steel. The results illustrate that the normal stress in each point of the cross section in the tangent area of the upper roller is equal to the yield strength.
It should be noted that the normal stress, σ X , is guided along the Ox-axis that is oriented along the neutral deformed surface of the plate, i.e., in each section, the Ox axis is “parallel” with the faces of the deformed plate. Furthermore, it was found that the distribution of stresses on a cross section depends on the mesh size when the calculation is in the elastoplastic domain.
The distribution of normal stress along the cross section in the tangent area of the upper roll is shown in Figure 6 for the S275JR steel. As can be seen, practically, the entire cross section of the plate is yielded, indicating that the deformation state corresponds to the “plastic hinge.”
Figure 7 shows the variation of the equivalent plastic strain in the cross section of the 10-mm plate in the tangent area of the upper roller. As can be seen, for almost the entire cross section, the plastic strain is greater than the σ Y /E ratio, which is 0.0013 for the S275JR steel.
Figure 8a shows the deformation of the 10-mm plate thickness at the end of stage 2. In order to estimate the final radius of the bending plate, the displacement of three points located on the inner surface of the plate was monitored during the FE analysis, as shown in Figure 8b. These points are approximately located on a circle with center ( x , y ) and radius r.
By using the least-square method, the center and the radius were determined by minimizing the sum
S ( x , y , r ) = i = 1 3 ( r r i ) 2 ,
in which r i , i = 1 ,   2 ,   3 , is the radius corresponding to each point in Figure 8b, i.e., the distance from the coordinates of each point ( x i , y i ) to the center ( x , y ) ,
r i = ( x x i ) 2 + ( y y i ) 2 .
Figure 9 shows the variation of the inner radius of the plate as a function of time for three vertical displacements. It can be seen that the inner radius decreases with increasing vertical displacement and tends to stabilize with increasing vertical displacement.
An approximate value for the inner radius related to a certain vertical displacement can be obtained by averaging the data over time. Therefore, the inner radius of the bending plate was found to be 218.6 ± 4.8 mm, 158.4 ± 2.9 mm, and 117.6 ± 1.5 mm for 30 mm, 40 mm, and 50 mm vertical displacement, respectively.

3. Analytical Modeling of the Three-Roller Bending Process

3.1. Radius of the Bending Plate

The radius of curvature of the bending plate (i.e., the radius of the neutral layer) can be determined based on the geometric compatibility of the three-roller bending system, shown in Figure 10, as follows:
R = R i + t 2 ,
where R i is the inner radius of the bending plate, and t is the plate thickness.
By invoking the geometric compatibilities of the bending system, taking into account that s i n ϕ 2 + c o s ϕ 2 = 1 , the following relation can be derived:
( L 2 R i + t + d 2 ) 2 + ( R i + t + y R i + t + d 2 ) 2 = 1 ,
where y is the distance between the horizontal axis of the lower rollers to the deformed outer layer of the plate and d is the diameter of the rollers.
After some mathematical manipulation and rearranging, the distance y reduces to
y = 1 ( L 2 R i + t + d 2 ) 2 ( R i + t + d 2 ) R i t .
Based on Figure 10, the vertical displacement, w, of the upper roller can be written as
w = d 2 y .
Substituting Equation (5) in Equation (6), the final expression for the vertical displacement as a function of the inner radius of the bending plate can be written as
w = d 2 1 ( L 2 R i + t + d 2 ) 2 ( R i + t + d 2 ) + R i + t .
The analytical model as in Equation (7), obtained only from geometrical compatibility and without taking into account the deformation mode of the plate or the spring-back, overestimates the inner radius of the bending plate as compared with the values predicted by the FE analysis. For the vertical displacement of 30 mm, 40 mm, and 50 mm, the relative error between the inner radius calculated by Equation (7) and FE values is 34%, 33%, and 33%, respectively. It should be pointed out that the greater the radius of the deformation, the greater the effect of the spring back. Thus, taking into account the deformation behavior from the FE analysis, which is similar to the deformation behavior of the plate during the bending process presented in [19], a new approach was proposed for calculating the inner radius as a function of the vertical displacement, as illustrated in Figure 11. It should be noted that the parameters marked by “*” stand for the correction according to the finding from the FE analysis.
Based on the new framework described in Figure 11, taking into account that ( s i n ϕ 1 * ) 2 + ( c o s ϕ 1 * ) 2 = 1 , the following relation can be derived:
[ L 2 ( R i * d 2 ) s i n ϕ * R i * + t + d 2 ] 2 + [ y * d 2 + ( R i * d 2 ) c o s ϕ * R i * + t + d 2 ] 2 = 1 ,
where
c o s ϕ * = d 2 d 2 + t ,   s i n ϕ * = 1 ( d 2 d 2 + t ) 2 .
After some mathematical manipulation and rearranging, the expression for y* takes the form
y * = 1 ( L 2 ( R i * d 2 ) s i n ϕ * R i * + t + d 2 ) 2 ( R i + t + d 2 ) d 2 ( R i * d 2 ) c o s ϕ * .
As before, the vertical displacement is
w * = d 2 y * .
Substituting Equation (10) into Equation (11), in term of the inner radius, following some manipulation of the resulting expression, the vertical displacement can be written as
w * = d 1 ( L 2 ( R i * d 2 ) s i n ϕ * R i * + t + d 2 ) 2 ( R i * + t + d 2 ) + ( R i * d 2 ) c o s ϕ * .
Equation (12) allows for a more accurate prediction of the inner radius of the bending plate as a function of vertical displacement.
Figure 12 compares the vertical displacement as a function of the inner radius of the bending plate corresponding to the geometric compatibility (as in Equation (7)), the deformation compatibility (as in Equation (12)), and the FE results, for the 10-mm plate thickness. As can be seen, the geometric compatibility model overestimates the vertical displacement, as compared with the deformation compatibility model, especially at the higher inner radius in which the effect of the spring-back is very important. The relative errors between the geometric model and deformation model are 45%, 35%, and 25% for 30 mm, 40 mm, and 50 mm vertical displacement, respectively. It can be seen that the higher the inner radius is, the higher the relative error. On the other hand, the vertical displacement of the upper roller can be predicted by Equation (12) relatively well. Compared with the FE results, the relative errors of the deformation compatibility model are 19%, 2%, and 12% for 30 mm, 40 mm, and 50 mm vertical displacement, respectively.

3.2. Bending Force

For estimating the bending force, in this study, the findings from the FE analysis were considered, namely, (i) the normal stress for each point of the cross section reaches the yield strength in the tangent area of the upper roller and (ii) the bending deformations of the plate beyond the tangent area are small. Therefore, in order to estimate the bending force, two scenarios were considered as given below.

3.2.1. Case 1 Scenario (C1)

The assumptions are as follows: (i) the plate has a finite thickness and width; (ii) the friction plays an important role in the three-roller bending process, i.e., the coefficient of friction between the rollers and plate is taken as 0.3; and (iii) the bending moment M p in the cross section of the plate situated in the tangent area of the upper roller corresponds to the so-called plastic hinge condition, i.e., the normal stress for each point of the cross-section becomes equal to the yield strength. Therefore, the bending moment can be expressed by
M p = σ Y l   ×   t 2 4 ,
where σ Y is the yielding strength, l is the width of the plate, and t is the thickness of the plate.
The initial configuration of the three-roller bending system is given in Figure 13a, while the configuration at a given time is presented in Figure 13b. From the notations given in Figure 13b, it follows that
t 1 = z   ×   c o s ϕ 1 d 2 t ,
and
L 2 s i n ϕ 1 + z c o s ϕ 1 d t = 0 ,
taking into account that
t a n α 1 = z L 2 .
Based on Equations (15) and (16), following some manipulation, it can be shown that
ϕ 1 + α 1 = a r c s i n [ d + t L 2   ×   c o s α 1 ] .
Furthermore, with the geometry of Figure 13b, the function b is given by
b = L 2 c o s ϕ 1 z s i n ϕ 1 ,
By taking into account the forces represented in Figure 13b, the bending moment M p , which corresponds to reaching the “plastic hinge” condition in the cross section, can be written as
M p = σ Y l   ×   t 2 4 = N   ×   b ,
considering the balance of forces in the vertical direction (in the y-axis direction), which implies that
F 1 = N c o s ϕ 1 + μ N s i n ϕ 1 ,
and the bending force is simply
F = 2   ×   F 1 .

3.2.2. Case 2 Scenario (C2)

The assumptions are as follows: (i) for estimating the deformed position of the plate, the thickness of the plate is negligible; (ii) the effect of friction is also negligible; (iii) the bending moment (Figure 14) corresponds to the high condition of plastic bending (i.e., the normal stress for each point of the cross-section is equal to the yield strength) and is given by Equation (13); and (iv) the moment in Section 3 is zero. It should be noted that, in Figure 14, the dotted lines represent the undeformed plate, while the solid line denotes the final position of the upper roller.
Considering the equilibrium of the forces that act on the three-roller bending system in Figure 14,
N   ×   b = σ Y l   ×   t 2 4 ,
the reaction force N is given by
N = σ Y l   ×   t 2 4 1 b ,
where b is the arm of the reaction force, and σ Y is the yield strength.
In the present context, the bending force can be calculated by
F = 2   ×   N c o s α .
Assuming that the vertical displacement of the upper roller can be approximated by
L 2 2 d 2 α = w ,
then, the angle α (in radians) is
α = L 2 w d .
Considering the three-roller bending system in Figure 14, the geometrical consideration leads to
2 d 2 s i n α + b c o s α = L 2 ,
and the arm b of the reaction force is given by
b = 1 c o s α [ L 2 d s i n α ] .
Substituting Equation (28) into Equation (23), the reaction force N becomes
N = σ Y l   ×   t 2 4 c o s α [ L 2 d s i n α ] .
The expression for N is then substituted into Equation (24) to obtain the bending force in the following form:
F = 2 σ Y l   ×   t 2 4 c o s 2 α [ L 2     d s i n α ] .

4. Predictive Models for Bending Force

In order to determine a predictive model for the bending force, FE simulations were carried out for different geometrical and material properties, as shown in Table 1. The FE simulations were performed with the model described in Section 2 for vertical displacement of the upper roller varying from 5 mm to 78 mm. The FE results for the bending force as a function of vertical displacement of the upper roller and the sheet thickness are shown in Figure 15 for S235JR and S275JR steels.
The results obtained from the FE simulations were used to determine a predictive model for the bending force as a function of two variables, i.e., vertical displacement of the upper roller and plate thickness. The FE results were fitted using Matlab curve fitting toolbox (version R2018a, The MathWorks, Inc., Natick, MA, USA) [20,21].
The regression models for the bending force (in Newton) for the two plates of steel are expressed by the following equations:
F ( w , t ) = 25.96 0.06384 w + 0.0542   e x p ( 0.001454 + 0.5549 t ) + 0.5711   w   e x p ( 1.852 + 0.2448 t )
for the S235JR steel, and
F ( w , t ) = 33.56 22.08 w + 0.008891   e x p ( 0.002553 + 0.7136 t ) + 1.123 w   e x p ( 2.909 + 0.0128 t )
for the S275JR steel.
The accuracy of the regression models was analyzed through the coefficient of determination (R2). The R2 coefficient for the S235JR and S275JR is 1.0% and 0.999%, respectively, indicating a very good correlation between the FE results and regression models. In addition, the adequacy of the derived models was also analyzed using the root-mean-square error (RMSE). The RMSE values for the S235JR and S275JR are 4.06 × 10−5% and 4.786 × 10−1%, respectively. These results show that the regression models are reliable and can be used to calculate the bending force.
Figure 16 compares the bending force as a function of vertical displacement of the upper roller for the two analytical models and the FE regression model (as in Equation (31)) for a plate thickness of 10 mm and S235JR steel. It can be seen that the C1 analytical model (frictional, as in Equation (21)) over predicts the bending force as compared with the C2 analytical model (frictionless, as in Equation (30)) and FE model. The relative errors for bending forces between the C1 and C2 models vary between 19% and 36%, and the higher the vertical displacement, the greater the relative error.
Compared with the FE regression model, the average relative error, the maximum relative error, and the RMSE for the C1 analytical model are 27%, 54%, and 18.2%, respectively, while the average relative error, the maximum relative error, and RMSE for the C2 model are 13.5%, 25%, and 12%, respectively. It should be noted that the maximum relative error of 54% for the C1 model and 25% for the C2 model occurs at 5 mm vertical displacement. For the C1 model, it appears that the presence of friction in the analytical formulation has a noticeable effect on the deformation behavior; thus, it predicts higher values as compared with the C2 model. However, the bending force values calculated based on the C2 model are in better agreement with the FE values than those calculated by the C1 model. Therefore, it can be concluded that the plastic hinge condition without friction is adequate to estimate the bending force with reasonable accuracy.

5. Conclusions

In this study, a hybrid approach based on the finite element (FE) analysis and analytical modeling of the three-roller bending process was proposed in order to predict the bending force and the curvature radius of the bending plate. The simulation of the three-roller bending process was carried out using the Ansys Workbench Static Structural program, under plane strain conditions. Starting from the observations following the FE analysis, analytical approaches for estimating the bending force and vertical displacement of the upper roller were derived. Based on the FE analysis and the analytical modeling, the following main findings and conclusions can be drawn:
(i) A 2D FE model was established to analyze the three-roller bending process for the S235JR and S275JR steels used in the naval industry. For the system under consideration, the FE results suggest that the cross section of the plate is practically a plastic hinge in the tangent area of the upper roller;
(ii) Taking into consideration the assumption of reaching the “plastic hinge” condition in the cross section of the plate, two analytical models for the bending force were derived based on the bent bar theory. These models can be used to first-hand estimate the bending force independently of FE analysis;
(iii) Using geometric and deformation compatibilities, simple analytical models for the vertical displacement of the upper roller as a function of the curvature radius of the plate were developed and verified against the FE results. However, in order to estimate the curvature radius of the bending plate using the analytical formulation, the deformation compatibilities must be considered;
(iv) FE-based regression models for estimation of the bending force for the S235JR and S275JR steels were formulated. The models take into consideration the plate thickness (8–12 mm) and vertical displacement of the upper roller (up to 78 mm). For other values of the yielding stress, the bending force can be obtained by interpolation (because the three-roller bending system has a diameter of 189 mm and the distance between the axis of the lower rollers of the system is 320 mm).
These results are of practical importance for the industry to estimate the setting parameters required in designing a three-roller bending machine. However, future research will address the experimental validation of the derived models.

Author Contributions

Conceptualization, D.B. and F.S.; methodology, I.G., D.B., and F.S.; software, I.G. and D.B.; validation, D.B. and I.G.; formal analysis, D.B. and F.S.; writing—original draft preparation, D.B. and F.S.; writing—review and editing, F.S.; funding acquisition, F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank Nicoleta-Violeta Stanciu for her support with the drawings.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roggendorff, S.; Haeusler, J. Plate bending: Three rolls and four rolls compared. Weld. Met. Fabr. 1979, 47, 353–357. [Google Scholar]
  2. Hua, M.; Lin, Y.H. Large deflection analysis of elastoplastic plate in steady continuous four-roll bending process. Int. J. Mech. Sci. 1999, 41, 1461–1483. [Google Scholar] [CrossRef]
  3. Shin, J.G.; Lee, J.H.; Kim, Y.I.; Yim, H. Mechanics-based determination of the center roller displacement in three-roll bending for smoothly curved rectangular plates. J. Mech. Sci. Technol. 2001, 1, 1655–1663. [Google Scholar] [CrossRef]
  4. Chudasama, M.K.; Raval, H.K. Analytical model for prediction of force during 3 Roller multipass conical bending and its experimental verification. Int. J. Mech. Eng. Robot. Res. 2012, 1, 91–105. [Google Scholar]
  5. Chudasama, M.K.; Rava, H.K. An approximate bending force prediction for 3-roller conical bending process. Int J. Mater. Form. 2013, 6, 303–314. [Google Scholar] [CrossRef]
  6. Padgan, N.P.; Deshpande, P.D.; Sakhale, C.N. Force Analysis of Metal Sheet in Bending Operation on Sheet Bending Machine. Int J. Eng. Res. Tech. 2015, 4, 267–270. [Google Scholar]
  7. Jadhav, C.S.; Talmale, P.S. Verification and Analysis of Stress-Strain Curve of Sheet Metal by using Three Roller Bending Machine. Int. J. Eng. Dev. Res. 2017, 5, 368–374. [Google Scholar]
  8. Tran, H.Q. Asymmetrical Roll Bending Process Study: Dynamic Finite Element Modeling and Experiments. Ph.D. Thesis, Université du Québec, Québec, QC, Canada, November 2014. [Google Scholar]
  9. Feng, Z.; Champliaud, H. Modeling and simulation of asymmetrical three-roll bending process. Simul. Model. Pract. Theory 2011, 19, 1913–1917. [Google Scholar] [CrossRef]
  10. Feng, Z.; Champliaud, H. Three stage process for improving roll bending quality. Simul. Model. Pract. Theory 2011, 19, 887–898. [Google Scholar] [CrossRef]
  11. Ktari, A.; Antar, Z.; Haddar, N.; Elleuch, K. Modeling and computation of the three-roller bending process of steel sheets. J. Mech. Sci. Technol. 2012, 26, 123–128. [Google Scholar] [CrossRef]
  12. Taylor, V.K.; Gandhim, A.H.; Moliya, R.D.; Raval, H.K. Finite element analysis of deformed geometry in three-roller plate bending process. In Proceedings of the ASME 2008 International Manufacturing Science and Engineering Conference, Illinois, IL, USA, 7–10 October 2008; Evanston: Illinois, IL, USA, 2008; pp. 1–7. [Google Scholar]
  13. Neto, D.M.; Martins, J.M.P.; Oliveira, M.C.; Menezes, L.F.; Alves, J.L. Evaluation of strain and stress states in the single point incremental forming process. Int J. Adv. Manuf. Tech. 2016, 85, 521–534. [Google Scholar] [CrossRef] [Green Version]
  14. Patel, B.N.; Pandit, D.; Srinivasan, S.M. Large elaso-plastic deflection of micro-beams using strain gradient plasticity theory. Procedia Eng. 2017, 173, 1064–1070. [Google Scholar] [CrossRef]
  15. Pandit, D.; Srinivasan, S.M. Large Elasto-Plastic Deflection of Thin Beams with Roller Support Contact, 11th International Symposium on Plasticity and Impact Mechanics. Procedia Eng. 2017, 173, 1079–1084. [Google Scholar] [CrossRef]
  16. Fu, Z.; Tian, X.; Hu, B.; Yao, X. Analytical modeling and numerical simulation for three-roll bending forming of sheet metal. Int. J. Adv. Manuf. Technol. 2013, 69, 1639–1647. [Google Scholar] [CrossRef]
  17. Ansys Workbench 2018, Version 19; Software For Technical Computation; Ansys, Inc.: Canonsburg, PA, USA, 2018.
  18. Gavrilescu, I.; Boazu, D. Simulation of roll bending with three rollers pyramid system using FEM. Ann. Univ. Dunarea de Jos of Galati Fascicle VI 2017, 35, 21–28. [Google Scholar]
  19. Hardt, D.E.; Constantine, E.; Wright, A. A model of a sequential bending process for manufacturing simulation. J. Eng. Ind. Trans. ASME 1992, 114/181, 181–187. [Google Scholar] [CrossRef]
  20. MATLAB 2018, Version R2018a; Software For Technical Computation; The MathWorks, Inc.: Natick, MA, USA, 2018.
  21. Kamiński, K.; Solecka, M. Optimization of the truss-type structures using the generalized perturbation-based Stochastic Finite Element Method. Finite Elem. Anal. Des. 2013, 63, 69–79. [Google Scholar] [CrossRef]
Figure 1. Flow chart for the hybrid numerical–analytical approach.
Figure 1. Flow chart for the hybrid numerical–analytical approach.
Materials 14 01204 g001
Figure 2. Schematic representation of the three-roller bending process (d—the outer diameter of the rollers, L—the distance between the lower rollers).
Figure 2. Schematic representation of the three-roller bending process (d—the outer diameter of the rollers, L—the distance between the lower rollers).
Materials 14 01204 g002
Figure 3. (a) A two-dimensional (2D) model for the three-roller bending process for stages 1 and 2 and (b) a 2D symmetric model for stage 1.
Figure 3. (a) A two-dimensional (2D) model for the three-roller bending process for stages 1 and 2 and (b) a 2D symmetric model for stage 1.
Materials 14 01204 g003aMaterials 14 01204 g003b
Figure 4. Discretization of the three-roller bending model.
Figure 4. Discretization of the three-roller bending model.
Materials 14 01204 g004
Figure 5. Distribution of the normal stress σ X (in Pa) for the 10-mm plate thickness (S235JR steel).
Figure 5. Distribution of the normal stress σ X (in Pa) for the 10-mm plate thickness (S235JR steel).
Materials 14 01204 g005
Figure 6. Distribution of the normal stress σ X (in Pa) in the tangent area of the upper roller (S275JR steel, 10-mm plate thickness).
Figure 6. Distribution of the normal stress σ X (in Pa) in the tangent area of the upper roller (S275JR steel, 10-mm plate thickness).
Materials 14 01204 g006
Figure 7. Equivalent plastic strain in the tangent area of the upper roller for the 1-mm plate thickness (S275JR steel).
Figure 7. Equivalent plastic strain in the tangent area of the upper roller for the 1-mm plate thickness (S275JR steel).
Materials 14 01204 g007
Figure 8. (a) Total deformation (in m) at the end of stage 2 and (b) location of the points for the estimation of the radius for the 10-mm plate thickness of S275JR steel.
Figure 8. (a) Total deformation (in m) at the end of stage 2 and (b) location of the points for the estimation of the radius for the 10-mm plate thickness of S275JR steel.
Materials 14 01204 g008
Figure 9. Variation of the inner radius of the bending plate as a function of time for different vertical displacements (10-mm plate thickness and S275JR steel).
Figure 9. Variation of the inner radius of the bending plate as a function of time for different vertical displacements (10-mm plate thickness and S275JR steel).
Materials 14 01204 g009
Figure 10. Schematic representation of the geometric compatibilities of the three-roller bending system.
Figure 10. Schematic representation of the geometric compatibilities of the three-roller bending system.
Materials 14 01204 g010
Figure 11. Schematic representation of the deformation compatibilities of the three-roller bending system based on FE observations.
Figure 11. Schematic representation of the deformation compatibilities of the three-roller bending system based on FE observations.
Materials 14 01204 g011
Figure 12. Vertical displacement of the upper roller as a function of the inner radius of the bending plate with 10 mm thickness (S275JR steel).
Figure 12. Vertical displacement of the upper roller as a function of the inner radius of the bending plate with 10 mm thickness (S275JR steel).
Materials 14 01204 g012
Figure 13. Configuration of the three-roller bending system for Case 1: (a) initial position of the rollers and undeformed plate and (b) position of the rollers and the plate at a given time and the diagram of forces applied on the plate.
Figure 13. Configuration of the three-roller bending system for Case 1: (a) initial position of the rollers and undeformed plate and (b) position of the rollers and the plate at a given time and the diagram of forces applied on the plate.
Materials 14 01204 g013
Figure 14. Schematic representation of the three-roller bending system and the forces applied to the sheet for Case 2.
Figure 14. Schematic representation of the three-roller bending system and the forces applied to the sheet for Case 2.
Materials 14 01204 g014
Figure 15. The bending force from FE analysis for (a) S235JR steel and (b) S275JR steel.
Figure 15. The bending force from FE analysis for (a) S235JR steel and (b) S275JR steel.
Materials 14 01204 g015
Figure 16. Bending force versus vertical displacement predicted based on the analytical models and the FE analysis for a 10-mm plate thickness and S235JR steel.
Figure 16. Bending force versus vertical displacement predicted based on the analytical models and the FE analysis for a 10-mm plate thickness and S235JR steel.
Materials 14 01204 g016
Table 1. Material properties and model parameters.
Table 1. Material properties and model parameters.
ParametersS235JRS275JR
Yield strength (MPa)235275
Young modulus (Pa)2.1 × 10112.1 × 1011
Tangent modulus (Pa)2.1 × 1092.1 × 109
Poisson ratio (-)0.30.3
Sheet thickness (mm)8, 10, 128, 10, 12
Maximum vertical displacement (mm)5050
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gavrilescu, I.; Boazu, D.; Stan, F. Estimating of Bending Force and Curvature of the Bending Plate in a Three-Roller Bending System Using Finite Element Simulation and Analytical Modeling. Materials 2021, 14, 1204. https://doi.org/10.3390/ma14051204

AMA Style

Gavrilescu I, Boazu D, Stan F. Estimating of Bending Force and Curvature of the Bending Plate in a Three-Roller Bending System Using Finite Element Simulation and Analytical Modeling. Materials. 2021; 14(5):1204. https://doi.org/10.3390/ma14051204

Chicago/Turabian Style

Gavrilescu, Ionel, Doina Boazu, and Felicia Stan. 2021. "Estimating of Bending Force and Curvature of the Bending Plate in a Three-Roller Bending System Using Finite Element Simulation and Analytical Modeling" Materials 14, no. 5: 1204. https://doi.org/10.3390/ma14051204

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