A New and Direct R-Value Measurement Method of Sheet Metal Based on Multi-Camera DIC System

: More and more attention has been given in the ﬁeld of mechanical engineering to a material’s R-value, a parameter that characterizes the ability of sheet metal to resist thickness strain. Conventional methods used to determine R-value are based on experiments and an assumption of constant volume. Because the thickness strain cannot be directly measured, the R-value is cur-rently determined using experimentally measured strains in the width, and loading directions in combination with the constant volume assumption, to determine the thickness strain indirectly. This paper provides an alternative method for determining the R-value without any assumptions. This method is based on the use of a multi-camera DIC system to measure strains in three directions simultaneously. Two sets of stereo-vision DIC measurement systems, each comprised of two GigE cameras, are placed on the front and back sides of the sample. Use of the double-sided calibration strategy uniﬁes the world coordinate system of the front and back DIC measurement systems to one coordinate system, allowing for the measurement of thickness strain and direct calculation of R-value. The Random Sample Consensus (RANSAC) algorithm is used to eliminate noise in the thickness strain data, resulting in a more accurate R-value measurement.


Introduction
The plastic strain ratio r is a parameter that indicates the ability of a sheet metal to resist thinning or thickening when subjected to either tensile or compressive forces in the plane of the sheet [1,2]. It has been used as an important parameter to evaluate the formability of automotive sheet metal. It is typically advantageous for the material to experience a minimal reduction in the area when subjected to a force; this means that it is a good drawing material with a high R-value. Therefore, an accurate measurement of R-value is of great significance to the study of the tensile and compressive properties of materials. Stickels [3] predicted the plastic-strain ratio of low-carbon steel sheets using the Young's Modulus, while Ghosh [4] calculated the R-value of Al-Mg-Si alloy sheets based on crystal plasticity models. However, since the crystal structure of polycrystalline materials is often difficult to model accurately, the accuracy of the current theoretical estimates is not high [5,6]. Some scholars [7][8][9] employed laser-ultrasound resonance spectroscopy to measure the texture, thickness, and plastic strain ratio on-line. Although the laser-ultrasonic method has the advantages of being fast and non-destructive, the detection accuracy is greatly affected by the vibration of the tested strip and changes in the environment temperature. In some traditional R-value measurement methods [10,11], the thickness strain is calculated using length strain and width strain. The stated measurement methods are based on the assumption of constant volume, which states that the volume of the material will not change during the deformation process. However, some materials do not always follow

Basic Theory
The stereo vision system shown in Figure 1, or a two-camera system with different perspectives, plays an important role in the DIC measurement system. A space point p is recorded simultaneously by two cameras on the image plane with a different perspective. The points p1 and p2 on the image plane of the two cameras are then matched using the DIC matching algorithm, and the 3D coordinates of point p can be recovered based on the geometry of the stereo-vision system.
The core idea of the DIC method is to measure displacement and strain information on an object's surface by matching and tracking feature points in the natural texture or artificial speckle pattern before and after deformation. The DIC matching process obtains the coordinate on the target image after deformation corresponding to a known coordinate position on the reference image taken before deformation. It is usually necessary to divide the image into multiple grids, or subsets, and the deformation of the subset is used to represent the deformation of the local area of the object. Figure 2 shows this principle. On the reference image, a rectangle with a size of (2M + 1) pixel × (2M + 1) pixel and point p(x 0 , y 0 ) at the center is taken as the reference subset. A correlation function is then used to find a subset on the target image with a similar distribution to that of the reference subset. The subset with the largest correlation coefficient is considered the deformed subset. The sum squared difference correlation criterion (SSD) is commonly used to evaluate the similarity between subsets [20]. Assuming that the intensity distribution of the two facets is represented as F(x, y) and G(x , y ), the SSD function has the following form: subset. The subset with the largest correlation coefficient is considered the deformed subset. The sum squared difference correlation criterion (SSD) is commonly used to evaluate the similarity between subsets [20]. Assuming that the intensity distribution of the two facets is represented as ( , ) and ( ′ , ′ ), the SSD function has the following form:

Double-Sided Calibration
Camera calibration plays an important role in DIC measurement, as its quality directly affects measurement accuracy. The principle of camera calibration is to establish the relationship between the position of the camera image pixel and the scene point position. According to the camera imaging model, the relationship between the image coordinate system and the world coordinate system is as follows [21]: where ̅ = [ , , 1] is the image homogeneous coordinate system, ̅ = [ , , , 1] is the world homogeneous coordinate system, is the scale factor, is the camera internal parameter matrix, and [ ] is the camera external parameter matrix. Both the internal and external parameters of the camera can be calculated by Zhang's calibration method [21].
The double-sided calibration strategy uses a double-sided calibration plate to calibrate the front and back dual-camera DIC systems simultaneously. The schematic of the double-sided calibration strategy is shown in Figure 3. Figure 3 shows two typical two- subset. The subset with the largest correlation coefficient is considered the deformed subset. The sum squared difference correlation criterion (SSD) is commonly used to evaluate the similarity between subsets [20]. Assuming that the intensity distribution of the two facets is represented as ( , ) and ( ′ , ′ ), the SSD function has the following form: (1)

Double-Sided Calibration
Camera calibration plays an important role in DIC measurement, as its quality directly affects measurement accuracy. The principle of camera calibration is to establish the relationship between the position of the camera image pixel and the scene point position. According to the camera imaging model, the relationship between the image coordinate system and the world coordinate system is as follows [21]: where ̅ = [ , , 1] is the image homogeneous coordinate system, ̅ = [ , , , 1] is the world homogeneous coordinate system, is the scale factor, is the camera internal parameter matrix, and [ ] is the camera external parameter matrix. Both the internal and external parameters of the camera can be calculated by Zhang's calibration method [21].
The double-sided calibration strategy uses a double-sided calibration plate to calibrate the front and back dual-camera DIC systems simultaneously. The schematic of the double-sided calibration strategy is shown in Figure 3. Figure 3 shows two typical two-

Double-Sided Calibration
Camera calibration plays an important role in DIC measurement, as its quality directly affects measurement accuracy. The principle of camera calibration is to establish the relationship between the position of the camera image pixel and the scene point position. According to the camera imaging model, the relationship between the image coordinate system and the world coordinate system is as follows [21]: where m = [u, v, 1] T is the image homogeneous coordinate system, M = [X, Y, Z, 1] T is the world homogeneous coordinate system, s is the scale factor, A is the camera internal parameter matrix, and [r t] is the camera external parameter matrix. Both the internal and external parameters of the camera can be calculated by Zhang's calibration method [21]. The double-sided calibration strategy uses a double-sided calibration plate to calibrate the front and back dual-camera DIC systems simultaneously. The schematic of the doublesided calibration strategy is shown in Figure 3. Figure 3 shows two typical two-camera stereo vision subsystems, and these two subsystems are linked through calibration. Figure 4 shows the design of the calibration plate. The role of this calibration is to connect these two subsystems into one single global system. The world coordinate system of the front camera is taken as the reference global coordinate, and the back camera coordinate system is transformed to the reference global system. This transformation can be expressed as follows: camera stereo vision subsystems, and these two subsystems are linked through calibration. Figure 4 shows the design of the calibration plate. The role of this calibration is to connect these two subsystems into one single global system. The world coordinate system of the front camera is taken as the reference global coordinate, and the back camera coordinate system is transformed to the reference global system. This transformation can be expressed as follows: In this equation, is the thickness of the double-sided calibration plate.

Optimized Thickness Strain Measurement
The measurement of the thickness strain is the core part of calculating the R-value. Li built a dual-camera DIC system on the front and back of the test piece to realize the non-destructive measurement of a sample's thickness strain [19]. However, as Li wrote in [19], the strain measurement results contain noise, so it is necessary to remove the noise from the strain data. Since the thickness of the test sample is usually a few millimeters, the noise is amplified during the thickness strain measurement. Figure 5 shows the DP980 camera stereo vision subsystems, and these two subsystems are linked through calibration. Figure 4 shows the design of the calibration plate. The role of this calibration is to connect these two subsystems into one single global system. The world coordinate system of the front camera is taken as the reference global coordinate, and the back camera coordinate system is transformed to the reference global system. This transformation can be expressed as follows: In this equation, is the thickness of the double-sided calibration plate.

Optimized Thickness Strain Measurement
The measurement of the thickness strain is the core part of calculating the R-value. Li built a dual-camera DIC system on the front and back of the test piece to realize the non-destructive measurement of a sample's thickness strain [19]. However, as Li wrote in [19], the strain measurement results contain noise, so it is necessary to remove the noise from the strain data. Since the thickness of the test sample is usually a few millimeters, the noise is amplified during the thickness strain measurement. Figure 5 shows the DP980 In this equation, d is the thickness of the double-sided calibration plate.

Optimized Thickness Strain Measurement
The measurement of the thickness strain is the core part of calculating the R-value. Li built a dual-camera DIC system on the front and back of the test piece to realize the non-destructive measurement of a sample's thickness strain [19]. However, as Li wrote in [19], the strain measurement results contain noise, so it is necessary to remove the noise from the strain data. Since the thickness of the test sample is usually a few millimeters, the noise is amplified during the thickness strain measurement. Figure 5 shows the DP980 thickness strain history, where the x-axis represents the number of photos taken and the yellow circle indicates the noise data. Two smoothing algorithms, namely the Least squares algorithm and Random sample consensus (RANSAC) algorithm, were used to fit the curve and reduce noise. A comparison of the results using these two algorithms follows. thickness strain history, where the x-axis represents the number of photos taken and the yellow circle indicates the noise data. Two smoothing algorithms, namely the Least squares algorithm and Random sample consensus (RANSAC) algorithm, were used to fit the curve and reduce noise. A comparison of the results using these two algorithms follows. Figure 5. DP980 thickness strain fitted by least squares algorithm.

Least Squares Algorithm
The least squares method is the most common method used to solve curve fitting problems. The basic idea of the method is to minimize the sum of squares of the error to find the best function match for the data. In general, we can model the expected value of y as an nth degree polynomial, yielding the general polynomial regression model: where = [ 1 , 2 , ⋯ , ] is a parameter to be determined. The objective function ( , ( , )) is minimized to find the optimal estimated value of the parameter in the function ( , ) for m given sets of data ( , ) ( = 1,2, ⋯ , ).
The result using least squares fitting is shown in Figure 5 as a red curve. The red curve represents the results of quadratic polynomial fitting. Observation of Figure 5 indicates that the fitted curve shifts relative to the expected data. Because of this, the algorithm's fitting results become inaccurate if the data contains a large amount of error or noise. This is the limitation of the least squares algorithm.

Random Sample Consensus (RANSAC) Algorithm
The RANSAC algorithm based on the iterative method was proposed by Fischler [22] to solve the problem of inaccurate solutions when using the least squares model on sample data with a large proportion of outliers. Song used the improved RANSAC algorithm to eliminate errors caused by high-temperature heat wave disturbances in DIC measurement

Least Squares Algorithm
The least squares method is the most common method used to solve curve fitting problems. The basic idea of the method is to minimize the sum of squares of the error to find the best function match for the data. In general, we can model the expected value of y as an nth degree polynomial, yielding the general polynomial regression model: where α = [α 1 , α 2 , · · · , α n ] is a parameter to be determined. The objective function L(y, f (x, α)) is minimized to find the optimal estimated value of the parameter α in the function f (x, α) for m given sets of data (x i , y i ) (i = 1, 2, · · · , m).
The result using least squares fitting is shown in Figure 5 as a red curve. The red curve represents the results of quadratic polynomial fitting. Observation of Figure 5 indicates that the fitted curve shifts relative to the expected data. Because of this, the algorithm's fitting results become inaccurate if the data contains a large amount of error or noise. This is the limitation of the least squares algorithm.

Random Sample Consensus (RANSAC) Algorithm
The RANSAC algorithm based on the iterative method was proposed by Fischler [22] to solve the problem of inaccurate solutions when using the least squares model on sample data with a large proportion of outliers. Song used the improved RANSAC algorithm to eliminate errors caused by high-temperature heat wave disturbances in DIC measurement [23]. This algorithm is able to classify data points as outliers or inliers and fit the mathematical model through inliers while ignoring outliers. The basic algorithm is summarized as follows:

1.
Randomly select several points in the thickness-strain data to calculate the polynomial thickness-strain history model.

2.
Set a preset tolerance value ε and calculate the number of data points that match the mathematical model in Step 1.

3.
If the proportion of data points that meet the mathematical model in Step 2 exceeds the preset threshold τ, recalculate the thickness-strain history model to use these inliers and terminate the algorithm. 4.
Otherwise, repeat Steps 1 to 3 K times.
The number of iterations K should be large enough to ensure that the probability p of at least one set of random samples not including outliers is greater than 0.99. Assuming that the probability that the data is selected as an inlier is a, then b = 1 − a is the probability of observing an outlier. The required number of iterations of the minimum number of points, denoted k, is calculated using the following function: Figure 6 shows that the RANSAC algorithm classifies the thickness strain data into inlier points and outlier points. The blue point set represent the inlier points, and the red circles represent the outlier point set. Using this method, the outlier points, or noisy data, can be ignored during fitting. The result of the RANSAC algorithm fitting is shown in Figure 7. It can be seen that the red fitting curve is fitted with inlier points as the data set. Compared with the least squares algorithm, it is hardly affected by noise in the data. This indicates that the RANSAC algorithm is effective and has important significance for the subsequent R-value calculation.
Metals 2021, 11, x FOR PEER REVIEW 6 of 16 [23]. This algorithm is able to classify data points as outliers or inliers and fit the mathematical model through inliers while ignoring outliers. The basic algorithm is summarized as follows: 1. Randomly select several points in the thickness-strain data to calculate the polynomial thickness-strain history model. 2. Set a preset tolerance value and calculate the number of data points that match the mathematical model in Step 1. 3. If the proportion of data points that meet the mathematical model in Step 2 exceeds the preset threshold , recalculate the thickness-strain history model to use these inliers and terminate the algorithm. 4. Otherwise, repeat Steps 1 to 3 K times.
The number of iterations K should be large enough to ensure that the probability p of at least one set of random samples not including outliers is greater than 0.99. Assuming that the probability that the data is selected as an inlier is a, then b = 1 − a is the probability of observing an outlier. The required number of iterations of the minimum number of points, denoted k, is calculated using the following function: Figure 6 shows that the RANSAC algorithm classifies the thickness strain data into inlier points and outlier points. The blue point set represent the inlier points, and the red circles represent the outlier point set. Using this method, the outlier points, or noisy data, can be ignored during fitting. The result of the RANSAC algorithm fitting is shown in Figure 7. It can be seen that the red fitting curve is fitted with inlier points as the data set. Compared with the least squares algorithm, it is hardly affected by noise in the data. This indicates that the RANSAC algorithm is effective and has important significance for the subsequent R-value calculation.

DP980 R-Value Determination
In this section, the R-value calculation method is reviewed first, followed by a discussion of the experimental setup and process. The composite material DP980 was used as the experimental sample. Finally, the experimental results are analyzed and discussed.

R-Value Calculation
The equation to calculate R-value is shown as Equation (8). For sheet metals, the R-value is usually determined from three different directions of in-plane loading (0°, 45°, 90° in the rolling direction) and the normal R-value is calculated using the following equation:

Experiment System
The experimental system is shown in Figures 8 and 9. Four GigE cameras (maximum resolution 2752 × 2200) were used to form two independent DIC binocular measurement systems. One DIC system was placed in front of the test sample, while the other DIC system was placed behind the test sample. Two pairs of 50 mm fixed focus lenses were used to achieve a 20 mm × 20 mm field of view. The spatial resolution was 0.08 mm/pixel and the acquisition rate was 10 Hz. Two powerful LED lights provided sufficient illumination for the front and back DIC measurement systems. A 50 KN MTS tensile test machine was used to do the tensile test, and a 1 mm thick dog bone shape sample made from DP980 sheet metal was prepared for the test. The tensile speed was set at 3 mm/min. The slow tensile rate used for the test enabled more images to be taken before the sample cracked, thus providing more data to analyze. Figure 10 is a schematic diagram of this sample.

DP980 R-Value Determination
In this section, the R-value calculation method is reviewed first, followed by a discussion of the experimental setup and process. The composite material DP980 was used as the experimental sample. Finally, the experimental results are analyzed and discussed.

R-Value Calculation
The equation to calculate R-value is shown as Equation (8).
where ε 2 , ε 3 represent the true width strain and true thickness strain, respectively. For sheet metals, the R-value is usually determined from three different directions of in-plane loading (0 • , 45 • , 90 • in the rolling direction) and the normal R-value is calculated using the following equation:

Experiment System
The experimental system is shown in Figures 8 and 9. Four GigE cameras (maximum resolution 2752 × 2200) were used to form two independent DIC binocular measurement systems. One DIC system was placed in front of the test sample, while the other DIC system was placed behind the test sample. Two pairs of 50 mm fixed focus lenses were used to achieve a 20 mm × 20 mm field of view. The spatial resolution was 0.08 mm/pixel and the acquisition rate was 10 Hz. Two powerful LED lights provided sufficient illumination for the front and back DIC measurement systems. A 50 KN MTS tensile test machine was used to do the tensile test, and a 1 mm thick dog bone shape sample made from DP980 sheet metal was prepared for the test. The tensile speed was set at 3 mm/min. The slow tensile rate used for the test enabled more images to be taken before the sample cracked, thus providing more data to analyze. Figure 10 is a schematic diagram of this sample. Current limitations with this system include the cost of the cameras used for the DIC systems and the large amount of space required to accommodate the front and rear camera groups. Current limitations with this system include the cost of the cameras used for the DIC systems and the large amount of space required to accommodate the front and rear camera groups.   Current limitations with this system include the cost of the cameras used for the DIC systems and the large amount of space required to accommodate the front and rear camera groups.

Experimental Results Analysis and Discussion
The 3D contour and strain distribution on the surfaces of the DP980 was obtained through DIC measurement. Figure 11 illustrates the simultaneous measurement results of length strain ( 1 ) map, width strain ( 2 ) map and thickness strain ( 3 ) map before fracture

Experimental Results Analysis and Discussion
The 3D contour and strain distribution on the surfaces of the DP980 was obtained through DIC measurement. Figure 11 illustrates the simultaneous measurement results of length strain (ε 1 ) map, width strain (ε 2 ) map and thickness strain (ε 3 ) map before fracture in the tensile test. The necking area, which is indicative of a strain concentration, is clearly visible.
R-value measurements were taken using two different methods, and these results were compared and analyzed. The first method is the traditional measurement, or indirect measurement, method based on the constant volume assumption. According to the assumption of constant volume, ε 1 + ε 2 + ε 3 = 0; thus, the thickness strain can be derived from measurements of the length and width strain: Based on the definition of R-value: The other method used is the direct measurement method, or the direct thickness strain measurement method proposed in this paper. The calculation method is shown in Equation (8). Data in both the non-necking and necking area was selected for analysis to better show the difference between the two measurement methods.

Experimental Results Analysis and Discussion
The 3D contour and strain distribution on the surfaces of the DP980 was obtained through DIC measurement. Figure 11 illustrates the simultaneous measurement results of length strain ( 1 ) map, width strain ( 2 ) map and thickness strain ( 3 ) map before fracture in the tensile test. The necking area, which is indicative of a strain concentration, is clearly visible.
R-value measurements were taken using two different methods, and these results were compared and analyzed. The first method is the traditional measurement, or indirect measurement, method based on the constant volume assumption. According to the assumption of constant volume, 1 + 2 + 3 = 0; thus, the thickness strain can be derived from measurements of the length and width strain: Based on the definition of R-value: The other method used is the direct measurement method, or the direct thickness strain measurement method proposed in this paper. The calculation method is shown in Equation (8). Data in both the non-necking and necking area was selected for analysis to better show the difference between the two measurement methods.

Non-Necking Area
A square area of 8 × 8 pixels, as shown in Figure 12a, was selected instead of a single point in the non-necking area to reduce the influence of noise on the data. Each pixel represents 0.2 mm. The calculation results for the traditional and proposed methods are shown in Figure 13. According to the definition of R-value, the R-value for DP980 is equal to the slope of these two curves. It can be seen from the figure that the relationship between thickness strain and width strain is approximately linear, and the amount of noise in the data is low. Therefore, the least squares principle can be used to calculate the slope of the two curves. The directly measured R-value is 0.8546, while the indirectly measured R-value is 0.8656. Based on previous research of the properties of DP980 material, its R-value generally ranges between 0.7-1.0, depending on the material's composition [24]. The results measured by both methods can be considered reasonable because they both fall within this range, and the difference in the results between the two methods is very small. Therefore, both methods are acceptable. point in the non-necking area to reduce the influence of noise on the data. Each pixel represents 0.2 mm. The calculation results for the traditional and proposed methods are shown in Figure 13. According to the definition of R-value, the R-value for DP980 is equal to the slope of these two curves. It can be seen from the figure that the relationship between thickness strain and width strain is approximately linear, and the amount of noise in the data is low. Therefore, the least squares principle can be used to calculate the slope of the two curves. The directly measured R-value is 0.8546, while the indirectly measured R-value is 0.8656. Based on previous research of the properties of DP980 material, its Rvalue generally ranges between 0.7-1.0, depending on the material's composition [24]. The results measured by both methods can be considered reasonable because they both fall within this range, and the difference in the results between the two methods is very small. Therefore, both methods are acceptable.   resents 0.2 mm. The calculation results for the traditional and proposed methods are shown in Figure 13. According to the definition of R-value, the R-value for DP980 is equal to the slope of these two curves. It can be seen from the figure that the relationship between thickness strain and width strain is approximately linear, and the amount of noise in the data is low. Therefore, the least squares principle can be used to calculate the slope of the two curves. The directly measured R-value is 0.8546, while the indirectly measured R-value is 0.8656. Based on previous research of the properties of DP980 material, its Rvalue generally ranges between 0.7-1.0, depending on the material's composition [24]. The results measured by both methods can be considered reasonable because they both fall within this range, and the difference in the results between the two methods is very small. Therefore, both methods are acceptable.

Necking Area
When necking occurs during a tensile test, a large amount of strain is disproportionately distributed in a central location of the material. As the deformation continues, the strain in the necking area continues to increase until the material ruptures. In the necking stage, the cross-sectional area of the material decreases drastically, while the strain in the thickness and width directions increases sharply. Figure 14 displays an example of the sharp strain increase that can occur during the necking stage.
When necking occurs during a tensile test, a large amount of strain is disproportionately distributed in a central location of the material. As the deformation continues, the strain in the necking area continues to increase until the material ruptures. In the necking stage, the cross-sectional area of the material decreases drastically, while the strain in the thickness and width directions increases sharply. Figure 14 displays an example of the sharp strain increase that can occur during the necking stage. An 8 × 8 pixel square area was selected in the necking area to obtain strain data, as shown in Figure 12b. Figure 15 illustrates the R-value measured by the two methods in the necking area. In the linear elastic deformation stage, the thickness strain is proportional to the width strain. The R-value at this time is a constant, and the results measured by the two methods are 0.8471 and 0.8526, respectively; these results are similar to those measured in the non-necked area. This shows that the material is uniformly distributed, stable in nature, and isotropic. When the width strain exceeds 60 mm/m, the strain value increases rapidly, and the material deformation enters the necking stage. It is worth noting that the direct method measurement curve no longer increases linearly; the growth rate reduces, and the R-value correspondingly becomes smaller. This shows that the ability of the material to resist strain in the thickness direction is reduced. Compared to the curve for the direct measurement method, the curve for the indirect measurement method remains linear in the early necking stage, and the R-value does not decrease significantly until the later necking stages. This result further confirms some researcher opinions that the material does not necessarily follow the constant volume assumption in the necking stage [12]. Therefore, the R-value measured by the direct measurement method is more reliable and accurate in the necking stage. It should be noted that some materials do follow the assumption of constant volume in the necking stage. The R-values for DP980 in the three rolling directions are displayed in Table 1. The R-values in the three rolling directions are very similar, indicating that the material has similar resistance to thickness deformation in different rolling directions. An 8 × 8 pixel square area was selected in the necking area to obtain strain data, as shown in Figure 12b. Figure 15 illustrates the R-value measured by the two methods in the necking area. In the linear elastic deformation stage, the thickness strain is proportional to the width strain. The R-value at this time is a constant, and the results measured by the two methods are 0.8471 and 0.8526, respectively; these results are similar to those measured in the non-necked area. This shows that the material is uniformly distributed, stable in nature, and isotropic. When the width strain exceeds 60 mm/m, the strain value increases rapidly, and the material deformation enters the necking stage. It is worth noting that the direct method measurement curve no longer increases linearly; the growth rate reduces, and the R-value correspondingly becomes smaller. This shows that the ability of the material to resist strain in the thickness direction is reduced. Compared to the curve for the direct measurement method, the curve for the indirect measurement method remains linear in the early necking stage, and the R-value does not decrease significantly until the later necking stages. This result further confirms some researcher opinions that the material does not necessarily follow the constant volume assumption in the necking stage [12]. Therefore, the R-value measured by the direct measurement method is more reliable and accurate in the necking stage. It should be noted that some materials do follow the assumption of constant volume in the necking stage. The R-values for DP980 in the three rolling directions are displayed in Table 1. The R-values in the three rolling directions are very similar, indicating that the material has similar resistance to thickness deformation in different rolling directions.

Additional Verification Tests
Two additional materials were used for uniaxial tensile tests to verify the repeatability of the method proposed in this paper: aluminum alloy 6061 and a polymer material.

Additional Verification Tests
Two additional materials were used for uniaxial tensile tests to verify the repeatabi ity of the method proposed in this paper: aluminum alloy 6061 and a polymer material.

Aluminum Alloy 6061
A 3 mm thick dog-bone shaped aluminum alloy 6061 sample was utilized for a ten sile test, and the test conditions were consistent with those used for DP980. The non-neck ing and necking areas were used to analyze the R-value of aluminum alloy 6061. Figure  16 and 17 show the R-value curves of the non-necking area and the necking area using th direct measurement method and the indirect measurement method, respectively. A thickness strain data have been filtered using the RANSAC algorithm before plotting th R-value curves. It can be seen from the figures that the two measurement curves of th non-necking area are very close. The R-values obtained by fitting the two curves with leas squares are 0.4895 and 0.4925, respectively. However, the data taken in the necking are shows the curves of the two measurement methods are different in the necking stage be fore fracture. The thickness strain of the direct measurement method increases more rap idly than that resulting from the indirect measurement method; thus, the R-value become smaller. The R-value curve change of this aluminum alloy 6061 based on direct measure ment is similar to that of DP980. The constant volume assumption is no longer followe in the necking stage. The possible reason for this is that the crystal structure of the meta material is dislocated during the necking process.  A 3 mm thick dog-bone shaped aluminum alloy 6061 sample was utilized for a tensile test, and the test conditions were consistent with those used for DP980. The nonnecking and necking areas were used to analyze the R-value of aluminum alloy 6061. Figures 16 and 17 show the R-value curves of the non-necking area and the necking area using the direct measurement method and the indirect measurement method, respectively. All thickness strain data have been filtered using the RANSAC algorithm before plotting the R-value curves. It can be seen from the figures that the two measurement curves of the non-necking area are very close. The R-values obtained by fitting the two curves with least squares are 0.4895 and 0.4925, respectively. However, the data taken in the necking area shows the curves of the two measurement methods are different in the necking stage before fracture. The thickness strain of the direct measurement method increases more rapidly than that resulting from the indirect measurement method; thus, the R-value becomes smaller. The R-value curve change of this aluminum alloy 6061 based on direct measurement is similar to that of DP980. The constant volume assumption is no longer followed in the necking stage. The possible reason for this is that the crystal structure of the metal material is dislocated during the necking process.

Polymer
The R-value of the polymer material is measured in this section. A 2.5 mm thick dogbone specimen was used for uniaxial tensile experiments. Since the polymer material is a typical brittle material, the material does not neck during the tensile test. Any rectangular area of 8 × 8 pixels on the material surface can be selected for R-value calculation. Figure  18 shows the R-value curves of the two measurement methods for this material. It can be found that the measurement results of the two methods are linearly increasing, but the growth rate is quite different. The R-value of the direct measurement is 0.3101, while the

Polymer
The R-value of the polymer material is measured in this section. A 2.5 mm thick dogbone specimen was used for uniaxial tensile experiments. Since the polymer material is a typical brittle material, the material does not neck during the tensile test. Any rectangular area of 8 × 8 pixels on the material surface can be selected for R-value calculation. Figure  18 shows the R-value curves of the two measurement methods for this material. It can be found that the measurement results of the two methods are linearly increasing, but the growth rate is quite different. The R-value of the direct measurement is 0.3101, while the

Polymer
The R-value of the polymer material is measured in this section. A 2.5 mm thick dogbone specimen was used for uniaxial tensile experiments. Since the polymer material is a typical brittle material, the material does not neck during the tensile test. Any rectangular area of 8 × 8 pixels on the material surface can be selected for R-value calculation. Figure 18 shows the R-value curves of the two measurement methods for this material. It can be found that the measurement results of the two methods are linearly increasing, but the growth rate is quite different. The R-value of the direct measurement is 0.3101, while the indirect measurement result is 0.2499. This states that the polymer material does not entirely follow the constant volume assumption in the tensile test. One possible explanation is the constant volume assumption is based on the crystal structure of metal materials, while polymers are composed of polymer compounds with complex structures and no longer follow the deformation conditions of metal materials. indirect measurement result is 0.2499. This states that the polymer material does not entirely follow the constant volume assumption in the tensile test. One possible explanation is the constant volume assumption is based on the crystal structure of metal materials, while polymers are composed of polymer compounds with complex structures and no longer follow the deformation conditions of metal materials.

Conclusions
A new R-value measurement method based on the use of two stereo-DIC systems has been presented in this paper. This method provides a new possibility for direct and simultaneous measurement of strains in the length, width, and thickness directions, thus allowing for direct measurement of the R-value. During the measurement process, it was found that DP980 and aluminum alloy 6061 do not follow the constant volume assumption in the necking stage. It was also found that a polymer material does not follow this assumption during the entire tensile test. The test results have certain significance for the research of different materials. Compared with the traditional R-value measurement method, this method has the advantages of direct measurement without assumption. The constant volume assumption has certain limitations, especially in the necking stage of material tensile deformation and polymer tensile deformation. Thus, the R-value obtained by the direct measurement method has more research value.
A main challenge to measuring strains in the third direction is that the strain data in this direction has more noise than other two strains due to the small value of the thickness. A standard DIC algorithm cannot be applied for such a measurement without using a special noise-removing algorithm, such as RANSAC presented in this paper. Some application limitations include the cost of four cameras and the large space requirement of the front and rear camera groups. Nevertheless, this paper provides an alternative method that directly measures R-value. The results for direct R-value measurement, as well as the comparison with the indirect R-value method, show that the method and the algorithm proposed is feasible and effective. This paper has directly and indirectly measured the R-value of different materials, compared the results of the two measurement methods, and shown the accuracy of the direct measurement method. However, these measurement results are based on uniaxial

Conclusions
A new R-value measurement method based on the use of two stereo-DIC systems has been presented in this paper. This method provides a new possibility for direct and simultaneous measurement of strains in the length, width, and thickness directions, thus allowing for direct measurement of the R-value. During the measurement process, it was found that DP980 and aluminum alloy 6061 do not follow the constant volume assumption in the necking stage. It was also found that a polymer material does not follow this assumption during the entire tensile test. The test results have certain significance for the research of different materials. Compared with the traditional R-value measurement method, this method has the advantages of direct measurement without assumption. The constant volume assumption has certain limitations, especially in the necking stage of material tensile deformation and polymer tensile deformation. Thus, the R-value obtained by the direct measurement method has more research value.
A main challenge to measuring strains in the third direction is that the strain data in this direction has more noise than other two strains due to the small value of the thickness. A standard DIC algorithm cannot be applied for such a measurement without using a special noise-removing algorithm, such as RANSAC presented in this paper. Some application limitations include the cost of four cameras and the large space requirement of the front and rear camera groups. Nevertheless, this paper provides an alternative method that directly measures R-value. The results for direct R-value measurement, as well as the comparison with the indirect R-value method, show that the method and the algorithm proposed is feasible and effective. This paper has directly and indirectly measured the R-value of different materials, compared the results of the two measurement methods, and shown the accuracy of the direct measurement method. However, these measurement results are based on uniaxial tensile tests. Thus, it cannot completely reveal the reasons and limitations of the constant volume assumption. Future work may be necessary to study the microstructure of different materials before and after uniaxial tensile tests and determine the scope of application of the constant volume assumption through the analysis of the material microstructure.