Acessibilidade / Reportar erro

Geotechnical site characterization by DMT and laboratory tests on an unsaturated tropical soil site for slope stability analysis

Abstract

The slope stability is an important topic because it presents risks of socio-economic losses caused by eventual ruptures. It is necessary to identify the site profile, as well as obtaining soil strength parameters for the slope stability analysis. This paper presents and discusses the use of the Flat Dilatometer Test (DMT) in the geotechnical site characterization for slope stability analysis in an unsaturated tropical soil site. Six DMTs were performed to define the stratigraphical profile and estimate geotechnical parameters. Shear strength parameters were determined in the laboratory using saturated and unsaturated triaxial compression as well as soil water retention curves (SWRC) to support DMT data interpretation. A commercial software was used to perform the slope stability analysis a cut slope with 6.6 m height and a gradient around 55° to illustrate the application of DMT and triaxial test data. The DMT allowed the representative site profile to be identified, as well as estimating the design parameters that compared well to those interpreted from unsaturated triaxial test data for the in-situ soil suction. The DMT can be used as logging test in the preliminary characterization of studied site specially to define the stratigraphical profile, site variability, select the regions to collect disturbed and undisturbed soil samples and as the first attempt to estimate the geotechnical design parameters via correlations. It is important to emphasize that laboratory tests on undisturbed soil samples are essential in the slope stability analysis of unsaturated tropical soil profiles.

Keywords:
Site characterization; Slope stability; In situ tests; DMT; Unsaturated soils; Suction

1. Introduction

Slopes can be natural or man-made, and their stability is a very important issue because of the risks of socio-economic losses caused by eventual failures. Slope failures can damage buildings, road infrastructure, as well as cause injury or death to people (Rahardjo et al., 2019Rahardjo, H., Kim, Y., & Satyanaga, A. (2019). Role of unsaturated soil mechanics in geotechnical engineering. International Journal of Geo-Engineering, 10(1), 8. http://dx.doi.org/10.1186/s40703-019-0104-8.
http://dx.doi.org/10.1186/s40703-019-010...
). One of the relevant aspects that leads to the complexity of slope stability is that it can be formed by rock, soil, or both, requiring a proper understanding of the stress state and failure mechanisms.

The need for stabilization and adequate preventive measures for slope protection has gained great relevance in the context of sustainable development of urban areas, especially in times of climate change. In the light of the effects of climate change, longer periods of drought and higher intensity and shorter duration rainfall are expected in the future, and it is necessary to prevent future rain-induced slope failures for slope protection against rainfall.

Changes in rainfall patterns, in particular, will influence the boundary condition of flow at the ground surface. Changes in groundwater hydrology associated with extreme events at the soil-atmosphere interface demonstrate the complex nature of time-dependent problems such as slope stability that can only be best analyzed in the context of Unsaturated Soil Mechanics.

The behavior of unsaturated soils depends on the soil suction (Alonso et al., 1990Alonso, E.E., Gens, A., & Josa, A. (1990). A constitutive model for partially saturated soils. Geotechnique, 40(3), 405-430. http://dx.doi.org/10.1680/geot.1990.40.3.405.
http://dx.doi.org/10.1680/geot.1990.40.3...
). Changes in the water content and suction in an unsaturated soil occurs due to climate variations (Blight, 2003Blight, G.E. (2003). The vadose zone soil-water balance and transpiration rates of vegetation. Geotechnique, 53(1), 55-64. http://dx.doi.org/10.1680/geot.2003.53.1.55.
http://dx.doi.org/10.1680/geot.2003.53.1...
; Cui et al., 2008Cui, Y.-J., Tang, A.-M., Mantho, A.T., & De Laure, E. (2008). Monitoring field soil suction using a miniature tensiometer. Geotechnical Testing Journal, 31(1), 100769. http://dx.doi.org/10.1520/GTJ100769.
http://dx.doi.org/10.1520/GTJ100769...
). Climate events (i.e., extreme precipitation and droughts) are dominated by soil-atmospheric interactions in which soil water content and soil suction can change seasonally (Rahardjo et al., 2019Rahardjo, H., Kim, Y., & Satyanaga, A. (2019). Role of unsaturated soil mechanics in geotechnical engineering. International Journal of Geo-Engineering, 10(1), 8. http://dx.doi.org/10.1186/s40703-019-0104-8.
http://dx.doi.org/10.1186/s40703-019-010...
). The changes in soil suction could reduce the shear strength of soil that may result in rainfall-induced slope failures and consequent destructions and deaths (Cho & Lee, 2002Cho, S.E., & Lee, S.R. (2002). Evaluation of surficial stability for homogeneous slopes considering rainfall characteristics. Journal of Geotechnical and Geoenvironmental Engineering, 128(9), 756-763. http://dx.doi.org/10.1061/(ASCE)1090-0241(2002)128:9(756).
http://dx.doi.org/10.1061/(ASCE)1090-024...
; Ng & Shi, 1998Ng, C.W., & Shi, Q. (1998). A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage. Computers and Geotechnics, 22(1), 1-28. http://dx.doi.org/10.1016/S0266-352X(97)00036-0.
http://dx.doi.org/10.1016/S0266-352X(97)...
; Tohari et al., 2007Tohari, A., Nishigaki, M., & Komatsu, M. (2007). Laboratory rainfall-induced slope failure with moisture content measurement. Journal of Geotechnical and Geoenvironmental Engineering, 133(5), 575-587. http://dx.doi.org/10.1061/(ASCE)1090-0241(2007)133:5(575).
http://dx.doi.org/10.1061/(ASCE)1090-024...
). Therefore, the effects of unsaturated soils should be considered in slope stability analyses for geotechnical designs (Alonso et al., 2003Alonso, E.E., Gens, A., & Delahaye, C.H. (2003). Influence of rainfall on the deformation and stability of a slope in overconsolidated clays: a case study. Hydrogeology Journal, 11(1), 174-192. http://dx.doi.org/10.1007/s10040-002-0245-1.
http://dx.doi.org/10.1007/s10040-002-024...
; Fredlund & Rahardjo, 1993Fredlund, D.G., & Rahardjo, H. (1993). Soil mechanics for unsaturated soils. John Wiley & Sons.; Ng et al., 2001Ng, C.W., Zhang, L.M., & Ho, K.K. (2001). Influence of laterally loaded sleeved piles and pile groups on slope stability. Canadian Geotechnical Journal, 38(3), 553-566. http://dx.doi.org/10.1139/t00-109.
http://dx.doi.org/10.1139/t00-109...
).

It is necessary to determine unsaturated soil properties for a proper unsaturated slope stability analysis, such as unsaturated shear strength, hydraulic conductivity, and soil water retention curve (SWRC). So, suction-controlled, and conventional triaxial tests as well as tests to determine the soil water characteristic curves (SWRC) (i.e., suction-plate, pressure-chamber, and filter-paper) should be carried out on undisturbed soil samples. It is important to recognize the significant progress that suction control tests have brought to the understanding of unsaturated soils, especially Hilf's (1956)Hilf, J.W. (1956). An investigation of pore-water pressure in compacted cohesive soils. Faculty of the Graduate School, University of Colorado. axis translation technique. It is important to mention, however, that suction control techniques are complex to use and expensive, which brings the need to search for more cost-effective alternatives.

In the slope stability analysis, it is necessary to identify the geomaterials in the site profile, their thicknesses, the ground water level position, as well as the strength parameters of each of these materials. Geotechnical site characterization can be carried out by means of in situ and laboratory tests. In situ tests, such as the Flat Dilatometer (DMT) can be used as an alternative to the traditional drilling, sampling, and laboratory testing, especially in sandy soils, where sampling is complex and does not always maintain the natural characteristics of the soil. Combining the use of a in situ profiling test with a specific determination of some indexes or properties is a modern approach to better characterize the site.

This paper presents and discusses the use of the Flat Dilatometer Test (DMT) in the geotechnical site characterization of an unsaturated tropical sandy soil profile, with emphasis on its application to the slope stability analysis. The classification of the soil profile and the estimative of strength parameters were done and compared with laboratory test data. In addition, the stability analysis on a cut slope in unsaturated sandy soil profile is presented and discussed.

2. Flat dilatometer test (DMT)

The Flat Dilatometer Test (DMT) is a simple and repeatable in situ test, and its data can be used to identify the soil type and estimate design parameters such as in situ coefficient of lateral earth pressure (K0), undrained shear strength (Su), friction angle (ϕ), overconsolidation ratio (OCR), and vertical drained constrained modulus (MDMT) (Marchetti, 1980Marchetti, S. (1980). In situ tests by flat dilatometer. Journal of the Geotechnical Engineering Division, 106(3), 299-321. http://dx.doi.org/10.1061/AJGEB6.0000934.
http://dx.doi.org/10.1061/AJGEB6.0000934...
). The test can also be used to predict settlements of shallow foundations, monitoring densification (i.e., soil compaction), liquefaction analysis as well as in the detecting slip surfaces in overconsolidated clays slopes. Moreover, dissipation tests can be performed to estimate the in situ consolidation (coefficient of consolidation, ch) and flow parameters (coefficient of permeability, kh) (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.; Marchetti & Monaco, 2018Marchetti, S., & Monaco, P. (2018). Recent improvements in the use, interpretation, and applications of DMT and SDMT in practice. Geotechnical Testing Journal, 41(5), 20170386. http://dx.doi.org/10.1520/GTJ20170386.
http://dx.doi.org/10.1520/GTJ20170386...
; Schmertmann, 1986Schmertmann, J.H. (1986). Suggested method for performing the flat dilatometer test. Geotechnical Testing Journal, 9(2), 93. http://dx.doi.org/10.1520/GTJ11035J.
http://dx.doi.org/10.1520/GTJ11035J...
).

The DMT has a flat stainless-steel blade with a 60 mm diameter steel membrane mounted flush on one side (Marchetti, 1980Marchetti, S. (1980). In situ tests by flat dilatometer. Journal of the Geotechnical Engineering Division, 106(3), 299-321. http://dx.doi.org/10.1061/AJGEB6.0000934.
http://dx.doi.org/10.1061/AJGEB6.0000934...
). The nominal dimensions of the blade are 95 mm width and 15 mm thickness, having a cutting edge angled between 24° and 32° to penetrate the soil (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.). The blade is designed to safely withstand up to 250 kN pushing force. The other components are the control unit, pneumatic-electrical cable, and gas tank (Figure 1).

Figure 1
Schematic representation of DMT [adapted from Marchetti et al. (2001)Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.].

The blade is advanced quasi statically or dynamically into the ground. The penetration rate is usually equal to 20 mm/s as in the cone penetration tests. The calibration procedure to obtain ΔA and ΔB pressures, necessary to overcome membrane stiffness, must be done before each profile. A-Pressure (just begin to move the membrane into the soil) and B-Pressure (move the membrane center 1.1 mm against the soil) are typically recorded every 200 mm during the test and p0 and p1 pressures were calculated. C-Pressure can also be recorded. The DMT interpretation starts with the intermediate parameters (ID, KD and ED) determination (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.):

I D = p 1 p 0 p 0 u 0 (1)
K D = p 0 u 0 σ ' v (2)
E D = 34.7. p 1 p 0 (3)

where u0 is the hydrostatic pressure, and the σ’v is the effective vertical stress.

ID is calculated to identify soil type. In general, ID provides an expressive profile of soil type and, in “normal” soils a reasonable soil description (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.). The KD provides the basis for several soil parameters correlations and is the key result of the dilatometer test (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.). ED in general should not be used as such, especially because it lacks information on stress history (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.).

3. Study site

The study site is the experimental research area at the São Paulo State University (Unesp) - Bauru Campus and a typical cut slope (the case study) that occurs in this region. The schematic representation of the study site in the city of Bauru, state of São Paulo, Brazil, is shown in Figure 2, as well as the location of the cut slope (case study).

Figure 2
Schematic representation of the Unesp research site and the cut slope in the city of Bauru, São Paulo State, Brazil.

Several site characterization programs including SPT, SPT-T, S-SPT, DMT, PMT, SCPT, CH and DH tests were carried out at the study site. Sample pits were also excavated to retrieve undisturbed and disturbed soil blocks at 1.5, 3.0, 5.0, 7.0, and 9.0 m depth. Soil specimens from these blocks were tested in laboratory for soil characterization and determination of mechanical properties and parameters.

The site profile is a red clayey fine to medium sand identified based on SPT. MCT Classification System (Nogami & Villibor, 1981Nogami, J.S., & Villibor, D.F. (1981). Uma nova classificação de solos para finalidades rodoviárias. In Anais do Simpósio Brasileiro de Solos Tropicais em Engenharia (pp. 30-41). Rio de Janeiro. ABMS.) classified the top 13 m as lateritic soil behavior (LA’) followed by a non-lateritic soil behavior (NA’), as discussed Giacheti et al. (1998)Giacheti, H.L., Coelho, V., & Carvalho, D. (1998). Caracterização geotécnica de dois perfis de solos tropicais com base em ensaios de laboratório. In Anais do XI Congresso Brasileiro de Mecânica dos Solos e Engenharia de Fundações – COBRAMSEG (pp. 195-202). São Paulo: ABMS.. The soil is classified in the Unified Soil Classification System as a SM Group soil. Some relevant geotechnical characteristics are summarized in Figure 3.

Figure 3
In situ and laboratory tests data for the experimental research site. The water level is deeper than 30 mat this location [adapted from Rocha & Giacheti (2018)Rocha, B.P., & Giacheti, H.L. (2018). Site characterization of a tropical soil by in situ tests. Dyna, 85(206), 211-219. http://dx.doi.org/10.15446/dyna.v85n206.67891.
http://dx.doi.org/10.15446/dyna.v85n206....
].

Figure 3a shows the typical site profile defined by SPT tests. Figure 3b shows that SPT N values increase almost linearly with depth. From CPT tests, the corrected cone resistance (qt) and the sleeve friction (fs) shows higher values at the top 1 m and tends to increase with depth leading to a friction ratio (Rf = (fs/qt)*100%) between 1.0 and 3.0%. Moreover, Figure 3 presents the variation of shear wave velocity (Vs) with depth (Figure 3e) measured by cross-hole (CH), SCPT and down-hole tests (DH 1 and DH2). The void ratio (e) at 1.0 m depth is equal to 0.72 and drops to about 0.60 at 16 m depth. The dry unit weight (γd) at 1.0 m depth is equal to 15.64 kN/m3 and increase with depth (Figure 3f). The groundwater level was deeper than 30 m depth at this location.

4. Results and analysis

Six DMTs were carried out at the site and Figure 4 presents the DMT profiles in terms of ID, KD and ED data plotted as an average profile as well as plus and minus one standard deviation (SD). The ID, KD and ED parameters were calculated using Equations 1, 2, and 3, respectively (Marchetti, 1980Marchetti, S. (1980). In situ tests by flat dilatometer. Journal of the Geotechnical Engineering Division, 106(3), 299-321. http://dx.doi.org/10.1061/AJGEB6.0000934.
http://dx.doi.org/10.1061/AJGEB6.0000934...
).

Figure 4
DMT data and interpretation for the study site. a) Corrected first reading, b) Corrected second reading, c) Material index, d) Horizontal stress index, e) Dilatometer modulus.

From the interpretation of the laboratory and in situ tests (Figure 3) it is observed that the site profile is relatively homogeneous, and that possible variations observed in the mechanical parameters of the soil may be associated with soil suction, as well as its variation over time. In this sense, the interpretation of in situ tests, as well as the design of geotechnical works (e.g., slope stability) should consider the effect of soil suction.

Since the DMT does not provide soil sample, the soil type can be defined based on the material index (ID) and by means of the chart proposed by Marchetti & Crapps (1981)Marchetti, S., & Crapps, D. K. (1981). Flat dilatometer manual. Internal Report of G.P.E. Inc. relating ID and ED . Geotechnical parameters can be estimated by using classical correlations and these parameters can be compared with laboratory-determined reference values obtained by triaxial tests on undisturbed soil samples.

4.1 Soil classification

The first step of the site characterization in the slope stability analysis consists of defining the representative site profile. The definition of the soil type was done by means of grain size distribution tests in the laboratory. Tropical soils normally present a macroaggregate structure, due to the cementing action of iron and aluminum oxides and hydroxides. It is recommended to carry out grain size distribution tests with and without the use of dispersant, to have a textural classification of the soil in these two conditions to better understand how this aggregation occurs.

Figure 5 shows the particle size distribution curves with and without the use of dispersant for samples collected at 0.25 m, 3.50 m, 5.50 m, and 8.00 m depth using a helical auger. It can be seen in this figure a significant difference between the clay and silt fractions for these two conditions, with implications for the soil classification. For instance, the sample collected at 8 m depth is classified as a clayey sand (78% Sand; 5% Silt; 17% Clay) with the use of dispersant and as a silty sand (65% Sand; 35% Silt) without the use of dispersant. It is due to the action of iron and aluminum oxides and hydroxides, present in tropical soils, which produce agglutination of the clay fraction, leading to a macroaggregate structure.

Figure 5
Grain size distribution curves with and without dispersant for the soil samples collected at 0.25, 3.50, 5.50, and 8.00 m depth [adapted from Rocha (2018)Rocha, B.P. (2018). Geotechnical Characterization of Unsaturated Tropical Soils by In situ tests [Doctoral thesis, University of São Paulo]. University of São Paulo’s repository (in Portuguese). https://doi.org/10.11606/T.18.2018.tde-03122018-103909
https://doi.org/10.11606/T.18.2018.tde-0...
].

The soil classification by means of the DMT uses the Material Index (ID) and the Dilatometric Modulus (ED) as proposed by Marchetti & Crapps (1981)Marchetti, S., & Crapps, D. K. (1981). Flat dilatometer manual. Internal Report of G.P.E. Inc.. The site profile is classified as a sandy silt by the DMT data, as shown in Figure 4 and Figure 6. This classification is not in agreement with the grain size of the soil obtained by the test with dispersant (Figure 5), which classifies the soil as a fine sand with little clay. However, the soil is classified as a silty sand when the grain size distribution curve is determined without the use of the dispersant. That classification is more consistent with the behavior identified by the DMT. It is important to note that the ID index is a parameter that reflects the mechanical behavior of the soil, as a stiffness index, and not the grain size composition of the soil (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.). Furthermore, the ID index indicates that mixtures of sand and clay are usually classified as silts, as observed for the study site (Marchetti et al., 2001Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.).

Figure 6
Soil classification for the soils from the study site in the DMT chart [adapted from Marchetti & Crapps (1981)Marchetti, S., & Crapps, D. K. (1981). Flat dilatometer manual. Internal Report of G.P.E. Inc.].

4.2 Geomechanical parameters

The geomechanical parameters estimative by means of the average profiles ± standard deviation of ID, KD and ED was done considering the six DMTs. It is important to emphasize that the investigated profile is unsaturated. Thus, the mechanical soil parameters were determined by triaxial tests both in saturated and unsaturated condition for a soil suction equal to 50 kPa (Fernandes, 2022Fernandes, J.B. (2022). Strength and deformability study of an unsaturated tropical soil profile [Doctoral thesis, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/217501
https://repositorio.unesp.br/handle/1144...
; Saab, 2016Saab, A.L. (2016). Influence of suction on the compressibility parameters of an unsaturated soil [Master's dissertation, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/bitstream/handle/11449/148696/saab_al_me_bauru.pdf?sequence=3
https://repositorio.unesp.br/bitstream/h...
) through the axis translation technique (Hilf, 1956Hilf, J.W. (1956). An investigation of pore-water pressure in compacted cohesive soils. Faculty of the Graduate School, University of Colorado.). This value was used in laboratory to represent the field condition as it was measured in situ by tensiometers and granular matrix sensors.

Figure 7 shows the profiles of total unit weight (γn), friction angle (ϕ), and cohesion intercept (c) estimated by the DMT. These values were compared with reference values determined in laboratory from undisturbed samples (Fagundes & Rodrigues, 2015Fagundes, L.S., & Rodrigues, R.A. (2015). Shear strength of a natural and compacted tropical soil. The Electronic Journal of Geotechnical Engineering, 20, 47-58.; Fernandes, 2022Fernandes, J.B. (2022). Strength and deformability study of an unsaturated tropical soil profile [Doctoral thesis, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/217501
https://repositorio.unesp.br/handle/1144...
). The γn values were estimated by Schmertmann (1986)Schmertmann, J.H. (1986). Suggested method for performing the flat dilatometer test. Geotechnical Testing Journal, 9(2), 93. http://dx.doi.org/10.1520/GTJ11035J.
http://dx.doi.org/10.1520/GTJ11035J...
’s equation (Equation 4) while friction angle was estimated by equation suggested by Marchetti (1997)Marchetti, S. (1997). The flat dilatometer: design applications. In Proceedings of the. Third International Geotechnical Engineering Conference - ICSMGE (pp. 421-448). Cairo.. The cohesion intercept values were calculated from Equation 6 proposed by Cruz & Viana Da Fonseca (2006)Cruz, N., & Viana Da Fonseca, A. (2006). Portuguese experience in residual soil characterization by DMT tests. In Proceedings of the Second International Conference on the Flat Dilatometer – DMT 2006 (pp. 359-364). Washington D.C.. It is worth noting that the proposed equation was developed for a residual silty sand from Granite.

γ n γ w = 1.12 . E D p a 0.1 . I D 0.05 (4)
ϕ = 28 + 14.6 . log K D 2.1 . log 2 K D (5)
c = 0 .376 . v O C R + 3 .08 (6)

γw is the specific weight of water, pa is the atmospheric reference pressure (100 kPa), and vOCR is the virtual overconsolidation ratio obtained from Marchetti & Crapps (1981)Marchetti, S., & Crapps, D. K. (1981). Flat dilatometer manual. Internal Report of G.P.E. Inc. approach.

Figure 7
Estimated parameters from DMT for the study site and reference values determined by triaxial tests.

Figure 7a shows that the density values estimated by Equation 4 are in good agreement with those determined from undisturbed soil samples, especially between 1 and 11 m depth. The ϕ value varied from 27.0º for 1 m depth to 31.2º for 9 m depth at the saturated condition and varied from 30.0º for 1 m depth to 29.5º for 9 m depth at the unsaturated condition (Figure 7b and Table 1). So, an average ϕ value equal to 30.9º in the saturated condition, and 31.2º in the unsaturated condition were calculated from the triaxial tests. The ϕ values determined from the average profiles of ID, KD, and ED agree relatively well with the reference values from triaxial tests below 2 m depth. It did not occur above 2 m depth, and it could be related to the effect of high soil suction values, which increased p0 and p1 values, mainly in the topsoil due to the soil-atmosphere interaction, which reflect in the estimative of geomechanical parameters by DMT (Giacheti et al., 2019Giacheti, H.L., Bezerra, R.C., Rocha, B.P., & Rodrigues, R.A. (2019). Seasonal influence on cone penetration test: an unsaturated soil site example. Journal of Rock Mechanics and Geotechnical Engineering, 11(2), 361-368. http://dx.doi.org/10.1016/j.jrmge.2018.10.005.
http://dx.doi.org/10.1016/j.jrmge.2018.1...
; Rocha, 2018Rocha, B.P. (2018). Geotechnical Characterization of Unsaturated Tropical Soils by In situ tests [Doctoral thesis, University of São Paulo]. University of São Paulo’s repository (in Portuguese). https://doi.org/10.11606/T.18.2018.tde-03122018-103909
https://doi.org/10.11606/T.18.2018.tde-0...
; Rocha et al., 2021Rocha, B.P., Rodrigues, R.A., & Giacheti, H.L. (2021). The flat dilatometer test in an unsaturated tropical soil site. Geotechnical and Geological Engineering, 39(8), 5957-5969. http://dx.doi.org/10.1007/s10706-021-01849-1.
http://dx.doi.org/10.1007/s10706-021-018...
).

Table 1
Strength parameters obtained by triaxial tests for different test depth for saturated and unsaturated conditions (Fernandes, 2022Fernandes, J.B. (2022). Strength and deformability study of an unsaturated tropical soil profile [Doctoral thesis, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/217501
https://repositorio.unesp.br/handle/1144...
; Fernandes et al., 2022Fernandes, J.B., Saab, A.L., Rocha, B.P., Rodrigues, R.A., Lodi, P.C., & Giacheti, H.L. (2022). Geomechanical parameters in the active zone of an unsaturated tropical soil site via laboratory tests. Soils and Rocks, 45(4), 1-15. http://dx.doi.org/10.28927/SR.2022.000422.
http://dx.doi.org/10.28927/SR.2022.00042...
).

Cohesion intercept values estimated from the DMT were compared with those determined in the laboratory by triaxial tests on undisturbed soil samples (Figure 7c). It should be mentioned that a soil suction value equal to 20 kPa to estimate the cohesion intercept from the DMT. This suction value was assumed from moisture content profiles and soil water retention curves available for the site and will be discussed latter on. The cohesion intercept values increased with depth and with soil suction (Table 1). The cohesion value was equal to 0 kPa in the saturated condition, and 3.0 kPa for 50 kPa suction for the sample collected at 1.5 m depth. The cohesion intercept was equal to 4.5 kPa in the saturated condition, and 22 kPa for 50 kPa suction for the sample collected at 9.0 m depth.

It is important to mention that the cohesion intercept for tropical sandy soils is mainly related to the cementing structures and the effect of the unsaturated condition and this parameter varies seasonally depending on the rainfall regime and water infiltration in the soil mass.

5. Slope stability analysis

The stability analysis of a cut slope in the city of Bauru, São Paulo state, which has a site profile and mechanical parameters like those found in the experimental research site at the Unesp, will be performed. DMT (Figure 4) and the triaxial tests (Table 1) previously presented will be used in these analyses.

The slope is 6.6 m high and has an approximate slope gradient of around 55°. The slope geometry was defined from planialtimetric surveys, and the groundwater was assumed at the slope foot, based on SPT data carried out at this location. No groundwater flow was considered in the stability analysis. Figure 8 shows the defined geometry and the assumed groundwater position for the slope stability analysis.

Figure 8
Slope geometry.

The GeoStudio software was used to perform the slope stability analysis using the Simplified Bishop's stability method. The shear strength criterion is that one proposed by Fredlund et al. (1978)Fredlund, D.G., Morgenstern, N.R., & Widger, R.A. (1978). The shear strength of unsaturated soils. Canadian Geotechnical Journal, 15(3), 313-321. http://dx.doi.org/10.1139/t78-029.
http://dx.doi.org/10.1139/t78-029...
for the unsaturated soil.

5.1 The analysis with DMT data

At the time the DMTs were carried out, disturbed soil samples were taken using a helical auger to determine the moisture content profiles (Figure 9a) to estimate the soil suction. The moisture content profile ranges from 4.0 to 11.2%, and the estimated soil suction varies from 4.2 to 50 kPa, with an average value equal to 20 kPa. These values were estimated from the soil water retention curves (SWRCs) determined from undisturbed samples collected at 1, 3, 5, and 7 m depth (Figure 9b) by Fernandes (2022)Fernandes, J.B. (2022). Strength and deformability study of an unsaturated tropical soil profile [Doctoral thesis, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/217501
https://repositorio.unesp.br/handle/1144...
. Table 2 presents the input parameters for the stability analyses by DMT and triaxial tests. In these analyses, the suction was assumed not to affect the internal friction angle (ϕ) of the soil.

Figure 9
Moisture content profiles (a) and the soil water retention curves (b) [adapted form Fernandes (2022)Fernandes, J.B. (2022). Strength and deformability study of an unsaturated tropical soil profile [Doctoral thesis, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/217501
https://repositorio.unesp.br/handle/1144...
].
Table 2
Average soil parameters assumed for slope analysis stability by DMT and triaxial tests for different soil suctions.

Figure 10 presents the safety factor (FS) equal to 0.943, as well as the possible rupture surfaces considering the soil parameters obtained by DMT. This value is lower than one, i.e., the slope is in an unstable condition for the assumed boundary conditions, stresses, and soil parameters. It is worth noting that the geomechanical parameters (c and ϕ) were obtained using empirical correlations, and they may not adequately represent the mechanical behavior of the studied unsaturated tropical soil.

Figure 10
Minimum safety factor (FS) and slip surfaces tested assuming shear strength parameters by DMT.

5.2 The analysis with triaxial test data

To compare the FS values calculated from the interpretation of DMT data, the cohesion intercept value determined using the triaxial test data for the same field suction value (i.e., 20 kPa) was estimated. For this, the model of Fredlund et al. (1978)Fredlund, D.G., Morgenstern, N.R., & Widger, R.A. (1978). The shear strength of unsaturated soils. Canadian Geotechnical Journal, 15(3), 313-321. http://dx.doi.org/10.1139/t78-029.
http://dx.doi.org/10.1139/t78-029...
was used to estimate the unsaturated cohesion for 20 kPa of suction (Equation 7), which resulted in a cohesion intercept value equal to 4.7 kPa.

In order to apply Equation 7 the values of cohesion intercept versus suction (Table 1) were reinterpreted to define ϕb that resulted in the value of 5.10o (tg ϕb = 0.089). Hence, the calculated FS was equal to 1.035, close to that found from the interpretation of the DMT data (Figure 11).

c = c s a t + u a u w . tan ϕ b (7)

where csat is the saturated cohesion, uauw is the soil suction and ϕb is the angle indicating the increase in shear strength due to soil suction.

Figure 11
Minimum safety factor and slip surfaces tested assuming shear strength parameters from the unsaturated triaxial tests with a suction value equal to 20 kPa.

Following the same approach, another stability analysis was performed considering an average cohesion value obtained for the suction condition equal to 50 kPa (Table 1). In this case, the calculated FS was 2.161 (Figure 12).

Figure 12
Minimum safety factor and slip surfaces tested assuming shear strength parameters from the unsaturated triaxial tests with a suction value equal to 50 kPa.

The results of this and other similar studies have led to a general consensus that unsaturated condition, with its associated soil suction, tend to increase slope stability due to increase in apparent cohesion (Dai et al., 2022Dai, G., Zhang, F., & Wang, Y. (2022). Stability analysis of layered slopes in unsaturated soils. Frontiers of Structural and Civil Engineering, 16(3), 378-387. http://dx.doi.org/10.1007/s11709-022-0808-2.
http://dx.doi.org/10.1007/s11709-022-080...
; Kang et al., 2020Kang, S., Lee, S.-R., & Cho, S.-E. (2020). Slope stability analysis of unsaturated soil slopes based on the site-specific characteristics: a case study of Hwangryeong Mountain, Busan, Korea. Sustainability, 12(7), 2839. http://dx.doi.org/10.3390/su12072839.
http://dx.doi.org/10.3390/su12072839...
; Sivakumar Babu & Murthy, 2005Sivakumar Babu, G.L., & Murthy, D.S. (2005). Reliability analysis of unsaturated soil slopes. Journal of Geotechnical and Geoenvironmental Engineering, 131(11), 1423-1428. http://dx.doi.org/10.1061/(ASCE)1090-0241(2005)131:11(1423).
http://dx.doi.org/10.1061/(ASCE)1090-024...
).

5.3 Slope failure back analysis

Vieira (2021)Vieira, R.H.M. (2021). Back-analysis of the slide of a planar slope of an unsaturated tropical sandy soil from Bauru-SP. [Master's dissertation, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/215781
https://repositorio.unesp.br/handle/1144...
performed the back analysis of similar shear slope slides by integrating data from laboratory tests and numerical modeling. The author analyzed a rupture that occurred after a 124.7 mm rainfall. Figure 13 shows the slope after rupture. The GeoStudio software (SLOPE/W and SEEP/W) was used to perform the numerical simulations in a time-dependent analysis. The approach adopted considered the Simplified Bishop's stability method, as well as the Mohr-Coulomb criterion that incorporates the effect of the saturated condition, as proposed by Fredlund et al. (1978)Fredlund, D.G., Morgenstern, N.R., & Widger, R.A. (1978). The shear strength of unsaturated soils. Canadian Geotechnical Journal, 15(3), 313-321. http://dx.doi.org/10.1139/t78-029.
http://dx.doi.org/10.1139/t78-029...
. The slope geometry and the failure surface were obtained from planialtimetric surveys. The rainfall history of the region during the slope failure period was obtained from a meteorological station near the slope. The slope failure back analysis was performed by combining the cohesion intercept and friction angle values to obtain them, considering the unitary safety factor (FS = 1) at the time of the 124.7 mm rain that caused the slope failure.

Figure 13
The slope after failure [adapted from Vieira (2021)Vieira, R.H.M. (2021). Back-analysis of the slide of a planar slope of an unsaturated tropical sandy soil from Bauru-SP. [Master's dissertation, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/215781
https://repositorio.unesp.br/handle/1144...
].

It is worth noting that slope stability analyses that use geomechanical parameters determined by triaxial compression tests with controlled suction in numerical simulations that incorporate the hydraulic properties of the soil, and flow conditions to simulate rainfall and its interaction with the soil mass are unusual in the practice of geotechnical jobs.

6. Conclusion

The main conclusions are:

  • The DMT allowed to define the site profile and to give a preliminary estimative of the strength parameters of the soils;

  • The natural specific weight profile (γn), friction angle (ϕ) and the cohesion intercept (c) presented reasonable agreement with the reference values;

  • The DMT was an interesting tool in defining the site profile for the case study as well as for estimating the strength parameters required in the slope stability analyses in the preliminary design phase;

  • The slope stability analyses performed considering the geometry, stresses and parameters estimated by the DMT and by the triaxial tests for a suction value of 20 kPa indicated that the slope is at the eminence of rupture (FS ≈ 1). In addition, the slope was found to be stable for a suction value of 50 kPa based on the triaxial test data, since the factor of safety in this case is greater than 2.0 (FS ≈ 2.2);

  • The back analysis of a slope with similar characteristics to those in this study confirms that the failure occurs due to the reduced contribution to soil shear strength by the effect of rainfall and its interaction with the soil mass;

  • The DMT can be used as logging test in the preliminary site characterization of unsaturated tropical soil sites specially to define the stratigraphical profile, site variability, select the regions to collect disturbed and undisturbed soil samples and as the first attempt to estimate the geotechnical design parameters via correlations;

  • The ideal condition is to characterize the site using a logging test like the SPT, CPT or DMT and refine it with triaxial compression tests with controlled suction to determine the shear strength parameters to have a better understanding of the suction effect during the lifetime of the slope considering the hydraulic properties of the soil and the influence of rainfall regime, infiltration and drainage conditions.

List of symbols

c Cohesion

ch Coefficient of consolidation

csat Saturated cohesion

e Void ratio

fs Sleeve friction

pa Atmospheric reference pressure (100 kPa)

p0 Corrected first reading

p1 Corrected second reading

qt Cone resistance

s Soil suction

A- Pressure required to just begin to move the membrane into the soil

B- Pressure required to move the membrane center 1.1 mm against the soil

C- Pressure taken by slowly deflating the membrane immediately after B-pressure is reached

CH Crosshole test

DH Downhole test

DMT Flat dilatometer tests

ED Dilatometer modulus

FS Safety factor

ID Material index

KD Horizontal stress index

Kh Coefficient of permeability

K0 In situ coefficient of lateral earth pressure

LA’ Lateritic soil behavior

MDMT Vertical drained constrained modulus

NA’ Non-lateritic soil behavior

OCR Overconsolidation ratio

PMT Pressuremeter tests

Rf Friction ratio

SCPT Seismic cone test

SM Silty sand

SPT Standard penetration tests

SPT-T Standard penetration tests with torque measurements

S-SPT Seismic SPT

SWRC Soil water retention curve

Su Undrained shear strength

Vs Shear wave velocity

ϕ Friction angle

ϕb Angle indicating the increase in shear strength due to soil suction

γd Dry unit weight

γn Total unit weight

γw Specific weight of water

u0 Hydrostatic pressure

ua – uw Soil suction

vOCR Virtual overconsolidation ratio

σ’v Effective vertical stress

ΔA Calibration parameters used to correct the A, B pressure

ΔB Calibration parameters used to correct the A, B pressure

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgements

The authors thank the São Paulo Research Foundation - FAPESP (Grants 2015/17260-0 and 2017/23174-5), the National Council for Scientific and Technological Development - CNPq (Grant 436478/2018-8), and the Coordination for the Improvement of Higher Education Personnel – CAPES, for supporting their research.

  • Discussion open until February 28, 2024.

References

  • Alonso, E.E., Gens, A., & Delahaye, C.H. (2003). Influence of rainfall on the deformation and stability of a slope in overconsolidated clays: a case study. Hydrogeology Journal, 11(1), 174-192. http://dx.doi.org/10.1007/s10040-002-0245-1
    » http://dx.doi.org/10.1007/s10040-002-0245-1
  • Alonso, E.E., Gens, A., & Josa, A. (1990). A constitutive model for partially saturated soils. Geotechnique, 40(3), 405-430. http://dx.doi.org/10.1680/geot.1990.40.3.405
    » http://dx.doi.org/10.1680/geot.1990.40.3.405
  • Blight, G.E. (2003). The vadose zone soil-water balance and transpiration rates of vegetation. Geotechnique, 53(1), 55-64. http://dx.doi.org/10.1680/geot.2003.53.1.55
    » http://dx.doi.org/10.1680/geot.2003.53.1.55
  • Cho, S.E., & Lee, S.R. (2002). Evaluation of surficial stability for homogeneous slopes considering rainfall characteristics. Journal of Geotechnical and Geoenvironmental Engineering, 128(9), 756-763. http://dx.doi.org/10.1061/(ASCE)1090-0241(2002)128:9(756)
    » http://dx.doi.org/10.1061/(ASCE)1090-0241(2002)128:9(756)
  • Cruz, N., & Viana Da Fonseca, A. (2006). Portuguese experience in residual soil characterization by DMT tests. In Proceedings of the Second International Conference on the Flat Dilatometer – DMT 2006 (pp. 359-364). Washington D.C.
  • Cui, Y.-J., Tang, A.-M., Mantho, A.T., & De Laure, E. (2008). Monitoring field soil suction using a miniature tensiometer. Geotechnical Testing Journal, 31(1), 100769. http://dx.doi.org/10.1520/GTJ100769
    » http://dx.doi.org/10.1520/GTJ100769
  • Dai, G., Zhang, F., & Wang, Y. (2022). Stability analysis of layered slopes in unsaturated soils. Frontiers of Structural and Civil Engineering, 16(3), 378-387. http://dx.doi.org/10.1007/s11709-022-0808-2
    » http://dx.doi.org/10.1007/s11709-022-0808-2
  • Fagundes, L.S., & Rodrigues, R.A. (2015). Shear strength of a natural and compacted tropical soil. The Electronic Journal of Geotechnical Engineering, 20, 47-58.
  • Fernandes, J.B. (2022). Strength and deformability study of an unsaturated tropical soil profile [Doctoral thesis, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/217501
    » https://repositorio.unesp.br/handle/11449/217501
  • Fernandes, J.B., Saab, A.L., Rocha, B.P., Rodrigues, R.A., Lodi, P.C., & Giacheti, H.L. (2022). Geomechanical parameters in the active zone of an unsaturated tropical soil site via laboratory tests. Soils and Rocks, 45(4), 1-15. http://dx.doi.org/10.28927/SR.2022.000422
    » http://dx.doi.org/10.28927/SR.2022.000422
  • Fredlund, D.G., & Rahardjo, H. (1993). Soil mechanics for unsaturated soils John Wiley & Sons.
  • Fredlund, D.G., Morgenstern, N.R., & Widger, R.A. (1978). The shear strength of unsaturated soils. Canadian Geotechnical Journal, 15(3), 313-321. http://dx.doi.org/10.1139/t78-029
    » http://dx.doi.org/10.1139/t78-029
  • Giacheti, H.L., Bezerra, R.C., Rocha, B.P., & Rodrigues, R.A. (2019). Seasonal influence on cone penetration test: an unsaturated soil site example. Journal of Rock Mechanics and Geotechnical Engineering, 11(2), 361-368. http://dx.doi.org/10.1016/j.jrmge.2018.10.005
    » http://dx.doi.org/10.1016/j.jrmge.2018.10.005
  • Giacheti, H.L., Coelho, V., & Carvalho, D. (1998). Caracterização geotécnica de dois perfis de solos tropicais com base em ensaios de laboratório. In Anais do XI Congresso Brasileiro de Mecânica dos Solos e Engenharia de Fundações – COBRAMSEG (pp. 195-202). São Paulo: ABMS.
  • Hilf, J.W. (1956). An investigation of pore-water pressure in compacted cohesive soils Faculty of the Graduate School, University of Colorado.
  • Kang, S., Lee, S.-R., & Cho, S.-E. (2020). Slope stability analysis of unsaturated soil slopes based on the site-specific characteristics: a case study of Hwangryeong Mountain, Busan, Korea. Sustainability, 12(7), 2839. http://dx.doi.org/10.3390/su12072839
    » http://dx.doi.org/10.3390/su12072839
  • Marchetti, S. (1980). In situ tests by flat dilatometer. Journal of the Geotechnical Engineering Division, 106(3), 299-321. http://dx.doi.org/10.1061/AJGEB6.0000934
    » http://dx.doi.org/10.1061/AJGEB6.0000934
  • Marchetti, S. (1997). The flat dilatometer: design applications. In Proceedings of the. Third International Geotechnical Engineering Conference - ICSMGE (pp. 421-448). Cairo.
  • Marchetti, S., & Crapps, D. K. (1981). Flat dilatometer manual Internal Report of G.P.E. Inc.
  • Marchetti, S., & Monaco, P. (2018). Recent improvements in the use, interpretation, and applications of DMT and SDMT in practice. Geotechnical Testing Journal, 41(5), 20170386. http://dx.doi.org/10.1520/GTJ20170386
    » http://dx.doi.org/10.1520/GTJ20170386
  • Marchetti, S., Monaco, P., Totani, G., & Calabrese, M. (2001). The Flat Dilatometer Test (DMT) in soil investigations. In Proceedings of the International Conference on in situ Measurement of Soil Properties – IN SITU 2001 (41 p.). Indonesia.
  • Ng, C.W., & Shi, Q. (1998). A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage. Computers and Geotechnics, 22(1), 1-28. http://dx.doi.org/10.1016/S0266-352X(97)00036-0
    » http://dx.doi.org/10.1016/S0266-352X(97)00036-0
  • Ng, C.W., Zhang, L.M., & Ho, K.K. (2001). Influence of laterally loaded sleeved piles and pile groups on slope stability. Canadian Geotechnical Journal, 38(3), 553-566. http://dx.doi.org/10.1139/t00-109
    » http://dx.doi.org/10.1139/t00-109
  • Nogami, J.S., & Villibor, D.F. (1981). Uma nova classificação de solos para finalidades rodoviárias. In Anais do Simpósio Brasileiro de Solos Tropicais em Engenharia (pp. 30-41). Rio de Janeiro. ABMS.
  • Rahardjo, H., Kim, Y., & Satyanaga, A. (2019). Role of unsaturated soil mechanics in geotechnical engineering. International Journal of Geo-Engineering, 10(1), 8. http://dx.doi.org/10.1186/s40703-019-0104-8
    » http://dx.doi.org/10.1186/s40703-019-0104-8
  • Rocha, B.P. (2018). Geotechnical Characterization of Unsaturated Tropical Soils by In situ tests [Doctoral thesis, University of São Paulo]. University of São Paulo’s repository (in Portuguese). https://doi.org/10.11606/T.18.2018.tde-03122018-103909
    » https://doi.org/10.11606/T.18.2018.tde-03122018-103909
  • Rocha, B.P., & Giacheti, H.L. (2018). Site characterization of a tropical soil by in situ tests. Dyna, 85(206), 211-219. http://dx.doi.org/10.15446/dyna.v85n206.67891
    » http://dx.doi.org/10.15446/dyna.v85n206.67891
  • Rocha, B.P., Rodrigues, R.A., & Giacheti, H.L. (2021). The flat dilatometer test in an unsaturated tropical soil site. Geotechnical and Geological Engineering, 39(8), 5957-5969. http://dx.doi.org/10.1007/s10706-021-01849-1
    » http://dx.doi.org/10.1007/s10706-021-01849-1
  • Saab, A.L. (2016). Influence of suction on the compressibility parameters of an unsaturated soil [Master's dissertation, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/bitstream/handle/11449/148696/saab_al_me_bauru.pdf?sequence=3
    » https://repositorio.unesp.br/bitstream/handle/11449/148696/saab_al_me_bauru.pdf?sequence=3
  • Schmertmann, J.H. (1986). Suggested method for performing the flat dilatometer test. Geotechnical Testing Journal, 9(2), 93. http://dx.doi.org/10.1520/GTJ11035J
    » http://dx.doi.org/10.1520/GTJ11035J
  • Sivakumar Babu, G.L., & Murthy, D.S. (2005). Reliability analysis of unsaturated soil slopes. Journal of Geotechnical and Geoenvironmental Engineering, 131(11), 1423-1428. http://dx.doi.org/10.1061/(ASCE)1090-0241(2005)131:11(1423)
    » http://dx.doi.org/10.1061/(ASCE)1090-0241(2005)131:11(1423)
  • Tohari, A., Nishigaki, M., & Komatsu, M. (2007). Laboratory rainfall-induced slope failure with moisture content measurement. Journal of Geotechnical and Geoenvironmental Engineering, 133(5), 575-587. http://dx.doi.org/10.1061/(ASCE)1090-0241(2007)133:5(575)
    » http://dx.doi.org/10.1061/(ASCE)1090-0241(2007)133:5(575)
  • Vieira, R.H.M. (2021). Back-analysis of the slide of a planar slope of an unsaturated tropical sandy soil from Bauru-SP [Master's dissertation, São Paulo State University]. São Paulo State University’s repository (in Portuguese). https://repositorio.unesp.br/handle/11449/215781
    » https://repositorio.unesp.br/handle/11449/215781

Publication Dates

  • Publication in this collection
    08 Dec 2023
  • Date of issue
    2023

History

  • Received
    15 May 2023
  • Accepted
    25 Oct 2023
Associação Brasileira de Mecânica dos Solos Av. Queiroz Filho, 1700 - Torre A, Sala 106, Cep: 05319-000, São Paulo - SP - Brasil, Tel: (11) 3833-0023 - São Paulo - SP - Brazil
E-mail: secretariat@soilsandrocks.com