Abstract

This paper presents experimental results from shaking table tests on two reduced-scale geogrid reinforced soil-retaining walls (RSRWs) constructed using standard soil, modular facing blocks, and uniaxial geogrid reinforcement to investigate the distribution of the geogrid strain and the mode of potential failure surface for dynamic loading conditions. Similitude relationships for shaking table tests in a 1 g gravitational field were used to scale the specimen geometry, applied characteristics of the earthquake motions. The lateral displacement of the top model is sufficiently large for the top-model block to fall down, and the RSRW is thus destroyed. The tensile strain at the lower part is greater than that at the upper part of the RSRW. The tensile strain in different layers for two-tiered RSRW is consistent with single-step RSRW. On comparing the measured maximum tensile strain lines of the geogrid with the result of the existing calculation method of the potential failure surface, it can be observed that the existing partial calculation method is conservative. Based on the calculation methods of various potential failure surfaces and the measured data, the use of a two-tiered fold-line failure surface is proposed for the two-tiered RSRW while taking into consideration the width of the platform. And it is advised that the failure surface calculation method of BS8006 be used as the calculation method for the potential failure surface of the single-step RSRW under dynamic motion.

1. Introduction

Reinforcement is an important part of a reinforced soil-retaining wall (RSRW). To ensure the overall stability of the RSRW, the reinforcement is required to have sufficient strength and sufficient friction with the backfill. The reinforcement stress condition is the key to the internal stability analysis of the RSRW. The damage to the RSRW is caused by the destruction of the reinforcement; for instance, a RSRW collapsed owing to reinforcement corrosion in the Chile 8.2 earthquake in April of 2014 [1].

The maximum point of the internal stress of each layer of the RSRW is connected to obtain the potential failure surface of the backfill Guler and Selek [2]. The shape and location of the potential failure surface directly affects the design length of the reinforcement and the internal stability analysis of the RSRW. The existing hypothesis regarding the shape and location of the potential failure surface of the RSRW comprises four categories, which are log spiral failure surface, linear failure surface, simplified failure surface of 0.3H, and broken line failure surface. These are presented in Figure 1.

The log spiral failure surface is a curved surface. Based on a centrifuge test, Lei [3] considered that the potential failure surface is the log spiral (Lei method) starting at the foot of the wall or the platform and reaching the rear of the top wall panel at approximately 0.169H.

The linear failure surface is the Rankine failure surface and comprises a self-defined fitting curve. The Rankine failure surface is a plane that passes through the wall foot and the horizontal plane at an angle of (45° + φ/2). The Chinese Technical Code for Application of Geosynthetics (Geosynthetics) [4] and US Department of Transportation Federal Highway Administration (FHWA) [5] use the Rankine failure surface to calculate the flexible reinforcement of RSRWs. Wang and Xu [6] believe that the failure surface of a two-tiered RSRW can be approximated as a straight line (H = −2.1452L + 14.766, Wang fitting curve) with an angle of (45° + φ/2) to the horizontal plane at the wall foot.

The simplified failure surface of 0.3H is used for the convenience of engineering applications. Chinese Specifications for Design of Highway Subgrades (Highway) [7] simplify the failure surface as the lower part as the wall foot (45° + φ/2 to the horizontal plane), the upper part as parallel to the wall, and the upper apex as (0.3H, H); the Chinese Code for Design on Retaining Structures of Railway Subgrades (Railway) [8] and FHWA method [5] (regarding rigid reinforcement) considering that the failure surface is simplified such that it is connected to the coordinate point (0.3H, 0.5H) through the wall and the upper part is parallel to the wall.

The broken line failure surface is a self-defined description of the test results and specifications presented by researchers. Huang et al. [9] performed a prototype observation and numerical simulation of the gabion RSRW. They inferred that the lower part of the failure surface passes through the wall foot and the position of the breaking point (0.2H, H/3) is connected to the upper point (0.3H, H): this failure surface was named as “the new broken-line-type failure surface.” Zhou et al. [10] proposed a simplified potential failure surface for an RSRW with different platform widths based on a finite element simulation; if the upper foot of the upper wall is located in the active zone of the lower wall, in the simplified potential failure surface, the lower part of the wall is the lower wall foot D (platform width) (45° + φ/2 with the horizontal plane). The upper part is D + 0.3 h from the wall (h is the height of the upper wall), upper part failure surface is parallel to the wall (Zhou surface). BS8006 [11] stipulates that the failure surface is the lower part passing through the wall foot, and the position of the broken line is (0.2H, 0.4H), which is connected with the upper point (0.3H, H).

The earthquake force on the RSRW is shared by the geogrid, the blocks, and the backfill (Cai et al. [12]). The acceleration has amplification effect along with the height inside the backfill (Li et al. [13]). The dynamic soil pressure increment is generated under the earthquake motion. The blocks are laterally displaced under the action of the earthquake. The geogrid generates strain due to the restraint of the soil and the blocks. Dynamic soil pressure and geogrid strain cause deformation, which is the external expressive from of RSRW. The limit equilibrium method is one of the applicable methods for calculating the lateral deformation of RSRW (Li et al. [14]). However, there is no unified understanding of the assumption of the failure surface. Few studies have been conducted on the internal stress and potential failure surface of reinforcement under dynamic motion. Thus, the stress and potential failure surface of the single-step RSRW and two-tiered RSRW were studied using the large-scale shaking table test. Based on the comparison of the experimental data with the results of the existing method of potential failure surface calculation, the forms of the potential failure surface suitable for the single-step RSRW and two-tiered RSRW under dynamic action are proposed. The result provides reference for the design and calculation of the geogrid RSRW.

2. Materials and Methods

2.1. Test Equipment

The walls were constructed in the Bidirectional Electrohydraulic Servo Seismic Simulation Shaking Table at the Civil Engineering Test Center of the Institute of Disaster Prevention. The main technical parameters of the shaking table are as follows: the table size is 3.0 m × 3.0 m; the bidirectional lateral seismic simulation is used; the maximum displacement is ±100 mm in the X direction and ±100 mm in the Y direction; the maximum acceleration is 2 g (full load) in the X direction and 2 g (full load) in the Y direction; and the maximum bearing is 20 t.

The acquisition system is mainly of two types: one is a domestic 128-channel dynamic acquisition system that can collect reinforcement strain data, the other is a domestic 16-channel acceleration acquisition system that can collect table acceleration data and soil acceleration data. The model box used for the test is fabricated using steel. The size of the model box: 3.0 m (length) × 2.0 m (height) × 1.5 m (width), as shown in Figure 2.

2.2. Similitude Relationship

To reflect the engineering characteristics of the actual project as realistically as possible, and based on the bearing capacity of the shaking table and the size of the model box, the scaling factor of two-tiered RSRW is 1 : 10. Furthermore, the scaling factor of the single-step RSRW is 1 : 4 and 1 : 2. However, scaling rules (Iai [15], Zheng et al. [16]) generally cannot satisfy the similarity rules for all the parameters. Therefore, it is necessary to pay attention to the main factors affecting the test results, which neglecting the secondary factors. The main scaling factors of the reduced-scale model are derived from the rules proposed by Iai [15], as listed in Table 1.

2.3. Soil

The test uses standard soil as the backfill, the relative density of which is 70%. The grain composition curve of the soil is shown in Figure 3, and the physical and mechanical parameters are listed in Table 2. To strictly control the relative density, the backfill is layered and filled, and its filling method is consistent with the method used by Wang et al. [17].

2.4. Facing Wall

The production of the modular block is mainly based on the prototype of the two-tiered RSRW. The modular-block size of the prototype project is 0.60 m (length) × 0.20 m (width) × 0.25 m (height) and is by the scaling rules of 1 : 10. The modular-block size should be 0.060 m (length) × 0.020 m (width) × 0.025 m (height). Concerning the Guler and Selek [2] test results, the model scale of the block did not affect maximum acceleration measured on the wall. Therefore, the size of the two modular blocks used in the actual project is selected (Wang [18]): 0.25 m (length) × 0.15 m (width) × 0.15 m (height) and 0.125 m (length) × 0.15 m (width) × 0.15 m (height), as shown in Figure 4. The connection between the geogrid and modular block is shown in Figure 5.

2.5. Reinforcement

A uniaxial high-density polyethene (HDPE) geogrid (Xuyu EG50#, Qingdao, China) was used as soil reinforcement. The length of the stretching unit is 22.5 cm. The rib spacing width is 2.22 cm. The rib thickness is 0.1 cm. Tensile tests were conducted on multi-rib specimens at a strain rate of 10%/min, according to ASTM D6637 (ASTM [19]). Results indicate that the geogrid has a tensile strength at 2% strain T2% = 17.4 kN/m and ultimate strength Tult = 50 kN/m in the machine direction. The same reinforcement was used in two reduced-scale models, and the stiffness of the geogrid in the prototype would be 100, 16, and 4 times in the scale model for 1/10, 1/4, and 1/2, respectively. Reinforcement was placed with each soil lift to give a vertical spacing  = 0.15 m. The reinforcement length of the two-tiered RSRW is 0.9 m, which length selected according to the prototype project. At the same time, the reinforcement length of the single-step RSRW is 1.26 m, which length selected from other literature (Ling et al. [20]; Guler and Selek, [2]). The model design of the RSRW, refer to the chosen measured-point-position of strain gauges, is shown in Figure 6.

2.6. Instrumentation

The model design diagrams of specimens are shown in Figure 6. Specimen 1 shows the model of the two-tiered RSRW, wherein the total height is 1.8 m, and the platform width is 0.18 m. A total of 104 strain gauges are placed to observe the dynamic strain of the geogrid. There are four layers on the East side, which are F1, F5, F6, and F7. Each layer is arranged with four positions, and each position is arranged on the upper side. A total of 11 layers are arranged on the West side, each layer is arranged with four positions, and each position has a strain gauge on each of the top and bottom, to correct for bending Bathurst et al. [21].

Specimen 2 shows the model of the single-step RSRW; with a total height of 1.8 m. A total of 88 strain gauges are arranged. To obtain more accurate measures data, the strain gauges are mounted on the top and bottom in the F1, F2, and F3 geogrids, and the top and bottom of the four measured points away from the panel in the F5, F6, and F7 layers.

2.7. Input Motions

To further understand the failure mechanism of the modular block RSRW and identify the failure phenomenon, and obtain the critical strains, the input acceleration is gradually increased until the retaining wall is damaged. Wolong wave (WL) was recorded in Wolong station during the Wenchuan earthquake in China in 2008, and the EL-Centro motion (EL) was recorded in EL-Centro station during the Imperial Valley earthquake in the USA in 1940. The peak acceleration of the WL is 1.0 g, and the time interval of the data points is 0.005 s. The peak acceleration of the EL is 1.0 g. and the time interval of the data points is 0.02 s. The duration time is 58.5 s. The time and peak acceleration may be adjusted and compressed. The acceleration time history and Fourier spectra are presented in Figure 7.

The tests are performed in the order of the peak acceleration from small to large, and the time compression ratio of specimen 2 is from high to low. Each time white noise is applied in between each seismic motion during specimen 1. White noise is input before and after each change in motion magnitude for specimen 2. The loading cases are listed in Tables 3 and 4.

After each loading case, data viewing, recording of the damage outside the specimens, photographing, and other related works are performed. All the acquisition channels are then zeroed, and the model is treated as a new model for the next case until the model is destroyed.

3. Results and Discussion

3.1. Model Damage Phenomena

As the damage mode of the specimens 1 and 2, the top model blocks falls, and the model is destroyed, as shown in Figure 8. The specimen 1 is destroyed in Case 9, as shown in Figure 8(a). The specimen 2 is destroyed in Case 32, as shown in Figure 8(b). It is a very well known facing stability issue; the falling of the top block has nothing to do with the stability of the wall. Therefore, the top two or three blocks are always connected rigidly both in practice and in research to avoid this type of top block failing Guler and Selek [2].

3.2. Geogrid Strain of Specimen 1

Figure 9 shows the strain time history of the F6 layer Y21 under Case 3. Part of the strain time history is selected and divided into two parts for analysis, which includes the horizontal same-layer comparison and vertical same position comparison.

3.2.1. Horizontal Same-Layer Comparison

Figure 10 shows the strain curve of longitudinal section L1, L2, L3, and L4 under different peak accelerations.

In general, the strain values of the strain points in each layer increase with the increase in the peak acceleration motion. When the peak acceleration is small, the strain values of different points are almost the same, and the strain trend approximates a straight line. The strain growth of each strain point is different, and the strain trend gradually changes into a curve with the acceleration increase. There is a difference under WL1.6 g that the strain values under WL 1.6 g are smaller than those under WL 1.2 g. The reason is the top blocks fall under WL 1.6 g, resulting in a reduction in the tensile force of the reinforcement. The strain trends of the L1 and L2 layers of the lower wall are one maximum (minimum) value and are consistent with the strain trend of the L3 and L4 layers of the upper wall. During the increase in the peak acceleration, the maximum strain of the L1 layer gradually approaches Y3, the maximum strain of the L3 layer gradually moves closer to Y28, and the maximum strain of the L4 layer gradually approaches Y44. It can be inferred that the maximum strain point of each layer expands toward the interior of the soil as the acceleration increases.

3.2.2. Vertical Same Position Comparison

Figure 11 shows the distribution of the geogrid strain values of the lower retaining wall transverse section T1 along with the height of the wall. As can be observed from the figure, the strain value increases with the acceleration increase, and the maximum value is observed at F2. Figure 12 shows the distribution of the geogrid strain values of the upper retaining wall T2 and T3 along with the height of the wall. It can be observed from Figure 12 that the strain value increases with the increase in the acceleration. On comparing the Figures 12(a) and 12(b), the distance from the retaining wall is different, and the distributions of the strain are different. The strain at T2 has maximum values at F7 (105 cm) and F9 (135 cm), and a minimum value at F8 (120 cm); the strain at T3 has maximum values at F7 (105 cm) and F10 (150 cm), and a minimum value at F9 (135 cm).

3.3. Geogrid Strain of Specimen 2

Figure 13 shows the strain time history of the Y53 under cases 13 and 14. Part of the strain time history is selected and divided into two parts for analysis, and these comprise the horizontal same-layer comparison and vertical same position comparison.

3.3.1. Horizontal Same-Layer Comparison

Figures 14 and 15 show the strain curve of longitudinal sections L5, L6, L7, and L8 under the WL and EL (model scale: 1/4).

As can be observed from Figures 14 and 15, the strain at each strain point of each layer increases as the acceleration increases. The strain values of each layer are different, but the trend of the same-layer strain is consistent under different peak accelerations. When the peak acceleration is small, the strain values of the strain points are almost the same, and the strain trend line approximates a straight line. The higher the geogrid layer in the retaining wall, the less the trend line of the strain will approximate a straight line. The strain growth of each strain point is different, and the trend line gradually changes into an oblique line or curve with the increase in the peak acceleration. The strain values of the various strain points in each layer are also different under different cases.

On comparing the reinforcement strain of different layers, it can be observed that the reinforcement strain value of the lower wall is greater than that of the upper wall. The strain trend is consistent with the trend of the specimen 1. This condition may be caused by the fact that the reinforcement of the different layers bears different vertical pressures.

3.3.2. Vertical Same Position Comparison

Figure 16 shows the distribution of the strain values at transverse section T4 and T5 under the WL and the EL (model scale: 1/4); Figure 17 shows the distribution of the strain values at T6 and T7 under the WL and EL (model scale: 1/4).

It can be observed from Figures 16 and 17, the strains values at the same position of different layers increase with the acceleration increase. When the similarity ratio and seismic motion are consistent, the distributions of the strain values at the same position (T4, T5, T6, and T7) are nearly identical. The strain law at T4 under WL 0.1 g—WL 1.6 g (model scale: 1/4) that there are two maximum points (F3 (45 cm) and F7 (105 cm)) and one minimum point (F5 (75 cm)). Moreover, the strain laws are nearly the same under different working cases (comparison of T4 with T5, and T6 with T7). The strain laws at T4 and T5 both have a maximum point (F3 (45 cm)) and a minimum point (F5 (75 cm)) under the EL (model scale: 1/4).

When the model scale is consistent, and the seismic waves are different, the strain laws at the same position are slightly different. When the model scale is 1 : 4, the strain trend has two maximum points (F3 (45 cm) and F7 (105 cm)) and one minimum point (F5 (75 cm)) under the WL, while the strain trend has just one maximum point (F3 (45 cm)) and one minimum point (F5 (75 cm)) under the EL. This may be caused by different types of seismic motions.

3.4. The Curve of the Potential Failure Surface

On taking the strain curve of each layer of the specimen 2 as an example, the maximum strain point of each layer of the specimen 2 under cases 31 and 32 are connected in a curve, and the curve is compared with the 0.3H simplified failure surface of the Railway Code [8]. It can be observed that the 0.3H simplified failure surface is conservative, as shown in Figure 18.

3.5. Potential Failure Surface Calculation Method for Specimen 1

The index fitting of the measured strain curve of the specimen 1 is shown in Figure 19. It can be observed that the position of the potential failure surface gradually expands toward the interior of the soil as the acceleration increases.

As for the design method of the multitiered RSRW, only the FHWA [5] and the National Concrete Masonry Association (NCMA) [22] guideline of the U.S are based on different assumptions, and only the design method of the two-tiered RSRW is simply specified. FHWA’s [5] design method is as follows: (1) For a small upper wall offset; D ≤ 1/20 (H1 + H2) (the definition of the wall heights: H1 is the upper wall height, H2 is the lower wall height, and D is the offset distance), it is assumed that the failure surface does not fundamentally change, and it is simply adjusted laterally using the offset distance D. The walls should be designed as a single wall of height H. (2) When D is greater than H2 tan (90° − φ), the walls are not considered to be superimposed and are independently designed from an internal stability viewpoint. (3) Only when (H1 + H2)/20 ≤ D ≤ H2 tan (90° − φ), the internal and external stability calculation of the lower wall are required to be performed by considering the influence of the upper wall load on the vertical pressure of the lower wall. However, the design method of the multitiered RSRW is not explained. The NCMA method stated that the interaction of the upper and lower adjacent walls is not considered in the design of the upper wall, and the upper wall should be equivalent to the top outside uniformly external load of the retaining wall when designing the lower wall.

According to the FHWA calculation method, the offset width of the specimen 1 is 0.18 m, which is in the stage of (H1 + H2)/20 ≤ D ≤ H2 tan (90° − φ). The various methods of calculating the potential failure surface are summarized in Figure 20. Although some methods (the FHWA method, Wang fitting curve, and Zhou simplify potential failure surface) can be used as the simplify potential failure surface at (H1 + H2)/20 ≤ D ≤ H2 tan (90° − φ) stage of specimen 1. Considering the engineering cost, the use of a two-tiered fold-line potential failure surface is proposed for the two-tiered RSRW at the (H1 + H2)/20 ≤ D ≤ H2 tan (90° − φ) stage.

For the two-tiered fold-line potential failure surface, it is assumed that the lower part passes through the wall foot, and the wall foot is used as the origin to establish the coordinate axis of the RSRW. The wall height is the independent variable y, and the failure surface position is the dependent variable x. The failure surface formula is as follows:where D is the platform width; d is the thickness of the wall, hl is the height of lower wall, and hu is the height of the upper wall.

The line type 9 is added in Figure 21; it can better react the failure surface of the two-tiered RSRW at (H1 + H2)/20 ≤ D ≤ H2 tan (90° − φ) stage.

3.6. Potential Failure Surface Calculation Method for Specimen 2

The various calculation methods for the potential failure surface are summarized in Figure 20. It can be observed that although the majority of the maximum measured strain points measured of the reinforcement are on both sides of the log spiral proposed by Lei [3], there are also a few points on the right side of the line. To ensure safety and stability and reduce engineering costs, line 5 (BS8006 surface) is more suitable for obtaining the failure surface formed by the measured maximum strain point of each layer. Therefore, it is suggested that the potential failure surface of the specimen 2 under dynamic motion be selected by BS8006 [11].

4. Discussions

Considering the suddenness and uncertainty of the earthquake, it is difficult to understand the behaviour of geogrid RSRW; the shaking table test is performed to study the actual prototype RSRW through the reduced-scale model. However, the test methods and scaling factors need to be considered. In this study, in measuring the location of the failure surface in model tests, other methods, such as inserting coloured sand layers in advance and carefully observing them after the test, would provide more clear and accurate information. It is impossible to record the failure surface during vibration. The author tried to record and observe the trend of the failure surface during the test, laid the layer of blue sand every 15 cm, but the result was not satisfactory. The test method for observing the position of the failure surface during vibration needs to be further studied.

Of the various results offered in the literature, it is difficult to failure of the reinforcement in the shaking table test. So the same geogrid was used in two reduced-scale models, and the geogrid was not scaled. The form of the failure surface for the RSRW is not only mentioned above but also a two-wedge failure surface. Due to the use of the location of maximum strain as a basis for judgment, the two-wedge failure surface is not mentioned. Also, the ultimate strain of the reinforcement is generally about 10%, and the measured strain by a strain gauge and an extensometer is within 2%, which is far from the limit. Therefore, the failure surface existing in these planes mentioned above remains to be further explored, and it can be used as a basis for calculation at present.

Specifications for Design of Highway Subgrades [7] (in Chinese) stipulates that the height of RSRW of expressways and first-class highways should not be greater than 12 m, the height of RSRW of second-class and below highway should not be greater than 20 m, and the height of multitiered RSRW should not exceed 10 m. It can be known that the specimens 1 and 2 have different application height ranges. Therefore, specimens 1 and 2 use different model scales and geogrid length. At present, the additional effect of the upper retaining wall to the lower retaining wall is unclear, and the design theory of the multitiered RSRW is not perfect. Therefore, the method of designing the length of the multitiered RSRW needs further take into account.

On comparing the reinforcement strain of different layers in the specimen 1, it can be observed that the reinforcement strain value of the lower wall is greater than that of the upper wall. This condition may be caused by the fact that the reinforcement of the lower retaining wall is required to bear its load and the upper wall load simultaneously, while the reinforcement of the upper wall is only subjected to its load. If the lower wall is regarded as the RSRW, and the upper wall and backfill are regarded as uniform loads, the research results are consistent with those of Wang et al. [17] comparing the same vertical position in the specimen 1. This law is contrary to the displacement distribution of two-tiered RSRW [12]: the strain value at the lower wall is large, the displacement of the lower wall is small, and the strain value at the upper wall is small, the displacement of the upper wall is large. It can be seen that the reinforcement and the modular-block panel cooperate to resist the seismic load. As for how to estimate the magnitude and influencing factors of each part of the earthquake, further research is needed.

Part of the calculation methods (Geoservice method and CTI method) [14] calculate the deformation of the RSRW using the strain state of the reinforcement under normal use conditions. However, there is no calculation method based on the strain state under seismic condition. Whether the deformation of the retaining wall can be inverted to obtain the strain state of the reinforcement under earthquake motion requires further study.

5. Conclusions

To understand the reinforcement stress and potential failure surface distribution of a geogrid RSRW under horizontal seismic loading, the calculation method for the potential failure surface under the dynamic motion of the geogrid RSRW is obtained using the large-scale shaking table tests. The following conclusions could be drawn:(1)The top model blocks of the specimens 1 and 2 both fall due to the excessive peak acceleration. The model block was damaged because the top two or three blocks are not connected rigidly.(2)The strain value of specimen 1 tends to increase with the acceleration increase, and the maximum point of the strain value extends into the interior of the soil. The strain values of different layers at the same position are different. The strain value of the lower retaining wall is greater than that strain of the upper retaining wall.(3)The strain value of the specimen 2 tends to increase as the acceleration increase. When comparing the same-layer strain, the strain law of each layer is consistent. When comparing the strains of different layers at the same location, the strain value law at the same position is nearly the same in all working cases. The strain value of the reinforcement at the lower part of the retaining wall is greater than that at the upper part of the retaining wall.(4)On comparing the measured maximum curve of the geogrid with the existing calculation method for the potential failure surface, the data show that the existing partial calculation method is conservative.(5)Based on the various calculation methods for the potential failure surfaces and measured data, a two-tiered fold-line failure surface is proposed for specimen 1 while taking into consideration the platform width. And it is advised that the failure surface calculation method of BS8006 be used as the calculation method for the potential failure surface of the specimen 2 under dynamic motion.

Data Availability

The data used to support the findings of this study are within the supplementary information files.

Disclosure

This paper and the article “Dynamic Behaviour of reinforced soil retaining wall under horizontal seismic loading” (Li et al. [13]) are the research results of different problems in the same experiment, so there are some similar text.

Conflicts of Interest

The author(s) declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (grant numbers 51778144), the Earthquake Technology Spark Program of China (grant numbers XH204402), and the Fundamental Research Funds for the Central Universities (grant numbers 20180208).

Supplementary Materials

Original data analysis is provided in the file named “Reinforcement strain and potential failure surface.” (Supplementary Materials)