University of Birmingham Reconfigurable multipoint forming using waffle-type elastic cushion and variable loading profile

There is an increasing demand for flexible, relatively inexpensive manufacturing techniques that can accommodate frequent changes to part design and production technologies, especially when limited batch sizes are required. Reconfigurable multi-point forming (MPF) is an advanced manufacturing technique which uses a reconfigurable die consisting of a set of moveable pins to shape sheet metal parts easily. This study investigates the use of a novel variable thickness waffle-type elastic cushion and a variable punch loading profile to either eliminate or minimise defects associated with MPF, namely wrinkling, thickness variation, shape deviation, and dimpling. Finite element modelling (FEM), analysis of variance (ANOVA), and the response surface methodology (RSM) were used to investigate the effect of process parameters pertaining to the cushion dimensions and type of loading profile on the aforementioned defects. The results of this study indicate that the most significant process parameters were maximum cushion thickness, cushion cut-out base radius, and cushion cut-out profile radius. The type of loading profile was found to be insignificant in all responses, but further investigation is required as the rate, and the thermal effects were not considered in the material modelling. Optimal process parameters were found to be a maximum cushion thickness of 3.01 mm, cushion cut-out base radius of 2.37 mm, cushion cut-out profile radius of 10 mm, and a “linear” loading profile. This yielded 0.50 mm, 0.00515 mm, 0.425 mm for peak shape deviation, thickness variation, and wrinkling, respectively.


Introduction
Mass customisation gains a continuous demand as a promising approach to combine personalization and flexibility of custom-produced parts with the low mass manufacturing costs.However, mass customisation conflicts with the current production lines, which produces identical parts with high quantities.
Additive manufacturing is a technology to manufacture highly customised components but with small quantities.The technology has achieved rapid worldwide popularity because of its ability to manufacture parts with high geometric freedom and material utilisation.The technology enables the processing of a wide range of materials, including, polymers [1], ceramics [2][3][4][5][6], metals [7], and composites [8] which promotes the adoption of this technology in healthcare [9], defence [10], energy [11], and aerospace [12,13].However, for large and curvilinear sheet metals surfaces as in the automotive industry, additive manufacturing is not state of the art.The poor surface roughness, the slow-building rate, especially with large parts, and the anisotropy properties of the fabricated parts along with the need for postprocessing steps are undesirable in automotive [14].
Traditional sheet metal manufacturing methods involve plastically deforming a metallic sheet using a set of complementary dies configured to a designated geometry.These methods are widely used in large-scale production to manufacture high-quality products quickly and inexpensively.However, due to the high tooling costs and time expense associated with traditional methods, they are suboptimal where limited batch sizes are required.In recent years, the requirement for flexible manufacturing processes which can accommodate frequent changes to the part design with minimal expense has risen drastically [15,16].Multi-point forming (MPF) is one such process; it replaces its solid dies for an ordered set of discrete pins which can generally move to the workpiece to construct a pseudo-die surface [17].
MPF has seen considerable progress with regards to its viability as a manufacturing process, as well as in the removal of defects associated with it, namely dimpling, wrinkling, and springback [18].Several investigations have studied the effect of pin tip and shape on surface quality [19][20][21].Schuh et al.
[21] recommended partial spherical pin tips based on machinability and formability but reported that hemi-ellipsoidal pin tips covered the largest range of contact angles.Walczyk and Hardt [22] recommended square-based pins as only they offered load path isolation.They also investigated other MPF design factors, including pin clamping and containment, forming force capacity, and die surface formation.Park et al. used design of experiment approach to investigate the effect of process parameters on stress distribution, forming force, and spring back [23].
Paunoiu et al. [24] developed bespoke FE models for MPF based on the pin contact points.They concluded that localised deformation is significant in MPF and is heavily dependent on contact points, and hence, recommended using an interpolator between the die and cushion to improve surface quality.
Similarly, Zhang et al. [25] proposed a variation of the established MPF method, namely multi-point sandwich forming (MPSF) to reduce process defects.Gorgi et al. [26] investigated the effect of inhomogeneities on the plastic strain of metal parts.The results showed that the presence of inhomogeneities enables an accurate estimation of localized necking.In MPSF, the punch pin matrix is replaced by a deformable die and a polyurethane interpolator, and FEM was used to model the stress distribution and springback.They concluded that the presence of a pliable interpolator could produce smooth surface quality and reduce the amount of dimpling seen.Zhong-qin et al. [27] proposed an optimisation algorithm for blank holder force (BHF); their algorithm would vary the force with punch stroke.
The variable BHF was found to improve the forming limit of the workpiece by circa 30%.Liu et al. [28] investigated the effect of a novel layered blank holder that deformed in tandem with the workpiece.They reported that the design eliminated wrinkling, improved stress and strain distribution homogeneity, and thickness variation across the workpiece.Qu et al. [29] implemented a segmented strip steel pad which was located between the dies and the elastic cushion.They analysed the impact to final part quality experimentally and with FEM.In both cases, they observed that the presence of the steel pad increased friction force and generated a surface compressive stress on the workpiece.This reduced the wrinkling, dimpling, springback, and straight edge defects in the final part.Quan et al. [25] and Zareh-Desari et al. [30] investigated the effect of the elastic cushion, with the former also studying the impact of cushion thickness.They both found that the presence of the elastic cushion is necessary to minimise dimpling defects and improve forming accuracy.experiments along with a response surface methodology to study several factors such as speed, feed rate, and coolant on the surface roughness, thickness reduction, and hardness of Aluminium sheet [33].Elgahwail et al.
[34] employed the response surface DOE and the analysis of variance to identify the optimsed process parameters of MPF process on the amount of springback.the effect of coefficient of friction, pin size, cushion thickness, and radius of curvature response surface method in order to minimise process defects and improve final part quality.To the best of the author's knowledge, although the effect of conventional and mesh-type elastic cushions have been investigated, no study currently exists for using a waffle-type cushion with a variable thickness profile.Furthermore, an investigation into different punch loading profiles has not yet been conducted.This work aims to employ FEM and the face-centred response surface method (RSM) to investigate the effect of maximum cushion thickness, cushion cut-out base radius, cushion cut-out profile radius, and punch loading profile on final part quality.The quality characteristics that will be considered are the thickness variation, peak shape deviation, and wrinkling.

Materials Properties
A steel sheet made of DC05 with a thickness of 1 mm was employed in this work.DC05 is popular non-alloy steel that is used for cold forming techniques of complex shapes parts such as deep drawing and incremental forming.The composition of the material is shown in Table 1.The material was assumed to be isotropic and homogeneous, and the elasticplastic model was used.Flow stress was assumed to obey a reduced Hollomon power law: Where  refers to true stress,  is the strain hardening exponent,  is the true strain, and  is the strength coefficient.
The parameters of the reduced Hollomon power law ( and ) were found by

Numerical Modelling of MPF
A finite element model was developed based on the MPF tool in Fig. 2   The compression results were compared with Mooney-Rivlin model according to Eq2.A good agreement was found between the two models, see Fig. 3.
Where  is the strain energy density,  ̅ and  ̅ are the first and second invariants of the deviatoric strain tensor,  is the elastic volume ratio for isotropic thermal expansion,  and  are coefficients relating to deviatoric response, and  is a coefficient relating to the volumetric response obtained from a uniaxial compression test conducted using a Shore hardness of 90.
Where,  and  refer to the initial bulk and shear moduli of the material, respectively.The material was assumed to be incompressible ( ≅ 0.5) due to lack of material data with evidence to the contrary.However, due to numerical stability constraints, true incompressibility cannot be directly modelled in ABAQUS/Explicit.Thus it was assumed that the material was almost incompressible such that  1 0 and  / 20.This corresponded to a Poisson's ratio of 0.475.
An overview of the FE model is shown in Fig. 4. A general contact algorithm was used to define interfacial contact.A friction coefficient of 0.1 was assumed between all bodies to model tangential behaviour, and this was achieved using penalty formulation [35][36][37].Normal behaviour was not considered in the model.The pin matrices, namely the punch and die, were modelled as rigid bodies and were meshed using R3D4 elements.In order to validate the model against previous work from our research group, the workpiece and cushion were defined as deformable and were meshed using C3D8R elements [34,35,38].Also, shell elements, S4R, were used to mesh the workpiece as their computational cost is small, and the thickness distribution could be gauged more easily.Mesh sensitivity analysis was conducted to identify the correct element size at which solution can be reached in a reasonable time without the model being mesh dependant.The number of elements in the punch, die, and cushion were 30900, 30900, and 11781, respectively.The workpiece consisted of 47586 and 15862 for continuum-solid and shell elements, respectively.Symmetric boundary conditions corresponding to the X and Z directions were applied to the workpiece and cushion, see Fig.In this work, "sigmoid" and "linear" loading profiles were studied; these are shown in Fig. 5 along with the overall shape and dimensions of the waffle-type cushion investigated.As can be seen in Fig. 5, the waffle-type cushion consisted of a set of ordered curved cut-outs at locations directly typical to the pin travel paths.These cut-outs were defined using two main dimensions, the spherical curvature of the cut out (henceforth referred to as the cut-out profile radius) and the width of the cut-out (henceforth referred to as the cut-out base radius).As stated above, ABAQUS explicit solver was employed, this was done to avoid convergence issues due to the non-linear deformation, a large number of elements, and the complex contact conditions involved in this problem.All simulations were performed on an Intel® CoreTM i5-7300HQ processor at 2.50 GHz.The analyses were performed using double precision to avoid round-off errors, and parallel processing was used due to the large node count.To reduce simulation time, a mass scaling factor of 10,000 was used for both the C3D8R and S4R element models; this reduced computation time significantly without sacrificing numerical accuracy.A is the maximum cushion thickness, B is the cushion cut-out base radius, C is the cushion cut-out profile radius.(b) Linear and sigmoid punch loading profiles.

Model Validation
The FE model was validated against experimentally published results obtained, in the same research group, by Abosaf et al. [38] where a flat cushion and linear punch loading profile were used.Figure 6 shows a comparison between the simulation results of the developed model and experimental results (target) obtained by Abosaf et al. [38].Fig. 6a shows the deformed workpiece and Figs.4b-d show the comparison of the forcedisplacement, shape, and thickness distribution profiles to their respective targets.Whereas Fig. 6e shows the energy history output of the FE model.As shown in Fig. 6b, it can be observed that the forming force increases gradually up until 40 mm displacement.At which point, all pins were in direct contact with the elastic cushion, and plastic deformation has commenced.
After this point, work-hardening of the material leads to a sharp rise to 60.1 kN until the end of motion [39].This results in a percentage error of 2% when compared to the 58.9 kN target.For the shape profile, Fig. 6c, the peak deflections observed occurred at the centre of the sheet and were found to be -11.82mm and -29.14 mm for axes AA and BB, respectively.Comparing these simulation outputs to the targets of -13.17 mm and -29.35 mm yields percentage errors of 10.3% and 0.7%, respectively.
In the case of thickness distribution, Fig. 6d, the workpiece is thinner nearer its centre and becomes thicker closer to the flange.This occurs as the centre of the sheet undergoes the most significant level of plastic deformation leading to sheet thinning and outward material flow.The peak thickness was found to be 1.013 mm and 1.023 mm for axes OA and OB, respectively.The relative change in values compared to the original workpiece thickness of 1 mm is very small, indicating that the normal to longitudinal plastic strain is approximately equal, which is ideal in sheet forming processes.Comparing these outputs to the targets of 1.005 mm and 1.013 mm yields percentage errors of 0.8% and 1.0% respectively.Overall, the developed model sees a good agreement with the experimental results obtained by Abosaf et al. [38] with a maximum error of 10.3%.To ensure reliability, it was necessary to check the stability of the solution given by the FE model.This was achieved by confirming that the workpiece deforms quasi-statically and element distortion via hourglassing is kept to a minimum.From Fig. 6e, the kinetic energy (KE) and artificial strain energy (AE) at the end of the analysis total to 2.1% and 10.2% of the internal energy (IE).As the KE is less than 5% of the IE, then, inertial forces can be considered small enough to not dominate the solution [40].For FE model to be reliable, the maximum KE of the deformed material and the maximum AE must both be less than 10% of the maximum IE [41].As the KE is 2.1% of the IE and the AE is approximately 10% of the IE [42], it was sufficient to conclude that artificial deformation had minimal impact on the solution and the FE model can be considered reliable.

Statistical Validation
Statistical methods, such as DOEs and ANOVA, have been used widely in manufacturing to investigate, predict, and optimise process response to a change in process parameters.A face-centred RSM was used to generate a set of experiments for analysing the effects of maximum cushion thickness, cut-out base radius, cut-out profile radius, and the type of punch loading profile on the defects seen in MPF.First-order orthogonal response-surface methods are generally used over a narrow set of process parameters; hence a second-order polynomial model was selected.This is given by the general expression [14]: Where  is the process response,  is any of the studied process factors,  refers to the polynomial coefficients, and  refers to the random error.The coefficients were derived using non-linear least-squares analysis.In this work, for each continuous parameter, three levels were tested: -1, 0, and 1.As the punch loading profile is a categoric factor, all continuous factor experiments were repeated for each level of that categoric factor, in this case, only two levels were tested.Table 3 summarises the process parameters (and their levels) which were used in the simulations.The response variables correspond to the quality characteristics (QCs) of the final part [31].In this work, only the quantitative QCs: thickness variation, wrinkling, and peak shape deviation, were considered in the DOE.Dimpling was considered a qualitative QC and was noted as being present when a nonuniform material distribution with highly localized strain was noticeably visible in the formed part.Wrinkling was defined as the normal deviation of the formed part from the target shape seen in Fig. 6c when measured at the sheet flange.It was quantified as a root-mean-square-error (RSME) using Eq. 5, Zi is a single deviation of the formed part from the target shape [38].To observe the trend more easily, only the wrinkling in the long edge of the deformed workpiece was used in the calculations.
Eq. 5 The sheet metal thickness was measured along the principal axes, and at the workpiece flanges, this was then quantified as a standard deviation,  (Eq.6 [29]), which was termed the thickness variation.Here,  denotes the number of points where the thickness was recorded,  is a single data point, and ̅ is the mean thickness in the data set.

𝑠 𝑥 𝑥̅ 𝑁
Eq. 6 Peak shape deviation was defined as the maximum normal distance between the target and formed part shapes, as stated previously, this would occur at the centre of the workpiece.

Results and Discussion
Table 4 shows the generated plan of 40 runs based on the DOE and the evaluated response of each QC in these simulations.The simulation results were then analysed using Design Expert 12, and an ANOVA study was conducted to identify statistically significant parameters.In this investigation, both peak shape deviation and wrinkling were fitted using standard response modelling, whereas thickness variation required a logarithmic Box-Cox transformation.The coefficient of determination, R 2 , was found to be 91.66 %, 90.46 %, and 97.50 % for peak shape deviation, thickness variation, and wrinkling, respectively.Similarly, the adjusted R 2 values were found to be 87.49%, 85.68 %, and 96.25 % and the residuals were all approximately normally distributed, indicating good agreement with the quadratic model.the change in punch speed which would be of a small consequence in a quasi-static process.Further to this, the literature finds that the strain rate sensitivity, , for DC05 steel is not insignificant, varying from 0.023 at low rates to 0.130 at high rates [44,45].This argument is also supported by the fact that the p-values for the interactions between the type of loading profile and cut-out base radius in thickness variation and wrinkling are very close to the 0.05 threshold and thus would likely become significant if the aforementioned effects were included.Hence, further investigation into the deformation behaviour of DC05 steel in MPF under various strain-rates is likely warranted.

Peak shape deviation
Fig. 7a shows the surface plot for the effect of maximum cushion thickness and cut-out base radius on the peak shape deviation.It can be observed that as the maximum cushion thickness increases, so too does the peak shape deviation.This is expected and is consistent with similar studies [34,38].
Increasing cushion thickness reduces local deformation and surface indentation.The larger material volume means that it can attenuate punch impact energy more effectively.This results in a reduced, more homogeneous pressure distribution as can be seen in Fig 7d, and an overall reduction in local sheet thinning.However, the lower stresses have the added drawback of leading to under-deformation of the workpiece and thus increased shape deviation [15,38].A minimum error from the ideal shape is achieved when the maximum cushion thickness is 3 mm.on the upper workpiece surface when using a 9 mm cushion.(e) Equivalent strain contour (S4R) on upper workpiece surface when using a 3 mm cushion with 5.13 mm cut-out base radius and 10 mm cut-out profile radius.
Although cut-out base radius is not statistically significant as a linear factor in peak shape deviation, it is as a quadratic one.It can be observed that the deviation is low when using base radii of 2.37 mm and 5.13 mm and is at its highest when an intermediate radius of 3.75 mm is used.At small radii, the variation in thickness across the cut-out is small, and so deformation behaviour is similar to that of a flat cushion.When the cushion is compressed, the surface adjacent to the pin matrices deforms, leading to a series of depressions at the point of contact [25].Flatter cushions will exhibit smoother deformation during this process as their geometry is simple.As the base radius increases, cushion material flow in these regions becomes increasingly more complex, which may lead to non-uniform deformation of the workpiece and greater deviation.However, it is also true that increasing the base radius reduces the local cushion thickness, which acts to increase the transmitted punch contact pressure.This, in effect, will lead to better deformation and a smaller deviation.Thus, a potential explanation for the trend observed in Fig 7a, is that below 3.75 mm radius, the former effect is dominant, and beyond this point is when the latter phenomenon becomes dominant.If the cushion is thin; however, the latter phenomenon can also lead to the formation of the dimpling defect, this can be observed from the discontinuous, highly localised equivalent strain regions seen in Fig. 7e.
Fig. 7b shows the surface plot for cushion thickness and cut-out profile radius on peak shape deviation.It can be observed that as the profile radius decreases so too does the peak shape deviation with the minimum error being achieved at 10 mm curvature.One explanation for this trend pertains how the profile radius affects the contact conditions between the pins and the cushion.
Fig. 8 shows the effect of using small and large profile radii on pin-to-cushion contact and the developed stress distribution.It can be observed that as the profile curvature decreases, the pin-to-cushion contact area is increased, with the maximum area being achieved when the pin curvature matches that of the cut-out.This achieves a similar effect to increasing cushion thickness, albeit to a much smaller degree, whereby the pressure distribution becomes more even.This is seen in the equivalent stress contours where the size of local low-stress concentrations decreases at smaller profile radii, with this effect being most observable at the sheet corners.These low-stress regions will elastically recover whilst the surrounding material, which is plastically deforming, will not, causing uneven workpiece deformation to take place.
Hence, the reduction in the size of these concentrations means that the stress is more uniform and local deformation is reduced.Additionally, decreasing the profile radius will not suppress macro-scale deformation of the entire workpiece like in the case of increasing maximum cushion thickness as there is no significant loss in the overall stress.Thus the overall deformation improves, thereby reducing deviation.

Sheet metal thickness variation
Fig. 9a shows the surface plot for the effect of cushion thickness and cut-out base radius on the thickness variation.It can be observed that as the maximum cushion thickness increases the sheet metal thickness variation decreases, with the minimum variation being achieved at 9-mm cushion thickness.As mentioned, thicker cushions lead to under-deformation; it follows that this results in reduced material flow outward from the workpiece centre, meaning the thickness across the sheet is more uniform overall.It should be noted that this differs from some existing findings in the literature [34].However, these can be attributed to some numerical modelling differences such as friction, and that these studies limited measuring thickness to only the principal axes of the workpiece, whilst in this work thickness at the workpiece flange was also considered.Moreover, the result found in this work is consistent with the results obtained by Abosaf et al. [15], so it can be considered to be reliable.
It is also observed that the sheet metal thickness variation increases as the cut-out base radius increases, with minimum thickness variation being achieved when a base radius of 2.37 mm is used.Figs.9c-d shows the thickness distribution contours when using small and large base radii.It can be observed that a larger base radius results in a larger sheet thinning region.This is attributed to the aforementioned larger contact pressures at larger base radii generating more sheet stretching.Fig. 9b shows the surface plot for the effect of cushion thickness and cut-out profile radius on the sheet metal thickness variation.As can be seen, the thickness variation decreases with decreasing profile radii curvature, with the minimum thickness variation being achieved at a 10 mm profile radius.As stated earlier, the smaller contact area and less uniform stress distribution at larger profile radii promote an increase in local deformation and cause the workpiece to deform more unevenly in regions where there are large differences in stress.This naturally leads to a more non-uniform workpiece thickness distribution as regions of high stress will exhibit more thinning than those of low stress.

Flange wrinkling
Wrinkling arises when in-plane tensile forces are insufficient, this can generate out-of-plane deformation in the form of wave-like perturbations.
These are due to local plastic deformation that occurs when some of the pins in the upper and lower dies starts to establish contact with the sheet during the deformation process.The force starts to increase rapidly when all pins establish contact with the sheet until the maximum plastic deformation is reached.Fig. 10a shows the surface plot for the effect of maximum cushion thickness and cut-out profile radius on wrinkling.It can be seen that wrinkling increases with cushion thickness, with maximum wrinkling being obtained at a thickness of 9 mm.This result agrees with findings in the literature [15].According to Abebe et al. [46], the punch contact pressure has to exceed the induced compressive instabilities that are generated during deformation to eliminate wrinkles.When a thick cushion expands due to compression, a greater volume of material is forced to accumulate near the workpiece flange [47], this means that contact pressure in this region is reduced more so than in thinner cushions, it follows that fewer in-plane compressive instabilities are suppressed and wrinkle wave amplitude increases as a result.Fig. 10b shows the surface plot for the effect of cushion thickness and cut-out base radius on wrinkling.It can be observed that wrinkling increases with decreasing cut-out base radius.This operates in the same manner as above, where the local reduction in cushion thickness and resulting increased pressure at larger base radii provide the necessary in-plane tensile force to counteract wrinkle formation.This is seen in Fig. 8b, where the fraction of the workpiece that corresponds to the sheet thickening region decreases in size when the base radius is increased from 2.37 mm to 5.13 mm.
It can also be seen from Fig. 10a that wrinkling increases with increasing cutout profile radius, with maximum wrinkling being obtained at 20 mm curvature.
A possible explanation for this phenomenon is that the steeper contact angles at larger profile radii cause cushion material to flow in the normal direction when the cushion is compressed during the punch stroke.At the flange, this can lead to larger undulations and cause the cushion to depart slightly from the workpiece.This reduces the tangential tensile force that is exerted by the cushion (due to friction) on the workpiece in this region and promotes material to flow out-of-plane to fill the resulting departure regions inevitably exacerbating wrinkling defects [48,49].This explanation is further supported when observing the deformed displacement contours shown in Figs.10c-d, where cushion wrinkling (and hence workpiece wrinkling) at the flange is seen to be relatively larger when using a 20 mm radius compared to 10 mm radius.
It is also observed that the increase in wrinkling is rather small, especially when compared to the change in wrinkling observed when changing the maximum cushion thickness and cut-out base radius.This is because although more material does flow in the normal direction at larger profile radii, the magnitude of material flow in the lateral direction is still much larger in comparison.Fig. 10e shows the plastic strain of the workpiece when two contact pins are used.

Optimisation of process parameters
Applying Eq. 4 to the studied process parameters gives the general form of the governing equation for each response.This can be defined such that: Eq. 9 In the case of punch loading profile, "linear" was assigned a normalised value of -1 and "sigmoid" was assigned a normalised coded value of +1.Table 6 shows the values of the polynomial coefficients   .Optimal parameters for minimising the studied responses were found using numerical optimisation, these are shown in Table 7.These optimal values were then validated using the same FE model.Table 8 compares the predicted responses, within a 95% confidence interval (±1.96SD), to the measured responses.It can be observed that the predicted values are underestimated in the case of peak shape deviation and thickness variation and overestimated in the case of wrinkling.However, all the measured values fall within the 95% confidence interval of their predictions, this seems reasonable given the complex deformation behaviour associated with a variable thickness waffle-type elastic cushion.

Conclusions
In this study, FEM was used in tandem with RSM and ANOVA to investigate the deformation of a DC05 workpiece when using a variable thickness waffletype elastic cushion and variable punch loading profile in MPF.This study demonstrated the following: 1.The maximum cushion thickness, cut-out base radius, and cut-out profile radius were all significant process parameters in their effects on peak shape deviation, wrinkling, and thickness variation.In terms of importance, maximum cushion thickness was the most significant followed by cut-out base radius and then by cut-out profile radius.The quadratic model was found to be the best fit for the response variables investigated.
2. The type of punch loading profile was deemed seemingly insignificant in all cases, but it is too early to completely discount it as rate and thermal effects were not considered in the FE model, thus further investigation is required.
3. Two-way interactions between process parameters were insignificant in all cases.
4. Maximum shape deviation was found to decrease with decreasing maximum cushion thickness and cut-out profile radius.For increasing the cut-out base radius, it was found to first increase then decrease.
5. Thickness variation was found to decrease with increasing maximum cushion thickness, decreasing cut-out base radius, and decreasing cut-out profile radius.
6. Wrinkling was found to decrease with decreasing maximum cushion thickness, decreasing cut-out profile radius, and increasing cut-out base radius.
7. In all cases, the results indicate that a waffle-type elastic cushion can be used to minimise the defects associated with MPF with optimal process parameters being found.However, further experimental investigations are still required.

fitting Eq. 1
to the stress-strain curve of DC05 sheet steel obtained by uniaxial tension test using a Zwick tensile testing machine.A comparison between the material model and the experimentally obtained tensile test properties is shown in Fig 1.

Figure 1 :
Figure 1: Material model and experimental stress strain curve of DC05 steel.
using ABAQUS CAE 2018.The setup is shown in Fig. 2 consists of a set of 30 x 20 pin matrices, two elastic cushions, and a workpiece.To reduce the computational cost, only a quarter of the setup was simulated because of geometrical symmetry in the X and Z directions.The pins had general dimensions of 10 mm x 10 mm x 10 mm, with the pin tips possessing 10 mm spherical curvature.Pin separation was set at 0.25 mm.The workpiece was modelled as a DC05 steel sheet of dimensions 153.5 mm x 102.5 mm x 1 mm and was set to have a final geometry of 400 mm spherical curvature.

Figure 4 :
Figure 4: Schematic of the initial FE model, including element type and boundary conditions (BCs).
4. A boundary condition was applied to the die to constrain it in all six-degrees-of-freedom (DOF).Similarly, a displacement boundary condition was used on the punch to constrain it in XYZ rotationally and XZ translationally, see Fig. 4. Punch displacement was set at 42.2 mm in the Y direction, and an amplitude operator was used to vary the loading profile.

Figure 5 :
Figure 5: (a) Overall geometry of the waffle-type elastic cushion with a variable thickness profile.The dimensions varied in this study are A, B, and C.

Figure 6 :
Figure 6: (a) Schematic overviewing principal axes of the deformed workpiece.(b) The force-displacement output of the C3D8R FE model compared to the target profile.(c) Deformed profile output of C3D8R FE model across principal axes compared to the target profile.(d) Thickness distribution output of S4R FE model across principal axes to centroid compared to the target distribution.(e) Energy history output of the C3D8R FE model.

Figure 7 :
Figure 7: (a) Surface plot for the effect of maximum cushion thickness and cut-out base radius on peak shape deviation.(b) Surface plot for the effect of maximum cushion thickness and cut-out profile radius on peak shape deviation.(c) Pressure distribution contour (S4R) on the upper workpiece surface when using a 3 mm cushion.(d) Pressure distribution contour (S4R)

Figure 8 :
Figure 8: Comparison of the pin-to-cushion contact points and equivalent stress contours (S4R) on the upper workpiece surface when using 10 mm profile radius and 20 mm profile radius.

Figure 9 :
Figure 9: (a) Surface plot for the effect of maximum cushion thickness and cut-out base radius on sheet metal thickness variation.(b) Surface plot for the effect of maximum cushion thickness and cut-out profile radius on sheet metal thickness variation.(c) Section thickness contour on the upper workpiece surface when using a cut-out base radius of 2.37 mm.(d) Section thickness contour on the upper workpiece surface when using a cut-out base radius of 5.13 mm.

Figure 10 :
Figure 10: (a) Surface plot for the effect of maximum cushion thickness and cut-out profile radius on wrinkling.(b) Surface plot for the effect of maximum cushion thickness and cut-out base radius on wrinkling.(c) Displacement contour (S4R) of deformed workpiece and cushion (at the flange) when using a cut-out base radius of 20 mm.(d) Displacement contour (S4R) of deformed workpiece and cushion (at the flange) when using a cut-out base radius of 10 mm, (e) plastic strain of the workpiece when when two contact pins are used.

Table 1 :
Chemical composition of DC05 steel, as supplied.

Table 2 :
Mechanical properties of DC05 steel.

Table 3 :
Process parameters and their corresponding levels.

Table 4 :
DOE results for shape deviation, thickness variation, wrinkling, and dimpling as obtained when using S4R elements.

Table 6 :
Polynomial coefficients for response variable equations.

Table 7 :
Optimum conditions for minimal defects.

Table 8 :
Comparison of predicted response variables from the optimal setting and the measured values.