<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Earth Sci.</journal-id>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-6463</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1078624</article-id>
<article-id pub-id-type="doi">10.3389/feart.2022.1078624</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Earth Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The strength law of coir fiber-reinforced soil based on modified Duncan-Chang model</article-title>
<alt-title alt-title-type="left-running-head">Yang et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/feart.2022.1078624">10.3389/feart.2022.1078624</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Hui</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Penghui</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2065742/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>E</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jiang</surname>
<given-names>Xueliang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Jiayu</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Civil Engineering and Engineering Management</institution>, <institution>Guangzhou Maritime University</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Civil Engineering</institution>, <institution>Central South University of Forestry and Technology</institution>, <addr-line>Changsha</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1223883/overview">Fanyu Zhang</ext-link>, Lanzhou University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2037276/overview">Kun Fang</ext-link>, China University of Geosciences Wuhan, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/583010/overview">Armstrong Ighodalo Omoregie</ext-link>, University of Technology Malaysia, Malaysia</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Penghui Li, <email>2015995276@qq.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Geohazards and Georisks, a section of the journal Frontiers in Earth Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>19</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>1078624</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>28</day>
<month>11</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Yang, Li, Yang, Jiang and Chen.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Yang, Li, Yang, Jiang and Chen</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>A large diameter triaxial sample of 61.9&#xa0;mm was made by adding coir fiber into red clay. The range of confining pressure was 50&#x2013;200&#xa0;kPa; fiber content was 0.1%&#x2013;0.4%; and fiber length was 10&#x2013;40&#xa0;mm. By varying the confining pressure, fiber content and fiber length, the unconsolidated and undrained triaxial tests were used to study the shear strength variation law of coir fiber-reinforced soil. The experimental data were processed to establish a linear model of the segmental elastic modulus, and linear analysis was used to determine the model fitting parameters and to improve the Duncan-Chang model by combining the concept of damage ratio. The modified Duncan-Chang model fits the stress-strain relationship of coir-reinforced soil. The results show a clear dividing line for the effect of fiber length and fiber content on the strength of the samples, which is about 30&#xa0;mm and 0.3%, respectively. At the same time, the modified model can fit the stress-strain relationship of coir fiber-reinforced soil and reflect the stress-strain curve characteristics of coir fiber-reinforced soil.</p>
</abstract>
<kwd-group>
<kwd>coir fiber-reinforced soil</kwd>
<kwd>strength law</kwd>
<kwd>elastic modulus</kwd>
<kwd>modified Duncan-Chang model</kwd>
<kwd>triaxial test</kwd>
<kwd>axial strain</kwd>
</kwd-group>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Science and Technology Program of Hunan Province<named-content content-type="fundref-id">10.13039/501100019081</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The ecological environment has always been closely related to human beings. In the field of civil engineering, a series of projects such as urban construction, slope reinforcement, and road planning require necessary soil reinforcement. Finding an ecologically beneficial soil reinforcement technology is significant and has broad application prospects. In recent years, the fiber-reinforced soil method has received more attention and research. Currently, the types of engineering materials are mainly classified as natural and synthetic fibers (<xref ref-type="bibr" rid="B14">Mahdi Hejazi et al., 2012</xref>). Natural fibers include plant fibers, mineral fibers, animal fibers, etc., and synthetic fibers contain carbon fibers, glass fibers, metal fibers, etc. Compared to synthetic fibers, natural fibers are not only easy to obtain but also do not require high cost. Coconut shell fiber is a plant fiber, and its price is only 1/20 of carbon fiber. At the same time, coir fiber has high strength, toughness, and renewable properties and is a natural environmental protection material. Fiber-reinforced red clay soil is a physical way to change the original mechanical properties of the soil (<xref ref-type="bibr" rid="B12">Liu et al., 2013</xref>). The use of coir fiber for the reinforcement of red clay soil is still relatively new, and its reinforcement mechanism and strength law need further research.</p>
<p>Currently, scholars have taken various forms of experimental studies on the mechanical properties of fiber-reinforced soils. Three main aspects are involved. The first is which aspects of the reinforcing effect of the soil after the addition of fibers are reflected. Compaction tests on soil samples after adding coir fibers to the soil revealed that the fibers affected physical properties such as maximum dry density and optimum moisture content of the soil (<xref ref-type="bibr" rid="B9">Jairaj et al., 2018</xref>). Direct shear tests on fiber-reinforced soils revealed that the soil was less brittle, cohesion increased by about 4.3&#x2013;8.8&#xa0;kPa, and internal friction angle increased by about 0.51&#xb0;&#x2013;4.36&#xb0; after fiber addition. The loss of soil strength after reaching its peak is also reduced (<xref ref-type="bibr" rid="B26">Zhang et al., 2001</xref>; <xref ref-type="bibr" rid="B23">Yetimoglu and Salbas, 2003</xref>; <xref ref-type="bibr" rid="B15">Qin et al., 2017</xref>). In addition, adding coir fibers will enhance parameters such as shear strength and modulus of elasticity of the soil. Adding coir fibers to expansive soils increases the strength of soil samples by about 1.1&#x2013;1.37 times (<xref ref-type="bibr" rid="B17">Sivakumar Babu et al., 2008</xref>; <xref ref-type="bibr" rid="B19">Widianti et al., 2021</xref>). Thus, what role do the fibers play in the soils? The fibers are combined with soil particles through friction and cohesion, and the fibers can interweave with each other to form a spatial mesh structure that increases the strength of the sample (<xref ref-type="bibr" rid="B25">Zhang et al., 2021</xref>). The researchers mentioned above mainly used the method of basic tests to study fiber-reinforced soils and less often used more accurate test protocols to test the fiber-reinforced soils. In this paper, a triaxial apparatus was used to study fiber-reinforced soils. The test results were more accurate, facilitating the study of fiber-reinforced soils using a model-fitting approach. The second one is the effect that the different properties of fiber length and content can have on the soil. The unconfined compressive strength of clay soils increases with increasing fiber content, increasing the compressive strength by about 150% compared to the clay strength (193&#xa0;kN/m<sup>2</sup>) (<xref ref-type="bibr" rid="B10">Kumar et al., 2005</xref>). The fiber content strengthens the soil more effectively than the fiber length (<xref ref-type="bibr" rid="B5">Consoli et al., 2007</xref>). The longer the fiber length and the more the incorporation at the same perimeter pressure, the better the reinforcement effect of the fiber (<xref ref-type="bibr" rid="B13">Liu et al., 2011</xref>). Both the tensile and compressive strength of fiber soil samples increased with time (<xref ref-type="bibr" rid="B2">Anggraini et al., 2015</xref>). The above studies only briefly described the relationship between fiber properties and reinforcement effectiveness and did not further explore the specific relationship between the two. For example, whether there is an optimal solution for the reinforcement effect of fiber length and fiber content on the soil will be the focus of the research in this paper. The third one is in the physical model test. Polyester fibers were mixed into the soil to simulate a sloping soil surface layer. The test results showed that the root system simulated by polyester fibers helped to reduce rainwater infiltration, delay the rise of the water table, and improve the soil shear strength (<xref ref-type="bibr" rid="B6">Eab et al., 2014</xref>). Using coir fibers to strengthen the foundation can reduce foundation settlement (<xref ref-type="bibr" rid="B11">Lal et al., 2017</xref>; <xref ref-type="bibr" rid="B7">Fang et al., 2023</xref>). Soils reinforced with vegetation were simulated with a mixture of clay and glass fibers, and the results showed that the deeper the roots, the deeper the sliding surface, with a difference of about twice the depth between the two (<xref ref-type="bibr" rid="B3">Catalina et al., 2022</xref>). These studies provide strong support for fiber-reinforced soils in terms of simulating engineering reality, etc., and are instructive for applying fiber-reinforced soils in engineering practice. It is also a direction for further research on coir fiber-reinforced soils.</p>
<p>Researchers have conducted numerous experimental studies on fiber-reinforced soils regarding fiber types, reinforcement effects, foundation simulation tests, and fiber-soil composition forms. However, the focus of the research is only at the level of experimental tests. The strength properties of fiber-reinforced soils have yet to be investigated by combining the test results with relevant theories. The Duncan-Chang model is essentially a linear stiffness softening the law, while the coir fiber-reinforced soil samples show a non-linear stiffness softening. Therefore, this model could be more accurate for describing the intrinsic relationship of non-linear strain softening. <xref ref-type="bibr" rid="B16">Shen and Zhang et al. (1988)</xref> proposed a theory of soil damage. According to the stress-strain relationship of the conventional triaxial test, the shear deformation of strong structural clay soil is divided into the undamaged structural stage, structural progressive damage stage, and structural complete damage stage (<xref ref-type="bibr" rid="B4">Zhang et al., 2022</xref>). Structural soils vary in their structural yield stresses with test conditions and require constructing a non-linear model of the segmental elastic modulus (<xref ref-type="bibr" rid="B18">Wang, 2006</xref>; <xref ref-type="bibr" rid="B22">Yang and Liang, 2014</xref>). The damage concept proposed by Shen Zhujiang was applied to the original Duncan-Chang model, and the model was improved. The improved hyperbolic function can simulate the stress-strain relationship of structural clay very well (<xref ref-type="bibr" rid="B27">Wang et al., 2004</xref>). Based on the improved Duncan-Chang model, the stress-strain relationship of the original red clay was simulated and calculated. The model could better reflect the development of soil stress-strain at different stages (<xref ref-type="bibr" rid="B8">Gu et al., 2020</xref>). These studies above have indicated clear directions for the improvement of the model, and this paper combines the above approaches to revise the model. In the field of fiber-reinforced soils, fewer studies combine with theoretical models. Based on the modified Duncan-Chang model, this paper introduces the stress-strain relationship data of coir fibers, fits the experimental data of fiber-reinforced soil by improving the Duncan-Chang model, and analyzes its strength law to explore the accuracy of the fitting results of the stress-strain curve of fiber-reinforced soil. These results will provide an essential basis for an in-depth understanding of the strength theory of coir fiber-reinforced soils and provide theoretical support for applying coir fiber-reinforced soils in engineering construction.</p>
</sec>
<sec id="s2">
<title>2 Experimental design</title>
<sec id="s2-1">
<title>2.1 Test materials</title>
<p>The red clay samples were taken from the Central South University of Forestry and Technology area in Changsha, Hunan province, China. The depth of the red clay samples was 50&#x2013;70&#xa0;cm, and the soil was brown-red. According to GB/T50123-2019, the physical properties of the soil were measured. The natural moisture content of undisturbed soil is 25.66%, the liquid limit is 49.0%, the plastic limit is 22.7%, the plastic index is 26.3, the dry density is 1.47&#xa0;g/cm<sup>3</sup>, the specific gravity is 2.74, the maximum dry density is 1.694&#xa0;g/cm<sup>3</sup>, the optimum moisture content is 20.2%. The performance parameters of coir fiber are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Performance parameters of coir fiber.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">No.</th>
<th align="left">Length (mm)</th>
<th align="left">Single ultimate tension (kN)</th>
<th align="left">Diameter (mm)</th>
<th align="left">Tensile strength (MPa)</th>
<th align="left">Ten ultimate tension (kN)</th>
<th align="left">Mean diameter (mm)</th>
<th align="left">Tensile strength (MPa)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="left">60</td>
<td align="left">5.3</td>
<td align="left">0.23</td>
<td align="left">127.56</td>
<td align="left">55.2</td>
<td align="left">0.23</td>
<td align="left">1328.60</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">60</td>
<td align="left">5.1</td>
<td align="left">0.22</td>
<td align="left">134.16</td>
<td align="left">52.7</td>
<td align="left">0.22</td>
<td align="left">1386.36</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-2">
<title>2.2 Sample preparation</title>
<p>The sample sizes commonly used in triaxial studies of fiber-reinforced soils are currently of two types, 39.1 and 50&#xa0;mm in diameter, under restricted test conditions. According to <xref ref-type="bibr" rid="B1">Ang and Erik Loehr (2003)</xref>, the sample diameter exceeds 70&#xa0;mm, and the sample is more likely to represent the strength of fiber-reinforced soil reasonably. Therefore, in order to avoid side effects and to consider the laboratory conditions, the tests in this paper used larger samples than previous sizes to investigate the strong pattern of coir fiber-reinforced soils. The sample diameter is 61.9&#xa0;mm, and the height is 125&#xa0;mm, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. The preparation process was carried out about the Standard for Geotechnical Test Methods (GB/T50123-2019). According to the calculation results, the prepared coir fibers with different contents were added to the soil particles in proportion to the mix well. It was put into the three-valve membrane, in turn, to compact it, thus forming a new soil structure. The distribution of coir fibers in the sample is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. The compaction method was used to make the sample, and the control compaction degree was 95% of the maximum dry density. To simulate the pore state of the soil in the actual project, the prepared samples were put into the vacuum saturator (Apparatus No. 21110155LQ, apparatus located in Room 301 of the Civil Engineering Experiment Center). After the saturation was finished, the samples were encapsulated and stored in cling film.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Coir fiber-reinforced soil with a diameter of 61.9&#xa0;mm and height of 125&#xa0;mm.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g001.tif"/>
</fig>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>The distribution of coir fiber in the sample.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g002.tif"/>
</fig>
</sec>
<sec id="s2-3">
<title>2.3 Test procedure</title>
<p>According to the experimental conditions, the stress-strain relationship of coir fiber-reinforced soil study by the unconsolidated and undrained triaxial shear test method. The experiment was carried out on a Nanjing TKA triaxial tester (Apparatus No. 202001007, apparatus located in Room 301 of the Civil Engineering Experiment Center), which consists of a computer, an axial compression system, a confining pressure system and a pressure chamber. It can automatically control the test process and data acquisition.</p>
<p>The test process is as follows: the coir fiber-reinforced soil sample covered with latex film is placed in the pressure chamber and filled with water. Apply specific confining pressure with a confining pressure system. When the confining pressure reaches a predetermined value, the sample is sheared by a computer-controlled axial compression system at a constant axial loading rate. According to the test code, the shear strain rate was calculated to be 0.625&#xa0;mm/min. A computer collects and processes the data of axial pressure, deflection stress, and axial deformation during loading. The maximum deviated stress value was used as the damage standard of the sample in the test. Four different circumferential pressures were selected for the tests, 50, 100, 150, and 200, respectively. Fiber content was chosen by mass percentages, 0, 0.1%, 0.2%, 0.3%, and 0.4%, respectively. The fiber lengths were 10, 20, 30, and 40&#xa0;mm. The test conditions for the specific samples are shown in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Test conditions.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">No.</th>
<th align="left">Fiber content (%)</th>
<th align="left">Fiber length (mm)</th>
<th align="left">Confining pressure (kPa)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="left">0</td>
<td align="left">0</td>
<td align="left">100</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">0.1, 0.2, 0.3, 0.4</td>
<td align="left">30</td>
<td align="left">100</td>
</tr>
<tr>
<td align="left">3</td>
<td align="left">0.3</td>
<td align="left">10, 20, 40</td>
<td align="left">100</td>
</tr>
<tr>
<td align="left">4</td>
<td align="left">0.3</td>
<td align="left">30</td>
<td align="left">50, 150, 200</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s3">
<title>3 The arrangement and analysis of test data</title>
<sec id="s3-1">
<title>3.1 Failure mode of coir fiber-reinforced soil</title>
<p>Influenced by the coir fiber, the samples underwent obvious plastic deformation during the test. As shown in <xref ref-type="fig" rid="F3">Figure 3A</xref>, without adding coir fibers, the soil samples showed an obvious diagonal crack during the damage, which exhibited brittle damage; as shown in <xref ref-type="fig" rid="F3">Figure 3B</xref>, after the addition of fibers, the damage showed bulging damage. Although a certain shear surface was also generated, the degree of shear was substantially reduced compared with that without the addition of coir fibers. Therefore, the damaged form of the coir fiber-reinforced soil is a mixture of ductile and brittle damage. This intuitively indicates that coir fiber can reduce the degree of damage to soil samples.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Mode of sample failure. <bold>(A)</bold> The red clay samples, <bold>(B)</bold> Coir fiber-reinforced soil.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g003.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>3.2 The variation law of deviating stress and strain</title>
<sec id="s3-2-1">
<title>3.2.1 Effect of fiber length</title>
<p>
<xref ref-type="fig" rid="F4">Figure 4</xref> shows the q-&#x3b5; (q &#x3d; &#x3c3;<sub>1</sub>&#x2013;&#x3c3;<sub>3</sub>, q is deviator stress) stress-strain curve of coir fiber at 100&#xa0;kPa confining pressure and 0.3% fiber content at different fiber lengths. During the test, it was found that the fiber length over 50&#xa0;mm is easy to twist into clumps, which will make the sample making more difficult and the strength of the sample is significantly reduced. Therefore, four different fiber lengths of 10, 20, 30, and 40&#xa0;mm were used in this paper. The values of the deviator stress in the stress-strain relationship curves differed significantly for the samples with different test parameters. When the fiber length increased in the range of 0&#x2013;30&#xa0;mm, the deviator stress of the samples gradually increased with the increasing length of fiber. When the fiber length exceeded 30&#xa0;mm, further increasing the length of coir fiber would not increase the deviator stress. When the fiber length was 40&#xa0;mm, the deviator stress of the sample at the same axial strain position was lower than that of the 20&#xa0;mm fiber length. This indicates that the effect of fiber length on the strength of the sample is pronounced, and the excessive fiber length causes the fibers themselves to associate with each other and form a dense web-like structure. This hinders the formation of a composite structure between the soil and the coir fibers, which results in the reinforcing effect of the fibers not being fully developed.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Effect of fiber length within 0&#x2013;40&#xa0;mm on stress-strain of samples at constant confining pressure and fiber content.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g004.tif"/>
</fig>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Effect of fiber content</title>
<p>
<xref ref-type="fig" rid="F5">Figure 5</xref> shows the stress-strain curves of different fiber contents under confining pressure of 100&#xa0;kPa and fiber length of 30&#xa0;mm. The percentage of fiber mass to soil sample mass was used as the fiber content to investigate the effect of fiber mass on the strength of the samples. When the fiber content reaches 0.5%, the excessive fibers cause them to aggregate with each other, and the produced samples easy form weak zones. Therefore, four different fiber contents of 0.1%, 0.2%, 0.3%, and 0.4% were used. As the fiber content grew from 0% to 0.3%, the deviator stress gradually increased. However, when the fiber content exceeded 0.3%, the deviator stress of the samples started to decrease, and then increasing the amount of coir fiber participation would not improve the strength of the samples.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Effect on stress-strain of samples at 0%&#x2013;0.4% of fiber content at constant confining pressure and fiber length.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g005.tif"/>
</fig>
<p>The content and length of coir fiber affect the strength of coir fiber, but they have different effects. When the fiber content increases to 0.4%, the deviator stress at the same axial strain position is lower than 0.3% fiber content but higher than that at 0.2% fiber content. This indicates that the excessive fiber content will also negatively affect the sample&#x2019;s strength, but the effect is weaker than the fiber length. Too much fiber also leads to fiber binding, which is not conducive to the connection between soil particles and fiber. Therefore, the effect of coir fiber on red clay is equivalent to the effect of reinforced material on cement, and there is an optimal content and length.</p>
<p>
<xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref> show that the strength of the sample changes with increasing fiber length and content. When no fiber is added, the stress-strain relationship is a softening curve. When fibers are added, the stress-strain relationship is a hardening type of curve, and the deviator stress increases with the development of axial strain. The coir fiber has a good anchoring effect. When the fiber action, the value of the deviator stress of the sample was significantly larger than the sample without the addition of fiber. When the axial strain did not reach 1.5%, both strengths increased with strain development. With the axial strain&#x2019;s further development, the fibre sample&#x2019;s deviator stress continued to grow. For the samples without added fiber, the value of deviator stress started to decrease, and the difference between the two gradually increased. As the strain grows, the deviator stress&#x2019;s differential value also increases. This indicates that the fibers are becoming more and more effective in reinforcing red clay soils.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Effect of stress-strain on samples within 50&#x2013;200&#xa0;kPa of confining pressure with constant fiber length and fiber content.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Relation between axial strain ratio deviator stress and axial strain.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g007.tif"/>
</fig>
</sec>
<sec id="s3-2-3">
<title>3.2.3 Effect of confining pressure</title>
<p>
<xref ref-type="fig" rid="F6">Figure 6</xref> shows the stress-strain curve of samples under different confining pressures with 0.3% coir fiber content. The depth of fiber soil cover in construction engineering is about 0&#x2013;15&#xa0;m, and the soil pressure is determined as the product of dry density (1.47&#xa0;g/cm<sup>3</sup>) and gravitational acceleration (9.8&#xa0;m/s<sup>2</sup>) and depth. Therefore, the range of soil pressure is 0&#x2013;220&#xa0;kPa, so four levels of 50, 100, 150, and 200&#xa0;kPa are used to simulate the actual working conditions of the fiber-reinforced soil. With the same fiber length and fiber content, the deviator stress increases with the step-by-step increase of the confining pressure.</p>
</sec>
</sec>
<sec id="s3-3">
<title>3.3 Discussion and analysis</title>
<p>After the above analysis, the strength properties of the fiber-reinforced soil samples exhibited the optimum situation at a fiber length of 30&#xa0;mm and a fiber content of 0.3%, and their stress-strain curves are shown in <xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref>. The stress-strain curves of the fiber-reinforced soil samples under different test parameters were similar in shape, and all of them retained the essential characteristics of hyperbolic curve. For the fiber-reinforced soil samples, the damage stage is significantly different from that of the red clay soil. After the red clay soil samples reached the damage stage, the stress decreased with the increase of axial strain. In contrast, after the fiber-reinforced soil reached the damage stage, the stress continued to increase with the growth of axial strain, and the damage form showed a hardening characteristic. After the addition of fibers to the soil, the combination of the two forms a new structure. The incorporation of fibers has a significant impact on the mechanical properties of the sample. Soil particles in contact with fibers form a new arrangement. Their structural composition is changed, manifested by a decrease in density and compaction, and a substantial increase in plastic deformation capabilities. The stress-strain characteristics are as follows: when the plastic deformation is small, the stress-strain curve rises faster, or the stress increase is more prominent. When the plastic deformation is larger, the deviator stress of some samples appears to fall and then increase. With the further growth of strain, the stress-strain curve gradually tends to flatten out, and the deviator stress will keep growing slowly. Therefore, adding fiber can effectively limit the deformation of the soil, improve the soil&#x2019;s strength index, and effectively change the soil&#x2019;s original stress-strain field, thus optimizing the engineering performance of the soil.</p>
</sec>
</sec>
<sec id="s4">
<title>4 A modified Duncan-Chang model of coir fiber-reinforced soil</title>
<sec id="s4-1">
<title>4.1 Modified Duncan-Chang model</title>
<p>The structural damage law promoted by <xref ref-type="bibr" rid="B16">Shen and Zhang (1998)</xref> regards the structural state of soil as two different forms: the non-damaged state and the remodeled state. The soil body is in a combination of both the non-damaged state and the remodeled state when the force deforms it, and with the gradual increase of the load, the state in which the soil sample is in can be expressed by the formula<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>In the formula: S is the mechanical parameter of the soil; S<sub>i</sub> is the strength property of the soil in the lossless state; S<sub>d</sub> is the strength property of the soil in the remodeled state; &#x3c9; is the damage ratio, the proportion of the soil in the remodeled state to the whole soil.</p>
<p>Eq. <xref ref-type="disp-formula" rid="e1">1</xref> is brought into the Duncan-Chang model to derive the modified hyperbolic function<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mfrac>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>In the formula: <inline-formula id="inf1">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf2">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf3">
<mml:math id="m5">
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>; E<sub>i1</sub> and E<sub>i2</sub> are the elastic modulus of the first and second stages, respectively; (&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)<sub>ult1</sub> and (&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)<sub>ult2</sub> are the deviator stress limit values of the fitted straight lines of the first and second stages, respectively; &#x3c9; is the damage ratio; &#x3b5;<sub>1u</sub> is the axial strain corresponding to the soil body at the time of damage.</p>
<p>By analyzing the triaxial test data, the undamaged structural stage and the structural gradually damaged stage of the reinforced soil sample are regarded as the first stage of strength deformation; then the structural completely damaged stage is made the second stage of strength deformation.</p>
<p>In <xref ref-type="fig" rid="F7">Figure 7</xref>, the stress-strain relationship curve &#x3b5;<sub>1</sub>/(&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>) &#x223c; &#x3b5;<sub>1</sub> consists of two stages, &#x3b5;<sub>1</sub> is at 0%&#x2013;0.5% and &#x3b5;<sub>1</sub> is greater than 0.5% after. Two straight lines can be fitted to these two stages respectively. Since the yield stress state of the soil changes with the change in loading conditions. As the stress changes, the sample gradually produces damage, and when the stress changes to a large enough, the damage will be more obvious, which is shown as the appearance of a turning point on the &#x3b5;<sub>1</sub>/(&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>) &#x223c; &#x3b5;<sub>1</sub> curve, and this turning point is regarded as the turning point of the two straight lines of the first and second stages. The turning point&#x2019;s appearance indicates that the sample&#x2019;s damage has developed from the first stage to the second stage, and its elastic modulus and ultimate deviator stress will also change significantly, which is presented in the coordinate system in <xref ref-type="fig" rid="F7">Figure 7</xref> as the change of intercept and slope.</p>
<p>As shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, each stage can be fitted linearly, and the corner mark k refers to two different steps in the hardening process. The formula is expressed as<disp-formula id="e3">
<mml:math id="m6">
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="normal">k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>The asymptotic values of elastic modulus and deviator stress at each stage are obtained from the following equations.<disp-formula id="e4">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
</sec>
<sec id="s4-2">
<title>4.2 Modified Duncan-Chang model parameter solution</title>
<p>In this paper, the stress-strain curve of coir fiber-reinforced soil is fitted using the hyperbolic function of Eq. <xref ref-type="disp-formula" rid="e2">2</xref>. First, a part of the data with abnormal values will appear at the beginning of the test. After removing these anomalous data, the relationship image of &#x3b5;<sub>1</sub>/(&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)&#x223c;&#x3b5;<sub>1</sub> was plotted by data processing software. The stress-strain relationship is shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, and there are two obvious phases in the image. Point P is the turning point. A straight-line fit was performed for each of the two phases, and the fit was above 97%. Then, the values of parameters E<sub>i1</sub>, E<sub>i2</sub>, (&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)<sub>ult1</sub>, and (&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)<sub>ult2</sub> are determined according to the Eqs <xref ref-type="disp-formula" rid="e3">3</xref>&#x2013;<xref ref-type="disp-formula" rid="e5">5</xref>. The same method was adopted to fit the data of each sample into a straight line to find out the values of the above four parameters. The calculated parameters are shown in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Stress-strain calculation parameters of coir fiber-reinforced soil.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="left">Model parameters</th>
<th colspan="5" align="left">Fiber content/%</th>
<th colspan="4" align="left">Fiber length/mm</th>
<th colspan="4" align="left">Confining pressure/kPa</th>
</tr>
<tr>
<th align="left">0</th>
<th align="left">0.1</th>
<th align="left">0.2</th>
<th align="left">0.3</th>
<th align="left">0.4</th>
<th align="left">10</th>
<th align="left">20</th>
<th align="left">30</th>
<th align="left">40</th>
<th align="left">50</th>
<th align="left">100</th>
<th align="left">150</th>
<th align="left">200</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">10<sup>3</sup>/E<sub>i1</sub>
</td>
<td align="left">2.6</td>
<td align="left">2.5</td>
<td align="left">2.0</td>
<td align="left">1.9</td>
<td align="left">2.0</td>
<td align="left">2.8</td>
<td align="left">2.3</td>
<td align="left">1.9</td>
<td align="left">2.3</td>
<td align="left">2.2</td>
<td align="left">1.9</td>
<td align="left">1.7</td>
<td align="left">1.4</td>
</tr>
<tr>
<td align="left">10<sup>3</sup>/(&#x3c3;<sub>1</sub>&#x2013;&#x3c3;<sub>3</sub>)<sub>ult1</sub>
</td>
<td align="left">7.1</td>
<td align="left">4.5</td>
<td align="left">4.9</td>
<td align="left">4.0</td>
<td align="left">4.3</td>
<td align="left">5.3</td>
<td align="left">4.2</td>
<td align="left">4.0</td>
<td align="left">5.1</td>
<td align="left">5.2</td>
<td align="left">4.0</td>
<td align="left">3.6</td>
<td align="left">2.8</td>
</tr>
<tr>
<td align="left">10<sup>3</sup>/E<sub>i2</sub>
</td>
<td align="left">2.2</td>
<td align="left">1.3</td>
<td align="left">1.2</td>
<td align="left">1.1</td>
<td align="left">1.3</td>
<td align="left">2.1</td>
<td align="left">1.9</td>
<td align="left">1.1</td>
<td align="left">1.6</td>
<td align="left">1.7</td>
<td align="left">1.1</td>
<td align="left">1.0</td>
<td align="left">0.8</td>
</tr>
<tr>
<td align="left">10<sup>3</sup>/(&#x3c3;<sub>1</sub>&#x2013;&#x3c3;<sub>3</sub>)<sub>ult2</sub>
</td>
<td align="left">11.8</td>
<td align="left">7.8</td>
<td align="left">6.9</td>
<td align="left">5.8</td>
<td align="left">6.5</td>
<td align="left">6.2</td>
<td align="left">5.7</td>
<td align="left">5.8</td>
<td align="left">6.2</td>
<td align="left">7.1</td>
<td align="left">5.8</td>
<td align="left">4.5</td>
<td align="left">3.7</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The elastic modulus of stage 1 and stage 2 are shown in <xref ref-type="table" rid="T3">Table 3</xref>, showing the relationship between elastic moduli and fiber content, fiber length, and confining pressure under different test parameters. With the increase in fiber content and length, the elastic modulus of coir fiber-reinforced soil shows a trend of increasing first and then decreasing. The elastic modulus reaches the maximum value when the fiber content and length reach 0.3% and 30&#xa0;mm, respectively. If the fiber content and length continue to increase, the elastic modulus decreases, and the ability to resist deformation of the sample decreases. For the same fiber content and length, the higher the confining pressure, the greater the elastic modulus. That is, the higher the resistance to deformation.</p>
<p>The values of ultimate deviator stress in stage 1 and stage 2 are shown in <xref ref-type="table" rid="T3">Table 3</xref>, showing the relationship between ultimate deviator stress and fiber content, fiber length, and confining pressure under different test parameters. The ultimate deviator stress value variation gradually increases when the fiber content is from 0.1% to 0.3%. Similarly, the ultimate deviator stress increases gradually when the fiber length is 10&#x2013;30&#xa0;mm. As the fiber content and fiber length continue to grow, the strength of the sample begins to decrease after reaching the optimum content and optimum length. The relationship between the modulus of elasticity and the confining pressure is shown by the fact that for the same sample with the same fiber content and fiber length, the value of the ultimate deviator stress increases with increasing confining pressure, and the sample exhibits higher strength. As can be seen from <xref ref-type="table" rid="T3">Table 3</xref>, the values of ultimate deviator stress (&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)<sub>ult1</sub> for the same group of samples are all higher than (&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)<sub>ult2</sub>, indicating that the samples reinforced with coir fiber are suitable for modified Duncan-Chang model.</p>
<p>Based on the model and the pattern exhibited by the experimental results, the sample deviator stress also increases when the fiber content increases from 0.1% to 0.3%. When the fiber content was 0.4%, the sample deviator stress still kept increasing, but the value was less than 0.3% of fiber content. It indicates that the fiber content significantly changes the reinforcement effect of red clay samples at 0.3%. Similarly, based on the test conditions and experimental results, it can be concluded that the fiber length is optimal when it reaches 30&#xa0;mm.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Discussion of the model fitting results</title>
<p>The experimental results and the fitted results are plotted on the same coordinate axes as in <xref ref-type="fig" rid="F8">Figure 8</xref>. The experimental data of the coir fiber-reinforced soil samples were consistent with the fitted data of the improved model. In the case of small strains (&#x3b5;<sub>1</sub> &#x3c; 1%), the fitted results are slightly larger compared to the measured values, but there is no significant deviation in general. The use of the modified model to fit the coir fiber-reinforced soil has the following advantages: its variation law is consistent with the actual situation, and the parameters of the modified model are obtained in a simple and clear sense, which provides theoretical guidance for the study of the strength law of coir fiber-reinforced soil.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Comparison of data at 0.3% fiber content with model fitting results.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g008.tif"/>
</fig>
<p>In order to explore the accuracy of the modified Duncan-Chang model in fitting the remaining coir fiber-reinforced soil samples, the test data of several groups of samples were fitted. The fit results are shown in <xref ref-type="fig" rid="F9">Figures 9</xref>&#x2013;<xref ref-type="fig" rid="F11">11</xref>.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Comparison of data with model fitting results for fiber content 0%&#x2013;0.4%.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Comparison of model fitting results for confining pressure 50&#x2013;200&#xa0;kPa.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Comparison of model fitting results for fiber length 0&#x2013;40&#xa0;mm.</p>
</caption>
<graphic xlink:href="feart-10-1078624-g011.tif"/>
</fig>
<p>In the axial deformation at about 1%, the fitted curve of the sample is higher than the curve drawn from the test data. The reason is that this part is in the junction area between the first and second stages of the &#x3b5;<sub>1</sub>/(&#x3c3;<sub>1</sub>-&#x3c3;<sub>3</sub>)&#x223c;&#x3b5;<sub>1</sub> curve, and this area contains the turning point of the two fitted straight lines in the modified model, which leads to a change in the values of the model parameters in this part, thus causing a difference between the fitted curve and the experimental data here. When the deformation of the sample develops from the first stage to the second stage, the solution of the fitted data changes from a straight line in the first stage to a straight line in the second stage. At this time, the intercept of the line becomes smaller and the slope increases. The modulus of elasticity increases while the ultimate deviator stress decreases. For the second stage, the fit between the model fitting curve and the test data curve is higher in the case of larger axial strain.</p>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>In this paper, larger triaxial samples are used. The experimental demonstration and model validation are provided for studying the strong law of coir fiber-reinforced soil by comparing the experimental data with the fitted curve of the modified Duncan-Chang model. At the same time, it provides a significant reference value for the strength theory research of coir fiber-reinforced soil. The following conclusions were drawn:<list list-type="simple">
<list-item>
<p>(1) At the beginning of loading, the eccentric stress increases rapidly when the axial strain is minor. With the increase of the axial strain, the curve tends to be smooth. The comparison of coir fiber-reinforced soil and red clay shows that the difference in deviator stress of different samples gradually increases, and the effect of coir fiber-reinforced soil is remarkable.</p>
</list-item>
<list-item>
<p>(2) The ultimate deviating stress and elastic modulus increase and then decrease with the increase of coir fiber content and length. The effect of fiber length and fiber content on the strength of coir fiber-reinforced soil has an obvious boundary, the cut-off length and the amount of doping are about 30&#xa0;mm and 0.3%, respectively. When it is below this limit, the elastic modulus and ultimate deviator stress of coir fiber-reinforced soil increase with the increase of fiber length and fiber content, and when it is above this limit, they begin to decrease.</p>
</list-item>
<list-item>
<p>(3) Under the condition of the same fiber length and content, the strength of coir fiber-reinforced soil under the same axial strain increases with the increase of confining pressure step by step. That is to say, the greater the confining pressure, the higher the sample&#x2019;s strength.</p>
</list-item>
<list-item>
<p>(4) From the point of view of structural damage, the stress-strain relationship of coir fiber-reinforced soil was divided into two stages, which the modified Duncan-Chang model fitted.</p>
</list-item>
<list-item>
<p>(5) The calculation of the formula obtains the exact parameters of the fit model. The modified Duncan-Chang model can fit the stress-strain relationship of coir fiber-reinforced soil, and the fit degree is high.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>Performed research, analyzed data, and wrote the paper, PL; designed research and revisions, HY; designed research and revisions, XJ; analyzed data, EY; analyzed data, JC.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>Five projects supported the study, which are the National Natural Science Foundation of China (31971727), the Science and Technology Innovation Program of Hunan Province (Grant No.: 2022NK2056), the Forest Science and Technology Innovation Program of Hunan Province (XLK202105-3), the Special Fund Project of Safety Production Prevention and Emergency Response of Hunan Province in 2021 (Grant No.: 2021-QYC-10008-24956) and the Hunan Provincial Natural Science Foundation Project (Grant No.: 2022JJ31005).</p>
</sec>
<ack>
<p>The authors are very grateful to the editor and reviewers for their valuable advice.</p>
</ack>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ang</surname>
<given-names>E. C.</given-names>
</name>
<name>
<surname>Erik Loehr</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Sample size effects for fiber-reinforced silty clay in unconfined compression[J]</article-title>. <source>Geotechnical Test. J.</source> <volume>26</volume> (<issue>2</issue>), <fpage>191</fpage>&#x2013;<lpage>200</lpage>. <pub-id pub-id-type="doi">10.1520/GTJ11320J</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anggraini</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Asadi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Huat</surname>
<given-names>B. B. K.</given-names>
</name>
<name>
<surname>Nahazanan</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Effects of coir fibers on tensile and compressive strength of lime treated soft soil</article-title>. <source>Measurement</source> <volume>59</volume>, <fpage>372</fpage>&#x2013;<lpage>381</lpage>. <pub-id pub-id-type="doi">10.1016/j.measurement.2014.09.059</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Catalina</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Mendoza</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Jose Vicente</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Physical and numerical modeling of clayey slopes reinforced with roots[J]</article-title>. <source>Int. J. Civ. Eng.</source> <volume>20</volume> (<issue>9</issue>), <fpage>1115</fpage>&#x2013;<lpage>1128</lpage>. <pub-id pub-id-type="doi">10.1007/S40999-022-00733-0</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Consoli</surname>
<given-names>N. C.</given-names>
</name>
<name>
<surname>Casagrande</surname>
<given-names>M. D. T.</given-names>
</name>
<name>
<surname>Coop</surname>
<given-names>M. R.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Performance of a fibre-reinforced sand at large shear strains</article-title>. <source>G&#xe9;otechnique</source> <volume>57</volume>, <fpage>751</fpage>&#x2013;<lpage>756</lpage>. <pub-id pub-id-type="doi">10.1680/geot.2007.57.9.751</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eab</surname>
<given-names>K. H.</given-names>
</name>
<name>
<surname>Takahashi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Likitlersuang</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Centrifuge modelling of root-reinforced soil slope subjected to rainfall infiltration</article-title>. <source>G&#xe9;otechnique Lett.</source> <volume>4</volume> (<issue>3</issue>), <fpage>211</fpage>&#x2013;<lpage>216</lpage>. <pub-id pub-id-type="doi">10.1680/geolett.14.00029</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Centrifuge modelling of landslides and landslide hazard mitigation: A review</article-title>. <source>Geosci. Front.</source> <volume>14</volume> (<issue>1</issue>), <fpage>101493</fpage>. <comment>ISSN 1674-9871</comment>. <pub-id pub-id-type="doi">10.1016/j.gsf.2022.101493</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Simulation of the stress-strain relationship of red clay based on the modified Duncan-Chang model [J]</article-title>. <source>J. Guilin Univ. Technol.</source> <volume>40</volume> (<issue>02</issue>), <fpage>351</fpage>&#x2013;<lpage>357</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1674-9057.2020.02.013</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jairaj</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Prathap Kumar</surname>
<given-names>M. T.</given-names>
</name>
<name>
<surname>Raghunandan</surname>
<given-names>M. E.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Compaction characteristics and strength of BC soil reinforced with untreated and treated coir fibers</article-title>. <source>Innov. Infrastruct. Solut.</source> <volume>3</volume> (<issue>1</issue>), <fpage>21</fpage>. <pub-id pub-id-type="doi">10.1007/s41062-017-0123-2</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumar</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Singh Walia</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Mohan</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Compressive strength of fiber reinforced highly compressible clay</article-title>. <source>Constr. Build. Mater.</source> <volume>20</volume> (<issue>10</issue>), <fpage>1063</fpage>&#x2013;<lpage>1068</lpage>. <pub-id pub-id-type="doi">10.1016/j.conbuildmat.2005.02.027</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lal</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Sankar</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chandrakaran</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Effect of reinforcement form on the behaviour of coir geotextile reinforced sand beds</article-title>. <source>Soils Found.</source> <volume>57</volume> (<issue>2</issue>), <fpage>227</fpage>&#x2013;<lpage>236</lpage>. <pub-id pub-id-type="doi">10.1016/j.sandf.2016.12.001</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Jian</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>etc. Research progress on engineering properties of fiber-reinforced soil</article-title>. <source>J. Eng. Geol.</source>, <fpage>540</fpage>&#x2013;<lpage>547</lpage>. <pub-id pub-id-type="doi">10.3969/J.ISSN.1004-9665.2013.04.009</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xiurun</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Triaxial test study on glass fiber-reinforced soil [J]</article-title>. <source>J. Shanghai Jiao Tong Univ.</source> <volume>45</volume> (<issue>05</issue>), <fpage>762</fpage>&#x2013;<lpage>766&#x2b;771</lpage>. <pub-id pub-id-type="doi">10.16183/J.CNKI.JSJTU.2011.05.028</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mahdi Hejazi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sheikhzadeh</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Abtahi</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>A simple review of soil reinforcement by using natural and synthetic fibers</article-title>. <source>Constr. Build. Mater.</source> <volume>30</volume>, <fpage>100</fpage>&#x2013;<lpage>116</lpage>. <pub-id pub-id-type="doi">10.1016/j.conbuildmat.2011.11.045</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qin</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Triaxial test study on Coir soil [J]</article-title>. <source>Sci. Technol. Eng.</source> <volume>17</volume> (<issue>10</issue>), <fpage>272</fpage>&#x2013;<lpage>276</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1671-1815.2017.10.047</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Shen</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>1988</year>). &#x201c;<article-title>Applying damage mechanics to soil mechanics [C]</article-title>,&#x201d; in <conf-name>3rd Proceedings of the National Symposium on Numerical Analysis and analytical methods of geotechnical mechanics</conf-name>, <fpage>595</fpage>&#x2013;<lpage>609</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sivakumar Babu</surname>
<given-names>G. L.</given-names>
</name>
<name>
<surname>Vasudevan</surname>
<given-names>A. K.</given-names>
</name>
<name>
<surname>Sayida</surname>
<given-names>M. K.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Use of coir fibers for improving the engineering properties of expansive soils</article-title>. <source>J. Nat. Fibers</source> <volume>5</volume> (<issue>1</issue>), <fpage>61</fpage>&#x2013;<lpage>75</lpage>. <pub-id pub-id-type="doi">10.1080/15440470801901522</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>L. Z.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C. Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L. L.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>A modified Duncan-Zhang model considering soil structural properties [J]</article-title>. <source>J. Water Resour.</source> (<issue>01</issue>), <fpage>83</fpage>&#x2013;<lpage>89</lpage>. <pub-id pub-id-type="doi">10.3321/j.issn:0559-9350.2004.01.016</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Study on piecewise tangent modulus of nonlinear soil model [J]</article-title>. <source>J. Hohai Univ. Sci.</source> (<issue>02</issue>), <fpage>204</fpage>&#x2013;<lpage>207</lpage>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Widianti</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Diana</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Alghifari</surname>
<given-names>M. R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Shear strength and elastic modulus behavior of coconut fiber-reinforced expansive soil</article-title>. <source>IOP Conf. Ser. Mat. Sci. Eng.</source> <volume>1144</volume> (<issue>1</issue>), <fpage>012043</fpage>. <pub-id pub-id-type="doi">10.1088/1757-899x/1144/1/012043</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Study on the stress-strain relationship of soft structural soil based on Duncan-Chang model [J]</article-title>. <source>Hydrogeology Eng. Geol.</source> <volume>41</volume> (<issue>04</issue>), <fpage>75</fpage>&#x2013;<lpage>79&#x2b;86</lpage>. <pub-id pub-id-type="doi">10.16030/J.CNKI.ISSN.1000-3665.2014.04.012</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yetimoglu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Salbas</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>A study on shear strength of sands reinforced with randomly distributed discrete fibers</article-title>. <source>Geotext. Geomembranes</source> <volume>21</volume> (<issue>2</issue>), <fpage>103</fpage>&#x2013;<lpage>110</lpage>. <pub-id pub-id-type="doi">10.1016/s0266-1144(03)00003-7</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Experimental study on strength characteristics of palm fiber-reinforced soil [J]</article-title>. <source>Mineral. Explor.</source> <volume>12</volume> (<issue>05</issue>), <fpage>1264</fpage>&#x2013;<lpage>1271</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1674-7801.2021.05.020</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Experimental study on strength characteristics of fibrous soil [J]</article-title>. <source>Roadbed Eng.</source> (<issue>01</issue>), <fpage>36</fpage>&#x2013;<lpage>38</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1003-8825.2001.01.012</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Q-M.</given-names>
</name>
<name>
<surname>Pei</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A modified Duncan-Chang model for strongly structural clays [J]</article-title>. <source>J. Water Conservancy Water Transp. Eng.</source> (<issue>02</issue>), <fpage>117</fpage>&#x2013;<lpage>125</lpage>. <pub-id pub-id-type="doi">10.12170/20210617003</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>