Next Article in Journal
Application of an Improved Link Prediction Algorithm Based on Complex Network in Industrial Structure Adjustment
Next Article in Special Issue
Customized Cost-Effective Cranioplasty for Large Asymmetrical Defects
Previous Article in Journal
Experimental Investigation on the Influence of Temperature on Coal and Gas Outbursts
Previous Article in Special Issue
Analyzing SPIF Product Characteristics Using Full Factorial Design-Integrated PCA Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mathematical Model for Force Prediction in Single Point Incremental Sheet Forming with Validation by Experiments and Simulation

1
Indian Institute of Technology (BHU), Varanasi 221005, Uttar Pradesh, India
2
Department of Mechanical Engineering, Vel Tech Rangarajan Sagunthala R & D Institute of Science and Technology, Avadi 600062, Tamil Nadu, India
3
Mechanical Engineering Department, Faculty of Engineering, Helwan University, Cairo 11732, Egypt
4
Department of Industrial Engineering, College of Engineering, King Saud University, P.O. Box 800, Riyadh 11421, Saudi Arabia
5
Industrial Engineering Department, College of Engineering, University of Houston, Houston, TX 77204, USA
*
Author to whom correspondence should be addressed.
Processes 2023, 11(6), 1688; https://doi.org/10.3390/pr11061688
Submission received: 26 March 2023 / Revised: 26 April 2023 / Accepted: 4 May 2023 / Published: 1 June 2023
(This article belongs to the Special Issue Computer-Aided Manufacturing Technologies in Mechanical Field)

Abstract

:
Incremental sheet forming (ISF) is an emerging technology that has shown great potential in forming customized three-dimensional (3D) parts without the use of product-specific dies. The forming force is reduced in ISF due to the localized nature of deformation and successive forming. Forming force plays an important role in modeling the process accurately, so it needs to be evaluated accurately. Some attempts have been made earlier to calculate the forming force; however, they are mostly limited to empirical formulae for evaluating the average forming force and its different components. The current work presents a mathematical model for force prediction during ISF in a 3D polar coordinate system. The model can be used to predict forces for axis-symmetric cones of different wall angles and also for incremental hole flanging. Axial force component, resultant force in the r-θ plane, and total force have been calculated using the developed mathematical model appearing at different forming depths. The cone with the same geometrical parameters and experimental conditions was modeled and simulated on ABAQUS, and finally, experiments were carried out using a six-axis industrial robot. The mathematical model can be used to calculate forces for any wall angle, but for comparison purposes, a 45° wall angle cone has been used for analytical, numerical, and experimental validation. The total force calculated from the mathematical model had a very high level of accuracy with the force measured experimentally, and the maximum error was 4.25%. The result obtained from the FEA model also had a good level of accuracy for calculating total force, and the maximum error was 4.89%.

Graphical Abstract

1. Introduction

Incremental Sheet Forming (ISF) is a flexible manufacturing process because complex three-dimensional (3D) parts can be fabricated by using only specific tool path programming on the existing setup without the use of dedicated dies. It is widely used for small batches and customized production. Patented in 1967, the process has attracted the eyes of aerospace, automotive [1], and biomedical sectors [2] for manufacturing complex sheet metal components. The process was used by Iseki et al. [3] for small batch production of non-symmetrical shallow shells using path-controlled spherical rollers. Various shapes were produced using the process of ISF, and a bulging test was used to predict the forming limit of ISF. In ISF, necking is either suppressed or delayed, due to which higher strains can be induced than in conventional stamping, deep drawing, and stretching processes [4].
Since the establishment of the process, there have been numerous attempts to establish the mechanism of ISF. Several FEA and theoretical analyses have been done to understand the underlying mechanism of the process. Initially, it was understood that the underlying mechanism behind the process is near biaxial stretching, as confirmed by Shim and Park [5] by FEA modeling. They concluded that Near Biaxial Stretching occurs in the closed loop, and near plane strain occurs along the straight line. They further concluded that the Forming Limit Curve (FLC), which is a plot of major strain vs. minor strain for ISF, is a straight line with a negative slope in the positive direction of minor strain and recommended the straight groove test as a standard test for plotting FLC for ISF [5]. Kim and Park predicted the same results by performing a series of experiments to produce various shapes [6].
Sheet thickness in the ISF process can be approximated by the famous sine law t = t 0 sin θ , where θ is the wall angle in the formed conical part. Since the thickness in the case of shear spinning is also given by the same sine law, the underlying deformation mechanism, as established by Avitzur and Yang [7] and demonstrated mathematically by Kim et al. [8], was initially considered to be pure shear through the thickness of the sheet and plane strain in the plane parallel to the undeformed sheet. Jackson et al. [9] found the presence of through-thickness shear along tool movement by measuring relative displacement in the sandwich panel. They took two pieces of 3 mm thick copper sheet, engraved their cross-section with square mesh, and brazed them. This method was adopted and suggested by Kalpakcioglu [10]. After deformation, they measured the mesh configuration and found the presence of through-thickness shear in the tool direction [11]. The presence of out-of-plane through-thickness stresses can lead to improved formability in the case of ISF, as suggested by Ma and Wang [12]. They further suggested that the conventional M-K model, which includes in-plane stresses only, needs to be extended and should also consider the through-thickness shear stresses. Allwood and Shoulder [13] proposed a generalized forming limit diagram (GFLD) that considered the out-of-plane stresses. It is evident that the stress state in ISF is complex, which was studied by Emmens and Boogard [14], and they presented an overview of the mechanisms responsible for higher strain in the single-point incremental forming (SPIF) process. They suggested that there are various factors responsible for higher strain in the ISF process. The nature of the complex stress state also plays a vital role in the higher strains in ISF [15]. Eyckens et al. [16] demonstrated, through finite element modeling, the presence of shear in the direction of tool movement but only a little shear in the perpendicular direction [8]. In addition to shear in the tool direction, Bending Under Tension (BUT) also plays a role in the enhancement of formability. Sawada et al. [17] performed a numerical study of stretch forming and ISF and found that bending and unbending happen along the meridional line, and shear happens along the circumferential direction. They also concluded that necking is induced on the part of the sheet which is in contact with the punch and stops after the punch surpasses the necking point. Recently, Li et al. [18] used the FE model to affirm that the deformation mechanism in ISF is a combination of bending, stretching, and shear.
Even after so many attempts, the underlying deformation mechanism of the process is still not clear due to the evolution of the complex stress state during the process. However, it can be said with affirmation that the stress state appearing during ISF consists of shear, normal, and bending stresses. Understanding of deformation mechanism is pivotal in developing an analytical model for force prediction in the process. Therefore, for developing such a model, all the stress components must be taken into consideration.
For in situ measurement of forces appearing during the ISF process, a dynamometer is fixed on the tool, and it measures the forces appearing during forming. The forces acting in ISF are (a) the axial force (Fz), acting in the direction along the tool axis, (b) force in the radial direction (Fr), and (c) force along the tangential direction (Fθ). The vector diagram of forces appearing during the process is given in Figure 1. Recently, Kumar et al. [19] presented a comprehensive review of the state of the arts of forming force and their effects on various forming outputs.
Duflou et al. [20] presented a model based on regression equations for force prediction during cone forming by SPIF as a function of step size, wall angle, tool diameter, and sheet thickness. Duflou et al. [21] performed a large set of experiments to form empirical relations for the calculation of force components for five different materials: AA3033, AA5754, DC01, AISI 304, and spring steel 65Cr2. Saidi et al. [22] conducted a set of experiments on sheets of AA1050 and 304 steel and validated their results, as predicted by Jeswiet et al. [23] and Aerens et al. [21]. Xiao et al. [24] measured forming forces by an in situ force monitoring system and studied the effect of forming temperature on the forming of an AA7075-T6 sheet. They found that the peak force was reduced to 1300 N from 1900 N when forming temperature was increased from 140 to 200 °C.
In addition to various experimental measuring methodologies, numerous attempts have been made to develop analytical and FEA models, considering various deformation mechanisms based on available literature for the calculation of forces and other parameters like strain, thickness, etc., during the process.
Flores et al. [25] used the FEA simulation of SPIF and studied the effects of constitutive laws on force prediction. They concluded that the selection of constitutive laws plays a vital role in force prediction, as the results were different for elastoplastic laws with isotropic and kinematic hardening. Henrard et al. [26] used FEA simulation to predict the forces for 20° and 60° cones using the Lagamine model, first-order brick element model for various simulation parameters, and inverse method for material data fit. They also conducted experiments for the validation of results obtained from FE modeling. Honarpisheh et al. [27] performed FEM analysis on ABAQUS to predict forces appearing during forming of Al/Cu bimetals and studied the effects of various parameters in ISF.
Iseki [28] used the equilibrium of the shell element to model bulge height, strain, and load during the process. He also used FEA analysis for validation of the results of the theoretical model. Silva et al. [29] proposed membrane analysis to consider the friction between the tool and sheet for analysis of the small plastic zone created during forming of symmetrical axis shapes in SPIF. They proposed that friction stress at the contact between tool and sheet can be assumed to be made of a couple of in-plane components, viz. meridional component ( μ Φ σ t ) and circumferential component ( μ θ σ t ). They considered the thickness, meridional, and circumferential direction as principal directions and gave the following equations for the corresponding stresses:
σ t = σ Y t ( r t o o l + t )
σ Φ = σ Y r t o o l ( r t o o l + t )
σ θ = σ Y ( r t o o l t ) 2 ( r t o o l + t )
Further, Chang et al. [30] used membrane analysis for the analytical modeling of the single-pass and multi-pass incremental sheet forming process (MPISF) and incremental hole flanging process (IHF). Bansal et al. [31] used the above equations for force calculation for single-stage and multi-stage ISF, modeling the contact area and validating with experiments and the formula given by Duflou et al. [20]. Similarly, Fang et al. [32] did analytical modeling of the process and validated it with numerical modeling, and also studied the fracture behavior of the sheet in ISF. Liu et al. [33] modeled the process in the 3D polar coordinate system and validated it experimentally. They further extended their work [34] and gave a comprehensive analytical model for the prediction of forming forces in three directions and studied the effect of step depth, tool diameter, and sheet thickness on forming forces occurring in ISF.
It can be concluded from the literature survey that although attempts have been made to develop a mathematical model of ISF, a comprehensive model for force prediction has not been established yet. Further, most of the models are empirical in nature, and they are used to calculate average forces during the process. Some other models consider only plane stresses in thickness, meridional, and circumferential direction and do not consider shear stresses. Very few works have been done considering all the six stress components to finally arrive at forces. In this current work, a simple mathematical model for the ISF process has been developed in a cylindrical coordinate system considering all the stress components. All force components have been calculated at different forming depths where the tool moves in different radii. The obtained results have been compared with the FEA model developed using the ABAQUS platform. Finally, experiments have been performed on aluminum alloy 6061 on a six-axis industrial robotic arm with the same experimental conditions as used for the analytical and FEA model. Different force components were measured using a dynamometer mounted on the robotic arm. The obtained results for various force components and total force have been compared.

2. Mathematical Model Development

For solving the mathematical equations, the cylindrical 3D coordinate has been chosen as tool motion to give a conical shape that can be easily replicated in this coordinate system. The center of the undeformed sheet has been chosen as the origin of the global coordinate system. The element has been taken in a local 3D polar co-ordinate system which is at a forming depth h and of length rdθ in the radial direction at an angle θ in the circumferential direction on the surface of the sheet in contact with the tool in the local polar coordinate system, as shown in Figure 2a,b. It is worth noting that the three directions—radial (r-direction), circumferential (θ-direction), and axial (z-direction)—may not necessarily be the principal directions. Hence all the six-stress component has been considered.
A few assumptions listed below have been made to develop the model:
  • All body forces have been neglected during the analysis;
  • τ θ z has been neglected throughout the analysis;
  • Normal stress in radial and z-directions are distributed in the same ratio as that of increments provided in this direction for forming a fixed wall angle cone. r z = tan α , σ r r σ z z = tan α , where α is the wall angle of the formed cone. The wall angle and step depths can be understood in Figure 3a;
  • For a given radius in which the tool is moving, τ r θ varies only along θ direction. If the considered element is under static equilibrium, the equilibrium equations in the differential form can be written directly [35]. The free-body diagram of the considered element is given in Figure 3b. Neglecting body forces, the equilibrium equations in the respective directions are:
In radial direction
σ r r r   +   1 r   τ r θ θ   +   σ r r σ z z r   +   τ r z z = 0
In circumferential direction
τ r θ r + 1 r σ θ θ θ + 2 r τ r θ = 0
In z direction
τ r z r + τ r z r + 1 r τ θ z θ + σ z z z = 0
The major shear component has been assumed to act only along the radial direction only. Hence τ θ z has been neglected throughout the model. Further, from membrane analysis, it can be said that
τ r z = f σ r r ,
Using the membrane analysis (4) and Equation (1) becomes
σ r r r   +   1 r   τ r θ θ   +   σ r r σ θ θ r   +   f σ r r z   =   0
Membrane analysis (4) and the assumption 2, Equation (3) becomes
f τ r r r + f r r r + σ z z z = 0
For any incremental sheet forming, the fabrication of 3D shapes is done by making the tool move over the sheet and providing increments to the tool so that it can perform fabrication in different passes. For fabricating a conical wall, the increments are provided in vertical and horizontal directions. In the current work, the increments are provided in radial and the vertically downward direction. These increments are generally small in comparison to r, and if the given assumption holds, then it can be assumed that, on neglecting spring back
z = r   tan α , where z is the increment in the vertical direction and r is the increment in radial direction. If these are small, then
r   tan α   =   z
Using assumption 3, it can be directly said that normal stress in radial and z-directions are distributed in the same ratio as that of increments.
σ r r tan α   =   σ z z
Taking differentials on both sides, it can be said that
σ r r tan θ   =   σ z z
From (7) and (8), it can be said that
σ r r r   =   σ z z z
Using Equation (9), Equation (6) further converges to
f σ r r r + f σ r r r + σ r r r = 0
σ r r r ( 1 + f ) = f σ r r r
Taking integrals on both sides of the above relation
σ r r σ r r = r r f 1 + f
l n σ r r = l n r f 1 + f + l n   C
σ r r = c r f 1 + f
From the above analysis, it is evident that for one complete cycle, the stress in the radial direction is constant. This analysis was performed on ABAQUS with the same parameters as used for the analytical model, and considering the element when the instantaneous radius of the undeformed part is 82.5 mm, the RMS value of the S11 was found to be approximately equal to 140 MPa as shown in Figure 3.
For f = 1 3 , C = 292.5
σ r r = 292.5 r 0.1667
Since σ r r depends on instantaneous circle radius only, hence
σ r r z = 0
From the equilibrium equation in the r-direction, we have
σ r r r + 1 r τ r θ z + σ r r σ θ θ r + σ r r z = 0
Using Equations (18) and (20)
r C r f 1 + f + 1 r τ r θ θ + C r ( r f 1 + f ) = σ θ θ r
We can directly say this using Equation (17)
σ r r 1 f f + 1 + τ r θ θ = σ θ θ
τ r θ θ   =   σ θ θ     σ r r ( 1 + f )
In the circumferential direction, using assumption 4, hence the equation in the circumferential direction converges to
1 r σ θ θ θ + 2 r τ r θ = 0
σ θ θ θ = 2 τ r θ
From Equations (23) and (25), we have
τ r θ σ θ θ = σ θ θ σ r r 1 + f 2 τ r θ
2 τ r θ τ r θ = σ θ θ σ θ θ + σ r r ( 1 + f ) σ θ θ
( τ r θ ) 2 = σ r r σ θ θ ( 1 + f ) σ θ θ 2 2 + k
If   σ θ θ = 0 ,   τ r θ = 0 ,   k = 0
τ r θ 2 = σ θ θ σ r r 1 + f σ θ θ 2
From Von-Mises criteria, we have
( σ r r σ θ θ ) 2 + ( σ r r σ z z ) 2 + ( σ θ θ σ z z ) 2 + 6 τ r θ 2 + τ r z 2 + τ θ z 2 = σ 2
The model is valid for all wall angles; however, for validation purposes, a 45° wall angle has been chosen because of its simplicity.
θ = 45 ° , σ r r = σ z z , τ θ z = 0
2 ( σ r r σ θ θ ) 2 + 6 τ r θ 2 + f 2 σ r r 2 = σ 2
σ r r 2 2 + 6 f 2 σ θ θ 2 + σ θ θ σ r r 2 4 f 1 + f = σ 2
If we consider the strain to be acting from bending and stretching only [34], we have from volume conservation
ε θ + ε r + ε f = 0
Total bending strain can be given by,   l n r t o o l R m
Whereas total tensile strain can be given as,   ε s t r e t c h i n g = l n t 0 t
Total strain can be given by
ε = l n r t o o l R m t 0 t
σ r r 2 2 + 6 f 2 σ θ θ 2 + σ θ θ σ r r 2 4 f 1 + f = σ 2
Taking the simple power law for annealed sheet, we have σ =   K ε n . From tensile test we have K = 205 MPa, and n = 0.2
σ r r 2 2 + 6 f 2 σ θ θ 2 + σ θ θ σ r r 2 4 f 1 + f = K 2 4 3 l n r t o o l t 0 R m t 2 n
From the analytical model developed, it can be said that the stress state appearing during the process is complex, but a few simple conclusions can be summarized:
(a)
The largest stress component is normal stress on the theta plane along the tool movement direction. Appreciable shear stress appears in the direction perpendicular to the tool motion along the radial and z planes;
(b)
From the Equation (28), it is clear that
σ r r 1 + f > σ θ θ 2

3. Contact Area Model

For contact area estimation, the approach of Bansal et al. [31] has been implemented. The contact area has been divided into two parts. The area directly below the tool in the region below the top of the undeformed sheet and the area above the undeformed sheet, as shown in Figure 4a,b. The area above the undeformed sheet leaves contact with the sheet tangentially in the circumferential direction. Let section XX′ be cut, and the top sectioned view is shown in Figure 4c.
From Figure 5b, it can be directly said that
cos φ = R Z R
The area of portion AB can be directly estimated by considering the concept of solid angle. The area can be given by
A A B = 2 π R Z
The area of the portion above the top of the undeformed surface can be given by integrating the perimeters of the strips of width d θ taken at an angle θ the length of which will be where r = R sin θ
A A B = P R d θ = R 2 β sin θ d θ
As θ   goes   from   φ   to   α , β   goes   from π 2 to 0.
Taking the variation of β to be linear with θ , we have ( β 0 ) = π 2 ( α φ ) ( θ α )
β = π 2 φ α θ α
From (39) and (40), we obtain
A 2 = R 2 π 2 φ α φ α ( θ α ) sin θ   d θ
A 2 = π R 2 2 φ α sin α + φ α cos φ sin φ
The total area of contact can be given by
A T = π R 2 2 φ α sin α + φ α cos φ sin φ + 2 π R Z
A T = π R 2 R + 3 z π R 2 2 α φ sin α sin φ
Putting R = 5 mm, z = 0.5 , α = π 4 and φ = 5 π 36 ,   A T = 19.05 .
Considering the whole contact area as a rectangle of l × b   area .
  l = R φ + α = 6.10 mm and b = 3.11 mm, the length of the portion above the top of the undeformed sheet is l C D = R α φ = 1.75 mm .
Total   projected   area   in   Z     plane   =   R sin α b   sin α = 7.775   mm 2
Total   projected   area   in   r plane = A r = Δ Z + R cos φ cos α b sin α = 3.28   mm 2
Total   projected   area   in   θ plane = A θ = t l = 4.52   mm 2
Different forces acting on the sheet in the contact region can be represented by the force diagram as given in Figure 1.
F r = σ r r × A r + z r θ A θ + z r z A z
F z = σ z z A z + z r z A r
F θ = σ θ θ A θ + z r θ A r
Once the relations for various force components were established, the calculation for forces using different developed equations was done on MATLAB for different instantaneous radii. The whole process was modeled on ABAQUS, and average forming forces were evaluated from the FEA model for the same instantaneous radii. The values obtained from the analytical model were compared with FEA and experimental results, as discussed in subsequent sections.

4. Simulation Study

The simulation analysis on the formed cone was carried out on ABAQUS. The finite element model was made on the software, which consisted of a rigid tool with a tool diameter of 10 mm, a deformable aluminum alloy 6061 sheet of thickness 1.05 mm, and a flange to hold the sheet. The geometry of the model used for simulation is given in Figure 5.
The properties of the material used in the simulation were obtained from the uniaxial tensile test and are given in Table 1.
For simulation, a true stress-true strain curve up to maximum load was obtained from the tensile test. Since a large strain is obtained in ISF for post-necking values of true stress and true strain, the following method was adopted. The log-log plot of true stress and the true strain was made till necking. Thereafter, a linear relation was assumed between the logarithm of true stress and the logarithm of true strain from necking to fracture. The strain hardening exponent in this region can be evaluated by finding the slope of this linear portion of the plot from necking to fracture. In this way, the true stress-true strain data was obtained and fed in ABAQUS till fracture. The engineering and true stress-strain curve obtained after tensile testing is given in Figure 6.
The isotropic hardening model was used for analysis as the sheet was annealed. Square meshing was done on the plate and tool with an average mesh size of 6.96 mm.
Loading and Boundary condition:
Penalty contact between the tool and sheet was taken with the coefficient of friction between the tool and the sheet to be 0.2, which was obtained from the wear test. The boundary of the sheet was encastred, and the blank was made rigid and encastred. The tip of the tool was chosen as a reference point and was given corresponding motions as per the conical geometry. The tool path was generated on MATLAB, and the corresponding positions were fed in ABAQUS to obtain the conical shape. The simulation was run, and dynamic analysis was done for stresses, strains, and tool forces arising during the tool motion, the results of which are given in Figure 7.
As can be seen from Figure 7a, the Von-Mises stress is maximum in the direct contact region between the undeformed sheet and the tool. The middle region of the cone is strained to the maximum value, as can be seen in Figure 7b. Further, it can also be seen from Figure 7c that the middle region of the cone has the minimum thickness showing that this region has undergone the maximum amount of thinning, which also affirms the distribution as seen in Figure 7b. A path along the central meridional plane was chosen in the cone after deformation and strain occurred along the path, and the thickness of the sheet along the path was plotted, which is shown in Figure 8a,b.
The simulation analysis was used to find out the forces appearing on the tool during the deformation. The tool force has been assumed to have two components. The axial components of the force Fzz and the resultant of in-plane X and Y components of forces is also equal to the resultant of tangential and radial force components.
As can be seen from Figure 9, the largest force component arising during the deformation is in the axial direction, as was reported by Duflou et al. [20] as well. Additionally, it was observed that the largest axial force appearing in the axial direction was found to be 2629.11 N. The mean force was found to be 1512.24 N. The variation of in-plane force, that is, the resultant of the forces appearing in the plane of the undeformed sheet with time, has been shown in Figure 9b, and the force diagram for the same has been given in Figure 9c. The largest value corresponding to in-plane forming force was found to be 2172.31 N, and the average force during forming was found to be 1109.3 N. Finally, the variation of total forming force with time has been given in Figure 9d. The maximum forming force appearing on the tool was found to be 2747.46 N, and the average tool force was found to be 1812 N.

5. Experimental Validation

Once the analytical model was developed and force calculations were done, a simulation was carried out for preliminary validation of obtained forces analytically. Finally, experiments were carried out on a six-axis industrial robotic arm. 1.05 mm thick aluminum 6061 alloy sheets were taken for experimentation. The material composition is given in Table 2.
Tensile tests and Erichsen cup tests were carried out for the determination of material properties before conducting the experiments. The sample was considered to be isotropic. The tensile test was carried out on the 100 kN INSTRON 8801 model. ASTM/E8/E8M standards were followed for the making of the specimens [36]. The sample after the test is shown in Figure 10.
Once the test was conducted, the results were calculated using the data obtained from the test, which are summarized in Table 3.
Before selecting a material, it is very crucial to have an idea of its formability. For the formability prediction of the sheet to be formed, the Erichsen cup test setup was used. The diameter of the used indenter was 20 mm. The main scale division was 1 mm, and the circular scale division was 50/5 MSD. The measured depth (HE) for all the domes was measured, which is given in Table 4. The domes are made as a result of indentation made by the indenter, as shown in Figure 11.
After obtaining the mechanical parameters from the uniaxial tensile test and Erichsen cup test, experiments were performed for which the chosen parameters are given in Table 5. The parameters were obtained using surface regression models. Dimension of the undeformed workpiece was taken as 250 mm × 250 mm × 1.05 mm. The sheet was deformed into a cone with a radius of the top circle 96 mm. The whole setup, the forming tool with the tool dynamometer, and the deformed axis-symmetric cone are shown in Figure 12a–c.

6. Results and Discussion

6.1. Thickness Distribution

The instantaneous thickness has a role to play in Equation (18). Hence the state of stress and the force are affected by the instantaneous sheet thickness, so the value of the sheet thickness was obtained from the FEA model and was validated by measuring the values after the experiments.
Once the experiments had been carried out, the thickness of the sheet was measured. The measurement was done by a micrometer with conical tips of least count 0.001 mm. The sheet was found to have the minimum thickness in the middle of the cone. The trend was also shown in previous work by Ambrogio et al. [37], who showed that the thickness prediction by the sine law is most accurate in the middle region of the deformed sheet. The conical region was divided into seven sub-regions (0–6), as shown in Figure 13a, and the thickness in the different regions was measured. Region-0 is the un-deformed region where thickness is 1.05 mm, and Region-6 is the region where the material gets accumulated because the tool moving along the sheet drags the material with it, and finally, it gets accumulated in the lowest region of the sheet; hence, the sheet thickness exceeds the original sheet thickness. However, Region-1–5 is the region in which sheet thinning occurs. The thickness plot for ISF in FEA and experiments are given in Figure 13c. For thickness prediction using FEA analysis, a path was chosen in the meridional direction, as shown in Figure 13b. The thickness variation obtained from FEA is shown in Figure 13c. All three thickness variations, obtained experimentally by FEA and by sine law [4], are plotted along the meridional path, as shown in Figure 13.

Force Measurement

For the measurement of forces, the cone was divided into seven zones. Zone 0 was taken as the zone of the upper undeformed cone, the middle region was taken as Zones 1–5 in which near uniform thinning occurs, and the sheet thinning followed the sine law, and Zone 6 was the region between the lower undeformed sheet and the cone, as shown in Figure 13a. The force was measured during deformation using a drill tool dynamometer. The setup is shown in Figure 12a,b. The capacity of the dynamometer is 500 kgf having a strain gauge based 350 Ω bridge sensor. The force was measured in X, Y, and Z directions, and the resultant of the forces in X and Y directions was calculated. Three forces were defined as Fz, Fp (resultant of forces in X and Y directions), and total force (Ft). Average values of Fz, Fp, and Ft were calculated for one complete cycle and were compared with the values obtained analytically. For comparison with the forces obtained from ABAQUS analysis, the average of forces for one complete circular cycle was calculated. Thus, all the components of forces were obtained from the analytical model, as well as the FEA model, using ABAQUS, and final validation was done by experimental results. The corresponding values of forces were plotted against the values of the instantaneous radius of the cone.

6.2. Variation of Axial Force

This force acts in the vertical direction. The force was calculated via an analytical model as well as the FEA model. Both were validated using final experimentation. The peak and the average value of the axial forces obtained in all three cases are given in Table 6.
The axial force increased with decreasing radius, as observed in the analytical model. The force in the axial direction is given by F z z = σ z z A z + z r z A r . The values of forces were calculated for different radii of the undeformed sheet using the formula, and the force was found to increase with the decreasing radius. The variation of the axial force with decreasing radii was almost linear. The same trend was seen in the FEA model, and experimentally, it was seen to increase with decreasing radii, but the shape of the plot departed from linearity slightly in the middle zone, as shown in Figure 14a,b. This can be due to sheet thinning, which can give rise to inaccurate sheet thickness at different points for a given radius. Nevertheless, the model gives a very good approximation of axial forming force during deformation. The error function for the error of the forces obtained analytically was defined as
Percentage   error = F analytical F experimental F experimental × 100
The maximum error was found to be 4.92%. The force trends, as observed in ABAQUS, also closely matched the data obtained experimentally. The maximum error obtained in this case was 5.52%. The error in the axial force with respect to the experimental force calculated at every value of radius is shown in Section 6.4.

6.3. Variation of Resultant Force in r-θ Plane

The force in the r-θ plane has two components. For the mathematical model, the two components of forces evaluated were the ones in the r and θ planes. However, the components of the forces measured by the dynamometer and the forces observed in the ABAQUS analysis are in X and Y directions. So, the resultant of Fr and Fθ was calculated, which was equal to the resultant of FX and FY. Finally, the force obtained analytically, and that obtained in FEA was compared with the experimental force. The peak and the average value of the resultant force in r-θ obtained in all three cases are given in Table 7.
The plot of forces obtained in the r-θ plane from the analytical model, FEA, and experiment is given in Figure 15. The force obtained analytically had a good level of accuracy with experimental force. The largest error was found to be 8.47%, and the force obtained on FEA was found to have a maximum deviation of 4.61% from the values measured experimentally. The error in the axial force with respect to the experimental force calculated at every value of radius is shown in Section 6.4.

6.4. Variation of Total Force

In the analytical model, the total force appearing on the forming tool-sheet interface was calculated by the formula F t = F p 2 + F z 2 . Once the force was calculated, it was compared with the total force obtained experimentally and the one obtained in the FEA model. The peak and the average value of the total force obtained in all three cases are given in Table 8.
The plot of total force obtained analytically, by the FEA model, and experimentally is given in Figure 16. The computed axial force from the mathematical model has a good level of accuracy as far as the total force is concerned. The maximum error in the calculated force from the total experimental force was found to be 4.25%.
The total force was also obtained from the FEA model on ABAQUS. The force obtained from the FEA model was also in close agreement with the total force obtained experimentally. The maximum error in obtained total force in the FEA model from the experimental model was found to be 4.89%, as shown in Figure 17c.

6.5. Error in Calculated and Experimentally Determined Forces

The forces calculated from the analytical model were compared with forces obtained by the FEA model and experimentally. For validation purposes, the true value of force was taken as the one which was obtained experimentally. The values of the forces were obtained from the analytical model for different instantaneous radii, and the corresponding values of the forces were also obtained from ABAQUS data. Both the values were compared with experimental values. So, validation of the mathematical model was done by comparing the values of the forces with the experimental values.
The   error   function   was   defined   as   percentage   error = F calculated F experimental F experimental × 100
For axial forces and the total forces, the error was maximized in Zone 2, i.e., the middle region of the cone. However, for the resultant force in the r-θ plane, the error was maximum in Region-1, which decreased, and then the error was increased in Region-2 before finally decreasing in Region-3.
As far as the FEA model is concerned, the validation was done by comparing the forces obtained experimentally. The error for all three forces followed the same trend. The percentage error vs. instantaneous radius followed a bell curve pattern. It was maximized in Region-2 and was minimal in Region-1 and Region-3. The error curves are given in Figure 17a–c.

7. Conclusions

This paper presents an analytical model to predict the forming forces during incremental sheet forming. The modeling was done in a polar 3D coordinate system. All six stress components were taken into consideration initially before neglecting τ θ z . Analytical formulae were developed for all the stress components. The contact area was modeled analytically before finally developing the equations for all forces. The three force components for which the equations were developed were Fr, Fθ, and Fz. The result of the three components is the total force Ft. The resultant of Fr and Fθ is the resultant force in the r-θ plane, represented as Fp. The material used in the mathematical model was AA6061. The process was modeled in FEA using the same experimental parameters, and the force components Fz and Fp were calculated before finally calculating the total force (Ft). Finally, the process was carried out on a six-axis industrial robot, and forces were measured using a tool dynamometer. The data obtained from the analytical model and FEA were validated by the experimental results. The model is for all wall angle cones; however, in this work, the validation has been done for a 45° cone. The results obtained from the analytical model and FEA were in good agreement with the experimental result. The maximum error in the analytically calculated total force was found to be 4.25%. Similarly, the results obtained from the FEA model were also in agreement with the experimental results. The maximum error in evaluated total force from the FEA model was found to be 4.89% with respect to the experimental results. It can be concluded that the mathematical model developed can be used for a force calculation to a good level of accuracy.

Author Contributions

Conceptualization, R.P.S. and M.M.; Methodology, R.P.S.; Software, R.P.S.; Validation, R.P.S., M.M., A.K.S. and E.S.A.N.; Formal analysis, R.P.S. and S.K.; Investigation, R.P.S. and M.M.; Resources, S.K.; Data curation, P.K.S.; Writing—original draft, R.P.S.; Writing—review & editing, S.K., P.K.S., S.S. and A.K.; Visualization, A.K.S.; Supervision, S.K.; Project administration, S.K.; Funding acquisition, S.S. and H.M.A.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [King Saud University] granted number [RSP2023R164] and [DST SERN NEWDELHI] granted number [MMER/2014/0068].

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to DST-SERB for sponsoring the work through project MMER/2014/0068 titled ‘Design, Development and Fabrication of an Incremental sheet hydroforming machine setup’. The authors extend their appreciation to King Saud University for funding this work through Researchers Supporting Project number (RSP2023R164), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

σ z z Normal stress in axial direction
σ r r Normal stress in radial direction
σ θ θ Normal   stress   in   circumferential   direction
τ r θ Shear   stress   in   radial   plane   in   circumferential   direction = τ θ r
τ r z Shear   stress   in   radial   plane   in   axial   direction = τ z r
τ z θ Shear   stress   in   axial   plane   in   circumferential   direction = τ θ z
f Coefficient   of   kinetic   friction
t o Initial   thickness   of   the   sheet
t Instantaneous   thickness   of   the   sheet
r Radius   of   current   pass   ( Instantaneous   radius )
R t Tool   radius
R m Radius   of   curvature   of   neutral   axis   of   the   taken   element
α Wall   angle   of   cone
F z Axial   force   component
F r Radial   force   component
F θ Tangential   force   component
F P Resultant   force   in   r θ   plane
σ t Normal stress in thickness direction
σ Φ Normal stress meridional direction
σ θ Normal stress circumferential direction
σ Y Normal stress perpendicular to the tool in the vertical direction

References

  1. Edward, L. Apparatus and Process for Incremental Dieless Forming. U.S. Patent 3,342,051, 19 September 1967. [Google Scholar]
  2. Jeswiet, J.; Micari, F.; Hirt, G.; Bramley, A.; Duflou, J.; Allwood, J. Asymmetric Single Point Incremental Forming of Sheet Metal. CIRP Ann. 2005, 54, 88–114. [Google Scholar] [CrossRef]
  3. Iseki, H. Flexible and incremental sheet metal bulging using a few spherical rollers. Trans. Jpn. Soc. Mech. Eng. 1993, 59, 2849. [Google Scholar] [CrossRef]
  4. Kim, T.; Yang, D.-Y. Improvement of formability for the incremental sheet metal forming process. Int. J. Mech. Sci. 2000, 42, 1271–1286. [Google Scholar] [CrossRef]
  5. Shim, M.-S.; Park, J.-J. The formability of aluminum sheet in incremental forming. J. Mater. Process. Technol. 2001, 113, 654–658. [Google Scholar] [CrossRef]
  6. Kim, Y.H.; Park, J.J. Effect of process parameters on formability in incremental forming of sheet metal. J. Mater. Process. Technol. 2002, 130, 42–46. [Google Scholar] [CrossRef]
  7. Avitzur, B.; Yang, C. Analysis of power spinning of cones. J. Eng. Ind. Aug. 1960, 82, 231–244. [Google Scholar] [CrossRef]
  8. Kim, J.H.; Park, J.H.; Kim, C. A Study on the mechanics of shear spinning of cones. J. Mech. Sci. Technol. 2006, 20, 806–818. [Google Scholar] [CrossRef]
  9. Jackson, K.; Allwood, J.; Landert, M. Incremental forming of sandwich panels. J. Mater. Process. Technol. 2008, 204, 290–303. [Google Scholar] [CrossRef]
  10. Kalpakcioglu, S. On the Mechanics of Shear Spinning. J. Eng. Ind. 1961, 83, 125–130. [Google Scholar] [CrossRef]
  11. Jackson, K.; Allwood, J. The mechanics of incremental sheet forming. J. Mater. Process. Technol. 2009, 209, 1158–1174. [Google Scholar] [CrossRef]
  12. Ma, L.; Wang, Z. The effects of through-thickness shear stress on the formability of sheet metal—A review. J. Manuf. Process. 2021, 71, 269–289. [Google Scholar] [CrossRef]
  13. Allwood, J.M.; Shouler, D.R. Generalised forming limit diagrams showing increased forming limits with non-planar stress states. Int. J. Plast. 2009, 25, 1207–1230. [Google Scholar] [CrossRef]
  14. Emmens, W.; Boogaard, A.V.D. An overview of stabilizing deformation mechanisms in incremental sheet forming. J. Mater. Process. Technol. 2009, 209, 3688–3695. [Google Scholar] [CrossRef]
  15. Emmens, W.; Boogaard, A.V.D. Incremental forming by continuous bending under tension—An experimental investigation. J. Mater. Process. Technol. 2009, 209, 5456–5463. [Google Scholar] [CrossRef]
  16. Eyckens, P.; Van Bael, A.; Aerens, R.; Duflou, J.; Van Houtte, P. Small-scale Finite Element Modelling of the Plastic Deformation Zone in the Incremental Forming Process. Int. J. Mater. Form. 2008, 1, 1159–1162. [Google Scholar] [CrossRef]
  17. Sawada, T.; Fukuhara, G.; Sakamoto, M. Deformation mechanism of sheet metal in stretch forming with computer numerical control machine tools. J. Jpn. Soc. Technol. Plast. 2001, 42, 1067–1069. [Google Scholar]
  18. Li, Y.; Daniel, W.J.T.; Meehan, P.A. Deformation analysis in single-point incremental forming through finite element simulation. Int. J. Adv. Manuf. Technol. 2017, 88, 255–267. [Google Scholar] [CrossRef]
  19. Kumar, A.; Gulati, V.; Kumar, P.; Singh, H. Forming force in incremental sheet forming: A comparative analysis of the state of the art. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 251. [Google Scholar] [CrossRef]
  20. Duflou, J.R.; Tunckol, Y.; Aerens, R. Force Analysis for Single Point Incremental Forming. Key Eng. Mater. 2007, 344, 543–550. [Google Scholar] [CrossRef]
  21. Aerens, R.; Eyckens, P.; Van Bael, A.; Duflou, J.R. Force prediction for single point incremental forming deduced from experimental and FEM observations. Int. J. Adv. Manuf. Technol. 2010, 46, 969–982. [Google Scholar] [CrossRef]
  22. Saidi, B.; Boulila, A.; Ayadi, M.; Nasri, R. Experimental force measurements in single point incremental sheet forming SPIF. Mech. Ind. 2015, 16, 410. [Google Scholar] [CrossRef]
  23. Jeswiet, J.; Duflou, J.R.; Szekeres, A. Forces in Single Point and Two Point Incremental Forming. Adv. Mater. Res. 2005, 6–8, 449–456. [Google Scholar] [CrossRef]
  24. Xiao, X.; Kim, C.-I.; Lv, X.-D.; Hwang, T.-S.; Kim, Y.-S. Formability and forming force in incremental sheet forming of AA7075-T6 at different temperatures. J. Mech. Sci. Technol. 2019, 33, 3795–3802. [Google Scholar] [CrossRef]
  25. Flores, P.; Duchêne, L.; Bouffioux, C.; Lelotte, T.; Henrard, C.; Pernin, N.; Van Bael, A.; He, S.; Duflou, J.; Habraken, A. Model identification and FE simulations: Effect of different yield loci and hardening laws in sheet forming. Int. J. Plast. 2007, 23, 420–449. [Google Scholar] [CrossRef]
  26. Henrard, C.; Bouffioux, C.; Eyckens, P.; Soly, H.; Duflou, J.; Van Houtte, P.; Van Bael, A.; Duchêne, L.; Habraken, A. Forming forces in single point incremental forming: Prediction by finite element simulations, validation and sensitivity. Comput. Mech. 2011, 47, 573–590. [Google Scholar] [CrossRef]
  27. Honarpisheh, M.; Keimasi, M.; Alinaghian, I. Numerical and experimental study on incremental forming process of Al/Cu bimetals: Influ-ence of process parameters on the forming force, dimensional accuracy and thickness variations. J. Mech. Mater. Struct. 2018, 13, 35–51. [Google Scholar] [CrossRef]
  28. Iseki, H. An approximate deformation analysis and FEM analysis for the incremental bulging of sheet metal using a spherical roller. J. Mater. Process. Technol. 2001, 111, 150–154. [Google Scholar] [CrossRef]
  29. Silva, M.; Skjoedt, M.; Martins, P.; Bay, N. Revisiting the fundamentals of single point incremental forming by means of membrane analysis. Int. J. Mach. Tools Manuf. 2008, 48, 73–83. [Google Scholar] [CrossRef]
  30. Chang, Z.; Li, M.; Chen, J. Analytical modeling and experimental validation of the forming force in several typical incremental sheet forming processes. Int. J. Mach. Tools Manuf. 2019, 140, 62–76. [Google Scholar] [CrossRef]
  31. Bansal, A.; Lingam, R.; Yadav, S.K.; Reddy, N.V. Prediction of forming forces in single point incremental forming. J. Manuf. Process. 2017, 28, 486–493. [Google Scholar] [CrossRef]
  32. Fang, Y.; Lu, B.; Chen, J.; Xu, D.; Ou, H. Analytical and experimental investigations on deformation mechanism and fracture behavior in single point incremental forming. J. Mater. Process. Technol. 2014, 214, 1503–1515. [Google Scholar] [CrossRef]
  33. Liu, F.; Li, Y.; Cheng, Z.; Wang, Z.; Li, F.; Li, J. Preliminary modelling of forming forces in three directions for incremental sheet forming process based on the contact area. Procedia Manuf. 2020, 50, 630–636. [Google Scholar] [CrossRef]
  34. Liu, F.; Li, Y.; Yuan, H.; Sun, R.; Li, F.; Li, J. Comprehensive modeling of forming forces in three directions for incremental sheet forming process based on the contact area. J. Manuf. Process. 2022, 84, 986–1000. [Google Scholar] [CrossRef]
  35. Dalir, M.A.; Seifi, R. Direct method for deriving equilibrium equations in solid continuous systems. Eng. Solid Mech. 2014, 2, 321–330. [Google Scholar] [CrossRef]
  36. E8/E8M-16a; Standard Test Methods for Tension Testing of Metallic Materials. ASTM International: West Conshohocken, PA, USA, 2016.
  37. Ambrogio, G.; Filice, L.; Gagliardi, F.; Micari, F. Sheet Thinning Prediction in Single Point Incremental Forming. Adv. Mater. Res. 2005, 6–8, 479–486. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic diagram of ISF (b) Forming forces in the plane of an undeformed sheet (c) Vector diagram of forming forces in the plane of an undeformed sheet.
Figure 1. (a) Schematic diagram of ISF (b) Forming forces in the plane of an undeformed sheet (c) Vector diagram of forming forces in the plane of an undeformed sheet.
Processes 11 01688 g001
Figure 2. (a) Definition of the 3D polar coordinate for considering an element, and (b) representation of the taken element. (c) Definition of wall angle in terms of step depth. (d) Stresses acting on the taken element.
Figure 2. (a) Definition of the 3D polar coordinate for considering an element, and (b) representation of the taken element. (c) Definition of wall angle in terms of step depth. (d) Stresses acting on the taken element.
Processes 11 01688 g002
Figure 3. The process simulated under the same condition on ABAQUS: (a) the radial stress distribution at a radius of 82.5 mm and (b) the chosen path on the radial line for stress distribution.
Figure 3. The process simulated under the same condition on ABAQUS: (a) the radial stress distribution at a radius of 82.5 mm and (b) the chosen path on the radial line for stress distribution.
Processes 11 01688 g003
Figure 4. (a) Three-dimensional view of the contact area, (b) contact area in rz plane, and (c) top view of (d) section XX′.
Figure 4. (a) Three-dimensional view of the contact area, (b) contact area in rz plane, and (c) top view of (d) section XX′.
Processes 11 01688 g004
Figure 5. Parts used in ABAQUS with dimensions in mm (a) tool holder (b) sheet (c) forming tool [32].
Figure 5. Parts used in ABAQUS with dimensions in mm (a) tool holder (b) sheet (c) forming tool [32].
Processes 11 01688 g005
Figure 6. Engineering and true stress and strain curve for Aluminum 6061 alloy.
Figure 6. Engineering and true stress and strain curve for Aluminum 6061 alloy.
Processes 11 01688 g006
Figure 7. Simulation results of formed part on ABAQUS: (a) stress distribution (MPa) in the deformed cone, (b) strain distribution in the deformed cone, and (c) thickness variation (mm) within the deformed cone.
Figure 7. Simulation results of formed part on ABAQUS: (a) stress distribution (MPa) in the deformed cone, (b) strain distribution in the deformed cone, and (c) thickness variation (mm) within the deformed cone.
Processes 11 01688 g007
Figure 8. (a) Strain vs. true distance along the path, (b) true thickness vs. true distance along the path, and (c) chosen path.
Figure 8. (a) Strain vs. true distance along the path, (b) true thickness vs. true distance along the path, and (c) chosen path.
Processes 11 01688 g008
Figure 9. (a) Axial force with a smoothened curve, (b) in plane force with a smoothened curve, (c) decomposition of plane forces, and (d) total force with a smoothened curve.
Figure 9. (a) Axial force with a smoothened curve, (b) in plane force with a smoothened curve, (c) decomposition of plane forces, and (d) total force with a smoothened curve.
Processes 11 01688 g009
Figure 10. (a) Setup for tensile testing, (b) sample under tension, and (c) tensile test sample.
Figure 10. (a) Setup for tensile testing, (b) sample under tension, and (c) tensile test sample.
Processes 11 01688 g010
Figure 11. Domes made after the Erichsen cup test.
Figure 11. Domes made after the Erichsen cup test.
Processes 11 01688 g011
Figure 12. (a) Setup with tool dynamometer, (b) tool used for sheet deformation, (c) deformed sheet.
Figure 12. (a) Setup with tool dynamometer, (b) tool used for sheet deformation, (c) deformed sheet.
Processes 11 01688 g012
Figure 13. (a) Different regions defined on the formed cone, (b) path chosen on formed cone in FEA model, and (c) thickness variation evaluated experimentally by FEA model and by sine law.
Figure 13. (a) Different regions defined on the formed cone, (b) path chosen on formed cone in FEA model, and (c) thickness variation evaluated experimentally by FEA model and by sine law.
Processes 11 01688 g013
Figure 14. (a) Variation of axial force and (b) fit curve of axial force with the instantaneous radius of the undeformed sheet as observed from the mathematical model, FEA, and experiments.
Figure 14. (a) Variation of axial force and (b) fit curve of axial force with the instantaneous radius of the undeformed sheet as observed from the mathematical model, FEA, and experiments.
Processes 11 01688 g014
Figure 15. (a) Variation of Fp and (b) fit curve of Fp with an instantaneous radius of the undeformed sheet as observed from the mathematical model, FEA, and experiments.
Figure 15. (a) Variation of Fp and (b) fit curve of Fp with an instantaneous radius of the undeformed sheet as observed from the mathematical model, FEA, and experiments.
Processes 11 01688 g015
Figure 16. (a) Variation of total force and (b) fit curve of total force with an instantaneous radius of the undeformed sheet as observed from the mathematical model, FEA, and experiments.
Figure 16. (a) Variation of total force and (b) fit curve of total force with an instantaneous radius of the undeformed sheet as observed from the mathematical model, FEA, and experiments.
Processes 11 01688 g016
Figure 17. Percentage error in (a) axial force, (b) resultant force, and (c) total force.
Figure 17. Percentage error in (a) axial force, (b) resultant force, and (c) total force.
Processes 11 01688 g017
Table 1. Details of parameters used for simulation on ABAQUS.
Table 1. Details of parameters used for simulation on ABAQUS.
Properties Used in SimulationCorresponding Values
0.2% offset yield strength (MPa)111
UTS(MPa)119
Total Elongation16.1%
Strength coefficient205 MPa
Strain hardening exponent0.2
Coefficient of Friction between tool and blank0.2
Hardening model usedIsotropic
Mesh Size6.96 mm
Mesh TypeSquare
No. of elements taken101 × 101
Simulation Run time28 h
Table 2. Chemical composition of 6061 Al alloy sheet used for experiments (wt%).
Table 2. Chemical composition of 6061 Al alloy sheet used for experiments (wt%).
ElementsAlSiFeZnCuCrMnMgTi
wt%97.3500.5100.4100.0600.4900.0510.0200.9500.050
Table 3. Tensile properties of 6061 Al sheet.
Table 3. Tensile properties of 6061 Al sheet.
PropertiesCorresponding Values
Young’s modulus (GPa)69.4
Yield strength (MPa)111.6
Ultimate tensile strength (MPa)118.2
Strain hardening exponent0.2
Strength coefficient (MPa)205
Table 4. Data obtained from the Erichsen cup test.
Table 4. Data obtained from the Erichsen cup test.
ParametersAl6061
Dome 1Dome 2Dome 3
Measured depth (HE)9.349.128.92
Dome height (Hd)9.028.828.61
Average value of HE9.13
Average value of Hd8.81
Table 5. Various parameters used for forming the axis-symmetric cone.
Table 5. Various parameters used for forming the axis-symmetric cone.
Parameters Corresponding Values
Tool speed200 mm/sProcesses 11 01688 i001
Tool diameter10 mm
Step depth0.5 mm
Initial circle radius96 mm
No. of cycles run100
Measured depth of cone47.6 mm
%ge spring back4.8%
Table 6. Comparison of peak forces and average forces obtained analytically, by FEA model, and determined experimentally.
Table 6. Comparison of peak forces and average forces obtained analytically, by FEA model, and determined experimentally.
ModelAverage Value (N)Peak Value (N)
Analytical1863.451957.15
FEA1876.691997.01
Experimentation1802.211963.41
Table 7. Comparison of peak and average forces resultant force in r-θ plane obtained analytically, by FEA model, and determined experimentally.
Table 7. Comparison of peak and average forces resultant force in r-θ plane obtained analytically, by FEA model, and determined experimentally.
ModelAverage Value (N)Peak Value (N)
Analytical 1538.631597.54
FEA1526.261570.94
Experimentation 1493.881509.53
Table 8. Comparison of peak and average value of total forces obtained analytically, by FEA model, and determined experimentally.
Table 8. Comparison of peak and average value of total forces obtained analytically, by FEA model, and determined experimentally.
ModelAverage Value (N)Peak Value (N)
Analytical 2416.702504.13
FEA1526.261570.94
Experimentation 1493.881509.53
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Singh, R.P.; Kumar, S.; Singh, P.K.; Meraz, M.; Srivastwa, A.K.; Salunkhe, S.; Hussein, H.M.A.; Nasr, E.S.A.; Kamrani, A. A Mathematical Model for Force Prediction in Single Point Incremental Sheet Forming with Validation by Experiments and Simulation. Processes 2023, 11, 1688. https://doi.org/10.3390/pr11061688

AMA Style

Singh RP, Kumar S, Singh PK, Meraz M, Srivastwa AK, Salunkhe S, Hussein HMA, Nasr ESA, Kamrani A. A Mathematical Model for Force Prediction in Single Point Incremental Sheet Forming with Validation by Experiments and Simulation. Processes. 2023; 11(6):1688. https://doi.org/10.3390/pr11061688

Chicago/Turabian Style

Singh, Ravi Prakash, Santosh Kumar, Pankaj Kumar Singh, Md. Meraz, Ashutosh Kumar Srivastwa, Sachin Salunkhe, H. M. A. Hussein, Emad S. Abouel Nasr, and Ali Kamrani. 2023. "A Mathematical Model for Force Prediction in Single Point Incremental Sheet Forming with Validation by Experiments and Simulation" Processes 11, no. 6: 1688. https://doi.org/10.3390/pr11061688

APA Style

Singh, R. P., Kumar, S., Singh, P. K., Meraz, M., Srivastwa, A. K., Salunkhe, S., Hussein, H. M. A., Nasr, E. S. A., & Kamrani, A. (2023). A Mathematical Model for Force Prediction in Single Point Incremental Sheet Forming with Validation by Experiments and Simulation. Processes, 11(6), 1688. https://doi.org/10.3390/pr11061688

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