Physically-Based Methodology for the Characterization of Wrinkling Limit Curve Validated by Yoshida Tests

: This paper presents a new experimental methodology for determining the formability limits by wrinkling in sheet metal forming using Yoshida buckling tests (YBTs). The YBT has been the most used formability test by the scientiﬁc community to determine the occurrence of wrinkling, on the one hand, due to its simple clamping conﬁguration directly related to the kinematics of tensile tests. On the other hand, the formability test replicates the occurrence of wrinkling in deformation regions under tension, similar to the side wall wrinkling of a deep drawn part. The new physically-based methodology is built upon the direct analysis of the in-plane minor strain and its strain rate for a point located at the center of the YBT specimen. Application of the methodology in AA1050–O aluminium alloy sheets enabled the experimental determination of the onset of wrinkling. Results also showed that in contrast to conventional techniques, where the wrinkling prediction is based in a speciﬁed percentage of the initial thickness of the blank (5–10%), considering a physical wrinkle-triggering factor enables one to identify the precise moment at which wrinkling occurs. The proposed physically-based methodology provides a new level of understanding of the YBT behavior and the onset of wrinkling determination for sheet metal forming with a higher degree of accuracy.


Introduction
The global need to reduce energy consumption is leading manufacturing industry to transform, and one of the strategies considered by the automotive, aerospace, and aeronautical industries is the use of lightweight construction [1].A solution to obtaining lightweight constructions is to use thinner sheet metals, which makes the material more prone to failure by tearing or wrinkling.
Thus, the identification of the failure limits by necking, fracture, and wrinkling is of major importance when designing the sheet metal forming process [2].The determination of the onset of necking and fracture in sheet metal forming has been addressed by the scientific community leading to the definition of general failure limits [3].However, wrinkling has been addressed by the scientific community in a case-to-case approach, and there is no generalized wrinkling limit that can be used independently of the sheet metal forming process [4].
The occurrence of wrinkling is of major relevance in the sheet forming process tools design due to its direct relation with the final geometry of the manufactured part, the processing parameters, and the active tool elements' damage/wear [5].
The experimental determination of the onset of wrinkling has been addressed through two specimens classification: Examples of the first class of tests include the deep drawing of cylindrical [6] or conical cups [7,8], the annular cup test [9], the spinning of a circular sheet blank [10], and in-plane roll-bending of a strip [11].
The second class of experimental wrinkling tests is directly related to the state of stress locally applied on the sheet that induces wrinkling [12].There are two major stress states that originate wrinkling, compression, and tension.The compression state of stress, Figure 1a, is directly related to the wrinkling flange area, and a non-uniform tensile state of stress, Figure 1b, to the side wall wrinkling of a deep drawn part.Additionally, in sheet metal forming of a complex geometry shear state of stress can also originate wrinkling, Figure 1c.

•
Simple geometries indirectly generate a stress state that induces wrinkling.
Examples of the first class of tests include the deep drawing of cylindrical [6] or conical cups [7,8], the annular cup test [9], the spinning of a circular sheet blank [10], and in-plane roll-bending of a strip [11].
The second class of experimental wrinkling tests is directly related to the state of stress locally applied on the sheet that induces wrinkling [12].There are two major stress states that originate wrinkling, compression, and tension.The compression state of stress, Figure 1a, is directly related to the wrinkling flange area, and a non-uniform tensile state of stress, Figure 1b, to the side wall wrinkling of a deep drawn part.Additionally, in sheet metal forming of a complex geometry shear state of stress can also originate wrinkling, Figure 1c.Experimental tests that allow the determination of the onset of wrinkling due to compressive states of stress are the rectangular test specimens proposed by Kasei et al. [13] and Magrinho et al. [14].The determination of the wrinkling initiation owing to a tensile stress state can be performed by the Yoshida Buckling Test (YBT) originally proposed by Yoshida et al. [15].The YBT has been the most used experimental test by the scientific community to determine the onset of wrinkling; Tang et al. [16] summarize the most common variants of the YBT proposed over the years.
In a previous work [14] the authors proposed a methodology to identify the onset of wrinkling from a rectangular specimen in compression.In the present work, the authors considered the original YBT (with an additional clamping area) and propose a new experimental methodology to determine the onset of wrinkling based on the in-plane minor true strain.The commonly used experimental methodology [17], which considers that wrinkling occurs when the wrinkle height reaches a predefined percentage of the initial thickness of the blank, is revisited and compared with the new proposed methodology.The new experimental methodology approach was shown as appropriate to determine the wrinkling limit curve (WLC) for a specific material.

Mechanical Characterization
The work was carried out on aluminium AA1050 sheets with 1.5 mm thickness annealed at 345 °C for 3 h before testing to obtain the O condition.The mechanical characterization of the material was performed by tensile testing conducted according to the standard ASTM E8/E8M-16a [18], where a series of tensile test specimens (standard specimen-sheet type) were cut from the aluminium sheet at an angle of 0, 45, and 90 degrees with respect to the rolling direction.The tensile testing of such specimens allowed the attainment of the mechanical properties present in Table 1, such as the modulus of elasticity, , the yield stress,  , the ultimate tensile stress,  , the elongation at break, , and lastly, the planar anisotropy coefficient, .Experimental tests that allow the determination of the onset of wrinkling due to compressive states of stress are the rectangular test specimens proposed by Kasei et al. [13] and Magrinho et al. [14].The determination of the wrinkling initiation owing to a tensile stress state can be performed by the Yoshida Buckling Test (YBT) originally proposed by Yoshida et al. [15].The YBT has been the most used experimental test by the scientific community to determine the onset of wrinkling; Tang et al. [16] summarize the most common variants of the YBT proposed over the years.
In a previous work [14] the authors proposed a methodology to identify the onset of wrinkling from a rectangular specimen in compression.In the present work, the authors considered the original YBT (with an additional clamping area) and propose a new experimental methodology to determine the onset of wrinkling based on the in-plane minor true strain.The commonly used experimental methodology [17], which considers that wrinkling occurs when the wrinkle height reaches a predefined percentage of the initial thickness of the blank, is revisited and compared with the new proposed methodology.The new experimental methodology approach was shown as appropriate to determine the wrinkling limit curve (WLC) for a specific material.

Mechanical Characterization
The work was carried out on aluminium AA1050 sheets with 1.5 mm thickness annealed at 345 • C for 3 h before testing to obtain the O condition.The mechanical characterization of the material was performed by tensile testing conducted according to the standard ASTM E8/E8M-16a [18], where a series of tensile test specimens (standard specimen-sheet type) were cut from the aluminium sheet at an angle of 0, 45, and 90 degrees with respect to the rolling direction.The tensile testing of such specimens allowed the attainment of the mechanical properties present in Table 1, such as the modulus of elasticity, E, the yield stress, σ y , the ultimate tensile stress, σ UTS , the elongation at break, A, and lastly, the planar anisotropy coefficient, r.
The average stress-strain response obtained from the tensile tests can be approximated by the following Ludwik-Hollomon's equation, The anisotropy coefficient obtained from the tensile tests for sheets aligned at 45 • with the rolling direction is considerably higher than those at 0 • , and 90 • , this phenomenon was identified previously by Sharma et al. [19] in sheets produced by accumulative roll bonding.

Experimental Plan
Table 2 presents a summary of the YBT specimens that were used in the characterization of the WLC of the aluminium sheet and at least ten repetitions were made for each dimension of the three dimensions.The YBTs were performed at room temperature and in quasi-static operating conditions, i.e., with a tensile load monotonously applied with a velocity of 5 mm/min.The specimens were cut out from the supplied sheets considering the rolling direction longitudinally in the specimens.The average stress-strain response obtained from the tensile tests can be approximated by the following Ludwik-Hollomon's equation,  147.38̅ . [MPa] The anisotropy coefficient obtained from the tensile tests for sheets aligned at 45° with the rolling direction is considerably higher than those at 0°, and 90°, this phenomenon was identified previously by Sharma et al. [19] in sheets produced by accumulative roll bonding.

Experimental Plan
Table 2 presents a summary of the YBT specimens that were used in the characterization of the WLC of the aluminium sheet and at least ten repetitions were made for each dimension of the three dimensions.The YBTs were performed at room temperature and in quasi-static operating conditions, i.e., with a tensile load monotonously applied with a velocity of 5 mm/min.The specimens were cut out from the supplied sheets considering the rolling direction longitudinally in the specimens.The YBTs were performed on a universal testing machine, Instron 4507 (INSTRON, Norwood, MA, EUA), due to the clamping configuration with the kinematics of the tensile tests (Figure 2b), the YBT specimen geometry induces a non-uniform tensile state of stress (Figure 1b) that causes compressive stress in the transversal direction (X direction) on the central area of the specimen.This compressive stress triggers the plastic instability in the Z direction of the specimen.The YBTs were performed on a universal testing machine, Instron 4507 (INSTRON, Norwood, MA, EUA), due to the clamping configuration with the kinematics of the tensile tests (Figure 2b), the YBT specimen geometry induces a non-uniform tensile state of stress (Figure 1b) that causes compressive stress in the transversal direction (X direction) on the central area of the specimen.This compressive stress triggers the plastic instability in the Z direction of the specimen.

Displacement and Strain Measurements
Three-dimensional displacement and strain measurement on the sheet surface of the YBT specimens that were needed to determine the onset of instability of the sheet material was performed by a DIC system, the Dantec Dynamics Model Q-400 3D (Figure 2a).
The surface of the specimens was sprayed with a stochastic black speckle pattern on a uniform background previously painted in matte white and the YBT specimens were illuminated with one spotlight.The DIC system is equipped with two cameras with 6

Displacement and Strain Measurements
Three-dimensional displacement and strain measurement on the sheet surface of the YBT specimens that were needed to determine the onset of instability of the sheet material was performed by a DIC system, the Dantec Dynamics Model Q-400 3D (Figure 2a).
The surface of the specimens was sprayed with a stochastic black speckle pattern on a uniform background previously painted in matte white and the YBT specimens were illuminated with one spotlight.The DIC system is equipped with two cameras with 6 megapixels of resolution and 50.2 mm of focal length with an aperture of f/11 to increase the depth of focus.The frequency of image acquisition was set to 10 frames per second and the correlation algorithm was performed with the INSTRA 4D software (version 4.6, Dantec Dynamics, Skovlunde, Denmark).A facet size of 17 pixels with a spacing grid of 13 pixels was considered.

Methodology
This paper proposes two physically-based methodologies to detect the onset of instability/wrinkling and to determine the formability limits by wrinkling in sheet metal forming through experimental YBTs.

Methodology A
Methodology A is based on a direct analysis of Z displacement (perpendicular direction to the sheet surface, Figure 2b) along the YBT specimen diagonal (section AA), as schematically represented in Figure 3a.The Z displacement analysis along section AA allows obtaining the amplitude of the plastic instability wave h that corresponds to the maximum Z displacement of section AA, Figure 3b.The evolution of the plastic instability wave h over time along the YBT test is schematically presented in Figure 3c and a critical amplitude plastic instability wave h crit can be defined, as 10% of the sheet thickness (Du et al. [17]).The instant of time when the amplitude of plastic instability wave ℎ is equal to the critical amplitude of plastic instability wave ℎ corresponds to the instant of time that allows for determining the critical strain pairs ( ,  ) from the strain loading path of a point localized in the central zone of section AA, as schematically represented in Figure 3d.The critical strain pairs for the YBT specimens with different dimensions are then used to create the WLC in the principal strain space.

Methodology B
Methodology B is based on the direct analysis of in-plane minor strain ( ) in the The instant of time when the amplitude of plastic instability wave h is equal to the critical amplitude of plastic instability wave h crit corresponds to the instant of time that allows for determining the critical strain pairs (ε crit 1 , ε crit 2 ) from the strain loading path of a point localized in the central zone of section AA, as schematically represented in Figure 3d.The critical strain pairs for the YBT specimens with different dimensions are then used to create the WLC in the principal strain space.

Methodology B
Methodology B is based on the direct analysis of in-plane minor strain (ε 2 ) in the central zone of the YBT specimen, as schematically represented in Figure 4a as the measuring area.The central zone of the YBT specimen, where the strains are measured by the DIC system, corresponds to the zone with maximum Z displacement of section AA, considering the wrinkling behavior presented schematically in Figure 4b.Preceding the onset of instability, the evolution of the in-plane minor strain (ε 2 ) in the central zone of the YBT specimen depends on the specimen geometry.However, after plastic instability occurs, the evolution of the in-plane minor strain (ε 2 ) changes direction.The evolution of the strain rate (firsttime derivate of the in-plane minor strain) allows for to detect of this phenomenon and the minimum point corresponds to the instant of time when buckling is triggered (critical time t crit ), as schematically represented in Figure 4c.phenomenon and the minimum point corresponds to the instant of time when buckling is triggered (critical time  ), as schematically represented in Figure 4c.
The critical time  corresponds to the instant of time that allows for determining the critical strain pairs ( ,  ) from the strain loading path of a point localized in the central zone of section AA, as represented in Figure 4d.The critical strain pairs for the YBT specimens with different geometries are then used to create the WLC in the principal strain space.
Methodology A is position-dependent because it relies on the analysis of displacements, while methodology B is time-dependent because it is based on the analysis of in-plane minor strain rates.Although the above-mentioned methodologies have different grounds, both are physically-based methodologies to detect the onset of plastic instability and determine the WLC of sheet metal forming processes.The critical time t crit corresponds to the instant of time that allows for determining the critical strain pairs (ε crit 1 , ε crit 2 ) from the strain loading path of a point localized in the central zone of section AA, as represented in Figure 4d.The critical strain pairs for the YBT specimens with different geometries are then used to create the WLC in the principal strain space.
Methodology A is position-dependent because it relies on the analysis of displacements, while methodology B is time-dependent because it is based on the analysis of in-plane minor strain rates.Although the above-mentioned methodologies have different grounds, both are physically-based methodologies to detect the onset of plastic instability and determine the WLC of sheet metal forming processes.

Methodology A
Figure 5a presents the force vs. displacement curves for specimens of the three dimensions considered experimentally, L = 90, 100, and 110 mm, identified as Y90, Y100, and Y110, respectively.Regardless of the dimension of the specimen, the force evolutions display the same tendency.For this reason, and clarity reasons, only one specimen for each dimension was included, however, all specimens registered similar results.Additionally, it should be highlighted that there is no change in the evolution of the force at the moment in which instability occurs.The DIC system was used to obtain the Z displacement, for methodology A, and in some cases, it was challenging to obtain results for the entire section AA of the YBT specimens.For this reason, it was needed to adapt the initial methodology presented before and used by Du et al. [17].Instead of considering all the section AA, only the center of each YBT specimen was considered, and the onset of wrinkling occurs when this particular point reaches a Z displacement equal to 5% of the thickness of the blank.Thus, The DIC system was used to obtain the Z displacement, for methodology A, and in some cases, it was challenging to obtain results for the entire section AA of the YBT specimens.For this reason, it was needed to adapt the initial methodology presented before and used by Du et al. [17].Instead of considering all the section AA, only the center of each YBT specimen was considered, and the onset of wrinkling occurs when this particular point reaches a Z displacement equal to 5% of the thickness of the blank.Thus, ensuring the evolution obtained is independent of the quality of the captured images by the DIC system in all YBT specimens.To sum up, it is an adaptation of the methodology used by Du et al. [17].
The evolution of the displacement in the Z direction (thickness direction) over time is presented in Figure 5b for two specimens of each set of specimen dimensions for clarity reasons, and since all the experimental tests for each specimen dimension revealed similar results.Contrary to the force-displacement curves, the increasing rate of the Z displacement, or the wrinkle height, is dependent on the specimen dimension, and, as expected, the z displacement increase is faster for the specimens with higher L [17,20].Thus, considering methodology A previously presented, the time instant for which the Z displacement in the central zone of the specimen, wrinkle height, reaches 5% of the initial thickness occurs is inversely proportional to the specimen dimension, i.e., the onset of wrinkling occurs earlier for the specimens with larger dimensions.
Figure 6 presents the principal strain space obtained by the DIC system at the center of the specimen for the three dimensions considered experimentally, and for clarity reasons only one specimen for each dimension was included.Despite the different dimensions, the strain path for all the specimens is approximately the same, corresponding to pure shear, β = ε 2 /ε 1 = −1.Even though the strain path remains with a slope of approximately −1 throughout the test, at the end of the test all specimens register a slight shift of the slope, as indicated in the work of Du et al. [21].In addition, the shift in the slope is higher for specimens with a bigger side length.The evolution of the displacement in the Z direction (thickness direction) over time is presented in Figure 5b for two specimens of each set of specimen dimensions for clarity reasons, and since all the experimental tests for each specimen dimension revealed similar results.Contrary to the force-displacement curves, the increasing rate of the Z displacement, or the wrinkle height, is dependent on the specimen dimension, and, as expected, the z displacement increase is faster for the specimens with higher L [17,20].Thus, considering methodology A previously presented, the time instant for which the Z displacement in the central zone of the specimen, wrinkle height, reaches 5% of the initial thickness occurs is inversely proportional to the specimen dimension, i.e., the onset of wrinkling occurs earlier for the specimens with larger dimensions.
Figure 6 presents the principal strain space obtained by the DIC system at the center of the specimen for the three dimensions considered experimentally, and for clarity reasons only one specimen for each dimension was included.Despite the different dimensions, the strain path for all the specimens is approximately the same, corresponding to pure shear,    ⁄ 1.Even though the strain path remains with a slope of approximately −1 throughout the test, at the end of the test all specimens register a slight shift of the slope, as indicated in the work of Du et al. [21].In addition, the shift in the slope is higher for specimens with a bigger side length.Table 3 presents, for each test specimen dimension, the average time instant for which the onset of instability occurs and the corresponding average principal strains obtained from the experimental work for all the specimens tested.It can be noted that the higher the L of the specimen, the earlier the time instant for which wrinkling occurs and the lower the corresponding principal strains.Table 3 presents, for each test specimen dimension, the average time instant for which the onset of instability occurs and the corresponding average principal strains obtained from the experimental work for all the specimens tested.It can be noted that the higher the L of the specimen, the earlier the time instant for which wrinkling occurs and the lower the corresponding principal strains.7 presents the evolutions of the minor principal strain, ε 2 obtained from the DIC and the corresponding strain rate over time, for the three YBT specimen dimensions: Y90, Y100, and Y110.The minor strain evolutions for all the specimens' geometry show a global decrease throughout the test, displaying no evident behavior that would allow identifying the onset of wrinkling.According to methodology B, the onset of wrinkling corresponds to a change in the direction of the strain rate evolution, or when the strain rate no longer decreases and starts to increase.At the instant where the strain rate evolution reaches a minimum, the principal strains are the ones corresponding to the onset of instability.Figure 7b,d,f show the minor strain distribution at the onset of wrinkling along the YBT specimens with a side length of 90, 100, and 110 mm, these results are representative of the results for the repetitions for each condition.It is possible to state that the YBT specimens with the smaller side lengths registered lower compressive values of minor strain in the center of the specimen at the instability instant.Considering all of the above-mentioned figures have the same color scale, Figure 7b has deeper tones of blue and purple at the center when compared to the other two geometries, thus demonstrating a more advanced state of compression at the instability instant than the other two geometries.
Similarly to Figure 6, Figure 8 presents the principal strain space obtained by the DIC system for a point located at the center of the specimen for the three specimens' dimensions considered experimentally, and, once again, for clarity reasons only one specimen for each dimension was included.The strain path for all the specimens is approximately the same, corresponding to pure shear stress state, β = ε 2 /ε 1 = −1.Table 4 presents for each YBT specimen geometry the time instant for which the onset of wrinkling occurs, and the corresponding average principal strains obtained from the experimental work.It can be noted that the higher the L of the specimen, the earlier the time instant for which wrinkling occurs and the lower the corresponding principal strains.The minor strain and minor strain rate behavior are independent of the side to which the wrinkle occurs.

Methodologies Comparison
This section presents a comparison between the results obtained using methodologies A and B to determine the onset of instability for Y100 specimens.In general, for all the specimens analyzed, the instability instant corresponding to the onset of wrinkling obtained from methodology B always occurs before the instability instant obtained from methodology A.
Figure 9 presents the Z displacement along the X position of section AA from Figure 3a, where each curve represents the section Z displacement evolution taken every 4 s from the beginning of the test up to a fully developed wrinkle.The solid black line corresponds to section AA for the instability instant obtained from methodology A and the black dotted Table 4 presents for each YBT specimen geometry the time instant for which the onset of wrinkling occurs, and the corresponding average principal strains obtained from the experimental work.It can be noted that the higher the L of the specimen, the earlier the time instant for which wrinkling occurs and the lower the corresponding principal strains.The minor strain and minor strain rate behavior are independent of the side to which the wrinkle occurs.

Methodologies Comparison
This section presents a comparison between the results obtained using methodologies A and B to determine the onset of instability for Y100 specimens.In general, for all the specimens analyzed, the instability instant corresponding to the onset of wrinkling obtained from methodology B always occurs before the instability instant obtained from methodology A. Figure 9 presents the Z displacement along the X position of section AA from Figure 3a, where each curve represents the section Z displacement evolution taken every 4 s from the beginning of the test up to a fully developed wrinkle.The solid black line corresponds to section AA for the instability instant obtained from methodology A and the black dotted line corresponds to section AA for the instability instant obtained from methodology B. The difference between the critical instant at the onset of wrinkling determined by methodologies A and B is not expressive.To compare both methodologies, the conversion of the Z displacement of the center of the specimen at the critical instant determined by methodology B to a percentage of the sheet metal thickness was performed to allow a direct comparison between the two methodologies.Results for Y90, Y100, and Y110 show a Z displacement of 2.70%, 2.67%, and 2.07%, respectively.These results allow a conclusion that by reducing the dimension of the YBT specimen the difference between the critical instant at the onset of wrinkling determined by methodologies A and B is smaller.
Figure 10 shows the evolution of the Z displacement over time at a point in the center of specimen Y100, and the critical instants determined as the onset of wrinkling obtained by methodologies A and B. Once again, the instability instant corresponding to the onset of wrinkling obtained from methodology B occurs before the instability instant obtained from methodology A.  The difference between the critical instant at the onset of wrinkling determined by methodologies A and B is not expressive.To compare both methodologies, the conversion of the Z displacement of the center of the specimen at the critical instant determined by methodology B to a percentage of the sheet metal thickness was performed to allow a direct comparison between the two methodologies.Results for Y90, Y100, and Y110 show a Z displacement of 2.70%, 2.67%, and 2.07%, respectively.These results allow a conclusion that by reducing the dimension of the YBT specimen the difference between the critical instant at the onset of wrinkling determined by methodologies A and B is smaller.
Figure 10 shows the evolution of the Z displacement over time at a point in the center of specimen Y100, and the critical instants determined as the onset of wrinkling obtained by methodologies A and B. Once again, the instability instant corresponding to the onset of wrinkling obtained from methodology B occurs before the instability instant obtained from methodology A. The difference between the critical instant at the onset of wrinkling determined by methodologies A and B is not expressive.To compare both methodologies, the conversion of the Z displacement of the center of the specimen at the critical instant determined by methodology B to a percentage of the sheet metal thickness was performed to allow a direct comparison between the two methodologies.Results for Y90, Y100, and Y110 show a Z displacement of 2.70%, 2.67%, and 2.07%, respectively.These results allow a conclusion that by reducing the dimension of the YBT specimen the difference between the critical instant at the onset of wrinkling determined by methodologies A and B is smaller.
Figure 10 shows the evolution of the Z displacement over time at a point in the center of specimen Y100, and the critical instants determined as the onset of wrinkling obtained by methodologies A and B. Once again, the instability instant corresponding to the onset of wrinkling obtained from methodology B occurs before the instability instant obtained from methodology A.  The methodologies A and B described in Sections 3.1 and 3.2 allowed the determination of the critical strain pairs (ε crit 1 , ε crit 2 ) at the onset of plastic instability for the YBT specimens listed in Table 2.The experimental plan considered the use of three different initial side lengths, L, of the YBT allowing obtaining the WLC for the AA1050-O aluminium sheet with 1.5 mm thickness.Figure 11 shows the critical strain pairs plotted in the principal strain space, as well as the WLC obtained with the two considered methodologies.
The methodologies A and B described in Sections 3.1 and 3.2 allowed the determination of the critical strain pairs ( ,  ) at the onset of plastic instability for the YBT specimens listed in Table 2.The experimental plan considered the use of three different initial side lengths, , of the YBT allowing obtaining the WLC for the AA1050-O aluminium sheet with 1.5 mm thickness.Figure 11 shows the critical strain pairs plotted in the principal strain space, as well as the WLC obtained with the two considered methodologies.
The WLC obtained by methodology A was built from the linear regression of the grey markers corresponding to the critical strain pairs obtained, and parallel the WLC obtained from methodology B was built from the linear regression of the black markers corresponding to the critical strain pairs obtained.Both WLCs present a slope near '−1′ similar to that of pure shear.As observed in Figures 9 and 10 the onset of wrinkling defined from methodology B occurs earlier than the one registered from methodology A. This fact is visible in the WLCs because the formability limit built from critical strain pairs obtained with methodology A presents larger strain values.
YBT specimens with larger dimensions () present an instability instant earlier and respective strain values lower than specimens with smaller dimensions.Therefore, the YBT must be performed with different geometries (different  ) in order to obtain a formability limit by wrinkling across a wider range of sheet-forming conditions.
Methodology B, originally proposed by the authors, revealed that considering a physical wrinkling trigger allows determining the exact onset of wrinkling, contrary to the conventional methodologies (e.g., Methodology A) for which the onset of wrinkling prediction considers a predefined percentage of the initial thickness of the blank.

Conclusions
The YBT has been the most used formability test by the scientific community to determine the formability limits by the wrinkling of sheet material.Its popularity is due to the simple clamping configuration with the kinematics of tensile tests.
A new experimental methodology to determine the onset of wrinkling based on the in-plane minor true strain is validated by comparison with the commonly used experimental methodology, which considers that wrinkling occurs when the wrinkle height reaches a predefined percentage of the initial thickness of the blank.The WLC obtained by methodology A was built from the linear regression of the grey markers corresponding to the critical strain pairs obtained, and parallel the WLC obtained from methodology B was built from the linear regression of the black markers corresponding to the critical strain pairs obtained.
Both WLCs present a slope near '−1' similar to that of pure shear.As observed in Figures 9 and 10 the onset of wrinkling defined from methodology B occurs earlier than the one registered from methodology A. This fact is visible in the WLCs because the formability limit built from critical strain pairs obtained with methodology A presents larger strain values.
YBT specimens with larger dimensions (L) present an instability instant earlier and respective strain values lower than specimens with smaller dimensions.Therefore, the YBT must be performed with different geometries (different L) in order to obtain a formability limit by wrinkling across a wider range of sheet-forming conditions.
Methodology B, originally proposed by the authors, revealed that considering a physical wrinkling trigger allows determining the exact onset of wrinkling, contrary to the conventional methodologies (e.g., Methodology A) for which the onset of wrinkling prediction considers a predefined percentage of the initial thickness of the blank.

Conclusions
The YBT has been the most used formability test by the scientific community to determine the formability limits by the wrinkling of sheet material.Its popularity is due to the simple clamping configuration with the kinematics of tensile tests.
A new experimental methodology to determine the onset of wrinkling based on the inplane minor true strain is validated by comparison with the commonly used experimental methodology, which considers that wrinkling occurs when the wrinkle height reaches a predefined percentage of the initial thickness of the blank.
The new experimental methodology revealed an appropriate and more accurate way to determine the WLC for sheet material than the conventional methodology.This is because methodology B is a physically-based methodology built upon the direct analysis of

Figure 1 .
Figure 1.States of stress that induce wrinkling in a metal sheet: (a) Compressive stress; (b) Nonuniform tensile stress; (c) Shear stress.

Figure 1 .
Figure 1.States of stress that induce wrinkling in a metal sheet: (a) Compressive stress; (b) Nonuniform tensile stress; (c) Shear stress.

Figure 2 .
Figure 2. Schematic representation of the setup utilized by the DIC system for measuring 3D displacement and strains of the YBT (a) before and (b) after deformation.

Figure 2 .
Figure 2. Schematic representation of the setup utilized by the DIC system for measuring 3D displacement and strains of the YBT (a) before and (b) after deformation.

Figure 3 .
Figure 3. Summary of methodology A to characterize the onset of instability: (a) schematic representation of YBT; (b) identification of wrinkling amplitude h; (c) evolution of the wrinkling amplitude h over time; (d) representation of the strain loading paths and identification of the critical strains in principal strain space.

Figure 3 .
Figure 3. Summary of methodology A to characterize the onset of instability: (a) schematic representation of YBT; (b) identification of wrinkling amplitude h; (c) evolution of the wrinkling amplitude h over time; (d) representation of the strain loading paths and identification of the critical strains in principal strain space.

Figure 4 . 4 .Figure 4 .
Figure 4. Summary of methodology B to characterize the onset of instability: (a) schematic representation of YBT; (b) identification of the displacement and strain measurement area; (c) evolution of the in-plane minor strain and strain rate; (d) representation of the strain loading paths and identification of the critical strains in principal strain space.4.Results and Discussion4.1.Methodology AFigure5apresents the force vs. displacement curves for specimens of the three dimensions considered experimentally,  = 90, 100, and 110 mm, identified as Y90, Y100,

Metals 2023 ,Figure 5 .
Figure 5. YBT specimens with a side length of 90, 100, and 110 mm: (a) force-displacement curves; (b) the z displacement of the center.

Figure 5 .
Figure 5. YBT specimens with a side length of 90, 100, and 110 mm: (a) force-displacement curves; (b) the z displacement of the center.

Figure 6 .
Figure 6.Strain loading path for the YBT specimens with a side length of 90, 100, and 110 mm with the critical strain pair for the onset of instability according to methodology A highlighted.

Figure 6 .
Figure 6.Strain loading path for the YBT specimens with a side length of 90, 100, and 110 mm with the critical strain pair for the onset of instability according to methodology A highlighted.

Figure 7 .
Figure 7. Minor strain and minor strain rate evolution over time for specimens: (a) Y90; (c) Y100; (e) Y110.Minor strain distribution at the onset of wrinkling obtained from the DIC according to methodology B for specimens: (b) Y90; (d) Y100; (f) Y110.

Figure 7 .
Figure 7. Minor strain and minor strain rate evolution over time for specimens: (a) Y90; (c) Y100; (e) Y110.Minor strain distribution at the onset of wrinkling obtained from the DIC according to methodology B for specimens: (b) Y90; (d) Y100; (f) Y110.

Metals 2023 , 15 Figure 8 .
Figure 8. Representative strain loading path for the modified Yoshida specimens with a side length of 90, 100, and 110 mm, with the onset of wrinkling according to methodology B highlighted.

Figure 8 .
Figure 8. Representative strain loading path for the modified Yoshida specimens with a side length of 90, 100, and 110 mm, with the onset of wrinkling according to methodology B highlighted.

Metals 2023 , 15 Figure 9 .
Figure 9. Z displacement of section AA of a Y100 specimen before wrinkling, at the onset of wrinkling according to methodologies A and B, and for a fully developed wrinkle.

Figure 10 .
Figure 10.Wrinkle height for a YBT specimen with a side length of 100 mm when the onset of wrinkling occurs for methodologies A and B.

Figure 9 .
Figure 9. Z displacement of section AA of a Y100 specimen before wrinkling, at the onset of wrinkling according to methodologies A and B, and for a fully developed wrinkle.

Metals 2023 , 15 Figure 9 .
Figure 9. Z displacement of section AA of a Y100 specimen before wrinkling, at the onset of wrinkling according to methodologies A and B, and for a fully developed wrinkle.

Figure 10 .
Figure 10.Wrinkle height for a YBT specimen with a side length of 100 mm when the onset of wrinkling occurs for methodologies A and B.

Figure 10 .
Figure 10.Wrinkle height for a YBT specimen with a side length of 100 mm when the onset of wrinkling occurs for methodologies A and B.

Figure 11 .
Figure 11.Formability limits by wrinkling (WLC) according to methodologies A and B of the AA1050-O aluminum with 1.5 mm thickness.

Figure 11 .
Figure 11.Formability limits by wrinkling (WLC) according to methodologies A and B of the AA1050-O aluminum with 1.5 mm thickness.

Table 1 .
Summary of the mechanical properties of aluminum AA1050-O sheet with 1.5 mm thickness.

Table 2 .
Summary of the experimental work plan for the Yoshida buckling tests.

Table 2 .
Summary of the experimental work plan for the Yoshida buckling tests.

Table 3 .
Average time instant at the onset of wrinkling for the three specimen dimensions determined by Methodology A.

Table 3 .
Average time instant at the onset of wrinkling for the three specimen dimensions determined by Methodology A.

Table 4 .
Average time instant at the onset of wrinkling for the three specimen dimensions determined by Methodology B.

Table 4 .
Average time instant at the onset of wrinkling for the three specimen dimensions determined by Methodology B.