Analysis of Springback Behaviour in Micro Flexible Rolling of Crystalline Materials

,is paper presents a constitutive modelling of the polycrystalline thin metal strip under a state of combined loading in microflexible rolling. ,e concept of grained inhomogeneity is incorporated into the classic Chaboche hardening model that accounts for the Bauschinger effect, in order to provide more precise description and analysis of the springback mechanism in the particular forming operation. ,e model is first implemented in the finite element program ABAQUS to numerically predict the stress-strain relationship of 304 stainless steel specimens over a range of average grain sizes. After validation of the developed model by comparison of predicted curves and actual stress-strain data points, it is further applied to predict the thickness directional springback in microflexible rolling of 304 stainless steel strips with initial thickness of 250 μm and reduction changing from 5 to 10%. ,e model predictions show a reasonable agreement with the experimental measurements and have proven to be more accurate than those obtained from the conventional multilinear isotropic hardeningmodel in combination with the Voronoi tessellation technique. In addition, the variation of thickness directional springback along with the scatter effect is compared and analysed in regard to the average grain size utilising both qualitative and quantitative approaches in respect of distinct types of data at different reductions.


Introduction
Rolling technology has been investigated extensively in the last few decades, during which theory and research have been formulated by using a combined analytical and numerical approach [1,2]. However, thanks to the increasing popularity of miniaturisation in modern industrial development over recent years, it has been promising to look into the flexible rolling process of micrometer thick materials, whereof the overall dimensions as well as the longitudinal thickness profiles of the workpieces may be adapted to the requirements of the microapparatus and systems. Pioneering investigations into the microflexible rolling process have been conducted with respect to the rolling force, the elastic recovery in thickness direction of the workpiece, the surface roughness evolution related to different process variables, and so forth [3][4][5][6]. Due to the thinness of the material, the thickness directional springback of the rolled workpiece, mainly caused by the elastic redistribution of the internal stresses upon removal of the rolling force, can lead to failure to produce a predefined thickness profile along rolling direction [7]. Consequently, an appropriate control for the springback is worthy of more attention and deeper examination.
One basic approach to the regulation of springback attempts to reduce or compensate for it through modification of process parameters or adjustment in tool design. For instance, Anggono et al. [8] developed a hybrid algorithm to rectify the springback error by iterative comparison between the deformed and target shapes, which provided a faster correction with a smaller deviation from the desired dimension during a U-channel forming case. Wang et al. [9] presented an incremental bending methodology to minimise the springback error in air bending process by utilising the data gathered from loading-unloading cycles, based on which the workpiece thickness, and material properties could be more accurately assessed and then used to calculate the final punch position to achieve the target bend angle. An irritation with a high-power diode laser was deployed by Gisario et al. [10] to control the extent of the springback by selectively heating the bent area in V-shape bending of aluminium alloy sheets. With regard to the flat rolling process, automatic gauge control technology has been widely adopted to create real-time adjustments of the roll gap in accordance with the feedback of the workpiece exit thickness, so as to correct the difference between the desired and actual thicknesses induced by springback [11][12][13]. Nevertheless, the cost and unpredictability of springback compensation can be reduced providing that the springback potential is recognised prior to the actual manufacturing [14].
As a consequence, springback prediction is becoming progressively more significant in metal forming processes and has been increasingly dependent upon numerical modelling and simulation techniques in recent years. For instance, Ankėnas and Barauskas [15] presented numerical evaluation of springback of sidewall and flange under various process conditions in U-shaped deep drawing with mild steel materials. Wagoner and Li [16] estimated the suitability of different numerical integration schemes as well as the number of through-thickness integration points in terms of springback accuracy in a beam bending model, which was reported to vary with material properties, R/t values, and sheet tensions. Sumikawa et al. [17] led a probe into the influence of material behaviour on the accuracy of springback analysis, and then developed a material model with a consideration of elastoplastic anisotropy and Bauschinger effect to further improve the accuracy of springback prediction during the hat-shaped forming using high-strength steel sheets. Furthermore, Wang et al. [18] proposed an analysis model based on the surface layer model and composite model to predict the unloading springback of copper alloy sheets in microscaled U-bending process, taking into account the interactive influence of geometry size, grain size, and punch radius. e pattern of the size effects affecting the springback was disclosed by the numerical simulations, physical experiments, and the comparison in-between. Diehl et al. [19] performed finite element simulations to reveal the variation of springback angle in relation to different material grain sizes and strain gradients in free bending of thin aluminium foils with thickness to the order of micrometers. Jiang and Chen [20] numerically investigated the grain size effect on the springback behaviour of microtubes in press bending process, and they discovered that the springback amount decreased with a decrease in thickness to average grain size (T/D) ratio for microtubes with a constant wall thickness and varying average grain size, whereas it decreased with an increase in T/D ratio for microtubes with constant average grain size. Since the springback phenomenon is closely related to the residual stress after the part has been released from the forming tools, a proper characterisation of residual stress may enhance the accuracy of springback prediction, which nonetheless highly depends on the modelling of material behaviour during the forming process [21].
For decades, numerous investigators have observed the constitutive behaviour of different materials subjected to an assortment of forming operations involving cyclic loading or complex loading conditions, under which the Bauschinger effect has widely been acknowledged, namely, the reduction of yield stress upon load reversals, which is usually accounted for by the release of dislocations piling up at the grain boundaries during initial deformation [22]. Prager [23] and Ziegler [24] first developed the linear kinematic hardening models to describe the Bauschinger effect, and the difference between their models was reflected by the translation direction of the yield surface. Later, Armstrong and Frederick [25] and Chaboche [26] modified these linear models into nonlinear ones by adding an extra term so as to capture the transient behaviour and reproduce the ratcheting in fatigue. Since then, increasing and diverse demands have encouraged improvements in either the Armstrong-Frederick model or the Chaboche model, such as to control the evolution of the decomposed kinematic hardening rules by Ohno and Wang [27,28], define effective quantities in stress and plastic strain rate by Chung and Richmond [29], or incorporate the permanent softening along with the Bauschinger effect and transient behaviour by Chun et al. [30,31]. Besides, two-surface models have been formulated by Krieg [32] and Dafalias and Popov [33] to define a continuous variation of hardening between these two surfaces, namely, that the evolution of the inner surface characterises the transient response of the material, whereas the evolution of the outer surface represents the long-term response of the material such as the Bauschinger and softening phenomena. McDowell [34] introduced a twosurface stress space model by incorporating the effects of changes in plastic strain range and nonproportionality of loading in the evolution of isotropic hardening in order to enhance the prediction accuracy of complex nonproportional deformation behaviours like cyclic ratcheting, mean stress relaxation, etc. D. Kim and J. Kim [35] developed a two-surface model for rate-dependent plasticity by combining both isotropic and kinematic hardening rules. is model was validated through both monotonic and cyclic loading cases using structural steels, where simulation results exhibited an excellent agreement with experimental values in terms of the maximum stress and shape of hysteresis. Hashiguchi [36] incorporated a subloading surface within the inner surface of the two-surface model so as to express the smooth elastic-plastic transition, avoiding an abrupt change of the strain rate/stress rate relationship when loaded from the stress state within the inner surface. Moreover, a multisurface model was proposed by Mroz [37] to combine properties of isotropic and kinematic workhardening via a configuration of surfaces of constant workhardening moduli, which specified an expression for the general features of the cyclic loading behaviour as well as a smooth transition from the elastic to plastic deformation. ese models have been set up and adopted primarily to describe the constitutive characteristics of materials at the macroscale; nevertheless, the so-called size effects, including geometry size, grain size, and grained inhomogeneity effects, cannot be negligible when constructing constitutive models will be applied on the microscale. For instance, Lu et al. [38] established a microscale constitutive model on the basis of grain size, shape, deformability, and specimen dimension to discuss the deformation behaviours of polycrystalline copper grain with different grain and feature sizes in microforming processes, like microcompression and microcross wedge rolling. According to the dislocation theory, Jiang et al. [39] proposed a constitutive model with consideration of forming temperature, Hall-Petch relationship and surface layer model to investigate the influences of temperature and grain size on the deformability of pure copper in microcompression process, where both factors were found to significantly affect the shape accuracy of the produced product, as well as the metal flow behaviour due to material heterogeneity. Liu et al. [40] developed a novel constitutive model to characterise the effects of grain size and boundary and geometry size on the flow stress of pure copper foils as well as on the microextrusion of CuZn30 alloy billets. Utilising the surface layer model, they divided the specimen into two portions, viz., surface layer and inner portions, while the grains of the inner portion were further divided into two portions, viz., grain interior and grain boundary. Likewise, Wang et al. [41] built a comprehensive constitutive model on the basis of surface layer model with consideration of the influences of grain size, material thickness, grain number through thickness, and surface property. e developed model was employed to predict the mechanical properties of thin sheet metal, and the results were compared with the physical experimental results, which confirmed that their model was particularly suitable for thin sheets with one or several grains across the thickness.
In this study, a nonlinear elastic-plastic constitutive model incorporating the two yield surface constitutive model as well as the Chaboche kinematic hardening rule is developed with the aim to give a more accurate mathematical description of the anisotropic workpiece responding to the complex loading during microflexible rolling process, thus enhancing the prediction accuracy of thickness directional springback after unloading. e newly constructed model is integrated into the ABAQUS geometric modelling, utilising the user-defined subroutine user material (UMAT) as an interface. Firstly, this model is validated by simulating the miniature tensile test and comparing results with experiments using 304 stainless steel specimens. Secondly, different constitutive/numerical models are employed to characterise the mechanical behaviour and predict the thickness directional springback during microflexible rolling of 250 µm thick 304 stainless steel strips with reduction of 5 to 10%. irdly, laboratory-scale microflexible rolling experiments are carried out to further verify and validate this novel model which produces more accurate predictions than the multilinear isotropic hardening model and the Voronoi tessellation based model. Finally, both simulation and experimental results of thickness directional springback are quantitatively revealed with respect to the average grain size, and the scatter effect caused by grained inhomogeneity is analysed synchronically in regard to the average grain size in a qualitative manner.

Elastic-Plastic Mixed-State Model.
It has been stated that the diversity of grains leads to the heterogeneous mechanical response of the whole granular material, viz., that grains with low yield strength may have suffered plastic deformation while those with higher yield strength are still in their elastic realm at a certain applied stress [5,6]. e yield anisotropy of this kind can be abstractly described using two yield surfaces metaphorically termed "initial yield surface" and "complete yield surface", which specify the transition from elastic deformation to mixed elasticplastic deformation and that from mixed elastic-plastic deformation to plastic deformation, respectively. As von Mises yield criterion has been adopted to analyse the constitutive behaviour of the rolled material, these two yield surfaces can be generally expressed by (for the planestress condition, , in which α ij and α ij ′ are the position tensors of the centres of the initial yield surface F 1 and complete yield surface F 2 in the stress space, respectively, associated with the deformation history of and the hardening rule for the material, whilst k 1 � 1/3σ 2 s (ε p ) and k 2 � 1/3σ ′2 s (ε p ), in which σ s and σ s ′ represent the sizes of these two yield surfaces, depending on the equivalent plastic strain where dε p ij denotes the components of the plastic strain increment. e relationship between the size of the yield surface and the effective plastic strain may be acquired by the uniaxial tensile test. e two surfaces are assumed to be of the identical form, having nearly parallel normal vectors at points σ ij and σ ij ′ on F 1 and F 2 , respectively.
In accordance with the applied stress, there are probably three stress states in the material, to wit, pure elastic, mixed elastic-plastic, and complete plastic. e constitutive laws for each stress state are described as follows: (a) For the pure elastic state: e sizes and centres of F 1 and F 2 remain unchanged for the applied stress, σ is within the initial yield surface F 1 , and only elastic strain occurs during this state. In the light of generalised Hooke's law, the relationship of stress to elastic strain is expressed as

Advances in Materials Science and Engineering
where D e is the constant elasticity tensor. Note that Equation (3) remains applicable for the elastic states of the material during mixed elastic-plastic loading and plastic loading. (b) For the mixed elastic-plastic state: During the mixed elastic-plastic loading, the size and centre of F 2 maintain their originality whereas F 1 keeps its size unaltered, but its centre evolves (i.e., the Bauschinger effect) until these two yield surfaces make contact with σ ij brought to coincide with σ ij ′ , which indicates the commencement of plastic loading. Based upon Ziegler's linear kinematic hardening rule, the evolution of the centre of F 1 (the back-stress α ij ) can be written as [24] dα where dμ is dependent on the material. Equation (4) states that F 1 moves in a translation in the direction of the vector connecting the two stress points σ ij and σ ij ′ . Due to the coexistence of both elastically and plastically deformed grains, a weighted heterogeneity coefficient ξ [5] is adopted to reflect the proportion, namely, that ξ/2 represents the ratio of strain occurring in the elastically deformed grains, whereas (1 − ξ/2) is formulated to represent the ratio of strain occurring in the grains that have entered the plastic region. From this point of view, the stress-strain relationship in the mixed elastic-plastic state follows that where ξ belongs to the interval [0, 2], ε e2 and ε p2 are the elastic and plastic parts of the strain occurring in the plastically deformed grain, respectively. (c) For the complete plastic state: After the material undergoes plastic deformation, it is postulated that the size of F 1 expands but its centre remains invariable, while both the size and centre of F 2 evolve following the Chaboche kinematic hardening rule, taking into account the complex loading conditions, such as compression and torsion, that may exist in the actual microflexible rolling process [43][44][45].
According to the basic form of the Chaboche model, the strain is partitioned into elastic and plastic parts, which, in the small strain hypothesis, can be formulated by [26] ε � ε e + ε p , (6) and Equation (3) still illustrates the relationship between the stress and elastic strain during the plastic loading. For yield surface F 2 , the kinematic hardening evolution of the back-stress α ′ is composed of three nonlinear terms as follows [46,47]: where a 1 , a 2 , a 3 , C 1 , C 2 , C 3 are material parameters, ε p and ε p are the plastic strain and effective plastic strain, respectively. e evolution of isotropic hardening takes the following form [46]: where Q defines the maximum change in the size of F 2 , and b indicates the rate of change in the size of F 2 as plastic deformation develops. en the change in the size of F 1 can be evaluated by [48] where k 1 ≪ k 2 , ensuring that no intersections appear between the two surfaces.
Firstly, it is assumed that the stress, back-stresses, elastic strain, and plastic strain are obtained as σ (n) , α (n) , α ′ (n) , ε (n) e , and ε (n) p at the nth step, respectively, as well as a strain increment Δε is given for the (n + 1)th step.
Neglecting the plastic part in the strain increment Δε, the initial value of the stress at the (n + 1)th step can be estimated by which is afterwards employed to determine the current stress state of the material: e stress state is pure elastic. Updating the stress and back-stress inside the initial yield surface gives e stress state is mixed elastic-plastic. e renewed stress is reckoned by e back-stress update is executed in the following manner: where Δμ is evaluated through the steps given below: (i) Adopting the von Mises yield criterion for F 1 produces where Σ x , Σ y , Σ z , T xy , T yz , T zx are the stress tensor components and where (ii) Substituting Equations (16) and (17) in Equation (15) yields where e updated back-stress α (n+1) is thus obtained by substituting Equation (20) in Equation (14).
e stress is in the complete plastic state. e stress renewal is achieved by ereby the updated back-stress and size of F 2 are, respectively, expressed as e update in the size of F 1 is subsequently made available from Figure 1 presents the schematic diagram of the constitutive model with a consideration of grained inhomogeneity.

Parametric Values of the Constitutive Model.
e 250 µm thick 304 stainless steel specimens with different average grain sizes have been utilised to assess the performance of the developed constitutive model. e mechanical properties as well as the weighted heterogeneity coefficient of each specimen have been determined by miniature tensile tests and microhardness tests, respectively. Note that the basic principles and practical procedures of both types of tests can be referred to Ref. [5]. Table 1 lists the average values of the properties identified for the experimental materials.
e parametric values in regard to the Chaboche kinematic hardening model involved in the constitutive modelling were adopted from Badnava et al.'s study [46] and are enumerated in Table 2. Note that in this study, the Chaboche model was utilised to describe the nonlinear Advances in Materials Science and Engineering kinematic hardening behaviour of the material with grains that had entered plastic deformation regime. As the microinhomogeneity e ect induced by di erent stress states of individual grains "disappears" at this stage, the material was considered uniform and treated equivalently with a development of plastic deformation at the macroscale.
erefore, macroscopic Chaboche model parameters have been adopted for numerical calculation of back-stress during plastic deformation. Another point worthy of note is that the six groups of materials were presumed to share the same parametric values involved in the Chaboche kinematic hardening model to reproduce the macroscopic behaviour in the regime of plastic activity; nevertheless, they were assigned with various elastic moduli and dissimilar weighted heterogeneity coe cients to distinguish them from one another for the time when they were subjected to the stress states of pure elastic and mixed elastic-plastic at the microscale, respectively.

Validation of the Constitutive Model.
In the developed constitutive model, size e ects were principally integrated into the mixed elastic-plastic regime by introducing a weighted heterogeneity coe cient to specify the proportions of elastically and plastically deformed grains, respectively. As can be read o from Table 1, variation in this    Figure 2 displays the stress-strain curves predicted by the elastic-plastic mixed-state model and the stress-strain data points obtained from miniature tensile tests for 304 stainless steel specimens S 1 -S 6 with average grain size of 22-74 µm. Generally, the ductility, the yield stress, and the ultimate tensile stress all increase with the decrease of the average grain size because of the enhancement of the grain boundary strengthening e ect as well as the change of the fracture mechanism [53]. It can be observed from Figure 2 that the predicted curves almost coincide with the experimental points except for several small deviations, which indicates that the developed constitutive model is valid and can be employed for the numerical prediction of springback behaviour in micro exible rolling.

Numerical Simulations of Micro exible
Rolling. e halfsymmetric 3D micro exible rolling nite element model has been established using the commercial FEA software package, ABAQUS/CAE. As exhibited in Figure 3, in order to achieve higher computational e ciency with less computational resource, the inner portion of the roll body has been removed as only its outer surface will be in contact with the material during the real forming process [54]. e basic geometric parameters of the FE model are as follows: outer diameter of the roll D 25 mm, length of the roll B 35 mm, length of the strip L 40 mm, width of the strip b 10 mm, and initial thickness of the strip H 250 µm.
One contact pair was set up in the nite element model, namely, the outer surface of the roll body and the top surface of the strip, which was de ned as frictional interface adopting a friction coe cient of 0.15, as well as applying the nite sliding formulation and surface-to-surface discretisation approach to obtain a stress accuracy improvement [5,54]. ree phases were involved in the simulation of micro exible rolling process, to wit, 5% thickness reduction for the rst rolling phase to form the thicker area, which was followed by the thickness transition area wherein thickness reduction increased swiftly from 5 to 10%, and then the third rolling phase with 10% thickness reduction to form the thinner area (i.e., the amount of thickness reduction is 12.5 μm for thicker area while 25 µm for thinner area) [5].
For comparison purposes, both the elastic-plastic mixedstate (EPMS) and multilinear isotropic hardening (MISO) models have been employed to describe the constitutive behaviour of the material and these simulation results are labelled as "Category 1" (obtained utilising the EPMS model) and "Category 2" (obtained utilising the MISO model), respectively. Furthermore, the 3D Voronoi tessellation technique has been exploited to further visualise the inhomogeneous microstructure of the workpiece which comprises grains that have dissimilar properties assigned to them. e basic concept of this method can be understood by distributing a set of Voronoi sites randomly within the workpiece volume, which partition the space into regions such that all locations within one certain region are closer to its generating site than to others [55].
is process was realised with the aid of MATLAB, which created a script le to help accomplish the Voronoi modelling in ABAQUS/CAE. Figure 4 displays an overview of the established 3D Voronoi tessellation in the mimicked 304 stainless steel strips, where di erent colours re ect di erent grain properties that were sorted and quanti ed with respect to the results produced by both the microhardness and miniature tensile tests [5]. Note that the MISO model has been adopted to represent the constitutive relations for each grain and the simulation results obtained using such a combined modelling approach (named VMISO model) are labelled as "Category 3".

Physical Experiments of Micro exible Rolling.
Physical experiments have been set up under the FEA simulation conditions and conducted on the self-developed prototype micro exible rolling mill. In addition, all the specimens were rolled at a speed of 20 cm/min and without lubrication [5]. Exemplarily, Figure 5 shows the overviews and magni ed pro les in thickness transition areas of the micro exible rolled 304 stainless steel specimens S 1 -S 6 , in which these areas were marked red, allowing them to be located easily under the microscope. Table 3 lists the experimental results for the microexible rolled 304 stainless steel specimens for both reductions of 5 and 10%, where the thickness directional springback of the specimen was estimated by the nal thickness minus the average value of the actual roll gap. Note that the measurement of the nal thickness was taken using a micrometer and was repeated at 10 di erent locations in the test area. ese results were afterwards averaged to ensure an accurate representation of workpiece thickness after rolling [5].

Simulation Results and Discussion
Taking 304 stainless steel specimens S 1 as an example, Figures 6-8 illustrate the distribution of e ective stress in the rolled workpiece using the aforementioned three categories of constitutive/numerical models, respectively.
As can be observed from these gures, the e ective stress distributions obtained from both the EPMS and MISO models generally exhibit a symmetric pattern with a smooth, gradual stress transition within the workpieces (for instance, see enlarged areas in Figures 6 and 7, respectively), whereas the VMISO model has introduced large inhomogeneities in the form of sudden changes in the effective stress distribution (for instance, see enlarged areas in Figure 8) and this irregularity can be accounted for by the di erent mechanical response of each grain to the external stress, which means some grains have slipped into the plastic state while others still remain in an elastic state with respect to the forces applied through the rolls; nonetheless, smaller e ective stress was generated at the centre than at both edges of the workpiece regardless of what type of constitutive/numerical model had been applied, mainly because the maximum de ection of the roll occurred at its centre, consequently exerting less compression at the centre of the workpiece [5].
is study has emphasised the investigation of thickness directional springback as it determines the dimensional accuracy as well as the supplementary amount of the reduction in micro exible rolling. e simulation result of thickness directional springback was calculated as the average distance between the upper and lower surface nodes of the workpiece after leaving the roll bite minus the minimum average distance between those in the roll bite. Four sets of data were collected and arranged, and each set was further sorted into six groups of fteen specimens each according to the average grain sizes summarised in Table 1, and specimens S 1 -S 6 were representative selections for they individually had the middle values in each group of one data set [5].    +1.744e + 02 Figure 6: Distribution of e ective stress in the micro exible rolled 304 stainless steel specimen S 1 utilising the EPMS model. Enlarged area I: thicker thickness area with 5% reduction; enlarged area II: thinner thickness area with 10% reduction. Figure 9 depicts the correlation between the thickness directional springback and average grain size such that the former basically exhibits uctuations with the latter, ascribable to the variability of the ratio of average rolling force to elastic modulus of the workpiece, which nevertheless stays approximately 4.50 µm at the reduction of 10%, nearly twice as much as that of around 2.50 µm at the reduction of 5%. Another point should be mentioned here is that only one type of grain size was obtained in each group, and a small o set value was deliberately given to set four types of data apart from each other so that the scatter of thickness directional springback can be seen obviously in this gure, which expands as the impact of each single grain becomes more pronounced with the decreasing number, as a consequence of the increasing size of grains within a constant strip thickness for each data set [5].
Additionally, it is observed from Figure 9 that the Category 2 of MISO model provides generally less satisfactory predictions of thickness directional springback at di erent average grain sizes, as the nonlinear loading involving the transition of the stress state is not likely to be represented adequately despite a continuous piecewise linear approximation to the ow curve of the material, while the Category 1 of EPMS model exhibits a better agreement with the measured data as it has captured the Bauschinger e ect, traced the evolution of the back-stress, as well as incorporated the mixed elastic-plastic stress state aroused by grained heterogeneity. Likewise, in spite of making use of the Finally, Table 4 illustrates a quantitative comparison of the thickness directional springback utilising the newly established EPMS model together with the conventional MISO model and the VMISO model containing a Voronoi structure involved in the workpiece. According to the calculated errors, the EPMS model has the capacity to provide for an almost 45% increase in the numerical accuracy compared with the VMISO model and an even larger increase of approximately 74% compared with the MISO model, which proves its feasibility and demonstrates its better performance on the quantitative prediction of the thickness directional springback in micro exible rolling.

Conclusions
is study looks into the springback behaviour in microexible rolling of thin crystalline metal strips with varying longitudinal thickness by means of a combined numerical and experimental approach. e following conclusions can be drawn from this study:  Figure 9: Correlation between thickness directional springback and average grain size for Groups S 1 -S 6 : (a) 5% reduction; (b) 10% reduction.
(i) Constitutive modelling in view of the Chaboche hardening model is conducted to account for the evolution of the back-stress under complicated loading patterns as well as the grain size and grained inhomogeneity effects owing to the decreasing specimen size. (ii) e novel model is implemented in ABAQUS/ Standard by using user-material (UMAT) subroutines based on an incremental deformation approach to predict the stress-strain curves for 304 stainless steel specimens over a range of average grain sizes. e coincidence between model predictions and experimental data demonstrates the correctness of the developed constitutive model. (iii) is model is further adopted to quantitatively evaluate the thickness directional springback during microflexible rolling of 250 µm thick 304 stainless steel strips with reduction of 5 to 10%. e simulation results are compared with the experimental data with an acceptable error of approximately 1.13%, which also serves as the verification and validation of the developed model. (iv) e newly established EPMS model significantly improves the springback prediction in microflexible rolling by capturing the Bauschinger effect and tracing the evolution of the back-stress with consideration of grained inhomogeneity, which increases the prediction accuracy by roughly 74% as compared with the conventional MISO model and by nearly 45% as compared with the VMISO model with a Voronoi structure constructed in the workpiece. (v) Quantitative analysis of the variation of thickness directional springback with average grain size is conducted with respect to both reductions for four types of data, and the scatter effect affiliated with the properties of each individual grain is qualitatively estimated in relation to the average grain size.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interests.