Effect of the Computational Model and Mesh Strategy on the Springback Prediction of the Sandwich Material

The effect of the computational model and mesh strategy on the springback prediction of the thin sandwich material made of micro-alloyed steel was investigated in this paper. To verify the chosen computational strategy, a comparison of the experimentally obtained specimen (U-bending) with the FEA result was performed. The Vegter yield criterion combined both with the isotropic and kinematic hardening law was used for the calculation. In addition, the effect of the deformation mesh element (surface and volume) on the accuracy of the springback prediction was investigated. It was concluded that the choice of the volume deformation mesh does not significantly improve the accuracy of the results. Moreover, it is quite a time-consuming approach. The much greater influence was monitored by concerning the selection of hardening law, where the anisotropic one was more suitable to be used on the springback prediction of a given sandwich material.


Introduction
In recent years, there has been a steady development in the range of materials ready to be used in various areas of industry. Such materials include, for example, so-called sandwich materials, which derive benefits from the advantageous utility properties of different material groups (e.g., combination metal-plastic). Nowadays, sandwich materials are mainly used either to reduce the overall weight of given products or due to their specific properties to reduce noise or vibrations.
The increasingly widespread utilisation of various types of sandwich materials, not only in the automotive industry, naturally also brings with it the necessity to know their specific deformation behaviour (generally material characteristics) under various loading conditions. For this reason, sandwich materials have been subjected to continuous testing and determination of their deformation behaviour, e.g., Fischer [1] performed static and dynamic tests (compression, shear, bending and impact tests) of aluminium foldcore specimens and compared the results with the finite element analysis (FEA). The effect of variable mechanical properties of the sandwich beam core in light of elastic loading was investigated by Grygorowicz et al. [2]. Moreover, other authors focused, for example, on the hardness measurement and thermal analysis of aluminium-based sandwich composites [3], mechanical, fatigue and thermal properties of the honeycomb sandwich structures [4][5][6][7][8][9][10][11][12] or properties of sandwich materials when joining technologies are taken into account. For example, Gladkovsky [13] researched a feature of interface structure at explosive welding and Nilsson [14] tested the fatigue properties of laser stake-welds. Generally, the joining of sandwich materials represents very important research and was investigated by many authors such as, e.g., [15][16][17]. Nevertheless, before the application of sandwich material, it is very important not only to verify these results but also to incorporate the chosen material and its material behaviour into the whole production process. For this reason, numerical the effect of anisotropy change at loading and springback recovery for AHSS materials. A combined isotropic-kinematic hardening model was modified in this case having a good matching with the real samples. Baara et al. [42] improved the accuracy of the finite element software (Abaqus) by introducing the new constitutive hardening model that can predict the final shape of a part that changes its dimensions after forming. The modification of the original Chord model and its combination with nonlinear isotropic-kinematic revealed a higher accuracy of this new model.
However, in light of the springback prediction, it is always important to consider the mutual effect between all major input parameters (elastic modulus, geometry, yield criterion and hardening laws). That is why quite a lot of papers deal with the effect of these parameters on the springback prediction. Mulidrán et al. [43,44] combined different yield criterions and hardening laws to determine the accuracy of springback prediction, in these cases for DP steel and aluminium alloy AA6451, respectively. Trzepiecinski and Lemu [45] studied the springback behaviour of steel DC04 under the different input parameters (geometry, conditions and material) for FEA (here, Abaqus). The work of Naofal et al. [46] presented calibrated plasticity models arising from the appropriate tests to predict springback during the roll forming process. In this case, the Yoshida-Uemori hardening law revealed a higher accuracy to predict springback than the isotropic hardening law. Sumikawa [47] developed and implemented an FEM software material model, which takes into account not only average Young s modulus and Bauschinger effect, but also elastic and plastic anisotropy. It was concluded that besides the Bauschinger effect and Young s modulus, there is also a great effect arising from elastic and plastic anisotropy on the results of springback predictions. The effects of the different yield criterions (Hill48 and Hill90), bending radii and directions on the aluminium alloys springback are described by Slota et al. [48].
In this paper, the springback prediction of the sandwich material was performed by using the FEA. Due to the previous results (see above), only the Vegter yield criterion was used and was combined with the two hardening laws (isotropic and kinematic). In addition to that, the effect of mesh strategy (selection of surface or volume elements) on the springback prediction accuracy was investigated as well. This quite complex combination (sandwich material-two hardening laws-type of mesh elements) of input parameters can be considered as the greatest novelty of this paper.

Materials and Methods
For the tests, the sandwich material CR270Y340T-LA -EG47/47-E-R-O (designation according to the standard VW 50065) was chosen. It is a micro-alloyed, electro-galvanized sheet designed for the production of car-body outer panels. Concerning the increased requirements for acoustic vibration absorption, this tested material was supplied as a sandwich with a layer of thickness as follows: 0.4-0.1-0.4 mm. The core material was identified by using the Fourier transform infrared spectroscopy (FTIR) on a device Nicolet iS10 manufactured by the company Thermo Scientific (Waltham, MA, USA). The method is based on the absorption of infrared radiation as it passes through the specimen, during which changes in the rotational and vibrational energy levels of the molecule occur as a function of the molecular dipole moment changes. The resulting infrared spectrum is a functional dependence of the energy, usually expressed in percent of transmittance or absorbance units, on the wavelength of incident radiation. For the present analysis, the attenuated total reflection (ATR) technique was used and the measured spectra of the individual materials (32 scans were used) were compared with the spectra in the instrument's database. From the results of the IR spectrometry, it could be concluded that it was an acrylate copolymer. The identification of core material was carried out for further material testing of the sandwich core mechanical properties.
The chemical composition analysis was performed on a Bruker-Quantron sparkoptical emission spectrometer Q4 Tasman. Such an analysis was determined based on 5 measurements. The results of the chemical composition analysis are given in Table 1. The subsequent structural analysis was carried out on a TESCAN MIRA 3 electron scanning microscope. From the results it is clear that the structure was a ferritic one with the dispersed particles containing an increased content of Mn, Si and Al micro-alloys. The analysed structure can be seen in Figure 1 (left) and detail of the particle is shown in Figure 2 (right). Finally, the EDX bulk analysis arising from the analysed structure is given in Figure 2.
The chemical composition analysis was performed on a Bruker-Quantron spark-optical emission spectrometer Q4 Tasman. Such an analysis was determined based on 5 measurements. The results of the chemical composition analysis are given in Table 1. The subsequent structural analysis was carried out on a TESCAN MIRA 3 electron scanning microscope. From the results it is clear that the structure was a ferritic one with the dispersed particles containing an increased content of Mn, Si and Al micro-alloys. The analysed structure can be seen in Figure 1 (left) and detail of the particle is shown in Figure  2 (right). Finally, the EDX bulk analysis arising from the analysed structure is given in Figure 2.   The chemical composition analysis was performed on a Bruker-Quantron sparkoptical emission spectrometer Q4 Tasman. Such an analysis was determined based on 5 measurements. The results of the chemical composition analysis are given in Table 1. The subsequent structural analysis was carried out on a TESCAN MIRA 3 electron scanning microscope. From the results it is clear that the structure was a ferritic one with the dispersed particles containing an increased content of Mn, Si and Al micro-alloys. The analysed structure can be seen in Figure 1 (left) and detail of the particle is shown in Figure  2 (right). Finally, the EDX bulk analysis arising from the analysed structure is given in Figure 2.

Static Tensile Test
A static tensile test was used to determine the basic material characteristics of the tested material. With a view to the future use of the experimentally obtained data for the definition of the material model in numerical simulations, these tests were already performed for specimens taken in the directions 0 • , 45 • and 90 • concerning the rolling direction. In order to eliminate any material anomalies, material characteristic values were determined each time from a set of 5 measured specimens. The shape and dimensions of the test specimens were chosen in accordance with the standard EN ISO 6892-1. In addition, the methodology of test execution and evaluation was in accordance with the standard above-mentioned EN ISO 6892-1. The static tensile test was carried out on a modernized device TIRA test 2300 using the software Labtest v.4 to evaluate all basic mechanical properties of the tested material (proof yield strength R p0.2 , ultimate strength R m , total ductility A 80mm , uniform ductility A g and Young s modulus E). The plastic strain ratio (Lankford parameter) r was determined in accordance with the standard ISO 10113 in the interval of true plastic strain values from 10% to 20%. The results of the individual material quantities are summarised in Table 2. Engineering stress R (MPa) vs. engineering strain ε (1), thus so-called engineering stress-strain curves (regarding the individual rolling directions) are shown in Figure 3. The evaluation of the mechanical values given in Table 2 was carried out by considering the measured force values to the whole cross-section, i.e., the measured values represent the average values of mechanical properties for the three layers of the tested sandwich material (steel-plastic-steel).

Static Tensile Test
A static tensile test was used to determine the basic material characteristics of the tested material. With a view to the future use of the experimentally obtained data for the definition of the material model in numerical simulations, these tests were already performed for specimens taken in the directions 0°, 45° and 90° concerning the rolling direction. In order to eliminate any material anomalies, material characteristic values were determined each time from a set of 5 measured specimens. The shape and dimensions of the test specimens were chosen in accordance with the standard EN ISO 6892-1. In addition, the methodology of test execution and evaluation was in accordance with the standard above-mentioned EN ISO 6892-1. The static tensile test was carried out on a modernized device TIRA test 2300 using the software Labtest v.4 to evaluate all basic mechanical properties of the tested material (proof yield strength Rp0.2, ultimate strength Rm, total ductility A80mm, uniform ductility Ag and Young´s modulus E). The plastic strain ratio (Lankford parameter) r was determined in accordance with the standard ISO 10113 in the interval of true plastic strain values from 10% to 20%. The results of the individual material quantities are summarised in Table 2. Engineering stress R (MPa) vs. engineering strain ε (1), thus so-called engineering stress-strain curves (regarding the individual rolling directions) are shown in Figure 3. The evaluation of the mechanical values given in Table 2 was carried out by considering the measured force values to the whole cross-section, i.e., the measured values represent the average values of mechanical properties for the three layers of the tested sandwich material (steel-plastic-steel).   A numerical simulation has to be a stress-strain curve (hardening curve) entered as a dependence of true stress (effective stress) vs. true plastic strain. Measured dependences of the true stress vs. true plastic strain for the individual rolling directions 0 • , 45 • and 90 • were subsequently approximated (fitted) using the Krupkowski law (power-law function) according to Equation (1). σ = K (ε pl + ε 0 ) n (1) where: σ-true stress (effective stress) (MPa) K-strength coefficient (MPa) ε pl -true plastic strain (1) ε 0 -offset true strain (pre-strain) (1) n-strain hardening exponent (1) An approximation of the stress-strain curve (hardening curve) was performed within the interval from 1% up to A g in the software OriginPro 9.0. The advantage of defining the hardening curve by using the Krupkowski hardening law is the simple possibility of integrating this function and thus it is easy to express the equivalent plastic work parametrically in dependence on the true strain and the approximation constants arising from such a function. An example of the approximation result can be seen in Figure 4. The results of the obtained approximation constants are presented in Table 3. where: σ-true stress (effective stress) (MPa) K-strength coefficient (MPa) εpl-true plastic strain (1) ε0-offset true strain (pre-strain) (1) n-strain hardening exponent (1) An approximation of the stress-strain curve (hardening curve) was performed within the interval from 1% up to Ag in the software OriginPro 9.0. The advantage of defining the hardening curve by using the Krupkowski hardening law is the simple possibility of integrating this function and thus it is easy to express the equivalent plastic work parametrically in dependence on the true strain and the approximation constants arising from such a function. An example of the approximation result can be seen in Figure 4. The results of the obtained approximation constants are presented in Table 3.

Hydraulic Bulge Test
When testing material by static tensile testing, the test specimen was deformed only by the uniaxial tension (uniaxial tensile stress state). To obtain information about the material deformation behaviour under the multi-axial loading, a hydraulic bulge test (HBT) was mostly performed. The principle of this test rests in stretching the specimen into the circular die by applying the pressure of fluid until a crack occurs. Similar to the static tensile test, the dependence of the true stress vs. true plastic strain (again true stress-strain curve) is evaluated. Figure 5 shows an element of the test specimen together with the parameters used to derive the stress in the test specimen during the HBT. was mostly performed. The principle of this test rests in stretching the specimen into the circular die by applying the pressure of fluid until a crack occurs. Similar to the static tensile test, the dependence of the true stress vs. true plastic strain (again true stress-strain curve) is evaluated. Figure 5 shows an element of the test specimen together with the parameters used to derive the stress in the test specimen during the HBT. Based on the right part of Figure 5, where the effects of fluid pressure on the specimen element are shown (equilibrium equation), can be acc. to membrane theory derived Equation (2) for the stress in the test specimen wall and so the effective stress σeff (arising from the Tresca's plastic flow criterion) can be finally expressed by relation (3). where: p-hydraulic pressure (MPa) R-radius of curvature (mm) t-actual thickness of specimen (mm) Finally, effective strain εeff can be calculated using principal strains ε1,2,3 according to Equation (4) where: ε1,2,3-principal strains (1) Actual thickness of the specimen can be subsequently determined through the principal strain ε3 (sometimes termed also as a thickness reduction) firstly via the constant volume law (5) and subsequently using Equation (6). Based on the right part of Figure 5, where the effects of fluid pressure on the specimen element are shown (equilibrium equation), can be acc. to membrane theory derived Equation (2) for the stress in the test specimen wall and so the effective stress σ eff (arising from the Tresca's plastic flow criterion) can be finally expressed by relation (3). where: p-hydraulic pressure (MPa) R-radius of curvature (mm) t-actual thickness of specimen (mm) Finally, effective strain ε eff can be calculated using principal strains ε 1,2,3 according to Equation (4) where: Actual thickness of the specimen can be subsequently determined through the principal strain ε 3 (sometimes termed also as a thickness reduction) firstly via the constant volume law (5) and subsequently using Equation (6). where: t-actual thickness (mm) t 0 -initial thickness (mm) To determine the hardening curve in the HBT, it is therefore necessary to measure the dependence of the test specimen deformation on the hydraulic pressure. That is why a contactless photogrammetry method using the software Mercury RT from the company Sobriety, Ltd. (Kuřim, Czech Republic) was used to measure the true strain values. A scheme of the measurement layout is shown in Figure 6.

t0-initial thickness (mm)
To determine the hardening curve in the HBT, it is therefore necessary to measure the dependence of the test specimen deformation on the hydraulic pressure. That is why a contactless photogrammetry method using the software Mercury RT from the company Sobriety, Ltd. (Kuřim, Czech Republic) was used to measure the true strain values. A scheme of the measurement layout is shown in Figure 6. Because stress relaxation occurs when the specimen is loaded, the pressure of the fluid was not increased continuously, but gradually at intervals of 5 s. The value of the pressure increase (pressure step) was chosen to be 0.25 MPa and the own magnitude of the required pressure was controlled by a precision hydraulic servo valve with PID control.
The achievable accuracy of the required pressure was 0.002 MPa. The system of hydraulic pressure activation and its control was solved by a separate hydraulic unit working independently of the hydraulic circuit of the given press. From the contactless deformation analysis using the system Mercury RT, the magnitudes of major and minor strain (ε1, ε2), as well as the radius of curvature R must be detected to subsequently determine the stress-strain curve. For this purpose, two surfaces have been defined on the measured area, in which the system calculates the desired quantities. An example of the analysis area is shown in Figure 7. The larger analysis area (area A) is used to find the radius of curvature R and the smaller area (area B) is used to calculate the major and minor strain, from which the actual thickness of the material can be calculated using Equation (5) and (6). An example of true strain distribution in the relevant area and evolution of the radius of curvature R is shown in Figure 8. Because stress relaxation occurs when the specimen is loaded, the pressure of the fluid was not increased continuously, but gradually at intervals of 5 s. The value of the pressure increase (pressure step) was chosen to be 0.25 MPa and the own magnitude of the required pressure was controlled by a precision hydraulic servo valve with PID control.
The achievable accuracy of the required pressure was 0.002 MPa. The system of hydraulic pressure activation and its control was solved by a separate hydraulic unit working independently of the hydraulic circuit of the given press. From the contactless deformation analysis using the system Mercury RT, the magnitudes of major and minor strain (ε 1 , ε 2 ), as well as the radius of curvature R must be detected to subsequently determine the stress-strain curve. For this purpose, two surfaces have been defined on the measured area, in which the system calculates the desired quantities. An example of the analysis area is shown in Figure 7. The larger analysis area (area A) is used to find the radius of curvature R and the smaller area (area B) is used to calculate the major and minor strain, from which the actual thickness of the material can be calculated using Equations (5) and (6). An example of true strain distribution in the relevant area and evolution of the radius of curvature R is shown in Figure 8.
Similar to the static tensile test, a hardening curve was obtained from the measured data, which was then approximated by the Krupkowski law (1). An example of this approximation of the hardening curve from HBT, so under the equi-biaxial stress state, can be seen in Figure 9. The average values of approximation constants obtained under equi-biaxial loading can be seen in Table 4. Table 4. Stress-strain curve approximation constants from the hydraulic bulge test.  Similar to the static tensile test, a hardening curve was obtained from the measured data, which was then approximated by the Krupkowski law (1). An example of this approximation of the hardening curve from HBT, so under the equi-biaxial stress state, can be seen in Figure 9. The average values of approximation constants obtained under equibiaxial loading can be seen in Table 4.   Similar to the static tensile test, a hardening curve was obtained from the measured data, which was then approximated by the Krupkowski law (1). An example of this approximation of the hardening curve from HBT, so under the equi-biaxial stress state, can be seen in Figure 9. The average values of approximation constants obtained under equibiaxial loading can be seen in Table 4.

Cyclic Test
In addition to the previous tests, we also carried out a cyclic test, which takes the effect of the change of loading direction (+, −), to define the material model termed as Yoshida-Uemori. Due to the compressive stress states, a performance of this test for sheet Figure 9. True stress-strain curve from HBT and application of the Krupkowski law.

Cyclic Test
In addition to the previous tests, we also carried out a cyclic test, which takes the effect of the change of loading direction (+, −), to define the material model termed as Yoshida-Uemori. Due to the compressive stress states, a performance of this test for sheet samples was quite demanding and there was a loss of stability resulting in sample buckling. Thus, a testing jig was proposed as an additional jig for the clamping grips. The testing jig consists of four subdivided supporting grips, which are hydraulically controlled and can prevent the specimen from buckling during the compression. The offset of measured curves from monitored cycles rests in the magnitude of the Bauschinger effect for the tested material. A total of 5 measurements were made under the data scanning frequency of 50 Hz, from which an average curve was determined. The layout of the cyclic test is shown in Figure 10 (left) and results are given in Figure 10 (right).

Cyclic Test
In addition to the previous tests, we also carried out a cyclic test, which takes the effect of the change of loading direction (+, −), to define the material model termed as Yoshida-Uemori. Due to the compressive stress states, a performance of this test for sheet samples was quite demanding and there was a loss of stability resulting in sample buckling. Thus, a testing jig was proposed as an additional jig for the clamping grips. The testing jig consists of four subdivided supporting grips, which are hydraulically controlled and can prevent the specimen from buckling during the compression. The offset of measured curves from monitored cycles rests in the magnitude of the Bauschinger effect for the tested material. A total of 5 measurements were made under the data scanning frequency of 50 Hz, from which an average curve was determined. The layout of the cyclic test is shown in Figure 10 (left) and results are given in Figure 10 (right).

U-Bending of the Specimens
For the experimental determination of the springback, it was suitable to choose such a test when a change of the stress state occurred in the bending area because this case can be fully developed the Bauschinger effect. Concerning the equipment of our laboratory, the U-bending of the workpiece (blank) in the tool with a fixed blank-holder was chosen for the springback analysis. This experiment simulates the process that occurs in a tool at the deep-drawing (stamping) process. During the test, the workpiece is two times bent in the opposite bending direction. The principle of the test is evident from Figure 11, where the scheme of this U-bending test is shown. The test conditions were as follows: distance between fixed blank-holder and die was 0.9 mm (i.e., draw clearance between tool and workpiece was 0 mm). The specimen was in contact with the blank-holder, but there was not force acting at the beginning of the test. The tool is equipped with a lower blank-holder that applies force against the bending punch and prevents the deflection of the sheet. The blank-holder force was developed by a coil spring with an initial force of 4 kN at the start of bending and 6 kN at the end of bending in the bottom dead centre of the bending punch. The die drawing radius was R5, punch drawing radius was R2.5, workpiece dimensions (rectangular shape) 30 mm × 140 mm, bending direction was perpendicular to the rolling direction and own bending was performed on the eccentric press LENP 40.
After finishing the bending test, a specimen was measured by the dimensional and shape analysis on the 3D coordinate measuring device SOMET XZY 464 by the relevant software TANGO1! 3D for own evaluation. As a result of the experimental measurement, there were coordinates of 100 points in *.stp format, which represented the actual shape of the specimen. The principle of measuring the resulting shape is shown in Figure 12.
Software CATIA V5 was used for further processing of the measured points. Through the measured points a curve (spline) was fitted, from which a 30 mm wide area was subsequently extruded. This area had to be further offset by 1.45 mm. This offset value (1.45 mm) was determined by the diameter of the measuring tip (2 mm) and the material thickness (0.9 mm). The offset plane then represented the midplane of the experimentally measured shape of the specimen. The procedure of obtaining the final shape of the experimental specimen in the Catia environment can be seen in Figure 13. The measured shape of the bent blank was subsequently used as a comparative criterion to verify the match between results of the real experiment and numerical simulation (in this case, software PAM-STAMP 2G was used).
for the springback analysis. This experiment simulates the process that occurs in a tool at the deep-drawing (stamping) process. During the test, the workpiece is two times bent in the opposite bending direction. The principle of the test is evident from Figure 11, where the scheme of this U-bending test is shown. The test conditions were as follows: distance between fixed blank-holder and die was 0.9 mm (i.e., draw clearance between tool and workpiece was 0 mm). The specimen was in contact with the blank-holder, but there was not force acting at the beginning of the test. The tool is equipped with a lower blank-holder that applies force against the bending punch and prevents the deflection of the sheet. The blank-holder force was developed by a coil spring with an initial force of 4 kN at the start of bending and 6 kN at the end of bending in the bottom dead centre of the bending punch. The die drawing radius was R5, punch drawing radius was R2.5, workpiece dimensions (rectangular shape) 30 mm × 140 mm, bending direction was perpendicular to the rolling direction and own bending was performed on the eccentric press LENP 40. After finishing the bending test, a specimen was measured by the dimensional and shape analysis on the 3D coordinate measuring device SOMET XZY 464 by the relevant software TANGO1! 3D for own evaluation. As a result of the experimental measurement, there were coordinates of 100 points in *.stp format, which represented the actual shape of the specimen. The principle of measuring the resulting shape is shown in Figure 12. Software CATIA V5 was used for further processing of the measured points. Through the measured points a curve (spline) was fitted, from which a 30 mm wide area was subsequently extruded. This area had to be further offset by 1.45 mm. This offset value (1.45 mm) was determined by the diameter of the measuring tip (2 mm) and the material thickness (0.9 mm). The offset plane then represented the midplane of the experimentally measured shape of the specimen. The procedure of obtaining the final shape of the experimental specimen in the Catia environment can be seen in Figure 13. The measured shape of the bent blank was subsequently used as a comparative criterion to verify the match between results of the real experiment and numerical simulation (in this case, software PAM-STAMP 2G was used). mm) was determined by the diameter of the measuring tip (2 mm) and the material thickness (0.9 mm). The offset plane then represented the midplane of the experimentally measured shape of the specimen. The procedure of obtaining the final shape of the experimental specimen in the Catia environment can be seen in Figure 13. The measured shape of the bent blank was subsequently used as a comparative criterion to verify the match between results of the real experiment and numerical simulation (in this case, software PAM-STAMP 2G was used).

Vegter Yield Criterion
The formulation of the Vegter yield criterion (function) is based on the principle of incorporating and using more experimentally measurable data in the definition of this yield criterion. This criterion is established for the general directions of the principal stresses, which do not need to be identical with the direction of sheet anisotropy. Using the suitable orientation of the coordinate system (principal axes), the shear stress components vanish and the stress state in the given zone is defined by using only the normal stress components (thus principal stresses) σ1 and σ2. The angle θ denotes the relative rotation between the principal stresses and chosen coordinate system. The principle of the Figure 13. Detection of the experimentally prepared specimen shape (CATIA V5).

Vegter Yield Criterion
The formulation of the Vegter yield criterion (function) is based on the principle of incorporating and using more experimentally measurable data in the definition of this yield criterion. This criterion is established for the general directions of the principal stresses, which do not need to be identical with the direction of sheet anisotropy. Using the suitable orientation of the coordinate system (principal axes), the shear stress components vanish and the stress state in the given zone is defined by using only the normal stress components (thus principal stresses) σ 1 and σ 2 . The angle θ denotes the relative rotation between the principal stresses and chosen coordinate system. The principle of the element planar loading and the individual stress components can be subsequently expressed by Equations (7) and (8).
The quadratic Bezier curve is used to mathematically describe the Vegter yield criterion passing through the experimentally determined reference points, as shown in Figure 14. This parametrically determined curve allows the construction of a smooth yield surface with smooth transitions at the reference points (via tangents intersecting in the hinge-points).
Points A and C are experimentally determined reference points defining the yield surface. The tangents to the Bezier curve passing through these reference points are defined by the strain vectors. The intersection of these tangents determines the final shape of the Bezier curve, and it is called a hinge-point [49]. Although the curve does not pass through this hinge-point, its shape strongly depends on this point. Coordinates of any point on the Bezier curve are given by the so-called Bezier parameter µ and position of the hinge-point B (strain vectors). Figure 15 Figure 15). However, it is quite difficult to determine the position of the point PS experimentally and for this reason the value α PS = 0.5 is very often assumed when defining the yield criterion. A similar problem is arising during the experimental determination of the point SH position characterising the pure shear stress. In this case, the presumption about the symmetry of yield criterion and parameter α SH having the value equal to 0.5 is used, as in the case of plane strain point. In the case of experimental determination of these points, no estimation of such extra input parameter α is necessary and the Vegter yield surface directly passes through these experimentally detected points. element planar loading and the individual stress components can be subsequently expressed by Equations (7) and (8).
The quadratic Bezier curve is used to mathematically describe the Vegter yield criterion passing through the experimentally determined reference points, as shown in Figure  14. This parametrically determined curve allows the construction of a smooth yield surface with smooth transitions at the reference points (via tangents intersecting in the hingepoints). Points A and C are experimentally determined reference points defining the yield surface. The tangents to the Bezier curve passing through these reference points are defined by the strain vectors. The intersection of these tangents determines the final shape of the Bezier curve, and it is called a hinge-point [49]. Although the curve does not pass through this hinge-point, its shape strongly depends on this point. Coordinates of any point on the Bezier curve are given by the so-called Bezier parameter μ and position of the hinge-point B (strain vectors). Figure 15 shows part of the curve defining the yield criterion, where coordinates of the point BI denote the yield strength at the equi-biaxial stress state. The position of the point UN denotes the transition between the elastic and plastic state (yield strength) under the uniaxial loading (uniaxial tensile test), and finally point PS defines the yield strength under the plane strain test (so-called plain strain point). Points a and b are the intersections of the tangent line passing through point PS with the tangent line passing through UN and BI, respectively. By introducing an extra input parameter α, which should have a value between 0 and 1, the position of the point PS can be determined (see Figure 15). However, it is quite difficult to determine the position of the point PS experimentally and for this reason the value αPS = 0.5 is very often assumed when  Using the Bezier curve to define the yield criterion, the functional relation can be subsequently expressed by the Equation (9): where: Using the Bezier curve to define the yield criterion, the functional relation can be subsequently expressed by the Equation (9): where: µ = Bezier parameter (0 ≤ µ ≤ 1) A, B, C, are vectors formulating in a two-dimensional stress plane σ 1 and σ 2 the plane stress state condition (σ 3 = 0). The yield criterion is completely dependent on the three parameters mentioned above but is not uniquely determined due to the interdependence of the parabolic functions. The yield surface is described as the individual parts of the Bezier curve affected by the hinge-point B and passing through the reference points A and C. The total yield surface is constructed through the four quadrants expressing the different stress states. Continuity of the first derivative is required between any part of the yield function. The strain vector is derived as follows: where: The independent variables in the above relations are the two principal stresses (σ 1 and σ 2 ) and angle θ between the material axes and directions of the principal stresses.

Kinematic Hardening Law
The Yoshida-Uemori model is based on the principle of two yield surfaces, where the inner surface behaves as a purely kinematic model (i.e., its shape does not change and only the origin of its coordinate system shifts due to the strain hardening). The second surface (outer one-bounding surface) behaves as a purely isotropic model (i.e., the origin of its coordinate system does not change and this surface is increasing due to the strain hardening). The outer yield surface (bounding surface) is thus the limiting factor for the displacement of the inner surface [51].
The chosen model allows us to properly describe all the phenomena occurring during the forming of thin sheets. These are mainly the so-called transient Bauschinger effect during the change of stress state and eventually the sharp yield point of the material. This model also allows expressing the change of Young's modulus depending on the plastic strain. As the most important input for proper setting of the required parameters defining the Yoshida-Uemori model, the true stress-strain curve (hardening curve) measured from the static tensile test and cyclic test as well was used. The change of the mechanical properties due to the Bauschinger effect can be seen in Figure 16.

Definition of the Material Model in the Software PAM-STAMP 2G
Three material models were created in the software PAM-STAMP 2G for the numerical simulation of the bending process. In all cases, the Vegter anisotropic yield criterion based on the measured material characteristics from static tensile tests and HBT were used. Such yield criterion was subsequently combined with different types of hardening laws and the selection of surface or volume elements in the FEA. Finally, three combinations of the mathematical computational model were tested in this paper.

Vegter Yield Criterion in Combination with the Isotropic Hardening Law
Vegter anisotropic yield criterion as defined in the software PAM-STAMP 2G allows taking into account the yield strengths in the individual testing directions (in this case in directions 0°, 45° and 90° with respect to the rolling direction), the yield strength under the equi-biaxial loading (HBT) and plastic strain ratios (Lankford parameters) both from the static tensile test and HBT. In addition to these values, Young's modulus, Poisson's ratio and density of the material to be tested must be entered into the material card. In the case of isotropic hardening law, the stress-strain curve is defined as the average curve from the measured values in the individual tested directions. A weighted average of the relevant quantity (acc. to standard EN ISO 10 113/2020) is calculated from Equation (17). The material card with the Vegter yield criterion in combination with the isotropic hardening law can be seen in Figure 17. For the calculation, a surface element of the deformation mesh, thus "surface blank", was chosen. where: x-weighted average of the monitored quantity x0, x45, x90-measured values in the relevant directions

Definition of the Material Model in the Software PAM-STAMP 2G
Three material models were created in the software PAM-STAMP 2G for the numerical simulation of the bending process. In all cases, the Vegter anisotropic yield criterion based on the measured material characteristics from static tensile tests and HBT were used. Such yield criterion was subsequently combined with different types of hardening laws and the selection of surface or volume elements in the FEA. Finally, three combinations of the mathematical computational model were tested in this paper.

Vegter Yield Criterion in Combination with the Isotropic Hardening Law
Vegter anisotropic yield criterion as defined in the software PAM-STAMP 2G allows taking into account the yield strengths in the individual testing directions (in this case in directions 0 • , 45 • and 90 • with respect to the rolling direction), the yield strength under the equi-biaxial loading (HBT) and plastic strain ratios (Lankford parameters) both from the static tensile test and HBT. In addition to these values, Young's modulus, Poisson's ratio and density of the material to be tested must be entered into the material card. In the case of isotropic hardening law, the stress-strain curve is defined as the average curve from the measured values in the individual tested directions. A weighted average of the relevant quantity (acc. to standard EN ISO 10 113/2020) is calculated from Equation (17). The material card with the Vegter yield criterion in combination with the isotropic hardening law can be seen in Figure 17. For the calculation, a surface element of the deformation mesh, thus "surface blank", was chosen. where: x-weighted average of the monitored quantity

Vegter Yield Criterion in Combination with the Kinematic Hardening Law
To determine the parameters for the model termed as the Yoshida-Uemori hardening law, it is necessary to calculate the parameters defining the kinematic hardening. For this purpose, the specialized software MatPara that has been developed and distributed by the authors of this model, was used. Input data for the calculation of the needed parameters are the data from tensile and cyclic tests. An example of fitting the measured curve to obtain material constants for the Yoshida-Uemori model can be seen in Figure 18. The material card for the Vegter yield criterion and kinematic hardening law is shown in Figure 19. For this material model, a surface element of mesh was chosen again.

Vegter Yield Criterion in Combination with the Kinematic Hardening Law
To determine the parameters for the model termed as the Yoshida-Uemori hardening law, it is necessary to calculate the parameters defining the kinematic hardening. For this purpose, the specialized software MatPara that has been developed and distributed by the authors of this model, was used. Input data for the calculation of the needed parameters are the data from tensile and cyclic tests. An example of fitting the measured curve to obtain material constants for the Yoshida-Uemori model can be seen in Figure 18. The material card for the Vegter yield criterion and kinematic hardening law is shown in Figure 19. For this material model, a surface element of mesh was chosen again.

Vegter Yield Criterion in Combination with the Kinematic Hardening Law
To determine the parameters for the model termed as the Yoshida-Uemori hardening law, it is necessary to calculate the parameters defining the kinematic hardening. For this purpose, the specialized software MatPara that has been developed and distributed by the authors of this model, was used. Input data for the calculation of the needed parameters are the data from tensile and cyclic tests. An example of fitting the measured curve to obtain material constants for the Yoshida-Uemori model can be seen in Figure 18. The material card for the Vegter yield criterion and kinematic hardening law is shown in Figure 19. For this material model, a surface element of mesh was chosen again.

Vegter Yield Criterion in Combination with the Kinematic Hardening Law: Volume Element of the Deformation Mesh
Compared to the previous mathematical model, the volume element of the deformation mesh was chosen for this calculation. In the definition of the bend workpiece, a sheet is entered directly as a sandwich material and three layers are already taken into account for the calculation, as in the real process. Compared to the previous definition, the material definition must take into account the fact that the micro-alloyed steel layer itself has higher mechanical properties than the sandwich as a whole. In the previous cases, the sandwich was taken as a whole and (as is mentioned above) all measured values were related to the total sheet thickness. In this case, every layer was defined to directly correspond to its mechanical properties. In light of the micro-alloyed steel layer, there was only the need to recalculate the approximation constants to a smaller thickness (0.8 mm compared to the original 0.9 mm). The shape of the hardening curve did not change, the curve only shifts towards the higher values, which is reflected by an increase in the approximation constant K in Equation (1). On the other hand, it was quite difficult to determine the mechanical properties of the sandwich core. It was necessary to remove the layers of micro-alloyed steels and to test the mechanical properties of the thin central foil. From the measured values, the hardening curve was again constructed, see Figure 20.
When creating the model of the sheet for bending, the option "Volume Blank-sandwich" was selected in the PAM-STAMP 2G. This option allows (as was mentioned above) us to define the formed sheet by using existing volume elements with the possibility to define any number of layers and also any thickness of these layers. This is quite an advanced approach to specifying a sheet forming operation. However, when using volume elements, it is necessary to follow certain rules when selecting the suitable mesh parameters to keep the stability of the subsequent calculation. Dimensions of the volume element edges should not exceed a ratio of 1:4. Simultaneously, it is necessary to describe individual layers in the thickness direction. In this case, layers had thicknesses of 0.4-0.1-0.4 mm. The size of the element in the sheet thickness direction was thus more or less determined by the thinnest layer in the sandwich material. Furthermore, the minimal number of elements in the sheet thickness direction (which is usually five layers) must be observed. In this case, nine layers in the thickness direction and element having a size of 0.4 mm in the

Vegter Yield Criterion in Combination with the Kinematic Hardening Law: Volume Element of the Deformation Mesh
Compared to the previous mathematical model, the volume element of the deformation mesh was chosen for this calculation. In the definition of the bend workpiece, a sheet is entered directly as a sandwich material and three layers are already taken into account for the calculation, as in the real process. Compared to the previous definition, the material definition must take into account the fact that the micro-alloyed steel layer itself has higher mechanical properties than the sandwich as a whole. In the previous cases, the sandwich was taken as a whole and (as is mentioned above) all measured values were related to the total sheet thickness. In this case, every layer was defined to directly correspond to its mechanical properties. In light of the micro-alloyed steel layer, there was only the need to recalculate the approximation constants to a smaller thickness (0.8 mm compared to the original 0.9 mm). The shape of the hardening curve did not change, the curve only shifts towards the higher values, which is reflected by an increase in the approximation constant K in Equation (1). On the other hand, it was quite difficult to determine the mechanical properties of the sandwich core. It was necessary to remove the layers of micro-alloyed steels and to test the mechanical properties of the thin central foil. From the measured values, the hardening curve was again constructed, see Figure 20.
When creating the model of the sheet for bending, the option "Volume Blank-sandwich" was selected in the PAM-STAMP 2G. This option allows (as was mentioned above) us to define the formed sheet by using existing volume elements with the possibility to define any number of layers and also any thickness of these layers. This is quite an advanced approach to specifying a sheet forming operation. However, when using volume elements, it is necessary to follow certain rules when selecting the suitable mesh parameters to keep the stability of the subsequent calculation. Dimensions of the volume element edges should not exceed a ratio of 1:4. Simultaneously, it is necessary to describe individual layers in the thickness direction. In this case, layers had thicknesses of 0.4-0.1-0.4 mm. The size of the element in the sheet thickness direction was thus more or less determined by the thinnest layer in the sandwich material. Furthermore, the minimal number of elements in the sheet thickness direction (which is usually five layers) must be observed. In this case, nine layers in the thickness direction and element having a size of 0.4 mm in the blank area were chosen so that both the minimal number of elements and the aspect ratio of 1:4 were fulfilled. The own definition of a blank with volume elements can be seen in Figure 21.
Machines 2022, 10, 114 19 of 28 blank area were chosen so that both the minimal number of elements and the aspect ratio of 1:4 were fulfilled. The own definition of a blank with volume elements can be seen in Figure 21.   blank area were chosen so that both the minimal number of elements and the aspect ratio of 1:4 were fulfilled. The own definition of a blank with volume elements can be seen in Figure 21. Figure 20. Stress-strain curve of the sandwich core from the acrylate copolymer. Figure 21. Definition of the blank using volume elements. Figure 21. Definition of the blank using volume elements.

Results from the Finite Element Analysis (FEA)
Regarding the chosen material models, three calculations of the sheet metal forming process by bending technology were performed. In all cases, the same boundary conditions for the calculation were used (rigid tools, constant friction coefficient, constant bending punch feed rate, the position of the bottom dead centre, own course of the lower blankholder force, etc.). For the evaluation of the numerical simulation accuracy, the criterion of matching the geometrical shape between the experimentally obtained shape and the shape obtained by the numerical simulations was chosen. When the springback computation was completed, the experimentally measured surface of the bent specimen was imported into the numerical simulation result and these two surfaces were subsequently compared. One of the post-processing options in the PAM-STAMP 2G is to express the dimensional deviations between two (or more) surfaces. This option was used, and the accuracy of the calculation was assessed based on this criterion. Due to the use of the volume element, the deviations were always related to the midplane of the sheet surfaces. Another possibility from the standard results evaluated in the sheet metal forming rested in the utilisation of the FLD diagram, which can be considered as a criterion for achieving the risk of rupture in the given part. For the chosen type of problem, any limit deformation was not achieved and the calculated deformations on the relevant bent part were well below the forming limit curve (FLC) for the tested sandwich material. The comparison of results between FEA and the real specimen can be seen in Figures 22-24.

Results from the Finite Element Analysis (FEA)
Regarding the chosen material models, three calculations of the sheet metal forming process by bending technology were performed. In all cases, the same boundary conditions for the calculation were used (rigid tools, constant friction coefficient, constant bending punch feed rate, the position of the bottom dead centre, own course of the lower blankholder force, etc.). For the evaluation of the numerical simulation accuracy, the criterion of matching the geometrical shape between the experimentally obtained shape and the shape obtained by the numerical simulations was chosen. When the springback computation was completed, the experimentally measured surface of the bent specimen was imported into the numerical simulation result and these two surfaces were subsequently compared. One of the post-processing options in the PAM-STAMP 2G is to express the dimensional deviations between two (or more) surfaces. This option was used, and the accuracy of the calculation was assessed based on this criterion. Due to the use of the volume element, the deviations were always related to the midplane of the sheet surfaces. Another possibility from the standard results evaluated in the sheet metal forming rested in the utilisation of the FLD diagram, which can be considered as a criterion for achieving the risk of rupture in the given part. For the chosen type of problem, any limit deformation was not achieved and the calculated deformations on the relevant bent part were well below the forming limit curve (FLC) for the tested sandwich material. The comparison of results between FEA and the real specimen can be seen in Figures 22-24.

Discussion
Numerical simulations already represent a standard tool for predicting technological processes in industrial practice. In terms of practically achievable results, where high requirements are set on the accuracy and calculation speed, it is obvious that the proper choice of the mathematical model and the quality of input material data are some of the key factors influencing the final result of FEA. Previous studies have shown that the basic formability calculations can be performed just by using fundamental yield criterions in combination with a simple isotropic hardening law. The definitions of these models can be completed by performing the basic material tests; therefore, just using static tensile testing, where the time required to carry out these tests is relatively short. However, in

Discussion
Numerical simulations already represent a standard tool for predicting technological processes in industrial practice. In terms of practically achievable results, where high requirements are set on the accuracy and calculation speed, it is obvious that the proper choice of the mathematical model and the quality of input material data are some of the key factors influencing the final result of FEA. Previous studies have shown that the basic formability calculations can be performed just by using fundamental yield criterions in combination with a simple isotropic hardening law. The definitions of these models can be completed by performing the basic material tests; therefore, just using static tensile testing, where the time required to carry out these tests is relatively short. However, in

Discussion
Numerical simulations already represent a standard tool for predicting technological processes in industrial practice. In terms of practically achievable results, where high requirements are set on the accuracy and calculation speed, it is obvious that the proper choice of the mathematical model and the quality of input material data are some of the key factors influencing the final result of FEA. Previous studies have shown that the basic formability calculations can be performed just by using fundamental yield criterions in combination with a simple isotropic hardening law. The definitions of these models can be completed by performing the basic material tests; therefore, just using static tensile testing, where the time required to carry out these tests is relatively short. However, in cases where the technological operations involve the localisation of deformation (e.g., hemming) and there is also an increased requirement for predicting the final shape due to the springback effect, it is recommended to use the advanced computational models (e.g., Hill90, Barlat91, Barlat2000, Vegter or Yoshida 6th) in combination with a kinematic hardening law (e.g., Lemaitre and Chaboche or Yoshida-Uemori). The proper definition of these computational models means performing more time-consuming material tests. In addition to the standard static tensile test, there are tests performing under biaxial loading as well as shear, compression and cyclic tests. A frequently discussed issue is the selection of surface or volume elements of the computational mesh. From this point of view, for thin sheets up to a thickness of about 2 mm, it is generally recommended to use a surface element of the mesh. This approach significantly reduces the number of mesh elements, and the accuracy of calculation is considered to be sufficient compared to volume elements, where higher accuracy is assumed.
In this paper, the effect of the chosen calculation strategy on the FEA result was carried out. From the results presented in Figures 22-24, it is evident that the choice of the Vegter yield criterion in combination with the isotropic hardening law provides the worst results in light of the springback prediction. The maximum measured distance between the midplanes of the experimentally measured specimen and results obtained from the numerical simulation is about 2.5 mm (see Figure 22). To assess the deviations of the actual surface (experimentally obtained one) and the FEA surface, a longitudinal section was made along the axis of symmetry through both surfaces. Such a section can be seen in Figure 25, where all tested combinations with the experimentally measured shape (the real one), which is shown in blue colour, are compared. From the mutual comparison of these sections it can be seen that there is a fairly good matching in the bending punch radii (lower parts of the bending punch). However, the results of sections on the bending die drawing edge (upper parts of the specimen) already revealed considerable variations at the use of isotropic hardening law. cases where the technological operations involve the localisation of deformation (e.g., hemming) and there is also an increased requirement for predicting the final shape due to the springback effect, it is recommended to use the advanced computational models (e.g., Hill90, Barlat91, Barlat2000, Vegter or Yoshida 6th) in combination with a kinematic hardening law (e.g., Lemaitre and Chaboche or Yoshida-Uemori). The proper definition of these computational models means performing more time-consuming material tests. In addition to the standard static tensile test, there are tests performing under biaxial loading as well as shear, compression and cyclic tests. A frequently discussed issue is the selection of surface or volume elements of the computational mesh. From this point of view, for thin sheets up to a thickness of about 2 mm, it is generally recommended to use a surface element of the mesh. This approach significantly reduces the number of mesh elements, and the accuracy of calculation is considered to be sufficient compared to volume elements, where higher accuracy is assumed. In this paper, the effect of the chosen calculation strategy on the FEA result was carried out. From the results presented in Figures 22-24, it is evident that the choice of the Vegter yield criterion in combination with the isotropic hardening law provides the worst results in light of the springback prediction. The maximum measured distance between the midplanes of the experimentally measured specimen and results obtained from the numerical simulation is about 2.5 mm (see Figure 22). To assess the deviations of the actual surface (experimentally obtained one) and the FEA surface, a longitudinal section was made along the axis of symmetry through both surfaces. Such a section can be seen in Figure 25, where all tested combinations with the experimentally measured shape (the real one), which is shown in blue colour, are compared. From the mutual comparison of these sections it can be seen that there is a fairly good matching in the bending punch radii (lower parts of the bending punch). However, the results of sections on the bending die drawing edge (upper parts of the specimen) already revealed considerable variations at the use of isotropic hardening law. As in the previous case, a section through the resulting surfaces (see Figure 23) was performed when the Vegter yield criterion in combination with the kinematic law and surface elements was chosen. The final sections can be seen also in Figure 25. These results show quite a good shape matching both in the lower parts of the bending punch and in As in the previous case, a section through the resulting surfaces (see Figure 23) was performed when the Vegter yield criterion in combination with the kinematic law and surface elements was chosen. The final sections can be seen also in Figure 25. These results show quite a good shape matching both in the lower parts of the bending punch and in the upper parts of the specimen at bending on the drawing edge of the bending die. The maximum measured deviation between the midplanes in this case was 1.59 mm. It is quite important that this deviation occurs only in the peripheral parts of the specimen and is largely influenced by the flange length. For the proper assessment of the shape matching, the most important from a practical point of view is the shape matching in the areas of rounding (radii area) and the walls between these radii. From Figures 22-24 can be seen that minimal differences could be found in these zones (order of 0.1 mm).
The last evaluated computational model combined the Vegter yield criterion with kinematic hardening law and volume elements. The FEA result was assessed again based upon Figure 24 and also on the performed section, which can be seen in the green colour again in Figure 25. In addition, in this case, the red colour is the shape of the experimentally obtained specimen and the grey colour is from the FEA. The results again show a very good shape matching both in the lower parts of the bending punch (radii) as well as in the upper parts of the specimen on the drawing edge of the bending die. The maximum measured deviation between the midplanes in this case was 1.47 mm, and in light of the shape, similar conclusions can be made as in the previous case where surface elements were used.
In light of the quantitative comparison of achieved results, histograms of deviations between experimentally and numerically determined surfaces were used. In the following the upper parts of the specimen at bending on the drawing edge of the bending die. The maximum measured deviation between the midplanes in this case was 1.59 mm. It is quite important that this deviation occurs only in the peripheral parts of the specimen and is largely influenced by the flange length. For the proper assessment of the shape matching, the most important from a practical point of view is the shape matching in the areas of rounding (radii area) and the walls between these radii. From Figures 22-24 can be seen that minimal differences could be found in these zones (order of 0.1 mm).
The last evaluated computational model combined the Vegter yield criterion with kinematic hardening law and volume elements. The FEA result was assessed again based upon Figure 24 and also on the performed section, which can be seen in the green colour again in Figure 25. In addition, in this case, the red colour is the shape of the experimentally obtained specimen and the grey colour is from the FEA. The results again show a very good shape matching both in the lower parts of the bending punch (radii) as well as in the upper parts of the specimen on the drawing edge of the bending die. The maximum measured deviation between the midplanes in this case was 1.47 mm, and in light of the shape, similar conclusions can be made as in the previous case where surface elements were used.
In light of the quantitative comparison of achieved results, histograms of deviations between experimentally and numerically determined surfaces were used. In the following    Figure 26 shows the distribution of deviations between real and computed surfaces by using the Vegter yield criterion, isotropic hardening law and surface elements.
The distribution of deviations between real and computed surfaces by using Vegter yield criterion, kinematic hardening law and surface elements is given in Figure 27.
Finally, Figure 28 illustrates almost the same situation as in the previous case; however, volume elements of the deformation mesh are used in this case.      Figure 26 shows the distribution of deviations between real and computed surfaces by using the Vegter yield criterion, isotropic hardening law and surface elements. From the previous figures it is obvious that in light of the deviations between real and computed surfaces, the utilisation of the isotropic hardening law provides by far the worst results from all tested combinations. Only about 50% of such deviations could be found up to 0.5 mm. In addition to that, quite a high frequency of 12.58% was determined for the maximal range of deviations (2-2.5 mm). On the other hand, results from the kinematic hardening laws revealed maximal values up to 1.5 mm. More precisely, the maximal deviation was 1.588 mm for surface elements (see Figure 23), but its proportion was zero due to rounding. Moreover, the vast majority (about 90%) of these values could be found in the first range of deviations (0-0.5 mm). Slightly better results were determined by using the volume elements.
When an isotropic hardening law is chosen, it is evident that this model is not able to describe accurately enough the change in the stress state on the drawing edge, which significantly affects the stress and strain in the formed material and thus also the subsequent springback of the used material. From the results, it can be seen that in the lower parts of the specimen on the bending punch radii, where any change of stress state does not occur, the accuracy is quite sufficient for all chosen material models. In the upper parts on the bending die drawing edge, where a change of the stress state already occurs, the FEA results are greatly influenced by the choice of the calculation model. The kinematic hardening law provides significantly more accurate results than the isotropic hardening law. When comparing the FEA results, which use the same computational model (Vegter yield criterion and kinematic hardening law), but differ in the used meshing strategy (surface element vs. volume element), it can be seen that both models provide very similar results. Both used models show quite a high shape matching with the experimentally measured ones. The model with surface elements slightly underestimates the springback magnitude, while the model with volume elements slightly overestimates the springback magnitude compared to the experiment results. However, these differences are minimal, and the calculation accuracy can be considered satisfactory in both cases. On the other hand, the choice of a volume element in the case of sheets having lower thicknesses significantly increases the total number of elements used in the computational mesh. The model with surface elements had 4200 elements (square elements selected with the edge of 1 mm) compared to 236, 250 volume elements. Such a significant increase in the number of elements leads to a hundred-fold increase in the computational time (minutes for surface elements vs. hours for volume elements).

Conclusions
The major aim of the research presented in this paper was to determine the influence of the computational material model (generally meshing strategy) on the result of the FEA springback prediction for a sandwich material made from a micro-alloyed steel and acrylate copolymer. Different material models describing the deformation behaviour of the tested sandwich material in the plastic state (isotropic hardening law vs. kinematic hardening law) were tested in the experimental part. Furthermore, the effect of the selection of the dimensional type of mesh element (surface element vs. volume element) on the result of the numerical simulation were investigated as well. The following conclusions can be stated from the measured and computed results: The isotropic hardening law cannot be used to correctly predict the springback of sandwich material in cases where the stress state changes during the forming process.

2.
The kinematic hardening law provides a more accurate springback prediction compared to the isotropic hardening model regardless of the surface or volume element selection for the computational mesh. 3.
The choice of the meshing strategy does not have any significant effect on the FEA result when the kinematic hardening law is used. The surface and volume elements give almost exactly comparable results for the springback prediction of the sandwich material.

4.
From a quantitative point of view (using histograms of surfaces' deviations, see , it was confirmed that the kinematic hardening law (regardless of the element type) has significantly higher accuracy in springback prediction than the isotropic hardening law. In addition to that, both kinematic hardening laws (surface and volume type of mesh) have almost 90% of the surfaces' deviations up to 0.5 mm compared to the isotropic hardening law, where only 46% can be found up to 0.5 mm.

5.
The definition of the sandwich material using layers of volume elements in the deformation mesh does not provide a significant improvement of the FEM result. 6.
In the numerical simulation of forming the sandwich material, the measured values of the mechanical quantities can be related to the entire sheet (sandwich) thickness, and it is not necessary to distinguish the different deformation and stress behaviour of the individual layers. 7.
From the calculation accuracy point of view, it does not make any sense to use volume elements of the deformation mesh for the thin sheets. Such an approach leads only to a significant increase in computational time.